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

: 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 ﬁlter is proposed by maximizing the mixture-correntropy-based cost. To improve numerical stability, the information ﬁlter 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


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,

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: e (e)de = 1 1−α logE e p α−1 e (e) = 1 1−α log(V α,σ (e)) (1) 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: Based on Expression (2), the correntropy can be defined as a generalized similarity measure between two random variables as follows: where κ(·, ·) is a shift-invariant kernel function. Without loss of generality, the Gaussian kernel function is used in this paper, which is given by: 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: 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: 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: 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.

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: 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: Multiple radar fusion scenario using Doppler measurements.
Assume that the state vector of the target is a variable in space ⊆ ℝ and can be denoted as is the position of the target, [˙,˙] are the speed, and the notation T represents the transpose operation.
The Doppler measurement is generated by each sensor in the scene, = 1, … , , 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 , y ].The multi-static Doppler frequency shift measured by two sensors is taken in space ⊆ ℝ , the measurement results are transmitted to the data fusion center, and the measurement equation shown in Formulas (3) and (4) are given as: xxxx H (9) and , is the measurement noise of sensor :

Mixture Correntropy Cost Function
Consider the discrete nonlinear state-space model: Assume that the state vector of the target is a variable in space X ⊆ R n x and can be denoted as 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: 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: and ε k,i is the measurement noise of sensor i:

Mixture Correntropy Cost Function
Consider the discrete nonlinear state-space model: 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 isx 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 2m . 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 statex k|k−1 and the predicted measurementŷ k|k−1,n can then be calculated as: Similarly, the predicted state covariance and the state-measurement cross-covariance matrix can then be calculated using the sigma points as: Define the measurement slope matrix as: H k,n = P −1 k|k−1 P xy,n T Then, the measurement equation for sensor n can be approximated by a statistical linear regression model as: 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: represent the augmented multi-source measurements and measurement matrix, respectively, and is the multi-source measurement noise vector with covariance and the mixture correntropy between predicted state and measurements can then be defined as: where e n,j is the jth 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 solutionx k can be obtained by minimizing the combined cost function based on correntropy according to the maximum correntropy criterion: where

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.
where , is the th element of , , is the selected kernel width for sensor according to the noise distribution characteristics, and is the dimension of the measurement for sensor .
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 ̂ can be obtained by minimizing the combined cost function based on correntropy according to the maximum correntropy criterion: where ∥ ∥ 2 = T .

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 −1 and information vector −1 at time − 1, and the state estimate can be represented as ̂− 1 = −1 −1 −1 and the associated covariance is −1 = −1 −1 . Then, the predicted state vector ̂| −1 and associated covariance | −1 are calculated using (12) and (13). Differentiating Equation (16) with respect to , the optimal estimate ̂ satisfies the equality: where , = , where , is the th element of , . The posterior information matrix is: 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 asx k−1 = V −1 k−1 v k−1 and the associated covariance is P k−1 = V −1 k−1 . Then, the predicted state vectorx 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 estimatex k satisfies the equality: where ζ k,n = R − 1 2 k,n ŷ k|k−1,n + H k,nxk − H k,nxk|k−1 − y k,n , and the term Λ −1 k|k−1 and Λ k,n are calculated as: where e j k,n is the jth element of ζ k,n . The posterior information matrix is: and the posterior information vector is: The detailed derivation is provided in Appendix A. Then, the posterior state estimate is: The state estimate x k is involved in the computation of Λ k,n , which implies that both sides of Equation (17) are interdependent withx 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 estimatex 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 MCCbased 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 1/2 k = QR. The Cholesky factor can also be computed using the QR factorization of P T/2 k−1 , 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 (2) Time update: Evaluate the sigma points χ i k−1 , and predict the points χ i k−1 and predicted state valueŝ x k|k−1 by using the process equation.
(5) Compute the term Λ k,y using (19). (6) Calculate the posterior information matrix: Calculate the posterior information matrix and the information vector: (7) Repeat steps 4-6 until convergence is achieved.

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: where m is the state dimension, · 2 is the 2-norm of matrix, h i,k,n is the ith row of H k,n , and v i,k is the element of the vector R According to (26), we obtain a function of kernel width: which satisfies ||g(x k )|| 1 ≤ φ(σ). Moreover, when σ n → ∞ , φ(σ) → 0, and when ∏ N n=1 σ 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 n=1 , and for any upper bound β of x 1 , we have: Differencing Equation (27) with respect to x k , when the kernel widths are large, we have: When the selected kernel widths σ n } N n=1 meet the conditions in (28) and (29) at the same time, the convergence of the fixed point iterative algorithm is guaranteed.

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: when the kernel width σ → ∞ , the term G σ n e i k,n can be approximated by the sum of zero order and first-order terms in Taylor expansion: Then, the cost function in (19) is as follows: which can be represented as follows: . 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 n=1 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.

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: for the CV model, and 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. proposed in this paper. In practical applications, the target signal is accumulated within the 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.
When the radar carrier frequency is 1 GHz, consider four kinds of non-Gaussian measurement noise: 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: and the covariance in the distribution can be represented as: When the radar carrier frequency is 1 GHz, consider four kinds of non-Gaussian measurement noise: In these equations, the first kind of noise is Gaussian noise with covariance 1000 Hz 2 ); the third is mixture Gaussian noise with covariance R 1 = diag(20 m 2 , 10 −5 rad 2 , 10 Hz 2 ), R 2 = diag 5 m 2 , 5 × 10 −6 rad 2 , 5 Hz 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. Figures 4 and 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. In these equations, the first kind of noise is Gaussian noise with covariance 1 = diag(20 m 2 , 10 −5 rad 2 , 10 Hz 2 ); the second is Gaussian noise with outliers, in which the outliers are generated with probability 0.1 and covariance 3 = diag(2000 m 2 , 10 −3 rad 2 , 1000 Hz 2 ); the third is mixture Gaussian noise with covariance 1 = diag(20 m 2 , 10 −5 rad 2 , 10 Hz 2 ) , 2 = diag(5 m 2 , 5 × 10 −6 rad 2 , 5 Hz 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. Figures 4 and 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 Figures 6  and 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.   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 Figures 6  and 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 correntropybased 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.  Evaluate the estimation performance with different kernel widths. Outliers generated with the parameters in case 2 are contained in the measurements. As shown in Figures 6  and 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.

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:   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.

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:  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.

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 Figures 9 and 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. In Figures 9 and 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.

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. In Figures 9 and 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.

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.

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.