贝塞尔函数释疑

数理方程中与贝塞尔函数有关的问题

据百度百科介绍：

贝塞尔（1784——1846）是德国天文学家，数学家，天体测量学的奠基人。20岁时发表了有关彗星轨道测量的论文。1810年任新建的柯尼斯堡天文台台长，直至逝世。1812年当选为柏林科学院院士。贝塞尔的主要贡献在天文学，以《天文学基础》（1818）为标志发展了实验天文学 ，还编制基本星表 ，测定恒星视差， 预言伴星的存在，导出用于天文计算的贝塞尔公式，较精确地计算出岁差常数等几个天文常数值，还编制大气折射表和大气折射公式，以修正其对天文观测的影响。他在数学研究中提出了贝塞尔函数，讨论了该函数的一系列性质及其求值方法，为解决物理学和天文学的有关问题提供了重要工具。此外，他在大地测量学方面也做出一定贡献，提出贝塞尔地球椭球体等观点。（图片来自维基百科）

一、 贝塞尔方程与贝塞尔函数 二、 贝塞尔方程与欧拉方程比较 三、 贝塞尔函数与伽马函数 四、 贝塞尔函数与几个常用函数的台劳级数比较 右图来自网页“维基百科——自由的百科全书”中贝塞尔

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

一、贝塞尔方程与贝塞尔函数

Bessel 方程是二阶线性变系数齐次常微分方程

0)(222

22

=-++y v x dx dy x dx

y d x 其中，v 是常数，称为Bessel 方程的阶（不一定是整数），可取任何实或复数。该方程

的解无法用初等函数表现。数理方程教科书采用第一类Bessel 函数和第二类Bessel 函数的线性组合表示方程的标准解函数。贝塞尔函数也被称为圆柱函数或圆柱谐波。通常所说的贝塞尔函数是指第一类Bessel 函数

m v m m v x

m v m x J 20)2

()1(!)1()(+∞

=∑++-=Γ

贝塞尔方程是在圆柱坐标或球坐标下使用分离变量法求解拉普拉斯方程和亥姆霍兹方程时得到的（在圆柱域问题中得到的是整阶形式；在球域问题中得到的是半奇数阶形式），因此贝塞尔函数在波动问题以及各种涉及有势场的问题中占有非常重要的地位，典型的问题有：在圆柱形波导中的电磁波传播问题；圆柱体中的热传导问题；圆形（或环形）薄膜的振动模态分析问题；在其他一些领域，贝塞尔函数也相当有用。如在信号处理中的调频合成（FM synthesis ）或凯泽窗（Kaiser window ）的定义中，都要用到贝塞尔函数。

在教科书中Bessel 方程来源

1． 在圆柱坐标系下解二维热传导方程；

⎪⎪

⎩

⎪⎪⎨⎧=+=<+=><++=2222

222222,0),,()0,,(0,),(R y x u R

y x y x y x u t R y x u u a u yy xx t ϕ 用分离变量法，令u (x ，y ，t ) = V (x ，y )T (t )，代入方程整得

λ-=+='V V V T

a T yy xx 2

由此得两个方程

0)()(2=+'t T a t T λ，0=++V V V yy xx λ

其中，一阶常微分方程的通解为

)ex p()(2t a A t T λ-=

而另一个是圆域上Laplace 算子的固有值问题，在极坐标系下

⎪⎩⎪⎨⎧=<<=+∂∂+∂∂+∂∂=,

0,0,0112

2222R V

R V V

V V ρρλθρρρρ

再一次使用分离变量法，令)()(),(θρθρΘP V =，代入方程整理得

μλρρρ='

'-=

+'+''Θ

ΘP

P

P P 22 由此得两个方程

0=+''ΘΘμ，0)(22=-+'+''P P P μλρρρ

第一个二阶常微分方程的通解为

θμθμθsin cos )(21C C +=Θ

引入周期边值条件)0()2(ΘΘ=π，得12cos =πμ。所以固有值

2n =μ，(n = 0，1，2，……)

固有函数系为

002

1

)(a =

θΘ，θθθn b n a n n n sin cos )(+=Θ，(n = 1，2，……) 将固有值代入第二个常微分方程，得

0)(222=-+'+''P n P P λρρρ

令ρλ=

x ，)/()(λx P x y =，则方程转化为标准的整数阶贝塞尔方程

0)(2

22

2

2=-++y n x dx dy x dx

y d x 2． 圆柱坐标系下解二维波动方程；

⎪⎪⎩

⎪⎪⎨⎧>=+=<+==><++=0,,0),,()0,,(),,()0,,(0

,),(2222222222t R y x u R y x y x y x u y x y x u t R y x u u a u t yy xx tt ψϕ 用分离变量法，令u (x ，y ，t ) = V (x ，y )T (t )，代入方程整得

λ-=+=''V V V T

a T yy

xx 2

由此得两个方程

0=++V V V yy xx λ，0)()(2=+''t T a t T λ

第一个是圆域上Laplace 算子的固有值问题，与热传导问题类似可得整数阶贝塞尔方程

0)(2

22

22

=-++y n x dx dy x dx

y d x 3．在圆柱坐标系下解三维拉普拉斯方程或亥姆霍夫方程。

圆域上亥姆霍兹方程边值问题

⎪⎩⎪⎨⎧≤≤=≤≤<<=+∂∂+∂∂+∂∂=π

θπθρθρρρρ

ρ20,020,0,01122

2222R V

R V k V

V V 用分离变量法，令)()(),(θρθρΘP V =，代入方程整理得

μρρρ='

'-=

+'+''Θ

ΘP

P

k P P 222 由此得两个方程

0=+''ΘΘμ，0)(222=-+'+''P k P P μρρρ

第一个二阶常微分方程的通解为

θμθμθsin cos )(21C C +=Θ

引入周期边值条件)0()2(ΘΘ=π，得12

cos =πμ。所以固有值

2n =μ，(n = 0，1，2，……)

固有函数系为

002

1

)(a =

θΘ，θθθn b n a n n n sin cos )(+=Θ，(n = 1，2，……) 将固有值代入第二个常微分方程，得

0)(2222=-+'+''P n k P P ρρρ

令ρk x =，)/()(k x P x y =，则方程转化为标准的整数阶贝塞尔方程

0)(2

22

2

2=-++y n x dx dy x dx

y d x 二、贝塞尔方程与欧拉方程比较

欧拉方程

02

22

=++y dx dy x dx y d x λ

也是一类二阶线性变系数齐次常微分方程。该方程的二阶导数项和一阶导数项表达式与

贝塞尔方程相同。

不同的是，贝塞尔方程中函数项系数为变系数，欧拉方程中函数项系数为常数。 贝塞尔方程只能求出级数形式的解，即使是零阶贝塞尔方程

022

22

=++y x dx dy x dx

y d x 欧拉方程可以通过自变量变换成为线性常系数常微分方程。作变换：)exp(t x =，即

x t ln =，未知函数的导数为