DE算法-matlab程序.doc

《DE算法-matlab程序.doc》由会员分享，可在线阅读，更多相关《DE算法-matlab程序.doc（12页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、Four short words sum up what has lifted most successful individuals above the crowd: a little bit more.-author-dateDE算法-matlab程序function DE(Gm,F0);function DE(Gm,F0);%差分进化算法程序基本程序%F是变异率 %Gm=1000; %最大迭代次数Gm=10000;F0=0.8;Np=100; %种群规模CR=0.9; %杂交参数G=1;%初始化代数N=10;%所求问题的维数ge=zeros(1,Gm);%各代的最优值bestx=zero

2、s(Gm,N);%各代的最优解%产生初始种群%xmin=-10;xmax=100;%带负数的下界xmin=-5.12;xmax=5.12;function y=f(XX); %y=sum(XX.4-16*XX.2+5*XX)/10+79.333; y=sum(XX.2-10.*cos(2.*pi.*XX)+10); endX0=(xmax-xmin)*rand(Np,N)+xmin; X=X0;%变异操作X1new=zeros(Np,N);%初始化X1_new=zeros(Np,N);%初始化X1=zeros(Np,N);%初始化value=zeros(1,Np);while Gxmin)&(

3、bonCR %利用二项分布来交叉 X1_new(i,:)=X(i,:); else X1_new(i,:)=X1new(i,:); end end %竞争操作 for i=1:Np if f(X1_new(i,:)f(X(i,:) X1(i,:)=X1_new(i,:); else X1(i,:)=X(i,:); end end %找出最小值 for i=1:Np value(i)=f(X1(i,:); end fmin,nmin=min(value); ge(G)=fmin; bestx(G,:)=X1(nmin,:); X=X1; trace(G,1)=G; trace(G,2)=ge(G

4、); G=G+1;end%ii=linspace(1,Np,Np);%plot(ii,ge)gmin,n=min(ge);bestvalue=gminbestsolution=bestx(n,:)%目标函数的倒数plot(trace(:,1),trace(:,2);end %求得结果%DE(1000,0.6)%ans =0.9993%best =-2.9035 -2.9035 -2.9035 -2.9035 -2.9035 -2.9035 -2.9035 -2.9035 -2.9035 -2.9035%DE(1000,0.8)%ans =0.8459%best =-2.9481 -2.9029 -3.0410 -2.9235 -2.9949 -2.9332 -3.0126 -2.9316 -2.9341 -2.8203-