IMGS方法的收敛性与比较性定理.pdf

2008年4月 重庆文理学院学报(自然科学版) Apr．，2008 Journalof ofArtsand Science voL27No,2 第27卷第2期 ChongqingUniversity Sciences(NaturalEdition) IMGS方法的收敛性与比较性定理 王福，袁东锦 (扬州大学数学科学学院，江苏扬州225002) [摘要]数学、物理、流体力学和经济学中的许多问题最终都可以归结为解一个或一些大型 稀疏矩阵的线性方程组．本文给出了一种IMGS方法，在理论上证明了当系数矩阵为胁一矩阵 此方法收敛，且其渐近收敛速度要快于基本的AOR迭代法，并用数值例子验证了本文所得的 主要结论． [关键词]AOR迭代法；IMGS方法；预条件算子；胁一矩阵 考虑线性方程组 Ax=6， (1) E R“1是未知向量，bR”1是已知向量．对任一分裂A= 其中A=’(a。)…ER““是非奇异矩阵，茗E 肘一Ⅳ，其中M是非奇异矩阵，则求解(1)式的基本迭代法通常定义为 髫(M’=M一1Nx(‘’+肘一。6，，k=0，l，2，… (2) 对／g×，l的矩阵A，不失一般性，可假设A=，一三一∥，其中一厶，一u分别是A的严格下三角和严 格上三角矩阵，那么基本的AOR迭代法的迭代矩阵为： (3) +L。。，=(J—yL)一[(1一∞)，+(03一T)L+oJU]． 为了加快线性方程组(1)的迭代法的收敛速度，人们通常采用预条件的方法．现对(1．)式两端左 乘P，使(1)式等价地转化为一个预条件的方程组， PAx=Pb． (4) 其中，非奇异的矩阵P通常称作是预条件子．则与(4)式相对应的基本迭代法为： 髫‘‘“’=M；1Nex‘‘’+M；1Pb，后=0，1，2，… (5) 其中，PA=Mp—Np是蹦的分裂，且肘，是非奇异的． 这里一个比较关键且令许多作者十分感兴趣的就是如何寻找预条件子尸．1997年，Kohno等人通 过在文献[1]中取P=，+提出了一种预条件的Gauss—Seidel迭代法，其中 0 一 o 。 O O ∞O 一 詈|船 O ¨ ～ ¨ 1●●●●●●●L (6) Sa= 一 … 0 0 0 一口n—I口n—I．n 0 0 O … O 由于此方法对Gauss--Seidel迭代法做了修正并改进，故在本文中，将此方法称作是IMGS方法． 是A。的严格下三角和严格上三角矩阵，那么可得IMGS方法的迭代矩阵为： ·[收稿B瓣1]200s—OI一02 [基金项目]国家自然科学基金资助项目 【作者简介】王福(1981一)，男，甘肃人，硕士研究生，主要从事数值代数研究．E—mail：waagf=a666@sohu．COffl