LINGO 教程
LINGO 是用来求解线性和非线性优化问题的简易工具。LINGO 内置了一种建立最优化模型
的语言,可以简便地表达大规模问题,利用 LINGO 高效的求解器可快速求解并分析结果。
§1
LINGO 快速入门
当你在 windows 下开始运行 LINGO 系统时,会得到类似下面的一个窗口:
外层是主框架窗口,包含了所有菜单命令和工具条,其它所有的窗口将被包含在主窗口之下。
在主窗口内的标题为 LINGO Model – LINGO1 的窗口是 LINGO 的默认模型窗口,建立的模型
都都要在该窗口内编码实现。下面举两个例子。
例 1.1 如何在 LINGO 中求解如下的 LP 问题:
3
2
x
1
min
..
ts
x
2
x
1
x
1
2
x
1
,
xx
1
2
x
2
x
2
350
100
600
0
在模型窗口中输入如下代码:
min=2*x1+3*x2;
x1+x2>=350;
x1>=100;
2*x1+x2<=600;
然后点击工具条上的按钮
即可。
例 1.2 使用 LINGO 软件计算 6 个发点 8 个收点的最小费用运输问题。产销单位运价如
下表。
单
位
销地
运
价
产地
B1
B2
B3
B4
B5
B6
B7
B8 产量
A1
A2
6
4
2
9
6
5
7
3
4
8
2
5
5
8
9
2
60
55
共 53 页 1
A3
A4
A5
A6
5
7
2
5
2
6
3
5
1
7
9
2
9
3
5
2
7
9
7
8
4
2
2
1
3
7
6
4
LINGO 教程
3
1
5
3
51
43
41
52
销量
35
37
22
32
41
32
43
38
使用 LINGO 软件,编制程序如下:
model:
!6 发点 8 收点运输问题;
sets:
warehouses/wh1..wh6/: capacity;
vendors/v1..v8/: demand;
links(warehouses,vendors): cost, volume;
endsets
!目标函数;
min=@sum(links: cost*volume);
!需求约束;
@for(vendors(J):
@sum(warehouses(I): volume(I,J))=demand(J));
!产量约束;
@for(warehouses(I):
@sum(vendors(J): volume(I,J))<=capacity(I));
!这里是数据;
data:
capacity=60 55 51 43 41 52;
demand=35 37 22 32 41 32 43 38;
cost=6 2 6 7 4 2 9 5
4 9 5 3 8 5 8 2
5 2 1 9 7 4 3 3
7 6 7 3 9 2 7 1
2 3 9 5 7 2 6 5
5 5 2 2 8 1 4 3;
enddata
end
然后点击工具条上的按钮
为了能够使用 LINGO 的强大功能,接着第二节的学习吧。
即可。
§2
LINGO 中的集
对实际问题建模的时候,总会遇到一群或多群相联系的对象,比如工厂、消费者群体、交通
工具和雇工等等。LINGO 允许把这些相联系的对象聚合成集(sets)。一旦把对象聚合成集,
就可以利用集来最大限度的发挥 LINGO 建模语言的优势。
现在我们将深入介绍如何创建集,并用数据初始化集的属性。学完本节后,你对基于建
模技术的集如何引入模型会有一个基本的理解。
2.1 为什么使用集
集是 LINGO 建模语言的基础,是程序设计最强有力的基本构件。借助于集,能够用一个
共 53 页 2
LINGO 教程
单一的、长的、简明的复合公式表示一系列相似的约束,从而可以快速方便地表达规模较大
的模型。
2.2 什么是集
集是一群相联系的对象,这些对象也称为集的成员。一个集可能是一系列产品、卡车或
雇员。每个集成员可能有一个或多个与之有关联的特征,我们把这些特征称为属性。属性值
可以预先给定,也可以是未知的,有待于 LINGO 求解。例如,产品集中的每个产品可以有一
个价格属性;卡车集中的每辆卡车可以有一个牵引力属性;雇员集中的每位雇员可以有一个
薪水属性,也可以有一个生日属性等等。
LINGO 有两种类型的集:原始集(primitive
一个原始集是由一些最基本的对象组成的。
一个派生集是用一个或多个其它集来定义的,也就是说,它的成员来自于其它已存在的
set)和派生集(derived set)。
集。
2.3 模型的集部分
集部分是 LINGO 模型的一个可选部分。在 LINGO 模型中使用集之前,必须在集部分事先
定义。集部分以关键字“sets:”开始,以“endsets”结束。一个模型可以没有集部分,或
有一个简单的集部分,或有多个集部分。一个集部分可以放置于模型的任何地方,但是一个
集及其属性在模型约束中被引用之前必须定义了它们。
2.3.1 定义原始集
为了定义一个原始集,必须详细声明:
·集的名字
·可选,集的成员
·可选,集成员的属性
定义一个原始集,用下面的语法:
setname[/member_list/][:attribute_list];
注意:用“[]”表示该部分内容可选。下同,不再赘述。
Setname 是你选择的来标记集的名字,最好具有较强的可读性。集名字必须严格符合标
准命名规则:以拉丁字母或下划线(_)为首字符,其后由拉丁字母(A—Z)、下划线、阿拉
伯数字(0,1,…,9)组成的总长度不超过 32 个字符的字符串,且不区分大小写。
注意:该命名规则同样适用于集成员名和属性名等的命名。
Member_list 是集成员列表。如果集成员放在集定义中,那么对它们可采取显式罗列和
隐式罗列两种方式。如果集成员不放在集定义中,那么可以在随后的数据部分定义它们。
① 当显式罗列成员时,必须为每个成员输入一个不同的名字,中间用空格或逗号搁开,
允许混合使用。
例 2.1 可以定义一个名为 students 的原始集,它具有成员 John、Jill、Rose 和 Mike,
属性有 sex 和 age:
sets:
students/John
Jill, Rose
Mike/: sex, age;
endsets
② 当隐式罗列成员时,不必罗列出每个集成员。可采用如下语法:
setname/member1..memberN/[: attribute_list];
这里的 member1 是集的第一个成员名,memberN 是集的最末一个成员名。LINGO 将自动产生
中间的所有成员名。LINGO 也接受一些特定的首成员名和末成员名,用于创建一些特殊的集。
列表如下:
隐式成员列表格式
示例
所产生集成员
1..n
1..5
1,2,3,4,5
StringM..StringN
Car2..car14
Car2,Car3,Car4,…,Car14
DayM..DayN
Mon..Fri
Mon,Tue,Wed,Thu,Fri
共 53 页 3
LINGO 教程
MonthM..MonthN
Oct..Jan
Oct,Nov,Dec,Jan
MonthYearM..MonthYearN
Oct2001..Jan2002
Oct2001,Nov2001,Dec2001,Jan2002
③ 集成员不放在集定义中,而在随后的数据部分来定义。
例 2.2
!集部分;
sets:
students:sex,age;
endsets
!数据部分;
data:
students,sex,age= John 1 16
Jill 0 14
Rose 0 17
Mike 1 13;
enddata
注意:开头用感叹号(!),末尾用分号(;)表示注释,可跨多行。
在集部分只定义了一个集 students,并未指定成员。在数据部分罗列了集成员 John、
Jill、Rose 和 Mike,并对属性 sex 和 age 分别给出了值。
集成员无论用何种字符标记,它的索引都是从 1 开始连续计数。在 attribute_ list 可
以指定一个或多个集成员的属性,属性之间必须用逗号隔开。
可以把集、集成员和集属性同 C 语言中的结构体作个类比。如下图:
集 ←→ 结构体
集成员 ←→ 结构体的域
集属性 ←→ 结构体实例
LINGO 内置的建模语言是一种描述性语言,用它可以描述现实世界中的一些问题,然后
再借助于 LINGO 求解器求解。因此,集属性的值一旦在模型中被确定,就不可能再更改。在
LINGO 中,只有在初始部分中给出的集属性值在以后的求解中可更改。这与前面并不矛盾,
初始部分是 LINGO 求解器的需要,并不是描述问题所必须的。
2.3.2 定义派生集
为了定义一个派生集,必须详细声明:
·集的名字
·父集的名字
·可选,集成员
·可选,集成员的属性
可用下面的语法定义一个派生集:
setname(parent_set_list)[/member_list/][:attribute_list];
setname 是集的名字。parent_set_list 是已定义的集的列表,多个时必须用逗号隔开。
如果没有指定成员列表,那么 LINGO 会自动创建父集成员的所有组合作为派生集的成员。派
生集的父集既可以是原始集,也可以是其它的派生集。
例 2.3
sets:
product/A B/;
machine/M N/;
week/1..2/;
allowed(product,machine,week):x;
endsets
LINGO 生成了三个父集的所有组合共八组作为 allowed 集的成员。列表如下:
编号
1
成员
(A,M,1)
共 53 页 4
2
3
4
5
6
7
8
2
3
4
5
6
7
8
LINGO 教程
(A,M,2)
(A,N,1)
(A,N,2)
(B,M,1)
(B,M,2)
(B,N,1)
(B,N,2)
成员列表被忽略时,派生集成员由父集成员所有的组合构成,这样的派生集成为稠密集。
如果限制派生集的成员,使它成为父集成员所有组合构成的集合的一个子集,这样的派生集
成为稀疏集。同原始集一样,派生集成员的声明也可以放在数据部分。一个派生集的成员列
表有两种方式生成:①显式罗列;②设置成员资格过滤器。当采用方式①时,必须显式罗列
出所有要包含在派生集中的成员,并且罗列的每个成员必须属于稠密集。使用前面的例子,
显式罗列派生集的成员:
allowed(product,machine,week)/A M 1,A N 2,B N 1/;
如果需要生成一个大的、稀疏的集,那么显式罗列就很讨厌。幸运地是许多稀疏集的成员都
满足一些条件以和非成员相区分。我们可以把这些逻辑条件看作过滤器,在 LINGO 生成派生
集的成员时把使逻辑条件为假的成员从稠密集中过滤掉。
例 2.4
sets:
!学生集:性别属性 sex,1 表示男性,0 表示女性;年龄属性 age.
students/John,Jill,Rose,Mike/:sex,age;
!男学生和女学生的联系集:友好程度属性 friend,[0,1]之间的数。 ;
linkmf(students,students)|sex(&1) #eq# 1 #and# sex(&2) #eq# 0: friend;
!男学生和女学生的友好程度大于 0.5 的集;
linkmf2(linkmf) | friend(&1,&2) #ge# 0.5 : x;
;
endsets
data:
sex,age = 1 16
0 14
0 17
0 13;
friend = 0.3 0.5 0.6;
enddata
用竖线(|)来标记一个成员资格过滤器的开始。#eq#是逻辑运算符,用来判断是否“相
等”,可参考§4. &1 可看作派生集的第 1 个原始父集的索引,它取遍该原始父集的所有成
员;&2 可看作派生集的第 2 个原始父集的索引,它取遍该原始父集的所有成员;&3,&4,……,
以此类推。注意如果派生集 B 的父集是另外的派生集 A,那么上面所说的原始父集是集 A 向
前回溯到最终的原始集,其顺序保持不变,并且派生集 A 的过滤器对派生集 B 仍然有效。因
此,派生集的索引个数是最终原始父集的个数,索引的取值是从原始父集到当前派生集所作
限制的总和。
总的来说,LINGO 可识别的集只有两种类型:原始集和派生集。
在一个模型中,原始集是基本的对象,不能再被拆分成更小的组分。原始集可以由显式
罗列和隐式罗列两种方式来定义。当用显式罗列方式时,需在集成员列表中逐个输入每个成
员。当用隐式罗列方式时,只需在集成员列表中输入首成员和末成员,而中间的成员由 LINGO
产生。
另一方面,派生集是由其它的集来创建。这些集被称为该派生集的父集(原始集或其它
的派生集)。一个派生集既可以是稀疏的,也可以是稠密的。稠密集包含了父集成员的所有
组合(有时也称为父集的笛卡尔乘积)。稀疏集仅包含了父集的笛卡尔乘积的一个子集,可
通过显式罗列和成员资格过滤器这两种方式来定义。显式罗列方法就是逐个罗列稀疏集的成
员。成员资格过滤器方法通过使用稀疏集成员必须满足的逻辑条件从稠密集成员中过滤出稀
疏集的成员。不同集类型的关系见下图。
共 53 页 5
LINGO 教程
集
派生集
原始集
稀疏集
稠密集
显式罗列
过滤器
LINGO 集类型
§3 模型的数据部分和初始部分
在处理模型的数据时,需要为集指派一些成员并且在 LINGO 求解模型之前为集的某些属
性指定值。为此,LINGO 为用户提供了两个可选部分:输入集成员和数据的数据部分(Data
Section)和为决策变量设置初始值的初始部分(Init Section)。
3.1 模型的数据部分
3.1.1 数据部分入门
数据部分提供了模型相对静止部分和数据分离的可能性。显然,这对模型的维护和维数
的缩放非常便利。
数据部分以关键字“data:”开始,以关键字“enddata”结束。在这里,可以指定集成
员、集的属性。其语法如下:
object_list = value_list;
对象列(object_list)包含要指定值的属性名、要设置集成员的集名,用逗号或空格
隔开。一个对象列中至多有一个集名,而属性名可以有任意多。如果对象列中有多个属性名,
那么它们的类型必须一致。如果对象列中有一个集名,那么对象列中所有的属性的类型就是
这个集。
数值列(value_list)包含要分配给对象列中的对象的值,用逗号或空格隔开。注意属
性值的个数必须等于集成员的个数。看下面的例子。
例 3.1
sets:
set1/A,B,C/: X,Y;
endsets
data:
X=1,2,3;
Y=4,5,6;
enddata
在集 set1 中定义了两个属性 X 和 Y。X 的三个值是 1、2 和 3,Y 的三个值是 4、5 和 6。
也可采用如下例子中的复合数据声明(data statement)实现同样的功能。
例 3.2
共 53 页 6
LINGO 教程
sets:
set1/A,B,C/: X,Y;
endsets
data:
X,Y=1 4
2 5
3 6;
enddata
看到这个例子,可能会认为 X 被指定了 1、4 和 2 三个值,因为它们是数值列中前三个,
而正确的答案是 1、2 和 3。假设对象列有 n 个对象,LINGO 在为对象指定值时,首先在 n
个对象的第 1 个索引处依次分配数值列中的前 n 个对象,然后在 n 个对象的第 2 个索引处依
次分配数值列中紧接着的 n 个对象,……,以此类推。
模型的所有数据——属性值和集成员——被单独放在数据部分,这可能是最规范的数据
输入方式。
3.1.2 参数
在数据部分也可以指定一些标量变量(scalar variables)。当一个标量变量在数据部
分确定时,称之为参数。看一例,假设模型中用利率 8.5%作为一个参数,就可以象下面一
样输入一个利率作为参数。
例 3.3
data:
interest_rate = .085;
enddata
也可以同时指定多个参数。
例 3.4
data:
interest_rate,inflation_rate = .085 .03;
enddata
3.1.3 实时数据处理
在某些情况,对于模型中的某些数据并不是定值。譬如模型中有一个通货膨胀率的参数,
我们想在 2%至 6%范围内,对不同的值求解模型,来观察模型的结果对通货膨胀的依赖有多
么敏感。我们把这种情况称为实时数据处理(what if analysis)。LINGO 有一个特征可方
便地做到这件事。
在本该放数的地方输入一个问号(?)。
例 3.5
data:
interest_rate,inflation_rate = .085
?;
enddata
每一次求解模型时,LINGO 都会提示为参数 inflation_rate 输入一个值。在 WINDOWS 操作
系统下,将会接收到一个类似下面的对话框:
直接输入一个值再点击 OK 按钮,LINGO 就会把输入的值指定给 inflation_rate,然后继续
求解模型。
共 53 页 7
除了参数之外,也可以实时输入集的属性值,但不允许实时输入集成员名。
3.1.4 指定属性为一个值
可以在数据声明的右边输入一个值来把所有的成员的该属性指定为一个值。看下面的例
LINGO 教程
子。
例 3.6
sets:
days /MO,TU,WE,TH,FR,SA,SU/:needs;
endsets
data:
needs = 20;
enddata
LINGO 将用 20 指定 days 集的所有成员的 needs 属性。对于多个属性的情形,见下例。
例 3.7
sets:
days /MO,TU,WE,TH,FR,SA,SU/:needs,cost;
endsets
data:
needs cost = 20 100;
enddata
3.1.5 数据部分的未知数值
有时只想为一个集的部分成员的某个属性指定值,而让其余成员的该属性保持未知,以
便让 LINGO 去求出它们的最优值。在数据声明中输入两个相连的逗号表示该位置对应的集成
员的属性值未知。两个逗号间可以有空格。
例 3.8
sets:
years/1..5/: capacity;
endsets
data:
capacity = ,34,20,,;
enddata
属性 capacity 的第 2 个和第 3 个值分别为 34 和 20,其余的未知。
3.2 模型的初始部分
初 始 部 分 是 LINGO 提 供 的 另 一 个 可 选 部 分 。 在 初 始 部 分 中 , 可 以 输 入 初 始 声 明
(initialization statement),和数据部分中的数据声明相同。对实际问题的建模时,初
始部分并不起到描述模型的作用,在初始部分输入的值仅被 LINGO 求解器当作初始点来用,
并且仅仅对非线性模型有用。和数据部分指定变量的值不同,LINGO 求解器可以自由改变初
始部分初始化的变量的值。
一个初始部分以“init:”开始,以“endinit”结束。初始部分的初始声明规则和数据
部分的数据声明规则相同。也就是说,我们可以在声明的左边同时初始化多个集属性,可以
把集属性初始化为一个值,可以用问号实现实时数据处理,还可以用逗号指定未知数值。
例 3.9
init:
X, Y = 0, .1;
endinit
Y=@log(X);
X^2+Y^2<=1;
好的初始点会减少模型的求解时间。
在这一节中,我们仅带大家接触了一些基本的数据输入和初始化概念,不过现在你应该
可以轻松的为自己的模型加入原始数据和初始部分啦。
共 53 页 8