Definitions

# Multilateration

Multilateration, also known as hyperbolic positioning, is the process of locating an object by accurately computing the time difference of arrival (TDOA) of a signal emitted from the object to three or more receivers. It also refers to the case of locating a receiver by measuring the TDOA of a signal transmitted from three or more synchronised transmitters.

Multilateration should not be confused with trilateration, which uses distances or absolute measurements of time-of-flight from three or more sites, or with triangulation, which uses a baseline and at least two angles measured e.g. with receiver antenna diversity and phase comparison.

## 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 series of emitter locations would give the same measurement of TDOA. Given two receiver locations and a known TDOA, the locus of possible emitter locations is a one half of a two-sheeted hyperboloid.

In simple terms, with two receivers at known locations, an emitter can be located onto a hyperboloid. 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 problem and solved using, for example, a least squares method or an 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 navigation system, developed during World War II, which used the 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.

## Derivation

Consider an emitter at unknown location $\left(x, y, z\right)$ which we wish to locate. Consider also a multilateration system comprising four receiver sites at known locations: a central site, C, a left site, L, a right site, R and a fourth site, Q.

The travel time (T) of pulses from the emitter at ($x, y, z$) to each of the receiver locations is simply the distance divided by the pulse propagation rate (c):

$T_L=frac\left\{1\right\}\left\{c\right\}left\left(sqrt\left\{\left(x-x_L\right)^2+\left(y-y_L\right)^2+\left(z-z_L\right)^2\right\}right\right)$
$T_R=frac\left\{1\right\}\left\{c\right\}left\left(sqrt\left\{\left(x-x_R\right)^2+\left(y-y_R\right)^2+\left(z-z_R\right)^2\right\}right\right)$
$T_Q=frac\left\{1\right\}\left\{c\right\}left\left(sqrt\left\{\left(x-x_Q\right)^2+\left(y-y_Q\right)^2+\left(z-z_Q\right)^2\right\}right\right)$
$T_C=frac\left\{1\right\}\left\{c\right\}left\left(sqrt\left\{\left(x-x_C\right)^2+\left(y-y_C\right)^2+\left(z-z_C\right)^2\right\}right\right)$
If the site C is taken to be at the coordinate system origin,
$T_C=frac\left\{1\right\}\left\{c\right\}left\left(sqrt\left\{x^2+y^2+z^2\right\}right\right)$

Then the time difference of arrival between pulses arriving directly at the central site and those coming via the side sites can be shown to be:

$tau_L=T_L-T_C=frac\left\{1\right\}\left\{c\right\}left\left(sqrt\left\{\left(x-x_L\right)^2+\left(y-y_L\right)^2+\left(z-z_L\right)^2\right\}-sqrt\left\{x^2+y^2+z^2\right\}right\right)$
$tau_R=T_R-T_C=frac\left\{1\right\}\left\{c\right\}left\left(sqrt\left\{\left(x-x_R\right)^2+\left(y-y_R\right)^2+\left(z-z_R\right)^2\right\}-sqrt\left\{x^2+y^2+z^2\right\}right\right)$
$tau_Q=T_Q-T_C=frac\left\{1\right\}\left\{c\right\}left\left(sqrt\left\{\left(x-x_Q\right)^2+\left(y-y_Q\right)^2+\left(z-z_Q\right)^2\right\}-sqrt\left\{x^2+y^2+z^2\right\}right\right)$
where $\left(x_L, y_L, z_L\right)$ is the location of the left receiver site, etc, and $c$ is the speed of propagation of the pulse, often the speed of light. Each equation defines a separate hyperboloid.

The multilateration system must then solve for the unknown target location $\left(x, y, z\right)$ in real time. All the other symbols are known.

Note that in civilian air traffic control multilateration systems, the unknown height, $z$, is often directly derived from the Mode C SSR transponder return. In this case, only three sites are required for a 3D solution.

## Multilateration 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 with known distances/sizes, mathematical a very sound system. In a triangle if the length of three sides is known, the angles can be calculated (congruence), however for calculating the other sides from angles at least one side (baseline) must be known additionally (otherwise just 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:

Accuracy of multilateration is a function of several variables, including:

• The antenna or sensor geometry of the receiver(s) and transmitter(s) for electronic or optical transmission.
• The timing accuracy of the receiver system, i.e. 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, 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.