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

relay1987

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

- _7 i7 S$ C9 A8 O; W4 [	甲 6 N: J; T( f: G# H0 t  F$ p	乙	丙 7 d1 \7 o1 r* G! ~5 r7 `  i6 L	丁
A   f) Z# w+ ]4 y  r' P	160	130	220 5 ]1 l8 M- @$ G3 E	170
B	140 $ [* ], C; G5 m6 w( d	130	190	150
C ; O' {. H% ?7 \2 \9 T* w& L	190	200	230	/ & M9 I6 m  L( r/ Y; Q9 g

: {  S+ T" _' L5 v* T  N% m3 V
7 \2 ^6 W: k# `+ u2 I4 {# `
; e  V4 |3 m" p5 f

         单位：元/千吨
" o7 }; j+ L2 q: g& @# {针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

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

mrx

【1】求解：
6 Z0 Z7 G9 h  c8 f) w设水库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
3 q1 M. F% |6 U% G6 ]st
9 F6 i+ l. I! M4 s* z+ fx34=0;
x11+x12+x13+x14<=50;
7 g* K) _( y3 v: wx21+x22+x23+x24<=60;
! f9 s. D" m* B  D5 f: |$ jx31+x32+x33+x34<=50;
7 n' q  I* B0 Kx11+x21+x31>=30;
! t: a; [. D6 G8 _x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
( j1 J" r: o5 Tx13+x23+x33>=10;
x13+x23+x33<=30;
( S, Q! w' C. Z. p5 i: nx14+x24+x34>=10;
x14+x24+x34<=50;
end
LINGO程序：
' Z' J4 X) h6 ]3 a' mmax=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;
x21+x22+x23+x24<=60;
8 ?& a; e* y: I8 bx31+x32+x33+x34<=50;
x11+x21+x31>=30;
* w% S2 S4 [9 e6 q1 _% Y# R+ _x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
" M; k, @; k9 j$ q. nx13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
x14+x24+x34<=50;
) w( `+ g# A, H, N  ^end
运行结果：
$ ^5 P* N' i. Q5 a- B  Global optimal solution found.
  Objective value:                              47600.00
0 L! _. j  @- p! H  Total solver iterations:                             7
7 q4 Z( w3 Z/ o2 V) X

2 S; J) J. a" P  ?Variable           Value        Reduced Cost
$ Z+ D5 M) l2 c! oX11        0.000000            30.00000
X12        50.00000            0.000000
X13        0.000000            50.00000
X14        0.000000            20.00000
X21        0.000000            10.00000
, {% `, y2 I- T. \2 d1 `X22        50.00000            0.000000
+ B* s8 }+ m/ oX23        0.000000            20.00000
) ]; |1 [' ~2 F0 hX24        10.00000            0.000000
3 l! p  E* E# A/ G- ^% sX31        40.00000            0.000000
" p! M/ j2 v5 ?4 @) f7 UX32        0.000000            10.00000
X33        10.00000            0.000000
6 E$ c/ ~/ a* UX34        0.000000            240.0000
- ]" k% s' r2 x
Row    Slack or Surplus      Dual Price
0 \8 R* e) c$ F4 v# M1        47600.00            1.000000
2        0.000000            320.0000
: k) r1 ~( P0 q& }0 }/ |- j" N" }3        0.000000            320.0000
4        0.000000            260.0000
( ^9 H+ r, n/ n5        10.00000            0.000000
6        40.00000            0.000000
  h  ^: _' |/ E. q' Z- U9 U9 y7        30.00000            0.000000
8        40.00000            0.000000
. h, ]" L( [# n9        0.000000           -40.00000
6 X$ D3 X) L" d7 u3 p10        20.00000            0.000000
11        0.000000           -20.00000
' b  _7 }# N  j4 l- S12        40.00000            0.000000

【2】求解：
2 d; }9 Q& ~% [3 u3 o- D' H3 _模型：
, \/ t- ]  B/ t% Vmax 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
st
9 Z6 [1 s' d* M. I2 q. a' ex34=0;
x11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
x11+x21+x31>=30;
x11+x21+x31<=80;
* w$ A( |; j' S9 \; r: n6 Vx12+x22+x32>=70;
% X) y$ H6 C6 Y. E( v- ^8 e& bx12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
1 p0 t4 V% f; A: nx14+x24+x34<=50;
end
4 @3 V& n9 R* i( y* ~LINGO程序：
. Q7 J1 E/ L: p; h! x! Ymax=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;
1 M5 H$ x9 F2 A* {x21+x22+x23+x24<=120;
, `; c4 h1 V: a2 e& Z& gx31+x32+x33+x34<=100;
8 A. r: b0 ^9 b! Y1 G, qx11+x21+x31>=30;
* [( e! B3 L6 Z/ Lx11+x21+x31<=80;
  H* }6 r% Z, ?$ u# `x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
* X! P# Y' ~+ R0 Dx14+x24+x34<=50;
) A8 G6 C- o- {  Hend
2 S, W7 R/ H$ b# }9 U运行结果：
  Global optimal solution found.
  Objective value:                              88700.00
  Total solver iterations:                             7
# ~4 G" o8 |; v2 y2 F/ `# v' m* G

Variable           Value        Reduced Cost
X11        0.000000            20.00000
! I8 M2 i" P  u: t% i# _% T" u( ~6 aX12        100.0000            0.000000
' s) r% X+ u: `' |: ~% n: V6 Z6 ?4 ]# H; QX13        0.000000            40.00000
. D3 G) L5 S. [) |2 H) q! qX14        0.000000            20.00000
: G$ m! u8 g+ x, E# GX21        30.00000            0.000000
# E, J* T8 {5 N6 s$ qX22        40.00000            0.000000
X23        0.000000            10.00000
" s3 X+ t7 k/ G" d. G+ c7 S7 f# MX24        50.00000            0.000000
2 D' H' b4 `4 m) t' S$ W+ `X31        50.00000            0.000000
% h" N) s, K+ I, f) \' i; L4 @; F8 oX32        0.000000            20.00000
X33        30.00000            0.000000
, p2 n/ Z$ q+ e0 YX34        0.000000            250.0000

Row    Slack or Surplus      Dual Price
1        88700.00            1.000000
( r' u2 @& q5 I" ?  B: D/ z2        0.000000            50.00000
/ D+ m2 w5 o' d& Q( F3        0.000000            50.00000
4        20.00000            0.000000
5        50.00000            0.000000
6        0.000000            260.0000
9 u: y  k; i$ w  ^2 f7        70.00000            0.000000
( F/ p, {- i: Y' }& ^8        0.000000            270.0000
$ I/ X* u) A$ p2 F9        20.00000            0.000000
: T$ q, \, s* Y10        0.000000            220.0000
11        40.00000            0.000000
12        0.000000            250.0000

alanloveai

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