整数阶复宗量贝塞尔函数的计算程序研究  

Program for calculating the integer order of Bessel functions with complex arguments

作　　者：任宏红 郭迎春 王兵兵[2] REN Hong-hong;GUO Ying-chun;WANG Bing-bing(School of Physics and Materials Science,East China Normal University,Shanghai 200241,China;Laboratory of Optical Physics,Beijing National Laboratory of Condensed Matter Physics,Institute of Physics,Chinese Academy of Sciences,Beijing 100190,China)

机构地区：[1]华东师范大学物理与材料科学学院,上海200241 [2]中国科学院物理研究所凝聚态物理国家重点实验室光物理实验室,北京100190

出　　处：《华东师范大学学报（自然科学版）》2019年第1期76-82,92,共8页Journal of East China Normal University(Natural Science)

基　　金：国家自然科学基金(61275128;11474348)

摘　　要：鉴于目前的算法程序集中没有现成的计算复宗量贝塞尔函数的程序,本文基于贝塞尔函数的逆向递推关系编写了计算整数阶复宗量第一类贝塞尔函数的Fortran程序源代码与Matlab软件的计算结果比较,两者至少有12位有效数字一致.接着运用此程序,分析了徐士良的《FORTRAN常用算法程序集》中的纯虚宗量的贝塞尔函数,即变形贝塞尔函数程序的准确度,发现其准确度为6位有效数字.最后,对基于实宗量贝塞尔函数和纯虚宗量贝塞尔函数相乘然后用无限求和来计算复宗量贝塞尔函数值的方法的准确性进行了探讨.证明其仅能对有限的贝塞尔函数进行准确计算.这是由于当求和项中有远大于最终的求和项时,会导致求和结果的有效数字减少甚至完全错误.Fortran source code for calculating the integer order of Bessel functions of the first kind with complex arguments is presented.The method is based on the backward recurrence relation of Bessel functions.Values of the Bessel function generated by our program are in-agreement with the values generated by Matlab to at least 12 significant digits.We use the program to calculate the integer order of Bessel functions of the first kind with pure imaginary arguments,provided by Xu Shiliang’s Fortran algorithm assembly.The results show that the first 6 significant digits are accurate.We also analyze the algorithm for calculating Bessel functions with complex arguments,which use the infinite sum of the product of the real arguments of the Bessel function and the pure imaginary arguments of the Bessel function,provided by Xu Shiliang’s algorithm assembly.The results show that this algorithm does not always get accurate values for Bessel functions with complex arguments.The reason lies with the fact that the term in the sum larger than the function value causes the loss of significant digits.

关 键 词：复宗量贝塞尔函数 Fortran源程序 逆递推计算 

分 类 号：O411.2[理学—理论物理]