Testing of the Conditioning of the GPS

Robert Argus, Aneesh Malhotra, Kathleen McLane

Satellites Spaced Apart

The fourth and fifth part of the project involved testing the conditioning system of the GPS problem. The three satellite positions were put into spherical coordinates such that \(A_i,B_i,C_i\) were now \begin{align*} A_i = \rho cos(\phi_i) cos(\theta_i)\\ B_i = \rho cos(\phi_i) sin(\theta_i)\\ C_i = \rho sin(\phi_i) \end{align*} where \(\rho\) = 26570 km, and \(\phi_i\) was four different values between 0 \(\leq \phi_i \leq \pi/2\) and \(\theta_i\) was four different values between 0 \(\leq \theta_i \leq 2\pi \). The four values used for \(\theta\) and \(\phi\) can be seen in the table below:

\(\phi \) Values \(\theta\) Values
\(\pi/2\) 0
\(\pi/2\) + 0.0785 0.0785
\(\pi/2\) -0.0785 0.0393
\(\pi/2\) +0.0393 0.05
It was then chosen that x=0, y=0, z=6370, and d=0.0001. These were then used to calculate the satellite ranges and travel times with the formulas, respectively: \begin{align*} R_i = \sqrt{A_i^2+B_i^2+(C_i+6370)^2}\\ t_i = d + R_i/c \end{align*}

Using this, we were able to define an error magnification factor for this situation using the formula provided in the book for Reality Check 4. It was stated that the atomic clocks aboard the satellites are correct up to 10 nanoseconds, so \(\Delta t_i = 10^{-8}\).
This allows the EMF to be defined as \begin{align*} emf = \frac{\|\Delta x, \Delta y, \Delta z\|_\inf}{c\|\Delta t_1,...,\Delta t_m\|_\inf} \end{align*}

From this, the \( t_i\) were changed between \(+10^{-8}\) and \(-10^{-8}\), and the error magnification factor was calculated. The codes used are H.m, and DH.m. The results can be seen in the table below.

EMF Maximum Position Error(km)
5.03009811062191 0.015079854765645
5.03009978100344 0.0150798597733228
1.15891466280988 0.00347433875376015
1.1589146589197 0.0034743387420977
1.42101635169679 0.00426009984933372
1.42101626925138 0.00426009960216862

Based on the table, it was seen that the max position error was 0.0150798597733228 kilometers, and the condition number was 5.03009978100344.

Satellites Tightly Bunched

The fifth part of the project was similar to the fourth, except the satellites were spaced within 5% of each other. This meant that all the \(\phi_i\) and \(\theta_i\) were all spaced within 5% of each other. The codes for this part of the project can be seen here I.m and DI.m. The error magnification factor was calculated for each of them and the results can be seen in the table below.

EMF Maximum Position Error(km)
433.097162600829 1.29839262928928
433.138427463415 1.29851633823512
957.010271932992 2.8690446175404
957.214585332525 2.86965713370289
714.193988739568 2.14109971373059
714.078079878503 2.14075222770697

The maximum position error was 2.86965713370289, with a condition number of 957.214585332525. It was found that the EMF decreased as the satellites moved further apart from each other. Thus it was better if the satellites are unbunched.