ベッセル関数のソースを表示

'''ベッセル関数'''（ベッセルかんすう、''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>が整数nに等しい場合がとくに重要である。

<math>\alpha</math>及び<math>-\alpha</math>はともに同一の微分方程式を与えるが、慣例としてこれら2つの異なる次数に対して異なるベッセル関数が定義される(例えば、<math>\alpha</math>の関数としてなるべく滑らかになるようにベッセル関数を定義する、など）。

そもそもベッセル関数は、惑星軌道の時間変化に関するケプラー方程式を、ベッセルが解析的に解いた際に導入された。

==応用==
ベッセル解は、[[ラプラス方程式]]または、[[ヘルムホルツ方程式]]の、円柱座標系または極座標系における分離解として見出される。したがって、ベッセル関数は、[[電波伝播]]や[[静電位差]]などの解を求めるときに重要なものとなっている。（円柱座標系においては、整数次数<math>\alpha=n</math>のベッセル関数が解として得られる。極座標系においては、半整数次数<math>\alpha=n+1/2</math>のベッセル関数が解として得られる。）</br>
例えば、
*円筒導波管における[[電磁波]]
*円柱物体の[[熱伝導]]
*薄い円(か環状の) 膜の振動のモード

ベッセル関数には、また信号処理のような他の問題のための有用な特性がある(例えば、[[FM合成]]、[[Kaiser窓]]や[[ベッセルフィルタ]]など)。

==定義式==
ベッセルの微分方程式は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> の整関数となる。

<gallery>
画像:BesselJ plot.svg|第1種ベッセル関数 
画像:BesselY plot.svg|第2種ベッセル関数
</gallery>

===超幾何級数との関係===
*ベッセル関数は[[超幾何級数]]（超幾何関数ともいう）によって、以下のように表現することができる。
:<math>J_\nu(z)=\frac{(z/2)^\nu}{\Gamma(\nu+1)}  \;_0F_1 (\nu+1; -z^2/4).</math>

===ハンケル関数===
*ベッセルの微分方程式に対する線形独立な2つの解を与える表式には、ハンケル関数''H''<sub>&alpha;</sub><sup>(1)</sup>(''x'') と ''H''<sub>&alpha;</sub><sup>(2)</sup>(''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つの線形独立な解を与えている。

<gallery>
Image:BesselI_plot.svg|第1種変形ベッセル関数
Image:BesselK_plot.svg|第2種変形ベッセル関数
</gallery>

変形ベッセル関数には以下の性質がある。
ここで、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''<sub>&alpha;</sub><sup>(1)</sup>(''x'') と ''h''<sub>&alpha;</sub><sup>(2)</sup>(''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 は正の整数またはゼロ。

== 関連項目 ==
* [[特殊関数]]
==参考文献==
* [http://www.math.sfu.ca/~cbm/aands/index.htm Handbook of Mathematical Functions], [[Abramowitz and Stegun]].
* [http://mathworld.wolfram.com/topics/BesselFunctions.html Bessel Functions], Weisstein, Eric W. "Modified Bessel Functions" From MathWorld--A Wolfram Web Resource. 
* [http://books.google.co.jp/books?id=Mlk3FrNoEVoC&lpg=PP1&pg=PP1#v=onepage&q=&f=false A treatise on the theory of Bessel functions], George Neville Watson, Cambridge University Press,(1995).

== 外部リンク ==
* [http://keisan.casio.jp/has10/Menu.cgi?path=08000000%2e%93%c1%8e%ea%8a%d6%90%94%2f04000100%2e%83x%83b%83Z%83%8b%8a%d6%90%94 ベッセル関数の計算]
{{デフォルトソート:へつせるかんすう}}
[[Category:特殊関数]]
[[Category:数学に関する記事]]