5.3 多元数量值函数的导数与微分

本节将把一元函数的导数与微分概念推广到多元数量值函数。我们将以二元函数为主进行讲解，然后推广到 $n$ 元函数。先介绍多元数量值函数的偏导数与全微分以及方向导数与梯度，再介绍高阶偏导数与高阶全微分以及复合函数的链式法则，最后介绍隐函数及其微分法。

3.1 偏导数

我们知道，一元函数在一点的导数表示函数在该点的变化率，它反映了在该点处函数值随自变量变化的快慢程度。对于二元函数 $z=f(x,y)$ 来说，当然也需要研究它的变化率问题。但是由于自变量多了一个，点 $(x,y)$ 在 $xOy$ 平面上变化，情况要复杂得多，因而因变量随自变量变化的情况也要比一元函数复杂得多。我们首先来研究二元函数关于它的一个自变量的变化率。例如，对于二元函数 $z=f(x,y)$，在点 $(x_0,y_0)$ 处，把 $y$ 固定在 $y=y_0$，令 $x$ 变化，这时候它就是 $x$ 的一元函数，该函数在 $x_0$ 处的导数，就称为函数 $f(x,y)$ 在点 $(x_0,y_0)$ 处对 $x$ 的偏导数。

定义 3.1 (偏导数)

设函数 $z=f(x,y)$ 在点 $(x_0,y_0)$ 的某一邻域 $U(x_0,y_0)$ 内有定义，当自变量 $y$ 固定在 $y=y_0$，而 $x$ 在 $x_0$ 处有改变量 $\Delta x$，$(x_0+\Delta x,y_0)\in U(x_0,y_0)$ 时，相应地，函数 $f$ 有改变量

f(x₀ + Δ x, y₀) - f(x₀, y₀).$$f(x_0+\Delta x,y_0)-f(x_0,y_0).$$

如果极限

lim_(Δ x → 0) (f(x₀ + Δ x, y₀) - f(x₀, y₀))/(Δ x)$$\underset{\Delta x\to 0}{\lim }\frac{f(x_0+\Delta x,y_0)-f(x_0,y_0)}{\Delta x}$$

存在，则称此极限值为函数 $z=f(x,y)$ 在点 $(x_0,y_0)$ 对 $x$ 的偏导数，记作

fₓ(x₀, y₀), (∂ f(x₀, y₀))/(∂ x), zₓ(x₀, y₀) 或 . (∂ z)/(∂ x) |_((x₀, y₀)),$$f_x(x_0,y_0),\frac{\partial f(x_0,y_0)}{\partial x},z_x(x_0,y_0)\text{或}{\frac{\partial z}{\partial x}|}_{(x_0,y_0)},$$

即

fₓ(x₀, y₀) = (∂ f(x₀, y₀))/(∂ x) = lim_(Δ x → 0) (f(x₀ + Δ x, y₀) - f(x₀, y₀))/(Δ x). (3.1)$$f_x(x_0,y_0)=\frac{\partial f(x_0,y_0)}{\partial x}=\underset{\Delta x\to 0}{\lim }\frac{f(x_0+\Delta x,y_0)-f(x_0,y_0)}{\Delta x}.\text{\hspace{1em}(3.1)}$$

想一想

AI 解答

试说明左边这段话的正确性。提示：设 $z=f(x,y)$ 在点 $(x_0,y_0)$ 处可偏导，令 $\phi (x)=f(x,y_0)$ (或 $\psi (y)=f(x_0,y)$)，证明 ${\phi }^{\prime }(x_0)=f_x(x_0,y_0)$ (或 ${\psi }^{\prime }(y_0)=f_y(x_0,y_0)$)。

类似地，可定义函数 $z=f(x,y)$ 在 $(x_0,y_0)$ 对 $y$ 的偏导数，记作

f_y(x₀, y₀), (∂ f(x₀, y₀))/(∂ y), z_y(x₀, y₀) 或 . (∂ z)/(∂ y) |_((x₀, y₀)),$$f_y(x_0,y_0),\frac{\partial f(x_0,y_0)}{\partial y},z_y(x_0,y_0)\text{或}{\frac{\partial z}{\partial y}|}_{(x_0,y_0)},$$

即

f_y(x₀, y₀) = (∂ f(x₀, y₀))/(∂ y) = lim_(Δ y → 0) (f(x₀, y₀ + Δ y) - f(x₀, y₀))/(Δ y). (3.2)$$f_y(x_0,y_0)=\frac{\partial f(x_0,y_0)}{\partial y}=\underset{\Delta y\to 0}{\lim }\frac{f(x_0,y_0+\Delta y)-f(x_0,y_0)}{\Delta y}.\text{\hspace{1em}(3.2)}$$

如果二元函数 $z=f(x,y)$ 在点 $(x_0,y_0)$ 对 $x$ 与对 $y$ 的偏导数均存在，那么称 $f(x,y)$ 在 $(x_0,y_0)$ 处可偏导。

由 (3.1) 式可见，$f(x,y)$ 在 $(x_0,y_0)$ 对 $x$ 的偏导数，实际上就是把 $y$ 固定在 $y_0$ 时，关于 $x$ 的一元函数 $f(x,y_0)$ 在 $x_0$ 处的导数；同理，由 (3.2) 式可见，偏导数 $\frac{\partial f}{\partial y}$ —— 就是把 $x$ 固定在 $x_0$ 时，关于 $y$ 的一元函数 $f(x_0,y)$ 在 $y_0$ 处的导数。

与一元函数类似，若 $f$ 在区域 $D\subseteq {\mathbf{R}}^2$ 内有定义，我们可以定义 $f$ 在区域 $D$ 内对 $x$ 及对 $y$ 的偏导函数分别为

(∂ f(x,y))/(∂ x) = lim_(Δ x → 0) (f(x + Δ x, y) - f(x, y))/(Δ x),$$\frac{\partial f(x,y)}{\partial x}=\underset{\Delta x\to 0}{\lim }\frac{f(x+\Delta x,y)-f(x,y)}{\Delta x},$$ (∂ f(x,y))/(∂ y) = lim_(Δ y → 0) (f(x, y + Δ y) - f(x, y))/(Δ y),$$\frac{\partial f(x,y)}{\partial y}=\underset{\Delta y\to 0}{\lim }\frac{f(x,y+\Delta y)-f(x,y)}{\Delta y},$$

其中 $(x,y)\in D,(x+\Delta x,y)\in D,(x,y+\Delta y)\in D$。函数 $z=f(x,y)$ 对 $x$ 的偏导函数可简记为

fₓ, (∂ f)/(∂ x), zₓ, 或 (∂ z)/(∂ x),$$f_x,\frac{\partial f}{\partial x},z_x,\text{或}\frac{\partial z}{\partial x},$$

$f$ 对 $y$ 的偏导函数可简记为

f_y, (∂ f)/(∂ y), z_y, 或 (∂ z)/(∂ y).$$f_y,\frac{\partial f}{\partial y},z_y,\text{或}\frac{\partial z}{\partial y}.$$

注意

由偏导数的定义可见，求函数 $z=f(x,y)$ 的偏导数，并不需要新的方法。因为求 $\frac{\partial z}{\partial x}$ 时，是把 $y$ 暂时看作常量而对 $x$ 求导数；求 $\frac{\partial z}{\partial y}$ 时，是把 $x$ 暂时看作常量而对 $y$ 求导数。所以，求偏导数的方法实际上仍旧是一元函数的求导方法。

在不致混淆时，偏导函数也简称为偏导数。由偏导函数的定义易知，$f(x,y)$ 在 $(x_0,y_0)$ 对 $x$ 的偏导数 $f_x(x_0,y_0)$ 就是偏导函数 $f_x(x,y)$ 在 $(x_0,y_0)$ 处的函数值，而 $f_y(x_0,y_0)$ 就是偏导函数 $f_y(x,y)$ 在 $(x_0,y_0)$ 处的函数值。

例 3.1 偏导数的计算与求值

设二元函数 $z=\arctan \frac{y}{x}$，求 $\frac{\partial z}{\partial x},\frac{\partial z}{\partial y}$ 及 ${\frac{\partial z}{\partial x}|}_{(1,1)}$。

解

把 $y$ 看作常数，对 $x$ 求导得

(∂ z)/(∂ x) = (1)/(1 + (y/x)²) · (-y)/(x²) = -(y)/(x² + y²).$$\frac{\partial z}{\partial x}=\frac{1}{1+{\left(\frac{y}{x}\right)}^2}\cdot \frac{-y}{x^2}=-\frac{y}{x^2+y^2}.$$

把 $x$ 看作常数，对 $y$ 求导得

(∂ z)/(∂ y) = (1)/(1 + (y/x)²) · 1/x = (x)/(x² + y²).$$\frac{\partial z}{\partial y}=\frac{1}{1+{\left(\frac{y}{x}\right)}^2}\cdot \frac{1}{x}=\frac{x}{x^2+y^2}.$$

由上述表达式得

. (∂ z)/(∂ x) |_((1,1)) = -(y)/(x² + y²) g|ₓ=1, y=1 = -1/2.$${\frac{\partial z}{\partial x}|}_{(1,1)}=-\frac{y}{x^2+y^2}{|}_{x=1,y=1}=-\frac{1}{2}.$$

或

. (∂ z)/(∂ x) |_((1,1)) = . (d)/(dx) ( arctan 1/x ) |ₓ=1 = . (1)/(1 + (1)/(x²)) ( -(1)/(x²) ) |ₓ=1 = -1/2.$${\frac{\partial z}{\partial x}|}_{(1,1)}={\frac{\mathrm{d}}{\mathrm{d}x}\left(\arctan \frac{1}{x}\right)|}_{x=1}={\frac{1}{1+\frac{1}{x^2}}\left(-\frac{1}{x^2}\right)|}_{x=1}=-\frac{1}{2}.$$

例 3.2 试由理想气体的状态方程

试由理想气体的状态方程

p = R T/V (R 为常数) (3.3)$$p=R\frac{T}{V}(R\text{ 为常数})\text{\hspace{1em}(3.3)}$$

求 $\frac{\partial p}{\partial T},\frac{\partial T}{\partial V},\frac{\partial V}{\partial p}$。

解

将 $V$ 看作常量（即等容变化过程），对 $T$ 求导得 $\frac{\partial p}{\partial T}=\frac{R}{V}$；由 (3.3) 式解得 $T=\frac{1}{R}pV$，将 $p$ 看作常量（即等压变化过程），对 $V$ 求导得 $\frac{\partial T}{\partial V}=\frac{p}{R}$；类似可求得 $\frac{\partial V}{\partial p}=-\frac{RT}{p^2}$。所以

(∂ p)/(∂ T) · (∂ T)/(∂ V) · (∂ V)/(∂ p) = R/V · p/R · ( -(RT)/(p²) ) = -(RT)/(Vp) = -1.$$\frac{\partial p}{\partial T}\cdot \frac{\partial T}{\partial V}\cdot \frac{\partial V}{\partial p}=\frac{R}{V}\cdot \frac{p}{R}\cdot \left(-\frac{RT}{p^2}\right)=-\frac{RT}{Vp}=-1.$$

注意

从例3.2 的结果可以看出，偏导数记号 $\frac{\partial z}{\partial x}$ 是一个整体记号，不能将它拆开，像一元函数那样将导数 $\frac{\mathrm{d}y}{\mathrm{d}x}$ 看作微商（微分之商）。

例 3.3 研究二元函数 $f(x,y)$ 在 $(0,0)$ 处是否可偏导

研究二元函数

f(x,y) = (xy)/(x² + y²), & x² + y² ≠ 0, 0, & x² + y² = 0$$f(x,y)=\{\begin{array}{ll}\frac{xy}{x^2+y^2},& x^2+y^2\ne 0,\\ 0,& x^2+y^2=0\end{array}$$

在点 $(0,0)$ 处是否可偏导。

解

由于

(f(0 + Δ x, 0) - f(0,0))/(Δ x) = 0, (f(0, 0 + Δ y) - f(0,0))/(Δ y) = 0,$$\frac{f(0+\Delta x,0)-f(0,0)}{\Delta x}=0,\frac{f(0,0+\Delta y)-f(0,0)}{\Delta y}=0,$$

故由偏导数的定义可得 $f(x,y)$ 在 $(0,0)$ 处可偏导（两个偏导数均存在），且 $f_x(0,0)=0,f_y(0,0)=0$。

二元函数偏导数的几何意义

根据定义3.1，二元函数 $z=f(x,y)$ 在点 $(x_0,y_0)$ 对 $x$ 的偏导数 $f_x(x_0,y_0)$ 就是一元函数 $z=f(x,y_0)$ 在 $x_0$ 处的导数，而二元函数 $z=f(x,y)$ 在几何上表示空间的一个曲面，函数 $z=f(x,y_0)$ 在几何上就表示曲面 $z=f(x,y)$ 与平面 $y=y_0$ 的交线，

Cₓ: z = f(x,y), y = y₀.$$C_x:\{\begin{array}{l}z=f(x,y),\\ y=y_0.\end{array}$$

于是由一元函数导数的几何意义立即可知 $f_x(x_0,y_0)$ 就是曲线 $C_x$ 在点 $M(x_0,y_0,f(x_0,y_0))$ 处切线 $T_x$ 的斜率（图 5.10），即

fₓ(x₀, y₀) = tan α,$$f_x(x_0,y_0)=\tan \alpha ,$$

同理可知，$f_y(x_0,y_0)$ 就是曲线 $C_y:\{\begin{array}{l}z=f(x,y),\\ x=x_0.\end{array}$ 在点 $M$ 处切线 $T_y$ 的斜率，即

f_y(x₀, y₀) = tan β.$$f_y(x_0,y_0)=\tan \beta .$$

图 5.10 偏导数的几何意义

应当指出，与一元函数不同，对于二元函数 $f$ 来说，它在一点处可偏导，并不能保证 $f$ 在此点连续。例如，例3.3 中的函数 $f$ 在点 $(0,0)$ 处可偏导，但在第二节中我们已经看到该函数在点 $(0,0)$ 却不连续。因此，一元函数在其可导点处一定连续的结论不能推广到二元函数。实际上，当函数 $f(x,y)$ 在 $(x_0,y_0)$ 处的两个偏导数均存在时，根据一元函数可导必连续的结论，仅能保证当 $(x,y)$ 分别沿平行于 $x$ 轴和 $y$ 轴两个特殊路径趋近于点 $(x_0,y_0)$ 时，$f(x,y)\to f(x_0,y_0)$，这当然不能保证当 $(x,y)$ 以任意路径和方式趋近于点 $(x_0,y_0)$ 时都有 $f(x,y)\to f(x_0,y_0)$，这就是不能保证 $f(x,y)$ 在点 $(x_0,y_0)$ 连续的原因。

仿照定义3.1，可以给出 $n$ 元函数的偏导数的定义。例如，设三元函数 $u=f(x,y,z)$ 定义在 $(x_0,y_0,z_0)$ 的某一邻域 $U(x_0,y_0,z_0)$ 内，$u$ 在 $(x_0,y_0,z_0)$ 对 $x$ 的偏导数的定义为

. (∂ u)/(∂ x) |_((x₀, y₀, z₀)) = lim_(Δ x → 0) (f(x₀ + Δ x, y₀, z₀) - f(x₀, y₀, z₀))/(Δ x).$${\frac{\partial u}{\partial x}|}_{(x_0,y_0,z_0)}=\underset{\Delta x\to 0}{\lim }\frac{f(x_0+\Delta x,y_0,z_0)-f(x_0,y_0,z_0)}{\Delta x}.$$

可以看出，与二元函数类似，它实际上是将 $y$ 固定在 $y_0$、将 $z$ 固定在 $z_0$ 时，一元函数 $f(x,y_0,z_0)$ 在 $x_0$ 处的导数，即

. (∂ u)/(∂ x) |_((x₀, y₀, z₀)) = . (d)/(dx) f(x, y₀, z₀) |ₓ=x₀.$${\frac{\partial u}{\partial x}|}_{(x_0,y_0,z_0)}={\frac{\mathrm{d}}{\mathrm{d}x}f(x,y_0,z_0)|}_{x=x_0}.$$

一般地，$n$ 元函数

u = f(x₁, x₂, …, xₙ)$$u=f(x_1,x_2,\dots ,x_n)$$

在 $x_0(x_{0,1},x_{0,2},\dots ,x_{0,n})$ 处对变量 $x_i$ ($i=1,2,\dots ,n$) 的偏导数也定义为将其余变量 $x_j$ ($j=1,\dots ,i-1,i+1,\dots ,n$) 固定时，关于 $x_i$ 的一元函数 $f(x_{0,1},\dots ,x_{0,i-1},x_i,x_{0,i+1},\dots ,x_{0,n})$ 在 $x_{0,i}$ 处的导数，即

(∂ f(x₀))/(∂ xᵢ) = . (d)/(dxᵢ) f(x₀,1, …, x₀,i-1, xᵢ, x₀,i+1, …, x₀,n) |ₓᵢ=x₀,i. (3.4)$$\frac{\partial f(x_0)}{\partial x_i}={\frac{\mathrm{d}}{\mathrm{d}{x}_i}f(x_{0,1},\dots ,x_{0,i-1},x_i,x_{0,i+1},\dots ,x_{0,n})|}_{x_i=x_{0,i}}.\text{\hspace{1em}(3.4)}$$

因此，同二元函数一样，我们在求 $n$ 元函数 $f(x_1,x_2,\dots ,x_n)$ 对某一自变量 $x_i$ 的偏导函数时，只要将 $f$ 看作 $x_i$ 的一元函数（将其他自变量都暂时看作常量），利用一元函数的求导法则，便可求得所要求的偏导数 $\frac{\partial f}{\partial x_i}$。

例 3.4 三元函数偏导数的计算与求值

设三元函数 $u(x,y,z)=\sqrt[z]{\frac{y}{x}}$，求 $u$ 对各个自变量的偏导数及 ${\frac{\partial u}{\partial x}|}_{(1,2,1)}$。

解

记 $z=f(x,y)=e^{\frac{1}{z}\ln (\frac{y}{x})}$，则

(∂ u)/(∂ x) = 1/z ( y/x )¹/z-1 ( -(y)/(x²) ), (∂ u)/(∂ y) = 1/z ( y/x )¹/z-1 1/x, (∂ u)/(∂ z) = ( y/x )¹/z ln( y/x ) ( -(1)/(z²) ),$$\frac{\partial u}{\partial x}=\frac{1}{z}{\left(\frac{y}{x}\right)}^{\frac{1}{z}-1}\left(-\frac{y}{x^2}\right),\frac{\partial u}{\partial y}=\frac{1}{z}{\left(\frac{y}{x}\right)}^{\frac{1}{z}-1}\frac{1}{x},\frac{\partial u}{\partial z}={\left(\frac{y}{x}\right)}^{\frac{1}{z}}\ln \left(\frac{y}{x}\right)\left(-\frac{1}{z^2}\right),$$. (∂ u)/(∂ x) |_((1,2,1)) = 1/1 ( 2/1 )¹/1-1 ( -(2)/(1²) ) g|_((1,2,1)) = -2.$${\frac{\partial u}{\partial x}|}_{(1,2,1)}=\frac{1}{1}{\left(\frac{2}{1}\right)}^{\frac{1}{1}-1}\left(-\frac{2}{1^2}\right){|}_{(1,2,1)}=-2.$$

或

. (∂ u)/(∂ x) |_((1,2,1)) = . (d)/(dx) u(x, 2, 1) |ₓ=1 = . (d)/(dx) ( 2/x ) |ₓ=1 = . -(2)/(x²) |ₓ=1 = -2.$${\frac{\partial u}{\partial x}|}_{(1,2,1)}={\frac{\mathrm{d}}{\mathrm{d}x}u(x,2,1)|}_{x=1}={\frac{\mathrm{d}}{\mathrm{d}x}\left(\frac{2}{x}\right)|}_{x=1}={-\frac{2}{x^2}|}_{x=1}=-2.$$

3.2 全微分

对于一元函数 $f:U(x_0)\subseteq \mathbf{R}\to \mathbf{R}$，若存在一个关于 $\Delta x$ 的线性函数 $L(\Delta x)=\alpha \Delta x$，使函数 $f$ 的改变量可表示为

f(x₀ + Δ x) - f(x₀) = α Δ x + o(Δ x),$$f(x_0+\Delta x)-f(x_0)=\alpha \Delta x+o(\Delta x),$$

其中常数 $\alpha$ 与 $\Delta x$ 无关，$o(\Delta x)$ 是当 $\Delta x\to 0$ 时关于 $\Delta x$ 的高阶无穷小，则称 $f$ 在点 $x_0$ 可微，且称函数 $f$ 变量的线性主部 $\alpha \Delta x$ 为 $f$ 在 $x_0$ 处的微分。

1. 二元函数的全微分

在实际应用中，对于多元函数，也需要讨论函数改变量的线性主部问题。例如，一个底面半径为 $r$ 高为 $h$ 的圆柱体受热（冷）后膨胀（收缩），底半径 $r$ 与高 $h$ 各自产生改变量 $\Delta r$ 与 $\Delta h$，需要求体积 $V$ 的改变量 $\Delta V$。由于

Δ V &= π(r + Δ r)²(h + Δ h) - π r² h &= 2π r h Δ r + π r² Δ h + 2π r Δ r Δ h + π h (Δ r)² + π (Δ r)² Δ h,$$\begin{array}{rl}\Delta V& =\pi (r+\Delta r{)}^2(h+\Delta h)-\pi r^{2}h\\ & =2\pi rh\Delta r+\pi r^2\Delta h+2\pi r\Delta r\Delta h+\pi h(\Delta r{)}^2+\pi (\Delta r{)}^2\Delta h,\end{array}$$

所以当 $(\Delta r,\Delta h)\to (0,0)$ 时，易见上式右端的后三项之和是 $\rho =\sqrt{(\Delta r{)}^2+(\Delta h{)}^2}$ 的高阶无穷小，因此此时 $\Delta V$ 的值可用其余部分（即关于 $\Delta r$ 与 $\Delta h$ 的线性部分） $2\pi rh\Delta r+\pi r^2\Delta h$ 近似代替，这个近似值就是体积 $V$ 的改变量 $\Delta V$ 的线性主部。

由此，类似于一元函数微分定义，可以给出二元函数全微分的定义如下。

定义 3.2 (全微分)

设二元函数 $z=f(x,y)$ 在点 $(x_0,y_0)$ 的某邻域 $U(x_0,y_0)$ 内有定义。如果对于

Δ z = f(x₀ + Δ x, y₀ + Δ y) - f(x₀, y₀)$$\Delta z=f(x_0+\Delta x,y_0+\Delta y)-f(x_0,y_0)$$

可以表示为

Δ z = a₁ Δ x + a₂ Δ y + o(ρ), (3.6)$$\Delta z=a_1\Delta x+a_2\Delta y+o(\rho ),\text{\hspace{1em}(3.6)}$$

其中 $a_1,a_2$ 是与 $\Delta x,\Delta y$ 无关的两个常数（但一般与点 $(x_0,y_0)$ 有关），$\rho =\sqrt{(\Delta x{)}^2+(\Delta y{)}^2}$，$o(\rho )$ 是当 $\rho \to 0$（即 $\Delta x\to 0,\Delta y\to 0$）时关于 $\rho$ 的高阶无穷小，则称函数 $f$ 在点 $(x_0,y_0)$ 处可微，并称 $a_1\Delta x+a_2\Delta y$ 为函数 $f$ 在点 $(x_0,y_0)$ 处的全微分，记作

dz|_((x₀, y₀)), 或 df(x₀, y₀), 即$$\mathrm{d}z{|}_{(x_0,y_0)},\text{或}\mathrm{d}f(x_0,y_0),\text{即}$$dz|_((x₀, y₀)) = a₁ Δ x + a₂ Δ y. (3.7)$$\mathrm{d}z{|}_{(x_0,y_0)}=a_1\Delta x+a_2\Delta y.\text{\hspace{1em}(3.7)}$$

习惯上，将自变量的改变量 $\Delta x$ 与 $\Delta y$ 分别写成 $\mathrm{d}x$ 与 $\mathrm{d}y$，并分别称为自变量 $x$ 与 $y$ 的微分，所以函数 $f$ 的全微分也常写成

dz|_((x₀, y₀)) = a₁ dx + a₂ dy. (3.8)$$\mathrm{d}z{|}_{(x_0,y_0)}=a_1\mathrm{d}x+a_2\mathrm{d}y.\text{\hspace{1em}(3.8)}$$

由上述定义可见，当 $\rho$ 充分小且 $a_1$ 与 $a_2$ 不全为零时，全微分 $\mathrm{d}z{|}_{(x_0,y_0)}$ 就是函数 $f$ 在 $(x_0,y_0)$ 处改变量的线性主部。

上面我们给出了可微与全微分的定义，但 $f$ 究竟在什么条件下才可微呢？又当 $f$ 可微时，定义中的 $a_1,a_2$ 究竟等于什么，即如何计算全微分？函数的激微性与连续性及可偏导之间又有什么关系？对于这些涉及全微分的基本问题，自然是我们应当关心的，下面我们来讨论它们。

定理 3.1 (可微的必要条件)

设函数 $z=f(x,y)$ 在点 $(x_0,y_0)$ 处可微，则 (1) $f$ 在 $(x_0,y_0)$ 处连续； (2) $f$ 在 $(x_0,y_0)$ 处的两个偏导数均存在，且 $a_1=f_x(x_0,y_0),a_2=f_y(x_0,y_0)$，即

df(x₀, y₀) = fₓ(x₀, y₀) dx + f_y(x₀, y₀) dy. (3.9)$$\mathrm{d}f(x_0,y_0)=f_x(x_0,y_0)\mathrm{d}x+f_y(x_0,y_0)\mathrm{d}y.\text{\hspace{1em}(3.9)}$$

证明

(1) 当 $f$ 在 $(x_0,y_0)$ 处可微时，(3.6) 式成立，在 (3.6) 式中令 $\rho \to 0$（即 $\Delta x\to 0,\Delta y\to 0$），得

lim_(ρ → 0) Δ z = 0,$$\underset{\rho \to 0}{\lim }\Delta z=0,$$

或

lim_(substackΔ x → 0 Δ y → 0) f(x₀ + Δ x, y₀ + Δ y) = f(x₀, y₀),$$\underset{\begin{array}{c}\Delta x\to 0\\ \Delta y\to 0\end{array}}{\lim }f(x_0+\Delta x,y_0+\Delta y)=f(x_0,y_0),$$

所以 $f$ 在 $(x_0,y_0)$ 处连续。

(2) 由可微的定义，函数 $f(x,y)$ 在 $(x_0,y_0)$ 处的改变量可表示为

f(x₀ + Δ x, y₀ + Δ y) - f(x₀, y₀) = a₁ Δ x + a₂ Δ y + o(ρ),$$f(x_0+\Delta x,y_0+\Delta y)-f(x_0,y_0)=a_1\Delta x+a_2\Delta y+o(\rho ),$$

取 $\Delta y=0$，则有 $\rho =|\Delta x|$，上式变为

f(x₀ + Δ x, y₀) - f(x₀, y₀) = a₁ Δ x + o(|Δ x|).$$f(x_0+\Delta x,y_0)-f(x_0,y_0)=a_1\Delta x+o(|\Delta x|).$$

所以

lim_(Δ x → 0) (f(x₀ + Δ x, y₀) - f(x₀, y₀))/(Δ x) = lim_(Δ x → 0) [ a₁ + (o(|Δ x|))/(Δ x) ] = a₁,$$\underset{\Delta x\to 0}{\lim }\frac{f(x_0+\Delta x,y_0)-f(x_0,y_0)}{\Delta x}=\underset{\Delta x\to 0}{\lim }\left[a_1+\frac{o(|\Delta x|)}{\Delta x}\right]=a_1,$$

即在点 $(x_0,y_0)$ 处 $f(x,y)$ 对 $x$ 的偏导数存在且 $f_x(x_0,y_0)=a_1$。

同理，当取 $\Delta x=0$ 时，可证该点 $f(x,y)$ 对 $y$ 的偏导数也存在且 $f_y(x_0,y_0)=a_2$。

综上可得

df(x₀, y₀) = a₁ dx + a₂ dy = fₓ(x₀, y₀) dx + f_y(x₀, y₀) dy.$$\mathrm{d}f(x_0,y_0)=a_1\mathrm{d}x+a_2\mathrm{d}y=f_x(x_0,y_0)\mathrm{d}x+f_y(x_0,y_0)\mathrm{d}y.$$

由定理3.1 可知，当 $z=f(x,y)$ 在 $x_0=(x_0,y_0)$ 处可微时，$f$ 在 $x_0$ 处的改变量必可表示为

Δ z = fₓ(x₀) Δ x + f_y(x₀) Δ y + o(ρ), (3.10)$$\Delta z=f_x(x_0)\Delta x+f_y(x_0)\Delta y+o(\rho ),\text{\hspace{1em}(3.10)}$$

其中 $o(\rho )$ 是当 $\rho \to 0$ 时关于 $\rho$ 的高阶无穷小；反之，若 $\Delta z$ 可表示为 (3.10) 式，则由全微分的定义可知 $f$ 在 $x_0$ 处可微。

(3.9) 式给出了全微分的计算公式，由此可见，与一元函数的微分类似，二元函数 $f(x,y)$ 在点 $(x_0,y_0)$ 处的全微分不仅与 $f$ 在点 $(x_0,y_0)$ 处的各个偏导数有关，还与 $\Delta x,\Delta y$ 有关。如果 $z=f(x,y)$ 在区域 $\Omega \subseteq {\mathbf{R}}^2$ 的每一点都可微，则称 $f$ 是 $\Omega$ 内的可微函数，此时，点 $(x,y)$ 处的全微分可记成 $\mathrm{d}f$ 或 $\mathrm{d}z$，且有

dz = fₓ dx + f_y dy.$$\mathrm{d}z=f_x\mathrm{d}x+f_y\mathrm{d}y.$$

在一元函数中，我们知道函数在一点可导与可微是等价的。对于二元函数，函数在一点的所有偏导数均存在能否保证函数在此点可微呢？例3.3 已经给出了否定的答案。在例3.3 中，函数 $f$ 在该点却不连续，而由定理3.1 函数在一点连续是该点可微的必要条件，所以 $O(0,0)$ 处的两个偏导数均存在，但 $f$ 在该点却不连续，而由定理 3.1 函数在一点连续是该点可微的必要条件，所以 $f$ 在 $O$ 点不可微。这表明，对于多元函数，在一点的所有偏导数都存在，并不能保证函数在该点处可微。下面的例子显示：即使函数 $f$ 在一点连续而且所有的偏导数均存在也不足以保证其在该点可微。

例 3.5 讨论函数 $f(x,y)$ 在 $(0,0)$ 处的连续性与可微性

讨论函数

f(x,y) = (xy)/(√(x² + y²)), & x² + y² ≠ 0, 0, & x² + y² = 0$$f(x,y)=\{\begin{array}{ll}\frac{xy}{\sqrt{x^2+y^2}},& x^2+y^2\ne 0,\\ 0,& x^2+y^2=0\end{array}$$

在点 $O(0,0)$ 处的连续性与可微性。

证

由偏导数的定义易得 $f_x(0,0)=f_y(0,0)=0$。因此有

Δ f - [fₓ(0,0) Δ x + f_y(0,0) Δ y] = f(Δ x, Δ y) - f(0,0) = (Δ x² + Δ y²) sin (1)/(Δ x² + Δ y²) = ρ² sin (1)/(ρ²) = o(ρ),$$\Delta f-[f_x(0,0)\Delta x+f_y(0,0)\Delta y]=f(\Delta x,\Delta y)-f(0,0)=(\Delta x^2+\Delta y^2)\sin \frac{1}{\Delta x^2+\Delta y^2}={\rho }^2\sin \frac{1}{{\rho }^2}=o(\rho ),$$

故再由 (3.10) 式知 $f$ 在点 $(0,0)$ 处可微。下面证明它在点 $(0,0)$ 处的偏导数不连续，事实上，当 $x^2+y^2\ne 0$ 时，有

fₓ(x,y) = 2x sin (1)/(x² + y²) - (2x)/(x² + y²) cos (1)/(x² + y²),$$f_x(x,y)=2x\sin \frac{1}{x^2+y^2}-\frac{2x}{x^2+y^2}\cos \frac{1}{x^2+y^2},$$

由于 ${\lim }_{\begin{array}{c}x\to 0\\ y\to 0\end{array}}2x\sin \frac{1}{x^2+y^2}=0$，而

lim_(substacky → 0 x → 0) (2x)/(x² + y²) cos (1)/(x² + y²) = lim_(substackx → 0 y → 0) cos (1)/(2x²)$$\underset{\begin{array}{c}y\to 0\\ x\to 0\end{array}}{\lim }\frac{2x}{x^2+y^2}\cos \frac{1}{x^2+y^2}=\underset{\begin{array}{c}x\to 0\\ y\to 0\end{array}}{\lim }\cos \frac{1}{2x^2}$$

不存在，所以 $f_x(x,y)$ 在点 $(0,0)$ 处间断。同法可证 $f_y(x,y)$ 在点 $(0,0)$ 处也间断。

数字资源 5.3.1二元函数连续、可偏导及可微的关系 ↗

2. $n$ 元函数的全微分

二元函数全微分的概念可直接推广到 $n$ 元函数。设 $n$ 元函数 $u=f(x)=f(x_1,\dots ,x_n)$ 在点 $x_0=(x_{0,1},\dots ,x_{0,n})\in {\mathbf{R}}^n$ 的邻域 $U(x_0)\subseteq {\mathbf{R}}^n$ 内有定义，如果 $\forall \mathbf{x}=x_0+\Delta \mathbf{x}\in U(x_0)$，存在一组与 $\Delta \mathbf{x}=(\Delta x_1,\dots ,\Delta x_n)$ 无关的常数 $a_1,\dots ,a_n$，使得函数 $f$ 在 $x_0$ 处的改变量

Δ u = f(x₀ + Δ x) - f(x₀)$$\Delta u=f(x_0+\Delta \mathbf{x})-f(x_0)$$

可以表示为

Δ u = a₁ Δ x₁ + … + aₙ Δ xₙ + o(ρ),$$\Delta u=a_1\Delta x_1+\cdots +a_n\Delta x_n+o(\rho ),$$

其中 $o(\rho )$ 是当 $\rho =\Vert \Delta \mathbf{x}\Vert \to 0$ 时关于 $\rho$ 的高阶无穷小，则称 $f$ 在点 $x_0$ 可微，且称关于 $\Delta x_1,\dots ,\Delta x_n$ 的线性函数

a₁ Δ x₁ + … + aₙ Δ xₙ$$a_1\Delta x_1+\cdots +a_n\Delta x_n$$

为 $f$ 在 $x_0$ 处的全微分，记为 $\mathrm{d}f(x_0)$，或 $\mathrm{d}u{|}_{x=x_0}$，即

df(x₀) = a₁ Δ x₁ + … + aₙ Δ xₙ.$$\mathrm{d}f(x_0)=a_1\Delta x_1+\cdots +a_n\Delta x_n.$$

同二元函数一样，常记 $\Delta x_i=\mathrm{d}{x}_i$，于是 $f$ 的全微分也常写成

df(x₀) = ∑ᵢ=1ⁿ aᵢ dxᵢ. (3.15)$$\mathrm{d}f(x_0)=\sum _{i=1}^n{a}_i\mathrm{d}{x}_i.\text{\hspace{1em}(3.15)}$$

此外，若 $f$ 在点 $x_0$ 处的所有偏导数均连续，则 $f$ 在点 $x_0$ 处一定可微。

例 3.8 求 $f(x,y,z)=\frac{1}{\sqrt{x^2 + y^2 + z^2}}$ 的全微分

求 $f(x,y,z)=\frac{1}{\sqrt{x^2+y^2+z^2}}$ 的全微分。

解

显然，在任意 $(x,y,z)\ne (0,0,0)$ 处，$f$ 的所有偏导数均存在且连续，因此由定理3.2 知 $f$ 可微，再由 (3.15) 式得

df = fₓ dx + f_y dy + f_z dz = -fracx dx + y dy + z dz(x² + y² + z²)³/2 (x² + y² + z² ≠ 0).$$\mathrm{d}f=f_x\mathrm{d}x+f_y\mathrm{d}y+f_z\mathrm{d}z=-\frac{x\mathrm{d}x+y\mathrm{d}y+z\mathrm{d}z}{(x^2+y^2+z^2{)}^{3/2}}(x^2+y^2+z^2\ne 0).$$

3. 全微分在近似计算和误差估计中的应用

我们知道，当 $n$ 元函数 $f$ 在点 $x_0$ 可微时，有

Δ f = f(x₀ + Δ x) - f(x₀) = df(x₀) + o(ρ),$$\Delta f=f(x_0+\Delta \mathbf{x})-f(x_0)=\mathrm{d}f(x_0)+o(\rho ),$$

从而当 $\rho =\Vert \Delta \mathbf{x}\Vert \ll 1$ 时，有

f(x₀ + Δ x) - f(x₀) ≈ df(x₀) = ∑ᵢ=1ⁿ fₓᵢ(x₀) Δ xᵢ,$$f(x_0+\Delta \mathbf{x})-f(x_0)\approx \mathrm{d}f(x_0)=\sum _{i=1}^n{f}_{x_i}(x_0)\Delta x_i,$$

或

f(x) ≈ f(x₀) + ∑ᵢ=1ⁿ fₓᵢ(x₀)(xᵢ - x₀,i). (3.16)$$f(\mathbf{x})\approx f(x_0)+\sum _{i=1}^n{f}_{x_i}(x_0)(x_i-x_{0,i}).\text{\hspace{1em}(3.16)}$$

可见与一元函数情况类似，当 $\rho \ll 1$ 时，用函数的全微分近似表示函数的改变量，就是在点 $x_0$ 的邻域内把一个非线性函数 $f(\mathbf{x})$ 线性化。 (3.16) 式右端称为 $f(\mathbf{x})$ 的局部线性化或局部线性逼近。这种局部线性逼近的思想可以用来解决以下两类问题：

(1) 函数数值的近似计算

例 3.9 求 $\sqrt{1.97^3 + 1.01^3}$ 的近似值

求 $\sqrt{1.97^3+1.01^3}$ 的近似值。

解

令 $f(x,y)=\sqrt{x^3+y^3}$，$(x_0,y_0)=(2,1)$，$\Delta x=-0.03,\Delta y=0.01$，则

fₓ(2,1) = (3x²)/(2√(x³ + y³)) g|_((2,1)) = 2, f_y(2,1) = (3y²)/(2√(x³ + y³)) g|_((2,1)) = 1/2,$$f_x(2,1)=\frac{3x^2}{2\sqrt{x^3+y^3}}{|}_{(2,1)}=2,f_y(2,1)=\frac{3y^2}{2\sqrt{x^3+y^3}}{|}_{(2,1)}=\frac{1}{2},$$√(1.97³ + 1.01³) = f(x₀ + Δ x, y₀ + Δ y) ≈ f(x₀, y₀) + df(x₀)$$\sqrt{1.97^3+1.01^3}=f(x_0+\Delta x,y_0+\Delta y)\approx f(x_0,y_0)+\mathrm{d}f(x_0)$$= f(2,1) + fₓ(2,1) Δ x + f_y(2,1) Δ y = 2.945.$$=f(2,1)+f_x(2,1)\Delta x+f_y(2,1)\Delta y=\mathrm{2.945.}$$

(2) 误差估计

设量 $z$ 由公式 $z=f(x,y)$ 确定，如果 $x$ 和 $y$ 的近似值 $x_0$ 和 $y_0$ 分别有绝对误差 ${\delta }_x$ 和 ${\delta }_y$，即 $|\Delta x|<{\delta }_x,|\Delta y|<{\delta }_y$，那么近似值 $x_0$ 和 $y_0$ 代入公式 $z=f(x,y)$ 所得的值 $z_0=f(x_0,y_0)$ 也是 $z$ 的近似值。由于 $|\Delta x|$ 和 $|\Delta y|$ 都很小，所以可用公式 $\Delta z\approx \mathrm{d}z$ 估计 $z_0$ 的绝对误差

|Δ z| &≈ |dz| = |fₓ(x₀, y₀) Δ x + f_y(x₀, y₀) Δ y| &≤ |fₓ(x₀, y₀)| |Δ x| + |f_y(x₀, y₀)| |Δ y| &< |fₓ(x₀, y₀)| δₓ + |f_y(x₀, y₀)| δ_y,$$\begin{array}{rl}|\Delta z|& \approx |\mathrm{d}z|=|f_x(x_0,y_0)\Delta x+f_y(x_0,y_0)\Delta y|\\ & \le |f_x(x_0,y_0)||\Delta x|+|f_y(x_0,y_0)||\Delta y|\\ & <|f_x(x_0,y_0)|{\delta }_x+|f_y(x_0,y_0)|{\delta }_y,\end{array}$$

从而可取 $z_0$ 的绝对误差为

δ_z = |fₓ(x₀, y₀)| δₓ + |f_y(x₀, y₀)| δ_y. (3.17)$${\delta }_z=|f_x(x_0,y_0)|{\delta }_x+|f_y(x_0,y_0)|{\delta }_y.\text{\hspace{1em}(3.17)}$$

由此得 $z_0$ 的相对误差为

(δ_z)/(|z₀|) = | (fₓ(x₀, y₀))/(f(x₀, y₀)) | δₓ + | (f_y(x₀, y₀))/(f(x₀, y₀)) | δ_y. (3.18)$$\frac{{\delta }_z}{|z_0|}=\left|\frac{f_x(x_0,y_0)}{f(x_0,y_0)}\right|{\delta }_x+\left|\frac{f_y(x_0,y_0)}{f(x_0,y_0)}\right|{\delta }_y.\text{\hspace{1em}(3.18)}$$

例 3.10 设 $z=xy$，求由测量值 $x_0, y_0$ 计算 $z$ 所产生的绝对误差与相对误差

设 $z=xy$，求由测量值 $x_0,y_0$ 计算 $z$ 所产生的绝对误差与相对误差。

解

将 $z_x=y,z_y=x$ 代入公式 (3.17) 与 (3.18)，得 绝对误差 ${\delta }_z=|y_0|{\delta }_x+|x_0|{\delta }_y$, 相对误差 $\frac{{\delta }_z}{|z_0|}=\frac{{\delta }_x}{|x_0|}+\frac{{\delta }_y}{|y_0|}$. 由此可见，乘积的相对误差等于各个因子的相对误差之和。同理可证，商的相对误差等于分子与分母的相对误差之差。

例 3.11 肾脏清除率误差估计

肾的一个重要功能是清除血液中的尿素。临床上用公式 $C=\frac{\sqrt{V}}{P}u$ 来计算尿素标准清除率，其中 $u$ 表示尿素中的尿素浓度（单位: $\text{mg/L}$），$V$ 表示每分钟排出的尿量（单位: $\text{mL/min}$），$P$ 表示血液中的尿素浓度（单位: $\text{mg/L}$）。某患者的实际测量值 $u=5000,V=1.44,P=200$，从而算得 $C=30$（正常人的 $C$ 约为 $54$）。如果该测量值 $u,V,P$ 的绝对误差分别为 $50,0.014,4.2$，试估算由测量值的误差对 $C$ 值所带来的绝对误差和相对误差。

解

因 $\frac{\partial C}{\partial u}=\frac{\sqrt{V}}{P},\frac{\partial C}{\partial V}=\frac{u}{2P\sqrt{V}},\frac{\partial C}{\partial P}=-\frac{u\sqrt{V}}{P^2}$，当 $u=5000,V=1.44,P=200$ 时，

(∂ C)/(∂ u) = 0.006, (∂ C)/(∂ V) = (125)/(2), (∂ C)/(∂ P) = -0.15,$$\frac{\partial C}{\partial u}=0.006,\frac{\partial C}{\partial V}=\frac{125}{2},\frac{\partial C}{\partial P}=-0.15,$$

3.3 方向导数与梯度

1. 方向导数

我们知道二元函数 $f(x,y)$ 在点 $(x_0,y_0)$ 处的两个偏导数，实际上是刻画函数 $f$ 在点 $(x_0,y_0)$ 处沿 $x$ 轴和 $y$ 轴平行的方向的变化率。然而，在许多实际问题中，仅仅研究函数沿这两个特殊方向的变化率还远远不够，还需要研究函数在点 $(x_0,y_0)$ 处沿各个不同方向的变化率。例如，在大气气象学中，为了确定某地的风速和风向，就要研究该地气压沿不同方向变化的快慢；在热传导问题中，需要研究某点温度沿不同方向变化的快慢等。这些研究函数在某点处沿不同方向变化率的问题就是方向导数问题。

设点 ${\mathbf{x}}_0\in {\mathbf{R}}^2$, $l$ 是平面上某一向量，其单位向量记为 ${\mathbf{e}}_l$，$f:U({\mathbf{x}}_0)\subseteq {\mathbf{R}}^2\to \mathbf{R}$ 是一个二元函数。我们来讨论 $f$ 在点 ${\mathbf{x}}_0$ 处沿 $l$ 方向的变化率（记作 ${\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}$）过点 ${\mathbf{x}}_0$ 作与 $l$ 平行的直线 $L$（图 5.11），它的方程为

x = x₀ + teₗ, t ∈ R. (3.19)$$\mathbf{x}={\mathbf{x}}_0+t{\mathbf{e}}_l,t\in \mathbf{R}.\text{\hspace{1em}(3.19)}$$

$f(\mathbf{x})$ 在点 ${\mathbf{x}}_0$ 处沿 $l$ 方向的变化率，就是当 $\mathbf{x}$ 在直线 $L$ 上变动时 $f(\mathbf{x})$ 在点 ${\mathbf{x}}_0$ 处的变化率。在点 ${\mathbf{x}}_0$ 与 ${\mathbf{e}}_l$ 固定的情况下，当点 $\mathbf{x}$ 在 $L$ 上变动时，函数

f(x) = f(x₀ + teₗ)$$f(\mathbf{x})=f({\mathbf{x}}_0+t{\mathbf{e}}_l)$$

实际上是自变量 $t$ 的一元函数，记作

F(t) = f(x₀ + teₗ).$$F(t)=f({\mathbf{x}}_0+t{\mathbf{e}}_l).$$

因此 $f(\mathbf{x})$ 在点 ${\mathbf{x}}_0$ 沿 $l$ 方向的变化率，也就是一元函数 $F(t)$ 在 $t=0$ 处的导数，即

. (∂ f)/(∂ l) |ₓ₀ = . (dF(t))/(dt) |ₜ=0 = limₜ → 0 (F(t) - F(0))/(t) = limₜ → 0 (f(x₀ + teₗ) - f(x₀))/(t).$${\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}={\frac{\mathrm{d}F(t)}{\mathrm{d}t}|}_{t=0}=\underset{t\to 0}{\lim }\frac{F(t)-F(0)}{t}=\underset{t\to 0}{\lim }\frac{f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)}{t}.$$

思路分析

多元函数 $f(\mathbf{x})$ 在点 ${\mathbf{x}}_0$ 沿 $l$ 方向的方向导数是通过引入参数 $t$ 将该函数转化为一元函数 $F(t)=f({\mathbf{x}}_0+t{\mathbf{e}}_l)$，用 $F(t)$ 在 $t=0$ 处的导数来定义的。这体现了“化多为一”的思想：用已知的（一元函数导数定义）去研究未知的（多元函数方向导数）。

定义 3.3 (方向导数)

设点 ${\mathbf{x}}_0\in {\mathbf{R}}^2$, $l$ 是平面上一向量，与 $l$ 同向的单位向量为 ${\mathbf{e}}_l$，二元函数 $f$ 定义在 ${\mathbf{x}}_0$ 的邻域 $U({\mathbf{x}}_0)\subseteq {\mathbf{R}}^2$ 内，在 $U({\mathbf{x}}_0)$ 内让自变量 $\mathbf{x}$ 由 ${\mathbf{x}}_0$ 沿与 ${\mathbf{e}}_l$ 平行的直线变到 ${\mathbf{x}}_0+t{\mathbf{e}}_l$，从而函数值有对应的改变量 $f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)$。若极限

limₜ → 0 (f(x₀ + teₗ) - f(x₀))/(t)$$\underset{t\to 0}{\lim }\frac{f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)}{t}$$

存在，则称此极限值为 $f$ 在点 ${\mathbf{x}}_0$ 处沿 $l$ 方向的方向导数，记作 ${\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}$ 或 $\frac{\partial f({\mathbf{x}}_0)}{\partial l}$，即

. (∂ f)/(∂ l) |ₓ₀ = limₜ → 0 (f(x₀ + teₗ) - f(x₀))/(t). (3.20)$${\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}=\underset{t\to 0}{\lim }\frac{f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)}{t}.\text{\hspace{1em}(3.20)}$$

我们对方向导数的定义再作如下说明：

(1) 定义中的 $t$ 的绝对值就是两点 ${\mathbf{x}}_0$ 与 $\mathbf{x}={\mathbf{x}}_0+t{\mathbf{e}}_l$ 之间的距离 $d$。事实上，注意到 ${\mathbf{e}}_l$ 为单位向量，是有

d = | (x₀ + teₗ) - x₀ | = | teₗ | = |t| |eₗ| = |t|.$$d=\Vert ({\mathbf{x}}_0+t{\mathbf{e}}_l)-{\mathbf{x}}_0\Vert =\Vert t{\mathbf{e}}_l\Vert =|t|\Vert {\mathbf{e}}_l\Vert =|t|.$$

(2) 方向导数实际上是函数 $f$ 在 ${\mathbf{x}}_0$ 处沿 $l$ 方向关于距离的变化率，它表示当 $\mathbf{x}$ 在 ${\mathbf{x}}_0$ 邻近沿 $l$ 方向变化（无论 $t>0$ 还是 $t<0$）时，函数 $f(\mathbf{x})$ 变化的快慢程度。事实上，当 $t>0$ 时，由点 ${\mathbf{x}}_0$ 到点 ${\mathbf{x}}_0+t{\mathbf{e}}_l$ 的向量 $t{\mathbf{e}}_l$ 与 $l$ 同向，所以此时 $f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)$ 就是 $f$ 在点 ${\mathbf{x}}_0$ 处沿 $l$ 方向变到 ${\mathbf{x}}_0+t{\mathbf{e}}_l$ 的改变量，又因 $t=d$，故

(f(x₀ + teₗ) - f(x₀))/(t) = (f(x₀ + teₗ) - f(x₀))/(d)$$\frac{f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)}{t}=\frac{f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)}{d}$$

表示当动点由 ${\mathbf{x}}_0$ 沿 $l$ 方向变到 ${\mathbf{x}}_0+t{\mathbf{e}}_l$ 时，函数 $f$ 关于距离的平均变化率，所以，方向导数 ${\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}$ 表示函数 $f$ 在点 ${\mathbf{x}}_0$ 处关于距离的变化率。而当 $t<0$ 时，由点 ${\mathbf{x}}_0+t{\mathbf{e}}_l$ 到点 ${\mathbf{x}}_0$ 的向量 $-t{\mathbf{e}}_l$ 与 $l$ 同向，又因 $d=-t$，于是由

(f(x₀ + teₗ) - f(x₀))/(t) = (f(x₀) - f(x₀ + teₗ))/(-t) = (f(x₀) - f(x₀ + teₗ))/(d)$$\frac{f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)}{t}=\frac{f({\mathbf{x}}_0)-f({\mathbf{x}}_0+t{\mathbf{e}}_l)}{-t}=\frac{f({\mathbf{x}}_0)-f({\mathbf{x}}_0+t{\mathbf{e}}_l)}{d}$$

可知方向导数 (3.20) 仍然是 $f$ 在 ${\mathbf{x}}_0$ 处沿 $l$ 方向关于距离的变化率。

(3) 由 (2) 可知，若方向导数 ${\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}>0$ ($<0$)，则 $f$ 在 ${\mathbf{x}}_0$ 处沿 $l$ 方向增加（减少）。而且，方向导数的绝对值愈大，则表明 $f$ 在 ${\mathbf{x}}_0$ 处沿 $l$ 方向增加（减少）得愈快。

2. 方向导数的几何意义

过直线 $L:\mathbf{x}={\mathbf{x}}_0+t{\mathbf{e}}_l$ 作平行于 $z$ 轴的平面 $\pi$（图 5.12），它与曲面 $z=f(x,y)$ 所交的曲线记作 $C$。在平面 $\pi$ 上考察曲线 $C$。容易看出，当 $t>0$ 时，$\overrightarrow{P_{0}P}$ 的方向与 $l$ 的方向相应，

(f(x₀ + teₗ) - f(x₀))/(t)$$\frac{f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)}{t}$$

表示曲线 $C$ 的割线向量 $\overrightarrow{P_{0}P}$ 与向量 $l$ 交角的正切值，即 $\overrightarrow{P_{0}P}$（关于 $l$ 方向）的斜率；当 $t<0$ 时，$P$ 在曲线 $C$ 上位于 $P_0$ 的另一侧，

(f(x₀ + teₗ) - f(x₀))/(t) = (f(x₀) - f(x₀ + teₗ))/(-t)$$\frac{f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)}{t}=\frac{f({\mathbf{x}}_0)-f({\mathbf{x}}_0+t{\mathbf{e}}_l)}{-t}$$

表示割线向量 $\overrightarrow{P{P}_0}$ 的斜率。当 $t\to 0$ 时，$\mathbf{x}\to {\mathbf{x}}_0$，割线转化为切线。

(3.20) 式中极限的存在，意味着当点 $\mathbf{x}$ 沿向量 $l$ 从 ${\mathbf{x}}_0$ 的两侧分别趋于 ${\mathbf{x}}_0$ 时，极限存在且相等，即曲线 $C$ 在点 $P_0$ 仅有唯一的切线 $T$，它关于 $l$ 方向的斜率就是方向导数 ${\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}$。

相应于一元函数的左（右）导数，对于方向导数，也可以按 $t\to 0^{+}$ 与 $t\to 0^{-}$ 讨论单侧方向导数。

图 5.12 方向导数的几何意义示意图

例 3.12 求二元函数的方向导数

设二元函数

f(x, y) = (xy²)/(x² + y⁴), & x² + y² ≠ 0, 0, & x² + y² = 0,$$f(x,y)=\{\begin{array}{ll}\frac{x{y}^2}{x^2+y^4},& x^2+y^2\ne 0,\\ 0,& x^2+y^2=0,\end{array}$$

求 $f$ 在点 $(0,0)$ 沿 ${\mathbf{e}}_l=(\cos \theta ,\sin \theta )$ 的方向导数。

解

当 $\cos \theta \ne 0$ 时，有

. (∂ f(0,0))/(∂ l) | = limₜ → 0 (f(t cos θ, t sin θ) - f(0,0))/(t) = limₜ → 0 (cos θ · sin² θ)/(cos² θ + t² sin⁴ θ) = (sin² θ)/(cos θ).$$\frac{\partial f(0,0)}{\partial l}|=\underset{t\to 0}{\lim }\frac{f(t\cos \theta ,t\sin \theta )-f(0,0)}{t}=\underset{t\to 0}{\lim }\frac{\cos \theta \cdot {\sin }^2\theta }{{\cos }^2\theta +t^2{\sin }^4\theta }=\frac{{\sin }^2\theta }{\cos \theta }.$$

当 $\cos \theta =0$ 时，由于 $f(t\cos \theta ,t\sin \theta )-f(0,0)=0$，从而 $\frac{\partial f(0,0)}{\partial l}|=0$。

由例3.12 可知，当 $\theta =\frac{\pi }{4}$ 时，$\frac{\partial f(0,0)}{\partial l}|=\frac{\sqrt{2}}{2}$；当 $\theta =\frac{\pi }{4}+\pi$ 时，$\frac{\partial f(0,0)}{\partial l}|=-\frac{\sqrt{2}}{2}$。可见，$f$ 在 $(0,0)$ 处沿 $\theta =\frac{\pi }{4}$ 的方向导数与沿 $\theta =\frac{\pi }{4}+\pi$ 的方向导数的绝对值相等但符号相反。一般地，由方向导数的定义容易看出，

. (∂ f)/(∂ (-l)) |ₓ₀ = -. (∂ f)/(∂ l) |ₓ₀.$${\frac{\partial f}{\partial (-l)}|}_{{\mathbf{x}}_0}=-{\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}.$$

即 $f$ 在 ${\mathbf{x}}_0$ 处沿 $l$ 方向的方向导数与沿 $-l$ 方向的方向导数只差一个符号。因此，若 $f$ 在 ${\mathbf{x}}_0$ 处沿 $l$ 方向是增加（减少）的，则沿 $-l$ 方向就是减少（增加）的。

数字资源 5.3.2方向导数与偏导数的区别与联系 ↗

3. 梯度

函数 $z=f(x,y)$ 在点 $(x_0,y_0)$ 处沿某方向的方向导数刻画了 $f(x,y)$ 在点 $(x_0,y_0)$ 处沿该方向的变化率。一般来说，在给定点沿不同方向的变化率（方向导数）是不相同的，而且在给定点可能存在无数个不同的方向。那么沿无数个方向的方向导数中，沿哪个方向的方向导数最大，怎样求这个最大的方向导数的值（如果存在的话），在科学技术中是一个非常重要的问题。例如，当热流从热源向周围扩散时，常常需要知道温度变化最快的方向以及沿此方向温度的变化率。为研究此类问题，就需要引入梯度的概念。

定义 3.4 (梯度)

设二元函数 $z=f(x,y)$ 定义在点 $(x_0,y_0)$ 的邻域中，如果存在一个向量，其方向为该函数在此点取得方向导数最大值的方向，其模等于该函数在此点方向导数的最大值，则称该向量为函数 $f(x,y)$ 在点 $(x_0,y_0)$ 处的梯度，记作 $\text{grad}f(x_0,y_0)$，其中 $\text{grad}$ 是英文 $\text{gradient}$ 的简称。

利用方向导数的计算公式 (3.21) 和梯度的定义，容易得到梯度存在的一个充分条件和梯度的计算公式。

定理 3.4

设二元函数 $f(x,y)$ 在点 $(x_0,y_0)$ 处可微，则 $f$ 在该点的梯度一定存在，并且

grad f(x₀, y₀) = fₓ(x₀, y₀)i + f_y(x₀, y₀)j. (3.23)$$\text{grad}f(x_0,y_0)=f_x(x_0,y_0)\mathbf{i}+f_y(x_0,y_0)\mathbf{j}.\text{\hspace{1em}(3.23)}$$

证明

设 $l$ 为任一向量，由于函数 $f(x,y)$ 在点 $(x_0,y_0)$ 处可微，根据定理3.3，$f$ 在 $(x_0,y_0)$ 处沿 ${\mathbf{e}}_l=(\cos \alpha ,\cos \beta )$ 的方向导数必定存在，并且

. (∂ f)/(∂ l) |_((x₀, y₀)) = fₓ(x₀, y₀) cos α + f_y(x₀, y₀) cos β.$${\frac{\partial f}{\partial l}|}_{(x_0,y_0)}=f_x(x_0,y_0)\cos \alpha +f_y(x_0,y_0)\cos \beta .$$

作向量 $\mathbf{g}=(f_x(x_0,y_0),f_y(x_0,y_0))$，下面证明 $\mathbf{g}$ 就是 $f$ 在点 $(x_0,y_0)$ 处的梯度。事实上，由于

. (∂ f)/(∂ l) |_((x₀, y₀)) = langle g, eₗ rangle = |g| cos(g, eₗ) ≤ |g|,$${\frac{\partial f}{\partial l}|}_{(x_0,y_0)}=⟨\mathbf{g},{\mathbf{e}}_l⟩=\Vert \mathbf{g}\Vert \cos (\mathbf{g},{\mathbf{e}}_l)\le \Vert \mathbf{g}\Vert ,$$

所以，当 $\cos (\mathbf{g},{\mathbf{e}}_l)=1$ 时，即 ${\mathbf{e}}_l$ 与 $\mathbf{g}$ 同向时，

. (∂ f)/(∂ l) |_((x₀, y₀)) = |g|.$${\frac{\partial f}{\partial l}|}_{(x_0,y_0)}=\Vert \mathbf{g}\Vert .$$

也就是说，向量 $\mathbf{g}$ 的方向是函数 $f$ 在点 $(x_0,y_0)$ 处取得方向导数最大值的方向，$\mathbf{g}$ 的模 $\Vert \mathbf{g}\Vert$ 就是 $f$ 在此点方向导数的最大值。因此，

g = fₓ(x₀, y₀)i + f_y(x₀, y₀)j = grad f(x₀, y₀).$$\mathbf{g}=f_x(x_0,y_0)\mathbf{i}+f_y(x_0,y_0)\mathbf{j}=\text{grad}f(x_0,y_0).$$

注意

上式左端表示 $f$ 在点 $(x_0,y_0)$ 处沿 $l$ 方向的方向导数等于向量 $\mathbf{g}$ 在 $l$ 方向的投影，即方向导数 ${\frac{\partial f}{\partial l}|}_{(x_0,y_0)}$ 等于梯度在 $l$ 方向的投影。

函数在 $(x_0,y_0)$ 处的梯度也可以用 $\text{nabla}$ 符号 $\nabla$（也称向量微分算子）来表示，记作 $\nabla f(x_0,y_0)$，其中符号 $\nabla =\left(\frac{\partial }{\partial x},\frac{\partial }{\partial y}\right)$，它本身并无意义，仅当 $\nabla$ 作用于函数 $f$ 后表示如下向量：

∇ f(x₀, y₀) = ( (∂ f(x₀, y₀))/(∂ x), (∂ f(x₀, y₀))/(∂ y) ).$$\nabla f(x_0,y_0)=\left(\frac{\partial f(x_0,y_0)}{\partial x},\frac{\partial f(x_0,y_0)}{\partial y}\right).$$

例 3.13 求梯度与方向导数

求二元函数 $u=x^2-xy+y^2$ 在点 $P(-1,1)$ 处沿方向 ${\mathbf{e}}_l=\frac{1}{\sqrt{5}}(2,1)$ 的方向导数。并指出 $u$ 在该点沿哪个方向的方向导数最大？这个最大的方向导数值是多少？$u$ 沿哪个方向减小得最快？沿哪个方向 $u$ 的值不变？

解

解此题的关键是求函数 $u$ 在 $(-1,1)$ 点处的梯度。由于

∇ u |_((-1,1)) = ( (∂ u)/(∂ x), (∂ u)/(∂ y) ) |_((-1,1)) = (2x - y, 2y - x) |_((-1,1)) = (-3, 3),$$\nabla u{|}_{(-1,1)}=\left(\frac{\partial u}{\partial x},\frac{\partial u}{\partial y}\right){|}_{(-1,1)}=(2x-y,2y-x){|}_{(-1,1)}=(-3,3),$$

所以

. (∂ u)/(∂ l) |_((-1,1)) = langle ∇ u|_((-1,1)), eₗ rangle = (1)/(√(5))(-6 + 3) = (-3)/(√(5)).$${\frac{\partial u}{\partial l}|}_{(-1,1)}=⟨\nabla u{|}_{(-1,1)},{\mathbf{e}}_l⟩=\frac{1}{\sqrt{5}}(-6+3)=\frac{-3}{\sqrt{5}}.$$

方向导数取得最大值的方向即梯度方向，其单位向量为 $\frac{1}{\sqrt{2}}(-1,1)$，方向导数的最大值 $\Vert \nabla u{|}_{(-1,1)}\Vert =3\sqrt{2}$。$u$ 沿梯度的负方向，即 $\frac{1}{\sqrt{2}}(1,-1)$ 的方向减小得最快。为使 $u$ 的值不变的方向，就是求使 $u$ 的变化率为零的方向，令 ${\mathbf{e}}_l=(\cos \theta ,\sin \theta )$，则

. (∂ u)/(∂ l) |_((-1,1)) = langle ∇ u|_((-1,1)), eₗ rangle = -3cos θ + 3sin θ = 3√(2) sin(θ - (π)/(4)).$${\frac{\partial u}{\partial l}|}_{(-1,1)}=⟨\nabla u{|}_{(-1,1)},{\mathbf{e}}_l⟩=-3\cos \theta +3\sin \theta =3\sqrt{2}\sin \left(\theta -\frac{\pi }{4}\right).$$

令 $\frac{\partial u}{\partial l}=0$，得 $\theta =\frac{\pi }{4}$ 或 $\theta =\pi +\frac{\pi }{4}$，故在点 $(-1,1)$ 处沿 $\theta =\frac{\pi }{4}$ 和 $\theta =\pi +\frac{\pi }{4}$ 的方向，函数 $u$ 的值不变。

图 5.13 等值线与梯度方向的关系

从图 5.13 的等值线来看，例3.13 的结论是显而易见的。事实上，函数 $u$ 是如图 5.13 所示的一簇椭圆。$u=x^2-xy+y^2$ 的等值线 $x^2-xy+y^2=C$ 是如图 5.13 所示的一簇椭圆。$u$ 在点 $P(-1,1)$ 处的梯度方向为 $\overrightarrow{P{P}_1}$，它垂直于等值线在点 $P$ 的切线，函数 $u$ 在点 $P$ 沿此方向的变化率显然最大，沿与 $\overrightarrow{P{P}_1}$ 相反的方向的变化率最小，而在点 $P$ 使 $u$ 的变化率为零的方向当然应是等值线在点 $P$ 的切线方向。

4. $n$ 元函数的方向导数与梯度

方向导数与梯度的概念均可直接推广到 $n$ 元函数。设 ${\mathbf{x}}_0$ 是 ${\mathbf{R}}^n$ 中一个点，$f$ 为定义在 ${\mathbf{x}}_0$ 的邻域 $U({\mathbf{x}}_0)$ 内的一个 $n$ 元函数，${\mathbf{e}}_l$ 为 ${\mathbf{R}}^n$ 中一个单位向量，则 $n$ 元函数 $f(\mathbf{x})$ 在点 ${\mathbf{x}}_0$ 沿 ${\mathbf{e}}_l$ 方向的方向导数定义与二元函数完全类似，即

. (∂ f)/(∂ l) |ₓ₀ = limₜ → 0 (f(x₀ + teₗ) - f(x₀))/(t).$${\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}=\underset{t\to 0}{\lim }\frac{f({\mathbf{x}}_0+t{\mathbf{e}}_l)-f({\mathbf{x}}_0)}{t}.$$

当 $f$ 在点 ${\mathbf{x}}_0$ 可微时，$f$ 在点 ${\mathbf{x}}_0$ 沿任意 $l$ 方向的方向导数都存在，并且有如下的计算公式：

. (∂ f)/(∂ l) |ₓ₀ = ∑ₖ=1ⁿ (∂ f(x₀))/(∂ xₖ) cos αₖ,$${\frac{\partial f}{\partial l}|}_{{\mathbf{x}}_0}=\sum _{k=1}^n\frac{\partial f({\mathbf{x}}_0)}{\partial x_k}\cos {\alpha }_k,$$

其中 $\cos {\alpha }_1,\cos {\alpha }_2,\dots ,\cos {\alpha }_n$ 为向量 $l$ 的方向余弦。

当 $f$ 在点 ${\mathbf{x}}_0$ 处可微时，$f$ 在 ${\mathbf{x}}_0$ 的梯度也存在，计算公式为

grad f(x₀) = ∇ f(x₀) = ( (∂ f)/(∂ x₁), (∂ f)/(∂ x₂), …, (∂ f)/(∂ xₙ) )ₓ₀.$$\text{grad}f({\mathbf{x}}_0)=\nabla f({\mathbf{x}}_0)={\left(\frac{\partial f}{\partial x_1},\frac{\partial f}{\partial x_2},\dots ,\frac{\partial f}{\partial x_n}\right)}_{{\mathbf{x}}_0}.$$

5. 梯度的运算法则

根据梯度的计算公式 (3.23) 可见，求函数 $u=f(\mathbf{x})$ 的梯度实际上就是求偏导数。故由已知的求导法则，可以得知梯度具有类似于求导法则的一些简单运算法则（其中的 $C_1,C_2$ 为任意常数，函数 $u,v$ 及 $f$ 均可微）：

(1) $\text{grad}(C_{1}u+C_{2}v)=C_1\text{grad }u+C_2\text{grad }v$，或 $\nabla (C_{1}u+C_{2}v)=C_1\nabla u+C_2\nabla v$;

(2) $\text{grad}(uv)=u\text{grad }v+v\text{grad }u$，或 $\nabla (uv)=u\nabla v+v\nabla u$;

(3) $\text{grad}\left(\frac{u}{v}\right)=\frac{1}{v^2}[v\text{grad }u-u\text{grad }v]$，或 $\nabla \left(\frac{u}{v}\right)=\frac{1}{v^2}[v\nabla u-u\nabla v](v\ne 0)$;

(4) $\text{grad }f(u)=f^{\prime }(u)\text{grad }u$，或 $\nabla f(u)=f^{\prime }(u)\nabla u$.

仅证 (4)，其余法则的证明留给读者。设 $u=u(\mathbf{x})=u(x_1,\dots ,x_n)$，由一元函数的链式法则，我们有

∇ f(u) = ( (∂ f(u))/(∂ x₁), …, (∂ f(u))/(∂ xₙ) ) = ( f'(u) (∂ u)/(∂ x₁), …, f'(u) (∂ u)/(∂ xₙ) ) = f'(u) ( (∂ u)/(∂ x₁), …, (∂ u)/(∂ xₙ) ) = f'(u) ∇ u.$$\nabla f(u)=\left(\frac{\partial f(u)}{\partial x_1},\dots ,\frac{\partial f(u)}{\partial x_n}\right)=\left(f^{\prime }(u)\frac{\partial u}{\partial x_1},\dots ,f^{\prime }(u)\frac{\partial u}{\partial x_n}\right)=f^{\prime }(u)\left(\frac{\partial u}{\partial x_1},\dots ,\frac{\partial u}{\partial x_n}\right)=f^{\prime }(u)\nabla u.$$

例 3.14 电场强度与电位

设点电荷 $q$ 放在坐标原点 $O(0,0,0)$ 处，则在其周围会产生电场，且任一点 $M(x,y,z)$ 的电位与电场强度分别为

u = (q)/(4πε r), E = (q)/(4πε r³) r (r ≠ 0),$$u=\frac{q}{4\pi \epsilon r},\mathbf{E}=\frac{q}{4\pi \epsilon r^3}\mathbf{r}(r\ne 0),$$

其中 $\epsilon$ 为介电系数，$r$ 为点 $M$ 的向径，$r=\Vert \mathbf{r}\Vert$，求电位函数 $u$ 的梯度。

解

解 $\nabla u=\nabla \left(\frac{q}{4\pi \epsilon r}\right)=\frac{q}{4\pi \epsilon }\nabla \left(\frac{1}{r}\right)=\frac{q}{4\pi \epsilon }{\left(\frac{1}{r}\right)}^{\prime }\nabla r=-\frac{q}{4\pi \epsilon r^2}\nabla r$,

而 $\nabla r=\left(\frac{\partial r}{\partial x},\frac{\partial r}{\partial y},\frac{\partial r}{\partial z}\right)=\left(\frac{x}{r},\frac{y}{r},\frac{z}{r}\right)=\frac{\mathbf{r}}{r}$，故

∇ u = -(q)/(4πε r²) r/r = -(q)/(4πε r³) r = -E.$$\nabla u=-\frac{q}{4\pi \epsilon r^2}\frac{\mathbf{r}}{r}=-\frac{q}{4\pi \epsilon r^3}\mathbf{r}=-\mathbf{E}.$$

因此，电位梯度的方向与电场强度 $\mathbf{E}$ 的方向相反，即沿径向的负方向电位 $u$ 增长最快。

例 3.15 鲨鱼追踪血腥味问题

一条鲨鱼在发现血腥味时，总是沿血腥味最浓的方向追寻。在海面上进行实验表明，如果把坐标原点取在血源处，在海平面上建立直角坐标系，那么点 $(x,y)$ 处血液的浓度 $C$（每百万分之一水中所含血的份数）的近似值为 $C=e^{-(x^2+2y^2)/10^4}$。求鲨鱼从点 $(x_0,y_0)$ 出发血源前进的路线。

解

设鲨鱼前进的路线为曲线 $\Gamma :y=f(x)$，我们首先来建立 $f(x)$ 应满足的方程。鲨鱼追踪最强的血腥味，所以每一瞬间它都将按血液浓度变化最快，即 $C$ 的梯度方向前进，由梯度的计算公式得

∇ C = ( (∂ C)/(∂ x), (∂ C)/(∂ y) ) = 10⁻4 e⁻(x² + 2y²)/10⁴ (-2x, -4y).$$\nabla C=\left(\frac{\partial C}{\partial x},\frac{\partial C}{\partial y}\right)=10^{-4}{e}^{-(x^2+2y^2)/10^4}(-2x,-4y).$$

取鲨鱼前进的方向为曲线 $\Gamma$ 的正向，相应方向的切线为正切线，正切线与 $x$ 轴正向的交角为 $\theta$（图 5.14），则在 $\Gamma$ 上点 $(x,y)$ 处 $\Gamma$ 的切线上的方向向量 $\mathbf{\tau }$ 可表示为 $(\cos \theta ,\sin \theta )$，或

(1, tan θ) = ( 1, (dy)/(dx) ),$$(1,\tan \theta )=\left(1,\frac{\mathrm{d}y}{\mathrm{d}x}\right),$$

从而也可表示为 $\mathbf{\tau }=(\mathrm{d}x,\mathrm{d}y)$。显然，$\mathbf{\tau }$ 必与 $\nabla C$ 平行同向，从而有

(dx)/(-2x) = (dy)/(-4y),$$\frac{\mathrm{d}x}{-2x}=\frac{\mathrm{d}y}{-4y},$$

于是可得 $y=f(x)$ 应满足微分方程

(dy)/(dx) = 2y/x. (3.24)$$\frac{\mathrm{d}y}{\mathrm{d}x}=2\frac{y}{x}.\text{\hspace{1em}(3.24)}$$

由于鲨鱼的初始位置为 $(x_0,y_0)$，即初值条件为 $y{|}_{x=x_0}=y_0$。方程 (3.24) 是一个一阶可分离变量的微分方程，容易求得其通解为

y = Ax²,$$y=A{x}^2,$$

其中 $A$ 为任意常数，代入初值条件 $y{|}_{x=x_0}=y_0$，得 $A=\frac{y_0}{x_0^2}$。于是鲨鱼从 $(x_0,y_0)$ 出发进攻的路线是

y = (y₀)/(x₀²) x².$$y=\frac{y_0}{x_0^2}{x}^2.$$

图 5.14 曲线切线方向与梯度方向的关系

3.4 高阶偏导数和高阶全微分

以二元函数为例，我们首先来介绍高阶偏导数的概念。若 $z=f(x,y)$ 在区域 $D$ 内的两个偏函数

zₓ = (∂ z)/(∂ x) = fₓ(x, y), z_y = (∂ z)/(∂ y) = f_y(x, y)$$z_x=\frac{\partial z}{\partial x}=f_x(x,y),z_y=\frac{\partial z}{\partial y}=f_y(x,y)$$

在 $D$ 内某点 $\mathbf{x}$ 处的偏导数仍存在，则这两个函数 $f_x(x,y)$ 与 $f_y(x,y)$ 的偏导数称为函数 $f(x,y)$ 的二阶偏导数。按照求导次序的不同，二元函数 $z$ 有以下四种不同的二阶偏导数，并分别记作

(∂)/(∂ x)((∂ z)/(∂ x)) = (∂² z)/(∂ x²) = zₓx = fₓx(x, y), (∂)/(∂ y)((∂ z)/(∂ x)) = (∂² z)/(∂ x ∂ y) = zₓy = fₓy(x, y),$$\frac{\partial }{\partial x}\left(\frac{\partial z}{\partial x}\right)=\frac{{\partial }^{2}z}{\partial x^2}=z_{xx}=f_{xx}(x,y),\frac{\partial }{\partial y}\left(\frac{\partial z}{\partial x}\right)=\frac{{\partial }^{2}z}{\partial x\partial y}=z_{xy}=f_{xy}(x,y),$$ (∂)/(∂ x)((∂ z)/(∂ y)) = (∂² z)/(∂ y ∂ x) = z_(yx) = f_(yx)(x, y), (∂)/(∂ y)((∂ z)/(∂ y)) = (∂² z)/(∂ y²) = z_(yy) = f_(yy)(x, y).$$\frac{\partial }{\partial x}\left(\frac{\partial z}{\partial y}\right)=\frac{{\partial }^{2}z}{\partial y\partial x}=z_{yx}=f_{yx}(x,y),\frac{\partial }{\partial y}\left(\frac{\partial z}{\partial y}\right)=\frac{{\partial }^{2}z}{\partial y^2}=z_{yy}=f_{yy}(x,y).$$

并称 $f_{xy}$ 和 $f_{yx}$ 为二阶混合偏导数。

类似地，可以定义 $n$ ($n\ge 3$) 阶偏导数。二阶及二阶以上的偏导数统称为高阶偏导数。

例 3.16 求二元函数的所有二阶偏导数

求二元函数 $z=x^y$ ($x>0$) 的所有二阶偏导数。

解

由 $\frac{\partial z}{\partial x}=y{x}^{y-1}$ 及 $\frac{\partial z}{\partial y}=x^y\ln x$，再分别关于变量 $x,y$ 求偏导数，得

(∂² z)/(∂ x²) = (∂)/(∂ x)((∂ z)/(∂ x)) = y(y-1)xʸ-2,$$\frac{{\partial }^{2}z}{\partial x^2}=\frac{\partial }{\partial x}\left(\frac{\partial z}{\partial x}\right)=y(y-1)x^{y-2},$$(∂² z)/(∂ x ∂ y) = (∂)/(∂ y)((∂ z)/(∂ x)) = xʸ-1 + yxʸ-1 ln x,$$\frac{{\partial }^{2}z}{\partial x\partial y}=\frac{\partial }{\partial y}\left(\frac{\partial z}{\partial x}\right)=x^{y-1}+y{x}^{y-1}\ln x,$$(∂² z)/(∂ y ∂ x) = (∂)/(∂ x)((∂ z)/(∂ y)) = yxʸ-1 ln x + xʸ-1,$$\frac{{\partial }^{2}z}{\partial y\partial x}=\frac{\partial }{\partial x}\left(\frac{\partial z}{\partial y}\right)=y{x}^{y-1}\ln x+x^{y-1},$$(∂² z)/(∂ y²) = (∂)/(∂ y)((∂ z)/(∂ y)) = xʸ (ln x)².$$\frac{{\partial }^{2}z}{\partial y^2}=\frac{\partial }{\partial y}\left(\frac{\partial z}{\partial y}\right)=x^y(\ln x{)}^2.$$

上例中混合偏导数 $z_{xy}=z_{yx}$，即在此例中，混合偏导数与求导的先后次序是无关的，但并非总是如此。

例 3.17 混合偏导数不相等的例子

设二元函数

f(x, y) = (xy(x² - y²))/(x² + y²), & x² + y² ≠ 0, 0, & x² + y² = 0,$$f(x,y)=\{\begin{array}{ll}\frac{xy(x^2-y^2)}{x^2+y^2},& x^2+y^2\ne 0,\\ 0,& x^2+y^2=0,\end{array}$$

证明：$f_{xy}(0,0)\ne f_{yx}(0,0)$。

证明

当 $x^2+y^2\ne 0$ 时利用求导法则，当 $x^2+y^2=0$ 时根据偏导数定义，可以求得

fₓ(x, y) = (y(x² - y²))/(x² + y²) + (xy · 4x² y²)/((x² + y²)²), & x² + y² ≠ 0, 0, & x² + y² = 0,$$f_x(x,y)=\{\begin{array}{ll}\frac{y(x^2-y^2)}{x^2+y^2}+\frac{xy\cdot 4x^2{y}^2}{(x^2+y^2{)}^2},& x^2+y^2\ne 0,\\ 0,& x^2+y^2=0,\end{array}$$f_y(x, y) = (x(x² - y²))/(x² + y²) - (xy · 4x² y²)/((x² + y²)²), & x² + y² ≠ 0, 0, & x² + y² = 0.$$f_y(x,y)=\{\begin{array}{ll}\frac{x(x^2-y^2)}{x^2+y^2}-\frac{xy\cdot 4x^2{y}^2}{(x^2+y^2{)}^2},& x^2+y^2\ne 0,\\ 0,& x^2+y^2=0.\end{array}$$

因此有 $f_x(0,y)=-y,f_y(x,0)=x$。再由偏导数定义可得

fₓy(0,0) = lim_(Δ y → 0) (fₓ(0, Δ y) - fₓ(0,0))/(Δ y) = lim_(Δ y → 0) (-Δ y)/(Δ y) = -1,$$f_{xy}(0,0)=\underset{\Delta y\to 0}{\lim }\frac{f_x(0,\Delta y)-f_x(0,0)}{\Delta y}=\underset{\Delta y\to 0}{\lim }\frac{-\Delta y}{\Delta y}=-1,$$f_(yx)(0,0) = lim_(Δ x → 0) (f_y(Δ x, 0) - f_y(0,0))/(Δ x) = lim_(Δ x → 0) (Δ x)/(Δ x) = 1.$$f_{yx}(0,0)=\underset{\Delta x\to 0}{\lim }\frac{f_y(\Delta x,0)-f_y(0,0)}{\Delta x}=\underset{\Delta x\to 0}{\lim }\frac{\Delta x}{\Delta x}=1.$$

所以 $f_{xy}(0,0)\ne f_{yx}(0,0)$。

我们看到，此例中的两个二阶混合偏导数不相等。一般地，$f_{xy}$ 与 $f_{yx}$ 是有区别的，前者是先对 $x$ 求导然后对 $y$ 求导，后者是先对 $y$ 求导再对 $x$ 求导，求导的次序是不同的。但是可以证明（证明略去）：当 $f_{xy}$ 和 $f_{yx}$ 都在点 $P$ 处连续时，则在点 $P$ 处有 $f_{xy}=f_{yx}$，即二阶混合偏导数与求导次序无关。

对于 $n$ 元函数的高阶偏导数，一般可以证明：如果所有的 $m$ 阶偏导数在点 $P$ 处连续，则在点 $P$ 处的 $m$ 阶偏导数与求导次序无关。例如，如果三元函数 $u=f(x,y,z)$ 在点 $P$ 的所有三阶偏导数连续，则在点 $P$ 有 $f_{xyy}=f_{yxy}=f_{yyx}$。

高阶全微分

我们知道，当二元函数 $u=f(x,y)$ 在区域 $\Omega \subseteq {\mathbf{R}}^2$ 内每一点均可微时，则 $\Omega$ 内 $u$ 的全分为

du = langle ∇ f, dx rangle = (∂ f)/(∂ x) dx + (∂ f)/(∂ y) dy.$$\mathrm{d}u=⟨\nabla f,\mathrm{d}\mathbf{x}⟩=\frac{\partial f}{\partial x}\mathrm{d}x+\frac{\partial f}{\partial y}\mathrm{d}y.$$

如果把 $\mathrm{d}x,\mathrm{d}y$ 看作固定不变，那么 $\mathrm{d}u$ 就是 $(x,y)$ 的函数。如果函数 $\mathrm{d}u$ 仍在 $\Omega$ 内可微，那么把这个函数 $\mathrm{d}u$ 再求全微分，其结果就称为 $u$ 的二阶全微分，记作 ${\mathrm{d}}^{2}u=\mathrm{d}(\mathrm{d}u)$。

设 $f$ 的各阶偏导数都存在且连续，对 $\mathrm{d}u=f_x\mathrm{d}x+f_y\mathrm{d}y$ 再求全微分可得

d² u = d(du) &= (∂)/(∂ x)(fₓ dx + f_y dy) dx + (∂)/(∂ y)(fₓ dx + f_y dy) dy &= fₓx dx² + 2fₓy dxdy + f_(yy) dy².$$\begin{array}{rl}{\mathrm{d}}^{2}u=\mathrm{d}(\mathrm{d}u)& =\frac{\partial }{\partial x}(f_x\mathrm{d}x+f_y\mathrm{d}y)\mathrm{d}x+\frac{\partial }{\partial y}(f_x\mathrm{d}x+f_y\mathrm{d}y)\mathrm{d}y\\ & =f_{xx}\mathrm{d}{x}^2+2f_{xy}\mathrm{d}x\mathrm{d}y+f_{yy}\mathrm{d}{y}^2.\end{array}$$

引入算符记号

(∂²)/(∂ x²) dx² + 2(∂²)/(∂ x ∂ y) dxdy + (∂²)/(∂ y²) dy² = ( (∂)/(∂ x) dx + (∂)/(∂ y) dy )²,$$\frac{{\partial }^2}{\partial x^2}\mathrm{d}{x}^2+2\frac{{\partial }^2}{\partial x\partial y}\mathrm{d}x\mathrm{d}y+\frac{{\partial }^2}{\partial y^2}\mathrm{d}{y}^2={\left(\frac{\partial }{\partial x}\mathrm{d}x+\frac{\partial }{\partial y}\mathrm{d}y\right)}^2,$$

上式左端相当于将右端按二项式展开，则二阶全微分也可简洁地写成

d² u = ( (∂)/(∂ x) dx + (∂)/(∂ y) dy )² f.$${\mathrm{d}}^{2}u={\left(\frac{\partial }{\partial x}\mathrm{d}x+\frac{\partial }{\partial y}\mathrm{d}y\right)}^{2}f.$$

类似地可以定义

d³ u = d(d² u), dⁿ u = d(dⁿ-1 u).$${\mathrm{d}}^{3}u=\mathrm{d}({\mathrm{d}}^{2}u),{\mathrm{d}}^{n}u=\mathrm{d}({\mathrm{d}}^{n-1}u).$$

而且可以验证

dⁿ u = ( (∂)/(∂ x) dx + (∂)/(∂ y) dy ) ( (∂)/(∂ x) dx + (∂)/(∂ y) dy )ⁿ-1 f = ( (∂)/(∂ x) dx + (∂)/(∂ y) dy )ⁿ f,$${\mathrm{d}}^{n}u=\left(\frac{\partial }{\partial x}\mathrm{d}x+\frac{\partial }{\partial y}\mathrm{d}y\right){\left(\frac{\partial }{\partial x}\mathrm{d}x+\frac{\partial }{\partial y}\mathrm{d}y\right)}^{n-1}f={\left(\frac{\partial }{\partial x}\mathrm{d}x+\frac{\partial }{\partial y}\mathrm{d}y\right)}^{n}f,$$

其中 ${\left(\frac{\partial }{\partial x}\mathrm{d}x+\frac{\partial }{\partial y}\mathrm{d}y\right)}^{n}f$ 表示将算符 $\left(\frac{\partial }{\partial x}\mathrm{d}x+\frac{\partial }{\partial y}\mathrm{d}y\right)$ 自乘 $n$ 次后再作用到 $f$ 上。

3.5 多元复合函数的偏导数和全微分

在一元函数的求导法中，复合函数的链式法则发挥了非常重要的作用。本段将把链式法则推广到多元函数。为了论述简洁，我们以由两个中间变量和两个自变量构成的复合函数 $z=f[u(x,y),v(x,y)]$ 为例来论述链式法则。

定理 3.5

若 $u=u(x,y)$ 和 $v=v(x,y)$ 均在点 $(x,y)$ 处可微，且函数 $z=f(u,v)$ 在对应的点 $(u,v)$ 处可微，则复合函数 $z=f[u(x,y),v(x,y)]$ 在点 $(x,y)$ 处也必可微，且其全微分为

dz = ( (∂ z)/(∂ u) (∂ u)/(∂ x) + (∂ z)/(∂ v) (∂ v)/(∂ x) ) dx + ( (∂ z)/(∂ u) (∂ u)/(∂ y) + (∂ z)/(∂ v) (∂ v)/(∂ y) ) dy. (3.25)$$\mathrm{d}z=\left(\frac{\partial z}{\partial u}\frac{\partial u}{\partial x}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial x}\right)\mathrm{d}x+\left(\frac{\partial z}{\partial u}\frac{\partial u}{\partial y}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial y}\right)\mathrm{d}y.\text{\hspace{1em}(3.25)}$$

注

在一元函数中，若函数 $u=\phi (x)$ 在 $x_0$ 处可导，而函数 $y=f(u)$ 在与 $x_0$ 对应的 $u_0$ 处可导，则复合函数 $y=f[\phi (x)]$ 在 $x_0$ 处必可导，且其导数为 ${\frac{\mathrm{d}y}{\mathrm{d}x}|}_{x_0}=f^{\prime }(u_0){\phi }^{\prime }(x_0)$。由于一元函数中可导与可微是等价的，所以定理3.5 可看作是一元复合函数求导法则在多元复合函数中的直接推广。

证

令自变量 $x,y$ 分别有改变量 $\Delta x,\Delta y$，相应地，函数 $u,v$ 分别有改变量 $\Delta u,\Delta v$，从而函数 $f$ 有改变量 $\Delta z$。由于 $u,v$ 均在点 $(x,y)$ 处可微，故有

Δ u = (∂ u)/(∂ x) Δ x + (∂ u)/(∂ y) Δ y + o₁(ρ), (3.26)$$\Delta u=\frac{\partial u}{\partial x}\Delta x+\frac{\partial u}{\partial y}\Delta y+o_1(\rho ),\text{\hspace{1em}(3.26)}$$Δ v = (∂ v)/(∂ x) Δ x + (∂ v)/(∂ y) Δ y + o₂(ρ), (3.27)$$\Delta v=\frac{\partial v}{\partial x}\Delta x+\frac{\partial v}{\partial y}\Delta y+o_2(\rho ),\text{\hspace{1em}(3.27)}$$

其中 $\rho =\sqrt{\Delta x^2+\Delta y^2}$，$o_i(\rho )(i=1,2)$ 是当 $\rho \to 0$ 时关于 $\rho$ 的高阶无穷小。又由于函数 $f$ 在 $(x,y)$ 所对应的 $(u,v)$ 处可微，故有

Δ z = (∂ z)/(∂ u) Δ u + (∂ z)/(∂ v) Δ v + o(√(Δ u² + Δ v²)). (3.28)$$\Delta z=\frac{\partial z}{\partial u}\Delta u+\frac{\partial z}{\partial v}\Delta v+o(\sqrt{\Delta u^2+\Delta v^2}).\text{\hspace{1em}(3.28)}$$

将 (3.26)、(3.27) 式代入 (3.28) 式并加以整理，则得复合函数 $z=f[u(x,y),v(x,y)]$ 的改变量为

Δ z = ( (∂ z)/(∂ u) (∂ u)/(∂ x) + (∂ z)/(∂ v) (∂ v)/(∂ x) ) Δ x + ( (∂ z)/(∂ u) (∂ u)/(∂ y) + (∂ z)/(∂ v) (∂ v)/(∂ y) ) Δ y + α, (3.29)$$\Delta z=\left(\frac{\partial z}{\partial u}\frac{\partial u}{\partial x}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial x}\right)\Delta x+\left(\frac{\partial z}{\partial u}\frac{\partial u}{\partial y}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial y}\right)\Delta y+\alpha ,\text{\hspace{1em}(3.29)}$$

其中

α = (∂ z)/(∂ u) o₁(ρ) + (∂ z)/(∂ v) o₂(ρ) + o(√(Δ u² + Δ v²)).$$\alpha =\frac{\partial z}{\partial u}{o}_1(\rho )+\frac{\partial z}{\partial v}{o}_2(\rho )+o(\sqrt{\Delta u^2+\Delta v^2}).$$

要证明 (3.25) 式，只需证明 (3.29) 式中的 $\alpha$ 为 $\rho$ 的高阶无穷小，即

lim_(ρ → 0) (α)/(ρ) = lim_(ρ → 0) [ (∂ z)/(∂ u) (o₁(ρ))/(ρ) + (∂ z)/(∂ v) (o₂(ρ))/(ρ) + (o(√(Δ u² + Δ v²)))/(ρ) ] = 0.$$\underset{\rho \to 0}{\lim }\frac{\alpha }{\rho }=\underset{\rho \to 0}{\lim }\left[\frac{\partial z}{\partial u}\frac{o_1(\rho )}{\rho }+\frac{\partial z}{\partial v}\frac{o_2(\rho )}{\rho }+\frac{o(\sqrt{\Delta u^2+\Delta v^2})}{\rho }\right]=0.$$

事实上，注意到 $\frac{\partial z}{\partial u},\frac{\partial z}{\partial v}$ 均与 $\rho$ 无关，以及 ${\lim }_{\rho \to 0}\frac{o_i(\rho )}{\rho }=0(i=1,2)$，从而有

lim_(ρ → 0) (1)/(ρ) [ (∂ z)/(∂ u) o₁(ρ) + (∂ z)/(∂ v) o₂(ρ) ] = 0,$$\underset{\rho \to 0}{\lim }\frac{1}{\rho }\left[\frac{\partial z}{\partial u}{o}_1(\rho )+\frac{\partial z}{\partial v}{o}_2(\rho )\right]=0,$$

因此，以下只需证明：

lim_(ρ → 0) (o(√(Δ u² + Δ v²)))/(ρ) = 0. (3.30)$$\underset{\rho \to 0}{\lim }\frac{o(\sqrt{\Delta u^2+\Delta v^2})}{\rho }=0.\text{\hspace{1em}(3.30)}$$

由于

(o(√(Δ u² + Δ v²)))/(ρ) = (o(√(Δ u² + Δ v²)))/(√(Δ u² + Δ v²)) · (√(Δ u² + Δ v²))/(ρ), (3.30)$$\frac{o(\sqrt{\Delta u^2+\Delta v^2})}{\rho }=\frac{o(\sqrt{\Delta u^2+\Delta v^2})}{\sqrt{\Delta u^2+\Delta v^2}}\cdot \frac{\sqrt{\Delta u^2+\Delta v^2}}{\rho },\text{\hspace{1em}(3.30)}$$

而当 $\rho$ 充分小时，由 (3.26) 式可见

| (Δ u)/(ρ) | ≤ | (∂ u)/(∂ x) | (|Δ x|)/(ρ) + | (∂ u)/(∂ y) | (|Δ y|)/(ρ) + (|o₁(ρ)|)/(ρ) < | (∂ u)/(∂ x) | + | (∂ u)/(∂ y) | + 1,$$\left|\frac{\Delta u}{\rho }\right|\le \left|\frac{\partial u}{\partial x}\right|\frac{|\Delta x|}{\rho }+\left|\frac{\partial u}{\partial y}\right|\frac{|\Delta y|}{\rho }+\frac{|o_1(\rho )|}{\rho }<\left|\frac{\partial u}{\partial x}\right|+\left|\frac{\partial u}{\partial y}\right|+1,$$

故 $\frac{\Delta u}{\rho }$ 有界，同理知 $\frac{\Delta v}{\rho }$ 也有界，于是知 $\frac{\sqrt{\Delta u^2+\Delta v^2}}{\rho }$ 有界。再由 $u,v$ 的可微性知 $u,v$ 在 $(x,y)$ 处连续，即当 $\rho \to 0$ 时 $\Delta u\to 0$ 及 $\Delta v\to 0$，所以

lim_(ρ → 0) (o(√(Δ u² + Δ v²)))/(√(Δ u² + Δ v²)) = 0,$$\underset{\rho \to 0}{\lim }\frac{o(\sqrt{\Delta u^2+\Delta v^2})}{\sqrt{\Delta u^2+\Delta v^2}}=0,$$

于是由 (3.30) 式知

lim_(ρ → 0) (o(√(Δ u² + Δ v²)))/(ρ) = 0.$$\underset{\rho \to 0}{\lim }\frac{o(\sqrt{\Delta u^2+\Delta v^2})}{\rho }=0.$$

当 $u=u(x,y),v=v(x,y)$ 与 $z=f(u,v)$ 均可微时，由 (3.25) 式可见，对复合函数 $z=f[u(x,y),v(x,y)]$ 有下列链式法则：

dfrac∂ z∂ x = dfrac∂ z∂ u dfrac∂ u∂ x + dfrac∂ z∂ v dfrac∂ v∂ x, dfrac∂ z∂ y = dfrac∂ z∂ u dfrac∂ u∂ y + dfrac∂ z∂ v dfrac∂ v∂ y, (3.31)$$\{\begin{array}{l}\frac{\partial z}{\partial x}=\frac{\partial z}{\partial u}\frac{\partial u}{\partial x}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial x},\\ \frac{\partial z}{\partial y}=\frac{\partial z}{\partial u}\frac{\partial u}{\partial y}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial y},\end{array}\text{\hspace{1em}(3.31)}$$

注意

公式 (3.31) 式的结构有如下特征：复合函数 $f$ 含有自变量 $x$ (或 $y$) 的中间变量有几个，则 $f$ 对自变量 $x$ (或 $y$) 的偏导数公式中就有几项之和，而且每一项的构成与一元函数的链式法则类似，即为“函数对中间变量的导数再乘以中间变量对自变量的导数”。

按照公式 (3.31)，我们可以将多元复合函数的求导法则推广到由 $m$ 个中间变量、$n$ 个自变量构成的一般复合函数中。设函数

y=f(u₁, u₂, …, uₘ) 及 uᵢ=uᵢ(x₁, x₂, …, xₙ), i=1, 2, …, m$$y=f(u_1,u_2,\cdots ,u_m)\text{及}{u}_i=u_i(x_1,x_2,\cdots ,x_n),i=1,2,\cdots ,m$$

都可微，则复合函数 $y=f(u_1(x),u_2(x),\cdots ,u_m(x))$ 也可微，其中 $\mathbf{x}=(x_1,x_2,\cdots ,x_n)$，且有

dy = (∂ y)/(∂ x₁) dx₁ + (∂ y)/(∂ x₂) dx₂ + … + (∂ y)/(∂ xₙ) dxₙ, (3.32)$$\mathrm{d}y=\frac{\partial y}{\partial x_1}\mathrm{d}{x}_1+\frac{\partial y}{\partial x_2}\mathrm{d}{x}_2+\cdots +\frac{\partial y}{\partial x_n}\mathrm{d}{x}_n,\text{\hspace{1em}(3.32)}$$

其中

(∂ y)/(∂ xⱼ) = (∂ y)/(∂ u₁) (∂ u₁)/(∂ xⱼ) + (∂ y)/(∂ u₂) (∂ u₂)/(∂ xⱼ) + … + (∂ y)/(∂ uₘ) (∂ uₘ)/(∂ xⱼ), j=1, 2, …, n. (3.33)$$\frac{\partial y}{\partial x_j}=\frac{\partial y}{\partial u_1}\frac{\partial u_1}{\partial x_j}+\frac{\partial y}{\partial u_2}\frac{\partial u_2}{\partial x_j}+\cdots +\frac{\partial y}{\partial u_m}\frac{\partial u_m}{\partial x_j},j=1,2,\cdots ,n.\text{\hspace{1em}(3.33)}$$

多元函数的复合可以有多种不同情况，读者应善于分析函数间的复合关系，弄清楚哪些是中间变量，哪些是自变量，把链式法则灵活地应用到各种场合。

复合函数求导实例

例 3.17 复合函数求导

(1) 设 $z=f(u,v),u=\phi (x),v=\psi (x)$ 均可微，则复合函数 $z=f[\phi (x),\psi (x)]$ 是 $x$ 的一元函数，应用公式 (3.32)，得

(dz)/(dx) = (∂ z)/(∂ u) (du)/(dx) + (∂ z)/(∂ v) (dv)/(dx).$$\frac{\mathrm{d}z}{\mathrm{d}x}=\frac{\partial z}{\partial u}\frac{\mathrm{d}u}{\mathrm{d}x}+\frac{\partial z}{\partial v}\frac{\mathrm{d}v}{\mathrm{d}x}.$$

它称为复合函数 $z$ 对 $x$ 的全导数。

(2) 设 $w=f(u),u=\phi (x,y,z)$ 均可微，则复合函数 $w=f[u(x,y,z)]$ 可微，它有一个中间变量，有三个自变量，应用公式 (3.33)，得

(∂ w)/(∂ x) = (dw)/(du) (∂ u)/(∂ x), (∂ w)/(∂ y) = (dw)/(du) (∂ u)/(∂ y), (∂ w)/(∂ z) = (dw)/(du) (∂ u)/(∂ z).$$\frac{\partial w}{\partial x}=\frac{\mathrm{d}w}{\mathrm{d}u}\frac{\partial u}{\partial x},\frac{\partial w}{\partial y}=\frac{\mathrm{d}w}{\mathrm{d}u}\frac{\partial u}{\partial y},\frac{\partial w}{\partial z}=\frac{\mathrm{d}w}{\mathrm{d}u}\frac{\partial u}{\partial z}.$$

(3) 设 $u=f(x,y,z),z=\phi (x,y)$ 均可微，则复合函数 $u=f[x,y,z(x,y)]$ 可微，它有三个中间变量，有两个自变量，应用公式 (3.33)，得

(∂ u)/(∂ x) = (∂ f)/(∂ x) + (∂ f)/(∂ z) (∂ z)/(∂ x), (∂ u)/(∂ y) = (∂ f)/(∂ y) + (∂ f)/(∂ z) (∂ z)/(∂ y).$$\frac{\partial u}{\partial x}=\frac{\partial f}{\partial x}+\frac{\partial f}{\partial z}\frac{\partial z}{\partial x},\frac{\partial u}{\partial y}=\frac{\partial f}{\partial y}+\frac{\partial f}{\partial z}\frac{\partial z}{\partial y}.$$

例 3.18 求偏导数

设 $z=f(x,xy)$，其中 $z=f(u,v)$ 可微，求 $\frac{\partial z}{\partial x},\frac{\partial z}{\partial y}$。

解

由于 $u=x$ 且 $v=xy$ 显然可微，故复合函数可微，应用公式 (3.31) 得

(∂ z)/(∂ x) = (∂ f)/(∂ x) + (∂ f)/(∂ v) (∂ v)/(∂ x) = (∂ f)/(∂ x) + (∂ f)/(∂ v) · y,$$\frac{\partial z}{\partial x}=\frac{\partial f}{\partial x}+\frac{\partial f}{\partial v}\frac{\partial v}{\partial x}=\frac{\partial f}{\partial x}+\frac{\partial f}{\partial v}\cdot y,$$(∂ z)/(∂ y) = (∂ f)/(∂ u) (∂ u)/(∂ y) + (∂ f)/(∂ v) (∂ v)/(∂ y) = (∂ f)/(∂ u) · 0 + (∂ f)/(∂ v) · x = x (∂ f)/(∂ v).$$\frac{\partial z}{\partial y}=\frac{\partial f}{\partial u}\frac{\partial u}{\partial y}+\frac{\partial f}{\partial v}\frac{\partial v}{\partial y}=\frac{\partial f}{\partial u}\cdot 0+\frac{\partial f}{\partial v}\cdot x=x\frac{\partial f}{\partial v}.$$

想一想

AI 解答

例3.18 的解中，左端的 $\frac{\partial z}{\partial x}$ 与右端的 $\frac{\partial f}{\partial x}$ 有什么不同？

把 $f(x,xy)$ 中的 $x$ 看作第一个变量，$xy$ 看作第二个变量，有时采用下面的记号来表示 $\frac{\partial z}{\partial x}$ 更为方便清晰：

(∂ z)/(∂ x) = f₁ + y f₂.$$\frac{\partial z}{\partial x}=f_1+y{f}_2.$$

其中 $f_1$ 表示 $f$ 对第一个变量的偏导数，$f_2$ 表示 $f$ 对第二个变量的偏导数。

例 3.19 证明恒等式

设 $u=\phi (x^2+y^2)$，其中 $\phi$ 可导，求证：$x\frac{\partial u}{\partial y}-y\frac{\partial u}{\partial x}=0$。

证

把 $u=\phi (x^2+y^2)$ 看作是由函数 $u=\phi (z)$ 及 $z=x^2+y^2$ 复合而成，分别对 $x$ 与 $y$ 求导得

(∂ u)/(∂ x) = φ'(z) · 2x, (∂ u)/(∂ y) = φ'(z) · 2y,$$\frac{\partial u}{\partial x}={\phi }^{\prime }(z)\cdot 2x,\frac{\partial u}{\partial y}={\phi }^{\prime }(z)\cdot 2y,$$

从而

x (∂ u)/(∂ y) - y (∂ u)/(∂ x) = x · 2y φ'(z) - y · 2x φ'(z) = 0.$$x\frac{\partial u}{\partial y}-y\frac{\partial u}{\partial x}=x\cdot 2y{\phi }^{\prime }(z)-y\cdot 2x{\phi }^{\prime }(z)=0.$$

例 3.20 求二阶偏导数

设 $z=f(u,x,y)$，其中 $f$ 具有对各变量的连续二阶偏导数，且 $u=x{e}^y$，求 $\frac{{\partial }^{2}z}{\partial x\partial y}$。

解

根据函数的复合结构（图 5.16）及复合函数的链式法则，得

(∂ z)/(∂ x) = (∂ f)/(∂ u) (∂ u)/(∂ x) + (∂ f)/(∂ x) = f₁ eʸ + f₂,$$\frac{\partial z}{\partial x}=\frac{\partial f}{\partial u}\frac{\partial u}{\partial x}+\frac{\partial f}{\partial x}=f_1{e}^y+f_2,$$

注意到 $f_1$ 和 $f_2$ 都是 $u,x,y$ 的三元函数，再由复合函数的链式法则得

(∂² z)/(∂ x ∂ y) = (∂)/(∂ y) ( (∂ z)/(∂ x) ) = (∂)/(∂ y) (f₁ eʸ + f₂) = (∂ f₁)/(∂ y) eʸ + f₁ eʸ + (∂ f₂)/(∂ y)$$\frac{{\partial }^{2}z}{\partial x\partial y}=\frac{\partial }{\partial y}\left(\frac{\partial z}{\partial x}\right)=\frac{\partial }{\partial y}(f_1{e}^y+f_2)=\frac{\partial f_1}{\partial y}{e}^y+f_1{e}^y+\frac{\partial f_2}{\partial y}$$= (f₁1 xeʸ + f₁3 eʸ) eʸ + f₁ eʸ + f₂1 xeʸ + f₂3,$$=(f_{11}x{e}^y+f_{13}{e}^y)e^y+f_1{e}^y+f_{21}x{e}^y+f_{23},$$

其中 $f_{ij}$ 表示先对第 $i$ 个变量求导，再对其第 $j$ 个变量求导的二阶偏导数。

图 5.16 复合函数结构示意图

坐标变换中的偏导数

例 3.21 极坐标变换下的偏导数

求 ${\left(\frac{\partial u}{\partial x}\right)}^2+{\left(\frac{\partial u}{\partial y}\right)}^2$ 与 $\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2}$ 在极坐标系中的表达式，其中 $u=F(x,y)$ 具有连续的二阶偏导数。

解

令 $x=\rho \cos \theta ,y=\rho \sin \theta$，从而

ρ = √(x²+y²), θ = arctan y/x (3.34)$$\rho =\sqrt{x^2+y^2},\theta =\arctan \frac{y}{x}\text{\hspace{1em}(3.34)}$$u = F(x, y) = F(ρ cos θ, ρ sin θ) stackreldef= F(ρ, θ) = overlineF(√(x²+y²), arctan y/x).$$u=F(x,y)=F(\rho \cos \theta ,\rho \sin \theta )\stackrel{\text{def}}{=}F(\rho ,\theta )=\overline{F}(\sqrt{x^2+y^2},\arctan \frac{y}{x}).$$

这样一来，就可以把 $u=F(x,y)$ 看作是由 $u=\overline{F}(\rho ,\theta )$ 与 $\rho =\sqrt{x^2+y^2},\theta =\arctan \frac{y}{x}$ 复合而成。显然，$\overline{F}(\rho ,\theta )$ 同样具有连续的二阶偏导数。应用链式法则得

(∂ u)/(∂ x) = (∂ u)/(∂ ρ) (∂ ρ)/(∂ x) + (∂ u)/(∂ θ) (∂ θ)/(∂ x), (∂ u)/(∂ y) = (∂ u)/(∂ ρ) (∂ ρ)/(∂ y) + (∂ u)/(∂ θ) (∂ θ)/(∂ y) (3.35)$$\frac{\partial u}{\partial x}=\frac{\partial u}{\partial \rho }\frac{\partial \rho }{\partial x}+\frac{\partial u}{\partial \theta }\frac{\partial \theta }{\partial x},\frac{\partial u}{\partial y}=\frac{\partial u}{\partial \rho }\frac{\partial \rho }{\partial y}+\frac{\partial u}{\partial \theta }\frac{\partial \theta }{\partial y}\text{\hspace{1em}(3.35)}$$

由 (3.34) 式有

(∂ ρ)/(∂ x) = (x)/(√(x²+y²)) = (x)/(ρ) = cos θ, (∂ ρ)/(∂ y) = (y)/(√(x²+y²)) = (y)/(ρ) = sin θ,$$\frac{\partial \rho }{\partial x}=\frac{x}{\sqrt{x^2+y^2}}=\frac{x}{\rho }=\cos \theta ,\frac{\partial \rho }{\partial y}=\frac{y}{\sqrt{x^2+y^2}}=\frac{y}{\rho }=\sin \theta ,$$(∂ θ)/(∂ x) = -(y)/(x²+y²) = -(sin θ)/(ρ), (∂ θ)/(∂ y) = (x)/(x²+y²) = (cos θ)/(ρ) (3.36)$$\frac{\partial \theta }{\partial x}=-\frac{y}{x^2+y^2}=-\frac{\sin \theta }{\rho },\frac{\partial \theta }{\partial y}=\frac{x}{x^2+y^2}=\frac{\cos \theta }{\rho }\text{\hspace{1em}(3.36)}$$

把它们代入 (3.35) 式得

(∂ u)/(∂ x) = (∂ u)/(∂ ρ) cos θ - (∂ u)/(∂ θ) (sin θ)/(ρ) (3.37)$$\frac{\partial u}{\partial x}=\frac{\partial u}{\partial \rho }\cos \theta -\frac{\partial u}{\partial \theta }\frac{\sin \theta }{\rho }\text{\hspace{1em}(3.37)}$$(∂ u)/(∂ y) = (∂ u)/(∂ ρ) sin θ + (∂ u)/(∂ θ) (cos θ)/(ρ) (3.38)$$\frac{\partial u}{\partial y}=\frac{\partial u}{\partial \rho }\sin \theta +\frac{\partial u}{\partial \theta }\frac{\cos \theta }{\rho }\text{\hspace{1em}(3.38)}$$

两式平方后相加得

( (∂ u)/(∂ x) )² + ( (∂ u)/(∂ y) )² = ( (∂ u)/(∂ ρ) )² + (1)/(ρ²) ( (∂ u)/(∂ θ) )².$${\left(\frac{\partial u}{\partial x}\right)}^2+{\left(\frac{\partial u}{\partial y}\right)}^2={\left(\frac{\partial u}{\partial \rho }\right)}^2+\frac{1}{{\rho }^2}{\left(\frac{\partial u}{\partial \theta }\right)}^2.$$

将 (3.37) 式两端再对 $x$ 求偏导数，应用链式法则，得

(∂² u)/(∂ x²) = (∂)/(∂ ρ) ( (∂ u)/(∂ ρ) cos θ - (∂ u)/(∂ θ) (sin θ)/(ρ) ) (∂ ρ)/(∂ x) + (∂)/(∂ θ) ( (∂ u)/(∂ ρ) cos θ - (∂ u)/(∂ θ) (sin θ)/(ρ) ) (∂ θ)/(∂ x)$$\frac{{\partial }^{2}u}{\partial x^2}=\frac{\partial }{\partial \rho }\left(\frac{\partial u}{\partial \rho }\cos \theta -\frac{\partial u}{\partial \theta }\frac{\sin \theta }{\rho }\right)\frac{\partial \rho }{\partial x}+\frac{\partial }{\partial \theta }\left(\frac{\partial u}{\partial \rho }\cos \theta -\frac{\partial u}{\partial \theta }\frac{\sin \theta }{\rho }\right)\frac{\partial \theta }{\partial x}$$

把 (3.36) 式中有关表达式代入并简化后得

(∂² u)/(∂ x²) = (∂² u)/(∂ ρ²) cos² θ - 2 (1)/(ρ) (∂² u)/(∂ ρ ∂ θ) sin θ cos θ + 2 (∂ u)/(∂ θ) (sin θ cos θ)/(ρ²) + (∂² u)/(∂ θ²) (sin² θ)/(ρ²) + (∂ u)/(∂ ρ) (sin² θ)/(ρ)$$\frac{{\partial }^{2}u}{\partial x^2}=\frac{{\partial }^{2}u}{\partial {\rho }^2}{\cos }^2\theta -2\frac{1}{\rho }\frac{{\partial }^{2}u}{\partial \rho \partial \theta }\sin \theta \cos \theta +2\frac{\partial u}{\partial \theta }\frac{\sin \theta \cos \theta }{{\rho }^2}+\frac{{\partial }^{2}u}{\partial {\theta }^2}\frac{{\sin }^2\theta }{{\rho }^2}+\frac{\partial u}{\partial \rho }\frac{{\sin }^2\theta }{\rho }$$

类似地，将 (3.38) 式两端再对 $y$ 求偏导并化简可得

(∂² u)/(∂ y²) = (∂² u)/(∂ ρ²) sin² θ + 2 (1)/(ρ) (∂² u)/(∂ ρ ∂ θ) sin θ cos θ - 2 (∂ u)/(∂ θ) (sin θ cos θ)/(ρ²) + (∂² u)/(∂ θ²) (cos² θ)/(ρ²) + (∂ u)/(∂ ρ) (cos² θ)/(ρ)$$\frac{{\partial }^{2}u}{\partial y^2}=\frac{{\partial }^{2}u}{\partial {\rho }^2}{\sin }^2\theta +2\frac{1}{\rho }\frac{{\partial }^{2}u}{\partial \rho \partial \theta }\sin \theta \cos \theta -2\frac{\partial u}{\partial \theta }\frac{\sin \theta \cos \theta }{{\rho }^2}+\frac{{\partial }^{2}u}{\partial {\theta }^2}\frac{{\cos }^2\theta }{{\rho }^2}+\frac{\partial u}{\partial \rho }\frac{{\cos }^2\theta }{\rho }$$

于是

(∂² u)/(∂ x²) + (∂² u)/(∂ y²) = (∂² u)/(∂ ρ²) + (1)/(ρ²) (∂² u)/(∂ θ²) + (1)/(ρ) (∂ u)/(∂ ρ).$$\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2}=\frac{{\partial }^{2}u}{\partial {\rho }^2}+\frac{1}{{\rho }^2}\frac{{\partial }^{2}u}{\partial {\theta }^2}+\frac{1}{\rho }\frac{\partial u}{\partial \rho }.$$

数字资源 5.3.4怎样求多元复合函数二阶偏导数 ↗

一阶全微分形式的不变性

在一元函数中，一阶微分具有形式不变性。下面以二元函数为例证明多元函数的一阶全微分也具有形式不变性。

设有函数 $z=f(u,v)$ 与 $u=\phi (x,y),v=\psi (x,y)$ 复合，若 $\phi ,\psi$ 在 $(x,y)$ 处可微，且 $f$ 在与 $(x,y)$ 相应的 $(u,v)$ 处可微，则复合函数的全微分为

dz = (∂ z)/(∂ x) dx + (∂ z)/(∂ y) dy.$$\mathrm{d}z=\frac{\partial z}{\partial x}\mathrm{d}x+\frac{\partial z}{\partial y}\mathrm{d}y.$$

由复合函数求导法则，

dfrac∂ z∂ x = dfrac∂ z∂ u dfrac∂ u∂ x + dfrac∂ z∂ v dfrac∂ v∂ x, dfrac∂ z∂ y = dfrac∂ z∂ u dfrac∂ u∂ y + dfrac∂ z∂ v dfrac∂ v∂ y,$$\{\begin{array}{l}\frac{\partial z}{\partial x}=\frac{\partial z}{\partial u}\frac{\partial u}{\partial x}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial x},\\ \frac{\partial z}{\partial y}=\frac{\partial z}{\partial u}\frac{\partial u}{\partial y}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial y},\end{array}$$

代入上式并整理得

dz &= ( (∂ z)/(∂ u) (∂ u)/(∂ x) + (∂ z)/(∂ v) (∂ v)/(∂ x) ) dx + ( (∂ z)/(∂ u) (∂ u)/(∂ y) + (∂ z)/(∂ v) (∂ v)/(∂ y) ) dy &= (∂ z)/(∂ u) ( (∂ u)/(∂ x) dx + (∂ u)/(∂ y) dy ) + (∂ z)/(∂ v) ( (∂ v)/(∂ x) dx + (∂ v)/(∂ y) dy ).$$\begin{array}{rl}\mathrm{d}z& =\left(\frac{\partial z}{\partial u}\frac{\partial u}{\partial x}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial x}\right)\mathrm{d}x+\left(\frac{\partial z}{\partial u}\frac{\partial u}{\partial y}+\frac{\partial z}{\partial v}\frac{\partial v}{\partial y}\right)\mathrm{d}y\\ & =\frac{\partial z}{\partial u}\left(\frac{\partial u}{\partial x}\mathrm{d}x+\frac{\partial u}{\partial y}\mathrm{d}y\right)+\frac{\partial z}{\partial v}\left(\frac{\partial v}{\partial x}\mathrm{d}x+\frac{\partial v}{\partial y}\mathrm{d}y\right).\end{array}$$

而 $\frac{\partial u}{\partial x}\mathrm{d}x+\frac{\partial u}{\partial y}\mathrm{d}y=\mathrm{d}u,\frac{\partial v}{\partial x}\mathrm{d}x+\frac{\partial v}{\partial y}\mathrm{d}y=\mathrm{d}v$，所以

dz = (∂ z)/(∂ u) du + (∂ z)/(∂ v) dv. (3.39)$$\mathrm{d}z=\frac{\partial z}{\partial u}\mathrm{d}u+\frac{\partial z}{\partial v}\mathrm{d}v.\text{\hspace{1em}(3.39)}$$

上面的全微分形式与把函数 $z=f(u,v)$ 中的中间变量 $u,v$ 看作自变量时的全微分形式是完全一样的，我们把这一性质称为一阶全微分形式不变性。换句话说，对于二元函数 $z=f(u,v)$ 来说，无论 $u,v$ 是自变量还是中间变量，它的一阶全微分（若存在的话）均可写成 (3.39) 式形式。

想一想

AI 解答

试由一阶全微分形式不变性证明全微分的运算法则：$\mathrm{d}(u\pm v)=\mathrm{d}u\pm \mathrm{d}v$。

由一阶全微分形式不变性，容易得到全微分的有理运算法则： (1) $\mathrm{d}(u\pm v)=\mathrm{d}u\pm \mathrm{d}v$; (2) $\mathrm{d}(uv)=v\mathrm{d}u+u\mathrm{d}v$; (3) $\mathrm{d}\left(\frac{u}{v}\right)=\frac{1}{v^2}(v\mathrm{d}u-u\mathrm{d}v),v\ne 0$.

例 3.22 利用不变性求偏导数

设 $f(u,v)$ 可微，求 $z=f\left(\frac{x}{y},\frac{y}{x}\right)$ 的偏导数。

解

利用一阶全微分形式不变性：

dz = f₁ d(x/y) + f₂ d(y/x) = f₁ (y dx - x dy)/(y²) + f₂ (x dy - y dx)/(x²)$$\mathrm{d}z=f_1\mathrm{d}\left(\frac{x}{y}\right)+f_2\mathrm{d}\left(\frac{y}{x}\right)=f_1\frac{y\mathrm{d}x-x\mathrm{d}y}{y^2}+f_2\frac{x\mathrm{d}y-y\mathrm{d}x}{x^2}$$= ( 1/y f₁ - (y)/(x²) f₂ ) dx + ( -(x)/(y²) f₁ + 1/x f₂ ) dy,$$=\left(\frac{1}{y}{f}_1-\frac{y}{x^2}{f}_2\right)\mathrm{d}x+\left(-\frac{x}{y^2}{f}_1+\frac{1}{x}{f}_2\right)\mathrm{d}y,$$

从而

(∂ z)/(∂ x) = 1/y f₁ - (y)/(x²) f₂, (∂ z)/(∂ y) = -(x)/(y²) f₁ + 1/x f₂.$$\frac{\partial z}{\partial x}=\frac{1}{y}{f}_1-\frac{y}{x^2}{f}_2,\frac{\partial z}{\partial y}=-\frac{x}{y^2}{f}_1+\frac{1}{x}{f}_2.$$

应当指出，高阶全微分不具有微分形式不变性。

3.6 由一个方程确定的隐函数的微分法

我们常常会碰到一些函数，其因变量与自变量的关系是以方程形式联系起来的。例如在球面方程

x² + y² + z² = 1$$x^2+y^2+z^2=1$$

中，如果把 $x$ 和 $y$ 看作是自变量，那么此球面方程在平面闭区域 $D=\{(x,y)\mid x^2+y^2\le 1\}$ 上确定了两个连续的二元函数

z = ± √(1 - (x² + y²)).$$z=\pm \sqrt{1-(x^2+y^2)}.$$

一般地，设有方程

F(x₁, …, xₙ, y) = 0, (3.40)$$F(x_1,\cdots ,x_n,y)=0,\text{\hspace{1em}(3.40)}$$

如果存在一个 $n$ 元函数 $y=\phi (\mathbf{x}),\mathbf{x}=(x_1,\cdots ,x_n)\in \Omega \subseteq {\mathbf{R}}^n$ ($\Omega$ 为一区域) 使得将 $y=\phi (\mathbf{x})$ 代入 (3.40) 后成为恒等式

F(x₁, …, xₙ, φ(x₁, …, xₙ)) ≡ 0,$$F(x_1,\cdots ,x_n,\phi (x_1,\cdots ,x_n))\equiv 0,$$

则称 $y=\phi (\mathbf{x})$ 是由方程 (3.40) 所确定的隐函数。

我们在上册中已经讲过如何不经过显化而直接由方程 $F(x,y)=0$ 求出由它所确定的隐函数的导数的方法。现在利用偏导数给出一个隐函数存在的充分条件和隐函数的导数公式。

定理 3.6 (隐函数存在定理)

如果二元函数 $F(x,y)$ 满足： (1) $F(x_0,y_0)=0$; (2) 在点 $(x_0,y_0)$ 的某邻域中有连续的偏导数; (3) $F_y(x_0,y_0)\ne 0$,

则方程 $F(x,y)=0$ 在点 $(x_0,y_0)$ 的某邻域中唯一确定了一个具有连续导数的函数 $y=f(x)$，它满足 $y_0=f(x_0)$ 及 $F(x,f(x))\equiv 0$，并且

(dy)/(dx) = -(Fₓ)/(F_y). (3.41)$$\frac{\mathrm{d}y}{\mathrm{d}x}=-\frac{F_x}{F_y}.\text{\hspace{1em}(3.41)}$$

这个定理的证明从略，下面仅在方程 $F(x,y)=0$ 已经确定了具有连续导数的函数 $y=f(x)$ 的假定下，来推出公式 (3.41)。将 $y=f(x)$ 代入 $F(x,y)=0$ 得

F(x, f(x)) ≡ 0,$$F(x,f(x))\equiv 0,$$

由链式法则，上式两端对 $x$ 求导得

Fₓ + F_y · (dy)/(dx) = 0.$$F_x+F_y\cdot \frac{\mathrm{d}y}{\mathrm{d}x}=0.$$

由于 $F$ 连续且 $F_y(x_0,y_0)\ne 0$，所以存在点 $(x_0,y_0)$ 的一个邻域，在这个邻域中 $F_y\ne 0$，于是

(dy)/(dx) = -(Fₓ)/(F_y).$$\frac{\mathrm{d}y}{\mathrm{d}x}=-\frac{F_x}{F_y}.$$

想一想

AI 解答

怎样由 $F_y$ 连续且 $F_y(x_0,y_0)\ne 0$ 推出存在 $(x_0,y_0)$ 的一个邻域，在此邻域中 $F_y\ne 0$?

如果 $F(x,y)$ 的二阶偏导数连续，则由链式法则，对 (3.41) 的两端关于 $x$ 再求导，就得到了二阶导数 $\frac{{\mathrm{d}}^{2}y}{\mathrm{d}{x}^2}$。

隐函数的求导方法可以推广到多元函数。例如，若一个三元方程

F(x, y, z) = 0$$F(x,y,z)=0$$

确定了一个二元函数 $z=f(x,y)$，则将 $z=f(x,y)$ 代入 (3.42)，得

F(x, y, f(x, y)) ≡ 0,$$F(x,y,f(x,y))\equiv 0,$$

应用链式法则，将上式两端分别对 $x$ 和 $y$ 求导，得

Fₓ + F_z dfrac∂ z∂ x = 0, F_y + F_z dfrac∂ z∂ y = 0,$$\{\begin{array}{l}{F}_x+F_z\frac{\partial z}{\partial x}=0,\\ F_y+F_z\frac{\partial z}{\partial y}=0,\end{array}$$

从而在 $F_z\ne 0$ 处有

dfrac∂ z∂ x = -dfracFₓF_z, dfrac∂ z∂ y = -dfracF_yF_z. (3.43)$$\{\begin{array}{l}\frac{\partial z}{\partial x}=-\frac{F_x}{F_z},\\ \frac{\partial z}{\partial y}=-\frac{F_y}{F_z}.\end{array}\text{\hspace{1em}(3.43)}$$

例 3.23 隐函数偏导数

设 $\phi (u,v)$ 具有连续的一阶偏导数，方程 $\phi (cx-az,cy-bz)=0$ 确定了函数 $z=z(x,y)$，求 $a{z}_x+b{z}_y$，其中 $a,b,c$ 均为常数。

解法一

令 $F(x,y,z)=\phi (cx-az,cy-bz)$，显然复合函数 $F(x,y,z)$ 有对 $x,y,z$ 的连续一阶偏导数。由函数求导公式 (3.43) 得

zₓ = -(Fₓ)/(F_z) = -(cφ₁)/(-aφ₁ - bφ₂) = (cφ₁)/(aφ₁ + bφ₂),$$z_x=-\frac{F_x}{F_z}=-\frac{c{\phi }_1}{-a{\phi }_1-b{\phi }_2}=\frac{c{\phi }_1}{a{\phi }_1+b{\phi }_2},$$z_y = -(F_y)/(F_z) = -(cφ₂)/(-aφ₁ - bφ₂) = (cφ₂)/(aφ₁ + bφ₂),$$z_y=-\frac{F_y}{F_z}=-\frac{c{\phi }_2}{-a{\phi }_1-b{\phi }_2}=\frac{c{\phi }_2}{a{\phi }_1+b{\phi }_2},$$

故

azₓ + bz_y = c.$$a{z}_x+b{z}_y=c.$$

解法二

利用一阶全微分形式不变性。由给定方程两端求全微分得

φ₁(cdx - adz) + φ₂(cdy - bdz) = 0,$${\phi }_1(c\mathrm{d}x-a\mathrm{d}z)+{\phi }_2(c\mathrm{d}y-b\mathrm{d}z)=0,$$

由此解得

dz = (cφ₁ dx + cφ₂ dy)/(aφ₁ + bφ₂),$$\mathrm{d}z=\frac{c{\phi }_1\mathrm{d}x+c{\phi }_2\mathrm{d}y}{a{\phi }_1+b{\phi }_2},$$

所以

zₓ = (cφ₁)/(aφ₁ + bφ₂), z_y = (cφ₂)/(aφ₁ + bφ₂),$$z_x=\frac{c{\phi }_1}{a{\phi }_1+b{\phi }_2},z_y=\frac{c{\phi }_2}{a{\phi }_1+b{\phi }_2},$$

从而得

azₓ + bz_y = c.$$a{z}_x+b{z}_y=c.$$

例 3.24 隐函数全微分

设方程 $xyz+\sqrt{x^2+y^2+z^2}=\sqrt{2}$ 确定了函数 $z=z(x,y)$，点 $(1,0,-1)$ 处的全微分 $\mathrm{d}z$。

解法一

利用一阶全微分形式不变性。由给定方程两端求全微分得

yzdx + xzdy + xydz + (xdx + ydy + zdz)/(√(x²+y²+z²)) = 0.$$yz\mathrm{d}x+xz\mathrm{d}y+xy\mathrm{d}z+\frac{x\mathrm{d}x+y\mathrm{d}y+z\mathrm{d}z}{\sqrt{x^2+y^2+z^2}}=0.$$

将 $x=1,y=0,z=-1$ 代入上式得

-dy + (1)/(√(2))(dx - dz) = 0.$$-\mathrm{d}y+\frac{1}{\sqrt{2}}(\mathrm{d}x-\mathrm{d}z)=0.$$

从而

dz|_((1,0,-1)) = dx - √(2)dy.$$\mathrm{d}z{|}_{(1,0,-1)}=\mathrm{d}x-\sqrt{2}\mathrm{d}y.$$

解法二

利用求导公式 (3.43) 得

zₓ = -( fracyz + (x)/(√(x²+y²+z²))xy + (z)/(√(x²+y²+z²)) ), z_y = -( fracxz + (y)/(√(x²+y²+z²))xy + (z)/(√(x²+y²+z²)) ).$$z_x=-\left(\frac{yz+\frac{x}{\sqrt{x^2+y^2+z^2}}}{xy+\frac{z}{\sqrt{x^2+y^2+z^2}}}\right),z_y=-\left(\frac{xz+\frac{y}{\sqrt{x^2+y^2+z^2}}}{xy+\frac{z}{\sqrt{x^2+y^2+z^2}}}\right).$$

在点 $(1,0,-1)$ 处，$z_x=1,z_y=-\sqrt{2}$，从而 $\mathrm{d}z=\mathrm{d}x-\sqrt{2}\mathrm{d}y$。

从以上几例可以看出，利用一阶全微分形式不变性求隐函数的全微分或偏导数有其显著的优点，它利用求导公式要方便，特别是在变量间的关系较复杂时，用此方法无需去弄清各变量间的关系，只是将各变量一律视为自变量去求全微分，这样既简化了问题，也不易出错。

例 3.25 求隐函数的二阶偏导数

求 $e^z+xyz=e$ 所确定函数 $z=z(x,y)$ 的一阶与二阶偏导数 $\frac{\partial z}{\partial x},\frac{\partial z}{\partial y},\frac{{\partial }^{2}z}{\partial x\partial y}$。

解

令 $F(x,y,z)=e^z+xyz-e$，则由公式 (3.43) 得

(∂ z)/(∂ x) = -(Fₓ)/(F_z) = -(yz)/(e^z + xy), (∂ z)/(∂ y) = -(F_y)/(F_z) = -(xz)/(e^z + xy).$$\frac{\partial z}{\partial x}=-\frac{F_x}{F_z}=-\frac{yz}{e^z+xy},\frac{\partial z}{\partial y}=-\frac{F_y}{F_z}=-\frac{xz}{e^z+xy}.$$

从而有

(∂² z)/(∂ x ∂ y) = (∂)/(∂ y) ( -(yz)/(e^z + xy) ) = -(( z + y (∂ z)/(∂ y) )(e^z + xy) - yz ( e^z (∂ z)/(∂ y) + x ))/((e^z + xy)²).$$\frac{{\partial }^{2}z}{\partial x\partial y}=\frac{\partial }{\partial y}\left(-\frac{yz}{e^z+xy}\right)=-\frac{\left(z+y\frac{\partial z}{\partial y}\right)(e^z+xy)-yz\left(e^z\frac{\partial z}{\partial y}+x\right)}{(e^z+xy{)}^2}.$$

将前面已经求得的 $\frac{\partial z}{\partial y}$ 代入上式得

(∂² z)/(∂ x ∂ y) = -(( z - (xyz)/(e^z + xy) )(e^z + xy) - yz ( -(xze^z)/(e^z + xy) + x ))/((e^z + xy)²) = -(ze² + xyz² e^z - x² y² z)/((e^z + xy)³).$$\frac{{\partial }^{2}z}{\partial x\partial y}=-\frac{\left(z-\frac{xyz}{e^z+xy}\right)(e^z+xy)-yz\left(-\frac{xz{e}^z}{e^z+xy}+x\right)}{(e^z+xy{)}^2}=-\frac{z{e}^2+xy{z}^2{e}^z-x^2{y}^{2}z}{(e^z+xy{)}^3}.$$

5.3 多元数量值函数的导数与微分