Chapter 13 柱坐标下的分离变量法 Bessel函数
- 格式:doc
- 大小:2.55 MB
- 文档页数:33
Chapter 13 柱坐标下的分离变量法 Bessel 函数Abstracts :以3+1D 实例通过柱坐标系下方程的变量分离与求解 引入各种柱函数(Bessel 函数、虚宗量Bessel 函数 和球Bessel 函数等)。
在分析这些函数性质的基础上, 表述相应定解问题的物理解。
一、柱坐标下的变量分离1. 柱坐标系下的稳定问题(Laplace 方程)222222110,u u uu zρρρρρϕ⎛⎫∂∂∂∂∇=++= ⎪∂∂∂∂⎝⎭ (1) 即:()22110.zz u u u u ρϕϕρρρρ∇=++= (2)令(,,)()()()u z R Z z ρϕρϕ=Φ,代入(2)得:()20.Z RZ R R Z ρρρΦ''''''+Φ+Φ= (3)2(3)R Zρ⨯Φ得:()2'.R Z RZρρρλ'''''Φ+=-=Φ(4) 分离变量得:()20.(5)'.(6)R Z R Zλρρρλ''Φ+Φ=⎧⎪'⎨''+=⎪⎩(5)与周期性边界条件(0)(2),(0)(2)ππ''Φ=ΦΦ=Φ构成本征值问题。
解得:2 (0,1,2,3,),m m m λ== (){cos ,sin }.m m m ϕϕϕΦ= (6)即为()22'R Z mRZρρρ'''+=⇒分离变量()22'.R m Z RZρμρρ'''-=-=- 得:()2220.0.Z Z R R m R μρρμρ''-=⎧⎪⎨'''++-=⎪⎩ 这两个方程,先求解哪一个以及μ如何取值,取决于哪一个可构成本征值问题,也就是取决于定解问题的边界条件。
如果()R ρ构成本征值问题,则()2220,R R m R ρρμρ'''++-=式中μ的取值范围不同,方程解的形式与性质不同。
题目:深入探讨分离变量法求解拉普拉斯方程和贝塞尔函数的应用在物理和工程领域中,求解微分方程是一项至关重要的任务。
而分离变量法作为一种常见的求解微分方程的方法,在求解拉普拉斯方程和应用贝塞尔函数中发挥着重要作用。
本文将深入探讨分离变量法在求解拉普拉斯方程和应用贝塞尔函数时的原理和应用,以及相关的物理和工程实际问题。
1. 分离变量法的基本原理我们需要了解分离变量法的基本原理。
对于一个多变量的微分方程,如果可以将变量分离,化为单变量的微分方程,那么就可以通过逐步求解单变量微分方程来求解原方程。
这种方法在求解偏微分方程中特别有用,因为它可以将原方程转化为一系列容易求解的常微分方程。
2. 分离变量法在求解拉普拉斯方程中的应用拉普拉斯方程是一种重要的二阶偏微分方程,它在电场、热传导和流体流动等问题中都有广泛的应用。
分离变量法正是一种常用的方法来解决拉普拉斯方程。
通过将方程中的变量分离,得到一系列常微分方程,并求解这些常微分方程,最终可以得到原拉普拉斯方程的解。
3. 贝塞尔函数在分离变量法中的应用在分离变量法中,贝塞尔函数是一种非常常见且重要的特殊函数。
它广泛地出现在圆形和圆柱形边界条件下的分离变量法中。
贝塞尔函数具有良好的性质,对于某些特定的问题有着特别方便的应用。
通过适当地选择边界条件和使用贝塞尔函数的性质,可以简化原方程的求解过程,从而得到更加简洁和优美的解析解。
4. 分离变量法的物理和工程应用我们将讨论分离变量法在物理和工程领域中的具体应用。
以电场分布、热传导问题和匹兹堡问题为例,我们将说明分离变量法是如何应用于这些实际问题中的。
通过分离变量法的应用,我们不仅可以求解这些问题中的微分方程,还可以得到这些问题的具体物理量和工程参数的解析表达式,为实际问题的分析和计算提供了重要的便利。
总结回顾通过本文的讨论,我们深入探讨了分离变量法在求解拉普拉斯方程和应用贝塞尔函数中的原理和应用。
我们分析了分离变量法的基本原理,探讨了其在求解拉普拉斯方程和应用贝塞尔函数时的具体方法,并讨论了其在物理和工程领域中的重要应用。
贝塞尔函数(Bessel functions)是数学上的一类特殊函数的总称。
一般贝塞尔函数是下列常微分方程(一般称为贝塞尔方程)的标准解函数y(x):这类方程的解是无法用初等函数系统地表示的。
贝塞尔函数的具体形式随上述方程中任意实数α变化而变化(相应地,α被称为其对应贝塞尔函数的阶数)。
实际应用中最常见的情形为α是整数n,对应解称为n阶贝塞尔函数。
尽管在上述微分方程中,α本身的正负号不改变方程的形式,但实际应用中仍习惯针对α和−α定义两种不同的贝塞尔函数(这样做能带来好处,比如消除了函数在α=0 点的不光滑性)。
历史贝塞尔函数的几个正整数阶特例早在18世纪中叶就由瑞士数学家丹尼尔·伯努利在研究悬链振动时提出了,当时引起了数学界的兴趣。
丹尼尔的叔叔雅各布·伯努利,欧拉、拉格朗日等数学大师对贝塞尔函数的研究作出过重要贡献。
1817年,德国数学家贝塞尔在研究开普勒提出的三体引力系统的运动问题时,第一次系统地提出了贝塞尔函数的总体理论框架,后人以他的名字来命名了这种函数[1][2]。
现实背景和应用范围贝塞尔方程是在柱坐标或球坐标下使用分离变量法求解拉普拉斯方程和亥姆霍兹方程时得到的(在圆柱域问题中得到的是整阶形式α = n;在球形域问题中得到的是半奇数阶形式α = n+½),因此贝塞尔函数在波动问题以及各种涉及有势场的问题中占有非常重要的地位,最典型的问题有:●在圆柱形波导中的电磁波传播问题;●圆柱体中的热传导问题;●圆形(或环形)薄膜的振动模态分析问题;在其他一些领域,贝塞尔函数也相当有用。
譬如在信号处理中的调频合成(FM synthesis)或凯泽窗(Kaiser window)的定义中,都要用到贝塞尔函数。
定义贝塞尔方程是一个二阶常微分方程,必然存在两个线性无关的解。
针对各种具体情况,人们提出了表示这些解的不同形式。
下面分别介绍这些不同类型的贝塞尔函数。
第一类贝塞尔函数图2 0阶、1阶和2阶第一类贝塞尔函数(贝塞尔J函数)曲线(在下文中,第一类贝塞尔函数有时会简称为“J函数”,敬请读者留意。
bessel函数的数值计算方法《Bessel函数的数值计算方法》一、Bessel函数简介Bessel函数,又称波索函数,是一类非平凡的复变函数,由德国数学家FriedrichBesel(1784-1846)于1817年提出。
它是圆坐标系下的一类无穷级数奇函数,是很多有关理论的的研究的基础,在理论物理学以及工程技术中都具有极为廣泛的应用。
Bessel函数的概念也可以用于解决类似拟辛普森中距离等特征,在椭圆坐标系和球坐标系中具有实际应用价值。
在物理学和数学上,Bessel函数也可以用于分析离散时间序列,例如金融市场分析。
二、Bessel函数的数值计算方法Bessel函数的数值计算有很多种方法,最简单的通过它的级数表示形式计算,例如Jn(x)可以表示为:Jn(x)=∑ (2n+1)!/[(n!) 2^(2n+1) x^n]同时,Bessel函数的数值计算还可以采用递推关系计算,例如Jn+1与Jn的递推关系:Jn+1(x) = Jn(x)(2n/x)-Jn-1(x)此外,Bessel函数可以采用切比雪夫多项式逼近法计算:Jn(x)≈∑A (2n-k)!/[(n-k)! 2^(2n-k) x^(n-k)]k=0,1,...,n在A有已知的值的情况下,可以采用以上的数值计算方法,在计算机科学和电子技术中,我们一般采用计算机软件来快速计算,采用数值积分的方法,具体的可以借助以下几种程序:MATLAB、MATEMATICA、MAPLE等软件。
三、Bessel函数的应用Bessel函数不仅仅有数学意义,它在实际应用中也十分重要,例如镜片折射,它可以用来求解折射方程,解决折射中空圆盘边界,例如折射台球中。
它还可以应用到角度投射视网膜上,用于精确投射光路,精确定位像素位置,从而形成精确的图像。
此外,Bessel函数也可以应用于放射和电磁场的分析,用来求解圆柱坐标放射介质和电荷分布的基础问题,以及球面放射和场的模拟,从而准确计算放射物体的空间分布。
柱坐标拉普拉斯方程分离变量拉普拉斯方程是数学物理中一种常见的偏微分方程,描述了二维和三维空间中的静态场的行为。
柱坐标是一种常见的坐标系,在处理与圆柱形物体相关的问题时非常有用。
本文将介绍柱坐标系下的拉普拉斯方程,以及如何利用分离变量的方法解决这个方程。
柱坐标系简介柱坐标系是一种三维坐标系,由径向r、极角$\\theta$和高度z三个坐标表示。
在柱坐标系中,点的位置可以用$(r, \\theta, z)$表示,其中$r\\geq0$,$0\\leq\\theta<2\\pi$,$-\\infty<z<+\\infty$。
通过极坐标系的转换,我们可以将直角坐标系中的点的坐标转换为柱坐标系中的坐标。
柱坐标系下的拉普拉斯方程在柱坐标系下,拉普拉斯方程可以写作:$$\\frac{1}{r}\\frac{\\partial}{\\partial r}\\left(r\\frac{\\partial\\Phi}{\\partial r}\\right) + \\frac{1}{r^2}\\frac{\\partial^2 \\Phi}{\\partial\\theta^2} + \\frac{\\partial^2 \\Phi}{\\partial z^2} = 0$$其中$\\Phi(r, \\theta, z)$是待求解的场量。
分离变量方法为了解决柱坐标系下的拉普拉斯方程,我们引入一个假设:$\\Phi(r, \\theta, z) = R(r)\\Theta(\\theta)Z(z)$。
将这个假设代入到拉普拉斯方程中,我们可以将方程分解为三个独立的普通微分方程。
首先考虑径向方程$\\frac{1}{r}\\frac{d}{dr}\\left(r\\frac{dR}{dr}\\right) +\\left(\\lambda - \\frac{m^2}{r^2}\\right)R = 0$,其中$\\lambda$和m是待定常数。
贝塞尔函数及其应用题目:贝塞尔函数及其应用摘要贝塞尔方程是在柱坐标或球坐标下使用分离变量法求解拉普拉斯方程时得到的,因此它在波动问题以及各种涉及有势场的问题的研究中占有非常重要的地位。
贝塞尔函数是贝塞尔方程的解。
它在物理和工程中,有着十分广泛的应用。
本文首先通过一个物理问题引入贝塞尔方程,并求出贝塞尔方程的解,即贝塞尔函数。
其次列出了贝塞尔函数的几个重要的结论,如递推公式,零点性质等,并对他们进行了深入的分析。
第二部分主要介绍了傅里叶-贝塞尔级数,通过matlab编程对函数按傅里叶-贝塞尔级数展开之后的图像进行分析,得到了它们的逼近情况。
最后一部分介绍了贝塞尔函数的几个重要应用,一个是在物理光学中的应用,着重分析了贝塞尔函数近似公式的误差;一个是在信号处理中调频制的应用,得到了特殊情况下的公式算法。
关键词:贝塞尔函数,傅里叶-贝塞尔级数,渐近公式目录一、起源1(一)贝塞尔函数的提出1(二)贝塞尔方程的引出1二、贝塞尔函数的基本概念4(一)贝塞尔函数的定义41.第一类贝塞尔函数52.第二类贝塞尔函数73.第三类贝塞尔函数104.虚宗量的贝塞尔函数10(二)贝塞尔函数的递推公式11(三)半奇数阶贝塞尔函数13(四)贝塞尔函数的零点14(五)贝塞尔函数的振荡特性16三、Fourier-Bessel级数16(一)傅里叶-贝塞尔级数的定义16(二)将函数按傅里叶-贝塞尔级数展开17四、贝塞尔函数的应用24(一)贝塞尔函数在光学中的应用24(二)贝塞尔函数在调频制中的应用26附录30一、起源(一)贝塞尔函数的提出随着科学技术的发展,数学的应用更为广泛。
在许多科技领域中,微积分及常微分方程已经不能够满足我们的需要,数学物理方程理论已经成为必须掌握的数学工具。
它们反映了未知函数关于时间的导数和关于空间变量的导数之间的制约关系,同时刻画了物理现象和过程的基本规律。
它的重要性,早在18世纪初就被人们认识。
在1715年,泰勒将弦线的横向振动问题归结为著名的弦振动方程。
bessel函数贝塞尔函数(Bessel Function)是一个在数学中具有广泛应用的特殊函数。
它得名于德国数学家弗里德里希·贝塞尔(Friedrich Bessel),他在19世纪早期首次引入了这个函数。
贝塞尔函数可以分为第一类贝塞尔函数(Bessel Function of the First Kind)和第二类贝塞尔函数(Bessel Function of the Second Kind),分别用J(x)和Y(x)表示。
这两个函数都是解贝塞尔微分方程而得到的,其方程形式为:x^2*y''(x)+x*y'(x)+(x^2-n^2)*y(x)=0其中,y(x)是贝塞尔函数,y'(x)和y''(x)分别表示y(x)的一阶和二阶导数,n是贝塞尔函数的阶数。
第一类贝塞尔函数J(x)在数学和物理学中应用非常广泛。
它在波动现象、电磁场理论、量子力学、光学等领域都有重要的作用。
贝塞尔函数具有周期性,其性质和三角函数类似。
当x趋近于无穷大时,贝塞尔函数的振幅会逐渐减小,并呈现振幅快速振荡的特点。
这种振荡现象在光学中有重要应用,例如描述光的衍射和干涉。
第二类贝塞尔函数Y(x)在数学和物理学中的应用较少,主要用于表示贝塞尔函数的通解形式。
贝塞尔函数的解可以表示为线性组合的形式,其中包括第一类和第二类贝塞尔函数。
第二类贝塞尔函数在x趋近于零的时候有发散的性质,因此在物理问题中较少使用。
除了第一类和第二类贝塞尔函数外,还存在修正贝塞尔函数(Modified Bessel Function),通常用I(x)和K(x)表示。
修正贝塞尔函数在数学分析中也有广泛的应用,特别是在处理边界值问题和椭圆型方程时会经常出现。
贝塞尔函数的计算通常使用数值方法进行,尤其是在高阶贝塞尔函数的计算中。
常用的数值计算方法包括泰勒展开法、渐进展开法、递推关系等。
此外,贝塞尔函数还有一系列的性质和恒等式,如递推关系、积分关系、级数展开等,这些性质可以用于简化贝塞尔函数的计算和分析。
Chapter 13 柱坐标下的分离变量法 Bessel 函数Abstracts :以3+1D 实例通过柱坐标系下方程的变量分离与求解 引入各种柱函数(Bessel 函数、虚宗量Bessel 函数 和球Bessel 函数等)。
在分析这些函数性质的基础上, 表述相应定解问题的物理解。
一、柱坐标下的变量分离1. 柱坐标系下的稳定问题(Laplace 方程)222222110,u u uu zρρρρρϕ⎛⎫∂∂∂∂∇=++= ⎪∂∂∂∂⎝⎭ (1) 即:()22110.zz u u u u ρϕϕρρρρ∇=++= (2)令(,,)()()()u z R Z z ρϕρϕ=Φ,代入(2)得:()20.Z RZ R R Z ρρρΦ''''''+Φ+Φ= (3)2(3)R Zρ⨯Φ得:()2'.R Z RZρρρλ'''''Φ+=-=Φ(4) 分离变量得:()20.(5)'.(6)R Z R Zλρρρλ''Φ+Φ=⎧⎪'⎨''+=⎪⎩(5)与周期性边界条件(0)(2),(0)(2)ππ''Φ=ΦΦ=Φ构成本征值问题。
解得:2 (0,1,2,3,),m m m λ== (){cos ,sin }.m m m ϕϕϕΦ= (6)即为()22'R Z mRZρρρ'''+=⇒分离变量()22'.R m Z RZρμρρ'''-=-=- 得:()2220.0.Z Z R R m R μρρμρ''-=⎧⎪⎨'''++-=⎪⎩这两个方程,先求解哪一个以及μ如何取值,取决于哪一个可构成本征值问题,也就是取决于定解问题的边界条件。
如果()R ρ构成本征值问题,则()2220,R R m R ρρμρ'''++-=式中μ的取值范围不同,方程解的形式与性质不同。
1)0:μ= 220,R R m R ρρ'''+-=即为Euler eq. 2)0:μ>)2220.R R m R ρρ⎡⎤'''++-=⎢⎥⎣⎦(7)记:()()x R y x ρ==⎪⎩则:)d d ()d ,d d d d d ,d d R y x x R x R R y ρρμρρ⎧''==⋅=⎪⎪⎨'⎪'''''====⎪⎩代入(7)得()2220,x y xy x m y '''++-=即为m 阶Bessel eq.3)0:μ< 令2k -=μ,代入()0222=-+'+''R m R R μρρρ得()22220.R R k m R ρρρ'''+-+= (8)记,()()k x R y x ρρ==,代入(8)得:()2220,x y xy x m y '''+-+=即为虚宗量Bessel eq. (9)令:,()()ix t y x t φ==代入(9)得()2220,t t t m φφφ'''++-=即为Bessel eq.2. 柱坐标系下的非稳定问题(振动、输运方程)2222(,)(,)0;(,)(,)0.tt t u r t a u r t u r t a u r t ⎧-∇=⎪⎨-∇=⎪⎩ 令(,,,)()(,,)u z t T t V z ρϕρϕ=,代入上式得:222222;.T Vk a T VT V k a TV ''⎧∇==-⎪⎪⎨'∇⎪==-⎪⎩分离变量得:2222T a k T T a k T ''⎧+=⎪⎨'+=⎪⎩ 和 220,V k V ∇+= 此为Helmholtz 方程,即:()22110.zz V V V k V ρϕϕρρρρ+++= 令(,,)()()()V z R Z z ρϕρϕ=Φ,代入上式得:()22220.0.0.Z Z R R k m R λμρρμρ⎧''Φ+Φ=⎪⎪''-=⎨⎪⎡⎤'''+++-=⎪⎣⎦⎩同样要求对2k μ+的符号()±加以讨论(下面的第二节为正,第三节为负—源于()Z z 的本征值问题)。
二、Bessel 函数 [(圆)柱函数] 1. Bessel 函数设,()(),x R y x ρ==则一般地(ν可以不为整数)()0222=-+'+''y x y x y x ν⇒解()J ()N ()y x A x B x νν=+,其中:()()201J ()!12kkk x x k k ννν+∞=-⎛⎫= ⎪⋅Γ++⎝⎭∑,J ()cos J ()N ()(integer)sin x x x ννννπνπν--=≠,J ()cos J ()N ()lim ( integer).sin n n x x x n ννννπννπ-→-== J ():x νν 阶(第一类)Bessel 函数; N ():x νν阶(第二类)Bessel 函数.J ()J ()J ()N ()N ():Norimann m m m x x m x x x νννν-≠==整数,和线性无关解;整数,和线性无关解,函数。
当=x J ()x ν和N ()x ν都是实函数,现在再引入两个复函数。
(1)H ()J ()N ()x x i x ννν=+,第一种Hankel 函数; (2)H ()J ()N ()x x i x ννν=-,第二种Hankel 函数,它们统称为ν阶(第三类)Bessel 函数,于是Bessel eq. 的解可以是以上四种函数中任何两个的线性组合。
类似于cos ,sin ,cos sin ,cos sin x x x i x x i x +-都是方程()()0y x y x ''+=的特解,其通解可以用以上四个函数中任何两个线性组合表示。
2. 各种柱函数的递推公式与渐近性质(1)递推公式()()11','.x Z x Z x Z x Z νννννννν---+⎧=⎪⎨=-⎪⎩11,.Z Z Z x Z Z Z x νννννννν-+⎧'+=⎪⎪⇒⎨⎪'-=-⎪⎩11112,2.Z Z Z Z Z Z x ννννννν-+-+'=-⎧⎪⇒⎨=+⎪⎩ Z ν代表(1)(2)J ,N ,H ,H νννν.证明:例如,()()201J ()!12kkk x x k k ννν+∞=-⎛⎫= ⎪Γ++⎝⎭∑, ()()()()()221121212001(22)1J 'J ,!12!112k k k k k kk k k x x x x x k k k k νννννννννννν+--+∞∞-+-+==-+-===Γ++Γ-++∑∑即:()1'.x Z x Z νννν-=同理又有:()1'x Z x Z νννν--+=-. 特例:01J'J =-()100J 1J ()xd x ξξ⇒=-⎰,()1-='ννννZx Z x ()1110J J ()xd x x ννννξξξ+++⇒=⎰.(2)渐近性质(A). x 很小(0)x →时,()()201J ()!12kkk x x k k ννν+∞=-⎛⎫= ⎪Γ++⎝⎭∑20J ()~1;21J ()~ (0).(1)2x x x x νννν⎧⎛⎫-⎪ ⎪⎝⎭⎪⇒⎨⎛⎫⎪≠ ⎪⎪Γ+⎝⎭⎩210221(1)!N ()ln J ()2!21(1)111111()!!222m nm m m n m nn mn m x m n x x C x n x n m n n m n πππ-+-=-+-∞=--⎛⎫⎛⎫=+- ⎪ ⎪⎝⎭⎝⎭-⎡⎤⎛⎫⎛⎫⎛⎫-+++++++ ⎪ ⎪ ⎪⎢⎥--⎝⎭⎝⎭⎝⎭⎣⎦∑∑22N ()~ln ~ln ;2()N ()~ (0).2x x x x x ννππννπ-⎧⎪⎪⇒⎨Γ⎛⎫⎪-≠ ⎪⎪⎝⎭⎩ (1)0(1)2H ()~ln ;()H ()~ (0).2i x x x x iννπννπ-⎧⎪⎪⎨Γ⎛⎫⎪-≠ ⎪⎪⎝⎭⎩(2)0(2)2H ()~ln ;()H ()~ (0).2i x x x x i ννπννπ-⎧-⎪⎪⎨Γ⎛⎫⎪≠ ⎪⎪⎝⎭⎩20J ()~12x x ⎛⎫- ⎪⎝⎭⇒0J (0)1=(上述特例积分时用过此).1J ()~ (0)(1)2x x νννν⎛⎫≠ ⎪Γ+⎝⎭⇒J (0)0 (0)νν=≠.可见0x =并非0J ()x 之零点,而是J (0)ν之ν阶零点。
22N ()~ln ~ln ;2N (0) (0).()N ()~ (0)2x x x x x νννππνννπ-⎧⎪⎪⇒=-∞≠⎨Γ⎛⎫⎪-≠ ⎪⎪⎝⎭⎩(1)(2)H (0) (0)νν=∞≠.(B). x 很大()x →∞时 [衰减式震荡函数,证明见教材§13.5]J ()~.24x x ννππ⎛⎫-- ⎪⎝⎭N ()~.24x x ννππ⎛⎫-- ⎪⎝⎭(1)24H ()~.i x x νππν⎛⎫-- ⎪⎝⎭(2)24H ()~.i x x νππν⎛⎫--- ⎪⎝⎭3. Bessel 函数J ()n x 的基本性质(这里仅仅讨论整数阶Bessel 函数)(1)生成函数(母函数)12J() (0).x z n z nn ex z z ⎛⎫∞- ⎪⎝⎭=-∞=<<∞∑特别地,令i z e θ=,有sin ().ix in nn eJx e θθ∞=-∞=∑证明:()120201,!21,!2l x z l l k k x zk k x e z l x e z k -∞=∞--=⎧⎛⎫=⎪ ⎪⎝⎭⎪⎨-⎪⎛⎫= ⎪⎪⎝⎭⎩∑∑则: ()()()()()()11222000012200011!!211!!2!!211!()!2!()!2kl kxx x z z z z l kl k kkl kl kl k l kk l k l k l kl nn kn ln nk n l n x ee ez l k x x z z l k l k x x z z k k n l l n --+∞∞---==++∞∞∞∞--====+-+-+∞∞∞-∞====--⎛⎫=⋅= ⎪⎝⎭--⎛⎫⎛⎫=+ ⎪ ⎪⎝⎭⎝⎭--⎛⎫⎛⎫=+ ⎪ ⎪+-⎝⎭⎝⎭=∑∑∑∑∑∑∑∑∑∑011J ()(1)J ()J ()(1)(1)J ()J().nnnnn n nn n n n n n n n n nn x z x z x z x z x z ∞-∞∞-∞----==-==-∞=-∞+-=+--=∑∑∑∑∑(1’)平面波按柱面波的展开cos 1()201J()J ()2J ()cos().i ikz ik x t ie t mtmx k m m m m e e ex t x k i m θρθρρθ∞=-==-∞∞=====+∑∑(2)加法公式 J ()J()J ().n kn k k x y x y ∞-=-∞+=∑证明:()12J(),x yz z n nn ex y z -+∞-=-∞=+∑又()()()()()()()111222J J J J .x yxyz z z z z z nln ln l n ln l n l eeex zy z x y z ---+∞∞∞∞---+=-∞=-∞=-∞=-∞=⋅==∑∑∑∑ 令l n k +=,则:()()()()()12J J J J .x yz z knnk nkn kn k n k ex y zx y z -+∞∞∞∞---=-∞=-∞=-∞=-∞==∑∑∑∑所以J ()J()J ().n kn k k x y x y ∞-=-∞+=∑(3)积分公式由sin J()ix in nn ex e θθ∞=-∞=∑得()()()sin sin cos cos 11J ()d d 2211cos sin d cos sin d 22d d .22ix in ix in n nnix in ix in x e e e x n n x i i ee ππθθθθππππππππθθθθππθθππθθθθθθππθθππ--+----+-+--===-=--==⎰⎰⎰⎰⎰⎰其中倒数第二等式的推导用了展开中的/2.θθπ→+ (4)J ()n x 的零点[方程J ()0n x =的根](A). J ()n x 的零点有无限多个,且0x ≠的零点都是一级零点()(1,2,3,)nx n ν=: J ()~24x x ννππ⎛⎫-- ⎪⎝⎭.0x =为J ()n x (0n ≠)的n 级零点:1J ()~ (0)(1)2x x νννν⎛⎫≠ ⎪Γ+⎝⎭.(B). J ()n x 的零点必正负成对:这是因为J ()n x 具有奇(偶)对称性,即J ()(1)J ()n n n x x -=-,因此可以只讨论正零点。