解对称线性方程组总体最小扰动方法.pdfVIP

维普资讯 南京大学学报数学半年刊 第22卷第2期 J。OURNALOFNAN。JINGUNIVERSITY Vo1．22，No．2 2005年 11月 MATHEMATICALBIQUARTERLY Nov．，2005 解对称线性方程组的总体最小扰动方法 李 欣 (黑龙江八一农垦大学数学系，大庆 163319) 戴 华 (南京航空航天大学理学院，南京 210016) 摘要 在利用 Lanczos方法求解大型对称线性方程组时，由于舍入误差的影 响，Lanczos过程易发生中断和数值不稳定．本文提出求懈对称线性方程组的总体 极小向后扰动 (TMINBACI()方法，新方法利用 Lanczos过程产生 Krylov子空间 km(A，TO)的一组基，并求 o+丘 (A，ro)中的近似解 ，使矩阵 A【，b】的向后扰动 范数 Jl[△4．，A~]IIF极小化．同时，为减少计算量和存储量，本文给出新算法的循环格 式．在迭代过程中，利用残量范数作为判断算法终止条件的缺点是，若近似值是精确 的，残量范数是小的，反之，不一定．本文利用总体向后扰动范数作为判断算法终止 条件，克服了范数作为判断算法终止条件的不足，提出了求解大型对称线性方程组的 循环总体极小向后扰动 (RTMINBACK)方法．数值实验表叽 新方法比一些旧的方 法求解大型对称线性方程组更有效，并且RTMINBACK方法适合求解病态线性方程 组． 关键词 对称线性方程组，Krylov子空间，MINRES方法，Lanczos方法，向后扰动， 病态线性方程组 中图法分类号 O241．6 黑龙江省 自然基金资助项 目，ZTA2005—30 收稿 日期：2004—11-04． e—mail：hljllxin@163．COIIl 维普资讯 316 南京大学学报数学半年刊 2005年 l1月 1 引 言 在科学计算和工程应用中，经常需要求解对称线性方程组 ．40= b 其中．4是Tt×扎的实对称矩阵，，b∈R“．求解对称正定线性方程组的最有效方法是共轭梯度 法 (ca)tlJ及其预处理技术．而求解对称不定线性方程组的有效方法有Paige和Saunders于 1975年提出的SYMMLQ方法和MINRES方法 2【】．另外，对称Lanczos方法 3【Ⅲ】也是解对称 不定线性方程组的有效方法．在理论上，Lanczos过程产生的向量组是正交向量组．但是，在实 际计算中，由于舍入误差的影响，Lanczos向量易失去正交性．为解决这个问题，同时为减少存 储量和运算量，人们采用重新开始的Lanczos方法，即循环Lanczos迭代法 【¨．一般迭代法的循 环格式是在迭代过程中用残最范数作为判断算法终止的准则．但是若近似解很接近精确解，残 量范数很小；而反过来不一定．为克服残量范数作为算法终止准则的不足，Kasenally提出了解 非对称线性方程组的GMBACK方法 [5J|曹志浩提出了解非对称线性方程组的总体GMBACK 方法 引【．本文给出一种解对称线性方程组的方法 ．总体最小扰动 (TMINBACK)方法，该方法 利用 Lanczos过程产生Krylov子空间km(A，r0)的一组基，并求近似解 Xrn∈XO+km(A，TO) 使矩阵 [A，6I的向后扰动范数 ll[△ ，△6】llF最小化． 本文第二节进行向后扰动分析，提出用总体向后扰动范数作为判断算法终止的准则； 第三节给出求解对称线性方程组的总体最小扰动 (TMINBACK)方法和循环总体最小扰动 (KTMINBACK)方法；最后给出一些数值例子，将 Lanczos方法与 TMINBACK方法、循 环Lanczos迭代法与RTMINBACK方法进行比较，并给出本文的重