格林公式(Green’s Theorem)是平面上的斯托克斯定理,核心思想是将闭曲线上的线积分转化为该曲线所围区域上的二重积分。
一、标准写法#
设 C 是平面上的分段光滑、简单闭曲线(正向为逆时针),D 是由 C 所围成的有界闭区域。函数 P(x,y) 和 Q(x,y) 在 D 上具有一阶连续偏导数,则:
∮CPdx+Qdy=∬D(∂x∂Q−∂y∂P)dxdy
二、等价的表达形式#
1. 向量形式(环量/旋度形式)#
设 F=(P,Q),T 为 C 的单位切向量,k 为 z 轴单位向量:
∮CF⋅Tds=∮CF⋅dr=∬D(∇×F)⋅kdA
其中 (∇×F)⋅k=∂x∂Q−∂y∂P 是 F 的旋度在 z 方向的分量。
2. 通量/散度形式(法向形式)#
设 n 为 C 的单位外法向量:
∮CF⋅nds=∬D∇⋅FdA=∬D(∂x∂P+∂y∂Q)dA
注意:此形式通常等价于对 (−Q,P) 应用标准格林公式。
3. 微分形式(斯托克斯定理的特例)#
∮∂Dω=∬Ddω
其中 ω=Pdx+Qdy,dω=(∂x∂Q−∂y∂P)dx∧dy。
三、常用的简化版与近似版#
1. 面积公式(最重要的推论)#
取 P=−2y,Q=2x,则 ∂x∂Q−∂y∂P=1,于是:
A=∬Ddxdy=21∮Cxdy−ydx
前提:C 是包围区域 D 的简单闭曲线。这是计算不规则区域面积的经典方法。
2. 保守场(积分与路径无关)#
若在单连通区域 D 内恒有:
∂x∂Q=∂y∂P
则对 D 内任意闭曲线 C:
∮CPdx+Qdy=0
前提:区域必须是单连通的(无洞)。如果区域有洞,即使旋度为零,沿包围洞的闭曲线积分也可能非零。
3. 复连通区域(带洞的区域)#
若 D 不是单连通的,边界由外边界 Cout(逆时针)和 n 条内边界 Cin(k)(顺时针,即洞的边界反向)组成:
∮CoutPdx+Qdy+∑k=1n∮Cin(k)Pdx+Qdy=∬D(∂x∂Q−∂y∂P)dA
前提:所有边界方向需协调(外逆内顺),且 P,Q 在 D 内(包括洞之间的区域)连续可微。
4. 中值定理型近似(小区域近似)#
若区域 D 很小,由积分中值定理,存在 (ξ,η)∈D 使得:
∮CPdx+Qdy≈(∂x∂Q−∂y∂P)(ξ,η)⋅Area(D)
前提:D 足够小,使得被积函数在区域内变化平缓。这在数值计算和物理场的局部分析中常用。
5. 矩形区域上的简化验证#
若 D=[a,b]×[c,d] 为矩形,可直接通过两次牛顿-莱布尼兹公式验证:
∬D∂x∂Qdxdy=∫cd[Q(b,y)−Q(a,y)]dy
∬D∂y∂Pdxdy=∫ab[P(x,d)−P(x,c)]dx
这既是证明格林公式的常用思路,也是数值离散化的基础。
总结:标准格林公式要求闭曲线 + 单连通 + 函数连续可微。最实用的简化是面积公式和保守场判据;遇到有洞区域时,必须使用复连通版本并注意边界的定向。
下面按标准格林公式逐符号拆解,并给出不同物理场景下的单位对照。
一、符号总表#
| 符号 | 数学含义 | 物理含义(以力场为例) | 量纲/单位 | 常量/变量 |
|---|
| x,y | 平面直角坐标 | 空间位置坐标 | 长度:m(米) | 自变量(变量) |
| dx,dy | 坐标微分 | 位移在坐标轴上的投影 | 长度:m | 变量(无穷小量) |
| P(x,y) | 被积函数的第一个分量 | 力场在 x 方向的分量 Fx | 取决于场景(见下文) | 场函数(变量) |
| Q(x,y) | 被积函数的第二个分量 | 力场在 y 方向的分量 Fy | 同 P | 场函数(变量) |
| C | 分段光滑的简单闭曲线 | 积分路径/边界回路 | 几何对象(无单位) | 通常固定已知 |
| ∮C | 沿闭曲线的线积分 | 沿闭合路径累积 | 由被积式决定 | 算符 |
| ds | 弧长微元 | 路径上的微小长度 | 长度:m | 变量 |
| T | 单位切向量 | 路径方向 | 无量纲(方向) | 沿 C 变化 |
| n | 单位外法向量 | 边界外向法线方向 | 无量纲(方向) | 沿 C 变化 |
| D | C 所围成的有界闭区域 | 积分区域/控制体截面 | 几何对象(无单位) | 通常固定已知 |
| ∬D | 在 D 上的二重积分 | 区域上的累积 | 由被积式决定 | 算符 |
| dA 或 dxdy | 面积微元 | 微小面积 | 面积:m2(平方米) | 变量 |
| ∂x∂Q | Q 对 x 的偏导数 | 场在 x 方向的变化率 | [Q]/[x] | 场函数的导出量 |
| ∂y∂P | P 对 y 的偏导数 | 场在 y 方向的变化率 | [P]/[y] | 场函数的导出量 |
| ∇×F | 向量场 F 的旋度 | 场的旋转强度 | [F]/[L] | 导出场 |
| ∇⋅F | 向量场 F 的散度 | 场的源汇强度 | [F]/[L] | 导出场 |
| k | z 轴单位向量 | 垂直于平面的方向 | 无量纲 | 常量 |
二、核心物理场景与单位对照#
格林公式是数学恒等式,本身无量纲;单位完全由 P,Q 的物理意义决定。
场景 A:力学(功与能)#
若 F=(P,Q) 是力场(如 F=(Fx,Fy)),则:
| 项 | 物理含义 | SI 单位 | 常用单位 |
|---|
| P,Q | 力的分量 | N(牛顿) | kgf、dyn |
| Pdx | 元功 | J(焦耳 = N⋅m) | erg、kgf⋅m |
| ∮CPdx+Qdy | 力沿闭路做的总功(环量) | J | — |
| ∂x∂Q−∂y∂P | 力场的”旋度密度” | N/m | — |
| 积分结果 | 总功 = 旋度在区域上的累积 | J | — |
场景 B:流体力学(平面流动)#
若 F=(P,Q) 是速度场 v=(vx,vy):
| 项 | 物理含义 | SI 单位 | 常用单位 |
|---|
| P,Q | 流速分量 | m/s | cm/s、km/h |
| ∮Cv⋅dr | 速度环量 Γ | m2/s | Stokes(St) |
| ∂x∂Q−∂y∂P | 涡量(Vorticity)ωz | s−1(1/s) | — |
| 通量形式 ∮Cv⋅nds | 体积流量(穿过边界的流量) | m2/s | L/s(若乘深度则为 m3/s) |
| ∇⋅v | 源汇强度 | s−1 | — |
场景 C:电磁学#
若 F 代表电场或磁场:
| 应用 | P,Q 代表 | 线积分 | 旋度项 | 单位 |
|---|
| 静电场 | E=(Ex,Ey) | 电动势(环量) | ∇×E | V(伏特)、V/m |
| 稳恒磁场 | H=(Hx,Hy) | 磁动势 | ∇×H=Jz | A(安培)、A/m |
场景 D:纯数学(面积公式)#
取 P=−y/2,Q=x/2:
| 符号 | 含义 | 单位 | 说明 |
|---|
| P=−y/2 | 坐标函数 | m | 此时 P 不是力,只是坐标 |
| Q=x/2 | 坐标函数 | m | 同上 |
| 被积式 | 面积微元的参数化 | m2 | 结果直接是面积 |
| 积分结果 | 区域 D 的面积 | m2 | 无物理量纲,纯几何 |
三、常量、变量与已知/待求关系#
1. 常量(在单次问题中不变)#
- 坐标轴方向、k:永远固定。
- 若 P,Q 是具体给定的场(如 P=x2y,Q=xy2),则表达式中的系数是常量。
- 区域 D 和边界 C 的形状:在具体问题中通常是给定的几何条件。
2. 变量(随位置或路径变化)#
- x,y:自变量,遍历整个平面或区域。
- P(x,y),Q(x,y):因变量,随位置变化。
- T,n:沿边界 C 的方向向量,随弧长参数变化。
- dx,dy,ds,dA:微元,随积分位置变化。
3. 已知量 vs 待求量(典型问题设定)#
| 问题类型 | 已知量 | 待求量 |
|---|
| 直接计算(最常见) | P,Q 的表达式;C 的参数方程 | 线积分 ∮C 的值,或验证与二重积分相等 |
| 保守场判定 | P,Q 的表达式 | 判断 ∂x∂Q=?∂y∂P,从而确定积分是否与路径无关 |
| 求原函数 | P,Q 满足 ∂xQ=∂yP | 势函数 u(x,y),使得 du=Pdx+Qdy |
| 面积计算 | 边界 C 的参数方程 | 区域 D 的面积 |
| 反问题(场重构) | 边界上的测量值(环量/通量) | 推断区域内部的旋度或散度分布 |
四、一个关键提醒:量纲一致性#
格林公式作为恒等式,左右两边量纲必须相同:
[P]⋅[L]∮CPdx=[Q]/[L]⋅[L]2=[Q]⋅[L]∬D(∂x∂Q)dA
因此,公式成立的前提是:
[P]=[Q](即 P 和 Q 具有相同的量纲)
例如:
- 力场中:P,Q 都是力(N),左边 N⋅m,右边 (N/m)⋅m2=N⋅m ✓
- 速度场中:P,Q 都是速度(m/s),左边 m2/s,右边 (s−1)⋅m2=m2/s ✓
如果 P 和 Q 量纲不同(例如 P 是力、Q 是速度),格林公式在物理上无意义,但在纯数学中仍可形式化计算。
格林公式的适用条件与范围#
一、严格数学条件(缺一不可)#
设标准格林公式为:
∮CPdx+Qdy=∬D(∂x∂Q−∂y∂P)dA
1. 对区域 D 的要求#
| 条件 | 具体说明 |
|---|
| 有界 | D 必须能被某个有限半径的圆包含 |
| 闭区域 | 包含边界 C 本身,即 D=D∘∪∂D |
| 可求面积 | 边界 C 必须是零面积集(不能是充满面积的”胖边界”) |
| 单连通或复连通 | 允许有有限个”洞”,但洞必须是闭区域 |
2. 对边界 C=∂D 的要求#
| 条件 | 具体说明 |
|---|
| 分段光滑 | 由有限条光滑曲线弧连接而成,允许有有限个角点(如矩形) |
| 简单闭曲线 | 不自交,且首尾相连形成封闭回路 |
| 正向(诱导定向) | 沿 C 行走时,D 始终在左侧;外边界逆时针,内边界(洞)顺时针 |
反例:若 C 自交(如”8”字形),则 C 围成多个区域,公式需对各区域分别应用。
3. 对函数 P,Q 的要求(最关键)#
| 条件 | 具体说明 |
|---|
| 在 D 上定义 | P,Q 在闭区域 D(含边界)上处处有定义 |
| 一阶偏导数存在 | ∂x∂P,∂y∂P,∂x∂Q,∂y∂Q 在 D 上存在 |
| 偏导数连续 | 上述四个偏导数在 D 上连续 |
⚠️ 重要:偏导数”存在”不足以保证公式成立,必须连续。这是格林公式与牛顿-莱布尼兹公式的重要区别之一。
二、适用范围与推广#
1. 单连通区域(标准情形)#
D 内无洞,边界只有一条外曲线 C。
典型例子:圆盘、椭圆内部、矩形、任意简单多边形内部。
2. 复连通区域(带洞的情形)#
D 内有有限个洞,边界由外边界 C0 和内边界 C1,C2,…,Cn 组成。
公式修正为:
∮C0Pdx+Qdy+∑k=1n∮CkPdx+Qdy=∬D(∂x∂Q−∂y∂P)dA
定向规则:
- C0:逆时针(正向)
- C1,…,Cn:顺时针(对洞而言是正向,即沿边界走时区域在左侧)
物理意义:外边界环量 + 各洞内边界环量 = 区域总涡量。
3. 非闭曲线的”补形”技巧#
若 C 不是闭曲线,但已知 ∂xQ=∂yP(保守场),则:
- 可任选一条从起点到终点的路径计算积分
- 或补上一段辅助线使其闭合,再用格林公式
三、失效情形与反例#
| 失效原因 | 反例 | 说明 |
|---|
| P,Q 在 D 内不连续 | P=x2+y2−y,Q=x2+y2x 在单位圆上 | 在原点 (0,0) 处无定义,偏导数不连续 |
| 偏导数在 D 内不连续 | 同上 | 虽然 P,Q 在除原点外处处光滑,但原点处偏导数无界 |
| C 不是简单曲线 | 自交的”8”字形 | 需分解为多个简单闭曲线分别处理 |
| C 不是分段光滑 | 科赫雪花曲线(处处不可微) | 属于分形边界,需用更一般的测度论工具 |
| D 无界 | 上半平面、全平面 | 需附加积分收敛条件,公式不自动成立 |
| 定向错误 | 外边界顺时针 | 结果差一个负号 |
经典反例详解:
∮x2+y2=1x2+y2−ydx+x2+y2xdy=2π=0
虽然计算得 ∂xQ−∂yP=0(在除原点外处处),但原点处 P,Q 无定义,格林公式不能直接用。正确做法是挖去原点处小圆盘,在复连通区域上应用,结果恰好是 2π。
四、与其他积分定理的关系#
| 定理 | 维度 | 关系 |
|---|
| 牛顿-莱布尼兹公式 | 1D | ∫abf′(x)dx=f(b)−f(a),格林公式是其 2D 推广 |
| 斯托克斯公式 | 3D | 格林公式是斯托克斯公式在平面上的特例 |
| 高斯公式(散度定理) | 3D | 格林公式的”法向形式” ∮CF⋅nds=∬D∇⋅FdA 是高斯公式在 2D 的对应 |
| 柯西积分公式 | 复分析 | 复变函数中的围道积分本质上依赖格林公式(或斯托克斯公式) |
五、判断”能否使用”的决策流程#
│否 → 格林公式失效,需用其他方法(如直接参数化、留数定理等)
六、一句话总结#
格林公式成立的核心:在有界闭区域上,分段光滑的简单闭曲线包围的区域中,若 P,Q 及其一阶偏导数处处连续,则线积分与二重积分相等。奇点、不光滑边界、错误定向是三大常见陷阱。