Bessel函数介绍.doc

贝塞尔函数（Bessel functions）是数学上的一类特殊函数的总称。一般贝塞尔函数是下列常微分方程（一般称为贝塞尔方程）的标准解函数y(x)：这类方程的解是无法用初等函数系统地表示的。贝塞尔函数的具体形式随上述方程中任意实数变化而变化（相应地，被称为其对应贝塞尔函数的阶数）。实际应用中最常见的情形为是整数n，对应解称为n 阶贝塞尔函数。尽管在上述微分方程中，本身的正负号不改变方程的形式，但实际应用中仍习惯针对和定义两种不同的贝塞尔函数（这样做能带来好处，比如消除了函数在=0 点的不光滑性）。历史贝塞尔函数的几个正整数阶特例早在18世纪中叶就由瑞士数学家丹尼尔伯努利在研究悬链振动时提出了，当时引起了数学界的兴趣。丹尼尔的叔叔雅各布伯努利，欧拉、拉格朗日等数学大师对贝塞尔函数的研究作出过重要贡献。1817年，德国数学家贝塞尔在研究开普勒提出的三体引力系统的运动问题时，第一次系统地提出了贝塞尔函数的总体理论框架，后人以他的名字来命名了这种函数 1 2。现实背景和应用范围贝塞尔方程是在柱坐标或球坐标下使用分离变量法求解拉普拉斯方程和亥姆霍兹方程时得到的（在圆柱域问题中得到的是整阶形式 = n；在球形域问题中得到的是半奇数阶形式 = n+），因此贝塞尔函数在波动问题以及各种涉及有势场的问题中占有非常重要的地位，最典型的问题有： 在圆柱形波导中的电磁波传播问题； 圆柱体中的热传导问题； 圆形（或环形）薄膜的振动模态分析问题；在其他一些领域，贝塞尔函数也相当有用。譬如在信号处理中的调频合成（FM synthesis）或凯泽窗（Kaiser window）的定义中，都要用到贝塞尔函数。定义 贝塞尔方程是一个二阶常微分方程，必然存在两个线性无关的解。针对各种具体情况，人们提出了表示这些解的不同形式。下面分别介绍这些不同类型的贝塞尔函数。第一类贝塞尔函数图2 0阶、1阶和2阶第一类贝塞尔函数（贝塞尔J函数）曲线（在下文中，第一类贝塞尔函数有时会简称为“J函数”，敬请读者留意。）第一类阶贝塞尔函数J(x)是贝塞尔方程当为整数或非负时的解，须满足在x = 0 时有限。这样选取和处理J的原因见本主题下面的性质介绍；另一种定义方法是通过它在x = 0 点的泰勒级数展开（或者更一般地通过幂级数展开，这适用于为非整数）：上式中(z)为函数（它可视为阶乘函数向非整型自变量的推广）。第一类贝塞尔函数的形状大致与按速率衰减的正弦或余弦函数类似（参见本页下面对它们渐进形式的介绍），但它们的零点并不是周期性的，另外随着x的增加，零点的间隔会越来越接近周期性。图2所示为0阶、1阶和2阶第一类贝塞尔函数J(x)的曲线（ = 0,1,2）。如果不为整数，则J(x)和J (x)线性无关，可以构成微分方程的一个解系。反之若是整数，那么上面两个函数之间满足如下关系：于是两函数之间已不满足线性无关条件。为寻找在此情况下微分方程与J(x)线性无关的另一解，需要定义第二类贝塞尔函数，定义过程将在后面的小节中给出。贝塞尔积分为整数时贝塞尔函数的另一种定义方法由下面的积分给出：（为任意实数时的表达式见参考文献2第360页）这个积分式就是贝塞尔当年提出的定义，而且他还从该定义中推出了函数的一些性质。另一种积分表达式为：和超几何级数的关系贝塞尔函数可以用超几何级数表示成下面的形式：第二类贝塞尔函数（诺依曼函数）图3 0阶、1阶和2阶第二类贝塞尔函数（贝塞尔Y 函数）曲线图（在下文中，第二类贝塞尔函数有时会简称为“Y函数”，敬请读者留意。）第二类贝塞尔函数也许比第一类更为常用。 这种函数通常用Y(x)表示，它们是贝塞尔方程的另一类解。x = 0 点是第二类贝塞尔函数的（无穷）奇点。Y(x)又被称为诺依曼函数（Neumann function），有时也记作N(x)。它和J(x)存在如下关系：若为整数（此时上式是0/0型未定式）则取右端的极限值。从前面对J(x)的定义可以知道，若不为整数时，定义Y是多余的（因为贝塞尔方程的两个线性无关解都已经用J函数表示出来了）。另一方面，若为整数，Y便可以和J构成贝塞尔方程的一个解系。与J函数类似，Y函数正负整数阶之间也存在如下关系：J(x)和Y(x)均为沿负实半轴割开的复平面内关于x的全纯函数。当为整数时，复平面内不存在贝塞尔函数的支点，所以J 和Y 均为x 的整函数。若将x 固定，则贝塞尔函数是的整函数。图3所示为0阶、1阶和2阶第二类贝塞尔函数Y(x)的曲线（ = 0,1,2）：汉开尔函数贝塞尔方程的另外一对重要的线性无关解称为汉开尔函数（Hankel functions）H(1)(x)和H(2)(x)，分别定义为：其中i 为虚数单位。以上的线性组合也成为第三类贝塞尔函数；它们描述了二维波动方程的内行柱面波解和外行柱面波解（行与在行动中同音）。利用前面推出的关系可将汉开尔函数表示成：若为整数，则须对等号右边取极限值。另外，无论是不是整数，下面的关系都成立：虚宗量的贝塞尔函数（修正贝塞尔函数）贝塞尔函数当宗量x 为复数时同样成立，并且当x 为纯虚数时能得到一类重要情形它们被称为第一类和第二类虚宗量的贝塞尔函数，或修正贝塞尔函数（有时还称为双曲型贝塞尔函数），定义为：以上形式保证了当宗量x 为实数时，函数值亦为实数。这两个函数构成了下列修正贝塞尔方程（与一般贝塞尔方程的差别仅在两个正负号）的一个相互线性无关的解系：修正贝塞尔函数与一般贝塞尔函数的差别在于：一般贝塞尔函数随实宗量是振荡型的，而修正贝塞尔函数I 和K则分别是指数增长和指数衰减型的。和第一类贝塞尔函数J一样，函数I当 0 时在x=0 点等于0，当=0时在x=0 点趋于有限值。类似地，K在x=0 点发散（趋于无穷）。图4-1 第一类修正贝塞尔函数I(x)对实自变量的曲线（ = 0,1,2）图4-2 第二类修正贝塞尔函数K(x)对实自变量的曲线（ = 0,1,2）复数宗量的贝塞尔函数之零值：J(x) = 0的解在-1的情况下都是实数；阶数-2-1的情况下，除了实数之外还有且仅有一对共轭的纯虚数解（G.N Watson 参考文献5）。球贝塞尔函数图5-1 第一类球贝塞尔函数jn(x)曲线（n = 0,1,2）图5-2 第二类球贝塞尔函数yn(x)曲线（n = 0,1,2）若使用分离变量法求解球坐标下的三维拉普拉斯方程，则可得到如下形式关于径向（r 方向）分量的常微分方程：关于上述方程的一对线性无关解称为球贝塞尔函数，分别用jn和yn表示（有时也记为nn）。这两个函数与一般贝塞尔函数Jn和Yn 存在关系：球贝塞尔函数也可写成：0阶第一类球贝塞尔函数j0(x)又称为sinc函数。头几阶整阶球贝塞尔函数的表达式分别为：第一类：-第二类：- -还可以依照前面构造汉开尔函数相同的步骤构造所谓 球汉开尔函数：事实上，所有半奇数阶贝塞尔函数都可以写成由三角函数组成的封闭形式的表达式，球贝塞尔函数也同样可以。特别地，对所有非负整数n，存在：而对实自变量x，hn(2)是上面hn(1)的复共轭（! 表示双阶乘）。由此我们可以通过得到h，再分离实部虚部，求出相应阶j 和h 的表达式，譬如j0(x) = sin(x)/x，y0(x) = -cos(x)/x，等等。黎卡提-贝塞尔函数黎卡提-贝塞尔函数（Riccati-Bessel functions）和球贝塞尔函数比较类似：该函数满足方程：这个方程以及相应的黎卡提-贝塞尔解是德国物理学家古斯塔夫米（Gustav Mie）于1908年研究电磁波在球状颗粒表面散射问题时提出的，后人将这种散射称为米氏散射（Mie scattering）。这个问题近几年的进展可参见文献 Du (2004)。后人有时会遵从德拜（Debye）在1909年的论文中的记法，用n,n 代替前面的Sn,Cn。渐近形式贝塞尔函数在非负时具有下面的渐近形式。当自变量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）相关，由阿贝尔恒等式（Abels 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/O177 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 0486604624. 5 G. N. Watson, A Treatise on the Theory of Bess