2022年运筹学安全工程案例 .pdf

运筹学典型题型案例集安全工程专业第一章线性规划1 生产计划问题（ （摘自王治祯环境应用数学309 页） ）某企业为了搞好综合利用，用三种废品生产三种副产品，生产情况和利润见下表，求最佳利润。副产品废品A B C 最大日产量 /件甲10 5 3 400 乙6 10 2 300 丙4 5 4 200 利润5 8 2 解：设 ABC 三种产品的产量为X1X2X3Max Z =5X1+8X2+2X3 10X1+5X2+3 X3=400 6X1+10X2+2 X3=400 4X1+5X2+4X3=200 经过求解X1=34.23 ，X2=8.19 X3=5.37 最大利润为274.41 2 投资问题解：用 Xij 表示第 i 年初（ i=1，2，3）给项目 j（A，B，C，D）的投资金额。第一年资金量：30 万，可投项目：A、B；故： X1A+X1B =30。第二年资金量：1.2*X1A ，可投项目： A、C；故： X2A+X2C =1.2*X1A 。第三年资金量：1.2*X2A+1.5*X1B，可投项目：A、B、D；故： X3A+X3B+X3D=1.2*X2A+1.5*X1B。其它条件： X1B =20；X2C =15； X3D =10。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 21 页 - - - - - - - - - 目标：第三年底收益最大。因投资X3B 在第 3 年底不能收回，故无收益。则目标函数为：f(x)=0.2* （X1A+ X2A + X3A）+0.5*X1B+0.6* X2C+0.4* X3D LINGO Model如下：max =0.2*(X1A+ X2A + X3A)+0.5*X1B+0.6* X2C+0.4* X3D; X1A+X1B=30; X2A+X2C=1.2*X1A; X3A+X3B+X3D=4; X8+ X9=4; X8+ X9+ X10=6; X9+ X10+ X11=6; X10+ X11+ X12=8;X11+ X12+ X13=8; X12+ X13+ X14=6; X13+ X14+ X15=6; X14+ X15+ X16=4; X15+ X16+ X17=4;X16+ X17+ X18=6; X17+ X18+ X19=6; X18+ X19+ X20=8; X19+ X20+ X21=8; 运行LINGO 软件得到问题的最优解（只列出非零变量）：最优目标函数值 =32.00000 X8=4.000000 X10=4.000000 X11=2.000000 X12=2.000000 X13=4.000000 X15=2.000000 X16=2.000000 X17=4.000000 X19=2.000000 X20=6.000000 根据运行结果，最优时间表确定如下，此时最少人数为32人时间8：00 9：00 10：00 11：00 12：00 13：00 14：00 15：00 16：00 17：00 18：00 19：00 20：00 21：00 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 21 页 - - - - - - - - - 开始工作人数4 0 4 2 2 4 0 2 2 4 0 2 6 0 需要人数4 4 6 6 8 8 6 6 4 4 6 6 8 8 （2）没有志愿者愿意在 12：00和18：00开始工作，即增加约束条件：X12=0 ； X18=0。LINGO Model 如下：min= X8 +X9+ X10+ X11+ X12+ X13+ X14+ X15+ X16+ X17+ X18+ X19+ X20+ X21;X8=4; X8+ X9=4; X8+ X9+ X10=6; X9+ X10+ X11=6; X10+ X11+ X12=8;X11+ X12+ X13=8; X12+ X13+ X14=6; X13+ X14+ X15=6; X14+ X15+ X16=4; X15+ X16+ X17=4;X16+ X17+ X18=6; X17+ X18+ X19=6; X18+ X19+ X20=8; X19+ X20+ X21=8; X12=0; X18=0;运行LINGO 软件得到问题的最优解（只列出非零变量）：最优目标函数值 =32.00000 X8=4.000000 X10=6.000000 X11=2.000000 X13=6.000000 X16=4.000000 X17=2.000000 X19=4.000000 X20=4.000000 根据运行结果，最优时间表确定如下，此时最少人数为32人时间8：00 9：00 10：00 11：00 12：00 13：00 14：00 15：00 16：00 17：00 18：00 19：00 20：00 21：00 开始工作人数4 0 6 2 0 6 0 0 4 2 0 4 4 0 需要人数4 4 6 6 8 8 6 6 4 4 6 6 8 8 第二章灵敏度分析名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 21 页 - - - - - - - - - 解： （1）设每天生产 A1产品用奶 x1桶，生产 A2产品用奶 x2桶。则LINGO Model如下：max= 24*3*x1+ 16*4*x2; x1+ x2 =50; 12*x1+ 8*x2 =480; 3*x1 =100;运行结果如下： Global optimal solution found. Objective value: 3360.000 Total solver iterations: 2 Variable Value Reduced Cost X1 20.00000 0.000000 X2 30.00000 0.000000 Row Slack or Surplus Dual Price 1 3360.000 1.000000 2 0.000000 48.00000 3 0.000000 2.000000 4 40.00000 0.000000 灵敏度分析结果如下：Ranges in which the basis is unchanged: 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 21 页 - - - - - - - - - Objective Coefficient Ranges Current Allowable Allowable Variable Coefficient Increase Decrease X1 72.00000 24.00000 8.000000 X2 64.00000 8.000000 16.00000 Righthand Side Ranges Row Current Allowable Allowable RHS Increase Decrease 2 50.00000 10.00000 6.666667 3 480.0000 53.33333 80.00000 4 100.0000 INFINITY 40.00000 生产计划：每天用 20桶奶生产 A1产品，用 30桶奶生产 A2产品获利最大，每天可获利 3360元。附加问题：由影子价格可知，原料增加1单位，利润增长 48元，成本为 35元，所以可以买。由灵敏度分析结果，每天最多再购买10桶牛奶。由影子价格可知，时间增加1单位，利润增长 2元，所以聘用临时工人的工资最多 2元/小时。由灵敏度分析可知， x1系数范围是 (72-8,72+24) ， 当A1产品获利增加到 30元/kg时，即 x1系数为 30*3=9072+24 ，在允许范围内，所以不应改变生产计划。2 课本81页2.12 林敏度分析第三章运输问题1 产销平衡的运输问题（摘自王治祯环境应用数学330 页）某城市有三个工厂,每个工厂生产都产出一定量的剩余物(通称为污染物),本着化害为宝的精神 ,需将各厂的废物分别输送到本市内其他单位搞综合利用,一直每厂的剩余物和各厂的需要量及运价表,试用表上作业法求满足现有条件的运费最少的分配方案. 需求地废品产出地A B C 产量甲10 30 27 30 乙16 19 22 50 丙20 15 10 80 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 21 页 - - - - - - - - - 需要量35 40 85 160 解:(中间过程略 )最优运送方案表为需求地废品产出地A B C 产量甲30 30 乙5 40 5 50 丙80 80 需要量35 40 85 160 此时总运费为2050 2 产销不平衡的运输问题解： （1）设 xij(i,j=1,2,3) 为产地 i 运往客户 j 的运量，列表如下（运量、运费） ：客户 1 客户 2 客户 3 发量产地 1 X11、10 X12、4 X13 、12 3000 产地 2 X21、8 X22、10 X23 、3 4000 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 21 页 - - - - - - - - - 产地 3（虚产地）X31、- X32、- X33 、- 1500 需求量2000 1500 5000 建立数学模型如下：min =10*x11+4*x12+12*x13+8*x21+10*x22+3*x23;目标x11+x12+x13=3000;约束条件x21+x22+x23=4000;约束条件x31+x32+x33=2000;约束条件x12+x22+x32=1500;约束条件x13+x23+x33=5000;约束条件运行结果如下：（只列出部分结果） Objective value: 33000.00 Variable Value X11 1500.000 X12 1500.000 X13 0.000000 X21 0.000000 X22 0.000000 X23 4000.000 X31 500.0000 X32 0.000000 X33 1000.000 运输方案：产地 1分别给客户 1、2发货1500单位；产地2给客户 3发货4000单位；产地3给客户 1发货500单位，给客户 3发货1000单位。最低运费为 33000。3 课本 104 页 3.9 飞行安全问题转化的运输问题4 人员分配问题某检测站在一定时间内要化验一批水样,共分析 16 个项目 ,其中重金属A4 项,有机物 B5项,物理分析C3 项,其他分析D4 项,在分析期间 ,一部分人外出采样,站立仅仅有4 人能够承担名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 21 页 - - - - - - - - - 化验任务 ,根据每人世纪工作情况,化验员甲乙丙丁能完成435项任务，每人完成不同任务所需要时间见下表所示，求需时最少的分配方案甲乙丙丁注释：划横线为不能胜任此项工作解：此例题为运输问题，经过表上作业法，结果为甲乙丙丁此时用时最少为第四章整数规划与分配问题（摘自北京工业大学薛毅编写的数学建模实验）1 监控摄像头的最优安装问题（0-1 整数规划）名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 21 页 - - - - - - - - - 2 0-1 整数规划课本 125 页例题 63 指派问题1 4 指派问题2 （摘自王治祯环境应用数学334 页）利用 5 种不同质量浓度的有机废水灌溉5 块土质相同的草地,各质量浓度的污水灌溉草地收获量见下表 ,求最优灌溉方案和最优分配下的最大生产量(单位公斤 ) 1 2 3 4 5 1 73 82 115 43 67 2 120 78 88 80 50 3 65 90 75 125 40 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 21 页 - - - - - - - - - 4 87 130 45 65 100 5 90 45 65 75 80 提示 : 将各行都剪去该行最大数,再列效率表 ,则求最大生产量问题转化为求最小用时问题. 1 2 3 4 5 1 42 33 0 72 48 2 0 42 32 40 70 3 60 35 50 0 85 4 43 0 85 65 30 5 0 45 25 15 10 以下仍按照表上作业指派问题的求解方法来做.最优指派方案为1 2 3 4 5 1 1 2 1 3 1 4 1 5 1 此时最佳收获量为120+130+115+125+80=570 第八章动态规划1 最短路径问题电器安全可靠性问题（课本216页8.9） 招聘问题（课本 216页8.12）第九章决策分析1 风险型决策的决策树法（课本317 页 11.511.611.8 类似）某厂因生产需要，考虑是否自行研制一个新的安全装置，首先，决定这个研制项目是否需要评审，如果需要评审，则需要评审费5000 元，不评审，则可省去这笔评审费用，如进行评审，通过概率为0.8，不通过概率为0.2。研制过程可以采取独立完成和外厂协作形式，如果研制成功，有6 万元的收益，若采用本厂独立完成形式，则研制费用为2.5 万元，成功概率为 0.7，失败概率为0.3；若采用外厂协作性质，成功概率为0.99，失败概率为0.01，但是需要支付4 万元的研制费用。试问该厂决策者应该如何决策？（用决策树法求解）摘自张景林主编安全系统工程149 页答案：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 21 页 - - - - - - - - - 故决策者应该参加项目评审，并且评审通过后应该采用写作完成的策略，此时预期的收益最大为 1.082 万元。2 不确定性决策317 页 11.1 河北科技大学环境科学专业第一章线性规划1 生产计划问题（摘自王治祯环境应用数学301 页）某厂生产甲乙两种产品，设生产1t 单位甲产品，排放污水2 个单位，排渣2 个单位，生产1t 单位乙产品，排放污水4 个单位，排渣1 个单位。甲产品单价位500 元，乙产品单价为400 元，按照环保部门要求， 为了使生产中每天排放的废水和废渣不超过20 个和 10 个单位，问甲乙两种产品各应生产多少，才能使产值最大？要求采用图解法和单纯性法求解，并与lingo 求解结果作比较？解：图解法略单纯性法略Lingo 程序max =5X1+4X2 2X1+4X2=20 2X1+X2=2.99 X1=1 X2=1 X3=1 计算结果X1=0 ，X2=0.859 X3=0.542 所需要的资金最少为132.977 万元3 城市安全供水问题（摘自王治祯环境应用数学312 页）一个 40 万人口的城市，由于工业及人口的不断发展，面临缺水问题，初步估计全市缺水约5000t/d，水源及有关资料见下表，同时要求水质中硫化物，氯化物和可溶性固体分别不超过 650， 200，100mg/L ，是决定最低成本的供水策略。A 河净化水B 河地下水可溶性固体mg/L 100 150 50 20 氯化物 mg/L 700 500 800 600 硫化物 mg/L 300 100 250 40 价格元/t 0.5 0.4 0.7 0.6 解：设 X1X2X3 X4位分别来自4 个水源地的最低成本下的供水量(单位万 t)。Min Z =0.5X1+0.4X2+0.7X3 +0.6X4 700X1+500X2+800X3 +600X4=650 （X1+X2+X3 +X4 ）300X1+100X2+250X3 +40X4=200 （X1+X2+X3 +X4 ）100X1+150X2+50X3 +20X4=0.5 经过求解X1=0.29，X2=0.14 X3=0, X4=0.11 最低成本为2.67 万元运筹学软件LINGO 的使用一LINGO 简介LINGO 是Linear Interactive and General Optimizer的缩写，即“ 交互式的线性和通用优化求解器” ，可以用于求解非线性规划，也可以用于一些线性和非线性方程组的求解等。其特色在于可以允许决策变量是整数(即整数规划，包括0-1 整数规划 )，方便灵活，而且执行速度非常快。二LINGO 在各类问题中的应用实例1 课本 27 页例题 5的lingo程序max 2x1+3x2 ST 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 21 页 - - - - - - - - - 2x1+2x2 =12 4x1 =16 5x2 =15 End lingo求解结果 Global optimal solution found. Objective value: 15.00000 Infeasibilities: 0.000000 Total solver iterations: 1 Model Class: LP Total variables: 2 Nonlinear variables: 0 Integer variables: 0 Total constraints: 4 Nonlinear constraints: 0 Total nonzeros: 6 Nonlinear nonzeros: 0 Variable Value Reduced Cost X1 3.000000 0.000000 X2 3.000000 0.000000 Row Slack or Surplus Dual Price 1 15.00000 1.000000 2 0.000000 1.000000 3 4.000000 0.000000 4 0.000000 0.2000000 2 数学建模基础216 页0-1 整数规划的 lingo程序max 3x1-2x2+5x3 ST x1+2x2-x3 =2 x1+4x2+x3 =4 x1+x3 =3 4x2+x3=1 x1+x7 =1 x3+x4+x5+x6=1 x6+x8 =1 end INT8 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 21 页 - - - - - - - - - lingo求解结果 Global optimal solution found. Objective value: 2.000000 Objective bound: 2.000000 Infeasibilities: 0.000000 Extended solver steps: 0 Total solver iterations: 0 Model Class: PILP Total variables: 8 Nonlinear variables: 0 Integer variables: 8 Total constraints: 5 Nonlinear constraints: 0 Total nonzeros: 18 Nonlinear nonzeros: 0 Variable Value Reduced Cost X1 1.000000 1.000000 X2 0.000000 1.000000 X3 0.000000 1.000000 X4 0.000000 1.000000 X5 0.000000 1.000000 X6 1.000000 1.000000 X7 0.000000 1.000000 X8 0.000000 1.000000 Row Slack or Surplus Dual Price 1 2.000000 -1.000000 2 0.000000 0.000000 3 0.000000 0.000000 4 0.000000 0.000000 5 0.000000 0.000000 4 数学建模基础215 页线性规划的lingo程序Min -3x1+x2+x3 ST x1-2x2 +x3=3 -2x1+x3=1 End 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 21 页 - - - - - - - - - lingo求解结果 Global optimal solution found. Objective value: -2.000000 Infeasibilities: 0.000000 Total solver iterations: 0 Model Class: LP Total variables: 3 Nonlinear variables: 0 Integer variables: 0 Total constraints: 4 Nonlinear constraints: 0 Total nonzeros: 11 Nonlinear nonzeros: 0 Variable Value Reduced Cost X1 4.000000 0.000000 X2 1.000000 0.000000 X3 9.000000 0.000000 Row Slack or Surplus Dual Price 1 -2.000000 -1.000000 2 0.000000 0.3333333 3 0.000000 -0.3333333 4 0.000000 -0.6666667 5 人力资源分配问题某中型商场经过统计对销售人员的需求如下周一至周五20,16,13,16，19,14,12 人，销售人员每周连续工作5 天，休息 2 天，试求应该如何安排销售人员的休息时间，使销售人员总数为最少？解：设某天开始休息的人数为xi，目标函数min= Xi St x2+x3+x4+x5+x6=20 x3+x4+x5+x6+x7 =16 x1+x4+x5+x6+x7 =13 x1+x2+x5+x6+x7=16 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 17 页，共 21 页 - - - - - - - - - x1+x2+x3+x6+x7=19 x1+x2+x3+x4+x7=14 x1+x2+x3+x4+x5=12 Xi=0, XiZ Lingo求解结果 Global optimal solution found. Objective value: 22.00000 Objective bound: 22.00000 Infeasibilities: 0.000000 Extended solver steps: 0 Total solver iterations: 6 Model Class: PILP Total variables: 7 Nonlinear variables: 0 Integer variables: 7 Total constraints: 8 Nonlinear constraints: 0 Total nonzeros: 42 Nonlinear nonzeros: 0 Variable Value Reduced Cost X1 0.000000 1.000000 X2 6.000000 1.000000 X3 3.000000 1.000000 X4 3.000000 1.000000 X5 0.000000 1.000000 X6 8.000000 1.000000 X7 2.000000 1.000000 Row Slack or Surplus Dual Price 1 22.00000 -1.000000 2 0.000000 0.000000 3 0.000000 0.000000 4 0.000000 0.000000 5 0.000000 0.000000 6 0.000000 0.000000 7 0.000000 0.000000 8 0.000000 0.000000 6 指派问题例 2 （选自运 P94 习题 2.4；整数规则）有四个工人，要分别指派他们完成四项不同的名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 18 页，共 21 页 - - - - - - - - - 工作， 每个人做各项工作所消耗的时间如表。问应该如何指派， 才能使总的消耗时间为最小？工人任务A B C D 甲15 18 21 24 乙19 23 22 18 丙26 17 16 19 丁19 21 23 17 这是一道典型的整数规则问题。我们记派第I 去做工作记为Xij 注意到每人只能做一项工作。每项工作一人做。我们得到目标函数和约束条件：min 15x11+19x21+26x31+19x41+18x12+23x22+17x32+21x42+24x13+22x23+16x33+23x43+24x14+18x24+19x34+17x44 ST 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 int 16 运算结果Global optimal solution found. Objective value: 70.00000 Objective bound: 70.00000 Infeasibilities: 0.000000 Extended solver steps: 0 Total solver iterations: 0 Model Class: PILP Total variables: 16 Nonlinear variables: 0 Integer variables: 16 Total constraints: 9 Nonlinear constraints: 0 Total nonzeros: 48 Nonlinear nonzeros: 0 Variable Value Reduced Cost 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 19 页，共 21 页 - - - - - - - - - X11 0.000000 15.00000 X21 1.000000 19.00000 X31 0.000000 26.00000 X41 0.000000 19.00000 X12 1.000000 18.00000 X22 0.000000 23.00000 X32 0.000000 17.00000 X42 0.000000 21.00000 X13 0.000000 24.00000 X23 0.000000 22.00000 X33 1.000000 16.00000 X43 0.000000 23.00000 X14 0.000000 24.00000 X24 0.000000 18.00000 X34 0.000000 19.00000 X44 1.000000 17.00000 Row Slack or Surplus Dual Price 1 70.00000 -1.000000 2 0.000000 0.000000 3 0.000000 0.000000 4 0.000000 0.000000 5 0.000000 0.000000 6 0.000000 0.000000 7 0.000000 0.000000 8 0.000000 0.000000 9 0.000000 0.000000 三 实践环节作业占总成绩 12% 。在实践成绩环节，学生可以根据所学的运筹学知识，结合所在专业寻找问题， 自己建模，通过计算机软件求解，通过对计算结果进行数据分析，结合实际问题给出可行性建议，最后以课程读书报告或者实践报告形式的形式以电子版提交。提交时间 15 周以前。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 20 页，共 21 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 21 页，共 21 页 - - - - - - - - -