资源预览内容
第1页 / 共24页
第2页 / 共24页
第3页 / 共24页
第4页 / 共24页
第5页 / 共24页
第6页 / 共24页
第7页 / 共24页
第8页 / 共24页
第9页 / 共24页
第10页 / 共24页
亲,该文档总共24页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
第2章 数独模型Lingo求解程序 本章信息指数理论性 应用性 趣味性 创新性 本章选择了LINGO建模语言,对数独模型编写了Lingo求解程序。2-1 Lingo简介 LINGO是用来求解线性和非线性优化问题的简易工具。LINGO内置了一种建立最优化模型的语言,可以简便地表达大规模问题,利用LINGO高效的求解器可快速求解并分析结果。 2.1.1 LINGO中的集 对实际问题建模的时候,总会遇到一群或多群相联系的对象,比如工厂、消费者群体、交通工具和雇工等等。LINGO允许把这些相联系的对象聚合成集(sets)。一旦把对象聚合成集,就可以利用集来最大限度的发挥LINGO建模语言的优势。集是LINGO建模语言的基础,是程序设计最强有力的基本构件。借助于集,能够用一个单一的、长的、简明的复合公式表示一系列相似的约束,从而可以快速方便地表达规模较大的模型。LINGO有两种类型的集:原始集(primitive set)和派生集(derived set)。 2.1.1.1 定义原始集 为了定义一个原始集,必须定义:集的名字、集的成员(可选)和集成员的属性(可选), 定义一个原始集,用下面的语法: setname/member_list/:attribute_list; 注意:用“”表示该部分内容可选。Setname是你选择的来标记集的名字, 最好具有较强的可读性。集名字必须严格符合标准命名规则:以拉丁字母或下划线(_)为首字符,其后由拉丁字母(AZ)、下划线、阿拉伯数字(0,1,9)组成的总长度不超过32个字符的字符串,且不区分大小写。 注意:该命名规则同样适用于集成员名和属性名等的命名。 比如number/1.9/;定义一个原始集number,由于在后续的程序中不直接使用该集,故省略了后面的属性,该集有9个成员,由于成员的规律性,可以采用隐式罗列法来表示,即不必罗列出每个集成员。具体语法如下: setname/member1.memberN/: attribute_list; 这里的member1是集的第一个成员名,memberN是集的最末一个成员名。LINGO将自动产生中间的所有成员名。 2.1.1.2 定义派生集 为了定义一个派生集,必须定义:集的名字、父集的名字、集成员(可选)和集成员的属性(可选)。定义一个派生集,用下面的语法: setname(parent_set_list)/member_list/:attribute_list;setname是集的名字。parent_set_list是已定义的集的列表,多个时必须用逗号隔开。如果没有指定成员列表,那么LINGO会自动创建父集成员的所有组合作为派生集的成员。派生集的父集既可以是原始集,也可以是其它的派生集。 比如shudu(line,col):a;定义了一个派生集,它没有指定成员列表,那么LINGO会自动创建父集line和col成员的所有组合作为派生集shudu的成员。A是该派生集的一个属性。 2.1.2 LINGO中的数据 在处理模型的数据时,需要为集指派一些成员并且在LINGO求解模型之前为集的某些属性指定值。为此,LINGO为用户提供了两个可选部分:输入集成员和数据的数据部分(Data Section)和为决策变量设置初始值的初始部分(Init Section)。比如:data:! 键盘输入;a=?;enddata 以data:开始,以enddata结束,语句前加了!后成为说明语句。LINGO规定,每个以分号“;”结束。 2.1.3 LINGO函数 LINGO有9种类型的函数:基本运算符、数学函数、金融函数、概率函数、变量界定函数、集操作函数、集循环函数、数据输入输出函数、辅助函数。 2.1.3.1 基本运算符 这些运算符是非常基本的,甚至可以不认为它们是一类函数。事实上,在LINGO中它们是非常重要的。 2.1.3.1.1 算术运算符 算术运算符是针对数值进行操作的。LINGO提供了5种二元运算符:乘方、乘、除、加、减。LINGO唯一的一元算术运算符是取反函数“”。这些运算符的优先级由高到底为:(取反)、。运算符的运算次序为从左到右按优先级高低来执行。运算的次序可以用圆括号“()”来改变。 2.1.3.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# 2.1.3.1.3 关系运算符 在LINGO中,关系运算符主要是被用在模型中,来指定一个表达式的左边是否等于、小于等于、或者大于等于右边,形成模型的一个约束条件。关系运算符与逻辑运算符#eq#、#le#、#ge#截然不同,前者是模型中该关系运算符所指定关系的为真描述,而后者仅仅判断一个该关系是否被满足:满足为真,不满足为假。 LINGO有三种关系运算符:“=”、“=”。LINGO中还能用“”表示大于等于关系。LINGO并不支持严格小于和严格大于关系运算符。 2.1.3.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的最大整数;当x0时,返回不低于x的最大整数。 smax(x1,x2,xn) 返回x1,x2,xn中的最大值 smin(x1,x2,xn) 返回x1,x2,xn中的最小值2.1.3.3 金融函数 LINGO提供了两个金融函数。 1fpa(I,n)返回如下情形的净现值:单位时段利率为I,连续n个时段支付,每个时段支付单位费用。若每个时段支付x单位的费用,则净现值可用x乘fpa(I,n)算得。fpa的计算公式为 。净现值就是在一定时期内为了获得一定收益在该时期初所支付的实际费用。 2fpl(I,n)返回如下情形的净现值:单位时段利率为I,第n个时段支付单位费用。fpl(I,n)的计算公式为 。 2.1.3.4 概率函数 1pbn(p,n,x)二项分布的累积分布函数。当n和(或)x不是整数时,用线性插值法进行计算。 2pcx(n,x)自由度为n的2分布的累积分布函数。 3peb(a,x)当到达负荷为a,服务系统有x个服务器且允许无穷排队时Erlang繁忙概率。 4pel(a,x)当到达负荷为a,服务系统有x个服务器且不允许排队时的Erlang繁忙概率。 5pfd(n,d,x)自由度为n和d的F分布的累积分布函数。 6pfs(a,x,c)当负荷上限为a,顾客数为c,平行服务器数量为x时,有限源的Poisson服务系统的等待或返修顾客数的期望值。a是顾客数乘以平均服务时间,再除以平均返修时间。当c和(或)x不是整数时,采用线性插值进行计算。 7phg(pop,g,n,x)超几何(Hypergeometric)分布的累积分布函数。pop表示产品总数,g是正品数。从所有产品中任意取出n(npop)件。pop,g,n和x都可以是非整数,这时采用线性插值进行计算。 8ppl(a,x)Poisson分布的线性损失函数,即返回max(0,z-x)的期望值,其中随机变量z服从均值为a的Poisson分布。 9pps(a,x)均值为a的Poisson分布的累积分布函数。当x不是整数时,采用线性插值进行计算。 10psl(x)单位正态线性损失函数,即返回max(0,z-x)的期望值,其中随机变量z服从标准正态分布。 11psn(x)标准正态分布的累积分布函数。 12ptd(n,x)自由度为n的t分布的累积分布函数。 13qrand(seed)产生服从(0,1)区间的拟随机数。qrand只允许在模型的数据部分使用,它将用拟随机数填满集属性。通常,声明一个mn的二维表,m表示运行实验的次数,n表示每次实验所需的随机数的个数。在行内,随机数是独立分布的;在行间,随机数是非常均匀的。这些随机数是用“分层取样”的方法产生的。 14rand(seed)返回0和1间的伪随机数,依赖于指定的种子。典型用法是U(I+1)=rand(U(I)。注意如果seed不变,那么产生的随机数也不变。 2.1.3.5 变量界定函数 变量界定函数实现对变量取值范围的附加限制,共4种: bin(x) 限制x为0或1 bnd(L,x,U) 限制LxU free(x) 取消对变量x的默认下界为0的限制,即x可以取任意实数 gin(x) 限制x为整数在默认情况下,LINGO规定变量是非负的,也就是说下界为0,上界为+。free取消了默认的下界为0的限制,使变量也可以取负值。bnd用于设定一个变量的上下界,它也可以取消默认下界为0的约束。 2.1.3.6集操作函数 LINGO提供了几个函数帮助处理集。 1in
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号