MathLabs

代数学

双线性型与二次型

对两个向量变量分别线性、或对一个变量为二次的函数,用于定义长度和角度。

直观用两个向量同时度量形状

熟悉的点积 u⋅v\mathbf{u}\cdot\mathbf{v} 接受两个向量并返回一个数,对每个变量都是线性的;双线性型正是把这个想法加以推广,B(u,v)=u⊤AvB(\mathbf{u},\mathbf{v}) = \mathbf{u}^\top A \mathbf{v},允许矩阵 AA 在组合向量之前对其进行扭曲和加权。把同一个向量代入两个位置就得到二次型 Q(x)=x⊤AxQ(\mathbf{x}) = \mathbf{x}^\top A \mathbf{x},其水平集 Q(x)=cQ(\mathbf{x})=c 依据 AA 的形状描绘出椭圆、双曲线或退化直线——正是中学几何中的圆锥曲线,如今推广到任意维度。

可旋转的三维曲面,展示二次型呈现为碗形、鞍形或山脊形。
双变量二次型的曲面 z=Q(x,y)z=Q(x,y):碗形表示正定,鞍形表示不定,山脊形表示存在退化方向。

中学从点积到一般的双线性型

定义: 双线性型与二次型

Rn\mathbb{R}^n 上的双线性型是对每个自变量分别线性的函数 B(u,v)B(\mathbf{u},\mathbf{v});在坐标下它总能用 n×nn \times n 矩阵 AA 写成 B(u,v)=u⊤AvB(\mathbf{u},\mathbf{v}) = \mathbf{u}^\top A \mathbf{v}。当 A⊤=AA^\top = A 时该形式对称,其伴随的二次型为 Q(x)=x⊤AxQ(\mathbf{x}) = \mathbf{x}^\top A \mathbf{x}。反过来,任何二次型都通过极化恒等式 B(u,v)=12(Q(u+v)−Q(u)−Q(v))B(\mathbf{u},\mathbf{v}) = \tfrac{1}{2}(Q(\mathbf{u}+\mathbf{v}) - Q(\mathbf{u}) - Q(\mathbf{v})) 唯一确定其对称双线性型,因此我们总可取 AA 为对称矩阵。

B(u,v)=u⊤Av,Q(x)=x⊤AxB(\mathbf{u},\mathbf{v}) = \mathbf{u}^\top A \mathbf{v}, \qquad Q(\mathbf{x}) = \mathbf{x}^\top A \mathbf{x}

通过 x=Py\mathbf{x} = P\mathbf{y}(PP 可逆)作基变换时,矩阵 AA 不是被相似变换 P−1APP^{-1}AP 替换,而是被合同变换 A′=P⊤APA' = P^\top A P 替换。在使 A′A' 为对角矩阵的基下,二次型化为纯粹的平方和 d1y12+⋯+dnyn2d_1 y_1^2 + \cdots + d_n y_n^2,正、负、零系数的个数构成惯性指数 (n+,n−,n0)(n_+, n_-, n_0)。

A′=P⊤AP,Q(Py)=d1y12+⋯+dnyn2A' = P^\top A P, \qquad Q(P\mathbf{y}) = d_1 y_1^2 + \cdots + d_n y_n^2
按惯性指数对实二次型分类
类型惯性指数 (n+,n−,n0)(n_+, n_-, n_0) 与几何
正定(当 x≠0\mathbf{x}\neq\mathbf{0} 时 Q(x)>0Q(\mathbf{x})>0)(n,0,0)(n,0,0);水平集为椭球面,图像为开口向上的碗
负定(当 x≠0\mathbf{x}\neq\mathbf{0} 时 Q(x)<0Q(\mathbf{x})<0)(0,n,0)(0,n,0);水平集为椭球面,图像为开口向下的穹顶
不定(同时取正负号)n+>0n_+>0 且 n−>0n_->0;水平集为双曲面,图像为马鞍面
半定 / 退化n0>0n_0>0;存在使 QQ 为零的平坦方向,水平集变成柱面

大学对角化与西尔维斯特惯性定理

每个实二次型 Q(x)=x⊤AxQ(\mathbf{x}) = \mathbf{x}^\top A \mathbf{x} 都可通过可逆线性变量代换 x=Py\mathbf{x}=P\mathbf{y} 化为纯对角形 d1y12+⋯+dnyn2d_1 y_1^2 + \cdots + d_n y_n^2。

为什么成立?

对角化二次型甚至不需要特征值:把中学熟悉的配方法逐个变量依次应用,就能把所有交叉项 xixjx_i x_j 全部吸收到平方项中。

证明

我们对变量个数 nn 用归纳法证明;n=1n=1 时已经是单项 a11x12a_{11}x_1^2。

当 n>1n > 1 时,若整个二次型为零则无需证明。否则,先设某个对角系数非零;重新编号变量后可设 a11≠0a_{11} \neq 0。把所有含 x1x_1 的项归并起来,得到 a11x12+2x1(a12x2+⋯+a1nxn)a_{11}x_1^2 + 2x_1(a_{12}x_2 + \cdots + a_{1n}x_n) 加上一个仅含 x2,…,xnx_2,\dots,x_n 的二次型。

令 y1=x1+a12a11x2+⋯+a1na11xny_1 = x_1 + \frac{a_{12}}{a_{11}}x_2 + \cdots + \frac{a_{1n}}{a_{11}}x_n,并对 k≥2k \ge 2 令 yk=xky_k = x_k,这是一个可逆线性变量代换,且 a11y12a_{11}y_1^2 与所有含 x1x_1 的项完全匹配,只差一个仅依赖于 y2,…,yny_2,\dots,y_n 的余项;因此 Q=a11y12+Q′(y2,…,yn)Q = a_{11}y_1^2 + Q'(y_2,\dots,y_n)。

反之,若所有对角元皆为零(对所有 ii 有 aii=0a_{ii}=0),则取一个非零交叉项 2a12x1x22a_{12}x_1 x_2,其中 a12≠0a_{12} \neq 0。可逆代换 x1=u1+u2x_1 = u_1 + u_2, x2=u1−u2x_2 = u_1 - u_2 把 2a12x1x22a_{12}x_1 x_2 变成 2a12(u12−u22)2a_{12}(u_1^2 - u_2^2),从而产生非零对角系数,归结为前一种情形。对 Q′(y2,…,yn)Q'(y_2,\dots,y_n) 应用归纳假设即完成对角化。

无论使用哪种可逆基变换将 Rn\mathbb{R}^n 上的实二次型 QQ 对角化,正系数的个数 n+n_+、负系数的个数 n−n_- 以及零系数的个数 n0n_0 总是相同的;三元组 (n+,n−,n0)(n_+, n_-, n_0) 是 QQ 的内蕴不变量。

为什么成立?

切换到新坐标系可以拉伸坐标轴并改变对角系数 did_i 各自的大小,但绝不可能在不经过平坦方向的情况下把向上弯曲的方向变成向下弯曲的方向——因此向上、向下与平坦坐标轴的数目被永久锁定。

证明

设 QQ 在两组基 {e1,…,en}\{\mathbf{e}_1,\dots,\mathbf{e}_n\} 与 {f1,…,fn}\{\mathbf{f}_1,\dots,\mathbf{f}_n\} 下被对角化,正、负、零系数的个数在第一组基下为 (n+,n−,n0)(n_+, n_-, n_0),在第二组基下为 (p+,p−,p0)(p_+, p_-, p_0)。对每组基排序,使正系数在前,负系数居中,零系数在后。

反设 n+>p+n_+ > p_+。令 V+=span(e1,…,en+)V_+ = \mathrm{span}(\mathbf{e}_1,\dots,\mathbf{e}_{n_+}),这是一个 n+n_+ 维子空间,对任意非零 x∈V+\mathbf{x} \in V_+ 都有 Q(x)>0Q(\mathbf{x}) > 0(因为在 V+V_+ 上只有 e\mathbf{e} 展开中的正平方项起作用)。

类似地,令 W−=span(fp++1,…,fn)W_- = \mathrm{span}(\mathbf{f}_{p_+ + 1},\dots,\mathbf{f}_n),这是一个 n−p+n - p_+ 维子空间,对任意 x∈W−\mathbf{x} \in W_- 都有 Q(x)≤0Q(\mathbf{x}) \le 0(因为在 W−W_- 上只有 f\mathbf{f} 展开中的负平方项和零项起作用)。

现在在 Rn\mathbb{R}^n 中计算维数:由于 n+>p+n_+ > p_+,有 dim⁡V++dim⁡W−=n++(n−p+)>n\dim V_+ + \dim W_- = n_+ + (n - p_+) > n。由子空间维数公式,维数之和大于 nn 的两个子空间不可能只有平凡交;因此存在非零向量 v∈V+∩W−\mathbf{v} \in V_+ \cap W_-。

由于 v∈V+\mathbf{v} \in V_+ 且 v≠0\mathbf{v} \neq \mathbf{0},必有 Q(v)>0Q(\mathbf{v}) > 0;又因 v∈W−\mathbf{v} \in W_-,同时必有 Q(v)≤0Q(\mathbf{v}) \le 0,直接矛盾。因此 n+≤p+n_+ \le p_+,由两组基的对称性又有 p+≤n+p_+ \le n_+,故 n+=p+n_+ = p_+。对 −Q-Q 应用完全相同的论证得 n−=p−n_- = p_-,最后 n0=n−n+−n−=p0n_0 = n - n_+ - n_- = p_0。

大学实际应用与典型例题

二次型及其惯性指数决定了多变量函数的临界点是极小、极大还是鞍点,并在狭义相对论中编码了时空的因果结构。

例题: 通过海森二次型作二阶导数判别

通过求二次型的惯性指数,判别 f(x,y)=x2+4xy+y2f(x,y) = x^2 + 4xy + y^2 的临界点 (0,0)(0,0) 的类型。

解答

在 (0,0)(0,0) 附近,该函数本身就是纯二次型 Q(x,y)=x2+4xy+y2Q(x,y) = x^2 + 4xy + y^2,其对称矩阵为 A=(1221)A = \begin{pmatrix} 1 & 2 \\ 2 & 1 \end{pmatrix}(注意非对角元是 xyxy 系数的一半)。

对 xx 配方:x2+4xy+y2=(x+2y)2−4y2+y2=(x+2y)2−3y2x^2 + 4xy + y^2 = (x + 2y)^2 - 4y^2 + y^2 = (x + 2y)^2 - 3y^2。在新坐标 u=x+2yu = x + 2y, v=yv = y 下,它化为 u2−3v2u^2 - 3v^2。

这里有一个正平方项(+u2+u^2)和一个负平方项(−3v2-3v^2),故惯性指数为 (1,1,0)(1,1,0)——该二次型是不定的。沿 v=0v=0 函数像 +u2+u^2 那样向上弯曲,而沿 u=0u=0 则像 −3v2-3v^2 那样向下弯曲,因此 (0,0)(0,0) 是鞍点,既非局部极小也非局部极大。

例题: 闵可夫斯基时空间隔与因果惯性指数

在狭义相对论中(取 c=1c=1),事件 (t,x,y,z)(t,x,y,z) 与原点之间的时空间隔是二次型 s2=t2−x2−y2−z2s^2 = t^2 - x^2 - y^2 - z^2。求其惯性指数,并利用西尔维斯特惯性定理解释为何所有惯性观测者对两个事件能否发生因果联系意见一致。

解答

二次型 s2=t2−x2−y2−z2s^2 = t^2 - x^2 - y^2 - z^2 已经是对角形,矩阵为 η=diag(1,−1,−1,−1)\eta = \mathrm{diag}(1,-1,-1,-1):一个 +1+1 和三个 −1-1,故其惯性指数为 (1,3,0)(1,3,0)——是不定的,而不像欧氏距离那样正定。

满足 s2>0s^2 > 0 的事件(类时间隔)位于光锥内部,可通过低于光速的信号相连;满足 s2<0s^2 < 0 的事件(类空间隔)位于光锥外部,彼此无法影响;满足 s2=0s^2 = 0 的事件(类光)只能由光线连接。

从一个惯性观测者切换到另一个惯性观测者是保持二次型 s2s^2 不变的线性坐标变换(洛伦兹变换),而根据西尔维斯特惯性定理,惯性指数 (1,3,0)(1,3,0)——一个时间维与三个空间维——不能被任何可逆坐标变换改变。因此 s2s^2 的符号是不变量,所有观测者对哪些事件对是类时、类空或类光完全一致。

与二次型 Q(x,y)=3x2−6xy+5y2Q(x,y) = 3x^2 - 6xy + 5y^2 相伴的对称矩阵 AA 是什么?

Q(x,y)=x2−2xy+y2Q(x,y) = x^2 - 2xy + y^2 的惯性指数 (n+,n−,n0)(n_+, n_-, n_0) 是什么?

当一个实二次型在两组不同基下被对角化时,西尔维斯特惯性定理断言什么?

若光滑函数 f(x,y)f(x,y) 在临界点处的海森矩阵二次型具有惯性指数 (2,0,0)(2,0,0),则该临界点是什么类型?

参考文献

  1. Roger A. Horn, Charles R. Johnson (2012). Matrix Analysis (2nd ed.)
  2. Gilbert Strang (2016). Introduction to Linear Algebra