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

relay1987

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

8 c9 N1 L1 l6 c8 F	甲 3 d# p3 q4 G* F  i. T0 n8 c  n	乙	丙 1 v: [9 ~1 k3 G  _	丁
A ' P% e( z$ Y- v* B/ N6 S	160 & d( l# H) x6 I	130	220 ! ]% {. b2 O; G8 E5 f; |	170
B	140 8 U$ V) L% n* f3 N, e$ K	130 8 \! y) V1 g- ^	190	150
C : b! }1 A' T8 n2 p, X% ]' F( H$ E# @) H	190	200 - M/ j" j8 _) \7 h/ g# F7 L	230 ! E0 ~7 K2 j/ d. }) D# J	/

) p% }* v+ t8 m: M& F' Z


9 T4 Y6 o# i$ e8 q$ P) u
         单位：元/千吨
针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

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

mrx

【1】求解：
7 C: K  Q% V( \; _$ [设水库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
) b. q; e* `2 U1 w  j& R4 Lst
& n3 a$ t( Q3 s9 Z- N' Fx34=0;
+ l% D5 o+ L5 {( l9 M3 J! @, Bx11+x12+x13+x14<=50;
9 F. n9 F( g2 @3 z) Ix21+x22+x23+x24<=60;
x31+x32+x33+x34<=50;
x11+x21+x31>=30;
- d2 v7 e$ e- I9 B  E. nx11+x21+x31<=80;
: X* K5 z0 B  t9 U8 `# Kx12+x22+x32>=70;
; s0 [7 t( ]4 r0 S+ n/ Xx12+x22+x32<=140;
x13+x23+x33>=10;
" N: S; j, S- U+ P. U  }x13+x23+x33<=30;
x14+x24+x34>=10;
& E0 `6 Y2 |9 c1 D" Zx14+x24+x34<=50;
end
+ Z' u4 h* X& [% t3 G, B6 ?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<=50;
x21+x22+x23+x24<=60;
  I# U# s+ ~) {) u. tx31+x32+x33+x34<=50;
  _/ H  Z5 A# c' Ix11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
$ F) C( {( q* l: c! C) _* l3 m- j( }x13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
x14+x24+x34<=50;
end
运行结果：
  Global optimal solution found.
  Objective value:                              47600.00
) `1 \4 v1 N5 U1 r  Total solver iterations:                             7
( X* {- `# A% a' Q
: L( C( m" u5 {) Q
Variable           Value        Reduced Cost
; Y3 I2 V3 I8 b4 F) `4 SX11        0.000000            30.00000
1 `" k. y3 T/ N% h& }: _X12        50.00000            0.000000
( \2 h' T* S9 H9 XX13        0.000000            50.00000
& }1 X; T, R1 u: b# |' r3 H5 PX14        0.000000            20.00000
X21        0.000000            10.00000
X22        50.00000            0.000000
& L/ F9 D, f: V& X0 C6 D$ EX23        0.000000            20.00000
, R" ^9 W- Y5 Y) Q# GX24        10.00000            0.000000
, {# _8 r0 {9 b* u0 M. p% n( IX31        40.00000            0.000000
X32        0.000000            10.00000
9 x/ ]6 [* P0 x4 W: W7 ?X33        10.00000            0.000000
7 K# \" I# M2 L, f' F: n6 aX34        0.000000            240.0000
0 r6 o) w3 P4 Z7 @0 i* S
Row    Slack or Surplus      Dual Price
1        47600.00            1.000000
0 [& m! _5 F" A) s6 e* y; T, X: ^2        0.000000            320.0000
3        0.000000            320.0000
4        0.000000            260.0000
5        10.00000            0.000000
6        40.00000            0.000000
7        30.00000            0.000000
8        40.00000            0.000000
% [7 i- n+ k1 I6 X9        0.000000           -40.00000
10        20.00000            0.000000
11        0.000000           -20.00000
0 a% K! P3 i7 C12        40.00000            0.000000

【2】求解：
6 s. L2 P  k; k: @模型：
0 m! s0 l1 e9 Z& wmax 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
/ M: Y$ W8 A. zst
x34=0;
  Y9 r$ [1 ?0 yx11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
+ k" a( a% W" R) R: mx11+x21+x31>=30;
# [3 `6 R* O# X5 Z& w% ?% Ox11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
  Y  J- ^' Z# N% E5 S2 i, `" e, B2 Gx13+x23+x33<=30;
8 @: B3 P$ r$ o. Q+ P& Z6 \x14+x24+x34>=10;
x14+x24+x34<=50;
. t' Y( I$ |/ iend
8 k. N: m$ J7 N+ |3 lLINGO程序：
max=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
  f2 G) W9 }8 }$ p2 T! M7 Wx11+x12+x13+x14<=100;
$ [  U& F! j/ O7 ^1 b1 ex21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
  u+ u( A9 S. \x11+x21+x31>=30;
5 n# Z& U" X) y& c. c& nx11+x21+x31<=80;
! a& s( C8 N9 N+ J1 E5 o) Ox12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
; e+ ~% p' P9 w$ X8 Qx14+x24+x34<=50;
end
' `  B9 e8 J) d0 y: h运行结果：
  Global optimal solution found.
/ _: ^2 ?' X+ e  Objective value:                              88700.00
  Total solver iterations:                             7
4 a) D3 U' O( l( B% f
/ ?, C6 M9 c% Y9 V6 R1 X6 b
- N0 k# f) o* h% F+ X7 m- s, B" uVariable           Value        Reduced Cost
4 n, ~5 \) J8 X1 FX11        0.000000            20.00000
X12        100.0000            0.000000
X13        0.000000            40.00000
6 f) V  y' K# `+ _) aX14        0.000000            20.00000
" B( u9 h+ x$ d1 V5 UX21        30.00000            0.000000
6 O+ v! ~$ @6 z  c& ?7 O9 FX22        40.00000            0.000000
X23        0.000000            10.00000
0 F3 n$ P& D6 R/ M. N2 M! }X24        50.00000            0.000000
: u' r8 H2 G; e2 ?6 I6 RX31        50.00000            0.000000
) y! n: I1 p4 A) R2 ]) oX32        0.000000            20.00000
X33        30.00000            0.000000
X34        0.000000            250.0000
0 [# G0 ]: z6 K8 |
Row    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
5        50.00000            0.000000
6        0.000000            260.0000
# |+ Q2 C+ R4 Q' ^7        70.00000            0.000000
' w! L6 v$ i% }& J! w5 A, x8        0.000000            270.0000
. ~5 j$ e8 {% b9        20.00000            0.000000
6 r6 \! H9 h  E10        0.000000            220.0000
11        40.00000            0.000000
12        0.000000            250.0000

alanloveai

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