求解对称矩阵特征问题的精化Arnoldi方法.pdfVIP

第 49卷 第 1期 厦 门大学学报 (自然科学版) Vo1．49 No．1 2010年 1月 JournalofXiamenUniversity(NaturalScience) Jan．2010 求解对称矩阵特征 问题的精化 Arnoldi方法 陈桂芝 ，叶莉瑛 (1．厦门大学数学科学学院，福建 厦门 361005；2．集美大学理学院，福建 厦门361021) 摘要 ：研究在有限精度下，如何用精化 Arnoldi方法求对称矩阵的一组正交程度可达到机器精度的近似特征向量组．首 先给出精化 Ritz向量的一个新的表达式。该表达式表明理论上对不同的近似特征值，一般地无法保证精化Arnoldi方法 所确定的精化 Ritz向量组是正交的．进一步，采用再正交化方法便可得到一组正交化程度可达到机器精度 的标准正交 近似特征向量组 ，最后的数值结果验证结论的准确性，同时再正交化后得到新的近似对的残量几乎是不变的． 关键词 ：有限精度；机器精度；特征问题；Arnoldi方法；精化Arnoldi方法 中图分类号 ：O241．6 文献标识码 ：A 文章编号 ：0438—0479(2010)一01—0011-05 精化 Arnoldi方法 [1]是一类 Krylov子空间方法 ， 不同小规模矩阵进行奇异值分解而得到 的，由此得到 其基本思想是对于指定的近似特征值，用在求解子空 的这些精化 Ritz向量之间似乎没有一个直接的联系． 问上达到最小残量的向量代替传统的Ritz向量作为 因此相对 Ritz向量 ，研究精化 Ritz向量之间的一些关 所求的特征 向量 的近似 ，新 向量称之为精化 Ritz向 系变得十分困难．本文借助一个中间量 ，给出所有精化 量．文献[2]的研究结果表明精化 Ritz向量的优越性． Ritz向量均满足 的表达式 ，通过该表达式研究精化 目前精化 Arnoldi方法 已经成为求解大规模矩 阵特征 Ritz向量组的正交性．最后的结论表 明：一般地 ，无法 问题的最有效方法之一． 保证精化 Arnoldi方法所确定的精化 Ritz向量是正交 尽管精化 Arnoldi方法的提出主要针对非对称矩 的．进一步本文同样使用再正交化方法将得到 的精化 阵的特征问题 ，但 因精化 Ritz向量 自身的最优性 ，该 Ritz向量组进行再正交．数值结果表明这样可得到一 方法对求解对称矩阵特征值 问题亦是非常有效的．本 组正交化程度可以达到机器精度 的标准正交近似特征 文主要研究用精化 Arnoldi方法求解对称矩阵特征 问 向量组 ，而且由此得到的新 的近似特征对 的残量和原 题时所得到的精化 Ritz向量的性质． 来的近似特征对的残量几乎是没有改变的． 传统的Ritz值和 Ritz向量是从求解一个小规模 1 Arnoldi方法和精化 Arnoldi方法 矩阵对应的特征对得到，因此用 Arnoldi方法求解对 称矩阵特征问题时 ，理论上如果 Arnoldi过程不 中断， 考虑矩阵特征问题 则可以得到一组不 同Ritz值所对应 的一组标准正交 A — ， 的Ritz向量组[3]．但在实际计算 中，由于舍入误差 的 其 中A为n阶对称矩阵，( ，)，：1，2，…，为A 的 影响 ，所求得不同Ritz值所对应的 Ritz向量往往数值 特征对 ll If。一1． 上并不是正交的，特别是当Ritz值 比较密集时，情况 Arnoldi方法嘲基本原理为用满足 可能会更 严重[4]．解决此类 问题 的方 法是将 得 到 的 f。f∈ (A，1) (1) Ritz向量组