Chapter9_1_特征值和乘幂法方法课件.pptVIP

第九章 矩阵特征值问题的数值方法 (Numerical methods for eigenvalue problems);引言/*Introduction*/; ;;A的特征值也称为A的特征根. ;§2 Hermite矩阵特征值问题;; 设 是Hermite矩阵A的n个特征值， 那么; 如果A的n个特征值为 其相 应的标准酉交的特征向量为 那么有 ;二、极值定理;矩阵特征值问题的性态是很复杂的，通常分别就 单个特征值或整体特征值给出条件数进行分 析. 对于Hermite矩阵，由于其特征值问题的特殊性 质，其特征值都是良态的.下面先证明Hermite矩 阵特征值的扰动定理. ; (扰动定理 ) 设矩阵A, E, A+E都是n阶Hermite矩阵，其特征值分别为 , , , 那么,;由定理; 设矩阵A和A’=A+E都是n阶Hermite矩 阵，其 特征值分别为 和 ， 那么 ;乘幂法是适用于求一般矩阵按模最大特征值及相应特征向量的算法.;任意取定初始向量x0;故当k→∞时， xk→λ1ka1v1.;对任意给定的初始向量x0;按模最大特征值λ1及其相应的特征向量v1的乘幂法的计算公式： ;试用幂法求按模最大的特征值和相应的特征向量;通过B求λ1-p 的速度快于通过A求λ1 的速度;求按模最小特征值及相应特征向量的反幂法,又称 为反迭代法.;一、λi的近似求法; “半迭代法”的命名也由此而得.;例2;由于A是对称矩阵，做一次Rayleigh商加速 ;理论上，实对称矩阵A正交相似于以A的特征值 为对角元 的 对角阵. 问题是如何构造这样的正交 矩阵呢? Jacobi方法就是通过构造特殊的正交矩 阵 序列，通过相似变换使A的非对角线元素逐次 零化来实现对角化的. ;一、平面旋转矩阵与相似约化;当 时;构造一个相似矩阵序列，使{Ak}收敛于一个对 角阵. 其中 Qk为平面旋转矩阵，其旋转角θk由 使Ak的绝对值最大元a(k)pq=a(k)qp=0 或按列依次 使A的非对角元 零化来确定. ; 设A是n阶实对称矩阵，那么由Jacobi方法产 生的相似矩阵序列{Ak}的非对角元收敛于0. 也就 是说，{Ak}收敛于以A的特征值为对角元的对角 ;又由假设,;循环Jacobi方法必须一次又一次扫描, 才能使{Ak} 收敛于对角阵 ，计算量很大. 在实际计算中，往 往用一些特殊方法来控制扫描次数, 减少计算量. ;当Ak+1的非对角元素充分小，Qk的第 j列qj可以看 成是近似特征值a(k+1)jj相应的特征向量了. ;理论上，一个实对称矩阵正交相似于一个以其特 征值为对角元的 对角阵. 但是，经典的结果告诉 我们，一个大于4次的多项式方程不可能用有限次 四则运算 求根. 因此，我们不可能期望只用有限 次相似变换将一个实对称矩阵约化为一个对角阵. 下面先介绍将一个实对称矩阵相似约化为实对称 三对角矩阵的方法，再讨论求其特征值的对分法. ;一、相似约化为实对称三对角矩阵;二、Sturm序列的性质; pi(λ)是i阶实对称矩阵的特征多项式， 因此, {pi(λ)} (i=1,2,…,n)的根全是实根. ;三、同号数和它的应用;四、求Hermite矩阵特征值的对分法