一文读懂Nabla算子

$\nabla$算子的介绍

$\nabla$称作Nabla算子或del算子（算子是一种映射，可以理解为特殊的函数），$\nabla$算子从功能上来看，被称作矢量微分算子，又因为它是被哈密顿引入的，因此也称哈密顿算子。其他名称还有劈形算子、倒三角算子等等，不过这都不重要。

那么这个$\nabla$算子到底是什么呢，我们先来看它在三维情况下的定义：
$$\nabla =\left(\frac{\partial }{\partial x},\frac{\partial }{\partial y},\frac{\partial }{\partial z}\right)={\vec{e}}_x\frac{\partial }{\partial x}+{\vec{e}}_y\frac{\partial }{\partial y}+{\vec{e}}_z\frac{\partial }{\partial z}$$

只看定义的话不难发现，$\nabla$算子本身似乎并没有什么明确的意义，它有偏导符号也有一个指向不明的方向，但我们似乎说不出它对谁求导以及它指向何处，就像微分算子$D=\frac{d}{dx}$一样，$\nabla$必须作用于函数才有意义。那么$\nabla$算子应当如何作用于其他函数呢？

梯度、散度和旋度

那么我们来观察一下$\nabla$的性质以确定其如何作用于函数，根据其定义我们可以发现，$\nabla$算子与矢量的形式$\vec{A}=\overrightarrow{e_x}a+\overrightarrow{e_y}b+\overrightarrow{e_z}c$也有一定的相似性，虽然它本身并不是矢量，但我们仍可以将其当做矢量进行使用，那么与矢量类似，$\nabla$算子也应当有三种乘法。

矢量$\vec{A}$的三种乘法 1. 与标量$a$相乘得$a\vec{A}$；2. 与另一个矢量$\vec{B}$点乘得 $\vec{A}\cdot \vec{B}$；3. 与另一个矢量叉乘得 $\vec{A}\times \vec{B}$。

所以我们将$\nabla$算子也作用于标量函数$u$上得 $\nabla u$，那么根据$\nabla$算子的定义有，$\nabla u=\left(\frac{\partial }{\partial x},\frac{\partial }{\partial y},\frac{\partial }{\partial z}\right)u=\frac{\partial u}{\partial x}\overrightarrow{e_x}+\frac{\partial u}{\partial y}\overrightarrow{e_y}+\frac{\partial u}{\partial z}\overrightarrow{e_z}$，回忆一下高数内容可以发现，这不就是梯度的表达式么，所以$\nabla$算子作用于标量函数$\nabla$即可得到标量函数的梯度，$\nabla u=gradu$。

接下来如法炮制将$\nabla$算子与矢量函数$\vec{f}$点乘，那么我们可以得到$\nabla \cdot \vec{f}=\left(\frac{\partial }{\partial x},\frac{\partial }{\partial y},\frac{\partial }{\partial z}\right)\cdot \left(f_x,f_y,f_z\right)=\frac{\partial f_x}{\partial x}\overrightarrow{e_x}+\frac{\partial f_y}{\partial y}\overrightarrow{e_y}+\frac{\partial f_z}{\partial z}\overrightarrow{e_z}$，这被称作$\vec{f}$的散度，高斯定律就是电场的散度为0，散度表征空间各点矢量场发散的强弱程度，$\nabla \cdot \vec{f}=divf$。

再将$\nabla$算子与矢量函数$\vec{A}$叉乘，则$\nabla \times \vec{A}=\left(\frac{\partial }{\partial x},\frac{\partial }{\partial y},\frac{\partial }{\partial z}\right)\times \left(A_x,A_y,A_z\right)=\left|\begin{array}{l}{\vec{e}}_x{\vec{e}}_y{\vec{e}}_z\\ \frac{\partial }{\partial x}\frac{\partial }{\partial y}\frac{\partial }{\partial z}\\ A_x{A}_y{A}_z\end{array}\right|$，这被称作$\vec{A}$的旋度，旋度表示向量场对某一点附近的微元造成的旋转程度，$\nabla \times \vec{A}=rot\vec{A}$。

散度反映的是无穷小曲面的通量和闭合曲面体积的比值，这里借用长尾科技的比喻体会散度的物理意义：在电场中，把电场线都想象成水流，然后拿一个非常轻的圆形橡皮筋放到这里，如果这个橡皮筋的面积变大，我们就说这个点的散度为正，反之为负。

旋度的定义就是无穷小非闭合曲面的环流和这个曲面的面积之比。旋度的物理意义：我们把一个小风车放在某一点上，如果这个风车能转起来，就说明这点的旋度不为0。

对于旋度、散度这两个概念来源不清楚的同学可以观看参考文献[2]的文章，想要更直观了解散度、旋度物理意义的同学可以在参考文献及视频[3]中观看3b1b的可视化视频，可以直接从14:11处开始。

不同坐标系下$\nabla$算子的形式

我们在定义中给出的是$\nabla$算子在三维直角坐标系下的定义 $\nabla =\left(\frac{\partial }{\partial x},\frac{\partial }{\partial y},\frac{\partial }{\partial z}\right)={\vec{e}}_x\frac{\partial }{\partial x}+{\vec{e}}_y\frac{\partial }{\partial y}+{\vec{e}}_z\frac{\partial }{\partial z}$。

向量微分算子的形式化定义为：$\nabla =\frac{d}{dr}$，在$n$维空间中，分母 $r$ 即为含$n$个分量的变量，在二维情况下$\nabla$算子即为 $\nabla =\left(\frac{\partial }{\partial x},\frac{\partial }{\partial y},\right)={\vec{e}}_x\frac{\partial }{\partial x}+{\vec{e}}_y\frac{\partial }{\partial y}+{\vec{e}}_z$，以此类推，可以得到任意高维的$\nabla$算子的形式。

在正交曲线坐标系下有
$$\begin{array}{rl}& \nabla ={\vec{e}}_1\frac{1}{h_1}\frac{\partial }{\partial x_1}+{\vec{e}}_2\frac{1}{h_2}\frac{\partial }{\partial x_2}+{\vec{e}}_3\frac{1}{h_3}\frac{\partial }{\partial x_3}\\ & \nabla \phi ={\vec{e}}_1\frac{1}{h_1}\frac{\partial \phi }{\partial x_1}+{\vec{e}}_2\frac{1}{h_2}\frac{\partial \phi }{\partial x_2}+{\vec{e}}_3\frac{1}{h_3}\frac{\partial \phi }{\partial x_3}\\ & \nabla \cdot \vec{A}=\frac{1}{h_1{h}_2{h}_3}\left[\frac{\partial }{\partial x_1}\left(h_2{h}_3{A}_1\right)+\frac{\partial }{\partial x_2}\left(h_3{h}_1{A}_2\right)+\frac{\partial }{\partial x_3}\left(h_2{h}_1{A}_3\right)\right]\end{array}$$
其中
$$h_i=\sqrt{{\left(\frac{\partial x}{\partial x_i}\right)}^2+{\left(\frac{\partial y}{\partial x_i}\right)}^2+{\left(\frac{\partial z}{\partial x_i}\right)}^2}(i=1,2,3)$$
称度量系数（或拉梅系数）, 正交坐标系可完全由三个拉梅系数 $h_1,h_2,h_3$ 来描述。

$\nabla$算子运算律

通过之前的分析，我们可以发现$\nabla$算子既是一个微分算符，同时也具有矢量形式，初学者往往不知道该在什么条件下使用何种性质，所以运算起来会比较复杂。

下面给出简单的运算律与相应的记忆方法

根据前面提到的具有的微分运算性质与具有矢量形式的性质可以将它的运算拆分为两部分：

• 微分运算，即按照乘积的微分(导数)公式：$\frac{\mathrm{d}}{\mathrm{d}x}(fg)=\frac{\mathrm{d}f}{\mathrm{d}x}g+f\frac{\mathrm{d}g}{\mathrm{d}x}$ 变形；
• 变形之后只要当作矢量去化简就好了，相当于已经考虑完了微分的性质，化简就是想办法利用标量矢量混合运算让 $\nabla$ 作用在该作用的地方。

下面我们对几个常用公式进行推导

首先简单体会一下微分运算的变形
$$\nabla (\phi \psi )=\nabla \left(\phi {\psi }_c\right)+\nabla \left({\phi }_c\psi \right)=\psi \nabla \phi +\phi \nabla \psi$$
中间部分就是微分的作用，下标 $c$ 的意思是 $\nabla$ 不作用于它，当它是常数。
然后就是第二步： $\nabla$ 是个矢量，$\phi ,\psi$ 都是标量所以就只要让 $\nabla$ 作用在该作用的部分即可。

引入向量后再次体会微分运算
$$\begin{array}{rl}& \nabla \cdot (\phi \vec{f})=\nabla \cdot \left(\phi {\vec{f}}_c\right)+\nabla \cdot \left({\phi }_c\vec{f}\right)=\nabla \phi \cdot \vec{f}+\phi \nabla \cdot \vec{f}\\ & \nabla \times (\phi \vec{f})=\nabla \times \left(\phi {\vec{f}}_c\right)+\nabla \times \left({\phi }_c\vec{f}\right)=\nabla \phi \times \vec{f}+\phi \nabla \times \vec{f}\end{array}$$
第一步，先让 $\nabla$ 分别作用，然后第二步，$\nabla$ 与 $\vec{f}$ 是矢量， $\phi$ 是标量，按照传统的运算手段整理它们即可。

接下来的例子为微分+矢量运算
$$\nabla \cdot (\vec{f}\times \vec{g})=\nabla \cdot \left(\vec{f}\times {\vec{g}}_c\right)+\nabla \cdot \left({\vec{f}}_c\times \vec{g}\right)=(\nabla \times \vec{f})\cdot \vec{g}-(\nabla \times \vec{g})\cdot \vec{f}$$
我们看第二步，这就是前面提到的想办法利用标量矢量混合运算让 $\nabla$ 作用在该作用的地方。大家都是矢量，那就要按照矢量的规则来，这是混合积，具有轮换对称性。

同样是微分+矢量运算，且需要配合矢量公式计算
$$\begin{array}{rl}\nabla \times (\vec{f}\times \vec{g})& =\nabla \times \left(\vec{f}\times {\vec{g}}_c\right)+\nabla \times \left({\vec{f}}_c\times \vec{g}\right)\\ & =(\vec{g}\cdot \nabla )\vec{f}-(\nabla \cdot \vec{f})\vec{g}+(\nabla \cdot \vec{g})\vec{f}-(\vec{f}\cdot \nabla )\vec{g}\end{array}$$
由矢量公式：$\vec{A}\times (\vec{B}\times \vec{C})=\vec{B}(\vec{A}\cdot \vec{C})-\vec{C}(\vec{A}\cdot \vec{B})$，先对上式进行微分变形，再使用矢量公式即可得到最后公式。

很好，继续下一个例子，仍是微分+矢量运算，并配合矢量公式
$$\begin{array}{rl}\nabla (\vec{f}\cdot \vec{g})& =\nabla \left(\vec{f}\cdot {\vec{g}}_c\right)+\nabla \left({\vec{f}}_c\cdot \vec{g}\right)\\ & =\left[\left({\vec{g}}_c\cdot \nabla \right)\vec{f}-{\vec{g}}_c\times (\vec{f}\times \nabla )\right]+\left[\left({\vec{f}}_c\cdot \nabla \right)\vec{g}-{\vec{f}}_c\times (\vec{g}\times \nabla )\right]\\ & =(\vec{g}\cdot \nabla )\vec{f}+\vec{g}\times (\nabla \times \vec{f})+(\vec{f}\cdot \nabla )\vec{g}+\vec{f}\times (\nabla \times \vec{g})\end{array}$$
所需矢量公式为：$\vec{C}(\vec{A}\cdot \vec{B})=\vec{B}(\vec{A}\cdot \vec{C})-\vec{A}\times (\vec{B}\times \vec{C})$，先对上式进行微分变形，然后使用矢量公式做进一步变形，对其进行整理则可得最后公式。

$\nabla$的微分运算性质并不只在$\frac{\mathrm{d}}{\mathrm{d}x}(fg)=\frac{\mathrm{d}f}{\mathrm{d}x}g+f\frac{\mathrm{d}g}{\mathrm{d}x}变$形时使用，例如变形 $(\vec{A}\cdot \nabla )\vec{A}=\nabla (\frac{{\vec{A}}^2}{2})$ 中就使用了其微分性质。

通过前面介绍的运算律，我们可以推导更多$\nabla$算子相关公式，接下来给出总表。

$\nabla$算子常用公式

$\nabla (\phi \psi )=\psi \nabla \phi +\phi \nabla \psi$；

$\nabla \cdot (\phi \vec{f})=\nabla \phi \cdot \vec{f}+\phi \nabla \cdot \vec{f}$；

$\nabla \times (\phi \vec{f})=\nabla \phi \times \vec{f}+\phi \nabla \times \vec{f}$；

$\nabla \cdot (\vec{f}\times \vec{g})=(\nabla \times \vec{f})\cdot \vec{g}-(\nabla \times \vec{g})\cdot \vec{f}$；

$\nabla \times (\vec{f}\times \vec{g})=(\vec{g}\cdot \nabla )\vec{f}-(\nabla \cdot \vec{f})\vec{g}+(\nabla \cdot \vec{g})\vec{f}-(\vec{f}\cdot \nabla )\vec{g}$；

$\nabla (\vec{f}\cdot \vec{g})=\vec{f}\times (\nabla \times \vec{g})+(\vec{f}\cdot \nabla )\vec{g}+\vec{g}\times (\nabla \times \vec{f})+(\vec{g}\cdot \nabla )\vec{f}$；

$\nabla \cdot (\vec{f}\vec{g})=(\nabla \cdot \vec{f})\vec{g}+(\vec{f}\cdot \nabla )\vec{g}$；

$\nabla \cdot \nabla \phi \equiv {\nabla }^2\phi$；

$\nabla \times (\nabla \times \vec{f})=\nabla (\nabla \cdot \vec{f})-{\nabla }^2\vec{f}$；

$\nabla \cdot (\vec{f}\vec{g})=\nabla \cdot \left(\vec{f}{\vec{g}}_c\right)+\nabla \cdot \left({\vec{f}}_c\vec{g}\right)=(\nabla \cdot \vec{f})\vec{g}+(\vec{f}\cdot \nabla )\vec{g}$；

$\nabla \phi =\nabla \cdot (\phi \overleftrightarrow{I})$；

$\nabla \cdot \overleftrightarrow{I}=\frac{\partial }{\partial x}\left({\vec{e}}_x\cdot \overleftrightarrow{I}\right)+\frac{\partial }{\partial y}\left({\vec{e}}_y\cdot \overleftrightarrow{I}\right)+\frac{\partial }{\partial z}\left({\vec{e}}_z\cdot \overleftrightarrow{I}\right)$；

$\overleftrightarrow{I}$ 是张量。

$\nabla$算子相关的常用公式

$\begin{array}{rl}& \nabla \times (\nabla \phi )\equiv 0\\ & \nabla \cdot (\nabla \times \vec{A})\equiv 0\\ & \nabla f(u)=(\nabla u)\frac{\mathrm{d}f}{\mathrm{d}u}\\ & \nabla \cdot \vec{A}(u)=(\nabla u)\cdot \frac{\mathrm{d}\vec{A}}{\mathrm{d}u}\\ & \nabla \times \vec{A}(u)=(\nabla u)\times \frac{\mathrm{d}\vec{A}}{\mathrm{d}u}\\ & {\nabla }^2\frac{1}{r}=-4\pi \delta (r)\end{array}$

最后一行也可是 $\nabla \cdot \left(\frac{\vec{r}}{r^3}\right)=4\pi \delta (r)$

补充内容

我们前面提到$\nabla$直接作用于标量函数$u$表示求其梯度，求梯度的结果为一个向量(梯度有方向嘛，所以结果自然是向量)，而$\nabla$也可以作用于向量函数$\vec{A}$，这种运算同样是叫做求向量函数$\vec{A}$的梯度，对向量求梯度的结果为二阶张量。

${\nabla }^2\equiv \frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}+\frac{{\partial }^2}{\partial z^2}=\nabla \cdot \nabla$称为拉普拉斯算符/算子，${\nabla }^{2}u=\nabla \cdot (\nabla u)$

在任意方向 $\vec{l}$ 上移动线元距离 $dl$，$\phi$ 的增量 $d\phi$ 称为方向微分，即 $d\phi =\frac{\partial \phi }{\partial l}dl=\nabla \phi \cdot d\vec{l}$。显然，任意两点 $\phi$ 值差为
$${\phi }_B-{\phi }_A={\int }_A^B\nabla \phi \cdot d\vec{l}$$

$(\vec{A}\cdot \nabla )\varphi$是什么？

$(\vec{A}\cdot \nabla )\varphi =\vec{A}\cdot (\nabla \varphi )$，是标量场$\varphi$在$\vec{A}$方向上的方向导数的$\mid \vec{A}\mid$倍。

$(\vec{A}\cdot \nabla )\vec{B}$是什么？

很多初学者会被$(\vec{A}\cdot \nabla )\vec{B}$这样的公式所迷惑，尤其是其中的$\vec{A}\cdot \nabla$，这可能会导致这样的误解，难道$\nabla$还可以作用于左边的函数吗？其实$\nabla$左乘矢量有结果，但结果就和$\nabla$一样是一种算子，必须作用在具体的函数上才有意义。

令$\vec{A}=(A_x,A_y,A_z),\nabla =\left(\frac{\partial }{\partial x},\frac{\partial }{\partial y},\frac{\partial }{\partial z}\right)$，因此两者点乘为$\vec{A}\cdot \nabla =\left(A_x\frac{\partial }{\partial x},A_y\frac{\partial }{\partial y},A_z\frac{\partial }{\partial z}\right)$，那么这就是新的算子$\vec{A}\cdot \nabla$，所以$(\vec{A}\cdot \nabla )\vec{B}=\left(A_x\frac{\partial B_x}{\partial x},A_y\frac{\partial B_y}{\partial y},A_z\frac{\partial B_z}{\partial z}\right)$。

同时在涉及到$\nabla$的运算中还需考虑到$\nabla$在对谁作用，这点至关重要。

例如在$\vec{A}\times (\nabla \times \vec{B})$，我们不能将其按照矢量公式直接变形为$\nabla (\vec{A}\cdot \vec{B})-(\vec{A}\cdot \nabla )\vec{B}$，而是要写为${\nabla }_B(\vec{A}\cdot \vec{B})-(\vec{A}\cdot \nabla )\vec{B}$，${\nabla }_B$表示$\nabla$算子只作用于$B$。
PS：这点是笔者临时补充进来的，暂时还没有严谨易懂的证明方法，读者中有大佬的可以在评论区内留言补充。

感谢您的观看，如果您觉得我的文章有帮助，可以点个关注哦，顺便可以留个赞​😁，如果对我的文章有什么建议，也欢迎在评论区留言讨论。

更多相关内容传送门：数学物理方法专栏阅读指南

参考文献及视频

[1] 长尾科技 最美的公式：你也能懂的麦克斯韦方程组（微分篇）

[2] 3Blue1Brown 散度与旋度：麦克斯韦方程组、流体等所用到的语言

[3] 東雲正樹 Nabla 算符 ∇ 的运算律以及常用公式