この章の目標
初期値問題を積分方程式に書き直し、リプシッツ条件のもとでの解の存在と一意性を、逐次近似と縮小写像の 2 通りの方法で証明できる
グロンウォールの不等式を証明し、一意性や解の誤差評価に使える
最大存在区間と解の延長定理を理解し、解の爆発を説明できる
解が初期値とパラメータに連続に依存することを証明できる
ペアノの存在定理をアスコリ–アルツェラの定理から証明できる
前提 :第1章 、微分積分学 第6章 、集合と位相 第5章 、第7章
第1章では、y ′ = y 2 y' = y^2 y ′ = y 2 の解が有限の時刻で爆発すること(例 1.10)、トリチェリの法則 h ′ = − k h h' = -k\sqrt{h} h ′ = − k h では解が一意でないこと(例 1.11)、y ′ = x 2 + y 2 y' = x^2 + y^2 y ′ = x 2 + y 2 は求積法で解けないこと(注意 1.6)を見た。本章では、解の式を書かずに次の問いに答える。解は存在するか。一意か。どこまで延長できるか。初期値を変えると解はどう変わるか。
以下 D ⊂ R × R n D \subset \mathbb{R} \times \mathbb{R}^n D ⊂ R × R n を開集合、f : D → R n f\colon D \to \mathbb{R}^n f : D → R n を連続写像とし、R n \mathbb{R}^n R n のノルム ∥ ⋅ ∥ \lVert \cdot \rVert ∥ ⋅ ∥ を一つ固定する(どのノルムも同値である)。n n n 階の方程式は 1 階の連立方程式に書き直せる(第1章 1.1 節)ので、x ′ = f ( t , x ) x' = f(t, x) x ′ = f ( t , x ) の形だけを考えれば十分である。ベクトル値関数の積分は成分ごとに定義し、∥ ∫ a b g ( s ) d s ∥ ≤ ∫ a b ∥ g ( s ) ∥ d s \left\lVert \int_a^b g(s)\ ds \right\rVert \leq \int_a^b \lVert g(s) \rVert\ ds ∫ a b g ( s ) d s ≤ ∫ a b ∥ g ( s )∥ d s (a ≤ b a \leq b a ≤ b )を使う(リーマン和に三角不等式を適用して極限をとればよい)。
2.1 初期値問題と積分方程式
( t 0 , x 0 ) ∈ D (t_0, x_0) \in D ( t 0 , x 0 ) ∈ D とし、初期値問題 x ′ = f ( t , x ) x' = f(t, x) x ′ = f ( t , x ) , x ( t 0 ) = x 0 x(t_0) = x_0 x ( t 0 ) = x 0 を考える。t 0 t_0 t 0 を含む区間 J J J 上の解 とは、微分可能な x : J → R n x\colon J \to \mathbb{R}^n x : J → R n で、グラフ { ( t , x ( t ) ) } \lbrace (t, x(t)) \rbrace {( t , x ( t ))} が D D D に含まれ、方程式と初期条件をみたすものである。
補題 2.1 x : J → R n x\colon J \to \mathbb{R}^n x : J → R n を連続でグラフが D D D に含まれる写像とする。x x x が初期値問題の解であるための必要十分条件は、すべての t ∈ J t \in J t ∈ J で
x ( t ) = x 0 + ∫ t 0 t f ( s , x ( s ) ) d s x(t) = x_0 + \int_{t_0}^t f(s, x(s))\,ds x ( t ) = x 0 + ∫ t 0 t f ( s , x ( s )) d s
が成り立つことである。このとき x x x は C 1 C^1 C 1 級である。
証明. 解ならば x ′ ( t ) = f ( t , x ( t ) ) x'(t) = f(t, x(t)) x ′ ( t ) = f ( t , x ( t )) は連続なので x x x は C 1 C^1 C 1 級で、微分積分学の基本定理から等式が従う。逆に等式が成り立てば、被積分関数は連続なので右辺は微分可能で x ′ = f ( t , x ) x' = f(t, x) x ′ = f ( t , x ) 、また x ( t 0 ) = x 0 x(t_0) = x_0 x ( t 0 ) = x 0 。□ \square □
この書き直しにより、微分が、一様収束と交換できる積分に置き換わる。
2.2 リプシッツ条件
定義 2.2 (リプシッツ条件, Lipschitz condition)定数 L ≥ 0 L \geq 0 L ≥ 0 があって、( t , x ) , ( t , y ) ∈ D (t, x), (t, y) \in D ( t , x ) , ( t , y ) ∈ D ならば常に
∥ f ( t , x ) − f ( t , y ) ∥ ≤ L ∥ x − y ∥ \lVert f(t, x) - f(t, y) \rVert \leq L\lVert x - y \rVert ∥ f ( t , x ) − f ( t , y )∥ ≤ L ∥ x − y ∥
が成り立つとき、f f f は D D D 上で x x x についてリプシッツ条件 をみたすといい、L L L をリプシッツ定数という。D D D の各点が、その上で f f f がリプシッツ条件をみたす近傍をもつとき、f f f は x x x について局所リプシッツ (locally Lipschitz) であるという。
命題 2.3 f f f と偏導関数 ∂ f i / ∂ x j \partial f_i/\partial x_j ∂ f i / ∂ x j が D D D 上で連続なら、f f f は x x x について局所リプシッツである。
証明. ( t 1 , x 1 ) ∈ D (t_1, x_1) \in D ( t 1 , x 1 ) ∈ D に対し、D D D は開集合なので C = [ t 1 − r , t 1 + r ] × B ‾ ( x 1 , r ) ⊂ D C = [t_1 - r, t_1 + r] \times \overline{B}(x_1, r) \subset D C = [ t 1 − r , t 1 + r ] × B ( x 1 , r ) ⊂ D となる r > 0 r > 0 r > 0 がある。C C C はコンパクトなので、ヤコビ行列 D x f D_xf D x f の作用素ノルムは C C C 上で最大値 L L L をとる。( t , x ) , ( t , y ) ∈ C (t, x), (t, y) \in C ( t , x ) , ( t , y ) ∈ C なら線分 y + θ ( x − y ) y + \theta(x - y) y + θ ( x − y ) (0 ≤ θ ≤ 1 0 \leq \theta \leq 1 0 ≤ θ ≤ 1 )は凸集合 B ‾ ( x 1 , r ) \overline{B}(x_1, r) B ( x 1 , r ) に含まれ、θ ↦ f ( t , y + θ ( x − y ) ) \theta \mapsto f(t, y + \theta(x - y)) θ ↦ f ( t , y + θ ( x − y )) に微分積分学の基本定理を使うと
∥ f ( t , x ) − f ( t , y ) ∥ = ∥ ∫ 0 1 D x f ( t , y + θ ( x − y ) ) ( x − y ) d θ ∥ ≤ L ∥ x − y ∥ \lVert f(t, x) - f(t, y) \rVert = \left\lVert \int_0^1 D_xf(t, y + \theta(x - y))(x - y)\,d\theta \right\rVert \leq L\lVert x - y \rVert ∥ f ( t , x ) − f ( t , y )∥ = ∫ 0 1 D x f ( t , y + θ ( x − y )) ( x − y ) d θ ≤ L ∥ x − y ∥
C C C は ( t 1 , x 1 ) (t_1, x_1) ( t 1 , x 1 ) の近傍なので主張が従う。□ \square □
例 2.4 (a) f ( x ) = ∣ x ∣ f(x) = \lvert x \rvert f ( x ) = ∣ x ∣ は微分可能でないがリプシッツ条件をみたす(L = 1 L = 1 L = 1 )。(b) f ( x ) = x 2 f(x) = x^2 f ( x ) = x 2 は局所リプシッツだが、∣ x 2 − y 2 ∣ = ∣ x + y ∣ ∣ x − y ∣ \lvert x^2 - y^2 \rvert = \lvert x + y \rvert\lvert x - y \rvert ∣ x 2 − y 2 ∣ = ∣ x + y ∣ ∣ x − y ∣ なので R \mathbb{R} R 全体ではリプシッツでない。(c) f ( x ) = ∣ x ∣ f(x) = \sqrt{\lvert x \rvert} f ( x ) = ∣ x ∣ は 0 の近傍でリプシッツでない。∣ f ( x ) − f ( 0 ) ∣ / ∣ x ∣ = ∣ x ∣ − 1 / 2 \lvert f(x) - f(0) \rvert/\lvert x \rvert = \lvert x \rvert^{-1/2} ∣ f ( x ) − f ( 0 )∣ / ∣ x ∣ = ∣ x ∣ − 1/2 が x → 0 x \to 0 x → 0 で有界でないからである。
補題 2.5 f f f が x x x について局所リプシッツなら、D D D の任意のコンパクト部分集合 K K K 上でリプシッツ条件をみたす。すなわち ( t , x ) , ( t , y ) ∈ K (t, x), (t, y) \in K ( t , x ) , ( t , y ) ∈ K なら ∥ f ( t , x ) − f ( t , y ) ∥ ≤ L K ∥ x − y ∥ \lVert f(t, x) - f(t, y) \rVert \leq L_K\lVert x - y \rVert ∥ f ( t , x ) − f ( t , y )∥ ≤ L K ∥ x − y ∥ となる定数 L K L_K L K がある。
証明. そうでないとすると、( t k , x k ) , ( t k , y k ) ∈ K (t_k, x_k), (t_k, y_k) \in K ( t k , x k ) , ( t k , y k ) ∈ K で ∥ f ( t k , x k ) − f ( t k , y k ) ∥ > k ∥ x k − y k ∥ \lVert f(t_k, x_k) - f(t_k, y_k) \rVert > k\lVert x_k - y_k \rVert ∥ f ( t k , x k ) − f ( t k , y k )∥ > k ∥ x k − y k ∥ となるものがある。M = max K ∥ f ∥ M = \max_K \lVert f \rVert M = max K ∥ f ∥ とすると ∥ x k − y k ∥ < 2 M / k → 0 \lVert x_k - y_k \rVert < 2M/k \to 0 ∥ x k − y k ∥ < 2 M / k → 0 。部分列をとって ( t k , x k ) → ( t ∗ , x ∗ ) ∈ K (t_k, x_k) \to (t^{\ast}, x^{\ast}) \in K ( t k , x k ) → ( t ∗ , x ∗ ) ∈ K としてよく、このとき y k → x ∗ y_k \to x^{\ast} y k → x ∗ 。( t ∗ , x ∗ ) (t^{\ast}, x^{\ast}) ( t ∗ , x ∗ ) の近傍で f f f がリプシッツ定数 L ′ L' L ′ をもつとすると、十分大きい k k k で k ∥ x k − y k ∥ < L ′ ∥ x k − y k ∥ k\lVert x_k - y_k \rVert < L'\lVert x_k - y_k \rVert k ∥ x k − y k ∥ < L ′ ∥ x k − y k ∥ となり(x k ≠ y k x_k \neq y_k x k = y k に注意)、k < L ′ k < L' k < L ′ に矛盾する。□ \square □
2.3 ピカール–リンデレフの定理
本節では次を仮定する。
仮定 (A) :a , b > 0 a, b > 0 a , b > 0 に対し R = { ( t , x ) ∣ ∣ t − t 0 ∣ ≤ a , ∥ x − x 0 ∥ ≤ b } ⊂ D R = \lbrace (t, x) \mid \lvert t - t_0 \rvert \leq a,\ \lVert x - x_0 \rVert \leq b \rbrace \subset D R = {( t , x ) ∣ ∣ t − t 0 ∣ ≤ a , ∥ x − x 0 ∥ ≤ b } ⊂ D とし、M = max R ∥ f ∥ M = \max_R \lVert f \rVert M = max R ∥ f ∥ 、h = min ( a , b / M ) h = \min(a, b/M) h = min ( a , b / M ) (M = 0 M = 0 M = 0 なら h = a h = a h = a )、I = [ t 0 − h , t 0 + h ] I = [t_0 - h, t_0 + h] I = [ t 0 − h , t 0 + h ] とおく。
補題 2.6 仮定 (A) のもとで、t 0 t_0 t 0 を含む区間 J ⊂ I J \subset I J ⊂ I 上の解 x x x は ∥ x ( t ) − x 0 ∥ ≤ M ∣ t − t 0 ∣ \lVert x(t) - x_0 \rVert \leq M\lvert t - t_0 \rvert ∥ x ( t ) − x 0 ∥ ≤ M ∣ t − t 0 ∣ をみたす。特にそのグラフは R R R に含まれる。
証明. t ≥ t 0 t \geq t_0 t ≥ t 0 の側を示す(逆側も同様)。S = { τ ∈ J ∣ τ ≥ t 0 , [ t 0 , τ ] 上で ∥ x ( t ) − x 0 ∥ ≤ b } S = \lbrace \tau \in J \mid \tau \geq t_0,\ [t_0, \tau] \text{ 上で } \lVert x(t) - x_0 \rVert \leq b \rbrace S = { τ ∈ J ∣ τ ≥ t 0 , [ t 0 , τ ] 上で ∥ x ( t ) − x 0 ∥ ≤ b } とおく。S S S は t 0 t_0 t 0 を含む区間で、連続性から閉である。τ ∈ S \tau \in S τ ∈ S なら [ t 0 , τ ] [t_0, \tau] [ t 0 , τ ] 上でグラフは R R R に含まれるので、補題 2.1 より ∥ x ( t ) − x 0 ∥ ≤ ∫ t 0 t ∥ f ( s , x ( s ) ) ∥ d s ≤ M ( t − t 0 ) \lVert x(t) - x_0 \rVert \leq \int_{t_0}^t \lVert f(s, x(s)) \rVert\ ds \leq M(t - t_0) ∥ x ( t ) − x 0 ∥ ≤ ∫ t 0 t ∥ f ( s , x ( s ))∥ d s ≤ M ( t − t 0 ) 。S S S が J ∩ [ t 0 , ∞ ) J \cap [t_0, \infty) J ∩ [ t 0 , ∞ ) 全体でないとすると、τ ∗ = max S \tau^{\ast} = \max S τ ∗ = max S の右にも J J J の点があるので τ ∗ < t 0 + h \tau^{\ast} < t_0 + h τ ∗ < t 0 + h であり、∥ x ( τ ∗ ) − x 0 ∥ ≤ M ( τ ∗ − t 0 ) < b \lVert x(\tau^{\ast}) - x_0 \rVert \leq M(\tau^{\ast} - t_0) < b ∥ x ( τ ∗ ) − x 0 ∥ ≤ M ( τ ∗ − t 0 ) < b (M > 0 M > 0 M > 0 なら M ( τ ∗ − t 0 ) < M h ≤ b M(\tau^{\ast} - t_0) < Mh \leq b M ( τ ∗ − t 0 ) < M h ≤ b 、M = 0 M = 0 M = 0 なら左辺は 0)。連続性から τ ∗ \tau^{\ast} τ ∗ の少し右でも ∥ x ( t ) − x 0 ∥ < b \lVert x(t) - x_0 \rVert < b ∥ x ( t ) − x 0 ∥ < b となり、τ ∗ \tau^{\ast} τ ∗ の最大性に矛盾する。□ \square □
定理 2.7 (ピカール–リンデレフの定理, Picard–Lindelöf theorem)仮定 (A) に加えて、f f f が R R R 上で x x x についてリプシッツ条件(定数 L L L )をみたすとする。このとき初期値問題 x ′ = f ( t , x ) x' = f(t, x) x ′ = f ( t , x ) , x ( t 0 ) = x 0 x(t_0) = x_0 x ( t 0 ) = x 0 は I I I 上の解 φ \varphi φ をもつ。また t 0 t_0 t 0 を含む区間 J ⊂ I J \subset I J ⊂ I 上の解は、J J J 上で φ \varphi φ と一致する。
証明 1(逐次近似). φ 0 ( t ) = x 0 \varphi_0(t) = x_0 φ 0 ( t ) = x 0 とし、t ∈ I t \in I t ∈ I に対して
φ k + 1 ( t ) = x 0 + ∫ t 0 t f ( s , φ k ( s ) ) d s \varphi_{k+1}(t) = x_0 + \int_{t_0}^t f(s, \varphi_k(s))\,ds φ k + 1 ( t ) = x 0 + ∫ t 0 t f ( s , φ k ( s )) d s
で逐次近似 (successive approximation) φ k \varphi_k φ k を定める。
(i) φ k \varphi_k φ k は連続でグラフは R R R に含まれる。帰納法による:φ k \varphi_k φ k がそうなら φ k + 1 \varphi_{k+1} φ k + 1 は定義できて連続で、∥ φ k + 1 ( t ) − x 0 ∥ ≤ M ∣ t − t 0 ∣ ≤ M h ≤ b \lVert \varphi_{k+1}(t) - x_0 \rVert \leq M\lvert t - t_0 \rvert \leq Mh \leq b ∥ φ k + 1 ( t ) − x 0 ∥ ≤ M ∣ t − t 0 ∣ ≤ M h ≤ b 。
(ii) ∥ φ k + 1 ( t ) − φ k ( t ) ∥ ≤ M L k ∣ t − t 0 ∣ k + 1 ( k + 1 ) ! \lVert \varphi_{k+1}(t) - \varphi_k(t) \rVert \leq \dfrac{ML^k\lvert t - t_0 \rvert^{k+1}}{(k+1)!} ∥ φ k + 1 ( t ) − φ k ( t )∥ ≤ ( k + 1 )! M L k ∣ t − t 0 ∣ k + 1 。k = 0 k = 0 k = 0 は (i) の評価そのものである。k k k で成り立てば、t ≥ t 0 t \geq t_0 t ≥ t 0 として
∥ φ k + 2 ( t ) − φ k + 1 ( t ) ∥ ≤ ∫ t 0 t L ∥ φ k + 1 ( s ) − φ k ( s ) ∥ d s ≤ ∫ t 0 t M L k + 1 ( s − t 0 ) k + 1 ( k + 1 ) ! d s = M L k + 1 ( t − t 0 ) k + 2 ( k + 2 ) ! \lVert \varphi_{k+2}(t) - \varphi_{k+1}(t) \rVert \leq \int_{t_0}^t L\lVert \varphi_{k+1}(s) - \varphi_k(s) \rVert\,ds \leq \int_{t_0}^t \frac{ML^{k+1}(s - t_0)^{k+1}}{(k+1)!}\,ds = \frac{ML^{k+1}(t - t_0)^{k+2}}{(k+2)!} ∥ φ k + 2 ( t ) − φ k + 1 ( t )∥ ≤ ∫ t 0 t L ∥ φ k + 1 ( s ) − φ k ( s )∥ d s ≤ ∫ t 0 t ( k + 1 )! M L k + 1 ( s − t 0 ) k + 1 d s = ( k + 2 )! M L k + 1 ( t − t 0 ) k + 2
(t ≤ t 0 t \leq t_0 t ≤ t 0 も同様)。
(iii) (ii) より ∑ k sup I ∥ φ k + 1 − φ k ∥ ≤ ∑ k M L k h k + 1 ( k + 1 ) ! < ∞ \sum_k \sup_I \lVert \varphi_{k+1} - \varphi_k \rVert \leq \sum_k \frac{ML^kh^{k+1}}{(k+1)!} < \infty ∑ k sup I ∥ φ k + 1 − φ k ∥ ≤ ∑ k ( k + 1 )! M L k h k + 1 < ∞ なので、ワイエルシュトラスの M 判定法により φ k = φ 0 + ∑ j < k ( φ j + 1 − φ j ) \varphi_k = \varphi_0 + \sum_{j < k}(\varphi_{j+1} - \varphi_j) φ k = φ 0 + ∑ j < k ( φ j + 1 − φ j ) は I I I 上で一様収束する。極限 φ \varphi φ は連続で、R R R は閉集合なのでグラフは R R R に含まれる。
(iv) ∥ f ( s , φ k ( s ) ) − f ( s , φ ( s ) ) ∥ ≤ L ∥ φ k ( s ) − φ ( s ) ∥ → 0 \lVert f(s, \varphi_k(s)) - f(s, \varphi(s)) \rVert \leq L\lVert \varphi_k(s) - \varphi(s) \rVert \to 0 ∥ f ( s , φ k ( s )) − f ( s , φ ( s ))∥ ≤ L ∥ φ k ( s ) − φ ( s )∥ → 0 (一様)なので、漸化式で k → ∞ k \to \infty k → ∞ とすると φ ( t ) = x 0 + ∫ t 0 t f ( s , φ ( s ) ) d s \varphi(t) = x_0 + \int_{t_0}^t f(s, \varphi(s))\ ds φ ( t ) = x 0 + ∫ t 0 t f ( s , φ ( s )) d s 。補題 2.1 より φ \varphi φ は解である。
(v) 一意性:ψ \psi ψ を J J J 上の解とすると、補題 2.6 よりグラフは R R R に含まれる。t 0 t_0 t 0 を含むコンパクト区間 J ′ ⊂ J J' \subset J J ′ ⊂ J 上で K = sup J ′ ∥ ψ − φ ∥ K = \sup_{J'}\lVert \psi - \varphi \rVert K = sup J ′ ∥ ψ − φ ∥ とおくと、∥ ψ ( t ) − φ ( t ) ∥ ≤ ∣ ∫ t 0 t L ∥ ψ − φ ∥ d s ∣ \lVert \psi(t) - \varphi(t) \rVert \leq \left\lvert \int_{t_0}^t L\lVert \psi - \varphi \rVert\ ds \right\rvert ∥ ψ ( t ) − φ ( t )∥ ≤ ∫ t 0 t L ∥ ψ − φ ∥ d s の繰り返しにより ∥ ψ ( t ) − φ ( t ) ∥ ≤ K L k ∣ t − t 0 ∣ k / k ! → 0 \lVert \psi(t) - \varphi(t) \rVert \leq KL^k\lvert t - t_0 \rvert^k/k! \to 0 ∥ ψ ( t ) − φ ( t )∥ ≤ K L k ∣ t − t 0 ∣ k / k ! → 0 。よって J J J 上で ψ = φ \psi = \varphi ψ = φ 。□ \square □
例 2.8 x ′ = x x' = x x ′ = x , x ( 0 ) = 1 x(0) = 1 x ( 0 ) = 1 の逐次近似は φ k ( t ) = ∑ j = 0 k t j / j ! \varphi_k(t) = \sum_{j=0}^k t^j/j! φ k ( t ) = ∑ j = 0 k t j / j ! であり、e t e^t e t のテイラー展開の部分和が現れる。
証明 2(縮小写像の原理). L > 0 L > 0 L > 0 としてよい。C ( I , R n ) C(I, \mathbb{R}^n) C ( I , R n ) の閉部分集合 X = { φ ∣ ∥ φ ( t ) − x 0 ∥ ≤ b ( t ∈ I ) } X = \lbrace \varphi \mid \lVert \varphi(t) - x_0 \rVert \leq b\ (t \in I) \rbrace X = { φ ∣ ∥ φ ( t ) − x 0 ∥ ≤ b ( t ∈ I )} に、重みつきの距離
d ( φ , ψ ) = sup t ∈ I e − 2 L ∣ t − t 0 ∣ ∥ φ ( t ) − ψ ( t ) ∥ d(\varphi, \psi) = \sup_{t \in I} e^{-2L\lvert t - t_0 \rvert}\lVert \varphi(t) - \psi(t) \rVert d ( φ , ψ ) = t ∈ I sup e − 2 L ∣ t − t 0 ∣ ∥ φ ( t ) − ψ ( t )∥
を入れる。e − 2 L h d ∞ ≤ d ≤ d ∞ e^{-2Lh}d_\infty \leq d \leq d_\infty e − 2 L h d ∞ ≤ d ≤ d ∞ (d ∞ d_\infty d ∞ は一様距離)なので d d d は一様距離と同値であり、C ( I , R n ) C(I, \mathbb{R}^n) C ( I , R n ) は一様距離で完備(集合と位相 定理 7.6 )だから ( X , d ) (X, d) ( X , d ) も完備である。( T φ ) ( t ) = x 0 + ∫ t 0 t f ( s , φ ( s ) ) d s (T\varphi)(t) = x_0 + \int_{t_0}^t f(s, \varphi(s))\ ds ( T φ ) ( t ) = x 0 + ∫ t 0 t f ( s , φ ( s )) d s とおくと、証明 1 の (i) と同様に T ( X ) ⊂ X T(X) \subset X T ( X ) ⊂ X 。t ≥ t 0 t \geq t_0 t ≥ t 0 では
∥ T φ ( t ) − T ψ ( t ) ∥ ≤ ∫ t 0 t L e 2 L ( s − t 0 ) d ( φ , ψ ) d s ≤ 1 2 e 2 L ( t − t 0 ) d ( φ , ψ ) \lVert T\varphi(t) - T\psi(t) \rVert \leq \int_{t_0}^t L e^{2L(s - t_0)} d(\varphi, \psi)\,ds \leq \frac{1}{2}e^{2L(t - t_0)}d(\varphi, \psi) ∥ T φ ( t ) − T ψ ( t )∥ ≤ ∫ t 0 t L e 2 L ( s − t 0 ) d ( φ , ψ ) d s ≤ 2 1 e 2 L ( t − t 0 ) d ( φ , ψ )
であり、t ≤ t 0 t \leq t_0 t ≤ t 0 でも同様なので d ( T φ , T ψ ) ≤ 1 2 d ( φ , ψ ) d(T\varphi, T\psi) \leq \frac{1}{2}d(\varphi, \psi) d ( T φ , T ψ ) ≤ 2 1 d ( φ , ψ ) 。縮小写像の原理(集合と位相 定理 7.14 )より T T T は X X X 内にただ一つの不動点 φ \varphi φ をもち、補題 2.1 より φ \varphi φ は解である。J ⊂ I J \subset I J ⊂ I 上の解 ψ \psi ψ は補題 2.6 より J J J 上の X X X に対応する空間に属し、J J J 上で同じ議論をすれば不動点の一意性から ψ = φ ∣ J \psi = \varphi\vert_J ψ = φ ∣ J 。□ \square □
2.4 グロンウォールの不等式
一意性や誤差の評価に繰り返し使う不等式を用意する。
定理 2.10 (グロンウォールの不等式, Gronwall's inequality)u , β u, \beta u , β を [ t 0 , t 1 ] [t_0, t_1] [ t 0 , t 1 ] 上の連続関数、β ≥ 0 \beta \geq 0 β ≥ 0 、α ∈ R \alpha \in \mathbb{R} α ∈ R とする。すべての t ∈ [ t 0 , t 1 ] t \in [t_0, t_1] t ∈ [ t 0 , t 1 ] で
u ( t ) ≤ α + ∫ t 0 t β ( s ) u ( s ) d s u(t) \leq \alpha + \int_{t_0}^t \beta(s)u(s)\,ds u ( t ) ≤ α + ∫ t 0 t β ( s ) u ( s ) d s
ならば、u ( t ) ≤ α exp ( ∫ t 0 t β ( s ) d s ) u(t) \leq \alpha\exp\left(\int_{t_0}^t \beta(s)\ ds\right) u ( t ) ≤ α exp ( ∫ t 0 t β ( s ) d s ) が成り立つ。
証明. 右辺を Φ ( t ) = α + ∫ t 0 t β u d s \Phi(t) = \alpha + \int_{t_0}^t \beta u\ ds Φ ( t ) = α + ∫ t 0 t β u d s 、B ( t ) = ∫ t 0 t β d s B(t) = \int_{t_0}^t \beta\ ds B ( t ) = ∫ t 0 t β d s とおく。β ≥ 0 \beta \geq 0 β ≥ 0 と u ≤ Φ u \leq \Phi u ≤ Φ より Φ ′ = β u ≤ β Φ \Phi' = \beta u \leq \beta\Phi Φ ′ = β u ≤ β Φ なので、( Φ e − B ) ′ = ( Φ ′ − β Φ ) e − B ≤ 0 (\Phi e^{-B})' = (\Phi' - \beta\Phi)e^{-B} \leq 0 ( Φ e − B ) ′ = ( Φ ′ − β Φ ) e − B ≤ 0 。よって Φ ( t ) e − B ( t ) ≤ Φ ( t 0 ) = α \Phi(t)e^{-B(t)} \leq \Phi(t_0) = \alpha Φ ( t ) e − B ( t ) ≤ Φ ( t 0 ) = α で、u ( t ) ≤ Φ ( t ) ≤ α e B ( t ) u(t) \leq \Phi(t) \leq \alpha e^{B(t)} u ( t ) ≤ Φ ( t ) ≤ α e B ( t ) 。□ \square □
t ≤ t 0 t \leq t_0 t ≤ t 0 の側(u ( t ) ≤ α + ∫ t t 0 β u d s u(t) \leq \alpha + \int_t^{t_0}\beta u\ ds u ( t ) ≤ α + ∫ t t 0 β u d s )でも、t t t を 2 t 0 − t 2t_0 - t 2 t 0 − t に置き換えれば同じ結論 u ( t ) ≤ α exp ∫ t t 0 β d s u(t) \leq \alpha\exp\int_t^{t_0}\beta\ ds u ( t ) ≤ α exp ∫ t t 0 β d s を得る。
系 2.11 (解の誤差評価)区間 J ∋ t 0 J \ni t_0 J ∋ t 0 上で、x x x は x ′ = f ( t , x ) x' = f(t, x) x ′ = f ( t , x ) の解、y y y は y ′ = g ( t , y ) y' = g(t, y) y ′ = g ( t , y ) の解とする。f f f は両者のグラフを含む集合上で x x x について定数 L > 0 L > 0 L > 0 のリプシッツ条件をみたし、y y y のグラフ上で ∥ f ( t , y ) − g ( t , y ) ∥ ≤ ε \lVert f(t, y) - g(t, y) \rVert \leq \varepsilon ∥ f ( t , y ) − g ( t , y )∥ ≤ ε とする。このとき t ∈ J t \in J t ∈ J で
∥ x ( t ) − y ( t ) ∥ ≤ ∥ x ( t 0 ) − y ( t 0 ) ∥ e L ∣ t − t 0 ∣ + ε L ( e L ∣ t − t 0 ∣ − 1 ) \lVert x(t) - y(t) \rVert \leq \lVert x(t_0) - y(t_0) \rVert e^{L\lvert t - t_0 \rvert} + \frac{\varepsilon}{L}\left(e^{L\lvert t - t_0 \rvert} - 1\right) ∥ x ( t ) − y ( t )∥ ≤ ∥ x ( t 0 ) − y ( t 0 )∥ e L ∣ t − t 0 ∣ + L ε ( e L ∣ t − t 0 ∣ − 1 )
特に g = f g = f g = f で x ( t 0 ) = y ( t 0 ) x(t_0) = y(t_0) x ( t 0 ) = y ( t 0 ) なら x = y x = y x = y である。
証明. t ≥ t 0 t \geq t_0 t ≥ t 0 とする。x ( t ) − y ( t ) = x ( t 0 ) − y ( t 0 ) + ∫ t 0 t ( f ( s , x ) − f ( s , y ) + f ( s , y ) − g ( s , y ) ) d s x(t) - y(t) = x(t_0) - y(t_0) + \int_{t_0}^t\bigl(f(s, x) - f(s, y) + f(s, y) - g(s, y)\bigr)ds x ( t ) − y ( t ) = x ( t 0 ) − y ( t 0 ) + ∫ t 0 t ( f ( s , x ) − f ( s , y ) + f ( s , y ) − g ( s , y ) ) d s より、u = ∥ x − y ∥ u = \lVert x - y \rVert u = ∥ x − y ∥ は u ( t ) ≤ u ( t 0 ) + ∫ t 0 t ( L u + ε ) d s u(t) \leq u(t_0) + \int_{t_0}^t (Lu + \varepsilon)\ ds u ( t ) ≤ u ( t 0 ) + ∫ t 0 t ( Lu + ε ) d s をみたす。w = u + ε / L w = u + \varepsilon/L w = u + ε / L とおくと w ( t ) ≤ w ( t 0 ) + ∫ t 0 t L w d s w(t) \leq w(t_0) + \int_{t_0}^t Lw\ ds w ( t ) ≤ w ( t 0 ) + ∫ t 0 t L w d s なので、定理 2.10 より w ( t ) ≤ w ( t 0 ) e L ( t − t 0 ) w(t) \leq w(t_0)e^{L(t - t_0)} w ( t ) ≤ w ( t 0 ) e L ( t − t 0 ) 。これが主張である。t ≤ t 0 t \leq t_0 t ≤ t 0 も同様。□ \square □
誤差は時間とともに高々指数関数的にしか増えないが、L L L が大きい系では初期値のわずかな誤差が急速に拡大しうる。
2.5 一意性が崩れる例
例 2.12 x ′ = ∣ x ∣ x' = \sqrt{\lvert x \rvert} x ′ = ∣ x ∣ , x ( 0 ) = 0 x(0) = 0 x ( 0 ) = 0 を考える。x ≡ 0 x \equiv 0 x ≡ 0 は解である。さらに 0 ≤ c ≤ ∞ 0 \leq c \leq \infty 0 ≤ c ≤ ∞ , 0 ≤ d ≤ ∞ 0 \leq d \leq \infty 0 ≤ d ≤ ∞ に対し
x ( t ) = { ( t − c ) 2 / 4 ( t > c ) 0 ( − d ≤ t ≤ c ) − ( t + d ) 2 / 4 ( t < − d ) x(t) = \begin{cases} (t - c)^2/4 & (t > c) \\ 0 & (-d \leq t \leq c) \\ -(t + d)^2/4 & (t < -d) \end{cases} x ( t ) = ⎩ ⎨ ⎧ ( t − c ) 2 /4 0 − ( t + d ) 2 /4 ( t > c ) ( − d ≤ t ≤ c ) ( t < − d )
も解である。実際 t > c t > c t > c では x ′ = ( t − c ) / 2 = x x' = (t - c)/2 = \sqrt{x} x ′ = ( t − c ) /2 = x 、t < − d t < -d t < − d では x ′ = − ( t + d ) / 2 = ∣ x ∣ x' = -(t + d)/2 = \sqrt{\lvert x \rvert} x ′ = − ( t + d ) /2 = ∣ x ∣ であり、つなぎ目の微分係数は 0 である。解は無数にあり、「いつ動き出すか」が決まらない。右辺は x = 0 x = 0 x = 0 でリプシッツ連続でない(例 2.4 (c))。x 0 > 0 x_0 > 0 x 0 > 0 から出発した解 ( t / 2 + x 0 ) 2 (t/2 + \sqrt{x_0})^2 ( t /2 + x 0 ) 2 も、過去にさかのぼって x = 0 x = 0 x = 0 に達すると枝分かれする。トリチェリの法則(例 1.11)と同じ現象である。
2.6 解の延長と最大存在区間
本節と次節では次を仮定する。
仮定 (B) :f : D → R n f\colon D \to \mathbb{R}^n f : D → R n は連続で、x x x について局所リプシッツである。
補題 2.14 (大域的な一意性)仮定 (B) のもとで、x 1 x_1 x 1 (区間 J 1 J_1 J 1 上)と x 2 x_2 x 2 (区間 J 2 J_2 J 2 上)が同じ初期値問題の解なら、J 1 ∩ J 2 J_1 \cap J_2 J 1 ∩ J 2 上で x 1 = x 2 x_1 = x_2 x 1 = x 2 。
証明. S = { t ∈ J 1 ∩ J 2 ∣ x 1 ( t ) = x 2 ( t ) } S = \lbrace t \in J_1 \cap J_2 \mid x_1(t) = x_2(t) \rbrace S = { t ∈ J 1 ∩ J 2 ∣ x 1 ( t ) = x 2 ( t )} は t 0 t_0 t 0 を含み、連続性から J 1 ∩ J 2 J_1 \cap J_2 J 1 ∩ J 2 で閉である。t 1 ∈ S t_1 \in S t 1 ∈ S なら、( t 1 , x 1 ( t 1 ) ) (t_1, x_1(t_1)) ( t 1 , x 1 ( t 1 )) の近傍の長方形で定理 2.7 を適用すると、t 1 t_1 t 1 の近くで x 1 x_1 x 1 と x 2 x_2 x 2 は一致するので、S S S は開でもある。区間 J 1 ∩ J 2 J_1 \cap J_2 J 1 ∩ J 2 は連結なので S = J 1 ∩ J 2 S = J_1 \cap J_2 S = J 1 ∩ J 2 。□ \square □
定理 2.15 (最大解)仮定 (B) のもとで、各 ( t 0 , x 0 ) ∈ D (t_0, x_0) \in D ( t 0 , x 0 ) ∈ D に対し、開区間 ( ω − , ω + ) ∋ t 0 (\omega_-, \omega_+) \ni t_0 ( ω − , ω + ) ∋ t 0 (− ∞ ≤ ω − < ω + ≤ ∞ -\infty \leq \omega_- < \omega_+ \leq \infty − ∞ ≤ ω − < ω + ≤ ∞ )上の解 x x x で、初期値問題の任意の解が x x x の制限になっているものがただ一つ存在する。
この x x x を最大解 (maximal solution)、( ω − , ω + ) (\omega_-, \omega_+) ( ω − , ω + ) を最大存在区間 (maximal interval of existence) という。
証明. 初期値問題のすべての解の定義域の和集合を J ∗ J^{\ast} J ∗ とすると、J ∗ J^{\ast} J ∗ は t 0 t_0 t 0 を含む区間である。t ∈ J ∗ t \in J^{\ast} t ∈ J ∗ で定義された解の t t t での値を x ( t ) x(t) x ( t ) とおくと、補題 2.14 によりこれは解の選び方によらない。x x x は各点の近くで何らかの解と一致するので J ∗ J^{\ast} J ∗ 上の解であり、どの解も x x x の制限である。J ∗ J^{\ast} J ∗ が右端 ω + \omega_+ ω + を含むとすると、( ω + , x ( ω + ) ) ∈ D (\omega_+, x(\omega_+)) \in D ( ω + , x ( ω + )) ∈ D で定理 2.7 を使って [ ω + , ω + + h ] [\omega_+, \omega_+ + h] [ ω + , ω + + h ] 上の解 y y y を作り、x x x とつなげた関数は積分方程式をみたす(補題 2.1)ので、J ∗ J^{\ast} J ∗ より長い区間上の解になり矛盾する。左端も同様で、J ∗ J^{\ast} J ∗ は開区間である。□ \square □
定理 2.16 (延長定理)仮定 (B) のもとで、x x x を最大解とする。ω + < ∞ \omega_+ < \infty ω + < ∞ ならば、任意のコンパクト集合 K ⊂ D K \subset D K ⊂ D に対し、ある t K < ω + t_K < \omega_+ t K < ω + があって、t K < t < ω + t_K < t < \omega_+ t K < t < ω + なら ( t , x ( t ) ) ∉ K (t, x(t)) \notin K ( t , x ( t )) ∈ / K となる。ω − \omega_- ω − についても同様である。
つまり、最大解のグラフは有限時刻で終わるなら、D D D のどのコンパクト集合からも最終的に出ていく。
証明. D D D は開、K K K はコンパクトなので、r > 0 r > 0 r > 0 を小さくとれば、K K K の点 ( s , y ) (s, y) ( s , y ) から ∣ t − s ∣ ≤ r \lvert t - s \rvert \leq r ∣ t − s ∣ ≤ r , ∥ x − y ∥ ≤ r \lVert x - y \rVert \leq r ∥ x − y ∥ ≤ r の範囲にある点全体 K r K_r K r は D D D に含まれるコンパクト集合になる(K r K_r K r は有界閉集合であり、K K K と D D D の補集合との距離は正だから)。M = max K r ∥ f ∥ M = \max_{K_r}\lVert f \rVert M = max K r ∥ f ∥ 、L L L を K r K_r K r 上のリプシッツ定数(補題 2.5)、h = min ( r , r / M ) h = \min(r, r/M) h = min ( r , r / M ) (M = 0 M = 0 M = 0 なら h = r h = r h = r )とする。各 ( s , y ) ∈ K (s, y) \in K ( s , y ) ∈ K で長方形 [ s − r , s + r ] × B ‾ ( y , r ) [s - r, s + r] \times \overline{B}(y, r) [ s − r , s + r ] × B ( y , r ) は K r K_r K r に含まれるので、定理 2.7 より ( s , y ) (s, y) ( s , y ) を通る解が [ s − h , s + h ] [s - h, s + h] [ s − h , s + h ] 上に存在する。h h h は ( s , y ) (s, y) ( s , y ) によらないことが要点である。いま t k → ω + t_k \to \omega_+ t k → ω + かつ ( t k , x ( t k ) ) ∈ K (t_k, x(t_k)) \in K ( t k , x ( t k )) ∈ K となる列があったとし、ω + − t k < h \omega_+ - t_k < h ω + − t k < h となる k k k をとる。( t k , x ( t k ) ) (t_k, x(t_k)) ( t k , x ( t k )) を通る [ t k − h , t k + h ] [t_k - h, t_k + h] [ t k − h , t k + h ] 上の解は補題 2.14 により x x x と共通部分で一致するので、x x x は t k + h > ω + t_k + h > \omega_+ t k + h > ω + まで延長でき、最大性に矛盾する。□ \square □
系 2.17 (爆発)D = J × R n D = J \times \mathbb{R}^n D = J × R n (J J J は開区間)で ω + < sup J \omega_+ < \sup J ω + < sup J ならば、t → ω + t \to \omega_+ t → ω + で ∥ x ( t ) ∥ → ∞ \lVert x(t) \rVert \to \infty ∥ x ( t )∥ → ∞ 。
証明. 任意の ρ > 0 \rho > 0 ρ > 0 に対し K = [ t 0 , ω + ] × B ‾ ( 0 , ρ ) K = [t_0, \omega_+] \times \overline{B}(0, \rho) K = [ t 0 , ω + ] × B ( 0 , ρ ) は D D D のコンパクト部分集合なので、定理 2.16 より ω + \omega_+ ω + の十分近くで ∥ x ( t ) ∥ > ρ \lVert x(t) \rVert > \rho ∥ x ( t )∥ > ρ 。□ \square □
例 2.18 (a) x ′ = x 2 x' = x^2 x ′ = x 2 , x ( 0 ) = x 0 > 0 x(0) = x_0 > 0 x ( 0 ) = x 0 > 0 の最大解は x = x 0 / ( 1 − x 0 t ) x = x_0/(1 - x_0t) x = x 0 / ( 1 − x 0 t ) で ω + = 1 / x 0 \omega_+ = 1/x_0 ω + = 1/ x 0 。初期値が大きいほど早く爆発する。(b) D = R × ( 0 , ∞ ) D = \mathbb{R} \times (0, \infty) D = R × ( 0 , ∞ ) 上の x ′ = − 1 / ( 2 x ) x' = -1/(2x) x ′ = − 1/ ( 2 x ) , x ( 0 ) = 1 x(0) = 1 x ( 0 ) = 1 の最大解は x = 1 − t x = \sqrt{1 - t} x = 1 − t で ω + = 1 \omega_+ = 1 ω + = 1 。解は有界だが、D D D の境界 x = 0 x = 0 x = 0 に近づくことで D D D のコンパクト集合から出ていく。(c) x ′ = − x 3 x' = -x^3 x ′ = − x 3 の解は ( x 2 ) ′ = − 2 x 4 ≤ 0 (x^2)' = -2x^4 \leq 0 ( x 2 ) ′ = − 2 x 4 ≤ 0 より t ≥ t 0 t \geq t_0 t ≥ t 0 で有界なので、系 2.17 より ω + = ∞ \omega_+ = \infty ω + = ∞ 。解を求めずに得られるこのような評価をアプリオリ評価 (a priori estimate) という。
系 2.19 (線形増大条件)D = J × R n D = J \times \mathbb{R}^n D = J × R n とし、仮定 (B) に加えて J J J 上の連続関数 α , β ≥ 0 \alpha, \beta \geq 0 α , β ≥ 0 により ∥ f ( t , x ) ∥ ≤ α ( t ) + β ( t ) ∥ x ∥ \lVert f(t, x) \rVert \leq \alpha(t) + \beta(t)\lVert x \rVert ∥ f ( t , x )∥ ≤ α ( t ) + β ( t ) ∥ x ∥ が成り立つとする。このとき最大解は J J J 全体で定義される。特に、連続な係数をもつ線形系 x ′ = A ( t ) x + b ( t ) x' = A(t)x + b(t) x ′ = A ( t ) x + b ( t ) の解は J J J 全体に延長できる。
証明. ω + < sup J \omega_+ < \sup J ω + < sup J とし、[ t 0 , ω + ] [t_0, \omega_+] [ t 0 , ω + ] 上で α ≤ A \alpha \leq A α ≤ A , β ≤ B \beta \leq B β ≤ B とする。t ∈ [ t 0 , ω + ) t \in [t_0, \omega_+) t ∈ [ t 0 , ω + ) で ∥ x ( t ) ∥ ≤ ∥ x 0 ∥ + A ( ω + − t 0 ) + ∫ t 0 t B ∥ x ( s ) ∥ d s \lVert x(t) \rVert \leq \lVert x_0 \rVert + A(\omega_+ - t_0) + \int_{t_0}^t B\lVert x(s) \rVert\ ds ∥ x ( t )∥ ≤ ∥ x 0 ∥ + A ( ω + − t 0 ) + ∫ t 0 t B ∥ x ( s )∥ d s なので、定理 2.10 より ∥ x ( t ) ∥ ≤ ( ∥ x 0 ∥ + A ( ω + − t 0 ) ) e B ( ω + − t 0 ) \lVert x(t) \rVert \leq (\lVert x_0 \rVert + A(\omega_+ - t_0))e^{B(\omega_+ - t_0)} ∥ x ( t )∥ ≤ (∥ x 0 ∥ + A ( ω + − t 0 )) e B ( ω + − t 0 ) と有界になり、系 2.17 に矛盾する。ω − \omega_- ω − も同様。線形系では ∥ A ( t ) x + b ( t ) ∥ ≤ ∥ b ( t ) ∥ + ∥ A ( t ) ∥ ∥ x ∥ \lVert A(t)x + b(t) \rVert \leq \lVert b(t) \rVert + \lVert A(t) \rVert\lVert x \rVert ∥ A ( t ) x + b ( t )∥ ≤ ∥ b ( t )∥ + ∥ A ( t )∥ ∥ x ∥ である。□ \square □
2.7 初期値とパラメータへの連続依存性
現実の初期値や係数(パラメータ)は誤差を含む。それらを少し変えたとき解が少ししか変わらないことは、微分方程式をモデルとして使うための最低限の条件である。
定理 2.20 (連続依存性)Λ ⊂ R m \Lambda \subset \mathbb{R}^m Λ ⊂ R m を開集合、f : D × Λ → R n f\colon D \times \Lambda \to \mathbb{R}^n f : D × Λ → R n を連続とし、D × Λ D \times \Lambda D × Λ の各点の近傍 U U U と定数 L L L で、( t , x , λ ) , ( t , y , λ ) ∈ U (t, x, \lambda), (t, y, \lambda) \in U ( t , x , λ ) , ( t , y , λ ) ∈ U なら ∥ f ( t , x , λ ) − f ( t , y , λ ) ∥ ≤ L ∥ x − y ∥ \lVert f(t, x, \lambda) - f(t, y, \lambda) \rVert \leq L\lVert x - y \rVert ∥ f ( t , x , λ ) − f ( t , y , λ )∥ ≤ L ∥ x − y ∥ となるものがあるとする。x ′ = f ( t , x , λ ) x' = f(t, x, \lambda) x ′ = f ( t , x , λ ) , x ( t 0 ) = ξ x(t_0) = \xi x ( t 0 ) = ξ の最大解を x ( t ; ξ , λ ) x(t; \xi, \lambda) x ( t ; ξ , λ ) と書く。x ∗ ( t ) = x ( t ; ξ 0 , λ 0 ) x^{\ast}(t) = x(t; \xi_0, \lambda_0) x ∗ ( t ) = x ( t ; ξ 0 , λ 0 ) がコンパクト区間 [ a , b ] ∋ t 0 [a, b] \ni t_0 [ a , b ] ∋ t 0 上で定義されているとする。このとき任意の ε > 0 \varepsilon > 0 ε > 0 に対し δ > 0 \delta > 0 δ > 0 があって、∥ ξ − ξ 0 ∥ < δ \lVert \xi - \xi_0 \rVert < \delta ∥ ξ − ξ 0 ∥ < δ , ∥ λ − λ 0 ∥ < δ \lVert \lambda - \lambda_0 \rVert < \delta ∥ λ − λ 0 ∥ < δ なら x ( t ; ξ , λ ) x(t; \xi, \lambda) x ( t ; ξ , λ ) は [ a , b ] [a, b] [ a , b ] 上で定義され
sup a ≤ t ≤ b ∥ x ( t ; ξ , λ ) − x ∗ ( t ) ∥ < ε \sup_{a \leq t \leq b}\lVert x(t; \xi, \lambda) - x^{\ast}(t) \rVert < \varepsilon a ≤ t ≤ b sup ∥ x ( t ; ξ , λ ) − x ∗ ( t )∥ < ε
証明. (i) Γ = { ( t , x ∗ ( t ) , λ 0 ) ∣ a ≤ t ≤ b } \Gamma = \lbrace (t, x^{\ast}(t), \lambda_0) \mid a \leq t \leq b \rbrace Γ = {( t , x ∗ ( t ) , λ 0 ) ∣ a ≤ t ≤ b } はコンパクトなので、0 < ρ < ε 0 < \rho < \varepsilon 0 < ρ < ε を小さくとれば
U = { ( t , x , λ ) ∣ a ≤ t ≤ b , ∥ x − x ∗ ( t ) ∥ ≤ ρ , ∥ λ − λ 0 ∥ ≤ ρ } U = \lbrace (t, x, \lambda) \mid a \leq t \leq b,\ \lVert x - x^{\ast}(t) \rVert \leq \rho,\ \lVert \lambda - \lambda_0 \rVert \leq \rho \rbrace U = {( t , x , λ ) ∣ a ≤ t ≤ b , ∥ x − x ∗ ( t )∥ ≤ ρ , ∥ λ − λ 0 ∥ ≤ ρ }
は D × Λ D \times \Lambda D × Λ に含まれるコンパクト集合になる。補題 2.5 と同じ議論で、U U U 上で x x x についてのリプシッツ定数 L > 0 L > 0 L > 0 がとれる。f f f は U U U 上で一様連続なので、η > 0 \eta > 0 η > 0 と δ ∈ ( 0 , ρ ] \delta \in (0, \rho] δ ∈ ( 0 , ρ ] を ( δ + η / L ) e L ( b − a ) < ρ (\delta + \eta/L)e^{L(b - a)} < \rho ( δ + η / L ) e L ( b − a ) < ρ となるように、かつ ∥ λ − λ 0 ∥ < δ \lVert \lambda - \lambda_0 \rVert < \delta ∥ λ − λ 0 ∥ < δ なら U U U 上で ∥ f ( t , x , λ ) − f ( t , x , λ 0 ) ∥ ≤ η \lVert f(t, x, \lambda) - f(t, x, \lambda_0) \rVert \leq \eta ∥ f ( t , x , λ ) − f ( t , x , λ 0 )∥ ≤ η となるようにとれる。
(ii) ∥ ξ − ξ 0 ∥ < δ \lVert \xi - \xi_0 \rVert < \delta ∥ ξ − ξ 0 ∥ < δ , ∥ λ − λ 0 ∥ < δ \lVert \lambda - \lambda_0 \rVert < \delta ∥ λ − λ 0 ∥ < δ とし、y ( t ) = x ( t ; ξ , λ ) y(t) = x(t; \xi, \lambda) y ( t ) = x ( t ; ξ , λ ) の最大存在区間を ( ω − , ω + ) (\omega_-, \omega_+) ( ω − , ω + ) とする。[ t 0 , min ( b , ω + ) ) [t_0, \min(b, \omega_+)) [ t 0 , min ( b , ω + )) のうち、[ t 0 , τ ] [t_0, \tau] [ t 0 , τ ] 上で ∥ y ( t ) − x ∗ ( t ) ∥ ≤ ρ \lVert y(t) - x^{\ast}(t) \rVert \leq \rho ∥ y ( t ) − x ∗ ( t )∥ ≤ ρ となる τ \tau τ の上限を T T T とする。[ t 0 , T ) [t_0, T) [ t 0 , T ) 上では ( t , y ( t ) , λ ) ∈ U (t, y(t), \lambda) \in U ( t , y ( t ) , λ ) ∈ U なので、系 2.11(g = f ( ⋅ , ⋅ , λ ) g = f(\cdot, \cdot, \lambda) g = f ( ⋅ , ⋅ , λ ) 、比較する方程式は f ( ⋅ , ⋅ , λ 0 ) f(\cdot, \cdot, \lambda_0) f ( ⋅ , ⋅ , λ 0 ) )より
∥ y ( t ) − x ∗ ( t ) ∥ ≤ ( δ + η L ) e L ( t − t 0 ) < ρ \lVert y(t) - x^{\ast}(t) \rVert \leq \left(\delta + \frac{\eta}{L}\right)e^{L(t - t_0)} < \rho ∥ y ( t ) − x ∗ ( t )∥ ≤ ( δ + L η ) e L ( t − t 0 ) < ρ
もし T < ω + T < \omega_+ T < ω + かつ T < b T < b T < b なら、連続性からこの不等式は T T T の少し先でも成り立ち、T T T の定義に矛盾する。もし T = ω + ≤ b T = \omega_+ \leq b T = ω + ≤ b なら、[ t 0 , ω + ) [t_0, \omega_+) [ t 0 , ω + ) 上で y y y のグラフはコンパクト集合 { ( t , x ) ∣ ( t , x , λ ) ∈ U } ⊂ D \lbrace (t, x) \mid (t, x, \lambda) \in U \rbrace \subset D {( t , x ) ∣ ( t , x , λ ) ∈ U } ⊂ D に留まり、定理 2.16 に矛盾する。よって T = b < ω + T = b < \omega_+ T = b < ω + で、評価は [ t 0 , b ] [t_0, b] [ t 0 , b ] で成り立つ。[ a , t 0 ] [a, t_0] [ a , t 0 ] も同様であり、ρ < ε \rho < \varepsilon ρ < ε から主張が従う。□ \square □
λ = λ 0 \lambda = \lambda_0 λ = λ 0 の場合、証明から sup [ a , b ] ∥ x ( t ; ξ , λ 0 ) − x ∗ ( t ) ∥ ≤ e L ( b − a ) ∥ ξ − ξ 0 ∥ \sup_{[a, b]}\lVert x(t; \xi, \lambda_0) - x^{\ast}(t) \rVert \leq e^{L(b - a)}\lVert \xi - \xi_0 \rVert sup [ a , b ] ∥ x ( t ; ξ , λ 0 ) − x ∗ ( t )∥ ≤ e L ( b − a ) ∥ ξ − ξ 0 ∥ というリプシッツ連続性 が得られる。同様に、( t , t 0 , ξ , λ ) ↦ x ( t ; t 0 , ξ , λ ) (t, t_0, \xi, \lambda) \mapsto x(t; t_0, \xi, \lambda) ( t , t 0 , ξ , λ ) ↦ x ( t ; t 0 , ξ , λ ) はその定義域(開集合になる)上で連続であることが示される(証明は省略)。これは第4章でフローを扱う基礎である。
さらに滑らかさを仮定すると、解は初期値について微分可能になる。
定理 2.22 (初期値についての微分可能性)f f f が ( x , λ ) (x, \lambda) ( x , λ ) について C 1 C^1 C 1 級で、その偏導関数が ( t , x , λ ) (t, x, \lambda) ( t , x , λ ) について連続なら、x ( t ; ξ , λ ) x(t; \xi, \lambda) x ( t ; ξ , λ ) は ( t , ξ , λ ) (t, \xi, \lambda) ( t , ξ , λ ) について C 1 C^1 C 1 級である。Y ( t ) = ∂ x / ∂ ξ ( t ; ξ 0 , λ 0 ) Y(t) = \partial x/\partial\xi\ (t; \xi_0, \lambda_0) Y ( t ) = ∂ x / ∂ ξ ( t ; ξ 0 , λ 0 ) は変分方程式 (variational equation)
Y ′ = D x f ( t , x ∗ ( t ) , λ 0 ) Y , Y ( t 0 ) = I Y' = D_xf(t, x^{\ast}(t), \lambda_0)\,Y, \qquad Y(t_0) = I Y ′ = D x f ( t , x ∗ ( t ) , λ 0 ) Y , Y ( t 0 ) = I
の解であり、Z = ∂ x / ∂ λ Z = \partial x/\partial\lambda Z = ∂ x / ∂ λ は Z ′ = D x f Z + D λ f Z' = D_xf\ Z + D_\lambda f Z ′ = D x f Z + D λ f , Z ( t 0 ) = O Z(t_0) = O Z ( t 0 ) = O の解である。
証明の方針. 方程式を形式的に ξ \xi ξ で微分すれば変分方程式が得られる。正当化には、差分商 q h ( t ) = ( x ( t ; ξ 0 + h v , λ 0 ) − x ∗ ( t ) ) / h q_h(t) = \bigl(x(t; \xi_0 + hv, \lambda_0) - x^{\ast}(t)\bigr)/h q h ( t ) = ( x ( t ; ξ 0 + h v , λ 0 ) − x ∗ ( t ) ) / h が、平均値の定理から q h ( t ) = v + ∫ t 0 t A h ( s ) q h ( s ) d s q_h(t) = v + \int_{t_0}^t A_h(s)q_h(s)\ ds q h ( t ) = v + ∫ t 0 t A h ( s ) q h ( s ) d s 、A h ( s ) = ∫ 0 1 D x f ( s , x ∗ + θ ( x h − x ∗ ) , λ 0 ) d θ A_h(s) = \int_0^1 D_xf\bigl(s, x^{\ast} + \theta(x_h - x^{\ast}), \lambda_0\bigr)d\theta A h ( s ) = ∫ 0 1 D x f ( s , x ∗ + θ ( x h − x ∗ ) , λ 0 ) d θ をみたすことに注目する。定理 2.20 から A h → D x f ( s , x ∗ ( s ) , λ 0 ) A_h \to D_xf(s, x^{\ast}(s), \lambda_0) A h → D x f ( s , x ∗ ( s ) , λ 0 ) が一様に成り立つので、系 2.11 と同様のグロンウォール評価によって q h q_h q h は変分方程式の解 Y ( t ) v Y(t)v Y ( t ) v に一様収束する。詳細は常微分方程式の標準的な教科書(README の参考文献)を参照のこと。
例 2.23 x ′ = x 2 x' = x^2 x ′ = x 2 , x ( 0 ) = ξ x(0) = \xi x ( 0 ) = ξ の解は x = ξ / ( 1 − ξ t ) x = \xi/(1 - \xi t) x = ξ / ( 1 − ξ t ) で、∂ x / ∂ ξ = ( 1 − ξ t ) − 2 \partial x/\partial\xi = (1 - \xi t)^{-2} ∂ x / ∂ ξ = ( 1 − ξ t ) − 2 である。一方、変分方程式 Y ′ = 2 x ( t ) Y Y' = 2x(t)Y Y ′ = 2 x ( t ) Y , Y ( 0 ) = 1 Y(0) = 1 Y ( 0 ) = 1 を解くと Y = exp ∫ 0 t 2 ξ 1 − ξ s d s = ( 1 − ξ t ) − 2 Y = \exp\int_0^t \frac{2\xi}{1 - \xi s}ds = (1 - \xi t)^{-2} Y = exp ∫ 0 t 1 − ξ s 2 ξ d s = ( 1 − ξ t ) − 2 となり、一致する。
2.8 ペアノの存在定理
リプシッツ条件がなくても、解の存在だけは保証される。
定理 2.24 (ペアノの存在定理, Peano existence theorem)f f f が R = [ t 0 − a , t 0 + a ] × B ‾ ( x 0 , b ) R = [t_0 - a, t_0 + a] \times \overline{B}(x_0, b) R = [ t 0 − a , t 0 + a ] × B ( x 0 , b ) 上で連続なら、M = max R ∥ f ∥ M = \max_R\lVert f \rVert M = max R ∥ f ∥ , h = min ( a , b / M ) h = \min(a, b/M) h = min ( a , b / M ) (M = 0 M = 0 M = 0 なら h = a h = a h = a )として、初期値問題 x ′ = f ( t , x ) x' = f(t, x) x ′ = f ( t , x ) , x ( t 0 ) = x 0 x(t_0) = x_0 x ( t 0 ) = x 0 は [ t 0 − h , t 0 + h ] [t_0 - h, t_0 + h] [ t 0 − h , t 0 + h ] 上に少なくとも一つの解をもつ。
証明. [ t 0 , t 0 + h ] [t_0, t_0 + h] [ t 0 , t 0 + h ] 上の解を作る(左側も同様に作り、t 0 t_0 t 0 でつなげれば、補題 2.1 の積分方程式をみたすので全体で解になる)。N ∈ N N \in \mathbb{N} N ∈ N に対し t k = t 0 + k h / N t_k = t_0 + kh/N t k = t 0 + k h / N (0 ≤ k ≤ N 0 \leq k \leq N 0 ≤ k ≤ N )とおき、折れ線 φ N \varphi_N φ N を
φ N ( t 0 ) = x 0 , φ N ( t ) = φ N ( t k ) + ( t − t k ) f ( t k , φ N ( t k ) ) ( t k ≤ t ≤ t k + 1 ) \varphi_N(t_0) = x_0, \qquad \varphi_N(t) = \varphi_N(t_k) + (t - t_k)f(t_k, \varphi_N(t_k)) \quad (t_k \leq t \leq t_{k+1}) φ N ( t 0 ) = x 0 , φ N ( t ) = φ N ( t k ) + ( t − t k ) f ( t k , φ N ( t k )) ( t k ≤ t ≤ t k + 1 )
で定める(オイラー折れ線 , Euler polygon)。
(i) k k k についての帰納法により、節点 ( t k , φ N ( t k ) ) (t_k, \varphi_N(t_k)) ( t k , φ N ( t k )) は R R R に含まれ、各線分の傾きのノルムは M M M 以下なので、∥ φ N ( t ) − x 0 ∥ ≤ M ( t − t 0 ) ≤ b \lVert \varphi_N(t) - x_0 \rVert \leq M(t - t_0) \leq b ∥ φ N ( t ) − x 0 ∥ ≤ M ( t − t 0 ) ≤ b である。よって φ N \varphi_N φ N のグラフは R R R に含まれる。
(ii) ∥ φ N ( t ) − φ N ( s ) ∥ ≤ M ∣ t − s ∣ \lVert \varphi_N(t) - \varphi_N(s) \rVert \leq M\lvert t - s \rvert ∥ φ N ( t ) − φ N ( s )∥ ≤ M ∣ t − s ∣ なので、{ φ N } \lbrace \varphi_N \rbrace { φ N } は同程度連続かつ一様有界である。アスコリ–アルツェラの定理(集合と位相 定理 5.38 ;各成分に順に適用して部分列を n n n 回とり直す)により、一様収束する部分列 φ N j → φ \varphi_{N_j} \to \varphi φ N j → φ がある。φ \varphi φ は連続で、グラフは R R R に含まれる。
(iii) s ∈ [ t k , t k + 1 ) s \in [t_k, t_{k+1}) s ∈ [ t k , t k + 1 ) に対し τ N ( s ) = t k \tau_N(s) = t_k τ N ( s ) = t k とおくと、φ N \varphi_N φ N は区分的に 1 次なので
φ N ( t ) = x 0 + ∫ t 0 t f ( τ N ( s ) , φ N ( τ N ( s ) ) ) d s \varphi_N(t) = x_0 + \int_{t_0}^t f\bigl(\tau_N(s), \varphi_N(\tau_N(s))\bigr)\,ds φ N ( t ) = x 0 + ∫ t 0 t f ( τ N ( s ) , φ N ( τ N ( s )) ) d s
である。ここで ∣ τ N ( s ) − s ∣ ≤ h / N \lvert \tau_N(s) - s \rvert \leq h/N ∣ τ N ( s ) − s ∣ ≤ h / N 、また ∥ φ N ( τ N ( s ) ) − φ ( s ) ∥ ≤ M h / N + sup ∥ φ N − φ ∥ \lVert \varphi_N(\tau_N(s)) - \varphi(s) \rVert \leq Mh/N + \sup\lVert \varphi_N - \varphi \rVert ∥ φ N ( τ N ( s )) − φ ( s )∥ ≤ M h / N + sup ∥ φ N − φ ∥ である。f f f はコンパクト集合 R R R 上で一様連続なので、N = N j → ∞ N = N_j \to \infty N = N j → ∞ のとき被積分関数は f ( s , φ ( s ) ) f(s, \varphi(s)) f ( s , φ ( s )) に一様収束する。極限をとれば φ ( t ) = x 0 + ∫ t 0 t f ( s , φ ( s ) ) d s \varphi(t) = x_0 + \int_{t_0}^t f(s, \varphi(s))\ ds φ ( t ) = x 0 + ∫ t 0 t f ( s , φ ( s )) d s で、補題 2.1 より φ \varphi φ は解である。□ \square □
オイラー折れ線は数値解法のオイラー法 そのものである。解が一意でなければ列全体が収束するとは限らないが、一意なら、どの部分列も一意な解へ収束する部分列をもつので列全体が収束する。
ペアノの定理は無限次元のバナッハ空間では成り立たない(ディユドネの反例)。証明の本質は、有限次元性に由来するコンパクト性(アスコリ–アルツェラの定理)にある。
まとめ
初期値問題は積分方程式 x ( t ) = x 0 + ∫ t 0 t f ( s , x ( s ) ) d s x(t) = x_0 + \int_{t_0}^t f(s, x(s))\ ds x ( t ) = x 0 + ∫ t 0 t f ( s , x ( s )) d s と同値であり、存在定理はこの形で証明する。
ピカール–リンデレフの定理:f f f が連続で x x x についてリプシッツなら、解は局所的に存在して一意である。証明は逐次近似の一様収束、または重みつき距離での縮小写像の原理による。
グロンウォールの不等式により、2 つの解の差は高々 e L ∣ t − t 0 ∣ e^{L\lvert t - t_0 \rvert} e L ∣ t − t 0 ∣ 倍にしか広がらない。一意性と誤差評価がここから従う。
局所リプシッツなら最大解が一意に存在し、有限時刻で終わる最大解は定義域のどのコンパクト集合からも出ていく。D = J × R n D = J \times \mathbb{R}^n D = J × R n ではこれは解の爆発を意味し、線形増大条件の下では爆発は起きない。
解は初期値とパラメータに連続に依存し、f f f が C 1 C^1 C 1 級なら微分可能である(変分方程式)。
ペアノの定理:f f f が連続なだけでも、オイラー折れ線とアスコリ–アルツェラの定理により解が存在する。
演習問題
問題 2.1 ★ x ′ = 2 t ( 1 + x ) x' = 2t(1 + x) x ′ = 2 t ( 1 + x ) , x ( 0 ) = 0 x(0) = 0 x ( 0 ) = 0 の逐次近似 φ k \varphi_k φ k を求め、その極限が解であることを確かめよ。
解答
φ 1 = ∫ 0 t 2 s d s = t 2 \varphi_1 = \int_0^t 2s\ ds = t^2 φ 1 = ∫ 0 t 2 s d s = t 2 、φ 2 = ∫ 0 t 2 s ( 1 + s 2 ) d s = t 2 + t 4 / 2 \varphi_2 = \int_0^t 2s(1 + s^2)\ ds = t^2 + t^4/2 φ 2 = ∫ 0 t 2 s ( 1 + s 2 ) d s = t 2 + t 4 /2 、φ 3 = t 2 + t 4 / 2 + t 6 / 6 \varphi_3 = t^2 + t^4/2 + t^6/6 φ 3 = t 2 + t 4 /2 + t 6 /6 。帰納法で φ k = ∑ j = 1 k t 2 j / j ! \varphi_k = \sum_{j=1}^k t^{2j}/j! φ k = ∑ j = 1 k t 2 j / j ! を示す:これを仮定すると φ k + 1 = ∫ 0 t 2 s ( 1 + ∑ j = 1 k s 2 j j ! ) d s = t 2 + ∑ j = 1 k t 2 j + 2 ( j + 1 ) ! \varphi_{k+1} = \int_0^t 2s\left(1 + \sum_{j=1}^k \frac{s^{2j}}{j!}\right)ds = t^2 + \sum_{j=1}^k \frac{t^{2j+2}}{(j+1)!} φ k + 1 = ∫ 0 t 2 s ( 1 + ∑ j = 1 k j ! s 2 j ) d s = t 2 + ∑ j = 1 k ( j + 1 )! t 2 j + 2 。極限は e t 2 − 1 e^{t^2} - 1 e t 2 − 1 で、x ′ = 2 t e t 2 = 2 t ( 1 + x ) x' = 2te^{t^2} = 2t(1 + x) x ′ = 2 t e t 2 = 2 t ( 1 + x ) 、x ( 0 ) = 0 x(0) = 0 x ( 0 ) = 0 をみたす。
問題 2.2 ★★ u , β u, \beta u , β を [ t 0 , t 1 ] [t_0, t_1] [ t 0 , t 1 ] 上の連続関数、β ≥ 0 \beta \geq 0 β ≥ 0 、α \alpha α を連続な単調増加関数とする。u ( t ) ≤ α ( t ) + ∫ t 0 t β ( s ) u ( s ) d s u(t) \leq \alpha(t) + \int_{t_0}^t \beta(s)u(s)\ ds u ( t ) ≤ α ( t ) + ∫ t 0 t β ( s ) u ( s ) d s ならば u ( t ) ≤ α ( t ) exp ∫ t 0 t β ( s ) d s u(t) \leq \alpha(t)\exp\int_{t_0}^t \beta(s)\ ds u ( t ) ≤ α ( t ) exp ∫ t 0 t β ( s ) d s であることを示せ。
解答
T ∈ [ t 0 , t 1 ] T \in [t_0, t_1] T ∈ [ t 0 , t 1 ] を固定する。t ∈ [ t 0 , T ] t \in [t_0, T] t ∈ [ t 0 , T ] では α ( t ) ≤ α ( T ) \alpha(t) \leq \alpha(T) α ( t ) ≤ α ( T ) なので u ( t ) ≤ α ( T ) + ∫ t 0 t β u d s u(t) \leq \alpha(T) + \int_{t_0}^t \beta u\ ds u ( t ) ≤ α ( T ) + ∫ t 0 t β u d s 。区間 [ t 0 , T ] [t_0, T] [ t 0 , T ] で定理 2.10 を定数 α ( T ) \alpha(T) α ( T ) に適用すると u ( t ) ≤ α ( T ) exp ∫ t 0 t β d s u(t) \leq \alpha(T)\exp\int_{t_0}^t \beta\ ds u ( t ) ≤ α ( T ) exp ∫ t 0 t β d s 。t = T t = T t = T とおけば、T T T は任意なので主張を得る。
問題 2.3 ★★ x ′ = 3 x 2 / 3 x' = 3x^{2/3} x ′ = 3 x 2/3 , x ( 0 ) = 0 x(0) = 0 x ( 0 ) = 0 のすべての解を求めよ(x 2 / 3 = ( x 1 / 3 ) 2 ≥ 0 x^{2/3} = (x^{1/3})^2 \geq 0 x 2/3 = ( x 1/3 ) 2 ≥ 0 とする)。
解答
0 ≤ c ≤ ∞ 0 \leq c \leq \infty 0 ≤ c ≤ ∞ , 0 ≤ d ≤ ∞ 0 \leq d \leq \infty 0 ≤ d ≤ ∞ に対し、t > c t > c t > c で x = ( t − c ) 3 x = (t - c)^3 x = ( t − c ) 3 、− d ≤ t ≤ c -d \leq t \leq c − d ≤ t ≤ c で x = 0 x = 0 x = 0 、t < − d t < -d t < − d で x = ( t + d ) 3 x = (t + d)^3 x = ( t + d ) 3 とおいたものが解である(x = ( t − c ) 3 x = (t - c)^3 x = ( t − c ) 3 なら x ′ = 3 ( t − c ) 2 = 3 x 2 / 3 x' = 3(t - c)^2 = 3x^{2/3} x ′ = 3 ( t − c ) 2 = 3 x 2/3 で、つなぎ目の微分係数は 0)。逆に任意の解は x ′ ≥ 0 x' \geq 0 x ′ ≥ 0 より単調増加で、x > 0 x > 0 x > 0 の範囲では右辺が局所リプシッツなので変数分離で x 1 / 3 = t − c x^{1/3} = t - c x 1/3 = t − c と一意に決まり、x < 0 x < 0 x < 0 の範囲も同様である。x = 0 x = 0 x = 0 となる t t t の集合は単調性から区間 [ − d , c ] [-d, c] [ − d , c ] であり、解は上の形に限る。
問題 2.4 ★★ f : R → R f\colon \mathbb{R} \to \mathbb{R} f : R → R が連続で単調減少なら、x ′ = f ( x ) x' = f(x) x ′ = f ( x ) , x ( t 0 ) = x 0 x(t_0) = x_0 x ( t 0 ) = x 0 の解は t ≥ t 0 t \geq t_0 t ≥ t 0 の範囲で一意であることを示せ。また f ( x ) = − x 1 / 3 f(x) = -x^{1/3} f ( x ) = − x 1/3 のとき、t ≤ t 0 t \leq t_0 t ≤ t 0 の範囲では一意性が成り立たないことを示せ。
解答
x , y x, y x , y を同じ初期値の解とすると、d d t ( x − y ) 2 = 2 ( x − y ) ( f ( x ) − f ( y ) ) ≤ 0 \frac{d}{dt}(x - y)^2 = 2(x - y)\bigl(f(x) - f(y)\bigr) \leq 0 d t d ( x − y ) 2 = 2 ( x − y ) ( f ( x ) − f ( y ) ) ≤ 0 (f f f は単調減少)なので、( x − y ) 2 (x - y)^2 ( x − y ) 2 は t ≥ t 0 t \geq t_0 t ≥ t 0 で増えず、t 0 t_0 t 0 で 0 だから恒等的に 0 である。
f ( x ) = − x 1 / 3 f(x) = -x^{1/3} f ( x ) = − x 1/3 で x ( t 0 ) = 0 x(t_0) = 0 x ( t 0 ) = 0 とする。任意の s ≤ t 0 s \leq t_0 s ≤ t 0 に対し、t ≤ s t \leq s t ≤ s で x s ( t ) = ( 2 3 ( s − t ) ) 3 / 2 x_s(t) = \left(\frac{2}{3}(s - t)\right)^{3/2} x s ( t ) = ( 3 2 ( s − t ) ) 3/2 、s ≤ t ≤ t 0 s \leq t \leq t_0 s ≤ t ≤ t 0 で x s ( t ) = 0 x_s(t) = 0 x s ( t ) = 0 とおく。t < s t < s t < s では x s ′ = − ( 2 3 ( s − t ) ) 1 / 2 = − x s 1 / 3 x_s' = -\left(\frac{2}{3}(s - t)\right)^{1/2} = -x_s^{1/3} x s ′ = − ( 3 2 ( s − t ) ) 1/2 = − x s 1/3 であり、t = s t = s t = s で左右の微分係数はともに 0 なので、x s x_s x s は ( − ∞ , t 0 ] (-\infty, t_0] ( − ∞ , t 0 ] 上の解で x s ( t 0 ) = 0 x_s(t_0) = 0 x s ( t 0 ) = 0 をみたす。s s s ごとに異なる解なので、過去向きの一意性は成り立たない(「いつ 0 に到達したか」が決まらない)。
問題 2.5 ★★ (1) x ′ ′ + x 3 = 0 x'' + x^3 = 0 x ′′ + x 3 = 0 のすべての解は R \mathbb{R} R 全体で定義されることを示せ。(2) x ′ ′ = x 3 x'' = x^3 x ′′ = x 3 には有限時刻で爆発する解があることを示せ。
解答
v = x ′ v = x' v = x ′ として 1 階の系 x ′ = v x' = v x ′ = v , v ′ = ∓ x 3 v' = \mp x^3 v ′ = ∓ x 3 を考える(右辺は C 1 C^1 C 1 級)。
(1) E = 1 2 v 2 + 1 4 x 4 E = \frac{1}{2}v^2 + \frac{1}{4}x^4 E = 2 1 v 2 + 4 1 x 4 は E ′ = v v ′ + x 3 x ′ = − v x 3 + x 3 v = 0 E' = vv' + x^3x' = -vx^3 + x^3v = 0 E ′ = v v ′ + x 3 x ′ = − v x 3 + x 3 v = 0 より保存される。よって ( x , v ) (x, v) ( x , v ) は有界な集合 { E ≤ E ( 0 ) } \lbrace E \leq E(0) \rbrace { E ≤ E ( 0 )} に留まり、系 2.17 より爆発しない。
(2) x ( 0 ) = 1 x(0) = 1 x ( 0 ) = 1 , x ′ ( 0 ) = 1 / 2 x'(0) = 1/\sqrt{2} x ′ ( 0 ) = 1/ 2 とすると 1 2 x ′ 2 − 1 4 x 4 = 0 \frac{1}{2}x'^2 - \frac{1}{4}x^4 = 0 2 1 x ′2 − 4 1 x 4 = 0 が保存されるので、x ′ = x 2 / 2 x' = x^2/\sqrt{2} x ′ = x 2 / 2 の解 x = 1 / ( 1 − t / 2 ) x = 1/(1 - t/\sqrt{2}) x = 1/ ( 1 − t / 2 ) を試すと、x ′ ′ = 2 x x ′ = x 3 x'' = \sqrt{2}\ xx' = x^3 x ′′ = 2 x x ′ = x 3 で初期条件もみたす。一意性からこれが解であり、t → 2 t \to \sqrt{2} t → 2 で爆発する。
問題 2.6 ★★★ f : R × R n → R n f\colon \mathbb{R} \times \mathbb{R}^n \to \mathbb{R}^n f : R × R n → R n が連続かつ有界なら、任意の初期値問題は R \mathbb{R} R 全体で定義された解をもつことを示せ。
解答
∥ f ∥ ≤ M \lVert f \rVert \leq M ∥ f ∥ ≤ M とする。刻み幅 1 / N 1/N 1/ N のオイラー折れ線 φ N \varphi_N φ N を t 0 t_0 t 0 から左右に R \mathbb{R} R 全体で定義すると(f f f が全空間で定義されているので、グラフが外に出る心配はない)、∥ φ N ( t ) − φ N ( s ) ∥ ≤ M ∣ t − s ∣ \lVert \varphi_N(t) - \varphi_N(s) \rVert \leq M\lvert t - s \rvert ∥ φ N ( t ) − φ N ( s )∥ ≤ M ∣ t − s ∣ 、∥ φ N ( t ) − x 0 ∥ ≤ M ∣ t − t 0 ∣ \lVert \varphi_N(t) - x_0 \rVert \leq M\lvert t - t_0 \rvert ∥ φ N ( t ) − x 0 ∥ ≤ M ∣ t − t 0 ∣ が成り立つ。各 m ∈ N m \in \mathbb{N} m ∈ N について、[ t 0 − m , t 0 + m ] [t_0 - m, t_0 + m] [ t 0 − m , t 0 + m ] 上でアスコリ–アルツェラの定理を適用し、m = 1 , 2 , … m = 1, 2, \dots m = 1 , 2 , … と順に部分列をとり直していく。m m m 番目の部分列の m m m 番目の項を並べた対角線列は、すべての [ t 0 − m , t 0 + m ] [t_0 - m, t_0 + m] [ t 0 − m , t 0 + m ] 上で一様収束する。その極限 φ \varphi φ は、各コンパクト区間上で定理 2.24 の証明 (iii) と同じ議論により積分方程式をみたすので、R \mathbb{R} R 全体での解である。