三重积分:空间区域中的累积 | 多变量微积分 I | 自在学三重积分:空间区域中的累积
二重积分把平面区域切成小矩形,再把每个小矩形上的量加起来。三重积分做的是同一件事,只是区域从一张平面变成一个空间实体,小块从小矩形变成小长方体。
如果 E 是空间中的实体,f(x,y,z) 是定义在实体内的量,那么
∭Ef(x,
y
,
z
)
d
V
表示把 E 内每一点附近的小体积所携带的量累加起来。这里的 dV 是体积微元。在直角坐标中,它对应一个很小的长方体体积。
三重积分的基本画面:空间区域、取样点、密度和小体积一起决定每一小块的贡献。
从小长方体开始
B=[a,b]×[c,d]×[e,h]. 把 [a,b] 分成 m 段,把 [c,d] 分成 n 段,把 [e,h] 分成 p 段。这样 B 被切成 mnp 个小长方体。每个小长方体的体积是
ΔV=ΔxΔyΔz. 在第 (i,j,k) 个小长方体中取一个点 (xijk∗,yijk∗,zijk∗),用函数值近似这一小块上的量。所有小块加起来,得到三重黎曼和
i=1∑mj=1∑nk=1∑pf(xijk∗,yijk∗,zijk∗)ΔV. 当三个方向的网格都不断变细,这个和的极限就是三重积分。
∭Bf(x,y,z)dV=maxΔV→0lim∑f(xijk∗,yijk∗,zijk∗)ΔV. 三重积分不是“多写一个积分号”这么简单。它表示累积对象从平面面积上的小片变成了空间体积中的小块。如果 f 是密度,fdV 是一小块质量;如果 f 是温度,∭EfdV 先得到温度对体积的总累积,再除以体积才是平均温度。
盒状区域与迭代积分
盒状区域最像一维定积分的直接推广。若 f 在盒状区域 B=[a,b]×[c,d]×[e,h] 上连续,就可以把三重积分写成迭代积分。
∭Bf(x,y,z)dV=∫ab∫cd∫ehf(x,y,z)dzdydx. 最内层的 dz 表示先固定 x 和 y,沿着竖直方向累加;得到一个关于 x,y 的函数后,再沿 y 累加;最后沿 x 累加。
盒状区域中,六种积分次序都可以使用。选择顺序时看计算是否更顺手。
例题:计算盒状区域上的三重积分
∭B(x+2y+3z)dV,B=[0,2]×[1,3]×[0,1]. 先把区域写成三个独立区间。这里 x 从 0 到 2,y 从 1 到 3,z 从 0 到 1,所以可以直接写成
∫02∫13∫0对 z 积分时,把 x 和 y 当常数。内层积分为
∫再对 y 积分:
∫13(x+2y最后对 x 积分:
∫02(2x+11)dx=这里的数值 26 不一定是体积。只有当被积函数恒为 1 时,三重积分才直接给出区域体积。
一般空间区域:先找影子,再找高度
真实题目中的空间区域很少刚好是盒子。更常见的情况是:实体在 xy 平面上有一个投影区域 D,而每个 (x,y) 上方有一段竖直线段。
E={(x,y,z):(x,y)∈D, g1(x,y)≤z≤g2(x,y)}, ∭Ef(x,y,z)dV=∬D(∫g1(x,y)g2(x,y)f(x,y,z)ddA. 这句话的读法很具体:先在平面上确定一点 (x,y),再沿着通过这一点的竖直线段从下表面积到上表面,最后把所有竖直线段在投影区域 D 上累加。
先找投影区域,再找穿过实体的线段上下端点,是写三重积分限的常用方法。
例题:第一卦限中平面下方的累积
2x+3y+z=6 ∭E2xdV. 第一卦限给出 x≥0、y≥0、z≥0。平面下方给出
0≤z≤6−2x−3y.要让上界非负,需要 2x+3y≤6。
把实体投影到 xy 平面,得到三角形区域
D={(x,y):0≤按照“先 z,再 y,再 x”写出积分:
先对 z 积分,得到竖直线段长度乘以 2x:
∫0继续计算外两层:
∫03∫0写一般空间区域的积分限时,最容易错的是把三维边界方程直接混在同一层。先决定最内层变量,再问“固定外层变量以后,它从哪里到哪里”,这样每一层的上下限才有清楚含义。
积分次序:先画切片,再写符号
dxdydz,dxdzdy,dydxdz,dydzdx,dzdxdy,dzdydx. 这些写法不是把字母随意排列。每一种次序都对应一种切片方式。若最内层是 dz,通常先看 xy 平面的投影;若最内层是 dx,通常先看 yz 平面的投影;若最内层是 dy,通常先看 xz 平面的投影。
以第一卦限内平面 x+y+z=4 下方的四面体为例。如果先沿 z 方向累加,可以写成
∫04∫04−x∫04−x−yf(x,y,z)dzdydx. 如果改成最内层先对 x 积分,就要看 yz 平面的投影。投影满足 y≥0、z≥0、y+z≤4,固定 (y,z) 后,
0≤x≤4−y−z. ∫04∫04−y∫04−y−zf(x,y,z)dxdzdy. 判断一个积分限是否合理,可以做两个快速检查。第一,最内层上下限只能含外层变量,不能含正在积分的变量。第二,每一层固定后,剩下的范围应该是一段线段、一个平面区域或一个区间,而不是还没有说清楚的三维实体。
体积、质量、平均值和质心
三重积分最直接的应用是体积。若 E 是空间实体,那么
V(E)=∭E1dV. 若实体密度为 ρ(x,y,z),质量是
M=∭Eρ(x,y,z)dV. 如果 f(x,y,z) 表示温度、浓度、电荷密度以外的某种场量,它在 E 上的平均值是
favg=V(E)1∭Ef(x,y,z)dV. 若 ρ(x,y,z) 是质量密度,质心坐标为
xˉ=M1∭Exρ(x,y,z)dV,yˉ=M1∭Eyρ(x,y,z)dV,zˉ=M1∭Ezρ(x,y,z)dV. 质心由质量分布决定。密度不均匀时,它不一定在几何中心。
例题:密度变化时质心会移动
E=[0,1]×[0,2]×[0,1] ρ(x,y,z)=1+x. 先计算质量。由于密度只依赖 x,y 和 z 方向只贡献区间长度:
M=∫01∫02∫01(1+x)dzdydx=3.计算 xˉ 的分子:
∭E由于密度不依赖 y 和 z,在这两个方向上仍然保持对称,所以
yˉ=1,质心是
(95,1,21
用截面理解空间区域
有些实体不好直接想象成“竖直线段堆起来”,但可以看成一张张截面叠起来。固定 z,得到截面面积 A(z),再把这些截面沿 z 方向累加,就得到体积:
V=∫mnA(z)dz. 这其实也是三重积分的切片思想。只是内层两重积分已经被概括成了截面面积。
固定一个变量得到横截面,是理解空间区域和估算体积的有效方式。
z=4−x2−y2 x2+y2≤4−z. 所以截面是半径 4−z 的圆,面积为
A(z)=π(4−z),0≤z≤4. V=∫04π(4−z)dz=8π. 这里没有把三重积分完全展开成 dxdydz 的形式,但思想相同:先把每一层的累积看清楚,再沿高度方向累加。
常见误区
把投影区域当成实体本身
投影区域 D 只是实体在某个坐标平面上的影子。它负责外两层变量的范围,不包含最内层变量的上下界。若区域写成
E={(x,y,z):(x,y)∈D, g1(x,y)≤z≤g2(x,y)}, 那么 D 只回答“哪些 (x,y) 会被用到”,不回答“z 到哪里为止”。
忘记被积函数的单位
若 ρ 是质量密度,ρdV 是质量;若 f 是温度,fdV 不是温度,而是“温度乘体积”的累积。要得到平均温度,必须除以体积。
机械改变积分次序
改变次序时不能只移动 dx,dy,dz。上下限也必须跟着新的投影和切片方向重写。对一般区域来说,积分次序改变后,上下限通常会变。
如果你发现某一层积分的上界含有正在积分的变量,例如 ∫04−zdz,几乎可以确定限写错了。上界和下界只能描述“当前变量的范围”,它们可以依赖外层已经固定的变量,但不能依赖当前变量本身。
练习
练习 1
∭BxyzdV,B=[0,2]×[0,1]×[1,3]. 因为区域是盒状区域,可以分离计算:
∭BxyzdV=(∫02xdx)(∫01ydy)(∫13zdz)=2⋅21⋅4=4.练习 2
设 E 位于曲面 z=3+x−y 下方、xy 平面上方,且投影区域为矩形 0≤x≤2、0≤y≤1。写出计算体积的三重积分,并求体积。
因为 0≤z≤3+x−y,体积为
V=∫02∫01∫03+x−y1dzdydx.先对 z 积分,得到
V=∫02∫01(3+x−练习 3
设 E 是第一卦限内平面 x+y+z=1 下方的四面体。用 dzdydx 的次序写出 ∭Ef(x,y,z)dV,并求 E 的体积。
投影到 xy 平面后有 0≤x≤1、0≤y≤1−x。固定 (x,y) 后,
0≤z≤1−x−y.所以
∭Ef(x,y,z)dV=∫体积取 f=1:
V=∫01∫01−x(1−练习 4
∫04∫04−x∫04−x−yf(x,y,z)dzdydx 原区域是第一卦限内 x+y+z≤4 的四面体。若最内层对 x 积分,就看 yz 平面投影:
y≥0,z≥0,y+z≤4.固定 (y,z) 后,
0≤x≤4−y−z.因此可以写成
∫04∫04−y∫练习 5
单位立方体 0≤x,y,z≤1 的密度为 ρ(x,y,z)=2+z。求质量和质心。
质量为
M=∫01∫01∫01(2+z)dzdydx=25.由于密度只依赖 z,所以
xˉ=21,yˉ=再计算
zˉ=M1∫质心为
(21,21,1581
(
x
+
2y+
3z)dzdydx.
0
1
(
x
+
2y+
3z)dz=
x+
2y+
23.
+
23
)
d
y
=
2x+
11.
26.
z
)
x
≤
3, 0≤
y≤
2−
32x}.
∭E
2
x
d
V
=
∫03∫02−32x∫06−2x−3y2xdzdydx.
6
−
2
x
−
3
y
2
x
d
z
=
2x(6−
2x−
3y).
2−32x
2
x
(
6
−
2x−
3y)dydx=
9.
x
(
1
+
x)dV=
2∫01(x+
x2)dx=
35.
xˉ=35/3=95.
zˉ
=
21.
)
.
几何中心是 (1/2,1,1/2)。因为密度随着 x 增大而增大,质心向 x 较大的一侧移动。
y)dydx=
7.
01
∫01−x
∫01−x−y
f
(
x
,
y
,
z
)
d
z
d
y
d
x
.
x−
y)dydx=
61.
04−y−z
f
(
x
,
y
,
z
)
d
x
d
z
d
y
.
21
.
01
∫01
∫01
z
(
2
+
z)dzdydx=
5/21(1+31)=
158.
)
.