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

relay1987

4 t5 s6 F* z7 E6 w1 M1 V: r6 P) [
有A、B、C三个水库往甲、乙、丙、丁四个小区供应自来水，其中C水库由于条件限制不能给丁小区供水，这三个水库每天的供水量分别为50，60，50（单位：千吨），这四个小区每天的基本用水量为30，70，10，10（单位：千吨），在保证基本用水量的前提下，每个小区还需要一定量的额外用水，这四个小区每天额外需求量分别为50，70，20，40（单位：千吨），小区只需按照900元/千吨的价格缴纳费用，水库往这四个小区的引水管理费如下表，且除引水管理费之外，每千吨水还有450元的其他费用。
$ ^4 G+ u7 j1 g$ o. ~* H

	甲	乙 * X9 H4 m3 m* i* X' o& y" `	丙	丁 6 m0 H9 A/ ~6 a- t/ T& A3 D7 f
A	160	130	220	170
B	140 ; R  k" U. W6 D  C	130	190	150 - Z$ P6 g. L' \+ S' e
C + L* B  k; E7 z9 V	190	200 : b8 ], x9 j: o/ r& J	230 ( M. j- F9 x6 U: |$ B5 z7 k- K	/

         单位：元/千吨
针对下面问题建立合适的数学模型，并写出程序（用Matlab、Mathematic、lingo均可）
- f0 R2 S/ L* o8 P: x1. 应如何分配水库供水量，公司才能获利最多？
2. 若水库供水量都提高一倍，公司利润可增加到多少？

relay1987

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

mrx

【1】求解：
; ~" V+ O: x' ~6 q0 C. v. z" G设水库A、B、C分别向甲、乙、丙、丁四个小区供水量为x11、x12、x13、x14、x21、x22、x23、x24、x31、x32、x33、x34，则有：
! x* ]3 `9 }0 o9 ~+ I+ L* e+ i# Vmax 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
6 j4 `' M( x7 V( wst
" Q+ F- |/ s( [4 T' Dx34=0;
3 i( J' g! T% l& qx11+x12+x13+x14<=50;
. F# c/ A9 K- w: a# [, Nx21+x22+x23+x24<=60;
x31+x32+x33+x34<=50;
& R+ l* z" z% _4 wx11+x21+x31>=30;
, [( s0 V/ z5 ~( p  S5 [( I; Ex11+x21+x31<=80;
x12+x22+x32>=70;
x12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
  R% q" c+ z0 G' H- \1 Yx14+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;
3 @3 t) g  ]) _/ M+ @; u' q. m# Vx11+x12+x13+x14<=50;
0 L5 |4 e: M( t2 I4 I+ wx21+x22+x23+x24<=60;
. T7 G) i1 ~. T& ]- G6 G( D! \4 `, ]$ mx31+x32+x33+x34<=50;
x11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
! r& x: p4 l& `. P! T2 p% Fx12+x22+x32<=140;
x13+x23+x33>=10;
x13+x23+x33<=30;
x14+x24+x34>=10;
: W& X7 ]) F( u( Rx14+x24+x34<=50;
, Q/ `: P0 O! S1 W- Q/ ~' hend
3 J( D, H) ]1 ^/ u( n运行结果：
  Global optimal solution found.
  Objective value:                              47600.00
# J$ `; T) G' \4 P: T8 o6 k  Total solver iterations:                             7


2 V0 R0 S& ?, W; Y% @6 Y  qVariable           Value        Reduced Cost
X11        0.000000            30.00000
X12        50.00000            0.000000
9 J& L  ]7 [0 a, z. H9 j% xX13        0.000000            50.00000
X14        0.000000            20.00000
X21        0.000000            10.00000
: M1 G5 V+ [4 L" M3 [X22        50.00000            0.000000
X23        0.000000            20.00000
4 f9 q" V/ Q/ TX24        10.00000            0.000000
! H$ U1 M  i, K* o" t( EX31        40.00000            0.000000
1 M  |- M/ x( o/ EX32        0.000000            10.00000
9 Q0 I" }3 U5 }0 e/ O# QX33        10.00000            0.000000
0 W3 R6 T( g; H) A7 i1 tX34        0.000000            240.0000
4 Q5 k$ z0 t+ p
Row    Slack or Surplus      Dual Price
1        47600.00            1.000000
# |0 o) \5 _' ^" D. r3 k* n- w+ G2        0.000000            320.0000
3        0.000000            320.0000
' u6 {4 A1 O" s4        0.000000            260.0000
5        10.00000            0.000000
6        40.00000            0.000000
7        30.00000            0.000000
5 A* x# o: A0 V6 r, U! h+ P8        40.00000            0.000000
# \% y2 Y# z  ^, [1 p9        0.000000           -40.00000
1 R( B' j1 ]# J2 ]4 B. Z: m: v" Y10        20.00000            0.000000
11        0.000000           -20.00000
. q0 T' t, U" Z; s5 ?12        40.00000            0.000000

【2】求解：
模型：
3 d4 V( H' @4 F6 }2 u7 Z; _max 290x11+320x12+230x13+280x14+310x21+320x22+260x23+300x24+260x31+250x32+220x33
+ y! v1 B( J1 A% F8 k% l. Q( kst
x34=0;
x11+x12+x13+x14<=100;
x21+x22+x23+x24<=120;
; n! O7 G  ^/ q$ u: ?$ x3 Lx31+x32+x33+x34<=100;
x11+x21+x31>=30;
x11+x21+x31<=80;
x12+x22+x32>=70;
0 z% C( w, r' a$ Hx12+x22+x32<=140;
" v" W2 _: F  }$ N* }  @x13+x23+x33>=10;
, ~8 }4 R! N6 B" k9 H3 k* Kx13+x23+x33<=30;
- N2 |+ T5 [3 ex14+x24+x34>=10;
x14+x24+x34<=50;
! A, t& z; Q2 d3 {end
' V- C1 q2 {! h# W& p( k2 TLINGO程序：
" @7 h! L5 u* a: tmax=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;
x21+x22+x23+x24<=120;
x31+x32+x33+x34<=100;
x11+x21+x31>=30;
0 Q0 z! m9 E; |. q! z9 H1 [x11+x21+x31<=80;
+ _5 z5 k$ m9 O4 c% _% mx12+x22+x32>=70;
x12+x22+x32<=140;
! }9 i3 Q* _# q8 t8 B: sx13+x23+x33>=10;
! x* G) J& O% m- b- `! ^& E1 `x13+x23+x33<=30;
3 i/ P9 B7 I0 G4 l# }x14+x24+x34>=10;
6 i. N" b* U" {1 e& t8 n0 gx14+x24+x34<=50;
end
运行结果：
# x; r( X& D: Z  Global optimal solution found.
6 T0 R: m: S! f  Objective value:                              88700.00
  Total solver iterations:                             7
- x1 w" f9 j7 p# c' S

( @7 j4 h- U( ^0 [' D0 @Variable           Value        Reduced Cost
: u' D2 ?& N, n+ W' G- eX11        0.000000            20.00000
) j% S/ X- D6 B* o9 n) gX12        100.0000            0.000000
X13        0.000000            40.00000
X14        0.000000            20.00000
/ K" Q' k- p' `2 n' _3 yX21        30.00000            0.000000
8 z8 j3 Y4 |, q2 l* {# n3 D5 y6 zX22        40.00000            0.000000
* O5 c  I7 a4 j: W0 H6 H& q  rX23        0.000000            10.00000
6 [4 p8 f. i* R$ |X24        50.00000            0.000000
, [: }; J# N) }- W  h6 ~9 w% iX31        50.00000            0.000000
X32        0.000000            20.00000
X33        30.00000            0.000000
X34        0.000000            250.0000
7 n2 a1 ]) Q+ p' i9 u' |
Row    Slack or Surplus      Dual Price
' \' l# E  @/ t5 \. U1        88700.00            1.000000
2        0.000000            50.00000
3        0.000000            50.00000
3 ~- f* ]. @' @4        20.00000            0.000000
7 a6 ~3 W) w5 n; z5        50.00000            0.000000
6        0.000000            260.0000
( {* S# Z2 h  d5 |8 o& C' P& D; S7        70.00000            0.000000
8        0.000000            270.0000
( K: Z# c& S6 g" U2 Z8 ?: N9        20.00000            0.000000
6 B0 x- S4 V9 J) f10        0.000000            220.0000
11        40.00000            0.000000
12        0.000000            250.0000

alanloveai

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