MathLabs

分析学

重积分

在二维或更高维区域上的积分,用于计算体积与高维总量。

直观从面积窄条到体积柱体

一元积分 ∫abf(x) dx\int_a^b f(x)\,dx 把区间切成宽度为 Δx\Delta x 的小段并累加矩形面积。二重积分 ∬Df(x,y) dA\iint_D f(x,y)\,dA 将完全相同的思想推广到二维区域 DD:用面积为 ΔA=Δx Δy\Delta A = \Delta x\,\Delta y 的微小矩形铺满 DD,在每个小格上竖起高为 f(xi,yj)f(x_i,y_j) 的细柱体,再令网格不断加密。所有柱体体积之和就给出曲面 z=f(x,y)z = f(x,y) 下方的有向体积,或者面密度为 f(x,y)f(x,y) 的薄板的总质量。

表示二维区域上二重积分所累积体积的交互式三维曲面图。
平面区域 DD 上的曲面 z=f(x,y)z = f(x,y):二重积分 ∬Df(x,y) dA\iint_D f(x,y)\,dA 累加曲面下方竖直柱体的体积,而富比尼定理通过逐片积分来计算它。

大学严格定义:二重与三重积分、富比尼累次化简与雅可比行列式

定义: 二重与三重黎曼积分

对于有界可测区域 DD 上的有界函数 f:D⊆R2→Rf : D \subseteq \mathbb{R}^2 \to \mathbb{R},二重积分 ∬Df(x,y) dA\iint_D f(x,y)\,dA 是当分割的网格直径趋于 00 时黎曼和 ∑i,jf(xi∗,yj∗) Δxi Δyj\sum_{i,j} f(x_i^*, y_j^*)\,\Delta x_i\,\Delta y_j 的极限。类似地,在立体区域 E⊆R3E \subseteq \mathbb{R}^3 上,三重积分 ∭Ef(x,y,z) dV\iiint_E f(x,y,z)\,dV 对体积为 ΔV=Δx Δy Δz\Delta V = \Delta x\,\Delta y\,\Delta z 的小长方体求和取极限。

∬Df(x,y) dA=∫ab(∫g1(x)g2(x)f(x,y) dy)dx=∫cd(∫h1(y)h2(y)f(x,y) dx)dy\iint_D f(x,y)\,dA = \int_a^b \left(\int_{g_1(x)}^{g_2(x)} f(x,y)\,dy\right) dx = \int_c^d \left(\int_{h_1(y)}^{h_2(y)} f(x,y)\,dx\right) dy

当积分区域呈圆形、柱形或球形时,光滑坐标变换 Φ(u,v)=(x(u,v),y(u,v))\Phi(u,v) = (x(u,v), y(u,v)) 能显著简化积分上下限。在映射 Φ\Phi 下,微小的 du×dvdu \times dv 矩形被拉伸并剪切成平行四边形,其面积缩放因子为雅可比行列式的绝对值 ∣det⁡DΦ∣=∣∂(x,y)∂(u,v)∣|\det D\Phi| = \left|\frac{\partial(x,y)}{\partial(u,v)}\right|。在极坐标 (x,y)=(rcos⁡θ,rsin⁡θ)(x,y) = (r\cos\theta, r\sin\theta) 中该因子为 rr,而在球坐标 (x,y,z)=(ρsin⁡φcos⁡θ,ρsin⁡φsin⁡θ,ρcos⁡φ)(x,y,z) = (\rho\sin\varphi\cos\theta, \rho\sin\varphi\sin\theta, \rho\cos\varphi) 中该因子为 ρ2sin⁡φ\rho^2\sin\varphi。

∬Φ(U)f(x,y) dx dy=∬Uf(Φ(u,v)) ∣det⁡DΦ(u,v)∣ du dv,∫−∞∞e−x2 dx=π\iint_{\Phi(U)} f(x,y)\,dx\,dy = \iint_U f(\Phi(u,v))\,|\det D\Phi(u,v)|\,du\,dv, \qquad \int_{-\infty}^{\infty} e^{-x^2}\,dx = \sqrt{\pi}
二维与三维坐标系及雅可比面积/体积微元
坐标系变换公式面积 / 体积微元
二维直角坐标(x,y)(x,y)dA=dx dydA = dx\,dy
二维极坐标x=rcos⁡θ,  y=rsin⁡θx = r\cos\theta,\; y = r\sin\thetadA=r dr dθdA = r\,dr\,d\theta
三维柱坐标x=rcos⁡θ,  y=rsin⁡θ,  z=zx = r\cos\theta,\; y = r\sin\theta,\; z = zdV=r dr dθ dzdV = r\,dr\,d\theta\,dz
三维球坐标x=ρsin⁡φcos⁡θ,  y=ρsin⁡φsin⁡θ,  z=ρcos⁡φx = \rho\sin\varphi\cos\theta,\; y = \rho\sin\varphi\sin\theta,\; z = \rho\cos\varphidV=ρ2sin⁡φ dρ dφ dθdV = \rho^2\sin\varphi\,d\rho\,d\varphi\,d\theta

大学核心定理:富比尼定理与雅可比换元公式

若 f(x,y)f(x,y) 在矩形 R=[a,b]×[c,d]R = [a,b] \times [c,d] 上连续,则二重积分等于两个一元累次积分:∬Rf(x,y) dA=∫ab(∫cdf(x,y) dy)dx=∫cd(∫abf(x,y) dx)dy\iint_R f(x,y)\,dA = \int_a^b \left(\int_c^d f(x,y)\,dy\right) dx = \int_c^d \left(\int_a^b f(x,y)\,dx\right) dy。更一般地,在纵线简单区域 D={(x,y):a≤x≤b,  g1(x)≤y≤g2(x)}D = \{(x,y) : a \le x \le b,\; g_1(x) \le y \le g_2(x)\} 上有 ∬Df(x,y) dA=∫ab∫g1(x)g2(x)f(x,y) dy dx\iint_D f(x,y)\,dA = \int_a^b \int_{g_1(x)}^{g_2(x)} f(x,y)\,dy\,dx。

为什么成立?

通过累加网格中的小柱体求体积时,无论是先沿 yy 方向把每一列加起来得到截面面积 A(x)A(x) 再将 A(x)A(x) 对 xx 积分,还是先垂直于 yy 轴切片再积分,所得结果完全一致。

证明

将 [a,b][a,b] 分为 mm 个宽为 Δx\Delta x 的子区间 [xi−1,xi][x_{i-1}, x_i],将 [c,d][c,d] 分为 nn 个宽为 Δy\Delta y 的子区间 [yj−1,yj][y_{j-1}, y_j]。对每个固定的 xx,定义截面积分 A(x)=∫cdf(x,y) dyA(x) = \int_c^d f(x,y)\,dy。由积分中值定理,在每个小段 [yj−1,yj][y_{j-1}, y_j] 上存在 yij∗∈[yj−1,yj]y_{ij}^* \in [y_{j-1}, y_j] 使得 ∫yj−1yjf(xi,y) dy=f(xi,yij∗) Δy\int_{y_{j-1}}^{y_j} f(x_i, y)\,dy = f(x_i, y_{ij}^*)\,\Delta y。

对 j=1,…,nj = 1, \dots, n 求和得 A(xi)=∑j=1nf(xi,yij∗) ΔyA(x_i) = \sum_{j=1}^n f(x_i, y_{ij}^*)\,\Delta y。再乘以 Δx\Delta x 并对 i=1,…,mi = 1, \dots, m 求和,得到 ∑i=1mA(xi) Δx=∑i=1m∑j=1nf(xi,yij∗) Δx Δy\sum_{i=1}^m A(x_i)\,\Delta x = \sum_{i=1}^m \sum_{j=1}^n f(x_i, y_{ij}^*)\,\Delta x\,\Delta y。因为 ff 在紧矩形 RR 上一致连续,令 Δx,Δy→0\Delta x, \Delta y \to 0,左边收敛到 ∫abA(x) dx=∫ab(∫cdf(x,y) dy)dx\int_a^b A(x)\,dx = \int_a^b \left(\int_c^d f(x,y)\,dy\right) dx,右边收敛到 ∬Rf(x,y) dA\iint_R f(x,y)\,dA。交换 xx 与 yy 的角色重复上述论证,即得与另一积分顺序的相等关系。

设 Φ:U→Φ(U)⊆R2\Phi : U \to \Phi(U) \subseteq \mathbb{R}^2 为 C1C^1 微分同胚,其雅可比矩阵为 DΦ(u,v)=(xuxvyuyv)D\Phi(u,v) = \begin{pmatrix} x_u & x_v \\ y_u & y_v \end{pmatrix}。对 Φ(U)\Phi(U) 上的任意可积函数 ff,有 ∬Φ(U)f(x,y) dx dy=∬Uf(x(u,v),y(u,v)) ∣det⁡DΦ(u,v)∣ du dv\iint_{\Phi(U)} f(x,y)\,dx\,dy = \iint_U f(x(u,v), y(u,v))\,|\det D\Phi(u,v)|\,du\,dv,其中 det⁡DΦ=xuyv−xvyu\det D\Phi = x_u y_v - x_v y_u。

为什么成立?

正如一元换元法中 dx=g′(u) dudx = g'(u)\,du 度量长度的伸缩一样,∣det⁡DΦ(u,v)∣|\det D\Phi(u,v)| 度量了微小 (u,v)(u,v) 矩形被映射为 (x,y)(x,y) 平面上的平行四边形时的局部面积缩放率。

证明

考虑 UU 中的小矩形 Rij=[ui,ui+Δu]×[vj,vj+Δv]R_{ij} = [u_i, u_i + \Delta u] \times [v_j, v_j + \Delta v]。利用 (ui,vj)(u_i, v_j) 处的一阶泰勒展开,两条边 (Δu,0)( \Delta u, 0 ) 与 (0,Δv)( 0, \Delta v ) 近似映射为切向量 a=Φu(ui,vj) Δu=(xu,yu) Δu\mathbf{a} = \Phi_u(u_i,v_j)\,\Delta u = (x_u, y_u)\,\Delta u 和 b=Φv(ui,vj) Δv=(xv,yv) Δv\mathbf{b} = \Phi_v(u_i,v_j)\,\Delta v = (x_v, y_v)\,\Delta v。在 R2\mathbb{R}^2 中由 a\mathbf{a} 与 b\mathbf{b} 张成的平行四边形面积等于 ∣xuyv−xvyu∣ Δu Δv=∣det⁡DΦ(ui,vj)∣ Δu Δv|x_u y_v - x_v y_u|\,\Delta u\,\Delta v = |\det D\Phi(u_i, v_j)|\,\Delta u\,\Delta v(相差高阶无穷小 o(Δu Δv)o(\Delta u\,\Delta v))。

将 ΔAij≈∣det⁡DΦ(ui,vj)∣ Δu Δv\Delta A_{ij} \approx |\det D\Phi(u_i, v_j)|\,\Delta u\,\Delta v 代入黎曼和 ∑i,jf(Φ(ui,vj)) ΔAij\sum_{i,j} f(\Phi(u_i, v_j))\,\Delta A_{ij},即得 ∑i,jf(Φ(ui,vj)) ∣det⁡DΦ(ui,vj)∣ Δu Δv\sum_{i,j} f(\Phi(u_i, v_j))\,|\det D\Phi(u_i, v_j)|\,\Delta u\,\Delta v。特别地,对极坐标 x=rcos⁡θx = r\cos\theta,y=rsin⁡θy = r\sin\theta,有 det⁡DΦ=(cos⁡θ)(rcos⁡θ)−(−rsin⁡θ)(sin⁡θ)=r(cos⁡2θ+sin⁡2θ)=r\det D\Phi = (\cos\theta)(r\cos\theta) - (-r\sin\theta)(\sin\theta) = r(\cos^2\theta + \sin^2\theta) = r,从而给出 dx dy=r dr dθdx\,dy = r\,dr\,d\theta。

大学实际应用与典型例题

重积分广泛出现在物理学、概率论与工程学中:在力学里,三重积分用于计算总质量 M=∭Eρ(x,y,z) dVM = \iiint_E \rho(x,y,z)\,dV、质心以及旋转机械的转动惯量 Iz=∭E(x2+y2)ρ dVI_z = \iiint_E (x^2+y^2)\rho\,dV;在概率论与数理统计中,正态分布的归一化依赖于高斯积分 ∫−∞∞e−x2 dx=π\int_{-\infty}^{\infty} e^{-x^2}\,dx = \sqrt{\pi},该结果正是通过平方化为极坐标下的二重积分证明的;而在天体物理与电磁学中,球坐标微元 ρ2sin⁡φ dρ dφ dθ\rho^2\sin\varphi\,d\rho\,d\varphi\,d\theta 被用来对恒星与行星上的引力场及电荷分布进行积分。

例题: 当内层积分没有初等原函数时交换积分次序

画出积分区域并利用富比尼定理交换积分次序,计算累次积分 I=∫01∫y1ex2 dx dyI = \int_0^1 \int_y^1 e^{x^2}\,dx\,dy。

解答

由于 ex2e^{x^2} 没有初等原函数,内层积分 ∫y1ex2 dx\int_y^1 e^{x^2}\,dx 无法直接求出解析表达式。但不等式 0≤y≤10 \le y \le 1 与 y≤x≤1y \le x \le 1 描述的是三角形区域 D={(x,y):0≤x≤1,  0≤y≤x}D = \{(x,y) : 0 \le x \le 1,\; 0 \le y \le x\}。

根据富比尼定理,先沿竖直方向切片可将积分改写为 I=∫01(∫0xex2 dy)dx=∫01xex2 dxI = \int_0^1 \left(\int_0^x e^{x^2}\,dy\right) dx = \int_0^1 x e^{x^2}\,dx。再令 u=x2u = x^2,du=2x dxdu = 2x\,dx,立得 I=[12ex2]01=e−12I = \left[\frac{1}{2}e^{x^2}\right]_0^1 = \frac{e - 1}{2}。

例题: 利用极坐标计算高斯积分

通过将 I2I^2 表示为 R2\mathbb{R}^2 上的二重积分并转换为极坐标,证明高斯积分公式 I=∫−∞∞e−x2 dx=πI = \int_{-\infty}^{\infty} e^{-x^2}\,dx = \sqrt{\pi}。

解答

用哑变量 xx 和 yy 将两个独立的 II 相乘:I2=(∫−∞∞e−x2 dx)(∫−∞∞e−y2 dy)=∬R2e−(x2+y2) dx dyI^2 = \left(\int_{-\infty}^{\infty} e^{-x^2}\,dx\right)\left(\int_{-\infty}^{\infty} e^{-y^2}\,dy\right) = \iint_{\mathbb{R}^2} e^{-(x^2+y^2)}\,dx\,dy。

换为极坐标 x=rcos⁡θx = r\cos\theta,y=rsin⁡θy = r\sin\theta,其中 0≤r<∞0 \le r < \infty,0≤θ≤2π0 \le \theta \le 2\pi,雅可比因子给出 dx dy=r dr dθdx\,dy = r\,dr\,d\theta。于是 I2=∫02π∫0∞e−r2r dr dθ=2π[−12e−r2]0∞=2π⋅12=πI^2 = \int_0^{2\pi} \int_0^{\infty} e^{-r^2} r\,dr\,d\theta = 2\pi \left[-\frac{1}{2}e^{-r^2}\right]_0^{\infty} = 2\pi \cdot \frac{1}{2} = \pi。因为 e−x2>0e^{-x^2} > 0,取算术平方根即得 I=πI = \sqrt{\pi}。

计算矩形区域 R=[0,2]×[0,1]R = [0, 2] \times [0, 1] 上的二重积分 ∬R(2x+4y) dA\iint_R (2x + 4y)\,dA。

交换累次积分 ∫02∫x24f(x,y) dy dx\int_0^2 \int_{x^2}^4 f(x,y)\,dy\,dx 的积分次序后,得到下列哪个累次积分?

在三维球坐标 (x,y,z)=(ρsin⁡φcos⁡θ,ρsin⁡φsin⁡θ,ρcos⁡φ)(x,y,z) = (\rho\sin\varphi\cos\theta, \rho\sin\varphi\sin\theta, \rho\cos\varphi) 中,体积微元 dVdV 是什么?

利用高斯积分 ∫−∞∞e−x2 dx=π\int_{-\infty}^{\infty} e^{-x^2}\,dx = \sqrt{\pi},二维概率积分 ∬R2e−(x2+y2)/2 dx dy\iint_{\mathbb{R}^2} e^{-(x^2+y^2)/2}\,dx\,dy 的值是多少?

参考文献

  1. Jerrold E. Marsden, Anthony J. Tromba (2012). Vector Calculus
  2. Tom M. Apostol (1974). Mathematical Analysis
  3. Tom M. Apostol (1969). Calculus, Vol. 2: Multi-Variable Calculus and Linear Algebra with Applications