跳转到内容
← 返回定理
向量分析中等约 19 分钟阅读

格林定理

Green's Theorem

green·1828
向量分析线积分面积分格林

定理陈述

格林定理（Green's Theorem）建立了平面区域上的二重积分与边界曲线积分之间的关系。

设 $D$ 是平面上的有界闭区域，边界 ∂D\partial D 分段光滑且取正向（逆时针方向）。设 $P(x,y)$ 和 $Q(x,y)$ 在包含 $D$ 的开集上连续可微，则

∮∂D(P dx+Q dy)=∬D(∂Q∂x−∂P∂y)dA\oint_{\partial D} (P\,dx + Q\,dy) = \iint_D \left(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\right) dA

两种特殊形式：

面积公式：取 $P = -y/2$，$Q = x/2$，则

Area(D)=12∮∂D(x dy−y dx)=∮∂Dx dy\text{Area}(D) = \frac{1}{2}\oint_{\partial D} (x\,dy - y\,dx) = \oint_{\partial D} x\,dy

直觉理解

格林定理说的是：沿边界的环流量等于内部旋度的总和。想象一个旋转的流体。如果你沿着边界走一圈，感受到的净推力（环量）等于内部所有微小旋转（旋度）的叠加。从面积的角度看：你可以通过测量边界的环量来计算内部的面积。这就是为什么格林定理在计算面积时非常有用——将面积分转化为线积分。

格林定理是微积分基本定理在二维的推广。FTC 说 ∫abf′(x)dx=f(b)−f(a)\int_a^b f'(x)dx = f(b) - f(a)——"内部的导数 = 边界上的值"。格林定理说的也是同样的话，只是维度提高了。

证明思路

证明概要

步骤 1：假设 $D$ 是 $x$-型区域（即 D={(x,y):a≤x≤b, g1(x)≤y≤g2(x)}D = \{(x,y) : a \leq x \leq b,\, g_1(x) \leq y \leq g_2(x)\}）。

步骤 2：计算 ∬D∂Q∂xdA\iint_D \frac{\partial Q}{\partial x} dA：

∬D∂Q∂xdA=∫ab∫g1(x)g2(x)∂Q∂xdy dx=∫ab[Q(x,g2(x))−Q(x,g1(x))]dx\iint_D \frac{\partial Q}{\partial x} dA = \int_a^b \int_{g_1(x)}^{g_2(x)} \frac{\partial Q}{\partial x} dy\,dx = \int_a^b [Q(x, g_2(x)) - Q(x, g_1(x))] dx

步骤 3：计算 ∮∂DQ dy\oint_{\partial D} Q\,dy。边界由下曲线 y=g1(x)y = g_1(x)（从左到右）和上曲线 y=g2(x)y = g_2(x)（从右到左）组成：

∮∂DQ dy=∫abQ(x,g2(x))dx−∫abQ(x,g1(x))dx\oint_{\partial D} Q\,dy = \int_a^b Q(x, g_2(x))dx - \int_a^b Q(x, g_1(x))dx

两式相等。类似地证明 $P$ 的部分。

步骤 4：一般区域可以分解为有限个 $x$-型和 $y$-型子区域，内部边界上的积分相互抵消。

历史背景

乔治·格林（George Green，1793—1841）是英国自学成才的数学家和物理学家——他白天经营一座磨坊，几乎没有受过正规学校教育，被称为"诺丁汉自学成才的磨坊主"。

他在 1828 年自费出版了论文《论数学分析在电磁学理论中的应用》，其中提出了格林函数，以及把体积分化为边界面积分的"格林恒等式"。这些是势论（potential theory，研究引力场、静电场等位势的数学分支）的核心工具。格林的论文在 1828 年发表时几乎无人注意——订阅者只有五十余位，多为诺丁汉的亲友。直到 1845 年初，刚从剑桥毕业的威廉·汤姆森（William Thomson，即后来的开尔文勋爵）经导师威廉·霍普金斯之手拿到三份副本，立刻意识到它对势论的价值。

当年夏天，汤姆森把论文带到巴黎，展示给刘维尔（Liouville）与斯图姆（Sturm）；随后它在 1850—1854 年间分三部分重刊于《克雷尔杂志》（Crelle's Journal）。至此，格林的工作才被广泛认可。格林函数后来成为数学物理中最重要的工具之一——用于求解偏微分方程。值得一提的是，我们今天课本里那个"平面环量 = 旋度二重积分"的格林定理，其实并未以这个形式出现在格林 1828 年的论文中。它的第一个印刷陈述出自柯西（Augustin-Louis Cauchy）1846 年的一篇论文，而第一个严格证明出自黎曼（Bernhard Riemann）1851 年的博士论文。

把这条平面定理冠以"格林"之名，是数学史上"以领袖人物命名"的惯例之一：格林真正的贡献是更一般的势论恒等式，平面特例的功劳更多属于柯西与黎曼。这也解释了为什么格林定理与复分析中的柯西积分定理、柯西-黎曼方程如此密不可分。

应用

  1. 面积计算：通过边界曲线的线积分计算平面区域的面积
  2. 流体力学：计算流体绕障碍物的环量和升力
  3. 电磁学：安培环路定理的二维形式
  4. 复分析：柯西积分定理的实变量版本——格林定理是复分析的基础
  5. 有限元方法：格林公式用于弱形式的推导

与其他定理的关系

  • 微积分基本定理：格林定理是 FTC 在二维的推广
  • 斯托克斯定理：格林定理是斯托克斯定理在平面上的特例
  • 散度定理：格林定理的散度形式——∮F⋅n ds=∬∇⋅F dA\oint \mathbf{F} \cdot \mathbf{n}\,ds = \iint \nabla \cdot \mathbf{F}\,dA
  • 柯西积分定理：复分析中 ∮f(z)dz=0\oint f(z)dz = 0——由格林定理推出

具体示例

设 $P = -y$，$Q = x$，$D$ 是单位圆盘 x2+y2≤1x^2 + y^2 \leq 1。格林定理给出：

∮∂D(−y dx+x dy)=∬D(1−(−1)) dA=2⋅π=2π\oint_{\partial D} (-y\,dx + x\,dy) = \iint_D (1 - (-1))\,dA = 2 \cdot \pi = 2\pi

左端是沿单位圆的线积分。参数化 x=cos⁡tx = \cos t，y=sin⁡ty = \sin t：

∫02π(−sin⁡t)(−sin⁡t) dt+(cos⁡t)(cos⁡t) dt=∫02π1 dt=2π\int_0^{2\pi} (-\sin t)(-\sin t)\,dt + (\cos t)(\cos t)\,dt = \int_0^{2\pi} 1\,dt = 2\pi

两边相等。利用面积公式 Area=12∮(x dy−y dx)=π\text{Area} = \frac{1}{2}\oint(x\,dy - y\,dx) = \pi。

常见误解与适用条件

很多人背下了公式，却跳过了那行前提："$P,Q$ 在包含 $D$ 的开集上连续可微"。这个前提不是装饰，去掉它结论就会出错。考虑涡旋场（vortex field，绕一点旋转的向量场）

F=(P,Q)=(−yx2+y2, xx2+y2)\mathbf{F} = (P, Q) = \left(\frac{-y}{x^2+y^2},\ \frac{x}{x^2+y^2}\right)

直接算它的旋度：

∂Q∂x−∂P∂y=y2−x2(x2+y2)2−y2−x2(x2+y2)2=0\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} = \frac{y^2 - x^2}{(x^2+y^2)^2} - \frac{y^2 - x^2}{(x^2+y^2)^2} = 0

旋度处处为零。一个粗心的"应用"会说：沿任意闭曲线的环量都等于 ∬D0 dA=0\iint_D 0\,dA = 0。但沿单位圆实际算一遍（x=cos⁡t, y=sin⁡tx = \cos t,\ y = \sin t）：

∮∂D(P dx+Q dy)=∫02π1 dt=2π≠0\oint_{\partial D} (P\,dx + Q\,dy) = \int_0^{2\pi} 1\,dt = 2\pi \neq 0

环量是 2π2\pi，不是 $0$。矛盾出在哪里？F\mathbf{F} 在原点处没有定义（分母为零），而原点恰好落在单位圆内部。区域里藏着一个"洞"，定理的连续可微前提在洞上被破坏，结论自然失效。这是最朴素也最深刻的一课：格林定理把环量算成"内部旋度之和"，前提是内部真的处处有定义。只要中间藏着奇点，环量就可能凭空多出来。

另一个常被忽略的点是取向（orientation）。"正向"的准确含义是：沿边界行走时，区域始终在你的左手边。对带洞的多连通区域，外边界逆时针、内边界顺时针，定理依然成立——这也是绕开奇点后重新让定理适用的标准办法。

进一步看，那个 2π2\pi 并非偶然。对任意一条绕原点 $n$ 圈的闭曲线，环量恰好是 2πn2\pi n，这里的整数 $n$ 就是环绕数（winding number，曲线绕原点转过的圈数）。环量在这里测量的是一个纯粹的拓扑量，与曲线的具体形状无关——这正是德拉姆上同调（de Rham cohomology，用微分形式刻画空间"洞"的理论）在最简单情形下的样子：挖去一点的平面有一个一维的"洞"，涡旋场就是探测这个洞的工具。

格林定理与复分析

格林定理在复分析中有重要应用。设 $f(z) = u + iv$ 是区域 $D$ 上的解析函数，∮∂Df(z) dz\oint_{\partial D} f(z)\,dz 可以分解为实部和虚部的线积分。利用柯西-黎曼方程 ∂u/∂x=∂v/∂y\partial u/\partial x = \partial v/\partial y 和 ∂u/∂y=−∂v/∂x\partial u/\partial y = -\partial v/\partial x，应用格林定理可得：

∮∂Df(z) dz=0\oint_{\partial D} f(z)\,dz = 0

这就是柯西积分定理——复分析的基石。格林定理为柯西积分定理提供了实变量的证明。更一般地，柯西积分公式 f(a)=12πi∮γf(z)z−adzf(a) = \frac{1}{2\pi i}\oint_\gamma \frac{f(z)}{z-a}dz 的证明也依赖于格林定理（通过将解析函数的积分化为实部和虚部的线积分）。

格林定理的物理应用

在流体力学中，格林定理的旋度形式给出：

∮∂Dv⋅dr=∬D(∇×v)⋅k dA\oint_{\partial D} \mathbf{v} \cdot d\mathbf{r} = \iint_D (\nabla \times \mathbf{v}) \cdot \mathbf{k}\,dA

即：沿边界的速度环量等于内部涡度的总和。这是开尔文环量定理的二维版本。在电磁学中，安培环路定理的二维形式 ∮B⋅dl=μ0I\oint \mathbf{B} \cdot d\mathbf{l} = \mu_0 I 本质上是格林定理的应用。

求积仪：把面积"走"出来

格林定理的面积公式 Area=∮∂Dx dy\text{Area} = \oint_{\partial D} x\,dy 不只是纸面技巧，它被造成了真实的仪器。

1854 年，瑞士数学家雅各布·阿姆斯勒（Jacob Amsler，1823—1912）发明了极式求积仪（polar planimeter）。你只需用仪器的描迹臂沿图形边界走一圈，一个测量轮就累计转动量，读数直接给出图形的面积。

它的工作原理正是格林定理：测量轮机械地积分了边界线积分 ∮x dy\oint x\,dy，把"沿边界的环绕"自动换算成"内部的面积"。无需把图形切成小方格去数，只要沿边界走一圈。阿姆斯勒把早期复杂的机构简化成两根用枢轴相连的臂，便宜又可靠。他在 1857 年开设工坊，此后六十年间在全世界卖出约五万台，被工程、地质、生物乃至医学领域广泛用来测量不规则图形（如地图上的湖泊、显微切片里的细胞）的面积。

格林第一与第二恒等式

格林定理可以推导出两个重要的恒等式。设 $u$ 和 $v$ 是区域 $D$ 上的二阶连续可微函数：

格林第一恒等式：

∬D(u∇2v+∇u⋅∇v) dA=∮∂Du∂v∂n ds\iint_D (u\nabla^2 v + \nabla u \cdot \nabla v)\,dA = \oint_{\partial D} u\frac{\partial v}{\partial n}\,ds

其中 ∂v∂n=∇v⋅n\frac{\partial v}{\partial n} = \nabla v \cdot \mathbf{n} 是法向导数。

格林第二恒等式（将第一恒等式中的 $u, v$ 互换后相减）：

∬D(u∇2v−v∇2u) dA=∮∂D(u∂v∂n−v∂u∂n)ds\iint_D (u\nabla^2 v - v\nabla^2 u)\,dA = \oint_{\partial D}\left(u\frac{\partial v}{\partial n} - v\frac{\partial u}{\partial n}\right)ds

这两个恒等式是求解拉普拉斯方程和泊松方程的基础。格林函数法正是利用第二恒等式，将偏微分方程的解表示为边界积分——这是边界元方法（BEM）的数学基础。

格林定理的离散化

在计算数学中，格林定理的离散化形式是有限体积法的核心。将区域划分为控制体积 ViV_i，对每个体积应用散度定理的二维形式：

∬Vi∇⋅F dA=∑edgesF⋅n Δs\iint_{V_i} \nabla \cdot \mathbf{F}\,dA = \sum_{\text{edges}} \mathbf{F} \cdot \mathbf{n}\,\Delta s

这种方法在计算流体力学（CFD）中被广泛使用，因为它天然地保持了物理守恒律——离散化后的方程仍然精确满足质量、动量和能量守恒。

跨域连接

  • 同调与上同调：绕奇点一周的环量等于一个固定常数乘环绕数，是只认"洞"的整数。推论是：把曲线连续变形而不越过奇点，环量丝毫不变——它测的是拓扑而非几何，所以可以拿来当洞的探测器。反过来说，若某区域上所有闭曲线的环量都为零，就能断定它没有被挖掉的点。
  • 遥感与地理信息系统：多边形面积用顶点坐标的鞋带公式算，正是边界线积分的离散化。推论是：只需边界顶点、不必填充内部，代价随顶点数而非面积增长；顶点顺序反向会给出负面积，这一点被直接用来标记多边形里的洞。
  • 计算几何：判断点是否在多边形内，可以数射线交点的奇偶，也可以算环绕数，两者在自相交多边形上给出不同答案。推论是：填充规则必须显式选择（奇偶还是非零环绕），这不是实现细节，而是定义问题。
  • 弹性力学：有限元的弱形式由分部积分得到，把二阶导数中的一阶转移到检验函数上。推论是：位移解只需一阶可微，于是分片线性单元就能处理带角点和材料界面的真实零件——光滑性要求降一级，正是有限元能算实物的关键。
  • 病理学：切片上不规则区域的面积，传统上靠沿边界描迹换算得到。推论是：误差只来自边界的定位，与区域内部的纹理无关——这也解释了为什么分割算法的评价对边界最敏感，内部像素几乎不贡献差异。另一个后果是：染色不均只影响边界判读，不改变面积的换算关系。

延伸阅读

  1. George Green, "An Essay on the Application of Mathematical Analysis to the Theories of Electricity and Magnetism" (1828).
  2. H.M. Schey, Div, Grad, Curl, and All That (4th ed., 2005).
  3. John B. Conway, Functions of One Complex Variable (2nd ed., 1978).
  4. 吉米多维奇, 《数学分析习题集》, 人民教育出版社.
  5. R.V. Churchill & J.W. Brown, Complex Variables and Applications (9th ed., 2013).

参考文献

  1. Victor J. Katz, "The History of Stokes' Theorem," Mathematics Magazine, 52(3), 1979, pp. 146–156（JSTOR 2690275）——格林、高斯、斯托克斯三定理的史料考订与归属问题（含柯西 1846、黎曼 1851）。
  2. MacTutor History of Mathematics Archive, "George Green (1793–1841)," University of St Andrews——格林生平、磨坊出身，以及汤姆森重新发现论文的经过。
  3. Oliver Knill, "Green's Theorem," lecture notes, Harvard University Department of Mathematics——定理陈述、涡旋场反例与环绕数。
  4. "Mathematical Treasure: Polar Planimeter Invented by Jacob Amsler," MAA Convergence, Mathematical Association of America——阿姆斯勒极式求积仪与格林定理面积公式的工程实现。

格林定理把平面区域 $D$ 上的二重积分与其边界曲线 ∂D\partial D 上的线积分联系起来：∮∂D(P dx+Q dy)=∬D(∂Q∂x−∂P∂y)dA\oint_{\partial D}(P\,dx+Q\,dy)=\iint_D\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)dA。它是斯托克斯定理在平面上的特例，被广泛用于计算平面区域面积、流体环量与电磁学中的环路积分。