单纯形法 之 出基入基由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“单纯形法怎么换基迭代”。
通过检验,初始可行解可能不是最优解。通过基变换得到一个新的可行基,具体做法是从可行基中换一个列向量,得到一个新的可行基,使得求解得到新的基本可行解,使目标函数值更优。为了换基就要确定换入变量与换出变量。(1)入基变量的确定
从最优解判别定理知道,当某个j0时,非基变量xj不取零值可以使目标函数值增大,故我们要选基检验数大于0的非基变量换到基变量中去。若有两个以上的,则为了是目标函数增加的更大一些,一般选j最大者的非基变量为入变量。
(2)出变量的确定
确定出基变量的方法如下。把已确定的入基变量在各约束方程中的系数除其所在约束方程中的常数项的值,把其中最小比值所在的约束方程中的原基变量确定为出基变量。
下面再进行检验其最优性,如果不是最优解还要继续进行基变换,直至找到最优解,或者能够判断出线性规划无最优解为止。
设p1,p2,....pm是一组线性独立的向量组,它们对应的基可行解是X(0)。将它代入约束方程组中pxb,x0,j1,2,...,n中得到 jjnjj1(0)xiPib
(1)i1m其他的向量Pm1,Pm2,...Pmt,...,Pn都可以用P1,P2,...Pm,线性表示,若确定非基变量pmt为换入变量,必然可以找到一组不全为0的数(i1,2,...,m)使得
Pmti,mtPi
i1或mPmti,mtPi0(2)
i1m在(2)式两边同乘一个正数,然后将它加到(1)式上,得到
m(0)xiPiPmti,mtPib i1i1(0)(xii,mt)PiP或mtb(3)i1mm当取适当值时,就能得到满足约束条件的一个可行解(即非零分量的数目不大于m个)。就应使(xi(0)i,mt)(i1,2,...,m)中的某一个为零,并保证其余的分量为非负。这个要求
x可以用以下的办法达到:比较各比值
(0)ii,mt(i1,2,...,m)。
又因必须是正数,xi(0)所以只选择i,mt0(i1,2,...,m)中比值最小的等于。以上描述用数学使表述为: xi(0)xi(0)min|i,mt0ii,mti,mt这时xi为换出变量。按最小比值确定值,称为最小比值规则。将便得到新的基可行解。
xi(0)i,mt带入X中,X1(0)(x1(0)
xi(0)l,mt(0)1,mt,...,0,...,xmxi(0)l,mtm,mt,0,...,xi(0)l,mt,0...,)
第l个分量第m+t分量 由此得到由X(0)转换到X(1)的各分量的转换公式
(0)X(1)ixi0xl,ili,mtx(0)l ,il1,mt(0)这里xi(0)是原基可行解X的各分量;Xi(1)是新基可行解X(1)的各分量;i,mt是换入向量Pmt的对应原来一组基向量的坐标。现在的问题是,这个新解X(1)的m个非零分量对应的列向量是否线性独立?事实上,因X(0)的第l个分量对应于X(1)的相应分量是零,即
Xl(0)l,mt0
其中Xl(0),均不为零,根据规则(最小比值),l,mt0。X(1)中的m个非零分量对应的m个列向量是Pj(j1,2,...,m,jl)和Pmt。若这组向量不是线性独立,则一定可以找到不全为零的数j,使
PmtjPj,jl(4)
j1m成立。又因
Pmtj,mtPj(5)
j1m将(5)式减(1)式得到
j1,j1m(j,mtj)Pjl,mtPl0
由于上式中至少有l,mt0,所以上式表明P1,P2,...,Pm是线性相关,这与假设相矛盾。由此可见,X(1)的m个非零分量对应的列向量Pj(j1,2,...,m)与Pmt是线性独立的,即经过基变换得到的解是基可行解。实际上,从一个基到另一个基可行解的变换, 就是进行一次基变换。从几何意义上讲, 就是从可行域的一个顶点转向另一个顶点。