雅克比矩阵知识介绍.pdfVIP

.- 雅可比矩阵(Jacobi 方法) Jacobi 方法 Jacobi 方法是求对称矩阵的全部特征值以及相应的特征向量的一种方法， 它是基于以下两个结论谢谢阅读 1) 任何实对称矩阵 A 可以通过正交相似变换成对角型，即存在正交矩阵 Q,感谢阅读 使得 QT AQ = diag(λ ,λ ,…,λ ) (3.1)感谢阅读 1 2 n 其中 λ (i=1,2,…,n)是 A 的特征值，Q 中各列为相应的特征向量。精品文档放心下载 i 2) 在正交相似变换下,矩阵元素的平方和不变。即设 A=(a ) ,Q 交矩阵,精品文档放心下载 ij n×n 记 B=QTAQ=(b ) , 则 ij n×n Jacobi 方法的基本思想是通过一次正交变换,将 A 中的一对非零的非对角 化成零并且使得非对角元素的平方和减小。反复进行上述过程,使变换后的矩 阵的非对角元素的平方和趋于零，从而使该矩阵近似为对角矩阵，得到全部特 征值和特征向量。谢谢阅读 1 矩阵的旋转变换 设 A 为 n 阶实对称矩阵，考虑矩阵 易见V (φ)是正交矩阵, 记 ij 注意到 B=V A 的第 i,j 行元素以及 的第 i,j 列元素为谢谢阅读 ij .- 可得 如果 a ≠0，取 φ 使得 则有 ij 对 A(1) 重复上述的过程，可得 A(2) ,这样继续下去, 得到一个矩阵序列 {A(k) }。 可以证明，虽然这种变换不一定能使矩阵中非对角元素零元素的个数单调增加， 但可以保证非对角元素的平方和递减，我们以 A 与 A(1) 为例进行讨论。感谢阅读 设 由式(3.4) 可得 这表明,在上述旋转变换下,非对角元素的平方和严格单调递减,因而由(3.2)可谢谢阅读 .- 知，对角元素的平方和单调增加。 2. Jacobi 方法 通过一系列旋转变换将 A 变成 A(k+1) ,求得 A 的全部特征值与特征向量 的方法称为 Jacobi 方法。计算过程如下谢谢阅读 1)令 k=0, A(k) =A 2) 求整数 i,j, 使得 3) 计算旋转矩阵 4) 计算 A(k+1) 5) 计算 6) 若 E(A(k+1))ε, 则 为特征值， .- QT = (V(0) V(1) …V(k+1))T 精品文档放心下载 的各列为相应的特 征 向量；否则，k+1=k 返回 2,重复上述过 程。 例 5 用 Jacobi 方 法求矩阵 的特征值和特征向 量。 一般地，Jacobi 法不能在有限步内将 A 化成对角阵，但有下面的定理。谢谢阅读 定理 3 设 A 为 n 阶使对称矩阵，对 A 用 Jacobi 法得到序列 {A(k) }, 其中 A(0) = A, 则感谢阅读 证明 由 Jacobi 法计算过程 故有 (3.5) 另一方面，有计算 A