数值分析讲义(东北大学) 第五章 矩阵特征值与特征向量的计算.pptVIP

所以 再取p=2,q=3,apq(1)=a23(1)=2.020190,类似地可得 以下依次有 * * 从而A的特征值可取为 ?1?2.125825, ?2?8.388761, ?3?4.485401 为了减少搜索非对角线绝对值最大元素时间, 对经典的Jacobi方法可作进一步改进. 1.循环Jacobi方法:按(1,2),(1,3),…,(1,n),(2,3), (2,4),…,(2,n),…,(n-1,n)的顺序, 对每个(p,q)的非零元素apq作Jacobi变换,使其零化,逐次重复扫描下去,直至?(A)?为止. 2.过关Jacobi方法: 取单调下降收敛于零的正数序列??k?,先以?1为关卡值,依照1中顺序,将绝对值超过?1的非对角元素零化,待所有非对角元素绝对值均不超过?1时,再换下一个关卡值?2 ,直到关卡值小于给定的精度? . * 练习题 第131页 习题5 5-1, 5-5, 5-7, 5-9(2), 5-10 Jacobi方法具有方法简单紧凑,精度高,收敛较快等优点, 是计算对称矩阵全部特征值和相应特征向量的有效方法,但计算量较大,一般适用于阶数不高的矩阵. 更多资料请访问：/u-459960.html * 第5章 矩阵特征值与特征向量的计算 n阶方阵A的特征值是特征方程 det(A-?E)=0 的根. Gerschgorin圆盘定理 设矩阵A=(aij)n?n ,记复平面上以aii为圆心,以ri= 为半径的n个圆盘为 Ri=??????aii??ri?,i=1,2,…,n A的特征向量是齐次线性方程组 (A-?E)x=0 的非零解. 则 (1)A的任一特征值至少位于其中一个圆盘内; (2)在m个圆盘相互连通(而与其余n-m个圆盘互不连通) 的区域内,恰有A的m个特征值(重特征值按重数记). 试讨论A的特征值的分布. 解 由A确定的3个圆盘分别为 所以 3??1?5 -2?2?2 -6??3-2 例1 设矩阵 R1=?????-4??1?, R2=???????2?, R3=?????+4??2? x y 0 -2 -4 -6 2 3 4 5 实际上， ?1=4.20308 , ?2=-0.442931 , ?3=-3.76010 * 适当选取非奇异对角矩阵D=diag(d1,d2,…,dn),则矩阵D-1AD与矩阵A有相同的特征值,且对角元素相同. 而矩阵D-1AD对应的n个圆盘为 如上例,取D=diag(2,1,1),则有 可见，R1?是独立的,所以可得3.5??1?4.5 . 对称矩阵A关于非零向量x的Rayleigh商定义为 设n阶实对称矩阵A的特征值为?1??2?…??n ,则有 * 可见,如果x是特征向量,则R(x)是对应的特征值. ?1? R(x)??n 如果x仅是一近似特征向量,则R(x)是对应的特征值的一个近似值. §1 乘幂法和反幂法 §1.1 乘幂法 乘幂法是用来求矩阵A按模最大的特征值和相应的特征向量的方法. 设A是单构矩阵, 即A有n个线性无关的特征向量. * A的n个特征值为 |?1? ? ??2? ??? ??n? 对应的特征向量为 x1,x2,…xn 线性无关. 我们要求?1 和 x1 . 乘幂法的基本思想是取初始向量v(0)?Rn,作迭代 v(k+1) =Av(k) =Ak+1v(0) , k=0,1,2,… 产生迭代序列?v(k)?. 由于x1,x2,…xn 线性无关, 从而有 v(0) =a1x1+a2x2+…+anxn 故有 v(k) = Akv(0) =a1?1kx1+a2?2kx2+…+an?nkxn (5.1) * 1. 设