Lingo基本用法总结(除集函数部分) - 图文(4)

2019-04-01 16:30

在数据部分也可以指定一些标量变量(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

16

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;

好的初始点会减少模型的求解时间。 在这一节中,我们仅带大家接触了一些基本的数据输入和初始化概念,不过现在你应该可以轻松的为自己的模型加入原始数据和初始部分啦。

17

§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具有9种逻辑运算符:

#not# 否定该操作数的逻辑值,#not#是一个一元运算符 #eq# 若两个运算数相等,则为true;否则为flase #ne# 若两个运算符不相等,则为true;否则为flase

#gt# 若左边的运算符严格大于右边的运算符,则为true;否则为flase #ge# 若左边的运算符大于或等于右边的运算符,则为true;否则为flase

18

#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中的最大值

19

@smin(x1,x2,?,xn) 返回x1,x2,?,xn中的最小值

例4.3 给定一个直角三角形,求包含该三角形的最小正方形。 解:如图所示。

求最小的正方形就相当于求如下的最优化问题:

C a x B b A

D E

LINGO代码如下:

model: sets:

object/1..3/: f; 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.3 金融函数

目前LINGO提供了两个金融函数。 1.@fpa(I,n)

返回如下情形的净现值:单位时段利率为I,连续n个时段支付,每个时段支付单位费用。若每个时段支付x单位的费用,则净现值可用x乘以@fpa(I,n)算得。@fpa的计算公式为

N??1??n?@fpa(I,N)??1/(1?I)?1??I ????1?I??n?1??N净现值就是在一定时期内为了获得一定收益在该时期初所支付的实际费用。

例4.4 贷款买房问题 贷款金额50000元,贷款年利率5.31%,采取分期付款方式(每年年末还固定金额,直至还清)。问拟贷款10年,每年需偿还多少元? LINGO代码如下:

50000 = x * @fpa(.0531,10);

答案是x=6573.069元。

2.@fpl(I,n)

返回如下情形的净现值:单位时段利率为I,第n个时段支付单位费用。@fpl(I,n)的计算公式为

20


Lingo基本用法总结(除集函数部分) - 图文(4).doc 将本文的Word文档下载到电脑 下载失败或者文档不完整,请联系客服人员解决!

下一篇:企业内部控制存在的问题及对策

相关阅读
本类排行
× 注册会员免费下载(下载后可以自由复制和排版)

马上注册会员

注:下载文档有可能“只有目录或者内容不全”等情况,请下载之前注意辨别,如果您已付费且无法下载或内容有问题,请联系我们协助你处理。
微信: QQ: