知识点

格林公式

定义

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

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

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

对复连通区域,格林公式仍然成立,但此时 LL 应理解为区域 DD 的全部边界,且外边界取正向(逆时针),内边界取负向(顺时针)。

使用格林公式前要检查两件事:PP、QQ 在 DD 上(包括 DD 的内部)具有一阶连续偏导数;LL 取正向。若 DD 内有使 PP 或 QQ 无定义的奇点,不能直接套公式,要在挖去奇点后的区域上使用,并把新增的边界计入 LL。

证明

下面给出 DD 同时为 xx 型与 yy 型区域时的证明;一般的分段光滑边界区域可以分割成有限个这样的区域,公共边界上的积分因方向相反而抵消。

设

D={(x,y):a≤x≤b, y1(x)≤y≤y2(x)}D=\{(x,y):a\le x\le b,\ y_1(x)\le y\le y_2(x)\}

也是 yy 型区域

D={(x,y):c≤y≤d, x1(y)≤x≤x2(y)}.D=\{(x,y):c\le y\le d,\ x_1(y)\le x\le x_2(y)\}.

先计算含 QQ 的项。固定 yy,对 xx 积分:

∬D∂Q∂x dx dy=∫cd[Q(x2(y),y)−Q(x1(y),y)] dy.\iint\limits_D\frac{\partial Q}{\partial x}\,\mathrm{d}x\,\mathrm{d}y =\int_c^d\Bigl[Q\bigl(x_2(y),y\bigr)-Q\bigl(x_1(y),y\bigr)\Bigr]\,\mathrm{d}y.

沿正向边界,右边界 x=x2(y)x=x_2(y) 从下往上走,左边界 x=x1(y)x=x_1(y) 从上往下走,所以

∮LQ dy=∫cdQ(x2(y),y) dy+∫dcQ(x1(y),y) dy=∫cd[Q(x2(y),y)−Q(x1(y),y)] dy.\oint_L Q\,\mathrm{d}y =\int_c^d Q\bigl(x_2(y),y\bigr)\,\mathrm{d}y +\int_d^c Q\bigl(x_1(y),y\bigr)\,\mathrm{d}y =\int_c^d\Bigl[Q\bigl(x_2(y),y\bigr)-Q\bigl(x_1(y),y\bigr)\Bigr]\,\mathrm{d}y.

两者相等。再计算含 PP 的项。固定 xx,对 yy 积分:

∬D−∂P∂y dx dy=∫ab[P(x,y1(x))−P(x,y2(x))] dx.\iint\limits_D-\frac{\partial P}{\partial y}\,\mathrm{d}x\,\mathrm{d}y =\int_a^b\Bigl[P\bigl(x,y_1(x)\bigr)-P\bigl(x,y_2(x)\bigr)\Bigr]\,\mathrm{d}x.

沿正向边界,下边界从 aa 到 bb,上边界从 bb 到 aa,所以

∮LP dx=∫abP(x,y1(x)) dx+∫baP(x,y2(x)) dx=∫ab[P(x,y1(x))−P(x,y2(x))] dx.\oint_L P\,\mathrm{d}x =\int_a^b P\bigl(x,y_1(x)\bigr)\,\mathrm{d}x +\int_b^a P\bigl(x,y_2(x)\bigr)\,\mathrm{d}x =\int_a^b\Bigl[P\bigl(x,y_1(x)\bigr)-P\bigl(x,y_2(x)\bigr)\Bigr]\,\mathrm{d}x.

两式相加即得

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

性质

计算平面图形面积

在格林公式中取 P=−yP=-y,Q=xQ=x,则

∂Q∂x−∂P∂y=1−(−1)=2.\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}=1-(-1)=2.

于是

2∬Ddx dy=∮Lx dy−y dx,2\iint\limits_D\mathrm{d}x\,\mathrm{d}y =\oint_L x\,\mathrm{d}y-y\,\mathrm{d}x,

即区域 DD 的面积为

A=∬Ddx dy=12∮Lx dy−y dx,A=\iint\limits_D\mathrm{d}x\,\mathrm{d}y =\frac12\oint_L x\,\mathrm{d}y-y\,\mathrm{d}x,

其中 LL 为 DD 的正向边界。

取 P=0P=0、Q=xQ=x 得 A=∮Lx dyA=\oint_L x\,\mathrm{d}y;取 P=−yP=-y、Q=0Q=0 得 A=−∮Ly dxA=-\oint_L y\,\mathrm{d}x。这两个公式与上面的对称形式等价,选择哪一个取决于边界方程用哪个变量作参数更方便。

边界写成 x=x(y)x=x(y) 时,∮Lx dy\oint_L x\,\mathrm{d}y 往往只需计算一个积分;边界写成 y=y(x)y=y(x) 时,−∮Ly dx-\oint_L y\,\mathrm{d}x 可能更省事。方向必须保持正向,否则结果差一个负号。

例子

用格林公式计算单位正方形的边界积分

设 D=[0,1]×[0,1]D=[0,1]\times[0,1],LL 为其正向边界。取 P=−yP=-y、Q=xQ=x,则

∮Lx dy−y dx=∬D2 dx dy=2⋅1=2.\oint_L x\,\mathrm{d}y-y\,\mathrm{d}x =\iint\limits_D 2\,\mathrm{d}x\,\mathrm{d}y =2\cdot 1=2.

这也正好是 DD 的面积的两倍。

椭圆的面积

椭圆 x2a2+y2b2=1\dfrac{x^2}{a^2}+\dfrac{y^2}{b^2}=1(a,b>0a,b>0)的正向边界可参数化为

x=acos⁡t,y=bsin⁡t,0≤t≤2π.x=a\cos t,\qquad y=b\sin t,\qquad 0\le t\le 2\pi.

于是

x dy−y dx=acos⁡t⋅bcos⁡t dt−bsin⁡t⋅(−asin⁡t) dt=ab dt.x\,\mathrm{d}y-y\,\mathrm{d}x =a\cos t\cdot b\cos t\,\mathrm{d}t -b\sin t\cdot(-a\sin t)\,\mathrm{d}t =ab\,\mathrm{d}t.

代入面积公式:

A=12∫02πab dt=πab.A=\frac12\int_0^{2\pi}ab\,\mathrm{d}t=\pi ab.

探索

边界积分与区域积分怎样对应?

用圆盘例子比较格林公式两侧。

正在准备演示…

题集

练习

  1. 使用格林公式时,下列哪组条件足够?

    选择答案
    参考解答

    格林公式要求 PP、QQ 在 DD 上有一阶连续偏导数,边界取正向;复连通区域也可以使用,只是要取全部边界。

  2. 设 D=[0,1]×[0,1]D=[0,1]\times[0,1],LL 为其正向边界,则 ∮Lx dy−y dx=\displaystyle\oint_L x\,\mathrm{d}y-y\,\mathrm{d}x= 。

    参考解答

    由格林公式,被积表达式对应 ∂Q∂x−∂P∂y=1−(−1)=2\dfrac{\partial Q}{\partial x}-\dfrac{\partial P}{\partial y}=1-(-1)=2,而 DD 的面积为 11,所以积分为 2⋅1=22\cdot 1=2。