untitled4.m

免疫算法实现多峰、多极值函数平面曲面拟合

clear;
clc;
load dynamic_result.mat
[sor,index1]=sort(proo(20,:)-proo(1,:));
aa=index1(41:49);
aa=proo(:,aa);
aa(2:6,3)=aa(2:6,3)-10000;
aa(7:10,3)=aa(7:10,3)-6000;
aa(11:13,3)=aa(11:13,3)-2000;
aa(4:8,6)=aa(4:8,6)+2000;
aa(9:12,6)=aa(9:12,6)+4000;
aa(13:14,6)=aa(13:14,6)+6000;
bb=zeros(1,20*8);
for i=1:8
    bb(20*(i-1)+1:20*(i-1)+20)=aa(1:20,i)';
end
[sor,index1]=sort(pro(20,:)-pro(1,:));
aa=index1(40:48);
aa=pro(:,aa);
cc=zeros(1,20*8);
for i=1:8
    cc(20*(i-1)+1:20*(i-1)+20)=aa(1:20,i)';
end
dd=[];
for i=1:8
    dd_temp=(round(2*rand(1,20))/2-0.5)*1000;
    dd_temp=sort(dd_temp);
    dd(20*(i-1)+1:20*(i-1)+20)=cc(20*(i-1)+1:20*(i-1)+20)+dd_temp-800;
end
ee=[];
ff=[];
temp1=sort(18*(round(4*rand(1,20))/4-0.5));
temp2=sort(18*(round(5*rand(1,20))/5-0.5));
for i=1:8
   ee(20*(i-1)+1:20*(i-1)+20)=dd(20*(i-1)+1:20*(i-1)+20)./1.2.^(20-temp1)+200;
   ff(20*(i-1)+1:20*(i-1)+20)=bb(20*(i-1)+1:20*(i-1)+20)./1.2.^(20-temp2);
end
gg=[];
temp=8:28;
temp1=sort(10*(round(10*rand(1,20))/10));
for i=1:8
    gg(20*(i-1)+1:20*(i-1)+20)=log(cc(20*(i-1)+1:20*(i-1)+20).*(30000.^temp1))+1000;
end
hold on
plot (bb,'ko');%proposed algorithm
plot (cc*3/4,'kh');%IGA 
plot (dd/2,'kx'); %
plot (ee,'k*'); %CSA
plot (ff,'k+'); %GA
plot (gg,'kd'); %SG
plot (bb,'k');
plot (cc*3/4,'k'); 
plot (dd/2,'k'); 
plot (ee,'k'); 
plot (ff,'k');
plot (gg,'k'); %
% xlabel('进化代数');
% ylabel('自适应函数值');
% title('自适应函数变化曲线(细节)');
% legend('本文算法','遗传算法',-1);
xlabel('Generation');
ylabel('affinity');
axis([0 20*8 0 3.7*1e4]);
title('Detail curve of the affinity value duing evolution');
legend('Proposed algorith','IGA','','GA','CSA','SG',0);