Logo

深度专题：高斯公式

高斯公式的故事

数学王子的贡献

卡尔·弗里德里希·高斯（Carl Friedrich Gauss，1777-1855），被誉为”数学王子”，是历史上最伟大的数学家之一。他在数论、代数、分析、几何、天文学、大地测量学等领域都有开创性贡献。

高斯公式（散度定理）最早由俄国数学家奥斯特罗格拉德斯基（Mikhail Ostrogradsky）在1826年正式发表，因此在俄语世界中被称为奥斯特罗格拉德斯基公式。但高斯在更早的研究中已经使用了类似的结果，因此通常以高斯的名字命名。

从”流量”到”散度”

高斯公式的物理直觉来自流体力学：如果流体在空间中流动，那么从一个封闭曲面流出的总流量，应该等于曲面内部所有”源”产生的流量减去所有”汇”消耗的流量。散度就是单位体积的”源强度”。

高斯公式的深入理解

公式的精确表述

定理1

设空间闭区域 Ω\Omega 由分片光滑的闭曲面 Σ\Sigma 围成，函数 P(x,y,z),Q(x,y,z),R(x,y,z)P(x,y,z), Q(x,y,z), R(x,y,z) 在 Ω\Omega 上具有一阶连续偏导数，则有

∯ΣP dy dz+Q dz dx+R dx dy=∭Ω(∂P∂x+∂Q∂y+∂R∂z)dv\oiint_\Sigma P \, dy\,dz + Q \, dz\,dx + R \, dx\,dy = \iiint_\Omega \left( \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} + \frac{\partial R}{\partial z} \right) dv

或用向量形式：

∯ΣF⃗⋅dS⃗=∭Ω(∇⋅F⃗) dv\oiint_\Sigma \vec{F} \cdot d\vec{S} = \iiint_\Omega (\nabla \cdot \vec{F}) \, dv

其中 Σ\Sigma 是 Ω\Omega 的整个边界曲面的外侧。

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

外侧的规定

散度的极限定义

散度可以用通量的极限来定义：

∇⋅F⃗(M)=lim⁡ΔV→01ΔV∯ΔΣF⃗⋅dS⃗\nabla \cdot \vec{F}(M) = \lim_{\Delta V \to 0} \frac{1}{\Delta V} \oiint_{\Delta\Sigma} \vec{F} \cdot d\vec{S}

其中 ΔV\Delta V 是包含点 MM 的小区域的体积，ΔΣ\Delta\Sigma 是其边界。这个定义揭示了散度的物理意义：单位体积的净流出通量。

证明思路

高斯公式的证明与格林公式类似，分为三步：

第一步：XY型区域

设 Ω\Omega 是 XYXY 型区域：(x,y)∈Dxy, z1(x,y)≤z≤z2(x,y)(x,y) \in D_{xy}, \, z_1(x,y) \leq z \leq z_2(x,y)。

∭Ω∂R∂zdv=∬Dxydx dy∫z1z2∂R∂zdz=∬Dxy[R(x,y,z2)−R(x,y,z1)]dx dy\iiint_\Omega \frac{\partial R}{\partial z} dv = \iint_{D_{xy}} dx\,dy \int_{z_1}^{z_2} \frac{\partial R}{\partial z} dz = \iint_{D_{xy}} [R(x,y,z_2) - R(x,y,z_1)] dx\,dy

而沿边界的曲面积分 ∯ΣR dx dy\oiint_\Sigma R \, dx\,dy 由三部分组成：上曲面（上侧）、下曲面（下侧）、侧面（法向量与 zz 轴垂直，投影为零）。计算结果与上式相等。

第二步：其他类型区域

类似地证明含 PP 和 QQ 的项。

第三步：一般区域

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

高斯公式的推广

非闭曲面的补面法

当曲面 Σ\Sigma 不是闭曲面时，可以补上一块曲面 Σ′\Sigma' 使其成为闭曲面，应用高斯公式后再减去补上的曲面积分：

∬ΣF⃗⋅dS⃗=∭Ω(∇⋅F⃗)dv−∬Σ′F⃗⋅dS⃗\iint_\Sigma \vec{F} \cdot d\vec{S} = \iiint_\Omega (\nabla \cdot \vec{F}) dv - \iint_{\Sigma'} \vec{F} \cdot d\vec{S}

高斯公式的微分形式

从积分形式出发，利用积分中值定理，令区域收缩到一点，可以得到高斯公式的微分形式：

∇⋅F⃗=lim⁡ΔV→01ΔV∯ΔΣF⃗⋅dS⃗\nabla \cdot \vec{F} = \lim_{\Delta V \to 0} \frac{1}{\Delta V} \oiint_{\Delta\Sigma} \vec{F} \cdot d\vec{S}

这就是散度的定义式。

高斯公式在物理学中的应用

1. 电磁学中的高斯定律

电场的高斯定律： ∯ΣE⃗⋅dS⃗=Q内ε0\oiint_\Sigma \vec{E} \cdot d\vec{S} = \frac{Q_{\text{内}}}{\varepsilon_0}

这是麦克斯韦方程组的第一个方程，表明通过任意闭曲面的电通量等于曲面内的总电荷量除以真空介电常数。利用高斯公式，其微分形式为：

∇⋅E⃗=ρε0\nabla \cdot \vec{E} = \frac{\rho}{\varepsilon_0}

其中 ρ\rho 是电荷体密度。

磁场的高斯定律： ∯ΣB⃗⋅dS⃗=0\oiint_\Sigma \vec{B} \cdot d\vec{S} = 0

其微分形式为 ∇⋅B⃗=0\nabla \cdot \vec{B} = 0，表明磁单极不存在。

2. 流体力学中的连续性方程

对于不可压缩流体，密度 ρ\rho 为常数，质量守恒定律给出：

∯Σρv⃗⋅dS⃗=0⇒∇⋅v⃗=0\oiint_\Sigma \rho \vec{v} \cdot d\vec{S} = 0 \Rightarrow \nabla \cdot \vec{v} = 0

对于可压缩流体，连续性方程为：

∂ρ∂t+∇⋅(ρv⃗)=0\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \vec{v}) = 0

3. 热传导中的傅里叶定律

热流密度 q⃗=−k∇T\vec{q} = -k\nabla T，能量守恒给出：

∂∂t(ρcT)=∇⋅(k∇T)\frac{\partial}{\partial t}(\rho c T) = \nabla \cdot (k\nabla T)

这就是热传导方程。

典型例题

例题1：点电荷的电通量

验证点电荷 qq 产生的电场 E⃗=q4πε0r3r⃗\vec{E} = \frac{q}{4\pi\varepsilon_0 r^3}\vec{r} 通过包围原点的任意闭曲面外侧的通量为 qε0\frac{q}{\varepsilon_0}。

参考答案(5 个标签)
高斯公式电场点电荷补面法奇点

解题思路： 电场在原点处有奇点，不能直接用高斯公式。用补面法：在闭曲面内作小球面，在两曲面之间的区域用高斯公式。

详细步骤：

  1. 在 r≠0r \neq 0 处，∇⋅E⃗=0\nabla \cdot \vec{E} = 0（已在调和场一节验证）。
  2. 设 Σ\Sigma 是包围原点的任意闭曲面（外侧），Σε\Sigma_\varepsilon 是半径为 ε\varepsilon 的小球面（内侧，即法向量指向原点），位于 Σ\Sigma 内部。
  3. 在 Σ\Sigma 和 Σε\Sigma_\varepsilon 之间的区域 Ω′\Omega' 上用高斯公式： ∯Σ∪ΣεE⃗⋅dS⃗=∭Ω′(∇⋅E⃗)dv=0\oiint_{\Sigma \cup \Sigma_\varepsilon} \vec{E} \cdot d\vec{S} = \iiint_{\Omega'} (\nabla \cdot \vec{E}) dv = 0
  4. 所以 ∬ΣE⃗⋅dS⃗=−∬ΣεE⃗⋅dS⃗\iint_\Sigma \vec{E} \cdot d\vec{S} = -\iint_{\Sigma_\varepsilon} \vec{E} \cdot d\vec{S}
  5. 在 Σε\Sigma_\varepsilon 上（外侧）：E⃗=q4πε0ε2e⃗r\vec{E} = \frac{q}{4\pi\varepsilon_0 \varepsilon^2}\vec{e}_r，dS⃗=e⃗rdSd\vec{S} = \vec{e}_r dS ∬Σε,外E⃗⋅dS⃗=q4πε0ε2⋅4πε2=qε0\iint_{\Sigma_{\varepsilon,外}} \vec{E} \cdot d\vec{S} = \frac{q}{4\pi\varepsilon_0 \varepsilon^2} \cdot 4\pi\varepsilon^2 = \frac{q}{\varepsilon_0}
  6. 因为 Σε\Sigma_\varepsilon 取内侧，所以 ∬ΣεE⃗⋅dS⃗=−qε0\iint_{\Sigma_\varepsilon} \vec{E} \cdot d\vec{S} = -\frac{q}{\varepsilon_0}
  7. 因此 ∬ΣE⃗⋅dS⃗=−(−qε0)=qε0\iint_\Sigma \vec{E} \cdot d\vec{S} = -(-\frac{q}{\varepsilon_0}) = \frac{q}{\varepsilon_0}

答案：通过任意包围点电荷的闭曲面的电通量为 qε0\frac{q}{\varepsilon_0}，与曲面形状无关。

（这就是高斯定律的验证。）

例题2：补面法计算曲面积分

计算 ∬Σ(x3 dy dz+y3 dz dx+z3 dx dy)\iint_\Sigma (x^3 \, dy\,dz + y^3 \, dz\,dx + z^3 \, dx\,dy)，其中 Σ\Sigma 是上半球面 z=R2−x2−y2z = \sqrt{R^2-x^2-y^2} 的上侧。

参考答案(4 个标签)
高斯公式补面法上半球面球坐标

解题思路： 曲面不是闭曲面，补上底面 z=0z=0 的圆面（下侧），用高斯公式。

详细步骤：

  1. 补 Σ′:z=0, x2+y2≤R2\Sigma': z=0, \, x^2+y^2 \leq R^2，取下侧。
  2. Σ∪Σ′\Sigma \cup \Sigma' 构成闭曲面（外侧），围成上半球体 Ω\Omega。
  3. ∇⋅F⃗=3x2+3y2+3z2=3(x2+y2+z2)\nabla \cdot \vec{F} = 3x^2 + 3y^2 + 3z^2 = 3(x^2+y^2+z^2)
  4. 由高斯公式： ∯Σ∪Σ′=∭Ω3(x2+y2+z2)dv=3∫02πdθ∫0π/2sin⁡ϕ dϕ∫0Rr2⋅r2dr=3⋅2π⋅1⋅R55=6πR55\oiint_{\Sigma \cup \Sigma'} = \iiint_\Omega 3(x^2+y^2+z^2) dv = 3\int_0^{2\pi} d\theta \int_0^{\pi/2} \sin\phi \, d\phi \int_0^R r^2 \cdot r^2 dr = 3 \cdot 2\pi \cdot 1 \cdot \frac{R^5}{5} = \frac{6\pi R^5}{5}
  5. 在 Σ′\Sigma' 上（下侧）：z=0, dz=0z=0, \, dz=0，dx dydx\,dy 投影为负 ∬Σ′=∬Σ′0 dx dy=0\iint_{\Sigma'} = \iint_{\Sigma'} 0 \, dx\,dy = 0（因为 z=0z=0，z3=0z^3=0；且 dy dz=dz dx=0dy\,dz=dz\,dx=0）
  6. 所以 ∬Σ=6πR55−0=6πR55\iint_\Sigma = \frac{6\pi R^5}{5} - 0 = \frac{6\pi R^5}{5}

答案：6πR55\frac{6\pi R^5}{5}

练习题

练习1

用高斯公式计算 ∯Σx dy dz+y dz dx+z dx dy\oiint_\Sigma x \, dy\,dz + y \, dz\,dx + z \, dx\,dy，其中 Σ\Sigma 是球面 x2+y2+z2=R2x^2+y^2+z^2 = R^2 的外侧。

参考答案(4 个标签)
高斯公式球面散度体积

解题思路： 散度为常数3。

详细步骤：

  1. ∇⋅F⃗=1+1+1=3\nabla \cdot \vec{F} = 1+1+1 = 3
  2. 积分 =∭Ω3 dv=3⋅43πR3=4πR3= \iiint_\Omega 3 \, dv = 3 \cdot \frac{4}{3}\pi R^3 = 4\pi R^3

答案：4πR34\pi R^3

练习2

计算 ∬Σ(x2+y2) dx dy\iint_\Sigma (x^2+y^2) \, dx\,dy，其中 Σ\Sigma 是锥面 z=x2+y2z = \sqrt{x^2+y^2} 被平面 z=1z=1 截下部分的下侧。

参考答案(4 个标签)
高斯公式补面法锥面柱坐标

解题思路： 补上顶面 z=1z=1 的圆面（上侧），用高斯公式。

详细步骤：

  1. 补 Σ′:z=1, x2+y2≤1\Sigma': z=1, \, x^2+y^2 \leq 1，取上侧。
  2. F⃗=(0,0,x2+y2)\vec{F} = (0, 0, x^2+y^2)，∇⋅F⃗=0\nabla \cdot \vec{F} = 0
  3. 由高斯公式：∯Σ∪Σ′=∭Ω0 dv=0\oiint_{\Sigma \cup \Sigma'} = \iiint_\Omega 0 \, dv = 0
  4. 所以 ∬Σ=−∬Σ′\iint_\Sigma = -\iint_{\Sigma'}
  5. 在 Σ′\Sigma' 上（上侧）：∬Σ′(x2+y2)dx dy=∫02πdθ∫01ρ2⋅ρ dρ=2π⋅14=π2\iint_{\Sigma'} (x^2+y^2) dx\,dy = \int_0^{2\pi} d\theta \int_0^1 \rho^2 \cdot \rho \, d\rho = 2\pi \cdot \frac{1}{4} = \frac{\pi}{2}
  6. 所以 ∬Σ=−π2\iint_\Sigma = -\frac{\pi}{2}

答案：−π2-\frac{\pi}{2}

总结

本文出现的符号

符号类型读音/说明在本文中的含义
Σ\Sigma闭曲面closed surface空间区域的边界曲面
Ω\Omega空间区域spatial domain高斯公式中的空间区域
∇⋅F⃗\nabla \cdot \vec{F}散度divergence向量场的散度
E⃗\vec{E}电场强度electric field电场强度向量
B⃗\vec{B}磁感应强度magnetic field磁感应强度向量
ρ\rho电荷密度/密度charge density / density电荷体密度或流体密度
ε0\varepsilon_0真空介电常数permittivity of free space真空介电常数
Q内Q_{\text{内}}内部电荷量enclosed charge闭曲面内的总电荷量

中英对照

中文术语英文术语音标说明
高斯公式Gauss’s theorem / divergence theorem/ɡaʊsɪz ˈθɪərəm/ /daɪˈvɜːrdʒəns ˈθɪərəm/闭曲面积分与三重积分的关系
散度定理divergence theorem/daɪˈvɜːrdʒəns ˈθɪərəm/即高斯公式
奥斯特罗格拉德斯基公式Ostrogradsky’s theorem/ˌɒstroʊˈɡrædski z ˈθɪərəm/高斯公式的俄语名称
通量flux/flʌks/向量场通过曲面的量
外侧outward normal/ˈaʊtwərd ˈnɔːrməl/闭曲面的外法向一侧
补面法surface completion method/ˈsɜːrfɪs kəmˈpliːʃən ˈmɛθəd/非闭曲面补成闭曲面的方法
高斯定律Gauss’s law/ɡaʊsɪz lɔː/电磁学中的通量定律
连续性方程continuity equation/ˌkɒntɪˈnjuːəti ɪˈkweɪʒən/质量守恒的微分形式
磁单极magnetic monopole/mæɡˈnɛtɪk ˈmɒnoʊpoʊl/单一磁极（不存在）
点源point source/pɔɪnt sɔːrs/集中在一点的源
上一章节深度专题：格林公式深入探讨格林公式的历史背景、证明思路、推广应用及其在向量场理论中的核心地位。
下一章节深度专题：斯托克斯公式深入探讨斯托克斯公式的历史背景、几何意义、与格林公式和高斯公式的统一关系，以及在电磁学中的应用。
多元函数积分学