Next Article in Journal
Philosophy and Application of High-Resolution Temperature Sensors for Stratified Waters
Next Article in Special Issue
Adaptive Fifth-Degree Cubature Information Filter for Multi-Sensor Bearings-Only Tracking
Previous Article in Journal
Developing Efficient Thin Film Temperature Sensors Utilizing Layered Carbon Nanotube Films
Previous Article in Special Issue
A Fusion Method for Combining Low-Cost IMU/Magnetometer Outputs for Use in Applications on Mobile Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Maximum Correntropy Based Unscented Particle Filter for Cooperative Navigation with Heavy-Tailed Measurement Noises

Department of Automation, Harbin Engineering University, Nantong Road No.145, Harbin 150001, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(10), 3183; https://doi.org/10.3390/s18103183
Submission received: 4 August 2018 / Revised: 17 September 2018 / Accepted: 18 September 2018 / Published: 20 September 2018
(This article belongs to the Special Issue Sensor Fusion and Novel Technologies in Positioning and Navigation)

Abstract

:
In this paper, a novel robust particle filter is proposed to address the measurement outliers occurring in the multiple autonomous underwater vehicles (AUVs) based cooperative navigation (CN). As compared with the classic unscented particle filter (UPF) based on Gaussian assumption of measurement noise, the proposed robust particle filter based on the maximum correntropy criterion (MCC) exhibits better robustness against heavy-tailed measurement noises that are often induced by measurement outliers in CN systems. Furthermore, the proposed robust particle filter is computationally much more efficient than existing robust UPF due to the use of a Kullback-Leibler distance-resampling to adjust the number of particles online. Experimental results based on actual lake trial show that the proposed maximum correntropy based unscented particle filter (MCUPF) has better estimation accuracy than existing state-of-the-art robust filters for CN systems with heavy-tailed measurement noises, and the proposed MCUPF has lower computational complexity than existing robust particle filters.

1. Introduction

Accurate navigation and localization of autonomous underwater vehicles (AUVs) are paramount for AUV autonomy. However, due to the existence of rapid attenuation in electro-magnetic, radio-frequency, and global position system (GPS) signals underwater, localization of AUVs has always been an intractable problem [1,2,3]. Inertial navigation system (INS) equipped high-accuracy inertial sensors is always employed in AUV navigation. However, the error of AUVs equipped with INS is accumulated over time and causes the localization error increase unboundedly [4,5]. With the aid of GPS signal, which can be obtained close to surface, the localization error can be alleviated, but it’s hard to realize in deep-water. Another method is to use acoustic baseline in AUVs, whose localization error is bounded [6,7,8]. Unfortunately, there is a high cost in equipment and limited working range for using static beacons. In this case, cooperative navigation (CN) based on acoustic range observation has been proposed [6,9,10,11]. High accuracy inertial sensors are installed in a small amount of AUVs in a CN system, which can assist the other AUVs that equipped with low precision dead-reckoning (DR) by using acoustic modems [5,11,12]. The CN of AUVs enables navigation information to be shared among each AUV, which makes the using scope more flexible and improves work efficiency of AUVs. For cooperative localization of AUVs, position estimation is of great importance and filtering technique based on state-space model has always been employed to solve this problem [4,9].
Particle filter (PF) is a popular method to address the state estimation problem in CN systems, in which the posterior probability density function (PDF) is approximated by a set of weighted particles based on a sequential Monte Carlo (MC) method [13,14,15]. A key component for PF is to select an appropriate importance PDF. In the conventional PF, the importance PDF is usually chosen as the prior state transition PDF for easy sampling and weight calculation [14,16]. However, state transition PDF is not the best choice for importance PDF since it does not include the latest measurement information. To solve this problem, an extended particle filter (EPF) has been proposed, where the extended Kalman filter (EKF) is employed to generate importance PDF [16,17,18]. The EPF contains not only prior information but also the latest measurement information so that it can match the optimal importance PDF (that is the true posterior PDF) better than conventional PF [14,17,18]. Based on a similar approach, an unscented particle filter (UPF) has also been proposed by utilizing the unscented Kalman filter (UKF) to generate the importance PDF. As compared with the EPF, the UPF can achieve better estimation accuracy [8,17,19,20].
Both EKF and UKF require the Gaussian assumption of measurement noise. However, there may be outliers of velocity and range in practical cooperative localization of AUVs, which results in a heavy-tailed non-Gaussian measurement noise [4,6,21]. Physically misalignment of body framework in Doppler velocity log (DVL) or water lock may induce outlier measurements of velocity [22]. Multiple acoustic propagation paths between receivers and source, which are caused by the reflection of sound wave that are caused by the changes of sound speed with depth and reflections off the sea bed and surface, may induce the outlier in range measurements [23]. In fact, such outliers lead to a bias from Gaussian distribution and can be modeled by a distribution that has heavier tails [21]. Outlier corrupted noises are generated according to:
v ~ { N ( 0 , R ) w . p .   ( 1 ε ) N ( 0 , 10 R ) w . p .   ε
where N ( β , Σ ) represents the Gaussian distribution with mean value β and covariance Σ ; ε is a perturbing parameter of outliers contamination. The PDF of Gaussian noise v 1 ( v 1 ~ N ( 0 , R ) , R = 10 ), heavy-tailed non-Gaussian noises v 2 (generated based on Equation (1) with ε = 0.1 ) and heavy-tailed non-Gaussian noises v 3 (generated based on Equation (1) with ε = 0.3 ) are shown in Figure 1. It can be seen that heavy-tailed non-Gaussian noise v 2 has more heavy tail than v 3 because that v 2 contains more outliers.
As a result, the estimation performance of EPF and UPF will degrade severely for such heavy-tailed non-Gaussian measurement noises. To address heavy-tailed measurement noises, many methods have been proposed, such as the Huber-based nonlinear Kalman filters (HNKFs) [6,16,24,25,26,27,28]. The HNKF is a kind of commonly used robust filter but its influence function doesn’t redescend, which results in deteriorating estimation performance [29]. By approximating the posterior density function (PDF) as a Student’s t distribution, a robust Student’s t based nonlinear filter (RSTNF) has been proposed [30,31,32,33]. However, such Student’s t approximation may be unreasonable in some engineering applications with slightly heavy-tailed measurement noises, which may result in deteriorating filtering performance [30]. To solve this problem, a novel robust Student’s t based Kalman filter (RSTKF) has been proposed based on variational-Bayesian (VB) method [26]. In the RSTKF, the posterior PDF is approximated as a Gaussian distribution. The RSTKF has good estimation accuracy but suffers from intensive computation. The correntropy, which is a local similarity measure in information theoretic, has been gaining more and more attention. The correntropy can capture higher-order statistical information of data directly not only the usual second-order statistical information, which has the potential to achieve better estimation performance [31,32]. Based on the maximum correntropy criterion (MCC), several robust filters including maximum correntropy Kalman filter (MCKF) and maximum correntropy unscented Kalman filter (MCUKF) have been proposed to suppress the impulsive noises [31,32,33,34,35].
In this paper, to further improve the positioning accuracy of UPF under non-Gaussian heavy-tailed measurement noise in CN of AUVs, a new maximum correntropy based unscented particle filter (MCUPF) is proposed, which modifies the update process of importance sampling of UPF. The proposed MCUPF exhibits good robustness to deviations from Gaussian distribution and has a recursive structure. Also, the Kullback-Leibler Divergence (KLD)-resampling is employed to further improve the computational efficiency of the proposed algorithm, in which the KLD is used to measure the approximate error of distribution represented by weighted particles to adjust the number of particle adaptively [36,37].
The remainder of this paper is organized as follows. In Section 2, the model of cooperative navigation and UPF algorithm are briefly reviewed. In Section 3, a new MCUPF is proposed. In Section 4, the proposed MCUPF and existing state-of-the-art robust filters are compared through an actual lake trial of AUVs, and simulation results are obtained. The conclusion is drawn in Section 5.

2. Problem Formulation

2.1. System Model

Now, we consider the model of cooperative navigation, which is based on a master-slave mode of AUVs. As the depth of AUVs can be measured accurately by pressure sensors, the three dimensional (3-D) model of CN can be transformed into two dimensional (2-D) model. The state space model of CN system based on acoustic range are described as follows [10]:
{ x k = x k 1 + Δ t ( c ^ k cos θ ^ k + s ^ k sin θ ^ k ) + w x , k 1 y k = y k 1 + Δ t ( c ^ k sin θ ^ k s ^ k cos θ ^ k ) + w y , k 1
z k = ( x k x k r ) 2 + ( y k y k r ) 2 + v k
where x k and y k are east and north positions of AUVs at time k , respectively; x k r and y k r are east and north positions of communication and navigation aid (CNA) at time k provided periodically by the acoustic modern, respectively; Δ t is the sampling time, s ^ k and c ^ k are starboard and forward velocities at time k , respectively, which are provided by DVL; θ ^ k is heading angle measured by compass.
Based on Equations (2) and (3), the state-space model of CN system can be formulated as follows:
{ x k = F x k 1 + μ k + w k 1 z k = h k ( x k ) + v k
where x k = [ x k , y k , θ k ] T denotes state vector, F = E 2 is state transition matrix and E 2 represents the 2-D identity matrix, w k 1 = [ w x , k 1 , w y , k 1 ] T denotes process noise vector, μ k denotes the control output and z k denotes relative distance between CNA and AUV. The process noise w k 1 and measurement noise v k are mutually uncorrelated Gaussian white noise processes. Here, w k 1 ~ N ( 0 , Q k ) and v k 1 ~ N ( 0 , R k ) , where Q k 1 denotes process noise covariance matrix and R k denotes measurement noise covariance. Besides, there may be correlations between process noise w x , k 1 and w y , k 1 because both of them rely on the measurement errors of compass and DVL.

2.2. Review of the Standard UPF Algorithm

Before deriving the proposed MCUPF, we first briefly introduce the standard UPF. As compared with conventional PF, UPF modifies the importance PDF through the UKF to include the latest observation information, which requires to assume that w k 1 and v k have Gaussian distributions. The UPF algorithm is summarized as follows [14,19,20].
Initialization: Inputs x ^ 0 | 0 ( j ) , P ^ 0 | 0 ( j ) and the number of particles N. Draw N particles x 0 | 0 ( j )   ( j = 1 ,   2 ,     ,   N ) from the known prior distribution p ( x 0 ) , with:
x ^ 0 | 0 ( j ) = E [ x 0 | 0 ( j ) ]
P ^ 0 | 0 ( j ) = E [ ( x 0 | 0 ( j ) x ^ 0 | 0 ) ( x 0 | 0 ( j ) x ^ 0 | 0 ) T ]
and the weights of all particles are set as f 0 ( j ) = 1 N .
For k = 1 : T
Importance sampling through UKF:
For j = 1 : N
1. Time Update
Calculate the sigma points:
X i , k 1 | k 1 = { x ^ k 1 | k 1 ( j ) , d i = κ λ + κ ,   i = 0 x ^ k 1 | k 1 ( j ) ± ( n + λ ) P ^ k 1 | k 1 ( j ) , d i = 1 2 ( λ + κ ) ,   i = 1   ,   2 ,     ,   2 n + 1
where d i are the weights of sigma points and n is the dimension of state vector; the free parameter of UKF is set as λ = 3 n .
Calculate the propagated sigma points:
X i ,   k | k 1 ( j ) = F X i , k 1 | k 1
Calculate the one-step prediction:
x ^ k | k 1 ( j ) = i = 0 2 n + 1 d i X i ,   k | k 1 ( j )
Calculate the one-step prediction error covariance matrix:
P ^ k | k 1 ( j ) = i = 0 2 n + 1 d i ( X i ,   k | k 1 ( j ) * x ^ k | k 1 ( j ) ) ( X i ,   k | k 1 ( j ) * x ^ k | k 1 ( j ) ) T + Q k 1   
2. Measurement Update
Calculate the sigma points:
X i , k | k 1 = { x ^ k | k 1 ( j ) , d i = κ n + κ i = 0 x ^ k | k 1 ( j ) ± ( n + κ ) P ^ k | k 1 ( j ) , d i = 1 2 ( n + κ ) i = 1   ,   2 ,     ,   2 n + 1
Calculate the propagated sigma points ( i = 1   ,   2 ,     ,   2 n + 1 ) :
Z i ,   k | k 1 ( j ) * = h ( X i ,   k | k 1 )
Calculate the predicted measurement vector ( i = 1   ,   2 ,     ,   2 n + 1 ) :
z ^ k | k 1 ( j ) = i = 0 2 n + 1 d i Z i ,   k | k 1 ( j ) *
Calculate the innovation covariance matrix, cross-covariance matrix and Kalman gain ( i = 1   ,   2 ,     ,   2 n + 1 ) , respectively.
P z z ,   k | k 1 ( j ) = i = 0 2 n + 1 d i ( Z i ,   k | k 1 ( j ) * z ^ k | k 1 ( j ) ) ( Z i ,   k | k 1 ( j ) * z ^ k | k 1 ( j ) ) T + R k
P x z ,   k | k 1 ( j ) = i = 0 2 n + 1 d i ( X i ,   k | k 1 ( j ) * x ^ k | k 1 ( j ) ) ( Z i ,   k | k 1 ( j ) * z ^ k | k 1 ( j ) ) T
W k ( j ) = P x z ,   k | k 1 ( j ) ( P z z ,   k | k 1 ( j ) ) 1
Calculate the state estimation of the j-th particle:
x ^ k | k ( j ) = x ^ k | k 1 ( j ) + W k ( j ) ( z k z ^ k | k 1 ( j ) )
Calculate the estimation error covariance matrix of the j-th particle:
P ^ k | k ( j ) = P ^ k | k 1 ( j ) W k ( j ) P z z , k | k 1 ( j ) ( W k ( j ) ) T
3. Sample from the importance PDF:
x k | k ( j ) ~ q ( x k | X 0 : k 1 j , Z 1 : k ) = N ( x ^ k | k ( j ) , P ^ k | k ( j ) )
4. Calculate weight of each particle:
f k ( j ) = f k 1 ( j ) p ( z k | x k ( j ) ) p ( x k ( j ) | x k 1 ( j ) ) q ( x k ( j ) | x k 1 ( j ) , z k )
5. Normalize the weights:
f k * ( j ) = f k ( j ) j = 1 N f k ( j )
Resampling: Multiply particles with high weights and suppress particles with low weights to generate a new particle set x ˜ k | k ( j )   ( j = 1 ,   2 ,     ,   N ) , which distributes according to P ( x 0 : k | z 1 : k ) approximately.
[ { x ˜ k | k ( j ) , f ˜ k ( j ) = 1 N } j = 1 N ] = Resampling [ { x k | k ( j ) , f k * ( j ) } j = 1 N ]
State Estimation:
x ^ k | k = j = 1 N f ˜ k ( j ) x ˜ k | k ( j )
P ^ k | k = j = 1 N f ˜ k ( j ) ( x ˜ k | k ( j ) x ^ k | k ) ( x ˜ k | k ( j ) x ^ k | k ) T
The standard UPF is summarized in Table 1.

3. Derivation of MCUPF

Although the UPF utilizes the latest measurement information, it can’t resist the measurement outliers due to the Gaussian assumption of measurement noise. Next, MCC is firstly used to robustify the UPF in importance sampling process, and then the KLD-resampling method is used to adjust the particle number online without decreasing the estimation accuracy.

3.1. Brief Introduction of MCC

Correntropy is a measure of similarity between two random variables X and Y , which is described as:
V ( X , Y ) = E [ κ ( X , Y ) ] = κ ( x , y ) d F X Y ( x , y )
where E ( ) is the expectation operator, κ ( ) is a shift-invariant Mercer kernel, and F X Y ( x , y ) denotes the joint distribution function of X and Y . In most cases, the joint distribution F X Y ( x , y ) is unknown and the only available data is a limited number of samples { x i , y i | i = 1 , , N } from F X Y ( x , y ) . Thus, based on the known data, the correntropy can be evaluated as follows:
V ^ ( X , Y ) = 1 N i = 1 N κ ( x i , y i ) .
The most common used Gaussian kernel is defined as:
κ ( x i , y i ) = G σ ( e i ) = exp ( ( e i ) 2 2 σ 2 ) .
where e i = x i y i and σ > 0 is the kernel bandwidth. In this case, the cost function of MCC can be written as:
J M C C = 1 N i = 1 N G σ ( e i ) .
Assume that the vector B ^ is the optimal solution based on MCC, the optimal solution can be obtained by maximizing the correntropy between desired signal y i and filter output x i :
B ^ = arg   max B   Ω 1 N i = 1 N G σ ( e i )
where Ω represents the feasible solution set of B ^ .
That is the maximum correntropy criterion, which can be incorporated in UPF to enhance robustness in non-Gaussian measurement noise.

3.2. Robustify the UPF

Considering that the correntropy is robust to outliers and can suppress negative impact of impulsive noises, MCC is incorporated into UKF to improve importance sampling in UPF.
Firstly, for measurement equation:
z k = h k ( x k ) + v k
decomposing the measurement noise covariance R k :
R k = a 2
Exploiting (31), Equation (30) can be further written as:
C k = g k ( x k ) + ζ k
with:
C k = 1 a · Z k
g k ( x k ) = 1 a · h k ( x k )
ζ k = 1 a · v k .
The main idea of MCC is to maximize the correntropy between C k and g k ( x k ) to obtain a minimized error ζ k ,   min , i.e.:
x ^ k | k = arg   min x k ( x k x ^ k | k 1 V k 2 + ( G σ ( 0 ) G σ ( ξ k , i ) ) )
where · 2 represents l 2 -norm of vector, V k denotes the inverse matrix of one-step predicted covariance matrix P k | k 1 and ξ k , i denotes the i-th component of ( C k g k ( x k ) ) . G σ ( ) is Gaussian kernel which is defined as follows [34]:
G σ ( e i ) = exp ( e i 2 2 σ 2 ) .
Considering only one fixed point iteration, the Equation (36) can be solved by using the modified measurement covariance matrix R ˜ k in the update process as [34,35]:
R ˜ k = R k / C k .
Therefore, based on MCC, robust UPF can be obtained by modifying update process of importance sampling of UPF. The proposed robust UPF is summarized in Table 2.
Remark 1.
The kernel bandwidth σ is an important parameter of robust UPF, which has significant effects on the estimation performance. In general, the algorithm will exhibit better robustness to outliers with a small σ . However, if σ is too small, it results in degraded performance and even diverge [33]. Furthermore, robust UPF is similar to UPF for a large σ , and particularly robust UPF will reduce to classic UPF when σ .

3.3. Modified Resampling Process

Based on the strong laws of large numbers, better estimation accuracy can be obtained when more particles are adopted [14]. Unfortunately, the greater the number of particles, the heavier the computational cost that is required. Therefore, it’s necessary to investigate some methods to improve the computational efficiency of PF. In the proposed MCUPF, the number of particles is adaptively adjusted online based on the KLD-resampling method, where KLD is used to determine the number of particles required for next iteration.
To ensure that the KLD between the posterior distribution before resampling and the posterior distribution after resampling is less than and equal to a pre-given error, the required number N K L D of particles can be approximated as follows [36]:
N K L D = k 1 2 ε { 1 2 9 ( k 1 ) + 2 9 ( k 1 ) z 1 δ } 3
where k is the number of bins with support and z 1 δ is the upper quartile of the standard Gaussian distribution. For a specific value of δ , the corresponding value of z 1 δ can be obtained from the normal distribution table.
KLD-resampling is composed of two parts as represented in Figure 2. Firstly, reselect particles based on their weights one by one until the required number N K L D is satisfied. Then, N K L D and bin size k are updated when a new sample is resampled [36].

3.4. The Proposed MCUPF

The proposed robust UPF is highly inefficient as a great deal of particles used for estimation. For making up the defect, the proposed MCUPF replace the resampling process with KLD-resampling, which can adjust the number of particles over time to determine the minimum amount of particles required to guarantee the estimation quality and reduce computing effort. The MCUPF, which is obtained based on MCC and KLD-resampling, is summarized in Table 3.
It can be seen from Table 3 that the importance PDF of the proposed MCUPF is generated by the robust UKF, which includes the latest observation information and can resist to measurement outliers. In addition, according to the KLD-resampling, the sample number of the proposed MCUPF can be adaptively adjusted in real time to improve the computational efficiency.

4. Lake-Water Filed Trial

The effectiveness of the proposed MCUPF is verified by employing the post processed date collected in an actual lake trial of AUVs. The filed trial was conducted in September 2014 in Lake Thai whose depth in most part varies seven to sixteen feet. In this test, two survey vessels acted as surface leaders known as CNAs, and only a single vessel surveyed as slaver AUV. All three vessels were equipped with an underwater acoustic modem ATM-885, and vessels can broadcast information to each other. For providing accurate position and time information to slaver AUV, the GPS OEMV-2RT-2 were installed in two CNAs. In addition, the DVL DS-99 for acquiring velocity information and a magnetic compass H/H HZ001 for obtaining a heading were also equipped on AUV. Figure 3 represents sensors and computer employed in this test and performance parameters and size of sensors used in this test are listed in Table 4 and Table 5, respectively. In this lake trial, acoustic data packets were sent from two leader CNAs to slaver AUV every 5 s. Range information and velocity information observed in this trial are presented in Figure 4 (velocity 1 and velocity 2 represent starboard velocity and forward velocity, respectively).
To further illustrate effect of outlier measurement on measurement noise, the measurement noise are obtained as follows:
v ^ k = z k ( x ^ k x k r ) 2 + ( y ^ k y k r ) 2
where ( x ^ k , y ^ k ) is the position of AUV at time k provided by GPS and v ^ k is approximate measurement noise value. The approximate measurement noise values can be obtained based on Equation (44) and PDF of measurement noise are shown in Figure 5. It can be seen from Figure 5 that the measurement noise have heavy-tailed distribution and Gaussian distribution can’t fit to the measurement noise value. Therefore, it’s essential to apply a robust filter algorithm to CN of AUVs.
The true north position and east position of two leaders and a slave AUV are presented in Figure 6, which are obtained using the information provided by GPSs installed on these three vessels. It’s clear to see from Figure 6 that AUV was between two leaders, and the observability of AUV can be improved by this way.
The performance of each filter is measured by position error (PE) and averaged position error (APE):
PE ( k ) = ( x ^ k x ^ k | k ) 2 + ( y ^ k y ^ k | k ) 2
APE = 1 T k = 1 T ( x ^ k x ^ k | k ) 2 + ( y ^ k y ^ k | k ) 2
where T = 1760 s is the total experimental time, ( x ^ k ,   y ^ k ) , ( x ^ k | k , y ^ k | k ) are reference position provided by GPS and estimated position at time k , respectively. The proposed MCUPF and existing state-of-the-art filters are coded with MATLAB and the simulations are run on a computer with Intel Pentium CPU at 2.90 GHz with 2.00 GB memory.

4.1. Comparisons of Different Filtering Methods

To save computation cost, the bins in KLD-resampling is usually divided in new state space which has a lower dimension than the origin state vector. The bin size is selected in terms of process noise covariance [36]. Following the advice of [36], bins in KLD-resampling of the proposed MCUPF is divided in 2-D (bin size is [ 10 , 10 ] ). The process noise covariance matrix is Q k = d i a g ( [ ( 10 m ) 2 ,   ( 10 m ) 2 ] ) and measurement noise covariance is R = 0.1 .
The existing PF [14], UPF [14], Huber unscented particle filter (HRUPF) [24], MCUKF [34], iterated VB based cubature Kalman filter (IVBCKF) [21], and the proposed MCUPF are compared, where the simple random resampling is employed in PF and UPF. In order to obtain a better performance of MCUKF, the kernel bandwidth in MCUKF and MCUPF are set as σ = 15 following the advice of [34,35]. Besides, the turning parameter and iterated times in HRUPF are commonly used γ = 1.345 and m = 3 [24,27], respectively. In IVBCKF, conjugate prior distributions for scale matrices of prior state and likelihood are inverse Wishart distribution. u k and t k denote degrees of freedom (dof) parameters of these two different scale matrices, respectively. Besides, dof parameters of prior state and likelihood all are Gamma distribution in IVBCKF. χ k and ϕ k denote shape parameter and rate parameter of dof parameters of prior state, respectively. ο k and ς k denote shape parameter and rate parameter of dof parameters of likelihood, respectively [21]. Based on reference [21], the prior parameters are set as χ k = ο k = 5 , ϕ k = ς k = 1 , u k = 8 , t k = 5 . In order to obtain a balance between computational burden and accuracy, the starting particle number are set as N s t a r t = 5000 . Besides, the maximum particle number must be less than starting number in MCUPF and the maximum particle number is set as N m a x = 2000 . Following the advice of [4,24,35,36], the parameters of above filters are presented in Table 6.
Experiment results are given as follows: Trajectory and positions taken by different filters are provided in Figure 7 and Figure 8, respectively, and APEs of different filters are shown in Figure 9, whose results and time consuming in a single step run are shown in Table 7.
It can be seen from Figure 9 and Table 7 that all filters exhibit robustness in the presence of measurement outliers. Particularly, the proposed MCUPF has the best performance among these filters. Moreover, due to the use of the KLD-resampling method, the proposed MCUPF has lower computational cost as compared with HRUPF.

4.2. Computational Complexity Analysis and Compares

In this subsection, the computational complexity of CKF, MCUKF, PF, UPF, HRUPF, IVBCKF, and the proposed MCUPF in terms of the floating point operations are analyzed. Based on [38], the total floating point operations of these filters are as follows:
S C K F = O ( 2 n 3 + m 3 ) + 2 n ( M f + M h ) + 6 n 3 + 10 n m 2 + 8 m n 2 + 6 n 2 + 2 m 2 + n m + m
S M C U K F = O ( 6 n 3 + m 3 ) + ( 2 n + 1 ) ( M f + M h ) + 6 n 3 + 4 n m 2 + 12 m n 2 + 20 m n + 20 n 2 + 4 m 2 + 5 m + 8 n
S P F = O ( n 3 + N p m 3 ) + ( M f + M h ) + N p ( 2 n 2 + m ) + S s i m p l e
S U P F = N p S H R U K F + N p [ O ( 3 n 3 + 2 m 3 ) + ( M f + M h ) + 8 n 2 + 2 m 2 + 2 m + 3 n 3 ] + S s a m p l e
S H R U P F = N p S H R U K F + N p [ O ( 3 n 3 + 2 m 3 ) + ( M f + M h ) + 8 n 2 + 2 m 2 + 2 m + 3 n 3 ] + S s a m p l e
S I V B C K F = O ( n 3 + T i 1 n 3 + m 3 + 3 T i 1 m 3 ) + [ M f + T i 1 ( 2 n + 1 ) M h ] + 3 n 2 + 12 T i 1 n m 2 + 8 T i 1 m n 2 + 4 T i 1 n 2 + 6 T i 1 m 2 + T i 1 m
S M C U P F = N p S M C U K F + N p [ O ( 3 n 3 + 2 m 3 ) + ( M f + M h ) + 8 n 2 + 2 m 2 + 2 m + 3 n 3 ] + S K L D
where:
S U K F = O ( 2 n 3 + m 3 ) + ( 2 n + 1 ) ( M f + M h ) + 6 n 3 + 10 n m 2 + 8 m n 2 +   13 n 2 + 3 m 2 + 7 n m + 4 m + 3 n
S H R U K F = O ( 2 n 3 + 2 T i 2 n 3 + m 3 ) + ( 2 n + 1 ) ( M f + M h ) + 6 n 3 + 4 n m 2 + 14 m n 2 + ( 15 + 8 T i 2 ) n 2 + 8 T i 2 m 2 + ( 7 + 16 T i 2 ) n m + ( 4 3 T i 2 ) n + ( 3 3 T i 2 ) m
where S C K F , S M C U K F , S P F , S U P F , S H R U P F , S I V B C K F , S M C U P F , S U K F , and S H R U K F denote the total floating point operations of CKF, MCUKF, PF, UPF, HRUPF, and IVBCKF, the proposed MCUPF, UKF, and Huber robust unscented Kalman filter(HRUKF), respectively; n and m denote the dimension of state vector and measurement vector, respectively; M f and M h denote the floating point operations of the computations of state equation and measurement equation, respectively; N p denotes particle number; S s i m p l e and S K L D denote the floating point operations of the simple-resampling and KLD-resampling, respectively; T i 1 and T i 2 denote iteration times of IVBCKF and HRUKF, respectively. In general, the iteration times of IVBCKF and HRUKF are more than 3 and the number of particles are tens of thousands. Thus, we can obtain the following equations based on (47)–(55):
S I V B C K F > S H R U K F > S M C U K F > S U K F > S C K F
S H R U P F > S M C U P F > S U P F > S P F > S I V B C K F > S M C U K F > S C K F .
It can be seen from (56) and (57) that four kind of PFs have higher computation complexity than others for using a large number of particles. Besides, the computation complexity of HRUPF is higher than the proposed MCUPF, due to the multi-iteration in HRUKF and the computation complexity of HRUKF is higher than the MCUKF.

5. Conclusions

In this paper, a novel MCUPF is proposed for CN systems of AUVs with heavy-tailed measurement noises, where a MCC based robust UKF is used to generate the importance PDF and KLD measure is utilized to adjust the number of particles in the resampling process. Experimental results demonstrate that the proposed MCUPF has better estimation accuracy than existing state-of-the-art robust filters and lower computational complexity than existing improved PFs.

Author Contributions

Y.F. proposed the main idea and finished the draft manuscript; Y.Z. conceived the experiment; G.W. analyzed the data; X.W. conducted the simulations; N.L. drew the figures and tables.

Funding

This research was funded by the National Natural Science Foundation of China under Grant No. 61773133, and the Natural Science Foundation of Heilongjiang Province under Grant No. F2016008.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vaulin, Y.V.; Dubrovin, F.S.; Scherbatyuk, A.F. Some Algorithms for Determining an Unknown Initial Position of AUV Using Information from a Single Beacon Navigation System. Gyroscopy Navig. 2017, 8, 209–216. [Google Scholar] [CrossRef]
  2. Scherbatyuk, A.F.; Dubrovin, F.S. About Accuracy Estimation of AUV Single-Beacon Mobile Navigation Using ASV, Equipped with DGPS. In Proceedings of the OCEANS MTS/IEEE Conference, Shanghai, China, 10–13 April 2016; pp. 1–4. [Google Scholar]
  3. Kebkal, K.G.; Mashoshin, A.I. AUV acoustic positioning methods. Gyroscopy Navig. 2017, 8, 80–89. [Google Scholar] [CrossRef]
  4. Huang, Y.L.; Zhang, Y.G.; Xu, B.; Wu, Z.M.; Chambers, J.A. A New Adaptive Extended Kalman Filter for Cooperative Localization. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 353–368. [Google Scholar] [CrossRef]
  5. Tan, H.P.; Diamant, R.; Seah, W.K.G.; Waldmeyer, M. A Survey of Techniques and Challenges in Underwater Localization. Ocean Eng. 2011, 38, 1663–1676. [Google Scholar] [CrossRef]
  6. Gao, W.; Liu, Y.L.; Xu, B. Robust Huber-based iterated divided difference filtering with application to cooperative localization of autonomous underwater vehicles. Sensors 2014, 4, 24523–24542. [Google Scholar] [CrossRef] [PubMed]
  7. Paull, L.; Saeedi, S.; Seto, M.; Li, H. AUV Navigation and Localization: A Review. IEEE J. Ocean. Eng. 2014, 39, 131–149. [Google Scholar] [CrossRef]
  8. Allotta, B.; Caiti, A.; Costanzi, R.; Fanelli, F.; Meli, E.; Ridolfi, A. Development and Online Validation of an Ukf-based Navigation Algorithm for AUVs. IFAC-PapersOnline 2016, 49, 69–74. [Google Scholar] [CrossRef]
  9. Bahr, A.; Leonard, J.J.; Fallon, M.F. Cooperative localization for autonomous underwater vehicles. Int. J. Robot. Res. 2009, 28, 714–728. [Google Scholar] [CrossRef]
  10. Fallon, M.F.; Papadopoulos, G.; Leonard, J.J. Cooperative AUV Navigation Using a Single Surface Craft; Massachusetts Institute of Technology: Cambridge, MA, USA, 2010; pp. 331–340. [Google Scholar]
  11. Papadopoulos, G.; Fallon, M.F.; Leonard, J.J.; Patrikalakis, N.M. Cooperative Localization of Marine Vehicles Using Nonlinear State Estimation. In Proceedings of the 2010 IEEE/RSJ International Conference on Intelligent Robots and Systems, Taipei, Taiwan, 18–22 October 2010; pp. 4874–4879. [Google Scholar]
  12. Paull, L.; Saeedi, S.; Seto, M.; Li, H. Sensor-driven Online Coverage Planning for Autonomous Underwater Vehicles. IEEE/ASME Trans. Mechatron. 2013, 18, 1827–1838. [Google Scholar] [CrossRef]
  13. Maurelli, F.; Krupinski, S.; Petillot, Y.; Salvi, J. A Particle Filter Approach for AUV Localization. Oceans 2008, 15, 1–7. [Google Scholar]
  14. Simon, D. Optimal State Estimation: Kalman, H and Nonlinear Approaches, 1st ed.; Wiley: Hoboken, NJ, USA, 2006. [Google Scholar]
  15. Bravo, F.G.; Vale, A.; Ribeiro, M.I.; Salvi, J. Navigation Strategies for Cooperative Localization Based on a Particle-filter Approach. Integr. Comput. Aided Eng. 2007, 14, 263–279. [Google Scholar] [CrossRef]
  16. Dan, C.; Arnaud, D. A survey of convergence results on particle filtering methods for practitioners. IEEE Trans. Signal Process. 2002, 50, 736–746. [Google Scholar] [Green Version]
  17. Aggarwal, P.; Gu, D.Q.; Nassar, S.; Syed, Z.; EI-sheimy, N. Extended Particle Filter (EPF) for INS/GPS Land Vehicle Navigation Applications. In Proceedings of the Institute of Navigation Satellite Division Technical Meeting, Fort Worth, TX, USA, 25–28 September 2007; pp. 2619–2626. [Google Scholar]
  18. Li, L.Q.; Ji, H.B.; Luo, J.H. The iterated extended Kalman particle filter. In Proceedings of the IEEE International Symposium on Communications and Information Technology, Beijing, China, 16 January 2005; pp. 1213–1216. [Google Scholar]
  19. Julier, S.; Uhlmann, J.; Durrant-Whyte, H.F. A new method for the nonlinear transformation of means and covariance in filters and estimations. IEEE Trans. Autom. Control 2000, 45, 477–482. [Google Scholar] [CrossRef]
  20. Carpenter, J.; Clifford, P.; Fearnhead, P. Improved particle filter for nonlinear problems. IEE Proc. Radar Sonar Navig. 1999, 146, 2–7. [Google Scholar] [CrossRef]
  21. Huang, Y.L.; Zhang, Y.G.; Xu, B.; Wu, Z.M.; Chambers, J.A. A new outlier-robust student’s t based gaussian approximate Filter for Cooperative Localization. IEEE/ASME Trans. Mechatron. 2017, 22, 2308–2386. [Google Scholar] [CrossRef]
  22. Vasilijevic, A.; Borovic, B.; Vukic, Z. Underwater Vehicle Localization with Complementary Filter: Performance Analysis in the Shallow Water Environment. J. Intell. Robot. Syst. 2012, 68, 373–386. [Google Scholar] [CrossRef]
  23. Vaganay, J.; Leonard, J.J.; Bellingham, J.G. Outlier Rejection for Autonomous Acoustic Navigation. In Proceedings of the IEEE International Conference on Robotics and Automation, Cambridge, MA, USA, 22–28 April 1996; pp. 2174–2181. [Google Scholar]
  24. Guo, R.X.; Gan, Q.; Zhang, J.W.; Guo, K.; Dong, J.K. Huber cubature particle filter and online state estimation. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2017, 231, 158–167. [Google Scholar] [CrossRef]
  25. Li, K.L.; Hu, B.Q.; Chang, L.B.; Li, Y. Robust square-root cubature Kalman filter based on Huber’s M-estimation methodology. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2015, 229, 1236–1245. [Google Scholar] [CrossRef]
  26. Tseng, C.H.; Lin, S.F.; Jwo, D.J. Robust Huber-based cubature Kalman filter for GPS navigation processing. J. Navig. 2016, 70, 527–546. [Google Scholar] [CrossRef]
  27. Chang, L.B.; Chang, G.B.; Li, A. Huber-based novel robust unscented Kalman filter. IET Sci. Meas. Technol. 2012, 6, 502–509. [Google Scholar] [CrossRef]
  28. Gandhi, M.A.; Mili, L. Robust Kalman Filter Based on a Generalized Maximum-Likelihood-Type Estimator. IEEE Trans. Signal Process. 2010, 58, 2509–2520. [Google Scholar] [CrossRef]
  29. Huang, Y.L.; Zhang, Y.G.; Li, N.; Chambers, J.A. Robust Student’s t Based Nonlinear Filter and Smoother. IEEE Trans. Aerosp. Electron. Syst. 2017, 52, 2586–2596. [Google Scholar] [CrossRef]
  30. Huang, Y.L.; Zhang, Y.G.; Wu, Z.M.; Li, N.; Chambers, J.A. A novel adaptive Kalman filter with inaccurate process and measurement noise covariance matrices. IEEE Trans. Autom. Control 2018, 63, 594–601. [Google Scholar] [CrossRef]
  31. Izanloo, R.; Fakoorian, S.A.; Yazdi, H.S.; Simon, D. Kalman filtering based on the maximum correntropy criterion in the presence of non-Gaussian noise. In Proceedings of the 2016 Annual Conference on Information Science and Systems (CISS), Princeton, NJ, USA, 16–18 March 2016; pp. 500–505. [Google Scholar]
  32. Liu, X.; Chen, B.D.; Xu, B.; Wu, Z.Z.; Honeine, P. Maximum correntropy unscented filter. Int. J. Syst. Sci. 2017, 48, 1607–1615. [Google Scholar] [CrossRef]
  33. Chen, B.D.; Wang, J.J.; Zhao, H.Q.; Zheng, N.N.; Principe, J.C. Convergence of a fixed-point algorithm under maximum correntropy criterion. IEEE Signal Process. Lett. 2015, 10, 1723–1727. [Google Scholar] [CrossRef]
  34. Liu, X.; Qu, H.; Zhao, J.H.; Yue, P.C.; Wang, M. Maximum Correntropy Unscented Kalman Filter for Spacecraft Relative State Estimation. Sensors 2016, 16, 1530. [Google Scholar] [CrossRef] [PubMed]
  35. Hou, B.W.; He, Z.M.; Li, D.; Zhou, H.Y.; Wang, J.Q. Maximum Correntropy Unscented Kalman Filter for Ballistic Missile Navigation System based on SINS/CNS Deeply Integrated Mode. Sensors 2018, 18, 1724. [Google Scholar] [CrossRef] [PubMed]
  36. Li, T.; Sun, S.; Sattar, T.P. Adapting sample size in particle filters through KLD-resampling. Electron. Lett. 2013, 49, 740–742. [Google Scholar] [CrossRef] [Green Version]
  37. Fox, D. Adapting the sample size in particle filters through KLD-sampling. Int. J. Robot. Res. 2003, 22, 985–1003. [Google Scholar] [CrossRef]
  38. Golub, G.H.; Van Loan, C.F. Matrix Computation, 2nd ed.; The Johns Hopkins University Press: Baltimore, MD, USA, 2013. [Google Scholar]
Figure 1. Gaussian distribution and heavy tailed non-Gaussian distribution induced by outlier interference.
Figure 1. Gaussian distribution and heavy tailed non-Gaussian distribution induced by outlier interference.
Sensors 18 03183 g001
Figure 2. Procedure of Kullback-Leibler Divergence (KLD)-resampling.
Figure 2. Procedure of Kullback-Leibler Divergence (KLD)-resampling.
Sensors 18 03183 g002
Figure 3. Sensors and computer employed in this test.
Figure 3. Sensors and computer employed in this test.
Sensors 18 03183 g003
Figure 4. (a) Range measurement; (b) Velocity measurement.
Figure 4. (a) Range measurement; (b) Velocity measurement.
Sensors 18 03183 g004
Figure 5. Probability density function (PDF) of measurement noise.
Figure 5. Probability density function (PDF) of measurement noise.
Sensors 18 03183 g005
Figure 6. True trajectory of two leaders and a slaver autonomous underwater vehicles (AUVs).
Figure 6. True trajectory of two leaders and a slaver autonomous underwater vehicles (AUVs).
Sensors 18 03183 g006
Figure 7. Paths taken by the proposed filter and existing filters.
Figure 7. Paths taken by the proposed filter and existing filters.
Sensors 18 03183 g007
Figure 8. Position of AUV taken by the proposed filter and existing filters.
Figure 8. Position of AUV taken by the proposed filter and existing filters.
Sensors 18 03183 g008
Figure 9. Position errors (PEs) of the proposed filter and existing filters in cooperative navigation (CN) of AUVs.
Figure 9. Position errors (PEs) of the proposed filter and existing filters in cooperative navigation (CN) of AUVs.
Sensors 18 03183 g009
Table 1. Algorithm of the standard unscented particle filter (UPF).
Table 1. Algorithm of the standard unscented particle filter (UPF).
1. Initialization by Equations (5) and (6).
2. Importance sampling by Equations (7)–(21):
Time update (Equations (7)–(10));
Measurement update (Equations (11)–(21)).
3. Resampling by Equation (22).
4. State estimation by Equations (23) and (24).
Table 2. Algorithm of the proposed robust UPF.
Table 2. Algorithm of the proposed robust UPF.
1. Initialization by Equations (5) and (6).
2. Importance sampling like UPF (Equations (7)–(21)) but modify the Equation (14) as follows:
P z z , k | k 1 ( j ) = i = 0 m d i ( Z i , k | k 1 ( j ) * z ^ k | k 1 ( j ) ) ( Z i , k | k 1 ( j ) * z ^ k | k 1 ( j ) ) T + R ˜ k (39)
3. Resampling by Equation (22).
4. State estimation by Equations (23) and (24).
Table 3. Algorithm of the proposed maximum correntropy based unscented particle filter (MCUPF).
Table 3. Algorithm of the proposed maximum correntropy based unscented particle filter (MCUPF).
1. Initialization by Equations (6) and (7).
2. Importance sampling as robust UPF
3. Resampling: KLD-resampling as shown in Figure 1.
[ { x ˜ k | k ( j ) , f ˜ k ( j ) = 1 N K L D } j = 1 N K L D ] = KLD resampling [ { x k | k ( j ) , f k * ( j ) } j = 1 N ] (41)
4. State estimation by N K L D resampled particles from last step.
x ^ k / k = j = 1 N K L D f ˜ k ( j ) x ˜ k / k ( j ) (42)
P ^ k | k = j = 1 N K L D f ˜ k ( j ) ( x k | k ( j ) x ^ k | k ) ( x k | k ( j ) x ^ k | k ) T (43)
Table 4. The parameters of employed sensors.
Table 4. The parameters of employed sensors.
SensorsMetricParameters
Acoustic modem: ATM-885Working range
Error rate
Up to 8 km
Less than 10 7
GPS: OEMV-2RT-2Position accuracy
Date update rate
1.8 m (CEP)
10 Hz
DVL:DS-99Working range
Measurement accuracy
−150 m/s–200 m/s
0.1–0.3%
Magnetic Compass:H/H HZ001Heading accuracy0.3°
Table 5. The size of employed sensors.
Table 5. The size of employed sensors.
SensorsSize
Acoustic modem: ATM-885 140   mm in diameter
850   mm long
GPS: OEMV-2RT-2 160 × 160 × 160   mm
DVL:DS-99(Transceiver) 350 × 300 × 30   mm
Magnetic Compass:H/H HZ001 300 × 300 × 280   mm
Table 6. Parameters of employed filters.
Table 6. Parameters of employed filters.
FiltersParameters
MCUKFKernel bandwidth σ = 15
PFThe number of particles N s t a r t = 5000
UPFThe number of particles N s t a r t = 5000
HRUPFTurning parameter γ = 1.345
Iterated times m = 3
IVBCKFPrior parameters χ k = ο k = 5 , ϕ k = ς k = 1 , u k = 8 , t k = 5
MCUPFKernel bandwidth σ = 15
Bin size [ 10 , 10 ]
Error bound ε = 0.15
Bound parameter δ = 0.01
Maximum particle number N m a x = 2000
Starting particle number N s t a r t = 5000
Table 7. Comparisons of averaged position errors (APEs) and computation time in a single step of the proposed filter and existing filters.
Table 7. Comparisons of averaged position errors (APEs) and computation time in a single step of the proposed filter and existing filters.
FiltersAPE (m)Time (s)
CKF [14]152.0 1.3644 × 10 4
MCUKF [34]18.2 1.6458 × 10 4
PF [14]14.1 2.6165 × 10 2
UPF [20]15.10.1725
HRUPF [24]12.40.3398
IVBCKF [21]12.1 7.9420 × 10 4
MCUPF8.60.2156

Share and Cite

MDPI and ACS Style

Fan, Y.; Zhang, Y.; Wang, G.; Wang, X.; Li, N. Maximum Correntropy Based Unscented Particle Filter for Cooperative Navigation with Heavy-Tailed Measurement Noises. Sensors 2018, 18, 3183. https://doi.org/10.3390/s18103183

AMA Style

Fan Y, Zhang Y, Wang G, Wang X, Li N. Maximum Correntropy Based Unscented Particle Filter for Cooperative Navigation with Heavy-Tailed Measurement Noises. Sensors. 2018; 18(10):3183. https://doi.org/10.3390/s18103183

Chicago/Turabian Style

Fan, Ying, Yonggang Zhang, Guoqing Wang, Xiaoyu Wang, and Ning Li. 2018. "Maximum Correntropy Based Unscented Particle Filter for Cooperative Navigation with Heavy-Tailed Measurement Noises" Sensors 18, no. 10: 3183. https://doi.org/10.3390/s18103183

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