PART41 » History » Version 16

COLIN, Tony, 03/19/2016 09:22 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 16 COLIN, Tony
These tables are extracted from *[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 7 COLIN, Tony
h3. a - Reminder about the 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 1 COLIN, Tony
*Figure 4.1 :* Pseudo-range measurement extracted from *[2]*
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 7 COLIN, Tony
or in matrix equation form,
57 7 COLIN, Tony
58 7 COLIN, Tony
p=. !Pos5.png!
59 7 COLIN, Tony
*Equation 5*
60 7 COLIN, Tony
61 7 COLIN, Tony
which can be expressed as :
62 7 COLIN, Tony
63 7 COLIN, Tony
p=. !Pos6.png!
64 7 COLIN, Tony
*Equation 6*
65 7 COLIN, Tony
66 7 COLIN, Tony
with the Least Square solution :
67 7 COLIN, Tony
68 7 COLIN, Tony
p=. !Pos7.png!
69 7 COLIN, Tony
*Equation 7*
70 7 COLIN, Tony
71 8 COLIN, Tony
Thus, it is possible to retrieve the receiver position.
72 8 COLIN, Tony
73 8 COLIN, Tony
_Note that all unknowns are depicted in red color._
74 5 COLIN, Tony
75 5 COLIN, Tony
---
76 5 COLIN, Tony
77 5 COLIN, Tony
*References :* 
78 5 COLIN, Tony
*[1]* K. Borre, D. M. Akos, N. Bertelsen, P. Rinder, S. H. Jensen, A software-defined GPS and GALILEO receiver
79 5 COLIN, Tony
*[2]* Position Estimation Workshop, March 2016