对角化是把矩阵通过相似变换化为对角矩阵的过程,可极大简化矩阵幂、矩阵指数等运算。
定义
矩阵A 通过相似变换矩阵P转化为对角矩阵D 的过程
P−1AP=D
条件
充分的线性无关特征向量:
- 对于一个 n×n 矩阵, 需要有 n 个线性无关 的特征向量.
代数重数=几何重数:
代数重数
几何重数
- 对每一个特征值,其几何重数必须等于代数重数;等价地,各特征空间维数之和为 n。特征向量不能通过“重复计算”增加数量,必须选出 n 个线性无关的特征向量。
步骤
设 A 为 n×n 矩阵:
- 计算特征值:解特征方程
A−λI=0
得到特征值 λ1,λ2,...,λn。 - 计算特征向量:对于每个 λi,求解
(A−λiI)xi=0
得到特征向量 xi。 - 构造相似变换矩阵
P=[x1,x2,...,xn]
- 相似变换得到对角矩阵
P−1AP=D
其中 D 为对角矩阵,对角线元素为 λ1,λ2,...,λn。
幂次
如果一个矩阵A可以被对角化,即存在可逆矩阵P和对角矩阵D,使得
A=PDP−1
那么A的幂次计算可以大大简化:对于任意正整数n,有
An=(PDP−1)n=PDnP−1
其中Dn是对角矩阵D的n次幂。如果D是一个对角矩阵,其形式为
D=λ10⋮00λ2⋮0⋯⋯⋱⋯00⋮λn
则Dn也是一个对角矩阵,其对角元素是λin,即
Dn=λ1n0⋮00λ2n⋮0⋯⋯⋱⋯00⋮λnn
因此,
An=PDnP−1=Pλ1n0⋮00λ2n⋮0⋯⋯⋱⋯00⋮λnnP−1
若 A 不可对角化,则无法走 PDnP−1 这条路,改用不可对角化矩阵的幂(凯莱-哈密顿定理带余除法或幂零矩阵拆分)。