Vote selection mechanisms and probabilistic data association-based mobile node localization algorithm in mixed LOS/NLOS environments

As one of the key technologies of wireless sensor networks (WSNs), the localization of mobile nodes (MN) is one of the most significant research topics in WSNs. When a line-of-sight (LOS) channel is available, accuracy localization result can be obtained. Motivated by the fact that the non-line-of-sight (NLOS) propagation of signal is ubiquitous and decreases the accuracy of localization, we propose a MN localization algorithm in mixed LOS/NLOS environments. Considering the characteristics of NLOS error, we propose a localization algorithm based on vote selection mechanisms to filter the distance measurements and reserve the reliable measurements. Then a modified probabilistic data association algorithm is proposed to fuse the multiple measurements reserved from the vote selection. The position of the MN is finally determined by a linear least squares algorithm based on reference nodes selection. This algorithm effectively mitigates various kinds of NLOS errors and largely improves the localization accuracy of the MN in mixed LOS/NLOS environments. The simulation and experiments results show that the proposed algorithm has better robustness and higher localization accuracy than other methods.


Introduction
Wireless sensor networks (WSNs) have attracted significant attention in recent years.Sensors have been used widely in different applications such as environment, health, military, national security and other fields.Recently, localization of mobile nodes (MNs) in WSNs has received considerable attention [1][2][3][4].Global Positioning System and Cellular Positioning System are the most common localization systems, but their performance deteriorate for indoor localization [5].
The common methods to estimate the distance between a transmitter and a receiver are received signal strength indicator, time of arrival, time difference of arrival and angle of arrival [6].Most of the papers assume that the signal propagation is line of sight (LOS) between the MN and the beacon nodes (BN).We can obtain an accurate distance estimation by using the filtering techniques [7].However, the direct propagation path between the MN and the BN may be blocked by obstacles in practical applications of WSNs, such as the indoor and urban environments.Non-line of sight (NLOS) signal propagation is often inevitable due to the reflection and diffraction.The distance measurement error includes the measurement noise and the NLOS propagation error.Therefore the NLOS error is one of the challenges for accurate mobile location estimation in WSNs.
Methods for mitigating negative impact from NLOS in existing localization algorithms fall into two major categories: Some researchers perform LOS reconstruction with the distance measurement.They mitigate the NLOS error of the distance measurements and make the revised measurements closer to the real distance measurements, then apply the localization algorithm to accomplish the locating.Wylie [8] finds the Nth order polynomial fit to smooth the data to correct the error.Based on the premise that standard deviation from the measured value is greater in NLOS than in LOS, this method determines whether the measured value contains NLOS error by comparing standard deviation of the actual measured value.Then the LOS reconstruction is performed on NLOS distance measurements.Communication will lead to increase in energy consumption because the method requires statistical properties of the error distribution.The proposed algorithm could effectively reduce the NLOS error, while the proposed algorithm needs the prior knowledge of the statistical characteristics of the system's standard measurement noise.Wann and Lin [9] proposes a method to fuse two different types of signals combining Bayes rules and fuzzy logic based on [8], which further improves localization accuracy.
Other researchers introduce weighted factor to mitigate the adverse effects from NLOS distance measurements on localization.In order to increase the robustness of the algorithm for NLOS error, the researchers give greater weight to less affected distance measurement and give less weight to greater affected distance measurement.Chan and Ho [10] presents the weighted linear least-squares estimation algorithm based on the Taylor series expansion.Although the algorithm has been used widely, it requires statistical information to determine the weights.Chen [11] proposes a residual weighting algorithm to mitigate the NLOS error.The proposed approach utilizes the sum of squared residuals obtained from least squares estimation as the indicator of the accuracy.Then they apply least squares approximation method and compute the estimated location.The method need not identify the distance measurement signal.Kirubarajan and Bar-Shalom [12] proposes a probabilistic data association with different probabilities to weight each distance measurement.The proposed algorithm could be effective for target tracking and is useful for a wide variety of real tracking problem.
In addition to the methods mentioned above, Wang et al. [13] presents a robust mobile target localization approach.The proposed algorithm utilizes the probability density function (PDF) estimator to identify the NLOS condition, then an improved biased Kalman filter (KF) algorithm is used to mitigate the NLOS error.The proposed algorithm needs the less prior information about communication environments, while has a good location performance even in severe NLOS situations.Mazuelas [14] utilizes known parameters and distributions of the NLOS error, computing the NLOS error from the distance measurements to effectively correct the measurements from NLOS propagation.The proposed algorithm could effectively correct the measurements from NLOS propagation while need the prior NLOS measurement.Chan et al. [5] proposes to utilize the residual test to identify the LOS base station, so the localization can proceed with the LOS base station.In [15], a linear programming approach is introduced to mitigate the NLOS error.By using the geometry of the base station, a constrained nonlinear optimization approach is presented to solve the NLOS localization [16].Wang et al. [17] treats localization in the LOS environment as unconstrained optimization problems.Therefore the NLOS situation localization problem is transformed into an optimization problem to solve.In [18], an NLOS localization algorithm based on semi-definite programming technique is presented, and the localization problem has been approximated by a convex optimization problem.Chen and co-author [19] proposes the interacting multiple model approach with the KF technique to interact the LOS distance measurement and the NLOS distance measurement.However, the algorithm needs to get an NLOS parameters as a priori information.In [20], a geometry-constrained location estimation (GLE) based on the two step least squares algorithm is proposed.The proposed algorithm possesses low complexity from the two-step LS algorithm and is suitable for the NLOS localization at the same time.Abumansoor et al. [21] propose a cooperative localization algorithm to mitigate the NLOS error.Yu and Dutkiewicz [22] proposes the high-frequency distance measurement data processing in base station localization.The proposed algorithm could significantly improve position accuracy.An NLOS identification and probability generation algorithm is proposed by Hammes and Zoubir [23].In this method, the M-estimate based robust KF is used to reduce the NLOS effect and the algorithm yields positioning accuracy similar to the EKF in the LOS environments and even significantly outperforms the REKF in the NLOS environments.Lin et al. [24] proposes a NLOS identification method to identify LOS measurement and NLOS measurement by using the statistics of measurements from the MN, but the algorithm focuses on the cellular network.Guvenc and Chong [1] presents the linear least square algorithm based on fixed BN.The disadvantage is that the localization accuracy will decline if the MN is far away from the BN without changing the BN deployment environment.Hence, localization accuracy is influenced by the selection of BN.Cheng et al. [25] proposes a likelihood matrix correction based mixed Kalman and H-infinity filter method.They propose the likelihood matrix based correction to correct the measurements, then using the mixed Kalman and H-infinity Filter method to improve the range measurement.This method does not need the prior information about the statistical properties of the NLOS errors.But the computational complexity of the method is high.Wang and Zekavat [26] proposes an omnidirectional MN localization technique in NLOS scenarios.They proposed an algorithm to identify if a MN is in the LOS of multiple BN or not.The method enables BN to determine and localize their shared reflection points, and to localize NLOS MN.
In this paper, we propose the voting weighted association reference (VWAR) algorithm and voting weighted association least (VWAL) algorithm.We firstly use data preprocessing method of measurements based on voting selection mechanisms and preserve the reliable data.Secondly, a weight value generation algorithm is proposed.The weight values of each distance measurements can be reasonably generated by the proposed algorithm.Probability data association algorithm is used to weight these distance measurements.Then, we use KF algorithm to filter the weighted distance data.Finally, we employ reference linear least square algorithm based on the selection of reference nodes to estimate the location of MN.
The main contribution of this paper is given as follows: (1) For the NLOS measurements, the vote selection mechanisms are proposed to filter the distance measurements.
Then we propose the probabilistic data association method to mitigate the NLOS error.(2) The proposed algorithm does not assume any statistical knowledge of the NLOS error.Therefore it is independent of the physical measurement ways.
The rest of the paper is organized as follows: In Sect.2, we introduce the architecture of the localization algorithm.Section 3 introduces our proposed strategy in detail.In Sect.4, some simulation and experiment results are presented.The conclusions are given in Sect. 5.

The proposed algorithm architecture
According to whether or not the NLOS error parameters are known, NLOS localization algorithm can be divided into two categories: localization algorithms with known parameters and localization algorithms with unknown parameters.If we can accurately get the NLOS error distribution model and its parameters as prior knowledge, we can mitigate the NLOS error and achieve higher localization performance.But the NLOS propagation model in the actual environment is not consistent with prior knowledge, which decrease the localization accuracy obviously.The unknown parameters localization algorithm does not need to obtain a NLOS distribution model and its parameters as prior knowledge, so this kind of algorithm will not have great influence when the NLOS propagation environments change.The deployment environment of wireless sensor networks is complicated, especially in NLOS propagation environment, thus the unknown parameters localization algorithm has higher research value and practical significance.The proposed algorithm is an unknown parameters localization algorithm and it can guarantee higher localization accuracy in different NLOS propagation environments.
Figure 1 shows the architecture of the proposed algorithm.The proposed algorithm is referenced from the high-frequency distance measurement data processing in base station localization.Ten distance measurements can be acquired by each beacon node in each localization process.In which, r i (t), i = 1, 2, . . ., N is original distance measurement vector of the ith beacon node at time t; ri (t) is the distance vector composed of the reserved data filtered through voting selection method which is initially the distance measurements of each beacon node; N v (t) is the number of the distance measurement through voting; θi , P i Fig. 1 Architecture of the proposed localization algorithm 123 are KF parameter; ri (t) is the distance by using associating probabilistic data association and KF; X(t) is the location vector of the MN through calculating.The proposed algorithm is described as follows: (1) Voting selection mechanisms: Selecting and processing the initial distance measurements of each BN by voting selection mechanisms, then we reserve the reliable distance measurements to form vector ri (t).

NLOS localization based on voting selection mechanisms 3.1 System model
In order to locate the MN, N BN are considered in this paper.The coordinate of the BN are represented as [x i , y i ] T , i = 1, 2, . . ., N ; the coordinate of the mobile node at time t is represented as [x(t), y(t)] T ; the measured distance can be acquired by multiplying the arrival time and the speed of signal.The real distance between the MN and the ith beacon node at the time t is : In the LOS environment, the measured distance between the mobile node and the ith beacon node at the time t is: where e L O S is the measurement noise error [27], which is modeled as zero-mean white Gaussian and standard deviation is σ g .In the NLOS propagation environment, the measured distance at the time t is modeled as: where e N L O S is the NLOS error and independent with e L O S .In different kinds of environments, e N L O S could obey Exponential distribution, Gaussian distribution or Uniform distribution [28].

Measurement data preprocessing method based on voting selection mechanisms
The differences of the mean values and standard deviations between LOS and NLOS distance measurements are used in the measurement data preprocessing method based on voting selection mechanisms.In order to provide the reliable remained data for following data processing, the unreliable part is eliminated by reasonable filtrating of the initial distance measurement.The precondition of using the proposed algorithm is that it is necessary to assume the sampling frequency of distance measurement is larger than 1, which is 10 Hz, and the localization frequency is set as 1 Hz in this paper.The idea is referenced from the high-frequency distance measurement data processing in base station localization [22].
Considering the temporal interval [0, K ] of the whole localization process, L distance measurements are executed by beacon node i. L distance measurements are shown as follows: Then, the set is divided as L − Z + 1 interlaced subsets according to the sampling time sequence.Each subset contains Z distance measurements, concluded within set ψ: where the parameter Z needs to consider the set R j contains enough samples for analyzing the statistical property.In the corresponding temporal interval with Z distance measurements in R j , the moving distance value of MN should not be too large for the robustness of the statistical analysis.Therefore the Z is set to be 10, which is accordant with the distance measurement frequency.We sort the distance measurements in each R j by value.Meanwhile, indexing label of each measurement is remained which is expressed as follows: where R a j is the sorted measurement, r 1 , r 2 , . . ., r Z ∈ , and Selecting the smallest l measurements in R a j , we calculate the mean value μ j (l) and standard deviation σ j (l) .The mean sequence μ a j and standard deviation sequence σ a j of R a j can be obtained: Comparing each estimated standard deviation in standard deviation sequence σ a j in each R j , we find the closest value of estimated standard deviation σm to the measurement error standard deviation σ g .The distance measurement, with which the estimated standard deviation is calculated, can be selected to get one vote.Finally, the voting weight of each distance measurement can be acquired: where α(k) is the voting weight value of distance measurement r k , c k is the final vote of r k , c Z is the probable maximum vote of r k , namely, c Z = Z .We treat the distance measurements with voting weight value α(k) ≥ 0.5 as reliable data, and that with voting weight value α(k) < 0.5 as unreliable data.Then we remain reliable data and eliminate unreliable data in following processing.In this paper, the distance measurement frequency has been set as 10 Hz and the localization frequency is 1 Hz in the proposed algorithm.Then in each localization process, 10 distance measurements can be acquired by each beacon node.Using all of the distance measurements to form an initial distance vector r(t).By voting selection, reserve the reliable distance measurements as vector r(t) for following calculation.
When the beacon node is under the LOS or NLOS situation, all the measurements in R j are more probable to be LOS measurements or NLOS measurements.When under LOS situation, the measurements in R j can be mostly remained for following calculating through the voting selection.When under NLOS situation, the NLOS standard deviation is usually larger than the LOS ones, the NLOS measurements with small errors can be selected for calculating in order to compare with the measure standard deviation.

Probabilistic data association algorithm based on voting selection mechanisms
Using measured data preprocessing method based on voting selection mechanisms, distance vector r(t) which contains reliable measured data can be obtained.After that, we process the distance in vector r(t) with weight value generation and probability data association.The result is filtered through KF.Since the number of BN is N , the number of KFs is also N to filter all data from BN.The state equation corresponding with the beacon node is defined as follows [29]: where θ (t) is the state vector of MN relative to ith beacon node T is white Gaussian noise; matrix A and G are defined as: where k is the period of localization.
The measurement equation of system is: where H is observation matrix, defined as H = [1, 0]; the vector v n (t) is random variable caused by measurement noise of sensors, follows Gaussian distribution whose mean value is 0 and standard deviation is The time update equation of KF can be expressed as follows: where Q is the variance of system processing noise.
Forming the reference of weight generation, and the residual variable is generated as: (16) where r (t|t − 1) is the reference of weight generation; rn (t) is the nth distance measurement in r(t); γ n (t) is the residual error caused by rn (t).
Equation ( 17) can be derived from Eqs. ( 12) and ( 16): where In each localization process, only one distance measurement r j (t) is deemed as the highest accurate distance measurement data.Other data in r(t) can be treated as less reliable clutter data since the existence of LOS or NLOS errors, which follows Uniform distribution in confidence interval.In order to explain the weight generation algorithm, we define as follows: The associated probability can be obtained [30]: Equation ( 20) can be inferred as: where j=1 ; R t is the accumulation set of distance measurements, namely, R t = {R(i)} t i=1 .Since the residual error variables γ j (t), j = 1, 2, . . ., N v (t) are mutually independent, applying Bayes theorem, Eq. ( 21) can be rewritten as: where, c is normalization constant; f (•) is the joint density function for each distance measurement under the condition of δ j (t), N v (t)andR t−1 .For δ j (t), f (•) is the product of the accurate data followed by Gaussian PDF and the clutter data in confidential interval followed by Uniform distribution probability density function.
Therefore, in Eq. ( 25), accurate data is deemed as following Gaussian distribution whose mean value is r (t|t − 1).On the other hand, the clutter data is deemed is modeled by a uniform probability distribution in the confidence interval of λ i (t).The less the distance measurement r j (t) differs from weight generation standard r (t|t − 1), the less the value of residual error γ j (t) is.Hence in order to have more effect in localization calculation with r j (t), it is probable to generate a large weight value w j (t).
The prior probability of Eq. ( 22) is: Multiplying Eq. ( 25) by ( 27) leads to: Normalizing the Eq. ( 28), the weight value r j (t) is processed as: Then weighting each distance measurement in r(t), and combining with KF algorithm.For calculating Kalman gain, firstly define update variance as [12]: And the Kalman gain is: The final state estimated value is obtained from using weighted distance measurement as: The covariance updating equation can be expressed as: where P L (t|t) is the posteriori covariance matrix of standard KF algorithm, and where Above is under the situation of N v (t) > 0, which is relatively normal.However, the situation of N v (t) = 0 is also probable.When N v (t) = 0, it is shown as follows: After calculating the posteriori state estimation vector θ (t|t), the distance value used in localization calculation can be obtained from:

Linear least square algorithm based on the selection of beacon node
We assume that the coordinates of BN are (x 1 , y 1 ), (x 2 , y 2 ), . . ., (x N , y N ), the coordinate of the MN is (x(t), y(t)).
According to the processed distances r1 (t), r2 (t), . . ., rN (t) in the 3.2 and the coordinates of both BN and the MN, the following linear equation can be concluded: where By using linear least square algorithm, the coordinate of the mobile node can be achieved as: The first beacon node is selected as the reference node in Eq. ( 40), thus the linear least square algorithm is based on fixed reference node.The linear least square algorithm based on the selection of reference node is adopted in this paper.When establishing linear equations in Eq. ( 40), the beacon node relevant to the smallest distance measurement is always to be chosen as reference node, which is shown as: Being similar with standard linear least square algorithm, the one based on the selection of reference node is one of LOS localization algorithms.From Eq. ( 42), the distance measurement from beacon node directly influences the selection of reference node.In the NLOS environment, the selection of reference node is effected by inaccurate NLOS distance measurement.However, the initial distance measurement is filtrated by voting selection mechanisms in the proposed algorithm.The remained reliable distance measurement is used in the final localization calculation through processing by probability data association algorithm.Moreover, the negative influence on localization result has been eliminated or at least weakened by filtering unreliable NLOS distance measurement.

Simulation results
In this section, we analyze the performance of the proposed algorithm described in Sect. 3 through simulation and experiment.Figure 2 shows the placement of seven BN in the 100 m × 100 m square area and one MN in this area.The acceleration of the MN is variable.We assume that the communication range of sensor node is 150 m.The measurement noise obeys Gaussian distribution whose mean value is zero and the standard deviation σ g is 3 m [31].The NLOS errors obey Exponential distribution whose mean value μ N L O S = 5 m [32].We compare the proposed VWAR algorithm and VWAL algorithm with the KF Linear Least Squares (KF-LLS) algorithm [29], the H-Infinity filter Linear Least Square (HH-LLS) algorithm [33] and the Particle filter Linear Least Square (PF-LLS) algorithm [34].The simulation results of the proposed algorithm are verified through MAT-LAB R2010a.In each simulation case, 50 Monte Carlo runs are performed with the same parameters.The performance of the proposed algorithm is measured by Root Mean Square Error (RMSE) and Cumulative Distribution Function: where (x j (t), y j (t)), ( x j (t), ŷ j (t)) describes the true and the estimated position of MN at time t for the jth Monte Carlo run, respectively.
N. Hu et al. Figure 3 shows the comparison of the root mean square error for five algorithms.The localization errors of KF-LLS algorithm, HH-LLS algorithm and PF-LLS algorithm are relatively larger than the proposed algorithms.KF-LLS algorithm, HH-LLS algorithm and PF-LLS algorithm are sensitive to NLOS error, these localization methods are easier to be affected by the NLOS error.In contrast, the localization performance of the VWAL algorithm and VWAR algorithm are stable and the localization accuracy improve greatly.It can be seen that VWAR algorithm and VWAL algorithm have strong inhibition for NLOS error, and VWAR algorithm has the highest localization accuracy.
Figure 4 shows the cumulative distribution function of localization error for the five algorithms.It can be observed that the 90-percentile of the KF-LLS algorithm is less than 7.1 m, the 90-percentile of the PF-LLS algorithm is less than 7.4 m and the 90-percentile of the HH-LLS algorithm is less than 7.5 m.The percentiles of the localization error of the VWAR algorithm is 90 % less than 3.7 m and the percentiles of the localization error of the VWAL algorithm is 90 % less than 4.9 m.Therefore, the localization performance of the VWAR algorithm and the VWAL algorithms are far superior to PF-LLS algorithm, KF-LLS algorithm and HH-LLS algorithm, and the VWAR algorithm has slightly higher localization accuracy than the VWAL algorithm.
In the following section, we evaluate the performance of our proposed algorithm when the NLOS errors obey the Uniform distribution, the Gaussian distribution and the Exponential distribution.

The NLOS errors obey uniform distribution
When the NLOS errors obey a Uniform distribution, the minimum and maximum values of parameters of Uniform distribution are 1 and U max .Figure 5 shows the relationship between U max and root mean square error.We observe that the localization errors of the five algorithms increase as the U max increases.When the U max = 2, the VWAR algorithm has higher localization accuracy than HH-LLS, PF-LLS KF-LLS and VWAL algorithms about 32.35, 22.30, 14.81 and 11.54 %, respectively.When the U max = 14, the VWAR algorithm is increased by 42.86, 41.18, 40.87 and 6.85 %.Thus by increasing U max , the localization accuracy of the HH-LLS, PF-LLS and KF-LLS algorithms decrease significantly.Therefore the localization performance of the proposed algorithms is stable and the VWAR algorithm has the higher localization accuracy than the VWAL algorithm.
Figure 6 illustrates the measurement noise standard deviation versus RMSE.We can observe that VWAR algorithm has

The NLOS errors obey Gaussian distribution
When the NLOS errors obey a Gaussian distribution, the performance of each algorithm is discussed in two cases.Figure 7 gives the mean of NLOS error versus RMSE.With the increase of the mean of NLOS error, the localization error increases.In contrast, VWAR algorithm always has the highest localization accuracy.VWAL algorithm is slightly worse than VWAR algorithm, but much higher than the KF-LLS algorithm, HH-LLS algorithm and PF-LLS algorithm.This indicates VWAR algorithm and VWAL algorithm have weak effect by NLOS. Figure 8 shows the relationship between the standard deviation of NLOS error and RMSE.It can be observed that with the increase of the standard deviation of NLOS error, the localization error of KF-LLS, HH-LLS and PF-LLS algorithms increase significantly, while the localization errors of VWAL and VWAR keep descending and finally flatten out.Therefore both VWAR and VWAL can obtain satisfactory localization effects when the standard deviations of NLOS errors are rather large.

The NLOS errors obey exponential distribution
When the NLOS errors obey an Exponential distribution, the performance of each algorithm is discussed in two cases.Figure 9 illustrates the mean of NLOS error versus RMSE.It shows all the localization accuracy increase with the mean of N. Hu et al.  Figure 10 shows the standard deviation of measurement noise versus root mean square error.We can observe when NLOS error obeys Exponential distribution, the localization performance of the KF-LLS algorithm, the HH-LLS algorithm and the PF-LLS algorithm are relatively stable, and less affected by the standard deviation variation of measurement noises, but the localization errors still maintain high.The localization accuracy of the VWAL algorithm and the VWAR algorithm decline with the standard deviation of measurement noise increase, but their accuracy are always much higher than the previous three algorithms.

Experiment results
So as to testify the localization performance of the proposed algorithm under practical environment, this paper designed field experiments to verify our methods.The experimental equipment used for field experiments is chirp spread spectrum (CSS) localization system, mainly consisting of CSS nodes.It is rather convenient to measure distance between nodes with CSS nodes.During locating experiments, BN and target node (i.e., mobile node) are established through switching working modes with CSS nodes.The physical map of CSS node is shown in Fig. 11.
Fig. 12 shows the deployment diagram of the field experiment As the figure shows, we deploy 8 BN in a 5 m × 7 m room, MN moves around a rectangle table in uniform velocity (long tables constitute the obstacles between nodes, which are not shown in this figure) following a rectangle trajectory.Through large quantities of distance measurement experiments, we find that under circumstances of LOS, the measurement error standard deviation of CSS node is around 0.7 m.Considering the experimental region size, the measurement error to this extent will apparently affect the localization performance of this algorithm.Thus, the distance measurement frequency of CSS nodes is set to 20 Hz, and 20 distance measurements are carried out at each sampling site in order to weaken the adverse impacts imposed to Figure 14 compares the root mean square error of five algorithms.We can observe the RMSE of the KF-LLS, HH-LLS and PF-LLS algorithms are large.By contrast, the localization accuracy of the VWAL and VWAR algorithms are much higher than the previous three algorithms, and the localization accuracy of the VWAR algorithm is slightly higher than the VWAL algorithm.
The cumulative distribution function of localization error in the field experiment is shown in Fig. 15.It shows that the localization error of the PF-LLS, KF-LLS, HH-LLS, VWAL algorithms are 3.3, 3.7, 3.4 and 1.7 m when the cumulative distribution function is close to 1, while the localization error To compare with the field experiment, we do the simulation in the same environment.In Figs.16 and 17, it can be obviously observed that the proposed VWAR and VWAL algorithms own higher localization accuracy in comparing with other methods.The trend of the simulation result is similar to the field experiment.However the localization accuracy of the simulation is higher than the field experiment.This is because the error is small in simulation, and the interference factors are varied in field experiment.So as to testify the localization performance of the proposed algorithm under other practical applied environments, we deploy the nodes in a 6.5 m × 6.5 m environment (square room) and a 2.3 m × 8.5 m environment (long corridor).
As the Figs.18 and 19 show, the proposed algorithm can get a better performance in compared with other methods in the 6.5 m×6.5 m environment.The result in the 6.5 m×6.5 m environment is close to the result in 5 m × 7 m environment.In Fig. 19, we show that the localization errors of the PF-LLS, KF-LLS, HH-LLS, VWAL and VWAR algorithms are 3.25, 3.5, 3.39 and 1.35 and 1.1m when the cumulative distribution function is close to 1.
Figures 20 and 21 show the performance of the methods in 2.3 m × 8.5 m environment (long corridor).As shown in Fig. 20, the localization error of the proposed methods are less than 1 m in most situations.Figure 21 shows that VWAR

Conclusion
This paper presents a novel localization algorithm for NLOS mobile node based on voting selection mechanisms.Firstly the algorithm makes use of a measurement data preprocessing method based on voting selection mechanisms in order to filter the initial distance measurement values, thus eliminating the adverse impacts on localization calculation caused by the NLOS distance measurements with large NLOS errors.Then the reserved reliable measurement data errors after voting selection are further mitigated through weighting, probabilistic data association and KF process.Finally, the localization calculation is finished with the filtered distance values using linear least squares method based on reference node selection.The simulation results and the field experiments show that the proposed algorithm has higher localization accuracy and is robust to the NLOS errors in different NLOS environments.

( 2 )
Weight generation: We use KF algorithm to predict state vectors, and utilize Bayes theorem to work out the weights corresponding to all the distance values of distance vector ri (t).(3) Probabilistic data association and update KF: The weight values of each distance measurement in ri (t) can be generated by the proposed algorithm.Then probability data association algorithm is used to weight these distance measurements.Using KF algorithm to filter the weighted distance data, we can obtain the filtered distance measurement ri (t) in real time.(4) Localization: Using the filtered distance measurement ri (t) from step (3) and linear least square algorithm based on the selection of BN to complete the localization calculation.

Fig. 2 Fig. 3
Fig. 2 The schematic diagram of the deployment environment

Fig. 4
Fig.4 The cumulative distribution function of localization error

Fig. 7 Fig. 8
Fig.7 The mean of NLOS errors versus root mean square error

Fig. 9 Fig. 10
Fig.9The mean of NLOS errors versus root mean square error

Fig. 11 Fig. 12
Fig. 11 The physical map of CSS node

Fig. 13 Fig. 14
Fig.13 The schematic diagram of the field experiment deployment environment

Fig. 15 Fig. 16
Fig.15 The cumulative distribution function of localization error in the field experiment

Fig. 17 Fig. 18
Fig.17 The cumulative distribution function of localization error

Fig. 19 Fig. 20
Fig.19 The cumulative distribution function of localization 6.5 m × 6.5 m environment

Fig. 21
Fig. 21 The cumulative distribution function of localization error in the 2.3 m × 8.5 m environment