4. 二重积分的计算

4.1 二重积分算的是什么

第 1.9 节把概率写成了曲顶柱体的体积。区域 $D$ 在 $xy$ 平面上,曲面 $z = f(x,y)$ 在它上方, 这个立体的体积记作

$$ \iint_D f(x,y)\,dA $$

其中 $dA$ 表示 $D$ 内的一小块面积。把 $D$ 分成许多小块,第 $i$ 块的面积记作 $\Delta A_i$, 在块内取一点 $(x_i, y_i)$,算出高度 $f(x_i, y_i)$,把这些高度与面积相乘再求和。 小块分得越细,和越接近一个确定的数,这个数就是二重积分。 所以二重积分就是曲面下方的体积。

极简例子:取 $f(x,y) = 1$。曲面是高度为 1 的水平面,柱体退化成一块薄板,体积等于底面积乘 1,所以 $\iint_D 1\,dA = D$ 的面积。矩形 $R = [0,2] \times [0,1]$ 上, $\iint_R 1\,dA = 2 \times 1 = 2$。
极简例子:把 $f$ 换成常数 3,高度是原来的 3 倍,体积也是 3 倍: $\iint_R 3\,dA = 3 \times 2 = 6$。

概率论里,如果 $f$ 是 $(X,Y)$ 的联合密度函数,那么 $\iint_D f\,dA = P\bigl((X,Y) \in D\bigr)$。 密度在全平面上的积分等于 1,就是整个曲顶柱体的体积等于 1。

4.2 化成一先一后

直接按定义求极限不现实,通常把二重积分化成两次一元积分,这一步依据的是 Fubini 定理。 做法是切片:固定一个变量,沿一个方向把区域切成细条,先算一条细条上的积分,再把所有细条累加。

图 1:三角形是底面 $G$,顶面是曲面 $z = xy$,里面的竖线标出每一点上的函数值。 这个立体的体积是 $\iint_G xy\,dA = 2$,计算过程见下面的极简例子。 4.4 节的例题 1 用同一块底面算了另一个被积函数 $x+y$,体积就变成了 4。

如果区域 $D$ 可以写成

$$ a \le x \le b, \qquad y_1(x) \le y \le y_2(x) $$

也就是竖直切下去,每条竖条的下端在 $y_1(x)$、上端在 $y_2(x)$,那么

$$ \iint_D f(x,y)\,dA = \int_a^b \left(\int_{y_1(x)}^{y_2(x)} f(x,y)\,dy\right) dx $$

内层对 $y$ 积分,积分限可以含 $x$;外层对 $x$ 积分,限是常数。这样的区域叫 X 型区域。

如果 $D$ 可以写成

$$ c \le y \le d, \qquad x_1(y) \le x \le x_2(y) $$

也就是水平切下去,每条横条的左端在 $x_1(y)$、右端在 $x_2(y)$,那么

$$ \iint_D f(x,y)\,dA = \int_c^d \left(\int_{x_1(y)}^{x_2(y)} f(x,y)\,dx\right) dy $$

这样的区域叫 Y 型区域。两式算的是同一个体积,区别只在切的方向。

写的时候注意两点。第一,内层积分的限可以含外层变量,外层积分的限必须是常数。 第二,内层的微分符号要和它的积分限对应,$dy$ 配 $y$ 的限,$dx$ 配 $x$ 的限。

同一个三角形的两种切片方式
图 2:区域 $D$ 由 $y = 0$、$y = x$、$x = 2$ 围成。 左图竖直切,$x$ 固定时 $y$ 从 0 到 $x$。右图水平切,$y$ 固定时 $x$ 从 $y$ 到 2。
极简例子:矩形 $R = [0,1] \times [0,2]$,被积函数 $f = xy$。两条限都是常数,可以直接拆成两个一元积分相乘: $$ \iint_R xy\,dA = \left(\int_0^1 x\,dx\right)\left(\int_0^2 y\,dy\right) = \frac{1}{2} \cdot 2 = 1 $$
极简例子:图 1 里的三角形 $T$:$0 \le y \le x \le 2$,被积函数 $f = xy$。 竖直切时内层的上限是变量: $$ \iint_T xy\,dA = \int_0^2 \left(\int_0^x xy\,dy\right) dx = \int_0^2 \frac{1}{2}x^3\,dx = 2 $$ 和图 1 显示的体积一致。

4.3 怎么确定积分限

定积分限分三步。

  1. 画图。把边界曲线画出来,标出交点。
  2. 写不等式。选一个变量先积,把区域写成「外层变量的范围」加上「内层的上下界」。
  3. 检查。外层限是常数,内层限只含外层变量。

以图 1 的三角形为例。竖直切时,$x$ 从 0 到 2,对每个 $x$,下边界是 $y = 0$,上边界是 $y = x$, 所以 $0 \le y \le x$。水平切时,$y$ 从 0 到 2,对每个 $y$,左边界是 $x = y$,右边界是 $x = 2$, 所以 $y \le x \le 2$。两种写法描述的是同一个区域。

4.4 直角坐标下的例题

例题 1

计算 $\iint_D (x+y)\,dA$,其中 $D$ 是由直线 $y = x$、$y = 0$ 和 $x = 2$ 围成的三角形区域。

解:

  1. 画图。三条边界围成三角形,顶点是 $(0,0)$、$(2,0)$、$(2,2)$。
  2. 定限。竖直切:$x$ 从 0 到 2。对每个 $x$,下边界是 $y = 0$,上边界是 $y = x$, 所以 $0 \le y \le x$。
  3. 计算: $$ \begin{align} \iint_D (x+y)\,dA &= \int_0^2 \int_0^x (x+y)\,dy\,dx \\ &= \int_0^2 \left[xy + \frac{1}{2}y^2\right]_0^x dx \\ &= \int_0^2 \left(x^2 + \frac{1}{2}x^2\right) dx \\ &= \int_0^2 \frac{3}{2}x^2\,dx \\ &= \left[\frac{1}{2}x^3\right]_0^2 \\ &= 4 \end{align} $$

例题 2

计算 $\iint_D (x^2+y^2)\,dA$,其中 $D$ 是以原点为中心、半径为 2 的圆形区域。

解:

被积函数在极坐标下等于 $r^2$,面积元素是 $r\,dr\,d\theta$。 圆形区域用极坐标表示最简单:$\theta$ 从 0 到 $2\pi$,$r$ 从 0 到 2。

$$ \begin{align} \iint_D (x^2+y^2)\,dA &= \int_0^{2\pi} \int_0^2 r^2 \cdot r\,dr\,d\theta \\ &= \int_0^{2\pi} \left[\frac{1}{4}r^4\right]_0^2 d\theta \\ &= \int_0^{2\pi} 4\,d\theta \\ &= 8\pi \end{align} $$

例题 3

计算 $\iint_D xy\,dA$,其中 $D$ 是由抛物线 $y = x^2$ 和直线 $y = 2x$ 围成的区域。

解:

先求交点:$x^2 = 2x$,解得 $x = 0$ 或 $x = 2$,交点是 $(0,0)$ 和 $(2,4)$。

两种切法都可以。

下面用水平切:

$$ \begin{align} \iint_D xy\,dA &= \int_0^4 \int_{\frac{y}{2}}^{\sqrt{y}} xy\,dx\,dy \\ &= \int_0^4 \left[\frac{1}{2}x^2y\right]_{\frac{y}{2}}^{\sqrt{y}} dy \\ &= \int_0^4 \left(\frac{1}{2}y^2 - \frac{1}{8}y^3\right) dy \\ &= \left[\frac{1}{6}y^3 - \frac{1}{32}y^4\right]_0^4 \\ &= \frac{32}{3} - 8 \\ &= \frac{8}{3} \end{align} $$

用竖直切验算,结果同样是 $\frac{8}{3}$。

4.5 极坐标下的计算

区域是圆、扇形或圆环时,用极坐标。令 $x = r\cos\theta$,$y = r\sin\theta$, 被积函数里的 $x$、$y$ 都要换掉:

$$ f(x,y) = f(r\cos\theta,\, r\sin\theta) $$

面积元素也要换。极坐标下的一小块面积约等于 $(r\,d\theta)(dr)$,所以

$$ dA = r\,dr\,d\theta $$

其中的 $r$ 不能漏。原因很直接:同样转过 $d\theta$,半径越大走过的路程越长,弧长是 $r\,d\theta$, 不是 $d\theta$。

极坐标下的面积元素
图 3:极坐标下的面积元素。径向宽度是 $dr$,弧长是 $r\,d\theta$, 两者的乘积就是面积 $r\,dr\,d\theta$。

于是

$$ \iint_D f(x,y)\,dA = \int_\alpha^\beta \int_{r_1(\theta)}^{r_2(\theta)} f(r\cos\theta,\, r\sin\theta)\, r\,dr\,d\theta $$

极简例子:单位圆盘 $D$:$0 \le r \le 1$,$0 \le \theta \le 2\pi$。 $$ \iint_D 1\,dA = \int_0^{2\pi}\int_0^1 r\,dr\,d\theta = 2\pi \cdot \frac{1}{2} = \pi $$ 和圆面积一致。漏掉 $r$ 会算出 $2\pi$,那个数是周长,不是面积。
极简例子:单位圆盘上算 $\iint_D (x^2+y^2)\,dA$。$\theta$ 从 0 到 $2\pi$,$r$ 从 0 到 1, 被积函数是 $r^2$,面积元素是 $r\,dr\,d\theta$,所以 $$ \iint_D (x^2+y^2)\,dA = \int_0^{2\pi}\int_0^1 r^3\,dr\,d\theta = 2\pi \cdot \frac{1}{4} = \frac{\pi}{2} $$

4.6 能简化的情形

有四种情形可以省事。每种都先写清条件。

4.6.1 被积函数是常数

如果 $f$ 与 $x$、$y$ 都无关,二重积分等于常数乘以区域的面积: $\iint_D c\,dA = c \cdot S_D$,其中 $S_D$ 是 $D$ 的面积。

例题 4

计算 $\iint_D 5\,dA$,其中 $D$ 是由 $y = x^2$ 和 $y = 2x$ 围成的区域。

解:

先算面积:

$$ \begin{align} S_D &= \int_0^2 (2x - x^2)\,dx \\ &= \left[x^2 - \frac{1}{3}x^3\right]_0^2 \\ &= 4 - \frac{8}{3} \\ &= \frac{4}{3} \end{align} $$

所以 $\iint_D 5\,dA = 5 \cdot \frac{4}{3} = \frac{20}{3}$。

4.6.2 被积函数可分离,且区域是矩形

条件有两条,缺一不可。被积函数要能写成 $f(x)g(y)$,区域要是矩形 $[a,b] \times [c,d]$。这时

$$ \iint_D f(x)g(y)\,dA = \left(\int_a^b f(x)\,dx\right)\left(\int_c^d g(y)\,dy\right) $$

区域不是矩形时不能这样拆。例如三角形的两条限会互相牵制,内层的上下界含外层变量, 两个积分拆不开。

例题 5

计算 $\iint_D e^x \sin y\,dA$,其中 $D = [0,1] \times [0,\pi]$。

解:

区域是矩形,被积函数是 $e^x$ 与 $\sin y$ 的乘积,两条条件都满足。

$$ \begin{align} \iint_D e^x \sin y\,dA &= \left(\int_0^1 e^x\,dx\right)\left(\int_0^\pi \sin y\,dy\right) \\ &= (e-1) \cdot 2 \\ &= 2(e-1) \end{align} $$

4.6.3 对称性

如果区域关于 $x$ 轴和 $y$ 轴都对称,被积函数关于 $x$、关于 $y$ 也都是偶函数, 可以先算第一象限里的四分之一,再乘 4。

例题 6

计算 $\iint_D x^2\,dA$,其中 $D$ 是圆 $x^2+y^2 \le 4$。

解:

区域关于两条轴都对称,$x^2$ 关于 $x$、关于 $y$ 都是偶函数,所以可以只算第一象限的四分之一圆, 再乘 4。

$$ 4 \int_0^2 \int_0^{\sqrt{4-x^2}} x^2\,dy\,dx = 4 \int_0^2 x^2\sqrt{4-x^2}\,dx $$

最后这个一元积分用 $x = 2\sin\theta$ 换元。$x$ 从 0 到 2 对应 $\theta$ 从 0 到 $\frac{\pi}{2}$, $dx = 2\cos\theta\,d\theta$:

$$ \begin{align} \int_0^2 x^2\sqrt{4-x^2}\,dx &= \int_0^{\frac{\pi}{2}} 4\sin^2\theta \cdot 2\cos\theta \cdot 2\cos\theta\,d\theta \\ &= 16 \int_0^{\frac{\pi}{2}} \sin^2\theta\cos^2\theta\,d\theta \end{align} $$

用 $\sin^2\theta\cos^2\theta = \frac{1}{8}(1-\cos 4\theta)$, $\int_0^{\frac{\pi}{2}} \frac{1}{8}(1-\cos 4\theta)\,d\theta = \frac{1}{8}\cdot\frac{\pi}{2} = \frac{\pi}{16}$。 所以一元积分等于 $16 \cdot \frac{\pi}{16} = \pi$,再乘 4 得 $4\pi$。

极简例子:单位圆盘上算 $\iint_D x^2\,dA$。由对称性,$\iint_D x^2\,dA = \iint_D y^2\,dA$, 两者相加等于 $\iint_D (x^2+y^2)\,dA = \frac{\pi}{2}$,所以各是 $\frac{\pi}{4}$。 不需要换元,也不需要例题 6 那一长串计算。

4.6.4 被积函数只含一个变量

如果被积函数只含 $y$,先对 $x$ 积分,内层积分就是横条的长度,外层只剩一次一元积分。 只含 $x$ 时同理,先对 $y$ 积分,内层是竖条的长度。

例题 7

计算 $\iint_D y^2\,dA$,其中 $D$ 是由 $y = x$、$y = 0$ 和 $x = 2$ 围成的三角形区域。

解:

被积函数只含 $y$,先对 $x$ 积分。$y$ 从 0 到 2,对每个 $y$,$x$ 从 $y$ 到 2。

$$ \begin{align} \iint_D y^2\,dA &= \int_0^2 y^2 \int_y^2 dx\,dy \\ &= \int_0^2 y^2(2-y)\,dy \\ &= \int_0^2 (2y^2 - y^3)\,dy \\ &= \left[\frac{2}{3}y^3 - \frac{1}{4}y^4\right]_0^2 \\ &= \frac{16}{3} - 4 \\ &= \frac{4}{3} \end{align} $$

内层 $\int_y^2 dx = 2-y$ 是这条横条的长度,外层把它乘以 $y^2$ 再按 $y$ 累加。 计算质心、转动惯量时经常这样处理。

4.7 回到概率论

联合密度为 $f(x,y)$ 时,随机向量落在区域 $G$ 内的概率是

$$ P\bigl((X,Y) \in G\bigr) = \iint_G f(x,y)\,dA $$

这就是 1.9 节那个曲顶柱体的体积。计算它,就是把 4.2 到 4.5 的方法用一遍。

极简例子:设 $(X,Y)$ 在矩形 $[0,2] \times [0,1]$ 上服从均匀分布,密度 $f = \frac{1}{2}$。 取 $G = [0,1] \times [0,1]$,则 $P = \iint_G \frac{1}{2}\,dA = \frac{1}{2} \times 1 = \frac{1}{2}$。 这个结果也等于面积之比,$G$ 的面积占整个矩形面积的一半。

例题 8

设 $(X,Y)$ 服从标准二元正态分布,密度 $f(x,y) = \frac{1}{2\pi}e^{-\frac{x^2+y^2}{2}}$。求 $(X,Y)$ 落在单位圆盘内的概率。

解:

圆形区域用极坐标。$x^2+y^2 = r^2$,面积元素是 $r\,dr\,d\theta$, $\theta$ 从 0 到 $2\pi$,$r$ 从 0 到 1。

$$ \begin{align} P &= \int_0^{2\pi} \int_0^1 \frac{1}{2\pi}e^{-\frac{r^2}{2}}\, r\,dr\,d\theta \\ &= \int_0^1 e^{-\frac{r^2}{2}}\, r\,dr \\ &= \left[-e^{-\frac{r^2}{2}}\right]_0^1 \\ &= 1 - e^{-\frac{1}{2}} \\ &\approx 0.393 \end{align} $$

这个数就是 1.9 节图 15 里把 $G$ 取成单位圆时显示的体积。换了坐标系,算的是同一件事。