文章目录
- 6 全微分
- 梯度的概念
- 全微分
前情提要 R语言微积分
- 极限
- π , e , γ \pi, e, \gamma π,e,γ
- 洛必达法则
- 连续性和导数
- 数值导数
- 差商与牛顿插值
- 方向导数
6 全微分
梯度的概念
对于任意函数 f ( x 0 , x 1 , ⋯ , x n ) f(x_0,x_1,\cdots,x_n) f(x0,x1,⋯,xn),其梯度为
∇ f = ( ∂ f ∂ x 0 , ∂ f ∂ x 0 , ⋯ , ∂ f ∂ x 0 ) \nabla f=(\frac{\partial f}{\partial x_0}, \frac{\partial f}{\partial x_0},\cdots, \frac{\partial f}{\partial x_0}) ∇f=(∂x0∂f,∂x0∂f,⋯,∂x0∂f)
由于在求导过程中,常数项被置零,所以 ∇ f \nabla f ∇f必然反映出 f = 0 f=0 f=0的某些几何特征。
z = 1 − x 2 − y 2 z=1-x^2-y^2 z=1−x2−y2对应的方程是 x 2 + y 2 = 0 x^2+y^2=0 x2+y2=0,是一个圆形,其梯度为 ( − 2 x , − 2 y ) (-2x,-2y) (−2x,−2y),做图如下
theta = seq(-pi,pi,0.1)
x = cos(theta)
y = sin(theta)
plot(x,y,type='l',col='red')
x1 = x*0
y1 = (x1-x)/(-2*x)*(-2*y)+y
for(i in 1:length(theta))
lines(c(x[i],x1[i]),c(y[i],y1[i]),col='green')
可见,梯度方向对应的是图形的法线方向。对于二维平面内的曲线而言,其法线方向与切线方向垂直。
关于切线方向,我们仅知道 y = f ( x ) y=f(x) y=f(x)的切线的斜率为 k = f ′ ( x ) k=f'(x) k=f′(x),对应的方向为 ( 1 , f ′ ( x ) ) (1,f'(x)) (1,f′(x))。或者更直观地写为 d y d x \frac{\text dy}{\text dx} dxdy。
而今却只知道 x , y x,y x,y的隐函数表达式 f ( x , y ) = 0 f(x,y)=0 f(x,y)=0,对此可以先写出 x , y x,y x,y的偏导数 f x ′ = ∂ f ∂ x , f y ′ = ∂ f ∂ y f'_x=\frac{\partial f}{\partial x},f'_y=\frac{\partial f}{\partial y} fx′=∂x∂f,fy′=∂y∂f。
全微分
回顾偏导数的定义
f x ′ ( x 0 , y 0 ) = lim Δ x → 0 f ( x 0 + Δ x , y 0 ) − f ( x 0 , y 0 ) Δ x f'_x(x_0,y_0)=\lim_{\Delta x\to0}\frac{f(x_0+\Delta x,y_0)-f(x_0,y_0)}{\Delta x} fx′(x0,y0)=Δx→0limΔxf(x0+Δx,y0)−f(x0,y0)
这里出现一个很吊诡的问题,既然满足 f ( x , y ) = 0 f(x,y)=0 f(x,y)=0,那么 f ( x , y ) = 0 f(x,y)=0 f(x,y)=0上的任意两点的差必为0,但显然在计算 f x ′ f'_x fx′的过程中,并没有遵守 f ( x , y ) = 0 f(x,y)=0 f(x,y)=0这个条件,即 ( x + Δ x , y ) (x+\Delta x,y) (x+Δx,y)可能并不在曲线 f ( x , y ) = 0 f(x,y)=0 f(x,y)=0上。
若想制造出一个完全位于 f ( x , y ) = 0 f(x,y)=0 f(x,y)=0上的变化关系,则可令
Δ f = f ( x 0 + Δ x , y 0 + Δ y ) − f ( x 0 , y 0 ) = f ( x 0 + Δ x , y 0 + Δ y ) − f ( x 0 , y 0 + Δ y ) + f ( x 0 , y 0 + Δ y ) − f ( x 0 , y 0 ) \begin{aligned} \Delta f&=f(x_0+\Delta x,y_0+\Delta y)-f(x_0,y_0)\\ &=f(x_0+\Delta x,y_0+\Delta y)-f(x_0,y_0+\Delta y)+f(x_0,y_0+\Delta y)-f(x_0,y_0) \end{aligned} Δf=f(x0+Δx,y0+Δy)−f(x0,y0)=f(x0+Δx,y0+Δy)−f(x0,y0+Δy)+f(x0,y0+Δy)−f(x0,y0)
则当 Δ x → 0 , Δ y → 0 \Delta x\to0,\Delta y\to0 Δx→0,Δy→0时,
Δ f = f x ′ ( x 0 , y 0 + Δ y ) Δ x + f y ′ ( x 0 , y 0 ) Δ y = f x ′ ( x 0 , y 0 ) Δ x + f y ′ ( x 0 , y 0 ) Δ y \begin{aligned} \Delta f&=f'_x(x_0,y_0+\Delta y)\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{aligned} Δf=fx′(x0,y0+Δy)Δx+fy′(x0,y0)Δy=fx′(x0,y0)Δx+fy′(x0,y0)Δy
由于满足 f ( x , y ) = 0 f(x,y)=0 f(x,y)=0,所以 Δ f = 0 \Delta f=0 Δf=0,且 Δ x , Δ y \Delta x,\Delta y Δx,Δy也会受到约束,二者同时趋近于0也可以表述为 Δ y \Delta y Δy在满足约束 f f f的情况下,随着 Δ → 0 \Delta\to0 Δ→0而趋近于0。故上式可写为
f x ′ ( x 0 , y 0 ) Δ x = − f y ′ ( x 0 , y 0 ) Δ y ↓ lim Δ x → 0 Δ y Δ x = − f x ′ f y ′ ↓ d y d x = − f x ′ f y ′ f'_x(x_0,y_0)\Delta x=-f'_y(x_0,y_0)\Delta y\\ \downarrow\\ \lim_{\Delta x\to0}\frac{\Delta y}{\Delta x}=-\frac{f'_x}{f'_y}\\ \downarrow\\ \frac{\text dy}{\text dx}=-\frac{f'_x}{f'_y} fx′(x0,y0)Δx=−fy′(x0,y0)Δy↓Δx→0limΔxΔy=−fy′fx′↓dxdy=−fy′fx′
故对于 f ( x , y ) = 0 f(x,y)=0 f(x,y)=0,其切线方向为 ( 1 , − f x ′ f y ′ ) (1,-\frac{f'_x}{f'_y}) (1,−fy′fx′)。
其法线方向为 ( f x ′ , f y ′ ) (f'_x,f'_y) (fx′,fy′)。
二者内积 ( 1 , − f x ′ f y ′ ) ⋅ ( f x ′ , f y ′ ) ≡ 0 (1,-\frac{f'_x}{f'_y})\cdot(f'_x,f'_y)\equiv0 (1,−fy′fx′)⋅(fx′,fy′)≡0。从而我们证明了 f ( x , y ) f(x,y) f(x,y)的梯度为 f ( x , y ) = 0 f(x,y)=0 f(x,y)=0这个曲线的法线。
而这个过程中,则得到了一种新的变化关系,即全微分,定义为
d f = f x ′ d x + f y ′ d y \text df=f'_x\text dx+f'_y\text dy df=fx′dx+fy′dy