第八章 矩阵特征值问题的求解.pdf

华北电力大学数理学院 School of mathematics physics §8.1 矩阵特征值问题 矩阵特征值问题： 设 A ∈R n×n ，求满足方程 Ax λx 或 (A −λI )x 0 的解(λ, x ). 应用背景丰富！例如系统稳定性分析。 MATLAB命令：D=eig(A)和[X,D]=eig(A) 华北电力大学数理学院 School of mathematics physics §8.2 矩阵特征值问题的幂法 幂法是求矩阵某指定特征值的一种简单方法！ 设矩阵A 的特征值： ，并有完全的 λ λ ≥≥ λ 1 2 n 特征向量系 x , x , , x . 则对任意非零向量v(0) ∈Rn 有 1 2 n (0) v αx =+αx ++α x 1 1 2 2 n n 华北电力大学数理学院 School of mathematics physics §8.2.1 求A 的按模最大特征值 对于 k 1, 2, 3, v(k) Av(k−1) k (0) A v k k k αλ x =+αλ x + +αλ x 1 1 1 2 2 2 n n n k ⎡ ⎛ ⎞λ k ⎛ ⎞λ k ⎤ ⎢ 2 n ⎥ λ αx =+⎜ ⎟ ⋅αx ++⎜ ⎟ ⋅α x 1 1 1 2 2 n n ⎢ λ1 λ1 ⎥ ⎣ ⎝ ⎠ ⎝ ⎠ ⎦ 当k充分大时，有 (k ) k k −1 (k −1) v ≈λ αx =λ λ α x ≈λv 1 1 1 1 ( 1 1 1 ) 1 Av(k ) ≈λv(k ) 1 华北电力大学数理学院 School of mathematics physics 矩阵按模最大特征值的幂法： k 1,2, (v(k-1) )T v(k) λ ≈lmta(k) 1 (v(k-1) )T v(k-1) (k )