帕累托最优前线

Question

我正在尝试获得两个健身功能的帕累托最优前沿。我通过使用在矩阵中为任何未指定的解决方案分配“1”的虚拟矩阵对未命名的解决方案进行了排序。当我绘制帕累托前线时，它会保留包含我知道不属于帕累托最优点的点。但是，我似乎无法找到这个问题的原因。任何帮助都会非常感激。

for  i = 1:1000
    f1(i) = x(i,1)^2;
    f2(i) = (x(i,1)-2)^2;
end
store = zeros(1000,1);
for i = 1:1000
    st = zeros(1000,1);
    for j = 1:1000
        if i == j
            st(j) = 1;                         
            continue;                           %Skip to next iteration.
        end
        if f1(i) > f1(j) && f2(i) > f2(j);      %Check for "x-dominated"
           continue;
        else st(j) = 1;                         %Dummy 1000x1 matrix
        end
    end
    if st == ones(1000,1)                       %Testing the dummy matrix for dominance
       store(i) = x(i);
    end
end

pareto = store(store ~= 0);                     
N = length(pareto);
for k = 1:N
    f3(k) = x(k,1)^2;
    f4(k) = (x(k,1)-2)^2;
end

Answer 1

不确定你做了什么，但这就是我用有限点绘制帕累托前线的方法。我认为这应该让你走上正轨：

t=1:10;
f1 = t.^2;
f2 = (t-2).^2;

ip = true(size(f1));

for k=1:numel(f1)
    if any(f1<f1(k)&(f2<f2(k)))
        ip(k) = false;
    end
end

plot(f1,f2)
hold all
plot(f1(ip),f2(ip),'ro')