1356 字
7 分钟
微积分 08|多重积分、换元与积分次序

数学总导航 · 本课程目录 · ← 上一章 · 下一章 →

一张人口密度图,怎样算全城人口?#

地图上的颜色表示每平方公里住多少人。只看最红的区域不能得到总人口,必须用每小块的“密度 × 面积”相加。

例如矩形城区长 2 公里、宽 1 公里,密度为 ρ(x,y)=1000+500x\rho(x,y)=1000+500x 人/平方公里。那么人口为

∫01∫02(1000+500x)dxdy=3000 人.\int_0^1\int_0^2(1000+500x)dxdy=3000\text{ 人}.

若把地图坐标改成圆心距离与角度,相同的距离格和角度格并不对应相同面积:越远离圆心,一格越宽。这就是极坐标积分要乘 rr 的原因;坐标格不是土地面积本身。

把人口统计真的拆成小格#

把城区分成四块,每块面积 0.5 平方公里。如果四块平均密度分别为 1000、1250、1500、1750,那么总人口估计为 0.5(1000+1250+1500+1750)=27500.5(1000+1250+1500+1750)=2750。

更精细的模型不依赖这四个估计,而是让网格不断变小,每格用当地密度乘面积。二维积分的本质仍是“速率式的量 × 对应小体积,再累加”,只是现在需要两个方向的位置。

非矩形区域为什么会出现随变量改变的上下限#

假设一片三角地界为 0≤y≤x≤10\le y\le x\le1。固定一个横坐标 x=0.2x=0.2,竖条长度是 0.2;固定 x=0.8x=0.8,长度是 0.8。所以先沿竖条积分时,上限必须是当前的 xx,不能固定成 1,否则会多算三角形外的土地。

换积分顺序时,先固定 yy,允许的横坐标变成 y≤x≤1y\le x\le1。因此同一片地可写为

∫01∫0xf(x,y)dydx=∫01∫y1f(x,y)dxdy,\int_0^1\int_0^x f(x,y)dydx =\int_0^1\int_y^1 f(x,y)dxdy,

这里先假设 ff 连续。画边界、固定一条截线、找到它穿过的端点,比机械交换 dx,dydx,dy 可靠。

为什么远处的一度比近处的一度覆盖更多土地#

在半径 rr 处,角宽 Δθ\Delta\theta 的小扇环沿圆周的宽度约为 rΔθr\Delta\theta,径向厚度为 Δr\Delta r,小面积约为 rΔrΔθr\Delta r\Delta\theta。这个多出来的 rr 不是换元公式的装饰,而是由实际土地宽度决定的。

同理,普通地图的一个经纬度小格在不同纬度也不代表相同地表面积。积分坐标可以任意选,但每格对应的实际权重必须正确。Jacobian 行列式是把这种局部体积修正统一写出的办法。

带你做:随机落点离圆心多远#

在单位圆盘内均匀随机选一点。离圆心的距离 RR 在 0 与 1 之间,但它不是均匀分布:外圈有更多面积。

半径不超过 rr 的概率等于小圆面积除以大圆面积,为 r2r^2。所以距离的密度是 2r2r,平均距离为

ER=∫01r⋅2r dr=23.\mathbb ER=\int_0^1r\cdot2r\,dr=\frac23.

若错误地把半径当成均匀分布,会得到 1/21/2,相当于过度偏爱靠近圆心的地方。用极坐标再算一次:均匀面积密度为 1/π1/\pi,平均距离是 π−1∫02π∫01r⋅r drdθ=2/3\pi^{-1}\int_0^{2\pi}\int_0^1r\cdot r\,drd\theta=2/3。

两个 rr 各有工作:一个是你要平均的距离,另一个是坐标换元后的面积权重。这道题把几何、Jacobian 和概率密度连在同一次计算里。

回头整理:定义、工具与练习#

区域上的累积#

对矩形 D=[a,b]×[c,d]D=[a,b]\times[c,d] 上连续的函数,可把区域分成小矩形,令面积和的极限定义二重积分。Fubini 定理在此情形给出

∬Df(x,y) dx dy=∫ab(∫cdf(x,y) dy)dx.\iint_Df(x,y)\,dx\,dy =\int_a^b\left(\int_c^df(x,y)\,dy\right)dx.

对非矩形区域,先画出边界。例如三角形 D={0≤x≤1,0≤y≤x}D=\{0\le x\le1,0\le y\le x\},积分 x+yx+y 得

∫01∫0x(x+y) dy dx=∫0132x2dx=12.\int_0^1\int_0^x(x+y)\,dy\,dx =\int_0^1\tfrac32x^2dx=\tfrac12.

换元为什么要乘行列式#

光滑一一坐标变换 x=Φ(u)x=\Phi(u) 在局部近似于线性映射 DΦD\Phi,小体积会放大 ∣det⁡DΦ∣|\det D\Phi| 倍。因此在满足换元定理的区域上,

∫Φ(U)f(x) dx=∫Uf(Φ(u))∣det⁡DΦ(u)∣ du.\int_{\Phi(U)}f(x)\,dx =\int_U f(\Phi(u))|\det D\Phi(u)|\,du.

通常要求在开区域内为 C1C^1 微分同胚;边界和可忽略的退化集合可另行处理。绝对值对应体积非负,而非有向面积。

极坐标 x=rcos⁡θ,y=rsin⁡θx=r\cos\theta,y=r\sin\theta 的 Jacobian 为 rr,于是单位圆盘面积为 ∫02π∫01r dr dθ=π\int_0^{2\pi}\int_0^1r\,dr\,d\theta=\pi。原点与角度接缝不满足全局一一性,但属于可忽略集合,可用分区处理。

概率密度也必须换体积#

若 X=Φ(U)X=\Phi(U),则在可逆、光滑且非退化的区域上

pX(x)=pU(Φ−1(x))∣det⁡DΦ−1(x)∣.p_X(x)=p_U(\Phi^{-1}(x))|\det D\Phi^{-1}(x)|.

这是“概率质量守恒”。例如二维标准正态密度转为极坐标后含 re−r2/2/(2π)r e^{-r^2/2}/(2\pi),积分时不能忘记 rr。

积分次序的准确边界#

在通常的 σ\sigma 有限测度空间中:Tonelli 定理允许对非负可测函数交换积分次序,即使结果是正无穷;Fubini 定理要求绝对可积,保证有符号积分可交换且几乎处处的截面可积。

非负并不等于积分有限;绝对可积也不能仅靠检查两个迭代积分恰好存在来代替。泛函分析中的积分算子和概率中的条件期望都会用到这个区别。

练习与解答#

  1. 求单位圆盘上 x2+y2x^2+y^2 的积分。
  2. X,YX,Y 独立且可积时,为何 E[XY]=EXEY\mathbb E[XY]=\mathbb EX\mathbb EY?
  3. 积分非负函数时,交换次序后得到无穷是否说明定理失败?

解答:1. ∫02π∫01r3drdθ=π/2\int_0^{2\pi}\int_0^1r^3drd\theta=\pi/2。2. 独立使联合分布为乘积;先由 Tonelli 得 E∣XY∣=E∣X∣E∣Y∣<∞\mathbb E|XY|=\mathbb E|X|\mathbb E|Y|<\infty,再由 Fubini 分离。3. 不失败,Tonelli 允许扩展实数值 +∞+\infty。


数学总导航 · 本课程目录 · ← 上一章 · 下一章 →

微积分 08|多重积分、换元与积分次序
https://f1yingwhite.github.io/posts/数学/微积分/08-多重积分、换元与积分次序/
作者
F1yingWhite
发布于
2026-11-03
许可协议
CC BY-NC-SA 4.0