clc
clear all
% 鐢诲浘
x=0:0.1:20;
y=0:0.1:20;
[x,y]=meshgrid(x,y);
z= (x .* sin(x) .* cos(2 * x) - 2 * x .* sin(3 * x)).*(y .* sin(y) .* cos(2 * y) - 2 * y .* sin(3 * y));
surf(x,y,z);
set(gcf,'position',[0 0 1920 1080]);%閾烘弧灞忓箷
colormap winter;%鏀逛竴鏇查潰鐨勯鑹?
hold on
%鍙傛暟鍒濆鍖?
c1 = 1.49445;
c2 = 1.49445;
maxgen = 100; %杩涘寲娆℃暟
sizepop =100; %绉嶇兢瑙勬ā
Vmax = 10; %閫熷害鏈?澶у??
Vmin =-10; %閫熷害鏈?灏忓??
popmax =20;
popmin =0;
%浜х敓鍒濆绮掑瓙鍜岄?熷害
%rand鍑芥暟鏄?0-1锛岄偅涔堝浣曚骇鐢熷埌锛坅锛宐锛夌殑鍖洪棿鐨勯殢鏈烘暟鍛€?俽and*锛坆-a锛?+a
for i=1:sizepop
pop(i,:)=rand(1,2)*20;
V(i,:)=rand(1,2)*(Vmax-Vmin)+Vmin;
fitness(i)=fun(pop(i,1),pop(i,2));
end
%%%鎵惧埌閫傚簲搴︽渶寮虹殑
gbest=pop; %%
fitnessgbest=fitness; %%涓綋鏈?寮?
[bestfitness,bestindex] = max(fitness);%%缇や綋鏈?寮?
fitnesszbest=fitness(bestindex);
zbest=pop(bestindex,:);
%%杩唬瀵讳紭
%鎻忕偣
temp = plot3(pop(:,1),pop(:,2),fitness,'y.','markersize',40);
temp2 = plot3(zbest(1),zbest(2),fitnesszbest,'r.','markersize',40);
hold on
for i=1:maxgen
for j=1:sizepop
V(j,:) = V(j,:) + c1*rand*(gbest(j,:)-pop(j,:)) + c2*rand*(zbest-pop(j,:));
%妫?鏌ラ?熷害鏄惁瓒婄晫
if(V(j,1)>Vmax) V(j,1)=Vmax;end
if(V(j,2)>Vmax) V(j,2)=Vmax;end
if(V(j,1)<Vmin) V(j,1)=Vmin;end
if(V(j,2)<Vmin) V(j,2)=Vmin;end
pop(j,:)=pop(j,:)+V(j,:);
%妫?鏌ュ彉閲忔槸鍚﹁秺鐣?
if(pop(j,1)>popmax) pop(j,1)=popmax;end
if(pop(j,2)>popmax) pop(j,2)=popmax;end
if(pop(j,1)<popmin) pop(j,1)=popmin;end
if(pop(j,2)<popmin) pop(j,2)=popmin;end
%閫傚簲搴︽洿鏂?
fitness(j)=fun(pop(j,1),pop(j,2));
%%鏇存柊涓綋鏈?浼?
if(fitness(j)>fitnessgbest(j))
fitnessgbest(j)=fitness(j);
gbest(j,:)=pop(j,:);
end
end
%%鏇存柊缇や綋鏈?浼?
[bestfitness,bestindex] = max(fitnessgbest);%%缇や綋鏈?寮?
fitnesszbest=fitnessgbest(bestindex);
zbest=gbest(bestindex,:);
%%鍔ㄦ?佺敾鍥?
title(strcat('第',num2str(i),'代'),'fontsize',40)%鍔ㄦ?佹弿缁樻爣棰?
set(temp,'visible','off');
set(temp2,'visible','off');%%鎶婁箣鍓嶇敾鐨勭偣闅愯棌
temp = plot3(pop(:,1),pop(:,2),fitness,'y.','markersize',50);
temp2 = plot3(zbest(1),zbest(2),fitnesszbest,'r.','markersize',50);%% 閲嶆柊鐢荤偣
hold on
pause(0.3);
end
