绘制已知散点数据的pareto前沿

        本代码用于美化图片，用软件Isight进行多目标优化后，软件自带的最大应力与质量的pareto图不美观，可导出数据后用MATLAB实现。

clear;clc;
data=xlsread('E:\Finalthesis\图\pareto.xlsx','D4:E9604')%导入数据
A=sortrows(data,[1,2])%从小到大排序
z=size(A,1);
for i=1:z-1
    for j=i+1:z
    if A(j,2)>A(i,2)&&A(j,1)~=A(i,1)
        A(j,2)=NaN;
    end
    end   
    A=A(all(~isnan(A),2),:);
    z=size(A,1);
end
scatter(data(:,1),data(:,2),5,[0 0 0],"filled")
hold on
scatter(A(:,1),A(:,2),20,[0 0 1])
hold on
scatter(135.68,252.01,30,[1 0 0],"filled")%初始值
xlabel("最大应力/Mpa")
ylabel("质量/kg")