最优化方法 第二章第二讲 单纯形法.ppt
《最优化方法 第二章第二讲 单纯形法.ppt》由会员分享,可在线阅读,更多相关《最优化方法 第二章第二讲 单纯形法.ppt(72页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、3 单纯形法单纯形法 p 单纯形法的一般原理单纯形法的一般原理 p 表格单纯形法表格单纯形法 p 借助人工变量求初始的基本可行解借助人工变量求初始的基本可行解1 单纯形法的一般步骤如下:单纯形法的一般步骤如下:(1 1)寻找一个寻找一个初始的基本可行解初始的基本可行解。(2 2)检查现行的基本可行解是否检查现行的基本可行解是否最优最优,如果为最优,如果为最优,则停止迭代,已找到最优解,否则转一步。则停止迭代,已找到最优解,否则转一步。(3 3)移至目标函数值有所改善的移至目标函数值有所改善的另一个基本可行解另一个基本可行解,然后转回到步骤(然后转回到步骤(2 2)。)。19471947年年,D
2、antzigDantzig提提出出的的单单纯纯形形法法把把寻寻优优的的目目标标集集中在所有基本可行解(即可行域顶点)中。中在所有基本可行解(即可行域顶点)中。2一、引例一、引例用单纯形方法求解下列线性规划:min z=6x14x2 2x1+3x2100 4x1+2x2120 x1、x20解:化为标准型 min z=6x14x2+0 x3+0 x4 2x1+3x2+x3 =100 4x1+2x2 +x4 =120 x1、x2,x3,x4 03基本解如下表基本解如下表 基基 基本解基本解 可行否可行否 目标值目标值 对应图中的点对应图中的点B1=(P1,P2)X1=(20,20,0,0)T 200
3、 BB2=(P1,P3)X2=(30,0,40,0)T 180 CB3=(P1,P4)X3=(50,0,0,80)T DB4=(P2,P3)X4=(0,60,80,0)T EB5=(P2,P4)X5=(0,100/3,0,160/3)T 400/3 AB6=(P3,P4)X6=(0,0,100,120)T 0 OABCDEOx1x22x1+3x2=1004x1+2x2 =120 min z=6x14x2 2x1+3x2+x3 =100 4x1+2x2 +x4 =120 x1、x2,x3,x4 04(1)找初始可行基:B1=(P3,P4)现成的初始可行基;此时,XB=(x3,x4)T,XN=(x
4、1,x2)T用XN表示Z和XB有:min z=6x14x2 x3 =1002x13x2 +x4 =1204x12x2令 XN=(0,0)T有 XB=(100,120)T则有:X(1)=(0,0,100,120)T为对应于基B1的基可行解。问:X(1)是否最优呢?否 称非基变量在目标函数中的系数为系数为检验数。min z=6x14x2 2x1+3x2+x3 =100 4x1+2x2 +x4 =120 x1、x2,x3,x4 0因为:x1和x2在目标函数中的系数为负,当x1,z;x2,z 。5(2)寻找可行)寻找可行基基B2,使其对应的基可行解使其对应的基可行解X(2)能使目标函数值增加。能使目标
5、函数值增加。选:x10则有:X(2)=(x1,0,x3,x4)Tmin z=6x14x2 x3 =1002x13x2 +x4 =1204x12x2要使为X(2)基可行解,x3,x4中必有一个为零,而另一个大等于零。只要取 x1=min100/2,120/4=30就有 x3=40,x4=0这样 X(2)=(30,0,40,0)T因为 P1,P3线性无关,因此,B2=(P1,P3)为一个基而 X(2)为对应于B2的基可行解此时 XB=(x1,x3)T,XN=(x2,x4)T用XN表示Z和XB有:min z=180 x2(3/2)x4 x1 =30(1/2)x2(1/4)x4 x3 =40 2 x2
6、 +(1/2)x4问:X(2)是否最优呢?否因为:x2在目标函数中的系数为负,当x2,z 。6(3)寻找可行寻找可行基基B3,使其对应的基可行解使其对应的基可行解X(3)能使目标函数值增加。能使目标函数值增加。选:x20则有:X(3)=(x1,x2,x3,0)T min z=180 x2(3/2)x4 x1 =30(1/2)x2(1/4)x4 x3 =40 2 x2 +(1/2)x4要使为X(3)基可行解,x1,x3中必有一个为零,而另一个大等于零。只要取 x2=min30/(1/2),40/2=20就有 x1=20,x3=0这样 X(3)=(20,20,0,0)T因为 P1,P2线性无关,因
7、此,B3=(P1,P2)为一个基而 X(3)为对应于B3的基可行解此时 XB=(x1,x2)T,XN=(x3,x4)T用XN表示Z和XB有:min z=200+(1/2)x3+(5/4)x4 x1 =20 +(1/4)x3 (3/8)x4 x2 =20 (1/2)x3+(1/4)x4问:X(3)是否最优呢?是,7从引例可以总结出求解过程:从引例可以总结出求解过程:(1)确定初始基及其基可行解;)确定初始基及其基可行解;(2)判断是否为最优解,是停止,否则)判断是否为最优解,是停止,否则转下一步;转下一步;(3)转换可行基,并求出相应的基可行)转换可行基,并求出相应的基可行解,使目标函数值有所改
8、进,转(解,使目标函数值有所改进,转(2)。)。8n 确定初始的基本可行解确定初始的基本可行解 确定初始的基本可行解等价于确定初始的可行基,一旦初确定初始的基本可行解等价于确定初始的可行基,一旦初始的可行基确定了,那么对应的初始基本可行解也就唯一确定始的可行基确定了,那么对应的初始基本可行解也就唯一确定.在在线线性性规规划划标标准准型型中中设设法法得得到到一一个个m m阶阶单单位位矩矩阵阵I I作作为为初初始始可行基。可行基。若在化标准形式前,若在化标准形式前,m m个约束方程都是个约束方程都是“”的形式,的形式,那么在化标准形时只需在一个约束不等式左端都加上一个那么在化标准形时只需在一个约束
9、不等式左端都加上一个松弛变量松弛变量x xn+in+i(i=12m)(i=12m)。9n判断现行的基本可行解是否最优判断现行的基本可行解是否最优 设有标准形式的线性规划问题:设有标准形式的线性规划问题:min z=cX;AX=b,X0;现假定现假定 A中存在一可行基中存在一可行基B,且且B为单位阵为单位阵这样这样 AX=b可以描述成如下形式,也就是用非基变量表示可以描述成如下形式,也就是用非基变量表示基变量基变量 x1 +a1,m+1xm+1+a1nxn=b1 x2 +a2,m+1xm+1+a2nxn=b2 xm+am,m+1xm+1+amnxn=bm即即 从上述约束方程中可以得到对应于从上述
10、约束方程中可以得到对应于基基B的基可行解的基可行解X=(b1,b2,bm,0,0)T10用非基变量表示目标函数有:令令 有有 式中式中 j为基可行解X的检验数检验数。更一般性:更一般性:11n向量形式向量形式 假如已求得一个基本可行解假如已求得一个基本可行解将这一基本可行解代入目标函数,可求得相应的目标函数值将这一基本可行解代入目标函数,可求得相应的目标函数值其中其中分别表示基变量和分别表示基变量和非基变量所对应的价值系数子向量。非基变量所对应的价值系数子向量。12要判定是否已经达到最小值,只需将要判定是否已经达到最小值,只需将代入目标函数,使目标函数用非代入目标函数,使目标函数用非基变量基变
11、量表示,即:表示,即:其中其中 称为非基变量称为非基变量N N的的检验向量检验向量,它的各个分量称为检验数。若,它的各个分量称为检验数。若N N的每一个检验数均小于的每一个检验数均小于等于等于0 0,即,即N N 0 0,那么现在的基本可行解就是最优解。,那么现在的基本可行解就是最优解。13定理定理1 1:最优解判别定理最优解判别定理 对于线性规划问题对于线性规划问题若某个基本可行解所对应的检验向量,若某个基本可行解所对应的检验向量,则这个基本可行解就是则这个基本可行解就是最优解最优解。定理定理2 2:无穷多最优解无穷多最优解判别定理判别定理 若是一个基本可行解,所对应的检验向量若是一个基本可
12、行解,所对应的检验向量,其中,其中存在一个检验数存在一个检验数m+km+k=0=0,则线性规划问题有则线性规划问题有无穷多最优解无穷多最优解。14 定理定理3 3:无最优解(无界)判别定理:无最优解(无界)判别定理 若若 是一个基本可行解,有一个检验数是一个基本可行解,有一个检验数 ,但是但是 ,则该线性规划问题无最优解。,则该线性规划问题无最优解。证:令证:令 ,则得新的可行解,则得新的可行解将上式代入将上式代入因为因为 ,故当故当+时,时,ZZ。15 如果现行的基本可行解不是最优解,即在检验向量如果现行的基本可行解不是最优解,即在检验向量 中中存存在在正正的的检检验验数数,则则需需在在原原
13、基基本本可可行行解解的的基基础上寻找一个新的基本可行解,并使目标函数值有所改善。础上寻找一个新的基本可行解,并使目标函数值有所改善。n 基本可行解的改进基本可行解的改进 具体做法是:具体做法是:先先从从检检验验数数为为正正的的非非基基变变量量中中确确定定一一个个换换入入变变量量,使使它它从从非非基基变量变成基变量(将它的值从零增至正值),变量变成基变量(将它的值从零增至正值),再再从从原原来来的的基基变变量量中中确确定定一一个个换换出出变变量量,使使它它从从基基变变量量变变成成非非基变量(将它的值从正值减至零)。基变量(将它的值从正值减至零)。由此可得一个新的基本可行解,由由此可得一个新的基本
14、可行解,由 可知,这样的变换一定能使可知,这样的变换一定能使目标函数值有所减少目标函数值有所减少。16l换入变量的确定换入变量的确定 最大减少原则最大减少原则 假设检验向量假设检验向量 ,选取选取最大正检验数最大正检验数所对应的非基变量为所对应的非基变量为换入变量换入变量,即若,即若 则选取对应的则选取对应的 为换入变量,为换入变量,由于且为最大,由于且为最大,因此当由零增至正值,因此当由零增至正值,可使目标函数值可使目标函数值 最大限度的最大限度的减少减少。17l如果确定为换入变量,方程如果确定为换入变量,方程其中为中与对应的系数列向量。其中为中与对应的系数列向量。现在需在现在需在 中确定一
15、个基变量为换中确定一个基变量为换出变量。出变量。l 换出变量的确定换出变量的确定 最小比值原则最小比值原则 l 当由零慢慢增加到某个值时,当由零慢慢增加到某个值时,的的非负性非负性可可能被打破。能被打破。为保持解的可行性,可以为保持解的可行性,可以按最小比值原则按最小比值原则确定换出变确定换出变量:量:若若则选取对应的基变量则选取对应的基变量 为换出变量。为换出变量。18单纯形表(1)19单纯形表(2)20计算步骤:(1)构造单纯形表1和2。(2)求 。(3)若 ,则此时基解就是最优解,否则转(4)。(4)若 ,则无最优解,否则转(5)。(5)求 。(6)以 为主元对初始单纯形表做换基运算得新
16、单纯形表,转(2)。214 单纯形表单纯形表例例1 求解求解LP问题问题 Z 0 1 -2 0 0 0 x1 x4 x5 1 -2 1 0 0 0 1 -3 1 0 0 1 -1 0 1 2 1 2x x1 1 x x2 2 x x3 3 x x4 4 x x5 5 RHSRHS22 Z 0 1 -2 0 0 0 x1 x4 x5 1 -2 1 0 0 0 1 -3 1 0 0 1 -1 0 1 2 1 2x x1 1 x x2 2 x x3 3 x x4 4 x x5 5 RHSRHS Z 0 0 1 -1 0 -1 x1 x2 x5 1 0 -5 2 0 0 1 -3 1 0 0 0 2
17、-1 1 4 1 1检验数检验数检验数检验数转轴元转轴元转轴元转轴元23 Z 0 0 1 -1 0 -1 x1 x2 x5 1 0 -5 2 0 0 1 -3 1 0 0 0 2 -1 1 4 1 1 Z 0 0 0 -0.5 -0.5-1.5 x1 x2 x3 1 0 0 -0.5 2.5 0 1 0 -0.5 1.5 0 0 1 -0.5 0.5 6.5 2.5 0.5 最优解:最优解:最优解:最优解:x x*=(6.5,2.5,0.5,0,0)=(6.5,2.5,0.5,0,0)T T最优值:最优值:最优值:最优值:z z*=-1.5=-1.524例例2:解解解解LPLP问题问题问题问题
18、4 单纯形表单纯形表25 Z 0 1 2 0 0 0 x1 x2 x3 1 0 0 -0.5 2.50 1 0 -0.5 1.5 0 0 1 -0.5 0.5 6.5 2.5 0.5x x1 1 x x2 2 x x3 3 x x4 4 x x5 5 Z 0 0 0 1.5 -2.5 -3.5 x1 x2 x3 1 0 0 -0.5 2.5 0 1 0 -0.5 1.5 0 0 1 -0.5 0.5 6.5 2.5 0.5 此问题无最优解此问题无最优解此问题无最优解此问题无最优解RHSRHS26例例3:解解解解LPLP问题问题问题问题4 单纯形表单纯形表27 Z 0 0 0 0 -1 18 x
19、1 x4 x2 1 0 1 0 0 0 0 3 1 -1 0 1 -3/2 0 1/2 4 6 3x x1 1 x x2 2 x x3 3 x x4 4 x x5 5 Z 0 0 0 0 -1 18 x1 x3 x2 1 0 0 -1/3 1/3 0 0 1 1/3 -1/3 0 1 0 1/2 0 2 2 6 此问题有多解。此问题有多解。此问题有多解。此问题有多解。RHSRHS28练习:1、用单纯形法求解 min z=6x14x2 2x1+3x2100 4x1+2x2120 x1、x20 2 2、用单纯形方法求解线性规划问题、用单纯形方法求解线性规划问题29因为非基变量因为非基变量x4的检验
20、数的检验数 ,由无穷多最优解判别定理,由无穷多最优解判别定理,本例的线性规划问题存在本例的线性规划问题存在无穷多最优解无穷多最优解X*=(20,20,0,0)T,Z*=200最优解最优值最优解最优值30n无最优解无最优解 无最优解与无可行解时两个不同的概念。无最优解与无可行解时两个不同的概念。无无可可行行解解是是指指原原规规划划不不存存在在可可行行解解,从从几几何何的的角角度度解释是指线性规划问题的可行域为空集;解释是指线性规划问题的可行域为空集;无无最最优优解解则则是是指指线线性性规规划划问问题题存存在在可可行行解解,但但是是可可行行解解的的目目标标函函数数达达不不到到最最优优值值,即即目目
21、标标函函数数在在可可行行域域内内可可以以趋趋于于无无穷穷小小(或或者者无无穷穷大大)。无无最最优优解解也也称称为为无有限最优解,或无界解。无有限最优解,或无界解。31p借助人工变量求初始的基本可行解借助人工变量求初始的基本可行解 假假定定我我们们下下面面研研究究的的线线性性规规划划都都是是非非退退化化的的,将将线线性性规规划划问问题题化化为为标标准准型型,如如果果约约束束方方程程组组中中包包含含有有一一个个单单位位矩阵矩阵I I,那么已经得到了一个那么已经得到了一个初始可行基初始可行基。否否则则在在约约束束方方程程组组的的左左边边加加上上若若干干个个非非负负的的人人工工变变量量,使使人人工工变
22、变量量对对应应的的系系数数列列向向量量与与其其它它变变量量的的系系数数列列向向量量共共同同构构成成一一个个单单位位矩矩阵阵。以以单单位位矩矩阵阵为为初初始始基基,即即可可求求得得一一个个初始的基本可行解初始的基本可行解。32 首先将原问题化为标准型首先将原问题化为标准型添加人工变量,并在目标函数中分别赋予添加人工变量,并在目标函数中分别赋予 例例4:4:线性规划问题线性规划问题:33加加入入人人工工变变量量后后的的约约束束方方程程组组与与原原约约束束方方程程组是组是不等价不等价的。的。加加上上人人工工变变量量以以后后,线线性性规规划划的的基基本本可可行行解解不不一一定定是是原线性规划的问题的基
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 最优化方法 第二章第二讲 单纯形法 优化 方法 第二 单纯
限制150内