矩阵特征值与特征向量的计算2

《矩阵特征值与特征向量的计算2》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第8章矩阵特征值和特征向量的计算很多工程计算中，会遇到特征值和特征向量的计算，如：机械、结构或电磁振动中的固有值问题；物理学中的各种临界值等。这些特征值的计算往往意义重大。求解线性方程组的迭代法，重要一点是判断迭代法的收敛性；判断方法之一就是看迭代矩阵的特征值的模是否都小于1。研究生学位课程数值分析1PA()是的高次的多项式，它的求根是很困难的。设法通过数值方法是求它的根。通常对某个特征值，可以用些针对性的方法来求其近似值。若要求所有的特征值，则可以对A做一系列的相似变换，“收敛”到对角阵或上（下）三角阵，从而求得所有

2、特征值的近似。n阶方阵A的特征值是特征方程PA()=det(A-E)=0的根.A的特征向量是齐次线性方程组(A-E)x=0的非零解.研究生学位课程数值分析2定理1：ARnn，1,…,n为A的特征值，则（2）A的行列式值等于全体特征值之积，即（1）A的迹数等于特征值之和，即特征根和特征向量的基本结论。定理2设为ARnn的特征值且Ax=x，其中x不为0,则（1）c为cA的特征值（c为常数且不为0);（2）-p为A-pI的特征值，即(A-pI)x=(-p)x;（3）k为Ak的特征值；（4）设A为非奇

3、异阵，那么且为特征值，即研究生学位课程数值分析3定义设矩阵A,BRnn，若有可逆阵P，使则称A与B相似。定理若矩阵A,BRnn且相似，则（1）A与B的特征值完全相同；（2）若x是B的特征向量，则Px便为A的特征向量。研究生学位课程数值分析48.1幂法和反幂法8.1.1幂法幂法是用来求矩阵A按模最大的特征值和相应的特征向量的方法.也称为主特征值和主特征向量。设A是单构矩阵,即A有n个线性无关的特征向量.A的n个特征值为

4、1>2n对应的特征向量为ξ1,ξ2,…,ξn线性无关.我们要求1和ξ1.幂

5、法的基本思想是取初始非零向量x0Rn,作迭代xk+1=Axk=Ak+1x(0),k=0,1,2,…产生迭代序列xk.由于ξ1,ξ2,…ξn线性无关,从而有x0=β1ξ1+β2ξ2+…+βnξn(8.3)研究生学位课程数值分析5故有xk=Akx0=β11kξ1+β22kξ2+…+βnnkξn设

6、1>2n,这时,上式可写成若β10,则对充分大的k有因而有从而特征向量ξ1xk.乘幂法的收敛速度取决于

7、2/1

8、的大小.研究生学位课程数值分析6实际计算时，常把每一步计算的迭代向量xk规范化。

9、对非零向量x,用max(x)表示x的按绝对值最大的分量,称向量y=x/max(x)为向量x的规范化向量.例如,设向量x=(2,1,-5,-1)T,则max(x)=5,y=(0.4,0.2,-1,-0.2)T.可见规范化向量y总满足‖y‖=1.幂法的规范化计算公式为:任取初始向量x0=y00,计算可得研究生学位课程数值分析7所以其收敛速度由比值

10、2/1

11、来确定.又由于所以因此,当k充分大时可取:1mk,ξ1xk.研究生学位课程数值分析8用乘幂法求A的按模最大的特征值和相应特征向量.例8.1设解取初值x0=y0=

12、(1,1,1)T,计算得kmkxk0123…101112107.26.5…6.0033526.0016756.000837(1,1,1)T(1,0.8,0.1)T(1,0.75,-0.111)T(1,0.730769,-0.188034)T…………………..(1,0.714405,-0.249579)T(1,0.714346,-0.249790)T(1,0.714316,-0.249895)T可取16.000837,ξ1(1,0.714316,-0.249895)T.实际上,A的3个特征值分别为1=6,2=3,

13、3=2.研究生学位课程数值分析98.1.2加速技术由于所以,乘幂法收敛速度取决于比值

14、2/1

15、,当

16、2/1

17、1时,收敛是很慢的.1.Aitken加速方法由上式可知可见,序列mk线性收敛于1.构造Aitken序列会达到加速收敛的目的.研究生学位课程数值分析102.原点位移法作矩阵B=A-pE,则B的特征值为qi=i-p(i=1,2,…,n),而且对应的特征向量相同.如果选取p,使q1仍然是B的按模最大特征值，且满足则对B应用乘幂法可达到加速收敛的目的。程序见P170，例8.2研究生学位课程数值分析11反幂法

18、是求矩阵按模最小的特征值和相应特征向量的方法.8.1.3反幂法设A是n阶非奇异矩阵,其特征值为

19、1

20、

21、2

22、…

23、n-1

24、>

25、n

26、>0对应的特征向量为ξ1,ξ2,…,ξn,则有A-1的特征值为对应的特征向量为ξn,ξn-1,…,ξ1.要想求n和ξn只需对A-1应用乘幂法,任取初始向量x0=y