Position estimation enhancement for a global navigation satellite system receiver

• US 8,508,406 B2
• Filed: 03/20/2008
• Issued: 08/13/2013
• Est. Priority Date: 03/20/2008
• Status: Active Grant

• Current Assignee: Telespazio SPA (Leonardo S.p.A.)
• Original Assignee: Telespazio SPA (Leonardo S.p.A.)
• Inventors: Rodriguez, Filippo, Plaia, Gian Paolo
• Primary Examiner(s): Keith, Jack W
• Assistant Examiner(s): Mull, Fred H
• Application Number:
  US12/933,279
• Publication Number:
  US 20110090117A1
• Time in Patent Office: 1,972 Days
• Field of Search: 342/357.61, 342/357.59
• US Class Current: 342/357.61
• CPC Class Codes: G01S 19/22 Multipath-related issues
• CPC Class Codes: G01S 19/428 using multipath or indirect path propagation signals in position determination

Abstract

A method of estimating a position of a satellite receiver, comprising computing a weight matrix, and computing an estimated position of the satellite receiver based on the weight matrix, wherein computing the weight matrix includes computing quantities indicative of degradations experienced by satellite signals and of multipath interference, and computing the weight matrix based on the computed quantities.

Claims

1. Method of estimating a position of a satellite receiver, comprising:

computing a weight matrix; and

computing, by a programmed processor, an estimated position of the satellite receiver based on the weight matrix;

characterized in that computing the weight matrix includes;

computing quantities indicative of degradations experienced by satellite signals and due to multipath interference;

computing the weight matrix based on the computed quantities; and

wherein computing the weight matrix includes;

classifying received signals into elevation classes based on elevations of the satellites that have emitted the signals;

for each elevation class, classifying the signals belonging to the elevation class into degradation classes based on the corresponding degradations, each degradation class having a corresponding cardinality defined by the number of signals with degradations falling within the degradation class;

for each elevation class, identifying different distribution functions which globally fit the distribution of the corresponding numbers of signal with degradations falling within the corresponding degradation classes;

for each elevation class, associating each degradation class with one of the corresponding distribution functions that satisfies a given criterion in relation to the number of signals with degradations falling within the corresponding degradation class; and

computing the weight matrix based on quantities related to parameters of the distribution functions associated with the degradation classes.

2. The method of claim 1, wherein the weight matrix includes weight entries each associated with a corresponding satellite, and wherein computing a quantity indicative of a degradation experienced by satellite signals from a corresponding satellite and due to multipath interference includes:

measuring a quantity indicative of an attenuation experienced by a satellite signal in a substantially multipath-free environment;

measuring a quantity indicative of an attenuation of a satellite signal in a multipath-affected environment; and

computing the quantity indicative of the degradation based on the measured quantities.

3. The method of claim 1, further including:

for each elevation class, identifying groups of consecutive degradation classes associated with the same distribution function.

4. The method of claim 1, wherein the weight matrix includes diagonal entries each associated with a corresponding satellite, and wherein computing the weight matrix based on quantities related to parameters of the distribution functions associated with the degradation classes includes:

computing each diagonal entry of the weight matrix based on a quantity related to a parameter of the distribution function associated with the degradation class containing the degradations of signals from the corresponding satellite, and the elevation class containing the elevation of the corresponding satellite.

5. The method of claim 1, wherein the parameter of the distribution function is the variance of the distribution function, and the quantity related to the parameter is the reciprocal of the variance.

6. The method of claim 1, wherein the different distribution functions include Rice, Rayleigh, and Loo distribution functions.

7. The method of claim 1, wherein the quantities indicative of degradations due to multipath interference are computed based on satellite signals emitted in a given time span.

8. The method of claim 1, wherein computing an estimated position of the satellite receiver based on the weight matrix includes:

computing iteratively the following equation system;

{circumflex over (X)}=({tilde over (G)}^T·

W·

{tilde over (G)})<sup>−

1</sup>·

{tilde over (G^T)}·

W·

{tilde over (y)}

wherein;

{circumflex over (X)} is a vector indicative of the estimated satellite receiver position;

W is the weight matrix;

{tilde over (G)} is a matrix with a number of rows equal to the number N of satellites the signals emitted thereby are received by the satellite receiver, each row being associated with a corresponding satellite and being in the form [cos El_i cos Az₁ cos El_i cos Az_i sin El_i 1], wherein El_i and Az_i are, respectively, the elevation and the azimuth of the corresponding satellite;

{tilde over (G)}^T is the transpose of the matrix {tilde over (G)}; and

{tilde over (y)} is a vector with N entries each associated with a corresponding satellite the signals emitted thereby are received by the satellite receiver, each entry being equal to the difference between a corresponding computed satellite pseudorange and a computed geometrical distance between the satellite receiver and the corresponding satellite.

9. A satellite navigation system including a satellite constellation, a ground local element configured to communicate with the satellites, and a satellite receiver configured to communicate with the satellites and the ground local element;

characterized by a system for estimating a position of the satellite receiver configured to implement the method according to claim 1.

10. The system of claim 9, wherein the system for estimating the position of the satellite receiver is comprised in the satellite receiver.

11. The system of claim 9, wherein the system for estimating the position of the satellite receiver is comprised in the ground local element.