如何在lingo中使用集合1
LINGO教程(精)
![LINGO教程(精)](https://img.taocdn.com/s3/m/991eeae69b89680203d8258d.png)
LINGO中的集集是一群相联系的对象,这些对象也称为集的成员。
一个集可能是一系列产品、卡车或雇员。
每个集成员可能有一个或多个与之有关联的特征,我们把这些特征称为属性。
属性值可以预先给定,也可以是未知的,有待于LINGO求解。
例如,产品集中的每个产品可以有一个价格属性;卡车集中的每辆卡车可以有一个牵引力属性;雇员集中的每位雇员可以有一个薪水属性,也可以有一个生日属性等等。
LINGO有两种类型的集:原始集(primitive set)和派生集(derived set)。
一个原始集是由一些最基本的对象组成的。
一个派生集是用一个或多个其它集来定义的,也就是说,它的成员来自于其它已存在的集。
1 模型的集部分集部分是LINGO模型的一个可选部分。
在LINGO 模型中使用集之前,必须在集部分事先定义。
集部分以关键字“sets:”开始,以“endsets”结束。
一个模型可以没有集部分,或有一个简单的集部分,或有多个集部分。
一个集部分可以放置于模型的任何地方,但是一个集及其属性在模型约束中被引用之前必须定义了它们。
2 定义原始集为了定义一个原始集,必须详细声明:·集的名字·可选,集的成员·可选,集成员的属性定义一个原始集,用下面的语法:setname[/member_list/][:attribute_list];注意:用“[]”表示该部分内容可选。
下同,不再赘述。
Setname是你选择的来标记集的名字,最好具有较强的可读性。
集名字必须严格符合标准命名规则:以拉丁字母或下划线(_)为首字符,其后由拉丁字母(A—Z)、下划线、阿拉伯数字(0,1,…,9)组成的总长度不超过32个字符的字符串,且不区分大小写。
注意:该命名规则同样适用于集成员名和属性名等的命名。
Member_list是集成员列表。
如果集成员放在集定义中,那么对它们可采取显式罗列和隐式罗列两种方式。
如果集成员不放在集定义中,那么可以在随后的数据部分定义它们。
lingo软件使用教程
![lingo软件使用教程](https://img.taocdn.com/s3/m/e4169aaca5e9856a561260b0.png)
lingo软件使用教程一般来说,一个优化模型将由以下三部分组成:1. 目标函数(Objective Function):要达到的目标。
2. 决策变量(Decision variables):每组决策变量的值代表一种方案。
在优化模型中需要确定决策变量的最优值,优化的目标就是找到决策变量的最优值使得目标函数取得最优。
3. 约束条件(Constraints):对于决策变量的一些约束,它限定决策变量可以取的值。
在写数学模型时,一般第一行是目标函数,接下来是约束条件,再接着是一些非负限制等。
在模型窗口输入如下代码:Max = 2*x1+3*x2;X1+2*x2<=8;4*x1<16;4*x2<12;注意:1.每一个lingo表达式最后要跟一个分号;2.多数电脑中没有符号,lingo中<=代替;为了方便可以用<代替小于等于,用>代替大于等于。
3.我们可以添加一些注释,增加程序的可读性。
注释以一个!(叹号必须在英文状态下输入,它会自动变为绿色)开始,以;(分号)结束。
4.Lingo中不区分变量名的大小写。
变量名必须以字母(A-Z)开头,后面的字符可以是字母、数字、下划线。
变量名不能超过32个字符。
Lingo程序的一些规则:1. 在Lingo中最开始都是“MAX=”或者“MIN=”开始表示求目标函数的最大或者最小值。
2. 变量和它前面的系数之间要用“*”连接,中间可以有空格。
3. 变量名不区分大小写,但必须以字母开始,不超过32个字符。
4. 数学表达式结束时要用分号“;”表示结束。
表达式可以写在多行上,但是表达式中间不能用分号。
5. 在电脑系统中一般没有“小于等于”符号,在Lingo采用“<=”来表示“小于等于”,用“>=”表示“大于等于”。
小于等于也可以用更简单的“<”表示,大于等于用“>”表示。
集合段:在我们已经得到的程序里有一些量没有定义,如WAREHOUSES( I),DEMAND( J), LINKS( I, J)。
Lingo的基本使用方法
![Lingo的基本使用方法](https://img.taocdn.com/s3/m/11ca5a9ef61fb7360b4c657d.png)
1. LINGO入门
内容提要
2.在LINGO中使用集合
3. 运算符和函数
4. LINGO的主要菜单命令
5. LINGO命令窗口
6.习题
1. LINGO入门
1. LINGO入门 2.在LINGO中使用集合 3. 运算符和函数 4. LINGO的主要菜单命令 5. LINGO命令窗口 6.习题
LINGO软件的主要特色
运行状态窗口
扩展 的求 解器 (求解 程序) 状态 框
使用的特殊求解程序 : B-and-B (分枝定界算法) Global (全局最优求解程序) Multistart(用多个初始点求解的程序)
目前为止找到的可行 解的最佳目标函数值
有效步数
目标函数值的界
特殊求解程序当前运行步数: 分枝数(对B-and-B程序); 子问题数(对Global程序); 初始点数(对Multistart程序)
INV (0) 10
加上变量的非负约束
注:LINGO中没有数组,只能对每个季度分别定义变量,如正常 产量就要有RP1,RP2,RP3,RP4 4个变量等。写起来就比较麻 烦,尤其是更多(如1000个季度)的时候。
记四个季度组成的集合QUARTERS={1,2,3,4},它们就是上 面数组的下标集合,而数组DEM,RP,OP, INV对集合QUARTERS中 的每个元素1,2,3,4分别对应于一个值。LINGO正是充分利用 了这种数组及其下标的关系,引入了“集合”及其“属性”的概 念,把QUARTERS={1,2,3,4}称为集合,把DEM,RP,OP, INV 称为该集合的属性(即定义在该集合上的属性)。
两种命令模式
Windows模式: 通过下拉式菜单命令驱动LINGO运行(多数
Lingo的基本使用方法
![Lingo的基本使用方法](https://img.taocdn.com/s3/m/4ed9160b581b6bd97f19ea33.png)
(1)
约
(2) 束
(3)
条 件
•无约束优化: 只有目标函数, 没有约束条件 •有约束优化: 目标函数和约束条件都有 •普通方程组:没有(1)(3), 只有(2) •不等式方程组:没有(1)(2), 只有(3)
优化类型
•连续优化 : 全部决策变量取值均 连续数值 (实数) •离散优化 : 部分或全部决策变量 只取离散数值
除“LG4”文件外, 另外几种格式的文件 都是普通的文本文件
.MPS:表示MPS(数学规划系统)格式的模,型可文以件用。任何文本编
辑器打开和编辑。
2.在LINGO中使用集合
1. LINGO入门 2.在LINGO中使用集合 3. 运算符和函数 4. LINGO的主要菜单命令 5. LINGO命令窗口 6.习题
LINGO
黑龙江八一农垦大学 文理学院数学系 高德宝
学习提纲
一. 优化模型与优化软件简介 二. LINGO软件的基本使用方法 三. LINGO程序流程控制和子模型
一. 优化模型与优化软件简介
优化问题的一般形式
min( or max) f ( x ) s.t. hi( x ) 0, i 1,...,m
运行状态窗口
扩展 的求 解器( 求解 程序) 状态 框
使用的特殊求解程序 : B-and-B (分枝定界算法) Global (全局最优求解程序) Multistart(用多个初始点求解的程序
)
目前为止找到的可行 解的最佳目标函数值
有效步数
目标函数值的界
特殊求解程序当前运行步数: 分枝数(对B-and-B程序); 子问题数(对Global程序); 初始点数(对Multistart程序)
Lingo的基本使用方法
![Lingo的基本使用方法](https://img.taocdn.com/s3/m/d83512e0227916888586d7e3.png)
建模时需要注意的几个基本问题
1、尽量使用实数优化,减少整数约束和整数变量 2、尽量使用光滑优化,减少非光滑约束的个数
如:尽量少使用绝对值、符号函数、多个变量求最大/最小值 、四舍五入、取整函数等 3、尽量使用线性模型,减少非线性约束和非线性变量的个数 (如x/y <5 改为x<5y)
4、合理设定变量上下界,尽可能给出变量初始值
目标函数是所有费用的和 约束条件主要有两个:
M I{ 4 N R 0(I) 0 P 4O 5(I) 0 P 2I0 N (I)V } I 1 ,2 ,3 ,4
1)能力限制: R(IP )4,I0 1 ,2 ,3 ,4
2)产品数量的平衡方程:
I( N I ) IV ( N I 1 ) R V ( I ) O P ( I ) D P ( I ) I E , 1 , 2 , 3 , 4 M
5、模型中使用的参数数量级要适当 (如小于103)
第四页,共149页。
1. LINGO入门
内容提要
2.在LINGO中使用集合
3. 运算符和函数
4. LINGO的主要菜单命令
5. LINGO命令窗口
6.习题
第五页,共149页。
1. LINGO入门
1. LINGO入门 2.在LINGO中使用集合 3. 运算符和函数 4. LINGO的主要菜单命令 5. LINGO命令窗口
第二十三页,共149页。
集合及其属性
QUARTERS集合的属性
DEM
RP
OP
INV
1
2
3
4
QUARTERS集合
第二十四页,共149页。
集合元素及集合的属性确定的所有变量
集合QUARTERS 1
LINGO软件的集合功能
![LINGO软件的集合功能](https://img.taocdn.com/s3/m/b81241e1bb4cf7ec4afed042.png)
二、集合循环函数
集合循环函数遍历整个集合进行操作。 其语法为: @function(setname[(set_index_list)[|c onditional_qualifier]]:expression_list); @function相应于下面罗列的五个集合循环 函数之一: @for, @sum, @prod, @min, @max
10
原始集合定义方法
隐式成员列表格式 示例
1..n 1..5
所产生集合成员
1, 2, 3, 4, 5
StringM..StringN
DayM..DayN MonthM..MonthN
Car2..Car14
Mon..Fri Oct..Jan
Car2, Car3, Car4, …, Car14
Mon, Tue, Wed, Thu, Fri Oct, Nov, Dec, Jan
MonthYearM..Month Oct2001..Jan YearN 2002
Oct2001, Nov2001, Dec2001, Jan2002
11
原始集合定义方法
3. 集合成员不放在集合定义中,而在随后的数 据部分来定义。例如:
12
原始集合定义方法
注意: 开头用感叹号(!),末尾用分号(;)表示 注释,可跨多行。 在集合部分只定义了一个集合students, 并未指定成员。在数据部分罗列了集合成员 John、Jill、Rose和Mike,并对属性sex和 age分别给出了值。 集合成员无论用何种字符标记,它的索引 都是从1开始连续计数。在attribute_ list可以 指定一个或多个集合成员的属性,属性之间必须 用逗号隔开。 13
派生集合定义方法
LINGO使用说明(比较简单)
![LINGO使用说明(比较简单)](https://img.taocdn.com/s3/m/c0f206565acfa1c7aa00cce3.png)
Lingo介绍Lingo是美国LINDO系统公司(Lindo Symtem Inc)开发的求解数学规划系列软件中的一个(其他软件为LINGDO,GINO,What’s Best等),它的主要功能是求解大型线性、非线性和整数规划问题,目前的版本是lingo11.0。
lingo分为Demo、solve suite、hyper、industrial、extended等六类不同版本,只有Demo版本是免费的,其他版本需要向LINDO系统公司(在中国的代理商)购买,Lingo的不同版本对模型的变量总数、非线性变量个数、整型变量个数和约束条件的数量做出不同的限制(其中extended版本无限制)。
Lingo的主要功能特色为:(1)既能求解线性规划,也有较强的求解非线性规划的能力;(2)输入模型简练直观;(3)运行速度快、计算能力强;(4)内置建模语言,提供几十种内部函数,从而能以较少语句,较直观的方式描述较大规模的优化模型;(5)将集合的概念引入编程语言,很容易将实际问题转换为Lingo语言;(6)能方便地与excel、数据库等其他软件交换数据。
学校图书馆40本《lingo和excel在数学建模中的应用》,袁新生、邵大宏、郁时炼主编,科学出版社Lingo程序设计简要说明在数学建模中会遇到如规划类的题型,在这种模型中总存在着一个目标,并希望这个目标的取值尽可能的大或小,同时与这个目标有关的一系列变量之间存在一些约束。
在构造出目标函数和约束条件的表达式后,我们需要对求出这个最值和各变量的取值。
一般我们用LINGO来对模型进行求解,本文将通过举一个简单的例子,围绕这个例子逐步学习LINGO 的使用。
LINGO只是一个求解工具,我们主要的任务还是模型的建立!当你在windows下开始运行LINGO系统时,会得到类似下面的一个窗口:外层是主框架窗口,包含了所有菜单命令和工具条,其它所有的窗口将被包含在主窗口之下。
在主窗口内的标题为LINGO Model –LINGO1的窗口是LINGO的默认模型窗口,建立的模型都都要在该窗口内编码实现。
lingo入门教程 --- 集合运用
![lingo入门教程 --- 集合运用](https://img.taocdn.com/s3/m/1b419cb2ccbff121dc368378.png)
lingo入门教程之二--- 集合运用lingo中的集合用法很多,这里主要通过几个例题来进行讲解对于每一个问题,都要先找到对应的目标函数,然后对相应值进行初始化,然后找到约束条件等进行求解例1:SAILCO公司需要决定下四个季度的帆船生产量。
下四个季度的帆船需求量分别是40条,60条,75条,25条,这些需求必须按时满足。
每个季度正常的生产能力是40条帆船,每条船的生产费用为400美元。
如果加班生产,每条船的生产费用为450美元。
每个季度末,每条船的库存费用为20美元。
假定生产提前期为0,初始库存为10条船。
如何安排生产可使总费用最小?分析:用DEM,RP,OP,INV分别表示需求量、正常生产的产量、加班生产的产量、库存量,则DEM,RP,OP,INV对每个季度都应该有一个对应的值,也就说他们都应该是一个由4个元素组成的数组,其中DEM是已知的,而RP,OP,INV是未知数接下里这里例子会讲到关于集合的派生问题,这个跟c++里面的继承与派生比较相像例2:建筑工地的位置(用平面坐标a,b表示,距离单位:公里)及水泥日用量d(吨)下表给出。
有两个临时料场位于P (5,1), Q (2, 7),日储量各有20吨。
从A, B两料场分别向各工地运送多少吨水泥,使总的吨公里数最小。
两个新的料场应建在何处,节省的吨公里数有多大?例3:(最短路问题) 在纵横交错的公路网中,货车司机希望找到一条从一个城市到另一个城市的最短路.下图表示的是公路网,节点表示货车可以停靠的城市,弧上的权表示两个城市之间的距离(百公里).那么,货车从城市S出发到达城市T,如何选择行驶路线,使所经过的路程最短?分析:假设从S到T的最优行驶路线P 经过城市C1, 则P中从S到C1的子路也一定是从S到C1的最优行驶路线;假设P 经过城市C2, 则P中从S到C2的子路也一定是从S到C2的最优行驶路线.因此, 为得到从S到T的最优行驶路线, 只需要先求出从S到Ck(k=1,2)的最优行驶路线,就可以方便地得到从S到T的最优行驶路线. 同样,为了求出从S到Ck(k=1,2)的最优行驶路线, 只需要先求出从S到Bj(j=1,2)的最优行驶路线;为了求出从S到Bj(j=1,2)的最优行驶路线, 只需要先求出从S到Ai(i=1,2,3)的最优行驶路线. 而S到Ai(i=1,2,3)的最优行驶路线是很容易得到的(实际上, 此例中S到Ai(i=1,2,3)只有唯一的道路) .此例中可把从S到T的行驶过程分成4个阶段,即S→Ai(i=1,2或3),Ai→Bj(j=1或2),Bj→Ck(k=1或2),Ck→T. 记d(Y,X)为城市Y与城市X之间的直接距离(若这两个城市之间没有道路直接相连,则可以认为直接距离为∞),用L(X)表示城市S到城市X的最优行驶路线的路长:。
第4讲LINGO基本使用方法续
![第4讲LINGO基本使用方法续](https://img.taocdn.com/s3/m/b56ea8af71fe910ef12df836.png)
优 优化建模 化建 优 化 建模 模
例2 某班8名同学准备分成4个调查队(每队两人)前往4个 地区进行社会调查。这8名同学两两之间组队的效率如 下表所示(由于对称性,只列出了严格上三角部分),问 如何组队可以使总效率最高?
学生 S1 S2 S3 S1 S2 S3 S4 S5 S6 S7 S8
56 6 8 A 2 3 6 7 3 A3 4
A1
优 优化建模 化建 优 化 建模 模
B1
7 B2
6 8 9
C1 C2
5 6
T
L B1 min L A1 6, L A2 8, L A3 7 10 L A3 7; L B2 min L A1 5, L A2 6, L A3 4 7 L A3 4; L C1 min L B1 6, L B2 8 15 L B2 8; L C2 min L B1 7, L B2 9 16 L B2 9; L T min L C1 5, L C2 6 20 L C1 5.
集合的使用小结
集合的不同类型及其关系 集合 派生集合 基本集合
稀疏集合
稠密集合
元素列表法 元素过滤法 直接列举法
2019/4/6
隐式列举法
上一页 下一页
新余学院
建模组
Xinyu University MCM
优 优化建模 化建 优 化 建模 模
基本集合的定义语法
基本集合的定义格式为(方括号“[ ]”中的内容是可选项, 可以没 有): setname [/member_list/] [: attribute_list]; 其中setname为定义的集合名,member_list为元素列表, attribute_list为属性列表。元素列表可以采用显式列举法(即直接 将所有元素全部列出,元素之间用逗号或空格分开),也可以采用隐 式列举法。隐式列举法可以有几种不同格式, 类型 隐式列举格式 示例 示例集合表示的元素
Lingo软件使用指南
![Lingo软件使用指南](https://img.taocdn.com/s3/m/fa9254c5192e45361166f598.png)
Lingo软件使用指南摘要:本文介绍了Lingo软件的基本使用方法。
从最基本的使用到复杂问题的解决,本文给出了比较详细的介绍。
Lingo软件是美国Lindo公司的产品,主要用来求解优化问题。
它是一个非常强大的软件,可以求解大部分优化问题,包括线性规划、二次规划、整数规划、运输问题等,是目前全球应用最广泛的优化软件之一。
这里我们简单介绍它的使用方法。
一进入Lingo如果你的计算机已经安装了Lingo,只需要在桌面上双击Lingo的快捷方式,就可以进入Lingo。
为了使自己的程序易于阅读,经常需要有一些注释,因此在编写程序中,每一行前面有感叹号的表示这一行是注释行,在程序运行中不起作用,希望初学者养成注释的好习惯。
二建立数学模型和 Lingo模型语言例1 在Lingo的命令窗口中输入下面的线性规划模型!目标函数;MAX = 100 * x1 + 150 * x2;!第一个约束;X1<= 100;!第二个约束;X2 <= 120;!第三个约束;X1 + 2 * x2<= 160;!end可有可无;end求解可得全局最优解:Objective value: 14500.00Variable ValueX1 100.0000X2 30.00000从这个例子可以看出,用Lingo软件求解一个简单的优化问题是非常容易的。
我们只需要输入优化问题的两个主要部分:目标函数和约束,就可以直接求解。
对于比较简单的问题,我们可以采取这种直接的方式去求解,但是,对于比较复杂的问题,用这种方式就不现实。
比如下面的例2,这就必须要使用Lingo的模型语言。
例2 一个运输问题假设WWW公司有6个仓库,储存着8个分厂生产所需要的原材料。
要求每一个仓库的供应量不能超过储存量,而且每一个分厂的需求必须得到满足。
问:如何组织运输,使总运输费用最小?已知从6个仓库到8个分厂的运输费用表。
表1 供应表2 需求表3 运输费用Wh5 2 3 9 5 7 2 6 5Wh6 5 5 2 2 8 1 4 3 这个问题是一个典型的优化问题,通常称为运输问题。
LINGO的使用方法说明大全
![LINGO的使用方法说明大全](https://img.taocdn.com/s3/m/5769f22e866fb84ae45c8d50.png)
LINGO的使用简介LINGO软件是美国的LINGO系统公司开发的一套专门用于求解最优化问题的软件包.LINGO除了能够用于求解线性规划和二次规划外,还可以用于非线性规划求解、以及一些线性和非线性方程(组)的求解等.LINGO软件的最大特色在于它允许优化模型中的决策变量为整数,即可以求解整数规划,而且执行速度快.LINGO是用来求解线性和非线性优化问题的简易工具.LINGO内置了一种建立最优化模型的语言,可以简便地表达大规模问题,利用LINGO高效的求解器可快速求解并分析结果.在这里仅简单介绍LINGO的使用方法.LINGO(Linear INteractive and General Optimizer )的基本含义是交互式的线性和通过优化求解器.它是美国芝加哥大学的 Linus Schrage 教授于1980年开发了一套用于求解最优化问题的工具包,后来经过完善成何扩充,并成立了LINDO系统公司.这套软件主要产品有:LINDO,LINGO,LINDO API和What’sBest.它们在求解最优化问题上,与同类软件相比有着绝对的优势.软件有演示版和正式版.正式版包括:求解包(solver suite)、高级版(super)、超级版(hyper)、工业版(industrial)、扩展版(extended).不同版本的LINGO对求解问题的规模有限制,如附表3-1所示.附表3-1 不同版本LINGO对求解规模的限制版本类型总变量数整数变量数非线性变量数约束数演示版 300 30 30 150求解包 500 50 50 250高级版 2000 200 200 1000超级版 8000 800 800 4000工业版 32000 3200 32000 16000扩展版无限无限无限无限3.1 LINGO程序框架LINGO可以求解线性规划、二次规划、非线性规划、整数规划、图论及网络最优化问题和最大最小求解问题,以及排队论模型中最优化等问题.一个LINGO程序一般会包括以下几个部分:(1) 集合段:集部分是LINGO模型的一个可选部分.在LINGO模型中使用集之前,必须在集部分事先定义.集部分以关键字“sets:”开始,以“endsets”结束.一个模型可以没有集部分,或有一个简单的集部分,或有多个集部分.一个集部分可以放置于模型的任何地方,但是一个集及其属性在模型约束中被引用之前必须先定义.(2) 数据段:在处理模型的数据时,需要为集部分定义的某些元素在LINGO求解模型之前为其指定值.数据部分以关键字“data:”开始,以关键字“enddata”结束.(3) 目标和约束段:这部分用来定义目标函数和约束条件等.该部分没有开始和结束的标记.主要是要用到LINGO的内部函数,尤其是与集合有关的求和与循环函数等.(4)初始段:这个部分要以关键字“INIT:”开始,以关键字“ENDINIT”结束,它的作用是对集合的属性定义一个初值.在一般的迭代算法中,如果可以给一个接近最优解的初始值,会大大减少程序运行的时间.(5) 数据预处理段:这一部分是以关键字“CALC:”开始,以关键字“ENDCALC”结束.它的作用是把原始数据处理成程序模型需要的数据,它的处理是在数据段输入完以后、开始正式求解模型之前进行的,程序语句是按顺序执行的.3.2 LINGO中集合的概念在对实际问题建模的时候,总会遇到一群或多群相联系的对象,比如工厂、消费者群体、交通工具和雇工等等.LINGO允许把这些相联系的对象聚合成集(sets).一旦把对象聚合成集,就可以利用集来最大限度地发挥LINGO建模语言的优势.现在将深入介绍如何创建集,并用数据初始化集的属性.3.2.1集的构成集是LINGO建模语言的基础,是程序设计最强有力的基本构件.借助于集能够用一个单一的、简明的复合公式表示一系列相似的约束,从而可以快速方便地表达规模较大的模型.集是一群相联系的对象,这些对象也称为集的元素.一个集可能是一系列产品、卡车或雇员.每个集的元素可能有一个或多个与之有关联的特征,把这些特征称为属性.属性值可以预先给定,也可以是未知的,有待于LINGO求解的.LINGO有两种类型的集:原始集(primitive set)和派生集(derived set).一个原始集是由一些最基本的对象组成的.一个派生集是用一个或多个其它集来定义的,也就是说,它的元素来自于其它已存在的集.3.2.2模型的集部分集部分在程序中又称为集合段,它是LINGO模型的一个可选部分.在LINGO模型中使用集之前,必须在集部分事先定义.集部分以关键字“sets:”开始,以“endsets”结束.一个模型可以没有集部分,或有一个简单的集部分,或有多个集部分.一个集部分可以放置于模型的任何地方,但是一个集及其属性在模型约束中被引用之前必须先定义.(1)原始集的定义为了定义一个原始集,必须详细说明集的名字,而集的元素和相应的属性是可选的.定义一个原始集,用下面的语法:setname[/member_list/][:attribute_list];注意:用“[]”表示该部分内容是可选的(下同).Setname是用来标记集的名字,最好具有较强的可读性.集名字必须严格符合标准命名规则:以拉丁字母或下划线为首字符,其后由拉丁字母、下划线、阿拉伯数字组成的总长度不超过32个字符的字符串,且不区分大小写.注意:该命名规则同样适用于集元素名和属性名等的命名.Member_list是集元素的列表.如果集元素放在集定义中,那么对它们可采取显式和隐式罗列两种方式.如果集元素不放在集定义中,那么可以在随后的数据部分定义.①当显式罗列元素时,必须为每个元素输入一个不同的名字,中间用空格或逗号隔开,允许混合使用.例3.1 定义一个名为friends的原始集,它具有元素John,Jill,Rose和Mike,其属性有sex和age:sets:friends/John Jill, Rose Mike/: sex, age;endsets②当隐式罗列元素时,不必罗列出每个集元素.可采用如下语法:setname/member1..member N/[: attribute_list];这里的member1是集的第一个元素名,member N是集的最后一个元素名.LINGO将自动产生中间的所有元素名.LINGO也接受一些特定的首元素名和末元素名,用于创建一些特殊的集.③集元素不放在集定义中,而在随后的数据部分来定义.例3.2!集部分;sets:friends:sex,age;endsets!数据部分;data:friends,sex,age=John,1,16 Jill,0,14 Rose,0,17 Mike,1,13;enddata注意:开头用感叹号(!),末尾用分号(;)表示注释,可跨多行.在集部分只定义了一个集friends,并未指定元素.在数据部分罗列了集元素John,Jill,Rose和Mike,并对属性sex和age分别给出了值.集元素无论用何种字符标记,它的索引都是从1开始连续计数.在attribute_ list可以指定一个或多个集元素的属性,属性之间必须用逗号隔开.LINGO内置的建模语言是一种描述性语言,用它可以描述现实世界中的一些问题,然后再借助于LINGO 求解器求解.因此,集属性的值一旦在模型中被确定,就不可能再更改.只有在初始部分中给出的集属性值在以后的求解中可更改.这与前面并不矛盾,初始部分是LINGO求解器的需要,并不是描述问题所必须的.(2) 定义派生集为了定义一个派生集,必须详细说明集的名字和父集的名字,而集元素和属性是可选的.可用下面的语法定义一个派生集:setname(parent_set_list)[/member_list/][:attribute_list];setname是集的名字.parent_set_list是已定义的集的列表,多个时要用逗号隔开.如果没有指定成员列表,那么LINGO会自动创建父集元素的所有组合作为派生集的元素.派生集的父集既可以是原始集,也可以是其它的派生集.例3.3sets:product/A,B/;machine/M,N/;week/1..2/;allowed(product,machine,week):x;endsetsLINGO生成了三个父集的所有组合共八组作为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生成派生集的元素时把使逻辑条件为假的元素从稠密集中过滤掉.例3.4sets:!学生集:性别属性sex,1表示男性,0表示女性;年龄属性age;students/John,Jill,Rose,Mike/:sex,age;!男学生和女学生的联系集:友好程度属性friend![0,1]之间的数;linkmf(students,students)|sex(&1)#eq#1#and#sex(&2)#eq#0: friend;!男学生和女学生的友好程度大于0.5的集;linkmf2(linkmf) | friend(&1,&2) #ge# 0.5 : x;data:sex,age =1 16,0 14,0 17,0 13;friend =0.3,0.5,0.6;enddata用竖线(|)来标记一个元素过滤器的开始.#eq#是逻辑运算符,用来判断是否“相等”. &1可看作派生集的第1个原始父集的索引,它取遍该原始父集的所有元素;&2可看作派生集的第2 个原始父集的索引,它取遍该原始父集的所有元素;&3,&4,…,依此类推.注意如果派生集B的父集是另外的派生集A,那么上面所说的原始父集是集A向前回溯到最终的原始集,其顺序保持不变,并且派生集A的过滤器对派生集B仍然有效.因此,派生集的索引个数是最终原始父集的个数,索引的取值是从原始父集到当前派生集所作限制的总和.3.3 LINGO数据部分和初始部分在处理模型的数据时,需要为集指定一些元素并且在LINGO求解模型之前为集的某些属性指定数值.为此,LINGO为用户提供了两个可选部分:输入集元素数值的数据部分(Data Section)和为决策变量设置初始值的初始部分(Init Section).3.3.1数据部分(1) 数据部分入门数据部分以关键字“data:”开始,“enddata”结束.在这里,可以指定集元素和集的属性.其语法如下:object_list = value_list;对象列(object_list)包含要指定值的属性名、要设置集元素的集名,用逗号或空格隔开.一个对象列中只能有一个集名,而属性名可以有任意多个.如果对象列中有多个属性名,那么它们的类型必须一致.数值列(value_list)包含要分配给对象列中对象的值,用逗号或空格隔开.注意属性值的个数必须等于集元素的个数.例3.5sets:SET0/A,B,C/: X,Y;endsetsdata:X=1,2,3;Y=4,5,6;enddata在集SET0中定义了两个属性X和Y.X的三个值是1,2,3,Y的三个值是4,5,6.也可采用如下例子中的复合数据说明(data statement)实现同样的功能.sets:SET0/A,B,C/: X,Y;endsetsdata:X,Y=1 4 2,5 3 6;enddata如果对象列中有n个对象,LINGO在为对象指定值时,首先在n个对象的第1个索引处依次分配数值列中的前n个对象,然后在n个对象的第2个索引处依次分配数值列中紧接着的n个对象,…,依此类推.(2) 参数输入在数据部分也可以指定一些标量变量(scalar variables).当一个标量变量在数据部分确定时,称之为参数.例如,假设模型中用利率9%作为一个参数,就可以输入一个利率作为参数.例3.7data:interest_rate = .09;enddata实际中也可以同时指定多个参数.如:data:interest_rate,inflation_rate = .09, .025;enddata(3) 实时数据处理在某些情况下,模型中的某些数据并不是定值.譬如模型中有一个参数在2%至6%范围内,对不同的值求解模型,观察模型的结果对参数依赖的程度,那么把这种情况称为实时数据处理.处理方法是在该语句的数值后面输入一个问号(?).例3.8data:interest_rate,inflation_rate = .09 ?;enddata在每一次求解模型时,LINGO都会提示为参数inflation_rate输入一个值.在WINDOWS操作系统下,将会看到一个如下面的对话框:直接输入一个值再点击OK按钮,LINGO就会把输入的值指定赋给inflation_rate,然后继续求解模型.除了参数之外,也可以实时输入集的属性值,但不允许实时输入集元素名.(4) 指定属性为一个值可以在数据定义的右边输入一个值来把所有的元素的该属性指定为一个值.如下面的例子.sets:days /MO,TU,WE,TH,FR,SA,SU/:needs;endsetsdata:needs = 40;enddataLINGO将用40指定days集的所有元素的needs属性.对于多个属性的情形如下:sets:days /MO,TU,WE,TH,FR,SA,SU/:needs,cost;endsetsdata:needs cost = 40 90;enddata(5) 数据部分的未知数值表示法有时候只需为一个集的部分元素的某个属性指定数值,而让其余元素的该属性是未知的,以便让LINGO 去求出它们的最优值.在数据定义中输入两个相连的逗号表示该位置对应元素的属性值未知,两个逗号间可以有空格.例3.10sets:years/1..6/: capacity;endsetsdata:capacity = ,24,40,,,;enddata属性capacity的第2个和第3个值分别为24和40,其余的未知.3.3.2初始部分初始部分是LINGO提供的另一个可选内容.在初始部分中,与数据部分中的数据定义相同,可以输入初始定义(initialization statement).在对实际问题的建模时,初始部分并不起到描述模型的作用,初始部分输入的值仅被LINGO求解器当作初始值来使用,并且仅仅对非线性模型有用.这与数据部分指定变量的值不同,LINGO求解器可以自由改变初始部分初始化变量的数值.一个初始部分以关键字“init:”开始,以关键字“endinit”结束.初始部分的初始定义规则和数据部分的数据定义规则相同.也就是说,可以在定义的左边同时初始化多个集属性,即可以把集属性初始化为一个数值,也可以用问号定义为实时数据,还可以用逗号指定为未知数值.例3.11X,Y = 1,0;endinitY=@log(X);X^2+Y^2<=1;3.4 LINGO函数3.4.1运算符及其优先级LINGO 中的运算符可以分为三类:算数运算符、逻辑运算符和关系运算符.(1) 算数运算符算数运算符分为5种: (加法), (减法), (乘法), (除法), (求幂).(2) 逻辑运算符逻辑运算符分为两类:#AND#(与),#OR#(或),#NOT#(非):这3个运算符是参与逻辑值之间的运算,其结果还是逻辑值.运算符#EQ#(等于),#NE#(不等于),#GT#(大于),#GE#(大于等于),#LT#(小于),#LE#(小于等于)是用于“数与数之间”的比较,其结果是实逻辑值.(3) 关系运算符LINGO中有3种关系运算符:<(小于等于),>(大于等于),=(等于).注意LINGO中优化模型的约束一般没有严格大于、严格小于,要和逻辑运算符区分开.运算符的优先等级如附表3-2所示.附表3-2 运算符的优先级优先级运算符高级#NOT# ,-(负号)^* ,/+,-#EQ#,#NE#,#GT#,#GE#,#LT#,#LE#, #AND#,#OR#最低< ,+ ,>3.4.2 LINGO数学函数(1) 基本数学函数LINGO中有相当丰富的数学函数,这些函数的用法简单.下面列表对各个函数的用法做简单的介绍,具体情况如附表3-3所示.(2) 集合循环函数集合循环是指对集合上的元素(下标)进行循环操作的函数,它的一般用法如下:@function(setname[(set_index_list)[|condition]]:expression_list);其中function是集合函数名,是FOR,MAX,MIN,PROD,SUM五种之一.setname是集合名;set_index_list 是集合索引列表(可以省略);condition是实用逻辑表达式描述的过滤条件(通常含有索引,可以省略);expression_list是一个表达式(对@FOR可以是一组表达式).下面对具体的集合函数作如下解释:@FOR(集合元素的循环函数):对集合setname的每个元素独立生成表达式,表达式由expression_list 描述.@MAX(集合属性的最大值):返回集合setname上的表达式的最大值.@MIN(集合属性的最小值) :返回集合setname上的表达式的最小值.@PROD(集合元素的乘积函数):返回集合setname上的表达式的积.@SUM(集合元素的求和函数) :返回集合setname上的表达式的和.表附3-3 基本数学函数函数调用格式含义@ABS(X) 返回X的绝对值@COS(X) 返回X的余弦值(X单位是弧度)@SIN(X) 返回X的正弦值(X单位是弧度)@FLOOR(X) 返回X的整数部分@LGM(X) 返回X的伽马(Gamma)函数的自然对数值@LOG(X) 返回X的自然对数值@MOD(X,Y) 返回X对Y取模的结果@POW(X,Y) 返回X Y的值@SIGN(X) 返回X的符号值@EXP(X) 返回e X的值@SMAX(LIST) 返回一列数的最大值@SMIN(LIST) 返回一列数的最小值@SQR(X) 返回X的平方@SQRT(X) 返回X的正的平方根值@TAN(X) 返回X的正切值(3) 集合操作函数集合操作函数是对集合进行操作的函数,主要有4种,下面分别介绍它们的一般用法.1)@INDEX([set_name,]primitive_set_element)这个函数给出元素primitive_set_element在集合set_name中的索引值(即按定义集合时元素出现顺序的位置编号).如果省略编号set_name,LINGO按模型中定义的集合顺序找到第一个含有元素primitive_set_element的集合,并返回索引值.通过下面例子解释函数的使用方法.例如,假设定义一个女孩的姓名集合和一个男孩的姓名集合:SETS:GIRLS/DEBBLE,SUE,ALICE/;BOYS/BOB,JOE,SUE,FRED/;ENDSETS注意到女孩集和男孩集中都有一个为SUE的元素,如果要调用此函数@INDEX(SUE),则得到返回索引值是2.因为集合GIRLS在集合BOYS之前,则索引函数只对集合GIRLS检索.如果想查找男孩集中的SUE,则应该使用@INDEX(BOYS,SUE),则此时得到的索引值是3.2)@IN(set_name,primitive_index_1[,primitive_index_2 …])这个函数用于判断一个集合中是否含有某个索引值.它的返回值是1(逻辑值“真”),或是0(逻辑值“假”).例3.12全集为I,B是I的一个子集,C是B的补集.sets:I/x1..x4/;B(I)/x2/;C(I)|#not#@in(B,&1):;endsets3)@wrap(index,limit)该函数返回j=index-k*limit,其中k是一个整数,取适当值保证j落在区间[1,limit]内.该函数相当于index模limit再加1.该函数在循环、多阶段计划编制中特别有用.4)@size(set_name)该函数返回集set_name的元素个数.在LINGO模型中,如果没有明确给出集的大小,则使用该函数能够使模型中的数据变化和集的大小改变更加方便.(4) 变量定界函数变量界定函数能够实现对变量取值范围的附加限制,共4种:1)@bin(x)表示限制就是x为0或1;2)@bnd(L,x,U)表示限制变量x满足;3)@free(x)表示取消对变量x的默认下界为0的限制,即x可以取任意实数;4)@gin(x)表示限制变量x为整数.在默认情况下,LINGO规定变量是非负的,即下界值为0,上界为+∞.@free取消了默认的下界为0的限制,使变量也可以取负值.@bnd用于设定一个变量的上下界,它也可以取消默认下界为0的约束.(5) 概率论中相关函数1)@pbn(p,n,x)二项分布的分布函数,当n和(或)x不是整数时,用线性插值法进行计算.2)@pcx(n,x)自由度为n的χ2分布的分布函数在x点的取值.3)@peb(load,x)当到达负荷(平均服务强度)为load,服务系统有x个服务台,且系统容量无限时的Erlang繁忙概率,多用于解决排队问题.4)@pel(load,x)当到达负荷(平均服务强度)为load,服务系统有x个服务台,系统容量为有限时的Erlang繁忙概率,多用于解决排队问题.5)@pfd(n,d,x)自由度为n和d的F分布的分布函数在x点的取值.6)@pfs(load,x,c)当负荷上限为load,顾客数为c,平行服务台数量为x时,顾客源有限的Poisson服务系统的等待或有返回顾客数的期望值.load是顾客数乘以平均服务时间,再除以平均返回时间.当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点的取值.当x不是整数时,采用线性插值进行计算.10)@psl(x)单位正态线性损失函数,即返回max(0,z-x)的期望值,其中随机变量z服从标准正态分布.11)@psn(x)标准正态分布的分布函数在x点的取值.12)@ptd(n,x)自由度为n的t分布的分布函数在x点的取值.13)@qrand(seed)产生(0,1)区间的拟随机数.@qrand只允许在模型的数据部分使用,它将用拟随机数填满集属性.通常定义一个m×n的二维表,m表示运行实验的次数,n表示每次实验所需的随机数的个数.在行内,随机数是独立分布的;在行间,随机数是非均匀的.这些随机数是用“分层取样”的方法产生的.目前LINGO提供了两个金融函数.1)@fpa(I,n)返回如下情形的净现值:单位时段利率为I,连续n个时段支付,每个时段支付单位费用.若每个时段支付x单位的费用,则净现值可用x乘以@fpa(I,n)得到.@fpa的计算公式为.净现值就是在一定时期内为了获得一定收益,在该时期初所支付的实际费用.2)@fpl(I,n)返回如下情形的净现值:单位时段利率为I,第n个时段支付单位费用.@fpl(I,n)的计算公式为.这两个函数间的关系:.(7)输入和输出函数输入和输出函数可以把模型与外部数据(如文本文件、数据库和电子表格等)连接起来.1)@file函数该函数用于从外部数据文件中输入数据,它可以放在模型中任何地方.该函数的语法格式为@file(’filename’).这里filename是文件名,可以采用相对路径和绝对路径两种表示方式.记录结束标记(~)之间的数据文件部分称为记录.如果数据文件中没有记录结束标记,那么整个文件被看作单个记录.除了记录结束标记外,从模型外部调用的文本和数据同在模型里是一样的.下面介绍一下在数据文件中的记录结束标记连同模型中@file函数调用是如何工作的.当在模型中第一次调用@file函数时,LINGO打开数据文件,然后读取第一个记录;第二次调用@file 函数时,LINGO读取第二个记录等等.文件的最后一条记录可以没有记录结束标记,当遇到文件结束标记时,LINGO会读取最后一条记录,然后关闭文件.如果最后一条记录也有记录结束标记,那么直到LINGO 求解完成模型后关闭该文件.注意,如果有多个文件同时保持打开状态,可能就会导致一些问题,LINGO允许同时打开文件的上限数是16.在LINGO中不允许嵌套调用@file函数.2)@text函数该函数被用在数据部分,用来把求解结果输出至文本文件中.它可以输出集元素和集属性值.其语法为@text([’filename’])这里filename是文件名,可以采用相对路径和绝对路径两种表示方式.如果忽略filename,那么数据就被输出到标准输出设备(大多数情形都是屏幕).@text函数仅能出现在模型数据部分的一条语句的左边,右边是集名(用来输出该集的所有元素名)或集属性名(用来输出该集属性的值).用接口函数产生输出的数据定义称为输出操作.输出操作仅当求解器求解完模型后才执行,执行次序取决于其在模型中出现的先后.@OLE是从EXCEL中引入或输出数据的接口函数,它是基于传输的OLE技术.OLE传输直接在内存中传输数据,并不借助于中间文件.当使用@OLE时,LINGO先装载EXCEL,再通知EXCEL装载指定的电子数据表,最后从电子数据表中获得Ranges.为了使用@OLE函数,必须有EXCEL5及其以上版本.@OLE函数可在数据部分和初始部分引入数据.@OLE可以同时读集元素和集属性,集元素最好使用文本格式,集属性最好使用数值格式.原始集每个集元素需要一个单元(cell),而对于n元的派生集每个集元素需要n个单元,这里第一行的n个单元对应派生集的第一个集元素,第二行的n个单元对应派生集的第二个集元素,依此类推.4)@ranged(variable_or_row_name)为了保持最优基不变,变量的费用系数或约束行的右端项允许减少的量.5)@rangeu(variable_or_row_name)为了保持最优基不变,变量的费用系数或约束行的右端项允许增加的量.6)@status()返回LINGO求解模型后的结束状态:0 --- Global Optimum(全局最优);1 --- Infeasible(不可行);2 --- Unbounded(无界);3 --- Undetermined(不确定);4 --- Feasible(可行);5 --- Infeasible or Unbounded(通常需要关闭“预处理”选项后重新求解模型,以确定模型究竟是不可行还是无界)6 --- Local Optimum(局部最优);7 --- Locally Infeasible(局部不可行,尽管可行解可能存在,但是LINGO并没有找到一个);8 --- Cutoff(目标函数的截断值被达到);9 --- Numeric Error(求解器因在某约束中遇到无定义的算术运算而停止).通常,如果返回值不是0,4或6时,那么解将不可信,几乎不能用.该函数仅被用在模型的数据部分来输出数据.7)@dual(variable_or_row_name)返回变量的判别数(检验数)或约束行的对偶(影子)价格(dual prices).(8) 辅助函数1)@if(logical_condition,true_result,false_result)@if函数将评价一个逻辑表达式logical_condition是否为真,如果为真,返回true_ result,否则返回false_result.2)@warn(’text’,logical_condition)如果逻辑条件logical_condition为真,则产生一个内容为’text’的信息框.。
LINGO软件中集合的使用
![LINGO软件中集合的使用](https://img.taocdn.com/s3/m/5ff6714e0b4e767f5bcfce11.png)
LINGO软件中集合的使用作者:陶思俊邹海林来源:《硅谷》2008年第24期[摘要]集合是LINGO软件中重要的组成要素。
介绍集合的定义,结合实例说明集合的基本用法。
对集合做分类,详细讲解派生集合中的稠密集合、稀疏集合的一般形式和使用方法。
[关键词]LINGO软件派生集合稠密集合稀疏集合中图分类号:O1-0文献标识码:A 文章编号:1671-7597(2008)1220074-01一、引言LINGO软件是由美国LINDO系统公司研制开发的求解大型数学规划问题的软件包,可以用来求解线性规划、整数规划、二次规划、非线性规划问题等。
在LINGO软件中,集合及其属性的使用占据了相当重要的作用。
集合能够实现只输入少数几行文字就能建立起含有大规模变量的目标函数或成千上万条约束条件的数学模型。
二、集合的定义和基本用法(一)集合的定义集合(set)是一组相关对象的全体,一个集合可以由一系列产品、任务或时间组成。
集合的定义包含以下几个部分:集合的名字(setname);集合的成员(member_list);集合中成员的属性(attribute_list)。
其表达式为以“SETS:”开始,以“ENDSETS”结束,各字母不区分大小写,集合成员之间可以以逗号或空格分开,各属性之间也以逗号分开。
一般形式为(以下语法中凡是在方括号“[]”中的内容,表示是可选项,即该项可以有也可以没有):SETS:setname[/member_list/][:attribute_list];ENDSETS(二)LINGO软件的数据段在集合的使用过程中,通常需要事先给出集合中的某些属性赋值以便LINGO求解,这些赋值构成LINGO的数据段。
因此,LINGO中的数据段一般与集合搭配使用。
数据段的表达式以“DATA:”开始,以“ENDDATA”结束。
常数列表(value_list)中的数据之间以逗号或空格分开,一般形式为DATA:attribute=value_list(常数列表);ENDDATA(三)集合循环函数集合循环函数是指对集合上的元素(一般指下标)进行循环操作的函数。
lingo教程
![lingo教程](https://img.taocdn.com/s3/m/da330f8c6529647d2728523a.png)
LINGO 教 程
一个简单的LINGO程序
LINGO的基本用法的几点注意事项
•LINGO中不区分大小写字母;变量和行名可以超过8个字符,但 不能超过32个字符,且必须以字母开头。 •用LINGO解优化模型时已假定所有变量非负(除非用限定变量取 值范围的函数@free或@sub或@slb另行说明)。 •变量可以放在约束条件的右端(同时数字也可放在约束条件的左 端)。但为了提高LINGO求解时的效率,应尽可能采用线性表达 式定义目标和约束(如果可能的话)。 •语句是组成LINGO模型的基本单位,每个语句都以分号结尾,编 写程序时应注意模型的可读性。例如:一行只写一个语句,按照 语句之间的嵌套关系对语句安排适当的缩进,增强层次感。 •以感叹号开始的是说明语句(说明语句也需要以分号结束))。
•.MPS:示MPS(数学规划系统)格式的模型文件。
LINGO 教 程
在LINGO中使用LINDO模型
在LINGO中可以直接使用LINDO语法编写的优化模型(即优化程序)。 作为一个最简单的例子,在名为EXAM0201.LTX的模型文件中保存了一个 LINDO模型,我们现在看看如何用LINGO把它打开。
• 除具有LINDO的全部功能外,还可用于求解非线
性规划问题,包括非线性整数规划问题;
• 内置建模语言,允许以简练、直观的方式描述较 大规模的优化问题,所需的数据可以以一定格式 保存在独立的文件中。
7
LINGO 教 程
在LINGO中使用LINDO模型
LINGO的界面
• LINGO软件的主窗口(用 户界面),所有其他窗口 都在这个窗口之内。
16
LINGO 教 程
LINGO早期版本对LINDO的兼容问题
在LINGO 9.0以前的版本中不能直接用File|Open命令打开LINDO 模型,但由FILE | IMPORT LINDO FILE (F12)命令可以直接把 LINDO的模型文件转化成LINGO模型。运行后屏幕上会显示一 个标准的“打开文件”的对话框,打开EXAM0201.LTX,在 LINGO主窗口中又打开了命令窗口(Command Window)显示原始 文件,名为“exam0201”的模型窗口显示的是等价的LINGO模型。 当前光标位于命令窗口。
LINGO中的集
![LINGO中的集](https://img.taocdn.com/s3/m/e148e91d227916888486d793.png)
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];注意:用“[]”表示该部分内容可选。
LINGO用户手册
![LINGO用户手册](https://img.taocdn.com/s3/m/ddc18919581b6bd97f19ea39.png)
原语集是一组唯一的对象不能被进一步减小构成。在无线小工具的例子(第27页),设置仓库,这是由六个仓库,是一种原始的集合。同样,集合八家厂商组成是一种原始的集合。 派生集合是使用一个或多个其他组中定义。换句话说,一个来自集源自其他预先存在的集的成员。此外,使用无线小工具例如,集合六个仓库和八家厂商(LINKS)之间的联系组成一个派生集。它源于它的成员来自独特的对仓库的成员和供应商集。虽然链接组被从原始集仅衍生,但也可以从其他来源的集建立来自集为好。请参阅下面的部分,定义导集,以了解更多信息
例如,准备一个仓库发货型号为100的仓库,如果你有明确写入每个约束(例如,“仓库1必须装运不超过其目前的库存较多,仓库2必须装运不超过其目前的库存,仓库更是繁琐3必须运输不超过其本库存......“,等等)。 LINGO让您表达最简单的形式公式你阅读和理解(例如,“每个仓库出货必须不超过其目前的库存比较”)。
模型的集科
套在一个LINGO模型的可选部分称为集部分定义。在您的LINGO模型中使用集,你必须在模型组部分定义它们。套节与关键字集开始:(包括结肠),并与关键字ENDSETS结束。模型可以具有没有套部中,单个集合部,或多组部分。一个集部分可能会出现在任何地方的模型。唯一的限制是你必须定义一组及其属性之前,他们在模型的约束引用。
当明确列出成员,为每个成员输入唯一的名称,可选用逗号隔开。与集名称,成员名称也必须符合标准命名约定。在无线小部件模型,我们可以用一个明确的成员名单在集合部分定义
该setname可以是你选择指定设置一个名称。这应该是一个描述性的名称,很容易记住。集名称必须符合标准LINGO命名约定。换句话说,该名称必须以字母字符,这可能是随后长达31个字母数字字符或下划线(_)。 LINGO不大写和小写字符之间的名称区分开。
如何在lingo中使用集合1
![如何在lingo中使用集合1](https://img.taocdn.com/s3/m/06873159f78a6529647d53cb.png)
例题1. 在lingo 中输入下列线性规划模型,并求解∑∈⋅=A j i j i x j i d z ),(),(),( mins.t. 1),1(≥∑∈Vj j x ,,},10,,2,1{,0),(x ,),(,1,1),(V V A V Vi i i j i x j j i x Vi ⨯==∈=>=∑∈Λ为非负实数所有的数值如下表:dd=0 8 5 9 12 14 12 16 17 228 0 9 15 16 8 11 18 14 225 9 0 7 9 11 7 12 12 179 15 7 0 3 17 10 7 15 1512 16 9 3 0 8 10 6 15 1514 8 11 17 8 0 9 14 8 1612 11 7 10 10 9 0 8 6 1116 18 12 7 6 14 8 0 11 1117 14 12 15 15 8 6 11 0 1022 22 17 15 15 16 11 11 10 0;分析:这个模型输入的难点,在于变量的数量太多,足足有100个。
约束条件也比较多,有没有什么方便的输入方法?下面介绍lingo 中集合的建立新建lingo 文件输入下面内容model :sets :V/1..10/;!创建集合V;A(V,V):d,x;!创建集合A 是V 乘V.而d,x 是与A 同结构的,即d ,x 分别是10*10矩阵;endsetsmin =sum (A(i,j):d(i,j)*x(i,j));!创建目标函数;sum (V(j):x(1,j))>=1; !第一个约束条件;for (V(j)|j#gt#1:!i#gt#1为逻辑判断语句表示i>1是返回真值,但这里不能直接写i>1,因为">"是关系运算符不是逻辑运算符;sum (V(i):x(i,j))=1;); !利用循环函数表达:当i>1(即i 从2到10)时,{x(i,j):j=1..10}的和等于1; for (V(i):x(i,i)=0;);data:d=0 8 5 9 12 14 12 16 17 228 0 9 15 16 8 11 18 14 225 9 0 7 9 11 7 12 12 179 15 7 0 3 17 10 7 15 1512 16 9 3 0 8 10 6 15 1514 8 11 17 8 0 9 14 8 1612 11 7 10 10 9 0 8 6 1116 18 12 7 6 14 8 0 11 1117 14 12 15 15 8 6 11 0 1022 22 17 15 15 16 11 11 10 0; !输入d的数值; enddataend这个模型如果按照常规输入是什么样子呢?Lingo也提供这样的功能点lingo菜单选generate点display model得到下面的界面(部分截图)可见,这个规划如果普通输入的话,工作量是很大的,而且容易出错。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
例题1. 在lingo 中输入下列线性规划模型,并求解
∑∈⋅=A j i j i x j i d z ),(),(),( min
s.t. 1),1(≥∑∈V
j j x ,
,
},10,,2,1{,0),(x ,),(,
1,1),(V V A V V
i i i j i x j j i x V
i ⨯==∈=>=∑∈ 为非负实数
所有
的数值如下表:d d=0 8 5 9 12 14 12 16 17 22
8 0 9 15 16 8 11 18 14 22
5 9 0 7 9 11 7 12 12 17
9 15 7 0 3 17 10 7 15 15
12 16 9 3 0 8 10 6 15 15
14 8 11 17 8 0 9 14 8 16
12 11 7 10 10 9 0 8 6 11
16 18 12 7 6 14 8 0 11 11
17 14 12 15 15 8 6 11 0 10
22 22 17 15 15 16 11 11 10 0;
分析:这个模型输入的难点,在于变量的数量太多,足足有100个。
约束条件也比较多,有没有什么方便的输入方法?下面介绍lingo 中集合的建立
新建lingo 文件
输入下面内容
model :
sets :
V/1..10/;!创建集合V;
A(V,V):d,x;!创建集合A 是V 乘V.而d,x 是与A 同结构的,即d ,x 分别是10*10矩阵; endsets
min =@sum (A(i,j):d(i,j)*x(i,j));!创建目标函数;
@sum (V(j):x(1,j))>=1; !第一个约束条件;
@for (V(j)|j#gt#1:
!i#gt#1为逻辑判断语句表示i>1是返回真值,但这里不能直接写i>1,因为">"是关系运算符不是逻辑运算符;
@sum (V(i):x(i,j))=1;); !利用循环函数表达:当i>1(即i 从2到10)时,
{x(i,j):j=1..10}的和等于1;
@for(V(i):x(i,i)=0;);
data:
d=0 8 5 9 12 14 12 16 17 22
8 0 9 15 16 8 11 18 14 22
5 9 0 7 9 11 7 12 12 17
9 15 7 0 3 17 10 7 15 15
12 16 9 3 0 8 10 6 15 15
14 8 11 17 8 0 9 14 8 16
12 11 7 10 10 9 0 8 6 11
16 18 12 7 6 14 8 0 11 11
17 14 12 15 15 8 6 11 0 10
22 22 17 15 15 16 11 11 10 0; !输入d的数值; enddata
end
这个模型如果按照常规输入是什么样子呢?
Lingo也提供这样的功能
点lingo菜单
选generate点display model
得到下面的界面(部分截图)
运行结果(部分截图)
这个报告很不简洁,因为我们只需要显示x中非零部分,其他信息都不是我们关心的。
有什么办法可以使得报告更简洁?
点lingo菜单,选solution,做如下设置,可以使得报告更加简洁
设置完后结果是
练习1. 求下列线性规划模型
∑∈=A
j i j i x j i d z ),(),(*),( max
s.t. V i i b j i x V
j ∈<∑∈),(),(
),(j i x 为非负实数,(i,j)∈A
这里V={1,2,3,…,7},V V A ⨯=.
b=49 38 49 44 58 50 47;
d=1 2 3 1 2 1 2
2 3 4 1 2 2 3
4 1 1 2 1
5 1
1 2 2 2 3 1 1
2 1 2
3 5 3 1
2 2 2 1
3 2 4
1 1 1 1
2 1 1;
参考答案
例题2.一家公司需要决定四个季度的帆船生产量,下四个季度的需求分别为40,60,75,25条,需要按时满足,每个季度生产能力是40条,每条生产费用400美元。
如果加班,每条船的生产费为450美元。
每个季度末,每条船的库存费用为20美元。
假定生产提前期为0,库存为10条,如何安排生产使得总费用最小?
创建lingo文件,输入
运行结果
注1:如果创建一个集合包含的元素较多,例如月份,只需
Sets:
Month/1..12/
endsets
而不需要把1到12都罗列出来.
注2:上面的程序由三部分组成(集合段,目标和约束段,数据段),但一般lingo程序可以由五部分组成(新增“初值设定“和”计算”两部分)
(1)集合段sets开始,endsets结束
(2)目标函数(min=或max=)和约束段(一些等式或不等式,或逻辑判断语句)(3)数据段data开始,enddata结束
(4)初始段,用于初值设定,init开始,endinit结束,对集合对属性(数组)定义初值(因为求解算法一般是迭代算法,用户如果能给出一个比较好的迭代初值,对提高算法的效果是有益处的)
(5)计算段,以calc开始,endcalc结束,作用在于对一些原始数据进行计算处理。
为什么要设计这个段呢?因为在实际问题中,可能出现输入的原始数据不能再模型中直
接使用,需要进行“预处理”。
练习2.用lingo求解下列问题
练习3.利用lingo的集合编写程序求解下列问题
练习4,利用lingo求解下列问题。