文章目錄
- 一、相似
- 二、若爾當型
- 1.1 認識若爾當型
- 1.2 凱萊-哈密頓定理 (Cayley-Hamilton Theorem)
一、相似
Every matrix C=B?1ABC = B^{-1}ABC=B?1AB has the same eigenvalues as A. These C’s are similar to A.
任意一個矩陣C,滿足 C=B?1ABC = B^{-1}ABC=B?1AB 都和A有著相同的特征值。C 和 A 是相似的。
相似具有:
- 自反性
- 對稱性
- 傳遞性
All the matrices A=B?1CBA = B^{-1}CBA=B?1CB are “similar”,They all share the eigenvalues of C.
所有滿足 A=B?1CBA = B^{-1}CBA=B?1CB 的矩陣都是相似的,他們的特征值相同。
證明:
如果Cx=λx,則BCB?1有特征值同樣為λ的特征向量Bx事實上,(BCB?1)(Bx)=λ(Bx)同理可證A的特征值都是C的特征值\begin{align} & 如果 Cx = \lambda x, 則 BCB^{-1} 有特征值同樣為\lambda 的特征向量 Bx \\ & 事實上,(BCB^{-1})(Bx) = \lambda (Bx) \\ & 同理可證 A 的特征值都是C的特征值 \end{align} ?如果Cx=λx,則BCB?1有特征值同樣為λ的特征向量Bx事實上,(BCB?1)(Bx)=λ(Bx)同理可證A的特征值都是C的特征值??
如果矩陣A 和 B 均可相似對角化,那么我們可以輕松判別二者是否相似。
那么不可對角化的情況呢?
二、若爾當型
1.1 認識若爾當型
(Jordan form) If A has s independent eigenvectors, it is similar to a matrix J that has s Jordan blocks J1,...,JsJ_1, ..., J_sJ1?,...,Js? on its diagonal. B contains “generalized eigenvectors”:
若爾當型(Jordan form):如果A 有s 個線性無關的特征向量,那么 A 相似于矩陣J,J的對角線上有 s 個 若爾當塊J1,...,JsJ_1, ..., J_sJ1?,...,Js?,矩陣B包含了**“廣義特征向量”**。
廣義特征向量定義如下:
一個向量 vvv 是對應于特征值 λ\lambdaλ 的廣義特征向量,如果它滿足: (A?λI)kv=0(A - \lambda I)^k v = 0(A?λI)kv=0 對于某個正整數 k≥1k \ge 1k≥1 成立。
當 k>1k>1k>1 時,我們得到了一系列向量,它們構成一個若爾當鏈 (Jordan Chain)。
一個長度為 kkk 的若爾當鏈由向量 v_1,v_2,…,v_k{v\_1, v\_2, \dots, v\_k}v_1,v_2,…,v_k 構成,它們滿足: (A?λI)vk=vk?1(A - \lambda I)v_k = v_{k-1}(A?λI)vk?=vk?1?
- 所有的特征向量都是廣義特征向量,但并非所有廣義特征向量都是特征向量。
- k 是使該零空間維度達到最大并穩定的最小整數(也就是 λ 對應的最大若爾當塊的階數)。
- 特征空間完全包含在廣義特征空間之中。
- 如何求解vkv_kvk? =>只需找到$k => v_k = N(A - \lambda I)^{k} \setminus N(A - \lambda I)^{k - 1} $
A 的若爾當型:
Jordan?form?of?AB?1AB=[J1?Js]=J\text { Jordan form of } A \quad B^{-1} A B=\left[\begin{array}{lll} J_{1} & & \\ & \ddots & \\ & & J_{s} \end{array}\right]=J ?Jordan?form?of?AB?1AB=?J1????Js???=J
每一個若爾當塊 JiJ_iJi? 主對角線上都是特征值 λi\lambda_iλi?,主對角線上方恰有1個1.
Jordan?block?in?JJi=[λi1???1λi].?\text { Jordan block in } J \quad J_{i}=\left[\begin{array}{cccc} \lambda_{i} & 1 & & \\ & \cdot & \cdot & \\ & & \cdot & 1 \\ & & & \lambda_{i} \end{array}\right] \text {. } ?Jordan?block?in?JJi?=?λi??1?????1λi???.?
一個 k×kk \times kk×k 的若爾當塊 J_k(λ)J\_k(\lambda)J_k(λ) 就對應一個長度為 kkk 的若爾當鏈。將這些若爾當鏈的所有向量集合在一起,就構成了若爾當基 (Jordan Basis),也就是相似變換矩陣 BBB 的列向量。
一些性質:
- 每一個矩陣都相似于一個若爾當型。
- 若爾當塊的數目等于線性無關特征向量的個數(幾何重數)。
如何求解若爾當型?
無需求解相似變換矩陣B
只需要知道每個特征值對應的若爾當塊的數量和大小。
步驟如下:
- 求解所有特征值,以及其代數重數(代數重數是指特征值作為特征方程根的重數)。
- 對每一個特征值 λi\lambda_iλi? 單獨分析,假設特征值 λ\lambdaλ 的代數重數是 mmm。這意味著與 λ\lambdaλ 相關的所有若爾當塊的階數之和必須為 mmm,關鍵在于確定塊的數量和各自的大小:
1. 若爾當塊的數量: 對應于特征值 λ\lambdaλ 的若爾當塊的數量,等于該特征值的幾何重數 (Geometric Multiplicity)(即線性無關特征向量的個數)。
2. k階若爾當塊的數量:rank((A?λI)k?1)+rank((A?λI)k?1)?2rank((A?λI)k)rank((A - \lambda I)^{k - 1}) + rank((A - \lambda I)^{k - 1}) - 2rank((A - \lambda I)^{k})rank((A?λI)k?1)+rank((A?λI)k?1)?2rank((A?λI)k)
至此,我們了解了若爾當型,并且得出 A 和 B 相似的充要條件是 A和B 有著相同的若爾當標準型。
1.2 凱萊-哈密頓定理 (Cayley-Hamilton Theorem)
凱萊-哈密頓定理 (Cayley-Hamilton Theorem):每一個n階方陣 A,都滿足其自身的特征方程。
記p(λ)=det(A?λI)展開:p(λ)=∑i=0nciλi定理稱:p(A)=∑i=0nciAi=0\begin{align} & 記 p(\lambda) = det(A - \lambda I) \\ & 展開:p(\lambda) = \sum_{i=0}^{n}c_i \lambda^i \\ & 定理稱:p(A) = \sum_{i=0}^{n}c_i A^i = 0 \end{align} ?記p(λ)=det(A?λI)展開:p(λ)=i=0∑n?ci?λi定理稱:p(A)=i=0∑n?ci?Ai=0??
證明:
設A=PJP?1(若爾當型)Ak=PJkP?1p(A)=∑i=0nciAi=∑i=0nciPJiP?1=P(∑i=0nciJi)P?1=Pp(J)P?1由于,p(λ)=det(A?λI)=det(J?λI)(相似矩陣具有相同的特征方程)可寫成:p(λ)=(λ?λ1)a1(λ?λ2)a2?(λ?λk)ak?,αk為代數重數則有:p(Jk)=(Jk?λ1I)a1?(Jk?λkI)ak?\begin{align} 設 A &= PJP^{-1}(若爾當型) \\ A^k &= PJ^kP^{-1} \\ p(A) &= \sum_{i=0}^{n}c_iA^i \\ &= \sum_{i=0}^{n}c_iPJ^{i}P^{-1} \\ &= P(\sum_{i=0}^{n} c_iJ^i)P^{-1} \\ &= Pp(J)P^{-1} \\ 由于,p(\lambda)&= det(A - \lambda I) = det(J - \lambda I)(相似矩陣具有相同的特征方程) \\ 可寫成&:p(\lambda) = (\lambda - \lambda_1)^{a_1} (\lambda - \lambda_2)^{a_2} \cdots (\lambda - \lambda_k)^{a_k} \cdots, \alpha_k 為代數重數 \\ 則有&: p(J_k) = (J_k - \lambda_1 I)^{a_1} \cdots (J_k - \lambda_k I)^{a_k} \cdots \\ \end{align} 設AAkp(A)由于,p(λ)可寫成則有?=PJP?1(若爾當型)=PJkP?1=i=0∑n?ci?Ai=i=0∑n?ci?PJiP?1=P(i=0∑n?ci?Ji)P?1=Pp(J)P?1=det(A?λI)=det(J?λI)(相似矩陣具有相同的特征方程):p(λ)=(λ?λ1?)a1?(λ?λ2?)a2??(λ?λk?)ak??,αk?為代數重數:p(Jk?)=(Jk??λ1?I)a1??(Jk??λk?I)ak????
我們來仔細分析其中的一個關鍵因子:(Jk?λkI)。Jk?λkI=(0100??100)這是一個冪零矩陣Nm×m,滿足Nm=0而m≤代數重數=α故p(Jk)=0,定理得證\begin{align} &我們來仔細分析其中的一個關鍵因子:(J_k - \lambda_k I)。\\ &J_k - \lambda_k I = \begin{pmatrix} 0 & 1 & & \mathbf{0} \\ & 0 & \ddots & \\ & & \ddots & 1 \\ \mathbf{0} & & & 0 \end{pmatrix}\\ & 這是一個冪零矩陣N_{m\times m}, 滿足 N^{m} = 0 \\ & 而 m \le 代數重數 = \alpha \\ & 故p(J_k) = 0, 定理得證 \end{align} ?我們來仔細分析其中的一個關鍵因子:(Jk??λk?I)。Jk??λk?I=?00?10????010??這是一個冪零矩陣Nm×m?,滿足Nm=0而m≤代數重數=α故p(Jk?)=0,定理得證??
一個比較實用的應用:求解矩陣的逆
p(λ)=λ2?6λ+5A2?6A+5I=05I=6A?A25A?1=6I?AA?1=6I?A5\begin{align} & p(\lambda) = \lambda^2 -6\lambda + 5 \\ & A^2?6A+5I=0 \\ & 5I = 6A - A^2 \\ & 5A^{-1} = 6I - A \\ & A^{-1} = \frac{6I - A}{5} \end{align} ?p(λ)=λ2?6λ+5A2?6A+5I=05I=6A?A25A?1=6I?AA?1=56I?A???