沿一个方向累计,再沿另一个方向累计
计算平面区域的面积时,我们常把它切成竖条,再把条带加起来。换成横条,区域没有改变,总面积似乎也不该变。但如果累计的量有正有负,或者某些方向的总量无穷,这个直觉需要条件才能成立。
这一章会把“长乘宽”推广到乘积测度,再说明 Tonelli 与 Fubini 分别允许怎样交换积分。计算之前,我建议你把区域画出来,同时写出可测性与非负性或绝对可积性;上下限和定理条件是两件都要做的事。
从可测矩形定义乘积空间
给定测度空间 ( X , Σ , μ ) (X,\Sigma,\mu) ( X , Σ , μ ) 与 ( Y , T , ν ) (Y,\mathcal T,\nu) ( Y , T , ν ) 。本章使用标准的 σ \sigma σ -有限 假设:每个空间都能写成可数个有限测度可测集的并。实数轴上的 Lebesgue 测度满足这一点,因为 R = ⋃ n [ − n , n ] \mathbb R=\bigcup_n[-n,n] R = ⋃ n [ − n , n ] ;自然数上的计数测度也满足,因为可用有限集合穷尽。
当 A ∈ Σ A\in\Sigma A ∈ Σ 、B ∈ T B\in\mathcal T B ∈ T 时,称 A × B A\times B A × B 为可测矩形。这里“矩形”并不要求两边是区间,A A A 或 B B B 可以是任意可测集。乘积 -代数定义为
Σ ⊗ T = σ ( { A × B : A ∈ Σ , B ∈ T } ) . \Sigma\otimes\mathcal T
=\sigma\bigl(\{A\times B:A\in\Sigma,\ B\in\mathcal T\}\bigr). Σ ⊗ T = σ ( { A × B : A ∈ Σ , B ∈ T } ) .
乘积测度应当保留的规则是
( μ × ν ) ( A × B ) = μ ( A ) ν ( B ) , (\mu\times\nu)(A\times B)=\mu(A)\nu(B), ( μ × ν ) ( A × B ) = μ ( A ) ν ( B ) ,
其中仍约定 0 ⋅ ∞ = 0 0\cdot\infty=0 0 ⋅ ∞ = 0 。例如一条竖直线即使向上无限延伸,二维 Lebesgue 测度仍为零。
这个测度怎样建立? 先把可测矩形的有限不交并收成一个集合代数 A \mathcal A A ,在其中定义
π ( ⨆ j = 1 m A j × B j ) = ∑ j = 1 m μ ( A j ) ν ( B j ) . \pi\left(\bigsqcup_{j=1}^m A_j\times B_j\right)
=\sum_{j=1}^m\mu(A_j)\nu(B_j). π ( j = 1 ⨆ m A j × B
两种表示可以交叉细分,有限可加性保证结果相同。但建立测度还需要可数可加,不只是有限可加。若一个矩形恰是可数个不交矩形的并,对每个固定 x x x ,对应的指示函数满足
1 A ( x ) 1 B ( y ) = ∑ j 1 A j ( x ) 1 B j ( y ) . \mathbf1_A(x)\mathbf1_B(y)
=\sum_j\mathbf1_{A_j}(x)\mathbf1_{B_j}(y). 1 A ( x ) 1 B ( y ) = j
先对 y y y 积分,用非负级数的单调收敛得到 1 A ( x ) ν ( B ) = ∑ j 1 A j ( x ) ν ( B j ) \mathbf1_A(x)\nu(B)=\sum_j\mathbf1_{A_j}(x)\nu(B_j) 1 A ( x ) ν ( B ) = ∑ j ;再对 用同样的步骤,就得到乘积数值的可数可加。对有限不交并细分后,结论推广到 ,所以 是前测度。
接下来重复第二章的覆盖思路:用 A \mathcal A A 中的集合覆盖任意 E ⊆ X × Y E\subseteq X\times Y E ⊆ X × Y ,取覆盖费用 ∑ j π ( R j ) \sum_j\pi(R_j) ∑ j π ( R 的下确界,得到外测度 。每个 都满足 Carathéodory 条件,因为覆盖可以按 与 切开,费用由有限可加性拆分。并且 :一个方向来自自覆盖,另一个方向把任意覆盖依次去掉重叠,变成 的不交分割,再用前测度的可数可加性。
因此 Carathéodory 可测集族包含 Σ ⊗ T \Sigma\otimes\mathcal T Σ ⊗ T ;把 π ∗ \pi^* π ∗ 限制在后者,就得到乘积测度。σ \sigma σ -有限性还保证满足矩形规则的测度在这个 σ \sigma σ -代数上唯一。它来自测度扩张的唯一性定理:在有限测度块上由生成代数确定测度,再用可数个有限块穷尽整个空间。
图中的网格说明有限细分的相容性。把定义延伸到所有乘积可测集,还要经过上面的外测度和可数运算步骤,不能只凭一张有限网格图断言完成。
固定一个坐标,得到一个切片
对 E ⊆ X × Y E\subseteq X\times Y E ⊆ X × Y ,定义
E x = { y : ( x , y ) ∈ E } , E y = { x : ( x , y ) ∈ E } . E_x=\{y:(x,y)\in E\},\qquad
E^y=\{x:(x,y)\in E\}. E x = { y : ( x , y ) ∈ E } , E
若 E ∈ Σ ⊗ T E\in\Sigma\otimes\mathcal T E ∈ Σ ⊗ T ,这些切片都可测。证明不需要把 E E E 画出来:矩形的切片是 B B B 或空集;取补集、可数并与取切片可以交换,所以“所有切片可测”的集合族是包含矩形的 σ \sigma σ -代数,必然包含整个乘积 σ \sigma σ -代数。
测度还能由切片累计:
( μ × ν ) ( E ) = ∫ X ν ( E x ) d μ ( x ) = ∫ Y μ ( E y ) d ν ( y ) . (\mu\times\nu)(E)
=\int_X\nu(E_x)\,d\mu(x)
=\int_Y\mu(E^y)\,d\nu(y). ( μ × ν ) ( E ) = ∫ X ν ( E
这包括一个不宜忽略的结论:切片测度随外层坐标变化形成的函数是可测的。对矩形,ν ( E x ) = ν ( B ) 1 A ( x ) \nu(E_x)=\nu(B)\mathbf1_A(x) ν ( E x ) = ν ( B ) 1 A ( x ) ,公式容易验证;一般集合要把这个事实延伸出去。
展开切片公式从矩形推广的论证 先假设 μ ( X ) , ν ( Y ) \mu(X),\nu(Y) μ ( X ) , ν ( Y ) 都有限。把同时满足切片测度可测及上述两个积分等式的集合收成族 C \mathcal C C 。它包含矩形代数:有限不交并的切片仍不交,可以直接相加。
若 E n ↑ E E_n\uparrow E E n ↑ E ,每个切片也递增。测度从下连续和 MCT 说明等式传给极限。若 E n ↓ E E_n\downarrow E ,由于两个空间测度有限,可以用测度从上连续;切片测度函数被常数 或 控制,再用 DCT。因此 对递增并、递减交都封闭,称为单调类。
对图中的 D = { ( x , y ) : 0 ≤ y ≤ x ≤ 1 } D=\{(x,y):0\le y\le x\le1\} D = {( x , y ) : 0 ≤ y ≤ x ≤ 1 } ,固定 x x x 得 y ∈ [ 0 , x ] y\in[0,x] ,长度为 ;固定 得 ,长度为 。所以 。换方向后,上下限也必须跟着换。
Tonelli:非负函数允许总量无穷
若 f : X × Y → [ 0 , ∞ ] f:X\times Y\to[0,\infty] f : X × Y → [ 0 , ∞ ] 对乘积 σ \sigma σ -代数可测,则
∫ X × Y f d ( μ × ν ) = ∫ X ( ∫ Y f ( x , y ) d ν ( y ) ) d μ ( x ) = ∫ Y ( ∫ X f ( x , y ) d μ ( x ) ) d ν ( y ) . \begin{aligned}
\int_{X\times Y}f\,d(\mu\times\nu)
&=\int_X\left(\int_Yf(x,y)\,d\nu(y)\right)d\mu(x)\\
&=\int_Y\left(\int_Xf(x,y)\,d\mu(x)\right)d\nu(y).
\end{aligned} ∫ X
三者都允许为 ∞ \infty ∞ 。证明顺着前面已经搭好的阶梯走:f = 1 E f=\mathbf1_E f = 1 E 时是切片公式;非负简单函数时是有限可加性;一般 f f f 取简单函数 s n ↑ f s_n\uparrow f s n ,对内层积分、外层积分和乘积积分分别用 MCT,便得到等式及内层积分函数的可测性。
例如对 [ 0 , 1 ] 2 [0,1]^2 [ 0 , 1 ] 2 上的 f ( x , y ) = 1 / ( x y ) f(x,y)=1/(xy) f ( x , y ) = 1/ ( x y ) ,在坐标轴上任取零值,它非负可测。两个方向的迭代积分都为 ∞ \infty ∞ ,Tonelli 仍成立。“可以换序”不等于“答案一定有限”。
Fubini:有符号函数检查绝对可积
若 f f f 为可测实值函数,且
∫ X × Y ∣ f ∣ d ( μ × ν ) < ∞ , \int_{X\times Y}|f|\,d(\mu\times\nu)<\infty, ∫ X × Y ∣ f ∣ d ( μ × ν ) < ∞ ,
那么几乎所有切片都是可积函数,两种迭代积分都存在且与乘积积分相等。这就是本章使用的 Fubini 定理。
为什么只保证“几乎所有”?对 ∣ f ∣ |f| ∣ f ∣ 使用 Tonelli,函数 x ↦ ∫ Y ∣ f ( x , y ) ∣ d ν ( y ) x\mapsto\int_Y|f(x,y)|\,d\nu(y) x ↦ ∫ Y ∣ f ( x , y ) ∣ d ν ( y ) 的积分有限,因此它只能在一个零测集上取无穷。去掉这些例外切片,其余内层积分都是有限数。在例外的外层坐标上把内层结果定义为零,不改变外层积分。
然后分别对 f + f^+ f + 、f − f^- f − 用 Tonelli。绝对可积性保证两部分的总积分都有限,所以可以相减,得到 Fubini 的等式。反过来,判断可积性时常直接算 ∣ f ∣ |f| ∣ f ∣ 的某一种迭代积分;由 Tonelli,它有限就意味着乘积空间上的绝对积分有限。
对于乘积测度的完备化,可测函数的切片可测性也一般只在几乎处处成立。平面上的 Lebesgue 积分通常采用这一完备版本。遇到这个版本时,不要额外声称每一条切线都一定没有例外。
例题:换序以前,把区域写对
计算
I = ∫ 0 1 ∫ 0 1 − x ( x + y ) d y d x . I=\int_0^1\int_0^{1-x}(x+y)\,dy\,dx. I = ∫ 0 1 ∫ 0 1 − x ( x +
区域是第一象限中 x + y ≤ 1 x+y\le1 x + y ≤ 1 的三角形。把函数在区域外补零,得到非负 Borel 可测函数,可用 Tonelli。它还被 1 1 1 控制在面积有限的区域上,所以也绝对可积。
固定 x x x ,内层积分为 x ( 1 − x ) + ( 1 − x ) 2 / 2 x(1-x)+(1-x)^2/2 x ( 1 − x ) + ( 1 − x ) 2 。把各项展开后, 项恰好抵消,得到 。
这里的对称性让计算很短。区域不对称时,更需要实际解不等式,不能只把 d x dx d x 与 d y dy d y 的位置调换。
实验:弯曲边界带来的平方根
取 D = { ( x , y ) : 0 ≤ x ≤ 1 , x 2 ≤ y ≤ x } D=\{(x,y):0\le x\le1,\ x^2\le y\le x\} D = {( x , y ) : 0 ≤ x ≤ 1 , x 2 ≤ y ≤ x 。竖切长度为 。若改为固定 ,由 与 得 ,横切长度便是 。因此
m 2 ( D ) = ∫ 0 1 ( x − x 2 ) d x = ∫ 0 1 ( y − y ) d y = 1 6 . m_2(D)=\int_0^1(x-x^2)\,dx
=\int_0^1(\sqrt y-y)\,dy=\frac16. m 2 ( D ) = ∫ 0 1 ( x −
拖动切片前,判断横切片在 y = 1 / 4 y=1/4 y = 1/4 还是 y = 1 / 2 y=1/2 y = 1/2 更长;再切换方向,用端点坐标核对自己的上下限。观察累积区域时留意:只累积到中途的横切和竖切,扫过的子区域通常不同,不能要求它们的部分面积相等。
横切长度 y − y \sqrt y-y y − y 在 y = 1 / 4 y=1/4 y = 1/4 达到最大值 1 / 4 1/4 1/4 ;它并不在图像高度的中点最宽。定量判断应来自边界方程,而不只是看图猜。
没有绝对可积性,次序可以改变答案
用一个离散例子就能看清问题。在 N × N \mathbb N\times\mathbb N N × N 上取计数测度,令
a i j = { 1 , j = i , − 1 , j = i + 1 , 0 , 其他 . a_{ij}=\begin{cases}
1,&j=i,\\
-1,&j=i+1,\\
0,&\text{其他}.
\end{cases} a ij = ⎩
每一行都有一正一负,因此先逐行求和再加起来,得到 0 0 0 。第一列只有一个 1 1 1 ,以后每列都有一正一负,因此先逐列求和再加起来,得到 1 1 1 。两种迭代和均有意义,却不相等。
为什么 Fubini 没有保护它?因为每行的绝对值和为 2 2 2 ,故 ∑ i , j ∣ a i j ∣ = ∞ \sum_{i,j}|a_{ij}|=\infty ∑ i , j ∣ a ij ∣ = ∞ 。正项、负项大小的总和都无穷,整个有符号乘积积分没有定义。
取前 N N N 行、前 M M M 列的有限矩形,正项有 min ( N , M ) \min(N,M) min ( N , M ) 个,负项有 min ( N , M − 1 ) \min(N,M-1) min ( N , M − 1 ) 个。因此
S N , M = { 1 , N ≥ M , 0 , N < M . S_{N,M}=\begin{cases}1,&N\ge M,\\0,&N<M.\end{cases} S N , M = { 1 , 0 ,
固定 N N N 让 M → ∞ M\to\infty M → ∞ ,结果为零;固定 M M M 让 N → ∞ N\to\infty N → ∞ ,结果为一。有限和仍可以任意换序,发生分歧的是把两个无穷极限先后取完的方式。
分别选择宽矩形、正方形和高矩形,预测哪个边缘留下了未抵消的 1 1 1 ,再扫描核对。切换到非负数组 a i j = 2 − i − j a_{ij}=2^{-i-j} a ij = 2 − i − j 时,总和为 1 1 1 ,两个方向都由 Tonelli 保证一致。别把一张有限矩阵扫描的先后,当成无限迭代和的定义。
练习
题 1。 写出区域 0 ≤ y ≤ x ≤ 1 0\le y\le x\le1 0 ≤ y ≤ x ≤ 1 两种方向的积分限,并求 ∬ D x d x d y \iint_D x\,dx\,dy ∬ D x d x d y 。
查看解答 竖切为 ∫ 0 1 ∫ 0 x x d y d x = ∫ 0 1 x 2 d x = 1 / 3 \int_0^1\int_0^x x\,dy\,dx=\int_0^1x^2\,dx=1/3 ∫ 0 1 ∫ 0 x x d y d x = 。横切为 。函数非负且区域可测,Tonelli 允许这两种累计方式。
题 2。 若已算出 ∫ X ∫ Y ∣ f ( x , y ) ∣ d ν d μ = 7 \int_X\int_Y|f(x,y)|\,d\nu\,d\mu=7 ∫ X ∫ Y ∣ f ( x , y ) ∣ d ν d μ = 7 ,还需要另算一遍乘积空间的绝对积分,才能用 Fubini 吗?说明假设。
查看解答 不需要重复计算。在两个空间 σ \sigma σ -有限、f f f 乘积可测的假设下,Tonelli 应用于 ∣ f ∣ |f| ∣ f ∣ ,直接说明乘积绝对积分也为 7 7 7 。因此 f f f 可积,Fubini 随即适用。必须先对绝对值算出的结果是有限数,而不只是某个有符号迭代积分有限。
3 对非负可测函数使用 Tonelli 时,哪种结果是允许的?
A. 两个迭代积分分别为0与1 B. 乘积积分与两个迭代积分都为∞ C. 只要图像连续就可不管可测空间条件 D. 任意有符号函数都可直接换序
题 4。 对本章的带符号数组,计算 S 3 , 5 S_{3,5} S 3 , 5 、S 5 , 3 S_{5,3} S 5 , 3 、S 4 , 4 S_{4,4} S 。这些有限和是否说明有限求和也依赖顺序?
查看解答 分别为 0 , 1 , 1 0,1,1 0 , 1 , 1 。它们对应不同形状的有限矩形,所包含的项不同。在任意固定的有限矩形中,逐行与逐列相加都得到同一个值。无限迭代和不相等,是因为两个方向趋于无穷的先后不同,而非有限加法交换律失效。