幂法和反幂法求矩阵特征值课程

题目具体内容要求

A幂和幂求阵征随机产生一对称矩阵，对不同的原点位移和初值至少取3)分别使用幂法求计算矩阵的主特征值及主特征向量，用反幂法求计算矩阵的按模最小特征值及特征向量，并比较不同的原点位移和初值说明收敛。1.认真读题，了解问题的数学原形；2.选择合适问题求解的数值计算方法；3.设计程序并进行计算；4.对结果进行解释说明；对于幂法和反幂法求解矩阵特征值和特征向量的问题将从问题分析，算法设计和流程图，理论依据，程序及结果进行阐述该问题。题求n阶方阵A的特征值和特征向量，是实际计算中常常碰到的问题，如：采

机械构或电磁振动中的固有值问题等于阶矩阵A存在

和用方法及结

n维向量x满足（1）则阵的特征值，x为相应的特征向量。由高等代数知识可知，特征值是代数方程|+an+…+a=0（2）n的根。从表面上看，矩阵特征值与特征向量的求解问题似乎很简单，只需果

求解方程（2）的根，就能得到特征

，再解齐次方程组说

I-A）x=0（3）明

的解，就可得到相应的特征向量。上述方法对于小时是可以的。但当大时，计算工作量将以惊人的速度增大且由于计算带有误差必是精确的特征方程，自然就不必说求解方程2与（）的困难了。幂法是一种计算矩阵主特征值（矩阵按模最大的特征值）及对应特征向量的迭代方法，特别是用于大型稀疏矩阵。反幂法是计算海森伯格阵或三角阵的对应一个给定近似特征值的特征向量的有效方法之一。A

A法程1、幂法算法（1）取初始向量u(0)例如取u(0)=(1,1,）,置精度要k=1.（2）计算v

(k)

=Au

(

,m=max(v

(k)

),u

(k)

=v

(k)

/m

（3）若|m=m

，则停止计算（m作为绝对值最大特征值，u1

(k)

作为相应的特征向量）否则置k=k+1，转（2）2、反幂法算法（1）取初始向量u

(0)

（例如取u

(0)

=(1,1,…1)

）,置精度要k=1.（2）对A作LU分解，即A=LU（3）解线性方程组()(()()（4）计算m=max(v

(k)

),u

(k)

=v

(k)

/m

（5）|m=m

，则停止计算1/m作为绝对值最小特征

n

()

作为相应的特征向量）;否则置k=k+1，转（3A

A幂法流程图：开始输入A；[m,u,index]=pow(A,1e-6)k=0;m1=0v=A*u[vmax,i]=max(abs(v))m=v(i);u=v/m

m1=m;k=k+1abs(m-m1)<1e-6index=1;break;输出：m，u，index结束A

A反幂法流程图开始输入A,u,index]=pow_inv(A,1e-6)k=0;m1=0v=invA*u[vmax,i]=max(abs(v))m=v(i);u=v/m

m1=m;k=k+1abs(m-m1)<1e-6index=1;break;输出：m，u，index输入A；结束A

A法据幂法是用来确定矩阵的主特征值的一种迭代方法，也即，绝对值最大的特征值。稍微修改该方法，也可以用来确定其他特征值。幂法的一个很有用的特性是它不仅可以生成特征值，而且可以生成相应的特征向量。实际上，幂法经常用来求通过其他方法确定的特征值的特征向量。1、幂法的迭代格式与收敛性质设n阶矩阵A的特征是绝对值大小编号的，1x(i=1,2,…,n)为对的特征向量，单根，即ii1||≥…1则计算最大特征值与特征向量的迭代格式为v

(k)

=Au

(

,m=max(vk

(k)

),u

(k)

=v

(k)

/m

（1）其中max(v

(k)

)表示向量v

(k)

绝对值的最大分量。2、对于幂法的定理按式（1）计算出m和u()满足klim=1

lim

()

=

1maxx1反幂法是用来计算绝对值最小的特征值忽然相应的特征向量的方法。是对幂法的修改，可以给出更快的收敛性。1、反幂法的迭代格式与收敛性质设A是非奇异矩阵，则零不是特征值，并设特征值为||≥…1

n则按A

的特征值绝对值的大小排序，有|

1

|>|

1

|≥…≥|

1

|

对A

实行幂法，就可得A

的绝对值最大的特征值1/

相应的特征向量，n即A的绝对值最小的特征值和相应的特征向量。A

kAk由于用A代替作幂法计算，因此该方法称为反幂法，反幂法的迭代格式为v

(k)

=A

u

(

,m=max(vk

(k)

),u

(k)

=v

(k)

/m

（2）2、对于反幂法的定理按式（2）计算出的m和u()满足：limm=,limu)=kmax()n在式（2）中，需要用到A这给计算带来很大的不方便，因此，把（）式的第一式改为求解线性方程组Av

(k)

=u

(

（3）但由于在反幂法中，每一步迭代都需求解线性方程组式，迭代做了大量的重复计算，为了节省工作量，可事先把矩阵作LU分解，即A=LU所以线性方程组（3）改为

(k)

(k

(k)

k)A

A法代幂法程序，在matlab中建立一个M文件并保存。在matlab输入面板，输入

产生一4维随机矩阵设立初始向量可省2个参数程序结束。在M文件中可以通过改变m0的值改变原点位移，从而达到原点位移加速。反幂法程序设计代码：在matlab中建立一个M文件并保存。%pow_inv.mA

A在matlab输入面板，输入

产生一4维随机矩阵设立初始向量可省2个参数程序结束。在M文件中可以通过改变m0的值改变原点位移，从而达到原点位移加速。A

AA

AA

AA

AA

AU[111]

m0

0.69340.56230.69340.56230.69350.56230.69340.5623

k[134]0.001[367]0.001

0.69340.56220.69340.56230.69370.56240.69340.56230.69370.5624

A

A反幂法结果显示：在m0为0时A

AM0=0.001U=[111]A

AM0=0.1u=[111]A

AM0=0357]A

AM0=0.1u=[137]A

AM0=0.5u=[137]A

AM0=0u=[235]A

AM0=0.1u=[235]A

AM0=0.7u=[235]A

A综上，反幂法结果如下：m[111]0.38470.38470.3847[157]0.38470.10.38470.3847[245]0.70.70910.10.38470.3847

[-0.89961.00000.2726-0.2364][-0.89961.00000.2726-0.2364][-0.89961.00000.2726-0.2364][-0.89951.00000.2726-0.2364][-0.89961.00000.2726-0.2364][-0.89961.00000.2726-0.2364][-0.6962[-0.89951.00000.2726-0.2364][-0.89961.00000.2726