Logo

深度专题：斯托克斯公式

斯托克斯公式的故事

从一道竞赛题到数学基石

1854年，爱尔兰数学家、物理学家乔治·斯托克斯（George Gabriel Stokes，1819-1903）在剑桥大学的史密斯奖考试中出了一道题，要求证明一个关于曲线积分和曲面积分的定理。这道题的结果后来以他的名字命名为斯托克斯公式。

然而，这个定理的思想最早可以追溯到开尔文勋爵（Lord Kelvin），他在1850年给斯托克斯的一封信中提到了类似的结果。斯托克斯将其整理并作为考题，因此获得了命名权。

斯托克斯是19世纪英国数学物理学派的代表人物之一，与开尔文、麦克斯韦等人共同推动了数学物理的发展。

“旋转”的数学表达

如果说高斯公式描述的是”发散”（通量与散度的关系），那么斯托克斯公式描述的就是”旋转”（环流量与旋度的关系）。它告诉我们：向量场沿闭曲线的环流量，等于旋度通过以该曲线为边界的任意曲面的通量。

斯托克斯公式的深入理解

公式的精确表述

定理1

设 Γ\Gamma 为分段光滑的空间有向闭曲线，Σ\Sigma 是以 Γ\Gamma 为边界的分片光滑的有向曲面，Γ\Gamma 的正向与 Σ\Sigma 的侧符合右手规则，函数 P,Q,RP, Q, R 在包含 Σ\Sigma 在内的一个空间区域内具有一阶连续偏导数，则有

∮ΓP dx+Q dy+R dz=∬Σ(∂R∂y−∂Q∂z)dy dz+(∂P∂z−∂R∂x)dz dx+(∂Q∂x−∂P∂y)dx dy\oint_\Gamma P \, dx + Q \, dy + R \, dz = \iint_\Sigma \left( \frac{\partial R}{\partial y} - \frac{\partial Q}{\partial z} \right) dy\,dz + \left( \frac{\partial P}{\partial z} - \frac{\partial R}{\partial x} \right) dz\,dx + \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) dx\,dy

或用向量形式：

∮ΓF⃗⋅dr⃗=∬Σ(∇×F⃗)⋅dS⃗\oint_\Gamma \vec{F} \cdot d\vec{r} = \iint_\Sigma (\nabla \times \vec{F}) \cdot d\vec{S}

几何解释
推论
证明
符号说明

右手规则

曲面的无关性

斯托克斯公式的一个重要性质是：积分值只依赖于边界曲线 Γ\Gamma，与以 Γ\Gamma 为边界的曲面 Σ\Sigma 的选择无关。这是因为对于两块具有相同边界的曲面 Σ1\Sigma_1 和 Σ2\Sigma_2，它们合起来构成一个闭曲面，而旋度的散度为零（∇⋅(∇×F⃗)=0\nabla \cdot (\nabla \times \vec{F}) = 0），由高斯公式知通过闭曲面的旋度通量为零，因此两块曲面的通量相等（符号相反，抵消后为零）。

证明思路

第一步：投影到坐标面

将曲面 Σ\Sigma 投影到三个坐标面，分别证明含 P,Q,RP, Q, R 的各项。以含 PP 的项为例：

∮ΓP dx=∬Σ∂P∂zdz dx−∂P∂ydx dy\oint_\Gamma P \, dx = \iint_\Sigma \frac{\partial P}{\partial z} dz\,dx - \frac{\partial P}{\partial y} dx\,dy

第二步：参数化曲面

设曲面 Σ\Sigma 由 z=z(x,y)z = z(x,y) 给出，利用复合函数求导法则和格林公式（在投影区域上）证明上式。

第三步：一般曲面

将一般曲面分割为若干个标准型曲面片，公共边界上的曲线积分相互抵消。

三大公式的统一

从维度看三大公式

广义斯托克斯定理

在微分形式的语言下，所有这些公式都统一为广义斯托克斯定理：

∫∂Mω=∫Mdω\int_{\partial M} \omega = \int_M d\omega

其中 MM 是一个定向流形，∂M\partial M 是其边界，ω\omega 是一个微分形式，dωd\omega 是其外微分。

  • 牛顿-莱布尼茨：MM 是区间，ω=f(x)\omega = f(x)（0-形式），dω=f′(x)dxd\omega = f'(x)dx（1-形式）
  • 格林公式：MM 是平面区域，ω=Pdx+Qdy\omega = Pdx + Qdy（1-形式），dω=(∂Q∂x−∂P∂y)dx∧dyd\omega = (\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y})dx\wedge dy（2-形式）
  • 高斯公式：MM 是空间区域，ω=Pdy∧dz+Qdz∧dx+Rdx∧dy\omega = Pdy\wedge dz + Qdz\wedge dx + Rdx\wedge dy（2-形式），dω=(∂P∂x+∂Q∂y+∂R∂z)dx∧dy∧dzd\omega = (\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z})dx\wedge dy\wedge dz（3-形式）
  • 斯托克斯公式：MM 是空间曲面，ω=Pdx+Qdy+Rdz\omega = Pdx + Qdy + Rdz（1-形式），dωd\omega 是2-形式

斯托克斯公式在物理学中的应用

1. 法拉第电磁感应定律

∮LE⃗⋅dr⃗=−ddt∬ΣB⃗⋅dS⃗\oint_L \vec{E} \cdot d\vec{r} = -\frac{d}{dt} \iint_\Sigma \vec{B} \cdot d\vec{S}

即：电场沿闭曲线的环流量等于磁通量变化率的负值。利用斯托克斯公式，其微分形式为：

∇×E⃗=−∂B⃗∂t\nabla \times \vec{E} = -\frac{\partial \vec{B}}{\partial t}

这是麦克斯韦方程组的第三个方程。

2. 安培-环路定理（含位移电流）

∮LB⃗⋅dr⃗=μ0∬Σ(J⃗+ε0∂E⃗∂t)⋅dS⃗\oint_L \vec{B} \cdot d\vec{r} = \mu_0 \iint_\Sigma \left( \vec{J} + \varepsilon_0 \frac{\partial \vec{E}}{\partial t} \right) \cdot d\vec{S}

其微分形式为：

∇×B⃗=μ0J⃗+μ0ε0∂E⃗∂t\nabla \times \vec{B} = \mu_0 \vec{J} + \mu_0\varepsilon_0 \frac{\partial \vec{E}}{\partial t}

这是麦克斯韦方程组的第四个方程，其中位移电流 ε0∂E⃗∂t\varepsilon_0 \frac{\partial \vec{E}}{\partial t} 是麦克斯韦的伟大贡献。

3. 流体力学中的环量

斯托克斯公式将流体的环流量与涡量（旋度）联系起来，是研究涡旋运动的基础。

典型例题

例题1：用斯托克斯公式计算空间曲线积分

计算 ∮Γ(y−z) dx+(z−x) dy+(x−y) dz\oint_\Gamma (y-z) \, dx + (z-x) \, dy + (x-y) \, dz，其中 Γ\Gamma 是圆柱面 x2+y2=a2x^2+y^2 = a^2 与平面 x+y+z=0x+y+z = 0 的交线，从 xx 轴正方向看去取逆时针方向。

参考答案(4 个标签)
斯托克斯公式空间曲线积分旋度椭圆

解题思路： 用斯托克斯公式转化为曲面积分，旋度为常数向量。

详细步骤：

  1. ∇×F⃗=(−2,−2,−2)\nabla \times \vec{F} = (-2, -2, -2)
  2. 取 Σ\Sigma 为平面 x+y+z=0x+y+z=0 上以 Γ\Gamma 为边界的椭圆面，法向量 n⃗=(1,1,1)3\vec{n} = \frac{(1,1,1)}{\sqrt{3}}（由右手定则确定）
  3. (∇×F⃗)⋅n⃗=−63=−23(\nabla \times \vec{F}) \cdot \vec{n} = \frac{-6}{\sqrt{3}} = -2\sqrt{3}
  4. 椭圆面积：投影到 xOyxOy 平面为圆 x2+y2≤a2x^2+y^2 \leq a^2，面积 πa2\pi a^2。cos⁡γ=13\cos\gamma = \frac{1}{\sqrt{3}}，椭圆面积 =πa2cos⁡γ=3πa2= \frac{\pi a^2}{\cos\gamma} = \sqrt{3}\pi a^2
  5. 积分 =−23×3πa2=−6πa2= -2\sqrt{3} \times \sqrt{3}\pi a^2 = -6\pi a^2

答案：−6πa2-6\pi a^2

例题2：旋度为零的场

验证向量场 F⃗=(yz,xz,xy)\vec{F} = (yz, xz, xy) 的旋度为零，并说明其曲线积分与路径无关。

参考答案(4 个标签)
斯托克斯公式无旋场有势场路径无关

解题思路： 计算旋度，利用斯托克斯公式说明沿任意闭曲线积分为零。

详细步骤：

  1. ∇×F⃗=(∂∂y(xy)−∂∂z(xz),∂∂z(yz)−∂∂x(xy),∂∂x(xz)−∂∂y(yz))=(x−x,y−y,z−z)=(0,0,0)\nabla \times \vec{F} = (\frac{\partial}{\partial y}(xy)-\frac{\partial}{\partial z}(xz), \frac{\partial}{\partial z}(yz)-\frac{\partial}{\partial x}(xy), \frac{\partial}{\partial x}(xz)-\frac{\partial}{\partial y}(yz)) = (x-x, y-y, z-z) = (0,0,0)
  2. 对任意闭曲线 Γ\Gamma，取以 Γ\Gamma 为边界的曲面 Σ\Sigma，由斯托克斯公式： ∮ΓF⃗⋅dr⃗=∬Σ(∇×F⃗)⋅dS⃗=∬Σ0⋅dS⃗=0\oint_\Gamma \vec{F} \cdot d\vec{r} = \iint_\Sigma (\nabla \times \vec{F}) \cdot d\vec{S} = \iint_\Sigma 0 \cdot d\vec{S} = 0
  3. 沿任意闭曲线积分为零，因此曲线积分与路径无关。
  4. 势函数：u=xyzu = xyz（因为 ∇u=(yz,xz,xy)=F⃗\nabla u = (yz, xz, xy) = \vec{F}）

答案：旋度为零，是有势场，势函数为 u=xyzu = xyz，曲线积分与路径无关。

练习题

练习1

用斯托克斯公式计算 ∮Γz dx+x dy+y dz\oint_\Gamma z \, dx + x \, dy + y \, dz，其中 Γ\Gamma 是平面 x+y+z=1x+y+z=1 被三个坐标面截成的三角形的边界，正向与平面上侧成右手规则。

参考答案(4 个标签)
斯托克斯公式三角形旋度平面

解题思路： 旋度为 (1,1,1)(1,1,1)，法向量为 (1,1,1)3\frac{(1,1,1)}{\sqrt{3}}。

详细步骤：

  1. ∇×F⃗=(1,1,1)\nabla \times \vec{F} = (1, 1, 1)
  2. 法向量 n⃗=(1,1,1)3\vec{n} = \frac{(1,1,1)}{\sqrt{3}}
  3. (∇×F⃗)⋅n⃗=3(\nabla \times \vec{F}) \cdot \vec{n} = \sqrt{3}
  4. 三角形面积：三个顶点 (1,0,0),(0,1,0),(0,0,1)(1,0,0), (0,1,0), (0,0,1)，边长 2\sqrt{2}，面积 =34(2)2=32= \frac{\sqrt{3}}{4}(\sqrt{2})^2 = \frac{\sqrt{3}}{2}
  5. 积分 =3×32=32= \sqrt{3} \times \frac{\sqrt{3}}{2} = \frac{3}{2}

答案：32\frac{3}{2}

练习2

计算 ∮Γ3y dx−xz dy+yz2 dz\oint_\Gamma 3y \, dx - xz \, dy + yz^2 \, dz，其中 Γ\Gamma 是圆周 x2+y2=2z, z=2x^2+y^2=2z, \, z=2，从 zz 轴正方向看去取逆时针方向。

参考答案(4 个标签)
斯托克斯公式圆周旋度平面圆面

解题思路： 取 z=2z=2 平面上的圆面，用斯托克斯公式。

详细步骤：

  1. ∇×F⃗=(∂∂y(yz2)−∂∂z(−xz),∂∂z(3y)−∂∂x(yz2),∂∂x(−xz)−∂∂y(3y))=(z2+x,0,−z−3)\nabla \times \vec{F} = (\frac{\partial}{\partial y}(yz^2)-\frac{\partial}{\partial z}(-xz), \frac{\partial}{\partial z}(3y)-\frac{\partial}{\partial x}(yz^2), \frac{\partial}{\partial x}(-xz)-\frac{\partial}{\partial y}(3y)) = (z^2+x, 0, -z-3)
  2. 在 z=2z=2 平面上：∇×F⃗=(4+x,0,−5)\nabla \times \vec{F} = (4+x, 0, -5)
  3. 取 Σ:z=2, x2+y2≤4\Sigma: z=2, \, x^2+y^2 \leq 4，上侧，法向量 (0,0,1)(0,0,1)
  4. (∇×F⃗)⋅n⃗=−5(\nabla \times \vec{F}) \cdot \vec{n} = -5
  5. 积分 =∬Σ(−5)dS=−5×π⋅22=−20π= \iint_\Sigma (-5) dS = -5 \times \pi \cdot 2^2 = -20\pi

答案：−20π-20\pi

总结

本文出现的符号

符号类型读音/说明在本文中的含义
Γ\Gamma空间闭曲线space closed curve斯托克斯公式中的边界曲线
Σ\Sigma有向曲面oriented surface以 Γ\Gamma 为边界的曲面
∇×F⃗\nabla \times \vec{F}旋度curl向量场的旋度
dr⃗d\vec{r}有向弧元素directed arc element曲线的有向弧微元
dS⃗d\vec{S}有向面积元素directed surface element曲面的有向面积微元
ω\omega微分形式differential form外微分形式
dωd\omega外微分exterior derivative微分形式的外微分
∂M\partial M流形边界boundary of manifold定向流形的边界

中英对照

中文术语英文术语音标说明
斯托克斯公式Stokes’ theorem/stoʊksɪz ˈθɪərəm/空间曲线积分与曲面积分的关系
右手规则right-hand rule/raɪt hænd ruːl/确定曲线正向与曲面侧的关系
环流量circulation/ˌsɜːrkjəˈleɪʃən/向量场沿闭曲线的积分
旋度curl / rotation/kɜːrl/ /roʊˈteɪʃən/向量场的旋转程度
广义斯托克斯定理generalized Stokes theorem/ˈdʒɛnrəlaɪzd stoʊks ˈθɪərəm/微分形式的统一积分定理
微分形式differential form/ˌdɪfəˈrɛnʃəl fɔːrm/外微分形式
外微分exterior derivative/ɪkˈstɪəriər dɪˈrɪvətɪv/微分形式的微分运算
法拉第定律Faraday’s law/ˈfærədeɪz lɔː/电磁感应定律
位移电流displacement current/dɪsˈpleɪsmənt ˈkɜːrənt/麦克斯韦引入的电流项
安培环路定理Ampère’s circuital law/æmˈpɛərz ˈsɜːrkjɪtəl lɔː/磁场环流量与电流的关系
涡量vorticity/vɔːrˈtɪsəti/流速场的旋度
上一章节深度专题：高斯公式深入探讨高斯公式（散度定理）的历史背景、物理意义、证明思路及其在电磁学和流体力学中的应用。
下一章节深度专题：雅可比行列式深入探讨雅可比行列式的几何意义、计算方法、在变量替换中的应用，以及与线性变换和微分几何的联系。
多元函数积分学