数学建模优化模型介绍.ppt
《数学建模优化模型介绍.ppt》由会员分享,可在线阅读,更多相关《数学建模优化模型介绍.ppt(134页珍藏版)》请在三一文库上搜索。
1、数学建模数学建模优化模型介绍http:/ -Francis Bacon 数学处于人类智能的中心领域数学处于人类智能的中心领域数学方数学方法渗透、支配着一切自然科学的理论分支法渗透、支配着一切自然科学的理论分支它已愈来愈成为衡量成就的主要标志。它已愈来愈成为衡量成就的主要标志。-von Neumann http:/ 一门科学只有当它达到能够成功地运用一门科学只有当它达到能够成功地运用 数学时,才算真正发展了。数学时,才算真正发展了。-Karl MarxGalileo:展现在我们眼前的宇宙像一本用数学语言写成的大书,如不掌握数学符号语言,就像在黑暗的迷宫里游荡,什么也认识不清。数学是一种语言,是一
2、切科学的共同语言数学是一种语言,是一切科学的共同语言http:/ 成的一种关键性的、可实现的技术二十世纪最伟大的数学家-二十世纪最伟大的物理学家-D.HilbertA.EinsteinGo back诺贝尔诺贝尔奖奖菲尔兹奖菲尔兹奖http:/ 数学建模数学建模其实并不是什么新东西,可以说有了数学并需要用数学去解决实际问题,就一定要用数学的语言、方法去近似地刻画该实际问题,这种刻划的数学表述的就是一个数学模型,其过程就是数学建模的过程数学模型一经提出,就要用一定的技术手段(计算、证明等)来求解并验证,其中大量的计算往往是必不可少的,高性能的计算机的出现使数学建模这一方法如虎添翼似的得到了飞速的发
3、展,掀起一个高潮数学建模将各种知识综合应用于解决实际问题中,是培养和提高同学们应用所学知识分析问题、解决问题的能力的必备手段之一.http:/ 规划模型的应用极其广泛,其作用已为越来规划模型的应用极其广泛,其作用已为越来来越急速地渗透于工农业生产、商业活动、军事来越急速地渗透于工农业生产、商业活动、军事行为行为 科学研究的各个方面,为社会节省的财富、科学研究的各个方面,为社会节省的财富、创造的价值无法估量创造的价值无法估量.特别是在数模竞赛过程中,规划模型是最常特别是在数模竞赛过程中,规划模型是最常见的一类数学模型见的一类数学模型.从从92-06年全国大学生数模竞年全国大学生数模竞越多的人所重
4、视越多的人所重视.随着计算机的逐渐普及,它越随着计算机的逐渐普及,它越赛试题的解题方法统计结果来看,规划模型共出赛试题的解题方法统计结果来看,规划模型共出现了现了15次,占到了次,占到了50%,也就是说每两道竞赛题,也就是说每两道竞赛题中就有一道涉及到利用规划理论来分析、求解中就有一道涉及到利用规划理论来分析、求解.http:/ min f(x)-目标函数目标函数 s.t.s.t.x S-约束集合,可行集约束集合,可行集其中,其中,S R Rn n,f:S R R,x S称(称(f S)的可行解的可行解n最优解最优解:x*S,满足满足f(x*)f(x),x S。则称则称 x*为为(f S)的全
5、局最优解的全局最优解(最优解最优解),),记记 g.opt.(global optimum),),简记简记 opt.n最优值最优值:x*为为(f S)的最优解的最优解,则称则称 f*=f(x*)为为 (f S)的最优值的最优值(最优目标函数值最优目标函数值)数学规划模型的一般形式数学规划模型的一般形式(f S)http:/ N(x*),使满足,使满足 f(x*)f(x),x S N(x*)。则称则称 x*为为(f S)的局部最的局部最优解优解,记记 l.opt.(local optimum)n在上述定义中,当在上述定义中,当x x*时有严格不等式成立,时有严格不等式成立,则分别则分别称称 x*
6、为为(f S)的严格全局最优解和严格局部最优解的严格全局最优解和严格局部最优解。严格严格l.opt.严格严格g.opt.l.opt.http:/ 优化优化线性规划线性规划非线性规划非线性规划二次规划二次规划连续优化连续优化整数规划整数规划问题求解的难度增加http:/ 线性规划线性规划LinearProgramminghttp:/ 1x x1 1+c+c2 2x x2 2+.+c.+cn nx xn ns.t.as.t.a1111x x1 1+a+a1212x x2 2+.+a.+a1n1nx xn n (=,(=,)b)b1 1 a a2121x x1 1+a+a2222x x2 2+.+a
7、a2n2nx xn n (=,(=,)b)b2 2 .a am1m1x x1 1+a+am2m2x x2 2+.+a.+amnmnx xn n (=,(=,)b)bm m x x1,1,x x2 2.x.xn n 0 0http:/ 线性规划问题线性规划问题有可行解有可行解(Feasible)无无可可行行解解(Infeasible)有有最最优优解解(Optimal)无无最最优优解解(Unbounded)http:/ x1+x23002x1+x2400 x2250 x1、x20该问题的最优解为该问题的最优解为x1=50=50;x2=250=250 x2z*=50 x1+100 x2=2750
8、0 x1+x2300 x1x22502x1+x2400z1=50 x1+100 x2=0BOACDz2=14000http:/ 1.模型:命令:x=linprog(c,A,b)2.模型:minz=cX 命令:x=linprog(c,A,b,Aeq,beq)注意:若没有不等式:存在,则令A=,b=.式中:linprog称为调用函数,C,A,b称为输入参数,全部由用户提供,必须按规定的位置放置在原括号内.http:/ VLBXVUB命令:1x=linprog(c,A,b,Aeq,beq,VLB,VUB)2x=linprog(c,A,b,Aeq,beq,VLB,VUB,X0)注意:1若没有等式约束:
9、则令Aeq=,beq=.2其中X0表示初始点,设置它有些情况下可以减少迭代工作量4.命令:x,fval=linprog()返回最优解及处的目标函数值fval.http:/ 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0;0 0.02 0 0 0.05 0;0 0 0.03 0 0 0.08;b=850;700;100;900;Aeq=;beq=;vlb=0;0;0;0;0;0;vub=;x,fval=linprog(c,A,b,Aeq,beq,vlb,vub)ToMATLAB(xxgh1)http:/ 3 4;A=0 1 0;b=50;Aeq=1 1 1
10、beq=120;vlb=30,0,20;vub=;x,fval=linprog(c,A,b,Aeq,beq,vlb,vub)ToMATLAB(xxgh2)http:/ 9 10 11 12 8;A=0.4 1.1 1 0 0 0 0 0 0 0.5 1.2 1.3;b=800;900;Aeq=1 0 0 1 0 0 0 1 0 0 1 0 0 0 1 0 0 1;beq=400 600 500;vlb=zeros(6,1);vub=;x,fval=linprog(f,A,b,Aeq,beq,vlb,vub)ToMATLAB(xxgh3)http:/ 600.0000 0.0000 400.0
11、000 0.0000 500.0000fval=1.3800e+004即在甲机床上加工600个工件2,在乙机床上加工400个工件1、500个工件3,可在满足条件的情况下使总加工费最小为13800.http:/ 36;A=-5-3;b=-45;Aeq=;beq=;vlb=zeros(2,1);vub=9 15;%调用linprog函数:x,fval=linprog(c,A,b,Aeq,beq,vlb,vub)ToMATLAB(xxgh4)http:/ 0.0000fval=360即只需聘用9个一级检验员.注:注:本问题应还有一个约束条件:本问题应还有一个约束条件:x1、x2取整数取整数.故它是一
12、个整数故它是一个整数线性规划问题线性规划问题.这里把它当成一个线性规划来解,求得其最优解刚好这里把它当成一个线性规划来解,求得其最优解刚好是整数:是整数:x1=9,x2=0,故它就是该整数规划的最优解,故它就是该整数规划的最优解.若用线性规划若用线性规划解法求得的最优解不是整数,将其取整后不一定是相应整数规划的解法求得的最优解不是整数,将其取整后不一定是相应整数规划的最优解,这样的整数规划应用专门的方法求解最优解,这样的整数规划应用专门的方法求解.返回http:/ x2原料供应原料供应劳动时间劳动时间加工能力加工能力决策变量决策变量 目标函数目标函数 每天获利每天获利约束条件约束条件非负约束非
13、负约束线性线性规划规划模型模型(LP)时间时间480小时小时 至多加工至多加工100kgA150桶牛奶桶牛奶每天每天基本基本模型模型http:/ 图解法图解法 x1x20ABCDl1l2l3l4l5约约束束条条件件目标目标函数函数 z=0z=2400z=3600z=c(常数常数)等值线等值线c在在B(20,30)点得到最优解点得到最优解目标函数和约束条件是线性函数目标函数和约束条件是线性函数可行域为直线段围成的凸多边形可行域为直线段围成的凸多边形目标函数的等值线为直线目标函数的等值线为直线最优解一定在凸多边最优解一定在凸多边形的某个顶点取得。形的某个顶点取得。http:/ 软件实现软件实现 L
14、INDOmax72x1+64x2st2)x1+x2503)12x1+8x24804)3x1100endOBJECTIVEFUNCTIONVALUE1)3360.000VARIABLEVALUEREDUCEDCOSTX120.0000000.000000X230.0000000.000000ROWSLACKORSURPLUSDUALPRICES2)0.00000048.0000003)0.0000002.0000004)40.0000000.000000DORANGE(SENSITIVITY)ANALYSIS?No20桶牛奶生产桶牛奶生产A1,30桶生产桶生产A2,利润,利润3360元。元。ht
15、tp:/ OBJECTIVEFUNCTIONVALUE1)3360.000VARIABLEVALUEREDUCEDCOSTX120.0000000.000000X230.0000000.000000ROWSLACKORSURPLUSDUALPRICES2)0.00000048.0000003)0.0000002.0000004)40.0000000.000000max72x1+64x2st2)x1+x2503)12x1+8x24804)3x1100end三三种种资资源源“资源资源”剩余为零的约束为紧约束(有效约束)剩余为零的约束为紧约束(有效约束)原料无剩余原料无剩余时间无剩余时间无剩余加工能
16、力剩余加工能力剩余40http:/ OBJECTIVEFUNCTIONVALUE1)3360.000VARIABLEVALUEREDUCEDCOSTX120.0000000.000000X230.0000000.000000ROWSLACKORSURPLUSDUALPRICES2)0.00000048.0000003)0.0000002.0000004)40.0000000.000000最优解下最优解下“资源资源”增增加加1单位时单位时“效益效益”的的增量增量35元可买到元可买到1桶牛奶,要买吗桶牛奶,要买吗?3548,应该买!应该买!聘用临时工人付出的工资最多每小时几元聘用临时工人付出的工资
17、最多每小时几元?2元!元!原料增加原料增加1单位单位,利润增长利润增长48时间增加时间增加1单位单位,利润增长利润增长2加工能力增长加工能力增长不不影响利润影响利润影子价格影子价格http:/ 3=72增至增至30 3=90”(或(或“=”(或(或“=”)功能相同)功能相同2.变量与系数间可有空格变量与系数间可有空格(甚至回车甚至回车),但无运算符但无运算符3.变量名以字母开头,不能超过变量名以字母开头,不能超过8个字符个字符4.变量名不区分大小写(包括变量名不区分大小写(包括LINDO中的关键字)中的关键字)5.目标函数所在行是第一行,第二行起为约束条件目标函数所在行是第一行,第二行起为约束
18、条件6.行号行号(行名行名)自动产生或人为定义自动产生或人为定义.行名以行名以“)”结束结束7.行中注有行中注有“!”符号的后面部分为注释符号的后面部分为注释.如如:!ItsComment.8.在模型的任何地方都可以用在模型的任何地方都可以用“TITLE”对模型命名对模型命名(最多(最多72个字符),如:个字符),如:TITLEThisModelisonlyanExamplehttp:/ x1+x2=345.5;x1=98;x1=98;2*x1+x2=600 2*x1+x2=345.5 x1+x2=345.5 x1=98 x1=98 2*x1+x2=600 2*x1+x2bj,其数学模型为:h
19、ttp:/ i bbj j;若若用用xi,n+1表表示示从从Ai到到Bn+1的的运运量量,可可令令ci,n+1=0或或等等于于第第Ai产产地地储储存存单单位位物物资资的的费费用用。因因为为xi,n+1实实际际上上表表示示Ai产产地地没没有有运运出出去去而而库库存存的的物物资资数数量量。经经处处理理后后,问问题题变变成成了了产产销销平平衡衡的的运输问题,其数学模型为:运输问题,其数学模型为:这样,这样,m个产地、个产地、n个销地的不平衡运输问题,转化成了个销地的不平衡运输问题,转化成了m个产地、个产地、n+1个销地的平衡运输问题,此时可用表上作业法求解。个销地的平衡运输问题,此时可用表上作业法求
20、解。http:/ bj,其数学模型为:http:/ 某公司有某公司有6个供货栈(仓库),库存货物总数分别为个供货栈(仓库),库存货物总数分别为60,55,51,43,41,52,现有,现有8个客户各要一批货数量分别为个客户各要一批货数量分别为35,37,22,32,41,32,43,38,各供货栈道,各供货栈道8个客户的单位货物运输价见表个客户的单位货物运输价见表供货栈到客户的单位货物运价供货栈到客户的单位货物运价客客户货栈V1V2V3V4V5V6V7V8W162674259W249538582W352197433W476739271W523957265W655228143试确定各货栈到各客户
21、处的货物调运数量,使总的运输费用最小。试确定各货栈到各客户处的货物调运数量,使总的运输费用最小。http:/ 引入决策变量引入决策变量代表从第代表从第个货栈到第个货栈到第个个客户的货物运量客户的货物运量.设设表示从第表示从第 个货栈个货栈到第 个客户的单位货物运价表示第表示第 个货栈的最大供货量个货栈的最大供货量,表示第表示第个客户的订货量个客户的订货量.目标函数是总运输费用最少目标函数是总运输费用最少.约束条件有三条:约束条件有三条:1.各货栈运出的货物总量不超过其库存数各货栈运出的货物总量不超过其库存数 2.各客户收到的货物总量等于其订货数量各客户收到的货物总量等于其订货数量3.决策变量决
22、策变量 非负非负数学模型为数学模型为http:/ 前面介绍的前面介绍的LinGO的基本用法,其优点是输入模型较直观,的基本用法,其优点是输入模型较直观,一般的数学表达式无需作大的变换即可直接输入一般的数学表达式无需作大的变换即可直接输入.对于规模较对于规模较小的的规划模型,用直接输入的方法是有利的,如果模型的小的的规划模型,用直接输入的方法是有利的,如果模型的变量和约束的条件个数比较多,若仍然用直接的输入方式,变量和约束的条件个数比较多,若仍然用直接的输入方式,虽然也能求解并得到结果,但这种做法有明显的不足之处。虽然也能求解并得到结果,但这种做法有明显的不足之处。模型的篇幅很长,不便于分析修改
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 建模 优化 模型 介绍
