US5621646A - Wide area differential GPS reference system and method - Google Patents
Wide area differential GPS reference system and method Download PDFInfo
- Publication number
- US5621646A US5621646A US08/681,087 US68108796A US5621646A US 5621646 A US5621646 A US 5621646A US 68108796 A US68108796 A US 68108796A US 5621646 A US5621646 A US 5621646A
- Authority
- US
- United States
- Prior art keywords
- gps satellite
- measurement
- sup
- reference stations
- ephemeris
- 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 - Fee Related
Links
Images
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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/05—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing aiding data
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
- G01S19/071—DGPS corrections
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
- G01S19/072—Ionosphere corrections
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
- G01S19/073—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections involving a network of fixed stations
- G01S19/074—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections involving a network of fixed stations providing integrity data, e.g. WAAS
Definitions
- the present invention relates generally to differential GPS reference systems and methods. More specifically, the present invention pertains to a wide area differential GPS (WDGPS) system and method that employs common view time transfer and minimum norm solution techniques to provide GPS satellite ephemeris and clock corrections to GPS users.
- WDGPS wide area differential GPS
- GPS global positioning system
- the system includes GPS satellites that each transmit L1 and L2 GPS signals.
- the L1 and L2 GPS signals comprise modulated carrier signals at the GPS L1 and L2 frequencies (1,575 and 1,227 MHz) which are each modulated by a coarse-acquisition pseudo-random (C/A) code, an encrypted precision pseudo-random (P) code, and navigation data.
- the navigation data includes the transmitting GPS satellite's orbit location and clock offset from GPS time.
- a user can determine the pseudo-range to each GPS satellite in view based on the time of arrival of the C/A code contained by the GPS signals transmitted by the GPS satellite. The user can compute a position fix based on the determined pseudo-ranges to the GPS satellites in its view and the locations and clock offsets for these GPS satellites contained in the navigation data of the GPS signals they transmit.
- the U.S. Department of Defense engages in selective availability and deliberately misrepresents the true location and clock offsets of the GPS satellites in the navigation data of the GPS signals they transmit.
- the ionosphere and the troposphere delay the arrival of these GPS signals at the user's location. Thus, these factors affect the accuracy and integrity of GPS for the user.
- GPS may be augmented by a network of ground-based reference stations.
- the reference stations monitor the health of the GPS signals transmitted by the GPS satellites in their view and generate corrections to these signals which are then transmitted to the user.
- the user can accurately determine his position based on the GPS signals he receives from the GPS satellites and the corrections to these signals he receives from the reference stations.
- This form of GPS is known as differential GPS (DGPS).
- position domain WDGPS position domain WDGPS and combines the information from multiple wide area reference stations (WRSs) in the position domain.
- WLSs wide area reference stations
- Such an approach is described in "Multi-Site Real-Time DGPS System Using Staffix Link: Operational Results", by D. Lapucha and M. Huff, Proceedings of the Fifth International Technical Meeting of the Satellite Division of the Institute of Navigation, Albuquerque, September 1992.
- each WRS acts as a local area DGPS (LDGPS) reference station. Based on its known fixed position, it computes for each GPS satellite in view a corresponding pseudo-range correction to the GPS satellite. Each WRS then transmits a stream of LDGPS pseudo-range corrections to the user. The user computes pseudo-ranges to the GPS satellites in his view based on the GPS signals received from these GPS satellites. The user then applies the pseudo-range corrections to the computed pseudo-ranges. If there are M WRSs, the user computes M position fixes and the corresponding covariance matrices. The user then uses the covariance matrices to form a single weighted position fix.
- LDGPS local area DGPS
- this approach does not involve the estimation of any of the underlying states of the GPS satellites, namely GPS satellite ephemeris and clock errors. Nor does it involve the estimation of the clock differences between the different reference stations. These type of errors behave quite differently and could be mitigated if they were separated. Since this approach allows these errors to be combined, it unfortunately does not allow a-priori information about each of the underlying states to be introduced into the estimation process. As a result, the accuracy of each computed pseudo-range correction decreases rapidly with age and user to reference station separation.
- this approach requires M independent data streams to be sent to the user if there are M different reference stations. Thus, it requires a rather large data bandwidth to implement.
- each WRS computes a pseudo-range correction to each GPS satellite in view.
- These pseudo-range corrections are sent to a wide area master station (WMS) rather than the user.
- WMS uses well known common view time transfer techniques to estimate the clock differences between the WRSs.
- the estimated clock differences are then removed from the pseudo-range corrections.
- these synchronized pseudo-range corrections are all referenced to a common master clock.
- the WMS then averages for each GPS satellite the synchronized pseudo-range corrections for the GPS satellite. These averaged pseudo-range corrections are then transmitted to the user and used by the user in the same way as described earlier for the position domain WDGPS approach.
- the measurement domain WDGPS approach requires less data bandwidth than the position domain WDGPS approach since a single averaged data stream for each GPS satellite is sent to the user rather than M different data streams.
- the pseudo-range corrections degrade rapidly with age and user to reference station separation since the underlying states of the GPS satellites are not estimated and therefore the user cannot leverage the known behavior of the underlying states to make accurate position fixes. This frailty is especially pronounced if the ephemeris errors associated with one or more of the GPS satellites is large or the ionosphere is introducing errors which vary spatially.
- state space WDGPS which uses the information from multiple WRSs to estimate the underlying states of the GPS satellites, the underlying state of the ionosphere, and the underlying clock differences between the WRSs.
- each WRS computes a pseudo-range residual (i.e., correction) for each GPS satellite in view.
- the computed pseudo-range residuals are sent to a WMS which uses one large estimator to simultaneously compute the clock differences between the WRSs and the ephemeris and clock corrections for the ephemeris and clock errors of each GPS satellite in view of the network.
- the estimator in this approach must solve a system of equations which contains 4K+(M-1) unknowns and MK knowns, where K is the number of satellites in view and M is the number of reference stations. This stems from the fact that each satellite has 3 unknown components of ephemeris error and 1 unknown clock error. In addition, there are M-1 unknown offsets between the clocks of the WRSs clocks. At the same time, the system of equations contains MK knowns corresponding to the K pseudo-range residuals computed by each of the M WRSs.
- the state space WDGPS approach separates the underlying states of the system, and as such, it enables the use of a-priori models and side information for these states.
- its complexity impedes the straight forward integration of such models.
- a GPS reference system and corresponding method utilizing a network of reference stations and a master station.
- Each reference station receives GPS signals from GPS satellites it observes and computes in response a pseudo-range residual for each GPS satellite it observes.
- the master station is in communication with the reference stations to receive the computed pseudo-range residuals and comprises a pseudo-range residual synchronizer, an ephemeris and clock correction estimator for each GPS satellite observed by the reference stations, and a transmitter.
- the pseudo-range residual synchronizer is responsive to the received pseudo-range residuals in order to compute clock differences between the reference stations and remove the clock differences from the received pseudo-range residuals so as to synchronize them.
- Each ephemeris and clock correction estimator computes an ephemeris correction and a clock correction for the corresponding GPS satellite in response to each synchronized pseudo-range residual for that GPS satellite.
- the transmitter transmits to users the computed ephemeris and clock corrections for each GPS satellite observed by the reference stations.
- FIG. 1 shows a wide area differential GPS (WDGPS) system in accordance with the present invention.
- WDGPS wide area differential GPS
- FIG. 2 provides a block diagram of a wide area reference station (WRS) used in the WDGPS system of FIG. 1.
- WRS wide area reference station
- FIG. 3 shows the operation of the WRS of FIG. 2.
- FIG. 4 provides a block diagram of a wide area master station (WMS) used in the WDGPS system of FIG. 1.
- WMS wide area master station
- FIG. 5 shows the operation of the WMS of FIG. 4.
- FIG. 6 demonstrates the situation where a GPS satellite is observed by only three WRSs in the WDGPS system of FIG. 1.
- FIG. 7 demonstrates the situation where a GPS satellite is observed by only two WRSs in the WDGPS system of FIG. 1.
- FIG. 8 demonstrates the situation where a GPS satellite is observed by only one WRSs in the WDGPS system of FIG. 1.
- FIG. 9 provides a block diagram of a user's GPS positioning unit used in the WDGPS system of FIG. 1.
- FIG. 10 shows the operation of the user's GPS positioning unit of FIG. 9.
- a wide area differential GPS (WDGPS) network 100 that includes a WDGPS reference or augmentation system 104.
- the WDGPS reference system includes a network of wide area reference stations (WRSs) 102.
- Each WRS receives L1 and L2 GPS signals from each of the GPS satellites 106 that are in its view.
- each WRS observes each of the GPS satellites in its view and computes a pseudo-range residual and an ionospheric delay estimate for each observed GPS satellite.
- the pseudo-range residuals and ionospheric delay estimates computed by each WRS are provided to the wide area master station (WMS) 108 of the WDGPS reference system via a data communication link 110.
- WMS wide area master station
- the WMS in response computes ephemeris and clock corrections for each GPS satellite observed by the WRSs and computes a grid of ionospheric corrections.
- the WMS then transmits the ephemeris, clock, and ionospheric corrections to the geostationary satellite 112.
- the geostationary satellite 112 re-transmits these corrections so that a mobile user's GPS positioning system 114 receives them along with the L1 GPS signals from each GPS satellite it observes.
- the user's GPS positioning system then computes its position based on the received GPS signals and the received ephemeris, clock, and ionospheric corrections.
- each WRS 102 includes a conventional survey quality dual frequency GPS receiver 116, a meteorological sensor 118, a data processor 120, and a data link interface 121.
- the GPS receiver includes a conventional GPS antenna 122 for receiving L 1 and L2 GPS signals.
- the data processor includes a central processing unit (CPU) 124 that runs the pseudo-range smoother 126, the navigation data processor 128, the tropospheric correction estimator 130, the pseudo-range residual estimator 132, and the clock steerer 134 which are all stored in the memory 136.
- the memory also stores the precisely surveyed location 138 of the WRS.
- the GPS antenna 122 of each WRS 102 receives L1 and L2 GPS signals from each of the GPS satellites 106 in view of the WRS and the GPS receiver 116 in response observes each of these GPS satellites.
- the GPS receiver computes at each of its measurement epochs L1 and L2 pseudo-ranges ⁇ m ,L1 k .sbsp.m and ⁇ m ,L2 k .sbsp.m (from the WRS to the corresponding GPS satellite) based on the L1/CA and L2/CA codes received from the GPS satellite, where k m represents the kth GPS satellite observed by the mth WRS.
- the GPS receiver makes at each epoch L1 and L2 carrier phase measurements ⁇ m ,L1 k .sbsp.m and ⁇ m ,L2 k .sbsp.m of the L1 and L2 carrier signals transmitted by each GPS satellite it observes.
- the measurement epochs occur at one second intervals (i.e., 1 Hz).
- the pseudo-range smoother uses the corresponding L1 and L2 carrier phase measurements and the corresponding L1 and L2 pseudo-ranges to produce at each epoch a corresponding ionospheric delay estimate l m k .sbsp.m and a corresponding smoothed pseudo-range ⁇ m k .sbsp.m from which the ionospheric delay estimate has been removed. This is done using the method described in "GPS Survey", by C. Goad, a chapter in "The Global Positioning System: Theory and Applications", edited by B. Parkinson, J. Spilker, P. Axelrod, and P.
- the GPS receiver 116 of each WRS 102 also extracts the navigation data contained by the L1 and L2 GPS signals it receives from each GPS satellite it observes.
- the navigation data processor 128 is provided with the navigation data for each observed GPS satellite and unpacks from it the GPS satellite's location vector r k .sbsp.m.sup.,B and the difference (or offset) B k .sbsp.m.sup.,B from GPS time of the GPS satellite's clock.
- the navigation data processor computes at each epoch a line of sight vector 1 m k .sbsp.m and a nominal range r k .sbsp.m.sup.,B ⁇ 1 m k .sbsp.m from the WRS to the GPS satellite based on the unpacked GPS satellite's location and the known reference station location 138 stored in the memory 130.
- the conventional meteorological sensor 118 of each WRS 102 senses surface pressure, humidity, and/or temperature and provides this sensed meteorological data to the tropospheric correction estimator 130.
- the tropospheric delay estimator 130 periodically computes at each epoch a tropospheric delay estimate T m k .sbsp.m for each of the GPS satellites 106 observed by the WRS in response to the sensed meteorological data and the line of sight to each GPS satellite. This is done using a modified Hopfield method as described in "Analysis of Tropospheric Delay Prediction Models: Comparisons with Ray-Tracing and Implications for GPS Relative Positioning", by H. Janes, R Langley, and S.
- the set of tropospheric delay estimates ⁇ T m k .sbsp.m ⁇ k .sbsb.m 1 K .sbsp.m are then provided to the pseudo-range residual estimator 132.
- Each smoothed pseudo-range ⁇ m k .sbsp.m received at each epoch by the pseudo-range residual estimator 132 may be modeled as:
- r k .sbsp.m is a vector that represents the actual location of the kth GPS satellite 106 observed by the WRS
- b m represents the actual difference from GPS time of the clock of the GPS receiver 116 of the mth WRS 102
- B k .sbsp.m represents the actual difference from GPS time of the clock of the kth GPS satellite observed by the mth WRS
- ⁇ m k .sbsp.m represents the noise associated with the kth GPS satellite and the mth WRS that observes it.
- the pseudo-range residual estimator 132 subtracts the corresponding nominal range r k .sbsp.m.sup.,B ⁇ 1 m k .sbsp.m, unpacked clock offset B k .sbsp.m.sup.,B, and tropospheric correction T m k .sbsp.m to produce a pseudo-range residual (i.e. correction) ⁇ m k .sbsp.m at each epoch for the error in the smoothed pseudo-range.
- This pseudo-range residual may be represented as:
- ⁇ r k .sbsp.m r k .sbsp.m -r k .sbsp.m.sup.,B and is a vector representing the correction for the ephemeris error of the kth GPS satellite observed by the mth WRS
- ⁇ B k .sbsp.m B k .sbsp.m -B k .sbsp.m.sup.,B and represents the correction for the clock error of the kth GPS satellite observed by the mth WRS.
- the clock difference estimate b m is subtracted by the pseudo-range residual estimator 132 from each pseudo-range residual ⁇ m k .sbsp.m to produce an adjusted pseudo-range residual ⁇ m k .sbsp.m given by:
- Each data communication link 110 may comprise a conventional modem/telephone line system, an RF or other frequency transmitter/receiver system, a fiber optic cable network, or other suitable transmission system.
- WMS Wide Area Master Station
- the WMS 108 includes data link interfaces 140, a data processor 142, and a transmitter 143.
- the data processor includes a pseudo-range residual synchronizer 144, ephemeris and clock correction estimators 146, and an ionospheric grid generator 148 which are run by the CPU 147 and stored in the memory 150.
- the WMS 108 includes the same number of data link interfaces 140 as there are WRSs 102 in the WDGPS reference system 104.
- the WMS may include a GPS receiver that receives the GPS signals transmitted by the GPS satellites 104 and then removes the navigation data contained in the GPS signals.
- the navigation data in the GPS signals received by the WRSs could be simply provided to the WMS via the communication links 110. In either case, based on the GPS satellite locations contained in the navigation data of the GPS signals and the known locations of the WRSs which would be stored in the memory 150 of the WMS, the WMS would then compute the lines of sight.
- the grid of ionospheric corrections is provided to the user's positioning system 114 in a correction message transmitted by the transmitter 143.
- the pseudo-range residual synchronizer 144 For each GPS satellite 106 in common view of the primary WRS and at least one of the secondary WRSs, the pseudo-range residual synchronizer 144 first subtracts at each epoch the pseudo-range residual ⁇ M k .sbsp.m,M computed for it by the primary WRS from each pseudo-range residual ⁇ m k .sbsp.m,M computed for it by those of the secondary WRSs that observe it, where M represents the primary WRS (in addition to representing the number of WRSs) and, for m ⁇ M,k m ,M represents the kth of the GPS satellites that is observed by the mth secondary WRS and by the primary WRS. This results in the difference ⁇ m k .sbsp.m,M at each epoch which is given as follows:
- ⁇ m ,M k .sbsp.m,M ⁇ m k .sbsp.m,M - ⁇ M k .sbsp.m,M.
- an estimate of the WRS clock error can be made by averaging the difference ⁇ m k .sbsp.m,M as follows: ##EQU1## where ⁇ b m - ⁇ b M represents the clock difference at each epoch between the primary WRS and the mth secondary WRS and K m ,M represents the total number of GPS satellites in common view of the primary WRS and the mth WRS.
- Equation (5) the expression in Equation (5) can be re-written as: ##EQU2## where h(t j ) represents the well known impulse response for a low pass filter.
- the secondary WRSs in order to compute a clock difference ⁇ b m - ⁇ b M at each epoch for each WRS 102, the secondary WRSs must be positioned with respect to the primary WRS so that each secondary WRS will always have at least one GPS satellite in common view with the primary WRS.
- the WRS that is closest to the WMS and most centrally located with respect to the more peripheral WRSs is the primary WRS.
- each synchronized pseudo-range residual can now be represented as ⁇ m k .sbsp.m and given by:
- ⁇ M k .sbsp.m ⁇ M k .sbsp.m
- ⁇ B k .sbsp.m ⁇ b M - ⁇ B k .sbsp.m and represents the clock correction for the kth GPS satellite observed by the mth WRS synchronized to the clock of the primary WRS.
- Equation (7) can be re-written as follows:
- k represents the kth of all of the GPS satellites 106 observed by the entire network of WRSs 102 and m k represents the mth of those of the WRSs that observe the kth GPS satellite.
- the WMS 108 can employ an ephemeris and clock correction estimator 146 for each of the GPS satellites 106 observed by the WRSs.
- the computational complexity of a single large estimator is avoided.
- each ephemeris and clock correction estimator 146 solves for the four dimensional state vector [ ⁇ r k , ⁇ B k ] T based on the following matrix system: ##EQU3##
- Equation (9) The system of Equation (9) can be more compactly re-written as follows:
- ⁇ k is the vector containing each synchronized pseudo-range residual ⁇ m k computed by the WRSs 102 for the kth GPS satellite 106
- x k is the four dimensional state vector containing the ephemeris and clock corrections ⁇ r k and ⁇ B k for the kth GPS satellite
- G W is the design matrix that contains, for each WRS that observes the kth GPS satellite, a line of sight from it to the kth GPS satellite and a matrix element (i.e., -1) that indicates that the corresponding synchronized pseudo-range residual in the vector ⁇ k should be adjusted by the clock correction ⁇ B k for the kth GPS satellite.
- the WDGPS reference system 104 includes at least four WRSs 102 (i.e., the primary WRS and at least three secondary WRSs) and the kth GPS satellite 106 is observed by exactly four of these WRSs, then the system of Equation (10) is exactly specified. Furthermore, when the WDGPS reference system 104 includes more than four WRSs 102 and the kth GPS satellite 106 is observed by more than four of these WRSs, then the system of Equation (10) is over specified. In either case, the corresponding estimator 146 first computes a pseudo-inverse (G W T G W ) -1 G W T of the design matrix G w .
- the pseudo-inverse reduces to the inverse G W -1 of the design matrix G W .
- the estimator then computes a minimum mean square error solution for the state vector x k as follows:
- the system of Equation (10) is under specified when the kth GPS satellite 106 is observed by less than four WRSs 102, regardless of whether the WDGPS reference system 104 includes less than four WRSs (i.e., the primary WRS and no more than two secondary WRSs) or at least four WRSs.
- the corresponding estimator 146 computes a pseudo-inverse G W T (G W G W T ) -1 of the design matrix G W and then computes a minimum norm solution for the state vector x k as follows:
- the minimum norm solution has the minimum length or energy among all of the satisfactory solutions.
- the transmitter transmits them as corrections to the user's GPS positioning system 114 via the geostationary satellite 112 or other suitable broadcast facility so that the user's positioning system can compute an accurate position fix.
- the computed ephemeris and clock corrections may be transmitted directly to the user's positioning system by the transmitter 143.
- FIGS. 6 illustrates the situations when the kth GPS satellite 106 is observed by only three WRSs 102 and the system of Equation (10) is under specified. Since the kth GPS satellite 106 is observed by only three WRSs, this geometry does not provide leverage for precisely observing each of the radial (i.e., orbit radius), along track (along satellite velocity), and cross track (i.e., orthogonal to radial and along track) components of the ephemeris error of the GPS satellite and the clock error of the GPS satellite.
- the radial i.e., orbit radius
- along track along satellite velocity
- cross track i.e., orthogonal to radial and along track
- Equation (12) does provide accurate ephemeris and clock corrections ⁇ r k and ⁇ B k for the ephemeris and clock errors since the components of these errors that effect the user's positioning system effect the three WRSs substantially in the same manner.
- FIG. 7 demonstrates the case when the kth GPS satellite 106 is observed by only two WRSs 102 and the system of Equation (10) is under specified. Similar to the situation described for FIG. 6, although this geometry does not provide leverage for precisely observing each of the components of the ephemeris and clock errors, the minimum norm solution of Equation (12) does provide accurate ephemeris and clock corrections ⁇ r k and ⁇ B k for when the user's positioning system 114 is between the two WRSs or adjacent to one of the WRSs. As in the situation of FIG. 6, this is due to the fact that the components of the ephemeris and clock errors that effect the user's positioning system effect the two WRSs substantially in the same manner.
- FIG. 8 shows the situation when the kth GPS satellite 106 is observed by only one WRS 102 and the system of Equation (10) is under specified.
- This geometry also does not provide the leverage for accurately observing the components of the ephemeris and clock errors.
- the minimum norm solution of Equation (12) provides accurate ephemeris and clock corrections ⁇ r k and ⁇ B k because the components of the ephemeris and clock errors that effect the user's positioning system effect the WRS in substantially the same manner.
- the WDGPS reference system 104 is capable of estimating applicable ephemeris and clock errors for any GPS satellite in view of at least one of the WRSs 102, provided that the user's positioning system 114 is between WRSs or adjacent to a WRSs.
- the WDGPS reference system can provide to users accurate and valid ephemeris and clock corrections for each GPS satellite from the time that it is first observed by one of the WRSs (i.e., appears on the horizon) until the time that it is no longer observed by any of the WRSs (i.e., disappears below the horizon).
- the high frequency energy in the clock correction ⁇ B k (resulting from the high frequency energy in the clock error due to selective availability) splatters into the ephemeris corrections ⁇ r k .
- each ephemeris and clock correction estimator 146 first subtracts each synchronized pseudo-range residual ⁇ m .sbsb.k k for the corresponding GPS satellite computed by the secondary WRSs from the synchronized pseudo-range residual ⁇ M .sbsb.k k for the corresponding GPS satellite computed by the primary WRS, where M k also represents the primary WRS (in addition to representing the total number of WRSs that observe the kth GPS satellite).
- Each difference ⁇ m .sbsb.k M .sbsb.k k can be expressed as:
- ⁇ k is a vector containing each computed difference ⁇ m .sbsb.k.sub.,M.sbsb.k k
- G W is the matrix that contains, for each secondary WRS that observes the kth GPS satellite, the difference between the line of sight 1 m .sbsb.k k from the corresponding secondary reference station to the kth GPS satellite and the line of sight 1 M .sbsb.k k from the primary reference station to the kth GPS satellite.
- the system of Equation (15) is exactly specified when the WDGPS reference system 104 includes at least four WRSs 102 (i.e., the primary WRS and at least three secondary WRSs) and the kth GPS satellite 106 is observed by exactly four of these WRSs. And, the system of Equation (15) is over specified when the WDGPS reference system includes more than four WRSs 102 and when the kth GPS satellite 106 is observed by more than four of these WRSs.
- the corresponding estimator 146 in this case computes a three dimensional minimum mean square error solution for the state vector ⁇ r k as follows:
- the system of Equation (15) is under specified if the kth GPS satellite 106 is observed by less than four WRSs 102 in either situation where the WDGPS reference system 104 includes less than four WRS (i.e., the primary WRS and no more than two secondary WRSs) or at least four WRSs.
- the corresponding estimator 146 computes a three dimensional minimum norm solution for the state vector ⁇ r k as follows:
- each estimator 146 then computes the clock correction ⁇ B k according to: ##EQU5##
- the three dimensional estimators 146 prevent the high frequency energy in the clock correction ⁇ B k from splattering into the computed ephemeris corrections ⁇ r k .
- the estimators 146 may utilize the three dimensional approach for each GPS satellite that is observed by the primary WRS and at least one secondary WRS and utilize the four dimensional approach for each GPS satellite that is observed by at least one WRS but is not observed by the primary WRS and at least one secondary WRS.
- the WDGPS reference system 104 is capable of estimating the ephemeris and clock corrections for any GPS satellite in view of at least one of the WRSs 102, the user's GPS positioning unit 114 can make accurate position fixes using the ephemeris and clock corrections transmitted by the WMS 108.
- the user's positioning unit 114 includes a correction message receiver 152, a conventional L1 frequency GPS receiver 154, a data processor 158, and a display 160.
- the data processor includes a CPU 166 that runs the correction message unpacker 168, the pseudo-range correction estimator 170, the pseudo-range estimator 172, the ionospheric correction estimator 174, the tropospheric correction estimator 176, and the position estimator 178 which are all stored in the memory 180.
- the receiver 152 includes an antenna 162 for receiving the correction message transmitted by the WMS 108 and re-transmitted by the geostationary satellite 112.
- the GPS receiver 154 includes a conventional GPS antenna 164 for receiving L1 GPS signals from GPS satellites in its view. Alternatively, the correction message receiver 152 and the GPS receiver 154 may be integrated as a single unit for receiving both the correction message and the GPS signals.
- the GPS antenna 164 of the GPS receiver 154 receives an L1 GPS signal from each of the GPS satellites 106 that is in view and the GPS receiver 154 in response observes each of these GPS satellites. For each observed GPS satellite, the GPS receiver 154 computes at each epoch a pseudo-range ⁇ u k .sbsp.u from the user's positioning unit 114 to the corresponding GPS satellite based on the L1/CA code received from the GPS satellite, where k u represents the kth GPS satellite observed by the user's positioning unit.
- the GPS receiver 154 also computes a line of sight 1 u k .sbsp.u from the user's positioning unit to each observed GPS satellite
- the tropospheric correction estimator 176 computes a tropospheric correction T u k .sbsp.u for each of the GPS satellites 106 observed by the user's positioning unit in response to the line of sight 1 u k .sbsp.u to each GPS satellite. This is done according to the modified Hopfield method discussed earlier.
- the user's positioning system could include a meteorological sensor to sense surface pressure, humidity, and/or temperature.
- the antenna 162 of the receiver 152 receives the correction message transmitted by the WMS 108 and the receiver in response provides the correction message to the correction message unpacker 168.
- the ionospheric delay estimator computes an ionospheric correction l u k .sbsp.u for each GPS satellite observed by the user's positioning unit.
- the user's positioning system could use a dual frequency GPS receiver to receive L1 and L2 GPS signals.
- the dual frequency GPS receiver could compute the ionospheric corrections directly without requiring the grid of ionospheric corrections transmitted by the WMS 108 to compute these corrections.
- the pseudo-range correction estimator 172 computes at each epoch a pseudo-range correction ⁇ u k .sbsp.u for the corresponding GPS satellite according to:
- Each computed pseudo-range correction ⁇ u k .sbsp.u, term (t-t j ) ⁇ k , tropospheric delay estimate T u k .sbsp.u, and ionospheric correction l u k .sbsp.u is subtracted from the corresponding pseudo-range ⁇ u k .sbsp.u by the pseudo-range estimator 172 to produce the corrected pseudo-range ⁇ u k .sbsp.u given by:
- the position estimator computes in response a position fix for the user's positioning unit 114 according to well known techniques.
- the computed position fix is provided to the display 160 which then displays the position fix for the user.
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Security & Cryptography (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
A GPS reference system and corresponding method utilizing a network of reference stations and a master station. Each reference station receives GPS signals from GPS satellites it observes and computes in response a pseudo-range residual for each GPS satellite it observes. The master station is in communication with the reference stations to receive the computed pseudo-range residuals and comprises a pseudo-range residual synchronizer, an ephemeris and clock correction estimator for each GPS satellite observed by the reference stations, and a transmitter. The pseudo-range residual synchronizer is responsive to the received pseudo-range residuals in order to compute clock differences between the reference stations and remove the clock differences from the received pseudo-range residuals so as to synchronize them. Each ephemeris and clock correction estimator computes an ephemeris correction and a clock correction for the corresponding GPS satellite in response to each synchronized pseudo-range residual for the corresponding GPS satellite. The transmitter transmits to users the computed ephemeris and clock corrections for each GPS satellite observed by the reference stations.
Description
This invention was supported in part by grant 93-G-004 from the Federal Aviation Administration. The U.S. Government may therefore have rights to this invention.
This is a continuation, of application Ser. No. 08/374,165 filed Jan. 17, 1995 now abandoned.
The present invention relates generally to differential GPS reference systems and methods. More specifically, the present invention pertains to a wide area differential GPS (WDGPS) system and method that employs common view time transfer and minimum norm solution techniques to provide GPS satellite ephemeris and clock corrections to GPS users.
The global positioning system (GPS) is a satellite-based navigation and time transfer system developed by the U.S. Department of Defense. It serves marine, airborne and terrestrial users.
The system includes GPS satellites that each transmit L1 and L2 GPS signals. The L1 and L2 GPS signals comprise modulated carrier signals at the GPS L1 and L2 frequencies (1,575 and 1,227 MHz) which are each modulated by a coarse-acquisition pseudo-random (C/A) code, an encrypted precision pseudo-random (P) code, and navigation data. The navigation data includes the transmitting GPS satellite's orbit location and clock offset from GPS time.
Since the P code is encrypted, it is unavailable to most users. However, a user can determine the pseudo-range to each GPS satellite in view based on the time of arrival of the C/A code contained by the GPS signals transmitted by the GPS satellite. The user can compute a position fix based on the determined pseudo-ranges to the GPS satellites in its view and the locations and clock offsets for these GPS satellites contained in the navigation data of the GPS signals they transmit.
However, the U.S. Department of Defense engages in selective availability and deliberately misrepresents the true location and clock offsets of the GPS satellites in the navigation data of the GPS signals they transmit. Moreover, the ionosphere and the troposphere delay the arrival of these GPS signals at the user's location. Thus, these factors affect the accuracy and integrity of GPS for the user.
In order to improve the accuracy and integrity of GPS, GPS may be augmented by a network of ground-based reference stations. The reference stations monitor the health of the GPS signals transmitted by the GPS satellites in their view and generate corrections to these signals which are then transmitted to the user. The user can accurately determine his position based on the GPS signals he receives from the GPS satellites and the corrections to these signals he receives from the reference stations. This form of GPS is known as differential GPS (DGPS).
A number of systems and methods have been developed for wide area DGPS (WDGPS). However, all of them suffer from various problems which render them undesirable.
One of these approaches is called position domain WDGPS and combines the information from multiple wide area reference stations (WRSs) in the position domain. Such an approach is described in "Multi-Site Real-Time DGPS System Using Staffix Link: Operational Results", by D. Lapucha and M. Huff, Proceedings of the Fifth International Technical Meeting of the Satellite Division of the Institute of Navigation, Albuquerque, September 1992.
In this approach, each WRS acts as a local area DGPS (LDGPS) reference station. Based on its known fixed position, it computes for each GPS satellite in view a corresponding pseudo-range correction to the GPS satellite. Each WRS then transmits a stream of LDGPS pseudo-range corrections to the user. The user computes pseudo-ranges to the GPS satellites in his view based on the GPS signals received from these GPS satellites. The user then applies the pseudo-range corrections to the computed pseudo-ranges. If there are M WRSs, the user computes M position fixes and the corresponding covariance matrices. The user then uses the covariance matrices to form a single weighted position fix.
However, this approach does not involve the estimation of any of the underlying states of the GPS satellites, namely GPS satellite ephemeris and clock errors. Nor does it involve the estimation of the clock differences between the different reference stations. These type of errors behave quite differently and could be mitigated if they were separated. Since this approach allows these errors to be combined, it unfortunately does not allow a-priori information about each of the underlying states to be introduced into the estimation process. As a result, the accuracy of each computed pseudo-range correction decreases rapidly with age and user to reference station separation.
Additionally, this approach requires M independent data streams to be sent to the user if there are M different reference stations. Thus, it requires a rather large data bandwidth to implement.
Another approach is called measurement domain WDGPS and combines the pseudo-range measurements from multiple WRSs. Such a system is described in "The FAA's WIBANDGPS Testbed and Recent Test Results", by M. Lage and B. Elrod, Proceedings of the Sixth International Technical Meeting of the Satellite Division of the Institute of Navigation, Salt Lake City, 1993, and in "GPS Wide Area Augmentation System (WAAS) Testbed Results--Phase ID Testbed Results", by F. Haas, M. Lage, and S. Kalinowski, Proceedings of the Annual Meeting of the Institute of Navigation, Colorado Springs, June 1994.
In this approach, each WRS computes a pseudo-range correction to each GPS satellite in view. These pseudo-range corrections are sent to a wide area master station (WMS) rather than the user. The WMS uses well known common view time transfer techniques to estimate the clock differences between the WRSs. The estimated clock differences are then removed from the pseudo-range corrections. As a result, these synchronized pseudo-range corrections are all referenced to a common master clock. The WMS then averages for each GPS satellite the synchronized pseudo-range corrections for the GPS satellite. These averaged pseudo-range corrections are then transmitted to the user and used by the user in the same way as described earlier for the position domain WDGPS approach.
The measurement domain WDGPS approach requires less data bandwidth than the position domain WDGPS approach since a single averaged data stream for each GPS satellite is sent to the user rather than M different data streams. However, in this approach as well, the pseudo-range corrections degrade rapidly with age and user to reference station separation since the underlying states of the GPS satellites are not estimated and therefore the user cannot leverage the known behavior of the underlying states to make accurate position fixes. This frailty is especially pronounced if the ephemeris errors associated with one or more of the GPS satellites is large or the ionosphere is introducing errors which vary spatially.
Still another approach is called state space WDGPS which uses the information from multiple WRSs to estimate the underlying states of the GPS satellites, the underlying state of the ionosphere, and the underlying clock differences between the WRSs. Such a system is described in "Algorithms and Implementation of Wide Area Differential GPS", by C. Kee and B. Parkinson, Proceedings of the Fifth International Technical Meeting of the Satellite Division of the Institute of Navigation, Albuquerque, September 1992, and also in "Wide Area Differential GPS (WADGPS)", by C. Kee, Stanford University, Ph.D Dissertation, December 1993.
In this approach as well, each WRS computes a pseudo-range residual (i.e., correction) for each GPS satellite in view. The computed pseudo-range residuals are sent to a WMS which uses one large estimator to simultaneously compute the clock differences between the WRSs and the ephemeris and clock corrections for the ephemeris and clock errors of each GPS satellite in view of the network.
The estimator in this approach must solve a system of equations which contains 4K+(M-1) unknowns and MK knowns, where K is the number of satellites in view and M is the number of reference stations. This stems from the fact that each satellite has 3 unknown components of ephemeris error and 1 unknown clock error. In addition, there are M-1 unknown offsets between the clocks of the WRSs clocks. At the same time, the system of equations contains MK knowns corresponding to the K pseudo-range residuals computed by each of the M WRSs.
When 4K+(M-1)>MK, the system is under determined and the WMS must use a minimum norm solution. When 4K+(M-1)=MK, then the system is exactly specified and can be readily inverted. When 4K+(M-1)<MK, then the system is over specified and the pseudo-inverse yields an estimate with minimum mean square error.
Unlike the previously discussed approaches, the state space WDGPS approach separates the underlying states of the system, and as such, it enables the use of a-priori models and side information for these states. However, its complexity impedes the straight forward integration of such models.
Moreover, it requires a greater data bandwidth than is necessary. This occurs, because the WRS network can be rather sparse and, as such, the system of equations used by the WMS is frequently under determined. Specifically, rising GPS satellites that have just appeared over the horizon are of great concern to WDGPS users because they can greatly improve the geometry and accuracy of the position fix. However, the system is always under determined in the case of a rising GPS satellite. Even though the minimum norm solution described above yields an accurate correction for a newly rising GPS satellite, the quickly varying components of the GPS satellite clock correction splatter into the GPS satellite ephemeris correction. Consequently, the ephemeris correction must be sent as frequently as the clock correction, thereby increasing the required data bandwidth.
A similar approach is described in U.S. Pat. No. 5,323,322 issued on Jun. 21, 1994 to Mueller, et. al. In this approach, ephemeris and clock corrections for a particular GPS satellite can be provided to a user only if the GPS satellite is observed by five or more WRSs. However, as indicated earlier, rising GPS satellites are very important to users since they provide added leverage for computing accurate position fixes. Thus, this approach is not capable of providing a user with ephemeris and clock corrections for newly rising GPS satellites.
The foregoing problems are solved by a GPS reference system and corresponding method utilizing a network of reference stations and a master station. Each reference station receives GPS signals from GPS satellites it observes and computes in response a pseudo-range residual for each GPS satellite it observes. The master station is in communication with the reference stations to receive the computed pseudo-range residuals and comprises a pseudo-range residual synchronizer, an ephemeris and clock correction estimator for each GPS satellite observed by the reference stations, and a transmitter. The pseudo-range residual synchronizer is responsive to the received pseudo-range residuals in order to compute clock differences between the reference stations and remove the clock differences from the received pseudo-range residuals so as to synchronize them. Each ephemeris and clock correction estimator computes an ephemeris correction and a clock correction for the corresponding GPS satellite in response to each synchronized pseudo-range residual for that GPS satellite. The transmitter transmits to users the computed ephemeris and clock corrections for each GPS satellite observed by the reference stations.
FIG. 1 shows a wide area differential GPS (WDGPS) system in accordance with the present invention.
FIG. 2 provides a block diagram of a wide area reference station (WRS) used in the WDGPS system of FIG. 1.
FIG. 3 shows the operation of the WRS of FIG. 2.
FIG. 4 provides a block diagram of a wide area master station (WMS) used in the WDGPS system of FIG. 1.
FIG. 5 shows the operation of the WMS of FIG. 4.
FIG. 6 demonstrates the situation where a GPS satellite is observed by only three WRSs in the WDGPS system of FIG. 1.
FIG. 7 demonstrates the situation where a GPS satellite is observed by only two WRSs in the WDGPS system of FIG. 1.
FIG. 8 demonstrates the situation where a GPS satellite is observed by only one WRSs in the WDGPS system of FIG. 1.
FIG. 9 provides a block diagram of a user's GPS positioning unit used in the WDGPS system of FIG. 1.
FIG. 10 shows the operation of the user's GPS positioning unit of FIG. 9.
Referring initially to FIG. 1, there is shown a wide area differential GPS (WDGPS) network 100 that includes a WDGPS reference or augmentation system 104. The WDGPS reference system includes a network of wide area reference stations (WRSs) 102. Each WRS receives L1 and L2 GPS signals from each of the GPS satellites 106 that are in its view. In response, each WRS observes each of the GPS satellites in its view and computes a pseudo-range residual and an ionospheric delay estimate for each observed GPS satellite. The pseudo-range residuals and ionospheric delay estimates computed by each WRS are provided to the wide area master station (WMS) 108 of the WDGPS reference system via a data communication link 110. The WMS in response computes ephemeris and clock corrections for each GPS satellite observed by the WRSs and computes a grid of ionospheric corrections. The WMS then transmits the ephemeris, clock, and ionospheric corrections to the geostationary satellite 112. The geostationary satellite 112 re-transmits these corrections so that a mobile user's GPS positioning system 114 receives them along with the L1 GPS signals from each GPS satellite it observes. The user's GPS positioning system then computes its position based on the received GPS signals and the received ephemeris, clock, and ionospheric corrections.
As shown in FIG. 2, each WRS 102 includes a conventional survey quality dual frequency GPS receiver 116, a meteorological sensor 118, a data processor 120, and a data link interface 121. The GPS receiver includes a conventional GPS antenna 122 for receiving L 1 and L2 GPS signals. The data processor includes a central processing unit (CPU) 124 that runs the pseudo-range smoother 126, the navigation data processor 128, the tropospheric correction estimator 130, the pseudo-range residual estimator 132, and the clock steerer 134 which are all stored in the memory 136. The memory also stores the precisely surveyed location 138 of the WRS.
Turning to FIGS. 1 and 3, the GPS antenna 122 of each WRS 102 receives L1 and L2 GPS signals from each of the GPS satellites 106 in view of the WRS and the GPS receiver 116 in response observes each of these GPS satellites. For each observed GPS satellite, the GPS receiver computes at each of its measurement epochs L1 and L2 pseudo-ranges ρm,L1k.sbsp.m and ρm,L2k.sbsp.m (from the WRS to the corresponding GPS satellite) based on the L1/CA and L2/CA codes received from the GPS satellite, where km represents the kth GPS satellite observed by the mth WRS. In addition, the GPS receiver makes at each epoch L1 and L2 carrier phase measurements φm,L1k.sbsp.m and φm,L2k.sbsp.m of the L1 and L2 carrier signals transmitted by each GPS satellite it observes. In the preferred embodiment, the measurement epochs occur at one second intervals (i.e., 1 Hz).
The sets of L1 and L2 pseudo-ranges {ρm,L1k.sbsp.m }k.sbsb.m=1 K.sbsp.m and {ρm,L2k.sbsp.m }k.sbsb.m=1 K.sbsp.m and the sets of L1 and L2 carrier phase measurements {φm,L1k.sbsp.m }k.sbsb.m=1 K.sbsp.m and {φm,L2k.sbsp.m }k.sbsb.m=1 K.sbsp.m at each epoch are then provided to the pseudo-range smoother 126 of the WRS, where Km represents the total number of GPS satellites observed by the mth WRS. For each GPS satellite observed by the WRS, the pseudo-range smoother uses the corresponding L1 and L2 carrier phase measurements and the corresponding L1 and L2 pseudo-ranges to produce at each epoch a corresponding ionospheric delay estimate lm k.sbsp.m and a corresponding smoothed pseudo-range ρm k.sbsp.m from which the ionospheric delay estimate has been removed. This is done using the method described in "GPS Survey", by C. Goad, a chapter in "The Global Positioning System: Theory and Applications", edited by B. Parkinson, J. Spilker, P. Axelrod, and P. Enge, American Institute of Aeronautics and Astronautics, 1995, which is hereby explicitly incorporated by reference. The set of smoothed pseudo-ranges {ρm k.sbsp.m }k.sbsb.m=1 K.sbsp.m at each epoch are sent to the pseudo-range residual estimator 132 and the set of ionospheric delay estimates {m k.sbsp.m }k.sbsb.m=1 k.sbsp.m at each epoch are sent to the WMS 108.
The GPS receiver 116 of each WRS 102 also extracts the navigation data contained by the L1 and L2 GPS signals it receives from each GPS satellite it observes. The navigation data processor 128 is provided with the navigation data for each observed GPS satellite and unpacks from it the GPS satellite's location vector rk.sbsp.m.sup.,B and the difference (or offset) Bk.sbsp.m.sup.,B from GPS time of the GPS satellite's clock. Then, for each GPS satellite observed by the WRS, the navigation data processor computes at each epoch a line of sight vector 1m k.sbsp.m and a nominal range rk.sbsp.m.sup.,B ·1m k.sbsp.m from the WRS to the GPS satellite based on the unpacked GPS satellite's location and the known reference station location 138 stored in the memory 130. The set of line of sights {1m k.sbsp.m }k.sbsb.m=1 K.sbsp.m are provided to the WMS 108 while the set of computed nominal ranges {rk,B ·1m k.sbsp.m }k.sbsb.m=1 K.sbsp.m and the set of unpacked clock differences {Bk.sbsp.m,B }k.sbsb.m=1 K.sbsp.m are provided to the pseudo-range residual estimator 132.
The conventional meteorological sensor 118 of each WRS 102 senses surface pressure, humidity, and/or temperature and provides this sensed meteorological data to the tropospheric correction estimator 130. The tropospheric delay estimator 130 periodically computes at each epoch a tropospheric delay estimate Tm k.sbsp.m for each of the GPS satellites 106 observed by the WRS in response to the sensed meteorological data and the line of sight to each GPS satellite. This is done using a modified Hopfield method as described in "Analysis of Tropospheric Delay Prediction Models: Comparisons with Ray-Tracing and Implications for GPS Relative Positioning", by H. Janes, R Langley, and S. Newby, Bulletin Geodesique, Vol 65, pp 151-161, 1991, which is hereby incorporated by reference. The set of tropospheric delay estimates {Tm k.sbsp.m }k.sbsb.m=1 K.sbsp.m are then provided to the pseudo-range residual estimator 132.
Each smoothed pseudo-range ρm k.sbsp.m received at each epoch by the pseudo-range residual estimator 132 may be modeled as:
ρ.sub.m.sup.k.sbsp.m =r.sup.k.sbsp.m ·1.sub.m.sup.k.sbsp.m +b.sub.m -B.sup.k.sbsp.m +T.sub.m.sup.k.sbsp.m +ν.sub.m.sup.k.sbsp.mEq. (1)
where rk.sbsp.m is a vector that represents the actual location of the kth GPS satellite 106 observed by the WRS, bm represents the actual difference from GPS time of the clock of the GPS receiver 116 of the mth WRS 102, Bk.sbsp.m represents the actual difference from GPS time of the clock of the kth GPS satellite observed by the mth WRS, and νm k.sbsp.m represents the noise associated with the kth GPS satellite and the mth WRS that observes it. For each smoothed pseudo-range ρm k.sbsp.m, the pseudo-range residual estimator 132 subtracts the corresponding nominal range rk.sbsp.m.sup.,B ·1m k.sbsp.m, unpacked clock offset Bk.sbsp.m.sup.,B, and tropospheric correction Tm k.sbsp.m to produce a pseudo-range residual (i.e. correction) Δρm k.sbsp.m at each epoch for the error in the smoothed pseudo-range. This pseudo-range residual may be represented as:
Δρ.sub.m.sup.k.sbsp.m =Δr.sup.k.sbsp.m ·1.sub.m.sup.k.sbsp.m +b.sub.m -ΔB.sup.k.sbsp.m +ν.sub.m.sup.K.sbsp.m Eq. (2)
where Δrk.sbsp.m =rk.sbsp.m -rk.sbsp.m.sup.,B and is a vector representing the correction for the ephemeris error of the kth GPS satellite observed by the mth WRS, and ΔBk.sbsp.m =Bk.sbsp.m -Bk.sbsp.m.sup.,B and represents the correction for the clock error of the kth GPS satellite observed by the mth WRS.
Since the clock difference bm dominates the expression in Equation (2), the reference station clock steerer 134 averages at each epoch all of the pseudo-range residuals {Δρm k.sbsp.m }k.sbsb.m=1 K.sbsp.m for the GPS satellites observed by the WRS to form a rough estimate bm of the actual clock difference bm, as is well known by those skilled in the art. The clock difference estimate bm is subtracted by the pseudo-range residual estimator 132 from each pseudo-range residual Δρm k.sbsp.m to produce an adjusted pseudo-range residual Δρm k.sbsp.m given by:
Δρ.sub.m.sup.k.sbsp.m =Δr.sup.k.sbsp.m ·1.sub.m.sup.k.sbsp.m +Δb.sub.m -ΔB.sup.k.sbsp.m +ν.sub.m.sup.k.sbsp.m Eq.(3)
where Δbm =bm -bm. The set of pseudo-range residuals {Δρm k.sbsp.m }k.sbsb.m=1 K.sbsp.m are then transmitted to the WMS 108 over one of the data communication links 110 by the data link interface 121. Each data communication link 110 may comprise a conventional modem/telephone line system, an RF or other frequency transmitter/receiver system, a fiber optic cable network, or other suitable transmission system.
Referring to FIG. 4, the WMS 108 includes data link interfaces 140, a data processor 142, and a transmitter 143. The data processor includes a pseudo-range residual synchronizer 144, ephemeris and clock correction estimators 146, and an ionospheric grid generator 148 which are run by the CPU 147 and stored in the memory 150.
Turning to FIGS. 1 and 5, the WMS 108 includes the same number of data link interfaces 140 as there are WRSs 102 in the WDGPS reference system 104. As a result, the interfaces 140 provide the pseudo-range residual synchronizer 144 at each epoch with the pseudo-range residuals {{Δρm k.sbsp.m }k.sbsb.m=1 K.sbsp.m }m=1 M, provide the ionospheric grid generator 148 with the ionospheric delay estimates {{lm k.sbsp.m }k.sbsb.m=1 K.sbsp.m }m=1 M, and provide the ephemeris and clock correction estimators 146 with the lines of sight {{1m k.sbsp.m }k.sbsb.m=1 K.sbsp.m }m=1 M (re-phrased as {{1m.sbsb.kk }k=1 K }m.sbsb.k=1 M.sbsp.k, as will be discussed later), where M represents the total number of WRSs.
As those skilled in the art will appreciate, the lines of sight {{1m k.sbsp.m }k.sbsb.m=1 K.sbsp.m }m=1 M could be computed by the WMS 108 rather than the WRSs 102. In this case, the WMS may include a GPS receiver that receives the GPS signals transmitted by the GPS satellites 104 and then removes the navigation data contained in the GPS signals. Or, the navigation data in the GPS signals received by the WRSs could be simply provided to the WMS via the communication links 110. In either case, based on the GPS satellite locations contained in the navigation data of the GPS signals and the known locations of the WRSs which would be stored in the memory 150 of the WMS, the WMS would then compute the lines of sight.
In response to the ionospheric delay estimates {{lm k.sbsp.m }k.sbsb.m=1 K.sbsp.m }m=1 M, the ionospheric grid generator 148 produces a grid of ionospheric corrections {ll }l=1 lat,lon over the area covered by the GPS satellites 106. This is done in the manner described in "Flight Trials of the Wide-Area Augmentation System (WAAS)", by Walter et. al., Proceedings of the International Technical Meeting of the Satellite Division of the Institute of Navigation, Salt Lake City, September 1994, which is hereby explicitly incorporated by reference. The grid of ionospheric corrections is provided to the user's positioning system 114 in a correction message transmitted by the transmitter 143.
In response to the pseudo-range residuals {{Δρm k.sbsp.m }k.sbsb.m=1 K.sbsp.m }m=1 M, the pseudo-range synchronizer 144 first performs at each epoch a common view time transfer to estimate the differences between the clocks of the WRSs 102 with respect to a master clock. The clock differences are then removed from the pseudo-range residuals so as to synchronize them.
In order to do so, one of the WRSs is designated as a primary WRS and the remaining WRSs are designated as secondary WRSs. For each GPS satellite 106 in common view of the primary WRS and at least one of the secondary WRSs, the pseudo-range residual synchronizer 144 first subtracts at each epoch the pseudo-range residual ΔρM k.sbsp.m,M computed for it by the primary WRS from each pseudo-range residual Δρm k.sbsp.m,M computed for it by those of the secondary WRSs that observe it, where M represents the primary WRS (in addition to representing the number of WRSs) and, for m≠M,km,M represents the kth of the GPS satellites that is observed by the mth secondary WRS and by the primary WRS. This results in the difference Δm k.sbsp.m,M at each epoch which is given as follows:
Δ.sub.m.sup.k.sbsp.m,M =Δρ.sub.M.sup.k.sbsp.m,M -Δρ.sub.M.sup.k.sbsp.m,M =Δr.sup.k.sbsp.m,M ·(1.sub.m.sup.k.sbsp.m,M -1.sub.M.sup.k.sbsp.m,M)+Δb.sub.m -Δb.sub.M +ν.sub.m,M.sup.k.sbsp.m,M Eq.(4)
where νm,Mk.sbsp.m,M =νm k.sbsp.m,M -νM k.sbsp.m,M. Assuming that the difference in line of sights (1m k.sbsp.m,M -1M k.sbsp.m,M) is small and/or the ephemeris correction Δrk.sbsp.m,M is small, then an estimate of the WRS clock error can be made by averaging the difference Δm k.sbsp.m,M as follows: ##EQU1## where Δbm -ΔbM represents the clock difference at each epoch between the primary WRS and the mth secondary WRS and Km,M represents the total number of GPS satellites in common view of the primary WRS and the mth WRS. Moreover, if the clocks of the GPS receivers 116 of the WRSs are high quality clocks, such as rubidium clocks, then the difference Δbm -ΔbM can be averaged over time to take advantage of the well behaved drift of the clocks. In this case, the expression in Equation (5) can be re-written as: ##EQU2## where h(tj) represents the well known impulse response for a low pass filter.
From the foregoing, it is evident that in order to compute a clock difference Δbm -ΔbM at each epoch for each WRS 102, the secondary WRSs must be positioned with respect to the primary WRS so that each secondary WRS will always have at least one GPS satellite in common view with the primary WRS. Preferably, the WRS that is closest to the WMS and most centrally located with respect to the more peripheral WRSs is the primary WRS.
The pseudo-range synchronizer 144 then at each epoch subtracts from the pseudo-range residuals {Δρm k.sbsp.m }k.sbsb.m=1 K.sbsp.m computed by each of the secondary WRSs the corresponding clock difference Δbm -ΔbM to produce the pseudo-range residuals {Δρm k.sbsp.m }k.sbsb.m=1 K.sbsp.m at each epoch. Thus, the pseudo-range residuals computed by the secondary WRSs {{Δρm k.sbsp.m }k.sbsb.m=1 K.sbsp.m are synchronized with the pseudo-range residuals {ΔρM k.sbsp.m }k.sbsb.M=1 K.sbsp.M computed by the primary WRS and referenced to the same master clock (i.e., the clock of the primary WRS). As a result, each synchronized pseudo-range residual can now be represented as Δρm k.sbsp.m and given by:
Δρ.sub.m.sup.k.sbsp.m =Δr.sup.k.sbsp.m ·1.sub.m.sup.k.sbsp.m +ΔB.sup.k.sbsp.m +ν.sub.m.sup.k.sbsp.m Eq.(7)
for m=1 to M, and where ΔρM k.sbsp.m =ΔρM k.sbsp.m, and ΔBk.sbsp.m =ΔbM -ΔBk.sbsp.m and represents the clock correction for the kth GPS satellite observed by the mth WRS synchronized to the clock of the primary WRS.
Furthermore, Equation (7) can be re-written as follows:
Δρ.sub.m.sbsb.k.sup.k =Δr.sup.k ·1.sub.m.sbsb.k.sup.k +ΔB.sup.k +ν.sub.m.sbsb.k.sup.kEq. (8)
where k represents the kth of all of the GPS satellites 106 observed by the entire network of WRSs 102 and mk represents the mth of those of the WRSs that observe the kth GPS satellite. The set of synchronized pseudo-range residuals {{Δρm.sbsb.kk }k=1 K }m.sbsb.k=1 M.sbsp.k and the corresponding set of lines of sight {{1m.sbsb.kk }k=1 K }m.sbsb.k=1 M.sbsp.k (re-phrased to correspond with the synchronized pseudo-range residuals) are then provided to the ephemeris and clock correction estimators 146 of the WMS 108, where Mk represents the total number of WRS that observe the kth GPS satellite.
Since the clock differences of the WRSs 102 have been removed from the set of synchronized pseudo-range residuals {{Δρm.sbsb.kk }k=1 K }m.sbsb.k=1 M.sbsp.k, the WMS 108 can employ an ephemeris and clock correction estimator 146 for each of the GPS satellites 106 observed by the WRSs. Each ephemeris and clock correction estimator then computes at each epoch the ephemeris correction vector Δrk and the clock correction ΔBk for the corresponding GPS satellite in response to the lines of sight {1m.sbsb.kk }m.sbsb.k=1 M.sbsp.k to the corresponding GPS satellite and the synchronized pseudo-range residuals {Δρm.sbsb.kk }m.sbsb.k=1 M.sbsp.k for the corresponding GPS satellite, independent of the synchronized pseudo-range residuals for the other GPS satellites. As will be evident from the following discussion, by employing separate ephemeris and clock correction estimators 146, the computational complexity of a single large estimator is avoided.
In one embodiment, each ephemeris and clock correction estimator 146 solves for the four dimensional state vector [Δrk,ΔBk ]T based on the following matrix system: ##EQU3##
The system of Equation (9) can be more compactly re-written as follows:
Δρ.sup.k =G.sub.W X.sup.k Eq. (10)
where Δρk is the vector containing each synchronized pseudo-range residual Δρm k computed by the WRSs 102 for the kth GPS satellite 106, xk is the four dimensional state vector containing the ephemeris and clock corrections Δrk and ΔBk for the kth GPS satellite, and GW is the design matrix that contains, for each WRS that observes the kth GPS satellite, a line of sight from it to the kth GPS satellite and a matrix element (i.e., -1) that indicates that the corresponding synchronized pseudo-range residual in the vector Δρk should be adjusted by the clock correction ΔBk for the kth GPS satellite.
When the WDGPS reference system 104 includes at least four WRSs 102 (i.e., the primary WRS and at least three secondary WRSs) and the kth GPS satellite 106 is observed by exactly four of these WRSs, then the system of Equation (10) is exactly specified. Furthermore, when the WDGPS reference system 104 includes more than four WRSs 102 and the kth GPS satellite 106 is observed by more than four of these WRSs, then the system of Equation (10) is over specified. In either case, the corresponding estimator 146 first computes a pseudo-inverse (GW T GW)-1 GW T of the design matrix Gw. In the case where the corresponding GPS satellite is observed by exactly four of the WRSs, the pseudo-inverse reduces to the inverse GW -1 of the design matrix GW. The estimator then computes a minimum mean square error solution for the state vector xk as follows:
X.sup.k =(G.sub.W.sup.T G.sub.W).sup.-1 G.sub.W.sup.T Δρ.sup.kEq. (11)
Moreover, the system of Equation (10) is under specified when the kth GPS satellite 106 is observed by less than four WRSs 102, regardless of whether the WDGPS reference system 104 includes less than four WRSs (i.e., the primary WRS and no more than two secondary WRSs) or at least four WRSs. In this situation, the corresponding estimator 146 computes a pseudo-inverse GW T (GW GW T)-1 of the design matrix GW and then computes a minimum norm solution for the state vector xk as follows:
x.sup.k =G.sub.W.sup.T (G.sub.W G.sub.W.sup.T).sup.-1 Δρ.sup.kEq. (12)
The minimum norm solution has the minimum length or energy among all of the satisfactory solutions.
The set of computed ephemeris and clock corrections {Δrk ; ΔBk }k=1 K for the GPS satellites observed by the WRSs 102 are then provided to the transmitter 144. The transmitter transmits them as corrections to the user's GPS positioning system 114 via the geostationary satellite 112 or other suitable broadcast facility so that the user's positioning system can compute an accurate position fix. Alternatively, as those skilled in the art will appreciate, the computed ephemeris and clock corrections may be transmitted directly to the user's positioning system by the transmitter 143.
FIGS. 6 illustrates the situations when the kth GPS satellite 106 is observed by only three WRSs 102 and the system of Equation (10) is under specified. Since the kth GPS satellite 106 is observed by only three WRSs, this geometry does not provide leverage for precisely observing each of the radial (i.e., orbit radius), along track (along satellite velocity), and cross track (i.e., orthogonal to radial and along track) components of the ephemeris error of the GPS satellite and the clock error of the GPS satellite. However, when the user's positioning system 114 is between the three WRSs or adjacent to one of the three WRSs, the minimum norm solution of Equation (12) does provide accurate ephemeris and clock corrections Δrk and ΔBk for the ephemeris and clock errors since the components of these errors that effect the user's positioning system effect the three WRSs substantially in the same manner.
FIG. 7 demonstrates the case when the kth GPS satellite 106 is observed by only two WRSs 102 and the system of Equation (10) is under specified. Similar to the situation described for FIG. 6, although this geometry does not provide leverage for precisely observing each of the components of the ephemeris and clock errors, the minimum norm solution of Equation (12) does provide accurate ephemeris and clock corrections Δrk and ΔBk for when the user's positioning system 114 is between the two WRSs or adjacent to one of the WRSs. As in the situation of FIG. 6, this is due to the fact that the components of the ephemeris and clock errors that effect the user's positioning system effect the two WRSs substantially in the same manner.
Similarly, FIG. 8 shows the situation when the kth GPS satellite 106 is observed by only one WRS 102 and the system of Equation (10) is under specified. This geometry also does not provide the leverage for accurately observing the components of the ephemeris and clock errors. But, when the user's positioning system 114 is adjacent to the WRS, the minimum norm solution of Equation (12) provides accurate ephemeris and clock corrections Δrk and ΔBk because the components of the ephemeris and clock errors that effect the user's positioning system effect the WRS in substantially the same manner.
From the foregoing, it is evident that the WDGPS reference system 104 is capable of estimating applicable ephemeris and clock errors for any GPS satellite in view of at least one of the WRSs 102, provided that the user's positioning system 114 is between WRSs or adjacent to a WRSs. As result, the WDGPS reference system can provide to users accurate and valid ephemeris and clock corrections for each GPS satellite from the time that it is first observed by one of the WRSs (i.e., appears on the horizon) until the time that it is no longer observed by any of the WRSs (i.e., disappears below the horizon).
However, in this embodiment, even though the four dimensional minimum norm solution for the state vector xk always generates corrections which are valid for users in between or adjacent the WRSs 102, the high frequency energy in the clock correction ΔBk (resulting from the high frequency energy in the clock error due to selective availability) splatters into the ephemeris corrections Δrk. As a consequence, the ephemeris corrections {Δrk }k=1 K must be transmitted by the WMS 108 in the correction message to users as frequently (i.e., at each epoch) as the clock corrections {ΔBk }k=1 K. This may be undesirable in an operational WDGPS reference system where the broadcast bandwidth is severely limited and only allows the clock corrections to be sent frequently. In such a case, three dimensional ephemeris and clock correction estimators may be used instead.
In this embodiment, each ephemeris and clock correction estimator 146 first subtracts each synchronized pseudo-range residual Δρm.sbsb.kk for the corresponding GPS satellite computed by the secondary WRSs from the synchronized pseudo-range residual ΔρM.sbsb.kk for the corresponding GPS satellite computed by the primary WRS, where Mk also represents the primary WRS (in addition to representing the total number of WRSs that observe the kth GPS satellite). Each difference Δm.sbsb.kM.sbsb.kk can be expressed as:
Δ.sub.m.sbsb.k.sub.,M.sbsb.k.sup.k =Δρ.sub.m.sbsb.k.sup.k -Δρ.sub.M.sbsb.k.sup.k =Δr.sub.m.sbsb.k.sub.,M.sbsb.k.sup.k ·(1.sub.m.sbsb.k.sup.k -1.sub.M.sbsb.k.sup.k)+ν.sub.m.sbsb.k.sub.,M.sbsb.k.sup.kEq. (13)
where νm.sbsb.k.sub.,M.sbsb.kk =νm.sbsb.kk -νM.sbsb.kk. This is done in order to separate the computation of the clock correction ΔBk from the computation of the ephemeris correction Δrk so that the high frequency energy of the clock correction does not splatter into the ephemeris correction.
Each estimator then solves the three dimensional ephemeris correction vector Δrk based on the following matrix system: ##EQU4## Similar to the system of Equation (9), the system of Equation (14) can be more compactly re-written as:
Δ.sup.k =G.sub.W Δr.sup.k Eq. (15)
where Δk is a vector containing each computed difference Δm.sbsb.k.sub.,M.sbsb.kk and GW is the matrix that contains, for each secondary WRS that observes the kth GPS satellite, the difference between the line of sight 1m.sbsb.kk from the corresponding secondary reference station to the kth GPS satellite and the line of sight 1M.sbsb.kk from the primary reference station to the kth GPS satellite.
As with the system of Equation (10), the system of Equation (15) is exactly specified when the WDGPS reference system 104 includes at least four WRSs 102 (i.e., the primary WRS and at least three secondary WRSs) and the kth GPS satellite 106 is observed by exactly four of these WRSs. And, the system of Equation (15) is over specified when the WDGPS reference system includes more than four WRSs 102 and when the kth GPS satellite 106 is observed by more than four of these WRSs. The corresponding estimator 146 in this case computes a three dimensional minimum mean square error solution for the state vector Δrk as follows:
Δr.sup.k =(G.sub.W.sup.T G.sub.W).sup.-1 G.sub.W.sup.T Δ.sup.kEq.(16)
In addition, the system of Equation (15) is under specified if the kth GPS satellite 106 is observed by less than four WRSs 102 in either situation where the WDGPS reference system 104 includes less than four WRS (i.e., the primary WRS and no more than two secondary WRSs) or at least four WRSs. In this case, the corresponding estimator 146 computes a three dimensional minimum norm solution for the state vector Δrk as follows:
Δr.sup.k =G.sub.W.sup.T (G.sub.W G.sub.W.sup.T).sup.-1 Δ.sup.kEq. (17)
For each of the situations described above, after computing the ephemeris correction Δrk, each estimator 146 then computes the clock correction ΔBk according to: ##EQU5##
The set of ephemeris and clock corrections {Δrk ;ΔBk }k=1 K computed by the estimators 146 are then provided to the transmitter 144 for transmission as ephemeris and clock corrections.
As was alluded to earlier, the three dimensional estimators 146 prevent the high frequency energy in the clock correction ΔBk from splattering into the computed ephemeris corrections Δrk. As a result, the set of computed ephemeris corrections {Δrk }k=1 K need not be transmitted as often as the set of computed clock corrections {ΔBk }k=1 K.
However, since newly rising GPS satellites are important to GPS users, a disadvantage to this three dimensional approach is that corrections for newly rising GPS satellites cannot be computed until the GPS satellite is viewed by the primary WRS and at least one secondary WRS. For this reason, the estimators 146 may utilize the three dimensional approach for each GPS satellite that is observed by the primary WRS and at least one secondary WRS and utilize the four dimensional approach for each GPS satellite that is observed by at least one WRS but is not observed by the primary WRS and at least one secondary WRS.
Referring to FIG. 1, since the WDGPS reference system 104 is capable of estimating the ephemeris and clock corrections for any GPS satellite in view of at least one of the WRSs 102, the user's GPS positioning unit 114 can make accurate position fixes using the ephemeris and clock corrections transmitted by the WMS 108.
As shown in FIG. 9, the user's positioning unit 114 includes a correction message receiver 152, a conventional L1 frequency GPS receiver 154, a data processor 158, and a display 160. The data processor includes a CPU 166 that runs the correction message unpacker 168, the pseudo-range correction estimator 170, the pseudo-range estimator 172, the ionospheric correction estimator 174, the tropospheric correction estimator 176, and the position estimator 178 which are all stored in the memory 180. The receiver 152 includes an antenna 162 for receiving the correction message transmitted by the WMS 108 and re-transmitted by the geostationary satellite 112. The GPS receiver 154 includes a conventional GPS antenna 164 for receiving L1 GPS signals from GPS satellites in its view. Alternatively, the correction message receiver 152 and the GPS receiver 154 may be integrated as a single unit for receiving both the correction message and the GPS signals.
Turning to FIGS. 1 and 10, the GPS antenna 164 of the GPS receiver 154 receives an L1 GPS signal from each of the GPS satellites 106 that is in view and the GPS receiver 154 in response observes each of these GPS satellites. For each observed GPS satellite, the GPS receiver 154 computes at each epoch a pseudo-range ρu k.sbsp.u from the user's positioning unit 114 to the corresponding GPS satellite based on the L1/CA code received from the GPS satellite, where ku represents the kth GPS satellite observed by the user's positioning unit. At each epoch, the GPS receiver 154 also computes a line of sight 1u k.sbsp.u from the user's positioning unit to each observed GPS satellite The set of pseudo-ranges {ρu k.sbsp.u }k.sbsb.u=1 K.sbsp.u at each epoch are then provided to the pseudo-range estimator 172 of the user's positioning unit and the set of lines of sight {1u k.sbsp.u }k.sbsb.u=1 K.sbsp.u at each epoch are provided to the tropospheric correction estimator 176, where Ku represents the total number of GPS satellites observed by the user's positioning unit.
The tropospheric correction estimator 176 computes a tropospheric correction Tu k.sbsp.u for each of the GPS satellites 106 observed by the user's positioning unit in response to the line of sight 1u k.sbsp.u to each GPS satellite. This is done according to the modified Hopfield method discussed earlier. The set of tropospheric corrections {Tu k.sbsp.u }k.sbsb.u=1 K.sbsp.u are then provided to the pseudo-range estimator 172.
Alternatively, the user's positioning system could include a meteorological sensor to sense surface pressure, humidity, and/or temperature. In response to this meteorological data, the tropospheric correction estimator 176 would compute the tropospheric corrections {Tu k.sbsp.u }k.sbsb.u=1 K.sbsp.u.
The antenna 162 of the receiver 152 receives the correction message transmitted by the WMS 108 and the receiver in response provides the correction message to the correction message unpacker 168. The correction message unpacker unpacks the sets of ephemeris and clock corrections {Δrk ; ΔBk }k=1 K and sends them to the pseudo-range estimator 172. In addition, it unpacks the grid of ionospheric corrections {ll }l=1 lat,lon and sends them to the ionospheric delay estimator 174.
In response to the grid of ionospheric corrections {ll }l=1 lat,lon, the ionospheric delay estimator computes an ionospheric correction lu k.sbsp.u for each GPS satellite observed by the user's positioning unit. The set of ionospheric corrections {lu k.sbsp.u }k.sbsb.u=1 K.sbsp.u are provided to the pseudo-range estimator 172.
Alternatively, the user's positioning system could use a dual frequency GPS receiver to receive L1 and L2 GPS signals. In this case, the dual frequency GPS receiver could compute the ionospheric corrections directly without requiring the grid of ionospheric corrections transmitted by the WMS 108 to compute these corrections.
In response to the unpacked ephemeris and clock corrections Δrk.sbsp.u and ΔBk.sbsp.u for each GPS satellite observed by the user's positioning system, the pseudo-range correction estimator 172 computes at each epoch a pseudo-range correction Δρu k.sbsp.u for the corresponding GPS satellite according to:
Δρ.sub.u.sup.k.sbsp.u =Δr.sub.u.sup.k.sbsp.u ·1.sub.u.sup.k.sbsp.u -ΔB.sup.k.sbsp.u Eq. (19)
and computes at each epoch the term (t-tj)Δρu k.sbsp.u, where Δρu k.sbsp.u is the pseudo-range correction velocity for the corresponding GPS satellite computed by differencing the computed pseudo-range corrections from the last two epochs as follows: ##EQU6##
Each computed pseudo-range correction Δρu k.sbsp.u, term (t-tj) Δρk, tropospheric delay estimate Tu k.sbsp.u, and ionospheric correction lu k.sbsp.u is subtracted from the corresponding pseudo-range ρu k.sbsp.u by the pseudo-range estimator 172 to produce the corrected pseudo-range ρu k.sbsp.u given by:
ρ.sub.u.sup.k.sbsp.u =ρ.sub.u.sup.k.sbsp.u -Δρ.sub.u.sup.k.sbsp.u -(t-t.sub.j)Δρ.sub.u.sup.k.sbsp.u -l.sub.u.sup.k.sbsp.u -T.sub.u.sup.k.sbsp.u Eq. (21)
The set of corrected pseudo-ranges {ρu k.sbsp.u }k.sbsb.u=1 K.sbsp.u are then provided to the position estimator 178. The position estimator computes in response a position fix for the user's positioning unit 114 according to well known techniques. The computed position fix is provided to the display 160 which then displays the position fix for the user.
While the present invention has been described with reference to a few specific embodiments, the description is illustrative of the invention and is not to be construed as limiting the invention. Various modifications may occur to those skilled in the art without departing from the true spirit and scope of the invention as defined by the appended claims.
Claims (15)
1. A GPS reference system for providing ephemeris and/or clock corrections for GPS satellites to users at measurement epochs of the GPS reference system, the GPS reference system comprising:
a network of reference stations, each of the reference stations receiving GPS signals from those of the GPS satellites it observes and computing in response at each of the measurement epochs a pseudo-range residual for each of the GPS satellites it observes at the measurement epoch; and
a master station in communication with the reference stations to receive the computed pseudo-range residuals at each of the measurement epochs, the master station comprising:
a pseudo-range residual synchronizer to remove at each of the measurement epochs clock differences between the reference stations from the received pseudo-range residuals at the measurement epoch so as to synchronize them in time;
for each of the GPS satellites observed at each of the measurement Epochs by at least one of the reference stations, a corresponding, independent ephemeris and clock correction estimator to compute at the measurement epoch an ephemeris correction and a clock correction for the corresponding GPS satellite in response to each of the synchronized pseudo-range residuals at the measurement for the corresponding GPS satellite and independent of the synchronized pseudo-range residuals at the measurement for the other GPS satellites observed by the reference stations; and
a transmitter to transmit to users at each of the measurement epochs the computed ephemeris and/or clock corrections for each of the GPS satellites observed at the measurement epoch by at least one of the reference stations;
wherein, for at least a subset of the measurement epochs, at least a subset of the GPS satellites for which the ephemeris and/or clock corrections are transmitted are observed by at least one but no more than three of the reference stations such that the transmitted ephemeris and/or clock corrections for each of the at least subset of GPS satellites at each of the at least subset of measurement epochs are computed in response to at least one but no more than three synchronized pseudo-range residuals for the GPS satellite at the measurement epoch.
2. A GPS reference system as recited in claim 1 wherein the network of reference stations comprises no more than three reference stations.
3. A GPS reference system as recited in claim 1 wherein the network of reference stations comprises at least four reference stations.
4. A GPS reference system as recited in claim 1 wherein:
the network of reference stations comprises no more than three reference stations;
each of the ephemeris and clock correction estimators computes at each of the measurement epochs the ephemeris and clock corrections for the corresponding GPS satellite when it is observed at the measurement epoch by at least one of the reference stations in accordance with:
x.sup.k =G.sub.W.sup.T (G.sub.W G.sub.W.sup.T).sup.-1 Δρ.sup.k
where:
k represents the corresponding GPS satellite;
xk is a state vector containing the ephemeris and clock corrections at the measurement epoch for the corresponding GPS satellite;
Δρk is a vector containing, for each of the reference stations that observes the corresponding GPS satellite at the measurement epoch, the corresponding synchronized pseudo-range residual at the measurement epoch for the corresponding GPS satellite; and
GW is a matrix containing, for each of the reference stations that observes the corresponding GPS satellite at the measurement epoch, a line of sight from the reference station to the corresponding GPS satellite and an element that indicates that the corresponding synchronized pseudo-range residual in the vector Δρk should be adjusted by the clock correction in the state vector xk for the corresponding GPS satellite.
5. A GPS reference system as recited in claim 1 wherein:
the network of reference stations comprises at least four reference stations;
each of the ephemeris and clock correction estimators computes at each of the measurement epochs the ephemeris and clock corrections for the corresponding GPS satellite when it is observed at the measurement epoch by at least one but no more than three of the reference stations in accordance with:
x.sup.k =G.sub.W.sup.T (G.sub.W G.sub.W.sup.T).sup.-1 Δρ.sup.k
and when it is observed at the measurement epoch by at least four of the reference stations in accordance with:
x.sup.k =(G.sub.W.sup.T G.sub.W).sup.-1 G.sub.W.sup.T Δρ.sup.k
where:
k represents the corresponding GPS satellite;
xk is a state vector containing the ephemeris and clock corrections at the measurement epoch for the corresponding GPS satellite;
Δρk is a vector containing, for each of the reference stations that observes the corresponding GPS satellite at the measurement epoch, the corresponding synchronized pseudo-range residual at the measurement epoch for the corresponding GPS satellite; and
GW is a matrix containing, for each of the reference stations that observes the corresponding GPS satellite at the measurement epoch, a line of sight from the reference station to the corresponding GPS satellite and an element that indicates that the corresponding synchronized pseudo-range residual in the vector Δρk should be adjusted by the clock correction in the state vector xk for the corresponding GPS satellite.
6. A GPS reference system for providing ephemeris and/or clock corrections for GPS satellites to users at measurement epochs of the GPS reference system, the GPS reference system comprising:
a network of reference stations including a primary reference station and secondary reference stations, each of the reference stations receiving GPS signals from those of the GPS satellites it observes and computing in response at each of the measurement epochs a pseudo-range residual for each of the GPS satellites it observes at the measurement epoch, each of the secondary reference stations positioned with respect to the primary reference station such that at least one of the GPS satellites it observes at each of the measurement epochs is simultaneously observed by the primary reference station; and
a master station in communication with the reference stations to receive the computed pseudo-range residuals at each of the measurement epochs, the master station comprising:
a pseudo-range residual synchronizer remove at each of the measurement epochs the clock differences from the received pseudo-range residuals at the measurement epoch so as to synchronize them in time in response to the received pseudo-range residuals at the measurement for each of the GPS satellites that is simultaneously observed at the measurement epoch by the primary reference station and at least one of the secondary reference stations;
for each of the GPS satellites observed at each of the measurements epoch by the primary reference station and at least one of the secondary reference stations, a corresponding, independent ephemeris and clock correction estimator to compute at the measurement epoch an ephemeris correction and a clock correction for the corresponding GPS satellite in response to each of the synchronized pseudo-range residuals at the measurement epoch for the corresponding GPS satellite and independent of the synchronized pseudo-range residuals at the measurement epoch for the other GPS satellites observed by the reference stations; and
a transmitter to transmit to users at each of the measurement epochs the computed ephemeris and/or clock corrections for each of the GPS satellites observed at the measurement epoch by the primary reference station and at least one of the secondary reference stations;
wherein, for at least a subset of the measurement epochs, at least a subset of the GPS satellites for which the ephemeris and/or clock corrections are transmitted are observed by the primary reference station and at least one but no more than two of the secondary reference stations such that the transmitted ephemeris and/or clock corrections for each of the at least subset of GPS satellites at each of the at least subset of measurement epochs are computed in response to at least two but no more than three synchronized pseudo-range residuals for the GPS satellite at the measurement epoch.
7. A GPS reference system as recited in claim 6 wherein the network of reference stations comprises no more than three reference stations.
8. A GPS reference system as recited in claim 6 wherein the network of reference stations comprises at least four reference stations.
9. A GPS reference system as recited in claim 6 wherein:
the network of reference stations comprises no more than three reference stations;
each of the ephemeris and clock correction estimators computes at each of the measurement epochs the ephemeris correction for the corresponding GPS satellite when it is observed at the measurement epoch by the primary reference station and at least one of the secondary reference stations in accordance with:
Δr.sup.k =G.sub.W.sup.T (G.sub.W G.sub.W.sup.T).sup.-1 Δ.sup.k
where:
k represents the corresponding GPS satellite;
Δrk is a vector containing the ephemeris correction at the measurement epoch for the corresponding GPS satellite;
Δk is a vector containing, for each of the secondary reference stations that observes the corresponding GPS satellite at the measurement epoch, the difference between the synchronized pseudo-range residual at the measurement epoch for the corresponding GPS satellite from the secondary reference station and the synchronized pseudo-range residual at the measurement epoch for the corresponding GPS satellite from the primary reference station;
GW is a matrix containing, for each of the secondary reference stations that observes the corresponding GPS satellite at the measurement epoch, the difference between a line of sight from the secondary reference station to the corresponding GPS satellite and a line sight from the primary reference station to the corresponding GPS satellite; and
each of the ephemeris and clock correction estimators computes the clock correction at each of the measurement epochs for the corresponding GPS satellite in response to the computed ephemeris correction at the measurement epoch for the corresponding GPS satellite.
10. A GPS reference system as recited in claim 6 wherein:
the network of reference stations comprises at least four reference stations;
each of the ephemeris and clock correction estimators computes the ephemeris correction at each of the measurement epochs for the corresponding GPS satellite when it is observed at the measurement epoch by the primary reference station and at least one but no more than two of the secondary reference stations in accordance with:
Δr.sup.k =G.sub.W.sup.T (G.sub.W G.sub.W.sup.T).sup.-1 Δ.sup.k
and when it is observed at the measurement epoch by the primary reference station and at least three of the secondary reference stations in accordance with: where:
Δr.sup.k =(G.sub.W.sup.T G.sub.W).sup.-1 G.sub.W.sup.T Δ.sup.k
k represents the corresponding GPS satellite;
Δrk is a vector containing the ephemeris correction at the measurement epoch for the corresponding GPS satellite;
Δk is a vector containing, for each of the secondary reference stations that observes the corresponding GPS satellite at the measurement epoch, the difference between the synchronized pseudo-range residual at the measurement epoch for the corresponding GPS satellite from the secondary reference station and the synchronized pseudo-range residual at the measurement epoch for the corresponding GPS satellite from the primary reference station; and
GW is a matrix containing, for each of the secondary reference stations that observes the corresponding GPS satellite at the measurement epoch, the difference between a line of sight from the secondary reference station to the corresponding GPS satellite and a line of sight from the primary reference station to the corresponding GPS satellite; and
each of the ephemeris and clock correction estimators computes the clock correction at each of the measurement epochs for the corresponding GPS satellite in response to the computed ephemeris correction at the measurement epoch for the corresponding GPS satellite.
11. A method of providing ephemeris and/or clock corrections for GPS satellites to users at measurement epochs of the method utilizing a network of reference stations and a master station, the method comprising the steps of:
with each of the reference stations:
receiving GPS signals from those of the GPS satellites it observes and in response computing at each of the measurement epochs a pseudo-range residual for each of the GPS satellites it observes at the measurement epoch;
providing the computed pseudo-range residuals at each of the measurement epochs to the master station;
with the master station:
removing at each of the measurement epochs clock differences from the provided pseudo-range residuals so as to synchronize them in response to the provided pseudo-range residuals at the measurement epoch;
for each of the GPS satellites observed at each of the measurement epochs by the reference stations, computing at the measurement epoch an ephemeris correction and a clock correction for the corresponding GPS satellite in response to the synchronized pseudo-range residuals at the measurement epoch for the corresponding GPS satellite and independent of the synchronized pseudo-range residuals at the measurement epoch for the other GPS satellites observed by the reference stations; and
transmitting at each of the measurement epochs the ephemeris and/or clock corrections for each of the GPS satellites observed at the measurement epoch by at least one of the reference stations;
wherein, for at least a subset of the measurement epochs, at least a subset of the GPS satellites for which the ephemeris and/or clock corrections are transmitted are observed by at least one but no more than three of the reference stations such that the transmitted ephemeris and/or clock corrections for each of the at least subset of GPS satellites at each of the at least subset of measurement epochs are computed in response to at least one but no more than three synchronized pseudo-range residuals for the GPS satellite at the measurement epoch.
12. A method as recited in claim 11 wherein the network of reference stations comprises no more than three reference stations.
13. A method as recited in claim 11 wherein the network of reference stations comprises at least four reference stations.
14. A method as recited in claim 11 wherein:
the network of reference stations comprises no more than three reference stations;
the ephemeris and clock corrections at each of the measurement epochs are computed in the ephemeris and clock correction computing step for the corresponding GPS satellite when it is observed at the measurement epoch by at least one of the reference stations in accordance with:
x.sup.k =G.sub.W.sup.T (G.sub.W G.sub.W.sup.T).sup.-1 Δρ.sup.k
where:
k represents the corresponding GPS satellite;
xk is a state vector containing the ephemeris and clock corrections at the measurement epoch for the corresponding GPS satellite;
Δρk is a vector containing, for each of the reference stations that observes the corresponding GPS satellite at the measurement epoch, the corresponding synchronized pseudo-range residual at the measurement epoch for the corresponding GPS satellite; and
GW is a matrix containing, for each of the reference stations that observes the corresponding GPS satellite at the measurement epoch, a line of sight from the reference station to the corresponding GPS satellite and an element that indicates that the corresponding synchronized pseudo-range residual in the vector Δρk should be adjusted by the clock correction in the state vector xk for the corresponding GPS satellite.
15. A method as recited in claim 11 wherein:
the network of reference stations comprises at least four reference stations;
the ephemeris and clock corrections at each of the measurement epochs are computed in the ephemeris and clock correction computing step for the corresponding GPS satellite when it is observed at the measurement epoch by at least one but no more than three of the reference stations in accordance with:
x.sup.k =G.sub.W.sup.T (G.sub.W G.sub.W.sup.T).sup.-1 Δρ.sup.k
and when it is observed at the measurement epoch by at least four of the reference stations in accordance with:
x.sup.k =(G.sub.W.sup.T G.sub.W).sup.-1 G.sub.W.sup.T Δρ.sup.k
where:
k represents the corresponding GPS satellite;
xk is a state vector containing the ephemeris and clock corrections at the measurement epoch for the corresponding GPS satellite;
Δpk is a vector containing, for each of the reference stations that observes the corresponding GPS satellite at the measurement epoch, the corresponding synchronized pseudo-range residual at the measurement epoch for the corresponding GPS satellite; and
GW is a matrix containing, for each of the reference stations that observes the corresponding GPS satellite at the measurement epoch, a line of sight from the reference station to the corresponding GPS satellite and an element that indicates that the corresponding synchronized pseudo-range residual in the vector Δρk should be adjusted by the clock correction in the state vector xk for the corresponding GPS satellite.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US08/681,087 US5621646A (en) | 1995-01-17 | 1996-07-22 | Wide area differential GPS reference system and method |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US37416595A | 1995-01-17 | 1995-01-17 | |
US08/681,087 US5621646A (en) | 1995-01-17 | 1996-07-22 | Wide area differential GPS reference system and method |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US37416595A Continuation | 1995-01-17 | 1995-01-17 |
Publications (1)
Publication Number | Publication Date |
---|---|
US5621646A true US5621646A (en) | 1997-04-15 |
Family
ID=23475593
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US08/681,087 Expired - Fee Related US5621646A (en) | 1995-01-17 | 1996-07-22 | Wide area differential GPS reference system and method |
Country Status (2)
Country | Link |
---|---|
US (1) | US5621646A (en) |
WO (1) | WO1996022546A1 (en) |
Cited By (102)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5760738A (en) * | 1995-09-29 | 1998-06-02 | Kabushiki Kaisha Toshiba | GPS overlay system for a geostationary satellite |
US5812961A (en) * | 1995-12-28 | 1998-09-22 | Trimble Navigation Limited | Method and reciever using a low earth orbiting satellite signal to augment the global positioning system |
US5877725A (en) * | 1997-03-06 | 1999-03-02 | Trimble Navigation Limited | Wide augmentation system retrofit receiver |
US5884220A (en) * | 1996-07-16 | 1999-03-16 | Trimble Navigation Limited | Method and apparatus to improve overall performance of a DGPS receiver |
WO1999018677A1 (en) * | 1997-10-02 | 1999-04-15 | Airsys Atm, Inc. | Gps augmentation system |
US6040798A (en) * | 1995-10-24 | 2000-03-21 | International Mobile Satellite Organization | Satellite radiodetermination |
WO2000020882A1 (en) * | 1998-09-28 | 2000-04-13 | Telefonaktiebolaget Lm Ericsson (Publ) | Improved use of global positioning system in locating a radio transmitter |
US6078286A (en) * | 1998-06-03 | 2000-06-20 | Motorola, Inc. | Method and apparatus for efficient acquisition and tracking of satellites |
US6104338A (en) * | 1998-05-04 | 2000-08-15 | Snaptrack, Inc. | Method and apparatus for operating a satellite positioning system receiver |
US6107939A (en) * | 1998-11-05 | 2000-08-22 | Trimble Navigation Limited | Lane change alarm for use in a highway vehicle |
US6133870A (en) * | 1999-06-29 | 2000-10-17 | Trw Inc. | Ephemeris determination using intersatellite ranging and ephemeris knowledge of other satellites |
WO2000065367A1 (en) * | 1999-04-23 | 2000-11-02 | Global Locate, Inc. | Method and apparatus for locating and providing services to mobile devices |
US6215441B1 (en) | 1997-04-15 | 2001-04-10 | Snaptrack, Inc. | Satellite positioning reference system and method |
US6230097B1 (en) | 1998-08-31 | 2001-05-08 | Trimble Navigation Limited | Accurate vehicle navigation |
US6278402B1 (en) * | 1998-09-28 | 2001-08-21 | Vxt, L.L.C. | Hybrid system and method for golf course distance determination using GPS navigation signals |
US6356232B1 (en) | 1999-12-17 | 2002-03-12 | University Corporation For Atmospheric Research | High resolution ionospheric technique for regional area high-accuracy global positioning system applications |
US6370475B1 (en) | 1997-10-22 | 2002-04-09 | Intelligent Technologies International Inc. | Accident avoidance system |
US6397147B1 (en) | 2000-06-06 | 2002-05-28 | Csi Wireless Inc. | Relative GPS positioning using a single GPS receiver with internally generated differential correction terms |
US6405132B1 (en) | 1997-10-22 | 2002-06-11 | Intelligent Technologies International, Inc. | Accident avoidance system |
US6411254B1 (en) * | 1997-04-15 | 2002-06-25 | Snaptrack, Inc. | Satellite positioning reference system and method |
US6445340B2 (en) | 2000-01-26 | 2002-09-03 | Novatel, Inc. | Method and apparatus for generating a time signal |
US20020172309A1 (en) * | 2001-05-15 | 2002-11-21 | International Business Machines Corporation | Universal clock reference |
US6539293B2 (en) * | 1999-04-01 | 2003-03-25 | Siemens Schweiz Ag | Method and device for monitoring bogies of multi-axle vehicles |
US6552680B1 (en) * | 2001-10-01 | 2003-04-22 | Garmin Ltd. | Method and system for minimizing storage and processing of ionospheric grid point correction information |
US6563765B1 (en) * | 1999-06-16 | 2003-05-13 | Matsushita Electric Industrial Co., Ltd. | Clock system |
US6567041B1 (en) | 2001-04-18 | 2003-05-20 | Sprint Spectrum, L.P. | Network system and method for a remote reference receiver system |
US6618671B2 (en) | 2000-09-20 | 2003-09-09 | Koninklijke Philips Electronics N.V. | Method of determining the position of a mobile unit |
WO2003096059A1 (en) * | 2002-05-13 | 2003-11-20 | Qualcomm Incorporated | System and method for reference data processing in network assisted position determination |
US6720920B2 (en) | 1997-10-22 | 2004-04-13 | Intelligent Technologies International Inc. | Method and arrangement for communicating between vehicles |
US20040140930A1 (en) * | 2001-03-29 | 2004-07-22 | Guy Harles | Ranging system for determining ranging information of a spacecraft |
US6799116B2 (en) | 2000-12-15 | 2004-09-28 | Trimble Navigation Limited | GPS correction methods, apparatus and signals |
US6829534B2 (en) | 1999-04-23 | 2004-12-07 | Global Locate, Inc. | Method and apparatus for performing timing synchronization |
US6839631B1 (en) | 2001-12-21 | 2005-01-04 | Garmin Ltd. | Space based augmentation system with hierarchy for determining geographical corrections source |
US20050024263A1 (en) * | 2003-07-30 | 2005-02-03 | Sharpe Richard T. | Method for generating clock corrections for a wide-area or global differential GPS system |
US20050219119A1 (en) * | 2004-02-26 | 2005-10-06 | Po-Hsin Hsu | Method and apparatus for signal code carrier coherency for wide area augmentation system |
US20060017610A1 (en) * | 2004-07-23 | 2006-01-26 | California Institute Of Technology | Generating high precision ionospheric ground-truth measurements |
US20060033658A1 (en) * | 2004-08-16 | 2006-02-16 | Camp William O Jr | Apparatus, methods and computer program products for GPS signal acquisition using an adaptive search engine |
US20060034354A1 (en) * | 2004-08-16 | 2006-02-16 | Camp William O Jr | Apparatus, methods and computer program products for positioning system signal processing using parallel computational techniques |
US20060064244A1 (en) * | 1994-01-03 | 2006-03-23 | Robbins James E | Differential GPS corrections using virtual stations |
US7031730B1 (en) * | 2001-10-01 | 2006-04-18 | Garmin Ltd. | Method and system for minimizing storage and processing of ionospheric grid point correction information in a wireless communications device |
FR2877437A1 (en) * | 2004-11-04 | 2006-05-05 | Thales Sa | POSITIONING DEVICE |
US7126527B1 (en) * | 2000-06-23 | 2006-10-24 | Intel Corporation | Method and apparatus for mobile device location via a network based local area augmentation system |
US20070096978A1 (en) * | 2005-10-27 | 2007-05-03 | Global Locate, Inc. | Method And Apparatus For Mitigating Multipath Effects At A Satellite Signal Receiver Using a Sequential Estimation Filter |
US20070149209A1 (en) * | 2005-12-22 | 2007-06-28 | Novariant Inc. | Wide area differential satellite positioning with local area updates |
US20070198185A1 (en) * | 2002-12-11 | 2007-08-23 | Mcclure John A | GNSS control system and method |
US20070241957A1 (en) * | 2001-10-05 | 2007-10-18 | Murphy Timothy A | Method and apparatus for providing an integrated communications, navigation and surveillance satellite system |
US20070299609A1 (en) * | 2006-06-23 | 2007-12-27 | Nemerix Sa | Method and system for ephemeris extension for GNSS applications |
US7327304B2 (en) | 2005-03-24 | 2008-02-05 | Agilent Technologies, Inc. | System and method for minimizing background noise in a microwave image using a programmable reflector array |
US20080167770A1 (en) * | 2007-01-05 | 2008-07-10 | Beeline Technologies Pty Ltd | Vehicle control system |
US20080175105A1 (en) * | 2006-07-31 | 2008-07-24 | Seiko Epson Corporation | Time Correction Device, Timepiece Having a Time Correction Device, and Time Correction Method |
US20080186228A1 (en) * | 2007-02-05 | 2008-08-07 | Nemerix Sa | Prediction refresh method for ephemeris extensions |
US20080205494A1 (en) * | 2007-02-27 | 2008-08-28 | Whitehead Michael L | Unbiased code phase discriminator |
US7453956B2 (en) | 2004-08-16 | 2008-11-18 | Sony Ericsson Mobile Communications Ab | Apparatus, methods and computer program products for signal acquisition using common demodulation templates |
US20090043435A1 (en) * | 2007-08-07 | 2009-02-12 | Quantum Engineering, Inc. | Methods and systems for making a gps signal vital |
US20090085815A1 (en) * | 2007-09-27 | 2009-04-02 | Jakab Andrew J | Tightly-coupled pcb gnss circuit and manufacturing method |
US20090099730A1 (en) * | 2003-03-20 | 2009-04-16 | Mcclure John A | Satellite based vehicle guidance control in straight and contour modes |
US20090121932A1 (en) * | 2003-03-20 | 2009-05-14 | Whitehead Michael L | Multi-antenna gnss positioning method and system |
US7535402B1 (en) * | 2004-04-19 | 2009-05-19 | Novariant, Inc. | Navigation with satellite communications |
US20090164067A1 (en) * | 2003-03-20 | 2009-06-25 | Whitehead Michael L | Multiple-antenna gnss control system and method |
US7574215B1 (en) | 2000-11-06 | 2009-08-11 | Trimble Navigation Limited | System and method for distribution of GPS satellite information |
US20090204281A1 (en) * | 2008-02-10 | 2009-08-13 | Hemisphere Gps Llc | Visual, gnss and gyro autosteering control |
US20090219206A1 (en) * | 2006-04-24 | 2009-09-03 | Nemerix Sa | Ephemeris extension method for gnss applications |
US20090251366A1 (en) * | 2008-04-08 | 2009-10-08 | Mcclure John A | Gnss-based mobile communication system and method |
US20100109944A1 (en) * | 2003-03-20 | 2010-05-06 | Whitehead Michael L | Gnss-based tracking of fixed or slow-moving structures |
US20100161179A1 (en) * | 2008-12-22 | 2010-06-24 | Mcclure John A | Integrated dead reckoning and gnss/ins positioning |
US20100176991A1 (en) * | 2008-12-11 | 2010-07-15 | Webber Mark R | Gnss superband asic with simultaneous multi-frequency down conversion |
US20100185366A1 (en) * | 2005-07-19 | 2010-07-22 | Heiniger Richard W | Adaptive machine control system and method |
US20100188285A1 (en) * | 2009-01-23 | 2010-07-29 | Her Majesty The Queen In Right Of Canada As Represented By The Minister Of Natural Resources | Decoupled clock model with ambiguity datum fixing |
US20100245172A1 (en) * | 2007-12-18 | 2010-09-30 | Thales Alenia Space Italia S.P.A | Method of Synchronising Nodes of a Network, and System and Device Therefor |
US20100312428A1 (en) * | 2003-03-20 | 2010-12-09 | Roberge Andre C | Gnss guidance and machine control |
US20100332058A1 (en) * | 2009-06-30 | 2010-12-30 | Quantum Engineering, Inc. | Vital speed profile to control a train moving along a track |
US20110015817A1 (en) * | 2009-07-17 | 2011-01-20 | Reeve David R | Optical tracking vehicle control system and method |
US20110022238A1 (en) * | 2009-07-22 | 2011-01-27 | Pollock Colin J | Gnss control system and method for irrigation and related applications |
US20110018765A1 (en) * | 2007-10-08 | 2011-01-27 | Whitehead Michael L | Gnss receiver and external storage device system and gnss data processing method |
US20110025555A1 (en) * | 2009-07-29 | 2011-02-03 | Whitehead Michael L | System and method for augmenting dgnss with internally-generated differential correction |
US20110054729A1 (en) * | 2004-03-19 | 2011-03-03 | Whitehead Michael L | Multi-antenna gnss control system and method |
US20110057834A1 (en) * | 2009-09-04 | 2011-03-10 | Miller Steven R | Multi-frequency gnss receiver baseband dsp |
US20110074628A1 (en) * | 2009-09-30 | 2011-03-31 | Astrium Limited | Positioning system |
US20110172887A1 (en) * | 2009-11-30 | 2011-07-14 | Reeve David R | Vehicle assembly control method for collaborative behavior |
US8085196B2 (en) | 2009-03-11 | 2011-12-27 | Hemisphere Gps Llc | Removing biases in dual frequency GNSS receivers using SBAS |
US8265826B2 (en) | 2003-03-20 | 2012-09-11 | Hemisphere GPS, LLC | Combined GNSS gyroscope 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 |
US8386129B2 (en) | 2009-01-17 | 2013-02-26 | Hemipshere GPS, LLC | Raster-based contour swathing for guidance and variable-rate chemical application |
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 |
US20130325323A1 (en) | 1998-10-22 | 2013-12-05 | American Vehicular Sciences | Vehicle software upgrade techniques |
US8649930B2 (en) | 2009-09-17 | 2014-02-11 | Agjunction Llc | GNSS integrated multi-sensor control system and method |
US20150022395A1 (en) * | 2013-07-16 | 2015-01-22 | O2Micro Inc. | Method and receiver for determining system time of a navigation system |
RU2542325C1 (en) * | 2013-07-24 | 2015-02-20 | Олег Иванович Завалишин | Aircraft location method |
RU2554090C1 (en) * | 2013-12-24 | 2015-06-27 | Акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва" | Method of protection of spacecraft command-and-measurement system |
US9443358B2 (en) | 1995-06-07 | 2016-09-13 | Automotive Vehicular Sciences LLC | Vehicle software upgrade techniques |
US9880562B2 (en) | 2003-03-20 | 2018-01-30 | Agjunction Llc | GNSS and optical guidance and machine control |
US9989644B2 (en) | 2016-02-09 | 2018-06-05 | Honeywell International Inc. | Use of wide area reference receiver network data to mitigate local area error sources |
RU2666554C1 (en) * | 2017-10-06 | 2018-09-11 | Олег Иванович Завалишин | Method for increasing safety of flight and landing of aircraft by local control and corrective station |
RU183917U1 (en) * | 2018-02-21 | 2018-10-09 | Федеральное государственное автономное образовательное учреждение высшего образования "Санкт-Петербургский государственный университет аэрокосмического приборостроения" | RADIO NAVIGATION SYSTEM SIGNAL FREQUENCY MONITORING DEVICE |
USRE47101E1 (en) | 2003-03-20 | 2018-10-30 | Agjunction Llc | Control for dispensing material from vehicle |
US10345448B2 (en) | 2016-01-21 | 2019-07-09 | Honeywell International Inc. | Using space based augmentation system (SBAS) ephemeris sigma information to reduce ground based augmentation systems (GBAS) ephemeris decorrelation parameter |
RU2697811C2 (en) * | 2018-01-22 | 2019-08-20 | Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации | Simulator of radio navigation signals |
RU2715059C1 (en) * | 2019-10-02 | 2020-02-25 | Акционерное общество научно-внедренческое предприятие "ПРОТЕК" | Method of determining coordinates of an aircraft in a satellite-pseudo-satellites multi-position surveillance system |
USRE48527E1 (en) | 2007-01-05 | 2021-04-20 | Agjunction Llc | Optical tracking vehicle control system and method |
CN115453582A (en) * | 2022-08-17 | 2022-12-09 | 武汉大学 | Clock model-based satellite clock error estimation method and system |
CN116908891A (en) * | 2023-02-03 | 2023-10-20 | 中国科学院国家天文台 | Frequency offset correction method and device applied to navigation signal of ground station |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19634529A1 (en) * | 1996-08-27 | 1998-03-12 | Rohde & Schwarz | DGPS tracking system |
GB9618067D0 (en) * | 1996-08-29 | 1996-10-09 | Philips Electronics Nv | Combination of a GPS receiver and a telecommunications apparatus |
TW526335B (en) * | 1998-11-12 | 2003-04-01 | Nippon Kouatsu Electric Co Ltd | Fault point location system |
JP4609095B2 (en) * | 2005-02-04 | 2011-01-12 | セイコーエプソン株式会社 | Positioning system |
FR2890265A1 (en) * | 2005-08-30 | 2007-03-02 | Alcatel Sa | Information broadcast consistency checking device for e.g. Galileo system, has processing unit generating alarm each time when combination between selected projection and distance representing difference between offsets is above threshold |
US8620306B2 (en) * | 2008-06-24 | 2013-12-31 | Qualcomm Incorporated | Location specific search for regional satellite vehicles |
CN102004258B (en) * | 2010-09-09 | 2012-10-10 | 中国计量科学研究院 | Time frequency transfer method and receiver based on multiple global navigation satellite system (GNSS) integration |
CN106533528B (en) * | 2016-09-27 | 2019-04-05 | 北京无线电计量测试研究所 | The two-way time compares compression and the calculation method of modem transmitted data amount |
CN106908817B (en) * | 2017-02-15 | 2019-06-25 | 深圳思凯微电子有限公司 | Assisting navigation localization method and system |
CN108931915B (en) * | 2018-05-08 | 2020-11-06 | 和芯星通科技(北京)有限公司 | Time service method and device using navigation satellite and computer readable storage medium |
CN117741716A (en) * | 2024-02-19 | 2024-03-22 | 中国电力科学研究院有限公司 | Carrier phase and pseudo-range jump exception handling method, module, equipment and medium |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4646096A (en) * | 1984-10-05 | 1987-02-24 | Litton Systems, Inc. | Enhanced global positioning system Delta-Range processing |
US4797677A (en) * | 1982-10-29 | 1989-01-10 | Istac, Incorporated | Method and apparatus for deriving pseudo range from earth-orbiting satellites |
US4912475A (en) * | 1987-03-20 | 1990-03-27 | Massachusetts Institute Of Technology | Techniques for determining orbital data |
US5093800A (en) * | 1989-09-14 | 1992-03-03 | Saburo Ifune | Global positioning system satellite signal simulator |
US5155490A (en) * | 1990-10-15 | 1992-10-13 | Gps Technology Corp. | Geodetic surveying system using multiple GPS base stations |
US5202829A (en) * | 1991-06-10 | 1993-04-13 | Trimble Navigation Limited | Exploration system and method for high-accuracy and high-confidence level relative position and velocity determinations |
US5225842A (en) * | 1991-05-09 | 1993-07-06 | Navsys Corporation | Vehicle tracking system employing global positioning system (gps) satellites |
US5323322A (en) * | 1992-03-05 | 1994-06-21 | Trimble Navigation Limited | Networked differential GPS system |
US5379224A (en) * | 1991-11-29 | 1995-01-03 | Navsys Corporation | GPS tracking system |
US5390124A (en) * | 1992-12-01 | 1995-02-14 | Caterpillar Inc. | Method and apparatus for improving the accuracy of position estimates in a satellite based navigation system |
-
1996
- 1996-01-17 WO PCT/US1996/000672 patent/WO1996022546A1/en active Application Filing
- 1996-07-22 US US08/681,087 patent/US5621646A/en not_active Expired - Fee Related
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4797677A (en) * | 1982-10-29 | 1989-01-10 | Istac, Incorporated | Method and apparatus for deriving pseudo range from earth-orbiting satellites |
US4646096A (en) * | 1984-10-05 | 1987-02-24 | Litton Systems, Inc. | Enhanced global positioning system Delta-Range processing |
US4912475A (en) * | 1987-03-20 | 1990-03-27 | Massachusetts Institute Of Technology | Techniques for determining orbital data |
US5093800A (en) * | 1989-09-14 | 1992-03-03 | Saburo Ifune | Global positioning system satellite signal simulator |
US5155490A (en) * | 1990-10-15 | 1992-10-13 | Gps Technology Corp. | Geodetic surveying system using multiple GPS base stations |
US5225842A (en) * | 1991-05-09 | 1993-07-06 | Navsys Corporation | Vehicle tracking system employing global positioning system (gps) satellites |
US5202829A (en) * | 1991-06-10 | 1993-04-13 | Trimble Navigation Limited | Exploration system and method for high-accuracy and high-confidence level relative position and velocity determinations |
US5379224A (en) * | 1991-11-29 | 1995-01-03 | Navsys Corporation | GPS tracking system |
US5323322A (en) * | 1992-03-05 | 1994-06-21 | Trimble Navigation Limited | Networked differential GPS system |
US5390124A (en) * | 1992-12-01 | 1995-02-14 | Caterpillar Inc. | Method and apparatus for improving the accuracy of position estimates in a satellite based navigation system |
Non-Patent Citations (16)
Title |
---|
C. Goad, "Surveying with GPS", The Global Positioning System, pp. 1-17. |
C. Goad, Surveying with GPS , The Global Positioning System, pp. 1 17. * |
C. Kee, "Wide Area Differential GPS(WADGPS)", Changdon Kee, Dec. 1993, pp. i-128. |
C. Kee, B. Parkinson, "Algorithms and Implementation of Wide Area Differential GPS", ION-GPS-92, pp. 1-8. |
C. Kee, B. Parkinson, Algorithms and Implementation of Wide Area Differential GPS , ION GPS 92, pp. 1 8. * |
C. Kee, Wide Area Differential GPS(WADGPS) , Changdon Kee, Dec. 1993, pp. i 128. * |
D. Lapucha, M. Huff, "Multi-Site Real-Time DGPS System Using Starfix Link; Operational Results", ION-GSP-92, pp. 581-588. |
D. Lapucha, M. Huff, Multi Site Real Time DGPS System Using Starfix Link; Operational Results , ION GSP 92, pp. 581 588. * |
F. Haas, M. Lage, S. Kalinowski, "GPS Wide Area Augmentation System (WAAS) Testbed Results -Phase 1D Testbed Results", ION Annual Meeting, pp. 239-248. |
F. Haas, M. Lage, S. Kalinowski, GPS Wide Area Augmentation System (WAAS) Testbed Results Phase 1D Testbed Results , ION Annual Meeting, pp. 239 248. * |
H. Janes, R. Langley, "Analysis of Tropospheric Delay Prediction Models: Comparisons with Ray-Tracing and Implications for GPS Relative Positioning, Bulletin Geodesiue", vol. 65, pp. 151-161, 1991. |
H. Janes, R. Langley, Analysis of Tropospheric Delay Prediction Models: Comparisons with Ray Tracing and Implications for GPS Relative Positioning, Bulletin Geodesiue , vol. 65, pp. 151 161, 1991. * |
M. Lage, B. Elrod, "The FAA's WIB/WDGPS Testbed and Recent Results", ION-GPS-93, pp. 487-493. |
M. Lage, B. Elrod, The FAA s WIB/WDGPS Testbed and Recent Results , ION GPS 93, pp. 487 493. * |
T. Walter, C. Kee, Y. C. Chao, Y. J. Tsai, U. Peled, J. Ceva, A. K. Barrows, E. Abbott, D. Powell, P. Enge, and B. Parkinson, "Flight Trials of the Wide Area Augmentation System (WAAS)", ION-GPS-94, pp. 1-9. |
T. Walter, C. Kee, Y. C. Chao, Y. J. Tsai, U. Peled, J. Ceva, A. K. Barrows, E. Abbott, D. Powell, P. Enge, and B. Parkinson, Flight Trials of the Wide Area Augmentation System (WAAS) , ION GPS 94, pp. 1 9. * |
Cited By (175)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7711480B2 (en) | 1994-01-03 | 2010-05-04 | Trimble Navigation Limited | Differential GPS corrections using virtual stations |
US20060064244A1 (en) * | 1994-01-03 | 2006-03-23 | Robbins James E | Differential GPS corrections using virtual stations |
US20060282216A1 (en) * | 1994-01-03 | 2006-12-14 | Robbins James E | Differential GPS corrections using virtual stations |
US9443358B2 (en) | 1995-06-07 | 2016-09-13 | Automotive Vehicular Sciences LLC | Vehicle software upgrade techniques |
US5760738A (en) * | 1995-09-29 | 1998-06-02 | Kabushiki Kaisha Toshiba | GPS overlay system for a geostationary satellite |
US6040798A (en) * | 1995-10-24 | 2000-03-21 | International Mobile Satellite Organization | Satellite radiodetermination |
US20040145517A1 (en) * | 1995-10-24 | 2004-07-29 | Inmarsat, Ltd. | Method and apparatus for a satellite payload and radiodetermination |
US5812961A (en) * | 1995-12-28 | 1998-09-22 | Trimble Navigation Limited | Method and reciever using a low earth orbiting satellite signal to augment the global positioning system |
US5884220A (en) * | 1996-07-16 | 1999-03-16 | Trimble Navigation Limited | Method and apparatus to improve overall performance of a DGPS receiver |
US5877725A (en) * | 1997-03-06 | 1999-03-02 | Trimble Navigation Limited | Wide augmentation system retrofit receiver |
US6215441B1 (en) | 1997-04-15 | 2001-04-10 | Snaptrack, Inc. | Satellite positioning reference system and method |
US6411254B1 (en) * | 1997-04-15 | 2002-06-25 | Snaptrack, Inc. | Satellite positioning reference system and method |
WO1999018677A1 (en) * | 1997-10-02 | 1999-04-15 | Airsys Atm, Inc. | Gps augmentation system |
US6370475B1 (en) | 1997-10-22 | 2002-04-09 | Intelligent Technologies International Inc. | Accident avoidance system |
US6720920B2 (en) | 1997-10-22 | 2004-04-13 | Intelligent Technologies International Inc. | Method and arrangement for communicating between vehicles |
US6405132B1 (en) | 1997-10-22 | 2002-06-11 | Intelligent Technologies International, Inc. | Accident avoidance system |
US6104338A (en) * | 1998-05-04 | 2000-08-15 | Snaptrack, Inc. | Method and apparatus for operating a satellite positioning system receiver |
US6078286A (en) * | 1998-06-03 | 2000-06-20 | Motorola, Inc. | Method and apparatus for efficient acquisition and tracking of satellites |
US6230097B1 (en) | 1998-08-31 | 2001-05-08 | Trimble Navigation Limited | Accurate vehicle navigation |
US6278402B1 (en) * | 1998-09-28 | 2001-08-21 | Vxt, L.L.C. | Hybrid system and method for golf course distance determination using GPS navigation signals |
US6266012B1 (en) | 1998-09-28 | 2001-07-24 | Telefonaktiebolaget Lm Ericsson | Use of global positioning system in locating a radio transmitter |
WO2000020882A1 (en) * | 1998-09-28 | 2000-04-13 | Telefonaktiebolaget Lm Ericsson (Publ) | Improved use of global positioning system in locating a radio transmitter |
EP1818680A1 (en) * | 1998-09-28 | 2007-08-15 | Telefonaktiebolaget LM Ericsson (publ) | Improved use of global positioning system in locating a radio transmitter |
US10240935B2 (en) | 1998-10-22 | 2019-03-26 | American Vehicular Sciences Llc | Vehicle software upgrade techniques |
US20130325323A1 (en) | 1998-10-22 | 2013-12-05 | American Vehicular Sciences | Vehicle software upgrade techniques |
US6107939A (en) * | 1998-11-05 | 2000-08-22 | Trimble Navigation Limited | Lane change alarm for use in a highway vehicle |
US6539293B2 (en) * | 1999-04-01 | 2003-03-25 | Siemens Schweiz Ag | Method and device for monitoring bogies of multi-axle vehicles |
US6487499B1 (en) * | 1999-04-23 | 2002-11-26 | Global Locate, Inc. | Method for adjusting a pseudo-range model |
US20070200753A1 (en) * | 1999-04-23 | 2007-08-30 | Global Locate, Inc. | Method and apparatus for locating position of a gps device |
WO2000065367A1 (en) * | 1999-04-23 | 2000-11-02 | Global Locate, Inc. | Method and apparatus for locating and providing services to mobile devices |
US6510387B2 (en) | 1999-04-23 | 2003-01-21 | Global Locate, Inc. | Correction of a pseudo-range model from a GPS almanac |
US6853916B2 (en) | 1999-04-23 | 2005-02-08 | Global Locate, Inc. | Method and apparatus for forming a pseudo-range model |
US7400974B2 (en) | 1999-04-23 | 2008-07-15 | Global Locate, Inc. | Method and apparatus for locating position of a GPS device |
US7158883B2 (en) | 1999-04-23 | 2007-01-02 | Global Locate, Inc | Method and apparatus for locating position of a GPS device |
US6829534B2 (en) | 1999-04-23 | 2004-12-07 | Global Locate, Inc. | Method and apparatus for performing timing synchronization |
US20070200754A1 (en) * | 1999-04-23 | 2007-08-30 | Global Locate, Inc. | Method and apparatus for locating position of a gps device |
US6484097B2 (en) | 1999-04-23 | 2002-11-19 | Global Locate, Inc. | Wide area inverse differential GPS |
US6795771B2 (en) | 1999-04-23 | 2004-09-21 | Global Locate, Inc. | Method and apparatus for correlating a satellite signal with a reference code |
US6453237B1 (en) | 1999-04-23 | 2002-09-17 | Global Locate, Inc. | Method and apparatus for locating and providing services to mobile devices |
US20050125152A1 (en) * | 1999-04-23 | 2005-06-09 | Global Locate, Inc. | Method and apparatus for locating position of a GPS device |
US6563765B1 (en) * | 1999-06-16 | 2003-05-13 | Matsushita Electric Industrial Co., Ltd. | Clock system |
US6133870A (en) * | 1999-06-29 | 2000-10-17 | Trw Inc. | Ephemeris determination using intersatellite ranging and ephemeris knowledge of other satellites |
US6356232B1 (en) | 1999-12-17 | 2002-03-12 | University Corporation For Atmospheric Research | High resolution ionospheric technique for regional area high-accuracy global positioning system applications |
US6445340B2 (en) | 2000-01-26 | 2002-09-03 | Novatel, Inc. | Method and apparatus for generating a time signal |
US6397147B1 (en) | 2000-06-06 | 2002-05-28 | Csi Wireless Inc. | Relative GPS positioning using a single GPS receiver with internally generated differential correction terms |
US7126527B1 (en) * | 2000-06-23 | 2006-10-24 | Intel Corporation | Method and apparatus for mobile device location via a network based local area augmentation system |
US7089112B2 (en) | 2000-09-20 | 2006-08-08 | Koninklijke Phillips Electronics N.V. | Method of determining the position of a mobile unit |
US6618671B2 (en) | 2000-09-20 | 2003-09-09 | Koninklijke Philips Electronics N.V. | Method of determining the position of a mobile unit |
US7574215B1 (en) | 2000-11-06 | 2009-08-11 | Trimble Navigation Limited | System and method for distribution of GPS satellite information |
US6862526B2 (en) | 2000-12-15 | 2005-03-01 | Trimble Navigation Limited | GPS correction methods, apparatus and signals |
US6799116B2 (en) | 2000-12-15 | 2004-09-28 | Trimble Navigation Limited | GPS correction methods, apparatus and signals |
US20040204852A1 (en) * | 2000-12-15 | 2004-10-14 | Robbins James E. | GPS correction methods, apparatus and signals |
US7512505B2 (en) | 2001-03-29 | 2009-03-31 | Ses Astra S.A. | Ranging system for determining ranging information of a spacecraft |
US20040140930A1 (en) * | 2001-03-29 | 2004-07-22 | Guy Harles | Ranging system for determining ranging information of a spacecraft |
US6567041B1 (en) | 2001-04-18 | 2003-05-20 | Sprint Spectrum, L.P. | Network system and method for a remote reference receiver system |
US20020172309A1 (en) * | 2001-05-15 | 2002-11-21 | International Business Machines Corporation | Universal clock reference |
US7031730B1 (en) * | 2001-10-01 | 2006-04-18 | Garmin Ltd. | Method and system for minimizing storage and processing of ionospheric grid point correction information in a wireless communications device |
US6646594B1 (en) | 2001-10-01 | 2003-11-11 | Garmin Ltd. | Method and system for minimizing storage and processing of ionospheric grid point correction information |
US6552680B1 (en) * | 2001-10-01 | 2003-04-22 | Garmin Ltd. | Method and system for minimizing storage and processing of ionospheric grid point correction information |
US7414573B2 (en) * | 2001-10-05 | 2008-08-19 | The Boeing Company | Method and apparatus for providing an integrated communications, navigation and surveillance satellite system |
US20070241957A1 (en) * | 2001-10-05 | 2007-10-18 | Murphy Timothy A | Method and apparatus for providing an integrated communications, navigation and surveillance satellite system |
US7110883B1 (en) | 2001-12-21 | 2006-09-19 | Garmin Ltd. | Space based augmentation system with hierarchy for determining geographical corrections source |
US6839631B1 (en) | 2001-12-21 | 2005-01-04 | Garmin Ltd. | Space based augmentation system with hierarchy for determining geographical corrections source |
WO2003096059A1 (en) * | 2002-05-13 | 2003-11-20 | Qualcomm Incorporated | System and method for reference data processing in network assisted position determination |
CN100465660C (en) * | 2002-05-13 | 2009-03-04 | 高通股份有限公司 | System and method for reference data processing in network assisted position determination |
US7522588B2 (en) | 2002-05-13 | 2009-04-21 | Qualcomm Incorporated | System and method for reference data processing in network assisted position determination |
US7885745B2 (en) | 2002-12-11 | 2011-02-08 | Hemisphere Gps Llc | GNSS control system and method |
US20070198185A1 (en) * | 2002-12-11 | 2007-08-23 | Mcclure John A | GNSS control system and method |
US20090099730A1 (en) * | 2003-03-20 | 2009-04-16 | Mcclure John A | Satellite based vehicle guidance control in straight and contour modes |
US8594879B2 (en) | 2003-03-20 | 2013-11-26 | Agjunction Llc | GNSS guidance and machine control |
US20100312428A1 (en) * | 2003-03-20 | 2010-12-09 | Roberge Andre C | Gnss guidance and machine control |
US8138970B2 (en) | 2003-03-20 | 2012-03-20 | Hemisphere Gps Llc | GNSS-based tracking of fixed or slow-moving structures |
US8140223B2 (en) | 2003-03-20 | 2012-03-20 | Hemisphere Gps Llc | Multiple-antenna GNSS control system and method |
US10168714B2 (en) | 2003-03-20 | 2019-01-01 | Agjunction Llc | GNSS and optical guidance and machine control |
US8190337B2 (en) | 2003-03-20 | 2012-05-29 | Hemisphere GPS, LLC | Satellite based vehicle guidance control in straight and contour modes |
USRE47101E1 (en) | 2003-03-20 | 2018-10-30 | Agjunction Llc | Control for dispensing material from vehicle |
US9886038B2 (en) | 2003-03-20 | 2018-02-06 | Agjunction Llc | GNSS and optical guidance and machine control |
US9880562B2 (en) | 2003-03-20 | 2018-01-30 | Agjunction Llc | GNSS and optical guidance and machine control |
US8265826B2 (en) | 2003-03-20 | 2012-09-11 | Hemisphere GPS, LLC | Combined GNSS gyroscope control system and method |
US20100109944A1 (en) * | 2003-03-20 | 2010-05-06 | Whitehead Michael L | Gnss-based tracking of fixed or slow-moving structures |
US8686900B2 (en) | 2003-03-20 | 2014-04-01 | Hemisphere GNSS, Inc. | Multi-antenna GNSS positioning method and system |
US20090164067A1 (en) * | 2003-03-20 | 2009-06-25 | Whitehead Michael L | Multiple-antenna gnss control system and method |
US20090121932A1 (en) * | 2003-03-20 | 2009-05-14 | Whitehead Michael L | Multi-antenna gnss positioning method and system |
WO2005043186A2 (en) * | 2003-07-30 | 2005-05-12 | Navcom Technology, Inc. | Method for generating clock corrections for a wide-area or global differential gps system |
US7117417B2 (en) | 2003-07-30 | 2006-10-03 | Navcom Technology, Inc. | Method for generating clock corrections for a wide-area or global differential GPS system |
US20050024263A1 (en) * | 2003-07-30 | 2005-02-03 | Sharpe Richard T. | Method for generating clock corrections for a wide-area or global differential GPS system |
WO2005043186A3 (en) * | 2003-07-30 | 2005-10-20 | Navcom Tech Inc | Method for generating clock corrections for a wide-area or global differential gps system |
US20050219119A1 (en) * | 2004-02-26 | 2005-10-06 | Po-Hsin Hsu | Method and apparatus for signal code carrier coherency for wide area augmentation system |
US8271194B2 (en) | 2004-03-19 | 2012-09-18 | Hemisphere Gps Llc | Method and system using GNSS phase measurements for relative positioning |
US8583315B2 (en) | 2004-03-19 | 2013-11-12 | Agjunction Llc | Multi-antenna GNSS control system and method |
US20110054729A1 (en) * | 2004-03-19 | 2011-03-03 | Whitehead Michael L | Multi-antenna gnss control system and method |
US7535402B1 (en) * | 2004-04-19 | 2009-05-19 | Novariant, Inc. | Navigation with satellite communications |
US7586438B1 (en) * | 2004-04-19 | 2009-09-08 | Novariant Inc. | Navigation with satellite communications |
US20060017610A1 (en) * | 2004-07-23 | 2006-01-26 | California Institute Of Technology | Generating high precision ionospheric ground-truth measurements |
US7289061B2 (en) * | 2004-07-23 | 2007-10-30 | California Institute Of Technology | Generating high precision ionospheric ground-truth measurements |
US7453956B2 (en) | 2004-08-16 | 2008-11-18 | Sony Ericsson Mobile Communications Ab | Apparatus, methods and computer program products for signal acquisition using common demodulation templates |
US7358897B2 (en) | 2004-08-16 | 2008-04-15 | Sony Ericsson Mobile Communicatios Ab | Apparatus, methods and computer program products for GPS signal acquisition using an adaptive search engine |
US20060033658A1 (en) * | 2004-08-16 | 2006-02-16 | Camp William O Jr | Apparatus, methods and computer program products for GPS signal acquisition using an adaptive search engine |
US20060034354A1 (en) * | 2004-08-16 | 2006-02-16 | Camp William O Jr | Apparatus, methods and computer program products for positioning system signal processing using parallel computational techniques |
FR2877437A1 (en) * | 2004-11-04 | 2006-05-05 | Thales Sa | POSITIONING DEVICE |
US7327304B2 (en) | 2005-03-24 | 2008-02-05 | Agilent Technologies, Inc. | System and method for minimizing background noise in a microwave image using a programmable reflector array |
US8214111B2 (en) | 2005-07-19 | 2012-07-03 | Hemisphere Gps Llc | Adaptive machine control system and method |
US20100185366A1 (en) * | 2005-07-19 | 2010-07-22 | Heiniger Richard W | Adaptive machine control system and method |
US8395543B2 (en) * | 2005-10-27 | 2013-03-12 | Global Locate, Inc. | Method and apparatus for mitigating multipath effects at a satellite signal receiver using a sequential estimation filter |
US7253768B2 (en) * | 2005-10-27 | 2007-08-07 | Global Locate, Inc | Method and apparatus for mitigating multipath effects at a satellite signal receiver using a sequential estimation filter |
US20100019960A1 (en) * | 2005-10-27 | 2010-01-28 | Global Locate, Inc. | Method and Apparatus for Mitigating Multipath Effects at a Satellite Signal Receiver Using a Sequential Estimation Filter |
US20070096978A1 (en) * | 2005-10-27 | 2007-05-03 | Global Locate, Inc. | Method And Apparatus For Mitigating Multipath Effects At A Satellite Signal Receiver Using a Sequential Estimation Filter |
US7653399B2 (en) | 2005-12-22 | 2010-01-26 | Novariant, Inc. | Wide area differential satellite positioning with local area updates |
US20070149209A1 (en) * | 2005-12-22 | 2007-06-28 | Novariant Inc. | Wide area differential satellite positioning with local area updates |
US9000979B2 (en) | 2006-04-24 | 2015-04-07 | Qualcomm Incorporated | Ephemeris extension method for GNSS applications |
US20090219206A1 (en) * | 2006-04-24 | 2009-09-03 | Nemerix Sa | Ephemeris extension method for gnss applications |
US8274430B2 (en) | 2006-04-24 | 2012-09-25 | Qualcomm Incorporated | Ephemeris extension method for GNSS applications |
US20070299609A1 (en) * | 2006-06-23 | 2007-12-27 | Nemerix Sa | Method and system for ephemeris extension for GNSS applications |
US20080175105A1 (en) * | 2006-07-31 | 2008-07-24 | Seiko Epson Corporation | Time Correction Device, Timepiece Having a Time Correction Device, and Time Correction Method |
US7474594B2 (en) * | 2006-07-31 | 2009-01-06 | Seiko Epson Corporation | Time correction device, timepiece having a time correction device, and time correction method |
USRE48527E1 (en) | 2007-01-05 | 2021-04-20 | Agjunction Llc | Optical tracking vehicle control system and method |
US20080167770A1 (en) * | 2007-01-05 | 2008-07-10 | Beeline Technologies Pty Ltd | Vehicle control system |
US20110118938A1 (en) * | 2007-01-05 | 2011-05-19 | Andrew John Macdonald | Vehicle control system |
US7835832B2 (en) | 2007-01-05 | 2010-11-16 | Hemisphere Gps Llc | Vehicle control system |
US20080186228A1 (en) * | 2007-02-05 | 2008-08-07 | Nemerix Sa | Prediction refresh method for ephemeris extensions |
US8497801B2 (en) * | 2007-02-05 | 2013-07-30 | Qualcomm Incorporated | Prediction refresh method for ephemeris extensions |
US8654009B2 (en) | 2007-02-05 | 2014-02-18 | Qualcomm Incorporated | Prediction refresh method for ephemeris extensions |
US20080205494A1 (en) * | 2007-02-27 | 2008-08-28 | Whitehead Michael L | Unbiased code phase discriminator |
US8000381B2 (en) | 2007-02-27 | 2011-08-16 | Hemisphere Gps Llc | Unbiased code phase discriminator |
US20090043435A1 (en) * | 2007-08-07 | 2009-02-12 | Quantum Engineering, Inc. | Methods and systems for making a gps signal vital |
US20090085815A1 (en) * | 2007-09-27 | 2009-04-02 | Jakab Andrew J | Tightly-coupled pcb gnss circuit and manufacturing method |
US7948769B2 (en) | 2007-09-27 | 2011-05-24 | Hemisphere Gps Llc | Tightly-coupled PCB GNSS circuit and manufacturing method |
US8456356B2 (en) | 2007-10-08 | 2013-06-04 | Hemisphere Gnss Inc. | GNSS receiver and external storage device system and GNSS data processing method |
US20110018765A1 (en) * | 2007-10-08 | 2011-01-27 | Whitehead Michael L | Gnss receiver and external storage device system and gnss data processing method |
US20100245172A1 (en) * | 2007-12-18 | 2010-09-30 | Thales Alenia Space Italia S.P.A | Method of Synchronising Nodes of a Network, and System and Device Therefor |
US8432850B2 (en) * | 2007-12-18 | 2013-04-30 | Thales Alenia Space Italia, S.p.A. | Method of synchronising nodes of a network, and system and device therefor |
US9002566B2 (en) | 2008-02-10 | 2015-04-07 | AgJunction, LLC | Visual, GNSS and gyro autosteering control |
US20090204281A1 (en) * | 2008-02-10 | 2009-08-13 | Hemisphere Gps Llc | Visual, gnss and gyro autosteering control |
US20090251366A1 (en) * | 2008-04-08 | 2009-10-08 | Mcclure John A | Gnss-based mobile communication system and method |
US8018376B2 (en) | 2008-04-08 | 2011-09-13 | 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 |
US20100176991A1 (en) * | 2008-12-11 | 2010-07-15 | Webber Mark R | Gnss superband asic with simultaneous multi-frequency down conversion |
US20100161179A1 (en) * | 2008-12-22 | 2010-06-24 | Mcclure John A | Integrated dead reckoning and gnss/ins positioning |
US8386129B2 (en) | 2009-01-17 | 2013-02-26 | Hemipshere GPS, LLC | Raster-based contour swathing for guidance and variable-rate chemical application |
USRE47055E1 (en) | 2009-01-17 | 2018-09-25 | Agjunction Llc | Raster-based contour swathing for guidance and variable-rate chemical application |
USRE48509E1 (en) | 2009-01-17 | 2021-04-13 | Agjunction Llc | Raster-based contour swathing for guidance and variable-rate chemical application |
US20100188285A1 (en) * | 2009-01-23 | 2010-07-29 | Her Majesty The Queen In Right Of Canada As Represented By The Minister Of Natural Resources | Decoupled clock model with ambiguity datum fixing |
US8018377B2 (en) | 2009-01-23 | 2011-09-13 | Her Majesty The Queen In Right Of Canada As Represented By The Minister Of Natural Resources | Decoupled clock model with ambiguity datum fixing |
US8085196B2 (en) | 2009-03-11 | 2011-12-27 | Hemisphere Gps Llc | Removing biases in dual frequency GNSS receivers using SBAS |
US8509970B2 (en) | 2009-06-30 | 2013-08-13 | Invensys Rail Corporation | Vital speed profile to control a train moving along a track |
US9168935B2 (en) | 2009-06-30 | 2015-10-27 | Siemens Industry, Inc. | Vital speed profile to control a train moving along a track |
US20100332058A1 (en) * | 2009-06-30 | 2010-12-30 | Quantum Engineering, Inc. | Vital speed profile to control a train moving along a track |
US20110015817A1 (en) * | 2009-07-17 | 2011-01-20 | Reeve David R | Optical tracking vehicle control system and method |
US8311696B2 (en) | 2009-07-17 | 2012-11-13 | Hemisphere Gps Llc | Optical tracking vehicle control system and method |
US8401704B2 (en) | 2009-07-22 | 2013-03-19 | Hemisphere GPS, LLC | GNSS control system and method for irrigation and related applications |
US20110022238A1 (en) * | 2009-07-22 | 2011-01-27 | Pollock Colin J | Gnss control system and method for irrigation and related applications |
US20110025555A1 (en) * | 2009-07-29 | 2011-02-03 | Whitehead Michael L | System and method for augmenting dgnss with internally-generated differential correction |
US8174437B2 (en) | 2009-07-29 | 2012-05-08 | Hemisphere Gps Llc | System and method for augmenting DGNSS with internally-generated differential correction |
US20110057834A1 (en) * | 2009-09-04 | 2011-03-10 | Miller Steven R | Multi-frequency gnss receiver baseband dsp |
US8334804B2 (en) | 2009-09-04 | 2012-12-18 | Hemisphere Gps Llc | Multi-frequency GNSS receiver baseband DSP |
USRE47648E1 (en) | 2009-09-17 | 2019-10-15 | Agjunction Llc | Integrated multi-sensor control system and method |
US8649930B2 (en) | 2009-09-17 | 2014-02-11 | Agjunction Llc | GNSS integrated multi-sensor control system and method |
US8339315B2 (en) | 2009-09-30 | 2012-12-25 | Astrium Limited | Positioning system |
EP2330433A1 (en) * | 2009-09-30 | 2011-06-08 | Astrium Limited | Positioning system |
WO2011039042A1 (en) * | 2009-09-30 | 2011-04-07 | Astrium Limited | Positioning system |
US20110074628A1 (en) * | 2009-09-30 | 2011-03-31 | Astrium Limited | Positioning system |
US8548649B2 (en) | 2009-10-19 | 2013-10-01 | Agjunction Llc | GNSS optimized aircraft control system and method |
US20110172887A1 (en) * | 2009-11-30 | 2011-07-14 | Reeve David R | Vehicle assembly control method for collaborative behavior |
US8583326B2 (en) | 2010-02-09 | 2013-11-12 | Agjunction Llc | GNSS contour guidance path selection |
US20150022395A1 (en) * | 2013-07-16 | 2015-01-22 | O2Micro Inc. | Method and receiver for determining system time of a navigation system |
RU2542325C1 (en) * | 2013-07-24 | 2015-02-20 | Олег Иванович Завалишин | Aircraft location method |
RU2554090C1 (en) * | 2013-12-24 | 2015-06-27 | Акционерное общество "Информационные спутниковые системы" имени академика М.Ф. Решетнёва" | Method of protection of spacecraft command-and-measurement system |
US10345448B2 (en) | 2016-01-21 | 2019-07-09 | Honeywell International Inc. | Using space based augmentation system (SBAS) ephemeris sigma information to reduce ground based augmentation systems (GBAS) ephemeris decorrelation parameter |
US9989644B2 (en) | 2016-02-09 | 2018-06-05 | Honeywell International Inc. | Use of wide area reference receiver network data to mitigate local area error sources |
RU2666554C1 (en) * | 2017-10-06 | 2018-09-11 | Олег Иванович Завалишин | Method for increasing safety of flight and landing of aircraft by local control and corrective station |
RU2697811C2 (en) * | 2018-01-22 | 2019-08-20 | Российская Федерация, от имени которой выступает Министерство обороны Российской Федерации | Simulator of radio navigation signals |
RU183917U1 (en) * | 2018-02-21 | 2018-10-09 | Федеральное государственное автономное образовательное учреждение высшего образования "Санкт-Петербургский государственный университет аэрокосмического приборостроения" | RADIO NAVIGATION SYSTEM SIGNAL FREQUENCY MONITORING DEVICE |
RU2715059C1 (en) * | 2019-10-02 | 2020-02-25 | Акционерное общество научно-внедренческое предприятие "ПРОТЕК" | Method of determining coordinates of an aircraft in a satellite-pseudo-satellites multi-position surveillance system |
CN115453582A (en) * | 2022-08-17 | 2022-12-09 | 武汉大学 | Clock model-based satellite clock error estimation method and system |
CN116908891A (en) * | 2023-02-03 | 2023-10-20 | 中国科学院国家天文台 | Frequency offset correction method and device applied to navigation signal of ground station |
Also Published As
Publication number | Publication date |
---|---|
WO1996022546A1 (en) | 1996-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5621646A (en) | Wide area differential GPS reference system and method | |
US10670734B2 (en) | Advanced navigation satellite system positioning method and system using delayed precise information | |
AU2012304928B2 (en) | Advanced timing and time transfer for satellite constellations using crosslink ranging and an accurate time source | |
US5757316A (en) | Attitude determination utilizing an inertial measurement unit and a plurality of satellite transmitters | |
US5467282A (en) | GPS and satellite navigation system | |
EP0904551B1 (en) | Spoofing detection system for a satellite positioning system | |
US5450448A (en) | Removal of signal errors for differential satellite positioning systems | |
US6317688B1 (en) | Method and apparatus for achieving sole means navigation from global navigation satelite systems | |
CN102540227B (en) | The method and system of aerogram target geo-location in search and rescue system | |
EP1248117A2 (en) | Method, apparatus and system for estimating user position with a satellite positioning system in poor signal conditions | |
US5563917A (en) | Compensation for multipath errors and ionospheric delays in differential satellite positioning systems | |
US8082099B2 (en) | Aircraft navigation using the global positioning system and an attitude and heading reference system | |
US5736960A (en) | Atomic clock augmented global positioning system receivers and global positioning system incorporating same | |
AU2009330687A1 (en) | Navigation receiver and method for combined use of a standard RTK system and a global carrier-phase differential positioning system | |
US5796773A (en) | Compensating for multipath errors, ionospheric time delay and noise errors in differential satellite positioning systems | |
EP2002279A1 (en) | Satellite ephemeris error | |
KR102314678B1 (en) | Method and Apparatus for GNSS multi-path mitigation | |
US20070085735A1 (en) | Spacecraft position monitoring and control | |
RU2181490C2 (en) | Device and process of radio positioning | |
WO2002069507A3 (en) | System and method for computing navigation information in the presence of interference | |
Kee et al. | Wide area differential GPS (WADGPS): Future navigation system | |
Kee et al. | Wide area differential GPS as a future navigation system in the US | |
US20040193373A1 (en) | Autonomous navigation error correction | |
US20060071849A1 (en) | Tactical all weather precision guidance and navigation system | |
US5914686A (en) | Utilization of exact solutions of the pseudorange equations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: PAT HLDR NO LONGER CLAIMS SMALL ENT STAT AS INDIV INVENTOR (ORIGINAL EVENT CODE: LSM1); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
FPAY | Fee payment |
Year of fee payment: 8 |
|
REMI | Maintenance fee reminder mailed | ||
LAPS | Lapse for failure to pay maintenance fees | ||
STCH | Information on status: patent discontinuation |
Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362 |
|
FP | Lapsed due to failure to pay maintenance fee |
Effective date: 20090415 |