多目標粒子群matlab代碼(編輯修改稿)

　

【文章內容簡介】 b,unfitx)。 end c = cmax (cmax cmin)*t/M。 w1=w()*t/M。 %c = cmax。 %c = cmax (cmax cmin)*mod(t,51)/50。 %w1=w()*mod(t,51)/50。 %獲得全局最優/////第5步///////////////////////////////////////// %if mod(t,3)==1 %[gbest,NPRule] = GetGlobalBest(NP,NPRule,Dnp)。 %end for i=1:N%更新粒子的位置和速度////第6步////////////////// %v(i,:)=w*v(i,:)+cmin*rand*(pbest(i,:)x(i,:))+cmax*rand*(gbestx(i,:))。 %[gbest,NPRule] = GetGlobalBest2(x(i,:),NP,NPRule)。%%%%%1217 [gbest,NPRule] = GetGlobalBest(NP,NPRule,Dnp)。 %w1=Inf。 %while(w1 || w10) % w1=+randn*。 %end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%模糊向導粒子群算法 if() nploc=x(i,:)。 npruleloc=[0,0,0]。 dnploc(1,1)=0。 tempv=zeros(10,D)。 tempx=zeros(10,D)。 for k=1:10 tempv(k,:)=w1*v(i,:)+c*rand*(pbest(i,:)x(i,:))+c*rand*(gbestx(i,:))。 for j=1:D if tempv(k,j)vmax(j) tempv(k,j)=vmax(j)。 elseif tempv(k,j)vmin(j) tempv(k,j)=vmin(j)。 end end tempx(k,:)=x(i,:)+tempv(k,:)。 faix = GetFai(x(i,:),funame,params)。 if faix=unfitx [nploc,npruleloc,dnploc] = pare(tempx(k,:),nploc,npruleloc,dnploc,10,funame,params)。 end end nnx=size(nploc,1)。 seltemp=randi([1,nnx])。 ttv=nploc(seltemp,:)x(i,:)。 x(i,:)=nploc(seltemp,:)。 if sum(abs(ttv))0 v(i,:)=ttv。 end else%%%%%%%%%%%%%%%%%%% v(i,:)=w1*v(i,:)+c*rand*(pbest(i,:)x(i,:))+c*rand*(gbestx(i,:))。 for j=1:D if v(i,j)vmax(j) v(i,j)=vmax(j)。 elseif v(i,j)vmin(j) v(i,j)=vmin(j)。 end end x(i,:)=x(i,:)+v(i,:)。 funame=upper(funame)。 if(strcmp(funame(1:3),39。CTP39。)) x(i,:)=x(i,:)+v(i,:)。 end end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%采取措施，避免粒子飛出空間////第7步///////////// %速度位置鉗制 if()%粒子隨機停留在邊界 for j=1:D if x(i,j)ub(j) if(randi([0,0],1)==0) x(i,j)=ub(j)。 v(i,j)=v(i,j)。 else x(i,j)=lb(j)+(ub(j)lb(j))*rand。 %隨機初始化位置 v(i,j)=(ub(j)lb(j))*rand*。 end end if x(i,j)lb(j) if(randi([0,0],1)==0) x(i,j)=lb(j)。 v(i,j)=v(i,j)。 else x(i,j)=lb(j)+(ub(j)lb(j))*rand。 %隨機初始化位置 v(i,j)=(ub(j)lb(j))*rand*。 end end end else%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% for j=1:D if x(i,j)ub(j)%粒子飛出上邊界 x(i,j)=x(i,j)v(i,j)。 v(i,j)=rand*(ub(j)x(i,j))。 x(i,j)=x(i,j)+v(i,j)。 elseif x(i,j)lb(j)%粒子飛出下邊界 x(i,j)=x(i,j)v(i,j)。 v(i,j)=rand*(x(i,j)lb(j))。 x(i,j)=x(i,j)v(i,j)。 end end end%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %每個粒子的目標向量//第8步/////////////// goals(t,i,:)=fitness(x(i,:),funame,params)。%調整自身//第10步///////////////// domiRel = DominateRel(pbest(i,:),x(i,:),funame,params)。%x,y的支配關系 if domiRel==1%pbest支配新解 continue。 else if domiRel==1%新解支配pbest pbest(i,:) = x(i,:)。 elseif(randi([0,1],1)==0)%新解與pbest互相不支配 pbest(i,:) = x(i,:)。 end%對NP進行更新和維護//第9步//////////////// faix = GetFai(x(i,:),funame,params)。 if faix=unfitx [NP,NPRule,Dnp] = pare(x(i,:),NP,NPRule,Dnp,Nnp,funame,params)。 end end endendnp = NP。%非劣解nprule=NPRule。dnp = Dnp。%非劣解之間的距離r=size(np,1)。fv=zeros(r,T)。for i=1:r fv(i,:)=fitness(np(i,:),funame,params)。endend%%%%%%%%%%%%%%%主函數結束%%%%%%%%%%%%%%%%%function [np_out,nprule_out,dnp_out] = pare(x,np,nprule,dnp,nnp,funame,params)%np:現有非劣解%x:需要比較的量Nnp = nnp。%非劣解集空間r=size(np,1)。%非劣解的個數np_out=np。%非劣解復本nprule_out = nprule。dnp_out = dnp。%非劣解集點之間距離if r==0 return。endfor i=r:1:1 domiRel=DominateRel(x,np(i,:),funame,params)。 if domiRel==1 %NP(i)被x支配 np_out(i,:)=[]。%非劣解剔除該解 nprule_out(i,:)=[]。 dnp_out(i,:)=[]。 if ~isempty(dnp_out) dnp_out(:,i)=[]。 end elseif domiRel==1 %x被NP(i)支配,返回不再比較 return。 endendr1=size(np_out,1)。%現有非劣解的行列np_out(r1+1,:)=x。%與所有非支配集粒子比較均占優或不可比較，則NP中加入xfaix = GetFai(x,funame,params)。 if faix 0 nprule_out(r1+1,:)=[0,faix,0]。else nprule_out(r1+1,:)=[0,0,0]。endif r1==0 dnp_out=0。endfor j=1:r1 dnp_out(r1+1,j)=GetDistance(np_out(j,:),x,funame,params)。 dnp_out(j,r1+1)=dnp_out(r1+1,j)。endif r1=Nnp %未達到非劣解種群極限 %移除密集距離最小的一個 densedis = GetDenseDis(dnp_out)。 n_min = find(min(densedis)==densedis)。%找出密度距離最小的一個 tempIndex = randi([1,length(n_min)],1)。 np_out(n_min(tempIndex),:)=[]。%非劣解剔除該解 nprule_out(n_min(tempIndex),:)=[]。 dnp_out(n_min(tempIndex),:)=[]。 if ~isempty(dnp_out) dnp_out(:,n_min(tempIndex))=[]。 endendend%%%%%%%%%%%%%%將粒子維護到外部種群%%%%%%%%%%%%%%%%%function dis=GetDistance