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

幂法求解矩阵主特征值的加速方法摘要：本论文主要研究的是幂法求解矩阵的主特征值和特征向量。物理、力学和工程技术中有许多需要我们求矩阵的按模最大的特征值（及称为主特征值）和特征向量。幂法是计算一个矩阵的模最大特征值和对应的特征向量的一种迭代方法。它最大的优点是方法简单，适合于大型稀疏矩阵的主特征值，但是收敛速度非常慢。所以我们要用加速的方法来加速收敛，加速方法包括原点平移加速、Rayleigh商加速和Aitken加速算法。关键词：幕法；原点平移加速；Rayleigh商加速；Aitken加速算法§1引言我们来介绍矩阵特征值和特征向量的计算方法，大家知道求一个矩阵的特征值的问题实质上是求一个多项式的根的问题，而数学上已经证明5阶以上的多项式的根一般不能用有限次运算求得。因此，矩阵特征值的计算方法本质上都是迭代，而对于大型的稀疏矩阵就需要用幕法求解最简单。但是由于收敛速度非常的慢所以我们需要用加速的方法加快收敛速度而本次论文也是针对加速问题来通过对几种方法的研究讨论。并且通过算法的实现来说明那种加速算法收敛得快，哪个计算量比较节省。其实本文主要讨论的问题是一个应用中常见的一类数值计算问题。§2加速算法的背景2.1幂法为了说明其基本思想我们先假设ACnn是可_—河南理工大学数学与信息科学学院本科毕业论文~~幂法是计算一个矩阵的模最大特征值和对应的特征向量的一种迭代方法。它适用于大型稀疏矩阵。对角化的，即A有如下分解其中diag非奇异，再假定现任取一向量u0.由于X的列向量构・成。的一组基，故u可表示为0由此可知1.Akulim0k1k乂，这里nAku0C.这样，我们有inAkxjjkXjjx11kx.j土这是A的一个很好的近似特k1征向量。然而，实际计算时，这是行不通的，其原因有二：一是我们事先这表明，当10而且k充分大时，向量u并不知道A的特征值1；二是对充分大的k计算Ak的工作量太大。所以找出一种工作量较小的方法，而幂法求解收敛速度很慢所以我们还要找出一种加快速度的算法。迭代格式：yAukk1k是y的模最大分量，jk其中U°Cn是任意给定的初始向量，通常要求||U』1.定理设ARnn有n个线性无关的特征向量，主特征值满足10,按下面构造的向量序，则有||||，则对任何非零初始向量vUTOC\o"1-5"\h\z1d1rnoo列u,v・••：kkvu0,vAu,kk1max(v),k1,2,kkvu〜,•kk(1)limu——X——,kkmaxx⑵limk（注：此定理证明参阅文献［5］）4140例1.计算矩阵A5130的主特征值。1022.2幂法的应用物理，力学和工程技术中的很多问题在数学上都归结为求矩阵特征值问题。例如，振动问题（大型桥梁或建筑物的振动，机械振动，电磁振荡等），物理学中某些临界值得确定，这些问题都归结为求矩阵的特征值的数学问题。而幕法求解实际应用矩阵特征值是十分有效方法之一，但是收敛速度太慢，所以在实际应用中它所需要的时间非常的长，而且计算过程中所消耗的时间造成了实际问题的完成进度。因而我们需要通过用加速算法来加快收敛速度，让实际问题提前或者按时完成。为了加快幕法求解矩阵主特征值的收敛速度，让幕法更有效广泛的运用在实际应用生活中，我们现在就来认识几种加速方法，如原点平移法、Rayleigh商加速、Aitken加速算法、一种改进的Aitken加速算法和一种新的改进的Aitken加速算法并且对他们进行比较，看哪种加速方法收敛得快，哪种计算量比较节省等。下面我们就来说说这几种加速方法。§3常见的几种加速算法3.1原点平移法定理设ACnn,p个互不相同的特征值满足|，并且模112n最大特征值1是半单的（即1的几何重数等于它的代数重数）。如果初始向量u0在1的特征子空间上的投影不为零，则定理（2.1.2）产生的向量序列uk收敛到1的一个特征向量X1,而且由定理（2.1.2）产生的数值序列k收敛到。（注：此定理证明参阅［1］）由定理（2.1.1）可知幕法的收敛速度主要取决于l/j的大小。在定理（2.1.1）的条件下，这个数总是小于1的，它越小收敛也就越快。当它接近于1时收敛是很慢的。所以为了加快幕法的收敛速度，通常用位移的方法，即应用幕法于入I上。如果适当选取u可使AI之模最大特征值与其他特征值之模的距离更大，就可起到加速的目的。首先我们引进矩阵BApI其中p为选择参数。