EURASIPJournal on Wireless Communications and Networking

RESEARCH Open Access

TOA localization for multipath and NLOS environment with virtual stations

Deliang Liu1, Ying Wang2*, Peng He1, You Zhai1 and Hongyun Wang1

Abstract

An effective geometrical method which can locate a source based on the time-of-arrival (TOA) measurements and the floor plan information is proposed in simple non-line-of-sight (NLOS) indoor environment where the multiple-bound scattering signal paths can be ignored. The novelty of our work resides in converting the NLOS problem into line-of-sight (LOS) problem based on the virtual stations (VSs), so that the localization algorithm for LOS conditions can be applied in simple NLOS conditions.

The performance bound of the proposed method is also analyzed. Simulation results demonstrate the effectiveness of the proposed method.

Keywords: Localization, Time-of-arrival, Non-line-of-sight, Virtual station

1 Introduction

Global positioning system (GPS) has played an important role for tracking and positioning in outdoors. However, it is limited in indoor environments because the satellite signal is seriously disturbed and the positioning accuracy cannot meet the requirements of indoor localization [1-4]. Therefore, indoor localization has received a high amount of attention over the last years. Among various indoor localization methods, which differ by the type of information and system parameters used, the representative ones are the time-of-arrival (TOA), time-difference-of-arrival (TDOA), direction-of-arrival (DOA), received-signal-strength (RSS), and their combinations.

Nevertheless, NLOS and multipath radio propagation between base station (BS) and mobile station (MS), which can cause biased range estimates in TOA measurements, will dramatically degrade the performance of conventional TOA localization methods. To mitigate the NLOS impact, a variety of techniques and algorithms have been proposed in the literature. Filtering-based method is presented in [5] which only applies to the tracking of mobile tags. The method in [6] struggles to identify the NLOS paths, from a set of mixed measurements including both LOS and NLOS paths. However,

* Correspondence: msyingwang@outlook.com

2Tianjin University of Technology and Education, 300222 Tianjin, China Full list of author information is available at the end of the article

this approach requires LOS measurements involved and cannot be expected to be efficient when signals only propagate through the NLOS paths. In [7], TOA scattering models are adopted to alleviate the NLOS error, whereas the distributions of NLOS error may not be available in the practical environment. A robust scheme in [8] utilizes one-bound scattering and discards multiple-bound scattering to improve the accuracy of NLOS localization. But it needs bidirectional estimations of both AOA and TOA measurements, which means more hardware costs. The VS method has been used to solve the NLOS and multipath problem in [9-11]. Using the known room geometry, the reflections can be mapped onto virtual anchors, i.e., the known positions of virtual sources that can be used for estimating the position of the receiver. However, the methods in [9, 10] only resolved the multipath problem and only considered the reflections off the walls, and methods in [11] ignored the reflections and diffractions caused by obstacles as well as the effect of people in the room.

In our previous work [12], we proposed a geometrical approach using the VS method based on TOA and AOA measurements. But measuring AOA needs array antennas which will increase the hardware costs. In this paper, we present a VS method only based on the TOA measurements for the simple NLOS environments where the signals undergo no more than one-bound scattering. We assume that the NLOS multipath signals

Springer Open

© The Author(s). 2017 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.C International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, anc reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

undergo one-bound scattering (reflection of diffraction) at most, because multiple-bound scattering can be ignored as they suffer from severe fading and may be too weak to be observed [9]. We first develop a geometric TOA model for NLOS environments, establish the VSs, and derive the path estimation algorithm to estimate the realistic signal propagation paths as well as the possible positions of the unknown MS by the TOA measurements between every BS and MS. Then, we determine the final position of the MS based on the least squares principle. Finally, the performance bound of the proposed algorithm is analyzed and the simulations are run to demonstrate the effectiveness of the proposed algorithm. The paper is organized as follows: In Section 2, the proposed VS method is described in detail. Section 3 will provide the analytic performance bound of the proposed method. Section 4 analyzes the performance of our model by simulations. Section 5 concludes the research efforts.

2 The TOA-based VS algorithm

2.1 Motivation

A typical TOA-based localization algorithm for LOS conditions consists of two kinds of parameters: one is distance parameters which contain the measured distances between every BS and the position of the unknown MS; the other is position parameters which include the positions of BSs. In the following analysis, we will take the least squares (LS) algorithm as an example of localization algorithm for LOS conditions.

To measure the time-of-arrival between a BS (transmitter) and a MS (receiver), we usually detect the peak of the received signal arriving on the first detectable path, and the time delay t of the peak is the measurement of the TOA. The estimated distance between the transmitter and the receiver is d = t x c, in which c is the speed of radio wave propagation.

Under LOS conditions, the direct path is always the strongest path so that the amplitude is higher than the threshold. Therefore, the first detectable path is the direct path. However, paths arriving in close vicinity of the direct path will shift the peak of the first path causing modest positive ranging errors [13]. So, the measured range ¿los between a BS and a MS can be expressed as [14]:

¿LOS = d + e (1)

where d is the Euclidean distance between the anchor and the tag and e represents the Gaussian random noise with zero-mean and standard deviation a.

But under NLOS conditions, the direct path is blocked by objects situated between the transmitter and the receiver and may not be detected at all. Therefore, the measured range dNLOS can be given by:

dNLOS = d + b NLOS + e (2)

where bNLOS is the propagation delay imposed by NLOS. Because e can be neglected in comparison with bNLOS, the measured range is always larger than the true distance under NLOS condition

dNLOS > d (3)

As a result, if we use dNLOS as the distance parameters of conventional LS algorithm, the localization error will be larger in NLOS conditions. So, many efforts have been made to estimate the Euclidean distance d from the measured distance ¿NLOS in order to reduce the difference between them [7].

But the other kind of parameters has not been attached importance yet, and the positions of the BSs. Instead of using the locations of real BSs, we employ the positions of VSs. We still use the measured distance in NLOS condition as the distance parameters of LS algorithm. Then, if we can estimate which path (direct path, reflection path, or diffraction path) is the first detectable path between every BS and the MS, the VSs can be established so that reflections and diffractions can be treated equivalently as signals coming from the VSs through the direct paths from the MS's perspective. Therefore, the NLOS problem has been converted into LOS problem, and the LS algorithm can be applied in NLOS conditions via the cooperation of several VSs.

2.2 Geometric TOA model for NLOS and establishment of VSs

Based on the above analysis, in order to estimate the distance between BS and MS, we measure the first detective path length l ¿, which can be expressed as

h = + e (4)

where represents the first detectable path length and i is the BS number, i = 1, 2,.., M.

Consider a MS with unknown position at X = [x, y]T receiving signal from a BS at known positions, namely, Xi = [xi, yi]T, where Tdenotes the transpose operator and i is the BS number. The TOA measurement between the BS and the unknown MS is l i .The VS is set up at

xf = [xf,yf'J^, where ki is the VS number of the ith

real BS, and the estimated virtual path length is l vi, where £ is the path type, and £ = 1, 2, or 3 represents the direct path, reflection path, and diffraction path, respectively.

In indoor environments with obstacles including people, there are mainly four kinds of signal paths: direct paths, reflection paths, diffraction paths, and penetration paths.

2.2.1 Direct paths

Assuming that the signal propagates from the BS at Xi to the MS at X through a direct path, lFPi in (4) can be expressed as

lFpt = ||X'-X||

where ||-||v represents the 2-norm.

For the consistency of the whole model, we also set VS Xi1 at the real BS position Xi, and the estimated virtual

path length lvi is equal to the measured path length li

L = h.

2.2.2 Penetration paths

When signal propagates through obstacles, it slows down and gets dispersed [15]. But these effects can be neglected in comparison with other NLOS errors. The only effect of penetration we consider in our work is attenuating the signal strength which may affect the first detective path.

Indoor environments are often populated with people. The experiment results in [16] revealed that no additional ranging error occurs due to the presence of the human body in TOA measurements with the proper choice of the detection threshold. The results also showed that the power of the direct path is strongly attenuated by the presence of the human body. Therefore, we treated human body effects equivalently as penetration effects.

2.2.3 Reflection paths

In our work, we assume that all reflections are specular reflections [8, 9]. Assuming that the signal propagates from the BS at Xi to the MS at X through a reflection path, the path length we measured with the TOA method is li. As depicted in Fig. 1a, we represent the reflector AjBj as a line

y = mjX + c

where j = 1, 2,..., N represents the number of obstacle.

The VS at Xk' and the real station X. are a pair of symmetry points of line AjBj, so the position of the VS can be expressed as

yJmj-2cj

2cj X'

lFPi = ||Xki-X||2 (9)

The virtual path length can be estimated as

lvi = l (10)

If the reflection happens, the estimated length of the first detectable path should be larger than the distance between the VS and the point P, which is the intersection point of line AjBj and line xxki.

IX -Xk' II2 >

It looks like that the signal transmits from x^ to X through a direct path. Therefore, lFPi in (4) can be given as

2.2.4 Diffraction paths

Assuming that the signal propagates from the BS at X to the MS at X through a diffraction path, the path length we measured is ¡¡. As depicted in Fig. 1b,

mj = ^

mj = ^

the diffraction happens at point A; and its position is considered as a known vector. So, in (4) can be given as

lFPi —

IX*'—XIL +

ixk'-X'i

3fli£*i, «2 e*2, • «m£*m

X a1,a2,-.-,aM — X*'

- - ' = 0

Hx*' —Xi

in (12) from the floor plan

We can get information.

We establish the VS at point A,-, so that it looks like the signal transmits from the VS to the MS at Xt through a direct path and the virtual path length can be estimated as

-„i — l'— ||x*'—X'

If the diffraction happens,

IIX—X*; II2 >

|X*; —X;

In this work, we assume that the NLOS environment is simple that the signals undergo no more than one-bound scattering (reflection or diffraction), and multi-bound scattering can be ignored as they suffer from severe fading and may be too weak to be observed [8].

In NLOS indoor environment, the floor plan can be easily obtained. Then, we can develop a group of VS for every BS based on the above analysis.

2.3 Path estimation algorithm

After considering all possible signal paths and establishing VSs for them, we choose one VS from every VS

group xi1 ; x22

; we can get the corresponding

-(1 <2

virtual path length lvi, lv2.

x11 ^ J22

,*1 J*M

A1 —yM

But if the errors do exist, then we try to minimize the left part of (16),

[flx, «2,— , am] = argmin {A-^-kM }

m f - -nr

D*1,*2,^*M — E (j|Xi*1'*2.-.*M —X;*' H2— ^ 2

Then, X'

«1 ;«2;" -;«M .

is the estimated MS position.

., ivM. Then, we can utilize the existing TOA localization algorithm, least squares

for example, to obtain a possible position X*1'*1' "'*M of the MS by (15), and so on; we can get all the possible positions of the MS, the total of which is K1x K2 x ... x KM.

*1'*2'-'*m — (ATA)—1ATb

3 Analytical performance bound

The Cramer-Rao lower bound (CRLB) is defined as the lower bound on the variance of the desired parameters of the estimator (in our case, this is the localization parameter x). The bound states that the variance of any unbiased estimator is at least as high as the inverse of the Fisher information. An unbiased estimator that achieves this bound is said to be an efficient estimator that provides the lowest possible minimum mean square error.

In indoor environment, there may be several paths between a BS and a MS because of the multipath problem. To get the CRLB of the proposed localization algorithm, we use a new parameter pki to represent the probability that the path number of the measured TOA from the ith BS is ki, ki = 1, 2, ..., K'.

When the indoor signal channel is stable and the measured path is the same path between the BS and the MS in every TOA measurements. There is only one path number q

P*' —

*'— q

*11)2—(xp)2+(y102—(y§02

:*)2—(4m) 2+(y102^yMM)2-

If the TOA measurements have no errors, then

When the indoor signal channel is complicated with relatively large flow of people, the measured path may be not the same in every TOA measurements, so pki needs to be practically measured.

The Fisher information matrix of the proposed method can be given as

Jo = E

- lnf (A|w)

- lnf (A|w)

K iXK 2X^XKn

lnf (A

lnf( A

pj=i = p1 x p2 x - x Pn Pj=2 = p2 x p2 x - x pn

pj=K i = P'k 1 x p2 x - x Pn

Pj=K 1 XK2X-Xkn = PK1 x PK x - x pnN

A=1 = [^1, A2, -, AN] A=2 = [A1, A1,-, AN]

Aj=K 1 = [AK1 x A2 x - x AN]

Aj=K1xK2X-xKN = [AK1 x AK2 x - x anN]

Aiki = \j (x-*k')2+(y-yfj2+«k'

Then, we can get

x x2 y y2

,- ,Hj=K 1 XK2 X-XKn

dvN dvN

x-xf1 y-yK1

•Vl dV1 Ki K 2

x x2 y y2

X xn y yN

tVN dVN K 1xK 2 X - X KN

Jo-1 = E Pj (H/ QHj)-1 j=1

K1XK2 X - xKn

= tr(Jo-1 )= E j=1

PjJxxj ^ J;

JxxjJyyj Jxyj

(21) (22)

Jxxj + Jyyj = ¡2 JxxjJyyj- Jxyj = E ( "¡j^T

2 N / -ykA 2

" (x-^') (y-yf) ^^¡dk

4 Simulation results

To evaluate the performance of the proposed VS method, simulations have been conducted. As shown in Fig. 2, three BSs were located at BS1 (6, 8), BS2 (14, 16), and BS3 (24, 5), respectively. The MS was located at three test points, A (13, 16), B (6, 12), and C (16, 1). The first detectable path between every anchor and every test point were assumed as the figure shows. We also set two people, P1 and P2, so that the LOS path may be not the first detectable path in some conditions. When the MS was located at point A, two paths were LOS and the other one is one-bound scattering path; when at point B, one path was LOS and the others were one-bound scattering paths; when at point C, all three paths were one-bound scattering paths.

Figures 3 and 4 illustrates the root mean square error (RMSE) comparison and the cumulative distribution function (CDF) comparison among the SDP (semidefinite programming) algorithm [17], ML [10], and TSWLS [10] estimators with and without using our VS method respectively with ai = 1 m. One thousand independent runs are simulated to obtain the RMSE. The CRLB in Fig. 3 is obtained in the situation which the indoor signal channel is complicated, and we assume that p1 = p2 = ...pki. As Figs. 3a and

Fig. 3 ALE comparison. a Point A. b Point B. c Point C

RMSE(m)

4 CDF comparison. a Point A. b Point B. c Point C

4a illustrated, when the MS was located at point A, where most paths were LOS paths, the performances of the five algorithms were almost the same. When the MS was located at point B where the NLOS problem is more serious than that at A, the algorithms which used our VS method had a better performance than others, as Figs. 3b and 4b show. When the MS was located at point C, where the LOS path

did not exist, the two algorithms which used our VS method are much more accurate than those that did not use our VS method, and the advantage of the VS method is more obvious.

5 Conclusions

In this paper, we have proposed a VS method for converting the NLOS problem into LOS problem by using virtual stations, so that the conventional localization algorithm

for LOS conditions can be applied to simple NLOS conditions where the multiple-bound scattering signal paths can be ignored. We utilized only TOA measurements and the floor plan of the indoor environment and established TOA model considering reflection, diffraction, penetration, and human body effect in NLOS environments, and the performance bound of the VS method was also analyzed. Simulation results confirmed the robustness of our VS method in terms of RMSE and ALE.

Acknowledgements

The experiments of this research were financially supported the by National Natural Science Foundation of China under Grant 61601494.

Authors' contributions

DL conceived and designed the study. YW carried out most of the analyses. DL and YW drafted the manuscript. PH and YZ performed the experiments. HW modified the English expressions. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author details

1Shijiazhuang Mechanical Engineering College, 050003 Shijiazhuang, China. 2Tianjin University of Technology and Education, 300222 Tianjin, China.

Received: 23 March 2017 Accepted: 29 May 2017 Published online: 09 June 2017

References

1. Z Li, K Liu, Y Zhao, Y Ma, MaPIT: an enhanced pending interest table for NDN with mapping bloom filter. IEEE Commun Lett 18(11), 1423-1426 (2014)

2. Z Li, L Song, H Shi, Approaching the capacity of K-user MIMO interference channel with interference counteraction scheme. Ad Hoc Netw 58,286-291 (2017)

3. Z Li, Y Chen, H Shi, K Liu, NDN-GSM-R: a novel high-speed railway communication system via named data networking. EURASIP J Wirel Commun Netw 48, 1-5 (2016)

4. X Liu, Z Li, P Yang, Y Dong, Information-centric mobile ad hoc networks and content routing: a survey. Ad Hoc Netw 58, 255-268 (2017)

5. K Yu, E Dutkiewicz, Improved kalman filtering algorithms for mobile tracking in NLOS scenarios, in Proc. IEEE WCNC, 2012, pp. 2390-2394

6. K Yu, Y Guo, Non-line-of-sight detection based on TOA and signal strength, in Proc. 19th IEEE Int. Symp. PIMRC, 2008, pp. 1-5

7. S Al-Jazzar, J Caffery, H You, Scattering-model-based methods for TOA location in NLOS environments. IEEE Trans Veh Technol 56(2), 583-593 (2007)

8. CK Seow, SYTan, Non-line-of-sight localization in multipath environments. IEEE Trans Mob Comput 7(5), 647-660 (2008)

9. P Setlur, GE Smith, F Ahmad et al., Target localization with a single sensor via multipath exploitation. IEEE Trans Aerosp Electron Syst 48(3), 1996-2014 (2012)

10. KWK Lui, HC So, Range-based source localization with pure reflector in presence of multipath propagation. IEEE Electron Lett 46(13), 957-958 (2010)

11. P Meissner, K Witrisal, Multipath-assisted single-anchor indoor localization in an office environment, in Proc. 19th IEEE IWSSIP, 2012, pp. 22-25

12. D Liu, K Liu, Y Ma et al., Joint TOA and DOA localization in indoor environment using virtual stations. IEEE Commun Lett 18(8), 1423-1426 (2014)

13. B Alavi, K Pahlavan, Modeling of the TOA based distance measurement error using UWB indoor radio measurements. IEEE Commun Lett 10(4), 275-277 (2006)

14. NA Alsindi, B Alavi, K Pahlavan, Measurement and modeling of ultrawideband TOA-based ranging in indoor multipath environments. IEEE Trans Veh Technol 58(3), 1046-1058 (2009)

15. N Iya, A Muqaibel, U Johar, MA Landolsi, Ultra-wideband characterization of obstructed propagation, in Proc 7th IWCMC, 2011, pp. 624-629

16. Y Kilic, AJ Ali, A Meijerink, W Scanlon et al., The effect of human-body shadowing on indoor UWB TOA-based ranging systems, in Proc 9th Workshop on Positioning Navigation and Communication (WPNC), 2012, pp. 126-130

17. CS Tai, SY Tan, CK Seow, Robust non-line-of-sight localization system in indoor environment. Electron Lett 46((8), 593-595 (2010)

Submit your manuscript to a SpringerOpen0 journal and benefit from:

► Convenient online submission

► Rigorous peer review

► Open access: articles freely available online

► High visibility within the field

► Retaining the copyright to your article

Submit your next manuscript at ► springeropen.com