急求该题答案，我试了10次也不会运行程序

relay1987

有A、B、C三个水库往甲、乙、丙、丁四个小区供应自来水，其中C水库由于条件限制不能给丁小区供水，这三个水库每天的供水量分别为50，60，50（单位：千吨），这四个小区每天的基本用水量为30，70，10，10（单位：千吨），在保证基本用水量的前提下，每个小区还需要一定量的额外用水，这四个小区每天额外需求量分别为50，70，20，40（单位：千吨），小区只需按照900元/千吨的价格缴纳费用，水库往这四个小区的引水管理费如下表，且除引水管理费之外，每千吨水还有450元的其他费用。
0 z/ V/ ~5 u! N! C

	甲	乙 9 A) }4 l' @5 T# R% I5 K	丙 / d6 C5 |$ X9 A- Z) K$ W	丁 # T: b* u1 ~% f$ ~
A	160	130 $ B1 g% N# e' V8 B	220	170 5 ]8 G* q: s- x8 i! s( p# d; E
B	140	130	190	150
C $ x$ i5 x6 d/ O' S2 n# e2 P- H	190 , m8 O" A, @# J7 A* I) h	200 ! I. K7 H1 q& R( ~  y8 J; i7 P6 z	230	/ 3 A; J, [. o9 u& X3 a' z' K6 X

1 u3 z( q" H; A. f3 q. I3 {& K

; R0 q! m8 d1 L7 t7 V
, r) R8 D( t" o2 `. G% Z7 |( K
+ O: {* y! `; \0 B, R# u         单位：元/千吨
针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
1. 应如何分配水库供水量，公司才能获利最多？
# t$ z9 \; X+ ~; M$ v2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

希望大家快些帮帮我，我参加学校比赛，老师病了，我不想打扰他，希望大家帮帮我，后天比赛，这些软件我学的不好，所以不怎么会用！今天看书写了个程序，可是不会运行，请高手们救救我吧，写好后发我的邮箱就行了relay1987@163.com，谢谢了！急急急，期待您的妙手回春！

mrx

【1】求解：
- V, E3 x8 Q: n9 g  y设水库A、B、C分别向甲、乙、丙、丁四个小区供水量为x11、x12、x13、x14、x21、x22、x23、x24、x31、x32、x33、x34，则有：
max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
st
2 u5 T4 ~7 e. r: o$ sx34=0;
x11+x12+x13+x14<=50;
$ K% f/ p1 ^$ t" X: _& A) d- Px21+x22+x23+x24<=60;
  q5 ]  z* d4 `3 Cx31+x32+x33+x34<=50;
x11+x21+x31>=30;
x11+x21+x31<=80;
7 h2 o/ _* R9 Q: Z. H6 Wx12+x22+x32>=70;
x12+x22+x32<=140;
7 |* p& g, g1 kx13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
+ P; R5 l1 m4 ]  m& ]* ]4 _x14+x24+x34<=50;
2 ~! b5 @& `8 h* i% Z9 a! P1 y3 oend
, o% F' i& U4 g3 n- GLINGO程序：
max=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
, G& z- C5 C, p6 Gx11+x12+x13+x14<=50;
( [9 Q3 Y0 Z; e* J8 mx21+x22+x23+x24<=60;
x31+x32+x33+x34<=50;
) w( \% |6 e( ]x11+x21+x31>=30;
5 \& p; r/ K/ C% k; s. z7 O2 \0 g$ [x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
4 K* @& W+ A" y! Jx13+x23+x33<=30;
* k' s/ R6 C  ]0 O3 ox14+x24+x34>=10;
2 C8 ^9 R/ ?0 I* rx14+x24+x34<=50;
end
运行结果：
  Global optimal solution found.
- ]. x4 e  t, B( X1 Q+ p  Objective value:                              47600.00
2 r! d% e2 e0 P0 t8 }' g7 G# N8 _; F% G  Total solver iterations:                             7

6 T2 t! ]; ]3 }" j* J
Variable           Value        Reduced Cost
X11        0.000000            30.00000
X12        50.00000            0.000000
X13        0.000000            50.00000
X14        0.000000            20.00000
X21        0.000000            10.00000
/ s+ n8 m( r1 r# \5 PX22        50.00000            0.000000
X23        0.000000            20.00000
X24        10.00000            0.000000
X31        40.00000            0.000000
; t7 z+ w- `" r# EX32        0.000000            10.00000
X33        10.00000            0.000000
X34        0.000000            240.0000

2 ]# }/ K/ q0 {" n* c" NRow    Slack or Surplus      Dual Price
6 C# v( j2 g6 k2 K2 m1        47600.00            1.000000
2        0.000000            320.0000
. ]4 _9 N9 [1 \! b3        0.000000            320.0000
4        0.000000            260.0000
5        10.00000            0.000000
9 O( C: I9 E5 i) n6        40.00000            0.000000
& I$ d: M% k4 a- S7        30.00000            0.000000
8        40.00000            0.000000
9        0.000000           -40.00000
. A) N5 R9 ?" E& `10        20.00000            0.000000
11        0.000000           -20.00000
12        40.00000            0.000000
5 H: w) X8 e7 k- Z) O* a
' m- d0 ~" Q7 G) s% t  D! m' u【2】求解：
% _( \( ~/ F6 F& {' P. W* ~模型：
6 D% i" W% L# w" q: imax 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
9 D+ h3 t) d% Y' g3 Tst
7 y, j6 K  I0 T$ O- U4 X* gx34=0;
x11+x12+x13+x14<=100;
  I* n# j7 I; _9 v* r: [3 u7 {+ zx21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
x11+x21+x31>=30;
" k" R2 x5 g$ Z2 Z4 ]2 @x11+x21+x31<=80;
. D+ O* D9 b" ^# [& Nx12+x22+x32>=70;
6 E" f& |5 I; X: Sx12+x22+x32<=140;
x13+x23+x33>=10;
& ~/ l. a6 a, b2 u! B& L7 p! C0 mx13+x23+x33<=30;
% |2 m# c) \( j/ Xx14+x24+x34>=10;
3 p/ h( `, k% Y; x6 h# nx14+x24+x34<=50;
( J; o& T# q4 F' [" M6 V8 z3 pend
' u; V; H9 x# ~+ n/ _LINGO程序：
max=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
; ^8 H7 J: R& h* d% Gx11+x12+x13+x14<=100;
7 A8 B1 s5 ~* l8 ]x21+x22+x23+x24<=120;
* m, y  R5 G2 r9 I  Px31+x32+x33+x34<=100;
5 A  l$ \$ ]7 c: Z6 ]x11+x21+x31>=30;
# H) L9 O; P$ M% s9 ]$ Hx11+x21+x31<=80;
1 m& F* S* _( ^7 D( R3 O0 \x12+x22+x32>=70;
x12+x22+x32<=140;
( B( @6 Y( E5 w% j) R- lx13+x23+x33>=10;
6 \/ j# T6 H4 r: L, P; o$ kx13+x23+x33<=30;
x14+x24+x34>=10;
x14+x24+x34<=50;
+ ]" ~$ R: ~- W( ?3 ?end
运行结果：
/ H$ Y& R$ `7 D- s# J  Global optimal solution found.
  Objective value:                              88700.00
  Total solver iterations:                             7

4 T! o( @0 g3 r5 E7 F5 @) P
' e: U4 p5 M- {4 i) b) MVariable           Value        Reduced Cost
" e  J# N. e8 u! |5 H8 v1 XX11        0.000000            20.00000
X12        100.0000            0.000000
" b. z. G. X+ UX13        0.000000            40.00000
X14        0.000000            20.00000
, k$ T! V, J$ @X21        30.00000            0.000000
# s8 t( }/ i  A# d( l9 ~+ t, w3 wX22        40.00000            0.000000
X23        0.000000            10.00000
X24        50.00000            0.000000
X31        50.00000            0.000000
X32        0.000000            20.00000
" O3 x# u- \5 r7 t. q) gX33        30.00000            0.000000
* |5 n8 U2 Y( I$ [9 f3 KX34        0.000000            250.0000

Row    Slack or Surplus      Dual Price
5 i" {& w2 Y+ J. a& O3 e1        88700.00            1.000000
  S% V" O* t0 B% x2 l$ m2        0.000000            50.00000
3        0.000000            50.00000
4        20.00000            0.000000
5        50.00000            0.000000
5 @- i% ]) }. F% f4 d9 z6        0.000000            260.0000
: m% h. z) @* w+ A, Q+ ~* ~7        70.00000            0.000000
8        0.000000            270.0000
# \3 d# C- X0 O' ]$ G3 L, E9        20.00000            0.000000
" [2 P6 ?- Z" J% B10        0.000000            220.0000
11        40.00000            0.000000
12        0.000000            250.0000

alanloveai

强············