以下计算方法适合于GPS L1 NAV星历 、BDII代 D1星历,其中:
是地心引力常数,
是地球自转角速率,其值的大小参考对应的ICD文档。
1、广播星历参数表
参考时间:
轨道长轴平方根:
偏心率:
近地点幅角:
卫星平均运动速率与计算值之差:
参考时刻平近点角:
参考时刻升交点赤经:
升交点赤经变化率:
参考时刻轨道倾角:
轨道倾角变化率:
轨道改正项参数:
2、计算卫星在ECEF坐标系下的位置坐标
(1)计算:
(2)计算卫星的平均角速率n:
(3)计算平近点角:
(4)计算偏近点角(迭代计算):
(5)计算真近点角:
(6)计算升交点角距:
(7)计算摄动校正项:
(8)计算摄动校正后的升交点角距:
(9)计算摄动校正后的矢径长度:
(10)计算摄动校正后的轨道倾角:
(11)计算卫星在轨道面上的位置:
(12)计算升交点赤经:
(13)计算卫星在ECEF坐标系下的位置:
3、计算卫星在ECEF坐标系下的速度
(1)计算平近点角对时间的一阶导数:
(2)计算偏近点角对时间的一阶导数:
(3)计算真近点角的一阶导数:
(4)计算升交点角距的一阶导数:
(5)计算摄动校正项的一阶导数:
(6)计算摄动校正后的升交点角距的一阶导数:
(7)计算摄动校正后的矢径长度的一阶导数:
(8)计算摄动校正后的轨道倾角的一阶导数:
(9)计算卫星在轨道面上的速度:
(10)计算升交点赤经的一阶导数:
(11)计算卫星在ECEF坐标系下的速度
4、计算卫星在ECEF坐标系下的加速度
(1)计算平近点角对时间的二阶导数:
(2)计算偏近点角对时间的二阶导数:
(3)计算真近点角的二阶导数:
(4)计算升交点角距的二阶导数:
(5)计算摄动校正项的二阶导数::
(6)计算摄动校正后的升交点角距的二阶导数:
(7)计算摄动校正后的矢径长度的二阶导数:
(8)计算摄动校正后的轨道倾角的二阶导数:
(9)计算卫星在轨道面上的加速度:
(10)计算升交点赤经的二阶导数:
(11)计算卫星在ECEF坐标系下的加速度
5、计算卫星在ECEF坐标系下的加加速度