卷二 · 数学之美13 分钟阅读

函数积分十:二重积分的计算

定理1

如果ff在I=[a,b]×[c,d]I = [a,b] \times [c,d]上可积,那么单变量函数φ(x)=∫‾cdf(x,y)dy\displaystyle \varphi(x) = \underline \int_c^d f(x, y) \mathrm{d} y和ψ(x)=∫‾cdf(x,y)dy\displaystyle \psi(x) = \overline \int_c^d f(x, y) \mathrm{d} y在区间[a,b][a,b]上可积,且

∫Ifdσ=∫abφ(x)dx=∫abψ(x)dx \int_I f \mathrm{d} \sigma = \int_a^b \varphi(x) \mathrm{d} x = \int_a^b \psi(x) \mathrm{d} x

证:分别对[a,b][a,b]和[c,d][c,d]的分割

πx:a=x0<x1<⋯<xn=b,πy:c=y0<y1<⋯<ym=d, \begin{aligned} \pi_x: a = x_0 < x_1 < \cdots < x_n = b, \\ \pi_y: c = y_0 < y_1 < \cdots < y_m = d, \end{aligned}

令
Ii=[xi−1,xi](i=1,2,⋯ ,n)Jj=[yj−1,yj](j=1,2,⋯ ,m) \begin{aligned} & I_i = [x_{i-1}, x_i] \quad (i=1,2,\cdots,n) \\ & J_j = [y_{j-1}, y_j] \quad (j=1,2,\cdots,m) \end{aligned}

则子矩形
Ii×Ij(i=1,2,⋯ ,n;j=1,2,⋯ ,m) I_i \times I_j \quad (i=1,2,\cdots,n; j=1,2,\cdots,m)

形成了矩形II的分割π=πx×πy\pi = \pi_x \times \pi_y。令A=∫Ifdσ\displaystyle A = \int_I f \mathrm{d} \sigma,由积分存在的定义可知,对任意的ε>0\varepsilon > 0,存在δ>0\delta > 0,当II的分割π\pi满足∥π∥<δ\Vert \pi \Vert < \delta时,必有
A−ε<∑i=1n∑j=1mf(ξi,ηj)ΔxiΔyj<A+ε A - \varepsilon < \sum_{i=1}^n \sum_{j=1}^m f(\xi_i, \eta_j) \Delta x_i \Delta y_j < A + \varepsilon

其中ξi∈Ii,ηj∈Jj(i=1,2,⋯ ,n;j=1,2,⋯ ,m)\xi_i \in I_i, \eta_j \in J_j (i=1,2,\cdots,n;j=1,2,\cdots,m),现取分割πx,πy\pi_x, \pi_y满足∥πx∥<δ/2,∥πy∥<δ/2\Vert \pi_x \Vert < \delta / \sqrt 2, \Vert \pi_y \Vert < \delta / \sqrt 2,那么∥π∥<δ\Vert \pi \Vert < \delta,从而上式成立,所以有
A−ε≤∑i=1n∑j=1minf⁡f(ξi,ηj)ΔxiΔyj≤∑i=1n∑j=1msup⁡f(ξi,ηj)ΔxiΔyj≤A+ε \begin{aligned} A - \varepsilon &\le \sum_{i=1}^n \sum_{j=1}^m \inf f(\xi_i, \eta_j) \Delta x_i \Delta y_j \\ & \le \sum_{i=1}^n \sum_{j=1}^m \sup f(\xi_i, \eta_j) \Delta x_i \Delta y_j \\ & \le A + \varepsilon \end{aligned}

而∑j=1minf⁡f(ξi,ηj)Δyj\displaystyle \sum_{j=1}^m \inf f(\xi_i, \eta_j) \Delta y_j表示函数f(ξi)f(\xi_i)在[c,d][c,d]上的下和,所以
∑j=1minf⁡f(ξi,ηj)Δyj≤∫‾cdf(ξi,y)dy=φ(ξi) \sum_{j=1}^m \inf f(\xi_i, \eta_j) \Delta y_j \le \underline \int_c^d f(\xi_i, y) \mathrm{d} y = \varphi(\xi_i)

同理
∑j=1msup⁡f(ξi,ηj)Δyj≤∫‾cdf(ξi,y)dy=ψ(ξi) \sum_{j=1}^m \sup f(\xi_i, \eta_j) \Delta y_j \le \overline \int_c^d f(\xi_i, y) \mathrm{d} y = \psi(\xi_i)

从而
A−ε≤∑i=1nφ(ξi)Δxi≤∑i=1nψ(ξi)Δxi≤A+ε A - \varepsilon \le \sum_{i=1}^n \varphi(\xi_i) \Delta x_i \le \sum_{i=1}^n \psi(\xi_i) \Delta x_i \le A + \varepsilon

即
lim⁡∥πx∥→0∑i=1nφ(ξi)Δxi=lim⁡∥π∥→0∑i=1nψ(ξi)Δxi=A \lim_{\Vert \pi_x \Vert \to 0} \sum_{i=1}^n \varphi(\xi_i) \Delta x_i = \lim_{\Vert \pi \Vert \to 0} \sum_{i=1}^n \psi(\xi_i) \Delta x_i = A

成立

Q.E.D.

定理2

设ff在I=[a,b]×[c,d]I = [a,b] \times [c, d]上可积,如果对每一个x∈[a,b]x \in [a,b],函数f(x,y)f(x, y)在[c,d][c, d]上可积,则

∫Ifdσ=∫ab(∫cdf(x,y)dy)dx \int_I f \mathrm{d} \sigma = \int_a^b \left(\int_c^d f(x, y) \mathrm{d} y \right) \mathrm{d} x

上面等式右边称为累次积分,也可以记为
∫abdx∫cdf(x,y)dy \int_a^b \mathrm{d} x \int_c^d f(x, y) \mathrm{d} y

同样,如果对于每一个y∈[c,d]y \in [c, d],函数f(x,y)f(x, y)在[a,b][a,b]上可积,那么有
∫Ifdσ=∫cd(∫abf(x,y)dx)dy \int_I f \mathrm{d} \sigma = \int_c^d \left(\int_a^b f(x, y) \mathrm{d} x \right) \mathrm{d} y

上面等式右边也可以记为
∫cddy∫abf(x,y)dx \int_c^d \mathrm{d} y \int_a^b f(x, y) \mathrm{d} x

证:由定理1可知,

φ(x)=ψ(x)=∫cdf(x,y)dy \varphi(x) = \psi(x) = \int_c^d f(x, y) \mathrm{d} y

所以
∫Ifdσ=∫abφ(x)dx=∫abdx∫cdf(x,y)dy \int_I f \mathrm{d} \sigma = \int_a^b \varphi(x) \mathrm{d} x = \int_a^b \mathrm{d} x \int_c^d f(x, y) \mathrm{d} y

后半部分同样的证明方法。

Q.E.D.

定理3

设ff是[a,b]×[c,d][a,b] \times [c, d]上的连续函数,则有

∫cddy∫abf(x,y)dx=∫abdx∫cdf(x,y)dy \int_c^d \mathrm{d} y \int_a^b f(x, y) \mathrm{d} x = \int_a^b \mathrm{d} x \int_c^d f(x, y) \mathrm{d} y

证:由于ff连续,从而φ(x)\varphi(x)与ψ(x)\psi(x)都可积,再由定理2易证。

Q.E.D.

定理4

设点集

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

其中函数y1,y2y_1,y_2在[a,b][a,b]上连续,函数ff在BB上可积。如果对任意的x∈[a,b]x \in [a,b],单变量积分
∫y1(x)y2(x)f(x,y)dy \int_{y_1(x)}^{y_2(x)} f(x, y) \mathrm{d} y

存在,那么
∫Bfdσ=∫abdx∫y1(x)y2(x)f(x,y)dy \int_B f \mathrm{d} \sigma = \int_a^b \mathrm{d} x \int_{y_1(x)}^{y_2(x)} f(x, y) \mathrm{d} y

证:令c=inf⁡y1([a,b]),d=sup⁡y2([a,b])c = \inf y_1([a, b]), d = \sup y_2([a, b]),从而I=[a,b]×[c,d]⊃BI = [a, b] \times [c, d] \supset B,由于ff在BB上可积,从而fBf_B在II上可积,而且

∫Bfdσ=∫IfBdσ \int_B f \mathrm{d} \sigma = \int_I f_B \mathrm{d} \sigma

显然易知,对每一个x∈[a,b]x \in [a,b],fB(x,y)f_B(x, y)在[c,d][c, d]上可积,所以由定理2可知,
∫IfBdσ=∫abdx∫cdfB(x,y)dy=∫abdx∫y1(x)y2(x)fBdy=∫abdx∫y1(x)y2(x)f(x,y)dy \int_I f_B \mathrm{d} \sigma = \int_a^b \mathrm{d} x \int_c^d f_B(x, y) \mathrm{d} y = \int_a^b \mathrm{d} x \int_{y_1(x)}^{y_2(x)} f_B \mathrm{d} y = \int_a^b \mathrm{d} x\int_{y_1(x)}^{y_2(x)} f(x, y) \mathrm{d} y

Q.E.D.

定义1:正则映射

设有界闭域D⊂R2D \subset \mathbb{R}^2,连续函数F:D→RF: D \to \mathbb{R},映射φ\boldsymbol{\varphi}由公式

x=x(u,v),y=y(u,v)((u,v)∈Δ) x = x(u, v), \quad y = y(u, v) \quad ((u, v) \in \Delta)

定义,其中Δ\Delta是uvuv平面上的有界闭区域,设映射φ\boldsymbol{\varphi}是正则的,即φ\boldsymbol{\varphi}是从Δ\Delta到DD上的一对一的映射,φ\boldsymbol{\varphi}在Δ\Delta上连续可导,并且
∂(x,y)∂(u,v)=∣∂x∂u∂x∂v∂y∂u∂y∂v∣≠0 \frac{\partial (x, y) }{\partial (u, v)} = \begin{vmatrix} \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} \end{vmatrix} \neq 0

在Δ\Delta上成立。

定理5

设正则映射φ\boldsymbol{\varphi}把R2\mathbb{R}^2中以

(u0,v0),(u0+h,v0),(u0+h,v0+h),(u0,v0+h) (u_0, v_0), \quad (u_0 + h, v_0), \quad (u_0+h, v_0 +h), \quad (u_0, v_0 +h)

为顶点的的矩形AhkA_{hk}一对一的映射为φ(Ahk)\boldsymbol{\varphi}(A_{hk}),那么
lim⁡h→0,k→0σ(φ(Ahk))σ(Ahk)=∣∂(x,y)∂(u,v)∣(u0,v0) \lim \limits_{ h\to 0, k \to 0} \frac{\sigma(\boldsymbol{\varphi}(A_{hk}))}{\sigma(A_{hk})} = \left| \frac{\partial(x, y)}{\partial(u, v)} \right|_{(u_0, v_0)}

证:由映射后的图像有

σ(φ(Ahk))≈∥(φ(u0+h,v0)−φ(u0,v0))×(φ(u0,v0+k)−φ(u0,v0))∥ \sigma(\boldsymbol{\varphi}(A_{hk})) \approx \Vert (\boldsymbol{\varphi}(u_0 + h, v_0) - \boldsymbol{\varphi}(u_0, v_0)) \times (\boldsymbol{\varphi}(u_0, v_0 + k) - \boldsymbol{\varphi}(u_0, v_0)) \Vert

由于
φ(u0+h,v0)−φ(u0,v0)=∂φ∂u(u0,v0)h+ξφ(u0,v0+h)−φ(u0,v0)=∂φ∂v(u0,v0)h+η \begin{aligned} \boldsymbol{\varphi}(u_0 + h, v_0) - \boldsymbol{\varphi}(u_0, v_0) = \frac{\partial \boldsymbol{\varphi}}{\partial u} (u_0, v_0) h + \boldsymbol{\xi} \\ \boldsymbol{\varphi}(u_0, v_0+h) - \boldsymbol{\varphi}(u_0, v_0) = \frac{\partial \boldsymbol{\varphi}}{\partial v} (u_0, v_0) h + \boldsymbol{\eta} \end{aligned}

其中
∥ξ∥=o(h),∥η∥=o(k) \Vert \boldsymbol{\xi} \Vert = o(h), \quad \Vert \boldsymbol{\eta} \Vert = o(k)

所以
σ(φ(Ahk))≈∥∂φ∂u(u0,v0)×∂φ∂v(u0,v0)∥hk+o(hk) \sigma(\boldsymbol{\varphi}(A_{hk})) \approx \Vert \frac{\partial \boldsymbol{\varphi}}{\partial u}(u_0, v_0) \times \frac{\partial \boldsymbol{\varphi}}{\partial v}(u_0, v_0) \Vert hk + o(hk)

而
hk=σ(Ahk) hk = \sigma(A_{hk})

所以上式可以写为
σ(φ(Ahk))σ(Ahk)≈∥∂φ∂u(u0,v0)×∂φ∂v(u0,v0)∥+o(1) \frac{\sigma (\boldsymbol{\varphi}(A_{hk}))}{\sigma(A_{hk})} \approx \Vert \frac{\partial \boldsymbol{\varphi}}{\partial u}(u_0, v_0) \times \frac{\partial \boldsymbol{\varphi}}{\partial v}(u_0, v_0) \Vert + o(1)

又因为
∂φ∂u×∂φ∂v=∣ijk∂x∂u∂x∂v0∂y∂u∂y∂v0∣=∂(x,y)∂(u,v)k \frac{\partial \boldsymbol{\varphi}}{\partial u} \times \frac{\partial \boldsymbol{\varphi}}{\partial v} = \left| \begin{matrix} \boldsymbol{i} & \boldsymbol{j} & \boldsymbol{k}\\ \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} & 0 \\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} & 0 \\ \end{matrix} \right| = \frac{\partial(x, y)}{\partial(u, v)} \boldsymbol{k}

所以
∥∂φ∂u×∂φ∂v∥=∣∂(x,y)∂(u,v)∣ \Vert \frac{\partial \boldsymbol{\varphi}}{\partial u} \times \frac{\partial \boldsymbol{\varphi}}{\partial v} \Vert = \left| \frac{\partial (x, y)}{\partial(u, v)} \right|

从而
σ(φ(Ahk))σ(Ahk)≈∣∂(x,y)∂(u,v)∣(u0,v0)+o(hk) \frac{\sigma (\boldsymbol{\varphi}(A_{hk}))}{\sigma(A_{hk})} \approx \left| \frac{\partial(x, y)}{\partial(u, v)} \right|_{(u_0, v_0)} + o(hk)

所以当h→0,k→0h \to 0, k \to 0时,有
σ(φ(Ahk))σ(Ahk)=∣∂(x,y)∂(u,v)∣(u0,v0) \frac{\sigma (\boldsymbol{\varphi}(A_{hk}))}{\sigma(A_{hk})} = \left| \frac{\partial(x, y)}{\partial(u, v)} \right|_{(u_0, v_0)}

笔者注:感觉该证明方法不太好,约等于记号不应该用在这里,就算使用也应该讲明原因。

Q.E.D.

定理6

设R2\mathbb{R}^2中的有界闭区域DD有面积,函数F:→RF: \to \mathbb{R},映射

φ:{x=x(u,v)y=y(u,v)(u,v)∈Δ \boldsymbol{\varphi}: \left\{ \begin{aligned} x = x(u, v) \\ y = y(u, v) \end{aligned} \right. \quad (u, v) \in \Delta

是从Δ\Delta到DD上的正则映射,那么
∬DF(x,y)dxdy=∬ΔF∘φ(u,v)∣∂(x,y)∂(u,v)∣dudv \iint \limits_D F(x, y) \mathrm{d} x \mathrm{d} y = \iint \limits_\Delta F \circ \boldsymbol{\varphi}(u, v) \left| \frac{\partial(x, y)}{\partial(u, v)} \right| \mathrm{d} u \mathrm{d} v

证:用矩形II把uvuv平面上的闭区域Δ\Delta覆盖起来,用两族平行直线u=ui(i=0,1,2,⋯ ,m)u =u_i (i=0,1,2, \cdots, m)和v=vj(j=0,1,2,⋯ ,n)v = v_j (j=0,1,2, \cdots,n)分割II,其中

u0<u1<⋯<um−1<um,v0<v1<⋯<vn−1<vn u_0 < u_1 < \cdots < u_{m-1} < u_m, \quad v_0 < v_1 < \cdots < v_{n-1} < v_n

令Δui=ui−ui−1(i=1,2,⋯ ,m),Δvj=vj−vj−1(j=1,2,⋯ ,n)\Delta u_i = u_i - u_{i-1}(i=1,2,\cdots,m), \Delta v_j = v_j - v_{j-1}(j=1,2,\cdots,n)。分割后可以得到mnmn个矩形,在映射φ\boldsymbol{\varphi}的作用下,它们变成了xyxy平面下的mnmn个曲边平行四边形,这时只需要考虑那些完全被包含在DD中的曲边平行四边形,将它们记为Di(i=1,2,⋯ ,k)D_i(i=1,2,\cdots,k),即Di=φ(Δi)D_i=\boldsymbol{\varphi}(\Delta_i),其中Δi(i=1,2,⋯ ,k)\Delta_i(i=1,2,\cdots,k)是完全被包含在Δ\Delta中的矩形,任取一点ηi∈Di\boldsymbol{\eta}_i \in D_i ,并设Δi\Delta_i中唯一的点ξi\xi_i,使得φ(ξi)=ηi(i=1,2,⋯ ,k)\boldsymbol{\varphi} (\xi_i) = \eta_i(i=1,2,\cdots,k),作积分和
∑i=1kF(ηi)σ(Di) \sum_{i=1}^k F(\boldsymbol{\eta}_i) \sigma(D_i)

由定理5可知
σ(Di)∼∣detJφ(ξi)∣σ(Δi)(i=1,2,⋯ ,k) \sigma(D_i) \sim |det J \boldsymbol{\varphi}(\xi_i)| \sigma(\Delta_i) \quad (i=1,2,\cdots, k)

所以得到
∑i=1kF(ηi)σ(Di)∼∑i=1kF∘φ(ξi)∣detJφ(ξi)∣σ(Δi) \sum_{i=1}^k F(\boldsymbol{\eta}_i) \sigma(D_i) \sim \sum_{i=1}^k F \circ \boldsymbol{\varphi} (\boldsymbol{\xi}_i) |det J \boldsymbol{\varphi} (\boldsymbol{\xi}_i)| \sigma(\Delta_i)

当分割无限细时,由函数积分九的定理7可得,
∫DFdσ=∫ΔF∘φ∣Jφ∣dσ \int_D F \mathrm{d} \sigma = \int_\Delta F \circ \boldsymbol{\varphi} |J \boldsymbol{\varphi}| \mathrm{d} \sigma

Q.E.D.

定理7:极坐标换元

令x=rcos⁡θ,y=rsin⁡θx = r \cos \theta, y= r \sin \theta,则

∬DF(x,y)dxdy=∬ΔF(rcos⁡θ,rsin⁡θ)rdrdθ \iint \limits_D F(x, y) \mathrm{d} x \mathrm{d} y = \iint \limits_\Delta F(r \cos \theta, r \sin \theta) r \mathrm{d} r \mathrm{d} \theta

证:这时

∂(x,y)∂(u,v)=∣cos⁡θ−rsin⁡θsin⁡θrcos⁡θ∣=r \frac{\partial(x, y)}{\partial (u, v)} = \left| \begin{matrix} \cos \theta & -r\sin \theta \\ \sin \theta & r \cos \theta \end{matrix} \right| = r

代入定理6的结果即证。

Q.E.D.

继续这个系列