pg572.m

包括了高教出版社出版的数学建模教材里部分习题的答案

k=0.19;
while k<0.23
    c=[0 0 0 0 0 1];
    A=[-0.05 -0.27 -0.19 -0.185 -0.185 0
        0 0.025 0 0 0 -1
        0 0 0.015 0 0 -1
        0 0 0 0.055 0 -1
        0 0 0 0 0.026 -1];
    b=[-1*k;0;0;0;0];
    Aeq=[1 1.01 1.02 1.045 1.065 0];
    beq=[1];
    vlb=[0 0 0 0 0 0];
    vub=[];
    [x,val]=linprog(c,A,b,Aeq,beq,vlb,vub)
    plot(k,val,'.');
    xlabel('收益'),ylabel('风险');
    %axis([0.1 0.5 0 0.5]);
    hold on
    k=k+0.001;
    pause
end