§E13.7連続依存性

最終更新

y′=yy'=yの二つの解について、時刻ttにおける解の差は、時刻00における差のete^t倍である。初期値の差から有限時間内の解の差をどのように評価するかが問題となる。

y′=yy'=yの例では解の差を直接計算することができるが、非線形方程式では一般に解の公式を得ることができない。それでも、解の差の増え方を右辺の Lipschitz 評価によって抑えれば、解そのものを求めることなく、二つの解の差の上界を得ることができる。この推論には Grönwall の不等式を用いる。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)

である。

証明. すべてのt∈[a,b]t\in[a,b]でu(t)≤A+∫at(ℓu+r)(s) dsu(t)\leq A+\int_a^t(\ell u+r)(s)\,dsが成り立つとし、

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は可積分であるから、§E13.2 定理 2.5により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

となる。§E13.2 定理 3.1 (1)と§E13.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

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

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を用いると、前半の評価を得る。

すべてのt∈[a,b]t\in[a,b]でu(t)≤A+∫tb(ℓu+r)(s) dsu(t)\leq A+\int_t^b(\ell u+r)(s)\,dsが成り立つとする。t∈[a,b]t\in[a,b]に対して

u~(t)=u(a+b−t),ℓ~(t)=ℓ(a+b−t),r~(t)=r(a+b−t)\widetilde u(t)=u(a+b-t),\qquad \widetilde\ell(t)=\ell(a+b-t),\qquad \widetilde r(t)=r(a+b-t)

とおく。反転s↦a+b−ss\mapsto a+b-sは Lebesgue 測度を保つから、ℓ~,r~∈L1([a,b])\widetilde\ell,\widetilde r\in L^1([a,b])は非負であり、u~\widetilde uは非負連続である。t∈[a,b]t\in[a,b]とt′=a+b−tt'=a+b-tについて、積分変数をs=a+b−σs=a+b-\sigmaと置き換えると

u~(t)=u(t′)≤A+∫t′b(ℓu+r)(s) ds=A+∫at(ℓ~ u~+r~)(σ) dσ\widetilde u(t)=u(t')\leq A+\int_{t'}^b(\ell u+r)(s)\,ds =A+\int_a^t(\widetilde\ell\,\widetilde u+\widetilde r)(\sigma)\,d\sigma

である。前半の評価をu~,ℓ~,r~\widetilde u,\widetilde\ell,\widetilde rに適用し、同じ置き換えで

∫atℓ~(σ) dσ=∫t′bℓ(s) ds,∫atexp⁡(∫σtℓ~(τ) dτ)r~(σ) dσ=∫t′bexp⁡(∫t′sℓ(τ) dτ)r(s) ds\int_a^t\widetilde\ell(\sigma)\,d\sigma=\int_{t'}^b\ell(s)\,ds, \qquad \int_a^t\exp\left(\int_\sigma^t\widetilde\ell(\tau)\,d\tau\right)\widetilde r(\sigma)\,d\sigma =\int_{t'}^b\exp\left(\int_{t'}^s\ell(\tau)\,d\tau\right)r(s)\,ds

を用いると、後半の評価がt′t'で成り立つ。ttが[a,b][a,b]を動くときt′t'も[a,b][a,b]を動くから、後半の評価はすべての点で成り立つ。▨

系 1.2.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)を連続関数とする。このとき、すべてのt∈[a,b]t\in[a,b]で

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

が成り立つならば、

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

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

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

注意 1.4. 各t∈[a,b]t\in[a,b]において、系 1.3の結論の係数eL(t−a)e^{L(t-a)}をそれより小さい数に置き換えた評価は成り立たない。実際、A>0A>0、L≥0L\geq0に対してu(t)=AeL(t−a)u(t)=Ae^{L(t-a)}は、すべてのt∈[a,b]t\in[a,b]でu(t)=A+L∫atu(s) dsu(t)=A+L\int_a^tu(s)\,dsを満たし、結論の評価も等号で成り立つ。

例 1.5.[a,b]=[0,1][a,b]=[0,1]とし、A=1A=1、ℓ≡−1\ell\equiv-1、r=0r=0、u≡1/2u\equiv1/2とおく。0≤t≤10\leq t\leq1で

A+∫0t(ℓ(s)u(s)+r(s)) ds=1−t2≥12=u(t)A+\int_0^t\bigl(\ell(s)u(s)+r(s)\bigr)\,ds=1-\frac t2\geq\frac12=u(t)

であるから、uuは定理 1.1の前半の仮定の不等式を満たす。一方、結論の右辺はt=1t=1でAe−1=e−1<1/2=u(1)Ae^{-1}=e^{-1}<1/2=u(1)である。従って、係数ℓ\ellの非負性を外すと定理 1.1の結論は成り立たない。

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は絶対連続であり、ほとんど至る所で

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

を満たすとする。y,z\boldsymbol y,\boldsymbol zのグラフはSSに含まれるとする。非負関数ℓ,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)

である。

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

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

と表される。s↦g(s,z(s))s\mapsto g(s,\boldsymbol z(s))とs↦f(s,y(s))s\mapsto f(s,\boldsymbol y(s))は、それぞれz′\boldsymbol z'とy′\boldsymbol y'にほとんど至る所で一致するから可積分である。共通零集合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)\|は非負連続であり、§E13.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:Ω→Rdf:\Omega\to\R^dを写像、g:Ω→Rdg:\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\}

とおく。K⊆ΩK\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)

が成り立つとする。t∗∈Jt_*\in J、η∈Rd\boldsymbol\eta\in\R^dが

δ=∥η−y(t∗)∥,Q=exp⁡(∫abℓ(s) ds)(δ+∫abr(s) ds)<ρ\delta=\|\boldsymbol\eta-\boldsymbol y(t_*)\|, \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(t∗)=η\boldsymbol z'=g(t,\boldsymbol z),\qquad \boldsymbol z(t_*)=\boldsymbol\eta

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

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

となる。この解はz(t∗)=η\boldsymbol z(t_*)=\boldsymbol\etaを満たすJJ上の絶対連続解の中で一意である。

証明.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)

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

主張 3.1.1.J′⊆JJ'\subseteq Jをt∗t_*を含む長さが正の閉区間とし、w:J′→Rd\boldsymbol w:J'\to\R^dをw′=g(t,w)\boldsymbol w'=g(t,\boldsymbol w)、w(t∗)=η\boldsymbol w(t_*)=\boldsymbol\etaの絶対連続解とする。このとき、すべてのt∈J′t\in J'で(t,w(t))∈K(t,\boldsymbol w(t))\in Kかつ∥w(t)−y(t)∥≤Q\|\boldsymbol w(t)-\boldsymbol y(t)\|\leq Qである。

証明.u(t)=∥w(t)−y(t)∥u(t)=\|\boldsymbol w(t)-\boldsymbol y(t)\|とおき、E={t∈J′:u(t)≥ρ}E=\{t\in J':u(t)\geq\rho\}が空でないと仮定する。uuは連続であるからEEはコンパクトであり、t↦∣t−t∗∣t\mapsto|t-t_*|はEE上で最小値をとる。最小値をとる点をτ\tauとする。u(t∗)=δ≤Q<ρu(t_*)=\delta\leq Q<\rhoであるからτ≠t∗\tau\neq t_*である。t∗t_*とτ\tauを端点とする閉区間をJ′′J''とする。J′′J''のτ\tau以外の点ttは∣t−t∗∣<∣τ−t∗∣|t-t_*|<|\tau-t_*|を満たすのでEEに属さず、u(t)<ρu(t)<\rhoである。uuの連続性からu(τ)≤ρu(\tau)\leq\rhoであり、τ∈E\tau\in Eと合わせてu(τ)=ρu(\tau)=\rhoである。従ってJ′′J''上で二つの解w,y\boldsymbol w,\boldsymbol yのグラフはKKに含まれる。定理 2.1を区間J′′J''、S=K∩(J′′×Rd)S=K\cap(J''\times\R^d)、t0=t∗t_0=t_*に適用すると

u(τ)≤exp⁡(∫J′′ℓ(s) ds)(δ+∫J′′r(s) ds)≤Q<ρu(\tau)\leq \exp\left(\int_{J''}\ell(s)\,ds\right) \left(\delta+\int_{J''}r(s)\,ds\right) \leq Q<\rho

となり、u(τ)=ρu(\tau)=\rhoと両立しない。従ってEEは空であり、J′J'上で両解のグラフはKKに含まれる。定理 2.1を区間J′J'、S=K∩(J′×Rd)S=K\cap(J'\times\R^d)、t0=t∗t_0=t_*に適用すると、J′J'上でu≤Qu\leq Qを得る。▨

δ≤Q<ρ\delta\leq Q<\rhoより(t∗,η)∈K⊆Ω(t_*,\boldsymbol\eta)\in K\subseteq\Omegaである。§E13.6 定理 2.1により、z(t∗)=η\boldsymbol z(t_*)=\boldsymbol\etaを満たすggの最大絶対連続解z:(α,β)→Rd\boldsymbol z:(\alpha,\beta)\to\R^dを取る。α<t∗<β\alpha<t_*<\betaである。

β≤b\beta\leq bと仮定する。任意のt∈[t∗,β)t\in[t_*,\beta)に対してt<c<βt<c<\betaとなるccを取ると、主張 3.1.1をJ′=[t∗,c]J'=[t_*,c]とz∣J′\boldsymbol z|_{J'}に適用して(t,z(t))∈K(t,\boldsymbol z(t))\in Kを得る。β\betaは有限で、K⊆ΩK\subseteq\Omegaはコンパクトであるから、これは§E13.6 系 4.3の右端点の結論に反する。従ってβ>b\beta>bである。α≥a\alpha\geq aと仮定すると、任意のt∈(α,t∗]t\in(\alpha,t_*]に対してα<c<t\alpha<c<tとなるccを取り、主張 3.1.1をJ′=[c,t∗]J'=[c,t_*]に適用して(t,z(t))∈K(t,\boldsymbol z(t))\in Kを得る。これは同じ系の左端点の結論に反するので、α<a\alpha<aである。従ってJ⊆(α,β)J\subseteq(\alpha,\beta)であり、主張 3.1.1をJ′=JJ'=Jとz∣J\boldsymbol z|_Jに適用すると結論の不等式を得る。

z~\widetilde{\boldsymbol z}を、z~(t∗)=η\widetilde{\boldsymbol z}(t_*)=\boldsymbol\etaを満たすJJ上の絶対連続解とする。主張 3.1.1により、z~\widetilde{\boldsymbol z}とz∣J\boldsymbol z|_JのグラフはKKに含まれる。定理 2.1をS=KS=K上の二つの右辺g,gg,g、二解z∣J,z~\boldsymbol z|_J,\widetilde{\boldsymbol z}、t0=t∗t_0=t_*、r=0r=0に適用する。初期値の差は零であるから、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、t∗=at_*=a、η=ηn\boldsymbol\eta=\boldsymbol\eta_nとして適用すると、共通区間JJ上の一意な解とQnQ_nによる評価を得る。▨

4 線形増大条件と区間全体での解

定理 4.1.d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\|\cdot\|を固定する。I⊆RI\subseteq\Rを内部が空でない区間とし、f:I×Rd→Rdf:I\times\R^d\to\R^dを写像とする。IIに含まれる長さが正の各有限閉区間JJとRd\R^dの各閉球BBについて、ffのJ×BJ\times Bへの制限は Carathéodory 条件を満たし、時間について可積分な Lipschitz 係数をもつとする。さらに、非負関数p,q∈Lloc1(I)p,q\in L^1_{\mathrm{loc}}(I)と Lebesgue 零集合N⊆IN\subseteq Iが存在し、すべてのt∈I∖Nt\in I\setminus Nとu∈Rd\boldsymbol u\in\R^dについて

∥f(t,u)∥≤p(t)+q(t)∥u∥\|f(t,\boldsymbol u)\|\leq p(t)+q(t)\|\boldsymbol u\|

が成り立つとする。このとき、任意のt∗∈It_*\in Iとη∈Rd\boldsymbol\eta\in\R^dに対して、IIのほとんど至る所でx′=f(t,x)\boldsymbol x'=f(t,\boldsymbol x)を満たし、x(t∗)=η\boldsymbol x(t_*)=\boldsymbol\etaを満たすx∈ACloc(I;Rd)\boldsymbol x\in AC_{\mathrm{loc}}(I;\R^d)が一意に存在する。この解はすべてのt∈It\in Iについて

∥x(t)∥≤(∥η∥+∣∫t∗tp(s) ds∣)exp⁡∣∫t∗tq(s) ds∣\|\boldsymbol x(t)\|\leq \left(\|\boldsymbol\eta\|+\left|\int_{t_*}^tp(s)\,ds\right|\right) \exp\left|\int_{t_*}^tq(s)\,ds\right|

を満たす。

証明.JJを、t∗t_*を含みIIに含まれる長さが正の有限閉区間とする。R×Rd\R\times\R^d上の写像fJf_Jを、t∈Jt\in JではfJ(t,u)=f(t,u)f_J(t,\boldsymbol u)=f(t,\boldsymbol u)、t∉Jt\notin JではfJ(t,u)=0f_J(t,\boldsymbol u)=\boldsymbol 0と定める。

長さが正の有限閉区間J′J'と閉球BBを取り、E=J′∩JE=J'\cap Jとおく。EEの長さが正ならば、EEはIIに含まれる長さが正の有限閉区間であるから、E×BE\times B上の Carathéodory 条件の零集合NEN_E、支配関数mEm_E、可積分な Lipschitz 係数ℓE\ell_Eが存在する。各u∈B\boldsymbol u\in Bについてt↦fJ(t,u)t\mapsto f_J(t,\boldsymbol u)はEE上でf(t,u)f(t,\boldsymbol u)、J′∖EJ'\setminus E上で0\boldsymbol 0であるから可測であり、J′×BJ'\times B上で零集合NEN_E、支配関数mE1Em_E\mathbf 1_E、Lipschitz 係数ℓE1E\ell_E\mathbf 1_Eについて条件が成り立つ。EEの長さが零ならば、EEは空集合または一点であり、J′∖EJ'\setminus E上でfJ=0f_J=\boldsymbol 0であるから、零集合EE、支配関数00、Lipschitz 係数00について条件が成り立つ。従ってfJf_JはR×Rd\R\times\R^d上で Carathéodory 条件と時間について可積分な Lipschitz 係数の条件を局所的に満たす。また、pJ=p1Jp_J=p\mathbf 1_J、qJ=q1Jq_J=q\mathbf 1_JはL1(R)L^1(\R)に属し、t∉Nt\notin Nならばすべてのu∈Rd\boldsymbol u\in\R^dについて∥fJ(t,u)∥≤pJ(t)+qJ(t)∥u∥\|f_J(t,\boldsymbol u)\|\leq p_J(t)+q_J(t)\|\boldsymbol u\|である。

§E13.6 定理 2.1により、x(t∗)=η\boldsymbol x(t_*)=\boldsymbol\etaを満たすfJf_Jの最大絶対連続解x:(α,β)→Rd\boldsymbol x:(\alpha,\beta)\to\R^dを取る。t∗<T<βt_*<T<\betaを取る。x\boldsymbol xは[t∗,T][t_*,T]上で絶対連続であり、ほとんど至る所でx′=fJ(t,x)\boldsymbol x'=f_J(t,\boldsymbol x)を満たすから、§E13.2 定理 2.5と§E13.4 系 3.1により、t∗≤t≤Tt_*\leq t\leq Tで

∥x(t)∥≤∥η∥+∫t∗t∥x′(s)∥ ds≤∥η∥+∫t∗t(qJ(s)∥x(s)∥+pJ(s)) ds\|\boldsymbol x(t)\|\leq\|\boldsymbol\eta\|+\int_{t_*}^t\|\boldsymbol x'(s)\|\,ds \leq\|\boldsymbol\eta\|+\int_{t_*}^t\bigl(q_J(s)\|\boldsymbol x(s)\|+p_J(s)\bigr)\,ds

である。定理 1.1の前半の評価を[t∗,T][t_*,T]上でA=∥η∥A=\|\boldsymbol\eta\|、ℓ=qJ\ell=q_J、r=pJr=p_Jとして用い、t∗≤s≤tt_*\leq s\leq tでexp⁡(∫stqJ)≤exp⁡(∫t∗tqJ)\exp(\int_s^tq_J)\leq\exp(\int_{t_*}^tq_J)であることを用いると、

∥x(t)∥≤(∥η∥+∫t∗tpJ(s) ds)exp⁡(∫t∗tqJ(s) ds)(t∗≤t≤T)\|\boldsymbol x(t)\|\leq \left(\|\boldsymbol\eta\|+\int_{t_*}^tp_J(s)\,ds\right) \exp\left(\int_{t_*}^tq_J(s)\,ds\right) \qquad(t_*\leq t\leq T)

を得る。TTは任意であるから、この評価はt∗≤t<βt_*\leq t<\betaで成り立つ。α<T<t∗\alpha<T<t_*を取り、T≤t≤t∗T\leq t\leq t_*で成り立つ

∥x(t)∥≤∥η∥+∫tt∗(qJ(s)∥x(s)∥+pJ(s)) ds\|\boldsymbol x(t)\|\leq\|\boldsymbol\eta\|+\int_t^{t_*}\bigl(q_J(s)\|\boldsymbol x(s)\|+p_J(s)\bigr)\,ds

に定理 1.1の後半の評価を用いると、α<t≤t∗\alpha<t\leq t_*で積分の範囲を[t,t∗][t,t_*]に替えた同じ形の評価を得る。従って、C=(∥η∥+∫Jp)exp⁡(∫Jq)C=(\|\boldsymbol\eta\|+\int_Jp)\exp(\int_Jq)とおくと、すべてのt∈(α,β)t\in(\alpha,\beta)で∥x(t)∥≤C\|\boldsymbol x(t)\|\leq Cである。

β<∞\beta<\inftyと仮定すると、[t∗,β)[t_*,\beta)上のx\boldsymbol xのグラフはコンパクト集合[t∗,β]×{u:∥u∥≤C}[t_*,\beta]\times\{\boldsymbol u:\|\boldsymbol u\|\leq C\}に含まれる。これは§E13.6 系 4.3の右端点の結論に反する。α>−∞\alpha>-\inftyと仮定すると、(α,t∗](\alpha,t_*]上のグラフはコンパクト集合[α,t∗]×{u:∥u∥≤C}[\alpha,t_*]\times\{\boldsymbol u:\|\boldsymbol u\|\leq C\}に含まれる。これは同じ系の左端点の結論に反する。従って(α,β)=R(\alpha,\beta)=\Rである。JJ上でfJ=ff_J=fであるから、xJ=x∣J\boldsymbol x_J=\boldsymbol x|_JはJJ上の絶対連続関数で、JJのほとんど至る所でxJ′=f(t,xJ)\boldsymbol x_J'=f(t,\boldsymbol x_J)、xJ(t∗)=η\boldsymbol x_J(t_*)=\boldsymbol\etaを満たし、t∈Jt\in Jで結論の評価を満たす。

w:J→Rd\boldsymbol w:J\to\R^dを、JJのほとんど至る所でw′=f(t,w)\boldsymbol w'=f(t,\boldsymbol w)を満たしw(t∗)=η\boldsymbol w(t_*)=\boldsymbol\etaを満たす絶対連続関数とする。xJ\boldsymbol x_Jとw\boldsymbol wは連続であるから、両者の値を含む原点中心の閉球BBが存在する。J×BJ\times B上の可積分な Lipschitz 係数ℓB\ell_Bと零集合NBN_Bを取り、定理 2.1をS=J×BS=J\times B上の二つの右辺f,ff,f、二解xJ,w\boldsymbol x_J,\boldsymbol w、t0=t∗t_0=t_*、ℓ=ℓB\ell=\ell_B、r=0r=0に適用する。初期値の差は零であるから、w=xJ\boldsymbol w=\boldsymbol x_Jである。

t∗t_*を含みIIに含まれる長さが正の有限閉区間の全体をJ\mathcal Jとする。J1,J2∈JJ_1,J_2\in\mathcal JならばJ1∪J2∈JJ_1\cup J_2\in\mathcal Jであり、xJ1∪J2\boldsymbol x_{J_1\cup J_2}のJ1J_1、J2J_2への制限は、前段の一意性によりそれぞれxJ1\boldsymbol x_{J_1}、xJ2\boldsymbol x_{J_2}に等しい。従って、t∈J∈Jt\in J\in\mathcal Jに対してx(t)=xJ(t)\boldsymbol x(t)=\boldsymbol x_J(t)と定めると、写像x:I→Rd\boldsymbol x:I\to\R^dがJJの取り方によらず定まる。IIに含まれる任意の有限閉区間はJ\mathcal Jのある元に含まれるから、x∈ACloc(I;Rd)\boldsymbol x\in AC_{\mathrm{loc}}(I;\R^d)であり、すべてのt∈It\in Iで結論の評価が成り立つ。IIはJ\mathcal Jの可算個の元の和であるから、x′=f(t,x)\boldsymbol x'=f(t,\boldsymbol x)はIIのほとんど至る所で成り立つ。同じ条件を満たすII上の関数は、各J∈JJ\in\mathcal Jへの制限が前段の一意性によりxJ\boldsymbol x_Jに等しいから、x\boldsymbol xに等しい。▨

系 4.2.K∈{R,C}\mathbb K\in\{\R,\C\}、d∈N≥1d\in\NNとし、Kd\mathbb K^dにノルム∥⋅∥\|\cdot\|を固定し、d×dd\times d行列にはその作用素ノルムを用いる。I⊆RI\subseteq\Rを内部が空でない区間とし、

A∈Lloc1(I;Kd×d),b∈Lloc1(I;Kd)A\in L^1_{\mathrm{loc}}(I;\mathbb K^{d\times d}),\qquad \boldsymbol b\in L^1_{\mathrm{loc}}(I;\mathbb K^d)

とする。このとき、任意のt∗∈It_*\in Iとη∈Kd\boldsymbol\eta\in\mathbb K^dに対して、IIのほとんど至る所でx′=A(t)x+b(t)\boldsymbol x'=A(t)\boldsymbol x+\boldsymbol b(t)を満たし、x(t∗)=η\boldsymbol x(t_*)=\boldsymbol\etaを満たすx∈ACloc(I;Kd)\boldsymbol x\in AC_{\mathrm{loc}}(I;\mathbb K^d)が一意に存在する。この解はすべてのt∈It\in Iについて

x(t)=η+∫t∗t(A(s)x(s)+b(s)) ds,∥x(t)∥≤(∥η∥+∣∫t∗t∥b(s)∥ ds∣)exp⁡∣∫t∗t∥A(s)∥ ds∣\boldsymbol x(t)=\boldsymbol\eta+\int_{t_*}^t\bigl(A(s)\boldsymbol x(s)+\boldsymbol b(s)\bigr)\,ds, \qquad \|\boldsymbol x(t)\|\leq \left(\|\boldsymbol\eta\|+\left|\int_{t_*}^t\|\boldsymbol b(s)\|\,ds\right|\right) \exp\left|\int_{t_*}^t\|A(s)\|\,ds\right|

を満たす。t<t∗t<t_*の積分は向き付き積分とする。

証明.K=C\mathbb K=\Cのときは、Cd\C^dを実部と虚部によりR2d\R^{2d}と同一視する。固定したノルムはR2d\R^{2d}のノルムであり、各ttについてu↦A(t)u\boldsymbol u\mapsto A(t)\boldsymbol uは実線形である。A,bA,\boldsymbol bの有限値の可測な代表元を取り、f(t,u)=A(t)u+b(t)f(t,\boldsymbol u)=A(t)\boldsymbol u+\boldsymbol b(t)とおく。作用素ノルムはKd×d\mathbb K^{d\times d}のノルムであるから、t↦∥A(t)∥t\mapsto\|A(t)\|は可測であり、§E10.1 命題 2.2によりIIに含まれる各有限閉区間上で可積分である。

IIに含まれる長さが正の有限閉区間JJと、中心c\boldsymbol c、半径κ\kappaの閉球BBを取る。各u\boldsymbol uについてt↦f(t,u)t\mapsto f(t,\boldsymbol u)は可測であり、各ttについてu↦f(t,u)\boldsymbol u\mapsto f(t,\boldsymbol u)は連続であって、u,v∈B\boldsymbol u,\boldsymbol v\in Bならば

∥f(t,u)∥≤∥A(t)∥(∥c∥+κ)+∥b(t)∥,∥f(t,u)−f(t,v)∥≤∥A(t)∥∥u−v∥\|f(t,\boldsymbol u)\|\leq\|A(t)\|(\|\boldsymbol c\|+\kappa)+\|\boldsymbol b(t)\|, \qquad \|f(t,\boldsymbol u)-f(t,\boldsymbol v)\|\leq\|A(t)\|\|\boldsymbol u-\boldsymbol v\|

である。右辺の係数はJJ上で可積分であるから、ffのJ×BJ\times Bへの制限は Carathéodory 条件を満たし、時間について可積分な Lipschitz 係数∥A∥\|A\|をもつ。すべてのt∈It\in Iとu\boldsymbol uについて∥f(t,u)∥≤∥b(t)∥+∥A(t)∥∥u∥\|f(t,\boldsymbol u)\|\leq\|\boldsymbol b(t)\|+\|A(t)\|\|\boldsymbol u\|である。定理 4.1をp=∥b∥p=\|\boldsymbol b\|、q=∥A∥q=\|A\|、N=∅N=\varnothingとして適用すると、解の存在、一意性、ノルムの評価を得る。x\boldsymbol xは局所絶対連続で、その導関数はほとんど至る所でAx+bA\boldsymbol x+\boldsymbol bに等しいから、§E13.2 系 2.6により積分表示が成り立つ。▨

例 4.3.f(t,u)=u2f(t,u)=u^2とおく。ffはR×R\R\times\R上で連続であり、状態変数について局所 Lipschitz 連続であるから、§E13.4 補題 4.3 (2)により、長さが正の各有限閉区間JJと長さが正の各閉区間BBについて、ffのJ×BJ\times Bへの制限は Carathéodory 条件を満たし、時間について可積分な Lipschitz 係数をもつ。一点BBへの制限は、BBを含む長さが正の閉区間への制限から同じ条件を受け継ぐ。一方、非負関数p,qp,qとttをどのように取っても、十分大きい∣u∣|u|ではu2>p(t)+q(t)∣u∣u^2>p(t)+q(t)|u|であるから、ffは定理 4.1の線形増大条件を満たさない。

R\R上の絶対連続解xxがx(0)=1x(0)=1を満たすと仮定する。§E13.4 補題 4.3 (3)によりxxは古典解であり、§E13.6 系 5.2 (2)によりxxはy′=y2y'=y^2、y(0)=1y(0)=1の最大解の制限である。§E13.6 例 6.1によりこの最大解の最大存在区間は(−∞,1)(-\infty,1)であるから、R⊆(−∞,1)\R\subseteq(-\infty,1)となり、両者は両立しない。従って、I=RI=\R上の解は存在せず、定理 4.1の結論は線形増大条件を外すと成り立たない。

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

系 5.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、L≥0L\ge0とし、

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

を連続写像とする。すべてのt∈Jt\in Jとu,v∈D\boldsymbol u,\boldsymbol v\in Dについて

∥g(t,u)−g(t,v)∥≤L∥u−v∥\lVert g(t,\boldsymbol u)-g(t,\boldsymbol v)\rVert \le L\lVert\boldsymbol u-\boldsymbol v\rVert

が成り立ち、

ε=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を、すべてのt∈Jt\in Jで

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)

である。

証明. 一階連続微分可能な両解は、§E13.2 系 2.7 (2)により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である。従って残りの二式も成り立つ。▨

系 5.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)の解とする。ρ>0\rho>0とし、

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

とおく。K⊆I×DK\subseteq I\times Dと仮定する。定数L,ε≥0L,\varepsilon\ge0が存在し、任意のt∈Jt\in Jとu,v∈D\boldsymbol u,\boldsymbol v\in Dに対して

∥g(t,u)−g(t,v)∥≤L∥u−v∥,∥f(t,u)−g(t,u)∥≤ε\lVert g(t,\boldsymbol u)-g(t,\boldsymbol v)\rVert \le L\lVert\boldsymbol u-\boldsymbol v\rVert,\qquad \lVert f(t,\boldsymbol u)-g(t,\boldsymbol u)\rVert \le\varepsilon

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

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

を満たすとする。このとき、初期値問題

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<ρ\max_{t\in J}\lVert\boldsymbol z(t)-\boldsymbol y(t)\rVert\le Q<\rho

が成り立つ。

証明.§E13.4 補題 4.3 (2)によりggはI×DI\times D上で Carathéodory 条件と時間について可積分な Lipschitz 係数の条件を局所的に満たす。§E13.2 系 2.7 (2)により古典解y\boldsymbol yはJJ上で絶対連続であり、絶対連続解である。定理 3.1にΩ=I×D\Omega=I\times D、t∗=at_*=a、ℓ(t)=L\ell(t)=L、誤差関数として定数ε\varepsilonを代入する。このとき

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

であるから、JJ上の絶対連続解と評価を得る。(a,η)∈K⊆I×D(a,\boldsymbol\eta)\in K\subseteq I\times Dであるから、得られた解は§E13.4 補題 4.3 (3)により古典解である。▨

系 5.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を代入する。§E13.2 系 2.7 (2)により古典解はJJ上で絶対連続であるから、同定理の評価が成り立つ。▨

系 5.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をyn′=fn(t,yn)\boldsymbol y_n'=f_n(t,\boldsymbol y_n)の解とし、y:J→D\boldsymbol y:J\to Dをy′=f(t,y)\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に対して系 5.1をg=fng=f_n、z=yn\boldsymbol z=\boldsymbol y_n、ε=εn\varepsilon=\varepsilon_nとして用いると、t∈Jt\in Jで

∥yn(t)−y(t)∥≤eL(t−a)δn+εn∫ateL(t−s) ds≤eL(b−a)(δn+(b−a)εn)\lVert\boldsymbol y_n(t)-\boldsymbol y(t)\rVert \leq e^{L(t-a)}\delta_n+\varepsilon_n\int_a^te^{L(t-s)}\,ds \leq e^{L(b-a)}\bigl(\delta_n+(b-a)\varepsilon_n\bigr)

である。右辺は零へ収束する。▨

例 5.5.系 5.4は、各解yn\boldsymbol y_nが区間JJ全体で定まりDDに値をとることを仮定している。この仮定はδn→0\delta_n\to0、εn→0\varepsilon_n\to0からは従わない。

J=[0,1]J=[0,1]、D=[0,∞)D=[0,\infty)、f(t,u)=0f(t,u)=0、fn(t,u)=−1/nf_n(t,u)=-1/nとし、y=0y=0とする。f,fnf,f_nは連続で、状態変数について共通の Lipschitz 定数00をもち、εn=1/n→0\varepsilon_n=1/n\to0である。一階連続微分可能なyn:J→Dy_n:J\to Dがyn′=fn(t,yn)y_n'=f_n(t,y_n)、yn(0)=0y_n(0)=0を満たすと仮定すると、平均値の定理によりyn(t)=−t/ny_n(t)=-t/nであり、yn(1)=−1/n∉Dy_n(1)=-1/n\notin Dとなってyn(J)⊆Dy_n(J)\subseteq Dと両立しない。従って、δn=0\delta_n=0となる初期値yn(0)=0y_n(0)=0に対して、JJ上でDDに値をとる解は存在しない。

この例では、どのρ>0\rho>0についてもyyの閉管J×[−ρ,ρ]J\times[-\rho,\rho]はJ×DJ\times Dに含まれない。閉管が右辺の定義域に含まれるという仮定の下では、系 3.2と系 5.2により、JJ全体で解が存在する。

6 計算例

例 6.1.T>0T>0、λ≥0\lambda\ge0、c,y0,z0∈Rc,y_0,z_0\in\Rとし、

y′=λy,z′=λz+cy'=\lambda y,\qquad z'=\lambda z+c

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

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

が与えられる。

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

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

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

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

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

となる。

例 6.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

であり、系 3.2の上界∫01rn(s) ds=1/n\int_0^1r_n(s)\,ds=1/nはこの例で等号となる。

例 6.3.f(t,u)=2∣u∣f(t,u)=2\sqrt{|u|}とおく。ffはR×R\R\times\R上で連続であるが、u≠0u\neq0について∣f(t,u)−f(t,0)∣/∣u∣=2/∣u∣|f(t,u)-f(t,0)|/|u|=2/\sqrt{|u|}であるから、u=0u=0を含むどの[0,1]×[−ρ,ρ][0,1]\times[-\rho,\rho]上でも状態変数について Lipschitz 連続でない。

y=0y=0はy(0)=0y(0)=0を満たす古典解である。ε>0\varepsilon>0に対してyε(t)=(ε+t)2y_\varepsilon(t)=(\sqrt\varepsilon+t)^2とおくと、[0,1][0,1]上でyε>0y_\varepsilon>0であり、

yε′(t)=2(ε+t)=2yε(t),yε(0)=εy_\varepsilon'(t)=2(\sqrt\varepsilon+t)=2\sqrt{y_\varepsilon(t)},\qquad y_\varepsilon(0)=\varepsilon

である。初期値の差はε\varepsilonであるが、

max⁡0≤t≤1∣yε(t)−y(t)∣=(1+ε)2>1\max_{0\leq t\leq1}|y_\varepsilon(t)-y(t)|=(1+\sqrt\varepsilon)^2>1

である。従って、ε→0\varepsilon\to0のときyεy_\varepsilonはyyへ一様収束せず、max⁡[0,1]∣yε−y∣≤Cε\max_{[0,1]}|y_\varepsilon-y|\leq C\varepsilonを満たす定数CCも存在しない。系 5.3の結論は、状態変数についての Lipschitz 条件を外すと成り立たない。なお、0≤t≤10\leq t\leq1で0≤yε(t)−t2=ε+2ε t≤ε+2ε0\leq y_\varepsilon(t)-t^2=\varepsilon+2\sqrt\varepsilon\,t\leq\varepsilon+2\sqrt\varepsilonであるから、yεy_\varepsilonは、y(0)=0y(0)=0を満たすもう一つの古典解t↦t2t\mapsto t^2へ一様収束する。

7 演習

問題 7.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

である。これは系 5.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であるから、tnt_nにおける差が11であることは、固定された[0,T][0,T]上での一様収束と矛盾しない。▨

問題 7.2 (Grönwall の不等式の変形と等号の例).a<ba<b、L≥0L\ge0とし、φ:[a,b]→[0,∞)\varphi:[a,b]\to[0,\infty)を単調非減少な連続関数、u:[a,b]→[0,∞)u:[a,b]\to[0,\infty)を連続関数とする。すべてのt∈[a,b]t\in[a,b]についてu(t)≤φ(t)+L∫atu(s) dsu(t)\le\varphi(t)+L\int_a^tu(s)\,dsが成り立つならば、u(t)≤φ(t)eL(t−a)u(t)\le\varphi(t)e^{L(t-a)}であることを示せ。また、a=0a=0、b=1b=1、L=1L=1、φ(t)=1−t\varphi(t)=1-t、u≡1u\equiv1が仮定の不等式を満たすが結論を満たさないことを示し、単調性の仮定を外すことができないことを確かめよ。

次に、λ∈R\lambda\in\R、T>0T>0とし、y,z:[0,T]→Ry,z:[0,T]\to\Rをy′=λyy'=\lambda yの二つの古典解で、δ=∣y(T)−z(T)∣>0\delta=|y(T)-z(T)|>0を満たすものとする。定理 2.1をt0=Tt_0=Tとして用い、0≤t≤T0\le t\le Tで∣y(t)−z(t)∣≤δe∣λ∣(T−t)|y(t)-z(t)|\le\delta e^{|\lambda|(T-t)}であることを導け。また、y−zy-zを求め、この評価がすべてのt∈[0,T)t\in[0,T)で等号となるための必要十分条件がλ≤0\lambda\le0であることを示せ。

解答.

t0∈(a,b]t_0\in(a,b]を固定する。a≤t≤t0a\le t\le t_0ではφ(t)≤φ(t0)\varphi(t)\le\varphi(t_0)であるから、u(t)≤φ(t0)+L∫atu(s) dsu(t)\le\varphi(t_0)+L\int_a^tu(s)\,dsが[a,t0][a,t_0]の上で成り立つ。系 1.3を区間[a,t0][a,t_0]とA=φ(t0)≥0A=\varphi(t_0)\ge0に用いると、u(t0)≤φ(t0)eL(t0−a)u(t_0)\le\varphi(t_0)e^{L(t_0-a)}である。t0=at_0=aでは仮定からu(a)≤φ(a)u(a)\le\varphi(a)である。ゆえに結論が成り立つ。

u≡1u\equiv1、φ(t)=1−t\varphi(t)=1-tとすると、φ(t)+∫0tu(s) ds=1−t+t=1=u(t)\varphi(t)+\int_0^tu(s)\,ds=1-t+t=1=u(t)であり、仮定の不等式は成り立つ。一方、t=1t=1ではφ(1)e1=0<1=u(1)\varphi(1)e^{1}=0<1=u(1)であるから結論は成り立たない。φ\varphiはφ(1)=0\varphi(1)=0を満たす非負の連続関数であるが、単調非減少でない。従って、単調性の仮定を外すと結論は成り立たない。

f(t,u)=λuf(t,u)=\lambda uはS=[0,T]×RS=[0,T]\times\R上ですべてのttについて∣f(t,u)−f(t,v)∣≤∣λ∣∣u−v∣|f(t,u)-f(t,v)|\le|\lambda||u-v|を満たす。§E13.2 系 2.7 (2)により古典解y,zy,zは[0,T][0,T]上で絶対連続であるから、定理 2.1をSS、g=fg=f、t0=Tt_0=T、ℓ=∣λ∣\ell=|\lambda|、r=0r=0に適用すると、0≤t≤T0\le t\le Tで

∣y(t)−z(t)∣≤δexp⁡(∫tT∣λ∣ ds)=δe∣λ∣(T−t)|y(t)-z(t)|\le\delta\exp\left(\int_t^T|\lambda|\,ds\right)=\delta e^{|\lambda|(T-t)}

である。w=y−zw=y-zはw′=λww'=\lambda wを満たすから、(e−λtw(t))′=e−λt(w′(t)−λw(t))=0(e^{-\lambda t}w(t))'=e^{-\lambda t}(w'(t)-\lambda w(t))=0であり、e−λtw(t)e^{-\lambda t}w(t)は定数である。従ってw(t)=w(T)e−λ(T−t)w(t)=w(T)e^{-\lambda(T-t)}、∣w(t)∣=δe−λ(T−t)|w(t)|=\delta e^{-\lambda(T-t)}である。0≤t<T0\le t<TではT−t>0T-t>0であるから、δe−λ(T−t)=δe∣λ∣(T−t)\delta e^{-\lambda(T-t)}=\delta e^{|\lambda|(T-t)}は−λ=∣λ∣-\lambda=|\lambda|、すなわちλ≤0\lambda\le0と同値である。ゆえに、評価がすべてのt∈[0,T)t\in[0,T)で等号となるための必要十分条件はλ≤0\lambda\le0である。λ>0\lambda>0のとき、時刻TTから過去へ遡ると二つの解の差は縮むが、上界はeλ(T−t)e^{\lambda(T-t)}倍に増える。▨

参考文献

  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 不等式、初期値に関する安定性、右辺の摂動評価を参考にした。

前提記事