Satellite navigation solution

Last updated June 19, 2024

Satellite navigation solution for the receiver's position (geopositioning) involves an algorithm. In essence, a GNSS receiver measures the transmitting time of GNSS signals emitted from four or more GNSS satellites (giving the pseudorange) and these measurements are used to obtain its position (i.e., spatial coordinates) and reception time.

The solution illustrated

Essentially, the solution shown in orange, $\scriptstyle ({\hat {\boldsymbol {r}}}_{\text{rec}},\,{\hat {t}}_{\text{rec}})$ , is the intersection of light cones.
The posterior distribution of the solution is derived from the product of the distribution of propagating spherical surfaces. (See animation.)

Calculation steps

A global-navigation-satellite-system (GNSS) receiver measures the apparent transmitting time, $\displaystyle {\tilde {t}}_{i}$ , or "phase", of GNSS signals emitted from four or more GNSS satellites ( $\displaystyle i\;=\;1,\,2,\,3,\,4,\,..,\,n$ ), simultaneously.^[1]
GNSS satellites broadcast the messages of satellites' ephemeris, $\displaystyle {\boldsymbol {r}}_{i}(t)$ , and intrinsic clock bias (i.e., clock advance), $\displaystyle \delta t_{{\text{clock,sv}},i}(t)$ ^{[ clarification needed ]} as the functions of (atomic) standard time, e.g., GPST.^[2]
The transmitting time of GNSS satellite signals, $\displaystyle t_{i}$ , is thus derived from the non-closed-form equations $\displaystyle {\tilde {t}}_{i}\;=\;t_{i}\,+\,\delta t_{{\text{clock}},i}(t_{i})$ and $\displaystyle \delta t_{{\text{clock}},i}(t_{i})\;=\;\delta t_{{\text{clock,sv}},i}(t_{i})\,+\,\delta t_{{\text{orbit-relativ}},\,i}({\boldsymbol {r}}_{i},\,{\dot {\boldsymbol {r}}}_{i})$ , where $\displaystyle \delta t_{{\text{orbit-relativ}},i}({\boldsymbol {r}}_{i},\,{\dot {\boldsymbol {r}}}_{i})$ is the relativistic clock bias, periodically risen from the satellite's orbital eccentricity and Earth's gravity field.^[2] The satellite's position and velocity are determined by $\displaystyle t_{i}$ as follows: $\displaystyle {\boldsymbol {r}}_{i}\;=\;{\boldsymbol {r}}_{i}(t_{i})$ and $\displaystyle {\dot {\boldsymbol {r}}}_{i}\;=\;{\dot {\boldsymbol {r}}}_{i}(t_{i})$ .
In the field of GNSS, "geometric range", $\displaystyle r({\boldsymbol {r}}_{A},\,{\boldsymbol {r}}_{B})$ , is defined as straight range, or 3-dimensional distance,^[3] from $\displaystyle {\boldsymbol {r}}_{A}$ to $\displaystyle {\boldsymbol {r}}_{B}$ in inertial frame (e.g., ECI one), not in rotating frame.^[2]
The receiver's position, $\displaystyle {\boldsymbol {r}}_{\text{rec}}$ , and reception time, $\displaystyle t_{\text{rec}}$ , satisfy the light-cone equation of $\displaystyle r({\boldsymbol {r}}_{i},\,{\boldsymbol {r}}_{\text{rec}})/c\,+\,(t_{i}-t_{\text{rec}})\;=\;0$ in inertial frame, where $\displaystyle c$ is the speed of light. The signal time of flight from satellite to receiver is $\displaystyle -(t_{i}\,-\,t_{\text{rec}})$ .
The above is extended to the satellite-navigation positioning equation, $\displaystyle r({\boldsymbol {r}}_{i},\,{\boldsymbol {r}}_{\text{rec}})/c\,+\,(t_{i}\,-\,t_{\text{rec}})\,+\,\delta t_{{\text{atmos}},i}\,-\,\delta t_{{\text{meas-err}},i}\;=\;0$ , where $\displaystyle \delta t_{{\text{atmos}},i}$ is atmospheric delay (= ionospheric delay + tropospheric delay) along signal path and $\displaystyle \delta t_{{\text{meas-err}},i}$ is the measurement error.
The Gauss–Newton method can be used to solve the nonlinear least-squares problem for the solution: $\displaystyle ({\hat {\boldsymbol {r}}}_{\text{rec}},\,{\hat {t}}_{\text{rec}})\;=\;\arg \min \phi ({\boldsymbol {r}}_{\text{rec}},\,t_{\text{rec}})$ , where $\displaystyle \phi ({\boldsymbol {r}}_{\text{rec}},\,t_{\text{rec}})\;=\;\sum _{i=1}^{n}(\delta t_{{\text{meas-err}},i}/\sigma _{\delta t_{{\text{meas-err}},i}})^{2}$ . Note that $\displaystyle \delta t_{{\text{meas-err}},i}$ should be regarded as a function of $\displaystyle {\boldsymbol {r}}_{\text{rec}}$ and $\displaystyle t_{\text{rec}}$ .
The posterior distribution of $\displaystyle {\boldsymbol {r}}_{\text{rec}}$ and $\displaystyle t_{\text{rec}}$ is proportional to $\displaystyle \exp(-{\frac {1}{2}}\phi ({\boldsymbol {r}}_{\text{rec}},\,t_{\text{rec}}))$ , whose mode is $\displaystyle ({\hat {\boldsymbol {r}}}_{\text{rec}},\,{\hat {t}}_{\text{rec}})$ . Their inference is formalized as maximum a posteriori estimation.
The posterior distribution of $\displaystyle {\boldsymbol {r}}_{\text{rec}}$ is proportional to $\displaystyle \int _{-\infty }^{\infty }\exp(-{\frac {1}{2}}\phi ({\boldsymbol {r}}_{\text{rec}},\,t_{\text{rec}}))\,dt_{\text{rec}}$ .

The GPS case

For Global Positioning System (GPS),^[2] the non-closed-form equations in step 3 result in

\scriptstyle {\begin{cases}\scriptstyle \Delta t_{i}(t_{i},\,E_{i})\;\triangleq \;t_{i}\,+\,\delta t_{{\text{clock}},i}(t_{i},\,E_{i})\,-\,{\tilde {t}}_{i}\;=\;0,\\\scriptstyle \Delta M_{i}(t_{i},\,E_{i})\;\triangleq \;M_{i}(t_{i})\,-\,(E_{i}\,-\,e_{i}\sin E_{i})\;=\;0,\end{cases}}

in which $\scriptstyle E_{i}$ is the orbital eccentric anomaly of satellite $i$ , $\scriptstyle M_{i}$ is the mean anomaly, $\scriptstyle e_{i}$ is the eccentricity, and $\scriptstyle \delta t_{{\text{clock}},i}(t_{i},\,E_{i})\;=\;\delta t_{{\text{clock,sv}},i}(t_{i})\,+\,\delta t_{{\text{orbit-relativ}},i}(E_{i})$ .

The above can be solved by using the bivariate Newton–Raphson method on $\scriptstyle t_{i}$ and $\scriptstyle E_{i}$ . Two times of iteration will be necessary and sufficient in most cases. Its iterative update will be described by using the approximated inverse of Jacobian matrix as follows:

$\scriptstyle {\begin{pmatrix}t_{i}\\E_{i}\\\end{pmatrix}}\leftarrow {\begin{pmatrix}t_{i}\\E_{i}\\\end{pmatrix}}-{\begin{pmatrix}1&&0\\{\frac {{\dot {M}}_{i}(t_{i})}{1-e_{i}\cos E_{i}}}&&-{\frac {1}{1-e_{i}\cos E_{i}}}\\\end{pmatrix}}{\begin{pmatrix}\Delta t_{i}\\\Delta M_{i}\\\end{pmatrix}}$

Tropospheric delay should not be ignored, while the Global Positioning System (GPS) specification^[2] doesn't provide its detailed description.

The GLONASS case

The GLONASS ephemerides don't provide clock biases $\scriptstyle \delta t_{{\text{clock,sv}},i}(t)$ , but $\scriptstyle \delta t_{{\text{clock}},i}(t)$ .

Notes

In the field of GNSS, $\scriptstyle {\tilde {r}}_{i}\;=\;-c({\tilde {t}}_{i}\,-\,{\tilde {t}}_{\text{rec}})$ is called pseudorange, where $\scriptstyle {\tilde {t}}_{\text{rec}}$ is a provisional reception time of the receiver. $\scriptstyle \delta t_{\text{clock,rec}}\;=\;{\tilde {t}}_{\text{rec}}\,-\,t_{\text{rec}}$ is called receiver's clock bias (i.e., clock advance).^[1]
Standard GNSS receivers output $\scriptstyle {\tilde {r}}_{i}$ and $\scriptstyle {\tilde {t}}_{\text{rec}}$ per an observation epoch.
The temporal variation in the relativistic clock bias of satellite is linear if its orbit is circular (and thus its velocity is uniform in inertial frame).
The signal time of flight from satellite to receiver is expressed as $\scriptstyle -(t_{i}-t_{\text{rec}})\;=\;{\tilde {r}}_{i}/c\,+\,\delta t_{{\text{clock}},i}\,-\,\delta t_{\text{clock,rec}}$ , whose right side is round-off-error resistive during calculation.
The geometric range is calculated as $\scriptstyle r({\boldsymbol {r}}_{i},\,{\boldsymbol {r}}_{\text{rec}})\;=\;|\Omega _{\text{E}}(t_{i}\,-\,t_{\text{rec}}){\boldsymbol {r}}_{i,{\text{ECEF}}}\,-\,{\boldsymbol {r}}_{\text{rec,ECEF}}|$ , where the Earth-centred, Earth-fixed (ECEF) rotating frame (e.g., WGS84 or ITRF) is used in the right side and $\scriptstyle \Omega _{\text{E}}$ is the Earth rotating matrix with the argument of the signal transit time.^[2] The matrix can be factorized as $\scriptstyle \Omega _{\text{E}}(t_{i}\,-\,t_{\text{rec}})\;=\;\Omega _{\text{E}}(\delta t_{\text{clock,rec}})\Omega _{\text{E}}(-{\tilde {r}}_{i}/c\,-\,\delta t_{{\text{clock}},i})$ .
The line-of-sight unit vector of satellite observed at $\scriptstyle {\boldsymbol {r}}_{\text{rec,ECEF}}$ is described as: $\scriptstyle {\boldsymbol {e}}_{i,{\text{rec,ECEF}}}\;=\;-{\frac {\partial r({\boldsymbol {r}}_{i},\,{\boldsymbol {r}}_{\text{rec}})}{\partial {\boldsymbol {r}}_{\text{rec,ECEF}}}}$ .
The satellite-navigation positioning equation may be expressed by using the variables $\scriptstyle {\boldsymbol {r}}_{\text{rec,ECEF}}$ and $\scriptstyle \delta t_{\text{clock,rec}}$ .
The nonlinearity of the vertical dependency of tropospheric delay degrades the convergence efficiency in the Gauss–Newton iterations in step 7.
The above notation is different from that in the Wikipedia articles, 'Position calculation introduction' and 'Position calculation advanced', of Global Positioning System (GPS).

Related Research Articles

In celestial mechanics, an orbit is the curved trajectory of an object such as the trajectory of a planet around a star, or of a natural satellite around a planet, or of an artificial satellite around an object or position in space such as a planet, moon, asteroid, or Lagrange point. Normally, orbit refers to a regularly repeating trajectory, although it may also refer to a non-repeating trajectory. To a close approximation, planets and satellites follow elliptic orbits, with the center of mass being orbited at a focal point of the ellipse, as described by Kepler's laws of planetary motion.

In physics, spacetime is a mathematical model that fuses the three dimensions of space and the one dimension of time into a single four-dimensional continuum. Spacetime diagrams are useful in visualizing and understanding relativistic effects, such as how different observers perceive where and when events occur.

A gyrocompass is a type of non-magnetic compass which is based on a fast-spinning disc and the rotation of the Earth to find geographical direction automatically. A gyrocompass makes use of one of the seven fundamental ways to determine the heading of a vehicle. A gyroscope is an essential component of a gyrocompass, but they are different devices; a gyrocompass is built to use the effect of gyroscopic precession, which is a distinctive aspect of the general gyroscopic effect. Gyrocompasses, such as the fibre optic gyrocompass are widely used to provide a heading for navigation on ships. This is because they have two significant advantages over magnetic compasses:

<span class="mw-page-title-main">Angular velocity</span> Direction and rate of rotation

In physics, angular velocity, also known as angular frequency vector, is a pseudovector representation of how the angular position or orientation of an object changes with time, i.e. how quickly an object rotates around an axis of rotation and how fast the axis itself changes direction.

<span class="mw-page-title-main">Moment of inertia</span> Scalar measure of the rotational inertia with respect to a fixed axis of rotation

The moment of inertia, otherwise known as the mass moment of inertia, angular/rotational mass, second moment of mass, or most accurately, rotational inertia, of a rigid body is a quantity that determines the torque needed for a desired angular acceleration about a rotational axis, akin to how mass determines the force needed for a desired acceleration. It depends on the body's mass distribution and the axis chosen, with larger moments requiring more torque to change the body's rate of rotation by a given amount.

In physics, a Langevin equation is a stochastic differential equation describing how a system evolves when subjected to a combination of deterministic and fluctuating ("random") forces. The dependent variables in a Langevin equation typically are collective (macroscopic) variables changing only slowly in comparison to the other (microscopic) variables of the system. The fast (microscopic) variables are responsible for the stochastic nature of the Langevin equation. One application is to Brownian motion, which models the fluctuating motion of a small particle in a fluid.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space $under the operation of composition.$

A directional derivative is a concept in multivariable calculus that measures the rate at which a function changes in a particular direction at a given point.

In physics, the Rabi cycle is the cyclic behaviour of a two-level quantum system in the presence of an oscillatory driving field. A great variety of physical processes belonging to the areas of quantum computing, condensed matter, atomic and molecular physics, and nuclear and particle physics can be conveniently studied in terms of two-level quantum mechanical systems, and exhibit Rabi flopping when coupled to an optical driving field. The effect is important in quantum optics, magnetic resonance and quantum computing, and is named after Isidor Isaac Rabi.

Dilution of precision (DOP), or geometric dilution of precision (GDOP), is a term used in satellite navigation and geomatics engineering to specify the error propagation as a mathematical effect of navigation satellite geometry on positional measurement precision.

Spacecraft flight dynamics is the application of mechanical dynamics to model how the external forces acting on a space vehicle or spacecraft determine its flight path. These forces are primarily of three types: propulsive force provided by the vehicle's engines; gravitational force exerted by the Earth and other celestial bodies; and aerodynamic lift and drag.

In geometry and linear algebra, a Cartesian tensor uses an orthonormal basis to represent a tensor in a Euclidean space in the form of components. Converting a tensor's components from one such basis to another is done through an orthogonal transformation.

<span class="mw-page-title-main">Flow (mathematics)</span> Motion of particles in a fluid

In mathematics, a flow formalizes the idea of the motion of particles in a fluid. Flows are ubiquitous in science, including engineering and physics. The notion of flow is basic to the study of ordinary differential equations. Informally, a flow may be viewed as a continuous motion of points over time. More formally, a flow is a group action of the real numbers on a set.

<span class="mw-page-title-main">Rotation around a fixed axis</span> Type of motion

Rotation around a fixed axis or axial rotation is a special case of rotational motion around an axis of rotation fixed, stationary, or static in three-dimensional space. This type of motion excludes the possibility of the instantaneous axis of rotation changing its orientation and cannot describe such phenomena as wobbling or precession. According to Euler's rotation theorem, simultaneous rotation along a number of stationary axes at the same time is impossible; if two rotations are forced at the same time, a new axis of rotation will result.

<span class="mw-page-title-main">Reaction–diffusion system</span> Type of mathematical model

Reaction–diffusion systems are mathematical models that correspond to several physical phenomena. The most common is the change in space and time of the concentration of one or more chemical substances: local chemical reactions in which the substances are transformed into each other, and diffusion which causes the substances to spread out over a surface in space.

In applied mathematics, discontinuous Galerkin methods (DG methods) form a class of numerical methods for solving differential equations. They combine features of the finite element and the finite volume framework and have been successfully applied to hyperbolic, elliptic, parabolic and mixed form problems arising from a wide range of applications. DG methods have in particular received considerable interest for problems with a dominant first-order part, e.g. in electrodynamics, fluid mechanics and plasma physics. Indeed, the solutions of such problems may involve strong gradients (and even discontinuities) so that classical finite element methods fail, while finite volume methods are restricted to low order approximations.

GNSS enhancement refers to techniques used to improve the accuracy of positioning information provided by the Global Positioning System or other global navigation satellite systems in general, a network of satellites used for navigation. Enhancement methods of improving accuracy rely on external information being integrated into the calculation process. There are many such systems in place and they are generally named or described based on how the GPS sensor receives the information. Some systems transmit additional information about sources of error, others provide direct measurements of how much the signal was off in the past, while a third group provides additional navigational or vehicle information to be integrated into the calculation process.

Orbit modeling is the process of creating mathematical models to simulate motion of a massive body as it moves in orbit around another massive body due to gravity. Other forces such as gravitational attraction from tertiary bodies, air resistance, solar pressure, or thrust from a propulsion system are typically modeled as secondary effects. Directly modeling an orbit can push the limits of machine precision due to the need to model small perturbations to very large orbits. Because of this, perturbation methods are often used to model the orbit in order to achieve better accuracy.

Beamforming is a signal processing technique used to spatially select propagating waves. In order to implement beamforming on digital hardware the received signals need to be discretized. This introduces quantization error, perturbing the array pattern. For this reason, the sample rate must be generally much greater than the Nyquist rate.

References

1 2 Misra, P. and Enge, P., Global Positioning System: Signals, Measurements, and Performance, 2nd, Ganga-Jamuna Press, 2006.
1 2 3 4 5 6 The interface specification of NAVSTAR GLOBAL POSITIONING SYSTEM
↑ 3-dimensional distance is given by $\displaystyle r({\boldsymbol {r}}_{A},\,{\boldsymbol {r}}_{B})=|{\boldsymbol {r}}_{A}-{\boldsymbol {r}}_{B}|={\sqrt {(x_{A}-x_{B})^{2}+(y_{A}-y_{B})^{2}+(z_{A}-z_{B})^{2}}}$ where $\displaystyle {\boldsymbol {r}}_{A}=(x_{A},y_{A},z_{A})$ and $\displaystyle {\boldsymbol {r}}_{B}=(x_{B},y_{B},z_{B})$ represented in inertial frame.

External links

PVT (Position, Velocity, Time): Calculation procedure in the open-source GNSS-SDR and the underlying RTKLIB

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[Misra_Enge-1] 1 2 Misra, P. and Enge, P., Global Positioning System: Signals, Measurements, and Performance, 2nd, Ganga-Jamuna Press, 2006.

[IS-GPS-2] 1 2 3 4 5 6 The interface specification of NAVSTAR GLOBAL POSITIONING SYSTEM

[3] 3-dimensional distance is given by $\displaystyle r({\boldsymbol {r}}_{A},\,{\boldsymbol {r}}_{B})=|{\boldsymbol {r}}_{A}-{\boldsymbol {r}}_{B}|={\sqrt {(x_{A}-x_{B})^{2}+(y_{A}-y_{B})^{2}+(z_{A}-z_{B})^{2}}}$ where $\displaystyle {\boldsymbol {r}}_{A}=(x_{A},y_{A},z_{A})$ and $\displaystyle {\boldsymbol {r}}_{B}=(x_{B},y_{B},z_{B})$ represented in inertial frame.

[1]

[2]

[3]