幂法求解矩阵主特征值的加速方法研究.docVIP

幂法求解矩阵主特征值的加速方法 摘要：本论文主要研究的是幂法求解矩阵的主特征值和特征向量。物理、力学和工程技术中有许多需要我们求矩阵的按模最大的特征值（及称为主特征值）和特征向量。幂法是计算一个矩阵的模最大特征值和对应的特征向量的一种迭代方法。它最大的优点是方法简单，适合于大型稀疏矩阵的主特征值，但是收敛速度非常慢。所以我们要用加速的方法来加速收敛，加速方法包括原点平移加速、Rayleigh商加速和Aitken加速算法。 关键词：幂法；原点平移加速；Rayleigh商加速；Aitken加速算法 §1 引言 我们来介绍矩阵特征值和特征向量的计算方法，大家知道求一个矩阵的特征值的问题实质上是求一个多项式的根的问题，而数学上已经证明5阶以上的多项式的根一般不能用有限次运算求得。因此，矩阵特征值的计算方法本质上都是迭代，而对于大型的稀疏矩阵就需要用幂法求解最简单。但是由于收敛速度非常的慢所以我们需要用加速的方法加快收敛速度而本次论文也是针对加速问题来通过对几种方法的研究讨论。并且通过算法的实现来说明那种加速算法收敛得快，哪个计算量比较节省。其实本文主要讨论的问题是一个应用中常见的一类数值计算问题。 §2 加速算法的背景 2.1幂法 幂法是计算一个矩阵的模最大特征值和对应的特征向量的一种迭代方法。它适用于大型稀疏矩阵。为了说明其基本思想我们先假设是可对角化的，即有如下分解 其中 非奇异，再