WO2001065271A1 - Methods and apparatus to position a mobile receiver using downlink signals - Google Patents

Methods and apparatus to position a mobile receiver using downlink signals Download PDF

Info

Publication number
WO2001065271A1
WO2001065271A1 PCT/CA2000/000224 CA0000224W WO0165271A1 WO 2001065271 A1 WO2001065271 A1 WO 2001065271A1 CA 0000224 W CA0000224 W CA 0000224W WO 0165271 A1 WO0165271 A1 WO 0165271A1
Authority
WO
WIPO (PCT)
Prior art keywords
location
bsts
aoa
signal
estimating
Prior art date
Application number
PCT/CA2000/000224
Other languages
French (fr)
Inventor
Michel Fattouche
Richard Klukas
Kevin Oler
Original Assignee
Cell-Loc 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
Priority claimed from US09/169,730 external-priority patent/US6204812B1/en
Priority claimed from US09/169,916 external-priority patent/US6208297B1/en
Priority claimed from US09/169,690 external-priority patent/US6266014B1/en
Application filed by Cell-Loc Inc. filed Critical Cell-Loc Inc.
Priority to AU29000/00A priority Critical patent/AU2900000A/en
Publication of WO2001065271A1 publication Critical patent/WO2001065271A1/en

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/0009—Transmission of position information to remote stations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-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/0246—Position-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 involving frequency difference of arrival or Doppler measurements
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/14—Relay systems
    • H04B7/15—Active relay systems
    • H04B7/185—Space-based or airborne stations; Stations for satellite systems
    • H04B7/1853—Satellite systems for providing telephony service to a mobile station, i.e. mobile satellite service
    • H04B7/18539—Arrangements for managing radio, resources, i.e. for establishing or releasing a connection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S2205/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S2205/001—Transmission of position information to remote stations
    • G01S2205/008—Transmission of position information to remote stations using a mobile telephone network
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-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/08—Position of single direction-finder fixed by determining direction of a plurality of spaced sources of known location
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-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/10—Position of receiver fixed by co-ordinating a plurality of position lines defined by path-difference measurements, e.g. omega or decca systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-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/12—Position-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 by co-ordinating position lines of different shape, e.g. hyperbolic, circular, elliptical or radial
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02—Position-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/14—Determining absolute distances from a plurality of spaced points of known location

Definitions

  • This invention relates to location finding and tracking of a Mobile Receiver (MR) using a Reference Receiver (RR) of known location.
  • radio signals that exist with worldwide coverage are satellite signals such as Global Positioning System (GPS), Iridium, Globalstar, etc.
  • GPS Global Positioning System
  • Examples of radio signals that can penetrate structures such as buildings are land-based systems such as cellular, PCS, paging, and ESMR, among others. Both satellite-based signals and terrestrial-based signals can be used to locate the MR by relying on their reception at both the MR and the RR.
  • a location system such as a handset-based Wireless Location System (WLS) (or Forward/downlink-based WLS), as opposed to a network-based WLS (or Reverse/uplink-based WLS), which locates a Mobile Transmitter (MT) using a plurality of receivers that receive the signal transmitted from the MT.
  • WLS Wireless Location System
  • MT Mobile Transmitter
  • An advantage of a network-based WLS is that it can be passive (in that it can locate MR's without modification or addition to wireless transmissions) and that no modifications to the hardware or operation of the MT are required for the WLS to locate the MT.
  • a handset-based WLS has the advantage of requiring no modification of existing network infrastructure.
  • the location in the handset-based WLS can be performed either at the MR, or at the RR, or at a Third Receiver (TR).
  • the RR is required to transmit the Reference Information to the MR.
  • the MR is required to transmit the Location Information to the RR.
  • both the MR and the RR are required to transmit their corresponding information to the TR.
  • Radio signals can interfere with one another unless multiple access techniques are employed.
  • Existing radio signals use any one of a number of multiple access techniques such as: Code Division Multiple Access (CDMA), Time Division Multiple Access (TDMA) and Frequency Division Multiple Access (FDMA).
  • CDMA is used in IS-95, the North American standard for CDMA CTs, as well as in GPS.
  • FDMA is used in AMPs, the North American standard for analogue CTs.
  • TDMA is used in IS- 136, the North American standard for TDMA CTs, and in GSM, the European standard for TDMA CTs.
  • GPS signals as an example of a satellite-based radio signal and IS-95 signals as an example of a land- based radio signal.
  • BST Base Station Transmitter
  • E91 1 91 1
  • RMS Root Mean Square
  • a handset-based Wireless Location System whereby an MR can be located using signals received by both the MR and an RR. More specifically, it is an object of the invention to estimate the static and kinematic positional information of a MR which receives signals transmitted from a number of BSTs.
  • a system that corrects for the sources of error that affect the different location technologies for a handset-based WLS in a unique and novel way.
  • a system that solves for the position of the MR using Least Squares where misclosures and standardized residuals are used to flag the observations that might have a blunder.
  • a system that uses Chaffee's method or Location On the Conic Axis (LOCA) to detect solution bifurcation.
  • LOCA Location On the Conic Axis
  • a system that uses LOCA and/or Plane Intersection to provide an initial position for Least Squares.
  • a system that uses individually, or in combination, TDOA observations, hybrid TDOA observations and AOA observations to locate the MR.
  • a system that resolves the AOA ambiguity at a cellular BS by using the TDOA information from two distant BSs or the AOA from several sectors at the same cellular BS.
  • a system that uses the range information from two BSTs or the AOA from several sectors at the same cellular base station to resolve the AOA ambiguity at such a BST.
  • WLS Wireless Location System
  • an MR can be located using signals received by both the MR and an RR. More specifically, it is an object of the invention to estimate the static and kinematic positional information of a MR which receives signals transmitted from a number of BSTs.
  • Yet another aspect of the invention is to estimate the Doppler shift ( ⁇ f,) in each received signal at the MR in order to estimate the speed and Direction of Travel (DOT) of the MR using Frequency Difference Of Arrival (FDOA). This applies to modifications: I and II.
  • Yet another aspect of the invention is to minimize the use of RRs in order to reduce the cost and complexity of a handset-based WLS. It is an object of the invention to produce a handset-based Wireless Location System (WLS) where an MR can be located without the need for an RR. This applies to all three modifications: 0, 1 and II.
  • a handset-based Wireless Location System whereby an MR can be located using signals received by both the MR and an RR. More specifically, it is an object of the invention to estimate the static and kinematic positional information of a MR which receives signals transmitted from a number of BSTs.
  • Another contribution of this invention is the way a CT is requested to act as an MR in a handset-based WLS where a transmitter in a cellular BS is considered to be a BST.
  • a CT to act as an MR, a certain level of CT modification might be required. Since there is a large number of existing CTs worldwide, it is important to keep in mind that some levels of modification to the CT might preclude existing CTs (which is not a desirable outcome). We attempt to reduce the amount of CT modification for a desired WLS performance level. Three levels of CT modification have been identified in this patent: modification 0 which requires no modification to the CT; modification I which requires only software modification to the existing CT; and modification II which requires both software and hardware modifications to the existing CT.
  • Yet another aspect of the invention is the use of the SR and Inverse SR algorithms over the correlation function (that is obtained by correlating the received signals at both the MR and RR with a clean version of themselves) in order to: resolve the correlation peaks beyond the Rayleigh resolution and to avoid using a threshold in selecting the correlation peak.
  • Yet another aspect of the invention is a system for locating an MR by estimating the Phase Of Arrival (POA) of existing or generated tones contained in the received signals.
  • the generated tones can be obtained by performing a nonlinear operation over the received signals at the MR.
  • the phases of the existing or generated tones can be extracted using SR algorithms in order to reduce the effect of multipath by selecting the tone corresponding to the first arrival in time. This applies to modifications: I and II.
  • Yet another aspect of the invention is use of the estimated POA to estimate either the TOA of the received signals or their AOA using Phase Difference Of Arrival (PDOA). This applies to modifications: I and II.
  • PDOA Phase Difference Of Arrival
  • Yet another aspect of the invention is provision of a system which estimates the Doppler shift ( ⁇ f) in each received signal at the MR in order to estimate the speed and Direction of Travel (DOT) of the MR using Frequency Difference Of Arrival (FDOA). This applies to modifications: I and II.
  • Yet another aspect of the invention is a system that minimizes the use of RRs in order to reduce the cost and complexity of a handset-based WLS.
  • WLS Handset-based Wireless Location System
  • a system that resolves the AOA ambiguity at a cellular BS by using the TDOA information from two distant BSs or the AOA from several sectors at the same cellular BS.
  • a system which uses the range information from two BSTs or the AOA from several sectors at the same cellular base station to resolve the AOA ambiguity at such a BST.
  • WLS Wireless Location System
  • an MR can be located using signals received by both the MR and an RR. More specifically, it is an object of the invention to estimate the static and kinematic positional information of a MR which receives signals transmitted from a number of BSTs.
  • a system that corrects for the sources of error that affect the different location technologies for a handset-based WLS in a unique and novel way.
  • a system that solves for the position of the MR using Least Squares where misclosures and standardized residuals are used to flag the observations that might have a blunder.
  • LOCA Location On the Conic Axis
  • a system that uses LOCA and/or Plane Intersection to provide an initial position for Least Squares.
  • a system that uses individually, or in combination, TDOA observations, hybrid TDOA observations and AOA observations to locate the MR.
  • a system that resolves the AOA ambiguity at a cellular BS by using the TDOA information from two distant BSs or the AOA from several sectors at the same cellular BS.
  • a system that uses the range information from two BSTs or the AOA from several sectors at the same cellular base station to resolve the AOA ambiguity at such a BST.
  • MR Mobile Receiver
  • MR Mobile Receiver
  • CTR Cellular Telephone
  • PCS Personal Communication Systems
  • PDA Personal Digital Assistant
  • RR Reference Receiver
  • RR is a device of known location with respect to the location of BSTs that are transmitting radio signals intended for reception by both the MR and the RR.
  • An RR can be portable or fixed.
  • the RR estimates the TOA, POA, SOA and/or FOA of signals transmitted by the BSTs. It either transmits the Reference Information to the MR or to the TR, or it receives the Location Information from the MR. In some unique situations when the BSTs are referenced to a common signal, the RR is not required.
  • Reference Information is the estimated TOA, POA, SOA and/or FOA of the signals transmitted by a number of BSTs at the RR as well as the location of such BSTs.
  • Location Information is the estimated TOA, POA, SOA and/or FOA of the signals transmitted by a number of BSTs at the MR.
  • TR Third Receiver
  • MSC Mobile Switching Center
  • BST Base Station Transmitter
  • a BST can be portable or fixed. It can consist of a transmitter alone or both a receiver and a transmitter. It can be a regular cellular Base Station (BS), a regular satellite transceiver, a GPS satellite, a Loran-C transmitter, a dedicated transmitter, a transponder, or any other type of transmitter/transceiver combination.
  • Host is a processor that processes the Location Information together with the Reference Information to estimate the location of the MR.
  • the function of the Host can be played by the MR, the RR or a TR.
  • the host can also be responsible for location services such as fleet management, location sensitive billing, etc.
  • Network-based WLS Form/Downlink-based WLS: locates a Mobile Receiver (MR) using a plurality of BSTs that transmit signals received by both the MR and the RR.
  • Handset-based WLS locates a Mobile Transmitter (MT) using a plurality of receivers that receive the signal transmitted from the MT.
  • Cell is a geographical area serviced by a cellular Base Station (BS).
  • BS Base Station
  • Sectorized cell is a cell that is made of several non-overlapping sectors. Each sector can be considered as an independent cell to be serviced by an independent cellular BS. However, sectors in a cell are usually serviced by the same cellular BS in order to minimize cost and complexity. We refer to such a BS as a sectorized BS. Pilot Channel: is the channel to be used by each sector to broadcast a (locally unique) pilot signal.
  • Pilot Signal is defined as a sequence known to the MR.
  • this sequence is a code repeated exactly 75 times every 2 seconds, and is known as the "pilot PN sequence” or the “short code”. Although all sectors broadcast the same pilot signal, the transmissions from each base station have different synchronization or phases, which distinguish the pilot signals. Phase increments are specified as multiples of 64 chips, where 1 chip is 1/1.2288 microseconds. The offset in units of 64 chips, is referred to as the pilot PN offset. Pilot increments in a network are typically 3 or 4 PN offsets. A PN offset of 4 corresponds to 256 chips, or 208.3 microseconds.
  • a SR algorithm transforms a frequency domain signal, S(f), to a time domain signal, s(t), with a better resolution than the resolution offered by the inverse Fourier transform, i.e. s(t) has a better resolution than the inverse Fourier transform F"' ⁇ S(f) ⁇ , of S(f)).
  • SR algorithms are well known in the literature and include: MUSIC/Root MUSIC, ESPRIT, Auto Regressive Moving Average, Minimum Variance, MUSIC using Higher Order Statistics, ESPRIT using Higher Order Statistics, Auto Regressive Moving Average using Higher Order Statistics, or Minimum Variance using Higher Order Statistics.
  • Frequency-domain Rayleigh Resolution is the frequency-domain resolution offered by the Fourier transform.
  • Time-domain Rayleigh Resolution is the time-domain resolution offered by the inverse Fourier transform.
  • Inverse SR algorithm is an algorithm which processes a time domain signal in order to improve its time-domain resolution over the conventional time-domain Rayleigh resolution.
  • a preferred embodiment of the inverse SR algorithm comprises a conventional time domain correlator, a time domain window, a Fourier Transform, a frequency domain window, a frequency domain equalizer, or a processor which performs a SR algorithm in order to resolve the TOAs in the received signal at a given MR.
  • an inverse SR algorithm can process a frequency domain signal in order to improve its frequency-domain resolution over the conventional frequency-domain Rayleigh resolution.
  • a preferred embodiment of the inverse SR algorithm comprises a frequency domain correlator, a frequency domain window, an inverse Fourier Transform, a time domain window, a time domain equalizer, or a processor which performs a SR algorithm in order to resolve the FOAs in the received signal at a given MR.
  • Sector Antennas are cellular antennas that exist in some BS that service more than one sector. Each sector has a dedicated antenna (or set of antennas if diversity is applied).
  • Fourier Transform-Based Filters are filters that: Fourier Transform the time domain signal, then window the transformed signal over a given band, and Inverse Fourier Transform the windowed signal.
  • the CT measures the quality of received pilot signals by computing the ratio between the received energy of a pilot (E c ) to the total received power by the mobile (I 0 ). This measurement is referred to as the E c /I 0 of the pilot. This measurement is closely related to the signal strength of the pilot.
  • An IS-95 CT uses the E c /I 0 measurement of a pilot to estimate the suitability of the communications link between itself (the CT) and that particular sector in a BS. For example, if a CT is in a mode in which it may only communicate with one sector at a time, it may choose to communicate with the sector which has the highest E c /I 0 measured by the CT.
  • An IS-95 CT may include a searcher element and multiple demodulator fingers as shown in Figure 12.
  • a searcher is used to search for pilot signals by measuring the E c /I 0 of a pilot signal at different delays.
  • a demodulator finger tracks a pilot signal, and may demodulate an information signal associated with said pilot signal.
  • a base station may also transmit one or more information signals associated with the pilot. These information-bearing signals are transmitted at the same radio frequency as their associated pilot signal and are synchronized with the pilot signal. Thus, an information signal will undergo propagation effects (attenuation, phase shift, time delays, multipath, frequency shift) nearly identical to those of its associated pilot signal.
  • a mobile then tracks the propagation effects of a pilot to aid in demodulating information from an associated signal, or code channel. As well, the mobile can estimate the reliability of information demodulated from a code channel from the E c /I 0 quality measurement of the associated pilot signal.
  • a CT searches for pilot signals using one or more searchers, and then based on the searcher results, assigns demodulator fingers to demodulate information from code channels associated with suitable pilots. This process is ongoing so that the finger assignment adapts to the changing conditions typical of a CT cellular environment.
  • the searchers and demodulator fingers in an IS-95 CT are usually assigned according to algorithms which optimize the communications performance of the phone, and the network as a whole. These assignment algorithms will limit the performance of estimation of the CT location, as they will often ignore pilot signals useful for the location estimation. We prefer to force the CT to enter into a state optimized for location.
  • a method of estimating the location of a Mobile Receiver comprising the steps of: receiving a plurality of signals from a plurality of Base Station Transmitters (BSTs) by the MR; and performing a nonlinear operation on the plurality of received signals at the MR in order to generate at least one tone per received signal. Additional refinements of this aspect of the invention are set forth in claims 2-5 attached hereto, which are hereby incorporated in this summary by reference.
  • a method of estimating the location of a Mobile Receiver comprising the steps of: receiving a plurality of signals from a plurality of Base Station Transmitters (BSTs) by the MR; and estimating the location of the MR by estimating the phases of existing tones in each received signal.
  • BSTs Base Station Transmitters
  • a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs where one or more antennas is used at the MR and in which location of the MR is achieved using a technique selected from the group consisting of:: AOA positioning (multi-angulation), TOA positioning (circular multilateration),
  • Range positioning (circular multilateration), TDOA positioning (hyperbolic multilateration), TDOA positioning (circular multilateration with TDOAs), AOA/TDOA positioning (multi-angulation/hyperbolic multilateration), AOA/TDOA positioning (multi-angulation/circular multilateration with TDOAs),
  • AOA/Range positioning multi-angulation/circular multilateration
  • AOA/Range positioning multi-angulation/circular multilateration
  • combining the signals from all antennas for location purposes Additional refinements of this aspect of the invention are set forth in claims 19, 20 and 22 attached hereto, which are hereby incorporated in this summary by reference
  • a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs using equation (29) and equation (30) (defined below).
  • a method of estimating the location of a MR wherein location of the MR is achieved based on a solution selected from the group consisting of: solution I, solution II, solution III, and solution IV (as described below).
  • a system to estimate the location of a Mobile Receiver comprising of: a receiver at the MR for receiving a plurality of signals from a plurality of Base Station Transmitters (BSTs); and a computing means connected to the receiver by a communication linkfor performing a nonlinear operation on the plurality of received signals in order to generate at least one tone per received signal.
  • a system to estimate the location of a Mobile Receiver comprising of: a receiver at the MR for receiving a plurality of signals from a plurality of Base Station Transmitters (BSTs); and a computing means connected by a communication link to the receiver for estimating the location the MR by estimating the phases of existing tones in each received signal.
  • BSTs Base Station Transmitters
  • Range positioning (circular multilateration), TDOA positioning (hyperbolic multilateration), TDOA positioning (circular multilateration with TDOAs), TOA/TDOA positioning (multi-angulation/hyperbolic multilateration),
  • AOA/TDOA positioning multi-angulation/circular multilateration with TDOAs
  • AOA/Range positioning multi-angulation/circular multilateration
  • AOA/Range positioning multi-angulation/circular multilateration
  • combining the signals from all antennas for location purposes Additional refinement of this aspect of the invention is set forth in claim 29 and 31 attached hereto, which are hereby incorporated in this summary by reference.
  • a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, based on any one of the following solutions: solution I, solution II, solution III, or solution IV.
  • a system to estimate the location of a Mobile Receiver comprising a Receiver that receives a plurality of signals transmitted from a plurality of BSTs; wherein the Receiver includes an IF-sampling Receiver comprising: a number of IF stages to convert a received RF signal r,(t) to an analog IF signal of IF frequency f i ; a pre-aliasing filter; an Analog-to-Digital (A/D) converter that converts the IF analog signal to an IF digital signal; and a processor for processing the IF digital signal.
  • MR Mobile Receiver
  • a system to estimate the location of a MR using a Reference Receiver having a LO comprising: an estimator for estimating the carrier offset (due to Doppler, ⁇ f backyard and due to LO offset, ⁇ f,) of a plurality of signals transmitted by a plurality of BSTs and received by the RR. Additional refinements of this aspect of the invention are set forth in claims 45- 54 and 60 attached hereto, which are hereby incorporated in this summary by reference.
  • a system to estimate the location of a MR having a LO comprising: an estimator for estimating the carrier offset (due to Doppler, ⁇ f,, and due to LO offset, ⁇ f,) of a plurality of signals transmitted by a plurality of BSTs and received by the MR.
  • a method of estimating the location of a MR comprising the steps of: receiving a plurality of signals transmitted from a plurality of BSTs at a Receiver; converting the received RF signal r,(t) to an IF signal of IF frequency fi in an
  • IF-sampling Receiver filtering the IF signal with a pre-aliasing filter; converting the filtered IF analog signal to an IF digital signal in an Analog-to- Digital converter; and processing the IF digital signal in a signal processor.
  • a method of estimating the location of a MR comprising the steps of: estimating the carrier offset (due to Doppler, Df,, and due to LO offset, Df,) of a plurality of signals transmitted by a plurality of BSTs and received by a RR.
  • a method of estimating the location of a MR comprising the steps of: estimating the carrier offset (due to Doppler, ⁇ f backyard and due to LO offset, ⁇ f,) of a plurality of signals transmitted by a plurality of BSTs and received by the MR.
  • a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs comprising the steps of: characterizing the relative group delay, ⁇ t, - ⁇ t j5 at a RR between a signal received from the i th BST and a signal received from the j th BST prior to operation; removing the characterized relative group delay, ⁇ t, - ⁇ t j , from the estimated Time Difference Of Arrival, ⁇ , - ⁇ ⁇ , at the MR between a signal received from the i th th BST and a signal received from the j BST to yield a corrected Time Difference of Arrival; and finding the location of the MR using the corrected Time Difference of Arrival.
  • a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs comprising: using Least Squares to solve for the positional information of the MR and using TDOA observations from distant BSTs to resolve for ambiguities in an AOA solution; using AOA observations from different sectors at the same cellular BS to resolve for ambiguities in an AOA solution; and using range observations from different BSTs to resolve for ambiguities in an AOA solution.
  • a system for estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs the system including a host, and the system comprising; the MR being adapted to receive some or all the following information from a remote source and to apply said information to the execution of a location state: identification of base station pilot signals for the mobile to search for pilot detection threshold; duration of time for searcher to accumulate pilot energy at each delay offset; search window time length or duration; spacing of successive sampling points within the search window; frequency at which to conduct the search for pilot signals; cellular system which the mobile should acquire; the number of times the mobile should repeat all or part of the location procedure; correlation interval; locations of base stations; and the MR being adapted to transmit one or more of the following information to the Host: times of arrival of pilot signals; identification of the pilot signals reported by said device; received signal qualities of pilot signals
  • a system for estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs comprising: first computing means for characterizing the relative group delay, ⁇ t, - ⁇ t j5 at a RR between a signal received from the i th BST and a signal received from the j th BST prior to operation; and second computing means for removing the characterized relative group delay, ⁇ t, - ⁇ t j , from the the estimated Time Difference Of Arrival, ⁇ , - ⁇ , at the MR between a signal received from the i th BST and a signal received from the j th BST.
  • a system for estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs where: Least Squares is used to solve for the positional information of the MR;
  • TDOA observations from distant BSTs are used to resolve for ambiguities in an AOA solution
  • AOA observations from different sectors at the same cellular BS are used to resolve for ambiguities in an AOA solution; and AOA observations from different sectors at the same cellular BS are used to resolve for ambiguities in an AOA solution.
  • a method to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs the MR receiving some or all the following information from a remote source and applies said information to the execution of a location state: identification of base station pilot signals for the mobile to search for pilot detection threshold; duration of time for searcher to accumulate pilot energy at each delay offset; search window time length or duration; spacing of successive sampling points within the search window; frequency at which to conduct the search for pilot signals; cellular system which the mobile should acquire; the number of times the mobile should repeat all or part of the location procedure; correlation interval; locations of base stations; and wherein some or all of the following information is transmitted from the MR to the Host: times of arrival of pilot signals; identification of the pilot signals reported by said device; received signal qualities of pilot signals; pilot search results; carrier frequencies of reported pilot signals.
  • a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs where misclosures and/or standardized residuals are used to flag observations that might contain a blunder.
  • MR can be located with as few as two AOA independent observations and one TDOA observation.
  • a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs wherein a combination of AOA and TOA observations are used such that the MR can be located with as little as one AOA observation and one TOA observation.
  • a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs where a combination of TOA observations are used such that the MR can be located with as little as two independent TOA observations.
  • a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs where a combination of AOA and TDOA observations are used such that the MR can be located with as few as two AOA independent observations and one TDOA observation.
  • a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs which quantifies the effect of geometry on position estimation, with both AOA and TDOA observations, using the design matrix in equation (46).
  • a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs which quantifies the effect of geometry on speed and direction of travel estimation using the design matrix in equation (51).
  • Figure 1 illustrates the broadcast of a signal, s,(t) (103), by the i th BST (102) where ⁇ 0, , is the Time of Transmission (TOT) of the signal.
  • a MR (101) receives the signal, together with other signals transmitted from a number of BSTs.
  • the signal, s,(t) (103) can be a pilot signal, a sync signal, a paging signal or a Forward Traffic signal.
  • s,(t) is a pilot signal, it consists of a single code, repeated exactly 75 times every 2 seconds, with a period of 32,768 chips.
  • the pilot code is known as the "pilot PN sequence" or the "short code”.
  • the pilot transmitted by each sector is time-offset by an integral multiple of 64 chips from any other pilot signal which could forseeably be simultaneously received by a MR. As the code has a period of 32,768 chips, and pilots are spaced by 64 chips, there are 512 unique pilot phases possible.
  • the signal, s,(t) (103), can be either a Coarse Acquisition (CA) signal or a Precision signal.
  • CA Coarse Acquisition
  • s,(t) is a CA signal, it consists of a unique Gold code per satellite (i.e. up to 24 distinct codes), repeated exactly 1000 times every 1 second, with a period of 1023 chips, with BPSK overhead messages at 50 bits/second.
  • the codes are known as the "CA-codes”.
  • the signal, s,(t) (103) can be a Forward Control Channel (FOCC) signal or a Forward Voice Channel (FVC) signal.
  • FVC Forward Voice Channel
  • DCC Digital Colour Codes
  • the time ⁇ 0 ,i is the time of transmission of a given pilot PN sequence while ⁇ i is the time of reception of such a PN sequence.
  • the time ⁇ 0 ,i is the time of transmission of a given DCC sequence while ⁇ j is the time of reception of such a sequence.
  • Figure 3 illustrates Angle of Arrival for several wavefronts from BSTs at a MR, and assumes that the MR is far from BST, and BST, with respect to the baseline between the two BSTs. Such an assumption implies that the transmitted wavefronts are planar, i.e. ⁇ y, where: • ⁇ s is the Angle Of Arrival (AOA), of signal s,(t) from the i th BST to the MR,
  • AOA Angle Of Arrival
  • ⁇ j is the AOA of signal s,(t) from the j th BST to the MR
  • BSTs (which we refer to as line 0 ,ij) in a clockwise manner from line 0 ,i,j to liney.
  • Figure 4 illustrates the Direction of Travel (DOT) ' ⁇ ' (401) relative to Northing (in a clockwise manner from Northing) and the speed v (402) of the MR (404) relative to the BSTs (403).
  • DOT Direction of Travel
  • v speed of Travel
  • ' ⁇ ' and v represent the velocity v of the MR with respect to the BSTs (403).
  • the BSTs are stationary which is not the case in a satellite- based system.
  • Figure 5 illustrates the two-dimensional (horizontal) Locus of Position (501) for TDOA2 ; l which is defined as
  • TDOA 2>1 ( ⁇ 2 - ⁇ 0j2 ) - ( ⁇ , - ⁇ 0 ⁇ l )
  • TDOA 2 1 ( ⁇ 2 - ⁇ 0 ⁇ 2 ) - ( ⁇ 1 - ⁇ 0 ⁇ l )
  • Figure 6 illustrates the two intersecting locii of Position for the two-dimensional coordinates (x,y) of the MR (603) based on TDOA2,l (601) and TDOA3 (602).
  • (x,y) it is possible to solve for (x,y) as the intersection between the two trajectories obtained using three Times of Arrival (TOAs) (after choosing the correct side of each hyperbola).
  • TOAs Times of Arrival
  • Figure 7a illustrates the Second Stage of the WLS where the signal, s,(t), is received by the MR (701).
  • the MR (701) processes the received signal, r,(t), to estimate its own coordinates (x,y), speed v and DOT ⁇ with the assistance of a Reference Receiver (RR) (702) of known location.
  • the function of the RR is to estimate the Reference Information regarding the BSTs and to transmit it to the MR (701).
  • Option I See Section I.
  • Figure 7b illustrates an alternative of the Second Stage of the Wireless location system where the signal, s,(t), is received by the MR (703).
  • the MR processes the received signal, r,(t), and transmits the Location information to the RR (704).
  • the RR estimates the coordinates (x,y) of the MR, its speed v and DOT ⁇ , based on the Location Information transmitted by the MR (703) and based on the Reference Information estimated by the RR (704).
  • Option II the Reference Information
  • Figure 7c illustrates yet another alternative of the Second Stage of the Wireless location system where the signal, s,(t), is received by the MR (705).
  • the MR processes the received signal, r,(t), and transmits the Location Information to a Third Receiver (TR) (706).
  • the RR (707) estimates the Reference Information regarding the BSTs and transmits it to the same Third Receiver (TR) (706).
  • the TR estimates the coordinates (x,y) of the MR, its speed v and DOT ⁇ , based on the Location Information transmitted by the MR (705) and based on the Reference Information transmitted by the RR (707).
  • Option III we will refer to such a system as Option III.
  • FIG. 8 illustrates the description of Design I for the MR.
  • the received RF signal r,(t) transmitted from the i th BST is initially filtered by an RF Band Pass Filter (BPF) (802), amplified by an RF amplifier (803), down-converted by one or more Intermediate Frequency (IF) stages (810).
  • BPF RF Band Pass Filter
  • IF Intermediate Frequency
  • Each IF stage consists of a mixer (804) with an LO (808) generated using an RF synthesizer (809), an IF BPF (805) and an IF amplifier (806).
  • the resulting analog IF signal (807) of IF frequency fi is filtered by a pre-aliasing filter (811), sampled by an A/D (812) and filtered again using digital filters (813) in order to generate a digital IF signal (814) of center frequency f 2 , where f 2 ⁇ fi .
  • An optional Direct Digital Converter (DDC) (815) is used to down-convert the digital IF signal (814) to Baseband.
  • DSP Digital Signal Processor
  • FOA Frequency Of Arrival
  • a function of the estimated FOA is fed- back to a Direct Digital Synthesizer (DDS) (819) which controls the RF synthesizer (809).
  • the DSP (817) also provides a reference signal (820) to the DDS (819).
  • a preferred source for the reference signal (820) is one that is common to all BSTs such as a GPS signal.
  • Figures 9a, 9b, 9c and 9d illustrate a flow diagram for the positioning strategy.
  • the positioning strategy consists of a number of algorithms (Least Squares, Location On the Conic Axis (LOCA), Plane Intersection) and numerous decisions. Given a set of observations, there are 4 possible outcomes: two LS position solutions (908), two LOCA position solutions (909), one LS position solution (923, 936), and no position solution (927).
  • LOCA Location On the Conic Axis
  • Figure 10 presents the preferred embodiment of an algorithm by which a MR can obtain TOA estimates of multiple IS-95 pilot signals.
  • a search is made for each pilot over multiple time lags. If a pilot is detected, then it is further observed to obtain a TOA estimate of refined accuracy.
  • Figure 11 shows correlation values which might be obtained when searching for a particular pilot over multiple time lags, or delay offsets.
  • Ec/Io is the ratio of average pilot chip energy to the total received power. A higher value of Ec/Io indicates greater likelihood of arrival of the pilot signal at a particular delay offset.
  • the delay offsets are expressed in IS-95 chips (approximately 0.814 ⁇ s).
  • the horizontal dotted line indicates a threshold (-15 dB in this case) for detection of the pilot arrival. An Ec/Io value greater than this threshold is considered to indicate an arrival of the pilot signal at the MR, at that particular time lag.
  • Figure 12 indicates a possible architecture for the embodiment of the MR as a CDMA CT, which may be used to implement the algorithm of Figure 10.
  • the Searcher searches for pilot signals at different time lags of the received signal of the MR.
  • the Demodulator Fingers are assigned to different pilot signals. When assigned to a pilot, a Demodulator Finger will slew its timing to track the pilot, and may demodulate data transmitted synchronously with that pilot.
  • This architecture is presented in U.S. Patent number 5,764,687 Mobile Demodulator Architecture for a Spread Spectrum Multiple Access Communications System.
  • the Low Pass (LP) equivalent signal, (t) transmitted from the i ⁇ BST can be modeled as
  • ⁇ (t) G j R e j( ⁇
  • r ⁇ (t) Re ⁇ e j2 ⁇ f ⁇ (t) ⁇ (6)
  • G represents the gain (real) due to the transmitting and receiving antennas: from the 1 th BST to the MR antenna (i.e. a function of the two antenna patterns),
  • n is a real number (usually between 2 and 4) that depends on the channel and .
  • R ⁇ is the range between the i ⁇ BST and the MR,
  • ⁇ j ⁇ -2 ⁇ ( + ⁇ + ⁇ fi ) (xj - ⁇ 0 ⁇ i + ⁇ t, ) ⁇ mod 2 ⁇ (7)
  • . ⁇ is the frequency offset due to the Doppler shift over the propagation channel (a function of speed, v, with respect to the BSTs and Direction Of Travel (DOT), ⁇ , with respect to Northing).
  • . X; TOA of a signal p;(t-x 0 f ) (a function of the range R,), . ⁇ t j is the overall group delay through the i m BST to the MR antenna (usually a function of frequency), . w ; (t) represents the Additive White Gaussian Noise (AWGN, complex) due to thermal noise at the MR (a function of temperature and bandwidth), .
  • MP j (t) represents all multipath components (complex) at MR antenna (a function of the environment and of the elevation of the antennas), and 00/00224
  • L (t) represents low-pass equivalent interference (both in-band and out-of- band) at the MR.
  • Total RF Reception In IS-95, a number of pilot signals, sync signals, paging signals, and forward traffic signals, are transmitted over the same frequency band in the forward link. In GPS, a number of signals are transmitted from the satellites over the same frequency band. In other words, in a CDMA system I(t) represents all CDMA signals excluding ⁇ (t).
  • ⁇ f 0 is the frequency error between the Local Oscillators (LOs) in the MR receiver and the carrier frequency f,.
  • ⁇ o is the carrier phase in the Local Oscillators (LOs) of the MR.
  • an Analog-to-Digital (A-to-D) Converter converts the baseband signal, ?(t), to a digital signal, r(m) , which is equivalent to performing the following operation:
  • • m is an integer; and • q(m) is the quantization noise which depends on the number of bits in the A-to-D.
  • _Rp the first source, Rp, can be used to estimate R;
  • Rp is unreliable due to the nature of the radio channel where RF shadowing and flat fading can deviate the value of "n" in "Rp” from 2 (for Line Of Sight (LOS)) to 4 (for an urban environment, as shown in Hata, M., "Empirical Formula for Radio Propagation Loss in Land Mobile Radio Services," IEEE Transactions on Vehicular Technology, Vol. VT-29, No. 3, August 1980), or even 6 in heavy urban environments.
  • LOS Line Of Sight
  • Rp can only offer an accuracy of a few kilometers for Rj.
  • G has to be known to the i m BST and "n" can be chosen to be 2.
  • the POA, ⁇ ; is related to the range, R,, between the i A BST, and the MR through equations (10) and (11):
  • RTD Round Trip Delay
  • ⁇ is the time of transmission from the MR to the BST and 2 ; is the TOA of the response of the BST to the MR.
  • equation (12) assumes that the BST is able to transmit a response back to the MR. This assumption makes the handset-based WLS depend on the functionality of the BSTs. This implies that the handset- based WLS is network dependent. Nonetheless, we will carry on using such a concept in order to be able to aid a handset-based solution with a network- based solution.
  • the amount of ambiguity is a function of the baseline between the two transmitting antennas, i.e. between the i m BST and the j ⁇ l BST with respect to the wavelength corresponding to the frequency: f c + ⁇ f ; + ⁇ f, in
  • n is the integer ambiguity
  • the CT frequency, ⁇ . + ⁇ can deviate by up to 0.05part per million (i.e. by up to 95Hz for a 1900MHz carrier frequency).
  • the carrier frequency, f c can take values around 800MHz (with a wavelength ⁇ of around 37.5cm) or values around 1.9GHz (with a wavelength ⁇ of around 15.8cm). This implies that for a velocity of lOOkm/hr ⁇ f corresponds to:
  • Solution I When the baseline, d, , between the antenna at the i" 1 BST and the antenna at the j" 1 BST is of the same order of magnitude as the wavelength, ⁇ , it is possible to generate a solution using either TDOA or PDOA.
  • the TDOA solution between the two antennas is a hyperbola while the PDOA solution between the two antennas is a line. Both solutions approximate the exact solution and coincide asymptotically, i.e. the PDOA line and the TDOA hyperbola coincide at infinity.
  • the baseline, cL j can be of the same order of magnitude as ⁇ only if the two antennas belong to different sectors of the same cellular base station.
  • the transmitted wavefront is planar, i.e. where:
  • Y J is the Angle Of Arrival (AOA) of r £ (t) from the antenna of the i ⁇ BST to the MR
  • Y j is the AOA of r-(t) from the antenna of the j m BST to the MR
  • the Cramer-Rao Lower Bound on the variance, van ⁇ ), of the estimated AOA, ⁇ , is equal to 0224
  • Equation (18) is based on the transmission of a narrowband tone from each element of the array at the same phase.
  • the transmitted signal Sj(t) is either a wideband signal (e.g. CDMA) or s ; (t) has a unique frequency f ⁇ (e.g. FDMA), (18) no longer applies directly, but can be made to apply after some mathematical adjustments to the received signals using the Reference Information estimated by the RR.
  • vanR j is the variance of the estimate, R ⁇ , of the range R;
  • • fj n is the frequency of the tone corresponding to ⁇ (i.e.
  • A is me amplitude of the tone at the MR from the 1 th BST
  • equation (19) is based on the transmission of a narrowband tone from each element of the array at the same phase.
  • the transmitted signal s ; (t) is either a wideband signal (e.g. CDMA) or s ⁇ (t) has a unique frequency ⁇ (e.g. FDMA)
  • (19) no longer applies directly, but can be made to apply after some mathematical adjustments to the received signals using the Reference Information estimated by the RR.
  • Phase offsets The phase of the tone f ⁇ is shifted by ⁇ ; .
  • ⁇ is a function of the frequency fj j , it has to be estimated otherwise the phase ⁇ j is distorted.
  • the propagation delay x is affected by ⁇ t j which represents the system delay through the antenna, cables, filters, amplifiers, etc.
  • ⁇ t j is a function of the frequency J ⁇ ,, it has to be estimated otherwise the propagation delay can vary significantly from one received signal to another.
  • multipath, MP,(f) In TDOA, the multipath MP j (t) is equivalent to extra delay over the propagation channel and has to be either estimated and removed, or mitigated. In the case where the phase estimate is used to solve for the AOA of the received signal r ; (t), the effect of multipath is to shift the AOA of r £ (t) by some amount depending on the AOA of MPj(t) and its magnitude.
  • Solution TJ The most common method for estimating x ; is to cross-correlate Pi(t - t j + ⁇ 0 j - ⁇ t j ) with P j (t), i.e. to cross-correlate r(t) with p £ (t).
  • P j (t) i.e. to cross-correlate r(t) with p £ (t).
  • peaks result from such a correlation particularly when receiving pilot signals in CDMA where all transmitted pilots are identical except for a relative delay between pilots.
  • the peak corresponding to x ; in p,(t - X j + x 0 j - ⁇ t j ) is located within a reasonable propagation delay, Xj - x o i , with respect to a clock derived from the signal, and can be distinguished from other pilots.
  • a possible method is to estimate the RTD between the MR transmission and the reception at the MR of the BST response to the MR transmission as shown in equation (12).
  • the resulting correlation peak corresponding to ⁇ ⁇ in P j (t - X j + x 0j - ⁇ t j ) is located within a reasonable propagation delay, ⁇ ⁇ - x oJ , with respect to a clock derived from the signal.
  • the TDOA, x, - ⁇ s can be used to estimate the relative delay, (x, - x 0 ,) - (x J - x OJ ) since the relative TOT, x 01 - x oj , between the i m pilot and the j" 1 pilot is known at the MR.
  • the received wavefront is planar, i.e. Y J ⁇ YJ ⁇ Y J J, where
  • Yj is the Angle Of Arrival (AOA) of ⁇ (t) from the antenna of the i ⁇ BST to the MR
  • Y j is the AOA of rj(t) from the antenna of the j* BST to the MR
  • Equation (21) relates the angle ⁇ - to the Time Difference of Arrival (TDOA), X; - X j , as follows
  • • dj j is the distance between the antenna at the i m BST and the antenna at thej ⁇ BSTj and
  • equation (22) is the fact that there are no ambiguities in equation (22) except for the ⁇ in the solution of Y j , while equation (23) can have a number of ambiguities depending on the value of the integer k.
  • Tins is equivalent to having a standard deviation for ⁇ of 7.07 x lO -3 radians, which corresponds to a range error of 7.11 meters for every kilometer range between the MR and the antennas.
  • p,(t - x, + x 0>1 - ⁇ t,) in (20) is suitable for TDOA positioning since in this case the baseline between antennas does not have to be small.
  • the lower-bound on the variance, va ⁇ R.)l c 2 vai j) , of the range R, is 2553meters 2 or equivalently ' 'p. ) ' P.
  • R j has a standard deviation of 50.53meters.
  • t in equation (25) to be lOdB for IS-95, the range Rj is
  • R s has a standard deviation of 12.28meters.
  • Solution TV When x ; is estimated indirectiy based on the phases of tones in Pi(t - Xj + x 0 j - ⁇ tj), an integer ambiguity can still exist in the solution unless the wavelength of the selected tones is large with respect to the range between the MR and the antenna at the i ⁇ BST.
  • the selected tones are discussed below where a solution for the integer ambiguity is presented.
  • BW is the RMS bandwidth of p,(t), .
  • (t) 2E, / N 0 ,
  • N 0 is the noise Power Spectral Density
  • the radio frequency (RF) channels have a bandwidth of approximately 1.25 MHz which is comparable to that of systems designed primarily for location such as Global Positioning Systems (GPS) with a BW of lMHz over Standard Positioning Services (SPS) channels and ISM-based location systems with a typical bandwidth of 10MHz.
  • GPS Global Positioning Systems
  • SPS Standard Positioning Services
  • ISM-based location systems with a typical bandwidth of 10MHz.
  • the location system may initially use a conventional sliding correlator to obtain a set of pseudo-ranges (one pseudo-range per satellite) (see e.g.
  • the Observation interval, T ; is directly related to the energy E ; in the received signal (t) at the MR since
  • BWj is the bandwidth of the filter that is processing (t).
  • the propagation delay ⁇ ⁇ is affected by ⁇ t j which represents the system delay through antenna, cables, filters, amplifiers, etc.
  • ⁇ tj represents the system delay through antenna, cables, filters, amplifiers, etc.
  • I j (f) can have a drastic effect on the accuracy of the estimated phase and may saturate the RF front end of the receiver. Its effects can be mitigated with both analog and digital hardware, as well as adequate software, whether it is in-band or out-of-band.
  • the multipath MPj(f) is equivalent to an extra delay over the propagation channel and has to be either estimated and removed, or mitigated.
  • the estimated time delay, ⁇ (> is used to solve for the TDOA between the received signal r,(t) and the received signal r,(t)
  • the difference in multipath, MPj(t) - MP j (f) is the distorting factor which has to be mitigated.
  • the cellular frequency band is between 824 and 894 MHz and the propagation characteristics at these UHF frequencies will have a significant impact on positioning by multi-lateration as shown in Parsons D., "The Mobile Radio Propagation Channel, " John Wiley & Sons, New York, 1992. That the ranges measured correspond to Line Of Sight (LOS) distances is a major assumption made when estimating position by multilateration. Although the dominant transmission mode in this band is LOS, reflections from natural and man-made objects as well as diffraction around said objects are also possibilities. Multipath and diffraction allow the cellular signal to propagate in heavily built up areas as well as indoors.
  • the measured ranges are longer than the true LOS distance which introduces error into the multilateration process.
  • the propagation distance at UHF is relatively short. This allows frequency reuse in the cellular system but limits the number of observables in the multilateration process. For instance, in a dense urban environment with a delay spread of 3 microseconds (as shown in Hata, M., "Empirical Formula for Radio Propagation Loss in Land Mobile Radio Services," IEEE Transactions on Vehicular Technology, Vol. VT-29, No. 3, August 1980) multipath causes the location accuracy to degrade to more than 1400m RMS.
  • the reason for this is that the correlation function from which the multipath may be estimated has a resolution which is limited to that of the Fourier transform which implies that any multipath within such a resolution is unresolvable using traditional methods. Further processing using an inverse SR algorithm often yields a result with higher multipath resolution as shown by Dumont, L.R., et al., " Super- resolution of Multipath Channels in a Spread Spectrum Location System," IEE Electronic Letters, Vol. 30, No. 19, pp. 1583-1584, September 15, 1994 and as shown by Fattouche et al., U.S. Patent No. 5,570,305 issued Oct., 1996, and as shown by Ziskind, I.
  • Tj(t) When the received signal, Tj(t), is periodic (e.g. the pilot signal in IS-95, the CA-code in GPS), it can be approximated using Fourier series, or equivalently, a number of tones can be used to approximate the pilot sequence:
  • the received signal, r ; (t), When the received signal, r ; (t), is cyclo-stationary, it generates a tone when a nonlinear operation is performed over it.
  • the phase of such a tone carries location information regarding the MR.
  • Most digital signals are cyclo-stationary and the generated tone is usually at the symbol rate (except in CDMA where the generated tone is the chipping rate).
  • the BSTs are AMPS TDMA/GSM-based
  • taking advantage of the cyclo-stationarity nature of the received signal is particularly important since in this case, it is not important to have a clean replica of the received signal.
  • a tone is generated of known (or estimated) frequency. This implies that the duration of observation can be extended as long as there is a signal to be received, i.e. as long as a call is in progress regarding a CT. This leads to an improved SNR for the received signal as explained in the CRLBs shown above.
  • a Treatment is the clockwise angle from Northing to the line formed by (x ⁇ y n ) and (x,y), and
  • Solution V we refer to such a solution as Solution V. Note that in Solution V the position of the MR (x,y) must be known (or estimated) prior to estimating the speed and velocity of the MR, in order to be able to know A ⁇ ⁇ * and Alic in equations (33).
  • is the variance of the estimate, ⁇ f, , of the frequency ⁇ f trip iTon •
  • SNR j ⁇ A? /(2 ⁇ 2 ) is the SNR of the tone at the MR from the antenna ofthe i ⁇ BST,
  • A is e amplitude of the tone at the MR from the antenna of the i" 1 BST
  • N j is the number of samples (which is directly related to the observation interval T ⁇ ).
  • I j (t) plays the same role in distorting the estimate of the frequency ⁇ as in distorting the estimate of the phase ⁇ , except that its effect can be reduced more significantly.
  • Ranges may be calculated by subtracting the known time of transmission x 0 of signal s(f) from the measured TOA, X j , of signal r ; (t). There are three possible solutions: 1. When the Time of Transmission, x 0 , is known, we refer to the positioning technique as Range positioning.
  • one way of deterrnining the time of transmission, x 0 is to use the RTD between a MR and a BST with respect to the round trip delay from the time the MR transmits a signal to the BST to the time it receives the response of the BST to the MR transmission.
  • a minimum of two independent equations is required to solve for the two unknowns x and y. If two ranges are estimated from TOA measurements of two different BST's, then there are two solutions for MR location at the intersections of two circles centered at the two BST's. This ambiguity can be resolved by including a third range measurement from another BST, or by (approximate) angle estimation from one or more sites by examining signal strength or phase difference of arrival.
  • X j is the Time Of Arrival of signal r ; (t) at the MR from the antenna of the i ⁇ BST, x 0 is the time of transmission of signal s(t) from the i ⁇
  • DOP Dilution Of Precision which is a measure of geometry.
  • HDOP Horizontal DOP
  • EDOP is the East DOP and is defined as the square root of the element in the 1 st row and l sl column of C ⁇ OA , 0224
  • 50 NDOP is the North DOP and is defined as the square root of the element in the 2 nd row and 2 nd column of -TOA'
  • C T0A is an unsealed matrix defined as
  • C is the unsealed measurement co variance matrix (the identity matrix of appropriate dimension)
  • A is the design matrix for the model of equation (35) and is defined as x-x, y-yi -c d. d, x-x 2 y- y 2 -c
  • dj is the best derived distance (range) between the MR and the antenna at the i ⁇ BST.
  • TOA observations thought to be more reliable may be weighted more heavily than those deemed less reliable. This is accomplished through the observation covariance matrix.
  • the inverse of C, is the weight matrix.
  • RSSI from each BST is one method of assigning weights to the TOAs.
  • a high RSSI from a BST implies a reliable TOA. This is due to two facts:
  • the RSSI from a BST usually consists of received signal power + received noise power.
  • the noise in the receiver is mainly thermal noise which is a function of bandwidth and temperature.
  • thermal noise is approximately the same in both BSTs.
  • a high RSSI implies a high received signal power, which in rum implies a high SNR. This is desirable.
  • a higher RSSI usually implies less shadowing than a lower RSSI, which in turn implies less multipath. This is also desirable.
  • Blunders are gross errors in the TOA observations. They may be caused by large signal level fluctuations due to either flat fading or sudden in-band interference. If unremoved, blunders cause large errors in the estimated position. It is possible to detect observations containing blunders by observing the misclosure of each observation during the iterative Least Squares process. Misclosure is defined as the value of the position model (35) given the best available position estimate (x,y). Observations containing blunders will generally have much larger misclosures than observations not containing blunders. When an observation is detected as having a blunder it may be removed from the position estimation process.
  • HDOP generally increases std(R j ).
  • blunder detection is crucial as part of the method of n inimizing positional error by maximizing the number of BSs (and hence reducing HDOP) without incurring a large penalty on std(R; ).
  • poor SNR at some BSs may cause the inclusion of measurements from those BSs to increase the std(range) more than they reduce HDOP.
  • the number of BSTs used to locate the MR is optimized such that stdfjposition) in (36) is minimized.
  • Blunders may also be detected by a statistical analysis of the observation residuals computed from the Least Squares process as shown by Vanicek, P., Krakiwsky, E., "Geodesy: The Concepts," North-Holland Publishing Company, Amsterdam, 1982.
  • the residual of each TOA observation may be standardized by its own estimated standard deviation such that the entire set of residuals is assumed to belong to the normal distribution with zero mean and unit standard deviation. If this hypothesis is correct, the standardized residuals should fall within some specified confidence region for the standard normal distribution. An observation whose residual is flagged is suspected of containing a blunder.
  • the 2-D positioning model for hyperbolic multilateration is
  • N in (40) is the number of BSTs. Note that the TOA from the antenna of the first BST is subtracted from all other TOAs.
  • the unsealed observation covariance matrix is
  • Weighting of the TDOA observations is possible. However, because the TDOAs are a function of two TOAs, the method of assigning weights is no longer straightforward.
  • Blunder detection may also be performed in TDOA positioning.
  • misclosures and residuals correspond to TDOAs. Therefore, a failing misclosure or residual may be due to a blunder in either of the TOA observations from which the TDOA is derived. It is not always possible to isolate the offending BST.
  • the probability of exact solution bifurcation is extremely small. It is, however, possible for bifurcation to exist for a subset of two TDOAs. In such a case, the second solution may affect the final solution obtained when using all available TDOAs.
  • the iterative Least Squares position estimation algorithm may converge to either of the solutions.
  • the solution converged to is a function of the initial starting position used to begin the iterative Least Squares process.
  • the initial position used to begin Least Squares must be relatively accurate. Given no a priori information about the location of the MR, a closed-form position estimation algorithm, using the TOA or TDOA observations, is the only choice.
  • Bancroft's method Bancroft, S., "An algebraic solution of the GPS equations," IEEE Transactions on Aerospace and Electronic Systems, Vol. AES-21, No. 7, Jan. 1985, pp. 56-59
  • Chan and Ho Chan, Y.T., et al., "A Simple and Efficient Estimator for Hyperbolic Location," IEEE
  • LOCA is used in the preferred embodiment of the invention.
  • LOCA is the mathematical dual of hyperbolic trilateration.
  • the fundamental theorem of LOCA states that TOA differences from three BSTs of known location yield a straight line of position. This straight line is the major axis of a conic.
  • the three BSTs lie on the conic and the MR, the location of which is being estimated, lies at one of the foci of the conic.
  • LOCA is expanded into Plane Intersection as shown by Schmidt, R.O., "A New Approach to Geometry of Range Difference Location," IEEE Transactions on Aerospace and Electronic Systems, Vol. AES-8, No. 6, Nov. 1972, pp. 821-835.
  • Plane Intersection is equally applicable to the case of redundancy in 2-D positioning.
  • the best available course of action may be to average the coordinates of the participating BSs and use that as the initial position for Least Squares.
  • Another alternative is to use the approximate coordinates of the BST with which the MR is communicating, as the initial position.
  • TDOAs In some cases it is desirable, or even necessary, to use TDOAs as opposed to TOAs. However, as explained above, it is not always possible to isolate blunders when using TDOAs. It is possible, however, to construct a positioning model that uses TDOAs but gives residuals and misclosures for individual TOAs.
  • one BST is chosen as the reference.
  • the assumption is then made that the TOA at the reference site is equal to the time of transmission, x 0 . From this point on, the development follows that of circular multilateration except that the TOA of the reference BST is subtracted from all TOAs (including that of the reference BST) instead subtracting the time of transmission.
  • the positioning model is then
  • R u is the range from the antenna of the reference (first) BST to the position of the MR.
  • the model of (42) has the same structure as that of circular multilateration (equation(35)) except that the third unknown is R u and not the time of transmission.
  • the model of (42) gives N observations, and therefore N residuals and misclosures, for N BSTs.
  • the N observations are a 0 and N-l TDOAs.
  • the structure of the equations is that of TOA positioning whereas the observations are that of TDOA positioning. Hence, the name Hybrid TDOA Positioning.
  • the design matrix for hybrid TDOA positioning is
  • the unweighted and unsealed observation covariance matrix, C, is the same as that for TOA positioning, the identity matrix of dimension N.
  • DOPs may be calculated, as before, from [AJCJ -1 A] .
  • the observations are, strictly speaking, TDOAs, they may be treated as observations corresponding to the individual BSTs. That is, the first observation, which is always zero, may be viewed as the observation of the reference (or first) BST.
  • the second observation consisting of the TOA of the reference BS subtracted from the TOA of the second BST, is treated as the observation of the second BST, and so on. Therefore, the N diagonal elements of C, may be weighted for the individual BSTs.
  • N independent misclosures and residuals are available for analysis. Blunder detection may then be performed as in the case of TOA positioning and blunders in individual TOA observations may be detected.
  • A is the clockwise angle from noisy to the line joining the MR to the antenna of the ith BST (similar to ⁇ j in equations (32) and (33); not to be confused with ⁇ , in equations (15), (16), (21), (22) and (23)).
  • Least Squares may be used to obtain a unique solution.
  • the relative positions of BSTs with respect to each other and the MR to be positioned are also important for AOA positioning.
  • DOP may again be used to quantify geometry. From the 2-D positioning model (44), the design matrix is found to be
  • the unsealed and unweighted observation covariance matrix, C, is the identity r _ j " i — i matrix of dimension N.
  • the DOPs can again be found from [A C, Aj .
  • the individual AOA observations may be weighted through the observation covariance matrix C,. Those AOA observations deemed more reliable are assigned a smaller variance, or conversely, a heavier weight.
  • AOA observations may contain blunders.
  • statistical testing of the observation residuals and misclosures may be used to detect observations which contain blunders.
  • the estimated AOA is the elevation angle.
  • Diversity observations may be combined by well- known methods including selection combining, maximal ratio combining, co-phasing combining, equal gain combining, or other methods.
  • Least Squares allows for the combination of different types of observations. In particular, it is possible to estimate 2-D position using both AOA observations and TDOA observations within a single Least Squares adjustment.
  • a combination of the two different observation types is particularly useful in situations where only two BSTs are available for estimating the 2-D position of a MR.
  • two BSTs are used to estimate the incoming signal AOA for these BSTs, numerous AOA ambiguities are possible as shown in equation (15). This is due both to the nature of AOA estimation with a 1-D linear array, and to the fact that the spacing between antennas can be greater than one wavelength. Therefore, with two independent AOA observations only, many position solutions are possible and there is no way of determining which of the many solutions is correct.
  • each AOA combination is combined with the TDOA observation in Least Squares. That combination of AOAs which results in the smallest sum of squares of residuals is chosen as correct. The corresponding position solution is used as the MR position estimate.
  • both the model for AOA positioning (44) and the model for TDOA positioning (39) are used.
  • the design matrix A for two independent AOA observations and one TDOA observation is cos ⁇ ! - sin ⁇ , cos A 2 - sin ⁇ 2 (46) x - x,
  • d is the best derived distance between the MR and the i ⁇ BST.
  • the first two rows correspond to the two AOA measurements whereas the third row corresponds to the TDOA measurement.
  • the two columns correspond to the two unknowns, x and y. Additional AOA and TDOA measurements may be included by adding appropriate rows to (46). Note that the 1/c factor seen in (40) is missing in the TDOA measurement row of (46). This is done such that the units throughout A are dimensionless.
  • the misclosure vector necessary for the Least Squares adjustment mechanism, merely consists of the misclosures of all observations.
  • the misclosure vector corresponding to the design matrix of (46) is
  • the unsealed observation covariance matrix is the identity matrix of appropriate dimension (number of AOA observations plus the number of TDOA observations).
  • DOP may again be used to quantify geometry.
  • the DOPs (HDOP, EDOP,
  • NDOP can be found from ATC ⁇ A immediately above.
  • Y where A and C, are defined Weighting of the AOA/TDOA Observations
  • the individual AOA and TDOA observations may be weighted through the observation covariance matrix C, . Those AOA and or TDOA observations deemed more reliable are assigned a smaller variance, or conversely, a heavier weight.
  • Both the TDOA and AOA observations may contain blunders.
  • Statistical testing of the observation residuals and misclosures may be used to detect observations which contain blunders.
  • the range, R j is the observed parameter.
  • AOA and range positioning may also be combined in order to estimate the position of the MR with three BSTs.
  • the design matrix A for two AOA observations and one range observations is
  • the misclosure vector corresponding to this design matrix is
  • the unsealed observation covariance matrix is the identity matrix of appropriate dimension (number of AOA observations plus the number of range observations).
  • the design matrix A for one AOA observations and two range observations is cos A, sin ⁇ .
  • the misclosure vector corresponding to this design matrix is
  • DOP may again be used to quantify geometiy.
  • the DOPs (HDOP, EDOP,
  • NDOP NDOP
  • the individual AOA and range observations may be weighted through the observation covariance matrix C, . Those AOA and/or range observations deemed more reliable are assigned a smaller variance, or conversely, a heavier weight.
  • Both the AOA and range observations may contain blunders.
  • Statistical testing of the observation residuals and misclosures may be used to detect observations which contain blunders.
  • the 3-D or 2-D location of the moving MR must be estimated at various epochs of time.
  • the doppler shift of the signal transmitted from each BST may be estimated as discussed earlier.
  • ⁇ f ( , the frequency offset for the i ⁇ BST, is assumed to be known by the RR.
  • the angle A in (52) may be calculated for each BST.
  • the knowns are f j , ⁇ , ⁇ f ⁇ and ⁇ j , (where f i is measured or observed), and v, ⁇ and ⁇ f 0 are the unknowns.
  • Three BSTs are required - the same number required to estimate the 2D MR position using TOA or TDOA positioning.
  • Geometry affects the estimation of speed and DOT as well. For instance, it is intuitively obvious that when the MR is traveling on a line perpendicular to the line connecting it and a BST, no information regarding the speed of the MR is available from observed data at that BST.
  • the design matrix for the model of (52) is The DOPs can again be found from A T C, -I AJ . In this case the DOPs will be speed DOP, direction of travel DOP, and frequency offset DOP.
  • C is the unsealed and unweighted identity matrix of dimension N.
  • a handset-based WLS it is desirable to be able to use a cellular BS as a BST and to request a CT to act as an MR in order to fulfill the FCC mandate for E911 as well as to offer location services to the current 60 million cellular subscribers in North America who are expected to grow to 120 million by the year 2001. This however might require a change to the CT as well as a change in the cellular standards. If the required changes preclude existing CTs, the FCC mandate for E911 will not be fulfilled using such a system. Nonetheless, a handset-based solution is attractive since it can offer an accuracy comparable to the one offered by systems dedicated for location such as GPS and Loran-C. This is why it is still acceptable to make changes to the CT regardless of the FCC mandate.
  • a hardware modification to the existing CT is required when a required change in the signal processing algorithms cannot be accomplished in software. This may be
  • Improvement 0 The most common method to estimate the TOAs, POAs, Strength of Arrival (SOAs), and FOAs of received signals at the MR consists of correlating a received signal with a clean version of itself.
  • Conventional methods cannot resolve correlation peaks beyond the conventional time- domain Rayleigh resolution.
  • conventional methods either choose the highest correlation peak above a certain threshold or the peak corresponding to the first arrival above a given threshold. Over a multipath channel, it is possible to assume that the first peak above a certain threshold corresponds to the correct peak, however, aliasing, noise, and interference can cause spurious peaks to appear in time before the correct one. It is an object of Improvement O to
  • Such an optimization is based on minimizing HDOP.
  • the optimized parameters: threshold and list of BSs, are then transferred to the CT whenever it is requested to enter into Location state 0.
  • the CT is forced to receive more signals from more BSTs in order to reduce HDOP.
  • the optimization process can be explained as follows. In order for the CT to receive more signals, its threshold has to be lowered. A lower threshold however forces the CT to accept weak signals. This produces low SNR signals that might suffer from multipath and interference.
  • the optimization process consists of solving for the position of the MR using Least Squares weighted by the individual RSSI from each received signal. 0224
  • Improvement I It consists of using the SR and Inverse SR algorithms over the correlation function (that is obtained by correlating the received signals at both the MR and RR with a clean version of themselves) in order to:
  • Location state I consists of selecting the correlation peak that produces the smallest residual when performing a weighted Least Squares operation on all possible correlation peaks resulting from super-resolving the correlation functions produced from the received signals using the inverse SR algorithm. It is also an object of Location state I to optimize the list of BSs considered for correlation based on minimizing HDOP.
  • modification H existing methods require a great deal of standard change as well as CT change. More specifically, it consists of receiving GPS signals partly at the CT and partly at the RR. Instead, several improvements over the existing methods are offered in modification H:
  • CT it allows the CT to choose between receiving cellular signals, or receiving non-cellular signals (such as GPS, Loran-C, etc.), or both;
  • CT it allows the CT to use short messaging cither for transmitting its Location Information or for receiving the Reference Information, so as not to be network-dependent;
  • an IS-95 receiver which measures the TOA of multiple pilot signals for use in estimating the location of said receiver.
  • This receiver may be obtained through modifications to the operation of a wide variety of CDMA receivers, including CDMA CTs.
  • An IS-95 CT is directed to enter location state 0 or Location State I by an order received from:
  • the CT assembles a list of pilots to search for (1011) as shown in Figure 10.
  • This list may be formed from a combination of the CT's Active set, neighbour set, remaining set, and a list of pilots specified in a message received by the CT from a remote source.
  • Block (1012) initializes P to the first pilot in the search list.
  • a searcher element is assigned to search for pilot P in the received signal (1013).
  • a window of PN offsets start and end, or center and width
  • the searcher then computes received energy of pilot P at multiple delay offsets within the search window. In the preferred embodiment, these delay offsets are spaced one half chip apart (0.4069 microseconds). For each delay offset, the searcher accumulates received pilot energy over a specified duration of time, (preferred embodiment uses 0.5 ms)
  • the searcher results are examined for evidence of reception of a pilot signal with sufficient quality (1014).
  • the ratio of the average pilot chip energy to the total received power (called EJI is computed at each delay offset.
  • Figure 11 illustrates a possible set of values for searcher results.
  • significant pilot energy appears near delay offset 7.
  • the pilot EJl 0 at each delay offset is compared to a threshold value (-15 dB in the preferred embodiment). If the fl ⁇ at one or more delay offsets exceed the threshold, then the earliest such delay offset is used as an initial estimate of the time of arrival of the pilot. If the threshold is not exceeded at any of the delay offsets, men it is determined that the pilot was not detected.
  • pilot P was detected, then a demodulator finger is assigned to pilot P at the initial estimate of the time of arrival of said pilot (1015). The finger is allowed to track pilot P for a period of time (preferred embodiment uses 2 ms), after which
  • search results are assembled (1019) for subsequent processing or transmission to a remote site.
  • Figure 10 illustrates steps performed sequentially. Some steps could be performed in parallel.
  • the searcher can begin to search for another pilot while a demodulator finger tracks a previously found pilot. If the CT's time reference is not sufficiently stable, then one demodulator finger can be assigned to one pilot while other pilots are searched and tracked by other fingers, in order to maintain a known time reference.
  • a demodulator finger may track a pilot using a delay-locked loop, as described in U.S. Patent 5,764,687. This results in the finger moving to the peak of the arrival of the pilot. Improved performance may be obtained by having a finger track towards a predetermined E/L 0 threshold on the rising edge of the pilot energy versus delay offset.
  • the location operation is performed at the Host which can be either located within the MR (Option I) or remotely from the MR (Options H and DT) (e.g. within the RR or a TR).
  • Option I the MR
  • Option H and DT the MR
  • TDOA a received signal
  • r;(t) that is transmitted from the i ⁇ BST
  • the required processing consists of correlating the received signal with a clean version of p ; (t) and sending the information regarding the correlation peaks to the Host.
  • the required processing consists of correlating the received signal with several clean versions of the p ; (t) at different frequencies and sending the information regarding the correlation peaks to the Host.
  • the correlation peaks between r ; (t) and a clean version of p ; (t) can be super- resolved using SR algorithms. It is sometimes preferable to perform the SR algorithms at the Host rather than performing them at the MR since it requires a number of computationally intensive operations that are easier to develop at the Host than at the MR. This implies that the correlation peak needs to be transferred to the Host which does not represent a large amount of information and can be easily transferred to the Host.
  • Procedure I In order to be able to correlate the received signal, r,(t), at the MR, a clean version of p,(t) must be available at the MR beforehand. In order to avoid correlating r,(t) with p,(t), one can use either one of the following two procedures: Procedure I:
  • TOA and AOA information are contained in the phase of the tone while the FOA information are contained in the frequency of the tone and can be extracted using Procedure HI (see below).
  • the nonlinear operation in Procedure H can be either a quadratic operation (e.g. squaring the signal), a delay product operation (i.e. multiplying the signal by a delayed version of itself) or any other suitable nonlinear operation.
  • Procedures I and H can be used both for acquiring the signal and for estimating its TOA. This is accomplished by stamping the time using a reliable clock at either the start of the generated tone or obtaining the time stamp at the end of the observation window (the latter is preferred).
  • Procedures I and H can be improved by resolving the existing or generated tones using SR and inverse SR algorithms. Moreover, the filtering required in step 1 of Procedure II can be carried out in the digital domain in order to reduce the effect of temperature and aging on group delay variation.
  • Procedures I and H can be improved by repeating them a number of times and then either averaging the Location Information or averaging the TDOA, PDOA or FDOA between pairs of BSTs.
  • An alternative method for improving Procedures I and ⁇ is by substantially increasing the duration of observation, T ; , of the received signal, r,(t), thereby increasing the received SNR.
  • T duration of observation
  • r received signal
  • p the received signal
  • p has to be available at the MR.
  • An important feature in Procedure II is the fact that the signal r,(t) transmitted by the i* BST does not have to be a deterministic signal, i.e. a clean replica of the signal does not have to be available at the MR.
  • Procedure II Equivalently, an important feature in Procedure II is the fact that the probability of location is 100%, unlike Procedure I which is below 100%.
  • T j5 can accommodate different situations. For example, in rural areas, where the cells are large, the SNR can be increased by increasing the duration of observation. Or equivalently, if the geometry is not good (i.e.
  • HDOP is large, e.g. highway
  • the duration of observation is thus preferably optimized with respect to the SNR and HDOP.
  • this also leads to reducing the effect of multipath, particularly in urban centers.
  • narrowband digital filters of a few hertz bandwidth are used (e.g. Fourier Transform-based filters, sliding Fourier Transform-based filters, ITR filters, FIR filters, etc.).
  • the Fourier Transform-based filter has a bandwidth which is a function of the observation duration, T ⁇ .
  • the ITR filter has a steady state response time which is a function of T ; .
  • the observation duration, T ; has to be selected in order to optimize the SNR (increase it) versus the bandwidth of the filter (decrease it) versus the response time of the wireless location system (minimize it).
  • MPj(t) In order to reduce the effect of the multipath, MPj(t), on the performance of a kinematic wireless location system, it is desirable to resolve the different Doppler frequencies ⁇ (which correspond to different TOAs and different AOAs) and to choose the frequency corresponding to the first arrival.
  • the method of resolving the different frequencies can be accomplished using SR algorithms.
  • DSP Digital Signal Processor
  • Initial rejection is achieved at RF using an RF Band Pass Filter (BPF) (802), followed by an RF amplifier (803) for good sensitivity.
  • BPF Band Pass Filter
  • Subsequent rejections are achieved at each TF stage which consist of a mixer (804) followed by an Intermediate Frequency (W) BPF filter (805) and an IF amplifier (806).
  • the sampling rate is intended to be below the Nyquist rate in order to force aliased copies of the IF signal (807) close to a lower digital IF, f 2 .
  • DDC Digital Down Converter
  • step 1 can be performed using a discrete Carrier Recovery Loop (CRL) described in Procedure IV (see below).
  • CTL Carrier Recovery Loop
  • step 2 can be performed using a Direct Digital Synthesiser (DDS) (809).
  • DDS Direct Digital Synthesiser
  • step 3 can be performed using a gradient-type search algorithm or any other suitable search algorithm.
  • step 4 can be performed using cither a DDS (809) only or a
  • step 1 in Procedure IV processes a new segment of the signal, r,(t), every iteration. If a DDS (809) and a DDC (809) are used is step 4, then step 1 in Procedure IV processes the same stored segment of the signal, r,(t), every iteration until the DDS (809) is used. It is preferred to use both the DDC and the DDS in order to reduce the effect of finite precision arithmetic on the estimation process.
  • the DSP (817) in step 1 of Procedure IV of the carrier recovery loop can consist of a digital notch filter centered at the desired IF frequency f 2 , whose power of response is to be minimized, i.e. f ⁇ - is the response of the digital notch filter.
  • the notch filter can consist of
  • the Band Pass filters can be implemented using (real) digital Biquadratic filters.
  • the DSP (817) in step 1 in Procedure IV of the carrier recovery loop can consist of a digital Band Pass Filter centered at the desired IF frequency f 2 , whose power of response is to be maximized, i.e. f ⁇ - is the response of the Band Pass Filter.
  • the BPF can be implemented as a (complex) 1 st order (or higher) digital HR filter.
  • the signal to process is a baseband signal (816) and the digital BPF turns into a simple (complex) LPF.
  • a simple (complex) LPF is an averaging filter which adds all real samples together and all imaginary samples together over the observation interval.
  • a common reference is used for all the LOs and for all the clocks.
  • a suitable common reference for the LOs and for the clocks can be derived from the GPS signal. Given that either all LOs in the BSTs have a common reference, which implies that
  • a handset-based Wireless Location System where an MR can be located without the need for an RR.
  • WLS Wireless Location System
  • the advantage is that such a sector is usually time- referenced to other sectors belonging to the same BST, i.e. sectors in the same BS are more or less synchronized to one another.
  • the MR can estimate the AOA of signals transmitted by such sectors without a need for an RR.
  • the MR can estimate its location using only two independent BSs based on the individual estimated AOA from each BS without the need for an RR_
  • a wireless location system has to be calibrated prior to operation (and occasionally during operation) in order to remove the effect of the overall Group Delay, ⁇ tj, at each RR. Moreover, temperature and aging can cause the overall Group Delay to change with time. Therefore, it is necessary to calibrate the RR on a regular basis.
  • a method to calibrate an RR is as follows: Procedure V:
  • Procedure V can be repeated a number of times in order to average the relative group delay over time.
  • an MR is stored in proximity to the RR. This allows for automatic calibration without a need for human intervention by forcing the MR to act as a transmitter.
  • SAW filters As the BPF (803, 805, 811) in Figure 8. It is also recommended to place the SAW filters (803, 805, 811) in a metallic block to minimize flexing effects due to temperature variations.
  • Figures 9a, 9b, 9c and 9d illustrate the preferred embodiment for adjusting the observed TOAs such that the sum of the squares of the residuals is a minimum.
  • Other methods for adjusting the observed according to some optimization criteria are available mcluding:
  • the positioning process begins with a set of TOA observations from a minimum of three BSs. If only three TOAs are available (901), solution bifurcation is tested for (902). In the event of solution bifurcation, the LOCA algorithm (903) is performed on the TOA observations yielding two solutions. These two solutions are each used as the initial position for Least Squares (904 and 905) resulting in two Least Squares solutions for the MR.
  • the Least Squares algorithm here, and in all other parts of Figure 9, is Hybrid TDOA positioning with received signal strength used to weight the individual BS observations. During every iteration of Least Squares, misclosures for each
  • both Least Squares solutions (904 and 905) converge (906) but are not equal (907)
  • both Least Squares solutions are reported as possible positions for the MR (908). If the either of the LS solutions diverge (906) or the convergent solutions are equal (907), the two LOCA solutions from (903) are reported as possible positions of the MR (909).
  • Plane Intersection (912) (or any other closed-form position estimation algorithm) is executed.
  • the resulting solution is used as the initial position for Least Squares (913).
  • a check is made as to whether Least Squares (913 or 911) converged (914). If Least Squares converges, another check for observational redundancy is made (915). If there is no redundancy, the Least Squares solution is reported as the position estimate of the MR (923). If there is redundancy, the normalized residuals are statistically tested for normality (916). The Least Squares solution is reported (923) should all residuals pass. If any residuals fail, the redundancy numbers of the failing observations are checked (917). The redundancy number of the i ft observation is defined as
  • C ? is the covariance matrix of the residuals and is defined as
  • Least Squares is executed again but with the average of participating BSTs coordinates used as the initial position (918). If Least Squares now converges (919) but there is no observation redundancy (920), the newly convergent Least Squares solution is reported as the position estimate of the MR (923). If there is redundancy (920), the standardized residuals are tested for normality (921). Should all residuals pass, the newly convergent Least Squares solution is reported (923). If some of the residuals fail but none of the corresponding observations have redundancy numbers greater than some threshold (922), the newly convergent Least Squares solution is reported (923). If observations with failing residuals have redundancy numbers larger than the threshold (922), that observation with redundancy number greater than the threshold and with the largest standardized residual is permanently removed from the observation set (928). The process then continues from (928) as described above.
  • the communication between the Host and the customer can be achieved over it.
  • a MR calls 911 for emergency
  • a powered-on MR is to be located by a customer
  • its positional information can be relayed from the Host to the customer through the internet.
  • MR location information residing at the Host can be integrated with third-party information (i.e. a map database, or a database consisting of geo-coded business addresses such as restaurants, towing companies, etc.) resid g in some other geographical location, perhaps even in another country.
  • third-party information i.e. a map database, or a database consisting of geo-coded business addresses such as restaurants, towing companies, etc.
  • the combination of this MR location information with the third party information can be transferred over the internet to allow customers to locate a particular MR wrt to either a map location (i.e.
  • the service provider will integrate together the various components including the MR location information derived from the wireless location system and the third-party information database.
  • the integration process of the different databases is transparent to the customer.
  • the customer will only know that the location-based service exists as such to bring all of the different components together to provide for a complete service which can be offered either on a regional, national or perhaps global basis.
  • these location-based services can now be offered on an economical basis to the customer. Examples of such services are fleet management, concierge services, roadside assistance, child find services, etc. Legahty and security aspects are a concern on the internet and in some cases a dedicated link may be necessary.
  • a clean replica referred to herein is a copy of a signal that is not noisy, as for example may be obtained by a stored version of the signal transmitted.
  • Low group delay variation preferably means a few 10s of nanoseconds of group delay.
  • An adequate objective function is a function known in the literature, for example minimizing the mean squared error.
  • a PSAP is a public safety answering point, for example a 911 station.
  • An RSSI is a receiver signal strength indicator.

Abstract

The invention consists of methods and apparatus to estimate the position and velocity of a Mobile Receiver (MR) using either the Time Of Arrival (TOA) of signals received by the MR, their Phase Of Arrival (POA), their Strength Of Arrival (SOA), their Frequency Of Arrival (FOA), or a combination thereof, with respect to a reference produced by a Reference Receiver (RR) of known location. In order to solve for the coordinates of the MR, the invention uses either hyperbolic multilateration based on Time Difference Of Arrival (TDOA), or linear multiangulation based on Phase Difference Of Arrival (PDOA), or both. In order to solve for the velocity of the MR, the patent uses FOA based on Frequency Difference Of Arrival (FDOA). An important contribution of this invention is the way the MR receives, processes and combines available signals for location purposes. Another important contribution of this invention is the way the RR receives, processes and combines available signals for reference purposes. Yet another important contribution of the invention is the application of Super-Resolution (SR) techniques at both the MR and the RR to increase the resolution of the estimated TOAs, POAs, SOAs, or FOAs.

Description

Methods and Apparatus to Position a Mobile Receiver using Downlink Signals
Field of the Invention: This invention relates to location finding and tracking of a Mobile Receiver (MR) using a Reference Receiver (RR) of known location.
Background of the Invention:
Whenever a number of signals are available for reception to the MR as well as to the RR, it is possible to locate the MR. Examples of radio signals that exist with worldwide coverage are satellite signals such as Global Positioning System (GPS), Iridium, Globalstar, etc. Examples of radio signals that can penetrate structures such as buildings are land-based systems such as cellular, PCS, paging, and ESMR, among others. Both satellite-based signals and terrestrial-based signals can be used to locate the MR by relying on their reception at both the MR and the RR. In this patent, we refer to such a location system as a handset-based Wireless Location System (WLS) (or Forward/downlink-based WLS), as opposed to a network-based WLS (or Reverse/uplink-based WLS), which locates a Mobile Transmitter (MT) using a plurality of receivers that receive the signal transmitted from the MT. An advantage of a network-based WLS is that it can be passive (in that it can locate MR's without modification or addition to wireless transmissions) and that no modifications to the hardware or operation of the MT are required for the WLS to locate the MT. On the other hand, a handset-based WLS has the advantage of requiring no modification of existing network infrastructure. The location in the handset-based WLS can be performed either at the MR, or at the RR, or at a Third Receiver (TR). When the location is performed at the MR, the RR is required to transmit the Reference Information to the MR. When the location is performed at the RR, the MR is required to transmit the Location Information to the RR. When the location is performed at the TR, both the MR and the RR are required to transmit their corresponding information to the TR.
Radio signals can interfere with one another unless multiple access techniques are employed. Existing radio signals use any one of a number of multiple access techniques such as: Code Division Multiple Access (CDMA), Time Division Multiple Access (TDMA) and Frequency Division Multiple Access (FDMA). CDMA is used in IS-95, the North American standard for CDMA CTs, as well as in GPS. FDMA is used in AMPs, the North American standard for analogue CTs. TDMA is used in IS- 136, the North American standard for TDMA CTs, and in GSM, the European standard for TDMA CTs. Without loss of generality, we will use in this patent GPS signals as an example of a satellite-based radio signal and IS-95 signals as an example of a land- based radio signal. This choice does not preclude any other standard or any other multiple access technique. In this patent, we refer to a transmitter as a Base Station Transmitter (BST). It is important to note that a RR is not required in this patent if the BSTs are referenced to the same reference signal.
The need for wireless location finding and tracking of MRs is ever increasing. Some of the possible services for locating MRs are:
1. Enhanced Emergency Services: In recent years pressure has been mounting for the development of technology to position Cellular Telephones (CTs). The primary driving force has been enhanced 91 1 (E91 1) services for wireless telephone subscribers. E911 services provide the 911 operator with information such as calling number, street address, and the primary subscriber's name. A ruling concerning E911 emergency calling systems was released by the FCC in October, 1996 (CC Docket No. 94-102) and clarified in December 1997. In the document, the FCC requires that by October 2001, a wireless 91 1 caller be located with a horizontal accuracy of 125 meters RMS (Root Mean Square). 2. Tracking of Fraudulent Calls.
3. Tracking of Stolen Vehicles: Conventionally, to track stolen vehicles requires installing an RF tag such as in Bird, U.S. Patent No. 5,418,537 issued May, 23, 1995, leaving it permanently on in each vehicle to be tracked and a new infrastructure for the BSs throughout the desired service area. An alternative solution is to either use existing CTs to take advantage of the existing cellular coverage and the broad availability of inexpensive CTs, or to use GPS receivers to take advantage of the availability and accuracy of satellite signals.
4. Fleet Management for Courier and Transportation Businesses: Once again, to manage a fleet requires installing an RF tag such as in Song, U.S. Patent No. 5,208,756 issued May, 4, 1993, and Sheffer et al., U.S. Patent No. 5,218,367 issued June 8, 1993. Instead, existing CTs can be used to take advantage of the existing cellular coverage and the broad availability of inexpensive CTs. Alternatively, GPS receivers can be used to take advantage of the availability and accuracy of satellite signals. 5. Location-sensitive billing.
Without loss of generality, we will assume throughout this document that the intention of the handset-based WLS is to estimate the horizontal position of the MR as well as its horizontal velocity. In case both the vertical and the horizontal position of a MR are to be estimated, an extra independent measurement is required in addition to the minimum number required for horizontal positioning.
Objects of the Invention:
According to an aspect of the invention, there is provided a handset-based Wireless Location System (WLS) whereby an MR can be located using signals received by both the MR and an RR. More specifically, it is an object of the invention to estimate the static and kinematic positional information of a MR which receives signals transmitted from a number of BSTs.
According to a further aspect of the invention, there is provided a system that corrects for the sources of error that affect the different location technologies for a handset-based WLS in a unique and novel way.
According to a further aspect of the invention, there is provided a system that solves for the position of the MR using Least Squares where misclosures and standardized residuals are used to flag the observations that might have a blunder. According to a further aspect of the invention, there is provided a system that uses Chaffee's method or Location On the Conic Axis (LOCA) to detect solution bifurcation.
According to a further aspect of the invention, there is provided a system that uses LOCA and/or Plane Intersection to provide an initial position for Least Squares. According to a further aspect of the invention, there is provided a system that uses individually, or in combination, TDOA observations, hybrid TDOA observations and AOA observations to locate the MR.
According to a further aspect of the invention, there is provided a system that resolves the AOA ambiguity at a cellular BS by using the TDOA information from two distant BSs or the AOA from several sectors at the same cellular BS.
According to a further aspect of the invention, there is provided a system that uses the range information from two BSTs or the AOA from several sectors at the same cellular base station to resolve the AOA ambiguity at such a BST. Objects of the Invention:
According to an aspect of the invention, there is provided a handset-based
Wireless Location System (WLS) whereby an MR can be located using signals received by both the MR and an RR. More specifically, it is an object of the invention to estimate the static and kinematic positional information of a MR which receives signals transmitted from a number of BSTs.
It is also an aspect of the invention to provide a system which corrects for the sources of error that affect the different location technologies for a handset-based WLS in a unique and novel way. More specifically, it is proposed to correct for frequency errors, clock offsets, overall group delays, multipath and interference.
Yet another aspect of the invention is to estimate the Doppler shift (δf,) in each received signal at the MR in order to estimate the speed and Direction of Travel (DOT) of the MR using Frequency Difference Of Arrival (FDOA). This applies to modifications: I and II. Yet another aspect of the invention is to minimize the use of RRs in order to reduce the cost and complexity of a handset-based WLS. It is an object of the invention to produce a handset-based Wireless Location System (WLS) where an MR can be located without the need for an RR. This applies to all three modifications: 0, 1 and II.
Objects of the Invention:
According to an aspect of the invention, there is provided a handset-based Wireless Location System (WLS) whereby an MR can be located using signals received by both the MR and an RR. More specifically, it is an object of the invention to estimate the static and kinematic positional information of a MR which receives signals transmitted from a number of BSTs.
Another contribution of this invention is the way a CT is requested to act as an MR in a handset-based WLS where a transmitter in a cellular BS is considered to be a BST. One should keep in mind that for a CT to act as an MR, a certain level of CT modification might be required. Since there is a large number of existing CTs worldwide, it is important to keep in mind that some levels of modification to the CT might preclude existing CTs (which is not a desirable outcome). We attempt to reduce the amount of CT modification for a desired WLS performance level. Three levels of CT modification have been identified in this patent: modification 0 which requires no modification to the CT; modification I which requires only software modification to the existing CT; and modification II which requires both software and hardware modifications to the existing CT.
It is an aspect of the invention to improve modification 0, 1 and II over existing methods.
It is an aspect of the invention to provide a system that optimizes the threshold used in selecting the correlation peak, and to optimize the list of signals considered for correlation, for all modifications: 0, I and II. It is a further aspect of the invention to optimize both the threshold as well as the list of signals considered for correlation based on minimizing HDOP.
Yet another aspect of the invention is the use of the SR and Inverse SR algorithms over the correlation function (that is obtained by correlating the received signals at both the MR and RR with a clean version of themselves) in order to: resolve the correlation peaks beyond the Rayleigh resolution and to avoid using a threshold in selecting the correlation peak.
According to another aspect of the invention, there is provided a system that optimizes the list of signals considered for correlation based once again on minimizing HDOP. This applies to modifications: I and II.
Yet another aspect of the invention is a system for locating an MR by estimating the Phase Of Arrival (POA) of existing or generated tones contained in the received signals. The generated tones can be obtained by performing a nonlinear operation over the received signals at the MR. The phases of the existing or generated tones can be extracted using SR algorithms in order to reduce the effect of multipath by selecting the tone corresponding to the first arrival in time. This applies to modifications: I and II.
Yet another aspect of the invention is use of the estimated POA to estimate either the TOA of the received signals or their AOA using Phase Difference Of Arrival (PDOA). This applies to modifications: I and II.
Yet another aspect of the invention is provision of a system which estimates the Doppler shift (δf) in each received signal at the MR in order to estimate the speed and Direction of Travel (DOT) of the MR using Frequency Difference Of Arrival (FDOA). This applies to modifications: I and II.
Yet another aspect of the invention is a system that minimizes the use of RRs in order to reduce the cost and complexity of a handset-based WLS. According a further aspect of the invention, there is provided a handset-based Wireless Location System (WLS) where an MR can be located without the need for an RR. This applies to all three modifications: 0, 1 and II.
According to a further aspect of the invention, there is provided a system that resolves the AOA ambiguity at a cellular BS by using the TDOA information from two distant BSs or the AOA from several sectors at the same cellular BS.
According to a further aspect of the invention there is provided a system which uses the range information from two BSTs or the AOA from several sectors at the same cellular base station to resolve the AOA ambiguity at such a BST.
Objects of the Invention:
According to an aspect of the invention, there is provided a handset-based
Wireless Location System (WLS) whereby an MR can be located using signals received by both the MR and an RR. More specifically, it is an object of the invention to estimate the static and kinematic positional information of a MR which receives signals transmitted from a number of BSTs.
According to a further aspect of the invention, there is provided a system that corrects for the sources of error that affect the different location technologies for a handset-based WLS in a unique and novel way. According to a further aspect of the invention, there is provided a system that solves for the position of the MR using Least Squares where misclosures and standardized residuals are used to flag the observations that might have a blunder.
According to a further aspect of the invention, there is provided a system that uses Chaffee's method or Location On the Conic Axis (LOCA) to detect solution bifurcation.
According to a further aspect of the invention, there is provided a system that uses LOCA and/or Plane Intersection to provide an initial position for Least Squares.
According to a further aspect of the invention, there is provided a system that uses individually, or in combination, TDOA observations, hybrid TDOA observations and AOA observations to locate the MR.
According to a further aspect of the invention, there is provided a system that resolves the AOA ambiguity at a cellular BS by using the TDOA information from two distant BSs or the AOA from several sectors at the same cellular BS. According to a further aspect of the invention, there is provided a system that uses the range information from two BSTs or the AOA from several sectors at the same cellular base station to resolve the AOA ambiguity at such a BST.
Definitions:
The following definitions take precedence over definitions for the same terms that can be found in the open literature.
Mobile Receiver (MR): is a device, which can be portable or fixed, that can consist of a receiver alone or both a receiver and a transmitter. It can be a regular Cellular Telephone (CT), a PCS (Personal Communication Systems) telephone, a cordless telephone, a Personal Digital Assistant (PDA), a GPS receiver, or a combination thereof. It can be a radio tag or a wireless telephone that does not contain the audio portion of the telephone. It can also be a transmitter that transmits periodically over a given channel, or a receiver that receives Radio signals, or both. Reference Receiver (RR): is a device of known location with respect to the location of BSTs that are transmitting radio signals intended for reception by both the MR and the RR. An RR can be portable or fixed. It can consist of a receiver alone or both a receiver and a transmitter. The RR estimates the TOA, POA, SOA and/or FOA of signals transmitted by the BSTs. It either transmits the Reference Information to the MR or to the TR, or it receives the Location Information from the MR. In some unique situations when the BSTs are referenced to a common signal, the RR is not required. Reference Information: is the estimated TOA, POA, SOA and/or FOA of the signals transmitted by a number of BSTs at the RR as well as the location of such BSTs. Location Information: is the estimated TOA, POA, SOA and/or FOA of the signals transmitted by a number of BSTs at the MR.
Third Receiver (TR): is a device that receives Location Information from the MR as well as Reference Information from the RR in order to estimate the location of the MR. It can be a receiver only, or a transmitter and a receiver. The Mobile Switching Center (MSC) can play the role of the TR. Base Station Transmitter (BST): is a device of known location with respect to other BSTs at a given time. A BST can be portable or fixed. It can consist of a transmitter alone or both a receiver and a transmitter. It can be a regular cellular Base Station (BS), a regular satellite transceiver, a GPS satellite, a Loran-C transmitter, a dedicated transmitter, a transponder, or any other type of transmitter/transceiver combination. Host: is a processor that processes the Location Information together with the Reference Information to estimate the location of the MR. The function of the Host can be played by the MR, the RR or a TR. The host can also be responsible for location services such as fleet management, location sensitive billing, etc. Network-based WLS (Forward/Downlink-based WLS): locates a Mobile Receiver (MR) using a plurality of BSTs that transmit signals received by both the MR and the RR.
Handset-based WLS (Reverse/uplink-based WLS): locates a Mobile Transmitter (MT) using a plurality of receivers that receive the signal transmitted from the MT. Cell: is a geographical area serviced by a cellular Base Station (BS).
Sectorized cell: is a cell that is made of several non-overlapping sectors. Each sector can be considered as an independent cell to be serviced by an independent cellular BS. However, sectors in a cell are usually serviced by the same cellular BS in order to minimize cost and complexity. We refer to such a BS as a sectorized BS. Pilot Channel: is the channel to be used by each sector to broadcast a (locally unique) pilot signal.
Pilot Signal: is defined as a sequence known to the MR.
As defined in the IS-95 standard, this sequence is a code repeated exactly 75 times every 2 seconds, and is known as the "pilot PN sequence" or the "short code". Although all sectors broadcast the same pilot signal, the transmissions from each base station have different synchronization or phases, which distinguish the pilot signals. Phase increments are specified as multiples of 64 chips, where 1 chip is 1/1.2288 microseconds. The offset in units of 64 chips, is referred to as the pilot PN offset. Pilot increments in a network are typically 3 or 4 PN offsets. A PN offset of 4 corresponds to 256 chips, or 208.3 microseconds.
As defined in the AMPs standard or in the IS-136/GSM standard, this sequence is a Digital Colour Code (DCC) that is repeated at specific intervals of time. Super-Resolution (SR) Algorithm: is an operation that transforms a time domain signal, s(t), to a frequency domain signal, S(f), in such a way that the frequency domain signal, S(f), has a better resolution than the resolution offered by the Fourier transform, i.e. S(f) has a better resolution than the Fourier transform F{s(t)}, of s(t). Vice-versa, a SR algorithm transforms a frequency domain signal, S(f), to a time domain signal, s(t), with a better resolution than the resolution offered by the inverse Fourier transform, i.e. s(t) has a better resolution than the inverse Fourier transform F"' {S(f)}, of S(f)). Examples of SR algorithms are well known in the literature and include: MUSIC/Root MUSIC, ESPRIT, Auto Regressive Moving Average, Minimum Variance, MUSIC using Higher Order Statistics, ESPRIT using Higher Order Statistics, Auto Regressive Moving Average using Higher Order Statistics, or Minimum Variance using Higher Order Statistics.
Frequency-domain Rayleigh Resolution: is the frequency-domain resolution offered by the Fourier transform. Time-domain Rayleigh Resolution: is the time-domain resolution offered by the inverse Fourier transform.
Inverse SR algorithm: is an algorithm which processes a time domain signal in order to improve its time-domain resolution over the conventional time-domain Rayleigh resolution. A preferred embodiment of the inverse SR algorithm comprises a conventional time domain correlator, a time domain window, a Fourier Transform, a frequency domain window, a frequency domain equalizer, or a processor which performs a SR algorithm in order to resolve the TOAs in the received signal at a given MR.
Alternatively, an inverse SR algorithm can process a frequency domain signal in order to improve its frequency-domain resolution over the conventional frequency- domain Rayleigh resolution. In this case, a preferred embodiment of the inverse SR algorithm comprises a frequency domain correlator, a frequency domain window, an inverse Fourier Transform, a time domain window, a time domain equalizer, or a processor which performs a SR algorithm in order to resolve the FOAs in the received signal at a given MR. Sector Antennas: are cellular antennas that exist in some BS that service more than one sector. Each sector has a dedicated antenna (or set of antennas if diversity is applied).
Fourier Transform-Based Filters: are filters that: Fourier Transform the time domain signal, then window the transformed signal over a given band, and Inverse Fourier Transform the windowed signal.
Description of an IS-95 CT receiving Pilot signals:
The CT measures the quality of received pilot signals by computing the ratio between the received energy of a pilot (Ec) to the total received power by the mobile (I0). This measurement is referred to as the Ec/I0 of the pilot. This measurement is closely related to the signal strength of the pilot. An IS-95 CT uses the Ec/I0 measurement of a pilot to estimate the suitability of the communications link between itself (the CT) and that particular sector in a BS. For example, if a CT is in a mode in which it may only communicate with one sector at a time, it may choose to communicate with the sector which has the highest Ec/I0 measured by the CT.
An IS-95 CT may include a searcher element and multiple demodulator fingers as shown in Figure 12. A searcher is used to search for pilot signals by measuring the Ec/I0 of a pilot signal at different delays. A demodulator finger tracks a pilot signal, and may demodulate an information signal associated with said pilot signal.
In addition to transmitting a pilot signal, a base station may also transmit one or more information signals associated with the pilot. These information-bearing signals are transmitted at the same radio frequency as their associated pilot signal and are synchronized with the pilot signal. Thus, an information signal will undergo propagation effects (attenuation, phase shift, time delays, multipath, frequency shift) nearly identical to those of its associated pilot signal. A mobile then tracks the propagation effects of a pilot to aid in demodulating information from an associated signal, or code channel. As well, the mobile can estimate the reliability of information demodulated from a code channel from the Ec/I0 quality measurement of the associated pilot signal.
In general, a CT searches for pilot signals using one or more searchers, and then based on the searcher results, assigns demodulator fingers to demodulate information from code channels associated with suitable pilots. This process is ongoing so that the finger assignment adapts to the changing conditions typical of a CT cellular environment.
During normal operation, the searchers and demodulator fingers in an IS-95 CT are usually assigned according to algorithms which optimize the communications performance of the phone, and the network as a whole. These assignment algorithms will limit the performance of estimation of the CT location, as they will often ignore pilot signals useful for the location estimation. We prefer to force the CT to enter into a state optimized for location.
Summary of the Invention In order to accomplish the aforementioned objects and intentions, there are provided inventive contributions which are claimed in this patent document.
Therefore, According to an aspect of the invention, there is provided a method of estimating the location of a Mobile Receiver (MR) comprising the steps of: receiving a plurality of signals from a plurality of Base Station Transmitters (BSTs) by the MR; and performing a nonlinear operation on the plurality of received signals at the MR in order to generate at least one tone per received signal. Additional refinements of this aspect of the invention are set forth in claims 2-5 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a method of estimating the location of a Mobile Receiver (MR) comprising the steps of: receiving a plurality of signals from a plurality of Base Station Transmitters (BSTs) by the MR; and estimating the location of the MR by estimating the phases of existing tones in each received signal.
Additional refinements of these aspects of the invention are set forth in claims 7- 17 attached hereto, which are hereby incorporated in this summary by reference. According to an aspect of the invention, there is provided a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where one or more antennas is used at the MR and in which location of the MR is achieved using a technique selected from the group consisting of:: AOA positioning (multi-angulation), TOA positioning (circular multilateration),
Range positioning (circular multilateration), TDOA positioning (hyperbolic multilateration), TDOA positioning (circular multilateration with TDOAs), AOA/TDOA positioning (multi-angulation/hyperbolic multilateration), AOA/TDOA positioning (multi-angulation/circular multilateration with TDOAs),
AOA/Range positioning (multi-angulation/circular multilateration), AOA/Range positioning (multi-angulation/circular multilateration), and combining the signals from all antennas for location purposes. Additional refinements of this aspect of the invention are set forth in claims 19, 20 and 22 attached hereto, which are hereby incorporated in this summary by reference According to an aspect of the invention, there is provided a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, using equation (29) and equation (30) (defined below).
According to an aspect of the invention, there is provided a method of estimating the location of a MR wherein location of the MR is achieved based on a solution selected from the group consisting of: solution I, solution II, solution III, and solution IV (as described below). According to an aspect of the invention, there is provided a system to estimate the location of a Mobile Receiver (MR) comprising of: a receiver at the MR for receiving a plurality of signals from a plurality of Base Station Transmitters (BSTs); and a computing means connected to the receiver by a communication linkfor performing a nonlinear operation on the plurality of received signals in order to generate at least one tone per received signal.
Additional refinements of this aspect of the invention are set forth in claims 26 and 27 attached hereto, which are hereby incorporated in this summary by reference According to an aspect of the invention, there is provided a system to estimate the location of a Mobile Receiver (MR) comprising of: a receiver at the MR for receiving a plurality of signals from a plurality of Base Station Transmitters (BSTs); and a computing means connected by a communication link to the receiver for estimating the location the MR by estimating the phases of existing tones in each received signal.
Additional refinements of this aspect of the invention are set forth in claims 2-5 attached hereto, which are hereby incorporated in this summary by reference According to an aspect of the invention, there is provided a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where one or more antennas are used at the MR and wherein the system comprises means for implementing a technique selected from the group consisting of: AOA positioning (multi-angulation), TOA positioning (circular multilateration),
Range positioning (circular multilateration), TDOA positioning (hyperbolic multilateration), TDOA positioning (circular multilateration with TDOAs), TOA/TDOA positioning (multi-angulation/hyperbolic multilateration),
AOA/TDOA positioning (multi-angulation/circular multilateration with TDOAs),
AOA/Range positioning (multi-angulation/circular multilateration), AOA/Range positioning (multi-angulation/circular multilateration), and combining the signals from all antennas for location purposes. Additional refinement of this aspect of the invention is set forth in claim 29 and 31 attached hereto, which are hereby incorporated in this summary by reference. According to an aspect of the invention, there is provided a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, using equation (29) and equation (30).
According to an aspect of the invention, there is provided a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, based on any one of the following solutions: solution I, solution II, solution III, or solution IV.
According to an aspect of the invention, there is provided a system to estimate the location of a Mobile Receiver (MR), the system comprising a Receiver that receives a plurality of signals transmitted from a plurality of BSTs; wherein the Receiver includes an IF-sampling Receiver comprising: a number of IF stages to convert a received RF signal r,(t) to an analog IF signal of IF frequency f i ; a pre-aliasing filter; an Analog-to-Digital (A/D) converter that converts the IF analog signal to an IF digital signal; and a processor for processing the IF digital signal.
Additional refinements of this aspect of the invention are set forth in claims 34- 44 and 49-54 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a system to estimate the location of a MR using a Reference Receiver having a LO, the system comprising: an estimator for estimating the carrier offset (due to Doppler, δf„ and due to LO offset, Δf,) of a plurality of signals transmitted by a plurality of BSTs and received by the RR. Additional refinements of this aspect of the invention are set forth in claims 45- 54 and 60 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a system to estimate the location of a MR having a LO, the system comprising: an estimator for estimating the carrier offset (due to Doppler, δf,, and due to LO offset, Δf,) of a plurality of signals transmitted by a plurality of BSTs and received by the MR.
Additional refinements of this aspect of the invention are set forth in claims 56-59 and 61-66 attached hereto, which are hereby incorporated in this summary by reference
According to an aspect of the invention, there is provided a method of estimating the location of a MR, the method comprising the steps of: receiving a plurality of signals transmitted from a plurality of BSTs at a Receiver; converting the received RF signal r,(t) to an IF signal of IF frequency fi in an
IF-sampling Receiver; filtering the IF signal with a pre-aliasing filter; converting the filtered IF analog signal to an IF digital signal in an Analog-to- Digital converter; and processing the IF digital signal in a signal processor.
Additional refinement of this aspect of the invention is set forth in claim 68 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a method of estimating the location of a MR, the method comprising the steps of: estimating the carrier offset (due to Doppler, Df,, and due to LO offset, Df,) of a plurality of signals transmitted by a plurality of BSTs and received by a RR.
According to an aspect of the invention, there is provided a method of estimating the location of a MR, the method comprising the steps of: estimating the carrier offset (due to Doppler, δf„ and due to LO offset, Δf,) of a plurality of signals transmitted by a plurality of BSTs and received by the MR.
According to an aspect of the invention, there is provided a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, the method comprising the steps of: characterizing the relative group delay, Δt, - Δtj5 at a RR between a signal received from the ith BST and a signal received from the jth BST prior to operation; removing the characterized relative group delay, Δt, - Δtj, from the estimated Time Difference Of Arrival, τ, - τ}, at the MR between a signal received from the ith th BST and a signal received from the j BST to yield a corrected Time Difference of Arrival; and finding the location of the MR using the corrected Time Difference of Arrival. Additional refinements of this aspect of the invention are set forth in claims 72- 74 attached hereto, which are hereby incorporated in this summary by reference. According to an aspect of the invention, there is provided a method of estimating the location of a MR in location state 0 or location state 1 from a plurality of signals transmitted from a plurality of BSTs, wherein the transfer of the positional information from the Host to a user is achieved over the internet.
Additional refinements of this aspect of the invention are set forth in claims 76- 78 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, comprising: using Least Squares to solve for the positional information of the MR and using TDOA observations from distant BSTs to resolve for ambiguities in an AOA solution; using AOA observations from different sectors at the same cellular BS to resolve for ambiguities in an AOA solution; and using range observations from different BSTs to resolve for ambiguities in an AOA solution.
Additional refinement of this aspect of the invention is set forth in claim 80 attached hereto, which are hereby incorporated in this summary by reference According to an aspect of the invention, there is provided a system for estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, the system including a host, and the system comprising; the MR being adapted to receive some or all the following information from a remote source and to apply said information to the execution of a location state: identification of base station pilot signals for the mobile to search for pilot detection threshold; duration of time for searcher to accumulate pilot energy at each delay offset; search window time length or duration; spacing of successive sampling points within the search window; frequency at which to conduct the search for pilot signals; cellular system which the mobile should acquire; the number of times the mobile should repeat all or part of the location procedure; correlation interval; locations of base stations; and the MR being adapted to transmit one or more of the following information to the Host: times of arrival of pilot signals; identification of the pilot signals reported by said device; received signal qualities of pilot signals; pilot search results; and carrier frequencies of reported pilot signals. According to an aspect of the invention, there is provided a system for estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, the system comprising: first computing means for characterizing the relative group delay, Δt, - Δtj5 at a RR between a signal received from the ith BST and a signal received from the jth BST prior to operation; and second computing means for removing the characterized relative group delay, Δt, - Δtj, from the the estimated Time Difference Of Arrival, τ, - τ,, at the MR between a signal received from the ith BST and a signal received from the jth BST. According to an aspect of the invention, there is provided a system for estimating the location of a MR in location state 0 or in location state 1 from a plurality of signals transmitted from a plurality of BSTs, wherein the transfer of the positional information from the Host to a user is achieved over the internet. Additional refinement of this aspect of the invention is set forth in claim 84 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a system for estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where: Least Squares is used to solve for the positional information of the MR;
TDOA observations from distant BSTs are used to resolve for ambiguities in an AOA solution;
AOA observations from different sectors at the same cellular BS are used to resolve for ambiguities in an AOA solution; and AOA observations from different sectors at the same cellular BS are used to resolve for ambiguities in an AOA solution.
Additional refinement of this aspect of the invention is set forth in claim 86 attached hereto, which are hereby incorporated in this summary by reference. According to an aspect of the invention, there is provided a method to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs; the MR receiving some or all the following information from a remote source and applies said information to the execution of a location state: identification of base station pilot signals for the mobile to search for pilot detection threshold; duration of time for searcher to accumulate pilot energy at each delay offset; search window time length or duration; spacing of successive sampling points within the search window; frequency at which to conduct the search for pilot signals; cellular system which the mobile should acquire; the number of times the mobile should repeat all or part of the location procedure; correlation interval; locations of base stations; and wherein some or all of the following information is transmitted from the MR to the Host: times of arrival of pilot signals; identification of the pilot signals reported by said device; received signal qualities of pilot signals; pilot search results; carrier frequencies of reported pilot signals.
According to an aspect of the invention, there is provided a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where misclosures and/or standardized residuals are used to flag observations that might contain a blunder.
Additional refinements of this aspect of the invention are set forth in claims 89-97, 106-112 and 114 attached hereto, which are hereby incorporated in this summary by reference. According to an aspect of the invention, there is provided a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of AOA and TDOA observations are used such that the
MR can be located with as few as two AOA independent observations and one TDOA observation.
Additional refinements of this aspect of the invention are set forth in claims 99 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, wherein a combination of AOA and TOA observations are used such that the MR can be located with as little as one AOA observation and one TOA observation.
Additional refinements of this aspect of the invention are set forth in claims 101- 102 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of TOA observations are used such that the MR can be located with as little as two independent TOA observations.
Additional refinements of this aspect of the invention are set forth in claims 104- 105 attached hereto, which are hereby incorporated in this summary by reference According to an aspect of the invention, there is provided a method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, which quantifies the effect of geometry on speed and direction of travel estimation using the design matrix in equation (51).
According to an aspect of the invention, there is provided a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where misclosures and/or standardized residuals are used to flag observations that might contain a blunder. Additional refinements of this aspect of the invention are set forth in claims 1 16-118 and 122-125 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of AOA and TDOA observations are used such that the MR can be located with as few as two AOA independent observations and one TDOA observation.
According to an aspect of the invention, there is provided a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of AOA and TOA observations are used such that the MR can be located with as little as one AOA observation and one TOA observation.
According to an aspect of the invention, there is provided a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of TOA observations are used such that the MR can be located with as little as two independent TOA observations.
Additional refinements of this aspect of the invention are set forth in claims 122- 123 attached hereto, which are hereby incorporated in this summary by reference.
According to an aspect of the invention, there is provided a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, which quantifies the effect of geometry on position estimation, with both AOA and TDOA observations, using the design matrix in equation (46).
According to an aspect of the invention, there is provided a system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, which quantifies the effect of geometry on speed and direction of travel estimation using the design matrix in equation (51).
Description of the Figures:
Preferred embodiments of the invention will be described with reference to the Figures, by way of example, and without intending to limit the invention to the specific embodiments disclosed, in which Figures like reference characters denote like elements, and in which:
Figure 1 illustrates the broadcast of a signal, s,(t) (103), by the ith BST (102) where τ0,, is the Time of Transmission (TOT) of the signal. A MR (101) receives the signal, together with other signals transmitted from a number of BSTs.
When the BST (102) corresponds to an IS-95-based sector in a BS, the signal, s,(t) (103), can be a pilot signal, a sync signal, a paging signal or a Forward Traffic signal. When s,(t) is a pilot signal, it consists of a single code, repeated exactly 75 times every 2 seconds, with a period of 32,768 chips. The pilot code is known as the "pilot PN sequence" or the "short code". Although in IS-95 all sectors broadcast the same short code, the pilot transmitted by each sector is time-offset by an integral multiple of 64 chips from any other pilot signal which could forseeably be simultaneously received by a MR. As the code has a period of 32,768 chips, and pilots are spaced by 64 chips, there are 512 unique pilot phases possible. When the BST (102) corresponds to a GPS-based satellite, the signal, s,(t)
(103), can be either a Coarse Acquisition (CA) signal or a Precision signal. When s,(t) is a CA signal, it consists of a unique Gold code per satellite (i.e. up to 24 distinct codes), repeated exactly 1000 times every 1 second, with a period of 1023 chips, with BPSK overhead messages at 50 bits/second. The codes are known as the "CA-codes". When the BST (102) corresponds to an AMPs-based sector or an IS-136/GSM- based sector in a BS, the signal, s,(t) (103), can be a Forward Control Channel (FOCC) signal or a Forward Voice Channel (FVC) signal. When s,(t) is a FOCC signal, it consists of Digital Colour Codes (DCCs) with overhead messages regarding the channel numbers of the paging channels. Each DCC consists of known data, unique to its corresponding servicing sector, interspersed with overhead messages.
Following Figure 1 , Figure 2 illustrates the broadcast of the signal s,(t) by the il BST (201) at Time 'τ0,i' and its reception by the MR (202) at Time of Arrival: 'τj', for i=l, 2, 3.
When the BST is IS-95-based or GPS-based, the time τ0,i is the time of transmission of a given pilot PN sequence while τi is the time of reception of such a PN sequence. When the BST is AMPs-based, the time τ0,i is the time of transmission of a given DCC sequence while τj is the time of reception of such a sequence. In order to solve for the horizontal coordinates (x,y) of the MR (202), the minimum number of BSTs required depends on the positioning technology used. Figure 3 illustrates Angle of Arrival for several wavefronts from BSTs at a MR, and assumes that the MR is far from BST, and BST, with respect to the baseline between the two BSTs. Such an assumption implies that the transmitted wavefronts are planar, i.e.
Figure imgf000025_0001
γy, where: • γs is the Angle Of Arrival (AOA), of signal s,(t) from the ith BST to the MR,
• γj is the AOA of signal s,(t) from the jth BST to the MR, and
• γi,j is defined as the angle formed between:
1. the line joining the ith BST to the jth BST (which we refer to as liney); and, 2. the line joining the MR (301) and the middle of liney between the two
BSTs (which we refer to as line0,ij) in a clockwise manner from line0,i,j to liney.
Figure 4 illustrates the Direction of Travel (DOT) 'φ' (401) relative to Northing (in a clockwise manner from Northing) and the speed v (402) of the MR (404) relative to the BSTs (403). 'φ' and v represent the velocity v of the MR with respect to the BSTs (403). In a cellular system, the BSTs are stationary which is not the case in a satellite- based system.
Without loss of generality, Figure 5 illustrates the two-dimensional (horizontal) Locus of Position (501) for TDOA2;l which is defined as
TDOA2>1 = (τ2 - τ0j2) - (τ, - τ0ιl )
= - V(χ2 - x)2 + (y2 - y)2 - - (χ, - χ)2 + (y, - y)2 c c
(1)
where c is the speed of propagation, (x, , yj , z, ) are the coordinates of the antenna of BSTι , (x2, y2 , z2) are the coordinates of the antenna at BST2, and (x,y) are the coordinates of the MR (502). This is achieved using TDOA-based Hyperbolic Multilateration (as shown in Turin, G. L. et al., "A Statistical Model of Urban Multipath Propagation," IEEE Transactions on Vehicular Technology, Vol. VT-21 , No. 1, February 1972, and as shown in Smith, J.O. et al., "Closed-Form Least-Squares Source Location Estimation from Range-Difference Measurements," IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-35, No. 12, December, 1987).
In the case where the three-dimensional coordinates of the MR are required, we have to solve for (x,y,z) using
TDOA2 1 = (τ2 - τ0ι2) - (τ1 - τ0ιl )
= - V(χ 2 - χ)2 + (y 2 - y)2 + (z2 - z)2 c
- -^ - x^ + Cy. - y^ + Cz. - z)2 c (2)
Without loss of generality, Figure 6 illustrates the two intersecting locii of Position for the two-dimensional coordinates (x,y) of the MR (603) based on TDOA2,l (601) and TDOA3 (602). In other words, it is possible to solve for (x,y) as the intersection between the two trajectories obtained using three Times of Arrival (TOAs) (after choosing the correct side of each hyperbola). In order to solve for (x,y,z) we require four TOAs. Figure 7a illustrates the Second Stage of the WLS where the signal, s,(t), is received by the MR (701). The MR (701) processes the received signal, r,(t), to estimate its own coordinates (x,y), speed v and DOT φ with the assistance of a Reference Receiver (RR) (702) of known location. The function of the RR is to estimate the Reference Information regarding the BSTs and to transmit it to the MR (701). We will refer to such a system as Option I.
Figure 7b illustrates an alternative of the Second Stage of the Wireless location system where the signal, s,(t), is received by the MR (703). The MR processes the received signal, r,(t), and transmits the Location information to the RR (704). In turn, the RR estimates the coordinates (x,y) of the MR, its speed v and DOT φ, based on the Location Information transmitted by the MR (703) and based on the Reference Information estimated by the RR (704). We will refer to such a system as Option II.
Figure 7c illustrates yet another alternative of the Second Stage of the Wireless location system where the signal, s,(t), is received by the MR (705). The MR processes the received signal, r,(t), and transmits the Location Information to a Third Receiver (TR) (706). In turn, the RR (707) estimates the Reference Information regarding the BSTs and transmits it to the same Third Receiver (TR) (706). The TR estimates the coordinates (x,y) of the MR, its speed v and DOT φ, based on the Location Information transmitted by the MR (705) and based on the Reference Information transmitted by the RR (707). We will refer to such a system as Option III.
Figure 8 illustrates the description of Design I for the MR. The received RF signal r,(t) transmitted from the ith BST is initially filtered by an RF Band Pass Filter (BPF) (802), amplified by an RF amplifier (803), down-converted by one or more Intermediate Frequency (IF) stages (810). Each IF stage consists of a mixer (804) with an LO (808) generated using an RF synthesizer (809), an IF BPF (805) and an IF amplifier (806). After the IF stages, the resulting analog IF signal (807) of IF frequency fi is filtered by a pre-aliasing filter (811), sampled by an A/D (812) and filtered again using digital filters (813) in order to generate a digital IF signal (814) of center frequency f2, where f2 < fi . An optional Direct Digital Converter (DDC) (815) is used to down-convert the digital IF signal (814) to Baseband. The resulting signal is then processed by a Digital Signal Processor (DSP) (817) in order to estimate the Frequency Of Arrival (FOA) of the received signal, r,(t). A function of the estimated FOA is fed- back to a Direct Digital Synthesizer (DDS) (819) which controls the RF synthesizer (809). The DSP (817) also provides a reference signal (820) to the DDS (819). A preferred source for the reference signal (820) is one that is common to all BSTs such as a GPS signal.
Figures 9a, 9b, 9c and 9d illustrate a flow diagram for the positioning strategy. The positioning strategy consists of a number of algorithms (Least Squares, Location On the Conic Axis (LOCA), Plane Intersection) and numerous decisions. Given a set of observations, there are 4 possible outcomes: two LS position solutions (908), two LOCA position solutions (909), one LS position solution (923, 936), and no position solution (927).
Figure 10 presents the preferred embodiment of an algorithm by which a MR can obtain TOA estimates of multiple IS-95 pilot signals. A search is made for each pilot over multiple time lags. If a pilot is detected, then it is further observed to obtain a TOA estimate of refined accuracy.
Figure 11 shows correlation values which might be obtained when searching for a particular pilot over multiple time lags, or delay offsets. Ec/Io is the ratio of average pilot chip energy to the total received power. A higher value of Ec/Io indicates greater likelihood of arrival of the pilot signal at a particular delay offset. In this figure, the delay offsets are expressed in IS-95 chips (approximately 0.814 μs). In this figure, the horizontal dotted line indicates a threshold (-15 dB in this case) for detection of the pilot arrival. An Ec/Io value greater than this threshold is considered to indicate an arrival of the pilot signal at the MR, at that particular time lag.
Figure 12 indicates a possible architecture for the embodiment of the MR as a CDMA CT, which may be used to implement the algorithm of Figure 10. Under the direction of the Controller, the Searcher searches for pilot signals at different time lags of the received signal of the MR. The Demodulator Fingers are assigned to different pilot signals. When assigned to a pilot, a Demodulator Finger will slew its timing to track the pilot, and may demodulate data transmitted synchronously with that pilot. This architecture is presented in U.S. Patent number 5,764,687 Mobile Demodulator Architecture for a Spread Spectrum Multiple Access Communications System.
DETAILED DESCRIPTION OF PREFERRED EMBODI ENTS_ THE INVENTION
Theory:
In a WLS, many factors affect the system performance:
1. RF shadowing and flat fading,
2. frequency offsets (including LOs drift and Doppler Shifts),
3. clock errors,
4. time delays,
5. noise,
6. multipath (selective fading),
7. interference;
8. geographical geometry of the BSTs relative to the intended MR.
Each factor degrades the estimated location of the MR depending on the technology employed for extracting the independent equations required for location.
RF Transmission: More specifically, the Low Pass (LP) equivalent signal, (t) , transmitted from the iΛ BST can be modeled as
s (t) = e-j2π(f<+Λfi)τ"- ej(2πΔf-t+ζi p,(t - τ0>1.) (3)
and the RF transmitted signal, s;(t), (see Figure 1) can be expressed as
Figure imgf000029_0001
where
• Re { . } denotes a real part operation;
• ξ. is the carrier frequency,
• Δf; is the frequency offset from fc (usually unknown to the MR since it depends on the BST; it is to be estimated by the RR),
• τ0 j = Time Of Transmission of signal p,(t) (also unknown to the MR),
• ζ; is the phase of the frequency fc + Δf, (also unknown to the MR),
• j = TT, and
• p,(t) is a band-limited baseband signal transmitted by the i,h BST. RF Reception: Then, the LP equivalent signal, r^(t) , received by the MR and corresponding to sj (t) can be modeled as
ξ(t) = Gj R ej(θ| + ) ej2π(δfi+Δfl) t PiCt - Tj + τo i - Δts) + Wi(t) + MPi(t) + ^(t)
(5)
and the RF received signal, r{(t), (see Figure 2) can be expressed as ri(t) = Re{ej2πf^(t)} (6) where
. Re { . } denotes a real part operation;
. G; represents the gain (real) due to the transmitting and receiving antennas: from the 1th BST to the MR antenna (i.e. a function of the two antenna patterns),
. represents the attenuation (real) due to the propagation channel, where
. n is a real number (usually between 2 and 4) that depends on the channel and . R^ is the range between the iώ BST and the MR,
. θj = {-2π( + δζ +Δfi ) (xj - τ0ιi + Δt, ) }mod 2π (7)
is the phase of the received RF signal at the MR, where:
. δζ is the frequency offset due to the Doppler shift over the propagation channel (a function of speed, v, with respect to the BSTs and Direction Of Travel (DOT), φ, with respect to Northing). . X; = TOA of a signal p;(t-x0 f ) (a function of the range R,), . Δtj is the overall group delay through the im BST to the MR antenna (usually a function of frequency), . w;(t) represents the Additive White Gaussian Noise (AWGN, complex) due to thermal noise at the MR (a function of temperature and bandwidth), . MPj(t) represents all multipath components (complex) at MR antenna (a function of the environment and of the elevation of the antennas), and 00/00224
29
. L (t) represents low-pass equivalent interference (both in-band and out-of- band) at the MR.
Total RF Reception: In IS-95, a number of pilot signals, sync signals, paging signals, and forward traffic signals, are transmitted over the same frequency band in the forward link. In GPS, a number of signals are transmitted from the satellites over the same frequency band. In other words, in a CDMA system I(t) represents all CDMA signals excluding ι (t).
Baseband Reception: Several IF stages down-convert the received RF signal, η(t), to an analog baseband signal, f;(t), which is equivalent to performing the following operation: t) = Re{ej2πf^ ιf(t) < «.+"•> e*- }= Rcføt) e-j2πΔf°' e*( (8)
where
• Δf0 is the frequency error between the Local Oscillators (LOs) in the MR receiver and the carrier frequency f,.; and
• ζo is the carrier phase in the Local Oscillators (LOs) of the MR.
From equations (5), (6) and (8), one can refer to:
• "ψ. = θ, - ζ, + j" as the Phase Of Arrival (POA) of the received signal, (t) , from the iΛ BST at the MR; and
• "f; = δfj - Δ + Δf " as the Frequency Of Arrival (FOA) of the received signal, ff (t) , from the im BST at the MR.
Digital Reception: Finally, an Analog-to-Digital (A-to-D) Converter converts the baseband signal, ?(t), to a digital signal, r(m) , which is equivalent to performing the following operation:
?( ) = f(t)|t=mΔT + q(m) (9)
where
• ΔT is the sampling interval in the MR;
• m is an integer; and • q(m) is the quantization noise which depends on the number of bits in the A-to-D.
Sources of Location Information in equation (5):
There are several sources of information where one can extract an independent equation from equation (5) regarding the static location of the MR, i.e. regarding the coordinates, (x,y), of the MR:
1. Rp,
Figure imgf000032_0001
3. Pi(t - τi+ x0> Δt .
In addition, there is one source of information where one can extract an independent equation from equation (5) regarding the kinematic location of the MR, i.e. regarding the velocity, v, namely:
4. ξ .
Starting with the static information, we have: _Rp: In the static case, the first source, Rp, can be used to estimate R;.
However, Rp is unreliable due to the nature of the radio channel where RF shadowing and flat fading can deviate the value of "n" in "Rp" from 2 (for Line Of Sight (LOS)) to 4 (for an urban environment, as shown in Hata, M., "Empirical Formula for Radio Propagation Loss in Land Mobile Radio Services," IEEE Transactions on Vehicular Technology, Vol. VT-29, No. 3, August 1980), or even 6 in heavy urban environments. Hence, unless LOS is guaranteed between the MR and the 1th BST, Rp can only offer an accuracy of a few kilometers for Rj. In the case when LOS is guaranteed (e.g. in a flat environment such as on the sea or over a lake), G; has to be known to the im BST and "n" can be chosen to be 2.
2. ψ;-. The second source, the POA
Ψi = { -2π( fc+ δf; +Δξ )(x; - τ0 i + Δt, ) + Δζ }mod 2π (10) can be used to estimate the TOA, τh where Δζ; = ζ, - o- The POA, ψ;, is related to the range, R,, between the iA BST, and the MR through equations (10) and (11):
Ri = ( τ0ιi) c (11)
where c is the speed of propagation.
Since the time of transmission, xo i , in (11) is unknown to the MR, one can either:
1. estimate it, then use the estimated value in TOA positioning, or
2. remove it using TDOA positioning.
Both TOA and TDOA positioning are explained later in the patent.
In TOA positioning, a possible method to estimate x0 is to estimate the Round Trip Delay (RTD) between:
• the transmission from the MR to an active BST and
• the reception by the MR of the response of the active BST. In this case,
. 2Ri = (2xi - xMR) c (12)
where ^ is the time of transmission from the MR to the BST and 2 ; is the TOA of the response of the BST to the MR. A drawback in estimating the RTD from equation (12) is that equation (12) assumes that the BST is able to transmit a response back to the MR. This assumption makes the handset-based WLS depend on the functionality of the BSTs. This implies that the handset- based WLS is network dependent. Nonetheless, we will carry on using such a concept in order to be able to aid a handset-based solution with a network- based solution.
Furthermore, in TOA positioning, a drawback in using θ; in equation (10) to estimate x; is the existence of a "modulo 2π" operation in (10) which implies that there can be a large integer ambiguity in the TOA estimation.
In TDOA positioning, a pair of TOAs, τ, and j, is needed per independent equation. In this case, the amount of ambiguity is a function of the baseline between the two transmitting antennas, i.e. between the im BST and the jϋl BST with respect to the wavelength corresponding to the frequency: fc + δf; +Δf, in
(10), as shown in equation (13):
τ, - Xj = (x0>i - x0J) - {(ψ; - ψj) - (Δζ; -ΔQ) + 2π n}/2πfc - (Δt, - ΔL) + extra term
(13) where
• X; corresponds to the TOA of (t) at the MR;
• Xj corresponds to the TOA of r,-(t) at the MR;
• x0 j corresponds to the TOT of p;(t- xoi);
• x0J- corresponds to the TOT of p;(t- x0j);
• n is the integer ambiguity; and
• the extra term in (13) is due to the existence of
• Δζ in (10) which depends on the error in the LOs for the im BS relative to £; and to
• δζ in (10) which depends on the velocity v of the MR and the wavelength λ corresponding to fj..
The term (xo i - x0j) is known to the MR and can be included in equation (13).
In the IS-95 and J-STD-008 standards, the CT frequency, ξ. + Δξ,, can deviate by up to 0.05part per million (i.e. by up to 95Hz for a 1900MHz carrier frequency). Furthermore, in cellular telephony in North America, the carrier frequency, fc, can take values around 800MHz (with a wavelength λ of around 37.5cm) or values around 1.9GHz (with a wavelength λ of around 15.8cm). This implies that for a velocity of lOOkm/hr δf corresponds to:
-75Hz > δfj >75Hz for ζ = 800Mhz, or to
-150Hz > δζ >150Hz for ξ = 1.9GHz.
In conclusion:
1. fc » δf, and fc » Δf0, i.e. the extra term in (13) is negligible; and
2. the amount of ambiguity in (13) can be potentially large unless the baseline between the two receiving antennas is small with respect to the wavelength, λ, corresponding to the frequency: fc. We refer to such a solution as Solution I. Solution I: When the baseline, d, , between the antenna at the i"1 BST and the antenna at the j"1 BST is of the same order of magnitude as the wavelength, λ, it is possible to generate a solution using either TDOA or PDOA. The TDOA solution between the two antennas is a hyperbola while the PDOA solution between the two antennas is a line. Both solutions approximate the exact solution and coincide asymptotically, i.e. the PDOA line and the TDOA hyperbola coincide at infinity.
Remark: From a practical point of view, the baseline, cLj, can be of the same order of magnitude as λ only if the two antennas belong to different sectors of the same cellular base station.
We assume that the MR is far from both antennas with respect to d^. Such an assumption implies that:
1. the transmitted wavefront is planar, i.e.
Figure imgf000035_0001
where:
• YJ is the Angle Of Arrival (AOA) of r£(t) from the antenna of the iώ BST to the MR,
• Yj is the AOA of r-(t) from the antenna of the jm BST to the MR, and
• YJ j is defined as the angle formed between:
• the line joining the antenna of the iΛ BST to the antenna of the jt βg-p (which we refer to as liney); and,
• the line joining the MR to the middle of lineu between the two antennas (which we refer to as lineoAJ) in a clockwise manner from line0ty to liney.
2. the following equation relates the Phase Difference of Arrival (PDOA), (ψi - ψj), to the angle y^:
— 2π^dϋ cos(γy ) = (ψj - ψj )mod 2π (14) where:
• ψ; is the Phase Of Arrival (POA) of (t) at the MR; and
• ψj is the POA of (t) at the MR.
The solution for γ;j in (14) is T A00/00224
34 λ γifj = ±cos~'<! (Ψi-ψj + 2πk) (15)
2πd >.J
where k in (15) is an integer that has to satisfy the following condition:
Figure imgf000036_0001
For example, when (Lj = λ, we have the following possible solutions for k:
when
Figure imgf000036_0002
= -1, we have k = 0,1 or 2;
2π
(ψj -ψ-)mod2π when = + 1, we have k = 0,-lor-2;
2π
(ψj-ψ:)mod2π when - = -1/2, we have k = 0orl; (17)
2π
(ψj -ψ-)mod2π when = + 1/2, wehave k = 0or-l;
2π
(ψj-ψi)mod2π when = 0, wehave k = 0,lor-l;
2π
Some of the solutions in (17) are trivial. As d becomes larger than λ, we have more nontrivial solutions. Even though having more than one solution implies ambiguity, it is possible to resolve the ambiguity using Least-Squares Techniques as explained later.
In processing ψ; to estimate either x; or YJ, we have the following sources of errors:
• The effect of noise, w,(t), on γt:
When using a linear array of antennas to estimate the AOA of an RF signal approximated as a planar wavefront (i.e. assuming a distant RF source), the Cramer-Rao Lower Bound on the variance, van^), of the estimated AOA, γ,, is equal to 0224
35
12c" vari (γ,)| Tone (18)
SNR 1 ; iTone x 4π2 x M x (M2 - 1) x d2 x sin2 γ, x f m 2 where
• var(γ,)| is the variance of the estimate γ, of the bearing
• SNRj e = Af /(2σ2),
• Ai is the amphtude of the tone,
• σn 2 is the noise variance,
• M is the number of elements in the antenna array,
• d is the distance between antenna elements,
• fn, is the frequency of the tone, and
• c is the speed of light.
Equation (18) is based on the transmission of a narrowband tone from each element of the array at the same phase. When the transmitted signal Sj(t) is either a wideband signal (e.g. CDMA) or s;(t) has a unique frequency f^ (e.g. FDMA), (18) no longer applies directly, but can be made to apply after some mathematical adjustments to the received signals using the Reference Information estimated by the RR.
The effect of noise, wt(t> on ψ,: The Cramer-Rao Lower Bound provides a lower bound on the variance, var(ψj)Lon , of the estimate,
Ψi , of the phase ψjfor a tone of frequency ^ due to AWGN :
Figure imgf000037_0001
where
• vanRj ) is the variance of the estimate, R { , of the range R;,
ITone
• λ is the wavelength of the tone,
• fjn is the frequency of the tone corresponding to λ (i.e.
• SNRi[rone = A? /(2σJ) where SNT^ is the Signal-to-Noise Ratio at the MR from the iΛ BST,
• A is me amplitude of the tone at the MR from the 1th BST,
• o2„ is the noise variance, • N, is the number of samples (which is directly related to the observation interval T, through T=N/fs where fs is the sampling frequency),
• c is the speed of propagation.
Once again equation (19) is based on the transmission of a narrowband tone from each element of the array at the same phase. When the transmitted signal s;(t) is either a wideband signal (e.g. CDMA) or s{(t) has a unique frequency ^(e.g. FDMA), (19) no longer applies directly, but can be made to apply after some mathematical adjustments to the received signals using the Reference Information estimated by the RR.
• Phase offsets: The phase of the tone f^ is shifted by Δζ;. When Δ is a function of the frequency fjj,, it has to be estimated otherwise the phase ψj is distorted.
• time delays, Δt;: The propagation delay x; is affected by Δtj which represents the system delay through the antenna, cables, filters, amplifiers, etc. Once again if Δtj is a function of the frequency Jζ,,, it has to be estimated otherwise the propagation delay can vary significantly from one received signal to another.
• multipath, MP,(f): In TDOA, the multipath MPj(t) is equivalent to extra delay over the propagation channel and has to be either estimated and removed, or mitigated. In the case where the phase estimate is used to solve for the AOA of the received signal r;(t), the effect of multipath is to shift the AOA of r£(t) by some amount depending on the AOA of MPj(t) and its magnitude.
• interference, I(t) : Depending on the level of interference, 1(f) can have a drastic effect on the accuracy of the estimated phase and may saturate the RF front end of the receiver. Its effects can be mitigated with both analog and digital hardware, as well as adequate software, whether it is in-band or out-of-band.
The practical description describes methods and apparatus to estimate θ, and to mitigate its sources of errors.
p,(t - x, + xo l - Δtj): The last (and most common) source of information one can use to extract an independent equation for the static location of the MR is P T/CA00/00224
37 P,(t - τ1 + τ0ιl - ΔL) (20)
which can be used to estimate either x, or γ,:
• When ; is estimated directly we refer to the solution as Solution II. In this case, either TOA positioning or TDOA positioning can be used.
• When Yj is estimated from x;, we refer to the solution as Solution HI. In this case, AOA positioning is used.
• When Xj is estimated indirectly based on the phases of tones in Pi(t - Xj + xo - Δtj), we refer to the solution as Solution IV. Once again, either TOA positioning or TDOA positioning can be used as discussed above.
Solution TJ: The most common method for estimating x; is to cross-correlate Pi(t - tj + τ0 j - Δtj) with Pj(t), i.e. to cross-correlate r(t) with p£(t). In general several peaks result from such a correlation particularly when receiving pilot signals in CDMA where all transmitted pilots are identical except for a relative delay between pilots. In this case, the peak corresponding to x; in p,(t - Xj + x0 j - Δtj) is located within a reasonable propagation delay, Xj - xo i, with respect to a clock derived from the signal, and can be distinguished from other pilots.
• In TOA positioning, the time of transmission, xo i , from the antenna of the iώ BST to the MR has to be estimated in order to estimate x; - x0 ; .
A possible method is to estimate the RTD between the MR transmission and the reception at the MR of the BST response to the MR transmission as shown in equation (12).
• In TDOA positioning, the time of transmission, x0 ; , is estimated using an extra independent equation. This is accomplished by estimating Xj. Once again, the most common method for estimating Xj is to cross- correlate pj(t - j + x0J- - Δtj) with Pj(t), i.e. to cross-correlate r,(t) with Pj(f). In CDMA, the resulting correlation peak corresponding to τ} in Pj(t - Xj + x0j - Δtj) is located within a reasonable propagation delay, τ} - xoJ, with respect to a clock derived from the signal. The TDOA, x, - ιs, can be used to estimate the relative delay, (x, - x0 ,) - (xJ - xOJ) since the relative TOT, x01- xoj, between the im pilot and the j"1 pilot is known at the MR. Solution III: In order to estimate γ, from the TDOA solution, x; - x we assume that the MR is far from both the antenna at the iUl BST and the antenna at the j"1 BST with respect to their baseline, cLj. As mentioned previously, such an assumption practically implies that:
1. Both antennas exist at the same BST.
2. The received wavefront is planar, i.e. YJ≤YJΞYJJ, where
• Yj is the Angle Of Arrival (AOA) of η(t) from the antenna of the i^ BST to the MR,
• Yj is the AOA of rj(t) from the antenna of the j* BST to the MR, and
• Yjj is defined as the angle formed between:
• the line joining the antenna of the iΛ BST to the antenna of the jm BST (which we refer to as lineu); and,
• the line joining the MR to the middle of liney between the two antennas (which we refer to as Hne0>y) in a clockwise manner from lineoAJ to liney.
3. Equation (21) relates the angle γ^- to the Time Difference of Arrival (TDOA), X; - Xj, as follows
— COS(Yj -) = Xj - X - (21) where
• djj is the distance between the antenna at the im BST and the antenna at thej^ BSTj and
• c is the speed of propagation.
The advantage in estimating the AOA, y{ j , using p;(t - τ{ + τo i - Δtj) in equation (20) from equation (22):
Figure imgf000040_0001
instead of using ψ; in (10) from equation (23): λ
Yi = ± cos"<j — — — (ψi - ψj + 2π k) (23)
2π d '.j
is the fact that there are no ambiguities in equation (22) except for the ± in the solution of Yj, while equation (23) can have a number of ambiguities depending on the value of the integer k.
The disadvantages in estimating the AOA, yi } , using p;(t - x; + x0>i - Δtj) in equation (20) instead of using ψ; in (10) are:
1. a small baseline is required between the two antennas; and
2. the noise affects the variance, var(xj)| , of x; substantially more than it affects the variance, var(ψj) , of ψj . The second disadvantage can be explained by comparing the Cramer-Rao Lower Bound for ψ; :
Figure imgf000041_0001
with the CRLB for XJ : (25)
Figure imgf000041_0002
Assuming that SNRj|Tone N( in (24) is equal to SNRj| in (25), then the difference between (24) and (25) is f^2 in (24) compared to BW2 in (25). In AMPs the BW is approximately equal to 30KHz while f^ can be equal to the carrier frequency fc which is equal to 800MHz. The ratio between the two values in dB is 44.26. In IS-95 and J-STD-008 the BW is approximately equal to 1.23MHz while f^ can be equal to the carrier frequency fc, i.e. 800MHz. The ratio between the two values in dB is 28.13dB. In both cases the difference is large. Nonetheless, pf(t - X; - Δtj) in (20) is suitable for TDOA positioning while ψ; in (10) is suitable for AOA positioning as explained here: • Using ψj in (10) to estimate the AOA, γ , we have a lower-bound of
5 x 10"5 radians2 for a SNR, SNR.I x N-, of 30dB and for a distance, d,
' 'Tone of one wavelength, i.e. 0.375m. Tins is equivalent to having a standard deviation for γ of 7.07 x lO-3 radians, which corresponds to a range error of 7.11 meters for every kilometer range between the MR and the antennas.
On the other hand, p,(t - x, + x0>1 - Δt,) in (20) is suitable for TDOA positioning since in this case the baseline between antennas does not have to be small. For example, assuming the SNR, SNR,| , in equation (25) to be 30dB for AMPs, the lower-bound on the variance, vaι R.)l =c2 vai j) , of the range R, is 2553meters2 or equivalently ' 'p. ) ' P.(
Rj has a standard deviation of 50.53meters. On the other hand, assuming the SNR, SNRj| t , in equation (25) to be lOdB for IS-95, the range Rj is
150.69meters2 or equivalently Rs has a standard deviation of 12.28meters.
Solution TV: When x; is estimated indirectiy based on the phases of tones in Pi(t - Xj + x0 j - Δtj), an integer ambiguity can still exist in the solution unless the wavelength of the selected tones is large with respect to the range between the MR and the antenna at the iώ BST. The selected tones are discussed below where a solution for the integer ambiguity is presented.
In processing p,(t - x; + x0 ; - Δtj) to estimate x; , we have the following sources of errors:
• noise, w,(t : The Cramer-Rao Lower Bound provides a lower bound on the effect of AWGN on the variance, var(x;)| (t) , of the estimate of the delay x,:
Figure imgf000042_0001
where
• var(R )| is the variance of the estimate, R, , of the range R.,
• BW is the RMS bandwidth of p,(t), . SNR,|p|(t) = 2E, / N0,
• E, is the energy of r,(t),
• N0 is the noise Power Spectral Density, and
• c is the speed of light. • The effect of Bandwidth, BW: Equation (26) demonstrates that the bandwidth of the signal plays an important role in the accuracy of the wireless location system. In IS-95, the radio frequency (RF) channels have a bandwidth of approximately 1.25 MHz which is comparable to that of systems designed primarily for location such as Global Positioning Systems (GPS) with a BW of lMHz over Standard Positioning Services (SPS) channels and ISM-based location systems with a typical bandwidth of 10MHz. In GPS, the location system may initially use a conventional sliding correlator to obtain a set of pseudo-ranges (one pseudo-range per satellite) (see e.g. Spilker, J.J., "GPS Signal Structure and Performance Characteristics," Global Positioning System, Volume I, The Institute of Navigation, Washington D.C., 1980). The pseudo-ranges are then used in multi-lateration to obtain a position fix of the GPS receiver. A typical accuracy for a commercial one point (i.e. no differential reception) GPS receiver with SPS is around 30m RMS without Selective Availability (SA). In direct proportions, an IS-95 land-based location system which initially uses a conventional sliding correlator to obtain a TOA estimate of the transmitted radio signal followed by a hyperbolic (differential) multi-lateration of all the TOA estimates should offer a location accuracy better than 30m RMS assuming no multipath. The reason is that the correlation function from which a TOA may be estimated has a resolution which is limited to that of the Fourier transform. The traditional resolution bound on Fourier-based methods is the Rayleigh resolution criterion as shown in Haykin, S., "Adaptive Filter Theory, " 2nd Edition, Prentice Hall, Englewood Cliffs, NJ, 1991. Further processing using SR and inverse SR algorithms often yields a result with higher resolution as shown by Dumont, L.R., et al., " Super-resolution of Multipath Channels in a Spread Spectrum Location System," IEE Electronic Letters, Vol. 30, No. 19, pp. 1583-1584, September 15, 1994 and as shown by Fattouche et al., U.S. Patent No. 5,570,305 issued Oct., 1996, and as shown by Ziskind, I. et al., "Maximum Likelihood Localization of Multiple Sources by Alternating Projection," IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-36, No. 10, October, 1988. It is proposed to use SR and inverse SR algorithms to improve the resolution of the WLS beyond the Rayleigh resolution, as opposed to
Dent, U.S. Patent No. 5,404,376 issued April 4, 1995; Dunn et al., U.S. Patent No. 5,600,706 issued Feb. 4, 1997; Khan et al., U.S. Patent No. 5,646,632 issued July 8, 1997; Ghosh et al., U.S. Patent No. 5,646,632 issued April 16, 1996;
Krasner, U.S. Patent No. 5,663,734 issued Sept. 2, 1997; and Krasner, U.S. Patent No. 5,781,156 issued July , 14, 1998; which do not further process the correlation peak.
• The Effect of the Observation Interval, Tt:
The Observation interval, T; , is directly related to the energy E; in the received signal (t) at the MR since
Figure imgf000044_0001
It is also related to the steady state response of a filter in such a way that the response of the filter reaches its steady state as long as
Figure imgf000044_0002
where BWj is the bandwidth of the filter that is processing (t).
• time delays. At,: The propagation delay τ{ is affected by Δtj which represents the system delay through antenna, cables, filters, amplifiers, etc. When Δtj changes from one BST to another independently from the RR, it has to be estimated otherwise the propagation delay can vary significantly from one received signal to another.
interference, 1(f) : Depending on the level of interference, Ij(f) can have a drastic effect on the accuracy of the estimated phase and may saturate the RF front end of the receiver. Its effects can be mitigated with both analog and digital hardware, as well as adequate software, whether it is in-band or out-of-band.
multipath, MP,(t): The multipath MPj(f) is equivalent to an extra delay over the propagation channel and has to be either estimated and removed, or mitigated. In the case the estimated time delay, χ(> is used to solve for the TDOA between the received signal r,(t) and the received signal r,(t), the difference in multipath, MPj(t) - MPj(f) , is the distorting factor which has to be mitigated.
The effect of Multipath, MP,(f):
When multipath is considered, the accuracy of the AMPs land-based WLS could potentially degrade even further. The cellular frequency band is between 824 and 894 MHz and the propagation characteristics at these UHF frequencies will have a significant impact on positioning by multi-lateration as shown in Parsons D., "The Mobile Radio Propagation Channel, " John Wiley & Sons, New York, 1992. That the ranges measured correspond to Line Of Sight (LOS) distances is a major assumption made when estimating position by multilateration. Although the dominant transmission mode in this band is LOS, reflections from natural and man-made objects as well as diffraction around said objects are also possibilities. Multipath and diffraction allow the cellular signal to propagate in heavily built up areas as well as indoors. However, they also cause the measured ranges to be longer than the true LOS distance which introduces error into the multilateration process. In addition, the propagation distance at UHF is relatively short. This allows frequency reuse in the cellular system but limits the number of observables in the multilateration process. For instance, in a dense urban environment with a delay spread of 3 microseconds (as shown in Hata, M., "Empirical Formula for Radio Propagation Loss in Land Mobile Radio Services," IEEE Transactions on Vehicular Technology, Vol. VT-29, No. 3, August 1980) multipath causes the location accuracy to degrade to more than 1400m RMS. Once again, the reason for this is that the correlation function from which the multipath may be estimated has a resolution which is limited to that of the Fourier transform which implies that any multipath within such a resolution is unresolvable using traditional methods. Further processing using an inverse SR algorithm often yields a result with higher multipath resolution as shown by Dumont, L.R., et al., " Super- resolution of Multipath Channels in a Spread Spectrum Location System," IEE Electronic Letters, Vol. 30, No. 19, pp. 1583-1584, September 15, 1994 and as shown by Fattouche et al., U.S. Patent No. 5,570,305 issued Oct., 1996, and as shown by Ziskind, I. et al., "Maximum Likelihood Localization of Multiple Sources by Alternating Projection," IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-36, No. 10, October, 1988. Another approach to multipath resolution is due to Morley, G.D. et al., "Improved Location Estimation with pulse-ranging in presence of shadowing and multipath excess-delay effects," Electronics Letters, Vol. 31, No. 18, pp. 1609-1610, 31st Aug. , 1995. It is proposed to use SR and inverse SR algorithms to better resolve the multipath components in the received radio signal, as opposed to
Dent, U.S. Patent No. 5,404,376 issued April 4, 1995;
Dunn et al., U.S. Patent No. 5,600,706 issued Feb. 4, 1997;
Khan et al., U.S. Patent No. 5,646,632 issued July 8, 1997; Ghosh et al., U.S. Patent No. 5,646,632 issued April 16, 1996;
Krasner, U.S. Patent No. 5,663,734 issued Sept. 2, 1997; and
Krasner, U.S. Patent No. 5,781,156 issued July , 14, 1998; which do not further process the correlation peak.
Choices of Tones for Solution IV:
There are existing tones and there are generated tones. Both types of tones carry the MR Location Information in their phases (mod 2π). Existing Tones in the periodic pilot sequence:
When the received signal, Tj(t), is periodic (e.g. the pilot signal in IS-95, the CA-code in GPS), it can be approximated using Fourier series, or equivalently, a number of tones can be used to approximate the pilot sequence:
Sj(t - Tj) = ∑^ exppπnf t - x,)] (29) n=-α>
where
• ζ,, is the fundamental frequency of the pilot sequence;
• X; is the propagation delay at the antenna of the im BST; and
• cn is the complex Fourier Coefficient of the pilot sequence. Using equation (30) below one can estimate the TOA, x„ as follows:
(t-x.) = {(phase of n* term in (29)) modulo 2π}/(2πnf (30) Generated Tones using a nonlinear operation on the pilot sequence:
When the received signal, r;(t), is cyclo-stationary, it generates a tone when a nonlinear operation is performed over it. The phase of such a tone carries location information regarding the MR. Most digital signals are cyclo-stationary and the generated tone is usually at the symbol rate (except in CDMA where the generated tone is the chipping rate).
In the case when the BSTs are AMPS TDMA/GSM-based, taking advantage of the cyclo-stationarity nature of the received signal is particularly important since in this case, it is not important to have a clean replica of the received signal. After performing the nonlinear operation on the received signal, a tone is generated of known (or estimated) frequency. This implies that the duration of observation can be extended as long as there is a signal to be received, i.e. as long as a call is in progress regarding a CT. This leads to an improved SNR for the received signal as explained in the CRLBs shown above.
The practical description introduces methods and apparatus to estimate x; and to mitigate its sources of errors.
4. fi: In the kinematic case where the MR moves with respect to the BSTs, the speed and Direction Of Travel (DOT) of the MR is of interest. The only source of information one can use to extract an independent equation for the speed and DOT ofthe MR is
Figure imgf000047_0001
which can be used to estimate the Doppler shift δfj . This is explained in Figure 4 which assumes a three antenna system: the antenna at the π BST with coordinates (xm,ym,zιn), the antenna at the kΛ BST with coordinates (x^y^,,) and the antenna of the n^ BST with coordinates (x^y^zj. In this case, the FOAs: fπ,, fk and are related
• to the Direction Of Travel (DOT), φ, relative to Northing (clockwise),
• to the relative speed of travel, v, with respect to the BSTs, and
• to the frequency offsets Δf„ as follows: f™ = £ + δfm - Δf0 + Δ = fc + v/λ cos(φ- Λ - Δf0 + Δfm (32a) fk = fc + δfk - Δf0 + Δfk = fc + v/λ cos(φ- Λk) - Δf0 + Δ k (32b) fn = fc + δf, - Δf0 + Δf; = fc + v/λ cos(φ- Λ - Δf0 + Δζ, (32c)
where
• A,,, is the clockwise angle from Northing to the line formed by (x^y^ and
(χ,y),
• A is the clockwise angle from Northing to the line formed by (x^yO and
(χ,y),
• A„ is the clockwise angle from Northing to the line formed by (x^yn) and (x,y), and
• (x,y) are the 2-D (horizontal) coordinates of the MR.
In other words, there are three equations (32a, b and c) with 6 unknowns: φ, v, Δfπ,, Δ , Δf„ and Δf^ This is a problem that can be resolved if the frequency offsets: Δ^, Δ , Δf„ are estimated by the RR and removed. Then we have two equations with two unknowns: v and φ after using Frequency Difference Of Arrival (FDOA):
fjπ - fk = v/λ cos(φ-ΛπL) - v/λ cosfø-Λ ) (33a) fk - = v/λ cos(φ-Λk) - v/λ cos(φ-ΛJ (33b)
We refer to such a solution as Solution V. Note that in Solution V the position of the MR (x,y) must be known (or estimated) prior to estimating the speed and velocity of the MR, in order to be able to know A^ Λ* and A„ in equations (33).
In processing f; to estimate δζ , we have the following sources of errors: • noise, w,(fl: The Cramer-Rao Lower Bound provides a lower bound on the variance, var(δf,)| , of the estimate of the frequency δζ of a tone due to fonc
AWGN:
var(δf,)| > ■ ; (34)
' hone SNRjL x N, χ (N,2 - l) χ 4π2
■ITone ' ' where vaτ(δf, )| is the variance of the estimate, δf, , of the frequency δf„ iTon • SNRj^ = A? /(2σ2 ) is the SNR of the tone at the MR from the antenna ofthe i^ BST,
• A is e amplitude of the tone at the MR from the antenna of the i"1 BST,
• σn 2 is the noise variance, and
• Nj is the number of samples (which is directly related to the observation interval T{).
• frequency offsets, Δf,: As mentioned above, there are three equations (32a- 32c) with 6 unknowns: φ, v, A^, Δfk, Δf„ and Δξ,. This problem can be resolved if the frequency offsets Δf, are estimated by the RR and removed.
• multipath, MP,ffl: The effect of the multipath MPj(t) in this case is to add some Doppler shift to δfj due to dynamic (non-stationary) reflectors such as cars, buses, trucks, etc.
• interference, 1(0 : In this case, Ij(t) plays the same role in distorting the estimate of the frequency δζ as in distorting the estimate of the phase θ, except that its effect can be reduced more significantly.
• Error in Estimating (x,y): In equations (33) it is assumed that the position, (x,y), of the MR is known prior to estimating its speed and DOT. This is usually not true and (x,y) needs to be estimated first. The estimation of (x,y) is imperfect, implying that it will contain errors that can affect the estimation of the speed of the MR and its DOT.
The practical description explains methods and apparatus to estimate δf-π, δfk and 6 and to mitigate sources of errors.
Positioning Technologies:
1. TO A/Range Positioning (Circular Multilateration)
In a 2-D (horizontal) location system it is possible to estimate the position of a MR from the range of the MR to at least two BSTs. The MR horizontal position estimate is simply the intersection of two horizontal circles with radii equal to the ranges and centered at the BSTs. This method of positioning may therefore be called circular multilateration. Ranges may be calculated by subtracting the known time of transmission x0 of signal s(f) from the measured TOA, Xj, of signal r;(t). There are three possible solutions: 1. When the Time of Transmission, x0, is known, we refer to the positioning technique as Range positioning.
2. When the time of transmission, x0, is unknown, and to be estimated, we refer to the positioning technique as TOA positioning.
3. When the time of transmission, x0 , is unknown, and to be eliminated using TDOA, we refer to the positioning technique as TDOA positioning (which is discussed in the following two sections).
In Range positioning, one way of deterrnining the time of transmission, x0, is to use the RTD between a MR and a BST with respect to the round trip delay from the time the MR transmits a signal to the BST to the time it receives the response of the BST to the MR transmission. In this case, a minimum of two independent equations is required to solve for the two unknowns x and y. If two ranges are estimated from TOA measurements of two different BST's, then there are two solutions for MR location at the intersections of two circles centered at the two BST's. This ambiguity can be resolved by including a third range measurement from another BST, or by (approximate) angle estimation from one or more sites by examining signal strength or phase difference of arrival.
In TOA positioning, there are now three unknowns: x, y and x0. TOAs from three BSs are required and the equation corresponding to the r* BST is
Figure imgf000050_0001
where
Xj is the Time Of Arrival of signal r;(t) at the MR from the antenna of the iώ BST, x0 is the time of transmission of signal s(t) from the iΛ
BST,
(x,y) is the best known 2-D position of the MR,
(xj. y, ) is the known 2-D position of the antenna at the i* BST where i=l, ..., N. 0224
49 When more than the minimum number of BSTs are available, as explained above, redundancy is said to exist. If a redundant set of measured TOAs contains errors, the TOAs must be adjusted in order to obtain a unique solution to the unknowns. The amount by which each TOA is adjusted is called the residual. Several methods of adjusting the observed according to some optimization criteria are available including:
• General Minimum Variance Unbiased Estimation,
• Best Linear Unbiased Estimation
• Maximum Likelihood Estimation,
• Least Squares Estimation,
• Method of Moments,
• General Bayesian Estimation,
• Linear Bayesian Estimation,
• Kalman Filtering, etc.
It is preferred to adjust the observed TOAs such that the sum of the squares of the residuals is a minimum. This is Least Squares Estimation.
Effect of Geometry on TOA Positioning
In a location system, geometry (that is the relative positions of the BSTs with respect to each other and the MR to be positioned) plays an important role as shown in equation (36)
std(ρosition) = DOP x std( R; ) (36) where
• std(Rj) is the standard deviation of the range estimate
Rj obtained from equations (19) or (26), and
• DOP is Dilution Of Precision which is a measure of geometry.
In a 2-D location system, it is common to refer to the Horizontal DOP (HDOP) which is defined as
HDOP = Λ/ED0P2 + ND0P2 (37)
where
• EDOP is the East DOP and is defined as the square root of the element in the 1st row and lsl column of CτOA, 0224
50 NDOP is the North DOP and is defined as the square root of the element in the 2nd row and 2nd column of -TOA'
CT0A is an unsealed matrix defined as
Figure imgf000052_0001
C, is the unsealed measurement co variance matrix (the identity matrix of appropriate dimension), A is the design matrix for the model of equation (35) and is defined as x-x, y-yi -c d. d, x-x2 y- y2 -c
A = (38) d* d2 x-xN y-yN -c
*N <*N
and dj is the best derived distance (range) between the MR and the antenna at the iΛ BST.
Weighting of the TOA Observations
In the position estimation process, not all TOA observations need carry the same weight. TOA observations thought to be more reliable may be weighted more heavily than those deemed less reliable. This is accomplished through the observation covariance matrix. The inverse of C, is the weight matrix.
Larger values on the diagonal of C,_1 correspond to heavier weighting for the corresponding TOA observations. In the context of cellular telephone positioning, RSSI from each BST is one method of assigning weights to the TOAs. A high RSSI from a BST implies a reliable TOA. This is due to two facts:
1. The RSSI from a BST usually consists of received signal power + received noise power. The noise in the receiver is mainly thermal noise which is a function of bandwidth and temperature. When two BSs have comparable temperatures and bandwidths, the received noise power is approximately the same in both BSTs. Thus, a high RSSI implies a high received signal power, which in rum implies a high SNR. This is desirable. 2. Furthermore, a higher RSSI usually implies less shadowing than a lower RSSI, which in turn implies less multipath. This is also desirable.
Blunder Detection in TOA Positioning
Blunders are gross errors in the TOA observations. They may be caused by large signal level fluctuations due to either flat fading or sudden in-band interference. If unremoved, blunders cause large errors in the estimated position. It is possible to detect observations containing blunders by observing the misclosure of each observation during the iterative Least Squares process. Misclosure is defined as the value of the position model (35) given the best available position estimate (x,y). Observations containing blunders will generally have much larger misclosures than observations not containing blunders. When an observation is detected as having a blunder it may be removed from the position estimation process.
We attempt to reduce the effect of geometry, i.e. reduce the value of HDOP in (37), by allowing the MR to monitor a large number of BSTs. In cellular communications, frequency reuse and fiat fading are common occurrences. Therefore, increasing the number of monitoring BSs (which probably reduces
HDOP) generally increases std(Rj). As a result blunder detection is crucial as part of the method of n inimizing positional error by maximizing the number of BSs (and hence reducing HDOP) without incurring a large penalty on std(R; ). In cases where blunders are not detected, poor SNR at some BSs may cause the inclusion of measurements from those BSs to increase the std(range) more than they reduce HDOP. With this trade-off in mind, the number of BSTs used to locate the MR is optimized such that stdfjposition) in (36) is minimized.
Blunders may also be detected by a statistical analysis of the observation residuals computed from the Least Squares process as shown by Vanicek, P., Krakiwsky, E., "Geodesy: The Concepts," North-Holland Publishing Company, Amsterdam, 1982. The residual of each TOA observation may be standardized by its own estimated standard deviation such that the entire set of residuals is assumed to belong to the normal distribution with zero mean and unit standard deviation. If this hypothesis is correct, the standardized residuals should fall within some specified confidence region for the standard normal distribution. An observation whose residual is flagged is suspected of containing a blunder.
2. TDOA Positioning (Hyperbolic Multilateration)
Instead of estimating the unknown time of transmission, it is possible to eliminate it. This is accomplished by differencing TOAs from two different BSTs. Since the time of transmission from both BSTs is known at both the RR and MR, it can be eliminated from the resulting TDOA (Time Difference Of Arrival). It can be shown that the locus of points for which a particular TDOA is valid corresponds to a hyperbola. The side of the hyperbola on which the MR must lie is known by the sign of the TDOA. Given TOAs from three BSs, two independent TDOAs may be formed. The intersection of the two corresponding hyperbolas estimates the position of the MR. This method is commonly referred to as hyperbolic multilateration.
The 2-D positioning model for hyperbolic multilateration is
Δxij -i (χ - i)2 + (y- yi)2 +1 /(χ -Xj)2 + (y- yj)2 = o (39)
where
• Δxij = xi - xj , i ≠ j, and
• Xj is the Time Of Arrival of signal r;(t) at the antenna of the iΛ BST.
Effect of Geometry on TDOA Positioning
Geometry affects TDOA positioning as well. The HDOP is again calculated from (37) where the design matrix is now
Figure imgf000055_0001
where N in (40) is the number of BSTs. Note that the TOA from the antenna of the first BST is subtracted from all other TOAs.
The unsealed observation covariance matrix is
2 1 •• • 1 l ■ • . l ;
C, = (41) i l • •. 1
1 • ■• 1 2
Weighting of the TDOA Observations
Weighting of the TDOA observations is possible. However, because the TDOAs are a function of two TOAs, the method of assigning weights is no longer straightforward.
Blunder Detection in TDOA Positioning
Blunder detection may also be performed in TDOA positioning. However, in this context, misclosures and residuals correspond to TDOAs. Therefore, a failing misclosure or residual may be due to a blunder in either of the TOA observations from which the TDOA is derived. It is not always possible to isolate the offending BST.
Multiple Solutions in TDOA Positioning
Two hyperbola halves, formed from two independent TDOAs, may intersect twice. This results in two mathematically correct solutions. This is particularly true in the case of positioning MRs where short distances and poor geometry are commonplace. Solution bifurcation (the existence of two solutions to the positioning equations) is most often a concern in the exactly determined case. For the exactly determined case, the existence of two solutions can be detected with the method given in Chaffee, J.W. et al., "Bifurcation of Pseudorange Equations,'' Proceedings of the 1993 National Technical Meeting, San Francisco, California, January 20-22, 1993, The Institute of Navigation. Although originally intended for the detection of bifurcation in the GPS (Global Positioning System), this method is equally applicable to the case of MR positioning.
When more than two TDOAs are available in the 2-D positioning case, the probability of exact solution bifurcation is extremely small. It is, however, possible for bifurcation to exist for a subset of two TDOAs. In such a case, the second solution may affect the final solution obtained when using all available TDOAs.
In any case, when solution bifurcation exists, the iterative Least Squares position estimation algorithm may converge to either of the solutions. The solution converged to is a function of the initial starting position used to begin the iterative Least Squares process. To converge to the solution corresponding to the actual position of the MR, the initial position used to begin Least Squares must be relatively accurate. Given no a priori information about the location of the MR, a closed-form position estimation algorithm, using the TOA or TDOA observations, is the only choice.
A number of closed-form positioning algorithms have been developed. Examples are, spherical interpolation (Smith, J.O., et al., "Closed-Form Least- Squares Source Location Estimation from Range-Difference Measurements," IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP- 35, No. 12, Dec. 1987, pp. 1661-1669), the method of Schau and Robinson (Schau, H.C., et al., "Passive source localization employing intersecting spherical surfaces from time-of-arrival differences," IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-35, No. 8, Aug. 1987, pp. 1223-1225), Bancroft's method (Bancroft, S., "An algebraic solution of the GPS equations," IEEE Transactions on Aerospace and Electronic Systems, Vol. AES-21, No. 7, Jan. 1985, pp. 56-59), the method of Chan and Ho (Chan, Y.T., et al., "A Simple and Efficient Estimator for Hyperbolic Location," IEEE
Transactions on Signal Processing, Vol. 42, No. 8, Aug. 1994, pp. 1905-1915), and LOCA (Location On the Conic Axis) by Schmidt, R.O. "A New Approach to Geometry of Range Difference Location," IEEE Transactions on Aerospace and Electronic Systems, Vol. AES-8, No. 6, Nov. 1972, pp. 821-835.
LOCA is used in the preferred embodiment of the invention. LOCA is the mathematical dual of hyperbolic trilateration. The fundamental theorem of LOCA states that TOA differences from three BSTs of known location yield a straight line of position. This straight line is the major axis of a conic. The three BSTs lie on the conic and the MR, the location of which is being estimated, lies at one of the foci of the conic. In the case of redundancy and 3- D positioning, LOCA is expanded into Plane Intersection as shown by Schmidt, R.O., "A New Approach to Geometry of Range Difference Location," IEEE Transactions on Aerospace and Electronic Systems, Vol. AES-8, No. 6, Nov. 1972, pp. 821-835.
Plane Intersection is equally applicable to the case of redundancy in 2-D positioning.
In LOCA, three conies are possible: an ellipse, hyperbola, or parabola. The conic of interest depends on the relative positions of the BSTs and the MR to be positioned. Each of these conies has two foci (for the parabola one focus is at infinity) and therefore two possible solutions. In the case, of an ellipse, the correct focus is determined from the signs of the TDOAs. In the case of the parabola, the solution is obvious. When the conic is a hyperbola, however, the correct focus cannot be determined from the given TDOAs. Chaffee's method of bifurcation detection in Chaffee, J.W. et al., "Bifurcation of Pseudorange Equations," Proceedings of the 1993 National technical Meeting, San Francisco, California, January 20-22, 1993, The Institute of Navigation, will indicate bifurcation for the case of the hyperbola and no bifurcation for the cases of the ellipse and parabola. With four BSs, a second conic is obtained. The intersection of the major axes of the two conies is the MR position estimate. Hence, there is no ambiguity and bifurcation is, therefore, not of concern when using Plane Intersection in the presence of redundancy. If the TOAs contain large errors, LOCA or Plane Intersection may give a very poor solution. When used as the initial position for Least Squares, this poor solution may cause Least Squares to diverge. In that case, the best available course of action may be to average the coordinates of the participating BSs and use that as the initial position for Least Squares. Another alternative is to use the approximate coordinates of the BST with which the MR is communicating, as the initial position.
3. Hybrid TDOA Positioning (Circular Multilateration with TDOAs)
In some cases it is desirable, or even necessary, to use TDOAs as opposed to TOAs. However, as explained above, it is not always possible to isolate blunders when using TDOAs. It is possible, however, to construct a positioning model that uses TDOAs but gives residuals and misclosures for individual TOAs.
To do this, one BST is chosen as the reference. The assumption is then made that the TOA at the reference site is equal to the time of transmission, x0. From this point on, the development follows that of circular multilateration except that the TOA of the reference BST is subtracted from all TOAs (including that of the reference BST) instead subtracting the time of transmission. The positioning model is then
Δτil +iRπ -- ('c-χi)2 +(y-yi)2 = 0 (42)
where
• ΔXj, = Xj - j , i = l,--,N and
• Ru is the range from the antenna of the reference (first) BST to the position of the MR.
The model of (42) has the same structure as that of circular multilateration (equation(35)) except that the third unknown is Ru and not the time of transmission.
The model of (42) gives N observations, and therefore N residuals and misclosures, for N BSTs. The N observations, however, are a 0 and N-l TDOAs. The structure of the equations is that of TOA positioning whereas the observations are that of TDOA positioning. Hence, the name Hybrid TDOA Positioning.
Effect of Geometry on Hybrid TDOA Positioning
The design matrix for hybrid TDOA positioning is
Figure imgf000059_0001
Note that (43) is the same as (38) except for the third column which corresponds to the third unknown.
The unweighted and unsealed observation covariance matrix, C, , is the same as that for TOA positioning, the identity matrix of dimension N.
DOPs may be calculated, as before, from [AJCJ-1 A] .
Weighting of the Hybrid TDOA Observations
Although the observations are, strictly speaking, TDOAs, they may be treated as observations corresponding to the individual BSTs. That is, the first observation, which is always zero, may be viewed as the observation of the reference (or first) BST. The second observation, consisting of the TOA of the reference BS subtracted from the TOA of the second BST, is treated as the observation of the second BST, and so on. Therefore, the N diagonal elements of C, may be weighted for the individual BSTs.
Blunder Detection in Hybrid TDOA Positioning
In like manner, N independent misclosures and residuals, one for each BST, are available for analysis. Blunder detection may then be performed as in the case of TOA positioning and blunders in individual TOA observations may be detected. Multiple Solutions in Hybrid TDOA Positioning
Since the observations used are actually TDOAs, this method of positioning is subject to solution bifurcation. The method of dealing with multiple solutions discussed under hyperbolic multilateration is equally applicable here.
4. AOA Positioning:
Given unambiguous AOA measurements, only two unknowns exist for position estimation - the 2-D coordinates of the MR to be positioned. Therefore, a minimum of two independent AOA measurements are required. Again let (x,y) be the 2-D coordinates of the MR and (xj.yj) the 2-D coordinates of the antenna of the i* BST. At the antenna of the iΛ BST, the AOA, Λ; , is measured, where:
• A; is the clockwise angle from Noithing to the line joining the MR to the antenna of the ith BST (similar to Λj in equations (32) and (33); not to be confused with γ, in equations (15), (16), (21), (22) and (23)).
The mathematical positioning model is then
(x - Xj)cosΛj -(y-yj)sinΛj = 0 (44)
When more than two independent AOA observations are available, Least Squares may be used to obtain a unique solution.
Effect of Geometry on AOA Positioning
The relative positions of BSTs with respect to each other and the MR to be positioned are also important for AOA positioning. DOP may again be used to quantify geometry. From the 2-D positioning model (44), the design matrix is found to be
cos A, -sin A, cos Λ2 - sin Λ2
A = (45)
cos ΛN - sin ΛN The unsealed and unweighted observation covariance matrix, C, , is the identity r _ j "i — i matrix of dimension N. The DOPs can again be found from [A C, Aj .
Weighting of the AOA Observations
The individual AOA observations may be weighted through the observation covariance matrix C,. Those AOA observations deemed more reliable are assigned a smaller variance, or conversely, a heavier weight.
Blunder Detection in AOA Positioning
As in the case of TOA or TDOA observations, AOA observations may contain blunders. When redundant observations are available, statistical testing of the observation residuals and misclosures may be used to detect observations which contain blunders.
Notes:
1. As mentioned above, the AOA solution in equation (15) has an ambiguity problem which can be resolved either by
• using TDOA as well as AOA as explained in the next section, or
• using more than two antennas at the same BS. This is possible in sectorized cells since each sector has usually two diversity antennas. In other words, in a three sector cell a total of six antennas could be available.
2. When the diversity antennas are separated only vertically, the estimated AOA is the elevation angle. In this case, it is possible to use both diversity antennas since they are independent from each other, i.e. they offer independent observations. Diversity observations may be combined by well- known methods including selection combining, maximal ratio combining, co-phasing combining, equal gain combining, or other methods.
5. AOA/TDOA Positioning:
Least Squares allows for the combination of different types of observations. In particular, it is possible to estimate 2-D position using both AOA observations and TDOA observations within a single Least Squares adjustment.
A combination of the two different observation types is particularly useful in situations where only two BSTs are available for estimating the 2-D position of a MR. When two BSTs are used to estimate the incoming signal AOA for these BSTs, numerous AOA ambiguities are possible as shown in equation (15). This is due both to the nature of AOA estimation with a 1-D linear array, and to the fact that the spacing between antennas can be greater than one wavelength. Therefore, with two independent AOA observations only, many position solutions are possible and there is no way of determining which of the many solutions is correct.
If, however, in addition to the two independent AOA observations, TOA is also measured at each one of two distant BSTs (maybe using p;(t - x, + x0 j - Δtj) in equation (20)), a TDOA can be calculated. This TDOA corresponds to a hyperbola side which, in the absence of error, will cross through the intersection point of two of the many bearing lines derived from the ambiguous AOAs.
When the AOA and TDOA observations include error, no two bearing lines and the measured hyperbola side will exactly intersect. For the measured TDOA and any two particular AOAs, Least Squares will give the position solution which minimizes the sum of the squares of the residuals. Observation residuals are available since there are two unknowns, the 2-D coordinates, with three independent observations (2 independent AOAs and 1 TDOA).
In order to separate the correct AOA pair from the ambiguities, each AOA combination is combined with the TDOA observation in Least Squares. That combination of AOAs which results in the smallest sum of squares of residuals is chosen as correct. The corresponding position solution is used as the MR position estimate.
When combining AOA and TDOA observations in Least Squares, both the model for AOA positioning (44) and the model for TDOA positioning (39) are used. Without loss of generality, the design matrix A for two independent AOA observations and one TDOA observation is cosΛ ! - sin Λ, cos A 2 - sin Λ2 (46) x - x,
+ y - y2 , y - y. d, d2 d,
where d; is the best derived distance between the MR and the iώ BST. Note that the first two rows correspond to the two AOA measurements whereas the third row corresponds to the TDOA measurement. The two columns correspond to the two unknowns, x and y. Additional AOA and TDOA measurements may be included by adding appropriate rows to (46). Note that the 1/c factor seen in (40) is missing in the TDOA measurement row of (46). This is done such that the units throughout A are dimensionless.
The misclosure vector, necessary for the Least Squares adjustment mechanism, merely consists of the misclosures of all observations. The misclosure vector corresponding to the design matrix of (46) is
(x - x,)cosΛ, -(y-y^sinΛ, w = (x - x2)cosΛ2 - (y- y2)sinΛ2 (47) c*Δx12-A/(χ-χ 1)2 +(y-y1)2 +V(χ-χ 2)2 +(y-y2)2J
where
• Δx12 = x, -x2. Since the AOA and TDOA observations are independent, the unsealed observation covariance matrix is the identity matrix of appropriate dimension (number of AOA observations plus the number of TDOA observations).
Effect of Geometry on AOA TDOA Positioning
DOP may again be used to quantify geometry. The DOPs (HDOP, EDOP,
NDOP) can be found from ATC^A immediately above. Y where A and C, are defined Weighting of the AOA/TDOA Observations
The individual AOA and TDOA observations may be weighted through the observation covariance matrix C, . Those AOA and or TDOA observations deemed more reliable are assigned a smaller variance, or conversely, a heavier weight.
Blunder Detection in AOA/TDOA Positioning
Both the TDOA and AOA observations may contain blunders. Statistical testing of the observation residuals and misclosures may be used to detect observations which contain blunders.
6. AOA/Range Positioning:
If the time of transmission or round-trip delay is known, the range, Rj , from the MR to the iώ BST is the observed parameter. In that case, AOA and range positioning may also be combined in order to estimate the position of the MR with three BSTs. Without loss of generality, the design matrix A for two AOA observations and one range observations is
Figure imgf000064_0001
The misclosure vector corresponding to this design matrix is
(49)
Figure imgf000064_0002
Since the AOA and range observations are independent, the unsealed observation covariance matrix is the identity matrix of appropriate dimension (number of AOA observations plus the number of range observations).
On the other hand, without loss of generality, the design matrix A for one AOA observations and two range observations is cos A, sin Λ.
A = y- y
(50) x - x. y- y2
The misclosure vector corresponding to this design matrix is
Figure imgf000065_0001
Effect of Geometry on AOA Range Positioning
DOP may again be used to quantify geometiy. The DOPs (HDOP, EDOP,
NDOP) can be found from [ATC,"1AJ where A and C, are defined immediately above.
Weighting of the AOA Range Observations
The individual AOA and range observations may be weighted through the observation covariance matrix C, . Those AOA and/or range observations deemed more reliable are assigned a smaller variance, or conversely, a heavier weight.
Blunder Detection in AOA/Range Positioning
Both the AOA and range observations may contain blunders. Statistical testing of the observation residuals and misclosures may be used to detect observations which contain blunders.
7. Speed and Direction of Travel Estimation:
In a kinematic location system, the 3-D or 2-D location of the moving MR must be estimated at various epochs of time. In addition, the doppler shift of the signal transmitted from each BST may be estimated as discussed earlier. The equations relating the frequency (including Doppler shift) of the arriving signal at the MR from the BST
1. to the relative speed of the MR with respect to the BSTs,
2. to its DOT and
3. to the frequency offset inherent in each BST are given in equations (35).
The estimation model for speed and DOT is, therefore,
fj - v/λ cos(φ - Λj ) - Δf j + Δf0 = 0 (52)
where Δf( , the frequency offset for the iώ BST, is assumed to be known by the RR. Given that at any particular epoch the position of the MR (x,y) is estimated by one of the above methods, the angle A; in (52) may be calculated for each BST. In equation (52) then, the knowns are fj , λ, Δf{ and Λj , (where f i is measured or observed), and v, φ and Δf0 are the unknowns. Three BSTs are required - the same number required to estimate the 2D MR position using TOA or TDOA positioning.
Effect of Geometry on Speed and Direction of Travel Estimation:
Geometry affects the estimation of speed and DOT as well. For instance, it is intuitively obvious that when the MR is traveling on a line perpendicular to the line connecting it and a BST, no information regarding the speed of the MR is available from observed data at that BST.
The design matrix for the model of (52) is
Figure imgf000066_0001
The DOPs can again be found from ATC,-IAJ . In this case the DOPs will be speed DOP, direction of travel DOP, and frequency offset DOP. C, is the unsealed and unweighted identity matrix of dimension N.
Practical Description of the Invention:
Forcing a CT to act as an MR:
In a handset-based WLS, it is desirable to be able to use a cellular BS as a BST and to request a CT to act as an MR in order to fulfill the FCC mandate for E911 as well as to offer location services to the current 60 million cellular subscribers in North America who are expected to grow to 120 million by the year 2001. This however might require a change to the CT as well as a change in the cellular standards. If the required changes preclude existing CTs, the FCC mandate for E911 will not be fulfilled using such a system. Nonetheless, a handset-based solution is attractive since it can offer an accuracy comparable to the one offered by systems dedicated for location such as GPS and Loran-C. This is why it is still acceptable to make changes to the CT regardless of the FCC mandate.
Some changes to the CT might require changes to the standard* We emphasize CT changes while attempting to reduce the changes to the standard. Three levels of CT modification are identified as follows:
1. modification 0 which requires no modifications to the CT;
2. modification I which requires only software modifications to the existing CT; and
3. modification H which requires both software and hardware modifications to the existing CT.
A hardware modification to the existing CT is required when a required change in the signal processing algorithms cannot be accomplished in software. This may be
• due to required changes (e.g. a change in the RF front end, in the antennas or in the mixers) which cannot be performed in software; or
• due to the fact that the digital hardware is not adequate to perform some tasks in software (e.g. when using ASIC chips instead of general purpose DSP chips in the existing CTs). It is proposed here to improve modification 0, 1 and II over existing methods. In this case, existing methods rely only on pilot signals transmitted from cellular sectors to locate the CT. More specifically, existing methods describe a typical operation of a CT, which searches for pilots and determines their TOAs. However, this normal mode of operation of a CT is designed to maximize the performance of the CT in receiving data signals, and to maximize the communications performance (i.e. capacity of the system, reliability of data transmission, etc.). Improvement 0 introduces Location state 0 which the CT enters into whenever it is requested to locate itself in modification 0. Location state 0 optimizes several design parameters in modification 0 as follows: Improvement 0: The most common method to estimate the TOAs, POAs, Strength of Arrival (SOAs), and FOAs of received signals at the MR consists of correlating a received signal with a clean version of itself. Conventional methods cannot resolve correlation peaks beyond the conventional time- domain Rayleigh resolution. Moreover, conventional methods either choose the highest correlation peak above a certain threshold or the peak corresponding to the first arrival above a given threshold. Over a multipath channel, it is possible to assume that the first peak above a certain threshold corresponds to the correct peak, however, aliasing, noise, and interference can cause spurious peaks to appear in time before the correct one. It is an object of Improvement O to
1. optimize the threshold used in selecting the correlation peak; and to
2. optimize the list of BSs considered for correlation.
Such an optimization is based on minimizing HDOP. The optimized parameters: threshold and list of BSs, are then transferred to the CT whenever it is requested to enter into Location state 0. In such a state, the CT is forced to receive more signals from more BSTs in order to reduce HDOP. The optimization process can be explained as follows. In order for the CT to receive more signals, its threshold has to be lowered. A lower threshold however forces the CT to accept weak signals. This produces low SNR signals that might suffer from multipath and interference. In order to minimize HDOP while at the same time reducing the effect of noise and multipath, the optimization process consists of solving for the position of the MR using Least Squares weighted by the individual RSSI from each received signal. 0224
67 In modification I, existing methods require a great deal of standard change. Instead, Improvement I does not require a standard change. It forces the CT to enter into Location state I:
Improvement I: It consists of using the SR and Inverse SR algorithms over the correlation function (that is obtained by correlating the received signals at both the MR and RR with a clean version of themselves) in order to:
1. resolve the correlation peaks beyond the conventional time-domain Rayleigh resolution; and to
2. avoid using a threshold in selecting the correlation peak.
An important feature in the Inverse SR algorithm is the fact that it generally produces a number of correlation peaks:
• a peak corresponding to the direct path signal as well as,
• a number of spurious peaks caused by noise, aliasing, multipath and interference, that can surround the correct peak.
Instead of using a threshold, Location state I consists of selecting the correlation peak that produces the smallest residual when performing a weighted Least Squares operation on all possible correlation peaks resulting from super-resolving the correlation functions produced from the received signals using the inverse SR algorithm. It is also an object of Location state I to optimize the list of BSs considered for correlation based on minimizing HDOP.
In modification H, existing methods require a great deal of standard change as well as CT change. More specifically, it consists of receiving GPS signals partly at the CT and partly at the RR. Instead, several improvements over the existing methods are offered in modification H:
• It allows the CT to choose between receiving cellular signals, or receiving non-cellular signals (such as GPS, Loran-C, etc.), or both;
• It allows the CT to use Improvement 0 and/or Improvement I, i.e. to enter into Location state 0 and/or Location state I;
• It allows the CT to use the carrier phase of the received signals for TDOA location;
• It allows the CT to use AOA for GPS Loran-C;
• It allows the CT to locate itself using only two independent BSs, with a minimum of two sectors per BS, when using AOA estimation;
• It allows the CT to use dual-frequency receivers for GPS;
• It allows the CT to use Pseudolite in GPS; It allows the CT to use cellular augmented GPS;
It allows the CT to use short messaging cither for transmitting its Location Information or for receiving the Reference Information, so as not to be network-dependent;
It allows the CT to combine location information from both GPS and cellular signals. GPS alone requires four satellites in order to solve for (x,y,z,) and the CT's clock error. Cellular alone require three BSs in order to solve for (x,y) and the CT's clock error. When combining GPS equations with cellular equations, we require four equations in total since we still have four unknowns: x, y, z, and the CT's clock error. A possible alternative is to replace z by the altitude of the BSs. However, this can introduce an error in the estimation for x and y;
Description of an IS-95 CT receiving Pilot Signals:
Here it is proposed to use an IS-95 receiver which measures the TOA of multiple pilot signals for use in estimating the location of said receiver. This receiver may be obtained through modifications to the operation of a wide variety of CDMA receivers, including CDMA CTs.
An IS-95 CT is directed to enter location state 0 or Location State I by an order received from:
1. a remote location,
2. an order received from the user of the CT, or
3. by other predetermined circumstances, such as initialization on power- up, or origination of a CT.
Upon entering location state 0 or I,
1. the CT assembles a list of pilots to search for (1011) as shown in Figure 10. This list may be formed from a combination of the CT's Active set, neighbour set, remaining set, and a list of pilots specified in a message received by the CT from a remote source.
2. The process of searching for pilots, and assigning a demodulator finger to obtain a refined TOA estimate, is executed for each pilot in the pilot search list. In blocks (1012), (1013), (1014), (1015), (1016), (1017), and (1018), P represents one of the pilots from the list. Block (1012) initializes P to the first pilot in the search list. 3. A searcher element is assigned to search for pilot P in the received signal (1013). A window of PN offsets (start and end, or center and width) is specified to the searcher. The searcher then computes received energy of pilot P at multiple delay offsets within the search window. In the preferred embodiment, these delay offsets are spaced one half chip apart (0.4069 microseconds). For each delay offset, the searcher accumulates received pilot energy over a specified duration of time, (preferred embodiment uses 0.5 ms)
4. The searcher results are examined for evidence of reception of a pilot signal with sufficient quality (1014). In the preferred embodiment, the ratio of the average pilot chip energy to the total received power (called EJI is computed at each delay offset. Figure 11 illustrates a possible set of values for searcher results. In this example, significant pilot energy appears near delay offset 7. The pilot EJl0 at each delay offset is compared to a threshold value (-15 dB in the preferred embodiment). If the fl^ at one or more delay offsets exceed the threshold, then the earliest such delay offset is used as an initial estimate of the time of arrival of the pilot. If the threshold is not exceeded at any of the delay offsets, men it is determined that the pilot was not detected.
5. If pilot P was detected, then a demodulator finger is assigned to pilot P at the initial estimate of the time of arrival of said pilot (1015). The finger is allowed to track pilot P for a period of time (preferred embodiment uses 2 ms), after which
6. the delay offset of the finger is retrieved and stored as the time of arrival of pilot P (1016).
7. After the searcher is finished searching for pilot P, P is assigned to the next pilot in the list (1017), (1018).
8. If the list is exhausted, then the search results are assembled (1019) for subsequent processing or transmission to a remote site.
Figure 10 illustrates steps performed sequentially. Some steps could be performed in parallel. Specifically, the searcher can begin to search for another pilot while a demodulator finger tracks a previously found pilot. If the CT's time reference is not sufficiently stable, then one demodulator finger can be assigned to one pilot while other pilots are searched and tracked by other fingers, in order to maintain a known time reference.
A demodulator finger may track a pilot using a delay-locked loop, as described in U.S. Patent 5,764,687. This results in the finger moving to the peak of the arrival of the pilot. Improved performance may be obtained by having a finger track towards a predetermined E/L0 threshold on the rising edge of the pilot energy versus delay offset.
Signal processing:
The location operation is performed at the Host which can be either located within the MR (Option I) or remotely from the MR (Options H and DT) (e.g. within the RR or a TR). In Options II and HI, it is preferable to perform as much of the processing as possible at the MR in order to reduce the amount of information transmitted from the MR to the Host. For example, if TDOA is employed for location, it is preferable to measure the TOA of a received signal, r;(t), that is transmitted from the iΛ BST, at the MR. The same is true for PDOA and FDOA. In the TOA case and in the POA case, the required processing consists of correlating the received signal with a clean version of p;(t) and sending the information regarding the correlation peaks to the Host. In the FOA case, the required processing consists of correlating the received signal with several clean versions of the p;(t) at different frequencies and sending the information regarding the correlation peaks to the Host.
The correlation peaks between r;(t) and a clean version of p;(t) can be super- resolved using SR algorithms. It is sometimes preferable to perform the SR algorithms at the Host rather than performing them at the MR since it requires a number of computationally intensive operations that are easier to develop at the Host than at the MR. This implies that the correlation peak needs to be transferred to the Host which does not represent a large amount of information and can be easily transferred to the Host.
In order to be able to correlate the received signal, r,(t), at the MR, a clean version of p,(t) must be available at the MR beforehand. In order to avoid correlating r,(t) with p,(t), one can use either one of the following two procedures: Procedure I:
1. Filter the entire received signal;
2. Extract the phase of an existing tone. In this case, TOA and AOA information are contained in the phase of the tone while the FOA information are contained in the frequency of the tone and can be extracted using Procedure HI (see below).
Procedure π:
1. Filter the entire received signal;
2. Perform a nonlinear operation on the filtered signal;
3. Extract the phase of a tone generated by the nonlinear operation in step 2. In this case, TOA and AOA information are contained in the phase of the tone while the FOA information are contained in the frequency of the tone and can be extracted using Procedure DH.
The nonlinear operation in Procedure H can be either a quadratic operation (e.g. squaring the signal), a delay product operation (i.e. multiplying the signal by a delayed version of itself) or any other suitable nonlinear operation.
Procedures I and H can be used both for acquiring the signal and for estimating its TOA. This is accomplished by stamping the time using a reliable clock at either the start of the generated tone or obtaining the time stamp at the end of the observation window (the latter is preferred).
Procedures I and H can be improved by resolving the existing or generated tones using SR and inverse SR algorithms. Moreover, the filtering required in step 1 of Procedure II can be carried out in the digital domain in order to reduce the effect of temperature and aging on group delay variation.
Procedures I and H can be improved by repeating them a number of times and then either averaging the Location Information or averaging the TDOA, PDOA or FDOA between pairs of BSTs.
An alternative method for improving Procedures I and π is by substantially increasing the duration of observation, T;, of the received signal, r,(t), thereby increasing the received SNR. In order to be able to increase T, at will, p;(t) has to be available at the MR. An important feature in Procedure II is the fact that the signal r,(t) transmitted by the i* BST does not have to be a deterministic signal, i.e. a clean replica of the signal does not have to be available at the MR.
Equivalently, an important feature in Procedure II is the fact that the probability of location is 100%, unlike Procedure I which is below 100%.
This leads to the fact that in Procedure H, increasing the duration of observation, Tj5 can accommodate different situations. For example, in rural areas, where the cells are large, the SNR can be increased by increasing the duration of observation. Or equivalently, if the geometry is not good (i.e.
HDOP is large, e.g. highway), one can improve it by increasing the duration of observation hence, increasing the number of BSTs that receive the signal Tj(t) at a reasonable level. The duration of observation is thus preferably optimized with respect to the SNR and HDOP. In addition to reducing HDOP by increasing the number of BSTs, this also leads to reducing the effect of multipath, particularly in urban centers.
In order to reduce the effect of in-band interference on the phase of the desired tone, narrowband digital filters of a few hertz bandwidth are used (e.g. Fourier Transform-based filters, sliding Fourier Transform-based filters, ITR filters, FIR filters, etc.). The Fourier Transform-based filter has a bandwidth which is a function of the observation duration, T{. The ITR filter has a steady state response time which is a function of T;. In other words, the observation duration, T;, has to be selected in order to optimize the SNR (increase it) versus the bandwidth of the filter (decrease it) versus the response time of the wireless location system (minimize it).
In order to reduce the effect of the multipath, MPj(t), on the performance of a kinematic wireless location system, it is desirable to resolve the different Doppler frequencies δξ (which correspond to different TOAs and different AOAs) and to choose the frequency corresponding to the first arrival. The method of resolving the different frequencies can be accomplished using SR algorithms.
Hardware
In order to improve over prior art for modification II, one should attempt to 1. reduce or avoid: — Gain Imbalance;
— Phase Imbalance;
— Carrier Feed-through; and
— DC offset,
2. be able to use the same RF front end for all standards (analog or digital) over a fixed band; and
3. reject out-of-band Interferers (e.g. paging, trunked radio, etc.) using digital filters in the Digital Signal Processor (DSP).
To achieve the three objectives above, it is proposed to use a linear IF- sampling receiver for the MR, or for the RR, or for both. Such an TF-sampling receiver is designed to have high rejection, low group delay variation and good sensitivity. In order to have high rejection with good sensitivity, the receiver is designed to have (see Figure 8): Design I:
1. A number of Intermediate Frequency (TF) stages (810) that convert the received RF signal Re{ r f)exp(j2πf;f)} (801) to an B? signal (807) of IF frequency f Initial rejection is achieved at RF using an RF Band Pass Filter (BPF) (802), followed by an RF amplifier (803) for good sensitivity. Subsequent rejections are achieved at each TF stage which consist of a mixer (804) followed by an Intermediate Frequency (W) BPF filter (805) and an IF amplifier (806).
2. A pre-aliasing filter (811) precedes the Analog-to-Digital converter (812) that converts the TF analog signal to an TF digital signal. The sampling rate is intended to be below the Nyquist rate in order to force aliased copies of the IF signal (807) close to a lower digital IF, f2.
3. Digital filters (813) which reject noise and interference outside of the digital TF bandwidth.
4. A Digital Down Converter (DDC) (804) which converts the (real) digital B? signal to a (complex) digital baseband signal (optional).
In order to have low group delay variation, the following procedure is suggested:
Procedure DTI:
1. Estimate the carrier offset (818), (due to Doppler, δf„ and due to LO offset, Δf,) of the signal r,(t) (801) using Digital Signal Processor
(DSP) (817). 2. Adjust the Local oscillators (808) in the receiver corresponding to the MR during reception of r;(t) (801) based on the estimated carrier offset (818), in step 2 above.
Explanation of Procedure HI: Since the overall group delay, Δtj, is a function of frequency, it can vary significantly from one frequency to another. Furthermore, according to the AMPs standard, Δfj can be as large as 2kHz. For this reason, it is important to force Δtj to be independent of the carrier offset, Δfj + δjζ. This is accomplished by estimating the carrier offset, Δξ + δf;, in step 1 of Procedure HI and removing it in step 2 of Procedure HI prior to processing it for location purposes. In other words, steps 1 and 2 force the overall group delay, Δtj, to remain constant, thereby reducing overall group delay variations due to frequency offsets. In order to reduce overall group delay variations due to temperature and aging in the RR, a calibration procedure has to be performed on a regular basis on the RR.
No tes:
1. In Procedure HI, step 1 can be performed using a discrete Carrier Recovery Loop (CRL) described in Procedure IV (see below).
2. In Procedure HI, step 2 can be performed using a Direct Digital Synthesiser (DDS) (809).
The discrete Carrier Recovery Loop consists of several steps that are repeated a number of times until a satisfactory level of performance is reached: Procedure IV:
1. Digitally Signal Process the sampled IF signal (814 or 816) using the DSP (817) in order to derive an adequate objective function, fobj.
2. If fobj is optimized to within a certain constraint, stop the loop, otherwise:
3. search for a new frequency offset (818),
4. adjust the LO (808) by the new frequency offset (818), and
5. go back to step 1.
In Procedure TV, step 3 can be performed using a gradient-type search algorithm or any other suitable search algorithm. In Procedure IV, step 4 can be performed using cither a DDS (809) only or a
DDS and a Digital Down-Converter (DDC) (815). If a DDS (815) only is used in step 4, then step 1 in Procedure IV processes a new segment of the signal, r,(t), every iteration. If a DDS (809) and a DDC (809) are used is step 4, then step 1 in Procedure IV processes the same stored segment of the signal, r,(t), every iteration until the DDS (809) is used. It is preferred to use both the DDC and the DDS in order to reduce the effect of finite precision arithmetic on the estimation process.
In the case when Design I does not contain step 4, the DSP (817) in step 1 of Procedure IV of the carrier recovery loop can consist of a digital notch filter centered at the desired IF frequency f2, whose power of response is to be minimized, i.e. f^- is the response of the digital notch filter. In this case, the notch filter can consist of
1. a first digital Band Pass Filter centered at the f2 + frequency of desired Tone; and
2. a second digital Band Pass Filter centered at the f2 - frequency of desired Tone.
The Band Pass filters can be implemented using (real) digital Biquadratic filters.
Alternatively, the DSP (817) in step 1 in Procedure IV of the carrier recovery loop can consist of a digital Band Pass Filter centered at the desired IF frequency f2, whose power of response is to be maximized, i.e. f^- is the response of the Band Pass Filter. In this case, the BPF can be implemented as a (complex) 1st order (or higher) digital HR filter.
In the case when Design I contains step 4, the signal to process is a baseband signal (816) and the digital BPF turns into a simple (complex) LPF. An example of such a LPF is an averaging filter which adds all real samples together and all imaginary samples together over the observation interval.
In order to reduce the effect of clock errors and LO offset (or drift), Δζ, a common reference is used for all the LOs and for all the clocks. A suitable common reference for the LOs and for the clocks can be derived from the GPS signal. Given that either all LOs in the BSTs have a common reference, which implies that
Δf, = Δf2 = Δf, = Δf
in equations (32), or that the RR estimates Δf,, Δf2, Δf3 in order for them to be removed, it is possible to estimate the speed, v, with respect to the BSTs and Direction Of Travel (DOT), φ, of the MR from the carrier offsets that are estimated from Procedure V at each BS using the discrete carrier recovery loop. This can be accomplished using the Frequency Difference of Arrival (FDOA) between each pair of BS as expressed in equations (33).
Calibration methods at the RR
It is preferred to minimize the use of the RRs in order to reduce cost and complexity of a handset-based WLS. It is desired to produce a handset-based Wireless Location System (WLS) where an MR can be located without the need for an RR. This is achieved by considering a transmitter in an existing cellular sector as a BST. The advantage is that such a sector is usually time- referenced to other sectors belonging to the same BST, i.e. sectors in the same BS are more or less synchronized to one another. This implies that the MR can estimate the AOA of signals transmitted by such sectors without a need for an RR. Moreover, the MR can estimate its location using only two independent BSs based on the individual estimated AOA from each BS without the need for an RR_
A wireless location system has to be calibrated prior to operation (and occasionally during operation) in order to remove the effect of the overall Group Delay, Δtj, at each RR. Moreover, temperature and aging can cause the overall Group Delay to change with time. Therefore, it is necessary to calibrate the RR on a regular basis. A method to calibrate an RR is as follows: Procedure V:
1. tiansmitting at precise time intervals signals from a transmitter of known location to the RR over several frequency bands. Both the RR and the transmitter have to be referenced to the same clock source;
2. characterizing the frequency response of the RR;
-SUBSTITWE 6H5&T (RW_E*6fr»« 3. estimating the time of travel between the transmitter and the
RR ;
4. estimating the group delay of the RR.
Procedure V can be repeated a number of times in order to average the relative group delay over time. In order to reduce deployment cost of a transmitter at known locations, an MR is stored in proximity to the RR. This allows for automatic calibration without a need for human intervention by forcing the MR to act as a transmitter.
In order to minimize the relative group delay variations between calibration intervals, it is recommended to use SAW filters as the BPF (803, 805, 811) in Figure 8. It is also recommended to place the SAW filters (803, 805, 811) in a metallic block to minimize flexing effects due to temperature variations.
Host Software
Given a number of positioning algorithms and a number of methods for blunder detection, positioning strategy is used as illustrated in Figures 9a, 9b, 9c and 9d. Figures 9 illustrate the preferred embodiment for adjusting the observed TOAs such that the sum of the squares of the residuals is a minimum. Other methods for adjusting the observed according to some optimization criteria are available mcluding:
• General Minimum Variance Unbiased Estimation,
• Best Linear Unbiased Estimation
• Maximum Likelihood Estimation,
• Method of Moments,
• General Bayesian Estimation,
• Linear Bayesian Estimation,
• Kalman Filtering, etc.
The positioning process begins with a set of TOA observations from a minimum of three BSs. If only three TOAs are available (901), solution bifurcation is tested for (902). In the event of solution bifurcation, the LOCA algorithm (903) is performed on the TOA observations yielding two solutions. These two solutions are each used as the initial position for Least Squares (904 and 905) resulting in two Least Squares solutions for the MR. The Least Squares algorithm here, and in all other parts of Figure 9, is Hybrid TDOA positioning with received signal strength used to weight the individual BS observations. During every iteration of Least Squares, misclosures for each
BST are calculated. If any misclosure exceeds some multiple of the RMS value of the entire set of misclosures, the corresponding BST, and its TOA, is no longer used within that particular execution of Least Squares.
If both Least Squares solutions (904 and 905) converge (906) but are not equal (907), it is assumed that Least Squares has converged to the two possible solutions and both Least Squares solutions are reported as possible positions for the MR (908). If the either of the LS solutions diverge (906) or the convergent solutions are equal (907), the two LOCA solutions from (903) are reported as possible positions of the MR (909).
If at (902) bifurcation is found to not exist, LOCA is again executed (910) but yields only one solution. This solution is used as the initial position in Least Squares (911).
When redundancy does exist (901), Plane Intersection (912) (or any other closed-form position estimation algorithm) is executed. The resulting solution is used as the initial position for Least Squares (913). At this point a check is made as to whether Least Squares (913 or 911) converged (914). If Least Squares converges, another check for observational redundancy is made (915). If there is no redundancy, the Least Squares solution is reported as the position estimate of the MR (923). If there is redundancy, the normalized residuals are statistically tested for normality (916). The Least Squares solution is reported (923) should all residuals pass. If any residuals fail, the redundancy numbers of the failing observations are checked (917). The redundancy number of the ift observation is defined as
Figure imgf000080_0001
where C? is the covariance matrix of the residuals and is defined as
CF = C, - A ATC, 'A A 1 (53) Should all failing observations have redundancy numbers less than some threshold (preferably, this threshold is 0.5), the Least Squares solution is reported as the MR position estimate (923).
If observations with failing residuals have redundancy numbers larger than the threshold (917), that observation with redundancy number greater than the threshold and with the largest standardized residual is permanently removed from the observation set (928). The initial position used for the previous execution of Least Squares is remembered (929) and used again in Least Squares with the truncated observation data set (930 or 931).
Should Least Squares now diverge (932), the previous Least Squares solution which did converge is reported as the position estimate of the MR (936). If Least Squares does converge (932), and there is no redundancy (933), the newly convergent Least Squares solution is reported (936). If there is redundancy (933), the standardized residuals are tested for normality (934). Should all the residuals pass, the newly convergent Least Squares solution is reported (936). Otherwise, if all failing observations have redundancy numbers less than some threshold (preferably, this threshold is 0.5), the newly convergent Least Squares solution is reported (936).
If observations with failing residuals have redundancy numbers larger than the threshold (935), that observation with redundancy number greater than the threshold and with the largest standardized residual is permanently removed from the observation set (928). The process then continues as described immediately above.
If at (914) Least Squares diverges, Least Squares is executed again but with the average of participating BSTs coordinates used as the initial position (918). If Least Squares now converges (919) but there is no observation redundancy (920), the newly convergent Least Squares solution is reported as the position estimate of the MR (923). If there is redundancy (920), the standardized residuals are tested for normality (921). Should all residuals pass, the newly convergent Least Squares solution is reported (923). If some of the residuals fail but none of the corresponding observations have redundancy numbers greater than some threshold (922), the newly convergent Least Squares solution is reported (923). If observations with failing residuals have redundancy numbers larger than the threshold (922), that observation with redundancy number greater than the threshold and with the largest standardized residual is permanently removed from the observation set (928). The process then continues from (928) as described above.
When Least Squares does not converge at (919), and there is no redundancy (924), either no solution is reported for this particular set of observation data or the solution from the last iteration before divergence is reported (927). If there is redundancy, but all of the standardized residuals pass the normality test (925), no solution is reported for this particular set of observation data or the solution from the last iteration before divergence is reported (927). Should some of the residuals fail, the observation with the largest standardized residual is permanently removed from the data set regardless of its redundancy number. The process then begins at point (912) as described above.
Figures 9 are also applicable for,
• AOA positioning,
• AOA TDOA positioning, and
• AOA/Range positioning.
Transferring Location Information of the MR over the Internet
Given that the Internet is global and inexpensive, the communication between the Host and the customer can be achieved over it. For example when a MR calls 911 for emergency, it is possible to relay the positional information of the MR from the Host to the PSAP over the internet. Similarly, in the case when a powered-on MR is to be located by a customer, its positional information can be relayed from the Host to the customer through the internet.
The utility of the internet in its capacity of transferring MR location information on a global scale, extends the application of the wireless location system beyond the cellular network for which the BSTs have been deployed. Through the latest internet technology such as Java, JavaBeans, as well as CORBA (Common Object Request Broker Architecture), MR location information residing at the Host can be integrated with third-party information (i.e. a map database, or a database consisting of geo-coded business addresses such as restaurants, towing companies, etc.) resid g in some other geographical location, perhaps even in another country. The combination of this MR location information with the third party information can be transferred over the internet to allow customers to locate a particular MR wrt to either a map location (i.e. street address) and/or a business location, provided that the MR is within network coverage for the wireless location system. Through this process the MR location information and the third party database information can be accessed or "served up" to the customer through efficient Java Internet technology processes. The service provider will integrate together the various components including the MR location information derived from the wireless location system and the third-party information database. The integration process of the different databases is transparent to the customer. The customer will only know that the location-based service exists as such to bring all of the different components together to provide for a complete service which can be offered either on a regional, national or perhaps global basis. Through the internet these location-based services can now be offered on an economical basis to the customer. Examples of such services are fleet management, concierge services, roadside assistance, child find services, etc. Legahty and security aspects are a concern on the internet and in some cases a dedicated link may be necessary.
A clean replica referred to herein is a copy of a signal that is not noisy, as for example may be obtained by a stored version of the signal transmitted. Low group delay variation preferably means a few 10s of nanoseconds of group delay. An adequate objective function is a function known in the literature, for example minimizing the mean squared error. A PSAP is a public safety answering point, for example a 911 station. An RSSI is a receiver signal strength indicator.
Every document referred to herein as containing useful information for implementing MR location is incorporated herein by reference.
Modifications to the preferred embodiments disclosed may occur to those skilled in the art without departing from the scope of the invention.

Claims

Wc Claim
1 A method of estimating the location of a Mobile Receiver (MR) comprising the steps of. receiving a plurality of signals from a plurality of Base Station Transmitters
(BSTs) by the MR; and performing a nonlinear operation on the plurality of received signals at the MR in order to generate at least one tone per received signal.
2. The method of claim 1 wherein the nonlinear operation is selected from the group consisting of a quadratic operation, a delay product operation; and an absolute value operation.
3. The method of claim 1 where the location operation is performed by estimating the phases of the generated at least one Tone.
4. The method of claim 1 wherein the received signals are filtered prior to performing the nonlinear operation.
5. The method of claim 4 wherein the received signals are correlated with a clean replica of themselves prior to performing the filtering operation.
6. A method of estimating the location of a Mobile Receiver (MR) comprising the steps of receiving a plurality of signals from a plurality of Base Station Transmitters
(BSTs) by the MR; and estimating the location of the MR by estimating the phases of existing tones in each received signal.
7 The method of claims 1 or 6 where the geneiated oi existing Tones are resolved using SR or Inverse SR algorithms
8 The method of estimating the location of the MR in claims 1 oi 6 by extracting the Tones using digital filteπng
9. The method of claim 8 whcicin the method of digital filtering compi iscs the steps of performing a Fouiier Transfor m (FT) on segments of the srgnal, windowing the desired frequency domain bands for each FT segment; and performing an inverse Fourier Transform (IFT) on the windowed bands
10. The method of claim 9 wherein the windowed frequency domain bands for each FT segment are transferred to a Host.
11. The method of claim 8 wherein the method of digital filtering performs an operation selected from the group consisting of:a sliding (or recursive) FT operation, an IIR filtering operation; and an FIR filtering operation.
12. The method of claim 10 where a time domain signal, r,(t), is constructed at the Host by appending all the received signal frequency domain segments; and performing an IFT operation.
13. The method of claims 9 or 12 where a plurality of tones in the reconstructed time domain signal, r,(t), are resolved using SR or Inverse SR algorithms.
14. The method of claim 13 where a tone is selected from the plurality of resolved tones such that it corresponds to the earliest Time-Of-Arrival among all resolved tones.
15. The methods of claims 7 or 13 where the SR or the inverse SR algorithms is selected from the group of algorithms consisting of MUSIC, ESPRIT, Auto Regressive Moving Average, Minrmum Var iance, MUSIC using Higher Older Statistics, ESPRIT using Higher Order Statistics, Auto Regressive Moving Aveiage using Higher Order Statistics, and Minimum Vaiiancc using Higher Older Statistics
16 The method of claim 1 5 further compπsing the method pei forming an Im ei se SR algonthm, the mveise SR algonthm compr ising a conventional time domain conelatoi , a time domain window, a Fourier Transform, a frequency domain window, a frequency domain equalizer, and a processor which performs a SR algorithm in order to resolve the TOAs in the received signal at a given BST.
17. The method of claim 15, further comprising the method performing an Inverse SR algorithm; the inverse SR algorithm comprising: a frequency domain correlator, a frequency domain window, an inverse Fourier Transform, a time domain window, a time domain equalizer, and a processor which performs a SR algorithm in order to resolve the FOAs in the received signal at a given BST.
18. A method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where one or more antennas is used at the MR and in which location of the MR is achieved using a technique selected from the group consisting of:
AOA positioning (multi-angulation), TOA positioning (circular multilateration),
Range positioning (circular multilateration), TDOA positioning (hyperbolic multilateration), TDOA positioning (circular multilateration with TDOAs), AOA/TDOA positioning (multi-angulation/hyperbolic multilateration), AOA/TDOA positioning (multi-angulation/circulai multilateration with
TDOAs),
AOA/Range posrtioning (multr-angulation/circular multilateration), AOA/Range positioning (multi-angulation/circulai multilateration), and combining the signals fi om all antennas foi location pui poses
19 The method of claim 18 in which moie than one antenna is used
20. The method of estimating the location of a MR in claim 18 where the combining of the signals from all antennas is selected from the group consisting of selection combining, maximal ratio combining, co-phasing combining, and equal gain combining.
21. A method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, using equation (29) and equation (30).
22. The method of claim 18 where location of the MR is achieved using a technique selected from the group consisting of:
Range positioning (circular multilateration) with Round Trip Delay, and AOA/Range positioning (multi-angulation/circular multilateration) with
Round Trip Delay.
23. A method of estimating the location of a MR wherein location of the MR is achieved based on a solution selected from the group consisting of: solution I, solution II, solution III, and solution IV.
24. A system to estimate the location of a Mobile Receiver (MR) comprising of a receiver at the MR for receiving a plurality of signals from a plurality of Base Station Transmitters (BSTs); and a computing means connected to the receiver by a communication linkfor performing a nonlinear operation on the plurality of leceived signals in order to generate at least one tone per received signal.
25 A system to estimate the location ol a Mobile Reccivei (MR) compnsing of a ieceivei at the MR for leceix mg a plurality of signals from a plurality ol Base Station Tiansmitteis (BSTs). and a computing means connected by a communication link to the receiver for estimating the location the MR by estimating the phases of existing tones in each leceived signal
26. The system of claim 24 further comprising a computing means for performing an Inverse SR algorithm; the inverse SR algorithm comprising: a conventional time domain correlator, a time domain window, a Fourier Transform, a frequency domain window, a frequency domain equalizer, and a processor which performs a SR algorithm in order to resolve the TOAs in the received signal at a given BST.
27. The system of claim 24 further comprising a computing means for performing an Inverse SR algorithm; the inverse SR algorithm comprising: a frequency domain correlator, a frequency domain window, an inverse Fourier Transform, a time domain window, a time domain equalizer, and a processor which performs a SR algorithm in order to resolve the FOAs in the received signal at a given BST.
28. A system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where one or more antennas are used at the MR and wherein the system comprises means for implementing a technique selected from the group consisting of:
AOA positioning (multi-angulation),
TOA positioning (circular multilateration),
Range positioning (circular multilateiation),
TDOA positioning (hypeibolic multilateration). TDOA positioning (circular multilateration with TDOAs),
AOA/TDOA positioning (multi-angulation/hypei bolic multilateiation), AOA/TDOA positioning (multi-angulation/cncular multilateiation with TDOAs),
AOA/Range positioning (multi-angulation/cncular multilateration), AOA/Range positioning (multi-angulation/cncular multilateration), and combining the signals from all antennas for location purposes
29. The system of claim 28 in which more than one antenna is used at the
MR
30. A system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, using equation (29) and equation (30).
31. The system of claim 28 further comprising means for implementing a technique selected from the group consisting of: Range positioning (circular multilateration) with Round Trip Delay, and
AOA/Range positioning (multi-angulation/circular multilateration) with Round Trip Delay.
32. A system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, based on any one of the following solutions: solution I, solution II, solution HI, or solution IV.
33 A system to estimate the location of a Mobile Receiver (MR), the system compiising a Receiver that receives a plurality of signals transmitted from a plurality of BSTs wheiein the Receiver includes an IF-samphng Receiver comprising a number of IF stages to conveit a received RF signal r,(t) to an analog IF signal of IF frequency fi ; a pre-aliasing filter; an Analog-to-Digital (A/D) converter that converts the IF analog signal lo an II- digital signal; and a piocessor for processing the IF digital signal
34. The system of claim 33 in which each IF stage of the IF-sampling receiver comprises: a mixer; an IF Bandpass Filter (BPF) filter; and an IF amplifier.
35. The system of claim 33 wherein the sampling rate of the A/D converter is below the Nyquist rate in order to force aliased copies of the digital IF signal close to a lower Intermediate Frequency f2.
36. The system of claim 33 wherein the processor performs functions selected from the group consisting of the functions: digital filtering to reject noise and interference outside of the digital IF bandwidth; digital down conversion to convert the digital IF signal to a complex digital baseband signal; and discrete carrier recovery in a carrier recovery loop.
37. The system of claim 33 wherein the IF-sampling Receiver has low group delay variations over a desired frequency band.
38. The system of claim 33 where the IF-sampling Receiver comprises one or more local oscillators (LOs) and clocks and has a common reference for all the LOs and for all the clocks.
39. The system of claim 38 wherein the common reference for the LOs and for the clocks is GPS-based.
40. The system of claim 34 wherein the pre-aliasing filter is a SAW filter and the BPF is a SAW filter.
41 . The system of claim 40 wherein the SAW filters arc placed in a metallic block to minimize flexing or temperature effects.
42. The system of claim 33 in which the Receiver is the MR.
43 The system of claim 33 in which the Receiver is a Reference Rcceivei distinct from the MR
44 A system to estimate the location of a MR using a Reference Receiver having a LO, the system comprising an estimator for estimating the carrier offset (due to Doppler, δf,, and due to LO offset, Δf,) of the plurality of signals transmitted by a plurality of BSTs and received by the RR
45. The system of claim 44 further comprising a discrete Carrier Recovery Loop.
46. The system of claim 45 wherein the carrier recovery loop comprises a processor configured to repeat the following steps until a satisfactory level of performance is reached:
1. process the sampled IF to derive an adequate objective function, fob,;
2. if fob is optimized to within a certain constraint, stop the loop, otherwise:
3. search for a new frequency offset, 4. adjust the LO by the new frequency offset, and
5. go back to step 1.
47. The system of claim 46 wherein the search for a new frequency offset in step 3 is performed using a gradient-type search algorithm
48 The system of claim 46 wherein adjusting the LO by the new frequency offset can be performed either using a DDS or a Digital Down-Converter or a combination thereof
49 The system of claim 36 or 46 wherein the piocessoi performs the function of a earner recovery loop by performing a digital notch filtei g operation centered at the desired IF frequency f2, and whose power of lesponse is an optimized objective function, fob,, derrved from the sampled IF
50 The system of claim 49 in which the digital notch filtei mg operation consists of a first digital Band Pass Filter centered at the f2 + frequency of desired Tone; and a second digital Band Pass Filter centered at the f2 - frequency of desired Tone.
51. The system of claim 36 or 46 wherein the processor is configured to perform the function of a carrier recovery loop by performing a digital Band Pass filtering operation centered at the desired IF frequency f2, and whose power of response is an optimized objective function, f0bj, derived from the sampled IF signal.
52. The system of claim 51 wherein the desired IF frequency f2 is equal to zero and the digital BPF is a LPF.
53. The system of claim 51 wherein both Band Pass filters are implemented using digital Biquadratic filters.
54. The system of claim 46 further comprising a Host and wherein the estimated carrier offset last found is transferred to the Host in order to aid in estimating the speed v of the MR and its Direction Of Travel. ************************************
55. A system to estimate the location of a MR having a LO, the system comprising: an estimator for estimating the carrier offset (due to Doppler, δfj, and due to LO offset, Δfj) of a plurality of signals transmitted by a plurality of BSTs and received by the MR.
56. The system of claim 55 in which the estimator comprises a discrete Carrier Recovery Loop.
57. The system of claim 56 wherein the carrier recovery loop comprises a processor configured to perform the following steps until a satisfactory level of performance is reached:
1 . process the sampled IF to derive an adequate objective function, f h,-
2. if f0h is optimized to within a certain constraint, stop the loop, otherwise. 3. search for a new frequency offset,
4. adjust the LO by the new frequency offset, and 5 go back to step 1
58. The system of claim 57 wherein the search for a new frequency offset in step 3 is performed using a gradient-type search algorithm
59. The system of claim 57 wherein adjusting the LO by the new frequency offset is performed either using a DDS or a Digital Down-Converter or a combination thereof.
60. The system of claim 46 wherein the processor is configured to perform the function of a carrier recovery loop by performing a digital notch filtering operation centered at the desired IF frequency f2, and whose power of response is an optimized objective function, f0bj, derived from the sampled If signal.
61. The system of claim 59 in which the digital notch filtering operation consists of: a first digital Band Pass Filter centered at the f2 + frequency of desired Tone; and a second digital Band Pass Filter centered at the f2 - frequency of desired Tone.
62. The system of claim 61 wherein the processor is configured to perform the function of a carrier recovery loop by performing a digital Band Pass filtering operation centered at the desired IF frequency f2, and whose power of response is an optimized objective function, f0b,, derived from the sampled If signal.
63. The system of claim 62 wherein the desired IF frequency f2 is equal to zero and the digital BPF is a LPF.
64. The system of claim 61 wheiein both Band Pass filters aie implemented using digital Biquadratic filters.
65 The system of claim 55 oi 56 fui thei compr ising a Host and wheiein the estimated carrier offset last found is transfer red to the Host in oidei to estimate the speed v of the MR and its Dnection Of Tiavel
66. The system of claim 65 wherein the speed of the CT and its direction of traveling is estimated as a function of the difference between the estimated carrier offsets.
67. A method of estimating the location of a MR, the method comprising the steps of: receiving a plurality of signals transmitted from a plurality of BSTs at a Receiver; converting the received RF signal r,(l) to an IF signal of IF frequency fi in an IF- sampling Receiver; filtering the IF signal with a pre-aliasing filter; converting the filtered IF analog signal to an IF digital signal in an Analog- to-Digital converter; and processing the IF digital signal in a signal processor.
68. The method of claim 67 in which the Receiver is the MR.
69. A method of estimating the location of a MR, the method comprising the steps of: estimating the carrier offset (due to Doppler, δfj, and due to LO offset, Δfj) of a plurality of signals transmitted by a plurality of BSTs and received by a RR.
70. A method of estimating the location of a MR, the method comprising the steps of: estimating the carrier offset (due to Doppler, δf„ and due to LO offset, Δf,) of a plurality of signals transmitted by a plurality of BSTs and received by the MR.
71 . A method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, the method comprising the steps of: characterizing the relative group delay, Δtj - Δtj, at a RR between a signal received from the i"1 BST and a signal received from the j' BST prior to operation; removing the characterized relative group delay, Δt, - Δtj, from the estimated Time Difference Of Arrival, τ, - τ,, at the MR between a signal received from the ilh BST and a signal received from the jlh BST to yield a corrected Time Difference of Arrival; and finding the location of the MR using the corrected Time Difference of Arrival
72 The method of estimating the location of a MR in claim 71 wherein characterizing the relative group delay, Δt, - Δt, is achieved by forcing a transmitter of known location to transmit at a known time of transmission (TOT), receiving the transmitted signal at the RR; measuring the time of arrival (TOA) of the signal at the RR; estimating the time of travel as the difference between the TOA and the TOT; and finding the difference between the estimated time of travel and the known time of travel, the difference being the estimated delay of the transmitted signal through the RR.
73. The method of claim 73 further comprising repeating the steps of claim 72 a plurality of times to yield a plurality of estimated delays, then averaging the estimated delays.
74. The method of claim 72 wherein the transmitter of known location is a CT.
75. A method of estimating the location of a MR in location state 0 or location state 1 from a plurality of signals transmitted from a plurality of BSTs, wherein the transfer of the positional information from the Host to a user is achieved over the internet
76. The method of claim 75 in which the MR is calling 91 1 (local emergency telephone number) and the information is transferred to a PSAP.
77. The method of claim 75 further comprising extending the capability of a handset-based wireless location system beyond the network into which it is deployed by using an internet technology as the basis to efficiently integrate third-party database information from various locations to provide holistic location-based services on a regional, national or global basis
78 The method of claim 77 wherein the internet technology is Java based
79 A method of estimating the location of a MR from a plurality of signals tiansmitted fi om a pluiality of BSTs, compi ising using Least Squares to solve for the positional information of the MR and using TDOA observations from distant BS fs to resolve for ambiguities in an AOA solution, using AOA observations from different sectors at the same cellulai BS to resolve for ambiguities in an AOA solution; and using range observations from different BSTs to resolve for ambiguities in an AOA solution
80. The method of claim 79 further comprising quantifying the effect of geometry on position estimation, with both AOA and range observations, using the design matrix in equation (48).
81. A system for estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, the system including a host, and the system comprising; the MR being adapted to receive some or all the following information from a remote source and to apply said information to the execution of a location state: identification of base station pilot signals for the mobile to search for pilot detection threshold; duration of time for searcher to accumulate pilot energy at each delay offset; search window time length or duration; spacing of successive sampling points within the search window; frequency at which to conduct the search for pilot signals; cellular system which the mobile should acquire; the number of times the mobile should repeat all or part of the location procedure, correlation interval; locations of base stations, and the MR being adapted to transmit one or more of the following information to the Host times of arrival of pilot signals, identification of the pilot signals reported by said device, leceived signal qualities of pilot signals, pilot search lesults, and earner ficquencies of reported pilot signals
82. A system for estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, the system comprising: first computing means for characterizing the relative group delay, Δt, - Δt,, at a RR between a signal received from the i" BST and a signal received from the jlh BST prior to operation; and second computing means for removing the characterized relative group delay, Δt, - ΔtJ; from the the estimated Time Difference Of Arrival, τ, - τ|7 at the MR between a signal received from the ith BST and a signal received from the jlh BST.
83. A system for estimating the location of a MR in location state 0 or in location state 1 from a plurality of signals transmitted from a plurality of BSTs, wherein the transfer of the positional information from the Host to a user is achieved over the internet.
84. The system of claim 83 in which positional data from a Host of a MR calling 91 1 is transferred over the internet.
85. A system for estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where: Least Squares is used to solve for the positional information of the MR;
TDOA observations from distant BSTs are used to resolve for ambiguities in an AOA solution;
AOA observations from different sectors at the same cellular BS are used to resolve for ambiguities in an AOA solution; and AOA observations from different sectors at the same cellular BS are used to resolve for ambiguities in an AOA solution.
86. The system of claim 85 further comprising a first computing means for quantifying the effect of geometry on position estimation, with both AOA and range observations, using the design matrix in equation (48).
87. A method to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs; the MR receiving some or all the following information from a remote source and applies said information to the execution of a location state: identification of base station pilot signals for the mobile to search for pilot detection threshold, duration of time for searcher to accumulate pilot energy at each delay offset, seaich window time length or duration, spacing of successive sampling points within the search window, frequency at which to conduct the search for pilot signals, cellular system which the mobile should acquire, the number of times the mobile should repeat all or part of the location procedure, correlation interval, locations of base stations, and wherein some or all of the following information is transmitted from the MR to the Host times of arrival of pilot signals, identification of the pilot signals reported by said device, received signal qualities of pilot signals, pilot search results, carrier frequencies of reported pilot signals
88 A method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where misclosures and/or standardized residuals are used to flag observations that might contain a blunder
89 The method of claim 88 in which location information from each BST is weighted based on its corresponding RSSI and wherein a method selected from the group consisting of the following methods is used to solve for the positional information of the MR
General Minimum Vanance Unbiased Estimation, Best Linear Unbiased Estimation Maximum Likelihood Estimation,
Least Squares Estimation, Method of Moments, General Bayesian Estimation, Linear Bayesian Estimation and Kalman Filtering
90. The method of claim 89 wherein a method selected from the group consisting of the following methods is used to solve for the positional information of the CT:
General Minimum Variance Unbiased Estimation, Best Linear Unbiased Estimation
Maximum Likelihood Estimation, Least Squares Estimation, Method of Moments, General Bayesian Estimation, Linear Bayesian Estimation, and
Kalman Filtering.
91. The method of claim 89 wherein the flagged observations are either discarded or kept in the LS solution based on their effect on the total residuals of the LS solution and their individual redundancy number.
92. The method of claim 91 where Least Squares is repeated without the flagged observations.
93. The method of claim 92 where newly flagged observations are removed and Least Squares repeated until Least Squares converges or no more redundant observations exist.
94. The method of claim 88 which uses a method selected from the group consiting of Chaffee's method and LOCA to detect solution bifurcation.
95. The method of claim 93 where two estimates of the MR position are provided in the event that solution bifurcation does exist and there is no observational redundancy.
96. The method of claim 88 in which the hybrid TDOA positioning model in equation (42), which yields misclosures and residuals for the individual BSTs, is used.
97. The method of claim 96 wherein a method selected from the group consisting of the following methods is used to solve for the positional information of the
MR: Gcneial Minimum Variance Unbiased Estimation, Best Linear Unbiased Estimation Maximum Likelihood Estimation, Least Squares Estimation, Method of Moments,
General Bayesian Estimation, Linear Bayesian Estimation, and Kalman Filtering.
98. A method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of AOA and TDOA observations are used such that the MR can be located with as few as two AOA independent observations and one TDOA observation.
99. The method of claim 98 wherein a method selected from the group consisting of the following methods is used to solve for the positional information of the MR:
General Minimum Variance Unbiased Estimation,
Best Linear Unbiased Estimation Maximum Likelihood Estimation,
Least Squares Estimation,
Method of Moments,
General Bayesian Estimation,
Linear Bayesian Estimation, and Kalman Filtering.
100. A method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, wherein a combination of AOA and TOA observations are used such that the MR can be located with as little as one AOA observation and one TOA observation
101 The method of claim 100 wherein a method selected fiom the group consisting of the following methods is used to solve for the positional infoi mation of the MR Geneial Minimum Vanance Unbiased Estimation,
Best Lineai Unbiased Estimation Maximum Likelihood Estimation, Least Squares Estimation, Method of Moments, General Bayesian Estimation, Linear Bayesian Estimation, and
Kalman Filtering.
102. The method of claim 100 where the TOA observations are based on estimating the RTD between the BS transmission and the reception of the MR's response to the BS transmission.
103. A method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of TOA observations are used such that the MR can be located with as little as two independent TOA observations.
104. The method of claim 103 wherein a method selected from the group consisting of the following methods is used to solve for the positional information of the MR:
General Minimum Variance Unbiased Estimation, Best Linear Unbiased Estimation
Maximum Likelihood Estimation,
Least Squares Estimation,
Method of Moments,
General Bayesian Estimation, Linear Bayesian Estimation, and
Kalman Filtering.
105. The method of claim 103 where the TOA observations are based on estimating the RTD between the BST transmission and the reception of the MR's response to the BST transmission.
106. The method of claim 89 in which Least Squares is used to solve for the positional information of the MR and where the number of BSTs used in locating the MR is optimized
107. The method of claim 106 where the optimization of the number of BSs used in locating the MR is accomplished based on minimizing HDOP while maximizing the average received RSSI from all BSTs.
108. The method of claim 89 in which Least Squares is used to solve for the positional information of the MR, the method comprising of a closed-form algorithm to provide an accurate initial position to start the Least Squares iterative process.
109. The method of claim 108 where the closed-form algorithm is selected from the group consisting of LOCA, Plane Intersection, Bancroft's method, spherical interpolation, Schau and Robinson's method and Chan and Ho's method to provide an accurate initial position to start the Least Squares iterative process.
1 10. The method of claim 88 in which the mean coordinates of the participating BSs are used as the initial position for Least Squares in the event that the closed-form causes Least Squares to diverge.
1 1 1. The method of claim 88 in which the effect of geometry on AOA positioning is quantified using the design matrix in equation (45).
1 12. The method of claim 88 in which the effect of geometry on position estimation, with both AOA and TDOA observations, is quantified using the design matrix in equation (46).
1 13. A method of estimating the location of a MR from a plurality of signals transmitted from a plurality of BSTs, which quantifies the effect of geometry on speed and direction of travel estimation using the design matrix in equation (51).
1 14. The methods of estimating the location of a MR in claim 90 where the observations are selected from the group consisting of TOA, POA, and FOA
1 15. A system to estimate the location of a MR Horn a plurality of signals transmitted from a plurality of BSTs, where misclosures and/or standardized residuals are used to flag observations that might contain a blunder.
1 16. The system of claim 1 15 in which location information from each BST is weighted based on its corresponding RSSI and wherein a method selected from group consisting of the following methods is used to solve for the positional information of the MR:
General Minimum Variance Unbiased Estimation, Best Linear Unbiased Estimation
Maximum Likelihood Estimation, Least Squares Estimation, Method of Moments, General Bayesian Estimation, Linear Bayesian Estimation, and
Kalman Filtering.
1 17. The system of claim 115 which uses a method selected from the group consisting of Chaffee's method and LOCA to detect solution bifurcation.
1 18. The system of claim 1 15 in which the hybrid TDOA positioning model in equation (42), which yields misclosures and residuals for the individual BSTs, is used.
1 19. A system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of AOA and TDOA observations are used such that the MR can be located with as few as two AOA independent observations and one TDOA observation.
120. A system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of AOA and TOA observations are used such that the MR can be located with as little as one AOA observation and one TOA observation.
12 1 . A system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, where a combination of TOA observations are used such that the MR can be located with as little as two independent TOA observations.
122. The system of claim 1 16 or 121 in which Least Squares is used to solve for the positional information of the MR and where the number of BSTs used in locating the MR is optimized.
123. The system of claim 1 16 or 121 in which Least Squares is used to solve for the positional information of the MR, the method comprising of a closed-form algorithm to provide an accurate initial position to start the Least Squares iterative process.
124. The system of claim 1 15 or 1 16 in which the mean coordinates of the participating BSs are used as the initial position for Least Squares in the event that a closed-form algorithm selected from the group consisting of LOCA, Plane Intersection, Bancroft's, spherical interpolation, Schau and Robinson's method and Chan and Ho's solution causes Least Squares to diverge.
125. The system of claim 1 15 or 116 in which the effect of geometry on AOA positioning is quantified using the design matrix in equation (45).
126. A system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, which quantifies the effect of geometry on position estimation, with both AOA and TDOA observations, using the design matrix in equation (46).
127. A system to estimate the location of a MR from a plurality of signals transmitted from a plurality of BSTs, which quantifies the effect of geometry on speed and direction of travel estimation using the design matrix in equation (5 I ).
PCT/CA2000/000224 1998-10-09 2000-03-03 Methods and apparatus to position a mobile receiver using downlink signals WO2001065271A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
AU29000/00A AU2900000A (en) 2000-03-03 2000-03-03 Methods and apparatus to position a mobile receiver using downlink signals

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US16985298A 1998-10-09 1998-10-09
US09/169,730 US6204812B1 (en) 1998-10-09 1998-10-09 Methods and apparatus to position a mobile receiver using downlink signals, part II
US09/169,916 US6208297B1 (en) 1998-10-09 1998-10-09 Methods and apparatus to position a mobile receiver using downlink signals, part I
US09/169,690 US6266014B1 (en) 1998-10-09 1998-10-09 Methods and apparatus to position a mobile receiver using downlink signals part IV

Publications (1)

Publication Number Publication Date
WO2001065271A1 true WO2001065271A1 (en) 2001-09-07

Family

ID=27496873

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CA2000/000224 WO2001065271A1 (en) 1998-10-09 2000-03-03 Methods and apparatus to position a mobile receiver using downlink signals

Country Status (1)

Country Link
WO (1) WO2001065271A1 (en)

Cited By (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003052451A3 (en) * 2001-12-14 2003-11-27 Qualcomm Inc Position determination system that uses a cellular communication system
EP1394562A1 (en) * 2002-08-28 2004-03-03 Cambridge Positioning Systems Limited Improvements in radio positioning systems
WO2006007282A1 (en) * 2004-06-23 2006-01-19 Intel Corporation Effective time-of-arrival estimation algorithm for multipath environment
EP1698188A2 (en) * 2003-12-08 2006-09-06 QUALCOMM Incorporated Excess delay estimation using total received power
WO2006110475A2 (en) 2005-04-08 2006-10-19 Trueposition, Inc. Augmentation of commercial wireless location system (wls) with moving and/or airborne sensors for enhanced location accuracy and use of real-time overhead imagery for identification of wireless device locations
WO2014117273A1 (en) * 2013-02-04 2014-08-07 Takemetuit Inc. Method of processing positioning signals in positioning systems to accurately determine a true arrival time of each signal
EP2957927A1 (en) * 2014-06-19 2015-12-23 Novatel, Inc. System and method for determining location of an interfering signal source
WO2016089266A1 (en) * 2014-12-04 2016-06-09 Telefonaktiebolaget Lm Ericsson (Publ) Improving accuracy when determining positions in a wireless network
US9372255B2 (en) 2011-03-25 2016-06-21 Teledyne Instruments, Inc. Determining a position of a submersible vehicle within a body of water
CN107907857A (en) * 2017-10-25 2018-04-13 深圳大学 A kind of real-time location method and positioner based on UWB
WO2018109435A1 (en) * 2016-12-14 2018-06-21 Bae Systems Plc A control system for controlling a projectile
EP3349045A1 (en) * 2017-01-12 2018-07-18 BAE SYSTEMS plc A control system for controlling a projectile
TWI637189B (en) * 2016-10-27 2018-10-01 臺灣高等法院檢察署 Indoor area positioning system and method thereof
CN111683344A (en) * 2020-06-02 2020-09-18 南京敏智达科技有限公司 Wireless indoor positioning method based on Wi-Fi
CN111815239A (en) * 2020-07-14 2020-10-23 上海燕汐软件信息科技有限公司 Error complement order statistics method, device and computer readable storage medium
CN112188387A (en) * 2020-08-24 2021-01-05 西安交通大学 BFOA-Fang high-precision positioning method based on bounding box model
CN112470023A (en) * 2018-07-10 2021-03-09 埃尔朗根-纽伦堡弗里德里希·亚历山大大学 Positioning method and positioning system for positioning at least one object using wave-based signals
CN113078937A (en) * 2021-03-19 2021-07-06 四川航天神坤科技有限公司 Method and system for capturing data link downlink under TDD system
WO2021154983A1 (en) * 2020-01-31 2021-08-05 Trakpoint Solutions, Inc. Distributed signal processing for radiofrequency indoor localization
WO2021209117A1 (en) * 2020-04-14 2021-10-21 Pozyx N.V. Improved method and system for positioning
CN113552535A (en) * 2021-07-21 2021-10-26 西安邮电大学 Monitoring method and device for radio equipment
WO2023044314A1 (en) * 2021-09-15 2023-03-23 Parsons Corporation Far-source position determination
CN116456460A (en) * 2023-06-14 2023-07-18 河北师范大学 Filtering pretreatment method for TDOA indoor positioning
WO2024067414A1 (en) * 2022-09-30 2024-04-04 华为技术有限公司 Communication method, communication apparatus, and communication system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5663734A (en) * 1995-10-09 1997-09-02 Precision Tracking, Inc. GPS receiver and method for processing GPS signals
US5701328A (en) * 1994-12-23 1997-12-23 Stanford Telecommunications, Inc. Chirped spread spectrum positioning system
US5734674A (en) * 1990-11-28 1998-03-31 Novatel Inc. Pseudorandom-noise receiver having automatic switching between regular and anti-jamming modes
US5890068A (en) * 1996-10-03 1999-03-30 Cell-Loc Inc. Wireless location system
DE19948556A1 (en) * 1998-10-09 2000-06-29 Cell Loc Inc Mobile receiver position estimation method involves using non-linear generation of tones from signals received from several transmitters with evaluation of tone phases

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5734674A (en) * 1990-11-28 1998-03-31 Novatel Inc. Pseudorandom-noise receiver having automatic switching between regular and anti-jamming modes
US5701328A (en) * 1994-12-23 1997-12-23 Stanford Telecommunications, Inc. Chirped spread spectrum positioning system
US5663734A (en) * 1995-10-09 1997-09-02 Precision Tracking, Inc. GPS receiver and method for processing GPS signals
US5890068A (en) * 1996-10-03 1999-03-30 Cell-Loc Inc. Wireless location system
DE19948556A1 (en) * 1998-10-09 2000-06-29 Cell Loc Inc Mobile receiver position estimation method involves using non-linear generation of tones from signals received from several transmitters with evaluation of tone phases

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KLUKAS R ET AL: "FIELD TESTS OF A CELLULAR TELEPHONE POSITIONING SYSTEM", IEEE VEHICULAR TECHNOLOGY CONFERENCE,US,NEW YORK, IEEE, vol. CONF. 47, 4 May 1997 (1997-05-04), pages 470 - 474, XP000736651, ISBN: 0-7803-3660-7 *

Cited By (47)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003052451A3 (en) * 2001-12-14 2003-11-27 Qualcomm Inc Position determination system that uses a cellular communication system
US7877100B2 (en) 2001-12-14 2011-01-25 Qualcomm Incorporated Position determination system that uses A cellular communication system
CN100409028C (en) * 2001-12-14 2008-08-06 高通股份有限公司 Position determination system that uses a cellular communication system
US7315745B2 (en) 2002-08-28 2008-01-01 Cambridge Positioning Systems Ltd. Radio positioning systems
EP1394562A1 (en) * 2002-08-28 2004-03-03 Cambridge Positioning Systems Limited Improvements in radio positioning systems
WO2004021033A1 (en) * 2002-08-28 2004-03-11 Cambridge Positioning Systems Limited Improvements in radio positioning systems
EP1698188A2 (en) * 2003-12-08 2006-09-06 QUALCOMM Incorporated Excess delay estimation using total received power
EP1698188A4 (en) * 2003-12-08 2011-06-15 Qualcomm Inc Excess delay estimation using total received power
US8725161B2 (en) 2003-12-08 2014-05-13 Qualcomm Incorporated Excess delay estimation using total received power
US7221913B2 (en) 2004-06-23 2007-05-22 Intel Corporation Effective time-of-arrival estimation algorithm for multipath environment
EP3512103A1 (en) * 2004-06-23 2019-07-17 Guangdong OPPO Mobile Telecommunications Corp., Ltd. Effective time-of-arrival estimation algorithm for multipath environment
EP3512104A1 (en) * 2004-06-23 2019-07-17 Guangdong OPPO Mobile Telecommunications Corp., Ltd. Effective time-of-arrival estimation algorithm for multipath environment
WO2006007282A1 (en) * 2004-06-23 2006-01-19 Intel Corporation Effective time-of-arrival estimation algorithm for multipath environment
WO2006110475A2 (en) 2005-04-08 2006-10-19 Trueposition, Inc. Augmentation of commercial wireless location system (wls) with moving and/or airborne sensors for enhanced location accuracy and use of real-time overhead imagery for identification of wireless device locations
US7427952B2 (en) 2005-04-08 2008-09-23 Trueposition, Inc. Augmentation of commercial wireless location system (WLS) with moving and/or airborne sensors for enhanced location accuracy and use of real-time overhead imagery for identification of wireless device locations
WO2006110475A3 (en) * 2005-04-08 2009-05-22 Trueposition Inc Augmentation of commercial wireless location system (wls) with moving and/or airborne sensors for enhanced location accuracy and use of real-time overhead imagery for identification of wireless device locations
US9372255B2 (en) 2011-03-25 2016-06-21 Teledyne Instruments, Inc. Determining a position of a submersible vehicle within a body of water
WO2014117273A1 (en) * 2013-02-04 2014-08-07 Takemetuit Inc. Method of processing positioning signals in positioning systems to accurately determine a true arrival time of each signal
GB2525110B (en) * 2013-02-04 2016-09-21 Takemetuit Inc Method of processing positioning signals in positioning systems to accurately determine a true arrival time of each signal
US9678195B2 (en) 2013-02-04 2017-06-13 Takemetuit Inc. Method of processing positioning signals in positioning systems to accurately determine a true arrival time of each signal
GB2525110A (en) * 2013-02-04 2015-10-14 Takemetuit Inc Method of processing positioning signals in positioning systems to accurately determine a true arrival time of each signal
US9766343B2 (en) 2014-06-19 2017-09-19 Novatel Inc. System and method for determining location of an interfering signal source
EP2957927A1 (en) * 2014-06-19 2015-12-23 Novatel, Inc. System and method for determining location of an interfering signal source
WO2016089266A1 (en) * 2014-12-04 2016-06-09 Telefonaktiebolaget Lm Ericsson (Publ) Improving accuracy when determining positions in a wireless network
US10222451B2 (en) 2014-12-04 2019-03-05 Telefonaktiebolaget Lm Ericsson (Publ) Improving accuracy when determining positions in a wireless network
TWI637189B (en) * 2016-10-27 2018-10-01 臺灣高等法院檢察署 Indoor area positioning system and method thereof
WO2018109435A1 (en) * 2016-12-14 2018-06-21 Bae Systems Plc A control system for controlling a projectile
AU2017374695B2 (en) * 2016-12-14 2022-08-11 Bae Systems Plc A control system for controlling a projectile
US11774544B2 (en) 2016-12-14 2023-10-03 Bae Systems Plc Control system for controlling a projectile
EP3349045A1 (en) * 2017-01-12 2018-07-18 BAE SYSTEMS plc A control system for controlling a projectile
CN107907857A (en) * 2017-10-25 2018-04-13 深圳大学 A kind of real-time location method and positioner based on UWB
CN107907857B (en) * 2017-10-25 2021-08-17 深圳大学 UWB-based real-time positioning method and positioning device
CN112470023B (en) * 2018-07-10 2024-01-16 埃尔朗根-纽伦堡弗里德里希·亚历山大大学 Positioning method and positioning system for positioning at least one object by using wave-based signals
CN112470023A (en) * 2018-07-10 2021-03-09 埃尔朗根-纽伦堡弗里德里希·亚历山大大学 Positioning method and positioning system for positioning at least one object using wave-based signals
WO2021154983A1 (en) * 2020-01-31 2021-08-05 Trakpoint Solutions, Inc. Distributed signal processing for radiofrequency indoor localization
WO2021209117A1 (en) * 2020-04-14 2021-10-21 Pozyx N.V. Improved method and system for positioning
CN111683344A (en) * 2020-06-02 2020-09-18 南京敏智达科技有限公司 Wireless indoor positioning method based on Wi-Fi
CN111815239A (en) * 2020-07-14 2020-10-23 上海燕汐软件信息科技有限公司 Error complement order statistics method, device and computer readable storage medium
CN111815239B (en) * 2020-07-14 2024-04-23 上海燕汐软件信息科技有限公司 Error complement order statistics method, device and computer readable storage medium
CN112188387A (en) * 2020-08-24 2021-01-05 西安交通大学 BFOA-Fang high-precision positioning method based on bounding box model
CN113078937A (en) * 2021-03-19 2021-07-06 四川航天神坤科技有限公司 Method and system for capturing data link downlink under TDD system
CN113552535A (en) * 2021-07-21 2021-10-26 西安邮电大学 Monitoring method and device for radio equipment
CN113552535B (en) * 2021-07-21 2024-04-09 西安邮电大学 Radio equipment monitoring method and device
WO2023044314A1 (en) * 2021-09-15 2023-03-23 Parsons Corporation Far-source position determination
WO2024067414A1 (en) * 2022-09-30 2024-04-04 华为技术有限公司 Communication method, communication apparatus, and communication system
CN116456460A (en) * 2023-06-14 2023-07-18 河北师范大学 Filtering pretreatment method for TDOA indoor positioning
CN116456460B (en) * 2023-06-14 2023-09-01 河北师范大学 Filtering pretreatment method for TDOA indoor positioning

Similar Documents

Publication Publication Date Title
US6266014B1 (en) Methods and apparatus to position a mobile receiver using downlink signals part IV
US6208297B1 (en) Methods and apparatus to position a mobile receiver using downlink signals, part I
US6204812B1 (en) Methods and apparatus to position a mobile receiver using downlink signals, part II
US6330452B1 (en) Network-based wireless location system to position AMPs (FDMA) cellular telephones, part I
WO2001065271A1 (en) Methods and apparatus to position a mobile receiver using downlink signals
US5890068A (en) Wireless location system
WO2000069198A1 (en) Wireless location system
US7593738B2 (en) GPS synchronization for wireless communications stations
Reed et al. An overview of the challenges and progress in meeting the E-911 requirement for location service
US7453961B1 (en) Methods and apparatus for detection of signal timing
US6765531B2 (en) System and method for interference cancellation in a location calculation, for use in a wireless location system
KR101083446B1 (en) Augmentation of commercial wireless location system(wls) with moving and/or airborne sensors for enhanced location accuracy and use of real-time overhead imagery for identification of wireless device locations
US6891500B2 (en) Method and apparatus for geolocating a wireless communications device
US6111538A (en) Cellular positioning system (CPS)
CA2701108C (en) Detection of time of arrival of cdma signals in a wireless location system
CA2821652C (en) Tdoa/gps hybrid wireless location system
US20120165037A1 (en) Robust Downlink Frame Synchronization Schemes in CDMA Wireless Networks for Geo-Location
KR20030045109A (en) Gps satellite signal acquisition assistance system and method in a wireless communications network
WO2001058195A1 (en) A network-based wireless location system to positon amps (fdma) cellular telephones
Klukas Cellular telephone positioning using GPS time synchronization
Klukas A superresolution based cellular positioning system using GPS time synchronization
Ma Techniques to improve ground-based wireless location performance using a cellular telephone network
AU2003200955B2 (en) Passive position determination using two low-earth orbit satellites
Mensing et al. GNSS positioning in critical scenarios: hybrid data fusion with communications signals
WO1996026579A1 (en) Novel cellular positioning system (cps)

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

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

AL Designated countries for regional patents

Kind code of ref document: A1

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

121 Ep: the epo has been informed by wipo that ep was designated in this application
REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

122 Ep: pct application non-entry in european phase