多目標粒子群matlab代碼-資料下載頁

　　

【正文】 數 case {39。ZDT139。,39。ZDT239。,39。ZDT339。,39。ZDT439。,39。ZDT539。,39。ZDT639。} fai=0。%CTP1 case 39。CTP139。 n=length(x)。 cx=41+sum(x(2:n).^210*cos(2*pi*x(2:n)))。 f1=x(1)。 f2=cx*exp(f1/cx)。 g1=*exp(*f1)f2。 g2=*exp(*f1)f2。 fai=max(g1,0)+max(g2,0)。%CTP2,CTP3,CTP4,CTP5,CTP6,CTP7 case {39。CTP239。,39。CTP339。,39。CTP439。,39。CTP539。,39。CTP639。,39。CTP739。} navars=[*pi,10,1,6,1。 *pi,10,1,1。 *pi,40,5,1,6,0。 *pi,10,2,1。 *pi,40,1,2,2。 *pi,10,1,1。]。 i=str2double(funame(4))1。%%選擇的是第幾個測試函數 sita=navars(i,1)。a=navars(i,2)。b=navars(i,3)。 c=navars(i,4)。d=navars(i,5)。e=navars(i,6)。 n=length(x)。 cx=41+sum(x(2:n).^210*cos(2*pi*x(2:n)))。 f1=x(1)。 f2=cx*(1f1/cx)。 g=a*abs(sin(b*pi*(sin(sita)*(f2e)+cos(sita)*f1)^c))^d(cos(sita)*(f2e)sin(sita)*f1)。 fai=max(g,0)。%CTP8,CONSTR case {39。CTP839。,39。CONSTR39。} g1=6x(2)9*x(1)。 g2=1+x(2)9*x(1)。 fai=max(g1,0)+max(g2,0)。%CTP9,SRN case {39。CTP939。,39。SRN39。} g1=x(1)^2+x(2)^2255。 g2=x(1)3*x(2)+10。 fai=max(g1,0)+max(g2,0)。%CTP10,TNK,MOPC4 case {39。CTP1039。,39。TNK39。,39。MOPC439。} g1=x(1)^2x(2)^2+1+*cos(16*atan(x(1)/x(2)))。 g2=(x(1))^2+(x(2))^。 fai=max(g1,0)+max(g2,0)。%MOPC1,BNH函數 case {39。MOPC139。,39。BNH39。} g1=(x(1)5)^2+x(2)^225。 g2=(x(1)8)^2(x(2)+3)。 fai=max(g1,0)+max(g2,0)。%MOPC2,OSY函數 case {39。MOPC239。,39。OSY39。} g1=2x(1)x(2)。 g2=6+x(1)+x(2)。 g3=2x(1)+x(2)。 g4=2+x(1)3*x(2)。 g5=4(x(3)3)^2+x(4)。 g6=4(x(5)3)^2x(6)。 fai = max(g1,0)+max(g2,0)+max(g3,0)+max(g4,0)+max(g5,0)+max(g6,0)。%MOPC3 case 39。MOPC339。 g1=x(2)+4*x(1)4。 g2=x(1)x(2)2。 fai = max(g1,0)+max(g2,0)。endendfunction fvout=GetLeastFunctionValue(fvin)%將外部種群中的非支配集剔除fvout=fvin。n=size(fvout,1)。i=1。while(i=n) j=i+1。 isdominated=false。 while(j=n) a=fvout(i,:)。b=fvout(j,:)。 if((a(1)b(1)amp。amp。a(2)=b(2))||(a(1)=b(1)amp。amp。a(2)b(2))) fvout(j,:)=[]。n=n1。 else if((b(1)a(1)amp。amp。b(2)=a(2))||(b(1)=a(1)amp。amp。b(2)a(2))) isdominated=true。 end j=j+1。 end end if isdominated fvout(i,:)=[]。n=n1。 else i=i+1。 endendendfunction [SP,SP1,MS,GD,GD2]=GetEvaluDis(funame,NP,stdfv)%計算所有評價指標的值if(iscell(NP)) n=length(NP)。 SP=Inf*ones(1,n)。 SP1=Inf*ones(1,n)。 MS=Inf*ones(1,n)。 GD=Inf*ones(1,n)。 GD2=Inf*ones(1,n)。 for i=1:n np=cell2mat(NP(i))。 [sp,sp1]=Spacing(funame,np)。 ms=MaximumSpread(funame,np)。 gd=GenerationalDistance(funame,np,stdfv)。 gd2=GenerationalDistance2(funame,np)。 SP(i)=sp。 SP1(i)=sp1。 MS(i)=ms。 GD(i)=gd。 GD2(i)=gd2。 endelse [SP,SP1]=Spacing(funame,NP)。 MS=MaximumSpread(funame,NP)。 GD=GenerationalDistance(funame,NP,stdfv)。 GD2=GenerationalDistance2(funame,NP)。endendfunction [SP,SP1]=Spacing(funame,np)%分散性(Spacing,SP)n=size(np,1)。m=size(fitness(np(1,:),funame),2)。f=zeros(n,m)。for i=1:n f(i,:)=fitness(np(i,:),funame)。endn=size(f,1)。%n:解的個數；dd=ones(1,n)*Inf。d=zeros(1,n)。for i=1:n k=0。 for j=1:n if j==i continue。 end k=k+1。 dd(i,k)=sum(abs(f(i,:)f(j,:)))。 end d(i)=min(dd(i,:))。enddaver=mean(d)。SP=sqrt(sum((ddaver).^2)/n)。SP1=sqrt(sum((ddaver).^2)/(n1))/daver。endfunction D=MaximumSpread(funame,np)%最大散布范圍(Maximum Spread,MS)n=size(np,1)。m=size(fitness(np(1,:),funame),2)。f=zeros(n,m)。if(n==1) D=0。 return。endfor i=1:n f(i,:)=fitness(np(i,:),funame)。endfmax=Inf*ones(1,m)。fmin=Inf*ones(1,m)。switch upper(funame)%ZDT1,39。ZDT239。,ZDT4 case {39。ZDT139。,39。ZDT239。,39。ZDT439。} fmax=[1,1。]。 fmin=[0,0。]。%ZDT3 case 39。ZDT339。 fmax=[,1]。 fmin=[0,]。%ZDT6 case 39。ZDT639。 fmax=[1,]。 fmin=[,0]。 case 39。CTP139。 fmax=[1,1]。 fmin=[0,]。 case 39。CTP239。 fmax=[,1]。 fmin=[0,]。 case 39。CTP339。 fmax=[,1]。 fmin=[0,]。 case 39。CTP439。 fmax=[1,]。 fmin=[0,0]。 case 39。CTP539。 case 39。CTP639。 case 39。CTP739。 case {39。CTP839。,39。CONSTR39。} case {39。CTP939。,39。SRN39。} case {39。CTP1039。,39。TNK39。,39。MOPC439。} fmax=[,]。 fmin=[,]。 case {39。MOPC239。,39。OSY39。} fmax=[,]。 fmin=[,]。 endfimax=min([max(f)。fmax])。fimin=max([min(f)。fmin])。fimin=min([fimax。fimin])。D=sqrt(sum((fimaxfimin)./(fmaxfmin))/m)。endfunction GD=GenerationalDistance(funame,np,stdfv)%世代距離GDn=size(np,1)。m=size(fitness(np(1,:),funame),2)。f=zeros(n,m)。if(n==1) GD=0。 return。endfor i=1:n f(i,:)=fitness(np(i,:),funame)。endn=size(f,1)。%n:解的個數；d=zeros(1,n)。for i=1:n d(i)=GetTrueDis(funame,f(i,:),stdfv)。endGD=sqrt(sum(d.^2))/n。endfunction d=GetTrueDis(funame,fv,stdfv)%點到真實Pareto前沿的距離x0=fv(1)。y0=fv(2)。dd=zeros(1,1001)。switch upper(funame)%ZDT1,ZDT4 case {39。ZDT139。,39。ZDT439。} x1=(1y0)^2。 for i=1:1001 x=x0+(x1x0)/1000*(i1)。 y=1sqrt(x)。 dd(i)=sqrt((xx0)^2+(yy0)^2)。 end%ZDT2 case {39。ZDT239。,39。ZDT639。} if(y0) x1=。 else x1=sqrt(1y0)。 end for i=1:1001 x=x0+(x1x0)/1000*(i1)。 y=1x^2。 dd(i)=sqrt((xx0)^2+(yy0)^2)。