LeoMath

线性代数基础

特征值

变换只是拉伸的方向。

约 7 分钟

从一个问题开始

问题

一个线性变换在大多数方向上既拉伸又旋转。有没有某些方向,它只做拉伸?

如果有,也许可以用这些方向重建坐标系,让变换在新坐标下变得极其简单。我们需要一个定义来抓住"方向不变"这件事,以及一种找到这些方向的方法。

观察

勾选"显示特征向量",然后按"拉伸":出现两条虚线,一条水平、一条竖直。这两个方向上的向量变换后方向不变,只是长度乘了一个数。再按"剪切":只剩一条虚线,水平方向。按"旋转":虚线消失了。

交互实验线性变换
A=(1001)A=\begin{pmatrix}\textcolor{#c2452d}{1}&\textcolor{#1f7a4d}{0}\\\textcolor{#c2452d}{0}&\textcolor{#1f7a4d}{1}\end{pmatrix}
行列式(有向面积) 1
特征值 1, 1
猜想

每个线性变换可能有一些"特殊方向",在这些方向上它只做拉伸。如果能找到足够多这样的方向,把它们当作新的基,变换在新基下就只是分别拉伸各个坐标:那是最简单的矩阵,对角矩阵。旋转没有这样的实方向,剪切只有一个。

定义

定义 5.1特征值与特征向量

设 AA 是 n×nn\times n 矩阵。若存在非零向量 vv 和数 λ\lambda 使 Av=λv,Av=\lambda v, 则称 λ\lambda 是 AA 的特征值,vv 是属于 λ\lambda 的特征向量。

"非零"是必要的:A0=λ0A0=\lambda0 对任何 λ\lambda 都成立,不含信息。

推导:怎么找

Av=λv  ⟺  (A−λI)v=0Av=\lambda v\iff(A-\lambda I)v=0。要有非零解 vv,矩阵 A−λIA-\lambda I 必须不可逆,即(见矩阵) det⁡(A−λI)=0.\det(A-\lambda I)=0.

定义 5.2特征多项式

pA(λ)=det⁡(A−λI)p_A(\lambda)=\det(A-\lambda I) 称为 AA 的特征多项式。对 2×22\times2 矩阵, pA(λ)=λ2−(tr⁡A)λ+det⁡A,tr⁡A=a+d.p_A(\lambda)=\lambda^2-(\operatorname{tr}A)\lambda+\det A,\qquad\operatorname{tr}A=a+d.

证明

det⁡(a−λbcd−λ)=(a−λ)(d−λ)−bc=λ2−(a+d)λ+(ad−bc)\det\begin{pmatrix}a-\lambda&b\\c&d-\lambda\end{pmatrix}=(a-\lambda)(d-\lambda)-bc=\lambda^2-(a+d)\lambda+(ad-bc)。

特征值就是特征多项式的根。求出 λ\lambda 后,解线性方程组 (A−λI)v=0(A-\lambda I)v=0 得特征向量。

例 5.1三种典型情形
  • 拉伸 (1.6000.7)\begin{pmatrix}1.6&0\\0&0.7\end{pmatrix}:p=(λ−1.6)(λ−0.7)p=(\lambda-1.6)(\lambda-0.7),特征值 1.6,0.71.6,0.7,特征向量 e1,e2\mathbf e_1,\mathbf e_2。两个独立方向。
  • 剪切 (1101)\begin{pmatrix}1&1\\0&1\end{pmatrix}:p=(λ−1)2p=(\lambda-1)^2,只有 λ=1\lambda=1。(A−I)v=0(A-I)v=0 迫使 v2=0v_2=0,特征向量只有 e1\mathbf e_1 的倍数。一个二重根只给了一个方向。
  • 旋转 (0−110)\begin{pmatrix}0&-1\\1&0\end{pmatrix}:p=λ2+1p=\lambda^2+1,无实根。没有实的特征方向;特征值是 ±i\pm i。

判别式 (tr⁡A)2−4det⁡A(\operatorname{tr}A)^2-4\det A 的符号决定了这三种情形。实验面板里显示的"复数"就是判别式为负。

定理与证明

定理 5.1不同特征值的特征向量线性无关

若 v1,…,vkv_1,\dots,v_k 是属于两两不同特征值 λ1,…,λk\lambda_1,\dots,\lambda_k 的特征向量,则它们线性无关。

证明

对 kk 归纳。k=1k=1 显然(v1≠0v_1\ne0)。设对 k−1k-1 成立。若 c1v1+⋯+ckvk=0c_1v_1+\cdots+c_kv_k=0,两边乘 AA 得 c1λ1v1+⋯+ckλkvk=0c_1\lambda_1v_1+\cdots+c_k\lambda_kv_k=0;原式乘 λk\lambda_k 再相减得 c1(λ1−λk)v1+⋯+ck−1(λk−1−λk)vk−1=0.c_1(\lambda_1-\lambda_k)v_1+\cdots+c_{k-1}(\lambda_{k-1}-\lambda_k)v_{k-1}=0. 由归纳假设各系数为零,而 λi≠λk\lambda_i\ne\lambda_k,故 c1=⋯=ck−1=0c_1=\cdots=c_{k-1}=0,进而 ckvk=0c_kv_k=0,ck=0c_k=0。

定理 5.2对角化

若 n×nn\times n 矩阵 AA 有 nn 个线性无关的特征向量 v1,…,vnv_1,\dots,v_n,令 P=[v1⋯vn]P=[v_1\cdots v_n],D=diag⁡(λ1,…,λn)D=\operatorname{diag}(\lambda_1,\dots,\lambda_n),则 A=PDP−1.A=PDP^{-1}.

证明

APAP 的第 jj 列是 Avj=λjvjAv_j=\lambda_jv_j;PDPD 的第 jj 列是 PP 的第 jj 列乘 λj\lambda_j,也是 λjvj\lambda_jv_j。故 AP=PDAP=PD。PP 的列线性无关,故可逆。

含义:换到以特征向量为基的坐标系,AA 就是对角矩阵 DD。"在特殊方向上只做拉伸"这句直觉,变成了一个精确等式。剪切不能对角化,因为它凑不齐两个独立方向。

常见错误

特征向量必须非零,但特征值可以为 00(此时 AA 有非平凡的核,行列式为零)。另外,二重特征值不保证两个线性无关的特征向量:剪切矩阵只有一个方向。代数重数与几何重数是两回事。

应用

应用矩阵的幂

Ak=PDkP−1A^k=PD^kP^{-1},而 Dk=diag⁡(λ1k,…,λnk)D^k=\operatorname{diag}(\lambda_1^k,\dots,\lambda_n^k)。Fibonacci 数列 Fn+1=Fn+Fn−1F_{n+1}=F_n+F_{n-1} 写成 (Fn+1Fn)=(1110)(FnFn−1)\begin{pmatrix}F_{n+1}\\F_n\end{pmatrix}=\begin{pmatrix}1&1\\1&0\end{pmatrix}\begin{pmatrix}F_n\\F_{n-1}\end{pmatrix},特征值 1±52\frac{1\pm\sqrt5}2,由此得到 Binet 公式。

应用二阶线性微分方程

y′′+py′+qy=0y''+py'+qy=0 化为一阶系统 (yy′)′=(01−q−p)(yy′)\begin{pmatrix}y\\y'\end{pmatrix}'=\begin{pmatrix}0&1\\-q&-p\end{pmatrix}\begin{pmatrix}y\\y'\end{pmatrix},其特征多项式是 λ2+pλ+q\lambda^2+p\lambda+q:正是教科书里的"特征方程"。这就是二阶线性方程一节的主题。

注

特征值回答的问题是:"这个变换在哪些方向上最简单?"对称矩阵永远有一整组正交的特征方向(谱定理),这是主成分分析、振动模态、量子力学可观测量的共同数学基础。

练习

01
求 (2112)\begin{pmatrix}2&1\\1&2\end{pmatrix} 的最大特征值。
02
剪切矩阵 (1101)\begin{pmatrix}1&1\\0&1\end{pmatrix} 的特征向量是?