pg573.m

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

s=0.92;
while s<=0.924
    c=[0.05*(s-1) 0.27*(s-1) 0.19*(s-1) 0.185*(s-1) 0.185*(s-1) s];
    A=[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 0-1];
   b=[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=[];
   s
   [x,val]=linprog(c,A,b,Aeq,beq,vlb,vub)
   plot(s,val,'.');
   xlabel('偏好系数'),ylabel('投资组合')
   hold on;
   s=s+0.00001;
   pause
end