1、第二章 线性规划问题及单纯形法,线性规划问题及其数学模型 图解法 单纯形法原理 单纯形法计算步骤 单纯形法的进一步讨论 数据包络分析,第一节 线性规划问题及其数学模型,线性规划在经营管理中,常常用来解决有限资源(人、财、物)的合理分配问题。在经营管理中,几乎一切问题都与有限资源的合理分配利用有关。线性规划为解决有限资源的合理分配利用提供了一个有效的数学工具。,建立线性规划数学模型是解决线性规划问题的一个重要步骤。建立的线性规划数学模型是否真正的反映客观实际,数学模型本身是否正确,都直接影响求解结果,从而影响决策结果,所以,建立正确的线性规划模型尤为重要。下面举例说明线性规划数学模型的建立。,一
2、、线性规划数学模型的建立,某厂利用A、B两种原料,生产甲、乙两种产品,有关数据如下:,例1:(产品组合问题),产品名称,甲 乙,单位产品消耗原料,原料名称,可供利用的原料数量(T/日),6 8,1 2 2 1,A B,产品售价 (千元/T),3 2,根据市场调查,有如下资料: 1.乙产品的需求量至多 2 T/日; 2.乙产品的需求量比甲产品的需求量至多大 1 T/日。 求该厂产值最大的生产方案。,提出三个问题大家考虑: 1.问题的未知数是什么? 设未知数 2.以什么准则进行决策? 目标函数 3.约束条件是什么? 约束方程,这里生产方案指的是如何安排甲、乙产品的产量。显然,产量是未知数。 设:甲
3、产品的产量为 x1 T/日乙产品的产量为 x2 T/日 决策准则是产值最大,用 Z 代表产值,则有: Z=3x1+2x2 Z 是x1、x2 的函数,称为目标函数,目标是求极大值,即:max Z= 3x1+2x2 约束条件(分三部分:资源限制、市场限制、非负限制)x1+2x262x1+x28x22x2 -x11x1,x20,约束条件,资源限制,市场限制,非负限制,整理得数学模型: 目标函数: min z = 1000 x1 + 800 x2 约束条件: s.t. x1 10.8 x1 + x2 1.6x1 2x2 1.4x1 0, x2 0,例3、配料问题(min, ),设 x1, x2分别代表
4、每粒胶丸中甲、乙两种原料的用量,某厂生产一种胶丸,已知如下资料:,例4、合理下料问题,用7.4m长的钢筋,分别截取2.9m、2.1m、1.5m各至少100根,要求用料最少。,设 xj 分别代表采用切割方案18所需7.4米的 钢筋的数量。,二、线性规划问题的共同特征, 每一个问题都用一组决策变量(x1,x2,xn)表示某一方案;这组决策变量的值就代表一个具体方案。一般这些变量取值都是非负的。 存在一定的约束条件,这些约束条件可以用一组线性等式或线性不等式来表示。 都有一个要求达到的目标,它可用决策变量的线性函数(称为目标函数)来表示。按问题的不同,要求目标函数实现最大化或最小化。,三、线性规划数
5、学模型的一般表示方式,求解线性规划问题的任务是:在满足约束条件的所有(x1,x2,xn)(可行解)中求出使目标函数达到最大(小)z 值的决策变量值(x1*,x2*,xn*)(最优解)。,1.和式,2.向量式,3.矩阵式,课堂作业:建立线性规划模型,某城市在一昼夜间,市内交通需要车辆数如图,对车辆的需求在昼夜间是变化的,车辆的工作制度是一天连续工作8小时,派车时间在各时间间隔的端点,一旦派出,就连续工作8小时。求保证需要的最小车辆数。,车辆数,时 间,0,4,7,12,16,20,24,4,8,12,4,8,12,10,8,4,派车时间在各时间间隔的端点,一旦派出,就连续工作8小时。设:各时间间
6、隔所派车辆数为xj j=1,2,6 则有:min Z=x1+x2+x3+x4+x5+x6x1+x64x1+x28x2+x3 10x3+x47x4+x512x5+x6 4x1,x2,x3,x4,x5,x6 0,第二章 线性规划问题及单纯形法,线性规划问题及其数学模型 图解法 单纯形法原理 单纯形法计算步骤 单纯形法的进一步讨论 数据包络分析,第二节 图解法,对模型中只含2个变量的线性规划问题,可以通过在平面上作图的方法求解。,一、图解法的步骤,1.等直线法,x1,x2,0,4,Q2(4,2),Q1,Q3,Q4,4x1=16,4x2=12,x1+2x2=8,3,Q2,1.建立平面直角坐标系;,4向
7、着目标函数的优化方向平移等值线,直至得到等值线与可行域 的最后交点,这种点就对应最优解。,2.找出表示每个约束的半平面,所有半平面的交集是可行域 (全体可行解的集合);,3.画出目标函数的等值线 ;,2.试算法,最优解在顶点达到: O点:X1=0, X2=0, Z=0 Q1: X1=4, X2=0, Z=8 Q2: X1=4, X2=2, Z=14 Q3: X1=2, X2=3, Z=10 Q4: X1=0, X2=3, Z=6,二、线性规划问题解的存在情况,1.存在唯一最优解,如例1,2.有无穷多最优解,若将例1目标函数变为 max z = 2x1+ 4x2,则问题变为存在无穷多最优解。如图
8、:,3. 有无界解,可行域可伸展到无穷,由此目标函数值也可增大至无穷。这种情况下问题的最优解无界。产生无界解的原因是由于在建立实际问题的数学模型时遗漏了某些必要的资源约束条件。,例如:max Z=2x1+2x2s.t. -2x1+x24x1-x2 2x1,x2 0,0,x1,x2,例如: min Z=60x1+50x22x1+4x2 803x1+2x2 60x1,x2 0,0,x1,x2,无界不一定无最有解,X1=10, x2 =15 Z=1350,模型的约束条件之间存在矛盾,建模时有错误。,4. 无可行解(可行域为空集),例如:max Z=x1+2x2-x1-x22 2x1+x24 x1,x
9、2 0,0,x1,x2,三、由图解法得到的启示,图解法虽只能用来求解只具有两个变量的线性规划问题,但它的解题思路和几何上直观得到的一些概念判断,对下面要讲的单纯形法有很大启示:1求解线性规划问题时,解的情况有:唯一最优解;无穷多最优解;无界解;无可行解。(见下页图示所示)2若线性规划问题的可行域存在,则可行域是一个凸集。3若线性规划问题的最优解存在,则最优解或最优解之一(如果有无穷多的话)一定是可行域的凸集的某个顶点。4解题思路是,先找出凸集的任一顶点,计算在顶点处的目标函数值。比较周围相邻点的目标函数值是否比这个值大,如果为否,则该顶点就是最优解的点或最优解的点之一,否则转到比这个点的目标函
10、数值更大的另一顶点,重复上述过程,一直到找出使目标函数值达到最大的顶点为止。,(d)可行域无界 (e)可行域无界 (f)可行域为空集多个最优解 目标函数无界 无可行解,(a)可行域有界 (b)可行域有界 (c)可行域无界唯一最优解 多个最优解 唯一最优解,四、线性规划问题的标准形式,(一)线性规划问题标准形式,为了使线性规划问题的解法标准,就要把一般形式化为标准形式。其一般形式如下所示:,线性规划的标准形式:,(二)线性规划问题的解法标准,1、目标函数为求极大值; 2、xj0 j=1,2,n; 3、bi0 i=1,2,m; 4、除非负约束外( xj0 ),其余 约束都为等式。,线性规划问题标准
11、形式的要求如下:,(三)标准形式的变换方法,目标函数为min型,价值系数一律反号。令 Z = -Z = -CX ,有 min Z = - max - Z = - max Z 第i 个约束的bi 为负值,则该行左右两端系数同时反号,同时不等号也要反向 第i 个约束为 型,在不等式左边增加一个非负的变量xn+i ,称为松弛变量;同时令 cn+i = 0,不等式变为等式。 第i 个约束为 型,在不等式左边减去一个非负的变量xn+i ,称为剩余变量;同时令 cn+i = 0,不等式变为等式。 若xj 0,令 xj= -xj ,代入非标准型,则有xj 0 若xj 不限,令 xj= xj - xj, xj
12、 0,xj 0,代入非标准型,(四)变换举例例1.将下述线性规划问题化为标准型:,令,其中,并按上述规则,该问题的标准形式为:,例2.将下述线性规划问题化为标准型,自己做一下练习,注意一下这几处,经过变换化为标准型如下:,x1 + x2 + x3 7x1 x2 + x3 2 3 x1+ x2 +2 x3 = 5 x1,x2 0,x3为无符号约束,例3将下述线性规划问题化为标准型min z = x1 +2x2 3x3,解:用x4 - x5 替换x3 ,令z = -z,x1 + x2 + (x4 - x5) + x6 = 7x1 x2 + (x4 - x5) - x7 = 2 3x1 + x2 +
13、2(x4 - x5) = 5x1,x2,x4,x5,x6,x7 0,max z= x1 2x2 + 3(x4 - x5)+0 x6+0 x7,用标准型求最优解后,再回到原变量。,线性代数基础知识补充与回顾,一、克莱姆规则,含有n个未知数x1,x2,xn的n个线性方程的方程组如下式所示:,克莱姆法则 如果上述线性方程组的系数行列式不等于零,即有:,那么,上述方程组有唯一解:,其中Dj(j=1,2,n)是把系数行列式D中的第j列的元素用方程组的常数项代替后得到的n阶行列式.,定理一:如果线性方程组得系数行列式D不等于零,则上述方程组一定有解,且解是唯一的。,定理二:如果上述方程组无解或有两个不同的
14、解,则它的系数行列式必为零。,二、矩阵的秩,定义1 在,矩阵A中,任取k行与k列,(K=m,k=n),位于这些行列交叉处的k的平方个元素,不改变他们在A中所处的位置次序而得到的k阶行列式,称为矩阵A的k阶子式。,定义二:设在矩阵中有一个不等于0的r阶子式D,并且所有的r+1阶子式(如果存在)全等于零,那么D称为矩阵A的最高阶非零子式,数r称为矩阵A的秩。,有了上述基本知识以后我们来看一下几个非常重要的概念,五、关于标准型解的若干基本概念,线性规划问题 :,可行解:满足上述约束条件(2.2),(2.3)的解 , 称为线性规划问题的可行解。全部可行解的集合称为可行域。 非可行解:满足约束条件(2.
15、2)但不满足非负条件(2.3)的解 X 称为非可行解 最优解:使目标函数(2.1)达到最大值的可行解称为最优解。基:设 A 为约束方程组(2.2)的 mn 阶系数矩阵,(设nm), 其秩为m,B是矩阵A中的一个mm阶的满秩子系数矩阵,称B是线性 规划问题的一个基。,不失一般性,设:,B中的每一个列向量Pj ( j1,m )称为基向量,与基向量Pj对应 的变量xj称为基变量。线性规划中除基变量以外的变量称为非基变量。基解: 在约束方程组(2.2)中,令所有非基变量xm+1xm+2xn0,又因为有 ,根据克莱姆规则,由m个 约束方程可解出m个基变量的唯一解 。将这个解 加上非基变量取0的值有 ,称
16、X为线性规划 问题的基解。显然在基解中变量取非零值的个数不大于方程数m,故基解 的总数不超过 个。 基可行解: 满足变量非负约束条件(2.3)的基解称为基可行解。可行基: 对应于基可行解的基称为可行基。退化解: 基础可行解的非零分量个数 m 时,称为退化解,例:找出下述线性规划问题的全部基解,指出其中的基可行解,并确定最优解。,解:该线性规划问题的全部基解见表l-4中的 -, 打者为基可行解,注*者为最优解,z* l9。,六、线性规划标准型问题解的关系,约束方程的 解空间,基础解,可行解,非可行解,基础 可行解,退化解,以(P3、P4、P5)作为基,令x1 = x2 =0,得到 X=(0,0,
17、8,16,12)T 为一个基可行解,对应图中O点;,2 x2 = 8x4 =16 4 x2 +x5 =12,以(P1、P2、P5)为基,令x3 = x4 =0,可得X=(4,2,0,0,4)T是基最优解,对应图中Q2点。,以(P2、P4、P5)作为基,令x1 = x3 =0,由,得X=(0,4,0,16,-4)T是个基解,不是基可行解,对应图中A点,某厂利用A、B两种原料,生产甲、乙两种产品,有关数据如下:,课堂作业:用图解法求解下列问题,产品名称,甲 乙,单位产品消耗原料,原料名称,可供利用的原料数量(T/日),6 8,2 2 1,A B,产品售价 (千元/T),3 2,根据市场调查,有如下
18、资料: 1.乙产品的需求量至多 2 T/日; 2.乙产品的需求量比甲产品的需求量至多大 1 T/日。 求该厂产值最大的生产方案。,max Z= 3x1+2x2 x1+2x262x1+x28x22x2 -x11x1,x20,0,x1,x2,X1=10/3,x2 =4/3 Z=12.67,第二章 线性规划问题及单纯形法,线性规划问题及其数学模型 图解法 单纯形法原理 单纯形法计算步骤 单纯形法的进一步讨论 数据包络分析,第三节 单纯形法原理,本节重点: 凸集与顶点 线性规划基本定理 检验数的概念和计算 最优性判别 基变换(换入变量和换出变量的确定),一、线性规划问题的几何意义,凸组合的概念 凸集的
19、概念 顶点 线性规划基本定理,二维空间,两点连线上的任何一点都是这两点的凸组合,1.凸组合 设 X1,X2,XmC,若存在 , 0 1,且 ,使 则称X 为X1,X2,Xm 的凸组合。,X2 X X1,令,2. 凸集,对简单的几何形体可以直观地判断其凹凸性,但在高维空间,只能给出点集的解析表达式,因此只能用数学解析式判断。凸集的概念为:如果集合C中任意两个点X1,X2,其连线上的所有点也都是集合C中的点,称C为凸集。由于X1,X2的连线可表示为,因此凸集定义用数学解析式可表为:对任何,有 则称C为凸集.,图中红粗线和红点是顶点。,3.顶点凸集C中满足下述条件的点X称为顶点。如果C中不存在任何两
20、个不同的点X1,X2,使X成为这两个点连线上的一个点,或者:对任何 ,不存在 ,则称X是凸集C的顶点。,4. 线性规划基本定理,定理1 若线性规划问题存在 可行解,则问题的可行域是凸集。,证 (方法1) 若满足线性规划约束条件 的所有点组成的几何图形C是凸集,根据凸集定义,C内任意两点Xl,X2连线上的点也必然在C内,下面给予证明。,设 为C内任意两点, 即 ,将X1,X2代入约束条件有,(2.4),X1,X2连线上任意一点可以表示为:,(2.5),将式(2.4)代入式(2.5)得:,所以 。由于集合中任意两点连线上的点均在集合内,所以C为凸集。,引理线性规划问题的可行解X=(x1,x2,xn
21、)T为基可行解的充分必要条件是X 的正分量所对应的系数列向量是线性无关的。,证明: (1)必要性 由基可行解的定义可知,X为基可行解 其正分量的系数列向量线性无关。 (2)充分性 若向量 线性独立,则必有km;当km时,它们恰好构成一个基, 从而 为相应的基可行解。当是km时,则一定可以从其余列向量中找出(m-k)个与 构成一个基,其对应的解恰为X,所以据定义它是基可行解。,定理2 线性规划问题的基本可行解 对应线性规划问题可行域(凸集)的顶点。,则问题可以转化为证明:的正分量对应的系数列向量线性相关 在可行 域内存在两点 、 , 可以用 、 的凸组合表示。,不是基本可行解 不是可行域的顶点。
22、 不是基本可行解 的正分量对应的系数列向量线性相关。不是可行域的顶点 在可行域内存在两点 、, 可以用 、 的凸组合表示。,证明思路: 利用反证法证明。,证明:本定理需要证明X是可行域顶点X是基可行解。下面采用 反证法,即X不是可行域的顶点X不是基可行解。分两步来证明。,(2.11),由引理知:X不是基可行解,所以 线性相关,即存在一组不全 为零的数 ,使得有:,(2.12),(2.12)式乘上一个不为零的数 得:,(2.13),(2.13)+ (2.11)得:(2.11)- (2.13)得:,令,又 可以这样来选取,( 取 ,则有 ) 使得对所有 有 ,且 。,由此 ,即X不是可行域的顶点。
23、,不失一般性,假设X的前m个分量为正,故有:,(1) X不是基可行解X不是可行域的顶点。,不失一般性,设 不是可行域的顶点,因而可以找到 可行 域内另外两个不同点Y和Z,有 , 或写为: 。,(2) X不是可行域的顶点 X不是基本可行解。,因 ,故当 时 ,必有 。,因有,(2.14) - (2.15)得,(2.15),(2.14),因 不全为零,故 线性相关,即 X不是基可行解。,定理3 若线性规划问题有最优解,一定存在一个基可行解是最优解。,情形2:当 线性相关时,存在不全为零的 ,使得。 已知,且至少有一个 ,取 ,则对任意 ,有 , ;存在某个 ,有或 。且 。,情形1:当 线性无关时
24、,由引理,线性规划问题的可行解 为基可行解的充分必要条件是 的正分量所对应的系数列向量是线性无关,定理结果显然成立。,证:设 是线性规划的一个最优解, 是目标函数的最大值。,(2-17),由可行解的定义知:,(2-16),(2.16)+ (2.17)得:(2.16)- (2.17)得:,反之,则可重复上述过程,直到剩下的大于零的分量 所对应的列向量组是线性无关为止。即最后一定可以找到一个基可行解,其目标函数值等于 。问题得证。,如果这些大于零的分量 所对应的列向量 是线性无关的,则问题转化为情形,因而定理的结论成立。,由此, 。即 。 、 是线性规划问题的最优可行解,且 、 中有一个解包含有
25、个分量 。,因 为目标函数的最大值,故有:,显然 、 是线性规划问题的可行解。分别代入目标函数可得:,令,结论:线性规划问题的可行域是凸集(凸多面体),有有限多个顶点。顶点对应基可行解。 当可行域有界时,必有顶点达到目标函数 最优值。,二、单纯型法的基本思路,由定理3可知,如果线性规划问题存在最优解,一定有一个基可行解是最优解。 因此单纯形法迭代的基本思路是:先找出一个基可行解,判断其是否为最优解,如果为否,则转换到相邻的基可行解,并使目标函数值不断增大,一直找到最优解为止。,三、单纯形法迭代原理,1确定初始基可行解 对标准型的线性规划问题 :,为说明其基本步骤,先考虑A中有单位矩阵 或其变形
26、的标准型问题。,不失一般性,在推导过程中,考虑有如下 约束方程组的问题:,Z,与当前非基变量的关系,由此可知,若存在,0 (m+1,j,n),,则有,x,j,0,其他非基变量仍为零的可行解,,其目标函数值为,这说明,当前解不是最优解。若所有,0,(m+1,j,n),,则,z,0,为可行解所,能取得的目标函数最大值,说明当前解是最优解。故称,为检验数。将基变量的检验数0也视为其检验数,可得:,注意:xj 的检验数 是当z 表示为非基变量的函数时 目标函数中xj 的系数。基变量的检验数为零。,最优性判别定理:若基可行解对应的检验数 0 ( j=1,2,,n) 则此解是最优解,否则不是最优解。,例1
27、0 中 z = 2x1+3x2 , x1 ,x2为非基变量,1=20, 2=30,X(0)不是最优解。,3.基变换,求一个改进的、“相邻”的可行基,一个基变量,将变成非基变量(换出),一个非基变量将变成 基变,量(换入)。,(1),换入变量的确定,一般,当,|,0=,s,k,,取,x,k,为换入变量。,例,10,中,,s,2,s,1,,可取,x,2,为换入变量。,第k列为主元列。,第2列为主元列。,(2) 换出变量的确定,在 中,令xk0 , 而xj =0(m+1 j n,j k),要保持xi 0 ( i=1,2,,m), 即,若所有 则xk 可取无穷大,问题无最优解。,必须 Xk,于是,当
28、为换出变量。,L行为主元行,alk为主元素,x2最多取值 = min8/2,-,12/4 = 3 =x2 =3,x5 =0,故第3个约束中的 x5 是换出变量.,新的基 B(1) = (P3 ,P4 , P2 ), 新的解 X(1)=(0,3,2,16,0)T,(1) 最优性判别定理若基可行解对应的检验数 0 ( j=1,2,,n),则此解是最优解,否则不是最优解。,4. 结论,(4)当所有的 0,又对某个非基变量 , 有 这表明可以找到另一顶点(基可行解)目标函数值也达到最大。由于该两点连线上的点也属可行域内的点,且目标函数值相等,即该线性规划问题有无穷多最优解。反之,当所有非基变量 的 O
29、时,线性规划问题具有唯一最优解。,(5) 如果存在某个 O,又Pj0,表明线性规划有无界解。,1. 系统中有j种活动,它们分享有限的资源bi; 2. 进行一个单位的第j种活动,需要第 i 种资源的量为aij; 3. 一个单位的第j种活动的产出以cj表示; 4. 第j项活动的量用xj表示。 5. 机会成本Zj:表示增加一个单位的xj所引起的目标函数的下降值。 6. 价值系数cj:表示增加一个单位的xj所引起的目标函数的增加值。 7. 判别数 =cj-zj :表示增加一个单位的xj所引起的目标函数的净增值。,四、线性规划问题的经济释义,课堂作业:,有如下线性规划:,1.变成标准型; 2.确定初始基
30、可行解; 3.确定换出变量; 4.确定换入变量; 5.说出主元行、主元列、和主元素。,1.标准型如下:,2. 初始基可行解:X(0)=(0,0,0,100,120); 3.换出变量:x2 4.确定换入变量:x4 5.说出主元行L=1;主元列k=2;主元素a12=3。,第二章 线性规划问题及单纯形法,线性规划问题及其数学模型 图解法 单纯形法原理 单纯形法计算步骤 单纯形法的进一步讨论 数据包络分析,第四节 单纯形法计算步骤,本节重点: 单纯形表(特别是检验数行) 单纯形法的计算步骤,一、单纯形表,用单纯形法求解线性规划时,设计了一种专门表格,称为单纯形表。迭代计算中每找出一个新的基可行解时,就
31、重画一张单纯形表。含初始基可行解的单纯形表称为初始单纯形表,含最优解的单纯形表称为最终单纯形表。,考虑系数矩阵中有单位矩阵的情况:,单纯形表,XB列基变量; CB列基变量的价值系数(目标函数系数); cj行价值系数; b列方程组右侧常数; 列确定换入变量时的比率计算值; 底行检验数; 中间约束方程系数。,二、计算步骤,(1)找出初始可行基,确定初始基可行解,建立初始单纯形表。,(2)检验各非基变量xj的检验数,若j 0,j=m+1,n;则 已得到最优解,可停止计算,否则转入下一步。,(3)在j 0,j=m+1,n中,若有某个k对应xk的系数列向量Pk 0, 则此问题是无界解,停止计算。否则,转
32、入下一步。,(4)根据max(j 0) =k,确定xk为换入变量,按 规则计算, =minbi/aikaik0,可确定第L行的基变量为换出变量。转入下一步。,2 3 0 0 0,1 2 1 0 04 0 0 1 00 4 0 0 1,0 2 3 0 0 0,0 0 0,8 16 12,x3 x4 x5,4 - 3,2 3 0 0 0,2 1 0 1 0 -1/2,9 2 0 0 0 -3/4,0 0 3,x3 x4 x2,2 4 -,( ),3 0 1 0 0 1/4,16 4 0 0 1 0,X(0)=(0,0,8,16,12)T, z0 =0,2 3 0 0 0,2 1 0 1 0 -1/
33、2,13 0 0 -2 0 1/4,2 0 3,x1 x4 x2,- 4 12,3 0 1 0 0 1/4,8 0 0 -4 1 2,X(1)=(0,3,2,16,0)T, z1 =9,2 3 0 0 0,4 1 0 0 1/4 0,14 0 0 -1.5 -1/8 0,2 0 3,x1 x5 x2,2 0 1 1/2 -1/8 0,4 0 0 -2 1/2 1,X(2)=(2,3,0,8,0)T, z2 =13,X(3)=(4,2,0,0, 4)T, z3 =14,第二章 线性规划问题及单纯形法,线性规划问题及其数学模型 图解法 单纯形法原理 单纯形法计算步骤 单纯形法的进一步讨论 数据包络
34、分析,第五节 单纯形法的进一步讨论,本节重点: 大M法 两阶段法 解的存在情况判别,由于所添加的剩余变量的技术系数为1,不能作为初始可行基变量,为此引入一个人为的变量(注意,此时约束条件已为“=”型),以便取得初始基变量,故称为人工变量。由于人工变量在原问题的解中是不能存在的,应尽快被迭代出去,因此人工变量在目标函数中对应的价值系数应具有惩罚性,称为罚系数。罚系数的取值视解法而定两种方法:大M法和二阶段法。,一、 人工变量的引入及其解法,1. 当约束条件为“”型,引入剩余变量和人工变量,2. 大M法的求解过程 例11,答:最优解为 x1=2, x2=2, x3=0, OBJ=36,例11的单纯
35、型表迭代过程,3.大M法的一些说明,(1)人工变量被迭代出去后一般就不会再成为基变量 (2)大M法实质上与原单纯型法一样,M 可看成一个很大的常数 (3)当检验数都满足最优条件,但基变量中仍有人工变量,说明原线性规划问题无可行解 (4)大M法手算很不方便因此提出了二阶段法。,4.二阶段法的求解过程,(1)第一阶段的任务是将人工变量尽快迭代出去,从而找到一个没有人工变量的基础可行解 (2)第二阶段以第一阶段得到的基础可行解为初始解,采用原单纯型法求解 (3)若第一阶段结束时,人工变量仍在基变量中,则原问题无解 (4)为了简化计算,在第一阶段重新定义价值系数如下:,用二阶段法求解例11的第一阶段,
36、用二阶段法求解例11的第二阶段,二、单纯型法的一些具体问题,1. 关于无界解问题 可行区域不闭合(约束条件有问题) 单纯型表中入变量 xk 对应的列中所有aik0,例12 的单纯型表及其迭代过程,2.关于退化问题 退化问题的原因很复杂,当原问题存在平衡约束时 当单纯型表中同时有多个基变量可选作出变量时 退化的严重性在于可能导致死循环,克服死循环的方法有“字典序”法,3 .关于多重解问题 多个基础可行解都是最优解,这些解在同一个超平面上,且该平面与目标函数等值面平行 最优单纯型表中有非基变量的检验数为0 最优解的线性组合仍是最优解,即 X=aX1+bX2,a+b=1,例14 的单纯型表及其迭代过
37、程,4 .关于无可行解问题,约束条件互相矛盾,无可行域 单纯型表达到最优解检验条件时,人工变量仍在基变量中,例15第一阶段的单纯型表,三、单纯行法小结 如下所示:,1. 线性规划模型及其变换,根据实际问题给出数学模型,列出初始单纯形表,进行标准化,见表,变,量,x,j,0,x,j,0,x,j,无约束,不需要处理,令,x,j,=,-,x,j,; x,j,0,令,x,j,=x,j,-,x,j,;,x,j,x,j,0,约,束,条,件,b0,b0,=,不需要处理,约束条件两端同乘,-,1,加人工变量,减去剩余,(,松弛,),变量,加人工变量,加松弛变量,目,标,函,数,Max z,Min z,加入变量
38、,的系数,松弛变量,人工变量,不需要处理,令,z,=,-,z,求,M,ax z,0,-,M,分别以每一个约束条件中松弛变量或人工变量为基,变量,列出初始单纯形表,2.线性规划解的情况 解除有唯一最优解的情况外,还有如下几种情况:,无界解:(1) 可行区域不闭合(约束条件有问题)(2)单纯型表中入变量 xk 对应的列中所有aik0退化解:(1) 退化问题的原因很复杂,当原问题存在平衡约束时(2)当单纯型表中同时有多个基变量可选作出变量时(3)退化的严重性在于可能导致死循环 无穷多解: (1)多个基础可行解都是最优解,这些解在同一个超平面上,且该平面与目标函数等值面平行(2)最优单纯型表中有非基变量的检验数为0(3)最优解的线性组合仍是最优解,即 X=aX1+bX2,a+b=1 无可行解: (1) 约束条件互相矛盾,无可行域(2)单纯型表达到最优解检验条件时,人工变量仍在基变量中,3.对目标函数求极大值标准型线性规划问题,单纯形法计算步骤的框图,