showmoeazdt3.m

基于GPS的多目标优化及动态多目标优化源代码 里面 包含MOP, DMOP的程序

function [f]=ShowMOEA(FileName)
% 画函数的点
    
    z1=0:0.001:1;
    z2=1:0.01:(1+29 *9 /29);
    
    for i=1:length(z1)
        y(i) = 100;
        for j=1: length(z2)    
            d = z2(j) * ( 1- sqrt(z1(i)/z2(j)) -z1(i)*sin(10*pi*z1(i))/z2(j));
            if (y(i) > d ) y(i) =d;    
            end;
        end;
    end;    
    hline2= plot(z1,y,'--r','LineWidth',1,...
            'MarkerEdgeColor','y',...
            'MarkerFaceColor','g',...
            'MarkerSize',1,'Marker','.');
    

hold on;
% 将结果的点集 显示出来

for fcount=1:length(FileName)
    t =char(FileName(fcount));
    fid = fopen(t,'rt');
    if fid ==-1 
        return ;
    end;
    s = fscanf(fid,'%g %g',[2,inf]);
    fclose(fid);
    s=s';
    k=1;
    for i=1:length(s)
        flag =0;
        for j=1:length(s)
            if  ((i~=j) &&(s(i,1) >=s(j,1) ) && (s(i,2) >=s(j,2))  )
                flag =1;
                break;
            end;
        end;
        if (flag==1) 
            continue;
        end;
        f(k,1) = s(i,1);
        f(k,2) = s(i,2);
        k=k+1;
    end;
    
    for i=1:length(f)
        x1(i)= f(i,1);
        y1(i)= f(i,2);
    end;
    % 开始画点
    if (fcount ==1)
        hline2= plot(x1,y1,'--r','LineWidth',1,...
            'MarkerEdgeColor','k',...
            'MarkerFaceColor','g',...
            'MarkerSize',2,'Marker','s');
    end;
    if (fcount ==2)
        hline2= plot(x1,y1,'--r','LineWidth',1,...
            'MarkerEdgeColor','b',...
            'MarkerFaceColor','g',...
            'MarkerSize',2,'Marker','+');
    end;
    if (fcount ==3)
        hline2= plot(x1,y1,'--r','LineWidth',1,...
            'MarkerEdgeColor','g',...
            'MarkerFaceColor','g',...
            'MarkerSize',2,'Marker','.');
    end;
    if (fcount ==4)
        hline2= plot(x1,y1,'--r','LineWidth',1,...
            'MarkerEdgeColor','y',...
            'MarkerFaceColor','g',...
            'MarkerSize',2,'Marker','d');
    end;
end;