奇异矩阵的Krylov方法的误差分析比较.pdf

维普资讯 第 41卷 第 5期 复 旦 学 报 (自然科学版) Vo1．4lNO．5 (…： l0，J II1，l1f1fI {}I】1【ni＼、【I『i{、(＼ l】ralS、【ic、I1(、、【1 0ct．2002 文章编号 ：0427—7104(2002)05—056604 奇异矩阵的Krylov方法的误差分析比较 周 峰 ，王 玫 (复旦大学 数学研究所 ，上海 200433) 摘 要 ：Krylov方 法是求解线性 方程组 A =b，A∈C ，b∈C 的一种迭 代方 法 ，当 A 非奇 异时 ，已有很 好 研 究．而 当A 奇异或接近奇异阵时 。在一定的假定条件下 ，Krylov方法的解与范教最小的最小二乘解 A b之间 的差是可 以估计 出来 的． 关键 词 ：Krylov方 法 ；奇异 阵；误差分 析 中图分类号 ：O2412 文献标识码 ：A Krylov方法是解线性方程组 Ar=b，A∈CN N，b∈C 的一种迭代方法 ，当 A 非奇异时，它的性质 已被很好地理解 [，而当A 奇异或接近奇异 阵时 ，Krylov方法未必会得 到解 (或最小二乘解)2． 然而在 一 定的假定条件下 ，Krylov方法的解与范数最小的最小二乘解 A b之间的差是可 以估计 出的．在文献 [1]中就 A是非奇异阵的情况作了详细的研究．本文在文献[1]的基础上将其结论推广到 A是奇异阵的 情况 ，并得到相应的结论 ．以下都假定 I5bII，=1． 1 Krylov方法 Krylor方法产生第 步迭代近似解 _Tm=VmYm Vm∈C m， Ym∈Cm， (1 其中 V = [l，…， ]构成 ，7z阶Krylov空间K (A，b) 表 1 各种 Krylov方法 比较 = span(b，A6，…，A b)的基 ， 是 在 中的系数 Tab．1 ThecomparisonofKrylovmethod 向量 ． 本文考 虑 4种 Krylov方 法 FOM ，GMRES。BiCG 和 QMR．按照构成 V 的方法 (Arnoldi和 Lanczos过程)和确 定系数 向量 Y，的方法 (Galerkin和 Minimization方法)，它 们之 间的关系见表 1HJ． 下面分析各个方法 的特点及性质 ． 先考虑求 的两种过程 ：Arnoldi过程和 Lanczos过程 它们都产生 K (A ，b)的基 和一个 ( + 1)×m 的上 Hessenberg矩阵 H ，并满足 AVm = Vm“Hm． (2) 对于 Arnoldi过程 ， l满足正交性 ，故有 ，：，+lAV =H ， 二AV =H ； (3) 其中H = [J，0]青 ．Lanczos过程在产生 V 的同时，还产生一组 K (AT，W1)的基 W 并且 V 和 收稿 日期 ：2001—06—25 基金项 目：国家 自然科学基金重点资助项 目；高等学校博士点基金资助项 目 作者简介 ：周 峰 (1978一)，男 ，硕士研究生 ． 维普资讯 第