前提:矩陣A必須可相似對角化!
充分條件:
- $A$ 是實對稱矩陣
- $A$ 有 $n$ 個互異特征值
- $A^{\wedge} 2=A $
- $\mathrm{A}^{\wedge} 2=\mathrm{E} $
- $ r(A)=1 且 \operatorname{tr}(A) !=0$
譜分解(Spectral Decomposition ),又稱特征分解,或相似標准形分解,是將矩陣分解為由其特征值和特征向量表示的矩陣之積的方法,需要注意只有對可對角化矩陣才可以施以特征分解。它體現了線性變換的旋轉和縮放的功效。
設 $A$ 為 $n$ 階實對稱陣,則必有正交陣 $P$ ,使 $A=P \Lambda P^{T} $ 其中 $ \Lambda $ 是以 $ A $ 的 $n$ 個特征值為對角元的對角陣, $P$ 是由 $A$ 的 $ n$ 個特征向量得到 的正交矩陣。
實對稱矩陣譜分解的步驟
設 $\boldsymbol{A} \in \boldsymbol{R}^{n \times n}$, $ \boldsymbol{A}^{\prime}=\boldsymbol{A} $
(i) 求出 $A$ 的所有不同的特征值: $\lambda_{1}, \lambda_{2}, \cdots, \lambda_{r} \in R $,其重數 $n_{1}, n_{2}, \cdots, n_{r} $ 必滿足 $\sum_{i=1}^{r} n_{i}=n $;
(ii) 對每個 $ \lambda_{i} $ ,解齊次線性方程組
$\left.\left(\lambda_{i}\right) E-A\right) X=0$
求出它的一個基礎解系: $\alpha_{i 1}, \alpha_{i 2}, \cdots, \alpha_{i n} $
把它們按 Schmidt 正交化過程化成 標准正交組
$\boldsymbol{\eta}_{i 1}, \boldsymbol{\eta}_{i 2}, \cdots, \boldsymbol{\eta}_{i n}$
(iii) 因為 $\lambda_{1}, \lambda_{2}, \cdots, \lambda_{r} $ 互不相同,所以
將 $\boldsymbol{\eta}_{11}, \boldsymbol{\eta}_{12}, \cdots, \boldsymbol{\eta}_{1 n_{1}}, \cdots, \boldsymbol{\eta}_{r 1}, \boldsymbol{\eta}_{r 2}, \cdots, \boldsymbol{\eta}_{r \boldsymbol{n}_{r}}$ 的分量依次作 矩陣 $P$的第 $1,2, \cdots, n$ 列, 則 $P$ 是正交矩陣, 且有 $\boldsymbol{A}=\boldsymbol{P} \boldsymbol{\Lambda} \boldsymbol{P}^{\boldsymbol{T}}$ .
例題
已知 $A=\left[\begin{array}{ccc}0 & -1 & 1 \\-1 & 0 & 1 \\1 & 1 & 0\end{array}\right]$ 求一個正交矩陣 $P$ , 使 $\boldsymbol{A}=\boldsymbol{P} \boldsymbol{\Lambda} \boldsymbol{P}^{\boldsymbol{T}}$ 為對角陣.
第1步: 求特征值.
$\begin{array}{c}|\lambda E-A|=\left|\begin{array}{ccc}\lambda & 1 & -1 \\1 & \lambda & -1 \\-1 & -1 & \lambda\end{array}\right| \underline{r_{1}-r_{2}}\left|\begin{array}{cccc}\lambda-1 & 1-\lambda & 0 \\1 & \lambda & -1 \\-1 & -1 & \lambda\end{array}\right| \\=(\lambda-1)\left|\begin{array}{ccc}1 & -1 & 0 \\1 & \lambda & -1 \\-1 & -1 & \lambda\end{array}\right|=(\lambda-1)\left|\begin{array}{ccc}1 & 0 & 0 \\1 & \lambda+1 & -1 \\-1 & -2 & \lambda\end{array}\right|=(\lambda+2)(\lambda-1)^{2} \\\lambda_{1}=-2, \lambda_{2}=\lambda_{3}=1\end{array}$
第2步: 求線性無關的特征向量.
對 $\lambda_{1}=-2$ , 解方程組 $(A+2 E) x=0$
$A+2 E=\left[\begin{array}{ccc}2 & -1 & 1 \\-1 & 2 & 1 \\1 & 1 & 2\end{array}\right] \rightarrow\left[\begin{array}{lll}1 & 0 & 1 \\0 & 1 & 1 \\0 & 0 & 0\end{array}\right]$
求得基礎解系(即最大無關特征向量)
$\alpha_{1}=\left[\begin{array}{c}-1 \\-1 \\1\end{array}\right]$
對 $\lambda_{2}=\lambda_{3}=1$ , 解方程組 $(A-E) x=0 $
$A-E=\left[\begin{array}{ccc}-1 & -1 & 1 \\-1 & -1 & 1 \\1 & 1 & -1\end{array}\right] \stackrel{r}{\rightarrow}\left[\begin{array}{ccc}1 & 1 & -1 \\0 & 0 & 0 \\0 & 0 & 0\end{array}\right]$
求得基礎解系(即最大無關特征向量)
$\begin{array}{l}\alpha_{2}=\left[\begin{array}{c}-1 \\1 \\0\end{array}\right], \alpha_{3}=\left[\begin{array}{c}1 \\0 \\1\end{array}\right] , \alpha_{1}=\left[\begin{array}{c}-1 \\-1 \\1\end{array}\right] \\{\left[\alpha_{1}, \alpha_{2}\right]=?\left[\alpha_{1}, \alpha_{3}\right]=? \quad\left[\alpha_{2}, \alpha_{3}\right]=0 ?}\end{array}$
第3步: 檢驗重特征值對應的特征向量是否正交, 如果不正交,
用施密特過程正交化; 再把正交的特征向量單位化.
$\begin{array}{l}\boldsymbol{\alpha}_{\mathbf{2}}=\left[\begin{array}{c}-\mathbf{1} \\\mathbf{1} \\\mathbf{0}\end{array}\right], \boldsymbol{\alpha}_{\mathbf{3}}=\left[\begin{array}{l}\mathbf{1} \\\mathbf{0} \\\mathbf{1}\end{array}\right] \\\boldsymbol{\beta}_{\mathbf{2}}=\boldsymbol{\alpha}_{\mathbf{2}} \\\boldsymbol{\beta}_{3}=\alpha_{3}-\frac{\left[\beta_{2}, \alpha_{3}\right]}{\left[\beta_{2}, \beta_{2}\right]} \boldsymbol{\beta}_{2}=\left[\begin{array}{l}1 \\0 \\1\end{array}\right]+\frac{1}{2}\left[\begin{array}{c}-1 \\1 \\0\end{array}\right]=\frac{1}{2}\left[\begin{array}{l}1 \\1 \\2\end{array}\right]\end{array}$
單位化:
$\xi_{1}=\frac{\alpha_{1}}{\left\|\alpha_{1}\right\|}=\frac{1}{\sqrt{3}}\left[\begin{array}{c}-1 \\-1 \\1\end{array}\right] \quad \xi_{2}=\frac{\beta_{2}}{\left\|\beta_{2}\right\|}=\frac{1}{\sqrt{2}}\left[\begin{array}{c}-1 \\1 \\0\end{array}\right] \quad \xi_{3}=\frac{\beta_{3}}{\left\|\beta_{3}\right\|}=\frac{1}{\sqrt{6}}\left[\begin{array}{l}1 \\1 \\2\end{array}\right]$
第4步: 把求得的規范正交特征向量拼成正交矩陣.
令 $ P=\left[\xi_{1}, \xi_{2}, \xi_{3}\right]=\left[\begin{array}{ccc}\frac{-1}{\sqrt{3}} & \frac{-1}{\sqrt{2}} & \frac{1}{\sqrt{6}} \\ \frac{-1}{\sqrt{3}} & \frac{1}{\sqrt{2}} & \frac{1}{\sqrt{6}} \\ \frac{1}{\sqrt{3}} & 0 & \frac{2}{\sqrt{6}}\end{array}\right] $
則 $A=P\left[\begin{array}{ccc}-2 & & \\ & 1 & \\ & & 1\end{array}\right] P^{T}$