Next Article in Journal
A High Precision Terahertz Wave Image Reconstruction Algorithm
Next Article in Special Issue
Gaussian Process Regression Plus Method for Localization Reliability Improvement
Previous Article in Journal
Design of a Direction-of-Arrival Estimation Method Used for an Automatic Bearing Tracking System
Previous Article in Special Issue
Distributed Information Compression for Target Tracking in Cluster-Based Wireless Sensor Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Letter

Localisation of Sensor Nodes with Hybrid Measurements in Wireless Sensor Networks  †

1
School of Electronic and Electrical Engineering, University of Leeds, Leeds LS2 9JT, UK
2
Department of Automatic Control and Systems Engineering, Sheffield University, Sheffield S1 3JD, UK
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in “Enhanced hybrid positioning in wireless networks II: AoA-RSS“. In Proceedings of the IEEE International Conference on Telecommunications and Multimedia (TEMU), Heraklion, Greece, 28–30 July 2014.
Sensors 2016, 16(7), 1143; https://doi.org/10.3390/s16071143
Submission received: 25 April 2016 / Revised: 29 June 2016 / Accepted: 13 July 2016 / Published: 22 July 2016
(This article belongs to the Special Issue Scalable Localization in Wireless Sensor Networks)

Abstract

:
Localisation in wireless networks faces challenges such as high levels of signal attenuation and unknown path-loss exponents, especially in urban environments. In response to these challenges, this paper proposes solutions to localisation problems in noisy environments. A new observation model for localisation of static nodes is developed based on hybrid measurements, namely angle of arrival and received signal strength data. An approach for localisation of sensor nodes is proposed as a weighted linear least squares algorithm. The unknown path-loss exponent associated with the received signal strength is estimated jointly with the coordinates of the sensor nodes via the generalised pattern search method. The algorithm’s performance validation is conducted both theoretically and by simulation. A theoretical mean square error expression is derived, followed by the derivation of the linear Cramer-Rao bound which serves as a benchmark for the proposed location estimators. Accurate results are demonstrated with 25%–30% improvement in estimation accuracy with a weighted linear least squares algorithm as compared to linear least squares solution.

1. Introduction

Localisation of wireless devices has become exceedingly important in many applications. These include logistics, robotics and surveillance [1]. Range based approaches are favoured for accurate localisation. Two main groups of techniques for estimating the range between sensor nodes are based on the time of arrival (ToA) and the received signal strength (RSS) approach. Location coordinates of nodes can also be estimated by utilising the angle of the impinging signal, this is known as the angle of arrival (AoA) technique. Individual analysis and optimisation of these techniques have been widely studied [2,3]. For range based models, the localisation problem can be solved via high complexity maximum likelihood techniques [4]. A low complexity linear least squares (LLS) approach has also been proposed for ToA based systems [5], its performance is analysed and enhanced in [6]. Similar approaches are adapted for RSS based localisation [3]. On the other hand, the AoA of the signal can be estimated using an array of antennas as in [7] or a rotating beam of radiation [8], and using techniques such as Multiple Signal Classification [9] or estimation of signal parameters via rotational invariance techniques [10].

Related work

With the increasing demand of high accuracy positioning, hybrid signal based localisation is becoming more and more popular. A highly celebrated hybrid AoA-ToA signal model is proposed in [11] which produces biased estimates of the location vector. An improved and unbiased version of which is presented in [12] where a weighted solution is presented. In [13], an angle based localisation model is presented, where the angles are derived from RSSs from different beacons. The algorithm works on a fingerprinting based approach and will fail to perform when the network scenario changes. In [14], the bearing measurements are utilised together with range-difference measurements obtained from time difference of arrival of the signal from multiple sensors. Together with geometric constraint on measurements errors an improved localisation algorithm is proposed for static target nodes. A two step range and angle based positioning is presented in [15], in which the range estimation in the first step is utilised with differential angle measurement obtained in the second step. A cooperative version of hybrid localisation based on ToA and AoA is proposed in [16], which achieves a very high accuracy at the cost of high computational load. Moreover, non-line of sight components of the signal are detected by an iterative algorithm which is based on the incoming hybrid signals in [17]. In [18], positioning and tracking of people is performed using the extended Kalman filter based on time difference of arrival and AoA. Positioning of people is an important application of WSN and has a vital significance in health care systems [19]. In contrast to the mentioned literature, our work is based on the noise covariance of the hybrid measurements and produces completely unbiased estimates of the unknown vector while considering a realistic assumption of unknown path-loss exponent (PLE) vector.
This paper proposes a LLS estimator based on a hybrid AoA-RSS measurement model which produces unbiased estimates of the unknown vector of location coordinates. Furthermore, in order to improve accuracy, the link quality between the anchor node (AN) and the target node (TN) is considered, which is provided by the covariance matrix. Thus, the covariance matrix is first derived and a weighted linear least squares (WLLS) estimator is proposed. Moreover it is noted that different combinations of ANs perform optimally in different sections of the network. As a result, the network can be decomposed into different zones where each zone corresponds to a unique combination of ANs. Hence, a technique that selects this optimal set of ANs for different zones is proposed. In case of ranging via RSS the correct knowledge of PLE associated with each link is required. In most studies the PLEs are assumed to be known, which is an oversimplification of real conditions. Some recent studies jointly estimate the location coordinates and the PLE for localisation [20,21] for RSS measurements only. However, these studies assume the same PLE for every AN-TN link, which is not a valid assumption for real data. In contrast with [20,21], in this paper, we assume an unknown and a different PLE value for each communication link and we propose a novel PLE estimator, based on the generalised pattern search algorithm.
The main contributions of this paper are as follows:
  • A new unbiased observation model for localisation of static nodes is developed based on hybrid measurements, namely angle-of-arrival and received-signal-strength data.
  • A WLLS framework based on the noise covariance of the signal is presented.
  • The mathematical derivation of unbiasness and unbias constant is given.
  • A two step AN selection technique is presented which further improves the performance.
  • Theoretical results for the mean square error (MSE) are derived.
  • Joint PLE and sensor node coordinates estimation is proposed via generalised pattern search (A dynamic version was presented in [22] for mobile nodes).
  • The linear Cramer-Rao bound (LCRB) is derived for the WLLS algorithm.
  • A more practical scenario for simulation is considered where the TNs are situated inside as well as outside the convex hull defined by ANs.
The rest of the paper is organized as follows: Section 2 presents the problem statement and the unbiased system model. The WLLS algorithm is proposed in Section 3. In Section 4 the two step AN selection strategy is presented and the theoretical MSE for LLS is derived. A PLE estimator via the generalised pattern search is proposed in Section 5. The LCRB is derived in Section 6. Finally, in Section 7, we discuss the simulation results which are followed by conclusions in Section 8.

2. System Model

The following notations are introduced: R n and Z n are the sets of n dimensional real numbers and integers respectively. Also, N ( μ , σ 2 ) and U U m i n , U m a x denotes the normal distribution with mean μ and variance σ 2 and uniform distribution between U m i n and U m a x , respectively. A two dimensional network is considered consisting N ANs with known locations i.e., u i = x i , y i T u i R 2 for i = 1 , , N and a TN which has unknown coordinates i.e., u = x , y T u R 2 . Unlike conventional trilateration, in hybrid systems the AN does not define a circle, but rather defines a line. At one end of the line the AN is situated with known position while the TN is situated at the opposite end for which the coordinates are to be estimated. If the slope (AoA) and the magnitude (RSS) information of this line is available, then the TN coordinates can be easily determined using trigonometric equations. The AN receives a signal with line of sight and non-line of sight components. The line of sight and non-line of sight detection is beyond the scope of this paper. Readers are referred to [23,24] for line of sight/non-line of sight detection and mitigation techniques. This section presents the angle of arrival-received signal strength measurement for localisation which serves as a base for the rest of the paper. Let ( x i , y i ) be the coordinates of ith AN then the x and y coordinates of the TN in the presence of both range and angle estimates are given by [25]
x ^ = x i + d ^ i cos θ ^ i δ i
y ^ = y i + d ^ i sin θ ^ i δ i
where d ^ i , θ ^ i and δ i represent range measurements, angle measurement and the unbiasing constant, respectively. The angle measurement θ ^ i is given by
θ ^ i = arctan y y i x x i + m i + ϕ i
where m i represents the zero mean Gaussian noise in angle estimate i.e., m i N 0 , σ m i 2 and ϕ i represents the angular spread caused by the non-line of sight signal which can be statistically described as a Gaussian random variable and can be calculated from experimental data [26]. In this paper, we restrict our attention to line of sight signals only. The non-line of sight detection/mitigation for angle estimation is studied in [27,28]. Equation (3) can be written in vector form as
θ ^ = f u + m
where θ ^ = θ ^ 1 , , θ ^ N T , f u = arctan y y 1 / x x 1 , , arctan y y N / x x N T , m = m 1 , , m N T is the noise vector and ( . ) T represents the transpose operator.
The range measurement, d ^ i , is extracted from the path-loss L i .
L i = L 0 + 10 α i log 10 d i + w i
where L 0 is the path-loss at reference distance d 0 , normally taken as 1m for indoor scenarios and is dependent on antenna characteristics, w i is the zero mean Gaussian random variable characterizing the shadowing effects i.e., w i N 0 , σ w i 2 , α i represents the PLE associated with ith AN with value range from 2–5, depending on the environment, The observed path-loss z i from d 0 to d i is given as L i L 0 , and can be represented as
z ^ i = γ α i ln d i + w i
where γ = 10 ln 10 . The distance estimates from Equation (6) can be obtained as [29]
d ^ i = d i exp w i γ α i κ i
where κ i is the unbiasing constant for RSS measurement only and is given by κ i = exp σ w i 2 / 2 γ α i 2 . Equation (7) can be represented in vector form as d ^ = d exp 1 γ α i w κ , where d ^ = d ^ 1 , , d ^ N T , d = d 1 , , d N T , κ = κ 1 , , κ N and w = w 1 , , w N T represents the shadowing component vector. The mathematical symbol ⊙ represents the Schur product. The unbiasing constant for AoA-RSS signal is given by
δ i = exp σ m i 2 2 σ w i 2 2 γ α i 2
which can be represented in vector form as δ = δ 1 , , δ N T .
Thus, Equations (1) and (2) can be written in matrix form as
b ^ = Au + q
where
A = diag e N , e N R 2 N × 2 , u = x , y T R 2 × 1 b ^ = b ^ ( x ) , b ^ ( y ) T R 2 N × 1
and e N is a column vector of N ones and q is the noise vector with zero mean vector and covariance C u . In Equation (10), b ^ x and b ^ y are given as
b ^ x = x 1 + d ^ 1 cos θ ^ 1 δ 1 , , x N + d ^ N cos θ ^ N δ N T R N × 1 b ^ y = y 1 + d ^ 1 sin θ ^ 1 δ 1 , , y N + d ^ N sin θ ^ N δ N T R N × 1
Then the LLS solution is given by
u ^ = A b ^
where A is the Moore–Penrose pseudoinverse of matrix A and is given by A = A T A 1 A T .
Justification For δ i :
Without considering the unbiasing constant, Equation (11) produces biased estimates of the unknown vector u . We now show mathematically that the unbiasing constant is imperative for unbias estimation. The bias of the LLS technique is given by
Bias = A ϵ x ϵ y T
where ϵ x = E b ^ x b x and ϵ y = E b ^ y b y . b x and b y represents noise free observation. Then the ith term of ϵ x and ϵ y is given by
ϵ x i = E m i , w i d i exp w i γ α i cos θ i + m i d i cos θ i
ϵ y i = E m i , w i d i exp w i γ α i sin θ i + m i d i sin θ i
where E ( . ) represents the mathematical expectation operation. Equations (13) and (14) are reduced to
ϵ x i = d i cos θ i exp σ m i 2 2 + σ w i 2 2 γ α i 2 1
ϵ y i = d i sin θ i exp σ m i 2 2 + σ w i 2 2 γ α i 2 1
Thus, we use δ i in Equations (1) and (2) to reduce Equations (15) and (16) to zero and consequently Equation (12) to zero. The proof is given in Appendix A.

3. Weighted Linear Least Squares Algorithm

The performance of LLS can be improved by utilising the communication link quality between ANs and the TN. Thus links with larger noise are given small weights as compared with links with small noise. This link quality is provided by the covariance matrix. In this section, we exploit the covariance matrix and propose a WLLS solution which is obtained by minimizing the cost function.
ε W L L S u ^ = b ^ A u T C 1 u b ^ A u
where C u is the covariance matrix given by C u = E m , w b ^ b b ^ b T . The matrix C u can be partitioned into sub-matrices as
C u = C x C x y C x y C y R 2 N × 2 N
The sub-matrices in Equation (18) are given as follows
C x = E m , w b ^ x b x b ^ x b x T R N × N
C y = E m , w b ^ y b y b ^ y b y T R N × N
C x y = E m , w b ^ x b x b ^ y b y T R N × N
Then for AoA-RSS measurement Equations (19)–(21) reduces to Equations (22)–(24) for i = j and to 0 for i j .
C x i i = d i 2 2 κ i + d i 2 2 cos 2 θ i κ ¯ i d i cos θ i 2
C y i i = d i 2 2 κ i d i 2 2 cos 2 θ i κ ¯ i d i sin θ i 2
C x y i i = d i 2 cos θ i sin θ i κ ¯ i 1
C x i j = 0 , C y i j = 0 , C x y i j = 0
where κ i = exp σ w i 2 γ α i 2 + σ m i 2 , κ ¯ i = exp σ w i 2 γ α i 2 σ m i 2 and the notation T i j refers to the element at the ith row and jth column of any matrix T .
The elements of C u in Equation (18) depends on the real values of distances and angles, which are not available. Thus, their estimated values are used to calculate the covariance matrix. Now the WLLS solution can be obtained as follows,
u ^ W L L S = A b
where A = A T C 1 u ^ A 1 A T and b ^ = C 1 u ^ b ^ .
Further performance improvement can be attained by optimal selection of ANs described in the next section.

4. Two Step Optimal AN Selection

Due to the unequal error associated with different ANs, some ANs may actually deteriorate the positioning accuracy. These ANs may be positioned at a large distance from the TN or they may receive signal through multiple paths or it may have a poor geometric dilution of precision. This scenario is more obvious in a network where some TNs are outside the convex hull defined by the ANs. Thus, for different TNs in a network there exists an optimal subset of ANs that will produce better estimates than estimates produced while using all ANs. In this section, we present a two step optimal subset selection scheme. A pre-processing step, called zone detection, selects different subsets of ANs for different TNs followed by localisation using the optimal subset of ANs.
Step I: Zone Detection: During this pre-processing step, the whole network is divided into a grid. The complexity of this step depends on the resolution of the grid and the total number of ANs. However, this step needs to be performed only once. Each point on the grid acts a pseudo-TN. For each of these pseudo-TN, the localisation error is calculated for all combinations of ANs using the theoretical MSE presented in the next subsection. The combination that shows the lowest MSE is selected as an optimal combination of anchor nodes for that point. Thus using this technique a particular combination is selected for different points on the grid. In this way the whole network is divided into different regions called zones, where each zone has its own optimal subset of ANs that shows that the minimum MSE during localisation in the next step.
The theoretical MSE of AoA-RSS signal model: The derivation of the theoretical MSE for LLS is performed as follows.
MSE u ^ = Tr E w , m u ^ u u ^ u T
where u ^ is the estimate of the location vector, u is A b is the ground truth and T r ( . ) represents the trace operator. Equation (26) can be simplified as follows
M S E ( u ) = T r E w , m A b ^ A b A b ^ A b T = T r E w , m A b ^ A b b ^ A T b A T = T r A E w , m b ^ b b ^ b A T = T r A C u ^ A T
Thus in the offline stage, the combination of ANs minimizing Equation (27) for a particular grid point is selected as the optimal AN subset.
Step II: Localisation with optimal combination of ANs: The second step is also two fold. Firstly, a rough estimate of the location of the TN is obtained using all ANs. This rough estimate is necessary to detect the zone where the TN belongs. Once the zone is detected, the location of this TN is refined by localising it again, this time using the optimal combination of ANs for its respective zone.

5. Estimation of Unknown PLE

In order to estimate the distance from the target node to the anchor node, the correct knowledge of PLE associated with each link is necessary. Most of the localisation techniques assume that the PLE is known and same for all links. However, even a small error in the PLE vector produces a significant error in the estimated location. In contrast with the common localisation techniques, in this paper we consider the case when the PLE is unknown and has a different value for every link. A new PLE-generalised pattern search algorithm is proposed. For the observation vector b ^ , given by Equation (10), the cost function Ψ u , α with unknown PLE vector and TN’s coordinates vector u is given by
Ψ u , α = A u b ^ 2
where α is the PLE vector given by, α = α 1 , , α N . In Equation (28), α and u are unknown. The LLS solution to u is given by Equation (11) and after replacing it in Equation (28) gives [30]
Ψ α = b x b y I 2 N A A b x b y T
Equation (29) has only one unknown i.e. the vector α . In Equation (29)
b x = exp z ^ 1 γ α 1 cos θ ^ 1 δ 1 , , exp z ^ N γ α N cos θ ^ N δ N T
and
b y = exp z ^ 1 γ α 1 sin θ ^ 1 δ 1 , , exp z ^ N γ α N sin θ ^ N δ N T
and I 2 N is an identity matrix of dimension 2 N . The solution to Equation (29) is given by
α ^ = arg min α Ψ α
Equation (32) can be solved by a brute force search method which is computationally expensive as the cost function has to be evaluated at all possible values of α. For a large number of nodes the brute force search method becomes impractical. Computationally more efficient the generalise pattern search method is therefore used to minimise Equation (32). The generalised pattern search for minimisation of Equation (32) is presented in the next subsection.

Generalised Pattern Search

The generalised pattern search belongs to a family of derivative-free optimisation techniques. Starting from an initial guess for α 0 [ 2 , 5 ] for most environments [31] and an initial step size Δ 0 , The generalised pattern search iteratively updates α k such that Ψ α k + 1 < Ψ α k , where α k represents the value at kth iteration. Each update evaluates the cost function Equation (32) at a point on the mesh, with the updated point closer to the minimum of Ψ α .
Each iteration consists of a search (optional) and a poll step. At each of these steps the cost function is evaluated on a mesh M k , centred at α k and defined by D , a finite set of direction which positively spans R N . The directions can be chosen by any strategy. However each direction d ¯ j j = 1 , , q must be a product of G z j , where q is the cardinality of D , G R N × N is a non-singular generating matrix which for the present problem is G = 1 ν I N for ν > 1 and z j Z N is an integer vector. These conditions are necessary to the convergence theory [32]. Let Z Z N × q denotes a matrix whose columns are z j j = 1 , , q , then D is represented as the product of GZ and the mesh centred at α k is given by
M k = α k + Δ k Dz : z Z q
At the kth poll, the cost function is evaluated at neighbouring poll points given by P k = α k + Δ k d ¯ , d ¯ D k . Thus at k + 1 th iteration if the cost function value i.e., Ψ α k + 1 is lower than Ψ α k then the step size is increased by Δ k + 1 = ξ Δ k for any scalar ξ > 1 and α k + 1 is accepted i.e., M k + 1 is centred at α k + 1 . Otherwise if Ψ α k + 1 > Ψ α k for all the poll points then the step size is decreased by Δ k + 1 = 1 ξ Δ k and α k + 1 = α k . The algorithm is repeated until a stopping condition is reached e.g., Ψ α k + 1 Ψ α k < τ , where τ is some small value. The generalised pattern search for PLE estimation is presented in Algorithm 1.
Algorithm 1: Generalised Pattern Search
for k = 1 ,
   i. Initialize α 0 2 5 , Δ 0 , τ , ξ , ν .
   ii. Evaluate Ψ α k + 1 with all poll points from poll set α k + Δ k d ¯ , d ¯ D .
   iii-a. If improved poll point is found, accept α k + 1 , set Δ k + 1 = ξ Δ k .
   iii-b. If improved poll point cannot be found, set α k + 1 = α k , set Δ k + 1 = Δ k ξ .
   Repeat until Ψ α k + 1 Ψ α k < τ .
end

6. Linear Cramer-Rao Bound

The Cramer-Rao bound characterizes the best possible accuracy of an unbiased estimator. The conventional localisation Cramer-Rao bound is based on individual readings from ANs. Conversely, the LLS and WLLS formulation is based on observation vector b ^ . In order to lower bound the performance of WLLS, we derive the LCRB in this section. The maximum accuracy of the two dimensional localisation is characterized by the MSE bound:
MSE u I u 11 + I u 22 det I u
where I u is the Fisher information matrix (FIM) whose elements are given by Equation (35) [30].
I u i j = μ u u i C 1 u μ u u j + 1 2 T r C 1 u C u u i C 1 u C u u j
where μ u = x 1 + d 1 cos θ 1 , , x N + d N cos θ N , , y 1 + d 1 sin θ 1 , , y N + d N sin θ N T is the mean of the observation vector. The derivatives are obtained as follows:
μ x = 1 1 , 1 2 , 1 N , 0 1 , 0 2 , 0 N T
μ y = 0 1 , 0 2 , 0 N , 1 1 , 1 2 , 1 N T
C x i i x = x x i κ i + y y i sin 2 θ i + x x i cos 2 θ i κ ¯ i 2 x x i
C x i i y = y y i κ i + y y i cos 2 θ i x x i sin 2 θ i κ ¯ i
C y i i x = x x i κ i y y i sin 2 θ i + x x i cos 2 θ i κ ¯ i
C y i i y = y y i κ i y y i cos 2 θ i x x i sin 2 θ i κ ¯ i 2 y y i
C x y i i x = sin 2 θ i x x i cos 2 θ i y y i κ ¯ i 1
C x y i i y = cos 2 θ i x x i + sin 2 θ i y y i κ ¯ i 1
The derivation of Equation (38) to Equation (43) are given in Appendix B.

7. Simulation Results

In this section, we report simulation results which evaluate the performance of the proposed techniques. A 200 m × 200 m network is considered with subsets of 8 ANs at fixed and known positions. To represent a realistic scenario the target nodes are taken at random locations. The parameters considered are given in Table 1 and the network deployment is shown in Figure 1.
In Figure 2, using Monte Carlo simulation, the Avg. RMSE is plotted obtained with LLS approach and WLLS approach for localisation using hybrid AoA-RSS measurements. The noise in angle estimates is kept fixed at σ m 2 = 4 while the shadowing noise is incremented gradually. Also, the PLE values are incremented from 2 to 4.5. The average RMSE is plotted across different values of shadowing noise variance and PLEs. It is evident from the Figure 2 that WLLS approach (based on the noise covariance matrix) has a better accuracy than LLS approach and showing around 25%–30% better results in terms of Avg. RMSE.
In Figure 3, the network is divided into different zone. Each zone have its own optimal subset of ANs. Each color represents a different combination of ANs. It is evident from Figure 3 that using all ANs for localisation does not show minimum error for all TNs. The optimal combinations of ANs for the zones shown in Figure 3 are given in Table 2.
For the configuration given in Figure 1, the network is divided into different zones, shown in Figure 3, using the theoretical MSE given by Equation (27). The blue and the red curves in Figure 4 represent the performance of the unbiased LLS algorithm using all available ANs and the subsets of ANs determined by the first step of the two step optimal AN selection scheme, respectively. In Figure 4, the Avg. RMSE is for all 30 TNs is plotted against shadowing noise variance on lower and angle noise variance on upper x-axis.
The theoretical MSE is used to demonstrate the system’s performance. In Figure 5, the curves represent the Avg. RMSE obtained by Monte Carlo simulation using different number of ANs. The markers represent the Avg. RMSE using the theoretical MSE expression derived in Section 4. Both performances are plotted across shadowing and angle noise variance which are taken across lower and upper x-axis. Since, the markers coincide with the curves, which shows the accurate prediction of the systems performance and which also proves that the mathematical derivation of the theoretical MSE expression in Section 4 is correct.
Figure 6 demonstrates the performance of the hybrid AoA-RSS signal model based on LLS approach when the PLE vector is estimated via the generalised pattern search. Each AN-TN link is associated with a different PLE, which is taken at random between 2–5. LLS using erroneous PLEs is also plotted for comparison, i.e., α ¯ i = α i + p i where α ¯ i is the erroneous PLE and p i is the error associated with true PLE α i . p i is considered to be zero mean Gaussianly distributed random variable with the standard deviation σ p i i.e., p i N 0 , σ p i 2 . For this simulation σ p i = 0 . 2 is considered. It is observed that even such a small error in the PLE vector produces considerable error in the final estimate of the location of TN, while localisation using estimated PLEs produce considerably better estimates.
Comparison of the LLS, WLLS with its corresponding LCRB for the AoA-RSS measurement model is given in Figure 7. In this case, the PLE is kept fixed at 2.5. The Avg. RMSE of all TNs is plotted across both noise variance. Again it can be seen from Figure 7 that the WLLS outperforms the LLS model and that LCRB tightly bounds WLLS.

8. Conclusions

This paper focuses on hybrid localisation using AoA-RSS measurements. An unbiased LLS estimator is proposed and improvements to the LLS model are achieved by proposing a WLLS algorithm and by developing a methodology for optimal AN selection. A PLE estimation technique using the generalised pattern search is also proposed and a closed form MSE expressions for LLS is derived. Further analysis is done and LCRB is derived that tightly bounds WLLS solution.
It is demonstrated via simulation that the WLLS algorithm performs with 25%–30% better accuracy than the LLS algorithm while the performance is further enhanced with the optimal AN selection. It is also observed that the MSE expression accurately predicts the performance of the LLS technique. Furthermore, the generalised pattern search algorithm considerably improves the performance by accurately estimating the PLEs. Finally, it is shown that the LCRB derived, tightly bound the performance of WLLS estimator. The proposed algorithms are suitable for a number of applications including first aid responders, logistics and fault and fire detection.

Acknowledgments

Lyudmila Mihaylova and Naveed Salman appreciate the support of the UK Engineering and Physical Sciences Research Council (EPSRC) via the project Bayesian Tracking and Reasoning over Time (BTaRoT) grant EP/K021516/1.

Author Contributions

Muhammad W. Khan and Naveed Salman derived the algorithms and wrote the paper. Simulation were performed by Muhammad W. Khan. Lyudmila Mihaylova and Andrew H. Kemp were involved in the theoretical development of the solutions, in editing and structuring the manuscript. All authors were involved in discussing potential solutions, the results and determining the test scenarios to be used and performance validation.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ANAnchor node
TNTarget Node
AoAAngle of Arrival
RSSReceived Signal Strength
ToATime of Arrival
MSEMean Squares Error
Avg. RMSEAverage Root Mean Squares Error
LLSLinear Least Squares
WLLSWeighted Linear Least Squares
PLEPath-loss Exponent
LCRBLinear Cramer-Rao Bound
FIMFisher Information Matrix

Appendix A. Derivation of Unbiasing Constant

The ith term of ϵ x i = E w i , m i b ^ x i b x i , for b ^ x i in which δ i is ignored, is given by,
ϵ x i = E w i , m i x i + d i exp w i γ α i cos θ i + m i x i d i cos θ i for i = 1 , , N
Using sum difference formula cos θ i + m i = cos θ i cos m i sin θ i sin m i .
ϵ x i = d i E w i exp w i γ α i cos θ i E m i cos m i sin θ i E m i sin m i d i cos θ i
Taking the expectations E w i exp w i γ α i = exp σ w i 2 γ α i , E m i cos m i = exp σ m i 2 2 and E m i sin m i = 0 , we obtain
ϵ x i = d i cos θ i exp 0 . 5 σ m i 2 + σ w i 2 2 γ α i 2 1
similarly the ith term of ϵ y i = E w i , m i b ^ y i b y i can be reduced to
ϵ y i = d i sin θ i exp 0 . 5 σ m i 2 + σ w i 2 2 γ α i 2 1

Appendix B. Derivation of FIM

Derivation of Equation (38): Taking derivative of Equation (22) with respect to x
C x i i x = x d i 2 2 κ i + x d i 2 2 cos 2 θ i κ ¯ i x d i cos θ i 2
Utilising product rule
C x i i x = x x i κ i + d i 2 2 x cos 2 θ i + cos 2 θ i x d i 2 2 κ ¯ i x d i cos θ i 2
= x x i κ i + d i 2 sin 2 θ i x θ i + cos 2 θ i x x i κ ¯ i x d i cos θ i 2
= x x i κ i + y y i sin 2 θ i + cos 2 θ i x x i κ ¯ i 2 x x i
This is the required solution and it is obtained from Equation (B2) after taking the following derivatives x θ i = y y i d i 2 , x d i 2 2 = x x i and x d i cos θ i 2 = 2 x x i .
Derivation of Equation (39): Taking derivative of Equation (22) with respect to y.
C x i i y = y d i 2 2 κ i + y d i 2 2 cos 2 θ i κ ¯ i y d i cos θ i 2
C x i i y = y y i κ i + d i 2 2 y cos 2 θ i + cos 2 θ i y d i 2 2 κ ¯ i y d i cos θ i 2
= y y i κ i + d i 2 sin 2 θ i y θ i + cos 2 θ i y y i κ ¯ i y d i cos θ i 2
= y y i κ i + cos 2 θ i y y i sin 2 θ i x x i κ ¯ i
Equation (B6) is obtained from Equation (B5) by using the following derivatives y θ i = x x i d i 2 and y d i cos θ i 2 = 0 .
Derivation of Equations (40) and (41) is similar, other than the fact that x variable is replaced by y.
Derivation of Equation (42): Taking derivative of Equation (24) with respect to x.
C x y i i x = x d i 2 cos θ i sin θ i κ ¯ i 1
which can written after using product rule as
x d i 2 cos θ i sin θ i = d i 2 cos θ i x sin θ i + sin θ i x d i 2 cos θ i = d i 2 cos θ i x sin θ i + d i 2 sin θ i x cos θ i + sin θ i cos θ i x d i 2 = d i 2 cos θ i cos θ i y y i d i 2 + d i 2 sin θ i sin θ i y y i d i 2 + sin θ i cos θ i 2 x x i
Replacing Equation (B8) in Equation (B7) we get Equation (42).
Derivation of Equation (43): Taking the derivative of Equation (24) with respect to y.
C x y i i y = y d i 2 cos θ i sin θ i κ ¯ i 1
The derivatives in Equation (B9) is given as
y d i 2 cos θ i sin θ i = d i 2 cos θ i y sin θ i + sin θ i y d i 2 cos θ i = cos 2 θ i x x i sin 2 θ i x x i + 2 sin θ i cos θ i y y i
Replacing Equation (B10) in Equation (B9) we get Equation (43).

References

  1. Patwari, N.; Ash, J.; Kyperountas, S.; Hero, A.; Moses, R.; Correal, N. Locating the nodes: Cooperative localization in wireless sensor networks. IEEE Signal Process. Mag. 2005, 22, 54–69. [Google Scholar] [CrossRef]
  2. Guvenc, I.; Chong, C.C. A Survey on TOA Based Wireless Localization and NLOS Mitigation Techniques. IEEE Commun. Surv. Tutor. 2009, 11, 107–124. [Google Scholar] [CrossRef]
  3. Salman, N.; Ghogho, M.; Kemp, A.H. Optimized Low Complexity Sensor Node Positioning in Wireless Sensor Networks. IEEE Sens. J. 2014, 14, 39–46. [Google Scholar] [CrossRef]
  4. Patwari, N.; Hero, A.; Perkins, M.; Correal, N.; O’Dea, R. Relative location estimation in wireless sensor networks. IEEE Trans. Signal Process. 2003, 51, 2137–2148. [Google Scholar] [CrossRef]
  5. Caffery, J.J. A new approach to the geometry of TOA location. In Proceedings of the 52nd Vehicular Technology Conference, Boston, MA, USA, 24–28 September 2000; Volume 4, pp. 1943–1949.
  6. Guvenc, I.; Gezici, S.; Sahinoglu, Z. Fundamental limits and improved algorithms for linear least-squares wireless position estimation. Wirel. Commun. Mob. Comput. 2012, 12, 1037–1052. [Google Scholar] [CrossRef]
  7. Chen, H.C.; Lin, T.H.; Kung, H.; Lin, C.K.; Gwon, Y. Determining RF angle of arrival using COTS antenna arrays: A field evaluation. In Proceedings of the IEEE Military Communications Conference, Orlando, FL, USA, 29 October–1 November 2012; pp. 1–6.
  8. Nasipuri, A.; Li, K. A Directionality Based Location Discovery Scheme for Wireless Sensor Networks. In Proceedings of the 1st ACM International Workshop on Wireless Sensor Networks and Applications, Atlanta, GA, USA, 28 September 2002; ACM: New York, NY, USA, 2002; pp. 105–111. [Google Scholar]
  9. Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
  10. Roy, R.; Kailath, T. ESPRIT-Estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef]
  11. Yu, K. 3-D localization error analysis in wireless networks. IEEE Trans. Wirel. Commun. 2007, 6, 3472–3481. [Google Scholar]
  12. Khan, M.; Salman, N.; Kemp, A.H. Enhanced hybrid positioning in wireless networks I: AoA-ToA. In Proceedings of the IEEE International Conference on Telecommunications and Multimedia (TEMU), Heraklion, Greece, 28–30 July 2014; pp. 86–91.
  13. Jiang, J.R.; Lin, C.M.; Lin, F.Y.; Huang, S.T. ALRD: AoA localization with RSSI differences of directional antennas for wireless sensor networks. In Proceedings of the International Conference on Information Society (i-Society), London, UK, 25–28 June 2012; pp. 304–309.
  14. Bishop, A.N.; Fidan, B.; Doğançay, K.; Anderson, B.D.O.; Pathirana, P.N. Exploiting Geometry for Improved Hybrid AOA/TDOA-based Localization. Signal Process. 2008, 88, 1775–1791. [Google Scholar] [CrossRef]
  15. Wang, Y.; Wiemeler, M.; Zheng, F.; Xiong, W.; Kaiser, T. Two-step hybrid self-localization using unsynchronized low-complexity anchors. In Proceedings of the International Conference on Localization and GNSS (ICL-GNSS), Turin, Italy, 25–27 June 2013; pp. 1–5.
  16. Khan, M.; Salman, N.; Kemp, A.H. Cooperative positioning using angle of arrival and time of arrival. In Proceedings of the Sensor Signal Processing for Defence (SSPD), Edinburgh, UK, 8–9 September 2014; pp. 1–5.
  17. Horiba, M.; Okamoto, E.; Shinohara, T.; Matsumura, K. An improved NLOS detection scheme for hybrid-TOA/AOA-based localization in indoor environments. In Proceedings of the IEEE International Conference on Ultra-Wideband (ICUWB), Sydney, Australia, 15–18 September 2013; pp. 37–42.
  18. Lategahn, J.; Muller, M.; Rohrig, C. TDoA and RSS Based Extended Kalman Filter for Indoor Person Localization. In Proceedings of the 78th IEEE Vehicular Technology Conference (VTC Fall), Las Vegas, NV, USA, 2–5 September 2013; pp. 1–5.
  19. Tennina, S.; Di Renzo, M.; Kartsakli, E.; Graziosi, F.; Lalos, A.S.; Antonopoulos, A.; Mekikis, P.V.; Alonso, L. WSN4QoL: A WSN-Oriented Healthcare System Architecture. Int. J. Distrib. Sens. Netw. 2014, 2014, 503417. [Google Scholar] [CrossRef]
  20. Salman, N.; Kemp, A.H.; Ghogho, M. Low Complexity Joint Estimation of Location and Path-Loss Exponent. IEEE Wirel. Commun. Lett. 2012, 1, 364–367. [Google Scholar] [CrossRef]
  21. Li, X. RSS-Based Location Estimation with Unknown Pathloss Model. IEEE Trans. Wirel. Commun. 2006, 5, 3626–3633. [Google Scholar] [CrossRef]
  22. Khan, M.W.; Kemp, A.H.; Salman, N.; Mihaylova, L.S. Tracking of wireless mobile nodes in the presence of unknown path-loss characteristics. In Proceedings of the 18th International Conference on Information Fusion (Fusion), Washington, DC, USA, 6–9 July 2015; pp. 104–111.
  23. Yu, K.; Guo, Y. Statistical NLOS Identification Based on AOA, TOA, and Signal Strength. IEEE Trans. Veh. Technol. 2009, 58, 274–286. [Google Scholar] [CrossRef]
  24. Guvenc, I.; Chong, C.C.; Watanabe, F. NLOS Identification and Mitigation for UWB Localization Systems. In Proceedings of the IEEE Wireless Communications and Networking Conference, Hong Kong, China, 11–15 March 2007; pp. 1571–1576.
  25. Salman, N.; Khan, M.W.; Kemp, A.H. Enhanced hybrid positioning in wireless networks II: AoA-RSS. In Proceedings of the IEEE International Conference on Telecommunications and Multimedia (TEMU), Heraklion, Greece, 28–30 July 2014; pp. 92–97.
  26. Tang, H.; Park, Y.; Qiu, T. A TOA-AOA-based NLOS Error Mitigation Method for Location Estimation. EURASIP J. Adv. Signal Process. 2008, 2008, 1–14. [Google Scholar] [CrossRef]
  27. Zhang, Y.; Brown, A.K.; Malik, W.Q.; Edwards, D.J. High Resolution 3-D Angle of Arrival Determination for Indoor UWB Multipath Propagation. IEEE Trans. Wirel. Commun. 2008, 7, 3047–3055. [Google Scholar] [CrossRef]
  28. Jiang, H.; Wang, S.-X. Azimuth and elevation estimation for multipath signals exploiting cyclostationarity and temporal smoothing technology. In Proceedings of the IEEE International Symposium on Microwave, Antenna, Propagation and EMC Technologies for Wireless Communications, Beijing, China, 8–12 August 2005; Volume 2, pp. 1066–1070.
  29. Salman, N.; Ghogho, M.; Kemp, A.H. On the Joint Estimation of the RSS-Based Location and Path-loss Exponent. IEEE Wirel. Commun. Lett. 2012, 1, 34–37. [Google Scholar] [CrossRef]
  30. Kay, S.M. Fundamentals of Statistical Signal Processing: Estimation Theory; Prentice Hall, Inc.: Upper Saddle River, NJ, USA, 1993. [Google Scholar]
  31. Rappaport, T. Wireless Communications: Principles and Practice; Prentice-Hal: Englewood Cliffs, NJ, USA, 1996. [Google Scholar]
  32. Audet, C.; Dennis, J.E. Analysis of Generalized Pattern Searches. SIAM J. Optim. 2000, 13, 889–903. [Google Scholar] [CrossRef]
Figure 1. Network deployment with 30 target nodes (TNs) positioned at random unknown locations and 8 anchor nodes (ANs) at fixed known locations.
Figure 1. Network deployment with 30 target nodes (TNs) positioned at random unknown locations and 8 anchor nodes (ANs) at fixed known locations.
Sensors 16 01143 g001
Figure 2. Performance comparison between linear least squares (LLS) and weighted linear least squares (WLLS) for hybrid angle of arrival (AoA)-received signal strength (RSS) measurement. σ m 2 = 4 0 , ANs = 1 8 , α i = 2 . 5 i , = 2500 .
Figure 2. Performance comparison between linear least squares (LLS) and weighted linear least squares (WLLS) for hybrid angle of arrival (AoA)-received signal strength (RSS) measurement. σ m 2 = 4 0 , ANs = 1 8 , α i = 2 . 5 i , = 2500 .
Sensors 16 01143 g002
Figure 3. Division of network into different zone based on the theoretical mean square error (MSE). ANs = 2 , 4 , 6 , 8 , α i = 2 . 5 i .
Figure 3. Division of network into different zone based on the theoretical mean square error (MSE). ANs = 2 , 4 , 6 , 8 , α i = 2 . 5 i .
Sensors 16 01143 g003
Figure 4. Performance comparison in terms of Avg. RMSE, using optimal subsets of ANs and using all ANs simultaneously. ANs = 2 , 4 , 6 , 8 , = 1000 , α i = 2 . 5 i .
Figure 4. Performance comparison in terms of Avg. RMSE, using optimal subsets of ANs and using all ANs simultaneously. ANs = 2 , 4 , 6 , 8 , = 1000 , α i = 2 . 5 i .
Sensors 16 01143 g004
Figure 5. Performance evaluation via theoretical MSE expression and simulation for LLS. ANs = 2 , 4 , 6 , 8 , 1 , 2 , 3 , 5 , 6 , 7 , 1 8 , = 1500 , α i = 2 . 5 i .
Figure 5. Performance evaluation via theoretical MSE expression and simulation for LLS. ANs = 2 , 4 , 6 , 8 , 1 , 2 , 3 , 5 , 6 , 7 , 1 8 , = 1500 , α i = 2 . 5 i .
Sensors 16 01143 g005
Figure 6. Avg. RMSE comparison using estimated PLEs and true PLE’s. ANs = 1 8 , = 2000 , τ = 1, ξ = 2 , Δ 0 = 0 . 5 , v = 10 , α i U 2 , 5 , α 0 U 2 , 5 , σ p = 0 . 2 .
Figure 6. Avg. RMSE comparison using estimated PLEs and true PLE’s. ANs = 1 8 , = 2000 , τ = 1, ξ = 2 , Δ 0 = 0 . 5 , v = 10 , α i U 2 , 5 , α 0 U 2 , 5 , σ p = 0 . 2 .
Sensors 16 01143 g006
Figure 7. Performance comparison between LLS, WLLS and LCRB using hybrid AoA-RSS measurements. ANs = 1 8 α i = 2 . 5 i , = 2000 .
Figure 7. Performance comparison between LLS, WLLS and LCRB using hybrid AoA-RSS measurements. ANs = 1 8 α i = 2 . 5 i , = 2000 .
Sensors 16 01143 g007
Table 1. Parameters description.
Table 1. Parameters description.
S.NoSymbolDescription
1 σ m 2 Angle noise variance
2 σ w 2 Shadowing noise variance
3 α i PLE associated with ith link
4 α 0 Initial PLE assumption (for initialising GenPS)
5 σ p Standard deviation of erroneous PLE
6 0 Initial step size in GenPS
7 k Step size at kth iteration
8ξStep size indicator in GenPS
9τStopping criteria for GenPS
10Number of iterations
Table 2. Optimal combinations of ANs for zones shown in Figure 3.
Table 2. Optimal combinations of ANs for zones shown in Figure 3.
ZonesOptimal AN Combination
Zone 1 i = 2 , 8
Zone 2 i = 2 , 4
Zone 3 i = 6 , 8
Zone 4 i = 4 , 6
Zone 5 i = 2 , 4 , 8
Zone 6 i = 2 , 6 , 8
Zone 7 i = 2 , 4 , 6
Zone 8 i = 4 , 6 , 8
Zone 9 i = 2 , 4 , 6 , 8

Share and Cite

MDPI and ACS Style

Khan, M.W.; Salman, N.; Kemp, A.H.; Mihaylova, L. Localisation of Sensor Nodes with Hybrid Measurements in Wireless Sensor Networks . Sensors 2016, 16, 1143. https://doi.org/10.3390/s16071143

AMA Style

Khan MW, Salman N, Kemp AH, Mihaylova L. Localisation of Sensor Nodes with Hybrid Measurements in Wireless Sensor Networks . Sensors. 2016; 16(7):1143. https://doi.org/10.3390/s16071143

Chicago/Turabian Style

Khan, Muhammad W., Naveed Salman, Andrew H. Kemp, and Lyudmila Mihaylova. 2016. "Localisation of Sensor Nodes with Hybrid Measurements in Wireless Sensor Networks " Sensors 16, no. 7: 1143. https://doi.org/10.3390/s16071143

APA Style

Khan, M. W., Salman, N., Kemp, A. H., & Mihaylova, L. (2016). Localisation of Sensor Nodes with Hybrid Measurements in Wireless Sensor Networks . Sensors, 16(7), 1143. https://doi.org/10.3390/s16071143

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