US6134484A - Method and apparatus for maintaining the integrity of spacecraft based time and position using GPS - Google Patents

Method and apparatus for maintaining the integrity of spacecraft based time and position using GPS Download PDF

Info

Publication number
US6134484A
US6134484A US09/493,482 US49348200A US6134484A US 6134484 A US6134484 A US 6134484A US 49348200 A US49348200 A US 49348200A US 6134484 A US6134484 A US 6134484A
Authority
US
United States
Prior art keywords
gps
kalman filter
integrity
fault vector
spacecraft
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
US09/493,482
Inventor
George Jeffrey Geier
Roger Charles Hart
Gaylin David Brenchley
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Dynamics Mission Systems Inc
Original Assignee
Motorola Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Motorola Inc filed Critical Motorola Inc
Priority to US09/493,482 priority Critical patent/US6134484A/en
Assigned to MOTOROLA, INC reassignment MOTOROLA, INC ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: BRENCHLEY, GAYLIN DAVID, GEIER, GEORGE JEFFREY, HART, ROGER CHARLES
Application granted granted Critical
Publication of US6134484A publication Critical patent/US6134484A/en
Assigned to GENERAL DYNAMICS DECISION SYSTEMS, INC. reassignment GENERAL DYNAMICS DECISION SYSTEMS, INC. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: MOTOROLA, INC.
Assigned to GENERAL DYNAMICS C4 SYSTEMS, INC. reassignment GENERAL DYNAMICS C4 SYSTEMS, INC. MERGER AND CHANGE OF NAME Assignors: GENERAL DYNAMICS DECISION SYSTEMS, INC.
Assigned to GENERAL DYNAMICS ADVANCED INFORMATION SYSTEMS, INC. reassignment GENERAL DYNAMICS ADVANCED INFORMATION SYSTEMS, INC. MERGER (SEE DOCUMENT FOR DETAILS). Assignors: GENERAL DYNAMICS C4 SYSTEMS, INC.
Assigned to GENERAL DYNAMICS MISSION SYSTEMS, INC. reassignment GENERAL DYNAMICS MISSION SYSTEMS, INC. MERGER (SEE DOCUMENT FOR DETAILS). Assignors: GENERAL DYNAMICS ADVANCED INFORMATION SYSTEMS, INC.
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/20Integrity monitoring, fault detection or fault isolation of space segment
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/14Relay systems
    • H04B7/15Active relay systems
    • H04B7/185Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851Systems using a satellite or space-based relay
    • H04B7/18519Operations control, administration or maintenance

Definitions

  • This invention relates to a method and apparatus for maintaining the integrity of precise time and position derived from GPS. More particularly, the present invention relates to method and apparatus for isolating GPS satellite failure and removing the satellite from the time and position computations more quickly and reliably.
  • GPS Global Positioning System
  • a spacecraft using GPS for navigation includes a GPS receiver with a plurality of channels (12 or more), each dedicated to tracking a different GPS satellite.
  • the spacecraft is moving in an orbit or path different than the satellites in the GPS constellation which means that the GPS satellites being tracked are continuously changing
  • SA Selective Availablity
  • FIG. 1 is a graphical representation of "along track position error" for nominal GPS conditions
  • FIG. 2 is a graphical representation of "along track position error" for a simulated GPS satellite failure
  • FIG. 3 is a block diagram illustrating major processing blocks of apparatus for maintaining the integrity of spacecraft based time and position using GPS in accordance with the present invention.
  • FIG. 4 is an orbital diagram illustrating the operation of the apparatus of FIG. 3.
  • GPS Global Positioning System
  • the GPS receiver signal processing function acquires and tracks all of the GPS satellites within view. For each GPS satellite which is tracked, pseudo range and carrier phase measurements are generated, in addition to the demodulated data which constitutes the ephemeris for each GPS Satellite and enables calculation of the position and velocity of the spacecraft.
  • the set of available pseudo ranges and the set of available delta ranges are input to a Receiver Autonomous Integrity Monitor (RAIM) processor, in addition to the satellite positions computed from the ephemeris.
  • RAIM algorithms allow the GPS receiver to automatically detect and perhaps isolate and remove failed GPS satellites from inclusion in the resulting GPS position and time solution.
  • Kalman filters developed in the 1960s, have been used extensively in navigation and are well known in the art.
  • the Kalman filter is a recursive algorithm that uses physical observations (Doppler, range, angles, etc.) to obtain an optimal estimate of desired navigation parameters.
  • the key attraction of the Kalman filter is that it depends only on the previous navigation update, requiring little processing time and storage at each update and is thus, well suited to real-time applications in embedded microprocessors.
  • a Kalman filter algorithm is used to reduce the nominal errors in the GPS measurements.
  • the Kalman filter utilizes a high precision model of the spacecraft orbital dynamics, which permits it to filter the long time constant error variations from SA. Because the resultant filtered accuracy is significantly improved relative to nominal, unfiltered GPS position and timing accuracy, ensuring that the requirements are kept when GPS satellites fail represents a challenge for the design of a RAIM algorithm.
  • the Kalman filter was designed to remove the errors induced by Selective Availability (SA) using a detailed model for the spacecraft orbital and maneuver dynamics and already provides a measure of isolation from GPS satellite failures.
  • SA Selective Availability
  • the Kalman filter residual test when implemented as the usual "three sigma" test, should be effective in screening measurements from GPS satellites which have failed, i.e. whose errors are already large relative to the normal expected error levels from SA.
  • the Kalman filter residual test may not be effective in detecting a failing satellite, i.e., one which is slowly drifting to a large error condition, as the Kalman filter may absorb some of the error in its estimates and so not be able to detect it in its residuals.
  • FIG. 1 a graphical representation of "along track position error" for nominal GPS conditions is illustrated. Position error in meters is represented along the vertical axis and time in seconds is represented along the horizontal axis.
  • the transient errors induced in the system illustrated at the beginning of the graphs of FIGS. 1 and 2, are caused by a spacecraft thermister firing (a very low level, but persistent ion thernister). The thermister firing occurs during the interval spanned by the two vertical lines, designated 9, on each of FIGS. 1 and 2.
  • FIG. 2 a graphical representation of "along track position error" for a simulated GPS satellite failure is illustrated.
  • the error transient in the along-track direction from the GPS satellite failure is ramping at a rate of 0.05 m/sec resulting in a peak error of nearly seventy meters.
  • a marriage of the RAIM processor and the Kalman filter, which will be explained in more detail presently, is used in the preferred embodiment of the present invention to reduce the magnitude of the transients.
  • FIG. 3 a block diagram is illustrated showing major processing blocks of apparatus 10 for maintaining the integrity of spacecraft based time and position using GPS in accordance with the present invention.
  • Apparatus 10 includes a standard space borne GPS receiver 11, which is connected to supply GPS time and position signals to a RAIM processor 12.
  • RAIM processor 12 supplies signals to a Kalman filter 14, which has an output 15 on which it supplies the navigation information.
  • a spacecraft acceleration model 16 supplies information to Kalman filter 14 and an open loop propagator receives processed information from the Kalman filter and supplies continuous position information to RAIM processor 12.
  • GPS receiver signal processing 11 acquires and tracks the GPS satellites within view in a well known manner. For each GPS satellite which is tracked, pseudo range and carrier phase measurements are generated, in addition to the demodulated data which constitutes the ephemeris from each GPS satellite and enables calculation of the position and velocity of the spacecraft.
  • the set of available pseudo ranges and the set of available delta ranges are input to RAIM processing 12, in addition to the GPS satellite positions computed from the ephemeris.
  • RAIM processing 12 performs two distinct sets of calculations.
  • a first set of calculations is designed to reduce the impact of a GPS satellite with a slowly developing failure condition (e.g. a GPS satellite clock which is drifting abnormally).
  • a second set of calculations is designed to isolate and remove the failed GPS satellite from further processing by Kalman filter 14.
  • fault vectors are used to compute the "a posteriori statistic of unit variance" or "unit variance" in the parlance of statistics.
  • This statistic is a measure of how the measurement set (whether pseudo ranges or delta ranges) conforms to its expected error statistics, a value less than or equal to unity implies that the measurement error variance is set conservatively. As the unit error variance becomes larger, however, there is greater evidence that off-nominal error conditions exist and it is , therefore, used to scale up the error variances assumed by Kalman filter 14. In this way, Kalman filter 14 is desensitized to off-nominal error conditions (i.e., including a failed satellite).
  • apparatus 10 is carried by a spacecraft at a position designated 20 in an orbit 21.
  • a spacecraft position designated 20 in an orbit 21 is carried by a spacecraft at a position designated 20 in an orbit 21.
  • an "open-loop" propagation of the estimated spacecraft position and velocity is performed.
  • These "open-loop" estimates are independent of the GPS pseudo range and carrier phase measurements generated by the GPS receiver during the time interval.
  • the second set of calculations in apparatus 10 forms an "extended" fault vector, formed using pseudo and delta range measurements generated by the receiver at spacecraft position 23, and the open-loop position and velocity.
  • the extended fault vector augments the GPS-only pseudo range and delta range measurement sets with position and velocity estimates propagated open-loop (i.e. without processing GPS measurement data) from Kalman filter 14 estimated position and velocity, using spacecraft acceleration model 16.
  • the period of the open-loop propagation i.e. the time interval between spacecraft positions 20 and 23 in FIG. 4
  • SA Selective Availability
  • SPS Standard Positioning Service
  • an open-loop propagation is initiated every five minutes from the current best estimate of the trajectory maintained by Kalman filter 14 and used to construct the extended fault vectors (one for the set of pseudo ranges and propagated position and a second for the set of delta ranges and propagated velocity).
  • the open-loop propagation process operates in parallel with the normal Kalman filter algorithm with the only interaction, other than the initialization, being through the action of the extended fault vector.
  • the extended fault vector is then used as a conventional fault vector to support failure isolation. Since the susceptibility of Kalman filter 14 to failures that have not been isolated is reduced by the deweighting in response to the unit variance calculation, the extended fault vector calculations can be done at a reduced rate.
  • GPS receiver 11, RAIM processor 12 and Kalman filter 14 continuously process GPS information to compensate for transients and the like, generally as illustrated in FIG. 1.
  • the combination of RAIM processor 12 and Kalman filter 14 generally isolates and removes failed GPS satellites from the process. However, this combination may still be slow in isolating and removing failing satellites which are slowly drifting to a large error condition. Therefore, in this preferred embodiment very accurate predicted position information is simultaneously provided using orbital acceleration information from the Kalman filter through the use of the extended fault vector.
  • the method includes two features: computation of an extended fault vector for use in the RAIM algorithm for satellite failure detection and isolation; and parallel use of the RAIM algorithm outputs to de-weight potentially faulty satellites in the Kalman filter solution as the satellites failure develops (i.e. prior to failure isolation).
  • the extended fault vector concept augments the usual (GPS only) fault vector with position information predicted using the orbital acceleration information from Kalman filter 14. This additional information enhances the failure detection properties of the RAIM algorithm so that GPS satellite failures are isolated and removed more quickly and reliably.
  • the following description of apparatus 10 first describes the use of the RAIM function in adjusting the assumed error variances in Kalman filter 14 and then describes the computation of the external fault vector.
  • the fault vector, f can be derived from a redundant set of equations in a procedure described in two articles: a first by Sturza, M. A. and entitled “NAVIGATION SYSTEM INTEGRITY MONITORING USING REDUNDANT MEASUREMENTS", in the Journal of the Institute of Navigation, Vol. 35, No. 4, Winter 1988-89; and a second by Geier, G. J., et al.
  • fault vector f is equal to the number of measurements. Also, it should be noted that separate fault vectors are computed, using equations (1) and (2) below, for the set of redundant pseudo ranges and the set of redundant delta ranges.
  • S is the fault vector covariance matrix
  • R is the pseudo (or delta) range measurement error covariance
  • H is the gradient matrix corresponding to the redundant measurement set
  • I is the identity matrix
  • T denotes vector transpose
  • n DOF is the number of degrees of freedom.
  • the number of degrees of freedom above is simply the number of measurements in excess of the number of unknowns (four for position and velocity determination with GPS). Given that a unit variance has been computed separately for the set of pseudo ranges (PR) and delta ranges (DR), the following tests are applied to each:
  • ⁇ 2 PR and ⁇ 2 DR are the variances used by the Kalman filter.
  • the thresholds which control the operation of the de-weighting, T PR and T DR can be set based on an allowable false alarm rate. Generally, low false alarm rates are desirable, subject to tolerable probabilities of missed detection. However, the penalty to performance associated with a false alarm is insignificant in this application, since the spacecraft acceleration model 16 used by Kalman filter 14 is highly accurate. Therefore, a relatively low threshold can be used and will minimize the impact of failing satellites.
  • the extended fault vector is constructed (at a lower rate). As described previously, every five minutes, an open-loop propagation from the Kalman filter estimated position and velocity is initiated. A highly numerical integration method (e.g. fourth order Runge-Kutta) is used to reduce the numerical error to insignificant levels. At the end of the propagation interval, three dimensional position and velocity estimates, denoted P orb and v orb , are available. Each is used to construct the appropriate extended fault vector, computed from an external measurement set Z ePR and Z eDR :
  • P is the Kalman filter estimated position
  • Z PR is the set of pseudo range residuals.
  • V is the Kalman filter estimated velocity
  • Z DR is the set of delta range residuals.
  • H PR is the pseudo range measurement gradient matrix
  • H DR is the delta range measurement gradient matrix
  • H orb has the following simple form: ##EQU1## where: I is a three dimensional identity matrix; and
  • 0 is a three dimensional zero vector.
  • P orbp is the covariance matrix associated with the open-loop propagation of position by the Kalman filter.
  • P orbv is the covariance matrix associated with the open-loop propagation of velocity by the Kalman filter.
  • Equation (15) represents factorization of the covariance matrix, where U is an upper diagonal matrix, and D is a diagonal matrix. This decomposition of the covariance matrix permits the square root and the inverse square root of the covariance matrices to be found as:
  • the apparatus includes a GPS receiver, a Receiver Autonomous Integrity Monitor (RAIM) processor coupled to receive signals from the GPS receiver, and a Kalman filter coupled to receive signals from the RAIM processor.
  • GPS signals from GPS satellites are processed in a normal fashion to provide navigation information including position, velocity, time, etc.
  • a spacecraft acceleration model is coupled to provide information to the Kalman filter to generate an extended fault vector and an open loop propagator is coupled to receive extended fault vector information from the Kalman filter and supply open-loop position and velocity errors to the RAIM processor.
  • the RAIM processor and Kalman filter use the open-loop position and velocity errors to isolate and remove failing GPS satellites from consideration. Accordingly, a method and apparatus are disclosed which are capable of quickly and reliably isolating and removing failed or failing GPS satellites from the GPS position and time calculations used in navigation of a spacecraft.

Abstract

A method and apparatus for maintaining the integrity of spacecraft based time and position using GPS (11) is disclosed. A GPS receiver is provided with a Receiver Autonomous Integrity Monitor (RAIM) processor (12) and a Kalman filter (14). GPS signals are received from at least four GPS satellites and processed in the GPS receiver to provide time and position information. Periodically, an extended fault vector is formed in the Kalman filter with propagated open-loop position and velocity measurements. The extended fault vector is used to isolate and remove failing GPS satellites.

Description

FIELD OF THE INVENTION
This invention relates to a method and apparatus for maintaining the integrity of precise time and position derived from GPS. More particularly, the present invention relates to method and apparatus for isolating GPS satellite failure and removing the satellite from the time and position computations more quickly and reliably.
BACKGROUND OF THE INVENTION
As more and more spacecraft and communications systems rely on the Global Positioning System (GPS) for their timing source, it is clear that the reliability of the spacecraft or communications link is a function of the entire GPS system. Even though GPS has demonstrated superb reliability and availability, a few operational glitches have occurred that affected the time and position solutions (i.e. the time and position output). At least one of these has contributed to a complete shutdown of a prototype TDMA system until the problem was cleared by the GPS Control Segment some 30 minutes later (the affected satellite was set unhealthy). A major and somewhat unadvertised flaw in the GPS system design is that there is no apparent way for the Control Segment to immediately remove a bad satellite from the constellation, the process time to do so is typically between 5 and 45 minutes.
Many spacecraft orbiting the earth use information derived from GPS for navigation and/or to maintain them in their specific orbit. A spacecraft using GPS for navigation includes a GPS receiver with a plurality of channels (12 or more), each dedicated to tracking a different GPS satellite. Generally, the spacecraft is moving in an orbit or path different than the satellites in the GPS constellation which means that the GPS satellites being tracked are continuously changing
The detection of failed GPS satellites is generally limited by the presence of Selective Availablity (SA). As is well known by experienced GPS users, SA, imposed on civil users by the Department of Defense, severely limits navigation accuracy. In addition, it also limits conventional methods of failure detection as, for example, referenced in Brown, "A Baseline RAIM Scheme and a Note on the Equivalence of Three RAIM Methods", Journal of the Institute of Navigation, Fall 1992. The fundamental limitation imposed by SA limits failure detection to error magnitudes which are large in comparison to SA. However, information provided by a failed or failing GPS satellite can seriously affect the navigation of a spacecraft relying on GPS information.
In U. S. Pat. No. 5,583,774, issued Dec. 10, 1996 and entitled "Assured-Integrity Monitored-Extrapolation Navigation Device", an integrity monitor is developed for GPS using sensor augmentation. In this device, an Inertial Navigation System (INS) is used to maintain a highly accurate trajectory for a period of time sufficiently long (e.g. 30 minutes) to enable isolation of failed satellites by the GPS Control Segment. A bank of Kalman filters, using the GPS-calibrated INS as a reference trajectory, is used to estimate a range bias for each satellite, which is used to confirm its health. Only satellites passing this health test can be used in the calibration. The major problem with this system is the need for external sensors and banks of Kalman filters, which greatly increases the size and complexity of the system.
Accordingly, it is highly desirable to provide a method and apparatus which is capable of quickly and reliably isolating and removing failed or failing GPS satellites from the GPS position and time calculations used in navigation of a spacecraft.
BRIEF DESCRIPTION OF THE DRAWINGS
Referring to the drawings:
FIG. 1 is a graphical representation of "along track position error" for nominal GPS conditions;
FIG. 2 is a graphical representation of "along track position error" for a simulated GPS satellite failure;
FIG. 3 is a block diagram illustrating major processing blocks of apparatus for maintaining the integrity of spacecraft based time and position using GPS in accordance with the present invention; and
FIG. 4 is an orbital diagram illustrating the operation of the apparatus of FIG. 3.
DESCRIPTION OF THE PREFERRED EMBODIMENT
Spacecraft relying on the Global Positioning System (GPS) for navigation (e.g. position and time) informnation include a GPS receiver having a plurality of channels, each of which tracks a different GPS satellite transmission. Generally, the GPS receiver signal processing function acquires and tracks all of the GPS satellites within view. For each GPS satellite which is tracked, pseudo range and carrier phase measurements are generated, in addition to the demodulated data which constitutes the ephemeris for each GPS Satellite and enables calculation of the position and velocity of the spacecraft.
The set of available pseudo ranges and the set of available delta ranges (computed from differenced carrier phase data) are input to a Receiver Autonomous Integrity Monitor (RAIM) processor, in addition to the satellite positions computed from the ephemeris. RAIM algorithms allow the GPS receiver to automatically detect and perhaps isolate and remove failed GPS satellites from inclusion in the resulting GPS position and time solution.
Kalman filters, developed in the 1960s, have been used extensively in navigation and are well known in the art. The Kalman filter is a recursive algorithm that uses physical observations (Doppler, range, angles, etc.) to obtain an optimal estimate of desired navigation parameters. The key attraction of the Kalman filter is that it depends only on the previous navigation update, requiring little processing time and storage at each update and is thus, well suited to real-time applications in embedded microprocessors. In maintaining precise position and time accuracy in orbiting applications, a Kalman filter algorithm is used to reduce the nominal errors in the GPS measurements. The Kalman filter utilizes a high precision model of the spacecraft orbital dynamics, which permits it to filter the long time constant error variations from SA. Because the resultant filtered accuracy is significantly improved relative to nominal, unfiltered GPS position and timing accuracy, ensuring that the requirements are kept when GPS satellites fail represents a challenge for the design of a RAIM algorithm.
In spacecraft such as that being discussed herein, the Kalman filter was designed to remove the errors induced by Selective Availability (SA) using a detailed model for the spacecraft orbital and maneuver dynamics and already provides a measure of isolation from GPS satellite failures. The Kalman filter residual test, when implemented as the usual "three sigma" test, should be effective in screening measurements from GPS satellites which have failed, i.e. whose errors are already large relative to the normal expected error levels from SA. However, the Kalman filter residual test may not be effective in detecting a failing satellite, i.e., one which is slowly drifting to a large error condition, as the Kalman filter may absorb some of the error in its estimates and so not be able to detect it in its residuals.
Simulations have indicated that the Kalman filter solution can exhibit relatively large error transients when a GPS satellite failure occurs during a maneuver period of the spacecraft. Referring to FIG. 1, a graphical representation of "along track position error" for nominal GPS conditions is illustrated. Position error in meters is represented along the vertical axis and time in seconds is represented along the horizontal axis. Here it can be seen that minor perturbations or transients in position, with properly operating GPS satellites, are quickly compensated by either the RAIM processor or the Kalman filter, resulting in a peak error of roughly 15 meters. The transient errors induced in the system, illustrated at the beginning of the graphs of FIGS. 1 and 2, are caused by a spacecraft thermister firing (a very low level, but persistent ion thernister). The thermister firing occurs during the interval spanned by the two vertical lines, designated 9, on each of FIGS. 1 and 2.
Referring to FIG. 2, a graphical representation of "along track position error" for a simulated GPS satellite failure is illustrated. In this example, the error transient in the along-track direction from the GPS satellite failure is ramping at a rate of 0.05 m/sec resulting in a peak error of nearly seventy meters. A marriage of the RAIM processor and the Kalman filter, which will be explained in more detail presently, is used in the preferred embodiment of the present invention to reduce the magnitude of the transients.
Turning now to FIG. 3, a block diagram is illustrated showing major processing blocks of apparatus 10 for maintaining the integrity of spacecraft based time and position using GPS in accordance with the present invention. Here it should be noted that while separate blocks are illustrated for the various functions, all of the processing functions can be performed in a single processor embedded in the GPS receiver. Apparatus 10 includes a standard space borne GPS receiver 11, which is connected to supply GPS time and position signals to a RAIM processor 12. RAIM processor 12 supplies signals to a Kalman filter 14, which has an output 15 on which it supplies the navigation information. Additionally, a spacecraft acceleration model 16 supplies information to Kalman filter 14 and an open loop propagator receives processed information from the Kalman filter and supplies continuous position information to RAIM processor 12.
Basically, GPS receiver signal processing 11 acquires and tracks the GPS satellites within view in a well known manner. For each GPS satellite which is tracked, pseudo range and carrier phase measurements are generated, in addition to the demodulated data which constitutes the ephemeris from each GPS satellite and enables calculation of the position and velocity of the spacecraft.
The set of available pseudo ranges and the set of available delta ranges (computed from differenced carrier phase data) are input to RAIM processing 12, in addition to the GPS satellite positions computed from the ephemeris. RAIM processing 12 performs two distinct sets of calculations. A first set of calculations is designed to reduce the impact of a GPS satellite with a slowly developing failure condition (e.g. a GPS satellite clock which is drifting abnormally). A second set of calculations is designed to isolate and remove the failed GPS satellite from further processing by Kalman filter 14. In the first set of calculations, fault vectors are used to compute the "a posteriori statistic of unit variance" or "unit variance" in the parlance of statistics. This statistic is a measure of how the measurement set (whether pseudo ranges or delta ranges) conforms to its expected error statistics, a value less than or equal to unity implies that the measurement error variance is set conservatively. As the unit error variance becomes larger, however, there is greater evidence that off-nominal error conditions exist and it is , therefore, used to scale up the error variances assumed by Kalman filter 14. In this way, Kalman filter 14 is desensitized to off-nominal error conditions (i.e., including a failed satellite).
Referring additionally to FIG. 4, an orbital diagram showing the operation of the apparatus of FIG. 3 is illustrated. In this specific example, apparatus 10 is carried by a spacecraft at a position designated 20 in an orbit 21. During the period of time (300 seconds) from spacecraft position 20 to spacecraft position 23, an "open-loop" propagation of the estimated spacecraft position and velocity is performed. These "open-loop" estimates are independent of the GPS pseudo range and carrier phase measurements generated by the GPS receiver during the time interval. The second set of calculations in apparatus 10 forms an "extended" fault vector, formed using pseudo and delta range measurements generated by the receiver at spacecraft position 23, and the open-loop position and velocity. Here it should be understood that FIG. 4 simply illustrates an instantaneous example and, as the spacecraft moves from one position (e.g. position 20) to the next position (e.g. position 23), the next open-loop position and extended fault vector is computed. The extended fault vector (and all of the following vectors) is used to isolate failed pseudo range and/or delta range measurements and remove them from further processing. During the open-loop propagation, the normal measurement processing steps of the Kalman filter, indicated by time ticks 22 (nominally occurring every 10-30 seconds), continue.
The extended fault vector augments the GPS-only pseudo range and delta range measurement sets with position and velocity estimates propagated open-loop (i.e. without processing GPS measurement data) from Kalman filter 14 estimated position and velocity, using spacecraft acceleration model 16. The period of the open-loop propagation (i.e. the time interval between spacecraft positions 20 and 23 in FIG. 4) is the minimum time interval required to effectively decorrelate the open-loop position and velocity errors from the set of pseudo range and delta range measurement errors. Since Selective Availability (SA) is the dominant error source for GPS Standard Positioning Service (SPS) receivers, an interval of roughly 300 seconds (three SA time constants) is appropriate. Thus, an open-loop propagation is initiated every five minutes from the current best estimate of the trajectory maintained by Kalman filter 14 and used to construct the extended fault vectors (one for the set of pseudo ranges and propagated position and a second for the set of delta ranges and propagated velocity). The open-loop propagation process operates in parallel with the normal Kalman filter algorithm with the only interaction, other than the initialization, being through the action of the extended fault vector. The extended fault vector is then used as a conventional fault vector to support failure isolation. Since the susceptibility of Kalman filter 14 to failures that have not been isolated is reduced by the deweighting in response to the unit variance calculation, the extended fault vector calculations can be done at a reduced rate.
In the operation of apparatus 10, GPS receiver 11, RAIM processor 12 and Kalman filter 14 continuously process GPS information to compensate for transients and the like, generally as illustrated in FIG. 1. The combination of RAIM processor 12 and Kalman filter 14 generally isolates and removes failed GPS satellites from the process. However, this combination may still be slow in isolating and removing failing satellites which are slowly drifting to a large error condition. Therefore, in this preferred embodiment very accurate predicted position information is simultaneously provided using orbital acceleration information from the Kalman filter through the use of the extended fault vector. The method includes two features: computation of an extended fault vector for use in the RAIM algorithm for satellite failure detection and isolation; and parallel use of the RAIM algorithm outputs to de-weight potentially faulty satellites in the Kalman filter solution as the satellites failure develops (i.e. prior to failure isolation). The extended fault vector concept augments the usual (GPS only) fault vector with position information predicted using the orbital acceleration information from Kalman filter 14. This additional information enhances the failure detection properties of the RAIM algorithm so that GPS satellite failures are isolated and removed more quickly and reliably. Also, contamination of the Kalman filter before a failing GPS satellite is removed from consideration by the RAIM algorithm or the Kalman filter residual test is prevented as the failure develops and is not yet isolated, by using the RAIM algorithm outputs to de-weight possible failed GPS satellites and scaling the error variances which Kalman filter 14 assigns to them.
The following description of apparatus 10 first describes the use of the RAIM function in adjusting the assumed error variances in Kalman filter 14 and then describes the computation of the external fault vector. The fault vector, f, can be derived from a redundant set of equations in a procedure described in two articles: a first by Sturza, M. A. and entitled "NAVIGATION SYSTEM INTEGRITY MONITORING USING REDUNDANT MEASUREMENTS", in the Journal of the Institute of Navigation, Vol. 35, No. 4, Winter 1988-89; and a second by Geier, G. J., et al. entitled "THE APPLICATION OF INTEGRITY MONITORING TO OFFSHORE SEISMIC SURVEY OPERATIONS" in the Journal ofthe Institute of Navigation, Vol. 39, No. 3, Fall 1992, both of which are incorporated herein by reference. The dimension of fault vector f is equal to the number of measurements. Also, it should be noted that separate fault vectors are computed, using equations (1) and (2) below, for the set of redundant pseudo ranges and the set of redundant delta ranges.
f=SR.sup.-1/2 z                                            (1)
S=I-R.sup.-1/2 H(H.sup.T R.sup.-1 H).sup.-1 H.sup.T R.sup.-1/2(2)
where:
S is the fault vector covariance matrix;
R is the pseudo (or delta) range measurement error covariance;
H is the gradient matrix corresponding to the redundant measurement set; and
I is the identity matrix.
The unit variance, u, is found directly from f:
u=f.sup.T f/n.sub.DOF                                      (3)
where:
T denotes vector transpose; and
nDOF is the number of degrees of freedom.
The number of degrees of freedom above is simply the number of measurements in excess of the number of unknowns (four for position and velocity determination with GPS). Given that a unit variance has been computed separately for the set of pseudo ranges (PR) and delta ranges (DR), the following tests are applied to each:
if(u.sub.PR >T.sub.PR)□.sup.2.sub.PR =u.sub.PR □.sup.2.sub.PR                                 (4)
if(u.sub.DR >T.sub.DR)□.sup.2.sub.DR =u.sub.DR □.sup.2.sub.DR                                 (5)
where:
2 PR and □2 DR are the variances used by the Kalman filter.
The thresholds which control the operation of the de-weighting, TPR and TDR, can be set based on an allowable false alarm rate. Generally, low false alarm rates are desirable, subject to tolerable probabilities of missed detection. However, the penalty to performance associated with a false alarm is insignificant in this application, since the spacecraft acceleration model 16 used by Kalman filter 14 is highly accurate. Therefore, a relatively low threshold can be used and will minimize the impact of failing satellites.
In addition to the de-weighting described above, the extended fault vector is constructed (at a lower rate). As described previously, every five minutes, an open-loop propagation from the Kalman filter estimated position and velocity is initiated. A highly numerical integration method (e.g. fourth order Runge-Kutta) is used to reduce the numerical error to insignificant levels. At the end of the propagation interval, three dimensional position and velocity estimates, denoted Porb and vorb, are available. Each is used to construct the appropriate extended fault vector, computed from an external measurement set ZePR and ZeDR :
Z.sub.ePR.sup.T =[Z.sub.PR.sup.T (P-P.sub.orb).sup.T ]     (6)
where:
P is the Kalman filter estimated position; and
ZPR is the set of pseudo range residuals.
Z.sub.eDR.sup.T =[Z.sub.DR.sup.T (V-V.sub.orb).sup.T ]     (7)
where:
V is the Kalman filter estimated velocity; and
ZDR is the set of delta range residuals.
Extended versions of the measurement gradient matrix H and the measurement variance matrix R are required:
H.sub.ePR.sup.T =[H.sub.PR.sup.T H.sub.orb ]               (8)
where HPR is the pseudo range measurement gradient matrix.
H.sub.eDR.sup.T =[H.sub.DR.sup.T H.sub.orb ]               (9)
where HDR is the delta range measurement gradient matrix.
Because position and velocity are directly estimated in the open-loop propagation, Horb has the following simple form: ##EQU1## where: I is a three dimensional identity matrix; and
0 is a three dimensional zero vector. ##EQU2## where: Porbp is the covariance matrix associated with the open-loop propagation of position by the Kalman filter. ##EQU3## where: Porbv is the covariance matrix associated with the open-loop propagation of velocity by the Kalman filter.
Note that the fault vector computation (eqs. 1 and 2) require that the R matrices (given by Eqs. 11 and 12) be inverted, and also that their square roots be taken. Because of the block-diagonal form of the extended matrices, this is relatively easy and straightforward, as given below: ##EQU4##
Inverting RPR and RDR is straightforward, since they are diagonal matrices. However, inverting Porbp and Porbv requires more work. The Kalman filter can be mechanized in the popular U-D mechanization, as explained in an article by Bierman, G. J., "Measurement Updating Using the U-D Factorization", Proceedings of the 1975 IEEE, Conference on Decision and Control, incorporated herein by reference. If the Kalman filter is mechanized as described, however, the inverted calculations have already been done, as indicated below:
P=UDU.sup.T                                                (15)
Equation (15) represents factorization of the covariance matrix, where U is an upper diagonal matrix, and D is a diagonal matrix. This decomposition of the covariance matrix permits the square root and the inverse square root of the covariance matrices to be found as:
P.sup.1/2 =UD.sup.1/2                                      (16)
P-.sup.1/2 =UD.sup.-1/2                                    (17)
Thus, a new and improved method and apparatus for maintaining the integrity of spacecraft based time and position using GPS is disclosed. In a preferred embodiment the apparatus includes a GPS receiver, a Receiver Autonomous Integrity Monitor (RAIM) processor coupled to receive signals from the GPS receiver, and a Kalman filter coupled to receive signals from the RAIM processor. GPS signals from GPS satellites are processed in a normal fashion to provide navigation information including position, velocity, time, etc. A spacecraft acceleration model is coupled to provide information to the Kalman filter to generate an extended fault vector and an open loop propagator is coupled to receive extended fault vector information from the Kalman filter and supply open-loop position and velocity errors to the RAIM processor. The RAIM processor and Kalman filter use the open-loop position and velocity errors to isolate and remove failing GPS satellites from consideration. Accordingly, a method and apparatus are disclosed which are capable of quickly and reliably isolating and removing failed or failing GPS satellites from the GPS position and time calculations used in navigation of a spacecraft.
While we have shown and described specific embodiments of the present invention, further modifications and improvements will occur to those skilled in the art. We desire it to be understood, therefore, that this invention is not limited to the particular forms shown and we intend in the appended claims to cover all modifications that do not depart from the spirit and scope of this invention.

Claims (10)

What is claimed is:
1. A method for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) comprising the steps of:
providing a GPS receiver coupled to a Receiver Autonomous Integrity Monitor (RAIM) processor and a Kalman filter;
receiving GPS signals from GPS satellites and processing the GPS signals in the GPS receiver to provide time and position information;
periodically forming an extended fault vector in the Kalman filter with propagated open-loop position and velocity measurements;
the step of periodically forming the extended fault vector includes inputting a spacecraft acceleration model into the Kalman filter; and
using the extended fault vector to isolate and remove failing GPS satellites.
2. A method for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) as claimed in claim 1 wherein the step of periodically forming the extended fault vector includes forming the extended fault vector at intervals which are long relative to Selective Availability (SA).
3. A method for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) as claimed in claim 1 wherein the step of periodically forming the extended fault vector includes forming the extended fault vector at intervals of approximately 300 seconds.
4. A method for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) comprising the steps of:
providing a GPS receiver coupled to a Receiver Autonomous Integrity Monitor (RAIM) processor and a Kalman filter;
receiving GPS signals from GPS satellites and processing the GPS signals in the GPS receiver to provide time and position information;
computing unit variance in the RAIM processor and using the unit variance to desensitize the Kalman filter to off-nominal conditions;
periodically forming an extended fault vector in the Kalman filter with propagated open-loop position and velocity measurements;
the step of periodically forming the extended fault vector includes inputting a spacecraft acceleration model into the Kalman filter; and
using the extended fault vector to isolate failures of GPS satellites and failing GPS satellites.
5. A method for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) as claimed in claim 4 wherein the step of periodically forming the extended fault vector includes forming the extended fault vector at intervals which are long relative to Selective Availability (SA).
6. A method for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) as claimed in claim 4, wherein the step of periodically forming the extended fault vector includes forming the extended fault vector at intervals of approximately 300 seconds.
7. A method for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) as claimed in claim 4 wherein the step of computing unit variance includes computing pseudo ranges from the received time and position information.
8. A method for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) as claimed in claim 4 wherein the step of receiving GPS signals includes receiving carrier phase data and the step of computing unit variance includes computing delta ranges from differenced carrier phase data.
9. An apparatus for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) comprising:
a GPS receiver;
a Receiver Autonomous Integrity Monitor (RAIM) processor coupled to receive signals from the GPS receiver;
a Kalman filter coupled to receive signals from the RAIM processor;
a spacecraft acceleration model coupled to provide information to the Kalman filter to generate an extended fault vector;
an open loop propagator coupled to receive extended fault vector information from the Kalman filter and supplying open-loop position and velocity measurements to the RAIM processor; and
the RAIM processor, Kalman filter, spacecraft acceleration model, and open loop propagator are all functions performed in a single processor embedded in the GPS receiver.
10. An apparatus for maintaining the integrity of spacecraft based time and position using the Global Positioning System (GPS) as claimed in claim 9 wherein the GPS receiver includes at least four channels for tracking GPS satellites.
US09/493,482 2000-01-28 2000-01-28 Method and apparatus for maintaining the integrity of spacecraft based time and position using GPS Expired - Lifetime US6134484A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US09/493,482 US6134484A (en) 2000-01-28 2000-01-28 Method and apparatus for maintaining the integrity of spacecraft based time and position using GPS

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US09/493,482 US6134484A (en) 2000-01-28 2000-01-28 Method and apparatus for maintaining the integrity of spacecraft based time and position using GPS

Publications (1)

Publication Number Publication Date
US6134484A true US6134484A (en) 2000-10-17

Family

ID=23960382

Family Applications (1)

Application Number Title Priority Date Filing Date
US09/493,482 Expired - Lifetime US6134484A (en) 2000-01-28 2000-01-28 Method and apparatus for maintaining the integrity of spacecraft based time and position using GPS

Country Status (1)

Country Link
US (1) US6134484A (en)

Cited By (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6205400B1 (en) * 1998-11-27 2001-03-20 Ching-Fang Lin Vehicle positioning and data integrating method and system thereof
US6445344B1 (en) 1999-11-16 2002-09-03 Northrop Grumman Corporation Local area positioning system
US6480789B2 (en) * 2000-12-04 2002-11-12 American Gnc Corporation Positioning and proximity warning method and system thereof for vehicle
US20030216865A1 (en) * 2002-04-12 2003-11-20 Ensco, Inc. Inertial navigation system for mobile objects with constraints
US20040220733A1 (en) * 2003-04-29 2004-11-04 United Parcel Service Of America, Inc. Systems and methods for fault detection and exclusion in navigational systems
US20050015680A1 (en) * 2002-09-25 2005-01-20 Rubin Arthur Leonard Methods and apparatus for evaluating operational integrity of a data processing system using moment bounding
US20050131636A1 (en) * 2003-12-16 2005-06-16 Alcatel Device for a mobile terminal for determining position by filtering integrity data from an augmentation device
US20050162794A1 (en) * 2003-10-30 2005-07-28 Donnelly David S. System for mitigating radiation effects on componentry in space
US20060074558A1 (en) * 2003-11-26 2006-04-06 Williamson Walton R Fault-tolerant system, apparatus and method
US20060195262A1 (en) * 2004-09-17 2006-08-31 Alexander Draganov GPS accumulated delta range processing for navigation applications
US20060208941A1 (en) * 2005-03-15 2006-09-21 Honeywell International Inc. Adjusting processor clock information using a clock drift estimate
US20070050102A1 (en) * 2005-08-29 2007-03-01 Johnson Gary N Space based navigation system
US20070115171A1 (en) * 2005-11-18 2007-05-24 Rahman Mohammad A Methods and apparatus to detect and correct integrity failures in satellite positioning system receivers
US20070156338A1 (en) * 2004-02-13 2007-07-05 Jacques Coatantiec Device for monitoring the integrity of information delivered by a hybrid ins/gnss system
US20070279279A1 (en) * 2006-06-06 2007-12-06 Honeywell International Inc. System and method for enhancing the performance of satellite navigation receivers
US20090091495A1 (en) * 2007-10-09 2009-04-09 Honeywell International Inc. Gps receiver raim with slaved precision clock
US20090146873A1 (en) * 2007-12-07 2009-06-11 Honeywell International Inc. Navigation system with apparatus for detecting accuracy failures
US20090182495A1 (en) * 2008-01-15 2009-07-16 Honeywell International, Inc. Navigation system with apparatus for detecting accuracy failures
US20090182493A1 (en) * 2008-01-15 2009-07-16 Honeywell International, Inc. Navigation system with apparatus for detecting accuracy failures
US20090182494A1 (en) * 2008-01-15 2009-07-16 Honeywell International, Inc. Navigation system with apparatus for detecting accuracy failures
US7783425B1 (en) * 2005-06-29 2010-08-24 Rockwell Collins, Inc. Integrity-optimized receiver autonomous integrity monitoring (RAIM)
US20110210890A1 (en) * 2008-11-05 2011-09-01 Robert Andre Laine System for monitoring a satellite constellation of a positioning system
US20110254729A1 (en) * 2009-09-29 2011-10-20 Texas Instruments Incorporated Cross coupled positioning engine (pe) architecture for sensor integration in global navigation satellite system (gnss)
CN101799524B (en) * 2009-07-10 2013-01-30 中国测绘科学研究院 Method for autonomously monitoring receiver integrity of global navigation satellite system
US20130088387A1 (en) * 2011-10-07 2013-04-11 Electronics And Telecommunications Research Institute Apparatus and method for monitoring malfunctioning state of global positioning system (gps) satellite
US8706319B2 (en) 2011-12-16 2014-04-22 The Boeing Company Space positioning system
US8922426B1 (en) * 2010-02-26 2014-12-30 Google Inc. System for geo-location
CN105487088A (en) * 2015-09-12 2016-04-13 北京大学 RAIM algorithm in satellite navigation system based on Kalman filtering
US9341718B2 (en) 2012-09-07 2016-05-17 Honeywell International Inc. Method and system for providing integrity for hybrid attitude and true heading
WO2018059532A1 (en) * 2016-09-30 2018-04-05 华为技术有限公司 Integrated navigational data calculation method of time-delay observation system, apparatus and navigational device
CN115032668A (en) * 2022-06-15 2022-09-09 北京讯腾智慧科技股份有限公司 Desensitization method and device for satellite navigation positioning reference station observation data

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5583744A (en) * 1994-09-06 1996-12-10 Citizen Watch Co., Ltd. Portable computer with detachable battery pack
US5825326A (en) * 1996-07-09 1998-10-20 Interstate Electronics Corporation Real-time high-accuracy determination of integer ambiguities in a kinematic GPS receiver
US5923286A (en) * 1996-10-23 1999-07-13 Honeywell Inc. GPS/IRS global position determination method and apparatus with integrity loss provisions

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5583744A (en) * 1994-09-06 1996-12-10 Citizen Watch Co., Ltd. Portable computer with detachable battery pack
US5825326A (en) * 1996-07-09 1998-10-20 Interstate Electronics Corporation Real-time high-accuracy determination of integer ambiguities in a kinematic GPS receiver
US5923286A (en) * 1996-10-23 1999-07-13 Honeywell Inc. GPS/IRS global position determination method and apparatus with integrity loss provisions

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
An article entitled "A Baseline GPS RAIM Scheme And a Note on the Equivalence of Three RAIM Methods", by R. Grover Brown, Iowa State University, Ames, Iowa, 1992 from the Navigation: Journal of The institute of Navigation.
An article entitled A Baseline GPS RAIM Scheme And a Note on the Equivalence of Three RAIM Methods , by R. Grover Brown, Iowa State University, Ames, Iowa, 1992 from the Navigation: Journal of The institute of Navigation. *

Cited By (54)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6205400B1 (en) * 1998-11-27 2001-03-20 Ching-Fang Lin Vehicle positioning and data integrating method and system thereof
US6445344B1 (en) 1999-11-16 2002-09-03 Northrop Grumman Corporation Local area positioning system
US6480789B2 (en) * 2000-12-04 2002-11-12 American Gnc Corporation Positioning and proximity warning method and system thereof for vehicle
US20030216865A1 (en) * 2002-04-12 2003-11-20 Ensco, Inc. Inertial navigation system for mobile objects with constraints
US6826478B2 (en) 2002-04-12 2004-11-30 Ensco, Inc. Inertial navigation system for mobile objects with constraints
US20050015680A1 (en) * 2002-09-25 2005-01-20 Rubin Arthur Leonard Methods and apparatus for evaluating operational integrity of a data processing system using moment bounding
US7089452B2 (en) * 2002-09-25 2006-08-08 Raytheon Company Methods and apparatus for evaluating operational integrity of a data processing system using moment bounding
US20040220733A1 (en) * 2003-04-29 2004-11-04 United Parcel Service Of America, Inc. Systems and methods for fault detection and exclusion in navigational systems
US6856905B2 (en) 2003-04-29 2005-02-15 Garmin At, Inc. Systems and methods for fault detection and exclusion in navigational systems
US20050096844A1 (en) * 2003-04-29 2005-05-05 Garmin At, Inc. An Oregon Corporation Systems and methods for fault detection and exclusion in navigational systems
US6944541B2 (en) 2003-04-29 2005-09-13 Garmin At, Inc. Systems and methods for fault detection and exclusion in navigational systems
US20050162794A1 (en) * 2003-10-30 2005-07-28 Donnelly David S. System for mitigating radiation effects on componentry in space
US20060074558A1 (en) * 2003-11-26 2006-04-06 Williamson Walton R Fault-tolerant system, apparatus and method
US6970785B2 (en) 2003-12-16 2005-11-29 Alcatel Device for a mobile terminal for determining position by filtering integrity data from an augmentation device
EP1544636A1 (en) * 2003-12-16 2005-06-22 Alcatel Determination device by filtering of integrity data in a differential positioning system, for a mobile terminal
FR2863714A1 (en) * 2003-12-16 2005-06-17 Cit Alcatel POSITION DETERMINING DEVICE BY FILTERING INTEGRITY DATA OF AN INCREASE SYSTEM FOR A MOBILE TERMINAL
US20050131636A1 (en) * 2003-12-16 2005-06-16 Alcatel Device for a mobile terminal for determining position by filtering integrity data from an augmentation device
US20070156338A1 (en) * 2004-02-13 2007-07-05 Jacques Coatantiec Device for monitoring the integrity of information delivered by a hybrid ins/gnss system
US7409289B2 (en) * 2004-02-13 2008-08-05 Thales Device for monitoring the integrity of information delivered by a hybrid INS/GNSS system
US20060195262A1 (en) * 2004-09-17 2006-08-31 Alexander Draganov GPS accumulated delta range processing for navigation applications
US7490008B2 (en) * 2004-09-17 2009-02-10 Itt Manufacturing Enterprises, Inc. GPS accumulated delta range processing for navigation applications
US20060208941A1 (en) * 2005-03-15 2006-09-21 Honeywell International Inc. Adjusting processor clock information using a clock drift estimate
US7889125B2 (en) * 2005-03-15 2011-02-15 Honeywell International Inc. Adjusting processor clock information using a clock drift estimate
US7783425B1 (en) * 2005-06-29 2010-08-24 Rockwell Collins, Inc. Integrity-optimized receiver autonomous integrity monitoring (RAIM)
US20070050102A1 (en) * 2005-08-29 2007-03-01 Johnson Gary N Space based navigation system
US20070115171A1 (en) * 2005-11-18 2007-05-24 Rahman Mohammad A Methods and apparatus to detect and correct integrity failures in satellite positioning system receivers
US7501981B2 (en) 2005-11-18 2009-03-10 Texas Instruments Incorporated Methods and apparatus to detect and correct integrity failures in satellite positioning system receivers
US20070279279A1 (en) * 2006-06-06 2007-12-06 Honeywell International Inc. System and method for enhancing the performance of satellite navigation receivers
US20090091495A1 (en) * 2007-10-09 2009-04-09 Honeywell International Inc. Gps receiver raim with slaved precision clock
US7667644B2 (en) * 2007-10-09 2010-02-23 Honeywell International Inc. GPS receiver RAIM with slaved precision clock
US20100149025A1 (en) * 2007-10-09 2010-06-17 Honeywell International Inc. Gps receiver raim with slaved precision clock
US20090146873A1 (en) * 2007-12-07 2009-06-11 Honeywell International Inc. Navigation system with apparatus for detecting accuracy failures
US20090150074A1 (en) * 2007-12-07 2009-06-11 Honeywell International Inc. Navigation system with apparatus for detecting accuracy failures
US8014948B2 (en) 2007-12-07 2011-09-06 Honeywell International Inc. Navigation system with apparatus for detecting accuracy failures
US8019539B2 (en) 2007-12-07 2011-09-13 Honeywell International Inc. Navigation system with apparatus for detecting accuracy failures
US20090182494A1 (en) * 2008-01-15 2009-07-16 Honeywell International, Inc. Navigation system with apparatus for detecting accuracy failures
US20090182493A1 (en) * 2008-01-15 2009-07-16 Honeywell International, Inc. Navigation system with apparatus for detecting accuracy failures
US20090182495A1 (en) * 2008-01-15 2009-07-16 Honeywell International, Inc. Navigation system with apparatus for detecting accuracy failures
US20110210890A1 (en) * 2008-11-05 2011-09-01 Robert Andre Laine System for monitoring a satellite constellation of a positioning system
US8717234B2 (en) * 2008-11-05 2014-05-06 Astrium Sas System for monitoring a satellite constellation of a positioning system
CN101799524B (en) * 2009-07-10 2013-01-30 中国测绘科学研究院 Method for autonomously monitoring receiver integrity of global navigation satellite system
US20110254729A1 (en) * 2009-09-29 2011-10-20 Texas Instruments Incorporated Cross coupled positioning engine (pe) architecture for sensor integration in global navigation satellite system (gnss)
US9030356B2 (en) * 2009-09-29 2015-05-12 Texas Instruments Incorporated Positioning system receiver sensor system coupled with measurement data output
US8922426B1 (en) * 2010-02-26 2014-12-30 Google Inc. System for geo-location
US20130088387A1 (en) * 2011-10-07 2013-04-11 Electronics And Telecommunications Research Institute Apparatus and method for monitoring malfunctioning state of global positioning system (gps) satellite
US9664794B2 (en) * 2011-10-07 2017-05-30 Electronics And Telecommunications Research Institute Apparatus and method for monitoring malfunctioning state of global positioning system (GPS) satellite
US8706319B2 (en) 2011-12-16 2014-04-22 The Boeing Company Space positioning system
US9341718B2 (en) 2012-09-07 2016-05-17 Honeywell International Inc. Method and system for providing integrity for hybrid attitude and true heading
CN105487088A (en) * 2015-09-12 2016-04-13 北京大学 RAIM algorithm in satellite navigation system based on Kalman filtering
WO2018059532A1 (en) * 2016-09-30 2018-04-05 华为技术有限公司 Integrated navigational data calculation method of time-delay observation system, apparatus and navigational device
CN107884800A (en) * 2016-09-30 2018-04-06 华为技术有限公司 Observe the integrated navigation data calculation method, apparatus and navigation equipment of time lag system
CN107884800B (en) * 2016-09-30 2020-06-26 华为技术有限公司 Combined navigation data resolving method and device for observation time-lag system and navigation equipment
CN115032668A (en) * 2022-06-15 2022-09-09 北京讯腾智慧科技股份有限公司 Desensitization method and device for satellite navigation positioning reference station observation data
CN115032668B (en) * 2022-06-15 2023-02-21 北京讯腾智慧科技股份有限公司 Desensitization method and device for satellite navigation positioning reference station observation data

Similar Documents

Publication Publication Date Title
US6134484A (en) Method and apparatus for maintaining the integrity of spacecraft based time and position using GPS
US6166683A (en) System and method for high-integrity detection and correction of cycle slip in a carrier phase-related system
EP0904551B1 (en) Spoofing detection system for a satellite positioning system
US5596328A (en) Fail-safe/fail-operational differential GPS ground station system
US7219013B1 (en) Method and system for fault detection and exclusion for multi-sensor navigation systems
US5917445A (en) GPS multipath detection method and system
US6424914B1 (en) Fully-coupled vehicle positioning method and system thereof
EP0950193B1 (en) Gps receiver fault detection method and system
US20020109628A1 (en) Integrated inertial/gps navigation system
US5841399A (en) Fault detection and exclusion used in a global positioning system GPS receiver
CN110596732B (en) GBAS ionosphere anomaly detection method based on LMS adaptive filtering
RU2634693C2 (en) Validity control method and combination/consolidation device with multiple processing modules
US20130090858A1 (en) System for measuring coseismic movements or vibrations of structures based on global navigation satellite systems-gnss and/or pseudolites
Nikiforov et al. Application of statistical fault detection algorithms to navigation systems monitoring
Da et al. A new failure detection approach and its application to GPS autonomous integrity monitoring
Da et al. Failure detection and isolation structure for global positioning system autonomous integrity monitoring
Harvey Development of a precision pointing system using an integrated multi-sensor approach
Lee et al. Detection of multiple faults in single-frequency differential GPS measurements
KR100438396B1 (en) Realtime multipath detection method in Global navigation satellite system and positioning apparatus utilizing the same
Yoo et al. Performance comparison of GPS fault detection and isolation via pseudorange prediction model based test statistics
Young et al. Fault detection and exclusion using normalized solution separation methods
Li et al. GPS receiver failure detection method in high dynamic environment
Riba et al. Multi-satellite cycle-slip detection and exclusion using the noise subspace of residual dynamics
Zhang et al. A Study on Fault-Tolerant Integrated Navigation Systems
Hunzinger et al. CDGPS RAIM: Algorithm and protection radius calculation

Legal Events

Date Code Title Description
AS Assignment

Owner name: MOTOROLA, INC, ILLINOIS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:GEIER, GEORGE JEFFREY;HART, ROGER CHARLES;BRENCHLEY, GAYLIN DAVID;REEL/FRAME:010542/0837

Effective date: 20000127

STCF Information on status: patent grant

Free format text: PATENTED CASE

AS Assignment

Owner name: GENERAL DYNAMICS DECISION SYSTEMS, INC., ARIZONA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:MOTOROLA, INC.;REEL/FRAME:012435/0219

Effective date: 20010928

FPAY Fee payment

Year of fee payment: 4

AS Assignment

Owner name: GENERAL DYNAMICS C4 SYSTEMS, INC., VIRGINIA

Free format text: MERGER AND CHANGE OF NAME;ASSIGNOR:GENERAL DYNAMICS DECISION SYSTEMS, INC.;REEL/FRAME:016996/0372

Effective date: 20050101

FPAY Fee payment

Year of fee payment: 8

FPAY Fee payment

Year of fee payment: 12

AS Assignment

Owner name: GENERAL DYNAMICS MISSION SYSTEMS, INC., VIRGINIA

Free format text: MERGER;ASSIGNOR:GENERAL DYNAMICS ADVANCED INFORMATION SYSTEMS, INC.;REEL/FRAME:039406/0107

Effective date: 20151209

Owner name: GENERAL DYNAMICS ADVANCED INFORMATION SYSTEMS, INC

Free format text: MERGER;ASSIGNOR:GENERAL DYNAMICS C4 SYSTEMS, INC.;REEL/FRAME:039405/0964

Effective date: 20151209