Method of estimating a line of sight signal propagation time using a reduced-multipath correlation function

• US 5,615,232 A
• Filed: 11/24/1993
• Issued: 03/25/1997
• Est. Priority Date: 11/24/1993
• Status: Expired due to Term

• Current Assignee: Novatel Incorporated (Hexagon AB), Technische Universiteit Delft
• Original Assignee: Novatel Communications Incorporated
• Inventors: Van Nee, D. J. R.
• Primary Examiner(s): Chin, Stephen
• Assistant Examiner(s): GLUCK, JEFFREY WILLIAM
• Application Number:
  US08/157,476
• Time in Patent Office: 1,217 Days
• Field of Search: 375/94, 375/99, 375/102, 375/340, 375/346, 375/349, 375/208, 375/209, 375/210, 455/52.3, 455/65
• US Class Current: 375/346
• CPC Class Codes: G01S 19/22 Multipath-related issues
• CPC Class Codes: H04B 1/709 Correlator structure
• CPC Class Codes: H04B 1/711 the interference being multi-path interference
• CPC Class Codes: H04B 1/712 Weighting of fingers for combining, e.g. amplitude control or phase rotation using an inner loop
• CPC Class Codes: H04B 2201/70715 with application-specific features
• CPC Class Codes: H04B 7/005 Control of transmission; Equalising

Abstract

A method and device for processing a signal are described, wherein an estimate of a multipath-induced contribution to a demodulated navigation signal is calculated and subtracted from said demodulated navigation signal to obtain an estimated line of sight contribution to said demodulated navigation signal, and a propagation time τ₀ is calculated from the thus calculated line of sight contribution to said demodulated navigation signal, such that a very accurate propagation time τ₀ of said navigation signal can be calculated.

Claims

1. In a signal processing system comprising a receiver, a method of processing a data carrying carrier modulated pseudo random noise signal received from a transmitter comprising the steps of:

receiving the data carrying signal;

providing a reference signal;

pertaining correlations on said received data carrying signal with said reference signal at a plurality of points along the correlation function R_x (τ

) of said received data carrying signal to provide a measure of the shape of the correlation function R_x (τ

);

estimating the multipath-induced contribution to the correlation function;

subtracting said estimated multipath-induced contribution from the correlation function to obtain an estimated line of sight contribution to the correlation function; and

calculating a propagation time τ

₀ for said signal propagating between said transmitter and said receiver, said propagation time τ

₀ being calculated from said estimated line of sight contribution to the correlation function.

2. Method according to claim 1, wherein said estimate calculating step further comprises the step of calculating parameters a_i, τ

_i and τ

_i as follows;

##EQU32## wherein ##EQU33## wherein R(τ

) is the reference correlation function;

wherein a_i, τ

_i and θ

_i are estimated values for a_i, τ

_i and θ

_i, respectively.

3. Method according to claim 2, wherein a reference correlation function is measured in the absence of noise and multipath, using a signal simulator, and wherein the shape of the measured reference correlation function is stored in a memory of the receiver, for instance in the form of a table or a function representing said shape.

4. Method according to claim 1, wherein said estimate calculating step further comprises the step of calculating parameters a_i, τ

_i and θ

_i as follows;

##EQU34## wherein ##EQU35## wherein R(τ

) is the reference correlation function;

a_i, τ

_i and θ

_i are estimated values for a_i, τ

_i and θ

_i, respectively.

5. Method according to claim 1, wherein said estimate calculating step further comprises performing the following steps:

a) calculating parameters τ

₀, a₀, θ

₀ as optimum values for the largest peak in R_x (τ

);

b) subtracting the contribution to the correlation function of said calculated parameters τ

₀, a₀, θ

₀ from R_x (τ

) to yield R₁ (τ

), according to the formula R₁ (τ

)=R_x (τ

) -a₀ ·

R(τ

-τ

₀) ·

exp(jθ

₀);

c) calculating parameters τ

₁, a₁, θ

₁ as optimum values for the largest peak in R₁ (τ

);

d) subtracting the contribution of said calculated parameters τ

₁, a₁, θ

₁ from R_x (τ

) to yield a new approximation of R₀ (τ

) according to the formula R₀ (τ

)=R_x (τ

)-a₁ -R(τ

-τ

₁) ·

exp(jθ

₁);

e) calculating new parameters τ

hd 0, a₀, θ

₀ as optimum values for the largest peak in the new approximation R₀ (τ

) as calculated in step (d); and

f) repeating the steps (b) , (c) , (d) , (e) until a suitable stop criterion is fulfilled.

6. Method according to claim 5, wherein the steps (b), (c), (d), (e) are repeated for a predetermined number of times.

7. Method according to claim 5, wherein the steps (b), (c), (d), (e) are repeated until at least one of the parameters τ

₀, τ

₀, θ

₀ and τ

₁, a₁, θ

₁ have converged according to predetermined convergence criteria.

8. Method according to claim 5, wherein, in order to calculate an estimated contribution for a third signal, the following steps are performed:

g) subtracting the contribution of the calculated parameters τ

₀, a₀, θ

₀ and τ

₁, a₁, θ

₁ from R_x (τ

) to yield R₂ (τ

), according to the formula R₂ (τ

)=R_x (τ

)-a₀ ·

R(τ

-τ

₀)·

exp (jθ

₀)-a_i ·

R(τ

-τ

₁) ·

exp(jθ

₁);

h) calculating parameters τ

₂, a₂, θ

₂ as optimum values for the largest peak in R₂ (τ

);

i) calculating new parameters τ

₀, τ

₀, θ

₀ as optimum values for the largest peak in R₁ (τ

)=R_x (τ

) -a₀ ·

R(τ

-τ

₀) ·

exp(jθ

₀)-a₂ ·

R (τ

-τ

₂)·

exp (jθ

₂), using the parameters τ

₀, a₀, θ

₀ as calculated in step (i);

k) calculating new parameters τ

₂, a₂, θ

₂ as optimum values for the largest peak in R₂ (τ

)=R_x (τ

)-a₀ ·

R(τ

-τ

₀) ·

exp(jθ

₀)-a₁ ·

R (τ

-τ

₁)·

exp(jθ

₁), using the parameters τ

₀, a₀, θ

₀ as calculated in step (i) and the parameters τ

₁, a₁, θ

₁ as calculated in step (j); and

     1) repeating steps (i), (j), (k), until a stop criterion is fulfilled.

9. Method according to claim 1, further comprising the step of estimating a phase θ

of a correlation peak using at least one sample R_i (τ

_max) of said peak where each sample is taken preferably but not necessarily at or close to a top of the peak of the correlation function R_x (τ

), wherein;

##EQU36##

10. Method according to claim 1 further comprising the step of calculating an inverse, function f^-1 (τ

_x) in advance, and storing said calculated inverse function in a memory of the receiver to facilitate calculating an estimated delay value τ

.

11. Method according to claim 10, wherein the function f(τ

_x) is defined as ##EQU37## wherein Δ

τ

is a sampling time interval.

12. Method according to claim 10, wherein the function f(τ

_x) is defined as ##EQU38## wherein Δ

τ

is a sampling interval, and τ

_c is the delay lying in the middle of two samples taken along the reference correlation function at times t₁, and t₂, respectively, and t₁ -{τ

₀ }τ

_c -1/2Δ

τ and

t₂ ={τ

₀ }τ

_c +1/2Δ

τ

.

13. Method according to claim 10, wherein the function f(τ

_x) is defined as ##EQU39## wherein N is a normalization parameter, and Δ

τ

is a sampling time interval.

14. Method according to claim 10, wherein the function f(τ

_x) is defined as ##EQU40## wherein N is a normalization parameter, and Δ

τ

is a sampling time interval.

15. Method according to claim 1, wherein said samples R_x (τ

) are obtained simultaneously.

16. In a positioning signal processing system comprising a receiver which receives a data carrying carrier modulated pseudo random noise signal from a transmitter, a method of processing the received data carrying signal to determine a propagation time τ

₀ indicative of the time it takes the received signal to propagate directly from the transmitter to the receiver, the method comprising the steps of;

providing a reference signal;

performing a correlation operation on the received data carrying signal with said reference signal in order to provide a measure of the shape of the correlation function R_x (τ

);

iteratively estimating the multipath-induced contribution to the correlation function and subtracting said estimated multipath-induced contribution from the correlation function to obtain an estimated line of sight contribution to the correlation function; and

calculating the propagation time τ

₀ from said estimated line of sight contribution to the correlation function.

17. In a positioning signal processing system comprising a receiver which receives a data carrying carrier modulated pseudo random noise signal from a transmitter, a method of processing the received data carrying signal to determine an estimated propagation time τ

₀ indicative of the time it takes the received signal to propagate directly from the transmitter to the receiver, the method comprising the steps of;

A. providing a reference signal;

B. generating a plurality of time delayed versions of said reference signal wherein each of said delayed reference signals has been time delayed by a different amount;

C. performing, for each said plurality of delayed reference signals, a correlation operation with said received data carrying signal and providing a correlation signal value indicative thereof;

D. estimating a multipath-induced contribution for each of said plurality of correlation signal values;

E. subtracting said estimated multipath-induced contribution from each of said plurality of correlation signal values to obtain an estimated line of sight signal contribution for each of said plurality of correlation signal values;

F. estimating the propagation time τ

₀ from said estimated line of sight contributions; and

G. repeating said steps D, E and F until a suitable stop criterion has been reached.

18. The method of claim 17, wherein said step of estimating said multi-path-induced contributions includes the step of using maximum likelihood estimation to provide said estimated multipath-induced contributions.

19. In a signal processing system comprising a receiver which receives a data carrying carrier modulated pseudo random noise signal from a transmitter, a method of processing the received data carrying signal to determine an estimated propagation time τ

₀ indicative of the time it takes the received data carrying signal to propagate directly from the transmitter to the receiver, the method comprising the steps of;

providing a plurality of delayed reference signals wherein each of said delayed reference signals has been delayed in time by a different amount;

performing, for each of said plurality of delayed reference signals, a correlation operation with said received data carrying signal and providing a correlation signal value indicative thereof;

estimating a multipath-induced correlation signal value for each of said plurality of correlation signal values;

subtracting each of said estimated multipath-induced correlation signal values from its associated one of said plurality of correlation signal values to obtain an estimated line of sight signal correlation contribution for each of said plurality of correlation signal values; and

calculating the estimated propagation time τ

₀ from said estimated line of sight signal correlation contributions.

20. The method of claim 19, further comprising the step of:

serially repeating said steps of estimating, subtracting and calculating until a predetermined stop criterion is satisfied.

21. The method of claim 19, further, comprising the steps of:

estimating a second multipath-induced correlation signal value for each of said correlation signal values;

subtracting each of said estimated first and second multipath-induced correlation signal values from its associated one of said plurality of correlation signal values to obtain a second estimated line of sight signal correlation contribution for each of said plurality of correlation signal values; and

calculating the estimated propagation time τ

₀ from said second estimated line of sight signal correlation contributions.