1 %----------------------------------------------------------------------

2 %程序功能：实现MOEAD算法，测试函数为ZDT1,ZDT2,ZDT3,ZDT4,ZDT6,DTLZ1,DTLZ23 %说明：遗传算子为模拟二进制交叉和多项式变异4 %作者：(晓风)5 %email: 18821709267@163.com6 %最初建立时间：2018.09.30

7 %最近修改时间：2018.10.08

8 %参考论文：9 %MOEA/D: A Multiobjective Evolutionary Algorithm Based on Decomposition10 %Qingfu Zhang, Senior Member, IEEE, andHui Li11 %IEEE TRANSACTIONS O12 %----------------------------------------------------------

13 clear all14 clc15 tic;16 %------------------------参数输入--------------------------

17 format long18 globalx_max x_min x_num f_num lamda z19 rand('state',sum(100*clock));20 N=300;%种群大小21 T=20;%邻居规模大小22 fun='DTLZ2';%

23 funfun;%测试函数24 lamda=genrate_lamda(N,f_num);%均匀分布的N个权重向量25 max_gen=250;%进化代数26 pc=1;%交叉概率27 pm=1/x_num;%变异概率28 yita1=2;%模拟二进制交叉参数229 yita2=5;%多项式变异参数530 %------------------------初始条件--------------------------

31 %%计算任意两个权重向量间的欧式距离，查找每个权向量最近的T个权重向量的索引32 B=look_neighbor(lamda,T);33 %%在可行空间均匀随机产生初始种群34 X=initialize(N,f_num,x_num,x_min,x_max,fun);35 %%初始化z36 for i=1:f_num37 z(i) = min(X(:,x_num+i));38 end39 %%初始化是否为非支配个体40 X=deterdomination(X,N,f_num,x_num);41 %%设置EP为初始种群里的非支配个体42 EP=[];43 for i=1:N44 if(X(i,x_num+f_num+1)==1)45 EP=[X(i,:);EP];46 end47 end48 %------------------------迭代更新--------------------------

49 for gen=1:max_gen50 for i=1:N51 %%基因重组，从B(i)中随机选取两个序列k，l52 index1 =randperm(T);53 parent1 = B(i,index1(1));54 parent2 = B(i,index1(2));55 off=cross_mutation(X(parent1,:),X(parent2,:),f_num,x_num,x_min,x_max,pc,pm,yita1,yita2,fun );56 %off=cross_mutation2(X(parent1,:),X(parent2,:),f_num,x_num,x_min,x_max,pc,pm,yita1,yita2,fun );57 %%更新z58 for j=1:f_num59 %%if(Zi

70 equal=0;%y'的目标函数值等于个体的目标函数值数目

71 greater=0;%y'的目标函数值大于个体的目标函数值数目

72 for mm=1:f_num73 if(off(:,x_num+mm)>EP(k,x_num+mm))74 greater=greater+1;75 elseif(off(:,x_num+mm)==EP(k,x_num+mm))76 equal=equal+1;77 else

78 less=less+1;79 end80 end81 %%%从EP中移除被y'支配的向量

82 if(greater==0 && equal~=f_num)%y'支配EP中的第K个个体

83 kk=[k kk];84 end85 %%%如果EP中没有支配y'的个体，将y'加入EP86 if(less==0 && equal~=f_num)%EP中的第K个个体支配y'87 temp=1;88 end89 end90 if(isempty(kk)==0)91 EP(kk,:)=[];92 end93 if(temp==0)94 EP=[EP;off];95 end96 end97 if mod(gen,10) ==098 fprintf('%d gen has completed!\n',gen);99 end100 % if f_num==2

101 % plot(EP(:,x_num+1),EP(:,x_num+2),'r*');102 % elseif f_num==3

103 % plot3( EP(:,x_num+1), EP(:,x_num+2),EP(:,x_num+3),'r*');104 % set(gca,'xdir','reverse'); set(gca,'ydir','reverse');105 %end106 %title(num2str(gen));107 %drawnow108 end109 filepath=pwd;110 cd('G:\xjl\learnPro\matlabLearn\临时数据\EP_DTLZ2');111 save solution5.txt EP -ASCII112 cd(filepath);113 toc;114 %------------------------画图对比--------------------------

115 if f_num==2

116 hold on117 plot(EP(:,x_num+1),EP(:,x_num+2),'r*');118 elseif f_num==3

119 hold on120 plot3( EP(:,x_num+1), EP(:,x_num+2),EP(:,x_num+3),'r*');121 set(gca,'xdir','reverse'); set(gca,'ydir','reverse');122 end123 %figure;124 % if f_num==2

125 %hold on126 % plot(X(:,x_num+1),X(:,x_num+2),'r*');127 % elseif f_num==3

128 %hold on129 % plot3( X(:,x_num+1), X(:,x_num+2),X(:,x_num+3),'r*');130 % set(gca,'xdir','reverse'); set(gca,'ydir','reverse');131 %end132 %--------------------Coverage(C-metric)---------------------

133 A=PP;B=EP(:,(x_num+1):(x_num+f_num));%%%%%%%%%%%%%%%%%%%%

134 [temp_A,~]=size(A);135 [temp_B,~]=size(B);136 number=0;137 for i=1:temp_B138 nn=0;139 for j=1:temp_A140 less=0;%当前个体的目标函数值小于多少个体的数目141 equal=0;%当前个体的目标函数值等于多少个体的数目142 for k=1:f_num143 if(B(i,k)

151 end152 end153 if(nn~=0)154 number=number+1;155 end156 end157 C_AB=number/temp_B;158 disp('C_AB:');159 disp(C_AB);160 %-----Distance from Representatives in the PF(D-metric)-----

161 A=EP(:,(x_num+1):(x_num+f_num));P=PP;%%%%%%%%%%%%%%%%%%%

162 [temp_A,~]=size(A);163 [temp_P,~]=size(P);164 min_d=0;165 for v=1:temp_P166 d_va=(A-repmat(P(v,:),temp_A,1)).^2;167 min_d=min_d+min(sqrt(sum(d_va,2)));168 end169 D_AP=(min_d/temp_P);170 disp('D_AP:');171 disp(D_AP);172 filepath=pwd;173 cd('G:\xjl\learnPro\matlabLearn\临时数据\EP_DTLZ2');174 save C_AB5.txt C_AB -ASCII175 save D_AP5.txt D_AP -ASCII176 cd(filepath);