线性规划问题的最优解
线性规划问题的最优解
引言
线性规划是运筹学的一个基本分支,其应用极其广泛,其作用以为越来越多的人所重视。线性规划主要就实际问题抽象成数学形式,即求一组变量的值,在满足一定的约束条件下,是某个目标达到最小或最大,而这些约束条件用可以用一组线性不等式或线性方程来表示。而求得目标函数的最优解尤为重要,本文就线性规划问题的最优解求解方法作出阐述,并举出实例加以强化,同时也指出了线性规划问题应用于生产与运作管理的重要性。
1.线性规划问题的最优解探讨
1.1线性规划问题的提出
考虑下面的线性规划问题的标准型: 目标函数:
minZCX (1)
约束条件:
AXb (2)
X0
其中,C(c1,c2,,cn),X(x1,x2,,xn)T,b(b1,b2,,bm)T,A(aij)mn阶矩阵。设B是A中m个线性无关的列向量构成的一个基,B(aij)mm 阶矩阵,这样将矩阵A分成两个部分,即A=(B,N),X=(XB,XN),C=CB,CN,XB,CB为基B对应的非基变量和系数,XN,XN为N对应的非基变量和系数,这样将线性规划问题改写为:
X
minZCB,CNB (3)
XB
约束条件:
XB(B,N)b
(4) XN
XX0BN
经过矩阵变换,得出关于基B的标准型如下:
minZCBB1+(CN-CBB1N)XN (5)
约束条件:
XBB1NXNB1b
(6)
X,X0BN
B1b(b'1,b2,,bm)T
B1N
a1m1a2m1
a1m2a1na2m2a2n
''
amm1amm2amn
将(5)(6)展开为:
minZcibi+
'
i1m
jm1
n
(cjciaij)xj (7)
'
i1
m
约束条件:
xi
jm1
a
n
'
ij
xjb'i ,i1,2,,m (8)
xj0 ,j1,2,,n (9)
令 Z0cibi , jcjciaij ,jm1,m2,,n ,称j为检验数。
'
'
i1
i1
m
m
1.2最优解判别准则
准则一:若 X(1)(b'1,b'2,,b'm,0,,0)T ,为对应于基B的基本可行解,且对于一切的 jm1,m2,,n ,j>0 ,则X 为线性规划问题的最优解。
证明:由(7)式可知,对任意一组可行解X(x1,x2,,xn),j>0 ,Zcjxj ,
'
T
n
j1
均有 ZZ0,但 X(1) 能使等式成立,即ZZ0 ,故 X(1) 为线性规划问题的最优解。
准则二:当j0,jm1,m2,,n ,有某一个j0,设
jm1 ,i1,2,,m ,a'im10 ,则该线性规划问题有第二个最优的基本可行解。
证明:构造一个行解 X(2) ,(8') 得:
xib'ia'im1xm1 i1,2,,m xm1 0 xj0 jm2,,n
根据 原则
b'ib'L'
min' |aim10 '
1imaim1aLm1
xm1
b'L ' , xL0
aLm1
将 X(2) 带入原目标函数(4)得:
b'L
Zcibi+(cm1-ciaim1+cLaLm1)'
aLm1i1i1,iL
m
'
m
'
'
由于 m1 cm1-ciaim10 ,故:Z
'
i1
m
i1,iL
cb
m
'
ii
+ cLbL cib'LZ0
'
i1
m
X(2)也是最优的基本可行解。
推论:若 X(1) 和 X(2) 均为最优的基本可行解,XX(1)(1)X(2) ,01 均为最优可行解。
准则三:当 j0 ,jm1,m2,,n ,有某一个 j0 ,对一切 i1,2,,m ,则该线性规划有无穷多个最优解。
证明:构造一个新解 X(3) ,由 (8')
xib'ia'im1xm1 i1,2,,m
xm1 0
xj0 jm2,,n
由于 a'im10 ,0 ,故 xi0 ,i1,2,,m 将X(1)代入原目标函数(4)得:
Z
m
i1,iL
cb
m
'
ii
+(cm1-cia'im1)
i1
m
由于:m1cm1-cia'im10
i1
故:
Z
i1,iL
cb
m
'
ii
'
+0 , cibi
i1
m
时,X(3)仍为可行解,得到无穷多可行解,而目标函数仍为 当
'
ZcibiZ0 ,即X(3)也是最优解。
i1m
以下举出一些实例,进一步说明线性规划最优解的具体求解方法:
2. 线性规划最优解的问题举例
2.1图解法求解线性规划问题
例15:求解下面的线性规划问题:
maxZ58x6,
x12x43x52x60,
x24x47x55x60, (1) x2x4x3x20,
456
3xj0,j1,2,,6.
显然 X*(x1,x2,,x6)T(0,0,20,0,0,0)T 是该线性规划问题(1)的一个最优解。
b'1b'i'
因c4c50 ,及 4min':ai40,i1,2,3'0 ,
a14ai4
'
'
b'2b'i'
5min':ai50,i1,2,3'0 ,
a25a5
可考虑如下线性规划问题:
maxWx4x5
x2x3x0145
(2)
x4x7x0452
xj0,j1,2,4,5.
易解得线性规划问题(2)的最优解为X''(x1,x2,x4,x5)T(0,0,0,0)T ,W(X'')0 , 于是可得 X*(x1,x2,,x6)T(0,0,20,0,0,0)T 是该线性规划问题(1)的唯一最优解。
例27:求解下面的线性规划问题:
maxZ5x6
x1x42x5x60
x24x442x52x625 (1) x2xx3x0
456
3xj0,j1,2,,6.
显然 X*(x1,x2,,x6)T(0,20,0,0,0,0)T 是该线性规划问题(1)的一个最优解。
b'1b'i'
因c4c50 ,及 4min':ai40,i1,2,3'0 ,
a14ai4
'
'
b'2b'i'
5min':ai50,i1,2,3'0 ,
a25a5
可考虑如下线性规划问题:
maxWx4x5xx2x0145
(2)
x2xx0453
xj0,j1,3,4,5.
易解得线性规划问题(2)有无界解,X''(x1,x3,x4,x5)T(0,3,2,1)T是该问题的一个可行解 ,W(X'')30 , 于是线性规划问题的最优解不唯一。只要取X*(x1,x2,,x6)T如下:
x10,x33s;
x225s(42421)0;x42s,x5s;x60.
那么 X'也是线性规划问题的最优解。例如,分别取s=0.5、0.25时,则(0,0,1.5,1,0.5,0)T和(0,12.5,0.75,0.5,0.25,0)T以及X*(0,25,0,0,0,0)T都是该线性规划问题(1)的最优解,其中,(0,25,0,0,0,0)T是一退化的基可行解,(0,0,1.5,1,0.5,0)T是一非退化的基可行解,而
(0,12.5,0.75,0.5,0.25,0)T解。
1TT
0,0,1.5,1,0.5,00,25,0,0,0,0是一可行解而不是基2
例34:要将两种大小不同的钢板结成A,B,C 三种规格,每张钢板可同时截得三种规格的小钢板的块数乳下表所示:
今需 A , B , C 三种规格的成品分别为 15 ,18 ,27 块,问:各截取这两种钢板多
少张可得所需三种规格成品,且使得所用钢板张数最少?
2xy15x2y18
解:需要截得第一种钢板X张,第二种钢板Y张,则x3y27 (*)
x0y0
作出可行区域图如下: 目标函数为zxy ,经过可行区域内的整点且与原点最近的直线是
xy12。它上面的整点有
(0,12)、(1,11)、(2,10)、(3,9)、(4,8)、(5,7)、(6,6)、(7,5)、(8,4)、(9,3)、(10,2)、(11,1)、(12,0),若逐一讨论其是否在可行域内比较麻烦时,只需先判断点A(
1839
,)附近的整数点是否满足条件,若满足条件,则再试附近的整55
数点;若不满足条件,则不需要再判断下去。
故此题,我们只需先判断A点附近的整数点(3,9)、(4,8),分别代入(*)式,易得它们都满足条件。故我们还需判断(2,10)、(5,7)两点,代入(*)式发现它们都不满足条件,则其余点不需要再判断。所以该题的最优解为(3,9)、(4,8)。
例43:maxS5x12x2,约束条件为:
x14
xx3
12,
5x12x210x1,x20
利用图解法求解如下:
此例中约束条件中存在和目标函数的系数成比例的约束条件,但是由于此约束条件在可行域的形成中没有发挥作用,所以此问题没有多个最优解。
图解法是求解含有两个变量的线性规划问题的一种很直观和有效的方法,所以在作出问题的可行域时,则可根据这个必要条件,若可行域中没有与目标函数平行的边界线时,则可直接判断出此问题一定没有多个最优解。
2.2单纯型法求解线性规划问题
例51:求解问题:
minZx22x3
s.t.x12x2x32
x23x3x41xxx2235xj0,j1,2,3,4,5
解:这里B(A1,A4,A5)是一个单位矩阵,且b(2,1,2)T0,故基B是可行基,
x1,x4,x5为基变量,x2,x3为非基变量,
基B对应的基本可行解为:x(2,0,0,1,2)T,其目标函数值Z00.方程组Axb已是典式,得到第一张单纯性表如下:
注意 ,第0行的元素应是将目标函数Zx22x3化成等价的方程Zx22x30后的相应元素。
检验数210 ,故当前解不是最优解,A2列中有两个元素a22,a32均为正数,取
bb12
min2,3min,1
11a22a32
故转轴元为a22,x2为进基变量,x4为出基变量。进行旋转变换后得下表:
它对应的基本可行解为x(4,1,0,0,1)T,其目标函数值为Z01.但310,仍不是最优解,此时a33为转轴元,进行旋转变换后得下表:
它对应的基本可行解为x(量0,故为最优解。
13513
,,,0,0)T,其目标函数值为Z0.此时检验数向2222
例61:接下面的LP问题:
minZx1x2x3
s.t.3xxx1123
x14x2x32x1,x2,x30
解:引进非负的剩余变量x40,x50将不等式约束化为等式约束
3x1x2x3x41
x14x2x3x52 x,x,x,x,x012345
若用原始单纯形法求解,需再引进两个非负的人工变量,然后利用两阶段法求解。由本例所具有的特点,我们只要将等式两端同乘以(-1),就直接得到原问题的一个基本(不可行)解和对偶问题的一个可行解(检验数向量0)其对应的单纯形表如下:
直接利用对偶单纯形法求解。b2
2b11,所以x5为离基变量,由以下比值决定进基变量。
1121 min,
41a224
因而x2为进基变量,以a22为转轴元,进行旋转变换后得下表:
3155显然x4为离基变量,计算min,,1 a1113
444确定x1为进基变量,以a11为转轴元,进行旋转变换后得下表:
此时,b0,故原问题的最优解为x(
由以上例题可见,在某些情况下使用对偶单纯形法比用原始单纯形法更具优越性。
279,,0,0,0)T,其最优解为。 131313
参考文献
[1] 胡运权等.运筹学基础及应用[M].第五版 北京:高等教育出版社,2008-06.
[2] 管梅谷,郑汉鼎.线性规划[M].济南:山东科学技术出版社,1983. [3] 张香云.线性规划[M].浙江大学出版社,2009-02.
[4] 利奥尼德.尼森.瓦泽斯坦,克里斯托弗.卡特里尔.伯恩.线性规划导论[M].机械工业出版社, 2005-06.
[5] 卢开澄,卢华明.线性规划[M].清华大学出版社,2009-02.
[6] 江道琪,何建坤,陈松华.编.使用线性规划方法及其支持系统[M].清华大学出版社,2006-04. [7] 卢向华,侯定丕,魏权龄.运筹学教程[M].北京高等教育出版社,1989.
致 谢
经过三四个月的忙碌和学习,本次毕业设计已经接近尾声,作为一个本科生的毕业设计,由于经验的匮乏,难免有许多考虑不周全的地方,如果没有导师的督促指导,以及一起工作的同学们的支持,想要完成这个设计是难以想象的。
在这里,首先要感谢我的导师马晓娜老师。马老师平日里工作繁多,但在我做毕业设计的每个阶段,从外出实习到查阅资料,设计草案的确定和修改,中期检查,后期工作等整个过程中,都给了我悉心的指导。她的治学严谨和科学研究的精神也是我永远学习的榜样,并且积极影响我以后的学习和工作。
其次还要感谢大学五年来所有的老师,是他们为我们打下了数学专业专业知识的基础;同时还要感谢所有的同学们,正是因为有了你们的支持和鼓励。此次毕业设计才会顺利完成。
最后感谢宿州学院五年来对我的大力栽培!