n × n 矩阵的特征方程含有 n 次多项式,如果考虑复根,方程恰好有 n 个根(重根重复计算)。本节将看到假如实矩阵 A 的特征方程有复根,那么这些复根将给我们提供有关 A 的关键信息。关键的问题是让 A 作用于 n 维复空间 ℂⁿ 中的 n 元复数。
我们对 ℂⁿ 感兴趣并不是为了把前几章的结果推广,尽管这可以开辟线性代数应用的新领域。然而,对复特征值的研究能使我们揭示各种实际生活问题中出现的某些实矩阵中“隐藏”的信息。这些问题包括很多蕴涵周期运动的实动力系统、振动或空间的某种旋转。
建立在 ℝⁿ 基础上的矩阵特征值-特征向量理论同样可以很好地应用到 ℂⁿ 。因此,一个复数 λ 满足 det(A - λ I) = 0 当且仅当在 ℂⁿ 中存在一个非零向量 x ,使得 Ax = λ x 。我们称这样的 λ 是(复)特征值, x 是对应于 λ 的(复)特征向量。
假设 A = 0 & -1 1 & 0 ,那么 ℝ² 上的线性变换 x mapsto Ax 将平面逆时针旋转 1/4 圈。 A 的作用是周期性的,因为在旋转 4 次 1/4 圈后,向量又回到了它的初始位置,很明显没有非零向量被映射成自身的数倍,所以 A 在 ℝ² 中没有特征向量,因此也没有实的特征值。实际上, A 的特征方程是
λ^ 2 + 1 = 0只有复根 λ = i 和 λ = -i 但是,如果我们让 A 作用在 ℂ² 上,那么
[ c c 0 & - 1 1 & 0 ] [ c 1 - i ] = [ c i 1 ] = i [ c 1 - i ][ c c 0 & - 1 1 & 0 ] [ l 1 i ] = [ c - i 1 ] = - i [ l 1 i ]因此i和-i是特征值, 1 -i 和 1 i 是对应的特征向量(在例2中讨论求复特征向量的方法).
下一个例子中的矩阵是本节的主要焦点所在.
设 A = 0.5 & -0.6 0.75 & 1.1 ,求 A 的特征值及每个特征空间的基.
解
A 的特征方程是
0 = det [ c c 0. 5 - λ & - 0. 6 0. 7 5 & 1. 1 - λ ] = (0. 5 - λ) (1. 1 - λ) - (- 0. 6) (0. 7 5) = λ^ 2 - 1. 6 λ + 1依求根公式,得 λ=1/2[1.6±sqrt(-1.6)²-4]=0.8±0.6i 。对特征值 λ=0.8-0.6i ,构造
r l A - (0. 8 - 0. 6 i) I & = [ c c 0. 5 & - 0. 6 0. 7 5 & 1. 1 ] - [ c c 0. 8 - 0. 6 i & 0 0 & 0. 8 - 0. 6 i ] & = [ c c - 0. 3 + 0. 6 i & - 0. 6 0. 7 5 & 0. 3 + 0. 6 i ] (1)由于有复数运算,手工对增广矩阵做行化简是相当慢的。不过,细心观察就能找到简化问题的方法:因为0.8-0.6i是特征值,所以方程组
r l (- 0. 3 + 0. 6 i) x _ 1 - & 0. 6 x _ 2 = 0 0. 7 5 x _ 1 + (0. 3 + 0. 6 i) x _ 2 = 0 (2)有非零解(这里的 x₁ 和 x₂ 可能是复数),因此,(2)中的2个方程确定的 x₁ 与 x₂ 之间的关系是同一关系,这样就可以通过其中的一个方程将某个变量用另一个来表示.
由(2)的第2个方程得
r l 0. 7 5 x _ 1 & = (- 0. 3 - 0. 6 i) x _ 2 x _ 1 & = (- 0. 4 - 0. 8 i) x _ 2为去掉小数,取 x₂ = 5 ,有 x₁ = -2 - 4i 。对应于 λ = 0.8 - 0.6i 的特征空间的基是
ν_ 1 = [ c - 2 - 4 i 5 ]对 λ = 0.8 + 0.6i 做同样的计算,可求出特征向量
boldsymbol v _ 2 = [ c - 2 + 4 i 5 ]为验证结果是否正确,可以计算
A ±b v _ 2 = [ c c 0. 5 & - 0. 6 0. 7 5 & 1. 1 ] [ c - 2 + 4 i 5 ] = [ c - 4 + 2 i 4 + 3 i ] = (0. 8 + 0. 6 i) ±b v _ 2例2 的矩阵 A 确定的变换 x mapsto Ax 实质上是一个旋转变换,画出一些适当的点后,这一事实会变得更明显.
例3 与例2中的矩阵 A 相乘会影响点,通过以下方法可以看到这种影响。先画出一个随机的初始点(如 x₀ = (2,0) ),然后画出重复乘 A 后该点的后续图像。即,画出
l boldsymbol x _ 1 = A boldsymbol x _ 0 = [ c c 0. 5 & - 0. 6 0. 7 5 & 1. 1 ] [ c 2 0 ] = [ c 1. 0 1. 5 ] boldsymbol x _ 2 = A boldsymbol x _ 1 = [ c c 0. 5 & - 0. 6 0. 7 5 & 1. 1 ] [ c 1. 0 1. 5 ] = [ c - 0. 4 2. 4 ] boldsymbol x _ 3 = A boldsymbol x _ 2 …在图 5-11 中,用粗点显示 x₀, …, x₈ ,用细点显示 x₉, …, x₁00 ,序列分布在椭圆形的轨道上.
当然,图 5-11 没有解释为什么会发生旋转,旋转的秘密隐藏在复特征向量的实部和虚部。向量的实部和虚部
ℂⁿ 中复向量 x 的共轭向量 overlinex 也是 ℂⁿ 中的向量,它的分量是 x 中对应分量的共轭复数,向量 operatornameRe x 和 operatornameIm x 称为复向量 x 的实部和虚部,分别由 x 的分量的实部和虚部组成.
假如 x = 3 - i i 2 + 5i = 3 0 2 +i -1 1 5 ,那么
operatornameRe boldsymbol x = [ l 3 0 2 ], operatornameIm boldsymbol x = [ c - 1 1 5 ], overline boldsymbol x = [ l 3 0 2 ] - i [ c - 1 1 5 ] = [ c 3 + i - i 2 - 5 i ]假设 B 是可能有复元素的 m × n 矩阵,那么,以 B 的元素的共轭复数为元素的矩阵记为 overlineB.复数的共轭运算性质对复矩阵代数亦成立:
overline r x = overline r overline x, overline B x = overline B overline x, overline B C = overline B overline C, overline r B = overline r overline B作用于 ℂⁿ 上的实矩阵的特征值和特征向量
设 A 为 n × n 的实矩阵,那么 overlineAx = overlineAoverlinex = Aoverlinex ,假如 λ 是 A 的特征值,x 是对应于 λ 的特征向量,那么
overline A x = overline A x = overline λ x = overline λ overline x故 barλ 同样是 A 的特征值,而 barx 是对应的特征向量。这表明,当 A 是实矩阵时,它的复特征值以共轭复数对出现(在这里或别的地方,我们用术语复特征值来表示形如 λ = a + bi (b ≠ 0) 的特征值)。
例2的实矩阵的特征值是共轭复数对 0.8 - 0.6i 和 0.8 + 0.6i ,对应的特征向量同样是共轭复向量:
±b v _ 1 = [ c - 2 - 4 i 5 ] text 和 ±b v _ 2 = [ c - 2 + 4 i 5 ] = overline ±b v _ 1为计算有复特征值的 2 × 2 实矩阵提供了基本模式.
设 C = a & -b b & a ,其中 a, b 为实数且不都等于零。那么 C 的特征值是 λ = a ± bi 。(见本节末尾的练习题。)同样,假如 r = |λ| = √(a² + b²) ,那么
C = r [ l l a / r & - b / r b / r & a / r ] = [ l l r & 0 0 & r ] [ l l cos φ & - sin φ sin φ & cos φ ]这里的 φ 是正 x 轴与 (0,0) 到 (a,b) 射线的夹角.参看图 5-12和附录B,角 φ 称为 λ = a± bi 的辐角.因此变换 xmapsto Cx 可看作由旋转 φ 变换复合而成(见图 5-13).|λ | 变换复合而成(见图5-13).
最后,我们准备讨论有复特征值的实矩阵中隐含的旋转.
例7 与例2一样, A = 0.5 & -0.6 0.75 & 1.1 ,λ = 0.8 - 0.6i,ν₁ = -2 - 4i 5 同时,设 P 是 2× 2 实矩阵,
P = [ operatornameRe boldsymbol v _ 1 operatornameIm boldsymbol v _ 1 ] = [ c c - 2 & - 4 5 & 0 ]并令
C = P ^ - 1 A P = frac 12 0 [ l l 0 & 4 - 5 & - 2 ] [ l l 0. 5 & - 0. 6 0. 7 5 & 1. 1 ] [ l l - 2 & - 4 5 & 0 ] = [ l l 0. 8 & - 0. 6 0. 6 & 0. 8 ]由例6,因为 |λ |² = (0.8)² +(0.6)² = 1 ,故 C 仅是旋转变换.由 C = P⁻1AP ,得
A = P C P ^ - 1 = P [ c c 0. 8 & - 0. 6 0. 6 & 0. 8 ] P ^ - 1A “含有”旋转!矩阵 P 提供变量代换,如 x = Pu 。 A 的作用相当于将 x 代换为 u ,再经过旋转,然后又代换回初始变量,见图 5-14。旋转产生一个椭圆,如图 5-11 所示,而不是圆,因为由 P 的列确定的坐标系不是长方形的,在两个轴上没有相等的单位长。
三、备忘录
下列定理表明例7 的计算适用于任意有复特征值 λ 的 2 × 2 实矩阵.证明利用了如下结论:
若 A 是实矩阵,则 A(operatornameRe x) = operatornameRe A x 和 A(operatornameIm x) = operatornameIm A x ,若 x 是对应于复特征值的特征向量,则 operatornameRe x 和 operatornameIm x 是线性无关的。(见习题25和26.)这里省略细节。
设 A 是 2 × 2 实矩阵,有复特征值 λ = a - bi (b ≠ 0) 及对应的 C² 中的复特征向量 v,那么
A = PCP⁻1 ,其中 P = [Rev Imv] , C = a & -b b & a
在更高维矩阵中亦存在例7中所示的现象.例如,若 A 是有复特征值的 3 × 3 矩阵,那么在 ℝ³ 中存在某个平面,A 对平面的作用是旋转 (可能还结合倍乘), 平面中的每个向量被旋转到该平面的另一点上,我们说平面在 A 的作用下是不变的.
矩阵 A = 0.8 & -0.6 & 0 0.6 & 0.8 & 0 0 & 0 & 1.07 有特征值 0.8 ± 0.6i 和1.07. x₁x₂ 平面(第3坐标为0)的任一向量 ±bw₀ 被 A 旋转到该平面的另一位置上.不在该平面的任一向量 ±bx₀ 的 x₃ 坐标乘1.07,图 5-15显示了点 ±bw₀ = (2,0,0) 和 ±bx₀ = (2,0,1) 乘 A 后的迭代结果.
证明:若 a 和 b 是实数,则 A = a & -b b & a 的特征值是 a ± bi ,对应的特征向量是 1 -i 和 1 i .
让习题 1~6 的每个矩阵作用于 C² ,求矩阵的特征值及对应于 C² 中特征空间的基.
[ c c 1 & - 2 1 & 3 ][ c c 1 & 5 - 2 & 3 ][ c c 5 & - 5 1 & 1 ][ c c 5 & - 2 1 & 3 ]5. [ c c 0 & 1 - 8 & 4 ][ c c 4 & 3 - 3 & 4 ]在习题 7~12 中,用例6 写出 A 的特征值。在每一例中,变换 x mapsto Ax 由旋转加倍乘复合而成。求旋转的角度 φ (-π < φ ≤slant π) 和倍乘因子 r。
- [ cc √(3) & -1 1 & √(3) ] 8. [ cc √(3) & 3 -3 & √(3) ] 9. [ cc -√(3) / 2 & 1 / 2 -1 / 2 & -√(3) / 2 ] 10. [ cc -5 & -5 5 & -5 ] 11. [ cc 0.1 & 0.1 -0.1 & 0.1 ] 12. [ cc 0 & 0.3 -0.3 & 0 ]
在习题 13~20 中,求可逆矩阵 P 和形如 a & -b b & a 的矩阵 C,把所给的矩阵表示为 A = PCP⁻1 的形式。对于习题 13~16,使用习题 1~4 中的信息。 13. 1 & -2 1 & 3 14. 5 & -5 1 & 1 15. 1 & 5 -2 & 3 16. 5 & -2 1 & 3 17. -1 & -0.8 4 & -2.2 18. 1 & -1 0.4 & 0.6 19. 1.52 & -0.7 0.56 & 0.4 20. -1.64 & -2.4 1.92 & 2.2
-
在例2 中,解(2)中的第 1 个方程,用 x₁ 表示 x₂ ,并由此求得 A 的特征向量 y=2 -1+2i 是例2 中 y 是例 2 中 v₁ 的(复)数倍.
-
设 A 是 n × n 的复(或实)矩阵,x 是复特征值对应于 Cⁿ 中的特征值 λ 的复特征向量。证明:对任意非零的复数 μ ,向量 μ x 是 A 的特征向量。
第 7 章将把重点放在具有性质 A^(T) = A 的矩阵 A 上,习题 23 和 24 证明了这种矩阵的特征值必定是实数。
- 设 n × n 实矩阵 A 有性质 A^(T) = A, x 是 ℂⁿ 中的向
量,令 q = overlinex^(T) A x 。下面的等式通过验证 overlineq = q 证明了 q 是实数。给出每一步的理由。 overlineq=overlineoverlinex^(T)Aoverlinex=overlinex^(T)Aoverlinex=overlinex^(T)Aoverlinex=(x^(T)Aoverlinex)^(T)=overlinex^(T)A^(T)x=q
-
设 n × n 实矩阵 A 有性质 A^(T) = A 。证明,若对 ±bx ∈ ℂⁿ 有 Ax = λ x ,则 λ 是实数,而 ±bx 的实部是 A 的特征向量。(提示:计算 overline±bx^(T) Ax 和利用习题23的结果,并检查 Ax 的实部和虚部。)
-
设 A 是 n × n 实矩阵, x ∈ Cⁿ ,证明 operatornameRe(Ax) = A(operatornameRex) 和 operatornameIm(Ax) = A(operatornameImx) .
-
设 A 是 2 × 2 实矩阵,有复特征值 λ = a - bi (b ≠ 0) 和对应的 C² 中的复特征向量 v, a.证明 A(Rev) = aRev + bImv 和 A(Imv) = -bRev + aImv . (提示: 记 v = Rev + iImv ,计算 Av.) b. 假设 P 和 C 是定理9 给出的矩阵,证明 AP = PC.
[M]在习题27和习题28中,求所给矩阵 A 的分解式 A = PCP⁻1 ,其中 C 是由形如例6矩阵的 2 × 2 子矩阵组成的分块对角矩阵。(对每一对共轭复特征值,利用属于 ℂ⁴ 的特征向量的实部和虚部来产生 P 的两列。)
A = [ r r r r 0. 7 & 1. 1 & 2. 0 & 1. 7 - 2. 0 & - 4. 0 & - 8. 6 & - 7. 4 0 & - 0. 5 & - 1. 0 & - 1. 0 1. 0 & 2. 8 & 6. 0 & 5. 3 ]A = [ c c c c - 1. 4 & - 2. 0 & - 2. 0 & - 2. 0 - 1. 3 & - 0. 8 & - 0. 1 & - 0. 6 0. 3 & - 1. 9 & - 1. 6 & - 1. 4 2. 0 & 3. 3 & 2. 3 & 2. 6 ]记住,要验证向量是否为特征向量是很容易的,并不需要检验特征方程.计算
A boldsymbol x = [ c c a & - b b & a ] [ l 1 - i ] = [ l a + b i b - a i ] = (a + b i) [ l 1 - i ]因此, 1 -i 是对应于 λ=a+bi 的特征向量,从本节讨论的内容知道, 1 i 一定是对应于 overlineλ=a-bi 的特征向量.