《第2章对偶理论与灵敏度分析(5h)16剖析》由会员分享，可在线阅读，更多相关《第2章对偶理论与灵敏度分析(5h)16剖析（54页珍藏版）》请在金锄头文库上搜索。
1、第二章 对偶理论与灵敏度分析,1线性规划的对偶问题 2对偶问题的基本性质 3影子价格 4对偶单纯形法 5现性规划问题的矩阵形式 6灵敏度分析,1线性规划的对偶问题,1.1 对偶问题的提出 1.2 对称形式下对偶问题的一般形式 1.3 非对称形式的原对偶问题关系 1.4 对偶问题的定义 1.5 对偶关系对应表,1.1 对偶问题的提出,1线性规划的对偶问题,如何组织生产，才能获得最大利润？假定有另一公司想把美佳公司的资源收买过来，它至少应付出多大代价，才能使美佳公司愿意放弃生产活动，出让自己的资源？,max Z=c1x1+ c2x2+cnxn,1.2 对称形式下对偶问题的一般形式,1线性规划的对偶
2、问题,min w=b1y1+ b2y2+bmym,1.3 对偶对应表,1线性规划的对偶问题,1线性规划的对偶问题,2对偶问题的基本性质,一、对称性 二、弱对偶性 三、无界性 四、可行解是最优解的性质 五、对偶定理 六、松紧定理,一、对称性 ： 对偶问题的对偶是原问题,2. 对偶问题的基本性质,二、弱对偶性：若x是原问题的可行解，y是对偶问题的可行解。则有 cxyb,推论（1）原问题任一可行解的目标函数值是其对偶问题目标函数值的下界，反之对偶问题任一可行解的目标函数值是其原问题目标函数值的上界。（2）若原问题（对偶问题）为无界解，则其对偶问题（原问题）无可行解。（3）若原问题有可行解而其对偶问题
3、无可行解，则原问题目标函数值无界，反之对偶问题有可行解而其原问题无可行解，则对偶问题的目标函数值无界。,三、最优性：设x是原问题的可行解，y是对偶问题的可行解.当cx=yb时 x, y是最优解。,四、 强对偶性（对偶定理）：若原问题及其对偶问题均具有可行解，则两者均具有最优解且它们最优解的目标函数值相等。,五、互补松弛性（松紧定理）在线性规划问题的最优解中，如果对应某一约束条件的对偶变量值为非零，则该约束条件取严格等式；反之如果约束条件取严格不等式，则其对应的对偶变量一定为零。,2. 对偶问题的基本性质,1.证明数学模型无最优解,2、考虑问题回答： 1 ）说明原问题和对偶问题都有最优解. 2
4、）求原问题和对偶问题的最优目标函数值的一个上界和下界.,2. 对偶问题的基本性质,1）解：原问题有可行解（0,0,8）,对偶问题有可行解（5,0） 强对偶性定理则有最优解,2. 对偶问题的基本性质,3、线性规划问题有最优解（0,0,4,4）最优值28。用互补松弛定理计算对偶问题的最优解。,解：原问题的对偶问题为,3影子价格,式中bi是线性规划原问题约束条件的右端项，它代表第i种资源的拥有量；对偶变量yi*的意义代表在资源最优利用条件下对单位第i仲资源的估价。这种估价不是资源的市场价格，而是根据资源在生产中作出的贡献而作的估价，为区别起见，称为影子价格(shadow price)。,这说明影子价
5、格的值相当于在资源得到最优利用的生产条件下，b每增加一个单位时目标函数Z的增量。,关于影子价格的几点说明： 1）影子价格是未知数,有赖于资源的利用情况。 2）影子价格是一种边际价格。,例1,A B C,15/2 0 0,资源,剩余,0 1/4 1/2,影子价格,3影子价格,所以第2种资源的影子价格为0.25，第3种资源的影子价格为0.5。,3影子价格,3）资源的影子价格实际上又是一种机会成本。 在纯市场经济条件下，可以根据影子价格与市场价格的关系确定资源的买进或卖出。 4）生产过程中如果某种资源末得到充分利用时，该种资源的影子价格为零；又当资源的影子价格不为零时，表明该种资源在生产中已耗费完毕
6、。,3影子价格,5）单纯形表中各个检验数的经济意义。当产品产值大于隐含成本时，表明生产该项产品有利，可在计划中安排，否则这些资源来生产别的产品更为有利，就不在计划中安排。这就是单纯形表中各个检验数的经济意义。,6）对线性规划问题的求解是确定资源的最优分配方案，而对于对偶问题的求解则是确定对资源的恰当估价。,对偶单纯形法的基本思路：先找出一个对偶问题的可行基，并保持对偶问题为可行解条件下，如不存在XB o，通过变换到一个相邻的目标函数值较小的基本解(因对偶问题是求目标函数极小化)，并循环进行，一直到原问题也为可行解(即XB o)，这时对偶问题与原问题均为可行解。 在迭代过程中保持原问题的检验数为
7、非正，逐步替换负基变量，从而得到最优解。,4对偶单纯形法,对偶单纯形法的计算步骤：,1. 列出初始单纯形表，且检验数非正。,4对偶单纯形法,4.以ars为主元素，进行迭代变换。,5 . 返 2，直到B 0为止。,原始单纯形法 按列选主元 对偶单纯形法 按行选主元,例5：用对偶单纯形法求解下列问题,4对偶单纯形法,4对偶单纯形法,4对偶单纯形法,4对偶单纯形法,已获得最优解： （x1, x2, x3, x4, x5, x6）=（5, 7, 6, 0, 0, 0） min z=35 对偶问题的最优解为： （w1, w2, w3, w4, w5, w6）=（-1, 5, 7, 0, 0, 0） ma
8、x y=35,4对偶单纯形法,从以上计算可以看出，用对偶单纯形法求解线性规划问题时，当约束条件为“”时，不必引进人工变量，使计算简化。但在初始单纯形表中其对偶问题应是基可行解这点，对多数线性规划问题很难实现。因此对偶单纯形法一般不单独使用，而主要应用于灵敏度分析及整数规划等有关章节中。,4对偶单纯形法,对偶单纯形法的优点： 1 简化计算 2 减少计算量 3 灵敏度分析 缺点： 第一个正则解很难找到,4对偶单纯形法,单纯形法计算的矩阵描述,对称形式线性规划问题的矩阵表达式加上松弛变量后为：,5单纯形法的矩阵描述,上式中XS为松弛变量,I 为mm单位矩阵,通过迭代后: X=XB + XN,相应地:
9、 C=CB + CN,系数矩阵 : A= B + N,参1-105,5单纯形法的矩阵描述,参1-105,迭代后新的单纯形表为：,初始单纯形表：,C- CBB-1A,5单纯形法的矩阵描述,可以看出：,(1)对应初始单纯形表中的单位矩阵I，迭代后的单纯形表中为B-1 ； (2)初始单纯形表中基变量Xsb，迭代后的表中 XB=B1b ； (3)初始单纯形表中约束系数矩阵为A，I B，N，I ，迭代后的表中约束系数矩阵为B-1A，B-1IB-1B，B-1N，B-1I I，B-1N，B-1 ； (4)若初始矩阵中变量Xi的系数向量为Pj，迭代后为Pj，则有Pj=B-1Pj 。,5单纯形法的矩阵描述,当B
10、为最优基时,因xB的检验数可写为 CB- CBI=0,CBB-1称为单纯形乘子，若令Y= CBB-1 ，则,所以,CN- CBB-1N 0 -CBB-1 0,C- CBB-1A 0 -CBB-1 0,5单纯形法的矩阵描述,6灵敏度分析,灵敏度分析是指对系统或事物因周围条件变化显示出来的敏感程度的分析。A，b，C 当这些参数中的一个或几个发生变化时，问题的最优解会有什么变化，或者这些参数在一个多大范围内变化时，问题的最优解不变。这就是灵敏度分析所要研究解决的问题。 单纯形法的迭代计算是从一组基向量变换为另一组基向量，表中每步迭代得到的数字只随基向量的不同选择而改变，因此有可能把个别参数的变化直接
11、在计算得到最优解的最终单纯形表上反映出来。,灵敏度分析的步骤 1将参数的改变计算反映到最终单纯形表上来 2检查原问题是否仍为可行解； 3检查对偶问题是否仍为可行解； 4按下表所列情况得以结论和决定继续计算的步骤。,6灵敏度分析,6灵敏度分析,(1)若家电I的利润降至1. 5元件，而家电的利润增至2元件时，美佳公司最优生产计划有何变化?,变量系数Cj的变化仅仅影响到检验数(Cj-Zj)约变化。,一、分析Cj的变化,6灵敏度分析,解 (1)将家电I，II的利润变化直接反映到最终单纯形表中。,以x4为换入变量,x3为换出变量，列新单纯形表。,所有检验数为非正，问题得到最优解：X（2，3，0，6，0）
12、,6灵敏度分析,(2)若家电I的利润不变，则家电II的利润在什么范围内变化时，则该公司的最优生产计划将不发生变化？ 设家电II的利润为(1十)元，反映到最终单纯形表中为：,为使表中的解仍为最优解，应有,即家电II的利润C2的变化范围应满足：,二、分析bi的变化,右端项bj的变化在实际问题中反映为可用资源数量的变化。,6灵敏度分析,解：因为：,(1)若设备A和调试工序的每天能力不变，而设备B每天的能力增加到32小时，分析公司最优计划的变化?,以X2为换出变量，X4为换入变量，列新单纯形表得：,问题的最优解：X（5，0，15，2，0）,6灵敏度分析,单纯形表中 b列数字为：,(2)若设备A和B每天
13、可用能力不变，则调试工序能力在什么范围内变化时，问题的最优基不变?,设调试工序每天可用能力为（5+）小时，因有,当b 0时问题的最优基不变，解得：-1 1 因此调试工序的能力应在4小时-6小时之间。,6灵敏度分析,三、增加一个变量xj的分析,增加一个变量在实际问题中反映为增加一种新的产品。 其分析步骤为：,3.若 ，原最优解不变，只需将计算得到的 和 直接写入最终单纯 形表中；若 ，则按单纯形法继续迭代计算找出最优。,6灵敏度分析,在美佳公司例子中，设该公司又计划推出新型号的家电III，生产一件所需设备A、B及调试工序的时间分别为3小时、4小时、2小时，该产品的预期盈利为3元件。试分析该种产品
14、是否值得投产，如投产，对该公司的最优生产汁划有何变化。,解：设该公司生产家电III 件，有,6灵敏度分析,以x2 为换出变量， x6 为换入变量，列新单纯形表得：,6灵敏度分析,问题的最优解：X（7/2，0，51/4，0，0，3/4）T。 所以，该公司新的生产计划应为每天生产家电I 7/2件，生产家电III 3/4件。,6灵敏度分析,四、分析参数aij的变化,aij的变化使线性规划的约束系数矩阵A发生变化。 若变量xj在最终单纯形表中为非基变量，其约束条件中系数aij的变化分析步骤可参照本节之三； 若变量xi在最终单纯形表中为基变量，则aij的变化将使相应的B和B-1 发生变化因此有可能出现原
15、问题和对偶问题均为非可行解的情况。 出现这种情况时，需引进入工变量将原问题的解转化为可行解，再用单纯形法求解。,6灵敏度分析,若家电II每件需设备A，B和调试工时变为8小时、4小时、1小时，该产品的利润变为3元/件，试重新确定该公司最优生产计划。,解：将生产工时变化后的新家电II看作是一种新产品，生产量为x2。,6灵敏度分析,因x2已变换为x2 ，故用单纯形算法将x2替换出基变量中的x2 ，并在单纯形表中不再保留x2。,6灵敏度分析,增加人工变量。,X3+ 4X4 - 24X5= -9,6灵敏度分析,美佳公司的最优生产计划为每天生产家电I 11/4件，新家电II 15/8件。,6灵敏度分析,五
16、、增加一个约束条件的分析,增加一个约束条件在实际问题中相当增添一道工序。从图形（二维）上分析相当于增加了一条直线。分析的方法是先将原问题最优解的变量值代入新增的约束条件。如满足，说明新增的约束未起到限制作用，原最优解不变。否则，将新增的约束直接反映到最终单纯形表中再进步分析。,6灵敏度分析,设家电I，II经调试后，还需经过道环境试验工序。家电I每件须环境试验3小时，家电II每件2小时，又环境试验工序每天生产能力为12小时试分析增加该工序后的美佳公司最优生产计划。 解：先将原问题的最优解x1 =7/2, x2 =3/2代入环境试验工序的约束条件： 3x1 +2 x2 12。因为37/2 + 23/2 12，故原问题最优解不是本例的最优解。,6灵敏度分析,
