形式陈述
设 U ⊂ R n 开,F : U → R n 为单射 C 1 映射,且 det D F ( x ) ≠ 0 对所有 x ∈ U 成立。由逆函数定理 理路 逆函数定理 Inverse function theorem 导数可逆的光滑映射在该点邻域内存在光滑局部逆。 ,V = F ( U ) 开,F : U → V 为 C 1 微分同胚。
令 S ⊂ U 是紧 Jordan 可测集:其边界可用有限个小盒子覆盖,使总体积任意小。令 f : F ( S ) → R 为有界 Riemann 可积 理路 Riemann 积分 Riemann integral 上下和或分割和在网格细化下共同收敛所定义的积分。 函数;在紧Jordan域上的积分可通过在外接矩形内补零来定义。则 F ( S ) 也是紧Jordan可测集,f ∘ F 在 S 上Riemann可积,并且
∫ F ( S ) f ( y ) d y = ∫ S f ( F ( x ) ) | det D F ( x ) | d x . 本页的 n 维Riemann积分沿用矩形分割、上下和;D F 由Jacobian矩阵 理路 Jacobian 矩阵 Jacobian matrix 多元映射各偏导数组成并表示其导数的矩阵。 表示。绝对值负责不带符号体积,不能省略。S 可以有洞或分片光滑边界,不要求它本身是一个开坐标域;光滑映射须在包含 S 的开集上定义。
更弱的版本只要求Jacobian在 S 上非零,或使用Lebesgue可积函数;这里证明的版本已经覆盖紧曲面片之间的换参数。含极点或接缝的常用坐标若不满足上述假设,须先去掉奇异集合或分片,再证明极限/零测集步骤,不能直接把公式当无条件代数规则。
直觉
换元做了两件事:f ( F ( x ) ) 把读数搬回参数坐标,| det D F ( x ) | 把每个小体积校正为真实大小。Jacobian矩阵一般会同时伸长、压缩、剪切、反射;它的行列式恰好汇总总体积倍率。映射若把同一区域走两遍,局部倍率不会替你删除多出来的一遍,所以单射或重数处理必不可少。
第一步:线性映射的体积倍率
可逆线性映射 A 把体积乘 | det A | 。可把 A 分解成初等矩阵:交换坐标不改体积,缩放某坐标乘缩放的绝对值,剪切只平移每张平行截面而不改其体积。对矩形先用Fubini定理 理路 Fubini 定理 Fubini's theorem 在适当可积条件下,多重积分等于任意次序的迭代积分。 逐截面积分,再用有限盒子逼近推广到Jordan集。这与行列式 理路 行列式 Determinant 交换含幺环上方阵的交替多线性标量不变量。 对上述三种初等操作的变化完全一致。
以下用最大范数 ‖ x ‖ ∞ 。若 H 在有界集 E 上是 L -Lipschitz,且 E , H ( E ) 都Jordan可测,则
vol ( H ( E ) ) ≤ L n vol ( E ) . 证明是用网格小立方体覆盖 E :边长 h 的每个交片在像中各坐标振幅至多 L h ,故可装入体积至多 ( L h ) n 的盒子。覆盖体积取极限即得不等式。同样的覆盖论证说明Lipschitz映射把零体积集送到零体积集。
第二步:局部近线性给出双边体积估计
在足够小、闭包位于 U 的立方体 Q 中,取 x 0 ∈ Q 、A = D F ( x 0 ) 。定义
H ( x ) = x 0 + A − 1 ( F ( x ) − F ( x 0 ) ) . 导数连续使 ‖ D H − I ‖ ∞ ≤ δ < 1 。沿 Q 内线段积分,得到
( 1 − δ ) ‖ x − y ‖ ∞ ≤ ‖ H ( x ) − H ( y ) ‖ ∞ ≤ ( 1 + δ ) ‖ x − y ‖ ∞ . 因此 H 与其像上的逆分别是 ( 1 + δ ) 与 1 / ( 1 − δ ) -Lipschitz。对任意Jordan可测 P ⊂ Q ,应用上一不等式两次,再乘线性倍率,可得
(1) | det A | ( 1 − δ ) n vol ( P ) ≤ vol ( F ( P ) ) ≤ | det A | ( 1 + δ ) n vol ( P ) . 这里所需的像可测性没有被暗中假定:紧Jordan集的边界为零体积,F 在紧集附近可由有限个有界导数的凸邻域控制,从而局部Lipschitz;同胚将边界送到像的边界,故后者仍为零体积。逆映射同理。
紧性保证对所有相交网格立方体可统一选足够细的网格,使 δ 任意小;D F 及其逆在相应紧邻域有界。式(1)就是把“微分近似”升级为“总体积误差可控”的关键。
第三步:从体积到一般积分
先设 f 连续。把 S 与细网格相交得到有限Jordan片 P j ,选 x j ∈ P j 。共有边界及其像均零体积,所以片积分可相加。F 与 f 在紧集上一致连续,故
∫ F ( S ) f = lim ∑ j f ( F ( x j ) ) vol ( F ( P j ) ) . 用式(1)替换每片体积,误差至多为某个趋零因子乘 sup | f | sup | det D F | vol ( S ) 。得到
lim ∑ j f ( F ( x j ) ) | det D F ( x j ) | vol ( P j ) , 它正是右侧积分。尤其取 f = 1 ,已证明每个紧Jordan子集 A ⊂ U 的体积公式 vol ( F ( A ) ) = ∫ A | det D F | 。
现在令 f 仅有界Riemann可积。在目标域外接矩形上补零,用Darboux上下阶梯函数 l ≤ f ≤ u 夹逼,使 ∫ F ( S ) ( u − l ) 任意小。阶梯函数的分界面是有限个坐标超平面;它们与 F ( S ) 的交及逆像边界都零体积,因为 F − 1 在紧邻域局部Lipschitz。因此各阶梯片的逆像Jordan可测。对这些片应用刚证明的体积公式,得到
∫ S l ( F ) | det D F | = ∫ F ( S ) l , ∫ S u ( F ) | det D F | = ∫ F ( S ) u . 正权重 | det D F | 保持夹逼。两边间隙任意小,便证明加权复合函数可积及换元等式。由于 | det D F | 在 S 上连续且有正下界,除以它又得 f ∘ F 可积。证明至此不只适用于连续被积函数。
例子与边界
极坐标环域:接缝怎么处理
F ( r , θ ) = ( r cos θ , r sin θ ) 有 det D F = r 。在 1 ≤ r ≤ 2 、ε ≤ θ ≤ 2 π − ε 的紧矩形上可选开邻域使 F 单射且 r > 0 ,所以面积为
∫ ε 2 π − ε ∫ 1 2 r d r d θ = 3 π − 3 ε . 让 ε ↓ 0 ,遗漏扇形在 0 < ε < π / 2 时位于 0 ≤ x ≤ 2 , | y | ≤ 2 sin ε 的窄矩形内,面积至多 8 sin ε → 0 ,得到完整环域面积 3 π 。也可分成几个不重叠扇区,仅让共有射线重叠;射线二维面积零。r = 0 的Jacobian为零,整圈角度也不单射,不能宣称闭极坐标矩形本身就是微分同胚。
反向与重复覆盖
反射 F ( x , y ) = ( − x , y ) 的行列式为 − 1 ,却保持正面积。换元若使用裸行列式,会把面积算成负数;有向微分形式积分才同时跟踪这种符号。
在极坐标中让 0 ≤ θ ≤ 4 π ,计算得到环域面积的两倍,这是两次覆盖。若映射可分成若干单射片,逐片换元后目标积分会带覆盖次数。不能在仍未去重的参数域上继续使用只算一次像集的公式。
更明显的一维例子是 x ↦ x 2 在 [ − 1 , 1 ] 上非单射:∫ − 1 1 2 | x | d x = 2 ,像集 [ 0 , 1 ] 长度为一。分成左右两支便看见重数二;不用绝对值则得到零,又变成了有向抵消,仍不是像集长度。
推论与应用
曲面面积 理路 曲面面积与通量积分 Surface area element · Scalar surface integral · Flux integral · 曲面积分 从切平行四边形导出面积密度,证明重参数不变性,并区分标量面积积分与依赖定向的通量。 在两张参数图之间换坐标时,叉积模恰好多出 | det D ϕ | ,本定理保证两张图的积分一致。长方形、带洞区域和分片光滑边界都可处理,只需保持局部可逆、重叠零测与边界可测条件。
概率密度换元还要指定原分布确有相应密度;一般推前测度 理路 推前测度 Pushforward measure · Image measure · 像测度 沿可测映射把源测度搬到目标空间,使目标集合的质量等于其原像的质量。 的积分恒等式对可测映射成立,却不自动给体积Jacobian。前者是换了测度,后者是把新测度相对Euclidean体积的倍率算出来,两项接口不可混同。
参考资料
Jiří Lebl,Basic Analysis II, §10.7 ,Proposition 10.7.1、Theorem 10.7.2及其证明:紧Jordan域、单射 C 1 映射与绝对Jacobian。该来源只需Jacobian在紧域上非零;本页采用在整个坐标开域非零的较强版本。
Michael Spivak,Calculus on Manifolds ,W. A. Benjamin,1965,Ch. 3 “Integration”:Riemann换元、零测边界与局部到整体分割。