【数学建模】数独模型Lingo求解程序（精品推荐 免费分享）
本章选择了LINGO建模语言,对数独模型编写了Lingo求解程序。
LINGO是用来求解线性和非线性优化问题的简易工具。LINGO内置了一种建立最优化模型的语言，可以简便地表达大规模问题，利用LINGO高效的求解器可快
速求解并分析结果。
对实际问题建模的时候，总会遇到一或多相联系的对象，比如工厂、消
费者体、交通工具和雇工等等。LINGO允许把这些相联系的对象聚合成集（sets）。一旦把对象聚合成集，就可以利用集来最大限度的发挥LINGO建模语言的优势。
集是LINGO建模语言的基础，是程序设计最强有力的基本构件。借助于集，
能够用一个单一的、长的、简明的复合公式表示一系列相似的约束，从而可以快
速方便地表达规模较大的模型。
LINGO有两种类型的集：原始集(primitive set)和派生集(derived set)。
为了定义一个原始集，必须定义：集的名字、集的成员（可选）和集成员的
属性（可选）， 定义一个原始集，用下面的语法：
setname[/member_list/][:attribute_list];
注意：用“[]”表示该部分内容可选。Setname是你选择的来标记集的名字，
最好具有较强的可读性。集名字必须严格符合标准命名规则：以拉丁字母或
下划线（_）为首字符，其后由拉丁字母（A—Z）、下划线、阿拉伯数字（0，1，„，9）组成的总长度不超过32个字符的字符串，且不区分大小写。
注意：该命名规则同样适用于集成员名和属性名等的命名。
比如number/1..9/;定义一个原始集number，由于在后续的程序中不直接使
用该集，故省略了后面的属性，该集有9个成员，由于成员的规律性，可以采用
隐式罗列法来表示，即不必罗列出每个集成员。具体语法如下：
berN/[: attribute_list];
这里的member1是集的第一个成员名，memberN是集的最末一个成员名。LINGO将自动产生中间的所有成员名。
为了定义一个派生集，必须定义：集的名字、父集的名字、集成员（可选）
和集成员的属性（可选）。定义一个派生集，用下面的语法：
setname(parent_set_list)[/member_list/][:attribute_list];
setname是集的名字。parent_set_list是已定义的集的列表，多个时必须用逗号
隔开。如果没有指定成员列表，那么LINGO会自动创建父集成员的所有组合作为
派生集的成员。派生集的父集既可以是原始集，也可以是其它的派生集。
比如shudu(line,col):a;定义了一个派生集，它没有指定成员列表，那么
LINGO会自动创建父集line和col成员的所有组合作为派生集shudu的成员。A是该派生集的一个属性。
在处理模型的数据时，需要为集指派一些成员并且在LINGO求解模型之前为集的某些属性指定值。为此，LINGO为用户提供了两个可选部分：输入集成员和数
据的数据部分（Data Section）和为决策变量设置初始值的初始部分（Init Section）。比如：
data:
! 键盘输入;
a=？;
enddata
以data:开始，以enddata结束，语句前加了!后成为说明语句。LINGO规定， 每个以分号“；”结束。
LINGO有9种类型的函数：基本运算符、数学函数、金融函数、概率函数、变 量界定函数、集操作函数、集循环函数、数据输入输出函数、辅助函数。
这些运算符是非常基本的，甚至可以不认为它们是一类函数。事实上，在LINGO 中它们是非常重要的。
算术运算符是针对数值进行操作的。LINGO提供了5种二元运算符：＾ 乘方、 ? 乘、／ 除、? 加、? 减。LINGO唯一的一元算术运算符是取反函数“?”。 这些运算符的优先级由高到底为：?（取反）、＾、?／、??。运算符的运算次 序为从左到右按优先级高低来执行。运算的次序可以用圆括号“（）”来改变。
在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#
在LINGO中，关系运算符主要是被用在模型中，来指定一个表达式的左边是否 等于、小于等于、或者大于等于右边，形成模型的一个约束条件。关系运算符与 逻辑运算符#eq#、#l
e#、#ge#截然不同，前者是模型中该关系运算符所指定关系 的为真描述，而后者仅仅判断一个该关系是否被满足：满足为真，不满足为假。
LINGO有三种关系运算符：“=”、“<=”和“>=”。LINGO中还能用“<”表 示小于等于关系，“>”表示大于等于关系。LINGO并不支持严格小于和严格大于关 系运算符。
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中的最小值
LINGO提供了两个金融函数。
1．@fpa(I,n)返回如下情形的净现值：单位时段利率为I，连续n个时段支付， 每个时段支付单位费用。若每个时段支付x单位的费用，则净现值可用x乘 @fpa(I,n)算得。@fpa的计算公式为
,nn11,(1，)I,,kI(1，)Ik,1 。
净现值就是在一定时期内为了获得一定收益在该时期初所支付的实际费用。
2．@fpl(I,n)返回如下情形的净现值：单位时段利率为I，第n个时段支付单
位费用。@fpl(I,n)的计算公式为 最强编程语言
,n(1，I)。
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表示每次实验所需的随机数的个数。在行内，随机数是独立分 布的；在行间，随机数是非常均匀的。这些随机数是用“分层取样”的方法产生 的。
14．@rand(seed)返回0和1间的伪随机数，依赖于指定的种子。典型用法是 U(I+1)=@rand(U(I))。注意如果seed不变，那么产生的随机数也不变。
变量界定函数实现对变量取值范围的附加限制，共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的约束。
LINGO提供了几个函数帮助处理集。
1．@in(set_name,primitive_index_1 [,primitive_index_2,„])如果元素
在指定集中，返回1；否则返回0。
2．@index([set_name,] primitive_set_element)该函数返回在集set_name 中原始集成员primitive_set_element的索引。如果set_name被忽略，那么LINGO 将返回与primitive_set_element匹配的第一个原始集成员的索引。如果不到， 则产生一个错误。

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。