Summary n 个变量的齐二次多项式,唯一对应一个对称矩阵,交叉项系数要减半。目标是通过可逆线性替换消去交叉项、化成平方和:合同变换保持秩与正负惯性指数,标准形有正交变换、配方法、初等变换三条路。惯性定理指出正、负惯性指数与秩由二次型本身决定,规范形唯一。正定性用特征值全正或顺序主子式全正判定,并延伸到半正定、负定。几何上按特征值符号把二次曲面分成椭球面、单叶与双叶双曲面等,应用见于多元极值的海森矩阵、协方差矩阵与最优化的凸性。
二次型 (Quadratic Form)是 n n n 个变量的 齐二次多项式 (每一项都是二次):
f ( x 1 , x 2 , … , x n ) = ∑ i = 1 n ∑ j = 1 n a i j x i x j f(x_1,x_2,\dots,x_n)=\sum_{i=1}^n\sum_{j=1}^na_{ij}x_ix_j f ( x 1 , x 2 , … , x n ) = i = 1 ∑ n j = 1 ∑ n a ij x i x j
它是中学里 a x 2 + b x y + c y 2 ax^2+bxy+cy^2 a x 2 + b x y + c y 2 向高维的推广。研究目标只有一个:通过 变量替换 ,把交叉项 x i x j x_ix_j x i x j 全部消掉,化成纯 平方和 形式,从而一眼看出它的符号性质和几何形状。
为什么值得专门研究?因为二次型无处不在:二次曲线/曲面的分类、多元函数极值的二阶判别(海森矩阵)、概率论里的协方差矩阵、最优化里的目标函数,背后都是同一套二次型理论。
任何二次型都能唯一地写成 对称矩阵 的形式:
f ( x ⃗ ) = x ⃗ T A x ⃗ , A T = A f(\vec x)=\vec x^TA\vec x,A^T=A f ( x ) = x T A x , A T = A
对应规则:a i i a_{ii} a ii 取 x i 2 x_i^2 x i 2 的系数;a i j = a j i a_{ij}=a_{ji} a ij = a j i 各取 x i x j x_ix_j x i x j 系数的 一半 (因为交叉项被 a i j a_{ij} a ij 和 a j i a_{ji} a j i 平摊)。
矩阵 A A A 称为二次型的 矩阵 ,其秩 r ( A ) r(A) r ( A ) 称为二次型的 秩 。
「系数减半」是新手最容易错的地方。例如 f = x 1 2 + 4 x 1 x 2 + 3 x 2 2 f=x_1^2+4x_1x_2+3x_2^2 f = x 1 2 + 4 x 1 x 2 + 3 x 2 2 ,对角填系数 1 1 1 和 3 3 3 ,交叉项 4 x 1 x 2 4x_1x_2 4 x 1 x 2 拆成 a 12 = a 21 = 2 a_{12}=a_{21}=2 a 12 = a 21 = 2 :
A = ( 1 2 2 3 ) A=\begin{pmatrix}
1 & 2 \\
2 & 3
\end{pmatrix} A = ( 1 2 2 3 ) 坚持用 对称矩阵 表示,是因为只有对称矩阵才一定能正交对角化,整套理论才走得通。
做可逆线性替换 x ⃗ = C y ⃗ \vec x=C\vec y x = C y (C C C 可逆),二次型变为:
f = x ⃗ T A x ⃗ = ( C y ⃗ ) T A ( C y ⃗ ) = y ⃗ T ( C T A C ) y ⃗ f=\vec x^TA\vec x=(C\vec y)^TA(C\vec y)=\vec y^T(C^TAC)\vec y f = x T A x = ( C y ) T A ( C y ) = y T ( C T A C ) y
于是新矩阵是 C T A C C^TAC C T A C 。若存在可逆 C C C 使 B = C T A C B=C^TAC B = C T A C ,称 A , B A,B A , B 合同 (Congruent),记 A ≃ B A\simeq B A ≃ B 。合同变换对应「换一组坐标看同一个二次型」,它保持二次型的本质(秩、正负惯性指数)不变。
关系 定义 几何意义 相似 P − 1 A P = B P^{-1}AP=B P − 1 A P = B 同一线性变换在不同基下 合同 C T A C = B C^TAC=B C T A C = B 同一二次型在不同坐标下
两者一般不同,但有一个重要交汇:当 C C C 是 正交矩阵 时 C T = C − 1 C^T=C^{-1} C T = C − 1 ,相似与合同 同时成立 。这正是用正交变换化标准形既保特征值(相似)又保二次型(合同)的原因。
只含平方项、没有交叉项的二次型称为 标准形 (Canonical Form):
f = d 1 y 1 2 + d 2 y 2 2 + ⋯ + d n y n 2 f=d_1y_1^2+d_2y_2^2+\dots+d_ny_n^2 f = d 1 y 1 2 + d 2 y 2 2 + ⋯ + d n y n 2
对应矩阵是对角阵 diag ( d 1 , … , d n ) \operatorname{diag}(d_1,\dots,d_n) diag ( d 1 , … , d n ) 。化标准形,就是给二次型的矩阵找一个合同的对角阵。
方法 思路 正交变换法 求 A A A 的特征值 λ i \lambda_i λ i 与正交矩阵 Q Q Q ,令 x ⃗ = Q y ⃗ \vec x=Q\vec y x = Q y ,得 f = λ 1 y 1 2 + ⋯ + λ n y n 2 f=\lambda_1y_1^2+\dots+\lambda_ny_n^2 f = λ 1 y 1 2 + ⋯ + λ n y n 2 配方法(拉格朗日) 逐个把含某变量的项凑成完全平方,再换元,反复直到无交叉项 初等变换法 对 ( A ∣ E ) (A\mid E) ( A ∣ E ) 做 成对的 行列初等变换,左边化对角时右边即变换矩阵
三种方法得到的对角元 不唯一 (系数大小可以不同),但下文的惯性指数是唯一的。
正交变换法 最特殊:因为 Q Q Q 正交,它在化简的同时 保持了几何形状 (不拉伸、不扭曲),所以对角元恰好是 特征值 ,对应各主轴方向的真实缩放。需要研究二次曲面形状时用它;只想快速判断正负性时,配方法 往往更省力。
例(有平方项):把 f = x 1 2 + 2 x 2 2 + 5 x 3 2 + 2 x 1 x 2 + 2 x 1 x 3 + 6 x 2 x 3 f=x_1^2+2x_2^2+5x_3^2+2x_1x_2+2x_1x_3+6x_2x_3 f = x 1 2 + 2 x 2 2 + 5 x 3 2 + 2 x 1 x 2 + 2 x 1 x 3 + 6 x 2 x 3 化标准形。思路是 先收齐含 x 1 x_1 x 1 的所有项凑成完全平方 ,再处理剩下的。含 x 1 x_1 x 1 的是 x 1 2 + 2 x 1 x 2 + 2 x 1 x 3 x_1^2+2x_1x_2+2x_1x_3 x 1 2 + 2 x 1 x 2 + 2 x 1 x 3 ,配成 ( x 1 + x 2 + x 3 ) 2 (x_1+x_2+x_3)^2 ( x 1 + x 2 + x 3 ) 2 会多出 x 2 2 + x 3 2 + 2 x 2 x 3 x_2^2+x_3^2+2x_2x_3 x 2 2 + x 3 2 + 2 x 2 x 3 ,扣掉:
f = ( x 1 + x 2 + x 3 ) 2 + ( 2 x 2 2 − x 2 2 ) + ( 5 x 3 2 − x 3 2 ) + ( 6 x 2 x 3 − 2 x 2 x 3 ) = ( x 1 + x 2 + x 3 ) 2 + x 2 2 + 4 x 3 2 + 4 x 2 x 3 f=(x_1+x_2+x_3)^2+(2x_2^2-x_2^2)+(5x_3^2-x_3^2)+(6x_2x_3-2x_2x_3)=(x_1+x_2+x_3)^2+x_2^2+4x_3^2+4x_2x_3 f = ( x 1 + x 2 + x 3 ) 2 + ( 2 x 2 2 − x 2 2 ) + ( 5 x 3 2 − x 3 2 ) + ( 6 x 2 x 3 − 2 x 2 x 3 ) = ( x 1 + x 2 + x 3 ) 2 + x 2 2 + 4 x 3 2 + 4 x 2 x 3
剩下的 x 2 2 + 4 x 2 x 3 + 4 x 3 2 x_2^2+4x_2x_3+4x_3^2 x 2 2 + 4 x 2 x 3 + 4 x 3 2 再对 x 2 x_2 x 2 配方,正好是 ( x 2 + 2 x 3 ) 2 (x_2+2x_3)^2 ( x 2 + 2 x 3 ) 2 :
f = ( x 1 + x 2 + x 3 ) 2 + ( x 2 + 2 x 3 ) 2 f=(x_1+x_2+x_3)^2+(x_2+2x_3)^2 f = ( x 1 + x 2 + x 3 ) 2 + ( x 2 + 2 x 3 ) 2
令 y 1 = x 1 + x 2 + x 3 y_1=x_1+x_2+x_3 y 1 = x 1 + x 2 + x 3 、y 2 = x 2 + 2 x 3 y_2=x_2+2x_3 y 2 = x 2 + 2 x 3 、y 3 = x 3 y_3=x_3 y 3 = x 3 (这是可逆线性替换),得标准形 f = y 1 2 + y 2 2 f=y_1^2+y_2^2 f = y 1 2 + y 2 2 。秩为 2 2 2 ,正惯性指数 p = 2 p=2 p = 2 ,负惯性指数 q = 0 q=0 q = 0 。
例(无平方项):把 f = 2 x 1 x 2 + 2 x 1 x 3 − 2 x 2 x 3 f=2x_1x_2+2x_1x_3-2x_2x_3 f = 2 x 1 x 2 + 2 x 1 x 3 − 2 x 2 x 3 化标准形。全是交叉项、没有平方项,配方无从下手,先做一个 预备替换 造出平方项:令 x 1 = y 1 + y 2 x_1=y_1+y_2 x 1 = y 1 + y 2 、x 2 = y 1 − y 2 x_2=y_1-y_2 x 2 = y 1 − y 2 、x 3 = y 3 x_3=y_3 x 3 = y 3 。代入,2 x 1 x 2 = 2 ( y 1 + y 2 ) ( y 1 − y 2 ) = 2 y 1 2 − 2 y 2 2 2x_1x_2=2(y_1+y_2)(y_1-y_2)=2y_1^2-2y_2^2 2 x 1 x 2 = 2 ( y 1 + y 2 ) ( y 1 − y 2 ) = 2 y 1 2 − 2 y 2 2 ,其余项 2 x 1 x 3 − 2 x 2 x 3 = 2 x 3 ( x 1 − x 2 ) = 2 y 3 ⋅ 2 y 2 = 4 y 2 y 3 2x_1x_3-2x_2x_3=2x_3(x_1-x_2)=2y_3\cdot2y_2=4y_2y_3 2 x 1 x 3 − 2 x 2 x 3 = 2 x 3 ( x 1 − x 2 ) = 2 y 3 ⋅ 2 y 2 = 4 y 2 y 3 ,于是:
f = 2 y 1 2 − 2 y 2 2 + 4 y 2 y 3 f=2y_1^2-2y_2^2+4y_2y_3 f = 2 y 1 2 − 2 y 2 2 + 4 y 2 y 3
现在有平方项了,对 y 2 y_2 y 2 配方:− 2 y 2 2 + 4 y 2 y 3 = − 2 ( y 2 2 − 2 y 2 y 3 ) = − 2 ( y 2 − y 3 ) 2 + 2 y 3 2 -2y_2^2+4y_2y_3=-2(y_2^2-2y_2y_3)=-2(y_2-y_3)^2+2y_3^2 − 2 y 2 2 + 4 y 2 y 3 = − 2 ( y 2 2 − 2 y 2 y 3 ) = − 2 ( y 2 − y 3 ) 2 + 2 y 3 2 ,故:
f = 2 y 1 2 − 2 ( y 2 − y 3 ) 2 + 2 y 3 2 f=2y_1^2-2(y_2-y_3)^2+2y_3^2 f = 2 y 1 2 − 2 ( y 2 − y 3 ) 2 + 2 y 3 2
再令 z 1 = y 1 z_1=y_1 z 1 = y 1 、z 2 = y 2 − y 3 z_2=y_2-y_3 z 2 = y 2 − y 3 、z 3 = y 3 z_3=y_3 z 3 = y 3 ,得标准形 f = 2 z 1 2 − 2 z 2 2 + 2 z 3 2 f=2z_1^2-2z_2^2+2z_3^2 f = 2 z 1 2 − 2 z 2 2 + 2 z 3 2 。秩为 3 3 3 ,正惯性指数 p = 2 p=2 p = 2 、负惯性指数 q = 1 q=1 q = 1 ,符号差 p − q = 1 p-q=1 p − q = 1 。「无平方项先凑 x i = y i + y j , x j = y i − y j x_i=y_i+y_j,x_j=y_i-y_j x i = y i + y j , x j = y i − y j 」是这类题的固定起手。
把标准形里每个非零平方项的系数再通过伸缩化成 ± 1 \pm 1 ± 1 ,得到 规范形 (Normal Form):
f = y 1 2 + ⋯ + y p 2 − y p + 1 2 − ⋯ − y p + q 2 f=y_1^2+\dots+y_p^2-y_{p+1}^2-\dots-y_{p+q}^2 f = y 1 2 + ⋯ + y p 2 − y p + 1 2 − ⋯ − y p + q 2
只剩 + 1 +1 + 1 、− 1 -1 − 1 、0 0 0 三种系数,是最「素」的形态。
一个二次型用 不同 的可逆线性变换化出的标准形,系数可以各不相同,但其中:
正系数的个数 p p p —— 正惯性指数 ;
负系数的个数 q q q —— 负惯性指数 ;
非零系数总数 r = p + q r=p+q r = p + q —— 二次型的 秩 ;
这三个数 完全由二次型本身决定 ,与采用何种变换无关。所以规范形 唯一 。p − q p-q p − q 称为 符号差 。
例:用正交变换把 f = 2 x 1 2 + 2 x 2 2 − 2 x 1 x 2 f=2x_1^2+2x_2^2-2x_1x_2 f = 2 x 1 2 + 2 x 2 2 − 2 x 1 x 2 化标准形并写出变换矩阵。二次型矩阵(交叉项系数减半):
A = ( 2 − 1 − 1 2 ) A=\begin{pmatrix}
2 & -1 \\
-1 & 2
\end{pmatrix} A = ( 2 − 1 − 1 2 )
特征方程 ∣ A − λ E ∣ = ( 2 − λ ) 2 − 1 = 0 |A-\lambda E|=(2-\lambda)^2-1=0 ∣ A − λ E ∣ = ( 2 − λ ) 2 − 1 = 0 给 λ 1 = 1 , λ 2 = 3 \lambda_1=1,\lambda_2=3 λ 1 = 1 , λ 2 = 3 。对 λ 1 = 1 \lambda_1=1 λ 1 = 1 :( A − E ) x ⃗ = ( 1 − 1 − 1 1 ) x ⃗ = 0 ⃗ (A-E)\vec x=\begin{pmatrix}1&-1\\-1&1\end{pmatrix}\vec x=\vec0 ( A − E ) x = ( 1 − 1 − 1 1 ) x = 0 ,得 ξ ⃗ 1 = ( 1 , 1 ) T \vec\xi_1=(1,1)^T ξ 1 = ( 1 , 1 ) T ;对 λ 2 = 3 \lambda_2=3 λ 2 = 3 :( − 1 − 1 − 1 − 1 ) x ⃗ = 0 ⃗ \begin{pmatrix}-1&-1\\-1&-1\end{pmatrix}\vec x=\vec0 ( − 1 − 1 − 1 − 1 ) x = 0 ,得 ξ ⃗ 2 = ( 1 , − 1 ) T \vec\xi_2=(1,-1)^T ξ 2 = ( 1 , − 1 ) T 。两特征值互异,特征向量已自动正交,只需单位化:
Q = 1 2 ( 1 1 1 − 1 ) Q=\frac{1}{\sqrt2}\begin{pmatrix}
1 & 1 \\
1 & -1
\end{pmatrix} Q = 2 1 ( 1 1 1 − 1 )
令 x ⃗ = Q y ⃗ \vec x=Q\vec y x = Q y ,则 f = λ 1 y 1 2 + λ 2 y 2 2 = y 1 2 + 3 y 2 2 f=\lambda_1y_1^2+\lambda_2y_2^2=y_1^2+3y_2^2 f = λ 1 y 1 2 + λ 2 y 2 2 = y 1 2 + 3 y 2 2 。对角元恰是特征值,Q Q Q 的列就是主轴方向。因为 Q Q Q 正交,这个化简同时保住了图形:f = 1 f=1 f = 1 是一个椭圆,长短半轴比由 3 : 1 \sqrt{3}:1 3 : 1 决定。
惯性定理是二次型的「身份证」:判断两个实二次型能否互相合同,只看它们的 秩 r r r 和正惯性指数 p p p 是否相同即可,不必真的去找变换矩阵。
设 f ( x ⃗ ) = x ⃗ T A x ⃗ f(\vec x)=\vec x^TA\vec x f ( x ) = x T A x (A A A 实对称),按它对所有非零 x ⃗ \vec x x 取值的符号分类:
类型 对所有 x ⃗ ≠ 0 ⃗ \vec x\ne\vec 0 x = 0 正定 (Positive Definite)f ( x ⃗ ) > 0 f(\vec x)>0 f ( x ) > 0 负定 (Negative Definite)f ( x ⃗ ) < 0 f(\vec x)<0 f ( x ) < 0 半正定 f ( x ⃗ ) ≥ 0 f(\vec x)\ge 0 f ( x ) ≥ 0 半负定 f ( x ⃗ ) ≤ 0 f(\vec x)\le 0 f ( x ) ≤ 0 不定 既能取正也能取负
正定的几何画面:图像是一个开口向上、最低点在原点的「碗」,从任何方向离开原点函数值都上升。
下列条件相互等价(A A A 为 n n n 阶实对称矩阵):
f f f 正定;
A A A 的 所有特征值 都 > 0 >0 > 0 ;
A A A 的 所有顺序主子式 都 > 0 >0 > 0 ;
正惯性指数 p = n p=n p = n ;
存在可逆矩阵 P P P 使 A = P T P A=P^TP A = P T P ;
A A A 合同于单位阵 E E E 。
挑哪条判定看场景:手头有特征值就用「特征值全正」;只有矩阵 A A A 就用 顺序主子式判别法(赫尔维茨判据) 最省事 ——:
a 11 > 0 , ∣ a 11 a 12 a 21 a 22 ∣ > 0 , … , ∣ A ∣ > 0 a_{11}>0,\begin{vmatrix}
a_{11} & a_{12} \\
a_{21} & a_{22}
\end{vmatrix}>0,\dots,|A|>0 a 11 > 0 , a 11 a 21 a 12 a 22 > 0 , … , ∣ A ∣ > 0 依次算 1 1 1 到 n n n 阶顺序主子式,全部 > 0 >0 > 0 即正定。注意是「顺序」主子式(从左上角起的那一串),不是任意主子式。
例(顺序主子式法):判定 f = 2 x 1 2 + 5 x 2 2 + 5 x 3 2 + 4 x 1 x 2 − 4 x 1 x 3 − 8 x 2 x 3 f=2x_1^2+5x_2^2+5x_3^2+4x_1x_2-4x_1x_3-8x_2x_3 f = 2 x 1 2 + 5 x 2 2 + 5 x 3 2 + 4 x 1 x 2 − 4 x 1 x 3 − 8 x 2 x 3 是否正定。先写矩阵(交叉项减半):
A = ( 2 2 − 2 2 5 − 4 − 2 − 4 5 ) A=\begin{pmatrix}
2 & 2 & -2 \\
2 & 5 & -4 \\
-2 & -4 & 5
\end{pmatrix} A = 2 2 − 2 2 5 − 4 − 2 − 4 5
依次算三个顺序主子式:
D 1 = 2 > 0 , D 2 = ∣ 2 2 2 5 ∣ = 10 − 4 = 6 > 0 , D 3 = ∣ A ∣ = 2 ( 25 − 16 ) − 2 ( 10 − 8 ) + ( − 2 ) ( − 8 + 10 ) = 18 − 4 − 4 = 10 > 0 D_1=2>0,D_2=\begin{vmatrix}
2 & 2 \\
2 & 5
\end{vmatrix}=10-4=6>0,D_3=|A|=2(25-16)-2(10-8)+(-2)(-8+10)=18-4-4=10>0 D 1 = 2 > 0 , D 2 = 2 2 2 5 = 10 − 4 = 6 > 0 , D 3 = ∣ A ∣ = 2 ( 25 − 16 ) − 2 ( 10 − 8 ) + ( − 2 ) ( − 8 + 10 ) = 18 − 4 − 4 = 10 > 0
三个全大于 0 0 0 ,故 f f f 正定 。只用矩阵本身、不必求特征值,这是手算判定正定最快的路子。
例(特征值法):判定 A = ( 3 1 1 3 ) A=\begin{pmatrix}3&1\\1&3\end{pmatrix} A = ( 3 1 1 3 ) 对应的二次型是否正定。特征值 λ = 3 ± 1 = 2 , 4 \lambda=3\pm1=2,4 λ = 3 ± 1 = 2 , 4 ,全为正,故正定。对照顺序主子式法:D 1 = 3 > 0 D_1=3>0 D 1 = 3 > 0 、D 2 = 9 − 1 = 8 > 0 D_2=9-1=8>0 D 2 = 9 − 1 = 8 > 0 ,同样判定正定,两法结论一致。若手头已经算出了特征值(如做完正交对角化),直接看符号最省事;只有矩阵时用顺序主子式。
半正定 :所有特征值 ≥ 0 \ge 0 ≥ 0 ,等价于所有 主子式 (含非顺序的)≥ 0 \ge 0 ≥ 0 (这里顺序主子式 ≥ 0 \ge 0 ≥ 0 不够,是常见陷阱)。
负定 :− A -A − A 正定 ⟺ \iff ⟺ 所有特征值 < 0 <0 < 0 ⟺ \iff ⟺ 顺序主子式 正负交替 ,奇数阶 < 0 <0 < 0 、偶数阶 > 0 >0 > 0 。
二次型 x ⃗ T A x ⃗ = 1 \vec x^TA\vec x=1 x T A x = 1 在不同正定性下,描绘出不同的二次曲线/曲面。正交变换把它转到主轴坐标系后,看特征值(即标准形系数)的符号就能定形状:
特征值符号(n = 3 n=3 n = 3 ) 曲面类型 全为正 椭球面(正定,封闭的「蛋」) 两正一负 单叶双曲面 一正两负 双叶双曲面 有零特征值 退化为柱面、抛物面等
二维同理:两个特征值同号得 椭圆 ,异号得 双曲线 ,有一个为零则退化。特征值的 绝对值 决定各主轴的长短,符号 决定曲线/曲面的开闭与朝向 —— 这就是「特征值即各方向缩放因子」在几何上的直接体现。
多元极值 :把 多元微积分 里的 海森矩阵 当成二次型矩阵,正定 ⟹ \implies ⟹ 极小值、负定 ⟹ \implies ⟹ 极大值、不定 ⟹ \implies ⟹ 鞍点。
二次曲线/曲面分类 :如上,靠特征值符号判别椭圆、双曲、抛物等。
概率论 :随机向量的 协方差矩阵 必半正定,正定则对应「没有退化方向」的非奇异分布。
最优化 :目标函数的二次近似正定时,问题是凸的,局部极小即全局极小。