GPS and Galileo Satellite Coordinates Computation - Navipedia

Download as pdf or txt
Download as pdf or txt
You are on page 1of 3

GPS and Galileo Satellite Coordinates Computation - Navipedia https://gssc.esa.int/navipedia/index.php/GPS_and_Galileo_Satellite_Co...

GPS and Galileo Satellite Coordinates


Computation

Table 1 provides the GPS or Galileo broadcast ephemeris Fundamentals


parameters to compute their satellite coordinates at any GPS and Galileo Satellite
Title
observation epoch. These parameters are periodically renewed Coordinates Computation
(typically every 2 hours for GPS and 3 hours for Galileo) and J. Sanz Subirana, J.M. Juan
Zornoza and M. Hernández-
must not be used out of the prescribed time (about four Author(s)
Pajares, Technical University of
hours), because the extrapolation error grows exponentially Catalonia, Spain.
beyond its validity period.
Level Intermediate

Year of
2011
The algorithm provided is from the [IS-GPS-200, table 20-IV] Publication
[footnotes 1].The Galileo satellites follow a similar scheme

Table 1: GPS and Galileo broadcast ephemeris and clock message parameters.

In order to compute satellite coordinates from navigation message, the algorithm provided as follows must be used. An

1 of 3 6/9/2023, 11:14 PM
GPS and Galileo Satellite Coordinates Computation - Navipedia https://gssc.esa.int/navipedia/index.php/GPS_and_Galileo_Satellite_Co...

accuracy of about 4 meters (RMS) is achieved for GPS satellites with S/A=0ff and several tens of meters with S/A=on
[footnotes 2]:

▪ Compute the time tk from the ephemerides reference epoch toe (t and toe are expressed in seconds in the GPS
week):

tk = t − toe
> 302 400 sec, subtract 604 800 sec from tk . If tk < −302 400 sec, add
If tk
604 800 sec.

▪ Compute the mean anomaly for tk ,

Mk = Mo + ( √
√μ
+ Δn) tk
a3

▪ Solve (iteratively) the Kepler equation for the eccentricity anomaly Ek :

Mk = Ek − e sin Ek

▪ Compute the true anomaly vk :

)
√1−e2 sin Ek
vk = arctan( cos Ek −e

▪ Compute the argument of latitude uk from the argument of perigee ω , true anomaly vk and corrections cuc and
cus :

uk = ω + vk + cuc cos 2 (ω + vk ) + cus sin 2 (ω + vk )

▪ Compute the radial distance rk , considering corrections crc and crs :

rk = a (1 − e cos Ek ) + crc cos 2 (ω + vk ) + crs sin 2 (ω + vk )

▪ Compute the inclination ik of the orbital plane from the inclination io at reference time toe , and corrections cic
and cis :

2 of 3 6/9/2023, 11:14 PM
GPS and Galileo Satellite Coordinates Computation - Navipedia https://gssc.esa.int/navipedia/index.php/GPS_and_Galileo_Satellite_Co...


ik = io + i tk + cic cos 2 (ω + vk ) + cis sin 2 (ω + vk )

▪ Compute the longitude of the ascending node λk (with respect to Greenwich). This calculation uses the right
ascension at the beginning of the current week (Ωo ), the correction from the apparent sidereal time variation in
Greenwich between the beginning of the week and reference time tk = t − toe, and the change in longitude of
the ascending node from the reference time toe :


λk = Ωo + (Ω −ωE ) tk − ωE toe

▪ Compute the coordinates in TRS frame, applying three rotations (around uk , ik and λk ):

⎡ Xk ⎤ ⎡ rk ⎤
⎢ Yk ⎥ = R3 (−λk ) R1 (−ik ) R3 (−uk ) ⎢ 0 ⎥
⎣Z ⎦ ⎣0⎦
k

where R1 and R3 are the rotation matrices defined in Transformation between Terrestrial
Frames.

Notes
1. ^ [IS-GPS-200], NAVSTAR GPS Space Segment/Navigation User Interfaces https://www.gps.gov/technical
/icwg/IS-GPS-200M.pdf
2. ^ Actually, the S/A was mainly applied to the satellite clocks and, apparently, not so often to the ephemeris.

Retrieved from "https://gssc.esa.int/navipedia/index.php?title=GPS_and_Galileo_Satellite_Coordinates_Computation&


oldid=15918"

This page was last edited on 10 September 2021, at 10:51.

3 of 3 6/9/2023, 11:14 PM

You might also like