关键词:
历表
数据结构
算法设计
伪波
基准序列
精度
摘要:
基于引潮位Doodson规格化展开的基本原理以及ELP/MPP02月球历表、Newcomb太阳历表的结构,设计了一个用以代表三角函数的数据结构;并由此自定义了三角函数的乘法与加法算法,将月球、太阳引潮位分别展开至5阶、3阶,振幅绝对值截断阈值为10-7,在展开过程中对"伪波"进行"滤波"处理后,最终得到包含4 686项展开式的引潮位展开表(其中振幅绝对值超过10-6的展开式有2 441项)。以德国BFO(Black Forest Observatory)测站为例,基于DE431历表,通过天球参考系变换计算得到1950—2050年间时间间隔为1 h的法向引潮力基准序列BFDE431;并根据各个引潮位展开表计算得到相应的法向引潮力序列;通过与BFDE431基准序列求差,得到各个差值序列的统计值。计算结果表明,文中给出的展开表对应的差值序列数值在±58×10^(-11)m·s^(-2)(nGal)以内,均方差为12.5×10^(-11)m·s^(-2),与XI89展开表的精度相当。但由于没有考虑行星及地球扁率的影响,仍未达到HW95和RATGP95展开表的精度水平。