第8章 矩阵特征值与特征向量的计算VIP免费

第8章矩阵特征值和特征向量的计算很多工程计算中，会遇到特征值和特征向量的计算，如：机械、结构或电磁振动中的固有值问题；物理学中的各种临界值等。这些特征值的计算往往意义重大。求解线性方程组的迭代法，重要一点是判断迭代法的收敛性；判断方法之一就是看迭代矩阵的特征值的模是否都小于1。PA()是的高次的多项式，它的求根是很困难的。设法通过数值方法是求它的根。通常对某个特征值，可以用些针对性的方法来求其近似值。若要求所有的特征值，则可以对A做一系列的相似变换，“收敛”到对角阵或上（下）三角阵，从而求得所有特征值的近似。n阶方阵A的特征值是特征方程PA()=det(A-E)=0的根.A的特征向量是齐次线性方程组(A-E)x=0的非零解.定理1：ARnn，1,…,n为A的特征值，则niiniiiaAtr11)(（2）A的行列式值等于全体特征值之积，即nA21)det(（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为非奇异阵，那么且为特征值，即011A.11xxA定义设矩阵A,BRnn，若有可逆阵P，使则称A与B相似。APPB1定理若矩阵A,BRnn且相似，则（1）A与B的特征值完全相同；（2）若x是B的特征向量，则Px便为A的特征向量。8.1幂法和反幂法8.1.1幂法幂法是用来求矩阵A按模最大的特征值和相应的特征向量的方法.也称为主特征值和主特征向量。设A是单构矩阵,即A有n个线性无关的特征向量.A的n个特征值为|1>2n对应的特征向量为ξ1,ξ2,…,ξn线性无关.我们要求1和ξ1.幂法的基本思想是取初始非零向量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)故有xk=Akx0=β11kξ1+β22kξ2+…+βnnkξn设|1>2n,这时,上式可写成])()([11222111nknkkknx若β10,则对充分大的k有111kkxkkkxx111111因而有nixxikik,,2,1)/()(11从而特征向量ξ1xk.乘幂法的收敛速度取决于|2/1|的大小.趋于零，时，〈范化。因为当需要对计算结果进行规在实际计算中kx1,1或上溢。从而计算时会出现下溢的非零分量趋于无穷。时当kx,11实际计算时，常把每一步计算的迭代向量xk规范化。对非零向量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,计算,3,2,1,/kmyxkkk1kkxyA)max(kkym可得)max(00xxxkkkAA))(max()(21121111niikiniikiii所以)max()max(lim111111kkx其收敛速度由比值|2/1|来确定.又由于])(max[])(max[2111211111niikiniikiii))max(max()max()max(0101xAxAAxymkkkkk所以1limkkm因此,当k充分大时可取:1mk,ξ1xk.)max(00xxxkkkAA))(max()(21121111niikiniikiii用乘幂法求A的按模最大的特征值和相应特征向量.例8.1设20101350144A解取初值x0=y0=(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,3=2.8.1.2加速技术由于)()max(121kkkoxm所以,乘幂法收敛速度取决于比值|2/1|,当|2/1|1时,收敛是很慢的.1.Aitken加速方法由上式可知0lim12111k...