1. Introduction
The Kalman filter and smoother [
1,
2,
3,
4,
5] have been the most commonly fundamental tools for filtering and smoothing in statistical time series analysis. Thus, the Kalman filter and smoother have been a standard choice and a beautiful reference for the state estimation and applied successfully for diverse engineering problems. However, due to their recursive formulations and infinite memory structure (IMS), the Kalman filter and smoother may exhibit performance degradation and even divergence in severe cases for mismodeling and temporary uncertainties.
Therefore, as an alternative to the Kalman filter, the finite memory structure (FMS) filter has been designed for state estimation and shown inherently to have BIBO stability and be more robust against temporary uncertainties [
6,
7,
8,
9]. This FMS filter has been applied successfully in various engineering fields [
10,
11,
12]. Likewise the recursive Kalman filter in [
1,
2,
3,
4,
5], the FMS filter in [
6,
7,
8,
9] is a causal filter that provides state estimates at given times based only on the relative past. Hence, as a noncausal filter, the FMS smoother has been also designed for estimation problems when there is a fixed delay
d between the original state and the estimated state [
13,
14,
15,
16]. The FMS smoother has been shown to be much less computationally complex and more robust against temporary uncertainties than the recursive Kalman smoother. Of course, noncausal filters such as the recursive Kalman smoother and the FMS smoother naturally yield more accurate estimates when there is a fixed delay
d between the original state and the estimated state. However, the FMS smoother could not be better than the FMS filter when there is no fixed delay between the original state and the estimated state. That is, the FMS filter should be applied when there is no fixed delay between signal generation and signal estimation.
Meanwhile, a target tracking problem in indoor positioning systems for wireless network environments has been an interesting research problem and thus applied successfully for locating and tracking people within a building and products stored within a warehouse using wireless sensors and mobile devices [
17,
18,
19]. Due to a variety of indoor wireless channel characteristics, an accurate tracking estimation of a moving target is required for indoor positioning systems. However, a couple of noises such as system noise and measurement noise can contaminate the estimated position and thus can cause deterioration of target tracking performance. Hence, these noises must be filtered for tracking a true path of a moving target. To provide the most accurate position and velocity of a moving target on the real-time indoor positioning in wireless network environments, some research adopted the Kalman filter to reduce the effect of noises [
20,
21,
22]. Actually, past measurement data might have a little information about the target’s current motion since the moving target changes its motion dynamics. In other words, the valid duration of the moving target’s dynamics might be limited to the recently finite time interval. In addition, estimation filters are typically designed with the state-space model called the constant velocity motion model. The constant velocity motion model assumes that targets move with constant velocity, that is, zero-acceleration within a short sampling time. However, in real situations, moving targets maneuver and change velocity and thus move temporarily with nonzero acceleration. Although this can be a temporary uncertainty and thus effects typically occur over a short time interval, the state estimation filter should be robust to diminish the effects of the temporary uncertainty. Therefore, the FMS filter and smoother might be very appropriate for target tracking approaches. The FMS filter was applied successfully for the target tracking in wireless network environments [
23,
24] while the FMS smoother has not been addressed.
In this paper, a state estimation problem is considered for a target tracking scheme in wireless communication environments. Firstly, an alternative FMS filtering and smoothing algorithm is derived for both cases with delay and without delay between the original state and the estimated state. This algorithm is thus called the unified algorithm of FMS filter and smoother. As shown in the terminology unified which can also mean united, combined, integrated, both FMS filter and smoother are derived by solving one optimization problem directly with incorporation of the unbiasedness constraint using only finite measurements and inputs on the most recent window . Thus, the proposed unified algorithm provides simultaneously state estimates at the current time i as well as at the lagged time , given finite measurements and inputs on the most recent window where is the window initial time and i is the current time. The proposed unified algorithm of FMS filtering and smoothing shows that there are some unique properties such as unbiasedness, deadbeat, time-invariance and intrinsic robustness while the recursive IMS filter such as Kalman filter does not have these properties. The on-line computational complexity of the proposed unified algorithm is discussed and compared with the IMS filter. Secondly, a target tracking in wireless network environments is considered as an application of the proposed unified algorithm of FMS filtering and smoothing. Through extensive computer simulations for moving target’s accelerations of various shapes, the proposed unified algorithm-based target tracking scheme is shown to provide estimates for position as well as acceleration of moving target in real time, while eliminating unwanted noise effects and maintaining desired moving positions. In particular, due to intrinsic robustness and deadbeat properties, the performance of the proposed unified algorithm-based scheme can be shown to be better than the existing IMS filtering-based scheme for suddenly changing acceleration.
This paper has the following structure. In
Section 2, a unified algorithm of FMS filtering and smoothing is proposed. In
Section 3, a target tracking scheme using the proposed unified algorithm is considered via extensive computer simulations. Then, concluding remarks are given in
Section 4.
2. Unified Algorithm of Finite Memory Structure Filter and Smoother
Consider the following linear discrete-time state-space model:
where
is the system state variable,
is the control input,
is the sensor measurement output. A couple of noises, the system noise
and the measurement noise
, are zero-mean white Gaussian. These noises are mutually uncorrelated and their covariances are denoted by positive definite matrices
Q and
R, respectively.
At the current time
i, finite measurements
and inputs
on the most recent window
can be expressed by the regression form with the current state
as follows:
where
and matrices
,
,
are defined by
Meanwhile, at the lagged time
, the lagged state
can be represented by
where matrices
and
is defined by
From (
2) and (
5), the current state
and the lagged state
can be represented by the regression form with the augmented state
as follows:
where
The noise term
in (
6) is zero-mean white Gaussian and its covariance given by the positive definite matrix
as follows:
where
Based on the approach of
best linear unbiased estimation in [
25], the unified algorithm of FMS filter and smoother is derived using the following series of equations. The unified algorithm of FMS filter and smoother
is assumed to be obtained from a matrix form using only finite measurements
and inputs
on the most recent window
as follows:
where
is the unified gain matrix. When taking the expectation both sides of (
9) as follows:
is unbiased, i.e.,
, with the following constraint:
Thus, the constraint (
11) can be called the
unbiasedness constraint for the unified algorithm of FMS filter and smoother
.
Subject to the unbiasedness constraint (
11), the objective is now to obtain the gain matrix
in order that the error of
has a minimum variance as follows:
Then, by solving the optimization problem directly with incorporation of the unbiasedness constraint, the unified algorithm of FMS filter and smoother
is obtained by the solution of (
12) as follows:
where
The proposed unified algorithm (
13) of FMS filtering and smoothing shows that there are some unique properties such as unbiasedness, deadbeat, time-invariance and intrinsic robustness. The gain matrix
in (
14) of the unified algorithm requires only once computation on the interval
and then is used for all windows, which means the time-invariance property. As shown in (
10) and (
11), the unified algorithm has the unbiasedness property
. In addition, the proposed unified algorithm
in (
13) for both the current state estimate
and the lagged state estimate
has the deadbeat property as following theorem.
Theorem 1. When , the proposed unified algorithm of FMS filter and smoother on the most recent window is exact for noise-free systems.
Proof of Theorem 1. When there are no noises on the most recent window
for the discrete-time state-space model (
1),
(
6) is represented by
Hence, the following is satisfied:
This completes the proof of the deadbeat property. □
The deadbeat property means that the unified algorithm of FMS filter and smoother
in (
13) tracks exactly the current state
and the lagged state
at every time for noise-free systems although the proposed unified algorithm has been designed assuming a couple noises such as
and
in the discrete-time state-space model (
1). Since the deadbeat property indicates finite convergence time and fast tracking ability, of the proposed unified algorithm of FMS filter and smoother. Thus, it can be expected that the proposed unified algorithm of FMS filter and smoother might be appropriate for a variety of applications requiring fast tracking capability such as fault detection and diagnosis for dynamic process systems as well as maneuver detection and target tracking for indoor positioning systems. Moreover, in contrast to the recursive IMS filter such as the Kalman filter, the proposed unified algorithm can have intrinsic robustness due to its finite memory structure. This means that the proposed unified algorithm can be more robust against round-off errors, mismodeling and temporary uncertainties.
The on-line computational complexity of the proposed unified algorithm is discussed. To simply compare the computational complexity of the proposed unified algorithm and the IMS filter, the measurement
is assumed as a scalar-valued (i.e.,
). In the case of the IMS filter such as the fixed-lag Kalman smoothing filter, the coefficient values for fixed-lag Kalman filtering and fixed-lag smoothing algorithms must be always computed for
d iterations by on-line computing before inputting any data to the fixed-lag Kalman smoothing filter algorithm stage. Thus, the IMS filter such as the fixed-lag Kalman smoothing filter with
nth dimension can be done with
flops as shown in [
26]. On the other hand, as mentioned before, the gain matrix (
14) of the unified algorithm requires only once computation on the interval
and then is used for all windows. Thus, FMS filter and smoother with
nth dimension in the proposed unified algorithm can be done respectively with
flops as shown in
Table 1. Thus, as the dimension grows, the on-line computational complexity of the proposed unified algorithm is shown to be much less than that of the fixed-lag Kalman smoothing filter. However, the proposed unified algorithm requires the memory with size of
M to store finite measurements on the most recent window
and the memory shift with
operations. On the other hand, the fixed-lag Kalman smoothing filter requires to store only intermediate estimates with size of
d. Since
, the fixed-lag Kalman smoothing filter can be better than the proposed unified algorithm in terms of the memory management.
3. Application for Target Tracking Problem
As an application of the proposed unified algorithm of FMS filtering and smoothing, a target tracking problem is considered through extensive computer simulations. In these days, a target tracking problem in indoor positioning systems for wireless network environments has been an interesting research problem and thus been applied successfully for locating and tracking people within a building and products stored within a warehouse using wireless sensors and mobile devices. Due to a variety of indoor wireless channel characteristics, an accurate tracking estimation of a moving target is required for indoor positioning systems. However, a couple of noises such as system noise and measurement noise can contaminate the estimated position and thus can cause deterioration of target tracking performance. These noises can be caused by various factors such as environmental interference, turbulence affecting the target’s movement, inaccuracies of sensor measurement data, and human’s navigation error in a perfect straight line. Hence, these noises must be filtered for tracking a true path of a moving target. To provide the most accurate position and velocity of a moving target on the real-time indoor positioning in wireless network environments, some research adopted both IMS and FMS filters to reduce the effect of noises [
20].
This paper considers only the one dimensional moving target for the
X direction for the simplicity. The state estimation filter for the target tracking is typically designed with the state-space model called the constant velocity motion model. The constant velocity motion model assumes that targets move with constant velocity, that is, zero-acceleration within a short sampling time. However, in real situations, moving targets maneuver and change velocity and thus move temporarily with nonzero acceleration. Thus, the moving target’s successive locations for the
X direction can be represented by the 3rd-order discrete-time state-space model with sampling time
T and acceleration term:
where
The state vector
at time
i consists of three state variables,
to represent the moving target’s random
position,
to represent the corresponding
velocity,
to represent the corresponding
accelerations. Typically, since the acceleration occurs over a short time interval, this can be considered as a temporary uncertainty, that is, an unknown input and thus considered as an unknown state term. Hence, since the target tracking problem in this paper does not need the control input term, the proposed unified algorithm (
13) is applied with assuming that the control input matrix is zero, that is,
. Because of random disturbances by fading and shadowing, the state vector
cannot be measured directly. Thus, the measurement vector
is modeled by only the position to take these effects into account filtering or smoothing estimates. As shown in [
23], target tracking applications using the constant velocity motion model often adopt the state-space model where the state equation is linear and only the measurement equation is nonlinear. Since little research is going on in the case of the FMS smoother for nonlinear systems, the proposed unified algorithm is currently difficult to use this type of state space model. Thus, this paper consider only the linear measurement equation. The application of the nonlinear measurement equation can be a future work of the current paper.
Using (
13) with (
15) and (
16), the proposed unified algorithm of a FMS filtering and smoothing-based target tracking scheme provides estimates for position as well as acceleration of moving target in real time, while eliminating unwanted noise effects and maintaining desired moving positions. Via extensive computer simulations, the performance of the proposed unified algorithm-based target tracking scheme is evaluated and compared with the existing IMS filtering-based scheme of [
20,
21,
22]. Computer simulations are performed for moving target’s accelerations of four kinds of shapes. The first and second scenarios consider that the target moves with suddenly changing acceleration such as step-type and ramp-type. The third case considers that the target moves with slowly changing acceleration such as triangle-type. The last case considers that the target moves with randomly changing acceleration. The sampling period is taken by
. The covariance of the system noise is
and the covariance of measurement noise is
. The window length is taken by
. The lagged length is taken by
. Simulations of 30 runs are performed using different system and measurement noises to make the comparison clearer. Each single simulation run lasts 800.
Figure 1,
Figure 2,
Figure 3 and
Figure 4 show simulation results for moving target’s accelerations of four kinds of shapes. The first plot of figures shows the root-mean-square(RMS) estimation error of the moving target’s random position for 30 simulations. The second plot of all figures shows the estimation error of the moving target’s random position for one of 30 simulations. The last plot of figures shows the estimate of unknown acceleration. In addition, time averaged values of RMS estimation errors are presented by
Table 2. As shown in simulation results, the proposed unified algorithm-based target tracking scheme can outperform the IMS filtering-based target tracking scheme when the target moves with suddenly changing acceleration. On the interval where the acceleration varies suddenly, the estimation error of the proposed unified algorithm-based scheme is remarkably smaller than that of the IMS filtering-based scheme. In addition, when the acceleration varies constantly, the convergence of the estimation error for the proposed unified algorithm-based scheme is much faster than that of the IMS filtering-based scheme. These observations for simulation results might come from the fast convergent time and the fast tracking ability due to intrinsic robustness and deadbeat properties of the proposed unified algorithm-based scheme.
4. Conclusions
This paper has dealt with the state estimation problem for the target tracking scheme in wireless communication environments. Firstly, the unified algorithm of FMS filtering and smoothing has been proposed. The unified algorithm has been derived by solving one optimization problem directly with incorporation of the unbiasedness constraint using only finite measurements and inputs on the most recent window. The proposed unified algorithm provides simultaneously the current state estimate as well as the lagged state estimate. The proposed unified algorithm has shown that there are some unique properties such as unbiasedness, deadbeat, time-invariance and intrinsic robustness while the recursive IMS filter such as Kalman filter does not have these properties. The on-line computational complexity of the proposed unified algorithm has been compared with the IMS filter such as the fixed-lag Kalman smoothing filter. Secondly, the target tracking in wireless network environment has been considered via extensive computer simulations as an application of the proposed unified algorithm of FMS filtering and smoothing. Through simulation results for moving targets’ accelerations of various shapes, the proposed unified algorithm-based target tracking scheme has been shown to provide estimates for position as well as acceleration of moving target in real time, while eliminating unwanted noise effects and maintaining desired moving positions. It has been shown that the performance of the proposed unified algorithm-based target tracking scheme can be better than the IMS filtering-based scheme for suddenly changing acceleration.
Actually, research on the FMS filter and smoother in nonlinear systems is relatively inactive compared to research in linear systems. There is some research, such as an adoption of nonlinear measurement model, on the FMS filter for nonlinear systems, but little research is going on in the case of the FMS smoother for nonlinear systems. Thus, for nonlinear systems, research on the FMS smoother should be preceded and then a unified algorithm of FMS filter and smoother can be researched as a future work. In addition, the estimation of target’s acceleration can be extended to identification and classification of moving targets, which could be another future work.