约束的奇异半定远线性方程组的迭代解法.docxVIP

约束的奇异半定远线性方程组的迭代解法 0 半收敛矩阵的解析确定 1965年，keell研究了难以置信的hetle函数。 Ax=b(1) 的迭代解法,证明了一个著名的结果:若A=M-N是P-正则分裂(即M非奇异,且M*+N为正定阵),则迭代 xj+1=M-1Nxj+M-1b(2) 是半收敛的(指迭代阵M-1N是半收敛阵)当且仅当A是半正定阵. 1997年,Ern与Giovangigli研究了求解约束方程组 Ax=b,x∈L(3) 的问题,其中A∈Rn×n是奇异的对称半正定阵,L为Rn的子空间,满足如下条件: b∈R(A) 与L?N(A)=Rn. (4) 据文报道,形如(3)的约束方程组的求解问题,有着多方面的应用背景. 文首先给出了Keller定理的一个很长的新证,然后在Keller定理的基础上,借助于迭代(2),提出了所谓的“投影迭代算法”(Projected iterative algorithms): y0=PL,N(A)x0,yj+1=PL,N(A)M-1Nyj+PL,N(A)M-1b,j=0,1,2,…, (5) 并证明了两个主要的结论: (1)ρ(PL,N(A)M-1N)1; (6) (2)y∞=limj→∞yj=limj→∞ΡL,Ν(A)xj=Ζb(7)(2)y∞=limj→∞yj=limj→∞PL,N(A)xj=Zb(7) 是(3)的惟一解,其中Z是A的对称半正定(1,2)-逆,且R(Z)=L. 本文的目标有两个:一是给出Keller定理的新证,它比文给出的新证浅显而简洁,并据之给出文关于投影迭代(5)的两个主要结论的简证;二是提出求解约束方程组(3)的两个十分简单的迭代格式.附带地,我们将给出半收敛阵的几个等价描述以及四分块三角阵半收敛的充要条件. 本文采用文中关于广义逆矩阵与投影算子的记号.另外,σ(A)表示A的谱,ρ(A)表示A的谱半径,在A仅有实特征值时,用λmax(A)与λmin(A)分别表示A的最大与最小特征值. 一如通常,若limj→∞Aj=0limj→∞Aj=0,则称A为收敛阵;若limj→∞Ajlimj→∞Aj存在(其值无需为零阵),则称A为半收敛阵. 引理0.1设A∈Cn×n.则下列提法是互相等价的: (1)A半收敛. (3)(I-A)#存在,且ρ(A(I-A)(I-A)#)1. (4)(I-A)#存在,且δ(A)≡max{|λ|:λ∈σ(A),λ≠1}1. (5)(I-A)#存在,且ρ(APL,N(I-A))1,其中L为N(I-A)的任一补子空间. 证明注意(I-A)#存在当且仅当I-A可逆或可表为,其中C可逆,并用关系式ρ(FG)=ρ(GF). (1)(I-G)#存在当且仅当如下两个条件成立: (a)A#与B#均存在,(b)(I-BB#)C(I-AA#)=0; (8) 其中W=B#2C(I-AA#)+(I-BB#)CA#2-B#CA#. (10) (2)G半收敛当且仅当如下两个条件成立: (a)I-A与I-B均半收敛,即 A#与B#均存在,ρ((I-A)AA#)1,且ρ((I-B)BB#)1; (11) (b) (I-BB#)C(I-AA#)=0, (12) 此时 其中Z=-(I-BB#)CA#-B#C(I-AA#). 证明(1) 熟知,(I-G)#存在?对某个矩阵X有I-G=(I-G)2X,即 [A0CB]=[A20CA+BCB2][X1X2X3X4][AC0B]=[A2CA+BC0B2][X1X3X2X4] ?{A=A2X1①,0=A2X2②,C=(CA+BC)X1+B2X3③?B=(CA+BC)X2+B2X4④ ①表明了A#存在.以A#左乘②的两边得AX2=0⑤ 以A#左乘①的两边得AX1=AA#⑥ 从③与⑥推出 C=CAA#+BCX1+B2X3⑦ 由⑦得C(I-AA#)=BCX1+B2X3,左乘以(I-BB#)得到(8)中所示的条件(b). 又,从⑦与⑤可得CX2=BCX1X2+B2X3X2,代入④,并注意⑤,可得 B=B2(CX1X2+BX3X2+X4). 此式即表明B#存在.以上证明了条件(8)的必要性. 反之,若条件(8)中的两个条件(a)与(b)成立,则可直接验明,(9)式右边的矩阵是(I-G)#.这证明了条件(8)的充分性. (2) 据引理0.1之(3)知,G半收敛?(I-G)#存在且ρ(G(I-G)(I-G)#)1. (1)中已证明(I-G)#存在?条件(8)成立,且此时(I-G)#的显式为(9). 不难算得 其中Q=-CAA#+(I-B)B#C(I-AA#)+(I-BB#)CA#. (14) 从(13)式易见 ρ(G(I-G)(I-G)#)1?ρ((I-A)AA#)1且ρ((I-B)BB