# Hybrid TOA/AOA Approximate Maximum Likelihood Mobile Localization

Authors: Mohamed Zhaounia; Mohamed Adnan Landolsi; Ridha Bouallegue
Journal: Journal of Electrical and Computer Engineering
 (2010)

---

## Abstract

This letter deals with a hybrid time-of-arrival/angle-of-arrival (TOA/AOA) approximate maximum likelihood (AML) wireless location algorithm. Thanks to the use of both TOA/AOA measurements, the proposed technique can rely on two base stations (BS) only and achieves better performance compared to the original approximate maximum likelihood (AML) method. The use of two BSs is an important advantage in wireless cellular communication systems because it avoids hearability problems and reduces network signaling burden. Simulation results show that, for certain scenarios, the proposed hybrid TOA/AOA AML with two BSs can outperform the AML with up to six BSs.


---

## Body

## 1. Introduction

Mobile positioning in wireless communication has gained a wide acceptance and trigged a large number of innovative applications [1, 2]. Various techniques have been proposed for mobile positioning, depending on the collected measurements by the wireless networking infrastructure. Broadly, these techniques can be classified as (i) time-based location that produces a circle centered at the base station (BS) and the mobile station (MS) location estimate is determined by the intersection of three circles at least, and (ii) angle-based location, which requires BSs equipped with antenna arrays, defines a line of bearing from the BS to the MS whose position is calculated from the intersection of a minimum of two BSs.Because of the noisy measurements, the mobile position is typically determined by some statistical estimation which processes a set of nonlinear equations constructed from the TOA and/or AOA statistics with knowledge of the BSs geometry. Among the different techniques proposed, one particularly interesting approach known as approximate maximum likelihood (AML) was developed in [3] and shown to achieve near-optimal performance without the complexity of “full” maximum likelihood estimation. One limitation of the original AML method in [3] is that it relies on TOA data only and needs measurements from several BSs (all assumed to have the same measurement error variance). However due to hearability problem, the signal received at BSs other than the serving BS may not be strong enough to provide reliable measurements [4]. Furthermore, the use of a large number of BSs results in a network overload by requiring complex messaging and synchronization. As such, obtaining reliable measurements from at least three BSs becomes difficult to achieve in real-world scenarios. In this letter, we present a new “hybrid” version of the AML algorithm by incorporating both TOA and AOA data and show that this new positioning method can have superior performance by using measurements from only 2 BSs.


## 2. The AML Approach

This section gives a brief description of the original AML approach in [3]. Assuming N≥3 BSs involved in the positioning process, the vector of measured distances is(1)Δ=[δ1⋯δN]T,
where δi=cti with ti being the measured TOA of the ith BS and c stands for the light velocity.The measured TOAdata are assumed to be corrupted by independent zero-mean Gaussian noise with a covariance matrix given by(2)Q=diag{varTOA1⋯varTOAN},
where the variances (varTOAi)1≤i≤N are not necessarily equals.The calculated vectorR is obtained from the estimated vector θ=[xy]T where (x,y) are the MS coordinates:(3)R(θ)=[r1⋯rN]T
where(4)ri2=s+ki-2xxi-2yyi,
with (xi,yi) being the coordinates of the ith BS and(5)s=x2+y2,(6)ki=xi2+yi2.Then, the probability density function (pdf) [5] of Δ given θ is expressed as(7)f(Δ∣θ)=(2π)-N/2det(Q)-1/2exp{-J2},
where(8)J=(Δ-R(θ))TQ-1(Δ-R(θ)).The ML estimation principle seeks to maximizef(Δ∣θ). Then, setting the gradient of J with respect to θ to zero gives two nonlinear ML equations, which are treated iteratively in a linear way in the simplified AML approach of [3]. Starting from some initial position, it first computes the coefficients of those equations expressed as a function of θ. Then applying (4) produces a quadratic in s. Next, a root selection routine (RSR) chooses the correct root. Finally, the AML selects the θ that gives the minimum J.


## 3. The Proposed HAML Approach

In the proposed approach the serving BS provides both the AOA and TOA measurements. By introducing the AOA information, we use statistics from two BSs only. Then, the vector of measured parameters becomes(9)ΔH=[δ1δ2β]T,
where β is the AOA at the home BS corrupted by an additive zero-mean Gaussian noise with variance varAOA independent from TOA noise. The covariance matrix becomes(10)QH=diag{varTOA1varTOA2varAOA}.
In the same way, the calculated vector includes the AOA measurement and is obtained as(11)RH(θ)=[r1r2α]T,
where α=arctan(y/x) and, without loss of generality, the serving BS is assumed at the origin (0.0). Then, the pdf of ΔH given θ is given by(12)f(ΔH∣θ)=(2π)-3/2det(QH)-1/2exp{-JH2},
where JH is expressed as(13)JH=∑i=12(δi-ri)2varTOAi+(β-α)2varAOA.Thus, the ML estimate is the vectorθ that minimizes JH.Setting the gradient ofJH with respect to θ to zero yields the following system of equations:(14)∑i=12(x-xi)(ri-δi)varTOAiri+(y-y1)(β-α)varAOAr12=0,∑i=12(y-yi)(ri-δi)varTOAiri-(x-x1)(β-α)varAOAr12=0.Substituting the quantityri-δi=(ri2-δi2)/(ri+δi), the system of equations above becomes(15)∑i=12pi(s+ki-δi2)=∑i=12pi(2xix+2yiy)-2w(y-y1),∑i=12qi(s+ki-δi2)=∑i=12qi(2xix+2yiy)+2w(x-x1),
where pi=(x-xi)/(varTOAi(ri+δi)), qi=(y-yi)/(varTOAi(ri+δi)) and w=1/2(β-α)/(varAOAr12).Then in matrix notation(16)2Aθ=b1s+b2,
where(17)A=[∑i=12pixi∑i=12piyi-w∑i=12qixi+w∑i=12qiyi],b1=[∑i=12pi∑i=12qi],b2=[∑i=12pi(ki-δi2)-2wy1∑i=12qi(ki-δi2)+2wx1].Performing the Least Squares (LS) estimator [5] leads to(18)x=θ(1)=12(A-1·b1)(1)s+12(A-1·b2)(1)=as+b,y=θ(2)=12(A-1·b1)(2)s+12(A-1·b2)(2)=a′s+b′.Substituting (18) into (5) a quadratic in s is obtained:(19)(as+b)2+(a′s+b′)2=s.Solving (19) raises 3 possible cases.Case 1.
The roots are real and there is only one positive root. Thens takes the value of this positive root and θ is estimated by applying (18).Case 2.
The roots are real and are both positive. Thenθ is calculated by applying (18) for each root and the corresponding value of JH is determined by performing (13). The correct root is the one that leads to the minimum value of JH and the estimated (x,y) is the corresponding θ.Case 3.
The roots are negative (or complex). In this case we take the absolute value (of the real part) and process as in the Case2.In order to calculate the initial estimate, we start with the noise-free case,(20)s+ki-δi2=2(xxi+yyi),(x-x1)sinβ-(y-y1)cosβ=0.
In matrix notation, this is equivalent to(21)Bθ=12c1s+c2,
where(22)B=[x1y1x2y2sinβ-cosβ],c1=[110],c2=[12(k1-δ12)12(k2-δ22)x1sinβ-y1cosβ].An initial position can be obtained from the LS estimate as(23)θ=12(BTB)-1BT{c1s+c2}.
Combining (23) and (5) provides a quadratic in s. Then, the adequate root is chosen similarly as discussed in (19).As can be seen and since it uses only 2 BSs, the proposed HAML algorithm is applicable regardless of the BSs’ geometry. In fact unlike the original AML algorithm, the HAML does not discriminate between cases where the BSs are on a straight line or not, which can be considered as an important advantage of the HAML regarding the implementation complexity of the AML in the case of 3 linear BSs [3].


## 4. Numerical Results

This section presents simulation results that illustrate the performance improvement of the HAML approach using 2 BSs only, compared to the original AML method using a larger number of BSs.In these simulations a hexagonal wireless cellular system layout with a radiusRc=4 Km is considered. TOA and AOA measurements are assumed to be corrupted by independent Gaussian noise with variances varTOA and varAOA, respectively. Due to range differences and geometrical distribution, we assume that the TOA standard deviation of the home BS is set to be 50 m and increases linearly with a step of 50 m for the other BSs without exceeding 200 m, that is, 100 m for the second BS, 150 m for the third BS, and 200 m for the other BSs. The AOA standard deviation is fixed to be 3∘. The statistical results are based on 5000 Monte-Carlo experiments.In the first simulation, we compare the performance of the proposed HAML technique used with 2 BSs to the original AML approach with 3 or more BSs. The MS position in m is at (1000,2000). Table1 lists the root mean squares (rms) of location error calculated as follows:(24)rms=(15000∑k=15000(x⌢k-x0)2+(y⌢k-y0)2)1/2,
where the superscript ()0 stands for the true value.Table 1
Performance comparison of 2-BS HAML and AML with various number of BSs.
HAMLAML# BS2 BSs3 BSs4 BSs5 BSs6 BSs7 BSsRMS (Km)0.09440.11170.10190.09910.09720.0937Although the AML technique uses a larger number of BSs, the HAML method shows better performance, as illustrated in Table1. The accuracy of AML surpasses that of HAML only if 7 BSs are involved in the positioning process which may not necessarily be available in a real-world scenario.In the second simulation, the mobile coordinates are varied on a linear trajectory, according to the expressionsx=1000+800*xMS and y=1000+800*yMS where xMS and yMS are two random variables uniformly distributed in [0,1]. The rms error is averaged over 100 ensemble runs. Figure 1 illustrates cumulative distribution function (CDF) of the error. As shown in this figure although the number of used BSs with the HAML technique is less than that of the AML approach, the performance of HAML is better than AML by providing a gain in accuracy of about 13%.Figure 1
Comparison between 2-BS HAML and 3-BS AML techniques.In the third simulation we study the influence of the variability of TOA variances between BSs, which is a real-world possibility in wireless cellular networks [6]. As shown in Figure 2, when we increase the difference in TOA variance between BSs, the AML performance becomes worse. However, the HAML is less influenced by the increase of this feature because of the contribution of the added AOA information, which in this case diminishes the influence of the TOA variances on the final accuracy by sharing the equations’ weighting coefficients.Figure 2
Influence of the variability of TOA variances between BSs on the position accuracy with 2-BS HAML and 3-BS AML techniques.The last simulation compares the performance of AML and HAML techniques against the AOA standard deviation. As depicted in Figure3, the accuracy of the mobile position estimate calculated with the HAML algorithm is enhanced when more accurate AOA measurements are used. However and as expected, no effect is introduced to the position estimate calculated with the AML technique since it does not use the AOA information. It is worth noticing that the performance of the proposed HAML algorithm still be better than that of the AML method using any AOA standard deviation.Figure 3
Influence of the AOA variances on the position accuracy with 2-BS HAML and 3-BS AML techniques.


## 5. Conclusion

In this letter a hybrid approximate maximum likelihood (HAML) approach for wireless location was proposed. The main contribution of this technique is that it can achieve good performance with TOA/AOA measurements from only two base stations, thereby avoiding problems of hearability and network overload. Simulation shows that the proposed technique is less influenced by the difference of noise powers between BSs. Moreover, in some geometric cases HAML results obtained from 2 BSs are as accurate as AML results obtained from 7 BSs.



---

*Source: 942657-2009-12-27.xml*