Jordan 标准型的实例
将学习到什么
练习一下如何把一个矩阵化为 Jordan 标准型.
将矩阵化为 Jordan 标准型需要三步:
- 第一步 求出矩阵 \(A \in M_n\) 全部的特征值 \(\lambda_1,\cdots,\lambda_t\), 假设有 \(t\) 个不同的特征值
- 第二步 Jordan 标准型定理 中的推论告诉我们:\(w_k(A,\lambda)-w_{k+1}(A,\lambda)\) 是以 \(\lambda\) 为特征值且阶恰好为 \(k\) 的 Jordan 块的个数. 我们就利用这个公式计算出以 \(\lambda\) 为特征值,阶为 \(\ell\) 的个数, \(\ell, \ell=1,2,\cdots\) 逐次计算. 以 \(\lambda\) 为特征值的 Jordan 块阶数之和等于特征值 \(\lambda\) 的代数重数,由此可知是否已经找出全部以 \(\lambda\) 为特征值的 Jordan 块
- 第三步 将所获得的 Jordan 块按任意次序排列成 Jordan 矩阵.
例 1
将矩阵
\begin{align}
A=\begin{bmatrix} 2 & 6 & -15 \\ 1 & 1 & -5 \\ 1 & 2 & -6 \end{bmatrix}
\end{align}
化为 Jordan 标准型.
第一步:求特征值
矩阵 \(A\) 的特征多项式为
\begin{align}
\lvert \lambda I-A \rvert =\begin{bmatrix} \lambda-2 & -6 & 15 \\ -1 & \lambda-1 & 5 \\ -1 & -2 & \lambda+6 \end{bmatrix} =(\lambda+1)^3
\end{align}
所以它只有一个特征值 \(\lambda_1=-1\), 代数重数为 3.
第二步:求 Jordan 块
对 \(\lambda_1=-1\), 令
\begin{align}
B=A-\lambda_1 I = A+I =\begin{bmatrix} 3 & 6 & -15 \\ 1 & 2 & -5 \\ 1 & 2 & -5 \end{bmatrix}, \qquad B^2=0
\end{align}
所以以 \(\lambda_1\) 为特征值阶为 1 的 Jordan 块的个数为
\begin{align}
w_1(A,\lambda_1)-w_2(A,\lambda_1)=[n-r_1(A,\lambda_1)] - [r_1(A,\lambda_1)-r_2(A,\lambda_1)] = [3-1]-[1-0]=1
\end{align}
其中 $r_k(A,\lambda)=\mathrm{rank} (A-\lambda I)^k, \quad r_0(A,\lambda):=n $, \(n\) 为方阵 \(A\) 的大小.
同理,以 \(\lambda_1\) 为特征值阶为 2 的 Jordan 块的个数为
\begin{align}
w_2(A,\lambda_1)-w_3(A,\lambda_1)=[r_1(A,\lambda_1)-r_2(A,\lambda_1)] - [r_2(A,\lambda_1)-r_3(A,\lambda_1)] = [1-0]-[0-0]=1
\end{align}
上面两个 Jordan 块阶数之和为 3,等于 \(\lambda_1\) 的重数,因而不再存在以 \(\lambda_1\) 为特征值的其它 Jordan 块. 因矩阵 \(A\) 没有其它特征值,故 Jordan 块求解完毕.
第三步:组成 Jordan 矩阵
只有一个重数为 3 的特征值 \(\lambda_1=-1\),一阶二阶各一个,所以矩阵 \(A\) 的 Jordan 标准型为
\begin{align}
J=\begin{bmatrix} -1 & 0 & 0 \\ 0 & -1 & 1 \\ 0 & 0 & -1\end{bmatrix}
\end{align}
例 2
将矩阵
\begin{align}
A=\begin{bmatrix} 3 & -4 & 0 & 2 \\ 4 & -5 & -2 & 4 \\ 0 & 0 & 3 & -2 \\ 0 & 0 & 2 & -1 \end{bmatrix}
\end{align}
化为 Jordan 标准型.
矩阵 \(A\) 的特征多项式为
\begin{align}
\lvert \lambda I-A \rvert =\begin{bmatrix} \lambda-3 & 4 & 0 & -2 \\ -4 & \lambda+5 & 2 & -4 \\ 0 & 0 & \lambda-3 & 2 \\ 0 & 0 & -2 & \lambda+1 \end{bmatrix}= (\lambda+1)^2(\lambda-1)^2
\end{align}
所以它有两个特征值 \(\lambda_1=-1\) 和 \(\lambda_2=1\), 代数重数都为 2.
第二步:求 Jordan 块
对 \(\lambda_1=-1\), 令
\begin{align}
B_1=A-\lambda_1 I = A+I &=\begin{bmatrix} 4 & -4 & 0 & 2 \\ 4 & -4 & -2 & 4 \\ 0 & 0 & 4 & -2 \\ 0 & 0 & 2 & 0 \end{bmatrix} \qquad \mathrm{rank}\,B_1=3 \\
B_1^2 &=\begin{bmatrix} 0 & 0 &12 & -8 \\ 0 & 0 & 8 & -4 \\ 0 & 0 & 12 & -8 \\ 0 & 0 & 8 & -4 \end{bmatrix} \qquad \mathrm{rank}\,B_1^2=2 \\
B_1^3 &=\begin{bmatrix} 0 & 0 & 32 & -24 \\ 0 & 0 & 24 & -16 \\ 0 & 0 & 32 & -24 \\ 0 & 0 & 24 & -16 \end{bmatrix} \qquad \mathrm{rank}\,B_1^3=2
\end{align}
所以以 \(\lambda_1\) 为特征值阶为 1 的 Jordan 块的个数为
\begin{align}
w_1(A,\lambda_1)-w_2(A,\lambda_1)=[n-r_1(A,\lambda_1)] - [r_1(A,\lambda_1)-r_2(A,\lambda_1)] = [4-3]-[3-2]=0
\end{align}
以 \(\lambda_1\) 为特征值阶为 2 的 Jordan 块的个数为
\begin{align}
w_2(A,\lambda_1)-w_3(A,\lambda_1)=[r_1(A,\lambda_1)-r_2(A,\lambda_1)] - [r_2(A,\lambda_1)-r_3(A,\lambda_1)] = [3-2]-[2-2]=1
\end{align}
上面第二个 Jordan 块阶数为 2,等于 \(\lambda_1\) 的重数,所以以 \(\lambda_1\) 为特征值的 Jordan 块求解完毕.
对 \(\lambda_2=1\), 令
\begin{align}
B_2=A-\lambda_2 I = A-I &=\begin{bmatrix} 2 & -4 & 0 & 2 \\ 4 & -6 & -2 & 4 \\ 0 & 0 & 2 & -2 \\ 0 & 0 & 2 & -2 \end{bmatrix} \qquad \mathrm{rank}\,B_2=3 \\
B_2^2 &=\begin{bmatrix} -12 & 16 & 12 & -16 \\ -16 & 20 & 16 & -20 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{bmatrix} \qquad \mathrm{rank}\,B_2^2=2 \\
B_2^3 &=\begin{bmatrix} 40 & -48 & -40 & 48 \\ 48 & -56 & -48 & 56 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{bmatrix} \qquad \mathrm{rank}\,B_2^3=2
\end{align}
所以以 \(\lambda_2\) 为特征值阶为 1 的 Jordan 块的个数为
\begin{align}
w_1(A,\lambda_2)-w_2(A,\lambda_2)=[n-r_1(A,\lambda_2)] - [r_1(A,\lambda_2)-r_2(A,\lambda_2)] = [4-3]-[3-2]=0
\end{align}
以 \(\lambda_2\) 为特征值阶为 2 的 Jordan 块的个数为
\begin{align}
w_2(A,\lambda_2)-w_3(A,\lambda_2)=[r_1(A,\lambda_2)-r_2(A,\lambda_2)] - [r_2(A,\lambda_2)-r_3(A,\lambda_2)] = [3-2]-[2-2]=1
\end{align}
上面第二个 Jordan 块阶数为 2,等于 \(\lambda_2\) 的重数,所以以 \(\lambda_2\) 为特征值的 Jordan 块求解完毕.
由 Jordan 标准型定理 的式 (6) 知,矩阵 \(B^k\), \(k>2\) 的秩不会再变化了,即次数大于特征值的最大的 Jordan 块的阶数时不再变化,最少为 \(n\) 减去 \(\lambda\) 的最大的 Jordan 块的阶数,这里也就是 2.
第三步:组成 Jordan 矩阵
以 \(\lambda_1=-1\) 和 \(\lambda_2=1\)为特征值的 Jordan 块各是一个二阶的,所以矩阵 \(A\) 的 Jordan 标准型为
\begin{align}
J=\begin{bmatrix} -1 & 1 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 0 & 1 & 1 \\ 0&0&0& 1\end{bmatrix}
\end{align}
Jordan 标准型的实例的更多相关文章
- 【线性代数】 06 - Jordan标准型
现在就来研究将空间分割为不变子空间的方法,最困难的是我们还不知道从哪里着手.你可能想到从循环子空间出发,一块一块地进行分割,但这个方案的存在性和唯一性都不能解决.不变子空间分割不仅要求每个子空间\(V ...
- Jordan 标准型的推论
将学习到什么 从 Jordan 标准型出发,能够获得非常有用的信息. Jordan 矩阵的构造 Jordan 矩阵 \begin{align} J=\begin{bmatrix} J_{n_1}( ...
- Jordan 标准型定理
将学习到什么 就算两个矩阵有相同的特征多项式,它们也有可能不相似,那么如何判断两个矩阵是相似的?答案是它们有一样的 Jordan 标准型. Jordan 标准型定理 这节目的:证明每个复矩阵都与一 ...
- 线性代数 | Jordan 标准型的笔记
内容概述: 把方阵 A 的特征多项式 \(c(λ)=|λE-A|\) 展开成 \(c(λ)=\sum_ia_i\lambda^i\) 的形式,然后使用神乎其技的证明,得到 \(c(A)=O\),特征多 ...
- [转载] $\mathrm{Jordan}$标准型的介绍
本文转载自陈洪葛的博客$,$ 而实际上来自xida博客朝花夕拾$,$ 可惜该博客已经失效 $\mathrm{Jordan}$ 标准形定理是线性代数中的基本定理$,$ 专门为它写一篇长文好像有点多余$: ...
- [Bilingual] Different proofs of Jordan cardinal form (Jordan标准型的几种证明)
- 实 Jordan 标准型和实 Weyr 标准型
将学习到什么 本节讨论关于实矩阵的实形式的 Jordan 标准型,也讨论关于复矩阵的另外一种形式的 Jordan 标准型,因为它在与交换性有关的问题中很有用. 实 Jordan 标准型 假设 \( ...
- Jordan 块的几何
设 $V$ 是复数域 $\mathbb{C}$ 上的 $n$ 维线性空间, $\varphi$ 是 $V$ 上的线性变换, $A\in M_n(\mathbb{C})$ 是 $\varphi$ 在某组 ...
- Jordan标准形
一.引入 前面已经指出,一切n阶矩阵A可以分成许多相似类.今要在与A相似的全体矩阵中,找出一个较简单的矩阵来作为相似类的标准形.当然以对角矩阵作为标准形最好,可惜不是每一个矩阵都能与对角矩阵相似.因此 ...
随机推荐
- HDU - 1715 - 大菲波数 - JAVA
http://acm.hdu.edu.cn/showproblem.php?pid=1715 import java.io.*; import java.util.*; import java.mat ...
- 百度地图API示例 JS
http://developer.baidu.com/map/jsdemo.htm#c2_2
- Window 10 Alt Tap 切换窗口,设置成Windows 7风格
升级了Windows 10 以后,切换窗口非常的难受.新版的窗口切换的图标变成了窗口预览了,这总让我一瞬间找不到要去的窗口,所以我打算切换成Windows 7 的小图标风格. Windows + R ...
- idea svn 问题
https://blog.csdn.net/liyantianmin/article/details/52837506
- c++中初始化列表的初始化变量顺序问题
例题来看:请问下面程序打印出的结果是什么? #include <iostream> #include <string> using namespace std; class b ...
- (转) Git
http://www.liaoxuefeng.com/wiki/0013739516305929606dd18361248578c67b8067c8c017b000
- 51Nod 1272 最大距离 (栈或贪心)
#include <cstdio> #include <queue> #include <cstring> #include <iostream> #i ...
- .Net Core应用框架Util介绍(二) 转
Util的开源地址 https://github.com/dotnetcore/util Util的开源协议 Util以MIT协议开源,这是目前最宽松的开源协议,你不仅可以用于商业项目,还能把Util ...
- Serervlet | 慕课课程实战 | 编写登录逻辑
Users.java package com.po; public class Users { private String username; private String password; pu ...
- [51Nod1952] 栈
Description 不支持后端删除的dequeue,每次操作后查询最大值. \(n\leq10^7\).时限1.5s,不用考虑读入/输出复杂度. Solution 首先考虑如果没有后端删除怎么做, ...