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

relay1987

' }  P: o, ?7 W
有A、B、C三个水库往甲、乙、丙、丁四个小区供应自来水，其中C水库由于条件限制不能给丁小区供水，这三个水库每天的供水量分别为50，60，50（单位：千吨），这四个小区每天的基本用水量为30，70，10，10（单位：千吨），在保证基本用水量的前提下，每个小区还需要一定量的额外用水，这四个小区每天额外需求量分别为50，70，20，40（单位：千吨），小区只需按照900元/千吨的价格缴纳费用，水库往这四个小区的引水管理费如下表，且除引水管理费之外，每千吨水还有450元的其他费用。
- n8 o) F7 R  d0 t: ^( C" F. `* U

1 k0 b1 F. G. V2 x	甲	乙 " T/ ^7 ?; ]: R& V	丙	丁 9 n0 H* g6 \3 ~8 X% p5 z: {
A   q8 M& W* ]. q6 ~1 @	160 # G" X3 ?, O( y1 u# H! c	130 / J9 S& o7 h2 l6 |6 ]) Y	220	170 . m) l2 ?& R- N, c) A7 D- h
B + J+ q3 {4 c; l/ v	140 # D1 u; y6 P" [% W4 `9 n	130	190 / D% a" }1 t( _2 o" b	150
C	190 5 x% t+ v  v, I+ Y# y6 Q	200	230	/ , u* T3 ?8 _& Y* U$ ~

1 x1 N+ @2 G- t
" d, K8 i8 `' |. P


8 R) L; R( E5 j: {1 Z% m4 h. \         单位：元/千吨
+ F3 H0 n7 U: |: T& O2 V针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
$ A# x, ~8 k/ [1. 应如何分配水库供水量，公司才能获利最多？
+ D) e+ j6 j5 w. u( q( D6 N2. 若水库供水量都提高一倍，公司利润可增加到多少？

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
5 U2 Y# \2 J; w9 }0 M. S  est
8 w% b' r8 W+ l/ Dx34=0;
x11+x12+x13+x14<=50;
x21+x22+x23+x24<=60;
x31+x32+x33+x34<=50;
x11+x21+x31>=30;
2 M5 Y: M: V# yx11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
0 U( j( r! b+ e; mx13+x23+x33>=10;
7 x, _9 i5 ]/ l. nx13+x23+x33<=30;
( x8 ]7 [% Y8 ^7 j5 l( _# e  Bx14+x24+x34>=10;
9 l# K; F# K5 z# U/ y* r4 }) }x14+x24+x34<=50;
/ J& |2 f$ n+ d) W6 a8 Bend
LINGO程序：
4 b3 R" C2 ~; Xmax=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
2 b" V1 f/ w& Z, t: Nx11+x12+x13+x14<=50;
x21+x22+x23+x24<=60;
- l% C5 J- B2 o7 r% E( Kx31+x32+x33+x34<=50;
+ F" a& H$ }) H5 W3 h$ S; cx11+x21+x31>=30;
; G/ N, n0 A6 I) H# ?. kx11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
& e1 g/ u, t/ j6 _3 A5 Gx13+x23+x33>=10;
x13+x23+x33<=30;
) O5 m. b% {8 r. n5 z' |x14+x24+x34>=10;
; D) B! W: ]* S7 h6 |) ~; Vx14+x24+x34<=50;
9 u, S0 \) x/ A: rend
运行结果：
  Global optimal solution found.
  Objective value:                              47600.00
4 v- C. D6 \# {0 s. H8 U1 q  Total solver iterations:                             7
1 n2 v" e0 q# U' _, s+ |) x/ q

Variable           Value        Reduced Cost
X11        0.000000            30.00000
X12        50.00000            0.000000
$ o5 z% ?! o% N/ a2 o) vX13        0.000000            50.00000
X14        0.000000            20.00000
X21        0.000000            10.00000
. Y4 U- d+ o/ SX22        50.00000            0.000000
X23        0.000000            20.00000
8 d& v& n0 B2 a2 s) K* y+ f2 ~X24        10.00000            0.000000
X31        40.00000            0.000000
7 C- @# e& X- L6 ^X32        0.000000            10.00000
X33        10.00000            0.000000
X34        0.000000            240.0000

Row    Slack or Surplus      Dual Price
1        47600.00            1.000000
1 i1 R- J% ^% B- q/ N2        0.000000            320.0000
/ a% K% ^! ?" X' W3        0.000000            320.0000
4        0.000000            260.0000
5        10.00000            0.000000
- W% R" Z& f6 G6 w/ V0 F6        40.00000            0.000000
7        30.00000            0.000000
/ @3 H, A6 R, b2 X% K8        40.00000            0.000000
9        0.000000           -40.00000
10        20.00000            0.000000
% d# k2 l% h0 s% Q11        0.000000           -20.00000
12        40.00000            0.000000
, Z4 U, _, X# u  u  ?
5 N( D# E! R: h. t  l/ w【2】求解：
模型：
max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
st
. b1 f/ |! D8 q: Tx34=0;
x11+x12+x13+x14<=100;
8 v: _0 k: z! Z' |/ O) `% I6 [$ b" K, Ix21+x22+x23+x24<=120;
/ h6 u- L! M9 e# I5 qx31+x32+x33+x34<=100;
5 r/ X" |$ K3 \8 f; Y8 H- sx11+x21+x31>=30;
$ h' C# K& s9 I" L8 ^x11+x21+x31<=80;
- o7 d1 u- J2 B. f5 tx12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
* r+ Q! z! Z. G# k+ o; k7 Sx14+x24+x34>=10;
2 V3 r* E' [( X1 G9 w/ ox14+x24+x34<=50;
' O- a: F# {  b! g+ `0 Qend
LINGO程序：
9 P; ]' C, ~- u" K. U. b5 umax=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
1 {# r, x$ a6 U. E' V+ Wx11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
7 }" L; o% a' N- ^1 r0 d# ]- Fx11+x21+x31>=30;
8 L; X4 {( n* d! _/ X( l5 }x11+x21+x31<=80;
% C7 N) n2 L; _1 d" H9 lx12+x22+x32>=70;
6 Z; m6 Z  i9 f( i1 L4 O/ X8 }x12+x22+x32<=140;
7 V. L9 y- G, U- p+ G7 T: wx13+x23+x33>=10;
; E$ ^+ Z7 m+ k) \, lx13+x23+x33<=30;
- }7 ]; {6 b: [1 sx14+x24+x34>=10;
0 d4 s( w8 I) O5 c7 B. Z1 ax14+x24+x34<=50;
2 M, r6 x. f  V9 J) ^3 Xend
运行结果：
  Global optimal solution found.
- Y$ L3 ?# N( |' p+ m3 e2 X3 Y  Objective value:                              88700.00
9 A8 }$ q+ B5 r  Total solver iterations:                             7

" N$ k9 E7 ~" Z7 y5 m! J, ~( T
. O! q7 t+ y5 h2 o6 K. ]+ `' xVariable           Value        Reduced Cost
X11        0.000000            20.00000
X12        100.0000            0.000000
! B/ r7 ]$ y* M8 |2 `* ?8 LX13        0.000000            40.00000
7 O- _% z7 u9 F  m5 Z4 e6 V3 Q/ u  Q6 JX14        0.000000            20.00000
X21        30.00000            0.000000
% k$ n& F, s/ NX22        40.00000            0.000000
/ {3 f8 {  \: G% YX23        0.000000            10.00000
4 x! ^" }; s. M8 C1 }9 lX24        50.00000            0.000000
X31        50.00000            0.000000
; J/ B4 }& V7 O% j9 r2 a0 A& ~X32        0.000000            20.00000
7 y/ i* H: H8 \6 q% n1 ^# ?( WX33        30.00000            0.000000
' T/ W5 D" ^# OX34        0.000000            250.0000

' w# L. k; _1 |8 g+ T6 p  {) s. ARow    Slack or Surplus      Dual Price
1        88700.00            1.000000
2        0.000000            50.00000
3        0.000000            50.00000
4        20.00000            0.000000
4 w1 s& m8 z# m: r5        50.00000            0.000000
! f+ ]9 [7 `6 I6        0.000000            260.0000
7        70.00000            0.000000
8        0.000000            270.0000
! H! G% i9 a: o) R& E$ ~9        20.00000            0.000000
10        0.000000            220.0000
5 x  k0 R: W/ E0 v3 v% m11        40.00000            0.000000
12        0.000000            250.0000

alanloveai

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