从广播星历计算卫星位置

从广播星历计算卫星位置：1.计算卫星运动的平均角速度n首先根据广播星历中给出的参数A计算出参考时刻TOE的平均角速度0n：30)(AGMn，式中，GM为万有引力常数G与地球总质量M之积。然后根据广播星历中给出的摄动参数n计算观测时刻卫星的平均角速度n：nnn0。2.计算观测瞬间卫星的平近点角M:)(0TOEtnMM式中，0M为参考时刻TOE时的平近点角，由广播星历给出。3.计算偏近点角E：EeMEsin解上述方程可用迭代法或微分方程改正法。4.计算真近点角fEeEefEeeEfcos1sin1sincos1coscos2式中，e为卫星轨道的偏心率，由广播星历给出。5.计算升交距角u：fu式中，为近地点角距，由广播星历给出。6.计算摄动改正项iru,,：广播星历中给出了下列6个摄动参数：isicrsrcusucCCCCCC,,,,,，据此可以求出由于2J项而引起的升交距角u的摄动改正项u、卫星矢径r的摄动改正项u和卫星轨道倾角i的摄动改正项i。计算公式如下：uCuCuCuCuCuCisicirsrcrusucu2sin2cos2sin2cos2sin2cos7.计算0,,iru进行摄动改正)()cos1(0TOEtdtdiiiEarruuirru式中：a为卫星轨道的长半径，2)(Aa，0i为TOE时刻的轨道倾角，由广播星历中的开普勒六参数给出，dtdi为i的变化率，由广播星历中的摄动九参数给出。8.计算卫星在轨道面坐标系中的位置在轨道平面直角坐标系中（坐标原点位于地心，X轴指向升交点）卫星的平面直角坐标为：uryurxsincos9.计算观测瞬间升交点的经度L若参考时刻TOE时升交点的赤经为TOE，升交点对时间的变化率，那么观测瞬间t的升交点赤经应为：)(TOEtTOE可从广播星历的摄动参数中给出。设本周开始时刻（星期日0时）格林尼治恒星时为weekGAST，则观测瞬间的格林尼治恒星时为：tGASTGASTeweek式中：e为地球自转角速度；t为本周内的时间)(s。这样就可求得观测瞬间升交点的经度值为：tTOEtGASTGASTLeweekTOE)(令weekTOEGAST0，则有：tTOEtLe)(010.计算卫星在瞬时地球坐标系中的位置已知升交点的大地经度L以及轨道平面的倾角i后，就可通过两次旋转方便地求得卫星在地固坐标系中的位置：LyLiyLxLiyLxyxiRLRZYXXZsinsincoscossincoscos0)()(11.计算卫星在协议地球坐标系中的位置观测瞬间卫星在协议地球坐标系中的位置ZYXyxyxZYXyRxRzyxpppppXPYCTS11001)()(