这是我写的程序和结果,请老师看看到底错哪里了?为什么结果不满足初始条件,比如第一个解应该在444和1110之间,得到的结果却是2135.,超出很多?谢谢各位
c=1/1000*[0 1.725 1.725 4.45 4.65 4.65 4.875 4.95 5.025];
A=[-1 0 -1 0 0 0 0 0 0; 0 0 0 0 1 1 1 1 1;...
0 0 0 0 0 0 0 1 1; 0 0 0 0 0 0 0 -1 1];
b=[-3700 55500 40800 0]';
aeq=[0 0 0 -2.5054 -2.9167 -2.5057 -2.0256 -0.2321 1.3702;...
0 0 0 6.2788 8.5069 6.2801 4.1187 0.4602 3.5878;...
1 1 1 1 1 1 1 1 1];
beq=[-74797.5; 361331.3; 88400];
lb=[444;6290;0;0;0;0;0;0;0];
ub=[1110;88400;3700;88400;55500;55500;55500;40800;40800];
[x,y]=linprog(-c,A,b,aeq,beq,lb,ub,[],[])
Exiting: One or more of the residuals, duality gap, or total relative error
has grown 100000 times greater than its minimum value so far:
the primal appears to be infeasible (and the dual unbounded).
(The dual residual < TolFun=1.00e-008.)
x =
1.0e+004 *
0.213598551439313
0.715574818077885
0.358399133790258
2.087357433992631
1.741695885227843
0.000000000002390
0.000000000279880
1.939537309900007
2.197992515813596
y =
-3.988585359189311e+002