线性代数(七)：Jordan标准形

Jordan标准形

定义7.1：称形如$$J(\lambda )={\left[\begin{array}{cccccc}\lambda & 1& & & & \\ & \lambda & 1& & & \\ & & \lambda & 1& & \\ & & & \ddots & \ddots & \\ & & & & \lambda & 1\\ & & & & & \lambda \end{array}\right]}_{r\times r}$$
的$r$阶方阵为为一个$r$阶 Jordan块。由Jordan块构成的准对角矩阵$$J=\left[\begin{array}{cccccc}{J}_1({\lambda }_1)& & & & & \\ & J_2({\lambda }_2)& & & & \\ & & J_3({\lambda }_3)& & & \\ & & & \ddots & & \\ & & & & J_{m-1}({\lambda }_{m-1})& \\ & & & & & J_m({\lambda }_m)\end{array}\right]$$则称为 Jordan矩阵。

Remark:
Jordan矩阵不仅是准对角矩阵，也是上三角矩阵，主对角元素为自身的特征值。Jordan矩阵非零元素只存在于主对角线及其上方元素 $a_{i,i+1}$,并且$a_{i,i+1}=0或1$。另外注意，不同Jordan块并不要求特征值不同。

我们知道并不是所有的$n$阶矩阵均可相似对角化，当且仅当其有$n$个线性无关的特征向量才可实现。但是我们可以得到如下结论：
定理7.1：在复数域上，任意$n$阶方阵$A$均相似于一个Jordan矩阵$J$，即存在$n$阶可逆矩阵$P$，使得：$$P^{-1}AP=J=\left[\begin{array}{cccccc}{J}_1({\lambda }_1)& & & & & \\ & J_2({\lambda }_2)& & & & \\ & & J_3({\lambda }_3)& & & \\ & & & \ddots & & \\ & & & & J_{m-1}({\lambda }_{m-1})& \\ & & & & & J_m({\lambda }_m)\end{array}\right]$$称$J$为$A$的Jordan标准形。

下面分析Jordan标准型的求解方法：

设$A_{n\times n}\sim J_A$，则 $J_A$ 与 $A$ 有相同的特征值，即 $A$ 的特征值为 $J_A$ 的主对角元素。

记 $$A的k_i重特征值为{\lambda }_i(i=1,2,\dots ,s;\sum _{i=1}^s{k}_i=n)$$那么将 $J_A$ 具有相同特征值的Jordan块排在一块有：$$J_A={\left[\begin{array}{cccc}{A}_1({\lambda }_1{)}_{k_1\times k_1}& & & \\ & A_2({\lambda }_2{)}_{k_2\times k_2}& & \\ & & \ddots & \\ & & & A_s({\lambda }_s{)}_{k_s\times k_s}\end{array}\right]}_{n\times n}其中,A_i({\lambda }_i)包含1\le t_i\le k_i个Jordan块,各Jordan块的阶数为n_{(i,j)}，则(\sum _{j=1}^{t_i}{n}_{(i,j)}=k_i)，即A_i({\lambda }_i)={\left[\begin{array}{cccc}{J}_1({\lambda }_i{)}_{n_{(i,1)}\times n_{(i,1)}}& & & \\ & J_2({\lambda }_i{)}_{n_{(i,2)}\times n_{(i,2)}}& & \\ & & \ddots & \\ & & & J_{t_i}({\lambda }_i{)}_{n_{(i,t_i)}\times n_{(i,t_i)}}\end{array}\right]}_{k_i\times k_i}(i=1,2,\dots ,s)$$
由$A_{n\times n}\sim J_A$ ，则$$P^{-1}AP=J_A⟹AP=P{J}_A⟹A[p_{n\times k_1}^1,p_{n\times k_2}^2,\dots ,p_{n\times k_s}^s]=[p_{n\times k_1}^1,p_{n\times k_2}^2,\dots ,p_{n\times k_s}^s]{\left[\begin{array}{cccc}{A}_1({\lambda }_1{)}_{k_1\times k_1}& & & \\ & A_2({\lambda }_2{)}_{k_2\times k_2}& & \\ & & \ddots & \\ & & & A_s({\lambda }_s{)}_{k_s\times k_s}\end{array}\right]}_{n\times n}=[p_{n\times k_1}^1{A}_1({\lambda }_1{)}_{k_1\times k_1},p_{n\times k_2}^2{A}_2({\lambda }_2{)}_{k_2\times k_2},\dots ,p_{n\times k_s}^s{A}_s({\lambda }_s{)}_{k_s\times k_s}]$$即：$$A{p}_{n\times k_r}^r=p_{n\times k_r}^r{A}_r({\lambda }_r{)}_{k_r\times k_r}(r=1,2,\dots ,s)$$
进一步，对$p_{n\times k_r}^r$进行列分块：$$A[{\tilde{P}}_{n\times n_{(r,1)}}^{r1},{\tilde{P}}_{n\times n_{(r,2)}}^{r2},\dots ,{\tilde{P}}_{n\times n_{(r,t_r)}}^{r{k}_r}]=[{\tilde{P}}_{n\times n_{(r,1)}}^{r1},{\tilde{P}}_{n\times n_{(r,2)}}^{r2},\dots ,{\tilde{P}}_{n\times n_{(r,t_r)}}^{r{k}_r}]{\left[\begin{array}{cccc}{J}_1({\lambda }_r{)}_{n_{(r,1)}\times n_{(r,1)}}& & & \\ & J_2({\lambda }_r{)}_{n_{(r,2)}\times n_{(r,2)}}& & \\ & & \ddots & \\ & & & J_{t_r}({\lambda }_r{)}_{n_{(r,t_r)}\times n_{(r,t_r})}\end{array}\right]}_{k_r\times k_r}=[{\tilde{P}}^{r1}{J}_1({\lambda }_r),{\tilde{P}}^{r2}{J}_2({\lambda }_r),\dots ,{\tilde{P}}^{r{k}_r}{J}_{t_r}({\lambda }_r)]$$即：$$A{\tilde{P}}_{n\times n_{(r,l)}}^{rl}={\tilde{P}}_{n\times n_{(r,l)}}^{rl}{J}_l({\lambda }_r{)}_{n_{(r,l)\times (r,l)}}(l=1,2,\dots ,t_r)$$
进一步分块有：$$A[\vec{\alpha },{\vec{\beta }}_2,\dots ,{\vec{\beta }}_{n_{(r,tr)}}]=[\vec{\alpha },{\vec{\beta }}_2,\dots ,{\vec{\beta }}_{n_{(r,tr)}}]{\left[\begin{array}{cccccc}\lambda & 1& & & & \\ & \lambda & 1& & & \\ & & \lambda & 1& & \\ & & & \ddots & \ddots & \\ & & & & \lambda & 1\\ & & & & & \lambda \end{array}\right]}_{n_{(r,l)\times (r,l)}}$$即：$$\{\begin{array}{l}(A-\lambda E)\vec{\alpha }=0\\ (A-\lambda E){\vec{\beta }}_2=\vec{\alpha }\\ (A-\lambda E){\vec{\beta }}_3={\vec{\beta }}_2\\ \dots \\ (A-\lambda E){\vec{\beta }}_{n_{(r,t_r)}}={\vec{\beta }}_{(n_{(r,t_r)}-1)}\end{array}(1)$$
将 { α ⃗ , β ⃗ 2 , … , β ⃗ n ( r , t r ) } \{\vec{\alpha},\vec{\beta}_2,\dots,\vec{\beta}_{n_{(r,tr)}}\} {α ,β ​2​,…,β ​n(r,tr)​​}称作Jotdan链，它们是同一特征值对应的某个Jordan块对应的可逆矩阵$P$的列

综上，求解可逆矩阵$P$与Jordan标准形$J_A$的步骤如下：
(1) 求解$A$的特征多项式，得到$k_i$重特征值${\lambda }_i(i=1,2,\dots ,s)$；
(2) 求解其中一个特征值${\lambda }_i$对应的所有线性无关的特征向量${\vec{\alpha }}_j(j=1,2,\dots ,t_i)$;
(3) 将每个特征向量${\vec{\alpha }}_j$分别代入方程$(1)$得到一个Jordan链；
(4) 对应排列Jordan链与相应的Jordan块得到P与Jordan标准形。