贝塞尔函数 - 维基百科,自由的百科全书
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
贝塞尔函数
维基百科,自由的百科全书
贝塞尔函数(Bessel functions),是数学上的一类特
殊函数的总称。通常单说的贝塞尔函数指第一类贝塞
尔函数(Bessel function of the first kind)。一
般贝塞尔函数是下列常微分方程(一般称为贝塞尔方
程)的标准解函数:
这类方程的解是无法用初等函数系统地表示。
由于贝塞尔微分方程是二阶常微分方程,需要由两个独
立的函数来表示其标准解函数。典型的是使用第一类贝
塞尔函数和第二类贝塞尔函数来表示标准解函数:图1 贝塞尔函数的一个实例:一个紧绷的鼓面在中心受到敲击后的二阶振动振型,其振幅沿半径方向上的分布就是一个贝塞尔函数(考虑正负号)。实际生
活中受敲击的鼓面的振动是各阶类似振
动形态的叠加。
注意,由于 在 x=0 时候是发散的(无穷),当取 x=0 时,相关系数 必须为0时,才能获得有物理意义的结果。
贝塞尔函数的具体形式随上述方程中任意实数或复数α变化而变化(相应地,α被称为其对应贝塞尔函数的阶数)。实际应用中最常见的情形为α是整数n ,对应解称为n 阶贝塞尔函数。
尽管在上述微分方程中,α本身的正负号不改变方程的形式,但实际应用中仍习惯针对α和−α定义两种不同的贝塞尔函数(这样做能带来好处,比如消除了函数在α=0 点的不光滑性)。贝塞尔函数也被称为柱谐函数、圆柱函数或圆柱谐波,因为他们是于拉普拉斯方程在圆柱坐标上的求解过程中被发现的。
目录
1 历史
2 现实背景和应用范围
3 定义
3.1 第一类贝塞尔函数
3.1.1 贝塞尔积分
3.1.2 和超几何级数的关系
3.2 第二类贝塞尔函数(诺依曼函数)
3.3 第三类贝塞尔函数(汉克尔函数)
3.4 修正贝塞尔函数
3.5 球贝塞尔函数
3.6 黎卡提-贝塞尔函数
4 渐近形式
5 性质
6 参考文献
7 外部连接
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
历史
贝塞尔函数的几个正整数阶特例早在18世纪中叶就由瑞士数学家丹尼尔·伯努利在研究悬链振动时提出了,当时引起了数学界的兴趣。丹尼尔的叔叔雅各布·伯努利,欧拉、拉格朗日等数学大师对贝塞尔函数的研究作出过重要贡献。1817年,德国数学家贝塞尔在研究开普勒提出的三体引力系统的运动问题时,第一次系统地提出了贝塞尔函数的总体理论框架,后人以他的名字来命名了这种函数 [1] (http://www.britannica.com/eb/article-9078932) [2] (http://www-
history.mcs.st-andrews.ac.uk/Biographies/Bessel.html)。
现实背景和应用范围
贝塞尔方程是在圆柱坐标或球坐标下使用分离变量法求解拉普拉斯方程和亥姆霍兹方程时得到的(在圆柱域问题中得到的是整阶形式 α = n ;在球形域问题中得到的是半奇数阶形式 α =
n +½),因此贝塞尔函数在波的传播问题以及各种涉及有势场的问题中占有非常重要的地位,最典型的问题有:
在圆柱形波导中的电磁波传播问题;
圆柱体中的热传导问题;
圆形(或环形)薄膜的振动模态分析问题;
在其他一些领域,贝塞尔函数也相当有用。譬如在信号处理中的调频合成(FM synthesis)或凯泽窗(Kaiser window)的定义中,都要用到贝塞尔函数。
定义
贝塞尔方程是一个二阶常微分方程,必然存在两个线性无关的解。针对各种具体情况,人们提出了表示这些解的不同形式。下面分别介绍这些不同类型的贝塞尔函数。
第一类贝塞尔函数
第一类贝塞尔函数(Bessel
function of the first kind),
又称贝塞尔函数(Bessel
function),下文中有时会简称
为J函数,记作J α。
第一类α阶贝塞尔函数J α(x )是贝
塞尔方程当α为整数或α非负时的
解,须满足在x = 0 时有限。这样
选取和处理J α的原因见本主题下
面的性质介绍;另一种定义方法是
通过它在x = 0 点的泰勒级数展开
(或者更一般地通过幂级数展开,
这适用于α为非整数):
图2 0阶、1阶和2阶第一类贝塞尔函数(贝塞尔J函数)曲线
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
上式中为Γ函数(它可视为阶乘函数向非整型自变量的推广)。第一类贝塞尔函数的形状大致与按速率衰减的正弦或余弦函数类似(参见本页下面对它们渐进形式的介绍),但它们的零点并不是周期性的,另外随着x 的增加,零点的间隔会越来越接近周期性。图2所示为0阶、1阶和2阶第一类贝塞尔函数的曲线()。如果α不为整数,则和线性无关,可以构成微分方程的一个解系
。反之若是整数,那么上面两个函数之间满足如下关系:
于是两函数之间已不满足线性无关条件。为寻找在此情况下微分方程与解,需要定义第二类贝塞尔函数,定义过程将在后面的小节中给出。
贝塞尔积分为整数时贝塞尔函数的另一种定义方法由下面的积分给出:
线性无关的另一(为任意实数时的表达式见参考文献[2]第360页)
这个积分式就是贝塞尔当年提出的定义,而且他还从该定义中推出了函数的一些性质。另一种积
分表达式为:
和超几何级数的关系
贝塞尔函数可以用超几何级数表示成下面的形式:
ɑ为整数。由于函数线性相关的特性(用了一个就少了一个,所以要再构造一个),才需定义如下详细介绍的第二类贝塞尔函数。
第二类贝塞尔函数(诺依曼函数)
第二类贝塞尔函数(Bessel function of the second kind),又称诺伊曼函数(Neumannfunction),下文中有时会简称为Y函数,记作Y α。
第二类贝塞尔函数也许比第一类更为常用。 这种函数通常用Y α(x )表示,它们是贝塞尔方程的另一类解。x = 0 点是第二类贝塞尔函数的(无穷)奇点。
14-5-6贝塞尔函数 - 维基百科,自由的百科全书Y α(x )又被称为诺依曼函数(Neumann function),有时也记作N α(x )。它和J α(x )存在如下关系:
若α为整数(此时上式是型未定式)则取右端的极限值。
从前面对J α(x )的定义可以知道,
若α不为整数时,定义Y α是多余
的(因为贝塞尔方程的两个线性无
关解都已经用J函数表示出来
了)。另一方面,若α为整
数,Y α便可以和J α构成贝塞尔方
程的一个解系。与J函数类似,Y函
数正负整数阶之间也存在如下关
系:
J α(x )和Y α(x )均为沿负实半轴割
开的复平面内关于x 的全纯函数。
图3 0阶、1阶和2阶第二类贝塞尔函数(贝塞尔Y 函数)曲线图当α为整数时,复平面内不存在贝
塞尔函数的支点,所以J 和Y 均
为x 的整函数。若将x 固定,则贝
塞尔函数是α的整函数。图3所示为0阶、1阶和2阶第二类贝塞尔函数的曲线(
):
第三类贝塞尔函数(汉克尔函数)
第三类贝塞尔函数(Bessel function of the third kind),又称汉克尔函数(Hankelfunction)。
贝塞尔方程的另外一对重要的线性无关解称为汉克尔函数(Hankel functions)H α(1)(x ) 和H α(2)(x ),分别定义为:
其中i 为虚数单位。以上的线性组合也成为第三类贝塞尔函数;它们描述了二维波动方程的内行柱面波解和外行柱面波解("行"与在"行动"中同音)。
利用前面推出的关系可将汉克尔函数表示成:
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
若α为整数,则须对等号右边取极限值。另外,无论α是不是整数,下面的关系都成立:
修正贝塞尔函数
贝塞尔函数当变量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 第二类修正贝塞尔函数自变量的曲线()对实
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
复数变量的贝塞尔函数之零值:的解在α≥-1的情况下都是实数;阶数-2>α>-1的情况下,除了实数之外还有且仅有一对共轭的纯虚数解(G.N Watson 参考文献[5])。
第二类修正贝塞尔函数有时候被称为第三类修正贝塞尔函数(modified Bessel function ofthe third kind)。
球贝塞尔函数
若使用分离变量法求解球坐标下的三维拉普拉斯
方程,则可得到如下形式关于径向(r 方向)分
量的常微分方程:
图5-1 第一类球贝塞尔函数
)曲线(
图5-2 第二类球贝塞尔函数
)曲线(
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
关于上述方程的一对线性无关解称为球贝塞尔函数,分别用j
n 和y n 表示(有时也记为n n )。这两个函数与一般贝塞尔函数J n 和Y n 存在关系:
球贝塞尔函数也可写成:
0阶第一类球贝塞尔函数
第一类:又称为sinc函数。头几阶整阶球贝塞尔函数的表达式分别为:
第二类:
还可以依照前面构造汉克尔函数相同的步骤构造所谓 球汉克尔函数:
事实上,所有半奇数阶贝塞尔函数都可以写成由三角函数组成的封闭形式的表达式,球贝塞尔函数也同样可以。特别地,对所有非负整数n ,存在:
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
而对实自变量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)和球贝塞尔函数比较类似:
该函数满足方程:
这个方程以及相应的黎卡提-贝塞尔解是德国物理学家古斯塔夫·米(Gustav Mie)于1908年研究电磁波在球状颗粒表面散射问题时提出的,后人将这种散射称为米氏散射(Mie
scattering)。这个问题近几年的进展可参见文献 Du (2004)。
后人有时会遵从德拜(Debye)在1909年的论文中的记法,用 代替前面的。渐近形式
贝塞尔函数在α非负时具有下面的渐近形式。当自变量x 为小量,即
有:时,
式中γ为欧拉-马歇罗尼常数(也叫欧拉常数,等于 0.5772156649...),Γ为Γ函数。对于很大的x ,即时,渐近形式为:
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
(α=1/2 时渐近号两边严格相等;参见前面对球贝塞尔函数的介绍)。其他形式贝塞尔函数的
渐近形式可以从上面的式子直接推得。譬如,对大自变量
,修正贝塞尔函数的渐近形式为:
对小自变量:
性质
整阶(α = n )第一类贝塞尔函数J n 常通过对其母函数
(generating function)的罗朗级数(Laurent series)展开来定义:
上式得左边即为整阶第一类贝塞尔函数的母函数,这是丹麦天文学家汉森于1843年提出的。(这种定义也可以通过路径积分或其他方法推广到非整数阶)。整阶函数的另一个重要性质是下列雅
可比-安格尔恒等式(Jacobi-Anger identity):
利用这一等式可以将平面波展开成一系列柱面波的叠加,或者将调频信号分解成傅里叶级数的叠加。
函数J α、Y α、H α(1)和
H α
(2)均满足递推关系:
其中Z 代表J , Y , H (1)或H (2)。(常将这两个恒等式联立推出其他关系)。从这组递推关系可以通过低阶贝塞尔函数(或它们的低阶导数)计算高阶贝塞尔函数(或它们的高阶导数)。特别地,有:
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
由于贝塞尔方程对应的作用算符除以x 后便是一个(自伴随的)厄米算符(Hermitian),所以它的解在适当的边界条件下须满足正交性关系。特别地,可推得:
其中α > -1,δm , n 为克罗内克δ,u α,m表示J α(x )的第m 级零点。这个正交性关系可用于计算傅里叶-贝塞尔级数中各项的系数,以利用该级数将任意函数写成α固定、m 变化的函数J α(x u α,m)的无穷叠加形式。(可以立即得到球贝塞尔函数相应的关系)。
另一个正交性关系是下列在α > -1/2时成立的“封闭方程”(closure equation):
其中δ为狄拉克δ函数。球贝塞尔函数的正交性条件为(当α > 0):
贝塞尔方程的另一个重要性质与其朗斯基行列式(Wronskian)相关,由阿贝尔恒等式(Abel'sidentity)得到:
其中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/O·177
[2] Milton Abramowitz and Irene A. Stegun, eds., Handbook of Mathematical
Functions with Formulas, Graphs, and Mathematical Tables (Dover: New York,
1972) (英文)
Chapter 9 (http://www.math.sfu.ca/~cbm/aands/page_355.htm) 整阶贝塞尔函数
Section 9.1 (http://www.math.sfu.ca/~cbm/aands/page_358.htm) J, Y
(韦伯) and H (汉开尔)
Section 9.6 (http://www.math.sfu.ca/~cbm/aands/page_374.htm) 修正贝塞
尔函数(I和K)
Section 9.9 (http://www.math.sfu.ca/~cbm/aands/page_379.htm) 开尔文函
14-5-6贝塞尔函数 - 维基百科,自由的百科全书
数
Chapter 10 (http://www.math.sfu.ca/~cbm/aands/page_435.htm) 分数阶贝塞尔函数
Section 10.1 (http://www.math.sfu.ca/~cbm/aands/page_437.htm) 球贝塞
尔函数(j、y和h)
Section 10.2 (http://www.math.sfu.ca/~cbm/aands/page_443.htm) 修正球
贝塞尔函数(I和K)
Section 10.3 (http://www.math.sfu.ca/~cbm/aands/page_445.htm) 黎卡提-
贝塞尔函数
Section 10.4 (http://www.math.sfu.ca/~cbm/aands/page_446.htm) 艾里函
数(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) ISBN0-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
(http://www.efunda.com/math/bessel/bessel.cfm)
取自“http://zh.wikipedia.org/w/index.php?title=贝塞尔函数&oldid=30961797”
本页面最后修订于2014年4月10日 (星期四) 16:04。
本站的全部文字在知识共享 署名-相同方式共享 3.0协议之条款下提供,附加条款亦可能应用。(请参阅使用条款)
Wikipedia和维基百科标志是维基媒体基金会的注册商标;维基是维基媒体基金会的商标。
维基媒体基金会是在美国佛罗里达州登记的501(c)(3)免税、非营利、慈善机构。
zh.wikipedia.org/wiki/贝塞尔函数11/11