《运筹学》复习参考资料.docx

1、运筹学复习参考资料 线性规划问题的求解重要算法：图解法、单纯形迭代、大M法单纯形迭代、对偶问题、表上作业法（找初始可行解：西北角法，最小元素法；最优性检验：闭回路法，位势法；）、目标规划：图解法、整数规划：分支定界法（次重点），匈牙利法（重点）、 动态规划问题的求解重要算法：图上标号法 网络分析问题的求解重要算法：破圈法、TP标号法、寻求网络最大流的标号法第一部分 线性规划问题的求解一、两个变量的线性规划问题的图解法：概念准备：定义：满足所有约束条件的解为可行解；可行解的全体称为可行（解）域。定义：达到目标的可行解为最优解。图解法：图解法采用直角坐标求解：x1横轴；x2竖轴。1、将约束条件（取

2、等号）用直线绘出；2、确定可行解域；3、绘出目标函数的图形（等值线），确定它向最优解的移动方向；注：求极大值沿价值系数向量的正向移动；求极小值沿价值系数向量的反向移动。4、确定最优解及目标函数值。参考例题：（只要求下面这些有唯一最优解的类型）例1：某厂生产甲、乙两种产品，这两种产品均需在A、B、C三种不同的设备上加工，每种产品在不同设备上加工所需的工时不同，这些产品销售后所能获得利润以及这三种加工设备因各种条件限制所能使用的有效加工总时数如下表所示：品产耗消备设 A B C利润（万元）甲乙3 5 99 5 37030有效总工时540 450 720问：该厂应如何组织生产，即生产多少甲、乙产品使

3、得该厂的总利润为最大？（此题也可用“”或化“”用大M法求解）解：设x1、x2为生产甲、乙产品的数量。 max z = 70x1+30x2 s.t. 、 可行解域为oabcd0，最优解为b点。由方程组 解出x1=75，x2=15X*=（75，15）Tmax z =Z*= 7075+3015=5700例2：用图解法求解 max z = 6x1+4x2 s.t. 、 可行解域为oabcd0，最优解为b点。由方程组 解出x1=2，x2=6X*=（2，6）Tmax z = 62+46=36例3：用图解法求解 min z =3x1+x2s.t. 、解：可行解域为bcdefb，最优解为b点。由方程组 解出x

4、1=4，x2=X*=（4，）Tmin z =34+=11二、标准型线性规划问题的单纯形解法：一般思路：1、用简单易行的方法获得初始基本可行解；2、对上述解进行检验，检验其是否为最优解，若是，停止迭代，否则转入3；3、根据L规则确定改进解的方向；4、根据可能改进的方向进行迭代得到新的解；5、根据检验规则对新解进行检验，若是最优解，则停止迭代，否则转入3，直至最优解。具体做法（可化归标准型的情况）：设已知max z = c1x1+ c2x2+ cnxns.t.对第i个方程加入松弛变量xn+i，i =1，2，m，得到列表计算，格式、算法如下：CBXBbc1c2cn+mLx1x2xn+mcn+1xn+

5、1b1a11a12a1 n+mc n+2xn+2b2a21a22a2 n++mxn+mbnam1am2am n+mz1z2zn+m12n+m注： zj =cn+1 a1j+ cn+2 a2j + cn+m amj=，（j=1，2，n+m）j =cjzj ,当j 0时，当前解最优。注：由maxj确定所对应的行的变量为“入基变量”；由L=确定所对应的行的变量为“出基变量”，行、列交叉处为主元素，迭代时要求将主元素变为1，此列其余元素变为0。例1：用单纯形法求解（本题即是本资料P2“”例1的单纯形解法；也可化“”求解）max z =70x1+30x2s.t.解：加入松弛变量x3，x4，x5，得到等效

6、的标准模型：max z =70x1+30x2+0 x3+0 x4+0 x5s.t.列表计算如下：CBXBb7030000Lx1x2x3x4x50x354039100540/3 =1800x445055010450/5 =900x5720（9）3001720/9 =800000070300000x33000810- 1/3300/8 =37.50x4500（10/3）01 - 5/950/10/3 =1570x1801 1/300 1/980/1/3 =2407070/30070/9020/30070/90x318000112/5130x215010 3/10- 1/670x175100- 1/

7、10 1/6570070300220/3000-220/3X*=（75，15，180，0，0）Tmax z =7075+3015=5700例2：用单纯形法求解max z =7x1+12x2s.t.解：加入松弛变量x3，x4，x5，得到等效的标准模型：max z =7x1+12x2+0 x3+0 x4+0 x5s.t.列表计算如下：CBXBb712000Lx1x2x3x4x50x336094100360/4 =900x420045010200/5 =400x53003（10）001300/10 =30000007120000x324078/10010- 2/5240/78/10 =2400/78

8、0x450（5/2）001- 1/250/5/2 =2012x2303/10100 1/1030/3/10 =10018/512006/517/50006/50x38400178/2529/257x1201002/5- 1/512x224010 3/254/28428712034/2511/3500034/2511/35X*=（20，24，84，0，0）Tmax z =720+1224=428三、非标准型线性规划问题的解法：1、一般地，对于约束条件组：若为“”，则加松弛变量，使方程成为“”；若为“”，则减松弛变量，使方程成为“”。我们在前面标准型中是规定目标函数求极大值。如果在实际问题中遇到的

9、是求极小值，则为非标准型。可作如下处理：由目标函数min z=变成等价的目标函数max（z）=令z=z/，min z=max z/2、等式约束大M法：通过加人工变量的方法，构造人造基，从而产生初始可行基。人工变量的价值系数为M，M是很大的正数，从原理上理解又称为“惩罚系数”。（课本P29）类型一：目标函数仍为max z，约束条件组与。例1：max z =3x1+5x2s.t.解：加入松弛变量x3，x4，得到等效的标准模型：max z =3x1+5x2s.t.其中第三个约束条件虽然是等式，但因无初始解，所以增加一个人工变量x5，得到： max z =3x1+5x2Mx5s.t. 单纯形表求解过程

10、如下：CBXBb3500MLx1x2x3x4x50x34（1）01004/1 =40x41202010Mx5183200118/3 =63M2M00M3M352M0003x14101000x4120201012/2 =6Mx560（2）3016/2 =332M33M0M0533M003x14101004/1 =40x4600（3）116/3 =25x23013/201/23/(2/3) =9/2359/205/2009/20M5/2305x121001/31/3x320011/31/3x260101/20363503/210003/2M1X*=（2，6，2，0）Tmax z =32+56=36

11、类型二：目标函数min z，约束条件组与。例2：用单纯形法求解min z =4x1+3x2s.t.解：减去松弛变量x3，x4，并化为等效的标准模型：max z/ =4x13x2s.t.增加人工变量x5、x6，得到：max z/ =4x13x2Mx5Mx6s.t单纯形表求解过程如下：CBXBb400MMLx1x2x3x4x5x6Mx5162（4）101016/4=4Mx61232010112/2=65M6MMMMM5M46M3MM003x241/211/401/404/1/2=8Mx64（2）01/211/214/2=22M3/233/4M/2MM/23/4M2M5/20M/23/4M3/43M

12、/203x23013/81/43/81/44x12101/41/21/41/217431/85/41/85/4001/85/4M1/8M5/4X*=（2，3，0，0）Tmin z =max z/ =（17）=17四、对偶问题的解法：什么是对偶问题？1、在资源一定的条件下，作出最大的贡献；2、完成给定的工作，所消耗的资源最少。引例（与本资料P2例1 “”、P7例1 “”同）：某工厂生产甲、乙两种产品，这些产品均需在A、B、C三种不同的设备上加工，每种产品在不同设备上加工时需要不同的工时，这些产品售后所能获得的利润值以及这三种加工设备因各种条件下所能使用的有效总工时数如下表：品产耗消备设 A B

13、C利润（万元）甲乙3 5 99 5 37030有效总工时540 450 720问：该厂应如何组织生产，即生产多少甲、乙产品使得该厂的总利润为最大？解：原问题设x1、x2为生产甲、乙产品的数量。max z = 70x1+30x2s.t.将这个原问题化为它的对偶问题设y1、y2、y2分别为设备A、B、C单位工时数的加工费。min w = 540y1+450y2+720y3s.t.用大M法，先化为等效的标准模型：max w/ =540y1450y2720y3s.t.增加人工变量y6、y7，得到：max z/ =540y1450y2720y3My6My7s.t大M法单纯形表求解过程如下：CBXBb54

14、045072000MMLy1y2y3y4y5y6y7My670359101070/3My730（9）53010130/9=10/312M10M12MMMMM12M54010M45012M720MM00My660010/3（8）11/311/360/8=2.5540y110/315/91/301/901/910/3/1/3=10-300+10/3M-8M180MM/3+60MM/3600-150+10/3M8M-540MM/3600M/3+60720y315/205/1211/81/241/81/2415/2/5/12=18540y15/61（5/12）01/241/81/241/85/6/5/

15、12=2540572720135/2475/12135/275/201250135/2475/12135/2M75/2M720450y320/31011/61/61/61/6y2212/5101/103/101/103/1057003604507207515751518000751575M15M该对偶问题的最优解是y*=（0，2，0，0）T最优目标函数值min w =（5700）=5700五、运输规划问题：运输规划问题的特殊解法“表上作业法”解题步骤：1、找出初始调运方案。即在(mn)产销平衡表上给出m+n-1个数字格。（最小元素法）2、（对空格）求检验数。判别是否达到最优解。如已是最优解，则

16、停止计算，否则转到下一步。（闭回路法）3、对方案进行改善，找出新的调运方案。（根据检验结果选择入基变量，用表上闭回路法调整即迭代计算，得新的基本可行解）4、重复2、3，再检验、再迭代，直到求得最优调运方案。类型一：供求平衡的运输规划问题（又称“供需平衡”、“产销平衡”）引例：某钢铁公司有三个铁矿和四个炼铁厂，铁矿的年产矿石量分别为100万吨、80万吨和50万吨，炼铁厂年需矿石量分别为50万吨、70万吨、80万吨和30万吨，这三个铁矿与四个炼铁厂的距离如下：矿铁厂铁离距炼 B1 B2 B3 B4A1A2A315 20 3 3070 8 14 2012 3 20 15问：该公司应如何组织运输，既满

17、足各炼铁厂需要，又使总的运输费用为最小（按吨.公里计）？解：用“表上作业法”求解。1 用最低费用法（最小元素法）求此问题的初始基础可行解：地产用费地销B1B2B3B4产量SiA1152067330651002080A270814442080302030A312533203325105050销量dj50708030 230230初始方案：B2203030B1B3A22080B1B3A1B250A3Z=1520+380+7030+820+2030+350=3550（吨.公里）2 的初始可行解进行迭代（表上闭回路法），求最优解：地产用费地销B1B2B3B4产量SiA11520143301210020

18、80A27053814920805030A312320202510503020销量dj50708030 230230用表上闭回路法调整后，从上表可看出，所有检验数=0，已得最优解。最优方案：3020B1B2A35030B2B4A22080B1B3A1Z=1520+380+850+2030+1230+320=1960（吨.公里）解法分析：如何求检验数并由此确定入基变量？有数字的空格称为“基格”、打的空格称为“空格”，标号为偶数的顶点称为偶点、标号为奇数的顶点称为奇点，出发点算0故为偶点。找出所有空格的闭回路后计算它们的检验数，必须0，才得到最优解。否则，应选所有中正的最大者对应的变量xj为入基变

19、量进行迭代（调整）。检验后调整运输方案的办法是：在空格的闭回路中所有的偶点均加上奇点中的最小运量，所有的奇点均减去奇点中的最小运量。重复以上两步，再检验、再调整，直到求得最优运输方案。求下面运输问题的最小值解：12341311310721923437410593656解：由最小元素法得到初始解：v1=2v2=9v3=3v4=101934u1=01311310743u2=-121923431u3=-53741059633656则：，最小值为-6，非基变量为，闭回路，最大调整量为1，得新解：，重新计算位势及影响系数，得下表：v1=8v2=9v3=3v4=101234u1=01311310752u2

20、=-721923431u3=-53741059633656，最小值为-5，非基变量为，闭回路，最大调整为2，得新解：重新计算位势及影响系数，得下表：v1=3v2=4v3=3v4=51234u1=01311310725u2=-221923413u3=03741059633656，此时，故当前解为最优解。最优解值为：。六、目标规划问题已知目标规划问题min z= p1 d1- + p2d2- + p3（5d3- + 3d4- ） + p4 d1+x1 +2x2 + d1- - d1+ =6x1 +2x2 + d2- - d2+ =9x1 -2x2 + d3- - d3+ =4x2 + d4- - d4+ =2x1 ，x20，di- , di+0 (i=1,2,3,4)分别用图解法和单纯形法求解x2d2+54d2-d4+d3-d1+3Ad4-d1-2d3+B10 2 4 6 8 x1七（1）、整数规划问题（分枝定界法）：用分枝定界法求解整数规划