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

relay1987

8 F" d  Q/ X' U' k. ?1 Y
( z; a) w2 `6 o0 a) K% t" S有A、B、C三个水库往甲、乙、丙、丁四个小区供应自来水，其中C水库由于条件限制不能给丁小区供水，这三个水库每天的供水量分别为50，60，50（单位：千吨），这四个小区每天的基本用水量为30，70，10，10（单位：千吨），在保证基本用水量的前提下，每个小区还需要一定量的额外用水，这四个小区每天额外需求量分别为50，70，20，40（单位：千吨），小区只需按照900元/千吨的价格缴纳费用，水库往这四个小区的引水管理费如下表，且除引水管理费之外，每千吨水还有450元的其他费用。

	甲 2 D9 R/ V4 P6 o& O$ R6 W8 X5 v! W5 w8 q	乙 9 f; w5 Y0 U: N. p# M	丙	丁
A 3 H5 k$ Q/ e4 ^8 I	160 3 X: z7 Q1 P' K" y) ~% y	130 & E# \3 C4 \; L$ ^/ Q/ X	220	170
B	140	130 & F; n& R2 D0 ]3 q	190 , J  H1 \+ b. v$ o	150 : J( x# g, n; @
C	190	200 : A4 p0 r& h, d+ R) x* x) U/ V$ E, @	230 ( A9 N7 k* _' O/ z4 D	/ 1 w/ R, O4 i. l2 r- j

) s$ ?4 I% F" V' W: o


4 G2 I0 {  q2 V! Y! C: q* {
         单位：元/千吨
针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

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

mrx

【1】求解：
/ n2 q1 W7 g# n& C0 S2 O( O设水库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
8 T" e: ^1 S& @& D$ mst
x34=0;
0 P5 a* y  M' V; y  L% yx11+x12+x13+x14<=50;
/ U# E9 F% |0 W: |' Ex21+x22+x23+x24<=60;
3 w, }$ i' ~8 e2 B$ g* {" j: }9 ^x31+x32+x33+x34<=50;
x11+x21+x31>=30;
x11+x21+x31<=80;
$ x& Q- v1 I2 ~; M4 ^* M6 o/ Tx12+x22+x32>=70;
7 J. M6 R* d9 v; b9 K  a) kx12+x22+x32<=140;
x13+x23+x33>=10;
" {4 q5 k5 [" a, J& r' J4 Gx13+x23+x33<=30;
7 O$ b4 b8 B0 t' G2 q! @) A+ [x14+x24+x34>=10;
7 h$ @/ V( n- a* {8 }9 Rx14+x24+x34<=50;
9 e; ]5 p6 _) T* D: V2 V) y" {end
0 S1 I2 E: ?( B' uLINGO程序：
/ Z8 _3 ?, B  ~6 l9 @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<=50;
8 Y3 K$ z8 C. s+ J" Sx21+x22+x23+x24<=60;
# e+ k7 z* l7 S+ Z+ c5 g; `0 `' vx31+x32+x33+x34<=50;
x11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
; K2 q: i; @/ T! G9 {& L) ?7 f1 r) mx13+x23+x33>=10;
! Y  G1 q+ w, U- r5 \x13+x23+x33<=30;
- b/ t- K# k: c1 ^* c' V( d- ?x14+x24+x34>=10;
x14+x24+x34<=50;
, Q7 {/ ~$ O" Xend
/ D1 f3 x' b3 U运行结果：
  Global optimal solution found.
  Objective value:                              47600.00
% G) d1 K, T( j$ O0 }  Total solver iterations:                             7

9 p7 X5 T* G" X+ a+ Q
Variable           Value        Reduced Cost
8 }/ }; g5 _# _/ c& F! yX11        0.000000            30.00000
X12        50.00000            0.000000
X13        0.000000            50.00000
X14        0.000000            20.00000
X21        0.000000            10.00000
" q/ n" `, `/ B7 a$ l. @! t' P& BX22        50.00000            0.000000
' N& w3 h( R1 E4 y: e% D% ZX23        0.000000            20.00000
X24        10.00000            0.000000
* ^0 N& H' W: LX31        40.00000            0.000000
% v% s- C* I8 c  rX32        0.000000            10.00000
X33        10.00000            0.000000
, h$ D, D0 j) zX34        0.000000            240.0000

3 n$ F# e4 ]7 W5 }Row    Slack or Surplus      Dual Price
# n+ k1 }- z; i9 B" O  m) {1        47600.00            1.000000
2        0.000000            320.0000
- }: M/ w8 C# @1 i+ v3        0.000000            320.0000
! G% S" U5 E% Z7 Z- w4        0.000000            260.0000
0 G/ E( {# I3 d* D+ G. ]1 U1 |, F5        10.00000            0.000000
; @+ R& `7 k3 M0 H! Q7 w# j6        40.00000            0.000000
5 R1 f1 a! ]8 S  U5 \3 w9 P! ?7        30.00000            0.000000
8        40.00000            0.000000
. x" {$ j  O9 I9        0.000000           -40.00000
) n$ ]3 J6 M8 l1 q8 y! A5 r10        20.00000            0.000000
11        0.000000           -20.00000
12        40.00000            0.000000
$ M; M) u' E# X, `: Q, B( o# p/ Q
7 `% i& `& p/ N4 p【2】求解：
模型：
" R2 v4 Y! V) x, Zmax 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
! B  n$ }) Q+ J3 jst
x34=0;
5 i9 s0 D  r+ e; J; m  Rx11+x12+x13+x14<=100;
" D4 C! u0 Q& _. x! W0 e$ q5 e3 t8 Nx21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
& k1 j7 w( p" X1 Hx11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
  \- c* E! \- N- rx13+x23+x33<=30;
8 I; [8 I1 h/ H/ lx14+x24+x34>=10;
x14+x24+x34<=50;
end
LINGO程序：
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;
x31+x32+x33+x34<=100;
) Q9 H/ w1 }0 T% ex11+x21+x31>=30;
: L. Y) v7 _0 [* c3 o& Px11+x21+x31<=80;
0 ]1 {6 t. p2 z6 T  L; g5 Ix12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
, N1 R- f. V$ m9 t; dx13+x23+x33<=30;
8 j% ], L4 q2 ]2 ^! U9 gx14+x24+x34>=10;
x14+x24+x34<=50;
end
运行结果：
$ J4 O) ]  ?% {' `- Q: W  Global optimal solution found.
  Objective value:                              88700.00
  Total solver iterations:                             7
, t- N' f( ^/ \% x' f8 \& @% e# b

Variable           Value        Reduced Cost
0 O) q. u/ s8 e' |X11        0.000000            20.00000
X12        100.0000            0.000000
X13        0.000000            40.00000
X14        0.000000            20.00000
/ I6 n; B8 j" u3 eX21        30.00000            0.000000
: N, C6 p; D+ U) DX22        40.00000            0.000000
X23        0.000000            10.00000
6 O1 R5 W7 E3 Z4 }/ qX24        50.00000            0.000000
% w2 U0 s2 {+ s  t5 L/ J0 U1 a6 H2 GX31        50.00000            0.000000
X32        0.000000            20.00000
X33        30.00000            0.000000
' A7 o  }/ I* i3 u& B" `X34        0.000000            250.0000

1 |: d7 V$ ^0 `1 Q1 V' SRow    Slack or Surplus      Dual Price
1        88700.00            1.000000
2        0.000000            50.00000
3        0.000000            50.00000
$ D+ x& s# ~+ B) e4 t  C% |4        20.00000            0.000000
( H# h! l7 t8 j8 V& ~5        50.00000            0.000000
: g. f0 K+ d) C; `  W0 l* ^6        0.000000            260.0000
# }# Y. ^) l5 s" F. W7        70.00000            0.000000
8        0.000000            270.0000
, |* Z9 x4 k  b3 u, s" c/ U! z9        20.00000            0.000000
10        0.000000            220.0000
11        40.00000            0.000000
, U% K" ]' g  ]3 N12        0.000000            250.0000

alanloveai

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