§E13.5Peano の存在定理

最終更新

常微分方程式の初期値問題で、右辺が連続であり、状態変数について局所 Lipschitz 連続である場合には、Picard の逐次近似によって局所解の存在と一意性が得られる。しかし、右辺が連続であっても Lipschitz 連続であるとは限らず、たとえばy′=∣y∣y'=\sqrt{\lvert y\rvert}、y(0)=0y(0)=0では解が一意でない。一方、右辺が連続であれば、Lipschitz 連続でない場合にも、Peano の存在定理により解が存在する。Peano の存在定理は、初期値問題に対する最も基本的な存在定理の一つである。Peano の存在定理は、時間について可測で状態について連続であり、大きさが時間について局所可積分な関数で支配される右辺に対する、絶対連続な解の存在定理へ拡張される。本記事では、この拡張を証明して Peano の存在定理を導き、仮定を外すと解が存在しない例、スカラー方程式の最大の解と最小の解、および微分不等式による比較を扱う。

1 Tonelli の近似解と局所解の存在

補題 1.1.d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\|\cdot\|を固定する。t0∈Rt_0\in\R、y0∈Rd\boldsymbol y_0\in\R^d、a,b>0a,b>0とし、

J=[t0−a,t0+a],B={u∈Rd:∥u−y0∥≤b}J=[t_0-a,t_0+a],\qquad B=\{\boldsymbol u\in\R^d:\|\boldsymbol u-\boldsymbol y_0\|\le b\}

とおく。写像f:J×B→Rdf:J\times B\to\R^dは、Lebesgue 零集合N⊆JN\subseteq Jと非負関数m∈L1(J)m\in L^1(J)について§E13.4 定義 2.1の Carathéodory 条件を満たすとする。0<h≤a0<h\le aは

∫t0−ht0m(s) ds≤b,∫t0t0+hm(s) ds≤b\int_{t_0-h}^{t_0}m(s)\,ds\le b,\qquad \int_{t_0}^{t_0+h}m(s)\,ds\le b

を満たすとし、Ih=[t0−h,t0+h]I_h=[t_0-h,t_0+h]とおく。r≥0r\ge0に対して

ω(r)=sup⁡{∫stm(τ) dτ:s,t∈J, s≤t≤s+r}\omega(r)=\sup\left\{\int_s^tm(\tau)\,d\tau:s,t\in J,\ s\le t\le s+r\right\}

とおく。n∈N≥1n\in\NNに対してδn=h/n\delta_n=h/nとし、σn:Ih→Ih\sigma_n:I_h\to I_hを

σn(t)={t−δn,t0+δn≤t≤t0+h,t0,∣t−t0∣≤δn,t+δn,t0−h≤t≤t0−δn\sigma_n(t)= \begin{cases} t-\delta_n,&t_0+\delta_n\le t\le t_0+h,\\ t_0,&\lvert t-t_0\rvert\le\delta_n,\\ t+\delta_n,&t_0-h\le t\le t_0-\delta_n \end{cases}

で定める。このとき、次が成り立つ。

  1. 各n∈N≥1n\in\NNについて、連続写像yn:Ih→B\boldsymbol y_n:I_h\to Bであって、s↦f(s,yn(s))s\mapsto f(s,\boldsymbol y_n(s))がIhI_h上で可積分であり、すべてのt∈Iht\in I_hに対して yn(t)=y0+∫t0σn(t)f(s,yn(s)) ds\boldsymbol y_n(t)=\boldsymbol y_0+\int_{t_0}^{\sigma_n(t)}f(s,\boldsymbol y_n(s))\,ds を満たすものが、ただ一つ存在する。
  2. すべてのn∈N≥1n\in\NNとs,t∈Ihs,t\in I_hに対して、∥yn(t)−yn(s)∥≤ω(∣t−s∣)\|\boldsymbol y_n(t)-\boldsymbol y_n(s)\|\le\omega(\lvert t-s\rvert)である。
  3. すべてのn∈N≥1n\in\NNとt∈Iht\in I_hに対して ∥yn(t)−y0−∫t0tf(s,yn(s)) ds∥≤ω(δn)\left\|\boldsymbol y_n(t)-\boldsymbol y_0-\int_{t_0}^tf(s,\boldsymbol y_n(s))\,ds\right\|\le\omega(\delta_n) である。
  4. ω\omegaは単調非減少であり、r↓0r\downarrow0のときω(r)→0\omega(r)\to0である。

証明.I′⊆JI'\subseteq Jをt0t_0を含む閉区間とし、x:I′→B\boldsymbol x:I'\to Bを連続写像とする。x\boldsymbol xをJ∖I′J\setminus I'上でI′I'の近い方の端点での値に等しく延長すると、JJからBBへの連続写像を得る。この延長に§E13.4 補題 2.3を適用すると、s↦f(s,x(s))s\mapsto f(s,\boldsymbol x(s))はI′I'上で可積分であり、s∈I′∖Ns\in I'\setminus Nで∥f(s,x(s))∥≤m(s)\|f(s,\boldsymbol x(s))\|\le m(s)である。

主張 (1)を示す。nnを固定し、k∈{1,…,n}k\in\{1,\ldots,n\}に対してI(k)=[t0−kδn,t0+kδn]I^{(k)}=[t_0-k\delta_n,t_0+k\delta_n]とおく。I(n)=IhI^{(n)}=I_hである。σn\sigma_nの定義から、各t∈Iht\in I_hについてσn(t)\sigma_n(t)はt0t_0とttの間にあり、∣σn(t)−t∣≤δn\lvert\sigma_n(t)-t\rvert\le\delta_nである。したがってσn(I(k))⊆I(k)\sigma_n(I^{(k)})\subseteq I^{(k)}であり、k<nk<nならばσn(I(k+1))⊆I(k)\sigma_n(I^{(k+1)})\subseteq I^{(k)}である。kkについての帰納法により、次の主張(Sk)(\mathrm S_k)を示す: 連続写像x:I(k)→B\boldsymbol x:I^{(k)}\to Bであって、すべてのt∈I(k)t\in I^{(k)}でx(t)=y0+∫t0σn(t)f(s,x(s)) ds\boldsymbol x(t)=\boldsymbol y_0+\int_{t_0}^{\sigma_n(t)}f(s,\boldsymbol x(s))\,dsを満たすものがただ一つ存在する。I(1)I^{(1)}上ではσn≡t0\sigma_n\equiv t_0であるから、この等式はx≡y0\boldsymbol x\equiv\boldsymbol y_0と同値であり、(S1)(\mathrm S_1)が成り立つ。

k<nk<nとし、(Sk)(\mathrm S_k)の写像をxk\boldsymbol x_k、gk(s)=f(s,xk(s))\boldsymbol g_k(s)=f(s,\boldsymbol x_k(s))(s∈I(k)s\in I^{(k)})とおく。gk\boldsymbol g_kはI(k)I^{(k)}上で可積分である。t∈I(k+1)t\in I^{(k+1)}に対して

x(t)=y0+∫t0σn(t)gk(s) ds\boldsymbol x(t)=\boldsymbol y_0+\int_{t_0}^{\sigma_n(t)}\boldsymbol g_k(s)\,ds

と定める。§E13.2 補題 1.3 (3)によりt′↦∫t0t′gkt'\mapsto\int_{t_0}^{t'}\boldsymbol g_kはI(k)I^{(k)}上で連続であり、σn\sigma_nは連続であるから、x\boldsymbol xは連続である。§E13.4 系 3.1により、∥x(t)−y0∥\|\boldsymbol x(t)-\boldsymbol y_0\|はt0t_0とσn(t)\sigma_n(t)を端点とする区間上のmmの積分以下であり、その区間はt0t_0とttを端点とする区間に含まれるから、hhの条件により∥x(t)−y0∥≤b\|\boldsymbol x(t)-\boldsymbol y_0\|\le bである。t∈I(k)t\in I^{(k)}では(Sk)(\mathrm S_k)の等式によりx(t)=xk(t)\boldsymbol x(t)=\boldsymbol x_k(t)である。したがってσn(I(k+1))⊆I(k)\sigma_n(I^{(k+1)})\subseteq I^{(k)}の上でf(s,x(s))=gk(s)f(s,\boldsymbol x(s))=\boldsymbol g_k(s)であり、x\boldsymbol xはI(k+1)I^{(k+1)}上で(Sk+1)(\mathrm S_{k+1})の等式を満たす。逆にx~\widetilde{\boldsymbol x}が(Sk+1)(\mathrm S_{k+1})の等式を満たす連続写像ならば、σn(I(k))⊆I(k)\sigma_n(I^{(k)})\subseteq I^{(k)}であるからx~∣I(k)\widetilde{\boldsymbol x}|_{I^{(k)}}は(Sk)(\mathrm S_k)の等式を満たし、x~∣I(k)=xk\widetilde{\boldsymbol x}|_{I^{(k)}}=\boldsymbol x_kである。このとき(Sk+1)(\mathrm S_{k+1})の等式の右辺はx(t)\boldsymbol x(t)の定義式に一致し、x~=x\widetilde{\boldsymbol x}=\boldsymbol xである。これで(Sk+1)(\mathrm S_{k+1})を得る。(Sn)(\mathrm S_n)が主張である。

主張 (2)を示す。s≤ts\le tとする。σn\sigma_nは単調非減少で、σn(t)−σn(s)≤t−s\sigma_n(t)-\sigma_n(s)\le t-sを満たす。§E13.4 系 3.1により

∥yn(t)−yn(s)∥=∥∫σn(s)σn(t)f(τ,yn(τ)) dτ∥≤∫σn(s)σn(t)m(τ) dτ≤ω(t−s)\|\boldsymbol y_n(t)-\boldsymbol y_n(s)\| =\left\|\int_{\sigma_n(s)}^{\sigma_n(t)}f(\tau,\boldsymbol y_n(\tau))\,d\tau\right\| \le\int_{\sigma_n(s)}^{\sigma_n(t)}m(\tau)\,d\tau \le\omega(t-s)

である。

主張 (3)を示す。yn(t)−y0−∫t0tf(s,yn(s)) ds=−∫σn(t)tf(s,yn(s)) ds\boldsymbol y_n(t)-\boldsymbol y_0-\int_{t_0}^tf(s,\boldsymbol y_n(s))\,ds=-\int_{\sigma_n(t)}^tf(s,\boldsymbol y_n(s))\,dsであり、σn(t)\sigma_n(t)とttを端点とする区間の長さはδn\delta_n以下である。§E13.4 系 3.1により、左辺のノルムはω(δn)\omega(\delta_n)以下である。

主張 (4)を示す。r≤r′r\le r'ならばω(r)\omega(r)の上限を取る集合はω(r′)\omega(r')の上限を取る集合に含まれるから、ω\omegaは単調非減少である。ε>0\ep>0とする。§E13.2 補題 1.3 (1)をmmに適用してδ>0\delta>0を取る。0≤r<δ0\le r<\deltaならば、s≤t≤s+rs\le t\le s+rを満たすs,t∈Js,t\in Jについてλ1([s,t])<δ\lambda_1([s,t])<\deltaであるから∫stm<ε\int_s^tm<\epであり、ω(r)≤ε\omega(r)\le\epである。▨

定理 1.2 (Carathéodory の存在定理).d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\|\cdot\|を固定する。t0∈Rt_0\in\R、y0∈Rd\boldsymbol y_0\in\R^d、a,b>0a,b>0とし、J=[t0−a,t0+a]J=[t_0-a,t_0+a]、B={u∈Rd:∥u−y0∥≤b}B=\{\boldsymbol u\in\R^d:\|\boldsymbol u-\boldsymbol y_0\|\le b\}とおく。写像f:J×B→Rdf:J\times B\to\R^dは、Lebesgue 零集合N⊆JN\subseteq Jと非負関数m∈L1(J)m\in L^1(J)について Carathéodory 条件を満たすとする。このとき、

∫t0−ht0m(s) ds≤b,∫t0t0+hm(s) ds≤b\int_{t_0-h}^{t_0}m(s)\,ds\le b,\qquad \int_{t_0}^{t_0+h}m(s)\,ds\le b

を満たす0<h≤a0<h\le aが存在する。そのようなhhを一つ取り、Ih=[t0−h,t0+h]I_h=[t_0-h,t_0+h]とおくと、IhI_h上の絶対連続関数y:Ih→Rd\boldsymbol y:I_h\to\R^dであって、すべてのt∈Iht\in I_hでy(t)∈B\boldsymbol y(t)\in Bかつ

y(t)=y0+∫t0tf(s,y(s)) ds\boldsymbol y(t)=\boldsymbol y_0+\int_{t_0}^tf(s,\boldsymbol y(s))\,ds

を満たすものが存在する。このy\boldsymbol yは(t0−h,t0+h)(t_0-h,t_0+h)のほとんど至る所でy′(t)=f(t,y(t))\boldsymbol y'(t)=f(t,\boldsymbol y(t))を満たす。

さらに、Ω⊆R×Rd\Omega\subseteq\R\times\R^dがJ×BJ\times Bを含む開集合であり、ffが写像F:Ω→RdF:\Omega\to\R^dのJ×BJ\times Bへの制限であるならば、y\boldsymbol yの(t0−h,t0+h)(t_0-h,t_0+h)への制限は、初期値問題y′=F(t,y)\boldsymbol y'=F(t,\boldsymbol y)、y(t0)=y0\boldsymbol y(t_0)=\boldsymbol y_0の絶対連続な局所解である。

証明.§E13.2 補題 1.3 (1)をJJ上の非負可積分関数mmに適用すると、あるδ>0\delta>0が存在して、JJに含まれる長さδ\delta未満の区間上のmmの積分はbb未満である。h=min⁡{a,δ/2}h=\min\{a,\delta/2\}は二つの不等式を満たす。

二つの不等式を満たす0<h≤a0<h\le aを一つ取る。補題 1.1のyn\boldsymbol y_n、ω\omega、δn\delta_nを取り、F={yn:n∈N≥1}\mathcal F=\{\boldsymbol y_n:n\in\NN\}とおく。Rd\R^dに距離(u,v)↦∥u−v∥(\boldsymbol u,\boldsymbol v)\mapsto\|\boldsymbol u-\boldsymbol v\|を入れた距離空間をYYと書く。§E10.1 系 2.3によりYYは完備であり、F⊆C(Ih,Y)\mathcal F\subseteq C(I_h,Y)である。

ε>0\ep>0に対して、補題 1.1 (4)によりω(r)<ε\omega(r)<\epとなるr>0r>0を取る。補題 1.1 (2)とω\omegaの単調性により、∣t−s∣<r\lvert t-s\rvert<rを満たすすべてのs,t∈Ihs,t\in I_hとすべてのnnについて∥yn(t)−yn(s)∥<ε\|\boldsymbol y_n(t)-\boldsymbol y_n(s)\|<\epである。したがってF\mathcal Fは同程度連続である。

∥⋅∥2\|\cdot\|_2を Euclid ノルムとする。§E10.1 命題 2.2により、定数c,C>0c,C>0が存在し、すべてのu∈Rd\boldsymbol u\in\R^dでc∥u∥2≤∥u∥≤C∥u∥2c\|\boldsymbol u\|_2\le\|\boldsymbol u\|\le C\|\boldsymbol u\|_2である。したがってBBは Euclid 距離について有界であり、u↦∥u−y0∥\boldsymbol u\mapsto\|\boldsymbol u-\boldsymbol y_0\|は Euclid 距離について連続であるからBBは閉集合である。§E2.9 定理 4.3によりBBは Euclid 距離についてコンパクトであり、§E2.9 系 3.3により全有界である。恒等写像(Rd,∥⋅∥2)→Y(\R^d,\|\cdot\|_2)\to YはCC-Lipschitz であるから、§E2.31 補題 2.2によりBBはYYの部分集合として全有界である。各t∈Iht\in I_hについて{yn(t):n∈N≥1}⊆B\{\boldsymbol y_n(t):n\in\NN\}\subseteq Bであるから、§E2.8 命題 3.1によりF\mathcal Fは各点で全有界である。

§E2.31 系 4.5により、部分列(ynk)k(\boldsymbol y_{n_k})_kとy∈C(Ih,Y)\boldsymbol y\in C(I_h,Y)が存在し、ynk\boldsymbol y_{n_k}はy\boldsymbol yへ一様収束する。BBはYYの閉集合であるから、すべてのt∈Iht\in I_hでy(t)∈B\boldsymbol y(t)\in Bである。

s∈Ihs\in I_hに対してgk(s)=f(s,ynk(s))\boldsymbol g_k(s)=f(s,\boldsymbol y_{n_k}(s))、g(s)=f(s,y(s))\boldsymbol g(s)=f(s,\boldsymbol y(s))とおく。y\boldsymbol yとynk\boldsymbol y_{n_k}をJJ上へ端点の値で定数に延長して§E13.4 補題 2.3を適用すると、gk\boldsymbol g_kとg\boldsymbol gはIhI_h上で可積分である。s∈Ih∖Ns\in I_h\setminus Nでは、§E13.4 定義 2.1 条件 (b)によりgk(s)→g(s)\boldsymbol g_k(s)\to\boldsymbol g(s)であり、∥gk(s)−g(s)∥≤2m(s)\|\boldsymbol g_k(s)-\boldsymbol g(s)\|\le2m(s)である。§E10.11 定理 3.1を可測関数列∥gk−g∥\|\boldsymbol g_k-\boldsymbol g\|に適用すると、∫Ih∥gk(s)−g(s)∥ ds→0\int_{I_h}\|\boldsymbol g_k(s)-\boldsymbol g(s)\|\,ds\to0である。t∈Iht\in I_hに対して、§E13.4 系 3.1と補題 1.1 (3)により

∥y(t)−y0−∫t0tg(s) ds∥≤∥y(t)−ynk(t)∥+∥ynk(t)−y0−∫t0tgk(s) ds∥+∥∫t0t(gk(s)−g(s)) ds∥≤max⁡τ∈Ih∥y(τ)−ynk(τ)∥+ω(δnk)+∫Ih∥gk(s)−g(s)∥ ds\begin{aligned} \left\|\boldsymbol y(t)-\boldsymbol y_0-\int_{t_0}^t\boldsymbol g(s)\,ds\right\| &\le\|\boldsymbol y(t)-\boldsymbol y_{n_k}(t)\| +\left\|\boldsymbol y_{n_k}(t)-\boldsymbol y_0-\int_{t_0}^t\boldsymbol g_k(s)\,ds\right\| +\left\|\int_{t_0}^t(\boldsymbol g_k(s)-\boldsymbol g(s))\,ds\right\|\\ &\le\max_{\tau\in I_h}\|\boldsymbol y(\tau)-\boldsymbol y_{n_k}(\tau)\|+\omega(\delta_{n_k})+\int_{I_h}\|\boldsymbol g_k(s)-\boldsymbol g(s)\|\,ds \end{aligned}

である。k→∞k\to\inftyのとき、右辺の第一項は一様収束により、第二項はδnk→0\delta_{n_k}\to0と補題 1.1 (4)により、第三項は上の優収束定理の適用により、いずれも00へ収束する。よってy\boldsymbol yは積分方程式を満たす。

g∈L1(Ih;Rd)\boldsymbol g\in L^1(I_h;\R^d)であるから、積分の加法性によりy(t)=y(t0−h)+∫t0−htg(s) ds\boldsymbol y(t)=\boldsymbol y(t_0-h)+\int_{t_0-h}^t\boldsymbol g(s)\,dsであり、§E13.2 補題 1.3 (3)によりy\boldsymbol yはIhI_h上で絶対連続である。§E13.2 定理 2.5により、(t0−h,t0+h)(t_0-h,t_0+h)のほとんど至る所でy′=g\boldsymbol y'=\boldsymbol gである。

ffがFFの制限である場合、y\boldsymbol yの(t0−h,t0+h)(t_0-h,t_0+h)への制限のグラフはJ×B⊆ΩJ\times B\subseteq\Omegaに含まれ、s↦F(s,y(s))=g(s)s\mapsto F(s,\boldsymbol y(s))=\boldsymbol g(s)は可積分である。§E13.4 補題 2.4により、この制限はt0t_0を含む開区間上の絶対連続解である。▨

系 1.3 (Peano の存在定理).d∈N≥1d\in\NNとし、Rd\R^dにノルム∥⋅∥\|\cdot\|を固定する。Ω⊆R×Rd\Omega\subseteq\R\times\R^dを開集合、f:Ω→Rdf:\Omega\to\R^dを連続写像、(t0,y0)∈Ω(t_0,\boldsymbol y_0)\in\Omegaとする。a,b>0a,b>0について

R=[t0−a,t0+a]×{u∈Rd:∥u−y0∥≤b}⊆ΩR=[t_0-a,t_0+a]\times\{\boldsymbol u\in\R^d:\|\boldsymbol u-\boldsymbol y_0\|\le b\}\subseteq\Omega

であるとし、M=max⁡(t,u)∈R∥f(t,u)∥M=\max_{(t,\boldsymbol u)\in R}\|f(t,\boldsymbol u)\|とおく(RRはコンパクトでありffは連続であるから、最大値が存在する)。M=0M=0ならばh=ah=a、M>0M>0ならばh=min⁡{a,b/M}h=\min\{a,b/M\}とする。このとき、初期値問題y′=f(t,y)\boldsymbol y'=f(t,\boldsymbol y)、y(t0)=y0\boldsymbol y(t_0)=\boldsymbol y_0は、[t0−h,t0+h][t_0-h,t_0+h]上の古典解であって、すべてのt∈[t0−h,t0+h]t\in[t_0-h,t_0+h]で(t,y(t))∈R(t,\boldsymbol y(t))\in Rを満たすものをもつ。

証明.f∣Rf|_Rは、N=∅N=\varnothingと定数関数m≡Mm\equiv Mについて Carathéodory 条件を満たす。実際、各u\boldsymbol uについてt↦f(t,u)t\mapsto f(t,\boldsymbol u)は連続であるから Lebesgue 可測であり、各ttについてu↦f(t,u)\boldsymbol u\mapsto f(t,\boldsymbol u)は連続であり、RR上で∥f∥≤M\|f\|\le Mである。hhの選び方から0<h≤a0<h\le aとMh≤bMh\le bが成り立つので、t0t_0の左右のmmの積分はともにbb以下である。定理 1.2により、[t0−h,t0+h][t_0-h,t_0+h]上の写像y\boldsymbol yであって、値が閉球に入り、全点で積分方程式を満たすものが存在する。y\boldsymbol yのグラフはR⊆ΩR\subseteq\Omegaに含まれるから、§E13.4 補題 2.5を区間[t0−h,t0+h][t_0-h,t_0+h]に適用すると、y\boldsymbol yは古典解である。▨

注意 1.4.定理 1.2の証明でRd\R^dが有限次元であることを用いるのは、閉球BBが全有界であることを確かめる段だけである。一方、無限次元の Banach 空間では、閉球は全有界でなく、Peano の存在定理そのものが成り立たない。実際、Godunov により、任意の無限次元の Banach 空間EEについて、連続写像f:R×E→Ef:\R\times E\to Eと(t0,y0)∈R×E(t_0,\boldsymbol y_0)\in\R\times Eであって、初期値問題y′=f(t,y)\boldsymbol y'=f(t,\boldsymbol y)、y(t0)=y0\boldsymbol y(t_0)=\boldsymbol y_0がt0t_0を含むどの開区間の上にも一階連続微分可能な解をもたないものが存在する。

注意 1.5. 格子点tk=t0+kh/nt_k=t_0+kh/n(k=0,1,…,nk=0,1,\ldots,n)を取り、y(t0)=y0\boldsymbol y(t_0)=\boldsymbol y_0から始めて[tk,tk+1][t_k,t_{k+1}]上でy(t)=y(tk)+(t−tk)f(tk,y(tk))\boldsymbol y(t)=\boldsymbol y(t_k)+(t-t_k)f(t_k,\boldsymbol y(t_k))と順に定める Euler 折れ線の定義に現れるのは、有限個の時刻tkt_kでのffの値だけである。一方、零集合N′⊆JN'\subseteq Jを取り、ffをN′×BN'\times Bの上だけで変更しても、Carathéodory 条件は零集合N∪N′N\cup N'について保たれる。連続写像x\boldsymbol xに対するs↦f(s,x(s))s\mapsto f(s,\boldsymbol x(s))はN′N'の外で変わらないので、その積分は変わらない。したがって補題 1.1 (1)の一意性により Tonelli の近似解yn\boldsymbol y_nは変わらない。値がBBに入る絶対連続関数y\boldsymbol yについても、t∉N′t\notin N'ではf(t,y(t))f(t,\boldsymbol y(t))が変わらないから、ほとんど至る所での等式y′(t)=f(t,y(t))\boldsymbol y'(t)=f(t,\boldsymbol y(t))は変更の前後で同値であり、値がBBに入る絶対連続解の全体も変わらない。

一方、Euler 折れ線は、ffをN′×BN'\times Bの上だけで変更すると変わることがある。d=1d=1とし、R2\R^2上のf(t,u)f(t,u)をttが有理数ならば11、無理数ならば00と定める。ffは任意の閉矩形の上でN=∅N=\varnothing、m≡1m\equiv1について Carathéodory 条件を満たす。ffは零集合Q×R\Q\times\Rの外で00に等しいから、t0=0t_0=0を含む区間上でy(0)=0y(0)=0を満たす絶対連続解yyはほとんど至る所でy′=0y'=0を満たし、§E13.2 系 2.6によりy≡0y\equiv0である。一方、t0=0t_0=0と有理数hhについて、すべてのtkt_kは有理数であり、[0,h][0,h]上の Euler 折れ線はy(t)=ty(t)=tである。したがって、y(t)=ty(t)=tは解でない。

注意 1.6. 他の文献では、支配の条件を、各u∈B\boldsymbol u\in Bごとに零集合NuN_{\boldsymbol u}を許して「t∉Nut\notin N_{\boldsymbol u}ならば∥f(t,u)∥≤m(t)\|f(t,\boldsymbol u)\|\le m(t)」と述べることもある。§E13.4 定義 2.1のJJ、BBと写像f:J×B→Rdf:J\times B\to\R^dが§E13.4 定義 2.1 条件 (a)と、零集合NNについての§E13.4 定義 2.1 条件 (b)を満たすとき、二つの形は同値である。実際、共通の零集合NNの外で支配が成り立つならば、すべてのu∈B\boldsymbol u\in BについてNu=NN_{\boldsymbol u}=Nとして状態ごとの零集合による形が成り立つ。逆に、BBの可算稠密集合QQを取るとN′=N∪⋃q∈QNqN'=N\cup\bigcup_{\boldsymbol q\in Q}N_{\boldsymbol q}は零集合であり、t∉N′t\notin N'ならばすべてのq∈Q\boldsymbol q\in Qで∥f(t,q)∥≤m(t)\|f(t,\boldsymbol q)\|\le m(t)であるから、f(t,⋅)f(t,\cdot)の連続性によりすべてのu∈B\boldsymbol u\in Bで∥f(t,u)∥≤m(t)\|f(t,\boldsymbol u)\|\le m(t)である。

2 仮定を外した場合の解の非存在

例 2.1.R2\R^2上の写像

f(t,u)={−1,u≥0,1,u<0f(t,u)= \begin{cases} -1,&u\ge0,\\ 1,&u<0 \end{cases}

と初期条件y(0)=0y(0)=0を考える。このとき、a,b>0a,b>0について、ffは[−a,a]×[−b,b][-a,a]\times[-b,b]上で§E13.4 定義 2.1 条件 (a)と、m≡1m\equiv1についての§E13.4 定義 2.1 条件 (c)を満たす。すべてのttでu↦f(t,u)u\mapsto f(t,u)はu=0u=0で連続でないから、ffは§E13.4 定義 2.1 条件 (b)を満たさない。任意のε>0\ep>0について、y∈AC([0,ε];R)y\in AC([0,\ep];\R)であってy(0)=0y(0)=0かつ(0,ε)(0,\ep)のほとんど至る所でy′(t)=f(t,y(t))y'(t)=f(t,y(t))を満たすものは存在しない。特に、原点を含む開区間上の絶対連続な局所解は存在しない。

実際、そのようなyyがあるとする。{y>0}\{y>0\}のほとんど至る所でy′=−1y'=-1であるから、§E13.2 補題 3.3 (1)により(0,ε)(0,\ep)のほとんど至る所で(y+)′=−1{y>0}≤0(y^+)'=-\mathbf 1_{\{y>0\}}\le0であり、§E13.2 定理 2.5によりy+(t)=y+(0)+∫0t(y+)′(s) ds≤0y^+(t)=y^+(0)+\int_0^t(y^+)'(s)\,ds\le0である。よって[0,ε][0,\ep]上でy+≡0y^+\equiv0である。{y<0}\{y<0\}のほとんど至る所でy′=1y'=1であるから、§E13.2 補題 3.3 (1)を−y-yに適用すると、ほとんど至る所で((−y)+)′=−1{y<0}y′=−1{y<0}≤0((-y)^+)'=-\mathbf 1_{\{y<0\}}y'=-\mathbf 1_{\{y<0\}}\le0であり、同じ積分表示から(−y)+≡0(-y)^+\equiv0である。したがって[0,ε][0,\ep]上でy≡0y\equiv0であり、すべてのt∈(0,ε)t\in(0,\ep)でy′(t)=0y'(t)=0である。一方、f(t,0)=−1f(t,0)=-1であるから、仮定により(0,ε)(0,\ep)のほとんど至る所でy′(t)=−1y'(t)=-1である。正の測度をもつ(0,ε)(0,\ep)の上で、すべての点でy′=0y'=0であることとほとんど至る所でy′=−1y'=-1であることは両立しない。

左側では、y(t)=−ty(t)=-t(t≤0t\le0)がy(0)=0y(0)=0とすべてのt<0t<0でのy′(t)=f(t,y(t))y'(t)=f(t,y(t))を満たす。

例 2.2.R2\R^2上の写像

f(t,u)={1/t,t>0,0,t≤0f(t,u)= \begin{cases} 1/t,&t>0,\\ 0,&t\le0 \end{cases}

を考える。このとき、a,b>0a,b>0について、ffは[−a,a]×[−b,b][-a,a]\times[-b,b]上でN=∅N=\varnothingについて§E13.4 定義 2.1 条件 (a)と§E13.4 定義 2.1 条件 (b)を満たす。非負関数mmと零集合N′N'がt∉N′t\notin N'で∣f(t,u)∣≤m(t)\lvert f(t,u)\rvert\le m(t)を満たすならば、0<δ<a0<\delta<aについて∫δam(t) dt≥∫δadt/t=log⁡(a/δ)\int_\delta^am(t)\,dt\ge\int_\delta^a dt/t=\log(a/\delta)であり、δ↓0\delta\downarrow0で右辺は非有界であるから、m∉L1([−a,a])m\notin L^1([-a,a])である。したがってffは§E13.4 定義 2.1 条件 (c)を満たさない。

任意のy0∈Ry_0\in\Rとε>0\ep>0について、y∈AC([0,ε];R)y\in AC([0,\ep];\R)であってy(0)=y0y(0)=y_0かつ(0,ε)(0,\ep)のほとんど至る所でy′(t)=f(t,y(t))y'(t)=f(t,y(t))を満たすものは存在しない。実際、そのようなyyがあるとすると、0<δ<ε0<\delta<\epについて§E13.2 定理 2.5を[δ,ε][\delta,\ep]に適用して

y(ε)−y(δ)=∫δεdtt=log⁡εδy(\ep)-y(\delta)=\int_\delta^\ep\frac{dt}{t}=\log\frac\ep\delta

を得る。yyは[0,ε][0,\ep]上で連続であるから、左辺の絶対値はδ\deltaによらず2max⁡[0,ε]∣y∣2\max_{[0,\ep]}\lvert y\rvert以下である。一方、右辺はδ↓0\delta\downarrow0で非有界である。特に、どの初期値についても、原点を含む開区間上の絶対連続な局所解は存在しない。左側では、定数関数y≡y0y\equiv y_0が(−∞,0](-\infty,0]上で方程式を満たす。

3 スカラー方程式の最大の解と最小の解

定義 3.1.t0,y0∈Rt_0,y_0\in\R、a,b>0a,b>0とし、J=[t0−a,t0+a]J=[t_0-a,t_0+a]、B=[y0−b,y0+b]B=[y_0-b,y_0+b]、f:J×B→Rf:J\times B\to\Rを写像とする。

  1. 0<h′≤a0<h'\le aとする。絶対連続関数y:[t0,t0+h′]→By:[t_0,t_0+h']\to Bがy(t0)=y0y(t_0)=y_0を満たし、(t0,t0+h′)(t_0,t_0+h')のほとんど至る所でy′(t)=f(t,y(t))y'(t)=f(t,y(t))を満たすとき、yyを[t0,t0+h′][t_0,t_0+h']上の値がBBに入る絶対連続解という。
  2. 0<h≤a0<h\le aとし、y∗,y∗y^*,y_*を[t0,t0+h][t_0,t_0+h]上の値がBBに入る絶対連続解とする。任意の0<h′≤h0<h'\le hと、[t0,t0+h′][t_0,t_0+h']上の値がBBに入る任意の絶対連続解yyについて、[t0,t0+h′][t_0,t_0+h']上でy≤y∗y\le y^*が成り立つとき、y∗y^*を[t0,t0+h][t_0,t_0+h]上の最大の解 (maximum solution) という。同じh′h'とyyのすべてについて[t0,t0+h′][t_0,t_0+h']上でy≥y∗y\ge y_*が成り立つとき、y∗y_*を[t0,t0+h][t_0,t_0+h]上の最小の解 (minimum solution) という。
  3. α<β\alpha<\betaを[α,β]⊆J[\alpha,\beta]\subseteq Jを満たす実数とする。絶対連続関数v:[α,β]→Bv:[\alpha,\beta]\to Bが(α,β)(\alpha,\beta)のほとんど至る所でv′(t)≤f(t,v(t))v'(t)\le f(t,v(t))を満たすとき、vvを一階の方程式y′=f(t,y)y'=f(t,y)の[α,β][\alpha,\beta]上の下解 (lower solution) という。(α,β)(\alpha,\beta)のほとんど至る所でv′(t)≥f(t,v(t))v'(t)\ge f(t,v(t))を満たすとき、vvを[α,β][\alpha,\beta]上の上解 (upper solution) という。

最大の解と最小の解は値の順序についての最大と最小であり、延長することができない解を指す最大解とは異なる概念である。

補題 3.2.J,B⊆RJ,B\subseteq\Rを区間、f:J×B→Rf:J\times B\to\Rを写像とし、α<β\alpha<\betaを[α,β]⊆J[\alpha,\beta]\subseteq Jを満たす実数、y1,y2:[α,β]→By_1,y_2:[\alpha,\beta]\to Bを絶対連続関数とする。このとき、(α,β)(\alpha,\beta)のほとんど至る所でy1′(t)≤f(t,y1(t))y_1'(t)\le f(t,y_1(t))とy2′(t)≤f(t,y2(t))y_2'(t)\le f(t,y_2(t))が成り立つならば、max⁡{y1,y2}\max\{y_1,y_2\}は[α,β][\alpha,\beta]からBBへの絶対連続関数であり、(α,β)(\alpha,\beta)のほとんど至る所でmax⁡{y1,y2}′(t)≤f(t,max⁡{y1,y2}(t))\max\{y_1,y_2\}'(t)\le f(t,\max\{y_1,y_2\}(t))を満たす。(α,β)(\alpha,\beta)のほとんど至る所でy1′(t)≥f(t,y1(t))y_1'(t)\ge f(t,y_1(t))とy2′(t)≥f(t,y2(t))y_2'(t)\ge f(t,y_2(t))が成り立つならば、min⁡{y1,y2}\min\{y_1,y_2\}は[α,β][\alpha,\beta]からBBへの絶対連続関数であり、(α,β)(\alpha,\beta)のほとんど至る所でmin⁡{y1,y2}′(t)≥f(t,min⁡{y1,y2}(t))\min\{y_1,y_2\}'(t)\ge f(t,\min\{y_1,y_2\}(t))を満たす。いずれの主張も、仮定と結論の不等号をすべて等号に替えて成り立つ。

証明.BBは区間であるから、w=max⁡{y1,y2}w=\max\{y_1,y_2\}とmin⁡{y1,y2}\min\{y_1,y_2\}の値はBBに入る。§E13.2 補題 3.3 (3)によりwwは絶対連続であり、ほとんど至る所でw′=1{y1>y2}y1′+1{y1≤y2}y2′w'=\mathbf 1_{\{y_1>y_2\}}y_1'+\mathbf 1_{\{y_1\le y_2\}}y_2'である。この等式と二つの仮定がすべて成り立つ点t∈(α,β)t\in(\alpha,\beta)を取る。そのような点の全体の補集合は零集合である。y1(t)>y2(t)y_1(t)>y_2(t)ならば、w(t)=y1(t)w(t)=y_1(t)でありw′(t)=y1′(t)≤f(t,y1(t))=f(t,w(t))w'(t)=y_1'(t)\le f(t,y_1(t))=f(t,w(t))である。y1(t)≤y2(t)y_1(t)\le y_2(t)ならば、w(t)=y2(t)w(t)=y_2(t)でありw′(t)=y2′(t)≤f(t,y2(t))=f(t,w(t))w'(t)=y_2'(t)\le f(t,y_2(t))=f(t,w(t))である。仮定が等号ならば、二つの場合の不等号は等号である。

min⁡{y1,y2}=−max⁡{−y1,−y2}\min\{y_1,y_2\}=-\max\{-y_1,-y_2\}に§E13.2 補題 3.3 (3)を適用すると、min⁡{y1,y2}\min\{y_1,y_2\}は絶対連続であり、ほとんど至る所でmin⁡{y1,y2}′=1{y1<y2}y1′+1{y1≥y2}y2′\min\{y_1,y_2\}'=\mathbf 1_{\{y_1<y_2\}}y_1'+\mathbf 1_{\{y_1\ge y_2\}}y_2'である。y1(t)<y2(t)y_1(t)<y_2(t)の点ではmin⁡{y1,y2}′(t)=y1′(t)≥f(t,y1(t))=f(t,min⁡{y1,y2}(t))\min\{y_1,y_2\}'(t)=y_1'(t)\ge f(t,y_1(t))=f(t,\min\{y_1,y_2\}(t))であり、y1(t)≥y2(t)y_1(t)\ge y_2(t)の点ではmin⁡{y1,y2}′(t)=y2′(t)≥f(t,y2(t))=f(t,min⁡{y1,y2}(t))\min\{y_1,y_2\}'(t)=y_2'(t)\ge f(t,y_2(t))=f(t,\min\{y_1,y_2\}(t))である。仮定が等号ならば、二つの場合の不等号は等号である。▨

補題 3.3.t0,y0∈Rt_0,y_0\in\R、a,b>0a,b>0とし、J=[t0−a,t0+a]J=[t_0-a,t_0+a]、B=[y0−b,y0+b]B=[y_0-b,y_0+b]、−B=[−y0−b,−y0+b]-B=[-y_0-b,-y_0+b]とおく。写像f:J×B→Rf:J\times B\to\Rに対して、f−:J×(−B)→Rf^-:J\times(-B)\to\Rをf−(t,u)=−f(t,−u)f^-(t,u)=-f(t,-u)で定め、f−f^-には初期値−y0-y_0と区間−B-Bについて定義 3.1を用いる。このとき、次が成り立つ。

  1. ffが Lebesgue 零集合N⊆JN\subseteq Jと非負関数m∈L1(J)m\in L^1(J)について Carathéodory 条件を満たすならば、f−f^-もNNとmmについて Carathéodory 条件を満たす。
  2. 0<h′≤a0<h'\le aとする。yyがffの[t0,t0+h′][t_0,t_0+h']上の値がBBに入る絶対連続解であることと、−y-yがf−f^-の[t0,t0+h′][t_0,t_0+h']上の値が−B-Bに入る絶対連続解であることは同値である。
  3. α<β\alpha<\betaを[α,β]⊆J[\alpha,\beta]\subseteq Jを満たす実数とする。vvがffの[α,β][\alpha,\beta]上の上解であることと、−v-vがf−f^-の[α,β][\alpha,\beta]上の下解であることは同値である。
  4. 0<h≤a0<h\le aとする。y^\hat yがf−f^-の[t0,t0+h][t_0,t_0+h]上の最大の解であることと、−y^-\hat yがffの[t0,t0+h][t_0,t_0+h]上の最小の解であることは同値である。

証明.主張 (1)を示す。各u∈−Bu\in-Bについてt↦f−(t,u)=−f(t,−u)t\mapsto f^-(t,u)=-f(t,-u)は可測である。t∈J∖Nt\in J\setminus Nならばu↦−f(t,−u)u\mapsto-f(t,-u)は−B-B上で連続であり、∣f−(t,u)∣=∣f(t,−u)∣≤m(t)\lvert f^-(t,u)\rvert=\lvert f(t,-u)\rvert\le m(t)である。

主張 (2)と主張 (3)を示す。§E13.2 補題 1.2 (2)により、yyが絶対連続であることと−y-yが絶対連続であることは同値である。yyの値がBBに入ることと−y-yの値が−B-Bに入ることは同値であり、y(t0)=y0y(t_0)=y_0と−y(t0)=−y0-y(t_0)=-y_0は同値である。yyが微分可能な点ttでは(−y)′(t)=−y′(t)(-y)'(t)=-y'(t)であり、y′(t)=f(t,y(t))y'(t)=f(t,y(t))と(−y)′(t)=f−(t,−y(t))(-y)'(t)=f^-(t,-y(t))は同値であり、y′(t)≥f(t,y(t))y'(t)\ge f(t,y(t))と(−y)′(t)≤f−(t,−y(t))(-y)'(t)\le f^-(t,-y(t))は同値である。

主張 (4)を示す。主張 (2)により、y↦−yy\mapsto-yは、ffの[t0,t0+h′][t_0,t_0+h']上の値がBBに入る絶対連続解の全体から、f−f^-の[t0,t0+h′][t_0,t_0+h']上の値が−B-Bに入る絶対連続解の全体への全単射である。−y≤y^-y\le\hat yとy≥−y^y\ge-\hat yは同値であるから、定義 3.1 (2)の最大の解の条件をf−f^-とy^\hat yについて述べたものは、最小の解の条件をffと−y^-\hat yについて述べたものと同値である。▨

定理 3.4.t0,y0∈Rt_0,y_0\in\R、a,b>0a,b>0とし、J=[t0−a,t0+a]J=[t_0-a,t_0+a]、B=[y0−b,y0+b]B=[y_0-b,y_0+b]とおく。写像f:J×B→Rf:J\times B\to\Rは、Lebesgue 零集合N⊆JN\subseteq Jと非負関数m∈L1(J)m\in L^1(J)について Carathéodory 条件を満たすとし、0<h≤a0<h\le aは∫t0−ht0m(s) ds≤b\int_{t_0-h}^{t_0}m(s)\,ds\le bと∫t0t0+hm(s) ds≤b\int_{t_0}^{t_0+h}m(s)\,ds\le bを満たすとする。このとき、[t0,t0+h][t_0,t_0+h]上の最大の解y∗y^*と最小の解y∗y_*がそれぞれただ一つ存在する。

さらに、ffが時間について可積分な Lipschitz 係数ℓ∈L1(J)\ell\in L^1(J)をもつならば、任意の0<h′≤h0<h'\le hについて、[t0,t0+h′][t_0,t_0+h']上の値がBBに入る絶対連続解はy∗y^*の[t0,t0+h′][t_0,t_0+h']への制限ただ一つである。特にy∗=y∗y^*=y_*である。

証明.[t0,t0+h][t_0,t_0+h]上の値がBBに入る絶対連続解の全体をSSとする。y1∗,y2∗y_1^*,y_2^*がともに最大の解ならば、定義 3.1 (2)をh′=hh'=hで適用してy1∗≤y2∗≤y1∗y_1^*\le y_2^*\le y_1^*を得る。y∗1,y∗2y_{*1},y_{*2}がともに最小の解ならば、同じ適用によりy∗1≥y∗2≥y∗1y_{*1}\ge y_{*2}\ge y_{*1}である。

0<h′≤h0<h'\le hとし、yyを[t0,t0+h′][t_0,t_0+h']上の値がBBに入る絶対連続解とする。§E13.2 定理 2.5とほとんど至る所での∣y′(r)∣=∣f(r,y(r))∣≤m(r)\lvert y'(r)\rvert=\lvert f(r,y(r))\rvert\le m(r)により、t0≤s≤t≤t0+h′t_0\le s\le t\le t_0+h'について

y(t)=y0+∫t0tf(r,y(r)) dry(t)=y_0+\int_{t_0}^tf(r,y(r))\,dr(3.4.1)\text{(3.4.1)}

と

∣y(t)−y(s)∣≤∫stm(r) dr\lvert y(t)-y(s)\rvert\le\int_s^tm(r)\,dr(3.4.2)\text{(3.4.2)}

が成り立つ。逆に、連続関数x:[t0,t0+h]→Bx:[t_0,t_0+h]\to Bがすべてのt∈[t0,t0+h]t\in[t_0,t_0+h]で式 (3.4.1)を満たすならば、s↦f(s,x(s))s\mapsto f(s,x(s))は§E13.4 補題 2.3により可積分であり、§E13.2 補題 1.3 (3)と§E13.2 定理 2.5によりx∈Sx\in Sである。

定理 1.2の解を[t0,t0+h][t_0,t_0+h]に制限するとSSの元を得るので、S≠∅S\ne\varnothingである。y1,y2∈Sy_1,y_2\in Sならばmax⁡{y1,y2}(t0)=y0\max\{y_1,y_2\}(t_0)=y_0であり、補題 3.2の等号の場合によりmax⁡{y1,y2}∈S\max\{y_1,y_2\}\in Sである。t∈[t0,t0+h]t\in[t_0,t_0+h]に対してy∗(t)=sup⁡y∈Sy(t)≤y0+by^*(t)=\sup_{y\in S}y(t)\le y_0+bとおく。t0≤s≤t≤t0+ht_0\le s\le t\le t_0+hとy∈Sy\in Sについて、式 (3.4.2)によりy(t)≤y∗(s)+∫stmy(t)\le y^*(s)+\int_s^tmとy(s)≤y∗(t)+∫stmy(s)\le y^*(t)+\int_s^tmであるから、yyについて上限を取ると∣y∗(t)−y∗(s)∣≤∫stm\lvert y^*(t)-y^*(s)\rvert\le\int_s^tmである。§E13.2 補題 1.3 (1)によりy∗y^*は連続である。

[t0,t0+h][t_0,t_0+h]の可算稠密集合{τj:j∈N≥1}\{\tau_j:j\in\NN\}を取る。j,k∈N≥1j,k\in\NNに対して、yj,k(τj)>y∗(τj)−1/ky_{j,k}(\tau_j)>y^*(\tau_j)-1/kを満たすyj,k∈Sy_{j,k}\in Sを取る。全単射N≥1→N≥1×N≥1\NN\to\NN\times\NNによって{yj,k}\{y_{j,k}\}をz1,z2,…z_1,z_2,\ldotsと並べ、wn=max⁡{z1,…,zn}w_n=\max\{z_1,\ldots,z_n\}とおく。wn∈Sw_n\in Sであり、wn≤wn+1≤y∗w_n\le w_{n+1}\le y^*である。各ttで単調有界な列(wn(t))(w_n(t))は極限w(t)≤y∗(t)w(t)\le y^*(t)をもつ。j,kj,kを固定すると、十分大きいnnについてwn(τj)≥yj,k(τj)>y∗(τj)−1/kw_n(\tau_j)\ge y_{j,k}(\tau_j)>y^*(\tau_j)-1/kであるから、w(τj)=y∗(τj)w(\tau_j)=y^*(\tau_j)である。wnw_nについての式 (3.4.2)でn→∞n\to\inftyとすると∣w(t)−w(s)∣≤∫stm\lvert w(t)-w(s)\rvert\le\int_s^tmであり、wwは連続である。連続関数wwとy∗y^*は稠密集合の上で一致するから、w=y∗w=y^*である。

s∈[t0,t0+h]∖Ns\in[t_0,t_0+h]\setminus Nでは、§E13.4 定義 2.1 条件 (b)によりf(s,wn(s))→f(s,y∗(s))f(s,w_n(s))\to f(s,y^*(s))であり、∣f(s,wn(s))−f(s,y∗(s))∣≤2m(s)\lvert f(s,w_n(s))-f(s,y^*(s))\rvert\le2m(s)である。§E10.11 定理 3.1により、wnw_nについての式 (3.4.1)でn→∞n\to\inftyとすると

y∗(t)=y0+∫t0tf(s,y∗(s)) ds(t∈[t0,t0+h])y^*(t)=y_0+\int_{t_0}^tf(s,y^*(s))\,ds\qquad(t\in[t_0,t_0+h])

を得る。BBは閉集合であるからy∗y^*の値はBBに入り、y∗∈Sy^*\in Sである。

0<h′≤h0<h'\le hとし、yyを[t0,t0+h′][t_0,t_0+h']上の値がBBに入る絶対連続解とする。h′=hh'=hならばy∈Sy\in Sであり、y≤y∗y\le y^*である。h′<hh'<hとし、t1=t0+h′t_1=t_0+h'、π(u)=min⁡{max⁡{u,y0−b},y0+b}\pi(u)=\min\{\max\{u,y_0-b\},y_0+b\}とおく。R2\R^2上の写像ggを、t∈Jt\in Jならばg(t,u)=f(t,π(u))g(t,u)=f(t,\pi(u))、t∉Jt\notin Jならばg(t,u)=0g(t,u)=0と定め、R\R上のm~\widetilde mをJJ上でmm、JJの外で00と定める。β=1+∫Jm\beta=1+\int_Jmとおく。ggは[t1−h,t1+h]×[y(t1)−β,y(t1)+β][t_1-h,t_1+h]\times[y(t_1)-\beta,y(t_1)+\beta]上で、零集合NNと支配関数m~\widetilde mについて Carathéodory 条件を満たし、t1t_1の左右のm~\widetilde mの積分はともにβ\beta以下である。定理 1.2により、[t1−h,t1+h][t_1-h,t_1+h]上の連続関数zzであって、すべてのttでz(t)=y(t1)+∫t1tg(s,z(s)) dsz(t)=y(t_1)+\int_{t_1}^tg(s,z(s))\,dsを満たすものが存在する。t∈[t1,t0+h]t\in[t_1,t_0+h]では[t1,t]⊆J[t_1,t]\subseteq Jであり、yyについての式 (3.4.2)と§E13.4 系 3.1により

∣z(t)−y0∣≤∣y(t1)−y0∣+∣z(t)−z(t1)∣≤∫t0t1m(s) ds+∫t1tm(s) ds≤b\lvert z(t)-y_0\rvert\le\lvert y(t_1)-y_0\rvert+\lvert z(t)-z(t_1)\rvert \le\int_{t_0}^{t_1}m(s)\,ds+\int_{t_1}^tm(s)\,ds\le b

である。したがって[t1,t0+h][t_1,t_0+h]上でz(t)∈Bz(t)\in B、g(t,z(t))=f(t,z(t))g(t,z(t))=f(t,z(t))である。[t0,t1][t_0,t_1]上でy^=y\hat y=y、[t1,t0+h][t_1,t_0+h]上でy^=z\hat y=zと定めると、y^\hat yは連続で値はBBに入る。t∈[t1,t0+h]t\in[t_1,t_0+h]について、yyについての式 (3.4.1)により

y^(t)=y(t1)+∫t1tf(s,z(s)) ds=y0+∫t0tf(s,y^(s)) ds\hat y(t)=y(t_1)+\int_{t_1}^tf(s,z(s))\,ds=y_0+\int_{t_0}^tf(s,\hat y(s))\,ds

であり、t∈[t0,t1]t\in[t_0,t_1]についても同じ等式が成り立つ。よってy^∈S\hat y\in Sであり、[t0,t1][t_0,t_1]上でy=y^≤y∗y=\hat y\le y^*である。以上より、y∗y^*は最大の解である。

補題 3.3のf−f^-を取る。補題 3.3 (1)によりf−f^-はNNとmmについて Carathéodory 条件を満たすので、前段をf−f^-に適用して[t0,t0+h][t_0,t_0+h]上の最大の解y^∗\hat y^*を取る。補題 3.3 (4)により、y∗=−y^∗y_*=-\hat y^*はffの最小の解である。

ffが可積分な Lipschitz 係数ℓ\ellをもつとする。R2\R^2上の写像FFを、J×BJ\times B上でff、J×BJ\times Bの外で00と定める。FFはJ×BJ\times B上でffに一致するから、開集合R2\R^2と閉矩形J×BJ\times B、支配関数mm、Lipschitz 係数ℓ\ellについて§E13.4 定理 3.5の仮定を満たす。JJの部分区間上の関数zzの値がBBに入るならばF(t,z(t))=f(t,z(t))F(t,z(t))=f(t,z(t))であるから、値がBBに入るFFの絶対連続解と値がBBに入るffの絶対連続解は一致する。§E13.4 定理 3.5 (1)により、Ih=[t0−h,t0+h]I_h=[t_0-h,t_0+h]上の絶対連続解yLy_Lで値がBBに入るものが存在する。0<h′≤h0<h'\le hとし、yyを[t0,t0+h′][t_0,t_0+h']上の値がBBに入る絶対連続解とする。I=[t0,t0+h′]⊆IhI=[t_0,t_0+h']\subseteq I_hはt0t_0を含み内部が空でないから、§E13.4 定理 3.5 (2)によりII上でy=yLy=y_Lである。y∗y^*の[t0,t0+h′][t_0,t_0+h']への制限も同じ条件を満たすからy∗=yLy^*=y_Lであり、[t0,t0+h′][t_0,t_0+h']上でy=y∗y=y^*である。h′=hh'=hとしてy∗y_*に適用するとy∗=y∗y_*=y^*である。▨

例 3.5.t0=0t_0=0、y0=0y_0=0、a=b=1a=b=1とし、J=B=[−1,1]J=B=[-1,1]上のf(t,u)=∣u∣f(t,u)=\sqrt{\lvert u\rvert}を考える。このとき、ffは連続であり、N=∅N=\varnothing、m≡1m\equiv1について Carathéodory 条件を満たし、h=1h=1は∫−10m=∫01m=1≤b\int_{-1}^0m=\int_0^1m=1\le bを満たす。0<h′≤10<h'\le1とすると、[0,h′][0,h']上の値がBBに入る絶対連続解は、§E13.4 例 5.3のycy_c(c≥0c\ge0)の[0,h′][0,h']への制限に限る。したがって[0,1][0,1]上の最大の解はy∗(t)=t2/4y^*(t)=t^2/4、最小の解はy∗≡0y_*\equiv0であり、すべてのc≥0c\ge0についてy∗≤yc≤y∗y_*\le y_c\le y^*が[0,1][0,1]上で成り立つ。系 1.3の解は一意とは限らない。

実際、yyを[0,h′][0,h']上の値がBBに入る絶対連続解とする。ほとんど至る所でy′=∣y∣≥0y'=\sqrt{\lvert y\rvert}\ge0であるから、§E13.2 定理 2.5によりyyは単調非減少であり、y≥y(0)=0y\ge y(0)=0である。τ=max⁡{t∈[0,h′]:y(t)=0}\tau=\max\{t\in[0,h']:y(t)=0\}とおくと、[0,τ][0,\tau]上でy=0y=0であり、(τ,h′](\tau,h']上でy>0y>0である。τ<h′\tau<h'ならば、§E13.2 定理 3.1 (1)を開集合(0,∞)(0,\infty)上の⋅\sqrt{\cdot}と区間(τ,h′](\tau,h']に適用すると、y\sqrt yは局所絶対連続であり、ほとんど至る所で(y)′=y′/(2y)=1/2(\sqrt y)'=y'/(2\sqrt y)=1/2である。§E13.2 系 2.6により、τ<s<t≤h′\tau<s<t\le h'についてy(t)=y(s)+(t−s)/2\sqrt{y(t)}=\sqrt{y(s)}+(t-s)/2であり、s↓τs\downarrow\tauとするとy(t)=(t−τ)2/4y(t)=(t-\tau)^2/4を得る。よってy=yτ∣[0,h′]y=y_\tau|_{[0,h']}である。逆に、各ycy_cは一階連続微分可能な解であり、[0,1][0,1]上で0≤yc(t)≤t2/4≤1/40\le y_c(t)\le t^2/4\le1/4である。§E13.2 系 2.7 (2)によりycy_cは[0,1][0,1]上で局所絶対連続であるから、その[0,h′][0,h']への制限は値がBBに入る絶対連続解である。

4 微分不等式による比較

定理 4.1.t0,y0∈Rt_0,y_0\in\R、a,b>0a,b>0とし、J=[t0−a,t0+a]J=[t_0-a,t_0+a]、B=[y0−b,y0+b]B=[y_0-b,y_0+b]とおく。写像f:J×B→Rf:J\times B\to\Rは、Lebesgue 零集合N⊆JN\subseteq Jと非負関数m∈L1(J)m\in L^1(J)について Carathéodory 条件を満たすとし、0<h≤a0<h\le aは∫t0−ht0m(s) ds≤b\int_{t_0-h}^{t_0}m(s)\,ds\le bと∫t0t0+hm(s) ds≤b\int_{t_0}^{t_0+h}m(s)\,ds\le bを満たすとする。y∗,y∗y^*,y_*を[t0,t0+h][t_0,t_0+h]上の最大の解と最小の解とし、0<h′≤h0<h'\le hとする。このとき、vvが[t0,t0+h′][t_0,t_0+h']上の下解でありv(t0)≤y0v(t_0)\le y_0を満たすならば、[t0,t0+h′][t_0,t_0+h']上でv≤y∗v\le y^*である。wwが[t0,t0+h′][t_0,t_0+h']上の上解でありw(t0)≥y0w(t_0)\ge y_0を満たすならば、[t0,t0+h′][t_0,t_0+h']上でw≥y∗w\ge y_*である。

証明.vvをJJ上へ延長したvˉ:J→B\bar v:J\to Bを、t<t0t<t_0でvˉ(t)=v(t0)\bar v(t)=v(t_0)、t>t0+h′t>t_0+h'でvˉ(t)=v(t0+h′)\bar v(t)=v(t_0+h')と定める。vˉ\bar vは連続である。(t,u)∈J×B(t,u)\in J\times Bに対してg(t,u)=f(t,max⁡{u,vˉ(t)})g(t,u)=f(t,\max\{u,\bar v(t)\})とおく。各u∈Bu\in Bについてt↦max⁡{u,vˉ(t)}t\mapsto\max\{u,\bar v(t)\}はJJからBBへの連続写像であるから、§E13.4 補題 2.3によりt↦g(t,u)t\mapsto g(t,u)は可測である。t∉Nt\notin Nならばu↦g(t,u)u\mapsto g(t,u)は連続であり、∣g(t,u)∣≤m(t)\lvert g(t,u)\rvert\le m(t)である。したがってggはNNとmmについて Carathéodory 条件を満たし、定理 1.2により、[t0−h,t0+h][t_0-h,t_0+h]上の絶対連続関数zzであって、値がBBに入り、z(t0)=y0z(t_0)=y_0とほとんど至る所でのz′(t)=g(t,z(t))z'(t)=g(t,z(t))を満たすものが存在する。

φ=v−z\varphi=v-zを[t0,t0+h′][t_0,t_0+h']上で考える。§E13.2 補題 1.2 (2)によりφ\varphiは絶対連続であり、φ(t0)=v(t0)−y0≤0\varphi(t_0)=v(t_0)-y_0\le0である。§E13.2 補題 3.3 (1)により、ほとんど至る所で(φ+)′=1{v>z}(v′−z′)(\varphi^+)'=\mathbf 1_{\{v>z\}}(v'-z')である。v(t)>z(t)v(t)>z(t)を満たすほとんどすべてのttで、max⁡{z(t),vˉ(t)}=v(t)\max\{z(t),\bar v(t)\}=v(t)であるからz′(t)=f(t,v(t))≥v′(t)z'(t)=f(t,v(t))\ge v'(t)である。したがってほとんど至る所で(φ+)′≤0(\varphi^+)'\le0であり、§E13.2 定理 2.5によりφ+(t)=φ+(t0)+∫t0t(φ+)′≤0\varphi^+(t)=\varphi^+(t_0)+\int_{t_0}^t(\varphi^+)'\le0である。よって[t0,t0+h′][t_0,t_0+h']上でv≤zv\le zであり、max⁡{z,vˉ}=z\max\{z,\bar v\}=zであるから、ほとんど至る所でz′(t)=f(t,z(t))z'(t)=f(t,z(t))である。z∣[t0,t0+h′]z|_{[t_0,t_0+h']}は値がBBに入る絶対連続解であり、最大の解の定義によりz≤y∗z\le y^*である。これでv≤z≤y∗v\le z\le y^*を得る。

補題 3.3のf−f^-を取る。補題 3.3 (1)によりf−f^-はNNとmmについて Carathéodory 条件を満たす。補題 3.3 (3)により−w-wはf−f^-の[t0,t0+h′][t_0,t_0+h']上の下解であり、−w(t0)≤−y0-w(t_0)\le-y_0である。補題 3.3 (4)により、−y∗-y_*はf−f^-の[t0,t0+h][t_0,t_0+h]上の最大の解である。前段をf−f^-に適用すると−w≤−y∗-w\le-y_*、すなわちw≥y∗w\ge y_*を得る。▨

補題 4.2.t0,y0∈Rt_0,y_0\in\R、a,b>0a,b>0とし、J=[t0−a,t0+a]J=[t_0-a,t_0+a]、B=[y0−b,y0+b]B=[y_0-b,y_0+b]とおく。写像f:J×B→Rf:J\times B\to\Rは、Lebesgue 零集合N⊆JN\subseteq Jと非負関数m∈L1(J)m\in L^1(J)について Carathéodory 条件を満たし、0<h≤a0<h\le aは∫t0−ht0m(s) ds≤b\int_{t_0-h}^{t_0}m(s)\,ds\le bと∫t0t0+hm(s) ds≤b\int_{t_0}^{t_0+h}m(s)\,ds\le bを満たすとする。τ(t)=2t0−t\tau(t)=2t_0-tとおき、(t,u)∈J×B(t,u)\in J\times Bに対してf^(t,u)=−f(τ(t),u)\hat f(t,u)=-f(\tau(t),u)と定める。0<h′≤h0<h'\le hについて、写像y:[t0−h′,t0]→By:[t_0-h',t_0]\to Bに対しy^=y∘τ:[t0,t0+h′]→B\hat y=y\circ\tau:[t_0,t_0+h']\to Bとおく。このとき、次が成り立つ。

  1. f^\hat fは零集合τ(N)\tau(N)と非負関数m∘τ∈L1(J)m\circ\tau\in L^1(J)について Carathéodory 条件を満たし、hhはm∘τm\circ\tauについても∫t0−ht0m(τ(s)) ds≤b\int_{t_0-h}^{t_0}m(\tau(s))\,ds\le bと∫t0t0+hm(τ(s)) ds≤b\int_{t_0}^{t_0+h}m(\tau(s))\,ds\le bを満たす。
  2. yyが絶対連続であることとy^\hat yが絶対連続であることは同値であり、y(t0)=y0y(t_0)=y_0とy^(t0)=y0\hat y(t_0)=y_0は同値である。(t0,t0+h′)(t_0,t_0+h')のほとんど至る所でのy′(τ(t))=f(τ(t),y(τ(t)))y'(\tau(t))=f(\tau(t),y(\tau(t)))とy^′(t)=f^(t,y^(t))\hat y'(t)=\hat f(t,\hat y(t))は同値であり、ほとんど至る所でのy′(τ(t))≥f(τ(t),y(τ(t)))y'(\tau(t))\ge f(\tau(t),y(\tau(t)))とy^′(t)≤f^(t,y^(t))\hat y'(t)\le\hat f(t,\hat y(t))は同値であり、ほとんど至る所でのy′(τ(t))≤f(τ(t),y(τ(t)))y'(\tau(t))\le f(\tau(t),y(\tau(t)))とy^′(t)≥f^(t,y^(t))\hat y'(t)\ge\hat f(t,\hat y(t))は同値である。

証明.τ\tauはJJからJJへの全単射であり、長さを保つ。したがって、τ\tauは Lebesgue 可測集合を Lebesgue 可測集合へ、零集合を零集合へ写す。また、Lebesgue 可測集合E⊆JE\subseteq Jと非負可積分関数c:J→Rc:J\to\Rに対して、∫Ec∘τ=∫τ(E)c\int_Ec\circ\tau=\int_{\tau(E)}cが成り立つ。

主張 (1)を示す。各u∈Bu\in Bについてt↦f^(t,u)=−f(τ(t),u)t\mapsto\hat f(t,u)=-f(\tau(t),u)は可測である。t∈J∖τ(N)t\in J\setminus\tau(N)ならばτ(t)∈J∖N\tau(t)\in J\setminus Nであるから、u↦f^(t,u)u\mapsto\hat f(t,u)は連続であり、∣f^(t,u)∣≤m(τ(t))\lvert\hat f(t,u)\rvert\le m(\tau(t))である。∫t0t0+hm(τ(s)) ds=∫t0−ht0m(s) ds≤b\int_{t_0}^{t_0+h}m(\tau(s))\,ds=\int_{t_0-h}^{t_0}m(s)\,ds\le bであり、同様に∫t0−ht0m(τ(s)) ds=∫t0t0+hm(s) ds≤b\int_{t_0-h}^{t_0}m(\tau(s))\,ds=\int_{t_0}^{t_0+h}m(s)\,ds\le bである。

主張 (2)を示す。τ\tauは[t0,t0+h′][t_0,t_0+h']と[t0−h′,t0][t_0-h',t_0]の間の等長写像であるから、yyの絶対連続性の定義のε\epとδ\deltaの組はy^\hat yについても同じ条件を満たし、逆も同様である。y^(t0)=y(t0)\hat y(t_0)=y(t_0)である。yyがτ(t)\tau(t)で微分可能であることとy^\hat yがttで微分可能であることは同値であり、そのときy^′(t)=−y′(τ(t))\hat y'(t)=-y'(\tau(t))である。またf^(t,y^(t))=−f(τ(t),y(τ(t)))\hat f(t,\hat y(t))=-f(\tau(t),y(\tau(t)))であるから、y′(τ(t))=f(τ(t),y(τ(t)))y'(\tau(t))=f(\tau(t),y(\tau(t)))とy^′(t)=f^(t,y^(t))\hat y'(t)=\hat f(t,\hat y(t))は同値であり、y′(τ(t))≥f(τ(t),y(τ(t)))y'(\tau(t))\ge f(\tau(t),y(\tau(t)))とy^′(t)≤f^(t,y^(t))\hat y'(t)\le\hat f(t,\hat y(t))は同値であり、y′(τ(t))≤f(τ(t),y(τ(t)))y'(\tau(t))\le f(\tau(t),y(\tau(t)))とy^′(t)≥f^(t,y^(t))\hat y'(t)\ge\hat f(t,\hat y(t))は同値である。▨

系 4.3.補題 4.2のt0t_0、y0y_0、aa、bb、JJ、BB、ff、NN、mm、hh、τ\tau、f^\hat fをそのまま用い、0<h′≤h0<h'\le hとする。このとき、f^\hat fの[t0,t0+h][t_0,t_0+h]上の最大の解y^∗\hat y^*と最小の解y^∗\hat y_*がそれぞれただ一つ存在し、次が成り立つ。

  1. v:[t0−h′,t0]→Bv:[t_0-h',t_0]\to Bが絶対連続関数であり、v(t0)≤y0v(t_0)\le y_0と(t0−h′,t0)(t_0-h',t_0)のほとんど至る所でのv′(t)≥f(t,v(t))v'(t)\ge f(t,v(t))を満たすならば、[t0−h′,t0][t_0-h',t_0]上でv(t)≤y^∗(τ(t))v(t)\le\hat y^*(\tau(t))である。
  2. w:[t0−h′,t0]→Bw:[t_0-h',t_0]\to Bが絶対連続関数であり、w(t0)≥y0w(t_0)\ge y_0と(t0−h′,t0)(t_0-h',t_0)のほとんど至る所でのw′(t)≤f(t,w(t))w'(t)\le f(t,w(t))を満たすならば、[t0−h′,t0][t_0-h',t_0]上でw(t)≥y^∗(τ(t))w(t)\ge\hat y_*(\tau(t))である。

証明.補題 4.2 (1)により、f^\hat fは定理 3.4の仮定を満たし、y^∗\hat y^*とy^∗\hat y_*がそれぞれただ一つ存在する。

主張 (1)を示す。v^=v∘τ\hat v=v\circ\tauは[t0,t0+h′][t_0,t_0+h']上の絶対連続関数であり、値はBBに入り、v^(t0)=v(t0)≤y0\hat v(t_0)=v(t_0)\le y_0である。補題 4.2 (2)により、ほとんど至る所でv^′(t)≤f^(t,v^(t))\hat v'(t)\le\hat f(t,\hat v(t))であるから、v^\hat vはf^\hat fの下解である。定理 4.1をf^\hat fとv^\hat vに適用するとv^≤y^∗\hat v\le\hat y^*であり、v(t)=v^(τ(t))≤y^∗(τ(t))v(t)=\hat v(\tau(t))\le\hat y^*(\tau(t))である。

主張 (2)を示す。w^=w∘τ\hat w=w\circ\tauとおくと、w^(t0)≥y0\hat w(t_0)\ge y_0であり、補題 4.2 (2)により、ほとんど至る所でw^′(t)≥f^(t,w^(t))\hat w'(t)\ge\hat f(t,\hat w(t))である。w^\hat wはf^\hat fの上解であり、定理 4.1によりw^≥y^∗\hat w\ge\hat y_*であり、w(t)=w^(τ(t))≥y^∗(τ(t))w(t)=\hat w(\tau(t))\ge\hat y_*(\tau(t))である。▨

注意 4.4.系 4.3の記号を用い、0<h′≤h0<h'\le hとし、y:[t0−h′,t0]→By:[t_0-h',t_0]\to Bを、y(t0)=y0y(t_0)=y_0と(t0−h′,t0)(t_0-h',t_0)のほとんど至る所でのy′(t)=f(t,y(t))y'(t)=f(t,y(t))を満たす絶対連続関数とする。yyは系 4.3 (1)のvvと系 4.3 (2)のwwの仮定をともに満たすから、[t0−h′,t0][t_0-h',t_0]上でy^∗(τ(t))≤y(t)≤y^∗(τ(t))\hat y_*(\tau(t))\le y(t)\le\hat y^*(\tau(t))である。補題 4.2 (2)により、y^∗∘τ\hat y^*\circ\tauとy^∗∘τ\hat y_*\circ\tauの[t0−h′,t0][t_0-h',t_0]への制限は、それ自身がyyの条件を満たす。したがって、[t0−h′,t0][t_0-h',t_0]上でyyの条件を満たす関数のうちで最大のものと最小のものが存在し、それぞれy^∗∘τ\hat y^*\circ\tauとy^∗∘τ\hat y_*\circ\tauの制限である。

5 演習

問題 5.1.t0=0t_0=0、y0=0y_0=0、a=b=1a=b=1、J=B=[−1,1]J=B=[-1,1]、f(t,u)=∣u∣f(t,u)=\sqrt{\lvert u\rvert}、m≡1m\equiv1、h=1h=1とする。補題 1.1の近似解yny_nが、すべてのn∈N≥1n\in\NNについて[−1,1][-1,1]上で恒等的に00であることを示せ。次に、0<h′≤10<h'\le1とし、y:[−h′,0]→By:[-h',0]\to Bをy(0)=0y(0)=0と(−h′,0)(-h',0)のほとんど至る所でのy′(t)=∣y(t)∣y'(t)=\sqrt{\lvert y(t)\rvert}を満たす絶対連続関数とするとき、[−h′,0][-h',0]上で−t2/4≤y(t)≤0-t^2/4\le y(t)\le0が成り立つことを示し、00と−t2/4-t^2/4がともにこの条件を満たすことを確かめよ。最後に、近似解の一様極限が、[0,1][0,1]上では最小の解であり、[−1,0][-1,0]上ではこの条件を満たす関数のうちで最大のものであることを結論せよ。

解答.

定数写像0:[−1,1]→B0:[-1,1]\to Bは連続であり、s↦f(s,0)=0s\mapsto f(s,0)=0は可積分であり、すべてのttについて0=0+∫0σn(t)∣0∣ ds0=0+\int_0^{\sigma_n(t)}\sqrt{\lvert0\rvert}\,dsである。補題 1.1 (1)の一意性によりyn≡0y_n\equiv0である。したがって近似解の列は00へ一様収束する。

yyを後半の条件を満たす関数とする。ほとんど至る所でy′≥0y'\ge0であるから、§E13.2 定理 2.5によりyyは単調非減少であり、[−h′,0][-h',0]上でy(t)≤y(0)=0y(t)\le y(0)=0である。σ=min⁡{t∈[−h′,0]:y(t)=0}\sigma=\min\{t\in[-h',0]:y(t)=0\}とおくと、[σ,0][\sigma,0]上でy=0y=0であり、[−h′,σ)[-h',\sigma)上でy<0y<0である。σ>−h′\sigma>-h'ならば、§E13.2 定理 3.1 (1)を開集合(−∞,0)(-\infty,0)上のu↦−uu\mapsto\sqrt{-u}と区間[−h′,σ)[-h',\sigma)に適用すると、−y\sqrt{-y}は局所絶対連続であり、ほとんど至る所で(−y)′=−y′/(2−y)=−1/2(\sqrt{-y})'=-y'/(2\sqrt{-y})=-1/2である。§E13.2 系 2.6により、−h′≤s<t<σ-h'\le s<t<\sigmaについて−y(t)=−y(s)−(t−s)/2\sqrt{-y(t)}=\sqrt{-y(s)}-(t-s)/2であり、t↑σt\uparrow\sigmaとするとy(s)=−(σ−s)2/4y(s)=-(\sigma-s)^2/4を得る。0≤σ−s≤−s0\le\sigma-s\le-sであるからy(s)≥−s2/4y(s)\ge-s^2/4である。[σ,0][\sigma,0]上でもy=0≥−t2/4y=0\ge-t^2/4である。よって[−h′,0][-h',0]上で−t2/4≤y(t)≤0-t^2/4\le y(t)\le0である。

00は条件を満たす。y(t)=−t2/4y(t)=-t^2/4は一階連続微分可能であり、t≤0t\le0でy′(t)=−t/2=t2/4=∣y(t)∣y'(t)=-t/2=\sqrt{t^2/4}=\sqrt{\lvert y(t)\rvert}、y(0)=0y(0)=0、y([−1,0])⊆[−1/4,0]⊆By([-1,0])\subseteq[-1/4,0]\subseteq Bである。§E13.2 系 2.7 (2)によりyyは[−1,0][-1,0]上で局所絶対連続であるから、その[−h′,0][-h',0]への制限は絶対連続であり、条件を満たす。

近似解の一様極限は00である。例 3.5により、00は[0,1][0,1]上の最小の解である。前段により、00は[−1,0][-1,0]の各部分区間[−h′,0][-h',0]の上で条件を満たす関数のうちで最大である。補題 4.2 (2)により、この条件を満たす関数yyの全体は、f^(t,u)=−∣u∣\hat f(t,u)=-\sqrt{\lvert u\rvert}の[0,h′][0,h']上の値がBBに入る絶対連続解y(−t)y(-t)の全体に対応し、各yyはy(−t)≤0y(-t)\le0を満たす。00はf^\hat fの解であるから、00はf^\hat fの[0,1][0,1]上の最大の解である。▨

問題 5.2.f(t,u)=t2+u2f(t,u)=t^2+u^2とする。t0=0t_0=0、y0=0y_0=0、a=b=1a=b=1、J=B=[−1,1]J=B=[-1,1]について、m≡2m\equiv2、h=1/2h=1/2、Lipschitz 係数ℓ≡2\ell\equiv2を取ることができることを確かめ、[0,1/2][0,1/2]上の値がBBに入る絶対連続解yyがただ一つ存在して[0,1/2][0,1/2]上でt3/3≤y(t)≤tan⁡tt^3/3\le y(t)\le\tan tを満たすことを示せ。次に、T>0T>0と、y(0)=0y(0)=0かつすべてのt∈[0,T)t\in[0,T)でy′(t)=t2+y(t)2y'(t)=t^2+y(t)^2を満たす一階連続微分可能な関数y:[0,T)→Ry:[0,T)\to\Rが与えられたとき、T≤1+π/2T\le1+\pi/2であることを示せ。

解答.

(t,u)∈[−1,1]2(t,u)\in[-1,1]^2ならば∣f(t,u)∣=t2+u2≤2\lvert f(t,u)\rvert=t^2+u^2\le2であり、u,v∈[−1,1]u,v\in[-1,1]ならば∣f(t,u)−f(t,v)∣=∣u+v∣∣u−v∣≤2∣u−v∣\lvert f(t,u)-f(t,v)\rvert=\lvert u+v\rvert\lvert u-v\rvert\le2\lvert u-v\rvertである。ffは連続であるから、N=∅N=\varnothing、m≡2m\equiv2について Carathéodory 条件を満たし、ℓ≡2\ell\equiv2は可積分な Lipschitz 係数である。h=1/2h=1/2についてt0t_0の左右のmmの積分はともに1=b1=bである。定理 3.4により、[0,1/2][0,1/2]上の値がBBに入る絶対連続解はただ一つであり、それは最大の解y∗y^*にも最小の解y∗y_*にも等しい。これをyyと書く。

v(t)=t3/3v(t)=t^3/3は[0,1/2][0,1/2]上で一階連続微分可能であり、値は[0,1/24]⊆B[0,1/24]\subseteq Bに入り、v′(t)=t2≤f(t,v(t))v'(t)=t^2\le f(t,v(t))、v(0)=0v(0)=0を満たす。§E13.2 系 2.7 (2)によりvvは有限閉区間[0,1/2][0,1/2]上で局所絶対連続であるから絶対連続であり、下解である。定理 4.1によりv≤y∗=yv\le y^*=yである。w(t)=tan⁡tw(t)=\tan tは[0,1/2][0,1/2]上で一階連続微分可能であり、0≤tan⁡t≤tan⁡(1/2)<tan⁡(π/4)=10\le\tan t\le\tan(1/2)<\tan(\pi/4)=1であるから値はBBに入る。t2≤1t^2\le1であるからw′(t)=1+w(t)2≥f(t,w(t))w'(t)=1+w(t)^2\ge f(t,w(t))であり、w(0)=0w(0)=0である。vvと同じく§E13.2 系 2.7 (2)によりwwは絶対連続であるから、wwは上解であり、定理 4.1によりw≥y∗=yw\ge y_*=yである。

y′≥0y'\ge0であるからyyは単調非減少であり、[0,T)[0,T)上でy≥0y\ge0である。T≤1T\le1ならばT≤1+π/2T\le1+\pi/2である。T>1T>1とする。arctan⁡y\arctan yは[0,T)[0,T)上で一階連続微分可能であり、t∈[1,T)t\in[1,T)ではy′(t)=t2+y(t)2≥1+y(t)2y'(t)=t^2+y(t)^2\ge1+y(t)^2であるから(arctan⁡y)′(t)=y′(t)/(1+y(t)2)≥1(\arctan y)'(t)=y'(t)/(1+y(t)^2)\ge1である。平均値の定理により、t∈[1,T)t\in[1,T)についてarctan⁡y(t)≥arctan⁡y(1)+(t−1)≥t−1\arctan y(t)\ge\arctan y(1)+(t-1)\ge t-1である。arctan⁡y(t)<π/2\arctan y(t)<\pi/2であるから、すべてのt∈[1,T)t\in[1,T)でt<1+π/2t<1+\pi/2であり、T≤1+π/2T\le1+\pi/2である。▨

参考文献

  1. Philip Hartman, Ordinary Differential Equations, 2nd ed., Classics in Applied Mathematics, Society for Industrial and Applied Mathematics, Philadelphia, 2002.Peano の存在定理、Carathéodory 条件の下の存在定理、最大の解と最小の解、微分不等式を参考にした。
  2. Wolfgang Walter, Ordinary Differential Equations, translated by Russell Thompson, Graduate Texts in Mathematics 182, Springer, 1998.Tonelli の近似による存在定理の証明と微分不等式による比較を参考にした。

前提記事

9 本の記事・単元を表示