球坐标变换与其他常用变换

球坐标变换

变换公式与雅可比

球坐标变换

{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) 是方位角。

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

体积元素:dv=dxdydz=Jdrdϕdθ=r2sinϕdrdϕdθdv = dx\,dy\,dz = |J| \, dr \, d\phi \, d\theta = r^2 \sin\phi \, dr \, d\phi \, d\theta

换元公式: Ω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

适用场景

广义极坐标变换(椭圆变换)

广义极坐标变换

{x=aρcosθy=bρsinθ\begin{cases} x = a\rho \cos\theta \\ y = b\rho \sin\theta \end{cases}

雅可比行列式:J=abρJ = ab\rho

面积元素:dσ=abρdρdθd\sigma = ab\rho \, d\rho \, d\theta

适用场景:积分区域是椭圆 x2a2+y2b21\frac{x^2}{a^2} + \frac{y^2}{b^2} \leq 1 或椭圆扇形。

椭圆区域在广义极坐标下变为:0ρ1,0θ2π0 \leq \rho \leq 1, \, 0 \leq \theta \leq 2\pi

推广:广义球坐标(椭球变换)

{x=arsinϕcosθy=brsinϕsinθz=crcosϕ\begin{cases} x = ar\sin\phi\cos\theta \\ y = br\sin\phi\sin\theta \\ z = cr\cos\phi \end{cases}

雅可比行列式:J=abcr2sinϕJ = abcr^2\sin\phi

适用场景:椭球区域 x2a2+y2b2+z2c21\frac{x^2}{a^2} + \frac{y^2}{b^2} + \frac{z^2}{c^2} \leq 1

双曲变换

双曲变换

{x=au2v22y=auv\begin{cases} x = a \cdot \frac{u^2 - v^2}{2} \\ y = a \cdot uv \end{cases}

雅可比行列式: J=auavavau=a2u2+a2v2=a2(u2+v2)J = \begin{vmatrix} au & -av \\ av & au \end{vmatrix} = a^2u^2 + a^2v^2 = a^2(u^2 + v^2)

面积元素:dσ=a2(u2+v2)dudvd\sigma = a^2(u^2+v^2) \, du \, dv

适用场景:积分区域由抛物线 y2=2a(x+a2)y^2 = 2a(x+\frac{a}{2})y2=2a(xa2)y^2 = -2a(x-\frac{a}{2}) 围成,或被积函数含有 x+x2+y2\sqrt{x+\sqrt{x^2+y^2}} 等形式。

指数变换(对数极坐标)

指数变换

{x=aeucosvy=aeusinv\begin{cases} x = ae^u \cos v \\ y = ae^u \sin v \end{cases}

雅可比行列式:J=a2e2uJ = a^2 e^{2u}

面积元素:dσ=a2e2ududvd\sigma = a^2 e^{2u} \, du \, dv

适用场景:被积函数含有 1(x2+y2)k\frac{1}{(x^2+y^2)^k} 且积分区域包含原点(此时极坐标下 ρ=0\rho=0 处可能发散,指数变换将原点映射到 u=u=-\infty,便于处理反常积分)。

线性变换

线性变换

{x=a11u+a12vy=a21u+a22v\begin{cases} x = a_{11}u + a_{12}v \\ y = a_{21}u + a_{22}v \end{cases}

雅可比行列式:J=a11a12a21a22=a11a22a12a21J = \begin{vmatrix} a_{11} & a_{12} \\ a_{21} & a_{22} \end{vmatrix} = a_{11}a_{22} - a_{12}a_{21}(常数)

面积元素:dσ=Jdudvd\sigma = |J| \, du \, dv

适用场景:积分区域是平行四边形,或被积函数是线性组合的函数。线性变换将平行四边形变为矩形,极大简化积分限。

常见线性变换

变换公式雅可比适用场景
平移x=u+a,y=v+bx = u + a, \, y = v + bJ=1J = 1区域中心不在原点
缩放x=au,y=bvx = au, \, y = bvJ=abJ = ab椭圆/矩形区域
旋转x=ucosαvsinα,y=usinα+vcosαx = u\cos\alpha - v\sin\alpha, \, y = u\sin\alpha + v\cos\alphaJ=1J = 1区域有倾斜对称轴
剪切x=u+cv,y=vx = u + cv, \, y = vJ=1J = 1平行四边形区域

典型例题

例题1:椭圆区域面积

用广义极坐标求椭圆 x2a2+y2b21\frac{x^2}{a^2} + \frac{y^2}{b^2} \leq 1 的面积。

参考答案(3 个标签)
广义极坐标椭圆面积变量替换

解题思路: 广义极坐标变换将椭圆区域变为单位圆区域。

详细步骤

  1. 变换:x=aρcosθ,y=bρsinθx = a\rho\cos\theta, \, y = b\rho\sin\theta,雅可比 J=abρJ = ab\rho
  2. 椭圆区域变为:0ρ1,0θ2π0 \leq \rho \leq 1, \, 0 \leq \theta \leq 2\pi
  3. 面积 S=Ddxdy=02πdθ01abρdρ=2πab12=πabS = \iint_D dx\,dy = \int_0^{2\pi} d\theta \int_0^1 ab\rho \, d\rho = 2\pi \cdot ab \cdot \frac{1}{2} = \pi ab

答案πab\pi ab

例题2:线性变换计算平行四边形区域积分

计算 D(x+y)dxdy\iint_D (x+y) \, dx\,dy,其中 DD 是由 x+y=0,x+y=1,xy=0,xy=1x+y=0, x+y=1, x-y=0, x-y=1 围成的平行四边形。

参考答案(3 个标签)
线性变换平行四边形变量替换

解题思路: 令 u=x+y,v=xyu = x+y, \, v = x-y,则区域变为矩形 0u1,0v10 \leq u \leq 1, \, 0 \leq v \leq 1

详细步骤

  1. 逆变换:x=u+v2,y=uv2x = \frac{u+v}{2}, \, y = \frac{u-v}{2}
  2. 雅可比:J=(x,y)(u,v)=1/21/21/21/2=12J = \frac{\partial(x,y)}{\partial(u,v)} = \begin{vmatrix} 1/2 & 1/2 \\ 1/2 & -1/2 \end{vmatrix} = -\frac{1}{2}J=12|J| = \frac{1}{2}
  3. 被积函数:x+y=ux+y = u
  4. D(x+y)dxdy=0101u12dudv=1201dv01udu=12112=14\iint_D (x+y) dx\,dy = \int_0^1 \int_0^1 u \cdot \frac{1}{2} du\,dv = \frac{1}{2} \int_0^1 dv \int_0^1 u \, du = \frac{1}{2} \cdot 1 \cdot \frac{1}{2} = \frac{1}{4}

答案14\frac{1}{4}

例题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

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

解题思路: 球壳区域,球坐标变换。

详细步骤

  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

用广义极坐标计算 Dx2a2+y2b2dxdy\iint_D \sqrt{\frac{x^2}{a^2} + \frac{y^2}{b^2}} \, dx\,dy,其中 D:x2a2+y2b21D: \frac{x^2}{a^2} + \frac{y^2}{b^2} \leq 1

参考答案(3 个标签)
广义极坐标椭圆区域变量替换

解题思路: 广义极坐标变换后,被积函数变为 ρ\rho

详细步骤

  1. 变换:x=aρcosθ,y=bρsinθx = a\rho\cos\theta, \, y = b\rho\sin\thetaJ=abρJ = ab\rho
  2. 区域:0ρ1,0θ2π0 \leq \rho \leq 1, \, 0 \leq \theta \leq 2\pi
  3. 被积函数:x2a2+y2b2=ρ\sqrt{\frac{x^2}{a^2}+\frac{y^2}{b^2}} = \rho
  4. Dρabρdρdθ=ab02πdθ01ρ2dρ=ab2π13=2πab3\iint_D \rho \cdot ab\rho \, d\rho \, d\theta = ab \int_0^{2\pi} d\theta \int_0^1 \rho^2 d\rho = ab \cdot 2\pi \cdot \frac{1}{3} = \frac{2\pi ab}{3}

答案2πab3\frac{2\pi ab}{3}

练习2

用线性变换计算 Dexyx+ydxdy\iint_D e^{\frac{x-y}{x+y}} \, dx\,dy,其中 DDx+y=1,x=0,y=0x+y=1, x=0, y=0 围成。

参考答案(3 个标签)
线性变换三角形区域变量替换

解题思路: 令 u=x+y,v=xyu = x+y, \, v = x-y,简化被积函数。

详细步骤

  1. 逆变换:x=u+v2,y=uv2x = \frac{u+v}{2}, \, y = \frac{u-v}{2}J=12|J| = \frac{1}{2}
  2. 区域变换:x0,y0u+v0,uv0uvux \geq 0, y \geq 0 \Rightarrow u+v \geq 0, u-v \geq 0 \Rightarrow -u \leq v \leq ux+y1u1x+y \leq 1 \Rightarrow u \leq 1
  3. 区域:0u1,uvu0 \leq u \leq 1, \, -u \leq v \leq u
  4. 被积函数:evue^{\frac{v}{u}}
  5. Dev/u12dudv=1201duuuev/udv=1201u[ev/u]uudu=1201u(ee1)du=ee1212=ee14\iint_D e^{v/u} \cdot \frac{1}{2} du\,dv = \frac{1}{2}\int_0^1 du \int_{-u}^u e^{v/u} dv = \frac{1}{2}\int_0^1 u[e^{v/u}]_{-u}^u du = \frac{1}{2}\int_0^1 u(e - e^{-1}) du = \frac{e-e^{-1}}{2} \cdot \frac{1}{2} = \frac{e-e^{-1}}{4}

答案ee14\frac{e-e^{-1}}{4}


总结

本文出现的符号

符号类型读音/说明在本文中的含义
rr径向距离radial distance球坐标中点到原点的距离
ϕ\phi极角polar angle球坐标中与 z 轴的夹角
θ\theta方位角azimuthal angle球坐标中在 xOy 平面上的投影角
a,ba, b半轴长semi-axes椭圆/椭球的半轴长度
u,v,wu, v, w新坐标new coordinates变量替换后的新自变量

中英对照

中文术语英文术语音标说明
球坐标变换spherical coordinate transformation/ˈsfɪrɪkəl koʊˈɔːrdɪneɪt ˌtrænsfərˈmeɪʃən/直角坐标到球坐标的变换
广义极坐标generalized polar coordinates/ˈdʒɛnrəlaɪzd ˈpoʊlər koʊˈɔːrdɪneɪts/椭圆变换,极坐标的推广
双曲变换hyperbolic transformation/ˌhaɪpərˈbɒlɪk ˌtrænsfərˈmeɪʃən/基于抛物线坐标的变换
线性变换linear transformation/ˈlɪniər ˌtrænsfərˈmeɪʃən/变换公式为线性的变量替换
椭球变换ellipsoidal transformation/ɪˈlɪpsɔɪdl ˌtrænsfərˈmeɪʃən/广义球坐标,椭球区域变换
指数变换exponential transformation/ˌɛkspəˈnɛnʃəl ˌtrænsfərˈmeɪʃən/对数极坐标变换
平移变换translation/trænsˈleɪʃən/区域平移
旋转变换rotation/roʊˈteɪʃən/区域旋转