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

relay1987

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

' w/ ?6 _  @% y9 I/ V	甲 ! F' i( W* @! o6 a! Y; u) Z, G	乙 & o/ o6 V+ s1 s; O, O+ P	丙	丁
A	160	130 ! x8 X- i1 N' z! ^% G	220	170 3 Y5 G7 F, _1 l
B 7 Y& d% V8 ^3 P' M0 j	140 ) V- n. V2 y" d' }	130	190 8 R1 v" J% G) S6 l+ i. K	150 + i4 p7 h! j, C5 s, I* E4 L7 u
C & z  W8 A9 \1 ^" R. c0 P, [	190 4 j3 c- j+ J% c0 u	200	230 ; Q7 u. R6 d9 }' a3 W7 l	/

! D8 [* D$ g$ G1 w$ b- l) P
1 F- T  N4 n2 t8 w2 Z+ r

# p; H  Y0 s/ M' t  G! g# x0 t         单位：元/千吨
针对下面问题建立合适的数学模型，并写出程序（用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
" [. |" N2 K' Z- Zst
8 G' [  A6 [, G- Y" Tx34=0;
x11+x12+x13+x14<=50;
x21+x22+x23+x24<=60;
7 Z9 C5 B3 {% T, Z3 C% E$ X% ax31+x32+x33+x34<=50;
# P1 w. }7 [, ix11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
$ l: u6 _6 r" O# i" p3 c# tx13+x23+x33<=30;
x14+x24+x34>=10;
0 Z5 T( q# D2 _x14+x24+x34<=50;
4 x. V& B$ a% E- p% z% O0 Fend
3 q& e: t  n2 uLINGO程序：
0 ?4 l2 K, S2 E& J) p: Z0 Z* \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;
6 H5 b6 Q' @2 O* v' h2 O4 B. Nx21+x22+x23+x24<=60;
* O& A" \& R5 y1 Lx31+x32+x33+x34<=50;
  {, C4 r# Y' i3 w- s6 vx11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
3 C0 ~6 o& ^7 D$ r6 Wx13+x23+x33>=10;
9 ]% \7 I) ]6 ^x13+x23+x33<=30;
x14+x24+x34>=10;
6 h3 y0 V9 s* X( }' }- lx14+x24+x34<=50;
; b+ L( ?; X+ L2 J+ F, i# h) eend
* j* ]; f- X, `: X运行结果：
7 O% S( j1 e8 g; |# E  Global optimal solution found.
  Objective value:                              47600.00
  Total solver iterations:                             7

  u! U7 p. z+ Z8 i: S- ^/ b
Variable           Value        Reduced Cost
X11        0.000000            30.00000
X12        50.00000            0.000000
3 ^1 ^1 E: w! x0 u' [, }( |X13        0.000000            50.00000
" Y) ~& Q$ p7 r6 a2 QX14        0.000000            20.00000
X21        0.000000            10.00000
7 E8 q8 ?& H- p* ]# vX22        50.00000            0.000000
X23        0.000000            20.00000
X24        10.00000            0.000000
X31        40.00000            0.000000
, Z  N# i' ?6 e' `6 V! X5 }X32        0.000000            10.00000
X33        10.00000            0.000000
  l5 I7 _" `8 N3 z8 h7 wX34        0.000000            240.0000
; x0 S( C4 F/ F
' f& Y, ]) d' ARow    Slack or Surplus      Dual Price
1        47600.00            1.000000
2        0.000000            320.0000
: U) v8 g; ]4 j% K3 H5 A6 y! b3        0.000000            320.0000
% |& @7 w" W; W- ?6 K9 P2 P4        0.000000            260.0000
% p- N. @  i5 t/ H8 o3 c5        10.00000            0.000000
6        40.00000            0.000000
# v# p( u  y& p5 }- X7        30.00000            0.000000
8        40.00000            0.000000
. L: Q& A  o' }9        0.000000           -40.00000
10        20.00000            0.000000
11        0.000000           -20.00000
$ `: r1 u2 F' Q/ ^' b  @5 _12        40.00000            0.000000
4 H0 V( D  m% i7 c! J- X$ Z& m- h
【2】求解：
3 R/ @1 g+ ^6 ~$ K8 X* ~1 u模型：
max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
3 T8 z% b0 W. ~9 M$ nst
x34=0;
9 o5 a' E4 M5 P. vx11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
: b# F( I0 |: o7 A( Nx11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
; Q% n% v3 V6 y) q/ cx13+x23+x33>=10;
# Q% M0 X, e( {2 j; T# px13+x23+x33<=30;
, f! Y3 Q9 A& K6 }) ?. h8 v1 lx14+x24+x34>=10;
x14+x24+x34<=50;
end
LINGO程序：
  {/ y" D6 ^4 F% k0 \" Dmax=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;
$ |: Z* f5 G4 E5 T. A3 a4 ^7 }x21+x22+x23+x24<=120;
! Z# [5 o- T, Y# vx31+x32+x33+x34<=100;
; s5 q- d3 E( Y$ |. D# B# ux11+x21+x31>=30;
. H: m, A% c# |$ X8 nx11+x21+x31<=80;
x12+x22+x32>=70;
  y2 \; ^, _: d4 \x12+x22+x32<=140;
8 t# C  P. j8 cx13+x23+x33>=10;
" `" C/ J4 ^+ d. P0 {, Q2 fx13+x23+x33<=30;
0 s' z3 ^$ `: {! w9 Tx14+x24+x34>=10;
x14+x24+x34<=50;
end
  q7 i: ]( [* h( f$ `) L4 t4 O运行结果：
6 `7 _8 |" a9 x( c0 q6 }' G/ s+ i  Global optimal solution found.
% j; B9 t4 R, C; n  Objective value:                              88700.00
' C0 l$ z1 i, t+ O2 @" O2 o  Total solver iterations:                             7


$ j7 p' N% D! _7 X2 XVariable           Value        Reduced Cost
* f% @) R4 V! g; Z" l' h4 C: |X11        0.000000            20.00000
X12        100.0000            0.000000
X13        0.000000            40.00000
X14        0.000000            20.00000
) d9 }$ d& O4 W/ R9 h* P4 h  K4 }X21        30.00000            0.000000
X22        40.00000            0.000000
X23        0.000000            10.00000
* K* `/ S, n- Q! C+ ]# jX24        50.00000            0.000000
- t$ G9 _8 |6 C; Y0 A7 EX31        50.00000            0.000000
$ X2 D0 N1 g! Q' r  DX32        0.000000            20.00000
X33        30.00000            0.000000
6 M* X, V  ]# q+ H8 N0 m  h7 d( ]& pX34        0.000000            250.0000

Row    Slack or Surplus      Dual Price
/ y9 u5 u  C3 b" L0 z2 `# H1        88700.00            1.000000
! y8 J8 X8 J9 h& D/ \  \& V. D2        0.000000            50.00000
3        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
5 L- O  m: O8 ^) H7 I& w' e, F9        20.00000            0.000000
% B; C& M6 ?. L" S10        0.000000            220.0000
11        40.00000            0.000000
12        0.000000            250.0000

alanloveai

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