ベッセル関数
ベッセル関数（ベッセルかんすう、英: Bessel function）とは、最初にスイスの数学者ダニエル・ベルヌーイによって定義され、フリードリヒ・ヴィルヘルム・ベッセルにちなんで名づけられた関数。円筒関数と呼ばれることもある。以下に示す、ベッセルの微分方程式における[math]y(x)[/math]の特殊解の1つである。
- [math]x^2 \frac{d^2y}{dx^2}+x \frac{dy}{dx}+(x^2-\alpha^2)y=0[/math]
上の式において、[math]\alpha[/math]は、任意の実数である（次数と呼ばれる）。[math]\alpha[/math]が整数[math]n[/math]に等しい場合がとくに重要である。
[math]\alpha[/math]及び[math]-\alpha[/math]はともに同一の微分方程式を与えるが、慣例としてこれら2つの異なる次数に対して異なるベッセル関数が定義される(例えば、[math]\alpha[/math]の関数としてなるべく滑らかになるようにベッセル関数を定義する、など）。
そもそもベッセル関数は、惑星軌道の時間変化に関するケプラー方程式を、ベッセルが解析的に解いた際に導入された。
Contents
応用
ベッセル解は、ラプラス方程式または、ヘルムホルツ方程式の、円柱座標系または極座標系における分離解として見出される。したがって、ベッセル関数は、電波伝播や静電位差などの解を求めるときに重要なものとなっている。（円柱座標系においては、整数次数[math]\alpha=n[/math]のベッセル関数が解として得られる。極座標系においては、半整数次数[math]\alpha=n+1/2[/math]のベッセル関数が解として得られる。）
例えば、
ベッセル関数には、また信号処理のような他の問題のための有用な特性がある(例えば、FM合成、カイザー窓やベッセルフィルタなど)。
定義式
ベッセルの微分方程式は2階の線形微分方程式であるので、線形独立な2つの解が存在するはずである。しかしながら、解を議論する状況に応じて解の様々な表現が便利に使われている。代表的ないくつかの解の表現について以下で説明する。
第1種及び第2種ベッセル関数
これらの関数が、ベッセル関数群としては、最も一般的な形式である。
- 第1種ベッセル関数
- 第1種ベッセル関数は [math]\displaystyle J_\alpha(x)[/math] と表記される。[math]\displaystyle J_\alpha(x)[/math] はベッセルの微分方程式の解であり、[math]\displaystyle \alpha[/math] が整数もしくは非負であるとき、[math]\displaystyle x=0[/math] で有限の値をとる。[math]\displaystyle J_\alpha[/math] における、特定解の選択及び正規化は定義された後に、後述する。第1種ベッセル関数はまた、[math]\displaystyle x=0[/math] のまわりでのテイラー展開（非整数の[math]\displaystyle \alpha[/math] に対しては、より一般にべき級数展開）によって定義することもできる。
- [math]\displaystyle J_\alpha(x)=\sum_{m=0}^\infty \frac{(-1)^m}{m!\Gamma(m+\alpha+1)} \left ( \frac{x}{2} \right)^{2m+\alpha}[/math]
- 非整数の [math]\displaystyle \alpha[/math] に対しては、[math]\displaystyle J_\alpha(x)[/math] と[math]\displaystyle J_{-\alpha}(x)[/math] とが、ベッセルの微分方程式に対する線形独立な2つの解を与える。整数の[math]\displaystyle \alpha[/math] に対してはこれらは線形独立な解を与えない。なぜなら、整数[math]\displaystyle n[/math] に対して、[math]\displaystyle J_n(x)[/math] と [math]\displaystyle J_{-n}(x)[/math] とのあいだには関係
- [math]\displaystyle J_{-n}(x)=(-1)^nJ_n(x)[/math]
- が成り立ち、両者は明らかに線形従属となるからである。整数次数に対して [math]\displaystyle J_n(x)[/math] と線形独立な第二の解は、以下の第2種ベッセル関数によって与えられる。
- 第2種ベッセル関数
- 第2種ベッセル関数は [math]\displaystyle Y_\alpha(x)[/math] と表記される。[math]\displaystyle Y_\alpha(x)[/math] はやはりベッセルの微分方程式の解であり、[math]\displaystyle x=0[/math] において特異性をもつ。また、[math]\displaystyle Y_\alpha(x)[/math] は、しばしばノイマン関数とも呼ばれ、[math]\displaystyle N_\alpha(x)[/math] とも表記される。第1種ベッセル関数[math]\displaystyle J_\alpha(x)[/math] との関係は、
- [math]\displaystyle Y_\alpha(x)=\frac{J_\alpha(x)\cos(\alpha \pi)-J_{-\alpha}(x)}{\sin(\alpha \pi)}[/math]
- で与えられる。ただし、[math]\displaystyle \alpha[/math] が整数のときは右辺は極限によって定義されるものとする。
- 非整数の [math]\displaystyle \alpha[/math] に対しては、[math]\displaystyle J_\alpha(x)[/math] と[math]\displaystyle J_{-\alpha}(x)[/math] とが線形独立な2つの解をすでに与えているので、[math]\displaystyle Y_\alpha(x)[/math] は解の表現としては冗長である。整数 [math]\displaystyle n[/math] に対しては、[math]\displaystyle Y_n(x)[/math]は[math]\displaystyle J_n(x)[/math] と線形独立な第二の解を与えている。
- 整数 [math]\displaystyle n[/math] に対して、[math]\displaystyle Y_n(x)[/math] と [math]\displaystyle Y_{-n}(x)[/math] とのあいだには関係
- [math]\displaystyle Y_{-n}=(-1)^nY_n(x)[/math]
- が成り立ち、したがって両者はやはり線形従属である。
- [math]\displaystyle J_\alpha(x)[/math] 及び [math]\displaystyle Y_\alpha(x)[/math] はどちらも、負の実軸を除く複素平面上で [math]\displaystyle x[/math] の解析的な関数（正則な関数）である。[math]\displaystyle \alpha[/math] が正の整数のとき、これらの関数は負の実軸上に分岐点を持たず、したがって [math]\displaystyle x[/math] の整関数となる。また、固定した [math]\displaystyle x[/math] に対して、ベッセル関数は [math]\displaystyle \alpha[/math] の整関数となる。
- BesselJ plot.svg
第1種ベッセル関数
- BesselY plot.svg
第2種ベッセル関数
超幾何級数との関係
- ベッセル関数は超幾何級数（超幾何関数ともいう）によって、以下のように表現することができる。
- [math]J_\nu(z)=\frac{(z/2)^\nu}{\Gamma(\nu+1)} \;_0F_1 (\nu+1; -z^2/4).[/math]
ハンケル関数
- ベッセルの微分方程式に対する線形独立な2つの解を与える表式には、ハンケル関数Hα(1)(x) と Hα(2)(x)があり、定義式は以下の通り。
- [math]H_\alpha^{(1)}(x) = J_\alpha(x) + i Y_\alpha(x)[/math]
- [math]H_\alpha^{(2)}(x) = J_\alpha(x) - i Y_\alpha(x)[/math]
ここで、[math]i[/math]は虚数単位である。[math]J_\alpha(x)[/math]と[math]Y_\alpha(x)[/math]との線形結合によって与えられるこれらの解の表現は、第三種ベッセル関数として知られている(ハンケル関数は円筒波の方程式における、内側もしくは外側への円筒波の伝播の解を表現する)。
変形ベッセル関数
ベッセル関数は[math]\displaystyle x[/math]の複素数値に対しても適切に定義されており、応用上は [math]\displaystyle x[/math]が純虚数の場合が特に重要である。この場合、ベッセルの微分方程式への解は第1種及び第2種の変形ベッセル関数と呼ばれ、以下のように定義される。
- [math]\displaystyle I_\alpha(x) = i^{-\alpha} J_\alpha(ix) = \sum_{m=0}^{\infty} \frac{1}{m! \Gamma(m+\alpha+1)}\left(\frac{x}{2}\right)^{2m+\alpha}[/math]
- [math]K_\alpha(x) = \frac{\pi}{2} \frac{I_{-\alpha} (x) - I_\alpha (x)}{\sin (\alpha \pi)} = \frac{\pi}{2} i^{\alpha+1} H_\alpha^{(1)}(ix).[/math]
これらの関数は、[math]\displaystyle x[/math]が実数のときに関数値が実数となるように定義されている。またこれらの関数は、変形されたベッセルの微分方程式
- [math]x^2 \frac{d^2 y}{dx^2} + x \frac{dy}{dx} - (x^2 + \alpha^2)y = 0.[/math]
に対する2つの線形独立な解を与えている。
- BesselI plot.svg
第1種変形ベッセル関数
- BesselK plot.svg
第2種変形ベッセル関数
変形ベッセル関数には以下の性質がある。 ここで、n は正の整数またはゼロ。
- [math]\begin{align} &I_{-n}(x) = I_{n}(x) \\ &K_{-\alpha}(x) = K_{\alpha}(x) \\ \end{align}[/math]
- [math]\begin{align} &K_{n+1/2}(x) = \left(\frac{\pi}{2x}\right)^{1/2}\exp(-x)\sum^n_{r=0}{\frac{(n+r)!}{r!(n-r)!}(2x)^{-r}} \end{align}[/math]
球ベッセル関数・球ノイマン関数
第1種及び第2種のベッセル関数から、球ベッセル関数(spherical Bessel functions)と球ノイマン関数(spherical Neumann functions)がそれぞれ以下のように定義される。
- [math]j_\alpha(x) = \left( \frac{\pi}{2x} \right)^{1/2} J_{\alpha+1/2}(x)[/math]
- [math]n_\alpha(x) = \left( \frac{\pi}{2x} \right)^{1/2} N_{\alpha+1/2}(x)[/math]
これらの関数は、球ベッセル微分方程式
- [math]\frac{d^2y}{dx^2} + \frac{2}{x}\frac{dy}{dx} + \left[ 1 - \frac{\alpha(\alpha+1)}{x^2} \right] y = 0[/math]
に対する2つの線形独立な解を与えている。
量子力学における3次元自由粒子のシュレーディンガー方程式の動径方向の解のうち、正則なものは球ベッセル関数で表され、正則でないものは球ノイマン関数で表される。
また3次元井戸型ポテンシャルのシュレディンガー方程式における、ポテンシャル内部の動径方向の解のうち、原点で発散しないものは球ベッセル関数で表され、原点で発散するものは球ノイマン関数で表される。
球ハンケル関数
- 球ベッセル微分方程式に対する線形独立な2つの解を与える表式には、球ハンケル関数hα(1)(x) と hα(2)(x)があり、定義式は以下の通り。
- [math]h_\alpha^{(1)}(x) = j_\alpha(x) + i y_\alpha(x) = \sqrt{\frac{\pi}{2x}}H^{(1)}_{n+1/2}(x)[/math]
- [math]h_\alpha^{(2)}(x) = j_\alpha(x) - i y_\alpha(x) = \sqrt{\frac{\pi}{2x}}H^{(2)}_{n+1/2}(x)[/math]
ここで、[math]i[/math] は虚数単位である。
また、非負の整数 n について:
- [math]h_n^{(1)}(x) = (-i)^{n+1} \frac{e^{ix}}{x} \sum_{m=0}^n \frac{i^m}{m!(2x)^m} \frac{(n+m)!}{(n-m)!}[/math]
[math]h_n^{(2)}[/math] は、実数xに関して[math]h_n^{(1)}(x)[/math]の複素共役となる。
量子力学では、3次元井戸型ポテンシャルのシュレディンガー方程式における、ポテンシャル外部の動径方向の解は、球ハンケル関数で表される。第一種球ハンケル関数は外向き、第二種球ハンケル関数は内向きを表す。
変形球ベッセル関数
第1種及び第2種の変形ベッセル関数から、変形球ベッセル関数(英:modified spherical Bessel functions)が以下のように定義される。
- [math]i_{\alpha}(x) = \left( \frac{\pi}{2x} \right)^{1/2} I_{\alpha+1/2}(x)[/math]
- [math]k_{\alpha}(x) = \left( \frac{2}{\pi x} \right)^{1/2} K_{\alpha+1/2}(x)[/math]
これらの関数は、変形球ベッセル微分方程式
- [math]x^2\frac{d^2y}{dx^2} + 2x\frac{dy}{dx} - \left(x^2 + \alpha(\alpha+1) \right)y = 0[/math]
に対する2つの線形独立な解を与えている。
変形球ベッセル関数には以下の性質がある。
- [math]\begin{align} k_{-n-1/2}(x) = k_{n+1/2}(x) \end{align}[/math]
ここで、n は正の整数またはゼロ。
関連項目
参考文献
- Handbook of Mathematical Functions, Abramowitz and Stegun.
- Bessel Functions, Weisstein, Eric W. "Modified Bessel Functions" From MathWorld--A Wolfram Web Resource.
- A treatise on the theory of Bessel functions, George Neville Watson, Cambridge University Press,(1995).