第二讲线性规划模型教材课件.ppt

《第二讲线性规划模型教材课件.ppt》由会员分享，可在线阅读，更多相关《第二讲线性规划模型教材课件.ppt（42页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第二讲第二讲 线性规划模型线性规划模型统计与应用数学系统计与应用数学系张耀峰张耀峰The The modelmodel of of linear programminglinear programming2/5/20231第第二二讲讲 线性规划模型线性规划模型1.1 优化的思想优化的思想1.2 什么是线性规划模型什么是线性规划模型1.3 如何使用如何使用Lingo软件求解软件求解 线性规划问题线性规划问题1.4 案案例解析例解析2/5/202321.1 优化的思想优化的思想2/5/20233烧水烧水 小明同学，烧一壶水要小明同学，烧一壶水要8分钟，灌开水要分钟，灌开水要1分钟，分钟，取牛奶和报

2、纸要取牛奶和报纸要5分钟（不能间断），整理书包要分钟（不能间断），整理书包要6分分钟（可间断），为了尽快做完这些事，怎样安排才能钟（可间断），为了尽快做完这些事，怎样安排才能使时间最少？最少需要几分钟？使时间最少？最少需要几分钟？例例1、如何安排早上的时间？、如何安排早上的时间？取牛奶和报纸取牛奶和报纸收拾书包收拾书包灌水灌水收拾书包收拾书包5891202/5/20234例例2、怎么排队才合理呢？、怎么排队才合理呢？码头上现在同时有码头上现在同时有3艘货船需要卸货，但是只艘货船需要卸货，但是只能一条一条地卸货，并且每艘船卸货所需的时间能一条一条地卸货，并且每艘船卸货所需的时间各不相同，分别为各

3、不相同，分别为4小时、小时、8小时和小时和1小时。按照怎小时。按照怎样的顺序卸货能使样的顺序卸货能使3艘货船等候的总时间最少呢？艘货船等候的总时间最少呢？2/5/20235方案卸货顺序船1的等候时间船2的等候时间船3的等候时间总的等候时间1船1船2船388+48+4+1332船1船3船288+1+48+1303船2船1船34+844+8+1294船2船3船14+1+844+1225船3船1船21+81+8+41236船3船2船11+4+81+41192/5/202361.2 什么是线性规划模型什么是线性规划模型2/5/20237例例3 3 运输问题运输问题2/5/20238解：设解：设A1,A

4、2调运到三个粮站的大米分别为调运到三个粮站的大米分别为x11，x12，x13，x21，x22，x23吨。吨。题设量可总到下表：题设量可总到下表：84库库存存量量x23x22x21A2542需要量需要量x13x12x11A1B3B2B1粮库粮库粮站粮站距离及运量距离及运量121224308242/5/20239结合存量限制和需量限制得数学模型结合存量限制和需量限制得数学模型：目标目标函数函数约束约束条件条件决策决策变量变量2/5/2023101.3如何使用如何使用Lingo软件求解线性规划问题软件求解线性规划问题2/5/202311程序编写程序编写model:min=12*x11+24*x12+

5、8*x13+30*x21+12*x22+24*x23;x11+x12+x134;x21+x22+x232;x12+x224;x13+x235;end2/5/202312运行结果运行结果 Global optimal solution found.Objective value:160.0000 Total solver iterations:5 Variable Value Reduced Cost X11 2.000000 0.000000 X12 0.000000 28.00000 X13 2.000000 0.000000 X21 0.000000 2.000000 X22 4.0000

6、00 0.000000 X23 3.000000 0.000000 Row Slack or Surplus Dual Price 1 160.0000 -1.000000 2 0.000000 16.00000 3 1.000000 0.000000 4 0.000000 -28.00000 5 0.000000 -12.00000 6 0.000000 -24.000002/5/202313 例例4 生产计划问题生产计划问题某工厂计划安排生产某工厂计划安排生产，两种产品，两种产品，已知每种单位产品的利润，生产单位产品所需设备台时已知每种单位产品的利润，生产单位产品所需设备台时及及A,B两种

7、原材料的消耗，现有原材料和设备台时的定额两种原材料的消耗，现有原材料和设备台时的定额如表所示，问：如表所示，问：）怎么安排生产使得工厂获利最大？）怎么安排生产使得工厂获利最大？）产品）产品的单位利润降低到的单位利润降低到1.8万元，要不要改变生产计万元，要不要改变生产计划，如果降低到划，如果降低到1万元呢？万元呢？）产品）产品的单位利润增大到的单位利润增大到5万，要不要改变生产计划万，要不要改变生产计划）如果产品）如果产品，的单位利润同时降低了的单位利润同时降低了1万元，要不万元，要不要改变生产计划？要改变生产计划？产品产品最大资源量设备128台时原材料A4016kg原材料B0412kg单位产

8、品利润232/5/2023142/5/202315程序编写程序编写model:title 生产计划问题生产计划问题;maxfmax=2*x1+3*x2;Ax1+2*x28;B4*x116;TIME4*x212;END2/5/202316运行结果运行结果 Model Title:生产计划问题生产计划问题 Variable Value Reduced Cost X1 4.000000 0.000000 X2 2.000000 0.000000 Row Slack or Surplus Dual Price MAXF 14.00000 1.000000 A 0.000000 1.500000 B 0

9、.000000 0.1250000 TIME 4.000000 0.000000 对问题对问题1，安排是生产产品，安排是生产产品4单位，产品单位，产品2单位，单位，最大盈利为最大盈利为14万元万元。2/5/202317 线性模型敏感性分析要使用敏感性分析要使用敏感性分析必须要在这里选择必须要在这里选择Prices&Ranges然后然后保存保存退出退出路径：路径：LINGOOptionsGeneral Solver(通用求解程序通用求解程序)选项卡选项卡2/5/202318要调出敏感性分析的结果，要调出敏感性分析的结果，必须必须先求解先求解后再后再在程序窗在程序窗口下口下点击点击LINGORan

10、ge，2/5/202319 Ranges in which the basis is unchanged:Objective Coefficient Ranges Current Allowable Allowable Variable Coefficient Increase Decrease X1 2.000000 INFINITY 0.5000000 X2 3.000000 1.000000 3.000000 Righthand Side Ranges Row Current Allowable Allowable RHS Increase Decrease A 8.000000 2.0

11、00000 4.000000 B 16.00000 16.00000 8.000000 TIME 12.00000 INFINITY 4.000000 当前变量系数允许增加量允许减少量2/5/202320 对问题对问题4，因为两个系数同时改变了，所以只，因为两个系数同时改变了，所以只有更有更 改程序的数据，重新运行得：不改变生产计改程序的数据，重新运行得：不改变生产计划，但是最大利润降低到划，但是最大利润降低到6万元万元.对问题对问题2，产品，产品的单位利润降低到的单位利润降低到1.8万元，在万元，在（1.5，）之间，所以不改变生产计划。如果降低）之间，所以不改变生产计划。如果降低到到1万元，

12、不在（万元，不在（1.5，）内，要改变生产计划。在）内，要改变生产计划。在程序中将目标函数的系数程序中将目标函数的系数“2”改为改为“1”，可得新，可得新的计划为安排是生产产品的计划为安排是生产产品2单位，产品单位，产品3单位，最单位，最大盈利为大盈利为11万元万元.对问题对问题3，要改变生产计划，更改程序得新，要改变生产计划，更改程序得新计划为生产产品计划为生产产品2单位，产品单位，产品3单位，最大盈利单位，最大盈利为为19万元万元.2/5/202321例例5 加工奶制品的生产计划加工奶制品的生产计划1桶牛奶 3公斤A1 12小时 8小时 4公斤A2 或获利24元/公斤 获利16元/公斤 5

13、0桶牛奶桶牛奶 时间时间480小时小时 至多加工至多加工100公斤公斤A1 制订生产计划，使每天获利最大制订生产计划，使每天获利最大 35元可买到元可买到1桶牛奶，买吗？若买，每天最多买多少桶牛奶，买吗？若买，每天最多买多少?可聘用临时工人，付出的工资最多是每小时几元可聘用临时工人，付出的工资最多是每小时几元?A1的获利增加到的获利增加到 30元元/公斤，应否改变生产计划？公斤，应否改变生产计划？每天：每天：2/5/2023221桶牛奶 3公斤A1 12小时 8小时 4公斤A2 或获利24元/公斤 获利16元/公斤 x1桶牛奶生产桶牛奶生产A1 x2桶牛奶生产桶牛奶生产A2 获利获利 243x

14、1 获利获利 164 x2 原料供应原料供应 劳动时间劳动时间 加工能力加工能力 决策变量决策变量 目标函数目标函数 每天获利每天获利约束条件约束条件非负约束非负约束 线性线性规划规划模型模型(LP)时间时间480小时小时 至多加工至多加工100公斤公斤A1 50桶牛奶桶牛奶 每天每天2/5/202323模型求解模型求解 OBJECTIVE FUNCTION VALUE 1)3360.000 VARIABLE VALUE REDUCED COST X1 20.000000 0.000000 X2 30.000000 0.000000 ROW SLACK OR SURPLUS DUAL PRIC

15、ES 2)0.000000 48.000000 3)0.000000 2.000000 4)40.000000 0.000000 NO.ITERATIONS=220桶牛奶生产桶牛奶生产A1,30桶生产桶生产A2，利润利润3360元。元。max=72*x1+64*x2;x1+x250;12*x1+8*x2480;3*x1100;2/5/202324模型求解模型求解 reduced cost值值表表示示当当该该非非基基变变量量增增加加一一个个单单位位时时（其其他他非非基基变变量量保保持持不不变变）目目标标函函数数减减少少的的量量(对对max型问题型问题)OBJECTIVE FUNCTION VAL

16、UE 1)3360.000 VARIABLE VALUE REDUCED COST X1 20.000000 0.000000 X2 30.000000 0.000000 ROW SLACK OR SURPLUS DUAL PRICES 2)0.000000 48.000000 3)0.000000 2.000000 4)40.000000 0.000000 NO.ITERATIONS=22/5/202325 OBJECTIVE FUNCTION VALUE 1)3360.000 VARIABLE VALUE REDUCED COST X1 20.000000 0.000000 X2 30.0

17、00000 0.000000 ROW SLACK OR SURPLUS DUAL PRICES 2)0.000000 48.000000 3)0.000000 2.000000 4)40.000000 0.000000原料无剩余原料无剩余时间无剩余时间无剩余加工能力剩余加工能力剩余40max 72x1+64x2st2）x1+x2503）12x1+8x24804）3x1100end三三种种资资源源结果解释结果解释 2/5/202326 OBJECTIVE FUNCTION VALUE 1)3360.000 VARIABLE VALUE REDUCED COST X1 20.000000 0.00

18、0000 X2 30.000000 0.000000 ROW SLACK OR SURPLUS DUAL PRICES 2)0.000000 48.000000 3)0.000000 2.000000 4)40.000000 0.000000结果解释结果解释 最优解下最优解下“资源资源”增增加加1单位时单位时“效益效益”的的增量增量 原料增原料增1单位单位,利润增利润增48 时间加时间加1单位单位,利润增利润增2 能力增减不影响利润能力增减不影响利润影子价格影子价格 35元可买到元可买到1桶牛奶，要买吗？桶牛奶，要买吗？35 sum(yuefen(i)|i#le#j:d);sum(yuefen

19、:x)=sum(yuefen:d);for(yuefen:xa);end 2/5/202333露天矿里铲位已分成矿石和岩石露天矿里铲位已分成矿石和岩石:平均铁含量不低于平均铁含量不低于25%的为矿石，否则为岩石。每个铲位的矿石、岩石数的为矿石，否则为岩石。每个铲位的矿石、岩石数量，以及矿石的平均铁含量（称为品位）都是已知的。量，以及矿石的平均铁含量（称为品位）都是已知的。每个铲位至多安置一台电铲，电铲平均装车时间每个铲位至多安置一台电铲，电铲平均装车时间5分钟分钟卡车在等待时所耗费的能量也是相当可观的，原则上在卡车在等待时所耗费的能量也是相当可观的，原则上在安排时不应发生卡车等待的情况。安排时

20、不应发生卡车等待的情况。例例7、露天矿生产的车辆安排露天矿生产的车辆安排(CUMCM-2003B)矿石卸点需要的铁含量要求都为矿石卸点需要的铁含量要求都为29.5%1%(品位限制）品位限制），搭配量在一个班次（，搭配量在一个班次（8小时）内满足品位限制即可。小时）内满足品位限制即可。卸点在一个班次内不变。卡车载重量为卸点在一个班次内不变。卡车载重量为154吨，平均时吨，平均时速速28km,平均卸车时间为平均卸车时间为3分钟。分钟。问题：出动几台电铲，分别在哪些铲位上；出动几辆问题：出动几台电铲，分别在哪些铲位上；出动几辆卡车，分别在哪些路线上各运输多少次卡车，分别在哪些路线上各运输多少次?2/

21、5/202334平面示意图2/5/202335问题数据问题数据 距离铲位1铲位2铲位3铲位4铲位5铲位6铲位7铲位8铲位9铲位10矿石漏5.265.194.214.002.952.742.461.900.641.27倒装1.900.991.901.131.272.251.482.043.093.51岩场5.895.615.614.563.513.652.462.461.060.57岩石漏0.641.761.271.832.742.604.213.725.056.10倒装4.423.863.723.162.252.810.781.621.270.50铲位1铲位2铲位3铲位4铲位5铲位6铲位7铲位

22、8铲位9铲位10矿石量095105100105110125105130135125岩石量125110135105115135105115135125铁含量30%28%29%32%31%33%32%31%33%31%2/5/202336问题分析问题分析 与典型的运输问题明显有以下不同：与典型的运输问题明显有以下不同：1.这是运输矿石与岩石两种物资的问题；这是运输矿石与岩石两种物资的问题；2.属于产量大于销量的不平衡运输问题；属于产量大于销量的不平衡运输问题；3.为了完成品位约束，矿石要搭配运输；为了完成品位约束，矿石要搭配运输；4.产地、销地均有单位时间的流量限制；产地、销地均有单位时间的流量限

23、制；5.运输车辆只有一种，每次满载运输，运输车辆只有一种，每次满载运输，154吨吨/车次；车次；6.铲位数多于铲车数意味着要最优的选择不多于铲位数多于铲车数意味着要最优的选择不多于7个个产地作为最后结果中的产地；产地作为最后结果中的产地；7.最后求出各条路线上的派出车辆数及安排。最后求出各条路线上的派出车辆数及安排。近似处理：近似处理：先求出产位、卸点每条线路上的运输量先求出产位、卸点每条线路上的运输量(MIP模型模型)然后求出各条路线上的派出车辆数及安排然后求出各条路线上的派出车辆数及安排2/5/202337模型假设模型假设v卡车在一个班次中不应发生等待或熄火后再启动卡车在一个班次中不应发生

24、等待或熄火后再启动的情况；的情况；v在铲位或卸点处由两条路线以上造成的冲突问题在铲位或卸点处由两条路线以上造成的冲突问题面前，我们认为只要平均时间能完成任务，就认面前，我们认为只要平均时间能完成任务，就认为不冲突。我们不排时地进行讨论；为不冲突。我们不排时地进行讨论；v空载与重载的速度都是空载与重载的速度都是28km/h，耗油相差很大；，耗油相差很大；v卡车可提前退出系统，等等。卡车可提前退出系统，等等。如理解为严格不等待，难以用数学规划模型来解如理解为严格不等待，难以用数学规划模型来解 个别参数队找到了可行解个别参数队找到了可行解（略）（略）2/5/202338符号符号vxij：从：从i铲位

25、到铲位到j号卸点的石料运量号卸点的石料运量（车）（车）单位：单位：吨；吨；vcij：从：从i号铲位到号铲位到j号卸点的距离号卸点的距离 公里；公里；vTij:从从i号铲位到号号铲位到号j卸点路线上运行一个周期平均时间卸点路线上运行一个周期平均时间 分；分；vAij：从号铲位到号卸点最多能同时运行的卡车数：从号铲位到号卸点最多能同时运行的卡车数 辆；辆；vBij：从号铲位到号卸点路线上一辆车最多可运行的次数：从号铲位到号卸点路线上一辆车最多可运行的次数 次；次；vpi：i号铲位的矿石铁含量号铲位的矿石铁含量 p=(30,28,29,32,31,33,32,31,33,31)%vqj:j号卸点任务

26、需求，号卸点任务需求，q=(1.2,1.3,1.3,1.9,1.3)*10000 吨吨vcki：i号铲位的铁矿石储量号铲位的铁矿石储量 万吨万吨vcyi：i号铲位的岩石储量号铲位的岩石储量 万吨万吨vfi:描述第描述第i号铲位是否使用的号铲位是否使用的0-1变量，取变量，取1为使用；为使用；0为关闭。为关闭。(近似近似)2/5/202339优化模型优化模型（1）道路能力）道路能力(卡车数卡车数)约束约束（2）电铲能力约束）电铲能力约束（3）卸点能力约束）卸点能力约束（4）铲位储量约束）铲位储量约束（5）产量任务约束）产量任务约束（6）铁含量约束）铁含量约束（7）电铲数量约束）电铲数量约束（8）

27、整数约束）整数约束.xij为非负整数为非负整数fi 为为0-1整数整数2/5/202340计算结果（计算结果（LINGOLINGO软件）软件）铲铲位位1 1铲铲位位2 2铲铲位位3 3铲铲位位4 4铲铲位位5 5铲铲位位6 6铲铲位位7 7铲铲位位8 8铲铲位位9 9铲铲位位1010矿矿漏漏131354541111倒倒42424343岩岩场场70701515岩漏岩漏81814343倒倒13132 27070铲铲位位1 1铲铲位位2 2铲铲位位3 3铲铲位位4 4铲铲位位5 5铲铲位位6 6铲铲位位7 7铲铲位位8 8铲铲位位9 9铲铲位位1010矿矿石漏石漏0.8671.8620.314倒倒场

28、场1.0771.162岩岩场场1.8920.326岩石漏岩石漏1.8411.229倒倒场场0.6840.11.489cumcm2003b1.lg4注注:LINGO8.0本来是可以得到最优解的，但有些本来是可以得到最优解的，但有些 LINGO8.0可能出现系统错误可能出现系统错误,可能是系统可能是系统BUG2/5/202341计算结果（派车）计算结果（派车）铲位1铲位2铲位3铲位4铲位5铲位6铲位7铲位8铲位9铲位10矿石漏1(29)倒场1(39)1(37)岩场1(37)岩石漏1(44)1(35)倒场1(47)结论：结论：铲位铲位1、2、3、4、8、9、10处各放置一台电铲。处各放置一台电铲。一共使用了一共使用了13辆卡车；总运量为辆卡车；总运量为85628.62吨公里；吨公里；岩石产量为岩石产量为32186吨；矿石产量为吨；矿石产量为38192吨。吨。此外：此外：6辆联合派车（方案略）辆联合派车（方案略）2/5/202342