运筹学—线性规划.ppt





《运筹学—线性规划.ppt》由会员分享,可在线阅读,更多相关《运筹学—线性规划.ppt(181页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、运筹学运筹学线性规划线性规划(5 5)根据根据 ,确定基变量,确定基变量 为换出变量。为换出变量。(6 6)用替代用替代 得新基得新基1 1,由变换公式,由变换公式计算新基的逆阵计算新基的逆阵1 1-1-1,求出新的基本可行解,求出新的基本可行解 其中为变换矩阵,构造方法是:其中为变换矩阵,构造方法是:从一个单位矩阵出发,把换出变量从一个单位矩阵出发,把换出变量 在基在基B B中的对应列的单位中的对应列的单位 向量,替换成换入变量向量,替换成换入变量 对应的系数列向量对应的系数列向量 ,并做如下变形,并做如下变形,主元素主元素 (应在主对角线上)取倒数,其它元素除以主元素(应在主对角线上)取倒
2、数,其它元素除以主元素 并取相反数。并取相反数。重复(重复(2 2)()(6 6)直至求得最优解。)直至求得最优解。换入换入无界解无界解换出换出最优解最优解初始基初始基新基新基 例例9 9、试用改进单纯形法求解、试用改进单纯形法求解解:解:()观察法确定()观察法确定,为基变量为基变量为非基变量为非基变量 (2 2)计算单纯行乘子)计算单纯行乘子 目标函数当前值目标函数当前值 (3 3)非基变量检验数)非基变量检验数(4(4)选择换入变量选择换入变量 故故 为换入变量。计算为换入变量。计算()确定换出变量确定换出变量确定确定 为换出变量,主元素为换出变量,主元素(6)(6)用用 代替代替 得新
3、可行基得新可行基 为基变量,为基变量,为非基变量,为非基变量,重复以上步骤,进入第二循环重复以上步骤,进入第二循环(1)(1)(2)(2)(3)(3)(4)(4)选择选择 换入变量换入变量(5)(5)选择选择 换出变量换出变量,主元素主元素(6)(6)用用 代替代替 得新可行基得新可行基 为基变量,为基变量,为非基变量,为非基变量,进入第三循环进入第三循环.(1)(1)(2)(2)(3)(3)非基变量对应的检验数全部非正,非基变量对应的检验数全部非正,故当前解故当前解 即为最优解,即为最优解,相应的目标函数最优值相应的目标函数最优值 。六、单纯形法总结六、单纯形法总结1、Min型单纯形表与型单
4、纯形表与Max型的区别仅在于:型的区别仅在于:令令k=minj0的的xk进基,当进基,当0时最优。时最优。2、解的几种情况及其在单纯形表上的体现(讨论、解的几种情况及其在单纯形表上的体现(讨论Max型)型)唯一唯一最优解最优解j0,非基非基0但但B-1Pk0无可行解无可行解用大用大M法法求解,最求解,最优基中含优基中含有有X人人退化解退化解最优解最优解中某基中某基变量为变量为0第三节 对偶问题与灵敏度分析一一、对偶问题及其模型、对偶问题及其模型例例7:回顾例:回顾例1这时有另一家厂商这时有另一家厂商提出要购买其煤、电、提出要购买其煤、电、油全部资源,并希望花油全部资源,并希望花费尽量少。试建立
5、购买费尽量少。试建立购买者的线性规划模型。者的线性规划模型。例例7称为例称为例1的对偶问题,记为的对偶问题,记为(D),例),例1称为例称为例7的原问题的原问题,记为(记为(P)。)。对偶模型的一般式对偶模型的一般式以例以例7为例,原问题为为例,原问题为(P)(D)这是最常见的对偶模型形式,称为对称式对偶模型。二者间这是最常见的对偶模型形式,称为对称式对偶模型。二者间具有十分对称的对应关系:具有十分对称的对应关系:原问题(原问题(P P)对偶问题对偶问题 (D D)目标目标max型型目标目标min型型有有n个变量(非负)个变量(非负)有有n个约束(大于等于)个约束(大于等于)有有m个约束个约束
6、(小于等于)(小于等于)有有m个变量(非负)个变量(非负)价格系数价格系数资源向量资源向量资源向量资源向量价格系数价格系数技术系数矩阵技术系数矩阵技术系数矩阵的转置技术系数矩阵的转置此外,还有一种情形原问题(原问题(P P)对偶问题对偶问题 (D D)第第j个变量为自由变量个变量为自由变量第第j个约束为等式约束个约束为等式约束第第i个约束为等式约束个约束为等式约束第第i个变量为自由变量个变量为自由变量例例8:写出下面线性规划的对偶规划模型:写出下面线性规划的对偶规划模型:例例8:写出下面线性规划的对偶规划模型:写出下面线性规划的对偶规划模型:minZ=-CXs.T-AX-bX0对称性定理:对称
7、性定理:对偶问题的对偶是原问题。对偶问题的对偶是原问题。minW=Ybs.TYACY0maxZ=Cs.TAXbX0对偶的定义对偶的定义对偶的定义对偶的定义maxW=-Ybs.TYACY0二、对偶问题的性质弱对偶原理:弱对偶原理:设设 和和 分别是问题(分别是问题(P P)和()和(D D)的)的可行解,则必有可行解,则必有推论推论:若若 和和 分别是问题(分别是问题(P P)和()和(D D)的可行)的可行解,则解,则 是(是(D D)的目标函数最小值的一)的目标函数最小值的一个下界;个下界;是(是(P P)的目标函数最大值的一)的目标函数最大值的一个上界。个上界。对偶问题的性质关于无界性有如
8、下结论:关于无界性有如下结论:问题无问题无界界无可无可行解行解无可无可行解行解问题无问题无界界对偶问对偶问题题原问题原问题无界无界如:如:(原)(原)无可无可行解行解(对)(对)对偶问题的性质推论推论:在一对对偶问题在一对对偶问题(P)(P)和和(D)(D)中,若其中一个问题可行但目中,若其中一个问题可行但目标函数无界,则另一个问题不标函数无界,则另一个问题不可行;反之不成立。这也是对可行;反之不成立。这也是对偶问题的无界性。偶问题的无界性。推推论论:在在一一对对对对偶偶问问题题(P P)和和(D D)中中,若若一一个个可可行行(如如P P),而而另另一一个个不不可可行行(如如D D),则则该
9、该可可行行的的问题无界。问题无界。例例1 1 试估计它们目标函数的界,并验证弱对偶性原理。试估计它们目标函数的界,并验证弱对偶性原理。(P P)对偶问题的性质解:解:(D D)由观察可知:由观察可知:,分别是(,分别是(P P)和)和(D D)的可行解。)的可行解。Z Z=10=10,W W=40=40,故有,故有 ,弱对偶定理成立。由推论弱对偶定理成立。由推论可知,可知,W W 的最小值不的最小值不能小于能小于1010,Z Z 的最大值不能超过的最大值不能超过4040。对偶问题的性质(P)(D)考虑考虑4.对偶定理 若(P)有最优解,则(D)也有最优解,且最优值相同。证:对(证:对(P)增加
10、松弛变量,)增加松弛变量,化为化为设其最优基为设其最优基为B,终表为,终表为其检验数为其检验数为问题:(1)由对偶定理可知,对偶问题最优解的表达式 Y*=?(2)求求Y*是否有必要重新求解(是否有必要重新求解(D)?)?CBB-1不必。可以从原问题(不必。可以从原问题(P)的单纯形终表获得。的单纯形终表获得。例如,已知例如,已知的终表为的终表为请指出其对偶问题的最优解和最优值。请指出其对偶问题的最优解和最优值。5.互补松弛定理y1yiymym+1ym+jyn+mx1xjxnxn+1xn+ixn+m对偶问题的变量对偶问题的变量对偶问题的松弛变量对偶问题的松弛变量原始问题的变量原始问题的变量原始问
11、题的松弛变量原始问题的松弛变量xjym+j=0yixn+i=0(i=1,2,m;j=1,2,n)在一对变量中,其中一个大于在一对变量中,其中一个大于0,另一个一定等于,另一个一定等于0直观上直观上定定义义:在在一一对对P和和D中中,若若P的的某某个个约约束束条条件件的的右右端端项项常常数数bi增增加加一一个个单单位位时时,所所引引起起的的目目标标函函数数最最优优值值Z*的的改改变变量量y*i称称为为第第i个个约约束条件的影子价格,又称为边际价格。束条件的影子价格,又称为边际价格。影子价格CCBCN0CBXBbXBXNXSCBXBB-1bIB-1NB-1ZCBB-1b0CNCBB-1NCBB-1
12、设:设:B是问题是问题P的最优基,由前式可知,的最优基,由前式可知,Z*=CBB-1b=Y*b=y*1b1+y*2b2+y*Ibi+y*mbm当当bi变为变为bi+1时(其余右端项不变,也不影响时(其余右端项不变,也不影响B)影子价格目标函数最优值变为:目标函数最优值变为:Z*=y*1b1+y*2b2+y*I(bi+1)+y*mbmZ*=Z*Z*=y*i也可以写成:也可以写成:即即y*i表示表示Z*对对bi的变化率。的变化率。其其经经济济意意义义是是:在在其其它它条条件件不不变变的的情情况况下下,单单位位资资源源变变化化所所引引起起的的目目标标函函数数的的最最优优值值的的变变化化。即即对对偶偶
13、变量变量yi就是第就是第i个约束条件的影子价格。个约束条件的影子价格。也也可可以以理理解解为为目目标标函函数数最最优优值值对对资资源源的的一一阶阶偏偏导导数(但问题中所有其它数据都保持不变)。数(但问题中所有其它数据都保持不变)。若若第第i种种资资源源的的单单位位市市场场价价格格为为mi,当当yimi时时,企企业业愿愿意意购购进进这这种种资资源源,单单位位纯纯利利为为yimi,则则有有利利可可图图;如如果果yi0。例11:在例1(煤电油例)中,其单纯形终表如下:(3)甲产品的价格在何范围内变化时,现最优解不变?)甲产品的价格在何范围内变化时,现最优解不变?解:甲产品的价格解:甲产品的价格c1是
14、基变量的价格系数。是基变量的价格系数。例例11:在例:在例1(煤电油例)中,其单纯形终表如下:(煤电油例)中,其单纯形终表如下:(4)若现又考虑一新产品丙,其资源单耗为)若现又考虑一新产品丙,其资源单耗为10,2,5,售价为售价为6.5,问该产品是否可投产?,问该产品是否可投产?故丙产品可以投产。故丙产品可以投产。首先将线性规划与经济问题联系起来的是首先将线性规划与经济问题联系起来的是T.G.Koopman(库普曼)和(库普曼)和L.V.Kamtorovich(康脱罗维奇)(康脱罗维奇)二人因此而共同分享了二人因此而共同分享了19751975年的第年的第7 7届诺贝尔经届诺贝尔经济学奖。济学奖
15、。求解线性规划的计算机软件举例LINDO、EXCELLINDO可以从下面的网址下载:可以从下面的网址下载:WWW.LLINDOLINDO由美国芝加哥大学开发,可求解线性规划和线性由美国芝加哥大学开发,可求解线性规划和线性整数规划等。其可按自然格式输入模型,使用方便。整数规划等。其可按自然格式输入模型,使用方便。输入例输入例:MAX2X+3Y?ST?4X+9Y9?7X+6Y13?END:GODORANGE(SENSITIVITY)ANALYSIS??Y可用可用HELP命命令令得到帮助。得到帮助。计算结果说明:计算结果说明:REDUCEDCOST为使该变量进基,其价格系数至少应为使该变量进基,其价
16、格系数至少应增加的数值;增加的数值;SLACKORSUPLUS松弛或剩余变量;松弛或剩余变量;DUALPRICES影子价格影子价格;ALLOWABLEINCREASE灵敏度分析中可使最优基不灵敏度分析中可使最优基不变的系数可增量之上界;变的系数可增量之上界;ALLOWABLEDECREASE灵敏度分析中可使最优基不灵敏度分析中可使最优基不变的系数可减量之上界;变的系数可减量之上界;使用EXCEL求解线性规划-进入进入EXCELEXCEL后先在表格的第一列输入变量、约束和目标的名后先在表格的第一列输入变量、约束和目标的名 称,在后面某一列对应位置输入约束和目标的表达式称,在后面某一列对应位置输入
17、约束和目标的表达式(前加前加 等号提示符等号提示符);-然后在工具中调用规划求解,按提示操作。然后在工具中调用规划求解,按提示操作。(可参看关于使用(可参看关于使用EXCELEXCEL的书,如谢国锋等,的书,如谢国锋等,EXCEL2000EXCEL2000中中文版入门与提高文版入门与提高,清华大学出版社,清华大学出版社,19991999)5 5 整数规划整数规划Integer Programming(Integer Programming(简称简称IP)IP)一、整数规划的一般模型一、整数规划的一般模型LP:max z=CX AX=b X0IP:max z=CX AX=b X0 X为整数例一、
18、合理下料问题例一、合理下料问题设用某型号的圆钢下零件设用某型号的圆钢下零件A A1 1,A A2 2,A,Am m 的毛坯。在一根的毛坯。在一根圆钢上下料的方式有圆钢上下料的方式有B B1 1,B,B2 2,B Bn n 种,每种下料方式可以得种,每种下料方式可以得到各种零件的毛坯数以及每种零件的需要量,如表所示。怎到各种零件的毛坯数以及每种零件的需要量,如表所示。怎样安排下料方式,使得即满足需要,所用的原材料又最少?样安排下料方式,使得即满足需要,所用的原材料又最少?零件 方 个数 式零件零 件毛坯数整数规划的模型设:设:xj表示用表示用Bj(j=1.2n)种方式下料根数种方式下料根数模型:
19、模型:整数规划的模型纯纯整整数数规规划划:所所有有决决策策变变量量要要求求取取非非负负整整数数(这这时时引引进进的的松松弛弛变变量量和和剩剩余余变变量量可可以以不不要要求求取整数)。取整数)。全全整整数数规规划划:除除所所有有决决策策变变量量要要求求取取非非负负整整数数外外,系系数数aij和和常常数数bi也也要要求求取取整整数数(这这时时引引进进的松弛变量和剩余变量也必须是整数)。的松弛变量和剩余变量也必须是整数)。混混合合整整数数规规划划:只只有有一一部部分分的的决决策策变变量量要要求求取取非负整数,另一部分可以取非负实数。非负整数,另一部分可以取非负实数。01整整数数规规划划:所所有有决决
20、策策变变量量只只能能取取0或或1两两个整数。个整数。整数规划的数学模型从从数数学学模模型型上上看看整整数数规规划划似似乎乎是是线线性性规规划划的的一一种种特特殊殊形形式式,求求解解只只需需在在线线性性规规划划的的基基础础上上,通通过过舍舍入入取取整整,寻寻求求满满足足整整数数要要求求的的解解即即可可。但但实实际际上上两两者者却却有有很很大大的的不不同同,通通过过舍舍入入得得到到的的解解(整整数数)也也不不一一定定就就是是最最优优解解,有有时时甚至不能保证所得倒的解是整数可行解。甚至不能保证所得倒的解是整数可行解。整数规划与线性规划的关系例:设整数规划问题如下例:设整数规划问题如下 首先不考虑整
21、数约束,得到线性规划问题(一般首先不考虑整数约束,得到线性规划问题(一般称为松弛问题)。称为松弛问题)。整数规划与线性规划的关系且为整数 用图解法求出最优解用图解法求出最优解x x1 13/2,3/2,x x2 2=10/3=10/3且有且有Z=29/6Z=29/6x1x233(3/2,10/3)现现求求整整数数解解(最最优优解解):如如用用“舍舍入入取取整整法法”可可得得到到4个个点点即即(1,3),(2,3),(1,4),(2,4)。显显然然,它它们们都都不不可可能能是是整整数数规规划划的的最优解。最优解。按按整整数数规规划划约约束束条条件件,其其可可行行解解肯肯定定在在线线性性规规划划问
22、问题题的的可可行行域域内内且且为为整整数数点点。故故整整数数规规划划问问题题的的可可行行解解集集是是一一个个有有限限集集,如图所示。如图所示。整数规划与线性规划的关系因因此此,可可将将集集合合内内的的整整数数点点一一一一找找出出,其其最最大目标函数的值为最优解,此法为完全枚举法。大目标函数的值为最优解,此法为完全枚举法。如如上上例例:其其中中(2 2,2 2)(3 3,1 1)点点为为最最大大值值,Z=4Z=4。目前,常用的求解整数规划的方法有:目前,常用的求解整数规划的方法有:分支定界法和割平面法;分支定界法和割平面法;对对于于特特别别的的0 01 1规规划划问问题题采采用用隐隐枚枚举举法法
23、和和匈匈牙利法。牙利法。整数规划与线性规划的关系考虑纯整数问题:考虑纯整数问题:整数问题的松弛问题:整数问题的松弛问题:分枝定界法1、先先不不考考虑虑整整数数约约束束,解解(IP)的的松松弛弛问问题题(LP),可可能得到以下情况之一:能得到以下情况之一:.若若(LP)没没有有可可行行解解,则则(IP)也也没没有有可可行行解解,停停止计算。止计算。.若若(LP)有有最最优优解解,并并符符合合(IP)的的整整数数条条件件,则则(LP)的最优解即为的最优解即为(IP)的最优解,停止计算。的最优解,停止计算。.若若(LP)有有最最优优解解,但但不不符符合合(IP)的的整整数数条条件件,转入下一步。为讨
24、论方便,设转入下一步。为讨论方便,设(LP)的最优解为:的最优解为:分枝定界法2、定界:、定界:记记(IP)的的目目标标函函数数最最优优值值为为Z*,以以Z(0)作作为为Z*的的上上界界,记记为为Z(0)。再再用用观观察察法法找找的的一一个个整整数数可可行行解解X,并并以以其其相相应应的的目目标标函函数数值值Z作作为为Z*的的下下界界,记记为为ZZ,也可以令,也可以令Z,则有,则有:Z Z*3、分枝:、分枝:在在(LP)的的最最优优解解X(0)中中,任任选选一一个个不不符符合合整整数数条条件件的的变变量量,例例如如xr=br(不不为为整整数数),以以br 表表示示不不超过超过br的最大整数。构
25、造两个约束条件的最大整数。构造两个约束条件xrbr 和和xrbr 1分枝定界法将将这这两两个个约约束束条条件件分分别别加加入入问问题题(IP),形形成成两两个个子子问问题题(IP1)和和(IP2),再再解解这这两两个个问问题题的的松松弛弛问问题题(LP1)和和(LP2)。4、修改上、下界:按照以下两点规则进行。、修改上、下界:按照以下两点规则进行。.在在各各分分枝枝问问题题中中,找找出出目目标标函函数数值值最最大大者者作作为为新的上界;新的上界;.从从已已符符合合整整数数条条件件的的分分枝枝中中,找找出出目目标标函函数数值值最大者作为新的下界。最大者作为新的下界。5、比较与剪枝、比较与剪枝:各
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 运筹学 线性规划

限制150内