MathLabs

分析学

偏导数与梯度

多元函数沿各坐标方向的变化率,汇总起来构成梯度向量。

直观弯曲地形上同时存在的两个坡度

想象站在一个山坡上,高度由二元函数 z=f(x,y)z = f(x,y) 给出,其中 xx 指向正东,yy 指向正北。与一元曲线不同,脚下并没有唯一的斜率:向正东迈步时高度按一种速率变化,向正北迈步时按另一种速率变化,而斜向迈步则会把两者结合起来。固定 yy、只让 xx 变化,相当于用一个东西向的竖直平面切开山坡,把问题化归为熟悉的一元导数。

展示坐标截面与切线斜率的交互式三维曲面图。
交互式三维曲面 z=f(x,y)z = f(x,y):沿固定的 yy 切开会看到斜率为 ∂f∂x\frac{\partial f}{\partial x} 的曲线,沿固定的 xx 切开则得到 ∂f∂y\frac{\partial f}{\partial y}。

大学严格定义:偏导数、切平面与梯度

定义: 偏导数与梯度向量

设 f:U⊆Rn→Rf : U \subseteq \mathbb{R}^n \to \mathbb{R} 定义在开集 UU 上。ff 在点 a∈Ua \in U 处关于 xix_i 的偏导数记作 ∂f∂xi(a)\frac{\partial f}{\partial x_i}(a) 或 fxi(a)f_{x_i}(a),是保持其余所有坐标 xjx_j (j≠ij \ne i) 不变时的极限。将全部 nn 个偏导数合为一个向量,就得到梯度 ∇f(a)=(∂f∂x1(a),…,∂f∂xn(a))\nabla f(a) = \left(\frac{\partial f}{\partial x_1}(a), \dots, \frac{\partial f}{\partial x_n}(a)\right)。

∂f∂xi(a)=lim⁡h→0f(a1,…,ai+h,…,an)−f(a1,…,an)h\frac{\partial f}{\partial x_i}(a) = \lim_{h \to 0} \frac{f(a_1, \dots, a_i + h, \dots, a_n) - f(a_1, \dots, a_n)}{h}

在 (x0,y0)(x_0,y_0) 处仅仅存在偏导数比全微分可微性要弱。若当 (h,k)→(0,0)(h,k) \to (0,0) 时,增量 Δf=f(x0+h,y0+k)−f(x0,y0)\Delta f = f(x_0+h, y_0+k) - f(x_0,y_0) 可以写成 fx(x0,y0)h+fy(x0,y0)k+o(h2+k2)f_x(x_0,y_0)h + f_y(x_0,y_0)k + o\left(\sqrt{h^2+k^2}\right),则称函数 f(x,y)f(x,y) 在 (x0,y0)(x_0,y_0) 处可微。当 ff 可微时,图像 z=f(x,y)z = f(x,y) 在 (x0,y0)(x_0,y_0) 处有良定义的切平面,且沿任意满足 ∥u∥=1\|\mathbf{u}\| = 1 的单位向量 u=(u1,u2)\mathbf{u} = (u_1, u_2) 的方向导数由点积 Duf=∇f⋅uD_{\mathbf{u}}f = \nabla f \cdot \mathbf{u} 给出。

z=f(x0,y0)+∂f∂x(x0,y0)(x−x0)+∂f∂y(x0,y0)(y−y0),Duf(x0,y0)=∇f(x0,y0)⋅uz = f(x_0,y_0) + \frac{\partial f}{\partial x}(x_0,y_0)(x - x_0) + \frac{\partial f}{\partial y}(x_0,y_0)(y - y_0), \qquad D_{\mathbf{u}}f(x_0,y_0) = \nabla f(x_0,y_0) \cdot \mathbf{u}
多元微分算子与二阶量汇总
对象公式几何 / 物理意义
偏导数∂f∂xi\frac{\partial f}{\partial x_i}沿坐标轴 xix_i 截面的斜率。
梯度向量∇f=(fx1,…,fxn)\nabla f = (f_{x_1}, \dots, f_{x_n})垂直于等值集并指向最陡上升方向;范数 ∥∇f∥\|\nabla f\| 为最大变化率。
方向导数Duf=∇f⋅uD_{\mathbf{u}}f = \nabla f \cdot \mathbf{u}以单位速度向量 u\mathbf{u} 移动时的瞬时变化率。
海森矩阵与判别式D=fxxfyy−fxy2D = f_{xx}f_{yy} - f_{xy}^2度量局部曲率,将驻点分类为极值点或鞍点。

大学核心定理:最陡上升、施瓦茨对称性与海森判别法

若 ff 在点 aa 可微且 ∇f(a)≠0\nabla f(a) \ne \mathbf{0},则对任意单位向量 u\mathbf{u}(满足 ∥u∥=1\|\mathbf{u}\| = 1),方向导数满足 Duf(a)=∇f(a)⋅u=∥∇f(a)∥cos⁡θD_{\mathbf{u}}f(a) = \nabla f(a) \cdot \mathbf{u} = \|\nabla f(a)\| \cos\theta,其中 θ\theta 是 ∇f(a)\nabla f(a) 与 u\mathbf{u} 之间的夹角。因此,当 u=∇f(a)∥∇f(a)∥\mathbf{u} = \frac{\nabla f(a)}{\|\nabla f(a)\|} 时,Duf(a)D_{\mathbf{u}}f(a) 取得最大值 ∥∇f(a)∥\|\nabla f(a)\|,在相反方向上取得最小值 −∥∇f(a)∥-\|\nabla f(a)\|。

为什么成立?

点积将梯度投影到所选的行进方向上:只有完全与梯度同向时才能获得梯度的全部模长,而垂直于梯度沿等高线行走时变化率为零。

证明

定义一元函数 g(t)=f(a+tu)g(t) = f(a + t\mathbf{u})。由方向导数的定义知 Duf(a)=g′(0)D_{\mathbf{u}}f(a) = g'(0)。因为 ff 在点 aa 可微,当 t→0t \to 0 时有 f(a+tu)−f(a)=∇f(a)⋅(tu)+o(∣t∣)f(a + t\mathbf{u}) - f(a) = \nabla f(a) \cdot (t\mathbf{u}) + o(|t|)。两边除以 tt 并取 t→0t \to 0 的极限,即得 Duf(a)=∇f(a)⋅uD_{\mathbf{u}}f(a) = \nabla f(a) \cdot \mathbf{u}。

由欧几里得内积的几何公式以及单位长度条件 ∥u∥=1\|\mathbf{u}\| = 1,可得 ∇f(a)⋅u=∥∇f(a)∥ ∥u∥cos⁡θ=∥∇f(a)∥cos⁡θ\nabla f(a) \cdot \mathbf{u} = \|\nabla f(a)\|\,\|\mathbf{u}\|\cos\theta = \|\nabla f(a)\|\cos\theta。由于 −1≤cos⁡θ≤1-1 \le \cos\theta \le 1,最大值 ∥∇f(a)∥\|\nabla f(a)\| 唯一地在 θ=0\theta = 0 即 u=∇f(a)∥∇f(a)∥\mathbf{u} = \frac{\nabla f(a)}{\|\nabla f(a)\|} 处取得,最小值 −∥∇f(a)∥-\|\nabla f(a)\| 在 θ=π\theta = \pi 处取得。

若 f(x,y)f(x,y) 在 (x0,y0)(x_0,y_0) 的某邻域内具有连续二阶偏导数,则 ∂2f∂x∂y(x0,y0)=∂2f∂y∂x(x0,y0)\frac{\partial^2 f}{\partial x\partial y}(x_0,y_0) = \frac{\partial^2 f}{\partial y\partial x}(x_0,y_0)。此外,若 ∇f(x0,y0)=0\nabla f(x_0,y_0) = \mathbf{0} 且记 D=fxx(x0,y0)fyy(x0,y0)−fxy(x0,y0)2D = f_{xx}(x_0,y_0)f_{yy}(x_0,y_0) - f_{xy}(x_0,y_0)^2,则当 D>0D > 0 且 fxx>0f_{xx} > 0 时 (x0,y0)(x_0,y_0) 为严格局部极小值点;当 D>0D > 0 且 fxx<0f_{xx} < 0 时为严格局部极大值点;当 D<0D < 0 时为鞍点。

为什么成立?

混合偏导数的对称性保证了海森矩阵是对称矩阵,对其二次型配方即可看出曲面在水平切平面附近是向上凹陷如碗、向下拱起如穹顶,还是像马鞍一样扭曲。

证明

对非零小量 h,kh, k,考虑二阶差分 Δ(h,k)=f(x0+h,y0+k)−f(x0+h,y0)−f(x0,y0+k)+f(x0,y0)\Delta(h,k) = f(x_0+h,y_0+k) - f(x_0+h,y_0) - f(x_0,y_0+k) + f(x_0,y_0)。先对 [x0,x0+h][x_0, x_0+h] 上的 g(x)=f(x,y0+k)−f(x,y0)g(x) = f(x,y_0+k) - f(x,y_0) 应用一元中值定理,再沿 yy 对 fxf_x 应用中值定理,得 Δ(h,k)=hk fyx(c1,d1)\Delta(h,k) = hk\,f_{yx}(c_1, d_1);按相反顺序应用则得 Δ(h,k)=hk fxy(c2,d2)\Delta(h,k) = hk\,f_{xy}(c_2, d_2),其中各中间点均收敛于 (x0,y0)(x_0,y_0)。令两者相等并利用连续性取极限 (h,k)→(0,0)(h,k) \to (0,0),即证得 ∂2f∂x∂y(x0,y0)=∂2f∂y∂x(x0,y0)\frac{\partial^2 f}{\partial x\partial y}(x_0,y_0) = \frac{\partial^2 f}{\partial y\partial x}(x_0,y_0)。

在满足 ∇f(x0,y0)=0\nabla f(x_0,y_0) = \mathbf{0} 的驻点处,泰勒公式给出 f(x0+h,y0+k)−f(x0,y0)=12Q(h,k)+o(h2+k2)f(x_0+h,y_0+k) - f(x_0,y_0) = \frac{1}{2}Q(h,k) + o(h^2+k^2),其中 Q(h,k)=fxxh2+2fxyhk+fyyk2Q(h,k) = f_{xx}h^2 + 2f_{xy}hk + f_{yy}k^2。当 fxx≠0f_{xx} \ne 0 时,配方可得 Q(h,k)=1fxx[(fxxh+fxyk)2+Dk2]Q(h,k) = \frac{1}{f_{xx}}\left[(f_{xx}h + f_{xy}k)^2 + Dk^2\right],其中 D=fxxfyy−fxy2D = f_{xx}f_{yy} - f_{xy}^2。若 D>0D > 0,方括号内的和在任意 (h,k)≠(0,0)(h,k) \ne (0,0) 处严格为正,故 Q(h,k)Q(h,k) 的符号与 fxxf_{xx} 一致(fxx>0f_{xx} > 0 时为严格极小,fxx<0f_{xx} < 0 时为严格极大)。若 D<0D < 0,则 Q(h,k)Q(h,k) 沿过原点的不同直线可正可负,从而形成鞍点。

大学实际应用与典型例题

偏导数与梯度驱动着现代科学与工程:在物理学中,傅里叶热传导定律指出热流密度正比于负温度梯度 −∇T-\nabla T;在机器学习中,训练神经网络通过不断沿负梯度 −∇L(θ)-\nabla L(\theta) 迈步来最小化损失函数 L(θ)L(\theta)(梯度下降法);而在经济学与力学中,海森判别法用于区分稳定平衡态(势能的局部极小值)与不稳定的鞍点。

例题: 温度场的方向导数与最陡上升方向

一块金属板在坐标 (x,y)(x,y)(单位:米)处的温度为 T(x,y)=x2+3xy+y3T(x,y) = x^2 + 3xy + y^3(单位:摄氏度)。求点 (1,2)(1,2) 处的梯度 ∇T(1,2)\nabla T(1,2)、沿向量 v=(3,4)\mathbf{v} = (3,4) 方向的方向导数,以及点 (1,2)(1,2) 处温度的最大增长率。

解答

首先,固定一个变量对另一个变量求导,计算偏导数:Tx(x,y)=2x+3yT_x(x,y) = 2x + 3y 和 Ty(x,y)=3x+3y2T_y(x,y) = 3x + 3y^2。代入点 (1,2)(1,2) 得 Tx(1,2)=2(1)+3(2)=8T_x(1,2) = 2(1) + 3(2) = 8 及 Ty(1,2)=3(1)+3(2)2=15T_y(1,2) = 3(1) + 3(2)^2 = 15,因此梯度向量为 ∇T(1,2)=(8,15)\nabla T(1,2) = (8, 15)。

接着,由于 ∥v∥=32+42=5\|\mathbf{v}\| = \sqrt{3^2+4^2} = 5,将 v=(3,4)\mathbf{v} = (3,4) 单位化得到单位方向向量 u=v∥v∥=(35,45)\mathbf{u} = \frac{\mathbf{v}}{\|\mathbf{v}\|} = \left(\frac{3}{5}, \frac{4}{5}\right)。方向导数为 DuT(1,2)=∇T(1,2)⋅u=8⋅35+15⋅45=845=16.8D_{\mathbf{u}}T(1,2) = \nabla T(1,2) \cdot \mathbf{u} = 8 \cdot \frac{3}{5} + 15 \cdot \frac{4}{5} = \frac{84}{5} = 16.8。点 (1,2)(1,2) 处的最大增长率沿 ∇T(1,2)\nabla T(1,2) 方向取得,等于 ∥∇T(1,2)∥=82+152=289=17\|\nabla T(1,2)\| = \sqrt{8^2 + 15^2} = \sqrt{289} = 17 摄氏度每米。

例题: 用海森行列式分类驻点

利用二阶导数海森判别法求出并分类函数 f(x,y)=x3+y3−3xyf(x,y) = x^3 + y^3 - 3xy 的所有驻点。

解答

令两个一阶偏导数同时为零:fx(x,y)=3x2−3y=0f_x(x,y) = 3x^2 - 3y = 0 且 fy(x,y)=3y2−3x=0f_y(x,y) = 3y^2 - 3x = 0。由第一式得 y=x2y = x^2,代入第二式得 3(x2)2−3x=3x(x3−1)=03(x^2)^2 - 3x = 3x(x^3 - 1) = 0,实数解为 x=0x = 0 和 x=1x = 1。因此两个驻点为 (0,0)(0,0) 和 (1,1)(1,1)。

接着计算二阶偏导数:fxx=6xf_{xx} = 6x,fyy=6yf_{yy} = 6y,fxy=−3f_{xy} = -3,从而海森判别式为 D(x,y)=fxxfyy−fxy2=36xy−9D(x,y) = f_{xx}f_{yy} - f_{xy}^2 = 36xy - 9。在 (0,0)(0,0) 处,D(0,0)=−9<0D(0,0) = -9 < 0,故 (0,0)(0,0) 是鞍点。在 (1,1)(1,1) 处,D(1,1)=36−9=27>0D(1,1) = 36 - 9 = 27 > 0 且 fxx(1,1)=6>0f_{xx}(1,1) = 6 > 0,故 (1,1)(1,1) 是严格局部极小值点,极小值为 f(1,1)=−1f(1,1) = -1。

设函数 f(x,y)=x3y2−4x+5yf(x,y) = x^3 y^2 - 4x + 5y,求梯度向量 ∇f(1,2)\nabla f(1, 2)。

若 ∇f(x0,y0)=(6,−8)\nabla f(x_0, y_0) = (6, -8),则在所有单位向量 u\mathbf{u} 上方向导数 Duf(x0,y0)D_{\mathbf{u}}f(x_0, y_0) 的可能最大值是多少?

在机器学习中,梯度下降法通过更新参数 θ\theta 来最小化可微损失函数 L(θ)L(\theta)。下列哪个更新公式使 θ\theta 沿 LL 局部下降最快的方向移动?

设 (x0,y0)(x_0, y_0) 是光滑函数 f(x,y)f(x,y) 的驻点,且 fxx(x0,y0)=−4f_{xx}(x_0,y_0) = -4,fyy(x0,y0)=−3f_{yy}(x_0,y_0) = -3,fxy(x0,y0)=2f_{xy}(x_0,y_0) = 2。点 (x0,y0)(x_0, y_0) 应如何分类?

参考文献

  1. Jerrold E. Marsden, Anthony J. Tromba (2012). Vector Calculus
  2. James Stewart (2015). Calculus: Early Transcendentals
  3. Augustin-Louis Cauchy (1847). Méthode générale pour la résolution des systèmes d'équations simultanées