第五章 - 贝塞尔函数

n阶第一类贝塞尔函数Jn(x)

第二类贝塞尔函数，或称Neumann函数Yn(x) 第三类贝塞尔函数汉克尔（Hankel）函数，Hn(1)(x) 第一类变形的贝塞尔函数In(x)

开尔文函数（或称汤姆孙函数）n阶第一类开尔文（Kelvin）

第五章 贝塞尔函数

在第二章中，用分离变量法求解了一些定解问题。从§2.3可以看出，当我们采用极坐标系后，经过分离变量就会出现变系数的线性常微分方程。在那里，由于只考虑圆盘在稳恒状态下的温度分布，所

以得到了欧拉方程。如果不是考虑稳恒状态而是考虑瞬时状态，就会得到一种特殊类型的常微分方程。本章将通过在柱坐标系中对定解问题进行分离变量，引出在§2.6中曾经指出过的贝塞尔方程，并讨论这个方程解的一些性质。下面将看到，在一般情况下，贝塞尔方程的解不能用初等函数表出，从而就导入一类特殊函数，称为贝塞尔函数。贝塞尔函数具有一系列性质，在求解数学物理问题时主要是引用正交完备性。

§5.1 贝塞尔方程的引出

下面以圆盘的瞬时温度分布为例推导出贝塞尔方程。设有半径为R的薄圆盘，其侧面绝缘，若圆盘边界上的温度恒保持为零摄氏度，且初始温度为已知，求圆盘内瞬时温度分布规律。

这个问题可以归结为求解下述定解问题：

2??u?2u22?u22??t?a(?x2??y2),x?y?R,t?0, (5.1)??222?ut?0??(x,y),x?y?R, (5.2) ??ux2?y2?R2?0, (5.3)??用分离变量法解这个问题，先令

u(x,y,t)?V(x,y)T(t)

代入方程（5.1）得

?2V?2VVT??a(2?2)T

?x?y2或

?2V?2V?T??x2?y2???? (??0) a2TV由此得到下面关于函数T(t)和V(x,y)的方程

（5.4） T??a2?T?0

?2V?2V ?2??V?0 （5.5）2?x?y从（5.4）得

T(t)?Ae?a?t

2方程（5.5）称为亥姆霍兹（Helmholtz）方程。为了求出这个方程满足条件

Vx2?y2?R2?0 （5.6）

的非零解，引用平面上的极坐标系，将方程（5.5）与条件（5.6）写成极坐标形式得

??2V1?v1?2V???2??????2??2??V?0,??R,0???2?, (5.7) ??V???R?0,0???2?, (5.8) 再令 V(?,?)?P(?)?(?)， 代入（5.7）并分离变量可得

???(?)???(?)?0

（5.9）

?2P??(?)??P?(?)?(??2??)P(?)?0 （5.10）

由于u(x,y,t)是单值函数，所以V(x,y)也必是单值得，因此?(?)应该是以2?为周期的周期函数，这就决定了?只能等于如下的数：

0,12,22,?,n2,?

对应于?n?n2，有

?0(?)?a0（为常数） 2?n(?)?ancosn??bnsinn?,(n?1,2,?)

以?n?n2代入（5.10）得

?2P??(?)??P?(?)?(??2?n2)P(?)?0 （5.11）

这个方程与（2.93）相比，仅仅是两者的自变量和函数记号有差别，所以，它是n阶贝塞尔方程。

若再作代换

r???，

并记

F(r)?P(r?)，

则得

r2F??(r)?rF?(r)?(r2?n2)F(r)?0.

这是n阶贝塞尔方程最常见的形式。

由条件（5.8）及温度u是有限的，分别可得

??P(R)?0 （5.12） ?P(0)?????因此，原定解问题的最后解决就归结为求贝塞尔方程（5.11）在条件（5.12）下的特征值与特征函数（（5.12中第一个条件是在??R处的第一类边界条件，第二个条件是在??0处的自然边界条件，由于。在下一k(?)??2在??0处为零，所以在这一点应加自然边界条件）节先讨论方程（5.11）的解法，然后在§5.5中再回过头来讨论这个特征值问题。

§5.2 贝塞尔方程的求解

在上一节中，从解决圆盘的瞬时温度分布问题引出了贝塞尔方程，本节来讨论这个方程的解法。按惯例，仍以x表示自变量，以y表示未知函数，则n阶贝塞尔方程为

d2ydyx?x?(x2?n2)y?0 （5.13） 2dxdx2其中n为任意实数或复数。我们仅限于n为任意实数，且由于方程中的系数出现n2的项，所以在讨论时，不妨先假定n?0。

设方程（5.13）有一个级数解，其形式为

y?x(a0?a1x?a2x???akx??)??akxc?k，a0?0 （5.14）

c2kk?0?其中常数c和ak(k?0,1,2,?)可以通过把y和它的导数y?,y??代入（5.13）来确定。

将（5.14）及其导数代入（5.13）后得

?{[(c?k)(c?k?1)?(c?k)?(xk?0?2?n2)]akxc?k}?0

化简后写成

(c?n)a0x?[(c?1)?n]a1x22c22c?1??{[(c?k)2?n2]ak?ak?2}xc?k?0

k?2?要上式为恒等式，必须各个x幂的系数全为零，从而得到下列各式： 1°a0(c2?n2)?0； 2°a1[(c?1)2?n2]?0；

3°[(c?k)2?n2]ak?ak?2?0(k?2,3,?)。

由1°得c??n，代入2°得a1?0。先暂取c?n，代入3°得

4°ak??ak?2。

k(2n?k)因为a1?0，由4°知a1?a3?a5?a7???0，而a2,a4,a6,?都可以用a0表示，即

a2??a0，

2(2n?2)a0，

2?4(2n?2)(2n?4)?a0，

2?4?6(2n?2)(2n?4)(2n?6)a4?a6??

a2m?(?1)m?a02?4?6?2m(2n?2)(2n?4)?(2n?2m)m(?1)a022mm!(n?1)(n?2)?(n?m).

由此知（5.14）的一般项为

(?1)ma0xn?2m

22mm!(n?1)(n?2)?(n?m)让a0取一个确定的值，就得（5.13）得一个特解。a0是一个任意常数，把a0取作

a0?1 n2?(n?1)这样选取a0可使一般项系数中2的次数与x的次数相同，并可以运用下列恒等式：

(n?m)(n?m?1)?(n?2)(n?1)?(n?1)??(n?m?1)

使分母简化，从而使（5.14）中一般项的系数变成

a2m?(?1)m1 （5.15） n?2m2m!?(n?m?1)这样就比较整齐、简单了。

以（5.15）代入（5.14）得到（5.13）的一个特解

xn?2my1??(?1)n?2m(n?0)

2m!?(n?m?1)m?0?m用级数的比率判别法（或称达朗贝尔判别法）可以判定这个级数在整个数轴上收敛。这个无穷级数所确定的函数，称为n阶第一类贝塞尔函数。记作

xn?2m Jn(x)??(?1)n?2m(n?0) （5.16）

2m!?(n?m?1)m?0?m至此，就求出了贝塞尔方程的一个特解Jn(x)。

当n为正整数或零时，?(n?m?1)?(n?m)!，故有

xn?2m Jn(x)??(?1)n?2m(n?0,1,2,?) （5.17）

2m!(n?m)!m?0?m取c??n时，用同样的方法可得（5.13）的另一特解

x?n?2m J?n(x)??(?1)?n?2m(n?1,2,?) （5.18）

2m!?(?n?m?1)!m?0?m比较（5.16）式与（5.18）式可见，只要在（5.16）右端把n换成?n，即可得到（5.18）式。因此不论n式正数还是负数，总可以用（5.16）统一地表达第一类贝塞尔函数。

当n不为整数时，这两个特解Jn(x)与J?n(x)是线性无关的，由齐次线性常微分方程的通解的结构定理知道，（5.13）的通解为

y?AJn(x)?BJ?n(x) （5.19）

其中A,B为两个任意常数。

当然，在n不为整数的情况，方程（5.13）的通解除了可以写成（5.19）式以外还可以写成其它的形式，只要能够找到该方程另一个与Jn(x)线性无关的特解，它与Jn(x)就可构成（5.13）的通解，这样的特解是容易找到的。例如，在（5.19）中取A?cotn?,B??cscn?，则得到（5.13）的一个特解

Yn(x)?cotn?Jn(x)?cscn?J?n(x)（5.20） Jn(x)cosn??J?n(x)?(n?整数)sinn?显然，Yn(x)与Jn(x)是线性无关的，因此，（5.13）的通解可以写成

y?AJn(x)?BYn(x) （5.21）

由（5.20）式所确定的函数Yn(x)称为第二类贝塞尔函数，或称Neumann函数。

§5.3 当n为整数时贝塞尔方程的通解

上一节说明，当n不为整数时，贝塞尔方程（5.13）的通解由（5.19）或（5.21）式确定，当n为整数时，（5.13）的通解应该是什么样子呢？ 首先，我们证明当n为整数时，Jn(x)与J?n(x)是线性相关的。事实上，不妨设n为正整数N（这不失一般性，因n为负整数时，会得到同样的结果），这在（5.18）中，

1当m?0,1,2,?,(N?1)时

?(?N?m?1)均为零，这时级数从m?N起才开始出现非零项。于是（5.18）可以写成

x?N?2mJ?N(x)??(?1)?n?2m2m!?(?N?m?1)!m?N?mxNxN?2xN?4 ?(?1){N?N?2?N?4??}

2N!2(N?1)!2(N?2)!2!N ?(?1)NJN(x)即JN(x)与J?N(x)线性相关，这时JN(x)与J?N(x)已不能构成贝塞尔方程的通解了。为了求出贝塞尔方程的通解，还要求出一个与JN(x)线性无关的特解。

取哪一个特解？自然我们想到第二类贝塞尔函数。不过当n为整数时（5.20）的右端没有意义，要想把整数阶贝塞尔方程的通解也写成（5.21）的形式，必须先修改第二类贝塞尔函数的定义。在n为整数的情况，我们定义第二类贝塞尔函数为

Yn(x)?limJ?(x)cos???J??(x)(n为整数) （5.22） ??nsin??由于当n为整数时，J?n(x)?(?1)nJn(x)?cosn?Jn(x)，所以上式右端的极限为“”形式的不定型的极限，应用洛必达法则并经过冗长的推导，最后得

x(?1)m()2mm?12x212Y0(x)?J0(x)(ln?c)??? ?2?m?0(m!)2k?0k?1?00x1n?1(n?m?1)!?x?Yn(x)?Jn(x)(ln?c)?????2?m?0m!?2? （5.23） mxn?2m(?1)()n?m?1m?11?112 ??(???),(n?1,2,3,?)?m?0m!(n?m)!k?0k?1k?0k?12??????ln)?0.5572?，称为欧拉常数。 其中c?lim(1n??12131n?n?2m根据这个函数的定义，它确是贝塞尔方程的一个特解，而且与

Jn(x)是线性无关的（因为当x?0时，Jn(x)为有限值，而Yn(x)为无穷

大）。

综上所述，不论n是否为整数，贝塞尔方程（5.13）的通解都可表示为

y?AJn(x)?BYn(x)

其中A,B为任意常数，n为任意实数。 §5.4贝塞尔函数的递推公式

不同阶的贝塞尔函数之间不是彼此鼓孤立的，而是有一定的联系，本节来建立反映这种联系的递推公式。

先考虑零阶与一阶贝塞尔函数之间的关系。 在（5.17）中令n?0及n?1得

x2x4x6x2kkJ0(x)?1?2?4?6???(?1)2k?? 22222(2!)2(3!)2(k!)xx3x5x7x2k?1kJ1(x)??3?5?7???(?1)2k?1??

22?2!2?2!?3!2?3!?4!2?k!?(k?1)!取出第一个级数的第k?2项求导数，得

2k?1dx2k?2x2k?1k?1k(2k?2)xk (?1)??(?1)2k?2??(?1)2k?12k?222dx2[(k?1)!]2[(k?1)!]2k!(k?1)!这个式子正好是J1(x)中含x2k?1这一项的负值，且知J0(x)的第一项导数为零，故得关系式

dJ0(x)??J1(x) （5.24） dx将J1(x)乘以x并求导数，又得

ddx2x4x2k?2k[xJ1(x)]?[?3???(?1)2k?1??]dxdx22?2!2k!?(k?1)!x3x2k?1k?x?2???(?1)2k??222(k!)x2x2kk?x[1?2???(?1)2k??]22(k!)2

即

d[xJ1(x)]?xJ0(x) （5.25） dx以上结果可以推广，现将Jn(x)乘以xn求导数，得

dnd?x2n?2mm[xJn(x)]??(?1)2n?2mm!?(n?m?1)dxdxm?0xn?2m?1?x?(?1)n?2m?12m!?(n?m)m?0n?m

?xnJn?1(x)即

dn[xJn(x)]?xnJn?1(x) （5.26） dx同理可得

d?n[xJn(x)]??x?nJn?1(x) （5.27） dx将（5.26）和（5.27）两式左端的导数求出来，并经过化简，这分别得

?(x)?nJn(x)?xJn?1(x) xJn?(x)?nJn(x)??xJn?1(x). 及 xJn将这两式相减及相加，分别得到

Jn?1(x)?Jn?1(x)?2nJn(x) （5.28） x?(x) （5.29） Jn?1(x)?Jn?1(x)?2Jn以上几式就是贝塞尔函数的递推公式，它们在有关贝塞尔函数的的分析运算中非常有用。特别值得一提的是，应用（5.28）式可以用

较低阶的贝塞尔函数把较高阶的贝塞尔函数表示出来，因此如果我们已有零阶与一阶贝塞尔函数表，这利用此表和（5.28），即可计算任意正整数阶的贝塞尔函数的数值。

第二类贝塞尔函数也具有与第一类贝塞尔函数相同的递推公式

?dnn[xY(x)]?xYn?1(x)n?dx??d[x?nY(x)]??x?nY(x)?nn?1 （5.30） ?dx?2n?Yn?1(x)?Yn?1(x)?Yn(x)x???Yn?1(x)?Yn?1(x)?2Yn?(x)作为递推公式的一个应用，考虑半奇数阶的贝塞尔函数，现计算

J12(x)，J?12(x)。由（5.16）可得

?2m(?1)mx1J12(x)??()2

3m?0m!?(?m)22?而

31?3?5???(2m?1)11?3?5???(2m?1)?(?m)??()?? 22m?122m?1从而

2?(?1)m2m?12 J12(x)?x?sinx （5.31）??xm?0(2m?1)!?x同理，可求得

J?12(x)?2cosx （5.32） ?x利用递推公式（5.28）得到

J32(x)? ?1J12(x)?J?12(x)x21(?cosx?sinx)?xx1dsinxx?()?xdxx2232

?? ???x(321dsinx)()xdxx同理可得

1231dcosxJ?32(x)?J12(x)?J?12(x)?x2()()

x?xdxx一般而言，有

J1(x)?(?1)2n2n??2xn?x12n?12(1dnsinx)() xdxxJ1?(n?)2(x)??(1dcosx )() （5.33）

xdxx这里为了方便起见，采用了微分算子(1d1dn)，它是算子连续作用

xdxxdx1dn1d2sinx1d1dsinx)()?[()]，例如(千万不能把它与nnn次的缩写，

xdxxxdxxdxxxdx混为一谈。

从（5.33）可以看出，半奇数阶的贝塞尔函数都是初等函数。 §5.5函数展成贝塞尔函数的级数

利用贝塞尔求解数学物理方程的定解问题，最终要把已知函数按贝塞尔方程的特征函数系进行展开。这一节我们先要所明贝塞尔方程的特征函数系是什么样的函数系，然后证明这个特征函数系是一个正交系。

5.5.1 贝塞尔函数的零点

在§5.1中，已经将求解圆盘的温度分布问题通过分离变量法转化成贝塞尔方程的特征值问题：

?r2P??(r)?rP?(r)?(?r2?n2)P(r)?0,0?r?R,(5.34)??P(r)r?R?0, (5.35) ??P(0)???(自然边界条件). (5.36)方程（5.34）的通解为

P(r)?AJn(?r)?BYn(?r)，

由条件（5.36）可得B?0，即

P(r)?AJn(?r)

利用条件（5.35）得

Jn(?R)?0 （5.37）

这就说明，为了求出上述特征值问题的特征值?必须要计算Jn(x)的零点。Jn(x)有没有实的零点？若存在实的零点，一共有多少个？关于这些问题，有以下结论：

1°Jn(x)有无穷多个单重实零点，且这无穷多个零点在x轴上关于原点实对称分布的，因而Jn(x)必有无穷多个正的零点。

2°Jn(x)的零点与Jn?1(x)的零点是彼此相间分布的，即Jn(x)的任意两个相邻零点之间必存在一个且仅有一个Jn?1(x)的零点。

(n)(n)(n)3°以?m表示Jn(x)的非负零点（m?1,2,?），则?m???1m当m??时无限地接近于?，即Jn(x)几乎是以2?为周期的函数。J0(x)与J1(x)的图形见图5.1。

为了便于工程技术上的应用，贝塞尔函数的正零点的数值已被详

细计算出来，并列成表格。下表给出了Jn(x)(n?0,1,2,?,5)的前9个正

(n)零点?m(m?1,2,?,9)的近似值：

利用上述关于贝塞尔函数零点的结论，方程（5.37）的解为

(n)（m?1,2,?） ?R??m即

?n?((n)?mR)2（m?1,2,?） （5.38）

与这些特征值相对应的特征函数为

Pm(r)?Jn((n)?mRr)（m?1,2,?） （5.39）

5.5.2 贝塞尔函数的正交性

(n)??m?现在来讨论特征函数系?Jn(我们将r)? (m?1,2,?)的正交性，

R??要证明

?R0?0,当m?k,??? （5.40） rJn(r)Jn(r)dr??R22R22(n)(n)RRJn?1(?m),当m?k,?Jn?1(?m)??22(n)m(n)k(n)??m?由于贝塞尔函数系?Jn(r)? (m?1,2,?)是特征值问题（5.34～

R??5.36）的特征函数系，所以它的正交性由§2.6中的施图姆－刘维尔理论可以直接推出。不过因为在那里我们并没有就一般情况证明这个

结论，因此，我们在这里把贝塞尔函数系的正交性详细证明一下，而且这个证明方法是富有启发性的，完全可以类似的步骤来证明§2.6中的结论3。下一章将要讲到的勒让德多项式的正交性，也是施图姆－刘维尔理论的另一个具体例子。

下面就来证明（5.40）。为了书写方便，令

F1(r)?Jn((n)?mRr)，

F2(r)?Jn(?r)(?为任意参数)，

按定义，F1(r)，F2(r)分别满足

(n)?mddF1(r)n22[r]?[()r?]F1(r)?0 drdrRrddF2(r)n22[r]?[?r?]F2(r)?0 drdrr以F2(r)乘第一个方程减去以F1(r)乘第二个方程，然后对r从0到R积分得

[((n)?mRR0)??]?rF1(r)F2(r)dr??F2(r)0022RR??F1(r)ddF2(r)[r]dr?0drdrddF1(r)[r]drdrdr

即

[((n)?mR)??]?rF1(r)F2(r)dr?{r[F2(r)F1?(r)?F1(r)F2?(r)]}0?0 220RR由此可得

?R0rF1(r)F2(r)dr??R[F2(R)F1?(R)?F1(R)F2?(R)]((n)?m

R)2??2(n)因F1(R)?Jn(?m)?0，故上式可写成

?R0(n)(n)?(?m?mJn(?R)Jn) （5.41） rJn(r)Jn(?r)dr????(n)(n)?m?R()2??2(m)2??2(n)?mRF2(R)F1?(R)RR若取???k(n)R,k?m，则

(n)Jn(?R)?Jn(?k)?0，

从而（5.41）的右端为零，即（5.40）中第一个式子已得证。

为了证明（5.40）中第二个式子，在（5.41）两端令??00(n)?mR，

此时（5.41）右端的极限是“”形式的不定型的极限，利用洛必达法则计算这个极限得

?R0(n)(n)?(?m?(?R)?RR2??mJn)Jn(n)2?(?mrJ(r)dr?lim?[Jn)] (n)?mR?2?2??2nR(n)?m由递推公式

?(x)?nJn(x)?xJn?1(x) xJn?(x)?nJn(x)??xJn?1(x) xJn(n)及Jn(?m)?0可知

(n)(n)(n)?(?mJn)?Jn?1(?m)??Jn?1(?m)

从而?0RR22R22(n)(n)rJ(r)dr?Jn?1(?m)?Jn?1(?m)，这就是（5.40）中第二

R222n(n)?m个式子。通常把定积分

?R0rJ(2n(n)?mRr)dr

的正平方根，称为贝塞尔函数Jn((n)?mRr)的模。

利用§2.6中关于特征函数系的完备性可知，任意在[0,R]上具有

一阶连续导数及分段连续的二阶导数的函数f(r)，只要它在r?0处有界，在r?R处等于零，则它必能展开成如下形式的绝对且一致收敛的级数

f(r)??AmJn(m?1?(n)?mR r) （5.42）

?k(n)Rr)，并

为了确定这个展开式的系数Am，在（5.42）两端同乘以rJn(对r从0到R积分，由正交关系式（5.40）得

?即

R0rf(r)Jn(?k(n)Rr)dr?Ak?rJ(0R2n?k(n)Rr)dr

Ak??0R22(n)Jn?1(?k)21Rrf(r)Jn(?k(n)R（5.43） r)dr

下一节将通过例子说明贝塞尔函数在求解定解问题时的用法。 §5.6贝塞尔函数应用举例

下面举两个例子，说明用贝塞尔函数求解定解问题的全过程。 例1 设有半径为1的薄均匀圆盘，边界上温度为零摄氏度，初始时刻圆盘内温度分布为1?r2，其中r是圆盘内任一点的极半径，求圆盘内温度分布规律。

解 由于是在圆盘内求解问题，故采用极坐标系较为方便，并考虑到定解条件与?无关，所以温度分布只能是r,t的函数，于是根据问题的要求，即可归结为求解下列定解问题：

2??u1?u2?u?a(?),0?r?1,t?0, (5.44)??t2?rr?r???ur?1?0,t?0, (5.45) ?2u?1?r,0?r?1, (5.46)t?0???此外，由物理意义，还有条件u???，且当t???时，u?0。令

u(r,t)?F(r)T(t)

代入方程（5.44）得

1FT??a2(F???F?)T

r或

1???F?FT?r???? 2aTF由此得

r2F???rF???r2F?0 （5.47） T??a2?T?0 （5.48）

方程（5.48）得解为

T(t)?Ce?a?t

2因为t???时，u?0。所以?只能大于零，令???2，则

T(t)?Ce?a2?2t

此时方程（5.47）的通解为

F(r)?C1J0(?r)?C2Y0(?r)

由u(r,t)的有界性，可知C2?0，再由（5.45）得J0(?)?0，即?是J0(x)的零点。以?n(0)表示J0(x)的正零点，则

(0)???n(n?1,2,3,?)

综合以上结果可得

(0)F(r)?J0(?nr)

T(t)?Ce?a2(0)2(?n)t

从而

un(r,t)??Cnen?1?(0)2?a2(?n)t(0)J0(?nr)

由条件（5.46）得

(0)1?r??CnJ0(?nr)

2n?1?从而

Cn??2?(?)][J02(0)J12(?n)(0)n2?1010(0)(1?r2)rJ0(?nr)dr1

(0)(0)[?rJ0(?nr)dr??r3J0(?nr)dr]0(0)(0)(0)(0)(0)因d[(?nr)J1(?nr)]?(?nr)[J0(?nr)d(?nr)]，即

d[(0)rJ1(?nr)(0)?n(0)]?rJ0(?nr)dr

故得

(0)(0)rJ1(?nr)1J1(?n) rJ(?r)dr??(0)(0)?000?n?n1(0)n另外

?10rJ0(?r)dr??rd[03(0)n12(0)rJ1(?nr)(0)?n(0)J1(?nr)12]?r?(0)(0)0?n?n3?10(0)r2J1(?nr)dr?J1(?)?(0)n(0)n?2(0)2(?n)(0)r2J2(?nr)10?J1(?)?(0)n(0)n?2J(?)(?)(0)2n(0)2n

从而

(0)4J2(?n)Cn?(0)22(0)

(?n)J1(?n)所以，所求定解问题的解为

(0)2(0)24J2(?n)(0)u(r,t)??(0)22(0)J0(?nr)e?a(?n)t (5.49)

n?1(?n)J1(?n)?其中?n(0)是J0(r)的正零点。

例2 求下列定解问题：

2??2u1?u2?u?a(?),0?r?R,t?0, (5.50)??t22?rr?r???u?r?R?0,ur?0???t?0, (5.51) ??r??ur2,0?r?R, (5.52)?ut?0?0,t?0?1?2?tR? 解 用分离变量法来解，令u(r,t)?F(r)T(t)，采用例1类似的运算，可以得到

F(r)?C1J0(?r)?C2Y0(?r) （5.53） T(t)?C3cosa?t?C4sina?t （5.54）

由u(r,t)在r?0处的有界性，可知C2?0，即

F(r)?C1J0(?r) （5.55）

再根据边界条件（5.51）中第一式，得

?(?R)?0 F?(R)?C1?J0因C1?不能为零，故有

?(?R)?0 J0利用贝塞尔函数的递推公式（5.24）可得

J1(?R)?0

(1)(1)即?R是J1(x)?0的非负零点，以?1(1),?2,?,?n,?表示J1(x)的所有正零

点，又因J1(0)?0，所以

(1)??0及?R??n(n?1,2,?)（5.56）

当??0时，由（5.47），（5.48）及（5.51）中第二个条件可知，方程（5.50）有一个特解

u0(r,t)?C0?D0t

其中C0,D0是待定常数。

当??(1)?nR,(n?1,2,?)时，由方程（5.55）及（5.54）得

(1)?nFn(r)?J0(Rr)

(1)(1)a?na?nT(t)?Cncost?Dnsint

RR即（5.50）由特解

(1)(1)(1)a?na?n?nun(r,t)?(Cncost?Dnsint)J0(r)

RRR其中Cn,Dn是待定常数。

利用叠加原理可得原定解问题的解为

(1)(1)(1)a?na?n?nu(r,t)?C0?D0t??(Cncost?Dnsint)J0(r)

RRRn?1?代入条件（5.52）得

C0??CnJ0(n?1?(1)?nR r)?0 （5.57）

(1)?na?r2(1)D0??Dn?nJ0(r)?1?2 （5.58）

Rn?1RR由（5.57）得Cn?0(n?0,1,2,?)，在（5.58）两端同乘以r并对r在[0,R]上积分得

D0?1?R0rdr?R0r21(1?2)rdr? R2由（5.58）并利用下面的结果（见习题五第14题）：如果?n(1)是J1(x)的正零点，则

?得到

R0R2R22(1)(1)(1)rJ(r)dr?J0(?n)J1?(?n)?J0(?n)

R2220(1)?n(1)R?n2r2Dn?(1)2(1)?(1?2)rJ0(r)dr0a?nRJ0(?n)RR?4RJ2(?)4R??(1)32(1)(1)3(1)a(?n)J0(?n)a(?n)J2(?n)(1)n

所以最后得到定解问题的解为

(1)(1)a?n?nt4R?1u(r,t)??sintJ0(r)。 ?(1)3(1)2an?1(?n)J2(?n)RR§5.7贝塞尔函数的其他类型

由于解决某些工程问题的需要，本节引入另外三种形式的贝塞尔函数。

5.7.1 第三类贝塞尔函数

第三类贝塞尔函数有名汉克尔（Hankel）函数，它是由下列公式来定义的：

(1)Hn(x)?Jn(x)?jYn(x) (2)Hn(x)?Jn(x)?jYn(x)

其中j??1，由于汉克尔函数是Jn(x)与Yn(x)的线性组合，所以，同样也具有第一类贝塞尔函数相同的递推形式：

dn(i)(i)[xHn(x)]?xnHn?1(x)， dxd?n(i)(i)[xHn(x)]??x?nHn?1(x)， (i?1,2) dx2n(i)(i)(i)HnHn(x)， ?1(x)?Hn?1(x)?xd(i)(i)(i)Hn(x)?H(x)?2Hn(x). ?1n?1dx在下一节将看到这种函数当x很大时有比较简单的渐近公式。

5.7.2 虚宗量的贝塞尔函数

当我们在圆柱形域内求解定解问题，如果圆柱上下两底的边界条件都是齐次的，侧面的边界条件是非齐次时，就会遇到形如

1n2y???y??(1?2)y?0 （5.60）

xx的方程，它和贝塞尔方程只有一项的符号有差别，若令x??jt就可将这个方程化成贝塞尔方程，因为

x??jt,dx??jdt

dy?1dyd2yd2y?,2??2 dxjdtdxdt代入（5.60）得到

d2y1dyn2??(1?2)y?0 dt2tdtt 因此方程（5.60）的通解为

y?AJn(jx)?BYn(jx)

这里

xn?2mJn(jx)?j?n?2m

2m!?(n?m?1)m?0n?将上式乘以j?n后，我们就定义它为第一类虚宗量的贝塞尔函数或称

第一类变形的贝塞尔函数，并记作

xn?2m （5.61） In(x)?jJn(jx)??n?2mm!?(n?m?1)m?02?n?特别，

x2x4x6I0(x)?1?2?4???

22(2!)226(3!)2关于第二类虚宗量贝塞尔函数Kn(x)定义如下：

当n是非整数时

1?[I?n(x)?In(x)]； Kn(x)?2sinn?当n是整数时

1?[I??(x)?I?(x)].（5.62） Kn(x)?lim2??nsin??所以方程（5.60）的通解又可写为

y?AIn(x)?BKn(x)

其中A,B为任意常数。

In(x)与Kn(x)不存在实的零点，所以它们的图形不是振荡型曲线，

这一点与Jn(x)及Yn(x)不同。

5.7.3 开尔文函数（或称汤姆孙函数）

n阶第一类开尔文（Kelvin）函数有两种形式，它们分别被定义

为Jn(x?j)的实部与虚部，记作bernx和beinx。在应用这种函数时，主要是用零阶和一阶的。由于

J0(x?j)?J0(xjj)??(?1)mm?0??jjx1()(m!)222m???m?0xjm()2m2(m!)2xxj2k()4k?j2k?1()4k?22?2???22k?0[(2k)!]k?0[(2k?1)!]kx4kkx4k?2?(?1)()?(?1)()2?j2???22k?0[(2k)!]k?0[(2k?1)!]

所以ber0x和bei0x分别为

x(?1)k()4k2 （5.63）ber0x?Re[x?j]?? 2[(2k)!]k?0?kx4k?2(?1)()?2bei0x?Im[x?j]?? （5.64） 2k?0[(2k?1)!]用类似的方法可以得到一阶开尔文函数：

xxx()3()5()71xber1(x)??[?2?2?2??] （5.65）

221!2!2!3!3!4!xxx()3()5()71xbei1(x)??[??2?2?2??] （5.67）

221!2!2!3!3!4!§5.8 贝塞尔函数的渐近公式

在应用贝塞尔函数解决工程技术问题时，常常需要求出这些函数当自变量x取很大值时的函数值，如果按照级数展开式来计算这些特定值，就要求计算级数很多项的和，这样作非常麻烦，因此想到要用另外的函数来代替收敛很慢的贝塞尔函数的级数表达式，这个函数既要能逼近贝塞尔函数，又要能节约计算的时间。

我们为寻找一些便于计算的公式，引入所谓贝塞尔函数的渐近公式，在这里我们不打算讨论这些公式是如何引出的，因为这些公式的引出是比较复杂的。下面只举处在应用最常见的渐近公式。

当x的值很大时，

Jn(x)?21111 [?n(x)cos(x???n?)??n(x)sin(x???n?)] （5.67）

?x424221111[?n(x)sin(x???n?)??n(x)cos(x???n?)] （5.68） ?x4242Yn(x)?1?n?)2j(x??H(x)?e42[?n(x)?j?nx] （5.69）

?x(1)n1?n?)2?j(x??42H(x)?e[?n(x)?j?nx] （5.70）

?x(2)n其中

(4n2?12)(4n2?32)(4n2?12)(4n2?32)(4n2?52)(4n2?72)?n(x)?1???? 242!(8x)4!(8x)4n2?12(4n2?12)(4n2?32)(4n2?52)?n(x)????3 1!(8x)3!(8x)只要x的值很大，用这些渐近展开式右端的前面几项来计算左端函数的近似值，就能达到比较满意的精确度，例如，我们常取

Jn(x)?211211cos(x???n?)，Yn(x)?sin(x???n?) ?x42?x42H(1)nnn???)???)2j(x?12?j(x?1(2)4242e?e，Hn? （5.71） ?x?x来作近似计算。

本文来源：https://www.bwwdw.com/article/lbow.html