多元微积分 · 二阶偏导数判别法

二阶偏导数判别法

在临界点,二阶偏导数组成的D = f_xx f_yy − f_xy²判断它是极大、极小还是鞍点;D = 0时判别法不下结论。转动过这一点的切片,看它怎样弯。

二阶偏导数判别法与本页画面

在函数f(x,y)f(x, y)的临界点P,两个一阶偏导数都等于0,切平面是水平的。二阶偏导数判别法用二阶导数A=fxxA = f_{xx}、B=fxyB = f_{xy}、C=fyyC = f_{yy}组成判别式D=AC−B2D = AC - B^2,也就是黑塞矩阵的行列式,由它判断这是什么样的点。若f在P附近有连续的二阶偏导数,则:

临界点处的判别法
{极小值,D>0, fxx>0极大值,D>0, fxx<0鞍点,D<0判别法无法判断,D=0\begin{cases} \text{极小值}, & D > 0,\ f_{xx} > 0 \\ \text{极大值}, & D > 0,\ f_{xx} < 0 \\ \text{鞍点}, & D < 0 \\ \text{判别法无法判断}, & D = 0 \end{cases}

页面用过临界点P的竖直切片,沿任意方向θ切开一个曲面。这个曲面可以是二阶导数由你设定的二次型,也可以是另外几个函数之一。切口在P处弯得多厉害,就是沿这个方向的二阶导数fuuf_{uu}。P周围的罗盘上,切片向上弯的方向涂满,向下弯的方向画斜线。判别法说的,其实是所有这些切片合在一起的情形。

各个方向上切片弯得多厉害
fuu(θ)=Acos⁡2θ+2Bcos⁡θsin⁡θ+Csin⁡2θ=m+scos⁡2(θ−θ0)f_{uu}(\theta) = A\cos^2\theta + 2B\cos\theta\sin\theta + C\sin^2\theta = m + s\cos 2(\theta - \theta_0)

θ转过半圈,弯曲程度绕着中线m=12(A+C)m = \tfrac12(A + C)走完一整个波,振幅是s,画在小图里。波的最高值λ1=m+s\lambda_1 = m + s和最低值λ2=m−s\lambda_2 = m - s出现在两个互相垂直的方向上,它们就是黑塞矩阵的特征值。乘开来得D=m2−s2=λ1λ2D = m^2 - s^2 = \lambda_1\lambda_2。所以D>0D > 0意味着波碰不到0。每个切片都朝同一边弯,P是极小值点或极大值点。这时随便看一个切片就知道是哪一种,fxxf_{xx}不过是沿x方向的弯曲。D<0D < 0意味着波穿过0:有的切片向上弯,有的向下弯,P是鞍点。

所以光有fxx>0f_{xx} > 0和fyy>0f_{yy} > 0还不够。交叉项fxyf_{xy}把波峰从坐标轴转开,让波变高,中线却不动。∣fxy∣|f_{xy}|大到一定程度,沿对角线的切片就会向下弯,尽管两条轴上的切片都向上弯,碗也就变成了鞍。

D<0D < 0时,恰好有两个方向的弯曲为0,就是波穿过0的地方,罗盘把它们画成横过圆盘的直线。它们是过P的那条等高线的切线,等高线在P处和自己相交。这就是鞍点在等高线图上画出的X形。

D=0D = 0时,波碰到0或者整条落在0上,所以在某个方向上二阶项为0,判别法说不准。x2−y4x^2 - y^4和x2+y4x^2 + y^4在原点的二阶导数完全相同,但沿y方向,前者像−t4-t^4那样下降,后者像t4t^4那样上升,一个是鞍点,一个是极小值点。猴鞍x3−3xy2x^3 - 3xy^2在原点的二阶导数全是0,曲面朝三个方向上升,朝另外三个方向下降。“判别法无法判断”不等于“没有极值”。

在临界点,一阶项都是0,所以在P附近,函数很接近它的二阶泰勒多项式,也就是f(P)加上一个二次型。“二次近似”把它画在圆盘上。D>0D > 0时,离P足够近的地方,二次项在每个方向上都压得过余项,所以P和二次近似一样,是极小值点或极大值点。D<0D < 0时,在二次项向上弯和向下弯的那些方向上,它仍然压得过余项,于是f沿一些方向上升、沿另一些方向下降:P是鞍点。

临界点处的二阶泰勒多项式
f(P+h)≈f(P)+12(Ahx2+2Bhxhy+Chy2)f(P + h) \approx f(P) + \tfrac12\big(A h_x^2 + 2B h_x h_y + C h_y^2\big)

值得注意

  • 临界点是局部的概念。x3−3x+y2x^3 - 3x + y^2在(1, 0)处的极小值并不是最低点:x越来越小时,函数会一直往下降。

操作

  • 转动切片有三种办法:拖罗盘上的红把手,转方向拨盘,或者在小图里拖动波。
  • θ栏可以填角度,也可以填π的倍数,比如π/4。

相关概念