logo

1계 상미분방정식의 초기값 문제에 대한 솔루션의 존재성과 유일성 📂상미분방정식

1계 상미분방정식의 초기값 문제에 대한 솔루션의 존재성과 유일성

정리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) 대신 립시츠 조건이 있어도 무방하다.


Part 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} ) |잡을 수 있다.


Part 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)supt[h,h]uk(t)ϕ0ε2\displaystyle \sup_{ t \in [ -h , h ] } | \mathbf{u}_{k} (t) - \phi_{0} | \le {{\varepsilon } \over {2 }} 를 만족시면서

{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 에서 연속함수로 정의될 수 있다.


Part 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 } | 의 상한을 계산할 것이다.

  • Case 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*}
  • Case 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*}

재귀적으로 풀어내면 Case 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) 을 만족해야한다.)


Part 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 ] 에 대해 주어진 초기값 문제의 솔루션으로써 존재한다.


Part 5. 유일성

u\mathbf{u}v\mathbf{v} 가 주어진 초기값 문제의 솔루션이라고 가정하자. u(t)v(t)| \mathbf{u} (t) - \mathbf{v} (t) | 가 가장 큰 값을 가지도록 하는 ttt0[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. William E. Boyce , Boyce's Elementary Differential Equations and Boundary Value Problems (11th Edition, 2017), p83-90 ↩︎