贝塞尔函数

贝塞尔函数

1. 贝塞尔方程及解:

令u (, ϕ, τ)=R(, )Φ(ϕ)Z(τ)为分离变量的解,则R (, )满足本征值问题的方程,

∂2R 1dy dR ⎛2m 2⎫ 2++ ω-2⎪R =0 (17.1.1) ∂dx d ⎝⎭

其中ω2是分量的本征值问题的本征值。 R() =R() =y(x);m=ν 则上面方程可以变换:若作变换x =ω(或x =λ) ;x

ω

x 2y //+x 2y /+(x-ν2) y =0 (17.1.1a )

当ν≠整数时,贝塞尔方程的通解为:

y (x)=AJ ν(x)+BJ -ν(x)

当ν=整数时,由于J -m =(-1) m J m (x),因此通解为

y (x)=AJ m (x)+BY m (x)

式中A 与B 为任意常数,J m (x)与Y m (x)分别定义为 m 阶第一类与m 阶第二类贝塞尔函数。

2. 贝塞尔方程的的级数解

二阶线性齐次常微分方程x 2y '' +xy ' +(x2-υ2) y =0,0≤x ≤b 为贝塞尔方程

现在x=0的领域求解贝塞尔方程的解

2.1级数解的形式

由1p(x)=x ν2,q(x)=1-2可见,x=0是p=(x )的一阶极点,是q(x)x

的二阶极点。因此,x=0是方程的正则奇点,方程的第一解具有形式;

y =x ∑C k x =∑C k x k +p 2.1.1 k =0k =0n ∞k ∞

2.2指标方程

将2.1.1代入贝塞尔方程可得:

22k +p k +ρ+3⎤∑⎡(k+ρ) -νC x +∑C x =0 2.1.2 k k ⎣⎦k =0k =0∞∞

由x 的最低次幂x ρ的系数为0,即得:

x ρ(ρ2-ν2)C 0=0

因C 0≠0,即得指标方程ρ2-ν2=0。由此得指标

ρ1=ν, ρ2=-ν

2.3.系数递推公式

为确定起见,令ν>0,并将ρ=ρ1=ν代入2.1.2中得到

22k +νk +ν+2⎤∑⎡(k+ν) -νC x +∑C x =0 k k ⎣⎦k =0k =0∞∞

改变第二项的求和指标,可得

k =0∑k(k+2ν) C k x k +ν+∑C k -2x k +ν=0 k =2∞∞

由x 的同次幂数之和为0,(1+2ν) C 1=0 k(k+2ν) k C k +C k -2=0

由此得

C 1=0

2.4. 推公式求系数得特解

………

将系数代入1.1中的贝塞尔方程的一个特解为 (-1) C k =C k -2 k(k+2ν)

(-1) n Γ(-ν+1)C 02n +νy 1(x)=∑2n x n =02n ! Γ(ν+n +1) ∞

2.5. 另一个特解

同理,令ρ=ρ2=-ν可得另一个特解为

(-1) n Γ(-ν+1)C 02n -νy 2(x)=∑2n x n =02n ! Γ(-ν+n +1) ∞

3. 第一类贝塞尔函数

第一类贝塞尔函数J ν(x)的级数形式为

dy x ν+2k J ν(x)=∑(-1) () k ! Γ(ν+κ+1) 2k =1∞k

m 经过证明可得:J -m (x)=(-1) J , m (x)

同理可得:J -m (x)=J , m (x)

m J (x)=(-1) J -m (x) 因此:, m

4. 第二类贝塞尔函数:

第二类贝塞尔函数是Weber 和Schlafli ,通常把它定义为

Y ν(x)

Y m (x)的级数形式为

Y m (x)=cos νπJ (x)-J -ν(x) sin νπ

2⎡x ⎤1m -1(m-k -1)! x -m +2k 1∞(-1) k x m +2κγ+ln J (x)-() -ϕ(k)+ϕ(m+κ) () {}∑∑m k ⎥π⎢2πk ! 2πk !(m+k) 2⎣⎦k =0k =0

式中γ=0.577216,而

ϕ(k)=∑1

n =1κn

当x 很小时,可得

Y 0≈

当x 很大时,

Y ν(x)≈

5. 第三类贝塞尔函数

通常定义为

(1) H ν=J ν(x)+iY ν(x)

(2) H ν=J ν(x)-iY ν(x) 2πlnx (ν=0) πνx (x--) (17.1.12) 42

则方程(17.1.1 a)的通解可以写成为

(1)(2)+B H ν(x) y(x)=AH ν

当x →∞时其渐进展开式为

H ν(1)3-i (x-ν2x ) =+o (x2) (17.1.14a )

3) --i (x-ν

2x -π

4=+o (x2) (17.1.14b ) H ν(2)

当x →0时其渐进展开式为

H ν(x)≈-i (ν-1)! 2() πx

2(2)H ν(x)≈-i ln x (ν>0) π

总结上述,ν阶贝塞尔方程

x 2y +xy /+(x2-ν2) y =0 的通解有三种形式:

(1)y(x)=AJ(x)+BJ (x) (ν≠0)

(2)y(x)=AJ(x)+BY ν(x) (ν可取任意整数)

(1) (3)y(x)=AH ν(x)+BH ν(2)(x) (ν可取任意整数)

其中A,B 为常数。

6贝塞尔函数的基本性质

6.1生成函数

如果一个函数的级数展开式系数是贝塞尔函数,则该函数为贝赛尔函数的生成函数,即如果有

n J (x)r f(x,r)=∑n n

则称f(x,r)为J n (x)的生成函数,其中r 为参数。生成函数有非常广泛的应用。

整数阶贝塞尔函数J n (x)的生成函数为

∞⎡x ⎛1⎫⎤n exp r -=J x r ⎪⎥∑n () ⎢2 r ⎭⎦n =-∞⎣⎝

6.2递推公式

阶数邻近的贝塞尔函数之间存在一定的关系式,称为递推公式。基本的递推公式是

d νν⎡⎤x J X =x J ν-1(X ) ()ν ⎣⎦dx

d -ν-ν⎡⎤x J X =-x J ν+1(X ) ()ν⎣⎦dx

贝塞尔函数还有如下两个常用的递推公式

1J ν(x)=⎡J ν-1(x )-J ν+1(x )⎤⎣⎦ 2'

2νJ ν-1(x )+J ν+1(x )=J ν(x ) x

N 阶贝塞尔函数的积分表示为:

当n=0时 J n (x )=π⎰1ππ⎰1π0cos (x sin θ-n θ)d θ

0c o (s x s θi )n θd

这是零阶贝塞尔函数积分表示

可以看出,贝塞尔函数的积分表示其实是用正弦和余弦函数积分形式表示贝塞尔函数。

6.3渐近公式

余弦函数是J 0(x )的渐近公式

J 0(

x )≈π-) 4

一般情况下,你n 阶贝塞尔函数J 0(x )有渐近公式

J n (

x )≈

πn π--) (n=0,1,2…. ) 42

贝塞尔函数

1. 贝塞尔方程及解:

令u (, ϕ, τ)=R(, )Φ(ϕ)Z(τ)为分离变量的解,则R (, )满足本征值问题的方程,

∂2R 1dy dR ⎛2m 2⎫ 2++ ω-2⎪R =0 (17.1.1) ∂dx d ⎝⎭

其中ω2是分量的本征值问题的本征值。 R() =R() =y(x);m=ν 则上面方程可以变换:若作变换x =ω(或x =λ) ;x

ω

x 2y //+x 2y /+(x-ν2) y =0 (17.1.1a )

当ν≠整数时,贝塞尔方程的通解为:

y (x)=AJ ν(x)+BJ -ν(x)

当ν=整数时,由于J -m =(-1) m J m (x),因此通解为

y (x)=AJ m (x)+BY m (x)

式中A 与B 为任意常数,J m (x)与Y m (x)分别定义为 m 阶第一类与m 阶第二类贝塞尔函数。

2. 贝塞尔方程的的级数解

二阶线性齐次常微分方程x 2y '' +xy ' +(x2-υ2) y =0,0≤x ≤b 为贝塞尔方程

现在x=0的领域求解贝塞尔方程的解

2.1级数解的形式

由1p(x)=x ν2,q(x)=1-2可见,x=0是p=(x )的一阶极点,是q(x)x

的二阶极点。因此,x=0是方程的正则奇点,方程的第一解具有形式;

y =x ∑C k x =∑C k x k +p 2.1.1 k =0k =0n ∞k ∞

2.2指标方程

将2.1.1代入贝塞尔方程可得:

22k +p k +ρ+3⎤∑⎡(k+ρ) -νC x +∑C x =0 2.1.2 k k ⎣⎦k =0k =0∞∞

由x 的最低次幂x ρ的系数为0,即得:

x ρ(ρ2-ν2)C 0=0

因C 0≠0,即得指标方程ρ2-ν2=0。由此得指标

ρ1=ν, ρ2=-ν

2.3.系数递推公式

为确定起见,令ν>0,并将ρ=ρ1=ν代入2.1.2中得到

22k +νk +ν+2⎤∑⎡(k+ν) -νC x +∑C x =0 k k ⎣⎦k =0k =0∞∞

改变第二项的求和指标,可得

k =0∑k(k+2ν) C k x k +ν+∑C k -2x k +ν=0 k =2∞∞

由x 的同次幂数之和为0,(1+2ν) C 1=0 k(k+2ν) k C k +C k -2=0

由此得

C 1=0

2.4. 推公式求系数得特解

………

将系数代入1.1中的贝塞尔方程的一个特解为 (-1) C k =C k -2 k(k+2ν)

(-1) n Γ(-ν+1)C 02n +νy 1(x)=∑2n x n =02n ! Γ(ν+n +1) ∞

2.5. 另一个特解

同理,令ρ=ρ2=-ν可得另一个特解为

(-1) n Γ(-ν+1)C 02n -νy 2(x)=∑2n x n =02n ! Γ(-ν+n +1) ∞

3. 第一类贝塞尔函数

第一类贝塞尔函数J ν(x)的级数形式为

dy x ν+2k J ν(x)=∑(-1) () k ! Γ(ν+κ+1) 2k =1∞k

m 经过证明可得:J -m (x)=(-1) J , m (x)

同理可得:J -m (x)=J , m (x)

m J (x)=(-1) J -m (x) 因此:, m

4. 第二类贝塞尔函数:

第二类贝塞尔函数是Weber 和Schlafli ,通常把它定义为

Y ν(x)

Y m (x)的级数形式为

Y m (x)=cos νπJ (x)-J -ν(x) sin νπ

2⎡x ⎤1m -1(m-k -1)! x -m +2k 1∞(-1) k x m +2κγ+ln J (x)-() -ϕ(k)+ϕ(m+κ) () {}∑∑m k ⎥π⎢2πk ! 2πk !(m+k) 2⎣⎦k =0k =0

式中γ=0.577216,而

ϕ(k)=∑1

n =1κn

当x 很小时,可得

Y 0≈

当x 很大时,

Y ν(x)≈

5. 第三类贝塞尔函数

通常定义为

(1) H ν=J ν(x)+iY ν(x)

(2) H ν=J ν(x)-iY ν(x) 2πlnx (ν=0) πνx (x--) (17.1.12) 42

则方程(17.1.1 a)的通解可以写成为

(1)(2)+B H ν(x) y(x)=AH ν

当x →∞时其渐进展开式为

H ν(1)3-i (x-ν2x ) =+o (x2) (17.1.14a )

3) --i (x-ν

2x -π

4=+o (x2) (17.1.14b ) H ν(2)

当x →0时其渐进展开式为

H ν(x)≈-i (ν-1)! 2() πx

2(2)H ν(x)≈-i ln x (ν>0) π

总结上述,ν阶贝塞尔方程

x 2y +xy /+(x2-ν2) y =0 的通解有三种形式:

(1)y(x)=AJ(x)+BJ (x) (ν≠0)

(2)y(x)=AJ(x)+BY ν(x) (ν可取任意整数)

(1) (3)y(x)=AH ν(x)+BH ν(2)(x) (ν可取任意整数)

其中A,B 为常数。

6贝塞尔函数的基本性质

6.1生成函数

如果一个函数的级数展开式系数是贝塞尔函数,则该函数为贝赛尔函数的生成函数,即如果有

n J (x)r f(x,r)=∑n n

则称f(x,r)为J n (x)的生成函数,其中r 为参数。生成函数有非常广泛的应用。

整数阶贝塞尔函数J n (x)的生成函数为

∞⎡x ⎛1⎫⎤n exp r -=J x r ⎪⎥∑n () ⎢2 r ⎭⎦n =-∞⎣⎝

6.2递推公式

阶数邻近的贝塞尔函数之间存在一定的关系式,称为递推公式。基本的递推公式是

d νν⎡⎤x J X =x J ν-1(X ) ()ν ⎣⎦dx

d -ν-ν⎡⎤x J X =-x J ν+1(X ) ()ν⎣⎦dx

贝塞尔函数还有如下两个常用的递推公式

1J ν(x)=⎡J ν-1(x )-J ν+1(x )⎤⎣⎦ 2'

2νJ ν-1(x )+J ν+1(x )=J ν(x ) x

N 阶贝塞尔函数的积分表示为:

当n=0时 J n (x )=π⎰1ππ⎰1π0cos (x sin θ-n θ)d θ

0c o (s x s θi )n θd

这是零阶贝塞尔函数积分表示

可以看出,贝塞尔函数的积分表示其实是用正弦和余弦函数积分形式表示贝塞尔函数。

6.3渐近公式

余弦函数是J 0(x )的渐近公式

J 0(

x )≈π-) 4

一般情况下,你n 阶贝塞尔函数J 0(x )有渐近公式

J n (

x )≈

πn π--) (n=0,1,2…. ) 42


相关文章

  • 贝塞尔函数的有关公式
  • C. 贝塞尔函数的有关公式 贝塞尔方程 的持解B p (z ) 为(柱) 贝塞尔函数.有 第一类柱贝塞尔函数J p (z ) p 为整数n 时,J -n =(-1) n J n : p 不为整数时,J p 与J -p 线性无关. 第二类柱贝 ...查看


  • 高等数学大纲(物理类)
  • <高等数学>教学大纲 课程名称:高等数学 适用层次.专业:理科.工科各专业 学 时:320学时 学 分:20学分 课程类型:通识教育平台课 课 程 性 质:必修课 一.课程的教学目标与任务 高等数学是理.工.管等相关专业的第一基 ...查看


  • Excel线性内插值计算实例
  • Excel 线性内插值计算实例 excel 插值,通常通过函数公式完成.如果手工计算插值,N 多烦琐. excel 插值函数有trend 和forecast 函数.通过这两个excel 插值函数来写相关的插值公式. 下面是两个相关excel ...查看


  • 数学物理方程公式总结
  • ===================== 无限长弦的一般强迫振动定解问题 ⎧u tt =a 2u xx +f (x , t )(x ∈R , t >0) ⎪ ⎨u t =0=ϕ(x ) ⎪ ⎩u t t =0=ψ(x ) 解u (x ...查看


  • 数学物理方程 课程论文
  • 数学物理方程在化工中的应用简介 摘要:数学是科学的基础, 是现代科学和工程技术的核心.虽然数学研究和数学教育的兴旺发达, 并不会自动地带来先进的科学技术和强大的化学工业, 但后者要先进和强大, 则离不开数学的繁荣与发展.数学在化学工程中的应 ...查看


  • Mie理论简介
  • 对流层大作业 --Mie理论的理解与介绍 学院:物理与光电工程学院 班级:071261班 学号:07126006 姓名:彭甜 指导老师:弓树宏 摘要: 本文首先推导了Mie 散射理论的基本公式,从而可以精确的计算出散射光强与各项参数之间的关 ...查看


  • 图形旋转和歪斜的椭圆
  • 图形旋转和歪斜的椭圆 介绍 窗口中的矩形,带圆角的矩形和椭圆只能由GDI 在轴向上绘制.假如有人希望在 Windows NT下绘制旋转或歪斜的图形,他可以使用世界坐标系变换.很不幸的是 在Windows 95/98下,是没有世界坐标系变换的 ...查看


  • 量子力学 第二版 第六章__散射 习题答案 周世勋
  • 第六章 散射 1.粒子受到势能为 U(r) ar 2 的场的散射,求S分波的微分散射截面. [解] 为了应用分波法,求微分散射截面,首先必须找出相角位移.注意到第l个分波的相角位移l是表示在辏力场中的矢径波函数Rl和在没有散射势时的矢径 ...查看


  • 实验二.傅里叶级数与傅里叶变换
  • 实验二.傅里叶级数与傅里叶变换 上机作业: 1. 在(0,T)周期上,锯齿波可表为f(x)=x/3,锯齿波可展开为傅里叶级数: TT63 试画出锯齿波的幅度频谱图. 解: %Fig2d2.m T=1;x=0.2;H=1; t=-0.5* ...查看


热门内容