球坐标下的三重积分计算

为什么需要球坐标

当积分区域是球体、球锥、半球等球形区域,或被积函数含有 x2+y2+z2x^2 + y^2 + z^2 时,球坐标可以极大简化计算。球坐标用”到原点的距离 + 两个角度”来定位空间中的点,是处理球对称问题的自然选择。

球坐标变换

球坐标变换公式

直角坐标与球坐标的关系: {x=rsin⁡ϕcos⁡θy=rsin⁡ϕsin⁡θz=rcos⁡ϕ\begin{cases} x = r \sin\phi \cos\theta \\ y = r \sin\phi \sin\theta \\ z = r \cos\phi \end{cases}

其中:

  • r≥0r \geq 0 是点到原点的距离(径向距离)
  • ϕ∈[0,π]\phi \in [0, \pi] 是点与原点连线和 zz 轴正方向的夹角(极角/余纬角)
  • θ∈[0,2π)\theta \in [0, 2\pi) 是点在 xOyxOy 平面上的投影与 xx 轴正方向的夹角(方位角)

雅可比行列式: J=∂(x,y,z)∂(r,ϕ,θ)=r2sin⁡ϕJ = \frac{\partial(x,y,z)}{\partial(r,\phi,\theta)} = r^2 \sin\phi

因此体积元素:dv=dx dy dz=r2sin⁡ϕ dr dϕ dθdv = dx\,dy\,dz = r^2 \sin\phi \, dr \, d\phi \, d\theta

角度范围的几何意义

  • ϕ=0\phi = 0:zz 轴正方向
  • ϕ=π2\phi = \frac{\pi}{2}:xOyxOy 平面
  • ϕ=π\phi = \pi:zz 轴负方向
  • θ=0\theta = 0:xOzxOz 平面的右半部分(y=0,x>0y=0, x>0)
  • θ=π2\theta = \frac{\pi}{2}:yOzyOz 平面的前半部分(x=0,y>0x=0, y>0)

球坐标下的三重积分公式

∭Ωf(x,y,z) dv=∭Ω′f(rsin⁡ϕcos⁡θ,rsin⁡ϕsin⁡θ,rcos⁡ϕ)⋅r2sin⁡ϕ dr dϕ dθ\iiint_\Omega f(x,y,z) \, dv = \iiint_{\Omega'} f(r\sin\phi\cos\theta, r\sin\phi\sin\theta, r\cos\phi) \cdot r^2 \sin\phi \, dr \, d\phi \, d\theta

常见球坐标区域

1. 球体(球心在原点)

球体

Ω:x2+y2+z2≤R2\Omega: x^2 + y^2 + z^2 \leq R^2

球坐标下: Ω′:0≤r≤R,0≤ϕ≤π,0≤θ≤2π\Omega': 0 \leq r \leq R, \quad 0 \leq \phi \leq \pi, \quad 0 \leq \theta \leq 2\pi

∭Ωf dv=∫02πdθ∫0πsin⁡ϕ dϕ∫0Rf(rsin⁡ϕcos⁡θ,…)r2 dr\iiint_\Omega f \, dv = \int_0^{2\pi} d\theta \int_0^\pi \sin\phi \, d\phi \int_0^R f(r\sin\phi\cos\theta, \ldots) r^2 \, dr

2. 上半球体

上半球体

Ω:x2+y2+z2≤R2,z≥0\Omega: x^2 + y^2 + z^2 \leq R^2, \quad z \geq 0

球坐标下: Ω′:0≤r≤R,0≤ϕ≤π2,0≤θ≤2π\Omega': 0 \leq r \leq R, \quad 0 \leq \phi \leq \frac{\pi}{2}, \quad 0 \leq \theta \leq 2\pi

3. 球锥区域

球锥区域

由球面 x2+y2+z2=R2x^2+y^2+z^2 = R^2 和圆锥面 z=x2+y2z = \sqrt{x^2+y^2} 围成的区域(圆锥内部、球面内部):

圆锥面 z=x2+y2z = \sqrt{x^2+y^2} 对应 ϕ=π4\phi = \frac{\pi}{4}(因为 tan⁡ϕ=x2+y2z=1\tan\phi = \frac{\sqrt{x^2+y^2}}{z} = 1)

球坐标下: Ω′:0≤r≤R,0≤ϕ≤π4,0≤θ≤2π\Omega': 0 \leq r \leq R, \quad 0 \leq \phi \leq \frac{\pi}{4}, \quad 0 \leq \theta \leq 2\pi

4. 球心不在原点的球

对于球心在 (0,0,a)(0,0,a)、半径为 aa 的球 x2+y2+z2≤2azx^2+y^2+z^2 \leq 2az:

代入球坐标:r2≤2arcos⁡ϕ⇒r≤2acos⁡ϕr^2 \leq 2ar\cos\phi \Rightarrow r \leq 2a\cos\phi

由 r≥0r \geq 0 得 cos⁡ϕ≥0\cos\phi \geq 0,即 0≤ϕ≤π20 \leq \phi \leq \frac{\pi}{2}

Ω′:0≤r≤2acos⁡ϕ,0≤ϕ≤π2,0≤θ≤2π\Omega': 0 \leq r \leq 2a\cos\phi, \quad 0 \leq \phi \leq \frac{\pi}{2}, \quad 0 \leq \theta \leq 2\pi

适用场景总结

典型例题

例题1:球体上的积分

计算 ∭Ω(x2+y2+z2) dv\iiint_\Omega (x^2+y^2+z^2) \, dv,其中 Ω:x2+y2+z2≤R2\Omega: x^2+y^2+z^2 \leq R^2。

参考答案
三重积分球坐标球体球对称

解题思路: 球体 + 被积函数含 x2+y2+z2x^2+y^2+z^2,球坐标是最佳选择。

详细步骤:

  1. 球坐标变换,体积元素 r2sin⁡ϕ dr dϕ dθr^2\sin\phi \, dr \, d\phi \, d\theta
  2. 积分区域:0≤r≤R, 0≤ϕ≤π, 0≤θ≤2π0 \leq r \leq R, \, 0 \leq \phi \leq \pi, \, 0 \leq \theta \leq 2\pi
  3. 被积函数:x2+y2+z2=r2x^2+y^2+z^2 = r^2
  4. ∭Ω(x2+y2+z2)dv=∫02πdθ∫0πsin⁡ϕ dϕ∫0Rr2⋅r2 dr\iiint_\Omega (x^2+y^2+z^2) dv = \int_0^{2\pi} d\theta \int_0^\pi \sin\phi \, d\phi \int_0^R r^2 \cdot r^2 \, dr =2π⋅[−cos⁡ϕ]0π⋅[r55]0R=2π⋅2⋅R55=4πR55= 2\pi \cdot [-\cos\phi]_0^\pi \cdot \left[\frac{r^5}{5}\right]_0^R = 2\pi \cdot 2 \cdot \frac{R^5}{5} = \frac{4\pi R^5}{5}

答案:4πR55\frac{4\pi R^5}{5}

例题2:球锥区域

计算 ∭Ωz dv\iiint_\Omega z \, dv,其中 Ω\Omega 由 x2+y2+z2≤2Rzx^2+y^2+z^2 \leq 2Rz 和 z≥x2+y2z \geq \sqrt{x^2+y^2} 围成。

参考答案
三重积分球坐标球锥区域偏心球

解题思路: 球面 x2+y2+z2=2Rzx^2+y^2+z^2=2Rz 即 x2+y2+(z−R)2=R2x^2+y^2+(z-R)^2=R^2(球心在 (0,0,R)(0,0,R),半径 RR)。圆锥面 z=x2+y2z=\sqrt{x^2+y^2} 对应 ϕ=π4\phi=\frac{\pi}{4}。

详细步骤:

  1. 球面的球坐标方程:r2=2Rrcos⁡ϕ⇒r=2Rcos⁡ϕr^2 = 2Rr\cos\phi \Rightarrow r = 2R\cos\phi
  2. 圆锥面:ϕ=π4\phi = \frac{\pi}{4}
  3. 积分区域:0≤r≤2Rcos⁡ϕ, 0≤ϕ≤π4, 0≤θ≤2π0 \leq r \leq 2R\cos\phi, \, 0 \leq \phi \leq \frac{\pi}{4}, \, 0 \leq \theta \leq 2\pi
  4. 被积函数:z=rcos⁡ϕz = r\cos\phi
  5. ∭Ωz dv=∫02πdθ∫0π/4dϕ∫02Rcos⁡ϕrcos⁡ϕ⋅r2sin⁡ϕ dr\iiint_\Omega z \, dv = \int_0^{2\pi} d\theta \int_0^{\pi/4} d\phi \int_0^{2R\cos\phi} r\cos\phi \cdot r^2\sin\phi \, dr =2π∫0π/4cos⁡ϕsin⁡ϕ dϕ∫02Rcos⁡ϕr3dr=2π∫0π/4cos⁡ϕsin⁡ϕ⋅(2Rcos⁡ϕ)44dϕ= 2\pi \int_0^{\pi/4} \cos\phi\sin\phi \, d\phi \int_0^{2R\cos\phi} r^3 dr = 2\pi \int_0^{\pi/4} \cos\phi\sin\phi \cdot \frac{(2R\cos\phi)^4}{4} d\phi =2π⋅4R4∫0π/4cos⁡5ϕsin⁡ϕ dϕ=8πR4⋅[−cos⁡6ϕ6]0π/4= 2\pi \cdot 4R^4 \int_0^{\pi/4} \cos^5\phi\sin\phi \, d\phi = 8\pi R^4 \cdot \left[-\frac{\cos^6\phi}{6}\right]_0^{\pi/4} =8πR4⋅16(1−18)=4πR43⋅78=7πR46= 8\pi R^4 \cdot \frac{1}{6}\left(1 - \frac{1}{8}\right) = \frac{4\pi R^4}{3} \cdot \frac{7}{8} = \frac{7\pi R^4}{6}

答案:7πR46\frac{7\pi R^4}{6}

例题3:球壳区域

计算 ∭Ω1x2+y2+z2 dv\iiint_\Omega \frac{1}{\sqrt{x^2+y^2+z^2}} \, dv,其中 Ω:1≤x2+y2+z2≤4\Omega: 1 \leq x^2+y^2+z^2 \leq 4。

参考答案
三重积分球坐标球壳反常积分

解题思路: 球壳区域,球坐标下积分限为常数,计算简单。

详细步骤:

  1. 球坐标:1≤r≤2, 0≤ϕ≤π, 0≤θ≤2π1 \leq r \leq 2, \, 0 \leq \phi \leq \pi, \, 0 \leq \theta \leq 2\pi
  2. 被积函数:1x2+y2+z2=1r\frac{1}{\sqrt{x^2+y^2+z^2}} = \frac{1}{r}
  3. ∭Ω1r⋅r2sin⁡ϕ dr dϕ dθ=∫02πdθ∫0πsin⁡ϕ dϕ∫12r dr=2π⋅2⋅[r22]12=4π⋅32=6π\iiint_\Omega \frac{1}{r} \cdot r^2\sin\phi \, dr \, d\phi \, d\theta = \int_0^{2\pi} d\theta \int_0^\pi \sin\phi \, d\phi \int_1^2 r \, dr = 2\pi \cdot 2 \cdot \left[\frac{r^2}{2}\right]_1^2 = 4\pi \cdot \frac{3}{2} = 6\pi

答案:6π6\pi


练习题

练习1

计算 ∭Ωx2+y2+z2 dv\iiint_\Omega \sqrt{x^2+y^2+z^2} \, dv,其中 Ω:x2+y2+z2≤1\Omega: x^2+y^2+z^2 \leq 1。

参考答案
三重积分球坐标单位球

解题思路: 单位球 + 被积函数含 x2+y2+z2\sqrt{x^2+y^2+z^2},球坐标。

详细步骤:

  1. 球坐标:0≤r≤1, 0≤ϕ≤π, 0≤θ≤2π0 \leq r \leq 1, \, 0 \leq \phi \leq \pi, \, 0 \leq \theta \leq 2\pi
  2. ∭Ωr⋅r2sin⁡ϕ dr dϕ dθ=∫02πdθ∫0πsin⁡ϕ dϕ∫01r3dr=2π⋅2⋅14=π\iiint_\Omega r \cdot r^2\sin\phi \, dr \, d\phi \, d\theta = \int_0^{2\pi} d\theta \int_0^\pi \sin\phi \, d\phi \int_0^1 r^3 dr = 2\pi \cdot 2 \cdot \frac{1}{4} = \pi

答案:π\pi

练习2

计算 ∭Ωz2 dv\iiint_\Omega z^2 \, dv,其中 Ω:x2+y2+z2≤R2\Omega: x^2+y^2+z^2 \leq R^2。

参考答案
三重积分球坐标球体对称性

解题思路: 球体上的积分,由对称性 ∭x2=∭y2=∭z2\iiint x^2 = \iiint y^2 = \iiint z^2,可以先求 ∭(x2+y2+z2)\iiint (x^2+y^2+z^2) 再除以3。

详细步骤:

  1. 由轮换对称性:∭Ωx2dv=∭Ωy2dv=∭Ωz2dv\iiint_\Omega x^2 dv = \iiint_\Omega y^2 dv = \iiint_\Omega z^2 dv
  2. ∭Ω(x2+y2+z2)dv=3∭Ωz2dv\iiint_\Omega (x^2+y^2+z^2) dv = 3\iiint_\Omega z^2 dv
  3. 由例题1知 ∭Ω(x2+y2+z2)dv=4πR55\iiint_\Omega (x^2+y^2+z^2) dv = \frac{4\pi R^5}{5}
  4. 因此 ∭Ωz2dv=13⋅4πR55=4πR515\iiint_\Omega z^2 dv = \frac{1}{3} \cdot \frac{4\pi R^5}{5} = \frac{4\pi R^5}{15}

答案:4πR515\frac{4\pi R^5}{15}


总结

本文出现的符号

符号类型读音/说明在本文中的含义
rr径向距离radial distance点到原点的距离,r≥0r \geq 0
ϕ\phi极角/余纬角polar angle与 z 轴正方向的夹角,ϕ∈[0,π]\phi \in [0,\pi]
θ\theta方位角azimuthal angle在 xOy 平面上与 x 轴的夹角
r2sin⁡ϕ dr dϕ dθr^2 \sin\phi \, dr \, d\phi \, d\theta球坐标体积元素spherical volume element球坐标系下的体积微元

中英对照

中文术语英文术语音标说明
球坐标spherical coordinates/ˈsfɪrɪkəl koʊˈɔːrdɪneɪts/用径向距离和两个角度表示点的坐标系
径向距离radial distance/ˈreɪdiəl ˈdɪstəns/点到原点的距离
极角polar angle/ˈpoʊlər ˈæŋɡəl/与 z 轴的夹角(余纬角)
方位角azimuthal angle/ˌæzɪˈmjuːθəl ˈæŋɡəl/在 xOy 平面上的投影角
球体sphere/sfɪr/球面围成的立体
球壳spherical shell/ˈsfɪrɪkəl ʃɛl/两个同心球面之间的区域
球锥spherical cone/ˈsfɪrɪkəl koʊn/球面与圆锥面围成的区域
雅可比行列式Jacobian determinant/dʒəˈkoʊbiən dɪˈtɜːrmɪnənt/坐标变换的偏导数行列式

发现内容有问题?

告诉我们哪里需要改进。