譜分解(SD)


前提:矩陣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}$

 

參考

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM