US20030004640A1 - Method and apparatus for providing accurate position estimates in instances of severe dilution of precision - Google Patents

Method and apparatus for providing accurate position estimates in instances of severe dilution of precision Download PDF

Info

Publication number
US20030004640A1
US20030004640A1 US09/922,472 US92247201A US2003004640A1 US 20030004640 A1 US20030004640 A1 US 20030004640A1 US 92247201 A US92247201 A US 92247201A US 2003004640 A1 US2003004640 A1 US 2003004640A1
Authority
US
United States
Prior art keywords
position estimate
convergence
device position
vector
update
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
US09/922,472
Other versions
US6505122B1 (en
Inventor
Alkinoos Vayanos
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qualcomm Inc
Original Assignee
Qualcomm Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qualcomm Inc filed Critical Qualcomm Inc
Priority to US09/922,472 priority Critical patent/US6505122B1/en
Assigned to QUALCOMM INCORPORATED reassignment QUALCOMM INCORPORATED ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: BENCHETRIT, URI
Priority to CA2451290A priority patent/CA2451290C/en
Priority to EP02737535A priority patent/EP1399756B1/en
Priority to JP2003507578A priority patent/JP4180505B2/en
Priority to IL15927802A priority patent/IL159278A0/en
Priority to CNB028165624A priority patent/CN100420958C/en
Priority to PCT/US2002/019406 priority patent/WO2003001236A2/en
Priority to MYPI20022343A priority patent/MY130562A/en
Priority to ARP020102381A priority patent/AR034627A1/en
Priority to TW091113877A priority patent/TWI229743B/en
Publication of US20030004640A1 publication Critical patent/US20030004640A1/en
Publication of US6505122B1 publication Critical patent/US6505122B1/en
Application granted granted Critical
Priority to IL159278A priority patent/IL159278A/en
Priority to HK05108326.3A priority patent/HK1076512A1/en
Assigned to QUALCOMM INCORPORATED reassignment QUALCOMM INCORPORATED CORRECTIVE ASSIGNMENT TO CORRECT THE INVENTOR NAME URI BENCHETRIT INCORRECT, SHOULD BE ALKINOOS H. VAYANOS PREVIOUSLY RECORDED ON REEL 012398 FRAME 0712. ASSIGNOR(S) HEREBY CONFIRMS THE NEW ASSIGNMENT FOR APPLICATION NO. 09/922,472. Assignors: VAYANOS, ALKINOOS H.
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • 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/14Determining absolute distances from a plurality of spaced points of known location
    • 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/45Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/46Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being of a radio-wave signal type
    • 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

Definitions

  • the present invention relates generally to position determination, and more specifically to a method and apparatus for improving estimating the location of a device under unfavorable dilution of precision (DOP) conditions.
  • DOP unfavorable dilution of precision
  • a common means by which to locate a device is to determine the amount of time required for signals transmitted from multiple sources at known locations to reach the device.
  • One system that provides signals from a plurality of transmitters of known locations is the well-known Global Positioning Satellite (GPS) system. Satellites in the GPS system are placed in precise orbits according to a GPS master plan. The position of the GPS satellites can be identified by different sets of information transmitted by the satellites themselves.
  • GPS Global Positioning Satellite
  • Another system that provides signals from transmitters (i.e., base stations) of known earth-bound locations is a wireless (cellular) communication system.
  • Signals from satellites and/or base stations may be used to determine the location of a device.
  • the amount of time required for the signals to travel from the transmitters to the device may be measured and used to compute the distances (or ranges) between the transmitters and the device.
  • the signals themselves may further include information indicative of the locations of the transmitters.
  • the device can “trilaterate” its position. For example, when range measurements are available, a measurement surface may be determined for each transmitter, with the surface originating at the transmitter's location and having a curvature determined by the range from the transmitter to the device. The intersection of three or more such measurement surfaces would then be the estimated position of the device.
  • the ranges to the transmitters and their locations are conventionally provided to an algorithm that operates on the data to provide a position estimate for the device.
  • This algorithm is typically a least mean square (LMS) algorithm that performs a number of iterations to arrive at a final solution for the device's position estimate. If provided with “good” range and transmitter location data, the LMS solution asymptotically converges toward a target “optimal” solution with each iteration, and the LMS algorithm may be terminated when the LMS solution is within a particular tolerance of the target solution.
  • LMS least mean square
  • the LMS algorithm does not converge toward the target solution.
  • DOP diffusion of precision
  • the measurement surfaces for these transmitters may intersect at a small angle (i.e., the surfaces are close to parallel as oppose to being more slanted toward perpendicular).
  • a small error in the range measurements would result in a large error in the point of intersection, and thus a corresponding large error in the estimated location of the device.
  • the measurement surfaces do not intersect, which then results in bad geometry and large errors at the optimal position estimate for the device.
  • Earth-bound measurements may be made from signals transmitted from base stations of a wireless communication system and received at the device or signals received at the base stations from the device. The distances from the device to the base stations are shorter than those to the satellites, and the curvature of the earth-bound measurement surfaces is larger, which then increases the likelihood that the measurement surfaces do not intersect.
  • aspects of the invention provide techniques to derive an improved position estimate for a receiver device (i.e., a “Device Position Estimate”) in cases of large “dilution of precision” (DOP).
  • the Device Position Estimate, “Actual Measurement Vectors” and “Transmitter Position Estimates” are provided to an iterative algorithm (e.g., an LMS algorithm) that is initially operated in a normal manner.
  • an iterative algorithm e.g., an LMS algorithm
  • the solutions from the algorithm may not converge toward a target solution, but may instead oscillate (i.e., overshoot) around the target solution.
  • a method and apparatus is disclosed herein to determine when convergence is unlikely and to adjust the algorithm to increase the likelihood of convergence. Accordingly, the disclosed method and apparatus frequently increases the accuracy of the final solutions calculated in cases of large DOP.
  • a determination that the LMS algorithm is unlikely to converge may be made based on: (1) the number of iterations performed without determining that the algorithm has converged to a final solution (hereafter referred to as the “target solution”), (2) phase reversal in an “Update Vector” used to update the Device Position Estimate, or (3) a combination of these parameters.
  • One embodiment of the disclosed method determines the Transmitter Position Estimates for a number of transmitters (e.g., GPS satellites and/or base stations) determines an “Actual Measurement Vector” by making measurements, and calculates a “Calculated Measurement Vector” related to the distance between a current device and each transmitter (e.g., the pseudo-ranges to the transmitters that would result from the transmission of signals between the current device and each transmitter). In addition, a “Residual Measurement Error Vector” is calculated. The Residual Measurement Error Vector represents the difference between the Actual Measurement Vector and the Calculated Measurement Vector.
  • An Update Vector for a current Device Position Estimate is then computed based on the Residual Measurement Error Vector, the Transmitter Position Estimates, and a previous Device Position Estimate. A determination is made as to whether convergence of the Device Position Estimate toward a target solution is likely. If convergence is not likely, then the Update Vector is adjusted to increase the likelihood of convergence toward the target solution.
  • the target solution is defined as that solution that results in the minimum value for a metric calculated as a function of the errors that are included in the Residual Measurement Error Vector.
  • the criteria are normally dependent on the particular algorithm being used to determine the Device Position Estimate.
  • the current Device Position Estimate is updated based on the Update Vector. It should be noted that in one embodiment of the disclosed method and apparatus, the Update Vector is adjusted prior to being used to update the Device Position Estimate.
  • the Update Vector may be adjusted by: (1) reducing the magnitude of the Update Vector based on a scaling factor, and/or (2) limiting the magnitude of the Update Vector based on the magnitude of a Residual Measurement Error Vector as computed when a determination has been made that convergence is unlikely.
  • the scaling factor may be increased (e.g., doubled) with each detected event indicative of non-convergence (possibly except for the first event).
  • FIG. 1 is a simplified illustration of a system capable of implementing various aspects of the disclosed method and apparatus
  • FIG. 2 is a simplified block diagram of a receiver device in accordance with one embodiment of the disclosed method and apparatus
  • FIG. 3 is a flow diagram illustrating a process to estimate the location of the receiver device, in accordance with one embodiment of the disclosed method and apparatus;
  • FIG. 4A is a flow diagram illustrating a process performed to estimate the location of the receiver device using an LMS algorithm, in accordance with one embodiment of the disclosed method and apparatus.
  • FIG. 4B is a flow diagram of a process to detect non-convergence and adjust parameters of the LMS algorithm, in accordance with one embodiment of the disclosed method and apparatus.
  • FIG. 1 is a simplified illustration of a system capable of implementing various aspects of the disclosed method and apparatus.
  • a receiver device 110 to be located receives signals transmitted from a plurality of transmitters 130 a - 130 z .
  • Transmitters 130 may be any type of transmitter having locations that are known or that can be ascertained.
  • transmitters 130 include satellites of a Global Positioning Satellite (GPS) system, such as transmitters 130 a - 130 d in FIG. 1.
  • Transmitters 130 may additionally or exclusively include terrestrial ground-based transmitters, such as transmitters 130 y and 130 z , which may be base stations of a wireless communication system.
  • Receiver device 110 may thus receive signals from any combination of satellites and/or ground-base transmitters.
  • GPS Global Positioning Satellite
  • receiver device 110 may be any device capable of determining the arrival times of received signals with respect to a reference time.
  • receiver device 110 is a cellular telephone capable of receiving signals from the plurality of transmitters 130 .
  • receiver device 110 may be a computer terminal having a wireless modem, a stand-alone GPS receiver, a receiver capable of receiving signals from ground-based transmitters, or any other type of receiver.
  • the position of receiver device 110 may be estimated based on signals received at the device (e.g., such as those transmitted by GPS satellites and/or base stations) and the locations of the transmitters from which the received signals originated.
  • the position of the receiver device may also be estimated based, in part, on a signal transmitted by the device and received at one or more base stations.
  • any combination of measurements for signals received at or transmitted by the receiver device may be used to estimate the location of the device.
  • various aspects and embodiments of the disclosed method and apparatus are described for cases in which the signals are transmitted by the transmitters (e.g., satellites and/or base stations) and received at the receiver device.
  • Estimates of the position of receiver device 110 may be derived by the receiver device itself, a Position Determining Equipment (PDE) in the wireless communication system, a base station, or some other entity.
  • PDE Position Determining Equipment
  • the entity performing the position estimate is provided with a set of measurements (hereafter referred to as an “Actual Measurement Vector”) and the locations of the transmitters (hereafter referred to as “Transmitter Position Estimates”) or the means to determine these locations.
  • the Transmitter Position Estimates for GPS satellites may be ascertained by processing the signals transmitted by the satellites. Each satellite transmits “Almanac” information, which includes information regarding the coarse location of all satellites in the “constellation”. Each satellite further transmits “Ephemeris” information, which includes a higher accuracy version of its own orbit, as tracked and reported by tracking stations on earth.
  • the Transmitter Position Estimates for the base stations may be made known to the entity performing the position estimate for the receiver device.
  • the receiver device may include a database of the Transmitter Position Estimates for the base stations.
  • these Transmitter Position Estimates may be provided by the PDE, or the PDE may calculate the Device Position Estimate and may have the Transmitter Position Estimates for the base stations.
  • the Transmitter Position Estimates for the base stations may also be transmitted via messages.
  • Transmitters 130 can be used as reference points to determine the Device Position Estimate. By accurately measuring the distances to three transmitters 130 at known locations, the Device Position Estimate can be determined by “trilateration”. The distance to each transmitter can be determined by measuring the time required for a signal to travel from the transmitter to the receiver device. If the time the signal is transmitted from the transmitter is known (e.g., stamped into the signal), then the travel time of the signal can be determined by observing the time the signal is received by the receiver device 110 based on the internal clock of the receiver device 110 . Typically however, the exact amount of time between transmission and reception cannot be exactly determined because of offsets in the clocks at the transmitter 130 and the receiver device 110 .
  • a “pseudo-range” is typically computed based on the difference between a reference time and the time that the signal is received.
  • the reference time may be any time, as long as the reference time is common to all pseudo-range measurements being made or the pseudo-range measurements can be adjusted to compensate for differences in the reference times used.
  • Various algorithms can be used to determine the Device Position Estimate based on knowledge of the Transmitter Position Estimates and the Actual Measurement Vector.
  • One such algorithm is a least mean square (LMS) algorithm that performs a number of iterations to arrive at an accurate Device Position Estimate. Each iteration should result in a more accurate Device Position Estimate. However, if the algorithm is not converging, this will not be the case.
  • LMS least mean square
  • DOP deciution of precision
  • the LMS algorithm asymptotically converges toward a target solution. However, assuming a large DOP around the target solution, the LMS algorithm may not converge if executed in the normal manner, resulting in an inaccurate solution.
  • the disclosed method and apparatus determines a Device Position Estimate in cases of large DOP.
  • an “Actual Measurement Vector” e.g., a set of measurements, such as measured pseudo-ranges
  • Transmitter Position Estimates are provided to an (e.g., LMS) algorithm that initially operates in a normal manner.
  • the solutions from the algorithm may not converge toward a target solution, but may instead overshoot the target solution, causing the algorithm to “oscillate”.
  • the presently disclosed method and apparatus detects cases in which convergence is unlikely and adjusts the operation of the algorithm to increase the likelihood of convergence toward the target solution. This results in a greater likelihood of convergence in cases of large DOP, and therefore a more accurate solution.
  • FIG. 2 is a simplified block diagram of receiver device 110 in accordance with one embodiment of the disclosed method and apparatus.
  • Receiver device 110 may be a component of a terminal in a wireless communication system, such as: (1) a cellular telephone, (2) a computer with a wireless modem, (3) a personal digital assistant, or (4) any other such portable device with a communication capability.
  • receiver device 110 may be a stand-alone position determining unit, such as a stand-alone GPS receiver, which has no communication capability, or a limited communication capability that is only available to assist the position location function of the device.
  • Receiver device 110 may be designed with the capability to process signals from satellites (such as GPS satellites or other position location satellites), terrestrial ground-base transmitters (e.g., base stations, Loran stations, etc.), or both.
  • receiver device 110 includes an antenna 210 , a terrestrial receiver unit 212 a , a GPS receiver unit 212 b , a processing unit 216 , a memory unit 218 , and a controller 220 .
  • Antenna 210 receives signals from transmitters 130 and couples the received signals to terrestrial and GPS receiver units 212 a and 212 b .
  • Terrestrial receiver unit 212 a includes front-end circuitry (e.g., radio frequency (RF) processing circuitry and/or other receiving circuitry) that processes signals transmitted from base stations to derive information used for position determination.
  • GPS receiver unit 212 b includes front-end circuitry that processes signals transmitted from GPS satellites to derive information used for position determination.
  • Such information from receiver units 212 a and 212 b may include timing information, the identities and locations of the transmitters whose signals are received, and possibly other information.
  • Processing unit 216 may initially determine an initial Device Position Estimate and receive from the GPS receiver unit 212 b an Actual Measurement Vector. Processing unit 216 then executes the LMS algorithm, which operates on the computed results to arrive at a final Device Position Estimate.
  • Memory unit 218 stores various data required for determining position.
  • memory unit 218 may store Almanac and/or Ephemeris information that indicates where each of a number of GPS satellites is located in the sky at any given time. Both the Almanac and Ephemeris are valid for a limited amount of time.
  • the Almanac information is considered to be accurate to approximately 3 kilometers for approximately one week from the time the Almanac is transmitted.
  • the Ephemeris provides information regarding the satellite orbit with an accuracy of approximately 1 meter for approximately 2 hours. The error in both the Almanac and Ephemeris grows as the information ages.
  • the satellite locations may be derived from the Almanac and Ephemeris or from other information obtained by GPS receiver unit 212 b from signals transmitted from the satellites or by terrestrial receiver unit 212 a from signals transmitted via the wireless link.
  • Controller 220 may receive signaling and/or intermediate results from processing unit 216 and may direct the operation of the processing unit. For example, controller 220 may select the particular types of signals to be processed (e.g., satellite, terrestrial, or both), the particular algorithm to be used (if more than one is available), the parameter values to be used for the selected algorithm, and so on, as described in further detail below.
  • controller 220 may select the particular types of signals to be processed (e.g., satellite, terrestrial, or both), the particular algorithm to be used (if more than one is available), the parameter values to be used for the selected algorithm, and so on, as described in further detail below.
  • receiver device 110 may communicate with a PDE, which may assist in determining the position estimate for the device.
  • the PDE may perform the computations to derive the receiver device position estimate, or may provide certain information used to determine the position estimate (e.g., the location of the GPS satellites and/or base stations).
  • An iterative algorithm may be used to determine the Device Position Estimate based on (1) Transmitter Position Estimates (e.g., information regarding the locations of GPS satellites and/or base stations) and (2) an Actual Measurement Vector (e.g., a set of pseudo-range measurements to the transmitters).
  • Transmitter Position Estimates e.g., information regarding the locations of GPS satellites and/or base stations
  • an Actual Measurement Vector e.g., a set of pseudo-range measurements to the transmitters.
  • the transmitter locations may be ascertained based on current Almanac and/or Ephemeris information.
  • the LMS iterative algorithm is specifically described below. However, other iterative algorithms may also be used to determine the Device Position Estimate in accordance with the disclosed method and apparatus.
  • Each Transmitter Position Estimate is an estimate of the location of one transmitter.
  • the Transmitter Position Estimate, ⁇ overscore (s) ⁇ i for the i-th transmitter is expressed as:
  • the b i is the clock bias expressed in the same dimension as x, y, and z and further includes the effects of the reference time.
  • ⁇ overscore (u) ⁇ k [x u,k y u,k z u,k b u,k ]
  • an initial Device Position Estimate ⁇ overscore (u) ⁇ 0
  • the value of the initial Device Position Estimate may be selected arbitrarily or based upon current or past knowledge about the position of the receiver device.
  • the initial Device Position Estimate, ⁇ overscore (u) ⁇ 0 may be: (1) the last Device Position Estimate (if one is available), (2) the center of the earth, (3) the center of the cell of which the receiver device is receiving the signal, or (4) any other location that provides an initial value to the LMS algorithm. The more accurate the initial Device Position Estimate, the faster the algorithm will convergence toward the target solution.
  • Each value of the Calculated Measurement Vector is a computed measurement associated with one of the transmitters used for position determination.
  • the values of the Calculated Measurement Vector are calculated from: (1) the Transmitter Position Estimates, ⁇ overscore (s) ⁇ i , (2) the Device Position Estimate, ⁇ overscore (u) ⁇ k , resulting from the k-th iteration of the iterative algorithm. Accordingly, it will be understood that there will be one Calculated Measurement Vector, ⁇ overscore ( ⁇ ) ⁇ k , associated with each iteration of the iterative algorithm.
  • ⁇ overscore (s) ⁇ i and ⁇ overscore (u) ⁇ k are not made, but are calculated for each iteration based on the values of ⁇ overscore (s) ⁇ i and ⁇ overscore (u) ⁇ k .
  • an Actual Measurement Vector, ⁇ overscore ( ⁇ ) ⁇ a is derived from actually measuring the received signals.
  • the elements of the Calculated Measurement Vector, ⁇ overscore ( ⁇ ) ⁇ k may be any type of measurement that is a function of the location of the receiver device and the clock bias associated with the transmitters.
  • the elements of the Calculated Measurement Vector may relate to ranges, pseudo-ranges as determined at the receiver device (i.e., r ⁇ b, where r is the range and b is the clock bias), reverse pseudo-ranges as determined at the base stations based on signals transmitted from the receiver device (i.e., r+b), range differences as determined based on “time difference of arrival” (TDOA), or a combination thereof.
  • the elements of the Calculated Measurement Vector may be computed by any one of a number of known algorithms.
  • the elements of the Calculated Measurement Vector, ⁇ overscore ( ⁇ ) ⁇ k are pseudo-ranges.
  • One way to compute the elements of the Calculated Measurement Vector is:
  • ⁇ k,i ⁇ square root ⁇ square root over (( x u,k ⁇ x i ) 2 +( y u,k ⁇ y i ) 2 +( z u,k ⁇ z i ) 2 ) ⁇ ( b u,k ⁇ b i ) Eq (1)
  • each element, ⁇ k,i , of the Calculated Measurement Vector, ⁇ overscore ( ⁇ ) ⁇ k is the pseudo-range between the Device Position Estimate, ⁇ overscore (u) ⁇ k , for the k-th iteration and the Transmitter Position Estimate, ⁇ overscore (s) ⁇ i , for the i-th transmitter.
  • a covariance matrix R of the measurements is defined as:
  • ⁇ overscore ( ⁇ ) ⁇ a [ ⁇ a,1 ⁇ a,2 . . . ⁇ a,n ] T
  • ⁇ overscore ( ⁇ ) ⁇ a is the Actual Measurement Vector and ⁇ a,n is the measurement element (e.g., pseudo-range measurement) attained from the signal received from the n-th transmitter.
  • E is the statistical expectation.
  • the superscript T denotes the complex conjugate transpose.
  • the covariance matrix R is typically estimated based on the particular design of the receiver device and the system (e.g., the signal bandwidth or the received signal power).
  • the Residual Measurement Error Vector, ⁇ overscore (e) ⁇ k is equal to the difference between the Actual Measurement Vector, ⁇ overscore ( ⁇ ) ⁇ a , and the Calculated Measurement Vector, ⁇ overscore ( ⁇ ) ⁇ k . It should be noted that one Residual Measurement Error Vector is calculated for each iteration of the algorithm.
  • the Residual Measurement Error Vector, ⁇ overscore (e) ⁇ k is calculated from the Device Position Estimate, ⁇ overscore (u) ⁇ k , determined on the k-th iteration.
  • the next Residual Measurement Error Vector, ⁇ overscore (e) ⁇ k+1 is then newly calculated for the next iteration from the Device Position Estimate, ⁇ overscore (u) ⁇ k+1 .
  • An Update Vector, ⁇ overscore (u) ⁇ k is computed for each iteration based on the geometry matrix H, the covariance matrix R, and the Residual Measurement Error Vector, ⁇ overscore (e) ⁇ k , for that iteration.
  • the Update Vector is defined as:
  • ⁇ ⁇ overscore (u) ⁇ k ( H T ⁇ R ⁇ 1 ⁇ H ) ⁇ 1 ⁇ H T ⁇ R ⁇ 1 ⁇ overscore (e) ⁇ k Eq (5)
  • This Update Vector represents an estimated error between the current Device Position Estimate, ⁇ overscore (u) ⁇ k , and an optimal position estimate.
  • an updated Device Position Estimate is expressed as:
  • Equations (1) through (6) comprise the computations for one iteration of the LMS algorithm. A number of iterations may be performed to derive a more accurate Device Position Estimate.
  • the LMS algorithm will most likely converge after a relatively small number of iterations. Accordingly, the Device Position Estimates provided by the LMS algorithm, ⁇ overscore (u) ⁇ k , asymptotically approach the target solution.
  • the target solution is defined as that solution that results in the minimum value for a metric calculated for the target solution.
  • the metric is a function of the errors that are included in the Residual Measurement Error Vector.
  • the disclosed method and apparatus (1) determines when the algorithm is unlikely to convergence (e.g., due to oscillation of the LMS solutions around the target solution) and (2) adjusts the algorithm to increase the likelihood of convergence.
  • the algorithm typically will convergence to the target solution after a predetermined number of iterations. This number is normally small (e.g., three). However, the typical number of iterations may be greater if the accuracy of the initial Device Position Estimate is poor. If the algorithm requires more iterations than what is typical, then there is less likelihood of convergence and greater likelihood of oscillation around the target solution.
  • Such oscillation can be identified by an inversion in the direction of two consecutive Update Vectors (i.e., 180° change in phase of the Update Vector). Whether the LMS algorithm is likely to converge can be determined by (1) the number of iterations performed without reaching the target solution and (2) changes in the phase of the Update Vector.
  • the magnitude of the Update Vector is reduced to (1) avoid significant overshoots and (2) adjust the rate at which the position estimate is updated, both of which help to ensure convergence.
  • FIG. 3 is a flow diagram illustrating a process 300 for determining the Device Position Estimate in accordance with one embodiment of the disclosed method and apparatus.
  • the information required to determine the Device Position Estimate is initially received, at step 312 .
  • This information typically includes the Transmitter Position Estimates and information used to determine the calculated Measurement Vector (e.g., pseudo-range measurements).
  • An iterative algorithm e.g., an LMS algorithm is then executed to determine the Device Position Estimate based on the received information, at step 314 .
  • a final Device Position Estimate is determined.
  • the algorithm may be terminated for any one of a number of reasons, such as the maximum number of iterations have been performed, the algorithm has sufficiently converged, etc.
  • a final Device Position Estimate is determined if the magnitude of the Update Vector is smaller than a predetermined threshold, indicating that the algorithm has sufficiently converged. If the algorithm has not terminated, then the process of FIG. 3 continues at step 314 .
  • Steps 314 through 322 may be viewed as comprising one iteration of the iterative algorithm, and any number of iterations may be performed until sufficiently accurate Device Position Estimate is reached.
  • a further discussion of how the likelihood of convergence is determined and the parameters are adjusted is provided below.
  • the disclosed method and apparatus is described below using an LMS algorithm as one example. However, it will be understood by those skilled in the art that other iterative algorithms could be used to calculate an estimate of the location of a receiver device.
  • FIG. 4A is a flow diagram illustrating a process 400 performed to determine the Device Position Estimate using the LMS algorithm in accordance with one embodiment of the disclosed method and apparatus.
  • the information required to determine the Device Position Estimate is first received at step 412 .
  • This information typically includes: (1) the Transmitter Position Estimates, ⁇ overscore (s) ⁇ t , (2) the initial Device Position Estimate, ⁇ overscore (u) ⁇ 0 ; and (3) the Actual Measurement Vector, ⁇ overscore ( 92 ) ⁇ a .
  • Variables used for the LMS algorithm are then initialized, at step 414 .
  • UV_modify False
  • k is the iteration index
  • g is a scaling factor used to adjust the magnitude of the Update Vector
  • ⁇ overscore (u) ⁇ k indicates whether or not to adjust the Update Vector.
  • Execution of the LMS algorithm then commences.
  • the elements of the Actual Measurement Vector represent pseudo-range measurements.
  • An initial Device Position Estimate is made.
  • a Calculated Measurement Vector is then determined from the initial Device Position Estimate.
  • the Calculated Measurement Vector, ⁇ overscore ( ⁇ ) ⁇ k is computed as shown in equation (1), at step 416 .
  • the geometry matrix H is updated based on the current Device Position Estimate, ⁇ overscore (u) ⁇ k , and the Transmitter Position Estimate, ⁇ overscore (s) ⁇ i , as shown in equation (2), at step 418 .
  • the measurement covariance matrix R is also derived, at step 418 .
  • the Residual Measurement Error Vector, ⁇ overscore (e) ⁇ k is then computed based on the Actual Measurement Vector and the Calculated Measurement Vector as shown in equation (4), at step 420 .
  • the Update Vector, ⁇ overscore (u) ⁇ k is next computed based on the matrices H and R and the Residual Measurement Error Vector, ⁇ overscore (e) ⁇ k , as shown in equation (5), at step 422 .
  • the new Device Position Estimate, ⁇ overscore (u) ⁇ k+1 , for the receiver device is computed based on the current Device Position Estimate, ⁇ overscore (u) ⁇ k , and the Update Vector, ⁇ overscore (u) ⁇ k , as shown in equation (6).
  • Termination of the algorithm may be triggered based on one or more criteria, such as: (1) the magnitude of the Update Vector, (2) the number of iterations already performed, or a combination thereof.
  • the algorithm is terminated if the magnitude of the Update Vector is smaller than a particular update threshold (i.e.,
  • the current Device Position Estimate is considered to be sufficiently accurate, at step 434 , and the process terminates.
  • the Device Position Estimate will in fact be accurate if the algorithm has converged. However the Device Position Estimate will not be as accurate as desired if the algorithm has not converged, but rather was terminated due to some other criterion (e.g., the number of iterations exceed the maximum number of iterations). If a termination criterion has not been satisfied, the iteration index k is incremented by one, at step 432 , and the process returns to step 416 to perform another iteration.
  • Steps 416 through 428 comprise one iteration of the LMS algorithm. Any number of iterations may be executed until the termination criterion or criteria are met.
  • FIG. 4B is a flow diagram of a process to determine whether convergence is likely. Also shown is a process for adjusting parameters of the LMS algorithm in accordance with one embodiment of the disclosed method and apparatus. This process may be used for step 426 in FIG. 4A.
  • a determination that convergence of the LMS algorithm is unlikely is in part based on the number of iterations performed without achieving convergence. Convergence is determined by attaining a sufficiently small norm of the Update Vector. As noted above, in normal circumstances in a wireless-assisted GPS position determination system, convergence typically occurs within a small number of iterations. Thus, a determination is made whether the number of iterations performed by the LMS algorithm for this collection of input data exceeds a particular threshold “k_threshold” (i.e., whether k>k_threshold), at step 448 .
  • k_threshold i.e., whether k>k_threshold
  • k_threshold is typically selected based on the number of iterations expected to be performed to converge for a typical application (which may be dependent on the accuracy of the initial device position estimate). Alternatively, any value may also be selected that can provide good performance. As a specific example, k threshold may be set to three. If the number of iterations already performed is greater than k_threshold, the likelihood that convergence by the algorithm will occur is assumed to be low and the process proceeds to step 450 where another test for convergence is performed. Otherwise, if the number of iterations k is less than or equal to k_threshold, a determination is made that the algorithm may yet converge.
  • a determination that the LMS algorithm is unlikely to convergence is made by examining the Update Vectors. Convergence is determined to be unlikely when the Device Position Estimates repeatedly overshoot the desired Device Position Estimate. Overshooting by the LMS algorithm causes the set of Device Position Estimates resulting from each iteration to “oscillate” around the desired Device Position Estimate. This oscillation can typically be seen by changes in the direction of the Update Vector essentially along the line that is coincident with the direction of the largest DOP. When such phase reversals occur, the phase difference ⁇ between two consecutive Update Vectors is close to 180°. This results in the cosine of the phase difference ⁇ being close to minus one. In one embodiment, reversal in the direction of the Update Vector is detected if the cosine of ⁇ is less than a particular negative cos_threshold. As a specific example, cos_threshold may be selected as ⁇ 0.95.
  • step 450 the cosine of the phase difference ⁇ of two consecutive Update Vectors is computed.
  • a variable UV_modify is set to True.
  • the factor g is multiplied by a positive factor greater than unity (e.g., two). This larger value for the factor g is then used to effectively divide the magnitude of the Update Vector by a corresponding amount.
  • the adjusted Update Vector may be smaller than optimal, and under normal operating conditions this smaller Update Vector would lead to a slower convergence rate. However, for unstable scenarios, the smaller Update Vector ensures that, given a sufficient number of iterations, the LMS algorithm has a greater likelihood of convergence.
  • the number of phase reversals does not correspond to the amount the Update Vector is reduced. Rather the Update Vector is reduced by a predetermined amount each time a reduction is made.
  • a parameter of the LMS algorithm may be continually (i.e., “adaptively”) adjusted based on the monitored performance of the algorithm.
  • such an “adaptive” adjustment is achieved based on the values of the variable d and the factor g.
  • the Update Vector is adjusted by the variable d and factor g, as described below.
  • the adjusted Update Vector retains the phase of the original Update Vector, but its norm is set to the magnitude of the Residual Measurement Error Vector at the time the last phase reversal was detected (as captured by the variable d) and further reduced by the factor g.
  • the adjusting ensures that the magnitude of the Update Vector is smaller each time the phase reversal is detected, since the amount of overshoot is less each time. In this manner, the step size can be made dependent on the monitored performance of the LMS algorithm. Moreover, by increasing (e.g., doubling) the factor g each time a phase reversal is detected, the size of the Update Vector is correspondingly reduced.
  • a phase reversal of the Update Vector was not detected at step 452 , a determination is made at to whether the variable UV_modify has been set to True, at step 464 . If the answer is no, indicating that the LMS may be taking longer to converge, but no oscillation in the Update Vector has been detected, then the process terminates. Otherwise, if UV-modify has already been set to True by a previous detected phase reversal, then the Update Vector is adjusted based on the current values for the variable d and factor g, at step 462 . In the embodiment shown in FIG. 4B, the Update Vector is adjusted by the variable d and factor g for each iteration, once UV-modify has been set to True, regardless of whether or not phase reversal of the Update Vector has been detected.
  • FIGS. 4A and 4B illustrate examples of the disclosed method for performing the LMS algorithm. However, these examples may also be applied to other algorithms used to determine the Device Position Estimate. For clarity, FIGS. 4A and 4B also describe a specific implementation for determining that convergence is unlikely based on the number of iterations and the Update Vectors. Also described is an example of one method by which the LMS algorithm can be adjusted to improve the likelihood of convergence. It should be noted that these methods are merely presented as examples of ways to determine whether convergence is likely and to adjust the parameters of the LMS algorithm. Other methods are clearly possible and would be understood by those skilled in the art from the examples provided herein.
  • the elements of receiver device 110 may be implemented as a digital signal processor (DSP), an application specific integrated circuit (ASIC), a processor, a microprocessor, a controller, a microcontroller, a field programmable gate array (FPGA), a programmable logic device, other electronic unit, or any combination thereof designed to perform the functions described herein.
  • DSP digital signal processor
  • ASIC application specific integrated circuit
  • FPGA field programmable gate array
  • Certain aspects of the disclosed method and apparatus may be implemented in hardware, software, or a combination of both.
  • the processing to determine the Device Position Estimate e.g., the LMS algorithm
  • the processing to determine the Device Position Estimate may be performed based on program codes stored within memory unit 218 and executed by a processor (processing unit 216 or controller 220 in FIG. 2).
  • such processing may be performed by dedicated hardware or some combination of hardware and software.

Abstract

Techniques for providing an improved position estimate for a device in instances of relatively large dilution of precision. In a method, position estimates for a number of transmitters (e.g., GPS satellites and/or base stations) and a set of initial measurements (e.g., pseudo-ranges) are initially received. An Update Vector for a current device position estimate is computed. Non-convergence of the current device position estimate toward a target position estimate is detected and, if true, the Update Vector is adjusted (e.g., reduced in magnitude by a scaling factor) to increase the likelihood of convergence to the target position estimate. The current device position estimate is then updated based on the (possibly adjusted) Update Vector. Non-convergence may be detected based on the number of iterations performed without reaching the target position estimate and/or phase reversal in two consecutive Update Vectors. The scaling factor may be increased with each detected event indicative of non-convergence.

Description

    RELATED APPLICATIONS
  • This application claims priority to U.S. Provisional Application No. 60/300,960, filed Jun. 25, 2001.[0001]
  • BACKGROUND
  • 1. Field [0002]
  • The present invention relates generally to position determination, and more specifically to a method and apparatus for improving estimating the location of a device under unfavorable dilution of precision (DOP) conditions. [0003]
  • 2. Background [0004]
  • A common means by which to locate a device is to determine the amount of time required for signals transmitted from multiple sources at known locations to reach the device. One system that provides signals from a plurality of transmitters of known locations is the well-known Global Positioning Satellite (GPS) system. Satellites in the GPS system are placed in precise orbits according to a GPS master plan. The position of the GPS satellites can be identified by different sets of information transmitted by the satellites themselves. Another system that provides signals from transmitters (i.e., base stations) of known earth-bound locations is a wireless (cellular) communication system. [0005]
  • Signals from satellites and/or base stations may be used to determine the location of a device. By receiving and processing the signals transmitted from these transmitters, the amount of time required for the signals to travel from the transmitters to the device may be measured and used to compute the distances (or ranges) between the transmitters and the device. The signals themselves may further include information indicative of the locations of the transmitters. By accurately measuring the distances from three or more transmitters at known locations, the device can “trilaterate” its position. For example, when range measurements are available, a measurement surface may be determined for each transmitter, with the surface originating at the transmitter's location and having a curvature determined by the range from the transmitter to the device. The intersection of three or more such measurement surfaces would then be the estimated position of the device. [0006]
  • The ranges to the transmitters and their locations are conventionally provided to an algorithm that operates on the data to provide a position estimate for the device. This algorithm is typically a least mean square (LMS) algorithm that performs a number of iterations to arrive at a final solution for the device's position estimate. If provided with “good” range and transmitter location data, the LMS solution asymptotically converges toward a target “optimal” solution with each iteration, and the LMS algorithm may be terminated when the LMS solution is within a particular tolerance of the target solution. [0007]
  • In certain circumstances, the LMS algorithm does not converge toward the target solution. These circumstances typically arise due to a large “dilution of precision” (DOP) around the device location, which may result from having inaccurate range data and/or poor geometry for the transmitter and device locations. As an example, if the transmitters are located in close proximity relative to the distance to the device, then the measurement surfaces for these transmitters may intersect at a small angle (i.e., the surfaces are close to parallel as oppose to being more slanted toward perpendicular). In this case, a small error in the range measurements would result in a large error in the point of intersection, and thus a corresponding large error in the estimated location of the device. In one especially case of very large DOP, the measurement surfaces do not intersect, which then results in bad geometry and large errors at the optimal position estimate for the device. [0008]
  • The likelihood of obtaining bad geometry increases when earth-bound measurements are used (exclusively or in conjunction with satellite measurements) to estimate the position of the device. Earth-bound measurements may be made from signals transmitted from base stations of a wireless communication system and received at the device or signals received at the base stations from the device. The distances from the device to the base stations are shorter than those to the satellites, and the curvature of the earth-bound measurement surfaces is larger, which then increases the likelihood that the measurement surfaces do not intersect. [0009]
  • Conventionally, solutions obtained in cases of bad geometry are discarded. This may be due to a combination of non-convergence by the LMS algorithm and/or reduced likelihood of obtaining the desired accuracy in the LMS solutions. However, in many instances, such as 911 emergency services, these solutions are beneficial and should be reported if they can be obtained. [0010]
  • There is therefore a need in the art for techniques to provide position estimate with improved likelihood of having the desired accuracy in cases of large DOP. [0011]
  • SUMMARY
  • Aspects of the invention provide techniques to derive an improved position estimate for a receiver device (i.e., a “Device Position Estimate”) in cases of large “dilution of precision” (DOP). The Device Position Estimate, “Actual Measurement Vectors” and “Transmitter Position Estimates” are provided to an iterative algorithm (e.g., an LMS algorithm) that is initially operated in a normal manner. In cases of large DOP, the solutions from the algorithm may not converge toward a target solution, but may instead oscillate (i.e., overshoot) around the target solution. A method and apparatus is disclosed herein to determine when convergence is unlikely and to adjust the algorithm to increase the likelihood of convergence. Accordingly, the disclosed method and apparatus frequently increases the accuracy of the final solutions calculated in cases of large DOP. [0012]
  • A determination that the LMS algorithm is unlikely to converge may be made based on: (1) the number of iterations performed without determining that the algorithm has converged to a final solution (hereafter referred to as the “target solution”), (2) phase reversal in an “Update Vector” used to update the Device Position Estimate, or (3) a combination of these parameters. [0013]
  • One embodiment of the disclosed method determines the Transmitter Position Estimates for a number of transmitters (e.g., GPS satellites and/or base stations) determines an “Actual Measurement Vector” by making measurements, and calculates a “Calculated Measurement Vector” related to the distance between a current device and each transmitter (e.g., the pseudo-ranges to the transmitters that would result from the transmission of signals between the current device and each transmitter). In addition, a “Residual Measurement Error Vector” is calculated. The Residual Measurement Error Vector represents the difference between the Actual Measurement Vector and the Calculated Measurement Vector. An Update Vector for a current Device Position Estimate is then computed based on the Residual Measurement Error Vector, the Transmitter Position Estimates, and a previous Device Position Estimate. A determination is made as to whether convergence of the Device Position Estimate toward a target solution is likely. If convergence is not likely, then the Update Vector is adjusted to increase the likelihood of convergence toward the target solution. [0014]
  • The target solution is defined as that solution that results in the minimum value for a metric calculated as a function of the errors that are included in the Residual Measurement Error Vector. [0015]
  • The criteria are normally dependent on the particular algorithm being used to determine the Device Position Estimate. The current Device Position Estimate is updated based on the Update Vector. It should be noted that in one embodiment of the disclosed method and apparatus, the Update Vector is adjusted prior to being used to update the Device Position Estimate. [0016]
  • The Update Vector may be adjusted by: (1) reducing the magnitude of the Update Vector based on a scaling factor, and/or (2) limiting the magnitude of the Update Vector based on the magnitude of a Residual Measurement Error Vector as computed when a determination has been made that convergence is unlikely. The scaling factor may be increased (e.g., doubled) with each detected event indicative of non-convergence (possibly except for the first event).[0017]
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • The features, nature, and advantages of the disclosed method and apparatus will become more apparent from the detailed description set forth below when taken in conjunction with the drawings in which like reference characters identify correspondingly throughout and wherein: [0018]
  • FIG. 1 is a simplified illustration of a system capable of implementing various aspects of the disclosed method and apparatus; [0019]
  • FIG. 2 is a simplified block diagram of a receiver device in accordance with one embodiment of the disclosed method and apparatus; [0020]
  • FIG. 3 is a flow diagram illustrating a process to estimate the location of the receiver device, in accordance with one embodiment of the disclosed method and apparatus; [0021]
  • FIG. 4A is a flow diagram illustrating a process performed to estimate the location of the receiver device using an LMS algorithm, in accordance with one embodiment of the disclosed method and apparatus; and [0022]
  • FIG. 4B is a flow diagram of a process to detect non-convergence and adjust parameters of the LMS algorithm, in accordance with one embodiment of the disclosed method and apparatus.[0023]
  • DETAILED DESCRIPTION
  • FIG. 1 is a simplified illustration of a system capable of implementing various aspects of the disclosed method and apparatus. A [0024] receiver device 110 to be located receives signals transmitted from a plurality of transmitters 130 a-130 z. Transmitters 130 may be any type of transmitter having locations that are known or that can be ascertained. In one embodiment, transmitters 130 include satellites of a Global Positioning Satellite (GPS) system, such as transmitters 130 a-130 d in FIG. 1. Transmitters 130 may additionally or exclusively include terrestrial ground-based transmitters, such as transmitters 130 y and 130 z, which may be base stations of a wireless communication system. Receiver device 110 may thus receive signals from any combination of satellites and/or ground-base transmitters.
  • In general, [0025] receiver device 110 may be any device capable of determining the arrival times of received signals with respect to a reference time. In one embodiment, receiver device 110 is a cellular telephone capable of receiving signals from the plurality of transmitters 130. In other embodiments, receiver device 110 may be a computer terminal having a wireless modem, a stand-alone GPS receiver, a receiver capable of receiving signals from ground-based transmitters, or any other type of receiver.
  • The position of [0026] receiver device 110 may be estimated based on signals received at the device (e.g., such as those transmitted by GPS satellites and/or base stations) and the locations of the transmitters from which the received signals originated. The position of the receiver device may also be estimated based, in part, on a signal transmitted by the device and received at one or more base stations. Thus, any combination of measurements for signals received at or transmitted by the receiver device may be used to estimate the location of the device. For clarity, various aspects and embodiments of the disclosed method and apparatus are described for cases in which the signals are transmitted by the transmitters (e.g., satellites and/or base stations) and received at the receiver device.
  • Estimates of the position of receiver device [0027] 110 (hereafter referred to as “Device Position Estimates”) may be derived by the receiver device itself, a Position Determining Equipment (PDE) in the wireless communication system, a base station, or some other entity. The entity performing the position estimate is provided with a set of measurements (hereafter referred to as an “Actual Measurement Vector”) and the locations of the transmitters (hereafter referred to as “Transmitter Position Estimates”) or the means to determine these locations.
  • The Transmitter Position Estimates for GPS satellites may be ascertained by processing the signals transmitted by the satellites. Each satellite transmits “Almanac” information, which includes information regarding the coarse location of all satellites in the “constellation”. Each satellite further transmits “Ephemeris” information, which includes a higher accuracy version of its own orbit, as tracked and reported by tracking stations on earth. The Transmitter Position Estimates for the base stations may be made known to the entity performing the position estimate for the receiver device. For example, the receiver device may include a database of the Transmitter Position Estimates for the base stations. Alternatively, these Transmitter Position Estimates may be provided by the PDE, or the PDE may calculate the Device Position Estimate and may have the Transmitter Position Estimates for the base stations. The Transmitter Position Estimates for the base stations may also be transmitted via messages. [0028]
  • Transmitters [0029] 130 can be used as reference points to determine the Device Position Estimate. By accurately measuring the distances to three transmitters 130 at known locations, the Device Position Estimate can be determined by “trilateration”. The distance to each transmitter can be determined by measuring the time required for a signal to travel from the transmitter to the receiver device. If the time the signal is transmitted from the transmitter is known (e.g., stamped into the signal), then the travel time of the signal can be determined by observing the time the signal is received by the receiver device 110 based on the internal clock of the receiver device 110. Typically however, the exact amount of time between transmission and reception cannot be exactly determined because of offsets in the clocks at the transmitter 130 and the receiver device 110. Thus, a “pseudo-range” is typically computed based on the difference between a reference time and the time that the signal is received. The reference time may be any time, as long as the reference time is common to all pseudo-range measurements being made or the pseudo-range measurements can be adjusted to compensate for differences in the reference times used.
  • Various algorithms can be used to determine the Device Position Estimate based on knowledge of the Transmitter Position Estimates and the Actual Measurement Vector. One such algorithm is a least mean square (LMS) algorithm that performs a number of iterations to arrive at an accurate Device Position Estimate. Each iteration should result in a more accurate Device Position Estimate. However, if the algorithm is not converging, this will not be the case. If provided with accurate range and Transmitter Position Estimates and assuming a relatively small “dilution of precision” (DOP), the LMS algorithm asymptotically converges toward a target solution. However, assuming a large DOP around the target solution, the LMS algorithm may not converge if executed in the normal manner, resulting in an inaccurate solution. [0030]
  • The disclosed method and apparatus determines a Device Position Estimate in cases of large DOP. To determine the Device Position Estimate, an “Actual Measurement Vector” (e.g., a set of measurements, such as measured pseudo-ranges) and Transmitter Position Estimates are provided to an (e.g., LMS) algorithm that initially operates in a normal manner. In cases of large DOP, the solutions from the algorithm may not converge toward a target solution, but may instead overshoot the target solution, causing the algorithm to “oscillate”. The presently disclosed method and apparatus detects cases in which convergence is unlikely and adjusts the operation of the algorithm to increase the likelihood of convergence toward the target solution. This results in a greater likelihood of convergence in cases of large DOP, and therefore a more accurate solution. [0031]
  • FIG. 2 is a simplified block diagram of [0032] receiver device 110 in accordance with one embodiment of the disclosed method and apparatus. Receiver device 110 may be a component of a terminal in a wireless communication system, such as: (1) a cellular telephone, (2) a computer with a wireless modem, (3) a personal digital assistant, or (4) any other such portable device with a communication capability. Alternatively, receiver device 110 may be a stand-alone position determining unit, such as a stand-alone GPS receiver, which has no communication capability, or a limited communication capability that is only available to assist the position location function of the device.
  • [0033] Receiver device 110 may be designed with the capability to process signals from satellites (such as GPS satellites or other position location satellites), terrestrial ground-base transmitters (e.g., base stations, Loran stations, etc.), or both. In the embodiment shown in FIG. 2, receiver device 110 includes an antenna 210, a terrestrial receiver unit 212 a, a GPS receiver unit 212 b, a processing unit 216, a memory unit 218, and a controller 220.
  • [0034] Antenna 210 receives signals from transmitters 130 and couples the received signals to terrestrial and GPS receiver units 212 a and 212 b. Terrestrial receiver unit 212 a includes front-end circuitry (e.g., radio frequency (RF) processing circuitry and/or other receiving circuitry) that processes signals transmitted from base stations to derive information used for position determination. Similarly, GPS receiver unit 212 b includes front-end circuitry that processes signals transmitted from GPS satellites to derive information used for position determination. Such information from receiver units 212 a and 212 b may include timing information, the identities and locations of the transmitters whose signals are received, and possibly other information.
  • The information from receiver units [0035] 212 is provided to processing unit 216. Processing unit 216 may initially determine an initial Device Position Estimate and receive from the GPS receiver unit 212 b an Actual Measurement Vector. Processing unit 216 then executes the LMS algorithm, which operates on the computed results to arrive at a final Device Position Estimate.
  • [0036] Memory unit 218 stores various data required for determining position. For example, memory unit 218 may store Almanac and/or Ephemeris information that indicates where each of a number of GPS satellites is located in the sky at any given time. Both the Almanac and Ephemeris are valid for a limited amount of time. The Almanac information is considered to be accurate to approximately 3 kilometers for approximately one week from the time the Almanac is transmitted. The Ephemeris provides information regarding the satellite orbit with an accuracy of approximately 1 meter for approximately 2 hours. The error in both the Almanac and Ephemeris grows as the information ages. The satellite locations may be derived from the Almanac and Ephemeris or from other information obtained by GPS receiver unit 212 b from signals transmitted from the satellites or by terrestrial receiver unit 212 a from signals transmitted via the wireless link.
  • [0037] Controller 220 may receive signaling and/or intermediate results from processing unit 216 and may direct the operation of the processing unit. For example, controller 220 may select the particular types of signals to be processed (e.g., satellite, terrestrial, or both), the particular algorithm to be used (if more than one is available), the parameter values to be used for the selected algorithm, and so on, as described in further detail below.
  • Although not shown in FIG. 2, [0038] receiver device 110 may communicate with a PDE, which may assist in determining the position estimate for the device. The PDE may perform the computations to derive the receiver device position estimate, or may provide certain information used to determine the position estimate (e.g., the location of the GPS satellites and/or base stations).
  • An iterative algorithm may be used to determine the Device Position Estimate based on (1) Transmitter Position Estimates (e.g., information regarding the locations of GPS satellites and/or base stations) and (2) an Actual Measurement Vector (e.g., a set of pseudo-range measurements to the transmitters). For GPS satellites, the transmitter locations may be ascertained based on current Almanac and/or Ephemeris information. [0039]
  • For clarity, the LMS iterative algorithm is specifically described below. However, other iterative algorithms may also be used to determine the Device Position Estimate in accordance with the disclosed method and apparatus. [0040]
  • Each Transmitter Position Estimate is an estimate of the location of one transmitter. The Transmitter Position Estimate, {overscore (s)}[0041] i, for the i-th transmitter is expressed as:
  • {overscore (s)} i =[x i y i z i] and b i
  • where the b[0042] i is the clock bias expressed in the same dimension as x, y, and z and further includes the effects of the reference time. The Transmitter Position Estimates, {overscore (s)}i for i=1, 2, . . . , n, are ascertained once prior to executing the LMS algorithm and are thereafter used by the algorithm to calculate the Device Position Estimate.
  • The LMS algorithm may be performed to determine a set of Device Position Estimates resulting in a final Device Position Estimate as follows. For the k-th iteration (where k=1 for the first iteration of the algorithm), the Device Position Estimate, {overscore (u)}[0043] k, is expressed as:
  • {overscore (u)} k =[x u,k y u,k z u,k b u,k]
  • Prior to the first iteration, an initial Device Position Estimate, {overscore (u)}[0044] 0, is assumed for the receiver device. The value of the initial Device Position Estimate may be selected arbitrarily or based upon current or past knowledge about the position of the receiver device. In one embodiment, the initial Device Position Estimate, {overscore (u)}0, may be: (1) the last Device Position Estimate (if one is available), (2) the center of the earth, (3) the center of the cell of which the receiver device is receiving the signal, or (4) any other location that provides an initial value to the LMS algorithm. The more accurate the initial Device Position Estimate, the faster the algorithm will convergence toward the target solution. A set of values make up a column vector, {overscore (ρ)}k, (hereafter referred to as the “Calculated Measurement Vector”). Each value of the Calculated Measurement Vector is a computed measurement associated with one of the transmitters used for position determination. The values of the Calculated Measurement Vector are calculated from: (1) the Transmitter Position Estimates, {overscore (s)}i, (2) the Device Position Estimate, {overscore (u)}k, resulting from the k-th iteration of the iterative algorithm. Accordingly, it will be understood that there will be one Calculated Measurement Vector, {overscore (ρ)}k, associated with each iteration of the iterative algorithm. It should be further understood that these “measurements” are not made, but are calculated for each iteration based on the values of {overscore (s)}i and {overscore (u)}k. In contrast, as will be more fully disclosed below, an Actual Measurement Vector, {overscore (ρ)}a, is derived from actually measuring the received signals. The elements of the Calculated Measurement Vector, {overscore (ρ)}k, may be any type of measurement that is a function of the location of the receiver device and the clock bias associated with the transmitters. For example, the elements of the Calculated Measurement Vector may relate to ranges, pseudo-ranges as determined at the receiver device (i.e., r−b, where r is the range and b is the clock bias), reverse pseudo-ranges as determined at the base stations based on signals transmitted from the receiver device (i.e., r+b), range differences as determined based on “time difference of arrival” (TDOA), or a combination thereof. Furthermore, the elements of the Calculated Measurement Vector may be computed by any one of a number of known algorithms. In one embodiment, the elements of the Calculated Measurement Vector, {overscore (ρ)}k, are pseudo-ranges. One way to compute the elements of the Calculated Measurement Vector is:
  • ρk,i={square root}{square root over ((x u,k −x i)2+(y u,k −y i)2+(z u,k −z i)2)}−(b u,k −b i)  Eq (1)
  • where each element, ρ[0045] k,i, of the Calculated Measurement Vector, {overscore (ρ)}k, is the pseudo-range between the Device Position Estimate, {overscore (u)}k, for the k-th iteration and the Transmitter Position Estimate, {overscore (s)}i, for the i-th transmitter.
  • A geometry matrix H for the Transmitter Position Estimates, {overscore (s)}[0046] i, and the Device Position Estimate, {overscore (u)}k, is defined as: H k = [ x u , k - x 1 u _ k - s _ 1 y u , k - y 1 u _ k - s _ 1 z u , k - z 1 u _ k - s _ 1 - 1 x u , k - x 2 u _ k - s _ 2 y u , k - y 2 u _ k - s _ 2 z u , k - z 2 u _ k - s _ 2 - 1 x u , k - x n u _ k - s _ n y u , k - y n u _ k - s _ n z u , k - z n u _ k - s _ n - 1 ] . Eq ( 2 )
    Figure US20030004640A1-20030102-M00001
  • A covariance matrix R of the measurements is defined as: [0047]
  • R=E{{overscore (ρ)} a·{overscore (ρ)}a T }−E{{overscore (ρ)} a }·E{{overscore (ρ)} a T}  Eq (3)
  • where: [0048]
  • {overscore (ρ)}a=[ρa,1ρa,2 . . . ρa,n]T
  • where {overscore (ρ)}[0049] a is the Actual Measurement Vector and ρa,n is the measurement element (e.g., pseudo-range measurement) attained from the signal received from the n-th transmitter. E is the statistical expectation. The superscript T denotes the complex conjugate transpose. The covariance matrix R is typically estimated based on the particular design of the receiver device and the system (e.g., the signal bandwidth or the received signal power).
  • “Residual Measurement Error Vector”, {overscore (e)}[0050] k, is then defined as:
  • {overscore (e)} k={overscore (ρ)}a{overscore (ρ)}k  Eq (4)
  • As shown in equation (4), the Residual Measurement Error Vector, {overscore (e)}[0051] k, is equal to the difference between the Actual Measurement Vector, {overscore (ρ)}a, and the Calculated Measurement Vector, {overscore (ρ)}k. It should be noted that one Residual Measurement Error Vector is calculated for each iteration of the algorithm. The Residual Measurement Error Vector, {overscore (e)}k, is calculated from the Device Position Estimate, {overscore (u)}k, determined on the k-th iteration. The next Residual Measurement Error Vector, {overscore (e)}k+1, is then newly calculated for the next iteration from the Device Position Estimate, {overscore (u)}k+1.
  • An Update Vector, Δ{overscore (u)}[0052] k, is computed for each iteration based on the geometry matrix H, the covariance matrix R, and the Residual Measurement Error Vector, {overscore (e)}k, for that iteration. The Update Vector is defined as:
  • Δ{overscore (u)} k=(H T ·R −1 ·H)−1 ·H T ·R −1 ·{overscore (e)} k  Eq (5)
  • This Update Vector represents an estimated error between the current Device Position Estimate, {overscore (u)}[0053] k, and an optimal position estimate. Thus, an updated Device Position Estimate is expressed as:
  • {overscore (u)} k+1 ={overscore (u)} k +Δ{overscore (u)} k  Eq (6)
  • Equations (1) through (6) comprise the computations for one iteration of the LMS algorithm. A number of iterations may be performed to derive a more accurate Device Position Estimate. [0054]
  • In cases of small DOP, if the Transmitter Position Estimates, {overscore (s)}[0055] i, and the Calculated Measurement Vector, {overscore (ρ)}k, are valid and accurate, then the LMS algorithm will most likely converge after a relatively small number of iterations. Accordingly, the Device Position Estimates provided by the LMS algorithm, {overscore (u)}k, asymptotically approach the target solution. The target solution is defined as that solution that results in the minimum value for a metric calculated for the target solution. The metric is a function of the errors that are included in the Residual Measurement Error Vector.
  • However, in cases of large DOP around the target solution, the Update Vector, Δ{overscore (u)}[0056] k, tends to become disproportionately large compared to the Residual Measurement Error Vector. This then leads to repeated overshoots and oscillation around the target solution. Therefore, the algorithm may not converge on the target solution.
  • The disclosed method and apparatus (1) determines when the algorithm is unlikely to convergence (e.g., due to oscillation of the LMS solutions around the target solution) and (2) adjusts the algorithm to increase the likelihood of convergence. It should be noted that in normal circumstances in a position determination system, the algorithm typically will convergence to the target solution after a predetermined number of iterations. This number is normally small (e.g., three). However, the typical number of iterations may be greater if the accuracy of the initial Device Position Estimate is poor. If the algorithm requires more iterations than what is typical, then there is less likelihood of convergence and greater likelihood of oscillation around the target solution. Such oscillation can be identified by an inversion in the direction of two consecutive Update Vectors (i.e., 180° change in phase of the Update Vector). Whether the LMS algorithm is likely to converge can be determined by (1) the number of iterations performed without reaching the target solution and (2) changes in the phase of the Update Vector. [0057]
  • To increase the likelihood of convergence, certain parameters of the LMS algorithm can be adjusted. In accordance with one embodiment of the disclosed method and apparatus, the magnitude of the Update Vector is reduced to (1) avoid significant overshoots and (2) adjust the rate at which the position estimate is updated, both of which help to ensure convergence. [0058]
  • FIG. 3 is a flow diagram illustrating a [0059] process 300 for determining the Device Position Estimate in accordance with one embodiment of the disclosed method and apparatus. The information required to determine the Device Position Estimate is initially received, at step 312. This information typically includes the Transmitter Position Estimates and information used to determine the calculated Measurement Vector (e.g., pseudo-range measurements). An iterative algorithm (e.g., an LMS algorithm) is then executed to determine the Device Position Estimate based on the received information, at step 314.
  • At [0060] step 318, a determination is made as to whether convergence is likely. In accordance with one embodiment, a determination that convergence is likely is based on one or more of the following tests. If convergence is unlikely, then one or more parameters of the algorithm may be adjusted to improve the likelihood of convergence, at step 320, as described below. The process then proceeds to step 322. Otherwise, if convergence is likely, as determined in step 318, the process simply proceeds to step 322.
  • At [0061] step 322, a determination is made as to whether or not the process should be terminated. In accordance with one embodiment of the presently disclosed method and apparatus, once the algorithm is terminated, a final Device Position Estimate is determined. The algorithm may be terminated for any one of a number of reasons, such as the maximum number of iterations have been performed, the algorithm has sufficiently converged, etc. In one such embodiment, a final Device Position Estimate is determined if the magnitude of the Update Vector is smaller than a predetermined threshold, indicating that the algorithm has sufficiently converged. If the algorithm has not terminated, then the process of FIG. 3 continues at step 314.
  • The process shown in FIG. 3 is applicable for any iterative algorithm that may be used to determine position. [0062] Steps 314 through 322 may be viewed as comprising one iteration of the iterative algorithm, and any number of iterations may be performed until sufficiently accurate Device Position Estimate is reached. A further discussion of how the likelihood of convergence is determined and the parameters are adjusted is provided below. For the sake of clarity, the disclosed method and apparatus is described below using an LMS algorithm as one example. However, it will be understood by those skilled in the art that other iterative algorithms could be used to calculate an estimate of the location of a receiver device.
  • FIG. 4A is a flow diagram illustrating a [0063] process 400 performed to determine the Device Position Estimate using the LMS algorithm in accordance with one embodiment of the disclosed method and apparatus. The information required to determine the Device Position Estimate is first received at step 412. This information typically includes: (1) the Transmitter Position Estimates, {overscore (s)}t, (2) the initial Device Position Estimate, {overscore (u)}0; and (3) the Actual Measurement Vector, {overscore (92 )}a. Variables used for the LMS algorithm are then initialized, at step 414. This may entail setting certain variables to their initial values as follows : k=1, g=1, and UV_modify=False, where k is the iteration index, g is a scaling factor used to adjust the magnitude of the Update Vector, Δ{overscore (u)}k, and UV_modify indicates whether or not to adjust the Update Vector.
  • Execution of the LMS algorithm then commences. In the embodiment shown in FIG. 4A, the elements of the Actual Measurement Vector represent pseudo-range measurements. An initial Device Position Estimate is made. A Calculated Measurement Vector is then determined from the initial Device Position Estimate. The Calculated Measurement Vector, {overscore (ρ)}[0064] k, is computed as shown in equation (1), at step 416. Next, the geometry matrix H is updated based on the current Device Position Estimate, {overscore (u)}k, and the Transmitter Position Estimate, {overscore (s)}i, as shown in equation (2), at step 418. The measurement covariance matrix R is also derived, at step 418. The Residual Measurement Error Vector, {overscore (e)}k, is then computed based on the Actual Measurement Vector and the Calculated Measurement Vector as shown in equation (4), at step 420. The Update Vector, Δ{overscore (u)}k, is next computed based on the matrices H and R and the Residual Measurement Error Vector, {overscore (e)}k, as shown in equation (5), at step 422.
  • At [0065] step 426, a determination is made as to whether the LMS algorithm is likely to converge. If not, then one or more parameters of the LMS algorithm and/or the Update Vector, Δ{overscore (u)}k, may be adjusted (e.g., reduced in magnitude) to increase the likelihood of convergence by the algorithm. Parameters used to adjust the Update Vector may also be updated in step 426. Methods for determining whether the algorithm is likely to convergence and adjusting the parameters of the LMS algorithm, including the Update Vector, are described below in FIG. 4B.
  • At step [0066] 428, the new Device Position Estimate, {overscore (u)}k+1, for the receiver device is computed based on the current Device Position Estimate, {overscore (u)}k, and the Update Vector, Δ{overscore (u)}k, as shown in equation (6).
  • A determination is then made as to whether or not the LMS algorithm should be terminated, at [0067] step 430. Termination of the algorithm may be triggered based on one or more criteria, such as: (1) the magnitude of the Update Vector, (2) the number of iterations already performed, or a combination thereof. In one embodiment, the algorithm is terminated if the magnitude of the Update Vector is smaller than a particular update threshold (i.e., |Δ{overscore (u)}k|<Δ{overscore (u)}k—threshold) or if the number of iterations already performed exceeds a particular maximum number of iterations (e.g., k>50). If either of these criteria are satisfied, the current Device Position Estimate is considered to be sufficiently accurate, at step 434, and the process terminates. The Device Position Estimate will in fact be accurate if the algorithm has converged. However the Device Position Estimate will not be as accurate as desired if the algorithm has not converged, but rather was terminated due to some other criterion (e.g., the number of iterations exceed the maximum number of iterations). If a termination criterion has not been satisfied, the iteration index k is incremented by one, at step 432, and the process returns to step 416 to perform another iteration.
  • [0068] Steps 416 through 428 comprise one iteration of the LMS algorithm. Any number of iterations may be executed until the termination criterion or criteria are met.
  • FIG. 4B is a flow diagram of a process to determine whether convergence is likely. Also shown is a process for adjusting parameters of the LMS algorithm in accordance with one embodiment of the disclosed method and apparatus. This process may be used for [0069] step 426 in FIG. 4A.
  • In one embodiment, a determination that convergence of the LMS algorithm is unlikely is in part based on the number of iterations performed without achieving convergence. Convergence is determined by attaining a sufficiently small norm of the Update Vector. As noted above, in normal circumstances in a wireless-assisted GPS position determination system, convergence typically occurs within a small number of iterations. Thus, a determination is made whether the number of iterations performed by the LMS algorithm for this collection of input data exceeds a particular threshold “k_threshold” (i.e., whether k>k_threshold), at [0070] step 448. The value of k_threshold is typically selected based on the number of iterations expected to be performed to converge for a typical application (which may be dependent on the accuracy of the initial device position estimate). Alternatively, any value may also be selected that can provide good performance. As a specific example, k threshold may be set to three. If the number of iterations already performed is greater than k_threshold, the likelihood that convergence by the algorithm will occur is assumed to be low and the process proceeds to step 450 where another test for convergence is performed. Otherwise, if the number of iterations k is less than or equal to k_threshold, a determination is made that the algorithm may yet converge.
  • In another embodiment, a determination that the LMS algorithm is unlikely to convergence is made by examining the Update Vectors. Convergence is determined to be unlikely when the Device Position Estimates repeatedly overshoot the desired Device Position Estimate. Overshooting by the LMS algorithm causes the set of Device Position Estimates resulting from each iteration to “oscillate” around the desired Device Position Estimate. This oscillation can typically be seen by changes in the direction of the Update Vector essentially along the line that is coincident with the direction of the largest DOP. When such phase reversals occur, the phase difference α between two consecutive Update Vectors is close to 180°. This results in the cosine of the phase difference α being close to minus one. In one embodiment, reversal in the direction of the Update Vector is detected if the cosine of α is less than a particular negative cos_threshold. As a specific example, cos_threshold may be selected as −0.95. [0071]
  • As part of step [0072] 450, the cosine of the phase difference α of two consecutive Update Vectors is computed. This cosine is computed by performing a dot product of the two most recent Update Vectors, which can be expressed as: cos ( α k ) = < Δ u ~ k , Δ u ~ k - 1 > Δ u ~ k · Δ u ~ k - 1 , Eq ( 7 )
    Figure US20030004640A1-20030102-M00002
  • where Δũ[0073] k represents the first three components of the Update Vector Δ{overscore (u)}k (i.e., Δũk=[Δxu,kΔyu,kΔzu,k]) and <Δũk, Δũk−1> denotes a dot product of the first three components of the two Update Vectors, which can be expressed as:
  • ũ k ,Δũ k−1>=(Δx i ·Δx i−1)+(Δy i ·Δy i−1)+(Δz i ·Δz i−1)
  • A determination is then made as to whether the cosine of the phase difference is less than the cos_threshold, (i.e., cos(α[0074] k)<cos_threshold), at step 452. If the answer is yes, then a phase reversal has been detected.
  • In one embodiment, the first time a phase reversal in the Update Vector is detected after the number of iterations exceeds k_threshold, a variable UV_modify is set to True. Each subsequent time the phase reversal is detected, the factor g is multiplied by a positive factor greater than unity (e.g., two). This larger value for the factor g is then used to effectively divide the magnitude of the Update Vector by a corresponding amount. The adjusted Update Vector may be smaller than optimal, and under normal operating conditions this smaller Update Vector would lead to a slower convergence rate. However, for unstable scenarios, the smaller Update Vector ensures that, given a sufficient number of iterations, the LMS algorithm has a greater likelihood of convergence. It should be noted that many other methods could be implemented for reducing the size of the Update Vector in response to such phase reversals. In one example of such an alternative, the number of phase reversals does not correspond to the amount the Update Vector is reduced. Rather the Update Vector is reduced by a predetermined amount each time a reduction is made. [0075]
  • In one embodiment, the norm of the Update Vector is set to the norm of the Residual Measurement Error Vector at the time a phase reversal is detected. This is achieved by setting a variable d to the magnitude of the Residual Measurement Error Vector (i.e., d=|{overscore (e)}[0076] k|) each time the phase reversal is detected. This variable d is thereafter used to determine the magnitude of the Update Vector, as described below.
  • For the embodiments described above, a parameter of the LMS algorithm (e.g., the Update Vector) may be continually (i.e., “adaptively”) adjusted based on the monitored performance of the algorithm. In one such embodiment, such an “adaptive” adjustment is achieved based on the values of the variable d and the factor g. Thus, in such an embodiment, once the variable UV_modify has been set to True and regardless of whether or not a phase reversal has been detected, the Update Vector is adjusted by the variable d and factor g, as described below. These actions are performed by the following steps. [0077]
  • If a phase reversal of the Update Vector was detected at [0078] step 452, a determination is made whether the variable UV_modify has been set to True, at step 454. If the answer is no, indicating that this is the first phase reversal detected, then the variable UV_modify is set to True, at step 456. Otherwise, if UV modify has already been set to True by a previously detected phase reversal, then the factor g is increased (e.g., doubled), at step 458. In either case, the variable d is set to the magnitude of the Residual Measurement Error Vector, at step 460. The Update Vector is then adjusted based on the variable d and factor g, at step 462.
  • In one embodiment, the Update Vector is adjusted as follows: [0079] Δ u _ k = Δ u _ k · d Δ u _ k · g . Eq ( 8 )
    Figure US20030004640A1-20030102-M00003
  • As shown in equation (8), the adjusted Update Vector retains the phase of the original Update Vector, but its norm is set to the magnitude of the Residual Measurement Error Vector at the time the last phase reversal was detected (as captured by the variable d) and further reduced by the factor g. The adjusting ensures that the magnitude of the Update Vector is smaller each time the phase reversal is detected, since the amount of overshoot is less each time. In this manner, the step size can be made dependent on the monitored performance of the LMS algorithm. Moreover, by increasing (e.g., doubling) the factor g each time a phase reversal is detected, the size of the Update Vector is correspondingly reduced. After adjusting the Update Vector at [0080] step 462, the process terminates.
  • If a phase reversal of the Update Vector was not detected at [0081] step 452, a determination is made at to whether the variable UV_modify has been set to True, at step 464. If the answer is no, indicating that the LMS may be taking longer to converge, but no oscillation in the Update Vector has been detected, then the process terminates. Otherwise, if UV-modify has already been set to True by a previous detected phase reversal, then the Update Vector is adjusted based on the current values for the variable d and factor g, at step 462. In the embodiment shown in FIG. 4B, the Update Vector is adjusted by the variable d and factor g for each iteration, once UV-modify has been set to True, regardless of whether or not phase reversal of the Update Vector has been detected.
  • For clarity, FIGS. 4A and 4B illustrate examples of the disclosed method for performing the LMS algorithm. However, these examples may also be applied to other algorithms used to determine the Device Position Estimate. For clarity, FIGS. 4A and 4B also describe a specific implementation for determining that convergence is unlikely based on the number of iterations and the Update Vectors. Also described is an example of one method by which the LMS algorithm can be adjusted to improve the likelihood of convergence. It should be noted that these methods are merely presented as examples of ways to determine whether convergence is likely and to adjust the parameters of the LMS algorithm. Other methods are clearly possible and would be understood by those skilled in the art from the examples provided herein. [0082]
  • In FIG. 2, the elements of receiver device [0083] 110 (e.g., processing unit 216 and controller 220) may be implemented as a digital signal processor (DSP), an application specific integrated circuit (ASIC), a processor, a microprocessor, a controller, a microcontroller, a field programmable gate array (FPGA), a programmable logic device, other electronic unit, or any combination thereof designed to perform the functions described herein. Certain aspects of the disclosed method and apparatus may be implemented in hardware, software, or a combination of both. For example, the processing to determine the Device Position Estimate (e.g., the LMS algorithm) may be performed based on program codes stored within memory unit 218 and executed by a processor (processing unit 216 or controller 220 in FIG. 2). Alternatively, such processing may be performed by dedicated hardware or some combination of hardware and software.
  • The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the presently disclosed method and apparatus. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.[0084]

Claims (33)

What is claimed is:
1. A method for determining an accurate device position estimate, comprising:
executing an iterative algorithm to provide a device position estimate;
determining whether convergence of the iterative algorithm is likely; and
if convergence of the iterative algorithm is unlikely, adjusting at least one parameter of the iterative algorithm to increase the likelihood of convergence.
2. The method of claim 1, wherein the iterative algorithm is based on an initial device position estimate, transmitter position estimates for a plurality of transmitters, and an actual measurement vector.
3. The method of claim 1, wherein convergence is determined to be unlikely if a particular number of iterations is performed without achieving convergence.
4. The method of claim 2, further comprising:
calculating an update vector for a current device position estimate based on a Residual Measurement Error Vector, the Transmitter Position Estimates, and the most recent Device Position Estimate calculated by the iterative algorithm;
calculating the norm of the update vector; and
wherein convergence is determined to be unlikely based on the value of the norm of the update vector.
5. The method of claim 4, wherein the metric is a function of the errors that are included in the Residual Measurement Error Vector.
6. The method of claim 5, wherein the function is linear.
7. The method of claim 5, wherein the function is quadratic.
8. The method of claim 1, wherein the iterative algorithm is a least mean square (LMS) algorithm.
9. A method for determining a device position estimate, comprising:
determining transmitter position estimates for a plurality of transmitters and measuring an actual measurement vector for the transmitters;
determining an initial device position estimate;
determining a calculated measurement vector based on the transmitter position estimates and the initial device position estimate;
revising the device position estimate based on the difference between the calculated measurement vector, and the actual measurement vector;
determining that convergence of the current device position estimate toward a final device position estimate is unlikely; and
if convergence is determined to be unlikely, adjusting the manner in which the revised device position estimate is determined to increase the likelihood of convergence.
10. The method of claim 9, wherein the revised device position estimate is calculated using a least mean square (LMS) algorithm.
11. The method of claim 9, wherein the difference between the calculated measurement vector, and the actual measurement vector is a residual measurement vector, and the manner for determining the revised device position estimate is adjusted by scaling a value derived from the residual measurement error vector.
12. The method of claim 11, wherein the value derived from the residual measurement error vector is an update vector.
13. The method of claim 9, further comprising:
determining a geometry matrix based on the current device position estimate and the transmitter position estimate; and
determining a measurement covariance matrix, and
wherein the adjustment is further based on the geometry matrix and the measurement covariance matrix.
14. The method of claim 9, further comprising:
repeating the steps of claim 9 iteratively until either the iterative algorithm converges to a final device position estimate, or a determination is made that convergence is unlikely.
15. The method of claim 9, wherein the step of adjusting the device position estimate includes:
determining a scaling factor which effects the amount of the adjustment to the device position estimate.
16. The method of claim 15, wherein the scaling factor is increased with each iteration during which convergence is determined to be unlikely.
17. The method of claim 16, wherein the convergence is determined to be unlikely based on phase reversals in the update vector.
18. The method of claim 16, wherein the scaling factor is not increased in response to the first predetermined number of iterations resulting in a determination that convergence is unlikely.
19. The method of claim 12, further including:
limiting the magnitude of the update vector based on the magnitude of a residual measurement error vector as computed when a determination has been made that convergence is unlikely.
20. The method of claim 14, wherein non-convergence is detected based in part on the number of iterations performed without reaching the target position estimate.
21. The method of claim 20, wherein non-convergence is further detected based on the update vector.
22. The method of claim 21, wherein non-convergence is detected if the number of iterations exceeds a first threshold and a phase difference between two consecutive updates for two iterations is less than a second threshold.
23. The method of claim 22, wherein the phase difference between two consecutive updates is derived by performing a dot product on the two consecutive updates.
24. The method of claim 14, further comprising:
providing the current device position estimate as a final position estimate for the device if one or more criteria for termination are met.
25. The method of claim 24, wherein the one or more criteria for termination are met if the number of iterations exceeds a first threshold.
26. The method of claim 24, wherein the one or more criteria for termination are met if magnitude of the update is within a second threshold.
27. The method of claim 9, wherein the actual measurement vector corresponds to pseudo-ranges indicative of distances between the device and the transmitters.
28. The method of claim 9, wherein the transmitters are satellites in a Global Positioning Satellite (GPS) system.
29. The method of claim 9, wherein the transmitters include base stations in a wireless communication system.
30. The method of claim 29, wherein the wireless communication system is a CDMA system.
31. A computer program product for determining the location of a device, comprising:
code for receiving position estimates for a plurality of transmitters and initial measurements for the transmitters;
code for computing an update for a current device position estimate based on the Transmitter Position Estimates, the current device position estimate, and the initial measurements, wherein the current device position estimate is initially set to an initial estimate for the location of the device;
code for detecting non-convergence of the current device position estimate toward a target position estimate;
code for, if non-convergence is detected, adjusting the update to increase likelihood of convergence to the target position estimate;
code for updating the current device position estimate based on the update; and
a data storage medium configured to store the codes.
32. A receiver device comprising:
at least one receiver unit configured to receive and process a plurality of signals from a plurality of transmitters to derive initial measurements for the transmitters; and
a processor coupled to the at least one receiver unit and configured to
receive the initial measurements and position estimates for the transmitters,
compute an update for a current device position estimate based on the transmitters position estimates, the current device position estimate, and the initial measurements, wherein the current device position estimate is initially set to an initial estimate for the location of the device,
detect non-convergence of the current device position estimate toward a target position estimate,
if non-convergence is detected, adjust the update to increase likelihood of convergence to the target position estimate,
update the current device position estimate based on the update, and
repeat the compute, detect, adjust, and update for a plurality of iterations
33. The receiver device of claim 32, wherein the processor is further configured to
scale the update by a scaling factor to reduce magnitude of the update,
wherein the scaling factor is increased based on detected events indicative of non-convergence.
US09/922,472 2001-06-25 2001-08-03 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision Expired - Lifetime US6505122B1 (en)

Priority Applications (12)

Application Number Priority Date Filing Date Title
US09/922,472 US6505122B1 (en) 2001-06-25 2001-08-03 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
CA2451290A CA2451290C (en) 2001-06-25 2002-06-17 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
EP02737535A EP1399756B1 (en) 2001-06-25 2002-06-17 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
JP2003507578A JP4180505B2 (en) 2001-06-25 2002-06-17 Method and apparatus for providing an accurate position estimate in a state where accuracy is significantly diluted
IL15927802A IL159278A0 (en) 2001-06-25 2002-06-17 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
CNB028165624A CN100420958C (en) 2001-06-25 2002-06-17 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
PCT/US2002/019406 WO2003001236A2 (en) 2001-06-25 2002-06-17 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
MYPI20022343A MY130562A (en) 2001-06-25 2002-06-21 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
ARP020102381A AR034627A1 (en) 2001-06-25 2002-06-24 METHOD AND APPLIANCE TO CONTRIBUTE PRECISE POSITION ESTIMATES IN CASES OF SEVERE DILUTION OF PRECISION
TW091113877A TWI229743B (en) 2001-06-25 2002-06-25 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
IL159278A IL159278A (en) 2001-06-25 2003-12-09 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
HK05108326.3A HK1076512A1 (en) 2001-06-25 2005-09-22 Method and apparatus for providing accurate position estimates in instances ; of severe dilution of precision

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US30096001P 2001-06-25 2001-06-25
US09/922,472 US6505122B1 (en) 2001-06-25 2001-08-03 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision

Publications (2)

Publication Number Publication Date
US20030004640A1 true US20030004640A1 (en) 2003-01-02
US6505122B1 US6505122B1 (en) 2003-01-07

Family

ID=26972069

Family Applications (1)

Application Number Title Priority Date Filing Date
US09/922,472 Expired - Lifetime US6505122B1 (en) 2001-06-25 2001-08-03 Method and apparatus for providing accurate position estimates in instances of severe dilution of precision

Country Status (11)

Country Link
US (1) US6505122B1 (en)
EP (1) EP1399756B1 (en)
JP (1) JP4180505B2 (en)
CN (1) CN100420958C (en)
AR (1) AR034627A1 (en)
CA (1) CA2451290C (en)
HK (1) HK1076512A1 (en)
IL (2) IL159278A0 (en)
MY (1) MY130562A (en)
TW (1) TWI229743B (en)
WO (1) WO2003001236A2 (en)

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007536510A (en) * 2004-05-07 2007-12-13 ナヴコム テクノロジー インコーポレイテッド GPS navigation using continuous difference of carrier-phase measurements
WO2008121178A1 (en) 2007-03-30 2008-10-09 Sony Corporation System and method for effectively performing enhanced mobile-device location procedures
US20090085802A1 (en) * 2006-04-18 2009-04-02 Andrew Corporation Method and apparatus for geolocation determination
CN101611328A (en) * 2007-01-22 2009-12-23 国立大学法人大阪教育大学 Position detecting device, method for detecting position, data judging device, data judging method, computer program and storage medium
US20100240389A1 (en) * 2009-03-18 2010-09-23 Andrew Llc System and method for locating mobile device in wireless communication network
US20100238070A1 (en) * 2009-03-18 2010-09-23 Andrew Llc System and method for concurrently determining locations of mobile device in wireless communication network
CN102175991A (en) * 2011-01-16 2011-09-07 哈尔滨工程大学 Target positioning method based on maximum positioning likelihood sensor configuration
US8538682B1 (en) * 2009-04-24 2013-09-17 Qualcomm Incorporated Systems and methods for satellite navigation using locally generated ephemeris data
US20140266911A1 (en) * 2013-03-15 2014-09-18 Nextnav, Llc Directional pruning of transmitters to improve position determination
US20140274125A1 (en) * 2013-03-15 2014-09-18 Nextnav, Llc Performance enhancements for local network of beacons
US20140365488A1 (en) * 2012-01-18 2014-12-11 Sensewhere Limited Positioning system
US8963764B1 (en) 2011-01-14 2015-02-24 Lockheed Martin Corporation Ship heading and pitch using satellite ephemerides and radar range measurement of satellite
US9459344B1 (en) * 2011-01-14 2016-10-04 Lockheed Martin Corporation Ship position and velocity using satellite ephemerides and radar range measurement of satellite
USRE46237E1 (en) * 2004-09-02 2016-12-13 Nariste Networks Pty Ltd Self monitoring GPS on mobile device
US9578459B2 (en) * 2015-06-05 2017-02-21 Qualcomm Incorporated Hybrid positioning techniques based on RTT and TOA/TDOA
US9885772B1 (en) 2014-08-26 2018-02-06 Vencore Labs, Inc. Geolocating wireless emitters
US10234566B2 (en) * 2015-02-11 2019-03-19 Trimble Inc. Fast GNSS receiver convergence by relative positioning
US10234565B2 (en) * 2015-02-11 2019-03-19 Trimble Inc. Global navigation satellite system receiver convergence selection
US10241213B2 (en) * 2015-02-11 2019-03-26 Trimble Inc. Timer initiated convergence of a GNSS receiver
CN110730502A (en) * 2019-10-23 2020-01-24 珠海优特电力科技股份有限公司 Positioning method and device
CN113219504A (en) * 2017-02-28 2021-08-06 荣耀终端有限公司 Positioning information determining method and device
WO2022077938A1 (en) * 2020-10-15 2022-04-21 华为技术有限公司 Ephemeris data transmission method and device, storage medium, chip, and system
US20220357464A1 (en) * 2021-05-10 2022-11-10 Qualcomm Incorporated Determining position information of mobile devices
CN117092588A (en) * 2023-10-20 2023-11-21 中国科学院深海科学与工程研究所 Method for estimating clock deviation of underwater sound positioning system

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8255149B2 (en) 1999-07-12 2012-08-28 Skybitz, Inc. System and method for dual-mode location determination
US20040143392A1 (en) * 1999-07-12 2004-07-22 Skybitz, Inc. System and method for fast acquisition reporting using communication satellite range measurement
US7660588B2 (en) * 2002-10-17 2010-02-09 Qualcomm Incorporated Method and apparatus for improving radio location accuracy with measurements
US7705775B2 (en) * 2005-12-30 2010-04-27 General Motors Llc Method of improving a vehicle emergency call network
CN100461969C (en) * 2006-09-22 2009-02-11 华为技术有限公司 A method to position the mobile station
KR20090076171A (en) * 2008-01-07 2009-07-13 삼성전자주식회사 The method for estimating position and the apparatus thereof
CN101726746B (en) * 2009-12-14 2012-11-14 北京航空航天大学 Intermediate frequency direct sequence spread spectrum receiver for satellite ranging
RU2490661C1 (en) * 2012-03-06 2013-08-20 Федеральное государственное унитарное предприятие "18 Центральный научно-исследовательский институт" Министерства обороны Российской Федерации Method of determining coordinates of short-wave radio-frequency source
CN102967870A (en) * 2012-11-02 2013-03-13 安凯 Fast iterative method of global position system (GPS) positioning
CN103399334B (en) * 2013-07-30 2015-01-07 中国科学院国家天文台 Method for improving positioning precision of satellite navigation system on basis of ultra-precise code
CN106569222A (en) * 2016-11-01 2017-04-19 湖北航天技术研究院总体设计所 Azimuth measurement method based on distance measurement principle
CN106959456B (en) * 2017-03-27 2019-04-26 中国电建集团西北勘测设计研究院有限公司 A kind of GNSS SURVEYING CONTROL NETWORK Accuracy Estimation
CN109765597A (en) * 2019-03-01 2019-05-17 广州达安临床检验中心有限公司 Medicine cold chain localization method, device, equipment and storage medium

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5648901A (en) * 1990-02-05 1997-07-15 Caterpillar Inc. System and method for generating paths in an autonomous vehicle
US5323322A (en) * 1992-03-05 1994-06-21 Trimble Navigation Limited Networked differential GPS system
DE4332945A1 (en) * 1993-09-28 1995-03-30 Bosch Gmbh Robert Positioning and navigation device with satellite support
US5848373A (en) * 1994-06-24 1998-12-08 Delorme Publishing Company Computer aided map location system
US5802492A (en) * 1994-06-24 1998-09-01 Delorme Publishing Company, Inc. Computer aided routing and positioning system
US5745868A (en) * 1995-12-26 1998-04-28 Motorola, Inc. Method for rapid recovery from dead reckoning system heading loss
US6133874A (en) * 1996-03-08 2000-10-17 Snaptrack, Inc. Method and apparatus for acquiring satellite positioning system signals
WO1998002762A2 (en) * 1996-07-12 1998-01-22 Eagle Eye Technologies, Inc. Method and apparatus for precision geolocation
US6246361B1 (en) * 1999-06-28 2001-06-12 Gary Sutton Method and apparatus for determining a geographical location of a mobile communication unit
US6052647A (en) * 1997-06-20 2000-04-18 Stanford University Method and system for automatic control of vehicles based on carrier phase differential GPS
US6009334A (en) * 1997-11-26 1999-12-28 Telefonaktiebolaget L M Ericsson Method and system for determining position of mobile radio terminals
US6208297B1 (en) * 1998-10-09 2001-03-27 Cell-Loc Inc. Methods and apparatus to position a mobile receiver using downlink signals, part I
US6331835B1 (en) * 1999-02-02 2001-12-18 The Charles Stark Draper Laboratory, Inc. Deeply-integrated adaptive GPS-based navigator with extended-range code tracking
US6351711B1 (en) * 2000-09-14 2002-02-26 Sirf Technology, Inc. GPS navigation system using neural networks

Cited By (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007536510A (en) * 2004-05-07 2007-12-13 ナヴコム テクノロジー インコーポレイテッド GPS navigation using continuous difference of carrier-phase measurements
USRE46237E1 (en) * 2004-09-02 2016-12-13 Nariste Networks Pty Ltd Self monitoring GPS on mobile device
USRE48206E1 (en) 2004-09-02 2020-09-15 Robert Anderson Malaney Reduced power use in mobile GPS-based technologies
US20090085802A1 (en) * 2006-04-18 2009-04-02 Andrew Corporation Method and apparatus for geolocation determination
US7696922B2 (en) * 2006-04-18 2010-04-13 Andrew Llc Method and apparatus for geolocation determination
US8396691B2 (en) * 2007-01-22 2013-03-12 Osaka Kyoiku University Position detector, position detecting method, data determination apparatus, data determination method, computer program, and storage medium
CN101611328A (en) * 2007-01-22 2009-12-23 国立大学法人大阪教育大学 Position detecting device, method for detecting position, data judging device, data judging method, computer program and storage medium
US20100019962A1 (en) * 2007-01-22 2010-01-28 Osaka Kyoiku University Position detector, position detecting method, data determination apparatus, data determination method, computer program, and storage medium
WO2008121178A1 (en) 2007-03-30 2008-10-09 Sony Corporation System and method for effectively performing enhanced mobile-device location procedures
EP2132529A4 (en) * 2007-03-30 2015-06-10 Sony Corp System and method for effectively performing enhanced mobile-device location procedures
WO2010107865A1 (en) * 2009-03-18 2010-09-23 Andrew Llc System and method for locating mobile device in wireless communication network
US8160610B2 (en) 2009-03-18 2012-04-17 Andrew Llc System and method for locating mobile device in wireless communication network
AU2010226679B2 (en) * 2009-03-18 2015-07-23 Maple Acquisition Llc System and method for locating mobile device in wireless communication network
US20100238070A1 (en) * 2009-03-18 2010-09-23 Andrew Llc System and method for concurrently determining locations of mobile device in wireless communication network
US9392521B2 (en) 2009-03-18 2016-07-12 Telecommunication Systems, Inc. System and method for concurrently determining locations of mobile device in wireless communication network
US20100240389A1 (en) * 2009-03-18 2010-09-23 Andrew Llc System and method for locating mobile device in wireless communication network
US8538682B1 (en) * 2009-04-24 2013-09-17 Qualcomm Incorporated Systems and methods for satellite navigation using locally generated ephemeris data
US9459344B1 (en) * 2011-01-14 2016-10-04 Lockheed Martin Corporation Ship position and velocity using satellite ephemerides and radar range measurement of satellite
US8963764B1 (en) 2011-01-14 2015-02-24 Lockheed Martin Corporation Ship heading and pitch using satellite ephemerides and radar range measurement of satellite
CN102175991A (en) * 2011-01-16 2011-09-07 哈尔滨工程大学 Target positioning method based on maximum positioning likelihood sensor configuration
US20140365488A1 (en) * 2012-01-18 2014-12-11 Sensewhere Limited Positioning system
US20170091223A1 (en) * 2012-01-18 2017-03-30 Sensewhere Limited Positioning system
US9639557B2 (en) * 2012-01-18 2017-05-02 Sensewhere Limited Positioning system
US20140266911A1 (en) * 2013-03-15 2014-09-18 Nextnav, Llc Directional pruning of transmitters to improve position determination
WO2014151164A1 (en) * 2013-03-15 2014-09-25 Nextnav, Llc Directional pruning of transmitters to improve position determination
US20140274125A1 (en) * 2013-03-15 2014-09-18 Nextnav, Llc Performance enhancements for local network of beacons
US9322900B2 (en) * 2013-03-15 2016-04-26 Nextnav, Llc Performance enhancements for local network of beacons
US9964647B2 (en) * 2013-03-15 2018-05-08 Nextnav, Llc Directional pruning of transmitters to improve position determination
US9885772B1 (en) 2014-08-26 2018-02-06 Vencore Labs, Inc. Geolocating wireless emitters
US10234566B2 (en) * 2015-02-11 2019-03-19 Trimble Inc. Fast GNSS receiver convergence by relative positioning
US10234565B2 (en) * 2015-02-11 2019-03-19 Trimble Inc. Global navigation satellite system receiver convergence selection
US10241213B2 (en) * 2015-02-11 2019-03-26 Trimble Inc. Timer initiated convergence of a GNSS receiver
US9578459B2 (en) * 2015-06-05 2017-02-21 Qualcomm Incorporated Hybrid positioning techniques based on RTT and TOA/TDOA
CN113219504A (en) * 2017-02-28 2021-08-06 荣耀终端有限公司 Positioning information determining method and device
CN110730502A (en) * 2019-10-23 2020-01-24 珠海优特电力科技股份有限公司 Positioning method and device
WO2022077938A1 (en) * 2020-10-15 2022-04-21 华为技术有限公司 Ephemeris data transmission method and device, storage medium, chip, and system
US20220357464A1 (en) * 2021-05-10 2022-11-10 Qualcomm Incorporated Determining position information of mobile devices
WO2022241346A1 (en) * 2021-05-10 2022-11-17 Qualcomm Incorporated Determining position information of mobile devices
CN117092588A (en) * 2023-10-20 2023-11-21 中国科学院深海科学与工程研究所 Method for estimating clock deviation of underwater sound positioning system

Also Published As

Publication number Publication date
CN1633604A (en) 2005-06-29
JP4180505B2 (en) 2008-11-12
AR034627A1 (en) 2004-03-03
IL159278A0 (en) 2004-06-01
MY130562A (en) 2007-06-29
WO2003001236A2 (en) 2003-01-03
EP1399756B1 (en) 2012-05-30
CN100420958C (en) 2008-09-24
TWI229743B (en) 2005-03-21
CA2451290C (en) 2012-02-07
JP2005521857A (en) 2005-07-21
EP1399756A2 (en) 2004-03-24
CA2451290A1 (en) 2003-01-03
WO2003001236A3 (en) 2003-04-24
IL159278A (en) 2009-08-03
US6505122B1 (en) 2003-01-07
HK1076512A1 (en) 2006-01-20

Similar Documents

Publication Publication Date Title
US6505122B1 (en) Method and apparatus for providing accurate position estimates in instances of severe dilution of precision
US6429809B1 (en) Method and apparatus for determining location using a coarse position estimate
US6832155B2 (en) Methods and apparatus for determining phase ambiguities in ranging and navigation systems
KR100924937B1 (en) Method and apparatus providing improved position estimate based on an initial coarse position estimate
RU2327303C2 (en) Positioning of wireless communication terminal device in mixed positioning system
US8255160B2 (en) Integrated mobile terminal navigation
US20220299652A1 (en) Gnss standard point positioning method based on spherical harmonics
US6587692B1 (en) Location determination using weighted ridge regression
US9291716B2 (en) Method and apparatus for determining physical characteristics of a receiving device in a navigation system
KR101761782B1 (en) System and method for determining the protection level
US6865478B2 (en) Determining the position of a receiver and/or the system time of a positioning system
US8358241B2 (en) Position calculating method and position calculating device
US20060063537A1 (en) Method and apparatus for determining position of mobile communication terminal
US20210325548A1 (en) Matching for gnss signals
CN115616623A (en) Satellite navigation time service method and device based on current model Kalman filtering
US20110181462A1 (en) System and Method for Positioning with GNSS Using Multiple Integer Candidates
US20030054759A1 (en) Integrity monitoring for geo-location systems
US6563461B1 (en) System, method, and software for non-iterative position estimation using range measurements
US20050080557A1 (en) Method and a system in positioning, and a device
JP3544887B2 (en) GPS receiving device and GPS receiving method applied to this device
Lerner Gauss-Newton Filter with Considered Covariance: Application to Geolocation

Legal Events

Date Code Title Description
AS Assignment

Owner name: QUALCOMM INCORPORATED, CALIFORNIA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:BENCHETRIT, URI;REEL/FRAME:012398/0712

Effective date: 20011022

STCF Information on status: patent grant

Free format text: PATENTED CASE

FPAY Fee payment

Year of fee payment: 4

AS Assignment

Owner name: QUALCOMM INCORPORATED, CALIFORNIA

Free format text: CORRECTIVE ASSIGNMENT TO CORRECT THE INVENTOR NAME URI BENCHETRIT INCORRECT, SHOULD BE ALKINOOS H. VAYANOS PREVIOUSLY RECORDED ON REEL 012398 FRAME 0712;ASSIGNOR:VAYANOS, ALKINOOS H.;REEL/FRAME:018645/0985

Effective date: 20011022

FPAY Fee payment

Year of fee payment: 8

FPAY Fee payment

Year of fee payment: 12