6优化Lingo
LINGO简介
例3 求解0—1整数线性规划问题 MAX f=-3x1+2x2-5x3 S.T x1+ x2- x3≤2 x1+4x2+x3≤4 x1+ x2 ≤3 4x2+x3≤6 xi (i=1,2,3)为0或1 程序如下:
MAX=-3*x1+2*x2-5*x3; x1+ x2-x3<=2; x1+4*x2+x3<=4; x1+x2<=3; 4*x2+x3<=6; @BIN(x1); @BIN(x2); @BIN(x3);
1、集合定义部分
这部分以“SETS:”开始,以 “ENDSETS”结束,作用在于定义必要的集合 变量(SET)及其元素(member,类似于数组 的下标)和属性(attribute,类似于数组)。 格式是: 集合名/1..n/:属性列表; 各属性之间用“,”分隔。属性的具体数值如 果是常量,则可在数据输入部分输入;如果是 未知量,则可在初始部分输入它的初值。
例 7、求解线性规划(LP)问题 min S.T
f (400xi 450yi 20zi )
i 1 4
zi zi 1 xi yi ri
xi 40 (i 1,,4)
(i 1,,4)
ri 40,60,75,25 (i 1,,4)
z0 10
三、LINGO模型的基本要素
当你遇到的是一个较大规模的优化问题时, 怎样才能做到输入几个语句就可以建立起含有 很多变量的目标函数和成千上万条约束的一个 复杂庞大的模型呢? 掌握LINGO提供的最优化模型语言是至关 重要的。 一个最复杂的LINGO模型由以下5个要素组 成(当然,如果实际问题不需要,则其中的某 些要素可以省略)。
优化软件LINGO使用入门
-实例1-
程序编写
model: title 生产计划问题; [maxf]max=2*x1+3*x2; [TIME]x1+2*x2<8; [A]4*x1<16; [B]4*x2<12; END
华中农业大学建模基地
-实例1-
运行结果 Model Title: 生产计划问题 Variable Value X1 4.000000 X2 2.000000 Row Slack or Surplus MAXF 14.00000 TIME 0.000000 A 0.000000 B 4.000000
更新时间间隔lingo实验华中农业大学建模基地求解报告窗口lingo实验华中农业大学建模基地2412302412问如何调运使运费最低如下公里单位距离两个粮库到三个粮站的大米分别为三个粮站至少需要两个粮库现存大米分别调运大米向三个粮站有两个粮库运输问题lingo实验华中农业大学建模基地1212243024lingo实验华中农业大学建模基地结合存量限制和需量限制得数学模型
华中农业大学建模基地
-界面介绍-
• LINGO软件的主窗口(用 软件的主窗口( 软件的主窗口 户界面), ),所有其他窗口 户界面),所有其他窗口 都在这个窗口之内。 都在这个窗口之内。
• 当前光标的 位置 • 模型窗口(Model 模型窗口( Window),用于输入 ),用于输入 ), LINGO优化模型(即 优化模型( 优化模型 LINGO程序)。 程序)。 程序
华中农业大学建模基地
-实例1-
生产计划问题 某工厂计划安排生产Ⅰ,Ⅱ两种产品, 已知每种单位产品的利润,生产单位产品所需设备台时 及A,B两种原材料的消耗,现有原材料和设备台时的定额 如表所示,问:怎么安排生产使得工厂获利最大?
2.Lingo优化软件的使用方法
然后再在求解的基础上,在菜单上选择lingo/range(ctr+R), 计算机会自动对每个变量和资源用有量进行范围分析:保 持最优基不变的情况下,变量或约束行的右端允许增加和 减少的量。注意:对整数规划和0-1规划灵敏度分析意义不 大。
Ranges in which the basis is unchanged:
数学模型如下:
max z [ t (i) r (i)]x (i) s(i) y(i)
i 1 i 1
3
3
s.t.
3 c ijx i b( j) j 1,2,3 i 1 x (i) My(i) i 1,2,3 ___ x (i) 0, x Z y(i) 0,1 i 1,2,3
产品1
资源A 资源B 资源C 2 2 1
产品2 产品3
4 3 2 8 4 3
资源量
500 300
单件可变费 固定费用 单件售价
4 100 8
5 150 10
6 200 12
分析问题和设置变量 x(i)表示第i种产品的产量;i=1,2,3 c(i,j)表示第i种产品对第j种资源的消耗量;i=1,23,j=1,23 b(j)表示第j种资源的拥有量; r(i)表示第i种产品的单件可变费用; s(i)表示第i种产品的固定费用; t(i)表示第i种产品的单位售价。 总收益=销售收入-固定费用-可变费用,由于不知道第i种产品 生产与否,所以,必须给定第i种产品的选择变量 y(i)表示第I种产品的生产选择,y(i)=0,表示不生产第i种产品, y(i)=1表示生产第i种产品。
#LT#(小于) #LE#(小于等于) 逻辑运算的结果只有“真”(TRUE)和“假”(FALES), Llingo用1表示True,其它的都是False。
6优化Lingo
@for(ren(i):@sum(yundong(j):x(i,j))=1);
@for(yundong(j):@sum(ren(i):x(i,j))=1);
end
例2 一家大建筑公司正在三个地点开掘。同 时又在其他四个地点建筑,这里需要土方的 填充。在1、2、3处挖掘产生的土方分别为每 天150,400,325立方码。建筑地点A、B、C、 D处需要的填充土方分别为175,125,225, 450立方码。也可以从地点4用每立方码5美元 的价格获得额外的填充土方。填充土方运输 的费用约为一货车容量每英里20美元。一辆 货车可以搬运10立方码的土方(每立方码土方 每英里运输费2美元)。表3-3给出了各地点间 距离的英里数。求使公司花费最少的运输计 划。
尽量少使用绝对值符号函数多个变量求最大最小值四舍五入取整函数等3尽量使用线性模型减少非线性约束和非线性变量的个数改为x5y4合理设定变量上下界尽可能给出变量初始值5模型中使用的参数数量级要适当如小于101lindo
数学模型与实验(六)
关于LINDO/LINGO软件
LINDO 公司软件产品简要介绍
美国芝加哥(Chicago)大学的Linus Schrage教授于1980 年前后开发, 后来成立 LINDO系统公司(LINDO Systems Inc.), 网址:
LINDO API: LINDO Application Programming Interface (V2.0)
演示(试用)版、学生版、高级版、超级版、工业版、 扩展版… (求解问题规模和选件不同)
建模时需要注意的几个基本问题
1、尽量使用实数优化,减少整数约束和整数变量
2、尽量使用光滑优化,减少非光滑约束的个数
练习7* 2005年D题: DVD在线租赁
2024年度LINGO软件
结果分析
对求解结果进行分析,验证模型的有 效性和可行性。
模型调整
根据结果分析,对模型进行调整和优 化,提高模型的实用性和准确性。
24
05 LINGO软件在各 个领域的应用案 例
2024/3/23
25
生产计划与调度问题
2024/3/23
生产线平衡
LINGO可以用于解决生产线平衡问题,通过优化生产线上 各个工位的任务分配,提高生产效率和资源利用率。
一些学术机构和研究团队也会分享他们的 LINGO使用经验和案例,为学习者提供更多 的学习资源和合作机会。
37
THANKS
感谢观看
2024/3/23
38
LINGO可用于金融市场预测和决策分析,通过建立预测模型和分析工具
,揭示金融市场的运行规律和趋势,为投资决策提供支持。
29
06 LINGO软件与其 他工具的集成与 应用
2024/3/23
30
与Excel的集成与应用
2024/3/23
数据交换
LINGO可以直接从Excel中导入数据,也可以将结果导出到Excel ,实现数据的无缝对接。
2024/3/23
28
金融工程与投资组合优化问题
01
投资组合优化
LINGO可以用于解决投资组合优化问题,通过优化投资组合中各个资产
的配置比例,实现风险和收益的平衡。
02
期权定价与风险管理
利用LINGO建立期权定价模型,可以准确计算期权的价值,为风险管理
提供决策依据。
2024/3/23
03
金融市场预测与决策分析
2024/3/23
整数规划算法
分支定界法、割平面法等适用于整数规划问 题的求解。
lingo软件 优化问题
2.2 LP模型在LINGO中的一个典型输入方式
集合定义部分从 (“SETS:”到 “ENDSETS” ): 定义集合及其属性
MIN
I 1,2,3,4
{400 RP ( I ) 450OP( I ) 20 INV ( I )}
RP( I ) 40, I 1, 2, 3, 4 INV ( I ) INV ( I 1) RP ( I ) OP ( I ) DEM ( I ), I 1, 2, 3, 4 INV (0) 10 INV ( I ) 0 OP ( I ) 0 RP ( I ) 0 I 1 4
Step5 回答问题
Global optimal solution found. Objective value: 78450.00 Total solver iterations: 2 Variable Value Reduced Cost RP1 40.00000 0.000000 RP2 40.00000 0.000000 RP3 40.00000 0.000000 RP4 25.00000 0.000000 OP1 0.000000 20.00000 OP2 10.00000 0.000000 OP3 35.00000 0.000000 OP4 0.000000 50.00000 INV1 10.00000 0.000000 INV2 0.000000 20.00000 INV3 0.000000 70.00000 INV4 0.000000 420.0000 Row Slack or Surplus Dual Price 1 78450.00 -1.000000 2 0.000000 30.00000 3 0.000000 50.00000 4 0.000000 50.00000 5 15.00000 0.000000 6 0.000000 430.0000 7 0.000000 450.0000 8 0.000000 450.0000 9 0.000000 400.0000
lingo最值和优化
ymax = -6.6080e-004
0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 -0.05
0
1
2
3
4
5
6
7
8
>> f = @(x)exp(x(2))*(4*x(1)^2+2*x(2)^2+4*x(1)*x(2)+2*x(2)+1); x0 = [-1, 1]; x=fminunc(f,x0); y=f(x) Warning: Gradient must be provided for trust-region method; using line-search method instead. > In fminunc at 243 Optimization terminated: relative infinity-norm of gradient less than options.TolFun. y = 1.4859e-013
fval = -11.0250
>> f=@(x)x.^3-2*x-5; z = fminbnd(f, 0, 2) x=0:0.1:2; f=@(x)x.^3-2*x-5; plot(x,f(x),'b',z,f(z),'rp') axis([0,2,-50,50]) legend('f(x)','极小点 极小点') 极小点 z =
>> f=@(x)(2*x(1)^2+4*x(2)^2-4*x(1)*x(2)-6*x(1)-3*x(2)); A=[1,1;4,1]; b=[3;9]; x0=[1;1] [x,fval]=fmincon(f,x0,A,b) x0 = 1 1 Warning: Large-scale (trust region) method does not currently solve this type of problem, switching to medium-scale (line search). > In fmincon at 274 Optimization terminated: first-order optimality measure less than options.TolFun and maximum constraint violation is less than options.TolCon. Active inequalities (to within options.TolCon = 1e-006): lower upper ineqlin ineqnonlin 1 x = 1.9500 1.0500
LINGO软件简介
LINGO 软件简介LINGO 软件是一个处理优化问题的专门软件,它尤其擅长求解线性规划、非线性规划、整数规划等问题.一个简单示例有如下一个混合非线性规划问题:⎪⎩⎪⎨⎧≥≤≤+++---+为整数2132121321322212121,;0,,210022..15023.027798max x x x x x x x x x x t s x x x x x x x .LINGO 程序模型:max =98x1+277x2-x1^2-0.3x1x2-2x2^2+150x3; x1+2x2+2x3<=100; x1<=2x2;gin x1;gin x2; Lingo 默认变量非负注意:binx 表示x 是0-1变量;ginx 表示x 是整数变量;bndL,x,U表示限制LxU ;freex 表示取消对x 的符号限制,即可正、可负.结果:Global optimal solution found.Objective value: 9561.200 Extended solver steps: 0 Total solver iterations: 45 Variable Value Reduced CostX1 6.000000 -76.70000X2 31.00000 -151.2000X3 16.00000 -150.0000Row Slack or Surplus Dual Price1 9561.200 1.0000002 0.000000 0.0000003 56.00000 0.000000———————— 非常简单在LINGO 中使用集合为了方便地表示大规模的规划问题,减少模型、数据表示的复杂程度,LINGO 引进了“集合”的用法,实现了变量、系数的数组化下标表示.例如:对⎪⎪⎩⎪⎪⎨⎧==-++-==≤++∑=.,,;10)0(;4,3,2,1),()())()1()(;4,3,2,1,20)(..)}(20)(450)(400{min4,3,2,1均非负INV OP RP INV I I DEM I OP I RP I INV I INV I I RP t s I INV I OP I RP I求解程序:model :sets :mark/1,2,3,4/:dem,rp,op,inv;也可以vmark/1..4/:dem,rp,op,inv;endsetsmin=sum mark:400rp+450op+20inv;也可以markI:400rpI+450opI+20invI;for markI: rpI<40;for markI|Igt1: invI=invI-1+rpI+opI-demI;inv1=10+rp1+op1-dem1;data:dem=40,60,75,35;enddataend上面程序在model…end之间有1集合定义、2数据输入和3其他三部分内容.集合定义部分从sets:到endsets:定义了一个指标集合mark可以理解为数组下标及其范围和其4个属性dem、rp、op、inv用此向量的数组变量.数据输入部分从data:到enddata依次给出常量dem的值.其他部分:给出优化目标及约束.一般而言,LINGO中建立优化模型的程序可以由五部分组成,或称为五段section:1集合段SETS:这部分以“SETS:”开始,以“ENDSETS”结束,作用在于定义必要的集合变量SET及其元素member,含义类似于数组的下标和属性attribute,含义类似于数组.2目标与约束段:这部分实际上定义了目标函数、约束条件等,但这部分没有段的开始和结束标记;该段一般常用到LINGO内部函数,尤其是和集合相关的求和函数SUM和循环函数FOR等.3数据段DATA:这部分以“DATA:”开始,以“ENDDATA”结束,作用在于对集合的属性数组输入必要的常数数据.格式为:attribute属性=value_list常数列表;常数列表中的数据之间可以用逗号、空格或回车符分隔.如果想要在运行时才对参数赋值,可以在数据段使用输入语句,其格式为“变量名=;”,但仅限对单个变量赋值,而不能用于属性变量数组的单个元素.4初始段INIT:这部分以“INIT:”开始,以“ENDINIT”结束,作用在于对集合的属性数组定义初值因为求解算法一般是迭代算法,提供一个较好的初值,能提高计算效果.定义初值的语句格式为:attribute属性=value_list常数列表;这与数据段中的用法类似.5计算段CALC:这部分以“CALC:”开始,以“ENDCALC”结束,作用在于对一些原始数据进行预处理加工,使其成为模型直接需要的数据.该段中通常是计算赋值语句.基本集合与派生集合为了处理二维数组变量等有多个下标的问题,LINGO引入了“派生集”的概念.我们把直接列出元素的指标集合叫“基本集合”,而基于其他集合派生出来的二维或多维指标集合称为“派生集”.派生集的定义格式为:派生集名原始集合1,原始集合2,…,原始集合n:属性变量列表;实际上就是笛卡儿积的意思,即:派生集={i1,i2, (i)n| i1集合1, i2集合2,…, in集合n}.1一个应用例子布局问题:某些建筑工地的位置用平面坐标a,b表示及水泥日用量d已知.现有A、B两临时料场位于P5,1、Q2,7,日储量20.问A、B两料场分别向各工地运输多少吨水泥,使总吨公里数最小若重新安排两料场的位置,应怎样安排才能使总吨公里数最小这样安排可节省多少吨公里设工地位置ai ,bi,水泥日用量为dii=1,2,…,6;料场位置xi,yi,日储量ej,j=1,2;从料场j向工地i运送量为cij.该问题的数学模型为:LINGO求解程序为:MODEL:sets:Imark/1..6/:a,b,d;Jmark/1,2/:x,y,e;IJmarkImark,Jmark:c;endsetsdata:Location for demand需求点位置;a=1.25,8.75,0.5,5.75,3,7.25;b=1.25,0.75,4.75,5,6.5,7.75;Quantities of the demand and supply供需量;d=3,5,4,7,6,11;e=20,20;enddatainit:Initial location for the supply初始点;x,y=5,1,2,7;endinitObjective function目标;OBJ min=sum IJmarki,j: ci,jxj-ai^2+yj-bi^2^1/2; demand contraints需求约束;for Imarki:DEMAND_CON SUM Jmarkj:ci,j=di;; supply constrains供给约束;for Jmarkj:SUPPLY_CON SUM Imarki:ci,j<=ej;;for Jmark: free x;free y;;2一个动态规划的例子:最短路问题从S城市到T城市之间找一条最短路径,道路情况如下:数学模型为:LINGO求解程序:model:sets:cities/s,a1,a2,a3,b1,b2,c1,c2,t/:L; 属性Li表示城市S到城市i的最优行驶路线的里程;roadscities,cities/ 派生集合roads表示的是网络中的道路;s,a1 s,a2 s,a3 由于并非所有城市间都有道路直接连接,所以将路具体列出;a1,b1 a1,b2 a2,b1 a2,b2 a3,b1 a3,b2b1,c1 b1,c2 b2,c1 b2,c2 属性Di,j是城市i到城市j的直接距离已知;c1,t c2,t/:D;endsetsD= 6 3 36 5 8 67 46 7 8 95 6;L=0,,,,,,,,; 因为Ls=0;enddatafor citiesi|igt index s: 这行中"indexs"可以直接写成"1";Li=min roadsj,i:Lj+Dj,i;; 这就是最短路关系式;endVariable ValueL S0.000000L A16.000000L A23.000000L A33.000000L B110.00000L B27.000000L C115.00000L C216.00000L T20.00000最短路径为: S-〉A3-〉B2-〉C1-〉T3指派问题设有6个人做6件事.其中cij表示第i人做第j事的收益;设第i人做第j事时xij =1,否则xij=0.该问题的规划模型:说明:其中“-”表示某人无法做该事.可令其为-表示绝对不行或0领薪不用干活LINGO求解程序:MODEL:sets:Imark/1..6/:i;Jmark/1..6/:j;IJmarkImark,Jmark:c,x;endsetsdata:第i人做第j事的收益;c=20,15,16,5,4,717,15,33,12,8,69,12,18,16,30,1312,8,11,27,19,14-99,7,10,21,10,32-99,-99,-99,6,11,13;enddataOBJ max=sum IJmarki,j: cx;每人做一项工作;for Imarki: SUM Jmarkj:xi,j=1;;每事一人做;for Jmarkj: SUM Imarki:xi,j=1;;for IJmark: bin x;本约束可以不要,因为有解时必为0或1; END4生产与销售计划问题某公司用两种原油A 和B 混合加工成两种汽油甲和乙.甲、乙两种汽油含原油A 的最低比例分别为50%和60%,每吨售价分别是4800元和5600元.该公司现有原油A 和B 的库存量分别为500吨和1000吨,还可以从市场上买到不超过1500吨的原油A.原油A 的市场价为:购买量不超500吨时单价为10000元/吨;购买量超过500吨但不超1000吨时,超过500吨部分单价为8000元/吨;购买量超过1000吨部分的单价是6000元/吨.该公司应如何安排原油的采购和加工以获得最大利润数学模型: 设原油A 用于生产甲、乙两种汽油的数量分别是x11和x12,原油B 用于生产甲、乙两种汽油的数量分别是x21和x22;购买原油A 的数量是x 吨,采购支出为cx 千元/吨.为了处理分段函数cx,将原油采购量x 分解为对应价格10千元/吨的采购量x1、对应对应价格8千元/吨的采购量x2和对应价格6千元/吨的采购量x3,它们应满足:0)500(21=-x x 表示要么x1=500要么x2=0,即x1的量不达到500时x2=00)500(32=-x x 表示要么x2=500要么x3=0,即x2的量不达到500时x3=0此时采购支出3216810)(x x x x c ++=模型改变为:LINGO 求解程序:model :init:x1=500;x2=500;x3=0;x12=1500;x22=1000;x11=0;x21=0;endinitmax=4.8x11+4.8x21+5.6x12+5.6x22-10x1-8x2-6x3; x11+x12<=x+500;x21+x22<=1000;0.5x11-0.5x21>=0;0.4x12-0.6x22>=0;x=x1+x2+x3;x1-500x2=0;x2-500x3=0;bnd0,x1,500;bnd0,x2,500;bnd0,x3,500;。
《LINGO简介》课件
某些复杂的数学表达可能无法直接在 LINGO中表示。用户可以通过混合编 程或使用其他建模语言(如GAMS) 来解决这一问题。
对于特定行业或领域的定制化需求, LINGO可能无法直接提供相应的功能 。在这种情况下,用户可以通过扩展 LINGO的API或与其他软件的集成来 实现定制化需求。
感谢您的观看
目标函数的设置
目标函数定义
在LINGO中,需要定义一个目标函数来描述决策变量 的优化目标。
目标函数类型
目标函数可以是最大化或最小化形式,根据实际问题 的需求进行选择。
目标函数编辑器
LINGO提供了一个目标函数编辑器,用户可以在其中 方便地定义和编辑目标函数。
求解操作
求解器选择
在LINGO中,可以选择不同的求解器 来求解模型,根据模型的规模和复杂
LINGO软件广泛应用于生产计划、资源分配 、工艺流程优化等方面。
物流运输
LINGO软件用于运输路线规划、车辆调度、 仓储优化等问题求解。
金融投资
LINGO软件用于投资组合优化、风险管理、 信贷决策等问题求解。
科研领域
LINGO软件在数学建模、统计分析、机器学 习等领域有广泛应用。
02
LINGO软件的基本操作
物流配送问题
总结词
物流配送问题是一个复杂的优化问题,LINGO软件能够通过建立有效的数学模型,优化配送路线和成本。
详细描述
物流配送问题涉及到如何合理规划配送路线、分配运输资源,以最小化运输成本并确保及时送达。LINGO软件通 过构建配送问题的数学模型,帮助企业找到最优的配送方案,降低运输成本、提高运输效率。
LINGO软件与其他软件的比较与选择
MATLAB
MATLAB在科学计算和数据分析领域具有广泛的应用,但 相比之下,LINGO在求解优化问题方面更加专业和高效。
优化模型与Lingo软件的使用
x 3 x 5 2 x 7 15
整数约束: xi 为整数 最优解:x2=12, x5=15, 其余为0; 最优值:27
模 式 1 2 3 4 5 6 7 需 求
4米 根数 4 3 2 1 1 0 0 50
6米 根数 0 1 0 2 1 3 0 20
下的最大值或最小值,其中
x
设计变量(决策变量) 目标函数 可行域
f (x)
x
二 优化模型的分类
1.根据是否存在约束条件 有约束问题和无约束问题。 2.根据设计变量的性质 静态问题和动态问题。 3.根据目标函数和约束条件表达式的性质 线性规划,非线性规划,二次规 划,多目标规划等。
模型的基本类型
可用下面的语法定义一个派生集:
setname(parent_set_list)[/member_list/][:attribute_list];
setname是集的名字。parent_set_list是已定义的集的列表,多个时必须用逗号隔开。 如果没有指定成员列表,那么LINGO会自动创建父集成员的所有组合作为派生集的 成员。派生集的父集既可以是原始集,也可以是其它的派生集。 例: sets: product/A B/; machine/M N/; week/1..2/; allowed(product,machine,week):x; endsets
g i ( x ) 0 ( g i ( x ) 0 ), i 1, 2 ,..., p .
例如:钢管下料问题
客户需求 4米50根
6米20根 原料钢管:每根19米 8米15根 节省的标准是什么? 5米10根
问题1. 如何下料最节省 ? 问题2. 客户增加需求:
LINGO的使用方法说明大全
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;endsetsdata: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)实现同样的功能.例3.6sets: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%范围内,对不同的值求解模型,观察模型的结果对参数依赖的程度,那么把这种情况称为实时数据处理.处理方法是在该语句的数值后面输入一个问号(?).data:interest_rate,inflation_rate = .09 ?;enddata在每一次求解模型时,LINGO都会提示为参数inflation_rate输入一个值.在WINDOWS操作系统下,将会看到一个如下面的对话框:直接输入一个值再点击OK按钮,LINGO就会把输入的值指定赋给inflation_rate,然后继续求解模型.除了参数之外,也可以实时输入集的属性值,但不允许实时输入集元素名.(4) 指定属性为一个值可以在数据定义的右边输入一个值来把所有的元素的该属性指定为一个值.如下面的例子.例3.9sets: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,,,;属性capacity的第2个和第3个值分别为24和40,其余的未知.3.3.2初始部分初始部分是LINGO提供的另一个可选内容.在初始部分中,与数据部分中的数据定义相同,可以输入初始定义(initialization statement).在对实际问题的建模时,初始部分并不起到描述模型的作用,初始部分输入的值仅被LINGO求解器当作初始值来使用,并且仅仅对非线性模型有用.这与数据部分指定变量的值不同,LINGO求解器可以自由改变初始部分初始化变量的数值.一个初始部分以关键字“init:”开始,以关键字“endinit”结束.初始部分的初始定义规则和数据部分的数据定义规则相同.也就是说,可以在定义的左边同时初始化多个集属性,即可以把集属性初始化为一个数值,也可以用问号定义为实时数据,还可以用逗号指定为未知数值.例3.11init:X,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 运算符的优先级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) 集合操作函数集合操作函数是对集合进行操作的函数,主要有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表示每次实验所需的随机数的个数.在行内,随机数是独立分布的;在行间,随机数是非均匀的.这些随机数是用“分层取样”的方法产生的.(6) 金融函数目前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(’’).这里是文件名,可以采用相对路径和绝对路径两种表示方式.记录结束标记(~)之间的数据文件部分称为记录.如果数据文件中没有记录结束标记,那么整个文件被看作单个记录.除了记录结束标记外,从模型外部调用的文本和数据同在模型里是一样的.下面介绍一下在数据文件中的记录结束标记连同模型中@file函数调用是如何工作的.当在模型中第一次调用@file函数时,LINGO打开数据文件,然后读取第一个记录;第二次调用@file 函数时,LINGO读取第二个记录等等.文件的最后一条记录可以没有记录结束标记,当遇到文件结束标记时,LINGO会读取最后一条记录,然后关闭文件.如果最后一条记录也有记录结束标记,那么直到LINGO 求解完成模型后关闭该文件.注意,如果有多个文件同时保持打开状态,可能就会导致一些问题,LINGO允许同时打开文件的上限数是16.在LINGO中不允许嵌套调用@file函数.2)@text函数该函数被用在数据部分,用来把求解结果输出至文本文件中.它可以输出集元素和集属性值.其语法为@text([’’])这里是文件名,可以采用相对路径和绝对路径两种表示方式.如果忽略,那么数据就被输出到标准输出设备(大多数情形都是屏幕).@text函数仅能出现在模型数据部分的一条语句的左边,右边是集名(用来输出该集的所有元素名)或集属性名(用来输出该集属性的值).用接口函数产生输出的数据定义称为输出操作.输出操作仅当求解器求解完模型后才执行,执行次序取决于其在模型中出现的先后.3)@ole函数@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’的信息框.3)@user(user_determined_arguments)该函数允许用户自己编写函数,可以用c语言等编写,返回值为用户函数计算的结果.3.5 LINGO程序出错信息在LINGO模型求解时,系统会对程序进行编译、求解或是执行于程序相关的命令,这都有可能出现一些语法或运行的错误.当出现时,系统会弹出一个出错报告框,显示错误代码,并且大致指出错误的所在位置.这些错误信息报告对于用户发现及改正程序中的错误有很大帮助.如附表3-4就出错提示信息,进行说明(没有说明的错误编号目前还没有使用).。
LINGO使用说明比较简单
LINGO使⽤说明⽐较简单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
集循环函数
集循环函数遍历整个集进行操作。其语法为: @function(setname[(set_index_list)[|conditional_qualifier]]:e xpression_list);
Lingo的四个集循环函数: @for(x()、@min():求最大值、最小值。
Lingo程序结构 Lingo程序结构
集合段,定义集合
数据段,输入数据
数据初始化段,赋 予数据初值
!目标函数; 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)); end
LINGO优化软件 LINGO优化软件
美国LINDO公司的主要产品: 公司的主要产品: 美国 公司的主要产品 LINDO: Linear INteractive and Discrete Optimizer (V6.1) LINGO: Linear INteractive General Optimizer (V10.0) LINDO API: LINDO Application Programming Interface (V5.0)
sets: product/A B/; machine/M N/; week/1..2/; allowed(product,machine,week):x; endsets 派生集成员:
定义了三个原 始集,allowed 为三个原始集 构成的派生集。 X为派生集的属 性。
(A,M,1) (A,M,2) (A,N,1) (A,N,2) 集成员按这个顺序排列。 (B,M,1) (B,M,2) (B,N,1) (B,N,2)
Lingo软件在求解数学优化问题的使用技巧
3*x3+2*y3+z3<=375;
END
得到的解如下:
X1=200,Y1=0,Z1=0;
X2=58.33333, Y2=312.5,Z2=0;
X3=0,Y3=187.5,Z3=0;
最大总净收益为253333.3元。
3.公司在各地有4项业务,选定了4位业务员去处理。由于业务能力、经验和其它情况不同,4业务员去处理4项业务的费用(单位:元)各不相同,见下表:
常见的集合函数如下:
@FOR(set_name:constraint_expressions)对集合(set_name)的每个元素独立地生成约束,约束由约束表达式(constraint_expressions)描述。
@MAX(set_name:expression)返回集合上的表达式(expression)的最大值。
=
返回如下情形下的净现值:单位时段利率为 ,第 个时段支付单位费用,即:
=
(5)概率函数
@PSN(X)标准正态分布的分布函数。
@PSL(X)单位正态线性损失函数(即返回 的期望值,其中Z为标准正态随机变量)
@PPS(A,X)均值为A的Possion分布的分布函数(当X不是整数时,采用线性插值进行计算)。
(4)初始化部分(INIT):这部分以“INIT:”开始,以“END INIT”结束。作用在于对集合的属性(数组)定义初值。格式为:attribute=value_list。由于非线性规划求解时,通常得到的是局部最优解,而局部最优解受输入的初值影响。通常可改变初值来得到不同的解,从而发现更好的解。
编写LINGO程序要注意的几点:
@PPL(X)Possion分布的线性损失函数(即返回 的期望值,其中Z为Possion分布随机变量)
LINGO课件
LINGO模型 LINGO模型 — 例:选址问题
某公司有6个建筑工地,位置坐标为 单位: 某公司有 个建筑工地,位置坐标为(ai, bi) (单位:公里 个建筑工地 单位 公里), 单位: 水泥日用量d 单位 水泥日用量 i (单位:吨)
i a b d 1 1 .2 5 1 .2 5 3 2 8 .7 5 0 .7 5 5 3 0 .5 4 .7 5 4 4 5 .7 5 5 7 5 3 6 .5 6 6 7 .2 5 7 .7 5 11
集是一群相联系的对象,这些对象也称为集的成员。 集是一群相联系的对象,这些对象也称为集的成员。一个集可能 是一系列产品、卡车或雇员。 是一系列产品、卡车或雇员。每个集成员可能有一个或多个与之 有关联的特征,我们把这些特征称为属性。属性值可以预先给定, 有关联的特征,我们把这些特征称为属性。属性值可以预先给定, 也可以是未知的,有待于LINGO求解。例如,产品集中的每个产 求解。 也可以是未知的,有待于 求解 例如, 品可以有一个价格属性; 品可以有一个价格属性;雇员集中的每位雇员可以有一个薪水属 也可以有一个生日属性等等。 性,也可以有一个生日属性等等。
•对实际问题建模的时候,总会遇到一群或多群相联系 对实际问题建模的时候, 对实际问题建模的时候 的对象,比如工厂、消费者群体、交通工具和雇工等等。 的对象,比如工厂、消费者群体、交通工具和雇工等等。 LINGO允许把这些相联系的对象聚合成集(sets)。 允许把这些相联系的对象聚合成集( 允许把这些相联系的对象聚合成集 )。
LINGO优化软件 优化软件 学习课件
LINGO软件简介 软件简介
LINGO模型的优点 模型的优点
•包含了 包含了LINDO的全部功能 包含了 的全部功能 •提供了灵活的编程语言(矩阵生成器) 提供了灵活的编程语言( 提供了灵活的编程语言 矩阵生成器)
什么是Lingo
LINGO目录简介步骤LINGO 综述编辑本段简介LINGO是Linear Interactive and General Optimizer的缩写,即“交互式的线性和通用优化求解器”,由美国LINDO系统公司(Lindo System Inc.)推出的,可以用于求解非线性规划,也可以用于一些线性和非线性方程组的求解等,功能十分强大,是求解优化模型的最佳选择。
其特色在于内置建模语言,提供十几个内部函数,可以允许决策变量是整数(即整数规划,包括0-1 整数规划),方便灵活,而且执行速度非常快。
能方便与EXCEL,数据库等其他软件交换数据。
编辑本段步骤一般地,使用LINGO 求解运筹学问题可以分为以下两个步骤来完成:1)根据实际问题,建立数学模型,即使用数学建模的方法建立优化模型;2)根据优化模型,利用LINGO 来求解模型。
主要是根据LINGO 软件,把数学模型转译成计算机语言,借助于计算机来求解。
例题:在线性规划中的应用max Z =5 X1+3 X2+6X3,s.t. X1 +2 X2 + X3 ≤182 X1 + X2 +3 X3 =16X1 + X2 + X3 =10X1 ,X2 ≥0 , X3 为自由变量应用LINGO 来求解该模型,只需要在lingo窗口中输入以下信息即可:max=5•x1 +3•x2 +6•x3 ;x1 +2•x2 + x3 <=18 ;2•x1 + x2+3•x3 =16 ;x1 + x2 + x3 =10 ;@free( x3) ;然后按运行按钮,得到模型最优解,具体如下:Objective value: 46.00000Variable Value Reduced Costx1 14.00000 0.000000x2 0.000000 1.000000x3 -4 .000000 0.000000由此可知,当x1 =14 , x2 =0 , x3 =-4 时,模型得到最优值,且最优值为46。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
每头奶牛需要1.5英亩土地,冬季需要付出100小 时劳动时间,夏季50小时,每年净收益为450美 元;相应地,每只母鸡不占用土地,冬季0.6小 时,夏季0.3小时,年净收益为3.5美元。养鸡房 最多容纳3000只母鸡,栅拦最多能容纳32头奶牛。
种植一英亩的大豆、玉米、燕麦分别需要冬季劳 动时间20、35、10小时,夏季劳动时间30、75、 40小时,年景收益分别为175、300、120美元。 建立数学模型,帮助该农户确定养殖计划,使得 年净收入最多。
d=5 2 6 10 4 5 7 5 7 6 4 4 9 10 6 2; enddata min=@sum(link(i,j):2*d(i,j)*x(i,j))+@su m(demand(i):5*x(4,i)); @for(demand(j):@sum(supply(i):x(i,j))>= a(j)); @for(supply(i)|i#LT#4:@sum(demand(j):x( i,j))<=b(i)); @for(link(i,j):x(i,j)>=0); end
数学模型与实验(六)
关于LINDO/LINGO软件
LINDO 公司软件产品简要介绍
美国芝加哥(Chicago)大学的Linus Schrage教授于1980 年前后开发, 后来成立 LINDO系统公司(LINDO Systems Inc.), 网址:
LINDO: Linear INteractive and Discrete Optimizer LINGO: Linear INteractive General Optimizer What’s Best!: (SpreadSheet e.g. EXCEL) (V6.1) (V8.0) (V7.0)
Lingo程序(初等形式)
model: min=99*x11+60*x12+59*x13+73*x14+79*x21+65*x22+ 93*x23+87*x24+67*x31+93*x32+63*x33+81*x34+56* x41+79*x42+86*x43+76*x44; x11+x12+x13+x14=1; x21+x22+x23+x24=1; x31+x32+x33+x34=1; x41+x42+x43+x44=1; x11+x21+x31+x41=1; x12+x22+x32+x42=1; x13+x23+x33+x43=1; x14+x24+x34+x44=1; end
@for(ren(i):@sum(yundong(j):x(i,j))=1);
@for(yundong(j):@sum(ren(i):x(i,j))=1);
end
例2 一家大建筑公司正在三个地点开掘。同 时又在其他四个地点建筑,这里需要土方的 填充。在1、2、3处挖掘产生的土方分别为每 天150,400,325立方码。建筑地点A、B、C、 D处需要的填充土方分别为175,125,225, 450立方码。也可以从地点4用每立方码5美元 的价格获得额外的填充土方。填充土方运输 的费用约为一货车容量每英里20美元。一辆 货车可以搬运10立方码的土方(每立方码土方 每英里运输费2美元)。表3-3给出了各地点间 距离的英里数。求使公司花费最少的运输计 划。
例1 混合泳接力赛由蛙泳、蝶泳、自由泳、仰泳
组成。如何根据 4位运动员的4种游泳竞赛成绩 安排混合泳接力队,以取得最佳成绩。 蛙泳 蝶泳 自由泳 仰泳 甲 99 乙 79 60 59 65 93 93
63
81
76
79 86
这就是简单的0-1规划。取决策变量 蛙泳 蝶泳 自由泳 仰泳 甲 X11 X12 X13 X14 乙 X21 X22 X23 X24 丙 X31 X32 X33 X34 丁 X41 X42 X43 X44 Xij=1 让第i人参加第j项比赛,否则 Xij=0
如:尽量少使用绝对值、符号函数、多个变量求 最大/最小值、四舍五入、取整函数等 3、尽量使用线性模型,减少非线性约束和非线性变 量的个数 (如x/y <5 改为x<5y)
4、合理设定变量上下界,尽可能给出变量初始值
5、模型中使用的参数数量级要适当 (如小于103)
需要掌握的几个重要方面
1、LINDO:
习题5
农作物问题
某农户有100英亩土地合5000美元可供投资。每年 冬季家庭成员可以贡献3500小时的劳动时间,而 夏季为4000小时。如果这些劳动时间有富裕,家 庭成员可以去附近农场打工,冬季每小时4.8美元, 夏季每小时5.1美元。 现金收入来源于3种农作物(大豆、玉米、燕麦) 以及2种家禽(奶牛、母鸡)。农作物不需要投资, 但每头奶牛需要400美元初始投资,每只母鸡需要 3美元初始投资。
习题 4
资源的最优配置策略
某工厂有1000台机器, 生产两种产品 A, B, 若投入 y 台机器生产A 产品, 则纯收入为 5y .若投入 y 台机器生产B 产品, 则纯收 入为 4y . 又知, 生产A 种产品机器的年折 损率为20%, , 生产B 种产品机器的年折损 率为10%, 问在5年内如何安排各年度的生产 计划, 才能使总收入最高.
挖掘与建筑地点间的距离(英里) 接收填充土方的地点 挖掘地点 A B C D 1 5 2 6 10 2 4 5 7 5 3 7 6 4 4 4 9 10 6 2
参量、变量 第 j 个工地需要填充土方量为 aj, 第 i 个挖掘点产生的土方量为 bi ,
第i个挖掘点到第 j 个工地的路程为 dij, 计划运输
LINDO API: LINDO Application Programming Interface (V2.0)
演示(试用)版、学生版、高级版、超级版、工业版、 扩展版… (求解问题规模和选件不同)
建模时需要注意的几个基本问题
1、尽量使用实数优化,减少整数约束和整数变量
2、尽量使用光滑优化,减少非光滑约束的个数
练习7* 2005年D题: DVD在线租赁
考虑如下的在线DVD租赁问题。顾客缴纳一定数量 的月费成为会员,订购DVD租赁服务。会员对哪些 DVD有兴趣,只要在线提交订单,网站就会通过快 递的方式尽可能满足要求。会员提交的订单包括多 张DVD,这些DVD是基于其偏爱程度排序的。 问题2 表2中列出了网站手上20种DVD的现有张数和 当前需要处理的100位会员的在线订单,如何对这些 DVD进行分配,才能使会员获得最大的满意度?请 具体列出前30位会员(即C0001~C0030)分别获得 哪些DVD。
的土方量为xij.
如何安排运输计划使总费用最小.
模型:
Min= ij (2 dij xij+ j 5 x4j) i xij = aj j=1,2,3,4. j xij = bi i=1,2,3.
xij >=0
Lingo程序(模式语言) model: sets: demand/1..4/:a; supply/1..4/:b; link(supply,demand):d,x; endsets data: a=175,125,225,450; b=150,400,325,1;
Lingo程序(模式语言)
model: sets: ren/1..4/:r; yundong/1..4/:y; link(ren, yundong):a,x; endsets
data: A=99 60 59 87 73 79 65 93
67 93 63
56 79 86
81
76;
enddata
min=@sum(link(i,j):a(i,j)*x(i,j)) ;
正确阅读求解报告(尤其要掌握敏感性分析)
2、LINGO: 掌握集合(SETS)的应用; 正确阅读求解报告; 正确理解求解状态窗口; 学会设置基本的求解选项(OPTIONS) ; 掌握与外部文件的基本接口方法
Lingo程序(初等形式)
用Lingo软件求解线性规划 Model: Max=110*x1+75*x2+60*x3; x1+x2+x3<=50; 1/2*x1+1/3*x2+1/4*x3<=20; end
练习6*:2004年D题 公务员招聘
经笔试科目:综合基础知识、专业知识和 “行政职业能力测验”三个部分, 每科满 分为100分。16人被选择进入第二阶段的面 试考核。面试考核从高到低分成A/B/C/D四 个等级。由招聘领导小组综合笔初试成绩、 面试成绩、以及各用人部门需求确定录用 名单,并分配到各用人部门。 问题 1. 如果不考虑应聘人员的意愿,择优 按需录用,试帮助招聘领导小组设计一种 录用分配方案。