格林公式

格林公式的故事

一位磨坊主的数学遗产

1793年,乔治·格林(George Green)出生在英国诺丁汉郡的一个磨坊主家庭。他几乎没有受过正规教育,大部分时间在父亲的磨坊里工作。然而,这位自学成才的数学天才在1828年自费出版了一本小册子——《数学分析在电磁理论中的应用》,其中包含了后来以他名字命名的格林公式和格林函数。

格林的工作在当时几乎无人问津,直到他去世后,开尔文勋爵发现了这本小册子的价值,将其推荐给数学界。如今,格林公式已成为向量分析、偏微分方程、电磁学等领域的基石。

从”边界”到”内部”的桥梁

格林公式揭示了一个深刻的数学思想:区域内部的积分可以用区域边界上的积分来表示。这一思想贯穿了整个多元积分学——格林公式(二维)、高斯公式(三维)、斯托克斯公式(三维曲面)都是这一思想的不同表现。

格林公式的深入理解

公式的精确表述

定理1

设闭区域 DD 由分段光滑的曲线 LL 围成,函数 P(x,y)P(x,y) 及 Q(x,y)Q(x,y) 在 DD 上具有一阶连续偏导数,则有

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

其中 LL 是 DD 的取正向的边界曲线。

正向的规定

两种形式的格林公式

格林公式有两种等价形式:

  1. 旋度形式(上述公式):∮LP dx+Q dy=∬D(∂Q∂x−∂P∂y)dx dy\oint_L P \, dx + Q \, dy = \iint_D (\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}) dx\,dy

  2. 散度形式(通量形式): ∮L(Pcos⁡α+Qcos⁡β)ds=∬D(∂P∂x+∂Q∂y)dx dy\oint_L (P\cos\alpha + Q\cos\beta) ds = \iint_D \left( \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} \right) dx\,dy 其中 (cos⁡α,cos⁡β)(\cos\alpha, \cos\beta) 是边界曲线的外法向量。这是高斯公式在二维的对应。

证明思路

格林公式的证明可以分为三步:

第一步:X型区域

设 DD 是 XX 型区域:a≤x≤b, φ1(x)≤y≤φ2(x)a \leq x \leq b, \, \varphi_1(x) \leq y \leq \varphi_2(x)。

∬D∂P∂ydx dy=∫abdx∫φ1(x)φ2(x)∂P∂ydy=∫ab[P(x,φ2(x))−P(x,φ1(x))]dx\iint_D \frac{\partial P}{\partial y} dx\,dy = \int_a^b dx \int_{\varphi_1(x)}^{\varphi_2(x)} \frac{\partial P}{\partial y} dy = \int_a^b [P(x,\varphi_2(x)) - P(x,\varphi_1(x))] dx

而沿边界的曲线积分: ∮LP dx=∫L1P dx+∫L2P dx=∫abP(x,φ1(x))dx−∫abP(x,φ2(x))dx\oint_L P \, dx = \int_{L_1} P \, dx + \int_{L_2} P \, dx = \int_a^b P(x,\varphi_1(x)) dx - \int_a^b P(x,\varphi_2(x)) dx (注意 L2L_2 的方向是从 bb 到 aa)

因此 −∬D∂P∂ydx dy=∮LP dx-\iint_D \frac{\partial P}{\partial y} dx\,dy = \oint_L P \, dx。

第二步:Y型区域

类似地,对 YY 型区域证明 ∬D∂Q∂xdx dy=∮LQ dy\iint_D \frac{\partial Q}{\partial x} dx\,dy = \oint_L Q \, dy。

第三步:一般区域

对于一般区域,可以将其分割为若干个 XX 型或 YY 型区域,在每个子区域上应用上述结果,相邻子区域的公共边界上的曲线积分相互抵消(方向相反),最终得到整个区域上的格林公式。

格林公式的推广

多连通区域

格林公式可以推广到多连通区域(有 nn 个”洞”的区域):

∮L0P dx+Q dy+∑i=1n∮LiP dx+Q dy=∬D(∂Q∂x−∂P∂y)dx dy\oint_{L_0} P \, dx + Q \, dy + \sum_{i=1}^n \oint_{L_i} P \, dx + Q \, dy = \iint_D \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) dx\,dy

其中 L0L_0 是外边界(逆时针),LiL_i 是内边界(顺时针)。

格林第一公式和第二公式

利用格林公式可以推导出偏微分方程中极为重要的格林公式:

格林第一公式: ∬Du∇2v dx dy=∮Lu∂v∂nds−∬D∇u⋅∇v dx dy\iint_D u \nabla^2 v \, dx\,dy = \oint_L u \frac{\partial v}{\partial n} ds - \iint_D \nabla u \cdot \nabla v \, dx\,dy

格林第二公式: ∬D(u∇2v−v∇2u)dx dy=∮L(u∂v∂n−v∂u∂n)ds\iint_D (u \nabla^2 v - v \nabla^2 u) dx\,dy = \oint_L \left( u \frac{\partial v}{\partial n} - v \frac{\partial u}{\partial n} \right) ds

其中 ∂∂n\frac{\partial}{\partial n} 是沿外法向的方向导数。这两个公式是研究拉普拉斯方程、泊松方程的基本工具。

格林公式的应用

1. 计算平面面积

取 P=−y, Q=xP = -y, \, Q = x,得: S=12∮Lx dy−y dxS = \frac{1}{2}\oint_L x \, dy - y \, dx

2. 计算曲线积分

当直接计算曲线积分困难时,可以用格林公式转化为二重积分。特别是当 ∂Q∂x−∂P∂y\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} 很简单时(如常数),效果显著。

3. 路径无关与全微分方程

格林公式是曲线积分与路径无关理论的基础。在单连通区域内,∂Q∂x=∂P∂y\frac{\partial Q}{\partial x} = \frac{\partial P}{\partial y} 等价于:

  • 曲线积分与路径无关
  • 沿任意闭曲线积分为零
  • Pdx+QdyPdx + Qdy 是某个函数的全微分
  • 向量场 (P,Q)(P,Q) 是有势场

4. 复变函数中的柯西定理

格林公式是复变函数中柯西积分定理的实分析基础。

典型例题

例题1:用格林公式计算面积

用格林公式计算椭圆 x2a2+y2b2=1\frac{x^2}{a^2} + \frac{y^2}{b^2} = 1 所围区域的面积。

参考答案
格林公式面积椭圆参数方程

解题思路: 用面积公式 S=12∮Lx dy−y dxS = \frac{1}{2}\oint_L x\,dy - y\,dx,椭圆参数方程 x=acos⁡t,y=bsin⁡tx=a\cos t, y=b\sin t。

详细步骤:

  1. dx=−asin⁡t dt, dy=bcos⁡t dtdx = -a\sin t \, dt, \, dy = b\cos t \, dt
  2. S=12∫02π[acos⁡t⋅bcos⁡t−bsin⁡t⋅(−asin⁡t)]dt=12∫02πab(cos⁡2t+sin⁡2t)dt=ab2∫02πdt=πabS = \frac{1}{2}\int_0^{2\pi} [a\cos t \cdot b\cos t - b\sin t \cdot (-a\sin t)] dt = \frac{1}{2}\int_0^{2\pi} ab(\cos^2 t + \sin^2 t) dt = \frac{ab}{2}\int_0^{2\pi} dt = \pi ab

答案:πab\pi ab

例题2:多连通区域的格林公式

计算 ∮Lx dy−y dxx2+y2\oint_L \frac{x \, dy - y \, dx}{x^2 + y^2},其中 LL 是不经过原点的任意正向简单闭曲线。

参考答案
格林公式多连通区域奇点挖洞法

解题思路: 被积函数在原点不连续,需分情况讨论。当 LL 不包围原点时直接用格林公式;当 LL 包围原点时用”挖洞法”。

详细步骤:

  1. P=−yx2+y2, Q=xx2+y2P = \frac{-y}{x^2+y^2}, \, Q = \frac{x}{x^2+y^2}
  2. ∂Q∂x=y2−x2(x2+y2)2=∂P∂y\frac{\partial Q}{\partial x} = \frac{y^2-x^2}{(x^2+y^2)^2} = \frac{\partial P}{\partial y}(在 x2+y2≠0x^2+y^2 \neq 0 处)
  3. 情况一:LL 不包围原点。在 LL 所围区域内 P,QP,Q 有连续偏导数,由格林公式,积分 =∬D0 dx dy=0= \iint_D 0 \, dx\,dy = 0。
  4. 情况二:LL 包围原点。在 LL 内作小圆周 l:x2+y2=ε2l: x^2+y^2 = \varepsilon^2(顺时针),在 LL 和 ll 之间的多连通区域上用格林公式: ∮L+∮l=∬D0=0\oint_L + \oint_l = \iint_D 0 = 0,所以 ∮L=−∮l\oint_L = -\oint_l
  5. 在 ll 上(逆时针):∮l逆x dy−y dxx2+y2=1ε2∮l逆x dy−y dx=1ε2⋅2πε2=2π\oint_{l_{逆}} \frac{x\,dy-y\,dx}{x^2+y^2} = \frac{1}{\varepsilon^2}\oint_{l_{逆}} x\,dy-y\,dx = \frac{1}{\varepsilon^2} \cdot 2\pi\varepsilon^2 = 2\pi
  6. 所以 ∮L=2π\oint_L = 2\pi

答案:当 LL 不包围原点时为 00;当 LL 包围原点时为 2π2\pi。

(这是复变函数中”绕数”的实分析版本。)


练习题

练习1

用格林公式计算 ∮L(x+y) dx−(x−y) dy\oint_L (x+y) \, dx - (x-y) \, dy,其中 LL 是椭圆 x2a2+y2b2=1\frac{x^2}{a^2}+\frac{y^2}{b^2}=1 的正向。

参考答案
格林公式椭圆曲线积分

解题思路: 用格林公式转化为二重积分。

详细步骤:

  1. P=x+y, Q=−(x−y)=y−xP = x+y, \, Q = -(x-y) = y-x
  2. ∂Q∂x−∂P∂y=−1−1=−2\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} = -1 - 1 = -2
  3. 积分 =∬D(−2)dx dy=−2⋅πab=−2πab= \iint_D (-2) dx\,dy = -2 \cdot \pi ab = -2\pi ab

答案:−2πab-2\pi ab

练习2

计算 ∮Ly dx−x dyx2+y2\oint_L \frac{y \, dx - x \, dy}{x^2+y^2},其中 LL 是圆周 x2+y2=a2x^2+y^2 = a^2 的正向。

参考答案
格林公式圆周奇点参数方程

解题思路: 在圆周上 x2+y2=a2x^2+y^2=a^2,直接用参数方程计算。

详细步骤:

  1. x=acos⁡t, y=asin⁡t, t:0→2πx = a\cos t, \, y = a\sin t, \, t: 0 \to 2\pi
  2. dx=−asin⁡t dt, dy=acos⁡t dtdx = -a\sin t \, dt, \, dy = a\cos t \, dt
  3. 积分 =∫02πasin⁡t(−asin⁡t)−acos⁡t(acos⁡t)a2dt=∫02π−a2(sin⁡2t+cos⁡2t)a2dt=∫02π(−1)dt=−2π= \int_0^{2\pi} \frac{a\sin t(-a\sin t) - a\cos t(a\cos t)}{a^2} dt = \int_0^{2\pi} \frac{-a^2(\sin^2 t+\cos^2 t)}{a^2} dt = \int_0^{2\pi} (-1) dt = -2\pi

答案:−2π-2\pi


总结

本文出现的符号

符号类型读音/说明在本文中的含义
LL边界曲线boundary curve区域 D 的正向边界曲线
DD平面区域planar domain格林公式中的平面区域
∂∂n\frac{\partial}{\partial n}法向导数normal derivative沿外法向的方向导数
∇2\nabla^2拉普拉斯算子Laplacian∂2∂x2+∂2∂y2\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}
∮L\oint_L闭曲线积分closed line integral沿闭曲线的积分

中英对照

中文术语英文术语音标说明
格林公式Green’s theorem/ɡriːnz ˈθɪərəm/平面曲线积分与二重积分的关系
格林函数Green’s function/ɡriːnz ˈfʌŋkʃən/偏微分方程的基本解
格林第一公式Green’s first identity/ɡriːnz fɜːrst aɪˈdɛntəti/含拉普拉斯算子的积分恒等式
格林第二公式Green’s second identity/ɡriːnz sɛkənd aɪˈdɛntəti/两个函数的格林公式
正向positive orientation/ˈpɒzətɪv ˌɔːriənˈteɪʃən/区域在左侧的边界方向
多连通区域multiply connected domain/ˈmʌltɪplaɪ kəˈnɛktɪd doʊˈmeɪn/有”洞”的区域
单连通区域simply connected domain/ˈsɪmpli kəˈnɛktɪd doʊˈmeɪn/没有”洞”的区域
挖洞法hole-digging method/hoʊl ˈdɪɡɪŋ ˈmɛθəd/处理奇点的方法
绕数winding number/ˈwaɪndɪŋ ˈnʌmbər/曲线绕点的次数
柯西定理Cauchy’s theorem/ˈkoʊʃiz ˈθɪərəm/复变函数的积分定理

发现内容有问题?

告诉我们哪里需要改进。