第4章 线性规划灵敏度分析精选文档.ppt
《第4章 线性规划灵敏度分析精选文档.ppt》由会员分享,可在线阅读,更多相关《第4章 线性规划灵敏度分析精选文档.ppt(95页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第4章 线性规划灵敏度分析1本讲稿第一页,共九十五页 线性规划解除有线性规划解除有唯一最优解唯一最优解的情况外,还有如的情况外,还有如下几种情况下几种情况 无可行解无可行解无可行解无可行解退化退化退化退化无穷多解无穷多解无穷多解无穷多解 无界解无界解无界解无界解人工人工变量变量不能不能从基从基底中底中换出换出基可行基可行解中非解中非零元素零元素个数小个数小于基变于基变量数量数检验数中检验数中零的个数零的个数多于基变多于基变量的个数量的个数检验数大于检验数大于零,但对应零,但对应列元素小于列元素小于等于零,无等于零,无换出变量换出变量2本讲稿第二页,共九十五页唯一最优解唯一最优解 否 否否 是是
2、是添添加加松松弛弛变变量量、人人工工变变量量列出初始单纯形表列出初始单纯形表计算非基变量计算非基变量各列的检验数各列的检验数j所有所有j 0基变量中基变量中有非零的有非零的人工变量人工变量某非基某非基变量检变量检验数为验数为零零无可行解无可行解无穷多最优解无穷多最优解对任一对任一j0有有aik0无界解无界解令令k=maxjxk为换入变量为换入变量对对 所所 有有 aik 0计计 算算i i=bi/aik令令l=min=mini i 第第l l个个基基变变量量为为换换出出变变量,量,alk为主元素为主元素 迭代运算迭代运算.用非基变量用非基变量xk替换换出变量替换换出变量.对主元素行对主元素行(
3、第第l行行)令令bl/alkbl;alj/alkajl对主元素列对主元素列(第第k列列)令令1alk;0;0其它其它元素表中其它行列元素元素表中其它行列元素令令aij-ali/alkaika aijij bi-bl/alkaikbi j-alj/alkk j否对目标函数求极大值标准型线性规划对目标函数求极大值标准型线性规划问题,单纯形法计算步骤的框图:问题,单纯形法计算步骤的框图:3本讲稿第三页,共九十五页第四章第四章 线性规划灵敏度分析线性规划灵敏度分析4.1灵敏度分析的基本原理灵敏度分析的基本原理4.2目标函数系数的灵敏度分析目标函数系数的灵敏度分析4.3右端常数的灵敏度分析右端常数的灵敏
4、度分析4.4技术系数的灵敏度分析技术系数的灵敏度分析*4.5参数线性规划参数线性规划4本讲稿第四页,共九十五页线线性性规规划划的的灵灵敏敏度度分分析析也也称称为为敏敏感感性性分分析析或或优优化化后后分分析析,它它是是研研究究和和分分析析参参数数(cj,bi,aij)的的波波动动对对最最优优解解的的影影响响程程度度,主主要要研研究究下面两个方面:下面两个方面:(1)参参数数在在什什么么范范围围内内变变化化时时,原原最最优优解解或或最最优优基基不不变变数数据的稳定区间;据的稳定区间;(2)当当参参数数超超出出(1)的的变变化化范范围围时时,最最优优解解或或最最优优基基有有何何变变化化如何求出新的最
5、优解和最优基。如何求出新的最优解和最优基。当当模模型型的的参参数数发发生生变变化化后后,可可以以不不必必对对线线性性规规划划问问题题重重新新求求解解,而而用用灵灵敏敏度度分分析析方方法法直直接接在在原原线线性性规规划划取取得得的的最最优优结结果果的的基基础础上上进进行行分分析析或或求求解解,既既可可减减少少计计算算量量,又又可可事事先先知知道道参参数数的的变变化化范范围,及时对原决策作出调整和修正。围,及时对原决策作出调整和修正。4.1灵敏度分析的基本原理灵敏度分析的基本原理5本讲稿第五页,共九十五页单纯形法:单纯形法:对应于基对应于基B的典则形式(典式)的典则形式(典式).Ax=b基变量用非
6、基变量表示:基变量用非基变量表示:代入目标函数:代入目标函数:6本讲稿第六页,共九十五页初始单纯形表初始单纯形表 c c1 c2 cm cm+1 cm+2 cncBxBb x1 x2 xm xm+1 xm+2 xnc1c2cmx1x2xmb1b2bm 100 a1m+1 a1m+2 a1n 010 a2m+1 a2m+2 a2n 001 amm+1 amm+2 amn-z(0)000m+1 m+2 n7本讲稿第七页,共九十五页分析分析 变化对最优解的影响。变化对最优解的影响。CCBCNCBXBbXBXNCBXBB-1 b IB-1 NZ-CB B-1 b0CN-CB B-1 b最优单纯形表最优
7、单纯形表8本讲稿第八页,共九十五页上表中上表中6个常数个常数a1,a2,a3,b,1,2取值在什么范围可使取值在什么范围可使1、现可行解最优,且唯一?何时不唯一?、现可行解最优,且唯一?何时不唯一?2、现基本解不可行;、现基本解不可行;3、问题无可行解;、问题无可行解;4、无有限最优解;、无有限最优解;5、现基本解可行,由、现基本解可行,由x1取代取代x6目标函数可改善。目标函数可改善。cjB-1bcBxBx1x2x3x4x5x6x34a110a20bx4-1-501-102x6a3-300-413j1200-309本讲稿第九页,共九十五页线性规划标准形式线性规划标准形式(1)、参数、参数A,
8、b,C在什么范围内变动,对当前方案无影响?在什么范围内变动,对当前方案无影响?(2)、参数、参数A,b,C中的一个中的一个(几个几个)变动,对当前方案影响?变动,对当前方案影响?(3)、如果最优方案改变,如何用简便方法求新方案?、如果最优方案改变,如何用简便方法求新方案?当线性规划问题中的一个或几个参数变化时,可以用单纯形法从头计当线性规划问题中的一个或几个参数变化时,可以用单纯形法从头计算,看最优解有无变化,但这样做既麻烦又没有必要。算,看最优解有无变化,但这样做既麻烦又没有必要。10本讲稿第十页,共九十五页4.1目标函数系数的灵敏度分析目标函数系数的灵敏度分析考虑检验数考虑检验数(1)若若
9、cj 是非基变量的系数:是非基变量的系数:11本讲稿第十一页,共九十五页解:最优单纯形表解:最优单纯形表 例例1试求试求 c3 在多大范围内变动时,原最优解保持不变。在多大范围内变动时,原最优解保持不变。cj-2-3-400B-1bcBxBx1x2x3x4x5-3x201-1/5-2/51/52/5-2x1107/5-1/5-2/511/5j00-9/5-8/5-1/5-28/512本讲稿第十二页,共九十五页从表中看到从表中看到c3=-4,3=-9/5可得到可得到c3-3=9/5时,即时,即c3-4+9/5=-11/5时原最优解不变。时原最优解不变。解:最优单纯形表解:最优单纯形表 cj-2-
10、3-400B-1bcBxBx1x2x3x4x5-3x201-1/5-2/51/52/5-2x1107/5-1/5-2/511/5j00-9/5-8/5-1/5-28/513本讲稿第十三页,共九十五页(2)若若cj是基变量的系数是基变量的系数14本讲稿第十四页,共九十五页以下分两种情况讨论:以下分两种情况讨论:1如果如果 cr 0上式才有可能不成立上式才有可能不成立,因此有:因此有:2如果如果 cr 0,只有只有arj0,则问题的最,则问题的最优解将发生变化,此时对原最终表适当修改后,应用单优解将发生变化,此时对原最终表适当修改后,应用单纯形法继续计算得到问题的最优解。纯形法继续计算得到问题的最
11、优解。16本讲稿第十六页,共九十五页例例2已知问题的已知问题的最优单纯形表,最优单纯形表,(1)求求c2在什么范围内变动时,在什么范围内变动时,原最优解保持不变;原最优解保持不变;(2)c2=5时,求新的最优解。时,求新的最优解。最优单纯形表最优单纯形表C i23000B-1bCBXBx1x2x3x4x52 x1 1001/4040 x5 00-21/2143 x2 011/2-1/802j00-3/2-1/8014最优解:最优解:最优值:最优值:17本讲稿第十七页,共九十五页C i23000B-1bCBXBx1x2x3x4x52 x1 1001/4040 x5 00-21/2143 x2 0
12、11/2-1/802j00-3/2-1/8014C i23+c2000B-1bCBXBx1x2x3x4x52 x1 1001/4040 x5 00-21/2143+c2 x2 011/2-1/802j00-3/2-c2/2-1/8+c2/8014+2c2(1)求求c2的变动范围,使原最优解保持不变;的变动范围,使原最优解保持不变;c2=c2+c218本讲稿第十八页,共九十五页从表中看到从表中看到可得到可得到-3c21时,原最优解不变。时,原最优解不变。C i23+c2000B-1bCBXBx1x2x3x4x52 x1 1001/4040 x5 00-21/2143+c2 x2 011/2-1/
13、802j00-3/2-c2/2-1/8+c2/8014+2c219本讲稿第十九页,共九十五页C i25000B-1bCBXBx1x2x3x4x52 x1 1001/404160 x5 00-21/21485 x2 011/2-1/802-j00-5/25/8018(2)c2=2,即,即c2由由3变为变为5时,求新的最优解时,求新的最优解C i25000B-1bCBXBx1x2x3x4x52 x1 1010-1/220 x4 00-41285 x2 010-1/8-1/43j00-20-7/419新的最优解:新的最优解:最优值:最优值:20本讲稿第二十页,共九十五页已知最优表如下已知最优表如下例
14、例3某企业利用两种资源生产三种产品的最优计划问题某企业利用两种资源生产三种产品的最优计划问题,归结为下列线性规划归结为下列线性规划cj54300CBXBbx1x2x3x4x545x2x140200110 232 1 11j26000 4 3 121本讲稿第二十一页,共九十五页cj54300CBXBbx1x2x3x4x545x2x140200110 232 1 11j26000 4 3 1最优计划是两种产品分别生产最优计划是两种产品分别生产40单位与单位与20单位,最大产值单位,最大产值z*=260单单位。位。(1)确定)确定x2的价值系数的价值系数c2的变化范围,使原最优解保持最优。的变化范围
15、,使原最优解保持最优。(2)c3在什么范围内变化,最优解不变?在什么范围内变化,最优解不变?(3)若)若c3从从3变为变为9,求新的最优计划。,求新的最优计划。22本讲稿第二十二页,共九十五页解解(1)因为)因为x2为基变量,为基变量,c2的变换范围:的变换范围:因此当因此当x2的价值系数在的价值系数在即在区间即在区间2.5,5变化时,最优解不变。变化时,最优解不变。23本讲稿第二十三页,共九十五页(2)因为)因为c3是是非基变量非基变量x3的价值系数,由公式得到的价值系数,由公式得到 c3+c3的变换范围:的变换范围:即当即当c3在在 c3+c33+4=7时,最优解不变。时,最优解不变。(3
16、)当)当c3从从3变为变为9时,原最优解失去最优性,在时,原最优解失去最优性,在表中作适当变动后(修改表中作适当变动后(修改c3的值,重新计算检验数)的值,重新计算检验数),用单纯形法容易求得新的最优表如下:,用单纯形法容易求得新的最优表如下:24本讲稿第二十四页,共九十五页cj54900CBXBbx1x2x3x4x545x2x140200110 232 1 11j260002 3 149x2x3160/320/32/31/310014/3 1/3 1/31/3j820/3 2/300 7/3 5/3故新的最优解为故新的最优解为X*=(0,160/3,20/3,0,0)最优值最优值z*=820
17、/3,即随着第三种产品价格的上升,总产值,即随着第三种产品价格的上升,总产值上升。上升。25本讲稿第二十五页,共九十五页设分量设分量br变化为变化为br+br,根据前面的讨论:,根据前面的讨论:最优解的基变量最优解的基变量XB=B-1b,那么只要保持,那么只要保持B-1(b+b)0则最优基不变,即基变量保持,只有值的变化;则最优基不变,即基变量保持,只有值的变化;否则,需要利用否则,需要利用对偶单纯形法对偶单纯形法继续计算。继续计算。4.2右端常数的灵敏度分析右端常数的灵敏度分析26本讲稿第二十六页,共九十五页27本讲稿第二十七页,共九十五页例例4已知前述例已知前述例2的最优解及最优单纯形表的
18、最优解及最优单纯形表28本讲稿第二十八页,共九十五页下表为最优单纯形表下表为最优单纯形表C i23000B-1bCBXBx1x2x3x4x52 x1 1001/4040 x5 00-21/2143 x2 011/2-1/802j00-3/2-1/801429本讲稿第二十九页,共九十五页由最优单纯形表得由最优单纯形表得:30本讲稿第三十页,共九十五页31本讲稿第三十一页,共九十五页不可行!用对偶单纯形法计算用对偶单纯形法计算将将b代入原最优单纯形表中,运用对偶单纯形法计算最优解。代入原最优单纯形表中,运用对偶单纯形法计算最优解。32本讲稿第三十二页,共九十五页将将b代入原最优单纯形表中,运用对偶
19、单纯形法计算最优解。代入原最优单纯形表中,运用对偶单纯形法计算最优解。经一次迭代后,求得新的最优解经一次迭代后,求得新的最优解:(4 3 2 0 0)TC i23000B-1bCBXBx1x2x3x4x52 x1 1001/4040 x5 00-21/21-43 x2 011/2-1/804j00-3/2-1/80143/42 x1 1001/4040 x3 001-1/4-1/223 x2 01001/43j000-1/2-3/41733本讲稿第三十三页,共九十五页(1)增加一个新变量增加一个新变量增加一个变新量,相当于系数矩阵增加一列。增加一个变新量,相当于系数矩阵增加一列。设增加变量设增
20、加变量xn+1,则有相应的,则有相应的Pn+1,cn+1。计算出计算出Pn+1=B-1Pn+1,n+1=cn+1-cB Pn+1填入最优单纯形表填入最优单纯形表,若若 n+10则则最优解不变;最优解不变;否则,进一步用单纯形法求解。否则,进一步用单纯形法求解。4.3技术系数的灵敏度分析技术系数的灵敏度分析34本讲稿第三十四页,共九十五页例例5 求当增加求当增加x6,P6=(2,6,3)T,c6=5时,原最优解是否保持时,原最优解是否保持不变,若变动求出新的最不变,若变动求出新的最优解。优解。解解:下表为最优单纯形表下表为最优单纯形表C j23000B-1bcBxBx1x2x3x4x52 x1
21、1001/4040 x5 00-21/2143 x2 011/2-1/802j00-3/2-1/801435本讲稿第三十五页,共九十五页36本讲稿第三十六页,共九十五页用单纯形法进一步求解,可得:用单纯形法进一步求解,可得:x*=(1,1.5,0,0,0,2)T z*=16.5C i230005B-1bCBXBx1x2x3x4x5x62 x1 1001/403/248/30 x5 00-21/212423 x2 011/2-1/801/428j00-3/2-1/805/4142 x1 103/2-1/8-3/4015 x6 00-11/41/2123 x2 013/4-3/16-1/403/2
22、j00-1/4-7/160033/237本讲稿第三十七页,共九十五页(2)增加一个新约束条件增加一个新约束条件增加一个新约束条件相当于在系数矩阵中增加一行。增加一个新约束条件相当于在系数矩阵中增加一行。增加一个约束条件之后,应把最优解带入新的约束,若满足增加一个约束条件之后,应把最优解带入新的约束,若满足则最优解不变,否则填入最优单纯形表作为新的一行,引则最优解不变,否则填入最优单纯形表作为新的一行,引入一个新的非负变量(原约束若是小于等于形式可引入非入一个新的非负变量(原约束若是小于等于形式可引入非负松弛变量,否则引入非负人工变量),并通过矩阵行变负松弛变量,否则引入非负人工变量),并通过矩
23、阵行变换把对应基变量的元素变为换把对应基变量的元素变为0,进一步用单纯形法或对偶单,进一步用单纯形法或对偶单纯形法求解。纯形法求解。38本讲稿第三十八页,共九十五页例例6 求当增加求当增加3x1+2x215时,原时,原最优解是否保持不变,若变最优解是否保持不变,若变动求出新的最优解。动求出新的最优解。解解:下表为最优单纯形表下表为最优单纯形表C i23000B-1bCBXBx1x2x3x4x52 x1 1001/4040 x5 00-21/2143 x2 011/2-1/802j00-3/2-1/801439本讲稿第三十九页,共九十五页将将3x1+2x2+x6=15代入原最优单纯形表。代入原最
24、优单纯形表。经经对偶单纯形法对偶单纯形法迭代一步,可得:迭代一步,可得:最优解为最优解为(3.5,2.25,0,0,3,2)T,最优值为,最优值为13.75C i230000B-1bCBXBx1x2x3x4x5x62 x1 1001/40040 x5 00-21/21043 x2 011/2-1/80020 x632000115j00-3/2-1/800142 x1 1001/40040 x5 00-21/21043 x2 011/2-1/80020 x600-1-1/201-1j00-3/2-1/8001440本讲稿第四十页,共九十五页(3)个别技术系数个别技术系数aij的改变的改变(计划生
25、产的产品工艺结构改变计划生产的产品工艺结构改变)非基变量非基变量xj工艺改变工艺改变只影响单纯形表只影响单纯形表Pj 列列,j.关键看关键看 j 0?还是还是0?.用增加新变量类似方法解决。用增加新变量类似方法解决。基变量基变量xj工艺改变,复杂,根据具体情况讨论。工艺改变,复杂,根据具体情况讨论。41本讲稿第四十一页,共九十五页分析分析:42本讲稿第四十二页,共九十五页注注:Y可由最优单纯形表查得可由最优单纯形表查得43本讲稿第四十三页,共九十五页最优单纯形表为最优单纯形表为求:(求:(1)P3由由(1,3)T改为改为(1,2)T;(2)P1由由(1,2)T改为改为(1,1)T,最优解的变化
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第4章 线性规划灵敏度分析精选文档 线性规划 灵敏度 分析 精选 文档
限制150内