第 28 章

相似矩阵与 Jordan 形:同一特征值下的不同家族

相似矩阵共享全部特征值,但不一定共享特征向量;重特征值会把矩阵分成多个相似类,Jordan 形为每一类选出最接近对角的标准代表。

Similar Matrices and Jordan Form: One Spectrum, Many Families

Similar matrices share all eigenvalues but not necessarily eigenvectors; repeated eigenvalues split matrices into separate similarity classes, and Jordan form picks the canonical representative of each.

首发于 2019.09.25中 / EN
系统课程linear-algebraeigenvaluesjordan-form
语言
字号
目录
篇内目录

开场:正定矩阵的收尾

这一讲先用十分钟补完正定矩阵,再进入主题。补的部分回答三个自然的问题。

正定矩阵的逆还是正定吗

是。理由用特征值最快:若 AA 的特征值是 λ1,…,λn>0\lambda_1,\dots,\lambda_n>0,则 A−1A^{-1} 的特征值是 1/λ1,…,1/λn1/\lambda_1,\dots,1/\lambda_n,仍然全为正。逆矩阵对称,因此也正定。

两个正定矩阵之和还是正定吗

是。不需要知道 A+BA+B 的特征值或主元,直接看二次型:

xT(A+B)x=xTAx+xTBx.x^{\mathsf T}(A+B)x=x^{\mathsf T}Ax+x^{\mathsf T}Bx.

右边是两个正数之和,对任何 x≠0x\neq 0 都为正。

这条证明也说明了一件事:遇到正定性问题,先试二次型,它往往绕开了特征值和主元。

ATAA^{\mathsf T}A 什么时候正定

设 AA 是 m×nm\times n 的长方矩阵。它本身通常不是方阵,谈不上正定;但 ATAA^{\mathsf T}A 是 n×nn\times n 对称矩阵,正是最小二乘法里的核心组合。

xTATAx=(Ax)T(Ax)=∥Ax∥2≥0.x^{\mathsf T}A^{\mathsf T}Ax=(Ax)^{\mathsf T}(Ax)=\lVert Ax\rVert^2\ge 0.

所以 ATAA^{\mathsf T}A 永远半正定。它正定的充要条件是 ∥Ax∥2>0\lVert Ax\rVert^2>0 对所有 x≠0x\neq 0 成立,也就是 Ax≠0Ax\neq 0,即 AA 的零空间只有零向量。

用秩来表述:

ATA 正定  ⟺  A 的列线性无关  ⟺  rank⁡A=n.A^{\mathsf T}A\text{ 正定}\iff A\text{ 的列线性无关}\iff \operatorname{rank}A=n.

这也解释了为什么最小二乘的正规方程在列无关时可以放心求解:ATAA^{\mathsf T}A 不仅可逆,而且正定。

补充:正定矩阵在数值计算中还有一个重要性质——消元时不需要换行。主元不会变成零或极小值,因此 A=LUA=LU 的分解稳定。这使正定矩阵成为数值线性代数里最「好算」的一类矩阵。

时间轴

时间内容
0:08相似矩阵的主题
1:23收尾:ATAA^{\mathsf T}A 的正定性
4:41长方矩阵与秩
6:15ATAA^{\mathsf T}A 正定当且仅当列向量独立
11:38预告 SVD
13:00相似的定义 B=M−1AMB=M^{-1}AM
14:47例子:AA 与对角矩阵 Λ\Lambda 相似
17:48构造另一个与 AA 相似的矩阵
22:27相似矩阵共享特征值:完整证明
25:34对角矩阵的相似类有多大
26:50坏矩阵:只有一个线性无关的特征向量
30:32不可对角化的例子
33:27Jordan 块的形状
37:35两个 Jordan 块与特征向量的不足
40:36特征值相同但矩阵不相似
43:11Jordan 形式给出相似类的完整分类

相似矩阵的定义

设 A,BA,B 都是 n×nn\times n 方阵。若存在可逆矩阵 MM,使得

B=M−1AM,B=M^{-1}AM,

则称 AA 与 BB 相似。

这个定义看起来像一次「换基」:同一个线性变换,在不同的基下写成不同的矩阵。MM 就是从旧基到新基的过渡矩阵。

最重要的相似:对角化

若 AA 有 nn 个线性无关的特征向量,把它们作为列组成 SS,则

S−1AS=Λ.S^{-1}AS=\Lambda.

所以 AA 相似于对角矩阵 Λ\Lambda。在所有与 AA 相似的矩阵里,Λ\Lambda 是形式最简单的一个。

但 SS 不是唯一的选择。任取一个可逆的 MM,M−1AMM^{-1}AM 都会得到另一个与 AA 相似的矩阵,只是通常不再是对角阵。

相似矩阵共享特征值

这是本讲的主定理:

相似矩阵有相同的特征值。

证明

设 Ax=λxAx=\lambda x,且 B=M−1AMB=M^{-1}AM。在 AA 与 xx 之间插入 MM−1MM^{-1}:

A(MM−1)x=λx.A(MM^{-1})x=\lambda x.

左边重新结合为 (AM)(M−1x)(AM)(M^{-1}x),再在两边左乘 M−1M^{-1}:

M−1AM (M−1x)=λ (M−1x).M^{-1}AM\,(M^{-1}x)=\lambda\,(M^{-1}x).

即

B(M−1x)=λ(M−1x).B(M^{-1}x)=\lambda(M^{-1}x).

所以 λ\lambda 也是 BB 的特征值,对应的特征向量是 M−1xM^{-1}x。

特征值不变,特征向量被 M−1M^{-1} 搬到了新基下。 这也符合直觉:换基不改变线性变换本身,只改变它的坐标表示。

一个可手算的例子

取

A=[2112],M=[1101],M−1=[1−101].A=\begin{bmatrix}2&1\\1&2\end{bmatrix},\qquad M=\begin{bmatrix}1&1\\0&1\end{bmatrix},\qquad M^{-1}=\begin{bmatrix}1&-1\\0&1\end{bmatrix}.

先算 AMAM:

AM=[2112][1101]=[2313].AM=\begin{bmatrix}2&1\\1&2\end{bmatrix} \begin{bmatrix}1&1\\0&1\end{bmatrix} =\begin{bmatrix}2&3\\1&3\end{bmatrix}.

再左乘 M−1M^{-1}:

B=M−1AM=[1−101][2313]=[1013].B=M^{-1}AM= \begin{bmatrix}1&-1\\0&1\end{bmatrix} \begin{bmatrix}2&3\\1&3\end{bmatrix} =\begin{bmatrix}1&0\\1&3\end{bmatrix}.

AA 的特征值是 33 与 11。BB 的迹是 44,行列式是 33,特征值同样是 33 与 11。

BB 看起来和 AA 很不一样,但它们属于同一个相似类。

相似类的图景

固定特征值 33 与 11,所有具有这两个特征值的 2×22\times2 矩阵构成一个家族,家族成员两两相似。家族的「代表」是对角矩阵

[3001].\begin{bmatrix}3&0\\0&1\end{bmatrix}.

除此之外还有上三角、下三角、各种非对称形式——只要特征值仍是 33 与 11,就都能通过某个 MM 与 AA 相连。

坏情形:特征值重复

若 λ1=λ2\lambda_1=\lambda_2,事情变得微妙:矩阵可能没有足够多的特征向量,因而不可对角化。

一个孤独的矩阵

考虑 4I=[4004]4I=\begin{bmatrix}4&0\\0&4\end{bmatrix}。对任何可逆 MM,

M−1(4I)M=4M−1M=4I.M^{-1}(4I)M=4M^{-1}M=4I.

它只与自己相似。在所有特征值为 4,44,4 的矩阵中,4I4I 自成一类。

另一个家族

矩阵

[4104]\begin{bmatrix}4&1\\0&4\end{bmatrix}

的特征值也是 4,44,4,但它只有一个线性无关的特征向量:零空间由 [10]\begin{bmatrix}1\\0\end{bmatrix} 张成,维数为 11。它不可对角化,因此不可能与 4I4I 相似。

所有其余特征值为 4,44,4 的矩阵都落在同一个相似类里,而 [4104]\begin{bmatrix}4&1\\0&4\end{bmatrix} 是这个类里最接近对角的形式。

更多同族成员

只要保持迹为 88、行列式为 1616,就得到同族的矩阵,例如

[51−13],[41704].\begin{bmatrix}5&1\\-1&3\end{bmatrix},\qquad \begin{bmatrix}4&17\\0&4\end{bmatrix}.

它们都只有一个特征向量,彼此相似。

数特征向量还不够

一个自然的猜想是:特征值相同、特征向量个数相同,就应该相似。这个猜想是错的。

考虑两个 4×44\times4 矩阵,特征值全是 00。

第一个:两个独立的 2×22\times2 块

J1=[0100000000010000].J_1=\begin{bmatrix} 0&1&0&0\\ 0&0&0&0\\ 0&0&0&1\\ 0&0&0&0 \end{bmatrix}.

它的秩为 22,零空间维数为 4−2=24-2=2,有两个特征向量。

第二个:一个 3×33\times3 块加一个 1×11\times1 块

J2=[0100001000000000].J_2=\begin{bmatrix} 0&1&0&0\\ 0&0&1&0\\ 0&0&0&0\\ 0&0&0&0 \end{bmatrix}.

它的秩同样是 22,零空间维数同样是 22,也有两个特征向量。

但 J1J_1 与 J2J_2 不相似。区别不在数量,而在块的大小:J1J_1 是 2+22+2,J2J_2 是 3+13+1。这说明要完整刻画相似类,必须知道特征向量的「归属结构」,而不只是总数。

Jordan 块与 Jordan 定理

Jordan 块

一个 Jordan 块 JiJ_i 形如

Ji=[λi1λi⋱⋱1λi],J_i=\begin{bmatrix} \lambda_i&1&&\\ &\lambda_i&\ddots&\\ &&\ddots&1\\ &&&\lambda_i \end{bmatrix},

对角线上是重复的特征值,超对角线上全是 11,其余为零。

关键性质:一个 Jordan 块只有一个特征向量。块越大,缺失的特征向量越多。

Jordan 定理

每个方阵 AA 都相似于某个 Jordan 矩阵 JJ,JJ 由若干 Jordan 块沿对角线拼成。

J=[J1⋱Jd].J=\begin{bmatrix} J_1&&\\ &\ddots&\\ &&J_d \end{bmatrix}.

两个要点:

  1. 块的数量等于线性无关特征向量的个数,因为每块贡献一个。
  2. 块的大小记录了特征向量的缺失方式,这正是相似类之间的区别。

若所有特征值互不相同,每个块都是 1×11\times1,JJ 就是对角矩阵,d=nd=n。所以对角化是 Jordan 形的特例。

为什么现在不强调 Jordan 形

Jordan 形在代数上很漂亮:它给出了相似类的完整分类。但它对数值极其敏感——只要特征值有微小扰动,重根就会分开,块结构随之改变,秩也会跳变。因此实际计算中很少真的去求 Jordan 形。

课堂的判断是:Jordan 形是理解相似类的正确语言,但不再是这门课的高潮;真正的主角是下一讲的奇异值分解。

本讲结论

  • 正定矩阵的逆、和仍然是正定;ATAA^{\mathsf T}A 恒为半正定,列无关时为正定。
  • 相似的定义是 B=M−1AMB=M^{-1}AM,等价于同一个线性变换换一组基。
  • 相似矩阵共享全部特征值,特征向量按 M−1M^{-1} 变换。
  • 特征值重复时,相似类会分裂;4I4I 自成一类,其余矩阵构成另一类。
  • 仅统计特征向量个数不足以判断相似,还要看 Jordan 块的结构。
  • 每个方阵都相似于一个 Jordan 矩阵;特征值互异时它退化为对角阵。

理解检验

1. 为什么 ATAA^{\mathsf T}A 总是半正定?它什么时候正定?

参考答案:对任意 xx,

xTATAx=(Ax)T(Ax)=∥Ax∥2≥0,x^{\mathsf T}A^{\mathsf T}Ax=(Ax)^{\mathsf T}(Ax)=\lVert Ax\rVert^2\ge0,

所以二次型从不为负,ATAA^{\mathsf T}A 半正定(它本身对称:(ATA)T=ATA(A^{\mathsf T}A)^{\mathsf T}=A^{\mathsf T}A)。要正定,需要 ∥Ax∥2>0\lVert Ax\rVert^2>0 对所有 x≠0x\ne0 成立,即 Ax≠0Ax\ne0 对所有非零 xx 成立,也就是 N(A)={0}N(A)=\{\mathbf{0}\}。这等价于 AA 的列线性无关,即 rank⁡A=n\operatorname{rank}A=n。列相关的方向会让 Ax=0Ax=0,二次型在该方向取零,此时只是半正定。

2. 证明相似矩阵有相同的特征值,并写出特征向量之间的对应关系。

参考答案:设 B=M−1AMB=M^{-1}AM,且 Ax=λxAx=\lambda x,x≠0x\ne0。在 AA 与 xx 之间插入 MM−1MM^{-1}:

A(MM−1)x=λx  ⟹  (AM)(M−1x)=λx.A(MM^{-1})x=\lambda x \;\Longrightarrow\; (AM)(M^{-1}x)=\lambda x .

两边左乘 M−1M^{-1}:

M−1AM (M−1x)=λ (M−1x),M^{-1}AM\,(M^{-1}x)=\lambda\,(M^{-1}x),

即 B(M−1x)=λ(M−1x)B(M^{-1}x)=\lambda(M^{-1}x)。因为 M−1M^{-1} 可逆且 x≠0x\ne0,所以 M−1x≠0M^{-1}x\ne0,它是 BB 的属于同一特征值 λ\lambda 的特征向量。因此特征值集合相同,特征向量按 x↦M−1xx\mapsto M^{-1}x 搬到新基下。

3. 为什么 4I4I 只与自己相似?

参考答案:对任何可逆 MM,

M−1(4I)M=4M−1M=4I.M^{-1}(4I)M=4M^{-1}M=4I .

所以与 4I4I 相似的所有矩阵都只能是 4I4I 本身,它的相似类只有一个成员。直观上,4I4I 在每个方向上都把向量拉伸 44 倍,换任何基都不改变这一表示;它的特征空间是整个平面,没有「特殊方向」可被换基改变。

4. 两个 4×44\times4 矩阵都有两个特征向量,为什么仍可能不相似?

参考答案:特征向量个数相同只是零空间维数相同,不足以刻画 Jordan 块的结构。正文的两个矩阵

J1=[0100000000010000],J2=[0100001000000000]J_1=\begin{bmatrix}0&1&0&0\\0&0&0&0\\0&0&0&1\\0&0&0&0\end{bmatrix}, \qquad J_2=\begin{bmatrix}0&1&0&0\\0&0&1&0\\0&0&0&0\\0&0&0&0\end{bmatrix}

秩都为 22,零空间维数都为 22,各有两列特征向量;但 J1J_1 是两块 2×22\times2,J2J_2 是一块 3×33\times3 加一块 1×11\times1。区分方法是看幂的秩:J12=0J_1^2=0,秩为 00;而 J22J_2^2 在第 (1,3)(1,3) 位仍为 11,秩为 11。相似变换保持 rank⁡((A−λI)k)\operatorname{rank}((A-\lambda I)^k),所以两者不可能相似。

5. 特征值互不相同时,Jordan 形是什么样子?为什么?

参考答案:此时 Jordan 形就是对角矩阵,每个块都是 1×11\times1,共有 nn 个块。原因是:不同特征值对应的特征向量线性无关,nn 个互不相同的特征值给出 nn 个线性无关的特征向量,矩阵可对角化;而 Jordan 块的大小由特征向量的缺失决定,没有缺失就没有大于 11 的块。于是 J=ΛJ=\Lambda,对角化只是 Jordan 形的特例。

与前后课程的关系

上一讲用二次型与海森矩阵刻画了正定矩阵;这一讲先补完它的运算性质,再进入相似与 Jordan 形。下一讲转向奇异值分解:它不需要矩阵对称,也不需要可对角化,却对所有矩阵都成立。

Finishing positive definiteness

The lecture spends its first ten minutes closing three loose ends.

Is the inverse of a positive definite matrix positive definite?

Yes. Eigenvalues give the shortest proof: if AA has eigenvalues λ1,…,λn>0\lambda_1,\dots,\lambda_n>0, then A−1A^{-1} has eigenvalues 1/λ1,…,1/λn1/\lambda_1,\dots,1/\lambda_n, still positive. The inverse is symmetric, so it is positive definite.

Is the sum of two positive definite matrices positive definite?

Yes. There is no need to know the eigenvalues or pivots of A+BA+B; use the quadratic form:

xT(A+B)x=xTAx+xTBx.x^{\mathsf T}(A+B)x=x^{\mathsf T}Ax+x^{\mathsf T}Bx.

The right side is a sum of two positive numbers for every x≠0x\neq 0.

The proof also carries a lesson: when a positive-definiteness question appears, try the quadratic form first. It often bypasses eigenvalues and pivots entirely.

When is ATAA^{\mathsf T}A positive definite?

Let AA be m×nm\times n. It is usually rectangular, so it is not itself positive definite. But ATAA^{\mathsf T}A is n×nn\times n and symmetric — the central object of least squares.

xTATAx=(Ax)T(Ax)=∥Ax∥2≥0.x^{\mathsf T}A^{\mathsf T}Ax=(Ax)^{\mathsf T}(Ax)=\lVert Ax\rVert^2\ge 0.

So ATAA^{\mathsf T}A is always positive semidefinite. It is positive definite exactly when ∥Ax∥2>0\lVert Ax\rVert^2>0 for every x≠0x\neq 0, that is, when Ax≠0Ax\neq 0, meaning the null space of AA is trivial.

In rank terms:

ATA positive definite  ⟺  A has independent columns  ⟺  rank⁡A=n.A^{\mathsf T}A\text{ positive definite}\iff A\text{ has independent columns}\iff \operatorname{rank}A=n.

This is why the normal equations are safe to solve when the columns are independent: ATAA^{\mathsf T}A is not merely invertible but positive definite.

Background: positive definite matrices also need no row exchanges during elimination. Pivots never become zero or dangerously small, so A=LUA=LU is numerically stable. It is the best-behaved class in numerical linear algebra.

Timeline

TimeTopic
0:08The theme: similar matrices
1:23Wrapping up: positive definiteness of ATAA^{\mathsf T}A
4:41Rectangular matrices and rank
6:15ATAA^{\mathsf T}A is positive definite iff the columns are independent
11:38A preview of the SVD
13:00Definition of similarity, B=M−1AMB=M^{-1}AM
14:47Example: AA similar to the diagonal matrix Λ\Lambda
17:48Constructing another matrix similar to AA
22:27Similar matrices share eigenvalues: the proof
25:34How large a similarity class can be
26:50A bad matrix with one independent eigenvector
30:32An example that cannot be diagonalized
33:27The shape of a Jordan block
37:35Two Jordan blocks and the eigenvector shortage
40:36Same eigenvalues yet not similar
43:11The Jordan form classifies similarity completely

The definition of similarity

Let AA and BB be n×nn\times n. They are similar if there is an invertible MM with

B=M−1AM.B=M^{-1}AM.

The definition is a change of basis: the same linear transformation written in a different basis, where MM is the change-of-basis matrix.

The most important similarity: diagonalization

If AA has nn independent eigenvectors, put them in the columns of SS. Then

S−1AS=Λ.S^{-1}AS=\Lambda.

So AA is similar to the diagonal matrix Λ\Lambda, the simplest member of its family.

But SS is not the only choice. Any invertible MM produces another matrix M−1AMM^{-1}AM similar to AA, usually no longer diagonal.

Similar matrices share eigenvalues

The main theorem:

Similar matrices have the same eigenvalues.

Proof

Suppose Ax=λxAx=\lambda x and B=M−1AMB=M^{-1}AM. Insert MM−1MM^{-1} between AA and xx:

A(MM−1)x=λx.A(MM^{-1})x=\lambda x.

Regroup the left side as (AM)(M−1x)(AM)(M^{-1}x) and multiply both sides on the left by M−1M^{-1}:

M−1AM (M−1x)=λ (M−1x).M^{-1}AM\,(M^{-1}x)=\lambda\,(M^{-1}x).

That is,

B(M−1x)=λ(M−1x).B(M^{-1}x)=\lambda(M^{-1}x).

So λ\lambda is also an eigenvalue of BB, with eigenvector M−1xM^{-1}x.

The eigenvalues stay fixed; the eigenvectors are moved by M−1M^{-1}. This matches intuition: changing basis does not change the transformation, only its coordinate representation.

A hand-checkable example

Take

A=[2112],M=[1101],M−1=[1−101].A=\begin{bmatrix}2&1\\1&2\end{bmatrix},\qquad M=\begin{bmatrix}1&1\\0&1\end{bmatrix},\qquad M^{-1}=\begin{bmatrix}1&-1\\0&1\end{bmatrix}.

First compute AMAM:

AM=[2112][1101]=[2313].AM=\begin{bmatrix}2&1\\1&2\end{bmatrix} \begin{bmatrix}1&1\\0&1\end{bmatrix} =\begin{bmatrix}2&3\\1&3\end{bmatrix}.

Then multiply on the left by M−1M^{-1}:

B=M−1AM=[1−101][2313]=[1013].B=M^{-1}AM= \begin{bmatrix}1&-1\\0&1\end{bmatrix} \begin{bmatrix}2&3\\1&3\end{bmatrix} =\begin{bmatrix}1&0\\1&3\end{bmatrix}.

The eigenvalues of AA are 33 and 11. The trace of BB is 44 and its determinant is 33, so its eigenvalues are again 33 and 11.

BB looks nothing like AA, yet they belong to the same similarity class.

The picture of a similarity class

Fix the eigenvalues 33 and 11. Every 2×22\times2 matrix with these eigenvalues forms one family, and any two members are similar. The representative is the diagonal matrix

[3001].\begin{bmatrix}3&0\\0&1\end{bmatrix}.

Upper triangular, lower triangular, and asymmetric forms all join the family as long as the spectrum stays {3,1}\{3,1\}.

The bad case: repeated eigenvalues

When λ1=λ2\lambda_1=\lambda_2, something subtler happens: the matrix may not have enough eigenvectors to diagonalize.

A matrix alone in its class

Consider 4I=[4004]4I=\begin{bmatrix}4&0\\0&4\end{bmatrix}. For any invertible MM,

M−1(4I)M=4M−1M=4I.M^{-1}(4I)M=4M^{-1}M=4I.

It is similar only to itself. Among all matrices with eigenvalues 4,44,4, it forms a class of one.

The other family

The matrix

[4104]\begin{bmatrix}4&1\\0&4\end{bmatrix}

also has eigenvalues 4,44,4, but only one independent eigenvector: the null space of A−4IA-4I is spanned by [10]\begin{bmatrix}1\\0\end{bmatrix}, dimension one. It is not diagonalizable, so it cannot be similar to 4I4I.

Every other matrix with eigenvalues 4,44,4 lies in one big class, and [4104]\begin{bmatrix}4&1\\0&4\end{bmatrix} is its closest-to-diagonal representative.

More members

Keep the trace at 88 and the determinant at 1616, and you stay in the same class, for example

[51−13],[41704].\begin{bmatrix}5&1\\-1&3\end{bmatrix},\qquad \begin{bmatrix}4&17\\0&4\end{bmatrix}.

Both have a single eigenvector and are similar to each other.

Counting eigenvectors is not enough

A natural guess is that equal eigenvalues plus equal eigenvector counts imply similarity. The guess is false.

Take two 4×44\times4 matrices whose eigenvalues are all zero.

First: two independent 2×22\times2 blocks,

J1=[0100000000010000].J_1=\begin{bmatrix} 0&1&0&0\\ 0&0&0&0\\ 0&0&0&1\\ 0&0&0&0 \end{bmatrix}.

Its rank is 22, so its nullity is 4−2=24-2=2: two eigenvectors.

Second: one 3×33\times3 block and one 1×11\times1 block,

J2=[0100001000000000].J_2=\begin{bmatrix} 0&1&0&0\\ 0&0&1&0\\ 0&0&0&0\\ 0&0&0&0 \end{bmatrix}.

Its rank is again 22, its nullity again 22, and it also has two eigenvectors.

Yet J1J_1 and J2J_2 are not similar. The difference is not the count but the block structure: J1J_1 is 2+22+2, J2J_2 is 3+13+1. A complete description of a similarity class needs to know how the eigenvectors are grouped, not just how many there are.

Jordan blocks and Jordan's theorem

The Jordan block

A Jordan block JiJ_i has the form

Ji=[λi1λi⋱⋱1λi],J_i=\begin{bmatrix} \lambda_i&1&&\\ &\lambda_i&\ddots&\\ &&\ddots&1\\ &&&\lambda_i \end{bmatrix},

with a repeated eigenvalue on the diagonal, ones on the superdiagonal, and zeros elsewhere.

The key property: a Jordan block has exactly one eigenvector. The larger the block, the more eigenvectors are missing.

The theorem

Every square matrix AA is similar to a Jordan matrix JJ, built from Jordan blocks along the diagonal.

J=[J1⋱Jd].J=\begin{bmatrix} J_1&&\\ &\ddots&\\ &&J_d \end{bmatrix}.

Two consequences:

  1. The number of blocks equals the number of independent eigenvectors, since each block contributes one.
  2. The block sizes record how eigenvectors are missing, and that is exactly what distinguishes similarity classes.

If all eigenvalues are distinct, every block is 1×11\times1, so JJ is diagonal and d=nd=n. Diagonalization is the special case.

Why Jordan form is no longer the climax

Jordan form is algebraically beautiful: it classifies similarity completely. But it is numerically fragile. A tiny perturbation splits repeated eigenvalues, changes the rank, and alters the block structure.

The lecture's verdict: Jordan form is the right language for understanding similarity classes, but it is no longer the centerpiece of the course. The next lecture's singular value decomposition is.

Takeaways

  • The inverse and the sum of positive definite matrices are positive definite; ATAA^{\mathsf T}A is always semidefinite and definite exactly when AA has independent columns.
  • Similarity means B=M−1AMB=M^{-1}AM, a change of basis for one linear transformation.
  • Similar matrices share all eigenvalues; eigenvectors transform by M−1M^{-1}.
  • Repeated eigenvalues split the spectrum into several similarity classes; 4I4I is a class of one.
  • Counting eigenvectors does not determine similarity; Jordan block structure does.
  • Every square matrix is similar to a Jordan matrix, which reduces to a diagonal matrix when eigenvalues are distinct.

Check your understanding

1. Why is ATAA^{\mathsf T}A always positive semidefinite, and when is it positive definite?

Reference answer: For every xx,

xTATAx=(Ax)T(Ax)=∥Ax∥2≥0,x^{\mathsf T}A^{\mathsf T}Ax=(Ax)^{\mathsf T}(Ax)=\lVert Ax\rVert^2\ge0,

so the quadratic form is never negative and ATAA^{\mathsf T}A is positive semidefinite (it is also symmetric, since (ATA)T=ATA(A^{\mathsf T}A)^{\mathsf T}=A^{\mathsf T}A). It is positive definite exactly when ∥Ax∥2>0\lVert Ax\rVert^2>0 for all x≠0x\ne0, i.e. Ax≠0Ax\ne0 for every nonzero xx, which says N(A)={0}N(A)=\{\mathbf{0}\}. Equivalently, the columns of AA are linearly independent, rank⁡A=n\operatorname{rank}A=n. A dependent direction makes Ax=0Ax=0 and the quadratic form vanish there, leaving only semidefiniteness.

2. Prove that similar matrices have the same eigenvalues and state how their eigenvectors correspond.

Reference answer: Let B=M−1AMB=M^{-1}AM and Ax=λxAx=\lambda x with x≠0x\ne0. Insert MM−1MM^{-1} between AA and xx:

A(MM−1)x=λx  ⟹  (AM)(M−1x)=λx.A(MM^{-1})x=\lambda x \;\Longrightarrow\; (AM)(M^{-1}x)=\lambda x .

Left-multiplying by M−1M^{-1},

M−1AM (M−1x)=λ (M−1x),M^{-1}AM\,(M^{-1}x)=\lambda\,(M^{-1}x),

i.e. B(M−1x)=λ(M−1x)B(M^{-1}x)=\lambda(M^{-1}x). Since M−1M^{-1} is invertible and x≠0x\ne0, we have M−1x≠0M^{-1}x\ne0, so it is an eigenvector of BB for the same λ\lambda. The eigenvalue set is unchanged and eigenvectors move by x↦M−1xx\mapsto M^{-1}x in the new basis.

3. Why is 4I4I similar only to itself?

Reference answer: For any invertible MM,

M−1(4I)M=4M−1M=4I.M^{-1}(4I)M=4M^{-1}M=4I .

So every matrix similar to 4I4I equals 4I4I itself; its similarity class has a single member. Intuitively, 4I4I stretches every direction by the same factor 44, and no change of basis alters that representation — its eigenspace is the whole space, so there is no distinguished direction for a change of basis to move.

4. Two 4×44\times4 matrices both have two eigenvectors. Why can they still fail to be similar?

Reference answer: Equal eigenvector counts only mean equal nullity; they do not capture the Jordan block structure. The two matrices from the text,

J1=[0100000000010000],J2=[0100001000000000],J_1=\begin{bmatrix}0&1&0&0\\0&0&0&0\\0&0&0&1\\0&0&0&0\end{bmatrix}, \qquad J_2=\begin{bmatrix}0&1&0&0\\0&0&1&0\\0&0&0&0\\0&0&0&0\end{bmatrix},

both have rank 22 and nullity 22, hence two eigenvectors each. But J1J_1 is two 2×22\times2 blocks while J2J_2 is one 3×33\times3 plus one 1×11\times1. The ranks of powers distinguish them: J12=0J_1^2=0 has rank 00, whereas J22J_2^2 has a 11 in position (1,3)(1,3) and rank 11. Similarity preserves rank⁡((A−λI)k)\operatorname{rank}((A-\lambda I)^k), so the two cannot be similar.

5. What does the Jordan form look like when all eigenvalues are distinct, and why?

Reference answer: It is diagonal, with every block 1×11\times1 — that is, nn blocks. Distinct eigenvalues have linearly independent eigenvectors, so nn distinct eigenvalues supply nn independent eigenvectors and the matrix is diagonalizable. Jordan block sizes measure how eigenvectors are missing; with no deficiency, no block is larger than 1×11\times1. Thus J=ΛJ=\Lambda, and diagonalization is the special case of Jordan form.

How this connects to the course

The previous lecture characterized positive definiteness through quadratic forms and Hessians; this one finishes its algebraic properties, then introduces similarity and Jordan form. The next lecture turns to the singular value decomposition, which requires neither symmetry nor diagonalizability and applies to every matrix.

来源