pg571.m

来自「包括了高教出版社出版的数学建模教材里部分习题的答案」· M 代码 · 共 19 行

a=0;
while a<0.1
    c=[-0.05 -0.27 -0.19 -0.185 -0.185];
    Aeq=[1 1.01 1.02 1.045 1.065];
    beq=[1];
    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;a;a;a];
    vlb=[0 0 0 0 0];
    vub=[];
    [x,val]=linprog(c,A,b,Aeq,beq,vlb,vub);
    a
    Q=-val
    plot(a,Q,'.')
    xlabel('风险'),ylabel('收益')
    axis([0 0.1 0 0.5])
    hold on
    a=a+0.001;
    pause
end