High dynamic global positioning system receiver

• US 4,578,678 A
• Filed: 11/14/1983
• Issued: 03/25/1986
• Est. Priority Date: 11/14/1983
• Status: Expired due to Fees

• Current Assignee: The United States of America As Represented By The Secretary of Agriculture
• Original Assignee: United States Administrator of National Aeronautics and Space Administration
• Inventors: Hurd, William J.
• Primary Examiner(s): Blum, Theodore M.
• Assistant Examiner(s): Issing, Gregory C.
• Application Number:
  US06/551,536
• Time in Patent Office: 862 Days
• Field of Search: 343/357, 343/356, 343/417, 343/418, 343/420, 343/422
• US Class Current: 342/357.25
• CPC Class Codes: G01S 19/29 carrier including Doppler, related G01S19/246 takes precedence

Abstract

A GPS receiver having a number of channels, one for each of a plurality of satellites, receives an aggregate of pseudorange code time division modulated signals from some or all satellites and converts the aggregate to baseband and then to digital form for separate processing in the separate channels, each of which includes a complex mixer for adjustment of the baseband signal to a smooth estimate of frequency, ω, for the next measurement interval, and an N-lag cross correlator which receives a replica of the pseudorange code adjusted for a smoothed estimate of the pseudorange (delay), τ. A fast Fourier transform processor computes the signal energy as a function of Doppler frequency for each correlation lag, and a range and frequency estimator computes estimates of pseudorange, τ, and frequency, ω. These raw estimates from all channels are used to estimate receiver position, velocity, clock offset and clock rate offset in a conventional navigation and control unit, and based on the total solution, that unit computes smoothed estimates, τ and ω, for the next measurement interval. The smoothed estimates are thus based on the total solution, and not just on the separate estimates of the satellites, all without the use of any phase-locked loop.

Claims

1. In a receiver for a navigation system using a number of satellites which synchronously transmit pseudorange codes that are code division multiplexed on a carrier frequency, said receiver having means for amplifying the aggregate of the signals from said satellites, and converting said aggregate to baseband, means for demodulating said baseband to inphase and quadrature signals, means for converting said inphase and quadrature signals from analog to digital form, a separate signal processing channel for each satellite signal received, and navigation means responsive to pseudorange and frequency estimates from all of said channels for computing smoothed estimates of receiver position, velocity, clock offset and clock rate offset, a method for navigation tracking with said channels in parallel comprising the steps of

computing in each channel the maximum likelihood estimation of pseudorange, τ

, to a satellite, and the frequency, ω

, of the satellite at baseband during a measurement time,

computing smoothed estimates of receiver position, velocity, clock offset and clock rate offset from the maximum likelihood estimates τ and

ω

, and data from said navigation means to predict pseudorange τ and

frequency ω

for the next measurement time for each satellite,

adjusting the digital baseband inphase and quadrature signals to a predicted frequency for the next measurement time in response to said smoothed estimate ω

of frequency, and

adjusting the phase of a code generator for generating a replica of the satellite pseudorange code with a predicted delay for the computation of the maximum likelihood estimate of pseudorange, τ

, in the next measurement time in response to said smoothed estimate τ

of pseudorange.

2. A method as defined in claim 1 wherein computing the maximum likelihood estimation of pseudorange, τ

, and frequency, ω

, is accomplished by

computing the cross correlation of the digital inphase and quadrature signal against a local replica of the code for each of N cross correlation lags,

computing the signal energy as a function of Doppler frequency for each correlation lag,

determining the point of maximum energy for a distinct lag and frequency, and

from said point of time lag and Doppler frequency that resulted in a maximum energy, and also from adjacent energy points arranged in terms of lag versus frequency, determining said maximum likelihood estimation of pseudorange τ and

frequency, ω

.

3. In a receiver for a navigation system using a number of satellites which synchronously transmit pseudorange codes that are code division multiplexed on a carrier frequency, said receiver having means for amplifying the aggregate of the signals from said satellites, and converting said aggregate to baseband, means for demodulating said baseband to inphase and quadrature signals, means for converting said inphase and quadrature signals from analog to digital form, a separate signal processing channel for each satellite signal received, and navigation means responsive to pseudorange and frequency estimates from all of said channels for computing smoothed estimates of receiver position, velocity, clock offset and clock rate offset for navigation tracking with said channels in parallel, an improvement comprising

means for computing in each channel the maximum likelihood estimation of pseudorange, τ

, to a satellite, and the frequency, ω

, of the satellite at baseband during a measurement time,

means for computing smoothed estimates of receiver position, velocity, clock offset and clock rate offset, means for predicting pseudorange τ and

frequency ω

for the next measurement time for each satellite from the maximum likelihood estimates τ and

ω

, and also from said position, velocity, clock offset and clock rate offset data,

means for adjusting the digital baseband inphase and quadrature signals to a predicted frequency for the next measurement time in response to said smoothed estimate ω

of frequency, and

means for adjusting the phase of a code generator for generating a replica of the satellite pseudorange code with a predicted delay equal to said smoothed estimate of pseudorange, τ

, for the next measurement time.

4. An improvement as defined in claim 3 wherein said means for computing the maximum likelihood estimation of pseudorange, τ

, and frequency, ω

, is comprised of

means for computing the cross correlation of the digital inphase and quadrature signal against a local replica of the code for each of N cross correlation lags,

means for computing the signal energy as a function of Doppler frequency for each correlation lag,

means for determining the point of maximum energy for a distinct lag and frequency, and

from the point of time lag and Doppler frequency that resulted in a maximum energy, and further from adjacent energy points arrayed in terms of lag versus frequency, means for determining said maximum likelihood estimation of pseudorange τ and

frequency ω

.

5. A GPS receiver having K channels, one for each of K satellites, means for receiving an aggregate of pseudorange code time division modulated signals from said satellites, and means for converting the aggregate signal to baseband, means for converting said baseband signal to digital form for separate processing in said K channels, each of which includes a complex mixer for adjustment of the digital baseband signal to a predicted frequency, ω

, for the next measurement interval, means for generating a replica of the satellite pseudorange code adjusted for a smoothed estimate of pseudorange delay, τ

, an N-lag cross correlator responsive to said smoothed pseudorange delay, τ

, means for computing the signal energy as a function of Doppler frequency for each correlation lag, means for computing maximum likelihood estimates of pseudorange, τ

, and frequency, ω

, from the maximum signal energy point and adjacent points, means for using these estimates from all channels to estimate receiver position, velocity, clock offset and clock rate offset in a conventional manner, and based on the total navigation solution computing said smoothed estimates, τ and

ω

, for the next measurement interval.