黑塞矩阵与极值理论深入

多元函数的泰勒展开

定理1

设函数 f(x,y)f(x, y) 在点 (x0,y0)(x_0, y_0) 的某邻域内具有直到 n+1n+1 阶连续偏导数,则对该邻域内任一点 (x0+h,y0+k)(x_0+h, y_0+k),有

f(x0+h,y0+k)=f(x0,y0)+(hx+ky)f(x0,y0)+12!(hx+ky)2f(x0,y0)++1n!(hx+ky)nf(x0,y0)+Rnf(x_0+h, y_0+k) = f(x_0,y_0) + \left(h\frac{\partial}{\partial x} + k\frac{\partial}{\partial y}\right)f(x_0,y_0) + \frac{1}{2!}\left(h\frac{\partial}{\partial x} + k\frac{\partial}{\partial y}\right)^2 f(x_0,y_0) + \cdots + \frac{1}{n!}\left(h\frac{\partial}{\partial x} + k\frac{\partial}{\partial y}\right)^n f(x_0,y_0) + R_n

其中二阶项展开为: 12(h2fxx+2hkfxy+k2fyy)=12(h,k)Hf(hk)\frac{1}{2}(h^2 f_{xx} + 2hk f_{xy} + k^2 f_{yy}) = \frac{1}{2}(h, k) H_f \begin{pmatrix} h \\ k \end{pmatrix}

这就是黑塞矩阵的来源——泰勒展开二阶项的系数矩阵。

几何解释
推论
证明
符号说明

黑塞矩阵与二次型

二次型与黑塞矩阵

在驻点 (x0,y0)(x_0, y_0) 处,一阶项为零,泰勒展开简化为:

f(x0+h,y0+k)f(x0,y0)=12(h,k)Hf(hk)+o(ρ2)f(x_0+h, y_0+k) - f(x_0,y_0) = \frac{1}{2}(h, k) H_f \begin{pmatrix} h \\ k \end{pmatrix} + o(\rho^2)

其中 Hf=(fxxfxyfyxfyy)H_f = \begin{pmatrix} f_{xx} & f_{xy} \\ f_{yx} & f_{yy} \end{pmatrix} 是黑塞矩阵。

函数增量的符号由二次型 Q(h,k)=(h,k)Hf(h,k)TQ(h,k) = (h,k)H_f(h,k)^T 决定:

  • Q>0Q > 0(正定)⟹ 极小值
  • Q<0Q < 0(负定)⟹ 极大值
  • QQ 可正可负(不定)⟹ 鞍点
  • Q0Q \geq 0Q0Q \leq 0(半定)⟹ 无法确定

正定性的判定方法

定理2

对于 nn 阶实对称矩阵 AA,以下条件等价:

  1. 正定

    • 对任意非零向量 x\vec{x}xTAx>0\vec{x}^T A \vec{x} > 0
    • 所有特征值 λi>0\lambda_i > 0
    • 所有顺序主子式 Dk>0D_k > 0k=1,,nk=1,\ldots,n
    • 存在可逆矩阵 CC 使 A=CTCA = C^T C
  2. 负定AA 负定 ⟺ A-A 正定 ⟺ 奇数阶顺序主子式 <0<0,偶数阶 >0>0

  3. 不定:既有正特征值又有负特征值

  4. 半正定/半负定:所有特征值 0\geq 0 / 0\leq 0,且至少有一个为零

几何解释
推论
证明
符号说明
二元函数的判别式与正定性

对于 2×22 \times 2 黑塞矩阵 H=(ABBC)H = \begin{pmatrix} A & B \\ B & C \end{pmatrix}

  • 一阶顺序主子式 D1=AD_1 = A
  • 二阶顺序主子式 D2=ACB2=ΔD_2 = AC - B^2 = \Delta
  • 正定 ⟺ A>0A > 0Δ>0\Delta > 0 ⟹ 极小值
  • 负定 ⟺ A<0A < 0Δ>0\Delta > 0 ⟹ 极大值
  • 不定 ⟺ Δ<0\Delta < 0 ⟹ 鞍点
  • 半定 ⟺ Δ=0\Delta = 0 ⟹ 无法确定

约束极值的二阶条件

定理3

对于条件极值问题 minf(x)\min f(x) s.t. gi(x)=0g_i(x)=0i=1,,mi=1,\ldots,m),拉格朗日函数 L=f+λigiL = f + \sum \lambda_i g_i

在驻点 xx^* 处,定义加边黑塞矩阵(bordered Hessian):

HB=(0JJTHL)H_B = \begin{pmatrix} 0 & J \\ J^T & H_L \end{pmatrix}

其中 JJ 是约束的雅可比矩阵,HLH_L 是拉格朗日函数的黑塞矩阵。

约束极值的二阶条件涉及 HLH_L 在约束切空间上的正定性,比无条件极值更复杂。在本科阶段,通常通过实际意义或代入法判断。

几何解释
推论
证明
符号说明

典型例题

例题1:用特征值判断极值

f(x,y)=x3+y33x3yf(x,y) = x^3 + y^3 - 3x - 3y 的极值,并用特征值验证。

参考答案(4 个标签)
黑塞矩阵特征值极值正定性
  1. fx=3x23,fy=3y23f_x = 3x^2-3, f_y = 3y^2-3,驻点 (±1,±1)(\pm1, \pm1),共4个
  2. H=(6x006y)H = \begin{pmatrix} 6x & 0 \\ 0 & 6y \end{pmatrix}
  3. (1,1)(1,1)H=(6006)H=\begin{pmatrix}6&0\\0&6\end{pmatrix},特征值 6,6>06,6>0,正定,极小值 f=4f=-4
  4. (1,1)(-1,-1)H=(6006)H=\begin{pmatrix}-6&0\\0&-6\end{pmatrix},特征值 6,6<0-6,-6<0,负定,极大值 f=4f=4
  5. (1,1)(1,-1)H=(6006)H=\begin{pmatrix}6&0\\0&-6\end{pmatrix},特征值 6,66,-6,不定,鞍点
  6. (1,1)(-1,1)H=(6006)H=\begin{pmatrix}-6&0\\0&6\end{pmatrix},特征值 6,6-6,6,不定,鞍点

答案:极小值 f(1,1)=4f(1,1)=-4,极大值 f(1,1)=4f(-1,-1)=4,其余为鞍点。

例题2:三元函数极值与顺序主子式

f(x,y,z)=x2+y2+z2+2x+4y6zf(x,y,z) = x^2 + y^2 + z^2 + 2x + 4y - 6z 的极值。

参考答案(4 个标签)
黑塞矩阵顺序主子式三元函数极值
  1. fx=2x+2,fy=2y+4,fz=2z6f_x=2x+2, f_y=2y+4, f_z=2z-6,驻点 (1,2,3)(-1,-2,3)
  2. H=(200020002)H = \begin{pmatrix} 2 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 2 \end{pmatrix}
  3. 顺序主子式:D1=2>0,D2=2002=4>0,D3=8>0D_1=2>0, D_2=\begin{vmatrix}2&0\\0&2\end{vmatrix}=4>0, D_3=8>0
  4. 黑塞矩阵正定,故 (1,2,3)(-1,-2,3) 是极小值点
  5. 极小值 f(1,2,3)=1+4+92818=14f(-1,-2,3) = 1+4+9-2-8-18 = -14

答案:极小值为 14-14


总结

本文出现的符号

符号类型读音/说明在本文中的含义
HfH_f黑塞矩阵Hessian of f二阶偏导数构成的矩阵
HBH_B加边黑塞矩阵bordered Hessian约束极值的加边黑塞矩阵
DkD_k顺序主子式leading principal minork阶顺序主子式
λi\lambda_i特征值eigenvalue矩阵的特征值
QQ二次型quadratic form由黑塞矩阵定义的二次型
JJ雅可比矩阵Jacobian约束函数的雅可比矩阵
RnR_n余项remainder泰勒展开的余项

中英对照

中文术语英文术语音标
黑塞矩阵Hessian matrix/ˈhɛsiən ˈmeɪtrɪks/
加边黑塞矩阵bordered Hessian/ˈbɔːrdərd ˈhɛsiən/
正定性positive definiteness/ˈpɒzətɪv ˌdiːˈfɪnətnəs/
负定性negative definiteness/ˈnɛɡətɪv ˌdiːˈfɪnətnəs/
不定性indefiniteness/ɪnˈdɛfɪnətnəs/
半正定positive semidefinite/ˈpɒzətɪv ˌsɛmiˈdɛfənət/
顺序主子式leading principal minor/ˈliːdɪŋ ˈprɪnsəpəl ˈmaɪnər/
特征值eigenvalue/ˈaɪɡənˌvæljuː/
二次型quadratic form/kwɒˈdrætɪk fɔːrm/
泰勒展开Taylor expansion/ˈteɪlər ɪkˈspænʃən/
雅可比矩阵Jacobian matrix/dʒəˈkoʊbiən ˈmeɪtrɪks/
二阶条件second-order condition/ˈsɛkənd ˈɔːrdər kənˈdɪʃən/