球坐标下的三重积分计算

为什么需要球坐标

当积分区域是球体、球锥、半球等球形区域,或被积函数含有 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}

其中:

  • r0r \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=dxdydz=r2sinϕdrdϕdθdv = dx\,dy\,dz = r^2 \sin\phi \, dr \, d\phi \, d\theta

角度范围的几何意义

  • ϕ=0\phi = 0zz 轴正方向
  • ϕ=π2\phi = \frac{\pi}{2}xOyxOy 平面
  • ϕ=π\phi = \pizz 轴负方向
  • θ=0\theta = 0xOzxOz 平面的右半部分(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ϕdrdϕ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+z2R2\Omega: x^2 + y^2 + z^2 \leq R^2

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

Ωfdv=02πdθ0πsinϕdϕ0Rf(rsinϕcosθ,)r2dr\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+z2R2,z0\Omega: x^2 + y^2 + z^2 \leq R^2, \quad z \geq 0

球坐标下: Ω:0rR,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

球坐标下: Ω:0rR,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+z22azx^2+y^2+z^2 \leq 2az

代入球坐标:r22arcosϕr2acosϕr^2 \leq 2ar\cos\phi \Rightarrow r \leq 2a\cos\phi

r0r \geq 0cosϕ0\cos\phi \geq 0,即 0ϕπ20 \leq \phi \leq \frac{\pi}{2}

Ω:0r2acosϕ,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+z2R2\Omega: x^2+y^2+z^2 \leq R^2

参考答案(4 个标签)
三重积分球坐标球体球对称

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

详细步骤

  1. 球坐标变换,体积元素 r2sinϕdrdϕdθr^2\sin\phi \, dr \, d\phi \, d\theta
  2. 积分区域:0rR,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ϕ0Rr2r2dr\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π2R55=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:球锥区域

计算 Ωzdv\iiint_\Omega z \, dv,其中 Ω\Omegax2+y2+z22Rzx^2+y^2+z^2 \leq 2Rzzx2+y2z \geq \sqrt{x^2+y^2} 围成。

参考答案(4 个标签)
三重积分球坐标球锥区域偏心球

解题思路: 球面 x2+y2+z2=2Rzx^2+y^2+z^2=2Rzx2+y2+(zR)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. 积分区域:0r2Rcosϕ,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. Ωzdv=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π4R40π/4cos5ϕsinϕdϕ=8πR4[cos6ϕ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πR416(118)=4πR4378=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+z2dv\iiint_\Omega \frac{1}{\sqrt{x^2+y^2+z^2}} \, dv,其中 Ω:1x2+y2+z24\Omega: 1 \leq x^2+y^2+z^2 \leq 4

参考答案(4 个标签)
三重积分球坐标球壳反常积分

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

详细步骤

  1. 球坐标:1r2,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. Ω1rr2sinϕdrdϕdθ=02πdθ0πsinϕdϕ12rdr=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+z2dv\iiint_\Omega \sqrt{x^2+y^2+z^2} \, dv,其中 Ω:x2+y2+z21\Omega: x^2+y^2+z^2 \leq 1

参考答案(3 个标签)
三重积分球坐标单位球

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

详细步骤

  1. 球坐标:0r1,0ϕπ,0θ2π0 \leq r \leq 1, \, 0 \leq \phi \leq \pi, \, 0 \leq \theta \leq 2\pi
  2. Ωrr2sinϕdrdϕdθ=02πdθ0πsinϕdϕ01r3dr=2π214=π\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

计算 Ωz2dv\iiint_\Omega z^2 \, dv,其中 Ω:x2+y2+z2R2\Omega: x^2+y^2+z^2 \leq R^2

参考答案(4 个标签)
三重积分球坐标球体对称性

解题思路: 球体上的积分,由对称性 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=134π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点到原点的距离,r0r \geq 0
ϕ\phi极角/余纬角polar angle与 z 轴正方向的夹角,ϕ[0,π]\phi \in [0,\pi]
θ\theta方位角azimuthal angle在 xOy 平面上与 x 轴的夹角
r2sinϕdrdϕ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/坐标变换的偏导数行列式