数理方程中与贝塞尔函数有关的问题据百度百科介绍：贝塞尔（1784——1846）是德国天文学家，数学家，天体测量学的奠基人。

20岁时发表了有关彗星轨道测量的论文。

1810年任新建的柯尼斯堡天文台台长，直至逝世。

1812年当选为柏林科学院院士。

贝塞尔的主要贡献在天文学，以《天文学基础》（1818）为标志发展了实验天文学 ，还编制基本星表 ，测定恒星视差， 预言伴星的存在，导出用于天文计算的贝塞尔公式，较精确地计算出岁差常数等几个天文常数值，还编制大气折射表和大气折射公式，以修正其对天文观测的影响。

他在数学研究中提出了贝塞尔函数，讨论了该函数的一系列性质及其求值方法，为解决物理学和天文学的有关问题提供了重要工具。

此外，他在大地测量学方面也做出一定贡献，提出贝塞尔地球椭球体等观点。

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

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

实际生活中受敲击的鼓面的振动是各阶类似振动形态的叠加一、贝塞尔方程与贝塞尔函数Bessel 方程是二阶线性变系数齐次常微分方程0)(22222=-++y v x dx dy x dxy d x 其中，v 是常数，称为Bessel 方程的阶（不一定是整数），可取任何实或复数。

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

数理方程教科书采用第一类Bessel 函数和第二类Bessel 函数的线性组合表示方程的标准解函数。

贝塞尔函数也被称为圆柱函数或圆柱谐波。

通常所说的贝塞尔函数是指第一类Bessel 函数m v m m v xm v m x J 20)2()1(!)1()(+∞=∑++-=Γ贝塞尔方程是在圆柱坐标或球坐标下使用分离变量法求解拉普拉斯方程和亥姆霍兹方程时得到的（在圆柱域问题中得到的是整阶形式；在球域问题中得到的是半奇数阶形式），因此贝塞尔函数在波动问题以及各种涉及有势场的问题中占有非常重要的地位，典型的问题有：在圆柱形波导中的电磁波传播问题；圆柱体中的热传导问题；圆形（或环形）薄膜的振动模态分析问题；在其他一些领域，贝塞尔函数也相当有用。

如在信号处理中的调频合成（FM synthesis ）或凯泽窗（Kaiser window ）的定义中，都要用到贝塞尔函数。

在教科书中Bessel 方程来源1． 在圆柱坐标系下解二维热传导方程；⎪⎪⎩⎪⎪⎨⎧=+=<+=><++=2222222222,0),,()0,,(0,),(R y x u Ry 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 Ta 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,01122222R VR V VV V ρρλθρρρρ再一次使用分离变量法，令)()(),(θρθρΘP V =，代入方程整理得μλρρρ=''-=+'+''ΘΘPPP P 22 由此得两个方程0=+''ΘΘμ，0)(22=-+'+''P P P μλρρρ第一个二阶常微分方程的通解为θμθμθsin cos )(21C C +=Θ引入周期边值条件)0()2(ΘΘ=π，得12cos =πμ。

所以固有值2n =μ，(n = 0，1，2，……)固有函数系为0021)(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)(22222=-++y n x dx dy x dxy 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 Ta T yyxx 2由此得两个方程0=++V V V yy xx λ，0)()(2=+''t T a t T λ第一个是圆域上Laplace 算子的固有值问题，与热传导问题类似可得整数阶贝塞尔方程0)(22222=-++y n x dx dy x dxy d x 3．在圆柱坐标系下解三维拉普拉斯方程或亥姆霍夫方程。

圆域上亥姆霍兹方程边值问题⎪⎩⎪⎨⎧≤≤=≤≤<<=+∂∂+∂∂+∂∂=πθπθρθρρρρρ20,020,0,011222222R VR V k VV V 用分离变量法，令)()(),(θρθρΘP V =，代入方程整理得μρρρ=''-=+'+''ΘΘPPk P P 222 由此得两个方程0=+''ΘΘμ，0)(222=-+'+''P k P P μρρρ第一个二阶常微分方程的通解为θμθμθsin cos )(21C C +=Θ引入周期边值条件)0()2(ΘΘ=π，得12cos =πμ。

所以固有值2n =μ，(n = 0，1，2，……)固有函数系为0021)(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)(22222=-++y n x dx dy x dxy d x 二、贝塞尔方程与欧拉方程比较欧拉方程0222=++y dx dy x dx y d x λ也是一类二阶线性变系数齐次常微分方程。

该方程的二阶导数项和一阶导数项表达式与贝塞尔方程相同。

不同的是，贝塞尔方程中函数项系数为变系数，欧拉方程中函数项系数为常数。

贝塞尔方程只能求出级数形式的解，即使是零阶贝塞尔方程02222=++y x dx dy x dxy d x 欧拉方程可以通过自变量变换成为线性常系数常微分方程。

作变换：)exp(t x =，即x t ln =，未知函数的导数为dtdyx dx dt dt dy dx dy 1== )(1)(11222222dt dydxy d x dt dy dx d x dt dy x dx y d -=+-=代入微分方程，得0)(22=++-y dt dydt dy dty d λ方程化简为：022=+y dtyd λ，该方程有初等函数表达式的通解。

三、贝塞尔函数与伽马函数 1．正整数阶贝塞尔函数贝塞尔函数的阶数v 不一定是整数。

引入伽马函数使表达式简化，但有一丝神秘m v m m v xm v m x J 20)2()1(!)1()(+∞=∑++-=Γ当阶数为正整数时，贝塞尔函数可写成m n m m n xm n m x J 20)2()!(!)1()(+∞=∑+-=零阶贝塞尔函数mm m x m x J 2020)2()!()1()(∑∞=-= 还有一种是积分形式（可用于数值计算实验）⎰-=πξξξπ20)sin cos(21)(d x v x J v2．负整数阶贝塞尔函数由于自变量为负值时，伽马函数的值趋于正无穷大，所以负整数阶贝塞尔函数m n m m n xm n m x J 20)2()1(!)1()(+-∞=-∑++--=Γ中对于m < n 的项为零，故m n n m m n xm n m x J 2)2()1(!)1()(+-∞=-∑++--=Γ令 k = m – n ，则有m = n + k 。

所以k n k k nk n k k n n x k n k x k k n x J 2020)2()1(!)1()1()2()1()!()1()1()(+∞=+∞=-∑∑++--=++--=ΓΓ对比 J n (x ) 的表达式，知J -n (x ) = (– 1)n J n (x )这说明两个整数阶贝塞尔函数线性相关。

3．伽马函数这一特殊函数以无穷积分的形式做为定义⎰+∞--=1)ex p()(dx x x s s Γ是正整数阶乘函数的推广。

其中，s 可以取正实数也可以取实部为正的复数。

几个简单性质如下：（1）．1)ex p()1(0=-=⎰+∞dx x Γ；（2）．)()1(n n n ΓΓ=+； 事实上⎰⎰∞-+∞+∞-+--=-=+0100)ex p()(ex p )ex p()1(dx x x n x x dx x x n n n n Γ)()ex p(01n n dx x x n n Γ=-=⎰∞-（3）．π=)2/1(Γ事实上⎰∞-=)ex p(1)2/1(dx x xΓ令2t x =，则上式化为概率积分π=-=-=⎰⎰+∞∞20)ex p(2)ex p(1)2/1(dt t dx x xΓ（4）．当伽马函数的自变量为负值时，无穷积分发散。

即⎰⎰⎰⎰+∞+∞+∞++---+-=-=-=-111111)exp(1)exp(1)exp(1)exp()(dxx x dx x xdx x xdx x xs s s s s Γ由于111111)exp(1ss s x sdx x dx x x-=≥-⎰⎰++∞+所以自变量为负值时，伽马函数的值趋于正无穷大。

四、贝塞尔函数与几个常用函数的台劳级数比较1．贝塞尔函数的级数收敛性 贝塞尔函数通常用级数表达式m n m m n xm n m x J 20)2()1(!)1()(+∞=∑++-=Γ利用交错级数的收敛判别法，用系数比值取极限0)1)(1(41lim )2()!1(4)1(!lim ||lim 222=+++=+++++=∞→∞→+∞→m n m m n m m n m a a m m m m m ΓΓ 所以，级数对任意自变量x 收敛。