(数值分析课件)第九章 矩阵特征值和特征向量.ppt

§3 Householder法 求一般实矩阵的全部特征值 Def 方阵B若满足：当ij+1时，bij=0,则称B为上Hessenberg阵(或准上三角阵)，即 i=j+1 i j+1 §2 平面旋转矩阵 ? 雅可比法的基本思想: 设法用一系列简单的正角阵Rk , 逐步地将 A 化为近似对角阵(非对角元近似化为0)。即选择Rk , 令 A的全部特征值 问题的关键：如何构造正交阵Rk ？ ? 平面旋转变换 ?雅可比算法: 设Ak-1 (k?1, A0 =A)未对角化，即非对角元中有较大的元素，设非对角元中按模最大的元素是 引入平面旋转矩阵 利用Rk(p,q)对Ak-1作旋转变换，使 中的非对角元 ? 应满足 常将?限制在 对Jacobi算法有几点说明： ? 1. 构造旋转矩阵时只需计算sin?,cos?,为了防止舍入误差扩大，sin?,cos?按下面公式计算： 否则， ? 2. 由于Ak 是对称阵，因此只要计算上三角（或下三角）元素即可，既节省计算量，有能保证Ak 严格对称。 ? 3. 的计算过程如下： ? 4. Ak 中经旋转变换化为零的元素，可能在Ak+1中又成为非零元素，因此不能期望通过有限次旋转变换将原矩阵A对角化，但可证 证明 Jacobi法的收敛性 由前面推论知 ? 5. 实际计算时