格林公式

格林公式的意义

格林公式建立了平面闭区域上的二重积分与沿区域边界的第二类曲线积分之间的联系。它是微积分基本定理(牛顿-莱布尼茨公式)在二维的推广,将”区域内部的变化”与”边界上的累积”联系起来。

格林公式

定理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 的取正向的边界曲线。

正向的规定

格林公式的物理意义

格林公式可以理解为:沿闭曲线的环流量 = 区域内旋度的积分。在二维情况下,向量场 (P,Q)(P,Q) 的旋度就是 ∂Q∂x−∂P∂y\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}。

格林公式的应用

应用一:计算平面图形的面积

在格林公式中取 P=−y, Q=xP = -y, \, Q = x,则 ∂Q∂x−∂P∂y=1−(−1)=2\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} = 1 - (-1) = 2,于是:

∮Lx dy−y dx=2∬D1 dx dy=2S\oint_L x \, dy - y \, dx = 2\iint_D 1 \, dx\,dy = 2S

因此区域 DD 的面积为:

S=12∮Lx dy−y dxS = \frac{1}{2}\oint_L x \, dy - y \, dx

也可以取 P=0,Q=xP = 0, Q = x,得 S=∮Lx dyS = \oint_L x \, dy;或取 P=−y,Q=0P = -y, Q = 0,得 S=∮L−y dxS = \oint_L -y \, dx。

应用二:简化曲线积分

当直接计算曲线积分困难时,可以利用格林公式将其转化为二重积分。特别是当 ∂Q∂x−∂P∂y\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} 很简单(如常数、零)时,二重积分往往更容易计算。

应用三:复连通区域

对于复连通区域 DD(有 nn 个洞),格林公式仍然成立,但边界曲线包括外边界和所有内边界,且各自取正向:

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

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

平面曲线积分与路径无关

定理2

设区域 GG 是单连通区域,函数 P(x,y),Q(x,y)P(x,y), Q(x,y) 在 GG 内具有一阶连续偏导数,则以下四个命题等价:

  1. 曲线积分 ∫LP dx+Q dy\int_L P \, dx + Q \, dy 在 GG 内与路径无关(只依赖起点和终点)

  2. 沿 GG 内任意闭曲线 CC 的曲线积分为零:∮CP dx+Q dy=0\oint_C P \, dx + Q \, dy = 0

  3. 在 GG 内存在函数 u(x,y)u(x,y),使得 du=P dx+Q dydu = P \, dx + Q \, dy(即 P=∂u∂x, Q=∂u∂yP = \frac{\partial u}{\partial x}, \, Q = \frac{\partial u}{\partial y})

  4. 在 GG 内恒有 ∂P∂y=∂Q∂x\frac{\partial P}{\partial y} = \frac{\partial Q}{\partial x}

原函数的求法

当曲线积分与路径无关时,Pdx+QdyPdx+Qdy 是某个函数 u(x,y)u(x,y) 的全微分,u(x,y)u(x,y) 称为原函数。原函数可以通过以下公式求得:

u(x,y)=∫(x0,y0)(x,y)P dx+Q dyu(x,y) = \int_{(x_0,y_0)}^{(x,y)} P \, dx + Q \, dy

通常选择折线路径(先水平后竖直,或先竖直后水平):

u(x,y)=∫x0xP(x,y0) dx+∫y0yQ(x,y) dyu(x,y) = \int_{x_0}^x P(x, y_0) \, dx + \int_{y_0}^y Q(x, y) \, dy

或

u(x,y)=∫y0yQ(x0,y) dy+∫x0xP(x,y) dxu(x,y) = \int_{y_0}^y Q(x_0, y) \, dy + \int_{x_0}^x P(x, y) \, dx

全微分方程

如果一阶微分方程 P(x,y)dx+Q(x,y)dy=0P(x,y)dx + Q(x,y)dy = 0 的左端是某个函数 u(x,y)u(x,y) 的全微分,即 du=Pdx+Qdydu = Pdx+Qdy,则称该方程为全微分方程。其通解为 u(x,y)=Cu(x,y) = C(CC 为常数)。

判断全微分方程的条件:∂P∂y=∂Q∂x\frac{\partial P}{\partial y} = \frac{\partial Q}{\partial x}(在单连通区域内)。

典型例题

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

用格林公式求椭圆 x2a2+y2b2≤1\frac{x^2}{a^2} + \frac{y^2}{b^2} \leq 1 的面积。

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

解题思路: 利用面积公式 S=12∮Lx dy−y dxS = \frac{1}{2}\oint_L x \, dy - y \, dx。

详细步骤:

  1. 椭圆参数方程:x=acos⁡t, y=bsin⁡t, t:0→2πx = a\cos t, \, y = b\sin t, \, t: 0 \to 2\pi(逆时针)
  2. dx=−asin⁡t dt, dy=bcos⁡t dtdx = -a\sin t \, dt, \, dy = b\cos t \, dt
  3. S=12∮Lx dy−y dx=12∫02π[acos⁡t⋅bcos⁡t−bsin⁡t⋅(−asin⁡t)]dt=12∫02πab(cos⁡2t+sin⁡2t)dt=ab2⋅2π=πabS = \frac{1}{2}\oint_L x \, dy - y \, dx = \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} \cdot 2\pi = \pi ab

答案:πab\pi ab

例题2:用格林公式计算曲线积分

计算 ∮L(x2−2y) dx+(3x+yey) dy\oint_L (x^2 - 2y) \, dx + (3x + ye^y) \, dy,其中 LL 是圆周 x2+y2=4x^2+y^2=4 取逆时针方向。

参考答案
格林公式圆周二重积分逆时针

解题思路: 用格林公式转化为二重积分,∂Q∂x−∂P∂y\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} 可能很简单。

详细步骤:

  1. P=x2−2y, Q=3x+yeyP = x^2 - 2y, \, Q = 3x + ye^y
  2. ∂Q∂x=3,∂P∂y=−2\frac{\partial Q}{\partial x} = 3, \quad \frac{\partial P}{\partial y} = -2
  3. ∂Q∂x−∂P∂y=3−(−2)=5\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} = 3 - (-2) = 5
  4. 由格林公式:∮LPdx+Qdy=∬D5 dx dy=5×圆面积=5×4π=20π\oint_L Pdx+Qdy = \iint_D 5 \, dx\,dy = 5 \times \text{圆面积} = 5 \times 4\pi = 20\pi

答案:20π20\pi

例题3:曲线积分与路径无关

验证曲线积分 ∫L(x+y) dx+(x−y) dy\int_L (x+y) \, dx + (x-y) \, dy 与路径无关,并求原函数。

参考答案
格林公式路径无关原函数全微分

解题思路: 验证 ∂P∂y=∂Q∂x\frac{\partial P}{\partial y} = \frac{\partial Q}{\partial x},然后用折线路径求原函数。

详细步骤:

  1. P=x+y, Q=x−yP = x+y, \, Q = x-y
  2. ∂P∂y=1,∂Q∂x=1\frac{\partial P}{\partial y} = 1, \quad \frac{\partial Q}{\partial x} = 1
  3. ∂P∂y=∂Q∂x\frac{\partial P}{\partial y} = \frac{\partial Q}{\partial x},在全平面(单连通)成立,故积分与路径无关。
  4. 求原函数,取 (x0,y0)=(0,0)(x_0,y_0) = (0,0),折线路径: u(x,y)=∫0xP(x,0)dx+∫0yQ(x,y)dy=∫0xx dx+∫0y(x−y)dy=x22+[xy−y22]0y=x22+xy−y22u(x,y) = \int_0^x P(x,0) dx + \int_0^y Q(x,y) dy = \int_0^x x \, dx + \int_0^y (x-y) dy = \frac{x^2}{2} + [xy - \frac{y^2}{2}]_0^y = \frac{x^2}{2} + xy - \frac{y^2}{2}
  5. 验证:∂u∂x=x+y=P,∂u∂y=x−y=Q\frac{\partial u}{\partial x} = x+y = P, \quad \frac{\partial u}{\partial y} = x-y = Q ✓

答案:原函数为 u(x,y)=x22+xy−y22+Cu(x,y) = \frac{x^2}{2} + xy - \frac{y^2}{2} + C


练习题

练习1

用格林公式计算 ∮L(2xy−x2) dx+(x+y2) dy\oint_L (2xy - x^2) \, dx + (x+y^2) \, dy,其中 LL 是由 y=x2y=x^2 和 x=y2x=y^2 围成的区域的正向边界。

参考答案
格林公式抛物线围成区域二重积分

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

详细步骤:

  1. P=2xy−x2, Q=x+y2P = 2xy - x^2, \, Q = x + y^2
  2. ∂Q∂x=1,∂P∂y=2x\frac{\partial Q}{\partial x} = 1, \quad \frac{\partial P}{\partial y} = 2x
  3. ∂Q∂x−∂P∂y=1−2x\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} = 1 - 2x
  4. 区域 D:0≤x≤1, x2≤y≤xD: 0 \leq x \leq 1, \, x^2 \leq y \leq \sqrt{x}
  5. ∬D(1−2x)dx dy=∫01dx∫x2x(1−2x)dy=∫01(1−2x)(x−x2)dx\iint_D (1-2x) dx\,dy = \int_0^1 dx \int_{x^2}^{\sqrt{x}} (1-2x) dy = \int_0^1 (1-2x)(\sqrt{x}-x^2) dx =∫01(x1/2−x2−2x3/2+2x3)dx=[23x3/2−x33−45x5/2+x42]01=23−13−45+12=130= \int_0^1 (x^{1/2} - x^2 - 2x^{3/2} + 2x^3) dx = [\frac{2}{3}x^{3/2} - \frac{x^3}{3} - \frac{4}{5}x^{5/2} + \frac{x^4}{2}]_0^1 = \frac{2}{3} - \frac{1}{3} - \frac{4}{5} + \frac{1}{2} = \frac{1}{30}

答案:130\frac{1}{30}

练习2

验证 ∫L(ey+x) dx+(xey−2y) dy\int_L (e^y + x) \, dx + (xe^y - 2y) \, dy 与路径无关,并计算从 (0,0)(0,0) 到 (1,2)(1,2) 的积分值。

参考答案
路径无关原函数全微分指数函数

解题思路: 验证 ∂P∂y=∂Q∂x\frac{\partial P}{\partial y} = \frac{\partial Q}{\partial x},求原函数后代入端点。

详细步骤:

  1. P=ey+x, Q=xey−2yP = e^y + x, \, Q = xe^y - 2y
  2. ∂P∂y=ey,∂Q∂x=ey\frac{\partial P}{\partial y} = e^y, \quad \frac{\partial Q}{\partial x} = e^y,相等,故与路径无关。
  3. 原函数:u(x,y)=∫0x(1+x)dx+∫0y(xey−2y)dy=x+x22+[xey−y2]0y=x+x22+xey−y2−x=x22+xey−y2u(x,y) = \int_0^x (1+x)dx + \int_0^y (xe^y-2y)dy = x + \frac{x^2}{2} + [xe^y - y^2]_0^y = x + \frac{x^2}{2} + xe^y - y^2 - x = \frac{x^2}{2} + xe^y - y^2
  4. 从 (0,0)(0,0) 到 (1,2)(1,2):u(1,2)−u(0,0)=(12+e2−4)−0=e2−72u(1,2) - u(0,0) = (\frac{1}{2} + e^2 - 4) - 0 = e^2 - \frac{7}{2}

答案:e2−72e^2 - \frac{7}{2}


总结

本文出现的符号

符号类型读音/说明在本文中的含义
∮L\oint_L闭曲线积分closed line integral沿闭曲线的积分
DD平面区域planar region闭曲线围成的平面区域
LL边界曲线boundary curve区域 D 的边界曲线
u(x,y)u(x,y)原函数potential function满足 du=Pdx+Qdydu = Pdx+Qdy 的函数
GG单连通区域simply connected domain没有”洞”的区域

中英对照

中文术语英文术语音标说明
格林公式Green’s theorem/ɡriːnz ˈθɪərəm/平面二重积分与边界曲线积分的关系
单连通区域simply connected domain/ˈsɪmpli kəˈnɛktɪd doʊˈmeɪn/区域内任意闭曲线可连续收缩为一点
复连通区域multiply connected domain/ˈmʌltɪplaɪ kəˈnɛktɪd doʊˈmeɪn/有”洞”的区域
正向positive orientation/ˈpɒzətɪv ˌɔːriɛnˈteɪʃən/区域在左侧的边界方向
与路径无关path independent/pæθ ˌɪndɪˈpɛndənt/积分值只依赖端点
原函数potential function/pəˈtɛnʃəl ˈfʌŋkʃən/全微分的原函数
全微分方程exact differential equation/ɪɡˈzækt ˌdɪfəˈrɛnʃəl ɪˈkweɪʒən/左端为全微分的微分方程
环流量circulation/ˌsɜːrkjəˈleɪʃən/沿闭曲线的第二类曲线积分
旋度curl/kɜːrl/向量场的旋转强度,二维为 ∂Q∂x−∂P∂y\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}

发现内容有问题?

告诉我们哪里需要改进。