注意,L±i>Li,i∈ J±如果wetake x±i→ ±∞, 我们得到了空的J±。将(206)代入(207),我们得到(K+m)未知量uα:Xβ的(K+m)方程组∈HbQαβuβ=yα(217),其中bqαβ≡ αβ+Xi∈eJbOhmiαbOhmiβξi(218)yα≡λXi∈eJbOhmiαξi[ρi- Liηi]+Xi∈J+x+ibOhmiα+Xi∈J-十、-ibOhmiα(219),所以我们有uα=Xβ∈HbQ-1αβyβ(220)式-这是BQ的倒数。注:给定ηi,J+,J,t(220)解uα-,eJ和J′。另一方面,(212)、(213)、(214)、(215)和(216)确定ηi,J+,J-,用uα表示eJ和J′。然后对整个系统进行迭代求解,在初始迭代时,取J(0)={1,…,N},因此J+(0),J-(0)和J′(0)为空,η(0)i=±1,i=1,N(221)而η(0)的先验值可以是二进制的,除非(K+m)<< N、 在某些情况下,可能会遇到收敛速度问题。然而,如果选择η(0)i=sign(ρi),i=1,N(222)那么迭代过程通常会收敛得相当快。以下技巧可以加快收敛速度。让我们这样做我∈ {1,…,N}:x-我≤ bx(s)i≤ x+i(223)NXi=1bx(s)iYia=0,a=1,m(224)设x(s+1)ibe为第(s+1)次迭代得到的解。该解满足线性约束,但可能不满足边界。乐其≡ x(s+1)i- bx(s)i(225)hi(t)≡ bx(s)i+t qi,t∈ [0,1](226)然后bx(s+1)i≡ 嗨(t)*) = bx(s)i+t*齐(227)其中t*是t的最大值,使得hi(t)满足边界。我们有:qi>0:pi≡ 闵x(s+1)i,x+i(228)气<0:pi≡ 最大值x(s+1)i,x-我(229)t*= 明皮- bx(s)iqiqi6=0,i=1,N(230)现在,在每个步骤中,我们可以定义J+和J,而不是(213)和(214)-via(J’isstill定义via(216))我∈ J+:bxi=x+i(231)我∈ J-: bxi=x-i(232),其中bxi是如上迭代计算的,我们可以取bx(0)i≡ 初始迭代时为0。