Skip to main content

解的存在唯一性与延拓

存在唯一性定理和延拓定理是 ODE 的理论基石:前者保证了解的存在性和唯一性(在一定条件下),后者说明了解可以延拓到多远。

一阶方程的存在唯一性

初值问题

y=f(x,y),y(x0)=y0y' = f(x, y), \quad y(x_0) = y_0

Picard-Lindelöf 定理

定理:设 f(x,y)f(x, y) 在矩形区域 R={(x,y)xx0a,yy0b}R = \{(x, y) \mid |x - x_0| \leq a, |y - y_0| \leq b\} 上连续,且对 yy 满足 Lipschitz 条件f(x,y1)f(x,y2)Ly1y2,(x,y1),(x,y2)R|f(x, y_1) - f(x, y_2)| \leq L |y_1 - y_2|, \quad \forall (x, y_1), (x, y_2) \in R

则存在 h>0h > 0,使初值问题在 [x0h,x0+h][x_0 - h, x_0 + h] 上存在唯一解。

Lipschitz 条件的判定

实际检验时常用以下充分条件:

fy\dfrac{\partial f}{\partial y}RR 上连续,则 ffyy 满足局部 Lipschitz 条件。

存在区间估计

h=min ⁣(a,bM)h = \min\!\left(a, \dfrac{b}{M}\right),其中 M=max(x,y)Rf(x,y)M = \max_{(x, y) \in R} |f(x, y)|

例子:非 Lipschitz \Rightarrow 解不唯一

y=y,y(0)=0y' = \sqrt{|y|}, \quad y(0) = 0

f(x,y)=yf(x, y) = \sqrt{|y|}y=0y = 0 处不满足 Lipschitz 条件。事实上,除零解 y0y \equiv 0 外,还有:

y(x)={0,xc(xc)24,x>c(c0)y(x) = \begin{cases} 0, & x \leq c \\ \dfrac{(x-c)^2}{4}, & x > c \end{cases} \quad (c \geq 0)

无穷多个解!这显示了 Lipschitz 条件对唯一性的必要性。

Picard 迭代法

Picard 存在唯一性定理的证明本身也是构造解的方法——Picard 逐次逼近法

y0(x)y0,yn+1(x)=y0+x0xf(t,yn(t))dty_0(x) \equiv y_0, \qquad y_{n+1}(x) = y_0 + \int_{x_0}^{x} f(t, y_n(t))\,dt

序列 {yn(x)}\{y_n(x)\} 一致收敛到真解 y(x)y(x)

解的延拓定理

"解的延拓不会在有限区间内无故中断——它能一直延拓,直到触及边界或趋于无穷。"

定理:设 f(x,y)f(x, y) 在区域 DR2D \subseteq \mathbb{R}^2 上连续且对 yy 满足局部 Lipschitz 条件,y=φ(x)y = \varphi(x) 是初值问题在区间 II 上的解。则 φ\varphi 可延拓到极大存在区间 (α,β)(\alpha, \beta),当 xα+x \to \alpha^+xβx \to \beta^- 时,(x,φ(x))(x, \varphi(x)) 必趋近于 DD 的边界,或 φ(x)\|\varphi(x)\| \to \infty(爆炸)。

高阶与方程组的推广

高阶方程

nn 阶方程 y(n)=f(x,y,y,,y(n1))y^{(n)} = f(x, y, y', \ldots, y^{(n-1)}) 可通过引入 y1=y,y2=y,,yn=y(n1)y_1 = y, y_2 = y', \ldots, y_n = y^{(n-1)} 化为一阶方程组。若 ff 对其后 nn 个变量满足 Lipschitz 条件,则存在唯一性定理成立。

线性方程组

对于线性方程组 y=A(x)y+b(x)\mathbf{y}' = A(x)\mathbf{y} + \mathbf{b}(x),若 A(x)A(x)b(x)\mathbf{b}(x) 在区间 II 上连续,则解在整个 II 上存在唯一(全局存在),无需延拓定理。

提示:线性 ODE 的一大优势——解在系数连续的整个区间上全局存在,不存在有限时间爆炸。

存在唯一性定理的重要意义

  1. 保证我们算出的解是唯一的,不同方法得到同一解
  2. 说明给定初值后,解完全确定(确定性)
  3. 延拓定理告诉我们解"最多能走多远"
  4. 为数值方法提供理论基础