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

relay1987

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

	甲	乙	丙	丁 1 T, g: ^: B- N1 T$ r- I  d7 j& ^, E
A $ v5 w9 `8 p' Z% F	160 ! |( ?# L4 O* t5 o; C	130	220	170 ' Z5 `2 K1 r5 g- m" O
B	140 4 ?, `/ M% J1 P5 B, }	130	190	150 & H8 r: t1 n! e/ J+ _! |1 P
C	190	200 ! I: k/ M+ Y% l* {$ I8 h	230 # G6 c1 \, W9 @4 i0 ~: j" o$ c% h	/ ; a- D" r. A" w

) A& [$ k( m' r+ h9 o! B
+ v5 p6 K0 f/ F8 j0 U
: z& {4 L5 ]- k2 L# O         单位：元/千吨
针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

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

mrx

【1】求解：
设水库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
* q: Y+ J  p2 ~: [* Z$ Vx34=0;
9 ]8 b3 y: R- Dx11+x12+x13+x14<=50;
! X3 K6 w& H! o2 ix21+x22+x23+x24<=60;
x31+x32+x33+x34<=50;
x11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
* ?6 K  A3 O" L( J* ~% z" o) c! ?x12+x22+x32<=140;
  Y# ?5 E1 N$ f- {2 nx13+x23+x33>=10;
5 u- B0 a& n: f2 P0 A- Ux13+x23+x33<=30;
# g7 ^6 _+ `/ Lx14+x24+x34>=10;
: `9 Q5 \( ^! R2 b- {$ c, H7 sx14+x24+x34<=50;
+ f3 i& V; ]* P2 h4 }end
LINGO程序：
" W2 E. l' d) ?! J. I( V4 E+ a, gmax=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
  \" z1 B, }( e! ]5 rx11+x12+x13+x14<=50;
% _* U, I1 J0 ?5 H1 n2 z& Mx21+x22+x23+x24<=60;
8 N# y) y! O' ~$ B4 b0 H2 ~x31+x32+x33+x34<=50;
  A* X/ I2 z0 `7 o2 Lx11+x21+x31>=30;
x11+x21+x31<=80;
) D" F# Q; |$ `7 M% ^- z0 ?x12+x22+x32>=70;
x12+x22+x32<=140;
# o# p1 b# v( j' z; cx13+x23+x33>=10;
x13+x23+x33<=30;
1 J1 y: ~* Y9 g1 @3 m0 E! Vx14+x24+x34>=10;
x14+x24+x34<=50;
end
运行结果：
  Global optimal solution found.
3 C( q0 D) L& `" S( j& F/ N- }- x  Objective value:                              47600.00
  Total solver iterations:                             7
- N- k# m/ w( P* L* I8 B9 ?

& t9 u$ t+ j$ D$ Y& ^  l5 o7 YVariable           Value        Reduced Cost
X11        0.000000            30.00000
X12        50.00000            0.000000
X13        0.000000            50.00000
4 m7 {: h  d" c8 cX14        0.000000            20.00000
X21        0.000000            10.00000
% L+ `# L6 y% i7 A( `; iX22        50.00000            0.000000
X23        0.000000            20.00000
! [- }, b- n* [2 q$ }; tX24        10.00000            0.000000
X31        40.00000            0.000000
X32        0.000000            10.00000
8 T7 ^9 c) z* M* nX33        10.00000            0.000000
2 s2 P3 U6 Q- A, mX34        0.000000            240.0000

( ~* ^# z; Q9 y- J2 T! wRow    Slack or Surplus      Dual Price
1        47600.00            1.000000
2        0.000000            320.0000
3        0.000000            320.0000
$ B2 s* Q2 B- {( V5 Z: y4        0.000000            260.0000
! t# Y! y8 ]" Z  S4 M& R- }; Z; z5        10.00000            0.000000
6        40.00000            0.000000
7        30.00000            0.000000
8        40.00000            0.000000
9        0.000000           -40.00000
* @* b8 T% d+ T$ l6 z9 q10        20.00000            0.000000
11        0.000000           -20.00000
  [. ^& k4 E9 b& j5 @6 q4 \: X12        40.00000            0.000000

3 c: _1 A/ Z1 z% T& a/ C【2】求解：
模型：
max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
st
  Z" {" ?8 E  O2 P% Dx34=0;
7 m; g- u0 G' o2 D# y+ ~1 P  Q2 ex11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
) n* F3 G. `$ ~9 y- |) b: Y6 Kx11+x21+x31>=30;
x11+x21+x31<=80;
" j" P. n$ t: E& C/ t- G: ^3 m2 lx12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
, i: I2 J9 d5 ~0 i( k5 Bx13+x23+x33<=30;
x14+x24+x34>=10;
x14+x24+x34<=50;
2 ~4 ~& \% ~. ~5 U% `7 pend
LINGO程序：
. E# J8 t, ~- f4 h( G" t+ ?* ~5 ?max=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
x11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
/ N' n; L/ D) _7 [! o6 sx31+x32+x33+x34<=100;
1 O- T% n0 @5 U! W) {8 ox11+x21+x31>=30;
9 n( t$ W$ q2 i0 H2 b7 Gx11+x21+x31<=80;
x12+x22+x32>=70;
" s) F! z1 C) d2 s4 U, _x12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
( u9 e2 ]2 Z8 N! |8 N0 Px14+x24+x34>=10;
0 U% E9 ~6 P" n2 Tx14+x24+x34<=50;
end
3 g; `+ [: Q; ~, J: U! ~6 K运行结果：
  Global optimal solution found.
9 t6 o0 E7 l8 s5 j  Objective value:                              88700.00
5 F0 h0 a8 r# ]  Total solver iterations:                             7
! Q3 j' Y- y' f: [: I5 J  }

0 C, R4 r6 h! M! ~9 D" M  @- {! L) DVariable           Value        Reduced Cost
X11        0.000000            20.00000
X12        100.0000            0.000000
4 a" z: R6 r" s/ i  E" i1 |X13        0.000000            40.00000
& Y' b( F; I5 N$ H) `X14        0.000000            20.00000
* ]. v& K! f  g" m- ~/ ]X21        30.00000            0.000000
( y- z$ B! b; N* rX22        40.00000            0.000000
X23        0.000000            10.00000
X24        50.00000            0.000000
X31        50.00000            0.000000
X32        0.000000            20.00000
- z8 ]: i+ G6 p. {1 \X33        30.00000            0.000000
/ M! B( Y, w2 }X34        0.000000            250.0000
* k' n/ \3 X5 s. ]. f$ e& q, _! y5 e
' Y0 r5 z8 Z" D% G: ORow    Slack or Surplus      Dual Price
$ Q6 I3 g: t# S) K1        88700.00            1.000000
# y, O# H: H* y7 j, p* n2        0.000000            50.00000
2 k9 J0 u) ~, x: a3        0.000000            50.00000
4        20.00000            0.000000
5        50.00000            0.000000
6        0.000000            260.0000
7        70.00000            0.000000
8        0.000000            270.0000
9        20.00000            0.000000
10        0.000000            220.0000
9 q2 u! E! s- ~) q4 U6 ^6 D11        40.00000            0.000000
12        0.000000            250.0000

alanloveai

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