贝索函数

王朝百科·作者佚名 2010-03-07

贝塞尔函数（Bessel functions）是数学上的一类特殊函数的总称。一般贝塞尔函数是下列常微分方程（一般称为'''贝塞尔方程'''）的标准解函数<math>y(x)</math>：

:<math>x^2 frac{d^2 y}{dx^2} + x frac{dy}{dx} + (x^2 - alpha^2)y = 0</math>

这类方程的解是无法用初等函数系统地表示的。

贝塞尔函数的具体形式随上述方程中任意实数α变化而变化（相应地，α被称为其对应贝塞尔函数的阶数）。实际应用中最常见的情形为α是整数''n''，对应解称为'''''n'' 阶贝塞尔函数'''。

尽管在上述微分方程中，α本身的正负号不改变方程的形式，但实际应用中仍习惯针对α和−α定义两种不同的贝塞尔函数（这样做能带来好处，比如消除了函数在α=0 点的不光滑性）。

==历史==

贝塞尔函数的几个正整数阶特例早在18世纪中叶就由瑞士数学家丹尼尔·伯努利在研究悬链振动时提出了，当时引起了数学界的兴趣。丹尼尔·伯努利|丹尼尔的叔叔雅各布·伯努利，莱昂哈德·欧拉|欧拉、约瑟夫·路易斯·拉格朗日|拉格朗日等数学大师对贝塞尔函数的研究作出过重要贡献。1817年，德国数学家弗里德里希·威廉·贝塞尔|贝塞尔在研究约翰内斯·开普勒|开普勒提出的三体万有引力|引力系统的运动问题时，第一次系统地提出了贝塞尔函数的总体理论框架，后人以他的名字来命名了这种函数 [http://www.britannica.com/eb/article-9078932] [http://www-history.mcs.st-andrews.ac.uk/Biographies/Bessel.html]。

==现实背景和应用范围==

贝塞尔方程是在柱坐标或球坐标下使用分离变量法求解拉普拉斯方程和亥姆霍兹方程时得到的（在圆柱域问题中得到的是整阶形式 α = ''n''；在球形域问题中得到的是半奇数阶形式 α = ''n''+½），因此贝塞尔函数在波动问题以及各种涉及有势场的问题中占有非常重要的地位，最典型的问题有：

来自“http://www.zhdbk.com”中华大百科

* 在圆柱形波导中的电磁波传播问题；

* 圆柱体中的热传导定律|热传导问题；

* 圆形（或环形）薄膜的振动模态分析问题；

在其他一些领域，贝塞尔函数也相当有用。譬如在信号处理中的调频合成（w:Frequency modulation synthesis|FM synthesis）或凯泽窗（w:Kaiser window|Kaiser window）的定义中，都要用到贝塞尔函数。

==定义==

贝塞尔方程是一个二阶常微分方程，必然存在两个矢量|线性无关的解。针对各种具体情况，人们提出了表示这些解的不同形式。下面分别介绍这些不同类型的贝塞尔函数。

===第一类贝塞尔函数===

image:BesselJ_plot.svg|right|thumb|400px|'''图2''' 0阶、1阶和2阶第一类贝塞尔函数（贝塞尔J函数）曲线

（''在下文中，第一类贝塞尔函数有时会简称为“J函数”，敬请读者留意。''）

第一类α阶贝塞尔函数''J''α(''x'')是贝塞尔方程当α为整数或α非负时的解，须满足在''x'' = 0 时有限。这样选取和处理''J''α的原因见本主题下面的贝塞尔函数#性质|性质介绍；另一种定义方法是通过它在''x'' = 0 点的泰勒级数展开（或者更一般地通过幂级数展开，这适用于α为非整数）：

:<math> J_alpha(x) = sum_{m=0}^infty frac{(-1)^m}{m! Gamma(m+alpha+1)} {left({frac{x}{2}}

ight)}^{2m+alpha} </math>

上式中<math>Gamma(z)</math>为Γ函数（它可视为阶乘|阶乘函数向非整型因变量和自变量|自变量的推广）。第一类贝塞尔函数的形状大致与按<math>1/sqrt x </math>速率衰减的正弦或三角函数|余弦函数类似（参见本页下面对它们渐进形式的介绍），但它们的零点并不是周期性的，另外随着''x''的增加，零点的间隔会越来越接近周期性。图2所示为0阶、1阶和2阶第一类贝塞尔函数<math>J_alpha (x)</math>的曲线（<math>alpha = 0, 1, 2</math>）。

如果α不为整数，则<math>J_alpha (x)</math>和<math>J_{-alpha} (x)</math>线性无关，可以构成微分方程的一个'''解系'''。反之若<math>alpha</math>是整数，那么上面两个函数之间满足如下关系：

:<math>J_{-alpha}(x) = (-1)^{alpha} J_{alpha}(x),</math>

于是两函数之间已不满足线性无关条件。为寻找在此情况下微分方程与<math>J_alpha (x)</math>线性无关的另一解，需要定义'''第二类贝塞尔函数'''，定义过程将在后面的小节中给出。

====贝塞尔积分====

<math>alpha</math>为整数时贝塞尔函数的另一种定义方法由下面的积分给出：

:<math>J_alpha (x) = frac{1}{2 pi} int_{0}^{2 pi} cos (alpha au - x sin au) dau.</math>

（<math>alpha</math>为任意实数时的表达式见贝塞尔函数#参考文献|参考文献[2]第360页）

这个积分式就是贝塞尔当年提出的定义，而且他还从该定义中推出了函数的一些性质。另一种积分表达式为：

:<math>J_alpha (x) = frac{1}{2 pi} int_{-pi}^{pi} e^{i(alpha au - x sin au)} dau</math>

====和超几何级数的关系====

贝塞尔函数可以用超几何级数表示成下面的形式：

:<math>J_alpha(z)=frac{(z/2)^alpha}{Gamma(alpha+1)} ;_0F_1 (alpha+1; -z^2/4).</math>

===第二类贝塞尔函数（诺依曼函数）===

image:BesselY_plot.svg|right|thumb|400px|'''图3''' 0阶、1阶和2阶第二类贝塞尔函数（贝塞尔''Y'' 函数）曲线图

（''在下文中，第二类贝塞尔函数有时会简称为“Y函数”，敬请读者留意。''）

'''第二类贝塞尔函数'''也许比第一类更为常用。

这种函数通常用''Y''α(''x'')表示，它们是贝塞尔方程的另一类解。''x'' = 0 点是第二类贝塞尔函数的（无穷）奇点。

''Y''α(''x'')又被称为'''诺依曼函数'''（Neumann function），有时也记作''N''α(''x'')。它和''J''α(''x'')存在如下关系：

:<math>Y_alpha(x) = frac{J_alpha(x) cos(alphapi) - J_{-alpha}(x)}{sin(alphapi)},</math>

若α为整数（此时上式是0/0型未定式）则取右端的极限值。

从前面对''J''α(''x'')的定义可以知道，若α不为整数时，定义''Y''α是多余的（因为贝塞尔方程的两个线性无关解都已经用J函数表示出来了）。另一方面，若α为整数，''Y''''α''便可以和''J''''α''构成贝塞尔方程的一个解系。与J函数类似，Y函数正负整数阶之间也存在如下关系：

:<math>Y_{-n}(x) = (-1)^n Y_n(x),</math>

''J''α(''x'')和''Y''α(''x'')均为沿负实半轴割开的复数 (数学)#.E5.A4.8D.E5.B9.B3.E9.9D.A2|复平面内关于''x''的全纯函数。当α为整数时，复平面内不存在贝塞尔函数的支点，所以''J'' 和''Y'' 均为''x'' 的整函数。若将''x'' 固定，则贝塞尔函数是α的整函数。图3所示为0阶、1阶和2阶第二类贝塞尔函数<math>Y_alpha (x)</math>的曲线（<math>alpha = 0, 1, 2</math>）：

===汉开尔函数===

贝塞尔方程的另外一对重要的线性无关解称为'''赫尔曼·汉开尔|汉开尔函数'''（Hankel functions）''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>

其中''i'' 为虚数单位<math>sqrt { - 1}</math>。以上的线性组合也成为'''第三类贝塞尔函数'''；它们描述了二维波动方程的内行柱面波解和外行柱面波解（"行"与在"行动"中同音）。

利用前面推出的关系可将汉开尔函数表示成：

:<math>H_{alpha}^{(1)} (x) = frac{J_{-alpha} (x) - e^{-alpha pi i} J_alpha (x)}{i sin (alpha pi)}</math>

:<math>H_{alpha}^{(2)} (x) = frac{J_{-alpha} (x) - e^{alpha pi i} J_alpha (x)}{- i sin (alpha pi)}</math>

若α为整数，则须对等号右边取极限值。另外，无论α是不是整数，下面的关系都成立：

:<math>H_{-alpha}^{(1)} (x)= e^{alpha pi i} H_{alpha}^{(1)} (x) </math>

:<math>H_{-alpha}^{(2)} (x)= e^{-alpha pi i} H_{alpha}^{(2)} (x) </math>

来自“http://www.zhdbk.com”中华大百科

===虚宗量的贝塞尔函数（修正贝塞尔函数）===

贝塞尔函数当宗量''x'' 为复数|复数时同样成立，并且当''x'' 为纯虚数时能得到一类重要情形——它们被称为第一类和第二类'''虚宗量的贝塞尔函数'''，或'''修正贝塞尔函数'''（有时还称为'''双曲型贝塞尔函数'''），定义为：

:<math>I_alpha(x) = i^{-alpha} J_alpha(ix) !</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>

以上形式保证了当宗量''x'' 为实数时，函数值亦为实数。这两个函数构成了下列'''修正贝塞尔方程'''（与一般贝塞尔方程的差别仅在两个正负号）的一个相互线性无关的解系：

:<math>x^2 frac{d^2 y}{dx^2} + x frac{dy}{dx} - (x^2 + alpha^2)y = 0.</math>

修正贝塞尔函数与一般贝塞尔函数的差别在于：一般贝塞尔函数随实宗量是振荡型的，而修正贝塞尔函数''I''α 和''K''α则分别是指数增长和指数衰减型的。和第一类贝塞尔函数''J''α一样，函数''I''α当α > 0 时在''x''=0 点等于0，当α=0时在''x''=0 点趋于有限值。类似地，''K''α在''x''=0 点发散（趋于无穷）。

{|align=center

''复数宗量的贝塞尔函数之零值''：<math>J_alpha (x) = 0</math>的解在α≥-1的情况下都是实数；阶数-2>α>-1的情况下，除了实数之外还有且仅有一对共轭的纯虚数解（G.N Watson 贝塞尔函数#参考文献|参考文献[5]）。

===球贝塞尔函数===

若使用分离变量法求解球坐标下的三维拉普拉斯方程，则可得到如下形式关于径向（''r'' 方向）分量的常微分方程：

:<math>x^2 frac{d^2 y}{dx^2} + 2x frac{dy}{dx} + [x^2 - n(n+1)]y = 0.</math>

关于上述方程的一对线性无关解称为'''球贝塞尔函数'''，分别用''j''''n''和''y''''n''表示（有时也记为''n''''n''）。这两个函数与一般贝塞尔函数''J''''n''和''Y''''n'' 存在关系：

:<math>j_n(x) = sqrt{frac{pi}{2x}} J_{n+1/2}(x),</math>

:<math>y_n(x) = sqrt{frac{pi}{2x}} Y_{n+1/2}(x) = (-1)^{n+1} sqrt{frac{pi}{2x}} J_{-n-1/2}(x).</math>

球贝塞尔函数也可写成：

:<math>j_n(x) = (-x)^n left(frac{1}{x}frac{d}{dx}

ight)^n,frac{sin x}{x} ,</math>

:<math>y_n(x) = -(-x)^n left(frac{1}{x}frac{d}{dx}

ight)^n,frac{cos x}{x}.</math>

0阶第一类球贝塞尔函数<math>j_0(x)</math>又称为sinc函数。头几阶整阶球贝塞尔函数的表达式分别为：

第一类：

:<math>j_0(x)=frac{sin x} {x}</math>

:<math>j_1(x)=frac{sin x} {x^2}- frac{cos x} {x}</math>

:<math>j_2(x)=left(frac{3} {x^2} - 1

ight)frac{sin x}{x} - frac{3cos x} {x^2}</math>

第二类：

:<math>y_0(x)=-j_{-1}(x)=-,frac{cos x} {x}</math>

:<math>y_1(x)=j_{-2}(x)=-,frac{cos x} {x^2}- frac{sin x} {x}</math>

:<math>y_2(x)=-j_{-3}(x)=left(-,frac{3}{x^2}+1

ight)frac{cos x}{x}- frac{3 sin x} {x^2}.</math>

还可以依照前面构造汉开尔函数相同的步骤构造所谓 球汉开尔函数：

:<math>h_n^{(1)}(x) = j_n(x) + i y_n(x)</math>

:<math>h_n^{(2)}(x) = j_n(x) - i y_n(x).</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>

而对实自变量''x''，''h''''n''(2)是上面''h''''n''(1)的复共轭（!! 表示'''双阶乘|阶乘'''）。由此我们可以通过得到''h''，再分离实部虚部，求出相应阶''j'' 和''h'' 的表达式，譬如''j''0(''x'') = sin(''x'')/''x''，''y''0(''x'') = -cos(''x'')/''x''，等等。

===黎卡提-贝塞尔函数===

黎卡提-贝塞尔函数（Riccati-Bessel functions）和球贝塞尔函数比较类似：

:<math>S_n(x)=x j_n(x)=sqrt{pi x/2}J_{n+1/2}(x)</math>

:<math>C_n(x)=-x y_n(x)=-sqrt{pi x/2}Y_{n+1/2}(x)</math>

:<math>zeta_n(x)=x h_n^{(2)}(x)=sqrt{pi x/2}H_{n+1/2}^{(2)}(x)=S_n(x)+iC_n(x)</math>

该函数满足方程：

:<math>x^2 frac{d^2 y}{dx^2} + [x^2 - n (n+1)] y = 0</math>

这个方程以及相应的黎卡提-贝塞尔解是德国物理学家古斯塔夫·米（w:Gustav Mie|Gustav Mie）于1908年研究电磁波在球状颗粒表面散射问题时提出的，后人将这种散射称为米氏散射（w:Mie theory|Mie scattering）。这个问题近几年的进展可参见文献 Du (2004)。

后人有时会遵从彼得·德拜|德拜（w:Peter Debye|Debye）在1909年的论文中的记法，用<math>psi_n,chi_n</math> 代替前面的<math>S_n,C_n</math>。

来自“http://www.zhdbk.com”中华大百科

==渐近形式==

贝塞尔函数在α非负时具有下面的渐近形式。当自变量''x'' 为小量，即<math>0 < x ll sqrt{alpha + 1}</math>时，有：

:<math>J_alpha(x)

ightarrow frac{1}{Gamma(alpha+1)} left( frac{x}{2}

ight) ^alpha </math>

:<math>Y_alpha(x)

ightarrow left{ egin{matrix}

frac{2}{pi} left[ ln (x/2) + gamma

ight] & mbox{if } alpha=0 \ \

-frac{Gamma(alpha)}{pi} left( frac{2}{x}

ight) ^alpha & mbox{if } alpha > 0

end{matrix}

ight.</math>

式中γ为欧拉-马歇罗尼常数（也叫欧拉常数，等于 0.5772156649...），Γ为Γ函数。对于很大的''x''，即<math>x gg |alpha^2 - 1/4|</math>时，渐近形式为：

:<math>J_alpha(x)

ightarrow sqrt{frac{2}{pi x}}

cos left( x-frac{alphapi}{2} - frac{pi}{4}

ight)</math>

:<math>Y_alpha(x)

ightarrow sqrt{frac{2}{pi x}}

sin left( x-frac{alphapi}{2} - frac{pi}{4}

ight).</math>

（α=1/2 时渐近号两边严格相等；参见前面对球贝塞尔函数的介绍）。其他形式贝塞尔函数的渐近形式可以从上面的式子直接推得。譬如，对大自变量<math>x gg |alpha^2 - 1/4|</math>，修正贝塞尔函数的渐近形式为：

:<math>I_alpha(x)

ightarrow frac{1}{sqrt{2pi x}} e^x,</math>

:<math>K_alpha(x)

ightarrow sqrt{frac{pi}{2x}} e^{-x}.</math>

对小自变量<math>0 < x ll sqrt{alpha + 1}</math>：

:<math>I_alpha(x)

ightarrow frac{1}{Gamma(alpha+1)} left( frac{x}{2}

ight) ^alpha </math>

:<math>K_alpha(x)

ightarrow left{ egin{matrix}

- ln (x/2) - gamma & mbox{if } alpha=0 \ \

frac{Gamma(alpha)}{2} left( frac{2}{x}

ight) ^alpha & mbox{if } alpha > 0

end{matrix}

ight.</math>

==性质==

整阶（α = ''n''）第一类贝塞尔函数''J''''n''常通过对其'''母函数'''（''generating function''）的罗朗级数（w:Laurent series|Laurent series）展开来定义：

:<math>e^{(x/2)(t-1/t)} = sum_{n=-infty}^infty J_n(x) t^n,</math>

上式得左边即为整阶第一类贝塞尔函数的母函数，这是丹麦天文学家w:Peter Andreas Hansen|汉森于1843年提出的。（这种定义也可以通过路径积分或其他方法推广到非整数阶）。整阶函数的另一个重要性质是下列'''雅可比-安格尔恒等式'''（''Jacobi-Anger identity''）：

:<math>e^{iz cos phi} = sum_{n=-infty}^infty i^n J_n(z) e^{inphi},</math>

利用这一等式可以将平面波展开成一系列柱面波的叠加，或者将频率调制|调频信号分解成傅里叶级数的叠加。

函数''J''α、''Y''α、''H''α(1)和''H''α(2)均满足递推关系：

:<math>Z_{alpha-1}(x) + Z_{alpha+1}(x) = frac{2alpha}{x} Z_alpha(x)</math>

:<math>Z_{alpha-1}(x) - Z_{alpha+1}(x) = 2frac{dZ_alpha}{dx}</math>

其中''Z''代表''J'', ''Y'', ''H''(1)或''H''(2)。（常将这两个恒等式联立推出其他关系）。从这组递推关系可以通过低阶贝塞尔函数（或它们的低阶导数）计算高阶贝塞尔函数（或它们的高阶导数）。特别地，有：

:<math>left( frac{d}{x dx}

ight)^m left[ x^alpha Z_{alpha} (x)

ight] = x^{alpha - m} Z_{alpha - m} (x)</math>

:<math>left( frac{d}{x dx}

ight)^m left[ frac{Z_alpha (x)}{x^alpha}

ight] = (-1)^m frac{Z_{alpha + m} (x)}{x^{alpha + m}}</math>

由于贝塞尔方程对应的作用算符除以''x'' 后便是一个（自伴随的）厄米算符（w:Hermitian|Hermitian），所以它的解在适当的边界条件下须满足正交性关系。特别地，可推得：

:<math>int_0^1 x J_alpha(x u_{alpha,m}) J_alpha(x u_{alpha,n}) dx = frac{delta_{m,n}}{2} J_{alpha+1}(u_{alpha,m})^2,</math>

其中α > -1，δ''m'',''n''为克罗内克尔δ，''u''α,m表示''J''α(''x'')的第''m'' 级零点。这个正交性关系可用于计算傅里叶-贝塞尔级数中各项的系数，以利用该级数将任意函数写成α固定、''m'' 变化的函数''J''α(''x'' ''u''α,m)的无穷叠加形式。（可以立即得到球贝塞尔函数相应的关系）。

另一个正交性关系是下列在α > -1/2时成立的“封闭方程”（''closure equation''）：

:<math>int_0^infty x J_alpha(ux) J_alpha(vx) dx = frac{1}{u} delta(u - v)</math>

其中δ为狄拉克δ函数。球贝塞尔函数的正交性条件为（当α > 0）：

:<math>int_0^infty x^2 j_alpha(ux) j_alpha(vx) dx = frac{pi}{2u^2} delta(u - v)</math>

贝塞尔方程的另一个重要性质与其朗斯基行列式（w:Wronskian|Wronskian）相关，由阿贝尔恒等式（w:Abel's identity|Abel's identity）得到：

:<math>A_alpha(x) frac{dB_alpha}{dx} - frac{dA_alpha}{dx} B_alpha(x) = frac{C_alpha}{x},</math>

其中''A''α 和''B''α是贝塞尔方程的任意两个解，''C''α是与''x'' 无关的常数（由α和贝塞尔函数的种类决定）。譬如，若''A''α = ''J''α、''B''α = ''Y''α，则''C''α is 2/π。该性质在修正贝塞尔函数中同样适用，譬如，若''A''α = ''I''α、''B''α = ''K''α，则''C''α为-1。

cs:Besselova funkce

de:Besselsche Differentialgleichung

en:Bessel function

es:Función de Bessel

fi:Besselin funktiot

fr:Fonction de Bessel

it:Funzioni di Bessel

ja:ベッセル関数

ko:베셀 함수

nl:Besselfunctie

pl:Funkcje Bessela

pt:Função de Bessel

ru:Функции Бесселя

sl:Besslova funkcija

sv:Besselfunktion

uk:Функція Неймана