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

ベッセルの微分方程式

ベッセル関数、それは難解な特殊関数の代名詞。
作成:2015/4/10

概要説明

次の形の微分方程式を「(l位の)ベッセルの微分方程式」と呼ぶ. 数式 lは実数である.lが正の場合と負の場合は同じ結果になるわけだが,独立な解が二つ得られるので,一方の解をlが正の場合に割り当て,他方をlが負の場合に割り当てるような定義を採用するのである.こうして決めた解のことを「第一種のベッセル関数」と呼ぶ.

ではlが 0 の場合にはどうなるかというのが気になるだろう.その場合にもちゃんと独立な解は二つ得られるのである.そのうちの一方は「第一種のベッセル関数」で表せるが,それとは独立な解を「第二種のベッセル関数」として別に定義する.

実はlが整数の場合はlが正の場合と負の場合の「第一種のベッセル関数」はそれぞれ独立な解にはなっていない.それで,lが 0 の場合と同じように,「第一種のベッセル関数」と「第二種のベッセル関数」の組み合わせが解になるのである.

「第一種のベッセル関数」をJ_l(x)という記号で表し,「第二種のベッセル関数」をY_l(x)という記号で表す.「第二種のベッセル関数」は「ノイマン関数」と呼ばれることもある.その場合にはN_l(x)という記号で表されることもあるが,この記事では混乱を避けるために使わない.

第二種のベッセル関数のことをノイマン関数と呼ぶ流儀の教科書では,第一種のベッセル関数のことを単にベッセル関数と呼ぶことが多い.

ベッセルの微分方程式は,物理ではラプラス方程式やヘルムホルツ方程式を円筒座標で解く時などに現れる.断面が円形の導波管内の電磁波や,円形膜の振動などである.


解法

(1) 式を級数解法で解くために両辺をx^2で割ってみると, 数式 という形になるのだが,級数解法が使えるための条件は各項のxの関数がテイラー展開できることだった.ところがこれらはx=0で特異点になっており,x=0を中心としたテイラー展開は無理である.中心をずらせば何とかならなくもないが,複雑になりそうである.

ところがだ,特別な場合に限ってはうまく行く計算法が知られているのである.それは「確定特異点での級数解」と呼ばれたりするものだが,その一般論を話していると長くなるので,ここではその知識をさりげなく使いながら解いてしまおう.

まず,この方程式の解y(x)が次のように展開できると仮定してみる. 数式 いつもの級数解法で使う仮定とは違って全体にx^rが掛けられており,まるで根拠がないかのように見えるが,こうすると結果的にうまく行くのである.rは整数だとは限らない.実数だと考えておこう.また,a_{0}≠0だとしておく.もしa_{0}=0ならrの値が1だけ減って次の項のa_{1}が代わりの役を演じることになるだけなので考えるだけ無駄だからである.これを微分したものは, 数式 となり,もう一度微分して, 数式 となる.これらを (2) 式に代入してやれば, 数式 となるが,括弧の外にある式を掛けてやれば,各行の先頭のxの次数が合うのではなかろうか. 数式 これならxの次数ごとにまとめてやることができて,右辺が 0 なので,その各項の係数が 0 に等しいということになる.例えばx^{r-2}の項ばかりを集めると, 数式 であり,a_{0}≠0であるからr=±lというのがその解だ.

このままでは分かりにくいので,とりあえずr=lの方だけを使って (3) 式の仮定からやり直そう. 数式 だと仮定する.これを微分することで, 数式 を得る.これらを (2) 式に代入すると, 数式 であり,もう少しまとめると,次のようだ. 数式 もう少しまとまりそうだ. 数式 まだ行けるな. 数式 ここで二つの項のxの次数を合わせておきたいので,少し小細工をする.第 2 項を次のように変えても意味は変わらないだろう. 数式 さらに第 1 項はn=0の時には 0 になるのでn=1から始めてもいい.第 2 項もそれに合わせてn=1から始めようとすると,もともと存在していなかったa_{-1}というものについて考えることになる.しかしa_{-1}=0だと定義しておけば問題ないのではなかろうか. 数式 それで,次のような式にまとめられる. 数式 要するに,次のような関係が成り立っていればいい. 数式 これを使って順次求めていけばいい.a_{1}を求めたければa_{-1}を使って導くことになるが,a_{-1}=0と定めたのでa_{1}も 0 になる.同様にa_{1}から求めるa_{3}も 0 だ.以後同様に,(4) 式の奇数番目の項は全て 0 になる.

では偶数番目の項はどうなるかというと,a_{0}から順に作っていけばいい. 数式 結局これは,次のようにまとめたらきれいなのではないだろうか. 数式 しかしまだ美しさに欠ける.分母には(l+k)から1つずつ値を減らしたものを掛けてゆく部分があって,それが(l+1)までで終わっている.もし,この続きとしてlの階乗があれば,この部分の全体を(l+k)!という形にすっきりまとめて書くことが出来るのではないだろうか.a_{0}の値を決める制限は出てきていないので,これを使ってそれを実現してやったらどうだろう.a_{0}=1/(l!)だと定めてやるのだ.それを (4) 式に当てはめると次のようになる. 数式 しかしもっとすっきりさせることができるだろう.ここではxの次数がx^{l+2k}と表されている.もし分母に2^lを掛けてやったら,2^{2k}と一緒にして2^{l+2k}と書ける.そうすればx^{l+2k}と次数が同じなのでまとめて表せるようになるのではないだろうか.そこで,先ほどの考えを変更してa_{0}=1/(2^l l!)だと定めてやることにする.そうすれば次のように書ける. 数式 ところがこれでは少し問題がある.階乗というのは整数にしか使えないのだったから,これではlが整数の場合にしか対応できていない.しかし世の中には階乗の性質をそのままに,整数以外にも使えるように拡張して作った「ガンマ関数」というものがあり,nが整数の場合にはガンマ関数との間に 数式 という関係が成り立っている.そこで,さらに考えを変更して,a_{0}=1/(2^l Γ(l+1))だと定めてやることにする.そうすれば, 数式 となり,lが実数の範囲であっても対応することができるのである.

さて,ここまではr=lの場合の話であったが,もう一つ,r=-lの場合もあったのだった.しかしそれはここまでの話と全く同じであり,l-lに置き換えるだけで実現できてしまう.結果の式も全く同じ形だ.

(5) 式を第 1 種の(l位の)ベッセル関数と呼んで,J_l(x)と表す.もう一度書いておこう. 数式



ベッセル関数は独立か

ベッセル関数のlが整数でない場合には,J_l(x)J_{-l}(x)は独立である.それはxの次数を見れば分かるが,x^{(l+2k)}x^{(-l+2k)}とは全く共通な点がないからだ.例えばl=1.3だとしたら,x^{1.3},x^{3.3},x^{5.3},x^{7.3},…となり,l=-1.3だとしたら,x^{-1.3},x^{0.7},x^{2.7},x^{4.7},…となり,共通項がない.

それで,lが整数以外の場合には任意定数をC_{1},C_{2}として, 数式 というのがベッセルの微分方程式の一般解となる.

ところがlが整数の場合にはJ_l(x)J_{-l}(x)は独立ではない.次のような関係になっていることが言えてしまうからだ. 数式 それを示すのは簡単である.ガンマ関数Γ(n)というのはnが 0 または負の整数になるところで値が+∞-∞のいずれかになってしまうという性質があり,(6) 式のlに負の整数値を入れると,最初の幾つかの項は 0 になってしまう.それで次のように変形できる. 数式 とにかくこういうわけだから,lが整数の場合には (7) 式はベッセルの微分方程式の一般解ではない.J_l(x)とは独立な基本解をもう一つ探さなくてはならない.


別の解を探す

lが整数の場合のJ_l(x)とは独立のもう一つの基本解を探すには,定数変化法を使う.(1) 式の解を 数式 だと仮定してやって (1) 式に代入してやってW(x)を求めるという方針で行く. 数式 だから,これを (1) 式に代入する.いや (2) 式を使おう. 数式 J_l(x)は (2) 式の解だから上の式の幾らかは消せる.斜線を引いた部分の合計は 0 であるというのが (2) 式から分かるからである.残った項を集めたのが次の式だ. 数式 いやはや,かなり単純な条件になったのではないかもう少しまとめてみよう. 数式 この両辺を積分すれば, 数式 であり,任意定数の部分は省いて考えてもいいので, 数式 となる.ここに (6) 式を代入して積分すれば一気に答えを得られそうだが,複雑になりすぎて扱いに困る.分母に無限級数が来るなんて,しかもその 2 乗を計算するなんて,さらにそれを積分するなんて,どう進めていいか分からない.

そこで,かなり遠回りになってしまうが別の道を探ることにしよう.まずは一度大雑把に式の形を見積もってみる.(6) 式を見返してもらうと分かるが,J_l(x) 数式 という形式になっている.それで, 数式 という形になるはずだ.この右辺の後ろの分数の方には特異点はないので普通にマクローリン展開ができるだろう.それで係数がどうなるかは不明だが,次の形に書ける. 数式 これを積分すれば, 数式 のようになる.Cは積分定数である.ややこしく見えるかも知れないが,ここで使っているb_{0},b_{1},b_{2},…は仮のものであり,全て不明なので,別の記号で書き換えればもっとすっきりする. 数式 右辺の最初にあるb_{2l}は 1 だと考えても構わない.どうせ後からこの全体に任意定数を掛けたものが解になるのである.全体をb_{2l}で割ってやれば跡形もなく消せるだろう.いちいち記号を変えるのがもったいないのでそのまま使うが,次の形になっているということだ. 数式 さて,後から説明するとわけが分からなくなるので今のうちに言っておくが,実はl=0の場合にはc_{0}は任意定数である.なぜならl=0の場合,一般解は二つの任意定数を使って, 数式 と表せると仮定しているからだ.この中にJ_{0}(x)に比例する項が 2 つもあり,どちらにも任意定数が付いている.c_{0}は幾つであっても最初の項に吸収されてしまうのだ.だから簡単にするためには,c_{0}=0と考えておいてやれば構わない.ただしこれはl=0の場合に限った話だ.

さて話を戻そう.(8) 式の仮定によれば,この (10) 式にJ_l(x)を掛けたものが解の候補なのだった.J_l(x)は (9) 式のように書けるのだから,また別の係数を使って表現すれば,解は次のような形をしていると言える. 数式 この第 2 項をu(x)と表現することにすると,要するに, 数式 であり,今度はこのu(x)を求めたいわけだ.この式を微分すると, 数式 であり,これを (2) 式に代入してやる. 数式 J_l(x)はもともと (2) 式の解だから消えてもらった.意外に単純な式が残る. 数式 今はlが整数だという仮定なのでJ_l(x)の定義で使っているガンマ関数の部分も階乗で表せて,2{J_l}'(x) /xの部分は 数式 となる.先ほど導いたことからすると,u(x) 数式 という形をしているのであり, 数式 である.これらをすべて (11) 式に代入すれば, 数式 となるが,左辺は割りと綺麗にまとめることが出来そうだ. 数式 この左辺は,a_{-l-1}=0,a_{-l-2}=0とすれば, 数式 とまとめることができるから, 数式 という関係が成り立っていることが分かる.これを解くのは大変そうだが,l=0の場合ならずっと簡単になるだろう.l=0の場合とそうでない場合とで問題の質も変わってくる.なぜならl=0なら負冪の項を考えなくて済むからだ.そこから試してみよう.問題は次のように比較的単純なものに変わる. 数式 ただしa_{-1}=a_{-2}=0である.右辺にはxの奇数乗の項はないから,左辺のnが奇数の時には 数式 が成り立っており,a_{-1}=0ゆえにa_{1}=0であり,続いてa_{3}=0が導かれ,この調子でずっと 0 である.一方nが偶数の時には (13) 式の左辺のnと右辺のkの間にn=2k-2の関係がある時にxの次数が一致するので, 数式 という関係が成り立っている.もう少し整理すると, 数式 であり,l=0の時にはa_{0}=0だと考えて良いという話を思い出して, 数式 となり,だんだんとパターンが見えてくる.これは結局は次のようにまとめることができる. 数式 以上より,探していたもう一つの基本解はl=0の場合には次のように表されることになる. 数式 もうかなり面倒なので検証を省くが,lが 0 でない正の整数の場合には次のようになる.(私自身も嫌になって確認ができていない) 数式 これらY_l(x)が「第二種ベッセル関数」である.


はるかにシンプルな表現

せっかくここまで苦労してlが整数の場合の「第二種ベッセル関数」の級数表現を求めたのだったが,実はこれをかなりシンプルに表すことができる式があるのである. 数式 驚くほどシンプルで,本当に同じものなのかという疑いたくなるほどだ.実は全く同じというわけではない.上で求めた級数解とは定数倍だけ違っていたり,J_l(x)に比例する項が余計に付いていたりというわずかな違いはあるのだが,それは許される範囲内の違いであり,大きな問題でもない.それより大きな問題がある.この式をよく見ると,lが整数の場合には分子も分母も 0 になってしまって値が不定なのではなかろうか.実は,lが整数の時には,lの値が整数に近づくときの極限での値を定義として採用するのである. 数式 このような手続きはかなり面倒だと思うかもしれないが,この場合には「ロピタルの定理」を当てはめることが出来る.それは次のようなものだった. 数式 つまり,次のような計算をして求めることができるのである. 数式 この結果が級数解と本当に同じになるかどうかの確認は省略しよう.(私もここまで書くのに疲れてしまって面倒になったのでやっていない.)先ほども書いたように,全く同じにはならない.

実はこの形式の解はlが整数でない場合にも使える.その場合には (14) 式をそのまま使えばいいのであって,ややこしい極限を考える必要もない.式の中にはJ_l(x)J_{-l}(x)とが含まれているから,これはJ_l(x)とは独立な解である.だからlが整数でない場合にも,J_l(x)Y_l(x)との組み合わせを基本解として採用することも可能なのである.

結局,lが整数であろうとなかろうと, 数式 というのをベッセルの微分方程式の一般解として採用することもできるというわけだ.もちろん,lが整数でなければ (7) 式を使った方がよっぽどシンプルだけれども.



趣味の物理学書店

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