本节将讨论:当方阵 A 满足什么条件时,能找到可逆的矩阵 P,使得 B = P⁻1AP 为对角矩阵。能表示成 B = P⁻1AP 的矩阵 B 称为 A 的相似矩阵。
一、相似矩阵
Section titled “一、相似矩阵”设 A, B 为 n 阶方阵。如果存在某个 n 阶可逆矩阵 P,使得
B = P⁻1AP,则称 B 是 A 的相似矩阵,或称 A 与 B 相似,记为 A ∼ B。
方阵之间的“相似”关系具有以下三条性质,因此是一种等价关系:
- 反身性:A ∼ A。即任意一个方阵都与自己相似。这是因为等式 A = E⁻1AE 显然成立;
- 对称性:若 A ∼ B,则 B ∼ A。事实上,如果有 B = P⁻1AP,则显然有 A = PBP⁻1 = Q⁻1BQ,其中 Q = P⁻1;
- 传递性:若 A ∼ B, B ∼ C,则 A ∼ C。事实上,如果有 B = P⁻1AP, C = Q⁻1BQ,则 C = Q⁻1BQ = Q⁻1P⁻1APQ = R⁻1AR,其中 R = PQ。
相似矩阵还有下面的性质:
如果 A ∼ B,则 A 与 B 有相同的特征多项式,因而有相同的特征值与行列式。
证
已知存在可逆矩阵 P,使得 B = P⁻1AP,所以
|λE - B| &= |λE - P⁻1AP| = |P⁻1(λE - A)P| &= |P⁻1| |λE - A| |P| = |λE - A|.如果 A ∼ B,则 f(A) ∼ f(B),其中
f(x) = aₘ x^m + aₘ-1x^(m-1) + … + a₁ x + a₀.证
已知存在可逆矩阵 P,使得 B = P⁻1AP,所以
B^k = (P⁻1AP)^k = P⁻1A^kP, k = 1, 2, …, m.f(B) &= aₘB^m + aₘ-1B^(m-1) + … + a₁B + a₀E &= aₘP⁻1A^mP + aₘ-1P⁻1A^(m-1)P + … + a₁P⁻1AP + a₀P⁻1EP &= P⁻1(aₘA^m + aₘ-1A^(m-1) + … + a₁A + a₀E)P = P⁻1f(A)P.如果 A ∼ B,A 可逆,则 B 可逆,且 A⁻1 ∼ B⁻1。
证
已知存在可逆矩阵 P,使得 B = P⁻1AP,且 |A| ≠ 0。所以 |B| = |A| ≠ 0,B 可逆,且 B⁻1 = (P⁻1AP)⁻1 = P⁻1A⁻1P,即 A⁻1 ∼ B⁻1。
下面讨论当方阵 A 满足什么条件时,能找到可逆的矩阵 P,使得 P⁻1AP 为对角矩阵,此时称 A 可对角化。
二、方阵可对角化的条件
Section titled “二、方阵可对角化的条件”n 阶方阵 A 相似于对角矩阵的充分必要条件是 A 有 n 个线性无关的特征向量。
证明
必要性:如果 A 与对角矩阵 Λ = diag(λ₁, λ₂, …, λₙ) 相似,即存在可逆矩阵 P,使 P⁻1AP = Λ。设 P = (p₁, p₂, …, pₙ),由 P⁻1AP = Λ 可得 AP = PΛ,即
A(p₁, p₂, …, pₙ) = (p₁, p₂, …, pₙ)Λ,(Ap₁, Ap₂, …, Apₙ) = (λ₁p₁, λ₂p₂, …, λₙpₙ).由此可得 Apₖ = λₖpₖ (k = 1, 2, …, n)。这说明 P 的 n 个列向量是 A 的特征向量。而 P 是可逆的,所以 p₁, p₂, …, pₙ 线性无关,即 A 有 n 个线性无关的特征向量。
充分性:设 A 有 n 个线性无关的特征向量 p₁, p₂, …, pₙ,且
Apⱼ = λⱼpⱼ, j = 1, 2, …, n.令 P = (p₁, p₂, …, pₙ),则 P 是 n 阶可逆矩阵,而且满足
AP &= A(p₁, p₂, …, pₙ) = (λ₁p₁, λ₂p₂, …, λₙpₙ) &= (p₁, p₂, …, pₙ)Λ,其中 Λ = diag(λ₁, λ₂, …, λₙ) 为对角矩阵,即 P⁻1AP = Λ 为对角矩阵。
定理充分性的证明过程告诉我们:只要找出 A 的 n 个线性无关的特征向量 p₁, p₂, …, pₙ,以它们为列向量构成 P = (p₁, p₂, …, pₙ),则 P⁻1AP 为对角矩阵。这个对角矩阵称为 A 的相似标准形,其对角元素恰好为 A 的 n 个特征值。
在求矩阵 A 的相似标准形 Λ 时,Λ 的对角元可以是按任意顺序排列的。但是 P 的各列的排列次序与 Λ 中各个对角元(即 A 的全体特征值)的排列次序必须互相对应,P 的第 k 列所对应的特征值就是 Λ 的第 k 个对角元素 (k = 1, 2, …, n)。
如果存在可逆矩阵 P,使得 P⁻1AP = Λ 为对角矩阵,则称对角矩阵 Λ 为 A 的相似标准形。
定理1 告诉我们,n 阶方阵 A 相似于对角矩阵的充要条件是 A 有 n 个线性无关的特征向量。那么,对于给定的 n 阶方阵 A,如何知道它是否有 n 个线性无关的特征向量呢?先证明下面的定理。
设 λ₁, λ₂, …, λₖ 是 n 阶方阵 A 的两两不同的特征值,pᵢ 是 A 的属于 λᵢ 的特征向量 (1 ≤slant i ≤slant k),则 p₁, p₂, …, pₖ 线性无关。
证明
用数学归纳法。k=1 时定理显然成立。假设 k=m 时定理成立。当 k=m+1 时,设
a₁p₁ + a₂p₂ + … + aₘ+1pₘ+1 = 0, (6.1)则
A(a₁p₁ + a₂p₂ + … + aₘ+1pₘ+1) = 0, a₁Ap₁ + a₂Ap₂ + … + aₘ+1Apₘ+1 = 0, a₁λ₁p₁ + a₂λ₂p₂ + … + aₘ+1λₘ+1pₘ+1 = 0. (6.2)将式 (6.1) 等式两边乘 λₘ+1,再与式 (6.2) 相减得
a₁(λₘ+1 - λ₁)p₁ + a₂(λₘ+1 - λ₂)p₂ + … + aₘ(λₘ+1 - λₘ)pₘ = 0.由归纳假设可知 p₁, p₂, …, pₘ 线性无关,所以
a₁(λₘ+1 - λ₁) = 0, a₂(λₘ+1 - λ₂) = 0, …, aₘ(λₘ+1 - λₘ) = 0.因为 λₘ+1 - λᵢ ≠ 0 (i = 1, 2, …, m),所以 a₁ = a₂ = … = aₘ = 0,代入式 (6.1) 得 aₘ+1pₘ+1 = 0。
因为 pₘ+1 ≠ 0,故 aₘ+1 = 0。这表明 p₁, p₂, …, pₘ+1 线性无关,定理得证。
结合定理1 与定理2,有下面的推论:
如果 n 阶方阵 A 有 n 个两两不同的特征值,则 A 可对角化。
由定理2 还可得到下面的推论:
设 λ₁, λ₂, …, λₘ 是 n 阶方阵 A 的两两不同的特征值,pᵢ1, pᵢ2, …, pᵢrᵢ 是 A 的属于 λᵢ 的线性无关的特征向量 (1 ≤slant i ≤slant m),则
p₁1, p₁2, …, p₁r₁, p₂1, p₂2, …, p₂r₂, …, pₘ1, pₘ2, …, pₘrₘ(共 r₁ + r₂ + … + rₘ 个向量)仍线性无关。
证*
设
a₁1p₁1 + … + a₁r₁p₁r₁ + a₂1p₂1 + … + a₂r₂p₂r₂ + … + aₘ1pₘ1 + … + aₘrₘpₘrₘ = 0,令 pᵢ = aᵢ1pᵢ1 + aᵢ2pᵢ2 + … + aᵢrᵢpᵢrᵢ (i = 1, 2, …, m),则
p₁ + p₂ + … + pₘ = 0.下面用反证法证明 p₁ = p₂ = … = pₘ = 0。否则,若 p₁, p₂, …, pₘ 中有非零向量,设 pₖ₁, pₖ₂, …, pₖₜ (1 ≤slant t ≤slant m) 为所有非零向量,则 pₖ₁ + pₖ₂ + … + pₖₜ = 0。由 6.1 节性质1 可知 pₖ₁, pₖ₂, …, pₖₜ 是 A,根据定理2,λₖ₁, λₖ₂, …, λₖₜ,根据定理 2,pₖ₁, pₖ₂, …, pₖₜ 是线性无关的,这与 pₖ₁ + pₖ₂ + … + pₖₜ = 0 矛盾。因此 p₁ = p₂ = … = pₘ = 0,即
pᵢ = aᵢ1pᵢ1 + aᵢ2pᵢ2 + … + aᵢrᵢpᵢrᵢ = 0, i = 1, 2, …, m.已知对任意 1 ≤slant i ≤slant m,pᵢ1, pᵢ2, …, pᵢrᵢ 是线性无关的,所以
aᵢ1 = aᵢ2 = … = aᵢrᵢ = 0, i = 1, 2, …, m,推论3 得证。
推论3 告诉我们:将属于不同特征值的线性无关的特征向量合并在一起,可以组成个数更多的线性无关的特征向量组。那么,对应于每个特征值,最多能找出多少个线性无关的特征向量呢?下面的定理3 给出了这个问题的答案。
设 λ₀ 是 n 阶方阵 A 的重数为 k 的特征值,对应于 λ₀,A 最多有 m 个线性无关的特征向量,则 m ≤slant k。
证*
用反证法。设 m > k,p₁, p₂, …, pₘ 为 A 的对应于 λ₀ 的线性无关的特征向量。因为 p₁, p₂, …, pₘ 线性无关,继续添加向量 pₘ+1, pₘ+2, …, pₙ 使 p₁, p₂, …, pₙ 成为复向量空间 ℂⁿ 的一组基。
设 Apᵢ = b₁ip₁ + b₂ip₂ + … + bₙipₙ (i = m+1, m+2, …, n),则
A(p₁, p₂, …, pₙ) = (p₁, p₂, …, pₙ) λ₀Eₘ & B₁ O & B₂ , (6.3)其中
B₁ = b₁, m+1 & b₁, m+2 & … & b₁n b₂, m+1 & b₂, m+2 & … & b₂n vdots & vdots & & vdots bₘ, m+1 & bₘ, m+2 & … & bₘn , B₂ = bₘ+1, m+1 & bₘ+1, m+2 & … & bₘ+1, n bₘ+2, m+1 & bₘ+2, m+2 & … & bₘ+2, n vdots & vdots & & vdots bₙ, m+1 & bₙ, m+2 & … & bₙn .令 P = (p₁, p₂, …, pₙ),则 P 可逆,由式 (6.3) 可得
P⁻1AP = λ₀Eₘ & B₁ O & B₂ .记 B = λ₀Eₘ & B₁ O & B₂,则 A 与 B 相似,所以 A 与 B 特征值相同。
|λE - B| = (λ - λ₀)Eₘ & -B₁ O & λEₙ-m - B₂ ,将等式右端的行列式依次按第 1 列,第 2 列,…,第 m 列展开可得
|λE - B| = (λ - λ₀)^m |λEₙ-m - B₂|.因为 |λEₙ-m - B₂| 是关于 λ 的 (n-m) 次多项式,所以 λ₀ 是 B 的特征值,重数至少为 m,而 A 与 B 的特征值是相同的,因此 λ₀ 也是 A 的特征值,重数至少为 m,但 m > k,这与已知重数为 k 矛盾。
定理3 说明,对应于每个特征值,矩阵能有的线性无关的特征向量个数不会超过这个特征值的重数。由定理1 与定理3 可得下面的结论:
方阵 A 可对角化的充分必要条件是对对应于 A 的任意特征值 λ,A 有 k 个线性无关的特征向量,其中 k 是 λ 作为特征值的重数。
三、相似对角化方法与典型例题
Section titled “三、相似对角化方法与典型例题”设 A = 3 & -1 & -2 2 & 0 & -2 2 & -1 & -1,讨论 A 是否可以对角化。
解
λ₁ = λ₂ = 1 是二重特征值。
考虑齐次线性方程组 (E - A)x = 0,其系数矩阵为
E - A = -2 & 1 & 2 -2 & 1 & 2 -2 & 1 & 2 , 显然 r(E - A) = 1,因此齐次线性方程组 (E - A)x = 0 的基础解系中有两个向量,也就是说,对应于二重特征值 λ₁ = λ₂ = 1,A 有两个线性无关的特征向量,所以 A 可以对角化。
对于单特征值,必有属于它的特征向量,无须对它进行讨论。
设 A = -1 & 1 & 0 -4 & 3 & 0 1 & 0 & 2,讨论 A 是否可以对角化。
解
λ₁ = λ₂ = 1 是二重特征值。
考虑齐次线性方程组 (E - A)x = 0,其系数矩阵为
E - A = 2 & -1 & 0 4 & -2 & 0 -1 & 0 & -1 , r(E - A) = 2,因此齐次线性方程组 (E - A)x = 0 的基础解系中只有一个向量,也就是说,对应于二重特征值 λ₁ = λ₂ = 1,A 只有一个线性无关的特征向量,所以 A 不能对角化。
设 A = 1 & 0 & 0 -1 & 2 & 0 1 & 1 & -1,求 A^k,k 为正整数。
解
A 为下三角矩阵,主对角元素就是其特征值,A 的特征值两两不同,所以 A 可以对角化。下面将 A 对角化。
先求出 A 的特征值和特征向量。A 的特征方程为
|λE - A| = λ - 1 & 0 & 0 1 & λ - 2 & 0 -1 & -1 & λ + 1 = (λ - 1)(λ - 2)(λ + 1) = 0.属于特征值 λ₁ = 1 的特征向量是方程组 (E - A)x = 0,即
0 & 0 & 0 1 & -1 & 0 -1 & -1 & 2 x₁ x₂ x₃ = 0 0 0的非零解向量。取特征向量 p₁ = (1, 1, 1)^T。
属于特征值 λ₂ = 2 的特征向量是方程组 (2E - A)x = 0,即
1 & 0 & 0 1 & 0 & 0 -1 & -1 & 3 x₁ x₂ x₃ = 0 0 0的非零解向量。取特征向量 p₂ = (0, 3, 1)^T。
属于特征值 λ₃ = -1 的特征向量是方程组 (-E - A)x = 0,即
-2 & 0 & 0 1 & -3 & 0 -1 & -1 & 0 x₁ x₂ x₃ = 0 0 0的非零解向量。取特征向量 p₃ = (0, 0, 1)^T。
令
P = (p₁, p₂, p₃) = 1 & 0 & 0 1 & 3 & 0 1 & 1 & 1 , 则 P⁻1AP = 1 & 0 & 0 0 & 2 & 0 0 & 0 & -1 = Λ, 即 A = PΛP⁻1.A^k &= (PΛP⁻1)^k = PΛ^kP⁻1 = 1 & 0 & 0 1 & 3 & 0 1 & 1 & 1 1 & 0 & 0 0 & 2^k & 0 0 & 0 & (-1)^k 1 & 0 & 0 -1/3 & 1/3 & 0 -2/3 & -1/3 & 1 &= 1 & 0 & 0 1 - 2^k & 2^k & 0 1 - (2^k)/(3) + frac2(-1)^(k+1)3 & frac2^k + (-1)^(k+1)3 & (-1)^k .设 A = 0 & -2 & -1 -2 & 3 & 2 -1 & 2 & 0,求可逆矩阵 P,使 P⁻1AP 为对角形。
解
A 的特征值为 λ₁ = 5, λ₂ = λ₃ = -1。
对应于特征值 λ₁ = 5 的特征向量是方程组 (5E - A)x = 0,即
5 & 2 & 1 2 & 2 & -2 1 & -2 & 5 x₁ x₂ x₃ = 0 0 0的非零解向量。方程组的基础解系为 ξ₁ = (-1, 2, 1)^T,取 p₁ = ξ₁。
对应于特征值 λ₂ = λ₃ = -1 的特征向量是方程组 (-E - A)x = 0,即
-1 & 2 & 1 2 & -4 & -2 1 & -2 & -1 x₁ x₂ x₃ = 0 0 0的非零解向量。方程组的基础解系为 ξ₂ = (1, 0, 1)^T, ξ₃ = (1, 1, -1)^T,取 p₂ = ξ₂, p₃ = ξ₃。
令
P = (p₁, p₂, p₃) = -1 & 1 & 1 2 & 0 & 1 1 & 1 & -1 , 则 P⁻1AP = 5 & 0 & 0 0 & -1 & 0 0 & 0 & -1 .这个例子中的矩阵 A 是实对称矩阵,这个矩阵 A 可以对角化并非偶然,实对称矩阵一定能相似于对角矩阵的,这个问题将在下一章中进行讨论。
已知三阶矩阵 A 的特征值与对应的特征向量分别为 λ₁ = 1, p₁ = (1, 2, 2)^T ; λ₂ = -1, p₂ = (0, 1, 1)^T ; λ₃ = 2, p₃ = (2, 1, -1)^T。求 A^T 的特征值与特征向量。
解
令
P = (p₁, p₂, p₃) = 1 & 0 & 2 2 & 1 & 1 2 & 1 & -1 , 则 P⁻1AP = 1 & 0 & 0 0 & -1 & 0 0 & 0 & 2 ,即
A = P 1 & 0 & 0 0 & -1 & 0 0 & 0 & 2 P⁻1.两边取转置,得
A^T = (P⁻1)^T 1 & 0 & 0 0 & -1 & 0 0 & 0 & 2 P^T = (P^T)⁻1 1 & 0 & 0 0 & -1 & 0 0 & 0 & 2 P^T,即
P^TA^T(P^T)⁻1 = 1 & 0 & 0 0 & -1 & 0 0 & 0 & 2 .记 Q = (P^T)⁻1,则
Q⁻1A^TQ = 1 & 0 & 0 0 & -1 & 0 0 & 0 & 2 .这说明 Q = (P^T)⁻1 的三个列向量就是 A^T 的特征向量。经计算得
Q = (P^T)⁻1 = (P⁻1)^T = 1/2 2 & -4 & 0 -2 & 5 & 1 2 & -3 & -1 .因此 A^T 的特征值与对应的特征向量分别为:
- λ₁ = 1, α₁ = k₁(1, -1, 1)^T (k₁ ≠ 0);
- λ₂ = -1, α₂ = k₂(-4, 5, -3)^T (k₂ ≠ 0);
- λ₃ = 2, α₃ = k₃(0, 1, -1)^T (k₃ ≠ 0)。
当矩阵 A 可对角化时,可逆矩阵 P 能使 P⁻1AP 为对角形的充分必要条件是 P 的列向量是 A 的特征向量。所以该例中并不要求求出 A^T,只要根据 Q⁻1A^TQ = P^TA^T(P^T)⁻1 为对角形,便可知 Q = (P^T)⁻1 的列向量就是 A^T 的特征向量。