CA2513154A1 - A method and system for locating a mobile radio receiver in a radio system with multiple transmitters - Google Patents
A method and system for locating a mobile radio receiver in a radio system with multiple transmitters Download PDFInfo
- Publication number
- CA2513154A1 CA2513154A1 CA002513154A CA2513154A CA2513154A1 CA 2513154 A1 CA2513154 A1 CA 2513154A1 CA 002513154 A CA002513154 A CA 002513154A CA 2513154 A CA2513154 A CA 2513154A CA 2513154 A1 CA2513154 A1 CA 2513154A1
- Authority
- CA
- Canada
- Prior art keywords
- receiver
- list
- function
- signals
- transmitter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W8/00—Network data management
- H04W8/02—Processing of mobility data, e.g. registration information at HLR [Home Location Register] or VLR [Visitor Location Register]; Transfer of mobility data, e.g. between HLR, VLR or external networks
- H04W8/08—Mobility data transfer
-
- 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/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
-
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W8/00—Network data management
- H04W8/02—Processing of mobility data, e.g. registration information at HLR [Home Location Register] or VLR [Visitor Location Register]; Transfer of mobility data, e.g. between HLR, VLR or external networks
Abstract
A method of determining the location of a mobile receiver is provided. The method works within a radio system in one or more dimensions, the radio system includes a network of transmitters the positions of which are known. The method comprises the steps of: compiling a list of the signals detected by the receiver; sending said list to a calculation device; constructing a position function which relates to the probability that the receiver is at a given location, which position function has components corresponding to each member of the list, each component including a predetermined probability function that depends on mobile receiver position; and evaluating said position function to determine a position that corresponds to the highest probability of location, and defining that position to be the location of the mobile receiver.
Description
A METHOD AND SYSTEM FOR LOCATING A MOBILE RADIO RECEIVER IN A
RADIO SYSTEM WITH MULTIPLE TRANSMITTERS
Technical field This invention relates to multiple transmitter radio systems and in particular to locating the position of a mobile radio receiver in a system where there is a plurality of transmitting stations each located at a different geographic location in order to provide the. area of interest with adequate radio coverage.
10. Description of Related Art Conventional systems for locating radio receivers . using a plurality of transmitters use different techniques for measuring the characteristics of the received signal to find the position of the receiver. These' can be, characterised as time referencing, phase referencing, or spatial referencing, or techniques which combine these three types.
In time referencing systems, for example the Loran C system set up by the US
coastguard, the radio transmissions take the form of radio pulses whose time of arrival at the mobile receiver can be accurately measured. The receiver is of a special type that can determine time differences between the reception of pulses from each of a pair of transmitters in the system. This information defines a curve connecting points of equal time difference on which the receiver must lie. By similar measurements of multiple transmitter pairs, further intersecting curves enable the receiver's position to be determined as the point of intersection of the curves which best satisfies all the measured time differences.
An extension of this system is the Global Positioning System (GPS) which uses the reception of signals from a number of satellite transmitters that are synchronised in time and which also transmit their positions to the receiver.
The receiver measures the time of arrival of a defined point in the received data, stream from each satellite with reference to the receiver's internal time source.
Each 30. measurement defines a sphere on which the receiver must lie and the radius of which is fixed by the time of flight of the radio waves from the satellite. Similar measurements using different satellites in different positions create multiple intersecting. spheres. The receiver's position can be determined, by solving the simultaneous equations for the distances to the satellites and the receiver's time offset, as the point of intersection of the spheres which best ,satisfies all of the measured time delays and the receiver clock offset.
Systems that use phase measurements include the Decca Navigation system and the Omega system. The i~nobile receiver measures the relative phase of the received . signal from special transmitters which have phase synchronised transmissions, either with respect, to each other or with respect to an internal phase reference. By knowing the positions of the transmitters and the predicted phase received at all points in the covered area, the position of the receiver consistent with a particular set of phase measurements may be determined.
All of the above systems require dedicated and expensive transmitter systems .
with specialised receivers: ~ .
Another prior art method, for example US-5859612-A uses the directional nature of the transmitted signals. This has been applied in a method which uses a beam of identifiable modulation, which synchronously sweeps a defined' area from a known start,position and start time. By knowing when the beam started.its sweep and the reception time of the defined modulation, the receiver's. position can be determined to be on a radial line of certain direction from the source transmitter.
Similar , receptions from multiple sources can fix the receiver's position in two dimensions by establishing the point that best satisfies all radial line bearings simultaneously.
Spatially-fixed beams of radio transmissions that have narrow transmission patterns from fixed antennas can also be used to determine the position of a receiver.
By knowing which directional pattern of the overall transmitted pattern is received, then the position of the receiver within a sector centred on the transmitting station can be determined. The reception of narrow beams from multiple overlapping transmitter patterns can similarly fix the receiver in space by searching for the point that best satisfies all the beam directions to the transmitting stations.
The above system also requires dedicated transmitters, but these are usually.
localised in particular areas of interest. The accuracy of the system diminishes rapidly as the beam-width, and therefore the exact bearing of the receiver, increases.
The receivers are generally less specialised than time- or phase-measuring receivers and can share their positioning ~ function with other functions such as weather report reception.
The above last method that uses the directional patterns of transmissions from fixed base transmitter stations can be applied to digital radio cellular systems. In such systems the area to be covered by the communication service is serviced by a plurality of fixed transmitting and receiving stations, usually co-located.
Tlie area covered by a single fixed transmitter site on one of its directional antenna patterns is .
RADIO SYSTEM WITH MULTIPLE TRANSMITTERS
Technical field This invention relates to multiple transmitter radio systems and in particular to locating the position of a mobile radio receiver in a system where there is a plurality of transmitting stations each located at a different geographic location in order to provide the. area of interest with adequate radio coverage.
10. Description of Related Art Conventional systems for locating radio receivers . using a plurality of transmitters use different techniques for measuring the characteristics of the received signal to find the position of the receiver. These' can be, characterised as time referencing, phase referencing, or spatial referencing, or techniques which combine these three types.
In time referencing systems, for example the Loran C system set up by the US
coastguard, the radio transmissions take the form of radio pulses whose time of arrival at the mobile receiver can be accurately measured. The receiver is of a special type that can determine time differences between the reception of pulses from each of a pair of transmitters in the system. This information defines a curve connecting points of equal time difference on which the receiver must lie. By similar measurements of multiple transmitter pairs, further intersecting curves enable the receiver's position to be determined as the point of intersection of the curves which best satisfies all the measured time differences.
An extension of this system is the Global Positioning System (GPS) which uses the reception of signals from a number of satellite transmitters that are synchronised in time and which also transmit their positions to the receiver.
The receiver measures the time of arrival of a defined point in the received data, stream from each satellite with reference to the receiver's internal time source.
Each 30. measurement defines a sphere on which the receiver must lie and the radius of which is fixed by the time of flight of the radio waves from the satellite. Similar measurements using different satellites in different positions create multiple intersecting. spheres. The receiver's position can be determined, by solving the simultaneous equations for the distances to the satellites and the receiver's time offset, as the point of intersection of the spheres which best ,satisfies all of the measured time delays and the receiver clock offset.
Systems that use phase measurements include the Decca Navigation system and the Omega system. The i~nobile receiver measures the relative phase of the received . signal from special transmitters which have phase synchronised transmissions, either with respect, to each other or with respect to an internal phase reference. By knowing the positions of the transmitters and the predicted phase received at all points in the covered area, the position of the receiver consistent with a particular set of phase measurements may be determined.
All of the above systems require dedicated and expensive transmitter systems .
with specialised receivers: ~ .
Another prior art method, for example US-5859612-A uses the directional nature of the transmitted signals. This has been applied in a method which uses a beam of identifiable modulation, which synchronously sweeps a defined' area from a known start,position and start time. By knowing when the beam started.its sweep and the reception time of the defined modulation, the receiver's. position can be determined to be on a radial line of certain direction from the source transmitter.
Similar , receptions from multiple sources can fix the receiver's position in two dimensions by establishing the point that best satisfies all radial line bearings simultaneously.
Spatially-fixed beams of radio transmissions that have narrow transmission patterns from fixed antennas can also be used to determine the position of a receiver.
By knowing which directional pattern of the overall transmitted pattern is received, then the position of the receiver within a sector centred on the transmitting station can be determined. The reception of narrow beams from multiple overlapping transmitter patterns can similarly fix the receiver in space by searching for the point that best satisfies all the beam directions to the transmitting stations.
The above system also requires dedicated transmitters, but these are usually.
localised in particular areas of interest. The accuracy of the system diminishes rapidly as the beam-width, and therefore the exact bearing of the receiver, increases.
The receivers are generally less specialised than time- or phase-measuring receivers and can share their positioning ~ function with other functions such as weather report reception.
The above last method that uses the directional patterns of transmissions from fixed base transmitter stations can be applied to digital radio cellular systems. In such systems the area to be covered by the communication service is serviced by a plurality of fixed transmitting and receiving stations, usually co-located.
Tlie area covered by a single fixed transmitter site on one of its directional antenna patterns is .
-2-called a sector or a cell. The combination of all of the cell positions has been designed to enable a large number of mobile transmitters and receivers in a .
large geographic area to be in continuous communication with at least one of these fixed .
transmitter-receivers (the serving base station) as the mobile station moves in any arbitrary direction. The fixed transmitter-receiver in this system is usually called a base station and the mobile receiver may receive the transmissions of neighbouring base stations because their antenna patterns overlap with that of~ the serving base station. In any particular case, the base station that the mobile receiver uses to .
connect to the underlying service network changes as the. receiver moves with respect to the fixed network. This process is called hand-off and it ensures that the mobile receiver's communications traffic is passed smoothly from one base station to another so that there is no lapse in service.
To enable the handofF process the cellular mobile receiver receives transmissions from a plurality of base stations. Related art shows how these transmissions made for telecommunication purposes may be also be used by the mobile receiver, or in .a processing device in communication with the receiver, to calculate the receiver's position. These prior art methods will now be described along with their disadvantages.
As described in US-6108553-A the mobile receiver. may measure the relative time delay between the signals received from two transmitted by time synchronised base stations. These measurements, when taken from several different pairs of base station transmissions, form an intersecting set of hyperbolic curves. The point that best satisfies the intersection of all of these curves predicts the position of the receiver.
The mobile receiver may also measure the reception of signals from each base station with reference to its own internal clock, see for example WO
97/11384 or WO 99/21028. This reception time may then be combined with measurements of fihe same transmission made at other external measurement receivers in the network.
The combined measurement information can then be combined in the handset or in other generalised computing means to provide a set of intersecting circles centred on the transmitters. The position of the handset is calculated to be at the point of intersection of the circles.
The preceding techniques require additional equipment to be installed and maintained at the transmitters in the network to enable timing measurements, and the handset must also be modified from its normal communications function to ~
make additional timing measurements.
large geographic area to be in continuous communication with at least one of these fixed .
transmitter-receivers (the serving base station) as the mobile station moves in any arbitrary direction. The fixed transmitter-receiver in this system is usually called a base station and the mobile receiver may receive the transmissions of neighbouring base stations because their antenna patterns overlap with that of~ the serving base station. In any particular case, the base station that the mobile receiver uses to .
connect to the underlying service network changes as the. receiver moves with respect to the fixed network. This process is called hand-off and it ensures that the mobile receiver's communications traffic is passed smoothly from one base station to another so that there is no lapse in service.
To enable the handofF process the cellular mobile receiver receives transmissions from a plurality of base stations. Related art shows how these transmissions made for telecommunication purposes may be also be used by the mobile receiver, or in .a processing device in communication with the receiver, to calculate the receiver's position. These prior art methods will now be described along with their disadvantages.
As described in US-6108553-A the mobile receiver. may measure the relative time delay between the signals received from two transmitted by time synchronised base stations. These measurements, when taken from several different pairs of base station transmissions, form an intersecting set of hyperbolic curves. The point that best satisfies the intersection of all of these curves predicts the position of the receiver.
The mobile receiver may also measure the reception of signals from each base station with reference to its own internal clock, see for example WO
97/11384 or WO 99/21028. This reception time may then be combined with measurements of fihe same transmission made at other external measurement receivers in the network.
The combined measurement information can then be combined in the handset or in other generalised computing means to provide a set of intersecting circles centred on the transmitters. The position of the handset is calculated to be at the point of intersection of the circles.
The preceding techniques require additional equipment to be installed and maintained at the transmitters in the network to enable timing measurements, and the handset must also be modified from its normal communications function to ~
make additional timing measurements.
-3-The prior art described in US-5293642-A also describes position estimating techniques using the spatial characteristics of radio waves applied to cellular systems. .In these,systems the base station receiver measures the radio propagation parameters of the signals transmitted by the mobile station transmitter. These may include signal strength and travel time. The parameters are then used, together with knowledge of the base station antenna pattern, to derive a probability density function (centred on each base station) of a mobile receiver giving the .predicted signal strength at any position. By combining this probability density function with similar ' probability density functions centred on other base stations, a maximum in the joint probability ~ density functions is derived which can be interpreted as being the .
maximum-likelihood position of the mobile station.
In other prior art, for example US-5613205-A, a mobile receiver measures the relative signal strengths of the control signals received from a plurality of base stations. The system then calculates the distance between the receiver and each base station as a function of the base station transmitted power, received power and the signal attenuation factors in the radio path from each base station to the mobile receiver.
Simpler systems for positioning mobile receivers have also been described, One of these takes the centroid of the serving cell as the handset's position.
Another measures the signal strength of the neighbouring cells and take the centroid of that with the biggest signal level as the best indication of the position of the receiver as in
maximum-likelihood position of the mobile station.
In other prior art, for example US-5613205-A, a mobile receiver measures the relative signal strengths of the control signals received from a plurality of base stations. The system then calculates the distance between the receiver and each base station as a function of the base station transmitted power, received power and the signal attenuation factors in the radio path from each base station to the mobile receiver.
Simpler systems for positioning mobile receivers have also been described, One of these takes the centroid of the serving cell as the handset's position.
Another measures the signal strength of the neighbouring cells and take the centroid of that with the biggest signal level as the best indication of the position of the receiver as in
4.
These particular prior art systems and methods are unreliable for estimating position since the centre of the serving cell may not be the closest approximation to the position of the receiver, and the cell is of indefinite shape and subject.
to anomalous propagation conditions. The signal level measurement of a single base station signal might also be erroneous for similar reasons. These techniques suffer from the highly unpredictable nature of the radio propagation path and the inaccuracies surrounding the mathematical models. In particular, estimated power levels do not usually conform to simple Gaussian_ distributions, especially when subject to multi-path fading. Of course, the actual power received at a number of points in the network may be determined a priori by employing a mobile test set which is calibrated for position and signal strength. However, such measurements are costly and time consuming to make, and can make no allowance for the changeable nature ~of real networks over time between measurements, The present invention is designed to overcome many of the shortcomings of the systems iri the prior art outlined above. In particular, the position of a mobile receiver in a communications network can be established with. relatively high . accuracy using only measurements made by the receiver without the need to provide additional equipment. in the network nor to carry out periodic calibrations of signal reception.
Summary of the invention According to the present invention, there is provided a method of determining . 10 the location of a mobile receiver in one or more dii~nensions in a radio system which includes a network of transmitters the positions of which are known, the method comprising the steps of: , , . .
a) compiling a list of the signals detected by the receiver;
b) sending said list to a calculation device;
~ c) constructing a position function. which relates to the probability that the receiver is at a given location, which position function has components corresponding to each member of the list, each component including a predetermined probability function that depends on mobile receiver position;
and d) evaluating said position function to determine a position that corresponds ~ to the highest probability of location, and defining that position to be the location of the mobile receiver.
The detection of signals by the receiver may include the reception of signals transmitted by a particular transmitter of the network and the decoding of information .
from the signals. It may also include identification of a particular signature, such as a synchronisation pattern or a pilot code, in signals received above a given threshold power, or it may simply be the reception of signal power on a ,particular frequency above a given threshold. . ~ .
The list of signals detected by the receiver may, for example, include decoded information such as the identifications of the transmitters, or it may consist of a list of frequencies on which signals above a given threshold power were received, or it may simply consist of a list of references to another list, known to the calculation device, which contains information about the transmitters The calculation device may be attached to the mobile receiver, or it may be at any location able to receive the list sent by the receiver.
These particular prior art systems and methods are unreliable for estimating position since the centre of the serving cell may not be the closest approximation to the position of the receiver, and the cell is of indefinite shape and subject.
to anomalous propagation conditions. The signal level measurement of a single base station signal might also be erroneous for similar reasons. These techniques suffer from the highly unpredictable nature of the radio propagation path and the inaccuracies surrounding the mathematical models. In particular, estimated power levels do not usually conform to simple Gaussian_ distributions, especially when subject to multi-path fading. Of course, the actual power received at a number of points in the network may be determined a priori by employing a mobile test set which is calibrated for position and signal strength. However, such measurements are costly and time consuming to make, and can make no allowance for the changeable nature ~of real networks over time between measurements, The present invention is designed to overcome many of the shortcomings of the systems iri the prior art outlined above. In particular, the position of a mobile receiver in a communications network can be established with. relatively high . accuracy using only measurements made by the receiver without the need to provide additional equipment. in the network nor to carry out periodic calibrations of signal reception.
Summary of the invention According to the present invention, there is provided a method of determining . 10 the location of a mobile receiver in one or more dii~nensions in a radio system which includes a network of transmitters the positions of which are known, the method comprising the steps of: , , . .
a) compiling a list of the signals detected by the receiver;
b) sending said list to a calculation device;
~ c) constructing a position function. which relates to the probability that the receiver is at a given location, which position function has components corresponding to each member of the list, each component including a predetermined probability function that depends on mobile receiver position;
and d) evaluating said position function to determine a position that corresponds ~ to the highest probability of location, and defining that position to be the location of the mobile receiver.
The detection of signals by the receiver may include the reception of signals transmitted by a particular transmitter of the network and the decoding of information .
from the signals. It may also include identification of a particular signature, such as a synchronisation pattern or a pilot code, in signals received above a given threshold power, or it may simply be the reception of signal power on a ,particular frequency above a given threshold. . ~ .
The list of signals detected by the receiver may, for example, include decoded information such as the identifications of the transmitters, or it may consist of a list of frequencies on which signals above a given threshold power were received, or it may simply consist of a list of references to another list, known to the calculation device, which contains information about the transmitters The calculation device may be attached to the mobile receiver, or it may be at any location able to receive the list sent by the receiver.
-5-The predetermined probability function may be based on a propagation model for the radio waves between each transmitter and the receiver assumed to beat fihat position, said model taking into account characteristics such as the transmitter power, transmitter antenna pattern, travel distance, and the variable transmission characteristics (such ~as fading) of signal propagation. Preferably, the predetermined probability function is proportional to the probability of detection of the' signals from that transmitter, and hence possible membership of the list. Alternatively, the said predetermined probability function may be based on a heuristic derived from detections obtained at calibration points. For example, it may be found in practice that the probability of detection of the signals from a transmitter at a particular point in the area covered by that transmitter can be approximated by a two-dimensional Gaussian function centred near the middle of the area. Use of such a heuristic can give very good results since it is based on real measurements rather than a theoretical model.
The list may be augmented, or a second list constructed,, to include information from signals that were not reported by the receiver, but which could have been detected had the signal propagation characteristics been favourable. In these cases, the, position function may be modified to take account of the probability of non-appearance in the list.
The position function may also be affected by the reliability~of the detection of each of the signals by the receiver. An estimate of the reliability could, for example, be proportional to the difference between the received power level and the detection threshold, or it may be a measure of the number of times a given signal is detected in a repeated series of measurements. The reliability measure could be transmitted in the messages communicated from the mobile receiver, either as part of the list or in a separate message.
In some applications, it may be an advantage for the receiver to sort the list before sending it. For example, the received signal power or the. reliability measure described previously could be used for sorting.
Digital communications 'systems often measure the round-trip time of flight of signals from the serving transmitter to the mobile receiver and back again.
.The measurement is used to advance the timing of the signals transmitted by the handset transmitter so that they arrive at a known time at the base station. Where it is available, this so-called timing advance value (TA) can preferably be included in the position.function, for example by assigning an additional probability that the mobile receiver is at the given location given the measured value,
The list may be augmented, or a second list constructed,, to include information from signals that were not reported by the receiver, but which could have been detected had the signal propagation characteristics been favourable. In these cases, the, position function may be modified to take account of the probability of non-appearance in the list.
The position function may also be affected by the reliability~of the detection of each of the signals by the receiver. An estimate of the reliability could, for example, be proportional to the difference between the received power level and the detection threshold, or it may be a measure of the number of times a given signal is detected in a repeated series of measurements. The reliability measure could be transmitted in the messages communicated from the mobile receiver, either as part of the list or in a separate message.
In some applications, it may be an advantage for the receiver to sort the list before sending it. For example, the received signal power or the. reliability measure described previously could be used for sorting.
Digital communications 'systems often measure the round-trip time of flight of signals from the serving transmitter to the mobile receiver and back again.
.The measurement is used to advance the timing of the signals transmitted by the handset transmitter so that they arrive at a known time at the base station. Where it is available, this so-called timing advance value (TA) can preferably be included in the position.function, for example by assigning an additional probability that the mobile receiver is at the given location given the measured value,
-6-The calculation of the position of the receiver using the position function to find the position corresponding to the highest probability can be carried out in a number of ways. Preferably, the position function is a 'cost function'.which describes the logarithm of the reciprocal of the probability as a function of the location of the receiver and membership of the list: This cost function has the lowest numerical value corresponding to the location which is most likely to give the measured list of, detected signals, including, if available; the list of signals not detected and the .
measured TA and any.reported ordering of the list. This minimum value can be found by,any one of a number of well-known general mathematical minimisation techniques including (a) analytical.analysis in which the position function can be manipulated to give the solution directly, (b) iterative methods which generally seek:~to move a,current estimate of the solution 'down the slope' in such a direction as to move closer to the minimum point, and (c) 'grid methods' in which the position function is evaluated at every point of a grid of positions with the position corresponding to the lowest value being taken as the result. These methods are described in many textbooks, for example; Numerical Recipes in C++, Second Edition by William H. Press' et al, Cambridge University Press, Chapter 10. pages 401-4.29.
The current invention therefore provides a method of obtaining the position in one or more dimensions of a mobile receiver by estimating its location in, for example, a geographic area that has a number of receivable and identifiable radio transmissions established within it. The method does not rely on a. detailed knowledge of the transmission modulation type or content or time delay measurements. It can therefore be applied to any radio system of any modulation type without knowledge of the system timing.
The mobile receiver at a given position receives signals from one or more of the transmitters in the network. Those that are detected can be called 'visible' and can be added to the list ~ maintained of visible stations. The position of the receiver affects which transmitting stations are visible to it. A systerri according to the present invention may also maintain a list of those transmitters that cannot be received and are called 'not visible', but are known to exist and be active, at the time of the position measurement. The receiver's position affects which transmitting base stations are not visible.
At any given location it is possible to predict theoretically which transmitters are visible and which are not visible, depending upon factors such as the distance from the transmitting station, the transmitted power of the transmitter, and the direction of the transmitter radiation pattern. According to one implementation of the
measured TA and any.reported ordering of the list. This minimum value can be found by,any one of a number of well-known general mathematical minimisation techniques including (a) analytical.analysis in which the position function can be manipulated to give the solution directly, (b) iterative methods which generally seek:~to move a,current estimate of the solution 'down the slope' in such a direction as to move closer to the minimum point, and (c) 'grid methods' in which the position function is evaluated at every point of a grid of positions with the position corresponding to the lowest value being taken as the result. These methods are described in many textbooks, for example; Numerical Recipes in C++, Second Edition by William H. Press' et al, Cambridge University Press, Chapter 10. pages 401-4.29.
The current invention therefore provides a method of obtaining the position in one or more dimensions of a mobile receiver by estimating its location in, for example, a geographic area that has a number of receivable and identifiable radio transmissions established within it. The method does not rely on a. detailed knowledge of the transmission modulation type or content or time delay measurements. It can therefore be applied to any radio system of any modulation type without knowledge of the system timing.
The mobile receiver at a given position receives signals from one or more of the transmitters in the network. Those that are detected can be called 'visible' and can be added to the list ~ maintained of visible stations. The position of the receiver affects which transmitting stations are visible to it. A systerri according to the present invention may also maintain a list of those transmitters that cannot be received and are called 'not visible', but are known to exist and be active, at the time of the position measurement. The receiver's position affects which transmitting base stations are not visible.
At any given location it is possible to predict theoretically which transmitters are visible and which are not visible, depending upon factors such as the distance from the transmitting station, the transmitted power of the transmitter, and the direction of the transmitter radiation pattern. According to one implementation of the
-7-invention, each transmitter in the visible and not-visible lists is assigned a 'cost' function. This cost function describes the combination of the transmitter power, transmitter radiation pattern, distance and bearing of the transmitter from the receiver, . radio propagation characteristics, and (where available) other characteristics such as timing advance or local topological data, which best models the probability ~
of detection of the signals by the mobile receiver. For mathematical convenience, the cost function combining all the characteristics affecting the visibility of the transmitters is ri~ade proportional to the negative of the logarithm of the' probability of detecting.
their signals (i.e. the logarithm of the reciprocal of the probability). Hence taking .one 10' transmitter in isolation, the cost function describing its visibility has a low value when the receiver is close to the transmitter position.and directly in line with the transmitter antenna's forward lobe as shown, for example, in Figure ~1. The cost function gradually increases as the receiver moves directly away from the base station, or moves tangentially at an angle from the direction of the forward lobe at fixed distance from the transmitter. The corresponding cost functions describing each transmitter are then combined as a sum to give an overall cost furiction that describes the combined visibility of the transmitters as shown for example in .Figure 2. The minimum of this overall cost function is the position of the receiver according to the invention.
According to the .present invention there is provided a mobile receiver device .
capable of detecting signals from a plurality of transmitters and arranged to compute a list of signals received by the mobile receiver and to send said list, in use, to a calculation device.
Furthermore, according to the present invention there is provided a calculation device for use in location determination of a mobile receiver, the calculation device being arranged to:
a) receive a list of the signals detected by a receiver, b) construct a position function which relates to the probability that .the receiver is at a given ~ location, which position , has components corresponding to each member of the list, each component including a predetermined probability function that depends on the mobile ,receiver position; and ' .
c) evaluate said position function to determine a position that corresponds to the highest probability of location, and defining that. position to be the location of the mobile receiver. w 35According to the present invention there is provided a system comprising: a mobile receiver device and a calculation device.
_g_ The calculation device.may be attached to the mobile receiver, or it maybe remote from, but in communication with, the mobile receiver.
Brief description of the drawings ~ Figure 1 shows an example plot of the position function ('Cost') associated with a single transmitter.
Figure 2 shows an example plot of the position function ('Cost') associated with multiple transmitters.
Figure 3 illustrates an example multiple transmitter system.
Figure 4 shows an example list of detected signals sent by the handset, and examples of lists maintained at the location calculation unit .of visible and not-visible .
signals. .
Figure 5 illustrates the variation of received power with distance from the transmitter.
Figure 6 shows an example fixed transmitter and a plot of the position function component associated .with timing ,advance.
Figure 7, illustrates a calculation grid. ..
Detailed descriptions of particular embodiments A mobile receiver positioning system has been implemented according to the current invention using an existing GSM mobile telephony network (see Figure 3).
The system 300 comprises a plurality of transceiver base stations 301, 302, 303, 304, which act as the transmitters described above, a handset which acts as the mobile receiver 310, a location calculation unit (LCU) 320, and a management database (MDB) 330. The base stations provide the signal transmissions which are detected' by the handset 310 and processed in the LCU 320 using the mathematical algorithms described below.
The handset 310 detects signals from a riumber of the nearby base stations 301-304, compiles a list of those it detects, and sends the list via a message over link 340 to the serving base station 304: The format and content of the message is shown in Figure 4 at 401. It is defined by the GSM standard and includes fields for the Broadcast Control Channel; BCCH, the Base Station Identity Code, BSIC, the Mobile Network Code, MNC, the Mobile Country Code, MCC, the Cell Identifier, CI, the Local . Area Code, LAC, the number of measurements taken, N, and an indicator of the measurement quality, Q. Message 401 includes some or all of these fields for each w ~ base station ~ signal it detects. The 'neighbour cell' information, i.e. ~
the data .9_ corresponding to the base stations other than the serving base station with which the handset is in two-way communication, are sorted in order of received signal power level. Communications links 340 & 341 are available as part of the telephony system.
Additional data are supplied to the LCU 320 from a Management Data Base 330 which include the identifications (CI, BSIC etc.) and locations of the fixed base stations, their corresponding antenna patterns and their transmitter powers.
These data are available as part of the standard configuration of the network contained within the management system database 330 and communicated to, the LCU by link 350. The LCU 320 uses the message 401 to create a list 402 of visible base stations.
The table shows the base station identifier, the quality value associated with the detection, the transmitting antenna pattern type, the pointing direction of the main lobe of the transmitter antenna, a flag indicating whether the detected signal was from the serving ,base station (S) or was a neighbour base station (N), and the position of the base station.
The LCU 320 also generates a list 403 of the base stations known to be active in the system using the management information, but which are not present in the list 402 of visible base stations.
The handset 310 is in time-synchronous communication with its serving base station 304 and the TA value is also communicated to the LCU 320 from the network via the link 341 in another message (not shown in Figures 3 and 4).
In order to determine the position of the handset 310, the LCU 320 uses the Fists of visible 402 and not-visible 403 base stations, and the TA value associated with the serving base station, to construct an function which describes the probability that the handset is at a given horizontal position within the GSM network, given the TA
. value and the membership of the visible and the not-visible lists. The position which corresponds to the maximum probability is returned as the position of the handset.
The calculations involved in the specific embodiment are now described.
In order to determine the position of the handset, the LCU 320 finds the position which minimises the function CTOT - CTA + CV + CN , ( 1 ) where CTOT is the function whose minimum indicates the most.likely position of the handset, CTA is the function associated with the timing advance measurement, Cv is the function associated with the base stations in the visible list, and CN is the function associated with the base stations which are in the not-visible list.
The functions C~ and CN are both similarly evaluated as a sum over all N base stations as ~%~ _ -~ ~09e(Pn )~ (2) n =1 where P" is the probability that the n~' base station is visible (in the case of C"). There are two preferred ways, of calculating this value, the Gaussian Version and the Power Version, as follows.
Gaussian Version In this version the probability of the handset detecting a signal from a transmitter is modelled by a Gaussian probability distribution. A two-dimensional Cartesian coordinate system is set up centred on the base station with the x-axis ~ oriented along the direction of the transmitting antenna's forward lobe. The Gaussian distribution is centred on the x-axis at position (xo,0): The probability, P, that the handset at vector position x detects.the signals from the transmitter at position (0,0) is given by .
P(x)=expj -2(x-xo)T ~A~(x-xo)J , ( ) where x = I x I x . = I x~ I A = I 1/ dX 0 I x and y are the co-ordinates of the LyJ ' ~ . LO J ' , . L 0 . 1/6YJ' receiver, and ( )T denotes the transpose of the matrix within the braces. That is, the distribution has its peak at the position xo on the x-axis, and has a standard deviation of 6 X in the x-direction and 6 y in the y-direction with no covariance terms.
It is convenient to work with a cost function, 1(x,y), which is proportional to the negative of the logarithm of P(x,y), i.e.
f(x~Y)=-210geP(X~Y)~ . ~ ('4') The definition of f ~x,y~ given above uses coordinates relative to the position and bearing of a particular base station. When combining the distributions for all the cells in the visible list, it is convenient to use a reference coordinate system which has an arbitrary origin somewhere close to the expected position of the handset and in which the x-axis is aligned to the east and the y-axis is aligned to the north, so that a coordinate transformation must be applied to each term in the summation. If a base station k has location ~xk,yk ~ relative to the reference coordinates and bearing ~k measured in radians anti-clockwise from east, then its contribution to the cost function will equal f(xk~yk)=(X-zo)r~B~(x-zo) , {5) where zo - . xk +Rr ~I xo ~, B=RT ~~ 1/aX 0 (~R, .and R=~ cos~k sink ~
yk L 0 ' 0 1 ~ 6y J ~ L-Sln ~k COS ~k J ~
The parameters xo, aX and 6y vary from cell to cell and can be derived using heuristic methods. In one implementation of the present invention, the results of measurements made at many known locations across the network were used to find the values which, gave the lowest overall error when the positions ,were re-calculated using the present method. Once calibrated in this fashion, the method was then used to locate handsets at unknown positions.
The formula for addition of functions of the form f ~xk, yk ~ is such 'that ~X-Z~~T ~B~ ~~X-Z~~+~X-Z2~T ~B2 ~~X-12~=~X-ZR~T ~BR ~~X-.zR)+e, (6) where BR = B~ + BZ and zR = BR ~ ~B~ ~ z~ + B2 ~ z2 ~ . The term a is a constant, and hence is irrelevant , since we are interested only iri the location of the function's minimum and not, its absolute value. At the minimum point, the gradient of the function is zero. Since the gradient of a constant is always zero, the term a makes no contribution.
From this it is clear that the total cost function, C, for all the base stations can be reduced, to a simple expression of the form C=(x-z$~T .'B$'(x-z$~
where the values of z$ and B$ are derived from the positions, and bearings of all of the base stations in the visible list by repeated use of equation 6. The function C has a unique minimum at the point where x = z$ , so there is no need in this case to perform a minimisation through the use of an iterative numerical technique.
Power Version In order to calculate P using this method it is first necessary to find the power received on average by the handset from a given base station. The average .is used here so thafi no account need be taken ~of the handset receiver's antenna pattern, which in real use points in any arbitrary direction. If the handset is at distance r from the base station and,at bearing angle 6 relative to the transmitter main lobe axis, then the average signal power. received by the handset is given by . , R(r,9)=G(9)INoyr ~, . (8) where R(r,9) is the average received signal power, G(B) is the transmitter antenna gain, lNo is the transmitted power, y is a log-normal random variable (see below), and ~3 is the exponent of power loss. The assumption is made that the handset has an omni-directional antenna pattern,, which is not true in each instance in practice, but is a good model for the average: The log-normal variable represents the variability of the signal power characteristics of the radio path (see below).
Taking the logarithm of both sides of,equation (8) produces loge (R(r,9)) = loge(G (9))+ loge VI/o + loge y- /3loge r . (9) The log normal variable can be modelled as Y = YoL'n' ~ ~ ~ (10) .
-.13-where ya is the median value and m is a normally-distributed random variable with a variance of 62. Substituting equation 10 into equation 9 gives, loge (R(r,9)) = loge(G (8))+ loge IiVo + loge yo + m - ~3loge r . , (11 ) As explained in the following paragraph, this indicates that the .logarithm of the received , power level is a Gaussian distributed variable, with a probability density function (PDF) given by PDF(IOge R ) _ (2~c~2 )'~ 2 exp[-(loge R -L ~ /(2~2 )] ~, (12) where L = IOge(G (~~)+ loge tNo + loge yo - /.i' loge r .
The handset detects a base station and lists it as 'visible' if the received power level is above the minimum detectable signal power, Rmj~, This is depicted diagrammatically in Figure 5. The graph' shows the logarithm of received power, log R, plotted against the logarithm of the .distance from a base-station transmitter, log r, on axes 500 and 501. respectively and the result is a straight line 502. An example point is shown by label ~A at a distance from the transmitter indicated by line 503 and with a predicted received power given by line 504. When many measurements of the received power are made at the same distance from the .transmitter, the results are found in practice to vary and in general follow the probability distribution shown by the curve 505. This is a normal distribution on a logarithmic plot; and hence this variability is best described by the log-normal variable, y, as discussed previously. , Line 506 is the minimum receive power, Rm,~, able to be detected. by the receiver. The probability, P, of detecting the signals at point A is therefore given by the area under the curve 505 above the threshold line 506. Mathematically this is given by the integral P = probablility that loge R > loge Rn,;~ , . ~ (13) i.e.
(~ya.2~ ~~2eXp~~(~O9eR _L~/(2~z)~ d(~o9eRJ. (14) IOge Rmin ~An integral of this sort is widely known to be related to the,'Error Function', Ert~(see for example, , in Numerical Recipes in C++ by William H. Press et al, Cambridge University Press, Chapter 6 page 225). Hence P is given by p _.(~..~- Erf ((L -loge Rmln )/ ~a-))/ 2 . ~ (15) As with the Gaussian version described previously, it is mathematically convenient to use a 'cost' function, C, which is the negative logarithm of P, i.e. .
(16) C =-lOge~P~~
This function has a minimum value at the position where the probability is a .
maximum. This cost function is plotted in Figure 1 under the assumption that the antenna pattern, G(~, is a cos29function. The x and y horizontal axes are the,position axes relative to the base station position. The vertical axis is the cost function value shown on an arbitrary scale of 0 (low cost) to 10 (high cost).
There are several distinctive features of the plotted surface in Figure 1. The first is that there is a low central region 101 around the transmitter base station position where the, cost of being detected by the receiver is a minimum. This is because the power' levels there are so high that even with an unusual degree of channel attenuation variation the received power level is still above the minimum.
The second is that that there is a,high cost 103 in the directions of the antenna nulls. .
When this ~ analysis is extended to include all the signals from all the BTSs, and combined into a single cost function, the shape is similar to that shown in Figure 2. The axes in this figure are centred on the handset's position 202.
In the preferred embodiment the corresponding cost function associated with the timing advance measurement, CTA, is as shown ~ diagrammatically in Figure 6.
This figure depicts a view of a base station 601 and a graphical representation of the timing advance distance, 602, which would, ideally, be calculated from a measurement of the TA value by multiplying by the speed of radio waves. In practice in a GSM system, the reported TA value measured at a particular radial distance, 602, is quantised by the measurement method into units of about 550 metres or half a GSM 'bit'. This means that the radial distance from the base station, 602 (DTA), has on average a plus/minus quarter-bit equivalent distance error, as ~ shown by the concentric circles ,603. The cost function 604-606 which is applied in the specific embodiment is shown superimposed above this diagram. In summary any positions which are closer to the base station than a quarter of a bit from the measured TA
distance (DTA - 225m) are given a cost value 604 proportional to the difference from this position. Position within the ambiguous band of plus/minus a quarter of a bif around the measured distance are given a zero cost value 605. Positions further.
away than a quarter bit beyond the TA distance (DTA + 225m) are given a value proportional to the square of the difference from this position 606.
~In the preferred embodiment the minimum of the combined cost function CTOT
is found as follows. The base station positions are projected onto a horizontal plane, and points, are identified forming an equal distance grid in the two horizontal directions within this horizontal plane common to all base stations. This is shown diagrammatically in Figure 7 with the grid 700 shown superimposed on the horizontal plane common to all base stations. An example point 701 is identified on the plane.
The total . cost function CTOT is evaluated at each . of these grid points and that corresponding to the lowest value of CTOr is found. In Figure 7 this is .shown as example point 701. This point is then taken as the centre of another square array of grid points at one tenth ~of the spacing 702. CTOr is again evaluated at each point of this sub-grid and the minimum found. This point is then used as the centre of yet another set of grid points at which CTOr is evaluated and so on, until the sub-grid spacing is equal to 10 m. The point corresponding to the lowest value of CTOr is used as the estimate of the position of the mobile receiver.
Trials according the specific embodiments described under'Gaussian Version' and 'Power Version' above have been carried out in several GSM networks. For example, in ohe of these the Gaussian Version was used. A total of 319 test points were selected which were scattered randomly over the test area, and at each the position of a handset measured according to the invention was compared with the true position as measured by more accurate means ,(the 'ground truth'). The results showed that two thirds of the positions were within 220 m of ground truth, and that nineteen out of twenty were within 380 m.
of detection of the signals by the mobile receiver. For mathematical convenience, the cost function combining all the characteristics affecting the visibility of the transmitters is ri~ade proportional to the negative of the logarithm of the' probability of detecting.
their signals (i.e. the logarithm of the reciprocal of the probability). Hence taking .one 10' transmitter in isolation, the cost function describing its visibility has a low value when the receiver is close to the transmitter position.and directly in line with the transmitter antenna's forward lobe as shown, for example, in Figure ~1. The cost function gradually increases as the receiver moves directly away from the base station, or moves tangentially at an angle from the direction of the forward lobe at fixed distance from the transmitter. The corresponding cost functions describing each transmitter are then combined as a sum to give an overall cost furiction that describes the combined visibility of the transmitters as shown for example in .Figure 2. The minimum of this overall cost function is the position of the receiver according to the invention.
According to the .present invention there is provided a mobile receiver device .
capable of detecting signals from a plurality of transmitters and arranged to compute a list of signals received by the mobile receiver and to send said list, in use, to a calculation device.
Furthermore, according to the present invention there is provided a calculation device for use in location determination of a mobile receiver, the calculation device being arranged to:
a) receive a list of the signals detected by a receiver, b) construct a position function which relates to the probability that .the receiver is at a given ~ location, which position , has components corresponding to each member of the list, each component including a predetermined probability function that depends on the mobile ,receiver position; and ' .
c) evaluate said position function to determine a position that corresponds to the highest probability of location, and defining that. position to be the location of the mobile receiver. w 35According to the present invention there is provided a system comprising: a mobile receiver device and a calculation device.
_g_ The calculation device.may be attached to the mobile receiver, or it maybe remote from, but in communication with, the mobile receiver.
Brief description of the drawings ~ Figure 1 shows an example plot of the position function ('Cost') associated with a single transmitter.
Figure 2 shows an example plot of the position function ('Cost') associated with multiple transmitters.
Figure 3 illustrates an example multiple transmitter system.
Figure 4 shows an example list of detected signals sent by the handset, and examples of lists maintained at the location calculation unit .of visible and not-visible .
signals. .
Figure 5 illustrates the variation of received power with distance from the transmitter.
Figure 6 shows an example fixed transmitter and a plot of the position function component associated .with timing ,advance.
Figure 7, illustrates a calculation grid. ..
Detailed descriptions of particular embodiments A mobile receiver positioning system has been implemented according to the current invention using an existing GSM mobile telephony network (see Figure 3).
The system 300 comprises a plurality of transceiver base stations 301, 302, 303, 304, which act as the transmitters described above, a handset which acts as the mobile receiver 310, a location calculation unit (LCU) 320, and a management database (MDB) 330. The base stations provide the signal transmissions which are detected' by the handset 310 and processed in the LCU 320 using the mathematical algorithms described below.
The handset 310 detects signals from a riumber of the nearby base stations 301-304, compiles a list of those it detects, and sends the list via a message over link 340 to the serving base station 304: The format and content of the message is shown in Figure 4 at 401. It is defined by the GSM standard and includes fields for the Broadcast Control Channel; BCCH, the Base Station Identity Code, BSIC, the Mobile Network Code, MNC, the Mobile Country Code, MCC, the Cell Identifier, CI, the Local . Area Code, LAC, the number of measurements taken, N, and an indicator of the measurement quality, Q. Message 401 includes some or all of these fields for each w ~ base station ~ signal it detects. The 'neighbour cell' information, i.e. ~
the data .9_ corresponding to the base stations other than the serving base station with which the handset is in two-way communication, are sorted in order of received signal power level. Communications links 340 & 341 are available as part of the telephony system.
Additional data are supplied to the LCU 320 from a Management Data Base 330 which include the identifications (CI, BSIC etc.) and locations of the fixed base stations, their corresponding antenna patterns and their transmitter powers.
These data are available as part of the standard configuration of the network contained within the management system database 330 and communicated to, the LCU by link 350. The LCU 320 uses the message 401 to create a list 402 of visible base stations.
The table shows the base station identifier, the quality value associated with the detection, the transmitting antenna pattern type, the pointing direction of the main lobe of the transmitter antenna, a flag indicating whether the detected signal was from the serving ,base station (S) or was a neighbour base station (N), and the position of the base station.
The LCU 320 also generates a list 403 of the base stations known to be active in the system using the management information, but which are not present in the list 402 of visible base stations.
The handset 310 is in time-synchronous communication with its serving base station 304 and the TA value is also communicated to the LCU 320 from the network via the link 341 in another message (not shown in Figures 3 and 4).
In order to determine the position of the handset 310, the LCU 320 uses the Fists of visible 402 and not-visible 403 base stations, and the TA value associated with the serving base station, to construct an function which describes the probability that the handset is at a given horizontal position within the GSM network, given the TA
. value and the membership of the visible and the not-visible lists. The position which corresponds to the maximum probability is returned as the position of the handset.
The calculations involved in the specific embodiment are now described.
In order to determine the position of the handset, the LCU 320 finds the position which minimises the function CTOT - CTA + CV + CN , ( 1 ) where CTOT is the function whose minimum indicates the most.likely position of the handset, CTA is the function associated with the timing advance measurement, Cv is the function associated with the base stations in the visible list, and CN is the function associated with the base stations which are in the not-visible list.
The functions C~ and CN are both similarly evaluated as a sum over all N base stations as ~%~ _ -~ ~09e(Pn )~ (2) n =1 where P" is the probability that the n~' base station is visible (in the case of C"). There are two preferred ways, of calculating this value, the Gaussian Version and the Power Version, as follows.
Gaussian Version In this version the probability of the handset detecting a signal from a transmitter is modelled by a Gaussian probability distribution. A two-dimensional Cartesian coordinate system is set up centred on the base station with the x-axis ~ oriented along the direction of the transmitting antenna's forward lobe. The Gaussian distribution is centred on the x-axis at position (xo,0): The probability, P, that the handset at vector position x detects.the signals from the transmitter at position (0,0) is given by .
P(x)=expj -2(x-xo)T ~A~(x-xo)J , ( ) where x = I x I x . = I x~ I A = I 1/ dX 0 I x and y are the co-ordinates of the LyJ ' ~ . LO J ' , . L 0 . 1/6YJ' receiver, and ( )T denotes the transpose of the matrix within the braces. That is, the distribution has its peak at the position xo on the x-axis, and has a standard deviation of 6 X in the x-direction and 6 y in the y-direction with no covariance terms.
It is convenient to work with a cost function, 1(x,y), which is proportional to the negative of the logarithm of P(x,y), i.e.
f(x~Y)=-210geP(X~Y)~ . ~ ('4') The definition of f ~x,y~ given above uses coordinates relative to the position and bearing of a particular base station. When combining the distributions for all the cells in the visible list, it is convenient to use a reference coordinate system which has an arbitrary origin somewhere close to the expected position of the handset and in which the x-axis is aligned to the east and the y-axis is aligned to the north, so that a coordinate transformation must be applied to each term in the summation. If a base station k has location ~xk,yk ~ relative to the reference coordinates and bearing ~k measured in radians anti-clockwise from east, then its contribution to the cost function will equal f(xk~yk)=(X-zo)r~B~(x-zo) , {5) where zo - . xk +Rr ~I xo ~, B=RT ~~ 1/aX 0 (~R, .and R=~ cos~k sink ~
yk L 0 ' 0 1 ~ 6y J ~ L-Sln ~k COS ~k J ~
The parameters xo, aX and 6y vary from cell to cell and can be derived using heuristic methods. In one implementation of the present invention, the results of measurements made at many known locations across the network were used to find the values which, gave the lowest overall error when the positions ,were re-calculated using the present method. Once calibrated in this fashion, the method was then used to locate handsets at unknown positions.
The formula for addition of functions of the form f ~xk, yk ~ is such 'that ~X-Z~~T ~B~ ~~X-Z~~+~X-Z2~T ~B2 ~~X-12~=~X-ZR~T ~BR ~~X-.zR)+e, (6) where BR = B~ + BZ and zR = BR ~ ~B~ ~ z~ + B2 ~ z2 ~ . The term a is a constant, and hence is irrelevant , since we are interested only iri the location of the function's minimum and not, its absolute value. At the minimum point, the gradient of the function is zero. Since the gradient of a constant is always zero, the term a makes no contribution.
From this it is clear that the total cost function, C, for all the base stations can be reduced, to a simple expression of the form C=(x-z$~T .'B$'(x-z$~
where the values of z$ and B$ are derived from the positions, and bearings of all of the base stations in the visible list by repeated use of equation 6. The function C has a unique minimum at the point where x = z$ , so there is no need in this case to perform a minimisation through the use of an iterative numerical technique.
Power Version In order to calculate P using this method it is first necessary to find the power received on average by the handset from a given base station. The average .is used here so thafi no account need be taken ~of the handset receiver's antenna pattern, which in real use points in any arbitrary direction. If the handset is at distance r from the base station and,at bearing angle 6 relative to the transmitter main lobe axis, then the average signal power. received by the handset is given by . , R(r,9)=G(9)INoyr ~, . (8) where R(r,9) is the average received signal power, G(B) is the transmitter antenna gain, lNo is the transmitted power, y is a log-normal random variable (see below), and ~3 is the exponent of power loss. The assumption is made that the handset has an omni-directional antenna pattern,, which is not true in each instance in practice, but is a good model for the average: The log-normal variable represents the variability of the signal power characteristics of the radio path (see below).
Taking the logarithm of both sides of,equation (8) produces loge (R(r,9)) = loge(G (9))+ loge VI/o + loge y- /3loge r . (9) The log normal variable can be modelled as Y = YoL'n' ~ ~ ~ (10) .
-.13-where ya is the median value and m is a normally-distributed random variable with a variance of 62. Substituting equation 10 into equation 9 gives, loge (R(r,9)) = loge(G (8))+ loge IiVo + loge yo + m - ~3loge r . , (11 ) As explained in the following paragraph, this indicates that the .logarithm of the received , power level is a Gaussian distributed variable, with a probability density function (PDF) given by PDF(IOge R ) _ (2~c~2 )'~ 2 exp[-(loge R -L ~ /(2~2 )] ~, (12) where L = IOge(G (~~)+ loge tNo + loge yo - /.i' loge r .
The handset detects a base station and lists it as 'visible' if the received power level is above the minimum detectable signal power, Rmj~, This is depicted diagrammatically in Figure 5. The graph' shows the logarithm of received power, log R, plotted against the logarithm of the .distance from a base-station transmitter, log r, on axes 500 and 501. respectively and the result is a straight line 502. An example point is shown by label ~A at a distance from the transmitter indicated by line 503 and with a predicted received power given by line 504. When many measurements of the received power are made at the same distance from the .transmitter, the results are found in practice to vary and in general follow the probability distribution shown by the curve 505. This is a normal distribution on a logarithmic plot; and hence this variability is best described by the log-normal variable, y, as discussed previously. , Line 506 is the minimum receive power, Rm,~, able to be detected. by the receiver. The probability, P, of detecting the signals at point A is therefore given by the area under the curve 505 above the threshold line 506. Mathematically this is given by the integral P = probablility that loge R > loge Rn,;~ , . ~ (13) i.e.
(~ya.2~ ~~2eXp~~(~O9eR _L~/(2~z)~ d(~o9eRJ. (14) IOge Rmin ~An integral of this sort is widely known to be related to the,'Error Function', Ert~(see for example, , in Numerical Recipes in C++ by William H. Press et al, Cambridge University Press, Chapter 6 page 225). Hence P is given by p _.(~..~- Erf ((L -loge Rmln )/ ~a-))/ 2 . ~ (15) As with the Gaussian version described previously, it is mathematically convenient to use a 'cost' function, C, which is the negative logarithm of P, i.e. .
(16) C =-lOge~P~~
This function has a minimum value at the position where the probability is a .
maximum. This cost function is plotted in Figure 1 under the assumption that the antenna pattern, G(~, is a cos29function. The x and y horizontal axes are the,position axes relative to the base station position. The vertical axis is the cost function value shown on an arbitrary scale of 0 (low cost) to 10 (high cost).
There are several distinctive features of the plotted surface in Figure 1. The first is that there is a low central region 101 around the transmitter base station position where the, cost of being detected by the receiver is a minimum. This is because the power' levels there are so high that even with an unusual degree of channel attenuation variation the received power level is still above the minimum.
The second is that that there is a,high cost 103 in the directions of the antenna nulls. .
When this ~ analysis is extended to include all the signals from all the BTSs, and combined into a single cost function, the shape is similar to that shown in Figure 2. The axes in this figure are centred on the handset's position 202.
In the preferred embodiment the corresponding cost function associated with the timing advance measurement, CTA, is as shown ~ diagrammatically in Figure 6.
This figure depicts a view of a base station 601 and a graphical representation of the timing advance distance, 602, which would, ideally, be calculated from a measurement of the TA value by multiplying by the speed of radio waves. In practice in a GSM system, the reported TA value measured at a particular radial distance, 602, is quantised by the measurement method into units of about 550 metres or half a GSM 'bit'. This means that the radial distance from the base station, 602 (DTA), has on average a plus/minus quarter-bit equivalent distance error, as ~ shown by the concentric circles ,603. The cost function 604-606 which is applied in the specific embodiment is shown superimposed above this diagram. In summary any positions which are closer to the base station than a quarter of a bit from the measured TA
distance (DTA - 225m) are given a cost value 604 proportional to the difference from this position. Position within the ambiguous band of plus/minus a quarter of a bif around the measured distance are given a zero cost value 605. Positions further.
away than a quarter bit beyond the TA distance (DTA + 225m) are given a value proportional to the square of the difference from this position 606.
~In the preferred embodiment the minimum of the combined cost function CTOT
is found as follows. The base station positions are projected onto a horizontal plane, and points, are identified forming an equal distance grid in the two horizontal directions within this horizontal plane common to all base stations. This is shown diagrammatically in Figure 7 with the grid 700 shown superimposed on the horizontal plane common to all base stations. An example point 701 is identified on the plane.
The total . cost function CTOT is evaluated at each . of these grid points and that corresponding to the lowest value of CTOr is found. In Figure 7 this is .shown as example point 701. This point is then taken as the centre of another square array of grid points at one tenth ~of the spacing 702. CTOr is again evaluated at each point of this sub-grid and the minimum found. This point is then used as the centre of yet another set of grid points at which CTOr is evaluated and so on, until the sub-grid spacing is equal to 10 m. The point corresponding to the lowest value of CTOr is used as the estimate of the position of the mobile receiver.
Trials according the specific embodiments described under'Gaussian Version' and 'Power Version' above have been carried out in several GSM networks. For example, in ohe of these the Gaussian Version was used. A total of 319 test points were selected which were scattered randomly over the test area, and at each the position of a handset measured according to the invention was compared with the true position as measured by more accurate means ,(the 'ground truth'). The results showed that two thirds of the positions were within 220 m of ground truth, and that nineteen out of twenty were within 380 m.
Claims (28)
1. A method of determining the location of a mobile receiver in one or more dimensions in a radio system which includes a network of transmitters the positions of which are known, the method comprising the steps of:
a) compiling a list of the signals detected by the receiver;
b) sending said list to a calculation device;
c) constructing a position function which relates to the probability that the receiver is at a given location, which position function has components corresponding to each member of the list, each component including a predetermined probability function that depends on mobile receiver position;
and d) evaluating said position function to determine a position that corresponds to the highest probability of location, and defining that position to be the location of the mobile receiver.
a) compiling a list of the signals detected by the receiver;
b) sending said list to a calculation device;
c) constructing a position function which relates to the probability that the receiver is at a given location, which position function has components corresponding to each member of the list, each component including a predetermined probability function that depends on mobile receiver position;
and d) evaluating said position function to determine a position that corresponds to the highest probability of location, and defining that position to be the location of the mobile receiver.
2. A method according to claim 1, wherein the detection of signals by the mobile receiver includes the reception of signals transmitted by a particular transmitter of the network on a particular frequency above a given threshold power.
3. A method according to claim 2, which includes the decoding of information from the signals.
4. A method according to claim 3, wherein the list of signals includes decoded information related to the identity of each transmitter.
5. A method according to any of the preceding claims, which includes the detection of a particular signature in,the signal.
6. A method according to any of the preceding claims, wherein the list of signals includes a list of frequencies on which the signals were received.
7. A method according to any of the preceding claims, wherein the list contains a list of references, each reference being representative of respective transmitter information, maintained by the calculation device.
8. A method according to any of the preceding claims, wherein the predetermined probability function is based on a propagation model for the radio waves between each transmitter and the receiver.
9. A method according to claim 8, wherein the propagation model is based on any or all of the transmitter power, transmitter antenna pattern, travel distance, bearing of the receiver from the transmitter, variable transmission characteristics of radio signal propagation, timing advance, or local topographical data.
10. A method according to any of the preceding claims, wherein the predetermined probability function is a Gaussian function with its origin at a pre-defined point from the transmitter.
11. A method according to any of the preceding claims, wherein the predetermined probability function is based on a heuristic derived from the probability of detections obtained at calibration points.
12. A method according to any of the preceding claims, wherein the predetermined probability function is proportional to the probability of the detection of the signals from a transmitter.
13. A method according to any of the preceding claims, in which the position function is a cost function which has its lowest numerical value at the position of maximum probability.
14. A method according to claim 13, wherein the cost function is proportional to minus the logarithm of the predetermined probability function.
15. A method according to any of the preceding claims, wherein a further list is created in the calculation device of information of transmitter signals not detected by the mobile receiver.
16. A method according to any of the preceding claims, wherein the position function is modified to take account of the list of signals not detected by the mobile receiver.
17. A method according to any of the preceding claims, wherein the timing advance is included in the position function.
18. A method according to any of the preceding claims, wherein position function is modified to take account of the reliability of the detection of each of the signals by the receiver.
19. A method according to claim 18, in which the reliability measurement is dependent upon the difference between the received power and the detection threshold.
20. A method according to claim 18 or 19, in which the reliability measurement includes the number of times the signal has been detected.
21. A method according to any of the preceding claims, in which the list sent from the receiver to the calculation device is sorted.
22. A method according to claim 21 in which the list sent from the receiver to the calculation device is sorted in terms of power received by the mobile receiver.
23. A method according claim 21, in which the list sent from the mobile receiver to the calculation device is sorted in terms of the reliability of signal reception.
24. A mobile receiver device capable of detecting signals from a plurality of transmitters and arranged to compute a list of signals received by the mobile receiver and to send said list, in use, to a calculation device.
25. A calculation device for use in location determination of a mobile receiver, the calculation device being arranged to:
a) receive a list of the signals detected by a receiver, b) construct a position function which relates to the probability that the receiver is at a given location, which position has components corresponding to each member of the list, each component including a predetermined probability function that depends on the mobile receiver position; and c) evaluate said position function to determine a position that corresponds to the highest probability of location, and defining that position to be the location, of the mobile receiver.
a) receive a list of the signals detected by a receiver, b) construct a position function which relates to the probability that the receiver is at a given location, which position has components corresponding to each member of the list, each component including a predetermined probability function that depends on the mobile receiver position; and c) evaluate said position function to determine a position that corresponds to the highest probability of location, and defining that position to be the location, of the mobile receiver.
26. A system comprising a mobile receiver device according to claim 24 and a calculation device according to claim 25.
27. A system according to claim 26, wherein the calculation device is attached to the mobile receiver.
28. A system according to claim 26, wherein the calculation device is remote from, but in communication with the mobile receiver.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP03250733A EP1445970B1 (en) | 2003-02-05 | 2003-02-05 | A method and system for locating a mobile radio receiver in a radio system with multiple tranmitters |
| EP03250733.7 | 2003-02-05 | ||
| PCT/GB2004/000281 WO2004071120A1 (en) | 2003-02-05 | 2004-01-23 | A method and system for locating a mobile radio receiver in a radio system with multiple transmitters |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CA2513154A1 true CA2513154A1 (en) | 2004-08-19 |
Family
ID=32605415
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002513154A Abandoned CA2513154A1 (en) | 2003-02-05 | 2004-01-23 | A method and system for locating a mobile radio receiver in a radio system with multiple transmitters |
Country Status (8)
| Country | Link |
|---|---|
| US (1) | US8630656B2 (en) |
| EP (1) | EP1445970B1 (en) |
| KR (1) | KR101108749B1 (en) |
| CA (1) | CA2513154A1 (en) |
| DE (1) | DE60326925D1 (en) |
| HK (1) | HK1064251A1 (en) |
| TW (1) | TWI345641B (en) |
| WO (1) | WO2004071120A1 (en) |
Families Citing this family (50)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7742776B2 (en) * | 2002-05-31 | 2010-06-22 | Telecom Italia S.P.A. | Method for locating mobile terminals, system and components therefor |
| US8971913B2 (en) * | 2003-06-27 | 2015-03-03 | Qualcomm Incorporated | Method and apparatus for wireless network hybrid positioning |
| US8483717B2 (en) | 2003-06-27 | 2013-07-09 | Qualcomm Incorporated | Local area network assisted positioning |
| US7904244B2 (en) * | 2003-11-18 | 2011-03-08 | Sarimo Technologies, LLC | Determining a location or position using information from multiple location and positioning technologies and applications using such a determined location or position |
| US7319878B2 (en) | 2004-06-18 | 2008-01-15 | Qualcomm Incorporated | Method and apparatus for determining location of a base station using a plurality of mobile stations in a wireless mobile network |
| US7881905B2 (en) | 2004-11-17 | 2011-02-01 | Qualcomm Incorporated | Method for ambiguity resolution in location determination |
| FR2879068A1 (en) * | 2004-12-06 | 2006-06-09 | France Telecom | DETERMINING A DISPLACEMENT STATE OF A MOBILE TERMINAL |
| FR2879067A1 (en) * | 2004-12-06 | 2006-06-09 | France Telecom | LOCATION OF A MOBILE TERMINAL IN RELATION TO A BUILDING |
| CN101180550A (en) | 2005-03-18 | 2008-05-14 | 探索无线公司 | Enhanced mobile location method and system |
| CN101171529A (en) * | 2005-03-18 | 2008-04-30 | 探索无线公司 | Enhanced mobile location |
| US8463285B2 (en) * | 2005-04-08 | 2013-06-11 | Wavemarket, Inc. | Systems and methods for mobile terminal location determination using profiles of radio signal parameter measurements |
| US7257413B2 (en) * | 2005-08-24 | 2007-08-14 | Qualcomm Incorporated | Dynamic location almanac for wireless base stations |
| EP1941758A4 (en) * | 2005-10-24 | 2012-04-25 | Wavemarket Inc D B A Location Labs | Detection in mobile service maintenance |
| CN101356765A (en) * | 2005-11-04 | 2009-01-28 | 探索无线公司 | Profile based communications service |
| US9042917B2 (en) * | 2005-11-07 | 2015-05-26 | Qualcomm Incorporated | Positioning for WLANS and other wireless networks |
| KR101071076B1 (en) * | 2005-11-07 | 2011-10-10 | 퀄컴 인코포레이티드 | Positioning for wlans and other wireless networks |
| US7626969B2 (en) * | 2006-10-04 | 2009-12-01 | Cisco Technology, Inc. | Relative location of a wireless node in a wireless network |
| FR2907620B1 (en) * | 2006-10-23 | 2009-02-06 | Inrets | METHOD AND DEVICE FOR LOCALIZATION BY DETECTING TRANSMITTERS OF RADIO SIGNALS. |
| WO2008115209A2 (en) * | 2006-10-25 | 2008-09-25 | Massachusetts Institute Of Technology | Cooperative localization for wireless networks |
| US9226257B2 (en) | 2006-11-04 | 2015-12-29 | Qualcomm Incorporated | Positioning for WLANs and other wireless networks |
| JP4814060B2 (en) * | 2006-11-09 | 2011-11-09 | 株式会社山武 | Position estimation method and position estimation system |
| US8798639B2 (en) * | 2007-01-17 | 2014-08-05 | Qualcomm Incorporated | Method and apparatus for using historic network information for determining approximate position |
| AU2008301216A1 (en) * | 2007-09-17 | 2009-03-26 | Seeker Wireless Pty Limited | Systems and methods for triggering location based voice and/or data communications to or from mobile radio terminals |
| US20090109885A1 (en) * | 2007-10-29 | 2009-04-30 | Motorola, Inc. | Method and Apparatus for Reducing Energy Consumption in Nodes by Adjusting Carrier Sensing Thresholds |
| US8737985B2 (en) | 2007-11-26 | 2014-05-27 | Wavemarket, Inc. | Methods and systems for zone creation and adaption |
| DE602007012912D1 (en) * | 2007-12-13 | 2011-04-14 | Swisscom Ag | System and method for finding a mobile communication terminal in a geographical area |
| EP2071355B1 (en) | 2007-12-13 | 2015-07-29 | Swisscom AG | System and method for determining a location area of a mobile user |
| US8803737B2 (en) | 2008-02-29 | 2014-08-12 | Apple Inc. | Location determination |
| WO2009124348A1 (en) * | 2008-04-07 | 2009-10-15 | Seeker Wireless Pty Limited | Efficient collection of wireless transmitter characteristics |
| US8213389B2 (en) | 2008-04-15 | 2012-07-03 | Apple Inc. | Location determination using formula |
| WO2011091062A1 (en) * | 2010-01-20 | 2011-07-28 | Worcester Polytechnic Institute | Precision location method and system |
| US9244175B2 (en) | 2010-03-17 | 2016-01-26 | Nokia Technologies Oy | Method and apparatus for testing received signals in a radio signal positioning system |
| US8618984B2 (en) * | 2010-03-19 | 2013-12-31 | Microsoft Corporation | Selecting beacons for location inference |
| US8244236B2 (en) | 2010-04-29 | 2012-08-14 | Wavemarket, Inc. | System and method for aggregating and disseminating mobile device tag data |
| US8504077B2 (en) | 2010-12-04 | 2013-08-06 | Wavemarket, Inc. | System and method for monitoring and disseminating mobile device location information |
| IT1404537B1 (en) | 2011-02-25 | 2013-11-22 | Sisvel Technology Srl | METHOD FOR ESTIMATING THE DISTANCE OF A RECEIVER FROM A RADIO TRANSMITTER, RELATED METHODS FOR CALCULATING THE POSITION OF A MOBILE TERMINAL, MOBILE TERMINAL AND DEVICE. |
| WO2013052651A2 (en) * | 2011-10-04 | 2013-04-11 | North Carolina State University | Receiver-based methods, systems, and computer readable media for controlling tcp sender behavior in cellular communications networks with large buffer sizes |
| US8797899B2 (en) * | 2011-12-22 | 2014-08-05 | Qualcomm Incorporated | System and method for probabilistic WLAN positioning |
| WO2013136121A1 (en) * | 2012-03-15 | 2013-09-19 | Nokia Corporation | Supporting location based operations |
| US9121922B2 (en) | 2012-06-26 | 2015-09-01 | Cambridge Silicon Radio Limited | Access point location identification methods and apparatus based on absolute and relative harvesting |
| US8909258B2 (en) | 2012-09-07 | 2014-12-09 | Cambridge Silicon Radio Limited | Context and map aiding for self-learning |
| US8705320B1 (en) * | 2012-09-13 | 2014-04-22 | Google Inc. | Audio localization using multilateration |
| CN104459617B (en) * | 2014-12-22 | 2017-08-04 | 上海斐讯数据通信技术有限公司 | The method of wireless location |
| CN104936286A (en) * | 2015-05-13 | 2015-09-23 | 中国矿业大学 | Tunnel three-channel difference range finding method based on RSSI |
| US10206193B2 (en) | 2016-08-31 | 2019-02-12 | At&T Intellectual Property, L.P. | Mobile device location based on detection of radiation sources |
| GB2563825A (en) * | 2017-06-19 | 2019-01-02 | Vodafone Ip Licensing Ltd | Localization of mobile devices |
| US10743141B2 (en) | 2018-06-05 | 2020-08-11 | Kenmar Corporation | Systems and methods for determining a location of an electronic device using bilateration |
| JP6640309B1 (en) * | 2018-11-09 | 2020-02-05 | モティーブリサーチ株式会社 | Position estimation device, program, position estimation method, information processing device, and information processing method |
| FR3107347B1 (en) * | 2020-02-13 | 2022-01-07 | Continental Automotive | Angle sensor for motor vehicle. |
| TWI774140B (en) * | 2020-11-26 | 2022-08-11 | 聚眾聯合科技股份有限公司 | Two-way signal positioning method and two-way signal positioning system thereof |
Family Cites Families (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US362783A (en) * | 1887-05-10 | Heney c | ||
| US5499032A (en) * | 1992-12-22 | 1996-03-12 | Terrapin Corporation | Navigation and positioning system and method using uncoordinated beacon signals |
| DE4225685A1 (en) * | 1992-08-04 | 1994-02-10 | Philips Patentverwaltung | Cellular mobile radio system |
| US5768686A (en) * | 1994-09-28 | 1998-06-16 | U S West, Inc. | Method and system for position determination using video dial tone |
| US5959580A (en) * | 1994-11-03 | 1999-09-28 | Ksi Inc. | Communications localization system |
| US5657487A (en) * | 1995-06-05 | 1997-08-12 | Airnet Communications Corporation | Mobile telephone location process making use of handoff data |
| GB9519087D0 (en) * | 1995-09-19 | 1995-11-22 | Cursor Positioning Sys Ltd | Navigation and tracking system |
| FI101445B1 (en) * | 1995-10-03 | 1998-06-15 | Nokia Mobile Phones Ltd | Location system of a mobile station |
| US5883598A (en) * | 1995-12-15 | 1999-03-16 | Signatron Technology Corporation | Position location system and method |
| GB2311697B (en) * | 1996-03-22 | 1999-07-28 | Matsushita Electric Ind Co Ltd | Wireless communication system and method and system for detection of position of radio mobile station |
| WO1998010307A1 (en) * | 1996-09-09 | 1998-03-12 | Dennis Jay Dupray | Location of a mobile station |
| FI106908B (en) | 1998-05-28 | 2001-04-30 | Nokia Networks Oy | Location procedure and radio system |
| DE19836778C2 (en) | 1998-08-13 | 2002-01-31 | Siemens Ag | Method for determining the position of a TDMA mobile radio mobile station |
| FR2782594B1 (en) * | 1998-08-20 | 2000-09-22 | France Telecom | METHOD FOR SIMULTANEOUS MEASUREMENT OF THE PROPAGATION CHARACTERISTICS OF A PLURALITY OF RADIOELECTRIC CHANNELS |
| US6496701B1 (en) * | 1998-08-25 | 2002-12-17 | Lucent Technologies Inc. | Pattern-recognition-based geolocation |
| JP3607516B2 (en) * | 1999-01-20 | 2005-01-05 | 松下電器産業株式会社 | Mobile map matching device |
| US6317599B1 (en) * | 1999-05-26 | 2001-11-13 | Wireless Valley Communications, Inc. | Method and system for automated optimization of antenna positioning in 3-D |
| US6493679B1 (en) * | 1999-05-26 | 2002-12-10 | Wireless Valley Communications, Inc. | Method and system for managing a real time bill of materials |
| US6263208B1 (en) * | 1999-05-28 | 2001-07-17 | Lucent Technologies Inc. | Geolocation estimation method for CDMA terminals based on pilot strength measurements |
| GB9915277D0 (en) | 1999-07-01 | 1999-09-01 | Aircom International Limited | Mobile telephone positioning system |
| US6499006B1 (en) * | 1999-07-14 | 2002-12-24 | Wireless Valley Communications, Inc. | System for the three-dimensional display of wireless communication system performance |
| US6556832B1 (en) * | 2000-02-04 | 2003-04-29 | Qualcomm Incorporated | Method and apparatus for evaluation of position location performance |
| US6625454B1 (en) * | 2000-08-04 | 2003-09-23 | Wireless Valley Communications, Inc. | Method and system for designing or deploying a communications network which considers frequency dependent effects |
| FI111901B (en) | 2000-12-29 | 2003-09-30 | Ekahau Oy | Estimation of position in wireless communication networks |
| US7224983B2 (en) * | 2003-01-08 | 2007-05-29 | Lucent Technologies Inc. | Method and apparatus for geolocation estimates in wireless networks |
-
2003
- 2003-02-05 EP EP03250733A patent/EP1445970B1/en not_active Expired - Fee Related
- 2003-02-05 DE DE60326925T patent/DE60326925D1/en not_active Expired - Lifetime
-
2004
- 2004-01-23 WO PCT/GB2004/000281 patent/WO2004071120A1/en active Application Filing
- 2004-01-23 CA CA002513154A patent/CA2513154A1/en not_active Abandoned
- 2004-01-23 US US10/544,072 patent/US8630656B2/en not_active Expired - Fee Related
- 2004-02-05 TW TW093102624A patent/TWI345641B/en not_active IP Right Cessation
- 2004-09-10 HK HK04106900.2A patent/HK1064251A1/en not_active IP Right Cessation
-
2005
- 2005-08-05 KR KR20057014524A patent/KR101108749B1/en not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| EP1445970A1 (en) | 2004-08-11 |
| HK1064251A1 (en) | 2005-01-21 |
| KR101108749B1 (en) | 2012-03-08 |
| KR20050098901A (en) | 2005-10-12 |
| WO2004071120A1 (en) | 2004-08-19 |
| EP1445970B1 (en) | 2009-04-01 |
| US8630656B2 (en) | 2014-01-14 |
| TW200428017A (en) | 2004-12-16 |
| US20070001867A1 (en) | 2007-01-04 |
| WO2004071120A8 (en) | 2005-09-01 |
| TWI345641B (en) | 2011-07-21 |
| DE60326925D1 (en) | 2009-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1445970B1 (en) | A method and system for locating a mobile radio receiver in a radio system with multiple tranmitters | |
| US20200142025A1 (en) | Observed time difference of arrival angle of arrival discriminator | |
| EP1197114B1 (en) | Selection of location measurement units for determining the position of a mobile communication station | |
| US9599696B2 (en) | Procedure to increase position location availability | |
| EP3278132B1 (en) | Method and apparatus for positioning a mobile terminal in a radio network | |
| EP1757138B1 (en) | Locating mobile terminals | |
| Yassin et al. | A survey of positioning techniques and location based services in wireless networks | |
| US9594149B2 (en) | Fingerprinting with radio channel related information | |
| US7933610B2 (en) | Method and apparatus to select an optimum site and/or sector to provide geo-location data | |
| US7570213B2 (en) | Method and apparatus for detecting and locating intrusion in a wireless network | |
| US20220271818A1 (en) | Non-Line-of-Sight Path Detection for User Equipment Positioning in Wireless Networks | |
| Kos et al. | Mobile user positioning in GSM/UMTS cellular networks | |
| KR20000062548A (en) | Method for combining multiple measurements to determine the position of a mobile transceiver | |
| KR20060096096A (en) | Estimating a position of a wireless mobile device with respect to one or more base station | |
| US20220221551A1 (en) | Method and system for geolocating an object using a mobile base station | |
| KR20120102277A (en) | Positioning method using the phase shift of base station | |
| Fischer | 5G NR positioning | |
| Fokin et al. | Topology Search Using Dilution of Precision Criterion for Enhanced 5G Positioning Service Area | |
| US20050176443A1 (en) | Cellular phone geolocation system | |
| Delgado et al. | Statistical beam information for mmW positioning | |
| Abugabal et al. | NOVEL position estimation using differential timing information for asynchronous LTE/NR networks | |
| EP4314857A1 (en) | Methods for positioning a wireless device, a related wireless node and a related location network node | |
| Buchanan et al. | Critical analysis and error determination of locating-finding techniques in GSM | |
| Singh et al. | Analysis of Cellular Positioning Techniques in UMTS Networks | |
| Singh et al. | Performance Analysis of Network Based Cellular Positioning Techniques in UMTS Networks |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FZDE | Discontinued |