实验-实验三详解.doc

《运筹学》课程实报告实验名称： 灵敏性报告与运输规划 任课教师： ：姓名： 学 号：完成日期： 一、实验目的 掌握在LINGO进行灵敏度分析的操作方法； 掌握运输规划的求解。 二、实验内容 实验1： 实验2： 实验3：略难 三、实验步骤及设计说明 实验1： 3.该报告说明： （1）运行 3 步找到全局最优解，目标函数值为 29000 ，变量值分别为x1= 100 ，x2= 30 ； （2）在“Slack or Surplus”（松弛或剩余）列，对于“=”不等式，右减左的差值为Slack，对于“=”不等式，左减右的差值为Surplus，从Row2和Row4的数值0，表示资源 1号生产线 和 劳动工时 在最优生产计划下已经用尽，从Row3的数值90，表示 2号生产线的资源使用剩余的资源还有90 ； （3）“Dual Price”是对偶价格（影子价格），从报告中可看出，生产线1、生产线2和劳动工时的影子价格分别是 50 、 0 、 150 。 ①若有2号生产线可以升级产量，升级成本是300，升级后的产量是200，问是否值得升级____不值得____? ②若通过加班来增加工人劳动力时间，那么付给工人的加班费每小时最多不能超过_150_______. 5.在模型窗口中，选择LINGO/Range选项，或按Ctrl+R，可得到该模型的灵敏度分析报告： 可以看到在最优基不改变的条件下， x1的目标函数系数c1有效范围是［150，+∞］， x2的目标函数系数c2有效范围是 [ 0，400 ] ； 对应1号生产线资源b1的有效范围是 [ 0,160 ] ； 对应2号生产线资源b2的有效范围是 [ 30，+∞ ] ； 对应劳动工时b3的有效范围是 [ 100,340 ] . 实验2： 实验代码： sets: col/1..3/:a; row/1..4/:b; aa(col,row):c,x; endsets data: a = 7,4,9; b = 3,6,5,6; c = 3,11,3,10, 1,9,2,8, 7,4,10,5; enddata min = @sum (aa:c*x); @for (col(i):@sum (row(j):x(i,j)) = a(i)); @for (row(j):@sum (col(i):x(i,j)) = b(j)); B1 B2 B3 B4 产量 A1 2 0 5 0 7 A2 1 0 0 3 4 A3 0 6 0 3 9 销量 3 6 5 6 85 实验3： 代码： max = 550*x1 + 600*x2 + 350*x3 + 400*x4 + 200*x5; 12*x1 + 20*x2 + 25*x4 + 15*x5 = 864; 10*x1 + 8*x2 + 16*x3 = 576; 20*x1 + 20*x2 + 20*x3 + 20*x4 + 20*x5 = 1152; x2 = 10; x2 = 100; x4 = 20; x4 = 150; 并结合模型结果报告和灵敏度分析报告回答下列问题： （1）写出正确的线性规划模型与生产计划。 解：设x(i)为一周内产品p(i)的生产量，i=1,2,3,4,5；由题得，各种资源的限额为： 研磨工时：9*6*2*8=864； 钻孔工时：6*6*2*8=576； 装配工时：24*6*8=1152； 又可得10≤x2≤100,20≤x4≤150； 所以原问题的线性规划模型如下： 目标函数：maxz=550*x1+600*x2+350*x3+400*x4+200*x5; 约束条件：12*x1+20*x2+25*x4+15*x5≤864; 10*x1+8*x2+16*x3≤576; 20*x1+20*x2+20*x3+20*x4+20*x5≤1152; x2≥10; x2≤100; x4≥20; x4≤150; x(i)≥0,i=1,2,3,4,5. 由lingo及联系实际解得：x1=14,x2=10,x3=14,x4=20,x5=0,maxz=26600 生产计划是生产产品p1,p2,p3,p4,p5分别为14,10,14,20,0件。 （2）这家工厂还有剩余的资源吗？若有，是哪种资源？有多少？ 答：由截图可以看出，该工厂还有13