CA2118013C - Navigation receiver with coupled signal-tracking channels - Google Patents
Navigation receiver with coupled signal-tracking channelsInfo
- Publication number
- CA2118013C CA2118013C CA002118013A CA2118013A CA2118013C CA 2118013 C CA2118013 C CA 2118013C CA 002118013 A CA002118013 A CA 002118013A CA 2118013 A CA2118013 A CA 2118013A CA 2118013 C CA2118013 C CA 2118013C
- Authority
- CA
- Canada
- Prior art keywords
- received signals
- platform
- signals
- estimates
- receiver
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 238000000034 method Methods 0.000 claims abstract description 87
- 238000005259 measurement Methods 0.000 claims abstract description 51
- 230000008569 process Effects 0.000 claims abstract description 46
- 230000008859 change Effects 0.000 claims abstract description 24
- 238000012937 correction Methods 0.000 claims description 54
- 239000011159 matrix material Substances 0.000 claims description 52
- 238000004364 calculation method Methods 0.000 claims description 21
- 230000006870 function Effects 0.000 claims description 11
- 239000000969 carrier Substances 0.000 claims description 8
- 238000004891 communication Methods 0.000 claims description 8
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 230000000644 propagated effect Effects 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims 12
- 230000033001 locomotion Effects 0.000 abstract description 9
- 239000013598 vector Substances 0.000 description 63
- 238000012545 processing Methods 0.000 description 23
- 230000005540 biological transmission Effects 0.000 description 10
- 238000005192 partition Methods 0.000 description 10
- 230000014509 gene expression Effects 0.000 description 6
- 230000007704 transition Effects 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 238000013519 translation Methods 0.000 description 5
- 230000014616 translation Effects 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000013213 extrapolation Methods 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000012886 linear function Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000000873 masking effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- PRPINYUDVPFIRX-UHFFFAOYSA-N 1-naphthaleneacetic acid Chemical compound C1=CC=C2C(CC(=O)O)=CC=CC2=C1 PRPINYUDVPFIRX-UHFFFAOYSA-N 0.000 description 1
- XUKUURHRXDUEBC-KAYWLYCHSA-N Atorvastatin Chemical compound C=1C=CC=CC=1C1=C(C=2C=CC(F)=CC=2)N(CC[C@@H](O)C[C@@H](O)CC(O)=O)C(C(C)C)=C1C(=O)NC1=CC=CC=C1 XUKUURHRXDUEBC-KAYWLYCHSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- 101100021570 Caenorhabditis elegans lon-3 gene Proteins 0.000 description 1
- 241000726409 Satellites Species 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000000872 buffer Substances 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 238000005183 dynamical system Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000005433 ionosphere Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/53—Determining attitude
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/53—Determining attitude
- G01S19/54—Determining attitude using carrier phase measurements; using long or short baseline interferometry
Abstract
The navigation receiver with coupled-tracking channels is intended for use in radio navigation systems based on trilatera-tion wherein the times of arrival of radio signals transmitted by a plurality of ground- or space-based terminals are measured by a user terminal and utilized by the user terminal in calculating its position and orientation. The coupled-tracking navigation re-ceiver periodically measures carrier phase, carrier frequency, modulation phase, and carrier amplitude for all of the signals arriving at the receiving ports of the receiver and periodically es-timates the present values of carrier phase, carrier frequency, and modulation phase for all of the received signals, the estimat-ing process utilizing for each parameter estimate the parameter measurements for a plurality of the received signals properly combined in a statistically appropriate manner by taking into ac-count the relative geometry of the line-of-sight paths, receiver clock time dynamics, and dynamics and motion constraints of the receiver platform, thereby obtaining better performance un-der poor signal reception conditions and more accurate esti-mates of carrier phase, carrier frequency, and modulation phase for each of the received signals than independent measurements alone could provide. It follows that these more accurate esti-mates of the basic signal parameters lead to more accurate esti-mates of platform position and attitude and the rates of change of these quantities.
Description
W O 93/22691 -PC~r/US93/04150 D E S C R I P T I O N
NAVIGATION RECEIVER WITH COUPLED SIGNAL-TRACKING CHANNELS
TECHNICAL FIELD
The invention relates generally to radio navigation systems and more particularly to systems based on trilateration wherein radio signals transmitted by a plurality of ground- or space-based terminals are received by a user terminal and utilized by the user terminal in calculating its position and orientation.
BACKGROUND ART
Trilateration navigation systems are based on the geometrical principle that the position of a point in space is uniquely determined by the distances of the point from three other suitably-positioned reference points whose locations in space are known. Radio navigation systems based on this principle utilize transmitting terminals at the reference points which transmit signals in which a common time reference is embedded. The user terminal, by measuring the time of arrival of the reference terminal signals with respect to the same common time reference, can determine its range to each of the refere~ce transmitting terminals and thereby its position.
If the us~er terminal does not have a clock of sufficient accuracy for independently keeping track of the common time reference, it can, by measuring the time of arrival of signals from four reference transmitters, also keep its cloc~
synchronized to the common time reference of the reference transmitters.
The principle can be extended to the measurement of the attitude of a platform by utilizing receiving ports at a plurality of points on the platform and measuring the locations of these points relative to the reference terminals.
3~ In general, the navigation data desired includes not onl~
position and attitude but also the rates of change of these parameters, i.e. linear velocity and angular velocity. To achieve these goals, not only must the carrier modulation be ~C
NAVIGATION RECEIVER WITH COUPLED SIGNAL-TRACKING CHANNELS
TECHNICAL FIELD
The invention relates generally to radio navigation systems and more particularly to systems based on trilateration wherein radio signals transmitted by a plurality of ground- or space-based terminals are received by a user terminal and utilized by the user terminal in calculating its position and orientation.
BACKGROUND ART
Trilateration navigation systems are based on the geometrical principle that the position of a point in space is uniquely determined by the distances of the point from three other suitably-positioned reference points whose locations in space are known. Radio navigation systems based on this principle utilize transmitting terminals at the reference points which transmit signals in which a common time reference is embedded. The user terminal, by measuring the time of arrival of the reference terminal signals with respect to the same common time reference, can determine its range to each of the refere~ce transmitting terminals and thereby its position.
If the us~er terminal does not have a clock of sufficient accuracy for independently keeping track of the common time reference, it can, by measuring the time of arrival of signals from four reference transmitters, also keep its cloc~
synchronized to the common time reference of the reference transmitters.
The principle can be extended to the measurement of the attitude of a platform by utilizing receiving ports at a plurality of points on the platform and measuring the locations of these points relative to the reference terminals.
3~ In general, the navigation data desired includes not onl~
position and attitude but also the rates of change of these parameters, i.e. linear velocity and angular velocity. To achieve these goals, not only must the carrier modulation be ~C
2 1 1 8 0 1 3 PCT/US93/04150 tracked but also the carrier phase and frequency for each of the received signals.
User receiver systems for trilateration radio navigation systems are comprised of an antenna, a receiver "front end"
for selecting and amplifying the frequency band of interest together with "down" converters for translating the radio frequency band of interest to an intermediate frequency band more suitable for subsequent signal processing, and signal processors for extracting the desired information from the received signals. Commonly today, the received signals are translated to baseband where in-phase (I) and quadrature (Q) signal samplers and digitizers followed by digital signal processors are used to e~tract range and range rate data from the signals.
The extraction of range and range rate data from the I
and Q samples of each received signal begins with a correlation process. The I and Q samples are multiplied by reference signal samples, the reference signal being a locally-generated replica of the transmitted signal, and the resulting products are accumulated with prior products to obtain smoothed I and Q values with greatly reduced noise and interference.
The smoothed I and Q samples are further processed to obtain information as to the offset, or misalignment, between the carrier phase, carrier frequency, and modulation phase of a received signal and the corresponding parameters of the corresponding locally-generated reference signal. Two additional functions performed at this stage are estimation of received signal amplitude and recovery of any au~iliary data modulation that may have been impressed upon the carrier.
Within the conventional receiver, all of these processing steps are carried out separately for each received signal.
Tracking loops are provided that provide estimated values of carrier phase, carrier frequency, and modulation phase.
After processing the current set of smoothed I and Q data, the tracking loops are updated resulting in updated estimates of carrier phase and frequency and modulation phase relative to the local clock in the receiver. At the same time, the WO93/22691 2 1 18 0 1~ PCT/US93/04150 auxiliary quantities of signal amplitude and carrier data bits are updated. Prior to the start of the ne~t correlation interval the updated estimates of carrier phase and frequenc and modulation phase are used to derive a corrective command signal that is used to adjust the carrier phase and frequency and the modulation phase of the reference signal generator in the correlation processor for the purpose of ensuring adequate signal correlation during the subsequent correlation interval.
An example of the type of navigation system described above is the NAVSTAR global positioning system (GPS) developed by the United States Government. The NAVSTAR system consists of a constellation of 18 to 24 orbiting satellites that transmit pseudo-random ranging signals from which users with appropriate equipment can obtain three-dimensional location, velocity, and timing information anywhere on or near the surface of the earth. Details concerning NAVSTAR/GPS are given in NAVIGATION: Journal of the Institution of Na-~igation, Volume 25, Number 2, December, 1978 (entire issue).
Conventional navigation receivers process the plurality of received signals independently. Thus, conventional navigation receivers fail to realize their theoretical performance potential by not taking advantage of the collective information contained in the plurality of received signals in the processing of each of the individual signals.
2~
DISCLOSURE OF INVENTION
The navigation receiver with coupled signal-tracking channels operates on a platform within line-of-sight range of a plurality of spatially-distributed transmitters with positions known to the receiver. The transmitters transmit modulated sinusoidal carriers that are particular functions of time as maintained by synchronized transmitter clocks. The navigation receiver, knowing the time associated with the transmission of the signals and measuring the reception time
User receiver systems for trilateration radio navigation systems are comprised of an antenna, a receiver "front end"
for selecting and amplifying the frequency band of interest together with "down" converters for translating the radio frequency band of interest to an intermediate frequency band more suitable for subsequent signal processing, and signal processors for extracting the desired information from the received signals. Commonly today, the received signals are translated to baseband where in-phase (I) and quadrature (Q) signal samplers and digitizers followed by digital signal processors are used to e~tract range and range rate data from the signals.
The extraction of range and range rate data from the I
and Q samples of each received signal begins with a correlation process. The I and Q samples are multiplied by reference signal samples, the reference signal being a locally-generated replica of the transmitted signal, and the resulting products are accumulated with prior products to obtain smoothed I and Q values with greatly reduced noise and interference.
The smoothed I and Q samples are further processed to obtain information as to the offset, or misalignment, between the carrier phase, carrier frequency, and modulation phase of a received signal and the corresponding parameters of the corresponding locally-generated reference signal. Two additional functions performed at this stage are estimation of received signal amplitude and recovery of any au~iliary data modulation that may have been impressed upon the carrier.
Within the conventional receiver, all of these processing steps are carried out separately for each received signal.
Tracking loops are provided that provide estimated values of carrier phase, carrier frequency, and modulation phase.
After processing the current set of smoothed I and Q data, the tracking loops are updated resulting in updated estimates of carrier phase and frequency and modulation phase relative to the local clock in the receiver. At the same time, the WO93/22691 2 1 18 0 1~ PCT/US93/04150 auxiliary quantities of signal amplitude and carrier data bits are updated. Prior to the start of the ne~t correlation interval the updated estimates of carrier phase and frequenc and modulation phase are used to derive a corrective command signal that is used to adjust the carrier phase and frequency and the modulation phase of the reference signal generator in the correlation processor for the purpose of ensuring adequate signal correlation during the subsequent correlation interval.
An example of the type of navigation system described above is the NAVSTAR global positioning system (GPS) developed by the United States Government. The NAVSTAR system consists of a constellation of 18 to 24 orbiting satellites that transmit pseudo-random ranging signals from which users with appropriate equipment can obtain three-dimensional location, velocity, and timing information anywhere on or near the surface of the earth. Details concerning NAVSTAR/GPS are given in NAVIGATION: Journal of the Institution of Na-~igation, Volume 25, Number 2, December, 1978 (entire issue).
Conventional navigation receivers process the plurality of received signals independently. Thus, conventional navigation receivers fail to realize their theoretical performance potential by not taking advantage of the collective information contained in the plurality of received signals in the processing of each of the individual signals.
2~
DISCLOSURE OF INVENTION
The navigation receiver with coupled signal-tracking channels operates on a platform within line-of-sight range of a plurality of spatially-distributed transmitters with positions known to the receiver. The transmitters transmit modulated sinusoidal carriers that are particular functions of time as maintained by synchronized transmitter clocks. The navigation receiver, knowing the time associated with the transmission of the signals and measuring the reception time
3~ by its own clock synchronized to the transmitter clocks, can determine its range to each transmitter. By measuring the Doppler shifts in carrier fre~uency, the receiver can determine the rate of change of range to each transmitter.
From this data the position and rate of change of position of the receiver and its platform can be determined. If the platform is equipped with multiple receiving ports, the receiver can also determine platform attitude and rate of change of attitude by including the additional line-of-sight paths in the measurement process.
The present invention pertains to the signal processor in the receiver which periodically measures carrier phase, carrier frequency, modulation phase, and carrier amplitude for all of the signals arriving at the platform receiving ports and periodically estimates the present values of carrier phase, carrier frequency, and modulation phase for all of the received signals, the estimating process utilizing for each parameter estimate the parameter measurements for a plurality of the received signals properly combined in a statistically appropriate manner by taking into account the relative geometry of the line-of-sight paths, receiver clock time dynamics, and dynamics and motion constraints of the receiver platform, thereby obtaining more accurate estimates of carrier phase, carrier frequency, and modulation phase for each of the received signals than independent measurements alone could provide. It follows that these more accurate estimates of the basic signal parameters lead to more accurate estimates of platform position and attitude and the rates of change of these quantities.
One object of the invention is to obtain more accurate and continuous updates of real-time position, velocit~, attitude, and attitude rate and the recovery of au~iliary data modulation down to a lower signal threshold by using a complete model for line-of-sight dynamics in such a manner that the carrier phase, carrier frequency, and modulation phase can be more reliably tracked during periods of intermittent signal attenuation, noise interference, multipath interference, and signal blockage and outage.
Another object of the invention is to provide impro~ed performance for post-processed precise carrier applications such as land, marine, or air survey by employing the above W O 93/22691 2 1 1 8 0 1 3 PC~r/US93/~4150 -described techniques of signal processing for combining wideband signal data from several locations.
BRIEF DESCRIPTION OF DRAWINGS
FIG. 1 is a functional block diagram of a signal processing apparatus for a conventional navigation receiver showing the decoupled signal-tracking structure characteristic of these receivers.
FIG. 2 is a functional block diagram of the sianal processing apparatus for the preferred embodiment of the present invention showing the coupled signal-tracking structure.
FIG. 3 is an illustration of the functional relationships of the software routines that control the operation of the coupled filter processor in the preferred embodiment of t}le invention.
FIG. 4 shows the timing of the signal processing operations performed in the signal processing apparatus of the preferred embodiment of the present invention.
FIG. 5 provides the definitions of the navigation state vector and the LOS state vector.
FIG. 6 describes how the base and incremental kOS states evolve with discrete time and defines the LOS state transition matrix.
FIG. 7 defines the input data vector to the coupled-filter signal processor and shows how the data vector can be modeled as a non-linear function of the incremental LOS state vector.
FIG. 8 further defines the non-linear function introduced in Fig. 7.
FIG. 9 continues the definition of the non-linear function introduced in Fig. 7 and relates the velocity aiding term supplied by the external velocity sensor to certain components of the incremental LOS state vector.
FIG. 10 details the operations carried out by the coupled-filter signal processor in performing the sianal amplitude estimator routine.
FIG. 11 identifies the operations carried out by the : 2118013 coupled-filter signal processor in performing the telemetry data demodulator routine.
FIG. 12 details the initial operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. 13 continues the detailing of the operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. 14 continues the ~etailing of the operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. 15 continues the detailing of the operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. 16 completes the detailing of the operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. l7 details the initial operations carried out by the coupled-filter signal processor in performing the IQ combiner routine.
FIG. 18 completes the detailing of the operations carried out by the coupled-filter signal processor in performing the IQ combiner routine.
FIG. l9 details the initial operations carried out by the coupled-filter signal processor in performing the vehicle position and attitude evaluator routine.
FIGS. 20, 21, 22, and 23 continue the detailing of the operations carried out by the coupled-filter signal processor in performing the vehicle position and attitude evaluator routine.
FIG. 24 completes the detailing of the operations carried out by the coupled-filter signal processor in performing the vehicle position and attitude evaluator routine.
FIG. 25 details the operations carried out by the 3~ coupled-fllter signal processor in performing the reference generator command setup routine.
BEST MODE FOR CARRYING OUT THE INVENTION
The preferred embodiment of the radio navigation receiver with coupled signal-tracking channels will be described in the context of its application in a satellite-referenced radio navigation system as exemplified by NAVSTAR/GPS.
A number of vsriants of a commonly used signal trac~ing technique used in GPS and other time-of-arrival navigation systems have been described in patents and other publications over the years. All such variants contain an antenna and receiver front end for selecting the general frequency band of the navigation system and down-converter stages followed by additional frequency- selective elements for further reducing noise and interference. In the case of multiple receiving antennas, parallel sets of these components are employed. Ir-J
present-day receiver art it is now common practice to sample in-phase and quadrature mix-down outputs at baseband and to carry out the remaining signal processing operations with digital hardware and microprocessors. For convenience in exposition, a similar digital implementation will form the basis for the description of the preferred embodiment of the present invention. It will be understood by one s~illed in the art that the principles of the invention can be applied with equal effectiveness to receiver implementations based on combinations of analog and digital techniques.
The demodulator portion of a conventional receiver for a single-antenna user terminal is shown in Fig. 1. The digital representations of the in-phase and quadrature signal samples 106 are routed to four correlator channels 110, 112, 11~, and 116, one for each of the four reference transmitter signals.
Each correlator channel consists of a correlator 118, a reference generator 120, and a loop filter 122. The samples of the composite received signal are multiplied in correlator 118 by samples of a reference signal supplied by the reference generator 120 and summed over the correlation interval to 3~ obtain smoothed I and Q values with greatly reduced noise and interference.
The reference signal is a replica of the signal broadcast by the particular reference transmitter that is associated 801~
with the particular correlator channel. In the case of GPS, the reference generator samples are generated by digital hardware consisting of spread spectrum code generators and carrier phase generators.
5The I and Q values resulting from signal correlation are further processed in loop filter 122 to obtain information as to the offset, or misalignment, between the incoming carrier phase and frequency and modulating code and the corresponding reference signal quantities. Two additional functions performed by the loop filter 122 are the estimation of received signal amplitude and recovery of any auxiliary data modulation that may be impressed upon the carrier. Within the conventional receiver, all of these processino steps are carried out separately for each path defined by a particular reference transmitter and a particular user antenna.
For carriers upon which the data has been impressed by phase-shift keying, the I and Q correlator outputs are typically processed by a loop filter 122 with a Costas loop.
If the carrier is further modulated at the transmitter by a spread-spectrum signal, the I and Q correlator outputs are typically also processed with a delay-lock loop or tau-dither loop. New I and Q data causes the loop filter 122 to update its estimates of carrier phase and frequency and modulating code, relative to the local clock in the receiver, and forward these values to the navigation estimator 150. At the same time, the auxiliary quantities of signal amplitude and carrier data bits are updated. Prior to the start of the ne~t correlation interval the updated estimates of carrier phase and frequency and modulating code offset are used by loop filter 122 to derive a corrective command signal to refere~ce generator 120 for the purpose of ensuring adequate signal correlation during the next correlation interval.
Further details concerning conventional navigation receivers of the type discussed above and the signal processing techniques employed therein can by found in U.S
patents 4,445,118, 4,468,793, 4,485,393, 4,578,678, 4,754,465,
From this data the position and rate of change of position of the receiver and its platform can be determined. If the platform is equipped with multiple receiving ports, the receiver can also determine platform attitude and rate of change of attitude by including the additional line-of-sight paths in the measurement process.
The present invention pertains to the signal processor in the receiver which periodically measures carrier phase, carrier frequency, modulation phase, and carrier amplitude for all of the signals arriving at the platform receiving ports and periodically estimates the present values of carrier phase, carrier frequency, and modulation phase for all of the received signals, the estimating process utilizing for each parameter estimate the parameter measurements for a plurality of the received signals properly combined in a statistically appropriate manner by taking into account the relative geometry of the line-of-sight paths, receiver clock time dynamics, and dynamics and motion constraints of the receiver platform, thereby obtaining more accurate estimates of carrier phase, carrier frequency, and modulation phase for each of the received signals than independent measurements alone could provide. It follows that these more accurate estimates of the basic signal parameters lead to more accurate estimates of platform position and attitude and the rates of change of these quantities.
One object of the invention is to obtain more accurate and continuous updates of real-time position, velocit~, attitude, and attitude rate and the recovery of au~iliary data modulation down to a lower signal threshold by using a complete model for line-of-sight dynamics in such a manner that the carrier phase, carrier frequency, and modulation phase can be more reliably tracked during periods of intermittent signal attenuation, noise interference, multipath interference, and signal blockage and outage.
Another object of the invention is to provide impro~ed performance for post-processed precise carrier applications such as land, marine, or air survey by employing the above W O 93/22691 2 1 1 8 0 1 3 PC~r/US93/~4150 -described techniques of signal processing for combining wideband signal data from several locations.
BRIEF DESCRIPTION OF DRAWINGS
FIG. 1 is a functional block diagram of a signal processing apparatus for a conventional navigation receiver showing the decoupled signal-tracking structure characteristic of these receivers.
FIG. 2 is a functional block diagram of the sianal processing apparatus for the preferred embodiment of the present invention showing the coupled signal-tracking structure.
FIG. 3 is an illustration of the functional relationships of the software routines that control the operation of the coupled filter processor in the preferred embodiment of t}le invention.
FIG. 4 shows the timing of the signal processing operations performed in the signal processing apparatus of the preferred embodiment of the present invention.
FIG. 5 provides the definitions of the navigation state vector and the LOS state vector.
FIG. 6 describes how the base and incremental kOS states evolve with discrete time and defines the LOS state transition matrix.
FIG. 7 defines the input data vector to the coupled-filter signal processor and shows how the data vector can be modeled as a non-linear function of the incremental LOS state vector.
FIG. 8 further defines the non-linear function introduced in Fig. 7.
FIG. 9 continues the definition of the non-linear function introduced in Fig. 7 and relates the velocity aiding term supplied by the external velocity sensor to certain components of the incremental LOS state vector.
FIG. 10 details the operations carried out by the coupled-filter signal processor in performing the sianal amplitude estimator routine.
FIG. 11 identifies the operations carried out by the : 2118013 coupled-filter signal processor in performing the telemetry data demodulator routine.
FIG. 12 details the initial operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. 13 continues the detailing of the operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. 14 continues the ~etailing of the operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. 15 continues the detailing of the operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. 16 completes the detailing of the operations carried out by the coupled-filter signal processor in performing the filter gain adjust routine.
FIG. l7 details the initial operations carried out by the coupled-filter signal processor in performing the IQ combiner routine.
FIG. 18 completes the detailing of the operations carried out by the coupled-filter signal processor in performing the IQ combiner routine.
FIG. l9 details the initial operations carried out by the coupled-filter signal processor in performing the vehicle position and attitude evaluator routine.
FIGS. 20, 21, 22, and 23 continue the detailing of the operations carried out by the coupled-filter signal processor in performing the vehicle position and attitude evaluator routine.
FIG. 24 completes the detailing of the operations carried out by the coupled-filter signal processor in performing the vehicle position and attitude evaluator routine.
FIG. 25 details the operations carried out by the 3~ coupled-fllter signal processor in performing the reference generator command setup routine.
BEST MODE FOR CARRYING OUT THE INVENTION
The preferred embodiment of the radio navigation receiver with coupled signal-tracking channels will be described in the context of its application in a satellite-referenced radio navigation system as exemplified by NAVSTAR/GPS.
A number of vsriants of a commonly used signal trac~ing technique used in GPS and other time-of-arrival navigation systems have been described in patents and other publications over the years. All such variants contain an antenna and receiver front end for selecting the general frequency band of the navigation system and down-converter stages followed by additional frequency- selective elements for further reducing noise and interference. In the case of multiple receiving antennas, parallel sets of these components are employed. Ir-J
present-day receiver art it is now common practice to sample in-phase and quadrature mix-down outputs at baseband and to carry out the remaining signal processing operations with digital hardware and microprocessors. For convenience in exposition, a similar digital implementation will form the basis for the description of the preferred embodiment of the present invention. It will be understood by one s~illed in the art that the principles of the invention can be applied with equal effectiveness to receiver implementations based on combinations of analog and digital techniques.
The demodulator portion of a conventional receiver for a single-antenna user terminal is shown in Fig. 1. The digital representations of the in-phase and quadrature signal samples 106 are routed to four correlator channels 110, 112, 11~, and 116, one for each of the four reference transmitter signals.
Each correlator channel consists of a correlator 118, a reference generator 120, and a loop filter 122. The samples of the composite received signal are multiplied in correlator 118 by samples of a reference signal supplied by the reference generator 120 and summed over the correlation interval to 3~ obtain smoothed I and Q values with greatly reduced noise and interference.
The reference signal is a replica of the signal broadcast by the particular reference transmitter that is associated 801~
with the particular correlator channel. In the case of GPS, the reference generator samples are generated by digital hardware consisting of spread spectrum code generators and carrier phase generators.
5The I and Q values resulting from signal correlation are further processed in loop filter 122 to obtain information as to the offset, or misalignment, between the incoming carrier phase and frequency and modulating code and the corresponding reference signal quantities. Two additional functions performed by the loop filter 122 are the estimation of received signal amplitude and recovery of any auxiliary data modulation that may be impressed upon the carrier. Within the conventional receiver, all of these processino steps are carried out separately for each path defined by a particular reference transmitter and a particular user antenna.
For carriers upon which the data has been impressed by phase-shift keying, the I and Q correlator outputs are typically processed by a loop filter 122 with a Costas loop.
If the carrier is further modulated at the transmitter by a spread-spectrum signal, the I and Q correlator outputs are typically also processed with a delay-lock loop or tau-dither loop. New I and Q data causes the loop filter 122 to update its estimates of carrier phase and frequency and modulating code, relative to the local clock in the receiver, and forward these values to the navigation estimator 150. At the same time, the auxiliary quantities of signal amplitude and carrier data bits are updated. Prior to the start of the ne~t correlation interval the updated estimates of carrier phase and frequency and modulating code offset are used by loop filter 122 to derive a corrective command signal to refere~ce generator 120 for the purpose of ensuring adequate signal correlation during the next correlation interval.
Further details concerning conventional navigation receivers of the type discussed above and the signal processing techniques employed therein can by found in U.S
patents 4,445,118, 4,468,793, 4,485,393, 4,578,678, 4,754,465,
4,785,463, 4,821,294, and 4,928,106. The above patents are incorporated in this specification by reference.
-It should be noted that in the prior art shown in Fig. 1 all of the processing steps are carried out in a decoupled fashion. There is no interaction between channels assigned to different transmitters. In essence, this type of structure
-It should be noted that in the prior art shown in Fig. 1 all of the processing steps are carried out in a decoupled fashion. There is no interaction between channels assigned to different transmitters. In essence, this type of structure
5 ignores statistical correlations in path dynamics and penalizes signal tracking performance, particularly during signal attenuations or while operating in the presence of strong interference. By coupling the signal-trac~ing channels together as shown in Fig. 2, it is possible to achieve superior signal-tracking performance under adverse reception conditions.
The demodulator portion of the preferred embodiment of a na~rigation receiver having coupled signal-tracking channels is shown in Fig. 2. The front end of the receiver is conventional and is full~ described in the patents incorporated by reference. As was mentioned previously, in present-day receiver art it is now common practice to sample in-phase and quadrature mix-down outputs at baseband and to carry out the remaining signal processing operations with digital hardware and microprocessors. These sampled in-phase and quadrature outputs from a conventional receiver front-end form the inputs 206 to the correlation system 208.
The correlator system 208 comprises the four correlator channels 210, 212, 214, and 216 where each of the correlator channels consists of a correlator 218 and a reference generator 220. The correlators 218 and the reference generators 220 are essentially the same as those used in the conventional decoupled receiver demodulators illustrated in Fig. 1. Details concerning the design of these receiver elements can also be found in the aforementioned patent references.
The essence of the invention is contained within the coupled filter processor 230 which in combination with the correlator system 208 produce navigational data more effectively as compared to the conventional uncoupled demodulator structure of Fig. 1.
The set of I and Q data from the correlator system 208 provides the basis for calculations performed by the coupled filter processor 230. In the GPS navigation conte~t, a fresh data set is typically made available on the data bus 232 at 20-millisecond intervals. The coupled filter processor, unlike a conventional phase-lock loop, Costas loop, or delay-lock loop, merges the data set in such a fashion that commandsfed back to each reference generator 220 over the control bus 234 are the result of the set of input bus data from all correlator channels during the previous correlation interval.
It should be noted that when multiple receiving antennas are employed, each antenna would feed its own receiver front end and its own correlation system 208. All correlation systems 208 would communicate with the coupled filter processor 230 by means of the data bus 232 and the control bus 234.
The velocity aiding system 235 provides measurements of platform velocity as projected onto the set of available LOS
paths between satellites and platform receiving antennas which are utilized by the coupled-filter signal processor 230 in cases of extreme signal blockage. In the case of a land vehicle these velocity components would typically be derived from a flux-gate compass and speedometer. In the case of a marine vehicle the velocity aiding information would typically be derived from a flux-gate compass and water speed indicator.
In the case of an aircraft these components might be derived from an inertial measurement unit (IMU). The velocity aiding system communicates with the signal processor by means of the data bus 232 and the control bus 234.
The preferred embodiment of the coupled filter processor 230 is one or more microprocessors of the Intel 8048~ or Motorola 68040 variety connected to the correlator system 208 by data bus 232 and control bus 234 and performing the software routines identified in Fig. 3.
Typically at 20-millisecond intervals the raw data vector Z_data from the correlator system 208 and the velocity aiding system 235 is sent over the data bus 232 for processinG. The format of this data is typically 8 bits for each component of the data vector. For each correlator channel there are W O 93/22691 2 118 0 13 PC~r/US93/041S0 typically six data components having to do with the measured correlations and a seventh component that is the projection of platform velocity onto the appropriate line-of-sight path.
The correlation components are known in the art as (1) punctual code I and Q outputs, (2) delay code I and Q outputs, and ~3) advanced code I and Q outputs, correspondina to the accumulated correlation values for in-phase and quadrature phase references and punctual, advance, and delay code references respectively.
The data vector Z_data is processed by the coupled filter signal processor 230 (Fig. 2) in accordance with the six routines identified in Fig. 3, the results of the processing being a set of commands conveyed to the reference generators 220 (Fig. 2), typically at 20-millisecond intervals, by means of the control bus 234. Each reference generator command is typically a 32-bit word comprised of values for reference signal phase and frequency and modulation phase at the beginning of the next correlation interval. The array of reference generator commands is denoted by the symbol sf.x_control.
Other inputs to the coupled-filter signal processor ~30 (Fig. 2) are the remotely-generated correction signals 236 from a differential correction station and the system initialization data 238.
The remotely-generated correction signals 236 are particularly important in compensating for uncertainties in phase of the received signals. In GPS such uncertainties are caused by deviations of the actual orbital model parameters from their nominal broadcast values, by drifts in the satellite clocks, and by unknown ionospheric and tropospheric delays. The deviations in carrier phase and frequency and modulation phase of the received signals are seen nearly identically by a fixed-site calibration station and the platform using the navigation system, and drift estimates may be broadcast to the platform where they can be removed by the coupled-filter signal processor. This mode of operation, known in the GPS art as differential GPS o.peration (DGPS) and described fully in the literature under industry standard W O 93/22691 21 1 80 1 3~ ` ` ; PC~r/US93/04150 format RTCM-104, release 2.0, (Radio Technical Commission for Maritime Services, P.O. Box 19087, Washington, D.C. 20036) is implemented in the present invention by applying the corrections directly within the signal tracking process rather than down-stream in the navigation calculations as is customary in the art.
The system initialization data 238 consists of time constants alpha_filt and process noise constants q_filt describing platform motionj parameters specifying measurement noise levels sf.r expected from the correlator system 208 (Fig. 2), and initial estimates of platform navigation states sf.part.
The outputs of the coupled-filter signal processor 230 (Fig. 2) are platform navigation data (or states) 242. Fol the case of a platform having a single antenna, these states are typically easting, northing, and altitude and the rates of change of these quantities. Additionally, the signal processor may estimate the offset between the internal system clock and the time reference of the navigation system. These outputs are of the same type as a conventional navigation system.
For the case of a platform having multiple antennas, platform orientation and motion states may be provided as outputs. Typically these are platform roll, pitch, and ya~-2~ angles with respect to the local east-north-up datum and the rates of change of these angles.
The scheduling of interactions between the correlator system 208 (Fig. 2) and the coupled-filter signal processor 230 (Fig. 2) is illustrated in Fig. 4. The figure depicts three frames of interactions: the (k-l)'th frame, the l~'th frame, and the (k+l)'th frame. Each frame begins at a time denoted by tl with an interrupt signal provided b~- the correlation system to the signal processor that signals the availabilily of data on the data bus 232 from one of the 3~ correlator channels 210, 212, 214, and 216 (Fig. 2). Tlle signal processor reads the data during the period tl-t2, performs the routines identified in Fig. 3 during the period t2-t3, and transfers the resulting command data to the WO93/22691 , PCT/US93/04150 - 2118013~
correlator system over control bus 234 (Fig. 2) during the interval t3-t4.
The command data consists of best estimates of carrier phase, carrier frequency offset, and modulating code offset to be used in adjusting the reference generator of the correlator channel in the correlator system. These estimates take into account satellite motion, DGPS biases, velocity aiding (if available), and dispersive propagation effects which are manifested by a difference in the LOS distance derived from the carrier phase measurement and the LOS distance derived from the modulating code offset measurement.
The next correlation interval begins at t4 and ends at the next tl when the next processing cycle begins.
Correlation intervals are normally loc~ed to the telemetry data bit transitions resulting in asynchrol~ous correlation intervals between correlator channels.
Specifically, during each system cycle, i.e. the time period required to service all correlator channels, the signal processor services the correlator channels on a "first come, first-served" basis. After servicing a correlator channel interrupt, the signal processor services the next correlator channel to generate an interrupt. The processing steps performed are identical for each correlator channel. ~ system cycle time for GPS is approximately 20 milliseconds.
The operations performed by the coupled-filter signal processor 230 (Fig. 2) and the routines in which they appear are as follows:
Step OPeration Routine 1 Read correl'ator data. sig. amp. est.
2 Estimate signal amplitude. sig. amp. est.
3 Recover telemetry data bit. telem. data demod.
4 Extrapolate LOS state from most filter oain recently updated channel. adjust WOg3/2269l 21 18 013 ~ pcT/uss3/o4lsn 5 Update residuals and evaluate signal filter gain quality. adjust
The demodulator portion of the preferred embodiment of a na~rigation receiver having coupled signal-tracking channels is shown in Fig. 2. The front end of the receiver is conventional and is full~ described in the patents incorporated by reference. As was mentioned previously, in present-day receiver art it is now common practice to sample in-phase and quadrature mix-down outputs at baseband and to carry out the remaining signal processing operations with digital hardware and microprocessors. These sampled in-phase and quadrature outputs from a conventional receiver front-end form the inputs 206 to the correlation system 208.
The correlator system 208 comprises the four correlator channels 210, 212, 214, and 216 where each of the correlator channels consists of a correlator 218 and a reference generator 220. The correlators 218 and the reference generators 220 are essentially the same as those used in the conventional decoupled receiver demodulators illustrated in Fig. 1. Details concerning the design of these receiver elements can also be found in the aforementioned patent references.
The essence of the invention is contained within the coupled filter processor 230 which in combination with the correlator system 208 produce navigational data more effectively as compared to the conventional uncoupled demodulator structure of Fig. 1.
The set of I and Q data from the correlator system 208 provides the basis for calculations performed by the coupled filter processor 230. In the GPS navigation conte~t, a fresh data set is typically made available on the data bus 232 at 20-millisecond intervals. The coupled filter processor, unlike a conventional phase-lock loop, Costas loop, or delay-lock loop, merges the data set in such a fashion that commandsfed back to each reference generator 220 over the control bus 234 are the result of the set of input bus data from all correlator channels during the previous correlation interval.
It should be noted that when multiple receiving antennas are employed, each antenna would feed its own receiver front end and its own correlation system 208. All correlation systems 208 would communicate with the coupled filter processor 230 by means of the data bus 232 and the control bus 234.
The velocity aiding system 235 provides measurements of platform velocity as projected onto the set of available LOS
paths between satellites and platform receiving antennas which are utilized by the coupled-filter signal processor 230 in cases of extreme signal blockage. In the case of a land vehicle these velocity components would typically be derived from a flux-gate compass and speedometer. In the case of a marine vehicle the velocity aiding information would typically be derived from a flux-gate compass and water speed indicator.
In the case of an aircraft these components might be derived from an inertial measurement unit (IMU). The velocity aiding system communicates with the signal processor by means of the data bus 232 and the control bus 234.
The preferred embodiment of the coupled filter processor 230 is one or more microprocessors of the Intel 8048~ or Motorola 68040 variety connected to the correlator system 208 by data bus 232 and control bus 234 and performing the software routines identified in Fig. 3.
Typically at 20-millisecond intervals the raw data vector Z_data from the correlator system 208 and the velocity aiding system 235 is sent over the data bus 232 for processinG. The format of this data is typically 8 bits for each component of the data vector. For each correlator channel there are W O 93/22691 2 118 0 13 PC~r/US93/041S0 typically six data components having to do with the measured correlations and a seventh component that is the projection of platform velocity onto the appropriate line-of-sight path.
The correlation components are known in the art as (1) punctual code I and Q outputs, (2) delay code I and Q outputs, and ~3) advanced code I and Q outputs, correspondina to the accumulated correlation values for in-phase and quadrature phase references and punctual, advance, and delay code references respectively.
The data vector Z_data is processed by the coupled filter signal processor 230 (Fig. 2) in accordance with the six routines identified in Fig. 3, the results of the processing being a set of commands conveyed to the reference generators 220 (Fig. 2), typically at 20-millisecond intervals, by means of the control bus 234. Each reference generator command is typically a 32-bit word comprised of values for reference signal phase and frequency and modulation phase at the beginning of the next correlation interval. The array of reference generator commands is denoted by the symbol sf.x_control.
Other inputs to the coupled-filter signal processor ~30 (Fig. 2) are the remotely-generated correction signals 236 from a differential correction station and the system initialization data 238.
The remotely-generated correction signals 236 are particularly important in compensating for uncertainties in phase of the received signals. In GPS such uncertainties are caused by deviations of the actual orbital model parameters from their nominal broadcast values, by drifts in the satellite clocks, and by unknown ionospheric and tropospheric delays. The deviations in carrier phase and frequency and modulation phase of the received signals are seen nearly identically by a fixed-site calibration station and the platform using the navigation system, and drift estimates may be broadcast to the platform where they can be removed by the coupled-filter signal processor. This mode of operation, known in the GPS art as differential GPS o.peration (DGPS) and described fully in the literature under industry standard W O 93/22691 21 1 80 1 3~ ` ` ; PC~r/US93/04150 format RTCM-104, release 2.0, (Radio Technical Commission for Maritime Services, P.O. Box 19087, Washington, D.C. 20036) is implemented in the present invention by applying the corrections directly within the signal tracking process rather than down-stream in the navigation calculations as is customary in the art.
The system initialization data 238 consists of time constants alpha_filt and process noise constants q_filt describing platform motionj parameters specifying measurement noise levels sf.r expected from the correlator system 208 (Fig. 2), and initial estimates of platform navigation states sf.part.
The outputs of the coupled-filter signal processor 230 (Fig. 2) are platform navigation data (or states) 242. Fol the case of a platform having a single antenna, these states are typically easting, northing, and altitude and the rates of change of these quantities. Additionally, the signal processor may estimate the offset between the internal system clock and the time reference of the navigation system. These outputs are of the same type as a conventional navigation system.
For the case of a platform having multiple antennas, platform orientation and motion states may be provided as outputs. Typically these are platform roll, pitch, and ya~-2~ angles with respect to the local east-north-up datum and the rates of change of these angles.
The scheduling of interactions between the correlator system 208 (Fig. 2) and the coupled-filter signal processor 230 (Fig. 2) is illustrated in Fig. 4. The figure depicts three frames of interactions: the (k-l)'th frame, the l~'th frame, and the (k+l)'th frame. Each frame begins at a time denoted by tl with an interrupt signal provided b~- the correlation system to the signal processor that signals the availabilily of data on the data bus 232 from one of the 3~ correlator channels 210, 212, 214, and 216 (Fig. 2). Tlle signal processor reads the data during the period tl-t2, performs the routines identified in Fig. 3 during the period t2-t3, and transfers the resulting command data to the WO93/22691 , PCT/US93/04150 - 2118013~
correlator system over control bus 234 (Fig. 2) during the interval t3-t4.
The command data consists of best estimates of carrier phase, carrier frequency offset, and modulating code offset to be used in adjusting the reference generator of the correlator channel in the correlator system. These estimates take into account satellite motion, DGPS biases, velocity aiding (if available), and dispersive propagation effects which are manifested by a difference in the LOS distance derived from the carrier phase measurement and the LOS distance derived from the modulating code offset measurement.
The next correlation interval begins at t4 and ends at the next tl when the next processing cycle begins.
Correlation intervals are normally loc~ed to the telemetry data bit transitions resulting in asynchrol~ous correlation intervals between correlator channels.
Specifically, during each system cycle, i.e. the time period required to service all correlator channels, the signal processor services the correlator channels on a "first come, first-served" basis. After servicing a correlator channel interrupt, the signal processor services the next correlator channel to generate an interrupt. The processing steps performed are identical for each correlator channel. ~ system cycle time for GPS is approximately 20 milliseconds.
The operations performed by the coupled-filter signal processor 230 (Fig. 2) and the routines in which they appear are as follows:
Step OPeration Routine 1 Read correl'ator data. sig. amp. est.
2 Estimate signal amplitude. sig. amp. est.
3 Recover telemetry data bit. telem. data demod.
4 Extrapolate LOS state from most filter oain recently updated channel. adjust WOg3/2269l 21 18 013 ~ pcT/uss3/o4lsn 5 Update residuals and evaluate signal filter gain quality. adjust
6 Update filter gains and covariance. filter gain adjust ~ Update LOS state vector for time t4. IQ combiner 8 Update current channel base vector IQ combiner for time t4.
9 Extrapolate LOS state vector to next IQ combiner time t4.
Extrapolate base and whole states to IQ combiner next time t4.
11 Update navigation state for next vehicle pos. &
time t4. att.
12 Find satellite positions and vehicle pos. &
velocities at time of transmission. att.
13 Estimate LOS phase and frequency for vehicle pos.
next correlation interval. att.
14 Correct LOS values with differential vehicle pos. &
correction. att.
15 Initialize new base values for phase vehicle pos. &
and frequency, this channel. att.
16 Update satellite status and vehicle pos. &
correlator channel assignments. att.
17 Compute commands for current ref. setup channel.
Each pass through the listed steps corresponds to the reading of the correlation data from one correlator channel and the writing of the commands to that correlator channel.
The first three steps performed by the coupled-filter signal processor 230 (Fig. 2) are encompassed by the telemetry data demodulator routine 244 (Fig. 3) and the signal amplitude estimator routine 246 and are concerned with reading data off the data bus, estimating signal amplitude, and reco~ering telemetry data bits. Steps 4 through 6 are performed in the filter gain adjustment routine 248. In this routine, W O 93/22691 PC~r/US93/04150 2118013 which must be performed promptly enough to enable the receiver to adapt to rapidly changing signal conditions, the statistically best weighting for measurement residuals is developed. Filter residuals are computed, signal quality is evaluated, and fresh gains are computed.
Steps 7 through 10 are performed in the IQ combiner routine 250. Unlike the conventional GPS receiver, where receiver channel measurement outputs (whole carrier phase and frequency and modulating code phase) are independently generated within each channel, without reference to line-of-sight geometery, the present invention processes the rah IQ
measurements of all channels in a statisticall~ optim~m fashion, using line-of-sight geometry and inter-path process correlations. Moreover, whenever a new set of IQ data is read 1~ from a iven correlator channel, the state estimates for all correlator channels are updated. The resulting estimates of correlator state are used to update whole estimates of LOS
quantities, and these are extrapolated to the start time for the next correlation interval. These whole values are fed to the vehicle position and attitude evaluator routine 252.
In the vehicle position and attitude evaluator routine 252 steps 11 through 16 are performed. The extrapolated LOS
values and linearized geometry are used to generate the least-squared-error estimate of the navigation states. These data may be sent to an external display or storage device, depending upon the application. Ne~t, the estimated satellit positions for signals received during the ne~t correlation interval are estimated, and the "best" estimates of transmission delay and rate of change of transmission dela~-are computed. These estimates are adjusted using differentialcorrection data, if such data is available. Finall~-, the house~eeping operations associated with scheduling the reception of signals from a newly-visible satellite alld searching for a satellite that was hidden from view for a lon 3~ period are performed.
The transmission delay and rate of change of transmission delay estimates are transformed in the reference generator command setup routine 254 (step 17) into carrier phase and W093/2269l PCT/US93/04150 2118013 ;I -frequency and modulating code phase adjustment commands for the reference generator 220 (Fig. 2) of the correlation channel being processed. In this transformation process, the ideal LOS transmission delay and delay rates are compensated for slowly changing biases that may exist between carrier phase and modulating code phase measurements. Also, the commands sent to the current correlator channel may be offset to a specified set-point.
Underlying the design of the coupled-filter signal processor 230 (Fig. 2) is the familiar state variable modeling approach for describing dynamical systems with random excitation. The structure of the model is defined in Figs.
through 9.
In the navigation platform frame of reference the quantities of interest are those comprising the navigation state vector sf.part which describe the platform's position, attitude, and motion in an earth-fixed coordinate system and those comprising the LOS state vector sf.los which describe the same quantities in terms of the distances and the rates of change of the distances of the individual platform recei~ing antennas to each of the satellites. A "state" corresponds to a particular set of values for the state variables.
The navigation state vector and the LOS state vector are defined in Fig. 5. The navigation state vector 256, in a fully configured system with multiple platform-mounted receiving antennas, represents a si~-degree-of-freedom model together with platform clock offset and offset rate resultillg in 14 state variables. Translations and rotations are defined for a body-fixed coordinate system referenced to an earth-fixed frame of reference, the latter typically being the east-north-up system. The navigation state drives the LOS state and a locally linear transformation is used to carry small changes observed in the navigation state into the LOS state and vice versa.
In the case of four satellites (or reference transmitters) and a platform having a single recei~ing antenna, the LOS state vector 258 consists of the four vector components 260 where the index i identifies a particular `~ 2118013~
satellite-to-platform path, di(t) is the line-of-sight distance to the i'th satellite at time t, ~c is the nominal carrier transmit frequency, C is the speed of light, d~(t) is the rate of change of di(t), carrier bias is the offset between the observed modulating code and the carrier phase, and carrier bias rate is the rate of change of carrier bias.
The quantity iq.phase_off[i] denotes the total path length of the i'th path at time t e~pressed in terms of phase of the carrier. The quantity iq.omega_off[i] denotes the rate of change of the path length of the i'th path. The last t~Yo terms of each vector component 260 are included to model biases between the free-space line-of-sight path length and the path length measured by accumulated carrier phase. In general, a carrier phase observation will be offset from the true path length due to the carrier cycle ambiguity as well as dispersive propagation phenomena caused by the ionosphere.
These slowly-varying states are estimated at the same time as the free-space quantities.
The LOS state vector X(k) 262 is e~pressèd as a sum of two parts, the base value Xbase(k) and the incremental part X(k). The base value is the value sent to a correlator channel at the start of its new correlation interval. The incremental value is the difference between the whole LOS
state and the specified base value. The LOS state vector 26~, in the case of four satellites and a single receiving antenna, consists of the four vector components 266 where the inde~; i identifies a particular satellite-platform path. Each of the vector components 266 is comprised of the free-space terms 260 together with the velocity aiding bias term iq.aid_bias[i].
As previously discussed, the signal processor accommodates ~;elocity aiding inputs from the velocit~- aiding system 23~
(Fig. 2). Aiding devices can have large biases ~hich are detrimental to system operation if uncompensated. ~'elocitv aiding input biases are estimated and compensated for in the IQ combiner routine 250 (Fig. 3).
The base LOS state evolves with discrete time in accordance with equation 268. The base LOS state vector 270, in the case of four satellites and a single receiving antenna, WO93/22691 ~ ; PCT/US93/04150 2118013 i`
consists of the four vector components 272. Each of the vector components 266 is comprised of the free-space terms 260 together with the zero velocity aiding bias term iq.aid_bias[i]. The state transition matrix ~ic~ will be discussed below in connection with Fig. 6.
The incremental LOS state vector 276, shown for the case of n line-of-sight paths, consists of the five vector components. Equation 278 follows from equations 262, 26~, 266, 268, 270, and 272. The quantity W(k) models the change on each LOS path occurring during one correlation interval.
These changes, observed as changes in carrier phase, carrier frequency, carrier bias, carrier bias rate, and aid bias, are described by the covariance matrix se.qe 281. The structure of se.qe plays a central role in coupling of the trac~ing channels and is defined in detail in Figs. 20, 21, and 22, to be discussed below.
The state transition matrix appearing in equations 268 and 274 is defined by expressions 282 and 284. Within each partition of the transition matrix are three constants. The first of these, Tc, is simply the correlation interval duration. The carrier bias time constant acarbias is deri~;ed from the process correlation time of the assumed carrier bias rate phenomenon, and the aid bias time constant aaidbia~ is established from a specified model for the drift rates of the external velocity sensor which is part of the velocity aiding system 235 (Fig. 2).
The structure of the observation model is given in Figs.
9 Extrapolate LOS state vector to next IQ combiner time t4.
Extrapolate base and whole states to IQ combiner next time t4.
11 Update navigation state for next vehicle pos. &
time t4. att.
12 Find satellite positions and vehicle pos. &
velocities at time of transmission. att.
13 Estimate LOS phase and frequency for vehicle pos.
next correlation interval. att.
14 Correct LOS values with differential vehicle pos. &
correction. att.
15 Initialize new base values for phase vehicle pos. &
and frequency, this channel. att.
16 Update satellite status and vehicle pos. &
correlator channel assignments. att.
17 Compute commands for current ref. setup channel.
Each pass through the listed steps corresponds to the reading of the correlation data from one correlator channel and the writing of the commands to that correlator channel.
The first three steps performed by the coupled-filter signal processor 230 (Fig. 2) are encompassed by the telemetry data demodulator routine 244 (Fig. 3) and the signal amplitude estimator routine 246 and are concerned with reading data off the data bus, estimating signal amplitude, and reco~ering telemetry data bits. Steps 4 through 6 are performed in the filter gain adjustment routine 248. In this routine, W O 93/22691 PC~r/US93/04150 2118013 which must be performed promptly enough to enable the receiver to adapt to rapidly changing signal conditions, the statistically best weighting for measurement residuals is developed. Filter residuals are computed, signal quality is evaluated, and fresh gains are computed.
Steps 7 through 10 are performed in the IQ combiner routine 250. Unlike the conventional GPS receiver, where receiver channel measurement outputs (whole carrier phase and frequency and modulating code phase) are independently generated within each channel, without reference to line-of-sight geometery, the present invention processes the rah IQ
measurements of all channels in a statisticall~ optim~m fashion, using line-of-sight geometry and inter-path process correlations. Moreover, whenever a new set of IQ data is read 1~ from a iven correlator channel, the state estimates for all correlator channels are updated. The resulting estimates of correlator state are used to update whole estimates of LOS
quantities, and these are extrapolated to the start time for the next correlation interval. These whole values are fed to the vehicle position and attitude evaluator routine 252.
In the vehicle position and attitude evaluator routine 252 steps 11 through 16 are performed. The extrapolated LOS
values and linearized geometry are used to generate the least-squared-error estimate of the navigation states. These data may be sent to an external display or storage device, depending upon the application. Ne~t, the estimated satellit positions for signals received during the ne~t correlation interval are estimated, and the "best" estimates of transmission delay and rate of change of transmission dela~-are computed. These estimates are adjusted using differentialcorrection data, if such data is available. Finall~-, the house~eeping operations associated with scheduling the reception of signals from a newly-visible satellite alld searching for a satellite that was hidden from view for a lon 3~ period are performed.
The transmission delay and rate of change of transmission delay estimates are transformed in the reference generator command setup routine 254 (step 17) into carrier phase and W093/2269l PCT/US93/04150 2118013 ;I -frequency and modulating code phase adjustment commands for the reference generator 220 (Fig. 2) of the correlation channel being processed. In this transformation process, the ideal LOS transmission delay and delay rates are compensated for slowly changing biases that may exist between carrier phase and modulating code phase measurements. Also, the commands sent to the current correlator channel may be offset to a specified set-point.
Underlying the design of the coupled-filter signal processor 230 (Fig. 2) is the familiar state variable modeling approach for describing dynamical systems with random excitation. The structure of the model is defined in Figs.
through 9.
In the navigation platform frame of reference the quantities of interest are those comprising the navigation state vector sf.part which describe the platform's position, attitude, and motion in an earth-fixed coordinate system and those comprising the LOS state vector sf.los which describe the same quantities in terms of the distances and the rates of change of the distances of the individual platform recei~ing antennas to each of the satellites. A "state" corresponds to a particular set of values for the state variables.
The navigation state vector and the LOS state vector are defined in Fig. 5. The navigation state vector 256, in a fully configured system with multiple platform-mounted receiving antennas, represents a si~-degree-of-freedom model together with platform clock offset and offset rate resultillg in 14 state variables. Translations and rotations are defined for a body-fixed coordinate system referenced to an earth-fixed frame of reference, the latter typically being the east-north-up system. The navigation state drives the LOS state and a locally linear transformation is used to carry small changes observed in the navigation state into the LOS state and vice versa.
In the case of four satellites (or reference transmitters) and a platform having a single recei~ing antenna, the LOS state vector 258 consists of the four vector components 260 where the index i identifies a particular `~ 2118013~
satellite-to-platform path, di(t) is the line-of-sight distance to the i'th satellite at time t, ~c is the nominal carrier transmit frequency, C is the speed of light, d~(t) is the rate of change of di(t), carrier bias is the offset between the observed modulating code and the carrier phase, and carrier bias rate is the rate of change of carrier bias.
The quantity iq.phase_off[i] denotes the total path length of the i'th path at time t e~pressed in terms of phase of the carrier. The quantity iq.omega_off[i] denotes the rate of change of the path length of the i'th path. The last t~Yo terms of each vector component 260 are included to model biases between the free-space line-of-sight path length and the path length measured by accumulated carrier phase. In general, a carrier phase observation will be offset from the true path length due to the carrier cycle ambiguity as well as dispersive propagation phenomena caused by the ionosphere.
These slowly-varying states are estimated at the same time as the free-space quantities.
The LOS state vector X(k) 262 is e~pressèd as a sum of two parts, the base value Xbase(k) and the incremental part X(k). The base value is the value sent to a correlator channel at the start of its new correlation interval. The incremental value is the difference between the whole LOS
state and the specified base value. The LOS state vector 26~, in the case of four satellites and a single receiving antenna, consists of the four vector components 266 where the inde~; i identifies a particular satellite-platform path. Each of the vector components 266 is comprised of the free-space terms 260 together with the velocity aiding bias term iq.aid_bias[i].
As previously discussed, the signal processor accommodates ~;elocity aiding inputs from the velocit~- aiding system 23~
(Fig. 2). Aiding devices can have large biases ~hich are detrimental to system operation if uncompensated. ~'elocitv aiding input biases are estimated and compensated for in the IQ combiner routine 250 (Fig. 3).
The base LOS state evolves with discrete time in accordance with equation 268. The base LOS state vector 270, in the case of four satellites and a single receiving antenna, WO93/22691 ~ ; PCT/US93/04150 2118013 i`
consists of the four vector components 272. Each of the vector components 266 is comprised of the free-space terms 260 together with the zero velocity aiding bias term iq.aid_bias[i]. The state transition matrix ~ic~ will be discussed below in connection with Fig. 6.
The incremental LOS state vector 276, shown for the case of n line-of-sight paths, consists of the five vector components. Equation 278 follows from equations 262, 26~, 266, 268, 270, and 272. The quantity W(k) models the change on each LOS path occurring during one correlation interval.
These changes, observed as changes in carrier phase, carrier frequency, carrier bias, carrier bias rate, and aid bias, are described by the covariance matrix se.qe 281. The structure of se.qe plays a central role in coupling of the trac~ing channels and is defined in detail in Figs. 20, 21, and 22, to be discussed below.
The state transition matrix appearing in equations 268 and 274 is defined by expressions 282 and 284. Within each partition of the transition matrix are three constants. The first of these, Tc, is simply the correlation interval duration. The carrier bias time constant acarbias is deri~;ed from the process correlation time of the assumed carrier bias rate phenomenon, and the aid bias time constant aaidbia~ is established from a specified model for the drift rates of the external velocity sensor which is part of the velocity aiding system 235 (Fig. 2).
The structure of the observation model is given in Figs.
7 through 9. In the material that follows it is assumed that four correlator channels are tracking four signals that were propagated over four LOS paths.
The input data vector to the coupled-filter si,nal processor 230 (Fig. 2) is defined by the expression 286 in Fig. 7. Each of the four vector components 288 consists of seven components. The first six are outputs from the associated correlator 218 (Fig. 2) and the seventh is the output from the velocity aiding system. The seven components are read from the data bus 232 (Figs. 2 and 3) b~- the signal processor 230.
WO93/22691 2 I I 8 0 I 3 PCT/USg3/04150 The data vector can be modeled as a nonlinear function 290 of the incremental LOS state vector X. This dependenc~-on only the incremental LOS state results from the signal processor command sent back to each channel at the start of its correlation interval which corrects the base value of the LOS state.
The quantity h~ 292 consists of four vector components 294. The seven components of each of the vector components 294 are defined in Figs. 8 and 9. The quantity ~k is the noise voltage that contaminates the correlation measurements.
The correlation measurement noise covariance matri~ 296 is initialized with values obtained from the hardware specification of the correlator system 208 (Fig. 2).
The parameters sf.setpoint_phase, sf.setpoint_freq, and sf.setpoint_code that appear in the equations shown in Figs.
The input data vector to the coupled-filter si,nal processor 230 (Fig. 2) is defined by the expression 286 in Fig. 7. Each of the four vector components 288 consists of seven components. The first six are outputs from the associated correlator 218 (Fig. 2) and the seventh is the output from the velocity aiding system. The seven components are read from the data bus 232 (Figs. 2 and 3) b~- the signal processor 230.
WO93/22691 2 I I 8 0 I 3 PCT/USg3/04150 The data vector can be modeled as a nonlinear function 290 of the incremental LOS state vector X. This dependenc~-on only the incremental LOS state results from the signal processor command sent back to each channel at the start of its correlation interval which corrects the base value of the LOS state.
The quantity h~ 292 consists of four vector components 294. The seven components of each of the vector components 294 are defined in Figs. 8 and 9. The quantity ~k is the noise voltage that contaminates the correlation measurements.
The correlation measurement noise covariance matri~ 296 is initialized with values obtained from the hardware specification of the correlator system 208 (Fig. 2).
The parameters sf.setpoint_phase, sf.setpoint_freq, and sf.setpoint_code that appear in the equations shown in Figs.
8 and 9 are the desired offsets between the reference signal phase, frequency, and code values and the actual received values at the start of correlation interval k+1. R() is the modulating code autocorrelation function, and _ denotes the offset delay between the on-time reference code of the GPS
correlator and its advanced and delayed code replicas. The quantity iq.am is the peak received signal amplitude. Time tag value t4[k-l,i] denotes the start of the k'th correlation interval, channel i. Time tag value tl[k,i] denotes the end time of the k'th correlation interval, channel i.
The velocity aiding input 298 (Fig. 9~, measured with respect to the base value of the LOS rate iq.omega_base, consists of the four vector components 300. The quantity aid_bias[i] is the bias in the velocity value supplied by the velocity aiding system 235 (Fig. 2).
The signal processing routines identified in Fig. 3 will now be described in the chronological order of performance.
The signal amplitude estimator routine 246 (Fig. 3) causes the coupled-filter signal processor 230 (Fig. 2) to 3~ read the raw correlator data vector Z_data and to produce estimates of signal amplitude iq.am_hat. In order for the signal processor to properly adapt to changing signal quality conditions as observed at the platform receiving antenna(s), W O 93/22691 PC~r/US93/04150 it is necessary that signal amplitude estimates be updated at a rapid rate, typically every 10 milliseconds. The two - processing steps are detailed in Fig. 10.
The first step 302 consists of reading the most recent data values from the correlator system 208 (Fig. 2) and the velocity aiding system 235 (Fig. 2) from the data bus 232 (Figs. 2 and 3).
The second step 304 consists of computing an updated amplitude estimate for each satellite line-of-sight path. To this end the I and Q values for the on-time code reference correlator are scaled, squared, and summed and the square root of the result is extracted as indicated by equation 306. The data elements being processed are identified by the expressions 308. The current signal processor estimate of the frequency state for the i'th correlator channel iq.omeae_h[i]
is obtained from the estimate sf.x_current_last of the incremental LOS state vector 278 (Fig. 6) for the channel currently being processed and referenced to time t~ of that channel. The estimate sf.x_current_last is made in the IQ
combiner routine 250 (Fig. 3).
The amplitude estimates are updated by passing the lambda values through a digital "infinite inpulse response" smoothing filter. The updated amplitude estimate for the i'th channel at cycle n is given by the equation 310.
In step 3 the signal processor performs the telemetry data demodulator routine 244 (Fig. 3) and produces one or more telemetry data bits. The operations performed are indicated in Fig. 11. In the case of phase-shift-keyed (PSK~ telemetry data modulation, the punctual I data components of the Z_data vector 286 (Fig. 7) are tested for sign. A plus sign is mapped into a binary one, a minus sign is mapped into a binary zero, and the raw Z_data vector 286 is then stripped of the data telemetry bits by inverting the sign of all component~ of a given correlator channel for which a binary zero ~as declared, prior to further processing of Z_data. !he recovered data telemetry bits, which for GPS are received at 50 baud, are collected in buffers and parity is checked. The data bits are tested for word and frame synchronization in a W O 93/22691 2 1 1 8 0 1 3 ~ PC~r/US93/04150 con~entional fashion as described for GPS in Interface Control Document ICD-200 tavailable from the NAVSTAR GPS Joint Program Office, United States Air Force Systems Command). .~lso, as described in this same document, the telemetry data is then used to compute satellite positions and drift correction parameters associated with satellite clock drift.
Periodically, the start and stop intervals t4 and tl of each correlator channel will need adjustment so as to remain in approximate alignment with the edges of the PS~ telemetry data. A standard bit synchronizer based upon the concept of the wide-sense cyclostationary process is incorporated for this purpose. (See Moghazy, Maral, and Blanchard, "Digilal PCM Bit Synchronizer and Detector", IEEE Transactions on Communication Systems, Vol. Com-28, August 1980, pp. 119l-1200.) In the present embodiment this synchronizationadjustment function is implemented by squaring values of the "on-time" I data stored in 288, and averaging these squared values over typically 100 correlation intervals to obtain a control signal which is responsive to the error in timing of correlation interval start relative to the time of data bit arrival at the receiver. This control signal is then used to periodically adjust the correlation interval start.
The Z_data vector produced by the correlator systenl 208 (Fig. 2) is used in the IQ combiner routine 250 (Fig. 3) to update the estimate sf.x of the true incremental LOS state vector _X 278 (Fig. 6). This key step can be performed using any one of a number of standard techniques, among them being the extended l~alman filter. A central part of any such method is the computation of the filter gain matrix which carries the measurement residuals into updated state estimates.
The signal processor accomplishes this function by means of the filter gain adjust routine 248 (Fig. 3~. It takes as input the signal amplitude estimates iq.am_hat 310 (Fig. 10), the nominal correlation measurement noise covariance matri~
3~ sf.r 296 (Fig. 7) that is user-supplied 238 (Fig. 3), and the path state noise covariance matrix se.qe 280 (Fig. 6) and produces the filter gain matrix sf.k. The detailed steps are shown in Figs. 12 through 16.
WO93/22691 21 18 0 i3 PCT/US93/04150 As discussed earlier, the signal processor services correlator channels on a "first come, first served" basis.
Before computing the residual errors in the LOS state estimates for the channel currently being served, the "best"
estimate of the LOS state based on the previously-served channel at its t4 time must be extrapolated to the t4 time of the current channel. It is important to note that the present invention, unlike conventional navigation receivers, performs an update of the LOS states of all correlator channels whenever fresh data becomes available from anv of its correlator channels.
The extrapolation of the LOS state from the most recently updated correlator channel (step 4) is described in Fig. 12 by equations 312. The vector sf.x_previous_last is the estimate of the LOS state following the most recently processed channel referenced to time t4 for that channel. The vector sf.x_current_last is the estimate of the LOS state for the channel currently being processed referenced to time t4 for that channel. The state transition matrix ~los was discussed earlier in connection with Fig. 6.
The filter state error covariance matrix sf.p is a product of the standard Kalman filter methodology. The error covariance matrix based on the previously-serviced channel at its t4 time is also extrapolated in step 4 to the t4 time of the current channel. The extrapolation process is specified in Fig. 12 by equation 314. The quantity sf.p_previous is the error covariance matrix at time t4[k-l,i-1] following the update of the LOS state at the last t4 time for the previous channel. The quantity sf.p_current is the error covariance matrix at time t4[k-l,i]. The error covariance matrix is extrapolated with the same time step as used in extrapolatin~
the LOS state 312. The quantity se.qe 280 (Fig. 6) was defined earlier.
Step 5 is detailed in Figs. 13, 14, and 15. The predicted data sf.z_hat 316 from the correlator system 208 (Fig. 2) and the velocity aiding system 23~ (Fig. 2) are subtracted from the current raw data Z_data to form the measurement residual vector sf.res 318 as indicated in Fig.
WO93/22691 2 1 18 ~ 13 PCT/US93/04150 13. At this time the seven residual components associated with the current channel are examined for reasonableness. A
channel with residuals above a rejection threshold is noted and its "bad residual count" is incremented.
There are two responses to a residual above threshold corresponding to whether the bad residual count has exceeded a limit or not:
(a) If the bad residual count has exceeded the limit, the current measurement is rejected which is most easily accomplished by increasing sf.r to sf.r_high, a quantity that is set to a large value, say 100,000. Then, this channel is removed from the list of active channels and is processed separately in a decoupled mode until residuals on this channel improve;
(b) If the bad residual count has not exceeded the limit, the current measurement is rejected as in the above case, but the channel is kept in the coupled mode.
The seven predicted data values 316 for the current LOS
path are defined by the equations 320 in Fig. 13 and all of the equations in Fig. 14. The quantities appearing in the equations are obtained from the estimate of the LOS state for the current channel sf.x_current_last 312 (Fig. 12) as shown in Fig. 15. The definition of the LOS state sf.x shown in Fig. 15 follows that given previously for the true LOS state 276, 278 (Fig.6).
The correlator system 208 (Fig. 2) is highly non-linear in the operating region of interest and is linearized about the current estimate of the LOS state vector in step 6.
kinearization is accomplished by ta~ing the gradient of the analytical expressions for the I and Q correlator outputs ~.-ith respect to the full set of LOS state variables. There are ~Ip to five LOS state variables defined for each line-of-sight path: phase and frequency error, modulating code offset and offset rate, and LOS velocity aiding input bias. For the i'th channel, the array of derivatives in the gradient are evaluated at sf.x_current_last and stored in array sf.h_i.
The gradient terms for the four-channel correlator system WO93/22691 2 1 1 8 0 1 3 ~ ` PCT/US93/~150 208 (Fig. 2) are defined by expressions 324, 326, and 328 in Fig. 16. The array 324 has dimensions 7x20 with three null partitions of dimensions 7x5.
The equation 330 for the Kalman filter gain sf.k_[i] for the i'th channel utilizes the current value for the correlator measurement noise covariance matrix, either sf.r_[i] if the residual is normal or sf.r_high_[i] if the residual is above the reasonableness threshold. The values in sf.r_[i] and sf.r_high_[i] are supplied as part of the initialization data 238 (Fig. 3). The overall LOS state covariance matri~
extrapolated to time t4[k-l,i], sf.p_current, in accordance with the standard Kalman filter design method.
The quantity sf.p_current is the filter state error covariance matrix extrapolated from the t4 time of the previously-updated correlator channel. The quantity sf.r_[i]
is the correlation measurement noise covariance matrix for the batch of seven measurement components for the current channel.
The dimensions of the arrays are 20x7 for sf.~_[i], 20x20 for sf.p_current, 7x20 with three zero partitions for sf.h_[i], and 7x7 for sf.r_[i].
Following the filter gain computation, the filter state error covariance matrix sf.p is updated 332 using the standard Kalman filter methodology. The updating reflects the improved estimation error for all correlator states following the processing of the current batch of correlator I and Q data together with the velocity aiding data. The quantity I in equation 332 is the identity matrix having dimensions 20x20.
It should be noted that the previously-calculated value of sf.p_previous is overwritten with the newly-calculated value.
It is also important to note that the new filter state error covariance matrix reflects improvement not onlv in the estimate for the LOS state for the current correlator channel but also in the estimates for the LOS states of all channels.
This result is in marked contrast to the traditional prior-art decoupled signal processing scheme.
The IQ combiner routine 250 (Fig. 3) shown in Figs. 1~
and 18 utilizes as inputs the residuals for the current channel sf.res_[i] and the filter gain matrix for the current WO93~22691 2 1 1 8 0 1 3 PCT/US93/041~0 channel sf.k_[i]. The output that results from the performance of the routine is the estimate of the LOS state sf.x for all LOS paths.
In step 7 the LOS state vector estimate sf.x (Fig. 15) referenced to the t4 time of the current channel is updated by multiplying the measurement residual vector sf.res by the filter gain matrix sf.k. in accordance with equation 334 shown in Fig. l7.
Following the updating of the LOS state vector estimate in step 7, the components of the LOS state base vector sf.base_last associated with the current channel are updated in step 8 by adding the updated LOS state values for the current channel to the LOS state base values for the current channel in accordance with equation 336. The term sf.x_base_last denotes the base values employed in each correlator channel for the k-l system cycle, the cycle that has just ended.
The LOS state values for the current channel, except for the velocity-aiding bias state, are then reset to zero in accordance with equation 338.
In step 9 the LOS state vector is next extrapolated to the start time t4[k,i] of the next correlation interval for the current channel, and the current channel is then defined to be the previous channel in accordance with equations 340 shown in Fig. 18. The vector sf.x_previous_next is the LOS
state estimate following the most recently completed channel referenced to the next t4 time for that channel. The vector sf.x_current_next is the LOS state estimate referenced to the next t4 time for the current channel.
In step l0 the base LOS values for the correlator system are then extrapolated to the current channel's ne~t correlation start t4[k,i] using the appropriate time step for each correlator channel in accordance with equations 3~2. The term sf.x_base_temp denotes extrapolations of sf.x_base to arbitrary time tags in the ~'th cycle.
The predicted whole values sf.x_whole_current for LOS
state estimates referenced to time t4[i,k] are formed b~-adding the updated and extrapolated LOS state vector to the W O 93/22691 PC~r/US93/04150 _base_temp values for each channel as indicated by equations 344.
The vehicle position and attitude evaluator routine 252 (Fig. 3) takes as initialization input: (1) the user-supplied initial value for platform position, velocity, attitude, and attitude rate sf.part; (2) the user-supplied initial values for the correlation measurement noise covariance matrix sf.r;
(3) the user-supplied values for platform velocity and platform rotation rate random process correlation times alpha_filt; and (4) the user-supplied values for platform velocity and platform rotation rate random process state noise disturbance power q_filt.
Quantities generated by the platform position and attitude evaluator routine 252 (Fig. 3) and utilized by other routines are: (1) the computed path state noise covariance matrix se.qe for the currently-defined satellite-to-platform antenna paths; and (2) the estimates of LOS whole phase and frequency values sf.x_base.
The details of the routine are shown in Figs. 19 through 24.
In step 11 the linearization of the LOS path oeometry is accomplished. The linearization relates small changes in platform navigation states (i.e. translation and rotation of the body-fixed coordinate system) to changes in the Euclidean distances and rates of change of distances measured between each satellite and the platform antenna(s). The linearization is found by taking the gradient of the set of distance vectors with respect to the navigation state vector and evaluating the resulting gradient matrix functions for the most recent estimate of the platform state vector se.state.
The line-of-sight linearization matrix 346 and the derivatives to be evaluated 348 are shown in Fig. 19.
As a prelude to the computation of the latest estimate of the platform na~igation vector, the offset vector iq.delt_x_los is obtained by subtracting the previous values of whole carrier phase and frequency from the current values for all channels, as indicated by equations 350 and 352.
Next, the pseudo-inverse matrix se.temp3 that relates W093/2269l 2118 013 - ~ ` PCT/US93/04150 _ changes in the line-of-sight quantities to the changes in navigation quantities is computed in accordance with equation 354. As is well known in the literature (see Gelb, APPlied OPtimal Estimation, MIT Press, 1974), the pseudo-inverse provides an optimal estimate in the least-squared-error sense for an overdetermined linear system of equations. This principle provides the basis for equation 356 which relates the current value sf.part_current of the platform navigation state to the previous value sf.part_previous. In accordance with equation 358, sf.part_current becomes sf.part_previous in anticipation of the next updating of sf.part.
As mentioned earlier, the central theme of the invention is the proper coupling of the signal tracking functions between all satellites being received. The quantity which permits this statistically optimum approach is the path state noise covariance matrix se.qe. The calculation of se.qe is accomplished by matrix operations involving the line-of-sight linearization matrix 346 and the platform motion state model noise covariance matrix se.qs. As wil be shown, the quantity se.qs is computed from the quantities alpha_filt and q_filt which are part of the user-supplied system initialization data 238 (Fig. 3).
The overall structure of se.qe 360 for n LOS paths is shown in Fig. 20. The QL terms contained in se.qe 360 describe the correlation in LOS paths due to user translations, rotations, and receiver clock bias and are derived from se.qe_los 362 which in turn is computed by means of equation 364 using the LOS linearization matrix previously discussed. The platform motion state model noise covariance matrix se.qs 366 is generated from input parameters as shown in Fig. 21. As indicated in the figure, the platform motion state noise covariance matrix 366 is partitioned into smaller arrays corresponding to the translation, cloc~ drift, and rotational motions of the platform. By way of an example, the 3~ partition 368 corresponding to the x-axis translation is shown. The elements of this 2-by-2 partition are computed from equations 370 where the various constants are obtained externally as components of the vectors alpha_filt 378 and WO93/22691 2 1 18 0 13 PCT/US93/~150 q_filt 380, discussed below in connection with Fig. 22. This computation represents the transformation between the continuous-time navigation model and the discrete-time navigation model.
The identical partitions QB of se.qe 360 (Fig. 20) are next computed. The QB partitions need only be computed once -at receiver initialization time. The 2-by-2 QB array 372 is shown in more detail in Fig. 22. The continuous-time state noise model describing the-carrier-code bias process is read in as components of the vectors alph_filt 378 and q_filt 380.
Then, the operations 374 needed to convert from the continuous-time model to the discrete-time model are performed resulting in the elements of QB.
Finally, the identical partitions QR of se.qe 360 (Fig.
l~ 20) are computed. The QR partitions also need onl~- be computed once - at receiver initialization time. The l-by-l QR partition 376 is defined in Fig. 22.
The computation of the new estimates of LOS whole phase and frequency values sf.x_base for the current channel requires the satellite positions and velocities to be updated.
Satellite position and velocity updating is accomplished conventionally in step 12 as described for GPS in Interface Control Document ICD-200 (available from the NAVSTAR GPS Joint Program Office, United States Air Force Systems Command).
First, the satellite location at its time of transmission corresponding to the beginning time t4[k,i] of the next correlation for the current channel is established. .~ signal observed at the platform at t4[k,i] will have been emitted at the satellite at an earlier time. The current estimate iq.phase_hat[i] of the line-of-sight distance from the satellite to the platform, converted from radian units to time units, is the appropriate time advance. The estimated time of transmission is computed using equation 382 shown in Fig. 23.
The quantity iq.phase_hat is obtained from sf.x_whole_current 3~ 344 (Fig. 18).
Next, the satellite position at the estimated time of transmission must be computed. For GPS the standard ICD-200 computation is employed to find the (x,y,z) earth-centered wo 93/226gl 2 1 18 0 13 PCT/US93/~150 . _ earth-fixed (ECEF) Cartesian coordinates corresponding to position ephem_position[i] 384 and to velocity - ephem_velocity[i] 386, the value of i corresponding to the satellite associated with the current channel.
The ECEF values obtained above are used in step 13 to compute the "best" ~stimates of LOS phase, ephem_LOS[i] 388, and frequency, ephem_LOS_dot[i] 390, for the current channel.
These values will be used for the next correlation interval.
If the externally-supplied differential corrections from a DGPS reference station are available, they are applied in step 14 to the above LOS values. The LOS and LOS rate values are also compensated for the user clock offset and offset rate, as previously estimated and stored in elements T and T_dot of the current navigation state vector sf.part_current.
1~
There are two possible system configurations for utilizing externally-generated differential corrections in obtaining solutions to the platform navigation problem. In one configuration the differential corrections for pseudorange and pseudorange-rate are generated in a conventional fashion at a reference station geographically separated from the platform using techniques and standards given in the ~GPS
RTCM-104 document and sent by radio link or other real-time communication medium to the platform where they are passed on to the coupled-filter signal processor 230 (Fig. 2) and utilized in the vehicle position and attitude evaluator routine 252 (Fig. 3).
For n LOS paths the psuedorange and pseudorange-rate differential corrections are defined by the expressions 392 and 394 respectively in Fig. 24. Although these differential corrections are generated by well-known techniques, their application directly in the signal tracking equations is unique to the present invention. By applying corrections in this fashion, the effects of satellite biases and propagation anomalies are removed from the phase and frequency estimates that are the products of the signal tracking equations. If this correction process were not performed at this stage, interpath line-of-sight correlations would be reduced thereby WO93/22691 2 1 18 0 1 3 ` PCT/US93/04150 causing the coupled-filter signal processor to be less effective.
There are certain situations where a real-time radio link between reference station and platform cannot be maintained, e.g. land surveying.
The technique employed for non-real-time operation is quite different from that employed for conventional DGPS.
Whereas in the conventional scheme pseudorange and carrier phase outputs from the reference station and platform are recorded for later merging in a signal processing scheme, the present invention requires that the sampled in-phase and quadrature outputs from the conventional receiver front-end (that form the inputs 206 (Fig. 2) to the correlation system 208 (Fig. 2)) be recorded by a high-speed digital recorder.
These recorded digital data streams together with recorded psueudorange and carrier phase outputs obtained from the reference station are then played at some convenient later time into the correlator system 208 (Fig. 2) and the coupled-filter signal processor 230 (Fig. 2). The computational process performed by the coupled-filter signal processor is the same as it was for the real-time situation. In this way the navigational history of a platform can be determined ir~ a manner fully-equivalent to the real-time process.
For either real-time or non-real-time operation the corrections are introduced as specified by equations 396 of Fig. 24. It should be noted that ephem_LOS_hat is expressed in radian units.
Finally, in step l5,the above corrected and compensated estimates are copied 398 into the first two components of sf.base. These base values, as well as the carrier bias values, replace the corresponding existing components in sf.x_base_next. The term sf.x_base_next denotes the base values to be employed in the h'th system cycle, the (k-l)'th cycle having just ended.
In step 16 satellite almanac data received over satellite telemetry links are used to evaluate on a regular basis the elevation angles of satellites currently being processed, as in conventional GPS receivers. The potential visibility of other satellites in the constellation not being tracked are also regularly checked using the same constellation almanac data. Satellites that are presently being tracked are taken out of the tracking mode .hen they fall below an elevation masking angle of typically five degrees. Fresh satellites are brought into the signal acquisition process if their predicted elevation angles are above the elevation masking angle. Once it is determined that a new satellite is to be acquired and tracked, an unused correlator channel is assigned to that satellite, and the signal acquisition process commences. Signal acquisition is essentially the same as for a conventional signal processor, as described in many references: e.g. M.~. Simon et al., Spread SPectrum Communications, Vol. III, Computer Science Press, Rockville, MD, 1985; C.R. Cahn and E.H. Martin, "Design Considerations for a Spread Spectrum Navigation Receiver,"
Proceedings of International Conference on Electronic Systems and Navigation Aids, Paris, France, November 1977. Note that during the initial acquisition process the coupled-filter signal processor 230 (Fig. 2) and the assigned correlator channel operate in a decoupled mode. Once code and frequency acquisition have been achieved, the correlator channel is controlled in the coupled mode.
In performing step 17 in the reference generator command setup routine 254 (Fig. 3), the signal processor 230 (Fig. 2) utilizes the data contained in sf.x_base_current for the current channel and generates commands to be sent to the reference generator 220 (Fig. 2) associated with the current channel which will cause the reference generator to generate a reference signal having specified values for carrier phase and frequency and modulating code phase during the next correlation interval. The reference generator command setup routine is detailed in Fig. 25.
Because of setup resolution limitations inherent in the 3~ design of the correlator system 208 (Fig. 2), the signal processor 230 (Fig.2) must provide for error setpoints which are non-zero. For example, in the case of GPS the code generator in the reference generator 220 (Fig. 2~ will only be WO93/22691 2 1 18 b i3 PCT/US93/04150 programmable in fractional code chip steps. Likewise, the phase and frequency of the signal produced by a numerically-controlled oscillator (NCO) in a typical GPS reference generator can only be adjusted in increments that are much coarser than the precision of the values maintained by the signal processor. In recognition of this situation, and to avoid any operating limitation postulated on zero LOS state values, provisions are made in the signal processor for non-zero set-points to be sent to the reference generators.
The commands sf.x_control to the reference generators are computed in accordance with equation 400. The base values come from the sf.x_base_next array (see equations 394 (Fig.
24) and discussion pertaining thereto).
The elements in the control vector correspond to NCO
phase and frequency setup and code generator setup for one correlator channel. The first component, which controls the NCO phase setup, is the LOS estimate iq.phase_base plus the carrier bias minus the desired NCO phase or setpoint. The second component, which controls the NCO frequency setup, is the estimated LOS rate plus the estimated carrier bias rate minus the desired NCO frequency offset or setpoint. The third component, which controls the code generator, is simply the LOS estimate iq.phase_base minus the desired code generator offset or setpoint.
The specified setpoints 402 are stored in read-only memory in the signal processor.
correlator and its advanced and delayed code replicas. The quantity iq.am is the peak received signal amplitude. Time tag value t4[k-l,i] denotes the start of the k'th correlation interval, channel i. Time tag value tl[k,i] denotes the end time of the k'th correlation interval, channel i.
The velocity aiding input 298 (Fig. 9~, measured with respect to the base value of the LOS rate iq.omega_base, consists of the four vector components 300. The quantity aid_bias[i] is the bias in the velocity value supplied by the velocity aiding system 235 (Fig. 2).
The signal processing routines identified in Fig. 3 will now be described in the chronological order of performance.
The signal amplitude estimator routine 246 (Fig. 3) causes the coupled-filter signal processor 230 (Fig. 2) to 3~ read the raw correlator data vector Z_data and to produce estimates of signal amplitude iq.am_hat. In order for the signal processor to properly adapt to changing signal quality conditions as observed at the platform receiving antenna(s), W O 93/22691 PC~r/US93/04150 it is necessary that signal amplitude estimates be updated at a rapid rate, typically every 10 milliseconds. The two - processing steps are detailed in Fig. 10.
The first step 302 consists of reading the most recent data values from the correlator system 208 (Fig. 2) and the velocity aiding system 235 (Fig. 2) from the data bus 232 (Figs. 2 and 3).
The second step 304 consists of computing an updated amplitude estimate for each satellite line-of-sight path. To this end the I and Q values for the on-time code reference correlator are scaled, squared, and summed and the square root of the result is extracted as indicated by equation 306. The data elements being processed are identified by the expressions 308. The current signal processor estimate of the frequency state for the i'th correlator channel iq.omeae_h[i]
is obtained from the estimate sf.x_current_last of the incremental LOS state vector 278 (Fig. 6) for the channel currently being processed and referenced to time t~ of that channel. The estimate sf.x_current_last is made in the IQ
combiner routine 250 (Fig. 3).
The amplitude estimates are updated by passing the lambda values through a digital "infinite inpulse response" smoothing filter. The updated amplitude estimate for the i'th channel at cycle n is given by the equation 310.
In step 3 the signal processor performs the telemetry data demodulator routine 244 (Fig. 3) and produces one or more telemetry data bits. The operations performed are indicated in Fig. 11. In the case of phase-shift-keyed (PSK~ telemetry data modulation, the punctual I data components of the Z_data vector 286 (Fig. 7) are tested for sign. A plus sign is mapped into a binary one, a minus sign is mapped into a binary zero, and the raw Z_data vector 286 is then stripped of the data telemetry bits by inverting the sign of all component~ of a given correlator channel for which a binary zero ~as declared, prior to further processing of Z_data. !he recovered data telemetry bits, which for GPS are received at 50 baud, are collected in buffers and parity is checked. The data bits are tested for word and frame synchronization in a W O 93/22691 2 1 1 8 0 1 3 ~ PC~r/US93/04150 con~entional fashion as described for GPS in Interface Control Document ICD-200 tavailable from the NAVSTAR GPS Joint Program Office, United States Air Force Systems Command). .~lso, as described in this same document, the telemetry data is then used to compute satellite positions and drift correction parameters associated with satellite clock drift.
Periodically, the start and stop intervals t4 and tl of each correlator channel will need adjustment so as to remain in approximate alignment with the edges of the PS~ telemetry data. A standard bit synchronizer based upon the concept of the wide-sense cyclostationary process is incorporated for this purpose. (See Moghazy, Maral, and Blanchard, "Digilal PCM Bit Synchronizer and Detector", IEEE Transactions on Communication Systems, Vol. Com-28, August 1980, pp. 119l-1200.) In the present embodiment this synchronizationadjustment function is implemented by squaring values of the "on-time" I data stored in 288, and averaging these squared values over typically 100 correlation intervals to obtain a control signal which is responsive to the error in timing of correlation interval start relative to the time of data bit arrival at the receiver. This control signal is then used to periodically adjust the correlation interval start.
The Z_data vector produced by the correlator systenl 208 (Fig. 2) is used in the IQ combiner routine 250 (Fig. 3) to update the estimate sf.x of the true incremental LOS state vector _X 278 (Fig. 6). This key step can be performed using any one of a number of standard techniques, among them being the extended l~alman filter. A central part of any such method is the computation of the filter gain matrix which carries the measurement residuals into updated state estimates.
The signal processor accomplishes this function by means of the filter gain adjust routine 248 (Fig. 3~. It takes as input the signal amplitude estimates iq.am_hat 310 (Fig. 10), the nominal correlation measurement noise covariance matri~
3~ sf.r 296 (Fig. 7) that is user-supplied 238 (Fig. 3), and the path state noise covariance matrix se.qe 280 (Fig. 6) and produces the filter gain matrix sf.k. The detailed steps are shown in Figs. 12 through 16.
WO93/22691 21 18 0 i3 PCT/US93/04150 As discussed earlier, the signal processor services correlator channels on a "first come, first served" basis.
Before computing the residual errors in the LOS state estimates for the channel currently being served, the "best"
estimate of the LOS state based on the previously-served channel at its t4 time must be extrapolated to the t4 time of the current channel. It is important to note that the present invention, unlike conventional navigation receivers, performs an update of the LOS states of all correlator channels whenever fresh data becomes available from anv of its correlator channels.
The extrapolation of the LOS state from the most recently updated correlator channel (step 4) is described in Fig. 12 by equations 312. The vector sf.x_previous_last is the estimate of the LOS state following the most recently processed channel referenced to time t4 for that channel. The vector sf.x_current_last is the estimate of the LOS state for the channel currently being processed referenced to time t4 for that channel. The state transition matrix ~los was discussed earlier in connection with Fig. 6.
The filter state error covariance matrix sf.p is a product of the standard Kalman filter methodology. The error covariance matrix based on the previously-serviced channel at its t4 time is also extrapolated in step 4 to the t4 time of the current channel. The extrapolation process is specified in Fig. 12 by equation 314. The quantity sf.p_previous is the error covariance matrix at time t4[k-l,i-1] following the update of the LOS state at the last t4 time for the previous channel. The quantity sf.p_current is the error covariance matrix at time t4[k-l,i]. The error covariance matrix is extrapolated with the same time step as used in extrapolatin~
the LOS state 312. The quantity se.qe 280 (Fig. 6) was defined earlier.
Step 5 is detailed in Figs. 13, 14, and 15. The predicted data sf.z_hat 316 from the correlator system 208 (Fig. 2) and the velocity aiding system 23~ (Fig. 2) are subtracted from the current raw data Z_data to form the measurement residual vector sf.res 318 as indicated in Fig.
WO93/22691 2 1 18 ~ 13 PCT/US93/04150 13. At this time the seven residual components associated with the current channel are examined for reasonableness. A
channel with residuals above a rejection threshold is noted and its "bad residual count" is incremented.
There are two responses to a residual above threshold corresponding to whether the bad residual count has exceeded a limit or not:
(a) If the bad residual count has exceeded the limit, the current measurement is rejected which is most easily accomplished by increasing sf.r to sf.r_high, a quantity that is set to a large value, say 100,000. Then, this channel is removed from the list of active channels and is processed separately in a decoupled mode until residuals on this channel improve;
(b) If the bad residual count has not exceeded the limit, the current measurement is rejected as in the above case, but the channel is kept in the coupled mode.
The seven predicted data values 316 for the current LOS
path are defined by the equations 320 in Fig. 13 and all of the equations in Fig. 14. The quantities appearing in the equations are obtained from the estimate of the LOS state for the current channel sf.x_current_last 312 (Fig. 12) as shown in Fig. 15. The definition of the LOS state sf.x shown in Fig. 15 follows that given previously for the true LOS state 276, 278 (Fig.6).
The correlator system 208 (Fig. 2) is highly non-linear in the operating region of interest and is linearized about the current estimate of the LOS state vector in step 6.
kinearization is accomplished by ta~ing the gradient of the analytical expressions for the I and Q correlator outputs ~.-ith respect to the full set of LOS state variables. There are ~Ip to five LOS state variables defined for each line-of-sight path: phase and frequency error, modulating code offset and offset rate, and LOS velocity aiding input bias. For the i'th channel, the array of derivatives in the gradient are evaluated at sf.x_current_last and stored in array sf.h_i.
The gradient terms for the four-channel correlator system WO93/22691 2 1 1 8 0 1 3 ~ ` PCT/US93/~150 208 (Fig. 2) are defined by expressions 324, 326, and 328 in Fig. 16. The array 324 has dimensions 7x20 with three null partitions of dimensions 7x5.
The equation 330 for the Kalman filter gain sf.k_[i] for the i'th channel utilizes the current value for the correlator measurement noise covariance matrix, either sf.r_[i] if the residual is normal or sf.r_high_[i] if the residual is above the reasonableness threshold. The values in sf.r_[i] and sf.r_high_[i] are supplied as part of the initialization data 238 (Fig. 3). The overall LOS state covariance matri~
extrapolated to time t4[k-l,i], sf.p_current, in accordance with the standard Kalman filter design method.
The quantity sf.p_current is the filter state error covariance matrix extrapolated from the t4 time of the previously-updated correlator channel. The quantity sf.r_[i]
is the correlation measurement noise covariance matrix for the batch of seven measurement components for the current channel.
The dimensions of the arrays are 20x7 for sf.~_[i], 20x20 for sf.p_current, 7x20 with three zero partitions for sf.h_[i], and 7x7 for sf.r_[i].
Following the filter gain computation, the filter state error covariance matrix sf.p is updated 332 using the standard Kalman filter methodology. The updating reflects the improved estimation error for all correlator states following the processing of the current batch of correlator I and Q data together with the velocity aiding data. The quantity I in equation 332 is the identity matrix having dimensions 20x20.
It should be noted that the previously-calculated value of sf.p_previous is overwritten with the newly-calculated value.
It is also important to note that the new filter state error covariance matrix reflects improvement not onlv in the estimate for the LOS state for the current correlator channel but also in the estimates for the LOS states of all channels.
This result is in marked contrast to the traditional prior-art decoupled signal processing scheme.
The IQ combiner routine 250 (Fig. 3) shown in Figs. 1~
and 18 utilizes as inputs the residuals for the current channel sf.res_[i] and the filter gain matrix for the current WO93~22691 2 1 1 8 0 1 3 PCT/US93/041~0 channel sf.k_[i]. The output that results from the performance of the routine is the estimate of the LOS state sf.x for all LOS paths.
In step 7 the LOS state vector estimate sf.x (Fig. 15) referenced to the t4 time of the current channel is updated by multiplying the measurement residual vector sf.res by the filter gain matrix sf.k. in accordance with equation 334 shown in Fig. l7.
Following the updating of the LOS state vector estimate in step 7, the components of the LOS state base vector sf.base_last associated with the current channel are updated in step 8 by adding the updated LOS state values for the current channel to the LOS state base values for the current channel in accordance with equation 336. The term sf.x_base_last denotes the base values employed in each correlator channel for the k-l system cycle, the cycle that has just ended.
The LOS state values for the current channel, except for the velocity-aiding bias state, are then reset to zero in accordance with equation 338.
In step 9 the LOS state vector is next extrapolated to the start time t4[k,i] of the next correlation interval for the current channel, and the current channel is then defined to be the previous channel in accordance with equations 340 shown in Fig. 18. The vector sf.x_previous_next is the LOS
state estimate following the most recently completed channel referenced to the next t4 time for that channel. The vector sf.x_current_next is the LOS state estimate referenced to the next t4 time for the current channel.
In step l0 the base LOS values for the correlator system are then extrapolated to the current channel's ne~t correlation start t4[k,i] using the appropriate time step for each correlator channel in accordance with equations 3~2. The term sf.x_base_temp denotes extrapolations of sf.x_base to arbitrary time tags in the ~'th cycle.
The predicted whole values sf.x_whole_current for LOS
state estimates referenced to time t4[i,k] are formed b~-adding the updated and extrapolated LOS state vector to the W O 93/22691 PC~r/US93/04150 _base_temp values for each channel as indicated by equations 344.
The vehicle position and attitude evaluator routine 252 (Fig. 3) takes as initialization input: (1) the user-supplied initial value for platform position, velocity, attitude, and attitude rate sf.part; (2) the user-supplied initial values for the correlation measurement noise covariance matrix sf.r;
(3) the user-supplied values for platform velocity and platform rotation rate random process correlation times alpha_filt; and (4) the user-supplied values for platform velocity and platform rotation rate random process state noise disturbance power q_filt.
Quantities generated by the platform position and attitude evaluator routine 252 (Fig. 3) and utilized by other routines are: (1) the computed path state noise covariance matrix se.qe for the currently-defined satellite-to-platform antenna paths; and (2) the estimates of LOS whole phase and frequency values sf.x_base.
The details of the routine are shown in Figs. 19 through 24.
In step 11 the linearization of the LOS path oeometry is accomplished. The linearization relates small changes in platform navigation states (i.e. translation and rotation of the body-fixed coordinate system) to changes in the Euclidean distances and rates of change of distances measured between each satellite and the platform antenna(s). The linearization is found by taking the gradient of the set of distance vectors with respect to the navigation state vector and evaluating the resulting gradient matrix functions for the most recent estimate of the platform state vector se.state.
The line-of-sight linearization matrix 346 and the derivatives to be evaluated 348 are shown in Fig. 19.
As a prelude to the computation of the latest estimate of the platform na~igation vector, the offset vector iq.delt_x_los is obtained by subtracting the previous values of whole carrier phase and frequency from the current values for all channels, as indicated by equations 350 and 352.
Next, the pseudo-inverse matrix se.temp3 that relates W093/2269l 2118 013 - ~ ` PCT/US93/04150 _ changes in the line-of-sight quantities to the changes in navigation quantities is computed in accordance with equation 354. As is well known in the literature (see Gelb, APPlied OPtimal Estimation, MIT Press, 1974), the pseudo-inverse provides an optimal estimate in the least-squared-error sense for an overdetermined linear system of equations. This principle provides the basis for equation 356 which relates the current value sf.part_current of the platform navigation state to the previous value sf.part_previous. In accordance with equation 358, sf.part_current becomes sf.part_previous in anticipation of the next updating of sf.part.
As mentioned earlier, the central theme of the invention is the proper coupling of the signal tracking functions between all satellites being received. The quantity which permits this statistically optimum approach is the path state noise covariance matrix se.qe. The calculation of se.qe is accomplished by matrix operations involving the line-of-sight linearization matrix 346 and the platform motion state model noise covariance matrix se.qs. As wil be shown, the quantity se.qs is computed from the quantities alpha_filt and q_filt which are part of the user-supplied system initialization data 238 (Fig. 3).
The overall structure of se.qe 360 for n LOS paths is shown in Fig. 20. The QL terms contained in se.qe 360 describe the correlation in LOS paths due to user translations, rotations, and receiver clock bias and are derived from se.qe_los 362 which in turn is computed by means of equation 364 using the LOS linearization matrix previously discussed. The platform motion state model noise covariance matrix se.qs 366 is generated from input parameters as shown in Fig. 21. As indicated in the figure, the platform motion state noise covariance matrix 366 is partitioned into smaller arrays corresponding to the translation, cloc~ drift, and rotational motions of the platform. By way of an example, the 3~ partition 368 corresponding to the x-axis translation is shown. The elements of this 2-by-2 partition are computed from equations 370 where the various constants are obtained externally as components of the vectors alpha_filt 378 and WO93/22691 2 1 18 0 13 PCT/US93/~150 q_filt 380, discussed below in connection with Fig. 22. This computation represents the transformation between the continuous-time navigation model and the discrete-time navigation model.
The identical partitions QB of se.qe 360 (Fig. 20) are next computed. The QB partitions need only be computed once -at receiver initialization time. The 2-by-2 QB array 372 is shown in more detail in Fig. 22. The continuous-time state noise model describing the-carrier-code bias process is read in as components of the vectors alph_filt 378 and q_filt 380.
Then, the operations 374 needed to convert from the continuous-time model to the discrete-time model are performed resulting in the elements of QB.
Finally, the identical partitions QR of se.qe 360 (Fig.
l~ 20) are computed. The QR partitions also need onl~- be computed once - at receiver initialization time. The l-by-l QR partition 376 is defined in Fig. 22.
The computation of the new estimates of LOS whole phase and frequency values sf.x_base for the current channel requires the satellite positions and velocities to be updated.
Satellite position and velocity updating is accomplished conventionally in step 12 as described for GPS in Interface Control Document ICD-200 (available from the NAVSTAR GPS Joint Program Office, United States Air Force Systems Command).
First, the satellite location at its time of transmission corresponding to the beginning time t4[k,i] of the next correlation for the current channel is established. .~ signal observed at the platform at t4[k,i] will have been emitted at the satellite at an earlier time. The current estimate iq.phase_hat[i] of the line-of-sight distance from the satellite to the platform, converted from radian units to time units, is the appropriate time advance. The estimated time of transmission is computed using equation 382 shown in Fig. 23.
The quantity iq.phase_hat is obtained from sf.x_whole_current 3~ 344 (Fig. 18).
Next, the satellite position at the estimated time of transmission must be computed. For GPS the standard ICD-200 computation is employed to find the (x,y,z) earth-centered wo 93/226gl 2 1 18 0 13 PCT/US93/~150 . _ earth-fixed (ECEF) Cartesian coordinates corresponding to position ephem_position[i] 384 and to velocity - ephem_velocity[i] 386, the value of i corresponding to the satellite associated with the current channel.
The ECEF values obtained above are used in step 13 to compute the "best" ~stimates of LOS phase, ephem_LOS[i] 388, and frequency, ephem_LOS_dot[i] 390, for the current channel.
These values will be used for the next correlation interval.
If the externally-supplied differential corrections from a DGPS reference station are available, they are applied in step 14 to the above LOS values. The LOS and LOS rate values are also compensated for the user clock offset and offset rate, as previously estimated and stored in elements T and T_dot of the current navigation state vector sf.part_current.
1~
There are two possible system configurations for utilizing externally-generated differential corrections in obtaining solutions to the platform navigation problem. In one configuration the differential corrections for pseudorange and pseudorange-rate are generated in a conventional fashion at a reference station geographically separated from the platform using techniques and standards given in the ~GPS
RTCM-104 document and sent by radio link or other real-time communication medium to the platform where they are passed on to the coupled-filter signal processor 230 (Fig. 2) and utilized in the vehicle position and attitude evaluator routine 252 (Fig. 3).
For n LOS paths the psuedorange and pseudorange-rate differential corrections are defined by the expressions 392 and 394 respectively in Fig. 24. Although these differential corrections are generated by well-known techniques, their application directly in the signal tracking equations is unique to the present invention. By applying corrections in this fashion, the effects of satellite biases and propagation anomalies are removed from the phase and frequency estimates that are the products of the signal tracking equations. If this correction process were not performed at this stage, interpath line-of-sight correlations would be reduced thereby WO93/22691 2 1 18 0 1 3 ` PCT/US93/04150 causing the coupled-filter signal processor to be less effective.
There are certain situations where a real-time radio link between reference station and platform cannot be maintained, e.g. land surveying.
The technique employed for non-real-time operation is quite different from that employed for conventional DGPS.
Whereas in the conventional scheme pseudorange and carrier phase outputs from the reference station and platform are recorded for later merging in a signal processing scheme, the present invention requires that the sampled in-phase and quadrature outputs from the conventional receiver front-end (that form the inputs 206 (Fig. 2) to the correlation system 208 (Fig. 2)) be recorded by a high-speed digital recorder.
These recorded digital data streams together with recorded psueudorange and carrier phase outputs obtained from the reference station are then played at some convenient later time into the correlator system 208 (Fig. 2) and the coupled-filter signal processor 230 (Fig. 2). The computational process performed by the coupled-filter signal processor is the same as it was for the real-time situation. In this way the navigational history of a platform can be determined ir~ a manner fully-equivalent to the real-time process.
For either real-time or non-real-time operation the corrections are introduced as specified by equations 396 of Fig. 24. It should be noted that ephem_LOS_hat is expressed in radian units.
Finally, in step l5,the above corrected and compensated estimates are copied 398 into the first two components of sf.base. These base values, as well as the carrier bias values, replace the corresponding existing components in sf.x_base_next. The term sf.x_base_next denotes the base values to be employed in the h'th system cycle, the (k-l)'th cycle having just ended.
In step 16 satellite almanac data received over satellite telemetry links are used to evaluate on a regular basis the elevation angles of satellites currently being processed, as in conventional GPS receivers. The potential visibility of other satellites in the constellation not being tracked are also regularly checked using the same constellation almanac data. Satellites that are presently being tracked are taken out of the tracking mode .hen they fall below an elevation masking angle of typically five degrees. Fresh satellites are brought into the signal acquisition process if their predicted elevation angles are above the elevation masking angle. Once it is determined that a new satellite is to be acquired and tracked, an unused correlator channel is assigned to that satellite, and the signal acquisition process commences. Signal acquisition is essentially the same as for a conventional signal processor, as described in many references: e.g. M.~. Simon et al., Spread SPectrum Communications, Vol. III, Computer Science Press, Rockville, MD, 1985; C.R. Cahn and E.H. Martin, "Design Considerations for a Spread Spectrum Navigation Receiver,"
Proceedings of International Conference on Electronic Systems and Navigation Aids, Paris, France, November 1977. Note that during the initial acquisition process the coupled-filter signal processor 230 (Fig. 2) and the assigned correlator channel operate in a decoupled mode. Once code and frequency acquisition have been achieved, the correlator channel is controlled in the coupled mode.
In performing step 17 in the reference generator command setup routine 254 (Fig. 3), the signal processor 230 (Fig. 2) utilizes the data contained in sf.x_base_current for the current channel and generates commands to be sent to the reference generator 220 (Fig. 2) associated with the current channel which will cause the reference generator to generate a reference signal having specified values for carrier phase and frequency and modulating code phase during the next correlation interval. The reference generator command setup routine is detailed in Fig. 25.
Because of setup resolution limitations inherent in the 3~ design of the correlator system 208 (Fig. 2), the signal processor 230 (Fig.2) must provide for error setpoints which are non-zero. For example, in the case of GPS the code generator in the reference generator 220 (Fig. 2~ will only be WO93/22691 2 1 18 b i3 PCT/US93/04150 programmable in fractional code chip steps. Likewise, the phase and frequency of the signal produced by a numerically-controlled oscillator (NCO) in a typical GPS reference generator can only be adjusted in increments that are much coarser than the precision of the values maintained by the signal processor. In recognition of this situation, and to avoid any operating limitation postulated on zero LOS state values, provisions are made in the signal processor for non-zero set-points to be sent to the reference generators.
The commands sf.x_control to the reference generators are computed in accordance with equation 400. The base values come from the sf.x_base_next array (see equations 394 (Fig.
24) and discussion pertaining thereto).
The elements in the control vector correspond to NCO
phase and frequency setup and code generator setup for one correlator channel. The first component, which controls the NCO phase setup, is the LOS estimate iq.phase_base plus the carrier bias minus the desired NCO phase or setpoint. The second component, which controls the NCO frequency setup, is the estimated LOS rate plus the estimated carrier bias rate minus the desired NCO frequency offset or setpoint. The third component, which controls the code generator, is simply the LOS estimate iq.phase_base minus the desired code generator offset or setpoint.
The specified setpoints 402 are stored in read-only memory in the signal processor.
Claims (36)
1. A navigation receiver for use on platforms (or vehicles) of a variety of types, some types of said platforms being equipped with platform velocity sensors that provide velocity inputs to said receiver, said velocity sensors being characterized by a velocity bias parameter, some types of said platforms being equipped with attitude reference systems that provide inputs to said receiver of heading in some cases, heading and elevation in other cases, and heading, elevation, and yaw in still other cases, all of said attitude data being referenced to an earth-fixed coordinate system, said platforms being characterized by platform dynamics characterization data comprising platform velocity random process correlation time, platform rotation rate random process correlation time, platform velocity random process state noise disturbance power, and platform rotation rate random process state noise disturbance power, said receiver being capable of receiving a plurality of signals, said signals being transmitted by a plurality of spatially-distributed transmitters with known positions and velocities and being propagated to said receiver over a plurality of line-of-sight paths, said transmitted signals being modulated sinusoidal carriers, said transmitted signals being particular functions of time, said time being maintained in synchronized transmitter clocks, said transmitted signals being received through at least one of a plurality of spatially-distributed receiver ports, said received signals being characterized by signal parameters comprising carrier phase, carrier frequency, modulation phase, modulation phase rate-of-change, and carrier amplitude, said receiver having a clock, said receiver clock being maintained in synchronism with said transmitter clocks, said navigation receiver comprising:
a means for periodically measuring the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by a correlation measurement noise covariance matrix; and a means for periodically estimating the values of a plurality of said signal parameters for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values for each of said received signals than could be obtained by direct measurements of said signal parameters for each of said received signals or by estimates of said signal parameters based only on information derived from the received signal to which each pertains.
a means for periodically measuring the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by a correlation measurement noise covariance matrix; and a means for periodically estimating the values of a plurality of said signal parameters for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values for each of said received signals than could be obtained by direct measurements of said signal parameters for each of said received signals or by estimates of said signal parameters based only on information derived from the received signal to which each pertains.
2. The navigation receiver of claim 1 wherein said measuring means is comprised of:
a means for generating said reference signals, said reference signals being replicas of said transmitted signals, said reference signals being characterized by a carrier phase, a carrier frequency, and a modulation phase, said received signals and said reference signals being grouped according to the transmitter with which said signals are associated; and a means for correlating each of said reference signals with the received signal in its group.
a means for generating said reference signals, said reference signals being replicas of said transmitted signals, said reference signals being characterized by a carrier phase, a carrier frequency, and a modulation phase, said received signals and said reference signals being grouped according to the transmitter with which said signals are associated; and a means for correlating each of said reference signals with the received signal in its group.
3. The navigation receiver of claim 2 further comprising:
a means for adjusting said reference signal carrier phase, said reference signal carrier frequency, and said reference signal modulation phase in accordance with commands supplied at an input port of said adjusting means;
a means for generating said adjustment commands from said signal parameter estimates and for supplying said adjustment commands to said adjusting means, said adjustment commands to result in each of said reference signals having a carrier phase, carrier frequency, and modulation phase within predetermined increments of said estimated values.
a means for adjusting said reference signal carrier phase, said reference signal carrier frequency, and said reference signal modulation phase in accordance with commands supplied at an input port of said adjusting means;
a means for generating said adjustment commands from said signal parameter estimates and for supplying said adjustment commands to said adjusting means, said adjustment commands to result in each of said reference signals having a carrier phase, carrier frequency, and modulation phase within predetermined increments of said estimated values.
4. The navigation receiver of claim 1 wherein said estimating means comprises:
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
5. The navigation receiver of claim 2 wherein said estimating means comprises:
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
a means for adjusting said reference signal carrier phase, said reference signal carrier frequency, and said reference signal modulation phase in accordance with commands supplied at an input port of said adjusting means;
a means for generating said adjustment commands from said signal parameter estimates and for supplying said adjustment commands to said adjusting means, said adjustment commands to result in each of said reference signals having a carrier phase, carrier frequency, and modulation phase within predetermined increments of said estimated values.
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
a means for adjusting said reference signal carrier phase, said reference signal carrier frequency, and said reference signal modulation phase in accordance with commands supplied at an input port of said adjusting means;
a means for generating said adjustment commands from said signal parameter estimates and for supplying said adjustment commands to said adjusting means, said adjustment commands to result in each of said reference signals having a carrier phase, carrier frequency, and modulation phase within predetermined increments of said estimated values.
6. The navigation receiver of claim 1 wherein said estimating means comprises:
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically estimating carrier amplitude for each of said received signals utilizing said measured correlations;
a means for entering said correlation measurement noise covariance matrix into said navigation receiver thereby making said covariance matrix available for calculations performed within said navigation receiver;
a means for entering said platform dynamics characterization data into said navigation receiver thereby making said characterization data available for calculations performed within said navigation receiver;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters, said set of weighting coefficients constituting a row of the filter gain matrix, said set of weighting coefficients to be used in estimating the present value of the associated signal parameter, the generation of said filter gain matrix utilizing user-supplied, measured, and calculated data comprising said carrier amplitude estimates, said correlation measurement noise covariance matrix, and said platform dynamics characterization data;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically estimating carrier amplitude for each of said received signals utilizing said measured correlations;
a means for entering said correlation measurement noise covariance matrix into said navigation receiver thereby making said covariance matrix available for calculations performed within said navigation receiver;
a means for entering said platform dynamics characterization data into said navigation receiver thereby making said characterization data available for calculations performed within said navigation receiver;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters, said set of weighting coefficients constituting a row of the filter gain matrix, said set of weighting coefficients to be used in estimating the present value of the associated signal parameter, the generation of said filter gain matrix utilizing user-supplied, measured, and calculated data comprising said carrier amplitude estimates, said correlation measurement noise covariance matrix, and said platform dynamics characterization data;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
7. The navigation receiver of claim 1 further comprising a means for entering velocity data from said platform velocity sensor and platform attitude data from said attitude reference system into said navigation receiver thereby making said platform velocity and platform attitude data available for calculations performed within said navigation receiver, said signal parameter estimating means also providing estimates of the projections of said platform velocity bias on said line-of-sight paths of said received signals, said signal parameter estimating means also utilizing the projections of said platform velocity on said line-of-sight paths in obtaining said estimated signal parameter values for each of said received signals.
8. The navigation receiver of claim 1 further comprising a means for entering differential correction data into said navigation receiver thereby making said data available for calculations performed within said navigation receiver, said differential correction data compensating for errors in predicted positions and velocities of said reference transmitters, propagation anomalies, and clock drifts, said differential correction data comprising corrections to said estimates of carrier phase and carrier frequency for said plurality of received signals, said differential corrections being provided by an external source, said estimating means obtaining corrected estimates of carrier phase and carrier frequency by adding said differential corrections to said uncorrected estimates.
9. The navigation receiver of claim 1 further comprising a means for obtaining position and velocity data of said reference transmitters, said estimating means also obtaining present estimates of the position of said platform and the rate of change of said position utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals and said position and velocity data of said reference transmitters.
10. The navigation receiver of claim 9 wherein said estimating means also obtains present estimates of the attitude of said receiver platform and the rate of change of said attitude utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals.
11. A tracking receiver for use on fixed or mobile platforms in the cooperative tracking of other platforms, said other platforms being known as tracked platforms, said mobile tracking platforms being equipped with independent navigation means, said tracked platforms being of a variety of types, some types of said tracked platforms being equipped with platform velocity sensors that provide velocity inputs to said tracking receiver by means of a communication link, said velocity sensors being characterized by a velocity bias parameter, some types of said tracked platforms being equipped with attitude reference systems that provide inputs to said tracking receiver by means of said communication link of heading in some cases, heading and elevation in other cases, and heading, elevation, and yaw in still other cases, all of said attitude data being referenced to an earth-fixed coordinate system, said tracked platforms being characterized by platform dynamics characterization data comprising platform velocity random process correlation time, platform rotation rate random process correlation time, platform velocity random process state noise disturbance power, and platform rotation rate random process state noise disturbance power, said tracking receiver being capable of receiving signals through a plurality of spatially-distributed receiver ports, said signals being transmitted from said tracked platform to said tracking platform over a plurality of line-of-sight paths, said transmitted signals being modulated sinusoidal carriers transmitted through at least one transmit port, said received signals being characterized by signal parameters comprising carrier phase, carrier frequency, modulation phase, modulation phase rate-of-change, and carrier amplitude, said navigation receiver comprising:
a means for periodically measuring the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by a correlation measurement noise covariance matrix; and a means for periodically estimating the values of a plurality of said signal parameters for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values for each of said received signals than could be obtained by direct measurements of said signal parameters for each of said received signals or by estimates of said signal parameters based only on information derived from the received signal to which each pertains.
a means for periodically measuring the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by a correlation measurement noise covariance matrix; and a means for periodically estimating the values of a plurality of said signal parameters for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values for each of said received signals than could be obtained by direct measurements of said signal parameters for each of said received signals or by estimates of said signal parameters based only on information derived from the received signal to which each pertains.
12. The tracking receiver of claim 11 wherein said measuring means is comprised of:
a means for generating said reference signals, said reference signals being replicas of said transmitted signals, said reference signals being characterized by a carrier phase, a carrier frequency, and a modulation phase, said received signals and said reference signals being grouped according to the transmitter with which said signals are associated; and a means for correlating each of said reference signals with the received signal in its group.
a means for generating said reference signals, said reference signals being replicas of said transmitted signals, said reference signals being characterized by a carrier phase, a carrier frequency, and a modulation phase, said received signals and said reference signals being grouped according to the transmitter with which said signals are associated; and a means for correlating each of said reference signals with the received signal in its group.
13. The tracking receiver of claim 12 further comprising:
a means for adjusting said reference signal carrier phase, said reference signal carrier frequency, and said reference signal modulation phase in accordance with commands supplied at an input port of said adjusting means;
a means for generating said adjustment commands from said signal parameter estimates and for supplying said adjustment commands to said adjusting means, said adjustment commands to result in each of said reference signals having a carrier phase, carrier frequency, and modulation phase within predetermined increments of said estimated values.
a means for adjusting said reference signal carrier phase, said reference signal carrier frequency, and said reference signal modulation phase in accordance with commands supplied at an input port of said adjusting means;
a means for generating said adjustment commands from said signal parameter estimates and for supplying said adjustment commands to said adjusting means, said adjustment commands to result in each of said reference signals having a carrier phase, carrier frequency, and modulation phase within predetermined increments of said estimated values.
14. The tracking receiver of claim 11 wherein said estimating means comprises:
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
15. The tracking receiver of claim 12 wherein said estimating means comprises:
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
a means for adjusting said reference signal carrier phase, said reference signal carrier frequency, and said reference signal modulation phase in accordance with commands supplied at an input port of said adjusting means;
a means for generating said adjustment commands from said signal parameter estimates and for supplying said adjustment commands to said adjusting means, said adjustment commands to result in each of said reference signals having a carrier phase, carrier frequency, and modulation phase within predetermined increments of said estimated values.
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
a means for adjusting said reference signal carrier phase, said reference signal carrier frequency, and said reference signal modulation phase in accordance with commands supplied at an input port of said adjusting means;
a means for generating said adjustment commands from said signal parameter estimates and for supplying said adjustment commands to said adjusting means, said adjustment commands to result in each of said reference signals having a carrier phase, carrier frequency, and modulation phase within predetermined increments of said estimated values.
16. The tracking receiver of claim 11 wherein said estimating means comprises:
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically estimating carrier amplitude for each of said received signals utilizing said measured correlations;
a means for entering said correlation measurement noise covariance matrix into said navigation receiver thereby making said covariance matrix available for calculations performed within said navigation receiver;
a means for entering said platform dynamics characterization data into said navigation receiver thereby making said characterization data available for calculations performed within said navigation receiver;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters, said set of weighting coefficients constituting a row of the filter gain matrix, said set of weighting coefficients to be used in estimating the present value of the associated signal parameter, the generation of said filter gain matrix utilizing user-supplied, measured and calculated data comprising said carrier amplitude estimates, said correlation measurement noise covariance matrix, and said platform dynamics characterization data;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
a means for periodically estimating the present values of the correlations of said received signals with said reference signals utilizing the most recent estimates of signal parameter values;
a means for periodically estimating carrier amplitude for each of said received signals utilizing said measured correlations;
a means for entering said correlation measurement noise covariance matrix into said navigation receiver thereby making said covariance matrix available for calculations performed within said navigation receiver;
a means for entering said platform dynamics characterization data into said navigation receiver thereby making said characterization data available for calculations performed within said navigation receiver;
a means for periodically generating for each of said received signals a set of weighting coefficients for each of said signal parameters, said set of weighting coefficients constituting a row of the filter gain matrix, said set of weighting coefficients to be used in estimating the present value of the associated signal parameter, the generation of said filter gain matrix utilizing user-supplied, measured and calculated data comprising said carrier amplitude estimates, said correlation measurement noise covariance matrix, and said platform dynamics characterization data;
a means for periodically calculating for each of said received signals the estimated present value of each of said signal parameters, said calculation consisting of taking the differences between said present-value estimates of the correlations and the present measured values of the correlations for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
17. The tracking receiver of claim 11 further comprising a means for entering velocity data from said tracked platform velocity sensor and attitude data from said tracked platform attitude reference system into said tracking receiver thereby making said tracked platform velocity and said tracked platform attitude data available for calculations performed within said tracking receiver, said signal parameter estimating means also providing estimates of the projections of said platform velocity bias on said line-of-sight paths of said received signals, said signal parameter estimating means also utilizing the projections of said platform velocity on said line-of-sight paths in obtaining said estimated signal parameter values for each of said received signals.
18. The tracking receiver of claim 11 further comprising a means for entering differential correction data into said tracking receiver thereby making said data available for calculations performed within said tracking receiver, said differential correction data compensating for errors in predicted positions and velocities of said tracking platforms, propagation anomalies, and clock drifts, said differential correction data comprising corrections to said estimates of carrier phase and carrier frequency for said plurality of received signals, said differential corrections being provided by an external source, said estimating means obtaining corrected estimates of carrier phase and carrier frequency by adding said differential corrections to said uncorrected estimates.
19. The tracking receiver of claim 11 further comprising a means for obtaining position and velocity data of said tracking platform, said estimating means also obtaining present estimates of the position of said tracked platform and the rate of change of said position utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals and said position and velocity data of said tracking platform.
20. The tracking receiver of claim 19 wherein said estimating means also obtains present estimates of the attitude of said tracked platform and the rate of change of said attitude utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals.
21. A method of using the navigation receiver of claim 8 comprising the steps:
recording said received signals prior to measuring said correlations;
playing back said recorded signals into said correlation measuring means in place of said real-time received signals, said correlation measuring means thereby periodically measuring the correlations of each of said plurality of recorded received signals with a plurality of reference signals.
recording said received signals prior to measuring said correlations;
playing back said recorded signals into said correlation measuring means in place of said real-time received signals, said correlation measuring means thereby periodically measuring the correlations of each of said plurality of recorded received signals with a plurality of reference signals.
22. The method of claim 21 comprising the additional steps:
recording said differential correction data;
playing back said recorded differential correction data into said receiver in place of said real-time differential correction data in synchronism with playback of said recorded received signals thereby obtaining corrected estimates of carrier phase and carrier frequency by the adding of said differential corrections to said uncorrected estimates in said receiver.
recording said differential correction data;
playing back said recorded differential correction data into said receiver in place of said real-time differential correction data in synchronism with playback of said recorded received signals thereby obtaining corrected estimates of carrier phase and carrier frequency by the adding of said differential corrections to said uncorrected estimates in said receiver.
23. A method of using the navigation receiver of claim 18 comprising the steps:
recording said received signals prior to measuring said correlations;
playing back said recorded signals into said correlation measuring means in place of said real-time received signals, said correlation measuring means thereby periodically measuring the correlations of each of said plurality of recorded received signals with a plurality of reference signals.
recording said received signals prior to measuring said correlations;
playing back said recorded signals into said correlation measuring means in place of said real-time received signals, said correlation measuring means thereby periodically measuring the correlations of each of said plurality of recorded received signals with a plurality of reference signals.
24. The method of claim 23 comprising the additional steps:
recording said differential correction data;
playing back said recorded differential correction data into said receiver in place of said real-time differential correction data in synchronism with playback of said recorded received signals thereby obtaining corrected estimates of carrier phase and carrier frequency by the adding of said differential corrections to said uncorrected estimates in said receiver.
recording said differential correction data;
playing back said recorded differential correction data into said receiver in place of said real-time differential correction data in synchronism with playback of said recorded received signals thereby obtaining corrected estimates of carrier phase and carrier frequency by the adding of said differential corrections to said uncorrected estimates in said receiver.
25. A method of using a navigation receiver on platforms (or vehicles) of a variety of types, some types of said platforms being equipped with platform velocity sensors that provide velocity inputs to said receiver, said velocity sensors being characterized by a velocity bias parameter, some types of said platforms being equipped with attitude reference systems that provide inputs to said receiver of heading in some cases, heading and elevation in other cases, and heading, elevation, and yaw in still other cases, all of said attitude data being referenced to an earth-fixed coordinate system, said platforms being characterized by platform dynamics characterization data comprising platform velocity random process correlation time, platform rotation rate random process correlation time, platform velocity random process state noise disturbance power, and platform rotation rate random process state noise disturbance power, said receiver being capable of receiving a plurality of signals, said signals being transmitted by a plurality of spatially-distributed transmitters with known positions and velocities and being propagated to said receiver over a plurality of line-of-sight paths, said transmitted signals being modulated sinusoidal carriers, said transmitted signals being particular functions of time, said time being maintained in synchronized transmitter clocks, said transmitted signals being received through at least one of a plurality of spatially-distributed receiver ports, said received signals being characterized by signal parameters comprising carrier phase, carrier frequency, modulation phase, modulation phase rate-of-change, and carrier amplitude, said receiver having a clock, said receiver clock being maintained in synchronism with said transmitter clocks, said navigation receiver comprising:
a means for periodically measuring the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by the correlation measurement noise covariance matrix; and a means for periodically estimating the values of a plurality of said signal parameters for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values for each of said received signals than could be obtained by direct measurements of said signal parameters for each of said received signals or by estimates of said signal parameters based only on information derived from the received signal to which each pertains;
said method of using comprising the steps:
recording said received signals prior to measuring said correlations;
playing back said recorded signals into said correlation measuring means in place of said real-time received signals, said correlation measuring means thereby periodically measuring the correlations of each of said plurality of recorded received signals with a plurality of reference signals.
a means for periodically measuring the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by the correlation measurement noise covariance matrix; and a means for periodically estimating the values of a plurality of said signal parameters for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values for each of said received signals than could be obtained by direct measurements of said signal parameters for each of said received signals or by estimates of said signal parameters based only on information derived from the received signal to which each pertains;
said method of using comprising the steps:
recording said received signals prior to measuring said correlations;
playing back said recorded signals into said correlation measuring means in place of said real-time received signals, said correlation measuring means thereby periodically measuring the correlations of each of said plurality of recorded received signals with a plurality of reference signals.
26. The method of claim 25, said navigation receiver further comprising a means for entering differential correction data into said navigation receiver thereby making said data available for calculations performed within said navigation receiver, said differential correction data compensating for errors in predicted positions and velocities of said reference transmitters, propagation anomalies, and clock drifts, said differential correction data comprising corrections to said estimates of carrier phase and carrier frequency for said plurality of received signals, said differential corrections being provided by an external source, said estimating means obtaining corrected estimates of carrier phase and carrier frequency by adding said differential corrections to said uncorrected estimates, said method comprising the additional steps:
recording said differential correction data;
playing back said recorded differential correction data into said receiver in place of said real-time differential correction data in synchronism with playback of said recorded received signals thereby obtaining corrected estimates of carrier phase and carrier frequency by the adding of said differential corrections to said uncorrected estimates in said receiver.
recording said differential correction data;
playing back said recorded differential correction data into said receiver in place of said real-time differential correction data in synchronism with playback of said recorded received signals thereby obtaining corrected estimates of carrier phase and carrier frequency by the adding of said differential corrections to said uncorrected estimates in said receiver.
27. A method of using a tracking receiver on fixed or mobile platforms in the cooperative tracking of other platforms, said other platforms being known as tracked platforms, said mobile tracking platforms being equipped with independent navigation means, said tracked platforms being of a variety of types, some types of said tracked platforms being equipped with platform velocity sensors that provide velocity inputs to said tracking receiver by means of a communiction link, said velocity sensors being characterized by a velocity bias parameter, some types of said tracked platforms being equipped with attitude reference systems that provide inputs to said receiver by means of said communication link of heading in some cases, heading and elevation in other cases, and heading, elevation, and yaw in still other cases, all of said attitude data being referenced to an earth-fixed coordinate system, said tracked platforms being characterized by platform dynamics characterization data comprising platform velocity random process correlation time, platform rotation rate random process correlation time, platform velocity random process state noise disturbance power, and platform rotation rate random process state noise disturbance power, said receiver being capable of receiving signals through a plurality of spatially-distributed receiver ports, said signals being transmitted from said tracked platform to said tracking platform over a plurality of line-of-sight paths, said transmitted signals being modulated sinusoidal carriers transmitted through at least one transmit port, said received signals being characterized by signal parameters comprising carrier phase, carrier frequency, modulation phase, modulation phase rate-of-change, and carrier amplitude, said tracking receiver comprising:
a means for periodically measuring the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by the correlation measurement noise covariance matrix; and a means for periodically estimating the values of a plurality of said signal parameters for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values for each of said received signals than could be obtained by direct measurements of said signal parameters for each of said received signals or by estimates of said signal parameters based only on information derived from the received signal to which each pertains;
said method of using comprising the steps:
recording said received signals prior to measuring said correlations;
playing back said recorded signals into said correlation measuring means in place of said real-time received signals, said correlation measuring means thereby periodically measuring the correlations of each of said plurality of recorded received signals with a plurality of reference signals.
a means for periodically measuring the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by the correlation measurement noise covariance matrix; and a means for periodically estimating the values of a plurality of said signal parameters for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values for each of said received signals than could be obtained by direct measurements of said signal parameters for each of said received signals or by estimates of said signal parameters based only on information derived from the received signal to which each pertains;
said method of using comprising the steps:
recording said received signals prior to measuring said correlations;
playing back said recorded signals into said correlation measuring means in place of said real-time received signals, said correlation measuring means thereby periodically measuring the correlations of each of said plurality of recorded received signals with a plurality of reference signals.
28. The method of claim 27, said tracking receiver further comprising a means for entering differential correction data into said tracking receiver thereby making said data available for calculations performed within said tracking receiver, said differential correction data compensating for errors in predicted positions and velocities of said tracking platforms, propagation anomalies, and clock drifts, said differential correction data comprising corrections to said estimates of carrier phase and carrier frequency for said plurality of received signals, said differential corrections being provided by an external source, said estimating means obtaining corrected estimates of carrier phase and carrier frequency by adding said differential corrections to said uncorrected estimates, said method comprising the additional steps:
recording said differential correction data;
playing back said recorded differential correction data into said receiver in place of said real-time differential correction data in synchronism with playback of said recorded received signals thereby obtaining corrected estimates of carrier phase and carrier frequency by the adding of said differential corrections to said uncorrected estimates in said receiver.
recording said differential correction data;
playing back said recorded differential correction data into said receiver in place of said real-time differential correction data in synchronism with playback of said recorded received signals thereby obtaining corrected estimates of carrier phase and carrier frequency by the adding of said differential corrections to said uncorrected estimates in said receiver.
29. The method practiced by the navigation receiver of claim 10.
30. The method practiced by the tracking receiver of claim 20.
31. A method for obtaining navigation data from signals received from a plurality of spatially-distributed reference transmitters with known positions and velocities, said method being used on platforms (or vehicles) of a variety of types, some types of said platforms being equipped with platform velocity sensors, said velocity sensors being characterized by a velocity bias parameter, said velocity and velocity bias parameter being used in practicing said method, some types of said platforms being equipped with attitude reference systems that provide heading in some cases, heading and elevation in other cases, and heading, elevation, and yaw in still other cases, said attitude data being used in practicing said method, all of said attitude data being referenced to an earth-fixed coordinate system, said platforms being characterized by platform dynamics characterization data comprising platform velocity random process correlation time, platform rotation rate random process correlation time, platform velocity random process state noise disturbance power, and platform rotation rate random process state noise disturbance power, said platform dynamics data being used in practicing said method, said transmitted signals being modulated sinusoidal carriers, said transmitted signals being particular functions of time, said time being maintained in synchronized transmitter clocks, said transmitted signals being propagated to said platforms over a plurality of line-of-sight paths and being received through at least one of a plurality of spatially-distributed receiving ports, said received signals being characterized by signal parameters comprising carrier phase, carrier frequency, modulation phase, modulation phase rate-of-change, and carrier amplitude, said method comprising the steps:
measuring periodically the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by the correlation measurement noise covariance matrix; and estimating periodically the values of a plurality of said signal parameters and the projection of said velocity sensor bias on the line-of-sight path for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals, the projections of said platform velocity on said line-of-sight paths, said correlation measurement noise covariance matrix, and said platform dynamics data in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values and velocity bias projection for each of said received signals than could be obtained by direct measurements of said signal parameters and velocity bias projection for each of said received signals or by estimates of said signal parameters and velocity bias projection based only on information derived from the received signal to which each pertains.
measuring periodically the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by the correlation measurement noise covariance matrix; and estimating periodically the values of a plurality of said signal parameters and the projection of said velocity sensor bias on the line-of-sight path for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals, the projections of said platform velocity on said line-of-sight paths, said correlation measurement noise covariance matrix, and said platform dynamics data in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values and velocity bias projection for each of said received signals than could be obtained by direct measurements of said signal parameters and velocity bias projection for each of said received signals or by estimates of said signal parameters and velocity bias projection based only on information derived from the received signal to which each pertains.
32. The method of claim 31 wherein said estimating step comprises the steps:
estimating the present values of said correlations of said received signals with said reference signals and said velocity projection utilizing the most recent estimates of signal parameter values;
generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
calculating for each of said received signals the estimated present value of each of said signal parameters and said velocity bias projection, said calculation consisting of taking the differences between said present-value estimates of the correlations and velocity projection and the present measured values of the correlations and velocity projection for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter or velocity bias projection being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
estimating the present values of said correlations of said received signals with said reference signals and said velocity projection utilizing the most recent estimates of signal parameter values;
generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
calculating for each of said received signals the estimated present value of each of said signal parameters and said velocity bias projection, said calculation consisting of taking the differences between said present-value estimates of the correlations and velocity projection and the present measured values of the correlations and velocity projection for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter or velocity bias projection being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
33. The method of claim 31 comprising the additional steps:
obtaining position and velocity data of said reference transmitters;
obtaining present estimates of the position of said platform and the rate of change of said position utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals and said position and velocity data of said reference transmitters;
obtaining present estimates of the attitude of said platform and the rate of change of said attitude utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals.
obtaining position and velocity data of said reference transmitters;
obtaining present estimates of the position of said platform and the rate of change of said position utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals and said position and velocity data of said reference transmitters;
obtaining present estimates of the attitude of said platform and the rate of change of said attitude utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals.
34. A method of cooperatively tracking mobile platforms known as tracked platforms, said tracked platforms transmitting at least one signal, said method being practiced on fixed or mobile platforms known as tracking platforms, said mobile tracking platforms being equipped with independent navigation means, said tracked platforms being of a variety of types, some types of said tracked platforms being equipped with platform velocity sensors that provide velocity inputs to said tracking receiver by means of a communication link, said velocity sensors being characterized by a velocity bias parameter, said velocity and velocity bias parameter being used in practicing said method, some types of said platforms being equipped with attitude reference systems that provide inputs by means of said communication link of heading in some cases, heading and elevation in other cases, and heading, elevation, and yaw in still other cases, said attitude data being used in practicing said method, all of said attitude data being referenced to an earth-fixed coordinate system, said tracked platforms being characterized by platform dynamics characterization data comprising platform velocity random process correlation time, platform rotation rate random process correlation time, platform velocity random process state noise disturbance power, and platform rotation rate random process state noise disturbance power, said platform dynamics data being used in practicing said method, said signals being transmitted from said tracked platform to said tracking platform over a plurality of line-of-sight paths, said signals being received through a plurality of spatially-distributed receiving ports, said transmitted signals being modulated sinusoidal carriers transmitted through at least one transmit port, said received signals being characterized by signal parameters comprising carrier phase, carrier frequency, modulation phase, modulation phase rate-of-change, and carrier amplitude, said method comprising the steps:
measuring periodically the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by the correlation measurement noise covariance matrix; and estimating periodically the values of a plurality of said signal parameters and the projection of said velocity sensor bias on the line-of-sight path for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals, the projections of said platform velocity on said line-of-sight paths, said correlation measurement noise covariance matrix, and said platform dynamics data in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values and velocity bias projection for each of said received signals than could be obtained by direct measurements of said signal parameters and velocity bias projection for each of said received signals or by estimates of said signal parameters and velocity bias projection based only on information derived from the received signal to which each pertains.
measuring periodically the correlations of each of said plurality of received signals with a plurality of reference signals, said measured correlations having noise components characterized by the correlation measurement noise covariance matrix; and estimating periodically the values of a plurality of said signal parameters and the projection of said velocity sensor bias on the line-of-sight path for each of said plurality of received signals, said estimating means utilizing said correlation measurements for a plurality of said received signals, the projections of said platform velocity on said line-of-sight paths, said correlation measurement noise covariance matrix, and said platform dynamics data in obtaining said estimated signal parameter values for each of said received signals thereby obtaining more reliable estimates of said signal parameter values and velocity bias projection for each of said received signals than could be obtained by direct measurements of said signal parameters and velocity bias projection for each of said received signals or by estimates of said signal parameters and velocity bias projection based only on information derived from the received signal to which each pertains.
35. The method of claim 34 wherein said estimating step comprises the steps:
estimating the present values of said correlations of said received signals with said reference signals and said velocity projection utilizing the most recent estimates of signal parameter values;
generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
calculating for each of said received signals the estimated present value of each of said signal parameters and said velocity bias projection, said calculation consisting of taking the differences between said present-value estimates of the correlations and velocity projection and the present measured values of the correlations and velocity projection for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter or velocity bias projection being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
estimating the present values of said correlations of said received signals with said reference signals and said velocity projection utilizing the most recent estimates of signal parameter values;
generating for each of said received signals a set of weighting coefficients for each of said signal parameters;
calculating for each of said received signals the estimated present value of each of said signal parameters and said velocity bias projection, said calculation consisting of taking the differences between said present-value estimates of the correlations and velocity projection and the present measured values of the correlations and velocity projection for a plurality of received signals, multiplying each of said differences by the appropriate weighting coefficient from the set associated with the particular received signal and the particular signal parameter or velocity bias projection being calculated, summing the weighted differences over the plurality of correlations, and adding the weighted sum to the prior estimate of the value of the signal parameter extrapolated to the present.
36. The method of claim 34 comprising the additional steps:
obtaining position and velocity data of said tracking platform;
obtaining present estimates of the position of said tracked platform and the rate of change of said position utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals and said position and velocity data of said tracking platform;
obtaining present estimates of the attitude of said tracked platform and the rate of change of said attitude utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals.
obtaining position and velocity data of said tracking platform;
obtaining present estimates of the position of said tracked platform and the rate of change of said position utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals and said position and velocity data of said tracking platform;
obtaining present estimates of the attitude of said tracked platform and the rate of change of said attitude utilizing said present estimates of carrier phase, carrier frequency, and modulation phase for a plurality of said received signals.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US07/880,135 | 1992-05-07 | ||
US07/880,135 US5343209A (en) | 1992-05-07 | 1992-05-07 | Navigation receiver with coupled signal-tracking channels |
Publications (1)
Publication Number | Publication Date |
---|---|
CA2118013C true CA2118013C (en) | 1996-04-09 |
Family
ID=25375580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CA002118013A Expired - Lifetime CA2118013C (en) | 1992-05-07 | 1993-05-03 | Navigation receiver with coupled signal-tracking channels |
Country Status (6)
Country | Link |
---|---|
US (1) | US5343209A (en) |
EP (1) | EP0639279B1 (en) |
JP (1) | JP3271977B2 (en) |
CA (1) | CA2118013C (en) |
DE (1) | DE69309050T2 (en) |
WO (1) | WO1993022691A1 (en) |
Families Citing this family (144)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5422813A (en) * | 1992-12-17 | 1995-06-06 | Stanford Telecommunications, Inc. | No-outage GPS/commercial RF positioning system |
US5406291A (en) * | 1993-10-05 | 1995-04-11 | Grumman Aerospace Corporation | Passive emitter location system |
US5477458A (en) * | 1994-01-03 | 1995-12-19 | Trimble Navigation Limited | Network for carrier phase differential GPS corrections |
US5506781A (en) * | 1994-06-03 | 1996-04-09 | Itt Corporation | RF link control of satellite clocks |
US5583774A (en) * | 1994-06-16 | 1996-12-10 | Litton Systems, Inc. | Assured-integrity monitored-extrapolation navigation apparatus |
US6334219B1 (en) | 1994-09-26 | 2001-12-25 | Adc Telecommunications Inc. | Channel selection for a hybrid fiber coax network |
US5519605A (en) * | 1994-10-24 | 1996-05-21 | Olin Corporation | Model predictive control apparatus and method |
DE69420418T2 (en) * | 1994-11-29 | 2000-05-25 | Zanavy Informatics Kk | Navigation system with switchover if a radio signal cannot be received |
US7280564B1 (en) | 1995-02-06 | 2007-10-09 | Adc Telecommunications, Inc. | Synchronization techniques in multipoint-to-point communication using orthgonal frequency division multiplexing |
USRE42236E1 (en) | 1995-02-06 | 2011-03-22 | Adc Telecommunications, Inc. | Multiuse subcarriers in multipoint-to-point communication using orthogonal frequency division multiplexing |
US7092369B2 (en) | 1995-11-17 | 2006-08-15 | Symbol Technologies, Inc. | Communications network with wireless gateways for mobile terminal access |
US5787384A (en) * | 1995-11-22 | 1998-07-28 | E-Systems, Inc. | Apparatus and method for determining velocity of a platform |
US5832247A (en) * | 1995-12-28 | 1998-11-03 | Trimble Navigation Limited | PCI card for receiving a GPS signal |
US5828694A (en) * | 1996-07-01 | 1998-10-27 | Trimble Navigation Limited | Determination of multipath tracking error |
US9134398B2 (en) | 1996-09-09 | 2015-09-15 | Tracbeam Llc | Wireless location using network centric location estimators |
CA2265875C (en) * | 1996-09-09 | 2007-01-16 | Dennis Jay Dupray | Location of a mobile station |
US5768319A (en) * | 1997-02-18 | 1998-06-16 | Motorola, Inc. | GPS data collection in high noise-to-signal ratio environments |
US5910788A (en) * | 1997-04-17 | 1999-06-08 | Honeywell, Inc. | Predictive approach integrity |
EP1793240A1 (en) * | 1997-04-21 | 2007-06-06 | Raytheon Company | Increased jamming immunity by optimizing processing gain for GPS/INS systems |
US6018784A (en) * | 1997-07-31 | 2000-01-25 | Trimble Navigation Limited | PCI card for receiving a GPS signal |
US6531982B1 (en) | 1997-09-30 | 2003-03-11 | Sirf Technology, Inc. | Field unit for use in a GPS system |
US6327471B1 (en) | 1998-02-19 | 2001-12-04 | Conexant Systems, Inc. | Method and an apparatus for positioning system assisted cellular radiotelephone handoff and dropoff |
US6366599B1 (en) | 1998-03-16 | 2002-04-02 | Trimble Navigation Limited | Fast acquisition of spread-spectrum signals by dynamically varying spacing of search bins |
US6348744B1 (en) | 1998-04-14 | 2002-02-19 | Conexant Systems, Inc. | Integrated power management module |
US7545854B1 (en) | 1998-09-01 | 2009-06-09 | Sirf Technology, Inc. | Doppler corrected spread spectrum matched filter |
US7711038B1 (en) | 1998-09-01 | 2010-05-04 | Sirf Technology, Inc. | System and method for despreading in a spread spectrum matched filter |
US6693953B2 (en) | 1998-09-30 | 2004-02-17 | Skyworks Solutions, Inc. | Adaptive wireless communication receiver |
US8135413B2 (en) | 1998-11-24 | 2012-03-13 | Tracbeam Llc | Platform and applications for wireless location and other complex services |
US6331835B1 (en) * | 1999-02-02 | 2001-12-18 | The Charles Stark Draper Laboratory, Inc. | Deeply-integrated adaptive GPS-based navigator with extended-range code tracking |
US6448925B1 (en) | 1999-02-04 | 2002-09-10 | Conexant Systems, Inc. | Jamming detection and blanking for GPS receivers |
US6606349B1 (en) | 1999-02-04 | 2003-08-12 | Sirf Technology, Inc. | Spread spectrum receiver performance improvement |
US6577271B1 (en) | 1999-03-30 | 2003-06-10 | Sirf Technology, Inc | Signal detector employing coherent integration |
US6304216B1 (en) * | 1999-03-30 | 2001-10-16 | Conexant Systems, Inc. | Signal detector employing correlation analysis of non-uniform and disjoint sample segments |
US6351486B1 (en) | 1999-05-25 | 2002-02-26 | Conexant Systems, Inc. | Accelerated selection of a base station in a wireless communication system |
GB9912329D0 (en) * | 1999-05-26 | 1999-07-28 | Symmetricon Inc | Positioning apparatus |
US6243649B1 (en) * | 1999-08-02 | 2001-06-05 | Trimble Navigation Limited | GPS true course guidance correction system for correcting antenna offset error |
US6278403B1 (en) | 1999-09-17 | 2001-08-21 | Sirf Technology, Inc. | Autonomous hardwired tracking loop coprocessor for GPS and WAAS receiver |
US6731237B2 (en) | 1999-11-09 | 2004-05-04 | The Charles Stark Draper Laboratory, Inc. | Deeply-integrated adaptive GPS-based navigator with extended-range code tracking |
US6526322B1 (en) | 1999-12-16 | 2003-02-25 | Sirf Technology, Inc. | Shared memory architecture in GPS signal processing |
US6952440B1 (en) | 2000-04-18 | 2005-10-04 | Sirf Technology, Inc. | Signal detector employing a Doppler phase correction system |
US6788655B1 (en) | 2000-04-18 | 2004-09-07 | Sirf Technology, Inc. | Personal communications device with ratio counter |
US6714158B1 (en) | 2000-04-18 | 2004-03-30 | Sirf Technology, Inc. | Method and system for data detection in a global positioning system satellite receiver |
US6931055B1 (en) | 2000-04-18 | 2005-08-16 | Sirf Technology, Inc. | Signal detector employing a doppler phase correction system |
US7885314B1 (en) | 2000-05-02 | 2011-02-08 | Kenneth Scott Walley | Cancellation system and method for a wireless positioning system |
US6462708B1 (en) | 2001-04-05 | 2002-10-08 | Sirf Technology, Inc. | GPS-based positioning system for mobile GPS terminals |
US6671620B1 (en) | 2000-05-18 | 2003-12-30 | Sirf Technology, Inc. | Method and apparatus for determining global position using almanac information |
US6389291B1 (en) | 2000-08-14 | 2002-05-14 | Sirf Technology | Multi-mode global positioning system for use with wireless networks |
US7970412B2 (en) | 2000-05-18 | 2011-06-28 | Sirf Technology, Inc. | Aided location communication system |
US6427120B1 (en) | 2000-08-14 | 2002-07-30 | Sirf Technology, Inc. | Information transfer in a multi-mode global positioning system used with wireless networks |
US7949362B2 (en) | 2000-05-18 | 2011-05-24 | Sirf Technology, Inc. | Satellite positioning aided communication system selection |
US8078189B2 (en) | 2000-08-14 | 2011-12-13 | Sirf Technology, Inc. | System and method for providing location based services over a network |
US7929928B2 (en) | 2000-05-18 | 2011-04-19 | Sirf Technology Inc. | Frequency phase correction system |
US8116976B2 (en) | 2000-05-18 | 2012-02-14 | Csr Technology Inc. | Satellite based positioning method and system for coarse location positioning |
US7970411B2 (en) | 2000-05-18 | 2011-06-28 | Sirf Technology, Inc. | Aided location communication system |
US6778136B2 (en) | 2001-12-13 | 2004-08-17 | Sirf Technology, Inc. | Fast acquisition of GPS signal |
US6738713B2 (en) | 2000-05-26 | 2004-05-18 | Parthus (Uk) Limited | Positioning apparatus and method |
US10684350B2 (en) | 2000-06-02 | 2020-06-16 | Tracbeam Llc | Services and applications for a communications network |
US10641861B2 (en) | 2000-06-02 | 2020-05-05 | Dennis J. Dupray | Services and applications for a communications network |
US6873288B2 (en) * | 2000-06-14 | 2005-03-29 | Stephen B. Heppe | Enhanced GNSS receiver |
US7616705B1 (en) | 2000-07-27 | 2009-11-10 | Sirf Technology Holdings, Inc. | Monolithic GPS RF front end integrated circuit |
US6856794B1 (en) * | 2000-07-27 | 2005-02-15 | Sirf Technology, Inc. | Monolithic GPS RF front end integrated circuit |
US7236883B2 (en) | 2000-08-14 | 2007-06-26 | Sirf Technology, Inc. | Aiding in a satellite positioning system |
US7680178B2 (en) | 2000-08-24 | 2010-03-16 | Sirf Technology, Inc. | Cross-correlation detection and elimination in a receiver |
US6529829B2 (en) * | 2000-08-24 | 2003-03-04 | Sirf Technology, Inc. | Dead reckoning system for reducing auto-correlation or cross-correlation in weak signals |
US6931233B1 (en) | 2000-08-31 | 2005-08-16 | Sirf Technology, Inc. | GPS RF front end IC with programmable frequency synthesizer for use in wireless phones |
US7047023B1 (en) | 2000-12-01 | 2006-05-16 | Sirf Technology, Inc. | GPS RF front end IC with frequency plan for improved integrability |
US7747236B1 (en) | 2000-12-11 | 2010-06-29 | Sirf Technology, Inc. | Method and apparatus for estimating local oscillator frequency for GPS receivers |
US7113552B1 (en) | 2000-12-21 | 2006-09-26 | Sirf Technology, Inc. | Phase sampling techniques using amplitude bits for digital receivers |
US7671489B1 (en) | 2001-01-26 | 2010-03-02 | Sirf Technology, Inc. | Method and apparatus for selectively maintaining circuit power when higher voltages are present |
US6680703B1 (en) | 2001-02-16 | 2004-01-20 | Sirf Technology, Inc. | Method and apparatus for optimally tuning a circularly polarized patch antenna after installation |
US6703971B2 (en) | 2001-02-21 | 2004-03-09 | Sirf Technologies, Inc. | Mode determination for mobile GPS terminals |
US7076256B1 (en) | 2001-04-16 | 2006-07-11 | Sirf Technology, Inc. | Method and apparatus for transmitting position data using control channels in wireless networks |
US8244271B2 (en) * | 2001-05-21 | 2012-08-14 | Csr Technology Inc. | Distributed data collection of satellite data |
US7877104B2 (en) * | 2001-05-21 | 2011-01-25 | Sirf Technology Inc. | Method for synchronizing a radio network using end user radio terminals |
US7668554B2 (en) | 2001-05-21 | 2010-02-23 | Sirf Technology, Inc. | Network system for aided GPS broadcast positioning |
US7925210B2 (en) | 2001-05-21 | 2011-04-12 | Sirf Technology, Inc. | Synchronizing a radio network with end user radio terminals |
US8082096B2 (en) | 2001-05-22 | 2011-12-20 | Tracbeam Llc | Wireless location routing applications and architecture therefor |
US6532251B1 (en) | 2001-08-16 | 2003-03-11 | Motorola, Inc. | Data message bit synchronization and local time correction methods and architectures |
US6775319B2 (en) | 2001-08-16 | 2004-08-10 | Motorola, Inc. | Spread spectrum receiver architectures and methods therefor |
US6965754B2 (en) * | 2001-10-09 | 2005-11-15 | Motorola, Inc. | Satellite positioning system receiver with reference oscillator circuit and methods therefor |
US7948769B2 (en) | 2007-09-27 | 2011-05-24 | Hemisphere Gps Llc | Tightly-coupled PCB GNSS circuit and manufacturing method |
US6944540B2 (en) | 2002-03-28 | 2005-09-13 | Motorola, Inc. | Time determination in satellite positioning system receivers and methods therefor |
DE10219216B4 (en) * | 2002-04-29 | 2005-04-07 | Eads Astrium Gmbh | Method and terminal for processing a navigation signal in a satellite navigation system |
US6661371B2 (en) | 2002-04-30 | 2003-12-09 | Motorola, Inc. | Oscillator frequency correction in GPS signal acquisition |
US7362799B1 (en) * | 2002-06-27 | 2008-04-22 | Arraycomm Llc | Method and apparatus for communication signal resolution |
US6889052B2 (en) * | 2002-08-30 | 2005-05-03 | Motorola, Inc. | Method and apparatus for generating time of arrival estimates for use in determining a location |
US7885745B2 (en) | 2002-12-11 | 2011-02-08 | Hemisphere Gps Llc | GNSS control system and method |
US7689354B2 (en) | 2003-03-20 | 2010-03-30 | Hemisphere Gps Llc | Adaptive guidance system and method |
US8140223B2 (en) | 2003-03-20 | 2012-03-20 | Hemisphere Gps Llc | Multiple-antenna GNSS control system and method |
US8265826B2 (en) * | 2003-03-20 | 2012-09-11 | Hemisphere GPS, LLC | Combined GNSS gyroscope control system and method |
US8594879B2 (en) | 2003-03-20 | 2013-11-26 | Agjunction Llc | GNSS guidance and machine control |
US8686900B2 (en) * | 2003-03-20 | 2014-04-01 | Hemisphere GNSS, Inc. | Multi-antenna GNSS positioning method and system |
US8634993B2 (en) | 2003-03-20 | 2014-01-21 | Agjunction Llc | GNSS based control for dispensing material from vehicle |
US8138970B2 (en) | 2003-03-20 | 2012-03-20 | Hemisphere Gps Llc | GNSS-based tracking of fixed or slow-moving structures |
US9002565B2 (en) | 2003-03-20 | 2015-04-07 | Agjunction Llc | GNSS and optical guidance and machine control |
US8214111B2 (en) | 2005-07-19 | 2012-07-03 | Hemisphere Gps Llc | Adaptive machine control system and method |
US8271194B2 (en) | 2004-03-19 | 2012-09-18 | Hemisphere Gps Llc | Method and system using GNSS phase measurements for relative positioning |
US8190337B2 (en) | 2003-03-20 | 2012-05-29 | Hemisphere GPS, LLC | Satellite based vehicle guidance control in straight and contour modes |
WO2008024123A2 (en) | 2005-10-28 | 2008-02-28 | Sirf Technology, Inc. | Global positioning system receiver timeline management |
KR20070012308A (en) | 2003-09-02 | 2007-01-25 | 서프 테크놀러지, 인코포레이티드 | Signal processing system for satellite positioning signals |
US8138972B2 (en) * | 2003-09-02 | 2012-03-20 | Csr Technology Inc. | Signal processing system for satellite positioning signals |
US20060074558A1 (en) * | 2003-11-26 | 2006-04-06 | Williamson Walton R | Fault-tolerant system, apparatus and method |
US20050114023A1 (en) * | 2003-11-26 | 2005-05-26 | Williamson Walton R. | Fault-tolerant system, apparatus and method |
US7248964B2 (en) * | 2003-12-05 | 2007-07-24 | Honeywell International Inc. | System and method for using multiple aiding sensors in a deeply integrated navigation system |
US7365680B2 (en) * | 2004-02-10 | 2008-04-29 | Sirf Technology, Inc. | Location services system that reduces auto-correlation or cross-correlation in weak signals |
US20050209762A1 (en) * | 2004-03-18 | 2005-09-22 | Ford Global Technologies, Llc | Method and apparatus for controlling a vehicle using an object detection system and brake-steer |
US8583315B2 (en) | 2004-03-19 | 2013-11-12 | Agjunction Llc | Multi-antenna GNSS control system and method |
US7061425B2 (en) * | 2004-11-12 | 2006-06-13 | Motorola Inc. | Frequency error tracking in satellite positioning system receivers |
US7986725B2 (en) * | 2005-03-04 | 2011-07-26 | Nokia Corporation | Spread spectrum transmission systems |
US8199049B2 (en) * | 2005-05-26 | 2012-06-12 | Trimble Navigation Limited | GNSS post positioning with selected precision |
US8358242B2 (en) * | 2005-05-26 | 2013-01-22 | Trimble Navigation Limited | GNSS post positioning with elongated dither sequence |
US7227496B2 (en) | 2005-05-26 | 2007-06-05 | Trimble Navigation Limited | GPS rover station having position dithering for controlling accuracy of secure positions |
US7202813B2 (en) * | 2005-05-27 | 2007-04-10 | Trimble Navigation Limited | GPS reference system providing a synthetic offset vector for dithering a rover station position |
KR100609905B1 (en) * | 2005-06-30 | 2006-08-08 | 삼성전자주식회사 | Image forming apparatus and light scanning apparatus |
US7388539B2 (en) | 2005-10-19 | 2008-06-17 | Hemisphere Gps Inc. | Carrier track loop for GNSS derived attitude |
US7328114B2 (en) * | 2005-12-09 | 2008-02-05 | General Electric Company | Methods and systems for measuring a rate of change of frequency |
US7835832B2 (en) | 2007-01-05 | 2010-11-16 | Hemisphere Gps Llc | Vehicle control system |
US8311696B2 (en) | 2009-07-17 | 2012-11-13 | Hemisphere Gps Llc | Optical tracking vehicle control system and method |
USRE48527E1 (en) | 2007-01-05 | 2021-04-20 | Agjunction Llc | Optical tracking vehicle control system and method |
US8000381B2 (en) | 2007-02-27 | 2011-08-16 | Hemisphere Gps Llc | Unbiased code phase discriminator |
US7808428B2 (en) | 2007-10-08 | 2010-10-05 | Hemisphere Gps Llc | GNSS receiver and external storage device system and GNSS data processing method |
US9002566B2 (en) | 2008-02-10 | 2015-04-07 | AgJunction, LLC | Visual, GNSS and gyro autosteering control |
WO2009126587A1 (en) | 2008-04-08 | 2009-10-15 | Hemisphere Gps Llc | Gnss-based mobile communication system and method |
US8217833B2 (en) | 2008-12-11 | 2012-07-10 | Hemisphere Gps Llc | GNSS superband ASIC with simultaneous multi-frequency down conversion |
US8386129B2 (en) | 2009-01-17 | 2013-02-26 | Hemipshere GPS, LLC | Raster-based contour swathing for guidance and variable-rate chemical application |
US8031115B2 (en) * | 2009-01-21 | 2011-10-04 | O2Micro Inc. | Navigation system with a signal processor |
US8085196B2 (en) | 2009-03-11 | 2011-12-27 | Hemisphere Gps Llc | Removing biases in dual frequency GNSS receivers using SBAS |
US8401704B2 (en) | 2009-07-22 | 2013-03-19 | Hemisphere GPS, LLC | GNSS control system and method for irrigation and related applications |
US8174437B2 (en) | 2009-07-29 | 2012-05-08 | Hemisphere Gps Llc | System and method for augmenting DGNSS with internally-generated differential correction |
US8334804B2 (en) | 2009-09-04 | 2012-12-18 | Hemisphere Gps Llc | Multi-frequency GNSS receiver baseband DSP |
US8649930B2 (en) | 2009-09-17 | 2014-02-11 | Agjunction Llc | GNSS integrated multi-sensor control system and method |
US8548649B2 (en) | 2009-10-19 | 2013-10-01 | Agjunction Llc | GNSS optimized aircraft control system and method |
US8583326B2 (en) | 2010-02-09 | 2013-11-12 | Agjunction Llc | GNSS contour guidance path selection |
US8248300B2 (en) | 2010-06-02 | 2012-08-21 | Honeywell International Inc. | System and method of compensating for micro-jump events |
EP2555017B1 (en) * | 2011-08-03 | 2017-10-04 | Harman Becker Automotive Systems GmbH | Vehicle navigation on the basis of satellite positioning data and vehicle sensor data |
US9482536B2 (en) * | 2012-05-31 | 2016-11-01 | International Electronic Machines Corp. | Pose estimation |
US9470798B2 (en) | 2012-11-13 | 2016-10-18 | Trimble Navigation Limited | GNSS rover having subscribed precision |
US9223026B2 (en) | 2012-11-13 | 2015-12-29 | Trimble Navigation Limited | GNSS reference for subscribed precisions |
US9223007B2 (en) * | 2012-11-21 | 2015-12-29 | Raytheon Company | Kalman filtering with indirect noise measurements |
CN104459743B (en) * | 2014-11-27 | 2017-01-04 | 西安空间无线电技术研究所 | Carrier phase offset determination methods between a kind of coherent multicarrier modulated signal component |
US9515697B2 (en) * | 2015-03-30 | 2016-12-06 | Honeywell International Inc. | Scanning correlator for global navigation satellite system signal tracking |
DE102015119308B4 (en) * | 2015-11-10 | 2019-02-21 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Method and apparatus for providing data for a satellite navigation-based automatic landing to an aircraft |
CN106353778B (en) * | 2016-08-23 | 2018-11-30 | 北京东方计量测试研究所 | A kind of precision test system and method for navigation terminal |
US11841442B2 (en) * | 2022-01-18 | 2023-12-12 | General Dynamics Mission Systems, Inc. | Doppler tracking for signals of opportunity using rate-line detection |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4485383A (en) * | 1980-12-01 | 1984-11-27 | Texas Instruments Incorporated | Global position system (GPS) multiplexed receiver |
US4445118A (en) * | 1981-05-22 | 1984-04-24 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Navigation system and method |
US4578678A (en) * | 1983-11-14 | 1986-03-25 | The United States Of America As Represented By The United States National Aeronautics And Space Administration | High dynamic global positioning system receiver |
US4754465A (en) * | 1984-05-07 | 1988-06-28 | Trimble Navigation, Inc. | Global positioning system course acquisition code receiver |
US4785463A (en) * | 1985-09-03 | 1988-11-15 | Motorola, Inc. | Digital global positioning system receiver |
US4807256A (en) * | 1985-12-23 | 1989-02-21 | Texas Instruments Incorporated | Global position system receiver |
US4821294A (en) * | 1987-07-08 | 1989-04-11 | California Institute Of Technology | Digital signal processor and processing method for GPS receivers |
US4928106A (en) * | 1988-07-14 | 1990-05-22 | Ashtech Telesis, Inc. | Global positioning system receiver with improved radio frequency and digital processing |
US5185610A (en) * | 1990-08-20 | 1993-02-09 | Texas Instruments Incorporated | GPS system and method for deriving pointing or attitude from a single GPS receiver |
-
1992
- 1992-05-07 US US07/880,135 patent/US5343209A/en not_active Expired - Lifetime
-
1993
- 1993-05-03 WO PCT/US1993/004150 patent/WO1993022691A1/en active IP Right Grant
- 1993-05-03 DE DE69309050T patent/DE69309050T2/en not_active Expired - Lifetime
- 1993-05-03 EP EP93911000A patent/EP0639279B1/en not_active Expired - Lifetime
- 1993-05-03 JP JP51957893A patent/JP3271977B2/en not_active Expired - Fee Related
- 1993-05-03 CA CA002118013A patent/CA2118013C/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
WO1993022691A1 (en) | 1993-11-11 |
JPH07506434A (en) | 1995-07-13 |
DE69309050T2 (en) | 1997-10-16 |
JP3271977B2 (en) | 2002-04-08 |
DE69309050D1 (en) | 1997-04-24 |
EP0639279B1 (en) | 1997-03-19 |
EP0639279A1 (en) | 1995-02-22 |
US5343209A (en) | 1994-08-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2118013C (en) | Navigation receiver with coupled signal-tracking channels | |
US7512505B2 (en) | Ranging system for determining ranging information of a spacecraft | |
US8364401B2 (en) | Highly integrated GPS, Galileo and inertial navigation system | |
JP3240393B2 (en) | Tracking system for world positioning system | |
CA2283904C (en) | A system using leo satellites for centimeter-level navigation | |
US6525688B2 (en) | Location-determination method and apparatus | |
US5543804A (en) | Navagation apparatus with improved attitude determination | |
US5185610A (en) | GPS system and method for deriving pointing or attitude from a single GPS receiver | |
US6683563B2 (en) | System for determining precise orbit of satellite and method thereof | |
US6630904B2 (en) | Deeply-integrated adaptive INS/GPS navigator with extended-range code tracking | |
US5969674A (en) | Method and system for determining a position of a target vehicle utilizing two-way ranging | |
USRE40642E1 (en) | Reduced-power GPS-based system for tracking multiple objects from a central location | |
JP3408600B2 (en) | Position calculation method in satellite navigation system | |
US5570096A (en) | Method and system for tracking satellites to locate unknown transmitting accurately | |
US20080309552A1 (en) | System and method for gps acquisition using advanced tight coupling | |
CA2172546A1 (en) | Navigation system using re-transmitted gps | |
US6417800B1 (en) | Method for determining reference time error and an electronic device | |
AU2005226022A1 (en) | Method for back-up dual-frequency navigation during brief periods when measurement data is unavailable on one of two frequencies | |
Weill | A high performance code and carrier tracking architecture for ground-based mobile GNSS receivers | |
Soloviev et al. | Deeply Integrated GPS/Low‐Cost IMU for Low CNR Signal Processing: Concept Description and In‐Flight Demonstration | |
US10534087B1 (en) | Differential vector phase locked loop GPS reception method | |
Betz | Fundamentals of Satellite‐Based Navigation and Timing | |
Sennott et al. | The use of satellite geometry for prevention of cycle slips in a GPS processor | |
US20220196851A1 (en) | Position, Navigation and Timing System Architecture Based on Signals of Opportunity | |
EP0664008A4 (en) | Method and apparatus for tracking the position and velocity of airborne instrumentation. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
EEER | Examination request | ||
MKLA | Lapsed |
Effective date: 20130503 |
|
MKEC | Expiry (correction) |
Effective date: 20131009 |