Next Article in Journal
A Comprehensive Study of ChatGPT: Advancements, Limitations, and Ethical Considerations in Natural Language Processing and Cybersecurity
Previous Article in Journal
InterviewBot: Real-Time End-to-End Dialogue System for Interviewing Students for College Admission
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Sensor Fusion Target Tracking Based on Maximum Mixture Correntropy in Non-Gaussian Noise Environments with Doppler Measurements

School of Astronautics and Aeronautics, Shanghai Jiao Tong University, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Information 2023, 14(8), 461; https://doi.org/10.3390/info14080461
Submission received: 10 July 2023 / Revised: 4 August 2023 / Accepted: 10 August 2023 / Published: 15 August 2023

Abstract

:
This paper addresses the multi-sensor fusion target tracking problem based on maximum mixture correntropy in non-Gaussian noise environments exclusively using Doppler measurements. As Doppler measurements are non-linear, a statistical linear regression model is constructed using the unscented transformation. Then, a centralized measurement model is developed, and the mixture correntropy is determined, which contains the high-order statistics of state prediction and the measurement error caused by noise. Then, a robust fusion filter is proposed by maximizing the mixture-correntropy-based cost. To improve numerical stability, the information filter and corresponding square root version are also derived. Furthermore, the performance of the proposed algorithm is analyzed, and the selection of the kernel width is discussed. Experiments are performed using simulated data and automatic driving software. The results show that the estimation performance of the proposed algorithm is better with respect to outliers and mixture Gaussian noise than that of traditional methods.

1. Introduction

Doppler tracking is a technique that is widely used in the fields of radar, wireless communication, and acoustics, and it is also used for the motion state tracking of targets. In this technique, by employing the Doppler effect, the speed and direction of the movement of a target are inferred by analyzing the frequency variation of the received signal. Doppler tracking plays an important and well-established role in radar systems: it can help a radar system to detect and track moving targets in real time. Doppler tracking offers the advantages of high-precision speed measurement, motion state estimation, real-time performance, anti-interference ability, and wide applicability. These properties make it an important technique for object tracking and motion analysis in many fields [1,2].
To solve a target-tracking problem using Doppler measurements, Ristic and Farina implemented Bernoulli PF for Doppler tracking in multi-static situations in [3]. These authors utilized fixed single-transmitter and Doppler measurements measured by multiple receivers in order to jointly detect and track targets. As the Doppler measurements were non-linear, the solution was obtained via particle implementation. Ristic and Farina also proposed an active receiver-based approach in [4] to improve the information collected from sensors, and to discard measurements from sensors aligned with the target and transmitter in the bistatic case. This study shows that the choice of receiver improves track-settling time and steady-state error performance. The observability and performance of dual/multistatic Doppler radars were analyzed by Xiao in [5]. For simplicity, the author analyzed the multi-transmitter-receiver case, and pointed out that this case can also be generalized to multi-transmitter and multi-receiver cases. In the case of a single receiver, observability is evaluated using the rank of the Fisher Information Matrix (FIM). If a target is completely observable, the FIM should have the highest rank. These Doppler-tracking filter methods and analyses for target tracking are all based on Gaussian assumptions.
A Kalman filter (KF) is an optimal estimator under the minimum mean square error criterion, for which a solution is achieved by recursively updating the estimated state and covariance matrix based on a state-space model. For nonlinear systems, various sub-optimal nonlinear extensions have been proposed. However, with respect to dealing with non-Gaussian noise, such as impulsive interference or outliers, their estimation performance is poor, as the minimum mean square error cannot indicate the bias of the predicted state. Accordingly, Robust filters, such as the H filter and Huber’s filter [6,7,8], are employed in this situation. In particular, the H filter ensures bounded estimation error but performs poorly in relation to Gaussian noise. Alternatively, Huber’s filter is suitable for dealing with data estimation problems involving outliers. Based on the theory of robust statistics, Huber’s filter has a certain robustness to outliers, and it can resist the interference of outliers. The Huber–Kalman filter can be used to obtain optimal results by minimizing a combined l 1 and l 2 norm. Huber estimation can provide better estimation results in the presence of outliers, and it has a better fitting effect on non-outliers.
Recently, a new cost function based on information theory was constructed, and a corresponding optimal criterion was proposed [9]. Correntropy is defined as a local similarity measure, and it is used as a cost for state estimation [10]. By maximizing the correntropy-based cost, the maximum correntropy criterion (MCC) was employed to derive robust solutions. As correntropy costs contain the high-order statistics of the error signal, the high-order characteristics of the error distribution can be captured, and the derived solutions are, therefore, more robust. MCC-based Kalman filters (MCC-KFs) can be used to solve linear and non-linear state estimation problems involving non-Gaussian noise, such as outliers and impulsive noise, and the performance is obviously better than traditional methods [11,12]. In fact, MCC-based filters can also be treated as a smoothed MAP estimator [13]. However, MCC-based filters cannot be obtained in closed forms; therefore, the solutions are usually obtained using a fixed-point iteration approach [14]. The selection of kernel width is also important in MCC-based algorithms, as it determines the relative weighting of different orders of norms in the loss function, thus potentially altering the estimation results. In [15], the algorithm performance under different kernel widths was shown and analyzed; better estimation performance was obtained using an appropriate kernel width.
This paper addresses the problem of multi-sensor fusion target tracking in non-Gaussian noise environments with Doppler measurements based on maximum mixture correntropy. A statistical linear regression model is constructed for nonlinear ranges, angles, and Doppler measurements by employing unscented transformation. Based on the centralized measurement model, a mixture-correntropy-based cost function is constructed to capture the high-order statistics of state prediction and the measurement error that are caused by non-Gaussian noise. Subsequently, a robust fusion filter is developed by maximizing the cost function based on the maximum correntropy criterion. The solution is obtained through an iterative approach and the use of square roots. To enhance numerical stability, the information filter version is also derived. This paper also analyzes the performance of the proposed algorithm and discusses the selection of kernel width. Experiments are performed using numerical simulation and automatic driving software; the results demonstrate that the proposed method is more robust against non-Gaussian noise than conventional methods.
The maximum correntropy criterion is introduced in Section 2. In Section 3, the mixture correntropy cost function is derived, and a robust fusion filter is developed. The simulation results are presented in Section 4. Finally, the conclusions of the study are summarized in Section 5.

2. Background

2.1. Maximum Correntropy Criterion

Define e = y x is the error between two variables x and y with probability distribution p e e . The information of the error e can be measured according to the definition of Renyi’s entropy, which is denoted as:
H α e 1 1 α l o g p e α e d e = 1 1 α l o g E e p e α 1 e = 1 1 α l o g V α , σ e
where E denotes the expectation operator. The variable α is the order of the function. Therefore, Expression (1) is called Renyi’s α entropy, and V α , σ e is α information potential estimator.
When α = 2 , the term H α e is quadratic Renyi’s entropy:
H 2 e = l o g E e p e = l o g V 2 e
Based on Expression (2), the correntropy can be defined as a generalized similarity measure between two random variables as follows:
H x , y = l o g κ x , y p x , y d x d y = l o g E X Y G σ x y
where κ , is a shift-invariant kernel function. Without loss of generality, the Gaussian kernel function is used in this paper, which is given by:
κ x , y = G σ e = e x p e 2 2 σ 2
where e = x y , and σ > 0 is the kernel width.
Usually, the joint distribution p x , y is unknown, and only limited numbers of data are available. In this situation, the PDF p x , y can be estimated from data using Parzen’s window estimator, and the expectation operator can be approximated with the sample mean. The nonparametric estimation of quadratic information potential can be calculated as:
V x , y = 1 N i = 1 N G σ e i
where G σ ( · ) represents the Gaussian kernel function. Based on the optimization theory, the optimal solution is obtained when the similarity of two variables is maximized with respect to e i , such as, for example:
x ^ = a r g m i n x H 2 e
This criterion is called the maximum correntropy criterion. Notice that the logarithm is dropped, and, therefore, the derivation of (6) is equal to maximize the term V e , such as, for example:
m i n x H 2 e = m a x x V e = m a x x 1 N i = 1 N G σ e i
It can be seen that the term G σ ( e i ) contains the high moments of the error. Thus, the MCC algorithm derives the optimal solution in high-dimensional space, and the results are more robust given non-Gaussian noise.

2.2. Doppler Measurement Equation

Consider the multi-radar fusion situational awareness background, assuming that the data fusion center knows the position of the transmitter and the sensor, the moving target detected by the sensor gives each Doppler measurement, and the sensor outputs the measured value to the data fusion center. The scene is illustrated in Figure 1:
Assume that the state vector of the target is a variable in space X R n x and can be denoted as x k = x k , y k , x ˙ k , y ˙ k T , in which x k , y k is the position of the target, x ˙ k , y ˙ k are the speed, and the notation T represents the transpose operation.
The Doppler measurement is generated by each sensor in the scene, i = 1 , , N s , and the motion of the target in the direction of the transmitter and receiver jointly cause the Doppler effect, which is located at the i-th of x i , y i .The multi-static Doppler frequency shift measured by two sensors is taken in space Z R n z , the measurement results are transmitted to the data fusion center, and the measurement equation shown in Formulas (3) and (4) are given as:
z k , i = h i x k + ε k , i
Among them, h i x k = d ˙ k t λ + d ˙ k , i λ is the Doppler frequency shift, λ is the wavelength of the transmitted signal, and ε k , i is the measurement noise of the sensor. The linearization matrix of h i x k can be expressed as:
H k , i = h i x k x k h i x k y k h i x k x ˙ k h i x k y ˙ k
and ε k , i is the measurement noise of sensor i :
ε k , i N ε ; 0 , σ ε 2

3. Robust Fusion Filter Based on Maximum Mixture Correntropy

3.1. Mixture Correntropy Cost Function

Consider the discrete nonlinear state-space model:
x k = f x k 1 + w k 1 y k , n = h n x k + v k , n
where f x and h n x are the nonlinear state transition equation and the nonlinear measurement equation with an independent non-linear process and a measurement noise of w k 1 and v k , n with the associated covariance Q k 1 and R k , n . Assume that at time k 1 , the prior state is x ^ k 1 with the associated covariance P k 1 , and the dimension of the state vector is m . To approximate the nonlinear of the state transition process and the measurement, a set of 2 m + 1 sigma points { χ i , k 1 } i = 1 2 m + 1 are generated based on the unscented transform [16].
According to the state transition equation and the nonlinear measurement equations, the predicted sigma points χ i , k | k 1 , γ i , k | k 1 , n can be obtained, and the predicted state x ^ k | k 1 and the predicted measurement y ^ k | k 1 , n can then be calculated as:
x ^ k | k 1 = 1 2 m + 1 i = 1 2 m + 1 χ i , k | k 1 y ^ k | k 1 , n = 1 2 m + 1 i = 1 2 m + 1 γ i , k | k 1 , n
Similarly, the predicted state covariance and the state-measurement cross-covariance matrix can then be calculated using the sigma points as:
P k | k 1 = i = 1 2 m + 1 χ i , k | k 1 χ i , k | k 1 T x ^ k | k 1 x ^ k | k 1 T + Q k 1 P x y , n = i = 1 2 m + 1 χ i , k | k 1 γ i , k | k 1 , n T x ^ k | k 1 y k | k 1 , n T
Define the measurement slope matrix as:
H k , n = P k | k 1 1 P x y , n T
Then, the measurement equation for sensor n can be approximated by a statistical linear regression model as:
y k , n y ~ k , n = y ^ k | k 1 , n + H k , n x k x ^ k | k 1 + v k , n
This equation is derived based on Kalman filtering equations [17], which provide the pseudo-measurement to approximate the nonlinear observation model for state update and filtering operations with unscented transformation. The pseudo-measurement y ~ k , n is an approximate measure obtained by predicting a state estimate and linearizing the observation model in order to transform the nonlinear observation model into a linear model. Pseudo-measurements can help us to deal with both nonlinear and non-Gaussian noise in estimation problems.
Then, the centralized measurement equation [16] can be constructed as:
Y k = H k x k + v k
where the vector Y k = y k , 1 T , y k , 2 T , , y k , n T T and matrix H k = H k , 1 T , H k , 2 T , , H k , n T T represent the augmented multi-source measurements and measurement matrix, respectively, and v k = v k , 1 T , v k , 2 T , , v k , n T T is the multi-source measurement noise vector with covariance R k = diag R 1 , R 2 , , R n .
Let ζ k , n = R k , n 1 Y k | k 1 + H k x k H k x ^ k | k 1 Y k , and the mixture correntropy between predicted state and measurements can then be defined as:
V x y = n = 1 N 1 L n j = 1 L n G σ n e n , j
where e n , j is the j th element of ζ k , n , σ n is the selected kernel width for sensor n according to the noise distribution characteristics, and L n is the dimension of the measurement for sensor n .
Assuming independence between the dynamic process noise and the measurement noise, the MMC criteria are applied to the state transition and the measurement models. The optimal solution x ^ k can be obtained by minimizing the combined cost function based on correntropy according to the maximum correntropy criterion:
x ^ k = a r g m i n x k G σ n ( x k x ^ k | k 1 ) n = 1 N 1 L n j = 1 L n G σ n e n , j
where x A 2 = x T A x .

3.2. Robust Fusion Information Filter Based on Maximum Correntropy Criterion

As shown in Figure 2, we provide the structural diagram of the proposed algorithm, which has a recursive framework. In this section, an iterative algorithm is proposed for multi-sensor fusion estimation with non-Gaussian noise. The mixture correntropy cost is constructed and recursively calculated as it related to the state estimates.
To improve the numerical stability, the fusion estimation algorithm within the information filter framework is proposed in this section.
Start with the Fisher information matrix V k 1 and information vector v k 1 at time k 1 , and the state estimate can be represented as x ^ k 1 = V k 1 1 v k 1 and the associated covariance is P k 1 = V k 1 1 . Then, the predicted state vector x ^ k | k 1 and associated covariance P k | k 1 are calculated using (12) and (13).
Differentiating Equation (16) with respect to x k , the optimal estimate x ^ k satisfies the equality:
Λ k | k 1 1 x ^ k x ^ k | k 1 = n = 1 N H k , n T R k , n T 2 Λ k , n R k , n 1 2 ζ k , n
where ζ k , n = R k , n 1 2 y ^ k | k 1 , n + H k , n x ^ k H k , n x ^ k | k 1 y k , n , and the term Λ k | k 1 1 and Λ k , n are calculated as:
Λ k = diag G σ x k x ^ k | k 1
Λ k , n = diag G σ e k , n 1 , G σ e k , n 2 , , G σ e k , n j
where e k , n j is the j th element of ζ k , n .
The posterior information matrix is:
V k = P k | k 1 1 + H k T Λ ~ k , y H k = V k | k 1 + n = 1 N H k , n T Λ ~ k , n H k , n
and the posterior information vector is:
v k P k 1 x ^ k = v k | k 1 + n = 1 N H k , n T Λ ~ k , n y k , n y ^ k | k 1 H k x ^ k | k 1
The detailed derivation is provided in Appendix A.
Then, the posterior state estimate is:
x ^ k = V k 1 v k
The state estimate x k is involved in the computation of Λ k , n , which implies that both sides of Equation (17) are interdependent with x ^ k . In order to update the estimates in the Kalman filter, a fixed-point iteration approach is employed where the innovation covariance matrix and the posterior information matrix and vector are recalculated iteratively until convergence is achieved. Using the posterior estimate x ^ k obtained from Equation (22), the terms Λ k and Λ k , n are recalculated, and the posterior information matrix and vector are updated iteratively.
The robust square root (SR) fusion information filter is a modification of the MCC-based fusion information filter using Cholesky factors. Then, re-formulate the filtering equations, and a more numerically stable algorithm will be derived than the original method.
The SR implementation of the fusion information filter uses QR factorization in each iteration step to update the corresponding Cholesky factors. This involves decomposing the information matrix into an orthogonal matrix Q and an upper triangular matrix R, such that V k 1 / 2 = Q R . The Cholesky factor can also be computed using the QR factorization of P k 1 T / 2 , where T denotes the transpose operation.
Algorithm: Maximum mixture correntropy fusion filter:
(1) Initialize V k 1 and vector v k 1 of the target state, and the posterior covariance is then P k 1 1 / 2 = V k 1 1 / 2 and the state vector is x ^ k 1 = P k 1 T / 2 v k 1 .
(2) Time update:
Evaluate the sigma points χ k 1 i , and predict the points χ k 1 i and predicted state values x ^ k | k 1 by using the process equation.
(3) Calculate the square root Λ k | k 1 1 / 2 using QR decomposition:
Λ k | k 1 1 / 2 , 0 T = Q R χ k k 1 , Q k 1 1 / 2 T
Calculate V k | k 1 1 / 2 = Λ k | k 1 1 / 2 and v k | k 1 1 / 2 = Λ k | k 1 1 / 2 x k | k 1 .
(4) Measurement update:
Calculate γ k + 1 k and the measurement slope matrix H k = χ k + 1 k 1 γ k + 1 k using the measurement equation to obtain the regression model (16).
(5) Compute the term Λ k , y using (19).
(6) Calculate the posterior information matrix:
V k 1 / 2 0 T = Q R V k | k 1 1 / 2 H k , 1 T Λ ~ k , 1 1 / 2 H k , n T Λ ~ k , n 1 / 2 T
where Λ ~ k , n T / 2 Λ ~ k , n 1 / 2 = Λ ~ k , n . Calculate the posterior information matrix and the information vector:
v k | k = v k | k 1 + n = 1 N H k , n T Λ ~ k , y y k , n y ^ k | k 1 H k , n x ^ k | k 1
(7) Repeat steps 4–6 until convergence is achieved.

3.3. Analysis of the Algorithm Convergence

The kernel width σ is a crucial parameter, and the selection of this parameter affects the convergence of both the algorithm and the fusion estimation performance.
In this section, we analyze the effect of parameter selection on the convergence of the iteration implementation of the proposed algorithm. Since the MCC solution cannot be obtained in closed form, iterative methods are used to derive the final results. The convergence of the fixed point iteration approach is guaranteed by the compression mapping theorem when the function g x k 1 and its derivative x k g x k 1 are bounded. Here, we define β as the upper bound of x ^ 1 , and we ensure that the function g x k 1 satisfies this inequality:
g x k 1 a Λ k | k 1 1 + H k T Λ ~ k , y H k 1 1 H k T Λ ~ k , y Y k Y ^ k 1 b Λ k | k 1 1 + H k T Λ ~ k , y H k 1 1 n = 1 N i = 1 L n v i , k × h i , k , n c m Λ k | k 1 1 + H k T Λ ~ k , y H k 1 2 n = 1 N i = 1 L n v i , k × h i , k , n = m λ m a x Λ k | k 1 1 + H k T Λ ~ k , y H k 1 n = 1 N i = 1 L n v i , k × h i , k , n d m n = 1 N i = 1 L n v i , k × h i , k , n λ m i n P k | k 1 1 + n = 1 N i = 1 L n ξ n h i , k , n T h i , k , n
where m is the state dimension, · 2 is the 2-norm of matrix, h i , k , n is the i th row of H k , n , and v i , k is the element of the vector R k , n 1 2 y k , n y ^ k | k 1 , n . The term ξ n = G σ n β h i , k , n 1 + d i , k , where d i , k is the element of the vector R k , n 1 2 y k , n y ^ k | k 1 , n H k , n x ^ k | k 1 . The inequality (a) and (c) hold according to the properties of the matrix norm, and the inequality (b) is established based on the triangle inequality. The inequality (d) holds as ξ n is a monotonically decreasing value of x .
According to (26), we obtain a function of kernel width:
ϕ σ = m n = 1 N i = 1 L n v i , k × h i , k , n λ m i n P k | k 1 1 + n = 1 N i = 1 L n ξ n h i , k , n T h i , k , n
which satisfies g x k 1 ϕ σ . Moreover, when σ n , ϕ σ 0 , and when n = 1 N σ n 0 , the function ϕ σ . The reason is that the function ϕ σ is continuous and involves the monotonic decreasing of σ n .
Therefore, the equation ϕ σ = β has the solution { σ n } n = 1 N , and for any upper bound β of x 1 , we have:
{ σ n } n = 1 N , s . t . g x 1 ϕ σ β
Differencing Equation (27) with respect to x k , when the kernel widths are large, we have:
x k g x k 1 α 1
When the selected kernel widths { σ n } n = 1 N meet the conditions in (28) and (29) at the same time, the convergence of the fixed point iterative algorithm is guaranteed.

3.4. Discussion of the Kernel Selection

Further, we analyze the influence of parameter selection on estimation performance. Extend the potential estimator around zero by the Taylor series yields:
n = 1 N 1 L n i = 1 L n G σ n e k , n , i = n = 1 N i = 1 L n 1 L n 2 π σ p = 0 1 p 2 p p ! e k , n , i 2 σ 2 p
when the kernel width σ , the term G σ n e k , n i can be approximated by the sum of zero order and first-order terms in Taylor expansion:
G σ n e k , n i 1 2 π σ 1 e k , n , j 2 2 σ n 2
Then, the cost function in (19) is as follows:
J F m i n x k x k | k 1 P k | k 1 1 2 n = 1 N 1 L n 2 π σ j = 1 L n e k , j 2
which can be represented as follows:
J F m i n x k x ^ k | k 1 P k | k 1 1 2 n = 1 N H k , n x k y k , n W k 2
where W = R k 1 L n 2 π σ . This equation is the same as the traditional MMSE-based fusion cost function, and the derived solution is similar to the conventional Kalman filter fusion algorithm.
The kernel widths σ x , { σ n } n = 1 N are crucial parameters in the maximum mixture correntropy (MMCC)-based fusion algorithm, and they can be set manually or optimized through trial-and-error in practical applications.

4. Simulation

4.1. Numerical Example

In this section, simulations are performed in order to demonstrate the effectiveness of the proposed algorithm, and they are then compared with traditional methods. In the scenario, consider linear constant velocity (CV) and nonlinear constant turn (CT) motion models with state transition equation:
x k = F k x k 1 + w k
where
F k = 1 Δ T 0 0 0 1 0 0 0 0 1 Δ T 0 0 1 0
for the CV model, and
F k = 1 s i n ω Δ T ω 0 1 c o s ω Δ T ω 0 c o s ω Δ T 0 s i n ω Δ T 0 1 c o s ω Δ T ω 1 s i n ω Δ T ω 0 s i n ω Δ T 0 c o s ω Δ T
for the CT model.
In the state transition matrix, the turn rate ω = 0.005 rad / s , and the dynamic state is initialed with x = p x , v x , p y , v y T = 0   m , 10   m / s , 0   m , 0   m / s T . The target track is shown in Figure 3.
In these equations, Δ T represents the interval between the target motion and the measurements, which also corresponds to the time step of the iterative estimation process proposed in this paper. In practical applications, the target signal is accumulated within the Δ T interval, and Doppler measurements are obtained by Fourier transform. Since this paper mainly studies the fusion estimation method with target Doppler measurement under non-Gaussian measurement conditions, the radar signal processing is not discussed in detail here.
The process noise w k follows mixture Gaussian distribution:
w k ~ N r 1 | 0 , Q 1 2 + 0.1 N r 3 | 0 , Q 3 2
and the covariance in the distribution can be represented as:
Q n = diag T 3 l n 3 T 2 l n 2 T 2 l n 2 T l n , T 3 l n 3 T 2 l n 2 T 2 l n 2 T l n
where l 1 = 10 3   m 2 / s 3 , l 2 = 10 3   m 2 / s 3 .
Four sensors receive Doppler measurements. As shown in Figure 1, the positions of the sensors are set at [−1500 m, −600 m], [−1500 m, 600 m], [500 m, −600 m], and [500 m, 600 m], respectively, and the measurement equation can be represented as:
h x = p x 2 + p y 2 arctan p y p x p x v x + p y v y p x 2 + p y 2
When the radar carrier frequency is 1 GHz, consider four kinds of non-Gaussian measurement noise:
r k , 1 N r 1 | 0 , R 1 2 r k , 2 N r 1 | 0 , R 1 2 + 0.1 N r 3 | 0 , R 3 2 r k , 3 N r 1 | 0 , R 1 2 + N r 2 | 0 , R 2 2 r k , 4 N r 1 | 0 , R 1 2 + N r 2 | 0 , R 2 2 + 0.1 N r 3 | 0 , R 3 2
In these equations, the first kind of noise is Gaussian noise with covariance R 1 = diag 20   m 2 , 10 5   rad 2 ,   10   H z 2 ; the second is Gaussian noise with outliers, in which the outliers are generated with probability 0.1 and covariance R 3 = diag 2000   m 2 , 10 3   rad 2 ,   1000   H z 2 ; the third is mixture Gaussian noise with covariance R 1 = diag 20   m 2 , 10 5   rad 2 ,   10   H z 2 , R 2 = diag 5   m 2 , 5 × 10 6   rad 2 ,   5   H z 2 ; and the fourth is mixture Gaussian noise with outliers.
In the simulation, the initial position and the velocity estimates are calculated by the Doppler cross approach. The termination threshold of the iteration implementation is set to 0.01. The proposed MMCC-based fusion method is compared with the traditional fusion method based on KF [17] and the robust Huber filter [8] in terms of the root mean squared error (RMSE) with 200 Monte Carlo trails. The selection of kernel widths is crucial for the performance here, and these parameters are chosen based on the error distribution through trial-and-error.
The fusion estimation results with Gaussian noise, mixture Gaussian noise, outliers, and complex noise are presented in Table 1, respectively. Appropriate kernel widths are selected in order to obtain the best performance for the correntropy-based fusion method. Figure 4 and Figure 5 show the estimation performance given mixture Gaussian noise and outliers. As shown in the figures, the MMCC-Fusion method outperforms the other algorithms. The reason for this is that the correntropy-based cost function in the MMCC-Fusion method captures higher-order statistics compared to the correction function in the Huber-Fusion algorithm, which allows it to accurately reflect the characteristics of the arbitrarily distributed process and the measurement noise. Therefore, the MMCC-Fusion-IF algorithm is able to effectively suppress different types of non-Gaussian noise under the MCC criterion. This experiment demonstrates the robustness and the adaptiveness of the correntropy-based method in complex non-Gaussian noise scenarios.
Table 1 provides the estimation performance given different noise patterns. It can be seen that KF is the optimal estimator, and that it performs best given Gaussian noise. However, the MMCC-Fusion method can still provide better results. The reason for this is that through select appropriate kernel width, relatively good results can be achieved under arbitrary noise conditions, and this also reflects the flexibility of the MMCC-based method. When the kernel width approaches infinity, the estimation results are close to the traditional Kalman filter fusion method.
Evaluate the estimation performance with different kernel widths. Outliers generated with the parameters in case 2 are contained in the measurements. As shown in Figure 6 and Figure 7 and also listed in Table 2, the estimation results are the best when the kernel widths are selected as σ = 9 for each Doppler radar, which leads to a desirable error distribution and optimal results. It can be seen that the performance decreases with inappropriate kernel widths because the characteristics cannot be accurate captured. With the increasing of kernel width, the second moment of the error has the largest proportion in the correntropy-based loss function. In this situation, the loss function is approximate to the minimum mean square error loss, and the estimation results will be close to the traditional Kalman filter fusion.
Furthermore, Table 2 shows that larger kernel widths result in faster convergence. With kernel widths larger than 10, the iterative implementation of the proposed algorithm converges to the optimal solution quickly in only two or three iterations. Therefore, very small kernel widths should be avoided to avoid slow convergence.

4.2. Autonomous Driving Data Simulation

This experiment provided the tracking of an autonomous driving vehicle. The data source is from the Udacity course self-driving car [18], and the Doppler measurements are obtained by four distributed radar sensors. The track of the car is shown in Figure 8:
In Figure 9 and Figure 10, the estimation results are shown. It can be seen that the KF-Fusion method performs worst, mainly due to the sensitivity of the MMSE criterion to outliers. The Huber-Fusion method yields better results by mitigating the impact of outliers using a correction function on the measurements. The MMCC-Fusion-IF method outperforms the others as non-Gaussian measurements are suppressed under the MCC criterion.

5. Conclusions

This paper addresses the multi-sensor fusion target tracking problem based on maximum mixture correntropy with the non-Gaussian process and measurement noise using Doppler measurements. A robust fusion filter is developed by maximizing the mixture correntropy-based cost, which contains the high-order statistics of state prediction and measurement error caused by noise. Furthermore, the selection of the kernel width is discussed. Simulations are performed given different non-Gaussian noise, such as outliers and mixture Gaussian noise. Experiments are performed using simulated data and automatic driving software. The results demonstrate the effectiveness of the proposed method.

Author Contributions

C.Y.: conceptualization, investigation, writing—review and editing. S.L.: writing—review and editing. M.L.: methodology, writing—original draft preparation. All authors have read and agreed to the published version of the manuscript.

Funding

This work is jointly supported by National Natural Science Foundation of China (Grant Nos. 61673262 and 61175028), China Postdoctoral Science Foundation(Grant No. 2019M651498).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Derivation of (20) and (21).
Let Λ ~ k , n = R k , n T 2 Λ k , n R k , n 1 2 , and Equation (17) can then be rewritten in a centralized form as:
Λ k | k 1 1 x ^ k x ^ k | k 1 = H k T Λ ~ k , y Y k Y ^ k + H x ^ k | k 1 H k x ^ k
where Λ ~ k , y = diag Λ ~ k , 1 , Λ ~ k , 2 , , Λ ~ k , n .
Change the form of (34) to obtain:
x ^ k = x ^ k | k 1 + K k Y k Y ^ k K k = Λ k | k 1 1 + H k T Λ ~ k , y H k 1 H k T Λ ~ k , y
Apply the matrix inversion lemma
A + B C D 1 = A 1 A 1 B C 1 + D A 1 B 1 D A 1
with the identifications
Λ k | k 1 A , H k T B , Λ ~ k , y C , H k D
and the innovation gain
K = Λ k | k 1 H k T Λ ~ k , y + H k P k | k 1 H k T 1
Then, the posterior covariance can be calculated as follows:
P k = I K H k P k | k 1 = P k | k 1 1 + H k T Λ ~ k , y H k 1
Therefore, the posterior information matrix is:
V k = P k | k 1 1 + H k T Λ ~ k , y H k = V k | k 1 + n = 1 N H k , n T Λ ~ k , n H k , n
and the posterior information vector is:
v k P k 1 x ^ k = P k | k 1 1 + H k T Λ ~ y , k H k x k | k 1 + P k 1 P k H k T Λ ~ k , y Y k Y ^ k = P k | k 1 1 x k | k 1 + H k T Λ ~ k , y Y k Y ^ k + H k x ^ k | k 1 = v k | k 1 + n = 1 N H k , n T Λ ~ k , n y k , n y ^ k | k 1 H k x ^ k | k 1

References

  1. Wu, W.; Jiang, J.; Liu, W.; Feng, X.; Gao, L.; Qin, X. Augmented state gm-phd filter with registration errors for multi-target tracking by doppler radars. Signal Process. 2016, 120, 117–128. [Google Scholar] [CrossRef]
  2. Zheng, J. Structure and Performance Analysis of Signal Acquisition and Doppler Tracking in LEO Augmented GNSS Receiver. Sensors 2021, 21, 525. [Google Scholar] [CrossRef] [PubMed]
  3. Ristic, B.; Farina, A. Joint detection and tracking using multi-static doppler-shift measurements. In Proceedings of the International Conference on Acoustics, Speech and Signal Processing (ICASSP) 2012, Kyoto, Japan, 25–30 March 2012; pp. 3881–3884. [Google Scholar]
  4. Ristic, B.; Farina, A. Target tracking via multi-static Doppler shifts. IET Radar Sonar Navig. 2013, 7, 508–516. [Google Scholar] [CrossRef]
  5. Vo, B.T.; Vo, B.N.; Cantoni, A. The Cardinality Balanced Multi-target Multi-Bernoulli Filter and Its Implementations. IEEE Trans. Signal Process. 2009, 57, 409–423. [Google Scholar]
  6. Karlgaard, C.D. Nonlinear Regression Huber-Kalman Filtering and Fixed-Interval Smoothing. J. Guid. Control Dyn. 2014, 38, 322–330. [Google Scholar] [CrossRef]
  7. Karlgaard, C.; Schaub, H. Huber-based divided difference filtering. J. Guid. Control Dyn. 2007, 30, 885–891. [Google Scholar] [CrossRef]
  8. Graham, M.C.; How, J.P.; Gustafson, D.E. Robust state estimation with sparse outliers. J. Guid. Control Dyn. 2015, 38, 1229–1240. [Google Scholar] [CrossRef] [Green Version]
  9. Principe, J.C. Information Theoretic Learning: Renyi’s Entropy and Kernel Perspectives; Springer: New York, NY, USA, 2010. [Google Scholar]
  10. Liu, W.; Pokharel, P.P.; Principe, J.C. Correntropy: Properties, and applications in non-Gaussian signal processing. IEEE Trans. Signal Process. 2007, 55, 5286–5298. [Google Scholar] [CrossRef]
  11. Chen, B.; Liu, X.; Zhao, H. Maximum correntropy Kalman filter. Automatic 2017, 76, 70–77. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, Y.; Zheng, W.; Sun, S.; Li, L. Robust Information Filter Based on Maximum Correntropy Criterion. J. Guid. Control Dyn. 2016, 39, 1126–1131. [Google Scholar] [CrossRef]
  13. Chen, B.; Principe, J.C. Maximum correntropy estimation is a smoothed MAP estimation. IEEE Signal Process. Lett. 2012, 19, 491–494. [Google Scholar] [CrossRef]
  14. Chen, B.; Wang, J.; Zhao, H. Convergence of a Fixed-Point Algorithm under Maximum Correntropy Criterion. IEEE Signal Process. Lett. 2015, 22, 1723–1727. [Google Scholar] [CrossRef]
  15. Wang, G.Q.; Gao, Z.X.; Zhang, Y.G.; Ma, B. Adaptive maximum correntropy Gaussian filter based on variational Bayes. Sensors 2018, 18, 1960. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Gan, Q.; Harris, C.J. Comparison of two measurement fusion methods for Kalman-filter-based multisensor data fusion. IEEE Trans. Aerosp. Electron. Syst. 2001, 37, 273–279. [Google Scholar] [CrossRef] [Green Version]
  17. Wang, G.; Li, N.; Zhang, Y. Maximum correntropy unscented Kalman and information filters for non-Gaussian measurement noise. J. Frankl. Inst. 2017, 354, 8659–8677. [Google Scholar] [CrossRef]
  18. Udacity’s Self-Driving Car Simulator, Udacity, Mountain View, CA, USA. 2017. Available online: https://github.com/udacity/self-drivingcar-sim (accessed on 2 August 2023).
Figure 1. Multiple radar fusion scenario using Doppler measurements.
Figure 1. Multiple radar fusion scenario using Doppler measurements.
Information 14 00461 g001
Figure 2. Algorithm framework.
Figure 2. Algorithm framework.
Information 14 00461 g002
Figure 3. Scenario 1.
Figure 3. Scenario 1.
Information 14 00461 g003
Figure 4. Comparison of position estimation performance with complex noise.
Figure 4. Comparison of position estimation performance with complex noise.
Information 14 00461 g004
Figure 5. Comparison of position estimation performance with complex noise.
Figure 5. Comparison of position estimation performance with complex noise.
Information 14 00461 g005
Figure 6. Position estimation performance with different kernel widths.
Figure 6. Position estimation performance with different kernel widths.
Information 14 00461 g006
Figure 7. Velocity estimation performance with different kernel widths.
Figure 7. Velocity estimation performance with different kernel widths.
Information 14 00461 g007
Figure 8. Scenario 2 of autonomous driving.
Figure 8. Scenario 2 of autonomous driving.
Information 14 00461 g008
Figure 9. Comparison of position estimation performance with outliers.
Figure 9. Comparison of position estimation performance with outliers.
Information 14 00461 g009
Figure 10. Comparison of velocity estimation performance with outliers.
Figure 10. Comparison of velocity estimation performance with outliers.
Information 14 00461 g010
Table 1. Fusion estimation RMSE in different non-Gaussian noise.
Table 1. Fusion estimation RMSE in different non-Gaussian noise.
AlgorithmCase 1Case 2Case 3Case 4
KF-Fusion1.4132.5101.9182.404
Huber-Fusion2.1652.0501.9072.226
MMCC-Fusion1.5551.7531.7251.890
Table 2. Velocity estimation performance with different kernel widths.
Table 2. Velocity estimation performance with different kernel widths.
AlgorithmKernel Widths RMSE   of   x ˙ RMSE   of   x ¨ Iteration Num
MMCC-Fusion-IF ( σ = 8 )1.8950.2493.85
( σ = 9 )1.8400.2342.98
( σ = 15 )2.1640.3062.36
KF-Fusion 2.1790.312
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yi, C.; Li, M.; Li, S. Multi-Sensor Fusion Target Tracking Based on Maximum Mixture Correntropy in Non-Gaussian Noise Environments with Doppler Measurements. Information 2023, 14, 461. https://doi.org/10.3390/info14080461

AMA Style

Yi C, Li M, Li S. Multi-Sensor Fusion Target Tracking Based on Maximum Mixture Correntropy in Non-Gaussian Noise Environments with Doppler Measurements. Information. 2023; 14(8):461. https://doi.org/10.3390/info14080461

Chicago/Turabian Style

Yi, Changyu, Minzhe Li, and Shuyi Li. 2023. "Multi-Sensor Fusion Target Tracking Based on Maximum Mixture Correntropy in Non-Gaussian Noise Environments with Doppler Measurements" Information 14, no. 8: 461. https://doi.org/10.3390/info14080461

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