7.3 Jacobi运算方法讲解.ppt

Jacobi方法 理论上，实对称矩阵A正交相似于以A的特征值为对角元 的 对角阵. 问题是如何构造这样的正交矩阵呢? Jacobi方法就是通过构造特殊的正交矩阵 序列，通过相似变换使A的非对角线元素逐次零化来实现对角化的. 设 是二阶实对称矩阵，即a21=a12，其特征值为λ1，λ2. 令 使得 记 容易验证BT=B, 且 解之得: 当 时 当 时 并规定 Jacobi方法 设A是实对称矩阵，记A1=A. Jacobi方法的基本思想是用迭代格式 Ak+1=QTkAkQk , k=1,2,… 构造一个相似矩阵序列，使｛Ak｝收敛于一个对角阵. 其中 Qk为平面旋转矩阵，其旋转角θk由使Ak的绝对值 最大元a(k)pq=a(k)qp=0 或按列依次使A的非对角元 零化来确定. 定理1 设A是n阶实对称矩阵，那么由Jacobi方法产生的相似矩阵序列｛Ak｝的非对角元收敛于0. 也就是说，｛Ak｝收敛于以A的特征值为对角元的对角阵. 记 其中Ek是Ak除主对角元