1926 字
10 分钟
9.3 三重积分

§\S9.3 三重积分#

  • abf(x)dx\displaystyle \int_a^b f(x)\mathrm dx:曲边梯形面积,铜棒质量
  • Df(x,y)dxdy\displaystyle \iint_D f(x,y)\mathrm dx \mathrm dy:曲顶柱体体积,薄片质量
  • Ωf(x,y,z)dxdydz\displaystyle \iiint_\Omega f(x,y,z) \mathrm dx \mathrm dy \mathrm dz:空间物体的质量

一、定义#

f(x,y,z)f(x,y,z) 为空间有界闭区域 Ω\Omega 的有界函数,任意分割 ΔV1ΔVn\Delta V_1 \cdots \Delta V_n,任意取点 (ai,bi,ci)ΔVi(a_i,b_i,c_i)\in \Delta V_i,称 limλ0i=1nf(ai,bi,ci)ΔVi\displaystyle \lim_{\lambda \to 0}\sum_{i=1}^nf(a_i,b_i,c_i)\cdot \Delta V_if(x,y,z)f(x,y,z)Ω\Omega 上的三重积分,记为

Ωf(x,y,z)dv=Ωf(x,y,z)dxdydz\begin{align*} \iiint_\Omega f(x,y,z)\mathrm dv=\iiint_\Omega f(x,y,z) \mathrm dx \mathrm dy \mathrm dz \end{align*}

f(x,y,z)=1f(x,y,z)=1 时,Ω1dV=ΩdV=V\displaystyle \iiint_\Omega 1\cdot \mathrm dV=\iiint_\Omega \mathrm dV = V

二、性质#

线性性质#

Ω(lf+kg)dV=lΩfdV+kΩgdV\begin{align*} \iiint_\Omega (lf+kg) \mathrm dV = l\iiint_\Omega f \mathrm dV + k\iiint_\Omega g \mathrm dV \end{align*}

可加性#

ΩfdV=Ω1fdV+Ω2fdV(Ω=Ω1+Ω2)\begin{align*} \iiint_\Omega f \mathrm dV &= \iiint_{\Omega_1} f \mathrm dV + \iiint_{\Omega_2} f\mathrm dV\\ (\Omega &= \Omega_1 + \Omega_2) \end{align*}

保序性#

f0ΩfdV0\begin{align*} f\ge 0 \Rightarrow \iiint_\Omega f \mathrm dV\ge 0 \end{align*}

估值定理#

ffΩ\Omega 连续,mfMm\le f\le M,则 mVΩfdVMV\displaystyle mV\le \iiint_\Omega f\cdot \mathrm dV\le MV

中值定理#

Ωf(x,y,z)dV=f(a,b,c)V\begin{align*} \iiint_\Omega f(x,y,z) \mathrm dV = f(a,b,c) V \end{align*}

对称性#

  • Ω\Omega 关于 xOyxOy 面对称,Ω\Omega_上Ω\Omega_下

    ΩfdV={2ΩfdVf(x,y,z)=f(x,y,z)0f(x,y,z)=f(x,y,z)\begin{align*} \iiint_\Omega f \mathrm dV= \begin{cases} 2\iiint_{\Omega_上} f \mathrm dV &\quad f(x,y,z)=f(x,y,-z)\\ 0 &\quad f(x,y,z)=-f(x,y,-z) \end{cases} \end{align*}
  • Ω\Omega 关于 yOzyOz 面对称,Ω\Omega_前Ω\Omega_后

    ΩfdV={2ΩfdVf(x,y,z)=f(x,y,z)0f(x,y,z)=f(x,y,z)\begin{align*} \iiint_\Omega f \mathrm dV= \begin{cases} 2\iiint_{\Omega_前} f \mathrm dV &\quad f(x,y,z)=f(-x,y,z)\\ 0 &\quad f(x,y,z)=-f(-x,y,z) \end{cases} \end{align*}
  • Ω\Omega 关于 xOzxOz 面对称,Ω\Omega_左Ω\Omega_右

    ΩfdV={2ΩfdVf(x,y,z)=f(x,y,z)0f(x,y,z)=f(x,y,z)\begin{align*} \iiint_\Omega f \mathrm dV= \begin{cases} 2\iiint_{\Omega_右} f \mathrm dV &\quad f(x,y,z)=f(x,-y,z)\\ 0 &\quad f(x,y,z)=-f(x,-y,z) \end{cases} \end{align*}
  • 轮换对称性:假设 x,y,zx,y,z 交换后 Ω\Omega 不变,则 Ωf(x,y,z)dV=Ωf(y,z,x)dV=Ωf(z,x,y)dV\displaystyle \iiint_\Omega f(x,y,z)\mathrm dV=\iiint_\Omega f(y,z,x)\mathrm dV=\iiint_\Omega f(z,x,y)\mathrm dV

例如:Ω:x2+y2+z2R2\Omega: x^2+y^2+z^2\le R^2,有 Ωx2dV=Ωy2dV=Ωz2dV=13Ω(x2+y2+z2)dV\displaystyle \iiint_\Omega x^2 \mathrm dV=\iiint_\Omega y^2 \mathrm dV=\iiint_\Omega z^2 \mathrm dV=\dfrac{1}{3}\iiint_\Omega (x^2+y^2+z^2) \mathrm dV

三、三重积分的计算方法#

直角坐标系#

  1. 投影法(“先一后二”)

    • 投影到 xOyxOy 面上:I=Dxydxdyϕ1(x,y)ϕ2(x,y)f(x,y,z)dz\displaystyle I=\iint_{D_{xy}} \mathrm dx \mathrm dy \int_{\phi_1(x,y)}^{\phi_2(x,y)} f(x,y,z) \mathrm dzDxyD_{xy}Ω\OmegaxOyxOy 上投影区域)
    • 投影到 yOzyOz 面上:I=Dyzdydzϕ1(y,z)ϕ2(y,z)f(x,y,z)dx\displaystyle I=\iint_{D_{yz}} \mathrm dy \mathrm dz \int_{\phi_1(y,z)}^{\phi_2(y,z)} f(x,y,z) \mathrm dxDyzD_{yz}Ω\OmegayOzyOz 上投影区域)
    • 投影到 xOzxOz 面上:I=Dxzdxdzϕ1(x,z)ϕ2(x,z)f(x,y,z)dy\displaystyle I=\iint_{D_{xz}} \mathrm dx \mathrm dz \int_{\phi_1(x,z)}^{\phi_2(x,z)} f(x,y,z) \mathrm dyDxzD_{xz}Ω\OmegaxOzxOz 上投影区域)

    投影法中的穿线法:垂直投影面,方向与坐标轴正向同向,先穿过积分下限面,再穿过积分上限面

  2. 截面法(“先二后一”)适用于被积函数为一元函数,且截面为三角形、矩形、圆的情况

    I=Ωg(x)dV=abdxDg(x)dydz=ab[g(x)Ddydz]dx=abg(x)δ(x)x=常数与Ω截面面积dx\begin{align*} I&=\iiint_\Omega g(x) \mathrm dV=\int_a^b \mathrm dx \iint_D g(x) \mathrm dy \mathrm dz \\ &=\int_a^b [g(x) \iint_D \mathrm dy \mathrm dz] \mathrm dx \\ &=\int_a^b g(x) \cdot \underset{x=常数与\Omega截面面积}{\delta(x)} \mathrm dx \end{align*}

例题#

  1. Ω:x2+y2+z2R2\Omega: x^2+y^2+z^2\le R^2

    I=Ωf(x,y,z)dxdydz=DxydxdyR2x2y2R2x2y2f(x,y,z)dz\begin{align*} I&=\iiint_\Omega f(x,y,z) \mathrm dx \mathrm dy \mathrm dz \\ &=\iint_{D_{xy}}\mathrm dx \mathrm dy \int_{-\sqrt{R^2-x^2-y^2}}^{\sqrt{R^2-x^2-y^2}}f(x,y,z) \mathrm dz \end{align*}

    Dxy:x2+y2R2D_{xy}:x^2+y^2\le R^2

  2. I=ΩfdVI=\iiint_\Omega f\cdot \mathrm dV,其中 Ω:x2+y2+(z2)24\Omega: x^2+y^2+(z-2)^2\le 4z=1z=1 所围成的包含球心的部分

    Ω\Omega 投影到 xOyxOy 平面上,Dxy:x2+y24D_{xy}:x^2+y^2\le 4

    {x2+y2+(z2)2=4z=2z{x2+y2=4z=0\begin{align*} \begin{cases} x^2+y^2+(z-2)^2=4\\ z=2 \end{cases} \overset{消z}{\Rightarrow } \begin{cases} x^2+y^2=4\\ z=0 \end{cases} \end{align*} {x2+y2+(z2)2=4z=1z{x2+y2=3z=0\begin{align*} \begin{cases} x^2+y^2+(z-2)^2=4\\ z=1 \end{cases} \overset{消z}{\Rightarrow } \begin{cases} x^2+y^2=3\\ z=0 \end{cases} \end{align*}

    I=D1dxdy12+4x2y2fdz+D2dxdy24x2y22+4x2y2fdz\displaystyle I=\iint_{D_1}\mathrm dx \mathrm dy\int_1^{2+\sqrt{4-x^2-y^2}} f \mathrm dz+\iint_{D_2}\mathrm dx \mathrm dy \int_{2-\sqrt{4-x^2-y^2}}^{2+\sqrt{4-x^2-y^2}} f \mathrm dz,其中 D1:z=1z=2+4x2y2D_1:z=1\to z=2+\sqrt{4-x^2-y^2}D2:z=24x2y2z=2+4x2y2D_2:z=2-\sqrt{4-x^2-y^2}\to z=2+\sqrt{4-x^2-y^2}

  3. Ωxdxdydz\iiint_\Omega x \mathrm dx \mathrm dy \mathrm dz,其中 Ω\Omegax+y+z=1x+y+z=1 与坐标面所围成的区域

    • 投影法

      • Ω\Omega 投影到 xOyxOy 面上,用平行 z 轴的射线穿过积分区域。Dxy:x+y1(x0,y0)D_{xy}:x+y\le 1(x\ge 0,y\ge 0)

        I=Dxydxdy01xyxdz=Dxyx(1xy)dxdy=01dx01xx(1xy)dy=124\begin{align*} I&=\iint_{D_{xy}} \mathrm dx \mathrm dy \int_0^{1-x-y}x \mathrm dz\\ &=\iint_{D_{xy}} x(1-x-y) \mathrm dx \mathrm dy\\ &=\int_0^1 \mathrm dx \int_0^{1-x} x(1-x-y) \mathrm dy\\ &=\dfrac{1}{24} \end{align*}
      • Ω\Omega 投影到 yOzyOz 面上,用平行 x 轴的射线穿过积分区域。Dyz:y+z1(y0,z0)D_{yz}:y+z\le 1(y\ge 0,z\ge 0)

        I=Dyzdydz01yzxdz=Dyzx(1yz)dydz=01dy01yx(1yz)dz=124\begin{align*} I&=\iint_{D_{yz}} \mathrm dy \mathrm dz \int_0^{1-y-z} x \mathrm dz\\ &=\iint_{D_{yz}} x(1-y-z) \mathrm dy \mathrm dz\\ &=\int_0^1 \mathrm dy \int_0^{1-y} x(1-y-z) \mathrm dz\\ &=\dfrac{1}{24} \end{align*}
    • 截面法

      • 用平行 xOyxOy 的面截 Ω\Omega

        I=01dzDzxdxdy=01dz01zdx01xzxdy=124\begin{align*} I&=\int_0^1 \mathrm dz \iint_{D_z} x \mathrm dx \mathrm dy\\ &=\int_0^1 \mathrm dz \int_0^{1-z} \mathrm dx \int_0^{1-x-z} x \mathrm dy\\ &=\dfrac{1}{24} \end{align*}
      • 用平行 yOzyOz 的面截 Ω\Omega

        I=01dxDxxdydz=01(xDxdydz)=01xδ(x)dx=1201x(1x)2dx=124\begin{align*} I&=\int_0^1 \mathrm dx \iint_{D_x} x \mathrm dy \mathrm dz\\ &=\int_0^1 (x\iint_{D_x} \mathrm dy \mathrm dz)\\ &=\int_0^1 x \delta (x)\mathrm dx\\ &=\dfrac{1}{2}\int_0^1 x(1-x)^2 \mathrm dx\\ &=\dfrac{1}{24} \end{align*}
  4. I=Ωz2dVI=\iiint_\Omega z^2 \mathrm dVΩ\Omegaz=x2+y2z=\sqrt{x^2+y^2}z=1z=1 所围成的闭区域

    截面法较易,用平行 xOyxOy 的面截 Ω\Omega

    I=01dzDzz2dxdy=01z2δ(z)dz=01z2πz2dz=15π\begin{align*} I&=\int_0^1 \mathrm dz \iint_{D_z}z^2 \mathrm dx \mathrm dy\\ &=\int_0^1 z^2 \delta (z) \mathrm dz \\ &=\int_0^1 z^2\pi z^2 \mathrm dz \\ &=\dfrac{1}{5}\pi \end{align*}
  5. I=Ωz2dVI=\iiint_\Omega z^2 \mathrm dVΩ\Omegaz=x2+y2z=x^2+y^2z=1z=1 所围成的闭区域

    I=01dzDzz2dxdy=01z2δ2(z)dz=01z2πzdz=π4\begin{align*} I&=\int_0^1 \mathrm dz \iint_{D_z}z^2 \mathrm dx \mathrm dy\\ &=\int_0^1 z^2\delta^2(z) \mathrm dz\\ &=\int_0^1 z^2 \pi z \mathrm dz\\ &=\dfrac{\pi}{4} \end{align*}
  6. I=Ωx2dVI=\iiint_\Omega x^2 \mathrm dVΩ\Omegaz=x2+y2z=\sqrt{x^2+y^2}z=1z=1 所围成的闭区域

    • 投影法:将 Ω\Omega 投到 xOyxOy 平面上

      I=Dxydxdyx2+y21x2dz(Dxy:x2+y21)=Dxyx2(1x2+y2)dxdy=02πdθρ2cos2θ(1ρ)ρdρ=π5\begin{align*} I&=\iint_{D_{xy}}\mathrm dx \mathrm dy \int_{\sqrt{x^2+y^2}}^1 x^2 \mathrm dz(D_{xy}:x^2+y^2\le 1)\\ &=\iint_{D_{xy}}x^2(1-\sqrt{x^2+y^2})\mathrm dx \mathrm dy\\ &=\int_0^{2\pi} \mathrm d\theta \rho^2 \cos^2\theta(1-\rho)\cdot \rho \mathrm d\rho\\ &=\dfrac{\pi}{5} \end{align*}
    • 截面法:用平行于 xOyxOy 的平面去截

      I=01dzDzx2dxdy(Dz:x2+y2z2,关于y=x对称)=01dz02πdθ0zρ2cos2θρdρ=π5\begin{align*} I&=\int_0^1 \mathrm dz \iint_{D_z}x^2 \mathrm dx \mathrm dy(D_z:x^2+y^2\le z^2,关于 y=x 对称)\\ &=\int_0^1 \mathrm dz \int_0^{2\pi} \mathrm d\theta \int_0^z \rho^2 \cos^2\theta \rho \mathrm d\rho\\ &=\dfrac{\pi}{5} \end{align*}
  7. I=Ω(x+z)dVI=\iiint_\Omega (x+z) \mathrm dVΩ\Omegaz=x2+y2z=\sqrt{x^2+y^2}z=1x2y2z=\sqrt{1-x^2-y^2} 所围成的闭区域

    投影到 xOyxOy 平面上:Dxy:x2+y212D_{xy}: x^2+y^2\le \dfrac{1}{2}

    I=Dxydxdyx2+y21x2y2(x+z)dz=Dxyx(1x2y2x2+y2)+Dxy12(1x2y2x2y2)dxdy=π8\begin{align*} I&=\iint_{D_{xy}}\mathrm dx \mathrm dy \int_{\sqrt{x^2+y^2}}^{\sqrt{1-x^2-y^2}}(x+z)\mathrm dz\\ &=\iint_{D_{xy}} x(\sqrt{1-x^2-y^2}-\sqrt{x^2+y^2})+\iint_{D_{xy}}\dfrac{1}{2}(1-x^2-y^2-x^2-y^2)\mathrm dx \mathrm dy\\ &=\dfrac{\pi}{8} \end{align*}
  8. 0adx0xdy0yf(z)dz\int_0^a \mathrm dx \int_0^x \mathrm dy \int_0^y f(z) \mathrm dz 化为定积分(交换积分次序为 xyzx\to y\to zyxzy\to x\to z

    • 法一:先确定积分区域 Ω\Omega(0xa,0yx,0zy)(0\le x\le a, 0\le y\le x, 0\le z\le y)

      • 截面法:z=z0z=z_0

        I=0adzDxyf(z)dxdy=0af(z)12(az)2dz=120af(z)(az)2dz\begin{align*} I&=\int_0^a \mathrm dz \iint_{D_{xy}} f(z)\mathrm dx \mathrm dy\\ &=\int_0^a f(z) \dfrac{1}{2}(a-z)^2 \mathrm dz\\ &=\dfrac{1}{2}\int_0^a f(z)(a-z)^2 \mathrm dz \end{align*}
      • 投影法:ΩyOz\Omega\to yOz 平面

        I=Dyzdydzyaf(z)dx=0adzzadyyaf(z)dx=120af(z)(az)2dz\begin{align*} I&=\iint_{D_{yz}}\mathrm dy \mathrm dz \int_y^a f(z)\mathrm dx\\ &=\int_0^a \mathrm dz \int_z^a \mathrm dy \int_y^a f(z)\mathrm dx\\ &=\dfrac{1}{2}\int_0^a f(z)(a-z)^2 \mathrm dz \end{align*}
    • 法二:直接交换

      I=0adx0xdy0yf(z)dz=0adxDzf(z)dydz=0adx0xdzzxf(z)dy=Dxzdxdzzxf(z)dy=0adzzadxzxdy=120af(z)(az)2dz\begin{align*} I&=\int_0^a \mathrm dx \int_0^x \mathrm dy \int_0^y f(z) \mathrm dz\\ &=\int_0^a \mathrm dx \iint_{D_z}f(z) \mathrm dy \mathrm dz\\ &=\int_0^a \mathrm dx\int_0^x \mathrm dz \int_z^x f(z)\mathrm dy\\ &=\iint_{D_{xz}}\mathrm dx \mathrm dz\int_z^x f(z)\mathrm dy\\ &=\int_0^a \mathrm dz\int_z^a \mathrm dx \int_z^x \mathrm dy\\ &=\dfrac{1}{2}\int_0^a f(z)(a-z)^2 \mathrm dz \end{align*}
  9. I=01dx0xdy0xyf(x,y,z)dzI=\int_0^1 \mathrm dx\int_0^x \mathrm dy\int_0^{xy}f(x,y,z)\mathrm dz,交换积分次序为 yxzy\to x\to z

    原次序:zyxz\to y\to x

    原式=01dx0x2dzzxxf(x,y,z)dy[交换yz]=01dzz1dxzxxf(x,y,z)dy[交换xz]\begin{align*} 原式 &=\int_0^1 \mathrm dx \int_0^{x^2}\mathrm dz \int_{\frac{z}{x}}^xf(x,y,z) \mathrm dy[交换y、z]\\ &=\int_0^1 \mathrm dz \int_{\sqrt{z}}^1 \mathrm dx \int_{\frac{z}{x}}^xf(x,y,z) \mathrm dy[交换x、z] \end{align*}

柱坐标系(~投影法+极坐标)#

适用于 x2+y2x^2+y^2yx\dfrac{y}{x}xy\dfrac{x}{y}Ω\OmegaxOyxOy 面上,且 DxyD_{xy} 与圆有关

{x=ρcosθy=ρsinθz=z,dV=ρdρdθdz\begin{align*} \begin{cases} x=\rho\cos\theta\\ y=\rho\sin\theta\\ z=z \end{cases} ,dV=\rho \mathrm d\rho \mathrm d\theta \mathrm dz \end{align*}

I=αβdθϕ1(θ)ϕ2(θ)ρdρψ1(ρ,θ)ψ2(ρ,θ)f(ρcosθ,ρsinθ,z)dz\displaystyle I=\int_\alpha^\beta \mathrm d\theta \int_{\phi_1(\theta)}^{\phi_2(\theta)}\rho \mathrm d\rho \int_{\psi_1(\rho,\theta)}^{\psi_2(\rho,\theta)}f(\rho\cos\theta,\rho\sin\theta,z)\mathrm dz

例题#

  1. z=1a(x2+y2)z=\dfrac{1}{a}(x^2+y^2)z=2ax2+y2z=2a-\sqrt{x^2+y^2} 所围成的体积

    • 二重积分法

      V=V上曲V下曲=Dxy(2ax2+y2)dxdyDxy1a(x2+y2)dxdy=Dxy[2ax2+y21a(x2+y2)]dxdy(Dxy:x2+y2a2)=D1dxdy0adz+D2dxdya2adz=πa0azdz+πa2a(2az)2dz=5π6a3\begin{align*} V &= V_{上曲}-V_{下曲}\\ &=\iint_{D_{xy}}(2a-\sqrt{x^2+y^2})\mathrm dx \mathrm dy - \iint_{D_{xy}}\dfrac{1}{a}(x^2+y^2) \mathrm dx \mathrm dy\\ &=\iint_{D_{xy}}[2a-\sqrt{x^2+y^2}-\dfrac{1}{a}(x^2+y^2)] \mathrm dx \mathrm dy(D_{xy}:x^2+y^2\le a^2)\\ &=\iint_{D_1} \mathrm dx \mathrm dy \int_0^a \mathrm dz + \iint_{D_2} \mathrm dx \mathrm dy \int_a^{2a}\mathrm dz\\ &=\pi a \int_0^a z \mathrm dz + \pi \int_a^{2a}(2a-z)^2 \mathrm dz \\ &=\dfrac{5\pi}{6}a^3 \end{align*}
    • 三重积分法:V=Ω1dVV=\iiint_\Omega 1\cdot \mathrm dVΩ:z=x2+y2az=2ax2+y2\Omega:z=\dfrac{x^2+y^2}{a} \rightarrow z=2a-\sqrt{x^2+y^2}

      • 投影:ΩxOy\Omega \to xOyDxy:x2+y2a2D_{xy}:x^2+y^2\le a^2

        I=Dxydxdyx2+y2a2ax2+y21dz=Dxy[2ax2+y2x2+y2a]dxdy=02πdθ0a(2aρρ2a)ρdρ=5π6a3\begin{align*} I&=\iint_{D_{xy}}\mathrm dx \mathrm dy \int_{\frac{x^2+y^2}{a}}^{2a-\sqrt{x^2+y^2}} 1\cdot \mathrm dz\\ &=\iint_{D_{xy}}[2a-\sqrt{x^2+y^2}-\dfrac{x^2+y^2}{a}]\mathrm dx \mathrm dy\\ &=\int_0^{2\pi}\mathrm d\theta \int_0^a(2a-\rho-\dfrac{\rho^2}{a})\rho \mathrm d\rho\\ &=\dfrac{5\pi}{6}a^3 \end{align*}
      • 截面:用平行于 xOyxOy 的平面截 Ω\Omega

        I=0adzDz11dxdy+a2adzDz21dxdy=0aπazdz+a2aπ(2az)2dz=56πa3\begin{align*} I&=\int_0^a \mathrm dz \iint_{D_{z_1}}1 \mathrm dx \mathrm dy+\int_a^{2a} \mathrm dz \iint_{D_{z_2}}1 \mathrm dx \mathrm dy\\ &=\int_0^a \pi az \mathrm dz + \int_a^{2a} \pi (2a-z)^2 \mathrm dz\\ &=\dfrac{5}{6}\pi a^3 \end{align*}
      • 柱坐标

        I=02πdθ0aρdρ1aρ22aρ1dz=2π0aρ(2aρ1aρ2)dρ=56πa3\begin{align*} I&=\int_0^{2\pi}\mathrm d\theta \int_0^a\rho \mathrm d\rho \int_{\frac{1}{a}\rho^2}^{2a-\rho}1\cdot \mathrm dz\\ &=2\pi\int_0^a \rho(2a-\rho-\dfrac{1}{a}\rho^2)\mathrm d\rho\\ &=\dfrac{5}{6}\pi a^3 \end{align*}
  2. I=Ωzx2+y2dVI=\iiint_\Omega z\sqrt{x^2+y^2} \mathrm dV,其中 Ω:x2+y2=2x\Omega: x^2+y^2=2xz=xz=xxOyxOy 面所围成的闭区域

    Dxy:x2+y22xD_{xy}:x^2+y^2\le 2x

    I=π2π2dθ02cosθρdρ0ρcosθzρdz=π2π2dθ02cosθρρ2cos2θ022ρdρ=π2π212cos2θ15(2cosθ)5dθ=π2π2165cos7θ=3250π2cos7θdθ=512175\begin{align*} I&=\int_{-\frac{\pi}{2}}^{\frac{\pi}{2}}\mathrm d\theta \int_0^{2\cos\theta}\rho \mathrm d\rho \int_0^{\rho\cos\theta} z\rho \mathrm dz\\ &=\int_{-\frac{\pi}{2}}^{\frac{\pi}{2}}\mathrm d\theta \int_0^{2\cos\theta}\rho \cdot \dfrac{\rho^2\cos^2\theta-0^2}{2} \rho \mathrm d\rho\\ &=\int_{-\frac{\pi}{2}}^{\frac{\pi}{2}}\dfrac{1}{2}\cos^2\theta \dfrac{1}{5}(2\cos \theta)^5 \mathrm d\theta\\ &=\int_{-\frac{\pi}{2}}^{\frac{\pi}{2}}\dfrac{16}{5}\cos^7\theta\\ &=\dfrac{32}{5}\int_0^{\frac{\pi}{2}}\cos^7\theta \mathrm d\theta\\ &=\dfrac{512}{175} \end{align*}
  3. Ωf(x,y,z)dV\iiint_\Omega f(x,y,z)\mathrm dV 化为柱坐标系下三次积分。Ω:z=4x2y2\Omega:z=\sqrt{4-x^2-y^2}z=13x2+y2z=\dfrac{1}{3}\sqrt{x^2+y^2} 所围成的闭区域

    {z=4x2y2z=13(x2+y2){x2+y23Dxyz=0\begin{cases}z=\sqrt{4-x^2-y^2}\\z=\dfrac{1}{3}(x^2+y^2)\end{cases}\rightarrow \begin{cases}x^2+y^2\le3 \rightarrow D_{xy}\\z=0\end{cases}

    θ:02π\theta:0\to 2\piρ:03\rho:0\to \sqrt{3}z:13ρ24ρ2z:\dfrac{1}{3}\rho^2\to \sqrt{4-\rho^2}

    I=02πdθ03ρdρ13ρ24ρ2f(x,y,z)dz\displaystyle I=\int_0^{2\pi}\mathrm d\theta\int_0^{\sqrt{3}}\rho \mathrm d\rho \int_{\frac{1}{3}\rho^2}^{\sqrt{4-\rho^2}}f(x,y,z) \mathrm dz

  4. Ωf(x,y,z)dV\iiint_\Omega f(x,y,z)\mathrm dV 化为柱坐标系下三次积分。Ω:x2+y2+z2=4\Omega:x^2+y^2+z^2=4z=1z=-1 所围成的包含球心的区域

    Dxy:x2+y24D_{xy}:x^2+y^2\le 4

    {x2+y2+z2=4z=1{x2+y2=3z=0\begin{cases}x^2+y^2+z^2=4\\z=-1\end{cases}\rightarrow \begin{cases}x^2+y^2=3\\z=0\end{cases}

    θ:02π\theta:0\to 2\piρ:02\rho:0\to 2z:14ρ2(Dxy:x2+y23),4ρ24ρ2(Dxy:3x2+y24)z:-1\to \sqrt{4-\rho^2}(D_{x'y}:x^2+y^2\le 3),-\sqrt{4-\rho^2}\to \sqrt{4-\rho^2}(D_{x''y}:3\le x^2+y^2\le 4)

    \displaystyle I=\int_0^{2\pi} \mathrm d\theta \int_0^\sqrt{3}\rho \mathrm d\rho \int_{-1}^{\sqrt{4-\rho^2}}f(x,y,z)\mathrm dz + \int_0^{2\pi}\mathrm d\theta \int_{\sqrt{3}}^2 \rho \mathrm d\rho \int_{-\sqrt{4-\rho^2}}^{\sqrt{4-\rho^2}}f(x,y,z)\mathrm dz

球坐标系#

适用于 Ω\Omega 与锥面或球面有关,被积函数 x2+y2+z2x^2+y^2+z^2

(ϕ=zOM,θ=xOM,r=OM){x=rsinϕcosθy=rsinϕsinθz=rcosϕ\begin{align*} (\phi=\angle zOM,\theta=\angle xOM',r=|\vec{OM}|)\\ \begin{cases} x=r\sin\phi\cos\theta\\ y=r\sin\phi\sin\theta\\ z=r\cos\phi \end{cases} \end{align*}

dV=r2sinϕdϕdθdr\mathrm dV=r^2\sin\phi \mathrm d\phi \mathrm d\theta \mathrm dr

I=Ωf(rsinϕcosθ,rsinϕsinθ,rcosϕ)r2sinϕdϕdθdr=αβdθa1(θ)a2(θ)sinϕdϕb1(ϕ,θ)b2(ϕ,θ)f(ϕ,θ,r)r2dr\begin{align*} I&=\iiint_\Omega f(r\sin\phi\cos\theta, r\sin\phi\sin\theta, r\cos\phi)r^2\sin\phi \mathrm d\phi \mathrm d\theta \mathrm dr\\ &=\int_\alpha^\beta \mathrm d\theta \int_{a_1(\theta)}^{a_2(\theta)}\sin\phi \mathrm d\phi\int_{b_1(\phi,\theta)}^{b_2(\phi,\theta)}f(\phi,\theta,r)r^2 \mathrm dr \end{align*}

例题#

  1. I=Ωx2+y2+z2dVI=\iiint_\Omega \sqrt{x^2+y^2+z^2}\mathrm dVΩ:x2+y2+z2=z\Omega: x^2+y^2+z^2=z 所围成的闭区域

    θ:02π\theta:0\to 2\piφ:0π2\varphi:0\to \dfrac{\pi}{2}r:0cosφr:0\to \cos\varphi

    I=02πdθ0π2sinφdφ0cosφrr2dr=π10\begin{align*} I&=\int_0^{2\pi} \mathrm d\theta \int_0^{\frac{\pi}{2}}\sin\varphi \mathrm d\varphi \int_0^{\cos\varphi}r\cdot r^2 \mathrm dr\\ &=\dfrac{\pi}{10} \end{align*}
  2. I=Ω(x+y+z)e(x2+y2+z2)dVI=\iiint_\Omega (x+y+z)e^{-(x^2+y^2+z^2)}\mathrm dV,其中 Ω:x2+y2+z2=1\Omega:x^2+y^2+z^2=1z0z\ge 0(上半球体)

    解:Ω\Omega 关于 xOzxOzyOzyOz 面对称

    I=Ωze(x2+y2+z2)dV=02πdθ0π2sinφdφ01rcosφer2r2dr=π(121e)\begin{align*} I&=\iiint_\Omega ze^{-(x^2+y^2+z^2)}\mathrm dV\\ &=\int_0^{2\pi}\mathrm d\theta \int_0^{\frac{\pi}{2}}\sin \varphi \mathrm d\varphi \int_0^1 r\cos\varphi e^{-r^2}r^2 \mathrm dr\\ &=\pi(\dfrac{1}{2}-\dfrac{1}{e}) \end{align*}
  3. I=Ωx2dVI=\iiint_\Omega x^2 \mathrm dV,其中 Ω:x2+y2+z21\Omega:x^2+y^2+z^2\le 1

解:Ω\Omega 具有轮换对称性

I=Ωy2dV=Ωz2dV=13Ω(x2+y2+z2)dV=13(02πdθ)(0πsinφdφ)(01r2r2dr)=132π215=4π15\begin{align*} I&=\iiint_\Omega y^2 \mathrm dV\\ &=\iiint_\Omega z^2 \mathrm dV\\ &=\dfrac{1}{3}\iiint_\Omega (x^2+y^2+z^2)\mathrm dV\\ &=\dfrac{1}{3}(\int_0^{2\pi}\mathrm d\theta)(\int_0^{\pi}\sin\varphi \mathrm d\varphi)(\int_0^1 r^2\cdot r^2 \mathrm dr)\\ &=\dfrac{1}{3}\cdot 2\pi \cdot 2\cdot \dfrac{1}{5}\\ &=\dfrac{4\pi}{15} \end{align*}

三重积分换元公式#

{x=x(u,v,w)y=y(u,v,w)z=z(u,v,w)\begin{align*} \begin{cases} x=x(u,v,w)\\ y=y(u,v,w)\\ z=z(u,v,w) \end{cases} \end{align*}I=Ωf(x,y,z)dxdydz=Ωf(x(u,v,w),y(u,v,w),z(u,v,w))Jdudvdw\begin{align*} I&=\iiint_\Omega f(x,y,z)\mathrm dx \mathrm dy \mathrm dz\\ &=\iiint_\Omega f(x(u,v,w),y(u,v,w),z(u,v,w))|J| \mathrm du \mathrm dv \mathrm dw \end{align*}

其中 J(u,v,w)=(x,y,z)(u,v,w)=xuxvxwyuyvywzuzvzwJ(u,v,w)=\dfrac{\partial (x,y,z)}{\partial (u,v,w)}=\begin{vmatrix}x_u & x_v & x_w\\y_u & y_v & y_w\\ z_u & z_v & z_w\end{vmatrix}

  • 广义柱坐标变换 {x=aρcosθy=bρsinθz=z\begin{cases}x=a\rho\cos\theta\\y=b\rho\sin\theta\\z=z\end{cases}dV=abρdρdθdz\mathrm dV=ab\rho \mathrm d\rho \mathrm d\theta \mathrm dzDxyD_{xy} 与椭圆有关,Ω\Omega 与椭球面有关)
  • 广义球坐标变换 {x=arsinφcosθy=brsinφsinθz=crcosφ\begin{cases}x=ar\sin\varphi\cos\theta\\y=br\sin\varphi\sin\theta\\z=cr\cos\varphi\end{cases}dV=abcr2sinφdφdθdr\mathrm dV=abcr^2\sin\varphi \mathrm d\varphi \mathrm d\theta \mathrm dr

例题#

x2a2+y2b2+z2c21\dfrac{x^2}{a^2}+\dfrac{y^2}{b^2}+\dfrac{z^2}{c^2}\le 1 的体积

解:

I=Ω1dV=02πdθ0πsinφdφ011abcr2dr=43πabc\begin{align*} I&=\iiint_\Omega 1\cdot \mathrm dV\\ &=\int_0^{2\pi} \mathrm d\theta \int_0^\pi \sin\varphi \mathrm d\varphi \int_0^1 1\cdot abc\cdot r^2 \mathrm dr\\ &=\dfrac{4}{3}\pi abc \end{align*}
9.3 三重积分
https://gitee.com/jason_ren/advanced-math-note
作者
Jason Ren
发布于
2025-12-03
许可协议
CC BY-SA 4.0

部分信息可能已经过时

封面
Music
Artist
封面
Music
Artist
0:00 / 0:00