可卷集
假设 F1 和 F2 是 Rn 中的两个闭集, 如果对任意的 R>0, 存在 R′>0 (可能依赖于 R) , 使得对任意的 x1∈F1, x2∈F2, 我们有∣x1+x2∣⩽R ⇒ ∣x1∣<R′,∣x2∣<R′,那么, 我们称 F1 与 F2 是可卷的.
更一般地, 假设 {Fi}i∈I 是 Rn 中的一族闭集, 如果对任意的 R>0, 存在 R′>0 (可能依赖于 R) , 使得对任意的可数子集 J⊂I, 对任意的 xj∈Fj, 其中 j∈J, 我们有∣j∈J∑xj∣⩽R ⇒ ∣xj∣<R′ 对所有j成立,那么, 我们称闭集族 {Fi}i∈I 是可卷的.
假设闭集 F1 与 F2 是可卷的, 那么, F1+F2 是闭集.
证明. 令
F=F2+F2 是闭集, 考虑
F 中的 Cauchy 列
{xi+yi}i⩾1⊂F, 其中,
xi∈F1,yi∈F2. 由于
{xi+yi}i⩾1 是 Cauchy 列, 所以是有界的, 即存在
R>0, 使得
∣xi+yi∣⩽R, 其中
i⩾1. 根据可卷集的性质, 存在
R′, 使得对每个
i,
∣xi∣⩽R′,
∣yi∣⩽R′. 从而, 存在
{xi}i⩾1 的子列
{xip}p⩾1 和
{yi}i⩾1 的子列
{yip}p⩾1, 它们都是 Cauchy 列, 从而, 存在
x∈F1,
y∈F2 (这是两个闭集) , 使得
p→∞limxip=x, p→∞limyip=y.所以,
x+y∈F 并且
i→∞limxi+yi=x+y.这说明
{xi+yi}i⩾1 的极限点仍然在
F 中, 即
F 是闭集.
如下的集合是可卷的:
1) | K1,⋯,Km⊂Rn 是紧集, K0 是闭集, 那么, {Ki}0⩽i⩽m 是可卷的. |
2) | 对 i=1,2,⋯,m, 我们令Fi=[xi,+∞).那么, {Fi}1⩽i⩽m 在 R 上是可卷的. |
3) | 我们在时空 R1+3=R1×R3 上考虑如下两个集合: |
• | 实心的未来光锥: C+={(t,x)∈R1+3∣∣t⩾∣x∣}. |
• | 平面 t=T 的未来: T∈R, 我们定义Rt⩾T1+3={(t,x)∈R1+3∣∣t⩾T}. |
那么, C+ 和 Rt⩾T1+3 是可卷的.
对于 u,v∈D′(Rn), 如果 supp(u) 和 supp(v) 是可卷的, 那么, 我们可以定义u∗v∈D′(Rn).实际上, 任意选取一组 χk∈C0∞(Rn), 使得对任意的 ∣x∣⩽k, χk(x)=1, 其中, k⩾1. 那么, 对任意的试验函数 φ∈D(Rn), 如下的极限定义了 u∗v: ⟨u∗v,φ⟩=k→∞lim⟨(χku)∗(χkv),φ⟩.
证明. 我们首先证明, 当给定了 {χk∈C0∞(Rn)}k⩾1 之后, 其中对任意的 ∣x∣⩽k, χk(x)=1, 极限k→∞lim⟨(χku)∗(χkv),φ⟩.存在. 实际上, 我们证明存在 N, 使得当 k,j⩾N 时, ⟨(χku)∗(χkv),φ⟩=⟨(χju)∗(χjv),φ⟩.考虑差⟨(χku)∗(χkv),φ⟩−⟨(χju)∗(χjv),φ⟩=I1⟨(χk−χj)u)∗(χkv),φ⟩+I2⟨(χju)∗(χk−χj)v),φ⟩.先考虑 I1, 我们假设它不是 0, 此时, supp(φ)∩{(supp(χk−χj)u)+supp(χkv)}=∅.特别地, supp(φ)∩{(supp(u)+supp(v)}=∅.假设 supp(φ) 在半径不超过 R 的球之内, 那么, 根据可卷集的性质, 存在 R′, 对于任意的 x+y∈supp(φ)∩{(supp(u)+supp(v)}, 其中, x∈supp(u), y∈supp(v), 我们有∣x∣⩽R′, ∣y∣⩽R′.我们取 N 足够大, 使得 χN∣∣BR′≡1. 那么, 对于 k,j⩾N, supp(φ)∩{(supp(χk−χj)u)+supp(χkv)}=∅与 ∣x∣⩽R′ 矛盾. 从而, I1≡0. 类似的, I2≡0.
上面的论证还表明⟨u∗v,φ⟩=k→∞lim⟨(χku)∗(χkv),φ⟩.不依赖于 {χk}k⩾1 的选取, 实际上, 我们只要令⟨u∗v,φ⟩=R→∞lim⟨(χRu)∗(χRv),φ⟩即可, 其中 χR∈C0∞(Rn) 在 BR 上恒为 1.
我们还要说明
u∗v 是分布. 实际上, 假设
supp(u)⊂BR, 我们有
⟨u∗v,φ⟩=k→∞lim⟨u,χk⋅((χkv)ˇ∗φ)⟩=⟨u,χN⋅((χNv)ˇ∗φ)⟩,其中,
N 如前述. 我们把分布定义所要证明的不等式留作作业.
对于 u,v,w∈D′(Rn), 假设 supp(u), supp(v) 和 supp(w) 是可卷的, 那么, 我们
1) | supp(u∗v)⊂supp(u)+supp(v); |
2) | u∗v=v∗u; |
3) | (u∗v)∗w=u∗(v∗w); |
4) | 对任意的多重指标 α, 我们有∂α(u∗v)=(∂αu)∗v=u∗(∂αv); |
证明. 我们只证明第四条, 其余的留作习题. 根据归纳法, 我们不妨假设
∣α∣=1. 首先观察到,
supp(∂αu) 和
supp(v) 还是可卷的. 我们选取试验函数
φ. 根据卷积的构造, 我们知道存在比较大的
k, 使得
⟨∂α(u∗v),φ⟩=(−1)∣α∣⟨u∗v,∂αφ⟩=(−1)∣α∣⟨χku,(χkv)ˇ∗∂αφ⟩=(−1)∣α∣⟨χku,∂α((χkv)ˇ∗φ)⟩=⟨∂α(χku),(χkv)ˇ∗φ⟩=⟨∂α(χk)u,(χkv)ˇ∗φ⟩+=⟨∂α(u)∗v,φ⟩⟨χk⋅∂α(u),(χkv)ˇ∗φ⟩=⟨(∂α(χk)u)∗(χkv),φ⟩+⟨∂α(u)∗v,φ⟩.上面第一项由于
∂α(χk)u 的支集在
∣x∣>k 上, 重复定理中的证明 (利用可卷性) , 我们知道这一项是
0, 所以,
⟨∂α(u∗v),φ⟩=⟨∂α(u)∗v,φ⟩.另外一个等号用交换律即可.
R3 中的波动方程
波动算子 □=−∂t2+△ 定义在 R1+3=R×R3 上的算子, 其中, 第一个坐标是时间 t 的坐标. 它作用在以 (t,x)∈R×Rn 为变量函数上的.
我们在
R1+3 中定义未来光锥
C+:
C+={(t,x)∈R1+3∣∣t=∣x∣}.除去
(0,0) 点外,
C+ 是
R4 中的一个光滑超曲面 (请证明这一点! ) , 我们用
dσ 表示
C+ 的曲面测度. 那么, 对于任意的
φ∈D(R1+3),
dσ 给出了一个
0 阶的分布:
⟨dσ,φ⟩=∫C+φ(t,x)∣∣C+dσ.由于
C+ 可以看作是函数图像 (
t=∣x∣) , 所以,
⟨dσ,φ⟩=2∫R3φ(∣x∣,x)dx. 分布 W=−4πt2+∣x∣2dσ 是 □ 的一个基本解, 即□(−4πt2+∣x∣2dσ)=δ0.
证明. 按照定义, 对任意的试验函数 φ∈D(R1+3), 我们有⟨□(t2+∣x∣2dσ),φ⟩=∫R3∣x∣(□φ)(∣x∣,x)dx.我们注意到, 右边是局部可积的.
根据上学期的作业 11, 在球坐标系下, 我们有
□=−∂t2+△=−∂t2+∂r2∂2+r2∂r∂+r21△S2,其中
r=∣x∣,
ϑ=∣x∣x∈S2 在单位球面上. 所以, 利用球坐标系, 我们有
∫R3∣x∣(□φ)(∣x∣,x)dx=∫R3∣x∣(□φ)(∣x∣,x)r2drdσS2=∫0∞∫S2(−∂t2φ+∂r2φ+2∂rφ+r1△S2φ)(r,r,ϑ)rdσSdr.根据上学期的作业 11 的结论 (球面上的散度公式) , 对于固定的
r>0, 我们有
∫S2(△S2φ)(r,r,ϑ)dσS=0.所以,
∫R3∣x∣(□φ)(∣x∣,x)dx=∫0∞∫S2(−∂t2φ+∂r2φ+2∂rφ)(r,r,ϑ)rdσS2dr=∫0∞∫S2(−∂t2φ+∂r2φ+2∂rφ)(r,r,ϑ)rdσS2dr=∫0∞∫S2−((∂t+∂r)(∂t−∂r)(rφ))(r,r,ϑ)dσS2dr.如果我们令
ψ=rφ, 那么,
⟨□(−t2+∣x∣2dσ),φ⟩=∫0∞∫S2((∂t+∂r)(∂t−∂r)(ψ))(r,r,ϑ)dσS2(ϑ)dr.我们令
L=∂t+∂r,
L=∂t−∂r. 对任意的光滑函数,
ϑ 固定, 我们计算
drd(f(r,r,ϑ))=(∂tf)(r,r,ϑ)+(∂rf)(r,r,ϑ).所以,
⟨□(−t2+∣x∣2dσ),φ⟩=∫0∞∫S2drd[(L(ψ))(r,r,ϑ)]dσS2(ϑ)dr=∫0∞drd[∫S2(L(ψ))(r,r,ϑ)dσS2(ϑ)]dr=ε→0lim∫S2(Lψ)(ε,ε,ϑ)dσS2(ϑ)另外,
(Lψ)(t,x)=φ(t,x)+∣x∣L(φ)(t,x).所以,
ε→0lim∫S2(Lψ)(ε,ε,ϑ)dσS2(ϑ)=∫S2φ(0,0)dσS2(ϑ)=4πφ(0,0).这就证明了结论.
给定
f∈D′(R1+3), 如果存在
T∈R, 使得
supp(f)⊂Rt⩾T1+3, 我们就说
f 的
过去是零.
假设 f∈D′(R1+3) 的过去是零, 那么, 存在唯一的过去为零的 u∈D′(R1+3), 使得□u=f.特别地, u 可以表示为u=−4πt2+∣x∣2dσ∗f.
证明. 我们注意到 Rt⩾T1+3 与 C+ 是两个可卷的闭集, 所以, u=−4πt2+∣x∣2dσ∗f是良好定义的 (因为 supp(dσ)⊂C+) . 由于Rt⩾T1+3+supp(dσ)⊂Rt⩾T1+3,所以, u 的过去为零.
现在证明唯一性: 假设
v∈D′(R1+3) 的过去为零, 并且
□v=f.那么, 我们有
v=v∗δ0=v∗□(−4πt2+∣x∣2dσ)=□(v)∗(−4πt2+∣x∣2dσ)=f∗(−4πt2+∣x∣2dσ)=u.这就证明了唯一性, 其中, 为了使得上面每一个式子都有定义, 我们用到了
v 的过去为零这个条件.
分布 −4πt2+∣x∣2dσ 是 □ 的唯一一个过去为零的基本解.
证明. 这是因为这样的基本解都满足
□u=δ0,其中,
δ0 的过去为零.
我们现在研究所谓的 Cauchy (初值) 问题, 这和常微分方程类似, 我们要在 t=0 这个时刻给定 u 的初始值 (此时, 由于方程对时间是 2 阶的, 我们需要给定 u(0,x) 和 (∂tu)(0,x)) , 然后解方程. 我们假设 u 是 R1+3 上的光滑函数, 它满足如下的波动方程: ⎩⎨⎧□u(t,x)=0, (t,x)∈R1+3;u∣∣t=0=u0(x), x∈R3;∂tu∣∣t=0=u1(x), x∈R3.我们要在 t⩾0 来研究这个问题: 给定了 t=0 的初始值, 我们想知道 u(t,x) 在未来 t⩾0 处的演化. 考虑过去为零的分布H(t)u(t,x)∈D′(R1+3).那么, 利用 Heaviside 函数 H(t) 的导数的计算, 我们有□(H(t)u(t,x))=D′−∂t(u(t,x)δ0(t)+H(t)∂tu)+H(t)△u=−∂t(u0(x)δ0(t)+H(t)∂tu)+H(t)△u=−∂t(u0(x)δ0(t))−u1(x)δ0(t)+H(t)□u.按照定义, u1(x)δ0(t) 是如下的分布: ⟨u1(x)δ0(t),φ(t,x)⟩=∫R3φ(0,x)u1(x)dx.这实际上是 t=0 所定义的曲面测度乘以 u1(x) 所定义的分布.
从而, □(H(t)u(t,x))=−∂t(u0(x)δ0(t))−u1(x)δ0(t).上式的右边是一个过去为零的分布, 所以, 对于 t⩾0, 我们有u(t,x)=−∂t(W∗[u0(x)δ0(t)])−W∗[u1(x)δ0(t)].
以上的公式是在假设波动方程有光滑解的情况下所给出的解的表达式! 不难看出, 我们只要假设 u(t,x) 是 C2 的, 上面的计算就成立.
下面我们把上面解的表达式显式地用微积分写清楚. 我们用
v 表示一个支集在
t=0 上的分布, 在应用的时候, 我们将会选取
v=−u0(x)δ0(t) 或者 −u1(x)δ0(t).为了一下计算明了, 我们这里不妨假设
u0 和
u1 都有紧支集 (否则, 我们将对下面的
vε 加一个
x 方向的截断函数, 请参考本次作业) . 我们再令
vε=χε(t)ui(x), i=0,1.其中
χε 是对
δ 的一个逼近. 很容易验证, 当
ε→0, 我们有
vε⟶D′v.从而,
W∗vε⟶D′W∗v.现在来计算
W∗vε. 按照定义, 我们有
W∗vε(x)=−⟨4π∣t′∣2+∣x′∣2dσ,χε(t−t′)ui(x−x′)⟩=−∫R34π∣x∣χε(t−∣x′∣)ui(x−x′)dx=−∫0∞∫S24πχε(t−r′)ui(x−x′)r′dσS2(ϑ)dr′=−∫0∞χε(t−r′)(r′∫S24πui(x−x′)dσS2(ϑ′))dr′⟶−t∫S24πui(x−tϑ′)dσS2(ϑ′).其中, 我们假设了
t⩾0 并且
x′=r′ϑ′. 在最后一步中, 我们用到了
χε(t−r′)⟶D′δt(r′).最终, 我们得到
W∗v(x)=−tH(t)∫S2ui(x−tϑ′)4πdσS2(ϑ′).特别地, 我们得到了波动方程的解的表达式: 对于
t⩾0, 我们有
u(t,x)=−t∫S2u1(x−tϑ′)4πdσS2(ϑ′)−∂t[t∫S2u0(x−tϑ′)4πdσS2(ϑ′)].在
u(t,x) 的表达式中, 由于
tϑ′ 是长度为
t 的向量, 所以, 积分项只与
u0 和
u1 在以
x 为中心以
t 为半径的圆上的值有关系.
我们还有一个更几何一点的表述,
u 在
(t,x) 处的值只和初始值在顶点在
(t,x) 处的倒向的光锥和
t=0 相截得到的球面
{x′∣∣∣x′−x∣=t} 上的值有关.