[发明专利]一种卫星位置解算方法

说明书

技术领域

本发明涉及卫星定位技术领域，特别是涉及一种卫星位置解算方法。

背景技术

研究GLONASS（GLOBAL NAVIGATION SATELLITE SYSTE，格洛纳斯）卫星系统定位，尤其是多系统组合定位有重要实用价值，如何根据GLONASS广播星历进行卫星位置计算是组合定位的重要一环。GLONASS每30min更换一组新的星历参数，采用数值积分的方法求解某一时刻t_e的卫星速度和位置时，要以最新星历的参考时刻t_b为参考时刻，并以该时刻卫星的坐标和速度为起算数据，对加速度和速度进行积分，其积分通式如（1-1）所示：

Vte=Vb+∫tbteadtXte=Xb+∫tbteVdt---(1-1)]]>

其中，V、X表示速度和位置，a表示加速度。目前GLONASS卫星轨道数值积分方法主要有龙格-库塔法（以下简称R-K法），其中以四阶龙格-库塔法最为常用。但是，龙格-库塔法在计算某一时刻卫星位置和速度时只用到前一步的信息，为了保证定位精度，需要多次重新计算多个点处的函数值，计算量较大。

使用R-K方法每进行一步积分，即由t_i时刻的速度和位置（V_i，X_i）计算计算t_i+1时刻的速度和位置坐标（V_i+1X_i+1），都要按照公式（1-2）进行计算，其中h＝t_i+1-t_i为步长，a₁,a₂,a₃,a₄为加速度中间变量值，V₁，V₂，V₃，V₄为速度中间变量值，为日月摄动加速度。