微分方程式の解に含まれるパラメータの数
上で述べたように、微分方程式には、微分方程式だけでは決まらないパラメータが必ず含まれる。それは微分方程式が局所的情報を表す式であることから必然的なのである。微分方程式を解く時にもこの点は大事なので、そのパラメータの数について考察しておこう。
dydx=yの解は y=Aexを例として考えよう。このときのパラメータAは
dydx=y(左辺にyを、右辺にxを集める)dyy=dx(両辺を積分)logy=x+Cの計算過程においてdyy=dxを積分するときの積分定数から現れたここで「両辺をyで割る」という計算をやっているがy=0の場合、これは許されない。ここでは暗黙のうちにy≠0を仮定している。と考えてよい。さらに、logy=x+Cは
y=ex+C=eCexとなるから、A=eCとすれば上の式と同じになる「この形だとAは負になれないのでは?」と心配する人もいるかもしれないが、Cがiπという虚部を持っていれば、Aは負にもなるので気にしなくてよい。。注釈に書いたように、ここまでの計算ではy≠0が仮定されていたが、幸いなことにy=0はこの一般解y=Aexの、A=0の場合に含まれているので、y≠0の条件は外してよいことになる。
左辺に積分定数をつけても、結果は同じなのである。もし左辺に積分定数をつけたとすれば、左辺の積分定数と右辺の積分定数は別の定数なのでそれをそれぞれA,Bとして時々、積分定数をどっちもCにしてlogy+C=x+Cとして両辺で打ち消してしまう、という計算をする人がいる(←ここは驚くか笑うかするところ)が、積分定数は左辺と右辺それぞれにおいて「任意の数」だから、両辺で一致する理由はない。
logy+A=x+Blogy=x+B−Aとなるが、AもBもまだ決まっていない数であり、しかも結果にはB−Aという組み合わせでしか出てこない。つまり、答えを出すためにはA,Bそれぞれを求める必要はなく、B−Aだけを求めればよいから、C=B−Aとおいて1つの積分定数と思えばよい。
非常に簡単な二階微分方程式(ddx)2f(x)=ddxf(x)を同様に解いてみよう。
(ddx)2f(x)=ddxf(x)(両辺を不定積分)ddxf(x)=f(x)+C(f(x)+C=yと置く)ddx(y−C)=y(ddx(−C)=0を使って、さらに変数分離)dyy=dx(もう一度積分)logy=x+D(eの肩に乗せて)y=ex+D(f(x)に戻して)f(x)=−C+ex+D積分を二度やったため、積分定数C,Dの二つが結果に現れる(パラメータの数は2)。
微分方程式を解くとは積分すること、と考えると「n階微分方程式なら不定積分をn回繰り返せば解ける」と言えて、結果はn個の積分変数を含む。上の具体例を見ると、確かに一階微分方程式の解は1個の、二階微分方程式の解は2個の積分定数を含んでいる。
結論として、n階微分方程式の解は常にn個の「微分方程式だけでは決まらないパラメータ」を含んでいることになるただし、1xの積分のところで示したように、途中で関数が定義できない点(この例の場合x=0)があると積分一つに対して二個の積分定数が出て来ることもあるので、そのような場合には注意が必要である。。
こういう考え方もできる。解析的な関数(つまり、テイラー展開できる関数)に限って考えると、「関数を決める」というのは、
f(x)=∞∑n=01n!f(n)(x0)(x−x0)nの係数f(n)(x0)を全て決めることである(少なくともテイラー展開の収束半径の内側ではこれで十分である)。
m階微分方程式は(適切な変更を行った後)
f(m)(x)=(f(m−1)(x),f(m−2)(x),⋯,f(1)(x),f(0)(x))の式のように書くことができる。さらにこれをどんどん微分することで、
f(m+1)(x)=(f(m)(x)⏟微分方程式を代入,f(m−1)(x),f(m−2)(x),⋯,f(1)(x),f(0)(x))の式=(f(m−1)(x),f(m−2)(x),⋯,f(1)(x),f(0)(x))の式のようにm階より高い階数の微係数も求めることができる。これらを使ってf(m)(x0)をそれより微分階数の低い係数と使って書き直すことができるから、f(x)の表現には、mより低い階数の微係数(f(m−1)(x0),f(m−2)(x0),⋯,f(1)(x0),f(0)(x0))だけが「決まらずに」残る。
たとえば一階微分方程式を満足する関数であればf(x0)のみが、二階微分方程式を満足する関数であればf(x0)とf′(x0)の二つだけが微分方程式だけでは決まらないパラメータとなる。
一階微分方程式で正規形の場合で、「決まらないパラメータ」の意味を考えておこう。
dydx=f(x,y)という式が与えられていると、x-y平面上である点(x,y)を指定すると、その点における関数のグラフの傾きdydxがわかる、という式になっている。つまり、dydx=f(x,y)という式は、各点各点において「グラフの線はどちらに伸ばすべきか」を与えている。
たとえば最初に考えた微分方程式はdydx=yと書けるが、その解のグラフは、各場所においてy座標と同じ傾きを持っている曲線である。
そういう曲線を次々と描いていくと、上のグラフにあるように全平面を埋め尽くしていく。
たとえばx=0の時y=1というふうに「出発点」を決めると、この場合はy=exという線(グラフでは1本だけ太い線で表現した)の上を進んでいくことになる。
一階微分方程式が指定するのは傾きのみであるから、出発点(上の例ではx=0から始めたが、実はどの場所でもよい)を指定すれば曲線は一つ決まる。別の点を出発点にすれば(たまたま同じ線上の2点を選ばない限り)また別の線が引ける。こうして、微分方程式だけからは決まらないパラメータが解には入っていることになる(後で、それを「初期条件」などで決めていく方法について述べる)。
二階微分方程式では、傾きではなく「曲がり具合」が微分方程式によって指定され、「場所」と「傾き」が微分方程式では決まらない量になる。
変数分離できる一階微分方程式
まず微分方程式がどういうものかに慣れることが必要だと思うので、以下では、微分方程式の中でも比較的簡単(でも応用範囲は広い)な「変数分離できる一階微分方程式」の具体例を考えていこう。
変数分離はいつでもできるとは限らないことに注意しよう。たとえばdydx(x)=x+yという簡単な場合でも左辺にyだけを集めることはできない(この微分方程式は解ける。つまり変数分離できなくても解ける時は解ける)。
以下この節では「変数分離できる場合」に限って話をする。
例として上げたdydx=yは変数分離できる微分方程式の例でもある。上でやったように、「変数分離した後で積分」という方法で解くことができた。
もう一つ、簡単な例を示そう。
dydx=−xyという式(前に図で考えた微分方程式で、答は円であった)は書き直すとydy=−xdxと変数分離できる。これを積分すると、
∫ydy=−∫xdxy22=−x22+Cとなる。Cは積分定数である。結果を整理すると、
x2+y2=2Cという円の式が出てくる(半径は√2Cである)。
次にdydx=yxを同様に解いてみよう(図で考えた時、この解は「原点を通る直線」であった)。まず変数分離して、
dyy=dxx(積分)logy=logx+C(両辺をeの肩へ)y=xeCとなり、確かに(傾きeCの)直線が解である(図を描いて考える方がすっとわかる)。