超高阶次Legendre函数的跨阶数递推算法  

A recursion arithmetic formula for Legendre functions of ultra-high degree and order on every other degree

作　　者：于锦海[1,2] 曾艳艳[1,2] 朱永超[1,2] 孟祥超[1,2] 

机构地区：[1]中国科学院计算地球动力学重点实验室,北京100049 [2]中国科学院大学地球科学学院,北京100049

出　　处：《地球物理学报》2015年第3期748-755,共8页Chinese Journal of Geophysics

基　　金：国家高技术研究发展计划(863计划)项目(2013AA122502-2);国家自然科学基金项目(41274034);CAS/CAFEA国际创新团队项目(KZZD-EW-TZ-19)联合资助

摘　　要：本文引入了Legendre函数的跨阶数递推算法,并利用该算法在双精度数范围内计算了按间隔为1°余纬从1°变化至89°对应的直到完整的20000阶次的归一化连带Legendre函数的值.为验证计算精度,通过多种途径对该算法的计算结果进行检验,结果表明:该算法算得的每个阶次连带Legendre函数的值至少具有10-10这样的绝对精度.此外还对该算法的计算用时进行了统计,结果为该算法的计算用时大约是Legendre函数计算中常用的按阶数递推算法用时的1.6倍.The most popular arithmetic for Legendre functions in the study of the gravity field is the increasing degree recursion method（IDR）.Although IDR has a simple expression in mathematics,it is not suitable to compute Legendre functions of ultra-high degree and order.For example,when degree reaches 1800,IDR cannot be used to compute Legendre functions because of under-flow in the double floating-point range.Hence,some modified versions for IDR are discussed.A typical modification is to extend the double floating-point range,and the result is that the run-time in computation increases rapidly too.In order to solve the problem in computing Legendre functions of ultra-high degree and order,a recursion arithmetic approach on every other degree for Legendre functions is presented in this paper.Our aim is to illustrate that this approach can be not only used to compute Legendre functions,but the run-time of computation is also saved.The main method is computation,that is,the values of Legendre functions are computed based on the recursion formulas of Legendre functions on every other degree,and then the computation accuracies and the run-time are assessed.The values of the fully normalized associated Legendre functions up to complete degree and order 20000 are computed from colatitude 1°to 89°with 1°interval in the double floating-point range.The computation accuracies are estimated according to the properties of Legendre functions.From our computation results,it can be acclaimed that the computation accuracies of Legendre functions up to complete degree and order 20000 can reach 10-10 at least if the recursion arithmetic on every other degree is applied.The run-time for computation is listed.From the statistical results,the run-time of the arithmetic in this work is approximately 1.6times that of IDR.The computation of Legendre functions of ultra-high degree and order plays an important role in refining the gravity field.The recursion arithmetic on every other degree is suitable to compute Legendre functions of ult

关 键 词：Legendre函数 递推算法 双精度数 

分 类 号：P312[天文地球—固体地球物理学]