showmoeazdt3.asv

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

function [f]=ShowMOEA(FileName)
% 将结果的点集 显示出来

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;

    gx=0;
    x1=0:0.001:1;
    x2=1:0.01:(1+29 *9 /29);
    
    for i=1:1000
        y(i) = 100;
        for j=1: length(x2)    
            d = x2(j) * ( 1- sqrt(x1(i)/x2(j)) -x1(i)*sin(10)
           if (y(i) > d ) y(i) =d;    
        end;
    end;
    
    
    
end;