这篇文章是 UncleBob 的 ODE (常微分方程) 中一阶 ODE 的基础核心理论部分笔记.
一阶ODE的基础核心理论
解的存在性 (Peano定理)
Euler折线法
以初值问题x˙=x;x(0)=1为例:
我们希望求得其解φ(t)在某个t0>0的取值, 我们去n很大, 记δ:=t0/n.
方向场在(0,1)放置的线段斜率为1, 在[0,δ]定义函数ϕn(t)=1+t, 在右端点处ϕn(δ)=1+δ.
方向场在(δ,1+δ)放置的线段斜率为1+δ, 在[δ,2δ]定义函数ϕn(t)=(1+δ)+(1+δ)(t−δ), 在右端点处ϕn(2δ)=(1+δ)2.
重复这个过程直到ϕn在整个[0,t0]定义好. 我们有
ϕn(t0)=ϕn(nδ)=(1+δ)n=(1+nt0)n
当n→+∞, 有ϕn(t0)→et0. 实际上在任何紧区间ϕn(t)⇉et.
这种逼近解的方式称作Euler折线法.
Peano定理
U⊂Rd+1开, f:U→Rd连续, (t0,x0)∈U, 则IVP{x˙=f(t,x)x(t0)=x0的解一定在t0的一个小领域存在.
定理 (Arzelà-Ascoli)
设{fn:n≥1}是紧区间I上一列在Rd中取值的向量值连续函数, 满足一致有界性质和等度连续性质, 则存在子列nk→∞以及I上的连续函数f, 使得fnk⇉f.
注
等度连续: ∀ϵ>0, ∃δ>0, s.t. ∣x−x′∣<δ时∀n, 有fn(x)−fn(x′)∣<ϵ.
定理 (Peano)
设f:U→Rd连续, 取ρ,η>0使得Cρ,η⊂U, 定义M:=max{∣f(t,x)∣:(t,x)∈Cρ,η}, 那么IVP{x˙=f(t,x)x(t0)=x0在区间[t0−ϵ,t0+ϵ]内至少有一个解, 这里ϵ=min{ρ,η/M}.
证明 这里只给出大纲, 具体过程略过
第一步: 利用Euler折线法构造一列逼近解ϕn(t);
第二步: 利用Arzelà-Ascoli定理找到目标解ϕ;
第三步: 证明ϕ是IVP在区间[t0,t0+ϵ]上的解.
解的唯一性(Picard迭代)
解的积分形式
设ϕ:I→Rd连续, 其图像Γϕ⊂U, 则ϕ为IVP{x˙=f(t,x)x(t0)=x0的解当且仅当其满足积分方程
ϕ(t)=x0+∫t0tf(τ,ϕ(τ))dτ
压缩映射原理
度量空间
定义
设X非空, d:X×X→[0,+∞)满足:
- 正定性: d(x,y)=0当且仅当x=y;
- 对称性: d(x,y)=d(y,x);
- 三角不等式: d(x,z)≤d(x,y)+d(y,z).
则称d时X上的一个度量.
定义
称{xn}⊂X为一个Cauchy列, 若d(xn,xm)→0,n,m→∞.
定义
称X为一个完备度量空间, 若{xn}⊂X是一个Cauchy列, 则∃xx∈X,xn→xx.
定义
(X,d),(Y,ρ)为两个度量空间, T:X→Y映射, 称T在x0∈X处连续, 若xn→x0⇒ρ(F(xn),F(x0))→0.
称T:X→Y连续, 若T在任何x∈X处连续.
定义
乘积度量:(X1,d1),(X2,d2)为两个度量空间, 令Y=X1×X2, 可以定义d:Y×Y→[0,+∞)为d((x1,x2),(y1,y2))=max(d1(x1,y1),d2(x2,y2)), 那么d是Y上的一个度量.
定义
开球 B(x0,r)={y∈X:d(x0,y)<r}
闭球 B(x0,r)={y∈X:d(x0,y)≤r}
性质
(X,d)是度量空间, (X×X,d^), 则
d:(X×X,d^)→(R,∣⋅∣)
是连续的.
定义
设V是一个实线性空间, 上面定义一个函数∣∣⋅∣∣:V→[0,+∞), 满足:
- 正定性: ∣∣v∣∣=0⇔v=0;
- 齐次性: ∣∣λv∣∣=∣λ∣∣∣v∣∣;
- 三角不等式: ∣∣u+v∣∣≤∣∣u∣∣+∣∣v∣∣.
则称∣∣⋅∣∣是V上的一个范数, (V,∣∣⋅∣∣)是一个赋范线性空间. (简写为n.l.s, 即normed linear space)
称∣∣v∣∣是向量v的长度或范数.
例
(Rn,∣⋅∣), ∣x∣=x12+⋯+xn2
∣∣x∣∣p:=(i=1∑n∣xi∣p)p1
(V,∣∣⋅∣∣)是n.l.s, 定义d:V×V→[0,∞)
d(v1,v2):=∣∣v1−v2∣∣
性质: d是V上的一个度量, 从而V在d下成为一个度量空间.
定义
称一个完备n.l.s.为Banach空间.
例
I:[a,b]是区间, V=C(I;Rn),
定义∣∣⋅∣∣∞:V→R+,
f↦∣∣f∣∣∞=t∈Isup∣f(t)∣=t∈Imax∣f(t)∣
性质: ∣∣⋅∣∣∞是C(I;Rn)上的一个范数, 称为无穷范数, 且(C(I;Rn),∣∣⋅∣∣∞)是Banach空间.
令K⊂Rn是一个紧集, 定义W=C(I;K)⊂V
性质: ∣∣⋅∣∣∞是C(I;K)上的一个范数, 称为无穷范数, 且(C(I;K),∣∣⋅∣∣∞)是Banach空间.
例
设V=C(I,R), I=[a,b],
∣∣f∣∣1=∫ab∣f(x)∣dx
∣∣⋅∣∣1是V上的范数, (C(I;R),∣∣⋅∣∣1)不是Banach空间.
例
设V=Mn(R), A∈Mn(R),
∣∣A∣∣2=i,j=1∑naij2
∣∣A∣∣op=∣v∣=1sup∣Av∣
则∣∣⋅∣∣op是Mn(R)上的一个范数, 且∃cn>1, s.t.
cn−1∣∣A∣∣op≤∣∣A∣∣2≤cn∣∣A∣∣op
压缩映射
定义
(X,d),T:X→X映射, 设0<c<1, 称T是一个c−压缩映射, 若d(T(x),T(y))≤c⋅d(x,y).
定理 (Banach不动点定理)
(X,d)是完备度量空间, T:X→X为c−压缩映射, c∈(0,1), 则T存在唯一的不动点x∗, T(x∗)=x∗. 进一步∀x0∈X,
d(Tn(x0),x∗)≤1−ccnd(x0,T(x0))
证明
x0∈X, 令xn=Tn(x0), {Tn(x0):n∈N}⊂X为Cauchy列,
d(Tn+1(x0),Tn(x0))≤d(T(Tn(x0)),T(Tn−1(x0)))≤c⋅d(Tn(x0),Tn−1(x0))≤⋯≤cnd(T(x0),x0)
d(Tn+k(x0),Tn(x0))≤d(Tn+k(x0),Tn+k−1(x0))+⋯+d(Tn+1(x0),Tn(x0))≤(cn+k−1+⋯+cn)d(T(x0),x0)≤1−ccnd(T(x0),x0)
由X完备, 存在唯一x∗∈X, s.t.
Tn(x0)→x∗,d(Tn(x0),x∗)→0
Lipschitz性质与解的唯一性
定义
U⊂Rd+1开, f:U→Rd连续, D⊂U, 称f∈Lipx(D), 若∃L=L(0)>0 s.t. ∀t,x,x′, (t,x),(t,x′)∈D, 有
∣f(t,x)−f(t,x′)∣≤L∣x−x′∣
定义
f:U→Rd连续, 称f∈Lipx,loc(U), 若∀K⊂U, 都存在L(K)>0, s.t.
f∈Lipx(K)
推论
U⊂Rd+1开, f:U→Rd为C1光滑 ⇒ f∈Lipx,loc(U).
定理(Picard-Lindelöf)
U⊂Rd+1开, f:U→Rd连续, 且f∈Lipx,loc(U). 设(t0,x0)∈U, 则IVP{x˙=f(t,x)x(t0)=x0
的解局部存在且唯一, 即存在δ>0以及φ:[t0−δ,t0+δ]→R s.t. φ为IVP的解且在[t0−δ,t0+δ]上是唯一解.
证明
记R=[t0−ϵ,t0+ϵ]×B(x0,ϵ),M=max(t,x)∈R∣f(t,x)∣,δ0=Mϵ,
第一步: 构造一个完备度量空间X:=C([t0−δ,t0+δ];B(x0,ϵ));
第二步: Picard映射T:X→C([t0−δ,t0+δ];Rd),ϕ↦Tϕ
(Tϕ)(t)=x0+∫t0tf(τ,ϕ(τ))dτ,t∈[t0−δ,t0+δ]
∣Tϕ(t)−x0∣=∣∫t0tf(τ,ϕ(τ))dτ∣≤∫t0t∣f(τ,ϕ(τ))∣dτ≤∣t−t0∣M≤δ0M≤ϵ
⇒Tϕ∈X;
第三步: T为压缩映射. ϕ,ψ∈X,
∣Tϕ(t)−Tψ(t)∣≤∫t0t∣f(τ,ϕ(τ))−f(τ,ψ(τ))∣dτ≤∫t0tL∣ϕ(τ)−ψ(τ)∣dτ≤∣∣ϕ−ψ∣∣∞,∀t∈[t0−δ,t0+δ]
取δ=min{δ0,2L1}
d(Tϕ,Tψ)=∣∣Tϕ−Tψ∣∣∞≤Lδ∣∣ϕ−ψ∣∣∞=21∣∣ϕ−ψ∣∣∞≤21d(ϕ,ψ)
即T为21-压缩映射.
由Banach, 存在唯一的φ∈X s.t. Tφ=φ, ∀t∈[t0−δ,to+δ],
φ(t)=Tφ(t)=∫t0tf(τ,ϕ(τ))dτ
在[t0−δ,t0+δ]上φ为IVP唯一解.
Picard迭代
引理 一个著名积分
n!(t1−t0)n=∫t0t1∫t0t2∫t0t3⋯∫t0tndtn+1dtn⋯dt2
取δ=δ0=Mϵ, T:X→X, ∀ϕ,ψ,t0≤t≤t0+δ0
∣Tnϕ(t)−Tnψ(t)∣=∣∫t0tf(t1,Tn−1ϕ(t1))−f(t1,Tn−1ψ(t1))dt1∣≤L∫t0t∣Tn−1ϕ(t1)−Tn−1ψ(t1)∣dt1≤L∫t0tL∫t0t1∣Tn−2ϕ(t2)−Tn−2ψ(t2)∣dt2dt1≤⋯≤Ln∫t0t∫t0t1⋯∫t0tn−1∣ϕ(tn)−ψ(tn)∣dtn⋯dt2dt1≤Ln∣∣ϕ−ψ∣∣∞n!δ0n
即
∣∣Tnϕ−Tnψ∣∣∞≤n!(Lδ0)n∣∣ϕ−ψ∣∣∞
断言: ∀ϕ∈X, {Tnϕ}⊂X为Cauchy列.
令ψ=Tkϕ, 有
∣∣Tn+kϕ−Tnϕ∣∣∞≤n!(Lδ0)n∣∣Tkϕ−ϕ∣∣∞≤2(∣x0∣+ϵ)n!(Lδ0)n→0,n→∞
∣∣ϕ∣∣∞≤∣x0∣+ϵ, 存在唯一的φ s.t. Tnϕ→φ,n→∞.
定理
I⊂R区间, f:I×Rd→Rd连续, 且∃L>0, s.t. ∀t∈I,∀x,x′∈Rd, 均有
∣f(t,x)−f(t,x′)∣≤L∣x−x′∣
那么∀t0∈I,∀x0∈Rd, IVP{x˙=f(t,x)x(t0)=x0的解在I上存在且唯一.
证明
任取J⊂I为紧区间, X=C(J;Rd), T:X→X,Tϕ(t)=x0+∫t0tf(τ,ϕ(τ))dτ, 证明同上, 得到单调上升Jn, I=∪Jn.
推论
I⊂R区间, A:I→Md(R), b:I→Rd连续, 则IVP{X˙=A(t)X+b(t)X(t0)=X0的解在I上存在且唯一.
证明
∀J⊂I为紧区间, A:J→Md(R), L=maxt∈J∣∣A(t)∣∣<∞, f(t,X):=A(t)X+b(t), 对于t∈J,X,Y∈Rd,
∣f(t,X)−f(t,Y)∣=∣A(t)(X−Y)∣≤∣∣A(t)∣∣∣X−Y∣≤L∣X−Y∣
极大解与极大积分路线
极大解的定义与存在性
定义
ω=(t0,x0), 称(ϕω,Iω)是I(ω)的一个极大解, 若以下两条成立:
- ϕω:Iω→Rd是I(ω)的解;
- 若(ψ,J)也是I(ω)的解, 且ψ为ϕω的延拓 (即Iω⊂J,ψ∣Iω=ϕω) , 则J=Iω,ψ=ϕω.
注:
若(ϕω,Iω)为极大解, 则Iω为开区间, 因为若Iω=(a,b], (b,ϕω(b))∈U, ϕω可以延拓到(a,b+ϵ)上, 从而与上述2) 矛盾.
从现在开始, 总假设∀ω=(t0,x0)∈U初值问题, I(ω)的解总是局部存在且唯一.
定理
U⊂Rd+1开, f:U→Rd连续且上述假设成立, 则∀(t0,ω0)∈U, I(ω)的极大解存在且唯一.
关键引理
(ϕ,I)与(ψ,J)均为I(ω)的解, 则ϕ∣I∩J=ψ∣I∩J.
证明
设I∩J=(a,b)∋t0, ϕ(t0)=ψ(t0)=x0.
反证: 不妨设∃t0<t1<b, s.t. ϕ(ti)=ψ(t1), ∃t0≤t∗<t1, s.t. ϕ∣[t0,t∗]=ψ∣[t0,t∗],ϕ(t∗)=ψ(t∗)=x∗, 但∀δ>0,ϕ∣(t∗,t∗+δ)=ψ∣(t∗,t∗+δ), 由假设, I(t∗,x∗)初值的解存在且唯一, 即∃δ0>0以及η:(t∗−δ0,t∗+δ0)→Rd为I(t∗,x∗)在(t∗−δ0,t∗+δ0)的唯一解, 矛盾.
称(ϕ,I)是I(ω)的一个局部解, 若I为开区间, 且ϕ:I→Rd为解.
记S={(ϕ,I):(ϕ,I)为I(ω)的局部解},
定义Imax=∪(ϕ,I)∈SI, 为开区间.
定义ϕmax:Imax→Rd,t↦ϕmax(t):=ϕt(t), ∃(ϕt,It), s.t. t∈It.
ϕmax′(t)=f(t,ϕmax(t)), ϕmax∣It=ϕt∣It, ϕmax′(t)=ϕt′(t)=f(t,ϕt(t)).
断言: (ϕmax,Imax)为I(ω)的一个极大解. 1) 成立; 对于2) , ψ:J→Rd为I(ω)的解, 且ψ是ϕmax的延拓⇒ J=Imax,ψ=ϕmax.
断言: J为开区间. 若已证J为开区间, (ψ,J)是一个局部解, J⊂Imax⇒Imax=J⇒ψ=ϕmax.
若ψ:J=(a∗,b∗]→Rd为解, ψ:J=(a∗,b∗+δ)→Rd为解, (a∗,b∗+δ)⊂Imax⊂J=(a∗,b∗], 矛盾. 从而成立.
唯一性: (ϕmax,Imax)为I(ω)的极大解, (ψmax,Jmax)也为极大解, 那么(ψmax,Jmax)是一个和局部解, 有Jmax⊂Imax, 同理有Imax⊂Jmax, 从而Jmax=Imax, 由引理知二者相等.
注:
- 假设(ϕmax,Imax)为I(ω)的极大解, 则称Imax为I(ω)的解的极大定义区间.
- 若没有假设, 用Zorn引理证明I(ω)的极大解必存在.
定义
若(ϕmax,Imax)是I(ω)的极大解, 称Γϕmax={(t,ϕmax(t)):t∈Imax}为极大积分曲线.
推论
过(t0,ω0)的极大积分曲线存在且唯一.
定义
称Γ是一条极大积分曲线, 若Γ是某个I(ω)的极大积分曲线.
极大解的几何性质
性质
设Γ1与Γ2是x˙=f(t,x)的两条极大积分曲线, 则要么Γ1∩Γ2=∅, 要么Γ1=Γ2.
引理
若Γ为I(t0,x0)的极大积分曲线, 且(t1,x1)∈Γ, 那么Γ也是I(t1,x1)的极大积分曲线.
证明
等价于证明I(t0,x0)的极大解(ϕmax,Imax)也是I(t1,x1)的极大解(ψmax,Jmax). 易证.
性质的证明
设Γ1是I(t0,x0)极大积分曲线, Γ2是I(t1,x1)极大积分曲线.
若Γ1∩Γ2=∅, 成立; 否则, ∃(t∗,x∗)∈Γ1∩Γ2, 从而Γ1也是I(t∗,x∗)极大积分曲线, Γ2也是I(t∗,x∗)极大积分曲线, 有Γ1=Γ2.
定理
U⊂Rd+1开, f:U→Rd连续, (t0,x0)∈U, 则过(t0,x0)的极大积分曲线“趋向于U的边界”, 即: 若(ϕmax,Imax)是I(t0,x0)的极大解, Imax=(a,b), 则∀K⊂D紧, 均存在a<t∗<t∗<b, s.t. 当t∈(a,t∗)∪(t∗,b)时, (t,ϕmax(t))∈/K.
证明
以t∗的存在性为例:
若b=+∞, 那么t∗存在是显然的;
若b<+∞, ∃t0≤t∗<b, s.t. 当t∈(t∗,b), (t,ϕmax(t))∈/K.
反证: ∃tn↑b, s.t. (tn,ϕ(tn))∈K紧. 不妨设ϕ(tn)→x∗, 取R=[b−2ϵ,b+2ϵ]×B(x∗,2ϵ)⊂U, 令M=max(t,x)∈R∣f(t,x)∣, 取tn, s.t. bMϵ<tn<b, 对(tn,ϕ(tn))使用Peano延拓定理: {x˙=f(t,x)x(tn)=ϕ(tn)解在[tn−Mϵ,tn+Mϵ]上存在, 设解为ψ(t), 定义
ψ(t)={ϕ(t),t≤tnψ(t),tn<t<tn+Mϵ
ψ(t)为解, 与ϕ为I(t0,x0)的极大解矛盾.
自治系统的极大解的性质
V⊂Rd为开集, g:V→Rd连续. 假设{x˙=g(x)x(t0)=x0的解总是局部存在且唯一.
f:R×V→Rd,(t,x)↦f(t,x):=g(x).
性质
x0∈V, (ϕmax,Imax)为I(0,x0)的极大解 ⇔ (ψmax,Jmax)为I(t0,x0)的极大解, 其中ψmax(t)=ϕmax(t−t0), Jmax=t0+Imax.
定义
{x˙=g(x)x(t0)=x0的极大解为(ϕmax,Imax), 定义I(t0,x0)的通过x0的极大相曲线为ϕmax的像, 即ϕmax(Imax).
推论
x˙=g(x), ∀x0∈V, 则
- 通过x0的极大相曲线唯一;
- 两条极大相曲线要么不相交, 要么完全重合.
性质
V⊂Rd, g:V→Rd, x˙=g(x), 任取K⊂V, K紧, ∀x0∈V, 设I(0,x0)的解为(ϕmax,Imax), Imax=(a,b), 则
- 要么b=+∞* (或者a=−∞) *;
- 要么b<+∞, 且存在0≤t∗<b, s.t. ϕ(t)∈/K, t∗<t<b* (或者a>−∞, 且存在a<t∗≤0, s.t. ϕ(t)∈/K, a<t<t∗) *.
证明
若b=+∞, 成立. 下设b<+∞.
f:Rd×V→Rd,(t,x)↦f(t,x)=g(x), 取K^=[0,b]×K⊂Rd+1紧, 由非自治系统结论知, ∃0≤t<b, s.t. (t,ϕ(t))∈/K^, 从而ϕ(t)∈/K.
解关于初值和参数的连续依赖性
对于{x˙=f(t,x)x(s)=z, 设解为ϕ(t,s,z).
Grönwall不等式
我们考虑对于
ϕ(t,s,z)=z+∫stf(τ,ϕ(τ,s,z))dτ
和
ϕ(t,s,z′)=z′+∫stf(τ,ϕ(τ,s,z′))dτ
定义
η(t):=∣ϕ(t,s,z)−ϕ(t,s,z′)∣≤∣z−z′∣+∫st∣f(τ,ϕ(τ,s,z))−f(τ,ϕ(τ,s,z′))∣dτ≤∣z−z′∣+L∫st∣ϕ(τ,s,z)−ϕ(τ,s,z′)∣dτ
即
0≤η(t)≤∣z−z′∣+L∫stη(τ)dτ
若η(t)=a+L∫stη(τ)dτ, 那么η′(t)=Lη(t), 有η(t)=aeL(t−s). 我们希望η(t)≤∣z−z′∣eL(t−s).
**定理 (Grönwall不等式) **
η:[t0,t1]→[0,∞)连续, 设a,b≥0, s.t.
η(t)≤a+b∫t0tη(τ)dτ,t∈[t0,t1]
则0≤η(t)≤aeb(t−t0).
证明
先设a>0, 令φ(t):=a+b∫t0tη(τ)dτ, 0≤η(t)≤φ(t), 那么φ∈C1, φ′(t)=bη(t), 从而
(lnφ(t))′=φ(t)φ′(t)=φ(t)bη(t)≤b
lnφ(t)−lnφ(t0)=∫t0t(lnϕ(τ))′dτ≤b(t−t0)
从而φ(t0)φ(t)≤eb(t−t0), 0≤η(t)≤φ(t)≤aeb(t−t0).
下设a=0. ∀a′>0,
0≤η(t)≤a′eb(t−t0)≤a′eb(t1−t0),t∈[t0,t1]
从而η(t)≡0.
解关于初值的连续依赖
引理
f:U→Rd连续, f∈Lipx,loc(U), 设(s,z)∈U, 设ϕ(t,s,z):[s−κ1,s+κ2]→Rd是I(s,z)在区间[s−κ1,s+κ2]上的解. 记Γ={(t,ϕ(t,s,z)):t∈I}, 令Γ(ϵ)={(t,x)∈U:d((t,x),Γ)≤ϵ}, 可知Γ(ϵ)⊂U为紧集. 设L=L(Γ(ϵ))>0, ∃δ=δ(L,κ1,κ2,ϵ)>0, s.t. ∀z′∈B(z,δ), I(s,z′)的解ϕ(t,s,z′)的定义区间至少也为I=[s−κ1,s+κ2], 且
∣ϕ(t,s,z)−ϕ(t,s,z′)∣≤∣z−z′∣eL∣t−s∣
证明
ϕ(t,s,z)=z+∫stf(τ,ϕ(τ,s,z))dτ
ϕ(t,s,z′)=z′+∫stf(τ,ϕ(τ,s,z′))dτ
η(t):=∣ϕ(t,s,z)−ϕ(t,s,z′)∣≤∣z−z′∣+∫st∣f(τ,ϕ(τ,s,z))−f(τ,ϕ(τ,s,z′))∣dτ≤∣z−z′∣+L∫st∣ϕ(τ,s,z)−ϕ(τ,s,z′)∣dτ≤∣z−z′∣+L∫stη(τ)dτ
∀z′∈B(z,δ), 设ϕ(t,s,z′)为I(s,z′)的解, 由极大积分曲线离开任何紧子集, 对Γ(s), 存在T>s, s.t. (t,ϕ(t,s,z′))∈Γ(ϵ), t∈[s,T]. ∃tn↓T, s.t.(tn,ϕ(tn,s,z))∈/Γ(ϵ).
断言: T≥s+κ2.
否则s<T<s+κ2, 则有d((T,ϕ(T,s,z′)),Γ)=ϵ, 对∀s≤t≤T, 有
η(t)≤∣z−z′∣+L∫stη(τ)dτ
由Grönwall
0≤η(t)≤∣z−z′∣eL(t−s)
那么
∣(T,ϕ(T,s,z))−(T,ϕ(T,s,z′))∣=∣ϕ(T,s,z)−ϕ(T,s,z′)∣=η(T)≤∣z−z′∣eL(T−s)≤δeLκ2=2ϵ
矛盾, 从而断言成立.
定理
对∀(s,z)∈U, 记I(s,z)的极大解为(ϕ(t,s,z),I(s,z)), 定义W={(t,s,z):(s,z)∈U,t∈I(s,z)}⊂Rd+2, 则W为开集, 且为ϕ的定义域, ϕ在W上连续.
证明
∀σ0=(t0,s0,z0)∈W, ∃δ>0,C>0, s.t. U(δ0):=(t0−δ,t0+δ)×(s0−δ,s0+δ)×B(z0,δ)⊂W, 且∀σ=(t,s,z)∈U(σ0), 有
∣ϕ(t,s,z)−ϕ(t0,s0,z0)∣≤C(∣t−t0∣+∣s−s0∣+∣z−z0∣)
记zs=ϕ(s,s0,z0), (s,zs)∈Γ(s,z0), 那么I(s,zs)的极大解ϕ(t,s,zs)=ϕ(t,s0,z0).
由t0∈I(s0,z0), ∃δ0>0, s.t. [t0−δ0,t0+δ0]⊂I(s0,z0), 记Γ={(t,ϕ(t,s0,z0)):t∈[s0−δ0,t0+δ0]}, M=max(t,x)∈Γ(ϵ)∣f(t,x)∣, 对(s,zs)与(s,z)用性质, 有∣z−zs∣<δ1=e−L(t0−s0+2δ0)ϵ/2, 故ϕ(t,s,z)的解也至少在[s0+δ0,t0+δ0]上存在且连续, 那么∣z−zs∣≤∣z−z0∣+∣z0−zs∣<δ+M∣s−s0∣<δ+Mδ<δ1. 当δ=δ1/(1+M)时∣s−s0∣<δ,∣z−z0∣<δ, 故ϕ(t,s,z)至少在[s0−δ0,t0+δ0]上有定义, (t0−δ,t0+δ)⊂[s0−δ0,t0+δ0]⊂I(s,t), (t0−δ,t0+δ)×(s0−δ,s0+δ)×B(z0,δ)⊂W, 从而W为开集.
∣ϕ(t,s,z)−ϕ(t0,s0,z0)∣≤∣ϕ(t,s,z)−ϕ(t,s0,z0)∣+∣ϕ(t,s0,z0)−ϕ(t0,s0,z0)∣≤∣ϕ(t,s,z)−ϕ(t,s,zs)∣+∣ϕ(t,s0,z0)−ϕ(t0,s0,z0)∣≤∣ϕ(t,s,z)−ϕ(t,s,zs)∣+∫t0t∣f(τ,s0,z0)∣dτ≤∣z−zs∣eLκ+M∣t−t0∣≤(∣z−z0∣+M∣s−s0∣)eLκ+M∣t−t0∣
解关于参数的连续依赖
U~⊂Rd+1+d′开, f:U~→Rd连续, 考虑IVP{x˙=f(t,x,λ)x(s)=z的解ϕ(t,s,z,λ).
定义
称f∈Lip(x,λ),loc(U~), 若∀K⊂U~为紧集, ∃L(K)>0, s.t.
∣f(t,x,λ)−f(t,x′,λ′)∣≤L(∣x−x′∣+∣λ−λ′∣)
性质
设U~⊂Rd+1+d′开, f:U~→Rd连续且f∈Lip(x,λ),loc(U~), ∀(s,z,λ)∈U~, 记I(s,z,λ)的极大解为(ϕ(t,s,z,λ),I(s,z,λ)), 则W~:={(t,s,z,λ):(s,z,λ)∈U~,t∈I(s,z,λ)}为Rd+1+d′中的开集, 其为ϕ的定义域, 进一步ϕ在W~上连续.
证明
参数变初值
记F:U~→Rd+d′,(t,x,y)↦(f(t,x,y),0), 那么F∈C(U~)∩Lip(x,y),loc(U~), 对于{X˙=F(t,X)X(s)=(z,λ), X=(x,y), 其解ϕ(t,s,z,λ), 那么ϕ(t,s,z,λ)在W~上连续, (x˙,y˙)=(f(t,x,y),0), {x˙=f(t,x,y)x(s)=zy˙=0y(s)=λ, 有y(t)≡λ, 带入上方即有结论正确.
解关于初值和参数的光滑依赖性
C1 光滑情形
f:U→Rd, f∈C1(U;Rd), 那么f∈C(U)∩Lipx,loc(U).
对于ϕ(t,s,z)=z+∫stf(τ,ϕ(τ,s,z))dτ,
∂z∂ϕ(t,s,z)=Id+∫st∂x∂f(τ,ϕ(τ,s,z))∂z∂ϕ(τ,s,z)dτ
记为Φ(t)=Id+∫stA(τ)Φ(τ)dτ, 则Φ′(t)=A(t)Φ(t).
考虑x˙=f(t,x)一阶变分方程{X˙=A(t)XX(s)=Id.
A(t,s,z):=∂x∂f(t,ϕ(t,s,z))连续, 下面讨论{X˙=A(t,s,z)X=F(t,X,s,z)X(s)=Id. 记ω=(s,z), 将s换为τ.
引理
I⊂R紧, K⊂Rm紧, A:I×K→Md(R),(t,ω)↦A(t,ω)连续, 取τ∈I, 取ω∈K, 取B∈Md(R), 则{X˙=A(t,ω)XX(τ)=B的解在I上有定义且唯一. 记解为Φ(t,τ,ω), 则Φ的定义域为I×I×K且Φ在I×I×K上连续.
证明
S=C(I×I×K;Md(R)), 定义P:S→S,ϕ(t,τ,ω)↦(Pϕ)(t,τ,ω), 其中
(Pϕ)(t,τ,ω):=B+∫τtA(u,ω)ϕ(u,τ,ω)du
可以证明{Pnϕ}是S中Cauchy列, 设Pnϕ→Φ, Φ(t,τ,ω)=B+∫τtA(u,ω)Φ(u,τ,ω)du.
性质
(t,ω)∈V⊂Rd+1, ∀ω∈Rd, 定义V的截口Vω={t∈R:(t,ω)∈V}, 设Vω均为开区间, A:V→MdR,(t,ω)↦A(t,ω), A连续. 固定B∈Md(R), ∀(τ,ω)∈V, 考虑初值问题{X˙=A(t,ω)XX(τ)=B, 则该初值问题的解在Vω上存在且唯一, 记该解为ϕ(t,τ,ω), 定义W={(t,τ,ω):(τ,ω)∈V,t∈Vm}, 则ϕ在W上连续.
证明
固定(τ,ω)∈V, 固定(τ0,ω0,t0), ω→Vω, 任取I, s.t. t0,τ0∈I⊂Vω紧, 取K=B(ω0,δ), 使用引理: 有解Φ(t,τ,ω)在I×I×B(ω0,δ)上连续.
V={(t,s,z):(s,z)∈U,t∈I(s,z)}, A:V→Md(R),(t,s,z)↦A(t,s,z)=∂x∂f(t,ϕ(t,s,z)), 初值问题{X˙=A(t,s,z)XX(τ)=Id的解存在且唯一, 记为Φ(t,τ,s,z), W={(t,τ,s,z):(s,z)∈U,τ,t∈I(s,z)}, Φ在W上连续.
初值问题{X˙=A(t,s,z)XX(s)=Id的解记为D(t,s,z)=Φ(t,τ,s,z). 定义域为V={(t,s,z):(s,z)∈U,t∈I(s,z)}.
有限增量定理
U⊂Rd开, f:U→Rd为C1光滑, [z0,z]⊂U, 则
∣f(z)−f(z0)−f′(z0)(z−z0)∣≤ξ∈[z0,z]sup∣∣f′(ξ)−f′(z0)∣∣∣z−z0∣
引理
U⊂Rd+1, f:U→Rd是C1的, (s,z)∈U, I(s,z)的解ϕ(t,s,z)的定义域为W, W={(t,s,z):(s,z)∈U,t∈I(s,z)}为开集, 则∂z∂ϕ在W上存在且连续.
证明
令Γ={(t,ϕ(t,s,z)):t∈[s,t0]}, M=max(t,x)∈Γ(ϵ0)∣∂x∂f(t,x)∣, ∂x∂f的连续模
δ(ϵ)=sup{∣∣∂x∂f(p)−∂x∂f(q)∣∣:p,q∈Γ(ϵ0),∣p−q∣<ϵ}
考虑ϕ(t0,s,z′)−ϕ(t0,s,z), 取δ>0, 当∣z′−z∣<δ, ϕ(t,s,z′)至少在[s,t0]上存在且位于Γ(ϵ0)内, 且
∣ϕ(t,s,z′)=ϕ(t,s,z)∣≤∣z′−z∣
D(t,s,z)是{X˙=∂x∂f(t,ϕ(t,s,z))XX(s)=Id的解, 而
D(t,s,z)=Id+∫st∂x∂f(τ,ϕ(τ,s,z))D(τ,s,z)dτ
ϕ(t,s,z′)−ϕ(t,s,z)=(z′−z)+∫stf(τ,ϕ(τ,s,z′))−f(τ,ϕ(τ,s,z))dτ
那么
η(t):=ϕ(t,s,z′)−ϕ(t,s,z)−D(t,s,z)(z′−z)=(z′−z)+∫stf(τ,ϕ(τ,s,z′))−f(τ,ϕ(τ,s,z))dτ−[(z′−z)+∫st∂x∂f(τ,ϕ(τ,s,z))D(τ,s,z)(z′−z)dτ]=∫st∂x∂f(τ,ϕ(τ,s,z))(ϕ(τ,s,z′)−ϕ(τ,s,z)−D(τ,s,z)(z′−z))+Δ(z′−z)dτ
后半部分
∣Δ(z′,z)∣≤∣ϕ(τ,s,z′)−ϕ(τ,s,z)∣≤δ(∣ϕ(τ,s,z)−ϕ(τ,s,z)∣)≤C∣z−z′∣δ(c∣z−z′∣)
从而
η(t)≤∫stMη(τ)dτ+∫st∣Δ(z′,z)∣dτ≤∫stMη(τ)dτ+c∣z′−z∣δ(c∣z′−z∣)
即
η(t)≤c∣z′−z∣δ(c∣z′−z∣)eM∣t−s∣=ο(∣z′−z∣)
这就证明了结论.
对于∂s∂ϕ的存在性和光滑性:
∂s∂ϕ(t,s,z)=−f(s,ϕ(s,s,z))+∫st∂x∂f(τ,ϕ(τ,s,z))∂s∂ϕ(τ,s,z)dτ
考察微分方程{X˙=∂x∂f(t,ϕ(t,s,z))XX(s)=−f(s,z)即可.
类似地推导可以得到:
性质
∂s∂ϕ也在W上存在且连续.
定理
I(s,z)的解ϕ(t,s,z)定义在W上, 且ϕ∈C1(W), ∂t∂ϕ=f(t,ϕ(t,s,z)).
证明
∂z∂ϕ,∂s∂ϕ,∂t∂ϕ∈C0(W)
推论
W∈Rd+d′+1, f:(t,x,λ)∋W^→Rd, f∈C1(W^), ∀(s,z,λ)∈W^, 记I(s,z,λ)的解为ϕ(t,s,z,λ), 则ϕ的定义域为W={(t,s,z,λ):(s,z,λ)∈W^,t∈I(s,z,λ)}, 且ϕ∈C(W).
Ck 光滑情形
定理
(1) f∈C(k)(U;Rd), 则ϕ(t,s,z;f)∈C(k)(W);
(2) F∈C(k)(W^;Rd), 则ϕ(t,s,z,λ;F)∈C(k)(W^).
这里k∈N∪{∞}.
证明
k=1时成立;
我们归纳法证明的路径是(1)k→(2)k→(1)k+1→⋯
假设k成立, U⊂Rd+1, 先设f:U→Rd, f∈C(k+1), 那么
ϕ∈C(k+1)(W)⇔∂z∂ϕ,∂s∂ϕ,∂t∂ϕ∈C(k)(W)
∂t∂ϕ(t,s,z)=f(t,ϕ(t,s,z))∈C(k)
∂z∂ϕ(t,s,z)=D(t,s,z)
{X˙=∂x∂f(t,ϕ(t,s,z))X=F(t,X,s,z)X(s)=Id断言F∈C(k), 有D(t,s,z)=Φ(t,s,z,s).
∂s∂ϕ(t,s,z)=D^(t,s,z)
{X˙=∂x∂f(t,ϕ(t,s,z))X=F(t,X,s,z)X(τ)=y, 有D^(t,s,z)=Ψ(t,s,z,s,−f(s,z)).
从而k+1时成立, 归纳有结论成立.