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
plays a critical role, but it is a highly uncertain design parameter [
12,
13,
14,
15]. The conventional EKF algorithm uses constant
values that are selected on the basis of an engineer’s knowledge and experience. However, inappropriately selected
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
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
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
and designing the TP. Without the need for the careful selection of
, 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:
where
c is the speed of light;
is the distance between transmitter and
i-th receiver;
are the 2D positions of the transmitter at time
k; and
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:
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
and velocities
as shown in
Figure 2.
State variables at time
k are represented as
and
. To simplify the notation of velocities, we used
instead of
. Thus, state vector
was constructed as
. Thus, the CV motion model [
14,
15] is represented as
where
T is the sampling interval, and
is the process-noise vector. We assumed that
was zero-mean white Gaussian noise with covariance matrix
, where
is a
identity matrix.
The measurement update process of the EKF uses a measurement model. The TDOA measurements are expressed as
where
, and
are TDOA measurements (in units of nanoseconds);
, and
are nonlinear measurement equations that produce TDOA measurements. The measurement vector and measurement equation vector were constructed as
and
, respectively. The measurement model can then be expressed as
where
is the Gaussian measurement noise vector with covariance
.
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:
where
is the estimation error covariance matrix;
is the Kalman gain; and superscripts − and + represent a priori and a posteriori, respectively.
is the Jacobian matrix defined as
and obtained as follows:
Process-noise covariance
in (8) had a significant effect on EKF performance.
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
is difficult. If an unsuitable
value is used, EKF accuracy degrades. Thus, we propose the SEKFB algorithm, which can provide accurate and reliable localization without selecting a suitable
.
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 hypotheses, which involves selecting q hypotheses because we assumed that . We constructed set of hypotheses , 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., ), 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
are obtained. The best state estimate is selected by minimizing Mahalanobis distance between actual TDOA measurement
and predicted measurements
, where
. The predicted measurement is computed using the state estimate as follows:
where
is the state estimate obtained from the
j-th EKF. Mahalanobis distance [
17] is computed as
If the
j-th EKF produces the minimal
,
is chosen as the output of the SEKFB at time
k. For the estimation of the next time step
, 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. |
|
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
m. Four receivers were installed at positions
,
,
, and
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
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
. We assumed that measurement-noise covariance was exactly known. However, process-noise covariance
was unknown. Thus, we roughly selected
q hypotheses as
. 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
where
is the final time step,
is the number of effective simulations, and
, where
is the number of localization failures. A case that produced the RTAMSE at a value greater than
was considered to be a localization failure.
Figure 5a shows the RMSPEs of the EKFs using the five
q hypotheses.
produced the largest RMSPE, that is, the worst performance.
and
produced much smaller RMSPEs compared with the other
q hypotheses. Because RMSPEs produced by
and
were similar and difficult to distinguish, we compared their RTAMSEs.
Table 2 compares the RTAMSEs of the filters;
and
produced RTAMSE values of
and
, respectively. Thus, the best EKF performance could be achieved with
.
In
Figure 5b, we compared the SEKFB and EKF with
. 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
, which was smaller than that of the EKF with
. 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
or
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
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
, where
,
,
,
, and
. 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,
and
, 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
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.