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轴标签