Picard–Lindelöfの定理(Picard–Lindelöf theorem)とは、初期値問題 $x'=f(t,x)$、$x(t_0)=x_0$ について、$f$ が連続で $x$ について Lipschitz 条件をみたせば、$t_0$ の近くで解がただ 1 つ存在するという定理である。$[t_0-T,t_0+T]\times\overline{B}_\delta(x_0)$ の上で $\lvert f\rvert\le M$ なら、解は少なくとも $\lvert t-t_0\rvert\le\min\{T,\delta/M\}$ で存在する。証明は積分方程式の右辺を連続関数の空間の縮小写像とみて不動点をとるもので、解は Picard の逐次近似の一様極限である。Lipschitz 条件を外すと一意性が崩れる($x'=3x^{2/3}$)。
前提知識: 常微分方程式, Lipschitz連続, 完備距離空間, 微分積分学の基本定理
Picard–Lindelöf の定理(Picard–Lindelöf theorem)は、常微分方程式の初期値問題
$$
x'=f(t,x),\qquad x(t_0)=x_0
$$
の解が、右辺 $f$ が連続で $x$ について Lipschitz 条件をみたすとき、$t_0$ の近くにただ 1 つ存在することを述べる定理である。証明の考え方は、微分方程式を積分方程式 $x(t)=x_0+\int_{t_0}^tf(s,x(s))\,ds$ に書き直し、右辺が定める写像 $K$ を連続関数の空間の上の縮小写像とみて、その不動点として解を得ることである。不動点は $x_0$ から始めて $K$ を繰り返し当てた列(Picard の逐次近似)の極限として得られ、その誤差も評価できる。この定理は、常微分方程式の理論のすべて(最大解、初期値への依存性、ベクトル場の流れ)の出発点になっている。
$\mathbb{R}^n$ のベクトルの長さを $\lvert x\rvert$(Euclid ノルム)で表し、$\overline{B}_\delta(x_0):=\{x\in\mathbb{R}^n\mid\lvert x-x_0\rvert\le\delta\}$ と書く。$U\subset\mathbb{R}\times\mathbb{R}^n$ を開集合、$f\colon U\to\mathbb{R}^n$ を連続写像とする。
$I$ を $t_0$ を含む区間(開区間でも閉区間でもよい)とする。$C^1$ 級の写像 $x\colon I\to\mathbb{R}^n$ で、すべての $t\in I$ について $(t,x(t))\in U$ かつ $x'(t)=f(t,x(t))$ となり、$x(t_0)=x_0$ をみたすものを、初期値問題 $x'=f(t,x)$、$x(t_0)=x_0$ の $I$ 上の解という。$I$ が端点を含むときは、端点での微分は片側微分係数とする。
$V\subset U$ とする。$f$ が $V$ の上で $x$ について Lipschitz 条件をみたすとは、定数 $L\ge0$ があって
$$
\lvert f(t,x)-f(t,y)\rvert\le L\,\lvert x-y\rvert\qquad\bigl((t,x),(t,y)\in V\bigr)
$$
となることをいう。$L$ を Lipschitz 定数という。$U$ の各点 $(t_0,x_0)$ に対し、$\eta>0$、$r>0$ で $[t_0-\eta,t_0+\eta]\times\overline{B}_r(x_0)\subset U$ となり、その上で $f$ が $x$ について Lipschitz 条件をみたすものがあるとき、$f$ は $x$ について局所 Lipschitz 条件をみたすという。
局所 Lipschitz 条件は常微分方程式の記事の定義「局所 Lipschitz 条件」と同じものである。$f$ が $U$ 上 $C^1$ 級なら、この条件はみたされる(Lipschitz連続の記事の命題「連続微分可能な写像の局所Lipschitz性」の議論を、$t$ を止めて $x$ について当てはめる)。
解 $x$ は連続なので $s\mapsto f(s,x(s))$ も連続であり、微分積分学の基本定理により、$x$ が解であることと、$x$ が連続で積分方程式
$$
x(t)=x_0+\int_{t_0}^{t}f\bigl(s,x(s)\bigr)\,ds\qquad(t\in I)
$$
をみたすことは同値である(常微分方程式の記事の補題「積分方程式への書き換え」。閉区間の端点でも、片側微分について同じ議論が成り立つ)。以下、この積分方程式を扱う。
$T>0$、$\delta>0$ を
$$
V:=[t_0-T,t_0+T]\times\overline{B}_\delta(x_0)\subset U
$$
となるようにとり、$f$ は $V$ の上で $x$ について Lipschitz 定数 $L$ の Lipschitz 条件をみたすとする。$M:=\max_{V}\lvert f\rvert$ とおき($V$ はコンパクトなので最大値がある)、
$$
T_0:=\min\Bigl\{T,\ \frac{\delta}{M}\Bigr\}\qquad(M=0\text{ のときは }T_0:=T),\qquad I_0:=[t_0-T_0,t_0+T_0]
$$
とする。このとき初期値問題 $x'=f(t,x)$、$x(t_0)=x_0$ の $I_0$ 上の解 $x$ がただ 1 つ存在し、$x(I_0)\subset\overline{B}_\delta(x_0)$ である。
連続性だけでは一意性は成り立たない。$x'=\sqrt{\lvert x\rvert}$、$x(0)=0$ は $x\equiv0$ と $x(t)=t^2/4$($t\ge0$)の二つの解をもち、$\sqrt{\lvert x\rvert}$ は $x=0$ の近くで Lipschitz 連続でない。Lipschitz 条件は一意性と、証明で使う縮小性の両方に効いている。また球の半径 $\delta$ と上界 $M$ から時間 $T_0=\min\{T,\delta/M\}$ を決めるのは、解が球から出ないことを保証するためである。
この形では、解の存在する時間 $T_0$ が $V$、$M$ だけで具体的に決まり、Lipschitz 定数 $L$ によらない。$L$ は一意性と証明の中の縮小の定数にだけ現れる。局所 Lipschitz 条件の下での形は次のとおりである。
$f$ が $x$ について局所 Lipschitz 条件をみたすなら、各 $(t_0,x_0)\in U$ に対してある $\varepsilon>0$ があり、初期値問題 $x'=f(t,x)$、$x(t_0)=x_0$ の解が開区間 $(t_0-\varepsilon,t_0+\varepsilon)$ 上に存在する。さらに、同じ初期値問題の 2 つの解は、定義域の共通部分の上で一致する。
局所 Lipschitz 条件の $\eta$、$r$ をとり、thm-pl-main を $T:=\eta$、$\delta:=r$ で使えば、$T_0>0$ について $[t_0-T_0,t_0+T_0]$ 上の解が得られる。$\varepsilon:=T_0$ としてその内部に制限すればよい。後半は常微分方程式の記事の系「解の一意性」である(そこでは、各点の近くで Grönwall の不等式から 2 つの解の差が $0$ であることを示し、一致する時刻の集合が開かつ閉であることを使っている)。$\square$
これで、常微分方程式の記事の定理「Picard–Lindelöf の定理」と、Lipschitz連続の記事の定理「初期値問題の解の存在と一意性」の存在の部分が証明されたことになる。
証明は、連続関数の空間に重み付きの距離を入れて、積分方程式の右辺を縮小写像にすることで行う。使う道具は完備距離空間の記事の定理「縮小写像の不動点定理」(空でない完備距離空間の縮小写像はただ 1 つの不動点をもつ)と、同じ記事の命題「有界関数と有界連続関数の空間の完備性」である。
thm-pl-main の記号で、$y$ を初期値問題 $x'=f(t,x)$、$x(t_0)=x_0$ の $I_0$ 上の任意の解とする。このときすべての $t\in I_0$ で $\lvert y(t)-x_0\rvert\le\delta$ である。
$t\ge t_0$ の側を示す($t\le t_0$ の側も同様)。$S:=\{s\in[t_0,t_0+T_0]\mid\text{すべての }r\in[t_0,s]\text{ で }\lvert y(r)-x_0\rvert\le\delta\}$ とおき、$s^\ast:=\sup S$ とする。$t_0\in S$ であり、$y$ の連続性から $s^\ast\in S$ である。$s^\ast< t_0+T_0$ と仮定する。$[t_0,s^\ast]$ の上では $(r,y(r))\in V$ なので $\lvert f(r,y(r))\rvert\le M$ であり、積分方程式から
$$
\lvert y(s^\ast)-x_0\rvert\le\int_{t_0}^{s^\ast}\bigl\lvert f(r,y(r))\bigr\rvert\,dr\le M(s^\ast-t_0)< MT_0\le\delta
$$
である($M=0$ なら左辺は $0<\delta$)。連続性から $s^\ast$ の少し先まで $\lvert y-x_0\rvert<\delta$ が続くので、$s^\ast=\sup S$ に反する。よって $s^\ast=t_0+T_0$ である。$\square$
$T_0\le\delta/M$ という制限は、ここで「速さ $M$ 以下で動く点は、時間 $\delta/M$ の間は半径 $\delta$ の球から出られない」という形で効いている。
段 1(空間).$L>0$ としてよい(Lipschitz 定数は大きくとり直してもよい)。$X:=C(I_0,\mathbb{R}^n)$ を $I_0$ 上の連続写像全体とし、
$$
\|x\|_w:=\max_{t\in I_0}e^{-2L\lvert t-t_0\rvert}\,\lvert x(t)\rvert,\qquad d_w(x,y):=\|x-y\|_w
$$
とおく。最大値ノルム $\|x\|_\infty:=\max_{I_0}\lvert x\rvert$ との間に $e^{-2LT_0}\|x\|_\infty\le\|x\|_w\le\|x\|_\infty$ が成り立つので、$d_w$ の Cauchy 列は $\|\cdot\|_\infty$ の Cauchy 列であり、収束も同じである。$(X,\|\cdot\|_\infty)$ は完備である(完備距離空間の記事の命題「有界関数と有界連続関数の空間の完備性」を各成分に使う)から、$(X,d_w)$ も完備である。
$$
C:=\bigl\{x\in X\ \big|\ \text{すべての }t\in I_0\text{ で }\lvert x(t)-x_0\rvert\le\delta\bigr\}
$$
は定数写像 $x_0$ を含み、$d_w$ について閉じている(各点での不等式は一様収束の極限で保たれる)。閉部分集合なので $(C,d_w)$ は空でない完備距離空間である。
段 2(写像 $K$).$x\in C$ なら $(s,x(s))\in V\subset U$ なので、
$$
K(x)(t):=x_0+\int_{t_0}^{t}f\bigl(s,x(s)\bigr)\,ds\qquad(t\in I_0)
$$
が定義でき、$K(x)$ は連続である。さらに $\lvert K(x)(t)-x_0\rvert\le M\lvert t-t_0\rvert\le MT_0\le\delta$ なので $K(x)\in C$ である。
段 3(縮小).$x,y\in C$、$t\ge t_0$ とすると、Lipschitz 条件と $\lvert x(s)-y(s)\rvert\le e^{2L(s-t_0)}\|x-y\|_w$ から
$$
\lvert K(x)(t)-K(y)(t)\rvert\le\int_{t_0}^{t}L\,\lvert x(s)-y(s)\rvert\,ds\le L\,\|x-y\|_w\int_{t_0}^{t}e^{2L(s-t_0)}\,ds\le\frac12\,e^{2L(t-t_0)}\,\|x-y\|_w
$$
である。$t\le t_0$ でも積分の向きを変えて同じ評価 $\frac12e^{2L\lvert t-t_0\rvert}\|x-y\|_w$ を得る。両辺に $e^{-2L\lvert t-t_0\rvert}$ を掛けて $t$ について最大をとると
$$
\|K(x)-K(y)\|_w\le\frac12\,\|x-y\|_w
$$
であり、$K\colon C\to C$ は定数 $\frac12$ の縮小写像である。
段 4(存在).縮小写像の不動点定理により、$K(x)=x$ となる $x\in C$ がただ 1 つある。$x$ は積分方程式をみたす連続写像なので、上で見た同値により $I_0$ 上の解であり、$x(I_0)\subset\overline{B}_\delta(x_0)$ である。
段 5(一意性).$y$ を $I_0$ 上の任意の解とする。lem-pl-stay により $y\in C$ であり、積分方程式から $K(y)=y$ である。不動点はただ 1 つなので $y=x$ である。$\square$
重み $e^{-2L\lvert t-t_0\rvert}$ を入れない最大値ノルムでは、段 3 の評価が $\|K(x)-K(y)\|_\infty\le LT_0\|x-y\|_\infty$ になり、縮小写像にするには $T_0<1/L$ という余分な制限が要る。重みは、時刻 $t_0$ から離れるほど誤差が指数関数的に広がりうることを見込んで、遠くの誤差を割り引いて測るためのものである。
縮小写像の不動点定理の証明は、任意の点から写像を繰り返し当てた列が不動点に収束することを示している。$K$ に当てはめたものが Picard の逐次近似である。
thm-pl-main の記号で、定数写像 $x_0(t):=x_0$ から始めて
$$
x_{m+1}(t):=x_0+\int_{t_0}^{t}f\bigl(s,x_m(s)\bigr)\,ds\qquad(t\in I_0,\ m\ge0)
$$
で定めた列 $(x_m)$ を Picard の逐次近似(Picard iteration)という。段 2 により各 $x_m$ は $C$ に属するので、この定義は意味をもつ。
$x$ を thm-pl-main の解とする。すべての $m\ge0$ と $t\in I_0$ について
$$
\lvert x(t)-x_m(t)\rvert\le\frac{M\,L^m\,\lvert t-t_0\rvert^{m+1}}{(m+1)!}
$$
が成り立つ。とくに $x_m$ は $I_0$ 上で $x$ に一様収束する。
$t\ge t_0$ の側を $m$ についての帰納法で示す($t\le t_0$ も同様)。$m=0$ では、積分方程式と $(s,x(s))\in V$ から $\lvert x(t)-x_0\rvert\le\int_{t_0}^t\lvert f(s,x(s))\rvert\,ds\le M(t-t_0)$ である。$m$ で成り立つとすると、$x$ と $x_m$ はともに $C$ に属するので Lipschitz 条件が使え、
$$
\lvert x(t)-x_{m+1}(t)\rvert\le\int_{t_0}^{t}L\,\lvert x(s)-x_m(s)\rvert\,ds\le\int_{t_0}^{t}\frac{M\,L^{m+1}(s-t_0)^{m+1}}{(m+1)!}\,ds=\frac{M\,L^{m+1}(t-t_0)^{m+2}}{(m+2)!}
$$
である。右辺の最大値 $ML^mT_0^{m+1}/(m+1)!$ は $m\to\infty$ で $0$ に近づく(階乗は指数関数より速く増える)ので、収束は一様である。$\square$
$x'=x$、$x(0)=1$($n=1$、$f(t,x)=x$、$L=1$)の逐次近似は
$$
x_0=1,\quad x_1(t)=1+t,\quad x_2(t)=1+t+\frac{t^2}{2},\quad\dots,\quad x_m(t)=\sum_{j=0}^{m}\frac{t^j}{j!}
$$
であり、$e^t$ の Taylor 多項式が順に現れる。実際、$x_m$ がこの形なら、項ごとに積分して
$$
x_{m+1}(t)=1+\int_0^t\sum_{j=0}^{m}\frac{s^j}{j!}\,ds=1+\sum_{j=0}^{m}\frac{t^{j+1}}{(j+1)!}=\sum_{j=0}^{m+1}\frac{t^j}{j!}
$$
となる。$V=[-T,T]\times[1-\delta,1+\delta]$ では $M=1+\delta$ なので、prop-pl-error は $\lvert e^t-x_m(t)\rvert\le(1+\delta)\lvert t\rvert^{m+1}/(m+1)!$($\lvert t\rvert\le T_0$)を与え、これは Taylor の定理の剰余の評価と同じ形である。
$x'=1+x^2$、$x(0)=0$ の解は $x(t)=\tan t$ であり、最大存在区間は $(-\pi/2,\pi/2)$ である($\tan'=1+\tan^2$)。thm-pl-main を $U=\mathbb{R}\times\mathbb{R}$、$V=[-T,T]\times[-\delta,\delta]$ で使うと $M=1+\delta^2$ なので、$T$ を大きくとっても保証される時間は
$$
T_0=\min\Bigl\{T,\ \frac{\delta}{1+\delta^2}\Bigr\}\le\frac12
$$
(最大は $\delta=1$ のとき)にとどまる。定理は解の存在する時間の下からの評価を与えるだけで、実際の最大存在区間はずっと長いことがある。区間を延ばすには、終点を新しい初期値として定理を当て直す(常微分方程式の記事の定理「最大解の存在と一意性」)。
$f$ がすべての $x$ で定義され、Lipschitz 定数が $x$ によらずにとれるなら、球 $\overline{B}_\delta(x_0)$ に閉じ込める必要がなくなり、解は時間の区間全体で存在する。
$J$ を区間、$t_0\in J$、$f\colon J\times\mathbb{R}^n\to\mathbb{R}^n$ を連続写像とし、$J$ に含まれる各有界閉区間 $I$ に対して定数 $L_I$ があって
$$
\lvert f(t,x)-f(t,y)\rvert\le L_I\,\lvert x-y\rvert\qquad(t\in I,\ x,y\in\mathbb{R}^n)
$$
となるとする。このとき各 $x_0\in\mathbb{R}^n$ について、初期値問題 $x'=f(t,x)$、$x(t_0)=x_0$ の $J$ 全体の上の解がただ 1 つ存在する。
$t_0$ を含む有界閉区間 $I\subset J$ を固定し、$L:=\max\{L_I,1\}$ とする。prf-pl-main と同じ重み付きノルム $\|x\|_w:=\max_{t\in I}e^{-2L\lvert t-t_0\rvert}\lvert x(t)\rvert$ を $X:=C(I,\mathbb{R}^n)$ に入れると、$(X,d_w)$ は完備である。今度は $f$ が $I\times\mathbb{R}^n$ 全体で定義されているので、$K(x)(t):=x_0+\int_{t_0}^tf(s,x(s))\,ds$ は $X$ 全体から $X$ への写像であり、段 3 と同じ計算で $\|K(x)-K(y)\|_w\le\frac12\|x-y\|_w$ である。よって $K$ はただ 1 つの不動点 $x_I\in X$ をもち、それが $I$ 上の解である。また $I$ 上の任意の解は $K$ の不動点なので $x_I$ に等しい。
$I\subset I'$ なら、$x_{I'}$ の $I$ への制限は $I$ 上の解なので $x_I$ に等しい。$J$ は $t_0$ を含む有界閉区間の増大列の和集合なので、$x(t):=x_I(t)$($t\in I$)は $I$ の選び方によらずに $J$ 上の写像を定め、各点の近くで解に一致するから $J$ 上の解である。$J$ 上の 2 つの解は各 $I$ の上で $x_I$ に一致するので等しい。$\square$
線形の方程式 $x'=A(t)x+b(t)$($A$、$b$ は $J$ 上連続)では、$I$ の上で $L_I:=\max_{t\in I}\|A(t)\|_{\mathrm{op}}$ とおけば仮定がみたされる。したがって線形の方程式の解は $J$ 全体で存在する。これは常微分方程式の記事の定理「線形方程式の解の大域存在と解空間」の (i) の別証明になっている(そちらでは爆発の二者択一と Grönwall の不等式を使っている)。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| $x$ についての Lipschitz 条件 | $x'=3x^{2/3}$、$x(0)=0$ | 解の一意性 |
| $f$ の $t$ についての連続性 | $x'=H(t)$、$x(0)=0$($H$ は下の関数) | 解の存在 |
| 大域的な Lipschitz 条件(局所的なものだけにする) | $x'=x^2$、$x(0)=1$ | 区間全体での存在 |
| 時間の制限 $T_0\le\delta/M$ | $x'=1$、$x(0)=0$ | 解が $\overline{B}_\delta(x_0)$ にとどまること |
1 行目:$f(x)=3x^{2/3}$($x^{2/3}:=(x^{1/3})^2$)は連続だが、$x>0$ で $\lvert f(x)-f(0)\rvert/\lvert x\rvert=3x^{-1/3}$ が $x\to+0$ で有界でないので、$0$ を含むどの $V$ の上でも Lipschitz 条件をみたさない(満たさない性質)。$x\equiv0$ と $x(t)=t^3$ がともに解なので、一意性が崩れる(常微分方程式の記事の例「反例:一意性が崩れる」)。破れる含意は「$f$ が連続なら初期値問題の解は一意である」であり、prf-pl-main の段 5 で使った不動点の一意性は、縮小の評価(段 3)がなければ得られない。
2 行目:$H(t):=0$($t\le0$)、$:=1$($t>0$)とし、$f(t,x):=H(t)$ とおく。$f$ は有界で、$x$ について Lipschitz 定数 $0$ の Lipschitz 条件をみたす(満たす性質)が、$t=0$ で連続でない(満たさない性質)。$\varepsilon>0$ と $[0,\varepsilon)$ を含む区間の上の解 $x$ があったとすると、$0< s< t<\varepsilon$ で平均値の定理から $x(t)-x(s)=t-s$ であり、$s\to+0$ として $x(t)=t$ を得る。すると $x$ の $0$ での右微分係数は $1$ だが、方程式は $x'(0)=H(0)=0$ を要求するので矛盾する。破れる含意は「$f$ が $x$ について Lipschitz 条件をみたせば解が存在する」であり、$f$ の連続性は Lipschitz 条件とは別に必要である。$x$ について不連続な例は常微分方程式の記事の例「反例:右辺が連続でないと解がない」にある(こちらは Lipschitz 条件も破る)。
3 行目:$f(x)=x^2$ は各有界閉区間で Lipschitz だが、$\lvert x^2-y^2\rvert=\lvert x+y\rvert\,\lvert x-y\rvert$ の係数 $\lvert x+y\rvert$ は有界でないので、thm-pl-global の仮定をみたさない(満たさない性質)。解 $1/(1-t)$ は $t\to1$ で発散し、$J=\mathbb{R}$ 全体では存在しない。破れる含意は「局所 Lipschitz なら解は区間全体で存在する」である。
4 行目:$f\equiv1$ では $M=1$、$L=0$ で、解 $x(t)=t$ は $\lvert t\rvert\le\delta=\delta/M$ でちょうど $\overline{B}_\delta(0)$ にとどまり、$\lvert t\rvert>\delta$ では球の外に出る。したがって $T>\delta$ としても、球の中の解は $[-\delta,\delta]$ より先へは延びない。thm-pl-main の $T_0=\min\{T,\delta/M\}$ は、「解が球にとどまる」という結論を保つには短くできない。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する