Next Article in Journal
Vibro-Shock Dynamics Analysis of a Tandem Low Frequency Resonator—High Frequency Piezoelectric Energy Harvester
Next Article in Special Issue
Flexible Fusion Structure-Based Performance Optimization Learning for Multisensor Target Tracking
Previous Article in Journal
Distributed Data Service for Data Management in Internet of Things Middleware
Previous Article in Special Issue
State Estimation Using Dependent Evidence Fusion: Application to Acoustic Resonance-Based Liquid Level Measurement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Auxiliary Truncated Unscented Kalman Filtering for Bearings-Only Maneuvering Target Tracking

Automatic Target Recognition Key Laboratory (ATR), Shenzhen University, Shenzhen 518060, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(5), 972; https://doi.org/10.3390/s17050972
Submission received: 13 March 2017 / Revised: 24 April 2017 / Accepted: 25 April 2017 / Published: 27 April 2017

Abstract

:
Novel auxiliary truncated unscented Kalman filtering (ATUKF) is proposed for bearings-only maneuvering target tracking in this paper. In the proposed algorithm, to deal with arbitrary changes in motion models, a modified prior probability density function (PDF) is derived based on some auxiliary target characteristics and current measurements. Then, the modified prior PDF is approximated as a Gaussian density by using the statistical linear regression (SLR) to estimate the mean and covariance. In order to track bearings-only maneuvering target, the posterior PDF is jointly estimated based on the prior probability density function and the modified prior probability density function, and a practical algorithm is developed. Finally, compared with other nonlinear filtering approaches, the experimental results of the proposed algorithm show a significant improvement for both the univariate nonstationary growth model (UNGM) case and bearings-only target tracking case.

1. Introduction

Bearings-only maneuvering target tracking has been widely researched for decades. It is important for many applications such as maritime surveillance, navigation and aerospace, wireless sensor networks (WSN), and infrared search and track (IRST) systems [1,2,3,4,5,6]. However, while implementing this technology in unlimited situations, there remain some challenging problems, such as multiple platform tracking, uncertainty of the target model and nonlinear non-Gaussian noise. To deal with the uncertainty of the motion model, such as abrupt target maneuver, heavy clutter measurements, highly nonlinearity of dynamic models and nonlinear non-Gaussian noise, etc., the interacting multiple model (IMM) [7] based on the nonlinear filtering algorithm is a promising approach. However, to model the uncertainty of the motion model, the performance of the IMM-type algorithm is directly proportional to the number of the motion models. Generally, the more motion models we produce, the greater accuracy of the estimated state we obtain. However, the computational complexity of the algorithm becomes larger with the increase of the numbers of motion models, particularly in heavily cluttered environments. Moreover, the nonlinear filtering has been studied extensively in bearings-only maneuvering target tracking.
As is well-known, the most widely used nonlinear filtering for bearing-only tracking is to employ an extended Kalman filter (EKF) [8,9]. However, when the nonlinearity of dynamic models becomes more severe, the performance of the EKF degrades sharply. In order to solve this problem, the unscented Kalman filter (UKF) [10] and the truncated unscented Kalman filtering (TUKF) were proposed [11,12]. Compared with other conventional Kalman filter-type approaches, the TUKF can achieve better performance in the conditions of the target tracking system, and can provide very informative nonlinear measurements compared to the prior. Moreover, to take into account the available additional information the state given by the constraint, Ondrej et al. [13] proposed a generic local filter for the inequality constrained estimation problem, and designed an efficient truncation technique based on the Monte Carlo integration method for the approximation of the state probability density function. Beatriz et al. [14] proposed a constrained dual state and parameter estimation algorithm using a dual Kalman filter (DKF) and a probability density function (PDF) truncation algorithm for analysis of lateral vehicle dynamics.
In recent years, particle filtering has been widely used for bearing-only tracking. In [15], Gordon proposed the first particle filtering algorithm based on the resampling step. The main idea is that the posterior distribution can be approximated by series of random samples with associated weights, and its parameter estimates can be computed by these samples and weights. Therefore, particle filtering can deal with nonlinear non-Gaussian problems in terms of the dynamics and measurements. Recently, many particle filtering methods have been proposed [16,17,18], for example, the extended Kalman particle filter (EKF-PF), unscented particle filter (UPF), and the multivariable feedback particle filter (GPF) [18]. Moreover, for the maneuvering target tracking problem, many particle filters have been proposed based on Markovian switching systems [19,20,21,22]. Boers et al. [19] proposed a interacting multiple model particle filter algorithm (IMM-PF) by combining a mixture of the interacting multiple model (IMM) filter with the particle filter. For the maneuvering target tracking problem in bearings-only wireless sensor networks (WSNs), Atiyeh et al. [20] proposed a interacting multiple model particle filter to estimate the state variables of the moving target. Li et al. [21] proposed a Rao–Blackwellized particle filter based on multiple model algorithm for maneuvering target tracking in a cluttered environment. Yu et al. [22] proposed a distributed particle filter by incorporating the curvature of the sensing region in the measurement model for bearings-only tracking of a moving target. In their method, to reduce the communication load, the transformation of the observations is approximated as Gaussian distribution, which the variance can be approximated using the average variance over all particles. However, abrupt target maneuvers, modeling uncertainty and the high nonlinearity of model function remain to be unsolved issues.
For achieving a successful tracking performance, the aforementioned methods require accurate motion models and adaptive nonlinear filtering methods. However, particularly in maneuvering target tracking, accurate motion modeling is almost impossible, and an adaptive nonlinear filtering needs to be used to handle abrupt maneuver of target. More importantly, these two challenges are not separate problems and should be considered simultaneously. In previous research [23,24], Ehsan et al. [23] proposed a new bearing-only bias estimation model based on triangulation using the associated measurement reports (AMR) or local bearing-only tracks from different sensor pairs for distributed tracking systems. Li et al. [24] proposed novel truncated quadrature Kalman filtering (TQKF) based on the Gauss-Hermite quadrature rule for bearings-only maneuvering target tracking. In order to avoid the requirement of the measurement function being bijective, the modified prior PDF of the TQKF algorithm can be approximately computed by the least square estimation approach. However, the most important limitation of the TQKF is the expensive computational burden, and it cannot be used for real time target tracking.
In this paper, novel auxiliary truncated Kalman filtering (ATUKF) is proposed for bearings-only maneuvering target tracking. Unlike the TUKF algorithm, to overcome the modeling uncertainty, a modified prior probability density function (PDF) is defined based on several auxiliary target characteristics and current measurements, which can effectively minimize the variance of the prior distribution. Moreover, to achieve the requirement of bijective measurement function, the statistical linear regression based on the unscented transformation is used to linearize the nonlinear measurement function, and the modified prior PDF is approximated as Gaussian. Finally, the posterior PDF can be approximately estimated based on the prior PDF and the modified prior PDF, and a practical algorithm is developed for bearings-only target tracking systems.
The rest of the paper is organized as follows. The proposed algorithm is given in Section 2. In Section 3, we provide the experimental results. Finally, some conclusions are given in Section 4.

2. Proposed Algorithm

In order to track the maneuvering target, accurate motion modeling and nonlinear filtering are two challenging problems that should not be separated. However, most research on maneuvering tracking investigates these problems separately. In this section, in the bearings-only maneuvering target tracking, novel auxiliary truncated unscented Kalman filtering is proposed. In Section 2.1, the joint prior distribution is approximately constructed. In Section 2.2, in order to track the bearings-only maneuvering target, the modified prior PDF is approximated based on statistical linear regression by introducing the target spatio-temporal information. Section 2.3 summarizes the proposed algorithm.

2.1. Joint Prior Distribution

Suppose the target dynamic system can be written as:
x n = f ( x n 1 ) + m n
z n = h ( x n ) + e n
where z n R n z denotes the observation vector at time n , x n R n x denotes the target state vector at time n , and f ( ) and h ( ) denote the nonlinear state transition function and observation function, respectively.
Suppose that r n denotes the set of target characteristics including c independent components r n = { r n 1 , r n 2 , , r n c } . In order to derive the proposed algorithm, there are two basic hypotheses, firstly, that the nonlinear function h n ( ) in (2) is a bijective, continuous function; and secondly, that the probability density function of the measurement noise e n has bounded, connected support.
p e n ( v n ) = 0 , e n I e n n z
where I e n is an n z -dimensional measurement validation region. Therefore, according to the second assumption, the measurement likelihood function can be defined as follows:
p ( z n | x n , r n ) = p e n ( z n h n ( x n ) ) χ I e n ( z n h n ( x n ) )
p ( z n | x n , r n ) = p e n ( z n h n ( x n ) ) χ I x n ( z n ) ( x n )
I x n ( z n ) = { x n | x n = h n 1 ( z n e n ) , e n I e n }
where χ I e n ( ) is the indicator function on the subset I e n . Therefore, the state posterior PDF can be defined as:
p ( x 0 : n | z 1 : n , r 1 : n ) = p ( z n | x 0 : n , z 1 : n 1 , r 1 : n ) p ( x n | x 1 : n 1 , z 1 : n 1 , r 1 : n ) p ( x 1 : n 1 | z 1 : n 1 , r 1 : n 1 ) p ( z n , r n | z 1 : n 1 , r 1 : n 1 ) = p e k ( z n h ( x n ) ) χ I x ( z n ) ( x n ) p ( x n | x n 1 , z 1 : n 1 , r 1 : n ) p ( x 1 : n 1 | z 1 : n 1 , r 1 : n 1 ) p ( z k , r k | z 1 : k 1 , r 1 : k 1 ) p ( z n | x n ) p 1 ( x n | z n , x n 1 , r 1 : n ) p ( x 1 : n 1 | z 1 : n 1 , r 1 : n 1 )
p 1 ( x n | z n , x n 1 , r 1 : n ) = p ( x n | x n 1 , z 1 : n 1 , r 1 : n ) χ I x n ( z n ) ( x n ) / ε 1
where ε 1 is a constant. From Equation (8), we can see that the modified prior PDF is defined by incorporating the current measurement information z n . According to the conclusions in [11], if the measurement noise is informative, the modified prior p 1 ( ) is not only the minimum variance of the prior p 0 ( ) , but also can improve the algorithm’s performance. Further, to deal with the uncertainty of motion models, the joint prior distribution p ( x n | x 0 : n 1 , z 1 : n , r 1 : n ) of the proposed algorithm can be defined as follows:
p ( x n | x 0 : n 1 , z 1 : n , r 1 : n ) = α n p 1 ( x n | z n , x n 1 , r 1 : n ) + ( 1 α n ) p 0 ( x n | x n 1 , z 1 : n 1 ) = α n p ( x n | x n 1 , r 1 : n ) χ I x n ( z n ) ( x n ) + ( 1 α n ) p 0 ( x n | x n 1 , z 1 : n 1 )
where α n [ 0 , 1 ] is a proper parameter. To approximately calculate the mean and covariance of the posterior distribution, we apply a UKF update to p 0 ( ) (UKF1), and another UKF update to p 1 ( ) (UKF2). Finally, the posterior estimates can be approximately calculated through merging both results obtained by UKF1 and UKF2.

2.2. Approximation of p 1 ( )

In the subsection, our object is to approximate the modified prior PDF p 1 ( ) as Gaussian. For this reason, we write the state vector as x n = [ a n T , b n T ] T , where a n n a denotes the position components of x n , b n n b denotes the velocity components of x n , and n x = n a + n b . The derivation of the mean x ^ p 1 , n and covariance P p 1 , n of p 1 ( ) is the same as in the truncated unscented Kalman filter (TUKF) [11], which can be shown as follows:
x ^ p 1 , n = [ μ a n , 1 μ b n , 1 ] , P p 1 , n = [ a n , 1     a n b n , 1 ( a n b n , 1 ) T    b n , 1 ]
Σ a n , 1 = H ˜ n 1 R n ( H ˜ n 1 ) T
where μ b n , 1 , Σ b n , 1 , Σ a n b n , 1 can be found in [11]. μ a n , 1 denotes the estimated mean of state a n , R n denotes the measurement noise covariance, and H ˜ n 1 = [ a n h n T ( a n ) ] T | a n = μ a n , 1 is the Jacobian of h n T ( a n ) evaluated at μ a n , 1 .
Now, how to calculate the estimated mean μ a n , 1 remains a key problem to be solved. For the passive sensor tracking system, the modeling of target dynamic system is a challenging problem when the target maneuvers, and some auxiliary target characteristics need to be used to deal with arbitrary changes in motion models. To achieve a high tracking performance, a statistical linear regression method (SLR) [25] is proposed to estimate the state mean μ a n , 1 .
Firstly, to evaluate the state mean μ a n , 1 , three approximations are used: (S1) the prior PDF p 0 ( a n ) is constant over the connected region I a n ( z n ) ; (S2) the nonlinear function h ( ) can be locally linearized; and (S3) the measurement noise satisfies uniform distribution e n ~ U I e n in the connected region I e n . According to S2, the nonlinear measurement Equation (2) can be approximated as a linear estimator of z n , z ^ n such that:
z ^ n = n a n + d n
where n denotes a linear measurement matrix, and d n denotes a noise vector, which are derived by minimizing the objective function defined as follows:
{ n , d n } = arg min E ( τ n T τ n )
where τ n is the linearization error, τ n = z n z ^ n .
Substituting τ n into (13), and setting the partial derivative of the objective function with d n to zero,
( 2 ) E ( z n n a n d n ) = 0 d n = z ¯ n n a ¯ n
where a ¯ n = E ( x n ) and z ¯ n = E ( z n ) . Substituting d n into (13)
τ T τ = [ ( z n z ¯ n ) n ( a n a ¯ n ) ] T [ ( z n z ¯ n ) n ( x n x ¯ n ) ]
Then, setting the gradient with respect to n to zero,
( 2 ) E { [ ( z n z ¯ n ) n ( a n a ¯ n ) ] [ z n z ¯ n ] ) } = 0
Solving (15) for n , we obtain:
n = P a n z n T P a n a n 1
where P a n a n = E [ ( a n a ¯ n ) ( a n a ¯ n ) T ] , P a n z n = E [ ( a n a ¯ n ) ( z n z ¯ n ) T ] . According to the measurement equation (12), the maximum likelihood position a ^ n ( z n ) of target state a n can be estimated as follows:
a ^ ( z n ) = n ( z n d n )
If the maximum likelihood estimate a ^ n ( z n ) is used to replace the estimated mean μ a n , 1 , the performance of the proposed algorithm will be improved because the current measurement information is incorporated. However, it cannot solve the uncertainty in motion models. In particular, when the target speed or the measurement sampling interval is large, the tracking performance degrades.
More recently, sophisticated techniques have been based on the target motion characteristics [6], which have been proposed to address the challenges in motion modeling. In their proposed method, three target characteristics, such as actual target speed v , time interval T of measurement and course angle θ of the target, are considered to improve the tracking performance. In (18), the relationship between three target characteristics and the state predicted error is given:
σ = ( T v ) 2 + ( T v n ) 2 2 T 2 v v n cos ( θ ) = T ( v ) 2 + ( v n ) 2 2 v v n cos ( θ )
where v n denotes the current estimated velocity, θ denotes the estimated error of course angle, and v denotes the actual target velocity. Supposing v = v v n , we can obtain
σ = T ( v ) 2 + ( v n ) 2 2 v v n cos ( θ ) = T v n 2 + 2 v ( v v n ) ( 1 cos ( θ ) )
From Equation (19), we can find that the predicted error σ increases monotonically with the increase of the parameters ( v n , T , θ ). In fact, when v n > 50 m / s and T > 20 s , the predicted error σ will larger than 1000 m. It shows that the predicted error becomes a major reason for the performance degradation. On the other hand, when the actual target speed is relatively small, the prediction error caused by T or θ is smaller than the measurement error. Therefore, to improve the performance of mean estimate μ a n , 1 , the maximum likelihood estimate a ^ n ( z n ) is considered as the latest observation, and the modified maximum likelihood estimates that can incorporate current target characteristics is defined as follows:
ϕ ^ ( z n ) = μ a n , 0 + K n ( a ^ ( z n ) n μ a n , 0 )
K n = ( T 2 v 2 σ v 2 ( n ) ) / ( λ σ e 2 ( k ) + T 2 v 2 σ v 2 ( n ) )
where λ is a constant, σ e 2 ( n ) denotes the variance of measurement noise, and σ v 2 ( n ) denotes the innovation variance. According to (20) and (16), the mean μ a n , 1 and the variance Σ a n , 1 in (10) can be approximated as:
μ a n , 1 = ϕ ^ ( z n ) , Σ a n , 1 = P a n a n
Finally, the modified prior PDF p 1 ( ) can be approximated as a Gaussian probability density function p 1 ( ) N ( x n , x ^ p 1 , n , P p 1 , n ) .

2.3. Summary of the Proposed Algorithm

According to the descriptions above, in order to describe clearly the proposed algorithm, the diagram of the ATUKF is shown in Figure 1. In Figure 1, it is shown that one cycle of the ATUKF algorithm consists of the following steps: (A) time update (predicted by using Kalman filtering); (B) the measurement updates based on the prior p 0 ( ) and the modified prior p 1 ( ) ; and (C) weight calculation and the joint state update. According to the derived results mentioned above and Figure 1, the detailed information of the new ATUKF algorithm can be summarized as follows.

Algorithm: Auxiliary Truncated Unscented Kalman Filtering (ATUKF)

ATUKF—Update based on the prior probability density function (PDF) p 0 ( )
  • Obtain N = 2 n x + 1 sigma points χ 0 1 , χ 0 2 ... , χ 0 N and the corresponding associated weights w 1 , w 2 , ... , w N using unscented transform (UT ) based on the mean x ^ n 1 | n 1 and covariance P n 1 | n 1 of the posterior PDF p 0 ( x n | x n 1 , z 1 : n 1 , r 1 ; n 1 ) , where n x denotes the dimension of state x . The predicted sigma points can be obtained by the nonlinear state function f ( ) :
    χ 0 , n | n 1 i = f ( χ 0 i ) , i = 1 , 2 , , N
  • Approximate the mean and covariance of the state-predicted prior PDF p 0 ( x n | x n 1 , z 1 : n 1 )
    x ^ p , 0 , n | n 1 = i = 1 N w i χ 0 , n | n 1 i
    P p , 0 , n | n 1 = Q n + i = 1 N w i ( χ 0 , n | n 1 i x ^ p , 0 , n | n 1 ) ( χ 0 , n | n 1 i x ^ p , 0 , n | n 1 ) T
  • Compute the predicted measurement z ^ 0 , n | n 1 based on the nonlinear measurement function h ( ) :
    z 0 , n | n 1 i = h ( χ 0 , n | n 1 i ) ,   i = 1 , 2 , , N
    z ^ 0 , n | n 1 = i = 1 N w i z 0 , n | n 1 j
  • The cross-covariance, innovation covariance and error covariance are estimated as follows:
    P x z , 0 , n | n 1 = i = 1 N w i ( χ 0 , n | n 1 i x ^ p , 0 , n | n 1 ) ( z 0 , n | n 1 i z ^ 0 , n | n 1 ) T , ( C r o s s   covariance )
    P z z , 0 , n | n 1 = R n + i = 1 N w i ( z 0 , n | n 1 i z ^ 0 , n | n 1 ) ( z 0 , n | n 1 i z ^ 0 , n | n 1 ) T ,   ( I n n o v a t i o n   covariance )
    P x x , 0 , n | n 1 = i = 1 N w i ( χ 0 , n | n 1 i x ^ p , 0 , n | n 1 ) ( χ 0 , n | n 1 i x ^ p , 0 , n | n 1 ) T ,   ( E r r o r   covariance )
  • Estimate the mean x ^ u , 0 , n | n and covariance P u , 0 , n | n using (30) and (31):
    x ^ u , 0 , n | n = x ^ p , 0 , n | n 1 + P x z , 0 , n | n 1 P z z , 0 , n | n 1 1 ( z n z ^ 0 , n | n 1 )
    P u , 0 , n | n = P p , 0 , n | n 1 P x z , 0 , n | n 1 P z z , 0 , n | 1 1 P x z , 0 , n | n 1 T
ATUKF—Update based on the modified prior PDF p 1 ( )
  • Calculation of the mean x ^ p , 1 , n | n 1 and covariance P p , 1 , n | n 1 of the prior p 1 ( )
    According to (14) and (16) in Section 2.2, the linear regression coefficients n and d n can be approximately computed by using Equations (27–29). The mean x ^ p , 1 , n | n 1 and covariance P p , 1 , n | n 1 of p 1 ( ) can be approximately estimated by (10) and (11), respectively.
  • Draw N new sigma points χ 1 , n | n 1 1 , χ 1 , n | n 1 2 , , χ 1 , n | n 1 N with the associated weights w 1 , w 2 , ...... w N by using the UT based on the mean x ^ p , 1 , n | n 1 and covariance P p , 1 , n | n 1 . The predicted measurements of new sigma points are estimated as follows:
    z 1 , n | n 1 i = h ( χ 1 , n | n 1 i )
  • Calculation of z ^ 1 , n | n 1 , P z z 1 , n | n 1 and P x z 1 , n | n 1
    z ^ 1 , n | n 1 = i = 1 N w i z 1 , n | n 1 j
    P z z , 1 , n | n 1 = R n + i = 1 N w i ( z 1 , n | n 1 i z ^ 1 , n | n 1 ) ( z 1 , n | n 1 i z ^ 1 , n | n 1 ) T ,   ( Innovation   covariance )
    P x z , 1 , n | n 1 = i = 1 N w i ( χ 1 , n | n 1 i x ^ p , 1 , n | n 1 ) ( z 1 , n | n 1 i z ^ 1 , n | n 1 ) T ,   ( Cross   covariance )
  • Estimate the mean x ^ u , 1 , n | n and covariance P u , 1 , n | n using (36) and (37):
    x ^ u , 1 , n | n = x ^ p , 1 , n | n 1 + P x z , 1 , n | n 1 P z z , 1 , n | n 1 1 ( z n z ^ 1 , n | n 1 )
    P u , 1 , n | n = P p , 1 , n | n P x z , 1 , n | n 1 P z z , 1 , n | n 1 1 P x z , 1 , n | n 1 T
ATUKF—Jointly update
  • Calculate the parameter α n using (38) and (39)
    μ 1 ( x ^ u , i , n ) = 1 / | P u , i , n | exp ( ( z n h n ( x ^ u , i , n ) ) / 2 ) , i = 0 , 1
    α n = μ 1 ( x ^ u , 1 , n ) / ( μ 0 ( x ^ u , 0 , n ) + μ 1 ( x ^ u , 1 , n ) )
  • Approximate the mean x ^ n and covariance P n of the posterior PDF p ( x n | z 1 : n ) using (23) and (24).
    x ^ n | n = α n x ^ u , 1 , n | n + ( 1 α n ) x ^ u , 0 , n | n
    P n | n = α n [ P u , 1 , n | n + ( x ^ u , 1 , n | n x ^ ) n | n ( x ^ u , 1 , n | n x ^ ) n | n T ] + ( 1 α n ) [ P u , 0 , n | n + ( x ^ u , 0 , n | n x ^ n | n ) ( x ^ u , 0 , n | n x ^ n | n ) T ]

3. Simulation Results

In this section, to evaluate the tracking performance of the ATUKF algorithm, two examples are employed. In Section 3.1, the univariate nonstationary growth model (UNGM) is discussed [11]. In Section 3.2, a bearings-only maneuvering tracking scenario [24], interested in defense applications, is discussed. In the first case, the EKF, UKF, the quadrature Kalman filtering (QKF) [25], the mixture truncated unscented Kalman filter (MTUKF, with three mixture components) [12] and particle filtering(PF) are utilized. In the second case, the TQKF, interacting multiple model extended Kalman filtering(IIMMEKF) and the interacting multiple model Rao–Blackwellized particle filter (IMMRBPF) are employed. In all the experiments, each simulation has been repeatedly performed 100 times.

3.1. Univariate Nonstationary Growth Model (UNGM)

In this section, due to the highly nonlinearity and non-stationarity of dynamic system, the univariate nonstationary growth model is considered. The discrete time system of this model can be written as:
x n = α x n 1 + β x n 1 1 + x n 1 2 + γ cos ( 1.2 ( n 1 ) ) + m n
z n = { φ 2 x n 2 + e n        n 30 φ 1 x n 3 2 + e n      n > 30
where the process noise m n is satisfied with Gaussian distribution with zero mean and variance one, and e n is satisfied with Gaussian distribution with zero mean and variance 0.01 . α = 0.5 , β = 25 , γ = 8 , φ 1 = 0.2 and φ 2 = 0.05 are known constants. In each Monte Carlo simulation, we assume that the initial distribution of state x 0 is uniform distribution in the interval [0 1]. The number of particles is 1000.
Figure 2 shows the root-mean-square position errors (RMSE) of the EKF, UKF, QKF, PF, MTUKF and ATUKF. It is obvious from Figure 2 that the performance of the ATUKF is much better than that of the EKF, UKF and QKF. In this case, the performance of the EKF is the poorest. A reason for the poor performance of the EKF, UKF and QKF is the increase of the approximate error arising from the high nonlinearity and the non-stationarity of the dynamic system. Figure 3 shows the RMS position errors of all filters with different noise variance σ e n ~ [ 0.1   5 ] . From Figure 3, it is seen that whenever the measurement is informative ( σ e n < 1 ) or the measurement is uninformative ( σ e n > 1 ), the ATUKF is robust in all situations, its performance is similar to the MTUKF’s, and it is very close to that of the PF. Moreover, among the EKF, UKF and QKF, the EKF has the poorest performance. In particular, when the measurement is very informative ( σ e n < 1 ), the EKF yields a divergent estimate.
Table 1 shows the computation time statistics for all algorithms. In this case, all the experiments are performed by using MATLAB programming on Intel-Core(TM)-i2-4030U processor (1.9 GHz) based on the Windows platform. It can be seen from Table 1 that the computational load of the PF is the largest than these of other filters, such as the EKF, UKF, QKF, MTUKF and ATUKF. The ATUKF is very close to the QKF. Furthermore, the computation time for the ATUKF is much lower than the MTUKF. The main reason is that the MTUKF approximates the posterior PDF as a Gaussian mixture, and it makes the computational burden increase.

3.2. Bearings-Only Maneuvering Tracking (BOT) Scenario

In this scenario, the target makes two circular turns with rectilinear segments connecting them. Figure 4 shows the true target trajectory. The speed modulus is kept constant throughout ( 0.3   km / s ) . The initial position is ( 2   km ,   8   km ,   1   km ) , and the initial velocity is ( 0.15   km / s ,   0.26   km / s ,   0.0   km / s ) . The segments are defined as follows:
First segment. Rectilinear flight until the plane is at ( 6.35   km ,   15.53   km ,   1   km ) (from t = 0 s to t = 30 s).
Second segment. Circular turn with turn rate 6 o / s (from t = 31 s to t = 50 s).
Third segment. Rectilinear flight until the plane is at ( 14.31   km ,   10.33   km ,   1   km ) (from t = 51 s to t = 70 s).
Fourth segment. Circular turn with turn rate 4.8 o / s (from t = 71 s to t = 95 s).
Fifth segment. Rectilinear flight until the plane is at ( 21.26   km ,   11.63   km ,   1   km ) (from t = 96s to t = 100 s).
The motion model of target is defined as follows:
x n s = F s x n 1 s + m n s
where x n s = ( x , y , z , x , y , z ) denotes the target state vector under model x , x , y and z denote the target position coordinates, x , y and z are the target speed in x , y and z directions, respectively, F s denotes the transition matrix, and s { 1 , 2 , , M } denotes the target model index. In the maneuvering target tracking scenario, only a constant velocity model is used for the ATUKF algorithm and TQKF algorithm. In the IMMEKF and IMMRBPF, there are both clockwise- and counterclockwise-coordinated turn models that are used to simulate the target maneuvering. The details of three target motion models are defined as follows:
Model 1: Constant Velocity Motion
The state transition matrix and the process noise covariance matrix are defined by:
F 1 = [ 1 0 0 T 0 0 0 1 0 0 T 0 0 0 1 0 0 T 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 ]
Q 1 = [ 1 4 T 4 0 0 1 2 T 3 0 0 0 1 4 T 4 0 0 1 2 T 3 0 0 0 1 4 T 4 0 0 1 2 T 3 1 2 T 3 0 0 T 2 0 0 0 1 2 T 3 0 0 T 2 0 0 0 1 2 T 3 0 0 T 2 ] σ n 2
where T denotes the sampling interval, in this paper, T is set to 1.
Model 2: Constant Turn Motion
The state transition matrix is:
F 2 = [ 1 0 0 sin ( w ) w cos ( w ) 1 w 0 0 1 0 1 cos ( w ) w sin ( w ) w 0 0 0 1 0 0 1 0 0 0 cos ( w ) sin ( w ) 0 0 0 0 sin ( w ) cos ( w ) 0 0 0 0 0 0 1 ]
where w is a constant angular rate. In this paper, w is set to 0.0175 . The process noise covariance matrix Q 2 is the same as in Model 1.
Model 3: w > 0 describes a clockwise turn, and Model 3 is its natural counterpart for a counterclockwise turn w < 0 .
Two passive sensors are deployed in ( 0 ,   5   km ,   0 ) and ( 0 ,   5   km ,   0 ) respectively. Using the detection fusion architecture [24], the azimuth and elevation angles of aircraft, α i and β i respectively, measured by sensor i , are transmitted to the fusion node. The measurement function is written as:
h ( x n ) = ( α i β i ) = ( arctan ( y s i , y x s i , x ) arctan ( z s i , z ( x s i , x ) 2 + ( y s i , y ) 2 ) )
where ( s i , x , s i , y , s i , z ) , i = 1 , 2 denote the positions of the stationary sensors. The measurement covariance can be defined as R n = [ 1 0 0 1 ] σ e n 2 .
The real initial position of the target is ( 2   km ,   8   km ,   1   km ) , and the initial velocity is ( 0.15   km / s ,   0.26   km / s ,   0.0   km / s ) . The prior PDF of state x 0 is assumed to be x 0 ~ N ( x ^ 0 | 0 , P ^ 0 | 0 ) , where x ^ 0 | 0 = [ 2.1   km   0 0.12   kms 1 7.95   km   0.23   kms 1   0.95   km   0   kms 1 ] T , P ^ 0 | 0 = d i a g [ 0.144   km 2   0.02 2   km 2 s 2   0.144   km 2   0.02. 2 km 2 s 2   0   0 ] T . The standard deviation of the process and the measurement noise are set to σ m n = 0.01 and σ e n = 0.001 , and these noises are zero-mean Gaussian-distributed and independent. The number of particles is set to 200.
Figure 5 shows the X RMSE, Y RMSE, Z RMSE and position RMSE of the ATUKF compared with IMMEKF, IMMRBPF and TQKF. From Figure 5, we can see that the RMSE of all algorithms abruptly increased from 30 s to 50 s, which is mainly due to the increase of the target predicted errors caused by the target maneuvering. However, in Figure 5a,c,d, it is shown that the performance of the ATUKF algorithm has outperformed the IMMEKF, IMMRBPF and TQKF. A key reason is that the proposed algorithm can incorporate the target characteristic information and current measurement information into the prior PDF, which can effectively degrade the variance of errors caused due to the maneuvering of the target. Moreover, because the flight height of target remained unchanged, from Figure 5c, the Z RMSE of all algorithms is very close.
Finally, the computation time statistics for all algorithms are given in Table 2. In this case, all the experiments are performed by using MATLAB programming on an Intel-Core(TM)-i2-4030U processor (1.9 GHz) based on the Windows platform. In Table 2, it is shown that the computational load of the IMMRBPF is much higher than these of the IMMEKF, TQKF and ATUKF. More importantly, the ATUKF requires much less of a computation time than the TQKF. However, the computational load of the ATUKF is nearly two times higher than that of the IMMEKF.

4. Conclusions

In this paper, we presented a bearings-only target tracking algorithm based on an auxiliary truncated unscented Kalman filtering (ATUKF) algorithm. Unlike the truncated unscented Kalman filtering, in the proposed algorithm, several target characteristics were introduced to construct the modified prior PDF, and the statistical linear regression was used to linearize the nonlinear non-bijective measurement function by using the sigma points. Moreover, we have developed a practical algorithm for a bearings-only target tracking system. Finally, in the simulation results, compared with the EKF, UKF, the quadrature Kalman filtering (QKF), the mixture truncated unscented Kalman filter (MTUKF) and the particle filter (PF), the ATUKF exhibits better performance. For the second case, compared with the IMMEKF algorithm, the IMMRBPF algorithm and the TQKF algorithm, the ATUKF algorithm not only improves the performance of the tracker, but significantly reduces the computation time.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant Nos. 61301074, 61271107, 61375015), Natural Science Foundation of the Guangdong Province of China (Grant No. S2012010009417), Defense Advance Research Fund Project (Grant No. 91400***501140C80340).

Author Contributions

Liang-Qun Li conceived and designed the experiments; Xiao-Li Wang performed the experiments; Zong-Xiang Liu analyzed the simulated results; Wei-Xin Xie contributed analysis tools; Liang-Qun Li and Xiao-Li Wang wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Battistini, S.; Shima, T. Differential games missile guidance with bearings-only measurements. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 2906–2915. [Google Scholar] [CrossRef]
  2. Li, W.L.; Jia, Y.M. An information theoretic approach to interacting multiple model estimation. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 1811–1825. [Google Scholar] [CrossRef]
  3. Zhang, Q.; Song, T.L. Improved Bearings-Only Multi-Target Tracking with GM-PHD Filtering. Sensors 2016, 16, 1469. [Google Scholar] [CrossRef] [PubMed]
  4. Zhao, Y.; Patwari, N. Robust estimators for variance-based device-free localization and tracking. IEEE Trans. Mob. Comput. 2015, 14, 2116–2129. [Google Scholar] [CrossRef]
  5. Viani, F.; Robol, F.; Polo, A.; Rocca, P.; Oliveri, G.; Massa, A. Wireless architectures for heterogeneous sensing in smart home applications-Concepts and real implementations. Proce. IEEE 2013, 101, 2381–2396. [Google Scholar] [CrossRef]
  6. Li, L.; Xie, W. A Novel Auxiliary Quadrature Particle Filtering Algorithm Based on Target Characteristic. Acta Electron. Sin. 2014, 42, 2069–2074. (In Chinese) [Google Scholar]
  7. Blom, H.A.P.; Bar-shalom, Y. The interacting multiple model algorithm for systems with markovian switching coefficients. IEEE Trans. Autom. Control 1988, 33, 780–783. [Google Scholar] [CrossRef]
  8. Syamantak, D.G.; Yu, J.Y.; Mahendra, M.; Mark, C.; Mark, M. Comparison of Angle-only Filtering Algorithms in 3D Using EKF, UKF, PF, PFF, and Ensemble KF. In Proceedings of the 18th International Conference on Information Fusion Washington, Washington, DC, USA, 6–9 July 2015; pp. 1649–1656. [Google Scholar]
  9. Gustafsson, F.; Hendeby, G. Some relations between extended and unscented Kalman filters. IEEE Trans. Signal Process. 2012, 60, 545–555. [Google Scholar] [CrossRef]
  10. Julier, S.J.; Uhlmann, J.K. Unscented filtering and nonlinear estimation. Proce. IEEE 2004, 92, 401–422. [Google Scholar] [CrossRef]
  11. Ángel, F.G.-F.; Mark, R.M.; Jesús, G. Truncated Unscented Kalman Filtering. IEEE Trans. Signal Process. 2012, 60, 3372–3386. [Google Scholar]
  12. García-Fernández, A.F.; Morelande, M.R.; Grajal, J. Mixture truncated unscented Kalman filtering. In Proceedings of the 15th International Conference on Information Fusion, Singapore, 9–12 July 2012; pp. 479–486. [Google Scholar]
  13. Straka, O.; Duník, J.; Šimandl, M. Truncation nonlinear filters for state estimation with nonlinear inequality constraints. Automatica 2012, 48, 273–286. [Google Scholar] [CrossRef]
  14. Boada, B.L.; Garcia-Pozuelo, D.; Boada, M.J.L.; Diaz, V. A Constrained Dual Kalman Filter Based on pdf Truncation for Estimation of Vehicle Parameters and Road Bank Angle: Analysis and Experimental Validation. IEEE Trans. Intell. Transp. Syst. 2017, 18, 1006–1016. [Google Scholar] [CrossRef]
  15. Gordon, N.J.; Salmond, D.J.; Smith, A.F.M. Novel approach to nonlinear/non-Gaussian Bayesian state estimation. IEEE Proceed. Radar Sonar Navig. 1993, 140, 107–113. [Google Scholar] [CrossRef]
  16. Sutharsan, S.; Kirubarajan, T.; Lang, T.; McDonald, M. An Optimization-Based Parallel Particle Filter for Multitarget Tracking. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 1601–1618. [Google Scholar] [CrossRef]
  17. Hlinka, O.; Hlawatsch, F.; Djuric, P.M. Distributed Particle Filtering in Agent Networks. IEEE Signal Process. Mag. 2013, 30, 61–81. [Google Scholar] [CrossRef]
  18. Yang, T.; Laugesen, R.S.; Mehta, P.G.; Meyn, S.P. Multivariable feedback particle filter. Automatica 2016, 71, 10–23. [Google Scholar] [CrossRef]
  19. Boers, Y.; Driessen, J.N. Interacting multiple model particle filter. IEEE Proc. Radar Sonar Navig. 2003, 150, 344–349. [Google Scholar] [CrossRef]
  20. Keshavarz-Mohammadiyan, A.; Khaloozadeh, H. Adaptive IMMPF for bearing-only maneuvering target tracking in Wireless Sensor networks. In Proceedings of the 4th International Conference on Control, Instrumentation, and Automation (ICCIA), Qazvin, Iran, 27–28 Janurary 2016; pp. 6–11. [Google Scholar]
  21. Li, L.; Xie, W.; Huang, J.; Huang, J. Multiple Model Rao–Blackwellized Particle Filter for Maneuvering Target Tracking. Int. J. Def. Sci. 2009, 59, 197–204. [Google Scholar]
  22. Yu, J.Y.; Coates, M.J.; Rabbat, M.G.; Blouin, S. A Distributed Particle Filter for Bearings-Only Tracking on Spherical Surfaces. IEEE Signal Process. Lett. 2016, 23, 326–330. [Google Scholar] [CrossRef]
  23. Taghavi, E.; Tharmarasa, R.; Kirubarajan, T.; McDonald, M. Multisensor-multitarget bearing-only sensor registration. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 1654–1666. [Google Scholar] [CrossRef]
  24. Li, L.; Xie, W.; Liu, Z. Bearings-only maneuvering target tracking based on truncated quadrature Kalman filtering. Int. J. Electron. Commun. 2015, 69, 281–289. [Google Scholar]
  25. Ienkaran, A, Simon, H. Discrete-Time Nonlinear Filtering Algorithms Using Gauss-Hermit Quadrature. Proc. IEEE 2007, 95, 953–976. [Google Scholar]
Figure 1. Auxiliary truncated unscented Kalman filtering.
Figure 1. Auxiliary truncated unscented Kalman filtering.
Sensors 17 00972 g001
Figure 2. Root-mean-square (RMS) position errors of the extended Kalman filter (EKF), unscented Kalman filter (UKF), quadrature Kalman filtering (QKF), PF, mixture truncated unscented Kalman filter (MTUKF) and auxiliary truncated unscented Kalman filtering (ATUKF).
Figure 2. Root-mean-square (RMS) position errors of the extended Kalman filter (EKF), unscented Kalman filter (UKF), quadrature Kalman filtering (QKF), PF, mixture truncated unscented Kalman filter (MTUKF) and auxiliary truncated unscented Kalman filtering (ATUKF).
Sensors 17 00972 g002
Figure 3. RMS position error for different noise variances.
Figure 3. RMS position error for different noise variances.
Sensors 17 00972 g003
Figure 4. True Trajectory.
Figure 4. True Trajectory.
Sensors 17 00972 g004
Figure 5. Root-mean-square error (RMSE) of the ATUKF, IMMEKF and IMMRBPF. (a) X RMSE; (b) Y RMSE; (c) Z RMSE; (d) position RMSE.
Figure 5. Root-mean-square error (RMSE) of the ATUKF, IMMEKF and IMMRBPF. (a) X RMSE; (b) Y RMSE; (c) Z RMSE; (d) position RMSE.
Sensors 17 00972 g005
Table 1. Comparison of the computation times of different filtering algorithms (s). UNGM: univariate nonstationary growth model.
Table 1. Comparison of the computation times of different filtering algorithms (s). UNGM: univariate nonstationary growth model.
CaseEKFUKFQKFPFMTUKF(3)ATUKF
UNGM1.1026.65015.264522.51943.14216.240
Table 2. Computation times for all algorithms (s). BOT: Bearings-only maneuvering tracking.
Table 2. Computation times for all algorithms (s). BOT: Bearings-only maneuvering tracking.
CaseIMMEKFIMMRBPFTQKFATUKF
BOT0.07414.4930.5530.150

Share and Cite

MDPI and ACS Style

Li, L.-Q.; Wang, X.-L.; Liu, Z.-X.; Xie, W.-X. Auxiliary Truncated Unscented Kalman Filtering for Bearings-Only Maneuvering Target Tracking. Sensors 2017, 17, 972. https://doi.org/10.3390/s17050972

AMA Style

Li L-Q, Wang X-L, Liu Z-X, Xie W-X. Auxiliary Truncated Unscented Kalman Filtering for Bearings-Only Maneuvering Target Tracking. Sensors. 2017; 17(5):972. https://doi.org/10.3390/s17050972

Chicago/Turabian Style

Li, Liang-Qun, Xiao-Li Wang, Zong-Xiang Liu, and Wei-Xin Xie. 2017. "Auxiliary Truncated Unscented Kalman Filtering for Bearings-Only Maneuvering Target Tracking" Sensors 17, no. 5: 972. https://doi.org/10.3390/s17050972

APA Style

Li, L. -Q., Wang, X. -L., Liu, Z. -X., & Xie, W. -X. (2017). Auxiliary Truncated Unscented Kalman Filtering for Bearings-Only Maneuvering Target Tracking. Sensors, 17(5), 972. https://doi.org/10.3390/s17050972

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