ベッセル関数

ベッセル関数（ベッセルかんすう、英: Bessel function）とは、最初にスイスの数学者ダニエル・ベルヌーイによって定義され、フリードリヒ・ヴィルヘルム・ベッセルにちなんで名づけられた関数。円筒関数と呼ばれることもある。以下に示す、ベッセルの微分方程式における y ( x ) {\displaystyle y(x)} の特殊解の1つである。 x 2 d 2 y d x 2 + x d y d x + ( x 2 − α 2 ) y = 0 {\displaystyle x^{2}{\frac {d^{2}y}{dx^{2}}}+x{\frac {dy}{dx}}+(x^{2}-\alpha ^{2})y=0} 上の式において、 α {\displaystyle \alpha } は、任意の実数である（次数と呼ばれる）。 α {\displaystyle \alpha } が整数 n {\displaystyle n} に等しい場合がとくに重要である。 α {\displaystyle \alpha } 及び − α {\displaystyle -\alpha } はともに同一の微分方程式を与えるが、慣例としてこれら2つの異なる次数に対して異なるベッセル関数が定義される(例えば、 α {\displaystyle \alpha } の関数としてなるべく滑らかになるようにベッセル関数を定義する、など）。 そもそもベッセル関数は、惑星の軌道運動に関するケプラー方程式をベッセルが解析的に解いた際に導入された。 ベッセル解はラプラス方程式またはヘルムホルツ方程式の円柱座標系および極座標系における分離解として見出される。従ってベッセル関数は、電波伝播や静電位差などの解を求める際に重要である（円柱座標系においては整数次数 α = n {\displaystyle \alpha =n} 、極座標系においては半整数次数 α = n + 1 / 2 {\displaystyle \alpha =n+1/2} のベッセル関数がそれぞれ解として得られる）。例えば、 円筒導波管における電磁波 円柱物体の熱伝導 薄い円（か環状の）膜の振動のモード など。 ベッセル関数はまた、信号処理のような問題で有用な特性を持つ（例えば、FM合成、カイザー窓やベッセルフィルタなど）。 ベッセルの微分方程式は2階の線形微分方程式であるので、線形独立な2つの解が存在するはずである。しかしながら、解を議論する状況に応じて解の様々な表現が便利に使われている。代表的ないくつかの解の表現について以下で説明する。 これらの関数がベッセル関数群としては最も一般的な形式である。 第1種ベッセル関数 第1種ベッセル関数は J α ( x ) {\displaystyle J_{\alpha }(x)} と表記される。 J α ( x ) {\displaystyle J_{\alpha }(x)} はベッセルの微分方程式の解であり、 α {\displaystyle \alpha } が整数もしくは非負であるとき、 x = 0 {\displaystyle x=0} で有限の値をとる。 J α {\displaystyle J_{\alpha }} における特定解の選択及び正規化は後述する。第1種ベッセル関数はまた、 x = 0 {\displaystyle x=0} のまわりでのテイラー展開（非整数の α {\displaystyle \alpha } に対しては、より一般にべき級数展開）によって定義することもできる。 J α ( x ) = ∑ m = 0 ∞ ( − 1 ) m m ! Γ ( m + α + 1 ) ( x 2 ) 2 m + α {\displaystyle J_{\alpha }(x)=\sum _{m=0}^{\infty }{\frac {(-1)^{m}}{m!\Gamma (m+\alpha +1)}}\left({\frac {x}{2}}\right)^{2m+\alpha }} 非整数の α {\displaystyle \alpha } に対しては、 J α ( x ) {\displaystyle J_{\alpha }(x)} と J − α ( x ) {\displaystyle J_{-\alpha }(x)} とがベッセルの微分方程式に対する線形独立な2つの解を与える。他方で α {\displaystyle \alpha } が整数の場合には、 J − n ( x ) = ( − 1 ) n J n ( x ) {\displaystyle J_{-n}(x)=(-1)^{n}J_{n}(x)} という関係が成り立つため、2つの解は線形従属となる。整数次数に対して J n ( x ) {\displaystyle J_{n}(x)} と線形独立な第2の解は、第2種ベッセル関数によって与えられる。 第2種ベッセル関数 ノイマン関数 第2種ベッセル関数 Y α ( x ) {\displaystyle Y_{\alpha }(x)} はベッセルの微分方程式の解であり x = 0 {\displaystyle x=0} において特異性を持つ。ベッセル関数はノイマン関数とも呼ばれ、 N α ( x ) {\displaystyle N_{\alpha }(x)} と表される。 第2種ベッセル関数と第1種ベッセル関数 J α ( x ) {\displaystyle J_{\alpha }(x)} は以下の関係を持つ。 Y α ( x ) = J α ( x ) cos ⁡ ( α π ) − J − α ( x ) sin ⁡ ( α π ) {\displaystyle Y_{\alpha }(x)={\frac {J_{\alpha }(x)\cos(\alpha \pi )-J_{-\alpha }(x)}{\sin(\alpha \pi )}}} ただし、 α {\displaystyle \alpha } が整数のときは右辺は極限によって定義されるものとする。 非整数の α {\displaystyle \alpha } に対しては、 J α ( x ) {\displaystyle J_{\alpha }(x)} と J − α ( x ) {\displaystyle J_{-\alpha }(x)} とが線形独立な2つの解を既に与えているので、 Y α ( x ) {\displaystyle Y_{\alpha }(x)} は解の表現としては冗長である。整数 n {\displaystyle n} に対しては、 Y n ( x ) {\displaystyle Y_{n}(x)} は J n ( x ) {\displaystyle J_{n}(x)} と線形独立な第2の解を与えている。 整数 n {\displaystyle n} に対して、 Y n ( x ) {\displaystyle Y_{n}(x)} と Y − n ( x ) {\displaystyle Y_{-n}(x)} の間に Y − n = ( − 1 ) n Y n ( x ) {\displaystyle Y_{-n}=(-1)^{n}Y_{n}(x)} という関係が成り立ち、従って両者は線形従属である。 J α ( x ) {\displaystyle J_{\alpha }(x)} 及び Y α ( x ) {\displaystyle Y_{\alpha }(x)} はどちらも、負の実軸を除く複素平面上で x {\displaystyle x} の解析的な関数（正則な関数）である。 α {\displaystyle \alpha } が正の整数のとき、これらの関数は負の実軸上に分岐点を持たず、したがって x {\displaystyle x} の整関数となる。また、固定した x {\displaystyle x} に対して、ベッセル関数は α {\displaystyle \alpha } の整関数となる。 第1種ベッセル関数 第2種ベッセル関数 ベッセル関数は超幾何級数（超幾何関数ともいう）によって、以下のように表現することができる。 J ν ( z ) = ( z / 2 ) ν Γ ( ν + 1 ) 0 F 1 ( ν + 1 ; − z 2 / 4 ) . {\displaystyle J_{\nu }(z)={\frac {(z/2)^{\nu }}{\Gamma (\nu +1)}}\;_{0}F_{1}(\nu +1;-z^{2}/4).} ベッセルの微分方程式に対する線形独立な2つの解を与える表式には、ハンケル関数Hα(1)(x) と Hα(2)(x)があり、定義式は以下の通り。 H α ( 1 ) ( x ) = J α ( x ) + i Y α ( x ) {\displaystyle H_{\alpha }^{(1)}(x)=J_{\alpha }(x)+iY_{\alpha }(x)} H α ( 2 ) ( x ) = J α ( x ) − i Y α ( x ) {\displaystyle H_{\alpha }^{(2)}(x)=J_{\alpha }(x)-iY_{\alpha }(x)} ここで、 i {\displaystyle i} は虚数単位である。 J α ( x ) {\displaystyle J_{\alpha }(x)} と Y α ( x ) {\displaystyle Y_{\alpha }(x)} との線形結合によって与えられるこれらの解の表現は、第三種ベッセル関数として知られている(ハンケル関数は円筒波の方程式における、内側もしくは外側への円筒波の伝播の解を表現する)。 ベッセル関数は x {\displaystyle \displaystyle x} の複素数値に対しても適切に定義されており、応用上は x {\displaystyle \displaystyle x} が純虚数の場合が特に重要である。この場合、ベッセルの微分方程式への解は第1種及び第2種の変形ベッセル関数と呼ばれ、以下のように定義される。 I α ( x ) = i − α J α ( i x ) = ∑ m = 0 ∞ 1 m ! Γ ( m + α + 1 ) ( x 2 ) 2 m + α {\displaystyle \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 }} K α ( x ) = π 2 I − α ( x ) − I α ( x ) sin ⁡ ( α π ) = π 2 i α + 1 H α ( 1 ) ( i x ) . {\displaystyle 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).} これらの関数は、 x {\displaystyle \displaystyle x} が実数のときに関数値が実数となるように定義されている。またこれらの関数は、変形されたベッセルの微分方程式 x 2 d 2 y d x 2 + x d y d x − ( x 2 + α 2 ) y = 0. {\displaystyle x^{2}{\frac {d^{2}y}{dx^{2}}}+x{\frac {dy}{dx}}-(x^{2}+\alpha ^{2})y=0.} に対する2つの線形独立な解を与えている。 第1種変形ベッセル関数 第2種変形ベッセル関数 変形ベッセル関数には以下の性質がある。 ここで、n は正の整数またはゼロ。 I − n ( x ) = I n ( x ) K − α ( x ) = K α ( x ) {\displaystyle {\begin{aligned}&I_{-n}(x)=I_{n}(x)\\&K_{-\alpha }(x)=K_{\alpha }(x)\\\end{aligned}}} K n + 1 / 2 ( x ) = ( π 2 x ) 1 / 2 exp ⁡ ( − x ) ∑ r = 0 n ( n + r ) ! r ! ( n − r ) ! ( 2 x ) − r {\displaystyle {\begin{aligned}&K_{n+1/2}(x)=\left({\frac {\pi }{2x}}\right)^{1/2}\exp(-x)\sum _{r=0}^{n}{{\frac {(n+r)!}{r!(n-r)!}}(2x)^{-r}}\end{aligned}}} 第1種及び第2種のベッセル関数から、球ベッセル関数(spherical Bessel functions)と球ノイマン関数(spherical Neumann functions)がそれぞれ以下のように定義される。 j α ( x ) = ( π 2 x ) 1 / 2 J α + 1 / 2 ( x ) {\displaystyle j_{\alpha }(x)=\left({\frac {\pi }{2x}}\right)^{1/2}J_{\alpha +1/2}(x)} n α ( x ) = ( π 2 x ) 1 / 2 N α + 1 / 2 ( x ) {\displaystyle n_{\alpha }(x)=\left({\frac {\pi }{2x}}\right)^{1/2}N_{\alpha +1/2}(x)} これらの関数は、球ベッセル微分方程式 d 2 y d x 2 + 2 x d y d x + [ 1 − α ( α + 1 ) x 2 ] y = 0 {\displaystyle {\frac {d^{2}y}{dx^{2}}}+{\frac {2}{x}}{\frac {dy}{dx}}+\left[1-{\frac {\alpha (\alpha +1)}{x^{2}}}\right]y=0} に対する2つの線形独立な解を与えている。 量子力学における3次元自由粒子のシュレーディンガー方程式の動径方向の解のうち、正則なものは球ベッセル関数で表され、正則でないものは球ノイマン関数で表される。 また3次元井戸型ポテンシャルのシュレディンガー方程式における、ポテンシャル内部の動径方向の解のうち、原点で発散しないものは球ベッセル関数で表され、原点で発散するものは球ノイマン関数で表される。 球ベッセル微分方程式に対する線形独立な2つの解を与える表式には、球ハンケル関数hα(1)(x) と hα(2)(x)があり、定義式は以下の通り。 h α ( 1 ) ( x ) = j α ( x ) + i y α ( x ) = π 2 x H n + 1 / 2 ( 1 ) ( x ) {\displaystyle h_{\alpha }^{(1)}(x)=j_{\alpha }(x)+iy_{\alpha }(x)={\sqrt {\frac {\pi }{2x}}}H_{n+1/2}^{(1)}(x)} h α ( 2 ) ( x ) = j α ( x ) − i y α ( x ) = π 2 x H n + 1 / 2 ( 2 ) ( x ) {\displaystyle h_{\alpha }^{(2)}(x)=j_{\alpha }(x)-iy_{\alpha }(x)={\sqrt {\frac {\pi }{2x}}}H_{n+1/2}^{(2)}(x)} ここで、 i {\displaystyle i} は虚数単位である。 また、非負の整数 n について: h n ( 1 ) ( x ) = ( − i ) n + 1 e i x x ∑ m = 0 n i m m ! ( 2 x ) m ( n + m ) ! ( n − m ) ! {\displaystyle 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)!}}} h n ( 2 ) {\displaystyle h_{n}^{(2)}} は、実数xに関して h n ( 1 ) ( x ) {\displaystyle h_{n}^{(1)}(x)} の複素共役となる。 量子力学では、3次元井戸型ポテンシャルのシュレディンガー方程式における、ポテンシャル外部の動径方向の解は、球ハンケル関数で表される。第一種球ハンケル関数は外向き、第二種球ハンケル関数は内向きを表す。 第1種及び第2種の変形ベッセル関数から、変形球ベッセル関数(英:modified spherical Bessel functions)が以下のように定義される。 i α ( x ) = ( π 2 x ) 1 / 2 I α + 1 / 2 ( x ) {\displaystyle i_{\alpha }(x)=\left({\frac {\pi }{2x}}\right)^{1/2}I_{\alpha +1/2}(x)} k α ( x ) = ( 2 π x ) 1 / 2 K α + 1 / 2 ( x ) {\displaystyle k_{\alpha }(x)=\left({\frac {2}{\pi x}}\right)^{1/2}K_{\alpha +1/2}(x)} これらの関数は、変形球ベッセル微分方程式 x 2 d 2 y d x 2 + 2 x d y d x − ( x 2 + α ( α + 1 ) ) y = 0 {\displaystyle x^{2}{\frac {d^{2}y}{dx^{2}}}+2x{\frac {dy}{dx}}-\left(x^{2}+\alpha (\alpha +1)\right)y=0} に対する2つの線形独立な解を与えている。 変形球ベッセル関数には以下の性質がある。 k − n − 1 / 2 ( x ) = k n + 1 / 2 ( x ) {\displaystyle {\begin{aligned}k_{-n-1/2}(x)=k_{n+1/2}(x)\end{aligned}}} ここで、n は正の整数またはゼロ。 Bessel の積分表示 J n ( z ) = 1 π ∫ 0 π cos ⁡ ( n θ − z sin ⁡ θ ) d θ = 1 2 π ∫ 0 2 π cos ⁡ ( n θ − z sin ⁡ θ ) d θ {\displaystyle J_{n}(z)={\frac {1}{\pi }}\int _{0}^{\pi }\cos(n\theta -z\sin \theta )d\theta ={\frac {1}{2\pi }}\int _{0}^{2\pi }\cos(n\theta -z\sin \theta )d\theta } Hansen の積分表示 J n ( z ) = 1 π i n ∫ 0 π e i z cos ⁡ θ cos ⁡ n θ d θ {\displaystyle J_{n}(z)={\frac {1}{\pi i^{n}}}\int _{0}^{\pi }e^{iz\cos \theta }\cos n\theta d\theta } Poisson の積分表示 J n ( z ) = ( z / 2 ) n π Γ ( n + 1 / 2 ) ∫ 0 π cos ⁡ ( z cos ⁡ θ ) sin 2 n ⁡ θ d θ {\displaystyle J_{n}(z)={\frac {(z/2)^{n}}{{\sqrt {\pi }}\Gamma (n+1/2)}}\int _{0}^{\pi }\cos(z\cos \theta )\sin ^{2n}\theta d\theta } Schläfli の積分表示 N ν ( z ) = 1 π ∫ 0 π sin ⁡ ( z sin ⁡ θ − ν θ ) d θ − 1 π ∫ 0 ∞ e − z sinh ⁡ t ( e ν t + e − ν t cos ⁡ ν π ) d t     ( R e ( z ) > 0 ) {\displaystyle N_{\nu }(z)={\frac {1}{\pi }}\int _{0}^{\pi }\sin(z\sin \theta -\nu \theta )d\theta -{\frac {1}{\pi }}\int _{0}^{\infty }e^{-z\sinh t}(e^{\nu t}+e^{-\nu t\cos \nu \pi })dt\ \ (\mathrm {Re} (z)>0)} Schafheitlin の積分表示 ただし複号は上が ι = 1 {\displaystyle \iota =1} , 下が ι = 2 {\displaystyle \iota =2} π Γ ( ν + 1 / 2 ) z ν H ν ( ι ) ( z ) = ∓ 2 ν + 1 i ∫ 0 π / 2 exp ⁡ { ± i [ z − ( ν − 1 2 ) θ ] − 2 z cot ⁡ θ } cos ν − 1 / 2 ⁡ θ c o s e c 2 ν + 1 θ d θ     [ R e ( ν ) > − 1 2 ,   R e ( z ) > 0 ] {\displaystyle {\frac {{\sqrt {\pi }}\Gamma (\nu +1/2)}{z^{\nu }}}H_{\nu }^{(\iota )}(z)=\mp 2^{\nu +1}i\int _{0}^{\pi /2}\exp \left\{\pm i\left[z-\left(\nu -{\frac {1}{2}}\right)\theta \right]-2z\cot \theta \right\}\,\cos ^{\nu -1/2}\theta \,\mathrm {cosec} ^{2\nu +1}\theta \,d\theta \ \ [\mathrm {Re} (\nu )>-{\frac {1}{2}},\ \mathrm {Re} (z)>0]} Heine の積分表示 ただし複号は上が ι = 1 {\displaystyle \iota =1} , 下が ι = 2 {\displaystyle \iota =2} H ν ( ι ) ( z ) = ∓ 2 i e ∓ ν π i / 2 π ∫ 0 ∞ e ± i z cosh ⁡ t cosh ⁡ ν t d t     [ | R e ( ν ) | < 1 , ±   I m ( z ) > 0 ] {\displaystyle H_{\nu }^{(\iota )}(z)={\frac {\mp 2ie^{\mp \nu \pi i/2}}{\pi }}\int _{0}^{\infty }e^{\pm iz\cosh t}\cosh \nu t\,dt\ \ [|\mathrm {Re} (\nu )|<1,\pm \ \mathrm {Im} (z)>0]} Whittaker の積分表示 ここに P n ( t ) {\displaystyle P_{n}(t)} はルジャンドル多項式 j n ( z ) = 1 2 i n ∫ − 1 1 e i z t P n ( t ) d t {\displaystyle j_{n}(z)={\frac {1}{2i^{n}}}\int _{-1}^{1}e^{izt}P_{n}(t)dt} | z | → ∞ {\displaystyle |z|\to \infty } のときベッセル関数は以下の漸近形を持つ。 J ν ( z ) ∼ 2 π z cos ⁡ ( z − 2 ν + 1 4 π ) {\displaystyle J_{\nu }(z)\sim {\sqrt {\frac {2}{\pi z}}}\cos \left(z-{\frac {2\nu +1}{4}}\pi \right)} N ν ( z ) ∼ 2 π z sin ⁡ ( z − 2 ν + 1 4 π ) {\displaystyle N_{\nu }(z)\sim {\sqrt {\frac {2}{\pi z}}}\sin \left(z-{\frac {2\nu +1}{4}}\pi \right)} H ν ( 1 ) ( z ) ∼ 2 π z exp ⁡ { i ( z − 2 ν + 1 4 π ) } {\displaystyle H_{\nu }^{(1)}(z)\sim {\sqrt {\frac {2}{\pi z}}}\exp \left\{i\left(z-{\frac {2\nu +1}{4}}\pi \right)\right\}} H ν ( 2 ) ( z ) ∼ 2 π z exp ⁡ { − i ( z − 2 ν + 1 4 π ) } {\displaystyle H_{\nu }^{(2)}(z)\sim {\sqrt {\frac {2}{\pi z}}}\exp \left\{-i\left(z-{\frac {2\nu +1}{4}}\pi \right)\right\}} j n ( z ) ∼ 1 z cos ⁡ ( z − n + 1 2 π ) {\displaystyle j_{n}(z)\sim {\frac {1}{z}}\cos \left(z-{\frac {n+1}{2}}\pi \right)} n n ( z ) ∼ 1 z sin ⁡ ( z − n + 1 2 π ) {\displaystyle n_{n}(z)\sim {\frac {1}{z}}\sin \left(z-{\frac {n+1}{2}}\pi \right)} h n ( 1 ) ( z ) ∼ ( − i ) n + 1 z e i z {\displaystyle h_{n}^{(1)}(z)\sim {\frac {(-i)^{n+1}}{z}}e^{iz}} h n ( 2 ) ( z ) ∼ i n + 1 z e − i z {\displaystyle h_{n}^{(2)}(z)\sim {\frac {i^{n+1}}{z}}e^{-iz}} [脚注の使い方] ^ “Bessel function”. Britannica. 2021年3月20日閲覧。 ^ a b c 岩波数学公式, p. 178. ^ 岩波数学公式, p. 182. ^ a b 岩波数学公式, p. 183. ^ 岩波数学公式, p. 185. ^ 岩波数学公式, pp. 154, 168. 特殊関数 Jackson の q-Bessel 関数（英語版） Hahn-Exton の q-Bessel 関数（英語版） 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). 森口繁一、宇田川銈久、一松信『岩波数学公式III 特殊函数』岩波書店、1987年。ISBN 4-00-005509-7。  ベッセル関数の計算 Bessel 関数の数値計算 古典的Bessel関数入門