#### 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.

#### 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 BSs involved in the positioning process, the vector of measured distances is where with being the measured TOA of the th BS and 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 where the variances are not necessarily equals.

The calculated vector is obtained from the estimated vector where are the MS coordinates: where with being the coordinates of the th BS and

Then, the probability density function (pdf) [5] of given is expressed as where

The ML estimation principle seeks to maximize . Then, setting the gradient of 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 . Next, a root selection routine (RSR) chooses the correct root. Finally, the AML selects the that gives the minimum .

#### 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 where is the AOA at the home BS corrupted by an additive zero-mean Gaussian noise with variance independent from TOA noise. The covariance matrix becomes In the same way, the calculated vector includes the AOA measurement and is obtained as where and, without loss of generality, the serving BS is assumed at the origin (0.0). Then, the pdf of given is given by where is expressed as

Thus, the ML estimate is the vector that minimizes .

Setting the gradient of with respect to to zero yields the following system of equations:

Substituting the quantity , the system of equations above becomes where , and .

Then in matrix notation where

Performing the Least Squares (LS) estimator [5] leads to

Substituting (18) into (5) a quadratic in is obtained:

Solving (19) raises 3 possible cases.

*Case 1. *The roots are real and there is only one positive root. Then 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 is determined by performing (13). The correct root is the one that leads to the minimum value of and the estimated 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 Case 2.

In order to calculate the initial estimate, we start with the noise-free case, In matrix notation, this is equivalent to where

An initial position can be obtained from the LS estimate as Combining (23) and (5) provides a quadratic in . 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 radius Km is considered. TOA and AOA measurements are assumed to be corrupted by independent Gaussian noise with variances and 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 . 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). Table 1 lists the root mean squares (rms) of location error calculated as follows: where the superscript stands for the true value.

Although the AML technique uses a larger number of BSs, the HAML method shows better performance, as illustrated in Table 1. 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 expressions and where and are two random variables uniformly distributed in . 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%.

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.

The last simulation compares the performance of AML and HAML techniques against the AOA standard deviation. As depicted in Figure 3, 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.

#### 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.

#### Acknowledgment

Mohamed Adnan Landolsi is duly acknowledges the support of King Fahd University of Petroleum & Minerals.