c/c++语言开发共享基于Matlab实现多目标粘液霉菌算法的示例代码

%%  (zdt3)  function f = zdt3 (x)  % 两个目标函数  % 变量数目是 30. x [0,1]  f = [];  n=length(x);  g=1+9*sum(x(2:n))/(n-1);  f(1)=x(1);  f(2)=1-sqrt(x(1)/g)-(x(1)/g)*sin(10*pi*x(1));

%%===mosma：多目标粘液霉菌算法=====     %% 目标函数  % 目标函数描述包含了有关目标函数的信息。m为目标空间的维数，d为决策变量空间的维数，  % lb和ub为决策变量空间中变量的取值范围。  % 用户必须使用决策变量定义目标函数。务必编辑功能‘评估_目标'，以适应您的需要。  clc  clear all  d = 30; % 决策变量数目  m = 2; % 目标函数的个数  k=m+d;  lb = ones(1, d).*0; %  lb - 表示每个决策变量的最小值的十进制值的向量。  ub = ones(1, d).*1; % ub - 决策变量的最大可能值向量。  gen = 200;  % 设置最大生成数  ecosize = 200;      % 种群数量 (np)  ishow = 10;  %% 开始进化过程  pareto = mosma(d,m,lb,ub,ecosize,gen,ishow);  obtained_pareto= pareto(:,d+1:d+m); % extract data to plot  obtained_pareto=sortrows(obtained_pareto,2);  true_pareto=load('zdt3.txt');  %% 可视化  %% =========两个目标函数时======================  if m == 2      plot(obtained_pareto(:,1),obtained_pareto(:,2),'o','linewidth',2,...          'markeredgecolor','r','markersize',2);      hold on      plot(true_pareto(:,1),true_pareto(:,2),'k');       title('采用mosma算法求解pareto最优解集');      legend('mosma');      xlabel('f_1');      ylabel('f_2');  %% =====3个目标函数时==========================  elseif m == 3         plot3(obtained_pareto(:,1),obtained_pareto(:,2),obtained_pareto(:,3),'o','linewidth',2,...          'markeredgecolor','r','markersize',2);      hold on      plot3(obtained_pareto(:,1),obtained_pareto(:,2),obtained_pareto(:,3),'.','linewidth',2,...          'markeredgecolor','k','markersize',6);      title('采用mosma算法求解pareto最优解集');      legend('mosma');      xlabel('f_1');      ylabel('f_2');      zlabel('f_3');  end  %%  metric value  m_igd=igd(obtained_pareto,true_pareto);  m_gd=gd(obtained_pareto,true_pareto);  m_hv=hv(obtained_pareto,true_pareto);  m_spacing=spacing(obtained_pareto,true_pareto);  m_spread=spread(obtained_pareto,true_pareto);  m_deltap=deltap(obtained_pareto,true_pareto);  display(['the igd metric obtained by mosma is     : ', num2str(m_igd)]);  display(['the gd metric obtained by mosma is      : ', num2str(m_gd)]);  display(['the hv metric obtained by mosma is      : ', num2str(m_hv)]);  display(['the spacing metric obtained by mosma is : ', num2str(m_spacing)]);  display(['the spread metric obtained by mosma is  : ', num2str(m_spread)]);  display(['the deltap metric obtained by mosma is  : ', num2str(m_deltap)]);