# 解的存在性 —— Peano 定理
对于一阶 ODE 的 IVP 问题 I0:f∈C(U;Rd)
⎩⎪⎨⎪⎧x˙=f(t,x)x(s)=z
# Euler 折线法
定义:
对于 I0,可以采用折线逼近的方法求解 t∈[s,s+t0] 中的一个解,称为 Euler 折线法。具体操作如下:
- 取 n∈N 均分 [s,s+t0],记线度为 δ=t0/n,记 s+kδ=sk
- 在 [s0,s1] 上构造
ϕn(t)=z+(t−s0)f(s,z)
- 在 [sk−1,sk] 上递归地构造(k=1,...,n)
ϕn(t)=ϕn(sk−1)+(t−sk−1)f(sk−1,ϕn(sk−1))
- 综上,在 [s,s+t0] 上定义 Euler 折线函数 ϕn(t)
性质:ϕ(t,n) 可以写成一个阶梯函数的积分。定义 ρn:[s,s+t0]→R 为
ρn(t):=f(sk,ϕ(sk)),t∈[sk,sk+1]
则有 Euler 折线的积分形式
ϕn(t)=z+∫stρn(τ)dτ
# Peano 定理
定理:设 I 是紧区间,设 {fn∈C(R;Rd)} 满足一致有界和等度连续性质,则存在子列 nk→∞ 以及 f∈C(I),使得
fnk⇉f
这是 Arzela-Ascoli 收敛定理。
在叙述 Peano 定理之前,给出一个圆柱体的定义:
Cρ,η:=[s−ρ,s+ρ]×B(z,η)⊂U
记速度上界:
M:=max{∣f(t,x)∣:(t,x)∈Cρ,η}
定理:对于 I0,其在 [s−ϵ,s+ϵ] 上至少存在一个解,这里
ϵ:=min{ρ,η/M}
这是 Peano 定理。
Proof
- 构造 Euler 折线函数 ϕn
- 利用 Arzela-Ascoli 定理说明一致收敛 ϕn⇉ϕ
- 验证极限函数 ϕ 是一个解,固定 t∗∈[s,s+t0],考虑 t→t∗ 时的情况
Δ(t,t∗):=∣∣∣∣∣t−t∗ϕ(t)−ϕ(t∗)−f(t∗,ϕ(t∗))∣∣∣∣∣,t=t∗
Remark
在构造 ϕn 时,也可以取一个收敛于 z 的初值序列 {zk} 开始按 Euler 折线法来定义 ϕn,这个方法还是有效的。
# 解的唯一性 —— Picard 迭代
# 积分算子与压缩映射
定义:对于 I0,设其解为 ϕ,那么
ϕ′(t)=f(t,ϕ(t))
两侧积分后就有 I0 的 解的积分方程
ϕ(t)=z+∫stf(τ,ϕ(τ))dτ
引理:设 ϕ∈C(I;Rd),且其图像 Γϕ⊂U。那么 ϕ 是 I0 的解当且仅当其满足 I0 的解的积分方程。
Remark
这个结果将微分方程的求解转化成为积分方程的求解,我们可以认为解的积分方程定义了一个算子,问题可以化为求解不动点。下面给出具体说明。
定义:对于 I0,定义关于方向场 f 的 积分算子 T,它将图像在 U 中的连续向量值函数 ϕ 映为另一个连续向量值函数 Tϕ
(Tϕ)(t):=z+∫stf(τ,ϕ(τ))dτ
度量空间、完备度量空间的概念参考本系列的其他文章。
定义:设 (X,d) 是一个度量空间。设 T:X→X 是一个映射。称 T 是一个 压缩映射,如果存在 c∈(0,1),使得
d(T(x),T(y))≤cd(x,y)
为了强调压缩程度,有时也称为 c - 压缩映射。
性质:压缩映射一定是连续映射。
定理:设 (X,d) 是一个完备度量空间,设 T:X→X 是一个 c - 压缩映射。则 T 有唯一的不动点 x∗。换言之,对任意 x0∈X,有
d(Tn(x0),x∗)≤1−ccnd(x0,T(x0))