数值分析课件 第8章1-2节.ppt

第8章 矩阵特征值问题计算 8.1 引 言 8 . 2 幂法及反幂法 8.2.2 加速方法 8.2.3 反幂法 应用幂法计算 的主特征值 及对应的特征向量时， 如果 （或 ）， 迭代向量 的各个不等于零 的分量将随 而趋向于无穷（或趋于零）. 这样,在计算机实现时就可能“溢出”. 为了克服这个缺点，就需要将迭代向量加以规范化. 设有一向量 ，将其规范化得到向量 其中 表示向量 的绝对值最大的分量， 即如果有 则 ，且 为所有绝对值最大的分量中的最小 下标. 主特征值为单特征值的条件下幂法可这样进行： 任取一初始向量 ， 构造向量序列 由（2.3）式 （2.8） （2.3） 这说明规范化向量序列收敛到主特征值对应的特征向量. 同理，可得到 收敛速度由比值 确定. 定理13 设 有 个线性无关的特征向量， 主特征值 满足 , 则有 （2.9） 非零初始向量 ， 按下述方法构造的向量序列 则对任意 例3 的主特征值和相应的特征向量. 用幂法计算 计算过程为 结果如表8-1. 和相应的特征向量的真值（8位数字）为 表8-1的结果是用8位浮点数字进行运算得到的， 的 分量值是舍入值. 于是得到 及相应的特征向量 原点平移法 一个补救的办法是采用加速收敛的方法. 引进矩阵 由前面讨论知, 应用幂法计算 的主特征值的收敛速度 主要由比值 来决定， 但当 接近于1时，收敛可能 很慢. 其中 为选择参数. 设 的特征值为 ， 则 的相应特征值为 而且 的特征向量相同. 如果要计算 的主特征值 ， 仍然是 的主特征值， 就要适当选择 使 且使 对 应用幂法，使得在计算 的主特征值 的过程中 得到加速. 这种方法通常称为原点平移法. 例4 设 有特征值 比值 . 则 的特征值为 做变换 应用幂法计算 的主特征值 的收敛速度的比值为 选择有利的 值，虽然能够使幂法得到加速，但问题 在于如何选择适当的参数 . 设 的特征值满足 （2.10） 则不管 如何， 的主特征值为 或 . 当希望计算 及 时，首先应选择 使 且使收敛速度的比值 显然，当 , 即 时 当 的特征值满足（2.10）且 能初步估计时， 就能确定 的近似值. 当希望计算 时，应选择 为最小， 这时收敛速度的比值为 使得应用幂法计算 得到加速. （2.10） 例5 的主特征值. 解 作变换 取 ， 计算矩阵 则 对 应用幂法，计算结果如表8-2. 由此得 的主特征值为 , 的主特征值 为 与例3结果比较，上述结果比例3迭代15次还好. 若迭代15次， （相应的 ). 原点位移的加速方法，是一个矩阵变换方法. 这种变换容易计算，又不破坏矩阵 的稀疏性，但 的选择依赖于对 的特征值分布的大致了解. 对应的特征向量满足 ，应用幂法计算 的主 特征值 ， 瑞利商加速 定理14 设 为对称矩阵，特征值满足 则规范化向量 的瑞利商给出 的较好的近似 证明 由（2.8）式及 （2.8） 得 （2.11） 反幂法用来计算矩阵按模最小的特征值及其特征向量， 设 为非奇异矩阵， 的特征值次序记为 相应的特征向量为 ， 则 的特征值为 对应的特征向量为 . 也可用来计算对应于一个给定近似特征值的特征向量. 物理、力学和工程技术的很多问题在数学上都归结为 求矩阵的特征值问题. 例如，振动问题（大型桥梁或建筑物的振动、机械的振动、电磁振荡等），物理学中某些临界值的确定. 这些问题都归结