线性规划和单纯形法幻灯片.ppt
《线性规划和单纯形法幻灯片.ppt》由会员分享,可在线阅读,更多相关《线性规划和单纯形法幻灯片.ppt(77页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、线性规划和单纯形法线性规划和单纯形法第1页,共77页,编辑于2022年,星期一目目 录录p线性规划实例与模型线性规划实例与模型p线性规划的图解法线性规划的图解法p单纯形法原理单纯形法原理p改进单纯形法改进单纯形法p应用应用第2页,共77页,编辑于2022年,星期一目目 录录p线性规划实例与模型线性规划实例与模型第3页,共77页,编辑于2022年,星期一实用举例实用举例 某公司通过市场调研,决定生产高中档新型拉杆箱。某分销商决定买进该公司3个月内的全部产品。拉杆箱生产需经过原材料剪裁、缝合、定型、检验和包装4过程。通过分析生产过程,得出:生产中档拉杆箱需要用7/10小时剪裁、5/10小时缝合、1
2、小时定型、1/10小时检验包装;生产高档拉杆箱则需用1小时剪裁、5/6小时缝合、2/3小时定型、1/4小时检验包装。由于公司生产能力有限,3月内各部的最大生产时间为剪裁部630小时、缝合部600小时、定型部708小时、检验包装部135小时。通过市场调研部和会计部的调查核算得出结论:生产中档拉杆箱的利润是10元,高档拉杆箱的利润是9元。公司应各生产多少中档和高档拉杆箱才能使公司利润最大?第4页,共77页,编辑于2022年,星期一实用举例实用举例 某公司通过市场调研,决定生产高中档新型拉杆箱。某分销商决定买进该公司3个月内的全部产品。拉杆箱生产需经过原材料剪裁、缝合、定型、检验和包装4过程。通过分
3、析生产过程,得出:生产中档拉杆箱需要用7/10小时剪裁、5/10小时缝合、1小时定型、1/10小时检验包装;生产高档拉杆箱则需用1小时剪裁、5/6小时缝合、2/3小时定型、1/4小时检验包装。由于公司生产能力有限,3月内各部的最大生产时间为剪裁部630小时、缝合部600小时、定型部708小时、检验包装部135小时。通过市场调研部和会计部的调查核算得出结论:生产中档拉杆箱的利润是10元,高档拉杆箱的利润是9元。公司应各生产多少中档和高档拉杆箱才能使公司利润最大?可以通过线性规划求解可以通过线性规划求解!第5页,共77页,编辑于2022年,星期一线性规划模型的建立线性规划模型的建立 设生产中、高档
4、拉杆箱数量分别为:称为决策变量。目标函数目标函数约束条件约束条件第6页,共77页,编辑于2022年,星期一一般线性规划模型Maxz=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=其中c=(c1,c2,cn),称为价值系数向量;第7页,共77页,编辑于2022年,星期一称为资源限制向量 X=(x1,x2,xn)T 称为决策变量向量
5、称为技术系数矩阵(也称消耗系数矩阵)一般线性规划模型第8页,共77页,编辑于2022年,星期一线性规划模型的标准形式线性规划模型的标准形式MaxZ=c1x1+c2x2+cnxnSubjectto(s.t.)a11x1+a12x2+a1nxn=b1a21x1+a22x2+a2nxn=b2am1x1+am2x2+amnxn=bmx x1 1 0,0,x x2 2 0,0,x xn n 0p为了讨论方便,把最大化、等式约束型的线性规划称为线性规划的标为了讨论方便,把最大化、等式约束型的线性规划称为线性规划的标准型,即准型,即第9页,共77页,编辑于2022年,星期一标准化标准化原问题原问题标准化方法
6、标准化方法目标函数目标函数Max f(x)Max f(x)Min f(x)Max f(x)约束条件约束条件引入松弛变量和人工变量引入松弛变量,不变变量变量 不变对某个对某个是任意的 第10页,共77页,编辑于2022年,星期一线性规划模型的标准化线性规划模型的标准化p例例:将如下线性规划模型标准化:将如下线性规划模型标准化:min 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-
7、3 x2+x3+x4 20 x1+2 x2+3 x3-x4 30 2x2+2 x3+3 x4-50 x1,x3,x4 0,x2取值无约束。取值无约束。目标优化标准化目标优化标准化第11页,共77页,编辑于2022年,星期一线性规划模型的标准化Maxz2=-x1-5y2+5y3+8x3+x4s.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
8、+x4 20 x1+2 x2+3 x3-x4 30 2x2+2 x3+3 x4-50 x1,x3,x4 0,x2取值无约束取值无约束 第12页,共77页,编辑于2022年,星期一线性规划模型的标准化Max z2=-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+
9、x4 20 -x1-2 y2+2y3-3 x3+x4 -30 -2y2+2y3-2 x3-3 x4 50 x1,x3,x4,y2,y3 0 第13页,共77页,编辑于2022年,星期一p可行解可行解:满足所有约束条件的解x=(x1,x2,.,xn),所有可行解的集合称为可行域可行域。p基基:设A是约束方程组的mn阶系数矩阵,秩为m,是A中阶非奇异子矩阵(即 ),则称是线性规划问题的一个基矩阵基矩阵,简称基基。B中的列向量称为基向量基向量,与基向量对应的变量x称为基变量基变量,其它变量称为非基变量非基变量。p基本解基本解:令非基变量值为0,由Ax=b可求出一个解,这个解x称为基本解基本解。p基本
10、可行解基本可行解:满足非负条件的基本解称为基本可行解基本可行解。p最优解最优解:使目标函数达到最优值的可行解称为最优解最优解。线性规划的解线性规划的解 第14页,共77页,编辑于2022年,星期一可行解、基本解、基本可行解的关系可行解可行解基本解基本解基本可行解基本可行解第15页,共77页,编辑于2022年,星期一目目 录录p线性规划实例与模型线性规划实例与模型p线性规划的图解法与基本性质线性规划的图解法与基本性质第16页,共77页,编辑于2022年,星期一线性规划的图解法线性规划的图解法 当只有两个决策变量时,当只有两个决策变量时,可用图解法求解!可用图解法求解!例例:用图解法求解以下线形规
11、划问题。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第17页,共77页,编辑于2022年,星期一解的特殊情况解的特殊情况多个最优解多个最优解第18页,共77页,编辑于2022年,星期一解的特殊情况解的特殊情况无可行解无可行解第19页,共77页,编辑于2022年,星期一解的特殊情况解的特殊情况无界解无界解第20页,共77页,编辑于2022年,星期一线性规划的基本性质线性规划的基本性质 X可行域内部的点 可行解?是是 最优解?不不 若线性规划有最若线性规
12、划有最优解,则最优解必在可行域的优解,则最优解必在可行域的顶点上达到。顶点上达到。第21页,共77页,编辑于2022年,星期一凸集的概念凸集的概念 p凸集是线性规划中一个很重要的概念,凸集的一般定义为:凸集是线性规划中一个很重要的概念,凸集的一般定义为:如果在集合如果在集合C C中任意取两个点中任意取两个点x x1 1,x x2 2,其连线上的所有点也,其连线上的所有点也都在集合都在集合C C中,则称集合中,则称集合C C为凸集合。用数学解析式表达为:为凸集合。用数学解析式表达为:对任意对任意 ,均有,均有 ,则称,则称C C是是凸集合。凸集合。非凸集非凸集非凸集非凸集凸集凸集第22页,共77
13、页,编辑于2022年,星期一线性规划的基本性质 定理定理2.12.1:线性规划的可行域:线性规划的可行域:是凸集(凸多面体)。是凸集(凸多面体)。引理引理2.12.1:线性规划的可行解线性规划的可行解 为基本可行解的为基本可行解的充分必要条件是充分必要条件是x x的正分量所对应的系数列向量是线性无关的,的正分量所对应的系数列向量是线性无关的,即每个正分量都是一个基变量即每个正分量都是一个基变量。定理定理2.22.2:线性规划问题的基本可行解线性规划问题的基本可行解x x对应于可行域的顶点对应于可行域的顶点 定理定理2.32.3:若线性规划有最优解,则最优解必在可行域的顶点上达到。若线性规划有最
14、优解,则最优解必在可行域的顶点上达到。第23页,共77页,编辑于2022年,星期一目目 录录p线性规划实例与模型线性规划实例与模型p线性规划的图解法线性规划的图解法p单纯形法原理单纯形法原理第24页,共77页,编辑于2022年,星期一一、初始基本可行解的确定一、初始基本可行解的确定p考虑如下形式的线性规划问题考虑如下形式的线性规划问题max z=c1x1+c2x2+.+cnxn s.t.a11x1+a12x2+a1nxnb1 a21x1+a22x2+a2nxnb2 (2.17).am1x1+am2x2+amnxnbm x1,x2,.,xn0 (2.18)对每个不等式约束引入一个非负松弛变量,得
15、标准形式:对每个不等式约束引入一个非负松弛变量,得标准形式: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+m0第25页,共77页,编辑于2022年,星期一是线性规划是线性规划(2.19)(2.19)的一个可行基。令的一个可行基。令 得得由此得到一个初始基本可行解为由此得到一个初始基本可行解为第26页,共77页,编辑于2022年,星期一二、最优性检验二、最优性检验 p定理定理2.42.4:在最大化问题中,对于
16、某个基本可行解,如果所有 的 ,则这个基本可行解是最优解。在最小化问题中,对于某个基本可行解,如果所有 的 ,则这个基本可行解是最优解。第27页,共77页,编辑于2022年,星期一单纯形法求解步骤单纯形法求解步骤p将所给问题化为标准形将所给问题化为标准形p找出一个初始可行基找出一个初始可行基,建立初始单纯形表建立初始单纯形表p检查所有检验数检查所有检验数(若全为非负若全为非负,则已得到最优解则已得到最优解,计算计算停止停止.否则继续下一步否则继续下一步)p考察是否无解考察是否无解(若是若是,计算停止计算停止,否则继续下一步否则继续下一步)p确定入基变量确定入基变量,出基变量出基变量p对初始单纯
17、形表进行单纯形变换对初始单纯形表进行单纯形变换第28页,共77页,编辑于2022年,星期一单纯形方法的一般步骤单纯形方法的一般步骤确定一个初始可行角点确定一个初始可行角点根据某种法则进行角点的根据某种法则进行角点的最优性判断最优性判断不是最优角点不是最优角点是最优角点是最优角点考察与当前角点相邻的考察与当前角点相邻的 “更好更好”的一个的一个角点,并置为当前角点。角点,并置为当前角点。根据某种法则进行根据某种法则进行LPLP的无界或不的无界或不可行判断可行判断无界或不可行无界或不可行还不能做出判断还不能做出判断求求解解结结束束第29页,共77页,编辑于2022年,星期一单纯形法举例单纯形法举例
18、解:引进松驰变量x3,x4,化为标准形得:第30页,共77页,编辑于2022年,星期一 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 0CBXBb 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/3第31页,共77页,编辑于2022年,星期一 4 3 0 0CBXBb x1 x2 x3 x434x2x146 0 1 3/5 -2/5 1 0
19、 -2/5 3/5 36 0 0 -1/5 -6/5表中最后一行的所有检验数均为非正数,表明目标函数已达到最大值,上述表为最优表。从表中可得到最优解为X=(x1,x2,x3,x4)=(6,4,0,0),最优值为Z=36 4 3 0 0CBXBb x1 x2 x3 x4 b/y04x3x120/326/3 0 5/3 1 -2/3 4 1 2/3 0 1/3 13 -104/3 0 1/3 0 -4/3第32页,共77页,编辑于2022年,星期一单纯形法的进一步讨论人工变量法单纯形法的进一步讨论人工变量法p人工变量法:p前面讨论了在标准型中系数矩阵有单位矩阵,很容易确定一组基可行解。在实际问题中
20、有些模型并不含有单位矩阵,为了得到一组基向量和初基可行解,在约束条件的等式左端加一组虚拟变量,得到一组基变量。这种人为加的变量称为人工变量,构成的可行基称为人工基,用大M法或两阶段法求解,这种用人工变量作桥梁的求解方法称为人工变量法。第33页,共77页,编辑于2022年,星期一单纯形法的进一步讨论人工变量法单纯形法的进一步讨论人工变量法p例1.10 用大M法解下列线性规划解:首先将数学模型化为标准形式解:首先将数学模型化为标准形式系数矩阵中不存在系数矩阵中不存在单位矩阵,无法建单位矩阵,无法建立初始单纯形表。立初始单纯形表。第34页,共77页,编辑于2022年,星期一单纯形法的进一步讨论人工变
21、量法单纯形法的进一步讨论人工变量法p故人为添加两个单位向量,得到人工变量单纯形法数学模故人为添加两个单位向量,得到人工变量单纯形法数学模型:型:其其中中:M是是一一个个很很大大的的抽抽象象的的数数,不不需需要要给给出出具具体体的的数数值值,可可以以理理解解为为它它能能大大于于给给定定的的任任何何一一个个确确定定数数值值;再再用用前前面面介介绍绍的的单单纯纯形形法法求解该模型,计算结果见下表。求解该模型,计算结果见下表。第35页,共77页,编辑于2022年,星期一单纯形法的进一步讨论人工变量法单纯形法的进一步讨论人工变量法cj32-100-M-MCBXBbx1x2x3x4x5x6x7i-Mx64
22、-431-101040 x5101-1201005-Mx712-21000113-2M2+M-1+2M-M0000 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/3第36页,共77页,编辑于2022年,星期一单纯形法的进一步讨论人工变量法单纯形法的进一步讨论人工变量法解的判别:解的判别:1)唯一最优解判别:最优表中所有非基变量
23、的检验数非零)唯一最优解判别:最优表中所有非基变量的检验数非零,则则线线 规划具有唯一最优解。(检验数为负数或规划具有唯一最优解。(检验数为负数或0)2)多重最优解判别:最优表中存在非基变量的检验数为零)多重最优解判别:最优表中存在非基变量的检验数为零,则线则则线则性规划具有无穷多最优解(检验数为负数或性规划具有无穷多最优解(检验数为负数或0)。)。3)无界解判别:某个检验数)无界解判别:某个检验数0且且aik(i=1,2,m)则线性)则线性规划具有无界解。规划具有无界解。4)无无可可行行解解的的判判断断:当当用用大大M单单纯纯形形法法计计算算得得到到最最优优解解并并且且存存在在Ri0时,则表
24、明原线性规划无可行解。时,则表明原线性规划无可行解。5)退化解的判别:存在某个基变量为零的基本可行解。)退化解的判别:存在某个基变量为零的基本可行解。第37页,共77页,编辑于2022年,星期一大大M M法法基本思想:基本思想:在约束条件中增加人工变量,同时修改目标函数,对目标函数加上具有很大正系数M的惩罚项,为使人工变量不影响目标函数取值,在迭代过程中,应把人工变量从基变量中退出,使其成为非基变量。其中,M是很大的正数,同时增加两个人工变量。不容易找到初始可行解第38页,共77页,编辑于2022年,星期一找初始可行解找初始可行解11-300MMbi/yibx1X2x3x4x5x6x70X41
25、11-21100011MX6321-40-1103/2MX71(1)0-2000111-3M1-M-3+6M0M00 要求最后得到的基变量中不含人工变量X1进基,x7出基11-300MMbi/yibx1x2x3x4x5x6x7-3X340011/3-2/32/3-5/31X210100-11-211X191002/3-4/34/3-7/30001/31/3M-1/3 M-2/3 可以从表中移去,然后继续求解 经过几次变换,得到基变量为x1,x2,x3第39页,共77页,编辑于2022年,星期一两阶段法原理两阶段法原理p两阶段法的第一阶段是用单纯形法 消去人工变量,即把人工变量都变为非基变量,求
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 线性规划 单纯 幻灯片
限制150内