主题
Search

特征分解


一个方阵 A矩阵分解,分解成所谓的特征值和特征向量,是非常重要的。这种分解通常被称为“矩阵对角化”。然而,这个名称并非最佳,因为所描述的过程实际上是将矩阵分解为其他三个矩阵的乘积,其中只有一个是对角矩阵,而且因为所有其他标准类型的矩阵分解都在其名称中使用术语“分解”,例如,Cholesky 分解Hessenberg 分解等等。因此,在这项工作中,将矩阵分解为由其特征向量和特征值组成的矩阵称为特征分解。

假设 A 具有非退化的特征值 lambda_1,lambda_2,...,lambda_k 和相应的线性独立的特征向量 X_1,X_2,...,X_k,可以表示为

 [x_(11); x_(12); |; x_(1k)],[x_(21); x_(22); |; x_(2k)],...[x_(k1); x_(k2); |; x_(kk)].
(1)

定义由特征向量组成的矩阵

P=[X_1 X_2 ... X_k]
(2)
=[x_(11) x_(21) ... x_(k1); x_(12) x_(22) ... x_(k2); | | ... |; x_(1k) x_(2k) ... x_(kk)]
(3)

和特征值

 D=[lambda_1 0 ... 0; 0 lambda_2 ... 0; | | ... |; 0 0 ... lambda_k],
(4)

其中 D 是一个对角矩阵。那么

AP=A[X_1 X_2 ... X_k]
(5)
=[AX_1 AX_2 ... AX_k]
(6)
=[lambda_1X_1 lambda_2X_2 ... lambda_kX_k]
(7)
=[lambda_1x_(11) lambda_2x_(21) ... lambda_kx_(k1); lambda_1x_(12) lambda_2x_(22) ... lambda_kx_(k2); | | ... |; lambda_1x_(1k) lambda_2x_(2k) ... lambda_kx_(kk)]
(8)
=[x_(11) x_(21) ... x_(k1); x_(12) x_(22) ... x_(k2); | | ... |; x_(1k) x_(2k) ... x_(kk)][lambda_1 0 ... 0; 0 lambda_2 ... 0; | | ... |; 0 0 ... lambda_k]
(9)
=PD,
(10)

给出 A 的惊人分解,分解成一个涉及 PD相似变换

 A=PDP^(-1).
(11)

只要 P 是一个方阵,对于一个方阵 A 来说,这种分解总是可能的,这在本文中被称为特征分解定理

此外,等式 (11) 两边平方得到

A^2=(PDP^(-1))(PDP^(-1))
(12)
=PD(P^(-1)P)DP^(-1)
(13)
=PD^2P^(-1).
(14)

通过归纳法,可以得出对于一般的正整数幂,

 A^n=PD^nP^(-1).
(15)

A 的逆矩阵是

A^(-1)=(PDP^(-1))^(-1)
(16)
=PD^(-1)P^(-1),
(17)

其中对角矩阵 D 的逆矩阵很容易由下式给出

 D^(-1)=[lambda_1^(-1) 0 ... 0; 0 lambda_2^(-1) ... 0; | | ... |; 0 0 ... lambda_k^(-1)].
(18)

因此,方程 (◇) 对于 n 以及正数都成立。

涉及矩阵 PD 的另一个显著结果来自矩阵指数的定义

e^(A)=sum_(n=0)^(infty)(A^n)/(n!)
(19)
=sum_(n=0)^(infty)(PD^nP^(-1))/(n!)
(20)
=P(sum_(n=0)^(infty)(D^n)/(n!))P^(-1)
(21)
=Pe^(D)P^(-1).
(22)

这是正确的,因为 D 是一个对角矩阵,并且

e^(D)=sum_(n=0)^(infty)(D^n)/(n!)
(23)
=sum_(n=0)^(infty)1/(n!)[lambda_1^n 0 ... 0; 0 lambda_2^n ... 0; | | ... |; 0 0 ... lambda_k^n]
(24)
=[sum_(n=0)^(infty)(lambda_1^n)/(n!) 0 ... 0; 0 sum_(n=0)^(infty)(lambda_2^n)/(n!) ... 0; | | ... |; 0 0 ... sum_(n=0)^(infty)(lambda_k^n)/(n!)]
(25)
=[e^(lambda_1) 0 ... 0; 0 e^(lambda_2) ... 0; | | ... |; 0 0 ... e^(lambda_k)],
(26)

因此,可以使用 D 找到 e^(A)


另请参阅

可对角化矩阵, 特征分解定理, 特征值, 特征向量, 矩阵分解, 奇异值分解

使用 Wolfram|Alpha 探索

请引用为

Weisstein, Eric W. “特征分解。” 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/EigenDecomposition.html

学科分类