Next Article in Journal
A Partial Ambiguity Resolution Algorithm Based on New-Breadth-First Lattice Search in High-Dimension Situations
Next Article in Special Issue
Adverse Weather Target Detection Algorithm Based on Adaptive Color Levels and Improved YOLOv5
Previous Article in Journal
A Comparison among Different Strategies to Detect Potential Unstable Behaviors in Postural Sway
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cooperative Location Method for Leader-Follower UAV Formation Based on Follower UAV’s Moving Vector

College of Automation Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(19), 7125; https://doi.org/10.3390/s22197125
Submission received: 22 July 2022 / Revised: 14 September 2022 / Accepted: 16 September 2022 / Published: 20 September 2022

Abstract

:
The traditional leader-follower Unmanned Aerial Vehicle (UAV) formation cooperative positioning (CP) algorithm, based on relative ranging, requires at least four leader UAV positions to be known accurately, using relative distance with leader UAVs to achieve the unknown position follower UAV’s high-precision positioning. When the number of the known position leader UAVs is limited, the traditional CP algorithm is not applicable. Aiming at the minimum cooperative unit, which consists of a known position leader UAV and an unknown position follower UAV, this paper proposes a CP method based on the follower UAV’s moving vector. Considering the follower UAV can only acquire the single distance with the leader UAV at each distance-sampling period, it is difficult to determine the follower UAV’s spatial location. The follower UAV’s moving vector is used to construct position observation of the follower UAV’s inertial navigation system (INS). High-precision positioning is achieved by combining the follower UAV’s moving vector. In the process of CP, the leader UAV obtains a high-precision position by an INS/Global Positioning System (GPS) loosely integrated navigation system and transmits its position information to the follower UAV. Based on accurate modeling of the follower UAV’s INS, the position, velocity and heading observation equation of the follower UAV’s INS are constructed. The improved extended Kalman filtering is designed to estimate the state vector to improve the follower UAV’s positioning accuracy. In addition, considering that the datalink system based on radio signals may be interfered with by the external environment, it is difficult for the follower UAV to obtain relative distance information from the leader UAV in real time. In this paper, the availability of the relative distance information is judged by a two-state Markov chain. Finally, a real flight test is conducted to validate the performance of the proposed algorithm.

1. Introduction

Over the past few decades, UAVs have been widely used in military and civilian fields. In the civilian field, they are mainly used for rescue missions [1], border patrol [2] and formation flying [3]. In the military field, they are mainly used for coordinated strikes and coordinated tracking [4]. However, due to the constraints of the volume and load of UAVs, a single UAV has low success rates in performing tasks. Multiple UAVs formation can make full use of their resource advantage, and the efficiency of performing tasks is much greater than a single UAV. However, multiple UAVs’ autonomous flight needs to solve two basic problems: firstly, solving the problem of UAVs’ real-time navigation and positioning. Secondly, solving the problem of UAVs’ real-time tracking control, mission and trajectory planning. The high-precision navigation information is the input of the UAV flight control system, so how to accurately acquire the UAVs’ high-precision navigation information is the key for a UAV formation to complete a flight mission.
INS is an autonomous, independent dead-reckoning system, that has been widely used in UAVs. However, the navigation accuracy is difficult to guarantee when only using the INS due to inertial measurement unit’s (IMU) error accumulation. When the satellite signals are available, the GPS receiver can provide the long-term reliable position information for the UAV. Thus, the INS/GPS loosely integrated navigation system is widely used in UAV formation [5,6,7]. In addition, the original information output from the satellite receiver, such as pseudo-range and carrier phase, can be used as system observation to correct INS’s long-term cumulative error. W Khalaf proposed an INS/GPS adaptive, tightly integrated method by fusing the GPS original pseudo-range, magnetometer and barometric altimeter information [8]. S Huang proposed an INS/GPS integrated navigation method based on GPS dual antennas. In this method, the two-difference calculation is performed by carrier phase measurement. The two-difference observation equation is established to improve the UAV’s attitude accuracy [9].
However, satellite signals may be interfered with by the external environment. Abundant external relative perception information among UAVs can effectively improve UAVs’ positioning accuracy. The airborne datalink system, as a communication system based on radio data transmission and receiving, is widely used in the field of collaborative navigation [10]. Under the condition of a UAV cooperative network, the high-precision time reference is established through a datalink time precision alignment method [11]. The relative distance, azimuth and pitch information between UAVs can be provided in real time by Time of Arrival (TOA) and Angle of Arrival (AOA) methods [12]. To the best of our knowledge, the datalink system can provide a relative ranging accuracy of about 30m within the 550 km communication range, under the condition of a time-synchronized cooperative network. Aiming at the relative ranging mode of the datalink system, many researchers have carried out a lot of works to realize multiple UAVs’ high-precision CP. Yang proposed a multi-UAVs collaborative navigation method based on relative distance and magnetic sensor measurement [13]. In this method, the UAV’s kinematic model is constructed combining the output of a magnetic sensor. The accumulated error of the dead-reckoning system is corrected using distance observation. Qu Y [14] proposed a CP method of low-cost UAV using relative range measurements in multi-UAVs flight. Under the condition that at least four leader UAV positions are known in the UAV formation, the relative distance with the leader UAVs is used as the system observation to improve the unknown position follower UAV’s positioning accuracy. In addition, Chen M [15] proposed a cooperative navigation method for a UAV swarm based on cooperative dilution of precision. The spatial configuration between leader UAVs and follower UAVs is considered, using relative distance with four known position leader UAVs in the optimal spatial geometry configuration as the observation method to further improve follower UAVs’ positioning accuracy. Through the above analysis, the current mainstream CP algorithm based on relative distance mainly adopts the multiple leader UAV CP framework, which is not applicable to the minimum cooperative unit composed of two UAVs.
Aiming at the minimum cooperative unit, which consists of a known position leader UAV and an unknown position follower UAV, there is little literature on only using the single relative distance to achieve follower UAV high-precision positioning, but there are a few applications in Autonomous Underwater Vehicles (AUVs). Considering the single relative distance between the AUV and a known position beacon, it is difficult to determine an AUV’s spatial position. Many researchers have constructed the virtual beacon arrays as known position beacons to achieve the AUV’s positioning [16,17]. However, the position of the virtual beacon is affected by the beacon calibration error, and the AUV’s velocity error and heading error affect the AUV’s positioning accuracy. In addition, by combining the navigation sensor mounted on the AUV, such as the pressure sensor, Doppler Velocity Log (DVL) and geomagnetic sensor, etc., many researchers have constructed the AUV’s motion model based on motion characteristics. The relative distance or signal propagation time between the known position beacon and AUV is used as the system observation to suppress positioning error divergence of the INS/DVL/pressure-sensor-integrated navigation system. However, due to the lack of sufficient position observation, it can inhibit positioning error accumulation of the integrated navigation system to a certain extent by using the above method. The positioning error of the integrated navigation system is still divergent under long-term conditions [18,19,20,21,22].
Thus, aiming at the minimum cooperative unit, this paper proposes a CP method based on the follower UAV’s moving vector. The error model of the follower UAV’s INS is constructed according to the IMU error characteristic. By combining the output of the speed sensor and magnetic sensor mounted on the follower UAV, the follower UAV’s moving vector is used to construct a position observation equation of the follower UAV’s INS. The improved extended Kalman filter is utilized to estimate the state vector to obtain a follower UAV high-precision navigation solution. In addition, considering that the airborne datalink system based on radio signals may be interfered with by the UAV formation-flying environment, the follower UAV cannot obtain relative distance information with the leader UAV in real time. In this paper, a two-state Markov chain is used to judge the availability of the distance information to improve the efficiency of the CP algorithm. The innovations of this paper are as follows:
  • Aiming at the minimum cooperative unit, which consists of a known position leader UAV and an unknown position follower UAV, this paper proposes a CP algorithm based on follower the UAV’s moving vector.
  • High-precision positioning of the follower UAV is achieved by using the single relative distance information and the follower UAV’s moving vector.
  • Introducing a two-state Markov chain to judge the availability of the distance information improves the efficiency of the cooperative navigation algorithm.
The rest of the paper is organized as follows: In the Section 2, the CP framework and the concept of the follower UAV’s moving vector are introduced. Observation equations based on the follower UAV’s moving vector are established in Section 3. The CP algorithm based on improved extended Kalman filtering is designed in Section 4. In Section 5, the comparative analysis is made between the proposed algorithm and traditional multiple leader UAV CP algorithm based on relative ranging. Finally, the conclusion of this study is summarized in Section 6.

2. Overview of the Proposed Cooperative Localization Algorithm

2.1. The Framework of Proposed Cooperative Localization Algorithm

The CP scenario of the minimum cooperative unit based on relative ranging is shown in Figure 1. In this paper, the geographic coordinate system is defined as the navigation coordinate system. In the navigation coordinate system, with the x-axis pointing in the east direction and the y-axis pointing in the north direction, the z-axis is determined by right-hand rule, and the origin is located at the centroid of the carrier. The front and right directions in the body coordinate system are pointing to the x-axis, y-axis, respectively, and the z-axis is determined by the right-hand criterion.
In the above CP scenario, the known position leader UAV flies in an area with available GPS signals and is equipped with a low-precision Micro-Electro-Mechanical System (MEMS) IMU, a GPS receiver in a single-point mode and a datalink communication device. The unknown position follower UAV flies in an area denied GPS signals and carries a low-precision MEMS-IMU, speed sensor, magnetic sensor, barometric altimeter and datalink communication device. The two UAVs perform time calibration by GPS system clock before cooperative flight to ensure each UAV’s time reference is synchronized. In the process of CP, the leader UAV obtains high-precision navigation information by fusing INS/GPS information. Meanwhile, the leader UAV broadcasts its position, speed and course information to the follower UAV node. According to the preset time interval period, the follower UAV receives navigation information transmitted by the leader UAV through the airborne datalink system. Under the condition that the clock of each UAV node in the minimum cooperative unit is synchronized, the relative distance between the leader and follower UAV is obtained in real time by the TOA method. Considering the positioning error of the follower UAV’s INS has accumulated and diverged for a long time, only relying on the positioning accuracy of the follower UAV’s INS is difficult for meeting the UAV cooperative flight requirement. By combining the altitude information obtained from the follower UAV’s barometric altimeter, three-dimensional navigation positioning of the follower UAV is transformed into two dimensions. The moving vector of the follower UAV at adjacent times is constructed using the information of the speed sensor and three-axis magnetic sensor. Based on the error model of the follower UAV’s INS being accurately constructed, the position observation equation of the follower UAV’s INS is constructed using the relative ranging constraint and the follower UAV’s moving vector. Meanwhile, combining the speed and heading information output from the speed and three-axis magnetic sensor, the velocity and heading observation equation of the follower UAV’s INS are constructed, respectively. In addition, considering that the datalink system based on radio signals may be interfered with by the external environment, this paper introduces a two-state Markov chain to judge the availability of the relative distance information. When the distance information among the UAVs is available, the position observation equation of follower the UAV’s INS is constructed by using the follower UAV’s moving vector. Otherwise, only the velocity and heading information obtained from the follower UAV’s airborne sensor are used as the system observation vector. Finally, the improved extended Kalman filtering is designed to estimate the state vector to improve the follower UAV’s positioning accuracy. The schematic diagram of the proposed method in this paper is shown in Figure 2.

2.2. Moving Vector of Follower UAV

Considering there is only one known position leader UAV in the minimum cooperative unit, the single relative distance with the leader UAV cannot determine the follower UAV’s spatial position. The moving vector of the follower UAV has to play an important role in constructing the position observation equation of the follower UAV’s INS. The follower UAV’s moving vector represents the relative displacement vector that the follower UAV removed from time t k to time t k + 1 . The schematic diagram of the follower UAV’s moving vector in the horizontal direction is illustrated in Figure 3. The relative position constraint of the follower UAV from time t k to time t k + 1 in the horizontal direction can be expressed as:
x k + 1 y k + 1 = 1 0 0 1 x k y k + 1 0 0 1 D x k , k + 1 D y k , k + 1 = 1 0 0 1 x k y k + v ^ cos θ ^ cos φ ^ Δ t v ^ cos θ ^ sin φ ^ Δ t
where D x k , k + 1 represents the moving vector of the follower UAV in X direction. D y k , k + 1 represents the moving vector of the follower UAV in Y direction. ( x k + 1 , y k + 1 ) is the horizontal position of the follower UAV at time t k + 1 . ( x k , y k ) is the horizontal position of the follower UAV at time t k . v ^ is the follower UAV’s velocity obtained from the speed sensor. θ ^ is the follower UAV’s pitch angle. φ ^ is the heading angle obtained from the follower UAV’s three-axis magnetic sensor. Δ t is the time interval from time t k to t k + 1 . Considering the existence of sensor measurement error, the follower UAV’s velocity, pitch angle and heading angle have the following form:
v ^ = v ¯ + Δ v θ ^ = θ ¯ + Δ θ φ ^ = φ ¯ + Δ φ
where ( v ¯ , θ ¯ , φ ¯ ) is the measured true value of velocity, pitch angle and heading angle, respectively, and ( Δ v , Δ θ , Δ φ ) is the velocity error, pitch angle error and heading error, respectively. Substituting Formula (2) into Formula (1) and ignoring the small error term, the relative displacement error of the follower UAV from time t k to t k + 1 can be expressed as:
Δ D x k , k + 1 Δ D y k , k + 1 = v ¯ Δ t ( cos θ ¯ sin φ ¯ Δ φ + sin θ ¯ cos φ ¯ Δ θ ) + Δ v Δ t ( cos θ ¯ cos φ ¯ cos θ ¯ sin φ ¯ Δ φ sin θ ¯ cos φ ¯ Δ θ ) v ¯ Δ t ( cos θ ¯ cos φ ¯ Δ φ sin θ ¯ sin φ ¯ Δ θ ) + Δ v Δ t ( cos θ ¯ sin φ ¯ + cos θ ¯ cos φ ¯ Δ φ sin θ ¯ sin φ ¯ Δ θ )
It can be seen from Formula (3) that the moving vector of the follower UAV from time t k to time t k + 1 in the horizontal direction is jointly affected by the follower UAV’s velocity error, pitch angle error and heading error. Considering the relatively small change of the pitch angle at adjacent times during UAV flight, the relative displacement error of the follower UAV at adjacent times is mainly affected by the velocity error and heading error. Generally, the airborne speed sensor can accurately measure the follower UAV’s velocity, while the three-axis magnetometer is vulnerable to interfere from the external environment. The three-axis magnetometer with temperature and magnetic compensation in advance can obtain higher heading measurement accuracy. The relative displacement accuracy of the follower UAV at adjacent times will affect the follower UAV’s positioning accuracy. The follower UAV’s positioning accuracy can be effectively improved by reducing the relative displacement error of the follower UAV at adjacent times. The smaller relative displacement error of the follower UAV at adjacent times, the higher the follower UAV’s positioning accuracy.

3. System Observation Equation Based on Follower UAV’s Moving Vector

Assuming that each UAV’s clock completes its time synchronization by GPS before UAV formation cooperative flight, in the process of CP, the leader UAV sends its pulse information and position broadcast through the airborne datalink system at time t k . Assuming that the ( x k m , y k m , z k m ) and ( x k , y k , z k ) represent the positions of the leader and follower UAV at time t k , respectively, the follower UAV receives a pulse signal from the leader UAV and measures relative distance d k . Considering the barometric altimeter mounted on the follower UAV can be used to measure altitude information in real time, based on the high-precision position information of the leader UAV and the height information of the follower UAV, the follower UAV’s position solution can be converted from three-dimensional to two-dimensional. The spatial position relationship in the horizontal direction between the leader and follower UAV at the adjacent times is illustrated in Figure 4. We can conclude from Figure 4 that at time t k , the follower UAV is located on a circle with ( x k m , y k m ) as the center and r k as the radius. r k can be expressed as:
r k = d k 2 ( z k m z k ) 2
The circular equation at time t k can be expressed as:
( x k x k m ) 2 + ( y k y k m ) 2 = r k 2
where d k is the three-dimensional relative distance between the leader and follower UAV at time t k . z k m is the altitude information at time t k obtained from the leader UAV’s INS/GPS loosely integrated navigation system. z k is the altitude information obtained from the follower UAV’s onboard barometric altimeter.
At time t k + 1 , the follower UAV is located on a circle with ( x k + 1 m , y k + 1 m ) as the center and r k + 1 as the radius. r k + 1 can be expressed as:
r k + 1 = d k + 1 2 ( z k + 1 m z k + 1 ) 2
The circular equation at time t k + 1 can be expressed as:
( x k + 1 x k + 1 m ) 2 + ( y k + 1 y k + 1 m ) 2 = r k + 1 2
where d k + 1 is the three-dimensional relative distance information between the leader UAV and the UAV at time t k + 1 . z k + 1 m is the altitude information at time t k + 1 obtained from the leader UAV’s INS/GPS loosely integrated navigation system. z k + 1 is the altitude information at time t k + 1 obtained from the follower UAV’s onboard barometric altimeter. The position solution of the follower UAV at time t k + 1 can be solved by combining the moving vector of the follower UAV from time t k to t k + 1 with the circular equation.
( x k + 1 D x k , k + 1 x k m ) 2 + ( y k + 1 D y k , k + 1 y k m ) 2 = r k 2 ( x k + 1 x k + 1 m ) 2 + ( y k + 1 y k + 1 m ) 2 = r k + 1 2
The first-order linearization Taylor expansion of Equation (8) is made and ignores the higher-order term above the second order, and by combining the altitude information output from the barometric altimeter, the position observation equation of the follower UAV’s INS at time t k + 1 can be obtained. The position observation equation of the follower UAV’s INS at time t k + 1 is constructed as follows:
                                                                        Z k + 1 p = H k + 1 p X k + 1 + v p                                                                     H k + 1 p = [ 0 3 × 6 H 3 × 3 ρ 0 3 × 9 ] 3 × 18 H 3 × 3 ρ = [ ( R N + h ) sin L ( cos λ f 11 + sin λ f 12 ) ( R N + h ) cos L ( sin λ f 11 cos λ f 12 ) cos L ( cos λ f 11 + sin λ f 12 ) ( R N + h ) sin L ( cos λ f 21 + sin λ f 22 ) ( R N + h ) cos L ( sin λ f 21 cos λ f 22 ) cos L ( cos λ f 21 + sin λ f 22 ) 0 0 1 ]
where Z k + 1 p is the position observation vector. H k + 1 p is the observation Jacobian matrix corresponding to the position measurement equation. X k + 1 is the state vector of the minimum cooperative unit’s positioning system. v p is the relative distance measurement white noise. ( L , λ , h ) represents the position information of the follower UAV in the navigation coordinate system. R N represents the curvature radius of the earth in prime vertical. f i 1 = f i x , f i 2 = f i y , i = 1 , 2 is the direction cosine of the position observation equation in the horizontal direction at time t k and t k + 1 , respectively.
The follower UAV’s onboard speed sensor can acquire the UAV’s forward velocity. The horizontal velocity can be obtained by decomposing the follower UAV’s forward velocity into the horizontal direction. The follower UAV’s horizontal velocity can be expressed as:
v x = v ^ cos φ ^ cos θ ^ v y = v ^ sin φ ^ cos θ ^
In the above formula, v x is the follower UAV’s velocity in x direction. v y is the follower UAV’s velocity in Y direction. Based on Formula (10), the horizontal velocity observation equation of the follower UAV’s INS can be written as follows:
Z k + 1 v = H k + 1 v X k + 1 + v v H k + 1 v = 0 2 × 3 1 0 0 0 1 0 0 2 × 12 2 × 18
where Z k + 1 v is the velocity observation vector. H k + 1 v is the Jacobian matrix corresponding to the velocity measurement equation. V v is the velocity measurement white noise.
At the same time, the follower UAV’s onboard three-axis magnetic sensor can provide heading information for the UAV through a sensitive geomagnetic field, which can be used as the heading observation to correct the error of the follower UAV’s INS. After temperature and magnetic compensation, the heading information output from three-axis magnetic sensor can be expressed as:
φ ^ m = φ ¯ m + v φ
where φ ^ m represents the measured value of the heading angle of the magnetic sensor. φ ¯ m represents the true value of heading angle measurement of the magnetic sensor. v φ is the magnetic sensor measurement white noise. By combining the heading information calculated by the follower UAV’s INS and the spatial coordinate system’s transformation relationship between the attitude angle and platform misalignment angle, the heading observation equation of the follower UAV’s INS can be constructed by Formula (11) as:
Z k + 1 φ = H k + 1 φ X k + 1 + v φ H k + 1 φ = H φ 0 1 × 15 1 × 18 H φ = 1 cos θ sin φ sin θ cos φ sin θ cos θ 1 × 3
where Z k + 1 φ is the heading observation vector. ( θ , φ ) represents the pitch angle and heading angle calculated by the follower UAV’s INS, respectively. H k + 1 φ is the Jacobian matrix corresponding to the heading observation equation. v φ is the magnetic sensor measurement white noise. Assuming that the relative distance measurement white noise, velocity measurement white noise and heading measurement white noise are independent of each other and satisfy the mean value being zero, the covariance is a fixed value.
E { v p ( k ) v v ( k ) } = 0 E { v φ ( k ) v v ( k ) } = 0 E { v p ( k ) v φ ( k ) } = 0 E { v φ ( k ) T v φ ( k ) } = R φ ( k ) E { v v ( k ) T v v ( k ) } = R v ( k ) E { v p ( k ) T v p ( k ) } = R p ( k )

4. Cooperative Positioning Algorithm Design Based on Improved Extended Kalman Filtering

According to the follower UAV’s INS state recursive equation, more details about the follower UAV’s INS state recursive equation can be found in [23]. The state propagation equation of the improved extended Kalman filtering is constructed as follow:
X ( k + 1 ) = f ( X ( k ) , k ) + w ( k )
where
X ( k ) = [ φ E ( k ) , φ N ( k ) , φ U ( k ) , δ v E ( k ) , δ v N ( k ) , δ v U ( k ) , δ L ( k ) , δ λ ( k ) , δ h ( k ) , ε b x ( k ) , ε b y ( k ) , ε b z ( k ) , ε r x ( k ) , ε r y ( k ) , ε r z ( k ) , a x ( k ) , a y ( k ) , a z ( k ) ] T
The system state vector X ( k ) can be, respectively, expressed as: platform misalignment angle error of the follower UAV’s INS, velocity error of the follower UAV’s INS, position error of the follower UAV’s INS, gyroscope random constant drift, gyroscope first-order Markov white noise and accelerometer first-order Markov white noise. w ( k ) is the system noise input matrix. Assuming that w ( k ) is the Gaussian white noise with a zero mean and satisfies E { w ( k ) w T ( k ) } = Q ( k ) , the state transition Jacobian matrix F ( k + 1 , k ) and system noise control input Jacobian matrix can be expressed as:
F ( k + 1 , k ) = f ( X ( k ) , k ) X | X = X ^ ( k | k ) W ( k + 1 , k ) = f ( X ( k ) , k ) u | X = X ^ ( k | k )
The airborne datalink system based on radio data transmitting and receiving is susceptible to radio signal interference. It is difficult for the follower UAV to receive communication data packets containing relative distance and the leader UAV’s position information at each sampling period. The position observation equation of the follower UAV’s INS has the following form:
Z k + 1 P = H k + 1 p X k + 1 + v p   if   communication   packet   received   successfully   0                   otherwise  
If the communication data packet from the leader UAV is successfully received at time t k , the position observation equation of the follower UAV’s INS can be expressed as Z k + 1 p = H k + 1 p X k + 1 + v p . Otherwise, the position observation equation of the follower UAV’s INS is zero. In this paper, in order to improve the follower UAV’s positioning accuracy, the random variable γ k { 0 , 1 } is introduced to describe whether the communication data packet arrived at time t k . γ k = 1 means that the communication data packet at time t k has been successfully received. In contrast, if γ k = 0 means that the communication data packet at time t k is not successfully received, assuming that the γ k follows the two-state Markov chain, the state transition matrix is
m = 1 q q p 1 p
where p , q > 0 , respectively, indicates the transmission failure and success rate of the communication data packet. The smaller the value of p or the larger the value of q , the better transmission performance of the airborne datalink system. Based on the above definition, the probability distribution of the follower UAV’s position measurement noise { v k , k 1 } is
Ρ ( v p | γ k ) ~ Ν ( 0 , R p ) , γ k = 1 Ν ( 0 , σ 2 I ) ,   γ k = 0
where σ 2 is a large positive number. The above formula indicated that if the follower UAV successfully receives the communication data packet at time t k , the measurement variance is R p . In contrast, the measurement variance is σ 2 I , σ . Due to the velocity and heading observation of the follower UAV’s INS being obtained from onboard sensors, there is no need for data transmission through the airborne datalink system. Therefore, the velocity and heading observation variance of the follower UAV’s INS at time t k are not affected by the datalink system. The position observation of the follower UAV’s INS before time t k is defined as follows:
Z k p = ( Z 1 p , Z 2 p , Z 3 p Z k p ) T γ k = ( γ 1 , γ 2 , γ 3 , γ k ) T
Under the condition of the position observation of the follower UAV’s INS being available, the covariance estimation equation of the traditional extended Kalman filtering can be expressed as:
X ^ k E { X k | Z k p , γ k }
P k E { ( X k X ^ k ) ( X k X ^ k ) T | Z k p , γ k }
P k + 1 / k E { ( X k + 1 X ^ k + 1 / k ) ( X k + 1 X ^ k + 1 / k ) T | Z k p , γ k }
Z ^ k + 1 / k p E { Z k + 1 p | Z k p , γ k }
However, under the condition of the datalink system being indirectly limited, the following equations are derived from the above formulas:
E { ( Z k + 1 p Z ^ k + 1 / k p ) ( X k + 1 p X ^ k + 1 / k ) T | Z k p , γ k + 1 } = H k + 1 p P k + 1 / k
E { ( Z k + 1 p Z ^ k + 1 / k p ) ( Z k + 1 p Z ^ k + 1 / k p ) T | Z k p , γ k + 1 } = H k + 1 p P k + 1 / k ( H k + 1 p ) T + γ k + 1 R k + 1 p + ( 1 γ k + 1 ) σ 2 I
Due to X k + 1 / k and Z k + 1 / k p being independent normal random vectors, its joint probability distribution under the conditions Z k p and γ k + 1 is still following normal distribution and satisfies:
E { X k + 1 / k , Z k + 1 / k p | Z k p , γ k + 1 } = X ^ k + 1 / k H k + 1 p X ^ k + 1 / k
cov { X k + 1 / k , Z k + 1 / k p | Z k p , γ k + 1 } = P k + 1 / k                                 P k + 1 / k ( H k + 1 p ) T   H k + 1 p P k + 1 / k H k + 1 p p k + 1 / k ( H k + 1 p ) T + γ k + 1 R k + 1 p + ( 1 γ k + 1 ) σ 2 I
Therefore, when the airborne datalink system is indirectly limited, the improved extended Kalman filtering equations are as follows:
X ^ k + 1 / k = F k X k
P k + 1 / k = F k P k F k T + Q k
X ^ k + 1 = X ^ k + 1 / k + γ k + 1 K k + 1 ( Z k + 1 p H k + 1 p X ^ k + 1 / k )
P k + 1 = P k + 1 / k γ k + 1 K k + 1 H k + 1 p P k + 1 / k
K k + 1 = P k + 1 / k ( H k + 1 p ) T ( H k + 1 p P k + 1 / k ( H k + 1 p ) T + R k + 1 ) 1

5. Multi-UAVs Flight Test and Results

An outdoor autonomous flight CP experiment involving two UAVs has been carried out on the campus playground. All sensors mounted on the two UAVs are shown in Figure 5. The performance parameters of airborne sensors are shown in Table 1. In the experiment, the leader UAV acquired high-precision position information by an INS/GPS loosely integrated navigation system. The forward velocity output from the follower UAV’s single-point GPS receiver was used to simulate the output of the speed sensor. The ultra-wideband (UWB) sensor was used to realize relative distance measurement in real time. In order to evaluate the follower UAV’s positioning performance, the real-time kinematic (RTK) navigation system mounted on the follower UAV was used as the reference system. The position information output from the RTK is only used as the position reference and does not participate in CP calculation. The three-axis magnetic sensor module with temperature and magnetic compensation and the barometric altimeter module mounted on the follower UAV, respectively, measure the UAV’s heading and altitude information. Before the experiment, the sampling frequency of the follower and leader UAV’s INS were set to 100Hz. The sampling frequencies of the other sensors were all set to 1Hz. In order to compare the positioning performance of the proposed algorithm and the traditional multi-leader UAV CP algorithm based on relative ranging, we made full use of five UWB base stations on the playground to replace the known position leader’s UAV nodes. The positions of the UWB base station were calibrated by GPS in a single-point mode before the experiment. During the experiment, the UWB tag on the follower UAV could obtain relative distance information relative to the UWB base station. Due to the influence of UAV’s motion, the UWB tag mounted on the follower UAV could not obtain relative distance information from the leader UAV at each sampling period. A two-state Markov chain was introduced to judge the availability of the relative distance. When the relative distance information with the leader UAV is available, the position observation equation of the follower UAV’s INS is constructed using the follower UAV’s moving vector. During the experiment, all sensor data were captured from the follower UAV’s processor. The feasibility of the proposed algorithm was verified by post-processing. The total flight time was about 109 s.
The autonomous flight CP experiment for the two UAVs is shown in Figure 6. The flight path of the two UAVs obtained from the leader and follower UAV navigation processors is shown in Figure 7. Figure 8 illustrates the positioning performance of the follower UAV generated by the traditional multi-leader UAV CP method based on relative ranging and the proposed method.
It can be seen from Figure 8 that during the 109 s flight time, the follower UAV’s positioning error was basically kept within 2 m by using the proposed method. Under the condition that the measurement accuracy of the sensors used in the experiment is consistent, the method proposed in this paper has higher positioning accuracy compared to the traditional multi-leader UAV CP method based on relative ranging.
In order to more visually compare the follower UAV’s positioning performance by using the proposed method and the traditional multi-leader UAV CP method based on relative ranging, the root mean square error (RMSE) value of the follower UAV’s positioning error is compared in Table 2. It can be seen from Table 2 that, under the condition that the follower UAV acquires relative distance information from the leader UAV at discontinuous times, the CP algorithm proposed in this paper has higher positioning accuracy than the traditional method. The follower UAV’s positioning accuracy is improved at least 50% compared to the traditional method. For the purpose of quantitatively comparing and analyzing the positioning performance of the proposed algorithm from various perspectives, Figure 9 gives the mean and variance of the follower UAV’s positioning error by using the traditional CP method and the proposed method. We can conclude from the Figure 9 that when there are not enough known position leader UAVs in the UAV cooperative network, the proposed algorithm makes full use of the historical cooperative information to construct position observations. The position closed-form solution of the follower UAV can be acquired in real time. Better results by using the proposed method can be obtained and positioning accuracy can meet the requirement of UAV formation cooperative flight. The results of the analysis based on mean and variance are consistent with the results of the RMSE.
The method proposed in this paper breaks through the limitation of the traditional CP algorithm on the number of known position leader UAV nodes, which lays a theoretical foundation for the future engineering application.

6. Conclusions

The traditional multi-leader UAV CP algorithm based on relative ranging has a requirement on the number of known position leader UAV nodes. When the number of leader UAVs and relative distance measurements are limited, the traditional CP algorithm is not applicable. Aimed at the minimum cooperative unit composed of a known position leader UAV and an unknown position follower UAV, this paper proposes a cooperative localization method based on the follower UAV’s moving vector. Compared to the traditional multi-leader UAV CP method based on relative ranging, only the single relative distance information with the leader UAV is required to obtain the follower UAV’s position solution. By combining the information of the follower UAV’s airborne speed sensor and three-axis magnetic sensor, the moving vector of the follower UAV at adjacent times is constructed. The position observation equation of the follower UAV’s INS is constructed by using a relative distance constraint and the follower UAV’s moving vector. Meanwhile, the velocity and heading observation equations of the follower UAV’s INS are constructed by using the information of the speed and magnetic sensor. The improved extended Kalman filtering is designed to estimate and compensate the system state vector to improve the follower UAV’s positioning accuracy. In addition, considering that the datalink system based on radio signals may be interfered with by the external environment, this paper introduces a two-state Markov chain to judge the availability of the relative distance information. The real autonomous flight CP experiment is conducted to validate the performance of the proposed algorithm. The results show that the proposed CP algorithm can achieve better positioning accuracy compared to the traditional multi-leader UAV CP algorithm.

Author Contributions

Conceptualization, S.C.; methodology, X.Z.; software, X.Z.; validation, X.Z.; formal analysis, J.L.; investigation, S.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 61973160.

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. Sun, J.; Li, B.; Jiang, Y.; Wen, C.-Y. A camera-based target detection and positioning UAV system for search and rescue (SAR) purposes. Sensors 2016, 16, 1778. [Google Scholar] [CrossRef] [PubMed]
  2. Liu, Y.; Liu, Z.; Shi, J.; Wu, G.; Chen, C. Optimization of base location and patrol routes for unmanned aerial vehicles in border intelligence, surveillance, and reconnaissance. J. Adv. Transp. 2019, 2019, 9063232. [Google Scholar] [CrossRef]
  3. Lázaro, M.T.; Castellanos, J.A. Localization of probabilistic robot formations in SLAM. In Proceedings of the International Robotics and Automation Conference, Anchorage, Alaska, 3–8 May 2010; pp. 3179–3184. [Google Scholar]
  4. Shaferman, V.; Shima, T. Unmanned aerial vehicles cooperative tracking of moving ground target in urban environments. J Guid Control Dyn. 2008, 31, 1360–1371. [Google Scholar] [CrossRef]
  5. Raju, V.A.; Vasundhara, P.; Reddy, V.C.K.; Aiswarya, A.S. Improvement of position and orientation of Unmanned Arial Vehicle (UAV) with INS/GPS. Int. J. Eng. Technol. 2018, 7, 642–646. [Google Scholar] [CrossRef]
  6. Guangcai, W.; Xu, X.; Zhang, T. M-M Estimation-Based Robust Cubature Kalman Filter for INS/GPS Integrated Navigation System. IEEE Trans. Instrum. Meas. 2020, 70, 1–11. [Google Scholar] [CrossRef]
  7. Outamazirt, F.; Fu, L.; Lin, Y.; Abdelkrim, N. A new SINS/GPS sensor fusion scheme for UAV localization problem using nonlinear SVSF with covariance derivation and an adaptive boundary layer. Chin. J. Aeronaut. 2016, 29, 424–440. [Google Scholar] [CrossRef]
  8. Khalaf, W.; Chouaib, I.; Wainakh, M. Novel adaptive UKF for tightly coupled INS/GPS integration with experimental validation on an UAV. Gyroscopy Navig. 2017, 8, 259–269. [Google Scholar] [CrossRef]
  9. Huang, S.; Huang, J.; Tang, D.; Chen, F. Research on UAV flight performance test method based on dual antenna GPS/ins integrated system. In Proceedings of the 3rd International Communication and Information Systems (ICCIS) Conference, Singapore, 28–30 September 2018; pp. 106–116. [Google Scholar]
  10. Cao, S.; Qin, H.; Cong, L.; Huang, Y. TDMA Datalink Cooperative Navigation Algorithm Based on INS/JTIDS/BA. Electronics 2021, 10, 782. [Google Scholar] [CrossRef]
  11. Duan, Q. A High Precision Time Synchronization Scheme for Avionics System. Telecommun. Eng. 2014, 54, 830–834. [Google Scholar]
  12. Ranger, J.F.O. Principles of JTIDS relative navigation. J. Navig. 1996, 49, 22–35. [Google Scholar] [CrossRef]
  13. Yang, C.; Strader, J.; Gu, Y.; Hypes, A.; Canciani, A.; Brink, K. Cooperative UAV navigation using Inter-Vehicle Ranging and Magnetic Anomaly Measurements. In Proceedings of the 2018 AIAA Guidance, Navigation and Control Conference, Kissimmee, FL, USA, 8–12 January 2018; p. 1595. [Google Scholar]
  14. Qu, Y.; Zhang, Y. Cooperative localization of low-cost UAV using relative range measurements in multi-UAV flight. In Proceedings of the AIAA Guidance, Navigation and Control Conference, Toronto, ON, Canada, 2–5 August 2010; p. 8187. [Google Scholar]
  15. 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. Rob. Syst. 2020, 17, 1729881420932717. [Google Scholar] [CrossRef]
  16. Jun, C.; Zheng, C.; Sun, D.; Zhang, D. AUV positioning based on single-beacon ranging in straight-line trajectory. J. Harbin Eng. Univ. 2017, 38, 733–738. [Google Scholar]
  17. Fangfang, C.; Weidong, L.; Juanli, L. Navigating and positioning based on EKF for unmanned underwater vehicle from a single beacon. Comput. Meas. Control 2011, 19, 2234–2237. [Google Scholar]
  18. Zhang, F.B.; Zhang, Y.Q. Correcting localization error with a single beacon for AUV. Torpedo Technol. 2012, 20, 38–41. [Google Scholar]
  19. Kepper, J.H., IV; Claus, B.C.; Kinsey, J.C. A Navigation Solution Using a MEMS IMU, Model-Based Dead-Reckoning, and One-Way-Travel-Time Acoustic Range Measurements for Autonomous Underwater Vehicles. IEEE J. Ocean. Eng. 2018, 44, 664–682. [Google Scholar] [CrossRef]
  20. Deng, Z.C.; Yu, X.; Qin, H.D.; Zhu, Z.B. Adaptive Kalman Filter-Based Single-Beacon Underwater Tracking with Unknown Effective Sound Velocity. Sensors 2018, 18, 4339. [Google Scholar] [CrossRef] [PubMed]
  21. Sun, J.; Hu, F.; Jin, W.; Wang, J.; Wang, X.; Luo, Y.; Yu, J.; Zhang, A. Model-aided localization and navigation for underwater gliders using single-beacon travel-time differences. Sensors 2020, 20, 893. [Google Scholar] [CrossRef] [PubMed]
  22. Qin, H.D.; Yu, X.; Zhu, Z.B.; Deng, Z.C. An expectation-maximization based single-beacon underwater navigation method with unknown ESV. Neurocomputing 2020, 378, 295–303. [Google Scholar] [CrossRef]
  23. Xu, J.; Xiong, Z.; Liu, J.; Wang, R. A dynamic vector-formed information-sharing algorithm based on two-state chi square detection in an adaptive federated filter. J. Navig. 2019, 72, 101–120. [Google Scholar] [CrossRef]
Figure 1. A known position leader UAV and an unknown position follower UAV cooperative localization scenario based on relative ranging.
Figure 1. A known position leader UAV and an unknown position follower UAV cooperative localization scenario based on relative ranging.
Sensors 22 07125 g001
Figure 2. The schematic diagram of proposed method in this paper.
Figure 2. The schematic diagram of proposed method in this paper.
Sensors 22 07125 g002
Figure 3. Schematic diagram of follower UAV’s moving vector in the horizontal direction.
Figure 3. Schematic diagram of follower UAV’s moving vector in the horizontal direction.
Sensors 22 07125 g003
Figure 4. The spatial position relationship in the horizontal direction between leader and follower UAV at adjacent times.
Figure 4. The spatial position relationship in the horizontal direction between leader and follower UAV at adjacent times.
Sensors 22 07125 g004
Figure 5. Sensor mounted on the UAVs in autonomous flight experiment.
Figure 5. Sensor mounted on the UAVs in autonomous flight experiment.
Sensors 22 07125 g005
Figure 6. Two UAVs autonomous flight cooperative localization experiment.
Figure 6. Two UAVs autonomous flight cooperative localization experiment.
Sensors 22 07125 g006
Figure 7. The flight path of leader and follower UAV.
Figure 7. The flight path of leader and follower UAV.
Sensors 22 07125 g007
Figure 8. Follower UAV’s positioning error curve comparison by using traditional multi-leader UAV CP method based on relative ranging and the proposed method.
Figure 8. Follower UAV’s positioning error curve comparison by using traditional multi-leader UAV CP method based on relative ranging and the proposed method.
Sensors 22 07125 g008
Figure 9. Mean and variance of follower UAV’s positioning error using traditional method and proposed method.
Figure 9. Mean and variance of follower UAV’s positioning error using traditional method and proposed method.
Sensors 22 07125 g009
Table 1. The performance parameters of airborne sensors in the experiment.
Table 1. The performance parameters of airborne sensors in the experiment.
Airborne Sensor TypesPerformance Parameters
Three-axis gyroscopeBias stabilityRandom work error
20 (°/h)0.2 ( ° / h )
Three-axis accelerometerBias stabilityRandom work error
45 µg 10 µg Hz
Three-axis magnetic sensor0.5°
Barometric altimeter0.5 m
UWB0.3 m
Speed sensor0.2 m/s
Table 2. Comparison of the RMSE of the follower UAV’s positioning error by using traditional CP method and the proposed method.
Table 2. Comparison of the RMSE of the follower UAV’s positioning error by using traditional CP method and the proposed method.
MethodLongitude (m)Latitude (m)Altitude (m)
Traditional method1.5481.2740.808
Proposed method0.5550.5050.164
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhu, X.; Lai, J.; Chen, S. Cooperative Location Method for Leader-Follower UAV Formation Based on Follower UAV’s Moving Vector. Sensors 2022, 22, 7125. https://doi.org/10.3390/s22197125

AMA Style

Zhu X, Lai J, Chen S. Cooperative Location Method for Leader-Follower UAV Formation Based on Follower UAV’s Moving Vector. Sensors. 2022; 22(19):7125. https://doi.org/10.3390/s22197125

Chicago/Turabian Style

Zhu, Xudong, Jizhou Lai, and Sheng Chen. 2022. "Cooperative Location Method for Leader-Follower UAV Formation Based on Follower UAV’s Moving Vector" Sensors 22, no. 19: 7125. https://doi.org/10.3390/s22197125

APA Style

Zhu, X., Lai, J., & Chen, S. (2022). Cooperative Location Method for Leader-Follower UAV Formation Based on Follower UAV’s Moving Vector. Sensors, 22(19), 7125. https://doi.org/10.3390/s22197125

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