Skip to content

定理Theorem

Peano 存在定理与 Euler 折线

Peano existence theorem · Cauchy–Peano theorem · 佩亚诺存在定理

在有限维连续向量场中构造不逃出的 Euler 折线,由一致积分残差与紧性取得经典解,并区分存在性、整列收敛和解误差。

形式陈述 ​

考虑有限维常微分方程初值问题

y′(t)=f(t,y(t)),y(t0)=y0∈Rd,1≤d<∞.

设 a,b>0,f 在紧圆柱

Q=[t0,t0+a]×B―(y0,b)

上联合连续;球和以下向量距离都用欧氏范数。记

M=max(t,z)∈Q‖f(t,z)‖,T={min{a,b/M},M>0,a,M=0.

则在 I=[t0,t0+T] 上至少有一条经典解,满足

y(t)=y0+∫t0tf(s,y(s))ds,‖y(t)−y0‖≤M(t−t0)≤b.

解在区间内部为 C1,导数连续延伸到端点,端点的方程以单侧导数解释。若 f 定义在 (t0,y0) 的开邻域并连续,可在左右各取这样的圆柱,得到初始时刻两侧的局部解。

定理没有要求对状态变量的 Lipschitz 常数。Picard–Lindelöf 定理中的状态 Lipschitz 条件还能保证局部唯一;这里的连续性只提供至少一条解。下面不仅证明这条解存在,也给出 Euler 近似满足积分方程到什么程度。

直觉

Euler 折线在每个网点读取一次向量场,然后沿这个方向走到下一个网点。它在一小段内使用了稍旧的时间和状态,因此通常不精确满足微分方程。然而,只要所有箭头都不长于 M,这些折线就都不能走得太快。选好总时间后,它们也不会离开我们掌握向量场的区域。

共同速度界提供紧性,保证至少能挑出一批越来越细、彼此一致接近的折线。还有一步同样重要:必须证明它们遗漏的方程量趋于零。这个量适合写成积分残差,因为折线在拐点未必有导数,而积分恒等式在整个区间都有意义。

Euler 折线的共同速度与积分残差

图中每段直线只使用左端点的斜率,细网格也仍是近似曲线。证明控制的是

Rh(t)=yh(t)−y0−∫t0tf(s,yh(s))ds,

而不是预先知道某条精确解后量出的距离。

第一步:构造每一段,并保证读数始终合法 ​

取任意有限分割

t0=τ0<τ1<⋯<τN=t0+T,h=maxj(τj+1−τj).

令 yh(t0)=y0,并递推定义

yh(t)=yh(τj)+(t−τj)f(τj,yh(τj)),τj≤t≤τj+1.

构造要先有依据才能继续:若前面的网点已定义且留在球中,当前斜率的范数就至多为 M。从初值累加所有已完成段与当前段的位移,得到

‖yh(t)−y0‖≤M∑i<j(τi+1−τi)+M(t−τj)=M(t−t0)≤b.

由归纳,每个新网点及整段都留在球内,下一次 f 求值因此合法。对任意两个时刻,把它们之间的路程按网点分开,同样得到

‖yh(t)−yh(s)‖≤M|t−s|.

如果从固定步长生成网格,就把最后一段截成终点剩余的长度。上面的证明使用实际段长,并不要求 T/h 为整数。M=0 时,圆柱上 f≡0,所有折线就是常值 y0,定理直接成立。

第二步:把冻结斜率的偏差变成一致残差 ​

在 Q 上使用距离 |(t−s)|+‖z−w‖。紧集上的连续函数一致连续,因此可定义一个趋零的连续模

ω(r)=sup{‖f(t,z)−f(s,w)‖:(t,z),(s,w)∈Q, |t−s|+‖z−w‖≤r}.

有限维紧性保证这个上确界有限,且 ω(r)→0。对 s∈[τj,τj+1],时间偏移至多为 h,状态偏移至多为 Mh,所以

‖f(τj,yh(τj))−f(s,yh(s))‖≤ω((1+M)h).

令 σh(s) 是所在网格段的左端点;在最后端点任意选相邻段,不影响积分。逐段积分折线的常斜率,得到严格恒等式

yh(t)−y0=∫t0tf(σh(s),yh(σh(s)))ds.

两个积分相减,便有

‖Rh‖∞≤Tω((1+M)h)⟶0.

这里只比较相近输入处的连续向量场,没有引入状态 Lipschitz 常数。

第三步:取得子列,并把方程传给极限 ​

取一列最大步长 hn→0 的分割。全部折线共同有界且有共同 Lipschitz 常数 M。对第一个坐标使用Arzelà–Ascoli 定理抽取一致收敛子列,再对第二个坐标从该子列继续抽取,依次进行有限次。最后一个子列的全部 d 个坐标同时一致收敛。由

‖v‖2≤dmax1≤i≤d|vi|

可知存在 yhnk→y 于向量值一致范数。一致极限保连续,而对闭球约束取极限,得到 y(t)∈B―(y0,b)。

记 ek=‖yhnk−y‖∞。同一连续模给出

sups∈I‖f(s,yhnk(s))−f(s,y(s))‖≤ω(ek)⟶0.

于是,对每个 t∈I,甚至对整个区间同时成立,

‖y(t)−y0−∫t0tf(s,y(s))ds‖≤ek+‖Rhnk‖∞+Tω(ek)⟶0.

因此积分方程精确成立。被积函数连续,逐坐标应用微积分基本定理,就得到经典微分方程及初值。证明没有对折线的导数直接取极限。

例子与边界

连续平方根右端的完整认证 ​

考虑

y′(t)=t+|y(t)|,y(0)=0,0≤t≤14.

在矩形 Q=[0,1/4]×[−1,1] 上,右端连续且

|t+|z||≤14+1=54=M.

取 b=1,则 MT=5/16<1,所以任意分割产生的 Euler 折线都留在矩形中;所有折线有共同速度界 5/4。事实上斜率非负,折线还满足 0≤yh(t)≤5/16,但存在性证明不需要这份更小的范围。

对任意实数 u,v,由绝对值的三角不等式和平方根差的估计,有

||u|−|v||≤||u|−|v||≤|u−v|.

所以在每个网格段,冻结斜率与当前右端的差至多为

h+54h.

积分长度最多为 1/4,得到可直接检查的证书

sup0≤t≤1/4|yh(t)−∫0t(s+|yh(s)|)ds|≤14(h+54h).

右侧趋于零,共同界又给一致收敛子列,所以极限满足这个具体方程。这里 f(t,⋅) 在零附近不 Lipschitz:正数 z 与零的差商是 1/z。本证明因而没有借用 Picard 压缩条件。

若想让这份积分残差证书不超过 1/100,可取 h≤1/1000。因为

54⋅11000=1800<391000,

其中平方后是 1250<1521,所以右侧严格小于 1/100。这只是一个方便的充分步长,不是最佳值,更不是“距某条指定解小于 1/100”的声明。

时间向左走时如何改写 ​

若连续圆柱位于 [t0−a,t0]×B―(y0,b),令 s=t0−t、z(s)=y(t0−s)。方程变成

z′(s)=−f(t0−s,z(s)),z(0)=y0.

新的右端仍连续,并有同一个高度界 M,因此正向证明完全适用。再把参数换回来,就得到向后的解。两侧圆柱都可用时,在 t0 拼接各自的解;两边导数都等于 f(t0,y0),拼接后仍是 C1。

紧性留下子列,不能代替唯一性或稳定性 ​

上面的结论是:每列越来越细的合法分割,都有一个子列一致趋于某条解。它没有指定抽取哪些下标,也没有控制整个序列到一条选定解的速度。连续右端可以有多解,具体等待解见Picard–Lindelöf 页的例子与边界。

若另外已经证明:在这个函数球中解唯一,那么整列折线便可由紧性升级为一致收敛。反证:若整列不趋于这条唯一解,就有一个子列始终距该解至少某个 ε>0;再对它抽取收敛子列,上面的残差证明使其极限也是解,与唯一性矛盾。这个推理仍不给显式速度。把残差换成解误差,需要另外的稳定性估计;Euler 方法中的全局一阶误差使用了状态 Lipschitz 和精确解的额外光滑性。

状态有限维也是实质条件。逐坐标抽取在有限次后能够结束;一般无限维 Banach 值域里,有界并不保证各点评值相对紧。常值函数 t↦en 就共同有界、共同零 Lipschitz,却彼此保持正距离,详见Ascoli 页的值域边界。不能把本页定理无条件换成 Banach 空间中的 ODE 存在定理。

推论与应用

迁移到耦合向量场 ​

例如初值为 (0,0) 的二维系统

x1′=t+|x2|,x2′=|x1|

在 [0,1/4]×B―(0,1) 上连续。欧氏范数下可以取 M=(5/4)2+1=41/4,且 MT=41/16<1。同样的折线构造与有限次坐标抽取给出经典向量解。

如果还需要显式残差,对同一网格段使用 ‖xh(s)−xh(τj)‖2≤Mh,便有

‖Rh‖∞≤14(h+Mh)2+Mh.

这里重新合成了两个坐标的误差,不能直接把标量例中的常数 5/4 当成整个向量场的欧氏范数界。

从这份证明读出存在性任务的三个关口 ​

给近似对象一个共同控制,使它们有合法极限;证明它们对原方程的残差一致消失;最后检查方程在所选极限下保持。Euler 更新只是产生近似对象的一种方式。Schauder 紧像不动点定理使用有限维凸包上的近似不动点,随后也要靠紧性与残差闭合极限。

可在连续方程的存在性终点任务中检查不足整步的末段、平方根残差预算及向量迁移;这些检查不需要先写出精确解。

参考资料
  • Gerald Teschl,Ordinary Differential Equations and Dynamical Systems, §2.7,AMS,2012 年作者预印本,印刷页 54–56,Theorem 2.19、式 (2.71)–(2.76):Euler 折线、共同速度界、积分残差和 Peano 存在证明。本文使用任意分割,并逐项处理零高度界、末段和反向时间。
  • 本文平方根右端、残差步长预算与二维耦合系统的常数均由正文不等式独立计算;Ascoli 和一致极限的完整证明见所链接条目。
关系图谱13 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系