大学精品课件:运筹学(二).ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《大学精品课件:运筹学(二).ppt》由用户(罗嗣辉)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 大学 精品 课件 运筹学
- 资源描述:
-
1、第二章第二章线性规划的线性规划的对偶理论与灵敏度分析对偶理论与灵敏度分析主要内容:主要内容:第一节第一节 线性规划的对偶问题线性规划的对偶问题第二节第二节 对偶问题的基本性质对偶问题的基本性质第三节第三节 影子价格影子价格第四节第四节 对偶单纯形法对偶单纯形法第五节第五节 灵敏度分析灵敏度分析第一节第一节线性规划的对偶问题线性规划的对偶问题一、对偶问题的提出一、对偶问题的提出 例例1 1(回顾第一章)(回顾第一章):美佳公司计划制造:美佳公司计划制造,两种家电产品。已知各制造一件两种家电产品。已知各制造一件时分别占用的设备时分别占用的设备A,BA,B的台时、调试时间、的台时、调试时间、调试工序
2、及每天可用于这两种家电的能调试工序及每天可用于这两种家电的能力、各售出一件的获利情况,如表力、各售出一件的获利情况,如表1-11-1所所示。问该公司应制造两种家电各多少件,示。问该公司应制造两种家电各多少件,使获取的利润为最大。使获取的利润为最大。表表1-1项目项目每天可用能力每天可用能力设备设备A(h)0515设备设备B(h)6224调试工序(调试工序(h)115利润(元)利润(元)21 设设x x1 1,x,x2 2 分别代表分别代表,两种家电两种家电的生产量,的生产量,此问题的数学模型为:此问题的数学模型为:目标函数目标函数 约束条件约束条件212maxxxz 0,52426155.21
3、21212xxxxxxxst 例例2 2:在例:在例1 1中,若有某个公司想把美佳公司的中,若有某个公司想把美佳公司的资源收购过来,它至少付出多大的代价,才能使资源收购过来,它至少付出多大的代价,才能使美佳公司愿意放弃生产活动,出让自己的资源美佳公司愿意放弃生产活动,出让自己的资源。但少生产一件但少生产一件产品,则产品,则丧失了丧失了1元的利润元的利润少生产一件少生产一件产品,则可以节省设产品,则可以节省设备备A、设备、设备B和调试工序和调试工序5、2、1个小时,把这些资源出租,每小时个小时,把这些资源出租,每小时租金分别为租金分别为y1,y2,y3,就可以获就可以获得租金得租金5y1+2y2
4、+y3?但少生产一件但少生产一件I产品,则产品,则丧失了丧失了2元的利润元的利润同样,少生产一件同样,少生产一件I产品,则可以产品,则可以节省设备节省设备A、设备、设备B和调试工序和调试工序0、6、1个小时,把这些资源出租,个小时,把这些资源出租,就可以获得租金就可以获得租金0y1+6y2+y3 所以,只有当所以,只有当 时,才可能出让。时,才可能出让。总的出让费总的出让费 ,最低出让费即为:最低出让费即为:2632 yy125321 yyy32152415minyyyw 32152415yyy 因此,该问题的数学模型为:因此,该问题的数学模型为:32152415minyyyw 2632 yy
5、125321 yyy0,321 yyy.t s二、对称形式下对偶问题的一般形式二、对称形式下对偶问题的一般形式 满足下列条件的线性规划问题称为具满足下列条件的线性规划问题称为具有有对称形式对称形式:其变量均具有非负约束,其约束条件其变量均具有非负约束,其约束条件当目标函数求极大时均取当目标函数求极大时均取“”,当目,当目标函数求极小时均取标函数求极小时均取“”。,n),(jxbxaxaxa bxaxaxabxaxaxas.t.jmnmnmmnnnn21022112222212111212111nnxcxcxcZ 2211maxmmybybybW2211min ,m),(iycyayaya cy
6、ayayacyayayas.t.inmmnnnmmmm21022112222211211221111原问题原问题对偶问题对偶问题CXZ maxbYWmin0bXAXts.0YCYAts.原问题原问题对偶问题对偶问题A A约束系数矩阵约束系数矩阵约束系数矩阵的转置约束系数矩阵的转置b b约束条件右端向量约束条件右端向量目标函数中的价值系数向量目标函数中的价值系数向量C C目标函数中的价值系数向量目标函数中的价值系数向量约束条件右端向量约束条件右端向量目标函数目标函数约束条件约束条件决策变量决策变量CXz maxbYwminbAXCYA0X0Y对称形式下原问题与对偶问题的对应关系:对称形式下原问题
7、与对偶问题的对应关系:即:即:对偶问题对偶问题212maxxxz 0,5 2426155 .2121212xxxxxxxst32152415minyyyw 2632 yy125321 yyy0,321 yyy.t s原问题原问题对偶问题对偶问题如:如:三、非对称形式原三、非对称形式原-对偶问题的关系对偶问题的关系332211maxxcxcxcz 约约束束无321333323213123232221211313212111,0,0.xxxbxaxaxabxaxaxabxaxaxast求下述线性规划问题的对偶问题:求下述线性规划问题的对偶问题::则则0,0,x x,x x,x x,x xx xx
8、x,x xx x令令3 33 32 23 33 33 32 22 2 33332211maxxcxcxcxcz 03213321333333323213123233232221211313313212111x,x,x,x)(bxaxaxaxa)(bxaxaxaxa)(bxaxaxaxast.约束(约束(2-2)和约束()和约束(3)两边同乘以)两边同乘以“-1”,则变换为:,则变换为:约束(约束(2)可以用以下两个约束来表示:)可以用以下两个约束来表示:2)-(2 1)-(2 23233232221212323323222121bxaxaxaxabxaxaxaxa 令各约束对应的对偶变量为:令
9、各约束对应的对偶变量为:,则其对偶问题为:,则其对偶问题为:33222211minybybybybw 0,.32213333223223113333322322311323322222221121331221221111yyyycyayayayacyayayayacyayayayacyayayayast3221,yyyy 33332211maxxcxcxcxcz 0,.33213333333232131232332322212123233232221211313313212111xxxxbxaxaxaxabxaxaxaxabxaxaxaxabxaxaxaxast,则有:,则有:令令33222,
10、yyyyy 332211minybybybw 0,0 .3213333223113333322311323322221121331221111yyycyayayacyayayacyayayacyayayast无无约约束束整理,得:332211minybybybw 0,0.321333322311323322221121331221111yyycyayayacyayayacyayayast无无约约束束因此,线性规划原问题与对偶问题的关系见下表:因此,线性规划原问题与对偶问题的关系见下表:原问题(或对偶问题)原问题(或对偶问题)对偶问题(或原问题)对偶问题(或原问题)目标函数目标函数 max目标函
11、数目标函数 min约约束束条条件件m个个m个个变变量量0000=无约束无约束变变量量n个个n个个约约束束条条件件0000无约束无约束=约束条件右端项约束条件右端项目标函数变量的系数目标函数变量的系数目标函数变量的系数目标函数变量的系数约束条件右端项约束条件右端项正常对正常,反常对反常正常对正常,反常对反常 在正常的情况下,线性规划的所有变量都是在正常的情况下,线性规划的所有变量都是0的,的,求最大化问题的约束是求最大化问题的约束是型的,而求最小化问题的约束是型的,而求最小化问题的约束是型的。型的。所谓所谓“正常对正常,反常对反常正常对正常,反常对反常”,是指,是指:只要原模型中的只要原模型中的
12、变量和约束变量和约束都符合这些规律,对偶模都符合这些规律,对偶模型中对应的型中对应的约束和变量约束和变量也一定符合这些规律;也一定符合这些规律;如果原模型中的变量和约束违反了这些规律,则对偶如果原模型中的变量和约束违反了这些规律,则对偶模型中的约束和变量也一定违反这些规律。模型中的约束和变量也一定违反这些规律。记忆口诀记忆口诀321362minxxxW无约束321321321321,0,060434803x 430 xxxxxxxxxxx例例3:写出下列线性规划问题的对偶问题:写出下列线性规划问题的对偶问题:正常正常正常正常反常反常反常反常321608030maxyyyZ0,034 363 4
13、24 321321321321yyyyyyyyyyyy无约束该问题的对偶问题为:该问题的对偶问题为:正常正常反常反常反常反常正常正常第二节第二节对偶问题的基本性质对偶问题的基本性质一、单纯形法的矩阵描述一、单纯形法的矩阵描述0bCXXAXtsz.max 考虑线性规划问题(有考虑线性规划问题(有n个变量,个变量,m个约束条件):个约束条件):引入松弛变量引入松弛变量TmnnnSxxxX),(21将其转化为标准形式,得到:将其转化为标准形式,得到:0,.0max SSSXXIXAXtsXz0bCX其中其中I为为mm阶单位矩阵。阶单位矩阵。设设A中存在一个可行基中存在一个可行基B,则变量,则变量X可
14、分为基变量可分为基变量XB和和非基变量非基变量XN,相应地,矩阵,相应地,矩阵A可以分为可以分为B和和N,价值系,价值系数数C可以分为可以分为CB和和CN,即:,即:),(NBXXX),(NBA),(NBCCC 则上述线性规划问题可写成如下形式:则上述线性规划问题可写成如下形式:0,.0max SSNBSNNBBXXIXNXBXtsXz0bXCXC列出初始单纯形表:列出初始单纯形表:CBCNC0BC0BxbSXbBXNXSXBNIzc BCNC0因为因为B为一个可行基,因此,可以通过迭代(上表中第三行为一个可行基,因此,可以通过迭代(上表中第三行左乘左乘B-1)得到另外一个单纯形表)得到另外一
15、个单纯形表:CBCNC0BCBCBxbBXbB1BXNXSXBB1NB11Bzc BBCCBB1NBCCBN11BCBCBCNC0BCBCBxbBXbB1BXNXSXBB1NB11Bzc BBCCBB1NBCCBN11BCB若若XB为最优基变量,则对应的目标函数值为为最优基变量,则对应的目标函数值为:bBCXzBSNNBB10 XCXC且对于上表中各检验数,有:且对于上表中各检验数,有:把上面三个式子的前两个式子合并,得到:把上面三个式子的前两个式子合并,得到:01 BBCCBB01 NBCCBN01 BCB01ABCCB01BCBbYWmin0YCYAts.可见,当原问题得到最优解时,其松弛
16、变量检验数的相反数可见,当原问题得到最优解时,其松弛变量检验数的相反数 是该问题的对偶问题的一个可行解。是该问题的对偶问题的一个可行解。该线性规划问题的对偶问题为:该线性规划问题的对偶问题为:1BCBCABCB101BCB212maxxxz0,5 2426155 .2121212xxxxxxxst32152415minyyyw2632 yy125321yyy0,321yyy.ts原问题原问题对偶问题对偶问题例:例:分别引入松弛变量和剩余变量,化为标准形式:分别引入松弛变量和剩余变量,化为标准形式:543210002maxxxxxxz0,5 2426155 .5432152142132xxxxx
17、xxxxxxxxst543210052415minyyyyyw26432yyy1255321yyyy0,54321yyyyy.ts原问题原问题对偶问题对偶问题0001210011500102624000 15015000012jjzcjcBCBxb4x2x3x4x5x3x1x5x021000015/20015/4-15/227/21001/4-1/213/2010-1/43/2000-1/4-1/2jjzc jcBCBxb1x2x3x4x5x3x1x2x原问题初始单纯形表:原问题初始单纯形表:原问题最终单纯形表:原问题最终单纯形表:110260501B6102103054*B4B2/34/10
18、2/14/102/154/51*1BBB2/32/72/152/1562/56052/15244/5151524152/34/102/14/102/154/511bB21000015/20015/4-15/227/21001/4-1/213/2010-1/43/2000-1/4-1/2jjzc jcBCBxb1x2x3x4x5x3x1x2x原问题最终单纯形表:原问题最终单纯形表:对偶问题最终单纯形表:对偶问题最终单纯形表:2100001/4-5/410-1/41/411/215/2011/2-3/215/2007/23/2jjzc jcBCBxb3y2y3y4y5y2y1y1.1.对称性对称性
19、2.2.弱对偶性弱对偶性3.3.最优性最优性4.4.强对偶性强对偶性5.5.互补松弛性互补松弛性二、对偶问题的基本性质二、对偶问题的基本性质对称性:对称性:对偶问题的问偶即原问题对偶问题的问偶即原问题证明证明:CXZ maxbYWmin0bXAXts.0YCYAts.对偶问题对偶问题原问题原问题对偶问题中若令对偶问题中若令w=-w,则该对偶问题可改写为:则该对偶问题可改写为:bYWmax0YCYAts.(1)bYWmax0YCYAts.CXzmin0XbAXts.对偶问题对偶问题原问题原问题CXz max0XbAXts.zz令即为原问题(即为原问题(1)把对偶问题作为原问题,写出其对偶问题如下
20、:把对偶问题作为原问题,写出其对偶问题如下:弱对偶性:弱对偶性:如果如果 是原问题的可行解,是原问题的可行解,是其对偶问题的可行是其对偶问题的可行解,则恒有:解,则恒有:XYbYXCCXZ maxbYWmin0bXAXts.0YCYAts.对偶问题对偶问题证明证明:因为因为 是原问题的可行解,故有:是原问题的可行解,故有:XbXA式子两边左乘式子两边左乘 (是对偶问题的可行解,因此是对偶问题的可行解,因此 ),得到:),得到:YY0YbYXAY 同理,因为同理,因为 是对偶问题的可行解,故有:是对偶问题的可行解,故有:YCYA式子两边右乘式子两边右乘 (是原问题的可行解,因此是原问题的可行解,
21、因此 ),得到:),得到:X0XXCXAYXCAY(1)(1)(2)(2)综合式子(综合式子(1)和()和(2),得到:),得到:bYXAYXC推论:推论:(1)原问题任一可行解的目标函数值是其对偶问题目标)原问题任一可行解的目标函数值是其对偶问题目标函数值的下界;反之对偶问题任一可行解的目标函数值函数值的下界;反之对偶问题任一可行解的目标函数值是其原问题目标函数值的上界。是其原问题目标函数值的上界。(2)在一对对偶问题中,若其中一个问题可行但目标函)在一对对偶问题中,若其中一个问题可行但目标函数无界,则另一个问题不可行;反之不成立。数无界,则另一个问题不可行;反之不成立。(3)在一对对偶问题
22、中,若一个问题可行,而另一个问)在一对对偶问题中,若一个问题可行,而另一个问题不可行,则该可行的问题的目标函数无界。题不可行,则该可行的问题的目标函数无界。因此因此 是原问题的最优解。是原问题的最优解。最优性:最优性:若若 和和 分别是原问题和对偶问题的可行解分别是原问题和对偶问题的可行解,且且 ,则则 和和 分别是原问题和对偶问题的最优解。分别是原问题和对偶问题的最优解。*X*YbYCX*X*Y证明证明:设:设 是原问题的任一可行解,由弱对偶性的推是原问题的任一可行解,由弱对偶性的推论论1,可得:,可得:XbYXC*又因为:又因为:,故有:,故有:bYCX*CXXC*X同理,可证同理,可证
23、是对偶问题最优解。是对偶问题最优解。*Y强对偶性:强对偶性:若一对对偶问题的原问题和对偶问题都有可行解,则它若一对对偶问题的原问题和对偶问题都有可行解,则它们都有最优解,且目标函数的最优值相等。们都有最优解,且目标函数的最优值相等。设设 是原问题的最优解,对应的最优基是是原问题的最优解,对应的最优基是B,我们已经证明,我们已经证明过过 是对偶问题的可行解。且对应的目标函数都为是对偶问题的可行解。且对应的目标函数都为 。因此,由最优性可知,因此,由最优性可知,是对偶问题的最优解。是对偶问题的最优解。结论得证。结论得证。*X1BCBbBCB11BCB证明证明:因为原问题和对偶问题都有可行解,由弱对
24、偶性推论:因为原问题和对偶问题都有可行解,由弱对偶性推论1可知,原问题和对偶问题都有界,即它们一定存在最优解。可知,原问题和对偶问题都有界,即它们一定存在最优解。综上所述,一对对偶问题的解,只能有下面三综上所述,一对对偶问题的解,只能有下面三种情况之一出现:种情况之一出现:(1)都有最优解,分别设为)都有最优解,分别设为 和和 ,且必,且必有有 ;(2)一个问题无界,另一个问题无可行解;)一个问题无界,另一个问题无可行解;(3)两个都无可行解。)两个都无可行解。*X*YbYCX*互补松弛性:互补松弛性:若若 和和 分别是原问题和对偶问题的可行解分别是原问题和对偶问题的可行解,为原问题的松弛变量
25、的值,为原问题的松弛变量的值,为对偶问题剩余变量为对偶问题剩余变量的值。的值。和和 分别是原问题和对偶问题最优解的分别是原问题和对偶问题最优解的充分必要条件是充分必要条件是 。*X*YSYSX0*XYXYSS*X*Y 若若 和和 分别是原问题和对偶问题的可行解。则有:分别是原问题和对偶问题的可行解。则有:左乘(左乘(1)证明:证明:(1)(1)充分性的证明充分性的证明 由上述所设,可知:由上述所设,可知:(1)*bIXAXS(2)*CIYAYS*X*Y(3)*bYXYAXYS(4)*CXXYAXYS*Y右乘(右乘(2)*X(3)-(4),得到:),得到:*CXbYXYXYSS若若 0*XYXY
展开阅读全文