examp081.m

来自「约束最优化」· M 代码 · 共 17 行

% max  z=3*x1+x2
% s.t. -x1+x2<=2
%      x1-2*x2<=2
%      3*x1+2*x2<=14
%      x1,x2>=0
c=-[3,1];a=[-1,1;1,-2;3,2];b=[2,2,14];
v1=[0 0];   
' ------  command one: x=linprog(c,a,b,[],[],v1) --------'
x=linprog(c,a,b,[],[],v1)
z=-c*x  %x= (4,1),max z= 13
pause
' ------  command two: [x,lag]=linprog(c,a,b,[],[],v1) --------'
[x,f,exitflag,output,lag]=linprog(c,a,b,[],[],v1)
z=-f
pause