求解单参数特征值题目的二维Arnoldi投影算法.pdf

摘要 本篇论文主要分为三个部分，讨论了求解大规模稀疏矩阵单参数特征值问题 的二维Amoldi投影算法． 第一部分包括第一章和第二章，主要对求解大规模稀疏矩阵的特征值问题和 广义特征值问题的Krylov子空间迭代法进行了回顾，并介绍了其中的核心部分 Arnoldi过程． 第二部分包括第三章和第四章．第三章详细介绍了新出现的所谓二维 Araoldi Arnoldi过程(Two-dimensional 详细算法，并给出了重正交化的二维Arnoldi过程和相应的数值实例．第四章 详细给出了如何利用二维Amoldi过程构造投影空间的一组标准正交基，并给 Arnoldi (Two-dimensionalProjection Method(TPM))．此外，还提出了基于上述 算法的两种不同形式的显式重开始策略．随后，将此求解单参数特征值的新方 法首次应用在系统无源性的检测和强制以及动力系统的分叉问题中出现的单参 数特征值问题中，通过详尽的数值例子分析了该方法的一些性质，并与已知的 Krylov子空间迭代方法进行了比较，给出了较好的结果． 第六章给出了与求解大规模稀疏矩阵的特征值问题相关的关于Sherman． Woodbury公式求解位移后近似奇异线性方程组，反迭代法仍然可以得到十分精 确的近似特征值和特征向量，并且当A，以y是稀疏矩阵时，所花费的时间少于 LU分解求解近似奇异线性方程组的时间． 关键词： 大规模稀疏矩阵特征值问题，迭代法，投影方法，Krylov子 中文摘要 ii Morri∞n-Woodbury公式，反迭代． 中图分类号：0241．6；0175．9 Abstract Inthisth∞iswedisctlS8theTwo-dimensionalArnoldi Method Projection its for (TPM)andapplicationssolvinglargesparsesingle eigenvalueproblems． In 1and introducetheiterativemethodsfor Chapter2，we solvinglargeparse thecore eigenvalueproblems，mainlyKrylov methods，and subspace algorithm Arnoldi process． In introducenew Two-dimensionalArnoldiProcess Chapter3，we developed the toconstructtheTAPfrom algorithm standardAmoldi (TAP)and process． Wealso the forTAP exRiil- give reorthogonalizationalgorithm andnumerical themethodtoconstructasetoforthonormalba- ples．InCh