赞
踩
我们继续积分学的第四个模块——三重积分。讲实话,这一块其实不太好理解,主要我空间想象能力不太行。不过最主要的考点应该还是计算,而且其实三重积分算到后面还是算定积分和二重积分,所以,前面的内容要先掌握好。
和二重积分很像,三重积分也是用来求解多元不规则量,下面从一个具体问题来引入。
设 Ω \Omega Ω 为几何体,其体密度为 ρ ( x , y , z ) \rho(x,y,z) ρ(x,y,z) ,求其质量 m m m 。如果每个点处的密度是一样的,那么, m = ρ V m=\rho V m=ρV 。但是现在各点处的密度不一,因此我们采用微元法的思想。
(1)将几何体划分为小体积元素 Δ v 1 , Δ v 2 , … , Δ v n . \Delta v_1,\Delta v_2,\dots,\Delta v_n. Δv1,Δv2,…,Δvn.
(2)在小体积元素中任取一点 ( ξ i , η i , ζ i ) ∈ Δ v i (\xi_i,\eta_i,\zeta_i) \in \Delta v_i (ξi,ηi,ζi)∈Δvi ,则每一小体积元素的质量为 Δ m i ≈ ρ ( ξ i , η i , ζ i ) Δ v i \Delta m_i \approx \rho (\xi_i,\eta_i,\zeta_i)\Delta v_i Δmi≈ρ(ξi,ηi,ζi)Δvi ,整个几何体的质量可近似表示为 m = ∑ i = 1 n ρ ( ξ i , η i , ζ i ) Δ v i ; m=\sum_{i=1}^n \rho (\xi_i,\eta_i,\zeta_i) \Delta v_i; m=i=1∑nρ(ξi,ηi,ζi)Δvi; (3)设 λ = m a x ( d i ) , d i \lambda =max(d_i),d_i λ=max(di),di 为小体积元素 Δ v i \Delta v_i Δvi 的直径,则该几何体精确的质量为 m = lim λ → 0 ∑ i = 1 n ρ ( ξ i , η i , ζ i ) Δ v i . m=\lim_{\lambda \to 0} \sum_{i=1}^n \rho (\xi_i,\eta_i,\zeta_i) \Delta v_i. m=λ→0limi=1∑nρ(ξi,ηi,ζi)Δvi.
将上述实际问题抽象出一个一般模型,可得到三重积分的定义如下:
设 Ω \Omega Ω 为几何体,三元函数 f ( x , y , z ) f(x,y,z) f(x,y,z) 在该几何体上有界。
(1)将几何体划分为小体积元素 Δ v 1 , Δ v 2 , … , Δ v n . \Delta v_1,\Delta v_2,\dots,\Delta v_n. Δv1,Δv2,…,Δvn.
(2)在小体积元素中任取一点 ( ξ i , η i , ζ i ) ∈ Δ v i (\xi_i,\eta_i,\zeta_i) \in \Delta v_i (ξi,ηi,ζi)∈Δvi ,作和 ∑ i = 1 n f ( ξ i , η i , ζ i ) Δ v i ; \sum_{i=1}^n f (\xi_i,\eta_i,\zeta_i) \Delta v_i; i=1∑nf(ξi,ηi,ζi)Δvi; (3)设 λ = m a x ( d i ) , d i \lambda =max(d_i),d_i λ=max(di),di 为小体积元素 Δ v i \Delta v_i Δvi 的直径,若极限 lim λ → 0 ∑ i = 1 n f ( ξ i , η i , ζ i ) Δ v i . \lim_{\lambda \to 0} \sum_{i=1}^n f (\xi_i,\eta_i,\zeta_i) \Delta v_i. λ→0limi=1∑nf(ξi,ηi,ζi)Δvi. 存在,则称函数 f ( x , y , z ) f(x,y,z) f(x,y,z) 在区域 Ω \Omega Ω 上可积。该极限的值称为函数 f ( x , y , z ) f(x,y,z) f(x,y,z) 在区域 Ω \Omega Ω 上的三重积分,记为 ∭ Ω f ( x , y , z ) d v . \iiint_{\Omega}f(x,y,z)dv. ∭Ωf(x,y,z)dv.
一般性质和二重积分类似,区域,常数,被积函数为 1 时,积分中值定理等等。
同样,三重积分也有奇偶性和对称性。我就列举一个,其余两个是同理的。
这个方法我认为相当直观,设区域 Ω \Omega Ω 表示为 { ( x , y , z ) ∣ ( x , y ) ∈ D z , c ≤ z ≤ d (x,y,z)|(x,y) \in D_z,c \leq z \leq d (x,y,z)∣(x,y)∈Dz,c≤z≤d} ,如下图所示。可以像二重积分那样,作一个向右的箭头。
则有如下计算方法:
∭
Ω
f
(
x
,
y
,
z
)
d
v
=
∫
c
d
d
z
∬
D
z
f
(
x
,
y
,
z
)
d
x
d
y
.
\iiint_{\Omega}f(x,y,z)dv=\int_c^ddz \iint_{D_z}f(x,y,z)dxdy.
∭Ωf(x,y,z)dv=∫cddz∬Dzf(x,y,z)dxdy. 举个例子:计算
∭
Ω
(
z
2
+
2
x
y
)
d
v
\iiint_{\Omega}(z^2+2xy)dv
∭Ω(z2+2xy)dv ,其中
Ω
\Omega
Ω 为锥面
z
=
x
2
+
y
2
z=\sqrt{x^2+y^2}
z=x2+y2
和
z
=
2
z=2
z=2 所围成的几何体,如下图所示。绿色的表示锥面,粉色的表示
z
=
2
z=2
z=2 平面。
首先,我们同样检查对称性,发现区域是关于
x
O
z
xOz
xOz 平面对称的,且
−
2
x
y
=
2
x
(
−
y
)
-2xy=2x(-y)
−2xy=2x(−y) 。故原积分可以拆成两部分,第一部分为
∭
Ω
z
2
d
v
\iiint_{\Omega}z^2dv
∭Ωz2dv 。第二部分为
∭
Ω
2
x
y
d
v
\iiint_{\Omega}2xydv
∭Ω2xydv ,因为对称性,第二部分为 0 。故只计算第一部分即可。
利用切片法,原区域可以表示为 { ( x , y , z ) ∣ 0 ≤ z ≤ 2 , ( x , y ) ∈ D z (x,y,z)|0 \leq z \leq 2,(x,y) \in D_z (x,y,z)∣0≤z≤2,(x,y)∈Dz} ,其中 D z D_z Dz 可表示为 { ( x , y ) ∣ x 2 + y 2 = z 2 (x,y)|x^2+y^2=z^2 (x,y)∣x2+y2=z2} 。平面区域的二重积分可采用极坐标进行简化计算,故原积分为 I = ∫ 0 2 d z ∬ D z z 2 d x d y = ∫ 0 2 d z ∫ 0 2 π d θ ∫ 0 z z 2 r d r = 32 π 5 I=\int_0^2dz\iint_{D_z}z^2dxdy=\int_0^2dz\int_0^{2\pi}d\theta \int_0^{z}z^2rdr=\frac{32\pi}{5} I=∫02dz∬Dzz2dxdy=∫02dz∫02πdθ∫0zz2rdr=532π
铅直投影法和切片法顺序是相反的,它先往 x O y xOy xOy 平面投影,先一后二。
设区域 Ω \Omega Ω 表示为 { ( x , y , z ) ∣ ( x , y ) ∈ D , φ 1 ( x , y ) ≤ z ≤ φ 2 ( x , y ) (x,y,z)|(x,y) \in D,\varphi_1(x,y) \leq z \leq \varphi_2(x,y) (x,y,z)∣(x,y)∈D,φ1(x,y)≤z≤φ2(x,y)} ,如下图所示。
则有如下计算方法:
∭
Ω
f
(
x
,
y
,
z
)
d
v
=
∬
D
d
x
d
y
∫
φ
1
(
x
,
y
)
φ
2
(
x
,
y
)
f
(
x
,
y
,
z
)
d
z
.
\iiint_{\Omega}f(x,y,z)dv= \iint_{D}dxdy\int_{\varphi_1(x,y) }^{\varphi_2(x,y)} f(x,y,z)dz.
∭Ωf(x,y,z)dv=∬Ddxdy∫φ1(x,y)φ2(x,y)f(x,y,z)dz. 我们用这个方法,做一做上面那道例题。用此方法,则区域
Ω
\Omega
Ω 可表示为 {
(
x
,
y
,
z
)
∣
(
x
,
y
)
∈
D
,
0
≤
z
≤
x
2
+
y
2
(x,y,z)|(x,y) \in D,0 \leq z \leq \sqrt{x^2+y^2}
(x,y,z)∣(x,y)∈D,0≤z≤x2+y2
},其中
D
D
D 可表示为 {
(
x
,
y
)
∣
x
2
+
y
2
=
2
2
(x,y)|x^2+y^2=2^2
(x,y)∣x2+y2=22} 。平面区域的二重积分可采用极坐标进行简化计算,故原积分为
I
=
∫
0
2
π
d
θ
∫
0
2
r
d
r
∫
x
2
+
y
2
2
z
2
d
z
=
32
π
5
I=\int_0^{2\pi}d\theta \int _0^{2}rdr\int_{\sqrt{x^2+y^2}}^{2}z^2dz=\frac{32\pi}{5}
I=∫02πdθ∫02rdr∫x2+y2
2z2dz=532π
书上还说了一种啥柱面坐标法,其实就是直角坐标加了极坐标法,上面的题目咱就已经用了这种方法了,所以就不多说。我们主要看看球坐标是怎么个一回事。
球坐标主要是对原坐标系进行了一定变换,如下: x = r c o s θ s i n φ , y = r s i n θ s i n φ , z = r c o s φ . x=rcos\theta sin\varphi,y=rsin\theta sin\varphi,z=rcos\varphi. x=rcosθsinφ,y=rsinθsinφ,z=rcosφ.
这样,原来的三重积分就变成了如下的形式: ∭ Ω f ( x , y , z ) d v = ∫ α β d θ ∫ θ 1 θ 2 d φ ∫ r 1 ( θ , φ ) r 2 ( θ , φ ) f ( x , y , z ) r 2 s i n φ d r . \iiint_{\Omega}f(x,y,z)dv=\int_{\alpha}^{\beta}d\theta \int_{\theta_1}^{\theta_2}d\varphi\int_{r_1(\theta,\varphi)}^{r_2(\theta,\varphi)}f(x,y,z)r^2sin\varphi dr. ∭Ωf(x,y,z)dv=∫αβdθ∫θ1θ2dφ∫r1(θ,φ)r2(θ,φ)f(x,y,z)r2sinφdr. 关键是如何找到三个参数的范围,而要找到三个参数的范围,就得了解这三个参数都代表什么意思。下图(图源网络)可以帮助我们理解。
r r r 是代表几何体上一点到原点距离,可以像极坐标那样,从原点画一条射线出来,看看范围。 θ \theta θ 表示的是 r r r 在 x O y xOy xOy 平面上的投影直线与 x x x 轴的夹角。 φ \varphi φ 表示 r r r 与 z z z 轴的夹角。我们用一个例题来试试。
计算
∭
Ω
(
x
2
+
y
2
)
d
v
\iiint_{\Omega}(x^2+y^2)dv
∭Ω(x2+y2)dv ,其中
Ω
=
\Omega=
Ω= {
(
x
,
y
,
z
)
∣
x
2
+
y
2
+
z
2
≤
1
,
z
≥
0
(x,y,z)|x^2+y^2+z^2 \leq 1,z \geq 0
(x,y,z)∣x2+y2+z2≤1,z≥0}. 我们很容易看出,该几何体是一个上半球,半径为 1 。
通过上图,我们可以找到三个参数的范围,如下: 0 ≤ θ ≤ 2 π , 0 ≤ φ ≤ π 2 , 0 ≤ r ≤ 1 0 \leq \theta \leq2\pi,0 \leq \varphi \leq \frac{\pi}{2},0 \leq r \leq 1 0≤θ≤2π,0≤φ≤2π,0≤r≤1 则原积分可计算为 I = ∫ 0 2 π d θ ∫ 0 π 2 d φ ∫ 0 1 r 2 s i n 2 φ ⋅ r 2 s i n φ d r = ∫ 0 2 π d θ ∫ 0 π 2 s i n 3 φ d φ ∫ 0 1 r 4 d r = 4 π 15 I=\int_0^{2\pi}d\theta \int_0^{\frac{\pi}{2}}d\varphi \int_0^1 r^2sin^2 \varphi \cdot r^2 sin\varphi dr=\int_0^{2\pi}d\theta \int_0^{\frac{\pi}{2}}sin^3 \varphi d\varphi \int_0^1r^4dr=\frac{4\pi}{15} I=∫02πdθ∫02πdφ∫01r2sin2φ⋅r2sinφdr=∫02πdθ∫02πsin3φdφ∫01r4dr=154π
应用我这里不细说,放到和前面的那些应用一起放到后面的空间几何来。三重积分应用主要有:
那到此,积分学的内容也将告一段落了,不过仅仅是理论部分,还需要大量的练习才能真正在考场上稳稳地拿下。
Copyright © 2003-2013 www.wpsshop.cn 版权所有,并保留所有权利。