数学建模线性和非线性规划

《数学建模线性和非线性规划》由会员分享，可在线阅读，更多相关《数学建模线性和非线性规划（61页珍藏版）》请在金锄头文库上搜索。

1、实验目的 1) 了解最优化问题的基本结构和基本建模方法; 2) 线性规划的求解方法; 3） 非线性规划的求解方法.一,优化问题的普遍性以及引例1，无处不在的优化每一个人，高致总统首相，总裁经理，平民百姓，无不在做决策：该做什么，该怎么做，才能有最好的效果?甚至自然中的动植物，也时刻面临这样的问题.类似的问题，还广泛的存在于无机世界中.一,优化问题的普遍性以及引例看看下面的例子分别属于哪一类？a)证券的投资组合；b)国家经济发展战略；c)产品规格、性能设计；d)球形的水滴；e)狼群的集体捕食；f)好的购物方案；g)物质分子结构； h)生物的身体构造；i)乘务组排班表; j)光传播路径:直线,反射

2、,折射课堂作业：和你的同桌讨论还有什么方面需要优化的。一,优化问题的普遍性以及引例2,一些成功的优化例子：“最优人员安排”为美国航空每年节约两千万美元. “改进的出货流程”每年为Yellow Freight 公司节约一千七百多万美元.“改进的卡车分派”为 Reynolds 公司每年节约七百万美元 .最优全局供应链为数字设备行业节约超过三亿美元.重建的 North America Operations, Proctor and Gamble减少 20%的工厂, 每年节约两亿美元.大阪的Hanshin高速的最优安排每年节约一千七百万人小时 . 为说明最优化的价值,建立了专门的网站,列举了哪些公司的

3、什么问题,运用最优化方法节约和增加了多少金额. 有可选的行业,考察的方面,受益的方式,希望同学们各选择其中的一个,提一份报告,以说明最优化的价值.一,优化问题的普遍性以及引例Google上相关搜索的结果：Search phrase number of hits（英文） 短语 点击数（中文）“optimize the supply chain”1,160,000 优化供应链414,000“optimize (the) return” 2,490,000 优化回报 453,000“optimal experience” 32,400,000 最优经历“optimal investment” 8,3

4、20,000 优化投资 8,250,000“optimal system” 84,200,000 优化系统 13,800,000“optimal decision” 28,800,000 最优决策 2,890,000“optimize your PC” 3,300,000 优化你的PC“optimal choice” 25,800,000 最优选择10,900,000“optimal design” 77,300,000 优化设计 1,270,000“optimal health” 31,900,000 优化健康还有如:优化产业结构 2,830,000 优化人员结构 3,110,000同学们有

5、没有发现,英文和中文短语间有很大的不同,原因可能是什么?一,优化问题的普遍性以及引例3,相关的几句格言：Waste neither time nor money, but make the best use of both. - Benjamin Franklin Obviously, the highest type of efficiency is that which can utilize existing material to the best advantage. -Jawaharlal NehruIt is more probable that the average man

6、could, with no injury to his health, increase his efficiency fifty percent.-Walter Scott请同学翻译上面的句子，你喜欢那一句？你有什么好的表述？引例1,动物饲料配置问题 美国一家公司以专门饲养并出售一种实验用的动物而闻名。这种动物的生长对饲料中的三种营养成分特别敏感，即蛋白质、矿物质和维生素。需要的营养量蛋白质：70克矿物质：3克维生素：9.1毫克 现有五种饲料，公司希望找出满足动物营养需要使成本达到最低的混合饲料配置。 饲料蛋白质(克)矿物质(克)维生素(毫克)1(x1)2(x2)3(x3)4(x4)5(x

7、5)需要量0.302.001.000.601.80700.100.050.020.200.0530.050.100.020.200.089.1每一种饲料每磅所含的营养成分每种饲料每磅的成本饲料12345成本(美元)0.020.070.040.030.05引例2:供应与选址 某公司有某公司有某公司有某公司有6 6个建筑工地要开工，每个工地的位个建筑工地要开工，每个工地的位个建筑工地要开工，每个工地的位个建筑工地要开工，每个工地的位置（用平面坐标置（用平面坐标置（用平面坐标置（用平面坐标a a，b b表示，距离单位：千米）及表示，距离单位：千米）及表示，距离单位：千米）及表示，距离单位：千米）及水

8、泥日用量水泥日用量水泥日用量水泥日用量d d 吨由下表给出。目前有两个临时料场吨由下表给出。目前有两个临时料场吨由下表给出。目前有两个临时料场吨由下表给出。目前有两个临时料场位于位于位于位于A(5,1), B(2,7)A(5,1), B(2,7)，日储量各有日储量各有日储量各有日储量各有2020吨。假设从料场吨。假设从料场吨。假设从料场吨。假设从料场到工地均有直线道路相连，（到工地均有直线道路相连，（到工地均有直线道路相连，（到工地均有直线道路相连，（1 1）试制定每天的供）试制定每天的供）试制定每天的供）试制定每天的供应计划，即从应计划，即从应计划，即从应计划，即从A A、B B 两料场分别

9、向各工地运送多少两料场分别向各工地运送多少两料场分别向各工地运送多少两料场分别向各工地运送多少吨水泥，使总的吨千米数最小。吨水泥，使总的吨千米数最小。吨水泥，使总的吨千米数最小。吨水泥，使总的吨千米数最小。a1.258.750.55.7537.25b1.250.754.7556.57.75d3547611二,优化问题建模的基本步骤介绍 在我们的生活中,始终有这样的问题:为了一定的目的做一些事情,我们可能要考虑有哪些重要的因素,这些因素和要完成的目标之间有什么样的关系.也就是说,我们在做一个决定时,会注意下面的三个要点:目的是什么?有哪些重要的因素?这些因素和你的目标之间有什么样的关系?二，优化

10、问题的表述目标函数 对应决策者而言,对其有利的程度必须定量的测度, 在商业应用中，有效性的测度经常是利润或者成本, 但对于政府，更经常的使用投入产出率来测度. 表示有效性测度的经常称为目标函数.目标函数要表出测度的有效性, 必须说明测度和导致测度改变的变量之间的关系. 系统变量分为决策变量和参数.决策变量是指能由决策者直接控制的变量. 而参数是指不能由决策者决定的量.实际上，数学模型很少有能表达变量和有效性测度之间的精确关系的. 实际上，运筹学分析者的任务就是找出对测度有最重要影响的变量 然后找出这些变量和测度之间的数学关系.这个数学关系也就是目标函数.二，优化问题的表述决策变量和参数 我们称

11、对应决策者可控的量称为决策变量,决策变量的取值确定了系统的最终性能,也是决策者采用决策的依据.在系统中还有一些量,它不能由决策者所控制,而是由系统所处的环境所决定,我们称之为参数.二，优化问题的表述约束条件 约束条件就是决策变量和参数之间的关系. 约束集界定决策变量可以取某些值而不能取其他的值.比如对应生产问题, 任何活动中,时间和物品不能为负数.当然,也有一些优化问题不带约束条件,我们称之为无约束优化问题.而在实际问题中,决策变量带有约束是普遍的. 三,优化问题的分类优化问题的分类可以从几个方面进行:1,从变量取值的连续和离散可以分成:连续优化,离散优化和混合优化2,从问题的线性非线性可以分

12、为:线性规划和非线性规划3,从变量是确定性和随机性可以分为:随机规划和确定性问题.John Von Neumann 约翰冯诺依曼(19031957),美藉匈牙利人.20世纪最杰出的数学家之一,被誉为”计算机之父”,”博弈论之父”.被认为是数学规划的三大创始人之一.以下的三个人物和线性规划的出现有重要的关系.George B. Dantzig George B. Dantzig(1914-2005),美国人,线性规划单纯形法的创始人,被誉为”线性规划之父”.美国科学院三院院士,美国军方数学顾问,教授.并以其名字设立Dantzig奖.数学规划的三大创始人之一.发现算法时非常年轻,以至到日本时,人们

13、以为”线性规划之父”是个老人,而对他无人问津.Leonid Vitalyevich Kantorovich Kantorovich(1912-1986)苏联人,著名数学家和经济学家,教授,年仅18岁获博士学位.因在经济学上提出稀缺资源的最优配置获诺贝尔奖.线性规划对偶理论的提出者,数学规划的三大创始人之一.非线性规划问题在实践中也是及其常见的.标志着这一学科的产生的奠基性工作由美国的数学家Tucker和Kuhn在1952年的一篇文章.该文章给出了非线性规划问题的必要条件和充分条件,后来成为Kuhn-Tucker条件.这为非线性规划问题的求解算法的提出提供了理论基础和算法的基本思路.相关的规划问

14、题,比如多目标规划,决策论等等. 美国一家公司以专门饲养并出售一种实验用的动物而闻名。这种动物的生长对饲料中的三种营养成分特别敏感，即蛋白质、矿物质和维生素。需要的营养量蛋白质：70克矿物质：3克维生素：9.1毫克 现有五种饲料，公司希望找出满足动物营养需要使成本达到最低的混合饲料配置。 四,线形规划问题的解法及举例饲料蛋白质(克)矿物质(克)维生素(毫克)1(x1)2(x2)3(x3)4(x4)5(x5)需要量0.302.001.000.601.80700.100.050.020.200.0530.050.100.020.200.089.1每一种饲料每磅所含的营养成分每种饲料每磅的成本饲料1

15、2345成本(美元)0.020.070.040.030.05建立数学模型决策变量：在混合饲料中，每天所需第j种饲料的磅数xj，j = 1，2，3，4，5； 约束条件： 蛋白质：0.30x1+2x2+x3+0.6x4+1.8x570 矿物质：0.10x1+0.05x2+0.02x3+0.2x4+0.05x53 维生素：0.05x1+0.1x2+0.02x3+0.2x4+0.08x510自然约束条件：xi0 确定目标：混合饲料的成本最低 0.02x1+0.07x2+0.04x3+0.03x4+0.05x5min完整的线性规划模型：min 0.02x1+0.07x2+0.04x3+0.03x4+0.

16、05x5s.t. 0.30x1+2x2+x3+0.6x4+1.8x570 0.10x1+0.05x2+0.02x3+0.2x4+0.05x53 0.05x1+0.1x2+0.02x3+0.2x4+0.08x510 xj0 j = 1,2,3,4,5;min cTxs.t. Axb x0归纳：返 回linprog min cTx s.t. Axb Aeqx beq lb x ubSolve a linear programming problem where c, x, b, beq, lb, and ub are vectors and A and Aeq are matrices.调用格式:

17、x = linprog(f,A,b,Aeq,beq)x = linprog(f,A,b,Aeq,beq,lb,ub)x = linprog(f,A,b,Aeq,beq,lb,ub,x0)x = linprog(f,A,b,Aeq,beq,lb,ub,x0,options)x,fval = linprog(.)x,fval,exitflag = linprog(.)x,fval,exitflag,output = linprog(.)x,fval,exitflag,output,lambda = linprog(.)原油生产计划原油生产计划原油类别原油类别 买入价买入价(元元/桶桶)买入量买入量

18、(桶桶/天天)辛烷值辛烷值(%)硫含量硫含量(%)A455000120.5B35500062.0C25500083.0汽油类别汽油类别卖出价卖出价(元元/桶桶)需求量需求量(桶桶/天天)辛烷值辛烷值(%)硫含量硫含量(%)甲甲703000101.0乙乙60200082.0丙丙50100061.01:1 加工费: 4元/桶 能力: = 14000桶/天I: 安排生产计划，在满足需求的条件下使利润最大 决策变量：目标：甲甲(3000)乙乙(2000)丙丙(1000)A/45X1X2X3B/35X4X5X6C/25X7X8X9约束：总利润最大 需求限制； 原料限制； 含量限制； 非负限制 含量限制非

19、负限制 原料限制需求限制约束总盈利：总盈利：126000元 c=454545353535252525;a1=100100100; 010010010; 001001001;a2=111000000; 000111000; 000000111; -1200-600-800; 0-1200-600-80; 00-1200-600-8; 0.500200300; 00.50020030; 0 00.5002003;b1=3000 20001000;b2=5000 50005000-30000-16000-6000300040001000;v1=zeros(1,9);x f=linprog(c,a2,

20、b2,a1,b1,v1)z=356000-f甲甲(3000)乙乙(2000)丙丙(1000)A/452400800800B/35000C/256001200200II: 通过广告增加销售(1元广告费：增加10桶销售)决策变量：目标：甲甲(3000+)乙乙(2000+)丙丙(1000+) A/45X1X2X3B/35X4X5X6C/25广告广告销售销售X7X103000+10X10X8X112000+10X11X9X121000+10X12约束：总利润最大需求限制； 原料限制； 产量限制；含量限制； 非负限制 含量限制非负限制 产量限制原料限制需求限制约束总盈利：总盈利：287750元 c=49

21、 49 49 39 39 39 29 2929-699-599-499;a1=1 0 0 100100 -100 0; 0 1 0 010010 0 -10 0; 0 0 1 001001 00 -10;a2=1110000 0 0 0 0 0; 0001110 0 0 0 0 0; 0000001 1 1 0 0 0; -1200-600-8 0 0 1000 0; 0-1200-600 -8 0080 0; 00-1200-60 0 -800 60; 0.5002003 0 0-100 0; 00.500200 3 00-20 0; 000.50020 0 300 -10 000000 0

22、 001 1 1;b1=300020001000;b2=500050005000-30000-16000-6000 3000 4000 1000800;v1=zeros(1,12);x f=linprog(c,a2,b2,a1,b1,v1)z=380000-f甲甲(3000)乙乙(2000)丙丙(1000)A/452121.82185.3692.9B/35695.54036.8267.6C/25广告广告182.703277.9750 39.40某公司有某公司有6个建筑工地，位置坐标为个建筑工地，位置坐标为(ai, bi) (单位：公单位：公里里),水泥日用量水泥日用量di (单位：吨）单位：吨

23、）假设：假设：料场料场和工地之间和工地之间有直线道路有直线道路五,非线形规划问题的解法及举例用例中数据计算，最优解为总吨公里数为总吨公里数为总吨公里数为总吨公里数为136.2136.2线性规划模型线性规划模型决策变量：决策变量：ci j (料场料场j到到工地工地i的的运量）运量）12维维Shili084lin.m选址问题：选址问题：NLPNLP2）改建两个新料场，需要确定新料场位置）改建两个新料场，需要确定新料场位置(xj,yj)和运量和运量cij ，在其它条件不变下使总吨公里数最小。在其它条件不变下使总吨公里数最小。决策变量：决策变量：ci j，(xj,yj)16维维非线性规划模型非线性规划

24、模型结果：结果：总吨公里数总吨公里数为为85.3，但局部最优，但局部最优解问题严重解问题严重Shili084.m: shili084fun.m决策变量：决策变量：ci ，(xj,yj)10维维 计算方法的改善计算方法的改善局部最优解问题有所改进+为工地, 数字为用量; *为新料场, 数字为供应量。约束非线性规划情形约束非线性规划情形标准模型标准模型fmincon语句的具体用法语句的具体用法x,fval,exitflag,output=fmincon(fun,x0,A,b,Aeq,beq,lb,ub,con) 建立建立m m文件函数文件函数 function f,G=fun(x)function

25、 f,G=fun(x) f=f(x); f=f(x); G=G G=G1 1(x),G(x),G2 2(x);(x); 选项选项optionsoptions 输出参数输出参数options(8)=options(8)=目标函数最优值；目标函数最优值； 输入参数输入参数options(13)=options(13)=等式约束的个数等式约束的个数(m)(m)； L L，U U是决策变量的下界和上界；是决策变量的下界和上界；2、fmincon语句的具体用法语句的具体用法max f(x) = x12+ x22-x1x2-2x1-5x2 s.t. -(x1 1)2+ x2 0 2 x1-3x2+60,

26、x0=0, 1例2fmincon语句的具体用法语句的具体用法转化成标准形转化成标准形min f(x) =- x12- x22+x1x2+2x1+5x2 s.t. (x1 1)2 - x2 0 -2 x1+3x2 60, x0=0, 1 function f=fun22(x) f=-x(1)2-x(2)2+x(1)*x(2)+2*x(1)+5*x(2); function G,Geq = cont2(x) G=(x(1)-1)2 - x(2); Geq=; x0=0 1; A=-2,3;b=6; Aeq=;beq=; lb=;ub=; fun22x,fval=fmincon(fun22,x0,A

27、,b,Aeq,beq,lb,ub,cont2) x = 1.0e+008*-0.0006 -2.7649 fval = -7.6432e+016fun22= -x(1)2-x(2)2+x(1)*x(2)+2*x(1)+5*x(2)fmincon语句的具体用法语句的具体用法例fmincon语句的具体用法语句的具体用法function fm=ex02(x)n=10;c=-6.089,-17.164,-34.054,-5.914, -24.721, -14.986,-24.1,-10.708, -26.662,-22.179;sx=0;for i=1:n sx=sx+x(i); %endfm=0;f

28、or i=1:n fm=fm+x(i)*(c(i)+log(x(i)/sx);end1) 建立目标函数建立目标函数 ex02.m供应与选址 某公司有某公司有某公司有某公司有6 6个建筑工地要开工，每个工地的位个建筑工地要开工，每个工地的位个建筑工地要开工，每个工地的位个建筑工地要开工，每个工地的位置（用平面坐标置（用平面坐标置（用平面坐标置（用平面坐标a a，b b表示，距离单位：千米）及表示，距离单位：千米）及表示，距离单位：千米）及表示，距离单位：千米）及水泥日用量水泥日用量水泥日用量水泥日用量d d 吨由下表给出。目前有两个临时料场吨由下表给出。目前有两个临时料场吨由下表给出。目前有两个

29、临时料场吨由下表给出。目前有两个临时料场位于位于位于位于A(5,1), B(2,7)A(5,1), B(2,7)，日储量各有日储量各有日储量各有日储量各有2020吨。假设从料场吨。假设从料场吨。假设从料场吨。假设从料场到工地均有直线道路相连，（到工地均有直线道路相连，（到工地均有直线道路相连，（到工地均有直线道路相连，（1 1）试制定每天的供）试制定每天的供）试制定每天的供）试制定每天的供应计划，即从应计划，即从应计划，即从应计划，即从A A、B B 两料场分别向各工地运送多少两料场分别向各工地运送多少两料场分别向各工地运送多少两料场分别向各工地运送多少吨水泥，使总的吨千米数最小。吨水泥，使总

30、的吨千米数最小。吨水泥，使总的吨千米数最小。吨水泥，使总的吨千米数最小。a1.258.750.55.7537.25b1.250.754.7556.57.75d3547611 （2 2）为进一步减少吨千米数，打算舍弃两个）为进一步减少吨千米数，打算舍弃两个）为进一步减少吨千米数，打算舍弃两个）为进一步减少吨千米数，打算舍弃两个临时料场，改建两个新的，日储量仍各为临时料场，改建两个新的，日储量仍各为临时料场，改建两个新的，日储量仍各为临时料场，改建两个新的，日储量仍各为2020吨，吨，吨，吨，问应建在何处，节省的吨千米数有多大？问应建在何处，节省的吨千米数有多大？问应建在何处，节省的吨千米数有多大

31、？问应建在何处，节省的吨千米数有多大？供应与选址范例：供应与选址(a1, b1), d1 = 3(a2, b2), d2 =5(a6, b6), d6= 11A(x1, y1)e1 = 20B(x2, y2)e2 = 20x11x12x16x21x26x22目标：使得吨公里min建立规划模型 记工地的位置为记工地的位置为记工地的位置为记工地的位置为( (a ai i, b, bi i) )，水泥日用量为水泥日用量为水泥日用量为水泥日用量为d di i, i = , i = 1,6,1,6,料场位置为料场位置为料场位置为料场位置为( (x xj j, , y yj j), ), 日储量为日储量为

32、日储量为日储量为e ej j, j=1,2, j=1,2；从料从料从料从料场场场场j j向工地运送量为向工地运送量为向工地运送量为向工地运送量为x xij ij（决策变量）。决策变量）。决策变量）。决策变量）。20吨吨 当料场位置为当料场位置为(xi, yi)为已知时为已知时, 上述模型上述模型为线性规划模型。从料场为线性规划模型。从料场i向工地运送量为向工地运送量为xij（决策变量）。决策变量）。 当料场位置为当料场位置为(xi, yi)为未知时为未知时, 上述模型上述模型为非线性规划模型。对决策变量为非线性规划模型。对决策变量(xi, yi)而言。而言。如何求解？如何求解？其中其中a，b，

33、d，e 都是常数向量（数组）。都是常数向量（数组）。 当料场位置为当料场位置为(xj, yj)为已知时为已知时, 上述模型为线性上述模型为线性规划模型。从料场规划模型。从料场j向工地运送量为向工地运送量为xij（决策变量）。决策变量）。思考： 矩阵形式 当料场位置为当料场位置为(xi, yi)为未知时为未知时, 上述模上述模型为非线性规划模型。对决策变量型为非线性规划模型。对决策变量(xi, yi)而言。而言。线性规划求解结果：线性规划求解结果：最优目标值最优目标值f = 136.23（吨千米）吨千米）j123456x1j(料场A)350701x2j(料场B)0040610程序程序liaocx

34、.mfunction f,g=liaoch(x)a=1.25 8.75 0.5 5.75 3 7.25 ;b=1.25 0.75 4.75 5 6.5 7.75 ;d=3 5 4 7 6 11;e=20 20; f1=0;for j=1:6s(j)=sqrt(x(13)-a(j)2+(x(14)-b(j)2);f1= f1 + s(j)*x(j);%A(x13,x14)到各工地的吨公里数endf2=0;for i=7:12s(i)=sqrt(x(15)-a(i-6)2+(x(16)-b(i-6)2);f2= f2 + s(i)*x(i);%B(x15,x16)到各工地的吨公里数endf=f1+

35、f2;% 总的目标函数for i=1:6 g(i)=x(i)+x(i+6)-d(i);%第一组约束条件endg(7)=sum(x(1:6)-e(1);g(8)=sum(x(7:12)-e(2); %第二组约束条件x0=ones(1,12);L=zeros(1,12);opt(13)=6;x,opt=constr(liaoch,x0,opt,L)f=opt(8),n=opt(10)MATLAB运行程序运行程序(liaochj.m)i123456ABCi1(料场A)0507085.906x3xCi2(料场B)3040635.073y6.5y计算结果：最优目标值计算结果：最优目标值f = 93.14

36、（吨千米）吨千米）新料场位置的改变，目标值比原来减少了新料场位置的改变，目标值比原来减少了43.1吨千米。吨千米。原最优目标值原最优目标值f = 136.23六,实验布置实验目的 1) 了解最优化问题的基本结构和基本建模方法; 2) 线性规划的求解方法; 3） 非线性规划的求解方法. 1,(糖果生产计划)某糖果公司用三种原料: 糖/果仁/巧克力生产两种糖果A和B.所有重量单位为公斤,所有价格单位为元 .假定原料混合过程中重量既不增加也不减少.三种原料分别只有100/20/30公斤,并要求糖果A中的果仁和巧克力含量不少于10%,糖果B中的果仁含量不少于20%.糖果A和B的单价分别为25元和20元

37、.问如何安排生产才能得到最高的收益.试建立相应的数学模型并进行求解. 2,(运输问题)从Toronto和Detroit两市分别有两批货物途径Chicago和Buffalo最终到达New York、Phila.和St.louis市.之间的路线表述如下图:ToFromChicagoBuffaloSupplyTorontoDetroit$4$ 5$ 7$ 7600500ToFromNew YorkPhila.St.louisChicagoBuffaloDemand$3$ 1450$ 2$ 3350$ 2$ 4300 其中Toronto和 Detroit 分别有600和500的货物需要运出,New Y

38、ork、Phila. 和St.louis的货物需求分别是450、350和300.每一段上的运输单价如下面两表： 问：如何进行运输安排使整个的运输费用最少？试建立问题的数学模型并求出最有解。3、求解非线性规划、求解非线性规划精品课件精品课件！精品课件精品课件！4、桃李花园服务中心选址、桃李花园服务中心选址初始点x0=20，20设设(ai, bi)(i=120)为第为第i 栋住宅楼的坐标；栋住宅楼的坐标；a = 29.74 4.9 69.32 65.0 98.3 55.27 40.0 19.8 62.5 73.3 37.58 0.98 41.98 75.37 79.38 92.0 84.47 36.77 62.08 73.13;b = 19.39 90.48 56.92 63.18 23.44 54.88 93.16 33.5 65.5 39.19 62.73 69.9 39.72 41.37 65.52 83.75 37.16 42.52 59.46 56.58;