PART41 » History » Version 22

COLIN, Tony, 03/19/2016 09:32 PM

1 3 COLIN, Tony
h1. PART 4 : Position Estimation.
2 2 COLIN, Tony
3 2 COLIN, Tony
{{toc}}
4 2 COLIN, Tony
5 13 COLIN, Tony
p(. Once the navigation bits from at least 4 satellites have been retrieved from the acquisition/tracking part, it is possible to estimate the desired position of the receiver.
6 2 COLIN, Tony
7 2 COLIN, Tony
---
8 2 COLIN, Tony
9 7 COLIN, Tony
h2. 1 - Ephemeris.
10 2 COLIN, Tony
11 16 COLIN, Tony
GPS uses a particular algorithm in order to characterise satellite position. In comparison with GLONASS, this method requires more parameters, but less complexity.
12 16 COLIN, Tony
13 14 COLIN, Tony
h3. a - GPS satellite ephemeris data.
14 9 COLIN, Tony
15 9 COLIN, Tony
p=. !Eph12min.png!
16 9 COLIN, Tony
17 14 COLIN, Tony
h3. b - GPS satellite position calculation algorithm.
18 9 COLIN, Tony
19 1 COLIN, Tony
p=. !Alg12min.png!
20 1 COLIN, Tony
21 18 COLIN, Tony
These tables are extracted from GPS Interface Control Document *[2]*
22 11 COLIN, Tony
23 7 COLIN, Tony
---
24 1 COLIN, Tony
25 7 COLIN, Tony
h2. 2 - Navigation computation.
26 7 COLIN, Tony
27 19 COLIN, Tony
h3. a - Reminder about the range impairments.
28 7 COLIN, Tony
29 4 COLIN, Tony
The following figure gives the impairments affecting the range in case of the GPS system as well as the correction process :
30 4 COLIN, Tony
31 4 COLIN, Tony
p=. !003.PNG!
32 20 COLIN, Tony
*Figure 4.1 :* Pseudo-range measurement extracted from *[3]*
33 7 COLIN, Tony
34 7 COLIN, Tony
h3. b - Demonstration of the Pseudo-ranges with Least Square method.
35 7 COLIN, Tony
36 7 COLIN, Tony
Starting from the fact that can determine most of the elements within the pseudo-range measurement PR_sat(i) from the information provided by each satellite, we have the equation : 
37 7 COLIN, Tony
38 7 COLIN, Tony
p=. !Pos1.png!
39 7 COLIN, Tony
*Equation 1*
40 7 COLIN, Tony
41 7 COLIN, Tony
or put in another way,
42 7 COLIN, Tony
43 7 COLIN, Tony
p=. !Pos2.png!
44 7 COLIN, Tony
*Equation 2*
45 7 COLIN, Tony
46 7 COLIN, Tony
Indeed 4 measurements are needed, providing 4 equations with 4 unknows which are the receiver coordinates and the clock bias of the receiver. As the equation is highly non-linear, it is important to proceed to a linearization such as the Taylor expansion :
47 7 COLIN, Tony
48 7 COLIN, Tony
p=. !Pos3.png!
49 7 COLIN, Tony
*Equation 3*
50 7 COLIN, Tony
51 7 COLIN, Tony
Hence,
52 7 COLIN, Tony
53 7 COLIN, Tony
p=. !Pos4.png!
54 7 COLIN, Tony
*Equation 4*
55 7 COLIN, Tony
56 22 COLIN, Tony
In practise, for a receiver located e.g. in France PR(t_0) can be described by Paris location as initialization for the algorithm.
57 21 COLIN, Tony
In vectorial form the equation becomes :
58 7 COLIN, Tony
59 7 COLIN, Tony
p=. !Pos5.png!
60 7 COLIN, Tony
*Equation 5*
61 7 COLIN, Tony
62 7 COLIN, Tony
which can be expressed as :
63 7 COLIN, Tony
64 7 COLIN, Tony
p=. !Pos6.png!
65 7 COLIN, Tony
*Equation 6*
66 7 COLIN, Tony
67 7 COLIN, Tony
with the Least Square solution :
68 7 COLIN, Tony
69 7 COLIN, Tony
p=. !Pos7.png!
70 7 COLIN, Tony
*Equation 7*
71 7 COLIN, Tony
72 8 COLIN, Tony
Thus, it is possible to retrieve the receiver position.
73 8 COLIN, Tony
74 8 COLIN, Tony
_Note that all unknowns are depicted in red color._
75 5 COLIN, Tony
76 5 COLIN, Tony
---
77 5 COLIN, Tony
78 5 COLIN, Tony
*References :* 
79 5 COLIN, Tony
*[1]* K. Borre, D. M. Akos, N. Bertelsen, P. Rinder, S. H. Jensen, A software-defined GPS and GALILEO receiver
80 17 COLIN, Tony
*[2]* GPS Interface Control Document under http://www.gps.gov/technical/icwg/IS-GPS-200H.pdf
81 17 COLIN, Tony
*[3]* Position Estimation Workshop, March 2016