1. 幂法求特征值
问题:求 A 的特征值。
1.1 做法
- 初始化:选择一个初始向量 x0(通常为随机向量)。
- 迭代:计算 xk+1=∥Axk∥Axk,其中 A 是给定矩阵,∥⋅∥ 是某种范数(通常为欧几里得范数)。
- 收敛:当 xk+1 与 xk 足够接近时,停止迭代。此时 xk+1 近似为矩阵 A 的主特征向量,而主特征值可以通过 λ≈xk+1Txk+1xk+1TAxk+1 计算得到。
1.2 幂法的适用条件
- 矩阵 A 必须有主特征值,即存在一个特征值 λ1,使得 ∣λ1∣>∣λ2∣≥∣λ3∣≥⋯≥∣λn∣,其中 λ2,λ3,…,λn 是其他特征值。
- 矩阵 A 必须是方阵,即 A∈Rn×n。
- 初始向量 x0 不能与主特征向量正交,否则迭代过程中无法收敛到主特征向量。
1.3 收敛速度
幂法的收敛速度与主特征值 λ1 和次主特征值 λ2 之间的差距有关。具体来说,收敛速度可以用以下公式表示:
k→∞lim∥xk−xk−1∥∥xk+1−xk∥=∣λ1∣∣λ2∣.