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

'''ベッセル関数'''（ベッセルかんすう、{{lang-en-short|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>の関数としてなるべく滑らかになるようにベッセル関数を定義する、など）。

そもそもベッセル関数は、[[惑星]]の[[軌道 (力学)|軌道]]運動に関する[[ケプラー方程式]]をベッセルが解析的に解いた際に導入された<ref>{{Cite web |url=https://www.britannica.com/science/Bessel-function |title=Bessel function |publisher=Britannica |date= |accessdate=2021-03-20}}</ref>。

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

ベッセル関数はまた、[[信号処理]]のような問題で有用な特性を持つ（例えば、[[FM合成]]、[[カイザー窓]]や[[ベッセルフィルタ]]など）。

== 定義 ==
ベッセルの微分方程式は2階の線形微分方程式であるので、線形独立な2つの解が存在するはずである。しかしながら、解を議論する状況に応じて解の様々な表現が便利に使われている。代表的ないくつかの解の表現について以下で説明する。

=== 第1種及び第2種ベッセル関数 ===

これらの関数がベッセル関数群としては最も一般的な形式である。
; 第1種ベッセル関数{{anchors|第1種ベッセル関数}} 
: 第1種ベッセル関数は <math>J_\alpha(x)</math> と表記される。<math>J_\alpha(x)</math> はベッセルの微分方程式の解であり、<math>\alpha</math> が整数もしくは非負であるとき、<math>x=0</math> で有限の値をとる。<math>J_\alpha</math> における特定解の選択及び正規化は後述する。第1種ベッセル関数はまた、<math>x=0</math> のまわりでの[[テイラー展開]]（非整数の <math>\alpha</math> に対しては、より一般にべき級数展開）によって定義することもできる。
:<math>J_\alpha(x)=\sum_{m=0}^\infty \frac{(-1)^m}{m!\Gamma(m+\alpha+1)} \left ( \frac{x}{2} \right)^{2m+\alpha}</math>

非整数の <math>\alpha</math> に対しては、<math>J_\alpha(x)</math> と <math>J_{-\alpha}(x)</math> とがベッセルの微分方程式に対する線形独立な2つの解を与える。他方で<math>\alpha</math> が整数の場合には、<math>J_{-n}(x)=(-1)^nJ_n(x)</math> という関係が成り立つため、2つの解は線形従属となる。整数次数に対して <math>J_n(x)</math> と線形独立な第2の解は、第2種ベッセル関数によって与えられる。

; 第2種ベッセル関数{{anchors|第2種ベッセル関数}}
; ノイマン関数{{anchors|ノイマン関数}}
: 第2種ベッセル関数 <math>Y_\alpha(x)</math> はベッセルの微分方程式の解であり <math>x=0</math> において特異性を持つ。ベッセル関数はノイマン関数とも呼ばれ、<math>N_\alpha(x)</math> と表される。
:第2種ベッセル関数と第1種ベッセル関数 <math>J_\alpha(x)</math> は以下の関係を持つ。
::<math>Y_\alpha(x)=\frac{J_\alpha(x)\cos(\alpha \pi)-J_{-\alpha}(x)}{\sin(\alpha \pi)}</math>
: ただし、<math>\alpha</math> が整数のときは右辺は極限によって定義されるものとする。

非整数の <math>\alpha</math> に対しては、<math>J_\alpha(x)</math> と <math>J_{-\alpha}(x)</math> とが線形独立な2つの解を既に与えているので、<math>Y_\alpha(x)</math> は解の表現としては冗長である。整数 <math>n</math> に対しては、<math>Y_n(x)</math> は <math>J_n(x)</math> と線形独立な第2の解を与えている。
整数 <math>n</math> に対して、<math>Y_n(x)</math> と <math>Y_{-n}(x)</math> の間に <math>Y_{-n}=(-1)^nY_n(x)</math> という関係が成り立ち、従って両者は線形従属である。

<math>J_\alpha(x)</math> 及び <math>Y_\alpha(x)</math> はどちらも、負の実軸を除く複素平面上で <math>x</math> の解析的な関数（正則な関数）である。<math>\alpha</math> が正の整数のとき、これらの関数は負の実軸上に[[分岐点_(数学)|分岐点]]を持たず、したがって <math>x</math> の[[整関数]]となる。また、固定した <math>x</math> に対して、ベッセル関数は <math>\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 は正の整数またはゼロ。

== 積分表示 ==
Bessel の積分表示<ref name="#1">岩波数学公式, p. 178.</ref>
{{Indent|<math>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</math>}}
Hansen の積分表示<ref name="#1"/>
{{Indent|<math>J_n ( z ) = \frac{ 1 }{ \pi i^n } \int_0^\pi e^{i z \cos \theta} \cos n \theta d\theta</math>}}
Poisson の積分表示<ref name="#1"/>
{{Indent|<math>J_n ( z ) = \frac{ ( z / 2 )^n }{ \sqrt{ \pi } \Gamma ( n + 1/2 ) } \int_0^\pi \cos ( z \cos \theta ) \sin^{2 n} \theta d\theta</math>}}
Schläfli の積分表示<ref>岩波数学公式, p. 182.</ref>
{{Indent|<math>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 )</math>}}
Schafheitlin の積分表示<ref name="#2">岩波数学公式, p. 183.</ref> ただし複号は上が <math>\iota = 1</math>, 下が <math>\iota = 2</math>
{{Indent|<math>\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] - 2 z \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]</math>}}
Heine の積分表示<ref name="#2"/> ただし複号は上が <math>\iota = 1</math>, 下が <math>\iota = 2</math>
{{Indent|<math>H_\nu^{(\iota)} ( z ) = \frac{ \mp 2 i e^{\mp \nu \pi i / 2} }{ \pi } \int_0^\infty e^{\pm i z \cosh t} \cosh \nu t \, dt \ \ [ |\mathrm{Re}(\nu)| < 1 , \pm \ \mathrm{Im}(z) > 0]</math>}}
Whittaker の積分表示<ref>岩波数学公式, p. 185.</ref> ここに <math>P_n ( t )</math> は[[ルジャンドル多項式]]
{{Indent|<math>j_n ( z ) = \frac{ 1 }{ 2 i^n } \int_{-1}^1 e^{i z t} P_n ( t ) dt</math>}}

== 漸近展開 ==
<math>| z | \to \infty</math> のときベッセル関数は以下の漸近形を持つ<ref>岩波数学公式, pp. 154, 168.</ref>。
{{Indent|<math>J_\nu ( z ) \sim \sqrt{ \frac{ 2 }{ \pi z } } \cos \left( z - \frac{ 2 \nu + 1 }{ 4 } \pi \right)</math>}}
{{Indent|<math>N_\nu ( z ) \sim \sqrt{ \frac{ 2 }{ \pi z } } \sin \left( z - \frac{ 2 \nu + 1 }{ 4 } \pi \right)</math>}}
{{Indent|<math>H_\nu^{(1)} ( z ) \sim \sqrt{ \frac{ 2 }{ \pi z } } \exp \left\{ i \left( z - \frac{ 2 \nu + 1 }{ 4 } \pi \right) \right\}</math>}}
{{Indent|<math>H_\nu^{(2)} ( z ) \sim \sqrt{ \frac{ 2 }{ \pi z } } \exp \left\{ - i \left( z - \frac{ 2 \nu + 1 }{ 4 } \pi \right) \right\}</math>}}
{{Indent|<math>j_n ( z ) \sim \frac{ 1 }{ z } \cos \left( z - \frac{ n + 1 }{ 2 } \pi \right)</math>}}
{{Indent|<math>n_n ( z ) \sim \frac{ 1 }{ z } \sin \left( z - \frac{ n + 1 }{ 2 } \pi \right)</math>}}
{{Indent|<math>h_n^{(1)} ( z ) \sim \frac{ ( - i )^{n + 1} }{ z } e^{i z}</math>}}
{{Indent|<math>h_n^{(2)} ( z ) \sim \frac{ i^{n + 1} }{ z } e^{- i z}</math>}}

== 脚注 ==
{{脚注ヘルプ}}
=== 出典 ===
{{Reflist|2}}

== 関連項目 ==
* [[特殊関数]]
* {{ill2|Jackson の q-Bessel 関数|en|Jackson q-Bessel function|label=Jackson の {{mvar|q}}-Bessel 関数}}
* {{ill2|Hahn–Exton の q-Bessel 関数|en|Hahn–Exton q-Bessel function|label=Hahn-Exton の {{mvar|q}}-Bessel 関数}}

==参考文献==
* [https://personal.math.ubc.ca/~cbm/aands/ 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. 
* [https://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).
* {{Cite book |和書 |author1=森口繁一 |author2=宇田川銈久 |author3=一松信 |title=岩波数学公式III 特殊函数 |publisher=岩波書店 |date=1987 |isbn=4-00-005509-7}}

== 外部リンク ==
* [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 ベッセル関数の計算]
* [http://nalab.mind.meiji.ac.jp/~mk/labo/text/computing-bessel-function/ Bessel 関数の数値計算]
* [http://math.univ-lyon1.fr/~iohara/Doc/Jap/Bessel.pdf 古典的Bessel関数入門] 

{{Normdaten}}

{{デフォルトソート:へつせるかんすう}}
[[Category:特殊関数]]
[[Category:フリードリヒ・ヴィルヘルム・ベッセル]]
[[Category:数学に関する記事]]
[[Category:数学のエポニム]]