1、运筹学课件,2,运筹学 课件,版权所有，未经准许，不得翻制,3,本课件配合运筹学教材，是总结我们多年教学中积累的教学课件做成的光盘。在这里为了支撑教师的教学和学生的学习，奉献给读者特别是教师。 其中的内容不是教材的简单复制，如此是为了扩大整个教学的信息量，作为附件仅供教师与其他读者参考。由于我们的水平，以及对工作投入的限制，课件存在许多不足和问题，诚挚欢迎提出宝贵的意见和建议。,运 筹 学,1.绪论2.线性规划建模及单纯形法3.线性规划问题的对偶与灵敏度分析4.运输问题5.动态规划6.排队论7.决策分析8.图与网络分析,5,第一章 绪 论,6,运筹学概况简述,运筹学(Operations Re
2、search)直译为“运作研究”。运筹学是运用科学的方法(如分析、试验、量化等)来决定如何最佳地运营和设计各种系统的一门学科。,7,运筹学概况简述,运筹学能够对经济管理系统中的人力、物力、财力等资源进行统筹安排，为决策者提供有依据的最优方案，以实现最有效的管理。通常以最优、最佳等作为决策目标，避开最劣的方案。,8,运筹学在工商管理中的应用,生产计划：生产作业的计划、日程表的编排、合理下料、配料问题、物料管理等。库存管理：多种物资库存量的管理，库存方式、库存量等。运输问题：确定最小成本的运输线路、物资的调拨、运输工具的调度以及建厂地址的选择等。,9,运筹学在工商管理中的应用,人事管理：对人员的需
3、求和使用的预测，确定人员编制、人员合理分配，建立人才评价体系等。市场营销：广告预算、媒介选择、定价、产品开发与销售计划制定等。,10,运筹学在工商管理中的应用,财务和会计：包括预测、贷款、成本分析、定价、证券管理、现金管理等。其他： 设备维修、更新，项目选择、评价，工程优化设计与管理等。,11,运筹学的产生和发展,运筹学思想的出现可以追溯到很早“田忌齐王赛马”（对策论）、孙子兵法等都体现了优化的思想。“运筹学”这一名词最早出现在第二次世界大战期间 美、英等国家的作战研究小组为了解决作战中所遇到的许多错综复杂的战略、战术问题而提出的。,12,运筹学的产生和发展,战后这些研究成果被应用到生产、经济
4、领域，并得到迅速发展有关理论和方法的研究、实践不断深入。1947年美国数学家丹捷格(G.B.Dantzig）提出了求解线性规划的有效方法单纯形法。,13,运筹学的产生和发展,数学对运筹学的作用是有关理论和方法的研究基础，是建立运筹学模型的工具。计算机的发展，促进运筹学的进一步发展高速、可靠的计算是运筹学解决问题的基本保障。,14,运筹学的分支,线性规划 非线性规划 整数规划 动态规划,多目标规划 随机规划 模糊规划等,15,运筹学的分支,图与网络理论 存储论 排队论 决策论,对策论 排序与统筹方法 可靠性理论等,16,运筹学方法使用情况(美1983),17,运筹学方法在中国使用情况 (随机抽样
5、),18,运筹学的推广应用前景,据美劳工局1992年统计预测:社会对运筹学应用分析人员的需求从1990年到2005年，其增长百分比预测为73%,增长速度排到各项职业的前三位。,19,运筹学的推广应用前景,结论: -运筹学在国内或国外的推广应用前景是非常广阔的。 -工商企业对运筹学应用的需求是很大的。 -在工商企业推广运筹学方面有大量的工作要做。,20,运筹学解决问题的过程,1）提出问题：认清问题。2）寻求可行方案：建模、求解。3）确定评估目标及方案的标准或方法、途径。4）评估各个方案：解的检验、灵敏性分析等。,21,运筹学解决问题的过程,5）选择最优方案：决策。6）方案实施：回到实践中。7）后
6、评估：考察问题是否得到完满解决。1）2）3）形成问题；4）5）分析问题：定性分析与定量分析相结合，构成决策。,22,如何学习运筹学课程,学习运筹学要把重点放在分析、理解有关的概念、思路上。在自学过程中，应该多向自己提问，例如一个方法的实质是什么，为什么这样进行，怎么进行等。自学时要掌握三个重要环节：,23,如何学习运筹学课程,1.认真阅读教材和参考资料，以指定教材为主，同时参考其他有关书籍。一般每一本运筹学教材都有自己的特点，但是基本原理、概念都是一致的。注意主从，参考资料会帮助你开阔思路，使学习深入。但是，把时间过多放在参考资料上，会导致思路分散，不利于学好。,24,2.要在理解了基本概念和
7、理论的基础上研究例题，注意例题是为了帮助理解概念、理论的。作业练习的主要作用也是这样，它同时还有让你自己检查自己学习的作用。因此，做题要有信心，要独立完成，不要怕出错。因为，整个课程是一个整体，各节内容有内在联系，只要学到一定程度，知识融会贯通起来，你自己就能够对所做题目的正确性作出判断。,如何学习运筹学课程,25,3、要学会做学习小结。每一节或一章学完后，必须学会用精炼的语言来概述该书所讲内容。这样，你才能够从较高的角度来看问题，更深刻地理解有关知识和内容。这就称作“把书读薄”，若能够结合相关参考文献并深入理解，把相关知识从更深入、广泛的角度进行论述，则称为“把书读厚”。,如何学习运筹学课程
8、,26,在建数学模型时，要结合实际应用。,如何学习运筹学课程,27,运筹学课件,28,第二章 线性规划建模及单纯形法,本章内容重点,线性规划模型与解的主要概念 线性规划的单纯形法，线性规划多解分析 线性规划应用建模,29,1.线性规划的概念,例2.1:某工厂拥有A、B、C三种类型的设备，生产甲、乙两种产品。每件产品在生产中需要占用的设备机时数，每件产品可以获得的利润以及三种设备可利用的时数如下表所示：,30,问题：工厂应如何安排生产可获得最大的总利润？解：设变量xi为第i种（甲、乙）产品的生产件数（i1，2）。根据题意，我们知道两种产品的生产受到设备能力（机时数）的限制。对设备A，两种产品生产
9、所占用的机时数不能超过65，于是我们可以得到不等式：3 x1 + 2 x2 65；对设备B，两种产品生产所占用的机时数不能超过40，于是我们可以得到不等式：2 x1 + x2 40；,1.线性规划的概念,31,对设备C，两种产品生产所占用的机时数不能超过75，于是我们可以得到不等式：3x2 75 ；另外，产品数不可能为负，即 x1 ,x2 0。同时，我们有一个追求目标，即获取最大利润。于是可写出目标函数z为相应的生产计划可以获得的总利润：z=1500x1+2500x2 。综合上述讨论，在加工时间以及利润与产品产量成线性关系的假设下，把目标函数和约束条件放在一起，可以建立如下的线性规划模型：,1
10、.线性规划的概念,32,目标函数 Max z =1500x1+2500x2 约束条件 s.t. 3x1+2x2 652x1+x2 403x2 75x1 ,x2 0,1.线性规划的概念,33,这是一个典型的利润最大化的生产计划问题。其中，“Max”是英文单词“Maximize”的缩写，含义为“最大化”；“s.t.”是“subject to”的缩写，表示“满足于”。因此，上述模型的含义是：在给定条件限制下，求使目标函数z达到最大的x1 ,x2 的取值。,1.线性规划的概念,34,一般形式 目标函数： Max(Min)z = c1x1 + c2x2 + + cnxn,约束条件： a11x1+a12x
11、2+a1nxn( =, )b1 a21x1+a22x2+a2nxn( =, )b2 . am1x1+am2x2 +amnxn( =, )bmx1 ，x2 ， ，xn 0,1.线性规划的概念,35,标准形式 目标函数： Max z = c1x1 + c2x2 + + cnxn,约束条件： a11x1 + a12x2 + + a1nxn = b1 a21x1 + a22x2 + + a2nxn = b2 . am1x1 + am2x2 + + amnxn = bmx1 ，x2 ， ，xn 0,1.线性规划的概念,36,可以看出，线性规划的标准形式有如下四个特点：目标最大化、约束为等式、决策变量均非
12、负、右端项非负。对于各种非标准形式的线性规划问题，我们总可以通过以下变换，将其转化为标准形式:,1.线性规划的概念,37,1.极小化目标函数的问题：设目标函数为Min f = c1x1 + c2x2 + + cnxn 则可以令z -f ，该极小化问 题与下面的极大化问题有相同的最优 解，即Max z = -c1x1 - c2x2 - - cnxn 但必须注意，尽管以上两个问题的最优解相同，但他们最优解的目标函数值却相差一个符号，即Min f - Max z,1.线性规划的概念,38,2、约束条件不是等式的问题：设约束条件为ai1 x1+ai2 x2+ +ain xn bi可以引进一个新的变量s
13、 ，使它等于约束右边与左边之差s=bi(ai1 x1 + ai2 x2 + + ain xn )显然，s 也具有非负约束，即s0，这时新的约束条件成为ai1 x1+ai2 x2+ +ain xn+s = bi,1.线性规划的概念,39,当约束条件为ai1 x1+ai2 x2+ +ain xn bi 时，类似地令s=(ai1 x1+ai2 x2+ +ain xn)- bi 显然，s 也具有非负约束，即s0，这时新的约束条件成为ai1 x1+ai2 x2+ +ain xn-s = bi,1.线性规划的概念,40,为了使约束由不等式成为等式而引进的变量s称为“松弛变量”。如果原问题中有若干个非等式约
14、束，则将其转化为标准形式时，必须对各个约束引进不同的松弛变量。,1.线性规划的概念,41,例2.2：将以下线性规划问题转化为标准形式 Min f = 3.6 x1 - 5.2 x2 + 1.8 x3s. t. 2.3 x1 + 5.2 x2 - 6.1 x3 15.74.1 x1 + 3.3 x3 8.9x1 + x2 + x3 = 38x1 , x2 , x3 0,1.线性规划的概念,解：首先,将目标函数转换成极大化： 令 z= -f = -3.6x1+5.2x2-1.8x3,42,其次考虑约束，有2个不等式约束，引进松弛变量x4，x5 0。 于是，我们可以得到以下标准形式的线性规划问题：
15、Max z = - 3.6 x1 + 5.2 x2 - 1.8 x3 s.t. 2.3x1+5.2x2-6.1x3+x4= 15.74.1x1+3.3x3-x5= 8.9x1+x2+x3= 38x1 ,x2 ,x3 ,x4 ,x5 0,1.线性规划的概念,43,3. 变量无符号限制的问题： 在标准形式中，必须每一个变量均有非负约束。当某一个变量xj没有非负约束时，可以令xj = xj- xj” 其中xj0，xj”0 即用两个非负变量之差来表示一个无符号限制的变量，当然xj的符号取决于xj和xj”的大小。,1.线性规划的概念,44,4.右端项有负值的问题： 在标准形式中，要求右端项必须每一个分量
16、非负。当某一个右端项系数为负时，如 bi0，则把该等式约束两端同时乘以-1，得到： -ai1 x1-ai2 x2- -ain xn = -bi 。,1.线性规划的概念,45,例2.3：将以下线性规划问题转化为标准形式 Min f= -3 x1 + 5 x2 + 8 x3 - 7 x4 s.t. 2 x1 - 3 x2 + 5 x3 + 6 x4 284 x1 + 2 x2 + 3 x3 - 9 x4 396 x2 + 2 x3 + 3 x4 - 58x1 , x3 , x4 0,1.线性规划的概念,46,解：首先，将目标函数转换成极大化： 令 z = -f = 3x15x28x3+7x4 ；其
17、次考虑约束，有3个不等式约束，引进松弛变量x5 ,x6 ,x7 0 ； 由于x2无非负限制，可令x2=x2-x2”,其中 x20，x2”0 ；由于第3个约束右端项系数为-58，于是把该式两端乘以-1 。于是，我们可以得到以下标准形式的线性规划问题：,1.线性规划的概念,47,Max z = 3x15x2+5x2”8x3 +7x4s.t. 2x13x2+3x2”+5x3+6x4+x5= 284x1+2x2-2x2”+3x3-9x4-x6= 39-6x2+6x2”-2x3-3x4-x7 = 58x1 ,x2,x2”,x3 ,x4 ,x5 ,x6 ,x7 0,1.线性规划的概念,48,2.线性规划的
18、图解法,线性规划的图解法(解的几何表示)对于只有两个变量的线性规划问题，可以二维直角坐标平面上作图表示线性规划问题的有关概念，并求解。图解法求解线性规划问题的步骤如下： (1)分别取决策变量x1 ,x2 为坐标向量建立直角坐标系。,49,2.线性规划的图解法,（2）对每个约束（包括非负约束）条件，先取其等式在坐标系中作出直线，通过判断确定不等式所决定的半平面。各约束半平面交出来的区域(存在或不存在），若存在，其中的点表示的解称为此线性规划的可行解。这些符合约束限制的点集合，称为可行集或可行域。然后进行（3）。否则该线性规划问题无可行解。,50,2.线性规划的图解法,（3）任意给定目标函数一个值
19、作一条目标函数的等值线，并确定该等值线平移后值增加的方向，平移此目标函数的等值线，使其达到既与可行域有交点又不可能使值再增加的位置(有时交于无穷远处，此时称无有限最优解）。若有交点时，此目标函数等值线与可行域的交点即最优解（一个或多个），此目标函数的值即最优值。,51,2.线性规划的图解法,例2.4:某工厂拥有A、B、C三种类型的设备，生产甲、乙两种产品。每件产品在生产中需要占用的设备机时数，每件产品可以获得的利润以及三种设备可利用的时数如下表所示：,52,2.线性规划的图解法,问题：工厂应如何安排生产可获得最大的总利润？用图解法求解。解：设变量xi为第i种（甲、乙）产品的生产件数（i1，2）
20、。根据前面分析，可以建立如下的线性规划模型：Max z = 1500 x1 + 2500 x2s.t. 3x1+2x2 65 (A)2x1+x2 40 (B)3x2 75 (C)x1 ,x2 0 (D, E),53,2.线性规划的图解法,按照图解法的步骤在以决策变量x1 ，x2 为坐标向量的平面直角坐标系上对每个约束（包括非负约束）条件作出直线，并通过判断确定不等式所决定的半平面。各约束半平面交出来的区域即可行集或可行域如下图阴影所示。,54,2、线性规划的图解法,图解法求解线性规划,55,2.线性规划的图解法,任意给定目标函数一个值作一条目标函数的等值线，并确定该等值线平移后值增加的方向，平
21、移此目标函数的等值线，使其达到既与可行域有交点又不可能使值再增加的位置，得到交点 (5,25)T ，此目标函数的值为70000。于是，我们得到这个线性规划的最优解x1=5、x2=25，最优值z = 70000。即最优方案为生产甲产品5件、乙产品25件，可获得最大利润为70000元。,56,2.线性规划的图解法,例2.5：在例2.1的线性规划模型中，如果目标函数变为：Max z = 1500 x1 + 1000 x2那么，最优情况下目标函数的等值线与直线（A）重合。这时，最优解有无穷多个，是从点 (5,25)T到点 (15,10)T 线段上的所有点，最优值为32500。如下图所示：,57,2.线
22、性规划的图解法,无穷多解的情况,58,2.线性规划的图解法,例2.6：在例2.1的线性规划模型中,如果约束条件（A）、（C）变为：3 x1 + 2 x2 65 (A）3 x2 75 (C）并且去掉（D、E）的非负限制。那么，可行域成为一个上无界的区域。这时，没有有限最优解，如下图所示：,59,2.线性规划的图解法,无有限解的情况,60,2、线性规划的图解法,例2.7：在例2.1的线性规划模型中，如果增加约束条件（F）为：x1 + x2 40 （F）那么，可行域成为空的区域。这时， 没有可行解，显然线性规划问题无解。如下图所示：,61,2.线性规划的图解法,无可行解的情况,62,根据以上例题，进
23、一步分析讨论可知线性规划的可行域和最优解有以下几种可能的情况1.可行域为封闭的有界区域(a)有唯一的最优解；(b)有无穷多个最优解；2.可行域为封闭的无界区域(c)有唯一的最优解；,2.线性规划的图解法,63,(d)有无穷多个最优解；(e)目标函数无界(即虽有可行解，但在可行域中，目标函数可以无限增大或无限减少)，因而没有有限最优解。3.可行域为空集(f)没有可行解，原问题无最优解,2.线性规划的图解法,64,以上几种情况的图示如下：,2.线性规划的图解法,可行域有界唯一最优解,可行域有界多个最优解,65,2.线性规划的图解法,可行域无界唯一最优解,可行域无界无穷多最优解,66,2.线性规划的
24、图解法,可行域无界目标函数无界,可行域为空集无可行解,67,可行解、可行解集（可行域） 最优解、最优值 基、基变量、非基变量 基本解、基本可行解 可行基、最优基,熟悉下列一些解的概念,2.线性规划解的概念,68,线性规划的基、基本解与基本可行解在一般情况下，由于图解法无法解决三个变量以上的线性规划问题，对于n个变量的线性规划问题，我们必须用解方程的办法来求得可行域的极点。再来进一步考察前例。例2.8 把例2.1的线性规划模型标准化，引入松驰变量 x3 ，x4 ，x5 0，得到,2.线性规划解的概念,69,Max z = 1500 x1 + 2500 x2 s.t. 3x1+2x2+x3= 65
25、 (A)2x1+x2+x4= 40 (B)3x2+x5= 75 (C)x1 ,x2 ,x3 ,x4 ,x5 0 用（D）（E）（F）（G）（H） 分别表示x1 = 0、x2 = 0、x3 = 0、 x4 = 0、x5 = 0 。 这里一共有8个约束条件，其中3个等式约束,2.线性规划解的概念,70,（一般情况下，等式约束的个数少于决策变量的个数），5个变量非负约束（与决策变量个数相同）。每5个方程若线性无关可解得一个点，我们可以看到前例图解法得到的区域中每两条直线的交点与此例的各个方程有如下关系：见下图。,2.线性规划解的概念,71,2.线性规划解的概念,平面上各不等式约束半平面得交点,72,
26、由上图可以看出：直线A、B的交点对应于约束条件(A)、（B）、（C）、（F）、（G）的解，即：x(1) = (15,10,0,0,45)T直线A、C的交点对应于约束条件(A)、（B）、（C）、（F）、（H）的解，即：x(2) = (5,25,0,5,0)T直线A、D的交点对应于约束条件(A)、（B）、（C）、（D）、（F）的解，即：x(3) = (0,32.5,0,7.5,-22.5)T,2.线性规划解的概念,73,直线A、E的交点对应于约束条件(A)、（B）、（C）、（E）、（F）的解，即：x(4) = (65/3,0,0,-10/3,75)T直线B、C的交点对应于约束条件(A)、（B）、（
27、C）、（G）、（H）的解，即：x(5) = (7.5,25,-7.5,0,0)T直线B、D的交点对应于约束条件(A)、（B）、（C）、（D）、（G）的解，即：x(6) = (0,40,-15,0,-45)T,2.线性规划解的概念,74,直线B、E的交点对应于约束条件(A)、（B）、（C）、（E）、（G）的解，即：x(7) = (20,0,5,0,75)T直线C、D的交点对应于约束条件(A)、（B）、（C）、（D）、（H）的解，即：x(8) = (0,25,15,15,0)T直线C、E无交点（C、E相互平行）直线D、E的交点对应于约束条件(A)、（B）、（C）、（D）、（E）的解，即：x(9)
28、= (0,0,65,40,75)T,2.线性规划解的概念,75,上图各约束直线的交点是由以下方法得到：在标准化的等式约束中，令其中某两个变量为零，得到其他变量的唯一解，这个解就是相应交点的坐标，如果某一交点的坐标 (x1 , x2 , x3 , x4 , x5 )T全为非负，则该交点就对应于线性规划可行域的一个极点（如A、B，A、C，B、E，C、D和D、E的交点）；如果某一交点的坐标中至少有一个分量为负值（如A、D，A、E，B、C和B、D的交点），则该交点不是可行域的极点。,2.线性规划解的概念,76,由上图可知，A、B交点对应于 x3 = 0， x4 = 0，在等式约束中令x3 = 0，x4
29、 = 0，得到x1 =15，x2 = 10，x5 = 45。即A、B交点对应于极点x= (x1 ，x2 ，x3 ，x4 ，x5)T =(15，10，0，0，45)T。由于所有分量都为非负，因此A、B交点是可行域的极点。又知，B、C交点对应于 x4= 0，x5= 0，在等式约束中令x4 = 0，x5 = 0，得到x1 =7.5，x2 = 25，x3 = -7.5。即B、C交点对应于点 x = (x1 ，x2 ，x3 ， x4, ， x5)T=(-7.5，25，-7.5，0，0)T。由于有负分量，因此B、C交点不是可行域的极点。我们同样可以讨论其他交点的情况。,2.线性规划解的概念,77,下面讨论
30、线性规划标准形式的基、基本解、基本可行解的概念。考虑线性规划标准形式的约束条件：Ax=b，x0其中A为mn的矩阵，nm，秩(A) = m，b Rm 。在约束等式中，令n维空间的解向量：x = (x1，x2，xn)T,2.线性规划解的概念,78,中n-m个变量为零，如果剩下的m个变量在线性方程组中有唯一解，则这n个变量的值组成的向量x就对应于n维空间Rn中若干个超平面的一个交点。当这n个变量的值都是非负时，这个交点就是线性规划可行域的一个极点。根据以上分析，我们建立以下概念：（1）线性规划的基：对于线性规划的约束条件Ax=b， x0,2.线性规划解的概念,79,设B是A矩阵中的一个非奇异（可逆）
31、的mm子矩阵，则称B为线性规划的一个基。用前文的记号，A=( p1 ，p2 ，pn ) ，其中 pj=( a1j ，a2j ，amj )T Rm ，任取A中的m个线性无关列向量 pj Rm 构成矩阵 B=( pj1 ，pj2 ，pjm )。那么B为线性规划的一个基。我们称对应于基B的变量xj1 ，xj2，xjm为基变量；而其他变量称为非基变量。,2.线性规划解的概念,80,可以用矩阵来描述这些概念。设B是线性规划的一个基，则A可以表示为A= B , N x也可相应地分成xBx=xN其中xB为m维列向量，它的各分量称为基变量，与基B的列向量对应；xN为n-m列向量，它的各分量称为非基变量，与非基
32、矩阵N的列向量对应。这时约束等式Ax=b可表示为,2.线性规划解的概念,81,xBB,N = bxN或BxB + NxN = b如果对非基变量xN取确定的值，则xB有唯一的值与之对应xB = B-1b - B-1NxN特别，当取xN = 0，这时有xB=B-1b。关于这类特别的解，有以下概念。,2.线性规划解的概念,82,（2）线性规划问题的基本解、基本可行解和可行基：对于线性规划问题，设矩阵B = ( pj1，pj2，pjm ) 为一个基，令所有非基变量为零，可以得到m个关于基变量xj1 ，xj2 ，xjm的线性方程，解这个线性方程组得到基变量的值。我们称这个解为一个基本解；若得到的基变量的
33、值均非负，则称为基本可行解，同时称这个基B为可行基。,2.线性规划解的概念,83,矩阵描述为，对于线性规划的解xB B-1bx= =xN 0称为线性规划与基B对应的基本解。若其中B-1b0，则称以上的基本解为一基本可行解，相应的基B称为可行基。,2.线性规划解的概念,84,我们可以证明以下结论：线性规划的基本可行解就是可行域的极点。这个结论被称为线性规划的基本定理，它的重要性在于把可行域的极点这一几何概念与基本可行解这一代数概念联系起来，因而可以通过求基本可行解的线性代数的方法来得到可行域的一切极点，从而有可能进一步获得最优极点。,2.线性规划解的概念,85,例2.9: 考虑例2.8的线性规划
34、模型Max z = 1500 x1 + 2500 x2s.t. 3 x1 + 2 x2 + x3 = 652 x1 + x2 + x4 = 403 x2 + x5 = 75x1 , x2 , x3 , x4 , x5 0 注意，线性规划的基本解、基本可行 解（极点）和可行基只与线性规划问题标准形式的约束条件有关。,2.线性规划解的概念,86,3 2 1 0 0 A = P1 ,P2 ,P3 ,P4 ,P5 = 2 1 0 1 00 3 0 0 1A矩阵包含以下10个33的子矩阵：B1=p1 ，p2 ，p3 B2=p1 ，p2 ，p4B3=p1 ，p2 ，p5 B4=p1 ，p3 ，p4 B5=
35、p1 ，p3 ，p5 B6=p1 ，p4 ，p5 B7=p2 ，p3 ，p4 B8=p2 ，p3 ，p5 B9=p2 ，p4 ，p5 B10=p3 ，p4 ，p5,2.线性规划解的概念,87,其中B4= 0，因而B4不是该线性规划问题的基。其余均为非奇异方阵，因此该问题共有9个基。对于基B3=p1 ，p2 ，p5，令非基变量x3 = 0， x4 = 0，在等式约束中令x3 = 0，x4 = 0，解线性方程组：3 x1 + 2 x2 + 0 x5 = 652 x1 + x2 + 0 x5 = 400 x1 + 3 x2 + x5 = 75得到x1 =15，x2 = 10，x5 = 45，对应的基
36、本可行解：x=(x1 ，x2 ，x3 ，x4 ，x5)T=(15，10，0，0，45)T。于是对应的基B3是一个可行基。,2.线性规划解的概念,88,类似可得到x(2) = (5,25,0,5,0)T （对应B2）x(7) = (20,0,5,0,75)T （对应B5）x(8) = (0,25,15,15,0)T （对应B7）x(9) = (0,0,65,40,75)T （对应B10）是基本可行解；而x(3)= (0,32.5,0,7.5,-22.5)T（对应B9）x(4)= (65/3,0,0,-10/3,75)T （对应B6）x(5)= (7.5,25,-7.5,0,0)T （对应B1）x
37、(6) = (0,40,-15,0,-45)T （对应B8）是基本解。,2.线性规划解的概念,89,因此，对应基本可行解（极点） 的B2 B3 B5 B7 B10都是可行基。这里指出了一种求解线性规划问题的可能途径，就是先确定线性规划问题的基，如果是可行基，则计算相应的基本可行解以及相应解的目标函数值。由于基的个数是有限的（最多个），因此必定可以从有限个基本可行解中找到最优解。,2.线性规划解的概念,90,利用求解线性规划问题基本可行解（极点）的方法来求解较大规模的问题是不可行的。单纯形法的基本思路是有选择地取基本可行解，即是从可行域的一个极点出发，沿着可行域的边界移到另一个相邻的极点，要求新
38、极点的目标函数值不比原目标函数值差。,3.单 纯 形 法,91,由上节的讨论可知，对于线性规划的一个基，当非基变量确定以后，基变量和目标函数的值也随之确定。因此，一个基本可行解向另一个基本可行解的移动，以及移动时基变量和目标函数值的变化，可以分别由基变量和目标函数用非基变量的表达式来表示。同时，当可行解从可行域的一个极点沿着可行域的边界移动到一个相邻的极点的过程中，所有非基变量中只有一个变量的值从0开始增加，而其他非基变量的值都保持0不变。,3.单 纯 形 法,92,3.单 纯 形 法,单纯形法的基本过程,93,考虑标准形式的线性规划问题： Max z = c1x1 + c2x2 + + cn
39、xn s.t. a11 x1 + a12 x2 + + a1n xn = b1 a21 x1 + a22 x2 + + a2n xn = b2 .am1 x1 + am2 x2 + + amn xn = bm x1 , x2 , , xn 0 x1 c1 b1 a11 a12a1nx2 c2 b2 a21 a22a2n x= . C= . B= . A= . . . . . . . .xn cn bn am1 am2amn,3.单 纯 形 法,94,这里，矩阵A表示为： A = ( p1 ，p2 ，pn ) ， 其中 pj = ( a1j ，a2j ，amj )T Rm。若找到一个可行基，无防
40、设 B = ( p1 ，p2 ，pm ) ，则m个基变量为 x1 , x2 , , xm，n-m个非基变量为 xm+1 ，xm+2 ，xn 。通过运算，所有的基变量都可以用非基变量来表示：,3.单 纯 形 法,95,3.单 纯 形 法,x1=b1-（a1m+1xm+1+a1m+2xm+2+a1nxn） x2=b2-（a2m+1xm+1+a2m+2xm+2+a2nxn）( 2-11 ).xm=bm-（amm+1xm+1+amm+2xm+2+amnxn） 把它们代入目标函数，得z = z+m+1xm+1+m+2xm+2+nxn ( 2-12 ) 其中j=cj-（c1a1j + c2a2j + +
41、cm amj）我们把由非基变量表示的目标函数形式称为基B相应的目标函数典式。,96,单纯形法的基本步骤可描述如下：（1）寻找一个初始的可行基和相应基本可行解（极点），确定基变量、非基变量以及基变量、非基变量（全部等于0）和目标函数的值，并将目标函数和基变量分别用非基变量表示；,3.单 纯 形 法,97,（2）在用非基变量表示的目标函数表达式（2-12）中，我们称非基变量xj的系数（或其负值）为检验数记为 j 。若 j 0，那么相应的非基变量xj，它的值从当前值0开始增加时，目标函数值随之增加。这个选定的非基变量xj称为“进基变量”，转（3）。如果任何一个非基变量的值增加都不能使目标函数值增加，
42、即所有 j 非正，则当前的基本可行解就是最优解，计算结束；,3.单 纯 形 法,98,（3）在用非基变量表示的基变量的表达式（2-11）中，观察进基变量增加时各基变量变化情况，确定基变量的值在进基变量增加过程中首先减少到0的变量xr ，满足， =minbi /aij aij 0 = br /arj 这个基变量xr称为“出基变量”。当进基变量的值增加到 时，出基变量xr的值降为0时，可行解就移动到了相邻的基本可行解（极点），转（4）。,3.单 纯 形 法,99,如果进基变量的值增加时，所有基变量的值都不减少，即所有aij 非正，则表示可行域是不封闭的，且目标函数值随进基变量的增加可以无限增加，此
43、时，不存在有限最优解，计算结束；（4）将进基变量作为新的基变量，出基变量作为新的非基变量，确定新的基、新的基本可行解和新的目标函数值。在新的基变量、非基变量的基础上重复（1）。,3.单 纯 形 法,100,例2.10：用单纯形法的基本思路解例2.8的线性规划问题Max z = 1500 x1 + 2500 x2s.t. 3 x1 + 2 x2 + x3 = 652 x1 + x2 + x4 = 403 x2 + x5 = 75x1 , x2 , x3 , x4 , x5 0,3.单 纯 形 法,101,第一次迭代： （1）取初始可行基B10= (p3 , p4 , p5)，那么x3 ，x4 ，
44、x5为基变量，x1 ，x2为非基变量。将基变量和目标函数用非基变量表示：z=1500x1+2500x2 x3 = 65 - 3 x1 - 2 x2 x4 = 40 - 2 x1 - x2 x5 = 75 - 3 x2 当非基变量x1，x2=0时，相应的基变量和目标函数值为x3=65，x4=40，x5= 75，z = 0，得到当前的基本可行解： x=(0,0,65,40,75)T，z = 0 。这个解对应于图2-7的D、E交点。,3.单 纯 形 法,102,（2）选择进基变量。在目标函数 z = 1500 x1 + 2500 x2中，非基变量x1，x2的系数都是正数，因此 x1 ，x2进基都可以
45、使目标函数z增大，但x2的系数为2500，绝对值比x1的系数1500大，因此把x2作为进基变量可以使目标函数z增加更快。选择x2为进基变量，使x2的值从0开始增加，另一个非基变量x1保持零值不变。,3.单 纯 形 法,103,（3）确定出基变量。在约束条件x3 = 65 - 3 x1 - 2 x2 x4 = 40 - 2 x1 - x2 x5 = 75 - 3 x2 中，由于进基变量x2在3个约束条件中的系数都是负数，当x2的值从0开始增加时，基变量x3 、x4 、x5的值分别从当前的值65、40和75开始减少，当x2增加到25时，x5首先下降为0成为非基变量。这时，新的基变量为x3 、x4
46、、x2 ，新的非基变量为x1 、x5 ，当前的基本可行解和目标函数值为： x = (0，25，15，15，0)T，z = 62500。这个解对应于图中的C、D交点。,3.单 纯 形 法,104,第二次迭代：（1）当前的可行基为B7 = (p2 , p3 , p4)，那么x2 ，x3 ，x4为基变量，x1 ，x5为非基变量。将基变量和目标函数用非基变量表示： z = 62500 + 1500 x1 (2500/3) x5 x2 = 25 (1/3) x5 x3 = 15 - 3 x1 + (2/3) x5 x4 = 15 - 2 x1 + (1/3) x5,3.单 纯 形 法,105,（2）选择
47、进基变量。在目标函数 z = 62500 + 1500 x1 (2500/3) x5 中，非基变量x1的系数是正数，因此 x1进基可以使目标函数z增大，于是选择x1进基，使x1的值从0开始增加, 另一个非基变量x5保持零值不变。 (3)确定出基变量。在约束条件 x2 = 25 (1/3) x5 x3 = 15 - 3 x1 + (2/3) x5 x4 = 15 - 2 x1 + (1/3) x5,3.单 纯 形 法,106,中，由于进基变量x1在两个约束条件中的系数都是负数，当x1的值从0开始增加时，基变量x3 、x4的值分别从当前的值15、15开始减少，当x1增加到5时，x3首先下降为0成为
48、非基变量。这时，新的基变量为x1 、x2 、x4 ，新的非基变量为x3 、x5 ，当前的基本可行解和目标函数值为： x = (5，25，0，5，0)T，z = 70000。这个解对应于图中的A、C交点。,3.单 纯 形 法,107,第三次迭代：（1）当前的可行基为B2 = (p1 , p2 , p4 )，那么x1 ，x2 ，x4为基变量，x3 ，x5为非基变量。将基变量和目标函数用非基变量表示： z = 70000 500 x3 500 x5 x1 = 5 (1/3) x3 + (2/9) x5 x2 = 25 (1/3) x5 x4 = 5 + (2/3) x3 (1/9) x5,3.单 纯 形 法,
