贝塞尔函数 柱函数

2 2 2
成柱面问题． （由于是二维问题，即退化为极坐标） 设 u ( x, y, t ) = u ( r , j , t ) = T (t )U ( r , j ) ，对泛定方程分离变量（取 l = k ）得
2
T ¢¢ + k 2 a 2 T = 0 （14.1.2） ì ¢¢ 1 ¢ 1 ¢¢ 2 ï U r + r U r + r 2 U j + k U = 0 í ïU | = 0 î r =l
(14.2.3) （14.2.4） （14.2.5） （14.2.6）
= ( -1) n å ( -1) l
l = 0
所以
J - n ( x ) = ( -1) n J n ( x )
同理可证
J - n ( x ) = J n (- x )
因此有重要关系
J n ( - x) = ( -1) n J n ( x )
（14.2.1）
J -n ( x ) = å ( -1) k
k = 0
式中 G ( x) 是伽马函数．当 n = n 整数时，上述的级数实际上是从 k = n 的项开始，即
¥
J n ( x ) = å (-1) k
k = 0
1 x ( ) n + 2 k , ( n ³ 0) k !( n + k )! 2
14.3 贝塞尔函数的基本性质
14.3.1 贝塞尔函数的递推公式
由贝塞尔函数的级数表达式（14.2.1）容易推出
d Jn ( x) J ( x ) [ n ] = - v +1 v dx x x d v [ x J v ( x)] = x v J v -1 ( x ) dx
（14.1.4）
d 2 y d y + x + ( x 2 -n 2 ) y = 0 2 dx dx y ( k r ) |r =l = y ( kl ) = 0 . 边界条件为 x2
（14.1.5）
方程（14.1.5）称为n 阶贝塞尔微分方程．这里n 和 x 可以为任意数．
14.1.2
) 与 是我们应该注意到：当 n = n 整数时，有 J - n ( x ) = ( -1) J n ( x ) ，故上述解中的 J n ( x J - n ( x ) 是线性相关的，所以(14.1.6)成为通解必须是n ¹ 整数. （2）当n 取任意值时： N ( x ) ，这样贝塞尔方程的通解可表示为 定义第二类贝塞尔函数 n y ( x) = AJn ( x) + BNn ( x )
(1) ì Hn ( x) = Jn ( x) + iNn ( x ) í (2) îHn ( x ) = Jn ( x) - iNn ( x) （14.1.8） (1) ( 2) 并分别将 Hn , Hn 称为第一种和第二种汉克尔函数． 最后，总结n 阶贝塞尔方程的通解通常有下列三种形式： n ¹ 整数） （i） y ( x) = AJn ( x ) + BJ -n ( x) (
第十四章
贝塞尔函数 柱函数
贝塞尔函数( 也称为圆柱函数) 是现代科学技术领域中经常遇到的一类特殊函数． 1732 年伯努利研究直悬链的摆动问题，以及 1764 年欧拉研究拉紧圆膜的振动问题时，都涉及到 这类函数．1824 年德国数学家贝塞尔（F.W.贝塞尔, 1784~1846）在研究天文学问题时又遇到了这类函数，并首次系统地研究了这类函数．因此 人们称这类 函数为贝塞尔函数，并被广泛应用到数学、物理、光通信和其它科 学技术领域之中． 在用分离变量法一章介绍了拉普拉斯方程在柱坐标系下分离变量得 到了一种特殊类型的常微分 方程:贝塞尔方程．通过幂级数解法得到了另一类特殊函数，称为贝 塞尔函数．贝塞尔函数具有 一系列性质，在求解数学物理问题时主要是引用贝塞尔函数的正交 完备性．
14.2.2 第二类贝塞尔函数
第二类贝塞尔函数（又称为诺依曼函数或第二类柱函数）定义为
Nn ( x) =
cos( n π)Jn ( x) - J -n ( x ) sin( n π) J n ( x ) 线性无关
（14.2.7）
这样定义的的理由是它既满足贝塞尔方程，又与
J ( x ) 线性独立； （i）当n ¹ 整数时，显然它与 n
(1)
(2)
14.2 三类贝塞尔函数的表示式及性质
14.2.1 第一类贝塞尔函数的表示式
第一类贝塞尔函数（也可直接简称为贝塞尔函数，或第一类柱函数） Jn ( x ) 的级数表示 式为
¥
Jn ( x ) = å ( -1) k
k = 0 ¥
1 x n + 2 k ( ) k ! G(n + k + 1) 2 1 x -n + 2 k ( ) k !G ( -n + k + 1) 2
d - v [ x Z v ( x)] = - x - v Z v +1 ( x ) dx
把两式左端展开, 又可改写为 (14.3.3) (14.3.4)
v Z v ( x ) = - Z v +1 ( x ) x v Zn¢ + Z v ( x) = Z v -1 ( x ) x 从(14.3.5)和(14.3.6)消去 Zn 或消去 Zn¢ 可得 Zn ¢ ( x) ¢ ( x) Z v +1 ( x ) = Z v -1 ( x ) - 2 Z v 2 v Z v +1 ( x) = - Z v -1 ( x) + Z v ( x ) x
再令 U ( r , j ) = R ( r ) F (j ) ，得到
（14.1.2）
F ¢¢ + n 2 F = 0 (14.1.3) 2 2 2 2 r R¢¢ + r R¢ + (k r -n ) R = 0
令 k r = x, R ( r ) = y ( x ) 于是(14.1.5)得到
例 14.3.1 求
（14.3.7） （14.3.8）
ò xJ
2
( x )dx
【解】 根据公式 (14.3.8) Z v -1 ( x) - Z v +1 ( x) = = J 0 ( x) - 2J1
ò xJ
2
¢ ( x)dx = xJ1 ( x) - 2[ xJ1 ( x ) - ò J1 ( x)dx ( x )dx = ò xJ 0 ( x)dx - 2ò xJ1 ]
(14.2.2)
而
¥
J - n ( x ) = å ( -1) k
k = n ¥
1 x - n + 2 k ( ) k !G( - n + k + 1) 2 1 x ( ) n + 2 l , (l = k - n ) l ! G( n + l + 1) 2
(14.3.9)
d dJ m m [r ] + {[k (j m ) ]2 - 2 }J m (k (j m ) r ) = 0 dr dr r
将(14.3.9)乘以 用分部积分法得到
2
(14.3.10)
J m (k
( m ) j
r ) ，将(14.3.10)乘以 J m ( k r ) ，然后两式相减，再积分，利
(14.3.5) (14.3.6)
x ) 和 Z v ( x ) 推算 Z v +1 ( x ) 的递推公式. 即为从 Z v -1 ( 上式也可以写成为 v Z v -1 ( x) + Z v +1 ( x ) = 2 Z v ( x ) x Z v -1 ( x) - Z v +1 ( x) = 2 Zn¢ ( x) ) 统称为柱函数． 任一满足一组递推关系的函数 Z v ( x
(14.3.1) (14.3.2)
) 和汉克尔函数也应该满足 以上两式都是贝塞尔函数的线性关系式. 诺伊曼函数 Nv ( x 上述递推关系．
) 代表 v 阶的第一或第二或第三类函数,总是有 若用 Z v ( x
d v [ x Z v ( x)] = x v Z v -1 ( x ) dx
( m ) i
{[ki( m ) ]2 - [k (j m ) ]2 }ò J m (ki( m ) r )J m (k (j m ) r ) r d r
0
r 0
故当
ki( m )
d d r = [ r J m (ki( m ) r ) J m (k (j m ) r ) - r J m (k (j m ) r ) J m (ki ( m ) r )] |0 = 0 dr d r ¹ k (j m ) 时
Hn ( x ) ，又
可以取任意数） （ii） y ( x) = AJn ( x ) + BNn ( x) (n
可以取任意数） （iii） y ( x) = AHn ( x ) + BHn ( x) (n 注明：第一、第二、第三类贝塞尔函数分别称为贝塞尔函数、诺依曼函数、汉克尔函数， 还可以分别简称为第一、第二、第三类柱函数．
（ii）当n = n （整数）时，可以用洛必达法则求极限的办法来证明它也与 独立，且其结果与用级数法寻找的第二个线性独立的解一致．
Jn ( x ) 线性
图 20.1
14.2.3 第三类贝塞尔函数
第三类贝塞尔函数（又称为汉克尔（Hankel）函数，或第三类柱函数） ，根据其定 义式即
(1) ì Hn ( x) = Jn ( x ) + i × Nn ( x ) í (2) îHn ( x) = Jn ( x ) - i × Nn ( x )
14.1 贝塞尔方程及其解
14.1.1 贝塞尔方程
拉普拉斯方程在柱坐标系下的分离变量得出了一般的贝塞尔方程，由于贝塞尔方程的 普遍性，我们还能从其它典型的数学物理定解问题来导出贝塞尔方程的一般形式． 考虑固定边界的圆膜振动，可以归结为下述定解问题
ìutt = a 2 (u xx + u yy ) (0 £ x 2 + y 2 < l 2 , t > 0) ï (t ³ 0) ï u |x + y =l = 0 í ) ï u ( x, y , t ) |t = 0 = j ( x, y ï u ( x, y, t ) | = y ( x, y ) î t t = 0 （14.1.1） 其中 l 为已知正数，j ( x, y ),y ( x, y ) 为已知函数．这个定解问题因宜于使用柱坐标，从而构
= xJ1 ( x) - 2[ xJ1 ( x) + ò J 0¢ ( x )dx] = - xJ1 ( x ) - 2J 0 ( x ) + c 14.3.2 贝塞尔函数正交性和模
1．正交性 对应不同本征值的本征函数分别满足
2 d dJ m m ( m) 2 [r ] + {[ki ] - 2 }J m (ki ( m ) r ) = 0 dr dr r
贝塞尔方程的解
通过数学物理方程的幂级数求解方法可以得出结论： （1）当n ¹ 整数时，贝塞尔方程(14.1.6)的通解为
y ( x) = AJn ( x) + BJ -n ( x)
n
（14.1.6）
其中 A, B 为任意常数， Jn ( x ) 定义为n 阶第一类贝塞尔函数（简称为贝塞尔函数） ．但
（14.1.7）
n ( x) ，本书之所以选取 Nn ( x ) ，是因为它又称为诺依曼函 （第二类贝塞尔函数也可写成 Y 数，取第一大写字母）下面我们会看到不管n 是否为整数，上式均成立． (3) 当n 取任意值时：
实际上由第一、二类贝塞尔函数还可以构成线性独立的第三类贝塞尔函数 称为汉克尔函数．
0
ò
0
r0
J m (ki( m ) r )J m (k (j m ) r )r dr = 0
(14.3.11)
2．贝塞尔函数的模
N
( m ) n
2 r 2 1 m ( m) 2 ( m ) 2 [ Nn ] = ( r 02 + 0 )[J m (k n r 0 )] 2 ln ln H