广义预条件对称-反对称分裂迭代法的收敛性分析.docx

广义预条件对称-反对称分裂迭代法的收敛性分析
广义预条件对称-反对称分裂迭代法是解线性方程组的常用方法之一，它的主要思想是将系数矩阵分裂成对称部分和反对称部分，然后将方程组分解成两个方程组分别对应这两个部分，再通过迭代的方式逐步逼近解。广义预条件对称-反对称分裂迭代法在实际应用中具有高效、稳定的优点，因此一直受到广泛关注。
对于线性方程组Ax=b，可以将系数矩阵A分解成两部分A=D-L-U，其中D是A的对角线元素构成的对角矩阵，L是A的下三角矩阵与D的和，U是A的上三角矩阵与D的和。则原方程可以变形为(D-L)x=Ux+b，和(D-U)x=Lx+b。
根据对称和反对称的定义，我们可以将矩阵D-L和D-U分别定义为对称部分S和反对称部分T，则方程组可以分解为Sx=Tx+b和Tx=Sx+b。同时，我们可以得到广义预条件对称-反对称分裂迭代法的迭代方程为：
x(k+1)=[D-(w/2)S]^-1[(D+(w/2)T)x(k)+wb]
其中，w是松弛因子，b为常数向量，x(k)为第k次迭代的解向量。
接下来我们对广义预条件对称-反对称分裂迭代法的收敛性进行分析。
首先，我们需要证明方程组Sx=Tx+b和Tx=Sx+b的解存在并唯一。由于矩阵S和T分别是对称和反对称矩阵，它们的特征值具有以下性质：
-S的特征值非负，T的特征值只有实部为0时才能取到
-S的特征向量正交归一，T的特征向量正交，但不一定归一
因此，矩阵S是正定矩阵，矩阵T的非零特征值的实部都为0，并且S和T的特征向量可以相互表示。综上所述，方程组Sx=Tx+b和Tx=Sx+b的解存在且唯一。
其次，我们需要证明广义预条件对称-反对称分裂迭代法的收敛性。设x*是方程组Ax=b的解向量，即Sx*=Tx*+b。我们将x(k+1)和x*的误差向量定义为e(k+1)=x(k+1)-x*，对迭代方程进行变形，可以得到：
e(k+1)=[D-(w/2)S]^-1[(D+(w/2)T)e(k)]
再对误差向量进行变形，可以得到：
||e(k+1)||≤||[D-(w/2)S]^-1||||(D+(w/2)T)||||e(k)||
因此，我们也需要证明||[D-(w/2)S]^-1||||(D+(w/2)T)||<1，才能保证广义预条件对称-反对称分裂迭代法的收敛性。
根据矩阵范数的性质和对称矩阵的对角化定理，我们可以得到：
||[D-(w/2)S]^-1||=||[D+(w/2)(T-S)]^-1||≤||[D+(w/2)(T-S)]||^-1
由于矩阵T-S是反对称矩阵，则它的特征值满足：
||T-S||=max{|λ|：λ是T-S的特征值}=max{|λ|：λ是T的特征值}
又因为T的特征值的实部为0，故有||T-S||<||T||。因此，可以得到：
||[D+(w/2)(T-S)]||≤||D||(1+w/2)+||T-S||≤||A||(1+w/2)
综合上述结论，可以得到：
||[D-(w/2)S]^-1||||(D+(w/2)T)||≤||A||(1+w/2)/||D+S||
接下来，我们需要找到w的取值范围使得上式小于1。由于D和S均是正定矩阵，则有：
||D+S||≥||D||+||S||=2||D||
因此，我们可以得到：
||A||(1+w/2)/(2||D||)<1
化简后，可以得到：
0<w<2/(1+||D^-1||||S||/||A||)
综上所述，广义预条件对称-反对称分裂迭代法的收敛性要求存在正定矩阵D，且松弛因子w的取值范围为(0,2/(1+||D^-1||||S||/||A||))。在这个取值范围内，迭代收敛，收敛速度较快，迭代次数不会受到较大影响。