运筹学OR1-Ch1-LP原理课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《运筹学OR1-Ch1-LP原理课件.ppt》由用户(三亚风情)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 运筹学 OR1 Ch1 LP 原理 课件
- 资源描述:
-
1、第一章第一章 线性规划原理线性规划原理Chapter1 the Principle ofLiner Programming1-1 LP问题的提出问题的提出1-2 LP模型及其解的概念模型及其解的概念1-3 LP图解法图解法1-4 LP的求解原理的求解原理1-5 LP在工商管理中的应用在工商管理中的应用1-1 LP问题的提出问题的提出例1:某工厂在计划期内要安排生产甲、乙两种产品,这两种产品分别需要在A、B、C三种不同的设备上加工。按工艺规定,产品甲、乙在各设备上所需的加工台时如表1-1。已知各设备在计划期内有效台时分别是8、16和12,该工厂每生产一件产品甲、乙可分别获得利润2元、3元。问应如
2、何安排生产计划,才能得到利润最多?设备 产品 A B C 单件利润 甲 1 4 0 2 乙 2 0 4 3 有效台时 8 16 12 2022年7月28日7时46分Page 2 of 43 LP问题的提出问题的提出此问题用数学语言描述为:v假设x1,x2分别表示在计划期内产品甲、乙的产量,则该企业的目标为利润z最大:z=2x1+3x2 v这个方程是实际问题追求的目标,因此称之为问题的目标方程。在此方程中,利润Z是x1,x2的单增函数,x1,x2越大,利润Z就越大。但是,由于受到有限资源的限制(此题中为设备台时的限制),x1,x2不可能任意增大,它要受到资源量的制约。比如,对于设备A而言,其总有
3、效台时为8,生产产品甲乙所需要的A设备的总台时不能超过8,即:x1+2x28v同样地,对于设备资源B、C,我们可以得到不等式:4x1164x212 2022年7月28日7时46分Page 3 of 43 LP问题的提出问题的提出v除了这些设备有效台时的限制条件之外,还有一个隐含的限制条件,就是产量x1,x2都必须是大于等于零的。v综上所述,得到问题的数学模型如下:Max z=2x1+3x2 x1+2x28 4x1 16 4x212 x1,x20例2:靠近某河流有两个化工厂(如图1-1),流经第一化工厂的河流流量为每天500万m3,在两个工厂之间有一条流量为每天200万m3的支流,第一化工厂每天
4、排放含有某种有害物质的工业污水2万m3,第二化工厂每天排放这种污水1.4万m3。从第一化工厂排出的工业污水流到第二化工厂以前有20%可自然2022年7月28日7时46分Page 4 of 43 LP问题的提出问题的提出净化。根据环保要求,河流中工业污水的含量应不大于0.2%。为此这两个工厂都需各自处理一部分工业污水以满足环保要求。第一化工厂处理工业污水的成本是1000元/万m3,第二化工厂是800元/万m3。现在要问在满足环保要求的条件下,每厂各应处理多少工业污水,使这两个工厂总的处理工业污水的费用最少。500 万3m 200 万3m O 工厂 1 O 工厂 2 例 2 示意图 设第一和第二化
5、工厂每天处理工业污水量分别为x1、x2万,则问题追求的目标为:费用 Min w=1000 x1+800 x2 且要满足环保条件:污水含量0.2%。2022年7月28日7时46分Page 5 of 43 LP问题的提出问题的提出为此,将河流分成两段考虑,第一段:工厂1工厂2之间,只有工厂1的排放污水,应满足的约束条件:(2x1)/5002/1000第二段:工厂2之后的河段,除了有工厂2的排放污水外,还有工厂1的残留污水,同时,水流量也是两个河流的流量和,即:0.8(2x1)+(1.4x2)/7002/1000除此之外,还应满足处理量不超过最大污水量和非负条件。综合以上分析得问题的数学模型为:21
6、8001000minxxz10002500)2(1 x10002700)4.1()2(8.021xx21x4.12x0,21xx2022年7月28日7时46分Page 6 of 43 1-2 LP模型及其解的概念模型及其解的概念1.线性规划模型的三要素:1)定义决策变量集合:x1,x2,xn,这组决策变量的每一组可能的值就代表问题的一个具体解决方案;2)定义目标函数:用决策变量的线性函数形势写出所要追求的目标,称之为目标函数:z=f(x1,x2,xn),根据问题的不同,要求目标函数实现最大化或最小化;3)定义约束条件集合:用一组决策变量的等式或不等式来表示在决策问题过程中所必须遵循的约束条件。
7、2.线性规划的一般模型 目标函数:约束条件:nnxcxcxcz2211max(min)11212111),(bxaxaxann22222121),(bxaxaxannmnmnmmbxaxaxa),(22110,21nxxx2022年7月28日7时46分Page 7 of 43 LP模型及其解的概念模型及其解的概念3.线性规划模型的矩阵描述形式0XbAXCXzmaxmn2m1mn22221n11211Tm21n21Tn21aaaaaaaaaAbbbbcccCxxxX约束系数矩阵资源列向量目标价值系数其中:决策变量集合),(),(),(n21j0 xm21ibxaxczjn1jijijn1jjj,
8、max还可以表示为:2022年7月28日7时46分Page 8 of 43 LP模型及其解的概念模型及其解的概念4.LP模型的标准型v概括线性规划在工农业生产方面的应用有两大类:第一类为在资源限制条件下使利益最大,第二类为在产出一定的前提下使费用最小。v从实际问题产生的线性规划模型,其目标函数可能要求最大或最小,约束条件可能是、或兼而有之,为探讨一种通用的解法,有必要将其标准化。vLP模型的标准型定义为:nnxcxcxcz2211max11212111bxaxaxann22222121bxaxaxannmnmnmmbxaxaxa22110,21nxxx2022年7月28日7时46分Page 9
9、 of 43 LP模型及其解的概念模型及其解的概念5。化一般LP模型为标准型v若问题的目标函数为最小化 ,则令 ,问题变成 ,得到最终结果 后反求 。v若约束条件为不等式,分两种情况,一种是约束,则将该不等式的左端加上一个非负的变量,使其变为等式,该变量称为松弛变量;另一种是约束,则将该不等式的左端减去一个非负的变量,使其变为等式,该变量称为剩余变量(也可统称为松弛变量)。v若某一决策变量 无非负约束,则用 代替,其中 。CXZ minCXZZCXZmaxZZZkxkkxx0,kkxx2022年7月28日7时46分Page 10 of 43 LP模型及其解的概念模型及其解的概念v例3:试将下面
10、线性规划问题化成标准型。32132minxxxz7321xxx2321xxx523321xxx321,0,xxx无符号要求 解:通过以下步骤:1.用 替代 ,其中 ;2.在第一个约束条件的号左端加入松弛变量 ;3.在第二个约束条件的号左端减去剩余变量 ;4.令 ,变求 为求 。)(33xx 3x0,33xx4x5xZZZminZmax54 332100)(32maxxxxxxxZ7)(4 3321xxxxx2)(5 3321xxxxx5)(23 3321xxxx71,0jxj2022年7月28日7时46分Page 11 of 43 LP模型及其解的概念模型及其解的概念6 LP解的概念v针对LP
11、的标准型.bAX .0X 可行解:满足所有约束条件包括非负条件的解 X=(x1,x2,xn)T 称为该问题的可行解。所有可行解的集合构成可行域。最优解:使目标函数达到最大值的可行解称为最优解。基:A是约束方程组的 mn 阶系数矩阵,设其秩为m,B矩阵是矩阵A中mm阶非奇异子矩阵(|B|0),则称B是线性规划问题的一个基。就是说,B是由m个线性独立的列向量组成。不失一般性,可以设:),(21212222111211mmmmmmmpppaaaaaaaaaBCXz max2022年7月28日7时46分Page 12 of 43 LP模型及其解的概念模型及其解的概念v称pj(j=1,2,m)为基向量,
12、与基向量pj对应的变量xj(j=1,2,m)为基变量,否则为非基变量。基本解:对于上述基B,其秩为m,一般mn有无数多个解,则上述式可变为:设XB 是对应于这个基B的基变量向量:若令 ,则可用高斯消元法,求出上式的一个解:,这个解非零分量的数目不大于方程数m,则称X为基本解。nnmmmmxpxpbxpxpxp112211TmBxxxX),(21021nmmxxxTmxxxX)0,0,(21基本可行解:满足非负条件的基本解称之为基本可行解。可行基:对应于基本可行解的基称为可行基。2022年7月28日7时46分Page 13 of 43 LP模型及其解的概念模型及其解的概念由上面分析可知,约束方程
13、组基本解的数目是有限的,最多为 个,一般情况下基本可行解的数目要小于基本解的数目,最多相等。另外,在基本解中,若非零分量的数目小于m,则称之为退化解。我们暂不讨论这种情况。非可行解 基本解 可行解 基本可行解 图 1-3:解的关系示意图 mnC2022年7月28日7时46分Page 14 of 43 1-3 LP的图解法现对上述例1进行图解,在以x1 x2 为坐标轴的直角坐标系中,原题中的每一个约束条件为一条直线(如图所示)。再来分析目标函数,把它反映在直角坐标系下,可以表示为以Z为参数的一族平行线:33212zxxx1x2152431234R*),(24Q16x41 12x42 8x2x21
14、 3zx32x12 图中粗线所围区域部分中(含边界)中每一点都满足约束条件,都是该问题的解,这个区域就是该问题的解集合,又称之为可行域。2022年7月28日7时46分Page 15 of 43 线性规划的图解法v注意:几种特殊情况:1.当目标方程直线与某一约束直线平行时,最优值不唯一。如,本例中若将目标函数改为max z=2x1+4x2。则其直线与约束x1+2x28的直线平行,则直线x1+2x28在可行域中的部分,其每一点都是问题的最优解。21maxxxz0 xx2xx4xx2212121 ,x1x2152431234R16x41 12x42 8x2x21 21x4x2z 2xx21 4xx2
15、21 2.有可行域,但无最优解。即目标函数的值z+如:2022年7月28日7时46分Page 16 of 43 线性规划的图解法3.无可行解。当约束条件出现相互矛盾时,则没有可行域。在实际问题中,当数学模型有错误时,才可能发生2、3种情况 通过以上图解,较直观地看到了线性规划模型的求解过程。对于复杂的问题是不可能用图解法求解的,它只能求解两个决策变量的情况。因此,还需要找到一种线性规划问题的一般的、通用的求解方法。线性规划的实用解法有:单纯形法;计算机软件求解。2022年7月28日7时46分Page 17 of 43 线性规划的图解法vLP图解法得到的启示1.若LP问题有最优解,则该解一定在可
16、行域的边界上(由目标求max或min);2.若LP问题有最优解,则一定在其可行域的某个顶点上达到最优;3.最优解一定存在于基本可行解中(理论定理证明P.1619)1)凸集定义:n维空间一个点集中的任意两点连线上的点都属于这个点集。2)顶点定义:凸集中不能用不同两点线性表示的点称为。3)定理1:LP问题的可行域是凸集。4)定理2:LP问题的基本可行解对应于其可行域的顶点。5)定理3:可行域有界的LP问题,其目标函数一定可以在其可行域的顶点上达到最优。2022年7月28日7时46分Page 18 of 43 1-4 LP的求解原理的求解原理1.解法思想:采用迭代的办法。首先找一个初始基本可行解,然
17、后判断是否是最优解,若是,则停止迭代;若不是,则以此解为基础,找比它更好的基本可行解。如此循环,直至找到最优解。其框图描述如图1-4。确定初始基本可行解 是最优解?寻找更好的基本可行解 停 是 否 图 1-4:LP 解法思想示意图 2022年7月28日7时46分Page 19 of 43 LP的求解原理的求解原理2.初始基本可行解的确定 1)若线性规划问题的约束方程均为“=”约束,即且从其系数列向量 中能观察到存在m个线性独立的单位向量,则由这些列向量构成初始基:,其对应的决策变量为初始基变量,其余变量为非基变量,且令非基变量等于零,则构成初始基本可行解。mibxpijnjj1,1)(n1jp
18、j100010001),(210mpppB2022年7月28日7时46分Page 20 of 43 LP的求解原理的求解原理2.初始基本可行解的确定 2)若原问题的约束条件均为“”约束,则给每一个约束条件的左端分别加上一个非负的松弛变量变成为等式约束。这m个松弛变量所对应的系数列向量必构成单位阵I,以此作为初始基:即以m个松弛变量为初始基变量,其余变量为非基变量并令其等于零,则得原问题的初始基本可行解。100010001),(210mnnnpppB2022年7月28日7时46分Page 21 of 43 LP的求解原理的求解原理2.初始基本可行解的确定 3)若原问题均为“=”约束,且观察不到存
19、在相互独立的m个单位列向量,则给每个约束左端加上一个非负的人工变量构成初始基,以此确定初始基本可行解。4)若原问题的约束条件均为“”约束,则按标准化的步骤给每个约束减去一个非负剩余变量,使其变为“=”;然后,再给每个约束的左端加上一个非负的人工变量,这m个人工变量所对应的系数列向量必构成单位阵,以此单位阵为初始基所形成的解为初始基本可行解。5)若约束条件中混合出现“”、“”、“=”,则综合运用上述原则,总可以得到m个相互独立的单位列向量构成为m阶单位阵,以此为初始基,求得初始基本可行解。注意:人工变量是为解决问题人为硬加入的不存在的量。因此,最终的解中这些人工变量必须成为非基变量,否则问题无解
20、。2022年7月28日7时46分Page 22 of 43 LP的求解原理的求解原理2.初始基本可行解的确定 例1初始解的确定51j0 x12xx416xx48xx2xx0 x0 x0 x3x2zj524132154321,max0 xx12x416x48x2xx3x2z21212121,max标准化100400100400121pppA521),(系数矩阵为:),(5430pppB基矩阵:)基向量:(543xxx,2514213x412xx416xx2x8x则有:0 x3x20Z0 xx2121得:令T01216800X),()(得到初始解:2022年7月28日7时46分Page 23 of
21、 43 LP的求解原理的求解原理3.最优性检验 按照线性规划的求解思路,在得到其初始可行解之后,就要检查一下,它是否是最优解,若是,则算法停止;若不是,继续迭代求解下一个基本可行解,并且每求得一个基本可行解都要检验它是否是最优解。那么,如何来检验呢?一般而言,经过迭代可得原问题的一个基,对于所有的约束,基变量留在等式的左端,所有非基变量移到右端,得到:m21ixabxjn1mjijii,jijm1iin1mjjim1iixaccbcz)(带入目标方程得:im1ii0bcz令:ijm1iijacz令:jjn1mjj0 xzczz)(得到:2022年7月28日7时46分Page 24 of 43
22、LP的求解原理的求解原理3.最优性检验 n1mjzcjjj,再令:jn1mjj0 xzz则得:从此式可以作出判断:若当前基所形成的j,对于所有的 j()都有j0,根据变量的非负性 (),有jxj0,所以只有当 所有 时,目标函数才能取得最大值。因此,称 为检验数。nmj1njxj1,0)(n1mj0 xjj2022年7月28日7时46分Page 25 of 43 LP的求解原理的求解原理3.最优性检验 1)最优性判定:若对应于基B的基本可行解为:且对于一切非基变量 ,有j0,则 为最优解。TmlbbbX)0,0,(21)(jxnmj1)(lX2)无穷多最优解判定:若对应于基B的基本可行解 ,均
展开阅读全文