关于战性方程组的m步多分裂预条件.pdfVIP

418 迦±主旦煎±墨堂查基金迨妻叁 关于战性方程组的m步多分裂预条件 刘仲云 (上海夫学数学．上海，200436) 摘要 关于求解线性方程组Axb的共轭梯度方法的基于A的一个多分裂的 预条件被考虑，其中A对称正定．本文中，多分裂的构造和理论研 究都是建立在A的正的非对角元的对角补偿约化基础上(参阅文【3】)； 预条件是针对一般对称正定矩阵而设计的，因此适合对称M一矩阵(参 阅文【41)；预条件的合理性被证明；并行操作的数值实验被给出。r 关键词 多分裂预条件对角补偿约化线性方程组 在高性能计算机上数值求解大型线性方程组 Ax--b (1) 预处理共轭梯度算法是强有力的计算工具，这里A对称正定。目前，预条件的选取 仍然是一个非常热门的研究课题。关于预条件的构造，有几种方法，其中一个常用的方法 由考虑矩阵A的一个分裂。 A=M-N f21 构成，这里．jIf非奇异且M1N的谱半径p(M1聊小于l。从该分裂可以通过应用A“的幂级 数的部分和构造预条件矩阵Pm，即 B产M(“H+…+旷1y1 (3) 这里H=M1N。它被称为基于近似逆的预条件。由(3)式可得 P21=(，一日”)A“(4) 由(4)式可推出 A=P,n-P酽l-Pm-Qk (5) 即分裂(5)被迭代矩阵}Im导出。因此Pm也被称为m步预条件。对于这种预条件，为了得 到巧1r(rJLj-给定向量)。即解方程组P2--v只需业=o做m步由分裂(2)引起的迭代 即始emlv。 如果预条件矩阵％对称正定，就说对应的m步预条件是合理的。 近来，一些基于多分裂的m步预条件被提出，参阅文献[24，7，8】。然而，所有这些预 条件的合理性都是建立在假设矩阵A为对称正定的M．矩阵或权矩阵为单位矩阵的倍数之 上。 本文给出了两种IN步多分裂预条件并证明了它们的合理性。所不同的是本文主要针 对一般对称正定矩阵且不需假设权矩阵为单位矩阵的倍数。同时，以于这两种预条件，本 文给出了数值比较。 1预备知识 我们首先回顾一些基本概念。矩阵A=(口“)∈R““称为玉矩阵，如果对于停寸，a口≤ 0。(如果)A是一非奇异二矩阵且A-1卸，则称矩阵A为胙矩阵。A的一个分裂A=M-N称 为弱正则分裂，如果M1卸且M1肥田。如果p(MIN)l，则称之为收敛分裂。 下面我们把一些后面要用的相关概念以定义的形式给出。 * k=l 非负对角矩阵，的三元矩阵组(Mk,NbEk)，k=l，…五称为A的一个多分裂。 定义2州A的一个多分裂(Mk,Nk,ED,k=1，…，l【其中地丑为方形块对角矩阵，分别表 阶和，l”的阶保持一致，则称之为块对角保形多分裂。 定义3明设A对称正定，矩阵A=曰十D被称为A的对角补偿约化矩阵。这里 B=D-R皿非负对称(约化矩阵)，D为对角矩阵(对角补偿矩阵)且满足对于任意正向量 u．Do=Ro。 2主要结果 设A对称正定，A为A的一个对角补偿约化(s如ltjes)矩阵。一旦A的一个多分裂 (肘t，ME)，k=1，…，K，被造出来时，A的相应的多分裂(蝎M最)，拓1．…，K，也随之产 生，它们之间的差别仅在于Nk=M。一A而N-蹦e-A。 k， 定理1设A对称正定，A为A的一个对角补偿约化(Stieltjes)矩阵。令A=MrN k=l，…点为弱正则分裂且^矗对称正定。假设A的多分裂(Mt，N^嘎)k=1，…，蜀，为 块对角保形。