• Jordan 标准型定理


    将学习到什么

    就算两个矩阵有相同的特征多项式,它们也有可能不相似,那么如何判断两个矩阵是相似的?答案是它们有一样的 Jordan 标准型.


    Jordan 标准型定理

    这节目的:证明**每个复矩阵都与一个本质上唯一的 Jordan 矩阵相似**. 分三步证明这个结论。其中前两步已经在其它章节中给出, - 第一步 [Schur 定理](http://www.cnblogs.com/zhoukui/p/7782688.html) 确保每个复矩阵都相似于一个上三角矩阵,这个上三角矩阵的特征值出现在其对角线上,且相等的特征值放在一起. - 第二步 [Schur 三角化定理推论](http://www.cnblogs.com/zhoukui/p/7786648.html) 中定理 1.3 确保第一步中所描述的那种形式的矩阵相似于一个分块对角的上三角矩阵, 其中每个对角分块都有相等的对角元素. - 第三步 在本节,要证明:**有相等对角元素的上三角矩阵相似于一个 Jordan 矩阵**.

    同时我们还对如下结论感兴趣:如果一个矩阵是实的,且只有实的特征值,那么它可以通过实相似化简为 Jordan 矩阵. 如果实矩阵 (A) 只有实的特征值,那么第一步与第二步确保存在一个相似矩阵 (S), 使得 (S^{-1}AS) 是一个实的分块对角的上三角矩阵. 于是,只要证明有相等主对角元素的实的上三角矩阵可以通过实相似化简成 Jordan 块的直和就足够了.

    下面给出一个有用的引理,它的证明完全是直接的计算. 特征值为零的 (k imes k) Jordan 块称为 幂零 Jordan 块.


      证明:如果 (n=1), 就有 (A=[0]), 故而结论是平凡的. 我们对 (n) 用归纳法. 假设 (n>1) 且结论对所有阶小于 (n) 的所有严格上三角矩阵都成立. 分划 (A=egin{bmatrix} 0 & a^T \ 0 & A_1 end{bmatrix}), 其中 (a in mathbb{C}^{n-1}), 且 (A_1 in M_{n-1}) 是严格上三角矩阵. 根据归纳假设,存在一个非奇异的 (S_1 in M_{n-1}), 使得 (S_1^{-1}A_1S_1) 有形式
    egin{align}
    S_1^{-1}A_1S_1= egin{bmatrix} J_{k_1} && \ & ddots & \ && J_{k_s} end{bmatrix} = egin{bmatrix} J_{k_1} & 0 \ 0 & J end{bmatrix}
    end{align}
    其中 (k_1 geqslant k_2 geqslant cdots geqslant k_s geqslant 1) , (k_1+k_2+cdots+k_s=n-1), (J_{k_i}=J_{k_i}(0)), 且 (J=J_{k_2}oplus cdots oplus J_{k_s} in M_{n-k_1-1}). (J) 中没有大于 (k_1) 的对角块,所以 (J^{k_1}=0). 计算表明
    egin{align} label{e1}
    egin{bmatrix} 1 & 0 \ 0 & S_1^{-1} end{bmatrix} A egin{bmatrix} 1 & 0 \ 0 & S_1 end{bmatrix} =egin{bmatrix} 0 & a^TS_1 \ 0 & S_1^{-1}A_1 S_1 end{bmatrix}
    end{align}
    分划 (a^TS_1=[a_1^T quad a_2^T]), 其中 (a_1 in mathbb{C}^{k_1}). 而 (a_2 in mathbb{C}^{n-k_1-1}), 并将 ef{e1} 写成
    egin{align}
    egin{bmatrix} 1 & 0 \ 0 & S_1^{-1} end{bmatrix} A egin{bmatrix} 1 & 0 \ 0 & S_1 end{bmatrix} =egin{bmatrix} 0 & a_1^T & a_2^T \ 0 & J_{k_1} & 0 \ 0 & 0 & J end{bmatrix}
    end{align}
    现在考虑相似性
    egin{align}
    & egin{bmatrix} 1 & -a_1^T J_{k_1}^T & 0 \ 0 & I & 0 \ 0 & 0 & I end{bmatrix} egin{bmatrix} 0 & a_1^T & a_2^T \ 0 & J_{k_1} & 0 \ 0 & 0 & J end{bmatrix} egin{bmatrix} 1 & a_1^TJ_{k_1}^T & 0 \ 0 & I & 0 \ 0 & 0 & I end{bmatrix} otag \ label{e22}
    &= egin{bmatrix} 0 & a_1^T(I-J_{k_1}^TJ_{k_1}) & a_2^T \ 0 & J_{k_1} & 0 \ 0 & 0 & J end{bmatrix}= egin{bmatrix} 0 & (a_1^Te_1)e_1^T & a_2^T \ 0 & J_{k_1} & 0 \ 0 & 0 & J end{bmatrix}
    end{align}
    其中用到恒等式 ([I-J_k^T J_k]x=(x^Te_1)e_1). 现在有两种可能性,它们要根据 (a_1^Te_1 eq 0) 还是 (a_1^Te_1 = 0) 而定.
    如果 (a_1^Te_1 eq 0), 那么
    egin{align}
    & egin{bmatrix} 1/a_1^Te_1 & 0 & 0 \ 0 & I & 0 \ 0 & 0 & (1/a_1^Te_1) I end{bmatrix} egin{bmatrix} 0 & (a_1^Te_1)e_1^T & a_2^T \ 0 & J_{k_1} & 0 \ 0 & 0 & J end{bmatrix}
    egin{bmatrix} a_1^Te_1 & 0 & 0 \ 0 & I & 0 \ 0 & 0 & a_1^Te_1 I end{bmatrix} otag \
    &= egin{bmatrix} 0 & e_1^T & a_2^T \ 0 & J_{k_1} & 0 \ 0 & 0 & J end{bmatrix}= egin{bmatrix} ilde{J} & e_1a_2^T \ 0 & J end{bmatrix} otag
    end{align}
    注意到 ( ilde{J}=egin{bmatrix} 0 & e_1^T \ 0 & J_{k_1} end{bmatrix}=J_{k_1+1}(0)). 由于 ( ilde{J}e_{i+1}=e_i)(对 (i=1,2,cdots,k_1)), 计算表明
    egin{align}
    egin{bmatrix} I & e_2a_2^T \ 0 & I end{bmatrix} egin{bmatrix} ilde{J} & e_1a_2^T \ 0 & J end{bmatrix} egin{bmatrix} I & -e_2a_2^T \ 0 & I end{bmatrix} =egin{bmatrix} ilde{J} & - ilde{J} e_2a_2^T+e_1a_2^T+e_2a_2^TJ \ 0 & J end{bmatrix} = egin{bmatrix} ilde{J} & e_2a_2^T J \ 0 & J end{bmatrix} otag
    end{align}
    我们可以用递推方式对 (i=2,3,cdots) 计算相似矩阵序列
    egin{align}
    egin{bmatrix} I & e_{i+1}a_2^T J^{i-1} \ 0 & I end{bmatrix} egin{bmatrix} ilde{J} & e_i a_2^T J^{i-1} \ 0 & J end{bmatrix} egin{bmatrix} I & -e_{i+1}a_2^T J^{i-1}\ 0 & I end{bmatrix} = egin{bmatrix} ilde{J} & e_{i+1}a_2^T J^i \ 0 & J end{bmatrix} otag
    end{align}
    由于 (J^{k_1} = 0),经过至多 (k_1) 步计算之后,这个相似矩阵序列中不在对角线上的元素最终都变为零. 我们得出结论: (A) 相似于 $ egin{bmatrix} ilde{J} & 0 \ 0 & J end{bmatrix} $, 它就是一个有所要求形状的严格上三角 Jordan 矩阵.
    如果 (a_1^Te_1 = 0), 则 ef{e22} 表明 (A) 相似于
    egin{align}
    egin{bmatrix} 0 & 0 & a_2^T \ 0 & J_{k_1} & 0 \ 0 & 0 & J end{bmatrix} otag
    end{align}
    它与
    egin{align} label{e3}
    egin{bmatrix} J_{k_1} & 0 & 0 \ 0 & 0 & a_2^T \ 0 & 0 & J end{bmatrix}
    end{align}
    置换相似. 根据归纳假设,存在一个非奇异的 (S_2 in M_{n-k_1}), 使得 (S_2^{-1}egin{bmatrix} 0 & a_2^T \ 0 & J end{bmatrix} S_2=hat{J} in M_{n-k_1}) 是主对角线为零的 Jordan 矩阵. 这样一来,矩阵 ef{e3} , (A) 本身都与 (egin{bmatrix} J_{k_1} & 0 \ 0 & hat{J} end{bmatrix}) 相似,这就是所要求形式的 Jordan 矩阵,除了对角 Jordan 块有可能不是按照其大小非增的次序排列. 如果需要的话,用分块置换相似就产生出所需要的形式.
    最后注意到,如果 (A) 是实的,那么这一证明中所有相似矩阵都是实的,所以 (A) 通过实相似与所要求形式的 Jordan 矩阵相似.

    由于一般情形是幂零情形的简单推论,所以定理 (1.1) 基本完成了第 3 步. 如果 (A in M_n) 是所有对角元素都为 (lambda) 的上三角矩阵,那么 (A_0=A-lambda I) 就是一个严格上三角矩阵. 如果 (S in M_n) 是非奇异的,且 (S^{-1}A_0S) 是幂零 Jordan 块 (J_{n_i}(0)) 的直和,那么 (S^{-1}AS=S^{-1}A_0S+lambda I) 就是特征值为 (lambda) 的 Jordan 块 (J_{n_i}(lambda)) 的直和. 下面给出 Jordan 标准型定理 中的存在性结论.

    定理 (1.2) 中的 Jordan 矩阵就是 (A)Jordan 标准型(相差不过直和项的排列). 矩阵 (J_k(lambda),lambda in mathbb{C},k=1,2,cdots)相似的标准分块. 两个事实为理解 Jordan 标准型定理中的唯一性提供了关键思路:(1) 如果两个矩阵都用同样的纯量矩阵来做平移,其相似性保持不变;(2) 秩是一个相似不变量.
     


    一些重要结论

    如果 (A,B,Sin M_n), (S) 是非奇异的,且 (A=SBS^{-1}), 那么对任何 (lambdain mathbb{C}), 有 (A-lambda I=SBS^{-1}-lambda SS^{-1}=S(B-lambda I)S^{-1}). 对每个 (k=1,2,cdots), 矩阵 ((A-lambda I)^k)((B-lambda I)^k) 相似,当然它们的秩相等. 我们着手于与 (A) 相似的 Jordan 矩阵 (B=J=J_{n_1}(lambda_1) oplus cdots oplus J_{n_q}(lambda_q)) 以及 (lambda)(A) 的一个特征值时的结论. 在对 (J) 的对角块用置换相似作排列之后。可以假设 (J=J_{m_1}(lambda) oplus cdots oplus J_{m_p}(lambda)oplus hat{J}), 其中的 Jordan 矩阵 (hat{J}) 是异于 (lambda) 的特征值对应的 Jordan 块的直和. 这样 (A-lambda I) 就相似于
    egin{align}
    J-lambda I&=(J_{m_1}(lambda)-lambda I) oplus cdots oplus (J_{m_p}(lambda)-lambda I)oplus( hat{J}-lambda I) otag \
    &=J_{m_1}(0) oplus cdots oplus J_{m_p}(0)oplus ( hat{J}-lambda I) otag
    end{align}
    它是 (p) 个不同阶的幂零 Jordan 块和一个非奇异 Jordan 矩阵 $ hat{J}-lambda I in M_m$ 的直和, 其中 (m=n-(m_1+cdots+m_p)). 此外,对每个 (k=1,2,cdots), ((A-lambda I)^k) 与 $(J-lambda I)^k=J_{m_1}(0)^k oplus cdots oplus J_{m_p}(0)^koplus ( hat{J}-lambda I)^k $ 相似. 由于直和的秩等于各个直和项的秩之和,故而对每个 (k=1,2,cdots), 我们有
    egin{align}
    mathrm{rank}\,(A-lambda I)^k &=mathrm{rank}\,(J-lambda I)^k= mathrm{rank}\, J_{m_1}(0)^k oplus cdots oplus mathrm{rank}\,J_{m_p}(0)^k oplus mathrm{rank}\, (hat{J}-lambda I)^k otag \ label{equ444}
    &=mathrm{rank}\, J_{m_1}(0)^k oplus cdots oplus mathrm{rank} \, J_{m_p}(0)^k +m
    end{align}
    幂零 Jordan 块的幂的秩等于什么?(J_{ell}(0)) 的第一列是零,而它的后 (ell-1) 列是线性无关的(第一条超对角元素中仅有的非零元素是 1),所以 (mathrm{rank}\, J_{ell}(0)=ell-1). (J_{ell}(0)^2) 中仅有的非零元素是在第二条超对角线中的 1,所以它的前两列都是零,它的后 (ell-2) 列是线性无关的,且 (mathrm{rank}\, J_{ell}(0)^2=ell-2). 随着幂每提高一次,这些 1 就向上移动一条超对角线(所以为零的列的个数就增加 1,而秩则减少 1),直到 (J_{ell}(0)^{ell-1}) 时恰好只有一个非零的元素(在位置 ((1,ell)) 处),且 (mathrm{rank}\,J_{ell}(0)^{ell-1}=1=ell-(ell-1)). 当然,对所有 (k=ell,ell+1,cdots), 都有 (J_{ell}(0)^k=0). 一般来说,对每个 (k=1,2,cdots), 我们有 (mathrm{rank}\,J_{ell}(0)^k=max {ell-k,0}), 所以
    egin{align} label{equ555}
    mathrm{rank}\,J_{ell}(0)^{k-1}-mathrm{rank}\,J_{ell}(0)^k= egin{cases} 1 quad ext{如果} k leqslant ell \ 0 quad ext{如果} k>ell end{cases} ,quad k=1,2,cdots
    end{align}
    其中我们标注 (mathrm{rank}\,J_{ell}(0)^0=ell).

    现在我们设 (Ain M_n,lambda in mathbb{C}), 令 (k) 是一个正整数,设
    egin{align}
    r_k(A,lambda)=mathrm{rank} \, (A-lambda I)^k, quad r_0(A,lambda):=n
    end{align}
    又定义
    egin{align}
    w_k(A,lambda)=r_{k-1}(A,lambda)-r_k(A,lambda), quad w_1(A,lambda)=n-r_1(A,lambda)
    end{align}
    如果 (Ain M_n), 且 (lambda in mathbb{C}) 不是 (A) 的特征值,对所有 (k=1,2,cdots), 有 (w_k(A,lambda)=0). 考虑 Jordan 矩阵
    egin{align} label{equ888}
    J=J_3(0)oplus J_3(0)oplus J_2(0) oplus J_2(0) oplus J_2(0) oplus J_1(0)
    end{align}
    容易验证 (r_1(J,0)=7), (r_2(J,0)=2) 以及 (r_3(J,0)=r_4 (J,0)=0), (w_1(J,0)=6) 是阶至少为 1 的分块的个数, (w_2(J,0)=5) 是阶至少为 2 的分块的个数, (w_3(J,0)=2) 是阶至少为 3 的分块的个数, (w_4(J,0)=0) 是阶至少为 4 的分块的个数。又注意到 (w_1(J,0)-w_2(J,0)=1) 是阶为 1 的分块的个数, (w_2(J,0)-w_3(J,0)=3) 是阶为 2 的分块的个数, (w_3(J,0)-w_4(J,0)=2) 是阶为 3 的分块的个数,这不是偶然的。
    利用式 ef{equ444} 和 ef{equ555} 来说明 (w_k(A,lambda)) 的代数意义:

    egin{align}
    w_k (A,lambda) &= Big ( mathrm{rank}\,J_{m_1}(0)^{k-1}-mathrm{rank}\,J_{m_1}(0)^k Big ) + cdots + Big ( mathrm{rank}\,J_{m_p}(0)^{k-1}-mathrm{rank}\,J_{m_p}(0)^k Big) otag \
    &=(1, \, ext{如果}\, m_1 geqslant k)+ cdots +(1, \, ext{如果} \, m_p geqslant k) otag \ label{equ666}
    &= ext{特征值为}\, lambda \,\, ext{且阶至少为} \, k \,\, ext{的分块的个数}
    end{align}

    特别地,(w_1(A,lambda))(A) 的以 (lambda) 为特征值的所有各阶的 Jordan 块的个数,由定义知:(w_1(A,lambda)) 也是 (lambda) 作为 (A) 的特征值的几何重数.
    利用特征刻画 ef{equ666},我们看出 (w_k(A,lambda)-w_{k+1}(A,lambda)) 是以 (lambda) 为特征值而阶至少为 (k) 且不含阶至少为 (k+1) 的 Jordan 块的个数,这也就是以 (lambda) 为特征值且阶恰好为 (k) 的 Jordan 块的个数.
    由秩恒等式 ef{equ444} 知,用 (q) 表示 (Ain M_n) 的以 (lambda in mathbb{C}) 为特征值的最大的 Jordan 块的阶,则当 (k geqslant q) 时,都有 (mathrm{rank} \, (A-lambda I)^k=mathrm{rank} \, (A-lambda I)^q), 这个整数 (q) 称为 (lambda) 作为 (A) 的特征值的指数. (w_1(A,lambda)+w_2(A,lambda)+cdots+w_q(A,lambda))(lambda) 的代数重数.

    Weyr 特征

    (Ain M_n) 的与 (lambda in mathbb{C}) 相关的 Weyr 特征定义为
    egin{align}
    w(A,lambda)=(w_1(A,lambda),cdots,w_q(A,lambda))
    end{align}
    其中 (q)(lambda) 作为 (A) 的特征值的指数.
    如前所述,与 (A) 相似的 Jordan 矩阵 (J) 的构造完全由 (A) 的与不同的特征值相关的 Weyr 特征所决定的,这就意味着两个本质上不同的 Jordan 矩阵不可能都与 (A) 相似,因为它们的 Weyr 特征必定是不相同的. 我们可以说:两个复方阵 (A,Bin M_n) 相似,当且仅当它们特征值相同,且每个特征值相关的 Weyr 特征相同.
     

    Segre 特征

    (A) 的每个不同的特征值 (lambda),列出 (A) 的以 (lambda) 为特征值的所有 Jordan 块的阶,则给定的 (Ain M_n) 的 Jordan 块构造就可以完全确定. 将 (A) 的以 (lambda) 为特征值的 Jordan 块的阶按照非增次序排列而成的表
    egin{align}
    s_1(A,lambda) geqslant s_2(A,lambda) geqslant cdots geqslant s_{w_1(A,lambda)}(A,lambda) >0= s_{w_1(A,lambda)+1}(A,lambda) =cdots
    end{align}
    称为 (A) 的特征值 (lambda) 想关的 (Segre) 特征. 注意到 (s_1(A,lambda))(lambda) 作为 (A) 的特征值的指数,而 (s_{w_1}(A,lambda))(A) 的以 (lambda) 为特征值的最小的 Jordan 块的阶. 如矩阵 ef{equ888} 的与零这个特征值相关的 Segre 特征是 (3,3,2,2,2,1), (s_1(J,0)=3) 以及 (s_6(J,0)=1). 如果已知 Segre 特征,那么 Weyr 特征就容易得出,反之亦然.


    应该知道什么

    • 每个复矩阵都与一个本质上唯一的 Jordan 矩阵相似
    • Weyr 特征与 Segre 特征及其元素代表的意义
  • 相关阅读:
    C语言博客I作业04
    C语言I博客作业03
    C语言I博客作业02
    C语言ll作业01
    C语言寒假大作战04
    C语言寒假大作战03
    C语言寒假大作战02
    C语言寒假大作战01
    C语言I作业12—学期总结
    C语言I博客作业11
  • 原文地址:https://www.cnblogs.com/zhoukui/p/7751703.html
Copyright © 2020-2023  润新知