差分进化算法matlab代码.pdfVIP

差分进化算法matlab代码

差分进化算法MATLAB代码

%差分进化算法（DE）matlab代码

%%初始化

clc

clear

%数据

fobj=@(x)3*sin(10*pi*x(1))+x(1)^2+x(2)^2;

%迭代次数

MaxGenerations=50;

%种群大小

nPop=20;

%变量范围

range=3;

xmin=-range;

xmax=range;

%变量个数

nVar=2;

%F

F=0.8;

%CR

-1-

CR=0.2;

%%主体程序

BestVal=zeros(1,MaxGenerations);

fori=1:MaxGenerations

%初始化种群

population=xmin+(xmax-xmin)*rand(nPop,nVar);

%评估

forj=1:nPop

val(j)=fobj(population(j,:));

end

%找出最小的变量值

[BestVal(i),index]=min(val);

%保留最优个体

best_value=population(index,:);

%产生重组种群

forj=1:nPop

%随机选取三个不同的个体

r=floor(rand(1,3)*nPop)+1;

whilelength(unique(r))3

r=floor(rand(1,3)*nPop)+1;

end

fork=1:nVar

-2-

%生成变异个体

v(k)=

population(r(1),k)+F*(population(r(2),k)-population(r(3),k)

);

%交叉

jcross=rand;

ifjcrossCR

u(k)=v(k);

else

u(k)=population(j,k);

end

end

%将变异后的个体放入新的种群

population(j,:)=u;

value(j)=fobj(u);

end

%更新种群中的最优值

forj=1:nPop

ifvalue(j)val(j)

val(j)=value(j);

population(j,:)=u;

end

-3-

end

%保存最优解

[tem,ind]=min(val);

BestVal(i)=tem;

end

%打印结果

disp([TheBestValis,num2str(best_value)]);

disp([TheminValis

,num2str(BestVal(MaxGenerations))]);

%结果绘图

figure;

plot(BestVal,LineWidth,3);

title(Thebestvalueofgeneration);

xlabel(Generations);

ylabel(value);

gridon;。

-4-