1、III每天可用能力每天可用能力设备设备A（h）设备设备B（h）调试工序（调试工序（h）06152115245利润（元）利润（元）21例1 0,52426155s.t.2max212121221xxxxxxxxxz设设y1,y2和和y3分别表示出让资源分别表示出让资源A，B和调试和调试工序的单价，则美佳公司同意出让的条件工序的单价，则美佳公司同意出让的条件将是将是同意出让生产产品同意出让生产产品I I的资源的资源同意出让生产产品同意出让生产产品IIII的资源的资源购买者希望用最少的代价获得这些资源购买者希望用最少的代价获得这些资源，因此因此2632 yy125321yyy32152415miny
2、yyz这样得到一个新的线性规划问题这样得到一个新的线性规划问题0,1252652415min32132132321yyyyyyyyyyyw称这一问题是原来的称这一问题是原来的LP问题的问题的对偶线性规对偶线性规划问题划问题或或对偶问题对偶问题，原来的，原来的LP问题也称为问题也称为原问题原问题。项目项目原原问题问题对偶问题对偶问题AbC目标函数目标函数约束条件约束条件决策变量决策变量约束条件系数矩阵约束条件系数矩阵约束条件右端项向量约束条件右端项向量目标函数系数向量目标函数系数向量max z=CX AXb X0约束条件系数矩阵转置约束条件系数矩阵转置目标函数的系数向量目标函数的系数向量约束条件
3、的右端项向量约束条件的右端项向量min w=YbAY CY 0原问题原问题max z对偶问题对偶问题min wn个决策变量个决策变量m个约束条件个约束条件n个约束条件个约束条件m个决策变量个决策变量约束条件约束条件“”型型决策变量决策变量0决策变量决策变量0约束条件约束条件“”型型对称形式的对应关系对偶问题的对偶是原问题，即对偶关系是对偶问题的对偶是原问题，即对偶关系是相互对称的关系相互对称的关系非对称形式下的对偶关系原问题原问题（对偶问题）（对偶问题）max z对偶问题对偶问题（原问题）（原问题）min wn个决策变量个决策变量m个约束条件个约束条件n个约束条件个约束条件m个决策变量个决策变
4、量约束条件约束条件“”型型约束条件约束条件“”型型约束条件约束条件“=”型型决策变量决策变量0决策变量决策变量0决策变量无约束决策变量无约束决策变量决策变量0决策变量决策变量0决策变量无约束决策变量无约束约束条件约束条件“”型型约束条件约束条件“”型型约束条件约束条件“=”型型0maxXbAXCXzNBNBCCCNBAXXX,00maxXbIXAXSXCXzsbBXBNXBIXbIXNXBXSNBSNB111原来添加松弛变量XS将XB的系数矩阵化为单位矩阵0, 00maxNBSNBXXCXCzXXbIXNXBXsNNBBCB CN 0XB XN XS0 XS bB N ICB CN 0CB C
5、N 0XB XN XSCB XB B-1bI B-1N B-10 CN CBB-1N CBB-1 初始单纯形表迭代后的单纯形表 在初始单纯形表中单位矩阵经过迭代后变为基矩阵B的逆 在初始单纯形表给出的解中基变量Xs=b,而在迭代后的表给出的解中基变量 XB=B-1b 系数矩阵的变化: A, IB-1A, I 在初始单纯形表中变量xj的系数为Pj经过迭代后变为Pj,并且Pj=B-1 Pj 若迭代后的单纯形表为最终表则该表也同时给出对偶问题的最优解项目原问题变量原问题松弛变量 x1 x2x3 x4 x5x3 15/2 0 01 5/4 15/2x1 7/2 1 00 1/4 -1/2x2 3/2
6、0 10 -1/4 3/2-j 0 00 1/4 1/2对偶问题剩余变量对偶问题变量 y4 y5y1 y2 y3项目对偶问题变量对偶问题剩余变量y1 y2 y3y4 y5y2 1/4-5/4 1 0-1/4 1/4 y3 1/215/2 0 1 1/2 -3/2j15/2 0 07/2 3/2原问题松弛变量原问题变量 x3 x4 x5x1 x2原问题最终单纯形表对偶问题最终单纯形表例1最大化问题检验数的相反数给出了对偶问题的解原本在对偶关系中，原问题的变量对应着对偶问题原本在对偶关系中，原问题的变量对应着对偶问题的约束条件，原问题的约束条件对应着对偶变量。的约束条件，原问题的约束条件对应着对偶
7、变量。但但在分别添加了松弛变量和剩余变量后，也可以建在分别添加了松弛变量和剩余变量后，也可以建立原问题变量与对偶问题变量之间的对应关系立原问题变量与对偶问题变量之间的对应关系原问题原问题对偶问题对偶问题第第i个约束条件中个约束条件中添加的松弛变量添加的松弛变量第第i个对偶变量个对偶变量第第j个变量个变量第第j个约束条件中个约束条件中添加的松弛变量添加的松弛变量注 上表中我们将松弛变量与剩余变量统称为松弛变量 弱对偶性原问题可行解的目标函数不超过对偶问题可行解的目标函数（1）原问题任一可行解的目标函数值是其对偶）原问题任一可行解的目标函数值是其对偶问题目标函数值的下界；反之对偶问题任一可行问题目
8、标函数值的下界；反之对偶问题任一可行解的目标函数值是原问题目标函数值的上界。解的目标函数值是原问题目标函数值的上界。（2）如原问题有可行解且目标函数无界（即原）如原问题有可行解且目标函数无界（即原问题为无界解），则对偶问题无可行解；反之对问题为无界解），则对偶问题无可行解；反之对偶问题有可行解且目标函数无界，则原问题无可偶问题有可行解且目标函数无界，则原问题无可行解。注意该推论的逆命题不成立。行解。注意该推论的逆命题不成立。（3）若原问题有可行解而对偶问题无可行解，）若原问题有可行解而对偶问题无可行解，则原问题目标函数无界；反之对偶问题有可行解则原问题目标函数无界；反之对偶问题有可行解而原问题
9、无可行解，则原问题目标函数无界。而原问题无可行解，则原问题目标函数无界。最优性若原问题一个可行解目标函数等于对偶问题的某个可行解的目标函数,则这两个可行解分别是原问题和对偶问题的最优解强对偶性若原问题和对偶问题都有可行解,则它们都有最优解,且最优解的目标函数值相等互补松弛性在线性规划问题的最优解中,如果对应某一约束条件的对偶变量值非零,则其对应的约束条件取等式;反之若一个约束条件为严格的不等式,则其对应的对偶变量为零在线性规划问题的最优解中,如果对应某一约束条件的对偶变量值非零,则该约束条件中松弛变量等于零;反之若一个约束条件中松弛变量非零,则其对应的对偶变量为零。例（p76.7）4 , 3
10、, 2 , 1, 096628342max321432214214321jxxxxxxxxxxxxxxxxzj4 , 3 , 2 , 1, 01143229668min314343214214321iyyyyyyyyyyyyyyyywi而由x1=20, 得22421yyy而由x2=20, 得434321yyyy而由x3=40, 得143 yy于是得到方程组0143224434321421yyyyyyyyyy得对偶问题最优解为0, 1,53,544321yyyy注：原问题与对偶问题最优目标函数值都是 z*=4+8+4=16式中bi是线性规划原问题约束条件的右端项，它代表第i种资源的拥有量；对偶变
11、量yi的意义代表在资源最优利用的条件下对第i种资源的估价。这种估价不是资源的市场价格，而是根据资源在生产中作出的贡献而作的估价，为区别起见，称为影子价格影子价格。*1*1*wybxczmiiinjjj设 和 分别是原问题和对偶问题的最优解，则由对偶性质，有), 1(*njxj), 1(*niyi 资源的影子价格随企业的生产任务、产品结构资源的影子价格随企业的生产任务、产品结构的改变而改变的改变而改变 影子价格是资源的影子价格是资源的边际价格边际价格 资源的影子价格也可视为一种资源的影子价格也可视为一种机会成本机会成本 在生产过程中若某种资源未得到充分利用则其在生产过程中若某种资源未得到充分利用
12、则其影子价格为零；只有在资源得到充分利用时，影子价格为零；只有在资源得到充分利用时，其影子价格才可能非零其影子价格才可能非零 利用影子价格可以说明：单纯形法中的检验数利用影子价格可以说明：单纯形法中的检验数可以看成生产某种产品的产值与隐含成本的差可以看成生产某种产品的产值与隐含成本的差 可以利用影子价格确定企业内部的核算价格，可以利用影子价格确定企业内部的核算价格，以便控制有限资源的使用和考核下属企业经营以便控制有限资源的使用和考核下属企业经营的好坏。的好坏。例1最优目标函数值的变化：8.5变到8.75，增加1/4资源的变化：设备B的可用时间从增加一小时参考文献：参考文献：李慧：资源影子价格分
13、析与经营管理决策，李慧：资源影子价格分析与经营管理决策，系统工程理论与实践系统工程理论与实践，2003年年4月号，月号，22-26单纯形法求解的基本思路单纯形法求解的基本思路基可行解基可行解检验数非正检验数非正保持解的可行性保持解的可行性对偶单纯形法的基本思路对偶单纯形法的基本思路对偶问题基可行解对偶问题基可行解（检验数非正）（检验数非正）原问题基可行原问题基可行解解保持对偶问题解的保持对偶问题解的可行性（检验数非可行性（检验数非正正（对偶问题可行解）（对偶问题可行解） 保持对偶问题有基可行解,而原问题只是基本解,通过迭代,使后者的负分量个数减少，一旦成为基可行解,则原问题与对偶问题同时实现最
14、优解. 适应于求解这样的适应于求解这样的LP问题：问题：标准化后不标准化后不含初始基变量，但将某些约束条件两端含初始基变量，但将某些约束条件两端乘以乘以“-1”后，即可找出初始基变量。后，即可找出初始基变量。要求：要求：初始单纯形表中的检验数满足最初始单纯形表中的检验数满足最优性条件优性条件对满足上述条件的对满足上述条件的LP问题，对偶单纯形法的问题，对偶单纯形法的步骤是：步骤是：旋转运算。然后回到第旋转运算。然后回到第2步。步。作出初始单纯形表（注意要求）作出初始单纯形表（注意要求）检查检查b列的数据是否非负，若是，表中已经列的数据是否非负，若是，表中已经给出最优解；否则转下一步给出最优解；
15、否则转下一步确定换出变量确定换出变量：取：取b列最小的数对应的变量列最小的数对应的变量为换出变量为换出变量确定换入变量确定换入变量：用检验数去除以换出变量行：用检验数去除以换出变量行的那些对应的负系数，在除得的商中选取其的那些对应的负系数，在除得的商中选取其中最小者对应的变量为换入变量中最小者对应的变量为换入变量例 用对偶单纯形法求解如下的LP问题0,1252652415min32132132321xxxxxxxxxxxz化成标准形式0,1252652415max3215321432321xxxxxxxxxxxxxw0,1252652415max3215321432321xxxxxxxxxxx
16、xxw将各约束条件两端同乘“-1”得用对偶单纯形法求解得注：通常很少直接使用对偶单纯形法求解线性规划问题。灵敏度分析 将讨论将讨论LP问题中的参数问题中的参数 中有一个或几个发生改变时问题的最优中有一个或几个发生改变时问题的最优解会有什么变化，或者这些参数在一个解会有什么变化，或者这些参数在一个多大的范围内变化时，问题的最优解不多大的范围内变化时，问题的最优解不变变ijijabc,研究的思路 将个别参数的变化直接在计算得到的最将个别参数的变化直接在计算得到的最终单纯形表中反映出来，这样就不需要终单纯形表中反映出来，这样就不需要从头计算，而从头计算，而直接检查在参数改变后最直接检查在参数改变后最
17、终表有什么改变终表有什么改变，若仍满足最终表的条，若仍满足最终表的条件，则表中仍给出最优解，否则从这个件，则表中仍给出最优解，否则从这个表开始进行迭代求改变以后的最优解。表开始进行迭代求改变以后的最优解。灵敏度分析的步骤灵敏度分析的步骤 将参数的改变计算反映到最终表上来。具体计将参数的改变计算反映到最终表上来。具体计算公式可以使用算公式可以使用 检查原问题是否仍为可行解检查原问题是否仍为可行解 检查对偶问题是否仍为可行解检查对偶问题是否仍为可行解 对检查情况按下表进行处理对检查情况按下表进行处理jBjjjjPBCcPBPbBb111原原问题问题对偶问题对偶问题结论或继续计算步骤结论或继续计算步
18、骤可行解可行解可行解可行解问题的最优解或最优基不变问题的最优解或最优基不变可行解可行解非非可行解可行解用用单纯形法继续迭代求最优单纯形法继续迭代求最优解解非非可行解可行解可行解可行解用用对偶单纯形法继续迭代求对偶单纯形法继续迭代求最优解最优解非非可行解可行解非非可行解可行解引进人工变量，编制新的单引进人工变量，编制新的单纯形表重新计算纯形表重新计算例：在第一章美佳公司的例例：在第一章美佳公司的例1中中（1）若产品）若产品I的利润降至的利润降至1.5元元/件，而产品件，而产品II的利润增至的利润增至2元元/件，美佳公司的最件，美佳公司的最优生产计划有何改变；优生产计划有何改变；（2）若产品）若产
19、品I的利润不变，则产品的利润不变，则产品II的利的利润在什么范围变化时，该公司的最优生润在什么范围变化时，该公司的最优生产计划不发生变化产计划不发生变化21000CbXbbx1x2x3x4x50 x37 1/2 0012 1/2-7 1/22x13 1/2 1001/4- 1/21x21 1/2 010- 1/41 1/2000- 1/4- 1/2原最终单纯形表原最终单纯形表（1）改变后）改变后0 x46004/51-61.5 x1210- 1/5012x23011/50000-1/100-1 1/2新的最优解为：0, 6, 0, 3, 254321xxxxx最优目标函数值为：9*z（2）改变
20、后）改变后为使表中的解仍为最优解必须为使表中的解仍为最优解必须0231, 02141aa因此产品因此产品II的利润变化范围为的利润变化范围为232 a例：在第一章美佳公司的例例：在第一章美佳公司的例1中中（1）若设备）若设备A与调试工序的每天能力不变，与调试工序的每天能力不变，而设备而设备B每天的能力增加到每天的能力增加到32小时，分析小时，分析公司最优计划的变化；公司最优计划的变化；（2）若设备）若设备A和和B每天可用能力不变，则每天可用能力不变，则调试工序能力在什么范围变化时，问题调试工序能力在什么范围变化时，问题的最优基不变的最优基不变（1）b由由(15，24，5)T 变为变为(15，3
21、2，5)T 后，相应地最后，相应地最终表中终表中b列的数据变为列的数据变为2/ 12/112/35532152/ 34/ 102/ 14/ 102/154/ 511bBb21000CbXbbx1x2x3x4x50 x317 1/20012 1/2-7 1/22 x15 1/21001/4- 1/21 x2- 1/2010- 1/41 1/2000- 1/4- 1/20 x312 1/2010107 1/22 x15110010 x420-401-60-100-2代入原最终表（2）设现在每天调试工序的时间为x，则最终表中b列的数变为xxxxbBb2362162154524152/ 34/ 102
22、/ 14/ 102/154/511故要使最优基不变必须6,4412602360216021545xxxxxxx利用Excle求解LP问题，以P45.7(2)为例变量，已经赋了初值目标函数值约束条件右端值其他专业软件：Lindo与Lingo，WinQSB例如Lingo，启动Lingo后，按图中的方式输入模型，然后点击求解的图标 。就可得到所需的最优解。24对偶问题为0, 06353322232max212121212121yyyyyyyyyyyyw由图解法可得对偶问题最优解为 .519*;51,5821wyy将该最优解代入对偶问题约束条件可知，第四个约束条件为严格不等式，因此在原问题最优解中 .
23、 04x而由于 0,21yy因此将原问题最优解代入原问题约束条件，它们成为等式。再由于原问题最优目标函数值等于对偶问题最优目标函数值。于是原问题最优解满足方程组 5/1953232232321321321xxxxxxxxx解方程组得原问题最优解： 0,51, 0,574321xxxx2.5对偶线性规划为 0, 0121222min321321321321321yyyyyyyyyyyyyyyw无约束，（2）直接观察可知对偶问题有解对应于该解的目标函数值 由弱对偶性，原问题的任何可行解的目标函数值都满足 0, 1, 0321yyy1w. 1 wz精品课件精品课件！精品课件精品课件！0, 0, 0321xxx. 0,01122min2121212121yyyyyyyyyyw0,21yy2.6 原问题显然有可行解。但对偶问题第一个约束条件与非负约束条件冲突，实际上，当时，第一个约束条件左端非正，因此不可能不小于正数1。这说明这些约束条件不能同时成立。故对偶问题无可行解。由弱对偶性，原问题目标函数无界。
