21.梯度场和势函数

1 判断梯度场

上一节对于向量场判断是否为梯度场已有较为全面的表述

本小节将从梯度场的性质出发,展示最常见且便捷的判断方法

如果F是梯度场,则F=f,M=fx,N=fy

fxy=fyx可知,梯度场F需要满足:My=Nx

  • F是梯度场,是My=Nx的充分不必要条件
  • My=NxF处处有定义可导的前提下,可推得F是梯度场
  • 在后续的第24节课中将提及,这一前提是在限定F定义域为单连通区域

以上一课时的第三小节例题说明:F=<y,x>

M=y,N=x,由于MyNx,势函数不存在,所以F不是梯度场

2 寻找势函数

寻找势函数的前提是My=Nx

方法一:计算线积分

易知cFdr=f(x,y)f(0,0),所以f(x,y)=cFdr+f(0,0)

考虑到梯度场路径独立的性质,路径c可拆分成路径c1和路径c2

所以f(x,y)=c1Fdr+c2Fdr+f(0,0)

  • 计算路径c1时,y=0,dy=0,0xx1
  • 计算路径c2时,x=x1,dx=0,0yy1
  • 借助路径独立的性质,计算过程得到了极大简化

方法二:求不定积分

势函数满足以下形式:(fx)dx=F(x)+g(y),(fy)dx=F(y)+g(x)

所以可以先根据fx=M(x,y)还原势函数中包含x的部分F(x),保留仅包含y的未知部分为g(y),再考虑根据fy=N(x,y)还原势函数中仅包含y的未知部分g(y)

此方法较难描述,可参考下一节的示例进行直观理解

3 寻找势函数的示例

例题:求解梯度场F=<4x2+8xy,3y2+4x2>的势函数

方法一(计算线积分): f(x,y)=c1Fdr+c2Fdr+f(0,0) =0x1(4x2+8x0)dx+0y1(3y2+4x12)dy+f(0,0) =43x13+y13+4x12y1+C 方法二(求不定积分):

(fx)dx=F(x)+g(y) =43x3+4x2y+g(y) fy=3y2+4x2=4x2+g(y),g(y)=y3+C f(x,y)=43x3+y3+4x2y+C

相比较来说,方法二更为简单常用 对于非线性复杂梯度场,求解势函数是困难的

4 旋度以及理解

在物理学中,常用向量场的旋度判断向量场是否为保守场

定义旋度的符号为curl,则旋度的计算方法定义如下: curl(F)=NxMy 容易发现,当Nx=My时,向量场F的旋度为0,说明F是保守场

  • My=Nx在单连通区域成立的前提下,可推得F是梯度场
  • 梯度场是保守场,所以保守场都是无旋的(旋度为0)

在速度场中,旋度可以衡量运动在各个点上的扭转程度大小,旋度的值等于运动在旋转分量上的角速度的两倍,而旋度的正/负则分别对应顺/逆时针旋转

在力场中,旋度可以衡量物体上任一点受到的扭矩(力用于转动的分量),类似于加速度=力/质量的关系,角加速度=扭矩/转动惯量

5 参考:

博客园-梯度场和势函数

往年同期文章