A satellite ephemeris correction calculating algorithm based on the Hill differential equation

The Hill differential equation is used to describe the change characteristic of the ephemeris prediction error in order to calculate the satellite broadcast ephemeris error and the error change rate. After the ECEF coordinate system converses the satellite inertial coordinate system, both the satellite ephemeris correction ΔR^j_S and the correction accuracy P^j_S,ocan be calculated by the Kalman filter. When the GEO, IGSO and MEO satellites move over the China, the calculation algorithms based on the weighted least square, the velocity model of the Kalman filter and the Hill differential equation are compared in the computing performance of the ephemeris corrections by using the reference date of IGS website. The simulation results show that: the ephemeris correction solver value based on the Hill differential equation can accurately reflect the ephemeris error characteristics at the current time. And it is able to correctly reflect the perturbation impact for the satellite ephemeris error. The Hill differential equation method is more suitable at the master station for engineering application about the MEO, IGSO and GEO satellites.