写在前面的话 " _6 E" ]% O0 Z2 O8 w) p2 Q 1 U" v8 L. g- L8 S在学习和工作之余写点东西,这种想法由来已久了。只是总是觉得很忙,抽不出时间做一点点自己喜欢做的事情。突然有一天发现,你必须牺牲某些东西才能得到某些东西。所以,就索性写上几句话聊以**吧。 $ R( t. p! @* f- S* P最开始,只想翻译一下GAMS的用户手册。但是又觉得这样写很枯燥,不光读者不喜欢看,连自己也要打瞌睡了。如果想催眠的话,还不如去听评书。但是本书还是以解释GAMS的用法为主。中间夹杂着一些自己做modelling 的体验,还可能还有一些毫不相关的东西。可能有很多错误,请读者们原谅。想学习GAMS的同学,想必一定是数学,运筹学或计算机相关专业的了。让大家失望了,其实我的专业是chemical engineering, 本科时候做的是量子化学。阴差阳错上了贼船。但是觉得自己用了这么多年了,总该有点体会,能过帮帮大家。) j% i4 m$ I, \ B
% W7 W/ A: C; V0 ?( l- b
进入正题之前,还想啰嗦几句题外的话。 * L, v% _- H. ~7 g$ Y9 F3 R8 E9 `能看到这篇文字的同学,恐怕大多数是在灯下苦读的学子。实际上学习某种工具并不重要。真正要学习的是知识,是思维的方法。就像很多计算机系的学生一样,把数据结构和算法精通了,任何工具都是一样的。所以,你用GAMS或者不用GAMS,都没有关系,要把握事物的本质,最关键的东西抓住了,任何语言就是一种输入法。另外,学习一定要有兴趣,要认识到他的beauty,不要把什么都换算**民币,那样你会觉得很累。除了,coding 以外,还要博览群书,文史哲,拿来翻翻,就当休息了。4 |% C$ \) g7 p7 P# E$ B
最后,男子汉们,别忘了,谈一场轰轰烈烈的恋爱!! 7 [! L8 A) ~1 f% |. [' w1 P - m3 A+ f' P! N: V7 J % G8 E/ m* E7 b- p* J1 ^% A / \# _& n; T7 l; S % e% C1 }2 l, }: \* }. Q ( o* S: Z w6 E7 Y9 w/ Z4 ? _8 P7 @
第一章,基本概况" x% s/ F( o# I
8 K6 d6 p: f0 S) Y9 y0 tGAMS是一种非常简单易学的运筹学建模语言 (www.gams.com),是由WORLD BANK 开发的。由于近几十年计算机软硬件的的高速发展,世界银行决定投资开发GAMS, 来提供:! I1 ~6 Y0 `, q8 u$ [
1. 建立复杂模型的高级语言 4 C, T8 g; N9 D0 Z6 G/ K! [9 J2. 简单易学4 ~/ r4 u- E2 b- ]5 l& {
3. 强大的描述代数及逻辑关系; j! r! x y: }0 x, q0 a, R: |
4. 模型可以独立于算法系统 7 Q. K, r0 b1 @6 WGAMS 可以让使用者专心致力于模型的建立而省略了解法方面的考虑因为GAMS包含了解决通用的LP,NLP,MILP,MINLP等各类模型所需要的solver. 另外,GAMS提供了和外界进行交互的输入输出借口。使得用户很方便的导入导出data。 $ v4 [. p8 k& P5 B下面当然是tutorial了,就想学习任何一种编程语言一样,最开始都要学习输出 hello world. 我们这里当然要复杂一点了,举一个最简单的例子,运输问题.. ?3 X4 B; O2 ?3 v2 ^8 w, g
这是一个经典的LP问题(不用解释什么是LP吧,呵呵)) r% E- m- R7 _
用点英语,呵呵, 5 M( N% i5 x- K+ j H! J' EIn the familiar transportation problem, we are given the supplies at several plants and the demands at several markets for a single commodity, and we are given the unit costs of shipping the commodity from plants to markets. The economic question is: how much shipment should there be between each plant and each market so as to minimize total transport cost?3 ]! u# N3 }, f ]' Y+ D0 x" N
就是说,我们已知有一些工厂,生产出的产品要运往一些市场来满足他们的需求,每一个工厂的生产能力(每个工厂能提供产品的数量)以及每个市场的需求(每个市场需要产品的数量)是已知的。另外,我们还知道从特定的工厂到特定的市场,运输一件产品的运费。那么我们怎么样设计从每个工厂到每个市场运送货物的数量,使得总运费最小?; l3 @3 E8 `; \# c; }" m
Indices:$ A1 z% O5 a# ?# @ o
i = plants5 L( k. G4 s3 U( B
j = markets & ]' O0 |% @1 m0 {; uGiven Data: ; _5 D: O5 ~( q! B. f2 Vai = supply of commodity of plant i (in cases) 1 b0 V0 U& a8 w5 M7 E5 Q9 f0 o! Zbj = demand for commodity at market j (cases) " `! O. f9 ` K Wcij = cost per unit shipment between plant i and market j ($/case) 5 ~- C+ L, \0 U8 t) [7 q1 c1 C+ [Decision Variables:: \/ O/ U# X/ h2 Q9 r
Xij = amount of commodity to ship from plant i to market j 8 {) Y& g, i7 V0 Z/ G在这里唯一的连续变量是Xi,j,4 p+ Q8 j" ?) V. o# {) Y2 ~
模型就不用我说了吧。。。。; f; j, w$ T! r! @2 Q( _ R7 E) V
4 c" V3 B; a \3 }! z这是一个典型的线性规划问题,他的GAMS code is' T7 Y' H! C' r# t* l/ {! \
7 A2 p S |+ s. t! pSets : g5 Z7 [7 w% p6 Wi canning plants / seattle, san-diego /7 [7 v/ s& H* H& H. ^5 N5 p9 f
j markets / new-york, chicago, topeka / ;: d) h" [2 ~: }" S& f
Parameters- u1 C5 O! E6 E6 _/ O/ ?" x& J" H
a(i) capacity of plant i in cases" D, n; b9 [3 U1 R
/ seattle 350% m5 O% I; c5 n+ d. N- X
san-diego 600 /5 m9 F9 Z6 i5 p
b(j) demand at market j in cases , v- d0 v+ m" G+ e. G7 B/ new-york 3259 E' Y, Q$ o# M& Q, y* v
chicago 300* ?( p7 n2 j2 u; E9 L) m2 |
topeka 275 / ; 6 E, P2 [- j; ~; b0 f, q- oTable d(i,j) distance in thousands of miles z1 p5 d" l! P5 A( O9 w8 c- znew-york chicago topeka 2 {' G; |$ @5 _% jseattle 2.5 1.7 1.8 & {# @2 p4 t# @" z9 Bsan-diego 2.5 1.8 1.4 ; & B+ o6 v% s8 dScalar f freight in dollars per case per thousand miles /90/ ;9 P( E! Z8 H7 I B+ s
Parameter c(i,j) transport cost in thousands of dollars per case ; 6 X7 E# E, K g, }2 C5 T' h6 S5 q4 Cc(i,j) = f * d(i,j) / 1000 ; ' I$ Y, z) e, |6 e) QVariables, P5 v# `4 {, R, E% r6 ?$ M
x(i,j) shipment quantities in cases ( l R$ X$ ?3 U9 | nz total transportation costs in thousands of dollars ; L( I) R! C1 t8 H( Z/ lPositive Variable x ; ! F; J( I7 r. q. w, sEquations 4 S. A% {/ I6 c# F% Mcost define objective function 3 D$ k' n N4 ]* N& xsupply(i) observe supply limit at plant i2 c6 a/ D ?9 A' i5 q9 I
demand(j) satisfy demand at market j ;1 M0 l( M' l+ F3 k7 d$ p0 N$ }
cost .. z =e= sum((i,j), c(i,j)*x(i,j)) ; ! q3 y# B3 Q! e4 Gsupply(i) .. sum(j, x(i,j)) =l= a(i) ; 9 h5 V5 |) n; _+ Cdemand(j) .. sum(i, x(i,j)) =g= b(j) ; 7 }( |* G- p: N7 e) {* mModel transport /all/ ; " [5 n. P W# @2 s% U7 y/ |Solve transport using lp minimizing z ;$ R8 n2 x {2 E/ v: W
Display x.l, x.m ; 9 E( Z' n6 j6 n% h& n: M这是一个非常典型的LP(linear programming)问题,所有的变量和约束都是线性相关的。