showkur.m

来自「基于GPS的多目标优化及动态多目标优化源代码 里面 包含MOP, DMOP的程序」· M 代码 · 共 47 行

function [f]=ShowKUR()
% 画函数的点
v1=-20:0.001:-14;
k1=-5:0.05:5;
k2=-5:0.05:5;
k3=-5:0.05:5;

for i=1:length(v1)
    v2(i)=0;
end;
for i=1:length(k1)
    for j=1:length(k2)
      for k=1:length(k3)
        dd =     -10 * exp(-0.2 * sqrt(k1(i).^2 + k2(j).^2))  +-10 * exp(-0.2 * sqrt(k3(k).^2 + k2(j).^2));
        if (dd <-20 ) || (dd >-14) 
            continue;
        end;
 
        td = abs(k1(i)).^0.8 + 5* sin(k1(i).^3) + abs(k2(j)).^0.8 + 5* sin(k2(j).^3) + abs(k3(k)).^0.8 + 5* sin(k3(k).^3) ;
%        if floor(dd*1000) +1 >length(v2)
%           disp(floor(dd*1000) +1);
%           disp(length(v1));
%        end;
        if v2(floor((dd+20)*1000) +1) > td || v2(floor((dd+20)*1000) +1)==0
            v2(floor((dd+20)*1000) +1) = td;
        end;
    end;
end;
end;

 for i=2:length(v2)
     if v2(i) > v2(i-1) 
         v2(i) = v2(i-1);
     end;
end;

hline2= plot(v1,v2,'--r','LineWidth',1,...
    'MarkerEdgeColor','y',...
    'MarkerFaceColor','g',...
    'MarkerSize',1,'Marker','.');