Analysis and improvement of the Bancroft algorithm for GNSS satellite orbit determination

Abstract

Satellite orbit information is crucial for ensuring that global navigation satellite systems (GNSSs) provide appropriate positioning, navigation and timing services. Typically, users can obtain access to orbit information of a specific accuracy level from navigation messages or precise ephemeris products. Without this information, a system will not be able to provide normal service. In response to this problem, initial orbit information of a certain level of precision must be obtained to support subsequent applications, such as broadcasting or precise ephemeris calculations, thereby ensuring the successful subsequent operation of the navigation system. One of two ways to calculate the initial orbit of a GNSS satellite is to utilize ground tracking stations to observe satellite vector information in the geocentric inertial system; the second way is to utilize GNSS range observations and known orbit information from other satellites. For the second approach, some researchers use the Bancroft algorithm combined with receiver clock offset to determine the initial orbit of GNSS satellites. Because this method requires an additional known receiver clock offset, we study the dependence of the Bancroft algorithm on clock offset in GNSS orbit determination. By assessing the impact of errors of different magnitude on the accuracy of the orbit results, we obtain experimental conclusions. After comprehensively analyzing various errors, we determine the accuracy level that the Bancroft algorithm can achieve for orbit determination without considering receiver clock correction. Dual-frequency and single-frequency pseudorange data from International GNSS Service stations are used in orbit determination experiments. When a small receiver clock offset is considered and no correction is made, the deviations in the calculated satellite positions in three dimensions are approximately 979.3 and 1118.1 m (dual and single frequency); with a satellite clock offset, these values are approximately 928.8 and 1062.7 m (dual and single frequency).