Next Article in Journal
The Application of Classical Control in the Design and Analysis of Power Amplifiers for Driving Piezoelectric Stack Actuators
Next Article in Special Issue
Indoor 3D Localization Scheme Based on BLE Signal Fingerprinting and 1D Convolutional Neural Network
Previous Article in Journal
Security Framework for IoT Based Real-Time Health Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Switching Extended Kalman Filter Bank for Indoor Localization Using Wireless Sensor Networks

Department of Electrical Engineering, Wonkwang University, Iksan 54538, Korea
Electronics 2021, 10(6), 718; https://doi.org/10.3390/electronics10060718
Submission received: 17 February 2021 / Revised: 13 March 2021 / Accepted: 17 March 2021 / Published: 18 March 2021
(This article belongs to the Special Issue Indoor Localization Using Wireless Sensor Networks)

Abstract

:
This paper presents a new filtering algorithm, switching extended Kalman filter bank (SEKFB), for indoor localization using wireless sensor networks. SEKFB overcomes the problem of uncertain process-noise covariance that arises when using the constant-velocity motion model for indoor localization. In the SEKFB algorithm, several extended Kalman filters (EKFs) run in parallel using a set of covariance hypotheses, and the most probable output obtained from the EKFs is selected using Mahalanobis distance evaluation. Simulations demonstrated that the SEKFB can provide accurate and reliable localization without the careful selection of process-noise covariance.

1. Introduction

Localization involves tracking locations of objects that are of interest, such as robots, humans, vehicles, and equipment [1,2,3]. Outdoor localization usually depends on a global positioning system, whereas indoor localization exploits various types of wireless sensor networks (WSNs) [4]. A typical WSN comprises several fixed nodes installed at designated locations and mobile nodes attached to target objects. Fixed and mobile nodes communicate with each other using wireless signals, and the locations of mobile nodes are obtained by analyzing the parameters of wireless signals. Time of arrival (TOA) [5], time difference of arrival (TDOA) [6], and angle of arrival [7] are typical wireless measurements used for indoor localization.
A wireless measurement is corrupted by noise, which is inevitable but can be mitigated using stochastic filters (also referred to as state estimators). To exploit stochastic filters, a localization system should be modeled in a state space form that comprises motion and measurement models. Because wireless measurements of WSNs are typically represented by nonlinear measurement models, nonlinear stochastic filters such as the extended Kalman filter (EKF) and the particle filter (PF) are often used for indoor localization [8,9,10,11]. In this study, the EKF was used for indoor localization because it has an advantage over the PF in terms of real-time processing.
In indoor localization, the constant-velocity (CV) motion model is commonly used to represent the motion of target objects. In the CV motion model, process-noise covariance Q plays a critical role, but it is a highly uncertain design parameter [12,13,14,15]. The conventional EKF algorithm uses constant Q values that are selected on the basis of an engineer’s knowledge and experience. However, inappropriately selected Q may degrade localization accuracy [12,13].
Uncertainty in the motion model is one of the oldest problems in the field of target tracking. To solve this problem, interacting multiple model (IMM) filtering [16] was developed. IMM filtering mixes models for which mixing probabilities are computed using the probabilities of transition between multiple models. Transitional probability (TP) is a key design parameter in IMM filtering. Guidelines to designing TP were presented for typical motion models used in target tracking, for example, CV motion and coordinated turn motion models [14,15]. However, TP design is still a cumbersome problem. In particular, it is difficult to design a TP between different Q values in the same CV motion model. Thus, a filtering algorithm that is simple and easy to implement is required for indoor localization.
This paper proposes a new filtering algorithm, the switching extended Kalman filter bank (SEKFB), to overcome the uncertain process-noise covariance problem. In the SEKFB algorithm, a filter bank consisting of several EKFs is constructed. EKFs in the filter bank use different Q hypotheses that are roughly selected without careful consideration or experience. In addition, the most probable among EKF outputs is selected using Mahalanobis distance evaluation [17]. SEKFB thereby alleviates the problem of selecting Q and designing the TP. Without the need for the careful selection of Q , the SEKFB algorithm performs accurate localization compared to the best achieved performance of an EKF. In addition, the SEKFB provides reliable localization, while the conventional EKF exhibits localization failures. Moreover, the SEKFB is easy to implement because of its simple algorithm. The excellent performance of the SEKFB is demonstrated via the simulation of indoor localization using a WSN.
The remainder of this paper is organized as follows. Section 2 discusses the indoor localization scheme utilizing the EKF, and proposes the SEKFB. Section 3 presents simulation results for the demonstration of SEKFB performance. Lastly, the conclusions of the study are presented in Section 4. The abbreviations used in this paper are listed in Table 1.

2. Indoor Localization Scheme and Proposed Algorithm

We considered an indoor localization system using a WSN described as follows. To simplify the problem, the indoor space was assumed to be a two-dimensional (2D) floor space. Figure 1 shows the configuration of a WSN using TDOA measurements.
Four receivers were installed at designated locations in the indoor space. A transmitter (mobile tag) was attached to a target object (e.g., vehicles, equipment, or human) in the space. The four receivers received signals from the transmitter and measured four TOAs that are defined as follows:
TOA i , k = 1 c d i , k , i = 1 , 2 , 3 , 4 ,
d i , k = ( x k x i ) 2 + ( y k y i ) 2 ,
where c is the speed of light; d i is the distance between transmitter and i-th receiver; ( x k , y k ) are the 2D positions of the transmitter at time k; and ( x i , y i ) are the fixed positions of the receivers. The receiver’s clocks were synchronized with each other. The measured TOAs were transmitted to a server computer, and the TDOA was computed as follows:
TDOA i , k = TOA 1 , k TOA i + 1 , k , i = 1 , 2 , 3 .
Because TDOA measurements are corrupted by noise, the EKF was used to estimate the target positions from noisy measurements. At each time step k, the EKF performs two processes, time and measurement updates. In the time-update process, state variables are updated using a motion model. We used the CV motion model for updating, and state variables to be estimated were 2D positions ( x , y ) and velocities ( v x , v y ) as shown in Figure 2.
State variables at time k are represented as ( x k , y k ) and ( v x , k , v y , k ) . To simplify the notation of velocities, we used ( x ˙ k , y ˙ k ) instead of ( v x , k , v y , k ) . Thus, state vector x k was constructed as x k [ x k y k x ˙ k y ˙ k ] T . Thus, the CV motion model [14,15] is represented as
x k + 1 = A x k + G w k , A = 1 0 T 0 0 1 0 T 0 0 1 0 0 0 0 1 , G = T 2 / 2 0 0 T 2 / 2 T 0 0 T ,
where T is the sampling interval, and w k is the process-noise vector. We assumed that w k was zero-mean white Gaussian noise with covariance matrix Q = q I 2 , where I 2 is a 2 × 2 identity matrix.
The measurement update process of the EKF uses a measurement model. The TDOA measurements are expressed as
z 1 , k z 2 , k z 3 , k = h 1 , k h 2 , k h 3 , k = 1 c d 1 d 2 d 1 d 3 d 1 d 4 ,
where z 1 , k , z 2 , k , and z 3 , k are TDOA measurements (in units of nanoseconds); h 1 , k , h 2 , k , and h 3 , k are nonlinear measurement equations that produce TDOA measurements. The measurement vector and measurement equation vector were constructed as z k = [ z 1 , k z 2 , k z 3 , k ] T and h k = [ h 1 , k h 2 , k h 3 , k ] T , respectively. The measurement model can then be expressed as
z k = h k ( x k ) + v k ,
where v k is the Gaussian measurement noise vector with covariance R .
Given the motion and measurement models, the EKF was used to estimate the 2D positions of the transmitter. EKF equations for the state-space models (4) and (6) are represented as follows:
x ^ k = A x ^ k 1 + ,
P k = A P k 1 + A T + G Q G T ,
K k = P k H T ( H P k H T + R ) 1 ,
x ^ k + = x ^ k + K k ( z k H x ^ k ) ,
P k + = ( I K k H ) P k ,
where P k is the estimation error covariance matrix; K k is the Kalman gain; and superscripts − and + represent a priori and a posteriori, respectively. H is the Jacobian matrix defined as
H k = h k x x ^ k ,
and obtained as follows:
H k = h 11 , k h 12 , k 0 0 h 21 , k h 22 , k 0 0 h 31 , k h 32 , k 0 0 , h 11 , k = 1 c x k x 1 d 1 x k x 2 d 2 , h 12 , k = 1 c y k y 1 d 1 y k y 2 d 2 , h 21 , k = 1 c x k x 1 d 1 x k x 3 d 3 , h 22 , k = 1 c y k y 1 d 1 y k y 3 d 3 , h 31 , k = 1 c x k x 1 d 1 x k x 4 d 4 , h 32 , k = 1 c y k y 1 d 1 y k y 4 d 4 .
Process-noise covariance Q in (8) had a significant effect on EKF performance. Q was derived from the CV motion model and is related to the movement rate of a target. Because the movement of a target (e.g., person) is unpredictable, selecting a suitable Q is difficult. If an unsuitable Q value is used, EKF accuracy degrades. Thus, we propose the SEKFB algorithm, which can provide accurate and reliable localization without selecting a suitable Q . Figure 3 illustrates the structure of SEKFB. EKFs using different process noise covariances constitute a filter bank. State estimates of EKFs are evaluated using Mahalanobis distance. The best estimate is selected as the output of SEKFB, and it is also used to reset the EKFs.
The first step of the SEKFB algorithm is to select Q hypotheses, which involves selecting q hypotheses because we assumed that Q = q I 2 . We constructed set of hypotheses { q 1 , q 2 , , q n } , where n is the number of hypotheses. Hypotheses do not need to be selected in a careful manner. Designers can roughly select the hypotheses (e.g., q = 1 , 10 , 100 , ), and the hypotheses perform suitably owing to the SEKFB algorithm.
In the next step, n EKFs using the n hypotheses are operated in parallel. At each time step, n state estimates { x ^ k + , 1 , x ^ k + , 2 , , x ^ k + , n } are obtained. The best state estimate is selected by minimizing Mahalanobis distance between actual TDOA measurement z ^ k and predicted measurements z ^ k j , where j = 1 , 2 , , n . The predicted measurement is computed using the state estimate as follows:
z ^ k j = h k ( x ^ k + , j ) ,
where x ^ k + , j is the state estimate obtained from the j-th EKF. Mahalanobis distance [17] is computed as
D k j = ( z k z ^ k j ) T R 1 ( z k z ^ k j ) .
If the j-th EKF produces the minimal D k j , x ^ k + , j is chosen as the output of the SEKFB at time k. For the estimation of the next time step k + 1 , the other EKFs except the j-th EKF are reset using the information of the most reliable EKF. Thus, EKFs can produce more reliable estimates than those estimated using previous algorithms. The overall algorithm of the SEKFB is summarized in Algorithm 1.
Algorithm 1: Filtering using SEKFB.
Electronics 10 00718 i001

3. Simulation

We simulated indoor localization using SEKFB. The simulation scenario was as follows. A person equipped with a mobile tag (transmitter) traveled along a square trajectory in an indoor space as shown in Figure 4. The size of the indoor space was 20 × 20 m. Four receivers were installed at positions ( 0 , 0 ) , ( 0 , 20 ) , ( 20 , 0 ) , and ( 20 , 20 ) m. The receivers obtained wireless signals from the transmitter, and three TDOA measurements were acquired as shown in Figure 1.
At each time step, the 2D positions of a person ( x k , y k ) were estimated using both the SEKFB and conventional EKFs. Motion and measurement models given in (4) and (6), respectively, were used for the filters. In the simulation, TDOA measurements were generated using the measurement model given in (6), and measurement-noise covariance R = 0 . 1 2 I 3 . We assumed that measurement-noise covariance was exactly known. However, process-noise covariance Q was unknown. Thus, we roughly selected q hypotheses as { 100 , 10 , 1 , 0.1 , 0.01 } . We operated the EKF using the five q hypotheses to find the best performing hypothesis. Next, we operated the SEKFB and compared its performance with the best achievable performance of the EKF.
Filter performance was evaluated by the root-mean-square position error (RMSPE) and root time-averaged mean square error (RTAMSE). We ran 100 Monte Carlo (MC) simulations. RMSPE at time k and RTAMSE were calculated as
RMSPE k = 1 N e m = 1 N e ( x k m x ^ k m ) 2 + ( y k m y ^ k m ) 2 ,
RTAMSE = 1 N e t max m = 1 N e k = 1 t final ( x k m x ^ k m ) 2 + ( y k m y ^ k m ) 2 ,
where t final is the final time step, N e is the number of effective simulations, and N e = 100 N f , where N f is the number of localization failures. A case that produced the RTAMSE at a value greater than 5 m was considered to be a localization failure.
Figure 5a shows the RMSPEs of the EKFs using the five q hypotheses. q = 0.01 produced the largest RMSPE, that is, the worst performance. q = 10 and q = 100 produced much smaller RMSPEs compared with the other q hypotheses. Because RMSPEs produced by q = 10 and q = 100 were similar and difficult to distinguish, we compared their RTAMSEs. Table 2 compares the RTAMSEs of the filters; q = 10 and q = 100 produced RTAMSE values of 0.0252 and 0.0260 , respectively. Thus, the best EKF performance could be achieved with q = 10 .
In Figure 5b, we compared the SEKFB and EKF with q = 10 . The two filters exhibited similar RMSPEs, as shown in Figure 5b, and we compared them in terms of the RTAMSE. As shown in Table 2, the RTAMSE of the SEKFB was 0.0249 , which was smaller than that of the EKF with q = 10 . The best choice for q was unknown in advance. The process-noise covariance of the CV model is a highly uncertain design parameter, and selecting an appropriate covariance value is difficult. EKF performance with q = 10 or q = 100 could be obtained when a selected covariance value is appropriate. However, the SEKFB exhibited accurate localization performance without selecting an appropriate covariance value.
We compared the filters in terms of the number of localization failures as shown in Table 2. EKFs using q = 10 and 100 produced RTAMSEs similar to those of the SEKFB, but they exhibited localization failures. However, the SEKFB did not fail, which means that the SEKFB was more reliable than EKFs using constant q values are. Thus, simulation results demonstrated that SEKFB is more accurate and reliable than a conventional EKF that uses constant covariance is.
Lastly, we compared the SEKFB with the IMM EKF under various signal-to-noise ratio (SNR) conditions. The IMM EKF algorithm ran five EKFs in parallel. Each EKF matched to the CV motion models using different process-noise covariances. The same q values were used for the IMM EKF as those used for the SEKFB. The IMM algorithm required elaborate TP design, but information on TP between the five motion models was completely unknown. Hence, we assumed that the five CV motion models had the same probabilities. Figure 6 compares the RTAMSEs of the SEKFB, IMM EKF, and EKF using the best q value. For various SNR conditions, we used five different measurement-noise covariances as R = r 2 I 3 , where r = 0.5 , 0.25 , 0.1 , 0.05 , and 0.025 . SNR increased as r decreased. In Figure 6, RTAMSEs tended to decrease as r decreased, which indicates that localization accuracy increased as SNR increased. The SEKFB exhibited smaller RTAMSEs than the IMM EKF did because the SEKFB selects the most suitable q at each time step on the basis of Mahalanobis distance, but the IMM EKF cannot. Under high SNR conditions, for example, r = 0.05 and 0.025 , the SEKFB was more accurate than the EKF using best constant q was because the adaptation of q in the SEKFB algorithm was based on Mahalanobis distance evaluation. When computing the Mahalanobis distance, predicted measurement z ^ k was obtained while ignoring measurement noise. Errors due to ignorance may be small under low-measurement-noise (i.e., high SNR) conditions. Thus, the SEKFB is more effective under high SNR conditions.

4. Conclusions

This paper proposed SEKFB for indoor localization using a WSN. The SEKFB switches multiple Q hypotheses in the CV motion model, whereas conventional algorithms use constant Q values. In the simulation, the SEKFB exhibited excellent localization accuracy without a priori information on the best covariance value. SEKFB accuracy was better than that of the IMM EKF. Moreover, the SEKFB exhibited better reliability than that of the conventional EKF in terms of localization failure. The SEKFB can thus provide more accurate and reliable performance than conventional EKFs can for indoor localization using WSNs.

Funding

This paper was supported by Wonkwang University in 2020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sinha, R.S.; Hwang, S.-H. Improved RSSI-based data augmentation technique for fingerprint indoor localization. Electronics 2020, 9, 851. [Google Scholar] [CrossRef]
  2. Hu, X.; Luo, A.; Jiang, W. AGV localization system based on ultra-wideband and vision guidance. Electronics 2020, 9, 448. [Google Scholar] [CrossRef] [Green Version]
  3. Pak, J.M.; Ahn, C.K.; Shmaliy, Y.S.; Shi, P.; Lim, M.T. Accurate and reliable human localization using composite particle/FIR filtering. IEEE Trans. Hum. Mach. Syst. 2017, 47, 332–342. [Google Scholar] [CrossRef]
  4. Liu, H.; Darabi, H.; Banerjee, P.; Liu, J. Survey of wireless indoor positioning technique and systems. IEEE Trans. Syst. Man Cybern. 2007, 37, 1067–1080. [Google Scholar] [CrossRef]
  5. Xu, C.; Ji, M.; Qi, Y.; Zhou, X. MCC-CKF: A distance constrained Kalman filter method for indoor TOA localization applications. Electronics 2019, 8, 478. [Google Scholar] [CrossRef] [Green Version]
  6. Su, C.; Liu, Y.; Liu, L.; Yang, M.; Zhao, H.; Yin, X. Experimental evaluation of multipath mitigation in TDOA-based indoor passive localization system using a beam steering circular polarization antenna. Electronics 2018, 7, 362. [Google Scholar] [CrossRef] [Green Version]
  7. Taponecco, L.; D’Amico, A.A.; Mengali, U. Joint TOA and AOA estimation for UWB localization applications. IEEE Trans. Wirel. Commun. 2011, 10, 2207–2217. [Google Scholar] [CrossRef]
  8. Martin, J.S.; Cortes, A.; Zamor-Cadenas, L.; Svensson, B.J. Precise positioning of autonomous vehicles combining UWB-ranging estimations with on-board sensors. Electronics 2020, 9, 1238. [Google Scholar] [CrossRef]
  9. Lipka, M.; Sippel, E.; Vossiek, M. An extended Kalman filter for direct, real-time, phase-based high precision indoor localization. IEEE Access 2019, 7, 25288–25297. [Google Scholar] [CrossRef]
  10. Pak, J.M.; Ahn, C.K.; Shi, P.; Shmaliy, Y.S.; Lim, M.T. Distributed hybrid particle/FIR filtering for mitigating NLOS effects in TOA-based localization using wireless sensor networks. IEEE Trans. Ind. Electron. 2017, 64, 5182–5191. [Google Scholar] [CrossRef]
  11. Yang, P.; Wu, W. Efficient particle filter localization algorithm in dense passive RFID tag environment. IEEE Trans. Ind. Electron. 2014, 61, 5641–5651. [Google Scholar] [CrossRef]
  12. Pak, J.M.; Ahn, C.K.; Shi, P.; Shmaliy, Y.S.; Lim, M.T. Self-recovering extended Kalman filtering algorithm based on model-based diagnosis and resetting using an FIR filter. Neurocomputing 2016, 173 Pt 3, 645–658. [Google Scholar] [CrossRef]
  13. Pak, J.M.; Kim, P.S.; You, S.H.; Lee, S.S.; Song, M.K. Extended least square unbiased FIR filter for target tracking using the constant velocity motion model. Int. J. Control Autom. Syst. 2017, 15, 947–951. [Google Scholar] [CrossRef]
  14. Blackman, S.; Popoli, R. Design and Analysis of Modern Tracking Systems; Arctech House: Norwood, MA, USA, 1999. [Google Scholar]
  15. Bar-Shalom, y.; Li, X.R.; Kirubarajan, T. Estimation with Applications to Tracking and Navigation; John Wiley and Sons: Hoboken, NJ, USA, 2001. [Google Scholar]
  16. Blom, H.A.P.; Bar-Shalom, Y.Y. The interacting multiple model algorithm for systems with Markovian switching coefficents. IEEE Trans. Autom. Control 1988, 33, 780–783. [Google Scholar] [CrossRef]
  17. Mahalanobis, P.C. On the generalised distance in statistics. In Proceedings of the National Institute of Sciences, Calcutta, India, 16 April 1936. [Google Scholar]
Figure 1. Schematic of indoor localization system using time difference of arrival (TDOA) measurements.
Figure 1. Schematic of indoor localization system using time difference of arrival (TDOA) measurements.
Electronics 10 00718 g001
Figure 2. State variables for constant-velocity (CV) motion model.
Figure 2. State variables for constant-velocity (CV) motion model.
Electronics 10 00718 g002
Figure 3. Structure of proposed switching extended Kalman filter bank (SEKFB).
Figure 3. Structure of proposed switching extended Kalman filter bank (SEKFB).
Electronics 10 00718 g003
Figure 4. Trajectory of person equipped with mobile tag.
Figure 4. Trajectory of person equipped with mobile tag.
Electronics 10 00718 g004
Figure 5. Root-mean-square position error (RMSPE) in localization simulation: (a) extended Kalman filters (EKFs) using five covariance hypotheses; (b) proposed filter and EKF using best covariance.
Figure 5. Root-mean-square position error (RMSPE) in localization simulation: (a) extended Kalman filters (EKFs) using five covariance hypotheses; (b) proposed filter and EKF using best covariance.
Electronics 10 00718 g005
Figure 6. Root time-averaged mean square error under various measurement-noise conditions.
Figure 6. Root time-averaged mean square error under various measurement-noise conditions.
Electronics 10 00718 g006
Table 1. List of abbreviations.
Table 1. List of abbreviations.
AbbreviationsExplanation
CVConstant velocity
EKFExtended Kalman filter
IMMInteractive multiple model
MCMonte Carlo
RMSPERoot-mean-square position error
RTAMSERoot-time-averaged mean square error
SEKFBSwitching extended Kalman filter bank
TDOATime difference of arrival
TOATime of arrival
TPTransition probability
WSNWireless sensor network
Table 2. Root time-averaged mean square error and number of localization failures.
Table 2. Root time-averaged mean square error and number of localization failures.
FilterRTAMSENumber of Localization Failures
SEKFB0.02490
EKF ( q = 100 )0.02606
EKF ( q = 10 )0.02521
EKF ( q = 1 )0.06260
EKF ( q = 0.1 )0.21990
EKF ( q = 0.01 )0.70450
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pak, J.M. Switching Extended Kalman Filter Bank for Indoor Localization Using Wireless Sensor Networks. Electronics 2021, 10, 718. https://doi.org/10.3390/electronics10060718

AMA Style

Pak JM. Switching Extended Kalman Filter Bank for Indoor Localization Using Wireless Sensor Networks. Electronics. 2021; 10(6):718. https://doi.org/10.3390/electronics10060718

Chicago/Turabian Style

Pak, Jung Min. 2021. "Switching Extended Kalman Filter Bank for Indoor Localization Using Wireless Sensor Networks" Electronics 10, no. 6: 718. https://doi.org/10.3390/electronics10060718

APA Style

Pak, J. M. (2021). Switching Extended Kalman Filter Bank for Indoor Localization Using Wireless Sensor Networks. Electronics, 10(6), 718. https://doi.org/10.3390/electronics10060718

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