一、三类柱函数
前面已经求出阶贝塞尔方程的通解,即:
或者是:
其中第一个通解不适用于整数阶贝塞尔方程,但是第二个对整数阶的情况照样适用。通常又取线性独立的两个函数:
这两个函数叫做第一种和第二种汉克尔函数,于是阶贝塞尔方程的通解又可以表示为:
贝塞尔函数、诺依曼函数、汉克尔函数分别叫做第一类、第二类、第三类柱函数。
前面提到过,当时,都不是有限的,由于圆柱内部的自然边界条件,应当舍弃诺依曼函数和负阶贝塞尔函数,只要零阶和正阶的贝塞尔函数。当时,根据柱函数的渐进公式(后面会证):
都趋向于,均不可舍弃。
下面推导三类柱函数的递推公式,以贝塞尔函数为例,根据其定义,得:
令,则:
仿照这种方法,可以推导出:
上面两式都是贝塞尔函数的线性关系式,由于诺依曼函数时正负阶贝塞尔函数的线性组合、汉克尔函数是贝塞尔函数和诺依曼函数的线性组合,所以上面的两式也适用于诺依曼函数和汉克尔函数。总之,用表示三类柱函数,总有:
将上面两式展开,并联立化简,可以得到递推公式:
二、贝塞尔方程
1、贝塞尔函数与本征值问题
对于圆柱的内部问题,如果柱侧面有齐次边界条件,则贝塞尔函数的情况应该排除,因为导出的虚宗量贝塞尔方程的解为,其中由于处有限,故,又由于恒正,故解恒不为,得到,为平凡解,除非。(对于拉普拉斯方程的虚宗量贝塞尔方程的通解为:)
这样,只需要考虑的情况,在这种情况下,应该是整数阶贝塞尔方程的解,由于圆柱轴上的自然边界条件,这两个方程的两个线性独立解之中,我们只要负阶贝塞尔函数:
柱侧的边界条件决定了的可能值,这些就是整数阶贝塞尔方程在所给边界条件的本征值,本征函数为,下面分别讨论不同类型的边界条件。
(1) 第一类齐次边界条件
第一类齐次边界条件,为圆柱的半径,代入到解中得到,这是一个超越方程,使用工具求解。本征值:
其中是的第个正零点,而表示在满足第一类齐次边界条件下方程:
的第个正根,故。为什么一定是正根呢,首先,如果,则原方程退化为简单常微分方程,其解为:
在第一类齐次边界条件下,只有平凡解,无意义。所以,不是本征值。
NOTE
贝塞尔函数的零点可以通过下面公式计算:
其中:
非常逆天。同时,特殊位置的贝塞尔函数存在下面关系:
即时,贝塞尔函数近似于三角函数。
由于:
则贝塞尔函数零点正负成对,绝对值相同。
由于:
则贝塞尔函数有无穷多个零点。
由于:
则的绝对值最小的零点比的绝对值最小的零点更接近于零。
由于:
则与的零点两两相间。
(2) 第二类齐次边界条件
第二类齐次边界条件,如果,则原方程退化为简单常微分方程,其解为:
在第二类齐次边界条件下,只有平凡解,无意义。所以,不是本征值。在的情况下,这个条件即为,所以本征值:
记,其中是在第二类齐次边界条件下方程:
的第个正根,即的第个正的零点。在的特例下:
也就是,其他情况下,根据递推公式:
得知的零点可以从和的零点得出。
NOTE
贝塞尔函数导数的零点可以使用下面公式计算:
其中:
非常逆天。
(3) 第三类齐次边界条件
第三类齐次边界条件形如,记,这个条件即为方程:
第三类齐次边界条件下的本征值为:
其中是方程的第个正根。这玩意你就去查表吧。
NOTE
后面统一使用表示的第个正零点,表示第类齐次边界条件下的第个正根或者本征值。不指明齐次边界条件种类的话,记本征值为或者。
2、贝塞尔函数的正交关系
贝塞尔函数是施图姆-刘维尔本征值问题正交关系的特例,对应于不同本征值的同节贝塞尔函数在区间上带权重正交:
3、贝塞尔函数的模
下面计算贝塞尔函数的模。
记,则:
然后将右边积分使用贝塞尔方程改写:
最后得到:
分别代入不同的齐次边界条件类型:
-
第一类齐次边界条件:
-
第二类齐次边界条件:
-
第三类齐次边界条件:
4、傅里叶-贝塞尔级数
根据施图姆-刘维尔本征值问题性质,本征函数组是完备的,可以作广义傅里叶级数展开的基。在区间上的函数的傅里叶-贝塞尔级数是:
TIP
计算时,可以使用下面的公式,分别是对之前结论的应用:
对于的情况,有傅里叶-贝塞尔积分:
5、母函数
将和分别展开为绝对收敛级数,逐项相乘得到:
上式前一个求和的内为阶贝塞尔函数,后面一个是,即,所以:
所以叫做贝塞尔函数的母函数,令,上式改写为:
又令:
再令:
上述公式都是等价的。
6、积分表达式
把母函数等式的右边看作是复数形式的傅里叶级数,那么就是的傅里叶系数,所以:
利用复数的三角函数展开,化简积分,得到:
这就是贝塞尔函数的积分表达式,使用变量代换可以得到其他形式。
7、诺依曼函数
研究圆柱内部问题时,圆柱轴上的自然边界条件排除了诺依曼函数,但是如果是空心圆柱,不涉及处的边界条件,就不能排除诺依曼函数了。例如对于定解问题:
首先利用,将边界条件化为齐次:
因为是当前是在极坐标系,所以不考虑,问题也与无关,所以,利用:
其中,是特征函数。代入边界条件:
该齐次线性方程组有非零解的条件时系数行列式为零,所以:
这个方程的根可以查表得出,然后求解得到的比值:
将上面的本征解叠加,得到:
上面的,尚未求出,为了求出,要代回初始条件:
将上式右边以为基展开(可以证明他是正交的),比较两边系数得出。算完之后反代得到的表达式。
8、汉克尔函数
波动方程在柱坐标系中分离变量形式的解为:
或
这里将和表示为了复数形式。在较大时,取时间因子为,并省略有关的因子,利用渐进公式,有:
上面第一个式子是朝着方向传播的波,即向外发散的波,第二个式子是朝着方向传播的波,即向内会聚的波。但是贝塞尔函数和诺依曼函数中是分离的,是驻波。
如果取时间因子为,则上面第一个式子是朝着方向传播的波,即向内会聚的波,第二个式子是朝着方向传播的波,即向外发散的波。贝塞尔函数和诺依曼函数仍然是驻波。
这样研究波的发散问题,使用汉克尔函数研究比较方便。例如一个半径为的长圆柱面,其径向速度分布为:
试求解柱面发射的恒定声振动的速度势,设远小于声波波长。该定解问题为:
该问题与无关,所以没有且,边界条件的时间因子,其实部是,能够得到问题的恒定解为:
且,表示波速,表示波数。由于这是一个声波发射问题,已知时间函数是,因此只取,舍弃,在这个例子中,只有一个值,所以无需叠加,故结果形式为:
将结果代入边界条件,由于,所以很小,根据汉克尔函数的表达式(前面找),当参数趋向于零时,利用泰勒展开,可以将上式化为:
即:
回代,得到答案:
当较大的时候使用渐进公式,并取实部:
这是一个振幅按照减小的柱面波。
三、虚宗量贝塞尔方程
前面研究的圆柱状区域的拉普拉斯方程定解问题都是柱侧面有齐次边界条件的,这种问题只需要考虑的情况,如果圆柱上下底面边界条件齐次,但是侧面边界条件非齐次,这时齐次方程和上下底面的齐次边界条件构成本征值问题,其中,应该考虑的解,的情况比较简单,当时,则引至虚宗量贝塞尔方程:
方程的解为阶虚宗量贝塞尔函数:
对于整数,,不是线性独立的,需要寻找另一个解,根据:
有:
继续变形,得到:
右边乘以,使其变为实值函数,记作,即虚宗量汉克尔函数:
于是阶虚宗量贝塞尔方程的两个线性解即为虚宗量贝塞尔函数和虚宗量汉克尔函数。令趋向于整数,有:
化简右边的式子,得到:
这个极限就是整数阶虚宗量汉克尔函数,记作,当时,第一个有限求和项不存在。
当时:
这样,如果所研究的区域包括圆柱轴线,即,那么在轴线上的自然边界条件应该排除虚宗量汉克尔函数。当时,引入渐渐公式,对于较大的,有:
所以如果研究区域伸向无穷远,应该排除。
四、球贝塞尔方程
使用球坐标对亥姆霍兹方程进行分离变量,得到球贝塞尔方程:
把自变数和函数分别变换为和:
则方程化为阶贝塞尔方程:
如果,则方程退化,两个线性独立解为,比较简单,下面讨论的情况。
1、线性独立解
阶贝塞尔方程有如下几种解
其中任取两个就组成阶贝塞尔方程的线性独立解。这样,球贝塞尔方程的线性独立解也就是下列五种之中任取的两种:
2、递推公式
用 代表球贝塞尔函数、球诺伊曼函数及球汉克尔函数,即:
利用之前的方法,可以得到如下几个递推公式:
3、初等函数表示式
根据前面得到的:
有:
于是反复使用递推公式,可以得到的初等函数表达式。
根据诺依曼函数定义:
于是得到:
反复使用递推公式,可以得到的初等函数表达式。
由球汉克尔函数的定义:
可以直接得到球汉克尔函数的初等函数表达式。
4、极限取值
当时,对于球贝塞尔函数:
把理解为非负整数,该级数只含有的正幂次项,故
对于球诺依曼函数:
既有的正幂次项,也有负幂次项,所以,当。
当时,使用渐进公式:
从波动方程分离出来的时间因子是,如果把时间因子附在上面的渐进公式中时,对应驻波,对应朝方向传播传播,即从球坐标系极点向外发散,对应朝方向传播传播,即从球坐标系极点向内会聚。如果把时间因子附在上面的渐进公式中时,对应驻波,对应朝方向传播传播,即从球坐标系极点向内会聚,对应朝方向传播传播,即从球坐标系极点向外发散。
5、球形区域内的本征值问题
球贝塞尔方程是施图姆-刘维尔型方程,即:
其中为本征值,为权重系数。方程在处有自然边界条件,所以取,舍弃,并且应当满足球面上的边界条件,这决定了本征值。对应不同本征值的本征函数在上带权重正交:
本征函数族是完备的,可以作为广义傅里叶级数展开的基:
系数
式中是本征函数的模: