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

積分因子

こんなやり方を紹介されても困っちゃうよね。
作成:2013/5/30

この方法は凡人には無理かも

微分方程式が完全形になっていなかったとしても,諦めるのは早いかも知れない.全体に何らかの関数M(x,y)を掛けることで完全形の微分方程式として解ける可能性が残されている. 数式 これが完全形であるための条件を調べてみよう. 数式 どうにも分かりにくい.こんなややこしい条件を満たす関数M(x,y)が見つけられるとは思えない.しかしたまたま運良く見つけた人,あるいは直感力に優れた人がこの手法を使って解いている場面に出くわすかも知れない.その時のために,こういう手法もあるのだと知っておくといいだろう.

この (1) 式の条件を満たすM(x,y)のことを「積分因子」と呼ぶ.


凡人にもチャンスあり

では運にも勘にも見放され,粘る気力も続かないという凡人はこの手法とは無縁なのだろうか.いや,二度だけ「スピードくじ」を引くことが許されている.それは外れることもあるのだが,その場合は潔く諦めよう.

Mxのみの関数であると仮定してみたらどうだろう.(1) 式の左辺第 1 項はなくなるので,もっと簡単になるかも知れない. 数式 Mxのみの関数だと仮定しているのだから,右辺にyが出てくるのはおかしい.右辺を計算してみて,たまたまxのみの関数になっていたら,それは当たりくじだ.それをr(x)とでも表現しよう. 数式 これは関数M(x)についての微分方程式になっている.もう少し変形してみよう. 数式 これは「1 階線形微分方程式」と呼ばれるものであり,次回説明するつもりでいたものが一足早く出てきてしまったしかしこれはその中でも特別簡単な形になっており,実は変数分離形でもあるから,すでに学んだやり方で解けてしまう.これを解くのは今回の目的ではないのでこれ以上はやめておこう.とにかくこのようにしてM(x)を導く方法があるということだ.

さて先ほど,くじは二度引けると言った.もう一つのチャンスはMyのみの関数だと仮定してみることである.やることの内容は先ほどと変わらないのだが,一方がダメでもどちらかがうまくいくことがあるし,どちらもうまくいくこともあるし,どちらもダメなこともある.(1) 式の右辺第 2 項が消えるので,次のように変形して行ける. 数式 もしこの右辺がyのみの関数になっていれば当たりくじである.右辺をr(y)とすると,次のような微分方程式を解けば良いだろう. 数式 これも「1 階線形微分方程式」であるが,変数分離形だと考えれば解けてしまう.変数がxではなくyになっているが,ここでのyxと同じただの変数と思えば良い.


具体例

次の微分方程式を解いてみる. 数式 これは変数分離形でも同次形でもない.完全形かも知れないと思って「第1項をy,第2項のカッコの中をxでそれぞれ偏微分する」ことを試してみても, 数式 となって条件を満たしていない.実は,両辺にxyを掛けてやれば完全形になるのである.そんなアイデアをどうして思い付くことが出来るのかと聞かれても分からない.(1) 式に合うものをたまたま思い付いたのである. 数式 試してみると, 数式 となっており,確かに完全形である.f(x,y)を探るために (3) 式の第 1 項をxで積分すると, 数式 と推定できる.また,(3) 式の第 2 項のカッコの中をyで積分することで, 数式 だと推定できる.見ての通り,互いに補うようなものはないからこのままでいい.任意の定数Cを使って,解は次の形に表せる. 数式 まとめた方が分かりやすいだろう. 数式 この例題の積分因子がxyであることをたまたま思いついたと書いたが,実はそんなものは全くの嘘であり,そうなるように (1) 式を分析して必死に例題を考えたのである.今回の具体例は解き方を教えることよりも,「積分因子を使ったときに何が起こるか」を試しに見てもらうことが主目的である.誰でも都合よく積分因子を思い付いて解けるものではないから,「諦め」を味わって欲しい.そのために,今回紹介した「凡人にも解ける方法」もわざわざ封じておいた.



趣味の物理学書店の案内バナー

EMAN物理note出張所の案内バナー