这些名字放在一起,很多人的第一反应是把高数下册和物理大地测量学教材的索引抄了一遍。我这两年帮人排查大地水准面解算代码时,见过太多次类似的翻车:拿着全球重力异常格网直接套Stokes积分,出来和已有模型差了快两米;或者把高数里“曲面积分转边界线积分”的Stokes公式,当成大地测量里用重力异常反演扰动位的Stokes积分,结果对不上号还以为是系数写错了。
如果你也正好被这组公式绕晕,那这篇文章就是写给你的。这里没有考证上的长篇大论,只有两件事:先讲清格林公式、高斯公式、Stokes公式在数学上的统一脉络,再讲清从Stokes积分、Hotine积分到Vening-Meinesz积分,它们分别解决哪类地球重力场问题、数据形式是什么、代码里最容易踩哪些坑。适合正在学物理大地测量的学生,也适合做GNSS高程转换、区域大地水准面精化的工程师。
1. 先把三个“高数公式”理顺:它们不是三兄弟,是同一个思想的三个投影
很多人把格林公式、高斯公式、Stokes公式当成三个需要分别背诵的结论,这也没错,考试时确实要分别写。但它们骨子里是同一件事:区域内部的积分,可以用边界的积分来表示。差别只在于“区域”是平面、立体还是曲面,以及被积函数是标量场还是向量场。
1.1 三个公式的字面区别
格林公式最常见的形式是平面闭曲线L和它围成的平面区域D之间建立联系:
- 格林公式:∮L P dx + Q dy = ∬D (∂Q/∂x - ∂P/∂y) dxdy
高斯公式把二维平面推广到三维空间闭区域,把空间闭曲面的通量积分和内部体积分联系起来:
- 高斯公式:∭V ∇·F dV = ∯∂V F·n dS
Stokes公式则是把空间中的有向曲面与其边界闭曲线联系起来:
- 斯托克斯公式:∮∂Σ F·dr = ∬Σ (∇×F)·n dS
从符号上很容易看花眼,但记的时候可以抓一句话:凡是“内部积分”能写成“边界积分”的,都源自微积分基本定理的高维版本。
1.2 用微积分基本定理统一理解
微积分基本定理在一维时说:函数在闭区间上的导数积分,等于两端点的函数值差。
∫ab f'(x) dx = f(b) - f(a)
这里“区间内部”是[a,b],“区间边界”就是两个端点a和b。如果把端点的值看成一种零维边界积分,这句话和格林公式、高斯公式、Sto
