lingo从入门到高手 Lingo教程 下载本文

5 5 (B,M,1) 6 6 (B,M,2) 7 7 (B,N,1) 8 8 (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产生。

页 第5

另一方面,派生集是由其它的集来创建。这些集被称为该派生集的父集(原始集或其它的派生集)。一个派生集既可以是稀疏的,也可以是稠密的。稠密集包含了父集成员的所有组合(有时也称为父集的笛卡尔乘积)。稀疏集仅包含了父集的笛卡尔乘积的一个子集,可通过显式罗列和成员资格过滤器这两种方式来定义。显式罗列方法就是逐个罗列稀疏集的成员。成员资格过滤器方法通过使用稀疏集成员必须满足的逻辑条件从稠密集成员中过滤出稀疏集的成员。不同集类型的关系见下图。

§3 模型的数据部分和初始部分

集 派生集 原始集 稀疏集 稠密集 显式罗列 过滤器 LINGO集类型

在处理模型的数据时,需要为集指派一些成员并且在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 sets:

页 第6

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,然后继续求解模型。

除了参数之外,也可以实时输入集的属性值,但不允许实时输入集成员名。 3.1.4 指定属性为一个值

可以在数据声明的右边输入一个值来把所有

页 第7

的成员的该属性指定为一个值。看下面的例子。 例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;

好的初始点会减少模型的求解时间。

在这一节中,我们仅带大家接触了一些基本的数据输入和初始化概念,不过现在你应该可以轻松的为自己的模型加入原始数据和初始部分啦。 §4 LINGO函数

页 第8