机械优化设计6线性规划.ppt
第六章第六章第六章第六章 线性规划线性规划线性规划线性规划一一.线性规划的基本概念线性规划的基本概念二二.求解线性规划的单纯形法求解线性规划的单纯形法三三.初始基本可行解初始基本可行解2021/9/231 某某厂厂生生产产甲甲、乙乙两两种种产产品品,已已知知:两两种种产产品品分分别别由由两两条条生生产产线线生生产产。第第一一条条生生产产甲甲,每每天天最最多多生生产产9 9件件,第第二二条条生生产产乙乙,每每天天最最多多生生产产7 7件件;该该厂厂仅仅有有工工人人2424名名,生生产产甲甲每每件件用用2 2工工日日,生生产产乙乙每每件件用用3 3工工日日;产产品品甲甲、乙乙的的单单件件利利润润分分别别为为4040元元和和8080元元。问问工工厂厂如如何何组组织织生生产产才才能能获获得得最最大大利润?利润?一)应用实例一)应用实例6-1 6-1 6-1 6-1 线性规划的基本概念线性规划的基本概念线性规划的基本概念线性规划的基本概念2021/9/232日利润最大日利润最大生产能力限制生产能力限制劳动力限制劳动力限制变量非负变量非负解解:设甲、乙两种产品的日产件数分别为设甲、乙两种产品的日产件数分别为s.t.2021/9/233二二)线性规划的一般形式线性规划的一般形式s.t.s.t.特点特点:1)1)为极小化问题为极小化问题;2);2)约束取等号约束取等号;3)3)限定系数非负限定系数非负;4);4)变量非负变量非负.式中,式中,价值系数;价值系数;结构系数结构系数 限定系数限定系数2021/9/234将数学模型化为标准型的方法将数学模型化为标准型的方法1 1)将极大化问题化为极小化问题)将极大化问题化为极小化问题 松弛变量松弛变量(开关变量)(开关变量)(两边乘(两边乘-1-1)4 4)将负的限定系数化为正值)将负的限定系数化为正值3 3)将任意变量化为非负变量)将任意变量化为非负变量2 2)将不等式约束变为等式约束:)将不等式约束变为等式约束:目标函数变号;目标函数变号;2021/9/235s.t.化为标准型化为标准型:2021/9/236三)线性规划的基本概念三)线性规划的基本概念s.t.1.1.线性规划的图解线性规划的图解x2x10F=0F*=620(1.5,7)2021/9/2372.2.线性规划的基本概念线性规划的基本概念1 1)可行解)可行解满足约束条件及非负条件的解。满足约束条件及非负条件的解。(D D内及其边界上的解)内及其边界上的解)2 2)基本解)基本解使使n-mn-m个变量等于个变量等于0 0,解约束方程,解约束方程组组(共有共有m m个约束方程个约束方程)所得的解。所得的解。基本解对应于约束边界的交点基本解对应于约束边界的交点.3 3)基本可行解)基本可行解可行域中的基本解(即可行域中的基本解(即D D的顶点)。的顶点)。4 4)基本变量与非基本变量)基本变量与非基本变量 预先取为零值的预先取为零值的n-mn-m个变量为非个变量为非基本变量,其余基本变量,其余m m个为基本变量。个为基本变量。x2x10F=0F*=-620(1.5,7)s.t.2021/9/238四)线性规划的基本性质四)线性规划的基本性质 1 1)可行域)可行域D D为凸集,每个基本可行解对应于为凸集,每个基本可行解对应于D D上的上的一个顶点;一个顶点;2 2)只要可行域存在且封闭,则起码有一个基本可)只要可行域存在且封闭,则起码有一个基本可行解为最优点;行解为最优点;*)若最优点所在的边界线与等值线平行,则该)若最优点所在的边界线与等值线平行,则该边界线上的点均为最优点;边界线上的点均为最优点;)若可行域不封闭,则可能有无界解。)若可行域不封闭,则可能有无界解。3)3)最优点可在最优点可在D D的顶点中寻找。的顶点中寻找。2021/9/2396-2 6-2 6-2 6-2 求解线性规划的单纯形法求解线性规划的单纯形法求解线性规划的单纯形法求解线性规划的单纯形法一一.基本思路基本思路 先先取取D D的的一一个个顶顶点点作作为为初初始始点点,由由此此出出发发朝朝可可使使目目标标函函数数降降低低最最快快的的方方向向依依次次经经过过一一系系列的基本可行解列的基本可行解,直至达到最优解直至达到最优解.*1)*1)需获得一个初始基本可行解需获得一个初始基本可行解;2)2)每次只更换一个非基本变量每次只更换一个非基本变量;3)3)保证下降性和可行性保证下降性和可行性.2021/9/2310二二.计算实例计算实例s.t.s.t.1.1.初始基本可行解初始基本可行解取取x x5 5,x,x6 6 为基本变量为基本变量,则有则有:0 0 0 0 4 5T2021/9/23112.2.第一次变换顶点第一次变换顶点(1)(1)选取进基变量选取进基变量原则原则:考虑下降性考虑下降性,且下降得最快且下降得最快判别数判别数:假定假定x x2 2进基进基,则有则有取取相应的目标函数变化量相应的目标函数变化量:即即2021/9/2312写成一般形式写成一般形式:最小最小,x,x3 3 应为进基变量应为进基变量推论推论:若若线线性性规规划划的的一一个个基基本本可可行行解解的的所所有有进进基基判判别别数数均均为非负为非负,则该解为最优解则该解为最优解.2021/9/2313(2)(2)确定离基变量确定离基变量原则原则:考虑可行性考虑可行性(该变量离基后该变量离基后,能使余下的基本变量为非负能使余下的基本变量为非负)判别数判别数:由于由于)若取若取 (离基离基),),则有则有 应取应取 为正且其值为最小者对应的基本变量离基为正且其值为最小者对应的基本变量离基.(可行可行)(不可行不可行)若取若取 (离基离基),),则有则有 2021/9/2314)推推论论:若若线线性性规规划划的的的的所所有有离离基基判别数均为负数时判别数均为负数时,则问题有无界解则问题有无界解.最小最小,x,x6 6 应为离基变量应为离基变量0 0 5/3 0 2/3 0T*)因因为为 ,故故 也也必必须须大大于于0,0,否否则则不不满满足足可行性要求可行性要求;2021/9/2315进基进基3.3.第二次变换顶点第二次变换顶点去掉了去掉了(1)(2)1)1)确定进基变量确定进基变量(3)(4)2021/9/23162)2)确定离基变量确定离基变量 离基离基(1)(1)(2)(2)0 0 8/5 1/5 0 0T(3)(3)(4)(4)2021/9/23174.4.第三次变换顶点第三次变换顶点1)1)确定进基变量确定进基变量 故故 为最优点为最优点,为最优值为最优值:0 0 8/5 1/5 0 0T2021/9/2318三三.用单纯形表求解线性规划用单纯形表求解线性规划例例.用初等变换法求解用初等变换法求解解解:增广矩阵增广矩阵:2021/9/2319s.t.离离基基判判别别数数进基判别数进基判别数 单单纯纯形形法法实实际际上上是是解解一一系系列列的的线线性性方方程程组组,也也可可用用初初等等变换方法列表求解变换方法列表求解.但需加入判别数的计算但需加入判别数的计算.421235基变量基变量x1x2x3x4x5x63x5112410425x612310155/3X0000045F037-4-11-20-15例例1 12021/9/232042123基变量基变量x1x2x3x4x5x63x51/3-1/3010/312/30.21x31/32/311/305/35X1005/302/30F111/38/37/3-25/3421235基变量基变量x1x2x3x4x5x63x5112410425x612310155/3X0000045F037-4-11-20-152021/9/232142123基变量基变量x1x2x3x4x5x63x51/3-1/3010/312/30.21x31/32/311/305/35X1005/302/30F111/38/37/3-25/34212基变量基变量x1x2x3x4x5x62x41/10-1/10010.21x33/107/10101.6X2001.60.200F223.51.5已获得最优解已获得最优解2021/9/2322-2-300基变量基变量x1x2x3x40 x3-1110330 x41-4014-1X00034F00-2-3-2-30基变量基变量x1x2x3x4-3x2-1103-30 x4-30116-16/3X103016F1-9-5s.t.例例2 2问题有无界解问题有无界解2021/9/23236-3 6-3 初始基本可行解初始基本可行解(1)(1)大大M M法法 引入一组人工变量引入一组人工变量,它们在目标函数中的它们在目标函数中的系数均是非常大的正数系数均是非常大的正数M;M;(2)(2)两相法两相法 引入一组人工变量引入一组人工变量,在人工变量未完全离在人工变量未完全离基前目标函数为各人工变量之和基前目标函数为各人工变量之和,当人工变当人工变量完全离基后恢复原目标函数。量完全离基后恢复原目标函数。当当A A内不包含单位矩阵时内不包含单位矩阵时,需引入由人工变量组需引入由人工变量组成的单位矩阵成的单位矩阵,以方便获得初始可行解以方便获得初始可行解.2021/9/2324一一.采用大采用大M M法获得初始基本可行解法获得初始基本可行解s.t.采用大采用大M M法法:s.t.原问题原问题:因因M M是比其他价值系数大得多的正数是比其他价值系数大得多的正数,且人工变量非负且人工变量非负,迭代迭代的结果会使人工变量趋于零的结果会使人工变量趋于零,而获得原问题的基本可行解而获得原问题的基本可行解.2021/9/2325s.t.411MM基变量基变量x1x2x3x4x5Mx42121042Mx53310131X000043F07M4-5M1-4M1-3M表表一一2021/9/2326411MM基变量基变量x1x2x3x4x5Mx42121042Mx53310131X000043F07M4-5M1-4M1-3M411M基变量基变量x1x2x3x4x5Mx40-14/3121.54x1111/3013X110020F14+2M-3+M-4/3-4M/3表表一一表表二二2021/9/2327411基变量基变量x1x2x3x4x51x30-3/411.5-24x115/400.50.4X20.501.5F23.50-13/40表表三三初始基本初始基本可行解可行解411M基变量基变量x1x2x3x4x5Mx40-14/3121.54x1111/3013X110020F14+2M-3+M-4/3-4M/3表表二二2021/9/2328411基变量基变量x1x2x3x4x51x30-3/411.5-24x115/400.50.4X20.501.5F23.50-13/40411基变量基变量x1x2x3x4x51x3011.81x2100.4X300.41.8F32.2表表三三表表四四初始基本初始基本可行解可行解最优解最优解2021/9/2329二二.采用两相法获得初始基本可行解采用两相法获得初始基本可行解大大M M法的法的M M是一个充分大的正数是一个充分大的正数,有时在计算机上不便处理有时在计算机上不便处理.s.t.原问题:原问题:s.t.相相1 1问题:问题:2021/9/233000011基变量基变量x1x2x3x4x51x421210421x53310131X0 00043F07-5-4-30001基变量基变量x1x2x3x4x51x40-14/3121.50 x1111/3013X1 10020F121-4/3表表一一表表二二2021/9/2331000基变量基变量x1x2x3x4x50 x30-3/411.5-20 x115/400.50.4X2 0.501.5F2 00表表三三初始基本初始基本可行解可行解411基变量基变量x1x2x3x4x51x30-3/411.5-24x115/400.50.4X00.501.5F03.5-13/4表表四四初始基本初始基本可行解可行解换换 相相2021/9/2332411基变量基变量x1x2x3x4x51x3011.81x2100.4X100.41.8F12.213/5表表五五最优解最优解2021/9/2333