摘要
基于珚(P_(nm)(cosθ))/(sinθ)m(>0)的非奇异递推公式,给出了基于球坐标的引力矢量和垂线偏差非奇异计算公式;针对极点λ可任意取值引起的地方指北坐标系方向的不确定性问题,证明了引力矢量在转换到地心直角坐标系后不随λ的变化而变化,即与λ的取值无关。最终的数值计算结果表明,直角坐标系下的非奇异计算公式与本文提出的球坐标下的非奇异计算公式所得计算结果绝对值差异小于10^(-16) m/s^2,证明了该非奇异公式的正确性。最后总结了所有引力位球函数一阶导、二阶导非奇异性计算的一般思路。
When computing gravity vector and vertical deviation using spherical harmonic function,singular problem exists in the formulas expressed by spherical coordinates.This will cause some errors in gravity vector and vertical deflection data and influence their application.This paper aims at proposing an alternative method to solve this problem.Based on the non-singular expression of （Pnm（cosθ））/（sinθ）m（〉0）,the paper gives the non-singular formulas expressed by spherical coordinates for computing gravity vector and vertical deviation.At North and South Poles,the paper proves that even values ofλare arbitrary,the values of gravity vector are sole when the values are transferred to Earth fixed rectangular coordinate system.In order to show the validity of our method,the paper computes gravity vectors at pointsθ=0,λ=i/3602π（i=0,1...359）using the former 100 degrees and orders of EGM2008.The absolute differences between the computing results by our method and the non-singular formula expressed in Cartesian coordinates are smaller than 10-16 m/s2,which show the validity of our method.The non-singular expression based on spherical function derived by the paper can make full use of the high accurate algorithms of Legendre function,so the proposed method has better generality ability compared with the non-singular formula expressed by Cartesian coordinates.Finally the methods for non-singular computing of all the first or second derivations of gravity field potential are summarized.The method of this paper can also be directly applied to the non-singular calculation of the spherical harmonic model of magnetic field,and the basic idea is similar to that of this paper.
作者
朱永超
万晓云
于锦海
ZHU Yongchao;WAN Xiaoyun;YU Jinhai(College of Earth Science, University of Chinese Academy of Sciences, Beijing 100049, China;Qian Xuesen Laboratory of Space Technology, Beijing 100049, China)
出处
《武汉大学学报（信息科学版）》
EI
CSCD
北大核心
2017年第12期1854-1860,共7页
Geomatics and Information Science of Wuhan University
基金
国家自然科学基金(41404019
41674026
41774089)
国家重点开发计划(2016YFB0501702)
中国科学院太空应用重点实验室开放基金(CSU-WX-A-KJ-2016-044)~~
关键词
引力
垂线偏差
极点
非奇异
球坐标
gravity
vertical deviation
North and South Poles
non-singular
spherical coordinates
