Bessel函数在物理学中应用.docVIP

Bessel函数在物理学中应用 　　摘要：本文首先介绍了贝塞尔函数的来源，然后介绍了其在物理学中的应用。贝塞尔函数在物理学中的应用是很广泛的，笔者主要以贝塞尔函数在热传导问题、量子力学和电动力学中的应用为例来说明其应用。 　　关键词：贝塞尔函数；热传导；球方势阱；圆柱形导体；圆柱形波导 　　中图分类号：G633.7 文献标识码：A 文章编号：1992-7711（2014）01-0119 　　一、引言 　　贝塞尔函数的几个正整数阶特例早在18世纪中叶就由瑞士数学家丹尼尔?伯努利在研究悬链振动时提出，丹尼尔的叔叔雅各布?伯努利，欧拉、拉格朗日等数学大师对贝塞尔函数的研究做出贡献。1817年，德国数学家贝塞尔在研究开普勒提出的三体引力系统的运动问题时，第一次系统地提出了贝塞尔函数的总体理论框架，后人以他的名字来命名了这种函数，现在贝塞尔函数广泛地应用于物理学中。贝塞尔方程是在柱坐标或球坐标下使用分离变量法求解拉普拉斯方程和亥姆霍兹方程时得到的（在圆柱域问题中得到的是整阶形式 α=n；在球形域问题中得到的是半奇数阶形式 α=n+1/2），因此贝塞尔函数在波动问题以及各种涉及有势场的问题中占有非常重要的地位，最典型的问题有： 　　在圆柱形波导中的电磁波传播问题；圆柱体中的热传导问题；圆形（或环形）薄膜的振动模态分析问题。在其他一些领域，贝塞尔函数也相当有用。譬如在信号处理中的调频合成（FM synthesis）或凯泽窗（Kaiser window）的定义中，都要用到贝塞尔函数。 　　二、贝塞尔方程和贝塞尔函数 　　1. 贝塞尔方程的来源 　　贝塞尔方程是在柱坐标或球坐标下使用分离变量法求解拉普拉斯方程和亥姆霍兹方程时得到的。柱坐标系中拉普拉斯算符的表达式为△=■■（ρ■）+■■+■，令U（ρ，φ，Z）=R（ρ）φ（φ）Z（z），通过对拉普拉斯方程 2 U=0分离变量得到关于变量ρ的方程： 　　ρ2R″（ρ）+ρR′（ρ）+[（λ-μ）ρ2-n2]R（ρ）=0 （1）或 　　ρ2R″（ρ）+ρR′（ρ）+[（-μ）ρ2-n2]R（ρ）=0 （2） 　　其中n=0，1，2……，若（λ-μ）或（-μ）≥0则k2=（λ-μ）或（-μ）此时上述二方程均变为： 　　ρ2R″（ρ）+ρR′（ρ）+[k2ρ2-n2]R（ρ）=0 （3） 　　称为n阶贝塞尔方程，令x=kρ，y（x）=R（ρ）则n阶贝塞尔方程又可表示为： 　　x2y″（x）+xy′（x）+（x2-n2）=0 （4） 　　（3）和（4）是整数阶贝塞尔方程，若将之进行推广即将n推广为实数υ则得到任意阶的贝塞尔方程： 　　x2y″（x）+xy′（x）+（x2-υ2）=0 （5） 　　2. 贝塞尔方程的通解 　　（1）υ和-υ阶的Bessel函数――第一类柱函数， 　　Jυ（x）=■■（■）2k+υ，Jυ（x）=■■（■）2k-υ均为Bessel方程的特解，当υ≠n（n=0，±1，±2……）时Bessel方程的通解是yc=CυJυ（x）+D-υ（x），当υ=n时J-n（x）=（-1）nJn（x）。 　　（2）第二类柱函数――Neumann函数 　　定义Nυ（x）=■为第二类柱函数，此时Bessel方程的通解表示为：yc=AυJυ（x）+BυNυ（x）。 　　（3）第三类柱函数――Hankel函数 　　定义Hυ（1）（x）=Jυ（x）+iNυ（x），Hυ（2）=Jυ（x）-iNυ（x）为第三类柱函数，于是贝塞尔函数的通解又可表示为：y（x）=C1Hυ（1）（x）+C2Hυ（2）（x）。 　　3. 球贝塞尔方程的来源 　　利用球坐标系拉普拉斯算符的表达式，可得球坐标系亥姆霍兹方程的表达式： 　　■■（r2■）+■■（sinθ■）+■■+k2υ=0 （6） 　　把变数r跟变数θ，φ分离开来，以υ（r，θ，φ）=R（r）Y（θ，φ）代入式（6），用■遍乘各项并适当移项得： 　　■■（r2■）+k2r2=■■（sinθ■）-■■ （7） 　　左边是r的函数，右边是θ，φ的函数，两边相等是不可能的，除非两边实际上是同一个常数，通常把这个常数记作l（l+1）。这就分解为两个方程： 　　■■（sinθ■）+■■+l（l+1）R=0 （8）■（r2■）+[k2r2-l（l+1）]R=0 （9） 　　其中式（9）亦即： 　　r2■+2r■+[k2r2-l（l+1）]R=0 （10） 　　这叫做l阶球贝塞尔方程。这是因为对于k0，可以把自变数r和函数R（r）分别换作x和x=krR（r）=■ ，则方程（10）式成为： 　　x2■+x■+[x2-（l+■）2]y=0 （11） 　　这正是l+■阶的贝塞尔方程。 　　4. 球贝塞尔方程的解 　　若