线性规划和单纯形法课件.ppt
《线性规划和单纯形法课件.ppt》由会员分享,可在线阅读,更多相关《线性规划和单纯形法课件.ppt(76页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、目目 录录p线性规划实例与模型线性规划实例与模型p线性规划的图解法线性规划的图解法p单纯形法原理单纯形法原理p改进单纯形法改进单纯形法p应用应用目目 录录p线性规划实例与模型线性规划实例与模型实用举例实用举例 某公司通过市场调研,决定生产高中档新型拉杆箱。某分销商决定买进该公司3个月内的全部产品。拉杆箱生产需经过原材料剪裁、缝合、定型、检验和包装4过程。 通过分析生产过程,得出:生产中档拉杆箱需要用7/10小时剪裁、5/10小时缝合、1小时定型、1/10小时检验包装;生产高档拉杆箱则需用1小时剪裁、5/6小时缝合、2/3小时定型、1/4小时检验包装。由于公司生产能力有限,3月内各部的最大生产时
2、间为剪裁部630小时、缝合部600小时、定型部708小时、检验包装部135小时。 通过市场调研部和会计部的调查核算得出结论:生产中档拉杆箱的利润是10元,高档拉杆箱的利润是9元。公司应各生产多少中档和高档拉杆箱才能使公司利润最大?实用举例实用举例 某公司通过市场调研,决定生产高中档新型拉杆箱。某分销商决定买进该公司3个月内的全部产品。拉杆箱生产需经过原材料剪裁、缝合、定型、检验和包装4过程。 通过分析生产过程,得出:生产中档拉杆箱需要用7/10小时剪裁、5/10小时缝合、1小时定型、1/10小时检验包装;生产高档拉杆箱则需用1小时剪裁、5/6小时缝合、2/3小时定型、1/4小时检验包装。由于公
3、司生产能力有限,3月内各部的最大生产时间为剪裁部630小时、缝合部600小时、定型部708小时、检验包装部135小时。 通过市场调研部和会计部的调查核算得出结论:生产中档拉杆箱的利润是10元,高档拉杆箱的利润是9元。公司应各生产多少中档和高档拉杆箱才能使公司利润最大?可以通过线性规划求解可以通过线性规划求解!线性规划模型的建立线性规划模型的建立 设生产中、高档拉杆箱数量分别为: 称为决策变量。0, 0 135 708 600 630 . .910 max21241110123212651212110721xxxxxxxxxxtsxx目标函数目标函数约束条件约束条件21,xx一般线性规划模型Ma
4、x z = c1x1 + c2x2 + + cnxns.t. a11x1+ a12x2 + + a1nxn = b1 ( 0) a21x1+ a22x2 + + a2nxn = b2 ( 0) : am1x1+ am2x2 + + amnxn= bm( 0) x1,x2 , ,xn 0s.t. 为约束限制(Subject to)的缩写(LP)x1.xnb1.bma11 a1nam1 amnx =b =A = c =0Xb),(AXt . scxzmin)max( 或其中c=(c1,c2,cn),称为价值系数向量;mbbbb21称为资源限制向量 X=(x1,x2,xn)T 称为决策变量向量aaa
5、aaaaaamnmmnnA212222111211称为技术系数矩阵(也称消耗系数矩阵)一般线性规划模型线性规划模型的标准形式线性规划模型的标准形式Max Z = c1x1 + c2x2 + + cnxnSubject to (s.t.)a11x1 + a12x2 + + a1nxn = b1 a21x1 + a22x2 + + a2nxn = b2 am1x1 + am2x2 + + amnxn= bm x x1 1 0, 0, x x2 2 0, 0, , , x xn n 0 为了讨论方便,把最大化、等式约束型的线性规划称为线为了讨论方便,把最大化、等式约束型的线性规划称为线性规划的标准型
6、,即性规划的标准型,即标准化标准化原问题原问题标准化方法标准化方法目标函数目标函数Max f(x)Max f(x)Min f(x)Max f(x)约束条件约束条件引入松弛变量和人工变量引入松弛变量, 不变变量变量 不变对某个njijijbxa1 i对某个njijijbxa1 i对某个njijijbxa1 人工变量松弛变量njijijbxa1 松弛变量0 iix对某个0 iix对某个0iiixxx,其中令0,iiiiixxxxx,其中令njijijbxa1 i对某个对某个ixi 是任意的 线性规划模型的标准化线性规划模型的标准化 例例: :将如下线性规划模型标准化:将如下线性规划模型标准化:mi
7、n z= x1 + 5 x2 - 8 x3 - x4s.t. 2 x1 - 3 x2 + x3 + x4 20 x1+ 2 x2 + 3 x3 - x4 30 2x2 + 2 x3 +3 x4- 50 x1 , x3 , x4 0,x2取值无约束。取值无约束。 max z1=-x1-5 x2+8 x3 +x4s.t. 2 x1 - 3 x2 + x3 + x4 20 x1+ 2 x2 + 3 x3 - x4 30 2x2 + 2 x3 +3 x4- 50 x1 , x3 , x4 0,x2取值无约束。取值无约束。 目标优化标准化目标优化标准化线性规划模型的标准化Max z2 =-x1-5y2+
8、5y3+8x3+x4 s.t. 2 x1 - 3 y2+3y3 + x3 + x4 20 -x1- 2 y2 +2y3 - 3 x3 + x4 -30 -2y2 +2y3 - 2 x3 -3 x4 50 x1 , x3 , x4 ,y2,y3 0 决策变量的标准化决策变量的标准化: y2 - y3 = x2max z1=-x1-5 x2+8 x3 +x4s.t. 2 x1 - 3 x2 + x3 + x4 20 x1+ 2 x2 + 3 x3 - x4 30 2x2 + 2 x3 +3 x4- 50 x1 , x3 , x4 0,x2取值无约束取值无约束 线性规划模型的标准化Max z2 =-
9、x1-5y2+5y3+8x3+x4 s.t. 2 x1 - 3 y2+3y3 + x3 + x4 +x5 =20 -x1- 2 y2 +2y3 - 3 x3 + x4 +x6 = -30 -2y2 +2y3 - 2 x3 -3 x4 +x7 = 50 x1 , x3 , x4 , x5, x6, x7, y2, y3 0约束关系标准化约束关系标准化Max z2 =-x1-5y2+5y3+8x3+x4 s.t. 2 x1 - 3 y2+3y3 + x3 + x4 20 -x1- 2 y2 +2y3 - 3 x3 + x4 -30 -2y2 +2y3 - 2 x3 -3 x4 50 x1 , x3
10、 , x4 ,y2,y3 0 p可行解可行解:满足所有约束条件的解x=(x1,x2,.,xn),所有可行解的集合称为可行域可行域。p基基:设A是约束方程组的mn阶系数矩阵,秩为m, 是A中阶非奇异子矩阵(即 ),则称是线性规划问题的一个基矩阵基矩阵,简称基基。B中的列向量称为基向量基向量,与基向量对应的变量x称为基变量基变量,其它变量称为非基变量非基变量。 p基本解基本解:令非基变量值为0,由Ax=b可求出一个解,这个解x称为基本解基本解。p基本可行解基本可行解:满足非负条件的基本解称为基本可行解基本可行解。p最优解最优解:使目标函数达到最优值的可行解称为最优解最优解。 线性规划的解线性规划的
11、解 ),(1mPPB| 0B 可行解、基本解、基本可行解的关系可行解可行解基本解基本解基本可行解基本可行解目目 录录p线性规划实例与模型线性规划实例与模型p线性规划的图解法与基本性质线性规划的图解法与基本性质线性规划的图解法线性规划的图解法 当只有两个决策变量当只有两个决策变量时,可用图解法求解!时,可用图解法求解!例例:用图解法求解以下线形规划问题。 max z=4x1+3x2 s.t. x16 2x28 2x1+3x218 x10,x20 x1 x2 L3:2x1+3x2=18可行域可行域L1:x1=6L2:x2=4最优解B4x1+3x2解的特殊情况解的特殊情况多个最优解多个最优解解的特殊
12、情况解的特殊情况无可行解无可行解解的特殊情况解的特殊情况无界解无界解线性规划的基本性质线性规划的基本性质 X可行域内部的点 可行解? 是是 最优解? 不不 若线性规划有最若线性规划有最优解,则最优解必在可优解,则最优解必在可行域的顶点上达到。行域的顶点上达到。凸集的概念凸集的概念 p凸集是线性规划中一个很重要的概念,凸集的一般定义为:凸集是线性规划中一个很重要的概念,凸集的一般定义为:如果在集合如果在集合C C中任意取两个点中任意取两个点x x1 1,x x2 2,其连线上的所有点也,其连线上的所有点也都在集合都在集合C C中,则称集合中,则称集合C C为凸集合。用数学解析式表达为:为凸集合。
13、用数学解析式表达为:对任意对任意 ,均有,均有 ,则称,则称C C是是凸集合。凸集合。Cxx21,Cxx21)1 () 10 (非凸集非凸集非凸集非凸集凸集凸集线性规划的基本性质 定理定理2.12.1:线性规划的可行域:线性规划的可行域: 是凸集(凸多面体)。是凸集(凸多面体)。0),(,|1nixxxxbAxxD 引理引理2.12.1:线性规划的可行解线性规划的可行解 为基本可行解的为基本可行解的充分必要条件是充分必要条件是x x的正分量所对应的系数列向量是线性无关的,的正分量所对应的系数列向量是线性无关的,即每个正分量都是一个基变量即每个正分量都是一个基变量。 Tnxxx),(1定理定理2
14、.22.2:线性规划问题的基本可行解线性规划问题的基本可行解x x对应于可行域的顶点对应于可行域的顶点 定理定理2.32.3:若线性规划有最优解,则最优解必在可行域的顶点上若线性规划有最优解,则最优解必在可行域的顶点上达到。达到。目目 录录p线性规划实例与模型线性规划实例与模型p线性规划的图解法线性规划的图解法p单纯形法原理单纯形法原理一、初始基本可行解的确定一、初始基本可行解的确定p考虑如下形式的线性规划问题考虑如下形式的线性规划问题max z=c1x1+c2x2+.+cnxn s.t. a11x1+a12x2+a1nxnb1 a21x1+a22x2+a2nxnb2 (2.17) . am1
15、x1+am2x2+ +amnxnbm x1,x2,.,xn0 (2.18)对每个不等式约束引入一个非负松弛变量,得标准形式:对每个不等式约束引入一个非负松弛变量,得标准形式:max z=c1x1+c2x2+.+cn xns.t. a11x1+a1nxn+xn+1 =b1 a21x1+a2nxn +xn+2 =b2 (2.19) . am1x1+amn xn +xn+m =bm x1,x2,.,xn,xn+1,xn+m0100010001),(1mnnPPB), 1(mibxiin是线性规划是线性规划(2.19)(2.19)的一个可行基。令的一个可行基。令 021nxxx得得由此得到一个初始基本
16、可行解为由此得到一个初始基本可行解为TmTmnnnbbxxxx), 0 , 0(), 0 , 0(121二、最优性检验二、最优性检验 p定理定理2.42.4: 在最大化问题中,对于某个基本可行解,如果所有 的 ,则这个基本可行解是最优解。 在最小化问题中,对于某个基本可行解,如果所有 的 ,则这个基本可行解是最优解。0jjzc0jjzc单纯形法求解步骤单纯形法求解步骤p将所给问题化为标准形将所给问题化为标准形p找出一个初始可行基找出一个初始可行基, ,建立初始单纯形表建立初始单纯形表p检查所有检验数检查所有检验数( (若全为非负若全为非负, ,则已得到最优解则已得到最优解, ,计计算停止算停止
17、. .否则继续下一步否则继续下一步) )p考察是否无解考察是否无解( (若是若是, ,计算停止计算停止, ,否则继续下一步否则继续下一步) )p确定入基变量确定入基变量, ,出基变量出基变量p对初始单纯形表进行单纯形变换对初始单纯形表进行单纯形变换单纯形方法的一般步骤单纯形方法的一般步骤确定一个初始可行角点确定一个初始可行角点根据某种法则进行角点根据某种法则进行角点的最优性判断的最优性判断不是最优角点不是最优角点是最优角点是最优角点考察与当前角点相邻的考察与当前角点相邻的 “ “更好更好”的的一个角点,并置为当前角点。一个角点,并置为当前角点。根据某种法则进行根据某种法则进行LPLP的无的无界
18、或不可行判断界或不可行判断无界或不可行无界或不可行还不能做出判断还不能做出判断求求解解结结束束单纯形法举例单纯形法举例0,26232432.34max21212121xxxxxxxxtsZ解:引进松驰变量x3, x4,化为标准形得:0,262324 32.0034max43214 213214321xxxxxxxxxxxxxxtsZ)26,24, 0 , 0(),(4321xxxx 4 3 0 0 CBXBb x1 x2 x3 x4 b/y00 x3x42426 2 3 1 0 24/2 3 2 0 1 26/3 0 4 3 0 04=4(03+02);3=3(02+03) 4 3 0 0CB
19、XBb x1 x2 x3 x4 b/y04x3x120/326/3 0 5/ 3 1 -2/3 1 2/3 0 1/3 -104/3 0 1/3 0 -4/3jjzcjjzcicic 4 3 0 0CBXBb x1 x2 x3 x434x2x146 0 1 3/5 -2/5 1 0 -2/5 3/5 36 0 0 -1/5 -6/5表中最后一行的所有检验数均为非正数,表明目标函数已达到最大值,上述表为最优表。从表中可得到最优解为X=(x1,x2,x3,x4)=(6,4,0,0),最优值为Z=36icjjzc 4 3 0 0CBXBb x1 x2 x3 x4 b/y04x3x120/326/3
20、0 5/ 3 1 -2/3 4 1 2/3 0 1/3 13 -104/3 0 1/3 0 -4/3icjjzc单纯形法的进一步讨论人工变量法单纯形法的进一步讨论人工变量法p人工变量法:p前面讨论了在标准型中系数矩阵有单位矩阵,很容易确定一组基可行解。在实际问题中有些模型并不含有单位矩阵,为了得到一组基向量和初基可行解,在约束条件的等式左端加一组虚拟变量,得到一组基变量。这种人为加的变量称为人工变量,构成的可行基称为人工基,用大M法或两阶段法求解,这种用人工变量作桥梁的求解方法称为人工变量法。单纯形法的进一步讨论人工变量法单纯形法的进一步讨论人工变量法p例1.10 用大M法解下列线性规划 01
21、2210243423max321321321321321xxxxxxxxxxxxxxxZ、解:首先将数学模型化为标准形式解:首先将数学模型化为标准形式 5 , 2 , 1, 012210243423max32153214321321jxxxxxxxxxxxxxxxZj系数矩阵中不存在系数矩阵中不存在单位矩阵,无法建单位矩阵,无法建立初始单纯形表。立初始单纯形表。单纯形法的进一步讨论人工变量法单纯形法的进一步讨论人工变量法p故人为添加两个单位向量,得到人工变量单纯形法数学模故人为添加两个单位向量,得到人工变量单纯形法数学模型:型:7 , 2 , 1, 012210243423max7321532
22、16432176321jxxxxxxxxxxxxxxMxMxxxxZj其中:其中:M是一个很大的抽象的数,不需要给出具体的数值,是一个很大的抽象的数,不需要给出具体的数值,可以理解为它能大于给定的任何一个确定数值;再用前面介可以理解为它能大于给定的任何一个确定数值;再用前面介绍的单纯形法求解该模型,计算结果见下表。绍的单纯形法求解该模型,计算结果见下表。 单纯形法的进一步讨论人工变量法单纯形法的进一步讨论人工变量法cj32-100-M-MCBXBbx1x2x3x4x5x6x7i-Mx64-431-101040 x5101-1201005-Mx712-21000113-2M2+M-1+2M-M0
23、000 x63-650-1013/5-Mx58-3300108/3-1x312-210005-6M5M0-M002x23/56/5101/50-Mx531/53/5003/5131/3-1x311/52/5012/505 00002x213010123x131/310015/3-1x319/300102/3000-5-25/3j j j j 单纯形法的进一步讨论人工变量法单纯形法的进一步讨论人工变量法解的判别:解的判别:1)唯一最优解判别:最优表中所有非基变量的检验数非零)唯一最优解判别:最优表中所有非基变量的检验数非零,则线则线 规划具有唯一最优解。(检验数为负数或规划具有唯一最优解。(检验
24、数为负数或0)2)多重最优解判别:最优表中存在非基变量的检验数为零)多重最优解判别:最优表中存在非基变量的检验数为零,则线则性规划具有无穷多最优解(检验数为负数或则线则性规划具有无穷多最优解(检验数为负数或0)。)。3)无界解判别:某个检验数)无界解判别:某个检验数0且且aik(i=1,2,m)则)则线性规划具有无界解。线性规划具有无界解。4)无可行解的判断:当用大)无可行解的判断:当用大M单纯形法计算得到最优解并单纯形法计算得到最优解并且存在且存在Ri0时,则表明原线性规划无可行解。时,则表明原线性规划无可行解。5)退化解的判别:)退化解的判别:存在某个基变量为零的基本可行解。存在某个基变量
25、为零的基本可行解。大大M M法法0,12 3421123min32131321321321xxxxxxxxxxxxxxz基本思想:基本思想:在约束条件中增加人工变量,同时修改目标函数,对目标函数加上具有很大正系数M的惩罚项,为使人工变量不影响目标函数取值,在迭代过程中,应把人工变量从基变量中退出,使其成为非基变量。0,1 2 3 4211 2)(3min72173165321432176321xxxxxxxxxxxxxxxxxMxxxz 其中,M是很大的正数,同时增加两个人工变量。 不容易找到初始可行解找初始可行解找初始可行解icBcBxjjzc11-300MMbi/yibx1X2x3x4x5
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 线性规划 单纯 课件
限制150内