梯度的应用

A. 误差

首先,我们介绍如何使用全微分来估计由已知误差的量计算出的某个量的误差。例如,质量 m 1 , m 2 在距离 r 时的万有引力由下式给出

m 1 m 2 r 2 = f ( m 1 , m 2 , r )

我们可能会问:如果我们测得 m 1 m 2 各为 1000  grams ± .01  g. ,且 r = 10  cm. ± .001  cm. ,那么用这些数值计算力时会产生多大的误差?由于我们只关心一个估计值,我们忽略 d f 表达式中的误差项,得到

d m 1 = + .01 = d m 2 , d r = .001 时,其绝对值最大。于是我们有

d f = 10 ( .01 ) + 10 ( .01 ) + 2 10 3 ( .001 ) = 2.2

因此我们对力的最大误差的估计为 ± 2.2 达因。

B. 极大值与极小值

我们可以对多元函数的相对极大值和极小值做一个简单的检验。我们将以二元函数的形式来讨论,但推广到更多变量的函数是显而易见的。点 ( x 0 , y 0 ) 称为 f ( x , y ) 相对极大值点,如果存在 ( x 0 , y 0 ) 的一个邻域,使得对于该邻域内的所有点 ( x , y )

f ( x , y ) f ( x 0 , y 0 )

相对极小值可类似定义。假设 f x f y 在相对极大值点 ( x 0 , y 0 ) 处存在。那么函数 g ( x ) = f ( x , y 0 ) x 0 处有相对极大值,所以 。同样地, h ( y ) = f ( x 0 , y ) y 0 处有相对极大值,所以 。因此我们得出结论:

定理 1

如果 f ( x , y ) ( x 0 , y 0 ) 处有相对极大值或极小值,那么

f x ( x 0 , y 0 ) = 0 = f y ( x 0 , y 0 )

其逆命题并不比一元函数的情形更成立。和以前一样,也有高阶导数检验法,但我们不在此讨论。

C. 曲面的切平面与法线

接下来考虑满足 f ( x , y , z ) = c 的点集 ( x , y , z ) ,其中 c 是某个常数,且 f ( x , y , z ) 具有连续的一阶偏导数。这样的点集称为函数 f ( x , y , z ) 的一个等值面。它完全类似于我们绘制等高线图时使用的二元函数 f ( x , y ) 的等值线。等值面的一个特例是二元函数 g ( x , y ) 的“图像”。在这种情况下我们绘制的点是 ( x , y , g ( x , y ) ) ,而这些点恰好是等值面 f ( x , y , z ) = 0 上的点,其中

f ( x , y , z ) = z g ( x , y )

一个不属于这种类型的例子是以原点为中心、半径为 5 的球面,它是等值面

x 2 + y 2 + z 2 = 25

我们希望求出这样一个曲面在曲面上一点 ( x 0 , y 0 , z 0 ) 处的法线切平面。如果我们能在某种意义下找到一个在 ( x 0 , y 0 , z 0 ) 处垂直于曲面的向量,那么问题就解决了。因为这个向量 A 将平行于法线且垂直于切平面。我们要证明的是,在 ( x 0 , y 0 , z 0 ) 处取值的 grad f 正是这样的向量 A

但是,一个向量在 ( x 0 , y 0 , z 0 ) 处垂直于曲面是什么意思呢?如果我们知道曲面上存在经过 ( x 0 , y 0 , z 0 ) 直线,我们可以说该向量应当垂直于所有这样的直线。但例如球面上并没有直线;所以接下来我们尝试找到一个在 ( x 0 , y 0 , z 0 ) 处垂直于曲面上所有经过 ( x 0 , y 0 , z 0 ) 曲线的向量。因此,设 X ( t ) = ( x ( t ) , y ( t ) , z ( t ) ) 是曲面上的一条曲线,并设 ( x 0 , y 0 , z 0 ) = X ( t 0 ) 。由于曲线位于曲面上,我们知道

f ( x ( t ) , y ( t ) , z ( t ) ) = c for all  t

因此

d d t ( f ( x ( t ) , y ( t ) , z ( t ) ) ) = d c d t = 0

但我们已经看到

所以

特别地,

是与曲线 X ( t ) t = t 0 处的切线平行的向量。由于 ( grad  f ) ( x 0 , y 0 , z 0 ) 垂直于 ,它垂直于曲线在 t 0 处的切线,并且这对曲面上经过 ( x 0 , y 0 , z 0 ) 的任何曲线都成立。因此 ( grad  f ) ( x 0 , y 0 , z 0 ) 就是所求的向量 A 。这为 grad f 的方向给出了一个几何解释:grad f 在 ( x 0 , y 0 , z 0 ) 处取值时,垂直于等值面 f ( x , y , z ) = f ( x 0 , y 0 , z 0 ) 在点 ( x 0 , y 0 , z 0 ) 。我们将在下一段中给出其长度的解释。

D. 方向导数

f ( x , y , z ) 是一个三元函数(尽管 n 元函数也同样适用)。设 A = ( a 1 , a 2 , a 3 ) 是一个给定的非零向量,并设 X 0 = ( x 0 , y 0 , z 0 ) 是一个给定的点。当我们在 X 0 处求 f x , f y , f z 时,我们求的是 f 关于从点 X 0 出发沿三个坐标轴方向之一所经过距离的变化率。现在假设改为将向量 A 以点 X 0 为起点放置并沿 A 行进,换句话说,我们从 X 0 沿 A 的方向行进,并且我们希望度量 f 关于所行进距离的变化率。如果我们从 X 0 走到 X 0 + d X ,那么“平均变化率”

f ( x 0 + d x , y 0 + d y , z 0 + d z ) f ( x 0 , y 0 , z 0 ) | d X |

现在向量 d X = ( X 0 + d X ) X 0 应与 A 具有相同的方向,所以 d X = h A ,其中某个 h > 0 。因此 | d X | = h | A | 。此外,

f ( x 0 + d x , y 0 + d y , z 0 + d z ) f ( x 0 , y 0 , z 0 ) = ( grad  f ) d X + error ,

其中

error | d X | 0 as  | d X | 0

“瞬时变化率”,我们称之为 f ( x , y , z ) X 0 处沿 A 方向的方向导数,将是平均变化率当 | d X | 0 时的极限,或者由于 | d X | 0 h 0 含义相同,即为

lim h 0 ( grad  f ) h A h | A | + error h | A |

但我们已经看到第二项趋于零,而第一项就是

( grad  f ) A | A |

这个结果为我们提供了一个求方向导数的简单公式。

现在注意到,如果 θ A grad f 之间的夹角,那么 ( grad  f ) A = | grad  f | | A | cos θ ,即方向导数就是 | grad  f | cos θ 。这当然在 cos θ = 1 时最大,即当 A grad f 方向相同时,此时它就是 | grad  f | 。因此我们得到了对 | grad  f | 的解释:它是 f 在任意方向上的方向导数的最大值。还要注意, f 的方向导数的最小值出现在方向与 grad f 相反时,且该最小值为 | grad  f |

习题

习题 1.

研究下列函数的相对极大值和极小值:
a) f ( x , y , z ) = 3 x 2 y 2 z 2
b) f ( x , y ) = 1 2 ( x 2 + y 2 )
c) f ( x , y ) = 1 2 ( x 2 y 2 )

习题 2.

f ( x , y , z ) = z e x sin y , P = ( log 3 , 3 π 2 , 3 )
求:
a) fP 处的 grad f
b) 经过 P f 的等值面在 P 处的法线;
c) 该曲面在 P 处的切平面;
d) f P 处沿向量 ( 1 , 2 , 2 ) 方向的方向导数。
e) f P 处方向导数的最大值和最小值。

习题 3.

考虑动能公式: K = 1 2 m v 2 。假设测得 m 20  grams ,最大误差为 ± .01  gram ,且测得 v 100  cm/sec ,最大误差为 ± 2  cm/sec 。估计计算动能时产生的最大误差,以及 K 百分比误差。

习题 4.

证明 f ( x , y , z ) 是常数当且仅当 grad f = 0