|
|
| Line 1: |
Line 1: |
| '''Multilateration''' is a [[navigation]] technique based on the measurement of the ''difference'' in distance to two stations at known locations that broadcast signals at known times. Unlike measurements of absolute distance or angle, measuring the difference in distance between two stations results in an infinite number of locations that satisfy the measurement. When these possible locations are plotted, they form a [[hyperbola|hyperbolic curve]]. To locate the exact location along that curve, multilateration relies on multiple measurements: a second measurement taken to a different pair of stations will produce a second curve, which intersects with the first. When the two curves are compared, a small number of possible locations are revealed, producing a "fix".
| | Emilia Shryock is my name but you can contact me something you like. One of the things he enjoys most is ice skating but he is struggling to discover time for it. Years ago we moved to North Dakota. Hiring is her day occupation now but she's usually needed her own business.<br><br>Here is my webpage: home std test kit [[http://www.cam4teens.com/blog/84472 Going to www.cam4teens.com]] |
| | |
| Multilateration is a common technique in [[radio navigation]] systems, where it is known as [[hyperbolic navigation]]. These systems are relatively easy to construct as there is no need for a common clock, and the difference in the signal timing can be measured visibly using an [[oscilloscope]]. This formed the basis of a number of widely used navigation systems starting in [[World War II]] with the British [[Gee (navigation)|Gee]] system and several similar systems introduced over the next few decades. The introduction of the [[microprocessor]] greatly simplified operation, greatly increasing popularity during the 1980s. The most popular hyperbolic navigation system was [[LORAN-C]], which was used around the world until the system was shut down in 2010. Other systems continue to be used, but the widespread use of [[satellite navigation]] systems like [[Global Positioning System|GPS]] (also Hyperbolic) have made these systems largely redundant.
| |
| | |
| Multilateration should not be confused with [[trilateration]], which uses distances or absolute measurements of [[Time of arrival|time-of-flight]] from three or more sites, or with [[triangulation]], which uses the measurement of absolute angles. Both of these systems are also commonly used with radio navigation systems; trilateration is the basis of GPS.
| |
| | |
| == Principle ==
| |
| Multilateration is commonly used in civil and military surveillance applications to accurately locate an aircraft, vehicle or stationary emitter by measuring the "Time Difference of Arrival" (TDOA) of a signal from the emitter at three or more receiver sites.
| |
| | |
| If a pulse is emitted from a platform, it will arrive at slightly different times at two spatially separated receiver sites, the TDOA being due to the different distances of each receiver from the platform. In fact, for given locations of the two receivers, a whole set of emitter locations would give the same measurement of TDOA. Given two receiver locations and a known TDOA, the [[locus (mathematics)|locus]] of possible emitter locations is a one half of a two-sheeted [[hyperboloid]].
| |
| | |
| [[File:HyperboloidOfTwoSheets.png|thumb|Fig1. A two-sheeted hyperboloid]]
| |
| | |
| In simple terms, with two receivers at known locations, an emitter can be located onto a hyperboloid.<ref>In other words, given two receivers at known locations, one can derive a three-dimensional surface (characterized as one half of a hyperboloid) for which any two points on said surface will have the same differential distance from said receivers, i.e., a signal transmitted from any point on the surface will have the same TDOA (measured by the receivers) as a signal transmitted from any other point on the surface. <br>Therefore, in practice, the TDOA corresponding to a (moving) transmitter is measured, a corresponding hyperbolic surface is derived, and the transmitter is said to be "located" somewhere on the that surface.</ref> Note that the receivers do not need to know the absolute time at which the pulse was transmitted – only the time difference is needed.
| |
| | |
| Consider now a third receiver at a third location. This would provide a second TDOA measurement and hence locate the emitter on a second hyperboloid. The intersection of these two hyperboloids describes a curve on which the emitter lies.
| |
| | |
| If a fourth receiver is now introduced, a third TDOA measurement is available and the intersection of the resulting third hyperboloid with the curve already found with the other three receivers defines a unique point in space. The emitter's location is therefore fully determined in 3D.
| |
| | |
| In practice, errors in the measurement of the time of arrival of pulses mean that enhanced accuracy can be obtained with more than four receivers. In general, ''N'' receivers provide ''N'' − 1 hyperboloids. When there are ''N'' > 4 receivers, the ''N'' − 1 hyperboloids should, assuming a perfect model and measurements, intersect on a single point. In reality, the surfaces rarely intersect, because of various errors. In this case, the location problem can be posed as an [[Optimization (mathematics)|optimization]] problem and solved using, for example, a [[least squares]] method or an [[Kalman filter#Extended Kalman filter|extended Kalman filter]].
| |
| | |
| Additionally, the TDOA of multiple transmitted pulses from the emitter can be averaged to improve accuracy.
| |
| | |
| === Reciprocal case: locating a receiver from multiple transmitter sites ===
| |
| Multilateration can also be used by a single receiver to locate itself, by measuring the TDOA of signals emitted from three or more synchronised transmitters at known locations. This can be used by navigation systems, an example being the British [[Decca Navigator System|DECCA]] navigation system, developed during World War II, which used the [[Phase (waves)|phase]]-difference of two transmitters, rather than the TDOA of a pulse, to define the hyperboloids. This allowed the transmitters to broadcast a continuous wave signal. Phase-difference and time-difference can be considered the same for narrow-band transmitters.
| |
| | |
| ==TDOA geometry==
| |
| | |
| [[Image:TDOA Geometry.png|thumb|Left|350px|alt=Geometry of spherical waves emanating from an emitter and passing through several receivers.|Fig 2. TDOA geometry.]]
| |
| | |
| Consider an emitter (E in Figure 2) at an unknown location vector
| |
| :E = (x, y, z)
| |
| | |
| which we wish to locate. The source is within range of N+1 receivers at known locations
| |
| :P<sub>0</sub>, P<sub>1</sub>, ..., P<sub>m</sub>, ..., P<sub>N</sub>.
| |
| The subscript m refers to any one of the receivers:
| |
| :P<sub>m</sub> = (x<sub>m</sub>, y<sub>m</sub>, z<sub>m</sub>)
| |
| :0 ≤ m ≤ N
| |
| | |
| The distance (R) from the emitter to one of the receivers in terms of the coordinates is
| |
| :{{NumBlk|:|
| |
| <math>
| |
| R_m = \left | \vec P_{m} - \vec E \right | = \sqrt{(x_m-x)^2+(y_m-y)^2+(z_m-z)^2}
| |
| </math>
| |
| |{{EquationRef|1}}}}
| |
| | |
| The math is made easier by placing the origin at one of the receivers (P<sub>0</sub>), which makes its distance to the emitter
| |
| :{{NumBlk|:|
| |
| <math>
| |
| R_0 = \sqrt{x^2+y^2+z^2}
| |
| </math>
| |
| |{{EquationRef|2}}}}
| |
| | |
| ==Measuring the time difference in a TDOA System==
| |
| | |
| {{multiple image
| |
| | width = 350
| |
| | direction = vertical
| |
| | footer = Examples of measuring time difference with cross correlation.
| |
| | image1 = CrossCorr_Pulse.png
| |
| | caption1 = Fig 3a. Pulse Signal
| |
| | image2 = CrossCorr_Wideband.png
| |
| | caption2 = Fig 3b. Wide-band Signal
| |
| | image3 = CrossCorr_Narrowband.png
| |
| | caption3 = Fig 3c. Narrow-band Signal
| |
| }}
| |
| | |
| The distance <math> R_m </math> in equation {{EquationNote|1}} is the wave speed (<math>v</math>) times transit time (<math>T_m</math>). A TDOA multilateration system measures the time difference (<math> \tau_m</math>) of a wavefront touching each receiver. The TDOA equation for receivers m and 0 is
| |
| :{{NumBlk|:|
| |
| <math>\begin{align}
| |
| v \tau_m & = v T_m - v T_0 \\
| |
| v \tau_m & = R_m - R_0
| |
| \end{align}</math>
| |
| |{{EquationRef|3}}}}
| |
| | |
| Figure 3a is a simulation of the a pulse waveform recorded by receivers <math>P_0</math> and <math>P_1</math>. The spacing between <math>E</math>, <math>P_1</math> and <math>P_0</math> is such that the pulse takes 5 time units longer to reach <math>P_1</math> than <math>P_0</math>. The units of time in Figure 3 are arbitrary. The following table gives approximate time scale units for recording different types of waves.
| |
| | |
| {| class="wikitable" style="text-align:left;
| |
| |-
| |
| ! Type of wave
| |
| ! Material
| |
| ! Time Units
| |
| |-
| |
| | Acoustic
| |
| | Air
| |
| | 1 millisecond
| |
| |-
| |
| | Acoustic
| |
| | Water
| |
| | 1/2 millisecond
| |
| |-
| |
| | Acoustic
| |
| | Rock
| |
| | 1/10 millisecond
| |
| |-
| |
| | Electromagnetic
| |
| | Vacuum, air
| |
| | 1 nanosecond
| |
| |}
| |
| | |
| The red curve in Figure 3a is the [[cross-correlation]] function <math>(P_1 \star P_0)</math>. The cross correlation function slides one curve in time across the other and returns a peak value when the curve shapes match. The peak at time = 5 is a measure of the time shift between the recorded waveforms, which is also the <math>\tau</math> value needed for Equation {{EquationNote|3}}.
| |
| | |
| Figure 3b is the same type of simulation for a wide-band waveform from the emitter. The time shift is 5 time units because the geometry and wave speed is the same as the Figure 3a example. Again, the peak in the cross correlation occurs at <math>\tau_1 = 5</math>.
| |
| | |
| Figure 3c is an example of a continuous, narrow-band waveform from the emitter. The cross correlation function shows an important factor when choosing the receiver geometry. There is a peak at Time = 5 plus every increment of the waveform period. To get one solution for the measured time difference, the largest space between any two receivers must be closer than one wavelength of the emitter signal. Some systems, such as the [[LORAN-C|LORAN C]] and [[Decca Navigator System|Decca]] mentioned at earlier (recall the same math works for moving receiver & multiple known transmitters), use spacing larger than 1 wavelength and include equipment, such as a [[Phase detector|Phase Detector]], to count the number of cycles that pass by as the emitter moves. This only works for continuous, narrow-band waveforms because of the relation between phase (<math>\theta</math>), frequency (f) and time (T)
| |
| :<math>\theta = 2 \pi f \cdot T</math>. | |
| The phase detector will see variations in frequency as measured [[phase noise]], which will be an [[Propagation of uncertainty|uncertainty that propagates]] into the calculated location. If the phase noise is large enough, the phase detector can become unstable.
| |
| | |
| ==3-D Solution==
| |
| Equation {{EquationNote|3}} is the hyperboloid described in the previous section, where 4 receivers (0 ≤ m ≤ 3) lead to 3 non-linear equations in 3 unknown values (x,y,z). The system must then solve for the unknown emitter location in real time.
| |
| | |
| Civilian air traffic control multilateration systems use the Mode C [[Secondary surveillance radar|SSR]] transponder return to find the altitude (z). Three or more receivers at known locations are used to find the other 2 dimensions (x, y).
| |
| | |
| R. Bucher and D. Misra show the detailed algebra to locate 1 transmitter with TDOA between 4 receivers.<ref name="Bucher Misra">[http://www.hindawi.com/journals/vlsi/2002/935925.abs.html A Synthesizable VHDL Model of the Exact Solution for Three-dimensional Hyperbolic Positioning System], Ralph Bucher and D. Misra, VLSI Design, Volume 15 (2002), Issue 2, Pages 507-520.</ref> Their solution is a set of linear equations to find (x, y) and a quadratic for (z).
| |
| | |
| Improving accuracy with a large number of receivers can be a problem for devices with small embedded processors because of the time required to solve several simultaneous, non-linear equations ({{EquationNote|1}}, {{EquationNote|2}} & {{EquationNote|3}}). The TDOA problem can be turned into a system of linear equations when there are 3 or more receivers, which can reduce the computation time. Starting with equation {{EquationNote|3}}, solve for R<sub>m</sub>, square both sides, collect terms and divide all terms by <math> v \tau_m </math>:
| |
| :{{NumBlk|:|
| |
| <math>\begin{align}
| |
| R_{m}^2 & = (v \tau_{m} + R_{0} )^2 \\
| |
| R_{m}^2 & = (v \tau_{m})^2 + 2 v \tau_{m} R_{0} + R_{0}^2 \\
| |
| 0 & = (v \tau_{m})^2 + 2 v \tau_{m} R_{0} + R_{0}^2 - R_{m}^2 \\
| |
| 0 & = (v \tau_{m}) + 2 R_{0} + \frac {(R_{0}^2 - R_{m}^2)} {v \tau_{m}}. \\
| |
| \end{align} </math>
| |
| |{{EquationRef|4}}}}
| |
| | |
| Removing the 2 R<sub>0</sub> term will eliminate all the square root terms. That is done by subtracting the TDOA equation of receiver m = 1 from each of the others (2 ≤ m ≤ N)
| |
| :{{NumBlk|:|
| |
| <math> \begin{align}
| |
| & 0 = v \tau_{m} + 2 R_{0} + \frac {(R_{0}^2 - R_{m}^2)} {v \tau_{m}} \\
| |
| & 0 = -v \tau_{1} - 2 R_{0} - \frac {(R_{0}^2 - R_{1}^2)} {v \tau_{1}} \\
| |
| | |
| \hline \\
| |
| | |
| & 0 = v \tau_{m} - v \tau_{1} + \frac {(R_{0}^2 - R_{m}^2)} {v \tau_{m}}
| |
| - \frac {(R_{0}^2 - R_{1}^2)} {v \tau_{1}}.
| |
| \end{align} </math>
| |
| |{{EquationRef|5}}}}
| |
| | |
| Focus for a moment on equation {{EquationNote|1}}. Square R<sub>m</sub>, group similar terms and use equation {{EquationNote|2}} to replace some of the terms with R<sub>0</sub>.
| |
| :{{NumBlk|:|
| |
| <math> \begin{align}
| |
| R_m^2 & = x_m^2+y_m^2+z_m^2 -x\ 2x_m - y\ 2y_m - z\ 2z_m +x^2 + y^2 + z^2 \\
| |
| & = x_m^2+y_m^2+z_m^2 -x\ 2x_m - y\ 2y_m - z\ 2z_m + R_0^2 \\
| |
| R_0^2 - R_m^2 & = -x_m^2-y_m^2-z_m^2 +x\ 2x_m + y\ 2y_m + z\ 2z_m.
| |
| \end{align} </math>
| |
| |{{EquationRef|6}}}}
| |
| | |
| Combine equations {{EquationNote|5}} and {{EquationNote|6}}, and write as a set of linear equations of the unknown emitter location x,y,z
| |
| :{{NumBlk|:|
| |
| <math> \begin{align}
| |
| 0 & = xA_m + yB_m + zC_m + D_m \\
| |
| & A_m = \frac {2x_m} {v \tau_m} - \frac {2x_1} {v \tau_1} \\
| |
| & B_m = \frac {2y_m} {v \tau_m} - \frac {2y_1} {v \tau_1} \\
| |
| & C_m = \frac {2z_m} {v \tau_m} - \frac {2z_1} {v \tau_1} \\
| |
| | |
| & D_m = v \tau_m - v \tau_1 - \frac {x_m^2 + y_m^2 + z_m^2} {v \tau_m}
| |
| + \frac {x_1^2 + y_1^2 + z_1^2} {v \tau_1}.
| |
| \end{align} </math>
| |
| |{{EquationRef|7}}}}
| |
| | |
| Use equation {{EquationNote|7}} to generate the four constants <math>A_m,B_m,C_m,D_m</math> from measured distances and time for each receiver 2 ≤ m ≤ N. This will be a set of N [[inhomogeneous linear equation]]s.
| |
| | |
| There are many robust linear algebra methods that can solve for the values of (x,y,z), such as [[Gaussian elimination|Gaussian Elimination]]. Chapter 15 in [[Numerical Recipes]]<ref>[http://www.numerical-recipes.com/ Numerical Recipes] official website</ref> describes several methods to solve linear equations and estimate the uncertainty of the resulting values.
| |
| | |
| == 2-D Solution ==
| |
| Finding the emitter location in a two dimensional geometry can use any of the methods used for the 3-D geometry. The coordinate frame is typically defined to make the z dimension zero or constant. Examples of 2-D multilateration are [[Short-wave radio|short wave radio]] long distance communications through the Earth's atmosphere, acoustic wave propagation in the [[SOFAR channel|sound fixing and ranging channel]] of the oceans and the [[Loran C|LORAN]] navigation system.
| |
| | |
| == Accuracy ==
| |
| | |
| For [[trilateration]] or multilateration, calculation is done based on distances, which requires the frequency and the wave count of a received transmission. For [[triangulation]] or multiangulation, calculation is done based on angles, which requires the phases of received transmission plus the wave count.
| |
| | |
| For lateration compared to angulation, the numerical problems compare, but the technical problem is more challenging with angular measurements, as angles require two measures per position when using optical or electronic means for measuring phase differences instead of counting wave cycles.
| |
| | |
| Trilateration in general is calculating with triangles of known distances/sizes, mathematically a very sound system. In a triangle, the [[angle]]s can be derived if one knows the length of all sides, (see [[Congruence_(geometry)#Congruence of triangles|congruence]]), but the length of the sides cannot be derived based on all of the angles, not without knowing the length of at least one of the sides (a baseline) (see [[Similarity_(geometry)#Angle/side similarities|similarity]]).
| |
| | |
| In 3D, when four or more angles are in play, locations can be calculated from ''n'' + 1 = 4 measured angles plus one known baseline or from just ''n'' + 1 = 4 measured sides.
| |
| | |
| Multilateration is, in general, far more accurate for locating an object than sparse approaches such as [[trilateration]], where with planar problems just three distances are known and computed. Multilateration serves for several aspects:
| |
| * over-determination of an ''n''-variable [[Quadratic equation|quadratic problem]] with (''n'' + 1) + ''m'' [[quadratic equation]]s
| |
| * [[Stochastics|stochastic]] errors prohibiting a [[Determinism|deterministic]] approach to solving the equations
| |
| * [[Cluster analysis|clustering]] needs to segregate members of various clusters contributing to various [[Model theory|models]] of solving, i.e. fixed locations, [[Oscillation|oscillating]] locations and moving locations
| |
| | |
| Accuracy of multilateration is a function of several variables, including:
| |
| * The [[Antenna (radio)|antenna]] or [[sensor]] [[geometry]] of the [[Receiver (radio)|receiver(s)]] and [[Transmitter|transmitter(s)]] for [[Electronics|electronic]] or [[Optics|optical]] transmission.
| |
| * The timing accuracy of the receiver system, i.e. [[Drift (telecommunication)|thermal stability]] of the clocking [[oscillators]].
| |
| * The accuracy of frequency synchronisation of the transmitter [[oscillators]] with the receiver oscillators.
| |
| * Phase synchronisation of the transmitted signal with the received signal, as propagation effects as e.g. diffraction or reflection changes the phase of the signal thus indication deviation from [[Line-of-sight propagation|line of sight]], i.e. multipath reflections.
| |
| * The bandwidth of the emitted pulse(s) and thus the rise-time of the pulses with pulse coded signals in transmission.
| |
| * Inaccuracies in the locations of the transmitters or receivers when used as a known location
| |
| The accuracy can be calculated by using the [[Cramér–Rao bound]] and taking account of the above factors in its formulation.
| |
| | |
| == Example applications ==
| |
| | |
| * [[Sound ranging]] – Using sound to locate artillery fire.
| |
| * [[Electronic Targets]] – Using the Mach wave of a bullet passing a sensor array to determine the point of arrival of the bullet on a firing range target.
| |
| * [[Decca Navigator System]] – A system used from the end of World War II to the year 2000, employing the phase-difference of multiple transmitters to locate on the intersection of hyperboloids
| |
| * [[OMEGA Navigation System]] – A worldwide system similar to Decca, shut down in 1997
| |
| * [[Radio navigation#Gee|GEE]] – British aircraft location technique from World War II, using accurate reference transmitters
| |
| * [[LORAN-C]] – navigation system using TDOA of signals from multiple synchronised transmitters
| |
| * Passive [[Electronic warfare support measures|ESM]] multilateration systems, including [[Kopáč passive sensor|Kopáč]], [[Ramona passive sensor|Ramona]], [[Tamara passive sensor|Tamara]], [[VERA passive sensor|VERA]] and possibly [[Kolchuga passive sensor|Kolchuga]] – location of a transmitter using multiple receivers
| |
| * [[Mobile phone tracking]] – using multiple base stations to estimate phone location (by either the phone itself, or the phone network)
| |
| * [[Reduced Vertical Separation Minima]] (RVSM) monitoring using [[Secondary Surveillance Radar]] – Mode C/S transponder replies to calculate the position of an aircraft. Application to RVSM was first demonstrated by [[Roke Manor Research Limited]] in 1989. <ref>{{cite web|title=Perfect Timing|url=http://www.roke.co.uk/resources/articles/Perfect-Timing-ATTI-2002.pdf|accessdate=31 August 2012|author=Air Traffic Technology International|year=2002}}</ref>
| |
| | |
| == Simplification ==
| |
| | |
| For applications where no need for absolute coordinates determination is assessed, the implementing of a more simple solution is advantageous. Compared to multilateration as the concept of ''crisp locating'', the other option is ''[[fuzzy locating]]'', where just one distance delivers the relation between detector and detected object. This most simple approach is unilateration. However, such unilateration approach never delivers the angular position with reference to the detector. Many solutions are available today.
| |
| Some of these vendors offer a position estimate based on combining several laterations. This approach is often not stable, when the wireless ambience is affected by metal or water masses. Other vendors offer room discrimination with a room-wise excitation, one vendor offers a position discrimination with a contiguity excitation.
| |
| | |
| == See also ==
| |
| | |
| * [[Hyperbolic navigation]]
| |
| * [[FDOA]] Frequency difference of arrival. Analogous to TDOA using differential doppler.
| |
| * [[Triangulation]] – Location by angular measurement on lines of bearing that intersect
| |
| * [[Trilateration]] – Location by distance (e.g. time-of-flight) measurement on coincident signals from multiple transmitters.
| |
| * [[Mobile phone tracking]] – used in [[GSM]] networks
| |
| * [[Multidimensional scaling]]
| |
| * [[Radiolocation]]
| |
| * [[Radio navigation]]
| |
| * [[Real-time locating]] – International standard for asset and staff tracking using wireless hardware and real-time software
| |
| * [[Real time location system]] – General techniques for asset and staff tracking using wireless hardware and real-time software
| |
| | |
| == Notes ==
| |
| {{reflist|2}}
| |
| | |
| == References ==
| |
| {{refbegin}}
| |
| * [http://www.multilateration.com The Multilateration Executive Reference Guide] is an easy-to-read reference for air traffic management, airport and airline professionals to learn more about this next-generation surveillance technology
| |
| {{refend}}
| |
| | |
| [[Category:Elementary geometry]]
| |
| [[Category:Euclidean geometry]]
| |
| [[Category:Radio navigation]]
| |
| [[Category:Geolocation]]
| |
| [[Category:Ubiquitous computing]]
| |
| [[Category:Wireless locating]]
| |