US20100190512A1 - Method for ambiguity resolution in location determination - Google Patents

Method for ambiguity resolution in location determination Download PDF

Info

Publication number
US20100190512A1
US20100190512A1 US12/754,360 US75436010A US2010190512A1 US 20100190512 A1 US20100190512 A1 US 20100190512A1 US 75436010 A US75436010 A US 75436010A US 2010190512 A1 US2010190512 A1 US 2010190512A1
Authority
US
United States
Prior art keywords
location determination
received power
ranging signals
power levels
solutions
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.)
Abandoned
Application number
US12/754,360
Inventor
Cristina A. Seibert
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 US12/754,360 priority Critical patent/US20100190512A1/en
Assigned to QUALCOMM INCORPORATED reassignment QUALCOMM INCORPORATED ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: SEIBERT, CRISTINA A
Publication of US20100190512A1 publication Critical patent/US20100190512A1/en
Abandoned 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/0257Hybrid positioning
    • G01S5/0268Hybrid positioning by deriving positions from different combinations of signals or of estimated positions in a single positioning system
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0252Radio frequency fingerprinting
    • G01S5/02521Radio frequency fingerprinting using a radio-map
    • G01S5/02524Creating or updating the radio-map
    • G01S5/02527Detecting or resolving anomalies in the radio frequency fingerprints of the radio-map
    • 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/0257Hybrid positioning
    • G01S5/0258Hybrid positioning by combining or switching between measurements derived from different systems
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01LSEMICONDUCTOR DEVICES NOT COVERED BY CLASS H10
    • H01L2224/00Indexing scheme for arrangements for connecting or disconnecting semiconductor or solid-state bodies and methods related thereto as covered by H01L24/00
    • H01L2224/01Means for bonding being attached to, or being formed on, the surface to be connected, e.g. chip-to-package, die-attach, "first-level" interconnects; Manufacturing methods related thereto
    • H01L2224/42Wire connectors; Manufacturing methods related thereto
    • H01L2224/47Structure, shape, material or disposition of the wire connectors after the connecting process
    • H01L2224/48Structure, shape, material or disposition of the wire connectors after the connecting process of an individual wire connector
    • H01L2224/4805Shape
    • H01L2224/4809Loop shape
    • H01L2224/48091Arched

Definitions

  • the disclosed subject matter relates generally to methods for location determination, and more particularly to methods using various heuristics techniques for resolving ambiguity in location determination in environments with or without noise.
  • range-based location determination systems measurements of ranging signals from a plurality of sources are converted to distance information associated with the source of each ranging signal. Distances to different sources with known locations are combined to solve for the unknown user location via geometric techniques known, for example as trilateration (a.k.a. triangulation). If delay of ranging signals cannot be known reliably (e.g. in asynchronous systems where the user clock is not synchronized to the network), location determination algorithms may treat user clock temporal bias as another unknown, to be solved for by the trilateration process, using an additional measurement.
  • trilateration a.k.a. triangulation
  • a location determination system is ambiguous if more than one location determination solution set of user coordinates and clock temporal bias is consistent with a set of distance measurements.
  • Location determination systems can produce ambiguous location determination solutions in three distinct ways: first, ambiguity can be caused by insufficient measurements; second, ambiguity may be introduced by the properties of the algorithm employed in location determination; and third, ambiguity may be introduced by the presence of noisy measurements.
  • a system has an insufficient number of measurements when the number of unknowns is greater than or equal to the number of independent measurements.
  • the unknowns are the two-dimensional user spatial coordinates and user clock temporal bias.
  • FIG. 1 There are three unknowns, namely the mobile station latitude, longitude and clock temporal bias.
  • Given three independent distance measurements there are two possible location determination solutions, depicted at the intersection of each set of circles.
  • the nature of the algorithm used for locating a user can also be a source of ambiguity.
  • a well known algorithm that is susceptible to ambiguity is described in the U.S. Pat. No. 6,289,280.
  • This algorithm solves for unknowns using a closed form system of equations. Because it solves for the user location algebraically, this algorithm runs efficiently, making it suitable to applications and devices with time or resource constraints.
  • the solution uses linear algebra manipulations to combine the measurements into a system of quadratic equations where the number of equations equals the number of unknowns. Two solutions are produced associated with the two roots of the quadratic equations. The two solutions form an ambiguous set of solutions which needs to be resolved by additional means.
  • the unknowns are the two-dimensional user spatial coordinates and user clock temporal bias.
  • the system can be said to have a sufficient number of measurements to unambiguously solve for the user location.
  • the algebraic method is used, the four measurements are combined into three “average” measurements and two solutions corresponding to these averages are identified, as shown in FIG. 2 .
  • noisy measurements can lead to error in the determination of user location.
  • a method for location determination in noisy environments by assuming the noise to be a discrete variable with known or computable statistical parameters.
  • a set of adjusted measurements and corresponding solution are generated for each assumed noise level.
  • Such location determination system is ambiguous, thus also warranting ambiguity resolution techniques.
  • FIG. 3 consider three noise levels, each 100 meters apart, associated with the measurement from base station BS 2 .
  • a set of circles is plotted as before, with a radius corresponding to the sum of distance measurement (in the case of BS 2 , this measurement is adjusted by the assumed noise level) and the clock temporal bias computed.
  • a method for resolving ambiguity in location determination with N ambiguous location determination solutions using clock temporal bias value includes the following steps: generating a probability distribution function (PDF) model for clock temporal bias, obtaining N clock temporal bias values for the N ambiguous location determination solutions, inserting each of the N clock temporal bias values into the PDF model, evaluating the PDF model to get N PDF values, setting N goodness metrics to the N PDF values, comparing the N goodness metrics, defining a maximum goodness metric as the largest of the N goodness metrics, and selecting a final location determination solution having the maximum goodness metric.
  • PDF probability distribution function
  • a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions using the order of receipt (a.k.a. time of arrival) of a plurality of ranging signals includes the following steps: ranking the plurality of ranging signals based on the order of actual time of arrival (TOA) from earliest to latest, ranking the plurality of ranging signals based on an expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions, and comparing the ranking of the ranging signals based on the order of the actual time of arrival (TOA) and based on the expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions.
  • TOA actual time of arrival
  • TOA actual time of arrival
  • a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions using the order of received power levels of a plurality of ranging signals includes the following steps: ranking the plurality of ranging signals based on the order of received power levels from strongest to weakest, ranking the plurality of ranging signals based on an expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions, and comparing the ranking of the ranging signals based on the order of received power levels and based on the expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions.
  • a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions includes the following steps: deriving a plurality of distances D 1 corresponding to each of the plurality of ambiguous location determination solutions for each of M sources, obtaining a plurality of original distances D 2 , comparing each of the plurality of distances D 1 with each of its corresponding plurality of original distances D 2 and computing a plurality of error measurements, and selecting a final location determination solution as having the lowest value of the plurality of error measurements.
  • a method for resolving ambiguity in location determination in a noisy environment includes the following steps: selecting L noise levels iteratively within a range [A, B] with an increment I for a first of Q ranging signals, repeating the selecting step Q minus 1 times for each of the rest of the Q ranging signals to produce a plurality of L noise levels, creating a plurality of ambiguous location determination solutions for each of the plurality of L noise levels based on a geometric technique, selecting the final location determination solution from the plurality of ambiguous location determination solutions based on a discriminator function.
  • FIG. 1 graphically illustrates an example of three measurements with unknown user spatial coordinates and clock temporal bias value with two ambiguous location determination solutions.
  • FIG. 2 graphically illustrates an example of four measurements with three unknowns with two ambiguous location determination solutions.
  • FIG. 3 graphically illustrates an example of noisy measurements with three ambiguous location determination solutions.
  • FIG. 4 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Manhattan dataset.
  • FIG. 5 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
  • CDF cumulative distribution function
  • FIG. 6 graphically illustrates the histogram of the clock temporal bias produced by each of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) for the Manhattan dataset.
  • FIG. 7 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Campbell dataset.
  • FIG. 8 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Campbell dataset.
  • CDF cumulative distribution function
  • FIG. 9 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Japan dataset.
  • FIG. 10 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset.
  • CDF cumulative distribution function
  • FIG. 11 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) and the chosen final location determination solution by resolving ambiguity through consistency in order of ranging signal arrival and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
  • CDF cumulative distribution function
  • FIG. 12 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) and the chosen final location determination solution by resolving ambiguity through consistency in order of ranging signal arrival and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset.
  • CDF cumulative distribution function
  • FIG. 13 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) and the chosen final location determination solution by resolving ambiguity through consistency of distance values and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
  • CDF cumulative distribution function
  • FIG. 14 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) and the chosen final location determination solution by resolving ambiguity through consistency of distance values and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset.
  • CDF cumulative distribution function
  • FIG. 15 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) with the best ambiguous location determination solutions with and without noise removal and are plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
  • CDF cumulative distribution function
  • Various heuristic approaches for ambiguity resolution are disclosed that exploit known or expected behavior of one or more unknowns or combinations thereof, consistency of various quantifiable parameters, and/or redundancy from measurements or different parameters within a measurement.
  • a main advantage of the heuristic approaches disclosed is that redundant measurements are not required for the purpose of location determination solution disambiguation.
  • the Manhattan database consists of 1316 fixes from about 25 locations across the city, with about 20 to 100 position fixes per location.
  • the Japan dataset consists of about 1912 fixes from all over Japan.
  • the Campbell dataset consists of a stationary location with about 2000 fixes.
  • the method for ambiguity resolution uses physical constraints on the behavior of some of the unknowns and combinations thereof.
  • asynchronous systems can benefit from a set of assumptions or information concerning clock temporal bias value, which can be exploited to disambiguate the user location.
  • user clock temporal bias value may include the time the reference timing signal takes to propagate from the source of the reference timing signal to the user, some multipath and receiver processing overhead and clock slewing error.
  • PDF probability distribution function
  • User clock temporal bias values may be characterized by known statistical distributions, with distribution parameters dictated by various factors, such as the receiver design, the user clock design characteristics, multipath characteristics, and, in a terrestrial system, by network deployment issues such as base station density.
  • the distribution parameters may be set statically or learned dynamically for a given location system.
  • the goodness metric of a final location determination solution can be derived from the value of the computed clock temporal bias value and the probability distribution function (PDF) of clock temporal bias value for that scenario.
  • PDF probability distribution function
  • the following steps are needed: First, generate a probability distribution function (PDF) model for the N clock temporal bias values.
  • the PDF model is parameterized by its mean and standard deviation.
  • Fourth, for each of the ambiguous location determination solutions evaluate the PDF model at the inserted clock temporal bias value to get a PDF value. The result should be N PDF values for the N clock temporal bias values.
  • N goodness metrics corresponding to the N PDF values which correspond to the N clock temporal bias values.
  • Sixth compare the N goodness metrics and defined a maximum goodness metric as the largest of the N goodness metrics.
  • Seventh select the final location determination solution from the ambiguous location determination solution with the maximum goodness metric.
  • the quantity N is defined as an integer greater than one.
  • clock temporal bias value is modeled to follow a Gaussian probability distribution function (PDF) with a mean and standard deviation of minus 100 meters and 100 meters respectively. If the computed clock temporal bias value (multiplied by the speed of light) is minus 200 meters, then goodness metric can be set to the value of the Gaussian probability distribution function one standard deviation away from the mean.
  • PDFs other than a Gaussian PDF (for example, a finite set of statistical parameters) may be used.
  • clock temporal bias value is assumed to have a static, simple uniform distribution, with permissible clock temporal bias value between some value A and another value B.
  • a and B are set to minus 1000 meter and 1000 meters respectively.
  • a and B are set to minus 2000 meters and 100 meters to better suit the Japan dataset as observed empirically.
  • FIG. 4 the relationship between clock temporal bias value and the horizontal error of each ambiguous location determination solution is shown in FIG. 4 .
  • FIG. 5 For large errors, it is observed that there is a strong, generally linear dependence of error on clock temporal bias value.
  • the disambiguation algorithm works well in that the final location determination solution statistically outperforms either of the two ambiguous location determination solutions across the entire dataset. It is noted that the ambiguous location determination solution labeled “Solution 2 ” generally outperforms the ambiguous location determination solution labeled as “Solution 1 ”. Such discrepancy is apparent and explainable by the behavior of clock temporal bias value as produced by each ambiguous location determination solution, and shown in FIG. 6 . One can see how the second ambiguous location determination solution has a “tight” distribution while the first ambiguous location determination solution produced highly-varying results.
  • FIG. 7 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Campbell dataset.
  • FIG. 8 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) and the chosen final location determination solution plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Campbell dataset.
  • CDF cumulative distribution function
  • FIG. 9 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions
  • FIG. 10 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution # 2 ) and the chosen final location determination solution plotted as a cumulative distribution function (CDF) of the horizontal errors in meters.
  • CDF cumulative distribution function
  • ranging signal order There are various ways for deriving ranging signal order, such as the time of arrival (TOA) of ranging signals and/or the power level of ranging signals.
  • TOA time of arrival
  • the method for ambiguity resolution uses the order of receipt (a.k.a. time of arrival) of ranging signals from a plurality of sources to resolve ambiguity.
  • the expected order corresponding to each solution is derived and compared to the actual order in which the ranging signals are received.
  • the chosen final location determination solution is the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
  • the following steps are needed: First, rank the ranging signals based on the order of actual time of arrival (TOA) from earliest to latest. Second, for each of the ambiguous location determination solution, rank the ranging signals based on the expected order of time of arrival corresponding to each ambiguous location determination solution. Third, compare the ranking of the ranging signals in the order of the actual time of arrival (TOA) from earliest to latest to the ranking of the ranging signals based on the expected order of arrival corresponding to each ambiguous location determination solution. Fourth, select the final location determination solution as the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
  • the a-posteriori ranging signal order is compared to the a-priori ranging signal order derived from the time-of-arrival (TOA) information.
  • TOA time-of-arrival
  • For each ranging signal its rankings in the a-priori and a-posteriori lists are compared and labeled a “1” in case of a match and a “0” in case of a mismatch.
  • the solution with the most matches (i.e., the most “1”s) across all ranging signals is picked as the location determination solution. In case of a tie, the average is used.
  • the results for this disambiguation technique applied to the Manhattan dataset are shown in FIG. 11 and for the Japan dataset are shown in FIG. 12 .
  • the final location determination solution fares similarly or statistically outperforms either of the two ambiguous location determination solutions.
  • the method for ambiguity resolution uses the received power level of ranging signals to discriminate among ambiguous location determination solutions.
  • To determine the final location determination solution from a set of ambiguous location determination solutions using the received power levels of the ranging signals the following steps are needed: First, rank the ranging signals based on the order of the received power level received, from strongest to weakest. Second, for each of the ambiguous location determination solution, rank the ranging signals based on the expected order of the received power levels corresponding to each ambiguous location determination solution. The expected order of the received power levels is assumed to be the same as the expected order of the time of arrival corresponding to each ambiguous location determination solution.
  • the method for ambiguity resolution uses the numerical values of the distances to M sources and combinations thereof to resolve ambiguity.
  • the quantity of sources is defined as M whereby M is an integer greater than one.
  • the ambiguous location determination solution with the closest match (as defined by the root-mean-square of D 1 minus D 2 or by the mean square of D 1 minus D 2 ), or some combination of a number of values with the closest matches is chosen as the final location determination solution.
  • the following steps are needed: First, for each of the M sources, derive the distance D 1 corresponding to each ambiguous location determination solution. Second, obtain the original distance D 2 (which is the original distance adjusted by clock temporal bias value). Third, compare the distance D 1 with the original distance D 2 and compute their error measurement. In one embodiment, the error measurement is the root-mean-square of D 1 -D 2 . In another embodiment, the error measurement is the root-sum-square of D 1 -D 2 . And, in yet another embodiment, the error measurement is the root-mean-square of D 1 -D 2 normalized by the size of the vector D 1 (or D 2 ). Fourth, select the final location determination solution as having the lowest value of error measurement.
  • noisy measurements can also lead to errors in location determination.
  • a subset of Q signals is selected by some method so that the algorithm may run efficiently.
  • the selection method may select the Q signals randomly or via a more systematic method, such as choosing those signals with the lowest power level (that are more likely to incur noise overhead).
  • L noise levels are postulated.
  • the location determination solutions correspond to each of L times Q combinations of noise levels.
  • a method for noise disambiguation involves the following steps. First, for the first of the Q ranging signals, select L noise levels iteratively within a range [A, B] with an increment I.
  • One skilled in the art would know the range [A, B] and the increment I (and hence, the quantity of L) based on the chosen location determination system.
  • Second repeat the above step Q minus 1 times for each of the rest of the Q minus 1 ranging signals to produce a plurality of L noise levels.
  • a discriminator function known to one skilled in the art, is used to select the final location determination solution from the plurality of ambiguous location determination solutions produced for each of the plurality of noise levels.
  • the discriminator function can use one or a combination of goodness metrics such as power ranking, distance values, clock temporal bias constraints, order of receipt of the ranging signals. In one embodiment, the discriminator function jointly maximizes the consistency of power ranking and distance values, while meeting the clock temporal bias constraint.
  • Other discriminator functions are known to one skilled in the art and may be used without departing from the spirit of the present invention.
  • N, A, B are set to 2, 0, and 500 meters respectively.
  • the increment I values is investigated in the set of ⁇ 100 meters, 50 meters ⁇ .
  • Disambiguation among solutions is done by means of a combination of power ranking, range consistency and clock temporal bias value.
  • the power ranking algorithm works as follows: Using the power level of incoming ranging signals, the a-priori order of ranging signals is computed. Then, for each solution corresponding to a noise setting, the a-posteriori order of ranging signals is computed using geographical distances. The a-priori and a-posteriori orders are compared for a given ranging signal.
  • the power ranking is computed as the sum of matches across all ranging signals, weighted for each ranging signal in the following way: if the ranging signal is the first ranging signal in the a-priori list (that is, for the base station which is assumed closest), then the weight is 1 ⁇ 2. Otherwise, the weight is 0.5 times (the number of ranging signals ⁇ 1). Note the weights are normalized, as they add up to 1. Also, the bias in the weighting allows the strongest ranging signal to be more heavily weighted which has been shown to give improved results.
  • the root-mean-square differences in the a-priori and a-posteriori ranges are also computed for each solution. Range consistency is calculated as 1—RRMS/max(RRMS).
  • the clock temporal bias value enforces a [ ⁇ 1000 meters to 1000 meters] range on the accepted solutions.
  • the final location determination solution is the one that meets the clock temporal bias value constraint while maximizing the product of the power ranking and range consistency. The results are shown in FIG. 15 . It is noted that an improvement is observed over the best algebraic method.

Abstract

A method using various heuristics techniques for resolving ambiguity in location determination in environments with or without noise. A final location determination solution may be determined from a set of ambiguous location determination solutions by using clock temporal bias value, by using consistency information of ranging signal order such as the time of arrival and/or the received power level of the ranging signals, by using the distances to the sources, and/or by using other discriminator functions to select the final location determination solution from a plurality of ambiguous location determination solutions. The main advantage of the heuristic approaches is that redundant measurements are not required for location determination solution disambiguation.

Description

    CROSS-REFERENCES TO RELATED APPLICATIONS
  • This application is a divisional application of U.S. patent application Ser. No. 11/034,705, filed on Jan. 12, 2005, which claims benefit of priority from U.S. Provisional Patent Application No. 60/629,011, filed Nov. 17, 2004, each of which is incorporated by reference for all purposes.
  • BACKGROUND OF THE INVENTION
  • 1. Field
  • The disclosed subject matter relates generally to methods for location determination, and more particularly to methods using various heuristics techniques for resolving ambiguity in location determination in environments with or without noise.
  • 2. Background
  • In range-based location determination systems, measurements of ranging signals from a plurality of sources are converted to distance information associated with the source of each ranging signal. Distances to different sources with known locations are combined to solve for the unknown user location via geometric techniques known, for example as trilateration (a.k.a. triangulation). If delay of ranging signals cannot be known reliably (e.g. in asynchronous systems where the user clock is not synchronized to the network), location determination algorithms may treat user clock temporal bias as another unknown, to be solved for by the trilateration process, using an additional measurement.
  • However, a location determination system is ambiguous if more than one location determination solution set of user coordinates and clock temporal bias is consistent with a set of distance measurements. Location determination systems can produce ambiguous location determination solutions in three distinct ways: first, ambiguity can be caused by insufficient measurements; second, ambiguity may be introduced by the properties of the algorithm employed in location determination; and third, ambiguity may be introduced by the presence of noisy measurements.
  • First, a system has an insufficient number of measurements when the number of unknowns is greater than or equal to the number of independent measurements. For example, consider the case where the unknowns are the two-dimensional user spatial coordinates and user clock temporal bias. Consider the case depicted in FIG. 1. There are three unknowns, namely the mobile station latitude, longitude and clock temporal bias. There are three base stations, namely BS1, BS2 and BS3, and three associated distance measurements. Circles are plotted centered at a particular base station, with radii given by the sum of the distance between the mobile station and the base station as measured at the mobile station, and the computed clock temporal bias corresponding to a fitting solution. Given three independent distance measurements, there are two possible location determination solutions, depicted at the intersection of each set of circles.
  • Second, the nature of the algorithm used for locating a user can also be a source of ambiguity. A well known algorithm that is susceptible to ambiguity is described in the U.S. Pat. No. 6,289,280. This algorithm solves for unknowns using a closed form system of equations. Because it solves for the user location algebraically, this algorithm runs efficiently, making it suitable to applications and devices with time or resource constraints. The solution uses linear algebra manipulations to combine the measurements into a system of quadratic equations where the number of equations equals the number of unknowns. Two solutions are produced associated with the two roots of the quadratic equations. The two solutions form an ambiguous set of solutions which needs to be resolved by additional means.
  • For example, consider the case where the unknowns are the two-dimensional user spatial coordinates and user clock temporal bias. With four measurements, the system can be said to have a sufficient number of measurements to unambiguously solve for the user location. Yet, when the algebraic method is used, the four measurements are combined into three “average” measurements and two solutions corresponding to these averages are identified, as shown in FIG. 2.
  • Third, noisy measurements can lead to error in the determination of user location. Consider a method for location determination in noisy environments by assuming the noise to be a discrete variable with known or computable statistical parameters. A set of adjusted measurements and corresponding solution are generated for each assumed noise level. Such location determination system is ambiguous, thus also warranting ambiguity resolution techniques. For example, in FIG. 3, consider three noise levels, each 100 meters apart, associated with the measurement from base station BS2. For each noise level, a set of circles is plotted as before, with a radius corresponding to the sum of distance measurement (in the case of BS2, this measurement is adjusted by the assumed noise level) and the clock temporal bias computed. There are three ambiguous solutions, associated with each noise level, shown by the intersection of circles.
  • Accordingly, since more than one possible solution is presented by these prior art algorithms, it would be desirable to provide a method for selecting the correct (a.k.a. final) location determination solution from a set of ambiguous location determination solutions.
  • BRIEF SUMMARY OF THE INVENTION
  • Disclosed are methods for selecting the correct (a.k.a. final) location determination solution from a set of ambiguous location determination solutions using various heuristics and/or noise removal.
  • According to one aspect, a method for resolving ambiguity in location determination with N ambiguous location determination solutions using clock temporal bias value includes the following steps: generating a probability distribution function (PDF) model for clock temporal bias, obtaining N clock temporal bias values for the N ambiguous location determination solutions, inserting each of the N clock temporal bias values into the PDF model, evaluating the PDF model to get N PDF values, setting N goodness metrics to the N PDF values, comparing the N goodness metrics, defining a maximum goodness metric as the largest of the N goodness metrics, and selecting a final location determination solution having the maximum goodness metric.
  • According to another aspect, a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions using the order of receipt (a.k.a. time of arrival) of a plurality of ranging signals includes the following steps: ranking the plurality of ranging signals based on the order of actual time of arrival (TOA) from earliest to latest, ranking the plurality of ranging signals based on an expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions, and comparing the ranking of the ranging signals based on the order of the actual time of arrival (TOA) and based on the expected order of time of arrival corresponding to each of the plurality of ambiguous location determination solutions.
  • According to yet another one aspect, a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions using the order of received power levels of a plurality of ranging signals includes the following steps: ranking the plurality of ranging signals based on the order of received power levels from strongest to weakest, ranking the plurality of ranging signals based on an expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions, and comparing the ranking of the ranging signals based on the order of received power levels and based on the expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions.
  • According to still another aspect, a method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions includes the following steps: deriving a plurality of distances D1 corresponding to each of the plurality of ambiguous location determination solutions for each of M sources, obtaining a plurality of original distances D2, comparing each of the plurality of distances D1 with each of its corresponding plurality of original distances D2 and computing a plurality of error measurements, and selecting a final location determination solution as having the lowest value of the plurality of error measurements.
  • According to yet another aspect, a method for resolving ambiguity in location determination in a noisy environment includes the following steps: selecting L noise levels iteratively within a range [A, B] with an increment I for a first of Q ranging signals, repeating the selecting step Q minus 1 times for each of the rest of the Q ranging signals to produce a plurality of L noise levels, creating a plurality of ambiguous location determination solutions for each of the plurality of L noise levels based on a geometric technique, selecting the final location determination solution from the plurality of ambiguous location determination solutions based on a discriminator function.
  • It is understood that other embodiments will become readily apparent to those skilled in the art from the following detailed description, wherein it is shown and described various embodiments by way of illustration. The drawings and detailed description are to be regarded as illustrative in nature and not as restrictive.
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • FIG. 1 graphically illustrates an example of three measurements with unknown user spatial coordinates and clock temporal bias value with two ambiguous location determination solutions.
  • FIG. 2 graphically illustrates an example of four measurements with three unknowns with two ambiguous location determination solutions.
  • FIG. 3 graphically illustrates an example of noisy measurements with three ambiguous location determination solutions.
  • FIG. 4 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Manhattan dataset.
  • FIG. 5 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
  • FIG. 6 graphically illustrates the histogram of the clock temporal bias produced by each of the two ambiguous location determination solutions (solution # 1 and solution #2) for the Manhattan dataset.
  • FIG. 7 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Campbell dataset.
  • FIG. 8 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Campbell dataset.
  • FIG. 9 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Japan dataset.
  • FIG. 10 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through clock temporal bias value constraints and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset.
  • FIG. 11 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency in order of ranging signal arrival and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
  • FIG. 12 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency in order of ranging signal arrival and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset.
  • FIG. 13 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency of distance values and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
  • FIG. 14 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution by resolving ambiguity through consistency of distance values and is plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Japan dataset.
  • FIG. 15 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) with the best ambiguous location determination solutions with and without noise removal and are plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Manhattan dataset.
  • DETAILED DESCRIPTION OF THE INVENTION
  • The description set forth below in connection with the appended drawings is intended as a description of various embodiments of the invention and is not intended to represent the only embodiments in which the invention may be practiced. Each embodiment is provided merely as an example or illustration, and should not necessarily be construed as preferred or advantageous over other embodiments. Specific details are used to provide an understanding of the invention. However, it will be apparent to those skilled in the art that the invention may be practiced without these specific details. Acronyms and other descriptive terminology may be used merely for convenience and clarity and are not intended to limit the scope of the invention.
  • Various heuristic approaches for ambiguity resolution are disclosed that exploit known or expected behavior of one or more unknowns or combinations thereof, consistency of various quantifiable parameters, and/or redundancy from measurements or different parameters within a measurement. A main advantage of the heuristic approaches disclosed is that redundant measurements are not required for the purpose of location determination solution disambiguation.
  • For the purpose of illustration only, simulations using real user data from Manhattan, Japan, and Campbell are presented. The Manhattan database consists of 1316 fixes from about 25 locations across the city, with about 20 to 100 position fixes per location. The Japan dataset consists of about 1912 fixes from all over Japan. The Campbell dataset consists of a stationary location with about 2000 fixes.
  • In one embodiment, the method for ambiguity resolution uses physical constraints on the behavior of some of the unknowns and combinations thereof. In particular, asynchronous systems can benefit from a set of assumptions or information concerning clock temporal bias value, which can be exploited to disambiguate the user location. In a system where the user clock timing is learned from some reference timing signal, user clock temporal bias value may include the time the reference timing signal takes to propagate from the source of the reference timing signal to the user, some multipath and receiver processing overhead and clock slewing error. A probability distribution function (PDF) of clock temporal bias value can be derived for that network and used to disambiguate the user location.
  • User clock temporal bias values may be characterized by known statistical distributions, with distribution parameters dictated by various factors, such as the receiver design, the user clock design characteristics, multipath characteristics, and, in a terrestrial system, by network deployment issues such as base station density. The distribution parameters may be set statically or learned dynamically for a given location system. The goodness metric of a final location determination solution can be derived from the value of the computed clock temporal bias value and the probability distribution function (PDF) of clock temporal bias value for that scenario.
  • To determine the final location determination solution from a set of N ambiguous location determination solutions using N clock temporal bias values, the following steps are needed: First, generate a probability distribution function (PDF) model for the N clock temporal bias values. In one embodiment, the PDF model is parameterized by its mean and standard deviation. Second, for each of the N ambiguous location determination solutions, obtain a clock temporal bias value, resulting in N clock temporal bias value. Third, for each of the ambiguous location determination solutions, insert its corresponding clock temporal bias value into the PDF model. Fourth, for each of the ambiguous location determination solutions, evaluate the PDF model at the inserted clock temporal bias value to get a PDF value. The result should be N PDF values for the N clock temporal bias values. Fifth, for each of the ambiguous location determination solutions, set the goodness metric of the ambiguous location determination solution to the PDF value. Again, there should be N goodness metrics corresponding to the N PDF values which correspond to the N clock temporal bias values. Sixth, compare the N goodness metrics and defined a maximum goodness metric as the largest of the N goodness metrics. Seventh, select the final location determination solution from the ambiguous location determination solution with the maximum goodness metric. The quantity N is defined as an integer greater than one.
  • In one embodiment, clock temporal bias value is modeled to follow a Gaussian probability distribution function (PDF) with a mean and standard deviation of minus 100 meters and 100 meters respectively. If the computed clock temporal bias value (multiplied by the speed of light) is minus 200 meters, then goodness metric can be set to the value of the Gaussian probability distribution function one standard deviation away from the mean. One skilled in the art would know that other PDFs other than a Gaussian PDF (for example, a finite set of statistical parameters) may be used.
  • In an example, clock temporal bias value is assumed to have a static, simple uniform distribution, with permissible clock temporal bias value between some value A and another value B. For Manhattan and Campbell datasets, A and B are set to minus 1000 meter and 1000 meters respectively. In the case of Japan, A and B are set to minus 2000 meters and 100 meters to better suit the Japan dataset as observed empirically. For the Manhattan dataset, the relationship between clock temporal bias value and the horizontal error of each ambiguous location determination solution is shown in FIG. 4. For large errors, it is observed that there is a strong, generally linear dependence of error on clock temporal bias value. For the Manhattan dataset, the performance results are shown in FIG. 5. It is observed that the disambiguation algorithm works well in that the final location determination solution statistically outperforms either of the two ambiguous location determination solutions across the entire dataset. It is noted that the ambiguous location determination solution labeled “Solution 2” generally outperforms the ambiguous location determination solution labeled as “Solution 1”. Such discrepancy is apparent and explainable by the behavior of clock temporal bias value as produced by each ambiguous location determination solution, and shown in FIG. 6. One can see how the second ambiguous location determination solution has a “tight” distribution while the first ambiguous location determination solution produced highly-varying results.
  • FIG. 7 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions for the Campbell dataset. FIG. 8 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution plotted as a cumulative distribution function (CDF) of the horizontal errors in meters for the Campbell dataset. For the Japan dataset, FIG. 9 graphically illustrates the relationship between the clock temporal bias value and the horizontal error for each of the ambiguous location determination solutions while FIG. 10 graphically illustrates the statistical performance of the two ambiguous location determination solutions (solution # 1 and solution #2) and the chosen final location determination solution plotted as a cumulative distribution function (CDF) of the horizontal errors in meters.
  • Alternative methods for ambiguity resolution use consistency of ranging signal order. There are various ways for deriving ranging signal order, such as the time of arrival (TOA) of ranging signals and/or the power level of ranging signals.
  • In one embodiment, the method for ambiguity resolution uses the order of receipt (a.k.a. time of arrival) of ranging signals from a plurality of sources to resolve ambiguity. The expected order corresponding to each solution is derived and compared to the actual order in which the ranging signals are received. The chosen final location determination solution is the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
  • To determine the final location determination solution from a set of ambiguous location determination solutions using the order of receipt (a.k.a. time of arrival) of ranging signals, the following steps are needed: First, rank the ranging signals based on the order of actual time of arrival (TOA) from earliest to latest. Second, for each of the ambiguous location determination solution, rank the ranging signals based on the expected order of time of arrival corresponding to each ambiguous location determination solution. Third, compare the ranking of the ranging signals in the order of the actual time of arrival (TOA) from earliest to latest to the ranking of the ranging signals based on the expected order of arrival corresponding to each ambiguous location determination solution. Fourth, select the final location determination solution as the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
  • In an example, the a-posteriori ranging signal order is compared to the a-priori ranging signal order derived from the time-of-arrival (TOA) information. For each ranging signal, its rankings in the a-priori and a-posteriori lists are compared and labeled a “1” in case of a match and a “0” in case of a mismatch. The solution with the most matches (i.e., the most “1”s) across all ranging signals is picked as the location determination solution. In case of a tie, the average is used. The results for this disambiguation technique applied to the Manhattan dataset are shown in FIG. 11 and for the Japan dataset are shown in FIG. 12. In general, the final location determination solution fares similarly or statistically outperforms either of the two ambiguous location determination solutions.
  • In another embodiment, the method for ambiguity resolution uses the received power level of ranging signals to discriminate among ambiguous location determination solutions. One can use received power level to discriminate among location determination solutions. To determine the final location determination solution from a set of ambiguous location determination solutions using the received power levels of the ranging signals, the following steps are needed: First, rank the ranging signals based on the order of the received power level received, from strongest to weakest. Second, for each of the ambiguous location determination solution, rank the ranging signals based on the expected order of the received power levels corresponding to each ambiguous location determination solution. The expected order of the received power levels is assumed to be the same as the expected order of the time of arrival corresponding to each ambiguous location determination solution. Third, compare the ranking of the ranging signals in the order of the received power level received from strongest to weakest to the ranking of the ranging signals based on the expected order of the received power levels corresponding to each ambiguous location determination solution. Fourth, select the final location determination solution as the one with the closest match in ordering, or a combination of a number of ambiguous location determination solutions with the closest match in ordering.
  • In another embodiment, the method for ambiguity resolution uses the numerical values of the distances to M sources and combinations thereof to resolve ambiguity. One can derive each of the distance D1 to each of the M sources corresponding to each ambiguous location determination solution and compare them to the original distance adjusted by the clock temporal bias value (labeled as original distance D2). The quantity of sources is defined as M whereby M is an integer greater than one. The ambiguous location determination solution with the closest match (as defined by the root-mean-square of D1 minus D2 or by the mean square of D1 minus D2), or some combination of a number of values with the closest matches is chosen as the final location determination solution. It will be understood by one skilled in the art that the present invention is not confined to the root-mean-square or mean-square of D1 minus D2 and that other error measures such as, but not limited to, root-sum-square (RSS) or sum-square (SS) may also be used without violating the spirit of the invention. The results for the Manhattan and Japan dataset are shown in FIG. 13 and FIG. 14 respectively. In both cases, the final location determination solution statistically outperforms either of the two ambiguous location determination solutions across the entire dataset.
  • To determine the final location determination solution from a set of ambiguous location determination solutions using numerical values of the distances to the M sources, the following steps are needed: First, for each of the M sources, derive the distance D1 corresponding to each ambiguous location determination solution. Second, obtain the original distance D2 (which is the original distance adjusted by clock temporal bias value). Third, compare the distance D1 with the original distance D2 and compute their error measurement. In one embodiment, the error measurement is the root-mean-square of D1-D2. In another embodiment, the error measurement is the root-sum-square of D1-D2. And, in yet another embodiment, the error measurement is the root-mean-square of D1-D2 normalized by the size of the vector D1 (or D2). Fourth, select the final location determination solution as having the lowest value of error measurement.
  • Noisy measurements can also lead to errors in location determination. Out of a set of P signals, a subset of Q signals is selected by some method so that the algorithm may run efficiently. The selection method may select the Q signals randomly or via a more systematic method, such as choosing those signals with the lowest power level (that are more likely to incur noise overhead). For each of the Q signals selected, L noise levels are postulated. The location determination solutions correspond to each of L times Q combinations of noise levels. A method for noise disambiguation involves the following steps. First, for the first of the Q ranging signals, select L noise levels iteratively within a range [A, B] with an increment I. One skilled in the art would know the range [A, B] and the increment I (and hence, the quantity of L) based on the chosen location determination system. Second, repeat the above step Q minus 1 times for each of the rest of the Q minus 1 ranging signals to produce a plurality of L noise levels. Third, create a plurality of ambiguous location determination solutions for each of the plurality of noise levels based on known geometric techniques such as trilateration. Fourth, a discriminator function, known to one skilled in the art, is used to select the final location determination solution from the plurality of ambiguous location determination solutions produced for each of the plurality of noise levels. The discriminator function can use one or a combination of goodness metrics such as power ranking, distance values, clock temporal bias constraints, order of receipt of the ranging signals. In one embodiment, the discriminator function jointly maximizes the consistency of power ranking and distance values, while meeting the clock temporal bias constraint. Other discriminator functions are known to one skilled in the art and may be used without departing from the spirit of the present invention.
  • In an example, the noise removal techniques for fixes in the Manhattan dataset. N, A, B are set to 2, 0, and 500 meters respectively. The increment I values is investigated in the set of {100 meters, 50 meters}. Disambiguation among solutions is done by means of a combination of power ranking, range consistency and clock temporal bias value. The power ranking algorithm works as follows: Using the power level of incoming ranging signals, the a-priori order of ranging signals is computed. Then, for each solution corresponding to a noise setting, the a-posteriori order of ranging signals is computed using geographical distances. The a-priori and a-posteriori orders are compared for a given ranging signal. Its ranking in the a-priori and a-posteriori lists are compared and a match is labeled a “1” while a mismatch returns a “0.” The power ranking is computed as the sum of matches across all ranging signals, weighted for each ranging signal in the following way: if the ranging signal is the first ranging signal in the a-priori list (that is, for the base station which is assumed closest), then the weight is ½. Otherwise, the weight is 0.5 times (the number of ranging signals −1). Note the weights are normalized, as they add up to 1. Also, the bias in the weighting allows the strongest ranging signal to be more heavily weighted which has been shown to give improved results. The root-mean-square differences in the a-priori and a-posteriori ranges (RRMS) are also computed for each solution. Range consistency is calculated as 1—RRMS/max(RRMS). The clock temporal bias value enforces a [−1000 meters to 1000 meters] range on the accepted solutions. The final location determination solution is the one that meets the clock temporal bias value constraint while maximizing the product of the power ranking and range consistency. The results are shown in FIG. 15. It is noted that an improvement is observed over the best algebraic method.
  • The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. 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.

Claims (14)

1. A method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions for a mobile station using an order of received power levels of a plurality of ranging signals from base stations in a cellular network, the method comprising:
receiving the plurality of ranging signals at the mobile station;
ranking the plurality of ranging signals based on the order of received power levels from strongest to weakest;
ranking the plurality of ranging signals based on an expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions;
comparing the ranking of the ranging signals based on the order of received power levels and based on the expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions; and
determining a final location determination solution based on a result of comparing the rankings of the ranging signals.
2. The method of claim 1, wherein determining the final location determination solution comprises determining the ambiguous location determination solution, from the plurality of ambiguous location determination solutions, for which the expected order of received power levels is a closest match to the order of received power levels.
3. The method of claim 1, wherein determining the final location determination solution comprises determining a combination of the plurality of ambiguous location determination solutions with a closest match in ordering.
4. A method of determining a final location determination solution for a mobile station from a plurality of candidate location determination solutions, the method comprising:
receiving ranging signals from a plurality of sources at the mobile station;
determining a received power level for each of the ranging signals;
ranking the ranging signals based on the received power level for each of the ranging signals;
determining an expected order of received power levels for each of the candidate location determination solutions;
determining a preferred candidate location determination solution from the plurality of candidate location determination solutions for which the expected order of received power levels is a closest match to the ranking of the ranging signals based on the received power level; and
determining the final location determination solution for the mobile station based on the preferred candidate location determination solution.
5. The method of claim 4, further comprising:
summing a number of matches between the ranging signals ranked based on the received power levels and the ranging signals ranked based on the expected order of the received power levels for each of the candidate location determination solutions.
6. The method of claim 5, wherein summing the number of matches comprises weighting each ranging signal according to a corresponding position in the expected order of the received power levels.
7. The method of claim 5, wherein the preferred candidate location determination solution has a greatest number of the matches among the plurality of candidate location determination solutions.
8. The method of claim 4, wherein each of the plurality of candidate location determination solutions corresponds to a distinct noise level.
9. A method for resolving ambiguity in location determination with a plurality of ambiguous location determination solutions for a mobile station using an order of received power levels of a plurality of ranging signals transmitted by base stations in a cellular network, the method comprising:
receiving the ranging signals at the mobile station from the base stations;
determining a received power level for each of the ranging signals;
ranking the ranging signals based on the order of received power levels from strongest to weakest;
determining an expected order of the received power levels for each of the ambiguous location determination solutions;
assigning a score to each of the ambiguous location determination solutions based on a number of matches between the ranging signals in the expected order of the received power levels corresponding to each of the ambiguous location determination solutions and the ranging signals ranked based on the order of the received power levels; and
determining a final location determination solution for the mobile stations based on the scores assigned to the ambiguous location determination solutions.
10. The method of claim 9, wherein assigning the score to each of the ambiguous location determination solutions comprises weighting each of the ranging signals according to a corresponding position in a corresponding one of the expected orders of the received power levels.
11. The method of claim 9, wherein each of the plurality of ambiguous location determination solutions corresponds to a distinct noise level.
12. A computer-readable medium having stored thereon a processor-readable program of instructions for causing a processor to resolve ambiguity in location determination with a plurality of ambiguous location determination solutions for a mobile station using an order of received power levels of a plurality of ranging signals, the program of instructions being configured to cause the processor to:
form a first ranking of the plurality of ranging signals based on the order of received power levels from strongest to weakest;
form a second ranking of the plurality of ranging signals based on an expected order of received power levels corresponding to each of the plurality of ambiguous location determination solutions;
compare the first ranking with each of the second rankings; and
determine a final location determination solution for the mobile station from the plurality of ambiguous location determination solutions based on a result of comparing the first ranking with the second rankings.
13. The computer readable medium of claim 13, wherein the program further comprises instructions configured to cause the processor to select as the final location determination solution an ambiguous location determination solution for which the expected order of the received power levels is a closest match to the order of the received power levels.
14. The computer readable medium of claim 13, wherein the program further comprises instructions configured to cause the processor to select as the final location determination solution a combination of the plurality of ambiguous location determination solutions with a closest match in ordering.
US12/754,360 2004-11-17 2010-04-05 Method for ambiguity resolution in location determination Abandoned US20100190512A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US12/754,360 US20100190512A1 (en) 2004-11-17 2010-04-05 Method for ambiguity resolution in location determination

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US62901104P 2004-11-17 2004-11-17
US11/034,705 US7881905B2 (en) 2004-11-17 2005-01-12 Method for ambiguity resolution in location determination
US12/754,360 US20100190512A1 (en) 2004-11-17 2010-04-05 Method for ambiguity resolution in location determination

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
US11/034,705 Division US7881905B2 (en) 2004-11-17 2005-01-12 Method for ambiguity resolution in location determination

Publications (1)

Publication Number Publication Date
US20100190512A1 true US20100190512A1 (en) 2010-07-29

Family

ID=35892584

Family Applications (3)

Application Number Title Priority Date Filing Date
US11/034,705 Expired - Fee Related US7881905B2 (en) 2004-11-17 2005-01-12 Method for ambiguity resolution in location determination
US12/754,348 Abandoned US20100190511A1 (en) 2004-11-17 2010-04-05 Method for ambiguity resolution in location determination
US12/754,360 Abandoned US20100190512A1 (en) 2004-11-17 2010-04-05 Method for ambiguity resolution in location determination

Family Applications Before (2)

Application Number Title Priority Date Filing Date
US11/034,705 Expired - Fee Related US7881905B2 (en) 2004-11-17 2005-01-12 Method for ambiguity resolution in location determination
US12/754,348 Abandoned US20100190511A1 (en) 2004-11-17 2010-04-05 Method for ambiguity resolution in location determination

Country Status (6)

Country Link
US (3) US7881905B2 (en)
EP (3) EP2211196A3 (en)
JP (2) JP5065040B2 (en)
CA (1) CA2587649A1 (en)
RU (1) RU2384859C2 (en)
WO (1) WO2006055790A2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100190511A1 (en) * 2004-11-17 2010-07-29 Qualcomm Incorporated Method for ambiguity resolution in location determination

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8626503B2 (en) 2005-07-14 2014-01-07 Erik Gosuinus Petrus Schuijers Audio encoding and decoding
EP2283641B1 (en) * 2008-06-06 2020-08-12 Skyhook Wireless, Inc. Method and system for determining location using a hybrid satellite and wlan positioning system by selecting the best wlan-ps solution
CN102368234A (en) * 2011-09-27 2012-03-07 邵阳神风动力制造有限责任公司 Method for converting PDF (portable document format) documents into DWG (drawing) documents
KR102051965B1 (en) * 2012-06-28 2019-12-04 한국전자통신연구원 Apparatus and method for measuring position
WO2014151118A1 (en) * 2013-03-15 2014-09-25 Nextnav, Llc Techniques to improve the performance of a fixed, timing-based radio positioning network using external assistance information
CN108352912B (en) * 2015-11-09 2020-07-07 智加系统公司 Method for synchronizing a plurality of devices and determining a position based on synchronized devices
CN110139303B (en) * 2019-04-23 2021-06-11 四川九洲电器集团有限责任公司 Rapid simulation method and device for equivalent signal level TOA measurement
JP7067541B2 (en) 2019-10-31 2022-05-16 Jfeスチール株式会社 Rolling mill control method and control device

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5293642A (en) * 1990-12-19 1994-03-08 Northern Telecom Limited Method of locating a mobile station
US5623414A (en) * 1995-01-24 1997-04-22 Massachusetts Inst Technology Clock-aided satellite navigation receiver system for enhanced position estimation and integrity monitoring
US5740048A (en) * 1992-08-07 1998-04-14 Abel; Jonathan S. Method and apparatus for GPS positioning, filtering and integration
US5878033A (en) * 1994-03-22 1999-03-02 Nortel Marta Cellular Process and equipment for broadcasting messages to mobile radio communication stations
US6114992A (en) * 1997-05-22 2000-09-05 Conexant Systems, Inc. Satellite acquisition and measurement system and process
US6249252B1 (en) * 1996-09-09 2001-06-19 Tracbeam Llc Wireless location using multiple location estimators
US6263208B1 (en) * 1999-05-28 2001-07-17 Lucent Technologies Inc. Geolocation estimation method for CDMA terminals based on pilot strength measurements
US6266537B1 (en) * 1998-03-27 2001-07-24 Nec Corporation Radio communication system
US20010022558A1 (en) * 1996-09-09 2001-09-20 Tracbeam Llc Wireless location using signal fingerprinting
US20020135462A1 (en) * 2001-03-26 2002-09-26 Nels Michael J. Method for proactive personalization and alerts to mobile devices
US6496489B1 (en) * 1993-11-01 2002-12-17 Telefonaktiebolaget Lm Ericsson Method and apparatus for distinguishing between a digital control channel and a digital traffic channel in a radiocommunication system
US6496701B1 (en) * 1998-08-25 2002-12-17 Lucent Technologies Inc. Pattern-recognition-based geolocation
US20030040300A1 (en) * 2001-08-27 2003-02-27 Alcatel System of interoperability between MMS messages and SMS/EMS messages and an associated exchange method
US20030146871A1 (en) * 1998-11-24 2003-08-07 Tracbeam Llc Wireless location using signal direction and time difference of arrival
US20100190511A1 (en) * 2004-11-17 2010-07-29 Qualcomm Incorporated Method for ambiguity resolution in location determination

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU1309747A1 (en) 1985-05-06 1999-09-10 Г.Н. Громов SYSTEM OF DETERMINATION OF THE LOCATION OF THE AIRCRAFT
US4963889A (en) 1989-09-26 1990-10-16 Magnavox Government And Industrial Electronics Company Method and apparatus for precision attitude determination and kinematic positioning
RU2060513C1 (en) 1992-10-06 1996-05-20 Олег Евгениевич Антонов Method for identifying object location
RU2098297C1 (en) 1996-01-05 1997-12-10 Герман Федорович Иоилев Method of and system for detection and identification of object
AU3119300A (en) * 1999-01-08 2000-07-24 Trueposition, Inc. Bandwidth synthesis for wireless location system
CA2296812A1 (en) * 1999-02-17 2000-08-17 Lucent Technologies Inc. Method for combining multiple measurements to determine the position of a mobile transceiver
JP2000244967A (en) * 1999-02-24 2000-09-08 Mitsubishi Electric Corp Mobile communication system, mobile unit and base station configuring the system and method for detecting position of the mobile unit in the system
GB2347571A (en) * 1999-03-03 2000-09-06 Secr Defence Locating system
US6564065B1 (en) * 1999-04-20 2003-05-13 Lucent Technologies Inc. Bayesian-update based location prediction method for CDMA systems
US6889053B1 (en) * 1999-07-26 2005-05-03 Lucent Technologies Inc. Likelihood-based geolocation prediction algorithms for CDMA systems using pilot strength measurements
RU2152050C1 (en) 1999-12-01 2000-06-27 Закрытое акционерное общество "НПО Космического Приборостроения" Satellite navigation system for detection of object position
US6289280B1 (en) * 1999-12-10 2001-09-11 Qualcomm Incorporated Method and apparatus for determining an algebraic solution to GPS terrestrial hybrid location system equations
EP1252533A2 (en) * 2000-01-24 2002-10-30 Integrinautics Corportion Multi-frequency pseudolites for carrier-based differential-position determination
RU2161317C1 (en) 2000-05-11 2000-12-27 Закрытое акционерное общество "НПО Космического Приборостроения" System for high-accuracy localization of objects- navigational information users by navigational radio signals with authorized access in regime of differential corrections
JP4396801B2 (en) * 2001-02-27 2010-01-13 パイオニア株式会社 POSITIONING DEVICE, POSITIONING METHOD, AND POSITIONING SYSTEM
JP2003075524A (en) * 2001-08-31 2003-03-12 Denso Corp Mobile communication terminal device, program and position server therefor
US7006834B2 (en) * 2001-10-29 2006-02-28 Qualcomm Incorporated Base station time calibration using position measurement data sent by mobile stations during regular position location sessions
AU2002360755A1 (en) 2001-12-27 2003-07-24 Qualcomm, Incorporated Creating and using base station almanac information in a wireless communication system having a position location capability
FI113092B (en) * 2002-05-31 2004-02-27 Ekahau Oy Measures of position differences and applications
US7257411B2 (en) * 2002-12-27 2007-08-14 Ntt Docomo, Inc. Selective fusion location estimation (SELFLOC) for wireless access technologies
EP1445970B1 (en) 2003-02-05 2009-04-01 Cambridge Positioning Systems Limited A method and system for locating a mobile radio receiver in a radio system with multiple tranmitters
JP4374021B2 (en) * 2003-05-23 2009-12-02 シンボル テクノロジーズ インコーポレイテッド Self-correction method of position search system by signal strength

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5293642A (en) * 1990-12-19 1994-03-08 Northern Telecom Limited Method of locating a mobile station
US5740048A (en) * 1992-08-07 1998-04-14 Abel; Jonathan S. Method and apparatus for GPS positioning, filtering and integration
US6496489B1 (en) * 1993-11-01 2002-12-17 Telefonaktiebolaget Lm Ericsson Method and apparatus for distinguishing between a digital control channel and a digital traffic channel in a radiocommunication system
US5878033A (en) * 1994-03-22 1999-03-02 Nortel Marta Cellular Process and equipment for broadcasting messages to mobile radio communication stations
US5623414A (en) * 1995-01-24 1997-04-22 Massachusetts Inst Technology Clock-aided satellite navigation receiver system for enhanced position estimation and integrity monitoring
US20010022558A1 (en) * 1996-09-09 2001-09-20 Tracbeam Llc Wireless location using signal fingerprinting
US6249252B1 (en) * 1996-09-09 2001-06-19 Tracbeam Llc Wireless location using multiple location estimators
US6114992A (en) * 1997-05-22 2000-09-05 Conexant Systems, Inc. Satellite acquisition and measurement system and process
US6266537B1 (en) * 1998-03-27 2001-07-24 Nec Corporation Radio communication system
US6496701B1 (en) * 1998-08-25 2002-12-17 Lucent Technologies Inc. Pattern-recognition-based geolocation
US20030146871A1 (en) * 1998-11-24 2003-08-07 Tracbeam Llc Wireless location using signal direction and time difference of arrival
US6263208B1 (en) * 1999-05-28 2001-07-17 Lucent Technologies Inc. Geolocation estimation method for CDMA terminals based on pilot strength measurements
US20020135462A1 (en) * 2001-03-26 2002-09-26 Nels Michael J. Method for proactive personalization and alerts to mobile devices
US20030040300A1 (en) * 2001-08-27 2003-02-27 Alcatel System of interoperability between MMS messages and SMS/EMS messages and an associated exchange method
US20100190511A1 (en) * 2004-11-17 2010-07-29 Qualcomm Incorporated Method for ambiguity resolution in location determination
US7881905B2 (en) * 2004-11-17 2011-02-01 Qualcomm Incorporated Method for ambiguity resolution in location determination

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100190511A1 (en) * 2004-11-17 2010-07-29 Qualcomm Incorporated Method for ambiguity resolution in location determination

Also Published As

Publication number Publication date
EP2211197A2 (en) 2010-07-28
RU2007122504A (en) 2008-12-27
CA2587649A1 (en) 2006-05-26
JP2008521009A (en) 2008-06-19
US20100190511A1 (en) 2010-07-29
WO2006055790A2 (en) 2006-05-26
US7881905B2 (en) 2011-02-01
US20060106573A1 (en) 2006-05-18
JP2011107158A (en) 2011-06-02
JP5065040B2 (en) 2012-10-31
EP2211197A3 (en) 2010-12-15
WO2006055790A3 (en) 2007-01-04
RU2384859C2 (en) 2010-03-20
EP2211196A3 (en) 2010-12-15
EP2211196A2 (en) 2010-07-28
EP1812806A2 (en) 2007-08-01
JP5129356B2 (en) 2013-01-30

Similar Documents

Publication Publication Date Title
US20100190512A1 (en) Method for ambiguity resolution in location determination
DK1514130T3 (en) Probabilistic model for a positioning technique
McGuire et al. Data fusion of power and time measurements for mobile terminal location
Shah et al. 3D weighted centroid algorithm & RSSI ranging model strategy for node localization in WSN based on smart devices
CN101536591B (en) Extended clustering for improved positioning
Tayebi et al. The application of ray-tracing to mobile localization using the direction of arrival and received signal strength in multipath indoor environments
Kul et al. IEEE 802.11 WLAN based real time indoor positioning: Literature survey and experimental investigations
CN102695269B (en) A kind of positioning correction method, relevant apparatus and system
EP2384549A1 (en) Method for position estimation using generalized error distributions
CN108882149B (en) NLOS compensation positioning method of distance correlation probability
US20100141529A1 (en) Set mode passive location in toa/tdoa modes
CN102325370A (en) High-precision three-dimensional positioner for wireless sensor network node
KR100714120B1 (en) Method and apparatus for locating a remote unit
Li et al. Cramer-rao lower bound analysis of data fusion for fingerprinting localization in non-line-of-sight environments
KR101258832B1 (en) Method for determining location using the access point and apparatus thereof
KR102118981B1 (en) Indoor positioning method based on beacon signal and fingerprint map and system having the method
Pradhan et al. Mathematical analysis of line intersection algorithm for TOA trilateration method
Zhao et al. Signal path reckoning localization method in multipath environment
Filípek et al. Indoor localization based on beacons and calculated by particle filter
Go et al. An efficient non-line-of-sight error mitigation method for TOA measurement in indoor environments
Zhang et al. Improved KNN algorithm with historical information fusion for indoor positioning
Başak et al. Accurate indoor localization with optimized fingerprinting algorithm
Wu et al. A Regional Positioning Algorithm for Non-line-of-sight (NOLS) Propagation
Shixun et al. Single base station hybrid TOA/AOD/AOA localization algorithms with the synchronization error in dense multipath environment
Konchenko et al. Multilateration surveillance system arrangement at the airport area

Legal Events

Date Code Title Description
AS Assignment

Owner name: QUALCOMM INCORPORATED, CALIFORNIA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SEIBERT, CRISTINA A;REEL/FRAME:024236/0904

Effective date: 20100414

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO PAY ISSUE FEE