整数阶贝塞尔函数算法探讨与研究.doc

整数阶贝塞尔函数算法探讨与研究摘要:贝塞尔函数在数学应用中十分广泛，通过对整数阶贝塞尔函数的性质及不同变量取值范围时的递推算法的计算分析，使计算结果更加准确的接近理论值，并通过仿真拟合验证递推结果的正确性。关键词：整数阶，贝塞尔函数，递推算法，大变量，小变量。The approach and research to the Bessel function of integral orderabstract: The Bessel function is used widely in mathematics, through analyzing the character of the Bessel function of integral order and the recursion in different variable dereferencing zone, we can get more exact numerical value ,which can approach to the theoretical value. Finally, the correctness of result is verified through fitted to theoretical value.1、引言贝塞尔函数是一种复杂的特殊函数，它最早应用于悬链振荡、长圆柱体冷却以及紧张膜振荡的物理问题中1,2。近年来随着轴频电场计算方法的不断完善与改进，在时谐源的建模求解电磁场强度表达式时经常会应用到贝塞尔函数的求解与换算，如时谐电偶极子及时谐电流元。为了进一步改进电场模型及其表达式，对贝塞尔函数的计算是十分必要的。本文从整数阶的贝塞尔函数入手利用收敛较快的递推公式对其进行计算，得出较为准确的近似解。2、整数阶贝塞尔函数整数阶贝塞尔函数的定义式为 (1)式中为整数。整数贝塞尔函数满足下述微分方程 (2)贝塞尔函数具有以下性质3,4：表1 贝塞尔函数性质3、 小变量时整数阶贝塞尔函数不论变量x是实数还是复数，式(1)对于任意x都是收敛的5，为了说明这个问题，以下取进行分析。令 (3)则 (4)不论阶数n的大小6，当时恒有。此时可通过计算级数式(4)的前若干项来逼近的值： (5)式中，当给定误差后，便可确定正整数的大小。知道的数值后，利用下面递推公式便可写出式(5)的递推计算式：递推初始值：递推关系：，其中计算结果：。以计算为例，取，经计算t=0.64,N=3。递推过程如下：，从而得出。实际上，递推值的真实误差为。可见，用以上递推公式计算小变量整数阶贝塞尔函数时收敛很快。4、大变量时整数阶贝塞尔函数以下分两种情况讨论时整数阶贝塞尔函数的计算方法。4.1 时的计算式由贝塞尔函数的性质 (6)假定是常数，当从n增加的方向递推上式时，其误差满足以下差分方程 (7)其特征方程为 (8)因,所以特征方程的解均是特征方程的单根，符合稳定性条件，这说明，当时用式(6)施行正向递推（即沿n增大的方向递推）可以计算。4.2 时的计算式由于时正向递推公式是不稳定的，所以可通过反向递推来完成，即 (9)对于给定的x，当时，所以可以取一大偶数，并令，应用式(9)即可得出计算结果。其表达式如下 (10)式中是整数阶贝塞尔函数的阶数，是计算精度的有效位数，符号表示取数值的整数部分。5、数值分析 应用本文提出的方法得到的数值结果同由其解析表达式求得的结果相对照7，如图1所示，所得到的数值结果与解析值完全吻合。 图 1 快速数值解和理论值的比较6、结论针对整数阶贝塞尔函数，提出了相对快速的收敛计算方法，由于贝塞尔函数的复杂性，分为不同情况下的递推计算公式，同时利用该方法的退到结果同解析表达式求出的结果进行拟合比较，结果完全相同，为以后的贝塞尔函数的计算研究提供了一定的理论依据。参考文献图 42 快速数值解和理论值的比较, ， f ()理论值数值解1 Watson G. N. A. Treatise on the theory of Bessel function, 2nd edM. Cambridge, England: Cambridge University Press,19662 莫平华. 一阶贝塞尔函数广义积分的数值计算J.数学理论与应用，2007,27(1):65-67.3 张爽，郭欣，宋立军.利用贝塞尔函数的级数形式进行数值计算的误差分析J.长春大学学报，2004，14(2):57-59.4 魏彦玉，宫玉彬，王文祥.任意阶复宗量贝塞尔函数的数值计算J.电子科技大学学报，1998,27(2):172-175.5 伍刚.贝塞尔函数的计算机仿真