从广播星历计算卫星位置.doc

从广播星历计算卫星位置：1. 计算卫星运动的平均角速度 n首先根据广播星历中给出的参数 计算出参考时刻 的平均角速度 ：ATOE0n，30)(GM式中， 为万有引力常数 与地球总质量 之积然后根据广播星历中给出的GM摄动参数 计算观测时刻卫星的平均角速度 ： nnn02. 计算观测瞬间卫星的平近点角 : )(0TOEt式中， 为参考时刻 时的平近点角，由广播星历给出0TE3. 计算偏近点角 ： eMsin解上述方程可用迭代法或微分方程改正法4. 计算真近点角 fEefcos1insinc2式中， 为卫星轨道的偏心率，由广播星历给出e5. 计算升交距角 ：u fu式中， 为近地点角距，由广播星历给出6. 计算摄动改正项 ： iru,广播星历中给出了下列 6 个摄动参数： ，据此可以求出iscrscuscCC,,,由于 项而引起的升交距角 的摄动改正项 、卫星矢径 的摄动改正项 和卫星轨道2J u倾角 的摄动改正项 计算公式如下：iiuCisici rrsuc2noi7. 计算 进行摄动改正0,iru)(cos10TOEtdiiarui r式中： 为卫星轨道的长半径， ， 为 时刻的轨道倾角，由广播a2A0i星历中的开普勒六参数给出， 为 的变化率，由广播星历中的摄动九参数给出。

dti8. 计算卫星在轨道面坐标系中的位置在轨道平面直角坐标系中（坐标原点位于地心，X 轴指向升交点）卫星的平面直角坐标为： uryxsinco9. 计算观测瞬间升交点的经度 L若参考时刻 时升交点的赤经为 ，升交点对时间的变化率 ，那么观测TOETOE&瞬间 的升交点赤经 应为：t )(tTE&可从广播星历的摄动参数中给出&设本周开始时刻（星期日 0 时）格林尼治恒星时为 ，则观测瞬间的格林wekGAST尼治恒星时为： tASTGewk式中： 为地球自转角速度； 为本周内的时间 et)(s这样就可求得观测瞬间升交点的经度值为： tTOEtASASTL ewekOE)(&令 ，则有：wekTOEG0 L010. 计算卫星在瞬时地球坐标系中的位置已知升交点的大地经度 以及轨道平面的倾角 后，就可通过两次旋转方便地求得i卫星在地固坐标系中的位置：  LyixyiRLZYXXsnco0)(11. 计算卫星在协议地球坐标系中的位置观测瞬间卫星在协议地球坐标系中的位置 ZYXyxZYXyRxzyx ppXPYCTS 10)()。