数学规划建模幻灯片.ppt
《数学规划建模幻灯片.ppt》由会员分享,可在线阅读,更多相关《数学规划建模幻灯片.ppt(115页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、数学规划建模数学规划建模第1页,共115页,编辑于2022年,星期六内容提要内容提要1什么是数学规划什么是数学规划2连续性线性规划连续性线性规划3整数线性规划整数线性规划4非线性规划非线性规划5多目标规划多目标规划6目标规划目标规划第2页,共115页,编辑于2022年,星期六最优化问题的最优化问题的数学模型的一般形式数学模型的一般形式为:为:(1)(2)三个要素:三个要素:决策变量决策变量decisionvariable,目标函数目标函数objectivefunction,约束条件约束条件constraints。什么是数学规划?什么是数学规划?第3页,共115页,编辑于2022年,星期六(2)
2、所确定的)所确定的x的范围称为的范围称为可行域可行域feasibleregion,满,满足(足(2)的解)的解x称为称为可行解可行解feasiblesolution,同时满足,同时满足(1)()(2)的解)的解x称为称为最优解最优解Optimalsolution,整个可,整个可行域上的最优解称为行域上的最优解称为全局最优解全局最优解globaloptimalsolution,可行域中某个领域上的最优解称为,可行域中某个领域上的最优解称为局部最优解局部最优解localoptimalsolution。最优解所对应的目标函数值称为。最优解所对应的目标函数值称为最优值最优值optimum。什么是数学规
3、划?什么是数学规划?第4页,共115页,编辑于2022年,星期六(一)(一)按有无约束条件按有无约束条件(2)可分为:)可分为:1.无约束优化无约束优化unconstrainedoptimization。2.约束优化约束优化constrainedoptimization。大部分实际问题都是约束优化问题。大部分实际问题都是约束优化问题。优化模型的分类优化模型的分类什么是数学规划?什么是数学规划?第5页,共115页,编辑于2022年,星期六(二)(二)按决策变量取值是否连续按决策变量取值是否连续可分为:可分为:1.数学规划或连续优化。数学规划或连续优化。可继续划分为可继续划分为线性规划线性规划(L
4、P)Linearprogramming和和非非线性规划线性规划(NLP)Nonlinearprogramming。在非线性规划。在非线性规划中有一种规划叫做中有一种规划叫做二次规划二次规划(QP)Quadraticprogramming,目标为二次函数,约束为线性函数。,目标为二次函数,约束为线性函数。2.离散优化或组合优化。离散优化或组合优化。包含:包含:整数规划整数规划(IP)Integerprogramming,整数规,整数规划中又包含很重要的一类规划:划中又包含很重要的一类规划:0-1(整数)规划(整数)规划Zero-oneprogramming,这类规划问题的决策变量只取,这类规划问
5、题的决策变量只取0或者或者1。什么是数学规划?什么是数学规划?第6页,共115页,编辑于2022年,星期六(三)(三)按目标的多少按目标的多少可分为:可分为:1.单目标规划。单目标规划。2.多目标规划。多目标规划。(四)(四)按模型中参数和变量是否具有不确定性按模型中参数和变量是否具有不确定性可分为:可分为:1.确定性规划。确定性规划。2.不确定性规划。不确定性规划。(五)(五)按问题求解的特性按问题求解的特性可分为:可分为:1.目标规划。目标规划。2.动态规划。动态规划。3.多层规划。多层规划。4.网络优化。网络优化。5.等等。等等。什么是数学规划?什么是数学规划?第7页,共115页,编辑于
6、2022年,星期六LINGO软件和软件和MATLAB软件。软件。求解优化问题常用的软件求解优化问题常用的软件什么是数学规划?什么是数学规划?第8页,共115页,编辑于2022年,星期六线性规划的一般形式:线性规划的一般形式:连续性线性规划连续性线性规划第9页,共115页,编辑于2022年,星期六一般线性规划问题都可以通过引入非负的松弛变量一般线性规划问题都可以通过引入非负的松弛变量slackvariable与非负的剩余变量与非负的剩余变量surplusvariable的方法化为的方法化为标准形标准形式式(约束全是等约束)。(约束全是等约束)。线性规划问题的可行域线性规划问题的可行域feasib
7、leregion是一个凸集是一个凸集convexset(任意两点的连线上的点都在区域内部,可以看(任意两点的连线上的点都在区域内部,可以看作是没有凹坑的凸多面体),所以最优解作是没有凹坑的凸多面体),所以最优解Optimalsolution/point在凸多面体的某个顶点上达到在凸多面体的某个顶点上达到求解方法:单纯形算法求解方法:单纯形算法simplexmethod。连续性线性规划连续性线性规划 第10页,共115页,编辑于2022年,星期六连续性线性规划连续性线性规划 例例1 运输问题运输问题第11页,共115页,编辑于2022年,星期六解解设设A1,A2调运到三个粮站的大米分别为调运到三
8、个粮站的大米分别为x1,x2,x3,x4,x5,x6吨。吨。题设量可总到下表:题设量可总到下表:连续性线性规划连续性线性规划 第12页,共115页,编辑于2022年,星期六结合存量限制和需量限制得数学模型结合存量限制和需量限制得数学模型:连续性线性规划连续性线性规划 第13页,共115页,编辑于2022年,星期六程序编写程序编写1model:min=12*x1+24*x2+8*x3+30*x4+12*x5+24*x6;x1+x2+x34;x4+x5+x62;x2+x54;x3+x65;end提示:课件中的程序提示:课件中的程序请先粘贴在记事本中请先粘贴在记事本中再再转帖于转帖于lingo软件中
9、软件中连续性线性规划连续性线性规划 第14页,共115页,编辑于2022年,星期六运行结果运行结果 Global optimal solution found.Objective value:160.0000 Total solver iterations:0 Variable Value Reduced Cost X1 2.000000 0.000000 X2 0.000000 28.00000 X3 2.000000 0.000000 X4 0.000000 2.000000 X5 4.000000 0.000000 X6 3.000000 0.000000 Row Slack or Su
10、rplus 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.00000连续性线性规划连续性线性规划 第15页,共115页,编辑于2022年,星期六模型修改为模型修改为 连续性线性规划连续性线性规划 第16页,共115页,编辑于2022年,星期六程序编写程序编写2MODEL:TITLE 调调运大米的运运大米的运输问题输问题程序程序2;!定义集合段定义集合段;SETS:HANG/1.5/:B
11、;!定义矩阵的行定义矩阵的行;LIE/1.6/:C,X;!定义矩阵的列以及变量定义矩阵的列以及变量;XISHU(HANG,LIE):A;!定义约束的系数矩阵定义约束的系数矩阵;ENDSETSDATA:A=1 1 1 0 0 0!系数矩阵赋值系数矩阵赋值;0 0 0 1 1 1 -1 0 0-1 0 0 0-1 0 0-1 0 0 0-1 0 0-1;连续性线性规划连续性线性规划 第17页,共115页,编辑于2022年,星期六C=12 24 8 30 12 24;!目标函数的系数;B=4 8-2-4-5;!约束的右端项;ENDDATA!标准形式的目标函数的矩阵形式;MIN=SUM(LIE:C*X
12、);FOR(HANG(I):SUM(LIE(J):A(I,J)*X(J)B(I);END连续性线性规划连续性线性规划 第18页,共115页,编辑于2022年,星期六84库库存存量量x23x22x21A2542需要量需要量x13x12x11A1B3B2B1粮库粮库粮站粮站距离及运量距离及运量12122430824变量更换为:变量更换为:连续性线性规划连续性线性规划 第19页,共115页,编辑于2022年,星期六模型模型:连续性线性规划连续性线性规划 第20页,共115页,编辑于2022年,星期六程序编写程序编写3MODEL:TITLE 调运大米的运输问题程序调运大米的运输问题程序3;!定义集合段
13、定义集合段;SETS:LIANGKU/1.2/:A;!定义粮库的集合定义粮库的集合;LIANGZHAN/1.3/:B;!定义粮站的集合定义粮站的集合;YULIANG(LIANGKU,LIANGZHAN):X,C;!定义运量和距离定义运量和距离;ENDSETSDATA:!粮库到粮站的距离粮库到粮站的距离;C=12 24 8 30 12 24;连续性线性规划连续性线性规划 第21页,共115页,编辑于2022年,星期六!粮粮库库的限量的限量;A=4 8;!粮站的限量粮站的限量;B=2 4 5;ENDDATAOBJMIN=SUM(YULIANG:C*X);!粮粮库库上限的上限的约约束束;FOR(LI
14、ANGKU(I):LK SUM(LIANGZHAN(J):X(I,J)B(J);END 连续性线性规划连续性线性规划 第22页,共115页,编辑于2022年,星期六程序的调试程序的调试1.直接点击运行,如果出错会弹出错误提示,根据提直接点击运行,如果出错会弹出错误提示,根据提示做相应的修改;示做相应的修改;2.可以用可以用“!”把约束变成说明语句,而把这条语句把约束变成说明语句,而把这条语句屏蔽掉,缩小寻找出错的范围;屏蔽掉,缩小寻找出错的范围;3.可以边写程序边运行,保证每行书写都是正确的程序;可以边写程序边运行,保证每行书写都是正确的程序;连续性线性规划连续性线性规划 第23页,共115页
15、,编辑于2022年,星期六例例2 阶段生产问题阶段生产问题(P31)某公司生产某产品某公司生产某产品,最大生产能力为最大生产能力为10000单位单位,每单每单位存储费位存储费2元元,预定的销售量与单位成本如下预定的销售量与单位成本如下:月份单位成本(元)销售量1234 70 6000 72 7000 80 12000 76 6000求一生产计划求一生产计划,使使1)满足需求满足需求;2)不超过生产能力不超过生产能力;3)成本成本(生产成本与存储费之和生产成本与存储费之和)最低最低.连续性线性规划连续性线性规划 第24页,共115页,编辑于2022年,星期六解解假定假定1月初无库存月初无库存,4
16、月底买完月底买完,当月生产的不库当月生产的不库存存,库存量无限制库存量无限制.连续性线性规划连续性线性规划 第25页,共115页,编辑于2022年,星期六model:title 生产计划程序生产计划程序1;Sets:yuefen/1.4/:c,x,e,d;endsetsdata:c=70 71 80 76;d=6000 7000 12000 6000;e=2 2 2 2;a=10000;enddatamin=sum(yuefen:c*x)+sum(yuefen(j)|j#lt#4:sum(yuefen(i)|i#le#j:x-d)*e(j+1);for(yuefen(j)|j#lt#4:sum
17、(yuefen(i)|i#le#j:x)sum(yuefen(i)|i#le#j:d);sum(yuefen:x)=sum(yuefen:d);for(yuefen:xa);end 连续性线性规划连续性线性规划 第26页,共115页,编辑于2022年,星期六连续性线性规划连续性线性规划 第27页,共115页,编辑于2022年,星期六Model:Title 生生产计产计划程序划程序2;Sets:yuefen/1.4/:c,x,e,d,s;endsetsdata:c=70 71 80 76;d=6000 7000 12000 6000;e=2 2 2 2;a=10000;enddatamin=su
18、m(yuefen:c*x+e*s);for(yuefen(i)|i#lt#4:s(i+1)=s(i)+x(i)-d(i);s(4)+x(4)-d(4)=0;s(1)=0;for(yuefen:xa);End连续性线性规划连续性线性规划 第28页,共115页,编辑于2022年,星期六月份单位成本(元)销售量1234 70 6000 72 7000 80 12000 76 6000连续性线性规划连续性线性规划 第29页,共115页,编辑于2022年,星期六76827676-80-7472-747270生产月10000100001000010000产量产量600041200070006000销量销量
19、4321321需求月费用cij连续性线性规划连续性线性规划 第30页,共115页,编辑于2022年,星期六连续性线性规划连续性线性规划 建立模型如下:建立模型如下:第31页,共115页,编辑于2022年,星期六model:title 生产计划程序生产计划程序3;sets:yuefen/1.4/:a,d,xx;!定义上三角矩阵定义上三角矩阵;link(yuefen,yuefen)|&2#ge#&1:c,x;endsetsdata:c=70 72 74 76 71 73 75 80 82 76;d=6000 7000 12000 6000;a=10000 10000 10000 10000;end
20、datamin=sum(link:c*x);for(yuefen(i):sum(yuefen(j)|j#ge#i:x(i,j)d(j););!得到每个月的生产量得到每个月的生产量;for(yuefen(i):xx=sum(yuefen(j)|j#ge#i:x(i,j);End连续性线性规划连续性线性规划 第32页,共115页,编辑于2022年,星期六Model Title:生生产计产计划程序划程序1 Variable Value Reduced Cost A 10000.00 0.000000 C(1)70.00000 0.000000 C(2)71.00000 0.000000 C(3)80
21、.00000 0.000000 C(4)76.00000 0.000000 X(1)10000.00 0.000000 X(2)10000.00 0.000000 X(3)5000.000 0.000000 X(4)6000.000 0.000000 E(1)2.000000 0.000000 E(2)2.000000 0.000000 E(3)2.000000 0.000000 E(4)2.000000 0.000000 D(1)6000.000 0.000000 D(2)7000.000 0.000000 D(3)12000.00 0.000000 D(4)6000.000 0.00000
22、0 连续性线性规划连续性线性规划 第33页,共115页,编辑于2022年,星期六连续投资连续投资10万元万元A:从第:从第1年年到第到第4年每年初要投资,次年末回收本年每年初要投资,次年末回收本利利1.15B:第第3年初投资,到第年初投资,到第5年末回收年末回收1.25,最大投资,最大投资4万万元元C:第第2年初投资,到第年初投资,到第5年末回收年末回收1.40,最大投资,最大投资3万元万元D:每年初投资,每年末回收每年初投资,每年末回收1.11。求:求:5年末总资本最大。年末总资本最大。练习一下吧练习一下吧 连续投资连续投资第34页,共115页,编辑于2022年,星期六 例例3 生产计划问题
23、生产计划问题某工厂计划安排生产某工厂计划安排生产,两种产品,已知每种单位产两种产品,已知每种单位产品的利润,生产单位产品所需设备台时及品的利润,生产单位产品所需设备台时及A,B两种原材料的消耗,现有原材两种原材料的消耗,现有原材料和设备台时的定额如表所示,问:料和设备台时的定额如表所示,问:)怎么安排生产使得工厂获利最大?)怎么安排生产使得工厂获利最大?)产品)产品的单位利润降低到的单位利润降低到1.8万元,要不要改变生产计划,如果降低到万元,要不要改变生产计划,如果降低到1万元呢?万元呢?)产品)产品的单位利润增大到的单位利润增大到5万元,要不要改变生产计划?万元,要不要改变生产计划?)如果
24、产品)如果产品,的单位利润同时降低了的单位利润同时降低了1万元,要不要改变生产计划?万元,要不要改变生产计划?产品产品最大资源量设备128台时原材料A4016kg原材料B0412kg单位产品利润23敏感性分析敏感性分析1 1 第35页,共115页,编辑于2022年,星期六敏感性分析敏感性分析1 1 第36页,共115页,编辑于2022年,星期六程序编写程序编写model:title 生产计划问题生产计划问题;maxfmax=2*x1+3*x2;TIMEx1+2*x28;A4*x116;B4*x212;END敏感性分析敏感性分析1 1 第37页,共115页,编辑于2022年,星期六运行结果运行结
25、果 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.000000 0.1250000 TIME 4.000000 0.000000 对问题对问题1,安排是生产产品,安排是生产产品4单位,产品单位,产品2单位,最大盈单位,最大盈利为利为14万元万元。敏感性分析敏感性分析1 1 第38页,共115页,编辑于2022
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 规划 建模 幻灯片
限制150内