第6章单纯形法的灵敏度分析与精选文档.ppt
《第6章单纯形法的灵敏度分析与精选文档.ppt》由会员分享,可在线阅读,更多相关《第6章单纯形法的灵敏度分析与精选文档.ppt(61页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第6章单纯形法的灵敏度分析与1本讲稿第一页,共六十一页上节回顾上节回顾4.退化问题(求出基变量过程中存在两个以上的最小比值,这样在下一次迭代就有一个或多个基变量等于零,造成的后果是:最优值在经过了一次或几次迭代而没有改善,降低了单纯型算法的效率)5.针对退化问题,讲解了一个迭代循环的例子,引出了勃兰特法则:按照决策变量、松弛变量、人工变量顺序排序。不管是检验数相等,还是比值相等,都选下标最小的为入基变量或出基变量。时间问题,习题没有做。2本讲稿第二页,共六十一页第六章第六章 单纯形法的灵敏度分析与对偶问题单纯形法的灵敏度分析与对偶问题11单纯形表的灵敏度分析单纯形表的灵敏度分析22线性规划的对
2、偶问题线性规划的对偶问题33对偶规划的基本性质对偶规划的基本性质44对偶单纯形法对偶单纯形法3本讲稿第三页,共六十一页单纯形表单纯形表4本讲稿第四页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析一、目标函数中变量系数Ck灵敏度分析(在什么范围内变化,最优解不变,与第二章,第三章联系起来)在线性规划的求解过程中,目标函数系数的变动将会影响检验数的取值,但是,当目标函数的系数的变动不破坏最优判别准则时,原最优解不变,否则,原最优解将发生变化,要设法求出新的最优解。下面我们具体的分析 1.在最终的单纯形表里,X k是非基变量 由于约束方程系数增广矩阵在迭代中只是其本身的行的初等变换与Ck没
3、有任何关系,所以当Ck变成Ck+Ck时,在最终单纯形表中其系数的增广矩阵不变,又因为Xk是非基变量,所以基变量的目标函数的系数不变,即CB不变,可知Zk也不变,只是Ck变成了Ck+Ck。这时 K=Ck-Zk就变成了Ck+Ck-Zk=K+Ck。要使原来的最优解仍为最优解,只要 K+Ck0即可,也就是Ck的增量 Ck-K。5本讲稿第五页,共六十一页2.在最终的单纯形表中,X k是基变量 当Ck变成Ck+Ck时,最终单纯形表中约束方程的增广矩阵不变,但是基变量的目标函数的系数CB变了,则ZJ(J=1,2,.,N)一般也变了,不妨设CB=(CB1,CB2。,Ck,,CBm),当CB变成=(CB1,CB
4、2。,Ck+Ck,CBm),则:ZJ=(CB1,CB2。,Ck,,CBm)(a1j,a2j,aKj ,amj)ZJ=(CB1,CB2。,Ck+Ck,,CBm)(a1j,a2j,aKj ,amj)=ZJ+Ck aKj 6本讲稿第六页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析根据上式可知 检验数 J(J=1,2,.,M)变成了 J,只要当J K时,有 J=CJ-ZJ=J-CK aKj。要使最优解不变,J 0,就可求出 的取值范围,也就是使得第K个约束条件的对偶价格不变的bk的变化范围。,21本讲稿第二十一页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析(k的取值固定的值,i
5、的取值是1到m)下面我们仍以第二章例1在最终单纯形表上对b1进行灵敏度分析。最终单纯形表如下所示:迭代次数基变量CBX1 X2 S1 S2 S3b50 100 0 0 02X1501 0 1 0 -150 S200 0 -2 1 150 X21000 1 0 0 1250 ZJ50 100 50 0 5027500CJ-ZJ0 0 -50 0 -5022本讲稿第二十二页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析 我们对b1进行灵敏度分析,因为在第一个约束方程中含有松弛变量S1,实际意义可以描述为:设备台时数在250与325之间变化,则设备台时数的对偶价格不变,都为每台设备台时50
6、元。23本讲稿第二十三页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析三、约束方程系数矩阵三、约束方程系数矩阵A A灵敏度分析(约束条件改变,对最优解的影响,比如有新产品要生灵敏度分析(约束条件改变,对最优解的影响,比如有新产品要生产、生产工艺改进等)产、生产工艺改进等)下面分两种情况讨论 1.在最终单纯形表上在最终单纯形表上Xk是非基变量:在初始单纯形表上的变量是非基变量:在初始单纯形表上的变量Xk的系数列的系数列Pk改变为改变为Pk经过迭代后,在最终单纯形表上经过迭代后,在最终单纯形表上Xk是非基变量。是非基变量。由于单纯形表的迭代是约束方程的增广矩阵的行变换,Pk变成Pk仅仅影
7、响最终单纯形表上第k列数据,包括Xk的系数列、Zk以及 k,这时最终单纯形表上的Xk的系数列就变成了B-1Pj,而Zk就变成CBB-1Pk,新的检验数 k=Ck-CBB-1Pk。若 k0,则原最优解仍然为最优解。若 k 0,则继续进行迭代以求出最优。例例 以第二章例1为基础,设该厂除了生产,种产品外,现在试制成一个新产品,已知生产产品,每件需要设备2台时,并消耗A原料0.5公斤。B原料1.5公斤,获利150元,问该厂应该生产该产品多少?解:这是一个增加新变量的问题。我们可以把它认为是一个改变变量X3在初始表上的系数列的问题,24本讲稿第二十四页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏
8、度分析迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 150X1501 0 1 0 -1 0.550 S200 0 -2 1 1 -250 X21000 1 0 0 1 1.5250 ZJ50 100 50 0 50 17527500CJ-ZJ0 0 -50 0 -50 -2525本讲稿第二十五页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析例 假设上例题中产品的工艺结构有了改进,这时生产1件产品需要使用1.5台设备,消耗原料A为2千克,原料B为1千克,每件产品的利润为160元,问该厂的生产计划是否要修改。解:首先求出X3在最终表上的系数列 迭代次数
9、基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 1602X1501 0 1 0 -1 0.55050/0.5 S200 0 -2 1 1 050 X21000 1 0 0 1 1250250/1 ZJ50 100 50 0 50 12527500CJ-ZJ0 0 -50 0 -50 3526本讲稿第二十六页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析接下来又可以有新的迭代S3进基,迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 1503X31602 0 2 0 -2 1100-S200 0 -2 1 1 05050/1 X
10、2100-20 1 -2 0 3 0150250/3 ZJ120 100 120 0 -20 16031000CJ-ZJ-70 0 -120 0 20 027本讲稿第二十七页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析接上页 可知此规划的最优解X1=0,X2=0,S1=0,S2=0,S3=50,X3=200,此时,最大目标函数为32000元。也就是说,该厂的新的生产计划为不生产、产品,生产产品200件,可获得最大利润32000元。迭代次数基变量CBX1 X2 S1 S2 S3 X3 b50 100 0 0 0 1504X31602 0 2 0 -2 1200-S300 0 -2 1
11、 1 05050/1 X2100-2 1 4 -3 0 00250/3 ZJ120 100 80 20 0 16032000CJ-ZJ-70 0 -80 -20 0 028本讲稿第二十八页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析 2.最终表上最终表上XK是基变量:在初始表上的变量是基变量:在初始表上的变量XK的系数的系数PK改变改变为为PK,经过迭代后,在最终表上,经过迭代后,在最终表上XK是基变量,(比如是基变量,(比如X1,X2,制作产品,制作产品1,产品,产品2的工艺有改善,在这种情况下原最优解的工艺有改善,在这种情况下原最优解的可行性和最优解都可能被破坏,问题十分复杂)
12、一般不去的可行性和最优解都可能被破坏,问题十分复杂)一般不去修改最终单纯形表,而是重新用单纯形法计算。修改最终单纯形表,而是重新用单纯形法计算。29本讲稿第二十九页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析四、增加一个约束条件的灵敏度分析(增加一个约束条件对最优解的影响)在原线性规划中增加一个约束条件时,(1)先将原问题的最优解的变量值代入新增的约束条件,如满足则说明新增的条件没有起到限制作用,故原最优解不变,停止计算。(2)否则将新增的约束添入原最终单纯形表上,进一步求解。下面仍以第三章例1为例来加以说明。例:假如该工厂除了在设备台时,原材料A、B上对该厂的生产有限制外,还有电
13、力供应上的限制。最高供应电量为5000度,而生产一个产品需要用电10度,而生产一个产品需要用电30度。试分析此时该厂获得最大利润的生产计划?30本讲稿第三十页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析 解:先将原问题的最优解=50=50,=250代入用电量的约束条件得:1050+30250=500+75005000,所以原题的最优解不是本题的最优解。在用电量的约束条件中加入松驰变量S4后得:把这个约束条件加入到原最终单纯形表上,其中S4为基变量,得表如下:迭代次数基变量b比值501000000501010-1050000-21105010001001025001030000150
14、00501005005002750000-500-50031本讲稿第三十一页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析 在上表中的X1,X2不是单位向量,故进行行的线性变换,得迭代次数基变量CBx1x2s1s2s3s4b比值501000000 x1501010-1050s2000-211050 x2100010010250s4000-100-201-3000zj501005005002750000-500-500把上表中的S4行的约束可以写为:上式两边乘以(-1),再加上人工变量a1得:用上式替换上表中的S4行,得下表:32本讲稿第三十二页,共六十一页11单纯形表的灵敏度分析单纯
15、形表的灵敏度分析迭代次数基变量x1x2s1s2s3s4a1b比值501000000-Mx1501010-10050s2000-21(1)0050 x21000100100250s4-M00-100-20113000zj5010050-10M050-20M0-M0010M-50020M-5000 x15010-11000100s3000-2110050 x2100012-1000200s4-M0050-200-112000zj50100150-50M20M-500M-M050M-15050-20M0-M0 x1501003/50-1/501/50140s300001/51-2/502/50130
16、 x2100010-1/502/50-2/50120s40001-2/50-1/501/5040zj5010001003-300-100-3-M+333本讲稿第三十三页,共六十一页11单纯形表的灵敏度分析单纯形表的灵敏度分析 由上表可知,最优解为:即该工厂在添加了用电限量以后的最优生产计划为产品生产140件,产品生产120件。34本讲稿第三十四页,共六十一页上节回顾上节回顾单纯形法的灵敏度分析1.目标函数系数Cj的灵敏度2.常数项bj的灵敏度分析3.系数矩阵A的灵敏度分析 在初始表上的变量XK的系数PK改变为PK,经过迭代后,在最终表上XK是基变量,(比如X1,X2,制作产品1,产品2的工艺有
17、改善,在这种情况下原最优解的可行性和最优解都可能被破坏,问题十分复杂)一般不去修改最终单纯形表,而是重新用单纯形法计算。4.增加一个约束条件的灵敏度分析(1)先将原问题的最优解的变量值代入新增的约束条件,如满足则说明新增的条件没有起到限制作用,故原最优解不变,停止计算。(2)否则将新增的约束添入原最终单纯形表上,进一步求解。35本讲稿第三十五页,共六十一页对偶单纯形法:对偶单纯形法也是解决线性规划问题的一种方法。对偶单纯形法是在保持原有问题的所有检验数都小于0的情况下,通过迭代使得所有的约束都大于等于0,最后求得最优解。优点:简化计算是对偶单纯形法的优点 例子:缺点:但是它在使用上有很大的局限
18、,这主要是 大多数线性规划问题很难找到初始解使得其所有检验数都小于0。适用目标函数的特点:A.max f=-M1*X1-M2*X2.的形式 B.用在灵敏度分析中(增加约束条件的)。上节回顾上节回顾36本讲稿第三十六页,共六十一页1.对偶问题的定义2.原问题与对偶问题的关系3.如何根据原问题的结果去找对偶问题的答案4.如何将原问题转化成对偶问题5.对偶规划的基本性质2 2 线性规划的对偶问题线性规划的对偶问题37本讲稿第三十七页,共六十一页任何一个求极大化的线性规划问题都有一个求极小化的线性规划问题与之对应,反之亦然,如果我们把任何一个求极大化的线性规划问题都有一个求极小化的线性规划问题与之对应
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 单纯 灵敏度 分析 精选 文档
限制150内