从广播星历计算卫星位置：
1. 计算卫星运动的平均角速度n
首先根据广播星历中给出的参数A计算出参考时刻TOE的平均角速度n0：
n0?GM(A)3，
式中，GM为万有引力常数G与地球总质量M之积。然后根据广播星历中给出的摄动参数?n计算观测时刻卫星的平均角速度n：n?n0??n。 2. 计算观测瞬间卫星的平近点角M:
M?M0?n(t?TOE)
式中，M0为参考时刻TOE时的平近点角，由广播星历给出。 3. 计算偏近点角E：
E?M?esinE
解上述方程可用迭代法或微分方程改正法。 4. 计算真近点角f
cosE?e?cosf???1?ecosE ?2?sinf?1?esinE?1?ecosE?
式中，e为卫星轨道的偏心率，由广播星历给出。 5. 计算升交距角u?：
u????f
式中，?为近地点角距，由广播星历给出。
6. 计算摄动改正项?u,?r,?i：
广播星历中给出了下列6个摄动参数：Cuc,Cus,Crc,Crs,Cic,Cis，据此可以求出由于J2项而引起的升交距角u的摄动改正项?u、卫星矢径r的摄动改正项?u和卫星轨道倾角
i的摄动改正项?i。计算公式如下：
??u?Cuccos2u??Cussin2u????r?Crccos2u??Crssin2u? ???Ccos2u??Csin2u?icis?i7. 计算u?,r?,i0进行摄动改正
?u?u???u???r?r???r?a(1?cosE)??r ?dii?i???(t?TOE)0i?dt? 式中：a为卫星轨道的长半径，a?(A)2，i0为TOE时刻的轨道倾角，由广播
星历中的开普勒六参数给出，
di为i的变化率，由广播星历中的摄动九参数给出。 dt 8. 计算卫星在轨道面坐标系中的位置
在轨道平面直角坐标系中（坐标原点位于地心，X轴指向升交点）卫星的平面直角坐标为：
?x?rcosu ?y?rsinu? 9. 计算观测瞬间升交点的经度L
?，那么观测瞬 若参考时刻TOE时升交点的赤经为?TOE，升交点对时间的变化率?间t的升交点赤经?应为：
?(t?TOE) ???TOE???可从广播星历的摄动参数中给出。 ? 设本周开始时刻（星期日0时）格林尼治恒星时为GASTweek，则观测瞬间的格林尼治恒星时为：
GAST?GASTweek??et
式中：?e为地球自转角速度；t为本周内的时间(s)。 这样就可求得观测瞬间升交点的经度值为：
?L???GAST??TOE?GASTweek??(t?TOE)??et
? 令?0??TOE?GASTweek，则有：L??0??(t?TOE)??et
10. 计算卫星在瞬时地球坐标系中的位置
已知升交点的大地经度L以及轨道平面的倾角i后，就可通过两次旋转方便地求得卫星在地固坐标系中的位置：
?X??x??xcosL?ycosisinL????????Y??RZ(?L)RX(?i)?y???xcosL?ycosisinL? ?Z??0???ysinL??????11. 计算卫星在协议地球坐标系中的位置
观测瞬间卫星在协议地球坐标系中的位置
?x??X??1??????RY(xP)RX(?yp)??Y???0?y??z??Z???x??CTS???p01ypxp??X?????yp??Y?
?Z?1????