《第五章-贝塞尔函数》由会员分享，可在线阅读，更多相关《第五章-贝塞尔函数（32页珍藏版）》请在金锄头文库上搜索。
1、n阶第一类贝塞尔函数Jn（x）第二类贝塞尔函数，或称Neumann0数Yn（x）第三类贝塞尔函数汉克尔（Hankel）函数，H：（x）第一类变形的贝塞尔函数In（x）开尔文函数（或称汤姆孙函数）n阶第一类开尔文（Kelvin）第五章贝塞尔函数f2盯co先0K=27iin，Jn(p)cosny,Jo*2厅ginne口鼠。一的=2亓，4(p)sinny.Jo在第二章中，用分离变量法求解了一些定解问题。从2.3可以看出，当我们采用极坐标系后，经过分离变量就会出现变系数的线性常微分方程。在那里，由于只考虑圆盘在稳恒状态下的温度分布，所以得到了欧拉方程。如果不是考虑稳恒状态而是考虑瞬时状态，就会得到一种
2、特殊类型的常微分方程。本章将通过在柱坐标系中对定解问题进行分离变量，引出在2.6中曾经指出过的贝塞尔方程，并讨论这个方程解的一些性质。下面将看到，在一般情况下，贝塞尔方程的解不能用初等函数表出，从而就导入一类特殊函数，称为贝塞尔函数。贝塞尔函数具有一系列性质，在求解数学物理问题时主要是引用正交完备性。5.1贝塞尔方程的引出下面以圆盘的瞬时温度分布为例推导出贝塞尔方程。设有半径为R的薄圆盘，其侧面绝缘，若圆盘边界上的温度恒保持为零摄氏度，且初始温度为已知，求圆盘内瞬时温度分布规律。这个问题可以归结为求解下述定解问题：22u2zuU222,a(),xyR,t0,(5.1)txy(5.2)(5.3)
3、V(x, y)T(t)2V2V 22c2ut0(x,y),xyR,ux2y2R20，用分离变量法解这个问题，先令u(x,y,t)代入方程(5.1)得VTa2(22xy精品资料2V2vTx2a2TV(0)由此得到下面关于函数T(t)和V(x,y)的方程2T a2 T 0(5.(4)从(5.4)得2V2V V 0 xy(5.(5)T(t) Ae a2 t方程(5.5)称为亥姆霍兹(Helmholtz)方程。为了求出这个方程满足条件22c20x y R(5.6 )的非零解，引用平面上的极坐标系，将方程(5.(6)(5.6 )写成极坐标形式得2V1 v-22V0,R,02 , (5.7)r 0,0(5
4、.8)再令V(P(代入(5.7)并分离变量可得(5.9)2P()P()()P()(5.10)应该是由于u(x,y,t)是单值函数，所以V(x,y)也必是单值得，因此(以2为周期的周期函数，这就决定了只能等于如下的数:2c2.2.0,1,2,L,n,L对应于n2,有 5.2贝塞尔方程的求解0()包(为常数)2n()ancosnbnsinn,(n1,2,L)以nn2代入(5.10)得2P()P()(2n2)P()0(5.11)这个方程与(2.93)相比，仅仅是两者的自变量和函数记号有差别，所以，它是n阶贝塞尔方程。若再作代换r1，并记rF(r)P(.,则得2一一22一r2F(r)rF(r)(r2n
5、2)F(r)0.这是n阶贝塞尔方程最常见的形式。P(R) 0P(0)由条件(5.8)及温度u是有限的，分别可得(5.12)因此，原定解问题的最后解决就归结为求贝塞尔方程(5.11)在条件(5.12)下的特征值与特征函数(5.12中第一个条件是在R处的第一类边界条件，第二个条件是在0处的自然边界条件，由于k()2在0处为零，所以在这一点应加自然边界条件)。在下一节先讨论方程(5.11)的解法，然后在5.5中再回过头来讨论这个特征值问题。精品资料在上一节中，从解决圆盘的瞬时温度分布问题引出了贝塞尔方程，本节来讨论这个方程的解法。按惯例，仍以x表示自变量，以y表示未知函数，则n阶贝塞尔方程为x2-y
6、x(x2n2)y0(5.13)dxdx其中n为任意实数或复数。我们仅限于n为任意实数，且由于方程中的系数出现n2的项，所以在讨论时，不妨先假定n0。设方程(5.13)有一个级数解，其形式为yxc(a0a1xa2x2LakxkL)akxck,a00(5.14)k0其中常数c和ak(k0,1,2,L)可以通过把y和它的导数y,y代入(5.13)来确定。将(5.14)及其导数代入(5.13)后得(ck)(ck1)(ck)(x2n2)akxck0k0化简后写成2_2c2_2c12_2ck(cn)ax(c1)nax(ck)nakak2x0k2要上式为恒等式，必须各个x哥的系数全为零，从而得到下列各式：1
7、a0(c2n2)0;2a1(c1)2n20；3(ck)2n2akak20(k2,3,L)。由1得cn,代入20得a10。先暂取cn,代入3。得4akak2。k(2nk)因为a10,由4知a1a3aa7L0,而a2,a4,a6,L都可以用a()表示，即a0a22(2n2)，aa0a4，2g4(2n2)(2n4)aa0a62g4g6(2n2)(2n4)(2n6)a2m(1)maO294g6L2m(2n2)(2n4)L(2n2m)m(1)a022mm!(n1)(n2)L(nm)由此知(5.14)的一般项为mn2m(1)a0X22mm!(n1)(n2)L(nm)a0是一个任意常数，让a0取一个确定的值
8、，就得(5.13)得一个特解。把a。取作a012n (n 1)这样选取a0可使一般项系数中2的次数与x的次数相同，并可以运用F列恒等式:(nm)(nm1)L(n2)(n1)(n1)(nm1)使分母简化，从而使(5.14)中一般项的系数变成(5.15)ma2m ( 1)1cn 2 m ；7772 m! (n m 1)这样就比较整齐、简单了。以(5.15 )代入(5.14)得到(5.13)的一个特解my1( 1) -m 02n 2mxn 2m ，/ m! (n m(n 0)1)用级数的比率判别法(或称达朗贝尔判别法)可以判定这个级数在整个数轴上收敛。这个无穷级数所确定的函数,称为n阶第一类贝塞尔精
9、品资料函数。记作Jn(x)( 1)mm 02n 2m xm! (n m 1)(n 0)(5.(16)至此，就求出了贝塞尔方程的一个特解Jn(x)。当n为正整数或零时，(n m 1) (n m)!,故有n 2mxJn(x)( 1)m2m-(n 0,1,2,L )n vimm 02 m!(n m)!(5.(17)取c n时，用同样的方法可得(5.13)的另一特解n 2mJn(x) m0( 1)m2n2mmx(n m 1)!(n 1,2,L(5.(18)比较(5.16)式与(5.18)式可见，只要在(5.16)右端把n换成n ,即可得到(5.18)式。因此不论n式正数还是负数，总可以用(5.16)统
10、一地表达第一类贝塞尔函数。当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)中取Acotn,Bcscn,则得到(5.13)的一个特解Yn(x)cotnJn(x)cscnJn(x)Jn(X)8snJ)(n整数)(
11、5.20)sinn显然，Yn(x)与Jn(x)是线性无关的，因此，(5.13)的通解可以写成yAJn(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)与Jn(x)是线性相关的。事实上，不妨设n为正整数N(这不失一般性，因n为负整数时，会得到同样的结果)，这在(5.18)中，1当m0,1,2,L,(N1)时(Nm1)均为
12、零，这时级数从mN起才开始出现非零项。于是(5.18)可以写成N 2mJ N (x)(1)mn 2m Xm n 2 m! ( N m 1)!NN 2( 1)NX2NN! 2 (N 1)!(1)NJN(x)N 4X N 4z2 (N 2)!2!L即Jn(x)与J N(X)线性相关,这时Jn(x)与J N(X)已不能构成贝塞尔方程的通解了。为了求出贝塞尔方程的通解，还要求出一个与Jn(x)线性无关的特解取哪一个特解？自然我们想到第二类贝塞尔函数。不过当n为整数时(5.20)的右端没有意义，要想把整数阶贝塞尔方程的通解也写成(5.21)的形式，必须先修改第二类贝塞尔函数的定义。在 n为整数的情况，我
13、们定义第二类贝塞尔函数为Yn(x) lim J (x)CosJ-(n为整数)(5.22 )n sin由于当n为整数时，J n(x) ( l)nJn(x) cosn Jn(x),所以上式右端的极限为“0 ”形式的不定型的极限，应用洛必达法则并经过冗长的推导，0最后得Y0(X)2X一J0(x)(lna c)m / X、2m2( 1) (2) m1 12m 0(m!) k 0 k 1K(X)2 x-Jn(x)(ln-、1 n 1 (n m 1)! xc)-m 0 m! 2n 2mmzX、n2m(15m!(nm)!(5.23)其中c lim(1n11L2ln)0.5572L，称为欧拉常数。根据这个函数
14、的定义，它确是贝塞尔方程的一个特解，而且与Jn(X)是线性无关的(因为当X0时，Jn(X)为有限值，而Yn(X)为无穷大)。综上所述，不论n是否为整数，贝塞尔方程(5.13)的通解都可表不为yAJn(x)BYn(x)其中A,B为任意常数，n为任意实数。 5.4 塞尔函数的递推公式不同阶的贝塞尔函数之间不是彼此鼓孤立的,而是有一定的联系,本节来建立反映这种联系的递推公式。先考虑零阶与一阶贝塞尔函数之间的关系在(5.17 )中令n 0及n 1得246J (x) 1 xx2224(2!) 226(3!)2(1)k2kx L2 (k!)2Ji(x)223g2!25 g2!g3!27g3!g4!2k 1L ( 1)kTTL22k1gk!g(k 1)!取出第一个级数的第k 2项求导数,
