4.2乘幂法和反乘幂法.pptVIP

二、幂法的加速 因为幂法的收敛速度是线性的，而且依赖 于比值 ，当比值接近于1时，幂法收敛很慢幂法加速有多种，以下介绍两种。 1.Aitken加速方法 由乘幂法的计算过程可知，存在正常数C,当K充分大时 从而得到 因此 是线性收敛 3. 对称矩阵的Rayleigh商加速法 定义 设A对称，x ? 0，则称 为x关于A的Rayleigh商 思想：A对称，特征值?1，?2，…，?n均为实数，且存在特征向量v1，v2，…，vn为标准正交基。 设 ，a1 ? 0，则 当k充分大时，M与k无关) 注: 此比Aitken加速中的(3.2-6)更快 公式 称为Rayleigh商加速法。 其中 反幂法是计算矩阵按模最小的特征值及特征向 量的方法，也是修正特征值、求相应特征向量的 最有效的方法。 三、反幂法 4.2 幂法与反幂法 乘幂法是求方阵的最大特征值及对应特征向量的一种迭代法。，若辅助以相应的收缩技巧，则可以逐步计算出该矩阵按模由大到小的全部特征值以及特征向量。 乘幂法的优点是计算简单，容易在计算机上实现，适合大型稀疏矩阵主特征值以及相应的特征向量。 2.2.1 乘幂法 设An有n个线性无关的特征向量v1，v2，…，vn，对应的特征值?1，?2，…，?n，满足 |?1| |?2| ? …? |?n| (4.11) 1. 基本思想 任意取一个非零的初始向量 由矩阵A 构造一个向量序列 由递推公式 n个线性无关的特征向量v1，v2，…，vn 为Rn的一组基 现在需要分情况进行讨论 1.如果A有唯一的主特征值，即|?1| |?2| ? …? |?n| 必有 主特征值 即相邻两个迭代向量分量的比值收敛到主特征值. 2.如果A有不唯一的主特征值， 即|?1| = |?2| |?2| ? …? |?n| 则可以分为三种情况 （i） ?1= ?2 (ii) ?1= -?2 (iii) （i） ?1= ?2 A的主特征值为二重根，根据 故k充分大的时候， 若 有 且序列 收敛到相应的特征向量 这种特征值可以推广到r重主特征值 故k充分大的时候， 若 因出现因子 则当k变化时， 是一个摆动序列 通过计算 可以计算 从而计算出主特征值 当k无穷大时，可以用 便可作相应的特征向量 对于情况（iii） 即A的主特征值为共轭复根时，因为A为是对称 即 设 对任意的 当k充分大时 则 类似的有 再利用三角函数运算性质以及 复数表示，不难验证 令 由方程组 确定p,q （通常用最小二乘法） 得到 由此可以求出对应的特征向量 令max(x)表示向量x各分量绝对值最大者，对任取 初始向量 记 定义 一般的 其中 ，相应地取