贝塞尔函数 柱函数

合集下载

3.4 特殊函数及其应用(柱函数)

3.4 特殊函数及其应用(柱函数)

阶贝塞尔方程
二 阶贝塞尔方程的解
• x 2 y xy ( x 2 v 2 ) y 0 两个线性无关的解:
y( x) C1 J ( x) C2 J ( x),
v 整数或半整数
• x 2 y xy ( x 2 m 2 ) y 0 两个线性无关的解:
d x n J n ( x) x n J n 1 ( x) dx d n x J n ( x) x n J n1 ( x) dx 2n J n 1 ( x) J n 1 ( x) J n ( x) x J n1 ( x) J n1 ( x) 2 J n ( x)
例1 求下列微积分
(1)
d n x N n ( x) x n N n 1 ( x) dx d n x N n ( x) x n N n 1 ( x) dx 2n N n 1 ( x) N n 1 ( x) N n ( x) x ' N n 1 ( x) N n 1 ( x) 2 N n ( x)
xJ n ( x) nJ n ( x) xJ n 1 ( x) xJ n ( x) nJ n ( x) xJ n 1 ( x) J n1 ( x) J n1 ( x) 2 J n ( x)
J n 1 ( x) J n 1 ( x)
2n J n ( x) x
n 1


J n ( x)
(1) x m 0 m! ( n m 1) 2
m
n2m
J ( x) cos J ( x) N n ( x) lim n sin
性质4 初值

数学物理方法2课件:第十三章-柱函数-渐近表示式

数学物理方法2课件:第十三章-柱函数-渐近表示式

ν
Jν ( x) + Jν′ ( x) = − Jν +1 ( x)
π
2
由此可以推出: 2 由前面的讨论,可以知道:当ν =1/ 2时,有: J1/2 ( x) = 这样可以得到: 2 这样,再由前面得到的递推关系,有: A0 = A1/2 = 所以最后得到 所以最后得到: 2 ⎛ νπ π ⎞ Jν ( x) ≈ sin ⎜ x − − ⎟ πx 2 4⎠ ⎝ 2 A1/2 = 2 2 2 π⎞ ⎛ sin x = cos ⎜ x − ⎟ πx πx 2⎠ ⎝ , θ1/2 = − Aν =A0 , θν = θ 0 −
d 2 y 1 dy d ⎛ ν2 ⎞ + + ⎜1 − 2 ⎟ y = 0 2 dx x dx ⎝ x ⎠
⎛ 1 − 4ν 2 ⎞ x → ∞ g ′′( x) + g ( x) = 0 = g ′′( x) + ⎜1 + g ( x ) 0 ⎟ x2 ⎠ ⎝
y ( x) ∼ 1 cos ( x + θν ) x
其中 θν 为相角,是待定的常数,但与 ν 有关。
Jν ( x) ∼
1 cos ( x + θν ) x
x 其中Aν 及θν 是两个待定的常数。将上式两边再对x求导,可以得到 x A ≈ − ν sin ( x + θν ) x 另一方面,由递推关系 Jν′ ( x) ≈ − Aν sin ( x + θν ) − 1 Aν cos ( x + θν ) 2x x
3. 贝塞尔函数的渐近表示式 主要是讨论当变量 x→ ∞ 时,柱函数的的渐近行为。先看一下贝塞尔 函数的渐近行为。令
y ( x) = g ( x) / x y ′ = g ′ / x1/ 2 − g / ( 2 x 3/2 ) y ′′ = g ′′ / x1/ 2 − g ′ / x 3/ 2 + 3 g / ( 4 x 5/2 )

贝塞尔函数

贝塞尔函数

贝塞尔函数
贝塞尔函数(Bessel functions),是数学上的一类特殊函数的总称。

通常单说的贝塞尔函数指第一类贝塞尔函数(Bessel function of the first kind)。

一般贝塞尔函数是下列常微分方程(一般称为贝塞尔方程)的标准解函数:
这类方程的解是无法用初等函数系统地表示。

由于贝塞尔微分方程是二阶常微分方程,因此需要两个独立的函数来表示其标准解函数。

通常,第一种贝塞尔函数和第二种贝塞尔函数用于表示标准解函数:
注意,由于在x=0 时候是发散的(无穷),当取x=0 时,相关系数必须为0时,才能获得有物理意义的结果。

贝塞尔函数的特定形式随上述方程式中的任何实数或复数α发生变化(相应地,α称为相应贝塞尔函数的阶数)。

实际应用中最常见的情况是α为整数n,相应的解称为n阶贝塞尔函数。

尽管在上述微分方程中,α本身的正负号不改变方程的形式,但实际应用中仍习惯针对α和−α定义两种不同的贝塞尔函数(这样做能带来好处,比如消除了函数在α=0 点的不光滑性)。

贝塞尔函数也被称为柱谐函数、圆柱函数或圆柱谐波,因为他们是于拉普拉斯方程在圆柱坐标上的求解过程中被发现的。

贝塞尔函数的几个正整数阶特例早在18世纪中叶就由瑞士数学家丹尼尔·伯努利在研究悬链振动时提出了,当时引起了数学界的兴趣。

丹尼尔的叔叔雅各布·伯努利,欧拉、拉格朗日等数学大师对贝塞尔函数的研究作出过重要贡献。

1817年,德国数学家贝塞尔在研究开普勒提出的三体引力系统的运动问题时,第一次系统地提出了贝塞尔函数的总体理论框架,后人以他的名字来命名了这种函数。

贝塞尔函数渐进表达式

贝塞尔函数渐进表达式

贝塞尔函数渐进表达式贝塞尔函数 (Bessel Function) 是一类重要的特殊函数,它在许多领域都有广泛的应用,如物理学、工程学、数学等等。

贝塞尔函数的渐进表达式在各种文献和教材中都有不同的表述和推导,其中最常见的是第一类和第二类柱函数的渐近表达式。

第一类柱函数,即贝塞尔函数 (Bessel Function)Jnu(x),在 x 趋近于正无穷时,其渐进表达式为:$$J_{u}(x)simfrac{1}{x}expleft(-x^{2}ight)left(frac{2x}{pi}ight)^{frac{u}{2}}$$其中 $u$ 是贝塞尔函数的阶数,$sim$ 表示当 $xightarrowinfty$ 时的极限。

第二类柱函数,即诺伊曼函数 (Neumann Function)Nnu(x),在 x 趋近于正无穷时,其渐进表达式为:$$N_{u}(x)simfrac{2}{pi x}expleft(-x^{2}ight)sum_{k=0}^{infty}frac{left(frac{x}{2}ight)^{2k}}{Gamma(u k+1)}$$其中 $u$ 是贝塞尔函数的阶数,$sim$ 表示当 $xightarrowinfty$ 时的极限。

第三类柱函数,即汉克尔函数 (Hankel Function)H${}_{u}$(x),在 x 趋近于正无穷时,其渐进表达式为:$$H_{u}(x)simfrac{1}{x}expleft(-x^{2}ight)left(frac{2x}{pi}ight)^{frac{u+1}{2}}$$其中 $u$ 是汉克尔函数的阶数,$sim$ 表示当 $xightarrowinfty$ 时的极限。

需要注意的是,这些渐近表达式仅仅是贝塞尔函数在一些特定的条件下的渐进表达式,并不是贝塞尔函数的所有渐进表达式。

此外,贝塞尔函数还有很多其他的定义和性质,都可以通过不同的数学工具和方法进行研究和探究。

第十六章 贝塞尔函数

第十六章 贝塞尔函数

(2) 当 不是整数时 ,得到一个特解为
x y1 (1) 2 m 2 k ! ( k 1) k 0
k
2 k
称之为 n 阶第一类贝塞尔函数,记作
1 x 2k J ( x ) ( 1) ( ) k ! ( k 1) 2 k 0
k

(3)当 v =正整数n或零时,有 (n m 1) (n m )!
2n (1) J n1 ( x ) J n ( x ) J n 1 ( x ) ; x 1 ' ( 2) J n ( x ) [ J n1 ( x ) J n1 ( x )] ; 2 ( 3) ( 4) x J n ( x ) nJ n ( x ) xJn1 ( x ) ; x J n ( x ) xJn1 ( x ) nJ n ( x ) ;
(16.11)
(16.13)
令x kr , y( x ) R( r ),则式(16.11)变为 d y dy 2 2 x x ( x ) y 0 2 dx dx
2 2
式(16.13)称为v阶Bessel微分方程
二、Bessel方程的解-Bessel函数
d y dy 2 2 x x ( x ) y 0 2 dx dx
(k 1)! (k 1)k!
(3 1)! 4! (3 1)3!
2( k 1) x 2 k 1 ( 1) 2 k 2 2 ( k 1)!( k 1)!
k 2 k 1 2 ( k 1 ) x ( 1)k 2 k 2 2 k! ( k 1)(k 1)!
k
n=1 ;m=0→∞ :
x x3 x5 x7 x 2 k 1 k J1 ( x ) 3 5 7 (1) 2k 1 2 2 2! 2 2! 3! 2 3! 4! 2 k!( k 1)!

贝塞尔函数 柱函数

贝塞尔函数 柱函数

成柱面问题.(由于是二维问题,即退化为极坐标)
设 u(x, y, t) = u(r,j, t) = T (t)U (r,j ) ,对泛定方程分离变量(取 l = k2 )得
T ¢¢ + k 2a2T = 0 (14.1.2)
ì ï í
U r¢¢
+
1 r
U r¢
+
1 r2
Uj¢¢
+
k 2U
=
0
ïîU |r=l = 0
ìíîHHnn((12))
( x) ( x)
= =
Jn Jn
(x) (x)
+ -
i i
× ×
Nn Nn
(x) (x)
14.3 贝塞尔函数的基本性质
14.3.1 贝塞尔函数的递推公式
由贝塞尔函数的级数表达式(14.2.1)容易推出
d [ Jn ( x) ] = - Jv+1 (x)
dx xn
xv
d dx
+ 0
u
yy
)
ï u(x, y,t) |t=0 = j (x, y)
ïî ut (x, y, t) |t=0 =y (x, y)
(0 £ x2 + y2 < l2 ,t > 0) (t ³ 0)
(14.1.1)
其中 l 为已知正数,j(x, y),y ( x, y) 为已知函数.这个定解问题因宜于使用柱坐标,从而构
14.1 贝塞尔方程及其解
14.1.1 贝塞尔方程
拉普拉斯方程在柱坐标系下的分离变量得出了一般的贝塞尔方程,由于贝塞尔方程的 普遍性,我们还能从其它典型的数学物理定解问题来导出贝塞尔方程的一般形式.

贝塞尔函数 柱函数

贝塞尔函数 柱函数

,令 m 1

d dx
[ x m1J m1 ( x)]
xm1Jm (x)
两边积分,故得到
xm1Jm1(x)
x 0
x m1J m
( x)dx
其中取 m 1 ,即为(22.3.18)式。
20.3.2 贝塞尔函数与本征问题
拉普拉斯方程在柱坐标系下的分离变量,得到了方程(14.6.7) 即
d2R
( x ) 2k 1) 2
式中 (x) 是伽马函数.满足关系
(20.2.1)
( k 1) ( k)( k 1)L ( 2)( 1)( 1)
当 为正整数或零时, ( k 1) ( k)!
当 取整数时 ( k 1) ,(k 0,1, 2,, 1)
所以当 n 整数时,上述的级数实际上是从 k n
将 (20.3.7)与(20.3.8)相加或相减消去 Z 1 或 Z 1 分别得到
Z
(
x)
x
Z
(
x)
Z
1
(
x)
(20.3.9)
Z
1
(
x)
x
Z
(
x)
Z
(
x)
(20.3.10)
将(20.3.9) 式中的 换成 1,得到
Z
(x)
Z1 ( x)
1 x Z
1 ( x)
(20.3.11)
将 (20.3.10)代入上式,立即得到 Z (x) 满足 阶贝塞尔方程.
Jm (k 0 ) 0 则
J m ( x) 0 就是决定本征值的方程.
若用
x(m) n
表征 Jm (x) 0
的第
n
个正根,于是本征值
(m) n

11.2 贝塞尔方程

11.2 贝塞尔方程
x
2
d R dx
2
2
x
dR dx
x m
2
2
R 0
(11.4.1)
数学物理方法
令 i x , y ( ) R ( x ) 代入上式,则得到贝塞尔方程
y y m
2 2 2
y0
(11.4.2)
令 i x , 即可得到虚宗量贝塞尔方程的解。 定义虚宗量贝塞尔方程的解具有如下形式
解:采用柱坐标系,极点在下底中心, z 轴沿圆柱的轴, 定解问题表为
2u 0 u 0, 0 u z 0 f 1 ( ), u
zL
f2 ( )
本例是圆柱内部的拉普拉斯方程定解问题, 柱侧是齐次的第 二类边界条件,故考虑 0 的情况。
况应舍弃。 故把特解叠加起来,有
v


Ap I0 (
p L
) sin
p z L
p 1
为确定系数,将上式代入柱侧的边界条件 q0 p p p z I '0 ( 0 ) sin Ap
p 1
L
L
L
k
数学物理方法
例 2 半径 0 ,高 L 的导体圆柱壳,用不导电的介质将柱壳 的上下底面和侧面隔离开,柱壳侧面电势为 u 0 z / L ,上底 面电势为 u 1 ,下底面接地,求柱壳外电势分布
v [ A J 0 ( ) B N 0 ( )]e
a t
2
A J 0 ( 1 ) B J 0 ( 1 ) 0 代入边界条件, ,从而解 AJ 0 ( 2 ) BJ 0 ( 2 ) 0
出本征值 ,从而定出相应系数,得解。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

6-2 贝塞尔函数柱函数在用分离变量法一章介绍了拉普拉斯方程在柱坐标系下分离变量得到了一种特殊类型的常微分方程:贝塞尔方程.通过幂级数解法得到了另一类特殊函数,称为贝塞尔函数.贝塞尔函数具有一系列性质,在求解数学物理问题时主要是引用贝塞尔函数的正交完备性.6.1 贝塞尔方程及其解6.1.1 贝塞尔方程拉普拉斯方程在柱坐标系下的分离变量得出了一般的贝塞尔方程。

考虑固定边界的圆膜振动,可以归结为下述定解问题222222200() (0,0)|0 (0)(,,)|(,)(,,)|(,)tt xx yy x y l t tt u a u u x y l t u t u x y t x y u x y t x y ϕψ+===⎧=+≤+<>⎪=≥⎪⎨=⎪⎪=⎩(6.1.1)其中l 为已知正数,(,),(,)x y x y ϕψ为已知函数.这个定解问题宜于使用柱坐标,从而构成柱面问题.(由于是二维问题,即退化为极坐标)设(,,)(,,)()(,)u x y t u t T t U ρϕρϕ==)得220a T =(6.1.2)22100 U U k U ρϕρ′′′++=(6.1.3)再令(,)()()U R ρϕρϕ=Φ,得到2ν′′Φ+Φ=(6.1.4)2222()0R R k R ρρρν′′++−=(6.1.5)于是(6.1.5)得到22d ()0d y x x y xν+−=(6.1.6)边界条件为()|()0l y k y kl ρρ===方程(6.1.6)称为ν阶贝塞尔微分方程.这里νx和可以为任意数.6.1.2贝塞尔方程的解通过数学物理方程的幂级数求解方法可以得出结论:(1)当ν≠整数时,贝塞尔方程(6.1.6)的通解为()J ()J ()y x A x B x νν−=+(6.1.7)其中,A B 为任意常数,J ()x ν定义为ν阶第一类贝塞尔函数但是当n ν=整数时,有J ()(1)J ()nn n x x −=−故上述解中的J ()n x 与J ()n x −是线性相关的,所以(6.1.7)成为通解必须是ν≠整数.(2)当ν取任意值时:定义第二类贝塞尔函数N ()x ν,这样贝塞尔方程的通解可表示为()J ()N ()y x A x B x νν=+(6.1.8)(3) 当ν取任意值时:由第一、二类贝塞尔函数还可以构成线性独立的第三类贝塞尔函数H ()x ν,又称为汉克尔函数.(1)(2)H ()J ()iN ()H ()J ()iN ()x x x x x x νννννν⎧=+⎨=−⎩(6.1.9)分别将(1)(2)H ,H νν称为第一种和第二种汉克尔函数.于是贝塞尔方程的通解又可以表示为(1)(2)(H ()H ()y x A x B x νν=+(6.1.10)最后,总结ν阶贝塞尔方程的通解通常有下列三种形式:(i )()J ()J () (y x A x B x ννν−=+≠整数)(ii )()J ()N ()(y x A x B x ννν=+可以取任意数)(iii )(1)(2)()H ()H ()(y x A x B x ννν=+可以取任意数)6.2 三类贝塞尔函数的表示式及性质6.2.1 第一类贝塞尔函数的表示式第一类贝塞尔函数J ()x ν的级数表示式为2201()1)2()!(1)2kkkk x x k k ννν∞+−+=Γ−++∑∑(6.2.1)伽马函数.满足关系()(1)(2)(1)(1)k k ννννν++−++Γ+"当ν为正整数或零时,(1)()!k k ννΓ++=+当ν取整数时(1),(0,1,2,,1)k k ννΓ−++=∞=⋅⋅⋅−所以当n ν=整数时,上述的级数实际上是从k n=的项开始,即, (0)n ≥(6.2.2)22011)()!(1)21(1)(), ()!(1)2kn kln lx k n k x l k n l n l −+∞+=Γ−++−=−Γ++∑(6.2.3)所以J ()(1)J ()nn n x x −=−(6.2.4)同理可证J ()J ()n n x x −=−(6.2.5)因此有重要关系J ()(1)J ()nx x −=−(6.2.6)贝塞尔函数表示式246223511()()()2(2!)2(3!)211()()2!22!3!2x x x x x −+−+−+−""当x 很小时(0)x →,保留级数中前几项,可得1J ()(),(1,2,3,)2(1)x x νννν≈≠−−−⋅⋅⋅Γ+(6.2.7)特别是0J (0)1,J (0)0 (=1,2,3,)n n ==⋅⋅⋅(6.2.8)32π)()π42o x x ν−−+(6.2.9)试证半奇阶贝塞尔函数122J ()s in πx x x=(6.2.1)有而13135(21)()π22k k k +⋅⋅⋅⋅+Γ+="(21)!k +2sin πx x122J ()cos πx xx−=121212220J ()(1)12!(1)2k kk k x x k k +∞+==−Γ++∑6.3 贝塞尔函数的基本性质6.3.1 贝塞尔函数的递推公式由贝塞尔函数的级数表达式(6.2.1)容易推出1J ()J ()d []v x x ν+=−(6.3.1))(6.3.2)以上两式都是贝塞尔函数的线性关系式.诺伊曼函数N()v x 也应该满足上述递推关系若用()v Z x 代表v阶的第一或第二或第三类函数,总是有1d [()]()d v vv v x Z x x Z x x−=(6.3.3)d [()]()v vx Z x x Z x −−=−(6.3.4)1)()()v v vx Z x Z x x+−=−(6.3.5)1()()v v vZ x Z x x−=(6.3.6)从(6.3.5)和(6.3.6)消去Z ν或消去Z ν′可得11()()2()v v v Z x Z x Z x +−′=−112()()()v v v vZ x Z x Z x x+−=−+即为从)(x Z 和)(x Z 推算)(1x Z v +的递推公式.1()2()v v vZ x Z x x++=(6.3.7)1)()2()v Z x Z x ν+′−=(6.3.8)任一满足一组递推关系的函数)(x Z v 统称为柱函数例6.3.1证明柱函数满足贝塞尔方程【证明】以满足(6.3.7)和(6.3.8)这一组递推公式来进行证明:将(6.3.7)与(6.3.8)相加或相减消去1Z ν+或1Z ν−分别得到()()()Z x Z x Z x ν′+=(6.3.9)x(6.3.10)换成1ν+,得到ν111()()x Z x x νν+++′+(6.3.11)将(6.3.10)代入上式,立即得到()Z x ν满足ν阶贝塞尔方程.例6.3.2 求2J ()d x x x∫【解】根据公式(6.3.8)11()()2()v v Z x Z x Z x ν−+′−=有201J ()J ()2J ()x x x ′=−20111111010J ()d J ()d 2J ()d J ()2[J ()J ()d ]J ()2[J ()J ()d ]J ()2J ()x x x x x x x x x x x x x x x x x x x x x x x x c′=−=−−′=−+=−−+∫∫∫∫∫例6.3.3 证明下式成立1110J ()d J ()xm m m m x x x x x +++=∫(6.3.17)特别是22J ()d J ()xx x x x x =∫(6.3.18)】利用递推公式(6.3.2)即1J ()vv x x −=,令1m ν=+则两边积分,故得到111d [J ()]J ()d m m m m x x x x x+++=1J ()d m x x∫1,即为(22.3.18)式。

6.3.2 贝塞尔函数与本征问题拉普拉斯方程在柱坐标系下的分离变量,得到了方程(14.6.7) 即2222d 1d ()0d d R R R νμρρρρ++−=(6.3.19)取整数,其它情况下ν分三种情况:0=μ,0>μ和0<μ进行讨论.方程(6.3.19)是欧拉方程;作代换x μρ=,则得到()22222d d 0 ()d d R R x x x R x x xνμρ++−==(6.3.21)即为ν阶贝塞尔(Bessel)方程.(3)0<μ.记20k μ−=>,以2kμ=−()22d 0d R x R xν−+=(6.3.22)虚宗量贝塞尔方程.如把贝塞尔方程(6.3.22)的宗量i x,就成了方程(6.3.21)贝塞尔方程本征值问题(即本征值0μ>的情况):1. 第一类边界条件的贝塞尔方程本征值问题220201d d []()()0(0)d d ()0 |(0)|R k R R R M νρρρρρρρρρ⎧+−=≤≤⎪⎨⎪=<⎩(6.3.23)()(2π+)ϕϕ=Φ0,1,2,3,="可进一步化为施—刘型本征值问题的形式20)()0(0)(0)|R k R R M ρρρρ+=≤≤<(6.3.24)相应于施-刘型方程中的22(), (),(),mk x x q x x x kxρλμ==−===故施-刘型本征值问题的结论对于贝塞尔方程的本征值问题也(6.3.24)的通解为)N ()m B μρμρ+(6.3.25)表征J ()0m x =的第n个正根,于是本征值())220][](1,2,3,)m m nxn ρ=="(6.3.26)代入边界条件决定本征值及本征函数.因为(0)R M<故B =又0()0R ρ=,要A ≠,则必须则=决定本征值的方程.施-刘型本征值问题的结论(1) 本征值存在,且都是非负的实数;<"(2) 本征值可编成单调递增的序列())()2222)()()m m nxxρρ<<<<""(6.3.27)()()()12J (),J (),,J (),m m m nm m m xxxρρρρρρ""(6.3.28)(3)对于每一个本征值有一个相应的本征函数()()2[]m m nnkλ=()J()m mnk ρ且本征函数()J ()m mnkρ在0[0,]ρ区间上有(1)n −个零点即()()()112000()()(),,,m m m n m m m nnnx x x xxxρρρ−"贝塞尔函数有无穷个零点.零点与的零点是彼此相间分布的,即1J ()m x +的任意两个相邻零点之间必有且仅有一个1J ()m x +1J (m x +(5)以()m nx表示J ()m x 的第n 个正零点,则()()1lim []πm m n nn xx+→+∞−=,即J ()m x 几乎是以2π为周期的周期函数.24612(1)3(4)15(4)105(4)C D EA A A ++++"其中221π(2),4,731,83982377922A m nB mC BD B B =−+==−=−+32694915385515857436277237E B B B =−+−)0=00)]J ()0m k k ρρρ=′==(6.3.30)对于不过,0k ≠,则本征值()()20(/)m m nnxλρ=(6.3.31) 其中()m nx是J ()m x ′的第n个零点.J ()m x ′的零点在一般的数学用表中并未列出.0m =的特例还是容易得到的:由公式(6.3.12)得到01J ()J ()x x ′=−这样,至于0J ()x ′的零点不过就是1J ()x 的零点,可从许多数学用表中查出.m ≠的情况,J ()m x ′的零点()m nx11[J ()J ()]m m x x −+−的零点可从曲线1J ()m x −和1J ()m x +对于0m ≠J ()mx ′的情况,的零点()m nx还可以用下面的公式计算:()m nx=35386(4)15(4)B C D A A A A +−−−−"2222),4,78292207530393537n B m C B B B B π==+−−+这个条件就是(3) 第三类齐次边界条件.0)()(00=′+ρρR H R ()()()00J ()J ()0m m m m nn m nkHk kρρ′+=记()000, /,m nx kh H ρρ==并引用(6.3.5)可将上式改写为,其中010J ()J ()m m x x x h m+=+所以本征值()()20(/)m m nnxμρ=()m nx贝塞尔函数正交性和模对应不同本征值的本征函数分别满足(6.3.35)乘以2()2()2dJ d []{[]}J ()0d d m m m i m i m k k ρρρρρ+−=(6.3.34)2()2()2dJ d []{[]}J ()0d d m m m j m j m k k ρρρρρ+−=(6.3.33))()m j ρ(J m i k ρ()J (m m ikρ两式相减,再积分,利用分部积分法得到0()()0()()()0J ()J ()d d d J ()J ()J ()]|0d d m m m i m jm m m m j m j m i k k k k k ρρρρρρρρρρρρ−=∫故当2.贝塞尔函数的模()()m m ij kk≠0()()0J ()J ()d 0m m m im jkk ρρρρρ=∫(6.3.36)()m nN为了用贝塞尔函数作基进行广义傅立叶级数展开,需要先()J ()m m nk ρ计算贝塞尔函数的模()m nN()2()2[][J ()]d m m nm nNkρρρρ=∫(6.3.37)注意对于()()()22[][]m m m nnnxkλρ==()0m n x >()0m x>把()m nkρ记为x2221[J ()]d [J ()]d()2x x m m nx x x x x λ=∫∫222001[J ()][J ()]J ()d x x m m m n x x x x x xλ′−∫(6.3.38)()222220011[][J ()][J ()J ()J ()]J ()d2x x m n m m m m m n n N x x x x x x m x x x λλ′′′=++−∫022222000dJ ()11[J ()][J ()(J )]d J dJ 2d x x x mm m mm m n n n x m x x x x x x x λλλ′′′=++−∫∫0220)][J ()]2xm nm x x λ−02222001)J ()][J ()]2x x m m nm x x x λ′+()22()20001)[J ()][J ()].2m m m n m n k k ρρρ′+第一类齐次边界条件则式(6.3.38)成为()0J()0,m mnkρ=()22()201[][J ()]m m Nk ρρ′=(6.3.39))代入上式,并且考虑到第一类齐次边界条件0)0,ρ=故得2()20101[J ()]2m m n k ρρ+(6.3.40)第二类齐次边界条件(6.3.38)成为()0J ()0,m m nkρ′=2()22()2001[]()[J ()]2m m nm n m Nk ρρλ=−(6.3.41))m H(6.3.38)成为222()2000)[J ()]m m n n n m k H ρρρλλ−+(6.3.42)6.3.3 广义傅立叶-贝塞尔级数按照施-刘型本征值问题的性质,本征函数族()J ()m m nkρ是完备的,可作为广义傅立叶级数展开的基.)(ρ可以展开为广义的傅立叶-贝塞尔级数()1()J ()m n m nn f f kρρ∞==∑(6.3.43)广义傅氏系数()20()J ()d ]m m nf kρρρρρ∫(6.3.44)例6.3.4在区间0[0,]ρ上,以(0)0J ()nk ρ为基,把函数0()f uρ=(常数)展开为傅里叶-贝塞尔级数.说明:其中(0)(0)2[]nnk λ=是本征函数(0)0J ()nkρ对应的本征值.【解】根据(6.3.43)和(6.3.44)则(0)001J ()n nn u f kρ∞==∑其中系数这里的0(0)00(0)21J ()d []n nnf u kNρρρρ=∫(0)nN由第一类边界条件所对应的模公式(6.3.40)给出.(0)22[]nxρ=是0阶贝塞尔函数0J ()x 的第n 个零点,贝塞尔函数表查出.这样令(0)2(0)200102J()d[J()]nnnu xfxρρρρρρ=∫(0)nxxρρ=,则(0)10(0)2(0)(0)12[J()]()]J()nxn n nu ux xx x x=(0)(0)(0)1102J()J()nn n nu xx xρρ∞=∑6.3.4 贝塞尔函数的母函数(生成函数)1. 母函数(生成函数)考虑解析函数)1(2),(zzxezxG−=在此处的x为参变数,不是复变数z的实部.∑∞==2!)2(k kkzx zk x e∞−−−)2(1l l z x x ∑∑∞=∞=−−=00)1)(!)2(!)2(k l l lk k zz l x z k x对于固定的z +∞<<z 0,以上两级数在内是可以相乘的,且可按任意方式并项.,,2,1,0,"±±==−n n l k 2000(1)(1)()[()]!!2()!!2l lk l k l l n n l n l x x z z k l n l l ∞∞∞∞+−+===−∞=−−=+∑∑∑∑称平面波用柱面波的形式展开(,)J()nnn G x z x z∞=−∞=∑(6.3.45))1(2z z x e−为贝塞尔函数的母函数(或生成函数).x kr=)的傅氏余弦展开式.i cos i i i 01J()i J ()[J ()i J ()i ]kr n n n n n n nn n n n ekr ekr kr ekr eθθθθ∞∞−−−=−∞===++∑∑1()cos n n kr n θ=∑(6.3.46)i cos kr eθ。

相关文档
最新文档