一、多目标粒子群优化算法MOPSO

二、部分代码

MOPSO求解46个多目标测试函数（ZDT1、ZDT2、ZDT3、ZDT4、ZDT6、DTLZ1-DTLZ7、WFG1-WFG10、UF1-UF10、CF1-CF10、Kursawe、Poloni、Viennet2、Viennet3）以及1个工程应用（盘式制动器设计），并采用四种评价指标（IGD，GD，HV，SP）评价。

close all;
clear ; 
clc;
%%
% TestProblem测试问题说明：1-47
%%
TestProblem=1;
MultiObj = GetFunInfo(TestProblem);
MultiObjFnc=MultiObj.name;%问题名
% Parameters
params.Np = 100;        % Population size
params.Nr = 200;        % Repository size
params.maxgen = 80;    % Maximum number of generations



% MOPSO
REP = MOPSO(params,MultiObj);
Xbest=REP.pos;
Fbest=REP.pos_fit;
%% 画结果图
figure(2)
if(size(Fbest,2)==2)
    h_rep = plot(Fbest(:,1),Fbest(:,2),'ok'); hold on;
       if(isfield(MultiObj,'truePF'))
            h_pf = plot(MultiObj.truePF(:,1),MultiObj.truePF(:,2),'.r'); hold on;
            legend('MOPSO','TruePF');
       else
           legend('MOPSO');
       end


        grid on; xlabel('f1'); ylabel('f2');
end
if(size(Fbest,2)==3)
    h_rep = plot3(Fbest(:,1),Fbest(:,2),Fbest(:,3),'ok'); hold on;
       if(isfield(MultiObj,'truePF'))
            h_pf = plot3(MultiObj.truePF(:,1),MultiObj.truePF(:,2),MultiObj.truePF(:,3),'.r'); hold on;
            legend('MOPSO','TruePF');
       else
           legend('MOPSO');
       end


        grid on; xlabel('f1'); ylabel('f2');zlabel('f3');
end


Obtained_Pareto=Fbest;
if(isfield(MultiObj,'truePF'))%判断是否有参考的PF
True_Pareto=MultiObj.truePF;
%%  Metric Value
% ResultData的值分别是IGD、GD、HV、Spacing  (HV越大越好，其他指标越小越好)
ResultData=[IGD(Obtained_Pareto,True_Pareto),GD(Obtained_Pareto,True_Pareto),HV(Obtained_Pareto,True_Pareto),Spacing(Obtained_Pareto)];
else
    %计算每个算法的Spacing，Spacing越小说明解集分布越均匀
    ResultData=Spacing(Obtained_Pareto);%计算的Spacing
end
% Display info
disp('Repository fitness values are stored in REP.pos_fit');
disp('Repository particles positions are store in REP.pos');




    

三、部分结果

四、完整MATLAB代码

文件夹包含MOPSO所有代码（含46个多目标测试函数+1个工程应用+4种评价指标）、原始文献以及多目标测试函数的TruePF。代码可点击main.m运行.