matlab线性优化

来源：百度知道 编辑：UC知道 时间：2024/05/15 01:24:33
已经编了一个，可是答案有点特殊得让人觉得不对
求解最大值Z=y1×5.2×105+(2.5+x1+x2-y1) ×4.3×105+y2×4.7 × 105+(8.2—x1+x3-y2) ×3.9×105+[0.25×y1+0.2×(2.5+x1+x2-y1)+0.23×y2+0.185×(8.2-x1+x2-y2)-2]×106-x1×2.0×105-x3×8.5×105-x2×10×105
约束条件：3.5-x2-x3≥0
8.2-x1+x3≥0
0≤y1≤2.5+x1+x2
0≤y2≤8.2-x1+x3
7.5y1+6.1 ×(2.5+x1+x2-y1)+9.0y2+7.35×(8.2-x1+x3-y2) ≤96.5
1.4y1+1.65y2≤7.5
0.25y1+0.2×(2.5+x1+x2-y1)+0.23y2+0.185×(8.2-x1+x3-y2) ≥2
x1≥0,x2≥0,x3≥0
谢谢了

clc;clear;
fmax=[-166.41 -557.69 -483 99.8 88.77];
fmin=-fmax;
A=[0 1 1 0 0;1 0 -1 0 0;-1 -1 0 1 0;1 0 -1 0 1;6.1-7.35 6.1 7.35 7.5-6.1 9.0-7.35;0 0 0 1.4 1.65;-0.2+0.185 -0.2 -0.185 -0.25+0.2 -0.23+0.185]
B=[3.5;8.2;2.5;8.2;96.5-75.52;7.5;-2+2.017]
LB=[0 0 0 0 0]%下限
[X,FVAL,EXITFLAG]=linprog(fmin,A,B,[],[],LB,[])

%关键是化简容易出错，化简过程如下：
%{
syms x1 x2 x3 y1 y2
f1=y1*5.2*105+(2.5+x1+x2-y1)*4.3*105+y2*4.7*105+(8.2-x1+x3-y2)*3.9*105+(0.25*y1+0.2*(2.5+x1+x2-y1)+0.23*y2+0.185*(8.2-x1+x2-y2)-2)*106-x1*2.0*105-x3*8.5*105-x2*10*105

f1 =
499/5*y1+1122113/250-16641/100*x1-55769/100*x2+8877/100*y2-483*x3

f2=vpa(f1,8)
f2 =
99.800000*y1+4488.4520-166.41000*x1-557.69000*x2+88.770000*y2-483.*x3
%}

%{
3.5-x2-x3≥0
8.2-x1+x3≥0
0≤y1≤2.5+x1+x2
0≤y2≤8.2-x1+x3
7.5y1+6.1 ×(2.5+x1+x2-y1)+9.0y2+7.35×(8.2-x1+x3-y2) ≤96.5

syms x1 x2 x3 y1 y2
ff=7.5*y1+6.1 *(2.5