logo

一階常微分方程式の初期値問題に対する解の存在性と一意性 📂微分方程式

一階常微分方程式の初期値問題に対する解の存在性と一意性

定理1

EERn\mathbb{R}^{n}開集合であり、fC1(E)f \in C^{1} (E)ϕ0E\phi_{0} \in Eに関する以下のような初期値問題が与えられたとする。

{ϕ˙=f(ϕ)ϕ(0)=ϕ0 \begin{cases} \dot{ \phi } = \mathbf{f} ( \phi ) \\ \phi (0) = \phi_{0} \end{cases}

すると、ある区間[h,h]R[-h,h] \subset \mathbb{R}で与えられた初期値問題の解ϕ(t)\phi (t)は一意に存在する。


証明

戦略:存在性を先に示した後に、一意性を示します。記法として、ユークリッド空間Rn\mathbb{R}^{n}ボールB(x0;d):={xRn:x0x<d}B[x0;d]:={xRn:x0xd} \begin{align*} B \left( \mathbf{x}_{0} ; d \right) :=& \left\{ \mathbf{x} \in \mathbb{R}^{n} : | \mathbf{x}_{0} - \mathbf{x} | < d \right\} \\ B \left[ \mathbf{x}_{0} ; d \right] :=& \left\{ \mathbf{x} \in \mathbb{R}^{n} : | \mathbf{x}_{0} - \mathbf{x} | \le d \right\} \end{align*} と表示します。当然ながら、fC1(E)f \in C^{1} (E)の代わりにリプシッツ条件を使うことも問題ありません。


パート1. ffは局所的にリプシッツである

局所的リプシッツ条件fC1(E)f \in C^{1}(E)のとき、ffEEで局所的リプシッツである。

fy\dfrac{\partial f}{\partial y}が連続と仮定されているので、定理によりffは局所的リプシッツである。よって、局所的リプシッツの定義により、全てのx,yB(ϕ0;ε)E\mathbf{x} , \mathbf{y} \in B \left( \phi_{0} ; \varepsilon \right) \subset Eに対して、以下の式を満たすε,K>0\varepsilon, K > 0が存在する。

f(t,y1)f(t,y2)Ky1y2 \left| f(t,y_{1}) - f(t,y_{2}) \right| \le K \left| y_{1} - y_{2} \right|

連続性とコンパクトの関係XXコンパクト距離空間、YY距離空間f:XYf:X\to Y連続であるとする。すると、f(X)f(X)はコンパクトである。

ffは連続であるので、コンパクトセットB:=B[ϕ0;ε2]B : = B \left[ \phi_{0} ; {{ \varepsilon } \over {2}} \right]ではバウンデッドであり、M:=supxBf(x)\displaystyle M : = \sup_{ \mathbf{x} \in B } | f ( \mathbf{x} ) |得ることができる。


パート2. ピカールの方法

EERn\mathbb{R}^{n} でオープンであり、fC1(E)f \in C^{1} (E)についての初期値問題

{ϕ˙=f(ϕ)ϕ(0)=ϕ0 \begin{cases} \dot{ \phi } = f ( \phi ) \\ \phi (0) = \phi_{0} \end{cases} があるとする。関数の {uk(t)}k=0\left\{ \mathbf{u}_{k} (t) \right\} _{ k =0}^{ \infty }

{u0(t)=ϕ0uk+1(t)=ϕ0+0tf(uk(s))ds \begin{cases} \mathbf{u}_{0} (t) = \phi_{0} \\ \displaystyle \mathbf{u}_{k+1} (t) = \phi_{0} + \int_{0}^{t} f \left( \mathbf{u}_{k} (s) \right) ds \end{cases}

のように定義すると、連続関数u(t):=limkuk(t)\displaystyle \mathbf{u} (t) := \lim_{k \to \infty} \mathbf{u}_{k} (t)は与えられた初期値問題の解である。

連続関数uk(t)\mathbf{u}_{k} (t)

{u0(t)=ϕ0uk+1(t)=ϕ0+0tf(uk(s))ds \begin{cases} \mathbf{u}_{0} (t) = \phi_{0} \\ \displaystyle \mathbf{u}_{k+1} (t) = \phi_{0} + \int_{0}^{t} f \left( \mathbf{u}_{k} (s) \right) ds \end{cases}

と定義されたと仮定する。

uk\mathbf{u}_{k}ff[h,h][-h,h]で連続なので、(fuk)( f \circ \mathbf{u}_{k} ) もまた連続である。

微積分学の基本定理:関数ffが閉区間[a,b][a,b]で連続であれば、関数F(x)=axf(t)dt\displaystyle F(x) = \int_{a}^{x} f(t) dt[a,b][a,b]で連続であり、(a,b)(a,b)微分可能である。 dF(x)dx=f(x){{dF(x)} \over {dx}} = f(x)

それにより、微積分学の基本定理によって、uk+1(t)=ϕ0+0tf(uk(s))ds\displaystyle \mathbf{u}_{k+1} (t) = \phi_{0} + \int_{0}^{t} f \left( \mathbf{u}_{k} (s) \right) dsもまた[h,h][-h, h]で連続である。少し整理して不等式を立てると、全てのt[h,h]t \in [-h, h]に対して、

uk+1(t)ϕ00tf(uk(s))ds0hf(uk(s))ds0hMdsMh \begin{align*} | \mathbf{u}_{k+1} (t) - \phi_{0} | \le & \int_{0}^{t} \left| f \left( \mathbf{u}_{k} (s) \right) \right| ds \\ \le & \int_{0}^{h} \left| f \left( \mathbf{u}_{k} (s) \right) \right| ds \\ \le & \int_{0}^{h} M ds \\ \le & Mh \end{align*}

つまり、h(0,ε2M]\displaystyle h \in \left( 0 , {{ \varepsilon } \over { 2M }} \right]を選ぶことによって、uk(t)\mathbf{u}_{k} (t)は全てのt[h,h]t \in [-h,h]k=1,2,3k = 1,2,3 \cdotsにおいて連続関数として定義することができる。


パート3. コーシー列 {uk}k=0\left\{ \mathbf{u}_{k} \right\}_{k=0}^{\infty}

t[h,h]t \in [-h , h ]についてuj+1uj| \mathbb{ u }_{ j + 1 } - \mathbb{ u }_{ j } |の上限を計算する。

  • ケース1. j=1j = 1
    u2(t)u1(t)0tf(u1(s))f(u0(s))dsK0tu1(s)u0(s)dsKhsupt[h,h]u1(t)ϕ0Khε2 \begin{align*} | \mathbb{ u }_{ 2 } ( t ) - \mathbb{ u }_{ 1 } ( t ) | \le & \int_{0}^{t} \left| f \left( \mathbf{u}_{1} (s) \right) - f \left( \mathbf{u}_{0} (s) \right) \right| ds \\ \le & K \int_{0}^{t} | \mathbb{ u }_{ 1 } ( s ) - \mathbb{ u }_{ 0 } ( s ) | ds \\ \le & K h \sup_{ t \in [ - h , h ] } | \mathbb{ u }_{ 1 } ( t ) - \phi_{ 0 } | \\ \le & {{ K h \varepsilon } \over { 2 }} \end{align*}
  • ケース2. j>1j > 1
    uj+1(t)uj(t)0tf(uj(s))f(uj1(s))dsK0tuj(s)uj1(s)dsKhsupt[h,h]uj(t)uj1(t) \begin{align*} | \mathbb{ u }_{ j+1 } ( t ) - \mathbb{ u }_{ j } ( t ) | \le & \int_{0}^{t} \left| f \left( \mathbf{u}_{ j } (s) \right) - f \left( \mathbf{u}_{ j - 1 } (s) \right) \right| ds \\ \le & K \int_{0}^{t} | \mathbb{ u }_{ j } ( s ) - \mathbb{ u }_{ j - 1 } ( s ) | ds \\ \le & K h \sup_{ t \in [ - h , h ] } | \mathbb{ u }_{ j } ( t ) - \mathbf{u}_{j-1} ( t ) \end{align*}

再帰的に解くと、**ケース1.**により

uj+1(t)uj(t)(Kh)jε2 | \mathbb{ u }_{ j+1 } ( t ) - \mathbb{ u }_{ j } ( t ) | \le {{ (Kh)^{j} \varepsilon } \over {2}}

m>k>Nm > k > Nh(0,1K)\displaystyle h \in \left( 0 , {{1} \over { K }} \right)とし、c:=Khc := K hとすると、

um(t)uk(t)j=km1uj+1(t)uj(t)j=Nuj+1(t)uj(t)j=N(Kh)jε2=cN1cε2 \begin{align*} | \mathbb{ u }_{ m } ( t ) - \mathbb{ u }_{ k } ( t ) | \le & \sum_{j = k}^{m-1} | \mathbf{u}_{j+1} (t) - \mathbf{u}_{j} (t) | \\ \le & \sum_{j = N}^{ \infty } | \mathbf{u}_{j+1} (t) - \mathbf{u}_{j} (t) | \\ \le & \sum_{j = N}^{ \infty } {{ ( K h )^{ j } \varepsilon } \over {2}} = {{ c^{N} } \over { 1 - c}} {{ \varepsilon } \over {2}} \end{align*}

c<1|c| < 1なので、NN \to \inftyのとき、cN1cε2\displaystyle {{ c^{N} } \over { 1 - c}} {{ \varepsilon } \over {2}}00に収束する。つまり、全てのε>0\varepsilon > 0に対して、

m,kN    umuk=supt[h,h]um(t)uk(t)<ε m , k \ge N \implies \| \mathbf{u}_{m} - \mathbf{u}_{k} \| = \sup_{ t \in [-h , h ] } | \mathbf{u}_{m} ( t ) - \mathbf{u}_{k} ( t ) | < \varepsilon

を満たすNNが存在する。これは{uk}k=0\left\{ \mathbf{u}_{k} \right\}_{k=0}^{\infty}C[h,h]C [ - h , h ]のコーシー列であることを意味する。(もちろん、選択されるhh0<h<min(bM,1K)\displaystyle 0 < h < \min \left( {{b} \over {M}} , {{1} \over {K}} \right)を満たさなければならない。)


パート4. バナッハ空間

C[h,h]C [ -h , h ]はバナッハ空間なので、u(t):=limkuk(t)\displaystyle \mathbf{u} (t) := \lim_{k \to \infty} \mathbf{u}_{k} (t)は連続関数である。u\mathbf{u}が連続であるので、(fu)( f \circ \mathbf{u} ) もまた連続であり、微積分学の基本定理によって

u˙(t)=(ϕ0+0tf(u(s))ds)=f(u(t)) \dot{\mathbf{u} } (t) = \left( \phi_{0} + \int_{0}^{t} f \left( u(s) \right) ds \right)' = f \left( \mathbf{u} (t) \right)

また、t=0t = 0ならばu(0)=ϕ0+00f(u(s))ds=ϕ0+0=ϕ0\displaystyle \mathbf{u} (0) = \phi_{0} + \int_{0}^{0} f \left( u(s) \right) ds = \phi_{0} + 0 = \phi_{0}従って、u\mathbf{u}は全てのt[h,h]t \in [ - h , h ]に対して、与えられた初期値問題の解として存在する。


パート5. 一意性

u\mathbf{u}v\mathbf{v}が与えられた初期値問題の解であるとする。u(t)v(t)| \mathbf{u} (t) - \mathbf{v} (t) |がいくつかのttに対して最大値を持つようにすると、それをt0[h,h]t_{0} \in [-h , h]とする、

uv=supt[h,h]u(t)v(t)=0t0[f(u(s))f(v(s))]ds0t0f(u(s))f(v(s))dsK0t0u(s)v(s)dsKhsups[h,h]u(s)v(s)Khuv \begin{align*} \| \mathbf{u} - \mathbf{v} \| =& \sup_{t \in [-h,h ] } | \mathbf{u} (t) - \mathbf{v} (t) | \\ =& \left| \int_{0}^{t_{0} } \left[ f \left( \mathbf{u} (s) \right) - f \left( \mathbf{v} (s) \right) \right] ds \right| \\ \le & \int_{0}^{t_{0} } \left| f \left( \mathbf{u} (s) \right) - f \left( \mathbf{v} (s) \right) \right| ds \\ \le & K \int_{0}^{t_{0} } \left|\mathbf{u} (s) - \mathbf{v} (s) \right| ds \\ \le & K h \sup_{s \in [-h, h ] } \left| \mathbf{u} (s) - \mathbf{v} (s) \right| \\ \le & K h \| \mathbf{u} - \mathbf{v} \| \end{align*}

要するに、uvKhuv\| \mathbf{u} - \mathbf{v} \| \le K h \| \mathbf{u} - \mathbf{v} \|であり、Kh<1Kh < 1があるため、uv=0\| \mathbf{u} - \mathbf{v} \| = 0でなければならない。従って、[h,h][-h, h ]で、u=v\| \mathbf{u} \| = \| \mathbf{v} \|でなければならない。


  1. ウィリアム・E・ボーイス, Boyce’s Elementary Differential Equations and Boundary Value Problems(第11版、2017年), p83-90 ↩︎