多目标粒子群matlab代码【精品文档】(完整版)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 改进的多目标粒子群算法，包括多个测试函数
% 对程序中的部分参数进行修改将更好地求解某些函数
%
ZDT1NP=cell(1,50);
ZDT1FV=cell(1,50);
ZDT1T=zeros(1,50);
for i=1:50
tic;
%[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('ZDT1',0.1,50,100,2.0,1.0,0.4,200,30,zer os(1,30),ones(1,30));%--ZDT1
elapsedTime=toc;
ZDT1NP(i)={np};
ZDT1FV(i)={fv};
ZDT1T(i)=elapsedTime;display(strcat('ZDT1',num2str(i)));
end
zdt1fv=cell2mat(ZDT1FV');
zdt1fv=GetLeastFunctionValue(zdt1fv);
ZDT2NP=cell(1,50);
ZDT2FV=cell(1,50);
ZDT2T=zeros(1,50);
for i=1:50
tic;
%[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('ZDT2',0.1,50,100,2.0,1.0,0.4,200,30,zer os(1,30),ones(1,30),[1,zeros(1,29)]);%--ZDT2
elapsedTime=toc;
ZDT2NP(i)={np};
ZDT2FV(i)={fv};
ZDT2T(i)=elapsedTime;display(strcat('ZDT2',num2str(i)));
end
zdt2fv=cell2mat(ZDT2FV');
zdt2fv=GetLeastFunctionValue(zdt2fv); %%%%%%%%%%%%%%%%%%%%%%%%%%%5
ZDT3NP=cell(1,50);
ZDT3FV=cell(1,50);
ZDT3T=zeros(1,50);
for i=1:50
tic;
%
[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('ZDT3',0.1,50,100,2.0,1.0,0.4,400,30,zeros(1,30 ),ones(1,30));%--ZDT3
elapsedTime=toc;
ZDT3NP(i)={np};
ZDT3FV(i)={fv};
ZDT3T(i)=elapsedTime;display(strcat('ZDT3',num2str(i)));
end
zdt3fv=cell2mat(ZDT3FV');
zdt3fv=GetLeastFunctionValue(zdt3fv);
ZDT4NP=cell(1,50);
ZDT4FV=cell(1,50);
ZDT4T=zeros(1,50);
for i=1:50
tic;
%
[np,nprule,dnp,fv,goals]=ParticleSwarmOpt('ZDT4',0.1,50,100,2.0,1.0,0.4,200,10,[0,-5,-5,-5,-5,-5, -5,-5,-5,-5],[1,5,5,5,5,5,5,5,5,5],[1,0,0,0,0,0,0,0,0,0]);%--ZDT4
elapsedTime=toc;
ZDT4NP(i)={np};
ZDT4FV(i)={fv};
ZDT4T(i)=elapsedTime;display(strcat('ZDT4',num2str(i)));
end
zdt4fv=cell2mat(ZDT4FV');
zdt4fv=GetLeastFunctionValue(zdt4fv); %%%%%%%%%%%%%%%%%%%%%%%%
ZDT6NP=cell(1,50);
ZDT6FV=cell(1,50);
ZDT6T=zeros(1,50);
for i=1:50
tic;
%[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('ZDT6',0.1,50,100,2.0,1.0,0.4,200,10,zer os(1,10),ones(1,10));%--ZDT6
elapsedTime=toc;
ZDT6NP(i)={np};
ZDT6FV(i)={fv};
ZDT6T(i)=elapsedTime;display(strcat('ZDT6',num2str(i)));
end
zdt6fv=cell2mat(ZDT6FV');
zdt6fv=GetLeastFunctionValue(zdt6fv);
CTP1NP=cell(1,50);
CTP1FV=cell(1,50);
CTP1T=zeros(1,50);
for i=1:50
tic;
[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('CTP1',0.1,50,100,2.0,1.0,0.4,1500,5,[0,-5,-5,-5, -5],[1,5,5,5,5],[],struct('isfmopso',false,'istargetdis',false,'stopatborder',true));%--CTP1 elapsedTime=toc;
CTP1NP(i)={np};
CTP1FV(i)={fv};
CTP1T(i)=elapsedTime;display(strcat('CTP1',num2str(i)));
end
ctp1fv=cell2mat(CTP1FV');
ctp1fv=GetLeastFunctionValue(ctp1fv);
CTP1fmNP=cell(1,50);
CTP1fmFV=cell(1,50);
CTP1fmT=zeros(1,50);
for i=1:50
tic;
[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('CTP1',0.1,50,100,2.0,1.0,0.4,400,5,[0,-5,-5,-5,-5],[1,5,5,5,5],[0 0 0 0 0],struct('isfmopso',true,'istargetdis',false,'stopatborder',true));%--CTP1 elapsedTime=toc;
CTP1fmNP(i)={np};
CTP1fmFV(i)={fv};
CTP1fmT(i)=elapsedTime;display(strcat('CTP1fm',num2str(i)));
end
ctp1fmfv=cell2mat(CTP1fmFV');
ctp1fmfv=GetLeastFunctionValue(ctp1fmfv);
CTP2NP=cell(1,50);
CTP2FV=cell(1,50);
CTP2T=zeros(1,50);
for i=1:50
tic;
[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('CTP2',0.1,50,100,2.0,1.0,0.4,1500,5,[0,-5,-5,-5, -5],[1,5,5,5,5],[],struct('isfmopso',false,'istargetdis',false,'stopatborder',true));%--CTP2 elapsedTime=toc;
CTP2NP(i)={np};
CTP2FV(i)={fv};
CTP2T(i)=elapsedTime;display(strcat('CTP2',num2str(i)));
end
ctp2fv=cell2mat(CTP2FV');
ctp2fv=GetLeastFunctionValue(ctp2fv);
CTP2fmNP=cell(1,50);
CTP2fmFV=cell(1,50);
CTP2fmT=zeros(1,50);
for i=1:50
tic;
[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('CTP2',0.1,50,100,2.0,1.0,0.4,400,5,[0,-5,-5,-5,-5],[1,5,5,5,5],[0 0 0 0 0],struct('isfmopso',true,'istargetdis',false,'stopatborder',true));%--CTP2 elapsedTime=toc;
CTP2fmNP(i)={np};
CTP2fmFV(i)={fv};
CTP2fmT(i)=elapsedTime;display(strcat('CTP2fm',num2str(i)));
end
ctp2fmfv=cell2mat(CTP2fmFV');
ctp2fmfv=GetLeastFunctionValue(ctp2fmfv);
CTP3NP=cell(1,50);
CTP3FV=cell(1,50);
CTP3T=zeros(1,50);
for i=1:50
tic;
[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('CTP3',0.1,50,100,2.0,1.0,0.4,1400,5,[0,-5,-5,-5, -5],[1,5,5,5,5],[],struct('isfmopso',false,'istargetdis',false,'stopatborder',true));%--CTP3 elapsedTime=toc;
CTP3NP(i)={np};
CTP3FV(i)={fv};
CTP3T(i)=elapsedTime;display(strcat('CTP3',num2str(i)));
end
ctp3fv=cell2mat(CTP3FV');
ctp3fv=GetLeastFunctionValue(ctp3fv);
CTP3fmNP=cell(1,50);
CTP3fmFV=cell(1,50);
CTP3fmT=zeros(1,50);
for i=1:50
tic;
[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('CTP3',0.1,50,100,2.0,1.0,0.4,400,5,[0,-5,-5,-5,-5],[1,5,5,5,5],[0 0 0 0 0],struct('isfmopso',true,'istargetdis',false,'stopatborder',true));%--CTP3 elapsedTime=toc;
CTP3fmNP(i)={np};
CTP3fmFV(i)={fv};
CTP3fmT(i)=elapsedTime;display(strcat('CTP3fm',num2str(i)));
end
ctp3fmfv=cell2mat(CTP3fmFV');
ctp3fmfv=GetLeastFunctionValue(ctp3fmfv);
CTP4NP=cell(1,50);
CTP4FV=cell(1,50);
CTP4T=zeros(1,50);
for i=1:50
tic;
[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('CTP4',0.1,50,100,2.0,1.0,0.4,1400,5,[0,-5,-5,-5, -5],[1,5,5,5,5],[],struct('isfmopso',false,'istargetdis',false,'stopatborder',true));%--CTP4 elapsedTime=toc;
CTP4NP(i)={np};
CTP4FV(i)={fv};
CTP4T(i)=elapsedTime;display(strcat('CTP4',num2str(i)));
end
ctp4fv=cell2mat(CTP4FV');
ctp4fv=GetLeastFunctionValue(ctp4fv);
CTP4fmNP=cell(1,50);
CTP4fmFV=cell(1,50);
CTP4fmT=zeros(1,50);
for i=1:50
tic;
[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('CTP4',0.1,50,100,2.0,1.0,0.4,400,5,[0,-5,-5,-5,-5],[1,5,5,5,5],[0 0 1 0 0],struct('isfmopso',true,'istargetdis',false,'stopatborder',true));%--CTP4 elapsedTime=toc;
CTP4fmNP(i)={np};
CTP4fmFV(i)={fv};
CTP4fmT(i)=elapsedTime;display(strcat('CTP4fm',num2str(i)));
end
ctp4fmfv=cell2mat(CTP4fmFV');
ctp4fmfv=GetLeastFunctionValue(ctp4fmfv);
CTP5NP=cell(1,50);
CTP5FV=cell(1,50);
CTP5T=zeros(1,50);
for i=1:50
tic;
%[np,nprule,dnp,fv,goals,pbest]=ParticleSwarmOpt('CTP5',0.1,50,100,2.0,1.0,0.4,200,5,[0,-5 ,-5,-5,-5],[1,5,5,5,5],[0 0 0 0 0]);%--CTP5
elapsedTime=toc;
CTP5NP(i)={np};
CTP5FV(i)={fv};