【正文】 gx=fitness(x,funame)。 np_out(n_min(tempIndex),:)=[]。 dnp_out(j,r1+1)=dnp_out(r1+1,j)。 if faix 0 nprule_out(r1+1,:)=[0,faix,0]。 end elseif domiRel==1 %x被NP(i)支配,返回不再比較 return。 if domiRel==1 %NP(i)被x支配 np_out(i,:)=[]。%非劣解復(fù)本nprule_out = nprule。for i=1:r fv(i,:)=fitness(np(i,:),funame,params)。%非劣解nprule=NPRule。 elseif(randi([0,1],1)==0)%新解與pbest互相不支配 pbest(i,:) = x(i,:)。 end end end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %每個(gè)粒子的目標(biāo)向量//第8步/////////////// goals(t,i,:)=fitness(x(i,:),funame,params)。 x(i,j)=x(i,j)+v(i,j)。 else x(i,j)=lb(j)+(ub(j)lb(j))*rand。 else x(i,j)=lb(j)+(ub(j)lb(j))*rand。CTP39。 elseif v(i,j)vmin(j) v(i,j)=vmin(j)。 x(i,:)=nploc(seltemp,:)。 if faix=unfitx [nploc,npruleloc,dnploc] = pare(tempx(k,:),nploc,npruleloc,dnploc,10,funame,params)。 for j=1:D if tempv(k,j)vmax(j) tempv(k,j)=vmax(j)。 dnploc(1,1)=0。 %w1=Inf。 %獲得全局最優(yōu)/////第5步///////////////////////////////////////// %if mod(t,3)==1 %[gbest,NPRule] = GetGlobalBest(NP,NPRule,Dnp)。 w1=w()*t/M。%自身最優(yōu)解%在確定每個(gè)粒子所對(duì)就的目標(biāo)方格//第4步///////////////////////////%進(jìn)入主要循環(huán)，按照公式依次迭代for t=1:M if mod(t,100)==0 unfitx = *(t+200)/M。%非劣解集參數(shù)Dnp(1,1)=0。 %隨機(jī)初始化速度 endend%計(jì)算目標(biāo)向量%速度控制vmax=(ublb)*。%記下N個(gè)粒子M次迭代T維目標(biāo)變化%初始化種群的個(gè)體///////第1步///////////////////////////////////x(1,:)=x0。,true)。,false,39。NP=[]。endosyfmfv=cell2mat(OSYfmFV39。 OSYfmT(i)=elapsedTime。,true))。,true,39。 [np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt(39。OSYfmNP=cell(1,50)。,num2str(i)))。 OSYFV(i)={fv}。stopatborder39。isfmopso39。for i=1:50 tic。bnhfv=GetLeastFunctionValue(bnhfv)。BNH39。 BNHNP(i)={np}。 %[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt(39。BNHNP=cell(1,50)。,num2str(i)))。 TNKfmFV(i)={fv}。stopatborder39。isfmopso39。for i=1:50 tic。tnkfv=GetLeastFunctionValue(tnkfv)。CTP1039。 TNKNP(i)={np}。,false,39。,50,100,1300,2,[0,0],[pi,pi],[],struct(39。TNKT=zeros(1,50)。)。display(strcat(39。%CTP9,SRN elapsedTime=toc。for i=1:50 tic。constrfv=GetLeastFunctionValue(constrfv)。CTP839。 CONSTRNP(i)={np}。 %[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt(39。CONSTRNP=cell(1,50)。,num2str(i)))。 CTP7FV(i)={fv}。CTP739。CTP7FV=cell(1,50)。endctp6fv=cell2mat(CTP6FV39。 CTP6T(i)=elapsedTime。,50,100,400,5,[0,5,5,5,5],[1,5,5,5,5],[0 0 0 0 0])。CTP6T=zeros(1,50)。)。display(strcat(39。%CTP5 elapsedTime=toc。for i=1:50 tic。ctp4fmfv=GetLeastFunctionValue(ctp4fmfv)。CTP4fm39。 CTP4fmNP(i)={np}。,false,39。,50,100,400,5,[0,5,5,5,5],[1,5,5,5,5],[0 0 1 0 0],struct(39。CTP4fmT=zeros(1,50)。)。display(strcat(39。%CTP4 elapsedTime=toc。istargetdis39。CTP439。CTP4FV=cell(1,50)。endctp3fmfv=cell2mat(CTP3fmFV39。 CTP3fmT(i)=elapsedTime。,true))。,true,39。 [np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt(39。CTP3fmNP=cell(1,50)。,num2str(i)))。 CTP3FV(i)={fv}。stopatborder39。isfmopso39。for i=1:50 tic。ctp2fmfv=GetLeastFunctionValue(ctp2fmfv)。CTP2fm39。 CTP2fmNP(i)={np}。,false,39。,50,100,400,5,[0,5,5,5,5],[1,5,5,5,5],[0 0 0 0 0],struct(39。CTP2fmT=zeros(1,50)。)。display(strcat(39。%CTP2 elapsedTime=toc。istargetdis39。CTP239。CTP2FV=cell(1,50)。endctp1fmfv=cell2mat(CTP1fmFV39。 CTP1fmT(i)=elapsedTime。,true))。,true,39。 [np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt(39。CTP1fmNP=cell(1,50)。,num2str(i)))。 CTP1FV(i)={fv}。stopatborder39。isfmopso39。for i=1:50 tic。zdt6fv=GetLeastFunctionValue(zdt6fv)。ZDT639。 ZDT6NP(i)={np}。 %[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt(39。%%%%%%%%%%%%%%%%%%%%%%%%ZDT6NP=cell(1,50)。,num2str(i)))。 ZDT4FV(i)={fv}。ZDT439。ZDT4FV=cell(1,50)。endzdt3fv=cell2mat(ZDT3FV39。 ZDT3T(i)=elapsedTime。,50,100,400,30,zeros(1,30),ones(1,30))。ZDT3T=zeros(1,50)。)。display(strcat(39。%ZDT2 elapsedTime=toc。for i=1:50 tic。zdt1fv=GetLeastFunctionValue(zdt1fv)。ZDT139。 ZDT1NP(i)={np}。 %[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt(39。.. . . ..%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 改進(jìn)的多目標(biāo)粒子群算法，包括多個(gè)測(cè)試函數(shù)% 對(duì)程序中的部分參數(shù)進(jìn)行修改將更好地求解某些函數(shù)% ZDT1NP=cell(1,50)。for i=1:50 tic。%ZDT1 elapsedTime=toc。display(strcat(39。)。ZDT2T=zeros(1,50)。,50,100,200,30,zeros(1,30),ones(1,30),[1,zeros(1,29)])。 ZDT2T(i)=elapsedTime。endzdt2fv=cell2mat(ZDT2FV39。ZDT3FV=cell(1,50)。ZDT339。 ZDT3FV(i)={fv}。,num2str(i)))。ZDT4NP=cell(1,50)。% [np,nprule,dnp,fv,goals]=ParticleSwarmOpt(39。 ZDT4NP(i)={np}。ZDT439。zdt4fv=GetLeastFunctionValue(zdt4fv)。for i=1:50 tic。%ZDT6 elapsedTime=toc。display(strcat(39。)。CTP1T=zeros(1,50)。,50,100,1500,5,[0,5,5,5,5],[1,5,5,5,5],[],struct(39。,false,39。 CTP1NP(i)={np}。CTP139。ctp1fv=GetLeastFunctionValue(ctp1fv)。for i=1:50 tic。isfmopso39。stopatborder39。 CTP1fmFV(i)={fv}。,num2str(i)))。CTP2NP=cell(1,50)。 [np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt(39。,false,39。,true))。 CTP2T(i)=elapsedTime。endctp2fv=cell2mat(CTP2FV39。CTP2fmFV=cell(1,50)。CTP239。istargetdis39。%CTP2 elapsedTime=toc。display(strcat(39。)。CTP3T=zeros(1,50)。,50,100,1400,5,[0,5,5,5,5],[1,5,5,5,5],[],struct(39。,false,39。 CTP3NP(i)={np}。CTP339。ctp3fv=GetLeastFunctionValue(ctp3fv)。for i=1:50 tic。isfmopso39。stopatborder39。 CTP3fmFV(i)={fv}。,num2str(i)))。CTP4NP=cell(1,50)。 [np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt(39。,false,39。,true))。 CTP4T(i)=elapsedTime。endctp4fv=cell2mat(CTP4FV39。CTP4fmFV=cell(