外观
Lesson 52 Fubini 定理
约 1715 字大约 6 分钟
2025-5-9
期末考试是 6.14,所以 5.17 要补一次很长时间的课.
Fubini 定理
/Theorem/
设 f 在 [a,b]×[c,d] 上可积,且
F(x)=∫cdf(x,y)dy
有定义 (∀x∈[a,b]),则 F 在 [a,b] 上可积,且
∬[a,b]×[c,d]fdxdy=∫abF(x)dx
/Proof/
我们需要操作一下 Darboux 上下和. 对于 [a,b] 剖分 P1:[a,b]=⋃iSi,估计 U(P1,F),L(P1,F),下面以上和为例.
U(P1,F)=∑∣Si∣x∈SisupF(x)
再对 [c,d] 取一个剖分 P2:[c,d]=⋃jTj,则由条件有
F(x)≤U(P2,f(x,y))
所以:
U(P1,F)=i∑∣Si∣x∈SisupF(x)≤i∑∣Si∣x∈Sisup(j∑∣Tj∣y∈Tjsupf(x,y))≤i∑∣Si∣(j∑∣Tj∣Si×Tsupf(x,y))≤U(P1×P2,f)
同理,L(P1,F)≥L(P1×P2,f).
从而,我们得到一个双侧的估计. 再利用 f(x,y) 可积这一条件,可证:
∬[a,b]×[c,d]f(x,y)dxdy=∫abF(x)dx
进一步写出:
∬If(x,y)dxdy=∫ab(∫cdf(x,y)dy)dx
为方便引用,称这个式子为“累次积分”. 由于证明过程中的对称性,当然也可以写成对称形式:
∬If(x,y)dxdy=∫cd(∫abf(x,y)dx)dy
更简便的写法是:
∬If(x,y)dxdy=∫abdx∫cdf(x,y)dy
有人说,这个定理固然很好,但是不是每一次积分都在一个矩体上. 因此我们对这个定理做一次推广:
/Theorem/
对一般的 D,设 D⊆R2 有界,且 f 在 D 上可积,则有
∬Df=∬IfχD=Fubini∫abdx∫cdf(x,y)χD(x,y)dy
这几乎是显然的. 但是实际上很多地方 χD 为 0,所以很多时候没有必要全部算在积分之中. 设 ∀x∈[a,b],{x}×[c,d] 与 D 的交线是有限个线段,则
∫cdf(x,y)χD(x,y)dy=j∑∫sjtjf(x,y)dy
但是写求和号会使得记号变得臃肿,我们简记为
=∫{x}×R∩Df(x,y)dy
当然这并不是一个 Riemann 积分的标准写法,但是一般不至于引起混淆.
有下面的几个特例:
(1) D 为 y=ϕ1(x),y=ϕ2(x) 和 x=a,b 围成的曲边四边形.
/Theorem/ (Fubini 定理的特例 1)
区域:
D={(x,y)a≤x≤bϕ1(x)≤y≤ϕ2(x)}
则积分为
∬Df=∫abdx∫ϕ1(x)ϕ2(x)f(x,y)dy
(2) D 为 x=ψ1(y),x=ψ2(y) 和 y=c,d 围成的曲边四边形.
/Theorem/
区域:
D={(x,y)c≤y≤dψ1(y)≤x≤ψ2(y)}
则积分为
∬Df=∫cddy∫ψ1(y)ψ2(y)f(x,y)dx
Fubini 定理保证了换序积分可以成立,同时对高维情况这个结论也是成立的.
以 3 维情况为例:
先积分 1 维再积分 2 维;
∭Ωfdxdydz=∬Ωxydxdy∫({(x,y)}×R)∩Ωf(x,y,z)dz
先积分 2 维再积分 1 维.
∭Ωfdxdydz=∫z1z2dz∬({z}×R2)∩Ωf(x,y,z)dxdy
对于 n 维情况,先积分 1 维再积分 n−1 维,等价于先积分 n−1 维再积分 1 维.
/Example/
Ω=⎩⎨⎧(x,y,z)x2+y2≤1x2+z2≤1y2+z2≤1⎭⎬⎫
(三个圆柱之交)
求:Vol(Ω).
积分,
Vol(Ω)=∭Ωdxdydz=∬x2+y2≤1dxdy∫z2≤max{1−x2,1−y2}dz=∬x2+y2≤1dxdy⋅2max{1−x2,1−y2}
讨论 max{1−x2,1−y2} 取值,可以分为两个区域再求积分.
V1=∬Adxdy⋅21−x2=⋯
最后 Vol(Ω)=V1+V2.
/Example/ (Feynman 参数化 & Schwinger trick)
提示
这个等式由两人同时发现,他们因此争论是谁先发现的.
令区域如下:
V={(x0,⋯,xn)xi≥0,∀ix0+⋯+xn≤1}
在量子场论中常用积分:
Fn(a0,⋯,an)=∫⋯∫V(a1x1+⋯+anxn+a0(1−x1−⋯−xn))n−1dx1⋯dxn
只能使用 Fubini 定理,积分最后一个维度,得到
Fn(a0,⋯,an)=−n(an−a0)1(Fn−1(an;a0,⋯,an−1)−Fn−1(a0,⋯,an−1))
只需要从低维结果猜出:
Fn(a0,⋯,an)=n!a0a1⋯an1
换元公式
Rn 中两个向量 a,b,张成一个平行四边形,得到 S□=∣a×b∣.
多重积分可以理解为加权的 volume,所以我们考虑用线性变换,换一个基底之后重新计算这一体积.
回忆线性代数中的余弦定理:
cosθ=∣a∣⋅∣b∣⟨a,b⟩
因此平行四边形面积公式是
S□=∣a∣2∣b∣2−⟨a,b⟩2
换一组基底,应该仍然可以用 Gram 矩阵表示 S□,
S□=det(a1b1a2b2)
如果推广到三维,有
S□=(a12+a22+a32)(b12+b22+b32)−(a1b1+a2b2+a3b3)2=(a1b3−a3b1)2+(a2b3−a3b2)2+(a1b2−a2b1)2=∣(a2b3−a3b2,a3b1−a1b3,a1b2−a2b1)∣
这就定义了 a×b,为了记忆这一公式,写成行列式:
deti^a1b1j^a2b2k^a3b3=a×b
叉乘的方向?考虑 a⋅(a×b):显然为 0;同样,b⋅(a×b)=0,所以叉乘的方向垂直于 a,b 张成的平面. 但是还是有两种选择,我们选择右手系.
对于一个三维空间 V,要和 R3 建立联系,需要一种同构,这里的映射就是将 R3 中的基底放到 V 中. 叉乘是在 R3 中定义的,因此我们需要建立这样的一种关系:
/Claim/
若 e^x,e^y,e^z 构成右手系,则对于 a,b,c∈R3,有 Φ(a),Φ(b),Φ(c) 构成右手系 ⟺
deta1b1c1a2b2c2a3b3c3>0
(讲义上有证明)
/Theorem/
P 的 k 维 volume =detG.