大学精品课件:第三章 对偶理论和灵敏度分析(第1-5节).ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《大学精品课件:第三章 对偶理论和灵敏度分析(第1-5节).ppt》由用户(罗嗣辉)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 大学精品课件:第三章 对偶理论和灵敏度分析第1-5节 大学 精品 课件 第三 对偶 理论 灵敏度 分析
- 资源描述:
-
1、第第1页页第第2页页第第3页页掌握用矩阵描述单纯形法计算过程的作用:掌握用矩阵描述单纯形法计算过程的作用:有助于加深对单纯形法的理解。有助于加深对单纯形法的理解。有助于下面对偶理论和灵敏度分析的学习。有助于下面对偶理论和灵敏度分析的学习。第第4页页0 .max BNBNNBBXbNXBXtsXCXCz从而线性规划问题的标准型从而线性规划问题的标准型0 max XbAXCXz第第5页页bNXBXNB 当前单纯形表中当前单纯形表中的右端常数的右端常数当前单纯形表中的非基当前单纯形表中的非基变量系数列向量变量系数列向量NBNXbBX NBNXBbBX11 当前单纯形表中的基变量系数列向量总是单位向量
2、当前单纯形表中的基变量系数列向量总是单位向量第第6页页max z=2x1+3x2 0,12 4 16 482 .212121xxxxxxts例:例:解:解:0,12 4 16 48 2 .543215241321xxxxxxxxxxxxtsmax z=2x1+3x2第第7页页 410004201B 08/12/112/1204/101B已知最终单纯形表中的基变量为已知最终单纯形表中的基变量为 xB=(x1,x5,x2),则可以则可以直接直接计算出最终单纯形表中的系数列向量:计算出最终单纯形表中的系数列向量:当前单纯形表中的非基变量系数列向量当前单纯形表中的非基变量系数列向量 8/12/12/1
3、24/1000100108/12/112/1204/101NB第第8页页 410004201B 08/12/112/1204/101B已知当前单纯形表中的基变量为已知当前单纯形表中的基变量为 xB=(x1,x5,x2),则可以则可以直接直接计算出当前单纯形表中的基变量的值:计算出当前单纯形表中的基变量的值:2441216808/12/112/1204/101bB当前单纯形表中的基变量的值(即右端常数)当前单纯形表中的基变量的值(即右端常数)第第9页页cj23000iCBXBbx1x2x3x4x52x141001/400 x5400-21/213x22011/2-1/80最终单纯形表的形式为:最
4、终单纯形表的形式为:第第10页页NNBBXCXCz 将将 XB 表达式带入目标函数:表达式带入目标函数:NBNXBbBX11 NNNBXCNXBbBC )(11NBNBXNBCCbBC)(11 最优目标函数值最优目标函数值第第11页页NBCCBN1 非基变量检验数非基变量检验数BBCCBB1 基变量检验数基变量检验数)0(ABCCB1 所有检验数:所有检验数:第第12页页 410004201B 08/12/112/1204/101B已知最终单纯形表中的基变量为已知最终单纯形表中的基变量为 xB=(x1,x5,x2),则可以则可以直接直接计算出最终单纯形表中的检验数:计算出最终单纯形表中的检验数
5、:)08/12/300(10040010040012108/12/112/1204/10)302()00032(1 ABCCB第第13页页cj23000iCBXBbx1x2x3x4x52x141001/400 x5400-21/213x22011/2-1/80c j z j00-3/2-1/8014最终最终单纯形表:单纯形表:bBCzB1 141216808/12/112/1204/10)302(第第14页页总总 结结从从任一单纯形表(不一定是初始单纯形表)任一单纯形表(不一定是初始单纯形表)出发,出发,只要能够知道当前单纯形表中的哪几个变量是基变只要能够知道当前单纯形表中的哪几个变量是基变量
6、,就可以直接计算出量,就可以直接计算出当前单纯形表当前单纯形表。ABA1 bBb1 ABCCB1 (1)约束条件系数矩阵:)约束条件系数矩阵:(2)约束条件右端常数:)约束条件右端常数:(3)检验数:)检验数:第第15页页只要能够知道当前单纯形表中的哪几个变量是基变只要能够知道当前单纯形表中的哪几个变量是基变量,就可以由初始单纯形表直接计算出当前单纯形量,就可以由初始单纯形表直接计算出当前单纯形表。表。要确定当前单纯形表中哪几个变量是基变量,必须要确定当前单纯形表中哪几个变量是基变量,必须要解决的问题:要解决的问题:(1)换入变量的确定。)换入变量的确定。(2)换出变量的确定。)换出变量的确定
7、。第第16页页(1)换入变量的确定)换入变量的确定ABCCB1 计算所有变量检验数:计算所有变量检验数:lklikikiiPBbBPBPBbB)()()0)()()(min11111 (2)换出变量的确定)换出变量的确定换入变量:换入变量:x k换出变量:换出变量:x l第第17页页基变量基变量XB非基变量非基变量 XN等式右边等式右边系数矩阵系数矩阵检验数检验数0IBB 1NB1 bB1 NBCCBN1 第第18页页由此可以看出,由此可以看出,B-1 的计算是整个问题的关键的计算是整个问题的关键第第19页页max z=2x1+3x2 0,12 4 16 482 .212121xxxxxxts
8、例:例:第第20页页cj23000iCBXBbx1x2x3x4x50 x381210040 x41640010-0 x5120 4 0013c j z j23000解解:(:(1)从初始单纯形表出发。)从初始单纯形表出发。0 x32 1 010-1/220 x4164001043x2301001/4-c j z j2000-3/41211400010201 B)(243xxx第第21页页cj23000iCBXBbx1x2x3x4x52x121010-1/2-0 x4800-41 2 43x2301001/412c j z j00-201/43411400014201 B11410004201
9、B2x141001/400 x5400-21/213x22011/2-1/80c j z j00-3/2-1/80)(241xxx)(251xxx第第22页页cj23000iCBXBbx1x2x3x4x50 x321010-1/220 x4164001043x2301001/4-c j z j2000-3/42x141001/40-0 x5400-21/213x22011/2-1/80c j z j00-3/2-1/8024(2)从任一单纯形表(非初始单纯形表)出发。)从任一单纯形表(非初始单纯形表)出发。1114/1000402/11 B 251xxx第第23页页1114/1000402/1
10、1 B 18/12/102/1204/10第第24页页 0 ,1 232 411 2 3min32131321321321xxxxxxxxxxxxxxz例:例:第第25页页 0,1 23 2 411 2 003 max76543217316532143217654321xxxxxxxxxxxxxxxxxxxMxMxxxxxxz解:加入解:加入 松弛变量:松弛变量:x4 剩余变量:剩余变量:x5 人工变量:人工变量:x6 和和 x7 得到:得到:第第26页页cj3-1-100-M-MiCBXBbx1x2x3x4x5x6x70 x4111-21100011-Mx63-4120-1103/2-Mx7
11、1-20 1 00011c j z j3-6MM-13M-10-M000 x4103-20100-1-Mx610 1 00-11-21-1x31-2010001-c j z j1M-100-M01-3M第第27页页cj3-1-100-M-MiCBXBbx1x2x3x4x5x6x70 x412 3 001-22-54-1x210100-11-2-1x31-2010001-c j z j1000-11-M-1-M3x141001/3-2/32/3-5/3-1x210100-11-2-1x390012/3-4/34/3-7/3c j z j000-1/3-1/31/3-M2/3-M第第28页页由此可
12、以看出,当前单纯形表中基变量所对应的由此可以看出,当前单纯形表中基变量所对应的B-1 等于在出发单纯形表中,系数列向量为单位向量的变等于在出发单纯形表中,系数列向量为单位向量的变量在当前单纯形表中所对应的列向量组成的矩阵。量在当前单纯形表中所对应的列向量组成的矩阵。第第29页页例例 已知某最大化线性规划问题用单纯形法求解时的初已知某最大化线性规划问题用单纯形法求解时的初始单纯形表及最终单纯形表如下表所示,求表中各括始单纯形表及最终单纯形表如下表所示,求表中各括号内未知数的值。号内未知数的值。第第30页页cj322000CBXBbx1x2x3x4x5x60 x4(b)1111000 x515(a
13、)120100 x6202(c)1001c j z j322000CBXBbx1x2x3x4x5x60 x45/400(d)(l)-1/4-1/43x125/410(e)03/4(i)2x25/201(f)0(h)1/2c j z j0(k)(g)0-5/4(j)第第31页页解:因为基变量所对应的系数列向量为单位向量解:因为基变量所对应的系数列向量为单位向量,故,故 l=1。显然,。显然,k=0。下面,关键是求下面,关键是求 B-1。2/104/304/14/111hiB第第32页页 2/54/254/520152/104/304/14/11bhi2/1 4/1 10 hib,2/54/254
14、/51015204/454/35hib第第33页页 2/12/104/14/304/14/111B第第34页页 010212/12/104/14/304/14/11a2 a 0102/112/14/34/12/1aaa第第35页页 100112/12/104/14/304/14/11c3 c 1002/12/14/14/34/14/3ccc第第36页页 fed1212/12/104/14/304/14/112/1 4/5 4/1 fed,fed2/14/54/1第第37页页最后可计算出:最后可计算出:g=-3/4。第第38页页例例 已知某线性规划问题,用单纯形法计算时得到中间已知某线性规划问题
15、,用单纯形法计算时得到中间某两步的计算表如下表,试将表中空白处数字填上。某两步的计算表如下表,试将表中空白处数字填上。第第39页页cj354000CBXBbx1x2x3x4x5x65x28/32/3101/3000 x514/3-4/305-2/3100 x620/35/304-2/301c j z j-1/304-5/300CBXBbx1x2x3x4x5x65x215/418/41-10/414x3-6/415/414/413x1-2/41-12/4114/41c j z j第第40页页解:解:41/1541/12041/441/5041/1041/811BCBXBbx1x2x3x4x5x6
16、5x215/418/41-10/414x3-6/415/414/413x1-2/41-12/4114/41c j z j10013/5403/4503/201 8/4150/4144/41010001第第41页页(1)确定初始基变量,并计算初始可行基的逆矩阵)确定初始基变量,并计算初始可行基的逆矩阵B-1;(2)检验各非基变量的检验数:)检验各非基变量的检验数:如果如果j 0,j=m+1,.,n,则当前解就是最优解,则当前解就是最优解,最有解为,最有解为 B-1b,停止计算:,停止计算:ABCCB1 第第42页页否则否则,转下步;,转下步;如果如果j 0,j=m+1,.,n 中,若某个非基变量
17、中,若某个非基变量 xj 的系数列向量的系数列向量 Pk 0,则此问题无界,停止计算,否,则此问题无界,停止计算,否则转下步;则转下步;(4)根据)根据max(j 0)=k,来确定,来确定 x k 为换入变量为换入变量(最大检验数有多个时,可任选一个);(最大检验数有多个时,可任选一个);第第44页页(5)根据下式来确定)根据下式来确定 x l 为换出变量:为换出变量:最小比值有多个时(下一步迭代会出现非基变量最小比值有多个时(下一步迭代会出现非基变量=0的情况,即退化现象),可任选一个(通常选择下的情况,即退化现象),可任选一个(通常选择下标较小的那个基变量标较小的那个基变量 x l 换出)
18、。换出)。lklikikiiPBbBPBPBbB)()()0)()()(min11111 (6)从而确定新的基变量,计算新的可行基的逆矩)从而确定新的基变量,计算新的可行基的逆矩阵阵 B-1。(7)重复(重复(2)(6),直到找到问题的最优解。),直到找到问题的最优解。第第45页页max z=2x1+3x2 0,12 4 16 482 .212121xxxxxxts例:例:第第46页页解:解:max z=2x1+3x2 0,12 4 16 48 2 .543215241321xxxxxxxxxxxxts(1)XB=(x3,x4,x5)100010001B 1000100011B第第47页页 3
19、412,28min40212168min40210001000112168100010001min 3240042110001000100032 2x5x第第48页页 2,416,12min0413162min0414/1000102/101121684/1000102/101min 4/321004014/1000102/10130002 1x3x(2)XB=(x3,x4,x2)400010201B 4/1000102/1011B第第49页页 44/13,28,min4/122/1382min1004/1002142/101121684/1002142/101min 4/121000014/
20、1002142/10130200 5x4x(3)XB=(x1,x4,x2)400014201B 4/1002142/1011B第第50页页 8/12/300100108/12/112/1204/1030200 (4)XB=(x1,x5,x2)410004201B 08/12/112/1204/101B 2441216808/12/112/1204/101251*bBxxxX 141216808/12/112/1204/103021*bBCzB第第51页页对偶是指对同一问题,从不同角度观察,有两种对对偶是指对同一问题,从不同角度观察,有两种对立的表述。立的表述。周长一定,面积最大的矩形是正方形。
21、周长一定,面积最大的矩形是正方形。面积一定,周长最小的矩形是正方形。面积一定,周长最小的矩形是正方形。第第52页页cyxts )(2 .xymax xcx21max221 maxxcx 0221)21(2 xcdxxcxdcx41 cy41 正方形正方形非线性规划非线性规划第第53页页sxyts .yx maxxsx max01)(2 xsdxxsxdsx sy 正方形正方形非线性规划非线性规划xsx max第第54页页问题问题1 某工厂有某工厂有 I、II 两种产品,已知生产单位产品两种产品,已知生产单位产品所需的设备台时数、原材料所需的设备台时数、原材料 A 的消耗量、原材料的消耗量、原材
22、料B的的消耗量,具体数据如表所示。消耗量,具体数据如表所示。设设 备备原材料原材料A原材料原材料BI140II2048台时台时16kg12kg每生产一件产品每生产一件产品 I 可获利可获利 2 元,每生产一件产品元,每生产一件产品 II 可可获得获得3元,问如果自己生产,产品元,问如果自己生产,产品 I 和和 II 各应生产多各应生产多少获利最大?少获利最大?第第55页页又知道,每生产一件产品又知道,每生产一件产品 I 可获利可获利 2 元,每生产一件元,每生产一件产品产品 II 可获得可获得 3 元,问应如何安排生产计划才能够使元,问应如何安排生产计划才能够使得工厂获利最多?得工厂获利最多?
23、0,12416482.32max21212121xxxxxxtsxxz第第56页页问题问题2 某工厂有某工厂有 I、II 两种产品,已知生产单位产品两种产品,已知生产单位产品所需的设备台时数、原材料所需的设备台时数、原材料 A 的消耗量、原材料的消耗量、原材料B的的消耗量,具体数据如表所示。消耗量,具体数据如表所示。设设 备备原材料原材料A原材料原材料BI140II2048台时台时16kg12kg每生产一件产品每生产一件产品 I 可获利可获利 2 元,每生产一件产品元,每生产一件产品 II 可可获得获得 3 元,问如不自己生产,而是将设备、原材料元,问如不自己生产,而是将设备、原材料 A、原材
24、料、原材料 B 等三种资源卖出,定价多少才能将资源卖等三种资源卖出,定价多少才能将资源卖出?出?第第57页页设设备台时数售价:设设备台时数售价:y1 资源资源 A 售价:售价:y2 资源资源 B 售价:售价:y3生产产品生产产品 I 所需的各种资源的售价总和,不应低所需的各种资源的售价总和,不应低于自己生产产品于自己生产产品 I 所带来的利润:所带来的利润:生产产品生产产品 II 所需的各种资源的售价总和,不应所需的各种资源的售价总和,不应低于自己生产产品低于自己生产产品 II 所带来的利润:所带来的利润:2421 yy34221 yy第第58页页 0,34224.12168 min32131
25、21321yyyyyyytsyyy 从工厂角度来看:从工厂角度来看:越大越好;越大越好;从接受者角度来看:从接受者角度来看:越小越好。越小越好。所以,工厂只能在满足大于等于自己生产产品所带来所以,工厂只能在满足大于等于自己生产产品所带来的利润的条件下,提出一个尽可能低的出售价格,才的利润的条件下,提出一个尽可能低的出售价格,才能实现其意愿,因此问题的目标函数是极小化。能实现其意愿,因此问题的目标函数是极小化。32112168 yyy 卖出的总收入为:卖出的总收入为:第第59页页 0,34224.12168 min3213121321yyyyyyytsyyy 0,12416482.32max21
展开阅读全文