LINGO函数.docx
《LINGO函数.docx》由会员分享,可在线阅读,更多相关《LINGO函数.docx(19页珍藏版)》请在冰豆网上搜索。
LINGO函数
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具有9种逻辑运算符:
#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返回-1;否则,返回1
@floor(x)返回x的整数部分。
当x>=0时,返回不超过x的最大整数;当x<0时,返回不低于x的最大整数。
@smax(x1,x2,…,xn)返回x1,x2,…,xn中的最大值
@smin(x1,x2,…,xn)返回x1,x2,…,xn中的最小值
endsets
data:
a,b=3,4;!
两个直角边长,修改很方便;
enddata
f
(1)=a*@sin(x);
f
(2)=b*@cos(x);
f(3)=a*@cos(x)+b*@sin(x);
min=@smax(f
(1),f
(2),f(3));
@bnd(0,x,1.57);
end
在上面的代码中用到了函数@bnd,详情请见4.5节。
4.4概率函数
1.@pbn(p,n,x)
二项分布的累积分布函数。
当n和(或)x不是整数时,用线性插值法进行计算。
2.@pcx(n,x)
自由度为n的χ2分布的累积分布函数。
3.@peb(a,x)
当到达负荷为a,服务系统有x个服务器且允许无穷排队时的Erlang繁忙概率。
4.@pel(a,x)
当到达负荷为a,服务系统有x个服务器且不允许排队时的Erlang繁忙概率。
5.@pfd(n,d,x)
自由度为n和d的F分布的累积分布函数。
6.@pfs(a,x,c)
当负荷上限为a,顾客数为c,平行服务器数量为x时,有限源的Poisson服务系统的等待或返修顾客数的期望值。
a是顾客数乘以平均服务时间,再除以平均返修时间。
当c和(或)x不是整数时,采用线性插值进行计算。
7.@phg(pop,g,n,x)
超几何(Hypergeometric)分布的累积分布函数。
pop表示产品总数,g是正品数。
从所有产品中任意取出n(n≤pop)件。
pop,g,n和x都可以是非整数,这时采用线性插值进行计算。
8.@ppl(a,x)
Poisson分布的线性损失函数,即返回max(0,z-x)的期望值,其中随机变量z服从均值为a的Poisson分布。
9.@pps(a,x)
均值为a的Poisson分布的累积分布函数。
当x不是整数时,采用线性插值进行计算。
10.@psl(x)
单位正态线性损失函数,即返回max(0,z-x)的期望值,其中随机变量z服从标准正态分布。
11.@psn(x)
标准正态分布的累积分布函数。
12.@ptd(n,x)
自由度为n的t分布的累积分布函数。
13.@qrand(seed)
产生服从(0,1)区间的拟随机数。
@qrand只允许在模型的数据部分使用,它将用拟随机数填满集属性。
通常,声明一个m×n的二维表,m表示运行实验的次数,n表示每次实验所需的随机数的个数。
在行内,随机数是独立分布的;在行间,随机数是非常均匀的。
这些随机数是用“分层取样”的方法产生的。
例4.5
model:
data:
M=4;N=2;seed=1234567;
enddata
sets:
rows/1..M/;
cols/1..N/;
table(rows,cols):
x;
endsets
data:
X=@qrand(seed);
enddata
end
如果没有为函数指定种子,那么LINGO将用系统时间构造种子。
14.@rand(seed)
返回0和1间的伪随机数,依赖于指定的种子。
典型用法是U(I+1)=@rand(U(I))。
注意如果seed不变,那么产生的随机数也不变。
例4.6利用@rand产生15个标准正态分布的随机数和自由度为2的t分布的随机数。
model:
!
产生一列正态分布和t分布的随机数;
sets:
series/1..15/:
u,znorm,zt;
endsets
!
第一个均匀分布随机数是任意的;
u
(1)=@rand(.1234);
!
产生其余的均匀分布的随机数;
@for(series(I)|I#GT#1:
u(I)=@rand(u(I-1))
);
@for(series(I):
!
正态分布随机数;
@psn(znorm(I))=u(I);
!
和自由度为2的t分布随机数;
@ptd(2,zt(I))=u(I);
!
ZNORM和ZT可以是负数;
@free(znorm(I));@free(zt(I));
);
end
4.5变量界定函数
变量界定函数实现对变量取值范围的附加限制,共4种:
@bin(x)限制x为0或1
@bnd(L,x,U)限制L≤x≤U
@free(x)取消对变量x的默认下界为0的限制,即x可以取任意实数
@gin(x)限制x为整数
在默认情况下,LINGO规定变量是非负的,也就是说下界为0,上界为+∞。
@free取消了默认的下界为0的限制,使变量也可以取负值。
@bnd用于设定一个变量的上下界,它也可以取消默认下界为0的约束。
4.6集操作函数
LINGO提供了几个函数帮助处理集。
1.@in(set_name,primitive_index_1[,primitive_index_2,…])
如果元素在指定集中,返回1;否则返回0。
例4.7全集为I,B是I的一个子集,C是B的补集。
sets:
I/x1..x4/;
B(I)/x2/;
C(I)|#not#@in(B,&1):
;
endsets
2.@index([set_name,]primitive_set_element)
该函数返回在集set_name中原始集成员primitive_set_element的索引。
如果set_name被忽略,那么LINGO将返回与primitive_set_element匹配的第一个原始集成员的索引。
如果找不到,则产生一个错误。
例4.8如何确定集成员(B,Y)属于派生集S3。
sets:
S1/ABC/;
S2/XYZ/;
S3(S1,S2)/AX,AZ,BY,CX/;
endsets
X=@in(S3,@index(S1,B),@index(S2,Y));
看下面的例子,表明有时为@index指定集是必要的。
例4.9
sets:
girls/debble,sue,alice/;
boys/bob,joe,sue,fred/;
endsets
I1=@index(sue);
I2=@index(boys,sue);
I1的值是2,I2的值是3。
我们建议在使用@index函数时最好指定集。
3.@wrap(index,limit)
该函数返回j=index-k*limit,其中k是一个整数,取适当值保证j落在区间[1,limit]内。
该函数相当于index模limit再加1。
该函数在循环、多阶段计划编制中特别有用。
4.@size(set_name)
该函数返回集set_name的成员个数。
在模型中明确给出集大小时最好使用该函数。
它的使用使模型更加数据中立,集大小改变时也更易维护。
4.7集循环函数
集循环函数遍历整个集进行操作。
其语法为
@function(setname[(set_index_list)[|conditional_qualifier]]:
expression_list);
@function相应于下面罗列的四个集循环函数之一;setname是要遍历的集;set_index_list是集索引列表;conditional_qualifier是用来限制集循环函数的范围,当集循环函数遍历集的每个成员时,LINGO都要对conditional_qualifier进行评价,若结果为真,则对该成员执行@function操作,否则跳过,继续执行下一次循环。
expression_list是被应用到每个集成员的表达式列表,当用的是@for函数时,expression_list可以包含多个表达式,其间用逗号隔开。
这些表达式将被作为约束加到模型中。
当使用其余的三个集循环函数时,expression_list只能有一个表达式。
如果省略set_index_list,那么在expression_list中引用的所有属性的类型都是setname集。
1.@for
该函数用来产生对集成员的约束。
基于建模语言的标量需要显式输入每个约束,不过@for函数允许只输入一个约束,然后LINGO自动产生每个集成员的约束。
例4.10产生序列{1,4,9,16,25}
model:
sets:
number/1..5/:
x;
endsets
@for(number(I):
x(I)=I^2);
end
2.@sum
该函数返回遍历指定的集成员的一个表达式的和。
例4.11求向量[5,1,3,4,6,10]前5个数的和。
model:
data:
N=6;
enddata
sets:
number/1..N/:
x;
endsets
data:
x=5134610;
enddata
s=@sum(number(I)|I#le#5:
x);
end
3.@min和@max
返回指定的集成员的一个表达式的最小值或最大值。
例4.12求向量[5,1,3,4,6,10]前5个数的最小值,后3个数的最大值。
model:
data:
N=6;
enddata
sets:
number/1..N/:
x;
endsets
data:
x=5134610;
enddata
minv=@min(number(I)|I#le#5:
x);
maxv=@max(number(I)|I#ge#N-2:
x);
end
下面看一个稍微复杂一点儿的例子。
例4.13职员时序安排模型一项工作一周7天都需要有人(比如护士工作),每天(周一至周日)所需的最少职员数为20、16、13、16、19、14和12,并要求每个职员一周连续工作5天,试求每周所需最少职员数,并给出安排。
注意这里我们考虑稳定后的情况。
model:
sets:
days/mon..sun/:
required,start;
endsets
data:
!
每天所需的最少职员数;
required=20161316191412;
enddata
!
最小化每周所需职员数;
min=@sum(days:
start);
@for(days(J):
@sum(days(I)|I#le#5:
start(@wrap(J+I+2,7)))>=required(J));
end
计算的部分结果为
Globaloptimalsolutionfoundatiteration:
0
Objectivevalue:
22.00000
VariableValueReducedCost
REQUIRED(MON)20.000000.000000
REQUIRED(TUE)16.000000.000000
REQUIRED(WED)13.000000.000000
REQUIRED(THU)16.000000.000000
REQUIRED(FRI)19.000000.000000
REQUIRED(SAT)14.000000.000000
REQUIRED(SUN)12.000000.000000
START(MON)8.0000000.000000
START(TUE)2.0000000.000000
START(WED)0.0000000.3333333
START(THU)6.0000000.000000
START(FRI)3.0000000.000000
START(SAT)3.0000000.000000
START(SUN)0.0000000.000000
从而解决方案是:
每周最少需要22个职员,周一安排8人,周二安排2人,周三无需安排人,周四安排6人,周五和周六都安排3人,周日无需安排人。
4.8输入和输出函数
输入和输出函数可以把模型和外部数据比如文本文件、数据库和电子表格等连接起来。
1.@file函数
该函数用从外部文件中输入数据,可以放在模型中任何地方。
该函数的语法格式为@file(’filename’)。
这里filename是文件名,可以采用相对路径和绝对路径两种表示方式。
@file函数对同一文件的两种表示方式的处理和对两个不同的文件处理是一样的,这一点必须注意。
例4.14以例1.2来讲解@file函数的用法。
注意到在例1.2的编码中有两处涉及到数据。
第一个地方是集部分的6个warehouses集成员和8个vendors集成员;第二个地方是数据部分的capacity,demand和cost数据。
为了使数据和我们的模型完全分开,我们把它们移到外部的文本文件中。
修改模型代码以便于用@file函数把数据从文本文件中拖到模型中来。
修改后(修改处代码黑体加粗)的模型代码如下:
model:
!
6发点8收点运输问题;
sets:
warehouses/@file('1_2.txt')/:
capacity;
vendors/@file('1_2.txt')/:
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=@file('1_2.txt');
demand=@file('1_2.txt');
cost=@file('1_2.txt');
enddata
end
模型的所有数据来自于1_2.txt文件。
其内容如下:
!
warehouses成员;
WH1WH2WH3WH4WH5WH6~
!
vendors成员;
V1V2V3V4V5V6V7V8~
!
产量;
605551434152~
!
销量;
3537223241324338~
!
单位运输费用矩阵;
62674259
49538582
52197433
76739271
23957265
55228143
把记录结束标记(~)之间的数据文件部分称为记录。
如果数据文件中没有记录结束标记,那么整个文件被看作单个记录。
注意到除了记录结束标记外,模型的文本和数据同它们直接放在模型里是一样的。
我们来看一下在数据文件中的记录结束标记连同模型中@file函数调用是如何工作的。
当在模型中第一次调用@file函数时,LINGO打开数据文件,然后读取第一个记录;第二次调用@file函数时,LINGO读取第二个记录等等。
文件的最后一条记录可以没有记录结束标记,当遇到文件结束标记时,LINGO会读取最后一条记录,然后关闭文件。
如果最后一条记录也有记录结束标记,那么直到LINGO求解完当前模型后才关闭该文件。
如果多个文件保持打开状态,可能就会导致一些问题,因为这会使同时打开的文件总数超过允许同时打开文件的上限16。
当使用@file函数时,可把记录的内容(除了一些记录结束标记外)看作是替代模型中@file(’filename’)位置的文本。
这也就是说,一条记录可以是声明的一部分,整个声明,或一系列声明。
在数据文件中注释被忽略。
注意在LINGO中不允许嵌套调用@file函数。
2.@text函数
该函数被用在数据部分用来把解输出至文本文件中。
它可以输出集成员和集属性值。
其语法为
@text([’filename’])
这里filename是文件名,可以采用相对路径和绝对路径两种表示方式。
如果忽略filename,那么数据就被输出到标准输出设备(大多数情形都是屏幕)。
@text函数仅能出现在模型数据部分的一条语句的左边,右边是集名(用来输出该集的所有成员名)或集属性名(用来输出该集属性的值)。
我们把用接口函数产生输出的数据声明称为输出操作。
输出操作仅当求解器求解完模型后才执行,执行次序取决于其在模型中出现的先后。
例4.15借用例4.12,说明@text的用法。
model:
sets:
days/mon..sun/:
required,start;
endsets
data:
!
每天所需的最少职员数;
required=20161316191412;
@text('d:
\out.txt')=days'至少需要的职员数为'start;
enddata
!
最小化每周所需职员数;
min=@sum(days:
start);
@for(days(J):
@sum(days(I)|I#le#5:
start(@wrap(J+I+2,7)))>=required(J));
end
3.@ole函数
@OLE是从EXCEL中引入或输出数据的接口函数,它是基于传输的OLE技术。
OLE传输直接在内存中传输数据,并不借助于中间文件。
当使用@OLE时,LINGO先装载EXCEL,再通知EXCEL装载指定的电子数据表,最后从电子数据表中获得Ranges。
为了使用OLE函数,必须有EXCEL5及其以上版本。
OLE函数可在数据部分和初始部分引入数据。
@OLE可以同时读集成员和集属性,集成员最好用文本格式,集属性最好用数值格式。
原始集每个集成员需要一个单元(cell),而对于n元的派生集每个集成员需要n个单元,这里第一行的n个单元对应派生集的第一个集成员,第二行的n个单元对应派生集的第二个集成员,依此类推。
@OLE只能读一维或二维的Ranges(在单个的EXCEL工作表(sheet)中),但不能读间断的或三维的Ranges。
Ranges是自左而右、自上而下来读。
例4.16
sets:
PRODUCT;!
产品;
MACHINE;!
机器;
WEEK;!
周;
ALLOWED(PRODUCT,MACHINE,WEEK):
x,y;!
允许组合及属性;
endsets
data:
rate=0.01;
PRODUCT,MACHINE,WEEK,ALLOWED,x,y=@OLE('D:
\IMPORT.XLS'