New Applications of Elliptic Functions and Integrals in GPS Inter-Satellite Communications with Account of General Relativity Theory

Abstract

During the last 15–20 years, the experimental methods for autonomous navigation and inter-satellite links have been developing rapidly in order to ensure navigation control and data processing without commands from Earth stations. Inter-satellite links are related to relative ranging between the satellites from one constellation or different constellations and measuring the distances between them with the precision of at least 1 $\mathsf{\mu }$m micrometer (=${10}^{-6}$ m), which should account for the bending of the light (radio or laser) signals due to the action of the Earth’s gravitational field. Thus, the theoretical calculation of the propagation time of a signal should be described in the framework of general relativity theory and the s.c. null cone equation. This review paper summarizes the latest achievements in calculating the propagation time of a signal, emitted by a GPS satellite, moving along a plane elliptical orbit or a space-oriented orbit, described by the full set of six Kepler parameters. It has been proved that for the case of plane elliptical orbit, the propagation time is expressed by a sum of elliptic integrals of the first, the second and the third kind, while for the second case (assuming that only the true anomaly angle is the dynamical parameter), the propagation time is expressed by a sum of elliptic integrals of the second and of the fourth order. For both cases, it has been proved that the propagation time represents a real-valued expression and not an imaginary one, as it should be. For the typical parameters of a GPS orbit, numerical calculations for the first case give acceptable values of the propagation time and, especially, the Shapiro delay term of the order of nanoseconds, thus confirming that this is a propagation time for the signal and not for the time of motion of the satellite. Theoretical arguments, related to general relativity and differential geometry have also been presented in favor of this conclusion. A new analytical method has been developed for transforming an elliptic integral in the Legendre form into an integral in the Weierstrass form. Two different representations have been found, one of them based on the method of four-dimensional uniformization, exposed in the monograph of Whittaker and Watson. The result of this approach is a new formulae for the Weierstrass invariants, depending in a complicated manner on the modulus parameter q of the elliptic integral in the Legendre form.