数值分析 Jacobi矩阵特征值和特征向量.pptVIP

阜师院数科院 * * 3.2 Jacobi 方法 J-方法用于求实对称阵的全部特征值、特征向量. 对于实对称阵 A，必有正交阵 U，使 U T A U = D。 其中 D 是对角阵，其主对角线元 ?i 是 A 的特征值. 正交阵 U 的第 j 列是 A 的属于 ?i 的特征向量. 原理：Jacobi 方法用平面旋转对矩阵 A 做相似变换，化 A 为对角阵，进而求出特征值与特征向量. 平面旋转矩阵 对于 p ≠ q, 下面定义的 n 阶矩阵 Upq 是平面旋转矩阵。 容易验证 Upq 是正交阵。对于向量 x，Upq x 相当于把坐标轴 Oxp和 Oxq 于所在的平面内旋转角度 ? . 变换过程: 在保证相似条件下，使主对角线外元素趋于零！ 记 n 阶方阵 A = [aij], 对 A 做下面的变换： A1= UpqT AUpq, (3.12) A1 仍然是实对称阵，因为，UpqT =Upq-1 ，知 A1与 A 的特征值相同. 下面，以 4 阶矩阵为例，来计算 (3.12) U 与 A 的行、列分块记法 如下 比较元素, 有 由此，乘开各式，我们得到下面的关系式： 由此见到，矩阵 A1 的第 p 行、列与第 q 行、列中的元素发生了变化，其它行、列中的元素不变。 特别，取旋转角 ? 满足下面的条件： 则可以得到 (3