求解非对称线性方程组的GMRES法的收敛性.pdfVIP

维普资讯 第37卷第l期 大 连 理 工 大 学 学 报 v01．37，No．1 l997年 1月 JournalofDalianUniversity0fTeelmology Jan．19 97 求解非对称线性方程组的GMRES法的收敛性 l一『 陈桂芝 廉庆荣 (二)z。 (大连理工大学应用数学系 116024) 摘要 对求解大型非对称线性方程组问题，saad提出了GMRES法．在理 论方面，Saad仅对系数阵可对角化时给出了收敛性分析．本文将取消这一 限制，对系数阵A为亏损的—般情况，建立了该方法的误差估计式，并由此 说明了该方法当A非亏损阵时亦是收敛的． 分类号；O收2敛41；称；性方程；多项式；篷堡生 ～统u7一兰‘ 孝”程’ 组一 ，， ～ 一 一 一 一 懒任，qH】iE 在科学工程计算中，许多问题最终都归结为大型线性方程组的求解．由Saad提出的 GMRES(广义极小残量法)是被公认为解决大型非对称线性方程组的较为有效的方法之一． 在文[1]中，Saad对该方法作了全面的介绍，并通过数值试验证实了该方法的有效性、实用 性．但在理论上Saad仅对系数阵可对角化时给出了其收敛性的证明． 1 GMRES法 考虑线性方程组 Ax—f (1) 其中：A是Ⅳ阶大型非对称的非奇异实矩阵，f6j ，z∈ ．对给定初始向量 ，令ro=f — A觚，取Krylov子空间墨 (ro，A)=span{to，Aro，…，A一。／0)，GMRF~法为求 ∈K。(ro，A) 使得 IIr—II Ilr0一A II一 Il 一 ll (2) 《 ． 将 = + 作为式(1)的近似解． 2 若干记号及引理 设A有M个相异特征值，按实部由大到小排列为 Re(A。)≥Re(屯)≥ … ≥Re( ) (3) 其中：^的重数为dl；i=1，2，…， ． 设A为亏损阵，不失一般性，再设A非退化(仅为证明简便)，则A有Jordan分解，设为 A~SJS ；其中S为 阶非奇异阵，J=dLag(-，l，J2，．．-，Ju)． 辩麓 訾 ∞ 陈枉芝：女．1962年生．讲师 维普资讯 大 连 理 工 大 学 学 报 第37卷 ^ 1 1 - J．= ● - i= 1，2，… ，M ‘× 设Q啊为次数不超过m的所有多项式的集合，则V ∈Q 有 (A)一s (．，)s ；其中 ， P (J)=d ( (．，，)，P (JD，… ， (．，))