MathLabs

代数学

特征值与对角化

揭示线性变换最简单形式的特殊标量与向量。

直观特征向量到底是什么

线性映射 AA 作用于大多数向量时,会把它们推向一个全新的方向,通常还伴随旋转与伸缩。但对少数特殊方向而言,该映射只沿着向量原来所在的那条直线拉伸或压缩,绝不旋转它。这些特殊方向就是 AA 的特征向量,沿每个方向的伸缩系数就是对应的特征值: Av=λvA\mathbf{v} = \lambda \mathbf{v}。可以想象一张被拉伸的橡皮膜:大多数标记点会被推向侧边,但膜上有几条线只是变长或变短,始终停留在同一条轴线上。这些轴线正是特征向量。

由二次型生成的三维碗状曲面,可旋转以显示其主轴。
对于对称矩阵(a12=a21a_{12}=a_{21}),两条绿色虚线特征向量轴相互正交:位于轴上的向量仅按 λ1,λ2\lambda_1, \lambda_2 伸缩而不发生旋转。

中学从几何图像到精确定义

定义: 特征值与特征向量

设 AA 是一个 n×nn \times n 矩阵。若存在标量 λ\lambda 使得 Av=λvA\mathbf{v} = \lambda \mathbf{v} 成立,则称非零向量 v\mathbf{v} 为 AA 的特征向量;标量 λ\lambda 是与 v\mathbf{v} 相伴的特征值。从几何上看,AA 把经过 v\mathbf{v} 的整条直线映到自身,只是按系数 λ\lambda 缩放。

Av=λvA\mathbf{v} = \lambda \mathbf{v}

把定义式改写为 (A−λI)v=0(A - \lambda I)\mathbf{v} = \mathbf{0} 可以看出,只有当 A−λIA - \lambda I 奇异、即其行列式为零时,才存在 v≠0\mathbf{v} \neq \mathbf{0}: det⁡(A−λI)=0\det(A - \lambda I) = 0。展开这个行列式得到一个关于 λ\lambda 的 nn 次多项式,称为特征多项式;它的根正是 AA 的特征值。

det⁡(A−λI)=0\det(A - \lambda I) = 0
特征值的类型及其几何意义
情形几何意义
nn 个互异实特征值nn 条独立的伸缩轴;AA 必可对角化
重特征值但特征空间维数足够仍可对角化:整个平面(或更高维)按相同系数伸缩
重特征值但特征空间维数不足(亏损矩阵)不可对角化:某些方向不仅伸缩,还发生剪切
共轭复特征值对不存在实的不变直线;映射一边旋转一边缩放,呈螺旋状

大学对角化与谱定理

若 λ1,…,λk\lambda_1, \dots, \lambda_k 是 AA 两两互异的特征值,对应特征向量为 v1,…,vk\mathbf{v}_1, \dots, \mathbf{v}_k,则集合 {v1,…,vk}\{\mathbf{v}_1, \dots, \mathbf{v}_k\} 线性无关。

为什么成立?

每个特征向量都标出了自己专属的不变方向,并配有专属的伸缩系数;如果某个方向能由其他方向组合而成,那么作用 AA 时就必须同时按所有这些不同系数缩放它,除非该向量为零,否则这是不可能的。

证明

我们对 kk 用归纳法证明。当 k=1k = 1 时结论显然成立:单个非零向量必然线性无关。

假设结论对 k−1k - 1 个互异特征值成立,反设 v1,…,vk\mathbf{v}_1, \dots, \mathbf{v}_k 线性相关。于是存在关系式 c1v1+⋯+ckvk=0c_1 \mathbf{v}_1 + \cdots + c_k \mathbf{v}_k = \mathbf{0},且并非所有 cic_i 都为零。

对两边作用 AA:由 Avi=λiviA\mathbf{v}_i = \lambda_i \mathbf{v}_i,得到 c1λ1v1+⋯+ckλkvk=0c_1 \lambda_1 \mathbf{v}_1 + \cdots + c_k \lambda_k \mathbf{v}_k = \mathbf{0}。将原关系式乘以 λk\lambda_k 后与这个新方程相减,含 vk\mathbf{v}_k 的项恰好抵消,剩下 c1(λ1−λk)v1+⋯+ck−1(λk−1−λk)vk−1=0c_1(\lambda_1 - \lambda_k)\mathbf{v}_1 + \cdots + c_{k-1}(\lambda_{k-1} - \lambda_k)\mathbf{v}_{k-1} = \mathbf{0}。

由归纳假设,v1,…,vk−1\mathbf{v}_1, \dots, \mathbf{v}_{k-1} 线性无关,故最后这个关系式中每个系数都必须为零:对所有 i<ki < k 有 ci(λi−λk)=0c_i(\lambda_i - \lambda_k) = 0。由于特征值两两互异,λi−λk≠0\lambda_i - \lambda_k \neq 0,因此对所有 i<ki < k 都有 ci=0c_i = 0。

代回原关系式,只剩下 ckvk=0c_k \mathbf{v}_k = \mathbf{0};由特征向量定义 vk≠0\mathbf{v}_k \neq \mathbf{0},故 ck=0c_k = 0 也成立。所有系数皆为零,与"并非所有 cic_i 为零"的假设矛盾。因此不存在这样的相关关系,特征向量线性无关。

若 AA 是 n×nn \times n 实对称矩阵,即 A⊤=AA^\top = A,则存在正交矩阵 QQ(即 Q⊤Q=IQ^\top Q = I)与实对角矩阵 Λ\Lambda,使得 A=QΛQ⊤A = Q \Lambda Q^\top。等价地,AA 有 nn 个实特征值,并存在一组由特征向量构成的标准正交基。

为什么成立?

对称矩阵无处不在——协方差矩阵、转动惯量、光滑函数的海森矩阵——这条定理保证我们总能旋转到一个坐标系,使得该矩阵仅沿相互垂直的坐标轴做纯粹的独立缩放,没有剪切,也没有任何复杂行为。

证明

我们先给出一个引理:实对称矩阵的特征值都是实数。设 Av=λvA\mathbf{v} = \lambda \mathbf{v},v≠0\mathbf{v} \neq \mathbf{0} 可能是复向量,考虑 v∗Av\mathbf{v}^{*}A\mathbf{v},其中 v∗\mathbf{v}^{*} 是共轭转置。由于 AA 是实对称矩阵,v∗Av=(v∗Av)∗\mathbf{v}^{*}A\mathbf{v} = (\mathbf{v}^{*}A\mathbf{v})^{*},故该量为实数;但它也等于 λ v∗v\lambda \, \mathbf{v}^{*}\mathbf{v},而 v∗v>0\mathbf{v}^{*}\mathbf{v} > 0 是正实数,于是 λ\lambda 本身必为实数。

下面对 nn 用归纳法证明定理。n=1n = 1 时结论平凡:任何 1×11 \times 1 矩阵已经是对角矩阵,取 Q=(1)Q = (1) 即可。

归纳步骤:假设定理对所有 n−1n - 1 阶实对称矩阵成立。由于 AA 是 n×nn \times n 实对称矩阵,由引理其特征多项式有实根 λ1\lambda_1;取相应的特征向量并归一化为单位向量 v1\mathbf{v}_1。

设 WW 是 v1\mathbf{v}_1 张成的直线的正交补,是一个 (n−1)(n-1) 维子空间。我们断言 WW 在 AA 下不变:对任意 w∈W\mathbf{w} \in W(即 w⊤v1=0\mathbf{w}^\top \mathbf{v}_1 = 0),计算 (Aw)⊤v1=w⊤A⊤v1=w⊤Av1=λ1w⊤v1=0(A\mathbf{w})^\top \mathbf{v}_1 = \mathbf{w}^\top A^\top \mathbf{v}_1 = \mathbf{w}^\top A \mathbf{v}_1 = \lambda_1 \mathbf{w}^\top \mathbf{v}_1 = 0,中间一步用到了 AA 的对称性。因此 AwA\mathbf{w} 仍与 v1\mathbf{v}_1 正交,即 Aw∈WA\mathbf{w} \in W。

取 WW 的一组标准正交基;在此基下,AA 限制在 WW 上由一个 (n−1)×(n−1)(n-1) \times (n-1) 矩阵 A′A' 表示,且 A′A' 也是对称的(AA 对称,把对称双线性型限制到子空间上,在标准正交基下仍保持对称)。由归纳假设,A′A' 在 WW 内有一组由特征向量 v2,…,vn\mathbf{v}_2, \dots, \mathbf{v}_n 构成的标准正交基,对应实特征值 λ2,…,λn\lambda_2, \dots, \lambda_n;由于 WW 是 AA-不变的,这些向量也是 AA 本身真正的特征向量。

将 v1,v2,…,vn\mathbf{v}_1, \mathbf{v}_2, \dots, \mathbf{v}_n 汇集起来,就得到了 Rn\mathbb{R}^n 的一组完全由 AA 的特征向量构成的标准正交基。把它们作为矩阵 QQ 的列,则 QQ 是正交矩阵,且 AQ=QΛAQ = Q\Lambda(其中 Λ=diag(λ1,…,λn)\Lambda = \mathrm{diag}(\lambda_1, \dots, \lambda_n));由于正交矩阵满足 Q−1=Q⊤Q^{-1} = Q^\top,这恰好整理为 A=QΛQ⊤A = Q \Lambda Q^\top,归纳完成。

大学实际应用与典型例题

特征值与特征向量是桥梁和建筑振动频率、网页搜索排名算法以及数据压缩所用主方向背后的隐藏变量。

例题: 双质量弹簧系统的固有频率

两个质量均为 mm 的质点在无摩擦直线上滑动;每个质点通过刚度为 kk 的弹簧连接到墙壁,两个质点之间又用刚度为 kk 的第三根弹簧相连。求振动的固有频率,并描述对应的振动模态。

解答

对位移 x1,x2x_1, x_2 应用牛顿第二定律,得 mx¨1=−2kx1+kx2m\ddot{x}_1 = -2kx_1 + kx_2 与 mx¨2=kx1−2kx2m\ddot{x}_2 = kx_1 - 2kx_2,即 mx¨=−Kxm\ddot{\mathbf{x}} = -K\mathbf{x},其中 K=(2k−k−k2k)K = \begin{pmatrix} 2k & -k \\ -k & 2k \end{pmatrix}。

寻求振动解 x(t)=vcos⁡(ωt)\mathbf{x}(t) = \mathbf{v}\cos(\omega t) 把系统化为特征值问题 Kv=mω2vK\mathbf{v} = m\omega^2 \mathbf{v},因此 mω2m\omega^2 必须是 KK 的特征值。特征方程 det⁡(K−μI)=(2k−μ)2−k2=0\det(K - \mu I) = (2k-\mu)^2 - k^2 = 0 给出 μ=3k\mu = 3k 或 μ=k\mu = k。

当 μ=3k\mu = 3k 时:(K−3kI)v=0(K - 3kI)\mathbf{v} = 0 给出 −kv1−kv2=0-k v_1 - k v_2 = 0,即 v=(1,−1)\mathbf{v} = (1,-1)——两质点反相摆动,ω1=3k/m\omega_1 = \sqrt{3k/m}。当 μ=k\mu = k 时:(K−kI)v=0(K - kI)\mathbf{v}=0 给出 kv1−kv2=0k v_1 - k v_2 = 0,即 v=(1,1)\mathbf{v} = (1,1)——两质点同相摆动,ω2=k/m\omega_2 = \sqrt{k/m}。

因此 KK 的两个特征向量正是两个简正振动模态,其特征值确定了两个固有频率——这正是该问题本质上是特征值问题的原因。

例题: 作为主特征向量的PageRank

一个由三个页面 A,B,CA, B, C 组成的小型网络,链接如下:AA 平均链接到 BB 和 CC;BB 只链接到 CC;CC 只链接到 AA。用列随机转移矩阵为随机浏览者建模,求每个页面长期占用时间的比例。

解答

列随机转移矩阵(每一列给出该页面的出链概率)为 P=(0010.5000.510)P = \begin{pmatrix} 0 & 0 & 1 \\ 0.5 & 0 & 0 \\ 0.5 & 1 & 0 \end{pmatrix}(行/列顺序为 A,B,CA,B,C)。稳态排名 π\boldsymbol{\pi} 是 PP 对应特征值 11 的特征向量,满足 Pπ=πP\boldsymbol{\pi} = \boldsymbol{\pi} 且 πA+πB+πC=1\pi_A + \pi_B + \pi_C = 1。

展开 Pπ=πP\boldsymbol{\pi} = \boldsymbol{\pi}:πC=πA\pi_C = \pi_A,0.5πA=πB0.5\pi_A = \pi_B,以及 0.5πA+πB=πC0.5\pi_A + \pi_B = \pi_C。前两个方程已给出 πC=πA\pi_C = \pi_A 和 πB=0.5πA\pi_B = 0.5\pi_A;代入第三个方程验证 0.5πA+0.5πA=πA=πC0.5\pi_A + 0.5\pi_A = \pi_A = \pi_C,一致。

令 πA=x\pi_A = x,则 πB=0.5x\pi_B = 0.5x,πC=x\pi_C = x,归一化 x+0.5x+x=2.5x=1x + 0.5x + x = 2.5x = 1 给出 x=0.4x = 0.4。稳态排名为 (πA,πB,πC)=(0.4,0.2,0.4)(\pi_A, \pi_B, \pi_C) = (0.4, 0.2, 0.4):页面 AA 与 CC 同样"最重要",这正符合预期,因为它们在一个紧密的循环中直接相互链接,而 BB 只从 AA 获得流量。

这正是 PageRank 背后的思想:重要性不是靠数链接数来衡量,而是通过求链接转移矩阵的主特征向量(特征值 11)来衡量。实践中通过反复乘以 PP(幂法)求得,无论初始猜测如何,总会收敛到同一个主特征向量。

A=(4123)A = \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix} 的特征值是什么?

在线性映射 AA 下,区分特征向量与普通向量的关键几何性质是什么?

在 PageRank 模型中,为什么网页的稳态重要性排名对应于一个特征向量?

为什么矩阵 (1101)\begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix} 尽管有重特征值 λ=1\lambda = 1,却不可对角化?

参考文献

  1. Gilbert Strang (2016). Introduction to Linear Algebra
  2. Sheldon Axler (2015). Linear Algebra Done Right