本文我们研究形如 y ′ = f ( x , y ) y' = f(x,y) y ′ = f ( x , y ) 的一阶显式常微分方程,分析其初值问题的解在局部到全局的存在性和唯一性。
{ y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0
我们称上述方程在 ( x 0 , y 0 ) ∈ D ⊂ R 2 (x_0, y_0) \in D \subset \mathbb{R}^2 ( x 0 , y 0 ) ∈ D ⊂ R 2 上有解意味着有一条解曲线 y ( x ) ⊂ D y(x) \subset D y ( x ) ⊂ D 即可,这不要求 y ( x ) y(x) y ( x ) 的定义域充满 D x D_x D x ,有唯一解同理叙述;而称上述方程在 x 0 ∈ I ⊂ R x_0 \in I \subset \mathbb{R} x 0 ∈ I ⊂ R 上有解意味着解曲线 y ( x ) y(x) y ( x ) 要在 I I I 上的每一点有定义且符合题意。
根据定义我们知道,上述方程的解 y ( x ) y(x) y ( x ) 一定要是连续可导的,且过点 ( x 0 , y 0 ) (x_0, y_0) ( x 0 , y 0 ) ,从而我们希望在 x y xy x y 平面上找到一个过点 ( x 0 , y 0 ) (x_0, y_0) ( x 0 , y 0 ) 且满足 y ′ = f ( x , y ) y' = f(x,y) y ′ = f ( x , y ) 的函数曲线。
但研究问题的时候我们通常不会直接在全局上找解,即把解的定义域直接限制在 R \mathbb{R} R 上,这意味着上述方程的解必须在 R \mathbb{R} R 上的每个点都有定义,这样的条件太苛刻了,故我们一般会先从 x 0 x_0 x 0 的一个小邻域内去找是否有解以及解是否唯一。
如果说在任何 x 0 x_0 x 0 的邻域内均没有解,那么显然在全局上就没有解,如果说在某个邻域有解,则我们可以找到一个最大解区间 I I I ,至于这个 I I I 是否可以取到 R \mathbb{R} R 是不确定的,从局部小邻域到最大解区间,这其中实际上蕴含了解的延拓思想。从而我们先希望在一个小区域 ( x 0 , y 0 ) ∈ D ⊂ R 2 (x_0, y_0) \in D \subset \mathbb{R}^2 ( x 0 , y 0 ) ∈ D ⊂ R 2 中所有过点 ( x 0 , y 0 ) (x_0, y_0) ( x 0 , y 0 ) 的连续可导函数里找到一个符合 y ′ = f ( x , y ) y' = f(x,y) y ′ = f ( x , y ) 的函数曲线,再看这个小区域里的解曲线是否是唯一的。
研究问题往往从特殊再到一般,我们不妨先研究 f f f 是二元连续函数,那么 y ′ y' y ′ 也是连续函数,从而可积,那么原来的初值问题等价于下述的积分方程:
y ( x ) = y 0 + ∫ x 0 x f ( t , , y ( t ) ) d t , x ∈ J y(x) = y_0 + \int_{x_0}^x f(t,, y(t))\,dt, \quad x \in J y ( x ) = y 0 + ∫ x 0 x f ( t ,, y ( t )) d t , x ∈ J
这里的定义域 J J J 先待定,或者可以理解为是一个很小的 x 0 x_0 x 0 的邻域。我们实际上就是解上述积分方程,不妨定义一个映射:
T : C ( J ) → C ( J ) , ( T y ) ( x ) = y 0 + ∫ x 0 x f ( t , , y ( t ) ) d t , x ∈ J T: C(J) \to C(J), \quad (Ty)(x) = y_0 + \int_{x_0}^x f(t,, y(t))\,dt, \quad x \in J T : C ( J ) → C ( J ) , ( T y ) ( x ) = y 0 + ∫ x 0 x f ( t ,, y ( t )) d t , x ∈ J
那么我们的目标就变成了找这个映射在 C ( J ) C(J) C ( J ) 上的不动点,这里的 C ( J ) C(J) C ( J ) 为 J J J 上的所有连续可导函数形成的空间,那么原初值问题在 J J J 上解的存在性和唯一性就变成了不动点的存在性和唯一性。值得注意的是这个映射 T T T 是将 J J J 上的任意连续可导函数都映射到了 J J J 上过点 ( x 0 , y 0 ) (x_0, y_0) ( x 0 , y 0 ) 的连续可导函数。
1. 解的唯一存在性
本节我们讨论更特殊的如下情形:设 D = [ x 0 − a , , x 0 + a ] × [ y 0 − b , , y 0 + b ] D = [x_0 - a,, x_0 + a] \times [y_0 - b,, y_0 + b] D = [ x 0 − a ,, x 0 + a ] × [ y 0 − b ,, y 0 + b ] ,a , b a, b a , b 均大于零,假设二元函数 f f f 在 D D D 上连续且满足关于 y y y 的 Lipschitz 条件,即存在常数 L L L ,使得
∣ f ( x , y 1 ) − f ( x , y 2 ) ∣ ≤ L ∣ y 1 − y 2 ∣ , ∀ ( x , y 1 ) , ( x , y 2 ) ∈ D |f(x,y_1) - f(x,y_2)| \leq L|y_1 - y_2|, \quad \forall\, (x,y_1), (x,y_2) \in D ∣ f ( x , y 1 ) − f ( x , y 2 ) ∣ ≤ L ∣ y 1 − y 2 ∣ , ∀ ( x , y 1 ) , ( x , y 2 ) ∈ D
我们接下来在区域 D D D 上研究不动点的存在性和唯一性,以及解曲线的定义域 J ⊂ [ x 0 − a , , x 0 + a ] J \subset [x_0 - a,, x_0 + a] J ⊂ [ x 0 − a ,, x 0 + a ] 的可取范围。在研究之前,我们先介绍压缩映射原理和解的局部延拓相关知识来得到关于解的局部存在唯一性的 Cauchy-Lipschitz 定理,最后通过全局延拓探究解在全局的存在唯一性。
1.1 压缩映射原理
设 B B B 为完备的赋范线性空间,即 Banach 空间,D ⊂ B D \subset B D ⊂ B 为非空闭集,若对于映射 T : D → B T: D \to B T : D → B ,存在 q ∈ ( 0 , 1 ) q \in (0,1) q ∈ ( 0 , 1 ) ,使得
∣ T ( x ) − T ( y ) ∣ ≤ q ∣ x − y ∣ , ∀ x , y ∈ D |T(x) - T(y)| \leq q|x - y|, \quad \forall\, x,y \in D ∣ T ( x ) − T ( y ) ∣ ≤ q ∣ x − y ∣ , ∀ x , y ∈ D
则我们有 T T T 在 D D D 上有唯一的不动点,即存在 x ˉ ∈ D \bar{x} \in D x ˉ ∈ D 使得 T ( x ˉ ) = x ˉ T(\bar{x}) = \bar{x} T ( x ˉ ) = x ˉ 成立。
更确切地说,任取 x 0 ∈ D x_0 \in D x 0 ∈ D ,令 x n + 1 = T ( x n ) x_{n+1} = T(x_n) x n + 1 = T ( x n ) ,∀ n ≥ 0 \forall\, n \geq 0 ∀ n ≥ 0 ,则序列 x n {x_n} x n 总会收敛于 x ˉ \bar{x} x ˉ ,且
∣ x n − x ˉ ∣ ≤ 1 1 − q ∣ x n + 1 − x n ∣ ≤ q n 1 − q ∣ x 1 − x 0 ∣ |x_n - \bar{x}| \leq \frac{1}{1-q}|x_{n+1} - x_n| \leq \frac{q^n}{1-q}|x_1 - x_0| ∣ x n − x ˉ ∣ ≤ 1 − q 1 ∣ x n + 1 − x n ∣ ≤ 1 − q q n ∣ x 1 − x 0 ∣
1.2 解的局部延拓
如果 ϕ 1 ( x ) \phi_1(x) ϕ 1 ( x ) 为初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 在 J 1 = [ x 0 , b ] J_1 = [x_0, b] J 1 = [ x 0 , b ] 上的唯一解,ϕ 2 ( x ) \phi_2(x) ϕ 2 ( x ) 为初值问题 { y ′ = f ( x , y ) y ( b ) = ϕ 1 ( b ) \begin{cases} y' = f(x,y) \\ y(b) = \phi_1(b) \end{cases} { y ′ = f ( x , y ) y ( b ) = ϕ 1 ( b ) 在 J 2 = [ b , c ] J_2 = [b,c] J 2 = [ b , c ] 上的唯一解,可以证明
ϕ ( x ) = { ϕ 1 ( x ) , x ∈ [ x 0 , b ] ϕ 2 ( x ) , x ∈ [ b , c ] \phi(x) = \begin{cases} \phi_1(x), & x \in [x_0, b] \\ \phi_2(x), & x \in [b, c] \end{cases} ϕ ( x ) = { ϕ 1 ( x ) , ϕ 2 ( x ) , x ∈ [ x 0 , b ] x ∈ [ b , c ]
为初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 在 J = [ x 0 , c ] J = [x_0, c] J = [ x 0 , c ] 上的唯一解,事实上这只需要证明在点 b b b 处的导数存在且满足导数等于 f ( x , y ) f(x,y) f ( x , y ) 即可。
根据上述结论,我们知道如果我们想研究 J = [ x 0 , x ˉ ] J = [x_0, \bar{x}] J = [ x 0 , x ˉ ] 上的唯一解,那么可以先研究 [ x 0 , x 1 ] [x_0, x_1] [ x 0 , x 1 ] 上的唯一解,再由这个解构建新的初值问题,研究在 [ x 1 , x 2 ] [x_1, x_2] [ x 1 , x 2 ] 上的唯一解,以此类推,这里的 x 1 , x 2 , ⋯ x_1, x_2, \cdots x 1 , x 2 , ⋯ 均小于等于 x n x_n x n ,其取值依赖于你能够求得出的初值问题的解区间,至于 x n x_n x n 是否能达到 x ˉ \bar{x} x ˉ 是不确定的。
更一般地,如果说 ϕ 1 ( x ) \phi_1(x) ϕ 1 ( x ) 为初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 在 J 1 J_1 J 1 上的唯一解,ϕ 2 ( x ) \phi_2(x) ϕ 2 ( x ) 为初值问题 { y ′ = f ( x , y ) y ( x 1 ) = y 1 \begin{cases} y' = f(x,y) \ y(x_1) = y_1 \end{cases} { y ′ = f ( x , y ) y ( x 1 ) = y 1 在 J 2 J_2 J 2 上的唯一解,且 J 1 ∩ J 2 ≠ ∅ J_1 \cap J_2 \neq \emptyset J 1 ∩ J 2 = ∅ ,ϕ 1 ( x ) = ϕ 2 ( x ) \phi_1(x) = \phi_2(x) ϕ 1 ( x ) = ϕ 2 ( x ) 对于某个 x ∈ J 1 ∩ J 2 x \in J_1 \cap J_2 x ∈ J 1 ∩ J 2 成立,那么初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 在 J 1 ∪ J 2 J_1 \cup J_2 J 1 ∪ J 2 上存在唯一解 ϕ ( x ) \phi(x) ϕ ( x ) ,使得 ϕ ∣ J n = ϕ n , ; n = 1 , 2 \phi|_{J_n} = \phi_n,; n = 1,2 ϕ ∣ J n = ϕ n ,; n = 1 , 2 成立。
1.3 Cauchy-Lipschitz 定理
根据压缩映射原理,自然我们会希望映射 T : C ( J ) → C ( J ) T: C(J) \to C(J) T : C ( J ) → C ( J ) 为压缩映射,这样它在 J J J 上就有唯一不动点,从而原来的初值问题在 J J J 上有唯一解,我们根据开头给出的条件,来分析一下映射 T T T 或者说找到一个合适的 J J J 使得 T T T 为压缩映射。任取 y , z ∈ C ( J ) y, z \in C(J) y , z ∈ C ( J ) ,
∣ T y ( x ) − T z ( x ) ∣ = ∣ ∫ x 0 x f ( t , y ( t ) ) − f ( t , z ( t ) ) d t ∣ ≤ ∫ x 0 x ∣ f ( t , y ( t ) ) − f ( t , z ( t ) ) ∣ d t ≤ L ∫ x 0 x ∣ y ( t ) − z ( t ) ∣ d t ≤ L ( x − x 0 ) ∣ y − z ∣ |T_y(x) - T_z(x)| = \left|\int_{x_0}^x f(t,y(t)) - f(t,z(t))\,dt\right| \leq \int_{x_0}^x |f(t,y(t)) - f(t,z(t))|\,dt \leq L\int_{x_0}^x |y(t) - z(t)|\,dt \leq L(x-x_0)|y-z| ∣ T y ( x ) − T z ( x ) ∣ = ∫ x 0 x f ( t , y ( t )) − f ( t , z ( t )) d t ≤ ∫ x 0 x ∣ f ( t , y ( t )) − f ( t , z ( t )) ∣ d t ≤ L ∫ x 0 x ∣ y ( t ) − z ( t ) ∣ d t ≤ L ( x − x 0 ) ∣ y − z ∣
这里空间 C ( J ) C(J) C ( J ) 的范数为 ∣ y ∣ = max x ∈ J ∣ y ( x ) ∣ |y| = \max_{x \in J}|y(x)| ∣ y ∣ = max x ∈ J ∣ y ( x ) ∣ ,上述过程为了使得第二个不等式成立,我们需要让
y ( t ) ∈ [ y 0 − b , y 0 + b ] , ∀ t ∈ [ x 0 , x ] ⟺ ∣ y ( t ) − y 0 ∣ ≤ b , ∀ t ∈ [ x 0 , x ] y(t) \in [y_0 - b, y_0 + b],\quad \forall\, t \in [x_0, x] \iff |y(t) - y_0| \leq b,\quad \forall\, t \in [x_0, x] y ( t ) ∈ [ y 0 − b , y 0 + b ] , ∀ t ∈ [ x 0 , x ] ⟺ ∣ y ( t ) − y 0 ∣ ≤ b , ∀ t ∈ [ x 0 , x ]
而根据微分中值定理,∣ y ( t ) − y ( x 0 ) ∣ = ∣ y ′ ( ξ ) ∣ ⋅ ∣ t − x 0 ∣ ≤ A ⋅ ∣ t − x 0 ∣ |y(t) - y(x_0)| = |y'(\xi)| \cdot |t - x_0| \leq A \cdot |t - x_0| ∣ y ( t ) − y ( x 0 ) ∣ = ∣ y ′ ( ξ ) ∣ ⋅ ∣ t − x 0 ∣ ≤ A ⋅ ∣ t − x 0 ∣ ,其中 A = max x , y ∈ D ∣ f ( x , y ) ∣ A = \max_{x,y \in D}|f(x,y)| A = max x , y ∈ D ∣ f ( x , y ) ∣ ,从而只需要 ∣ t − x 0 ∣ < b A |t - x_0| < \dfrac{b}{A} ∣ t − x 0 ∣ < A b 即可,那么也就是说在区间 [ x 0 − α , , x 0 + α ] [x_0 - \alpha,, x_0 + \alpha] [ x 0 − α ,, x 0 + α ] 内(其中 α = min { a , , b A } \alpha = \min\{a,, \dfrac{b}{A}\} α = min { a ,, A b } ),有
∣ T ( y ) − T ( z ) ∣ ≤ L α ∣ y − z ∣ |T(y) - T(z)| \leq L\alpha|y - z| ∣ T ( y ) − T ( z ) ∣ ≤ Lα ∣ y − z ∣
如果 L α < 1 L\alpha < 1 Lα < 1 ,那么在区间 [ x 0 − α , , x 0 + α ] [x_0 - \alpha,, x_0 + \alpha] [ x 0 − α ,, x 0 + α ] 上初值问题有唯一解;
如果 L α ≥ 1 L\alpha \geq 1 Lα ≥ 1 ,我们不妨先在一个更小的区间 [ x 0 , , x 0 + b ] [x_0,, x_0 + b] [ x 0 ,, x 0 + b ] 分析,其中 b = L α n < 1 b = \dfrac{L\alpha}{n} < 1 b = n Lα < 1 ,那么在该小区间 T T T 为压缩映射,原初值问题有唯一解 y 1 ( x ) y_1(x) y 1 ( x ) ,同理在区间 [ x 0 + b , x 0 + 2 b ] [x_0 + b, x_0 + 2b] [ x 0 + b , x 0 + 2 b ] 上,T ′ ( y ) = y 1 ( x 0 + b ) + ∫ x 0 x f ( t , y ( t ) ) d t T'(y) = y_1(x_0 + b) + \int_{x_0}^x f(t, y(t))\,dt T ′ ( y ) = y 1 ( x 0 + b ) + ∫ x 0 x f ( t , y ( t )) d t 仍为压缩映射,从而我们找到初值问题 { y ′ = f ( x , y ) y ( x 0 + b ) = y 1 ( x 0 + b ) \begin{cases} y' = f(x,y) \\ y(x_0 + b) = y_1(x_0 + b) \end{cases} { y ′ = f ( x , y ) y ( x 0 + b ) = y 1 ( x 0 + b ) 也有唯一解,再由解的延拓中的结论,得到在区间 [ x 0 , x 0 + 2 b ] [x_0, x_0 + 2b] [ x 0 , x 0 + 2 b ] 上原初值问题有唯一解,以此类推,得到在区间 [ x 0 , , x 0 + α ] [x_0,, x_0 + \alpha] [ x 0 ,, x 0 + α ] 上有唯一解,同理 [ x 0 − α , , x 0 ] [x_0 - \alpha,, x_0] [ x 0 − α ,, x 0 ] 上有唯一解,最终就得到 [ x 0 − α , , x 0 + α ] [x_0 - \alpha,, x_0 + \alpha] [ x 0 − α ,, x 0 + α ] 上初值问题有唯一解。
实际上我们根据压缩映射原理,可以有计算初值问题解的数值方法,在满足条件的定义域内构造逐次逼近序列
y k + 1 ( x ) = y 0 + ∫ x 0 x f ( t , , y k ( t ) ) d t , ∀ k y_{k+1}(x) = y_0 + \int_{x_0}^x f(t,, y_k(t))\,dt, \quad \forall\, k y k + 1 ( x ) = y 0 + ∫ x 0 x f ( t ,, y k ( t )) d t , ∀ k
取 y 0 ( x ) = y 0 y_0(x) = y_0 y 0 ( x ) = y 0 ,例如上述如果 L α ≥ 1 L\alpha \geq 1 Lα ≥ 1 ,我们则需要分区间进行逼近。
综上我们得到如下的 Cauchy-Lipschitz 定理 :设 D = [ x 0 − a , , x 0 + a ] × [ y 0 − b , , y 0 + b ] D = [x_0 - a,, x_0 + a] \times [y_0 - b,, y_0 + b] D = [ x 0 − a ,, x 0 + a ] × [ y 0 − b ,, y 0 + b ] ,a , b a,b a , b 均大于零,假设二元函数 f f f 在 D D D 上连续且满足关于 y y y 的 Lipschitz 条件,则初值问题
{ y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0
在 D D D 上有唯一解,且解的存在区间 J J J 至少包含区间 [ x 0 − α , , x 0 + α ] [x_0 - \alpha,, x_0 + \alpha] [ x 0 − α ,, x 0 + α ] ,其中
α = min { a , b A } , A = max x , y ∈ D ∣ f ( x , y ) ∣ \alpha = \min\{a, \frac{b}{A}\}, \quad A = \max_{x,y \in D}|f(x,y)| α = min { a , A b } , A = x , y ∈ D max ∣ f ( x , y ) ∣
1.4 解的全局延拓
根据上述分析,我们已经知道如果二元连续函数 f f f 在平面上每个点都满足关于 y y y 的局部 Lipschitz 条件,那么 y ′ = f ( x , y ) y' = f(x,y) y ′ = f ( x , y ) 关于每一点 ( x 0 , y 0 ) (x_0, y_0) ( x 0 , y 0 ) 的初值问题在 x 0 x_0 x 0 的一个小邻域一定有唯一解,或者说在平面上一定唯一解曲线。那么一个自然的想法就是不断进行局部延拓,是否可以得到全局的解?
具体思路如下:我们首先根据 y ′ = f ( x , y ) y' = f(x,y) y ′ = f ( x , y ) 关于点 ( x 0 , y 0 ) (x_0, y_0) ( x 0 , y 0 ) 的初值问题得到唯一解区间 [ x 0 − a 1 , , x 0 + a 1 ] [x_0 - a_1,, x_0 + a_1] [ x 0 − a 1 ,, x 0 + a 1 ] 和对应的解函数 y 1 ( x ) y_1(x) y 1 ( x ) ,然后根据 y ′ = f ( x , y ) y' = f(x,y) y ′ = f ( x , y ) 关于点 ( x 0 + a 1 , , y 1 ( x 0 + a 1 ) ) (x_0 + a_1,, y_1(x_0 + a_1)) ( x 0 + a 1 ,, y 1 ( x 0 + a 1 )) 的初值问题得到唯一解区间 [ x 0 + a 1 − a 2 , , x 0 + a 1 + a 2 ] [x_0 + a_1 - a_2,, x_0 + a_1 + a_2] [ x 0 + a 1 − a 2 ,, x 0 + a 1 + a 2 ] ,从而得到原初值问题在 [ x 0 , , x 0 + a 1 + a 2 ] [x_0,, x_0 + a_1 + a_2] [ x 0 ,, x 0 + a 1 + a 2 ] 上有唯一解,左侧同理,如此一直进行下去,记 S n = ∑ k = 1 n a k S_n = \sum_{k=1}^n a_k S n = ∑ k = 1 n a k ,如果说 lim n → ∞ S n = ∞ \lim_{n \to \infty} S_n = \infty lim n → ∞ S n = ∞ ,那么原初值问题自然在全局有唯一存在性,否则,便无法延拓到全局。
同样以右侧为例,如果初值问题的解无法延拓到正无穷,那么它一定存在一个最大解区间,记为 [ x 0 , b ) [x_0, b) [ x 0 , b ) ,如果可以取到 b b b ,那么自然可以继续向右延拓,又因为倘若这个解区间对应的解函数 y ( x ) y(x) y ( x ) 在其中有界,则 ( x , y ( x ) ) (x, y(x)) ( x , y ( x )) 有界,则 f ( x , y ( x ) ) f(x, y(x)) f ( x , y ( x )) 有界,则 ∣ y ′ ∣ |y'| ∣ y ′ ∣ 在这个解区间有界,从而 y ( x ) y(x) y ( x ) 在 [ x 0 , b ) [x_0, b) [ x 0 , b ) 上是一致连续的,从而 b b b 处的左极限是存在且有限的,从而根据积分方程令 x → b − x \to b^- x → b − 可得解区间可延拓到 [ x 0 , b ] [x_0, b] [ x 0 , b ] ,从而如果初值问题的解无法延拓到正无穷,则 lim x → b − y ( x ) = ∞ \lim_{x \to b^-} y(x) = \infty lim x → b − y ( x ) = ∞ 。
如果我们假设的是二元连续函数 f f f 在开集 D ⊂ R 2 D \subset \mathbb{R}^2 D ⊂ R 2 上满足每个点的关于 y y y 的局部 Lipschitz 条件,同样任给 ( x 0 , y 0 ) ∈ D (x_0, y_0) \in D ( x 0 , y 0 ) ∈ D ,易知初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 在 D D D 上存在唯一解 ϕ ( x ) \phi(x) ϕ ( x ) ,并且 ϕ ( x ) \phi(x) ϕ ( x ) 对应的右侧最大解区间形如 [ x 0 , b ) [x_0, b) [ x 0 , b ) ,并且我们有下列三种情形之一必成立:
b = ∞ b = \infty b = ∞ ,即全局有唯一解
b < ∞ b < \infty b < ∞ ,且 lim x → b − y ( x ) = ∞ \lim_{x \to b^-} y(x) = \infty lim x → b − y ( x ) = ∞ ,即可以一直延拓,且解曲线一直在 D D D 内,但无法延拓到全局
b < ∞ b < \infty b < ∞ ,且 lim x → b − d ∂ D ( x , , ϕ ( x ) ) = 0 \lim_{x \to b^-} d_{\partial D}(x,, \phi(x)) = 0 lim x → b − d ∂ D ( x ,, ϕ ( x )) = 0 ,即解曲线到达 D D D 的边界,而 D D D 之外不一定满足 Lipschitz 条件,从而无法向外延拓
事实上上述的三个情形合起来就等价于 G = { ( x , ϕ ( x ) ) ∣ x ≥ x 0 } G = \{(x,\phi(x)) \mid x \geq x_0\} G = {( x , ϕ ( x )) ∣ x ≥ x 0 } 为非有界闭集,即解函数在最大解区间对应的像为非有界闭集。事实上如果将无穷远处也看做 D D D 的边界,那么就相当于唯一解一定能延拓到 D D D 的边界。
当然如果 f f f 满足的是全局 Lipschitz 条件,那么解曲线就不会在有限时间爆炸,从而解可以向左或向右延拓到边界,即上述的第二个情形不会成立。
1.5 Osgood 条件
设函数 f ( x , y ) f(x,y) f ( x , y ) 在区域 G G G 内连续,而且满足不等式 ∣ f ( x , y 1 ) − f ( x , y 2 ) ∣ ≤ F ( ∣ y 1 − y 2 ∣ ) |f(x,y_1) - f(x,y_2)| \leq F(|y_1 - y_2|) ∣ f ( x , y 1 ) − f ( x , y 2 ) ∣ ≤ F ( ∣ y 1 − y 2 ∣ ) ,其中 F ( r ) > 0 F(r) > 0 F ( r ) > 0 是 r > 0 r > 0 r > 0 上的连续函数,并且瑕积分 ∫ 0 r 1 d r F ( r ) = ∞ \int_0^{r_1} \dfrac{dr}{F(r)} = \infty ∫ 0 r 1 F ( r ) d r = ∞ ,其中 r 1 > 0 r_1 > 0 r 1 > 0 为常数,则称 f ( x , y ) f(x,y) f ( x , y ) 在 G G G 内对 y y y 满足 Osgood 条件。
容易证明关于 y y y 的 Lipschitz 条件是 Osgood 条件的特例。
并且我们有更一般的满足对 y y y 的 Osgood 条件的解的唯一存在性的结论:若 f ( x , y ) f(x,y) f ( x , y ) 在区域 G G G 内满足对 y y y 的 Osgood 条件,那么初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 在 G G G 内有唯一解,且可以延拓到 G G G 的边界。解的存在性可以根据下面的 Peano 定理得出,唯一性采取反证法即可。
2. 解的存在性
前面我们对二元连续函数 f f f 施加了局部关于 y y y 的 Lipschitz 条件来研究初值问题的解的唯一存在性以及最大解区间的范围,那么我们现在丢掉这个 Lipschitz 条件来分析解是否存在,即用一个弱一点的条件来研究弱一点的结论。
我们先来观察这个初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 ,分析 y ′ = f ( x , y ) y' = f(x,y) y ′ = f ( x , y ) ,不难发现我们可以理解为在二维平面上,每一点都唯一对应一个斜率,而方程的解就是过点 ( x 0 , y 0 ) (x_0, y_0) ( x 0 , y 0 ) 的且每点斜率吻合的连续可导曲线,于是我们自然会有一个想法,从给定的点 ( x 0 , y 0 ) (x_0, y_0) ( x 0 , y 0 ) 出发,该点的斜率我们知道,我们就可以向两边作极小的线性延伸,然后我们会得到一个新的点 ( x 1 , y 1 ) (x_1, y_1) ( x 1 , y 1 ) ,根据这个点我们又得到该点的斜率,再作延伸,无限进行下去,似乎我们就可以得到一个近似解,只要不断减小每次延伸的长度,似乎就能提高近似的精度,这在理论上其实需要证明这一列折线函数在每次延伸的长度趋于零时是收敛的,并且收敛到初值问题的解。
事实上上述的思想其实就是欧拉折线,下面我们更严谨地来定义欧拉折线。
2.1 欧拉折线
对于初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 ,设 f ( x , y ) f(x,y) f ( x , y ) 是在矩形区域 D = [ x 0 − a , x 0 + a ] × [ y 0 − b , y 0 + b ] D = [x_0 - a, x_0 + a] \times [y_0 - b, y_0 + b] D = [ x 0 − a , x 0 + a ] × [ y 0 − b , y 0 + b ] 内给定的连续函数,则初值问题同样等价于前面的积分方程,且为了保证解曲线一定落在 D D D 内部,我们加以一定的限制,即令 A = max x , y ∈ D ∣ f ( x , y ) ∣ A = \max_{x,y \in D}|f(x,y)| A = max x , y ∈ D ∣ f ( x , y ) ∣ ,h = min { a , b M } h = \min\{a, \frac{b}{M}\} h = min { a , M b } ,则在区间 [ x 0 − h , , x 0 + h ] [x_0 - h,, x_0 + h] [ x 0 − h ,, x 0 + h ] 上,原初值问题若有解则解曲线一定落在 D D D 内。接着,我们将区间 ∣ x − x 0 ∣ ≤ h |x - x_0| \leq h ∣ x − x 0 ∣ ≤ h 分成 2 n 2n 2 n 等份,则每份的长度为 h n = h n h_n = \dfrac{h}{n} h n = n h ,而 2 n + 1 2n+1 2 n + 1 个分点为 x k = x 0 + k h n x_k = x_0 + kh_n x k = x 0 + k h n (k = 0 , ± 1 , ⋯ , ± n k = 0, \pm1, \cdots, \pm n k = 0 , ± 1 , ⋯ , ± n ),从而 x − n x_{-n} x − n 和 x n x_n x n 为区间 [ x 0 − h , , x 0 + h ] [x_0 - h,, x_0 + h] [ x 0 − h ,, x 0 + h ] 的两个端点,根据上述的想法,我们令欧拉折线的表达式为:
φ n ( x ) = y 0 + ∑ k = 0 s − 1 f ( x k , y k ) ( x k + 1 − x k ) + f ( x s , y s ) ( x − x s ) , ∀ x s < x ≤ x s + 1 ; ( 0 ≤ s ≤ n − 1 ) \varphi_n(x) = y_0 + \sum_{k=0}^{s-1} f(x_k, y_k)(x_{k+1} - x_k) + f(x_s, y_s)(x - x_s), \quad \forall\, x_s < x \leq x_{s+1}\,; (0 \leq s \leq n-1) φ n ( x ) = y 0 + k = 0 ∑ s − 1 f ( x k , y k ) ( x k + 1 − x k ) + f ( x s , y s ) ( x − x s ) , ∀ x s < x ≤ x s + 1 ; ( 0 ≤ s ≤ n − 1 )
这就给出了当 x 0 < x ≤ x 0 + h x_0 < x \leq x_0 + h x 0 < x ≤ x 0 + h 的欧拉折线的表达式,对于左边区域同理可得。
于是我们相当于需要证明欧拉折线 φ n ( x ) \varphi_n(x) φ n ( x ) 关于 n n n 在区间 ∣ x − x 0 ∣ ≤ h |x - x_0| \leq h ∣ x − x 0 ∣ ≤ h 上是收敛的,或者说有一个收敛的子列,更强一点,我们希望序列是一致收敛的,这样可以保证连续性,不仅如此我们最终要说明它是可导的,且满足 y ′ = f ( x , y ) y' = f(x,y) y ′ = f ( x , y ) 这一条件。为了解决这个收敛性问题,我们引入如下的 Ascoli 引理。
2.2 Ascoli 引理
我们先介绍两个概念,分别是函数列的一致有界和等度连续。
设在区间 I I I 上给定一个函数序列 f n ( x ) {f_n(x)} f n ( x ) ,如果存在常数 K > 0 K > 0 K > 0 ,使得 ∣ f n ( x ) ∣ < K |f_n(x)| < K ∣ f n ( x ) ∣ < K 对任意的 x ∈ I x \in I x ∈ I 和 n n n 都成立,则称上述的函数列在 I I I 上是一致有界 的。
如果对任意的正数 ε \varepsilon ε ,存在正数 δ \delta δ ,使得 ∀ , x 1 , x 2 ∈ I \forall, x_1, x_2 \in I ∀ , x 1 , x 2 ∈ I 且满足 ∣ x 1 − x 2 ∣ < δ |x_1 - x_2| < \delta ∣ x 1 − x 2 ∣ < δ 时,有 ∣ f n ( x 1 ) − f n ( x 2 ) ∣ |f_n(x_1) - f_n(x_2)| ∣ f n ( x 1 ) − f n ( x 2 ) ∣ 对于所有的 n n n 成立,则称上述的函数列在 I I I 上是等度(一致)连续 的。
并且我们可以证明如果一函数列在闭区间上收敛,且等度连续,那么它一致收敛,证明方法用定义加闭区间的无限开覆盖存在有限子覆盖加等度连续,根据这个我们容易证明如下的 Ascoli 引理 :
设 f n ( x ) {f_n(x)} f n ( x ) 在有限闭区间 I I I 上是一致有界和等度连续的,则可以选取它的一个子序列使得,该子序列在 I I I 上是一致收敛的。具体证明通过不断从某点收敛子列中找另一个收敛子列,再得到最终逐点收敛的对角序列,然后根据上一个结论易得。
2.3 Peano 存在性定理
接下来要证明的 Peano 存在性定理给出了当 f f f 仅连续时初值问题的解的存在性,在此之前,我们先来看两个引理,都是关于上述的欧拉折线的。
首先我们证明上述的欧拉序列在区间 ∣ x − x 0 ∣ ≤ h |x - x_0| \leq h ∣ x − x 0 ∣ ≤ h 至少有一个一致收敛的子序列,根据 Ascoli 引理,我们需要这一列欧拉折线满足一致有界和等度连续即可,一致有界是显然的都处于 D D D 的内部,且由于欧拉折线的各个直线段的斜率绝对值都小于等于 A A A ,从而容易证明每个欧拉折线的任何割线的斜率的绝对值也小于等于 A A A ,即有 ∣ φ n ( x 1 ) − φ n ( x 2 ) ∣ ≤ A ∣ x 1 − x 2 ∣ |\varphi_n(x_1) - \varphi_n(x_2)| \leq A|x_1 - x_2| ∣ φ n ( x 1 ) − φ n ( x 2 ) ∣ ≤ A ∣ x 1 − x 2 ∣ 对于任意 n n n 成立,从而满足等度连续,从而得证。
现在我们知道这一列欧拉折线存在子列一致收敛,收敛的函数有什么性质以及是否为初值问题的解还不确定,于是我们先证明下述的一个结论:欧拉折线 y = φ n ( x ) y = \varphi_n(x) y = φ n ( x ) 在 ∣ x − x 0 ∣ ≤ h |x - x_0| \leq h ∣ x − x 0 ∣ ≤ h 上满足
φ n ( x ) = y 0 + ∫ x 0 x f ( x , , φ n ( x ) ) d x + δ n ( x ) \varphi_n(x) = y_0 + \int_{x_0}^x f(x,, \varphi_n(x))\,dx + \delta_n(x) φ n ( x ) = y 0 + ∫ x 0 x f ( x ,, φ n ( x )) d x + δ n ( x )
其中 lim n → ∞ δ n ( x ) = 0 \lim_{n \to \infty} \delta_n(x) = 0 lim n → ∞ δ n ( x ) = 0 对任意 ∣ x − x 0 ∣ ≤ h |x - x_0| \leq h ∣ x − x 0 ∣ ≤ h ,事实上我们有
δ n ( x ) = ∑ i = 0 s − 1 ∫ x i x i + 1 [ f ( x i , y i ) − f ( x , , φ n ( x ) ) ] d x + ∫ x s x [ f ( x s , y s ) − f ( x , , φ n ( x ) ) ] d x \delta_n(x) = \sum_{i=0}^{s-1}\int_{x_i}^{x_{i+1}}\bigl[f(x_i, y_i) - f(x,, \varphi_n(x))\bigr]\,dx + \int_{x_s}^x \bigl[f(x_s, y_s) - f(x,, \varphi_n(x))\bigr]\,dx δ n ( x ) = i = 0 ∑ s − 1 ∫ x i x i + 1 [ f ( x i , y i ) − f ( x ,, φ n ( x )) ] d x + ∫ x s x [ f ( x s , y s ) − f ( x ,, φ n ( x )) ] d x
又因为在区间 [ x i , x i + 1 ] [x_i, x_{i+1}] [ x i , x i + 1 ] 上,我们有
∣ x − x i ∣ ≤ h n , ∣ φ n ( x ) − y i ∣ ≤ M ∣ x − x i ∣ ≤ M h n |x - x_i| \leq \frac{h}{n}, \quad |\varphi_n(x) - y_i| \leq M|x - x_i| \leq \frac{Mh}{n} ∣ x − x i ∣ ≤ n h , ∣ φ n ( x ) − y i ∣ ≤ M ∣ x − x i ∣ ≤ n M h
从而根据 f ( x , y ) f(x,y) f ( x , y ) 的连续性可推出 lim n → ∞ δ n ( x ) = 0 \lim_{n \to \infty} \delta_n(x) = 0 lim n → ∞ δ n ( x ) = 0 。
根据上述的两个结论,我们知道我们可以选取欧拉折线序列的一个子序列 φ n k ( x ) {\varphi_{n_k}(x)} φ n k ( x ) ,使其在区间 ∣ x − x 0 ∣ ≤ h |x - x_0| \leq h ∣ x − x 0 ∣ ≤ h 上一致收敛,从而极限函数 φ ( x ) = lim n → ∞ φ n k ( x ) \varphi(x) = \lim_{n \to \infty} \varphi_{n_k}(x) φ ( x ) = lim n → ∞ φ n k ( x ) 在该区间是连续的,又因为我们有
φ n k ( x ) = y 0 + ∫ x 0 x f ( x , φ n k ( x ) ) d x + δ n k ( x ) \varphi_{n_k}(x) = y_0 + \int_{x_0}^x f(x, \varphi_{n_k}(x))\,dx + \delta_{n_k}(x) φ n k ( x ) = y 0 + ∫ x 0 x f ( x , φ n k ( x )) d x + δ n k ( x )
再令 k → ∞ k \to \infty k → ∞ ,由一致收敛性和 f f f 的连续性,我们推出
lim k → ∞ ∫ x 0 x f ( x , φ n k ( x ) ) d x = ∫ x 0 x f ( x , φ ( x ) ) d x \lim_{k \to \infty}\int_{x_0}^x f(x, \varphi_{n_k}(x))\,dx = \int_{x_0}^x f(x, \varphi(x))\,dx k → ∞ lim ∫ x 0 x f ( x , φ n k ( x )) d x = ∫ x 0 x f ( x , φ ( x )) d x
从而有
φ ( x ) = y 0 + ∫ x 0 x f ( x , φ ( x ) ) d x \varphi(x) = y_0 + \int_{x_0}^x f(x, \varphi(x))\,dx φ ( x ) = y 0 + ∫ x 0 x f ( x , φ ( x )) d x
这说明 y = φ ( x ) y = \varphi(x) y = φ ( x ) 即为原初值问题在区间 ∣ x − x 0 ∣ ≤ h |x - x_0| \leq h ∣ x − x 0 ∣ ≤ h 的一个解。
2.4 解的延拓补充
前面我们讲述了对于初值问题唯一解的延拓,实际上一句话概括就是如果二元连续函数 f f f 在区域 G G G 上每个点都满足关于 y y y 的局部 Lipschitz 条件,那么区域 G ⊂ R 2 G \subset \mathbb{R}^2 G ⊂ R 2 上的初值问题的解曲线一定可以唯一延拓到 G G G 的边界,这个边界可以为 x , y x, y x , y 无穷远处。同样,该定理适用于 f f f 仅在 G G G 内连续的情况,只不过不是唯一延拓。
2.5 Schauder 不动点定理
设 X X X 为 Banach 空间,D ⊂ X D \subset X D ⊂ X 为闭凸子集,若映射 T : D → X T: D \to X T : D → X 为连续紧算子,且 T ( D ) ⊂ D T(D) \subset D T ( D ) ⊂ D ,则 T T T 存在不动点 x ∈ D x \in D x ∈ D 。
上述的连续算子和紧算子的概念如下:设 ( X , ∣ ⋅ ∣ X ) (X, |\cdot|_X) ( X , ∣ ⋅ ∣ X ) ,( Y , ∣ ⋅ ∣ Y ) (Y, |\cdot|_Y) ( Y , ∣ ⋅ ∣ Y ) 为赋范线性空间,
若映射 T : X → Y T: X \to Y T : X → Y 满足 X X X 中任意收敛于 x x x 的序列 x n {x_n} x n 有 T x n {T_{x_n}} T x n 收敛于 T x T_x T x ,则称 T T T 为连续算子
若映射 T : X → Y T: X \to Y T : X → Y 满足对于任意有界集 D ⊂ X D \subset X D ⊂ X ,有 T ( D ) T(D) T ( D ) 为 Y Y Y 中的相对紧集,则称 T T T 为紧算子
为了更好地理解这个定理的相关概念,我们引入拓扑学中度量空间中若干于紧性相关的概念,设 A A A 为度量空间 ( X , d ) (X, d) ( X , d ) 中的子集,则:
若 A A A 的任何开覆盖存在有限子覆盖,则 A A A 称为紧集
若 A A A 的闭包为紧集,则 A A A 称为相对紧集
若 A A A 中任意序列均有收敛子列且收敛于 A A A 中元,则 A A A 称为序列紧
若任给 ε > 0 \varepsilon > 0 ε > 0 ,存在有限的 B ( x 1 , ε ) , ⋯ , B ( x N , ε ) B(x_1,\varepsilon), \cdots, B(x_N,\varepsilon) B ( x 1 , ε ) , ⋯ , B ( x N , ε ) ,x i ∈ A x_i \in A x i ∈ A ,使得 A ⊂ ⋃ 1 ≤ i ≤ N B ( x i , ε ) A \subset \bigcup_{1 \leq i \leq N}B(x_i,\varepsilon) A ⊂ ⋃ 1 ≤ i ≤ N B ( x i , ε ) ,则称 A A A 完全有界
若 A A A 中任意 Cauchy 序列均收敛于 A A A 中元,则称 A A A 是完备的
有了上述定义我们有如下说法等价:
A A A 为紧集
A A A 为序列紧
A A A 完备且完全有界
以及如下说法等价:
A A A 为相对紧集
A A A 中任意序列均存在收敛子列
A A A 完全有界
根据 Schauder 不动点定理,我们也可以得出前述的解的存在性定理,这里省略证明。
3. 一阶显式常微分方程组与欧拉解
前面我们讨论了一阶显式常微分方程初值问题解的存在性和唯一性,实际其思想方法以及结论可以推广到一阶显式常微分方程组形如:
{ y 1 ′ = f 1 ( x , y 1 , ⋯ , y n ) ⋮ y n ′ = f n ( x , y 1 , ⋯ , y N ) \begin{cases} y_1' = f_1(x, y_1, \cdots, y_n) \\ \vdots \\ y_n' = f_n(x, y_1, \cdots, y_N) \end{cases} ⎩ ⎨ ⎧ y 1 ′ = f 1 ( x , y 1 , ⋯ , y n ) ⋮ y n ′ = f n ( x , y 1 , ⋯ , y N )
其中 n n n 个函数 f 1 , ⋯ , f n f_1, \cdots, f_n f 1 , ⋯ , f n 定义于 n + 1 n+1 n + 1 维空间 R n + 1 \mathbb{R}^{n+1} R n + 1 的子集 D D D 上,一个向量值函数 ( y 1 ( x ) , ⋯ , y n ( x ) ) T (y_1(x), \cdots, y_n(x))^T ( y 1 ( x ) , ⋯ , y n ( x ) ) T 称作上述方程组在区间 J J J 上的解是指每一 y i ( x ) y_i(x) y i ( x ) 在 J J J 上可微,且满足上述方程,自然地,我们要求任意 x ∈ J x \in J x ∈ J ,有 ( x , y 1 ( x ) , ⋯ , y n ( x ) ) ∈ D (x, y_1(x), \cdots, y_n(x)) \in D ( x , y 1 ( x ) , ⋯ , y n ( x )) ∈ D ,定义
y ( x ) = ( y 1 ( x ) ⋮ y n ( x ) ) , f ( x , y ) = ( f 1 ( x , y ) ⋮ f n ( x , y ) ) y(x) = \begin{pmatrix} y_1(x) \\\ \vdots \\\ y_n(x) \end{pmatrix}, \quad f(x,y) = \begin{pmatrix} f_1(x,y) \\ \vdots \ \\ f_n(x,y) \end{pmatrix} y ( x ) = y 1 ( x ) ⋮ y n ( x ) , f ( x , y ) = f 1 ( x , y ) ⋮ f n ( x , y )
则原方程组可简记为 y ′ ( x ) = f ( x , y ) y'(x) = f(x,y) y ′ ( x ) = f ( x , y ) ,相应初值问题记为 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 ,于一阶方程是一致的。
设 f : [ a , b ] × R n → R n f: [a,b] \times \mathbb{R}^n \to \mathbb{R}^n f : [ a , b ] × R n → R n ,我们考虑下面一阶常微分方程组的初值问题 { x ′ ( t ) = f ( t , x ( t ) ) x ( a ) = x 0 \begin{cases} x'(t) = f(t, x(t)) \\ x(a) = x_0 \end{cases} { x ′ ( t ) = f ( t , x ( t )) x ( a ) = x 0 ,同样我们采用欧拉折线法来研究这个方程组初值问题解的存在性,这里的 x ( t ) : [ a , b ] → R n x(t): [a,b] \to \mathbb{R}^n x ( t ) : [ a , b ] → R n 为一元向量值函数。
设 π : a = t 0 < t 1 < ⋯ < t n = b \pi: a = t_0 < t_1 < \cdots < t_n = b π : a = t 0 < t 1 < ⋯ < t n = b 为区间 [ a , b ] [a,b] [ a , b ] 的一个划分,同样由初值点 ( a , x 0 ) (a, x_0) ( a , x 0 ) 以及该点处的切向量 x ′ ( a ) = f ( a , x 0 ) x'(a) = f(a, x_0) x ′ ( a ) = f ( a , x 0 ) ,我们可以构造局部近似解,同前面所述一致,我们可以得到欧拉折线表达式为:
x π ( t ) = x 0 + ∑ k = 0 s − 1 f ( t k , x k ) ( t k + 1 − t k ) + f ( t s , x s ) ( t − t s ) , ∀ t s < t ≤ t s + 1 ; ( 0 ≤ s ≤ n − 1 ) x_\pi(t) = x_0 + \sum_{k=0}^{s-1} f(t_k, x_k)(t_{k+1} - t_k) + f(t_s, x_s)(t - t_s), \quad \forall\, t_s < t \leq t_{s+1}\,; (0 \leq s \leq n-1) x π ( t ) = x 0 + k = 0 ∑ s − 1 f ( t k , x k ) ( t k + 1 − t k ) + f ( t s , x s ) ( t − t s ) , ∀ t s < t ≤ t s + 1 ; ( 0 ≤ s ≤ n − 1 )
实际上之前的欧拉折线的划分也不一定要为等距划分,只需要限制分割的模就行。
我们这里先减弱一下 f f f 的要求,我们不要求它连续只要求其满足:存在 γ , c > 0 \gamma, c > 0 γ , c > 0 ,使得 ∣ f ( t , x ) ∣ ≤ γ ∣ x ∣ + c |f(t,x)| \leq \gamma|x| + c ∣ f ( t , x ) ∣ ≤ γ ∣ x ∣ + c ,∀ , ( t , x ) ∈ [ a , b ] × R n \forall, (t,x) \in [a,b] \times \mathbb{R}^n ∀ , ( t , x ) ∈ [ a , b ] × R n ,我们同样可以推出 Ascoli 定理成立的两个条件,从而证明一列分割的模收敛于零的欧拉折线,其必然存在子列一致收敛,设其收敛于 x ( t ) x(t) x ( t ) ,我们称这个 x ( t ) x(t) x ( t ) 为原初值问题的一个欧拉解 ,注意欧拉解未必是原方程的解,因为我们还没有验证其是否满足积分方程,只是一个近似解,且我们有其任意欧拉解 x ( t ) x(t) x ( t ) 满足
∣ x ( t ) − x ( a ) ∣ ≤ ( t − a ) e γ ( t − a ) + γ ∣ x ( a ) ∣ , ∀ t ∈ [ a , b ] |x(t) - x(a)| \leq (t-a)e^{\gamma(t-a)} + \gamma|x(a)|, \quad \forall\, t \in [a,b] ∣ x ( t ) − x ( a ) ∣ ≤ ( t − a ) e γ ( t − a ) + γ ∣ x ( a ) ∣ , ∀ t ∈ [ a , b ]
当然如果 f f f 是连续的,则欧拉解为原初值问题的解。
而如果设 D ⊂ R × R n D \subset \mathbb{R} \times \mathbb{R}^n D ⊂ R × R n 为开区域,f : D → R n f: D \to \mathbb{R}^n f : D → R n 满足全局的关于 y y y 的 Lipschitz 条件,则 D D D 上每一点的初值问题都存在唯一解,且不会出现前述的情形二有限时间爆炸,从而解可以向左和向右延拓到 D D D 的边界。对于上述方程组问题的讨论中我们没有去严格证明,虽然大同小异但还是有一定区别,这里不赘述。
4. 比较定理
4.1 第一比较定理
如果 f ( x , y ) f(x,y) f ( x , y ) 与 F ( x , y ) F(x,y) F ( x , y ) 在平面区域 G G G 内都连续且满足 f ( x , y ) < F ( x , y ) f(x,y) < F(x,y) f ( x , y ) < F ( x , y ) 在 G G G 上恒成立,那么对于这两个二元函数形成的初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 与 { y ′ = F ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = F(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = F ( x , y ) y ( x 0 ) = y 0 ,我们由 Peano 存在性定理易知都是有解的,不妨设 [ a , b ] [a,b] [ a , b ] 属于二者的公共解存在区间,令 ϕ 1 , ϕ 2 \phi_1, \phi_2 ϕ 1 , ϕ 2 分别为上述的初值问题在 [ a , b ] [a,b] [ a , b ] 上的解。
那么可以证明有
{ ϕ 1 ( x ) < ϕ 2 ( x ) , x 0 < x < b ϕ 1 ( x ) > ϕ 2 ( x ) , a < x < x 0 \begin{cases} \phi_1(x) < \phi_2(x), & x_0 < x < b \\ \phi_1(x) > \phi_2(x), & a < x < x_0 \end{cases} { ϕ 1 ( x ) < ϕ 2 ( x ) , ϕ 1 ( x ) > ϕ 2 ( x ) , x 0 < x < b a < x < x 0
我们称之为第一比较定理 。需要注意的是这里都不能取等。
实际上我们可以知道只要 ϕ 1 \phi_1 ϕ 1 和 ϕ 2 \phi_2 ϕ 2 的函数值相等,就必然会有此点的斜率前者小于后者,而由于初始点函数值已经相等,从而若右侧其他点函数值相等,必然会有前者斜率大于后者,从而矛盾得出上述结论,从几何意义来说就是斜率小的曲线向右不可能从斜率大的曲线的下方穿越到上方。
4.2 极大解与极小解
考虑初值问题 E : { y ′ = f ( x , y ) y ( x 0 ) = y 0 E: \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} E : { y ′ = f ( x , y ) y ( x 0 ) = y 0 ,f f f 在矩形区域 D = [ x 0 − a , , x 0 + a ] × [ y 0 − b , , y 0 + b ] D = [x_0 - a,, x_0 + a] \times [y_0 - b,, y_0 + b] D = [ x 0 − a ,, x 0 + a ] × [ y 0 − b ,, y 0 + b ] 上连续,令 A = max x , y ∈ D ∣ f ( x , y ) ∣ A = \max_{x,y \in D}|f(x,y)| A = max x , y ∈ D ∣ f ( x , y ) ∣ ,h = min { a , b A } h = \min\{a, \frac{b}{A}\} h = min { a , A b } ,如果在区间 ∣ x − x 0 ∣ ≤ h |x - x_0| \leq h ∣ x − x 0 ∣ ≤ h 上初值问题有两个解 y = z ( x ) y = z(x) y = z ( x ) 和 y = w ( x ) y = w(x) y = w ( x ) 使得初值问题的任何解 y = y ( x ) y = y(x) y = y ( x ) 都满足 w ( x ) ≤ y ( x ) ≤ z ( x ) w(x) \leq y(x) \leq z(x) w ( x ) ≤ y ( x ) ≤ z ( x ) ,则称 y = w ( x ) y = w(x) y = w ( x ) 和 y = z ( x ) y = z(x) y = z ( x ) 分别为初值问题的极小解 和极大解 。容易证明大区间的极大解和极小解一定是由小区间的极大解和极小解延拓来的,也就是说不可能出现极大解或极小解在唯一延拓后不是极大解或极小解的情况。
关于极大解和极小解的存在性,我们下面给出证明其一定存在正数 σ < h \sigma < h σ < h ,使得在区间 ∣ x − x 0 ∣ < σ |x - x_0| < \sigma ∣ x − x 0 ∣ < σ 上,初值问题一定有极大解和极小解。
我们考虑与上述初值问题相联系的初值问题 E m : { y ′ = f ( x , y ) + ε m y ( x 0 ) = y 0 E_m: \begin{cases} y' = f(x,y) + \varepsilon_m \\ y(x_0) = y_0 \end{cases} E m : { y ′ = f ( x , y ) + ε m y ( x 0 ) = y 0 ,其中 ε m > 0 \varepsilon_m > 0 ε m > 0 且单调递减趋于零,由 Peano 存在性定理,存在 h > h m > 0 h > h_m > 0 h > h m > 0 ,使得 E m E_m E m 在区间 ∣ x − x 0 ∣ < h m |x - x_0| < h_m ∣ x − x 0 ∣ < h m 上有解 y = ϕ m ( x ) y = \phi_m(x) y = ϕ m ( x ) ,且易知 lim m → ∞ h m = h \lim_{m \to \infty} h_m = h lim m → ∞ h m = h ,并且
ϕ m ( x ) = y 0 + ∫ x 0 x [ f ( x , ϕ m ( x ) ) + ε m ] d x \phi_m(x) = y_0 + \int_{x_0}^x \bigl[f(x, \phi_m(x)) + \varepsilon_m\bigr]\,dx ϕ m ( x ) = y 0 + ∫ x 0 x [ f ( x , ϕ m ( x )) + ε m ] d x
从而可知存在 σ > 0 \sigma > 0 σ > 0 ,使得初值问题 E E E 和 E m E_m E m 在区间 I : ∣ x − x 0 ∣ < σ I: |x - x_0| < \sigma I : ∣ x − x 0 ∣ < σ 上解存在。
又因为任意 m m m 和任意 x ∈ I x \in I x ∈ I ,我们有 ∣ ϕ m ( x ) − y 0 ∣ ≤ b |\phi_m(x) - y_0| \leq b ∣ ϕ m ( x ) − y 0 ∣ ≤ b ,且 ∣ ϕ m ( x 1 ) − ϕ m ( x 2 ) ∣ ≤ ( A + ε 1 ) ∣ x 1 − x 2 ∣ |\phi_m(x_1) - \phi_m(x_2)| \leq (A + \varepsilon_1)|x_1 - x_2| ∣ ϕ m ( x 1 ) − ϕ m ( x 2 ) ∣ ≤ ( A + ε 1 ) ∣ x 1 − x 2 ∣ ,从而 ϕ m ( x ) \phi_m(x) ϕ m ( x ) 在 I I I 上一致有界与等度连续,利用 Ascoli 引理可知其存在一致收敛的子列,设其收敛于 φ ( x ) \varphi(x) φ ( x ) ,从而对积分方程取极限可知 φ ( x ) \varphi(x) φ ( x ) 为 E E E 在 I I I 上的一个解,由第一比较定理易知 y = φ ( x ) y = \varphi(x) y = φ ( x ) 是 I I I 上的右行最大解和左行最小解。
对上述过程用 − ε m -\varepsilon_m − ε m 替换 ε m \varepsilon_m ε m ,可证 I I I 上 E E E 有左行极大解和右行极小解。我们知道初值问题的所有解在交点处必然相切,于是我们可以根据上述得到的两个解拼接出 I I I 上的极大解和极小解,事实上我们在进行上述分析的时候可以不用限制定义域,类似于解的延拓对于每个初值问题我们将其解延拓至边界即可,尽管每个解的定义域可能不一致,但如此我们也能得到矩形区域 D D D 乃至平面区域 G G G 上的极小解和极大解。
从而我们知道上述初值问题的解唯一等价于它的极大解和极小解是恒同的。
4.3 第二比较定理
我们容易得到如下的第二比较定理 :设 f ( x , y ) f(x,y) f ( x , y ) 与 F ( x , y ) F(x,y) F ( x , y ) 在 G G G 上连续且满足 f ( x , y ) ≤ F ( x , y ) f(x,y) \leq F(x,y) f ( x , y ) ≤ F ( x , y ) 在 G G G 上恒成立,不妨设 [ a , b ] [a,b] [ a , b ] 属于二者的公共解存在区间,令 y 1 ( x ) , y 2 ( x ) y_1(x), y_2(x) y 1 ( x ) , y 2 ( x ) 分别为初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 与 { y ′ = F ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = F(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = F ( x , y ) y ( x 0 ) = y 0 在 [ a , b ] [a,b] [ a , b ] 上的解,且 y 1 ( x ) y_1(x) y 1 ( x ) 为右行极小解和左行极大解(或者 y 2 ( x ) y_2(x) y 2 ( x ) 为右行极大解和左行极小解),则有
{ y 1 ( x ) ≤ y 2 ( x ) , x 0 < x < b y 1 ( x ) ≥ y 2 ( x ) , a < x < x 0 \begin{cases} y_1(x) \leq y_2(x), & x_0 < x < b \\ y_1(x) \geq y_2(x), & a < x < x_0 \end{cases} { y 1 ( x ) ≤ y 2 ( x ) , y 1 ( x ) ≥ y 2 ( x ) , x 0 < x < b a < x < x 0
这里是对第一比较定理无法取等的补充。
对于上述定理我们常常运用下面的推论,即常常把比较对象定为唯一可解的常微分方程的初值问题:假设 f f f 在 G G G 上连续,易知初值问题 { y ′ = f ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = f(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = f ( x , y ) y ( x 0 ) = y 0 一定有解,如果初值问题 { y ′ = F 1 ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = F_1(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = F 1 ( x , y ) y ( x 0 ) = y 0 与 { y ′ = F 2 ( x , y ) y ( x 0 ) = y 0 \begin{cases} y' = F_2(x,y) \\ y(x_0) = y_0 \end{cases} { y ′ = F 2 ( x , y ) y ( x 0 ) = y 0 在 G G G 上有唯一解,分别为 y 1 , y 2 y_1, y_2 y 1 , y 2 ,且 F 2 ( x , y ) ≤ f ( x , y ) ≤ F 1 ( x , y ) F_2(x,y) \leq f(x,y) \leq F_1(x,y) F 2 ( x , y ) ≤ f ( x , y ) ≤ F 1 ( x , y ) 在 G G G 上恒成立,那么对于 f f f 的初值问题的任意解 y ( x ) y(x) y ( x ) ,以右侧为例,y 2 ( x ) ≤ y ( x ) ≤ y 1 ( x ) y_2(x) \leq y(x) \leq y_1(x) y 2 ( x ) ≤ y ( x ) ≤ y 1 ( x ) 在两两之间的公共定义域内恒成立,也就是说如果 y 1 y_1 y 1 在 G G G 上没有爆炸,那么任意 y y y 都不会爆炸,如果 y 2 y_2 y 2 在 G G G 上爆炸,那么任意 y y y 在 G G G 上都爆炸。
实际上我们也可以这么理解,对于矩形区域 D = [ x 0 , a ] × R D = [x_0, a] \times \mathbb{R} D = [ x 0 , a ] × R 上解的延拓如果没有延拓到 D D D 的右边界,说明解在有限时间爆炸,可以认为爆炸点右侧的区域的函数值都是无穷大,从而 y 2 ( x ) ≤ y ( x ) ≤ y 1 ( x ) y_2(x) \leq y(x) \leq y_1(x) y 2 ( x ) ≤ y ( x ) ≤ y 1 ( x ) 可以认为在 [ x 0 , a ] [x_0, a] [ x 0 , a ] 上成立。
5. 奇解与包络
前面我们都在探究形如 y ′ = f ( x , y ) y' = f(x,y) y ′ = f ( x , y ) 的一阶显式常微分方程的初值问题的解,下面我们直接探究形如 F ( x , y , y ′ ) = 0 F(x, y, y') = 0 F ( x , y , y ′ ) = 0 的一阶常微分方程的解,首先我们需要了解几个概念,对于一阶常微分方程的通解形如 y = φ ( x , C ) y = \varphi(x, C) y = φ ( x , C ) ,其中 C C C 为任意常数时 y y y 均为原方程的解,而特解是满足微分方程的具体解,可以是通解中给定固定的常数 C C C 得到的解,也可以是通解之外的解,如奇解或一些其他的孤立解。
那么什么是奇解呢?首先它是不属于通解族的特解,其次该解曲线上的每一点,在其任何邻域内原微分方程都有一个不同于该特解的解与之在该点相切,则称该特解为奇解 。
根据上述定义,我们知道微分方程如果有奇解,那么奇解上的每一点都至少有两个解,也就是解不唯一,如果 F ( x , y , p ) F(x,y,p) F ( x , y , p ) 是 C 1 C^1 C 1 的,且 y = φ ( x ) y = \varphi(x) y = φ ( x ) 为微分方程的奇解,则必然有 F ( x , φ ( x ) , φ ′ ( x ) ) = 0 F(x, \varphi(x), \varphi'(x)) = 0 F ( x , φ ( x ) , φ ′ ( x )) = 0 ,且 F p ( x , φ ( x ) , φ ′ ( x ) ) = 0 F_p(x, \varphi(x), \varphi'(x)) = 0 F p ( x , φ ( x ) , φ ′ ( x )) = 0 ,第一个等式是显然的,对于第二个等式,假设某点不成立,那么该点由隐映射定理 F ( x , y , p ) = 0 F(x,y,p) = 0 F ( x , y , p ) = 0 局部可以写为 p = f ( x , y ) p = f(x,y) p = f ( x , y ) 参数曲面形式,且
f y ( x , y ) = − F y ( x , y , f ( x , y ) ) F p ( x , y , f ( x , y ) ) f_y(x,y) = -\frac{F_y(x, y, f(x,y))}{F_p(x, y, f(x,y))} f y ( x , y ) = − F p ( x , y , f ( x , y )) F y ( x , y , f ( x , y ))
实际上 f f f 也是 C 1 C^1 C 1 的,从而 f y ( x , y ) f_y(x,y) f y ( x , y ) 是连续的,从而可知在 ( x 0 , φ ( x 0 ) ) (x_0, \varphi(x_0)) ( x 0 , φ ( x 0 )) 的某个邻域中解是唯一的,这与该点的两条不同解曲线矛盾,故原等式成立。
值得注意的是,如果反过来,满足上述两个等式的函数则不一定是奇解。
事实上有时候我们可以根据上述的两个方程通过消去 p p p ,得到一个函数 y = φ ( x ) y = \varphi(x) y = φ ( x ) ,当然此时这个函数不一定是解,但是如果是解,即满足 F ( x , φ ( x ) , φ ′ ( x ) ) = 0 F(x, \varphi(x), \varphi'(x)) = 0 F ( x , φ ( x ) , φ ′ ( x )) = 0 ,那么它还要满足什么条件才是奇解。值得注意的是,此时 φ ( x ) \varphi(x) φ ( x ) 不一定满足 F p ( x , φ ( x ) , φ ′ ( x ) ) = 0 F_p(x, \varphi(x), \varphi'(x)) = 0 F p ( x , φ ( x ) , φ ′ ( x )) = 0 ,例如考虑微分方程 y = 2 x + y ′ − 1 3 ( y ′ ) 3 y = 2x + y' - \dfrac{1}{3}(y')^3 y = 2 x + y ′ − 3 1 ( y ′ ) 3 ,所以我们首先要加约束条件 F p ( x , φ ( x ) , φ ′ ( x ) ) = 0 F_p(x, \varphi(x), \varphi'(x)) = 0 F p ( x , φ ( x ) , φ ′ ( x )) = 0 ,其次我们需要要求 F F F 是 C 2 C^2 C 2 的,且 F y ( x , φ ( x ) , φ ′ ( x ) ) ≠ 0 F_y(x, \varphi(x), \varphi'(x)) \neq 0 F y ( x , φ ( x ) , φ ′ ( x )) = 0 ,以及 F p p ( x , φ ( x ) , φ ′ ( x ) ) ≠ 0 F_{pp}(x, \varphi(x), \varphi'(x)) \neq 0 F pp ( x , φ ( x ) , φ ′ ( x )) = 0 ,如此 φ ( x ) \varphi(x) φ ( x ) 一定是奇解。
下面我们先来讲述一般曲线族的包络,再探究包络与奇解之间的关系。
设单参数 C C C 的曲线族 K ( C ) : V ( x , y , C ) = 0 K(C): V(x,y,C) = 0 K ( C ) : V ( x , y , C ) = 0 ,如果在平面上有一条连续可导的曲线 Γ \Gamma Γ ,对于该曲线上每一点,在 K ( C ) K(C) K ( C ) 中均有一条在该点某邻域不同于 Γ \Gamma Γ 的曲线通过该点且在该点与 Γ \Gamma Γ 相切,则称 Γ \Gamma Γ 为 K ( C ) K(C) K ( C ) 的一支包络 。
并不是所有的曲线族都有包络,如同心圆族。
对于 F ( x , y , y ′ ) = 0 F(x,y,y') = 0 F ( x , y , y ′ ) = 0 ,我们假设其有通积分 U ( x , y , C ) = 0 U(x,y,C) = 0 U ( x , y , C ) = 0 ,且该曲线族有包络 Γ : y = φ ( x ) \Gamma: y = \varphi(x) Γ : y = φ ( x ) ,则我们可以证明 y = φ ( x ) y = \varphi(x) y = φ ( x ) 为微分方程的奇解。实际上这只需要证明 y = φ ( x ) y = \varphi(x) y = φ ( x ) 为微分方程的解即可,由于该曲线上每一点及其斜率与某个通解的一致,易知满足微分方程。
如果曲线族 K ( C ) : V ( x , y , C ) = 0 K(C): V(x,y,C) = 0 K ( C ) : V ( x , y , C ) = 0 中 V V V 是 C 1 C^1 C 1 的,且存在一支包络 Γ \Gamma Γ ,则包络上的每一点 ( x 0 , y 0 ) (x_0, y_0) ( x 0 , y 0 ) ,存在 C 0 C_0 C 0 ,使得 V ( x 0 , y 0 , C 0 ) = 0 V(x_0, y_0, C_0) = 0 V ( x 0 , y 0 , C 0 ) = 0 ,V C ( x 0 , y 0 , C 0 ) = 0 V_C(x_0, y_0, C_0) = 0 V C ( x 0 , y 0 , C 0 ) = 0 同时成立。
同样我们可以通过上述两个方程消去 C C C 得到一个函数 y = φ ( x ) y = \varphi(x) y = φ ( x ) ,这个函数也不一定是曲线族的包络,但如果满足 V V V 是 C 2 C^2 C 2 的,且
∣ V x V y V C x V C y ∣ ≠ 0 , V C C ≠ 0 \begin{vmatrix} V_x & V_y \\ V_{Cx} & V_{Cy} \end{vmatrix} \neq 0, \quad V_{CC} \neq 0 V x V C x V y V C y = 0 , V C C = 0
则 φ ( x ) \varphi(x) φ ( x ) 为曲线族的包络。