shili082.m

来自「matlab讲义」· M 代码 · 共 15 行

% P225
r=-1.06;s=-1.15;
c=-[0,0,0,1.4,0,0,1.25,0,1.15,0,1.06];
A=[1,1,0,0,0,0,0,0,0,0,0;
0,r,1,1,1,0,0,0,0,0,0;
s,0,0,0,r,1,1,1,0,0,0;
0,0,s,0,0,0,0,r,1,1,0;
0,0,0,0,0,s,0,0,0,r,1;
0,0,0,1,0,0,0,0,0,0,0;
0,0,0,0,0,0,1,0,0,0,0];
b=[10,0,0,0,0,3,4];
[x,lag,how]=lp(c,A,b,zeros(length(c),1),[],[],5)
z=-c*x