GPS卫星位置的计算

一、用广播星历计算卫星位置

1. 计算卫星运动的平均角速度n
   

2. 计算观测瞬间卫星的平近点角M
   
   广播星历每2h更新一次，将参考时刻设在中央时刻，外推间隔小于等于1h，而卫星的运行周期为12h左右，采用卫星过近地点时刻来计算时，外推间隔最大有可能达6h，用参考时刻来取代卫星过近地点时刻后，外推间隔将大大减小，也能获得精度较高的结果。

3. 计算偏近点角
   

4. 计算真近点角f
   
   式中，e为卫星轨道的偏心率，由广播星历给出，所以
   

5. 计算升交角距
   

6. 计算摄动改正项
   

7. 对u、r、i进行摄动改正
   

8. 计算卫星在轨道面坐标系中的位置
   在轨道平面直角坐标系中（坐标原点位于地心，X轴指向升交点）卫星的平面直角坐标为：
   {x=rcosu, y=rsinu

9. 计算观测瞬间升交点的经度
   

10. 计算卫星在瞬时地球坐标系中的位置
    

11. 计算卫星在协议地球坐标系中的位置
    
    

二、用精密星历计算卫星位置

精密星历是按一定的时间间隔来给出卫星在空间的三维坐标、三维运动速度及卫星钟改正数等信息。观测瞬间的卫星位置及运动速度可采用内插法求得，其中拉格朗日多项式内插法被广泛采用