矩阵对角化

**可對角化矩陣**是可化簡為對角矩陣的方陣。可對角化大幅化簡了矩陣的計算難度，因為對角矩陣在線性代數中特別容易處理，比如其行列式就是對角線上所有數字的乘積，而那些對角線上的數字正好就是其特徵值。

可對角化也使該線性變換的幾何意義更直觀，因為每個線性變換都可以對應到一個矩陣，所以將矩陣對角化等價於找到一組基底，使的線性變換的作用僅僅是伸縮基底向量而已。類似的，若用對角矩陣表示差分方程組或者微分方程組的係數的話，這樣每條等式只含有一個未知函數，這樣也大幅度了化簡了方程式的難度。

若爾當 - 謝瓦萊分解表達一個算子為它的對角部分與它的冪零部分的和。

## 正式定義

### 可對角化的線性映射

### 特徵化

關於可對角化映射和矩陣的基本事實可表達為如下:

* 在域 _F_ 上的 _n_ × _n_ 矩陣 _A_ 是可對角化的，若且唯若它的特徵空間的和的維度等於 _n_，它為真若且唯若存在由 _A_ 的特徵向量組成的 _F_n 的基。如果找到了這樣的基，可以形成有基向量作為縱列的矩陣 _P_，而 _P_-1_AP_ 將是對角矩陣。這個矩陣的對角元素是 _A_ 的特徵值。
* 線性映射 _T_ : _V_ → _V_ 是可對角化的，若且唯若它的特徵空間的維度等於 dim (_V_)，它為真若且唯若存在由 _T_ 的特徵向量組成的 _V_ 的基。_T_ 關於這個基將表示為對角矩陣。這個矩陣的對角元素是 _T_ 的特徵值。

另一個特徵化：矩陣或線性映射在域 _F_ 上可對角化的，若且唯若它的極小多項式在 _F_ 上有不同的線性因子。

下列充分 (但非必要) 條件經常是有用的。

* _n_ × _n_ 矩陣 _A_ 只在域 _F_ 上可對角化的，如果它在 _F_ 中有 _n_ 個不同的特徵值，就是說，如果它的特徵多項式在 _F_ 中有 _n_ 個不同的根。
* 線性映射 _T_ : _V_ → _V_ 帶有 _n_=dim(_V_) 是可對角化的，如果它有 _n_ 個不同的特徵值，就是說它的特徵多項式在 _F_ 中有 _n_ 個不同的根。

作為經驗規則，在複數域 **C** 上幾乎所有矩陣都是可對角化的。更精確地說：在 **C** 上不可對角化的複數 _n_ × _n_ 矩陣的集合被當作 **C**n×_n_ 的子集，它是關於勒貝格測度的零集。也可以說可對角化矩陣形成了關於 扎里斯基拓撲的稠密子集 : 補位於特徵多項式的判別式變為零的集合內，後者是超平面。從中得出的還有在平常的 (強拓撲) 中密度由範數給出。

對於 **R** 域就不是這樣了。隨著 _n_ 增長，隨機選擇的實數矩陣是在 **R** 上可對角化的可能性越來越小。

## 例子

### 可對角化矩陣

* 對合在實數上 (甚至特徵不是 2 的任何域) 是可對角化的，帶有 1 和 -1 在對角線上。
* 有限階自同態 (包括對合) 是在複數，或域的特徵不整除自同態的階的任何代數閉合域 (因為單位一的根是不同的) 是可對角化的，帶有單位根在對角線上。這是循環群的表示理論的一部分。
* 投影是可對角化的，帶有 0 和 1 在對角線上。

### 非可對角化的矩陣

某些矩陣在任何域上都是不可對角化的，最著名的是冪零矩陣。如果特徵值的幾何重次和代數重次不一致，這會更一般的出現。例如考慮

　
: $C={\begin{bmatrix}0&1\\0&0\end{bmatrix}}$

這個矩陣是不可對角化的：沒有矩陣 _U_ 使得 $U^{-1}CU$ 是對角矩陣。實際上，_C_ 有一個特徵值 (就是零) 而這個特徵值有代數重次 2 和幾何重次 1。

某些實數矩陣在實數上是不可對角化的。例如考慮

　
: $B={\begin{bmatrix}0&1\\-1&0\end{bmatrix}}$

矩陣 _B_ 沒有任何實數特徵值，所以沒有實數矩陣 _Q_ 使得 $Q^{-1}BQ$ 是對角矩陣。但是如果允許複數的話 ， $B$ 仍可以對角化。實際上，如果我們取

　
: $Q={\begin{bmatrix}1&{\textrm {i}}\\{\textrm {i}}&1\end{bmatrix}}$

則 $Q^{-1}BQ$ 是對角的。

### 矩陣對角化的方法

考慮矩陣

　
: $A={\begin{bmatrix}1&2&0\\0&3&0\\2&-4&2\end{bmatrix}}$

這個矩陣有特徵值

　
: $\lambda _{1}=3,\quad \lambda _{2}=2,\quad \lambda _{3}=1$

所以 _A_ 是有三個不同特徵值的 3 × 3 矩陣，所以它是可對角化的。

如果我們要對角化 _A_，我們需要計算對應的特徵向量。它們是

　
: $v_{1}={\begin{bmatrix}-1\\-1\\2\end{bmatrix}}\quad v_{2}={\begin{bmatrix}0\\0\\1\end{bmatrix}}\quad v_{3}={\begin{bmatrix}-1\\0\\2\end{bmatrix}}$

我們可以輕易的驗證 $Av_{k}=\lambda _{k}v_{k}$ 。

現在，設 _P_ 是由這些特徵向量作為縱列的矩陣:

　
: $P={\begin{bmatrix}-1&0&-1\\-1&0&0\\2&1&2\end{bmatrix}}$

則 _P_ 對角化了 _A_，簡單的計算可驗證:

　
: $P^{-1}AP={\begin{bmatrix}0&-1&0\\2&0&1\\-1&1&0\end{bmatrix}}{\begin{bmatrix}1&2&0\\0&3&0\\2&-4&2\end{bmatrix}}{\begin{bmatrix}-1&0&-1\\-1&0&0\\2&1&2\end{bmatrix}}={\begin{bmatrix}3&0&0\\0&2&0\\0&0&1\end{bmatrix}}$

注意特徵值 $\lambda _{k}$ 出現在對角矩陣中。

## 應用

對角化可被用來有效的計算矩陣 _A_ 的冪，假如矩陣是可對角化的。比如我們找到了

　
: $P^{-1}AP=D\ $

是對角矩陣，因為矩陣的積是結合的，

　
: ${\begin{aligned}A^{k}&=(PDP^{-1})^{k}=(PDP^{-1})\cdot (PDP^{-1})\cdots (PDP^{-1})\\&=PD(P^{-1}P)D(P^{-1}P)\cdots (P^{-1}P)DP^{-1}=PD^{k}P^{-1}\end{aligned}}$

而後者容易計算，因為它只涉及對角矩陣的冪。

在找到線性遞歸序列比如斐波那契數列的項的閉合形式的表達中這是非常有用的。

### 特定應用

例如，考慮下列矩陣:

　
: $M={\begin{bmatrix}a&b-a\\0&b\end{bmatrix}}$

計算 _M_ 個各次冪揭示了一個驚人的模式:

　
: $M^{2}={\begin{bmatrix}a^{2}&b^{2}-a^{2}\\0&b^{2}\end{bmatrix}},\quad M^{3}={\begin{bmatrix}a^{3}&b^{3}-a^{3}\\0&b^{3}\end{bmatrix}},\quad M^{4}={\begin{bmatrix}a^{4}&b^{4}-a^{4}\\0&b^{4}\end{bmatrix}},\quad \ldots $

上面的現象可以通過對角化 _M_ 來解釋。要如此我們需要由 _M_ 的特徵向量組成的 **R**2 的基。一個這樣的特徵向量基給出自

　
: $\mathbf {u} ={\begin{bmatrix}1\\0\end{bmatrix}}=\mathbf {e} _{1},\quad \mathbf {v} ={\begin{bmatrix}1\\1\end{bmatrix}}=\mathbf {e} _{1}+\mathbf {e} _{2}$

這裡的 **e**i 指示 **R**n 的標準基。
逆的基變更給出自

　
: $\mathbf {e} _{1}=\mathbf {u} ,\qquad \mathbf {e} _{2}=\mathbf {v} -\mathbf {u} $

直接計算證實

　
: $M\mathbf {u} =a\mathbf {u} ,\qquad M\mathbf {v} =b\mathbf {v} $

所以，_a_ 和 _b_ 是分別是對應於 **u** 和 **v** 的特徵值。
根據矩陣乘法的線性，我們有

　
: $M^{n}\mathbf {u} =a^{n}\,\mathbf {u} ,\qquad M^{n}\mathbf {v} =b^{n}\,\mathbf {v} $

切換回標準基，我們有

　
: $M^{n}\mathbf {e} _{1}=M^{n}\mathbf {u} =a^{n}\mathbf {e} _{1}$

　
: $M^{n}\mathbf {e} _{2}=M^{n}(\mathbf {v} -\mathbf {u} )=b^{n}\mathbf {v} -a^{n}\mathbf {u} =(b^{n}-a^{n})\mathbf {e} _{1}+b^{n}\mathbf {e} _{2}$

前面的關係用矩陣形式表達為

　
: $M^{n}={\begin{bmatrix}a^{n}&b^{n}-a^{n}\\0&b^{n}\end{bmatrix}}$

因此解釋了上述現象。

## 參見

* 若爾當標準型
* 縮放
* 三角矩陣

## 外部連結

* Diagonalization. PlanetMath.

## 引用

* Roger A. Horn and Charles R. Johnson, _Matrix Analysis_, Chapter 1, Cambridge University Press, 1985. ISBN 0-521-30586-1 (hardback), ISBN 0-521-38632-2 (paperback).