showsph150.m

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

function [f]=ShowSPH150()
% 画函数的点
v1=0:0.001:2.0;
k1=-0:0.001:1;
k2=-0:0.001:1;
k3=0:0.001:0.001;

v2(length(v1))=0;
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 =( (k1(i) - 1) *  (k1(i) - 1) +(k2(j)) * (k2(j)) +k3(k)); 
            td =k1(i) * k1(i)+ (k2(j) - 1) * (k2(j) - 1) +k3(k); 
            if dd>=2.0 
                continue;
            end;

%               disp(floor(dd*1000) +1);
%               disp(length(v1));
%               disp(dd);
            if (floor(dd*1000) +1 > length(v2))
                continue;
            end;

            if v2(floor(dd*1000) +1) > td || v2(floor(dd*1000) +1)==0
                v2(floor(dd*1000) +1) = td;
            end;
        end;
    end;
end;

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