Wireless Communications and Mobile Computing

Volume 2019, Article ID 2707469, 8 pages

https://doi.org/10.1155/2019/2707469

## Probability Weighting Localization Algorithm Based on NLOS Identification in Wireless Network

Information Science and Engineering College, Chongqing Jiaotong University, No. 66 Xuefu Road, Nan’an Dist, Chongqing 400074, China

Correspondence should be addressed to Shixun Wu; moc.361@333nuxihsuw

Received 22 November 2018; Revised 17 February 2019; Accepted 4 March 2019; Published 28 March 2019

Guest Editor: Mohamed Laaraiedh

Copyright © 2019 Shixun Wu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In this paper, a localization scenario that the home base station (BS) measures time of arrival (TOA) and angle of arrival (AOA) while the neighboring BSs only measure TOA is investigated. In order to reduce the effect of non-line of sight (NLOS) propagation, the probability weighting localization algorithm based on NLOS identification is proposed. The proposed algorithm divides these range and angle measurements into different combinations. For each combination, a statistic whose distribution is chi-square in LOS propagation is constructed, and the corresponding theoretic threshold is derived to identify each combination whether it is LOS or NLOS propagation. Further, if those combinations are decided as LOS propagation, the corresponding probabilities are derived to weigh the accepted combinations. Simulation results demonstrate that our proposed algorithm can provide better performance than conventional algorithms in different NLOS environments. In addition, computational complexity of our proposed algorithm is analyzed and compared.

#### 1. Introduction

Wireless localization which can determine the position of mobile station (MS) in wireless network has received considerable attention over the past years, especially the application of the location based services (LBSs). The existing wireless localization techniques such as received signal strength (RSS) [1], time of arrival (TOA) [2, 3], time difference of arrival (TDOA) [4, 5], angle of arrival (AOA) [6–9], and the combination of the above one are often used in wireless network. The non-line of sight (NLOS) propagation is one of the dominant factors to affect the localization accuracy of MS, and it happens when the direct signal path between MS and base station (BS) is blocked. Comparing with the line of sight (LOS) propagation, the signal travels extra distance, inducing power loss and angle bias in NLOS propagations.

There are two ways to cope with the NLOS condition. The first way localizes with all NLOS and LOS measurements, but provides weighting, nonlinear optimization, or scaling to minimize the effects of the NLOS error. Residual weighting algorithm (Rwgh) [10] is very effective in reducing the NLOS error for TOA-based localization system. Its main idea is to divide the range measurements into different combinations, each combination obtains the intermediate position estimate of MS with nonlinear least square (NLS) algorithm, and the final position estimate of MS is weighted by the intermediate position estimate and the corresponding normalized residual. However, it does not discard any combination which may be corrupted greatly by the NLOS propagation and also has high computational complexity when the number of involved BSs is big. As we know, if the home BS is equipped with antenna array, AOA measurement is obtainable, and it is helpful to improve the localization accuracy of TOA-based wireless network. In [11], Geometric Dilution of Precision (GDOP) is introduced into hybrid TOA/AOA measurements to propose GDOP-weighted localization algorithm. A nonlinear constrained optimization algorithm whose constraints on range and angle are inferred from geometry with hybrid TOA/AOA measurements is proposed in [12]. Based on a single bounce scattering environment, a joint TOA/AOA constrained minimization method which incorporates the unknown scatterers into the nonlinear optimization model is proposed in [13, 14]. By introducing scale factors to build the relation between the true distances and measured distances, the work in [15] proposes a geometric method to locate MS with only two BSs. Taylor series least square (TS-LS) algorithm which is developed for TOA-based systems to incorporate AOA measurements is proposed in [16]. However, all the nonlinear optimization algorithms or scaling algorithms with hybrid TOA/AOA measurements have a good localization accuracy at the cost of computational complexity. The second way attempts to identify and localize with the LOS BSs. NLOS identification is done with a time-history based hypothesis test [17, 18], the feature of channel statistics [19–21], or a residual test which compares the residuals of a group against a predetermined threshold [22, 23]. However, the selection of a predetermined threshold is obtained by experience.

In this paper, we investigate hybrid TOA/AOA NLOS identification with a residual test and the weighting localization approach to minimize the effect of NLOS error. Different from the residual test in [22, 23], the selection of a predetermined threshold is based on theoretical analysis. Moreover, different from the Rwgh algorithm in [10], AOA measurement from home BS is introduced, and the weight of each combination is the corresponding probability rather than the corresponding normalized residual. Specifically, we divide the range and angle measurements into different combinations. For each combination, the position estimate of MS and its corresponding covariance matrix are computed by linear least square algorithm. Then, the range errors from different BSs are easily obtained with the position estimate of MS, and a statistic is derived and constructed by utilizing the first-order Taylor series. If the combination is derived from LOS measurements, this statistic obeys the chi-square distribution; otherwise it is not. Thus, given a false alarm probability, a theoretical threshold whose value can be computed from chi-square probability density function is derived to identify the combination whether it is LOS or NLOS. If the statistic is larger than the threshold, it is rejected; otherwise it is retained. Next, for those retained combinations, the corresponding probabilities are used to weigh the intermediate position estimates. Simulation results show that the proposed probability weighting localization algorithm based on NLOS identification has better performance than the existing algorithms in different NLOS environments.

The rest of the paper is described as follows. In Section 2, the system model is presented. In Section 3, the probability weighting localization algorithm is proposed. Section 4 presents the simulation results and computational complexity. Finally, Section 5 provides some conclusions.

#### 2. System Model

There are BSs available to localize MS in wireless network, without loss of generality; we suppose that is the home BS which can obtain range and angle measurements, while the neighboring only has range measurement [10, 12]. The system model is described aswhere is the measured TOA between MS and , is the speed of light, the position of MS is , and is the position of . and are the NLOS error and measurement noise, respectively. Measurement noise is a Gaussian distribution . If has a LOS path to MS, then . is the range measurement from i-th BS. is the angle measurement in home BS. is the function of inverse tangent and the angle measurement error is the sum of the angle measurement noise and NLOS angle deviation. If the home BS experiences LOS propagation, the angle measurement noise is a Gaussian distribution . Due to the obtainable parameters about range and angle measurement noise in [24], we assume that the variances of range and angle measurement noise are known, whereas the NLOS errors are unknown in this paper.

#### 3. Probability Weighting Localization Algorithm

In this section, we present the proposed NLOS identification and probability weighting localization algorithm in terms of system model shown in Section 2. The proposed algorithm contains two steps: NLOS identification and probability weighting.

##### 3.1. NLOS Identification

As shown in Section 2, range measurements and one angle measurement are available to localize MS. Generally speaking, two range measurements and one angle measurement can provide the position estimate of MS. Moreover, it was explained that the measurements in the serving BS are more reliable than ones in the neighboring BSs [12]. Thus, we can divide the range and angle measurements into different combinations; each combination must contain the measurements in serving BS. BSs have combinations. In LOS environment, the model shown in (1) can be described as

For AOA equation in (2), using the fact that when , we have the following geometrical relationship [16, 25]:

Squaring the range equations in (2), two noise terms are present. Since the measurement noise is relatively small, the square term of noise is assumed to be negligible in comparison with the first-order term of noise. Thus, we have the following approximate equations:

By fixing the first equation in (4) as the reference, subtracting it from the rest of equations and combining them with equation (3), we can obtain the following linear equations:where

The position estimate of MS can be determined by solving (5) using least square algorithm,where , , is the transpose of matrix , and is the inverse matrix of . Due to the unknown actual distance , we can use the range measurement to replace them for estimating .

We choose as the initial estimate and use the first-order Taylor Series expansion to approximate the nonlinear range equations in (2). Then we obtainwhere

Putting (7) into (8) can obtain the error vector as follows:where .

Due to the Gaussian distribution of and in , the error vector is a Gaussian distribution with mean , covariance matrix .where

In order to validate (11), we define the following hypotheses and alternatives:

The hypothesis holds true, if the BSs from the combination are LOS. The alternative holds if at least one BS from the combination is corrupted by NLOS errors. If the error vector is computed from LOS BSs, we can construct the test statistic and compare it with the threshold , based on preset false alarm probability under assumption (13). The distribution of under is the chi-square distribution with N degree of freedom [26]. If is larger than the threshold , the hypothesis is rejected; otherwise it is retained. Then, the false alarm probability expressed as the probability to decide NLOS combination if this combination is derived from LOS measurements is defined as follows:where is the chi-square probability density function with N degrees of freedom and is the gamma function.

In the above discussion, we only consider the combination that all the BSs are involved. It is easily extended to other combinations. For example, if the combination contains two range measurements and one angle measurement, there are combinations with BSs. These combinations have the same threshold whose value can be computed from (14) based on the chi-square probability density function with 2 degrees of freedom. By constructing the corresponding matrix or vectors , , , , , , and shown above with two range measurements and one angle measurement, each combination can obtain the position estimate of MS, the error vector , and covariance matrix from (7), (8), and (11), respectively. Then a test statistic is computed and compared with the threshold to decide whether this combination is LOS or NLOS.

##### 3.2. Probability Weighting

For our system model in Section 2, there are combinations. The NLOS identification is performed for each combination and the accepted ones are weighted with different probabilities. Note that the accepted combinations are labeled as , where is the number of the accepted combinations. We assume that the corresponding error vector, covariance matrix, and test statistic of each accepted combination are denoted as , , and , respectively. As we know, the smaller the value of is, the bigger probability the combination is LOS. Therefore, the probabilities of each accepted combination can be obtained approximately as . To ensure that the sum of probabilities is one, we normalize as

The final position estimate of MS is weighted aswhere is the intermediate position estimate of the j-th accepted combination.

In extreme circumstances, none of these combinations is accepted in the tests; the proposed algorithm will not output a valid position estimate of MS. If this situation happened, it means that the range measurements in neighboring BSs deteriorate significantly. The localization accuracy will be degraded if they are combined with the range and angle measurements in home BS. Thus, only home BS is reliable to provide the position estimate of MS. With the assumption of LOS propagation, the position estimate of MS is easily obtained as .

#### 4. Simulation Results

In this section, we carry out some simulations to prove the performance of the proposed NLOS identification and probability weighting localization algorithm. Three BSs with a hexagonal layout shown in Figure 1 are deployed. Without loss of generality, we assume that the position of is , being the home BS. Because the radius of hexagon is 1000m, the position of and can be easily obtained as and , respectively. The position of MS is .