3.2_Jordan标准形:若干论断和应用

3.2 Jordan标准形:若干论断和应用

3.2.1 Jordan 矩阵的结构 Jordan 矩阵

J=[Jn1(λ1)0Jn2(λ2)0Jnk(λk)],n1n2nk=n(3.2.1.1)J = \left[ \begin{array}{c c c} J _ {n _ {1}} (\lambda_ {1}) & & 0 \\ & J _ {n _ {2}} (\lambda_ {2}) & \\ & & \ddots \\ 0 & & J _ {n _ {k}} (\lambda_ {k}) \end{array} \right], \quad n _ {1} \mid n _ {2} \mid \dots \mid n _ {k} = n \tag {3.2.1.1}

有一个确定的结构,它使得该矩阵以及与其相似的任一矩阵的某些基本性质变得更加明显。

  1. Jordan 块的个数 kk (相同的子块计重复出现的次数) 是 JJ 的线性无关特征值向量的个数.

  2. 矩阵 JJ 可对角化,当且仅当 k=nk = n .

  3. 相应于一个已知特征值的 Jordan 块的个数是该特征值的几何重数。它是相应的特征空间的维数。相应于一个已知特征值的所有 Jordan 块的阶数之和是该特征值的代数重数。

  4. 知道了诸特征值以及它们的代数重数和几何重数,一般不能完全确定一个 Jordan 矩阵。我们还需知道相应于每个特征值的各 Jordan 块的阶数,相应于一个特征值 λ\lambda 的最大 Jordan 块的阶数是 λ\lambda 作为其极小多项式的一个根的重数(3.3.6)。

  5. 相应于某个特征值的诸 Jordan 块的阶数可以通过计算某些幂的秩来确定。例如,如果

j=[21002100022102210022].j = \left[ \begin{array}{c c c c c c c} 2 & 1 & 0 & & & & \\ 0 & 2 & 1 & & & 0 & \\ 0 & 0 & 2 & & & \\ & & & 2 & 1 & & \\ & & & 0 & 2 & & \\ & & & & & 2 & 1 \\ & & 0 & & & 0 & 2 \\ & & & & & & 2 \end{array} \right].

经计算

J2I=[01000100000100010000000000].J - 2 I = \left[ \begin{array}{c c c c c c c c} 0 & 1 & 0 & & & & \\ 0 & 0 & 1 & & & & \\ 0 & 0 & 0 & & & & 0 \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ & & 0 & 1 & & & \\ & & 0 & 0 & & & \\ & & \vdots & \vdots & \vdots & \vdots & \vdots \\ & & & 0 & 1 & \\ 0 & 0 & 0 & 0 & 0 & \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ 0 & 0 & 0 & 0 & 0 \end{array} \right].
(J2I)2=[00100000000000000000](J - 2 I) ^ {2} = \left[ \begin{array}{c c c c c c c c c} 0 & 0 & 1 & & & & & \\ 0 & 0 & 0 & & & & & \\ 0 & 0 & 0 & \vdots & & & 0 & \\ \dots & \dots & \dots & \dots & \dots & \vdots & \\ & & 0 & 0 & \vdots & & \\ & & 0 & 0 & \vdots & \\ & & & \vdots & 0 & 0 \\ & 0 & & \vdots & 0 & 0 \\ & & \vdots & \vdots & \vdots & \vdots & 0 \end{array} \right]

(J2I)1=0(J - 2I)^{-1} = 0 ,因此,我们知道

(J2I)4=0,(J - 2 I) ^ {4} = 0,
rank(J2I)2=1,\operatorname {r a n k} (J - 2 I) ^ {2} = 1,
rank(J2I)=4,\operatorname {r a n k} (J - 2 I) = 4,

(J2I)(J - 2I)8×88 \times 8 矩阵.

这组数足以确定 JJ 的分块结构。事实上, (J2I)3=0(J - 2I)^3 = 0 说明最大子块的阶是3。 (J2I)2(J - 2I)^2 的秩就是3阶子块的个数,因而只有一个3阶子块。 (J2I)(J - 2I) 的秩是3阶子块的个数的2倍加上2阶子块的个数,因而有两个2阶子块。1阶子块的个数是 8(2×2)3=18 - (2 \times 2) - 3 = 1 。可以把同样的手续应用于具有任意阶数的若干Jordan块的直和,只要它们是相应于同一个特征值的所有子块。如果 JJ 是相应子特征值 λ\lambda 的这样一个直和,那么使 (JλI)k1=0(J - \lambda I)^{k_1} = 0 的最小整数 k1k_1 便是最大子块的阶数。 (JλI)k1(J - \lambda I)^{k_1} 的秩是 k1k_1 阶子块的个数, (JλI)k1(J - \lambda I)^{k_1} 的秩是 k1k_1 阶子块的个数的2倍再加上 k11k_1 - 1 阶子块的个数,如此等等。对于 i=0,1,2,,k11i = 0, 1, 2, \dots, k_1 - 1(JλI)k1i(J - \lambda I)^{k_1 - i} 的秩序列递归地确定出 JJ 中所有子块的阶数。

  1. 在一般的 Jordan 矩阵 (3.1.1.1) 中的所有 Jordan 子块的阶数可以通过分析某些幂的秩来确定。如果 λ1\lambda_{1} 是 Jordan 矩阵 JMnJ \in M_{n} 的特征值,那么,当我们取幂 (Jλ1I)(J - \lambda_{1}I)(Jλ1I)2(J - \lambda_{1}I)^{2} ,…时,只有相应子 λ1\lambda_{1} 的诸 Jordan 块会变成零,因为 Jλ1IJ - \lambda_{1}I 中的其他子块都有非零对角元。最终, (Jλ1I)k(J - \lambda_{1}I)^{k} 的秩将停止下降(不必考察任一 k>nk > n );使 (Jλ1I)k(J - \lambda_{1}I)^{k} 的秩达到它的极小值的那个最小的 kk 值是相应于 λ1\lambda_{1} 的最大子块的阶。这个极小值称为特征值 λ1\lambda_{1} 的指标。对 (Jλ1I)(J - \lambda_{1}I) 的各次幂组成的序列的秩进行分析,就足以确定相应于 λ1\lambda_{1} 的诸 Jordan 块的阶数与个数,然后再对 λ2,λ3\lambda_{2}, \lambda_{3} ,等等作同样的分析,因而确定出所有的子块。

131

虽然上述各个论断是针对Jordan矩阵 JJ 作出的,然而,每一个论断也同样适用于任一相似于 JJ 的矩阵。因此,如果已知矩阵 AMnA \in M_n ,那么 AA 的Jordan标准形(而不是把 AA 变换成Jordan标准形的相似)可以通过下述步骤来确定:

  1. 求出 AA 的所有不同的特征值,或许这可以通过求特征多项式的根来办到。

  2. AA 的每个不同的特征值 λi\lambda_{i} ,作 (AλiI)k(A - \lambda_{i}I)^{k} ,其中 k=1,2,,nk = 1, 2, \dots, n ,然后分析这些矩阵的秩所组成的序列便可看出 AA 的相应于特征值 λi\lambda_{i} 的所有 Jordan 块的阶数与个数。

这个算法对于手算一些简单形式的小矩阵可能是有用的,但是用计算机计算就可能失效,因为确定一个矩阵的秩的过程本来就是一个不稳定的过程。例如, Λε=[100ε]\Lambda_{\varepsilon} = \begin{bmatrix} 1 & 0 \\ 0 & \varepsilon \end{bmatrix} 就可以使我们看清这一点;对所有 ε0\varepsilon \neq 0AεA_{\varepsilon} 有秩2。但对 ε=0\varepsilon = 0 ,它有秩1。

例如,这个算法在下述情形是有用的,假定有这样一个问题,要确定 Jordan 块 Jk(0)MkJ_{k}(0) \in M_{k} 的平方的 Jordan 标准形:

A=[010100]2=[00101000].A = \left[ \begin{array}{c c c c} 0 & 1 & & 0 \\ & & \ddots & \\ & & \ddots & 1 \\ 0 & & & 0 \end{array} \right] ^ {2} = \left[ \begin{array}{c c c c c} 0 & 0 & 1 & \dots & 0 \\ & & & \ddots & \vdots \\ & & \ddots & & 1 \\ & & \ddots & \ddots & \\ & 0 & & \ddots & 0 \\ & & & & 0 \end{array} \right].

Λ\Lambda 的各特征值都是零,对 m[(k+1)/2]=(k+1)/2m - [(k + 1) / 2] = (k + 1) / 2 的整数部分, Am=0A^m = 0 ,而如果 p=1,2,,m1p = 1, 2, \dots, m - 1 ,则 Ap0A^p \neq 0 。对于 p1,2,,m2p - 1, 2, \dots, m - 2 ,每个幂 Ap+1A^{p+1} 的秩比它的前一个幂 Λp\Lambda^p 的秩少 2。如果 kk 是偶数,则 Λm1\Lambda^{m-1} 有秩 2,而如果 kk 为奇数,则 Am1A^{m-1} 有秩 1。因此, A=Jk2(0)A = J_k^2(0) 的 Jordan 标准形是

Jm(0)00Jm(0)],k=2m是 偶 数;\begin{array}{r} {\left. \begin{array}{l l} {- J _ {m} (0)} & {0} \\ {- 0} & {J _ {m} (0)} \end{array} \right], \quad k = 2 m \text {是 偶 数};} \end{array}

以及

[Jm(0)00Jm1(0)],k2m1是 奇 数.\left[ \begin{array}{l l} J _ {m} (0) & 0 \\ 0 & J _ {m - 1} (0) \end{array} \right], \quad k - 2 m - 1 \text {是 奇 数}.

这一结果对于确定一个已知矩阵是否有平方根是有用的。比如,它说明不存在使 A2[0100]A^2 - \left[ \begin{array}{ll}0 & 1\\ 0 & 0 \end{array} \right] 的矩阵 AM2A \in M_2

3.2.2 常系数线性微分方程组 在理论上具有重要意义的 Jordan 标准形的一个应用是分析常系数微分方程组的解。设 AMnA \in M_n 是已知的,且考虑一阶初值问题

λ(t)Ax(t).x(0)=x0是 已 知 的,(3.2.2.1)\begin{array}{l} \lambda^ {\prime} (t) \quad A x (t). \\ x (0) = x _ {0} \text {是 已 知 的}, \tag {3.2.2.1} \\ \end{array}

其中 x(t)=[x1(t),x2(t),,xn(t)]Tx(t) = [x_1(t), x_2(t), \dots, x_n(t)]^T ,且撇号 ()(^{\prime}) 表示对 tt 微分。如果 AA 不是对角矩阵,则这个方程组是耦合的;即 xi(t)x_i(t) 不仅与 xj(t)x_j(t) 有关,而且还与向量 x(t)x(t) 的其他分量有关。这种耦合情形使问题变得难以解决,但是,如果可以把 AA 变成对角(或近乎对角)的形式,那么耦合的影响可以减弱,或者甚至可以消除,因而可以很容易地求解。如果 A=SJS1A = SJS^{-1} ,其中 JJAA 的Jordan标准形,那么(3.2.2.1)变成

y(t)Jy(t),y(0)y0是 已 知 的,(3.2.2.2)\begin{array}{l} y ^ {\prime} (t) - J y (t), \\ y (0) - y _ {0} \text {是 已 知 的}, \tag {3.2.2.2} \\ \end{array}

其中 x(t)=Sy(t)x(t) = S y(t) ,且 y0=S1x0y_{0} = S^{-1} x_{0} 。如果可以解出问题(3.2.2.2),那么(3.2.2.1)的解的每个分量就是(3.2.2.2)的解的诸分量的一个线性组合,且这个线性组合是由 SS 给出的。

如果 Λ\pmb{\Lambda} 是可对角化的,那么 JJ 是对角矩阵,且(3.2.2.2)正好是一个形如 yk(t)=λkyk(t)y_{k}^{\prime}(t) = \lambda_{k}y_{k}(t) 的非耦合方程组,它有解 yk(t)yk(0)eλkty_{k}(t) - y_{k}(0)e^{\lambda_{k}t} 。如果特征值 λk\lambda_{k} 是实的,这是一个简单的幂,而如果 λk=ak+ibk\lambda_{k} = a_{k} + ib_{k} 是复数,它则是一个振动项 yk(t)=yk(0)eakt[cos(bkt)+isin(bkt)]y_{k}(t) = y_{k}(0)e^{a_{k}t}\left[\cos (b_{k}t) + i\sin (b_{k}t)\right]

如果 JJ 是非对角的,解就更复杂了。相应于 JJ 中的不同Jordan块的 y(t)y(t) 的诸分量是非耦合的,因此只要考虑 JJ 是单独的Jordan块

J=[λ1010λ]MmJ = \left[ \begin{array}{c c c c} \lambda & 1 & & 0 \\ & & \ddots & \\ & \ddots & & 1 \\ 0 & & & \lambda \end{array} \right] \in M _ {m}

的情形就可以了。方程组(3.2.2.2)是

y1(t)=λy1(t)+y2(t),ym1(t)=λym1(t)+ym(t),ym(t)=λym(t),\begin{array}{l} y _ {1} ^ {\prime} (t) = \lambda y _ {1} (t) + y _ {2} (t), \\ \vdots \qquad \vdots \qquad \vdots \\ y _ {m - 1} ^ {\prime} (t) = \lambda y _ {m - 1} (t) + y _ {m} (t), \\ y _ {m} ^ {\prime} (t) = \lambda y _ {m} (t), \\ \end{array}

[133

它可以由下往上直接求解。从最后一个方程开始,得到

ym(t)=ym(0)ea,y _ {m} (t) = y _ {m} (0) e ^ {a},

于是

ym1(t)=λym1(t)+ym(0)eλt.y _ {m - 1} ^ {\prime} (t) = \lambda y _ {m - 1} (t) + y _ {m} (0) e ^ {\lambda t}.

它有解

ym1(t)=ym(0)teλt+ym1(0)eλt,y _ {m - 1} (t) = y _ {m} (0) t e ^ {\lambda t} + y _ {m - 1} (0) e ^ {\lambda t},

于是可以把它用于下一个方程,它变成

ym2(t)=λym2(t)+ym(0)teit+ym1(0)eit.y _ {m - 2} ^ {\prime} (t) = \lambda y _ {m - 2} (t) + y _ {m} (0) t e ^ {i t} + y _ {m - 1} (0) e ^ {i t}.

它有解

ym2(t)=ym(0)t22eu+ym1(0)teu+ym2(0)eu,y _ {m - 2} (t) = y _ {m} (0) \frac {t ^ {2}}{2} e ^ {u} + y _ {m - 1} (0) t e ^ {u} + y _ {m - 2} (0) e ^ {u},

如此等等。显然解的每一个分量具有形式

yk(t)=eλtqk(t),y _ {k} (t) = e ^ {\lambda t} q _ {k} (t),

其中 qk(t)q_{k}(t) 是次数至多为 m1m - 1 的多项式.

从上述分析可以知道,对于任一给定的初始条件 x0x_0 ,问题(3.2.2.1)的解 x(t)x(t) 的各个分量有形式

xt(t)=p1(t)eλ1t+p2(t)eλ2t++pk(t)eλkt,x _ {t} (t) = p _ {1} (t) e ^ {\lambda_ {1} t} + p _ {2} (t) e ^ {\lambda_ {2} t} + \dots + p _ {k} (t) e ^ {\lambda_ {k} t},

其中 λ1,λ2,,λt\lambda_1, \lambda_2, \dots, \lambda_tAA 的不同的特征值,而每个 pi(t)p_i(t) 是多项式,它的次数严格小于相应于 λi\lambda_i 的最大 Jordan 块的阶数。实特征值产生纯指数项,而复特征值可以产生混合指数和振动项。

3.2.3 矩阵与它的转置的相似性 每个 Jordan 块置换相似于它的转置,因为从相似性可以看出

[01110][λ1010λ][01110]=[λ0101λ].\left[ \begin{array}{c c c c} 0 & & & 1 \\ & & \ddots & \\ & 1 & & \\ 1 & & & 0 \end{array} \right] \left[ \begin{array}{c c c c} \lambda & 1 & & 0 \\ & & \ddots & \\ & & \ddots & 1 \\ 0 & & & \lambda \end{array} \right] \left[ \begin{array}{c c c c} 0 & & & 1 \\ & & \ddots & \\ & 1 & & \\ 1 & & & 0 \end{array} \right] = \left[ \begin{array}{c c c c} \lambda & & & 0 \\ 1 & & \ddots & \\ & \ddots & & \\ 0 & & 1 & \lambda \end{array} \right].

因此,如果 AMnA \in M_{n} 是已知的,且 A=SJS1A = SJS^{1} 是它的Jordan标准形,那么, AA 相似于 JJ ,而 JJ 相似于 JTJ^{T}JTJ^{T} 又相似于 AT=(ST)1JT(ST)A^{T} = (S^{T})^{-1}J^{T}(S^{T}) 。结论是,每个复矩阵相似于它的转置。由此推出,一个复矩阵的行秩(线性无关行的最大个数)等于它的列秩(线性无关列的最大个数),这是因为秩是相似不变量。同时也推出, AAATA^{T} 有相同的特征值,不过,直接作出所有这些推论也是比较容易的。

对任意的域 F\mathbf{F} ,情形也是这样, Mn(F)M_{n}(\mathbf{F}) 中的每个矩阵都可经 Mn(F)M_{n}(\mathbf{F}) 中的一个矩阵相似于它的转置;且不必假定 F=C\mathbf{F} = \mathbf{C} 。实际上,相似矩阵可以取为对称矩阵。

[134]

3.2.4 交换矩阵和非减次矩阵 如果 p(t)p(t) 是多项式,而 AMnA \in M_n 是已知矩阵,显然 p(A)p(A)AA 可交换,这是一个有用的事实。可是反过来又如何呢?如果 A,BMnA, B \in M_n 是已知的,且 AABB 可交换, BB 必定是关于 AA 的多项式吗?当然不是这样,因为如果取 A=IA = I ,则 AA 与每个矩阵都可交换, p(I)=p(1)Ip(I) = p(1)I 不可能产生任何非纯量矩阵。问题是 AA 这一形式使它能与许多矩阵可交换,但又只允许它产生少数形如 p(A)p(A) 的矩阵。为了得到这方面的任何结果,必须寻求协调这两种能力之间的关系的某个折中办法。

3.2.4.1 定义 矩阵 AMnA \in M_{n} 叫做非减次矩阵,是指 AA 的每个特征值有几何重数 1.

因为 Jordan 矩阵的一个特征值的几何重数等于相应于该特征值的 Jordan 块的个数,所以,一个矩阵是非减次的,当且仅当相应于每个不同的特征值恰好有 - 个 Jordan 块。例如,如果 AMnA \in M_nnn 个不同的特征值,或者 AA 只有一个特征值,而该特征值有几何重数 1,那么它是非减次的。纯量矩阵不是非减次矩阵。

3.2.4.2 定理 设 AMnA \in M_{n} 是给定的非减次矩阵,矩阵 BMnB \in M_{n}AA 可交换,当且仅当存在一个次数至多为 n1n - 1 的多项式 p()p(\cdot) 使得 B=p(A)B = p(A) .

证明:如果 B=p(A)B = p(A) ,那么 BB 一定与 BB 可交换.关于逆命题,设 A=SJS1A = SJS^{-1}AA 的Jordan标准形.如果 BA=ABBA = AB ,那么 BSJS1=SJS1BBSJS^{-1} = SJS^{-1}B(S1BS)J=J(S1BS)(S^{-1}BS)J = J(S^{-1}BS) ,如果能够证明 S1BS=p(J)S^{-1}BS = p(J) ,那么 BSp(J)S1=p(SJS1)=p(A)B - Sp(J)S^{-1} = p(SJS^{-1}) = p(A) ,因此,只要假定 AA 本身就是一个Jordan矩阵就可以了.因为 AA 是非减次的,所以可以假定

A=[Jn1(λ1)0Jn2(λ2)0Jnk(λk)],A = \left[ \begin{array}{c c c c} J _ {n _ {1}} (\lambda_ {1}) & & & 0 \\ & J _ {n _ {2}} (\lambda_ {2}) & & \\ & & \ddots & \\ 0 & & & J _ {n _ {k}} (\lambda_ {k}) \end{array} \right],

其中 λ1,λ2,,λk\lambda_1, \lambda_2, \dots, \lambda_kAAkk 个不同的特征值,如果把 BB 写成与 AA 的上述形式相同的分块形式 B=(Bn)B = (B_n) ,那么 ABBAAB - BA 的相应非对角子块具有形式

Jn(λi)BijBijJnj(λj)=0,ij.J _ {n} \left(\lambda_ {i}\right) B _ {i j} - B _ {i j} J _ {n _ {j}} \left(\lambda_ {j}\right) = 0, \quad i \neq j.

因为特征值 λ1\lambda_{1}λj\lambda_{j} 是不同的,由此可以得出[见(2.4)节习题9], Bij=0B_{ij} = 0 是这些方程的唯一解,因此矩阵 BB 必须是分块对角矩阵

B=[B10B20Bk],B = \left[ \begin{array}{c c c c} B _ {1} & & & 0 \\ & B _ {2} & & \\ & & \ddots & \\ 0 & & & B _ {k} \end{array} \right],

其中每个 BiMniB_{i} \in M_{n_{i}} 。交换性假设是说,对于所有 i=1,2,,ki = 1, 2, \dots, k ,有 BiJni(λi)=Jni(λi)BiB_{i} J_{n_{i}}(\lambda_{i}) = J_{n_{i}}(\lambda_{i}) B_{i} 。如果记 Jni(λi)=λiI+NiJ_{n_{i}}(\lambda_{i}) = \lambda_{i} I + N_{i} ,其中

Ni=[010100]Mn,N _ {i} = \left[ \begin{array}{c c c} 0 & 1 & 0 \\ & \ddots & \\ & \ddots & 1 \\ 0 & & 0 \end{array} \right] \in M _ {n},

那么这些恒等式变成 BiNi=NiBiB_{i}N_{i} = N_{i}B_{i}i=1i = 1 ,2,…,k.直接计算可知,由于 NiN_{i} 的特殊形式每个 BiB_{i} 必须是具有Toeplitz形式(0.9.7)的上三角矩阵,即

Bt=[b1(t)b2(t)bn1(t)b1(t)0b2(t)b1(t)],(3.2.4.3)B _ {t} = \left[ \begin{array}{c c c c} b _ {1} ^ {(t)} & b _ {2} ^ {(t)} & \dots & b _ {n _ {1}} ^ {(t)} \\ & b _ {1} ^ {(t)} & \ddots & \vdots \\ 0 & & \ddots & b _ {2} ^ {(t)} \\ & & & b _ {1} ^ {(t)} \end{array} \right], \tag {3.2.4.3}

136

其中沿着各对角线的诸元素取常值。如果可以构造出一些次数至多是 n1n - 1 的多项式 pi(t)p_{i}(t) ,且具有性质:对所有 iji \neq jpi(Jκi(λj))=0p_{i}(J_{\kappa_{i}}(\lambda_{j})) = 0pi(Jκi(λi))=Bip_{i}(J_{\kappa_{i}}(\lambda_{i})) = B_{i} ,那么

p(t)=p1(t)++pk(t)p (t) = p _ {1} (t) + \dots + p _ {k} (t)

将满足定理的论断. 定义

qt(t)=i=1j=1k(tλj)ni,qt(t)的 次 数=nni,q _ {t} (t) = \prod_ {\substack {i = 1 \\ j = 1}} ^ {k} (t - \lambda_ {j}) ^ {n _ {i}}, \quad q _ {t} (t) \text {的 次 数} = n - n _ {i},

由此可见,对所有 iji \neq jqi(Jnj(λj))=0q_{i}(J_{n_{j}}(\lambda_{j})) = 0 ,这是因为 (Jnj(λj)λjI)nj=0(J_{n_{j}}(\lambda_{j}) - \lambda_{j}I)^{n_{j}} = 0 。显然 qi(Jni(λi))q_{i}(J_{n_{i}}(\lambda_{i})) 不一定等于 BiB_{i} ,可是它是非奇异矩阵(因为诸 λi\lambda_{i} 是各不相同的),且像关于 Jni(λi)J_{n_{i}}(\lambda_{i}) 的任一个多项式那样,具有形式(3.2.4.3).

因为形如(3.2.4.3)的任一非奇异矩阵的逆具有同一形式,又因为这种形式的任意两个矩阵的乘积还是这同一形式,所以矩阵

[qi(Jn(λi))]1Bi\left[ q _ {i} \left(J _ {n} (\lambda_ {i})\right) \right] ^ {- 1} B _ {i}

是具有Toeplitz形式(3.2.4.3)的上三角矩阵,每个这样的矩阵可以写成关于 Jni(λi)J_{n_i}(\lambda_i) 的一个多项式.例如

Btb1(1)(Jnt(λt)λtI)c+b2(1)(Jnt(λt)λtI)1+bnt(1)(Jnt(λt)λtI)nt.B _ {t} - b _ {1} ^ {(1)} \left(J _ {n _ {t}} (\lambda_ {t}) - \lambda_ {t} I\right) ^ {c} + b _ {2} ^ {(1)} \left(J _ {n _ {t}} (\lambda_ {t}) \quad \lambda_ {t} I\right) ^ {1} + \dots - b _ {n _ {t}} ^ {(1)} \left(J _ {n _ {t}} (\lambda_ {t}) \quad \lambda_ {t} I\right) ^ {n _ {t}}.

因此,存在一个次数至多是 ni1n_{i} - 1 的多项式 ri(t)r_i(t) ,使得

[qi(Jn1(λi))]1Bi=ri(Jni(λi)).\left[ q _ {i} \left(J _ {n _ {1}} (\lambda_ {i})\right) \right] ^ {1} B _ {i} = r _ {i} \left(J _ {n _ {i}} (\lambda_ {i})\right).

如果现在令 pi(t)=qi(t)ri(t)p_{i}(t) = q_{i}(t)r_{i}(t)pi(t)p_{i}(t) 的次数至多是 n1n - 1

pr(Jni(λj))qr(Jnj(λj))rr(Jnj(λj))=0rr(Jnj(λj))=0,p _ {r} \left(J _ {n _ {i}} \left(\lambda_ {j}\right)\right) - q _ {r} \left(J _ {n _ {j}} \left(\lambda_ {j}\right)\right) r _ {r} \left(J _ {n _ {j}} \left(\lambda_ {j}\right)\right) = 0 r _ {r} \left(J _ {n _ {j}} \left(\lambda_ {j}\right)\right) = 0,

如果 iji \neq j ,则 pi(Jni(λi))qi(Jni(λi))ri(Jni(λi))=Bip_{i}(J_{n_{i}}(\lambda_{i})) - q_{i}(J_{n_{i}}(\lambda_{i}))r_{i}(J_{n_{i}}(\lambda_{i})) = B_{i}

定理的逆命题也同样成立,因而导出非减次矩阵的一个特征:矩阵 AMnA \in M_{n} 是非减次的,当且仅当每个与 AA 可交换的矩阵是 AA 的一个多项式.

3.2.5 收敛矩阵 矩阵 AMnA \in M_{n} 如果有具性质:当 mm \to \infty 时, AmA^{m} 的所有元素都趋于零,就称 AA 为收敛矩阵。这样的矩阵在数值线性代数的算法分析中起着重要的作用。如果 AA 是对角矩阵,那么显然有: AA 收敛,当且仅当 AA 的所有特征值的模都严格小于 1,并且同样的结果可以推广到可对角化矩阵。

[137]

由于涉及极限运算,如何运用扰动理论把这一结果推广到未必是对角化矩阵的一般情形是不明显的。但是,可以利用Jordan标准形,如果 A=SJS1A = SJS^{-1}AA 的Jordan标准形,那么 Am=SJmS1A^{m} = SJ^{m}S^{-1} ,因而,当 mm \to \inftyAm0A^{m} \to 0 当且仅当 mm \to \inftyJm0J^{m} \to 0 。因为 JJ 是诸Jordan块的直和,

所以只要考虑Jordan块

Jk(λ)=[λ1010λ]=[λ00λ]+[010100]=λINkMk,其 中Nk=Jk(0)\begin{array}{l} J _ {k} (\lambda) = \left[ \begin{array}{c c c} \lambda & 1 & 0 \\ & \ddots & \\ & \ddots & 1 \\ 0 & & \lambda \end{array} \right] = \left[ \begin{array}{c c c} \lambda & & 0 \\ & \ddots & \\ 0 & & \lambda \end{array} \right] + \left[ \begin{array}{c c c} 0 & 1 & 0 \\ & & \\ & \ddots & \\ & \ddots & 1 \\ 0 & & 0 \end{array} \right] \\ = \lambda I \cdot N _ {k} \in M _ {k}, \quad \text {其 中} N _ {k} = J _ {k} (0) \\ \end{array}

的各次幂的性质就可以了,因为,对所有 mkm \geq kNkm=0\mathbf{N}_k^m = 0 ,有

[Jk(λ)]m=(λI+Nk)m=i=1m(mi)λiNkmi=i=mm(mi)λiNkmi.\left[ J _ {k} (\lambda) \right] ^ {m} = \left(\lambda I + N _ {k}\right) ^ {m} = \sum_ {i = 1} ^ {m} \binom {m} {i} \lambda^ {i} N _ {k} ^ {m - i} = \sum_ {i = m} ^ {m} \binom {m} {i} \lambda^ {i} N _ {k} ^ {m - i}.

对所有的 mkm \geqslant k 成立,因为各对角元都为 λm\lambda^m ,如果 Jm0J^{m} \to 0 ,则必定有 λm0\lambda^m \to 0 ,这说明 λ<1|\lambda| < 1 。反过来,如果 λ<1|\lambda| < 1 ,我们想证明,对每个 j=0,1,2,,k1j = 0, 1, 2, \dots, k-1

m:×;,(mmj)λmj0.\text {当} m \rightarrow : \times ; \text {时}, \quad \binom {m} {m - j} \lambda^ {m - j} \rightarrow 0.

然而,

(mmj)λm1=m(m1)(m2)(mj+1)λmj!λjmλmj!λj,\left| \binom {m} {m - j} \lambda^ {m - 1} \right| = \left| \frac {m (m - 1) (m - 2) \cdots (m - j + 1) \lambda^ {m}}{j ! \lambda^ {j}} \right| \leqslant \left| \frac {m ^ {\prime} \lambda^ {m}}{j ! \lambda^ {j}} \right|,

因此,只要证明当 mm \to \infty 时有 mλm0m^{\prime} \mid \lambda^{m} \to 0 就够了。看出这个关系的一个简单办法是取对数,并且注意到,因为 logλ<0\log |\lambda| < 0 ,又根据L'Hopital法则,当 xx \to \infty(logx)/x0(\log x) / x \to 0 ,所以当 mm \to \infty 时,

jlogmmlogλ.j \log m \mid m \log | \lambda | \rightarrow - \infty .

上述证法主要是利用A的Jordan标准形来证明,当 m+m \rightarrow +\infty 时有 Am0A^{m} \rightarrow 0 。当且仅当A的所有特征值有小于1的模。另一个证明将在(5.6.12)中给出,它根本不依赖Jordan标准形。

3.2.6 几何重数-代数重数不等式 矩阵 AMnA \in M_{n} 的一个特征值的几何重数是 A 的相应于该特征值的 Jordan 块的个数。这个数小于或等于相应于该特征值的所有 Jordan 块的阶数的和。这个和是代数重数。因此一个特征值的几何重数不大于它的代数重数;试与(1.4.9)比较。

3.2.7 可对角化矩阵和幂零矩阵 一个矩阵 AMnA \in M_{n} 是幂零的,是指对某个正整数 kk ,有 Ak=0A^{k} = 0 。任一Jordan块 Jk(λ)J_{k}(\lambda) 可以写成 Jk(λ)=λI+NkJ_{k}(\lambda) = \lambda I + N_{k} ,其中, (Nk)k=0(N_{k})^{k} = 0 。因此,任一Jordan块是一个对角矩阵矩与一个幂零矩阵的和。

更一般地,一个Jordan矩阵(3.2.1.1)可以写成 J=D+NJ = D + N ,其中 DD 是对角矩阵,它的主对角线与 JJ 的主对角线相同,且 N=JDN = J - D 。矩阵 NN 是幂零的,且只要 kkJJ 的最大Jordan块的阶,就有 Nk=0N^{k} = 0

最后,如果 AMnA \in M_{n} 是任一给定的矩阵,且 A=SJS1A = SJS^{-1} 是Jordan标准形,那么 A=SDS1+SNS1=AD+ANA = SDS^{-1} + SNS^{-1} = A_{D} + A_{N} ,其中 ADA_{D} 可对角化, ANA_{N} 是幂零的,且 ADAN=ANADA_{D}A_{N} = A_{N}A_{D} ,这是因为 DDNN 都是具有相应的同阶子块的分块对角矩阵,且 DD 的诸子块是纯量矩阵。

我们得出,任一 AMnA \in M_{n} 可以表示成一个可对角化矩阵与一个幂零矩阵的和,且这两个被

加矩阵可交换.

习题

  1. rˉ={Aα:αf}Mn\bar{r} = \{A_{\alpha} : \alpha \in f\} \subset M_{n} 是给定的矩阵族,且附以指标集 yˉ\bar{y} ,又假定存在非减次矩阵 A1yˉA_{1} \in \bar{y} ,使得对所有 αf\alpha \in f ,有 AαAα=AαAαA_{\alpha} A_{\alpha} = A_{\alpha} A_{\alpha} ,证明,对每个 αyˉ\alpha \in \bar{y} ,都存在一个次数至多为 n1n - 1 的多项式 pα(t)p_{\alpha}(t) 使得 Aα=pα(Aα)A_{\alpha} = p_{\alpha}(A_{\alpha}) ,因而 yˉ\bar{y} 是一个交换族。

  2. AMnA \in M_{n} 是已知的, λi\lambda_{i}AA 的一个特征值。证明 AA 的相应于特征值 λi\lambda_{i} 的最大 Jordan 块的阶数 ( λi\lambda_{i} 的指标) 是使 rank(AλiI)k=rank(AλiI)k+1\operatorname{rank}(A - \lambda_{i}I)^{k} = \operatorname{rank}(A - \lambda_{i}I)^{k + 1}k=1,2,,n1k = 1, 2, \dots, n - 1 的最小值。

  3. 如果对某个 knk \geq n , AMnA \in M_{n}Ak=0A^{k} = 0 , 证明, 对某个 rnr \leq n , 有 Ar=0A^{r} = 0 , 因此每个幂零矩阵有一个次数不大于其阶数的等于零的幂. 提示: 证明 0 是 AA 仅有的特征值. 想一想 AA 的 Jordan 标准形是什么样子? 再取幂.

  4. Jk(0)J_{k}(0) 是一个已知的 Jordan 块。试利用(3.2.1)末尾的证法确定矩阵 Jk3(0)J_{k}^{3}(0) 的三种可能的 Jordan 标准形。

  5. AMnA \in M_{n} 是幂零矩阵,因而对某个 kkAk=0A^{k} = 0 。证明, AA 的特征多项式是 pA(t)=tnp_{A}(t) = t^{n}

  6. 作用在次数至多为3的全体多项式所组成的向量空间上的线性变换 d/dt\mathrm{d} / \mathrm{d}tp(t)p(t)p(t) \rightarrow p'(t) ,在基 β={1,t,t2,t3}\beta = \{1, t, t^2, t^3\} 下有基表示

[0100002000030000].\left[ \begin{array}{c c c c} 0 & 1 & 0 & 0 \\ 0 & 0 & 2 & 0 \\ 0 & 0 & 0 & 3 \\ 0 & 0 & 0 & 0 \end{array} \right].

这个矩阵的Jordan标准形是什么?

  1. 使 Aδ=IA^{\delta} = I 的矩阵 AMnA \in M_{n} 的各种可能的 Jordan 形是什么?

  2. 只有特征多项式 p1(t)=(t+3)t(t4)2p_1(t) = (t + 3)^t (t - 4)^2 的矩阵 AM0A \in M_0 的各种可能的 Jordan 标准形是什么?

  3. 试用(3.2.1)中所叙述的方法确定

A=[i11i]A = \left[ \begin{array}{c c} i & 1 \\ 1 & - i \end{array} \right]

的Jordan标准形,

  1. 试证定理(3.2.4.2)的证明中的论断:形如(3.2.4.3)的两个矩阵的乘积有同一形式,再证明一个形如(3.2.4.3)的非奇异矩阵的逆有同一形式。提示: AA 的逆是一个关于 AA 的多项式。

  2. A,BMnA, B \in M_{n} , 试用定理(1.3.20)的证明中的恒等式证明, ABABBABA 的 Jordan 形中的各非奇异子块是恒等的. 这能说明 ABABBABA 相似吗? 如果 ABABBABA 不相似, 讨论这与它们能够相似有多大差别.

  3. 假定 A1,,AkA_{1}, \cdots, A_{k} 是已知矩阵,其中, AiMnA_{i} \in M_{n}i=1,2,,ki = 1, 2, \cdots, k ,又假定 J1,,JkJ_{1}, \cdots, J_{k} 是它们的相应 Jordan 标准形。证明直和

[A10A20Ak]Mn,n1+n2++nk=n\left[ \begin{array}{c c c c} A _ {1} & & & 0 \\ & A _ {2} & & \\ & & \ddots & \\ 0 & & & A _ {k} \end{array} \right] \in M _ {n}, \quad n _ {1} + n _ {2} + \dots + n _ {k} = n

(除了差诸对角子块的一个排列以外)有Jordan标准形

[J10J20Jk].\left[ \begin{array}{c c c c} J _ {1} & & & 0 \\ & J _ {2} & & \\ 0 & & \ddots & \\ & & & J _ {k} \end{array} \right].
  1. AMnA \in M_{n}B,CMmB, C \in M_{m} 是已知的,证明,直和 [A00B]Mn+m\left[ \begin{array}{cc} A & 0 \\ 0 & B \end{array} \right] \in M_{n+m} 相似于直和 [A00C]\left[ \begin{array}{cc} A & 0 \\ 0 & C \end{array} \right] ,当且仅当 BB 相似于 CC

  2. B,CMmB, C \in M_{m} 是已知的,证明, kk 次直和

[B0B0B]Mkn[C0C0C]Mkn,k1\left[ \begin{array}{c c c c} {B} & & & {0} \\ & {B} & & \\ {0} & & {\ddots} & \\ & & & {B} \end{array} \right] \in M _ {k n} \text {与} \left[ \begin{array}{c c c c} {C} & & & {0} \\ & {C} & & \\ {0} & & {\ddots} & \\ & & & {C} \end{array} \right] \in M _ {k n}, k \geqslant 1

相似,当且仅当 BBCC 相似

  1. AMnA \in M_{n}B,CMmB, C \in M_{m} 是已知的,证明,直和

[A0B0B]Mn+kn[A0C0C]Mnkn,k1\left[ \begin{array}{c c c c} {A} & & & {0} \\ & {B} & & \\ {0} & & {\ddots} & \\ & & & {B} \end{array} \right] \in M _ {n + k n} \text {与} \left[ \begin{array}{c c c c} {A} & & & {0} \\ & {C} & & \\ {0} & & {\ddots} & \\ & & & {C} \end{array} \right] \in M _ {n \cdot k n}, \quad k \geqslant 1

相似,当且仅当 BBCC 相似.

  1. AMnA \in M_{n} 有 Jordan 标准形 Jn1(λ1)Jnk(λk)J_{n_{1}}(\lambda_{1}) \oplus \cdots \oplus J_{n_{k}}(\lambda_{k}) 。如果 AA 是非奇异形矩阵,证明 A2A^{2} 的 Jordan 标准形是 Jn1(λ12)Jnk(λk2)J_{n_{1}}(\lambda_{1}^{2}) \oplus \cdots \oplus J_{n_{k}}(\lambda_{k}^{2}) ;即 A2A^{2} 的 Jordan 标准形恰好由与 AA 相同的一组 Jordan 块组成,只是相应的特征值要取平方。与此类似的结论对所有幂 Ak(k2)A^{k} (k \geqslant 2) 成立吗?举一个 2×22 \times 2 的例子说明,如果 AA 是奇异矩阵,上述结论不成立。提示:如果 A0A \neq 0 ,试证 Jk2(λ)J_{k}^{2}(\lambda) 的 Jordan 标准形(一个 Jordan 块)是 Jk(λ2)J_{k}(\lambda^{2}) 。再证明,如果 λ0\lambda \neq 0

rank[Jk(λ)λI]m=rank[Jk2(λ)λ2I]m,m=1,2,,k\operatorname {r a n k} \left[ J _ {k} (\lambda) - \lambda I \right] ^ {m} = \operatorname {r a n k} \left[ J _ {k} ^ {2} (\lambda) - \lambda^ {2} I \right] ^ {m}, \quad m = 1, 2, \dots , k
  1. 如果 AMnA \in M_{n} , 证明, rankA=rankA\operatorname{rank} A = \operatorname{rank} A^{*} , 当且仅当特征值 λ=0\lambda = 0 的几何重数与代数重数相等, 即 AA 的 Jordan 标准形中相应于 λ=0\lambda = 0 的所有 Jordan 块 (如果有的话) 是 1×11 \times 1 的.

进一步阅读一个给定的矩阵与其转置间的相似变换矩阵总可以取为对称矩阵。该事实的证明可参看0. Taussky and H. Zassenhaus, “On the Similarity Transformation between a Matrix and Its Transpose,” Pacific J. Math. 9(1959), 893-896. 在(3.2.4)节末提到了定理(3.2.4.2)的逆命题也成立,其证明可见[HJ].

3.2_Jordan标准形:若干论断和应用 - 矩阵分析 | OpenTech