logo
导航

三重积分

三重积分的定义

三重积分的概念

设函数 f(x,y,z)f(x, y, z) 在有界闭区域 Ω\Omega 上有定义,将 Ω\Omega 任意分割成 nn 个小区域 ΔV1,ΔV2,,ΔVn\Delta V_1, \Delta V_2, \ldots, \Delta V_n,在每个小区域 ΔVi\Delta V_i 上任取一点 (ξi,ηi,ζi)(\xi_i, \eta_i, \zeta_i),作和式:

i=1nf(ξi,ηi,ζi)ΔVi\sum_{i=1}^n f(\xi_i, \eta_i, \zeta_i) \Delta V_i

其中 ΔVi\Delta V_i 表示小区域 ΔVi\Delta V_i 的体积。当分割越来越细,最大直径趋于零时,如果这个和式的极限存在,则称此极限为函数 f(x,y,z)f(x, y, z) 在区域 Ω\Omega 上的三重积分,记作:

Ωf(x,y,z)dV=limλ0i=1nf(ξi,ηi,ζi)ΔVi\iiint_\Omega f(x, y, z) dV = \lim_{\lambda \to 0} \sum_{i=1}^n f(\xi_i, \eta_i, \zeta_i) \Delta V_i

其中 λ\lambda 表示所有小区域的最大直径。

三重积分的几何意义

  • f(x,y,z)0f(x, y, z) \geq 0 时,三重积分表示以 Ω\Omega 为底,以曲面 w=f(x,y,z)w = f(x, y, z) 为顶的四维体积
  • f(x,y,z)<0f(x, y, z) < 0 时,三重积分的绝对值表示相应四维体积

三重积分的性质

线性性质

Ω[af(x,y,z)+bg(x,y,z)]dV=aΩf(x,y,z)dV+bΩg(x,y,z)dV\iiint_\Omega [af(x, y, z) + bg(x, y, z)] dV = a\iiint_\Omega f(x, y, z) dV + b\iiint_\Omega g(x, y, z) dV

区域可加性

如果区域 Ω\Omega 可以分解为两个不相重叠的子区域 Ω1\Omega_1Ω2\Omega_2,则:

Ωf(x,y,z)dV=Ω1f(x,y,z)dV+Ω2f(x,y,z)dV\iiint_\Omega f(x, y, z) dV = \iiint_{\Omega_1} f(x, y, z) dV + \iiint_{\Omega_2} f(x, y, z) dV

保号性

如果 f(x,y,z)0f(x, y, z) \geq 0Ω\Omega 上成立,则:

Ωf(x,y,z)dV0\iiint_\Omega f(x, y, z) dV \geq 0

单调性

如果 f(x,y,z)g(x,y,z)f(x, y, z) \geq g(x, y, z)Ω\Omega 上成立,则:

Ωf(x,y,z)dVΩg(x,y,z)dV\iiint_\Omega f(x, y, z) dV \geq \iiint_\Omega g(x, y, z) dV

中值定理

如果 f(x,y,z)f(x, y, z) 在闭区域 Ω\Omega 上连续,则存在点 (ξ,η,ζ)Ω(\xi, \eta, \zeta) \in \Omega,使得:

Ωf(x,y,z)dV=f(ξ,η,ζ)V(Ω)\iiint_\Omega f(x, y, z) dV = f(\xi, \eta, \zeta) \cdot V(\Omega)

其中 V(Ω)V(\Omega) 表示区域 Ω\Omega 的体积。

三重积分的计算

直角坐标系下的计算

累次积分法

如果区域 Ω\Omega 可以表示为:

Ω={(x,y,z)axb,ϕ1(x)yϕ2(x),ψ1(x,y)zψ2(x,y)}\Omega = \{(x, y, z) | a \leq x \leq b, \phi_1(x) \leq y \leq \phi_2(x), \psi_1(x, y) \leq z \leq \psi_2(x, y)\}

则:

Ωf(x,y,z)dV=abϕ1(x)ϕ2(x)ψ1(x,y)ψ2(x,y)f(x,y,z)dzdydx\iiint_\Omega f(x, y, z) dV = \int_a^b \int_{\phi_1(x)}^{\phi_2(x)} \int_{\psi_1(x, y)}^{\psi_2(x, y)} f(x, y, z) dz dy dx

柱面坐标系下的计算

在柱面坐标系下,三重积分的计算公式为:

Ωf(x,y,z)dV=Ωf(rcosθ,rsinθ,z)rdrdθdz\iiint_\Omega f(x, y, z) dV = \iiint_\Omega f(r\cos\theta, r\sin\theta, z) r dr d\theta dz

其中 rr 为极径,θ\theta 为极角,zz 为高度。

柱面坐标下的累次积分

如果区域 Ω\Omega 在柱面坐标系下可以表示为:

Ω={(r,θ,z)αθβ,r1(θ)rr2(θ),z1(r,θ)zz2(r,θ)}\Omega = \{(r, \theta, z) | \alpha \leq \theta \leq \beta, r_1(\theta) \leq r \leq r_2(\theta), z_1(r, \theta) \leq z \leq z_2(r, \theta)\}

则:

Ωf(x,y,z)dV=αβr1(θ)r2(θ)z1(r,θ)z2(r,θ)f(rcosθ,rsinθ,z)rdzdrdθ\iiint_\Omega f(x, y, z) dV = \int_\alpha^\beta \int_{r_1(\theta)}^{r_2(\theta)} \int_{z_1(r, \theta)}^{z_2(r, \theta)} f(r\cos\theta, r\sin\theta, z) r dz dr d\theta

球面坐标系下的计算

在球面坐标系下,三重积分的计算公式为:

Ωf(x,y,z)dV=Ωf(ρsinϕcosθ,ρsinϕsinθ,ρcosϕ)ρ2sinϕdρdϕdθ\iiint_\Omega f(x, y, z) dV = \iiint_\Omega f(\rho\sin\phi\cos\theta, \rho\sin\phi\sin\theta, \rho\cos\phi) \rho^2\sin\phi d\rho d\phi d\theta

其中 ρ\rho 为球半径,ϕ\phi 为天顶角,θ\theta 为方位角。

球面坐标下的累次积分

如果区域 Ω\Omega 在球面坐标系下可以表示为:

Ω={(ρ,ϕ,θ)αθβ,ϕ1(θ)ϕϕ2(θ),ρ1(ϕ,θ)ρρ2(ϕ,θ)}\Omega = \{(\rho, \phi, \theta) | \alpha \leq \theta \leq \beta, \phi_1(\theta) \leq \phi \leq \phi_2(\theta), \rho_1(\phi, \theta) \leq \rho \leq \rho_2(\phi, \theta)\}

则:

Ωf(x,y,z)dV=αβϕ1(θ)ϕ2(θ)ρ1(ϕ,θ)ρ2(ϕ,θ)f(ρsinϕcosθ,ρsinϕsinθ,ρcosϕ)ρ2sinϕdρdϕdθ\iiint_\Omega f(x, y, z) dV = \int_\alpha^\beta \int_{\phi_1(\theta)}^{\phi_2(\theta)} \int_{\rho_1(\phi, \theta)}^{\rho_2(\phi, \theta)} f(\rho\sin\phi\cos\theta, \rho\sin\phi\sin\theta, \rho\cos\phi) \rho^2\sin\phi d\rho d\phi d\theta

三重积分的应用

体积计算

区域 Ω\Omega 的体积为:

V(Ω)=ΩdVV(\Omega) = \iiint_\Omega dV

质量计算

如果 ρ(x,y,z)\rho(x, y, z) 表示体密度,则区域 Ω\Omega 的质量为:

M=Ωρ(x,y,z)dVM = \iiint_\Omega \rho(x, y, z) dV

重心计算

区域 Ω\Omega 的重心坐标为:

xˉ=1V(Ω)ΩxdV,yˉ=1V(Ω)ΩydV,zˉ=1V(Ω)ΩzdV\bar{x} = \frac{1}{V(\Omega)} \iiint_\Omega x dV, \quad \bar{y} = \frac{1}{V(\Omega)} \iiint_\Omega y dV, \quad \bar{z} = \frac{1}{V(\Omega)} \iiint_\Omega z dV

转动惯量计算

区域 Ω\Omega 关于 zz 轴的转动惯量为:

Iz=Ω(x2+y2)ρ(x,y,z)dVI_z = \iiint_\Omega (x^2 + y^2) \rho(x, y, z) dV

例题

例 1:直角坐标系下的三重积分

计算 ΩzdV\iiint_\Omega z dV,其中 Ω\Omega 为长方体 0x1,0y1,0z20 \leq x \leq 1, 0 \leq y \leq 1, 0 \leq z \leq 2

ΩzdV=010102zdzdydx\iiint_\Omega z dV = \int_0^1 \int_0^1 \int_0^2 z dz dy dx

=0101[12z2]02dydx= \int_0^1 \int_0^1 \left[\frac{1}{2}z^2\right]_0^2 dy dx

=01012dydx= \int_0^1 \int_0^1 2 dy dx

=012dx=2= \int_0^1 2 dx = 2

例 2:柱面坐标系下的三重积分

计算 Ωx2+y2dV\iiint_\Omega \sqrt{x^2 + y^2} dV,其中 Ω\Omega 为圆柱体 x2+y21,0z1x^2 + y^2 \leq 1, 0 \leq z \leq 1

: 在柱面坐标系下,Ω={(r,θ,z)0r1,0θ2π,0z1}\Omega = \{(r, \theta, z) | 0 \leq r \leq 1, 0 \leq \theta \leq 2\pi, 0 \leq z \leq 1\}

Ωx2+y2dV=02π0101rrdzdrdθ\iiint_\Omega \sqrt{x^2 + y^2} dV = \int_0^{2\pi} \int_0^1 \int_0^1 r \cdot r dz dr d\theta

=02π01r2drdθ= \int_0^{2\pi} \int_0^1 r^2 dr d\theta

=2π01r2dr=2π×13=2π3= 2\pi \int_0^1 r^2 dr = 2\pi \times \frac{1}{3} = \frac{2\pi}{3}

例 3:球面坐标系下的三重积分

计算 Ω(x2+y2+z2)dV\iiint_\Omega (x^2 + y^2 + z^2) dV,其中 Ω\Omega 为单位球 x2+y2+z21x^2 + y^2 + z^2 \leq 1

: 在球面坐标系下,Ω={(ρ,ϕ,θ)0ρ1,0ϕπ,0θ2π}\Omega = \{(\rho, \phi, \theta) | 0 \leq \rho \leq 1, 0 \leq \phi \leq \pi, 0 \leq \theta \leq 2\pi\}

Ω(x2+y2+z2)dV=02π0π01ρ2ρ2sinϕdρdϕdθ\iiint_\Omega (x^2 + y^2 + z^2) dV = \int_0^{2\pi} \int_0^\pi \int_0^1 \rho^2 \cdot \rho^2\sin\phi d\rho d\phi d\theta

=02π0π01ρ4sinϕdρdϕdθ= \int_0^{2\pi} \int_0^\pi \int_0^1 \rho^4\sin\phi d\rho d\phi d\theta

=02π0π[15ρ5]01sinϕdϕdθ= \int_0^{2\pi} \int_0^\pi \left[\frac{1}{5}\rho^5\right]_0^1\sin\phi d\phi d\theta

=02π0π15sinϕdϕdθ= \int_0^{2\pi} \int_0^\pi \frac{1}{5}\sin\phi d\phi d\theta

=02π15[cosϕ]0πdθ= \int_0^{2\pi} \frac{1}{5} \left[-\cos\phi\right]_0^\pi d\theta

=02π25dθ=4π5= \int_0^{2\pi} \frac{2}{5} d\theta = \frac{4\pi}{5}

练习题

练习 1

计算 ΩxyzdV\iiint_\Omega xyz dV,其中 Ω\Omega0x1,0y1,0z10 \leq x \leq 1, 0 \leq y \leq 1, 0 \leq z \leq 1

练习 2

计算 Ω(x+y+z)dV\iiint_\Omega (x + y + z) dV,其中 Ω\Omegax2+y21,0z1x^2 + y^2 \leq 1, 0 \leq z \leq 1

练习 3

计算 Ωex2+y2+z2dV\iiint_\Omega e^{x^2 + y^2 + z^2} dV,其中 Ω\Omegax2+y2+z21x^2 + y^2 + z^2 \leq 1

练习 4

求球体 x2+y2+z2R2x^2 + y^2 + z^2 \leq R^2 的体积。

练习 5

计算 ΩzdV\iiint_\Omega z dV,其中 Ω\Omega 为锥体 x2+y2z2,0z1x^2 + y^2 \leq z^2, 0 \leq z \leq 1

参考答案

1. 计算 ΩxyzdV\iiint_\Omega xyz dV,其中 Ω\Omega0x1,0y1,0z10 \leq x \leq 1, 0 \leq y \leq 1, 0 \leq z \leq 1

ΩxyzdV=010101xyzdzdydx\iiint_\Omega xyz dV = \int_0^1 \int_0^1 \int_0^1 xyz dz dy dx

=0101xy[12z2]01dydx= \int_0^1 \int_0^1 xy \left[\frac{1}{2}z^2\right]_0^1 dy dx

=0101xy×12dydx= \int_0^1 \int_0^1 xy \times \frac{1}{2} dy dx

=1201x[12y2]01dx= \frac{1}{2} \int_0^1 x \left[\frac{1}{2}y^2\right]_0^1 dx

=1401xdx=18= \frac{1}{4} \int_0^1 x dx = \frac{1}{8}


2. 计算 Ω(x+y+z)dV\iiint_\Omega (x + y + z) dV,其中 Ω\Omegax2+y21,0z1x^2 + y^2 \leq 1, 0 \leq z \leq 1

在柱面坐标系下:Ω={(r,θ,z)0r1,0θ2π,0z1}\Omega = \{(r, \theta, z) | 0 \leq r \leq 1, 0 \leq \theta \leq 2\pi, 0 \leq z \leq 1\}

Ω(x+y+z)dV=02π0101(rcosθ+rsinθ+z)rdzdrdθ\iiint_\Omega (x + y + z) dV = \int_0^{2\pi} \int_0^1 \int_0^1 (r\cos\theta + r\sin\theta + z) r dz dr d\theta

=02π0101(r2cosθ+r2sinθ+rz)dzdrdθ= \int_0^{2\pi} \int_0^1 \int_0^1 (r^2\cos\theta + r^2\sin\theta + rz) dz dr d\theta

=02π01(r2cosθ+r2sinθ+12r)drdθ= \int_0^{2\pi} \int_0^1 (r^2\cos\theta + r^2\sin\theta + \frac{1}{2}r) dr d\theta

=02π(13cosθ+13sinθ+14)dθ= \int_0^{2\pi} \left(\frac{1}{3}\cos\theta + \frac{1}{3}\sin\theta + \frac{1}{4}\right) d\theta

=02π14dθ=π2= \int_0^{2\pi} \frac{1}{4} d\theta = \frac{\pi}{2}


3. 计算 Ωex2+y2+z2dV\iiint_\Omega e^{x^2 + y^2 + z^2} dV,其中 Ω\Omegax2+y2+z21x^2 + y^2 + z^2 \leq 1

在球面坐标系下:Ω={(ρ,ϕ,θ)0ρ1,0ϕπ,0θ2π}\Omega = \{(\rho, \phi, \theta) | 0 \leq \rho \leq 1, 0 \leq \phi \leq \pi, 0 \leq \theta \leq 2\pi\}

Ωex2+y2+z2dV=02π0π01eρ2ρ2sinϕdρdϕdθ\iiint_\Omega e^{x^2 + y^2 + z^2} dV = \int_0^{2\pi} \int_0^\pi \int_0^1 e^{\rho^2} \rho^2\sin\phi d\rho d\phi d\theta

=2π0π01eρ2ρ2sinϕdρdϕ= 2\pi \int_0^\pi \int_0^1 e^{\rho^2} \rho^2\sin\phi d\rho d\phi

=2π0πsinϕdϕ01eρ2ρ2dρ= 2\pi \int_0^\pi \sin\phi d\phi \int_0^1 e^{\rho^2} \rho^2 d\rho

=2π×2×01eρ2ρ2dρ= 2\pi \times 2 \times \int_0^1 e^{\rho^2} \rho^2 d\rho

=4π×1201eρ2d(ρ2)=2π(e1)= 4\pi \times \frac{1}{2} \int_0^1 e^{\rho^2} d(\rho^2) = 2\pi(e - 1)


4. 求球体 x2+y2+z2R2x^2 + y^2 + z^2 \leq R^2 的体积

在球面坐标系下:Ω={(ρ,ϕ,θ)0ρR,0ϕπ,0θ2π}\Omega = \{(\rho, \phi, \theta) | 0 \leq \rho \leq R, 0 \leq \phi \leq \pi, 0 \leq \theta \leq 2\pi\}

V=ΩdV=02π0π0Rρ2sinϕdρdϕdθV = \iiint_\Omega dV = \int_0^{2\pi} \int_0^\pi \int_0^R \rho^2\sin\phi d\rho d\phi d\theta

=2π0πsinϕdϕ0Rρ2dρ= 2\pi \int_0^\pi \sin\phi d\phi \int_0^R \rho^2 d\rho

=2π×2×R33=4πR33= 2\pi \times 2 \times \frac{R^3}{3} = \frac{4\pi R^3}{3}


5. 计算 ΩzdV\iiint_\Omega z dV,其中 Ω\Omega 为锥体 x2+y2z2,0z1x^2 + y^2 \leq z^2, 0 \leq z \leq 1

在柱面坐标系下:Ω={(r,θ,z)0rz,0θ2π,0z1}\Omega = \{(r, \theta, z) | 0 \leq r \leq z, 0 \leq \theta \leq 2\pi, 0 \leq z \leq 1\}

ΩzdV=02π010zzrdrdzdθ\iiint_\Omega z dV = \int_0^{2\pi} \int_0^1 \int_0^z z \cdot r dr dz d\theta

=2π01z0zrdrdz= 2\pi \int_0^1 z \int_0^z r dr dz

=2π01z[12r2]0zdz= 2\pi \int_0^1 z \left[\frac{1}{2}r^2\right]_0^z dz

=2π01z×12z2dz= 2\pi \int_0^1 z \times \frac{1}{2}z^2 dz

=π01z3dz=π×14=π4= \pi \int_0^1 z^3 dz = \pi \times \frac{1}{4} = \frac{\pi}{4}

常见错误与注意事项

  1. 坐标变换:不同坐标系下的积分元素不同
  2. 积分次序:要注意积分变量的上下限对应关系
  3. 区域边界:要正确确定积分区域的边界
  4. 对称性:利用对称性可以简化计算
  5. 体积元素:柱面坐标下为 rdrdθdzr dr d\theta dz,球面坐标下为 ρ2sinϕdρdϕdθ\rho^2\sin\phi d\rho d\phi d\theta

提示:三重积分的计算要特别注意坐标变换和积分区域的确定,建议多做练习来熟练掌握各种计算方法。