Next Article in Journal
A Mobility Management Using Follow-Me Cloud-Cloudlet in Fog-Computing-Based RANs for Smart Cities
Next Article in Special Issue
Investigation of the Temperature Fluctuation of Single-Phase Fluid Based Microchannel Heat Sink
Previous Article in Journal
Off-Line Evaluation of Mobile-Centric Indoor Positioning Systems: The Experiences from the 2017 IPIN Competition
Previous Article in Special Issue
A Comprehensive Study of a Micro-Channel Heat Sink Using Integrated Thin-Film Temperature Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Sensor Optimal Data Fusion Based on the Adaptive Fading Unscented Kalman Filter

1
School of Automatics, Northwestern Polytechnical University, Xi’an 710072, China
2
School of Engineering, RMIT University, Bundoora, VIC 3083, Australia
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(2), 488; https://doi.org/10.3390/s18020488
Submission received: 4 January 2018 / Revised: 1 February 2018 / Accepted: 3 February 2018 / Published: 6 February 2018
(This article belongs to the Special Issue Integrated Sensors)

Abstract

:
This paper presents a new optimal data fusion methodology based on the adaptive fading unscented Kalman filter for multi-sensor nonlinear stochastic systems. This methodology has a two-level fusion structure: at the bottom level, an adaptive fading unscented Kalman filter based on the Mahalanobis distance is developed and serves as local filters to improve the adaptability and robustness of local state estimations against process-modeling error; at the top level, an unscented transformation-based multi-sensor optimal data fusion for the case of N local filters is established according to the principle of linear minimum variance to calculate globally optimal state estimation by fusion of local estimations. The proposed methodology effectively refrains from the influence of process-modeling error on the fusion solution, leading to improved adaptability and robustness of data fusion for multi-sensor nonlinear stochastic systems. It also achieves globally optimal fusion results based on the principle of linear minimum variance. Simulation and experimental results demonstrate the efficacy of the proposed methodology for INS/GNSS/CNS (inertial navigation system/global navigation satellite system/celestial navigation system) integrated navigation.

1. Introduction

With the rapid development of electronics technologies, various sensors have been developed and applied to many engineering fields such as integrated navigation, target tracking, signal processing, and networked communications [1,2,3,4]. Since the use of multiple sensors provides more accurate and reliable information than that of a single sensor alone, multi-sensor data fusion has received considerable attention in recent years. The basic idea of multi-sensor data fusion is to combine the data obtained from multiple sensors and further associate them with database information to achieve improved estimation accuracy compared to the use of a single sensor [5,6].
The existing studies on multi-sensor data fusion can be divided into two categories [4,6]. One is the centralized fusion method, which can achieve a globally optimal state estimation by directly combining local measurement models to form an augmented measurement model [7,8]. This method has minimal information loss. However, it causes large computation and communication burdens in the fusion center due to high-dimension computations and large data memories. Further, it has a poor robustness and reliability when one sensor is faulty [1,9]. The other is the decentralized fusion method, which can yield globally optimal or suboptimal state estimations according to certain criteria [9,10]. This method has small computation and communication burdens in the fusion center due to the decentralized structure. Furthermore, the decentralized structure also enables easy fault detection and isolation. Due to its superiority over the centralized fusion method, the decentralized fusion method has been extensively used in engineering fields.
The federated Kalman filter (FKF) is a typical example of the decentralized fusion method. It employs the principle of information sharing for local and global filters, and eliminates the correlation between local estimations by using the technique of covariance upper bound [11,12]. However, FKF is only suitable for multi-sensor linear stochastic systems [13,14]. Its fusion accuracy is sharply deteriorated when applied to multi-sensor nonlinear stochastic systems. In addition, due to the use of process noise covariance’s upper bound rather than process noise covariance itself, the resultant local state estimations are at suboptimal level, making the global state estimation also at suboptimal level. Furthermore, FKF also requires the local estimations resulting from local filters to be independent initially, which is difficult to implement in practical engineering [4,14,15].
The unscented Kalman filter (UKF) is a promising filtering method to estimate the state of a nonlinear stochastic system [16,17]. This method can approximate the posterior mean and covariance of any Gaussian random variable in third-order accuracy by using unscented transformation (UT). It has advantages of high estimation accuracy, high convergence rate and simple implementation compared to other nonlinear filtering methods [17,18]. Due to these merits, Hu and Huang [13] presented an UKF-based FKF (UKF-FKF) by combining UKF with FKF to address the data fusion problem in multi-sensor nonlinear stochastic systems. This method adopts UKFs as local filters to calculate local state estimations, avoiding the linearization error of the system model involved in FKF and improving the estimation accuracy of local filters. However, it still suffers from the problem due to the use of covariance upper bound in FKF. Further, since UKF is designed based on the condition that the nonlinear system can be exactly modeled [19,20,21], UKF-FKF also requires the system model to be accurate. If the system model of a multi-sensor nonlinear stochastic system involves modeling error, the performances of local filters in UKF-FKF will be degraded or even divergent, leading to the poor fusion accuracy.
Currently, the existing multi-sensor data fusion methods are mainly dominated by inhibiting the influence of covariance upper bound on the fusion accuracy. Chang et al. [22] studied an information matrix-based distributed fusion Kalman filter. The authors’ group also reported a random weighting estimation-based multi-sensor optimal data fusion approach and a matrix weighted multi-sensor data fusion method [23,24]. However, the above studies are mainly focused on multi-sensor linear stochastic systems. Just recently, the authors’ group reported a UKF-based multi-sensor optimal data fusion method (UKF-MODF) for multi-sensor nonlinear stochastic systems to inhibit the influence of covariance upper bound on the fusion solution [14]. This method develops an optimal data fusion scheme according to the principle of linear minimum variance to overcome the problem due to the use of covariance upper bound in FKF. It can achieve the globally optimal state estimation for multi-sensor nonlinear stochastic systems. However, since UKF still serves as local filters to calculate local state estimations, the estimation accuracy of UKF-MODF also relies on the accuracy of the system model. Furthermore, the global state estimation process of UKF-MODF is only suitable for the case with two local filters. It is necessary to provide a process of global state estimation for general cases with N local filters.
The system model of a nonlinear stochastic system consists of process and measurement models. In practical applications, the measurement model’s accuracy can be guaranteed by using high-precision measurement equipment together with plenty of available measurement data [25,26]. However, since the process model is a theoretical approximation to the real-world dynamic system [26], it inevitably involves error. In some applications with the preference of computational simplicity and real-time performance, such an approximation is intended given the complexity of the dynamic environment [25]. Therefore, the process-modeling error is generally considered to be the main source of model error for a nonlinear stochastic system [25,26,27,28].
Adaptive UKF is a strategy to improve the UKF adaptability and robustness against process-modeling error by adjusting the Kalman gain matrix for state estimation. Song and Han presented an adaptive UKF to inhibit the process-modeling error by updating the process noise covariance via minimization of the difference between computed and actual innovation covariances [29]. However, this method requires the calculation of partial derivatives, causing a relatively large computation burden. Meng et al. reported a covariance matching-based adaptive UKF to make online estimates of and adjust process noise covariance using innovation and residual sequences [30]. However, this covariance matching technique yields a steady-state estimation error, leading to a limited improvement in the filtering accuracy [20]. Compared to the above two adaptive UKFs, the adaptive fading UKF (AFUKF) can effectively inhibit the influence of process-modeling error on the filtering solution under a modest computation load [27,28]. In this method, when process-modeling error is detected, an adaptive fading factor is constructed based on innovation vector and its corresponding statistical information. By using the adaptive fading factor to scale process noise covariance or predicted state covariance, historical information is less weighted while current measurement information is more weighted to tune UKF recursively to alleviate the influence of process-modeling error on the estimation solution [19,28]. Therefore, the introduction of AFUKF into multi-sensor optimal data fusion provides a solution to improve the adaptability and robustness of UKF-MODF.
This paper focuses on improvement of our previous work (UKF-MODF) for multi-sensor nonlinear stochastic systems. It presents a novel multi-sensor optimal data fusion methodology based on AFUKF to address the UKF-MODF problems, leading to improved adaptability and robustness of data fusion for multi-sensor nonlinear stochastic systems. This methodology is in a two-level fusion structure: at the bottom level, an AFUKF based on the Mahalanobis distance is developed and serves as local filters to improve the adaptability and robustness of local state estimations against process-modeling error; at the top level, an UT-based multi-sensor optimal data fusion for the case of N local filters is established according to the principle of linear minimum variance to calculate globally optimal state estimation via fusion of local estimations. The proposed methodology not only improves the robustness of multi-sensor data fusion against process-modeling error, but it also achieves globally optimal fusion results based on the principle of linear minimum variance. Simulations and practical experiments as well as comparison analysis with FKF, UKF-FKF and UKF-MODF have been conducted to comprehensively evaluate the performance of the proposed multi-sensor optimal fusion methodology for INS/GNSS/CNS (inertial navigation system/global navigation satellite system/celestial navigation system) integrated navigation system.

2. Adaptive Fading UKF Based Multi-Sensor Optimal Data Fusion

The framework of the proposed AFUKF-based multi-sensor optimal data fusion methodology (AFUKF-MODF) is shown in Figure 1. It has a two-level fusion structure: at the bottom level, an AFUKF based on the Mahalanobis distance is developed and serves as local filters to improve the robustness of local state estimations against process-modeling error; at the top level, according to the principle of linear minimum variance, the UT-based multi-sensor optimal data fusion for the case of N local filters is established to calculate the globally optimal state estimation by fusion of local estimations.

2.1. Local State Estimation

Consider the i th ( i = 1 , 2 , , N ) local filter with the following dynamic model
{ x ( k ) = f ( x ( k 1 ) ) + w ( k ) z i ( k ) = h i ( x ( k ) ) + v i ( k )
where x ( k ) R n is the system state vector; f ( ) is the nonlinear state function; w ( k ) is the process noise which is commonly assumed as a zero-mean Gaussian white noise with covariance Q 0 ; z i ( k ) R m i is the measurement of the i th local filter; h i ( ) is the nonlinear measurement function of the i th local filter; and v i ( k ) is the measurement noise of the i th local filter and is commonly assumed as a zero-mean Gaussian white noise with covariance R i > 0 .

2.1.1. Classical UKF

In order to describe the Mahalanobis distance-based AFUKF clearly, let us briefly review the classical UKF at first. For the nonlinear system model given by (1), the procedure of the classical UKF can be summarized as:
Step 1:Initialization
{ x ^ i ( 0 ) = E [ x i ( 0 ) ] P i ( 0 ) = E [ ( x i ( 0 ) x ^ i ( 0 ) ) ( x i ( 0 ) x ^ i ( 0 ) ) T ]
Step 2:Time Update
Assume that the state estimate x ^ i ( k 1 ) and the corresponding error covariance matrix P i ( k 1 ) are given, the sigma points are selected as
{ χ j , k 1 ( i ) = x ^ i ( k 1 ) j = 0 χ j , k 1 ( i ) = x ^ i ( k 1 ) + a ( n P i ( k 1 ) ) j j = 1 , 2 , , n χ j , k 1 ( i ) = x ^ i ( k 1 ) a ( n P i ( k 1 ) ) j n j = n + 1 , n + 2 , , 2 n
where a is a tuning parameter to determine the spread of the sigma points around x ^ i ( k 1 ) and it is commonly set to a small positive value; and ( n P i ( k 1 ) ) j denotes the jth column of the square root of the matrix n P i ( k 1 ) .
Calculate the predicted state mean and covariance as
χ j , k / k 1 ( i ) = f ( χ j , k 1 ( i ) )   ( j = 0 , 1 , , 2 n )
x ^ i ( k / k 1 ) = j = 0 2 n ω j χ j , k / k 1 ( i )
P i ( k / k 1 ) = j = 0 2 n ω j ( χ j , k / k 1 ( i ) x ^ i ( k / k 1 ) ) ( χ j , k / k 1 ( i ) x ^ i ( k / k 1 ) ) T + Q
where { ω j = 1 1 a 2 , j = 0 ω j = 1 2 n a 2 , j = 1 , 2 , , 2 n .
Step 3:Sigma Point Update
Select a set of new sigma points with the mean x ^ i ( k / k 1 ) and the covariance P i ( k / k 1 )
{ χ j , k / k 1 ( i ) = x ^ i ( k / k 1 ) , j = 0 χ j , k / k 1 ( i ) = x ^ i ( k / k 1 ) + ( a n P i ( k / k 1 ) ) j , j = 1 , 2 , , n χ j , k / k 1 ( i ) = x ^ i ( k / k 1 ) ( a n P i ( k / k 1 ) ) j n , j = n + 1 , n + 2 , , 2 n
Step 4:Measurement Update
The weighted mean and covariance of the predicted measurement are computed by
γ j , k / k 1 ( i ) = h ( χ j , k / k 1 ( i ) )
z ^ i ( k / k 1 ) = j = 0 2 n ω j γ j , k / k 1 ( i ) = j = 0 2 n ω j h ( χ j , k / k 1 ( i ) )
P i ( z ^ i ( k / k 1 ) ) = j = 0 2 n ω j ( γ j , k / k 1 ( i ) z ^ i ( k / k 1 ) ) ( γ j , k / k 1 ( i ) z ^ i ( k / k 1 ) ) T + R i
Subsequently, the state estimate and associated error covariance matrix are updated as
P i ( x ^ i ( k / k 1 ) z ^ i ( k / k 1 ) ) = j = 0 2 n ω j ( χ j , k / k 1 ( i ) x ^ i ( k / k 1 ) ) ( γ j , k / k 1 ( i ) z ^ i ( k / k 1 ) ) T
K i ( k ) = P i ( x ^ i ( k / k 1 ) z ^ i ( k / k 1 ) ) P i 1 ( z ^ i ( k / k 1 ) )
x ^ i ( k ) = x ^ i ( k / k 1 ) + K i ( k ) ( z i ( k ) z ^ i ( k / k 1 ) )
P i ( k ) = P i ( k / k 1 ) K i ( k ) P i ( z ^ i ( k / k 1 ) ) K i T ( k )
Step 5:Get back to Step 2 for the next sample until all samples are processed.

2.1.2. Mahalanobis Distance Based Adaptive Fading UKF

The Mahalanobis distance of innovation vector, which is based on hypothesis testing theory, is commonly used as a measure to identify system modeling error for Gaussian systems [27,31]. This paper adopts the Mahalanobis distance into AFUKF and further develop a Mahalanobis distance-based AFUKF to improve the adaptability and robustness of the classical UKF against process-modeling error for multi-sensor nonlinear stochastic systems.
Define the innovation vector of the ith local filter as
z ˜ i ( k ) = z i ( k ) z ^ i ( k / k 1 )
For the nonlinear Gaussian system given by (1), z ˜ i ( k ) should obey the zero-mean Gaussian distribution with the covariance
P i ( z ^ i ( k / k 1 ) ) = j = 0 2 n ω j ( γ j , k / k 1 ( i ) z ^ i ( k / k 1 ) ) ( γ j , k / k 1 ( i ) z ^ i ( k / k 1 ) ) T + R i
Thus, the square of the Mahalanobis Distance of the innovation vector should obey the chi-square distribution with m i degrees of freedom [27,31], i.e.,
M i 2 ( k ) = z ˜ i T ( k ) ( P i ( z ^ i ( k / k 1 ) ) ) 1 z ˜ i ( k ) χ m i 2
where m i is the dimension of the measurement in the i th local filter.
According to the hypothesis testing theory, for a given significance level α , we have
P ( M i 2 ( k ) χ m i , α 2 ) = 1 α   ( 0 < α 1 )
where P ( ) represents the probability of a random event.
If (18) holds, this means the system works under the optimal condition, i.e., the nonlinear multi-sensor system described by (1) does not have process-modeling error. However, if (18) does not hold, it can be concluded with high probability that there exists process-modeling error in the multi-sensor system (1). In this case, AFUKF incorporates a time-varying adaptive fading factor into the predicted state covariance matrix to refrain from the influence of prior knowledge of the current state estimate. Thus, the predicted state covariance matrix in AFUKF is modified as
P i ( k / k 1 ) = λ k i ( j = 0 2 n ω j ( χ j , k / k 1 ( i ) x ^ i ( k / k 1 ) ) ( χ j , k / k 1 ( i ) x ^ i ( k / k 1 ) ) T + Q )
Introduce the adaptive fading factor λ k i to the predicted state covariance matrix such that (18) holds. Thus, we have the following equation
g i ( λ k i ) = z ˜ i T ( k ) ( P i ( z ^ i ( k / k 1 ) ) ) 1 z ˜ i ( k ) χ m i , α 2 = 0
where P i ( z ^ i ( k / k 1 ) ) is the measurement prediction covariance matrix calculated by using P i ( k / k 1 ) .
It can be seen from (20) that the determination of the adaptive fading factor λ k i is a problem of solving the nonlinear equation. Since the derivative of g i ( λ k i ) with respect to λ k i is difficult to calculate, the traditional Newton’s method cannot be used to solve the nonlinear Equation (20). In this paper, a chord secant method [32,33] is adopted to iteratively determine the adaptive fading factor λ k i by solving the nonlinear Equation (20). Thus, we have
λ k i ( t + 1 ) = λ k i ( t ) g i ( λ k i ( t ) ) [ λ k i ( t ) λ k i ( t 1 ) ] g i ( λ k i ( t ) ) g i [ λ k i ( t 1 ) ] ,   ( t = 1 , 2 , )
where t represents the iteration time.
Substituting (20) into (21) yields
λ k i ( t + 1 ) = λ k i ( t ) [ M ¯ k i ( λ k i ( t ) ) χ m i , α 2 ] [ λ k i ( t ) λ k i ( t 1 ) ] M ¯ k i ( λ k i ( t ) ) M ¯ k i ( λ k i ( t 1 ) ) ,   ( t = 1 , 2 , )
where M ¯ k i ( λ k i ( t ) ) = z ˜ i T ( k ) ( P i ( z ^ i ( k / k 1 ) ) ) 1 z ˜ i ( k ) .
The above iterative process to determine the adaptive fading factor λ k i is initialized as λ k i ( 1 ) = 1 and λ k i ( 0 ) = 0 . It will be terminated when the criterion M ¯ k i [ λ k i ( t ) ] χ m i , α 2 is satisfied. The calculation process of the Mahalanobis distance-based AFUKF can be summarized as:
Step 1:
Initialization
Initiate the Local filters by presetting the x ^ i ( 0 ) and P i ( 0 ) as (2).
Step 2:
Local Filtering
(i)
Conduct the classical UKF procedures (3)–(10) and calculate the M ¯ k i [ λ k i ( t ) ] .
(ii)
If M ¯ k i [ λ k i ( t ) ] χ m i , α 2 ,
  • Perform the classical UKF procedures (11)–(14) to compute the local state estimations.
Else,
  • Determine the adaptive fading factor λ k i by iteratively computing (22) until the criterion M ¯ k i [ λ k i ( t ) ] χ m i , α 2 is satisfied.
  • Then, replace the predicted state covariance matrix P i ( k / k 1 ) of the classical UKF with the modified type P i ( k / k 1 ) as described by (19).
  • Complete the classical UKF procedures (7)–(14) to update the local state estimations.
(iii)
Repeat Steps (i) and (ii) for the next sample.
Step 3:
Complete the local state estimations until all samples are processed.
Each local filter is performed in a parallel manner to generate the local state estimation x ^ i ( k ) and the corresponding error covariance P i ( k ) ( i = 1 , 2 , , N ) by using the Mahalanobis distance-based AFUKF. After this, these local state estimations will be further fused by the UT-based data fusion method, which is to be described in the next section, to obtain the globally optimal state estimation.

2.2. Global State Estimation

Suppose that the local state estimation achieved by the ith ( i = 1 , 2 , , N ) local filter is x ^ i ( k ) . According to the linear weighting fusion criterion, the globally optimal state estimation can be described as the following form
x ^ ( k ) = i = 1 N α i x ^ i ( k )
i = 1 N α i = I n
where α i ( i = 1 , 2 , , N ) is the weighting matrix to be determined; and I n represents an n-dimensional identity matrix.
Based on the principle of linear minimum variance, the globally optimal state estimation x ^ ( k ) should satisfy the following conditions [14,24]:
Condition 1.
x ^ ( k ) must be the unbiased estimation of x ( k ) , i.e., E { x ^ ( k ) x ( k ) } = 0 .
Condition 2.
x ^ ( k ) makes t r { P ( k ) } minimum, where P ( k ) is the error covariance matrix of x ^ ( k ) , and t r { } is the trace operator.
Theorem 1.
Denote the local state estimation achieved by the ith ( i = 1 , 2 , , N ) local filter by x ^ i ( k ) , and its corresponding error covariance matrix by P i i ( k ) (i.e., P i ( k ) , i = 1 , 2 , , N in the local filters). The cross-covariance matrix between the state estimation errors of x ^ i ( k ) and x ^ j ( k ) is P i j ( k ) ( i , j = 1 , 2 , , N ; i j ). Define ( Σ ( k ) ) i j = P i j ( k ) R n × n ( i , j = 1 , 2 , , N ), E ¯ = [ I n , I n , , I n ] T R n N × n and α ¯ = [ α 1 , α 2 , , α N ] T R n N × n . Based on the principle of linear minimum variance, the globally optimal state estimation x ^ ( k ) is obtained as
x ^ ( k ) = i = 1 N α i x ^ i ( k )
where α i ( i = 1 , 2 , , N ) is given by
α ¯ = Σ 1 ( k ) E ¯ [ E ¯ T Σ 1 ( k ) E ¯ ] 1
Proof. 
As E { x ^ i ( k ) x ( k ) } = 0 , taking the expectation of (25) and considering the relation i = 1 N α i = I n , we can obtain
E ( x ^ ( k ) ) = i = 1 N α i E ( x ^ i ( k ) ) = x ( k )
It can be seen from (27), x ^ ( k ) obtained by (25) is the unbiased estimation of x ( k ) .
Now let us consider Condition 2 to derive the globally optimal state estimation x ^ ( k ) . Define the estimation error of the local state estimation x ^ i ( k ) ( i = 1 , 2 , , N ) as
x ˜ i ( k ) = x ( k ) x ^ i ( k )
Thus, the estimation error of x ^ ( k ) is described as
x ˜ ( k ) = x ( k ) x ^ ( k ) = i = 1 N α i x ˜ i ( k )
According to (29), the error covariance matrix of x ^ ( k ) is
P ( k ) = α ¯ T Σ ( k ) α ¯
Applying the Lagrange multiplier to solve the minimum of t r { P ( k ) } , the Lagrange function is defined as
𝕃 = t r { P ( k ) } + t r { Λ [ α ¯ T E ¯ I n ] }
where Λ R n × n is the Lagrange multiplier.
Taking the partial derivative of (31) with respect to α ¯ and letting it to be zero, we have
Σ ( k ) α ¯ + 1 2 E ¯ Λ = 0
Combining i = 1 N α i = I n with (32), we readily have
[ Σ ( k ) E ¯ E ¯ T 0 ] [ α ¯ 1 2 Λ ] = [ 0 I n ]
By solving (33), it is verified that
[ α ¯ 1 2 Λ ] = [ Σ ( k ) E ¯ E ¯ T 0 ] 1 [ 0 I n ] = [ Σ 1 ( k ) E ¯ [ E ¯ T Σ 1 ( k ) E ¯ ] 1 [ E ¯ T Σ 1 ( k ) E ¯ ] 1 ]
Form (34), we have
α ¯ = Σ 1 ( k ) E ¯ [ E ¯ T Σ 1 ( k ) E ¯ ] 1
The proof of Theorem 1 is completed. ☐
It can be seen from (25) and (26), in order to calculate the global optimal state estimation x ^ ( k ) , it is necessary to predetermine the matrix Σ ( k ) . In the matrix Σ ( k ) , P i i ( k ) ( i = 1 , 2 , , N ) can be directly determined by the error covariance matrix of the state estimation in the ith local filter. However, it is difficult to calculate the cross-covariance matrix P i j ( k ) between the state estimation errors of x ^ i ( k ) and x ^ j ( k ) . This paper adopts the UT concept to approximate P i j ( k ) , and the derivation process is shown in Theorem 2.
Theorem 2.
For the multi-sensor system given by (1), the cross-covariance matrix P i j ( k ) between the estimation errors of local state estimations x ^ i ( k ) and x ^ j ( k ) ( i , j = 1 , 2 , , N ; i j ) is computed as
P i j ( k ) = s = 0 2 n ω s ( χ s , k / k 1 ( i ) x ^ i ( k / k 1 ) ) ( χ s , k / k 1 ( j ) x ^ j ( k / k 1 ) ) T [ s = 0 2 n ω s ( χ s , k / k 1 ( i ) x ^ i ( k / k 1 ) ) ( γ s , k / k 1 ( j ) z ^ j ( k / k 1 ) ) T ] K j T ( k ) K i ( k ) [ s = 0 2 n ω s ( γ s , k / k 1 ( i ) z ^ i ( k / k 1 ) ) ( χ s , k / k 1 ( j ) x ^ j ( k / k 1 ) ) T ] + K i ( k ) [ s = 0 2 n ω s ( γ s , k / k 1 ( i ) z ^ i ( k / k 1 ) ) ( γ s , k / k 1 ( j ) z ^ j ( k / k 1 ) ) T ] K j T ( k )
where χ s , k / k 1 ( i ) denotes the sigma point transformed by the nonlinear function f ( ) in (4) for the ith local filter; γ s , k / k 1 ( i ) denotes the sigma point transformed by the nonlinear function h i ( ) in (8) for the ith local filter; and s = 0 , 1 , , 2 n represents the order of the transformed sigma points.
Proof. 
Define the estimation error of the ith local filter as
x ( k ) x ^ i ( k ) = x ( k ) [ x ^ i ( k / k 1 ) + K i ( k ) ( z i ( k ) z ^ i ( k / k 1 ) ) ] = ( x ( k ) x ^ i ( k / k 1 ) ) K i ( k ) ( z i ( k ) z ^ i ( k / k 1 ) )
Then, it is verified from (37)
P i j ( k ) = E { [ x ( k ) x ^ i ( k ) ] [ x ( k ) x ^ j ( k ) ] T } = E { [ ( x ( k ) x ^ i ( k / k 1 ) ) K i ( k ) ( z i ( k ) z ^ i ( k / k 1 ) ) ]     [ ( x ( k ) x ^ j ( k / k 1 ) ) K j ( k ) ( z j ( k ) z ^ j ( k / k 1 ) ) ] T } = E { [ ( x ( k ) x ^ i ( k / k 1 ) ) ] [ ( x ( k ) x ^ j ( k / k 1 ) ) ] T } E { [ ( x ( k ) x ^ i ( k / k 1 ) ) ] [ ( z j ( k ) z ^ j ( k / k 1 ) ) ] T } K j T ( k ) K i ( k ) E { [ ( z i ( k ) z ^ i ( k / k 1 ) ) ] [ ( x ( k ) x ^ j ( k / k 1 ) ) ] T } + K i ( k ) E { [ ( z i ( k ) z ^ i ( k / k 1 ) ) ] [ ( z j ( k ) z ^ j ( k / k 1 ) ) ] T } K j T ( k )
Based on the tranformated sigma points χ s , k / k 1 ( i ) and γ s , k / k 1 ( i ) ( i = 1 , 2 , , N ), P i j ( k ) can be approximated by using UT
P i j ( k ) = s = 0 2 n ω s ( χ s , k / k 1 ( i ) x ^ i ( k / k 1 ) ) ( χ s , k / k 1 ( j ) x ^ j ( k / k 1 ) ) T [ s = 0 2 n ω s ( χ s , k / k 1 ( i ) x ^ i ( k / k 1 ) ) ( γ s , k / k 1 ( j ) z ^ j ( k / k 1 ) ) T ] K j T ( k ) K i ( k ) [ s = 0 2 n ω s ( γ s , k / k 1 ( i ) z ^ i ( k / k 1 ) ) ( χ s , k / k 1 ( j ) x ^ j ( k / k 1 ) ) T ] + K i ( k ) [ s = 0 2 n ω s ( γ s , k / k 1 ( i ) z ^ i ( k / k 1 ) ) ( γ s , k / k 1 ( j ) z ^ j ( k / k 1 ) ) T ] K j T ( k )
The proof of Theorem 2 is completed. ☐

3. Performance Evaluation and Discussion

A prototype system of INS/GNSS/CNS integration was implemented using the proposed AFUKF-MODF. Simulations and experiments as well as comparison analysis with FKF, UKF-FKF [13], and UKF-MODF [14] were conducted to comprehensively evaluate the performance of the proposed AFUKF-MODF.

3.1. System Model of INS/GNSS/CNS Integration

Fundamentally, INS/GNSS/CNS integration generates navigation solutions by utilizing the high-precision GNSS position and velocity to correct the INS velocity and position errors and utilizing the high-precision CNS attitude to correct the INS attitude error.

3.1.1. Process Model

The process model of the INS/GNSS/CNS integrated navigation system is established by combining the INS error equations with the inertial measurement unit (IMU) error equations.
The navigation frame (n-frame) is selected as the E-N-U (East-North-Up) geography frame (g-frame). Denote the inertial frame by i, the earth frame e, the body frame b and the INS simulated actual platform frame n . The system state vector is defined as
x ( t ) = [ ϕ E , ϕ N , ϕ U , δ v E , δ v N , δ v U , δ L , δ λ , δ h , ε x b , ε y b , ε z b , x b , y b , z b ] T
where ( ϕ E , ϕ N , ϕ U ) is the attitude error, ( δ v E , δ v N , δ v U ) the velocity error, ( δ L , δ λ , δ h ) the position error, ( ε x b , ε y b , ε z b ) the gyro constant drift, and ( x b , y b , z b ) the accelerometer zero-bias.
The nonlinear attitude and velocity error equations are given as [34,35]
{ ϕ ˙ = C ω 1 [ ( I C n n ) ω ^ i n n + C n n δ ω i n n C b n δ ω i b b ] δ v ˙ n = δ g n + [ I ( C n n ) T ] C b n f ^ b + C b n δ f b ( 2 ω ^ i e n + ω ^ e n n ) × δ v n ( 2 δ ω i e n + δ ω e n n ) × v n
where ϕ = ( ϕ E , ϕ N , ϕ U ) T and δ v n = ( δ v E , δ v N , δ v U ) T ; v n = ( v E , v N , v U ) T is the velocity of the vehicle in n-frame; C n n , C b n and C b n are the rotation matrices; δ g n is the gravity error; f ^ b is the measured specific force in the b-frame, and it is composed of accelerometer zero-bias b and its white noise ω a b ; δ f b is the corresponding error; δ ω i b b is the measurement error of the gyro, which is composed of gyro constant drift ε b and its white noise ω g b ; ω i e n is the rotational angular velocity of the earth; ω e n n is the angular velocity of the vehicle relative to the earth; ω i n n = ω i e n + ω e n n is the relative rotational angular velocity between the n-frame and i-frame; ω ^ i e n , ω ^ e n n and ω ^ i n n are the actual values of ω i e n , ω e n n and ω i n n in the n -frame; and δ ω i e n , δ ω e n n and δ ω i n n represent the corresponding errors. The above parameters can be calculated as
{ ω i e n = [ 0 ω i e cos L ω i e sin L ] T δ ω i e n = [ 0 δ L ω i e sin L δ L ω i e cos L ] T ω e n n = [ v N R M + h v E R N + h v E tan L R N + h ] T δ ω e n n = [ δ v N R M + h + δ h v N ( R M + h ) 2 δ v E R N + h δ h v E ( R N + h ) 2 δ v E tan L R N + h + δ L v x sec L 2 R N + h δ h v E sec L ( R N + h ) 2 ]
where L and h represent the latitude and altitude of the vehicle; and R M and R N are the median radius and normal radius.
C ω 1 is computed as
C ω 1 = 1 cos ϕ E [ cos ϕ N cos ϕ E 0 sin ϕ N cos ϕ E sin ϕ N sin ϕ E cos ϕ E cos ϕ N sin ϕ E sin ϕ N 0 cos ϕ N ]
The position error equation of INS is described by [36]
{ δ L ˙ = δ v N R M + h δ h v N ( R M + h ) 2 δ λ ˙ = δ v E sec L R N + h + δ L v E tan L sec L R N + h δ h v E sec L ( R N + h ) 2 δ h ˙ = δ v U
The gyro constant drift ε b and accelerometer zero-bias b are commonly described as random constants [14,30], i.e.,
ε ˙ i b = 0   ( i = x , y , z )
˙ i b = 0   ( i = x , y , z )
According to the selected system state vector, the process model of the INS/GNSS/CNS integration can be established by combining (41)–(46)
x ˙ ( t ) = f ¯ ( x ( t ) ) + w ( t )
where f ¯ ( ) is a nonlinear function describing the system state equation in continuous form; and w ( t ) = [ ( C ω 1 C b n ω a b ) T , ( C b n ω g b ) T , 0 1 × 9 ] T is the process noise vector.
By discretizing (47) with the improved Euler formulation [37], the discrete-time process model of the INS/GNSS/CNS integration is obtained as
x ( k ) = f ( x ( k 1 ) ) + w ( k )
where f ( ) is a nonlinear function describing the system state equation in discrete form; and w ( k ) is the discrete-time process noise vector.

3.1.2. Measurement Model of INS/GNSS Subsystem

Take the difference between INS and GNSS in terms of velocity and position as the measurement of the INS/GNSS subsystem, i.e.,
z 1 ( k ) = [ v E I v E G v N I v N G v U I v U G L I L G λ I λ G h I h G ] T
where ( v E I , v N I , v U I ) T and ( L I , λ I , h I ) T are the velocity and position output by INS; and ( v E G , v N G , v U G ) T and ( L G , λ G , h G ) T are the velocity and position achieved by GNSS.
Then, the measurement model of the INS/GNSS subsystem can be described as [11,24]
z 1 ( k ) = H 1 ( k ) x ( k ) + v 1 ( k ) = [ H v ( k ) H P ( k ) ] x ( k ) + [ v v ( k ) v P ( k ) ]
where H v ( k ) = [ 0 3 × 3 , I 3 × 3 , 0 3 × 9 ] , H P ( k ) = [ 0 3 × 6 , d i a g ( R M , R N cos L , 1 ) , 0 3 × 6 ] ; and v v ( k ) and v P ( k ) are the measurement noises, which correspond to the velocity and position errors of GNSS, respectively.

3.1.3. Measurement Model of INS/CNS Subsystem

CNS provides high-precise attitude information for a vehicle through the star sensor. Take the difference between INS and CNS in terms of attitude as the measurement of the INS/CNS subsystem, i.e.,
z 2 ( k ) = [ ϕ E I ϕ E C ϕ N I ϕ N C ϕ U I ϕ U C ]
where ( ϕ E I , ϕ N I , ϕ U I ) T and ( ϕ E C , ϕ N C , ϕ U C ) T are the attitude information obtained by INS and CNS, respectively.
Accordingly, the measurement model of the INS/CNS subsystem is constructed as [11,24]
z 2 ( k ) = H 2 ( k ) x ( k ) + v 2 ( k )
where H 2 ( k ) = [ I 3 × 3 , 0 3 × 12 ] , and v 2 ( k ) is the measurement noise corresponding to the measurement error of the star sensor.
The INS/GNSS/CNS integrated navigation system is a multi-sensor nonlinear system with 2 local filters (N = 2). The nonlinear system model, which is generally described by (1), is now specified by (48), (50) and (52) for the specific INS/GNSS/CNS integrated navigation system. The system models of the INS/GNSS and INS/CNS subsystems are described by (48) and (50), and (48) and (52), respectively. Based on these, the AFUKF based on Mahalanobis distance described in Section 2.1.2 serves as the local filters to calculate the local state estimations of the INS/GNSS and INS/CNS subsystems. Subsequently, both local state estimations of the INS/GNSS and INS/CNS subsystems are fused by the UT-based multi-sensor optimal data fusion described in Section 2.2 to calculate the globally optimal state estimation for the INS/GNSS/CNS integration.

3.2. Simulations and Analysis

Monte Carlo simulations were conducted to comprehensively evaluate the performance of the proposed AFUKF-MODF for the dynamic flight of an unmanned aerial vehicle (UAV) using the INS/GNSS/CNS integration for navigation and positioning. The flight trajectory, which involves various maneuvers such as climbing, pitching, rolling and turning, is shown in Figure 2. The simulation parameters are listed in Table 1. Both filtering period and fusion period were 1 s. The Monte Carlo simulations were carried out 100 times.
The initial state error covariances and process noise covariance for the local filters are set as
P i ( 0 ) = d i a g [ ( 1 ) 2 , ( 1 ) 2 , ( 1.5 ) 2 , ( 0.4   m / s ) 2 I 3 × 3 , ( 10   m ) 2 , ( 10   m ) 2 , ( 15   m ) 2 , ( 0.1 / h ) 2 I 3 × 3 , ( 1 × 10 3   g ) 2 I 3 × 3 ]   ( i = 1 , 2 )
Q = d i a g [ ( 0.05 / h ) 2 I 3 × 3 , ( 1 × 10 4   g s ) 2 I 3 × 3 , 0 9 × 9 ]
However, P i ( 0 ) ( i = 1 , 2 ) and Q were enlarged to the twice of their initial values for both FKF and UKF-FKF to eliminate the correlation between the two local state estimations by the use of covariance upper bound. Moreover, the process model of the INS/GNSS/CNS integration should be linearized in local fitering processes for the FKF.
The measurement noise covariances are set as
R 1 = d i a g [ ( 0.05   m / s ) 2 I 3 × 3 , ( 5   m ) 2 , ( 5   m ) 2 , ( 8   m ) 2 ]
R 2 = ( 5 ) 2 I 3 × 3
In order to evaluate the performance of the proposed AFUKF-MODF in terms of process-modeling error, the following process-modeling error is introduced to the process model of the INS/GNSS/CNS integration during the time interval from 400 s to 600 s
Δ x = [ 0 1 × 3 , 0.02   m / s , 0.02   m / s , 0.02   m / s , ( 2 × 10 6 ) , ( 2 × 10 6 ) , 5   m , 0 1 × 6 ] T
Accordingly, the process model used in the local filters is
{ x ( k ) = f ( x ( k 1 ) ) + w ( k ) ,   other   time   intervals x ( k ) = f ( x ( k 1 ) ) + Δ x + w ( k ) , ( 400   s , 600   s )
To identify the above process-modeling error, χ m 1 , α 2 and χ m 2 , α 2 used in the proposed AFUKF-MODF for the INS/GNSS subsystem and INS/CNS subsystem were chosen as 12.592 and 7.815, respectively. These values are resulted from the χ 2 distribution when the reliability level is 95% ( α = 0.05 ) and the degrees of freedom are 6 and 3, respectively.
For comparison analysis, simulation trials were conducted at the same conditions by using FKF, UKF-FKF, UKF-MODF and AFUKF-MODF, respectively. Moreover, the overall estimation error is adopted to evaluate the navigation accuracy in the simulation analysis for the four fusion methods. It is defined as the norm of the navigation parameters estimation error [38]
Δ x = Δ x E 2 + Δ x N 2 + Δ x U 2
where Δ x E , Δ x N and Δ x U are the components of Δ x in East, North and Up, respectively.
Figure 3, Figure 4 and Figure 5 depict the root mean squared errors (RMSEs) of overall attitude errors, velocity errors and position errors by the above four fusion methods, respectively. During the time intervals (0 s, 400 s) and (600 s, 1000 s) without process-modeling error involved, FKF has poor navigation accuracy compared to UKF-FKF, UKF-MODF and AFUKF-MODF. Its estimation RMSEs in attitude, velocity and position are around 0.2010 , 0.2084 m/s and 15.1129 m. This is because the use of covariance upper bound and the linearization of the system model cause suboptimal fusion performance. UKF-FKF adopts UKF to calculate the local state estimations, thus effectively refraining from the error caused by the linearization of system model and leading to the improved fusion accuracy. Its estimation RMSEs in attitude, velocity and position are around 0.1844 , 0.1647 m/s and 13.1980 m. However, due to the use of covariance upper bound, the fusion results of UKF-FKF are still suboptimal. Different to FKF and UKF-FKF, both UKF-MODF and AFUKF-MODF are directly derived based on the principle of linear minimum variance, avoiding the use of covariance upper bound in the local filters. Thus, the fusion results obtained by both UKF-MODF and AFUKF-MODF are globally optimal, which are much more accurate than those by FKF and UKF-FKF. The estimation RMSEs in attitude, velocity and position are around 0.1612 , 0.1155 m/s and 10.6099 m for UKF-MODF, and 0.1643 , 0.1176 m/s and 10.7212 m for AFUKF-MODF. Moreover, since there is no process-modeling error identified, the proposed AFUKF-MODF does not incorporate the adaptive fading factor into the local estimation processes. Accordingly, its navigation accuracy is close to that of UKF-MODF during these time intervals.
During the time interval (400 s, 600 s), due to the influence of the introduced process-modeling error, the navigation accuracy of FKF, UKF-FKF and UKF-MODF degrades seriously. This is because these three methods do not have the ability to resist the influence of process-modeling error on the fusion solution. The estimation RMSEs in attitude, velocity and position are around 0.2319 , 0.2471 m/s and 17.8467 m for FKF; 0.2078 , 0.1917 m/s and 15.2703 m for UKF-FKF; and 0.1887 , 0.1472 m/s and 13.1225 m for UKF-MODF. In contrast, the proposed AFUKF-MODF has higher navigation accuracy than FKF, UKF-FKF and UKF-MODF during this time interval. Its estimation RMSEs in attitude, velocity and position are around 0.1702 , 0.1280 m/s and 11.5798 m. This is because the proposed AFUKF-MODF can identify the process-modeling error according to the hypothesis testing theory and effectively refrain from the influence of process-modeling error on the fusion solution by adjusting the time-varying adaptive fading factor incorporated in the predicted state covariance matrix, leading to a strong adaptability and robustness.
The mean RMSEs of the over attitude errors, velocity errors and position errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the time interval (400 s, 600 s) and the other time intervals are listed in Table 2. The results in Table 2 also verify that the proposed AFUKF-MODF has a better adaptability and robustness than the other methods, thus leading to improved navigation accuracy for the INS/GNSS/CNS integration.
The above simulations and analysis demonstrate that the proposed AFUKF-MODF can overcome the limitation of using covariance upper bound and provide globally optimal fusion results. It can also enhance the adaptability and robustness against process-modeling error to overcome the limitation of UKF-MODF, thus leading to higher navigation accuracy than FKF, UKF-FKF and UKF-MODF in the presence of process-modeling error.

3.3. Experiments and Analysis

Practical experiments were also conducted to evaluate the performance of the proposed AFUKF-MODF by observing the flight of an UAV. As shown in Figure 6, the UAV uses an INS/GNSS/CNS integration system for navigation. This navigation system includes a MTi-100 IMU, a Hemisphere P307 BDS/GNSS receiver and a SODERN SED26 star sensor. The main parameters of the above three devices are listed in Table 3, Table 4 and Table 5. Moreover, another Hemisphere P370 BDS/GNSS receiver, which was placed at a local reference station (around 1km from the initial UAV position), was used along with the one mounted on the UAV to provide the differential GPS (DGPS) data. The maximal distance between the UAV and local reference station was less than 60 km to achieve the position accuracy of less than 0.1 m from the DGPS via post difference processing. The DGPS data were used as the reference values to evaluate the positioning error of the INS/GNSS/CNS integrated system.
The UAV flight test was conducted at the City of Yanliang in Shaanxi, China. The experimental navigation data were selected from the UAV flight test within a continuous time period of 1000 s, where different maneuvers were involved. Within the selected UAV flight test, the UAV initial position was at East longitude 109.217 ° , North latitude 34.647 ° and altitude 3525 m. The UAV initial velocity was 180 m/s, 150 m/s and 50 m/s in East, North and Up, respectively. The other initial parameters were set identically to those in the simulation case. The filtering periods of the local filters, i.e., the INS/GNSS and INS/CNS subsystems, were 1 s. The period for fusion of local state estimations was also 1 s.
Due to the complexity of the dynamic flight environment, the process model of the INS/GNSS/CNS integration involves modeling error during the UAV flight process. Figure 7 illustrates the UAV position errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF, respectively. As shown in Figure 7, the position error achieved by FKF is relatively large due to the influence of process-modeling error, the use of the covariance upper bound and the linearization of the process model. During the time period (100 s, 1000 s), the position errors in longitude, latitude and altitude for FKF are within (−23.0174 m, 20.7557 m), (−21.7712 m, 20.6464 m) and (−26.7975 m, 27.7306 m), respectively. UKF-FKF improves the fusion accuracy of FKF by adopting UKF to calculate the local state estimations, leading to the longitude, latitude and altitude errors within (−15.1817 m, 16.5516 m), (−16.5740 m, 15.9417 m) and (−21.8423 m, 24.0462 m). UKF-MODF can provide globally optimal fusion results based on the principle of linear minimum variance. Thus, it has higher navigation accuracy than FKF and UKF-FKF, leading to the longitude, latitude and altitude errors within (−12.0942 m, 11.8042 m), (−12.0311 m, 12.2480 m) and (−16.9458 m, 15.9384 m). However, its position error curve still has pronounced oscillations due to the influence of the process-modeling error. In contrast, the position errors in longitude, latitude and altitude for the proposed AFUKF-MODF are within (−7.1247 m, 7.8715 m), (−7.0576 m, 7.5495 m) and (−11.8677 m, 11.2610 m), which are much smaller than those by FKF, UKF-FKF and UKF-MODF. This is because the proposed AFUKF-MODF has the capability to inhibit the disturbances of process-modeling error.
Table 6 lists the mean absolute errors (MAEs) and standard deviations (STDs) of the position errors achieved by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF. It can be seen that the MAE and STD of the position errors by the proposed AFUKF-MODF are also much smaller than those of the other three methods.
The above experimental results and analysis demonstrate that the proposed AFUKF-MODF enhances the adaptability and robustness of multi-sensor data fusion in presence of process-modeling errors, leading to improved navigation accuracy for INS/GNSS/CNS integration in comparison with FKF, UKF-FKF and UKF-MODF.

4. Conclusions

This paper presents a novel AFUKF-MODF to improve the fusion adaptability and robustness for multi-sensor nonlinear stochastic systems against process-modeling error. The contributions of this paper are: (i) an AFUKF based on the Mahalanobis distance is developed to improve the adaptability and robustness of local state estimations against process-modeling error; and (ii) an UT-based multi-sensor optimal data fusion for the case of N local filters is also established to achieve the globally optimal state estimation via fusion of local state estimations. Simulations and practical experiments as well as comparison analysis demonstrate that the proposed AFUKF-MODF not only resists the disturbance of process-modeling error on the fusion solution, leading to improved adaptability and robustness for multi-sensor data fusion, but it also provides globally optimal data fusion results in the sense of linear minimum variance. The achieved fusion accuracy is much higher than that of FKF, UKF-FKF and UKF-MODF for INS/GNSS/CNS integration.
Future research work will focus on two aspects. One is the improvement of the proposed AFUKF-MODF. It is expected to combine the proposed AFUKF-MODF with artificial intelligence technologies such as machine learning, pattern recognition and neural networking to intelligently identify and compensate for system modeling error. The other is on the applications of the proposed AFUKF-MODF to address the multi-sensor data fusion problems in other fields such as target tracking, fault diagnosis and signal processing.

Acknowledgments

The work of this paper was supported by the National Natural Science Foundation of China (Project Numbers: 61174193&41704016) and the Specialized Research Fund for the Doctoral Program of Higher Education (Project Number: 20136102110036).

Author Contributions

Bingbing Gao and Gaoge Hu conceived and designed this paper; Gaoge Hu and Chengfan Gu performed the experiments and analyzed the data; Bingbing Gao wrote the paper; Shesheng Gao and Yongmin Zhong reviewed and edited the manuscript. All authors read and approved this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

FKFFederated Kalman filter
UKFUnscented Kalman filter
UTUnscented transformation
UKF-FKFUnscented Kalman filter-based federated Kalman filter
UKF-MODFUnscented Kalman filter-based multi-sensor optimal data fusion
AFUKFAdaptive fading unscented Kalman filter
AFUKF-MODFAdaptive fading unscented Kalman filter-based multi-sensor optimal data fusion
INSInertial navigation system
GNSSGlobal navigation Satellite System
CNSCelestial navigation system
IMUInertial measurement unit
E-N-UEast-north-up
UAVUnmanned aerial vehicle
RMS/RMSERoot mean squared error
BDSBeidou system
GPSGlobal positioning system
DGPSDifferential global positioning system
FOVField of view
MAEMean absolute error
STDStandard deviation

References

  1. Gravina, R.; Alinia, P.; Ghasemzadeh, H.; Fortinoa, G. Multi-sensor fusion in body sensor networks: State-of-the-art and research challenges. Inf. Fusion 2017, 35, 68–80. [Google Scholar] [CrossRef]
  2. Gao, S.S.; Zhong, Y.M.; Zhang, X.Y.; Shirinzadeh, B. Multi-sensor optimal data fusion for INS/GPS/SAR integrated navigation system. Aerosp. Sci. Technol. 2009, 13, 232–237. [Google Scholar] [CrossRef]
  3. Khaleghi, B.; Khamis, A.; Karray, F.O.; Razavi, S.N. Multisensor data fusion: A review of the state-of-the-art. Inf. Fusion 2013, 14, 28–44. [Google Scholar] [CrossRef]
  4. Gao, J.B.; Harris, C. Some remarks on Kalman filters for the multisensor fusion. Inf. Fusion 2002, 3, 191–201. [Google Scholar] [CrossRef]
  5. Qi, W.J.; Zhang, P.; Nie, G.H.; Deng, Z.L. Robust weighted fusion Kalman predictors with uncertain noise variances. Digit. Signal Process. 2014, 30, 37–54. [Google Scholar] [CrossRef]
  6. Caballero-Águila, R.; García-Garrido, I.; Linares-Pérez, J. Information fusion algorithms for state estimation in multi-sensor systems with correlated missing measurements. Appl. Math. Comput. 2014, 226, 548–563. [Google Scholar] [CrossRef]
  7. Gao, Y.; Krakiwsky, E.J.; Abousalem, M.A.; Mclellan, J.F. Comparison and analysis of centralized, decentralized, and federated Filters. Navigation 1993, 40, 69–86. [Google Scholar] [CrossRef]
  8. Ma, J.; Sun, S.L. Centralized fusion estimators for multisensor systems with random sensor delays, multiple packet dropouts and uncertain observations. IEEE Sens. J. 2013, 13, 1228–1235. [Google Scholar] [CrossRef]
  9. Mahmoud, M.S.; Khalid, H.M. Distributed Kalman filtering: A bibliographic review. IET Control Theory Appl. 2013, 7, 483–501. [Google Scholar] [CrossRef]
  10. Feng, J.X.; Wang, Z.D.; Zeng, M. Distributed weighted robust Kalman filter fusion for uncertain systems with autocorrelated and cross-correlated noises. Inf. Fusion 2013, 14, 78–86. [Google Scholar] [CrossRef]
  11. Hu, G.G.; Gao, S.S.; Zhong, Y.M.; Gao, B.B.; Subic, A. Modified federated Kalman filter for INS/GNSS/CNS integration. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2016, 230, 30–44. [Google Scholar] [CrossRef]
  12. Carlson, N.A.; Berarducci, M.P. Federated Kalman filter simulation results. Navigation 1994, 41, 297–322. [Google Scholar] [CrossRef]
  13. Hu, H.D.; Huang, X.L. SINS/CNS/GPS integrated navigation algorithm based on UKF. J. Syst. Eng. Electron. 2010, 21, 102–109. [Google Scholar] [CrossRef]
  14. Gao, B.B.; Hu, G.G.; Gao, S.S.; Zhong, Y.M.; Gu, C.F. Multi-sensor optimal data fusion for INS/GNSS/CNS integration based on unscented Kalman filter. Int. J. Control Autom. Syst. 2018, 16, 129–140. [Google Scholar] [CrossRef]
  15. Sun, S.L.; Deng, Z.L. Multi-sensor optimal information fusion Kalman filter. Automatica 2004, 40, 1017–1023. [Google Scholar] [CrossRef]
  16. Hu, G.G.; Gao, S.S.; Zhong, Y.M. A derivative UKF for tightly coupled INS/GPS integrated navigation. ISA Trans. 2015, 56, 135–144. [Google Scholar] [CrossRef] [PubMed]
  17. Hu, G.G.; Gao, S.S.; Zhong, Y.M.; Gao, B.B. Stochastic stability of the derivative unscented Kalman filter. Chin. Phys. B 2015, 24, 070202. [Google Scholar] [CrossRef]
  18. Gao, Z.H.; Mu, D.J.; Gao, S.S.; Zhong, Y.M.; Gu, C.F. Adaptive unscented Kalman filter based on maximum posterior and random weighting. Aerosp. Sci. Technol. 2017, 71, 12–24. [Google Scholar] [CrossRef]
  19. Gao, B.B.; Gao, S.S.; Zhong, Y.M.; Hu, G.G.; Gu, C.F. Interacting multiple model estimation-based adaptive robust unscented Kalman filter. Int. J. Control Autom. Syst. 2017, 15, 2013–2025. [Google Scholar] [CrossRef]
  20. Gao, B.B.; Gao, S.S.; Hu, G.G.; Zhong, Y.M.; Gu, C.F. Maximum likelihood principle and moving horizon estimation-based adaptive unscented Kalman filter. Aerosp. Sci. Technol. 2018, 73, 184–196. [Google Scholar] [CrossRef]
  21. Gao, S.S.; Hu, G.G.; Zhong, Y.M. Windowing and random weighting-based adaptive unscented Kalman filter. Int. J. Adapt. Control Signal Process. 2015, 29, 201–223. [Google Scholar] [CrossRef]
  22. Chang, K.C.; Zhi, T.; Saha, R.K. Performance evaluation of track fusion with information matrix filter. IEEE Trans. Aerosp. Electron. Syst. 2002, 38, 455–466. [Google Scholar] [CrossRef]
  23. Gao, S.S.; Zhong, Y.M.; Shirinzadeh, B. Random weighting estimation for fusion of multi-dimensional position data. Inf. Sci. 2010, 180, 4999–5007. [Google Scholar] [CrossRef]
  24. Hu, G.G.; Gao, S.S.; Zhong, Y.M.; Gao, B.B.; Subic, A. Matrix weighted multi-sensor data fusion for INS/GNSS/CNS integration. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2016, 230, 1011–1026. [Google Scholar] [CrossRef]
  25. Hu, G.G.; Gao, S.S.; Zhong, Y.M.; Gao, B.B.; Subic, A. Modified strong tracking unscented Kalman filter for nonlinear state estimation with process model uncertainty. Int. J. Adapt. Control Signal Process. 2015, 29, 1561–1577. [Google Scholar] [CrossRef]
  26. Chang, L.B.; Li, K.L.; Hu, B.Q. Huber’s M-estimation-based Process Uncertainty Robust Filter for Integrated INS/GPS. IEEE Sens. J. 2015, 15, 3367–3374. [Google Scholar] [CrossRef]
  27. Chang, G.B.; Liu, M. An adaptive fading Kalman filter based on Mahalanobis distance. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2015, 229, 1114–1123. [Google Scholar] [CrossRef]
  28. Soken, H.E.; Hajiyev, C. Adaptive fading UKF with Q-adaptation: Application to picosatellite attitude estimation. J. Aerosp. Eng. 2011, 26, 628–636. [Google Scholar] [CrossRef]
  29. Song, Q.; Han, J.D. An adaptive UKF algorithm for the state and parameter estimation of a mobile robot. Acta Autom. Sin. 2008, 34, 72–79. [Google Scholar] [CrossRef]
  30. Meng, Y.; Gao, S.S.; Zhong, Y.M.; Hu, G.G.; Subic, A. Covariance matching-based adaptive unscented Kalman filter for direct filtering in INS/GNSS integration. Acta Astronaut. 2016, 120, 171–181. [Google Scholar] [CrossRef]
  31. Chang, G.B. Robust Kalman filtering based on Mahalanobis distance as outlier judging criterion. J. Geodesy 2014, 88, 391–401. [Google Scholar] [CrossRef]
  32. Allen, M.B.; Isaacson, E.L. Numerical Analysis for Applied Science; John Wiley & Sons: New York, NY, USA, 1998. [Google Scholar]
  33. Kelley, C.T. Solving Nonlinear Equations with Newton's Method; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2003. [Google Scholar]
  34. Cui, B.B.; Chen, X.Y.; Tang, X.H. Improved Cubature Kalman Filter for GNSS/INS Based on Transformation of Posterior Sigma-Points Error. IEEE Trans. Signal Process. 2017, 65, 2975–2987. [Google Scholar] [CrossRef]
  35. Sun, F.; Tang, L.J. INS/GPS integrated navigation filter algorithm based on cubature Kalman filter. Control Decis. 2012, 27, 1032–1036. [Google Scholar]
  36. Qin, Y.Y. Inertial Navigation; Science Press: Beijing, China, 2006. [Google Scholar]
  37. Xiong, K.; Liu, L.D.; Zhang, H.Y. Modified unscented Kalman filtering and its application in autonomous satellite navigation. Aerosp. Sci. Technol. 2009, 13, 238–246. [Google Scholar] [CrossRef]
  38. Wang, D.J.; Lv, H.F.; Wu, J. Augmented cubature Kalman filter for nonlinear RTK/MIMU integrated navigation with non-additive noise. Measurement 2016, 97, 111–125. [Google Scholar] [CrossRef]
Figure 1. The framework of the proposed AFUKF-based multi-sensor optimal data fusion.
Figure 1. The framework of the proposed AFUKF-based multi-sensor optimal data fusion.
Sensors 18 00488 g001
Figure 2. Flight trajectory of the UAV.
Figure 2. Flight trajectory of the UAV.
Sensors 18 00488 g002
Figure 3. The RMSEs of overall attitude errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the simulation case.
Figure 3. The RMSEs of overall attitude errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the simulation case.
Sensors 18 00488 g003
Figure 4. The RMSEs of overall velocity errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the simulation case.
Figure 4. The RMSEs of overall velocity errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the simulation case.
Sensors 18 00488 g004
Figure 5. The RMSEs of overall position errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the simulation case.
Figure 5. The RMSEs of overall position errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the simulation case.
Sensors 18 00488 g005
Figure 6. Navigation setup of the UAV.
Figure 6. Navigation setup of the UAV.
Sensors 18 00488 g006
Figure 7. The position errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the UAV experiment test: (a) Longitude error; (b) Latitude error; (c) Altitude error.
Figure 7. The position errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the UAV experiment test: (a) Longitude error; (b) Latitude error; (c) Altitude error.
Sensors 18 00488 g007aSensors 18 00488 g007b
Table 1. Simulation parameters.
Table 1. Simulation parameters.
ParametersValues
Initial parametersInitial position (longitude-latitude-altitude)(108.997°, 34.246°, 5000 m)
Initial velocity (east-north-up)(0 m/s, 150 m/s, 0 m/s)
Initial attitude (pitch-roll-yaw)(0°, 0°, 0°)
Initial parameter errorsInitial position error (longitude-latitude-altitude)(10 m, 10 m, 15 m)
Initial velocity error (east-north-up)(0.4 m/s, 0.4 m/s, 0.4 m/s)
Initial attitude error (pitch-roll-yaw)(1′, 1′, 1.5′)
INS parametersGyro parametersConstant drift0.1 °/h
Random walk coefficient0.05 °/ h
Sampling frequency20 Hz
Accelerometer parametersZero-bias1 × 10−3 g
Random walk coefficient1 × 10−4 g × s
Sampling frequency20 Hz
GNSS parametersHorizontal position error (RMS)5 m
Altitude error (RMS)8 m
Velocity error (RMS)0.05 m/s
Data update rate1 Hz
CNS parametersAttitude error (RMS)5″
Data update rate1 Hz
Simulation time1000 s
Table 2. Mean RMSEs of the overall estimation errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the simulation case.
Table 2. Mean RMSEs of the overall estimation errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the simulation case.
Data fusion MethodsNavigation ErrorsMean RMSE
(400 s, 600 s)The Other Time Intervals
FKFAttitude error (′)0.23190.2010
Velocity error (m/s)0.24710.2084
Position error (m)17.846715.1129
UKF-FKFAttitude error (′)0.20780.1844
Velocity error (m/s)0.19170.1647
Position error (m)15.270313.1980
UKF-MODFAttitude error (′)0.18870.1612
Velocity error (m/s)0.14720.1155
Position error (m)13.122510.6099
AFUKF-MODFAttitude error (′)0.17020.1643
Velocity error (m/s)0.12800.1176
Position error (m)11.579810.7212
Table 3. Noise parameters of MTi-100 IMU.
Table 3. Noise parameters of MTi-100 IMU.
Noise ParametersValues
GyroConstant drift 10 ° / h
Random walk coefficient 0.6 / h
AccelerometerZero-bias 40   μ g
Random walk coefficient 80   μ g × s
Table 4. Main parameters of Hemisphere P307 BDS/GNSS receiver.
Table 4. Main parameters of Hemisphere P307 BDS/GNSS receiver.
Feature ParametersValues
Satellite signalsBDS(B1, B2, B3), GPS(L1, L2), GLONASS(G1, G2)
Horizontal position error (RMS)1.2 m
Altitude error (RMS)3 m
Velocity error (RMS)0.02 m/s
Data update rate20 Hz
Table 5. Main parameters of SODERN SED26 star sensor.
Table 5. Main parameters of SODERN SED26 star sensor.
Feature ParametersValues
Field of view (FOV)15.437° × 15.437° (Wide FOV)
Observable star number 6
Attitude error (RMS) 5
Data update rate10 Hz
Table 6. MAEs and STDs of the position errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the UAV experiment test.
Table 6. MAEs and STDs of the position errors obtained by FKF, UKF-FKF, UKF-MODF and AFUKF-MODF for the UAV experiment test.
Data Fusion MethodsPosition
LongitudeLatitudeAltitude
FKFMAE (m)7.27197.37079.0803
STD (m)8.83258.915011.1273
UKF-FKFMAE (m)5.45305.59248.0972
STD (m)6.58796.84479.8784
UKF-MODFMAE (m)3.88723.89966.1051
STD (m)4.80514.81427.3485
AFUKF-MODFMAE (m)2.35032.36104.0656
STD (m)2.90152.92254.9706

Share and Cite

MDPI and ACS Style

Gao, B.; Hu, G.; Gao, S.; Zhong, Y.; Gu, C. Multi-Sensor Optimal Data Fusion Based on the Adaptive Fading Unscented Kalman Filter. Sensors 2018, 18, 488. https://doi.org/10.3390/s18020488

AMA Style

Gao B, Hu G, Gao S, Zhong Y, Gu C. Multi-Sensor Optimal Data Fusion Based on the Adaptive Fading Unscented Kalman Filter. Sensors. 2018; 18(2):488. https://doi.org/10.3390/s18020488

Chicago/Turabian Style

Gao, Bingbing, Gaoge Hu, Shesheng Gao, Yongmin Zhong, and Chengfan Gu. 2018. "Multi-Sensor Optimal Data Fusion Based on the Adaptive Fading Unscented Kalman Filter" Sensors 18, no. 2: 488. https://doi.org/10.3390/s18020488

APA Style

Gao, B., Hu, G., Gao, S., Zhong, Y., & Gu, C. (2018). Multi-Sensor Optimal Data Fusion Based on the Adaptive Fading Unscented Kalman Filter. Sensors, 18(2), 488. https://doi.org/10.3390/s18020488

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