GAMS中文教程

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

二GAMS简明教程由Richard E. Rosenthal编著

(译者:该翻译由华中科技大学能源与动力工程学院小海完成,由于译者水平有限,其中肯定会有些问题,欢迎大

家阅读并批评指正!)

2.1 简介

这本书的简介部分以一个详细的例子结束,这个例子是一个用GAMS进行表达、解决和分析的小而简单的优化问题。

由来自蒙特利尔市Naval研究生院的Richard E.Rownthal编写。通过这个例子我们可以对GAMS有一个快速而全面的

了解。很多地方引用到本书其他章节的知识,方便大家进一步查阅;而阅读这一部分即使不看本书其他章节也能看

懂并且有所收获。

例子是一个关于运输的线性规划问题,这个问题在最优化技术发展过程中曾经被当做实验对象。用这个例子展示

GAMS这类代数建模语言的强大功能是很好的选择,因为不管手头要处理的问题有多大,运输问题本身具有一个简

单的、可用的代数结构。你可以发现,如果要解决一个非常大的运输问题,我们将呈现的用于解决小规模运输问题

的GAMS输入文件中的内容并不需要改变多少。

在熟悉的运输问题中,我们考虑几个工厂的供应和几个市场的需求的商品,我们也给出从工厂运输商品到市场的单

位花费。这其中的经济学问题是:怎样安排运输使得我们的总运输成本最小?

这个问题的代数表达常常采用下列方式:

Indices:

i = plants

j = markets

Given Data:

i

a= supply of commodity of plant i (in cases)

ij

b= demand for commodity at market j

ij

c = cost per unit shipment between plant i an

d market j ($/case)

Decision Variables:

ij

x = amount of commodity to ship from plant i to market j (cases),

Where

ij

x ≥0, for all i, j

Constraints:

Observe supply limit at plant i:

j

ij j

a x≤∑for all i (cases)

Satisfy demand at market j:

ij j

i

x b ≥∑for all j (cases)

Objective Function: Minimize

ij ij

ijcx ∑∑($K)

注意这个例子显示了一些一般情况下我们认为是好习惯的建模方式,这些在GAMS里面被继承了。首先,模型程序

中的各组成部分都按类型分组定义了。其次,GAMS程序各组成部分先后次序已经定好了,因此没有标识符能在被

定义之前使用。第三,各组成部分都有特定的单位。第四,选择的单位要使优化过程中得到的数值具有相对较小的

绝对数量级。(例如标识$K表示dollar的千倍)

各组成部分类型的名称在不同的模型中可能不一样。例如,经济学家分别用外生变量和内生变量来表示已知数据和

决策变量。(译者:在经济模型中,内生变量是指该模型所要决定的变量。外生变量指由模型以外的因素所决定的

已知变量,它是模型据以建立的外部条件。内生变量可以在模型体系内得到说明,外生变量决定内生变量,而外生

变量本身不能在模型体系中得到说明。参数通常是由模型以外的因素决定的,因此也往往被看成外生变量。例:

P=a+bQ,表示价格与数量的关系,则a、b是参数,都是外生变量;P、Q是模型要决定的变量,所以是内生变量。

除此之外,譬如相关商品的价格,人们的收入等其他于模型有关的变量,都是外生变量)在GAMS中,被采用的相关术语是:sets表示指数下标,parameters表示已知数,variables 表示决策变量,equations

表示约束方程和目标方程。

运输问题的GAMS语言表述紧密的联系了上述几个部分。最主要的区别在于GAMS表述可以被电脑读取和运行。

作为运输问题的例子,假设有两个罐头厂和三个市场,已知数据如表2.1所示。运输距离的单位是千英里,运输成

本是$90.00每箱每千英里。这个例子的GAMS表述是:

Sets

i canning plants / seattle, san-diego /

j markets / new-york, chicago, topeka / ;

Parameters

a(i) capacity of plant i in cases

/ seattle 350

san-diego 600 /

b(j) demand at market j in cases

/ new-york 325

chicago 300

topeka 275 / ;

Table d(i,j) distance in thousands of miles

new-york chicago topeka

seattle 2.5 1.7 1.8

san-diego 2.5 1.8 1.4 ;

Scalar f freight in dollars per case per thousand miles /90/ ;

Parameter c(i,j) transport cost in thousands of dollars per case ;

c(i,j) = f * d(i,j) / 1000 ;

Variables

x(i,j) shipment quantities in cases

z total transportation costs in thousands of dollars ;

Positive Variable x ;

Equations

cost define objective function

supply(i) observe supply limit at plant i

demand(j) satisfy demand at market j ;

cost .. z =e= sum((i,j), c(i,j)*x(i,j)) ;

supply(i) .. sum(j, x(i,j)) =l= a(i) ;

demand(j) .. sum(i, x(i,j)) =g= b(j) ;

Model transport /all/ ;

Solve transport using lp minimizing z ;

Display x.l, x.m ;

如果你在GAMS中建立一个文件,将以上内容输入进去,运输模型就可以被建立和进行计算了。要使GAMS在不同

的计算机上运行需要改变一些细节,但是最简单的方法(不提供非必要服务的方法)是在输入文件的名字后面加上

GAMS这几个字母。在程序运行过程中,你将看到一些精练的描述GAMS运行过程的字符行,包括了写入输出结果

相关文档
最新文档