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

relay1987

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

$ S" h4 R" X% X1 }3 k2 k	甲	乙	丙	丁 0 M0 B4 ]! B2 j! A5 c/ Z
A	160	130 ) s$ c; A7 {) a6 c/ J	220 % Y2 f1 m  A* u	170 $ P2 F7 k7 t3 x- b. G
B ; x) H. ]) J8 K# U	140   I1 e1 {% A1 n6 Z+ Z5 u- j	130	190 ; a# [# A" T& X; l	150 + N( Q$ p# j* n6 }3 A- o
C	190	200	230 # k0 m; t3 u; W3 l+ F: Y+ O& W% O2 Y	/ + T0 `- G  A* f% u$ }3 m  h; P

& T# M5 B. S8 I) b* q. g4 n
! b3 E# g7 {7 t4 F
9 Y+ }& C5 l2 b0 q         单位：元/千吨
7 d  u. u) K3 O' h+ t+ `. W! Q% q针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

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

mrx

【1】求解：
1 [4 ~; z2 x/ |' ^设水库A、B、C分别向甲、乙、丙、丁四个小区供水量为x11、x12、x13、x14、x21、x22、x23、x24、x31、x32、x33、x34，则有：
2 A6 g6 j% Z' [max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
7 x2 H' f  ~$ ~  B' H* r; gst
$ b7 U" T; g7 a% |3 D; Xx34=0;
x11+x12+x13+x14<=50;
& l/ X6 i+ V! ]) p3 z5 I  P- tx21+x22+x23+x24<=60;
x31+x32+x33+x34<=50;
x11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
$ F  [, v' a' n1 C% tx13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
" Y- N; _7 }( z7 z7 E+ B$ Lx14+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;
% E. h: y* ?- M, W, s0 Yx11+x12+x13+x14<=50;
x21+x22+x23+x24<=60;
x31+x32+x33+x34<=50;
x11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
/ G$ C3 B1 }. r/ L6 J& d7 Ux12+x22+x32<=140;
" t/ v: w% Q5 ax13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
) u3 R' a( Z( l: {* vx14+x24+x34<=50;
end
! B/ N  ^8 W9 q运行结果：
  Global optimal solution found.
  Objective value:                              47600.00
  Total solver iterations:                             7

5 Y% W) Z# R+ r! x" F7 m
; b) C% ~; X- rVariable           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
% x) O9 H3 o* b4 r% |! Y% yX22        50.00000            0.000000
: u$ c5 ?2 d9 s! P# v4 oX23        0.000000            20.00000
' ]" A/ ~" p1 s/ h/ `, b* O" d* oX24        10.00000            0.000000
X31        40.00000            0.000000
X32        0.000000            10.00000
- `+ u5 |( h  x1 m: H, Z( DX33        10.00000            0.000000
X34        0.000000            240.0000

0 [2 T9 F; _6 z9 ~' ORow    Slack or Surplus      Dual Price
6 T& z4 H  w1 t- y7 H+ o1        47600.00            1.000000
8 `( Y4 A' `9 h7 k5 ?2        0.000000            320.0000
! K( d& i* p+ e0 {5 h* b4 l; C3        0.000000            320.0000
# w6 ^# Y( p5 p4 T4 l& ]% r4        0.000000            260.0000
5        10.00000            0.000000
# R, O7 R" |% K% w+ b6        40.00000            0.000000
8 a- m- u$ N8 ~3 F7        30.00000            0.000000
% Q& v0 s7 z/ j  a1 s; [4 u8        40.00000            0.000000
9        0.000000           -40.00000
4 h8 a# C) y$ l" f$ r, \: o10        20.00000            0.000000
11        0.000000           -20.00000
12        40.00000            0.000000
) o, W/ Q1 W( F7 G( j
  n$ b& _4 ?/ `+ R2 S$ t7 s【2】求解：
模型：
max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
st
. U  X, y. v( P+ }* @2 ?! ex34=0;
x11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
; P9 k; r; Y! P0 j! q6 c: {. B4 T' bx11+x21+x31>=30;
2 T7 r" S- o, D( dx11+x21+x31<=80;
/ X! q: D; t# I4 ax12+x22+x32>=70;
4 m  D" _9 f  j, |; ]9 U% H( Zx12+x22+x32<=140;
x13+x23+x33>=10;
  ~0 T$ j! i9 Kx13+x23+x33<=30;
x14+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;
8 W+ ~/ L; ?7 V' e. b' Rx11+x12+x13+x14<=100;
6 O' A8 b( j, Ax21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
8 e6 U- i& q* T1 a8 dx11+x21+x31>=30;
( I8 |& w; k, L! @# }x11+x21+x31<=80;
0 l9 s  h. \6 P% A5 ~x12+x22+x32>=70;
x12+x22+x32<=140;
  t, `& [$ h# M3 Px13+x23+x33>=10;
. C* E& M) k3 w, h8 vx13+x23+x33<=30;
+ Z( T1 k0 ]# m. y, Z/ o0 B& a& e" ]x14+x24+x34>=10;
x14+x24+x34<=50;
end
0 J# q5 x; \. Z运行结果：
/ ?! H% x5 E$ y+ |* \% M: M  Global optimal solution found.
2 {& L8 A* z# b$ N+ G6 o  Objective value:                              88700.00
  Total solver iterations:                             7
  i* `4 Q) H6 e) ~) ~2 i% d
# a0 Y5 Q- u: A
0 t2 T) y$ o) y9 j- f$ eVariable           Value        Reduced Cost
0 x3 [% l: h# m8 D' c; u- C0 eX11        0.000000            20.00000
8 S+ K  C8 R8 H, d# |+ m$ SX12        100.0000            0.000000
" G, ~- U6 Z. O4 G, `! T0 DX13        0.000000            40.00000
X14        0.000000            20.00000
2 B3 a  q1 O  _0 s/ g6 X  xX21        30.00000            0.000000
X22        40.00000            0.000000
X23        0.000000            10.00000
! L9 \, m% r- g( Y7 n1 r) x! D1 Q) U, y+ qX24        50.00000            0.000000
X31        50.00000            0.000000
) S+ {1 k& |( n, GX32        0.000000            20.00000
* V3 S; N% p, @. S  [X33        30.00000            0.000000
X34        0.000000            250.0000

2 e7 C: d: A) E3 Y0 L5 L) i  BRow    Slack or Surplus      Dual Price
1        88700.00            1.000000
! Z. c. \, t1 n7 G, n7 t2        0.000000            50.00000
3        0.000000            50.00000
9 h& M$ H7 |* A4 d6 b. a4        20.00000            0.000000
4 I, O% t5 N0 M) C8 T+ W5        50.00000            0.000000
6        0.000000            260.0000
9 K5 H! h; n: }' {7        70.00000            0.000000
' G/ F/ @0 E' P, [% j' M8        0.000000            270.0000
$ [9 n2 D$ D$ Y% c; `2 T" D) e) G9        20.00000            0.000000
10        0.000000            220.0000
11        40.00000            0.000000
12        0.000000            250.0000

alanloveai

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