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

relay1987

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

' b5 ?( k! \6 [$ }7 v) _	甲 - D& @% B2 f: E5 y% M	乙	丙 : h# P' l1 C* R' B2 u	丁 + P9 S8 ~( K9 ]& ?  l) V  Y
A	160	130 8 L! k, I! q/ n# S/ e	220 , m3 e& f  L3 t, K' J* X( r	170
B 5 `2 _% G2 y! ~7 P* Q' O, o7 H- s	140 % Y/ |( U0 \2 E2 w  h7 M, d	130	190	150
C   l; H7 t6 M/ U1 D6 `	190	200	230	/ 8 w: l' ]# k8 E

9 P! d: n7 u! ]9 K% {6 D


0 W# ^) z3 s3 E, e% U6 F2 ]         单位：元/千吨
' g; S+ R$ A  m1 @8 w" T针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
. c  T# j8 l/ w/ Z1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

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

mrx

【1】求解：
+ ~& c' s3 X! g4 @' Q) c/ v; k设水库A、B、C分别向甲、乙、丙、丁四个小区供水量为x11、x12、x13、x14、x21、x22、x23、x24、x31、x32、x33、x34，则有：
; j$ z, t& h/ Bmax 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
, |" R# @8 x' e4 V9 nst
% t; ^) X" e0 ^. |3 [/ nx34=0;
x11+x12+x13+x14<=50;
" l' B% H' k& p& sx21+x22+x23+x24<=60;
x31+x32+x33+x34<=50;
2 a9 K8 O0 F$ Bx11+x21+x31>=30;
" U/ a- c: G  Q& W, F0 zx11+x21+x31<=80;
- A) ~8 S+ b& p3 T4 T# Xx12+x22+x32>=70;
& n& r% C0 v" _7 ]x12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
1 w1 J+ ~7 I3 a( kx14+x24+x34<=50;
end
" ^0 v* r+ c( cLINGO程序：
! Z1 ^( r; e9 Y( H  fmax=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;
: f  g% a" w4 v# M* G) L6 ax31+x32+x33+x34<=50;
7 g" q' F7 c$ W4 U: _( ]x11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
" }' ~2 k: n) |) ex13+x23+x33>=10;
; c# g2 n  J' Z4 Kx13+x23+x33<=30;
x14+x24+x34>=10;
x14+x24+x34<=50;
end
运行结果：
8 @: k+ q* v" {" I) w5 V  Global optimal solution found.
9 P, ?) M  I- a5 {7 Y0 ?  Objective value:                              47600.00
  Total solver iterations:                             7
+ S/ r7 }  Q% h6 o9 [% }+ l  L
" H9 X) T9 B0 V' \3 w: l* N. \
Variable           Value        Reduced Cost
; O4 @7 }  I& Z/ M: I. zX11        0.000000            30.00000
$ [- G1 `; f" N  gX12        50.00000            0.000000
X13        0.000000            50.00000
X14        0.000000            20.00000
X21        0.000000            10.00000
X22        50.00000            0.000000
& I- D/ e7 {) S$ k0 Y$ k' wX23        0.000000            20.00000
8 E* O" M. H0 |, dX24        10.00000            0.000000
# ]2 {8 p- i: l3 c; q( W# T% pX31        40.00000            0.000000
' g1 G7 \  Z! k5 YX32        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
- ]/ S) ]: W* i  J3 x4 |2        0.000000            320.0000
% r7 ^' O0 i' P- X* @3        0.000000            320.0000
0 X, g, s$ |+ W0 y4        0.000000            260.0000
5        10.00000            0.000000
; ^& o. h. g) R8 C( }! a6 v% ]. ?6        40.00000            0.000000
7        30.00000            0.000000
. Z: _) f' c+ D: P8        40.00000            0.000000
* U6 v7 Q! ^) \9 Z6 `6 Y! X( [9        0.000000           -40.00000
5 i( v0 Z9 l* }' {) j10        20.00000            0.000000
11        0.000000           -20.00000
12        40.00000            0.000000
9 e6 ]' Q) ~; F8 n$ U* L% G6 e
【2】求解：
1 w% x/ v' J' U, J模型：
max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
st
x34=0;
2 N; z4 f* h6 @! G4 o: u6 f# Px11+x12+x13+x14<=100;
; e" Z& _; T  v3 e* j: P/ C& Mx21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
x11+x21+x31>=30;
9 S5 k. L! [2 j4 }" Ex11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
0 b! g0 R  f' `' l" L3 k, o* Nx13+x23+x33<=30;
7 F9 }. @# G( M$ @6 U1 p8 k. Bx14+x24+x34>=10;
; Z* _4 [  |- o) `; \  D- ix14+x24+x34<=50;
end
7 F: _+ k# T. w& T- v. SLINGO程序：
max=290*x11+320*x12+230*x13+280*x14+310*x21+320*x22+260*x23+300*x24+260*x31+250*x32+220*x33;
- o: \2 m, ^# ux11+x12+x13+x14<=100;
8 i# C5 W" m" r+ sx21+x22+x23+x24<=120;
: b& w& B0 c' V1 [$ _- Dx31+x32+x33+x34<=100;
x11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
' `0 G7 q. m( z! zx12+x22+x32<=140;
x13+x23+x33>=10;
' g) E/ t- y7 B6 p6 U0 ]. g- ix13+x23+x33<=30;
x14+x24+x34>=10;
# |5 }: r8 x7 G& J0 y0 Q# bx14+x24+x34<=50;
end
  e' }/ u) B5 S5 z, O5 S; N7 H" q运行结果：
  z4 Q* X# w+ j  Global optimal solution found.
  Objective value:                              88700.00
  Total solver iterations:                             7


+ F8 D2 N5 `6 ?0 Z2 k0 r- |Variable           Value        Reduced Cost
5 I& a4 E3 L9 c# h. tX11        0.000000            20.00000
X12        100.0000            0.000000
X13        0.000000            40.00000
( S' e0 r8 H1 W' ?X14        0.000000            20.00000
: p: O1 R, f+ I$ b+ L6 |# f: YX21        30.00000            0.000000
X22        40.00000            0.000000
X23        0.000000            10.00000
X24        50.00000            0.000000
. M0 d8 ]- L/ N  {# E  wX31        50.00000            0.000000
5 }/ W* {0 s7 x1 d3 DX32        0.000000            20.00000
X33        30.00000            0.000000
X34        0.000000            250.0000
& B/ @0 S( @* F. [! n  u! |
Row    Slack or Surplus      Dual Price
! s1 J2 ~) E, `4 }* |1        88700.00            1.000000
2        0.000000            50.00000
3        0.000000            50.00000
$ a. s& n$ ?4 ~3 P& b6 N$ i9 ^( o4        20.00000            0.000000
5        50.00000            0.000000
6        0.000000            260.0000
& k# d3 e2 j& y+ Z) g" A: L" c8 \; a/ E7        70.00000            0.000000
8        0.000000            270.0000
) d. k$ t( Z4 @8 g3 |0 |6 a9 w9        20.00000            0.000000
9 e$ |' k7 Z4 m! l* l1 G7 R6 R10        0.000000            220.0000
) g2 A" e: t7 x# B4 @11        40.00000            0.000000
* m! N7 D6 L9 X9 b# A12        0.000000            250.0000

alanloveai

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