§E10.6Grönwall の不等式と連続依存性

最終更新

同じ増加率に従う二つの量でも、出発時のわずかな差は時間とともに拡大する。たとえば、y′=yy'=yに従う二つの解について、初期値の差から有限時間内の解の差をどのように評価するかが問題となる。

この例では差を直接計算することができるが、非線形方程式では一般に解の公式を得ることができない。それでも、差の増え方を右辺の Lipschitz 評価によって抑えれば、解そのものを求めずに誤差の上界を与えることができる。Grönwall の不等式はこの推論を可能にし、初期値だけでなく方程式の右辺を摂動したときにも、解の連続依存性を定量化する。

本記事では、Grönwall の不等式を導入し、有限の時間区間における初期値と右辺の摂動に対する解の連続依存性を明らかにする。

1 積分形の Grönwall 不等式

定理 1.1 (Grönwall の不等式).a<ba<b、A≥0A\geq0とし、u:[a,b]→[0,∞)u:[a,b]\to[0,\infty)を連続関数、ℓ,r∈L1([a,b])\ell,r\in L^1([a,b])を非負関数とする。すべてのt∈[a,b]t\in[a,b]で

u(t)≤A+∫at(ℓ(s)u(s)+r(s)) dsu(t)\leq A+\int_a^t\bigl(\ell(s)u(s)+r(s)\bigr)\,ds

が成り立つならば、

u(t)≤Aexp⁡(∫atℓ(s) ds)+∫atexp⁡(∫stℓ(τ) dτ)r(s) ds(a≤t≤b)u(t)\leq A\exp\left(\int_a^t\ell(s)\,ds\right) +\int_a^t\exp\left(\int_s^t\ell(\tau)\,d\tau\right)r(s)\,ds \qquad(a\leq t\leq b)

である。

また、すべてのt∈[a,b]t\in[a,b]で

u(t)≤A+∫tb(ℓ(s)u(s)+r(s)) dsu(t)\leq A+\int_t^b\bigl(\ell(s)u(s)+r(s)\bigr)\,ds

が成り立つならば、

u(t)≤Aexp⁡(∫tbℓ(s) ds)+∫tbexp⁡(∫tsℓ(τ) dτ)r(s) ds(a≤t≤b)u(t)\leq A\exp\left(\int_t^b\ell(s)\,ds\right) +\int_t^b\exp\left(\int_t^s\ell(\tau)\,d\tau\right)r(s)\,ds \qquad(a\leq t\leq b)

である。

証明.

v(t)=A+∫at(ℓ(s)u(s)+r(s)) ds,H(t)=∫atℓ(s) dsv(t)=A+\int_a^t\bigl(\ell(s)u(s)+r(s)\bigr)\,ds, \qquad H(t)=\int_a^t\ell(s)\,ds

とおく。uuは有界であり、ℓu+r\ell u+rは可積分であるから、§E10.2 定理 2.4によりv,Hv,Hは絶対連続である。右向きの不等式を仮定すると、0≤u≤v0\leq u\leq vであり、ほとんど至る所で

v′=ℓu+r≤ℓv+r,H′=ℓv'=\ell u+r\leq\ell v+r,\qquad H'=\ell

となる。§E10.2 定理 3.1 (1)と§E10.2 定理 3.1 (2)によりe−Hve^{-H}vは絶対連続であり、ほとんど至る所で

(e−Hv)′=e−H(v′−ℓv)≤e−Hr(e^{-H}v)'=e^{-H}(v'-\ell v)\leq e^{-H}r

が成り立つ。§E10.2 定理 2.4の全点積分表示から

e−H(t)v(t)≤A+∫ate−H(s)r(s) dse^{-H(t)}v(t)\leq A+\int_a^t e^{-H(s)}r(s)\,ds

を得る。eH(t)e^{H(t)}を掛け、u(t)≤v(t)u(t)\leq v(t)とH(t)−H(s)=∫stℓ(τ) dτH(t)-H(s)=\int_s^t\ell(\tau)\,d\tauを用いると、右向きの評価を得る。

左向きの不等式を仮定し、

w(t)=A+∫tb(ℓ(s)u(s)+r(s)) dsw(t)=A+\int_t^b\bigl(\ell(s)u(s)+r(s)\bigr)\,ds

とおく。wwは絶対連続で、0≤u≤w0\leq u\leq wである。ほとんど至る所で

w′=−ℓu−r≥−ℓw−r,(eHw)′=eH(w′+ℓw)≥−eHrw'=-\ell u-r\geq-\ell w-r,\qquad (e^Hw)'=e^H(w'+\ell w)\geq-e^Hr

となる。ttからbbまで積分し、w(b)=Aw(b)=Aを用いると

eH(t)w(t)≤eH(b)A+∫tbeH(s)r(s) dse^{H(t)}w(t)\leq e^{H(b)}A+\int_t^b e^{H(s)}r(s)\,ds

である。e−H(t)e^{-H(t)}を掛け、u(t)≤w(t)u(t)\leq w(t)を用いて左向きの評価を得る。▨

系 1.2 (Grönwall の不等式).a<ba<bとし、u,r:[a,b]→[0,∞)u,r:[a,b]\to[0,\infty)を連続関数、A,L≥0A,L\ge0とする。任意のt∈[a,b]t\in[a,b]について

u(t)≤A+∫at(Lu(s)+r(s)) dsu(t)\le A+\int_a^t\bigl(Lu(s)+r(s)\bigr)\,ds

が成り立つならば、

u(t)≤AeL(t−a)+∫ateL(t−s)r(s) ds(a≤t≤b)u(t)\le A e^{L(t-a)} +\int_a^t e^{L(t-s)}r(s)\,ds \qquad(a\le t\le b)

である。

証明.定理 1.1にℓ(t)=L\ell(t)=Lを代入する。連続関数rrは有限閉区間上で可積分であり、∫stℓ(τ) dτ=L(t−s)\int_s^t\ell(\tau)\,d\tau=L(t-s)であるから、結論を得る。▨

系 1.3.a<ba<b、A,L≥0A,L\geq0とし、u:[a,b]→[0,∞)u:[a,b]\to[0,\infty)が連続で、

u(t)≤A+L∫atu(s) dsu(t)\le A+L\int_a^t u(s)\,ds

を満たすならば、

u(t)≤AeL(t−a)u(t)\le A e^{L(t-a)}

である。とくにA=0A=0ならば、uuは[a,b][a,b]上で恒等的に零である。

証明.系 1.2でr=0r=0と置けばよい。▨

指数因子は単なる証明上の上界ではない。L≥0L\geq0、z(a)≥0z(a)\geq0に対して、方程式z′=Lzz'=Lzの解z(t)=z(a)eL(t−a)z(t)=z(a)e^{L(t-a)}では等号が成り立つため、一般にこの増大率をより小さい形へ置き換えることはできない。

2 二つの絶対連続解の比較

定理 2.1.d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\|\cdot\|を固定する。a<ba<b、J=[a,b]J=[a,b]、t0∈Jt_0\in J、S⊆J×RdS\subseteq J\times\R^dとし、f,g:S→Rdf,g:S\to\R^dを写像とする。y,z:J→Rd\boldsymbol y,\boldsymbol z:J\to\R^dは絶対連続で、そのグラフはSSに含まれ、ほとんど至る所で

y′=f(t,y),z′=g(t,z)\boldsymbol y'=f(t,\boldsymbol y),\qquad \boldsymbol z'=g(t,\boldsymbol z)

を満たすとする。非負関数ℓ,r∈L1(J)\ell,r\in L^1(J)と Lebesgue 零集合N⊆JN\subseteq Jが存在し、すべてのt∈J∖Nt\in J\setminus Nと、(t,u),(t,v)∈S(t,\boldsymbol u),(t,\boldsymbol v)\in Sを満たすすべての状態対に対して

∥g(t,u)−g(t,v)∥≤ℓ(t)∥u−v∥,∥g(t,y(t))−f(t,y(t))∥≤r(t)\|g(t,\boldsymbol u)-g(t,\boldsymbol v)\| \leq\ell(t)\|\boldsymbol u-\boldsymbol v\|, \qquad \|g(t,\boldsymbol y(t))-f(t,\boldsymbol y(t))\|\leq r(t)

が成り立つと仮定する。δ=∥z(t0)−y(t0)∥\delta=\|\boldsymbol z(t_0)-\boldsymbol y(t_0)\|とおくと、t0≤t≤bt_0\leq t\leq bで

∥z(t)−y(t)∥≤δexp⁡(∫t0tℓ(s) ds)+∫t0texp⁡(∫stℓ(τ) dτ)r(s) ds,\|\boldsymbol z(t)-\boldsymbol y(t)\| \leq\delta\exp\left(\int_{t_0}^t\ell(s)\,ds\right) +\int_{t_0}^t\exp\left(\int_s^t\ell(\tau)\,d\tau\right)r(s)\,ds,

a≤t≤t0a\leq t\leq t_0で

∥z(t)−y(t)∥≤δexp⁡(∫tt0ℓ(s) ds)+∫tt0exp⁡(∫tsℓ(τ) dτ)r(s) ds\|\boldsymbol z(t)-\boldsymbol y(t)\| \leq\delta\exp\left(\int_t^{t_0}\ell(s)\,ds\right) +\int_t^{t_0}\exp\left(\int_t^s\ell(\tau)\,d\tau\right)r(s)\,ds

となる。特に

max⁡t∈J∥z(t)−y(t)∥≤exp⁡(∫abℓ(s) ds)(δ+∫abr(s) ds)\max_{t\in J}\|\boldsymbol z(t)-\boldsymbol y(t)\| \leq\exp\left(\int_a^b\ell(s)\,ds\right) \left(\delta+\int_a^b r(s)\,ds\right)

である。

証明.§E10.2 定理 2.4により、二つの解の差は全点で

z(t)−y(t)=z(t0)−y(t0)+∫t0t(g(s,z(s))−f(s,y(s))) ds\boldsymbol z(t)-\boldsymbol y(t) =\boldsymbol z(t_0)-\boldsymbol y(t_0) +\int_{t_0}^t\bigl(g(s,\boldsymbol z(s))-f(s,\boldsymbol y(s))\bigr)\,ds

と表される。両合成は可積分な導関数とほとんど至る所で一致するから可積分である。共通零集合NNの外で

∥g(s,z(s))−f(s,y(s))∥≤∥g(s,z(s))−g(s,y(s))∥+∥g(s,y(s))−f(s,y(s))∥≤ℓ(s)∥z(s)−y(s)∥+r(s)\begin{aligned} \|g(s,\boldsymbol z(s))-f(s,\boldsymbol y(s))\| &\leq\|g(s,\boldsymbol z(s))-g(s,\boldsymbol y(s))\| +\|g(s,\boldsymbol y(s))-f(s,\boldsymbol y(s))\|\\ &\leq\ell(s)\|\boldsymbol z(s)-\boldsymbol y(s)\|+r(s) \end{aligned}

である。u(t)=∥z(t)−y(t)∥u(t)=\|\boldsymbol z(t)-\boldsymbol y(t)\|は非負連続であり、§E10.4 補題 3.1により

u(t)≤δ+∫t0t(ℓu+r)(s) ds(t≥t0),u(t)≤δ+∫tt0(ℓu+r)(s) ds(t≤t0)u(t)\leq\delta+\int_{t_0}^t(\ell u+r)(s)\,ds\quad(t\geq t_0), \qquad u(t)\leq\delta+\int_t^{t_0}(\ell u+r)(s)\,ds\quad(t\leq t_0)

を満たす。定理 1.1を[t0,b][t_0,b]と[a,t0][a,t_0]のうち長さが正の区間へ適用して、二つの評価を得る。t=t0t=t_0の評価はu(t0)=δu(t_0)=\deltaによる。ℓ,r≥0\ell,r\geq0なので、各部分区間の積分を[a,b][a,b]の積分で抑えれば最後の評価が従う。▨

3 閉管内の共通存在区間

定理 3.1.d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\|\cdot\|を固定する。Ω⊆R×Rd\Omega\subseteq\R\times\R^dを開集合とし、f,g:Ω→Rdf,g:\Omega\to\R^dはそれぞれ Carathéodory 条件と時間について可積分な Lipschitz 係数の条件を局所的に満たすとする。a<ba<b、J=[a,b]J=[a,b]とし、y:J→Rd\boldsymbol y:J\to\R^dをy′=f(t,y)\boldsymbol y'=f(t,\boldsymbol y)の絶対連続解とする。ρ>0\rho>0に対して

K={(t,u)∈J×Rd:∥u−y(t)∥≤ρ}⊆ΩK=\{(t,\boldsymbol u)\in J\times\R^d: \|\boldsymbol u-\boldsymbol y(t)\|\leq\rho\}\subseteq\Omega

と仮定する。

非負関数ℓ,r∈L1(J)\ell,r\in L^1(J)と Lebesgue 零集合N⊆JN\subseteq Jが存在し、すべてのt∈J∖Nt\in J\setminus Nと、(t,u),(t,v)∈K(t,\boldsymbol u),(t,\boldsymbol v)\in Kを満たすすべての状態対に対して

∥g(t,u)−g(t,v)∥≤ℓ(t)∥u−v∥,∥g(t,y(t))−f(t,y(t))∥≤r(t)\|g(t,\boldsymbol u)-g(t,\boldsymbol v)\| \leq\ell(t)\|\boldsymbol u-\boldsymbol v\|, \qquad \|g(t,\boldsymbol y(t))-f(t,\boldsymbol y(t))\|\leq r(t)

が成り立つとする。η∈Rd\boldsymbol\eta\in\R^dに対して

δ=∥η−y(a)∥,Q=exp⁡(∫abℓ(s) ds)(δ+∫abr(s) ds)<ρ\delta=\|\boldsymbol\eta-\boldsymbol y(a)\|, \qquad Q=\exp\left(\int_a^b\ell(s)\,ds\right) \left(\delta+\int_a^b r(s)\,ds\right)<\rho

ならば、初期値問題

z′=g(t,z),z(a)=η\boldsymbol z'=g(t,\boldsymbol z),\qquad \boldsymbol z(a)=\boldsymbol\eta

はJJ全体に絶対連続解をもち、

max⁡t∈J∥z(t)−y(t)∥≤Q<ρ\max_{t\in J}\|\boldsymbol z(t)-\boldsymbol y(t)\|\leq Q<\rho

となる。この解は同じ初期値をもつJJ上の絶対連続解の中で一意である。

証明.δ≤Q<ρ\delta\leq Q<\rhoより(a,η)∈K⊆Ω(a,\boldsymbol\eta)\in K\subseteq\Omegaである。§E10.5 定理 2.1により、z(a)=η\boldsymbol z(a)=\boldsymbol\etaを満たす最大絶対連続解z:(α,β)→Rd\boldsymbol z:(\alpha,\beta)\to\R^dを取る。α<a<β\alpha<a<\betaである。また、KKは連続写像

J×{v:∥v∥≤ρ}⟶R×Rd,(t,v)⟼(t,y(t)+v)J\times\{\boldsymbol v:\|\boldsymbol v\|\leq\rho\} \longrightarrow\R\times\R^d, \qquad (t,\boldsymbol v)\longmapsto(t,\boldsymbol y(t)+\boldsymbol v)

によるコンパクト集合の像なのでコンパクトである。

a<c≤ba<c\leq b、c<βc<\betaを満たすccを取り、u(t)=∥z(t)−y(t)∥u(t)=\|\boldsymbol z(t)-\boldsymbol y(t)\|とおく。[a,c][a,c]にu(t)≥ρu(t)\geq\rhoとなる点が存在すると仮定する。uuの連続性により、集合{t∈[a,c]:u(t)≥ρ}\{t\in[a,c]:u(t)\geq\rho\}は最小元τ\tauをもち、u(a)=δ<ρu(a)=\delta<\rhoからa<τa<\tau、u(τ)=ρu(\tau)=\rhoである。a≤t≤τa\leq t\leq\tauでは二つの解のグラフはKKに含まれる。定理 2.1を[a,τ][a,\tau]とS=K∩([a,τ]×Rd)S=K\cap([a,\tau]\times\R^d)に適用すると

u(τ)≤exp⁡(∫aτℓ(s) ds)(δ+∫aτr(s) ds)≤Q<ρu(\tau)\leq \exp\left(\int_a^\tau\ell(s)\,ds\right) \left(\delta+\int_a^\tau r(s)\,ds\right) \leq Q<\rho

となり、u(τ)=ρu(\tau)=\rhoと両立しない。従って[a,c][a,c]全体でu<ρu<\rhoである。両グラフがKKに含まれるので、同じ比較定理を[a,c][a,c]全体へ適用してu(t)≤Qu(t)\leq Qを得る。

β≤b\beta\leq bと仮定する。任意のt∈[a,β)t\in[a,\beta)に対してt<c<βt<c<\betaとなるccを取ると、直前の評価から(t,z(t))∈K(t,\boldsymbol z(t))\in Kである。β\betaは有限で、K⊆ΩK\subseteq\Omegaはコンパクトであるから、§E10.5 定理 4.1によりz\boldsymbol zはβ\betaを越えて延長することができる。最大性に反するのでβ>b\beta>bである。c=bc=bとして得られる評価が結論の不等式である。

同じ初期値をもつ別のJJ上の絶対連続解z~\widetilde{\boldsymbol z}を取る。z~\widetilde{\boldsymbol z}に対しても、c=bc=bとして最初の離脱時刻を取る議論を適用すると、そのグラフはKKに含まれる。定理 2.1をS=KS=K上の二つの右辺g,gg,gと二解z,z~\boldsymbol z,\widetilde{\boldsymbol z}に適用する。初期値の差と右辺の差はともに零なので、z~=z∣J\widetilde{\boldsymbol z}=\boldsymbol z|_Jである。▨

系 3.2.d∈N≥1d\in\NNとし、Rd\R^dにノルムを固定する。Ω⊆R×Rd\Omega\subseteq\R\times\R^dを開集合とし、f:Ω→Rdf:\Omega\to\R^dと各n∈N≥1n\in\NNに対するfn:Ω→Rdf_n:\Omega\to\R^dは、それぞれ Carathéodory 条件と時間について可積分な Lipschitz 係数の条件を局所的に満たすとする。a<ba<b、J=[a,b]J=[a,b]、ρ>0\rho>0とし、y:J→Rd\boldsymbol y:J\to\R^dをffに対する絶対連続解とする。その閉管

K={(t,u):t∈J, ∥u−y(t)∥≤ρ}K=\{(t,\boldsymbol u):t\in J,\ \|\boldsymbol u-\boldsymbol y(t)\|\leq\rho\}

はΩ\Omegaに含まれるとする。

非負関数ℓ∈L1(J)\ell\in L^1(J)、各n∈N≥1n\in\NNに対する非負関数rn∈L1(J)r_n\in L^1(J)、Lebesgue 零集合N⊆JN\subseteq Jが存在し、すべてのn∈N≥1n\in\NN、t∈J∖Nt\in J\setminus Nと、(t,u),(t,v)∈K(t,\boldsymbol u),(t,\boldsymbol v)\in Kを満たすすべての状態対に対して

∥fn(t,u)−fn(t,v)∥≤ℓ(t)∥u−v∥,∥fn(t,y(t))−f(t,y(t))∥≤rn(t)\|f_n(t,\boldsymbol u)-f_n(t,\boldsymbol v)\| \leq\ell(t)\|\boldsymbol u-\boldsymbol v\|, \qquad \|f_n(t,\boldsymbol y(t))-f(t,\boldsymbol y(t))\|\leq r_n(t)

が成り立つと仮定する。初期値列ηn∈Rd\boldsymbol\eta_n\in\R^dが

δn=∥ηn−y(a)∥,δn+∫abrn(s) ds⟶0\delta_n=\|\boldsymbol\eta_n-\boldsymbol y(a)\|, \qquad \delta_n+\int_a^b r_n(s)\,ds\longrightarrow0

を満たすならば、十分大きいすべてのnnに対して

yn′=fn(t,yn),yn(a)=ηn\boldsymbol y_n'=f_n(t,\boldsymbol y_n),\qquad \boldsymbol y_n(a)=\boldsymbol\eta_n

の絶対連続解yn\boldsymbol y_nがJJ全体に一意に存在し、

max⁡t∈J∥yn(t)−y(t)∥≤exp⁡(∫abℓ(s) ds)(δn+∫abrn(s) ds)⟶0\max_{t\in J}\|\boldsymbol y_n(t)-\boldsymbol y(t)\| \leq\exp\left(\int_a^b\ell(s)\,ds\right) \left(\delta_n+\int_a^b r_n(s)\,ds\right)\longrightarrow0

となる。

証明.

Qn=exp⁡(∫abℓ(s) ds)(δn+∫abrn(s) ds)Q_n=\exp\left(\int_a^b\ell(s)\,ds\right) \left(\delta_n+\int_a^b r_n(s)\,ds\right)

とおく。指数因子はnnに依存しない有限な定数なので、Qn→0Q_n\to0である。従って、あるn0∈N≥1n_0\in\NNが存在し、n≥n0n\geq n_0ならばQn<ρQ_n<\rhoとなる。各n≥n0n\geq n_0に定理 3.1をg=fng=f_n、r=rnr=r_n、η=ηn\boldsymbol\eta=\boldsymbol\eta_nとして適用すると、共通区間JJ上の一意な解とQnQ_nによる評価を得る。▨

4 連続な右辺と一様誤差の場合

系 4.1.d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\lVert\cdot\rVertを固定する。a<ba<b、J=[a,b]J=[a,b]、D⊆RdD\subseteq\R^dとし、

f,g:J×D→Rdf,g:J\times D\to\R^d

を連続とする。ある共通の定数L≥0L\ge0に対して

∥f(t,u)−f(t,v)∥≤L∥u−v∥,∥g(t,u)−g(t,v)∥≤L∥u−v∥\begin{aligned} \lVert f(t,\boldsymbol u)-f(t,\boldsymbol v)\rVert &\le L\lVert\boldsymbol u-\boldsymbol v\rVert,\\ \lVert g(t,\boldsymbol u)-g(t,\boldsymbol v)\rVert &\le L\lVert\boldsymbol u-\boldsymbol v\rVert \end{aligned}

が、t∈Jt\in Jとu,v∈D\boldsymbol u,\boldsymbol v\in Dについて成り立つと仮定する。また、

ε=sup⁡(t,u)∈J×D∥f(t,u)−g(t,u)∥<∞\varepsilon =\sup_{(t,\boldsymbol u)\in J\times D} \lVert f(t,\boldsymbol u)-g(t,\boldsymbol u)\rVert <\infty

とする。

y,z:J→D\boldsymbol y,\boldsymbol z:J\to Dがそれぞれ

y′=f(t,y),z′=g(t,z)\boldsymbol y'=f(t,\boldsymbol y),\qquad \boldsymbol z'=g(t,\boldsymbol z)

を満たすならば、t∈Jt\in Jに対して

∥y(t)−z(t)∥≤eL(t−a)∥y(a)−z(a)∥+ε∫ateL(t−s) ds\lVert\boldsymbol y(t)-\boldsymbol z(t)\rVert \le e^{L(t-a)} \lVert\boldsymbol y(a)-\boldsymbol z(a)\rVert +\varepsilon\int_a^t e^{L(t-s)}\,ds

である。すなわち、L>0L>0ならば

∥y(t)−z(t)∥≤eL(t−a)∥y(a)−z(a)∥+εL(eL(t−a)−1),\lVert\boldsymbol y(t)-\boldsymbol z(t)\rVert \le e^{L(t-a)} \lVert\boldsymbol y(a)-\boldsymbol z(a)\rVert +\frac{\varepsilon}{L}\bigl(e^{L(t-a)}-1\bigr),

L=0L=0ならば

∥y(t)−z(t)∥≤∥y(a)−z(a)∥+ε(t−a)\lVert\boldsymbol y(t)-\boldsymbol z(t)\rVert \le \lVert\boldsymbol y(a)-\boldsymbol z(a)\rVert +\varepsilon(t-a)

である。

証明. 二つの解の導関数は連続であるから、両解はJJ上で絶対連続である。定理 2.1にS=J×DS=J\times D、t0=at_0=a、ℓ=L\ell=L、r=εr=\varepsilonを代入すると、最初の評価を得る。L>0L>0ならば

∫ateL(t−s) ds=eL(t−a)−1L,\int_a^t e^{L(t-s)}\,ds=\frac{e^{L(t-a)}-1}{L},

L=0L=0ならば積分はt−at-aである。従って残りの二式も成り立つ。▨

系 4.2.d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\lVert\cdot\rVertを固定する。I⊆RI\subseteq\Rを開区間、D⊆RdD\subseteq\R^dを開集合とし、a<ba<b、J=[a,b]⊆IJ=[a,b]\subseteq Iとする。f,g:I×D→Rdf,g:I\times D\to\R^dは連続であり、状態変数について局所 Lipschitz 連続であるとする。y:J→D\boldsymbol y:J\to Dをy′=f(t,y)\boldsymbol y'=f(t,\boldsymbol y)の解とし、あるr>0r>0に対して

K={(t,u)∈J×Rd:∥u−y(t)∥≤r}K=\{(t,\boldsymbol u)\in J\times\R^d: \lVert\boldsymbol u-\boldsymbol y(t)\rVert\le r\}

がI×DI\times Dに含まれるコンパクト集合であると仮定する。定数L,ε≥0L,\varepsilon\ge0が存在し、任意のt∈Jt\in Jとu,v∈D\boldsymbol u,\boldsymbol v\in Dに対して

∥f(t,u)−f(t,v)∥≤L∥u−v∥,∥g(t,u)−g(t,v)∥≤L∥u−v∥,∥f(t,u)−g(t,u)∥≤ε\begin{aligned} \lVert f(t,\boldsymbol u)-f(t,\boldsymbol v)\rVert &\le L\lVert\boldsymbol u-\boldsymbol v\rVert,\\ \lVert g(t,\boldsymbol u)-g(t,\boldsymbol v)\rVert &\le L\lVert\boldsymbol u-\boldsymbol v\rVert,\\ \lVert f(t,\boldsymbol u)-g(t,\boldsymbol u)\rVert &\le\varepsilon \end{aligned}

が成り立つとする。η∈Rd\boldsymbol\eta\in\R^dに対して

δ=∥η−y(a)∥,Q=eL(b−a)(δ+(b−a)ε)<r\delta=\lVert\boldsymbol\eta-\boldsymbol y(a)\rVert, \qquad Q=e^{L(b-a)}\bigl(\delta+(b-a)\varepsilon\bigr)<r

ならば、初期値問題

z′=g(t,z),z(a)=η\boldsymbol z'=g(t,\boldsymbol z),\qquad \boldsymbol z(a)=\boldsymbol\eta

はJJ上に解z:J→D\boldsymbol z:J\to Dをもち、

max⁡t∈J∥z(t)−y(t)∥≤Q<r\max_{t\in J}\lVert\boldsymbol z(t)-\boldsymbol y(t)\rVert\le Q<r

が成り立つ。

証明.§E10.5 補題 5.1によりf,gf,gは局所的な Carathéodory 条件と可積分 Lipschitz 係数の条件を満たす。定理 3.1にΩ=I×D\Omega=I\times D、ρ=r\rho=r、ℓ(t)=L\ell(t)=L、誤差関数として定数ε\varepsilonを代入する。このとき

exp⁡(∫abℓ(s) ds)(δ+∫abε ds)=eL(b−a)(δ+(b−a)ε)=Q<r\exp\left(\int_a^b\ell(s)\,ds\right) \left(\delta+\int_a^b\varepsilon\,ds\right) =e^{L(b-a)}\bigl(\delta+(b-a)\varepsilon\bigr)=Q<r

であるから、JJ上の絶対連続解と評価を得る。ggは連続であるため、得られた解は§E10.5 補題 5.1により古典解である。▨

系 4.3.d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\lVert\cdot\rVertを固定する。a<ba<b、J=[a,b]J=[a,b]、D⊆RdD\subseteq\R^d、L≥0L\geq0とする。写像f:J×D→Rdf:J\times D\to\R^dが状態変数について Lipschitz 定数LLをもち、y,z:J→D\boldsymbol y,\boldsymbol z:J\to Dが同じ方程式

u′=f(t,u)\boldsymbol u'=f(t,\boldsymbol u)

の古典解ならば、

∥y(t)−z(t)∥≤eL(t−a)∥y(a)−z(a)∥(t∈J)\lVert\boldsymbol y(t)-\boldsymbol z(t)\rVert \le e^{L(t-a)} \lVert\boldsymbol y(a)-\boldsymbol z(a)\rVert \qquad(t\in J)

である。

証明.定理 2.1にS=J×DS=J\times D、g=fg=f、t0=at_0=a、ℓ=L\ell=L、r=0r=0を代入する。古典解は絶対連続であるから、同定理の評価が成り立つ。▨

系 4.4.d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\lVert\cdot\rVertを固定する。a<ba<b、J=[a,b]J=[a,b]、D⊆RdD\subseteq\R^dとする。連続写像f:J×D→Rdf:J\times D\to\R^dと、各n∈N≥0n\in\Nに対する連続写像

fn:J×D→Rdf_n:J\times D\to\R^d

が、状態変数について共通の Lipschitz 定数L≥0L\geq0をもつと仮定する。また、各n∈N≥0n\in\Nに対してyn:J→D\boldsymbol y_n:J\to Dをfnf_nが定める方程式の解とし、y:J→D\boldsymbol y:J\to Dをffが定める方程式の解とする。すなわち、

yn′=fn(t,yn),y′=f(t,y)\boldsymbol y_n'=f_n(t,\boldsymbol y_n),\qquad \boldsymbol y'=f(t,\boldsymbol y)

を満たす解であるとする。初期値誤差と右辺の一様誤差を

δn=∥yn(a)−y(a)∥,εn=sup⁡(t,u)∈J×D∥fn(t,u)−f(t,u)∥\delta_n=\lVert\boldsymbol y_n(a)-\boldsymbol y(a)\rVert, \qquad \varepsilon_n=\sup_{(t,\boldsymbol u)\in J\times D} \lVert f_n(t,\boldsymbol u)-f(t,\boldsymbol u)\rVert

と置く。各nnについてεn<∞\varepsilon_n<\inftyであり、

δn⟶0,εn⟶0\delta_n\longrightarrow0, \qquad \varepsilon_n\longrightarrow0

ならば、

max⁡t∈J∥yn(t)−y(t)∥≤eL(b−a)(δn+(b−a)εn)⟶0\max_{t\in J}\lVert\boldsymbol y_n(t)-\boldsymbol y(t)\rVert \leq e^{L(b-a)}\bigl(\delta_n+(b-a)\varepsilon_n\bigr) \longrightarrow0

である。したがって、yn\boldsymbol y_nはJJ上でy\boldsymbol yへ一様収束する。特に、fn=ff_n=fとすれば初期値だけを変える場合が得られ、yn(a)=y(a)\boldsymbol y_n(a)=\boldsymbol y(a)とすれば右辺だけを変える場合が得られる。

証明. 各nnに対して定理 2.1にg=fng=f_n、z=yn\boldsymbol z=\boldsymbol y_n、S=J×DS=J\times D、t0=at_0=a、ℓ=L\ell=L、r=εnr=\varepsilon_nを代入すると、

max⁡t∈J∥yn(t)−y(t)∥≤eL(b−a)(δn+(b−a)εn)\max_{t\in J}\lVert\boldsymbol y_n(t)-\boldsymbol y(t)\rVert \leq e^{L(b-a)}\bigl(\delta_n+(b-a)\varepsilon_n\bigr)

を得る。右辺は零へ収束する。LLによる除算を含まないため、L=0L=0の場合もこの評価に含まれる。▨

ここで仮定した共通区間上の存在は、初期値や右辺の収束だけからは従わない。既知解の閉管が開定義域に含まれ、右辺が局所存在条件を満たす場合には、系 3.2が十分大きい添字に対する共通区間上の存在を与える。

5 計算例

例 5.1 (線形方程式における評価の鋭さ).T>0T>0、λ≥0\lambda\ge0、η,y0,z0∈R\eta,y_0,z_0\in\Rとし、

y′=λy,z′=λz+ηy'=\lambda y,\qquad z'=\lambda z+\eta

を[0,T][0,T]上で考える。二つの初期値をy(0)=y0y(0)=y_0、z(0)=z0z(0)=z_0とする。右辺の状態変数に関する共通 Lipschitz 定数はL=λL=\lambdaであり、同じ状態における右辺の差は∣η∣|\eta|である。系 4.1は、λ>0\lambda>0のとき

∣y(t)−z(t)∣≤eλt∣y0−z0∣+∣η∣λ(eλt−1)|y(t)-z(t)| \le e^{\lambda t}|y_0-z_0| +\frac{|\eta|}{\lambda}(e^{\lambda t}-1)

を与える。

λ>0\lambda>0のとき、実際の解は

y(t)=eλty0,z(t)=eλtz0+ηλ(eλt−1)y(t)=e^{\lambda t}y_0,\qquad z(t)=e^{\lambda t}z_0+\frac{\eta}{\lambda}(e^{\lambda t}-1)

である。y0−z0y_0-z_0と−η-\etaの符号が同じ場合には三角不等式で損失がなく、上の評価で等号が成り立つ。したがって、指数因子と右辺摂動の積分項はいずれも一般に必要である。

λ=0\lambda=0のときはy(t)=y0y(t)=y_0、z(t)=z0+ηtz(t)=z_0+\eta tであり、

∣y(t)−z(t)∣≤∣y0−z0∣+∣η∣t|y(t)-z(t)|\leq|y_0-z_0|+|\eta|t

となる。

例 5.2.n∈N≥1n\in\NNに対して、R\R上の関数を

rn(t)=n1(0,1/n2)(t)r_n(t)=n\mathbf1_{(0,1/n^2)}(t)

とおく。f(t,u)=0f(t,u)=0、fn(t,u)=rn(t)f_n(t,u)=r_n(t)はR×R\R\times\R上で局所的な Carathéodory 条件を満たし、状態についての Lipschitz 係数はℓ=0\ell=0である。J=[0,1]J=[0,1]、基準解y=0y=0、ρ=1\rho=1、初期値ηn=0\eta_n=0とすると、閉管は[0,1]×[−1,1][0,1]\times[-1,1]であり、

sup⁡t∈Jrn(t)=n,∫01rn(t) dt=1n\sup_{t\in J}r_n(t)=n,\qquad \int_0^1r_n(t)\,dt=\frac1n

となる。従って系 3.2は、外力が一様には零へ収束しなくても、解がJJ上で零へ一様収束することを与える。

実際、

yn(t)=∫0trn(s) ds=nmin⁡{t,1/n2}(0≤t≤1)y_n(t)=\int_0^t r_n(s)\,ds =n\min\{t,1/n^2\}\qquad(0\leq t\leq1)

は絶対連続で、ほとんど至る所でyn′=rny_n'=r_n、yn(0)=0y_n(0)=0を満たす。従って

max⁡0≤t≤1∣yn(t)∣=1n\max_{0\leq t\leq1}|y_n(t)|=\frac1n

であり、積分誤差による上界はこの例で等号となる。

6 演習

問題 6.1 (固定した区間と増大する時刻). 各n∈N≥1n\in\NNに対して、初期値問題

yn′=yn,yn(0)=1ny_n'=y_n,\qquad y_n(0)=\frac1n

の解をyn(t)=et/ny_n(t)=e^t/nとし、同じ方程式の零解をy(t)=0y(t)=0とする。任意の固定したT>0T>0に対して、yny_nが[0,T][0,T]上でyyへ一様収束することを示せ。一方、tn=log⁡nt_n=\log nにおける差を計算し、[0,∞)[0,\infty)上での一様収束は成り立たないことを示せ。

解答.

yn′(t)=et/n=yn(t)y_n'(t)=e^t/n=y_n(t)、yn(0)=1/ny_n(0)=1/nであるから、yny_nは指定された初期値問題の解である。零関数もy′=yy'=yの解である。ete^tは増加するので、任意の固定したT>0T>0について

max⁡0≤t≤T∣yn(t)−y(t)∣=eTn⟶0\max_{0\le t\le T}|y_n(t)-y(t)|=\frac{e^T}{n}\longrightarrow0

である。これは系 4.4においてJ=[0,T]J=[0,T]、D=RD=\R、fn(t,u)=f(t,u)=uf_n(t,u)=f(t,u)=u、L=1L=1、δn=1/n\delta_n=1/n、εn=0\varepsilon_n=0とした評価にも一致する。

n≥1n\ge1よりtn=log⁡n≥0t_n=\log n\ge0であり、

∣yn(tn)−y(tn)∣=elog⁡nn=1|y_n(t_n)-y(t_n)|=\frac{e^{\log n}}{n}=1

である。したがって、どのnnに対しても[0,∞)[0,\infty)上のすべての時刻で差を1/21/2未満にすることはできず、この区間上での一様収束は成り立たない。tn→∞t_n\to\inftyであるから、固定された[0,T][0,T]上での一様収束とは両立する。▨

問題 6.2 (Grönwall の証明と摂動評価).a<ba<bとし、連続関数u:[a,b]→[0,∞)u:[a,b]\to[0,\infty)が

u(t)≤A+∫at(Lu(s)+r(s)) dsu(t)\le A+\int_a^t\bigl(Lu(s)+r(s)\bigr)\,ds

を満たすとする。ここでA,L≥0A,L\ge0であり、r:[a,b]→[0,∞)r:[a,b]\to[0,\infty)は連続である。補助関数

v(t)=A+∫at(Lu(s)+r(s)) dsv(t)=A+\int_a^t\bigl(Lu(s)+r(s)\bigr)\,ds

と積分因子を用いて、積分形の Grönwall の不等式を再証明せよ。

次に、d∈N≥1d\in\NNとし、Rd\R^dに一つのノルムを固定する。連続な右辺に対するy′=f(t,y)\boldsymbol y'=f(t,\boldsymbol y)とz′=g(t,z)\boldsymbol z'=g(t,\boldsymbol z)の解が共通区間J=[a,b]J=[a,b]上で共通集合DDに値を取り、ffとggがDD上で共通の Lipschitz 定数LLをもつとする。さらに

sup⁡(t,u)∈J×D∥f(t,u)−g(t,u)∥≤ε\sup_{(t,\boldsymbol u)\in J\times D} \lVert f(t,\boldsymbol u)-g(t,\boldsymbol u)\rVert\le\varepsilon

とする。二つの解の差を評価し、L=0L=0の場合も明記せよ。

解答.

仮定から0≤u(t)≤v(t)0\le u(t)\le v(t)であり、

v′(t)=Lu(t)+r(t)≤Lv(t)+r(t)v'(t)=Lu(t)+r(t)\le Lv(t)+r(t)

である。したがって、

ddt(e−L(t−a)v(t))≤e−L(t−a)r(t)\frac{d}{dt}\left(e^{-L(t-a)}v(t)\right) \le e^{-L(t-a)}r(t)

となる。aaからttまで積分してv(a)=Av(a)=Aを用いると、

u(t)≤v(t)≤AeL(t−a)+∫ateL(t−s)r(s) dsu(t)\le v(t) \le A e^{L(t-a)} +\int_a^t e^{L(t-s)}r(s)\,ds

を得る。

二つの積分方程式を差し引き、積分内の差を

f(s,y(s))−g(s,z(s))=f(s,y(s))−f(s,z(s))+f(s,z(s))−g(s,z(s))f(s,\boldsymbol y(s))-g(s,\boldsymbol z(s)) =f(s,\boldsymbol y(s))-f(s,\boldsymbol z(s)) +f(s,\boldsymbol z(s))-g(s,\boldsymbol z(s))

と分ける。Grönwall の不等式をu(t)=∥y(t)−z(t)∥u(t)=\lVert\boldsymbol y(t)-\boldsymbol z(t)\rVert、A=∥y(a)−z(a)∥A=\lVert\boldsymbol y(a)-\boldsymbol z(a)\rVert、r(t)=εr(t)=\varepsilonに適用すると、

∥y(t)−z(t)∥≤eL(t−a)∥y(a)−z(a)∥+ε∫ateL(t−s) ds\lVert\boldsymbol y(t)-\boldsymbol z(t)\rVert \le e^{L(t-a)}\lVert\boldsymbol y(a)-\boldsymbol z(a)\rVert +\varepsilon\int_a^t e^{L(t-s)}\,ds

を得る。L>0L>0ならば右辺の第二項はε(eL(t−a)−1)/L\varepsilon(e^{L(t-a)}-1)/Lであり、L=0L=0ならば評価は

∥y(t)−z(t)∥≤∥y(a)−z(a)∥+ε(t−a)\lVert\boldsymbol y(t)-\boldsymbol z(t)\rVert \le\lVert\boldsymbol y(a)-\boldsymbol z(a)\rVert+\varepsilon(t-a)

となる。共通区間、共通集合、共通ノルム、および共通の Lipschitz 定数を用いることが、この比較には必要である。▨

参考文献

  1. Wolfgang Walter, Ordinary Differential Equations, translated by Russell Thompson, Graduate Texts in Mathematics 182, Springer, 1998.Grönwall の積分不等式と解の連続依存性の評価を参考にした。
  2. Gerald Teschl, Ordinary Differential Equations and Dynamical Systems, American Mathematical Society, Providence, Rhode Island, 2012.Grönwall 不等式、初期値に関する安定性、右辺の摂動評価を参考にした。

前提記事