物理を解説 ♪
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン
ツイッター用のシェアボタン フェイスブック用のシェアボタン はてなブックマーク用のシェアボタン ライン用のシェアボタン

線形微分方程式の級数解法

教科書の後の方で出てくるからと言って、難しいわけじゃないんだよ。
作成:2013/6/6

原理は簡単

非線形の微分方程式は解けないことの方が多いが,線形ならば何とかして必ず解ける・・・などと匂わせる説明をしたのだった.前回は 1 階の場合の解き方を説明しただけだから,2 階以上でもちゃんと解けることを示す必要がある.

しかしそう単純ではなく,いつでも頼れる唯一の方法があるわけではない.方程式の形が単純な場合には,それなりに楽な解法が幾つも存在しているのである.そういうものを一つ一つゆっくり紹介していくといつまで経ってもスッキリしなくて疲れてしまうので,必ず解ける最終手段をいきなり説明してしまおう.

線形微分方程式は次のような形をしているのだった. 数式 ここでy^{(n)}というのは未知関数y(x)n階微分であり,y^{(0)}は微分していないままのy(x)のことである.またp_{1}(x),…,p_{n}(x)q(x)xのみの既知の関数で,これらの関数の形によって問題が決まるのである.

予備説明はこれで終わり.ではそろそろ始めよう.テイラー展開についてはすでに学んでいるだろう.まず,未知関数y(x)を次のように表せると仮定する. 数式 すると,y'y''も次のように計算すれば良いことになる. 数式 これ以上の微分を繰り返しても同じような形に書けるはずだ.既知の関数p_{1}(x),…,p_{n}(x)q(x)についてもテイラー展開できる. 数式 これらの全てを (1) 式へと代入すれば・・・.うーん,その結果を正確には書きたくないな.あまりに複雑で,本当に言いたいことが分かりにくくなってしまうからだ.とにかく全て展開して,x^kの次数が同じ項を集めて整理すれば次のような形に書けるはずである. 数式 この式がxの値に関係なく成り立つようにするためには,それぞれの(…)の中身がどれも 0 でなくてはならない.その中身が具体的にどんな感じの式になるかを思い描けることが大事なのだが,いきなりn階の場合について考えるのは酷なので,後でとりあえず 2 階の場合でやってみよう.

それらは連立方程式になるわけだが,ある程度のパターンを持ったものになっている.幾つかの任意定数を仮定することで,順々に全てのa_kの値を定めて行くことができる.要するに漸化式の形になっているのでそれを解くことになるわけだ.

そのような手順で求めたa_kを (2) 式に当てはめることで,y(x)が級数の形で導けることになる.

ただし,この方法はp_{1}(x),…,p_{n}(x)q(x)がテイラー展開できないような関数である場合には使えない.線形なら必ず解けるようなことを言ったのはそのような特殊な場合については目をつぶってのことである.

また,今の説明ではx=0の周りでのテイラー展開を行っているが,問題によってはp_k(x)x=0の周りではテイラー展開できないということもあるだろう.その場合はx=dの周りで展開して(x-d)^kの冪で表せばいい.今回の話でx^kとなっている部分が(x-d)^kになるだけで,全く同じ手続きである.


確認作業

2 階から始めようと思ったが,思ったよりも面倒だったので 1 階から試してみることにする.1 階の線形微分方程式は前回解いたので,それが今回の方法だとどうなるのかを比較するのも悪くない.

1 階の線形微分方程式は次の形をしている. 数式 これに (2) (3) (5) (6) 式を代入する. 数式 2 行目が無限級数どうしの積を行うことになるので一番大変だと思う.しかしそれぞれの係数の添字の合計がx^kkと等しくなることに気付けば,列挙するだけだからそう難しくもない.とにかく,この結果をx^kkが同じものどうしで集めて整理すると,次のような多数の等式が出来上がる. 数式 まずk=0の場合を見ると,b_{0}c_{0}は既知なので,結局は未知の変数であるa_{0}a_{1}の関係を表した式だということになる.これは一方の値,例えばa_{0}の方を任意定数として定めてやってa_{0}=Cとでも置いてやれば,そこからa_{1}の値が決められるようになる.

次のk=2の式を見ると,新たに登場する未知変数はa_{2}であり,他の値はすでに定まっているからそれを元にこの値も決められる.k=3の場合はa_{3}の値が決まるし,以下も同様である.


推して知るべし

2 階の微分方程式の場合にはどうなるか,やらなくても想像が付くだろう.k=0の式でいきなりa_{0}a_{1}a_{2}の関係が出てくるだろうから,この中の二つを任意定数として定めてやれば,残る一つの値も定まる.それ以降の式ではそれを元にして一つ一つ値を決めていけばいいのである.

先ほど,漸化式を解いてやれば良いなどと書いたが,先へ進むほど項の数が爆発的に増えるので,とても考える気になれない.しかしp(x)q(x)が無限に続く級数ではなくてかなり単純な形をしていれば・・・,例えば (5) 式や (6) 式のようではなく, 数式 のような形をしていれば,既知の係数はb_{1}=2c_{0}=3以外は全て 0 であり,項の数はいつも一定になるから,何か単純な漸化式のようなものが作れることだろう.しかし今回は「原理的に必ず解ける」ということだけ言いたかったのでこれくらいにしておこう.今回の手法を具体的に利用する話は,また後の方で「特殊関数」の話題と絡めてするつもりである.

次回からは,方程式が単純な形をしている時にだけ使える,実用的な解法を色々と紹介していこう.



趣味の物理学書店

趣味で量子力学2の広告バナー