解析学基礎/常微分方程式

常微分方程式とは何ぞや

微分方程式とは、独立変数xと、xの関数y(x)、およびその何階かの導関数を含む方程式である。一般化すれば、微分方程式は

$f(x,y,y',\cdots ,y^{(n)})=0$

の形に書くことのできる方程式である。そして、この方程式に含まれる導関数のうちもっとも高階の導関数が$y^{(n)}$であるとき、これをn階微分方程式と呼び、この方程式を満たすような関数を求める操作を、微分方程式を解く、という。

微分方程式は、大きく分けて常微分方程式と偏微分方程式に分かれる。常微分方程式とは、一変数関数とその導関数からなる方程式である。一方、偏微分方程式とは、多変数関数とその偏導関数との方程式である。ここでは、常微分方程式の解き方について記述することにし、本書では特に断りのない場合「微分方程式」は常微分方程式をさしているものとする。

解の種類

微分方程式は微分された関数が含まれた方程式であるから、その解を求めるためには多くの場合積分操作が必要であり、解には積分定数が含まれる。n階微分方程式であればn個の任意の積分定数が含まれる。このような積分定数を含む形の解を一般解と呼ぶ。

一般解のうち、積分定数にある値を与えた解を特殊解と呼ぶ。

さらに、微分方程式の解のなかには、方程式の解であるにもかかわらず、積分定数にどのような値を代入しても表すことのできない解も存在する。このような解を特異解と呼ぶ。

問題の種類

微分方程式を解くとき、一般解が重要になる場面はそう多くない。現実にはある独立変数xにおける従属変数yの値が定まっていて、その条件を満たすような特殊解が必要になる場合がほとんどである。

たとえば時刻$t=0$において、ある関数$y(t)$の値が$y_{0}$と分かっている時に$y(t)$に関する微分方程式を解くような場合である。

このような、ある初期条件

$y(x_{0})=y_{0}$

を満たすような微分方程式

$f(x,y,y',\cdots ,y^{(n)})$

の特殊解を求める問題を初期値問題といい、これらをみたす特殊解を求めることを初期値問題を解くという。

また、例えば位置$x=0$と$x=L$で常に$y=0$となるような波（固定端）の変位$y(x)$に関する微分方程式を解くという状況もある。

このような、ある境界条件

$y(x_{i})=y_{i}(i=0,1,\cdots ,n)$

を満たすような微分方程式

$f(x,y,y',\cdots ,y^{(n)})$

の特殊解を求める問題を境界値問題といい、これらをみたす特殊解を求めることを境界値問題を解くという。

初等解法

微分方程式を、有限回の式変形や変数変換や積分によって解く方法を初等解法と呼ぶ。はじめに、微分方程式の初等解法について解説する。なお、どのような微分方程式であっても初等解法によって解くことができるとは限らず、この手法が適用できる場合は限られてくることに注意されたい。

まずは、1階の微分方程式について考えることにする。

変数分離形

一般に、n階微分方程式が

$y^{(n)}=f(x,y,y',\cdots ,y^{(n-1)})$

の形で書き表されるとき、これを正規形と呼ぶ。

1階微分方程式の正規形

$y'=f(x,y)$

において、右辺の式が

$f(x,y)=X(x)Y(y)$

のようにxのみの関数とyのみの関数との積の形に変形できるとき、これを変数分離形の微分方程式と呼ぶ。この場合、微分方程式は

$y'=X(x)Y(y)$

の形になっているから、$Y(y)\neq 0$と仮定して両辺を$Y(y)$で割ることにより

${\frac {y'}{Y(y)}}=X(x)$

と変形して、左辺がyとその導関数のみの式、右辺はxのみの式となるように分離することができる。

もし$Y(y)=0$を満たすyの値が存在すれば、その値を$y=a$とすると、もとの微分方程式に代入して

$y'(x)=0$

を得る。一方、いま置いた$y=a$も$y'=0$を満たす関数である。すなわち、微分方程式の解は

$y(x)=a$

と簡単に求めることができる。これは微分方程式の特殊解である。

では、$Y(y)\neq 0$として変数を分離した式に戻ろう。 分離した式の両辺をxで積分して、

$\int {\frac {1}{Y(y(x))}}y'(x)dx=\int X(x)dx$

左辺は置換積分の公式より、$\int {\frac {1}{Y(y(x))}}y'(x)dx=\int {\frac {1}{Y(y)}}dy$であるので、

$\int {\frac {1}{Y(y)}}dy=\int X(x)dx$

を得る。これで両辺の不定積分が計算できれば、微分方程式の解が求まることになる。これは微分方程式の一般解である。

例題

微分方程式$y'=xy$を解く。

これは変数分離形の1階微分方程式である。$y=0$のとき$y'=0$となって、これは微分方程式を満たす。

$y\neq 0$と仮定して両辺をyで割ると、

${\frac {y'}{y}}=x$

であるから、両辺をxで積分して、

$\int {\frac {1}{y}}dy=\int xdx+C$

となる。両辺の不定積分を計算すれば、

$\log |y|={\frac {1}{2}}x^{2}+C$

となるから、これより

$y=e^{{\frac {1}{2}}x^{2}+C}=Ae^{{\frac {1}{2}}x^{2}}$（Aは任意の定数）

とすることができる。これは微分方程式の一般解である。

先に求めた$y=0$は、一般解で$A=0$とした場合であるから、微分方程式の特殊解である。したがって、微分方程式の解は

$y(x)=Ae^{{\frac {1}{2}}x^{2}}$

である。

同次形

一見変数分離形でないように見える微分方程式であっても、適切な変数変換によって変数分離形へ持ち込むことのできる微分方程式が存在する。

1階微分方程式の正規形

$y'=f(x,y)$

において、右辺の式が

$f(x,y)=g\left({\frac {y}{x}}\right)$

のように${\frac {y}{x}}$の関数として記述できるとき、これを同次形の微分方程式と呼ぶ。このとき微分方程式は

$y'=g\left({\frac {y}{x}}\right)$

の形をしている。

$z(x)={\frac {y(x)}{x}}$とおく。このとき$y=xz$であるから、

$y'=(xz)'=z+xz'$

が成り立つ。これを元の微分方程式に代入すると、

$z+xz'=g(z)$

となる。これを$z'$について解くと、

$z'={\frac {g(z)-z}{x}}$

となって、変数分離形の微分方程式となる。

変数分離形の方程式の解き方にしたがってこれを解くと、

${\begin{aligned}{\frac {z'}{g(z)-z}}&={\frac {1}{x}}\\\int {\frac {1}{g(z)-z}}dz&=\int {\frac {1}{x}}dx+C\\\int {\frac {1}{g(z)-z}}dz&=\log |x|+C\end{aligned}}$

となる。これで左辺の不定積分を計算し、$z={\frac {y}{x}}$を代入し直せば微分方程式の解が得られる。

例題

微分方程式$y'={\frac {y}{x}}+{\frac {x}{y}}$を解く。

これは同次形の1階微分方程式である。$z={\frac {y}{x}}$とおくと、$y=xz$であるからこの微分方程式は

${\begin{aligned}z+xz'&=z+{\frac {1}{z}}\\xz'&={\frac {1}{z}}\end{aligned}}$

と書き直すことができる。これは変数分離形の微分方程式である。$z\neq 0$に注意して変数分離を行うと

$zz'={\frac {1}{x}}$

であるから、両辺をxで積分して式変形を行うと、

${\begin{aligned}\int zdz&=\int {\frac {1}{x}}dx\\{\frac {1}{2}}z^{2}&=\log |x|+C\\\log |x|&={\frac {1}{2}}z^{2}+C\\x&=e^{{\frac {1}{2}}z^{2}+C}=Ae^{{\frac {1}{2}}z^{2}}\end{aligned}}$

となる。ここで$z={\frac {y}{x}}$を代入しなおすと、

$x=Ae^{\frac {y^{2}}{2x^{2}}}$

となる。これが求める微分方程式の一般解である。

同次形の応用

正規形の1階微分方程式

$y'=f(x,y)$

について、右辺がxとyの有理関数になっている場合、すなわち

$f(x,y)={\frac {h(x,y)}{g(x,y)}}$

の場合を考える。このとき、$g(x,y)$および$h(x,y)$が特定の形をしている場合は、上手な式変形や変数変換によって同次形の解法を適用することができることが知られている。ここでは、いくつかの例題を用いてそれらの解法を紹介することにする。

例題1

微分方程式$y'={\frac {2x^{2}+3xy+y^{2}}{x^{2}-4xy+2y^{2}}}$を解く。

これは、$g(x,y)$と$h(x,y)$がともにすべての項でx,yについて同次であるような場合である。例えばこのような場合には、右辺の分子と分母を$x^{2}$で割ることで

$y'={\frac {2+3(y/x)+(y/x)^{2}}{1-4(y/x)+2(y/x)^{2}}}$

となって、容易に同次形の微分方程式へ持ち込むことができる。あとは同次形の解法に従って解けばよい。

例題2

微分方程式$y'={\frac {2x+3y-8}{x-y+1}}$を解く。

これは、$g(x,y)$と$h(x,y)$がともにx,yの1次式になっている場合である。例えばこのような場合は、次の手順で解くことができることが知られている。

はじめに、連立方程式

${\begin{cases}2x+3y-8=0\\x-y+1=0\end{cases}}$

を解く。これを解くと、解は$(x,y)=(1,2)$である。この解を用いて、

${\begin{cases}x=u+1\\y=v+2\end{cases}}$

とおく。これをもとの微分方程式へ代入すると、

${\begin{aligned}{\frac {dy}{dx}}&={\frac {2(u+1)+3(v+2)-8}{(u+1)-(v+2)+1}}\\{\frac {dv}{du}}&={\frac {2u+3v}{u-v}}\end{aligned}}$

となる。ここで、

${\frac {dy}{dx}}={\frac {d}{dx}}(v+2)={\frac {dv}{dx}}={\frac {dv}{du}}{\frac {du}{dx}}={\frac {dv}{du}}{\frac {d}{dx}}(x-1)={\frac {dv}{du}}$

を用いた。

このようにx,yからu,vへの変数変換を施すと、例題1で見た形の方程式となり、右辺の分母分子をuで割ることによって同次形の微分方程式として扱うことができる。そしてu,vの式として同次形の微分方程式を解いた後、変数をu,vからx,yに戻せば、求めるべき微分方程式の解が得られる。

例題3

微分方程式$y'={\frac {2x+3y-4}{4x+6y-3}}$を解く。

これは、例題2のようにして連立方程式を解こうとしても、連立方程式の解が存在しないか、あるいは1つに定まらない場合である。このような場合は、右辺の分母をzと置くことによって一般解を求めることができる。

この問題では、

$z=4x+6y-3$

とおくと、分子は

$2x+3y-4={\frac {1}{2}}(4x+6y-3)-{\frac {5}{2}}={\frac {1}{2}}z-{\frac {5}{2}}$

である。また、

$z'=4+6y'$

であるから、この微分方程式は変数分離形へと変形することができて

${\begin{aligned}{\frac {z'-4}{6}}&={\frac {{\frac {1}{2}}z-{\frac {5}{2}}}{z}}\\{\frac {z}{7z-15}}z'&=1\\\int {\frac {z}{7z-15}}dz&=x+C\end{aligned}}$

と変形できる。この左辺に$z=4x+6y-3$を代入すれば求めるべき微分方程式の一般解が求まる。

例題4

微分方程式$y'={\frac {2x+3y-4}{3}}$を解く。

これは、例題3と同様に連立方程式を解こうとしても解が一意に定まらず、かつ、分母が定数になっている場合である。この場合は右辺をzと置けばよい。

$z={\frac {2x+3y-4}{3}}$

とおくと、

$z'={\frac {2+3y'}{3}}={\frac {2+3z}{3}}$

となる。これは変数分離形の微分方程式であるから、その方法に従って解いた後で$z={\frac {2x+3y-4}{3}}$を代入すれば求めるべき一般解が得られる。

例題5

微分方程式$y'={\frac {x^{2}y}{x^{3}+y}}$を解く。

これは、同次形をさらに一般化させた微分方程式である。同次形では、正規形の微分方程式

$y'=f(x,y)$

の右辺$f(x,y)$について、$\lambda$を定数として

$f(\lambda x,\lambda y)=f(x,y)$

が成り立つ。この例題は同次形ではないため、微分方程式の右辺を$f(x,y)$と置いてもこれは成り立たない。しかし、

$f(\lambda x,\lambda ^{3}y)={\frac {(\lambda x)^{2}\lambda ^{3}y}{(\lambda x)^{3}+\lambda ^{3}y}}=\lambda ^{2}{\frac {x^{2}y}{x^{3}+y}}=\lambda ^{2}f(x,y)$

が成り立っている。

一般に、$f(x,y)$について$n\neq 0$として

$f(\lambda x,\lambda ^{n}y)=\lambda ^{n-1}f(x,y)$

が成り立つとき、

$z={\frac {y}{x^{n}}}$

とおいて変数変換を施し、式変形を行うことで変数分離形へ持ち込むことができることが知られている。

この例題では、$z={\frac {y}{x^{3}}}$と置くと、$y=x^{3}z$であるから

$y'=3x^{2}z+x^{3}z'$

これをもとの微分方程式へ代入すると、

${\begin{aligned}3x^{2}z+x^{3}z'&={\frac {x^{2}x^{3}z}{x^{3}+x^{3}z}}\\&={\frac {x^{2}z}{1+z}}\\z'&=-{\frac {3z^{2}+2z}{x(1+z)}}\end{aligned}}$

となる。これは変数分離形の微分方程式であるから、変数分離形の解法に従って解き、最後に$z={\frac {y}{x^{3}}}$を代入すればよい。

1階線型微分方程式

1階微分方程式が線型であるとは、与えられた微分方程式が

$y'+f(x)y=g(x)$

と書けることである。このように書けない1階微分方程式は1階非線型微分方程式という。

斉次1階線型微分方程式

斉次1階線型微分方程式とは、1階線型微分方程式であって、特に$g(x)=0$であるものをいい、この時この微分方程式は斉次であるという。 $g(x)\neq 0$の場合は非斉次であるという。また、斉次は「同次」とも呼ばれることがあるが、本書では前者で統一することにする。

まずは斉次1階線型微分方程式を解いてみよう。 簡単な微分積分法しか知らない我々は、これ程までに限定してやっと解けるようになるのである。

今解こうとしているのは、次の微分方程式である。

$y'+f(x)y=0$

これは変数分離形の微分方程式である。まず$y\neq 0$を仮定して、この式を同値変型する。

${{y'} \over {y}}=-f(x)$

両辺を積分して

${\rm {ln}}|y|=\int {-f(x)dx+C_{0}}$

両辺をeの肩に掛けて、

$|y|=e^{\int {-f(x)dx+C_{0}}}$

右は常に正なので、$e^{C_{0}}=C$として、

$y=Ce^{\int {-f(x)dx}}$

この解法を変数分離法といい、得られた結果がこの斉次方程式の一般解である。

一般解はこのようにして求められたが、$y(x_{0})=y_{0}$となるときの特殊解yを求めなければならないときもある。斉次1階微分方程式の初期値問題について考えてみよう。

初期値問題

$y'+f(x)y=0;y(x_{0})=y_{0}$

を解く。

はじめに微分方程式を解くと、先に導いたように一般解

$y=Ce^{\int {-f(x)dx}}$

を得る。この式の両辺に$(x,y)=(x_{0},y_{0})$を代入すれば、積分定数Cの値が求められるため、改めてそれをこの式に代入しなおすことで特殊解が得られる。

あるいは、微分方程式を解く際に不定積分ではなく$x_{0}$から$x$までの定積分を求めることによって初期値問題を解くこともできる。多少厄介だが、積分記号を外せないときにも解を求めることができる。

変数分離を施した形

${\frac {y'}{y}}=-f(x)$

より、両辺を$x_{0}$から$x$まで定積分する。

${\begin{aligned}&{\int _{x_{0}}^{x}{{y'} \over {y}}\ dx'}={\int _{x_{0}}^{x}-f(x')\ dx'}\\\iff &{\rm {ln}}y-{\rm {ln}}y_{0}={\int _{x_{0}}^{x}-f(x')\ dx'}\\\iff &{{y} \over {y_{0}}}=e^{\int _{x_{0}}^{x}-f(x')\ dx'}\\\iff &y=y_{0}e^{\int _{x_{0}}^{x}-f(x')\ dx'}\\\end{aligned}}$

結局、一般解における積分定数Cが$y_{0}$に、不定積分が定積分になっただけであった。

例題1

微分方程式$y'-4xy=0$を解く。

上の解説の通り、両辺をyで割り変数分離法によって計算する。この微分方程式の一般解は

$y=Ce^{\int {-(-4x)dx}}=Ce^{2x^{2}}$

である。

例題2

次の微分方程式の初期値問題を解け。

$y'+y\sin x=0;y(0)={3 \over 2}$

この微分方程式の一般解として、変数分離法によって

$y=Ce^{\cos x}$

が求められる。この式に$(x,y)=(0,3/2)$を代入すれば、

$C={\frac {3/2}{e^{\cos 0}}}={\frac {3}{2e}}$

したがって求めるべき特殊解は

$y={\frac {3}{2e}}e^{\cos x}={\frac {3}{2}}e^{\cos x-1}$

あるいは、不定積分の代わりに定積分を行うことにより、

$y={3 \over 2}e^{-\int _{0}^{x}\sin t\ dt}={3 \over 2}e^{\cos x-1}$

が導かれる。

非斉次1階線型微分方程式

次に、非斉次1階線型微分方程式

$y'+f(x)y=g(x)$　　　　(1.1)

の解き方を考えてみよう。しかし今、我々にできる事は二つしかない。それは、斉次微分方程式を解くことと、各種式変形を行うことである。これを最大限駆使して解くしかない。具体的には、なんとかして(1.1)を斉次微分方程式

$z'=\nu (x)$　　　　　★

の形に式変形して、これを解くのである。

天下り式であるが、(1.1)にある関数$h(x)$をかけて

$h(x)y'+h(x)f(x)y=h(x)g(x)$　　　　(1.2)

とする。ここで$h(x)$が

$\{h(x)y\}'=h(x)y'+h(x)f(x)y$　　　　(1.3)

をみたすような関数であるとすると、$z=h(x)y,\nu (x)=h(x)g(x)$とすれば★の形に変形できる。

ではそのような$h(x)$は存在するのだろうか。具体的に求めてみる。

$\{h(x)y\}'=h(x)y'+h'(x)y$であるから、これを(1.3)に代入すると

${\begin{aligned}h(x)y'+h'(x)y&=h(x)y'+h(x)f(x)y\\h'(x)&=h(x)f(x)\end{aligned}}$

を得る。$h(x)$についてはこの変数分離形の微分方程式を解けばよい。定数倍は関係ないので、

$h(x)=e^{\int f(x)dx}$　　　　　(1.4)

としてよい。この$h(x)$は積分因子と呼ばれる。

さて、(1.2)と(1.3)より、

$\{h(x)y\}'=h(x)g(x)$　　　　　(1.5)

を得る。これを変形すると、

${\begin{aligned}h(x)y&=\int h(x)g(x)dx+C\\y&={1 \over {h(x)}}\left(\int h(x)g(x)dx+C\right)\end{aligned}}$

あとはこれに(1.4)を代入すると、一般解

$y={\frac {1}{e^{\int f(x)dx}}}\left(\int e^{\int f(x)dx}g(x)dx+C\right)$

を得る。

初期値問題$y'+f(x)y=g(x);y(x_{0})=y_{0}$を解くには、(1.5)の両辺を積分する際に定積分とすれば、

$\int _{x_{0}}^{x}\{h(x)y\}'dx=\int _{x_{0}}^{x}h(x)g(x)dx$

を得る。あとはこれをyについて解けばよい。

以上、非斉次微分方程式の解法を述べた。手順をまとめると、

1. 積分因子$h(x)=e^{\int f(x)dx}$を求める。
2. $g(x)$に積分因子を掛け積分する。
3. それを$h(x)$で割って一般解とする。

となる。

非斉次1階線型微分方程式の別の解法として、定数変化法と呼ばれる方法を紹介する。

非斉次な微分方程式

$y'+f(x)y=g(x)$

を解くのが最終的な目標であるが、ひとまず、右辺を$g(x)=0$とおいて、斉次な微分方程式

$y'+f(x)y=0$

を解くことにする。この形ならば、前々節で見た方法によって、一般解

$y_{h}(x)=Ce^{\int {-f(x)dx}}$

を得ることができる。ここで、非斉次な場合は積分定数のCがxの関数になると考えて、仮に非斉次微分方程式の解を

$y=C(x)y_{h}(x)$

とおく。これを解くべき微分方程式へ代入すると、

${\begin{aligned}\left\{C(x)y_{h}(x)\right\}'+f(x)C(x)y_{h}(x)&=g(x)\\C'y_{h}+Cy_{h}'+f(x)Cy_{h}&=g(x)\\C'y_{h}+C(y_{h}'+f(x)y_{h})&=g(x)\end{aligned}}$

となるが、ここで$y_{h}$が斉次微分方程式$y'+f(x)y=0$の解であることから、

$C'(x)y_{h}(x)=g(x)$

が得られる。この中で未知関数は$C'(x)$のみであるから、両辺を$y_{h}(x)$で割ってxで積分すると、

${\begin{aligned}C(x)&=\int {\frac {g(x)}{y_{h}(x)}}dx+C\\&=\int {\frac {g(x)}{e^{\int {-f(x)dx}}}}dx+C\\&=\int g(x)e^{\int f(x)dx}dx+C\end{aligned}}$

したがって、求めるべき非斉次微分方程式の一般解は、

$y=C(x)y_{h}(x)=e^{\int -f(x)dx}\left\{\int g(x)e^{\int f(x)dx}dx+C\right\}$

となる。これは積分因子を用いて求めた一般解と等しい。

例題1

微分方程式$y'-2xy=x$を解く。

$f(x)=-2x$より、積分因子$h(x)$は、

$h(x)=e^{\int -2xdx}=e^{-x^{2}}$

である。これを与式右辺（$g(x)$）に掛けて積分すると、

$\int h(x)g(x)dx=\int e^{-x^{2}}xdx=-{\frac {1}{2}}e^{-x^{2}}+C$ (1.6)

したがって、微分方程式の一般解は

$y={\frac {-{\frac {1}{2}}e^{-x^{2}}+C}{e^{-x^{2}}}}=-{1 \over 2}+Ce^{x^{2}}$

となる。

あるいは、定数変化法によって求めることもできる。仮に斉次な微分方程式

$y'-2xy=0$

を解くと、この一般解は

$y=Ce^{x^{2}}$

となる。これより、仮に求めるべき微分方程式の解を

$y=C(x)e^{x^{2}}$

と置いて元の微分方程式に代入すると、

$C'(x)e^{x^{2}}=x$

が得られる。これより、

$C(x)=\int {\frac {x}{e^{x^{2}}}}dx=\int e^{-x^{2}}xdx=-{\frac {1}{2}}e^{-x^{2}}+C$

となるから、求める一般解は

$y=e^{x^{2}}\left(-{\frac {1}{2}}e^{-x^{2}}+C\right)=-{\frac {1}{2}}+Ce^{x^{2}}$

である。

例題2

初期値問題$y'-2xy=x;y(1)=2$を解く。

例題1で(1.6)を積分するときに定積分にする。

${\begin{aligned}\int _{1}^{x}\{ye^{-t^{2}}\}'dt&=\int _{1}^{x}e^{-t^{2}}tdt\\\left[ye^{-t^{2}}\right]_{1}^{x}&=\left[-{1 \over 2}e^{-t^{2}}\right]_{1}^{x}\\ye^{-x^{2}}-2e^{-1}&=-{1 \over 2}e^{-x^{2}}+{1 \over 2}e^{-1}\\ye^{-x^{2}}&=-{1 \over 2}e^{-x^{2}}+{5 \over 2}e^{-1}\end{aligned}}$

したがって求める特殊解は

$y(x)=-{1 \over 2}+{5 \over 2}e^{x^{2}-1}$

あるいは、例題1で求めた一般解に$(x,y)=(1,2)$を代入することによってCの値を求めてもよい。

ベルヌーイの微分方程式

1階微分方程式のなかでも、特に

$y'+f(x)y=g(x)y^{n}$

の形の微分方程式をベルヌーイ（Bernoulli）の微分方程式と呼ぶ。$n=0,1$であれば上で見た非斉次1階微分方程式あるいは斉次1階微分方程式の形となり、これらの解法が適用できるが、それ以外の場合でも適切な式変形によって線型微分方程式へ帰着できることが知られている。

ベルヌーイの1階微分方程式

$y'+f(x)y=g(x)y^{n},(n\neq 0,1)$

の両辺に$(1-n)y^{-n}$をかけると、

$(1-n)y^{-n}y'+f(x)(1-n)y^{1-n}=g(x)(1-n)$

となるから、ここで$z=y^{1-n}$とおくと、$z'=(1-n)y^{-n}y'$なので、

$z'+f(x)(1-n)z=g(x)(1-n)$

となる。これはzに関する1階線型微分方程式であるから、定数変化法あるいは積分因子を用いる方法によって計算することができて、一般解

$z=e^{-(1-n)\int f(x)dx}\left\{(1-n)\int g(x)e^{(1-n)\int f(x)dx}dx+C\right\}$

を得る。これに$z=y^{1-n}$を代入しなおすと、

${\begin{aligned}y^{1-n}&=e^{-(1-n)\int f(x)dx}\left\{(1-n)\int g(x)e^{(1-n)\int f(x)dx}dx+C\right\}\\y&=e^{-\int f(x)dx}\left\{(1-n)\int g(x)e^{(1-n)\int f(x)dx}dx+C\right\}^{\frac {1}{1-n}}\end{aligned}}$

を得る。

リッカチの微分方程式

1階微分方程式のなかでも、特に

$y'=f(x)y^{2}+g(x)y+h(x)$

の形に書くことのできる微分法定式をリッカチ（Riccati）の微分方程式と呼ぶ。この形の方程式は初等解法によって一般解を求めることはできない。しかし、なにか1つの特殊解$y_{0}$が見つかれば、それを元にして一般解を求めることができる。

リッカチの微分方程式

$y'=f(x)y^{2}+g(x)y+h(x)$

について、ある特殊解$y_{0}$が与えられているとする。この時、$z=y-y_{0}$とおいて元の微分方程式へ代入すると、

${\begin{aligned}z'+y_{0}'&=f(x)(z+y_{0})^{2}+g(x)(z+y_{0})+h(x)\\z'&=f(x)z^{2}+\left\{2f(x)y_{0}+g(x)\right\}z+\left\{f(x)y_{0}^{2}+g(x)y_{0}+h(x)-y_{0}'\right\}\end{aligned}}$

となる。ここで$y_{0}$がこの微分方程式の特殊解であることから

$y_{0}'=f(x)y_{0}^{2}+g(x)y_{0}+g(x)$

が成り立っているので、

$z'=f(x)z^{2}+\left\{2f(x)y_{0}+g(x)\right\}z$

となる。これはベルヌーイの微分方程式で$n=2$の場合であるから、前節で見た方法で解くことができる。両辺に$-z^{-2}$をかけて

$-z^{-2}z'=-f(x)-\left\{2f(x)y_{0}+g(x)\right\}z^{-1}$

さらに$u=z^{-1}$とおくと$u'=-z^{-2}z'$であるから

$u'=-f(x)-\left\{2f(x)y_{0}+g(x)\right\}u$

となって、1階線型微分方程式に帰着する。この一般解は、前節で見た式から

$z=e^{\int \left\{2f(x)y_{0}+g(x)\right\}dx}\left\{-\int f(x)e^{\int \left\{2f(x)y_{0}+g(x)\right\}dx}dx+C\right\}^{-1}$

となり、求めるべき微分方程式の一般解も

${\begin{aligned}y=z+y_{0}&=e^{\int \left\{2f(x)y_{0}+g(x)\right\}dx}\left\{-\int f(x)e^{\int \left\{2f(x)y_{0}+g(x)\right\}dx}dx+C\right\}^{-1}+y_{0}\\&={\frac {e^{\int \left\{2f(x)y_{0}+g(x)\right\}dx}}{-\int f(x)e^{\int \left\{2f(x)y_{0}+g(x)\right\}dx}dx+C}}+y_{0}\end{aligned}}$

と求まる。

演習

次の方程式を解け

1. $f'(x)+f(x)\cos x=0$
2. $f'(x)+f(x){\sqrt {x}}\sin x=0$
3. $f'(x)+f(x){{2x} \over {x^{2}+1}}={1 \over {x^{2}+1}}$
4. $f'(x)+f(x)=xe^{x}$
5. $f'(x)+f(x)x^{2}=1$
6. $f'(x)+f(x){x \over {x^{2}+1}}=1-{x^{3} \over {x^{4}+1}}f(x)$
7. $f'(x)+f(x){\sqrt {x^{2}+1}}=0,f(0)={\sqrt {5}}$
8. $f'(x)+f(x){\sqrt {x^{2}+1}}e^{-}x=0,f(0)=1$
9. $f'(x)+f(x){\sqrt {x^{2}+1}}e^{-}x=0,f(0)=0$
10. $f'(x)-2xf(x)=x,f(0)=1$
11. $f'(x)+xf(x)=x+1,f({3 \over 2})=0$
12. $f'(x)+2f(x)={1 \over {x^{2}+1}},f(1)=2$
13. $f'(x)-2xf(x)=x,f(0)=1$
14. $(x^{2}+1)f'(x)+xf(x)=(x^{2}+1)^{5 \over 2}$
15. $(x^{2}+1)f'(x)+3xf(x)=(x^{2}+1)^{5 \over 2},f(1)={1 \over 3}$

原子核の崩壊速度

線型微分方程式のひとつの応用例として、原子核の崩壊に関するものを見てみよう。

物理学者ラザフォードは、放射性元素の原子核は不安定で、一定の割合で崩壊する事を示した。つまり、原子核の数をyという関数で表すことにすれば

y'=-λy　　　　　　　　　　　　　　　　　　　(5.1)

という関係式が成り立つ。ここで比例定数λは崩壊定数と呼ばれる正数である。 この関係式は、まさに一階線形常微分方程式となっているので、これまでに述べた方法で解くことができる。 y(x₀)=y₀とすれば、(5.1)は

$y=y_{0}e^{-\lambda (x-x_{0})}$　　　　　　　　　　　　　　　　　　　(5.2)

と解ける。適当に文字を置き換えると、高等学校理科 物理II 原子と原子核の1.2.3で述べた式が導かれたことになる。

一階定数係数連立線型常微分方程式と高階定数係数線型常微分方程式

連立線型常微分方程式と行列の指数関数

上の節では一階の線型常微分方程式の解法を述べた。その中でも最もやさしい定数係数の方程式

$y'=ay$

の解は、変数分離法により簡単に求まり、

$y=Ce^{xa}$

であった。ただし、C=y(0)である。

次に、n本の一階定数係数線型常微分方程式を連立させた方程式

${\begin{cases}y_{1}'=a_{11}y_{1}+a_{12}y_{2}+\cdots a_{1n}y_{n}\\y_{2}'=a_{21}y_{1}+a_{22}y_{2}+\cdots a_{2n}y_{n}\\\vdots \\y_{n}'=a_{n1}y_{1}+a_{n2}y_{2}+\cdots a_{nn}y_{n}\end{cases}}$

を考えよう。この方程式は、行列を用いて

$\mathbf {y} '=A\mathbf {y} \cdots \bigstar$

と表すことができる。ただし $\mathbf {y} ={\begin{pmatrix}y_{1}\\y_{2}\\\vdots \\y_{n}\end{pmatrix}},A={\begin{pmatrix}a_{11}&a_{12}&\cdots &a_{1n}\\a_{21}&a_{22}&\cdots &a_{2n}\\\vdots &&&\vdots \\a_{n1}&a_{n2}&\cdots &a_{nn}\\\end{pmatrix}}$である。

方程式が1本のときの例から類推すれば、この連立方程式の解は

$e^{xA}$

のようなものが定義できれば、それを用いて表せそうである。しかし、行列の指数関数をどうやって定義すればよいだろうか？そのために、そもそも実数上の関数としての指数関数がどのように定義されるかを考えてみると、次のようにしてTaylor展開で定義できることが思い出される。

$e^{x}=\sum _{k=0}^{\infty }{\frac {x^{k}}{k!}}$

行列であっても、この式に代入することは可能そうである。すなわち、次のように定義する。

定義　正方行列Aに対して、$e^{xA}:=\sum _{k=0}^{\infty }{\frac {(xA)^{k}}{k!}}$

この級数が収束するのか、またどの程度よい収束をするのかが問題だが、結論から言えば一様絶対収束する。詳しい証明は省くが、ゆえにこの級数を項別微分することができ、

$(e^{xA})'=Ae^{xA}$

が成り立つ。

このことから、連立線型微分方程式は初期条件を与えると次のように解けることがわかる。

定理

$\mathbf {y} =e^{xA}\mathbf {y} (0)$

は、方程式$\bigstar$の初期値$\mathbf {y} (0)={\begin{pmatrix}c_{1}\\c_{2}\\\cdots \\c_{n}\end{pmatrix}}$における解になっている。

実際に解になっていることは代入によって確かめることができる。

高階定数係数線型常微分方程式

次に、n階の定数係数線型常微分方程式

$y^{(n)}=a_{n-1}y^{(n-1)}+\cdots +a_{1}y'+a_{0}y$

を考える。この方程式は、実は次のようにして連立常微分方程式とみなして行列を使って表せる。

${\frac {d}{dx}}{\begin{pmatrix}y\\y'\\\vdots \\y^{(n-2)}\\y^{(n-1)}\end{pmatrix}}={\begin{pmatrix}0&1&0&\cdots &0\\0&0&1&\cdots &0\\\vdots &&&&\vdots \\0&0&0&\cdots &1\\a_{0}&a_{1}&a_{2}&\cdots &a_{n-1}\end{pmatrix}}{\begin{pmatrix}y\\y'\\\vdots \\y^{(n-2)}\\y^{(n-1)}\end{pmatrix}}$

よって、上の節で述べた方法により初期値問題を解くことができる。

具体的な行列に対する計算法

では、具体的な係数行列が与えられたとき、どのようにすれば行列の指数関数が計算できるかを見てみよう。

対角行列の場合

対角行列

$D={\begin{pmatrix}c_{1}&0&\cdots &0\\0&c_{2}&\cdots &0\\\vdots &&&\vdots \\0&0&\cdots &c_{n}\end{pmatrix}}$

に対して$e^{xD}$を計算してみよう。

すぐにわかるように、

$D^{k}={\begin{pmatrix}c_{1}^{k}&0&\cdots &0\\0&c_{2}^{k}&\cdots &0\\\vdots &&&\vdots \\0&0&\cdots &c_{n}^{k}\end{pmatrix}}$

である。よって、各成分ごとの計算から $e^{xD}={\begin{pmatrix}e^{c_{1}x}&0&\cdots &0\\0&e^{c_{2}x}&\cdots &0\\\vdots &&&\vdots \\0&0&\cdots &e^{c_{n}x}\end{pmatrix}}$ である。

対角化可能な行列の場合

行列Aが$P^{-1}AP=D$と対角化可能な場合も行列の指数関数は容易に計算できる。なぜならば、

$A^{k}=(PDP^{-1})^{k}=PD^{k}P^{-1}$

なので、これを代入することで

$e^{xA}=Pe^{xD}P^{-1}$

となり、対角行列の指数関数は容易に計算できるからである。

対角化不可能な行列の場合

係数行列が対角化不可能なときは上記のようにはいかず、一般にはJordan標準形を用いることになる。しかし、特殊な場合にはそこまでの計算をする必要はない。たとえば、固有値がすべて等しい場合には次のようにして計算することができる。

n次正方行列Aのn個の固有値がすべて$\lambda$のとき、この行列の固有多項式は$(t-\lambda )^{n}$なので、Cayley-Hamiltonの定理より

$(A-\lambda I)^{n}=O$

である。このことを用いると、

${\begin{aligned}e^{xA}&=e^{\lambda xI+x(A-\lambda I)}\\&=e^{\lambda x}\sum _{k=0}^{\infty }{\frac {x^{k}(A-\lambda I)^{k}}{k!}}\\&=e^{\lambda x}\sum _{k=0}^{n-1}{\frac {x^{k}(A-\lambda I)^{k}}{k!}}\\\end{aligned}}$

と有限回の計算で指数関数を計算することができる。

具体例

二階の線型常微分方程式の具体例として、ばねにつながれた物体の運動を記述してみよう。ばねにつながれた物体の時刻xにおける変位をyとする。このとき、ばねから物体が受ける力は（負の比例定数で）変位に比例することが知られている。このことを用いて物体の運動方程式を記述すると、

$y''={\frac {-k}{m}}y$

となる。ただしkはばね定数と呼ばれる正の数、mは物体の質量である。

この方程式を行列を用いて書き直すと、

${\frac {d}{dx}}{\begin{pmatrix}y\\y'\end{pmatrix}}={\begin{pmatrix}0&1\\{\frac {-k}{m}}&0\end{pmatrix}}{\begin{pmatrix}y\\y'\end{pmatrix}}$

と表せる。$A={\begin{pmatrix}0&1\\{\frac {-k}{m}}&0\end{pmatrix}}$とする。この行列は対角化できるので、指数関数が計算できて、

$e^{xA}={\begin{pmatrix}\cosh(i{\sqrt {\frac {k}{m}}}x)&-i{\sqrt {\frac {m}{k}}}\sinh(i{\sqrt {\frac {k}{m}}}x)\\i{\sqrt {\frac {k}{m}}}\sinh(i{\sqrt {\frac {k}{m}}}x)&\cosh(i{\sqrt {\frac {k}{m}}}x)\\\end{pmatrix}}={\begin{pmatrix}\cos({\sqrt {\frac {k}{m}}}x)&{\sqrt {\frac {m}{k}}}\sin({\sqrt {\frac {k}{m}}}x)\\-{\sqrt {\frac {k}{m}}}\sin({\sqrt {\frac {k}{m}}}x)&\cos({\sqrt {\frac {k}{m}}}x)\\\end{pmatrix}}$

である。初期条件を ${\begin{pmatrix}y(0)\\y'(0)\end{pmatrix}}={\begin{pmatrix}y_{0}\\v_{0}\end{pmatrix}}$ で定めると、解は

${\begin{pmatrix}y\\y'\end{pmatrix}}={\begin{pmatrix}\cos({\sqrt {\frac {k}{m}}}x)&{\sqrt {\frac {m}{k}}}\sin({\sqrt {\frac {k}{m}}}x)\\-{\sqrt {\frac {k}{m}}}\sin({\sqrt {\frac {k}{m}}}x)&\cos({\sqrt {\frac {k}{m}}}x)\\\end{pmatrix}}{\begin{pmatrix}y_{0}\\v_{0}\end{pmatrix}}={\begin{pmatrix}y_{0}\cos({\sqrt {\frac {k}{m}}}x)+v_{0}{\sqrt {\frac {m}{k}}}\sin({\sqrt {\frac {k}{m}}}x)\\-y_{0}{\sqrt {\frac {k}{m}}}\sin({\sqrt {\frac {k}{m}}}x)+v_{0}\cos({\sqrt {\frac {k}{m}}}x)\end{pmatrix}}$

と求められた。これがばねによって振動する物体の時刻xにおける変位と速度である。

常微分方程式の初期値問題の解の存在と一意性

前節まででみたように、いくつかの微分方程式については積分計算によって解を具体的に求めることができるが、一方でそのような求積法の存在しない常微分方程式も多い。だが、そのような方程式についても、ある条件を満たせば解の存在や一意性が保証されることがある。ここではそのような例を見ていこう。

もし解の存在や一意性が保証されるならば、簡単に求積できない微分方程式でも少しは調べやすくなる。一意性が保証されるということは、まぐれやあてずっぽうであっても解をひとつみつけさえすれば、解けたのと同じになるからだ。また、ここで扱う存在と一意性に関する定理は、その解を（ある関数列の極限として）具体的に構成する方法を含んでおり、その意味であてずっぽうではなく解を見つける方法を提供してくれてもいるのである。

本節では、独立変数xの関数yについての1階常微分方程式

$y'=f(x,y),\ y(x_{0})=y_{0}$…(*)

について考える。関数yが(*)を満たすことは、

$y(x)=y_{0}+\int _{x_{0}}^{x}f(t,y(t))dt$…(*)'

を満たすことと同値であることも注意しておく。2変数関数$f(x,y)$に対していくつかの仮定を課したときに、この方程式の解がどのように構成されるかを見ていく。

冪級数による解法

本節では、fが次の仮定(H1)を満たすとする。

(H1) $f(x,y)$は点$(x_{0},y_{0})$の近傍で解析的（すなわち冪級数展開可能）であり、$f(x,y)=\sum _{j,l=0}^{\infty }f_{j,l}(x-x_{0})^{j}(y-y_{0})^{l}$と表される。

このとき、次が成り立つ。

定理5.1.1 $f(x,y)$が仮定(H1)を満たすとき、(*)を満たす$x=x_{0}$の近傍で解析的な関数yがただひとつ存在する。

これを証明したい。ただ、冪級数の中心が一般の形だと計算が煩雑になるので、ここでは次の形の定理を証明することにする。

定理5.1.1' $f(x,y)$が原点の近傍で解析的であり、$f(x,y)=\sum _{j,l=0}^{\infty }f_{j,l}x^{j}y^{l}$と表されるとき、常微分方程式

$y'=f(x,y),\ y(0)=0$…(☆)

を満たす$x=0$の近傍で解析的な関数yがただひとつ存在する。

いくつかの補題に分けて証明しよう。

補題5.1.2 冪級数$y=\sum _{j=0}^{\infty }y_{j}x^{j}$であって(☆)を満たすものがあるならば、その係数$y_{j}$は一意に定まる。

（証明） $y_{0}=0$である。$j\geq 1$のときは$y=\sum _{j=0}^{\infty }y_{j}x^{j}$を(☆)に代入すると、

$(lhs)=y_{1}+2y_{2}x+3y_{3}x^{2}+\cdots$

${\begin{aligned}(rhs)&=f_{0,0}+f_{0,1}(y_{1}x+y_{2}x^{2}+\cdots )+f_{0,2}(y_{1}x+y_{2}x^{2}+\cdots )^{2}+f_{1,0}x+f_{1,1}x(y_{1}x+y_{2}x^{2}+\cdots )+f_{2,0}x^{2}+\cdots \\&=f_{0,0}+(f_{0,1}y_{1}+f_{1,0})x+(f_{0,1}y_{2}+f_{0,2}y_{1}^{2}+f_{1,1}y_{1}+f_{2,0})x^{2}+\cdots \\\end{aligned}}$

であり、次数の低い方から係数を比較することで、係数$y_{j}$が順に決まっていくことがわかる。//

補題5.1.3 $f(x,y)$の冪級数展開の優級数$F(x,y)$があるとき、常微分方程式

$y'=F(x,y),\ y(0)=0$…(☆)'

の冪級数解は、補題5.1.2で定まる(☆)の冪級数解の優級数である。

（証明） $f(x,y)=\sum _{j,l=0}^{\infty }f_{j,l}x^{j}y^{l},F(x,y)=\sum _{j,l=0}^{\infty }F_{j,l}x^{j}y^{l}$とし、(☆)の解を$y=\sum _{j=0}^{\infty }y_{j}x^{j}$、(☆)'の解を$Y=\sum _{j=0}^{\infty }Y_{j}x^{j}$とする。 ただし、$y,Y$の冪級数表示は現時点では収束性については何も分かっていない、形式的冪級数である。すべてのj,lについて$|f_{j,l}|\leq |F_{j,l}|$が成り立つならばすべてのjについて$|y_{j}|\leq |Y_{j}|$であることを数学的帰納法で証明する。$y_{0}=Y_{0}=0,y_{1}=f_{0,0},Y_{1}=F_{0,0}$なので、$j=0,1$のときは成り立つ。 $j\leq m$なるすべてのjで成り立つと仮定する。補題5.1.2の証明から、$y_{m}$は$f_{j,l}(j+l\leq m-1),y_{j}(j\leq m-1)$に関する多項式の値であり、その係数は非負である。$Y_{m}$も同様に、同じ非負係数多項式に$F_{j,l}(j+l\leq m-1),Y_{j}(j\leq m-1)$を代入した値である。よって、帰納法の仮定より、$|y_{m}|\leq |Y_{m}|$が成り立つ。よって、すべての自然数jについて$|y_{j}|\leq |Y_{j}|$が成り立つ。//

補題5.1.4 $f(x,y)$が原点の近傍$|x|\leq r,|y|\leq \rho$において$|f(x,y)|\leq M$を満たすとき、$\sum _{j,l=0}^{\infty }{\frac {M}{r^{j}\rho ^{l}}}x^{j}y^{l}$は$f(x,y)$の冪級数展開の優級数である。

（証明） $|f_{j,l}|\leq {\frac {M}{r^{j}\rho ^{l}}}$を示せばよい。$f(x,y)$の定義域を複素変数に拡張してw:コーシーの積分公式を用いると、$|x|<r,|y|<\rho$のとき $f(x,y)=\int _{|\zeta |=r}{\frac {d\zeta }{2i\pi }}\int _{|\xi |=\rho }{\frac {d\xi }{2i\pi }}{\frac {f(\zeta ,\xi )}{(\zeta -x)(\xi -y)}}=-{\frac {1}{4\pi ^{2}}}\int _{|\zeta |=r}d\zeta \int _{|\xi =\rho }d\xi f(\zeta ,\xi )\left(\sum _{j=0}^{\infty }{\frac {x^{j}}{\zeta ^{j+1}}}\right)\left(\sum _{l=0}^{\infty }{\frac {y^{l}}{\xi ^{l+1}}}\right)=-\sum _{j,l=0}^{\infty }{\frac {1}{4\pi ^{2}}}\int _{|\zeta |=r}d\zeta \int _{|\xi |=\rho }d\xi {\frac {f(\zeta ,\xi )}{\zeta ^{j+1}\xi ^{l+1}}}x^{j}y^{l}$であるから、$|f_{j,l}|\leq {\frac {1}{4\pi ^{2}}}\int _{|\zeta |=r}d\zeta \int _{|\xi |=\rho }d\xi {\frac {|f(\zeta ,\xi )|}{|\zeta |^{j+1}|\xi |^{l+1}}}\leq {\frac {M}{r^{j}\rho ^{l}}}$である。//

補題5.1.5 補題5.1.3の微分方程式で$F(x,y)=\sum _{j,l=0}^{\infty }{\frac {M}{r^{j}\rho ^{l}}}x^{j}y^{l}$としたものの解は、$x=0$の近傍で解析的な関数であり、収束する冪級数で表される。

（証明）

$y'=\sum _{j,l=0}^{\infty }{\frac {M}{r^{j}\rho ^{l}}}x^{j}y^{l}={\frac {M}{\left(1-{\frac {x}{r}}\right)\left(1-{\frac {y}{\rho }}\right)}}$

は変数分離形なので解を求めることができて、

$\left(1-{\frac {y}{\rho }}\right)dy={\frac {M}{\left(1-{\frac {x}{r}}\right)}}dx$

$-{\frac {\rho }{2}}\left(1-{\frac {y}{\rho }}\right)^{2}=-Mr\log \left(1-{\frac {x}{r}}\right)+C$

であり、$y(0)=0$より$C=-{\frac {\rho }{2}}$であることに注意して整理すると、

$y=\rho -{\sqrt {\rho \left(\rho +2Mr\log \left(1-{\frac {x}{r}}\right)\right)}}$

である。これは確かに$|x|<r(1-e^{-{\frac {\rho }{2Mr}}})$で解析的な関数である。//

（定理5.1.1’の証明） 補題5.1.3,5.1.4,5.1.5より、補題5.1.2の冪級数は収束する優級数をもち、したがって自身も収束する。よって、この冪級数の極限として、解が一意的に存在することがわかる。//

ピカールの逐次近似法

本節では、fが次の仮定(H2)を満たすとする。

(H2) $f(x,y)$は点$(x_{0},y_{0})$の近傍$D=\{|x-x_{0}|\leq r,|y-y_{0}|\leq \rho \}$でリプシッツ連続である、すなわちあるKに対して$|f(x_{1},y_{1})-f(x_{2},y_{2})|\leq K{\sqrt {(x_{1}-x_{2})^{2}+(y_{1}-y_{2})^{2}}}$が成り立つ。

このとき、解は次のようにして構成される。

定義5.2.1 $f(x,y)$が仮定(H2)を満たすとき、漸化式$y_{j+1}(x)=y_{0}+\int _{x_{0}}^{x}f(t,y_{j}(t))dx,\ y_{0}(x)=y_{0}$で定まる関数列$y_{j}$をピカールの逐次近似列という。

定理5.2.2 $f(x,y)$が仮定(H2)を満たすとき、$M=\max _{(x,y)\in D}|f(x,y)|,\delta =\min \left\{{\frac {\rho }{M}},r\right\}$とする。閉区間$[x_{0}-\delta ,x_{0}+\delta ]$において(*)を満たす関数yがただひとつ存在し、それはピカールの逐次近似列$y_{j}$の$j\to \infty$における極限として定義される。

これをいくつかの補題に分けて証明しよう。

補題5.2.3 $x_{0}-\delta \leq x\leq x_{0}+\delta$のとき、$\left|\int _{x_{0}}^{x}f(t,y(t))dt\right|\leq M|x-x_{0}|\leq \rho$である。

（証明） $\left|\int _{x_{0}}^{x}f(t,y(t))dt\right|\leq M|x-x_{0}|\leq M\delta \leq \rho$//

補題5.2.3を帰納的に用いることで、任意のjについて$y_{j}$の値域が$|y-y_{0}|\leq \rho$に含まれ、したがって関数列$y_{j}$がwell-definedであることが従う。

次に、解の一意性を先に示しておこう。

補題5.2.4 $f(x,y)$が仮定(H2)を満たすとき、閉区間$[x_{0}-\delta ,x_{0}+\delta ]$において(*)を満たす関数はただ一つである。

（証明） $y(x),{\tilde {y}}(x)$がともに(*)'を満たすとすると、fがリプシッツ連続であることから

$|y(x)-{\tilde {y}}(x)|=\left|\int _{x_{0}}^{x}f(t,y(t))-f(t,{\tilde {y}}(t))dt\right|\leq \left|\int _{x_{0}}^{x}|f(t,y(t))-f(t,{\tilde {y}}(t))|dt\right|\leq K\left|\int _{x_{0}}^{x}|y(t)-{\tilde {y}}(t)|dt\right|$

である。一方、補題5.2.3より、

$|y(x)-{\tilde {y}}(x)|\leq |y(x)-y_{0}|+|{\tilde {y}}(x)-y_{0}|=\left|\int _{x_{0}}^{x}f(t,y(t))dt\right|+\left|\int _{x_{0}}^{x}f(t,{\tilde {y}}(t))dt\right|\leq 2\rho$

なので、

$|y(x)-{\tilde {y}}(x)|\leq 2\rho K|x-x_{0}|$

であり、よって

$|y(x)-{\tilde {y}}(x)|\leq \left|\int _{x_{0}}^{x}2\rho K|t-x_{0}|dt\right|=2\rho {\frac {(K|x-x_{0}|)^{2}}{2}}$

である。同様に繰り返すことで、任意の自然数lに対して

$|y(x)-{\tilde {y}}(x)|\leq 2\rho {\frac {(K|x-x_{0}|)^{l}}{l!}}\leq 2\rho {\frac {(K\delta )^{l}}{l!}}$

であることがわかるが、$\lim _{l\to \infty }{\frac {(K\delta )^{l}}{l!}}=0$なので、$y(x)={\tilde {y}}(x)$である。//

補題5.2.6 関数列$y_{j}$は一様収束する。

（証明） 補題5.2.3とfがリプシッツ連続であることより、

$|y_{1}(x)-y_{0}|\leq M|x-x_{0}|$

$|y_{2}(x)-y_{1}(x)|\leq \left|\int _{x_{0}}^{x}|f(t,y_{1}(t))-f(t,y_{0}(t))|dt\right|\leq K\left|\int _{x_{0}}^{x}|y_{1}(t)-y_{0}(t)|dt\right|\leq K\int _{x_{0}}^{x}M|t-x_{0}|dt\leq KM{\frac {|x-x_{0}|^{2}}{2}}$

以下同様に繰り返して、

$|y_{j}(x)-y_{j-1}(x)|\leq {\frac {M}{K}}{\frac {(K|x-x_{0}|)^{l}}{l!}}\leq {\frac {M}{K}}{\frac {(K\delta )^{l}}{l!}}$

である。よって、$j<l$のとき

$|y_{j}(x)-y_{l}(x)|\leq \sum _{i=j+1}^{l}{\frac {M}{K}}{\frac {(K\delta )^{i}}{i!}}$

であるが、右辺は$j\to \infty$において0に収束するので、$y_{j}$は一様収束する。//

以上で定理5.2.2を示す準備は整った。

（定理5.2.2の証明） $y_{j}$が一様収束することに注意して$y_{j+1}(x)=y_{0}+\int _{x_{0}}^{x}f(t,y_{j}(t))dx$の両辺の$j\to \infty$における極限を考えると、

$\lim _{j\to \infty }y_{j}(x)=y_{0}+\int _{x_{0}}^{x}f(t,\lim _{j\to \infty }y_{j}(t))dx$

である。これは$\lim _{j\to \infty }y_{j}(x)$が(*)'の解であることを示している。

コーシーの折れ線とペアノの定理

本節では、fが次の仮定(H3)を満たすとする。

(H3) $f(x,y)$は点$(x_{0},y_{0})$の近傍$D=\{|x-x_{0}|\leq r,|y-y_{0}|\leq \rho \}$で連続である。

このとき、解は次のようにして構成される。

定義5.3.1 自然数jに対し、$x_{0}-r\leq x\leq x_{0}+r$における関数$y_{j}$を次のように定める。まず、$y_{j}(x_{0})=y_{0}$とする。次に、非負整数iに対して$x_{i}=x_{0}+{\frac {ir}{j}}$と定め、$x_{i}<x\leq x_{i+1}$のときには

$y_{j}(x)=y_{0}+{\frac {r}{j}}\sum _{d=0}^{i-1}f(x_{d},y_{j}(x_{d}))+(x-x_{i})f(x_{i},y_{j}(x_{i}))$…(#)

とする。$x<0$についても同様にする。このようにして定まる関数$y_{j}$のグラフは連続な折れ線になり、これをコーシーの折れ線という。

前節までに見た「解析的」や「リプシッツ連続」と比べ、「連続」はとても弱い仮定であり、より多くの関数が満たしている。だがその分本節では解の一意性は失われ、存在しか示すことができない。すなわち、次が成り立つのみである。

定理5.3.2 コーシーの折れ線$y_{j}$は一様収束する部分列$y_{j_{l}}$を持ち、$\lim _{l\to \infty }y_{j_{l}}$は方程式(*)の解である。

まず、次の補題を示す。

補題5.3.3 $M=\max _{(x,y)\in D}|f(x,y)|,\delta =\min \left\{{\frac {\rho }{M}},r\right\}$とする。関数列$y_{j}$は$x_{0}-\delta \leq x\leq x_{0}+\delta$において一様有界かつ同程度連続である。

（証明） 定義より

$|y_{j}(x_{1})-y_{j}(x_{2})|\leq M|x_{1}-x_{2}|$

なので、同程度連続である。また、この式で$x_{2}=0$とすると

$|y_{j}(x)|\leq |y_{0}|+M\delta \leq |y_{0}|+\rho$

なので、一様有界である。//

補題5.3.3とアスコリ＝アルツェラの定理により、$y_{j}$が一様収束する部分列を持つことがわかるので、あとはこの部分列の極限が解になっていることを示せばよい。

（定理5.3.2の証明） $y(x)=\lim _{l\to \infty }y_{j_{l}}$が(*)'の解になっていることを示したい。(*)と(#)を辺々引いた式

$y(x)-y_{j,l}(x)+\sum _{d=0}^{i-1}\int _{x_{d}}^{x_{d+1}}(f(x,y(x))-f(x_{d},y_{j}(x_{d})))dx+\int _{t_{i}}^{t}(f(x,y(x))-f(x_{i},y_{j}(x_{i})))dx=0$

が成り立つことを示せばよい。左辺を$A_{j_{l}}$とする。$A_{j_{l}}=0$を示したい。まず、任意の$\varepsilon >0$に対してあるNが存在して、$l>N$ならば$|y(x)-y_{j_{l}}(x)|<\varepsilon$である。次に有界閉集合D上の連続関数fは一様連続なので、任意の$\varepsilon '>0$に対して、$|x_{1}-x_{2}|+|y_{1}-y_{2}|<\delta '$ならば$|f(x_{1},y_{1})-f(x_{2},y_{2})|<\varepsilon '$となるように、$\delta '>0$をとることができる。この$\delta$に対して$\varepsilon <{\frac {\delta '}{2}}$を満たすように$\varepsilon >0$をとり、この$\varepsilon$に対して$|y(x)-y_{j_{l}}(x)|<\varepsilon$かつ$(M+1){\frac {r}{j_{l}}}<{\frac {\delta '}{2}}$を満たすようにlをとれば、$x_{d}<x\leq x_{d+1}$のときには

$|x-x_{d}|+|y(x)-y_{j_{l}}(x_{d})|\leq |x-x_{d}|+|y(x)-y_{j_{l}}(x)|+|y_{j_{l}}(x)-y_{j_{l}}(x_{d})|\leq {\frac {r}{j_{l}}}+\varepsilon +M{\frac {r}{j_{l}}}<\delta '$

なので、$|f(x,y(x))-f(x_{d},y_{j}(x_{d}))|<\varepsilon '$である。よって、

$|A_{j_{l}}|<\varepsilon +(i+1)\int _{x_{d}}^{x_{d+1}}\varepsilon '=\varepsilon +{\frac {(i+1)r\varepsilon '}{j_{l}}}$

である。$\varepsilon ,\varepsilon '$は任意なので、$A_{j_{l}}=0$である。

陰関数型の1階常微分方程式

陰関数型の1階常微分方程式

$\displaystyle F{\Bigl (}k+\ell x+my,\;\;{\frac {dy}{dx}}\;{\Bigr )}=0$

は求積法で一般解を表示することができる。ここに、$F$ は任意の既知関数であり、$k,\;\ell ,\;m$ は任意定数である。 この陰関数型1階常微分方程式の一般解は、次に示す三通りの式で与えられる。

$x=\!\int \!\left\{{\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;d\phi (t)\;}{dt}}\;\right\}dt+C,$

$k+\ell x+my=\phi (t),$

$F{\bigl (}\phi (t),\;\psi (t){\bigr )}\equiv 0.$

ここに、$t$ は媒介変数であり、$\phi (t)$ と $\psi (t)$ は $t$ の関数で、 $F{\bigl (}\phi (t),\;\psi (t){\bigr )}\equiv 0$ は $\phi (t),\;\psi (t)$ に関する恒等式である。なお $C$ は積分定数である。 以下で、その解法を示す。

与えられた常微分方程式 $\displaystyle F{\Bigl (}k+\ell x+my,\;\;{\frac {dy}{dx}}\;{\Bigr )}=0$ に対して、$t$ を媒介変数とする任意関数 $\phi (t),\;\psi (t)$ を導入し、

$k+\ell x+my=\phi (t),$

$\displaystyle {\frac {dy}{dx}}=\psi (t)$

と置く。ただし、$m\neq 0$ とする。 上式 $k+\ell x+my=\phi (t)$ の両辺を$x$で微分すると、

$\displaystyle \ell +m{\frac {dy}{dx}}={\frac {d\phi (t)}{dt}}\cdot {\frac {dt}{dx}}$

となる。 ここで、${\frac {dy}{dx}}=\psi (t)$ と $\ell +m{\frac {dy}{dx}}={\frac {d\phi (t)}{dt}}\cdot {\frac {dt}{dx}}$ から ${\frac {dy}{dx}}$ を消去すると、

$\displaystyle \ell +m\psi (t)={\frac {d\phi (t)}{dt}}\cdot {\frac {dt}{dx}}$

を得る。この式を変形すると、

${\frac {dx}{dt}}={\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {d\phi (t)}{dt}}$

となる。上式は変数分離形であるから積分すると、

$x=\!\int \!\left\{{\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}\,\right\}dt+C,$ ・・・・・（▲）

となり、$x$ が $t$ の関数として表示された。 これを用いれば $y$ は $k+\ell x+my=\phi (t)$ と上式 $x=\int \left\{{\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {d\phi (t)}{dt}}\right\}dt+C$ により $t$ の関数として与えられる。なお $C$ は積分定数である。

途中の計算は省略して $y$ の式のみを以下に示しておく。

$\displaystyle y={\frac {1}{\;m\;}}\left[\phi (t)-k-\ell \left\{\int \left(\!{\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}\,\right)dt+C\right\}\right].$ ・・・・・（▲▲）

上記の（▲）式と（▲▲）式を用いて、

$\displaystyle {\frac {dy}{\;dx\;}}={\frac {\displaystyle {\frac {dy}{\;dt\;}}}{\displaystyle {\frac {dx}{\;dt\;}}}}$

を計算すると、

$\displaystyle {\frac {dy}{\;dx\;}}=\psi (t)$

を得る。この $\displaystyle {\frac {dy}{\;dx\;}}=\psi (t)$ は、 $\displaystyle F{\Bigl (}k+\ell x+my,\;\;{\frac {dy}{dx}}\;{\Bigr )}=0$ を解く時、 始めに仮定した式と同一である。 次に、${\frac {\;{d}x\;}{dt}}$ と ${\frac {\;{d}y\;}{dt}}$ から ${\frac {\;{d}y\;}{dx}}=\psi (t)$ に至る計算の詳細を下に示しておく。

（▲）式を $t$ で微分した式は前掲の、 $\displaystyle {\frac {\;{d}x\;}{dt}}={\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}$ である。一方、（▲▲）式を $t$ で微分すると、

$\displaystyle {\frac {\;{d}y\;}{dt}}={\frac {1}{\;m\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}-{\frac {\ell }{\;m\;}}\cdot \!{\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}$

となる。従って、

$\displaystyle {\frac {\displaystyle \;{\frac {dy}{\;dt\;}}\;}{\displaystyle {\frac {dx}{\;dt\;}}}}=\displaystyle {\frac {\;\displaystyle {\frac {1}{\;m\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}-{\frac {\ell }{\;m\;}}\cdot \!{\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}\;}{\displaystyle {\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}}}$

である。更に計算してゆくと、

$\displaystyle {\frac {\displaystyle \;{\frac {dy}{\;dt\;}}\;}{\displaystyle {\frac {dx}{\;dt\;}}}}=\displaystyle {\frac {\;\displaystyle {\frac {1}{\;m\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}\;}{\displaystyle {\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}}}-{\frac {\;\displaystyle {\frac {\ell }{\;m\;}}\cdot \!{\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}\;}{\displaystyle {\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;{d}\phi (t)\;}{dt}}}}$

となるから、これを計算して整理すると、

$\displaystyle {\frac {\displaystyle \;{\frac {dy}{\;dt\;}}\;}{\displaystyle {\frac {dx}{\;dt\;}}}}=\displaystyle {\frac {\;\displaystyle {\frac {1}{\;m\;}}\;}{\displaystyle {\frac {1}{\;\ell +m\psi (t)\;}}}}-{\frac {\;\displaystyle {\frac {\ell }{\;m\;}}\;}{1}}=\displaystyle {\frac {\;\ell +m\psi (t)\;}{\;m\;}}-\displaystyle {\frac {\ell }{\;m\;}}=\displaystyle {\frac {\ell }{\;m\;}}+\psi (t)-\displaystyle {\frac {\ell }{\;m\;}}=\psi (t),$

である。すなわち，

$\displaystyle {\frac {dy}{\;dx\;}}={\frac {\displaystyle {\frac {dy}{\;dt\;}}}{\displaystyle {\frac {dx}{\;dt\;}}}}$

であるから、

$\displaystyle {\frac {dy}{\;dx\;}}=\psi (t),$

が得られたことになる。

例題1

陰関数型の関数 $F$ が、 $F{\bigl (}\phi (t),\;\psi (t){\bigr )}=\phi (t)-\psi (t)=0$ のとき、

$x=\!\int \!\left\{{\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;d\phi (t)\;}{dt}}\,\right\}dt+C,$

$k+\ell x+my=\phi (t)$

から、一般解を求めよ。 解きかたは、$\psi (t)=\phi (t)$ の関係を $x=\!\int \!\left\{{\frac {1}{\;\ell +m\psi (t)\;}}\cdot {\frac {\;d\phi (t)\;}{dt}}\;\right\}dt+C$ に適用すればよい。 すなわち、 $x=\!\int \!\left\{{\frac {1}{\;\ell +m\phi (t)\;}}\cdot {\frac {\;d\phi (t)\;}{dt}}\;\right\}dt+C$ として、積分すると、

$x={\frac {1}{\;m\;}}\cdot \ln {\bigl (}\,\ell +m\phi (t)\,{\bigr )}+C,$・・・・・（•）

を得る。 $x$ が $\phi (t)$ の関数として表示されたので、$y$ は、 $k+\ell x+my=\phi (t)$ を用いて $\phi (t)$ の関数として求めることができる。

計算の途中は省略して、求めた $y$ の数式のみを記述しておく。

$\displaystyle y={\frac {1}{\;m\;}}\left[\phi (t)-\left\{k+\ell \left({\frac {1}{\;m\;}}\cdot \ln {\bigl (}\,\ell +m\phi (t)\,{\bigr )}+C\right)\right\}\right].$・・・・・（••）

一般解 $x$ と $y$ が $\phi (t)$ の関数として表示されたので、次に、

$\displaystyle {\frac {dy}{\;dx\;}}={\frac {\displaystyle {\frac {dy}{\;dt\;}}}{\displaystyle {\frac {dx}{\;dt\;}}}}$

を確かめておく必要がある。（•）式と（••）式の両辺を $t$ で微分すると、

${\frac {dx}{\;dt\;}}={\frac {1}{\;\ell +m\phi (t)\;}}\cdot {\frac {\;d\phi (t)\;}{dt}},$

${\frac {dy}{\;dt\;}}={\frac {\phi (t)}{\;\ell +m\phi (t)\;}}\cdot {\frac {\;d\phi (t)\;}{dt}}$

が得られる。この、上記の 2式により、$\displaystyle {\frac {dy}{\;dx\;}}$ は、

$\displaystyle {\frac {dy}{\;dx\;}}={\frac {\displaystyle \;{\frac {dy}{\;dt\;}}}{\displaystyle {\frac {dx}{\;dt\;}}}}={\frac {\displaystyle {\frac {\phi (t)}{\;\ell +m\phi (t)\;}}\cdot {\frac {\;d\phi (t)\;}{dt}}}{\displaystyle \;{\frac {1}{\;\ell +m\phi (t)\;}}\cdot {\frac {\;d\phi (t)\;}{dt}}\;}}=\phi (t)$

となる。すなわち，

$\displaystyle {\frac {dy}{\;dx\;}}=\phi (t)$

が得られた。 例題1 は、$\psi (t)=\phi (t)$ の場合であるから、 上式 $\displaystyle {\frac {dy}{\;dx\;}}=\phi (t)$ は、常微分方程式 $\displaystyle F{\Bigl (}k+\ell x+my,\;\;{\frac {dy}{dx}}\;{\Bigr )}=0$ を解くとき、 始めに仮定した式 $\displaystyle {\frac {dy}{\;dx\;}}=\psi (t)$ と同一である。 次に、（•）式と（••）式から、 $\phi (t)$ を消去した式を示しておく。

$\displaystyle y={\frac {1}{\;m\;}}\left[{\frac {1}{\;m\;}}\{\exp(mx-mC\;)-\ell \;\}-k-\ell x\;\right].$ ・・・・・（•••）

上記の（•••）式は、$k+\ell x+my={\frac {dy}{dx}}$ の一般解である。なお、$C$ は積分定数である。

例題2

陰関数型の関数 $F$ が、 $F{\bigl (}\phi (t),\;\psi (t){\bigr )}=\phi (t)-\psi (t)=0$ の場合は、$\phi (t)=\psi (t)$ であるから、 最初の仮定 $k+\ell x+my=\phi (t)$ と ${\frac {dy}{dx}}=\psi (t)$ により、

$k+\ell x+my={\frac {dy}{dx}},$

が成り立つ。この式は、1階線型常微分方程式であるため求積法で解ける。 その一般解と、例題1 の一般解とが一致することを確かめよ。 （注：積分定数が 例題1 の場合とは等しくなく、異なる形をしている。非常に興味深い。研究する価値がある。（••••）と（★★）。）

上記の式 $k+\ell x+my={\frac {dy}{dx}}$ の一般解は、公式を用いて解くと、

$\displaystyle y=C\exp(mx)-{\frac {1}{\;m\;}}\left\{k+\ell \left(x+{\frac {1}{\;m\;}}\right)\right\},$・・・・・（★）

となる。なお、$C$ は積分定数である。 以下で、例題1 の（•••）式と （★）式を比較する。 まず、例題1 の（•••）式を展開すると、

$\displaystyle y={\frac {\;\exp(-mC)\;}{m^{2}}}\cdot \exp(mx)-{\frac {\;\ell x\;}{m\;}}-{\frac {k}{\;m\;}}-{\frac {\ell }{\;m^{2}\;}},$・・・・・（••••）

を得る。 次に、例題2 の（★）式を展開すると、

$\displaystyle y=C\exp(mx)-{\frac {\;\ell x\;}{m}}-{\frac {k}{\;m\;}}-{\frac {\ell }{\;m^{2}\;}},$・・・・・（★★）

となる。 上記（••••）式、および（★★）式は、いずれも、 $k+\ell x+my={\frac {dy}{dx}}$ を満たす。