|
前两个不等式定义了(F+2)未知量vA、u和eu的J和J′。在本小节的其余部分中,我们假设:≥ 0.我们在第5.4小节中讨论u<0情况。将(19 3)代入(192)、(169)和(170),我们得到(F+2)未知数vA、u和eu的(F+2)方程组:FXB=1QABvB+aAu+bAeu=0(197)FXA=1aAvA+cu+d eu=-1(198)FXA=1bAvA+du+eu=-eP(199)其中QAB≡ δAB+Xi∈日本脑炎Ohm依斯克拉OhmiBξi(200)aA≡xi∈JηieOhmiAξi(201)bA≡xi∈JαieOhmiAξi(202)c≡xi∈Jξi(203)d≡xi∈Jαiηiξi(204)e≡xi∈Jαiξi(205)qab,a,b=1,(F+2)是以下对称的(F+2)×(F+2)matr ix:eQAB≡ QAB(206)eQA(F+1)=eQ(F+1),A≡ aA(207)eQA,(F+2)=eQ(F+2),A≡ bA(208)等式(F+1),(F+1)≡ c(209)等式(F+1),(F+2)=等式(F+2),(F+1)≡ d(210)等式(F+2),(F+2)≡ e(211)同样,让xa和ya,a=1,(F+2)是以下(F+2)-向量:xA≡ vA(212)xF+1≡ u(213)xF+2≡ eu(214)yA=0(215)yF+1=-1(216)yF+2=-eP(217)那么我们有xA=F+2Xb=1bQabyb(218),其中bQab是矩阵的逆toeQab:bQ≡情商-1.注意,(218)求解给定ηi,J a和J′的vA,u和eu。另一方面,(194)、(195)和(196)根据vA、u和eu确定ηi、J和J′。然后对整个系统进行迭代求解,在初始迭代中,取J(0)={1,…,N},因此J′(0)为空,η(0)i=±1,i=1,N(219)而η(0)的值可以是任意的,除非F<< N、 在某些情况下,可能会遇到收敛问题。然而,如果选择η(0)i=sign(αi),i=1,N(220)那么迭代过程通常会收敛得相当快。此外,请注意,该解决方案实际上是精确的,即收敛标准由给出(根据第5.5小节,我们有u的收敛标准)≥ 这就产生了全局最优解J(s+1)=J(s)(221)我∈ J(s+1):η(s+1)i=η(s)i(222)u(s+1)=u(s)(223)eu(s+1)=eu(s)(224)A.∈ {1, . . .
|