Next Article in Journal
A Penetration Method for UAV Based on Distributed Reinforcement Learning and Demonstrations
Next Article in Special Issue
Multi-UAV Collaborative Absolute Vision Positioning and Navigation: A Survey and Discussion
Previous Article in Journal
High Precision Mesh-Based Drone Image Stitching Based on Salient Structure Preservation and Regular Boundaries
Previous Article in Special Issue
Multi-UAV Cooperative Trajectory Planning Based on FDS-ADEA in Complex Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A UAV Formation Control Method Based on Sliding-Mode Control under Communication Constraints

1
Coast Guard Academy, Naval Aviation University, Yantai 264000, China
2
Weapons Academy, Naval Engineering University, Wuhan 430000, China
*
Author to whom correspondence should be addressed.
Drones 2023, 7(4), 231; https://doi.org/10.3390/drones7040231
Submission received: 28 February 2023 / Revised: 13 March 2023 / Accepted: 22 March 2023 / Published: 27 March 2023
(This article belongs to the Special Issue A UAV Platform for Flight Dynamics and Control System)

Abstract

:
The problem of vision-based fixed-wing UAV formation control under communication limitations and the presence of measurement errors was investigated. In the first part of this paper, the single UAV motion model and the process of estimating the neighboring UAV states using the Extended Kalman Filter are introduced. The second part describes how we designed a sliding mode controller considering the sensor measurement errors and discusses the sufficient conditions for the stability and formation system in the presence of state transfer time delays in the formation. The main motivation of this paper was to develop a hierarchical, globally stable sliding mode controller that could enable the considered vision-based multiple fixed-wing UAVs to achieve and maintain formation in the presence of measurement errors. To this end, the selected problem was first transformed into a state-tracking problem for UAVs in the neighborhood, and then the stability criterion was established using the Lyapunov stability theory. Finally, the effectiveness of the proposed control method was illustrated using three numerical arithmetic examples.

1. Introduction

In recent years, with the development of technology, UAV formations are able to integrate more sensors compared to single UAVs to achieve diverse tasks in complex mission environments, so UAV formations are widely used in activities such as search [1], reconnaissance and target tracking [2], aerial refueling [3], and cloud computing [4]. In ref. [5], Alsamhi et al. provide a thorough review of the application of UAVs in search. Additionally, considering that although rotary-wing UAVs have the advantages of being mobile and flexible and able to hover at will, but their disadvantages, such as limited load capacity and short flight duration, cannot be ignored, this paper selects fixed-wing UAV formations capable of performing various complex tasks with long flight times as the object of study.
In previous studies, UAV formation control with communication links has become a research topic for researchers because of stable inter-aircraft information transmission, easy implementation, and the ability to adapt to most mission situations to achieve a variety of formation control strategies. According to the control strategy, it can be classified as Leader-Follower [6,7,8], Virtual Structure [9,10], Behavior-Based [11,12,13], and Consensus-Based [14,15]. A distributed control method was proposed in [16] to effectively reduce the requirement of the follower UAV speed and achieve formation control under the follower speed no greater than the leader speed limit. Additionally, in [17], to reduce the formation error due to parameter uncertainty, the feedback linearization approach was used to design the formation controller to achieve a zero dynamic error in formation. Unlike the previous two, Zhao et al. in [18] designed controllers with full consideration of UAV input constraints as well as inter-aircraft collision avoidance constraints and implemented simulation experiments in a 3D environment. Kim et al. [19] designed control laws for the dynamic compensation of UAVs during turns and implemented algorithmic simulations of hardware-in-the-loop. Arumugam et al. in [20] derived sufficient conditions for nonlinear intelligence consistency based on matrix linear inequalities and later in [21] studied multi-intelligence consistency for semi-Markov jump in the perturbed case.
However, in recent years, with the complexity of the mission environment and various types of communication sabotage facilities for UAV formations, it was difficult for UAV formations to form effective communication links within the mission area. Therefore, formation control methods that rely only on visual information transfer are gradually being emphasized by researchers and are widely used in various mission environments. Refs. [22,23,24] use adaptive primitives to achieve control enhancement without modeling the dynamics and interference, while ref. [25] builds on the former to design the adaptive control of the guiding law line of sight to achieve formation control on a pair of UAVs. A control law based on the inverse height between the optical center of the camera and the feature point was proposed in [26] to achieve formation control under visibility constraints. The ref. [27], on the other hand, designs adaptive control laws for formation control in the absence of a distance state using position navigator velocity and vision-related parameters. To reduce the constraint of information transfer direction, a pilotless formation control method was proposed in [28] and achieved obstacle avoidance and reorganization of the formation. A controller achieved large-scale cluster formation in a measurable error environment, which was presented by Chen et al. in [29] and realized comparative simulations in multiple cases.
Although the above-mentioned article implements vision-based formation control, it makes high demands on the controller design due to the measurement error of the sensor in estimating the state of the leader UAV and the need to form a stable form in a short period. The sliding mode control, on the other hand, is not sensitive to error and has a good control effect, and is widely used by researchers to realize the formation control of UAV formation. To solve the drawbacks of jitter in the sliding mode control method and the difficulty in setting the control gain, adapting the sliding mode control was applied to formation control [30]. The refs. [31,32] integrated sliding-mode controllers into feedback linearization to achieve formation control in the case of unknown leader UAV speed. A bilateral control law using linear quadratic adjustment of the outer loop control and sliding mode control achieved inner loop attitude control and was proposed by Ghamry et al. in [33] to solve the problems of the long computation time of the sliding mode controller and system jitter. To achieve the formation control of multiple aircraft, ref. [34] used a higher-order sliding mode control to achieve the extension of UAV formation. The UAV formation control method is applicable to the model-free case and was proposed in the literature [35], while Mobayen et al. designed two adaptive laws [36] to deal with uncertainties and disturbances with unknown upper bounds, and demonstrated the effectiveness of the algorithm in practical applications. The work in [37], on the other hand, used an integral sliding mode control to achieve the formation control of three UAVs with little change in the speed of the leader.
The above research results show that the sliding mode control algorithm has the advantages of being insensitive to state estimation errors and fast tracking, but there are problems with large initial input values and system jitter. To solve the above problems, this paper adopts the two-layer global stable sliding mode control to design the controller. While ensuring the fast convergence of the system, the sliding mode control input amplitude was limited to ensure that the UAV formation could be controlled under the input limitation and the existence of state estimation delay. The contributions of this paper can be summarized as follows:
(1)
The formation control problem of fixed-wing UAVs in an uncommunicated situation was studied. UAVs use only onboard vision sensors to obtain the status of UAVs in the neighborhood and to achieve formation control.
(2)
In this paper, we considered sensor measurements and used Extended Kalman Filtering to reduce errors and designed a sliding mode controller to reduce sensitivity to measurement errors.
(3)
In this paper, the formation control was divided into a two-layer control structure of position control and attitude control, and a suitable Lyapunov function was constructed to prove the system’s stability.
(4)
The UAV state transfer delay was considered, the UAV error state equation was given, and a suitable Lyapunov- Krasovskii generalized function was constructed to derive sufficient conditions for the stability of the delayed system, and finally, the obtained theoretical results were illustrated by three numerical simulations.
The rest of the article is outlined below. Section 2 describes the process of estimating the neighborhood UAV state using the Extended Kalman Filter. Section 3 describes the controller development process and sufficient conditions for system stability in the presence of time-delayed states. Section 4 provides numerical simulations to determine the effectiveness of the proposed control method. Finally, then Section 5 gives conclusions.

2. Extended Kalman Filter-Based State Estimation

2.1. Equation of Motion for Single Machine

In order to introduce the UAV formation movement and modeling process, we introduced three coordinate frames I, Vi, and Li as the inertial coordinate frame, the follower velocity frame of the UAV i , and the line of sight frame. The origin of the follower velocity frame and the line of sight frame is the mass point of the UAV. The relative relationship between the coordinate frames is shown in Figure 1.
The inter–aircraft relative relationship diagram (shown in Figure 2) can be obtained from Figure 1, where l i j , χ i j , and γ i j are the relative distance, relative yaw angle, and relative pitch angle between the follower and the leader UAV in its neighborhood.
From [29], the equation of the state of motion of a single UAV is given by:
x ˙ i = [ x ˙ i y ˙ i z ˙ i v ˙ i χ ˙ i γ ˙ i ] = [ v i cos γ i cos χ i v i cos γ i sin χ i v i sin γ i a x i a z i v i - a y i v i ]
where: a i = [ a x i V i a y i V i a z i V i ] T 3 is the UAV control input. Considering the actual flight of the UAV, the input is not able to be chosen arbitrarily, considering the actual state of the UAV, and the input is restricted to a certain action area. Therefore, in this paper, we assumed that the UAV input satisfied:
a i = [ | a x i | | a y i | | a z i | ] [ a x max a y max a z max ]
[ a x max a y max a z max ] T in Equation (2) are known constants. The position error between the leader UAV and the follower UAV can be expressed as:
e x = [ e x e y e z ] = [ x i x l y i y l z i z l ]

2.2. Leader UAV Status Estimation

In this paper, considering the measurement error in the binocular camera, this paper used the Extended Kalman Filter to estimate the state of the leader UAV. The Extended Kalman Filter for the continuous measurement of the UAV state equation can be expressed as:
Ψ l ( k + 1 ) = A Ψ l ( k ) + w ( k )
z ( k ) = H ( Ψ l ( k ) Ψ i ( k ) ) + v ( k )
where, Ψ i ( k ) = [ x i y i z i v x i v y i v z i ] is the position and velocity state vector of Leader UAV. z ( k ) = [ l i j χ i j γ i j ] T is the relative information vector between the navigator and the follower where:
A = [ 1 0 0 T 0 0 0 1 0 0 T 0 0 0 1 0 0 T 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 ] ,   H = [ l i j x i l i j y i l i j z i l i j v x i l i j v y i l i j v z i χ i j x i χ i j y i χ i j z i χ i j v x i χ i j v y i χ i j v z i γ i j x i γ i j y i γ i j z i γ i j v x i γ i j v y i γ i j v z i ]
The estimation process of the Extended Kalman Filter can be divided into the following estimation parts as well as the gain update part,
Ψ ^ l ( k + 1 | k ) = A Ψ ^ l ( k | k ) P ( k + 1 | k ) = A P ( k | k ) A T + Q ( k )
where Ψ ^ l is the leader state estimation vector, P is the covariance matrix of the state errors, and is the covariance matrix of noise due to sensor errors in the state variables. After obtaining the state estimation vector, the state estimation vector was corrected by the following equation, and the Kalman gain was updated by Equation (8).
Ψ ^ l ( k + 1 | k + 1 ) = Ψ ^ l ( k + 1 | k ) + K ( k ) [ z ( k ) H ( Ψ l ( k + 1 | k ) Ψ i ( k ) ) ] P ( k + 1 | k ) = P ( k + 1 | k ) K ( k ) B P ( k + 1 | k )
K ( k ) = P ( k + 1 | k ) B T ( B P ( k + 1 | k ) B T + R ( k ) ) 1
where R ( k ) is the measurement noise caused by the measurement and B is the output matrix, which can be expressed in the following form.
B = I 6 × 6
By applying the extended Kalman filter, the leader state vector Ψ l ( k ) = [ x l y l z l v x l v y l v z l ] could be obtained, and the leader state could be obtained by calculating:
p l = [ x y z ] , v l = [ x ˙ y ˙ z ˙ ] = [ v x l v y l v z l ]
Through the above state estimation for UAVs in the neighborhood, the follower UAV can obtain the reference state information of UAVs in the neighborhood.

3. Controller Design

The controller designed in this paper controlled the inner and outer loops of the UAV separately to realize the position control of the outer loop and the attitude control of the inner loop, respectively. The ideal state of the inner-loop control was generated while realizing the outer-loop control, which then passed to the inner-loop, and then the controller design of the inner-loop realized the tracking of the passed signal and finally realized the closed-loop control of the whole system. In this paper, U A V l was used to represent the reference UAV in the neighborhood of U A V i . The flow chart of the algorithm is shown in Figure 3.
One note here: in the controller design of this paper, we will use the following second-order differential controller to implement the derivation of the intermediate ideal state differential generated by the outer loop.
{ x ˙ 1 = x 2 x ˙ 2 = 2 R 2 ( x 1 n ( t ) ) R x x y = x x
where n ( t ) is the input signal and is the component to be differentiated, x 1 achieves the tracking of the signal, and x 2 is an estimate of the first-order derivative of the signal n ( t ) .The initial states of x 1 and x 2 are defined as x 1 ( 0 ) = 0 , x 2 ( 0 ) = 0 .

3.1. Controller Design

For the UAV to have a smooth trajectory, this also meant that it needed to be a smooth curve to make the input to the UAV continuous in the 3D environment. From Equation (1), we obtained:
v i = x ˙ i 2 ( t ) + y ˙ i 2 ( t ) + z ˙ i 2 ( t )
χ i = arctan 1 ( y ˙ ( t ) x ˙ ( t ) )
γ i = arctan 1 ( z ˙ ( t ) x ˙ 2 ( t ) + y ˙ 2 ( t ) )
Together with Equations (1), (12)–(14), we can obtain the acceleration of the UAV as:
a x i = v ˙ i = x ˙ i ( t ) x ¨ i ( t ) + y ˙ i ( t ) y ¨ i ( t ) + z ˙ i ( t ) z ¨ i ( t ) x ˙ i 2 ( t ) + y ˙ i 2 ( t ) + z ˙ i 2 ( t )
a y i = v i γ ˙ i = z ¨ i ( t ) [ x ˙ i 2 ( t ) + y ˙ i 2 ( t ) ] z ˙ i ( t ) [ x ˙ i ( t ) x ¨ i ( t ) + y ˙ i ( t ) y ¨ i ( t ) ] x ˙ i 2 ( t ) + y ˙ i 2 ( t ) + z ˙ i 2 ( t ) x ˙ i 2 ( t ) + y ˙ i 2 ( t )
a z i = v i χ ˙ i = x ˙ i 2 ( t ) + y ˙ i 2 ( t ) + z ˙ i 2 ( t ) x ˙ i ( t ) y ¨ i ( t ) x ¨ i ( t ) y ˙ i ( t ) x ˙ i 2 ( t ) + y ˙ i 2 ( t )
From the above equation, the UAVs could achieve formation control in the case of multiple UAVs with smooth paths under continuous trajectories. It was assumed that the desired trajectory of the UAV in the 3D environment was:
[ x d ( t ) y d ( t ) z d ( t ) ] = [ f d ( t ) g d ( x d ( t ) ) h d ( x d ( t ) ) ]
Considering that the UAVs’ input is constrained, the following assumptions were made on the trajectory function.
| d g d ( x ) d x | λ 1 , | d 2 g d ( x ) d 2 x | λ 2 , x [ 0 , )
| d h d ( x ) d x | μ 1 , | d 2 h d ( x ) d 2 x | μ 2 , x [ 0 , )
where λ 1 , λ 2 , μ 1 , μ 2 is the determining constant. By substituting Equation (18) into Equations (15)–(17), which yielded:
a x i = f ˙ d 2 ( t ) ( d 2 g d ( x ) d 2 x + d 2 h d ( x ) d 2 x ) + f ¨ d ( t ) ( d g d ( x ) d x + d h d ( x ) d x ) 1 + ( d g d ( x ) d x ) 2 + ( d h d ( x ) d x ) 2
a y i = f ˙ d ( t ) d h d ( x ) d x ( 1 + ( d g d ( x ) d x ) 2 d g d ( x ) d x d 2 g d ( x ) d 2 x ) 1 + ( d g d ( x ) d x ) 2 + ( d h d ( x ) d x ) 2 1 + ( d g d ( x ) d x ) 2
a z i = 1 + ( d g d ( x ) d x ) 2 + ( d h d ( x ) d x ) 2 f ˙ d ( t ) d 2 g d ( x ) d 2 x 1 + ( d g d ( x ) d x ) 2
By applying Equations (21)–(23) to Equation (2), we obtained:
| f ˙ d 2 ( t ) ( d 2 g d ( x ) d 2 x + d 2 h d ( x ) d 2 x ) + f ¨ d ( t ) ( d g d ( x ) d x + d h d ( x ) d x ) 1 + ( d g d ( x ) d x ) 2 + ( d h d ( x ) d x ) 2 | a x max
| f ˙ d ( t ) d h d ( x ) d x ( 1 + ( d g d ( x ) d x ) 2 d g d ( x ) d x d 2 g d ( x ) d 2 x ) 1 + ( d g d ( x ) d x ) 2 + ( d h d ( x ) d x ) 2 1 + ( d g d ( x ) d x ) 2 | a y max
| 1 + ( d g d ( x ) d x ) 2 + ( d h d ( x ) d x ) 2 f ˙ d ( t ) d 2 g d ( x ) d 2 x 1 + ( d g d ( x ) d x ) 2 | a z max
Considering Equations (19) and (20), it was obtained that:
| f ˙ d 2 ( t ) ( λ 2 + μ 2 ) + f ¨ d ( t ) ( λ 1 + μ 1 ) 1 + λ 1 2 + μ 1 2 | a x max
| f ˙ d ( t ) μ 1 ( 1 + λ 1 2 λ 1 λ 2 ) 1 + λ 1 2 + μ 1 2 1 + λ 1 2 | a y max
| 1 + λ 1 2 + μ 1 2 λ 2 f ˙ d ( t ) 1 + λ 1 2 | a z max
From Equations (27)–(29), we obtained:
| f ˙ d ( t ) | λ 3 , λ 3 = min ( a x max , a y max 1 + λ 1 2 + μ 1 2 1 + λ 1 2 1 + λ 1 2 λ 1 λ 2 , a z max ( 1 + λ 1 2 ) λ 2 1 + λ 1 2 + μ 1 2 )
| f ¨ d ( t ) | a x max 1 + λ 1 2 + μ 1 2 λ 1 + μ 1
It can be derived from the above equation that when the system input was bounded, the boundedness of the trajectory could be guaranteed. The design of the outer loop controller is discussed below. From Equation (3), the system error tracking equation can be obtained as:
e ˙ i = [ x ˙ i x ˙ l y ˙ i y ˙ l z ˙ i z ˙ l ] = [ v i cos γ i cos χ i v x l v i cos γ i sin χ i v y l v i sin γ i v z l ]
The outer-loop virtual input can be taken as:
u i = [ u 1 u 2 u 3 ] = [ v i cos γ i cos χ i v i cos γ i sin χ i v i sin γ i ] = [ v x l a 1 tanh ( ρ 1 e x ) v y l a 2 tanh ( ρ 2 e y ) v z l a 3 tanh ( ρ 3 e z ) ]
To achieve the inner loop control when the input is u i , the ideal tracking state of the inner loop can be calculated according to Equation (33) as:
v d = u 3 sin γ d χ d = arctan ( u 2 u 1 ) γ d = arctan ( u 3 u 1 2 + u 1 2 )
Take the inner-loop sliding mode functions as:
s v = v i v d s χ = χ i χ d s γ = γ i γ d
The system input can be taken as:
a x i = v ˙ d k 1 s v τ 1 sign s v b 1 tanh ρ 4 s v
a y i = v d ( γ ˙ d k 2 s γ τ 2 sign s γ b 2 tanh ρ 5 s γ )
a z i = v d ( χ ˙ d k 3 s χ τ 3 sign s χ b 3 tanh ρ 6 s χ )
when the system input is (36)–(38), and the system can ensure global stability (Stability proofs are made in the next section).

3.2. System Stability Proof

With the input of Equation (33) for the outer loop system, the tracking error system (32) can become:
e ˙ i = [ x ˙ i x ˙ l y ˙ i y ˙ l z ˙ i z ˙ l ] = [ a 1 tanh ( ρ 1 e x ) a 2 tanh ( ρ 2 e y ) a 3 tanh ( ρ 3 e z ) ]
To prove the stability of the outer-loop system, Theorem 1 is given.
Theorem 1. 
When the system has a tame solution, as shown in equation:
η ˙ = α tanh ( k η )
for any constant α . η satisfying α . η > 0 , the system converges to 0 in any initial state η ( 0 ) .
The proof procedure is shown in [38] and is not repeated here. For the inner-loop system, the Lyapunov function of the inner-loop tracking error system can be taken as:
V 1 = 1 2 s v 2 + 1 2 s χ 2 + 1 2 s γ 2
Taking the derivative of Equation (41) and substituting it into Equations (1) and (35), we obtained:
V ˙ 1 = s v s ˙ v + s χ s ˙ χ + s γ s ˙ γ = s v ( v ˙ i v ˙ d ) + s χ ( χ ˙ i χ ˙ d ) + s γ ( γ ˙ i γ ˙ d ) = s v ( a x i v ˙ d ) + s χ ( a z i v d χ ˙ d ) + s γ ( a y i v d γ ˙ d )
By substituting into Equations (36)–(38), Equation (42) was simplified to obtain:
V ˙ 1 = s v ( a x i v ˙ d ) + s χ ( a z i v d χ ˙ d ) + s γ ( a y i v d γ ˙ d ) = ( k 1 s v 2 + k 2 s χ 2 + k 3 s γ 2 ) ( τ 1 s v sign s v + τ 2 s χ sign s χ + τ 3 s γ sign s γ ) ( b 1 s v tanh ρ 4 s v + b 2 s χ tanh ρ 5 s χ + b 3 s γ tanh ρ 6 s γ ) < ( b 1 s v tanh ρ 4 s v + b 2 s χ tanh ρ 5 s χ + b 3 s γ tanh ρ 6 s γ )
Based on Equation (43) and Theorem 1, it was obtained that the inner loop of the system was convergent under the action of the input, and when the parameter values were appropriate, the inner loop converged faster than the outer loop, and the inner and outer loops converged equally when and only when the error was zero. The global stability of the system is analyzed below. Assuming that there exists a pitch angle as well as a deflection angle to satisfy the follower for the leader, the position state part of Equation (1) can be written as:
{ x ˙ i = v d cos γ d cos χ d + cos γ d cos χ d ( v i v d ) + v i cos γ i ( cos χ i cos χ d ) + v i cos χ d ( cos γ i cos γ d ) y ˙ i = v d cos γ d sin χ d + cos γ d sin χ d ( v i v d ) + v i cos γ i ( sin χ i sin χ d ) + v i sin χ d ( cos γ i cos γ d ) z ˙ i = v d sin γ d sin γ d ( v i v d ) v i ( sin γ i sin γ d )
From Equation (44), the state tracking system produces corresponding deviations when there are errors in χ i , χ d as well as γ i , γ d . Considering that the system inputs (36) to (38) were bounded, the system state quantity v i , χ i , γ i satisfies the global Lipschitz condition. Substituting the outer-loop reference control and inputting them into Equation (44) yielded:
{ e ˙ x = a 1 tanh ( ρ 1 e x ) + cos γ d cos χ d ( v i v d ) + v i cos γ i ( cos χ i cos χ d ) + v i cos χ d ( cos γ i cos γ d ) e ˙ y = a 2 tanh ( ρ 2 e y ) + cos γ d sin χ d ( v i v d ) + v i cos γ i ( sin χ i sin χ d ) + v i sin χ d ( cos γ i cos γ d ) e ˙ z = a 3 tanh ( ρ 3 e z ) sin γ d ( v i v d ) v i ( sin γ i sin γ d )
The following is an example of analyzing the stability of the system with e x . the system Lyapunov function can be taken as:
V = 3 a 1 ln ( cosh ρ 1 e x ) + 1 2 ρ 1 e x 2
where α , ρ 1 > 0 . Taking the derivative of Equation (46) leads to,
V ˙ = 3 a 1 ρ 1 e ˙ x tanh ρ 1 e x + ρ 1 e x e ˙ x
Let m 1 = a 1 tanh ( ρ 1 e x ) , m 2 = cos γ d cos χ d ( v i v d ) , m 3 = v i cos γ i ( cos χ i cos χ d ) , m 4 = v i cos χ d ( cos γ i cos γ d ) , then e ˙ x can be expressed as: e ˙ x = m 1 + m 2 + m 3 + m 4 , and, when substituted into Equation (47), gives:
V ˙ e x = 3 ρ 1 m 1 ( m 1 + m 2 + m 3 + m 4 ) + ρ 1 e x ( m 1 + m 2 + m 3 + m 4 ) = ρ 1 [ m 1 2 3 m 1 m 2 + ( 3 2 m 2 ) 2 ] ρ 1 [ m 1 2 3 m 1 m 2 + ( 3 2 m 2 ) 2 ] ρ 1 [ m 1 2 3 m 1 m 2 + ( 3 2 m 2 ) 2 ] ρ 1 e x m 1 + 9 4 ρ 1 m 2 ( m 2 + e x ) + 9 4 ρ 1 m 3 ( m 3 + e x ) + 9 4 ρ 1 m 3 ( m 3 + e x )
Additionally, because of a 1 tanh ( ρ 1 e x ) > 0 , therefore:
V ˙ e x < ρ 1 [ m 1 2 3 m 1 m 2 + ( 3 2 m 2 ) 2 ] ρ 1 [ m 1 2 3 m 1 m 2 + ( 3 2 m 2 ) 2 ] ρ 1 [ m 1 2 3 m 1 m 2 + ( 3 2 m 2 ) 2 ] + 9 4 ρ 1 m 2 ( m 2 + e x ) + 9 4 ρ 1 m 3 ( m 3 + e x ) + 9 4 ρ 1 m 3 ( m 3 + e x ) < ρ 1 [ ( m 1 3 2 m 2 ) 2 + ( m 1 3 2 m 3 ) 2 + ( m 1 3 2 m 4 ) 2 ] + 9 4 ρ 1 m 2 ( m 2 + e x ) + 9 4 ρ 1 m 3 ( m 3 + e x ) + 9 4 ρ 1 m 3 ( m 3 + e x )
Additionally, the following can be considered:
| cos χ i cos χ d | = | 2 sin χ i + χ d 2 sin χ i χ d 2 | 2 | sin χ i χ d 2 | | χ i χ d |
The same reasoning leads to:
| cos γ i cos γ d | | γ i γ d |
From Equation (43), v i , γ i , χ i is exponentially convergent. For ε > 0 , t v i > 0 , when t > t v i , | v i v d | < ε ; similarly, this applies when t > t χ i , | χ i χ d | < ε and when t > t γ i , | γ i γ d | < ε . So, there exists a finite time t e x , when t e x > max { t v i , t χ i , t γ i } , V ˙ e x < 0 and when t , m i 0 , ( i = 2 , 3 , 4 ) . So when t , e x 0 , Similarly, it can be proved that when t , e y 0 , e z 0 , and the global stability of the system is proved.

3.3. Delayed System Stability Proof

When the navigator exists within the field of view of the follower UAV, the system can achieve tracking in a time-delay-free state. However, when the navigator is not present in the field of view of the follower UAV, there is a time delay in the state transfer, and Theorem II is given below and demonstrates that the delayed system can achieve convergence.
Theorem 2. 
When the time delay of system  lim t x l 4 = x r is  h , satisfying h > 0 , h ( 0 , h ¯ ] , h ¯ appears as the upper limit of delay, then when:
Γ ( h ) = [ Q 1 0 0 P 1 A 12 r 0 Q 1 0 0 0 0 h Z 1 0 A 12 r T P 1 0 0 A 12 r T h Z 1 A 12 r + P 2 A 22 r + A 22 r T P 2 + Q 2 + A 22 r T h Z 2 A 22 r A 12 d T P 0 0 A 12 d T h Z 1 A 12 r + A 22 d T P 2 + A 22 d T h Z 2 A 22 r 0 0 0 0 0 0 0 B 2 r T P 2 + B 2 r T h Z 2 A 22 r 0 0 0 B 2 d T P 2 + B 2 d T h Z 2 A 22 r P 1 A 12 d 0 0 0 0 0 0 0 0 0 0 0 A 12 r T h Z 1 A 12 d + P 2 A 22 d + A 22 r T h Z 2 A 22 d 0 P 2 B 2 r + A 22 r T h Z 2 B 2 r P 2 B 2 d + A 22 r T h Z 2 B 2 d A 12 d T h Z 1 A 12 d Q 2 + A 22 d T h Z 2 A 22 d 0 A 22 d T h Z 2 B 2 r A 22 d T h Z 2 B 2 d 0 h Z 2 0 0 B 2 r T h Z 2 A 22 d 0 R + B 2 r T h Z 2 B 2 r B 2 r T h Z 2 B 2 d B 2 d T h Z 2 A 22 d 0 B 2 d T h Z 2 B 2 r B 2 d T h Z 2 B 2 d R ] < 0
this indicates delayed system stability. Among others, P 1 , Z 1 , P 2 , Q , R , Z 2 > 0 .
Proof of Theorem 2. 
First, linearizing the error state Equation (32) results in:
e ˙ i = A e i + B u i
where:
A = [ A 11 A 12 A 13 A 14 ] = [ 0 0 0 cos γ r cos χ r v r cos γ r sin χ r v r sin γ r cos χ r 0 0 0 cos γ r sin χ r v r cos γ r cos χ r v r sin γ r sin χ r 0 0 0 sin γ r 0 v r cos γ r 0 0 0 0 0 0 0 0 0 1 / v r 2 0 0 0 0 0 1 / v r 2 0 0 ] ,   B = [ B 11 B 21 ] = [ 0 0 0 0 0 0 0 0 0 1 0 0 0 1 / v r 0 0 0 1 / v r ]
namely,
e ˙ 1 ( t ) = A 12 r e 2 ( t ) + A 12 d e 2 ( t h ) e ˙ 2 ( t ) = A 22 r e 2 ( t ) + A 22 d e 2 ( t h ) + B 2 r u ( t ) + B 2 d u ( t h )
where A 12 r , A 22 r , B 2 r is the ideal state error matrix and A 12 d , A 22 d , B 2 d is the error between the ideal state matrix and the actual matrix. The Lyapunov–Krasovskii generalized function can be taken as follows.
V 3 ( e 2 , t ) = e 1 T ( t ) P 1 e 1 ( t ) V 4 ( e 2 , t ) = t h t e 1 T ( α ) Q 1 e 1 ( α ) d α V 5 ( e 2 , t ) = h 0 t h t e ˙ 1 T ( t ) Z 1 e ˙ 1 ( t ) d α d β V 6 ( e 2 , t ) = e 2 T ( t ) P 2 e 2 ( t ) V 7 ( e 2 , t ) = t h t e 2 T ( α ) Q 2 e 2 ( α ) d α V 8 ( u 2 , t ) = t h t u 2 T ( t ) R u 2 ( t ) d α V 9 ( e 2 , t ) = h 0 t h t e ˙ 2 T ( t ) Z 2 e ˙ 2 ( t ) d α d β
where P 1 , Z 1 , P 2 , Q , R , Z 2 > 0 , and it is clear that the above equation satisfies this.
V ( e , t ) = i = 3 8 V i ( e , t ) > 0
From the system (52), the derivative of Equation (53) can be obtained as:
V ˙ 3 ( e 1 , t ) = 2 e 1 T ( t ) P 1 e ˙ 1 ( t ) = 2 e 1 T ( t ) P 1 A 12 r e 2 ( t ) + 2 e 1 T ( t ) P 1 A 12 d e 2 ( t h ) + 2 e 1 T ( t ) P 1 w 1 ( t ) = 1 h t h t [ 2 e 1 T ( t ) P 1 A 12 r e 2 ( t ) + 2 e 1 T ( t ) P 1 A 12 d e 2 ( t h ) + 2 e 1 T ( t ) P 1 w 1 ( t ) ] d α
V ˙ 4 ( e 1 , t ) = 1 h t h t [ e 1 T ( t ) Q 1 e 1 ( t ) e 1 T ( t h ) Q 1 e 1 ( t h ) ] d α
V ˙ 5 ( e 1 , t ) = h 0 t h t e ˙ 1 T ( t ) Z 1 e ˙ 1 ( t ) d α d β = t h t [ e ˙ 1 T ( t ) Z 1 e ˙ 1 ( t ) e ˙ 1 T ( t + β ) Z 1 e ˙ 1 ( t + β ) ] d β = t + β = α t h t [ e ˙ 1 T ( t ) Z 1 e ˙ 1 ( t ) e ˙ 1 T ( α ) Z 1 e ˙ 1 ( α ) ] d α = 1 h t h t h [ ( A 12 r e 2 ( t ) + A 12 d e 2 ( t h ) + w 1 ( t ) ) T Z 1 ( A 12 r e 2 ( t ) + A 12 d e 2 ( t h ) + w 1 ( t ) ) e ˙ 1 T ( α ) Z 1 e ˙ 1 ( α ) ] d α
V ˙ 6 ( e 2 , t ) = 2 e 2 T ( t ) P 2 e ˙ 2 ( t ) = 2 e 2 T ( t ) P 2 A 22 r e 2 ( t ) + 2 e 2 T ( t ) P 2 A 22 d e 2 ( t h ) + 2 e 2 T ( t ) P 2 B 2 r u ( t ) + 2 e 2 T ( t ) P 2 B 2 d u ( t h ) = 1 h t h t [ 2 e 2 T ( t ) P 2 A 22 r e 2 ( t ) + 2 e 2 T ( t ) P 2 A 22 d e 2 ( t h ) + 2 e 2 T ( t ) P 2 B 2 r u ( t ) + 2 e 2 T ( t ) P 2 B 2 d u ( t h ) ] d α
V ˙ 7 ( e 2 , t ) = 1 h t h t [ e 2 T ( t ) Q 2 e 2 ( t ) e 2 T ( t h ) Q 2 e 2 ( t h ) ] d α
V ˙ 8 ( u 2 , t ) = 1 h t h t [ u 2 T ( t ) R u 2 ( t ) u 2 T ( t h ) R u 2 ( t h ) ] d α
V ˙ 9 ( e 2 , t ) = h 0 t h t e ˙ 2 T ( t ) Z 2 e ˙ 2 ( t ) d α d β = t h t [ e ˙ 2 T ( t ) Z 2 e ˙ 2 ( t ) e ˙ 2 T ( t + β ) Z 2 e ˙ 2 ( t + β ) ] d β = t + β = α t h t [ e ˙ 2 T ( t ) Z 2 e ˙ 2 ( t ) e ˙ 2 T ( α ) Z 2 e ˙ 2 ( α ) ] d α = 1 h t h t h [ ( A 22 r e 2 ( t ) + A 22 d e 2 ( t h ) + B 2 r u ( t ) + B 2 d u ( t h ) ) T h Z 2 ( A 22 r e 2 ( t ) + A 22 d e 2 ( t h ) + B 2 r u ( t ) + B 2 d u ( t h ) ) e ˙ 2 T ( α ) Z 2 e ˙ 2 ( α ) ] d α
From Equations (55)–(61), it can be obtained that:
V ˙ ( e , t ) = 1 h t h t [ Ξ T ( t , α ) Γ ( h ) Ξ ( t , α ) ] d α
where:
Ξ ( t , α ) = [ e 1 ( t ) e 1 ( t h ) e ˙ 1 ( α ) e 2 ( t ) e 2 ( t h ) e ˙ 2 ( α ) u ( t ) u ( t h ) ] ,
So, when Γ ( h ) < 0 , V ˙ ( e , t ) < 0 , theorem II can be proved. □

4. Simulation and Test Results

In order to verify the algorithm proposed in this paper, the control algorithm proposed was simulated with MATLAB, as outlined in this section. We assumed that the UAVs used for simulation carried onboard sensors that were capable of measuring the relative distance, relative pitch, and relative yaw angles. The simulation parameters were set as shown in Table 1. The initial states of the drones are shown in Table 2.

4.1. Case 1

In case 1, we first simulated the case of four UAVs in formation in a simple curve environment. The simulation results are shown in Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9. We set the desired trajectory as:
[ x d y d z d ] T = [ t sin ( 0.5 × x d ) + 0.5 × x d + 0.5 x d ] T
As can be seen in Figure 4, the UAV is able to achieve the expected formation control and maintain the stability of the formation until the end of the simulation, verifying the effectiveness of the algorithm in simple curved paths. Additionally, when the leader is not in the follower’s field of view, and there is a certain delay in state transfer, No. 4 is able to follow the leader effectively. In Figure 5, at the first sampling moment, the solution time of the UAV controller of No. 1 was significantly higher than that of the other UAVs, but it was less than 0.01 s, verifying that the controller could complete the solution within the sampling time to avoid the errors caused by an untimely solution. Figure 6 shows the acceleration input of the UAV during the actual motion, and it can be seen that the algorithm of this controller avoided the jitter of the input. Figure 7 gives the virtual input of the system calculated based on the UAV position error, from which it can be seen that the virtual input was able to generate a continuous inner-loop control signal (as shown in Figure 8 for the inner-loop ideal tracking state). Regardless of the UAV state volume (as shown in Figure 8) or the position error (as shown in Figure 9), the following UAV was able to achieve finite time convergence, while the inner loop convergence rate was greater than the outer loop convergence rate, ensuring the convergence of the system.

4.2. Case 2

In case 2, we simulated the formation control and retention problem of four UAVs in a complex path, and the simulation results are shown in Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13. We set the desired trajectory as:
[ x d y d z d ] T = [ t 5 × sin ( 0.5 × x d ) 5 × cos ( 0.5 × x d ) + 0.5 × x d ] T
In Figure 10 of case 2, it can be seen that the four UAVs were still able to achieve the formation and maintenance of the formation under complex path conditions, but the tracking effect was not particularly good in the initial phase when the UAVs started tracking, and further improvement is needed in the next step of work. From Figure 11, it can be seen that the No. 1 UAV and No. 2 UAV had a long controller calculation time in the first sampling moment, and the controller solution time within the cluster was less than 0.01 s. In Figure 12, it can be seen that the inputs of the four UAVs were still able to achieve formation control and acceleration convergence in a short time in a complex environment. Additionally, as can be seen in Figure 13, when the path was continuous, the generated virtual inputs were also continuous and could provide continuous desired signals to the inner loop control. In Figure 14 and Figure 15, it can be seen that the state, as well as the position error of the UAV, can achieve fast convergence under the action of the algorithm, and the inner loop converged faster than the outer loop converges, which ensured the formation and maintenance of the formation.

4.3. Case 3

To further verify the effectiveness of the algorithm for the formation control of large-scale clusters, the formation control of eight UAVs was simulated in a complex environment with the initial states of the UAVs shown in Table 3, and the simulation results are shown in Figure 14–18. In case 3, we set the desired trajectory as:
[ x d y d z d ] T = [ t 5 × sin ( 0.5 × x d ) 5 × cos ( 0.5 × x d ) + 0.5 × x d ] T
In Figure 16, it can be seen that the eight UAVs were able to effectively achieve and maintain an ideal formation under the action of the algorithm, but there was also a situation where some of the UAVs were poorly tracked during the initial phase, which required correcting in the next step to ensure that the system could achieve stability from the initial state under various paths. In Figure 17, it is shown that among the cluster of 8 UAVs, there existed 3 UAVs (UAVs No. 1,3,5) whose first solving time is much longer than the other sorties, and all other numbered UAVs were able to satisfy the solving of the controller within one sampling moment. In Figure 18, it can be seen that the acceleration of each UAV experienced large fluctuations in the initial state, and after the system error started to converge, the system input gradually converged, and the amplitude decreased, ensuring the formation control at low acceleration. Figure 19 then shows that the virtual inputs of each UAV gradually converged after stabilization and were able to provide a stable and continuous tracking signal to the inner loop (as shown in Figure 20). In Figure 20 and Figure 21, the UAV inner-loop state error, as well as the outer-loop position error are given, respectively; from this, it can be seen that the large-scale cluster could achieve the convergence of the state error as well as formation control in the complex path case.

4.4. Discussion

From the above three numerical simulation examples, it can be seen that:
  • The designed controller can be effectively applied to the formation control of fixed-wing UAVs, even in the presence of measurement errors as well as state transfer time delays in multiple UAV formations.
  • It can be seen from the computation times of the three examples that, since the controller was distributed, it did not increase the solution time of individual UAVs as the number of UAVs increased. However, there was a long initial moment at the beginning of each phase in the time of each numerical simulation, which would be the next step to improve.
  • When there were multiple UAVs in the neighborhood, the selection of the follower UAV for the desired leader UAV became an issue. In the next step of the study, rules need to be developed or the research method updated to ensure the efficiency and accuracy of the follower UAV’s selection of the desired leader UAV.

5. Conclusions

In this paper, a global sliding mode control algorithm for the vision-based control of multiple UAV formations was proposed. In the presence of sensor errors, the follower UAV uses extended Kalman filtering to achieve an estimate of the position as well as the state of the leader UAV in the field of view, and the follower UAV uses the estimated position as a reference value for the local outer loop control. After that, the continuous signal of the inner loop control reference is generated using the outer loop virtual control, which is then passed to the inner loop. Later, the inner-loop controller was designed to ensure the tracking of the inner-loop reference signal, and it was mathematically proven that the inner-loop controller converged faster than the outer loop while the stability of the global algorithm was demonstrated. Finally, the stability of the formation in the presence of time delay was given considering the time delay in the transmission of the ideal state of the UAV due to the limited range of action of the UAV sensors. In a future study, this algorithm will be further investigated in its initial state for the tracking of the reference value, the control of a larger number of clusters, and the further improvement of the algorithm in the case of poor initial state tracking.

Author Contributions

Conceptualization, Q.C. and Y.J.; methodology, Q.C.; validation, Q.C., T.W. and B.Q.; formal analysis, Q.C. and T.W.; investigation, Y.W.; resources, Y.W.; data curation, Y.W. and B.Q.; writing—original draft preparation, Q.C. and Y.J.; writing—review and editing, Y.W. and B.Q.; visualization, T.W.; supervision, T.W. and Y.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data used in the experiments are available in the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alsamhi, S.H.; Shvetsov, A.V.; Kumar, S.; Shvetsova, S.V.; Alhartomi, M.A.; Hawbani, A.; Rajput, N.S.; Srivastava, S.; Saif, A.; Nyangaresi, V.O. UAV Computing-Assisted Search and Rescue Mission Framework for Disaster and Harsh Environment Mitigation. Drones 2022, 6, 154. [Google Scholar] [CrossRef]
  2. Radmanesh, M.; Kumar, M. Flight formation of UAVs in presence of moving obstacles using fast-dynamic mixed integer linear programming. Aerosp. Sci. Technol. 2016, 50, 149–160. [Google Scholar] [CrossRef] [Green Version]
  3. Madyastha, V.K.; Caliset, A.J. An adaptive filtering approach to target tracking. In Proceedings of the American Control Conference, Portland, OR, USA, 8–10 June 2005. [Google Scholar]
  4. Alsamhi, S.H.; Shvetsov, A.V.; Kumar, S.; Hassan, J.; Alhartomi, M.A.; Shvetsova, S.V.; Sahal, R.; Hawbani, A. Computing in the Sky: A Survey on Intelligent Ubiquitous Computing for UAV-Assisted 6G Networks and Industry 4.0/5.0. Drones 2022, 7, 177. [Google Scholar] [CrossRef]
  5. Alsamhi, S.H.; Ma, O.; Ansari, M.S.; Gupta, S.K. Collaboration of Drone and Internet of Public Safety Things in Smart Cities: An Overview of QoS and Network Performance Optimization. Drones 2019, 3, 13. [Google Scholar] [CrossRef] [Green Version]
  6. Desai, J.P.; Ostrowski, J.P.; Kumar, V. Modeling and control of formations of nonholonomic mobile robots. IEEE Trans. Robot. Autom. 2001, 17, 905–908. [Google Scholar] [CrossRef] [Green Version]
  7. Roldão, V.; Cunha, R.; Cabecinhas, D.; Silvestre, C.; Oliveira, P. A leader-following trajectory generator with application to quadrotor formation flight. Robot. Auton. Syst. 2014, 62, 1597–1609. [Google Scholar] [CrossRef]
  8. Fierro, R.; Belta, C.; Desai, J.P.; Kumar, V. On controlling aircraft formations. In Proceedings of the 40th IEEE Conference on Decision and Control (Cat. No.01CH37228), Orlando, Fl, USA, 4–7 December 2001; Volume 2, pp. 1065–1070. [Google Scholar]
  9. Eren, T. Formation shape control based on bearing rigidity. Int. J. Control. 2012, 85, 1361–1379. [Google Scholar] [CrossRef]
  10. Lewis, M.A.; Tan, K.H. High precision formation control of mobile robots using virtual structures. Auton. Robots. Auton. Robot. 1997, 4, 387–403. [Google Scholar] [CrossRef]
  11. Sharma, R.; Ghose, D. Collision avoidance between UA V clusters using swarm intelligence techniques. Int. J. Syst. Sci. 2009, 5, 521–538. [Google Scholar] [CrossRef]
  12. Giulietti, F.; Innocenti, M.; Pollini, L. Formation flight control—A behavioral approach. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Montreal, QC, Canada, 6–9 August 2001. [Google Scholar]
  13. Park, C.; Tahk, M.; Bang, H. Multiple Aerial Vehicle Formation Using Swarm Intelligence. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, Austin, TX, USA, 11–14 August 2003. [Google Scholar]
  14. Ren, W. Consensus based formation control strategies for multi-vehicle systems. In Proceedings of the 2006 American Control Conference, Minneapolis, MN, USA, 14–16 June 2006. [Google Scholar]
  15. Jadbabaie, A.; Lin, J.; Morse, A.S. Coordination of groups of mobile autonomous agents using nearest neighbor rules. IEEE Trans. Autom. 2003, 6, 988–1001. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, X.; Yu, Y.; Li, Z. Distributed sliding mode control for leader-follower formation flight of fixed-wing unmanned aerial vehicles subject to velocity constraints. Int. J. Robust Nonlinear Control 2021, 31, 2110–2125. [Google Scholar] [CrossRef]
  17. Liu, S.; Tan, D.; Liu, G. Robust Leader-follower Formation Control of Mobile Robots Based on a Second Order Kinematics Model. Acta Autom. Sin. 2007, 33, 947–955. [Google Scholar] [CrossRef] [Green Version]
  18. Zhao, C.; Dai, S.; Zhao, G.; Gao, C.; Liu, S. UAV formation control based on distributed model predictive control. Control. Decis. 2021, 1, 1–9. [Google Scholar]
  19. Kim, S.; Cho, H.; Jung, D. Evaluation of Cooperative Guidance for Formation Flight of Fixed-wing UAVs using Mesh Network. In Proceedings of the 2020 International Conference on Unmanned Aircraft Systems, Athens, Greece, 1–4 September 2020. [Google Scholar]
  20. Parivallal, A.; Sakthivel, R.; Amsaveni, R.; Alzahrani, F.; Alshomrani, A.S. Observer-based memory consensus for nonlinear multi-agent systems with output quantization and Markov switching topologies. Phys. A Stat. Mech. Its Appl. 2020, 551, 123949. [Google Scholar] [CrossRef]
  21. Sakthivel, R.; Parivallal, A.; Huy Tuan, N.; Manickavalli, S. Nonfragile control design for consensus of semi-Markov jumping multiagent systems with disturbances. Int. J. Adapt. Control. Signal Process. 2021, 35, 1039–1061. [Google Scholar] [CrossRef]
  22. Polloni, L.; Mati, R.; Innocenti, M.; Campa, G.; Napolitano, M. A Synthetic Environment for the Simulation of Visionbased Formation Flight. In Proceedings of the AIAA Modeling and Simulation Technologies Conference, Austin, TX, USA, 11–14 August 2003. [Google Scholar]
  23. Ray, R.; Cobleigh, B.; Vachon, M.; St John, C. Flight Test Techniques used to Evaluate Performance Benefits During Formation Flight. In Proceedings of the AIAA Atmospheric Flight Mechanics Conference and Exhibit, Monterey, CA, USA, 5–8 August 2002. [Google Scholar]
  24. Seanor, B.; Campa, G.; Gu, Y.; Napolitano, M.; Rowe, L.; Perhinschi, M. Formation Flight Test Results for UAV Research Aircraft Models. In Proceedings of the AIAA 1st Intelligent Systems Technical Conference, Chicago, IL, USA, 20–22 September 2004. [Google Scholar]
  25. Sattigeri, R.; Calise, A.; Kim, B.S.; Volyanskyy, K.; Kim, N. 6-DOF Nonlinear Simulation of Vision-Based Formation Flight. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, San Francisco, CA, USA, 15–18 August 2005. [Google Scholar]
  26. Lin, J.; Miao, Z.; Zhong, H.; Peng, W.; Wang, Y.; Fierro, R. Adaptive Image-Based Leader-Follower Formation Control of Mobile Robots with Visibility Constraints. IEEE Trans. Ind. Electron. 2021, 68, 6010–6019. [Google Scholar] [CrossRef]
  27. Chen, X.; Jia, Y. Adaptive leader-follower formation control of nonholonomic mobile robots using active vision. IET Control. Theory Appl. 2015, 8, 1302–1311. [Google Scholar] [CrossRef]
  28. Sattigeri, R.; Calise, A.J.; Evers, J.H. An adaptive vision-based approach to decentralized formation control. J. Aerosp. Comput. Inf. Commun. 2004, 1, 502–525. [Google Scholar] [CrossRef] [Green Version]
  29. Chen, Q.; Jin, Y.; Wang, T.; Wang, Y.; Yan, T.; Long, Y. UAV Formation Control Under Communication Constraints Based on Distributed Model Predictive Control. IEEE Access 2022, 2022, 126494–126507. [Google Scholar]
  30. Qin, J.; Zhang, G.; Zheng, W.X.; Kang, Y. Adaptive Sliding Mode Consensus Tracking for Second-Order Nonlinear Multiagent Systems with Actuator Faults. IEEE Trans. Cybern. 2019, 49, 1605–1615. [Google Scholar] [CrossRef]
  31. Desai, J.P.; Ostrowski, J.; Kumar, V. Controlling formations of multiple mobile robots. In Proceedings of the IEEE Transactions on Robotics & Automation, Leuven, Belgium, 20 May 2010; Volume 6, pp. 905–908. [Google Scholar]
  32. Das, A.K.; Fierro, R.; Kumar, V.; Ostrowski, J.P.; Spletzer, J.; Taylor, C.J. A vision-based formation control framework. IEEE Trans. Robot. Autom. 2002, 18, 813–825. [Google Scholar] [CrossRef] [Green Version]
  33. Mariottini, G.L.; Morbidi, F.; Prattichizzo, D.; Vander Valk, N.; Michael, N.; Pappas, G.; Daniilidis, K. Vision-Based Localization for Leader-Follower Formation Control. IEEE Trans. Robot. 2009, 25, 1431–1438. [Google Scholar] [CrossRef]
  34. Galzi, D.; Shtessel, Y. UAV formations control using high order sliding modes. In Proceedings of the American Control Conference, Minneapolis, MN, USA, 14–16 June 2006. [Google Scholar]
  35. Kadungoth Sreeraj, A.R. A Model-Free Control Algorithm Based on the Sliding Mode Control Method with Applications to Unmanned Aircraft Systems. Ph.D. Thesis, Rochester Institute of Technology, Rochester, NY, USA, 2019. [Google Scholar]
  36. Mobayen, S. Adaptive Global Terminal Sliding Mode Control Scheme with Improved Dynamic Surface for Uncertain Nonlinear Systems. Int. J. Control. Autom. Syst. 2018, 16, 1692–1700. [Google Scholar] [CrossRef]
  37. Dehghani, M.A.; Menhaj, M.B. Integral sliding mode formation control of fixed-wing unmanned aircraft using seeker as a relative measurement system. Aerosp. Sci. Technol. 2016, 58, 318–327. [Google Scholar] [CrossRef]
  38. Ailon, A.; Zohar, I. Controllers for trajectory tracking and string-like formation in Wheeled Mobile Robots with bounded inputs. In Proceedings of the Melecon 2010–2010 15th IEEE Mediterranean Electrotechnical Conference, Valletta, Malta, 26–28 April 2010. [Google Scholar]
Figure 1. Coordinate frames.
Figure 1. Coordinate frames.
Drones 07 00231 g001
Figure 2. The relative kinematics of the leader-follower.
Figure 2. The relative kinematics of the leader-follower.
Drones 07 00231 g002
Figure 3. Schematic diagram of algorithm principle.
Figure 3. Schematic diagram of algorithm principle.
Drones 07 00231 g003
Figure 4. UAV trajectory.
Figure 4. UAV trajectory.
Drones 07 00231 g004
Figure 5. Single UAV solving time.
Figure 5. Single UAV solving time.
Drones 07 00231 g005
Figure 6. Actual UAV input.
Figure 6. Actual UAV input.
Drones 07 00231 g006
Figure 7. UAV virtual input.
Figure 7. UAV virtual input.
Drones 07 00231 g007
Figure 8. Attitude error.
Figure 8. Attitude error.
Drones 07 00231 g008
Figure 9. Position error.
Figure 9. Position error.
Drones 07 00231 g009
Figure 10. UAV trajectory.
Figure 10. UAV trajectory.
Drones 07 00231 g010
Figure 11. Single UAV solving time.
Figure 11. Single UAV solving time.
Drones 07 00231 g011
Figure 12. Actual UAV input.
Figure 12. Actual UAV input.
Drones 07 00231 g012
Figure 13. UAV virtual input.
Figure 13. UAV virtual input.
Drones 07 00231 g013
Figure 14. Attitude error.
Figure 14. Attitude error.
Drones 07 00231 g014
Figure 15. Position error.
Figure 15. Position error.
Drones 07 00231 g015
Figure 16. UAV trajectory.
Figure 16. UAV trajectory.
Drones 07 00231 g016
Figure 17. Single UAV solving time.
Figure 17. Single UAV solving time.
Drones 07 00231 g017
Figure 18. Actual UAV input.
Figure 18. Actual UAV input.
Drones 07 00231 g018
Figure 19. UAV virtual input.
Figure 19. UAV virtual input.
Drones 07 00231 g019
Figure 20. Attitude error.
Figure 20. Attitude error.
Drones 07 00231 g020
Figure 21. Position error.
Figure 21. Position error.
Drones 07 00231 g021
Table 1. Simulation parameters table.
Table 1. Simulation parameters table.
Parameter NameSymbolsNumerical Value
Sampling interval T k 0.01
Sampling time T 30
Number of drones N v 4
Intermachine delay h 0.004
Table 2. UAV initial parameter table.
Table 2. UAV initial parameter table.
UAV NumberInitial Location v ( m / s ) Yaw Angle
1 [ 2 , 2 , 0 ] T 0 π 3
2 [ 2 , 3 , 0 ] T 0 π 3
3 [ 2 , 1 , 0 ] T 0 π 3
4 [ 2 , 1 , 0 ] T 0 π 3
Table 3. UAV initial parameter table.
Table 3. UAV initial parameter table.
UAV NumberInitial Location v ( m / s ) Yaw Angle
1 [ 2 , 2 , 0 ] T 0 π 3
2 [ 2 , 3 , 0 ] T 0 π 3
3 [ 2 , 1 , 0 ] T 0 π 3
4 [ 2 , 1 , 0 ] T 0 π 3
5 [ 2 , 2 , 0 ] T 0 π 3
6 [ 2 , 3 , 0 ] T 0 π 3
7 [ 2 , 1 , 0 ] T 0 π 3
8 [ 2 , 2 , 0 ] T 0 π 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, Q.; Wang, T.; Jin, Y.; Wang, Y.; Qian, B. A UAV Formation Control Method Based on Sliding-Mode Control under Communication Constraints. Drones 2023, 7, 231. https://doi.org/10.3390/drones7040231

AMA Style

Chen Q, Wang T, Jin Y, Wang Y, Qian B. A UAV Formation Control Method Based on Sliding-Mode Control under Communication Constraints. Drones. 2023; 7(4):231. https://doi.org/10.3390/drones7040231

Chicago/Turabian Style

Chen, Qijie, Taoyu Wang, Yuqiang Jin, Yao Wang, and Bei Qian. 2023. "A UAV Formation Control Method Based on Sliding-Mode Control under Communication Constraints" Drones 7, no. 4: 231. https://doi.org/10.3390/drones7040231

APA Style

Chen, Q., Wang, T., Jin, Y., Wang, Y., & Qian, B. (2023). A UAV Formation Control Method Based on Sliding-Mode Control under Communication Constraints. Drones, 7(4), 231. https://doi.org/10.3390/drones7040231

Article Metrics

Back to TopTop