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

定数係数線形非同次微分方程式

似たようなタイトルが続くね。
作成:2013/7/9

非同次形で言えること

前々回,前回と,定数係数の同次形を考えてきた.今回は非同次形を考えてみよう.例えば 2 階の定数係数非同次形というのは次の形で書かれる. 数式 右辺が 0 ではない辺りが今までとは違っている.また,右辺まで定数というわけではなく,広い範囲で考えてxの関数だとしておこう.

右辺のq(x)が邪魔なせいで,今までやったような解を代入しても成り立たない.しかし,この式を成り立たせる解z(x)が何か一つ,どうにかしてうまい具合に見つかったとしよう.すると,q(x)=0とした場合に成り立つ解と,z(x)とを足してやったものは,ことごとくこの方程式を満たすことになる.(1) 式の一般解は次のような形になると言っているのである. 数式 z(x)には任意定数が付いていないし,付けると成り立たない.(1) 式の右辺にあるq(x)のことはz(x)だけに任せた形になるわけだ.

というわけで,何とかして一つでいいからz(x)を求めることが必要になる.


係数変化法

手順は簡単である.まずは (1) 式の右辺を 0 と置いてやって,前々回の方法で二つの基本解y_{1}(x)y_{2}(x)を求めてやる.次に 数式 という形の解を仮定してやる.この段階でなぜかと問うてはならない.この方法でうまく行くのを見てから,もしまだ納得が行かなければ考えて欲しい.(3) 式を 1 階微分してやろう. 数式 ここで,もし 1 行目のカッコの中身について, 数式 という条件が満たされていれば,(4) 式は 2 行目だけになって, 数式 となるから,これをさらに微分してやる. 数式 ここで,もし 1 行目のカッコの中身について, 数式 という条件が満たされていれば,(7) 式は 数式 となる.ここでz(x)が (1) 式の解であるというので (3) (6) (9) 式を (1) 式に代入してやると, 数式 となって,確かにz(x)は (1) 式を満たしていることが分かる.なぜこんなことが成り立ったかと言えば,途中で使った (5) 式と (8) 式の条件のお陰である.もう一度並べて書いてみよう. 数式 これを満たすようなA(x)B(x)を求めてやれば良いのである.未知関数はA'(x)B'(x)だけであり,これはただの連立方程式だ.

A'(x)B'(x)とが求まったら積分してやってA(x)B(x)が求まるから,それを (3) 式に代入してやればいい.A(x)B(x)は積分して求めるので任意定数が付いた形になっているが,それを残したまま (3) 式に代入してやれば勝手に (2) 式と同じ形式のものが出来上がるので心配は要らない.


階数を増やす拡張

今は 2 階の場合をやったわけだが,階数が増えても全く同じ手法が通用するだろう.

まずq(x)=0と置いてやって,前回の記事の通りにn個の同次形の基本解y_{n}(x)を求めてやり,次のようなn個の条件式を作る. 数式 階数が増えてもA_{n}(x)の 1 階微分しか出てこないのは本当に有難い.これは微分しては,条件式を作って 0 と置いておいてやって,また微分ということを繰り返すからそうなるのである.

このような式の列挙ではごちゃごちゃして分かりにくいかも知れないから行列を使って表現してみよう. 数式 この連立方程式を解いてA'_{1}(x),…,A'_{n}(x)を求めてやり,それを積分して求めたA_{1}(x),…,A_{n}(x)を任意定数の付いたままで, 数式 に代入してやれば,それが一般解である.



趣味の物理学書店

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