Next Article in Journal
A New Approach of Oil Spill Detection Using Time-Resolved LIF Combined with Parallel Factors Analysis for Laser Remote Sensing
Next Article in Special Issue
On Reliable and Efficient Data Gathering Based Routing in Underwater Wireless Sensor Networks
Previous Article in Journal
A Wideband Circularly Polarized Pixelated Dielectric Resonator Antenna
Previous Article in Special Issue
A Wireless Sensor Network with Soft Computing Localization Techniques for Track Cycling Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Practical Performance Analysis for Multiple Information Fusion Based Scalable Localization System Using Wireless Sensor Networks

1
Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences, Shenzhen 518055, China
2
Shenzhen Institute of Radio Testing & Tech., Shenzhen 518000, China
3
Department of Mathematics and Computer Science, Freie Universität Berlin, Berlin 14195, Germany
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(9), 1346; https://doi.org/10.3390/s16091346
Submission received: 27 June 2016 / Revised: 1 August 2016 / Accepted: 9 August 2016 / Published: 23 August 2016
(This article belongs to the Special Issue Scalable Localization in Wireless Sensor Networks)

Abstract

:
In practical localization system design, researchers need to consider several aspects to make the positioning efficiently and effectively, e.g., the available auxiliary information, sensing devices, equipment deployment and the environment. Then, these practical concerns turn out to be the technical problems, e.g., the sequential position state propagation, the target-anchor geometry effect, the Non-line-of-sight (NLOS) identification and the related prior information. It is necessary to construct an efficient framework that can exploit multiple available information and guide the system design. In this paper, we propose a scalable method to analyze system performance based on the Cramér–Rao lower bound (CRLB), which can fuse all of the information adaptively. Firstly, we use an abstract function to represent all of the wireless localization system model. Then, the unknown vector of the CRLB consists of two parts: the first part is the estimated vector, and the second part is the auxiliary vector, which helps improve the estimation accuracy. Accordingly, the Fisher information matrix is divided into two parts: the state matrix and the auxiliary matrix. Unlike the theoretical analysis, our CRLB can be a practical fundamental limit to denote the system that fuses multiple information in the complicated environment, e.g., recursive Bayesian estimation based on the hidden Markov model, the map matching method and the NLOS identification and mitigation methods. Thus, the theoretical results are approaching the real case more. In addition, our method is more adaptable than other CRLBs when considering more unknown important factors. We use the proposed method to analyze the wireless sensor network-based indoor localization system. The influence of the hybrid LOS/NLOS channels, the building layout information and the relative height differences between the target and anchors are analyzed. It is demonstrated that our method exploits all of the available information for the indoor localization systems and serves as an indicator for practical system evaluation.

Graphical Abstract

1. Introduction

Locating a target using a wireless sensor network (WSN) is an efficient way to support multiple Internet of Things (IoT) applications, and many measurement and sensing techniques are proposed [1]. The techniques of measurement, such as angle-of-arrival (AOA), time-of-arrival (TOA), time-difference-of-arrival (TDOA) and received-signal-strength (RSS), play important roles in many applications, e.g., navigation, localization, target tracking and location-based service for mobile communication [2,3]. Just like other wireless applications, localization systems are sensitive to the signal noise. The major sources that influence the measurements result from the multi-path effect and non-line-of-sight (NLOS) during the propagation of the wireless signal due to the complicated infrastructures or environment. In addition, motivated by cloud computing and urban sensing technologies, multiple sensors are deployed widely, and much information in the complicated environment is also useful to calibrate the target position, e.g., map information and the trajectory records.
The main purpose of the localization system is the position estimation accuracy for a better location-based service. Thus, many research efforts focus on developing accurate location estimation algorithms [4]. The key idea is to fuse multiple pieces of information, derive the relationship to the target position and eliminate the error. The most popular tools are the nonlinear filters based on the recursive Bayesian estimation model, e.g., extended Kalman filters, particle filters and Gaussian filters, which fuse the prior information of the previous state and the current measurement data to derive positions [5,6,7]. One important factor that may influence the target position estimation is the NLOS signals. Thus, many algorithms try to identify the NLOS measurements and mitigate or adapt them based on the prior knowledge to improve the position estimation [8,9]. A hybrid method that combines the recursive Bayesian model and NLOS information is proposed according to the jump Markov model estimation [10]. Multiple information fusion methods are also studied, which fuse both TOA and RSS measurements to improve the estimation [11,12]. The map or building layout information is also modeled for the wireless localization system and is applied for estimation [13,14]. It can be found that with the development of the sensing technique, more and more information will be integrated into the sensors and used for localization in the future. Thus, the localization system requires a scalable architecture to fuse such information. Although the indoor localization system is designed via employing multiple pieces of information, a theoretical analysis for such heterogeneous information is not provided yet.
Before the practical implementation of the WSN localization system, theoretical analysis is also required to evaluate the performance of the heterogeneous information fusion. The Cramér–Rao lower bound (CRLB) as the optimal performance indicator for the unbiased estimator is widely applied in localization and positioning systems. Theoretical investigations have been researched as the nonlinear problem for wireless localization systems. Tichavsky et al. provided the formulation of recursive posterior CRLB for nonlinear filters based on the Bayesian framework [15]. Zuo et al. proposed a conditional CRLB, which considered that the posterior probability is conditioned on the prior probability [16]. For range-based wireless localization system, many research works have provided CRLB results for different scenarios. A generalized CRLB (G-CRLB) of the wireless system is proposed for the NLOS environment [17]. The hybrid LOS/NLOS environment is analyzed, and the authors indicated that with prior knowledge of the wireless transmission channel, the estimation performance can be improved [17]. Shen et al. defined an equivalent CRLB (E-CRLB) for a general framework of the wideband wireless network [18]. The multi-path and NLOS effect are both considered, and the CRLB with or without prior information is compared to the E-CRLB [18]. A linear CRLB (L-CRLB) is proposed, which considers the linearized effect and provided the lower bound for such an estimator [19]. Other similar works also give CRLB for different ranging techniques [20,21]. Although some other methods can be used for performance analysis [22], CRLB is still popular for wireless localization researchers due to its simplicity and general expression. Although the above-mentioned CRLBs try to provide the general fundamental limits of the localization systems, they assume only parts of the factors in the localization problem, which limit their applications. In this case, these CRLBs still cannot be analyzed precisely since more information makes the environment complicated, and they are influenced by many unknown factors. Therefore, a general scalable framework of CRLB is required to further collect more information and evaluate the practical environment.
In this paper, we propose a general analysis method to evaluate the practical WSN localization system, which can fuse multiple pieces of information. The first major contribution is that we construct a scalable framework to model the multiple information fusion in the localization problem generally. We derive the Fisher information matrix (FIM) based on the proposed abstract function of all of the wireless localization system model instead of using a specific wireless propagation model. Furthermore, we employ an extendable estimated random vector θ , which contains all of the unknown parameters, e.g., the previous estimated state, the current state and the unknown indeterministic parameters that may influence the estimation accuracy. Then, we divide θ into two parts: the estimated state vector, which indicates the final position estimation, and the auxiliary vector, which helps improve the estimation accuracy. When θ only contains the state vectors or together with NLOS indicators, the our method is equivalent to the above-mentioned CRLBs. However, when a large number of parameters is involved in the estimation, e.g., the map information, the previous state, the relative height difference or the prior information of the unknown random parameters, our method can still provide practical optimal performance of the localization systems. The major advantage is that it is suitable for complicated and dynamic environment and fully considers the prior information, hybrid unknown factors and the recursive feature of the tracking algorithms.
The second contribution is that we employ the proposed CRLB framework to analyze the TOA range-based indoor localization system as a case study. The simulation environment considers all of the possible factors, e.g., the target-anchor geometry effect, the building layout, the relative height differences between the target and anchors, the NLOS transmission channel, the related prior information and the recursive feature of the tracking algorithm. The impact of each factor for the estimation accuracy is illustrated in the simulation. We use the spatial position error distribution (SPED) as a metric to evaluate the performance. The SPED indicates not only the target-anchor geometry effect, but also illustrates the impacts of height differences between the target and anchors, NLOS transmission and building layout information for the indoor localization. In addition, the impacts of the related prior information are also evaluated via SPED. To indicate which factor is important to the location estimation, we numerically evaluate the NLOS ranging, the height difference and the prior information in multiple scenarios. For the dynamic moving trajectories, the recursive form of the CRLB is applied. Finally, the estimation improvement using multiple anchors is also presented in the simulation. The results indicate that the NLOS ranging measurement mainly influences the estimation accuracy, and the prior information of the NLOS channel and the target position play important roles for improving the estimation. The relative height differences also degrade the estimation accuracy if we ignore them. However, with reliable prior information, we can make the estimation accuracy approach the location performance without any relative height differences. In general, our proposed method is suitable to exploit all of the available information to analyze the performance of the WSN localization system.
The rest parts of this paper is organized as follows: Section 2 provides the WSN localization system model; Section 3 introduces the scalable CRLB framework; then, we use the CRLB to derive the formulation of a practical TOA system in Section 4; the simulation evaluations and analysis are illustrated in Section 5; and Section 6 concludes the whole paper. To make the content more clear, we list all of the fundamental notations in Table 1 for the mathematical formulation.

2. System Model

In the WSN localization system, the mobile device with an unknown position is called the target, such as a mobile sensor node or robot. The position state of the target is denoted by x t = [ p t X p t Y ] T , where p t X and p t Y are the coordinates in the two-dimensional positioning system, and T is the transpose operator. The wireless sensor devices with known positions, which measure the ranges (or distances) to the target, are called anchors. For each anchor, the position is denoted by a j = [ a j X a j Y ] T , where a j X and a j Y are the coordinates. In this paper, we assume the WSN localization system uses the time-of-arrival (TOA) method to measure the distance between the target and anchor. Thus, we will have the relative factors that may influence the TOA and the final location estimation.

2.1. Time-Of-Arrival Ranging

In the TOA measurement method, the distance between the target and anchor is calculated according to the wireless propagation time. Consider a synchronous wireless communication, where clocks at the target and anchors are strictly synchronized. Anchors send the ranging measurements periodically via a time-division multiplexing (TDM) method based on the related WSN protocol. The target receives a radio signal transmitted from one anchor via a single propagation path. Let τ t j be the time delay of the received signal from anchor j at time t:
τ t j = 1 c | | x t a j | | + l t j
where c = 3 × 10 8 m / s is the propagation speed of the signal, a j = [ a j X a j Y ] T is the anchor position and | | · | | denotes the distance between two positions; l t j 0 is the range drift, which is caused by the NLOS effect. The indicator l t j = 0 for the LOS propagation, whereas l t j > 0 for the NLOS propagation. For many indoor systems, the TOA ranging measurement is obtained through the packet transmission time based on the network protocol; thus, the TOA is denoted as the time observation instead of the received waveform, as depicted in [17,18]:
z t j = c τ t j + v t j
where z t j is the range measurement for anchor j, the measurement function h t j ( ( d ( x t , k ) , l t j ) ) = c τ t j and v t j is the measurement noise for anchor j. The measurement noise v t j follows the zero-mean Gaussian distribution v t j N ( 0 , R j ) , where R j is the variance of the ranging noise.

2.2. TOA Noise

The localization performance mainly relies on the measurement noise. Thus, it is necessary to list the factors that are related to the TOA noise. TOA relies on the quality of detecting the direct path (DP) signal. Thus, the TOA ranging not only tolerates the attenuation of the signal strength, but also depends on the peak of all received multipath signals related to the bandwidth. Thus, the radio interference, the multipath effect and the scenario variability affect the ranging error. In addition, the TOA ranging noise also depends on the system noise. The system noise comes from the unsynchronized signals and the disordered received signals, which are transmitted from different sensors simultaneously. Apart from the NLOS effect, the comprehensive impacts on the ranging noise can be generally modeled as the normal distribution [23].

2.3. Relative Height

Another important factor that should be considered in the real-world system is that the localization happens in the 3D world. In many theoretical research works, it is assumed that the anchors and targets are on the same plane in many real localization applications. The goal is to calculate 2D positions, the X Y coordinates of the target, which is already implemented in our daily apps. In this case, the height difference between anchor and target is ignored. However, the existing height difference in the real world actually affects the accuracy in the 2D position estimation. Here, we define the height difference between the anchor and target as relative height, which is a positive variable in the Z axis. Then, we use k t 0 to denote the coordinate of the target in the Z axis and a j Z to denote the coordinate of the anchor. Thus, the relative height is k t a j z at time t. If the relative height is zero or assumed to be zero in the simulation, we define the range measurement as 2D-ranging. The distance of 2D-ranging for each anchor is formulated as:
d j ( x t , k ) = ( p t X a j X ) 2 + ( p t Y a j Y ) 2
where d j ( x t , k t ) denotes the distance function from anchor j to the target. If the relative height between the anchor and target is not zero, which is always applicable in the real case, the measurement depends on 3D coordinates, then we define the range measurement as 3D-ranging. The 3D-ranging for each anchor is formulated as:
d j ( x t , k t ) = ( p t X a j X ) 2 + ( p t Y a j Y ) 2 + ( k t a j Z ) 2
Figure 1 illustrates the difference between 2D-ranging and 3D-ranging. Suppose a person or a robot is carrying a mobile device and walking in the building. Anchor 1 and Anchor 2 are deployed on the roof, and the distance to the target depends on the target position on the ground and the relative height. Even if the target is just below the anchor, the measurement is still not zero. Anchor 3 is deployed on the same plane of the target. In this case, the range measurement is 2D-ranging. However, 2D-ranging is an ideal case, and the anchors cannot be always on the same plane of the targets. No matter whether 2D-ranging or 3D-ranging, the position estimation is still for 2D in the playing field.

3. Scalable Framework of FIM

We attempt to construct the scalable analytical framework according to the Bayesian estimation process. The Bayesian estimation is an efficient way to fuse multiple pieces of information. Firstly, we build the predict-update model of the Bayesian estimation framework. Then, we derive the scalable Fisher information matrix.

3.1. Bayesian Model

According to the Bayesian estimation framework, the relationship between the estimated state x t and the measurement z t follows:
x t = f t ( x t 1 ) + q t
where (5) is the prediction function. In (5), the target’s movement is based on the transition function f t ( ) , and q t is the prediction noise, which follows the normal distribution N ( 0 , Q t ) , where Q t indicates the covariance matrix.
In order to generally denote the wireless localization model, we propose an abstract measurement function of the unknown vectors:
z t = h t ( d ( x t , k ) , l ) + v t
where (6) is the abstract measurement function, which is a general expression based on several range measurement techniques. In (6), z t = [ z t 1 z t j z t N ] T is the measurement vector, and N denotes the number of anchors. Note that z t can represent the RSS vector or TOA vector, and (6) can be rewritten according to the different measurement techniques. Then, h t ( ) = [ h t 1 ( ) h t j ( ) h t N ( ) ] T is the nonlinear observation function, which relates to the actual received signal waveforms at the target from the anchors; v t = [ v t 1 v t j v t N ] T is the ranging noise, which is assumed as independent noise; d ( ) = [ d 1 ( ) d j ( ) d N ( ) ] T represents the distance vector between the target and anchors. The auxiliary parameters are defined as the optional unknown factors that may influence the estimation. In the position estimation, the auxiliary parameters are not necessarily calculated. However, the information of such parameters can improve the estimation accuracy. We define two kinds of auxiliary parameters that may influence z t in (6). The first one is the nonlinear auxiliary vector k = [ k 1 k j k N k ] T , which affects the actual distance together with the position state x t , and N k denotes the number of elements of k . The second one is the linear auxiliary vector l = [ l 1 l j l N l ] T , which affects the observed measurement z t and is independent of x t , where N l N indicates the number of the parameters. Since (6) is a general expression for several range measurement techniques, it should be specified and rewritten in a particular system. In this paper, (6) is simplified as the TOA formulation, which is discussed in the following sections. According to the Bayesian theorem, the posterior probability of x t is expressed as p ( x t | z t , x t 1 ) = p ( x t | x t 1 ) p ( z t | x t ) , where t 1 indicates the previous time interval, and p ( x t | x t 1 ) is the prior probability [24].

3.2. Recursive FIM

Our analysis fully considers all of the possible unknown random factors that may influence the position estimation; hence, the parameter vector includes: the current state x t , the previous state x t 1 and auxiliary parameter vectors k and l . Thus, θ is expressed as:
θ x t T x t 1 T k T l T T
The CRLB, which is given by the inverse of the Fisher information matrix (FIM), sets the lower limit for the variance (or covariance matrix) of any unbiased estimators of an unknown parameter (or unknown parameters) [25]. If p ( θ , z t ) denotes the joint probability density function (PDF) of observations z t and the state θ , then the score function is defined as the gradient of its log-likelihood:
U ( x t ) θ ln p ( θ , z t )
where θ = ` 1 , , ` N T is the operator of first order partial derivatives. The FIM, J ( θ ) , is the covariance of the score function:
J ( θ ) E θ ln p ( θ , z t ) θ ln p ( θ , z t ) T
where E · indicates the expectation operator. Additionally, the CRLB is just the inverse of FIM, and the estimation covariance cannot be lower than it:
Cov θ ( θ ˜ ) { J ( θ ) } 1
where “ A B ” should be interpreted as matrix A B is non-negative definite.
Since p ( θ , z t ) = p ( z t | θ ) p ( θ ) is based on the Bayesian theorem, it is easily seen that J ( θ ) can be decomposed into two parts:
J ( θ ) = J D ( θ ) + J P ( θ )
where J D ( θ ) represents the information obtained from measurement data and J P ( θ ) represents the prior information.
Firstly, we use the notations h = h t ( d ( x t , k ) , l ) , h j = h t j ( d ( x t , k ) , l ) and decompose J D using the chain rule as:
J D ( θ ) = H · J h · H T
where H = [ θ h ] and J h is the FIM conditioned on h :
J h = E h ln p ( z t | θ ) h ln p ( z t | θ ) T
The matrix H is further decomposed into four components:
H = H t H t 1 K L T
where H t = [ x t h ] 2 × N , H t 1 = [ x t 1 h ] 2 × N , K = [ k h ] N k × N and L = [ l h ] N l × N . Since d is independent of the previous state x t 1 , H t 1 = 0 . Then, H is formulated as:
H = H t 0 K L T
For J h , we can use diagonal matrices of order N to represent it:
J h = Λ = diag ( λ 1 , , λ j , , λ N )
where the diagonal term λ j depends on h t j ( ) and will be derived based on the typical system later. Then, J D , which is a ( 4 + N k + N l ) × ( 4 + N k + N l ) matrix, is written as:
J D = H · Λ · H T = D 11 0 D 13 D 14 0 0 0 0 D 13 T 0 D 33 D 34 D 14 T 0 D 34 T D 44
where:
D 11 = H t Λ H t T D 13 = H t Λ K T D 14 = H t Λ L T D 33 = K Λ K T D 34 = K Λ L T D 44 = L Λ L T
Without prior information, J D cannot indicate the CRLB directly, because det ( J D ) = 0 and J D is not reversible. However, the CRLB can also be attained through a calculation rule. The calculation rule is illustrated in the next subsection. Next, we derive J P . The prior probability for θ is extended as p ( θ ) = p ( x t | x t 1 ) p ( k ) p ( l ) , then the prior information is written as:
ln p ( θ ) = ln p ( x t | x t 1 ) + ln p ( k ) + ln p ( l )
where p ( k ) and p ( l ) are the independent prior information of x t and x t 1 . If we decompose θ into two sub-vectors: the state vector [ x t x t 1 ] T and the auxiliary vector [ k l ] T . Then, J P can be formulated as:
J P = E θ ln p ( θ ) θ ln p ( θ ) T = J P 11 J P 12 J P 12 T J P 22
where J P 11 is the recursive form of x t and x t 1 , which is formulated by Tichavsky et al. [15]:
J P 11 = M 11 M 12 M 12 T M 22 + J ( x t 1 )
where:
M 11 = Q t 1 M 12 = x t 1 f t ( x t 1 ) Q t 1 M 22 = x t 1 f t ( x t 1 ) Q t 1 x t 1 f t ( x t 1 ) T
where J ( x t 1 ) is the previous FIM of x t 1 . Additionally, J P 12 are 0 matrices, since p ( k ) and p ( l ) are independent of x t and x t 1 :
J P 12 = 0 0 0 0
The prior distribution p ( x t | x t 1 ) is also independent of l and k ; thus, J P 21 = J P 12 T = 0 . Finally, the last element J P 22 is expressed as:
J P 22 = J K 0 0 J L
where J K and J L are the FIMs conditioned on k and l , respectively:
J K = E k ln p ( k ) k ln p ( k ) T J L = E l ln p ( l ) l ln p ( l ) T
The formulations of J K and J L are based on the features of k and l , which are discussed in the later sections. Then, substitute (17) and (20) into (11); we obtain:
J ( θ ) = M 11 + D 11 M 12 D 13 D 14 M 12 T M 22 + J ( x t 1 ) 0 0 D 13 T 0 D 33 + J K D 34 D 14 T 0 D 34 T D 44 + J L
However, J ( θ ) is a matrix of high dimensions, while only a small submatrix J 1 ( x t ) [ J 1 ( θ ) ] 2 × 2 is of interest. Using the form of the Schur complement of the submatrix [26], we decompose the FIM into two parts: the state matrix J S and the auxiliary matrix J A :
J ( x t ) = J S J A
where:
J S = M 11 + D 11 M 12 M 22 + J ( x t 1 ) 1 M 12 T J A = D 13 D 14 D 33 + J K D 34 D 34 T D 44 + J L 1 D 13 D 14 T
Additionally, the formulation of each element can be found in (18), (22) and (25).

3.3. Calculation Rule

Equation (27) only holds when all of the elements in θ are to be estimated and the prior information for the whole θ is available. If the CRLB is used to analyze a specific localization system, not all of the elements are necessary for θ , and some vectors are absent sometimes. For instance, for the non-recursive scenario, the system does not consider x t 1 . In addition, when the system has a deterministic value of the auxiliary vectors, k and l are not estimated and useless for J ( x t ) . Thus, the calculation principle for our CRLB is that: when any vector in θ is absent, the related matrix in (27) turns into 0 , and we will treat such a 0 matrix as the empty matrix, then we mitigate the empty matrix for further calculation.
Here, we consider two examples. The first example is that the estimated vector θ does not include x t 1 and k , then M 12 , M 22 , J ( x t 1 ) , J K , D 13 , D 33 and D 34 are 0 and mitigated from (27):
J ( x t ) = M 11 + D 11 D 14 D 44 + J L 1 D 14 T
which has the same form as the formulation of E-CRLB in the wideband NLOS environment [18] and L-CRLB for the linear formulation [19].
The second example is that the estimated vector θ does not include k and l , and J A is empty, then only J S remains, which is the form for the recursive nonlinear filter [15]:
J ( x t ) = M 11 + D 11 M 12 M 22 + J ( x t 1 ) 1 M 12 T

4. Application to a TOA WSN Localization System

In this section, we apply the proposed CRLB to analyze the time-of-arrival (TOA)-based WSN localization systems. Note that this paper provides a scalable framework for theoretical analysis rather than a particular analysis. Thus, our goal is not to find the major impacts or factors via mathematical derivation, but through the numerical simulations. In addition, according to the formulation of (27), many factors are involved in the calculation. It is not easy to analyze them via a single formulation. Thus, after enumerating the important factors in the TOA system, we will evaluate the performance via numerical analytical simulations. The impacts of the practical conditions, e.g., the real 3D anchor deployment and the NLOS transmission based on the building layout, are formulated. The related parts in (26) are derived for both cases with and without prior knowledge.

4.1. Relative Height

As mentioned before, no matter whether 2D-ranging or 3D-ranging, the position estimation is still for 2D in the playing field, which means we just want to obtain x t = [ p t X p t Y ] T and not k t . However, for CRLB analysis, we still consider if k t is involved in the calculation and check whether the estimation accuracy is influenced. Then, H t is obtained:
H t = d t 1 ( x t , k t ) p t X d t j ( x t , k t ) p t X d t N ( x t , k t ) p t X d t 1 ( x t , k t ) p t Y d t j ( x t , k t ) p t Y d t N ( x t , k t ) p t Y
where:
d t j ( x t , k t ) p t X = p t X a j X ( p t X a j X ) 2 + ( p t Y a j Y ) 2 + ( k t a j Z ) 2 d t j ( x t , k t ) p t Y = p t Y a j Y ( p t X a j X ) 2 + ( p t Y a j Y ) 2 + ( k t a j Z ) 2
Then, K = [ k t d ( x t , k t ) ] 1 × N is expressed as:
K = d t 1 ( x t , k t ) k t d t j ( x t , k t ) k t d t N ( x t , k t ) k t
where:
d t j ( x t , k t ) k t = k t a j Z ( p t X a j X ) 2 + ( p t Y a j Y ) 2 + ( k t a j Z ) 2
Since x t and k t are within the same function d ( x t , k t ) , H t and K are both derived from d ( x t , k ) . In order to analyze the impacts of x t and k t , we decompose the matrix and obtain (31) and (33). For a 3D positioning case, we can combine H t and K together to indicate the 3D position state. However, it is more convenient to separate H t and K to analyze the estimation performance of x t in a 2D positioning case.
If any prior information is unknown to the system, substitute (31) and (33) into (27), then we obtain:
J ( x t ) = H t Λ H t T H t Λ K T ( K Λ K T ) K Λ H t T
For the prior information of k t , we assume that the target is always above the ground, which is k t 0 . Thus, we apply the Gamma distribution to indicate the potential distribution of k t , where k t G ( α k , β k ) ( k t ) = ( β k ) α k Γ ( α k ) k t α k 1 exp ( β k k t ) , α k is the shape parameter and β k is the rate parameter. For the Gamma distribution, J k is complicated. To obtain an analytical expression, we assume α k > 2 for simplicity. Then, the Gamma function is Γ ( α k ) = 0 + exp ( x ) x α k 1 d x . Thus, J k = E k G ( α k , β k ) ( k ) 2 is derived as:
J k = β k 2 2 ( α k 1 ) β k Γ ( α k ) 0 + β k α k k t α k 2 exp ( β k k t ) d k t + ( α k 1 ) 2 Γ ( α k ) 0 + β k α k k t α k 3 exp ( β k k t ) d k t
Use the property Γ ( α k ) = α k Γ ( α k 1 ) , and substitute it into (36); we obtain:
J k = β k 2 α k 2

4.2. NLOS Impact

We use the vector l to represent the NLOS indicator. In the TOA system, the NLOS delays the wireless packet propagation linearly. Thus, l is the linear auxiliary vector. We assume there are N l N NLOS measurements, and the drift for each measurement is independent of the others, then L = [ l h ] N l × N is formulated as:
L = h t 1 l 1 h t N l l 1 h t N l 1 h t 1 l N l h t N l l N l h t N l N l = I N l 0
where I N l is the identity matrix of order N l , and the rest is a N l × ( N N l ) zero matrix due to the independent condition for the LOS measurement. As mentioned before, the zero matrix is mitigated during the calculation, then L = I N l . Note that L indicates the sub-matrix of the FIM, which means that the NLOS condition is unknown to the system. Then, the system should also estimate the NLOS parameters together with the position state in the real environment. In this case, the estimated error is propagated and degrades the position estimation.
Since the range drift for the NLOS is also nonnegative, we still use the Gamma distribution as the prior information l m G ( a m , b m ) ( l m ) = ( b m ) a m Γ ( a m ) l m a m 1 exp ( b m l m ) according to [17], where a m 2 is the shape parameter, b m is the rate parameter and m is the m-th NLOS measurement. Similar to (37), we obtain J L :
J L = diag ( b 1 2 a 1 2 , , b m 2 a m 2 , , b N l 2 a N l 2 )
where J L is the prior matrix for L , which illustrates that the system can use the prior information to determine whether a measurement is LOS or NLOS if this information is unknown. In this case, the estimation error can be reduced to some extent.

5. Simulation

We set up several WSN localization simulations to evaluate the analytical performance. In each simulation, we consider several different factors, e.g., the recursive process during the target tracking, estimations with and without considering l and k t . To make the results clear, we mark the CRLBs for different situations by adding superscripts and subscripts, which can be depicted as CRLB . The subscripts indicate the considered vectors, including the state vector and the auxiliary vector. This means that the system contains other factors in the ranging measurement, e.g., relative height or NLOS measurement or both. The superscripts indicate the available prior information of the related vectors. For instance, if we want to simulate the estimation with the NLOS range drift, the results of the CRLB are marked by CRLB x t , l . Additionally, if the prior information of x t is attained in the simulation, the results are marked by CRLB x t , l x t . For the recursive estimation, we use the notation CRLB x t , l x t 1 , l to indicate the results, which considers the prior state x t 1 . Note that, the superscript of x t in CRLB x t , indicates that we have arbitrary prior knowledge of x t . This information is provided directly by the system instead of recursive estimation; while the superscript of x t 1 in CRLB x t 1 , also means the prior knowledge of x t . However, such information is obtained based on the previous estimation of x t 1 in the recursive target tracking. Thus, we use CRLB x t , to analyze the static positioning scenarios and use CRLB x t 1 , to analyze the dynamic continuous target tracking scenarios. Some CRLB notations that appear in the following sections are listed below.
  • CRLB x t : Classical CRLB without any prior information, relative height or NLOS measurement
  • CRLB x t , k t : CRLB when the relative height exists
  • CRLB x t , l : CRLB with the NLOS measurement
  • CRLB x t , k t , l : CRLB both contains the relative height and NLOS measurement
  • CRLB x t , k t k t : CRLB with the prior information of the relative height
  • CRLB x t , l l : CRLB with the prior information of NLOS measurement
  • CRLB x t , x t , : CRLB with the arbitrary prior information of x t in the static scenario
  • CRLB x t , x t 1 , : CRLB with the prior information based on x t 1 in the dynamic scenario
To approach the real environment, we initially set the related parameters according to the test building in [27]. Since the localization system can be affected by many factors, we tune the parameters to provide a comprehensive analysis in the following simulations, e.g., k t , l and the number of anchors. To better understand the performances in multiple environments, we uniformly locate the anchors in the SPED evaluation. Additionally, we evaluate the performances considering the relative height, NLOS and recursive estimation based on randomly-deployed anchors and Monte Carlo simulation to draw a general performance. In the real experiment, the anchors are randomly deployed throughout the whole building.

5.1. Spatial Position Error Distribution

In the first simulation, a 100 × 100 m 2 playing field with four anchors is constructed. All of the anchors are deployed on the roof being 2.5 m high, and the target is about 0.5∼1 m. To approach the real applications, we use several statistical results according to [27]. We set the relative height as the constant value 1.5 m. The range error for each anchor follows zero-mean Gaussian distribution v t j N ( 0 , R t j ) , where R t j is the variance of v t j and is set to 5 2 . The range drift for the NLOS measurements is set to 2 m. For the prior information, k t G ( 2 . 5 , 2 ) ( k ) . The prior information of the NLOS range drift l m is l m G ( 3 . 5 , 1 . 8 ) ( l m ) . For the position state prior information, we assume that the prediction function is a linear static identity matrix with the zero-mean Gaussian prediction noise q t N ( 0 , Q t ) , where Q t = diag ( σ x 2 , σ y 2 ) is the covariance of q t . We assume σ x = σ y = 2 m.
We apply the CRLB to indicate the optimal squared error, which is tr ( J 1 ( x t ) ) . To illustrate the target-anchor geometry effect for the 2D localization system in the playing field, we employ tr ( J 1 ( x t ) ) to depict the spatial position error distribution (SPED) [28]. The SPED is defined as the distribution of the position error for every possible target position, which estimates the positions point by point in the whole playing field and draws statistical results. Thus, every position in the playing field is assumed as the target position during the simulation. The SPED is derived according to all of the statistical results of the estimation error of all of the positions in the playing field. It illustrates how the performance changes according to the target-anchor (RX-TX) geometry effect. When the anchor positions and the error model change, the SPED changes accordingly, which helps to understand the relationship between the anchor deployment and the algorithm. The SPED results are represented by the contours in Figure 2, where the four anchor positions are marked by triangles.
In Figure 2a, the SPED is drawn based on the classical CRLB ( CRLB x t ) in which the unknown vector is only x t , and no relative height or NLOS propagation is introduced. When both k t and l are introduced in the simulation, the dimension of the estimated vector is increased again. Then, a special pattern of the contours is drawn in Figure 2b, which demonstrates that the estimation error is propagated and increased if more unknown factors are introduced. Each anchor seems to be a center of an independent contour area, and the estimation error is extremely larger than Figure 2a, which indicates the strong uncertainty and a special geometric relationship between the target position and the anchor positions. Since more factors cannot be ignored in the complicated environment, the localization problem turns out to be a high dimensional estimation problem. It is quite possible that the squared error of the high dimensional estimation is larger than the low dimensional estimation, as the probability of the wrong estimation increases if more unknown parameters appear. Additionally, it also imposes new error on the original x t estimation error. However, such a high dimensional localization problem cannot happen in the real world. On the one hand, the calculation complexity is increased dramatically in the real system, and only x t is useful, which is a waste of computational resources. On the other hand, the algorithm designers ignore some unimportant factors, simplify the calculation and increase the estimation accuracy based on the prior knowledge.
Then, the prior information of x t , k t and l is used in the CRLB to improve the estimation accuracy. The SPEDs with prior information are depicted in Figure 2c,d. When the prior state information of x t is introduced in Figure 2c, the estimation error is reduced effectively from more than 5 m in Figure 2a to 1.4 m below. In addition, the position error is almost the same everywhere, although the geometric pattern is similar to Figure 2a, which indicates that the geometric impact is reduced by using the prior information of x t . In Figure 2d, we assume that all of the prior information is available, then the estimation error is reduced to below 1.4 m, just as Figure 2c. Although Figure 2d still has a geometric pattern, which is similar to Figure 2c, the errors in different positions are almost the same. Thus, the geometric effect is actually reduced. In addition, since the estimation error is not further reduced, the results in Figure 2d can be the limits based on the prior information. It is also demonstrated that the prior information of x t significantly reduces the estimation error and is the most important prior information especially in the recursive estimation. For the other two parts of the prior information, they are used for the improvement of the impacts of the relative height and the NLOS effects.

5.2. Impact of Relative Height

In this simulation, we evaluate the impact of the relative height independently. For the indoor localization, the relative height cannot be too high unless the room is tall enough. Thus, we choose the value of k t between 0.5 m and 10 m. The simulation evaluates CRLB x t , k t and also CRLB x t as a comparison. To illustrate the general relationship between the relative height and the measurement noise, we also tune the variance of the noise from 0 . 5 2 to 5 2 and k t from 0.5 m and 10 m. In this simulation, we run 1000 Monte Carlo experiments. In each experiment, we choose a random position for analysis. The averaged results are listed in Figure 3.
It is clearly observed that the estimation error rises according to the increase of k t . However, the increased value is quite small. Take Figure 3c for instance; the increased error is only 0.06 m when k t is tuned from 0.5 m to 10 m, and the estimation error is 3.43 m by that time. Therefore, the value of k t does not affect the estimation accuracy too much. On the other hand, the existence of k t does affect the performance no matter what the value of CRLB x t is. Compare the results in Figure 3c; k t affects the estimation with a typical certain value. When k t is small, the average RMSE of CRLB x t is 1.4390 m, but the average RMSE of CRLB x t , k t = 0 . 5 is 1.6856 m and CRLB x t , k t = 10 is 1.7137 m. According to the analytical results, when the relative height is introduced into the simulation, there is a gap between 2D localization without the relative height. Additionally, the gap becomes larger with the increased value of the measurement noise. Take v t j N ( 0 , 0 . 5 2 ) for instance; when k t = 0 , the average RMSE is only 0.84 m. When v t j N ( 0 , 5 2 ) , the average RMSE without the relative height is about 8.3 m, which is almost 10-times that of the low noise. Thus, the relative height really affects the estimation no matter the value.
In addition, we evaluate the impact of the prior knowledge of k t in the multiple noise environment. The simulation results are presented in Figure 4. It is clearly observed that the prior knowledge can significantly improve the RMSE. If v t j N ( 0 , 5 2 ) , the prior information of k t can even reduce 1 m RMSE. Additionally, if the error is small, the improvement is only a little bit. However, with the increased value of the relative height, the RMSE rises accordingly even with the prior information. For the real application, the relative height can be within 10 m, in which the improvement can effectively reduce the RMSE to a reasonable range. Therefore, even if the impact of the relative height is limited, it is still necessary to employ the prior information.

5.3. Impact of NLOS

The analysis of NLOS measurements has been mentioned in several literature works [17,18]. Here, we evaluate the impact of the number of NLOS measurements through simulations. The playing field is still 100 × 100 m 2 with 16 randomly-deployed anchors. We set up 1000 Monte Carlo runs, and a random set of anchors is chosen to have NLOS measurements. We tune the number of NLOS measurements from 1 to 16 and evaluate the performance of CRLB x t and CRLB x t , l l . The averaged results are depicted in Figure 5. It is observed that the RMSE of CRLB x t , l increases according to the rise of the NLOS measurement number. It is also possible to have useless estimation when all of the measurements are NLOS, just as mentioned in [17]. When there is only one NLOS measurement, the estimation error reaches 3 m. Additionally, the error is increased to 5 m if the number of NLOS measurements is changed to 16. In real applications, the number of NLOS measurements is unknown, and the system has to calculate all of the parameters from all of the measurements. Thus, the estimation error is quite high without the help of prior information. In addition, with the prior information of l , the estimation error is almost a constant value. Thus, the prior information of the NLOS measurement can successfully reduce the NLOS effect, which is applicable for real applications. Such prior information can be a combination of the building layout information and the wireless propagation model, e.g., map-matching-based algorithms [13]. Additionally, several literature works have proposed many algorithms that either use pre-assumptions or estimate the parameters to obtain the NLOS information. The NLOS-related methods are beyond the scope of this paper; please refer to other research works in [20,23,29].

5.4. Map Assist Localization System

In the following simulation, we evaluate a practical scenario where the layout of the building is involved. The playing field is still 100 × 100 m 2 . There are four big rooms located at four corners of the playing field. The area for each room is 40 × 40 m 2 . The rest of the playing field is the hallways. Four anchors are deployed inside each room and placed at the four corners. Thus, 16 anchors are uniformly deployed in the playing field, which is depicted in Figure 6a. The triangles mark the positions of anchors. The LOS measurements can only be obtained in one room with four associate anchors, and others measurements are NLOS. For the positions in the hallways, all of the measurements are NLOS rangings.
All of the anchors are deployed on the roof with relative height 1∼2 m to the target. The distributions of the range error and the related prior information are the same in Figure 2. The range drift for the NLOS measurements depends on the signal transmission from the anchor to the target. According to our previous research, the average positive NLOS drift of the signal through one wall is set to 2 m and to 5 m for the signal through two or three walls if the signal can be detected [30]. The prior information of the NLOS range drift l m depends on the signal transmission. Here, we use the conclusions based on the experimental models in [27]: the distribution for the signal transmission through one wall is l m G ( 3 . 5 , 1 . 8 ) ( l m ) ; the distribution for the signal transmission through more than two walls is l m G ( 4 . 5 , 2 . 2 ) ( l m ) . In this simulation, we evaluate the SPED of CRLB x t , k t , l , CRLB x t , k t , l k t , l and CRLB x t , k t , l x t , k t , l , which means that every position is calculated, and the error contours are drawn to indicate the target-anchor geometry effect. The results are depicted in Figure 6.
For numerical comparison, the RMSE in Figure 6b is higher than in Figure 6c,d, which is more than 3.39 m in the central area. Additionally, the error becomes higher and higher when the position is approaching the corner, which is more than 8 m. Due to the lack of prior information, the geometric shape does not have special characteristics, which are related to rooms or corridors. The contours are almost like rectangles located in the center of the playing field. When the prior information of k t and l is introduced, the accuracy is significantly improved, which is reduced to 2.55 m on average. The geometric shapes of the contours are different in the rooms and hallways. This indicates that the localization algorithms using the prior knowledge of NLOS conditions based on the building layout information and the NLOS identification and mitigation methods can reasonably improve the estimation performance. Thus, the layout information in the building map is an important information source for localizations. When the prior information of x t is introduced in the estimation as indicated in Figure 6d, the RMSE is further reduced, which is 1.235 m in almost all of the playing field where the target-anchor geometry effect is reduced effectively.

5.5. Bayesian-Based Target Tracking Estimation

In this simulation, we evaluate the performance of the recursive Bayesian estimation for target tracking. Since the impacts of k t and l are extensively analyzed above, this section mainly focuses on analyzing the impact of x t 1 . The playing field is still 100 × 100 m 2 , and 16 anchors are deployed randomly in the field. We evaluate the CRLB in three scenarios: The first one considers NLOS and relative height measurement; the second scenario considers only the relative height; the third scenario only has NLOS measurement. In each scenario, both CRLB without prior information and the recursive CRLB with prior information based on (27) are analyzed. Therefore, we evaluate CRLB x t , k t , l and CRLB x t , k t , l x t 1 , k t , l in the first scenario, CRLB x t , k t and CRLB x t , k t x t 1 , k t in the second scenario and CRLB x t , l and CRLB x t , l x t 1 , l in the third scenario. We assume that the prior information is unknown initially in the recursive estimation. This information can only be obtained recursively after the initial estimation is obtained. Then, the recursive estimation is applied.
We run 1000 Monte Carlo simulations, and the target moves in a separate random path in each simulation. In addition, the target can also be static. Since x t 1 can also be estimated in the static scenario and be used for recursive estimation, the analysis results are the same as the dynamic target tracking scenarios. The estimation results are averaged and represented by the RMSE in Figure 7. There are three solid parallel straight lines, which indicate the estimations without prior information: CRLB x t , k t , l , CRLB x t , k t and CRLB x t , l . The three other dashed curves illustrate the recursive estimations according to time steps, which are CRLB x t , k t , l x t 1 , k t , l , CRLB x t , k t x t 1 , k t and CRLB x t , l x t 1 , l .
The Bayesian recursive estimation method with related prior information effectively reduced the estimation as indicated in Figure 7. The RMSEs of the three curves gradually converge to low values according to time steps. The impacts of the relative heights and the NLOS drifts still degrade the estimation performance. Even with the recursive estimation, the estimation error cannot be further reduced, where the CRLB x t , k t , l x t 1 , k t , l is 0.5 m larger than CRLB x t , k t x t 1 , k t when t = 20 .

5.6. Multiple Anchors

The relationship between the CRLB and the number of anchors is illustrated in Figure 8. We randomly deploy multiple anchors in the playing field. The number of anchors is adapted from four to 30. Both k t and l are considered in this simulation. We assume that all of the measurements are NLOS. In this simulation, we evaluate CRLB x t , k t , l , CRLB x t , k t , l k t , l , CRLB x t , k t , l x t 1 and CRLB x t , k t , l x t 1 , k t , l . With a few anchors, CRLB x t , k t , l contains large error, which achieves more than 15 m. When the number of anchors is four, the prior information of k t and l can reduce much of the error as indicated by the curve of CRLB x t , k t , l k t , l . Furthermore, the recursive estimation based on CRLB x t , k t , l x t 1 can reduce much of the error, which is even smaller than CRLB x t , k t , l k t , l . The error of the recursive estimation based on CRLB x t , k t , l x t 1 , k t , l is only a little lower than CRLB x t , k t , l x t 1 . It is demonstrated that compared to the prior information of k t and l , x t 1 is the dominant factor for reducing the estimation error. However, when the number of anchors is increased, all of the curves converge gradually to a low value. In this case, the prior information does not improve much for the estimation. Thus, even if the prior information is not available, using more observations still improves the estimation accuracy effectively.

5.7. Practical Evaluation

For practical usage, we employ the CRLB to evaluate a reference system. In this system, we deployed 17 wireless sensor nodes either along the corridor or in the offices of the research building. A robot carrying a sensor node as a target moved along the corridor of the building with constant speed while recording its own positions [31]. All sensors are integrated with the nanoPAN 5375 RF (Nanotron Tech. GmbH, Berlin, Germany) module with a 2.4-GHz transceiver and a 1-Mbps data rate for range measurement, the LPC 2387 as the micro-controller (Nanotron Tech. GmbH, Berlin, Germany) and the CC1101 900-MHz transceiver (Nanotron Tech. GmbH, Berlin, Germany) as the radio transceiver for communication. The data collected from sensor nodes are also range measurement values, which are based on TOA. Figure 9 depicts the map of our experimental building. The triangles, which are randomly deployed, mark the sensor node positions.
According to the collected data and map information, we construct the distribution model for both LOS and NLOS measurements, which are a zero-mean Gaussian distribution v t j N ( 0 , R t j ) ; R t j is set to 5 2 for LOS measurement; and the NLOS measurements drift is set to 2 m. Then, the SPEDs with or without prior knowledge are depicted in Figure 10. As illustrated in Figure 10a, the overall estimation error can be extremely high due to the measurement error. In addition, the anchor-target geometry shape is seriously distorted due to the randomly-deployed anchor. In the corners of the building, the RMSEs are changing rapidly. However, with the help of prior information, we can effectively reduce the error to about 1.25 m.

5.8. Discussion

In this section, we use our proposed method to evaluate the optimal performance in multiple environments. It can be observed that our method can adapt the estimated parameters to fit the practical environment. In addition, the scalable architecture can effectively fuse multiple pieces of information, which outperforms other methods, which only consider parts of the information. In Table 2, we compare our method to three mainly used CRLBs, which are generalized CRLB [17], conditional CRLB [16] and equivalent CRLB [18]. The first row in Table 2 indicates the considered information, and the word “prior” means the related prior information. As illustrated in the table, the generalized CRLB can only analyze the NLOS effect. The conditional CRLB provides the formulation of the state vector instead of other factors. The equivalent CRLB exploits more information. However, it is a closed form formulation, which is not scalable and cannot fuse the relative height information.

6. Conclusions

In this paper, we propose a scalable analyzing method for a WSN localization system, which can fuse multiple heterogeneous information to indicate the optimal performance. Theoretically, we divide the estimated vector θ into three parts: the estimated state vector and two auxiliary vectors. The recursive formulation of FIM is provided, which fully considers all of the possible factors that may influence the estimation accuracy, and it exploits all of the available information to derive the fundamental limits. It is a suitable tool to indicate the optimal estimation bound of practical systems.
We employ our theoretical contribution to analyze the TOA range-based WSN localization system. The impacts of the height k t and NLOS range drift l are considered as the auxiliary vectors. The target-anchor geometry, prior information and the recursive form are extensively analyzed in the simulation. In addition, we employ a real test-bed and practical data to evaluate the overall performance. In the simulation, we find that both the relative height and NLOS impact can heavily degrade the estimation performance. However, many pieces of available information can improve it, e.g., the prior state distribution, the prior knowledge of the building and multiple anchors.
According to the simulation and experimental demonstration, the proposed CRLB is a general framework for analyzing the WSN localization systems, and it is not restricted to any specific technique. Future work will use the scalable CRLB to exploit other localization techniques and in other complicated environments to find potential factors that may influence the performance.

Acknowledgments

This work is partially supported by the Major Project of the Chinese National Programs for Fundamental Research and Development (973 Program, No. 2015CB352400), the National Nature Science Foundation of China (Grant No. 61501443, U1401258), the Science and Technology Planning Project of Guangdong Province (2015B010129011), the Shenzhen Overseas High-Caliber Personnel Innovation Funds (KQCX20140521115045446), the Shenzhen-Hong Kong Innovative Technology Cooperation Funding (Grant No. SGLH20131009154139588), the Shenzhen Fundamental Research Project (Grant No. JCYJ20150630151329298, JCY20151117161854942), Shenzhen Peacock Innovation Team (Grant No. KQTD2015071715073798), the Leading Talents of Guangdong Province Program (Grant No. 00201510).

Author Contributions

Yubin Zhao proposed the scalable framework of the CRLB and derived the formulation. Xiaofan Li built the simulation and processed the data. Sha Zhang processed the practical data and obtained the fitting results of the prior information, e.g., NLOS information and relative height information. Tianhui Meng constructed the real test-bed evaluation. Yiwen Zhang executed the simulation of the relative height part and polished the whole paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, H.; Darabi, H.; Banerjee, P.; Liu, J. Survey of wireless indoor positioning techniques and systems. IEEE Trans. Syst. Man Cybern. Part C Appl. Rev. 2007, 37, 1067–1080. [Google Scholar] [CrossRef]
  2. Axel, K. Location-Based Services: Fundamentals and Operation; John Wiely & Sons: Hoboken, NJ, USA, 2005. [Google Scholar]
  3. Mao, G.; Fidan, B.; Anderson, B. Wireless sensor network localization techniques. Comput. Netw. 2007, 51, 2529–2553. [Google Scholar] [CrossRef]
  4. Patwari, N.; Ash, J.N.; Kyperountas, S.; Hero, A.O., III; Moses, R.L.; Correal, N.S. Locating the Nodes: Cooperative Localization in Wireless Sensor Networks. IEEE Signal Process. Mag. 2005, 22, 54–69. [Google Scholar] [CrossRef]
  5. Chen, B.S.; Yang, C.Y.; Liao, F.K.; Liao, J.F. Mobile location estimator in a rough wireless environment using extended Kalman-Based IMM and Data Fusion. IEEE Trans. Veh. Technol. 2009, 58, 1157–1169. [Google Scholar] [CrossRef]
  6. Farahmand, S.; Roumeliotis, S.; Giannakis, G. Particle filter adaptation for distributed sensors via set membership. In Proceedings of the 2010 IEEE International Conference on Acoustics Speech and Signal Processing (ICASSP), Dallas, TX, USA, 14–19 March 2010; pp. 3374–3377.
  7. Kotecha, J.; Djuric, P. Gaussian particle filtering. IEEE Trans. Signal Process. 2003, 51, 2592–2601. [Google Scholar] [CrossRef]
  8. Mazuelas, S.; Lago, F.; Blas, J.; Bahillo, A.; Fernandez, P.; Lorenzo, R.; Abril, E. Prior NLOS measurement correction for positioning in cellular wireless networks. IEEE Trans. Veh. Technol. 2009, 58, 2585–2591. [Google Scholar] [CrossRef]
  9. Prieto, J.; Bahillo, A.; Mazuelas, S.; Lorenzo, R.; Blas, J.; Fernandez, P. NLOS mitigation based on range estimation error characterization in an RTT-based IEEE 802.11 indoor location system. In Proceedings of the IEEE International Symposium on Intelligent Signal Processing, WISP 2009, Budapest, Hungary, 26–28 August 2009; pp. 61–66.
  10. Nicoli, M.; Morelli, C.; Rampa, V. A jump markov particle filter for localization of moving terminals in multipath indoor scenarios. IEEE Trans. Signal Process. 2008, 56, 3801–3809. [Google Scholar] [CrossRef]
  11. Song, Y.; Yu, H. A new hybrid TOA/RSS location tracking algorithm for wireless sensor network. In Prceedings of the 9th International Conference on Signal Processing, ICSP 2008, Beijing, China, 26–29 October 2008; pp. 2645–2648.
  12. Prieto, J.; Mazuelas, S.; Bahillo, A.; Fernandez, P.; Lorenzo, R.M.; Abril, E.J. Adaptive data fusion for wireless localization in harsh environments. IEEE Trans. Signal Process. 2012, 60, 1585–1596. [Google Scholar] [CrossRef]
  13. Montorsi, F.; Pancaldi, F.; Vitetta, G. Map-aware models for indoor wireless localization systems: An experimental study. IEEE Trans. Wirel. Commun. 2014, 13, 2850–2862. [Google Scholar]
  14. Haque, I.; Assi, C. Profiling-based indoor localization schemes. IEEE Syst. J. 2015, 9, 76–85. [Google Scholar] [CrossRef]
  15. Tichavsky, P.; Muravchik, C.; Nehorai, A. Posterior Cramer-rao bounds for discrete-time nonlinear filtering. IEEE Trans. Signal Process. 1998, 46, 1386–1396. [Google Scholar] [CrossRef]
  16. Zuo, L.; Niu, R.; Varshney, P.K. Conditional posterior Cramér–rao lower bounds for nonlinear sequential bayesian estimation. IEEE Trans. Signal Process. 2011, 59, 1–14. [Google Scholar] [CrossRef]
  17. Qi, Y. Wireless Geolocation in a Non-Line-of-Sight Environment. Ph.D. Thesis, Princeton University, Princeton, NJ, USA, 2003. [Google Scholar]
  18. Shen, Y.; Win, M.Z. Fundamental limits of wideband localization Part I: A general framework. IEEE Trans. Inf. Theory 2010, 56, 4956–4980. [Google Scholar] [CrossRef]
  19. Tseng, P.H.; Feng, K.T. Geometry-assisted localization algorithms for wireless networks. IEEE Trans. Mob. Comput. 2013, 12, 774–789. [Google Scholar] [CrossRef]
  20. Arias-de Reyna, E.; Djuric, P.M. Indoor localization with range-based measurements and little prior information. IEEE Sens. J. 2013, 13, 1979–1987. [Google Scholar] [CrossRef]
  21. Yang, Z.; Liu, Y. Quality of trilateration: Confidence-based iterative localization. IEEE Trans. Parallel Distrib. Syst. 2010, 21, 631–640. [Google Scholar] [CrossRef]
  22. Wang, J.; Urriza, P.; Han, Y.; Cabric, D. Weighted centroid localization algorithm: Theoretical analysis and distributed implementation. IEEE Trans. Wirel. Commun. 2011, 10, 3403–3413. [Google Scholar] [CrossRef]
  23. Alsindi, N.A.; Alavi, B.; Pahlavan, K. Measurement and Modeling of Ultrawideband TOA-based ranging in Indoor multipath environments. IEEE Trans. Veh. Technol. 2009, 58, 1046–1058. [Google Scholar] [CrossRef]
  24. Arulampalam, M.; Maskell, S.; Gordon, N.; Clapp, T. A tutorial on particle filters for online Nonlinear/non-Gaussian Bayesian tracking. IEEE Trans. Signal Process. 2002, 50, 174–188. [Google Scholar] [CrossRef]
  25. Van Trees, H.L. Optimum Array Processing (Detection, Estimation, and Modulation Theory, Part IV); Wiley-Interscience: Princeton, NJ, USA, 2002; p. 100. [Google Scholar]
  26. Horn, R.A.; Johnson, C.R. Matrix Analysis; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  27. Adler, S.; Schmitt, S.; Kyas, M. Path loss and multipath effects in a real world indoor localization scenario. In Proceedings of the 2014 11th Workshop on Positioning, Navigation and Communication (WPNC), Dresden, Germany, 12–13 March 2014; pp. 1–7.
  28. Hillebrandt, T.; Will, H.; Kyas, M. Quantitative and spatial evaluation of distance-based localization algorithms. In Progress in Location-Based Services; Springer: Berlin, Germany, 2013; pp. 173–194. [Google Scholar]
  29. Guvenc, I.; Chong, C.C.; Watanabe, F. Analysis of a linear least-squares localization technique in LOS and NLOS environments. In Proceedings of the IEEE 65th Vehicular Technology Conference (VTC2007-Spring), Dublin, Republic of Ireland, 22–25 April 2007; pp. 1886–1890.
  30. Yang, Y.; Zhao, Y.; Kyas, M. A non-parametric modeling of Time-of-flight ranging error for indoor network localization. In Proceedings of the 2013 IEEE Global Communications Conference (GLOBECOM), Atlanta, GA, USA, 9–13 December 2013; pp. 189–194.
  31. Schmitt, S.; Will, H.; Aschenbrenner, B.; Hillebrandt, T.; Kyas, M. A reference system for indoor localization testbeds. In Proceedings of the Internatiuonal Conference on Indoor Positioning and Indoor Navigation, IPIN 2012, Sydney, Australia, 13–15 November 2012; pp. 1–4.
Figure 1. Demonstration of 3D-ranging.
Figure 1. Demonstration of 3D-ranging.
Sensors 16 01346 g001
Figure 2. The spatial position error distribution (SPED) results of the CRLB with and without prior information. (a) The SPED of CRLB x t ; (b) the SPED of CRLB x t , k t , l ; (c) the SPED of CRLB x t x t ; (d) the SPED of CRLB x t , k t , l x t , k t , l .
Figure 2. The spatial position error distribution (SPED) results of the CRLB with and without prior information. (a) The SPED of CRLB x t ; (b) the SPED of CRLB x t , k t , l ; (c) the SPED of CRLB x t x t ; (d) the SPED of CRLB x t , k t , l x t , k t , l .
Sensors 16 01346 g002
Figure 3. The evaluation of the impact of k t . (a) v t j N ( 0 , 0 . 5 2 ) ; (b) v t j N ( 0 , 1 2 ) ; (c) v t j N ( 0 , 2 2 ) ; (d) v t j N ( 0 , 3 2 ) ; (e) v t j N ( 0 , 4 2 ) ; (f) v t j N ( 0 , 5 2 ) .
Figure 3. The evaluation of the impact of k t . (a) v t j N ( 0 , 0 . 5 2 ) ; (b) v t j N ( 0 , 1 2 ) ; (c) v t j N ( 0 , 2 2 ) ; (d) v t j N ( 0 , 3 2 ) ; (e) v t j N ( 0 , 4 2 ) ; (f) v t j N ( 0 , 5 2 ) .
Sensors 16 01346 g003
Figure 4. The evaluation of the impact of k t with the multiple noise environment and prior knowledge of k. (a) v t j N ( 0 , 0 . 5 2 ) ; (b) v t j N ( 0 , 1 2 ) ; (c) v t j N ( 0 , 3 2 ) ; (d) v t j N ( 0 , 5 2 ) .
Figure 4. The evaluation of the impact of k t with the multiple noise environment and prior knowledge of k. (a) v t j N ( 0 , 0 . 5 2 ) ; (b) v t j N ( 0 , 1 2 ) ; (c) v t j N ( 0 , 3 2 ) ; (d) v t j N ( 0 , 5 2 ) .
Sensors 16 01346 g004
Figure 5. The impact of the number of l.
Figure 5. The impact of the number of l.
Sensors 16 01346 g005
Figure 6. The simulation of a building layout. (a) Anchor deployment of the playing field; (b) the SPED of CRLB x t , k t , l ; (c) the SPED of CRLB x t , k t , l k t , l ; (d) the SPED of CRLB x t , k t , l x t , k t , l .
Figure 6. The simulation of a building layout. (a) Anchor deployment of the playing field; (b) the SPED of CRLB x t , k t , l ; (c) the SPED of CRLB x t , k t , l k t , l ; (d) the SPED of CRLB x t , k t , l x t , k t , l .
Sensors 16 01346 g006
Figure 7. Sequential estimation lower bound.
Figure 7. Sequential estimation lower bound.
Sensors 16 01346 g007
Figure 8. Cramér–Rao lower bound with different number of anchors.
Figure 8. Cramér–Rao lower bound with different number of anchors.
Sensors 16 01346 g008
Figure 9. Building layout for the indoor localization experiment and the robot trajectory. The triangles mark the positions of sensor nodes, which are placed either in the offices or along the corridor.
Figure 9. Building layout for the indoor localization experiment and the robot trajectory. The triangles mark the positions of sensor nodes, which are placed either in the offices or along the corridor.
Sensors 16 01346 g009
Figure 10. The real system evaluation. (a) SPED without prior knowledge; (b) SPED with prior knowledge.
Figure 10. The real system evaluation. (a) SPED without prior knowledge; (b) SPED with prior knowledge.
Sensors 16 01346 g010
Table 1. Notations and explanations.
Table 1. Notations and explanations.
x t The target position state h t The measurement function
a j The target position vector θ Unknown state vector
k t The relative height vector J ( θ ) Fisher information matrix
l The NLOS drift J P The prior matrix
z t The measurement vector J D The measurement matrix
q t The prediction noise J K Nonlinear auxiliary matrix
v t The measurement noise J L Linear auxiliary matrix
d j ( ) The distance function J S State matrix
f t ( ) The prediction function J A Auxiliary matrix
Table 2. Algorithm comparisons.
Table 2. Algorithm comparisons.
AlgorithmStatePrior StateNLOSPrior NLOSRelative HeightPrior Height
Generalized CRLB
Conditional CRLB
Equivalent CRLB
Our Method

Share and Cite

MDPI and ACS Style

Zhao, Y.; Li, X.; Zhang, S.; Meng, T.; Zhang, Y. Practical Performance Analysis for Multiple Information Fusion Based Scalable Localization System Using Wireless Sensor Networks. Sensors 2016, 16, 1346. https://doi.org/10.3390/s16091346

AMA Style

Zhao Y, Li X, Zhang S, Meng T, Zhang Y. Practical Performance Analysis for Multiple Information Fusion Based Scalable Localization System Using Wireless Sensor Networks. Sensors. 2016; 16(9):1346. https://doi.org/10.3390/s16091346

Chicago/Turabian Style

Zhao, Yubin, Xiaofan Li, Sha Zhang, Tianhui Meng, and Yiwen Zhang. 2016. "Practical Performance Analysis for Multiple Information Fusion Based Scalable Localization System Using Wireless Sensor Networks" Sensors 16, no. 9: 1346. https://doi.org/10.3390/s16091346

APA Style

Zhao, Y., Li, X., Zhang, S., Meng, T., & Zhang, Y. (2016). Practical Performance Analysis for Multiple Information Fusion Based Scalable Localization System Using Wireless Sensor Networks. Sensors, 16(9), 1346. https://doi.org/10.3390/s16091346

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop