Lingo基本用法总结 本文关键词：用法，Lingo
Lingo基本用法总结 本文简介：Lingo基本用法总结（除集函数部分）LINGO是用来求解线性和非线性优化问题的简易工具。Lingo免费版可以支持30个未知数，lingo破解版可以支持几万个未知数、几万个约束条件。当你在windows下开始运行LINGO系统时，会得到类似下面的一个窗口：外层是主框架窗口，包含了所有菜单命令和工具条
Lingo基本用法总结 本文内容：
Lingo基本用法总结（除集函数部分）
LINGO是用来求解线性和非线性优化问题的简易工具。Lingo免费版可以支持30个未知数，lingo破解版可以支持几万个未知数、几万个约束条件。
当你在windows下开始运行LINGO系统时，会得到类似下面的一个窗口：
外层是主框架窗口，包含了所有菜单命令和工具条，其它所有的窗口将被包含在主窗口之下。在主窗口内的标题为LINGO
Model
–
LINGO1的窗口是LINGO的默认模型窗口，建立的模型都都要在该窗口内编码实现。下面举两个例子。
例1.1
如何在LINGO中求解如下的LP问题：
在模型窗口中输入如下代码：
min=2*x1+3*x2;
x1+x2>=350;
x1>=100;
2*x1+x2=350，x1>=100，2*x1+x2=350;
x1>=100;
2*x1+x2=100，2*x1+x2=350;
x1>=100;
2*x1+x2<=600;
然后点击工具条上的按钮即可。
例1.2
使用LINGO软件计算6个发点8个收点的最小费用运输问题。产销单位运价如下表。
单
位
销地
运
价
产地
B1
B2
B3
B4
B5
B6
B7
B8
产量
A1
6
2
6
7
4
2
5
9
60
A2
4
9
5
3
8
5
8
2
55
A3
5
2
1
9
7
4
3
3
51
A4
7
6
7
3
9
2
7
1
43
A5
2
3
9
5
7
2
6
5
41
A6
5
5
2
2
8
1
4
3
52
销量
35
37
22
32
41
32
43
38
使用LINGO软件，编制程序如下：
model:
!6发点8收点运输问题;
sets:
warehouses/wh1wh6/:
capacity;
vendors/v1v8/:
demand;
links(warehouses,vendors):
cost,volume;
endsets
!目标函数;
[email protected](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建模语言的基础，是程序设计最强有力的基本构件。借助于集，能够用一个单一的、长的、简明的复合公式表示一系列相似的约束，从而可以快速方便地表达规模较大的模型。
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/member1memberN/[:
attribute_list];
这里的member1是集的第一个成员名，memberN是集的最末一个成员名。LINGO将自动产生中间的所有成员名。LINGO也接受一些特定的首成员名和末成员名，用于创建一些特殊的集。列表如下：
隐式成员列表格式
示例
所产生集成员
1n
15
1,2,3,4,5
StringMStringN
Car2car14
Car2,Car3,Car4,…,Car14
DayMDayN
MonFri
Mon,Tue,Wed,Thu,Fri
MonthMMonthN
OctJan
Oct,Nov,Dec,Jan
MonthYearMMonthYearN
Oct2001Jan2002
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/12/;
allowed(product,machine,week):x;
endsets
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)
成员列表被忽略时，派生集成员由父集成员所有的组合构成，这样的派生集成为稠密集。如果限制派生集的成员，使它成为父集成员所有组合构成的集合的一个子集，这样的派生集成为稀疏集。同原始集一样，派生集成员的声明也可以放在数据部分。一个派生集的成员列表有两种方式生成：①显式罗列；②设置成员资格过滤器。当采用方式①时，必须显式罗列出所有要包含在派生集中的成员，并且罗列的每个成员必须属于稠密集。使用前面的例子，显式罗列派生集的成员：
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(
!男学生和女学生的友好程度大于0.5的集;
linkmf2(linkmf)
|
friend(
endsets
data:
sex,age
=
1
16
0
14
0
17
0
13;
friend
=
0.3
0.5
0.6;
enddata
用竖线（|）来标记一个成员资格过滤器的开始。#eq#是逻辑运算符，用来判断是否“相等”，可参考4.
对象列（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:
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
指定属性为一个值
可以在数据声明的右边输入一个值来把所有的成员的该属性指定为一个值。看下面的例子。
例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/15/:
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
[email protected](X);
X^2+Y^2<=1;
好的初始点会减少模型的求解时间。
在这一节中，我们仅带大家接触了一些基本的数据输入和初始化概念，不过现在你应该可以轻松的为自己的模型加入原始数据和初始部分啦。
4
LINGO函数
有了前几节的基础知识，再加上本节的内容，你就能够借助于LINGO建立并求解复杂的优化模型了。
LINGO有9种类型的函数：
1．
基本运算符：包括算术运算符、逻辑运算符和关系运算符
2．
数学函数：三角函数和常规的数学函数
3．
金融函数：LINGO提供的两种金融函数
4．
概率函数：LINGO提供了大量概率相关的函数
5．
变量界定函数：这类函数用来定义变量的取值范围
6．
集操作函数：这类函数为对集的操作提供帮助
7．
集循环函数：遍历集的元素，执行一定的操作的函数
8．
数据输入输出函数：这类函数允许模型和外部数据源相联系，进行数据的输入输出
9．
辅助函数：各种杂类函数
4.1
基本运算符
这些运算符是非常基本的，甚至可以不认为它们是一类函数。事实上，在LINGO中它们是非常重要的。
4.1.1
算术运算符
算术运算符是针对数值进行操作的。LINGO提供了5种二元运算符：
＾
乘方
﹡
乘
／
除
﹢
加
﹣
减
LINGO唯一的一元算术运算符是取反函数“﹣”。
这些运算符的优先级由高到底为：
高
﹣（取反）
＾
﹡／
低
﹢﹣
运算符的运算次序为从左到右按优先级高低来执行。运算的次序可以用圆括号“（）”来改变。
例4.1
算术运算符示例。
2﹣5／3，(2﹢4)／5等等。
4.1.2
逻辑运算符
在LINGO中，逻辑运算符主要用于集循环函数的条件表达式中，来控制在函数中哪些集成员被包含，哪些被排斥。在创建稀疏集时用在成员资格过滤器中。
LINGO具有９种逻辑运算符：
#not#
否定该操作数的逻辑值，＃not＃是一个一元运算符
#eq#
若两个运算数相等，则为true；否则为flase
#ne#
若两个运算符不相等，则为true；否则为flase
#gt#
若左边的运算符严格大于右边的运算符，则为true；否则为flase
#ge#
若左边的运算符大于或等于右边的运算符，则为true；否则为flase
#lt#
若左边的运算符严格小于右边的运算符，则为true；否则为flase
#le#
若左边的运算符小于或等于右边的运算符，则为true；否则为flase
#and#
仅当两个参数都为true时，结果为true；否则为flase
#or#
仅当两个参数都为false时，结果为false；否则为true
这些运算符的优先级由高到低为：
高
#not#
#eq#
#ne#
#gt#
#ge#
#lt#
#le#
低
#and#
#or#
例4.2
逻辑运算符示例
2
#gt#
3
#and#
4
#gt#
2，其结果为假（0）。
4.1.3
关系运算符
在LINGO中，关系运算符主要是被用在模型中，来指定一个表达式的左边是否等于、小于等于、或者大于等于右边，形成模型的一个约束条件。关系运算符与逻辑运算符#eq#、#le#、#ge#截然不同，前者是模型中该关系运算符所指定关系的为真描述，而后者仅仅判断一个该关系是否被满足：满足为真，不满足为假。
LINGO有三种关系运算符：“=”、“=”。LINGO中还能用“”表示大于等于关系。LINGO并不支持严格小于和严格大于关系运算符。然而，如果需要严格小于和严格大于关系，比如让A严格小于B：
A 那么可以把它变成如下的小于等于表达式： A+ε<=B， 这里ε是一个小的正数，它的值依赖于模型中A小于B多少才算不等。 下面给出以上三类操作符的优先级： 高 #not# ﹣（取反） ＾ ﹡ ／ ﹢﹣ #eq# #ne# #gt# #ge# #lt# #le# #and# #or# 低 = 4.2 数学函数 LINGO提供了大量的标准数学函数： @abs(x) 返回x的绝对值 @sin(x) 返回x的正弦值，x采用弧度制 @cos(x) 返回x的余弦值 @tan(x) 返回x的正切值 @exp(x) 返回常数e的x次方 @log(x) 返回x的自然对数 @lgm(x) 返回x的gamma函数的自然对数 @sign(x) 如果x=0时，返回不超过x的最大整数；当x<0时，返回不低于x的最大整数。 @smax(x1,x2,…,xn) 返回x1，x2，…，xn中的最大值 @smin(x1,x2,…,xn) 返回x1，x2，…，xn中的最小值 例4.3 给定一个直角三角形，求包含该三角形的最小正方形。 解：如图所示。 A B C D DA E a b x 求最小的正方形就相当于求如下的最优化问题： LINGO代码如下： model: sets: object/1