LINGO的使用方法说明大全Word文档下载推荐.docx

1、3.2 LINGO中集合的概念在对实际问题建模的时候，总会遇到一群或多群相联系的对象，比如工厂、消费者群体、交通工具和雇工等等LINGO允许把这些相联系的对象聚合成集（sets）一旦把对象聚合成集，就可以利用集来最大限度地发挥LINGO建模语言的优势现在将深入介绍如何创建集，并用数据初始化集的属性3.2.1集的构成集是LINGO建模语言的基础，是程序设计最强有力的基本构件借助于集能够用一个单一的、简明的复合公式表示一系列相似的约束，从而可以快速方便地表达规模较大的模型集是一群相联系的对象，这些对象也称为集的元素一个集可能是一系列产品、卡车或雇员每个集的元素可能有一个或多个与之有关联的特征，把这

2、些特征称为属性属性值可以预先给定，也可以是未知的，有待于LINGO求解的LINGO有两种类型的集：原始集（primitiveset）和派生集（derived set）一个原始集是由一些最基本的对象组成的一个派生集是用一个或多个其它集来定义的，也就是说，它的元素来自于其它已存在的集3.2.2模型的集部分集部分在程序中又称为集合段，它是LINGO模型的一个可选部分在LINGO模型中使用集之前，必须在集部分事先定义集部分以关键字“sets:（1）原始集的定义为了定义一个原始集，必须详细说明集的名字，而集的元素和相应的属性是可选的定义一个原始集，用下面的语法：setname/member_list/:

3、attribute_list;注意：用“”表示该部分内容是可选的（下同）Setname是用来标记集的名字，最好具有较强的可读性集名字必须严格符合标准命名规则：以拉丁字母或下划线为首字符，其后由拉丁字母、下划线、阿拉伯数字组成的总长度不超过32个字符的字符串，且不区分大小写该命名规则同样适用于集元素名和属性名等的命名Member_list是集元素的列表如果集元素放在集定义中，那么对它们可采取显式和隐式罗列两种方式如果集元素不放在集定义中，那么可以在随后的数据部分定义 当显式罗列元素时，必须为每个元素输入一个不同的名字，中间用空格或逗号隔开，允许混合使用例3.1 定义一个名为friends的原始集

4、，它具有元素John，Jill，Rose和Mike，其属性有sex和age：sets:friends/John Jill, Rose Mike/: sex, age;endsets 当隐式罗列元素时，不必罗列出每个集元素可采用如下语法：setname/member1.memberN/: attribute_list;这里的member1是集的第一个元素名，memberN是集的最后一个元素名LINGO将自动产生中间的所有元素名LINGO也接受一些特定的首元素名和末元素名，用于创建一些特殊的集 集元素不放在集定义中，而在随后的数据部分来定义例3.2!集部分;friends:sex,age;数据部分

5、;data:friends,sex,age=John,1,16 Jill,0,14 Rose,0,17 Mike,1,13;enddata开头用感叹号（!），末尾用分号（;）表示注释，可跨多行在集部分只定义了一个集friends，并未指定元素在数据部分罗列了集元素John，Jill，Rose和Mike，并对属性sex和age分别给出了值集元素无论用何种字符标记,它的索引都是从1开始连续计数在attribute_ list可以指定一个或多个集元素的属性，属性之间必须用逗号隔开LINGO内置的建模语言是一种描述性语言，用它可以描述现实世界中的一些问题，然后再借助于LINGO求解器求解因此，集属性的

6、值一旦在模型中被确定，就不可能再更改只有在初始部分中给出的集属性值在以后的求解中可更改这与前面并不矛盾，初始部分是LINGO求解器的需要，并不是描述问题所必须的（2） 定义派生集为了定义一个派生集，必须详细说明集的名字和父集的名字，而集元素和属性是可选的可用下面的语法定义一个派生集：setname（parent_set_list）/member_list/:setname是集的名字parent_set_list是已定义的集的列表，多个时要用逗号隔开如果没有指定成员列表，那么LINGO会自动创建父集元素的所有组合作为派生集的元素派生集的父集既可以是原始集，也可以是其它的派生集例3.3produc

7、t/A,B/;machine/M,N/;week/1.2/;allowed（product,machine,week）:x;LINGO生成了三个父集的所有组合共八组作为allowed集的元素，列表如下：编号 元素 1 （A,M,1） 2 （A,M,2） 3 （A,N,1） 4 （A,N,2） 5 （B,M,1） 6 （B,M,2） 7 （B,N,1） 8 （B,N,2）元素列表被忽略时，派生集成员由父集成员所有的组合构成，这样的派生集成为稠密集如果限制派生集的成员，使它成为父集成员所有组合构成的集合的一个子集，这样的派生集成为稀疏集同原始集一样，派生集元素的说明也可以放在数据部分一个派生集的元

8、素列表有两种方式生成：显式罗列；设置元素选择的过滤器当采用方式时，必须显式罗列出所有要包含在派生集中的元素，并且罗列的每个元素要属于稠密集使用前面的例子，显式罗列派生集的元素，如：allowed（product,machine,week）/A M 1,A N 2,B N 1/;如果需要生成一个大的、稀疏的集，那么显式罗列就十分麻烦但是许多稀疏集的元素都满足一些条件，可以把这些逻辑条件看作过滤器，在LINGO生成派生集的元素时把使逻辑条件为假的元素从稠密集中过滤掉例3.4学生集：性别属性sex，1表示男性，0表示女性；年龄属性age;students/John,Jill,Rose,Mike/:男

9、学生和女学生的联系集：友好程度属性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;sex,age =1 16,0 14,0 17,0 13;friend =0.3,0.5,0.6;用竖线（|）来标记一个元素过滤器的开始#eq#是逻辑运算符，用来判断是否“相等”. &1可看作派生集的第1个原始父集的索引，它取遍该原始父集的所有元素；&2可看作派生集的第2 个原

10、始父集的索引，它取遍该原始父集的所有元素；3，&4，依此类推注意如果派生集B的父集是另外的派生集A，那么上面所说的原始父集是集A向前回溯到最终的原始集，其顺序保持不变，并且派生集A的过滤器对派生集B仍然有效因此，派生集的索引个数是最终原始父集的个数，索引的取值是从原始父集到当前派生集所作限制的总和3.3 LINGO数据部分和初始部分在处理模型的数据时，需要为集指定一些元素并且在LINGO求解模型之前为集的某些属性指定数值为此，LINGO为用户提供了两个可选部分：输入集元素数值的数据部分（Data Section）和为决策变量设置初始值的初始部分（Init Section）3.3.1数据部分（1

11、） 数据部分入门数据部分以关键字“data:”开始，“enddata”结束在这里，可以指定集元素和集的属性其语法如下：object_list = value_list;对象列（object_list）包含要指定值的属性名、要设置集元素的集名，用逗号或空格隔开一个对象列中只能有一个集名，而属性名可以有任意多个如果对象列中有多个属性名，那么它们的类型必须一致数值列（value_list）包含要分配给对象列中对象的值，用逗号或空格隔开注意属性值的个数必须等于集元素的个数例3.5SET0/A,B,C/: X,Y;X=1,2,3;Y=4,5,6;在集SET0中定义了两个属性X和YX的三个值是1，2，3，

12、Y的三个值是4，5，6也可采用如下例子中的复合数据说明（data statement）实现同样的功能例3.6X,Y=1 4 2,5 3 6;如果对象列中有n个对象，LINGO在为对象指定值时，首先在n个对象的第1个索引处依次分配数值列中的前n个对象，然后在n个对象的第2个索引处依次分配数值列中紧接着的n个对象，依此类推 （2） 参数输入在数据部分也可以指定一些标量变量（scalar variables）当一个标量变量在数据部分确定时，称之为参数例如，假设模型中用利率9%作为一个参数，就可以输入一个利率作为参数例3.7interest_rate = .09;实际中也可以同时指定多个参数如：int

13、erest_rate,inflation_rate = .09, .025;（3） 实时数据处理在某些情况下，模型中的某些数据并不是定值譬如模型中有一个参数在2%至6%范围内，对不同的值求解模型，观察模型的结果对参数依赖的程度，那么把这种情况称为实时数据处理处理方法是在该语句的数值后面输入一个问号（?）例3.8interest_rate,inflation_rate = .09 ?;在每一次求解模型时，LINGO都会提示为参数inflation_rate输入一个值在WINDOWS操作系统下，将会看到一个如下面的对话框：直接输入一个值再点击OK按钮，LINGO就会把输入的值指定赋给inflati

14、on_rate，然后继续求解模型除了参数之外，也可以实时输入集的属性值，但不允许实时输入集元素名（4） 指定属性为一个值可以在数据定义的右边输入一个值来把所有的元素的该属性指定为一个值如下面的例子例3.9days /MO,TU,WE,TH,FR,SA,SU/:needs;needs = 40;LINGO将用40指定days集的所有元素的needs属性对于多个属性的情形如下：needs,cost;needs cost = 40 90;（5） 数据部分的未知数值表示法有时候只需为一个集的部分元素的某个属性指定数值，而让其余元素的该属性是未知的，以便让LINGO去求出它们的最优值在数据定义中输入两个

15、相连的逗号表示该位置对应元素的属性值未知，两个逗号间可以有空格例3.10years/1.6/: capacity;capacity = ,24,40,;属性capacity的第2个和第3个值分别为24和40，其余的未知3.3.2初始部分初始部分是LINGO提供的另一个可选内容在初始部分中，与数据部分中的数据定义相同，可以输入初始定义（initialization statement）在对实际问题的建模时，初始部分并不起到描述模型的作用，初始部分输入的值仅被LINGO求解器当作初始值来使用，并且仅仅对非线性模型有用这与数据部分指定变量的值不同，LINGO求解器可以自由改变初始部分初始化变量的数值

16、一个初始部分以关键字“init:”开始，以关键字“endinit”结束初始部分的初始定义规则和数据部分的数据定义规则相同也就是说，可以在定义的左边同时初始化多个集属性，即可以把集属性初始化为一个数值，也可以用问号定义为实时数据，还可以用逗号指定为未知数值例3.11init:X,Y = 1,0;endinitY=log（X）;X2+Y2=1;3.4 LINGO函数3.4.1运算符及其优先级LINGO 中的运算符可以分为三类：算数运算符、逻辑运算符和关系运算符（1） 算数运算符算数运算符分为5种： （加法）， （减法）， （乘法）， （除法）， （求幂）（2） 逻辑运算符逻辑运算符分为两类：#AN

17、D#（与），#OR#（或），#NOT#（非）：这3个运算符是参与逻辑值之间的运算，其结果还是逻辑值运算符#EQ#（等于），#NE#（不等于），#GT#（大于），#GE#（大于等于），#LT#（小于），#LE#（小于等于）是用于“数与数之间”的比较，其结果是实逻辑值（3） 关系运算符LINGO中有3种关系运算符：（大于等于），=（等于）注意LINGO中优化模型的约束一般没有严格大于、严格小于，要和逻辑运算符区分开运算符的优先等级如附表3-2所示附表3-2 运算符的优先级优先级运算符高级#NOT# ,（负号）* ,/+,#EQ#,#NE#,#GT#,#GE#,#LT#,#LE#,#AND#,#OR

18、#最低 ,+ ,3.4.2 LINGO数学函数（1） 基本数学函数LINGO中有相当丰富的数学函数，这些函数的用法简单下面列表对各个函数的用法做简单的介绍，具体情况如附表3-3所示（2） 集合循环函数集合循环是指对集合上的元素（下标）进行循环操作的函数，它的一般用法如下：function（setname（set_index_list）|condition:expression_list）；其中function是集合函数名，是FOR,MAX,MIN,PROD,SUM五种之一setname是集合名；set_index_list是集合索引列表（可以省略）；condition是实用逻辑表达式描述的过滤

19、条件（通常含有索引，可以省略）；expression_list是一个表达式（对FOR可以是一组表达式）下面对具体的集合函数作如下解释：FOR（集合元素的循环函数）：对集合setname的每个元素独立生成表达式，表达式由expression_list描述MAX（集合属性的最大值）：返回集合setname上的表达式的最大值MIN（集合属性的最小值） ：返回集合setname上的表达式的最小值PROD（集合元素的乘积函数）：返回集合setname上的表达式的积SUM（集合元素的求和函数） ：返回集合setname上的表达式的和表附3-3 基本数学函数函数调用格式含义ABS（X）返回X的绝对值COS（

20、X）返回X的余弦值（X单位是弧度）SIN（X）返回X的正弦值（X单位是弧度）FLOOR（X）返回X的整数部分LGM（X）返回X的伽马（Gamma）函数的自然对数值LOG（X）返回X的自然对数值MOD（X,Y）返回X对Y取模的结果POW（X,Y）返回XY的值SIGN（X）返回X的符号值EXP（X）返回eX的值SMAX（LIST）返回一列数的最大值SMIN（LIST）返回一列数的最小值SQR（X）返回X的平方SQRT（X）返回X的正的平方根值TAN（X）返回X的正切值 （3） 集合操作函数集合操作函数是对集合进行操作的函数，主要有4种，下面分别介绍它们的一般用法1）INDEX（set_name,p

21、rimitive_set_element）这个函数给出元素primitive_set_element在集合set_name中的索引值（即按定义集合时元素出现顺序的位置编号）如果省略编号set_name，LINGO按模型中定义的集合顺序找到第一个含有元素primitive_set_element的集合，并返回索引值通过下面例子解释函数的使用方法例如，假设定义一个女孩的姓名集合和一个男孩的姓名集合：SETS:GIRLS/DEBBLE,SUE,ALICE/;BOYS/BOB,JOE,SUE,FRED/;ENDSETS注意到女孩集和男孩集中都有一个为SUE的元素，如果要调用此函数INDEX（SUE），

22、则得到返回索引值是2因为集合GIRLS在集合BOYS之前，则索引函数只对集合GIRLS检索如果想查找男孩集中的SUE，则应该使用INDEX（BOYS,SUE），则此时得到的索引值是32）IN（set_name,primitive_index_1,primitive_index_2 ）这个函数用于判断一个集合中是否含有某个索引值它的返回值是1（逻辑值“真”），或是0（逻辑值“假”）例3.12全集为I，B是I的一个子集，C是B的补集I/x1.x4/;B（I）/x2/;C（I）|#not#in（B,&1）:3）wrap（index,limit）该函数返回j=index-k*limit，其中k是一个整

23、数，取适当值保证j落在区间1，limit内该函数相当于index模limit再加1该函数在循环、多阶段计划编制中特别有用4）size（set_name）该函数返回集set_name的元素个数在LINGO模型中，如果没有明确给出集的大小，则使用该函数能够使模型中的数据变化和集的大小改变更加方便（4） 变量定界函数变量界定函数能够实现对变量取值范围的附加限制，共4种：1）bin（x）表示限制就是x为0或1；2）bnd（L,x,U）表示限制变量x满足；3）free（x）表示取消对变量x的默认下界为0的限制，即x可以取任意实数；4）gin（x）表示限制变量x为整数在默认情况下，LINGO规定变量是非负

24、的，即下界值为0，上界为+free取消了默认的下界为0的限制，使变量也可以取负值bnd用于设定一个变量的上下界,它也可以取消默认下界为0的约束（5） 概率论中相关函数1）pbn（p,n,x）二项分布的分布函数，当n和（或）x不是整数时，用线性插值法进行计算2）pcx（n,x）自由度为n的2分布的分布函数在x点的取值3）peb（load,x）当到达负荷（平均服务强度）为load，服务系统有x个服务台，且系统容量无限时的Erlang繁忙概率，多用于解决排队问题4）pel（load,x）当到达负荷（平均服务强度）为load，服务系统有x个服务台，系统容量为有限时的Erlang繁忙概率，多用于解决排队问题5）pfd（n,d,x）自由度为n和d的F分布的分布函数在x点的取值6）pfs（load,x,c）当负荷上限为load，顾客数为c，平行服务台数量为x时，顾客源有限的Poisson服务系统的等待或有返回顾客数的期望值load是顾客数乘以平均服务时间，再除以平均返回时间当c和（或）x不是整数时，采用线性插值进行计算7）phg（pop,g,n,x）超几何（Hyperge