Next Article in Journal
Investigation on the Model-Based Control Performance in Vehicle Safety Critical Scenarios with Varying Tyre Limits
Previous Article in Journal
Design and Research of All-Terrain Wheel-Legged Robot
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

UAV Swarm Navigation Using Dynamic Adaptive Kalman Filter and Network Navigation

School of Instrument Science and Optoelectronics Engineering, Beihang University, XueYuan Road No. 37, HaiDian District, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(16), 5374; https://doi.org/10.3390/s21165374
Submission received: 13 July 2021 / Revised: 4 August 2021 / Accepted: 6 August 2021 / Published: 9 August 2021
(This article belongs to the Topic GNSS Measurement Technique in Aerial Navigation)

Abstract

:
Aiming to improve the positioning accuracy of an unmanned aerial vehicle (UAV) swarm under different scenarios, a two-case navigation scheme is proposed and simulated. First, when the Global Navigation Satellite System (GNSS) is available, the inertial navigation system (INS)/GNSS-integrated system based on the Kalman Filter (KF) plays a key role for each UAV in accurate navigation. Considering that Kalman filter’s process noise covariance matrix Q and observation noise covariance matrix R affect the navigation accuracy, this paper proposes a dynamic adaptive Kalman filter (DAKF) which introduces ensemble empirical mode decomposition (EEMD) to determine R and adjust Q adaptively, avoiding the degradation and divergence caused by an unknown or inaccurate noise model. Second, a network navigation algorithm (NNA) is employed when GNSS outages happen and the INS/GNSS-integrated system is not available. Distance information among all UAVs in the swarm is adopted to compensate the INS position errors. Finally, simulations are conducted to validate the effectiveness of the proposed method, results showing that DAKF improves the positioning accuracy of a single UAV by 30–50%, and NNA increases the positioning accuracy of a swarm by 93%.

1. Introduction

In recent years, the unmanned aerial vehicle (UAV) swarm has been given more and more attention due to the tremendous application prospects. Compared with a single UAV, the UAV swarm is more efficient, reliable, flexible and intensive [1,2]. As a prerequisite for a UAV to complete various tasks, the navigation and positioning technology has been extensively studied by scholars. At present, the most popular navigation system usually consists of an inertial navigation system (INS) and a Global Navigation Satellite System (GNSS) [3]. This system is very effective for improving navigation accuracy; however, its disadvantage of a fragile robustness due to the easy loss or disturbance of the GNSS signal limits its applications. When the GNSS signal is available, the Kalman filter (KF) is the most popular INS/GNSS data fusion algorithm. This paper aims to improve the fusion accuracy of the traditional KF and solve the problem of the UAV swarm navigation in the case of GNSS outage.
The Kalman filter uses the next observation to update the estimates of the current state variables, along with their uncertainties. KF corrects the noisy sensor data to a more realistic sensor output (with less noise) [4]. However, the performance of KF is affected by the stochastic model describing the process noise and observation noise and the dynamic model describing the dynamics of the system over time [5]. A variety of adaptive Kalman filters (AKF) has been proposed by scholars to improve the fusion precision of traditional Kalman filters when the stochastic model is inaccurate [6]. The AKF based on the attenuation factor reduces the influence of historical data on the current results by adjusting the weight of the one-step prediction covariance matrix. In fact, the problem of the inaccurate noise model has not been solved [7]. In most adaptive Kalman filtering algorithms, one of Q and R is usually fixed to adjust the other adaptively. For example, Li et al. proposed an adaptive Kalman filter based on the ensemble empirical mode (EMD) for X-ray pulsar navigation on the premise that the model has an accurate process noise covariance matrix Q [8]. The adaptive Kalman filter algorithm based on an innovation sequence and pseudo-measurement vector method proposed by Zhang [9] and the improving adaptive Kalman estimation algorithm based on a covariance-matching technology mentioned in [10] are implemented on the assumption that R is completely known. The multi-model adaptive Kalman filter can estimate the unknown probability of measurement loss, but its premise is that both Q and R are clear and the calculation of the algorithm requires multiple Kalman filters [11]. However, these three premises are difficult to achieve in practical applications.
Unfortunately, in case of GNSS outage, the INS/GNSS-integrated system does not work. The UAV position estimation based solely on the INS will drift significantly over time [12]. At present, the artificial neural network (ANN) is the mainstream solution to deal with this situation. The introduction of the ANN was essentially to establish a nonlinear relationship between different inputs and outputs, and then estimate and correct inertial navigation errors when the GNSS is unavailable [13]. For example, the author uses the position information output by the INS and the output of the interactive multi-model extended Kalman filter (IMM-EKF) to train the extreme learning machine when the GNSS is available. This training model can compensate for divergent position errors during GNSS outage [14]. Wei et al. trained a wavelet neural network based on a random forest regression to provide the observation information for the AKF during a GNSS outage [15]. Chen et al. introduced a wavelet neural network (WNN) to establish the time-varying characteristic model of INS error [16]. In this system, historical error data were used to predict the current error when the GNSS signal was blocked. Abdolkarimi et al. [17] introduced a neural network based on an extreme learning machine to compensate for INS errors during a GNSS outage, and improve the signal-to-noise ratio of the sensor measurements through wavelet transform. An ensemble learning method based on the multilayer perceptron (MLP) is proposed to provide a continuous position estimation for ships when GNSS signal is unavailable [18]. However, as we all know, the prediction result of the neural network depends heavily on its own structure and training data. In addition, a lot of training time is also necessary to maintain a satisfactory accuracy.
The purpose of this paper is to provide a navigation solution that meets multiple situations (with or without GNSS). The main contents are as follows:
First, aiming at the situation that R and Q will affect each other when the adaptive Kalman filter works [10], the dynamic adaptive Kalman filter (DAKF) based on ensemble empirical mode decomposition (EEMD) is designed. This method applies EEMD to the noise extraction of multi-channel GNSS signals and estimates the noise variance to dynamically adjust the observation noise covariance matrix R. Then, through the comparison of the actual and theoretical innovation sequences, Q is adaptively modified [10]. Through the above work, the influence of the uncertain stochastic model on the performance of KF will be eliminated. A simulation will show that the DAKF has an effective improvement compared with the traditional KF.
Second, a network navigation algorithm (NNA) based on the distances among all UAVs is used. If a GNSS outage occur, the UAV swarm will automatically adjust from integrated navigation to network navigation. The NNA applies the translation–rotation approach between the measurement network with measured distances between nodes as edges and the inertial position network with inertial distances calculated from an INS output as edges to estimate the inertial positioning errors of each node in the swarm [19]. The simulation will verify that NNA can significantly improve the UAV swarm positioning accuracy.
The main contributions of this paper are as follows:
Compared with the adaptive Kalman filter algorithm mentioned in [7,8,9,10,11], the algorithm proposed in this paper can simultaneously estimate Q and R. This estimate changes dynamically with time and has a strong robustness to changing noise. The proposed algorithm also guarantees the positive semi-definiteness of the estimated state covariance matrix by using its Joseph’s form. In addition, compared with the EMD-based adaptive filtering method mentioned in [8], the proposed algorithm has the following advantages: (1) EMD suffers from the problem of mode mixing, which causes insufficient decomposition of the original signal. The EEMD algorithm introduced in this paper is not bothered by this problem [20]. This paper may be the first approach to introduce EEMD into the parameter solving of the Kalman filter in the field of integrated navigation. (2) All historical data were used for noise extraction in [8], which will lead to that R cannot effectively follow the changes of environmental noise. The proposed algorithm adds a sliding window for noise extraction so that R can better respond to the changes of real noise. (3) The noise is calculated by the fixed number of the intrinsic mode function (IMF) components in [8], which is unrealistic. The proposed algorithm uses the power spectral density (PSD) feature of the IMF to select different components to reconstruct the noise signal in different environments, which is adaptive [21]. Compared with the multi-model adaptive Kalman filter in [11], the algorithm proposed in this paper only needs one Kalman filter, which will save resources and time for calculation.
Compared with the neural network-based algorithm mentioned in [14,15,16,17,18], this paper introduces the NNA to compensate the divergent position error of the UAV swarm during a GNSS outage. The proposed algorithm does not require a large number of network parameters and training resources, and the compensation result is not affected by historical data. In addition, this paper designs a simulation to prove the superiority of the proposed algorithm compared with the neural network-based algorithm.
The rest of the paper is organized as follows. Section 2 introduces the EEMD and the adaptive Kalman filter. Then, a dynamic adaptive Kalman filter algorithm is derived. Section 3 describes the principles of the networked navigation algorithm. Section 4 provides the test results and analysis to show the superiority of the proposed methods. The conclusion is given in Section 5.

2. Dynamic Adaptive Kalman Filter

When the observation noise covariance matrix matches the real noise, the adjustment of the process noise covariance matrix can help the KF perform better. However, in practical applications, it is extremely difficult to obtain the real noise which will vary with the measurement environment and sensor characteristics in the measurement process [22]. This situation not only damages the performance of the Kalman filter itself, but also leads to an erroneous adjustment of Q, which makes the performance of the AKF inferior to the ordinary KF.
In order to solve this problem, a new adaptive algorithm that can simultaneously adjust Q and R in real time is introduced. When the GNSS signal changes under the environment influence or intentional jamming, the observation noise possesses the characteristics of non-linearity and non-uniformity. General signal processing methods, such as the Fourier transform, have difficulty dealing with this problem. Therefore, this paper introduces EEMD which is self-adaptive into the AKF to estimate the observation noise.

2.1. Ensemble Empirical Mode Decomposition

Huang invented a signal processing algorithm called EMD that can decompose a nonlinear signal into a series of near-orthogonal intrinsic mode functions (IMFs) [23]. EMD is an adaptive time–frequency data processing algorithm, which is widely used to extract signals of interest from non-linear and non-stationary processes containing noise [24,25]. However, it is well known that EMD suffers from a mode mixing problem, which is defined as signals of different frequencies contained in a single IMF, that is, the decomposition of the original signal is not sufficient [26]. To solve this problem, EEMD based on noise assistance is proposed [20].
The principle of EEMD is very complete.
Adding white noise of limited amplitude to the original data will help the work of EMD. Signals of different scales can be automatically projected onto appropriate IMF components [20]. This phenomenon depends on the statistical characteristics of the white noise itself [27].
Since the added white noise in each experiment is random, after averaging the results of enough experiments, the influence of the added white noise on the original signal will be eliminated [20].
By combining the original signal with white noise of limited amplitude, we can obtain a clearer IMF component [26]. The pseudo-code implementation of the EEMD is listed in Algorithm 1.
Algorithm 1 EEMD
Input:
o ( t ) : The original signal;
begin
 (1) o ( t ) = o ( t ) + n ( t ) //Add Gaussian noise series to the original signal o ( t ) ; n ( t ) is the Gaussian noise series with amplitude of σ n ;
 (2) I M F s ( t ) = E M D ( o ( t ) ) //Decompose noise-added signal o ( t ) into I M F s ( t ) using EMD;
 (3) Repeat (1) and (2) M times with different Gaussian noise series which have the same amplitude;
 (4) Get the means of corresponding I M F s ( t ) of the decompositions as the final I M F s ( t ) ;
end
As a result, the original signal is decomposed into the following form:
o ( t ) = f = 1 F I M F f ( t ) + r ( t )
where o ( t ) is the original signal, I M F f ( t ) is the fth IMF signal and r ( t ) is the residual signal.

2.2. Acquisition of Observation Noise Covariance Matrix R

The general steps of the acquisition of R are as follows:
The EEMD algorithm decomposes the GNSS signal within a fixed window size w into a series of IMF components (the introduction of the window makes the calculation of R pay more attention to the current noise situation instead of obtaining R from the data of the entire historical period, which wastes a lot of resources and is inaccurate). Then, the noise components are distinguished by comparing the power spectral density (PSD) of each component [21] and the original noise is reconstructed by the summation of the noise components. At last, the observed noise covariance matrix R is estimated by the standard deviation of the reconstructed noise.
Since EEMD requires a suitable window size w for the noise extraction, the setting of R is determined by the empirical value in the first w seconds, but the inaccuracy of R within a short period of time will not have a significant impact on the final result.
The calculation process of the observation noise covariance matrix R using EEMD is listed in Algorithm 2.
Algorithm 2 f E E M D ( ) : observation noise covariance matrix by EEMD
Input:
Z G N S S : a numeric sequence of GNSS signal;
w : a window size used for computing R;
Output:
R k : observation noise covariance matrix at time k;
begin
 (1) I M F s : = E E M D ( Z G N S S [ k w : k ] ) //Decompose input signal using EEMD;
 (2) I M F n i o s e : = S e l e c t ( I M F s ) // Use power spectral density (PSD) to select IMF [21];
 (3) Z n i o s e : = S u m ( I M F n i o s e ) // Reconstruct the noise signal by summing I M F n i o s e ;
 (4) R k : = C a l c u l a t e V A R ( Z n i o s e ) // Calculate the variance of Z n i o s e ;
 (5) Return R k // return the estimated value of R at time k;
end
The inputs of this function are the GNSS signal Z G N S S and the window size w for filtering, and the output is R at time k which changes with the actual noise.

2.3. INS Error Model

In this section, we will describe the INS error model in the east–north–up navigation reference frame. The error equation is as follows:
(1)
Attitude error equation:
φ ˙ n = δ ω i e n + δ ω e n n + ( ω i e n + ω e n n ) × φ n δ ω i b n
(2)
Velocity error equation:
δ v ˙ n = φ n × f n ( 2 δ ω i e n + δ ω e n n ) × v n ( 2 ω i e n + ω e n n ) × δ v n + δ f n
(3)
Position error equation:
δ L ˙ = δ v N R M + h δ λ ˙ = δ v E R N + h sec L + v E R N + h sec L tan L δ L δ h ˙ = δ v U
Navigation information errors include: the attitude error shown in the navigation frame φ n = ( φ E   φ N   φ U ) the velocity error δ v n = ( δ v E   δ v N   δ v U ) and the position error δ P = ( δ L   δ λ   δ h ) . f n is the output of the accelerometer. R M and R N are the radii of the earth in the latitude and longitude directions, respectively. ω i e n and ω e n n are the rotational angular velocities of the earth and navigation frame, respectively. δ ω i b n and δ f n are the errors of the gyro and accelerometer, respectively. δ ω i b n consists of the gyro bias ε n and white noise ω g n . δ f n is composed of the accelerometer bias n and white noise ω a n .

2.4. The Filtering Algorithm of Dynamic Adaptive Kalman Filter

Consider the following multivariable linear discrete system [28]:
X k = Φ k / k 1 X k 1 + W k 1 Z k = H k X k + V k
where X k is the state vector ( X k = [ δ P   δ v n   φ n   ε n   n ] ), Φ k / k 1 is the transition matrix formed by the error equation described in Section 2.3, Z k is the observation vector formed by the position difference of the GNSS receiver and the INS and H k is the observation matrix. W k 1 and V k are the Gaussian white noise sequences satisfying the following conditions:
{ E [ W k ] = 0 , E [ W k W j T ] = Q k j = k E [ V k ] = 0 , E [ V k V j T ] = R k j = k E [ W k V j T ] = 0
where Q k is the process noise covariance matrix and R k is the observation noise covariance matrix.
The prediction equations are:
X ^ k / k 1 = Φ k / k 1 X ^ k 1 P ^ k / k 1 = Φ k / k 1 P ^ k 1 Φ k / k 1 T + Q k
where X ^ k / k 1 is the predicted state vector; X ^ k 1 is the estimated state vector; P ^ k / k 1 is the predicted state covariance matrix; P ^ k 1 is the estimated state covariance matrix.
The update equations are:
K k = P ^ k / k 1 H k T ( H k P ^ k / k 1 H k T + R k ) 1 X ^ k = X ^ k / k 1 + K k ( Z k H k X ^ k / k 1 ) P ^ k = ( I K k H k ) P ^ k / k 1 ( I K k H k ) T + K k R k K k T
where K k is the Kalman gain and H k is the observation matrix. It is worth mentioning that in order to ensure the positive semi-definiteness of P ^ k , we used its Joseph’s form [28].
In order to settle the uncertainty of the process noise, Ding et al. proposed an improving AKF based on the covariance matching method under the assumption that the observation noise is completely known and rewritten (7) as follows [10].
X ^ k / k 1 = Φ k / k 1 X ^ k 1 P ^ k / k 1 = Φ k / k 1 P ^ k 1 Φ k / k 1 T + s k Q k 1
where s k is the scaling factor expressed as follows:
s k = t r a c e { 1 m e = 0 m 1 d k e d k e T R k } t r a c e { H k ( Φ k / k 1 P k 1 Φ k / k 1 T + Q k 1 ) H k T }
The innovation sequence d k is defined as:
d k = Z k H k X ^ k / k 1
When R k is unknown or inaccurate, the calculation of the scaling factor will be affected. The adjustment of Q will become worse. When R k is accurately calculated by the proposed method, the solution of the scaling factor will be accurate. Therefore, Q will be adjusted more appropriately.
After substituting the R k calculated by Algorithm 2 into (10) used to solve the scaling factor, the complete process of DAKF is shown in Figure 1.

3. Network Navigation Algorithm

When the UAV swarm cannot receive the GNSS signal, the three-dimensional network navigation technology based on the distance information among all UAVs is introduced to reduce the error of the INS [19].
The specific process of the algorithm is as follows:
By linearizing the difference between the inertial distances calculated from the INS output and the measured distance between nodes, the following Equation can be obtained [29]:
l i j = d i j m d i j I x i I x j I d i j I Δ x i I y i I y j I d i j I Δ y i I z i I z j I d i j I Δ z i I + x i I x j I d i j I Δ x j I + y i I y j I d i j I Δ y j I + z i I z j I d i j I Δ z j I + Δ i j
where d i j m is the measured distance between nodes i and j, d i j I is the inertial distance between nodes i and j, l i j is the difference between these two distances, { x i I , y i I , z i I } is the INS output position of node i and { Δ x i I , Δ y i I , Δ z i I } is the position error of node i; it is assumed that the error of distance difference l i j obeys a normal distribution. This assumption is based on the fact that the noise of the INS and the error of the measured distance between nodes obey a normal distribution. Δ i j is the error of distance difference l i j , Δ i j ~ N ( 0 ,   σ r 2 ) and σ r 2 is the variance of the normal distribution [29].
The matrix form of (12) is [30]:
{ L d = T S + Δ E ( L d ) = 0 L d = σ 0 2 P c o 1
where L d is the p-dimensional measurement vector, T is the p × 3a dimensional coefficient matrix (a is the number of nodes), S is the 3a-dimensional position error vector, Δ is the p-dimensional measurement error vector, σ 0 2 is the unit weight variance and P c o is the p × p dimensional measurement error covariance matrix of L d .
According to the principle of least squares, Δ and S are replaced with their respective estimated values V n and S ^ . The model of the network navigation algorithm is established as follows [29]:
{ V n = T S ^ L d V n T P c o V n = min
The above formula can be rewritten as follows:
S ^ = ( T T P c o T ) 1 T T P c o L d
where S ^ is the estimated position error vector.
However, the above formula can only prove that the shapes of the measurement network with measured distances between nodes as edges and the inertial position network with inertial distances calculated from the INS output as edges are the same.
To make the absolute positions of the two networks consistent in three-dimensional space, it is necessary to rotate and translate the two networks through the principle of the free network theory with rank deficiency mentioned in [29].
According to this principle, the above rotation and translation process is equivalent to adding the following new solution conditions [29]:
G T P X S ^ = 0
where P X is a 3a × 3a dimensional weight matrix determined by the accuracy level of each node’s inertial navigation system. The expression of P X is as follows:
P X = [ p x 1 0 0 0 0 0 0 p y 1 0 0 0 0 0 0 p z 1 0 0 0 0 0 0 p x n 0 0 0 0 0 0 p y n 0 0 0 0 0 0 p z n ]
where G T is a matrix composed of six linearly independent eigenvectors corresponding to the zero eigenvalue of T T P c o T . G T satisfies the following conditions:
{ T G = 0 T T P c o T G = 0 rank ( G T ) = 6
The first item in (18) can ensure that the added solution conditions (16) and the error Equation (15) are independent of each other. It can be seen from the second item in (18) that G is composed of the zero eigenvectors of T T P c o T . The third item in (18) means that G supplements the error Equation (15) with six constraints caused by rotation and translation [29].
The final estimated error of the node’s INS output can be expressed as:
S ^ = ( T T P c o T + P X G G T P X ) 1 T T P c o L d
According to the position S output by the INS and the error S ^ estimated by the network navigation algorithm, the corrected value S ˜ of the node position can be obtained:
S ˜ = S S ^
The theory of the network navigation algorithm has been well established, and the specific details can be found in [30].

4. Simulation Evaluation

4.1. Simulation Configuration

The simulation uses four UAVs carrying the GNSS and INS, and employs the data link method to obtain the distances among all UAVs [31].
At the beginning, four UAVs formed a square swarm with sides of 60 km. Subsequently, the swarm moved north and then east at a speed of 200 m/s. The true trajectory of the UAV swarm is shown in Figure 2.
The velocity of the four UAVs in the east–north–up navigation reference frame is shown in Figure 3. (the change at about 220 s was caused by the movement of the UAV).
The initial attitude error was ( 0.5   0.5   1.5 ) , the initial velocity error was ( 0.5   m / s   0.5   m / s   0.5   m / s ) and initial position error was ( 20   m   20   m   20   m ) . The bias and standard deviation of the white noise of the accelerometer were 800   μ g and 100   μ g s . The bias and standard deviation of white noise of the gyro were 3 ° / h and 1 ° / h . The standard deviation of the measurement error vector was 20 m. The process noise covariance matrix Q was set as:
Q = d i a g [ 0 3 × 3 , ( 100   μ g s ) 2 · I 3 × 3 , ( 1 / h ) 2 · I 3 × 3 , 0 6 × 6 ]

4.2. Evaluation of Dynamic Adaptive Kalman Filter

In order to highlight the superiority of the proposed algorithm, two simulations were conducted.
In the first simulation, the noise of the GNSS signal was designed as a step from 10 m to 30 m with the purpose of verifying the adaptability of the proposed algorithm. The parameters of EEMD were set as follows: M = 20 and σ n = 0.2 .
It can be seen from Figure 4 that the proposed method could follow the changes of noise in three direction channels well. It is worth mentioning that in order to verify the effect of the proposed algorithm on noise reconstruction, this paper compares the standard deviation of the reconstructed noise and the real noise in the same window. Although the standard deviation of the real noise was set to three constant values, the real noise was actually a random value fluctuating around this constant value, due to the randomness of the noise itself [32]. Therefore, the green line is actually the standard deviation of the real noise calculated by the window at different positions. In addition, any digital signal processing method requires a window size for the accurate extraction of noise; therefore, errors in the first four seconds are allowed ( w = 4   s ).
The average difference between the two curves in Figure 4 is ( X :   1.450   m   Y :   0.574   m   Z :   0.568   m )
The second simulation compared the classic Kalman filter with the proposed algorithm to prove the effectiveness of the DAKF. This simulation used the trajectory of UAV1 in Figure 2, and this method was also verified on the other three trajectories.
Experimental parameters are shown in Table 1. In the first four seconds, the observation noise matrix R of DAKF was consistent with the KF, due to the introduction of the window in the proposed algorithm; in the following time, R of DAKF was dynamically obtained by EEMD and R of the KF was still a constant empirical value. In order to verify the universality of the proposed algorithm, setting the real noise of the three channels of GNSS included three stages: the real noise less than the model setting, the real noise greater than the model setting, and the real noise equal to the model setting. To ensure the effectiveness of the simulation, in Table 1, the real noise of the DAKF and KF were consistent. In other words, the data in the second row belong to the second and third columns together in Table 1.
Figure 5 shows the change of the scaling factor. The fluctuation of the scaling factor at the initial moment is a normal phenomenon caused by the initial parameter setting which has been proven [10]. At 200 s and 400 s, it was used to prove the normal operation of the adaptive part that the scaling factor following the change of the GNSS noise [33]. Except for the above three transitional periods, the scaling factor was well maintained around one.
Figure 6 shows the simulation results of the KF and DAKF with dynamic observation noise. When the model noise did not match the real noise, the DAKF performed significantly better than the KF. When equal, the effects of the DAKF and KF were similar. A further explanation of Figure 6 is given by Table 2. Table 2 shows the average position error of the three channels, which proves the effective improvement of the DAKF.

4.3. Evaluation of Networked Navigation Algorithm

In order to fully verify the performance of the NNA, the UAV movement during the GNSS interruption should include two modes: straight and turning. The total time of the UAV swarm flying was designed to be 520 s while there was a GNSS outage of 150 s start at the 100th second.
Since the swarm had no displacement in the vertical direction, Figure 7 shows the swarm trajectory in the horizontal direction in order to better show the effect of the NNA. The purple line is the navigation results of the INS. The yellow line is the trajectory after the position compensation of the NNA. The gray line is the output of the UAV’s INS/GNSS integrated system, which was basically consistent with the true trajectory, when the GNSS was available. The black line is the true trajectory shown in Figure 2. It can be easily obtained from Figure 7 that the NNA could obviously improve the navigation accuracy of the overall UAV swarm, no matter whether in the straight or turning mode.
Table 3 describes the improvement in positioning accuracy of the UAV swarm after using the NNA. Δ p represents the position error of four UAV’s inertial navigation system and Δ p ˜ represents the position error after the compensation of the NNA. In this simulation, the above position error was calculated from the last moment of the GNSS signal being unavailable. The expression is as follows:
Δ p i = ( Δ x i I ) 2 + ( Δ y i I ) 2 + ( Δ z i I ) 2 Δ p ˜ i = ( Δ x ˜ i I ) 2 + ( Δ y ˜ i I ) 2 + ( Δ z ˜ i I ) 2 i = 1 , 2 , 3 , 4
where Δ x ˜ i I , Δ y ˜ i I and Δ z ˜ i I are the position errors in three directions after the NNA compensation of the ith aircraft.
To define the improvement of positioning accuracy as the ratio of Δ p and Δ p ˜ :
P i m p = Δ p Δ p ˜
Table 3 shows that the NNA could increase the positioning accuracy of the UAV swarm by 1.93 times. In addition, the purpose of the NNA was to improve the positioning accuracy of the entire swarm. The deterioration of the positioning accuracy of a single node (UAV3) due to the change of the swarm’s shape which was caused by the divergence of the INS navigation errors with time was allowed [19,34].
In order to verify the superiority of the proposed algorithm compared with the neural network-based algorithm, the NNA and MLP-based Bagging method mentioned in [18] were compared on UAV1 (the simulation result was similar on the other three UAVs). The parameter settings and experimental results of the NNA were consistent with the above simulation. The topology and transfer function of the MLP were 12 (input layer) × 6 (hidden layer) × 1 (output layer) and logsig-linear, respectively. The network number of Bagging was set to 50. The training data of the neural network were UAV information during the first 100 s when the GNSS signal was available. The output of the MLP was the increment of the UAV’s position during the GNSS outage.
Figure 8 shows the comparison result between the NNA and MLP-based Bagging method. The pink line is the prediction result of the MLP-based Bagging method.
The positioning error of the MLP-based Bagging was 6385.014 km, which was approximately 41 times larger than the positioning error of the NNA. One of the disadvantages of the neural network-based algorithm is that it relies heavily on training data. If the training sample and the test sample were significantly different, the above situation would be an inevitable phenomenon [18]. In contrast, the NNA was only related to the distance between nodes at the current moment, so it would not be affected by historical data. It can be seen from the simulation results that the NNA was more robust than the neural network-based methods.

5. Conclusions

This paper proposed a new navigation scheme for the UAV swarm. The proposed algorithm included the following two cases: (1) By adaptively acquiring the observation noise covariance matrix using EEMD and adjusting the process noise covariance matrix, a dynamic adaptive Kalman filter which is superior to the traditional Kalman filter was designed when the GNSS signal was available. The simulation shows that the DAKF could effectively improve the accuracy of the Kalman filter. (2) In the phase of GNSS signal outage, the positioning error of the UAV swarm was reduced through the use of the networked navigation algorithm. The simulation demonstrated that the NNA could improve the positioning accuracy of the swarm by 93% during the GNSS interruption of 150 s and better than the neural network-based algorithm.
Future work mainly includes the following research directions. First, the colored noise and non-Gaussian noise were not considered in this simulation. It would determine whether the proposed algorithm has a sufficient generality in all environments. Second, the window size of the DAKF was an empirical value obtained from a large number of experiments. How to decrease its impact on the algorithm will be considered in the next step. Finally, the proposed algorithm has not been verified on an actual UAV swarm, which will be carried out in future work.

Author Contributions

Conceptualization, J.Z., X.W. and W.Z.; methodology, X.W. and W.Z.; validation, X.W. and W.Z.; writing—original draft preparation, W.Z.; writing—review and editing, X.W. and W.Z.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Beijing Natural Science Funds (grant number: 4204103) and the Aeronautical Science Fund (grant number: 2019ZC051009).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, M.; Xiong, Z.; Liu, J.; Wang, R.; Xiong, J. Cooperative navigation of unmanned aerial vehicle swarm based on cooperative dilution of precision. Int. J. Adv. Robot. Syst. 2020, 17, 1729881420932717. [Google Scholar] [CrossRef]
  2. Wang, R.; Du, J.; Xiong, Z.; Chen, X.; Liu, J. Hierarchical Collaborative Navigation Method for UAV Swarm. J. Aerosp. Eng. 2021, 34, 4020097. [Google Scholar] [CrossRef]
  3. Sabzevari, D.; Chatraei, A. INS/GPS Sensor Fusion based on Adaptive Fuzzy EKF with Sensitivity to Disturbances. IET Radar Sonar Navig. 2021. [Google Scholar] [CrossRef]
  4. Valade, A.; Acco, P.; Grabolosa, P.; Fourniols, J.-Y. A study about Kalman filters applied to embedded sensors. Sensors 2017, 17, 2810. [Google Scholar] [CrossRef] [Green Version]
  5. Pan, C.; Gao, J.; Li, Z.; Qian, N.; Li, F. Multiple fading factors-based strong tracking variational bayesian adaptive Kalman filter. Measurement 2021, 176, 109139. [Google Scholar] [CrossRef]
  6. Huang, Y.; Zhang, Y.; Wu, Z.; Li, N.; Chambers, J. A novel adaptive Kalman filter with inaccurate process and measurement noise covariance matrices. IEEE Trans. Autom. Control 2017, 63, 594–601. [Google Scholar] [CrossRef] [Green Version]
  7. Liu, Y.; Fan, X.; Lv, C.; Wu, J.; Li, L.; Ding, D. An innovative information fusion method with adaptive Kalman filter for integrated INS/GPS navigation of autonomous vehicles. Mech. Syst. Signal Process. 2018, 100, 605–616. [Google Scholar] [CrossRef] [Green Version]
  8. Li, N.; Kang, Z.W.; Liu, J.; Xu, X.M. An adaptive filtering method based on EMD for X-ray pulsar navigation with uncertain measurement noise. MATEC Web Conf. EDP Sci. 2017, 114, 4017. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, L.; Wang, S.; Selezneva, M.S.; Neusypin, K.A. A new adaptive Kalman filter for navigation systems of carrier-based aircraft. Chin. J. Aeronaut. 2021. [Google Scholar] [CrossRef]
  10. Ding, W.; Wang, J.; Rizos, C.; Kinlyside, D. Improving adaptive Kalman estimation in GPS/INS integration. J. Navig. 2007, 60, 517. [Google Scholar] [CrossRef] [Green Version]
  11. Youn, W.; Ko, N.Y.; Gadsden, S.A.; Myung, H. A novel multiple-model adaptive Kalman filter for an unknown measurement loss probability. IEEE Trans. Instrum. Meas. 2020, 70, 1–11. [Google Scholar] [CrossRef]
  12. Nourmohammadi, H.; Keighobadi, J. Fuzzy adaptive integration scheme for low-cost SINS/GPS navigation system. Mech. Syst. Signal Process. 2018, 99, 434–449. [Google Scholar] [CrossRef]
  13. Zhou, Y.; Wan, J.; Li, Z.; Song, Z. GPS/INS integrated navigation with BP neural network and Kalman filter. In Proceedings of the 2017 IEEE International Conference on Robotics and Biomimetics (ROBIO), Macau, Macao, 5–8 December 2017; pp. 2515–2520. [Google Scholar]
  14. Li, D.; Jia, X.; Zhao, J. A novel hybrid fusion algorithm for low-cost GPS/INS integrated navigation system during GPS outages. IEEE Access. 2020, 8, 53984–53996. [Google Scholar] [CrossRef]
  15. Wei, X.; Li, J.; Feng, K.; Zhang, D.; Li, P.; Zhao, L.; Jiao, Y. A Mixed Optimization Method Based on Adaptive Kalman Filter and Wavelet Neural Network for INS/GPS During GPS Outages. IEEE Access. 2021, 9, 47875–47886. [Google Scholar] [CrossRef]
  16. Chen, X.; Shen, C.; Zhang, W.; Tomizuka, M.; Xu, Y.; Chiu, K. Novel hybrid of strong tracking Kalman filter and wavelet neural network for GPS/INS during GPS outages. Measurement 2013, 46, 3847–3854. [Google Scholar] [CrossRef]
  17. Abdolkarimi, E.S.; Abaei, G.; Mosavi, M.R. A wavelet-extreme learning machine for low-cost INS/GPS navigation system in high-speed applications. GPS Solut. 2018, 22, 15. [Google Scholar] [CrossRef]
  18. Zhang, C.; Guo, C.; Guo, M.-Z. Information fusion based on artificial intelligence method for SINS/GPS integrated navigation of marine vessel. J. Electr. Eng. Technol. 2020, 15, 1345–1356. [Google Scholar] [CrossRef]
  19. Peide, L.J.Z.J. Swarming aircraft collaborative localization based on mutual rangings. J. Beijing Univ. Aeronaut. Astronaut. 2012, 38, 541–545. [Google Scholar]
  20. Wu, Z.; Huang, N.E. Ensemble empirical mode decomposition: A noise-assisted data analysis method. Adv. Adapt. Data Anal. 2009, 1, 1–41. [Google Scholar] [CrossRef]
  21. Liu, C.; Yang, Z.; Shi, Z.; Ma, J.; Cao, J. A gyroscope signal denoising method based on empirical mode decomposition and signal reconstruction. Sensors 2019, 19, 5064. [Google Scholar] [CrossRef] [Green Version]
  22. Park, S.; Gil, M.-S.; Im, H.; Moon, Y.-S. Measurement noise recommendation for efficient Kalman filtering over a large amount of sensor data. Sensors 2019, 19, 1168. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.; Yen, N.-C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  24. Zhao, Q.; Gao, W.; Xiang, W.; Shi, R.; Liu, C.; Zhai, T.; Huang, H.A.; Gumley, L.E.; Strabala, K. Analysis of air quality variability in Shanghai using AOD and API data in the recent decade. Front. Earth Sci. 2013, 7, 159–168. [Google Scholar] [CrossRef]
  25. Zhou, Y.; Zhou, S.-T.; Zhong, Z.-Y.; Li, H.-G. A de-illumination scheme for face recognition based on fast decomposition and detail feature fusion. Opt. Express 2013, 21, 11294–11308. [Google Scholar] [CrossRef] [PubMed]
  26. Lei, Y.; He, Z.; Zi, Y. EEMD method and WNN for fault diagnosis of locomotive roller bearings. Expert Syst. Appl. 2011, 38, 7334–7341. [Google Scholar] [CrossRef]
  27. Flandrin, P.; Rilling, G.; Goncalves, P. Empirical mode decomposition as a filter bank. IEEE Signal Process. Lett. 2004, 11, 112–114. [Google Scholar] [CrossRef] [Green Version]
  28. Grewal, M.S.; Andrews, A.P. Kalman Filtering: Theory and Practice with MATLAB; John Wiley & Sons: New York, NY, USA, 2014; ISBN 1118984919. [Google Scholar]
  29. Cui, X. General Adjustment; Wuhan University Press: Wuhan, China, 2009; pp. 20–190. (In Chinese) [Google Scholar]
  30. Zhang, J.; Lei, H.; Wang, X.; Qin, F. Research on Navigation Error Estimation Technology Based on Distance Measurement Information of Aircraft Three-Dimensional Network. Air Space Def. 2020, 3, 124–130. (In Chinese) [Google Scholar]
  31. Wang, X.G.; Guo, J.F.; Cui, N.G. Cooperative localization approach to intelligent missile based on data link. J. Chinese Inert. Technol. 2009, 17, 319–323. [Google Scholar]
  32. Rice, S.O. Mathematical analysis of random noise. Bell Syst. Tech. J. 1944, 23, 282–332. [Google Scholar] [CrossRef]
  33. Xu, T.-L.; You, W.-H.; Cui, P.-Y. Research on GPS/INS integrated navigation system based on fuzzy adaptive Kalman filtering. Yuhang Xuebao/J. Astronaut. 2005, 26, 571–575. (In Chinese) [Google Scholar]
  34. Gao, J.; Li, K.; Chen, Y. Study on integration of FOG single-axis rotational INS and odometer for land vehicle. IEEE Sens. J. 2017, 18, 752–763. [Google Scholar] [CrossRef]
Figure 1. Process of Dynamic Adaptive Kalman Filter.
Figure 1. Process of Dynamic Adaptive Kalman Filter.
Sensors 21 05374 g001
Figure 2. True trajectory of the UAVs.
Figure 2. True trajectory of the UAVs.
Sensors 21 05374 g002
Figure 3. Velocity of UAV.
Figure 3. Velocity of UAV.
Sensors 21 05374 g003
Figure 4. Proof of algorithm following effect. The green curve represents standard deviation (STD) of real noise, the blue curve represents STD estimated by EEMD.
Figure 4. Proof of algorithm following effect. The green curve represents standard deviation (STD) of real noise, the blue curve represents STD estimated by EEMD.
Sensors 21 05374 g004
Figure 5. The change of scaling factor.
Figure 5. The change of scaling factor.
Sensors 21 05374 g005
Figure 6. Position error of the two algorithms.
Figure 6. Position error of the two algorithms.
Sensors 21 05374 g006
Figure 7. Swarm trajectory in the horizontal direction.
Figure 7. Swarm trajectory in the horizontal direction.
Sensors 21 05374 g007
Figure 8. Comparison results between NNA and MLP-based Bagging method.
Figure 8. Comparison results between NNA and MLP-based Bagging method.
Sensors 21 05374 g008
Table 1. Experimental parameters.
Table 1. Experimental parameters.
DAKFKF
R m o d e l / m Dynamic(30,30,30)
R e a l   n o i s e / m (2,2,2)→(100,100,100)→(30,30,30)
Table 2. Average position error.
Table 2. Average position error.
DAKFKFImprovement
X/m24.170252.423554%
Y/m48.825174.259434%
Z/m49.288374.155134%
Table 3. Improvement in positioning accuracy of NNA.
Table 3. Improvement in positioning accuracy of NNA.
Δ p ˜ / k m Δ p / k m P i m p
UAV 1157.142460.4432.93
UAV 2102.436224.7732.20
UAV 3127.55596.9230.76
UAV 4177.985307.2561.73
Average141.280272.3491.93
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, J.; Zhou, W.; Wang, X. UAV Swarm Navigation Using Dynamic Adaptive Kalman Filter and Network Navigation. Sensors 2021, 21, 5374. https://doi.org/10.3390/s21165374

AMA Style

Zhang J, Zhou W, Wang X. UAV Swarm Navigation Using Dynamic Adaptive Kalman Filter and Network Navigation. Sensors. 2021; 21(16):5374. https://doi.org/10.3390/s21165374

Chicago/Turabian Style

Zhang, Jingjuan, Wenxiang Zhou, and Xueyun Wang. 2021. "UAV Swarm Navigation Using Dynamic Adaptive Kalman Filter and Network Navigation" Sensors 21, no. 16: 5374. https://doi.org/10.3390/s21165374

APA Style

Zhang, J., Zhou, W., & Wang, X. (2021). UAV Swarm Navigation Using Dynamic Adaptive Kalman Filter and Network Navigation. Sensors, 21(16), 5374. https://doi.org/10.3390/s21165374

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