fbpx
维基百科

贝塞尔函数

贝塞尔函数Bessel functions),是数学上的一类特殊函数的总称。通常单说的贝塞尔函数第一类贝塞尔函数Bessel function of the first kind)。一般贝塞尔函数是下列常微分方程(一般称为贝塞尔方程)的标准解函数

图1 贝塞尔函数的一个实例:一个紧绷的鼓面在中心受到敲击后的二阶振动振型,其振幅沿半径方向上的分布就是一个贝塞尔函数(考虑正负号)。实际生活中受敲击的鼓面的振动是各阶类似振动形态的叠加。

该方程的通解无法用初等函数表示。

由於贝塞尔微分方程是二階常微分方程,需要由兩個獨立的函數來表示其标准解函数。典型的是使用第一类贝塞尔函数第二类贝塞尔函数來表示标准解函数:

注意,由於 在 x=0 時候是發散的(無窮),當取 x=0 時,相關係數 必須為0時,才能獲得有物理意義的結果。

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

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

贝塞尔函數也被稱為柱諧函數圓柱函數圓柱諧波,因為他們是於拉普拉斯方程圓柱坐標上的求解過程中被發現的。

历史

贝塞尔函数的几个正整数阶特例早在18世纪中叶就由瑞士数学家丹尼尔·伯努利在研究悬链振动时提出了,当时引起了数学界的兴趣。丹尼尔的叔叔雅各布·伯努利欧拉拉格朗日等数学大师对贝塞尔函数的研究作出过重要贡献。1817年,德国数学家贝塞尔在研究开普勒提出的三体引力系统的运动问题时,第一次系统地提出了贝塞尔函数的总体理论框架,后人以他的名字来命名了这种函数 [2](页面存档备份,存于互联网档案馆)。

现实背景和应用范围

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

在其他一些领域,贝塞尔函数也相当有用。譬如在ttf字体文件压缩,信号处理中的调频合成英语Frequency modulation synthesis凯泽窗的定义中,都要用到贝塞尔函数。

定义

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

第一类贝塞尔函数

 
图2 0阶、1阶和2阶第一类贝塞尔函数(贝塞尔J函数)曲线

第一类贝塞尔函数Bessel function of the first kind),又称贝塞尔函数Bessel function),下文中有时会简称为J函数,記作Jα

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

 

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

如果α不为整数,则  线性无关,可以构成微分方程的一个解系。反之若 是整数,那么上面两个函数之间满足如下关系:

 

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

贝塞尔积分

 为整数时贝塞尔函数的另一种定义方法由下面的积分给出:

 

 为任意实数时的表达式见参考文献[2]第360页)

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

 

和超几何级数的关系

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

 

ɑ為整數。由於函數線性相關的特性(用了一個就少了一個,所以要再構造一個),才需定義如下詳細介紹的第二類貝塞爾函數。

第二类贝塞尔函数(诺依曼函数)

 
图3 0阶、1阶和2阶第二类贝塞尔函数(贝塞尔Y 函数)曲线图

第二类贝塞尔函数Bessel function of the second kind),又称诺伊曼函数Neumann function),下文中有时会简称为Y函数,記作Yα

第二类贝塞尔函数也许比第一类更为常用。 这种函数通常用Yα(x)表示,它们是贝塞尔方程的另一类解。x = 0 点是第二类贝塞尔函数的(无穷)奇点。

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

 

若α为整数(此时上式是 未定式)则取右端的极限值。

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

 

Jα(x)和Yα(x)均为沿负实半轴割开的复平面内关于x全纯函数。当α为整数时,复平面内不存在贝塞尔函数的支点,所以JY 均为x整函数。若将x 固定,则贝塞尔函数是α的整函数。图3所示为0阶、1阶和2阶第二类贝塞尔函数 的曲线( ):

第三类贝塞尔函数(汉克尔函数)

第三类贝塞尔函数Bessel function of the third kind),又称汉克尔函数Hankel function)。

贝塞尔方程的另外一对重要的线性无关解称为汉克尔函数(Hankel functions)Hα(1)(x)和Hα(2)(x),分别定义为:

 
 

其中i虚数单位 。以上的线性组合也成为第三类贝塞尔函数;它们描述了二维波动方程外向行柱面波解和内向行柱面波解("行"与在"行动"中同音)。

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

 
 

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

 
 

修正贝塞尔函数

贝塞尔函数当变量x复数时同样成立,并且当x 为纯虚数时能得到一类重要情形——它们被称为第一类修正贝塞尔函数modified Bessel function of the first kind)和第二类修正贝塞尔函数modified Bessel function of the second kind),或虚变量的贝塞尔函数(有时还称为双曲型贝塞尔函数),定义为:

 
 

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

 

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

 
图4-1 第一类修正贝塞尔函数 对实自变量的曲线( 
 
图4-2 第二类修正贝塞尔函数 对实自变量的曲线( 


复数变量的贝塞尔函数之零值 的解在α≥-1的情况下都是实数;阶数-2>α>-1的情况下,除了实数之外还有且仅有一对共轭的纯虚数解(G.N Watson 参考文献[5])。

第二类修正贝塞尔函数有时候被称为第三类修正贝塞尔函数modified Bessel function of the third kind)。

球贝塞尔函数

 
图5-1 第一类球贝塞尔函数 曲线( 
 
图5-2 第二类球贝塞尔函数 曲线( 

若使用分离变量法求解球坐标下的三维亥姆霍兹方程,则可得到如下形式关于径向(r 方向)分量的常微分方程

 

关于上述方程的一对线性无关解称为球贝塞尔函数,分别用jnyn表示(有时也记为nn)。这两个函数与一般贝塞尔函数JnYn 存在关系:

 
 

球贝塞尔函数也可写成:

 
 

0阶第一类球贝塞尔函数 又称为sinc函数。头几阶整阶球贝塞尔函数的表达式分别为:

第一类:

 
 
 

第二类:

 
 
 

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

 
 

事实上,所有半奇数阶贝塞尔函数都可以写成由三角函数组成的封闭形式的表达式,球贝塞尔函数也同样可以。特别地,对所有非负整数n,存在:

 

而对实自变量xhn(2)是上面hn(1)的复共轭(!! 表示阶乘)。由此我们可以通过得到h,再分离实部虚部,求出相应阶jh 的表达式,譬如j0(x) = sin(x)/xy0(x) = -cos(x)/x,等等。

球贝塞尔函数的生成函数为:

 
 

黎卡提-贝塞尔函数:Sn, Cn, ξn, ζn

黎卡提-贝塞尔函数(Riccati-Bessel functions)和球贝塞尔函数比较类似:

 
 
 

该函数满足方程:

 

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

后人有时会遵从德拜(Debye)在1909年的论文中的记法,用  代替前面的 

渐近形式

贝塞尔函数在α非负时具有下面的渐近形式。当自变量x 为小量,即 时,有:

 
 

式中γ为歐拉-馬歇羅尼常數(也叫歐拉常數,等于 0.5772156649...),Γ为Γ函数。对于很大的x,即 时,渐近形式为:

 
 

(α=1/2 时渐近号两边严格相等;参见前面对球贝塞尔函数的介绍)。其他形式贝塞尔函数的渐近形式可以从上面的式子直接推得。譬如,对大自变量 ,修正贝塞尔函数的渐近形式为:

 
 

对小自变量 

 
 

性质

整阶(α = n)第一类贝塞尔函数Jn常通过对其母函数generating function)的罗朗级数(Laurent series)展开来定义:

 

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

 

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

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

 
 

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

 
 

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

 

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

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

 

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

 

贝塞尔方程的另一个重要性质与其朗斯基行列式(Wronskian)相关,由阿贝尔恒等式(Abel's identity)得到:

 

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

參見

参考文献

  • [1] 严镇军编,《数学物理方程》,第二版,中国科学技术大学出版社,合肥,2002,第82页~第123页,ISBN 7-312-00799-6
  • [2] Milton Abramowitz and Irene A. Stegun, eds., Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (Dover: New York, 1972) (英文)
    • Chapter 9(页面存档备份,存于互联网档案馆) 整阶贝塞尔函数
      • Section 9.1(页面存档备份,存于互联网档案馆) J, Y (韦伯) and H (汉开尔)
      • Section 9.6(页面存档备份,存于互联网档案馆) 修正贝塞尔函数(I和K)
      • Section 9.9(页面存档备份,存于互联网档案馆) 开尔文函数
    • Chapter 10(页面存档备份,存于互联网档案馆) 分数阶贝塞尔函数
      • Section 10.1(页面存档备份,存于互联网档案馆) 球贝塞尔函数(j、y和h)
      • Section 10.2(页面存档备份,存于互联网档案馆) 修正球贝塞尔函数(I和K)
      • Section 10.3(页面存档备份,存于互联网档案馆) 黎卡提-贝塞尔函数
      • Section 10.4(页面存档备份,存于互联网档案馆) 艾里函数(Airy functions)
  • [3] George B. Arfken and Hans J. Weber, Mathematical Methods for Physicists (Harcourt: San Diego, 2001).
  • [4] Frank Bowman, Introduction to Bessel Functions (Dover: New York, 1958) ISBN 0-486-60462-4.
  • [5] G. N. Watson, A Treatise on the Theory of Bessel Functions, Second Edition, (1966) Cambridge University Press.
  • [6] G. Mie, "Beiträge zur Optik trüber Medien, speziell kolloidaler Metallösungen", Ann. Phys. Leipzig 25(1908), p.377.
  • [7] Hong Du, "Mie-scattering calculation," Applied Optics 43 (9), 1951-1956 (2004).

外部链接

  • Engineering Fundamentals - Bessel Function(页面存档备份,存于互联网档案馆

贝塞尔函数, bessel, functions, 是数学上的一类特殊函数的总称, 通常单说的指第一类, bessel, function, first, kind, 一般是下列常微分方程, 一般称为贝塞尔方程, 的标准解函数y, displaystyle, 图1, 的一个实例, 一个紧绷的鼓面在中心受到敲击后的二阶振动振型, 其振幅沿半径方向上的分布就是一个, 考虑正负号, 实际生活中受敲击的鼓面的振动是各阶类似振动形态的叠加, displaystyle, frac, frac, alpha, 该方程的通解无法用. 贝塞尔函数 Bessel functions 是数学上的一类特殊函数的总称 通常单说的贝塞尔函数指第一类贝塞尔函数 Bessel function of the first kind 一般贝塞尔函数是下列常微分方程 一般称为贝塞尔方程 的标准解函数y x displaystyle y x 图1 贝塞尔函数的一个实例 一个紧绷的鼓面在中心受到敲击后的二阶振动振型 其振幅沿半径方向上的分布就是一个贝塞尔函数 考虑正负号 实际生活中受敲击的鼓面的振动是各阶类似振动形态的叠加 x 2 d 2 y d x 2 x d y d x x 2 a 2 y 0 displaystyle x 2 frac d 2 y dx 2 x frac dy dx x 2 alpha 2 y 0 该方程的通解无法用初等函数表示 由於贝塞尔微分方程是二階常微分方程 需要由兩個獨立的函數來表示其标准解函数 典型的是使用第一类贝塞尔函数和第二类贝塞尔函数來表示标准解函数 y x c 1 J a x c 2 Y a x displaystyle y x c 1 J alpha x c 2 Y alpha x 注意 由於 Y a x displaystyle Y alpha x 在 x 0 時候是發散的 無窮 當取 x 0 時 相關係數 c 2 displaystyle c 2 必須為0時 才能獲得有物理意義的結果 贝塞尔函数的具体形式随上述方程中任意实数或複數a变化而变化 相应地 a被称为其对应贝塞尔函数的阶数 实际应用中最常见的情形为a是整数n 对应解称为n 阶贝塞尔函数 尽管在上述微分方程中 a本身的正负号不改变方程的形式 但实际应用中仍习惯针对a和 a定义两种不同的贝塞尔函数 这样做能带来好处 比如消除了函数在a 0 点的不光滑性 贝塞尔函數也被稱為柱諧函數 圓柱函數或圓柱諧波 因為他們是於拉普拉斯方程在圓柱坐標上的求解過程中被發現的 目录 1 历史 2 现实背景和应用范围 3 定义 3 1 第一类贝塞尔函数 3 1 1 贝塞尔积分 3 1 2 和超几何级数的关系 3 2 第二类贝塞尔函数 诺依曼函数 3 3 第三类贝塞尔函数 汉克尔函数 3 4 修正贝塞尔函数 3 5 球贝塞尔函数 3 6 黎卡提 贝塞尔函数 Sn Cn 3n zn 4 渐近形式 5 性质 6 參見 7 参考文献 8 外部链接历史 编辑贝塞尔函数的几个正整数阶特例早在18世纪中叶就由瑞士数学家丹尼尔 伯努利在研究悬链振动时提出了 当时引起了数学界的兴趣 丹尼尔的叔叔雅各布 伯努利 欧拉 拉格朗日等数学大师对贝塞尔函数的研究作出过重要贡献 1817年 德国数学家贝塞尔在研究开普勒提出的三体引力系统的运动问题时 第一次系统地提出了贝塞尔函数的总体理论框架 后人以他的名字来命名了这种函数 1 2 页面存档备份 存于互联网档案馆 现实背景和应用范围 编辑贝塞尔方程是在圆柱坐标或球坐标下使用分离变量法求解拉普拉斯方程和亥姆霍兹方程时得到的 在圆柱域问题中得到的是整阶形式 a n 在球形域问题中得到的是半奇数阶形式 a n 因此贝塞尔函数在波的传播问题以及各种涉及有势场的问题中占有非常重要的地位 最典型的问题有 在圆柱形波导中的电磁波传播问题 圆柱体中的热传导问题 圆形 或环形 薄膜的振动模态分析问题 在其他一些领域 贝塞尔函数也相当有用 譬如在ttf字体文件压缩 信号处理中的调频合成 英语 Frequency modulation synthesis 或凯泽窗的定义中 都要用到贝塞尔函数 定义 编辑贝塞尔方程是一个二阶常微分方程 必然存在两个线性无关的解 针对各种具体情况 人们提出了表示这些解的不同形式 下面分别介绍这些不同类型的贝塞尔函数 第一类贝塞尔函数 编辑 图2 0阶 1阶和2阶第一类贝塞尔函数 贝塞尔J函数 曲线 第一类贝塞尔函数 Bessel function of the first kind 又称贝塞尔函数 Bessel function 下文中有时会简称为J函数 記作Ja 第一类a阶贝塞尔函数Ja x 是贝塞尔方程当a为整数或a非负时的解 须满足在x 0 时有限 这样选取和处理Ja的原因见本主题下面的性质介绍 另一种定义方法是通过它在x 0 点的泰勒级数展开 或者更一般地通过幂级数展开 这适用于a为非整数 J a x m 0 1 m m G m a 1 x 2 2 m a displaystyle J alpha x sum m 0 infty frac 1 m m Gamma m alpha 1 left frac x 2 right 2m alpha 上式中G z displaystyle Gamma z 为G函数 它可视为阶乘函数向非整型自变量的推广 第一类贝塞尔函数的形状大致与按1 x displaystyle 1 sqrt x 速率衰减的正弦或余弦函数类似 参见本页下面对它们渐近形式的介绍 但它们的零点并不是周期性的 另外随着x的增加 零点的间隔会越来越接近周期性 图2所示为0阶 1阶和2阶第一类贝塞尔函数J a x displaystyle J alpha x 的曲线 a 0 1 2 displaystyle alpha 0 1 2 如果a不为整数 则J a x displaystyle J alpha x 和J a x displaystyle J alpha x 线性无关 可以构成微分方程的一个解系 反之若a displaystyle alpha 是整数 那么上面两个函数之间满足如下关系 J a x 1 a J a x displaystyle J alpha x 1 alpha J alpha x 于是两函数之间已不满足线性无关条件 为寻找在此情况下微分方程与J a x displaystyle J alpha x 线性无关的另一解 需要定义第二类贝塞尔函数 定义过程将在后面的小节中给出 贝塞尔积分 编辑 a displaystyle alpha 为整数时贝塞尔函数的另一种定义方法由下面的积分给出 J a x 1 2 p 0 2 p cos a t x sin t d t displaystyle J alpha x frac 1 2 pi int 0 2 pi cos alpha tau x sin tau d tau a displaystyle alpha 为任意实数时的表达式见参考文献 2 第360页 这个积分式就是贝塞尔当年提出的定义 而且他还从该定义中推出了函数的一些性质 另一种积分表达式为 J a x 1 2 p p p e i a t x sin t d t displaystyle J alpha x frac 1 2 pi int pi pi e i alpha tau x sin tau d tau 和超几何级数的关系 编辑 贝塞尔函数可以用超几何级数表示成下面的形式 J a z z 2 a G a 1 0 F 1 a 1 z 2 4 displaystyle J alpha z frac z 2 alpha Gamma alpha 1 0 F 1 alpha 1 z 2 4 ɑ為整數 由於函數線性相關的特性 用了一個就少了一個 所以要再構造一個 才需定義如下詳細介紹的第二類貝塞爾函數 第二类贝塞尔函数 诺依曼函数 编辑 图3 0阶 1阶和2阶第二类贝塞尔函数 贝塞尔Y 函数 曲线图 第二类贝塞尔函数 Bessel function of the second kind 又称诺伊曼函数 Neumann function 下文中有时会简称为Y函数 記作Ya 第二类贝塞尔函数也许比第一类更为常用 这种函数通常用Ya x 表示 它们是贝塞尔方程的另一类解 x 0 点是第二类贝塞尔函数的 无穷 奇点 Ya x 又被称为诺依曼函数 Neumann function 有时也记作Na x 它和Ja x 存在如下关系 Y a x J a x cos a p J a x sin a p displaystyle Y alpha x frac J alpha x cos alpha pi J alpha x sin alpha pi 若a为整数 此时上式是0 0 displaystyle frac 0 0 型未定式 则取右端的极限值 从前面对Ja x 的定义可以知道 若a不为整数时 定义Ya是多余的 因为贝塞尔方程的两个线性无关解都已经用J函数表示出来了 另一方面 若a为整数 Ya便可以和Ja构成贝塞尔方程的一个解系 与J函数类似 Y函数正负整数阶之间也存在如下关系 Y n x 1 n Y n x displaystyle Y n x 1 n Y n x Ja x 和Ya x 均为沿负实半轴割开的复平面内关于x的全纯函数 当a为整数时 复平面内不存在贝塞尔函数的支点 所以J 和Y 均为x 的整函数 若将x 固定 则贝塞尔函数是a的整函数 图3所示为0阶 1阶和2阶第二类贝塞尔函数Y a x displaystyle Y alpha x 的曲线 a 0 1 2 displaystyle alpha 0 1 2 第三类贝塞尔函数 汉克尔函数 编辑 第三类贝塞尔函数 Bessel function of the third kind 又称汉克尔函数 Hankel function 贝塞尔方程的另外一对重要的线性无关解称为汉克尔函数 Hankel functions Ha 1 x 和Ha 2 x 分别定义为 H a 1 x J a x i Y a x displaystyle H alpha 1 x J alpha x iY alpha x H a 2 x J a x i Y a x displaystyle H alpha 2 x J alpha x iY alpha x 其中i 为虚数单位 1 displaystyle sqrt 1 以上的线性组合也成为第三类贝塞尔函数 它们描述了二维波动方程的外向行柱面波解和内向行柱面波解 行 与在 行动 中同音 利用前面推出的关系可将汉克尔函数表示成 H a 1 x J a x e a p i J a x i sin a p displaystyle H alpha 1 x frac J alpha x e alpha pi i J alpha x i sin alpha pi H a 2 x J a x e a p i J a x i sin a p displaystyle H alpha 2 x frac J alpha x e alpha pi i J alpha x i sin alpha pi 若a为整数 则须对等号右边取极限值 另外 无论a是不是整数 下面的关系都成立 H a 1 x e a p i H a 1 x displaystyle H alpha 1 x e alpha pi i H alpha 1 x H a 2 x e a p i H a 2 x displaystyle H alpha 2 x e alpha pi i H alpha 2 x 修正贝塞尔函数 编辑 贝塞尔函数当变量x 为复数时同样成立 并且当x 为纯虚数时能得到一类重要情形 它们被称为第一类修正贝塞尔函数 modified Bessel function of the first kind 和第二类修正贝塞尔函数 modified Bessel function of the second kind 或虚变量的贝塞尔函数 有时还称为双曲型贝塞尔函数 定义为 I a x i a J a i x displaystyle I alpha x i alpha J alpha ix K a x p 2 I a x I a x sin a p p 2 i a 1 H a 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 为实数时 函数值亦为实数 这两个函数构成了下列修正贝塞尔方程 与一般贝塞尔方程的差别仅在两个正负号 的一个相互线性无关的解系 x 2 d 2 y d x 2 x d y d x x 2 a 2 y 0 displaystyle x 2 frac d 2 y dx 2 x frac dy dx x 2 alpha 2 y 0 修正贝塞尔函数与一般贝塞尔函数的差别在于 一般贝塞尔函数随实变量是振荡型的 而修正贝塞尔函数Ia 和Ka则分别是指数增长和指数衰减型的 和第一类贝塞尔函数Ja一样 函数Ia当a gt 0 时在x 0 点等于0 当a 0时在x 0 点趋于有限值 类似地 Ka在x 0 点发散 趋于无穷 图4 1 第一类修正贝塞尔函数I a x displaystyle I alpha x 对实自变量的曲线 a 0 1 2 displaystyle alpha 0 1 2 图4 2 第二类修正贝塞尔函数K a x displaystyle K alpha x 对实自变量的曲线 a 0 1 2 displaystyle alpha 0 1 2 复数变量的贝塞尔函数之零值 J a x 0 displaystyle J alpha x 0 的解在a 1的情况下都是实数 阶数 2 gt a gt 1的情况下 除了实数之外还有且仅有一对共轭的纯虚数解 G N Watson 参考文献 5 第二类修正贝塞尔函数有时候被称为第三类修正贝塞尔函数 modified Bessel function of the third kind 球贝塞尔函数 编辑 图5 1 第一类球贝塞尔函数j n x displaystyle j n x 曲线 n 0 1 2 displaystyle n 0 1 2 图5 2 第二类球贝塞尔函数y n x displaystyle y n x 曲线 n 0 1 2 displaystyle n 0 1 2 若使用分离变量法求解球坐标下的三维亥姆霍兹方程 则可得到如下形式关于径向 r 方向 分量的常微分方程 x 2 d 2 y d x 2 2 x d y d x x 2 n n 1 y 0 displaystyle x 2 frac d 2 y dx 2 2x frac dy dx x 2 n n 1 y 0 关于上述方程的一对线性无关解称为球贝塞尔函数 分别用jn和yn表示 有时也记为nn 这两个函数与一般贝塞尔函数Jn和Yn 存在关系 j n x p 2 x J n 1 2 x displaystyle j n x sqrt frac pi 2x J n 1 2 x y n x p 2 x Y n 1 2 x 1 n 1 p 2 x J n 1 2 x displaystyle y n x sqrt frac pi 2x Y n 1 2 x 1 n 1 sqrt frac pi 2x J n 1 2 x 球贝塞尔函数也可写成 j n x x n 1 x d d x n sin x x displaystyle j n x x n left frac 1 x frac d dx right n frac sin x x y n x x n 1 x d d x n cos x x displaystyle y n x x n left frac 1 x frac d dx right n frac cos x x 0阶第一类球贝塞尔函数j 0 x displaystyle j 0 x 又称为sinc函数 头几阶整阶球贝塞尔函数的表达式分别为 第一类 j 0 x sin x x displaystyle j 0 x frac sin x x j 1 x sin x x 2 cos x x displaystyle j 1 x frac sin x x 2 frac cos x x j 2 x 3 x 2 1 sin x x 3 cos x x 2 displaystyle j 2 x left frac 3 x 2 1 right frac sin x x frac 3 cos x x 2 第二类 y 0 x j 1 x cos x x displaystyle y 0 x j 1 x frac cos x x y 1 x j 2 x cos x x 2 sin x x displaystyle y 1 x j 2 x frac cos x x 2 frac sin x x y 2 x j 3 x 3 x 2 1 cos x x 3 sin x x 2 displaystyle y 2 x j 3 x left frac 3 x 2 1 right frac cos x x frac 3 sin x x 2 还可以依照前面构造汉克尔函数相同的步骤构造所谓球汉克尔函数 h n 1 x j n x i y n x displaystyle h n 1 x j n x iy n x h n 2 x j n x i y n x displaystyle h n 2 x j n x iy n x 事实上 所有半奇数阶贝塞尔函数都可以写成由三角函数组成的封闭形式的表达式 球贝塞尔函数也同样可以 特别地 对所有非负整数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 而对实自变量x hn 2 是上面hn 1 的复共轭 表示双阶乘 由此我们可以通过得到h 再分离实部虚部 求出相应阶j 和h 的表达式 譬如j0 x sin x x y0 x cos x x 等等 球贝塞尔函数的生成函数为 1 z cos z 2 2 z t n 0 t n n j n 1 z displaystyle frac 1 z cos sqrt z 2 2zt sum n 0 infty frac t n n j n 1 z 1 z sin z 2 2 z t n 0 t n n y n 1 z displaystyle frac 1 z sin sqrt z 2 2zt sum n 0 infty frac t n n y n 1 z 黎卡提 贝塞尔函数 Sn Cn 3n zn 编辑 黎卡提 贝塞尔函数 Riccati Bessel functions 和球贝塞尔函数比较类似 S n x x j n x p x 2 J n 1 2 x displaystyle S n x xj n x sqrt pi x 2 J n 1 2 x C n x x y n x p x 2 Y n 1 2 x displaystyle C n x xy n x sqrt pi x 2 Y n 1 2 x z n x x h n 2 x p x 2 H n 1 2 2 x S n x i C n x displaystyle zeta n x xh n 2 x sqrt pi x 2 H n 1 2 2 x S n x iC n x 该函数满足方程 x 2 d 2 y d x 2 x 2 n n 1 y 0 displaystyle x 2 frac d 2 y dx 2 x 2 n n 1 y 0 这个方程以及相应的黎卡提 贝塞尔解是德国物理学家古斯塔夫 米 Gustav Mie 于1908年研究电磁波在球状颗粒表面散射问题时提出的 后人将这种散射称为米氏散射 Mie scattering 这个问题近几年的进展可参见文献 Du 2004 后人有时会遵从德拜 Debye 在1909年的论文中的记法 用ps n x n displaystyle psi n chi n 代替前面的S n C n displaystyle S n C n 渐近形式 编辑贝塞尔函数在a非负时具有下面的渐近形式 当自变量x 为小量 即0 lt x a 1 displaystyle 0 lt x ll sqrt alpha 1 时 有 J a x 1 G a 1 x 2 a displaystyle J alpha x rightarrow frac 1 Gamma alpha 1 left frac x 2 right alpha Y a x 2 p ln x 2 g if a 0 G a p 2 x a if a gt 0 displaystyle Y alpha x rightarrow left begin matrix frac 2 pi left ln x 2 gamma right amp mbox if alpha 0 frac Gamma alpha pi left frac 2 x right alpha amp mbox if alpha gt 0 end matrix right 式中g为歐拉 馬歇羅尼常數 也叫歐拉常數 等于 0 5772156649 G为G函数 对于很大的x 即x a 2 1 4 displaystyle x gg alpha 2 1 4 时 渐近形式为 J a x 2 p x cos x a p 2 p 4 displaystyle J alpha x rightarrow sqrt frac 2 pi x cos left x frac alpha pi 2 frac pi 4 right Y a x 2 p x sin x a p 2 p 4 displaystyle Y alpha x rightarrow sqrt frac 2 pi x sin left x frac alpha pi 2 frac pi 4 right a 1 2 时渐近号两边严格相等 参见前面对球贝塞尔函数的介绍 其他形式贝塞尔函数的渐近形式可以从上面的式子直接推得 譬如 对大自变量x a 2 1 4 displaystyle x gg alpha 2 1 4 修正贝塞尔函数的渐近形式为 I a x 1 2 p x e x displaystyle I alpha x rightarrow frac 1 sqrt 2 pi x e x K a x p 2 x e x displaystyle K alpha x rightarrow sqrt frac pi 2x e x 对小自变量0 lt x a 1 displaystyle 0 lt x ll sqrt alpha 1 I a x 1 G a 1 x 2 a displaystyle I alpha x rightarrow frac 1 Gamma alpha 1 left frac x 2 right alpha K a x ln x 2 g if a 0 G a 2 2 x a if a gt 0 displaystyle K alpha x rightarrow left begin matrix ln x 2 gamma amp mbox if alpha 0 frac Gamma alpha 2 left frac 2 x right alpha amp mbox if alpha gt 0 end matrix right 性质 编辑整阶 a n 第一类贝塞尔函数Jn常通过对其母函数 generating function 的罗朗级数 Laurent series 展开来定义 e x 2 t 1 t n J n x t n displaystyle e x 2 t 1 t sum n infty infty J n x t n 上式得左边即为整阶第一类贝塞尔函数的母函数 这是丹麦天文学家汉森于1843年提出的 这种定义也可以通过路径积分或其他方法推广到非整数阶 整阶函数的另一个重要性质是下列雅可比 安格尔恒等式 Jacobi Anger identity e i z cos ϕ n i n J n z e i n ϕ displaystyle e iz cos phi sum n infty infty i n J n z e in phi 利用这一等式可以将平面波展开成一系列柱面波的叠加 或者将调频信号分解成傅里叶级数的叠加 函数Ja Ya Ha 1 和Ha 2 均满足递推关系 Z a 1 x Z a 1 x 2 a x Z a x displaystyle Z alpha 1 x Z alpha 1 x frac 2 alpha x Z alpha x Z a 1 x Z a 1 x 2 d Z a d x displaystyle Z alpha 1 x Z alpha 1 x 2 frac dZ alpha dx 其中Z代表J Y H 1 或H 2 常将这两个恒等式联立推出其他关系 从这组递推关系可以通过低阶贝塞尔函数 或它们的低阶导数 计算高阶贝塞尔函数 或它们的高阶导数 特别地 有 d x d x m x a Z a x x a m Z a m x displaystyle left frac d xdx right m left x alpha Z alpha x right x alpha m Z alpha m x d x d x m Z a x x a 1 m Z a m x x a m displaystyle left frac d xdx right m left frac Z alpha x x alpha right 1 m frac Z alpha m x x alpha m 由于贝塞尔方程对应的作用算符除以x 后便是一个 自伴随的 厄米算符 Hermitian 所以它的解在适当的边界条件下须满足正交性关系 特别地 可推得 0 1 x J a x u a m J a x u a n d x d m n 2 J a 1 u a m 2 displaystyle int 0 1 xJ alpha xu alpha m J alpha xu alpha n dx frac delta m n 2 J alpha 1 u alpha m 2 其中a gt 1 dm n为克罗内克d ua m表示Ja x 的第m 级零点 这个正交性关系可用于计算傅里叶 贝塞尔级数中各项的系数 以利用该级数将任意函数写成a固定 m 变化的函数Ja x ua m 的无穷叠加形式 可以立即得到球贝塞尔函数相应的关系 另一个正交性关系是下列在a gt 1 2时成立的 封闭方程 closure equation 0 x J a u x J a v x d x 1 u d u v displaystyle int 0 infty xJ alpha ux J alpha vx dx frac 1 u delta u v 其中d为狄拉克d函数 球贝塞尔函数的正交性条件为 当a gt 0 0 x 2 j a u x j a v x d x p 2 u 2 d u v displaystyle int 0 infty x 2 j alpha ux j alpha vx dx frac pi 2u 2 delta u v 贝塞尔方程的另一个重要性质与其朗斯基行列式 Wronskian 相关 由阿贝尔恒等式 Abel s identity 得到 A a x d B a d x d A a d x B a x C a x displaystyle A alpha x frac dB alpha dx frac dA alpha dx B alpha x frac C alpha x 其中Aa 和Ba是贝塞尔方程的任意两个解 Ca是与x 无关的常数 由a和贝塞尔函数的种类决定 譬如 若Aa Ja Ba Ya 则Ca is 2 p 该性质在修正贝塞尔函数中同样适用 譬如 若Aa Ia Ba Ka 则Ca为 1 參見 编辑Hankel變換 以貝塞爾函數作展開 参考文献 编辑 1 严镇军编 数学物理方程 第二版 中国科学技术大学出版社 合肥 2002 第82页 第123页 ISBN 7 312 00799 6 2 Milton Abramowitz and Irene A Stegun eds Handbook of Mathematical Functions with Formulas Graphs and Mathematical Tables Dover New York 1972 英文 Chapter 9 页面存档备份 存于互联网档案馆 整阶贝塞尔函数 Section 9 1 页面存档备份 存于互联网档案馆 J Y 韦伯 and H 汉开尔 Section 9 6 页面存档备份 存于互联网档案馆 修正贝塞尔函数 I和K Section 9 9 页面存档备份 存于互联网档案馆 开尔文函数 Chapter 10 页面存档备份 存于互联网档案馆 分数阶贝塞尔函数 Section 10 1 页面存档备份 存于互联网档案馆 球贝塞尔函数 j y和h Section 10 2 页面存档备份 存于互联网档案馆 修正球贝塞尔函数 I和K Section 10 3 页面存档备份 存于互联网档案馆 黎卡提 贝塞尔函数 Section 10 4 页面存档备份 存于互联网档案馆 艾里函数 Airy functions 3 George B Arfken and Hans J Weber Mathematical Methods for Physicists Harcourt San Diego 2001 4 Frank Bowman Introduction to Bessel Functions Dover New York 1958 ISBN 0 486 60462 4 5 G N Watson A Treatise on the Theory of Bessel Functions Second Edition 1966 Cambridge University Press 6 G Mie Beitrage zur Optik truber Medien speziell kolloidaler Metallosungen Ann Phys Leipzig 25 1908 p 377 7 Hong Du Mie scattering calculation Applied Optics 43 9 1951 1956 2004 外部链接 编辑Engineering Fundamentals Bessel Function 页面存档备份 存于互联网档案馆 取自 https zh wikipedia org w index php title 贝塞尔函数 amp oldid 71894371, 维基百科,wiki,书籍,书籍,图书馆,

文章

,阅读,下载,免费,免费下载,mp3,视频,mp4,3gp, jpg,jpeg,gif,png,图片,音乐,歌曲,电影,书籍,游戏,游戏。