运筹学实验报告-lingo软件的使用-习题代码
LINGO软件学习入门实验报告
LINGO 实验报告一.实验目的1、熟悉LINGO 软件的使用方法、功能；2、学会用LINGO 软件求解一般的线性规划问题。
二.实验内容1、求解线性规划：12121212max z x 2x 2x 5x 12s.t.x 2x 8x ,x 0=++≥⎧⎪+≤⎨⎪≥⎩2、求解线性规划：12121212min z 20x 10x 5x 4x 24s.t.2x 5x 5x ,x 0=++≤⎧⎪+≥⎨⎪≥⎩3、假设现在一个计算机厂商要生产两种型号的PC ：标准型(standard)和增强型(turbo)，由于生产线和劳动力工作时间的约束，使得标准型PC 最多生产100台。
增强型PC 最多生产120台；一共耗时劳动力时间不能超过160小时。
已知每台标准型PC 可获利润$100，耗掉1小时劳动力工作时间；每台增强型PC 可获利润$150，耗掉2小时劳动力工作时间。
请问：该如何规划这两种计算机的生产量才能够使得最后获利最大？三. 模型建立1、求解线性规划：12121212max z x 2x 2x 5x 12s.t.x 2x 8x ,x 0=++≥⎧⎪+≤⎨⎪≥⎩2、求解线性规划：12121212min z 20x 10x 5x 4x 24s.t.2x 5x 5x ,x 0=++≤⎧⎪+≥⎨⎪≥⎩3、设生产标准型为1x 台；生产增强型2x 台，则可建立线性规划问题数学模型为12121212max z 100x 150x x 100x 120s.t.x 2x 160x ,x 0=+≤⎧⎪≤⎪⎨+≤⎪⎪≥⎩四. 模型求解（含经调试后正确的源程序）1、求解线性规划：model:max=x1+2*x2;2*x1+5*x2>12;x1+2*x2<8;end结果显示：2、求解线性规划：model:min=20*x1+10*x2;5*x1+4*x2<24;2*x1+5*x2>5;End结果显示：3、求解线性规划：model:mAX=100*x1+150*x2; x1+2*x2<160;x1<100;x2<120;end结果显示：五．结果分析对于第一题中我们得出最优解：x1=0;x2=4;最优值max=8;第二题中我们得出最优解：x1=0;x2=1;最优值min=10;第三题中我们得出最优解：x1=1000;x2=30;最优值max=14500;即：生产标准型100台，生产增强型30台时，使得最后获利达最大，为14500。
运筹学上机实践报告LINGO软件
Southwest university of science and technology实验报告LINGO软件在线性规划中的运用学院名称环境与资源学院专业名称采矿工程学生姓名学号____________________________________ 指导教师陈星明教授二◦一五年十一月实验LINGO软件在线性规划中的运用实验目的掌握LINGO软件求解线性规划问题的基本步骤，了解LINGO软件解决线性规划问题的基本原理，熟悉常用的线性规划计算代码，理解线性规划问题的迭代关系。
实验仪器、设备或软件电脑，LINGO软件实验内容1. LINGO软件求解线性规划问题的基本原理；2•编写并调试LINGO软件求解线性规划问题的计算代码；实验步骤1•使用LINGO计算并求解线性规划问题；2 •写出实验报告，并浅谈学习心得体会（线性规划的基本求解思路与方法及求解过程中出现的问题及解决方法）。
实验过程有一艘货轮，分为前、中、后三个舱位，它们的容积与允许载重量如下表所示。
现有三种商品待运，已知有关数据列于下表中。
又为了航运安全，要求前、中、后舱在实际载重量上大体保持各舱最大允许载重量的比例关系。
具体要求前、后舱分别与中舱之间的载重量比例偏差不超过15%,前、后舱之间不超过10%。
问货轮首先分析问题，建立数学模型:确定决策变量假设i=1,2,3分别代表商品A、B C, 8用j=1,2,3分别代表前、中、后舱，设决策变量X ij为装于j舱位的第i种商品的数量（件）。
确定目标函数商品A的件数为：x11- x12x13商品B的件数为：x21x22x23商品A的件数为：X31 - X32 - X33为使运费最高，目标函数为：确定约束条件前、中、后舱位载重限制为：前、中、后舱位体积限制为：A、B、C三种商品数量的限制条件：各舱最大允许载重量的比例关系构成的约束条件：且决策变量要求非负，即X j > 0,i=1,2；j=1,2,3。
运筹学实验
《运筹学》上机实验报告学 院 机电工程学院 专 业 工业工程 指导教师 吴小东 班 级 工业10-1班 学生姓名 林 金 铎 学生学号实验时间 2012-2013学年第一学期实验一 使用LINGO 求解线性规划问题班级：工业10-1班 姓名：林金铎 学号： 评阅成绩： 已知如下线性规划模型：123max 303540z x x x =++1231231231233251823412229,,0x x x x x x x x x x x x ++≤⎧⎪++≤⎪⎨++≤⎪⎪≥⎩ 一、利用集的方法编写上述线性规划模型的LINGO 程序。
在LINGO 软件模型中编写本题的程序如下图1-1所示所示。
图1-1 LINGO 模型窗口截图点击LINGO 菜单下的Solve 选项，LINGO 软件求解所输入的模型，得到LINGO 运行状态窗口如图1-2所示图1-2 LINGO运行状态窗口截图运行结束后，关闭LINGO运行状态窗口，获得LINGO软件的结果报告窗口，如图1-3、1-4所示。
图1-3 LINGO结果报告窗口截图（一）图1-4 LINGO结果报告窗口截图（二）二、根据编写的程序，回答以下问题：1、哪些是原始集答：var(j), const(i)是原始集2、哪个是派生集该派生集是稠密集还是稀疏集该派生集有多少个成员答：A(i,j)是派生集,属于稠密集合，共有9个成员3、属性值“5”是属于成员(b1,x3)还是(b3,x1)的属性值答：属于成员(b1,x3)的属性值三、根据程序的运行结果，回答以下问题：1、全局最优值是否已经找到该值是多少答：已经找到，最优值为1652、该模型求解一共迭代了多少次答：共迭代了2次3、在求解结果的界面中，Variable、Value、Reduced Cost、Row、Slack or Surplus 和Dual Price分别表示什么答：Variable表示运算时各定义变量的取值；Value表示给出最优解中各变量的值；Reduced Cost表示列出最优单纯形表中判别数所在行的变量的系数，表示当变量有微小变动时, 目标函数的变化率；Row表示行数；Slack or Surplus 表示给出松驰变量的值；Dual Price表示当对应约束有微小变动时, 目标函数的变化率。
lingo软件使用
一、基本语法--- Lingo的内部函数
◆ 关系运算符 “=” 表示等于关系 “<=” 表示小于等于关系 “>=” 表示大于等于关系
◆ 数学函数 @abs(x) 返回x的绝对值
@sin(x) 返回x的正弦值，x采用弧度制
7 i 1
7
5
154 x ij s ij x ij q j j 1, 2 , , 5 i 1, 2 , , 7 i 1, 2 , , 7 j 1, 2 , , 5 j 1, 2 , , 5
i 1 j 1
Hale Waihona Puke x i 1 x i 2 x i 5 ks i x i 3 x i 4 ys i
一、基本语法---程序结构的三个部分
◆ 目标函数与约束条件部分 按模型中的目标函数及约束条件逐句写入 格式：max = …… 或 min =……
利用lingo中的内部函数组织目标函数及约束
条件的具体内容
例 sets: factory/1..60/:cost, volume; endsets min=@sum(factory(i): cost(i)*volume(i));
返回x1，x2，…，xn中的最大值
@smin(x1,x2,…,xn)
返回x1，x2，…，xn中的最小值
一、基本语法--- Lingo的内部函数
◆ 集循环函数 集循环函数遍历整个集进行操作。其语法为 @function(setname[(set_index_list)
集循环函数名 集合名 循环控制变量
◆ 输入和输出函数 输入和输出函数可以把程序和外部数据文件， 如文本文件、电子表格等连接起来。 @file 格式：@file(’文件名’) 功能：从指定的外部文本文件中输入数据 在数据文件中，以”~”作为一批数据结束的标 记，同一文件中，可以含多批数据，系统依次 打开数据文件时，依次读取各批数据
运筹学实验报告lingo
二. 实验题目
1、求解线性规划：
max
z x 1 2x
2
2x 1 5x 2 12 s.t. x 1 2x 2 8 x , x 0 2 1
并对价值系数、右端常量进行灵敏度分析。
2、已知某工厂计划生产I，II，III三种产品，各 产品需要在A、B、C设备上加工，有关数据如下：
Allowable Decrease：允许减少量
Current RHS ：当前右边常数项
结论1：
该线性规划问题的最优解为：X*=(35,10,0)T 最优值为：z*=215
结论2：
c1＝5 c1在（4,8）内原最优解不变，但最优值是要变的 c2＝4 c2在（2.7,5）内原最优解不变，但最优值是要变的 c3＝3 c3在（-∞ ,7）内原最优解，最优值都是不变的 b1=45 b1在（40, 50）内原最优基不变，但最优解和最优值是要变的 b2=80 b2在（67.5, 90）内原最优基不变，但最优解和最优值是要变的 b3=90 b3在（65, ∞ ）内原最优基不变，但最优解和最优值是要变的
Row 1 2 3 4 Slack or Surplus 215.0000 0.000000 0.000000 25.00000 Dual Price 1.000000 3.000000 1.000000 0.000000
激活灵敏度计算功能
法一：打开command window，输入range；
法二：LINGO——options —— General Solver —Dual Computations——Prices & Ranges
LINGO
Outline
一.熟悉LINDO软件的灵敏度分析功能
运筹学实验报告-lingo软件的使用-习题代码
运筹学实验报告姓名：学号：班级：相关问题说明：一、实验性质和教学目的本实验是运筹学课内安排的上机操作实验。
目的在于了解、熟悉计算机Lingo软件在运筹学模型求解中的作用，激发学习兴趣，提高学习效果，增强自身的动手能力，提高实际应用能力。
二、实验基本要求要求学生：1. 实验前认真做好理论准备，仔细阅读实验指导书；2. 遵从教师指导，认真完成实验任务，按时按质提交实验报告。
三、主要参考资料1．LINGO软件2. LINGO8.0及其在环境系统优化中的应用，天津大学出版社，20053. 优化建模与LINDO/LINGO软件，清华大学出版社，20054．运筹学编写组主编，运筹学（修订版），清华大学出版社，19905．蓝伯雄主编，管理数学（下）—运筹学，清华大学出版社，19976．胡运权主编，运筹学习题集（修订版），清华大学出版社，19957．胡运权主编，运筹学教程（第二版），清华大学出版社，2003实验内容1、线性规划问题：⎪⎪⎩⎪⎪⎨⎧≥≤+≤+≤++=0,13119241171289..68max 2121212121x x x x x x x x t s x x z (1) 给出原始代码；(2) 计算结果(包括灵敏度分析，求解结果粘贴)；(3) 回答下列问题(手写)：a ) 最优解及最优目标函数值是多少；b ) 资源的对偶价格各为多少，并说明对偶价格的含义；c ) 为了使目标函数值增加最多，让你选择一个约束条件，将它的常数项增加一个单位，你将选择哪一个约束条件？这时目标函数值将是多少？d ) 对x 2的目标函数系数进行灵敏度分析；e ) 对第2个约束的约束右端项进行灵敏度分析；f ) 结合本题的结果解释“Reduced Cost ”的含义。
对偶价格就是说 约束方程右端变量增加1对目标函数值的影响 答案： （1）代码max =8*x1+6*x2; 9*x1+8*x2<=12; 7*x1+11*x2<=24; 9*x1+11*x2<=13; x1>=0; x2>=0;（2）计算结果Global optimal solution found.Objective value: 10.66667 Total solver iterations: 2Variable Value Reduced Cost X1 1.333333 0.000000 X2 0.000000 1.111111Row Slack or Surplus Dual Price 1 10.66667 1.000000 2 0.000000 0.8888889 3 14.66667 0.000000 4 1.000000 0.000000 5 1.333333 0.0000006 0.000000 0.000000 Ranges in which the basis is unchanged:Objective Coefficient RangesCurrent Allowable Allowable Variable Coefficient Increase Decrease X1 8.000000 INFINITY 1.250000 X2 6.000000 1.111111 INFINITYRighthand Side RangesRow Current Allowable Allowable RHS Increase Decrease 2 12.00000 1.000000 12.00000 3 24.00000 INFINITY 14.66667 4 13.00000 INFINITY 1.000000 5 0.0 1.333333 INFINITY 6 0.0 0.0 INFINITY(3)a) b) c) d) e) f)2、运输问题：(1) 给出原始代码；(2) 计算结果(决策变量求解结果粘贴)Min Z = Cij Xij∑=61i Xij <=bj （j=1...8） 销量约束∑∑==6181i j∑=81j Xij = ai （i=1...6） 产量约束Xij ≥ 0（i=1...6；j=1...8）代码：model :!6发点8 model :!6发点8收点运输问题; sets :warehouses/wh1..wh6/: capacity; vendors/v1..v8/: demand;links(warehouses,vendors): cost, volume; endsetsmin =@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)); !产量约束; !这里是数据; data :capacity=55 47 42 52 41 32; demand=60 55 51 43 41 52 43 38; cost=6 2 9 7 4 2 5 9 4 5 5 3 8 5 3 25 2 1 3 7 4 8 3 767 9 9 2 7 1 2 3 6 5 7 2 6 5 5 9 2 28 1 4 3; enddata end答案Global optimal solution found.Objective value: 473.0000 Infeasibilities: 0.000000 Total solver iterations: 9Model Class: LPTotal variables: 48 Nonlinear variables: 0 Integer variables: 0Total constraints: 15Nonlinear constraints: 0Total nonzeros: 144Nonlinear nonzeros: 0Variable Value Reduced Cost CAPACITY( WH1) 55.00000 0.000000 CAPACITY( WH2) 47.00000 0.000000 CAPACITY( WH3) 42.00000 0.000000 CAPACITY( WH4) 52.00000 0.000000 CAPACITY( WH5) 41.00000 0.000000 CAPACITY( WH6) 32.00000 0.000000 DEMAND( V1) 60.00000 0.000000 DEMAND( V2) 55.00000 0.000000 DEMAND( V3) 51.00000 0.000000 DEMAND( V4) 43.00000 0.000000 DEMAND( V5) 41.00000 0.000000 DEMAND( V6) 52.00000 0.000000 DEMAND( V7) 43.00000 0.000000 DEMAND( V8) 38.00000 0.000000 COST( WH1, V1) 6.000000 0.000000 COST( WH1, V2) 2.000000 0.000000 COST( WH1, V3) 9.000000 0.000000 COST( WH1, V4) 7.000000 0.000000 COST( WH1, V5) 4.000000 0.000000 COST( WH1, V6) 2.000000 0.000000 COST( WH1, V7) 5.000000 0.000000 COST( WH1, V8) 9.000000 0.000000 COST( WH2, V1) 4.000000 0.000000 COST( WH2, V2) 5.000000 0.000000 COST( WH2, V3) 5.000000 0.000000 COST( WH2, V4) 3.000000 0.000000 COST( WH2, V5) 8.000000 0.000000 COST( WH2, V6) 5.000000 0.000000 COST( WH2, V7) 3.000000 0.000000 COST( WH2, V8) 2.000000 0.000000 COST( WH3, V1) 5.000000 0.000000 COST( WH3, V2) 2.000000 0.000000 COST( WH3, V3) 1.000000 0.000000 COST( WH3, V4) 3.000000 0.000000 COST( WH3, V5) 7.000000 0.000000 COST( WH3, V6) 4.000000 0.000000COST( WH3, V8) 3.000000 0.000000 COST( WH4, V1) 7.000000 0.000000 COST( WH4, V2) 6.000000 0.000000 COST( WH4, V3) 7.000000 0.000000 COST( WH4, V4) 9.000000 0.000000 COST( WH4, V5) 9.000000 0.000000 COST( WH4, V6) 2.000000 0.000000 COST( WH4, V7) 7.000000 0.000000 COST( WH4, V8) 1.000000 0.000000 COST( WH5, V1) 2.000000 0.000000 COST( WH5, V2) 3.000000 0.000000 COST( WH5, V3) 6.000000 0.000000 COST( WH5, V4) 5.000000 0.000000 COST( WH5, V5) 7.000000 0.000000 COST( WH5, V6) 2.000000 0.000000 COST( WH5, V7) 6.000000 0.000000 COST( WH5, V8) 5.000000 0.000000 COST( WH6, V1) 5.000000 0.000000 COST( WH6, V2) 9.000000 0.000000 COST( WH6, V3) 2.000000 0.000000 COST( WH6, V4) 2.000000 0.000000 COST( WH6, V5) 8.000000 0.000000 COST( WH6, V6) 1.000000 0.000000 COST( WH6, V7) 4.000000 0.000000 COST( WH6, V8) 3.000000 0.000000 VOLUME( WH1, V1) 0.000000 4.000000 VOLUME( WH1, V2) 55.00000 0.000000 VOLUME( WH1, V3) 0.000000 7.000000 VOLUME( WH1, V4) 0.000000 5.000000 VOLUME( WH1, V5) 0.000000 2.000000 VOLUME( WH1, V6) 0.000000 0.000000 VOLUME( WH1, V7) 0.000000 3.000000 VOLUME( WH1, V8) 0.000000 8.000000 VOLUME( WH2, V1) 0.000000 1.000000 VOLUME( WH2, V2) 0.000000 2.000000 VOLUME( WH2, V3) 0.000000 2.000000 VOLUME( WH2, V4) 43.00000 0.000000 VOLUME( WH2, V5) 0.000000 5.000000 VOLUME( WH2, V6) 0.000000 2.000000 VOLUME( WH2, V7) 4.000000 0.000000 VOLUME( WH2, V8) 0.000000 0.000000 VOLUME( WH3, V1) 0.000000 4.000000 VOLUME( WH3, V2) 0.000000 1.000000VOLUME( WH3, V4) 0.000000 2.000000 VOLUME( WH3, V5) 0.000000 6.000000 VOLUME( WH3, V6) 0.000000 3.000000 VOLUME( WH3, V7) 0.000000 7.000000 VOLUME( WH3, V8) 0.000000 3.000000 VOLUME( WH4, V1) 0.000000 5.000000 VOLUME( WH4, V2) 0.000000 4.000000 VOLUME( WH4, V3) 0.000000 5.000000 VOLUME( WH4, V4) 0.000000 7.000000 VOLUME( WH4, V5) 0.000000 7.000000 VOLUME( WH4, V6) 14.00000 0.000000 VOLUME( WH4, V7) 0.000000 5.000000 VOLUME( WH4, V8) 38.00000 0.000000 VOLUME( WH5, V1) 41.00000 0.000000 VOLUME( WH5, V2) 0.000000 1.000000 VOLUME( WH5, V3) 0.000000 4.000000 VOLUME( WH5, V4) 0.000000 3.000000 VOLUME( WH5, V5) 0.000000 5.000000 VOLUME( WH5, V6) 0.000000 0.000000 VOLUME( WH5, V7) 0.000000 4.000000 VOLUME( WH5, V8) 0.000000 4.000000 VOLUME( WH6, V1) 0.000000 4.000000 VOLUME( WH6, V2) 0.000000 8.000000 VOLUME( WH6, V3) 0.000000 1.000000 VOLUME( WH6, V4) 0.000000 1.000000 VOLUME( WH6, V5) 0.000000 7.000000 VOLUME( WH6, V6) 32.00000 0.000000 VOLUME( WH6, V7) 0.000000 3.000000 VOLUME( WH6, V8) 0.000000 3.000000Row Slack or Surplus Dual Price1 473.0000 -1.0000002 19.00000 0.0000003 0.000000 0.0000004 9.000000 0.0000005 0.000000 0.0000006 41.00000 0.0000007 6.000000 0.0000008 39.00000 0.0000009 0.000000 1.00000010 0.000000 -2.00000011 0.000000 -3.00000012 0.000000 -1.00000013 0.000000 -2.00000014 0.000000 -2.00000015 0.000000 -1.0000003、一般整数规划问题：某服务部门各时段（每2h为一时段）需要的服务员人数见下表。
lingo实验报告
lingo实验报告《lingo 实验报告》一、实验目的本次 lingo 实验旨在深入了解和掌握 Lingo 软件在优化问题求解中的应用，通过实际操作和案例分析，提高对数学建模和优化算法的理解与运用能力。
二、实验环境本次实验使用的计算机配置为：处理器_____，内存_____，操作系统_____。
Lingo 软件版本为_____。
三、实验内容（一）线性规划问题1、问题描述考虑一个生产计划问题，某工厂生产两种产品 A 和 B，生产 A 产品每个需要消耗 2 个单位的原材料和 3 个单位的工时，生产 B 产品每个需要消耗 3 个单位的原材料和 2 个单位的工时。
工厂现有原材料 100 个单位，工时 80 个单位。
A 产品的单位利润为 5 元，B 产品的单位利润为 4 元。
求如何安排生产计划，使得总利润最大。
2、数学模型设生产 A 产品 x 个，生产 B 产品 y 个，则目标函数为：Max Z ＝5x ＋ 4y约束条件为：2x ＋ 3y ＜＝ 1003x ＋ 2y ＜＝ 80x ＞＝ 0y ＞＝ 03、 Lingo 代码及求解结果｀｀｀lingomodel:max ＝ 5x ＋ 4y;2x ＋ 3y ＜＝ 100;3x ＋ 2y ＜＝ 80;x ＞＝ 0;y ＞＝ 0;end｀｀｀求解结果为：x ＝ 20，y ＝ 20，最大利润为 180 元。
（二）整数规划问题1、问题描述某物流公司需要安排车辆运输货物，有两种车型可供选择，大型车每辆可载货10 吨，小型车每辆可载货5 吨。
共有货物80 吨需要运输，且大型车的数量不能超过 8 辆，小型车的数量不能超过 12 辆。
大型车每辆的运输成本为 100 元，小型车每辆的运输成本为 60 元。
求如何安排车辆，使得运输成本最低。
2、数学模型设安排大型车 x 辆，小型车 y 辆，则目标函数为：Min C ＝ 100x ＋60y约束条件为：10x ＋ 5y ＞＝ 80x ＜＝ 8y ＜＝ 12x ＞＝ 0 且为整数y ＞＝ 0 且为整数3、 Lingo 代码及求解结果｀｀｀lingomodel:min ＝ 100x ＋ 60y;10x ＋ 5y ＞＝ 80;x ＜＝ 8;y ＜＝ 12;＠gin(x)；＠gin(y)；end｀｀｀求解结果为：x ＝ 6，y ＝ 4，最低运输成本为 840 元。
LINGO实验项目
实验项目一线性规划实验学时：2实验目的：线性规划（Linear Programming，简写LP）是运筹学中最成熟的一个分枝，而且是应用最为广泛的一个运筹学分枝，是解决最优化问题的重要工具。
而目前 Lindo/lingo 是求解线性规划比较成熟的一个软件，通过本实验，掌握线性规划模型在 Lindo/lingo 中的求解，并能达到灵活运用。
实验要求：1.掌握线性规划的建模步骤及方法；2.掌握Lindo/lingo 的初步使用；3.掌握线性规划模型在Lindo/lingo 建模及求解；4.掌握线性规划的灵敏度分析实验内容及步骤：例：美佳公司计划制造I、II 两种家电产品。
已知各制造一件时分别占用设备A、B 的台时、调试时间、调试工序每天可用于这种家电的能力、各售出一件时的获利情况，如表1-1 所示。
1.问该公司应制造两种家电各多少件，使其获取的利润最大。
2. 如果资源出租，资源出租的最低价格至少是多少（即每种资源的影子价格是多少）。
3.若家电I 的利润不变，家电II 的利润在什么范围内变化时，则该公司的最优生产计划将不发生变化。
4. 若设备A 和B 每天可用能力不变，则调试工序能力在什么范围内变化时，问题的最优基不变。
解：设x1表示产品I 的生产量; x2表示产品II 的生产量,所在该线性规划的模型为：从此线性规划的模型中可以看出，第一个小问是典型的生产计划问题，第二小问是相应资源的影子价格，第三和第四个小问则是此问题的灵敏度分析。
现在我们利用lingo8.0 来教你求解线性规划问题。
第一步，启动lingo 进入初始界面如下图1-1 和图1-2 所示：第二步，在进行线性规划模型求解时，先要对初始求解方法及参数要进行设置，首先选择lingo 菜单下的Option 菜单项，并切换在general solver（通用求解器）页面下，如下图1-3 所示：general solver 选项卡上的各项设置意义如下表格1-1 所示：表格1-1 general solver 选项卡上的各项设置意义接下来再对Linear Solver（线性求解器）选项卡进行设置，切换界面如所示：其各项设置意义如下表格1-2 所示：表格1-2 Linear Solver 选项卡各项设置意义因为这个线性规划模型较为简单，数字也是比较小的，而且需要进行灵敏度分析，所以对general solver 选项卡上的Dual Computations（对偶计算）项设为“Prices and Ranges（计算对偶价格并分析敏感性）”。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
运筹学实验报告姓名：学号：班级：相关问题说明：一、实验性质和教学目的本实验是运筹学课安排的上机操作实验。
目的在于了解、熟悉计算机Lingo软件在运筹学模型求解中的作用，激发学习兴趣，提高学习效果，增强自身的动手能力，提高实际应用能力。
二、实验基本要求要求学生：1. 实验前认真做好理论准备，仔细阅读实验指导书；2. 遵从教师指导，认真完成实验任务，按时按质提交实验报告。
三、主要参考资料1．LINGO软件2. LINGO8.0及其在环境系统优化中的应用，大学，20053. 优化建模与LINDO/LINGO软件，清华大学，20054．运筹学编写组主编，运筹学（修订版），清华大学，19905．蓝伯雄主编，管理数学（下）—运筹学，清华大学，19976．胡运权主编，运筹学习题集（修订版），清华大学，19957．胡运权主编，运筹学教程（第二版），清华大学，2003实验容1、线性规划问题：⎪⎪⎩⎪⎪⎨⎧≥≤+≤+≤++=0,13119241171289..68max 2121212121x x x x x x x x t s x x z (1) 给出原始代码；(2) 计算结果(包括灵敏度分析，求解结果粘贴)；(3) 回答下列问题(手写)：a ) 最优解及最优目标函数值是多少；b ) 资源的对偶价格各为多少，并说明对偶价格的含义；c ) 为了使目标函数值增加最多，让你选择一个约束条件，将它的常数项增加一个单位，你将选择哪一个约束条件？这时目标函数值将是多少？d ) 对x 2的目标函数系数进行灵敏度分析；e ) 对第2个约束的约束右端项进行灵敏度分析；f ) 结合本题的结果解释“Reduced Cost ”的含义。
对偶价格就是说 约束方程右端变量增加1对目标函数值的影响答案：（1）代码max =8*x1+6*x2;9*x1+8*x2<=12;7*x1+11*x2<=24;9*x1+11*x2<=13;x1>=0;x2>=0;（2）计算结果Global optimal solution found.Objective value: 10.66667Total solver iterations: 2Variable Value Reduced Cost X1 1.333333 0.000000 X2 0.000000 1.111111Row Slack or Surplus Dual Price 1 10.66667 1.000000 2 0.000000 0.8888889 3 14.66667 0.000000 4 1.000000 0.0000006 0.000000 0.000000Ranges in which the basis is unchanged:Objective Coefficient RangesCurrent Allowable Allowable Variable Coefficient Increase Decrease X1 8.000000 INFINITY 1.250000 X2 6.000000 1.111111 INFINITYRighthand Side RangesRow Current Allowable Allowable RHS Increase Decrease 2 12.00000 1.000000 12.00000 3 24.00000 INFINITY 14.66667 4 13.00000 INFINITY 1.000000 5 0.0 1.333333 INFINITY 6 0.0 0.0 INFINITY(3)a)b)c)d)e)f)2、运输问题：(1) 给出原始代码；(2) 计算结果(决策变量求解结果粘贴)Min Z = Cij Xij∑=61i Xij <=bj （j=1...8） 销量约束∑∑==6181i j∑=81j Xij = ai （i=1...6） 产量约束Xij ≥ 0（i=1...6；j=1...8）代码：model :!6发点8 model :!6发点8收点运输问题;sets :warehouses/wh1..wh6/: capacity;vendors/v1..v8/: demand;links(warehouses,vendors): cost, volume;endsetsmin =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)); !产量约束;!这里是数据;data :capacity=55 47 42 52 41 32;demand=60 55 51 43 41 52 43 38;cost=6 2 9 7 4 2 5 94 5 5 3 8 5 3 25 2 1 3 7 4 8 37 6 7 9 9 2 7 12 3 6 5 7 2 6 55 9 2 2 8 1 4 3;enddataend答案Global optimal solution found.Objective value: 473.0000Infeasibilities: 0.000000Total solver iterations: 9Model Class: LPTotal variables: 48Nonlinear variables: 0Total constraints: 15Nonlinear constraints: 0Total nonzeros: 144Nonlinear nonzeros: 0Variable Value Reduced Cost CAPACITY( WH1) 55.00000 0.000000 CAPACITY( WH2) 47.00000 0.000000 CAPACITY( WH3) 42.00000 0.000000 CAPACITY( WH4) 52.00000 0.000000 CAPACITY( WH5) 41.00000 0.000000 CAPACITY( WH6) 32.00000 0.000000 DEMAND( V1) 60.00000 0.000000 DEMAND( V2) 55.00000 0.000000 DEMAND( V3) 51.00000 0.000000 DEMAND( V4) 43.00000 0.000000 DEMAND( V5) 41.00000 0.000000 DEMAND( V6) 52.00000 0.000000 DEMAND( V7) 43.00000 0.000000 DEMAND( V8) 38.00000 0.000000 COST( WH1, V1) 6.000000 0.000000 COST( WH1, V2) 2.000000 0.000000 COST( WH1, V3) 9.000000 0.000000 COST( WH1, V4) 7.000000 0.000000 COST( WH1, V5) 4.000000 0.000000 COST( WH1, V6) 2.000000 0.000000 COST( WH1, V7) 5.000000 0.000000 COST( WH1, V8) 9.000000 0.000000 COST( WH2, V1) 4.000000 0.000000 COST( WH2, V2) 5.000000 0.000000 COST( WH2, V3) 5.000000 0.000000 COST( WH2, V4) 3.000000 0.000000 COST( WH2, V5) 8.000000 0.000000 COST( WH2, V6) 5.000000 0.000000 COST( WH2, V7) 3.000000 0.000000 COST( WH2, V8) 2.000000 0.000000 COST( WH3, V1) 5.000000 0.000000 COST( WH3, V2) 2.000000 0.000000 COST( WH3, V3) 1.000000 0.000000 COST( WH3, V4) 3.000000 0.000000 COST( WH3, V5) 7.000000 0.000000 COST( WH3, V6) 4.000000 0.000000COST( WH4, V1) 7.000000 0.000000 COST( WH4, V2) 6.000000 0.000000 COST( WH4, V3) 7.000000 0.000000 COST( WH4, V4) 9.000000 0.000000 COST( WH4, V5) 9.000000 0.000000 COST( WH4, V6) 2.000000 0.000000 COST( WH4, V7) 7.000000 0.000000 COST( WH4, V8) 1.000000 0.000000 COST( WH5, V1) 2.000000 0.000000 COST( WH5, V2) 3.000000 0.000000 COST( WH5, V3) 6.000000 0.000000 COST( WH5, V4) 5.000000 0.000000 COST( WH5, V5) 7.000000 0.000000 COST( WH5, V6) 2.000000 0.000000 COST( WH5, V7) 6.000000 0.000000 COST( WH5, V8) 5.000000 0.000000 COST( WH6, V1) 5.000000 0.000000 COST( WH6, V2) 9.000000 0.000000 COST( WH6, V3) 2.000000 0.000000 COST( WH6, V4) 2.000000 0.000000 COST( WH6, V5) 8.000000 0.000000 COST( WH6, V6) 1.000000 0.000000 COST( WH6, V7) 4.000000 0.000000 COST( WH6, V8) 3.000000 0.000000 VOLUME( WH1, V1) 0.000000 4.000000 VOLUME( WH1, V2) 55.00000 0.000000 VOLUME( WH1, V3) 0.000000 7.000000 VOLUME( WH1, V4) 0.000000 5.000000 VOLUME( WH1, V5) 0.000000 2.000000 VOLUME( WH1, V6) 0.000000 0.000000 VOLUME( WH1, V7) 0.000000 3.000000 VOLUME( WH1, V8) 0.000000 8.000000 VOLUME( WH2, V1) 0.000000 1.000000 VOLUME( WH2, V2) 0.000000 2.000000 VOLUME( WH2, V3) 0.000000 2.000000 VOLUME( WH2, V4) 43.00000 0.000000 VOLUME( WH2, V5) 0.000000 5.000000 VOLUME( WH2, V6) 0.000000 2.000000 VOLUME( WH2, V7) 4.000000 0.000000 VOLUME( WH2, V8) 0.000000 0.000000 VOLUME( WH3, V1) 0.000000 4.000000 VOLUME( WH3, V2) 0.000000 1.000000 VOLUME( WH3, V3) 42.00000 0.000000VOLUME( WH3, V6) 0.000000 3.000000 VOLUME( WH3, V7) 0.000000 7.000000 VOLUME( WH3, V8) 0.000000 3.000000 VOLUME( WH4, V1) 0.000000 5.000000 VOLUME( WH4, V2) 0.000000 4.000000 VOLUME( WH4, V3) 0.000000 5.000000 VOLUME( WH4, V4) 0.000000 7.000000 VOLUME( WH4, V5) 0.000000 7.000000 VOLUME( WH4, V6) 14.00000 0.000000 VOLUME( WH4, V7) 0.000000 5.000000 VOLUME( WH4, V8) 38.00000 0.000000 VOLUME( WH5, V1) 41.00000 0.000000 VOLUME( WH5, V2) 0.000000 1.000000 VOLUME( WH5, V3) 0.000000 4.000000 VOLUME( WH5, V4) 0.000000 3.000000 VOLUME( WH5, V5) 0.000000 5.000000 VOLUME( WH5, V6) 0.000000 0.000000 VOLUME( WH5, V7) 0.000000 4.000000 VOLUME( WH5, V8) 0.000000 4.000000 VOLUME( WH6, V1) 0.000000 4.000000 VOLUME( WH6, V2) 0.000000 8.000000 VOLUME( WH6, V3) 0.000000 1.000000 VOLUME( WH6, V4) 0.000000 1.000000 VOLUME( WH6, V5) 0.000000 7.000000 VOLUME( WH6, V6) 32.00000 0.000000 VOLUME( WH6, V7) 0.000000 3.000000 VOLUME( WH6, V8) 0.000000 3.000000Row Slack or Surplus Dual Price1 473.0000 -1.0000002 19.00000 0.0000003 0.000000 0.0000004 9.000000 0.0000005 0.000000 0.0000006 41.00000 0.0000007 6.000000 0.0000008 39.00000 0.0000009 0.000000 1.00000010 0.000000 -2.00000011 0.000000 -3.00000012 0.000000 -1.00000013 0.000000 -2.00000014 0.000000 -2.0000003、一般整数规划问题：某服务部门各时段（每2h为一时段）需要的服务员人数见下表。