showpol.m

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

function [f]=ShowPOL()
% 画函数的点
v1=0:0.001:25;
k1=-pi:0.001:pi;
k2=-pi:0.001:pi;
for i=1:length(v1)
    v2(i)=0;
end;
for i=1:length(k1)
    for j=1:length(k2)
        a1 = 0.5 * sin(1)     - 2* cos(1)     +    sin(2)   - 1.5 *cos(2);
        b1 = 0.5 * sin(k1(i)) - 2* cos(k1(i)) + sin(k2(j))  - 1.5 *cos(k2(j));
        a2 = 1.5 * sin(1)     -    cos(1)   + 2* sin(2)     - 0.5 *cos(2);
        b2 = 1.5 * sin(k1(i)) -  cos(k1(i)) + 2*sin(k2(j)) - 0.5 *cos(k2(j));
      
        dd =     1 + (a1-b1)*(a1-b1) +(a2-b2)*(a2-b2);    
        if (dd <=0 ) || (dd >=25) 
            continue;
        end;
 
        td = (k1(i) +3)*(k1(i) +3) + (k2(j) +1)*(k2(j) +1);
%        if floor(dd*1000) +1 >length(v2)
%           disp(floor(dd*1000) +1);
%           disp(length(v1));
%        end;
        if v2(floor(dd*1000) +1) > td || v2(floor(dd*1000) +1)==0
            v2(floor(dd*1000) +1) = td;
        end;
    end;
end;

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