自动积分步长的GLONASS卫星轨道龙格库塔积分法

摘要

为了快速准确地计算GLONASS卫星坐标,通过对不同的卫星状态方程进行比较分析,给出一个准确的GLONASS卫星运动状态方程.以该状态方程为基础,提出了古典形式、龙格-库塔、基尔公式3种不同形式的龙格-库塔轨道积分方法,比较分析得出3种积分方法精度相当,但基尔公式的舍入误差较小.然后,以基尔公式为基础,给出了定步长和自动选择步长2种不同的积分方法的详细步骤.讨论了定步长的选择以及自动选择步长收敛值的取值,得到其最佳值在20~30之间.最后,利用GLONASS的广播星历对基于基尔公式的自动选择积分步长的龙格-库塔法进行精度分析,积分时间间隔为60min,积分误差小于3m.