文章目录
- 一、相似
- 二、若尔当型
- 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=λi1⋅⋅⋅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∑nciλi定理称:p(A)=i=0∑nciAi=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∑nciAi=i=0∑nciPJiP−1=P(i=0∑nciJi)P−1=Pp(J)P−1=det(A−λI)=det(J−λI)(相似矩阵具有相同的特征方程):p(λ)=(λ−λ1)a1(λ−λ2)a2⋯(λ−λk)ak⋯,αk为代数重数:p(Jk)=(Jk−λ1I)a1⋯(Jk−λkI)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−λkI)。Jk−λkI=0010⋱⋱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