MCM线性规划（二）

Matlab code:

---

clc,clear
a=0;hold on     %保持原图并接受此后绘制的新的曲线
while a<0.05
    f=[-0.05,-0.27,-0.19,-0.185,-0.185];
    A=[0 0.025 0 0 0 ;0 0 0.015 0 0 ;0 0 0 0.055 0;0 0 0 0 0.026];
    b=a*ones(4,1);
    aeq=[1 1.01 1.02 1.045 1.046];
    beq=1;
    lb=zeros(5,1);
    [x,y]=linprog(f,A,b,aeq,beq,lb,[]);
    y=-y;plot(a,y,'*k');%第三个参数是LineSpec曲线属性,见下图
    a=a+0.001;
end
xlabel('a'),ylabel('y')%添加x轴，y轴标签