WO2002046788A2 - Method and apparatus for determining an error estimate in a hybrid position determination system - Google Patents

Method and apparatus for determining an error estimate in a hybrid position determination system Download PDF

Info

Publication number
WO2002046788A2
WO2002046788A2 PCT/US2001/050896 US0150896W WO0246788A2 WO 2002046788 A2 WO2002046788 A2 WO 2002046788A2 US 0150896 W US0150896 W US 0150896W WO 0246788 A2 WO0246788 A2 WO 0246788A2
Authority
WO
WIPO (PCT)
Prior art keywords
terminal
error
position determination
sub
pseudo
Prior art date
Application number
PCT/US2001/050896
Other languages
French (fr)
Other versions
WO2002046788A3 (en
Inventor
Alkinoos Vayanos
Original Assignee
Qualcomm Incorporated
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
Priority to AU3973602A priority Critical patent/AU3973602A/en
Priority to EP01987533A priority patent/EP1328822B1/en
Priority to DE60139493T priority patent/DE60139493D1/en
Priority to AU2002239736A priority patent/AU2002239736B2/en
Priority to IL15563501A priority patent/IL155635A0/en
Priority to CA002426341A priority patent/CA2426341A1/en
Application filed by Qualcomm Incorporated filed Critical Qualcomm Incorporated
Priority to AT01987533T priority patent/ATE438868T1/en
Priority to KR1020037005790A priority patent/KR100829666B1/en
Priority to JP2002548471A priority patent/JP4414136B2/en
Publication of WO2002046788A2 publication Critical patent/WO2002046788A2/en
Publication of WO2002046788A3 publication Critical patent/WO2002046788A3/en
Priority to IL155635A priority patent/IL155635A/en

Links

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
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0205Details
    • G01S5/021Calibration, monitoring or correction
    • 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
    • G01S1/00Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith
    • G01S1/02Beacons or beacon systems transmitting signals having a characteristic or characteristics capable of being detected by non-directional receivers and defining directions, positions, or position lines fixed relatively to the beacon transmitters; Receivers co-operating therewith using radio waves
    • 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/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • 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
    • 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/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/48Determining position by combining or switching between position solutions derived from the satellite radio beacon positioning system and position solutions derived from a further system
    • 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
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0205Details
    • G01S5/0244Accuracy or reliability of position solution or of measurements contributing thereto
    • 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
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0252Radio frequency fingerprinting

Definitions

  • the current invention relates to position determination. More specifically, the present invention relates to a method and apparatus for improving error estimates of a position determination measurement made when determining the position of a device.
  • GPS Global Positioning System
  • GPS is a "constellation" of 24 well-spaced satellites that orbit the Earth.
  • the accuracy with which GPS can determine the position of a GPS receiver is anywhere from 100 to 10 meters for most receivers.
  • Each satellite within the constellation of GPS satellites transmits signals encoded with information. The information allows receivers on earth to measure the time of arrival of the received signals relative to an arbitrary point in time. This relative time of arrival measurement is generally referred to as a "pseudo- range" measurement.
  • GPS is owned and operated by the U.S. Department of Defense, but is available for general use around the world. Briefly, GPS includes 21 "regular" satellites and three spare satellites in orbit at 10,600 miles above the Earth. The satellites are spaced so that from any point on Earth, at least four satellites will be above the horizon. Each satellite contains a computer, an atomic clock, and a radio. With an understanding of its own orbit and the clock, each satellite continually broadcasts its changing position and time. Once a day, each satellite checks its own sense of time and position with ground stations and corrects the information as necessary. On the ground, each GPS receiver contains a computer that "triangulates" its own position by getting bearings from three satellites for a two dimensional solution.
  • the result is provided in the form of a geographic position.
  • This position is typically in the form of longitude and latitude.
  • the accuracy of the position determination is typically within 100 meters.
  • the receiver is also equipped with a display screen that shows a map, the position can be shown on the map.
  • the receiver/computer can figure out the altitude as well as the geographic position. If the receiver is moving, that receiver may also be able to calculate the speed and direction of travel of the receiver and give an estimated time of arrival to specified destinations.
  • GPS receivers Unfortunately, signals from GPS satellites are received at very low power levels due to the relatively large distances between the transmitting satellites and the receivers. , Therefore, minimal obstructions in the signal path that either block or disperse the signal make it impossible for receivers to receive the signals. For example, most GPS receivers have great difficulty receiving signals inside a building, under dense foliage, in urban settings in which tall buildings block much of the sky, etc. Accordingly, other techniques are used in place of, or to supplement, GPS. One such system is commonly referred to as a "hybrid position determination" system.
  • a hybrid position determination system includes a position determination terminal that includes both a GPS receiver and a communication system receiver.
  • the communication system receiver is a cellular telephone receiver.
  • a position determination beacon within the communication system communicates with the hybrid position determination terminal.
  • Signals from GPS satellites are received when available by the hybrid position determination terminal via the GPS receiver.
  • "Aiding information” is received from the position determination beacon by the hybrid position determination terminal via the communication system receiver.
  • the aiding information includes information that allows GPS satellite signals to be rapidly located in frequency and time.
  • the communication system signals can also be used to determine pseudo-ranges to base stations, one or more of which may be a position determination beacon. The pseudo- ranges to the base stations are used together with the pseudo-ranges to the satellites to calculate the position of the receiver.
  • base stations provide a time reference to the position determination receiver within the hybrid position determination terminal.
  • the time reference provided to the receiver by the communication system is GPS time.
  • the GPS time that is provided is offset by the amount of time required for the signal communicating the GPS time to propagate from the position determination beacon to the position determination receiver. This offset can be determined by measuring the propagation delay encountered by a signal that is transmitted on a "round trip" from the communication system receiver to the position determination beacon and back to the communication system receiver. The offset is then equal to one half of the total round trip delay (RTD).
  • RTD total round trip delay
  • Multipathing occurs when a signal takes an indirect path between the transmitter (i.e., the position determination beacon) and the receiver (i.e., the position determination terminal).
  • An indirect path is defined as a path that is longer than the shortest distance between the transmitter and receiver.
  • the word multipathing implies that more than one signal path will be traversed by the signal between the transmitter and receiver. However, for the purposes of this discussion, a signal would still be considered to be a multipath signal, even if the signal takes only one indirect path between the transmitter and the receiver.
  • Multipathing increases the amount of time required for the signal to traverse the distance between the position determination beacon and the position determination terminal. This increase is due to the longer distance traveled by the signal as a consequence of reflections off obstacles, such as buildings.
  • the increase in the amount of time required for the signal to arrive at the receiver results in an error in the pseudo- range measurement. This pseudo-range measurement error is then translated into an error in the position that is calculated from the pseudo-range measurements.
  • Multipathing can be a problem in GPS signals. However, it is easier to mitigate the effects of multipathing in GPS signals, since it is likely that the signal will still arrive at the position determination terminal via the direct path. That is, the signal between the GPS satellite and the position determination terminal is likely to take more than one path. However, one of those paths is likely to be the direct path. Accordingly, the direct path is assumed to be the one that is first to arrive. In addition, the direct path typically will have greater signal strength. In contrast, communication signals transmitted from a position determination beacon are more likely to take only indirect paths.
  • the following description discloses a method and apparatus for determining an estimate of the amount of error that is present in pseudo-range measurements made in a hybrid position determination system.
  • the presently disclosed method and apparatus allows correlations between a selected parameter and the error in a pseudo range measurement to be exploited.
  • a database is established in which the amount of error estimated for particular pseudo range measurements to a beacon is maintained.
  • Clusters are defined. Each cluster is associated with a range of values for the selected parameter. Pseudo range measurements are then associated with a particular cluster based upon the value of the selected parameter at the time (or proximate to the time) the pseudo range measurement was taken. As more estimates of the pseudo range measurements are made, the size of the clusters (i.e., the range of values of the selected parameter) can be reduced. Due to the correlation between the selected parameter and the errors in the pseudo range measurements, reducing the size of the clusters reduces the variance of the error estimates. The mean value of the error estimates is used to correct the error in future pseudo range measurements.
  • the position of a terminal that measures the pseudo range to a beacon is the selected parameter.
  • any other correlated parameter, such as the power level of the beacon signal may be the selected parameter.
  • the size of the cluster is initially relatively large, since the database will have relatively few error estimates in any particular geographic region. However, as the number of error estimates increases, the size of the clusters can be decreased, thus reducing the variance of the error estimates within the smaller clusters with respect to the larger clusters.
  • error estimates are made by first calculating what the pseudo range to a particular beacon should be. This calculation is made by determining (using a highly accurate first position determination sub-system) the current position of the terminal used to make the pseudo range measurements to the beacon.
  • the pseudo range measurements to a beacon can easily be calculated, assuming that the location of the beacon is known.
  • the pseudo range to the beacon from the terminal is then measured using a less accurate second position determination sub-system.
  • the difference between the pseudo range that is calculated based on the more accurate first position determination sub-system and the pseudo range measurement made by the less accurate second position determination sub-system is determined. This difference is assumed to be due to the error in the measurement made by the less accurate second position determination sub-system.
  • the database includes information that allows pseudo range measurements made by the less accurate second position determination sub-system to be corrected when the more accurate first position determination sub-system is not available.
  • the database is self generating, in that the information required in the database is taken during the operation of the terminal based upon the availability of the more accurate first position determination sub-system.
  • the more points the terminal can locate using the more accurate first position determination sub-system the smaller the clusters in the database.
  • the variance in the error estimates that are maintained in the database for each cluster will be reduced.
  • the presently disclosed method and apparatus could be used with position determination systems other than hybrid position determination systems if there is some other means by which to determine the position of the terminal and that means is available at some times or locations, but not at others. In that case, the position of the terminal would be used as a reference to determine the amount of error in the pseudo range measurements in the same way as the position determination from the more accurate sub-system described above.
  • an iterative approach is used when the selected parameter is the location of the terminal. Assuming that the more accurate first position determination sub-system is not available, and that a sufficient number of initial error estimates have been made, the iterative approach uses corrected measurements of the pseudo range based upon a relatively large cluster to determine the position of the terminal. Once the location of the terminal has been determined in this way, the corrections to the pseudo ranges can be recalculated based upon a much smaller cluster, assuming that a statistically valid number of error estimates has been made for the smaller cluster.
  • FIG. 1 is an illustration of a hybrid position determination system in accordance with the disclosed method and apparatus.
  • FIG. 2 illustrates the disclosed method in which there is a relation between the approximate position of a terminal and the expected error in pseudo range measurements made by the terminal.
  • Fig. 3 is an illustration of an entire region shown as a single cluster.
  • Fig. 4 is an illustration of four smaller clusters that have been subdivided from one larger cluster.
  • Fig. 5 is an illustration of an example in which some of the clusters of Fig. 4 have been combined to form a new cluster.
  • Fig. 6 illustrates the process for building an MSD in the case in which there is no assumption regarding the errors in the satellite measurements.
  • FIG. 7 is a simplified block diagram of one embodiment of an iterative process for correcting errors in position determinations.
  • FIG. 8 is a simplified block diagram of one embodiment of a terminal used to implement the presently disclosed method.
  • FIG. 1 is an illustration of a hybrid position determination system 100.
  • a hybrid position determination system is defined as a system in which the position of a terminal 105 can be determined from a first position determination sub-system alone or in combination with a second sub-system.
  • the first position determination sub-system is preferably capable of determining the position of the terminal 105 "independent" of the second sub-system, the present context, the phrase "independent of the second sub-system" means that the first sub-system may be used to determine the position of a terminal at times when some or all of the components of the second sub-system are not available.
  • components of the second sub-system may be a necessary part of the first sub-system.
  • the second sub-system may or may not be capable of performing a position determination without the use of components of, or information from, the first subsystem.
  • At least one "parameter” is related to the magnitude of the errors in position determination measurements of one of the position determination sub-systems.
  • a parameter is any variable that is related to the magnitude of the errors in position determination measurements of one of the position determination sub-systems.
  • a parameter for a position determination terminal 105 can be: (1) the location of the position determination terminal 105, (2) the amount of power in a position determination signal received by the position determination terminal 105, (3) the number of base stations that are present in the geographic region in which the position determination terminal 105 is presently located, (4) the type of buildings in the geographic region in which the position determination terminal 105 is presently located, (5) the density of the structures in the geographic region in which the position determination terminal 105 is presently located, etc.
  • a parameter which is related to the magnitude of the errors in measurements made by the second sub-system is preferably unrelated (and alternatively, weakly related) to the magnitude of the errors in measurements made by the first position determination sub-system.
  • each of the two sub-systems includes a set of transmitters.
  • the transmitters of each sub-system transmit position determination signals.
  • four satellites 101, three position determination beacons 103, and one position determination terminal 105 are shown.
  • the satellites 101 are a first of the two types of transmitters of position determination signals and are associated with the first sub-system.
  • the satellites 101 provide signals that can be received by the terminal 105.
  • the received signals enable the terminal to determine the position of the terminal 105 Without the need to receive position determination signals from any of the beacons 103, assuming that signals from a sufficient number of satellites signals are received and that the information received on those signals can be decoded. Nonetheless, in some embodiments of the disclosed method and apparatus, the terminal 105 may require communication with a beacon 103 for acquisition and processing assistance. Whether acquisition assistance is desirable or not, once a sufficient number of satellites 101 are acquired, the information provided by the satellites 101 is sufficient to calculate the position of the terminal 105.
  • the beacons 103 are a part of a second of the two position determination subsystems. Like the satellites 101, the beacons 103 provide signals that can be received by the terminal 105. In accordance with one embodiment of the disclosed method and apparatus, the signals enable the terminal 105 to determine its own position without the need to receive position determination signals from any of the satellites 101. However, in an alternative embodiment, the second sub-system requires at least one transmitter from the first sub-system to provide information in order to determine the position of the terminal 105. Once the information is received from each of the required transmitters, a determination of the position can be performed by the terminal 105, either alone or with acquisition and processing assistance from one of the beacons 103.
  • the position of the terminal 105 can be determined by a device (not explicitly shown) that is remote from the terminal 105.
  • the terminal 105 transmits any information necessary to determine the position of the terminal 105 to the remote device.
  • the remote device is one of the beacons 103 or a sub-system 104 within one of the beacons 103.
  • the remote device is shown in FIG. 1 as "optional" to highlight the fact that the remote device may be located either in the beacon 103 or external to the beacon 103. It will be understood by those skilled in the art that there may be more or less satellites 101 or beacons 103 than are shown in FIG. 1, and that there may be more than one terminal 105.
  • the satellites 101 are GPS satellites.
  • the satellites 101 may be any other type of transmitter that enables a relatively independent and relatively accurate position determination to be made.
  • a terrestrial position determination system such as LORAIN
  • LORAIN that is capable of providing an accurate position determination may be used instead of satellites.
  • the beacons 103 are cellular base stations, including (1) base station transceiver sub-systems (BTSs), (2) base station controllers (BSCs), and (3) position determination equipment (PDEs).
  • BTSs base station transceiver sub-systems
  • BSCs base station controllers
  • PDEs position determination equipment
  • the beacons 103 may include a subset of these components.
  • the beacons 103 may be any other transmitters capable of transmitting position determination signals to enable the position of a terminal 105 to be determined. It should also be understood by those skilled in the ait that in a beacon 103 that includes a PDE, the PDE would likely serve as the remote device 104 shown in FIG. 1.
  • information received from the first type of transmitters should enable a substantially more accurate position determination to be calculated than can be determined from the information received from the second type of transmitters.
  • the accuracy of the position determination made using measurements to the satellites is typically greater than the accuracy of the position determinations made using measurements to the cellular base stations.
  • the presently disclosed method and apparatus will be described in the context of a hybrid position determination system in which the satellites are GPS satellites and the beacon is a code division multiple access (CDMA) cellular base station.
  • CDMA code division multiple access
  • the present invention does not rely upon any of the unique features of these systems. Accordingly, the invention may be implemented using any other systems for determining position.
  • the satellites are GPS satellites and the beacons are CDMA cellular base stations (which may or may not include a BSC and PDE)
  • "pseudo range" measurements psi - Ps are made with respect to the satellites 101 and similarly, pseudo range measurement pt > i - pb3 are made with respect to the beacons 103.
  • Pseudo range measurements represent the relative distance between the receiving terminal 105 and the sources of position determination signals. It should be noted that the particular manner in which position determinations are made is not relevant to the presently disclosed method and apparatus. However, a discussion of pseudo range based position determination techniques is provided as an example of one embodiment of the disclosed method and apparatus. Other means for measuring pseudo ranges are well known and could be used to perform the disclosed method and apparatus. Furthermore, other means for determining position without the need for pseudo range measurements are well known and could be used to perform the disclosed method and apparatus.
  • Pseudo range measurements typically represent a time difference between an arbitrary point in time and the time at which a signal arrives at a position determination terminal 105.
  • pseudo ranges are typically presented in units of meters. The amount of time can be translated into a distance in meters by multiplying the time difference by the speed of light in meters per second.
  • the position determination terminal 103 is a CDMA base station, the position determination terminal notes when the beginning of a series of bits that make up a spreading code for the CDMA signal has been received with respect to a prior arbitrary point in time. This series of bits is commonly referred to as a "pseudo-random noise (PN) spreading code”.
  • PN pseudo-random noise
  • the signals transmitted by each of the base stations within a CDMA communications system are encoded with the same PN spreading code.
  • these offsets must be taken into account before comparing the relative timing of the signals received from each of the beacons 103 with respect to one another.
  • these offsets are easily determined and can be subtracted out. In fact, these offsets are typically used to identify the particular base station from which the signal originated.
  • the duration of the clock (assumed for this example to be one millisecond) must be such that any ambiguity can be resolved. That is, the duration of the clock must be such that the values of the clock that are associated with the received signals should uniquely identify only one relationship between the received signals. Accordingly, the duration of the clock must be much greater than twice the difference between the time at which the signals will be received from each beacon 103, taking into account the offsets in the timing of the codes transmitted by each beacon 103. It should be noted that making pseudo range measurements is very well known in the art.
  • the pseudo range measurements can be used to determine the position of the terminal 105.
  • the pseudo range measurements and the position of the beacons 103 are applied to a well-known process commonly referred to as the "Least Mean Squares" (LMS) process.
  • LMS LMS
  • the inverse is true. That is, if the position of the terminal 105 and three beacons 103 are known, the distance between the terminal 105 and each of the three beacons 103 can be used to determine the pseudo range values from the position of the terminal 105 to the beacons 103.
  • the presently disclosed method and apparatus takes advantage of the fact that a hybrid position determination system will frequently have two independently derived sets of pseudo range measurements.
  • One set of measurements is typically going to be more accurate than the other. Therefore, if the more accurate set of measurements is sufficient to determine the position of the terminal 105, then the more accurate set of measurements can be used to determine the error in the less accurate set of measurements. Furthermore, it has been determined that there is a predictable relationship between particular parameters (such as position of the terminal 105) and the amount of error in the less accurate set of measurements. Therefore, by nowing both the value of the parameter and the amount of error associated with that parameter, an estimate as to the amount of error in the less accurate measurements can be made. For example, by knowing the approximate position of the terminal 105 and the relationship between the position and the amount of error, the amount of error in the less accurate measurement can be estimated.
  • a "Measurement Statistics Database” is generated in accordance with one embodiment of the presently disclosed method and apparatus. It should be noted that means other than a database may be used for associating the values of a parameter with the amount of error in the measurements made by the less accurate position determination sub-system.
  • the MSD provides an effective means by which a parameter (such as the position of terminal 105) can be associated with a correction factor to be applied to the less accurate pseudo range measurements. Accordingly, when the value of the parameter is known, the associated correction factor can be determined from the MSD and applied to the pseudo range measurement.
  • the parameter at issue is the position of the terminal 105 as determined from essentially uncorrected pseudo range measurements made by the terminal 105.
  • the pseudo range measurements are completely uncorrected.
  • the pseudo range measurements may be corrected for time offsets due to differences between the time at which the same code is transmitted from different beacons.
  • the parameter at issue may be any one of a large number of other parameters. It should be noted that it is advantageous for the terminal 105 or the beacon 103 to be able to directly measure the value of the parameter at issue. Direct measurement of the parameter by the terminal 105 or the beacon 103 allows the disclosed method and apparatus to operate without input from an external source. However, it is possible for the value of the parameter to be input by an external source, such as the terminal operator or signals received from a remote source.
  • Possible parameters include, but are not limited to: (1) the amount of power with which the position determination signals have been received by the terminal 105; (2) the type of buildings that are proximate to the terminal 105; (3) the amount of urban development in the general vicinity of the terminal 105; (4) the distance from the beacons 103 to the terminal 105 that will be receiving the position determination signals; and (5) the shape of the correlation peak determined from the correlation of the received signal with a known PN spreading code. It should be clear that this list merely provides a small sample of the many types of parameters that could be used to predict the amount of error that is likely to be present in pseudo range measurements. The present invention should not be limited to the types of parameters that are listed herein.
  • the scope of the present disclosure should be considered to include using any parameter that is correlated with a first set of pseudo range measurements and not (or alternatively, to a lesser degree) with a second set of pseudo range measurements to predict the amount of error in the first set of measurements.
  • FIG. 2 particularly relates to the example in which there is a relation between the approximate position of a terminal 105 and the expected error in pseudo range measurements made by the terminal 105.
  • the process of generating the MSD requires that the more accurate position determination capability be available. If not available, then the process of generating an MSD can go no further until the more accurate position determination capability becomes available. Nonetheless, it should be understood that in this case, the MSD, if already sufficiently generated, can be used to estimate the error in the less accurate sub-system.
  • the terminal 105 has the ability to determine position based upon either pseudo range measurements to satellites 101 or pseudo range measurements to beacons 103.
  • the pseudo range measurements to the satellites 101 tend to be more accurate. Therefore, in such an embodiment, the position of the terminal 105 is determined using pseudo ranges to the satellites 101, if possible.
  • the beacon pseudo ranges are used only to supplement the satellite pseudo range measurements when an insufficient number of satellites 101 are available. However, similar to the case in which the beacon pseudo ranges can be used with a satellite to determine the position of the terminal 105, using satellite pseudo ranges is preferred over using beacon pseudo ranges.
  • the position of the terminal 105 is determined from pseudo ranges to the satellites 101 in essentially the same way as described above with respect to the determination of pseudo ranges from the beacons 103. Greater accuracy in the pseudo range measurements to the satellites 101 than can typically be had from the beacons 103 comes from several factors, including the fact that the satellites 101 are overhead.
  • the probability that the signals received from the satellites 101 will reach the terminal 105 on a direct path is greater than the probability that the signals from the beacons 103 will reach the terminal 105 directly.
  • the failure of the signals from the beacons 103 to reach the terminal 105 directly adds to the distance along the path propagated by the signals from the beacon 103 to the terminal 105.
  • the additional distance causes an error in calculating the position of the terminal 105 when using the pseudo ranges from the beacons 103.
  • the pseudo ranges to the beacons 103 can be computed (STEP 207). It can be seen that the expected pseudo range to a particular beacon 103 can be easily determined from knowledge of the position of the terminal 105 and from knowledge of the position of the beacon 103.
  • the pseudo ranges are provided to the beacon 103, which determines the expected pseudo ranges for a particular terminal 105 to a particular beacon 103.
  • the terminal 105 does the determination and calculation, then the presently disclosed method and apparatus assumes that the terminal 105 has access to knowledge of the position of the beacons 103.
  • the base stations provide the terminal 105 with information regarding their own position.
  • the terminal 105 maintains a database that provides the terminal 105 with the position of the beacon 103 based upon an identification indication received as part of the transmissions from the beacon 103 to the terminal 105.
  • the terminal 105 sends information to one or more of the beacons 103 to process.
  • the beacon 103 knows or has access to the location of the beacons in the area.
  • the terminal 105 measures the pseudo ranges to each beacon 103 (STEP 209). Each measured pseudo range is associated with one beacon 103 and is compared to a pseudo range associated with that beacon 103 as computed in STEP 207 (STEP 211). It is assumed that the difference between a computed and measured pseudo range is an error in the measured pseudo range. It should be clear that the pseudo ranges to the both beacons 103 and the satellites 101 can be measured and sent to the beacon 103 for processing prior to the calculation in STEP 207 of the beacon pseudo ranges.
  • a cluster in which the terminal 105 is presently located is determined (STEP 213).
  • a cluster can be defined in a number of ways. In accordance with the preferred embodiment, a cluster is defined as a set of contiguous values of the parameter. For example, if the parameter is the position of the terminal 105, then the cluster would be a contiguous geographic region.
  • the set of contiguous parameter values is preferably large enough to include several error measurements. That is, preferably enough error estimates have been made within the cluster to calculate with reasonable accuracy, the mean of the pseudo range errors made by the less accurate sub- system at any position within the cluster.
  • the size of the cluster is increased to include more error estimates.
  • the cluster will be considered to be immature until a statistically sufficient number of error estimates are made.
  • the cluster is divided into two. Each cluster should have a sufficient number of error estimates to allow the calculated mean of the error associated with that cluster to reflect with reasonable accuracy the mean that would be calculated from an essentially infinite number of error estimates taken from within the cluster. If the mean of the errors in contiguous clusters are the same, then those clusters may be combined, even though there are sufficient number of error estimates to support more than one cluster. However, in accordance with one embodiment of the disclosed method and apparatus, the cluster may be divided differently if such a division would result in two clusters, each having a sufficient number of error estimates and different error estimate means.
  • the desired result is to associate each cluster with a value, or set of values, of a parameter. If the value of the parameter detected by the terminal 105 is associated with the cluster, then the terminal 105 is assumed to be in the associated cluster. For example, if the parameter is the position of the terminal 105, then the cluster would be associated with a range of geographical positions. If the terminal 105 is detected to be located within a cluster, then the terminal 105 is considered to be within that cluster for the purposes of the disclosed method and apparatus. Alternatively, assume that the parameter of interest is the power level of the received signal. If the terminal 105 detects that the received signals are within the range of power levels associated with a particular cluster, then that terminal 105 would be considered to be within the associated cluster.
  • the parameter is selected such that there is a correlation between the parameter and the amount of error in the pseudo range measurements made by the less accurate position determination sub-system. Accordingly, each value of the parameter is associated with a particular value in the error in the position determinations made by the less accurate position determination sub-system. That is, whenever the terminal 105 notes that the parameter has a particular value (or is within a particular range of values), the error in the determination of the position made by the less accurate position determination sub-system will be a particular value (or within a particular range of values).
  • each position of the terminal 105 would be associated with a particular amount of error in the position determination measurement made by the terminal 105 using the less accurate position determination sub-system. Accordingly, to select position as an appropriate parameter, each position of the terminal 105 must be related to a particular amount of error in the position determination measurement made by the less accurate position determination sub-system. Accordingly, for a terminal at position X the amount of error would be Y. This error value would be essentially constant and predictable, within a range of tolerable uncertainty.
  • FIGs. 3 and 4 A region 300 is shown in which Eleven positions 301 are identified. Pseudo range error estimates were made at each of these positions 301. The error estimates are based upon the difference between pseudo range measurements taken at those locations 301 and pseudo ranges calculated from a knowledge of the position based on the more accurate position determination sub-system.
  • FIG. 3 is an illustration of the entire region 300 shown as a single cluster. This is so because initially there are only relatively few (eleven) error estimates made in the region 300. It should be noted that the particular number of estimates that is required to constitute a sufficient number depends on the particular application. Therefore, the number may be much greater or less than eleven depending upon the application. It should also be noted that while the beacons 103 are shown in FIG. 3 to be outside the region 300, the positional relationship between the beacons 103 and the region is not relevant to the present invention other than noting that it should be possible to receive the signals from the beacons 103 from within the region 300.
  • FIG. 4 is an illustration of four smaller clusters 401, 403, 405, and 407 which have been subdivided from one larger cluster 300.
  • FIG. 5 is an illustration of an example in which clusters 401, 403, and 407 of FIG. 4 have been combined to form yet another new cluster 501 including positions 500 from which error estimates have been made. Accordingly, clusters having relatively complex shapes can be formed using clusters having relatively simple shapes. Having more error estimates taken at different positions within a cluster provides more opportunities to reshape the clusters.
  • the measurements taken within a cluster will typically have differing amounts of error.
  • the set of error values form an array associated with a particular beacon 103 within a particular cluster.
  • the values within the array will vary about a mean value within some range.
  • one way to characterize the size of the range in this variation is to calculate the "variance" of the array.
  • the size of the range can be characterized by the "deviation" of the array.
  • the boundaries of a cluster can be determined based upon the value of the error estimates of the pseudo ranges measured. In this case, the variation in the error estimates included in the array is minimized.
  • the area within a cluster is defined such that variations in the error estimates that have been made by position determination terminals 105 within the cluster are within a predetermined range.
  • the shape of the cluster is altered to remove the offending position. This results in greater accuracy in determining the correction to be applied to measurements taken by the less accurate position determination sub-system. This is particularly valuable when it is necessary to rely on at least one measurement from the less accurate sub-system in order to determine the position of the terminal 105.
  • the new variance and mean replace the old (STEP 217).
  • the error itself would be the mean and the variance would be equal to zero.
  • Each newer error value associated with that cluster would then be placed into the array of older values.
  • the entire array is maintained within the MSD.
  • only the new variance, mean, and number of elements in the array are maintained. It should be noted that in order to update the variance and mean of the error associated with that cluster the entire array must be available. Alternatively, the mean, variance, and the number of elements in the array must be known.
  • the MSD is dynamically generated.
  • the values of all of the elements in the array are stored in the MSD.
  • the values for the mean and variance can be calculated as needed from the stored array elements, rather than being calculated after each measurement is made.
  • the mean and variance can be weighted by giving more significance to certain elements in the array.
  • the elements that are given more weight are those elements that are known to be more reliable for one reason or another. For example, if a terminal 105 makes a pseudo range measurement based on a relatively weak received signal, the weight of such a pseudo range measurement might be reduced with respect to other measurements made from stronger signals.
  • the satellite errors i.e., the errors in the position determined using the satellite pseudo range measurements
  • the satellite errors will typically not bias the mean of the errors in a cluster.
  • no such assumption is made.
  • FIG. 6 illustrates the process for building an MSD in the case in which there is no assumption regarding the errors in the satellite measurements, as described below.
  • Pseudo range measurements to the satellites 101 are made by the terminal 105 (STEP 601).
  • the terminal notes the variance in the satellite measurements (STEP 603). That is, since the measurements will have a certain amount of error, the measurements will vary about a mean within some range. This variation is characterized by calculating the variance in the error in the satellite measurements. Alternatively, the errors can be characterized by their deviation.
  • the variance (or deviation) and the mean in the errors in the satellite measurements are typically known, hi general, the variance and means are provided in a covariance matrix associated with the each satellite 101, as is well known in the art.
  • the satellite measurements and the error statistics associated with the satellite measurements are applied to a well- known LMS process (STEP 604).
  • the satellite geometry is also taken into account by the LMS process.
  • the result of applying the satellite measurements, the error statistics associated with those measurements, and the satellite geometry to the LMS process is a position determination solution.
  • the position determination solution provides the position of the terminal 105 from which the pseudo range measurements were taken.
  • the LMS process results in a solution covariance matrix that indicates the error statistics regarding the amount of error in the solution.
  • solution covariance matrices are well know in the art. The means are assumed to be zero due to the fact that there is assumed to be no multipathing errors in the signals received from the satellites.
  • the solution is presented in the form x, y, z, and b, where x, y, and z are the Cartesian coordinates of the terminal 105 with respect to a selected reference, such as the center of the earth, and b is the bias of the clock in the terminal 105 with respect to true GPS time.
  • x, y, and z are the Cartesian coordinates of the terminal 105 with respect to a selected reference, such as the center of the earth, and b is the bias of the clock in the terminal 105 with respect to true GPS time.
  • a beacon covariance matrix can be calculated from the solution covariance matrix (STEP 610).
  • the beacon covariance matrix characterizes the amount of error in each of the calculated error estimates based upon the amount of error in the satellite pseudo range measurements, as determined from the satellite covariance matrix. That is, estimates of the error in the beacon pseudo range measurements are based on the assumption that the calculated pseudo ranges are perfectly accurate. However, the statistics of the error is known for the satellite measurements used to derive the calculated pseudo ranges. Therefore, there is an error in the estimates of the amount of error in the beacon pseudo ranges.
  • the statistics of the error in the error estimates can be determined from the solution covariance matrix, which is in turn calculated from the satellite covariance matrix.
  • the MSD is generated such that each entry to the MSD is associated with a cluster, a group of beacons, the mean of the array of pseudo range measurement error estimates to each beacon, and the beacon covariance matrix that characterizes the error in each measurement of the pseudo range due to the errors introduced by multipathing (STEP 612).
  • each entry includes identification of an associated cluster, identification of an associated beacon, an array of pseudo range measurement error estimates associated with the beacon, and beacon covariance matrix for the associated beacon. It should be noted that it may not be possible to determine the covariance between beacons associated with different entries.
  • the error estimates of the array are weighted according to the reliability of each such estimate.
  • these values can be used. In particular, they are used to estimate the error in pseudo ranges to the particular beacons 103 when the pseudo range measurements are being made by a terminal 105 within the cluster. The estimates of the error in the beacon pseudo range measurements can then be used to correct the measured beacon pseudo ranges. This is particularly useful when there are too few satellite pseudo range measurements. Accordingly, if there are not enough satellites available to make a position determination based on satellite measurements alone, then corrected pseudo ranges to beacons 103 can be used. It should be understood that the presently disclosed method and apparatus assumes that the MSD will be generated and stored within the terminal 105.
  • the MSD may be generated and stored within a beacon or other component remote from the terminal 105.
  • the terminal 105 may be responsible only for taking pseudo range measurements. These measurements are then transmitted to the beacon 103.
  • the beacon 103 then processes the pseudo range measurements as described above in a device 104 within the beacon 103.
  • the device 104 is remote from the beacon 103 and the beacon 103 communicates the pseudo range measurements to the remote device 104.
  • the remote device 104 then processes the pseudo range measurements. It should be understood that the processing may be distributed such that some of the processing is performed in one device and other processing is performed in other devices. Nonetheless, without regard for which device is responsible for performing which process, the disclosed method will be performed essentially as defined above.
  • the present method and apparatus provides a means by which a position determination from a less accurate position determination sub-system can be corrected in a single correction step.
  • an iterative process can be performed as shown in FIG. 7 and described below.
  • the position of a terminal 105 is determined using a less accurate subsystem (STEP 701). That determination is then corrected based on a correction factor associated with a first relatively large cluster (STEP 703).
  • the first cluster is intentionally made to be relatively large to account for the fact that the error that is present in the less accurate sub-system's initial estimate of the position may be great enough to make it difficult to accurately determine which cluster the terminal 105 is in if the cluster is made to small.
  • the position of the terminal 105 will be known to a greater accuracy. Smaller clusters are then used to determine a more precise correction factor (STEP 705), since the improvement in accuracy provided by the first correction will make it possible to determine in which smaller cluster the terminal 105 current resides.
  • the presently disclosed method and apparatus could be used with position determination systems other than hybrid position determination systems if there is some other means by which to determine the position of the terminal is available at some times or locations, but not at others.
  • the position of the terminal would be used as a reference to determine the amount of error in the pseudo range measurements in the same way that the position determination from the more accurate sub-system described above is used.
  • FIG. 8 is a simplified block diagram of one embodiment of a terminal 105 used to implement the presently disclosed method.
  • the terminal 105 includes a receiver 501, a decoder 503, a correlator 505, a reference code generator 507, a processor 509, a clock 511, and a memory 513.
  • the receiver 501 typically includes a conventional radio frequency front end section, including antenna, down converter, filters, amplifiers, etc., all of which are well known and for simplicity sake are not shown. Signals, including position determination signals, are received by the receiver 510.
  • the received signals are down-converted, filtered, and amplified (not necessarily in that order) as required for input to the decoder 503. It will be understood by those skilled in the art that other processing not disclosed herein may occur in the receiver 501, but is essentially immaterial to the novel aspects of the presently disclosed method and apparatus.
  • the output from the receiver 501 is coupled to the decoder 503.
  • the decoder 503 decodes the intelligence from the received signal. When position determination signals are received and decoded, the output from the decoder will be a code sequence.
  • the decoded intelligence is coupled to the correlator 505.
  • the correlator 505 is also coupled to the reference code generator 507.
  • the reference code generator 507 provides a reference code sequence to the correlator 505 to allow the correlator 505 to identify a predetermined point in the received code sequence.
  • the correlator 505 outputs a signal to the processor 509 indicating when the reference code sequence and the received code sequence are synchronized.
  • the reference code generator 507 is also coupled to the processor 509. Once the reference code sequence and the received code sequence are synchronized, the output from the reference code generator 507 to the processor 509 allows the processor 509 to determine the timing of the received code sequence.
  • the clock 511 is also coupled to the processor 509.
  • the clock 511 may be either free running or synchronized to an external clock reference. If the clock is synchronized to an external clock reference, then the timing of the received code sequence can be determined with respect to the external clock reference.
  • the processor 509 is also coupled to the receiver 501.
  • the receiver 501 provides the processor 509 with information as to which frequency the receiver 501 is be tuned.
  • the processor 509 receives information from the timing of the received code sequence that, together with the frequency on which the signal was received, is indicative of the particular transmitter 101, 103 from which the code sequence originated. It should be noted that multiple correlators 505 may be used to concurrently correlate multiple signals.
  • the processor 509 takes all of the information that is provided to it and determines the relative timing (i.e., the pseudo ranges) of all of the signals received. In addition, the processor 509 performs all of the pseudo range processing functions described above.
  • a sector coverage map is a map that shows which geographic regions are being served by a particular sector of a beacon 103. Beacons 103 are said to be sectorized (i.e., serve a plurality of sectors) when separate antennas are used to direct the signals transmitted by the server to a particular geographic sector and to receive signals from that particular geographic sector. Typically, the sectors at issue are assumed to be essentially pie shaped wedges emanating out from the beacon 103.
  • a sector coverage map allows both a terminal 105 and a beacon 103 to determine which sector of the beacon 103 is servicing the terminal 105.
  • a sector coverage map is generated in accordance with the presently disclosed method and apparatus by associating (1) the sector in which the beacon 103 receives signals from the terminal 105, (2) the position of the terminal 105 as determined by the more accurate position determination sub-system, and (3) the position of the terminal 105 as determined by the less accurate position determination sub-system. Accordingly, a sector coverage map generated as indicated above can be used together with knowledge of the position of the terminal, as determined by the less accurate subsystem, to more reliably attain the sector of the beacon 105 that is servicing the terminal. This is particularly useful when more than one sector receives signals from a terminal 105. This is even more useful when more than one sector is reporting service to the terminal 105 and the position determined by the less accurate position determination sub-system determines that the terminal 105 is not in any of the sectors that are reporting that they are serving the terminal 105.

Abstract

A method and apparatus that allows correlations between a selected parameter and the error in a pseudo range measurement to be exploited. A database is established in which the amount of error estimated for particular pseudo range measurements to a beacon is maintained. Clusters are defined. Each cluster is associated with a range of values for the selected parameter. Pseudo range measurements are then associated with a particular cluster based upon the value of the selected parameter at the time (or proximate to the time) the pseudo range measurement was taken. As more estimates of the pseudo range measurements are made, the size of the clusters (i.e., the range of values of the selected parameter) can be reduced. Due to the correlation between the selected parameter and the errors in the pseudo range measurements, reducing the size of the clusters reduces the variance of the error estimates. The mean value of the error estimates is used to correct the error in future pseudo range measurements.

Description

METHOD AND APPARATUS FOR DETERMINING AN ERROR ESTIMATE IN A HYBRID POSITION DETERMINATION
SYSTEM
BACKGROUND OF THE INVENTION
I. Field of the Invention
The current invention relates to position determination. More specifically, the present invention relates to a method and apparatus for improving error estimates of a position determination measurement made when determining the position of a device.
II. Description of the Related Art
It has always been desirable for man to know his geographic location. Devices ranging from compasses, maps, sextants, surveying equipment, etc. have been used to determine a person's location. Today, we enjoy the benefits of a system of satellites that orbit about the earth and that provide information to receivers on earth. Each such receiver can use the information provided by the satellites to determine its position. One such system is the well-known Global Positioning System (GPS). GPS is a "constellation" of 24 well-spaced satellites that orbit the Earth. The accuracy with which GPS can determine the position of a GPS receiver is anywhere from 100 to 10 meters for most receivers. Each satellite within the constellation of GPS satellites transmits signals encoded with information. The information allows receivers on earth to measure the time of arrival of the received signals relative to an arbitrary point in time. This relative time of arrival measurement is generally referred to as a "pseudo- range" measurement.
GPS is owned and operated by the U.S. Department of Defense, but is available for general use around the world. Briefly, GPS includes 21 "regular" satellites and three spare satellites in orbit at 10,600 miles above the Earth. The satellites are spaced so that from any point on Earth, at least four satellites will be above the horizon. Each satellite contains a computer, an atomic clock, and a radio. With an understanding of its own orbit and the clock, each satellite continually broadcasts its changing position and time. Once a day, each satellite checks its own sense of time and position with ground stations and corrects the information as necessary. On the ground, each GPS receiver contains a computer that "triangulates" its own position by getting bearings from three satellites for a two dimensional solution. The result is provided in the form of a geographic position. This position is typically in the form of longitude and latitude. The accuracy of the position determination is typically within 100 meters. If the receiver is also equipped with a display screen that shows a map, the position can be shown on the map. If a fourth satellite can be received, the receiver/computer can figure out the altitude as well as the geographic position. If the receiver is moving, that receiver may also be able to calculate the speed and direction of travel of the receiver and give an estimated time of arrival to specified destinations.
Unfortunately, signals from GPS satellites are received at very low power levels due to the relatively large distances between the transmitting satellites and the receivers. , Therefore, minimal obstructions in the signal path that either block or disperse the signal make it impossible for receivers to receive the signals. For example, most GPS receivers have great difficulty receiving signals inside a building, under dense foliage, in urban settings in which tall buildings block much of the sky, etc. Accordingly, other techniques are used in place of, or to supplement, GPS. One such system is commonly referred to as a "hybrid position determination" system.
A hybrid position determination system includes a position determination terminal that includes both a GPS receiver and a communication system receiver. In one example of such a hybrid position determination system, the communication system receiver is a cellular telephone receiver. A position determination beacon within the communication system communicates with the hybrid position determination terminal.
Signals from GPS satellites are received when available by the hybrid position determination terminal via the GPS receiver. "Aiding information" is received from the position determination beacon by the hybrid position determination terminal via the communication system receiver. The aiding information includes information that allows GPS satellite signals to be rapidly located in frequency and time. In addition, the communication system signals can also be used to determine pseudo-ranges to base stations, one or more of which may be a position determination beacon. The pseudo- ranges to the base stations are used together with the pseudo-ranges to the satellites to calculate the position of the receiver.
In addition, base stations provide a time reference to the position determination receiver within the hybrid position determination terminal. In one particular hybrid system, the time reference provided to the receiver by the communication system is GPS time. However, the GPS time that is provided is offset by the amount of time required for the signal communicating the GPS time to propagate from the position determination beacon to the position determination receiver. This offset can be determined by measuring the propagation delay encountered by a signal that is transmitted on a "round trip" from the communication system receiver to the position determination beacon and back to the communication system receiver. The offset is then equal to one half of the total round trip delay (RTD). However, it should be noted that there is a delay that is added to the RTD by the internal delays associated with the reception and retransmission of the signal at the position determination beacon. Therefore, in order to get an accurate GPS time transfer from the position determination beacon to the position determination terminal, these internal delays must be determined and subtracted from the measured RTD. This is often referred to as "calibrating" the position determination beacon. Calibrating the position determination beacon requires measuring the amount of delay internal to the position determination beacon. Calibrating the position determination beacons is a time consuming and difficult task. Accordingly, it would be advantageous to provide a method and apparatus that would allow the position of a hybrid position determination terminal to be determined without requiring calibration of position determination beacons.
Even after having calibrated each of the position determination beacons within a communication system, the accuracy of the pseudo-range measurements that are made between the position determination terminal and the position determination beacon are not necessarily accurate. This is due to a phenomenon known as "multipathing". Multipathing occurs when a signal takes an indirect path between the transmitter (i.e., the position determination beacon) and the receiver (i.e., the position determination terminal). An indirect path is defined as a path that is longer than the shortest distance between the transmitter and receiver. The word multipathing implies that more than one signal path will be traversed by the signal between the transmitter and receiver. However, for the purposes of this discussion, a signal would still be considered to be a multipath signal, even if the signal takes only one indirect path between the transmitter and the receiver.
Multipathing increases the amount of time required for the signal to traverse the distance between the position determination beacon and the position determination terminal. This increase is due to the longer distance traveled by the signal as a consequence of reflections off obstacles, such as buildings. The increase in the amount of time required for the signal to arrive at the receiver results in an error in the pseudo- range measurement. This pseudo-range measurement error is then translated into an error in the position that is calculated from the pseudo-range measurements.
Multipathing can be a problem in GPS signals. However, it is easier to mitigate the effects of multipathing in GPS signals, since it is likely that the signal will still arrive at the position determination terminal via the direct path. That is, the signal between the GPS satellite and the position determination terminal is likely to take more than one path. However, one of those paths is likely to be the direct path. Accordingly, the direct path is assumed to be the one that is first to arrive. In addition, the direct path typically will have greater signal strength. In contrast, communication signals transmitted from a position determination beacon are more likely to take only indirect paths.
Therefore, there is a need to determine the error that is introduced by multipathing. The following description discloses a method and apparatus for determining an estimate of the amount of error that is present in pseudo-range measurements made in a hybrid position determination system.
SUMMARY OF THE INVENTION
The presently disclosed method and apparatus allows correlations between a selected parameter and the error in a pseudo range measurement to be exploited. A database is established in which the amount of error estimated for particular pseudo range measurements to a beacon is maintained. Clusters are defined. Each cluster is associated with a range of values for the selected parameter. Pseudo range measurements are then associated with a particular cluster based upon the value of the selected parameter at the time (or proximate to the time) the pseudo range measurement was taken. As more estimates of the pseudo range measurements are made, the size of the clusters (i.e., the range of values of the selected parameter) can be reduced. Due to the correlation between the selected parameter and the errors in the pseudo range measurements, reducing the size of the clusters reduces the variance of the error estimates. The mean value of the error estimates is used to correct the error in future pseudo range measurements.
In one embodiment of the disclosed method and apparatus, the position of a terminal that measures the pseudo range to a beacon is the selected parameter. Alternatively, any other correlated parameter, such as the power level of the beacon signal, may be the selected parameter. The size of the cluster is initially relatively large, since the database will have relatively few error estimates in any particular geographic region. However, as the number of error estimates increases, the size of the clusters can be decreased, thus reducing the variance of the error estimates within the smaller clusters with respect to the larger clusters. hi accordance with one embodiment, error estimates are made by first calculating what the pseudo range to a particular beacon should be. This calculation is made by determining (using a highly accurate first position determination sub-system) the current position of the terminal used to make the pseudo range measurements to the beacon. Once the position of the terminal is known, the pseudo range measurements to a beacon can easily be calculated, assuming that the location of the beacon is known. The pseudo range to the beacon from the terminal is then measured using a less accurate second position determination sub-system. The difference between the pseudo range that is calculated based on the more accurate first position determination sub-system and the pseudo range measurement made by the less accurate second position determination sub-system is determined. This difference is assumed to be due to the error in the measurement made by the less accurate second position determination sub-system.
Accordingly, the database includes information that allows pseudo range measurements made by the less accurate second position determination sub-system to be corrected when the more accurate first position determination sub-system is not available. The database is self generating, in that the information required in the database is taken during the operation of the terminal based upon the availability of the more accurate first position determination sub-system. The more points the terminal can locate using the more accurate first position determination sub-system, the smaller the clusters in the database. As a consequence of smaller clusters, the variance in the error estimates that are maintained in the database for each cluster will be reduced. It should be understood that the presently disclosed method and apparatus could be used with position determination systems other than hybrid position determination systems if there is some other means by which to determine the position of the terminal and that means is available at some times or locations, but not at others. In that case, the position of the terminal would be used as a reference to determine the amount of error in the pseudo range measurements in the same way as the position determination from the more accurate sub-system described above.
In accordance with one embodiment of the disclosed method and apparatus, an iterative approach is used when the selected parameter is the location of the terminal. Assuming that the more accurate first position determination sub-system is not available, and that a sufficient number of initial error estimates have been made, the iterative approach uses corrected measurements of the pseudo range based upon a relatively large cluster to determine the position of the terminal. Once the location of the terminal has been determined in this way, the corrections to the pseudo ranges can be recalculated based upon a much smaller cluster, assuming that a statistically valid number of error estimates has been made for the smaller cluster.
BRIEF DESCRIPTION OF THE DRAWINGS
The features, objects, and advantages of the present invention will become more apparent from the detailed description set forth below when taken in conjunction with the drawings, in which like reference characters identify like elements in the figures throughout and wherein:
FIG. 1 is an illustration of a hybrid position determination system in accordance with the disclosed method and apparatus. FIG. 2 illustrates the disclosed method in which there is a relation between the approximate position of a terminal and the expected error in pseudo range measurements made by the terminal. Fig. 3 is an illustration of an entire region shown as a single cluster.
Fig. 4 is an illustration of four smaller clusters that have been subdivided from one larger cluster.
Fig. 5 is an illustration of an example in which some of the clusters of Fig. 4 have been combined to form a new cluster.
Fig. 6 illustrates the process for building an MSD in the case in which there is no assumption regarding the errors in the satellite measurements.
FIG. 7 is a simplified block diagram of one embodiment of an iterative process for correcting errors in position determinations.
FIG. 8 is a simplified block diagram of one embodiment of a terminal used to implement the presently disclosed method.
DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS
FIG. 1 is an illustration of a hybrid position determination system 100. For the purpose of this disclosure, a hybrid position determination system is defined as a system in which the position of a terminal 105 can be determined from a first position determination sub-system alone or in combination with a second sub-system. The first position determination sub-system is preferably capable of determining the position of the terminal 105 "independent" of the second sub-system, the present context, the phrase "independent of the second sub-system" means that the first sub-system may be used to determine the position of a terminal at times when some or all of the components of the second sub-system are not available. However, in accordance with one embodiment of the presently disclosed method and apparatus, components of the second sub-system may be a necessary part of the first sub-system. It should be noted that the second sub-system may or may not be capable of performing a position determination without the use of components of, or information from, the first subsystem.
Furthermore, in a hybrid position determination system in accordance with the presently disclosed method and apparatus, at least one "parameter" is related to the magnitude of the errors in position determination measurements of one of the position determination sub-systems. For the purpose of this disclosure, a parameter is any variable that is related to the magnitude of the errors in position determination measurements of one of the position determination sub-systems. For example, a parameter for a position determination terminal 105 can be: (1) the location of the position determination terminal 105, (2) the amount of power in a position determination signal received by the position determination terminal 105, (3) the number of base stations that are present in the geographic region in which the position determination terminal 105 is presently located, (4) the type of buildings in the geographic region in which the position determination terminal 105 is presently located, (5) the density of the structures in the geographic region in which the position determination terminal 105 is presently located, etc. Each of these variables will effect the accuracy of the measurements made by a position determination terminal within a terrestrial position determination system, and can therefore be considered to be a parameter for the position determination terminal 105. However, a parameter which is related to the magnitude of the errors in measurements made by the second sub-system is preferably unrelated (and alternatively, weakly related) to the magnitude of the errors in measurements made by the first position determination sub-system.
In one example of a hybrid position determination system including two subsystems, each of the two sub-systems includes a set of transmitters. The transmitters of each sub-system transmit position determination signals. In such a hybrid system, it is not necessary for transmissions from both sub-systems to be available at the same time. In the embodiment depicted in FIG.l, four satellites 101, three position determination beacons 103, and one position determination terminal 105 are shown. The satellites 101 are a first of the two types of transmitters of position determination signals and are associated with the first sub-system. The satellites 101 provide signals that can be received by the terminal 105. The received signals enable the terminal to determine the position of the terminal 105 Without the need to receive position determination signals from any of the beacons 103, assuming that signals from a sufficient number of satellites signals are received and that the information received on those signals can be decoded. Nonetheless, in some embodiments of the disclosed method and apparatus, the terminal 105 may require communication with a beacon 103 for acquisition and processing assistance. Whether acquisition assistance is desirable or not, once a sufficient number of satellites 101 are acquired, the information provided by the satellites 101 is sufficient to calculate the position of the terminal 105.
The beacons 103 are a part of a second of the two position determination subsystems. Like the satellites 101, the beacons 103 provide signals that can be received by the terminal 105. In accordance with one embodiment of the disclosed method and apparatus, the signals enable the terminal 105 to determine its own position without the need to receive position determination signals from any of the satellites 101. However, in an alternative embodiment, the second sub-system requires at least one transmitter from the first sub-system to provide information in order to determine the position of the terminal 105. Once the information is received from each of the required transmitters, a determination of the position can be performed by the terminal 105, either alone or with acquisition and processing assistance from one of the beacons 103. Alternatively, the position of the terminal 105 can be determined by a device (not explicitly shown) that is remote from the terminal 105. In this case, the terminal 105 transmits any information necessary to determine the position of the terminal 105 to the remote device. In one embodiment of the disclosed method and apparatus, the remote device is one of the beacons 103 or a sub-system 104 within one of the beacons 103. The remote device is shown in FIG. 1 as "optional" to highlight the fact that the remote device may be located either in the beacon 103 or external to the beacon 103. It will be understood by those skilled in the art that there may be more or less satellites 101 or beacons 103 than are shown in FIG. 1, and that there may be more than one terminal 105. In one embodiment of the hybrid position determination system 100, the satellites 101 are GPS satellites. However, in an alternative embodiment of the hybrid position determination system 100, the satellites 101 may be any other type of transmitter that enables a relatively independent and relatively accurate position determination to be made. For example, a terrestrial position determination system (such as LORAIN) that is capable of providing an accurate position determination may be used instead of satellites.
In one embodiment of the hybrid position determination system 100, the beacons 103 are cellular base stations, including (1) base station transceiver sub-systems (BTSs), (2) base station controllers (BSCs), and (3) position determination equipment (PDEs). However, it will be understood that in other embodiments, the beacons 103 may include a subset of these components. Furthermore, the beacons 103 may be any other transmitters capable of transmitting position determination signals to enable the position of a terminal 105 to be determined. It should also be understood by those skilled in the ait that in a beacon 103 that includes a PDE, the PDE would likely serve as the remote device 104 shown in FIG. 1. In accordance with the presently disclosed method and apparatus, information received from the first type of transmitters should enable a substantially more accurate position determination to be calculated than can be determined from the information received from the second type of transmitters. For example, in the case of GPS satellites and cellular base stations, the accuracy of the position determination made using measurements to the satellites is typically greater than the accuracy of the position determinations made using measurements to the cellular base stations.
For the purpose of clarity, the presently disclosed method and apparatus will be described in the context of a hybrid position determination system in which the satellites are GPS satellites and the beacon is a code division multiple access (CDMA) cellular base station. However, as noted above, the present invention does not rely upon any of the unique features of these systems. Accordingly, the invention may be implemented using any other systems for determining position.
In the embodiment in which the satellites are GPS satellites and the beacons are CDMA cellular base stations (which may or may not include a BSC and PDE) "pseudo range" measurements psi - Ps are made with respect to the satellites 101 and similarly, pseudo range measurement pt>i - pb3 are made with respect to the beacons 103. Pseudo range measurements represent the relative distance between the receiving terminal 105 and the sources of position determination signals. It should be noted that the particular manner in which position determinations are made is not relevant to the presently disclosed method and apparatus. However, a discussion of pseudo range based position determination techniques is provided as an example of one embodiment of the disclosed method and apparatus. Other means for measuring pseudo ranges are well known and could be used to perform the disclosed method and apparatus. Furthermore, other means for determining position without the need for pseudo range measurements are well known and could be used to perform the disclosed method and apparatus.
Pseudo range measurements typically represent a time difference between an arbitrary point in time and the time at which a signal arrives at a position determination terminal 105. However, pseudo ranges are typically presented in units of meters. The amount of time can be translated into a distance in meters by multiplying the time difference by the speed of light in meters per second.
In one embodiment of the disclosed method and apparatus in which the beacon
103 is a CDMA base station, the position determination terminal notes when the beginning of a series of bits that make up a spreading code for the CDMA signal has been received with respect to a prior arbitrary point in time. This series of bits is commonly referred to as a "pseudo-random noise (PN) spreading code".
It should be noted that the signals transmitted by each of the base stations within a CDMA communications system are encoded with the same PN spreading code. However, there are offsets in time between the start of the code transmitted from one base station with respect to the start of the code transmitted from each other base station. Therefore, these offsets must be taken into account before comparing the relative timing of the signals received from each of the beacons 103 with respect to one another. As is well known, these offsets are easily determined and can be subtracted out. In fact, these offsets are typically used to identify the particular base station from which the signal originated.
Furthermore, it should be noted that the duration of the clock (assumed for this example to be one millisecond) must be such that any ambiguity can be resolved. That is, the duration of the clock must be such that the values of the clock that are associated with the received signals should uniquely identify only one relationship between the received signals. Accordingly, the duration of the clock must be much greater than twice the difference between the time at which the signals will be received from each beacon 103, taking into account the offsets in the timing of the codes transmitted by each beacon 103. It should be noted that making pseudo range measurements is very well known in the art.
If the positions are known for at least three of the beacons 103 from which pseudo range measurements are taken, then the pseudo range measurements can be used to determine the position of the terminal 105. The pseudo range measurements and the position of the beacons 103 are applied to a well-known process commonly referred to as the "Least Mean Squares" (LMS) process. Likewise, the inverse is true. That is, if the position of the terminal 105 and three beacons 103 are known, the distance between the terminal 105 and each of the three beacons 103 can be used to determine the pseudo range values from the position of the terminal 105 to the beacons 103. The presently disclosed method and apparatus takes advantage of the fact that a hybrid position determination system will frequently have two independently derived sets of pseudo range measurements. One set of measurements is typically going to be more accurate than the other. Therefore, if the more accurate set of measurements is sufficient to determine the position of the terminal 105, then the more accurate set of measurements can be used to determine the error in the less accurate set of measurements. Furthermore, it has been determined that there is a predictable relationship between particular parameters (such as position of the terminal 105) and the amount of error in the less accurate set of measurements. Therefore, by nowing both the value of the parameter and the amount of error associated with that parameter, an estimate as to the amount of error in the less accurate measurements can be made. For example, by knowing the approximate position of the terminal 105 and the relationship between the position and the amount of error, the amount of error in the less accurate measurement can be estimated. In order to take advantage of the relationship between a particular parameter and the amount of error in the less accurate pseudo range measurements, a "Measurement Statistics Database" (MSD) is generated in accordance with one embodiment of the presently disclosed method and apparatus. It should be noted that means other than a database may be used for associating the values of a parameter with the amount of error in the measurements made by the less accurate position determination sub-system.
However, the MSD provides an effective means by which a parameter (such as the position of terminal 105) can be associated with a correction factor to be applied to the less accurate pseudo range measurements. Accordingly, when the value of the parameter is known, the associated correction factor can be determined from the MSD and applied to the pseudo range measurement.
The following example describes the presently disclosed method and apparatus in the case in which the parameter at issue is the position of the terminal 105 as determined from essentially uncorrected pseudo range measurements made by the terminal 105. In one embodiment of the presently disclosed method and apparatus, the pseudo range measurements are completely uncorrected. However, it will be clear to those skilled in the art that some corrections may be made to the pseudo range measurements without deviating from the scope of the presently described method and apparatus. For example, the pseudo range measurements may be corrected for time offsets due to differences between the time at which the same code is transmitted from different beacons.
It will also be clear to those skilled in the art that the parameter at issue may be any one of a large number of other parameters. It should be noted that it is advantageous for the terminal 105 or the beacon 103 to be able to directly measure the value of the parameter at issue. Direct measurement of the parameter by the terminal 105 or the beacon 103 allows the disclosed method and apparatus to operate without input from an external source. However, it is possible for the value of the parameter to be input by an external source, such as the terminal operator or signals received from a remote source.
Possible parameters include, but are not limited to: (1) the amount of power with which the position determination signals have been received by the terminal 105; (2) the type of buildings that are proximate to the terminal 105; (3) the amount of urban development in the general vicinity of the terminal 105; (4) the distance from the beacons 103 to the terminal 105 that will be receiving the position determination signals; and (5) the shape of the correlation peak determined from the correlation of the received signal with a known PN spreading code. It should be clear that this list merely provides a small sample of the many types of parameters that could be used to predict the amount of error that is likely to be present in pseudo range measurements. The present invention should not be limited to the types of parameters that are listed herein.
Rather, the scope of the present disclosure should be considered to include using any parameter that is correlated with a first set of pseudo range measurements and not (or alternatively, to a lesser degree) with a second set of pseudo range measurements to predict the amount of error in the first set of measurements.
Generation of an MSD
The process for generating an MSD is illustrated in FIG. 2. FIG. 2 particularly relates to the example in which there is a relation between the approximate position of a terminal 105 and the expected error in pseudo range measurements made by the terminal 105.
Initially, a determination is made as to whether the more accurate position determination sub-system is available (STEP 201). If so, then the position of the terminal 105 is determined by using only measurements from the more accurate position determination sub-system (STEP 203). The process of generating the MSD requires that the more accurate position determination capability be available. If not available, then the process of generating an MSD can go no further until the more accurate position determination capability becomes available. Nonetheless, it should be understood that in this case, the MSD, if already sufficiently generated, can be used to estimate the error in the less accurate sub-system.
In one embodiment of the presently disclosed method and apparatus, the terminal 105 has the ability to determine position based upon either pseudo range measurements to satellites 101 or pseudo range measurements to beacons 103. The pseudo range measurements to the satellites 101 tend to be more accurate. Therefore, in such an embodiment, the position of the terminal 105 is determined using pseudo ranges to the satellites 101, if possible. In an alternative embodiment, the beacon pseudo ranges are used only to supplement the satellite pseudo range measurements when an insufficient number of satellites 101 are available. However, similar to the case in which the beacon pseudo ranges can be used with a satellite to determine the position of the terminal 105, using satellite pseudo ranges is preferred over using beacon pseudo ranges.
The position of the terminal 105 is determined from pseudo ranges to the satellites 101 in essentially the same way as described above with respect to the determination of pseudo ranges from the beacons 103. Greater accuracy in the pseudo range measurements to the satellites 101 than can typically be had from the beacons 103 comes from several factors, including the fact that the satellites 101 are overhead. The probability that the signals received from the satellites 101 will reach the terminal 105 on a direct path is greater than the probability that the signals from the beacons 103 will reach the terminal 105 directly. The failure of the signals from the beacons 103 to reach the terminal 105 directly adds to the distance along the path propagated by the signals from the beacon 103 to the terminal 105. The additional distance causes an error in calculating the position of the terminal 105 when using the pseudo ranges from the beacons 103. Once the position of the terminal 105 has been determined using the more accurate position determination sub-system, the pseudo ranges to the beacons 103 can be computed (STEP 207). It can be seen that the expected pseudo range to a particular beacon 103 can be easily determined from knowledge of the position of the terminal 105 and from knowledge of the position of the beacon 103. In one embodiment of the disclosed method and apparatus, the pseudo ranges are provided to the beacon 103, which determines the expected pseudo ranges for a particular terminal 105 to a particular beacon 103. If the terminal does the determination and calculation, then the presently disclosed method and apparatus assumes that the terminal 105 has access to knowledge of the position of the beacons 103. For example, in a system in which the beacons 103 are CDMA base stations, the base stations provide the terminal 105 with information regarding their own position. Alternatively, the terminal 105 maintains a database that provides the terminal 105 with the position of the beacon 103 based upon an identification indication received as part of the transmissions from the beacon 103 to the terminal 105. In yet another embodiment of the disclosed method and apparatus, the terminal 105 sends information to one or more of the beacons 103 to process. The beacon 103 knows or has access to the location of the beacons in the area.
The terminal 105 then measures the pseudo ranges to each beacon 103 (STEP 209). Each measured pseudo range is associated with one beacon 103 and is compared to a pseudo range associated with that beacon 103 as computed in STEP 207 (STEP 211). It is assumed that the difference between a computed and measured pseudo range is an error in the measured pseudo range. It should be clear that the pseudo ranges to the both beacons 103 and the satellites 101 can be measured and sent to the beacon 103 for processing prior to the calculation in STEP 207 of the beacon pseudo ranges.
Once the error in the measured pseudo ranges to the beacons 103 has been determined, the "cluster" in which the terminal 105 is presently located is determined (STEP 213). A cluster can be defined in a number of ways. In accordance with the preferred embodiment, a cluster is defined as a set of contiguous values of the parameter. For example, if the parameter is the position of the terminal 105, then the cluster would be a contiguous geographic region. The set of contiguous parameter values is preferably large enough to include several error measurements. That is, preferably enough error estimates have been made within the cluster to calculate with reasonable accuracy, the mean of the pseudo range errors made by the less accurate sub- system at any position within the cluster.
If there are an insufficient number of error estimates, then the size of the cluster is increased to include more error estimates. Alternatively, the cluster will be considered to be immature until a statistically sufficient number of error estimates are made. At the point at which enough error estimates have been made to support two clusters, the cluster is divided into two. Each cluster should have a sufficient number of error estimates to allow the calculated mean of the error associated with that cluster to reflect with reasonable accuracy the mean that would be calculated from an essentially infinite number of error estimates taken from within the cluster. If the mean of the errors in contiguous clusters are the same, then those clusters may be combined, even though there are sufficient number of error estimates to support more than one cluster. However, in accordance with one embodiment of the disclosed method and apparatus, the cluster may be divided differently if such a division would result in two clusters, each having a sufficient number of error estimates and different error estimate means.
The desired result is to associate each cluster with a value, or set of values, of a parameter. If the value of the parameter detected by the terminal 105 is associated with the cluster, then the terminal 105 is assumed to be in the associated cluster. For example, if the parameter is the position of the terminal 105, then the cluster would be associated with a range of geographical positions. If the terminal 105 is detected to be located within a cluster, then the terminal 105 is considered to be within that cluster for the purposes of the disclosed method and apparatus. Alternatively, assume that the parameter of interest is the power level of the received signal. If the terminal 105 detects that the received signals are within the range of power levels associated with a particular cluster, then that terminal 105 would be considered to be within the associated cluster.
The parameter is selected such that there is a correlation between the parameter and the amount of error in the pseudo range measurements made by the less accurate position determination sub-system. Accordingly, each value of the parameter is associated with a particular value in the error in the position determinations made by the less accurate position determination sub-system. That is, whenever the terminal 105 notes that the parameter has a particular value (or is within a particular range of values), the error in the determination of the position made by the less accurate position determination sub-system will be a particular value (or within a particular range of values).
Refer again to the example in which the position of the terminal 105 is the parameter of interest. Each position of the terminal 105 would be associated with a particular amount of error in the position determination measurement made by the terminal 105 using the less accurate position determination sub-system. Accordingly, to select position as an appropriate parameter, each position of the terminal 105 must be related to a particular amount of error in the position determination measurement made by the less accurate position determination sub-system. Accordingly, for a terminal at position X the amount of error would be Y. This error value would be essentially constant and predictable, within a range of tolerable uncertainty.
As a more explicit example, refer to FIGs. 3 and 4. A region 300 is shown in which Eleven positions 301 are identified. Pseudo range error estimates were made at each of these positions 301. The error estimates are based upon the difference between pseudo range measurements taken at those locations 301 and pseudo ranges calculated from a knowledge of the position based on the more accurate position determination sub-system.
FIG. 3 is an illustration of the entire region 300 shown as a single cluster. This is so because initially there are only relatively few (eleven) error estimates made in the region 300. It should be noted that the particular number of estimates that is required to constitute a sufficient number depends on the particular application. Therefore, the number may be much greater or less than eleven depending upon the application. It should also be noted that while the beacons 103 are shown in FIG. 3 to be outside the region 300, the positional relationship between the beacons 103 and the region is not relevant to the present invention other than noting that it should be possible to receive the signals from the beacons 103 from within the region 300.
Over time, additional error estimates are made at various locations 302 within the region 300. As these additional error estimates are made at locations 400, a sufficient number of error estimates may be available in portions of the cluster to allow the cluster to be subdivided into smaller clusters. FIG. 4 is an illustration of four smaller clusters 401, 403, 405, and 407 which have been subdivided from one larger cluster 300.
The mean of the errors in the pseudo range measurements made in each small cluster 401, 403, 405, and 407 is then determined. If the mean of two or more of the smaller clusters 303, 305, and 307 is sufficiently close in value, then those clusters may be put back together to form one cluster. FIG. 5 is an illustration of an example in which clusters 401, 403, and 407 of FIG. 4 have been combined to form yet another new cluster 501 including positions 500 from which error estimates have been made. Accordingly, clusters having relatively complex shapes can be formed using clusters having relatively simple shapes. Having more error estimates taken at different positions within a cluster provides more opportunities to reshape the clusters.
The measurements taken within a cluster will typically have differing amounts of error. The set of error values form an array associated with a particular beacon 103 within a particular cluster. The values within the array will vary about a mean value within some range. In the terms of statistics, one way to characterize the size of the range in this variation is to calculate the "variance" of the array. Alternatively, the size of the range can be characterized by the "deviation" of the array. Alternatively, the boundaries of a cluster can be determined based upon the value of the error estimates of the pseudo ranges measured. In this case, the variation in the error estimates included in the array is minimized. The area within a cluster is defined such that variations in the error estimates that have been made by position determination terminals 105 within the cluster are within a predetermined range. If a measurement is outside of the predetermined range, then the shape of the cluster is altered to remove the offending position. This results in greater accuracy in determining the correction to be applied to measurements taken by the less accurate position determination sub-system. This is particularly valuable when it is necessary to rely on at least one measurement from the less accurate sub-system in order to determine the position of the terminal 105.
Once a new variance and mean are calculated for the array associated with the particular cluster in which the terminal 105 is located, the new variance and mean replace the old (STEP 217). In the case in which no previous values for the variance and mean had been calculated, the error itself would be the mean and the variance would be equal to zero. Each newer error value associated with that cluster would then be placed into the array of older values. In one embodiment of the presently disclosed method and apparatus, the entire array is maintained within the MSD. Alternatively, only the new variance, mean, and number of elements in the array are maintained. It should be noted that in order to update the variance and mean of the error associated with that cluster the entire array must be available. Alternatively, the mean, variance, and the number of elements in the array must be known. By updating the mean, variance, and number of elements in the array each time a terminal 105 makes a measurement, the MSD is dynamically generated. Alternatively, the values of all of the elements in the array are stored in the MSD. The values for the mean and variance can be calculated as needed from the stored array elements, rather than being calculated after each measurement is made. It should be noted that in accordance with one embodiment of the disclosed method and apparatus, the mean and variance can be weighted by giving more significance to certain elements in the array. In accordance with one such embodiment, the elements that are given more weight are those elements that are known to be more reliable for one reason or another. For example, if a terminal 105 makes a pseudo range measurement based on a relatively weak received signal, the weight of such a pseudo range measurement might be reduced with respect to other measurements made from stronger signals.
Furthermore, it should be noted that in accordance with one embodiment of the presently disclosed method and apparatus, it can be assumed that there is essentially no correlation between the errors in the satellite pseudo range measurements and the errors in the beacon pseudo range measurements. Therefore, in light of this assumption, the satellite errors (i.e., the errors in the position determined using the satellite pseudo range measurements) will typically not bias the mean of the errors in a cluster. However, in accordance with another embodiment of the presently disclosed method and apparatus, no such assumption is made.
FIG. 6 illustrates the process for building an MSD in the case in which there is no assumption regarding the errors in the satellite measurements, as described below.
Pseudo range measurements to the satellites 101 are made by the terminal 105 (STEP 601). In addition, the terminal notes the variance in the satellite measurements (STEP 603). That is, since the measurements will have a certain amount of error, the measurements will vary about a mean within some range. This variation is characterized by calculating the variance in the error in the satellite measurements. Alternatively, the errors can be characterized by their deviation. The variance (or deviation) and the mean in the errors in the satellite measurements are typically known, hi general, the variance and means are provided in a covariance matrix associated with the each satellite 101, as is well known in the art. The satellite measurements and the error statistics associated with the satellite measurements (typically in the form of the covariance matrix) are applied to a well- known LMS process (STEP 604). In addition, in one embodiment of the disclosed method, the satellite geometry is also taken into account by the LMS process. The result of applying the satellite measurements, the error statistics associated with those measurements, and the satellite geometry to the LMS process is a position determination solution. The position determination solution provides the position of the terminal 105 from which the pseudo range measurements were taken. In addition, the LMS process results in a solution covariance matrix that indicates the error statistics regarding the amount of error in the solution. Such solution covariance matrices are well know in the art. The means are assumed to be zero due to the fact that there is assumed to be no multipathing errors in the signals received from the satellites.
The solution is presented in the form x, y, z, and b, where x, y, and z are the Cartesian coordinates of the terminal 105 with respect to a selected reference, such as the center of the earth, and b is the bias of the clock in the terminal 105 with respect to true GPS time. From the position determination solution, an estimate can be made of the pseudo range measurements to each of the beacons 103 from which the terminal 105 is receiving position determination signals (STEP 606). The pseudo range estimates can then be compared with actual pseudo range measurements taken by the terminal 105 to each beacon 103 (STEP 608) to determine the difference between them. This difference is assumed to be an estimate of the error in the measured pseudo ranges. In addition, a beacon covariance matrix can be calculated from the solution covariance matrix (STEP 610). The beacon covariance matrix characterizes the amount of error in each of the calculated error estimates based upon the amount of error in the satellite pseudo range measurements, as determined from the satellite covariance matrix. That is, estimates of the error in the beacon pseudo range measurements are based on the assumption that the calculated pseudo ranges are perfectly accurate. However, the statistics of the error is known for the satellite measurements used to derive the calculated pseudo ranges. Therefore, there is an error in the estimates of the amount of error in the beacon pseudo ranges. The statistics of the error in the error estimates can be determined from the solution covariance matrix, which is in turn calculated from the satellite covariance matrix. This can be done in accordance with methods well known to those skilled in the art. In accordance with this embodiment of the presently disclosed method and apparatus, the MSD is generated such that each entry to the MSD is associated with a cluster, a group of beacons, the mean of the array of pseudo range measurement error estimates to each beacon, and the beacon covariance matrix that characterizes the error in each measurement of the pseudo range due to the errors introduced by multipathing (STEP 612). Alternatively, each entry includes identification of an associated cluster, identification of an associated beacon, an array of pseudo range measurement error estimates associated with the beacon, and beacon covariance matrix for the associated beacon. It should be noted that it may not be possible to determine the covariance between beacons associated with different entries. Accordingly, if there is only one beacon associated with an entry, then it may not be possible to determine the covariance matrix for that entry. As noted above, in accordance with one embodiment of the presently disclosed method and apparatus, the error estimates of the array are weighted according to the reliability of each such estimate.
It should be understood that by selecting clusters that cover a larger amount of area, the variance within each cluster will typically become greater. This is due to the more diverse conditions that would exist in a larger cluster. However, the advantage of having relatively large clusters is that less data is required (i.e., fewer points per square kilometer) to accurately determine the variance and mean of the error within the cluster.
Once the variance and mean of the errors associated with a particular beacon within a cluster are known, these values can be used. In particular, they are used to estimate the error in pseudo ranges to the particular beacons 103 when the pseudo range measurements are being made by a terminal 105 within the cluster. The estimates of the error in the beacon pseudo range measurements can then be used to correct the measured beacon pseudo ranges. This is particularly useful when there are too few satellite pseudo range measurements. Accordingly, if there are not enough satellites available to make a position determination based on satellite measurements alone, then corrected pseudo ranges to beacons 103 can be used. It should be understood that the presently disclosed method and apparatus assumes that the MSD will be generated and stored within the terminal 105. However, the MSD may be generated and stored within a beacon or other component remote from the terminal 105. For example, it will be clear to those skilled in the art that the terminal 105 may be responsible only for taking pseudo range measurements. These measurements are then transmitted to the beacon 103. The beacon 103 then processes the pseudo range measurements as described above in a device 104 within the beacon 103. Alternatively, the device 104 is remote from the beacon 103 and the beacon 103 communicates the pseudo range measurements to the remote device 104. The remote device 104 then processes the pseudo range measurements. It should be understood that the processing may be distributed such that some of the processing is performed in one device and other processing is performed in other devices. Nonetheless, without regard for which device is responsible for performing which process, the disclosed method will be performed essentially as defined above.
As contemplated above, the present method and apparatus provides a means by which a position determination from a less accurate position determination sub-system can be corrected in a single correction step. However, in an alternative embodiment of the process described above, an iterative process can be performed as shown in FIG. 7 and described below.
Initially, the position of a terminal 105 is determined using a less accurate subsystem (STEP 701). That determination is then corrected based on a correction factor associated with a first relatively large cluster (STEP 703). The first cluster is intentionally made to be relatively large to account for the fact that the error that is present in the less accurate sub-system's initial estimate of the position may be great enough to make it difficult to accurately determine which cluster the terminal 105 is in if the cluster is made to small. Once the correction to the position is made using the correction factor from the first relatively large cluster, the position of the terminal 105 will be known to a greater accuracy. Smaller clusters are then used to determine a more precise correction factor (STEP 705), since the improvement in accuracy provided by the first correction will make it possible to determine in which smaller cluster the terminal 105 current resides.
It should be understood that the presently disclosed method and apparatus could be used with position determination systems other than hybrid position determination systems if there is some other means by which to determine the position of the terminal is available at some times or locations, but not at others. In this case, the position of the terminal would be used as a reference to determine the amount of error in the pseudo range measurements in the same way that the position determination from the more accurate sub-system described above is used.
Apparatus Used to Implement one Embodiment of the Presently Disclosed Method FIG. 8 is a simplified block diagram of one embodiment of a terminal 105 used to implement the presently disclosed method. The terminal 105 includes a receiver 501, a decoder 503, a correlator 505, a reference code generator 507, a processor 509, a clock 511, and a memory 513. The receiver 501 typically includes a conventional radio frequency front end section, including antenna, down converter, filters, amplifiers, etc., all of which are well known and for simplicity sake are not shown. Signals, including position determination signals, are received by the receiver 510. The received signals are down-converted, filtered, and amplified (not necessarily in that order) as required for input to the decoder 503. It will be understood by those skilled in the art that other processing not disclosed herein may occur in the receiver 501, but is essentially immaterial to the novel aspects of the presently disclosed method and apparatus.
The output from the receiver 501 is coupled to the decoder 503. The decoder 503 decodes the intelligence from the received signal. When position determination signals are received and decoded, the output from the decoder will be a code sequence. The decoded intelligence is coupled to the correlator 505. The correlator 505 is also coupled to the reference code generator 507. The reference code generator 507 provides a reference code sequence to the correlator 505 to allow the correlator 505 to identify a predetermined point in the received code sequence. The correlator 505 outputs a signal to the processor 509 indicating when the reference code sequence and the received code sequence are synchronized. The reference code generator 507 is also coupled to the processor 509. Once the reference code sequence and the received code sequence are synchronized, the output from the reference code generator 507 to the processor 509 allows the processor 509 to determine the timing of the received code sequence.
The clock 511 is also coupled to the processor 509. The clock 511 may be either free running or synchronized to an external clock reference. If the clock is synchronized to an external clock reference, then the timing of the received code sequence can be determined with respect to the external clock reference. The processor 509 is also coupled to the receiver 501. The receiver 501 provides the processor 509 with information as to which frequency the receiver 501 is be tuned. In accordance with one type of received signal, the processor 509 receives information from the timing of the received code sequence that, together with the frequency on which the signal was received, is indicative of the particular transmitter 101, 103 from which the code sequence originated. It should be noted that multiple correlators 505 may be used to concurrently correlate multiple signals. The processor 509 takes all of the information that is provided to it and determines the relative timing (i.e., the pseudo ranges) of all of the signals received. In addition, the processor 509 performs all of the pseudo range processing functions described above. One application for the presently disclosed method and apparatus allows more accurate production of a "sector coverage map". A sector coverage map is a map that shows which geographic regions are being served by a particular sector of a beacon 103. Beacons 103 are said to be sectorized (i.e., serve a plurality of sectors) when separate antennas are used to direct the signals transmitted by the server to a particular geographic sector and to receive signals from that particular geographic sector. Typically, the sectors at issue are assumed to be essentially pie shaped wedges emanating out from the beacon 103. It is common for cellular base stations to be constructed with three such sectors, each sector forming a non-overlapping pie shaped wedge with a width of 120 degrees. However, due to several factors, sectors of a beacon 103 may, in fact, receive signals from, and transmit signals to, areas that would not typically be assumed to be within the sector. Accordingly, a sector coverage map allows both a terminal 105 and a beacon 103 to determine which sector of the beacon 103 is servicing the terminal 105.
A sector coverage map is generated in accordance with the presently disclosed method and apparatus by associating (1) the sector in which the beacon 103 receives signals from the terminal 105, (2) the position of the terminal 105 as determined by the more accurate position determination sub-system, and (3) the position of the terminal 105 as determined by the less accurate position determination sub-system. Accordingly, a sector coverage map generated as indicated above can be used together with knowledge of the position of the terminal, as determined by the less accurate subsystem, to more reliably attain the sector of the beacon 105 that is servicing the terminal. This is particularly useful when more than one sector receives signals from a terminal 105. This is even more useful when more than one sector is reporting service to the terminal 105 and the position determined by the less accurate position determination sub-system determines that the terminal 105 is not in any of the sectors that are reporting that they are serving the terminal 105.
The previous description of the preferred embodiments of the disclosed method and apparatus is provided to enable any person skilled in the art to make or use the present invention recited in the claims presented below. Various modifications to these embodiments can be made without extending beyond the scope of the invention claimed below. Accordingly, it should be recognized that the principles disclosed herein may be applied to other undisclosed embodiments of the disclosed method and apparatus without extending beyond the scope of the invention recited in the claims provided below. The present invention is not intended to be limited to the embodiments shown herein, but is instead intended to be accorded the widest scope consistent with claims provided below.

Claims

WHAT IS CLAIMED IS:
1. A method for determining an error estimate using information provided by a first sub-system, the error estimate being used to correct errors in position determinations made by a second sub-system, including the steps of: a) determining the position of a terminal using the first sub-system; b) calculating pseudo ranges to at least one transmitter within the second sub- system using the determined position; c) measuring pseudo ranges to the at least one transmitter within the second sub- system; and d) comparing the calculated and measured pseudo ranges for the at least one transmitter to determine an error estimate.
2. The method of Claim 1, wherein the first position determination sub-system typically determines the position of a terminal with more accuracy than the second sub-system.
3. A method for using an aggregate value to correct a pseudo range measurement made by a position determination sub-system, including the steps of: a) determining a set of error estimates associated with a cluster, each error estimate within the set representing the error in a pseudo range to at least one transmitter, and each pseudo range being measured from a position within the cluster; b) calculating an aggregate value for the set of error estimates; and c) using the aggregate value to correct for errors in additional pseudo ranges measured within the cluster.
4. The method of Claim 3, wherein the aggregate value is calculated by calculating the mean of the set of error estimates.
5. The method of Claim 4, wherein the set of error estimates is weighted before calculating the mean.
6. The method of Claim 5, wherein the weighting is done based upon the relative reliability of each error estimate.
7. The method of Claim 3, further including the steps of: a) associating a cluster with a range of values of a parameter; and b) determining that an error estimate is associated with the cluster if a value of the parameter determined by a terminal to which a pseudo range associated with the error estimate is taken is within the range.
8. The method of Claim 1, further including the steps of: a) determining error estimates associated with pseudo ranges to at least one transmitter from several positions; b) defining the scope of a cluster such that variations in the error estimates that have been made by position determination terminals within the cluster are within a predetermined range.
9. A method for determining an error estimate in a first sub-system using information provided by a second sub-system, including the steps of: a) determining the position of a position determination terminal using the second sub-system; b) determining second sub-system error statistics regarding the amount of error in the determined position; c) calculating first sub-system pseudo ranges to at least one transmitter using the determined position; d) determining from the second sub-system error statistics, error statistics for the calculated first sub-system pseudo ranges; e) measuring first sub-system pseudo ranges; and f) comparing the calculated and measured first sub-system pseudo ranges to determine an error estimate.
10. The method of Claim 9, wherein the first position determination sub-system typically determines the position of a terminal with more accurate than the second sub-system.
11. The method of Claim 9, further including the steps of: a) determining at least one set of error estimates from several positions within a cluster in accordance with the method recited in Claim 9, each error estimate of the same set being associated with the same beacon and having been taken by a terminal within the same cluster; b) calculating an aggregate value for at least one set of error estimates; and c) using the aggregate value to correct for errors in subsequent measurements of pseudo ranges taken within the cluster and associated with the same beacon to which the aggregate value is associated.
12. The method of Claim 11, wherein the aggregate value is calculated by calculating the mean of the set of error estimates.
13. The method of Claim 12, wherein the set of error estimates is weighted before calculating the mean.
14. The method of Claim 13, wherein the weighting is done based upon the relative reliability of the calculated pseudo ranges used to determine the position of the position determination terminal.
15. The method of Claim 11, further including the steps of: a) associating a cluster with a range of values of a parameter; and b) determining that a position determination terminal is within a cluster if the value of the parameter for the position determination terminal is within the range of values associated with that cluster.
16. The method of Claim 11, further including the steps of: a) determining error estimates associated with pseudo ranges to at least one transmitter from several positions; and b) defining the scope of a cluster such that variations in the error estimates that have been made by position determination terminals within the cluster are within a predetermined range.
17. A method for correcting errors in a pseudo range measurement, including the steps of: a) measuring pseudo ranges to beacons; b) measuring the value of a parameter; c) determining the amount of correction to apply to the measured pseudo ranges to each beacon based upon a predetermined relationship between the value of the parameter and the correction factor to be applied; and d) applying the determined amount of correct to each of the measured pseudo ranges.
18. A method for correcting errors in a pseudo range measurement, including the steps of: a) measuring pseudo ranges to a particular beacon; b) measuring the value of a parameter; c) determining the amount of correction to apply to the measured pseudo ranges based upon the particular beacon, a predetermined relationship between the value of the parameter and the correction factor to be applied for the particular beacon; and d) applying the determined amount of correct to each of the measured pseudo ranges.
19. The method of Claim 18, wherein the parameter is a position of the terminal as determined from uncorrected pseudo range measurements.
20. A method for correcting errors in a pseudo range measurement, including the steps of: a) measuring pseudo ranges to a particular beacon; b) measuring the value of a parameter; c) determining which cluster a terminal current resides in based upon the value of the measured parameter; d) determining the amount of correction to apply to the measured pseudo ranges to the particular beacon based upon the cluster in which the terminal resides; and e) applying the determined amount of correct to each of the measured pseudo ranges to the particular beacon.
21. A method for generating a sector coverage map, including the steps of: a) determining the position of a terminal using a first position determination subsystem, the first position determination sub-system being independent of the any pseudo range measurements to the beacons; b) determining which sectors of a beacon are servicing the terminal; and c) for each sector, creating a sector coverage map that indicates each of the positions, determined by the satellite position determination sub-system, from which a sector can provide service to a terminal.
22. The method of Claim 21, further including the steps of: a) determining the position of the terminal using a second position determination sub-system while the terminal remains at the position from which the position was determined by the first position determination sub-system; and b) associating the location determined by the second position determination system with the location determined by the first position determination sub-system.
23. A computer readable medium encoded with a data structure for determining in which particular sector of a beacon a terminal resides, including: a) a first list of positions associated with a particular sector of a beacon, each such position having been determined by a first position determination sub-system, from which a terminal may be serviced by the associated sector; and b) a second list of positions associated with the particular sector, each such position having been determined by a second position determination sub-system, from which a terminal may be serviced by a particular sector associated with the list of positions.
24. A computer readable medium encoded with a data structure for determining a correction factor to be applied to a pseudo range measurement taken by a first position determination sub-system, including: a) a set of error estimates to be used to correct a pseudo range measurement to an associated beacon; and b) an indication associated with the set of error estimates for identifying the associated beacon.
25. A computer readable medium encoded with a data structure, for determining a correction factor to be applied to a pseudo range measurement taken by a first position determination sub-system, including: a) an aggregate error estimate to be used to correct a pseudo range measurement to an associated beacon, the aggregate error estimate having been calculated from a set of error estimates, each error estimate within the set of error estimates having been calculated from the difference between a measured and a calculated pseudo range to the associated beacon; and b) an indication associated with the aggregate error estimate for identifying the associated beacon.
26. The computer readable medium of Claim 25, further including statistics for characterizing the amount of error in the error estimates.
27. The computer readable medium of Claim 26, wherein the statistics include a covariance matrix.
28. The computer readable medium of Claim 26, wherein the statistics include a mean of the set of error estimates.
29. The computer readable medium of Claim 28, wherein the statistics include a variance of the set of error estimates.
PCT/US2001/050896 2000-10-26 2001-10-24 Method and apparatus for determining an error estimate in a hybrid position determination system WO2002046788A2 (en)

Priority Applications (10)

Application Number Priority Date Filing Date Title
EP01987533A EP1328822B1 (en) 2000-10-26 2001-10-24 Method and apparatus for determining an error estimate in a hybrid position determination system
DE60139493T DE60139493D1 (en) 2000-10-26 2001-10-24 METHOD AND DEVICE FOR ERROR EVALUATION IN A HYBRID LOCATION SYSTEM
AU2002239736A AU2002239736B2 (en) 2000-10-26 2001-10-24 Method and apparatus for determining an error estimate in a hybrid position determination system
IL15563501A IL155635A0 (en) 2000-10-26 2001-10-24 Method and apparatus for determining an error estimate in a hybrid position determination system
CA002426341A CA2426341A1 (en) 2000-10-26 2001-10-24 Method and apparatus for determining an error estimate in a hybrid position determination system
AU3973602A AU3973602A (en) 2000-10-26 2001-10-24 Method and apparatus for determining an error estimate in a hybrid position determination system
AT01987533T ATE438868T1 (en) 2000-10-26 2001-10-24 METHOD AND DEVICE FOR ERROR ASSESSMENT IN A HYBRID LOCATION SYSTEM
KR1020037005790A KR100829666B1 (en) 2000-10-26 2001-10-24 Method and apparatus for determining an error estimate in a hybrid position determination system
JP2002548471A JP4414136B2 (en) 2000-10-26 2001-10-24 Method and apparatus for determining error estimates in a hybrid position determination system
IL155635A IL155635A (en) 2000-10-26 2003-04-28 Method and apparatus for determining an error estimate in a hybrid position determination system

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US09/697,781 US6420999B1 (en) 2000-10-26 2000-10-26 Method and apparatus for determining an error estimate in a hybrid position determination system
US09/697,781 2000-10-26

Publications (2)

Publication Number Publication Date
WO2002046788A2 true WO2002046788A2 (en) 2002-06-13
WO2002046788A3 WO2002046788A3 (en) 2003-01-09

Family

ID=24802506

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2001/050896 WO2002046788A2 (en) 2000-10-26 2001-10-24 Method and apparatus for determining an error estimate in a hybrid position determination system

Country Status (11)

Country Link
US (1) US6420999B1 (en)
EP (2) EP1328822B1 (en)
JP (4) JP4414136B2 (en)
KR (1) KR100829666B1 (en)
CN (4) CN101101326A (en)
AT (1) ATE438868T1 (en)
AU (2) AU2002239736B2 (en)
CA (1) CA2426341A1 (en)
DE (1) DE60139493D1 (en)
IL (2) IL155635A0 (en)
WO (1) WO2002046788A2 (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1294773C (en) * 2001-09-13 2007-01-10 艾尔塞奇公司 System and method for providing traffic information using operational data of a wireless network
WO2012064860A1 (en) * 2010-11-09 2012-05-18 Apple Inc. Beacon-based geofencing
US8504061B2 (en) 2010-04-07 2013-08-06 Apple Inc. Multi-tier geofence detection
US8971913B2 (en) 2003-06-27 2015-03-03 Qualcomm Incorporated Method and apparatus for wireless network hybrid positioning
US9014719B2 (en) 2002-10-17 2015-04-21 Qualcomm Incorporated Method and apparatus for improving radio location accuracy with measurements
US9042917B2 (en) 2005-11-07 2015-05-26 Qualcomm Incorporated Positioning for WLANS and other wireless networks
USRE45808E1 (en) 2004-06-18 2015-11-17 Qualcomm Incorporated Method and apparatus for determining location of a base station using a plurality of mobile stations in a wireless mobile network
US9226257B2 (en) 2006-11-04 2015-12-29 Qualcomm Incorporated Positioning for WLANs and other wireless networks
EP2901737A4 (en) * 2012-09-27 2016-06-15 Ericsson Telefon Ab L M Detecting multipath and determining positioning measurement uncertainty
US9749876B2 (en) 2003-06-27 2017-08-29 Qualcomm Incorporated Local area network assisted positioning

Families Citing this family (63)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040143392A1 (en) 1999-07-12 2004-07-22 Skybitz, Inc. System and method for fast acquisition reporting using communication satellite range measurement
US8255149B2 (en) 1999-07-12 2012-08-28 Skybitz, Inc. System and method for dual-mode location determination
US20070200752A1 (en) * 2001-06-06 2007-08-30 Global Locate, Inc. Method and apparatus for maintaining integrity of long-term orbits in a remote receiver
US6429809B1 (en) * 2001-01-30 2002-08-06 Qualcomm Incorporated Method and apparatus for determining location using a coarse position estimate
US6865394B2 (en) * 2001-01-31 2005-03-08 Hitachi, Ltd Location detection method, location detection system and location detection program
JP4417583B2 (en) * 2001-05-08 2010-02-17 パイオニア株式会社 Navigation device
US8358245B2 (en) * 2001-06-06 2013-01-22 Broadcom Corporation Method and system for extending the usability period of long term orbit (LTO)
US7499710B2 (en) * 2001-08-30 2009-03-03 Alcatel-Lucent Usa Inc. Integrity monitoring for geo-location systems
US7623871B2 (en) * 2002-04-24 2009-11-24 Qualcomm Incorporated Position determination for a wireless terminal in a hybrid position determination system
US20040203989A1 (en) * 2002-09-12 2004-10-14 Broadcom Corporation Using location information to control transmission signal levels of wireless devices
US20040080454A1 (en) * 2002-10-23 2004-04-29 Camp William O. Methods and systems for determining the position of a mobile terminal using digital television signals
US6839020B2 (en) * 2003-06-02 2005-01-04 Motorola, Inc. Aiding location determinations in satellite positioning system receivers
BR122018004395B1 (en) * 2003-06-27 2020-11-10 Qualcomm Incorporated method and equipment for hybrid wireless network positioning
JP2005035532A (en) * 2003-06-30 2005-02-10 Komatsu Ltd Crawler tension adjustment device
FI115167B (en) * 2003-09-10 2005-03-15 Nokia Corp Method and system for positioning and device
US7490008B2 (en) * 2004-09-17 2009-02-10 Itt Manufacturing Enterprises, Inc. GPS accumulated delta range processing for navigation applications
WO2006069597A1 (en) * 2004-12-27 2006-07-06 Telecom Italia S.P.A. Hybrid locating method and system for locating a mobile terminal in a wireless communications network
US7479922B2 (en) * 2005-03-31 2009-01-20 Deere & Company Method and system for determining the location of a vehicle
US7653483B2 (en) 2005-03-31 2010-01-26 Deere & Company System and method for determining a position of a vehicle
US7720598B2 (en) 2005-03-31 2010-05-18 Deere & Company System and method for determining a position of a vehicle with compensation for noise or measurement error
US7647177B2 (en) * 2005-03-31 2010-01-12 Deere & Company System and method for determining a position of a vehicle
US7593811B2 (en) * 2005-03-31 2009-09-22 Deere & Company Method and system for following a lead vehicle
US7705775B2 (en) * 2005-12-30 2010-04-27 General Motors Llc Method of improving a vehicle emergency call network
US7925278B2 (en) * 2006-06-27 2011-04-12 Motorola Mobility, Inc. Method and system for locating a wireless device in a wireless communication network
JP2009543074A (en) 2006-07-07 2009-12-03 スカイフック ワイヤレス,インク. System and method for collecting information from a WLAN-enabled access point to estimate the location of a WLAN positioning device
DE102006033744A1 (en) * 2006-07-21 2008-01-24 Deutsche Telekom Ag Method and device for merging traffic data with incomplete information
US7856234B2 (en) 2006-11-07 2010-12-21 Skyhook Wireless, Inc. System and method for estimating positioning error within a WLAN-based positioning system
CN101234169B (en) * 2007-02-01 2012-05-23 河北以岭医药研究院有限公司 Application of Chinese medicine composition in preparing medicament for adjusting cardiac muscle cell potassium ion channel
US20090231192A1 (en) * 2008-03-14 2009-09-17 Van Diggelen Frank Method and system for generating temporary ephemeris
EP2283641B1 (en) * 2008-06-06 2020-08-12 Skyhook Wireless, Inc. Method and system for determining location using a hybrid satellite and wlan positioning system by selecting the best wlan-ps solution
AU2012200417B2 (en) * 2008-06-06 2014-05-01 Skyhook Wireless, Inc. Method and system for determining location using a hybrid satellite and WLAN positioning system by selecting the best WLAN-PS solution
US8155666B2 (en) 2008-06-16 2012-04-10 Skyhook Wireless, Inc. Methods and systems for determining location using a cellular and WLAN positioning system by selecting the best cellular positioning system solution
JP2010151629A (en) * 2008-12-25 2010-07-08 Ntt Docomo Inc Positioning system and positioning method
US8917206B2 (en) * 2009-06-09 2014-12-23 Qualcomm Incorporated Mobile-based positioning with non-conforming use of assistance data
US8022877B2 (en) 2009-07-16 2011-09-20 Skyhook Wireless, Inc. Systems and methods for using a satellite positioning system to detect moved WLAN access points
US8063820B2 (en) * 2009-07-16 2011-11-22 Skyhook Wireless, Inc. Methods and systems for determining location using a hybrid satellite and WLAN positioning system by selecting the best SPS measurements
DE102010011982A1 (en) * 2009-09-29 2011-04-07 Siemens Aktiengesellschaft Method for the computer-aided creation and / or updating of a reference map for a satellite-based location of an object
US8638256B2 (en) 2009-09-29 2014-01-28 Skyhook Wireless, Inc. Accuracy and performance of a hybrid positioning system
US20110080318A1 (en) * 2009-10-02 2011-04-07 Skyhook Wireless, Inc. Determining A Dilution of Precision Metric Using Two or Three GPS Satellites
US8279114B2 (en) * 2009-10-02 2012-10-02 Skyhook Wireless, Inc. Method of determining position in a hybrid positioning system using a dilution of precision metric
US8259010B2 (en) * 2009-10-14 2012-09-04 Qualcomm Incorporated Qualifying coarse position injection in position determination systems
KR101162585B1 (en) * 2010-01-08 2012-07-04 동명대학교산학협력단 Trianglulation, system, and storage using bias modelling of measured distance
US8531332B2 (en) 2010-03-22 2013-09-10 Qualcomm Incorporated Anti-spoofing detection system
KR101972606B1 (en) 2010-11-03 2019-04-25 스카이후크 와이어리스, 인크. Method of system for increasing the reliability and accuracy of location estimation in a hybrid positioning system
US8547870B2 (en) 2011-06-07 2013-10-01 Qualcomm Incorporated Hybrid positioning mechanism for wireless communication devices
US8509809B2 (en) 2011-06-10 2013-08-13 Qualcomm Incorporated Third party device location estimation in wireless communication networks
US8909244B2 (en) 2011-06-28 2014-12-09 Qualcomm Incorporated Distributed positioning mechanism for wireless communication devices
US8489114B2 (en) 2011-09-19 2013-07-16 Qualcomm Incorporated Time difference of arrival based positioning system
US8521181B2 (en) 2011-09-19 2013-08-27 Qualcomm Incorporated Time of arrival based positioning system
US8457655B2 (en) 2011-09-19 2013-06-04 Qualcomm Incorporated Hybrid time of arrival based positioning system
JP5871536B2 (en) * 2011-09-27 2016-03-01 三菱重工業株式会社 Position specifying device, position specifying method, and program
US8755304B2 (en) 2011-10-21 2014-06-17 Qualcomm Incorporated Time of arrival based positioning for wireless communication systems
US8824325B2 (en) * 2011-12-08 2014-09-02 Qualcomm Incorporated Positioning technique for wireless communication system
US9606217B2 (en) * 2012-05-01 2017-03-28 5D Robotics, Inc. Collaborative spatial positioning
JP5932603B2 (en) * 2012-10-23 2016-06-08 株式会社日立製作所 Extracting regional trends in life prediction
US20140258201A1 (en) * 2013-03-05 2014-09-11 Qualcomm Incorporated Generating a geofence via an analysis of a gps fix utilization distribution
WO2014151118A1 (en) * 2013-03-15 2014-09-25 Nextnav, Llc Techniques to improve the performance of a fixed, timing-based radio positioning network using external assistance information
EP3193187A1 (en) * 2016-01-13 2017-07-19 ST SPORTSERVICE GmbH Method for calibrating a local positioning system based on time-difference-of-arrival measurements
US11740365B2 (en) 2018-02-19 2023-08-29 Mitsubishi Heavy Industries Machinery Systems, Ltd. Positional error prediction device, prediction model generation device, positional error prediction method, prediction model generation method, and program
WO2020230645A1 (en) * 2019-05-13 2020-11-19 日本電気株式会社 Position estimation system, position estimation method, program, and recording medium
US11501245B2 (en) 2019-09-16 2022-11-15 P44, Llc Systems and methods for imputation of shipment milestones
US11378697B2 (en) * 2019-10-30 2022-07-05 Qualcomm Incorporated Multipath mitigation for multiband GNSS receiver
CN113466910A (en) * 2021-05-12 2021-10-01 广西科技大学 Signal positioning system based on WiFi + GPS algorithm

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5717406A (en) * 1995-06-07 1998-02-10 Sanconix Inc. Enhanced position calculation
US5899957A (en) * 1994-01-03 1999-05-04 Trimble Navigation, Ltd. Carrier phase differential GPS corrections network
WO1999061934A1 (en) * 1998-05-28 1999-12-02 Ericsson Inc. Location system combining ranging measurements from gps and cellular networks
US5999124A (en) * 1998-04-22 1999-12-07 Snaptrack, Inc, Satellite positioning system augmentation with wireless communication signals
EP1148344A1 (en) * 2000-04-20 2001-10-24 Lucent Technologies Inc. Positioning of a wireless terminal with satellite positioning signals or base station signals

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0693650B2 (en) * 1988-08-05 1994-11-16 日本電信電話株式会社 Moving body position detection method
JPH0675036A (en) * 1992-08-27 1994-03-18 Oki Electric Ind Co Ltd Position detector
JPH0743446A (en) * 1993-08-02 1995-02-14 Aisin Seiki Co Ltd Verifying apparatus for gps satellite data
JPH07146351A (en) * 1993-11-24 1995-06-06 Sumitomo Electric Ind Ltd Position detector
JP3439827B2 (en) * 1994-04-22 2003-08-25 三菱電機株式会社 Positioning device and moving object with positioning device
US5602903A (en) * 1994-09-28 1997-02-11 Us West Technologies, Inc. Positioning system and method
US5913170A (en) * 1994-11-16 1999-06-15 Highwaymaster Communications, Inc. Locating system and method using a mobile communications network
JPH0974585A (en) * 1995-09-06 1997-03-18 Nippon Motorola Ltd Location detector for mobile station
JP3557024B2 (en) * 1995-12-18 2004-08-25 日本無線株式会社 Positioning device
US5828336A (en) * 1996-03-29 1998-10-27 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Robust real-time wide-area differential GPS navigation
US5995043A (en) * 1996-12-06 1999-11-30 The Boeing Company Aircraft satellite navigation precision-approach system including CDMA datalink
JPH10210548A (en) * 1997-01-20 1998-08-07 Matsushita Electric Ind Co Ltd Mobile communication equipment
US6151551A (en) * 1997-03-03 2000-11-21 Motorola, Inc. Method and apparatus for generating an indication of loss of positioning integrity in emergency call systems
FI973569A (en) * 1997-09-01 1999-03-02 Nokia Telecommunications Oy Subscriber information management system
JPH11186955A (en) * 1997-12-18 1999-07-09 Canon Inc Portable communication terminal
US6075987A (en) * 1998-02-27 2000-06-13 Ericsson Inc. Stand alone global positioning system (GPS) and method with high sensitivity
US6014102A (en) * 1998-04-17 2000-01-11 Motorola, Inc. Method and apparatus for calibrating location finding equipment within a communication system
JP2000244967A (en) * 1999-02-24 2000-09-08 Mitsubishi Electric Corp Mobile communication system, mobile unit and base station configuring the system and method for detecting position of the mobile unit in the system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5899957A (en) * 1994-01-03 1999-05-04 Trimble Navigation, Ltd. Carrier phase differential GPS corrections network
US5717406A (en) * 1995-06-07 1998-02-10 Sanconix Inc. Enhanced position calculation
US5999124A (en) * 1998-04-22 1999-12-07 Snaptrack, Inc, Satellite positioning system augmentation with wireless communication signals
WO1999061934A1 (en) * 1998-05-28 1999-12-02 Ericsson Inc. Location system combining ranging measurements from gps and cellular networks
EP1148344A1 (en) * 2000-04-20 2001-10-24 Lucent Technologies Inc. Positioning of a wireless terminal with satellite positioning signals or base station signals

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1294773C (en) * 2001-09-13 2007-01-10 艾尔塞奇公司 System and method for providing traffic information using operational data of a wireless network
US9014719B2 (en) 2002-10-17 2015-04-21 Qualcomm Incorporated Method and apparatus for improving radio location accuracy with measurements
US9778372B2 (en) 2003-06-27 2017-10-03 Qualcomm Incorporated Wireless network hybrid positioning
US10841892B2 (en) 2003-06-27 2020-11-17 Qualcomm Incorporated Local area network assisted positioning
US9814016B2 (en) 2003-06-27 2017-11-07 Qualcomm Incorporated Local area network assisted positioning
US8971913B2 (en) 2003-06-27 2015-03-03 Qualcomm Incorporated Method and apparatus for wireless network hybrid positioning
US10895648B2 (en) 2003-06-27 2021-01-19 Qualcomm Incorporated Method and apparatus for wireless network hybrid positioning
US10849092B2 (en) 2003-06-27 2020-11-24 Qualcomm Incorporated Local area network assisted positioning
US9749876B2 (en) 2003-06-27 2017-08-29 Qualcomm Incorporated Local area network assisted positioning
US9810761B2 (en) 2003-06-27 2017-11-07 Qualcomm Incorporated Local area network assisted positioning
US9335419B2 (en) 2003-06-27 2016-05-10 Qualcomm Incorporated Wireless network hybrid positioning
USRE45808E1 (en) 2004-06-18 2015-11-17 Qualcomm Incorporated Method and apparatus for determining location of a base station using a plurality of mobile stations in a wireless mobile network
US9042917B2 (en) 2005-11-07 2015-05-26 Qualcomm Incorporated Positioning for WLANS and other wireless networks
US9226257B2 (en) 2006-11-04 2015-12-29 Qualcomm Incorporated Positioning for WLANs and other wireless networks
US10568062B2 (en) 2006-11-04 2020-02-18 Qualcomm Incorporated Positioning for WLANs and other wireless networks
US8504061B2 (en) 2010-04-07 2013-08-06 Apple Inc. Multi-tier geofence detection
US8868104B2 (en) 2010-04-07 2014-10-21 Apple Inc. Multi-tier geofence detection
US10064003B2 (en) 2010-11-09 2018-08-28 Apple Inc. Beacon-based geofencing
WO2012064860A1 (en) * 2010-11-09 2012-05-18 Apple Inc. Beacon-based geofencing
US8396485B2 (en) 2010-11-09 2013-03-12 Apple Inc. Beacon-based geofencing
US8996030B2 (en) 2010-11-09 2015-03-31 Apple Inc. Beacon-based geofencing
US9467803B2 (en) 2012-09-27 2016-10-11 Telefonaktiebolaget Lm Ericsson (Publ) Detecting multipath and determining positioning measurement uncertainty
EP2901737A4 (en) * 2012-09-27 2016-06-15 Ericsson Telefon Ab L M Detecting multipath and determining positioning measurement uncertainty

Also Published As

Publication number Publication date
EP1328822B1 (en) 2009-08-05
EP1328822A2 (en) 2003-07-23
IL155635A0 (en) 2003-11-23
JP4414136B2 (en) 2010-02-10
JP4550872B2 (en) 2010-09-22
CN1483148A (en) 2004-03-17
ATE438868T1 (en) 2009-08-15
EP2000820A2 (en) 2008-12-10
AU2002239736B2 (en) 2006-05-18
JP2008032730A (en) 2008-02-14
JP2004515767A (en) 2004-05-27
IL155635A (en) 2009-05-04
KR100829666B1 (en) 2008-05-16
JP2008032731A (en) 2008-02-14
CN101101326A (en) 2008-01-09
EP2000820A3 (en) 2009-01-14
AU3973602A (en) 2002-06-18
WO2002046788A3 (en) 2003-01-09
US6420999B1 (en) 2002-07-16
CN101187700A (en) 2008-05-28
KR20030045142A (en) 2003-06-09
DE60139493D1 (en) 2009-09-17
JP2008039786A (en) 2008-02-21
CA2426341A1 (en) 2002-06-13
CN100401093C (en) 2008-07-09
CN101101327A (en) 2008-01-09

Similar Documents

Publication Publication Date Title
US6420999B1 (en) Method and apparatus for determining an error estimate in a hybrid position determination system
AU2002239736A1 (en) Method and apparatus for determining an error estimate in a hybrid position determination system
KR100684541B1 (en) Method and apparatus for determining an algebraic solution to gps terrestrial hybrid location system equations
RU2327303C2 (en) Positioning of wireless communication terminal device in mixed positioning system
EP1552323B1 (en) Method and apparatus for improving radio location accuracy
CN102209906B (en) Method and system for localizing mobile communications terminals
US6570529B2 (en) Autonomous calibration of a wireless-global positioning system
US7646338B2 (en) Method of optimization of processing of location data in the presence of a plurality of satellite positioning constellations
CN1318752A (en) Acquisition of pilot phase offset delay parameter for radio terminal
AU2005291960A1 (en) A method for finding the location of a mobile terminal in a cellular radio system
KR20030062344A (en) Apparatus for processing of gps signals
EP1417508A2 (en) Method and apparatus for determining location using a coarse position estimate
KR100573951B1 (en) Method for correcting gps location using environmental information
Lee et al. Weighted Positioning Algorithm for Considering Geometrical Errors of Relay-based Positioning System
JP2002006027A (en) Method for obtaining reciprocating delay time (rtd) for ratio terminal of wireless network global positioning integrating (wgp) system

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A2

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NO NZ PH PL PT RO RU SD SE SG SI SK SL TJ TM TR TT TZ UA UG UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: A2

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2426341

Country of ref document: CA

Ref document number: 1020037005790

Country of ref document: KR

Ref document number: 628/CHENP/2003

Country of ref document: IN

WWE Wipo information: entry into national phase

Ref document number: 155635

Country of ref document: IL

Ref document number: 2002548471

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 2001987533

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2002239736

Country of ref document: AU

WWP Wipo information: published in national office

Ref document number: 1020037005790

Country of ref document: KR

WWE Wipo information: entry into national phase

Ref document number: 018213464

Country of ref document: CN

WWP Wipo information: published in national office

Ref document number: 2001987533

Country of ref document: EP

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642