Adaptive Robust Trajectory Tracking Control of Multiple Quad-Rotor UAVs with Parametric Uncertainties and Disturbances

Recently, formation flying of multiple unmanned aerial vehicles (UAVs) found numerous applications in various areas such as surveillance, industrial automation and disaster management. The accuracy and reliability for performing group tasks by multiple UAVs is highly dependent on the applied control strategy. The formation and trajectories of multiple UAVs are governed by two separate controllers, namely formation and trajectory tracking controllers respectively. In presence of environmental effects, disturbances due to wind and parametric uncertainties, the controller design process is a challenging task. This article proposes a robust adaptive formation and trajectory tacking control of multiple quad-rotor UAVs using super twisting sliding mode control method. In the proposed design, Lyapunov function-based adaptive disturbance estimators are used to compensate for the effects of external disturbances and parametric uncertainties. The stability of the proposed controllers is guaranteed using Lyapunov theorems. Two variants of the control schemes, namely fixed gain super twisting SMC (STSMC) and adaptive super twisting SMC (ASTSMC) are tested using numerical simulations performed in MATLAB/Simulink. From the results presented, it is verified that in presence of disturbances, the proposed ASTSMC controller exhibits enhanced robustness as compared to the fixed gain STSMC.


Introduction
A flight in which more than one quad-rotors fly and maintain the relative distance among each other is called formation flight. Recently, the interest in the formation control of quad-rotors has attracted a lot of attention. This trend is due to its potential applications in the defense industry, aerial mapping, search and rescue operations, oil fields monitoring, agriculture and transportation of suspended loads [1]. It is expected that by 2027 the payload market value of global UAV may reach USD 3 billion due to its anticipated usefulness [2]. Multiple quad-rotors increase the capacity for equipping sensors, provides larger payload capacity and a wider surveillance area as compared to single quad-rotor [3][4][5][6][7]. However, controlling the formation of multiple quad-rotors in the presence of uncertainties is a challenging task. Moreover, the derivation of formation dynamic model for multiple UAVs in presence of external perturbations has also become an important topic . The transnational and rotational dynamics of a quad rotor is modeled as six degree of freedom nonlinear differential equations [8][9][10][11]. For multiple UAVs, different formation geometries exist depending upon the number of quad-rotors and the purpose of flight. These include V shape geometry and finger four geometry. V shape geometry is used in this paper for flight formation of three quad-rotors.The quad-rotors formation can be stepped up, stepped down, and leveled based on the altitude position between leader and follower. A leveled formation is used in this paper. Since a quad-rotor is a complex system with under-actuated multi variable non-linear model hence its formation control problem is more difficult to control. To ensure robust formation control, the modeling uncertainties and the disturbance due to wind gust must be compensated using appropriate control schemes.
To solve the formation control problem of the quad rotor, many research efforts have been made. In [12] the authors proposed a leader-follower formation control using classical proportional derivative scheme and fuzzy logic system for the formation pattern. However the above controller did not take uncertainties into considerations. In [13], a prescribed performance controller is proposed for formation control of multiple UAVs. The prescribed performance controller ensures robust formation pattern and trajectory. In [14], a classical PID control scheme with a sliding mode controller (SMC) is proposed for multiple quad-rotors. However the above mentioned PID-SMC controller ignores the disturbances and communication delays between multiple UAVs. In [15][16][17], classical SMC method is proposed for formation control problems of multiple UAVs, however the classical SMC method offers high frequency chattering in the excitation signal. Chattering phenomena degrades the life of the actuators.A control scheme for the circular formation scheme of multiple UAVs is presented in [18]. A classical PI control-based synchronization control for the formation of two UAVs is presented in [19]. In [20] a distributed controller is presented to compensate the communication delays in multiple UAVs formation. Similarly a nonlinear distributed controller is proposed for formation of micro UAVs [21]. In [22], the authors proposed a cohesive formation controller multiple UAVs.A back-stepping control scheme is proposed for the formation control of multiple UAVs [23]. Similarly a model predictive control scheme is proposed for multiple UAVs using adaptive gain tuning method in [24,25]. Velocity tracking and formation control of quad-rotors is achieved by the design of prior-bounded intermediary adaptive controller which gives the reference orientation and bounded control thrust [26]. A guidance algorithm based on Lyapunov function is used for the formation control of quad-rotors with attached slung load, where the quad-rotors are controlled using linear quadratic tracking controller [27]. Leader follower formation controller is designed for two parrot drones in [28], where a proportional derivative controller is implemented in the respective models. The problems related to formation and tracking control of quad-rotors in leader follower formation are addressed and a formation controller is designed to avoid collision in swarm [29]. Adaptive law for the formation control of swarm UAVs in leader follower mechanism influenced by motion constraints and and unknown external disturbances is discussed in [30]. For a swarm of three quad-rotors, a model reference adaptive control algorithm is presented in [31]. The controller gains are tuned online, by which the algorithm allows the system to adapt to unexpected disturbances. However, in this method, no robust controller is investigated. A semi physical platform for formation control of multiple fixed wings UAV is proposed in [32]. In [33], a detailed survey on low cost UAV platforms for infra structure monitoring is proposed.
The above cited work is specifically focused on the formation control of multiple UAVs. It is also necessary to describe the back ground of robust control system due to its utmost importance in control community. Robust control is designed for uncertain systems in which the uncertainty belongs to bounded set [34]. Robust controllers are designed both in frequency and time domains. A widely used frequency domain robust controller is the H ∞ method and it was first reported in [34]. Later on several variants of H ∞ control were reported in the literature such as loop shaping in [35], optimal H ∞ control using Riccati equations [36] and Linear Matrix Inequalities (LMIs)-based design in [37]. In frequency domain, the performance and stability of the control system is measured in terms of gain and phase margins, percent overshoot, rise, delay and settling times of the signals. Apart from frequency domain, a modern approach for designing robust controllers is the state space frame work. Sliding mode-based control (SMC) system variants are the most widely used methods and it find numerous applications in all areas of science and technology [38]. Classical SMC has several disadvantages such as high frequency chattering and asymptotic convergence property. These shortcoming are addressed by introducing new variants of SMC such as global sliding mode [39], LMI-based SMC [40], Higher order SMC [41,42], Lyapunov-based adaptive SMC and Non singular terminal SMC in [43]. A widely and important criteria to ensure the stability of the SMC controllers is the Luyapunov theorem [44]. The proposed theorems ensure global stability of the nonlinear systems and control using the condition that the disturbances are bounded.
Considering the aforementioned literature review, this paper proposes adaptive robust formation and trajectory tracking of multiple UAVs using super-twisting sliding mode control method.The proposed controller compensates for the disturbances using adaptive control laws derived by Lyapunov function method.System stability is ensured using Lyapunov theorem. Furthermore, the formation flight between multiple UAVs are also controller using super twisting sliding mode control methods.
The rest of the paper is organized into four parts. First part is related to "system description and mathematical modeling". The objective of this part is the modelling of a single quad-rotor and transnational dynamics of multiple UAVs. The second part includes the trajectory and the formation controller formulation. Then the simulation results are presented and comparative analysis for different controllers is done. At last, conclusions are made. Following specific contributions are highlighted:

1.
In presence of external disturbances, robust formation and trajectory tracking of multiple UAVs is achieved using adaptive super twisting sliding mode control method.

2.
The adaptive laws are derived using Lyapunov theorem and implemented using projection operators. Figure 1 shows a quad-rotor UAV in earth's reference coordinates (X, Y, Z). Apart from inertial frame of reference, the body coordinates of the UAV are given as: (X B Y B Z B ). To derive the model, the following assumptions are made. Assumption 1. It is assumed that the UAVs are represented by a symmetrical rigid body configurations with masses m.

System Description and Mathematical Modeling
Assumption 2. The external disturbances affect the X and Y accelerations components of each UAV.

Assumption 3.
It is assumed that the disturbances are affecting the leader and followers UAV uniformly.Ẍ  Referring to Figure 2, the transnational dynamics of the UAVs are expressed as follows: where V Xi and V Yi represents the velocities in X and Y directions of the inertial frame. As shown in Figure 2, let the follower UAVs maintain d Xj and d Yj distances in X and Y planes respectively with respect to the leader UAV, so d Xj and d Yj are expressed as follows: The error in ψ dynamics as defined as follows: e ψ = ψ j − ψ L . By taking the first derivatives of Equations (10) and (11) with respect to time and combining the resultant expressions with Equations (7) and (8) yields the following expressions: where V Xj , V Yj , V XL and V YL represent the longitudinal and lateral velocities of the fol-lower1, follower2 and leader UAVs respectively. By defining errors in the longitudinal and lateral dynamics of Equations (12) and (13), the error state equation is represented as follows: Equation (14) are explained as follows: Also the terms G(χ) and F(χ) are expressed as follows: In Equation (17), e ψ is already defined, while c represents cos while s is a sin function. Also from Equation (15) UAVs are expressed as follows:

Trajectory and Formation Controllers Formulation
In this section, as a first step, the derivations of the attitude, altitude and position controllers are formulated for the leader UAV. As a second step, the formation controller is derived and based on it, new references are calculated for follower1 and follower2 UAVs. In the last step, the trajectory and attitude controllers of the leader UAV are generalized for follower UAVs. Before deriving the control schemes, the following assumptions are made: Assumption 4. It is assumed that the following condition is true for the uncertainty terms: represent the upper bound of the mentioned uncertainties.

Leader UAV Control Formulation
In this subsection, the attitude, altitude and position controllers are derived for leader UAV using adaptive super twisting sliding mode control method.

Attitude Control
Attitude controllers regulate the roll, yaw and pitch angles of the UAV. Let the reference Euler angle commands for the leader UAV are set as φ dL , θ dL, ψ dL , then the desired sliding manifold is chosen as follows: where S φ L represents the sliding surface for φ L loop, k 1 , k 2 are the design constants and the φ L loop error dynamics are expressed as follows, i.e., By taking the first time derivative of Equation (19), the following expression is obtained : Equations (4) and (20) are combined and expressed as follows: In Equation (21), the coefficients are defined as follows, i.e., a 1L = then the equivalent control law for φ L loop is derived as follows: Using super twisting algorithm, the switching control law is derived as follows: Referring to Equation (22) and (23), the total control action is the sum of equivalent and switching control parts, i.e., U 2L = U 2Leq + U 2Lsw . Similar procedure is adopted to derive the pitch and yaw controllers. The sliding surfaces for θ L and ψ L loops are defined as follows i.e., S θ L = k 3 e θ L + k 4ėθ L and S ψ L = k 5 e ψ L + k 6ėψ L , then the θ L and ψ L loops controllers are formulated as follows: From Equations (19)- (27), constant parameters k 1 , k 2 , k 3 , k 4 , k 5 , k 6 , k d1 , k d2 , k d3 , k d4 , k d5 , k d6 represent controllers and sliding surface gains. The coefficients are defined as follows: a 3L = I zL −I xL I yL , a 4L = J rL I yL , b 2L = lL I yL , a 5L = I xL −I yL I zL and b 3L = lL I xL , Moreover, S θ L and S ψ L represent the sliding surfaces for θ L and ψ L loops. The corresponding error dynamics for θ L and ψ L loops are expressed as follows, i.e., (4)-(6),satisfying assumptions 1-3, then under the proposed controllers of Equations (22)-(25), states of the attitude dynamics will converge to the origin in finite time [45].

Theorem 1. Consider the nonlinear system presented in Equations
Proof. Proof of Theorem 1 [45] The stability proof is only derived for φ L loop only. Similar procedures can be be adopted for the other two loops of attitude dynamics. Equation (23) is modified as follows: sgn(S φ L ). By combining the above terms with Equations (21) and (22),Ṡ φL is expressed as follows: Let the Lyapunov function for φ loop dynamics is chosen as follows: −τ 1 2 and then the Lyapunov function is expressed as follows: The time derivative of the Lyapunov function along (28) yields the following relation [45]: where the new matrices are represented as follows: Q φL = τ 1 2 2τ 2 + τ 2 1 −τ 1 −τ 1 1 and Applying the uncertainty bounds mentioned in Assumption 4, expression (29) is simplified as follows [45]: where (30) is negative definite only ifQ φL > 0. If the gains satisfy the following criteria Remark 1. The proof of finite time convergence property can be derived by using the procedures adopted in [45].

Altitude and Position Control
This section formulates the altitude and position control system for the leader UAV expressed in Equations (1)- (6). First the altitude control system is derived and then using the transformation matrix, the position controllers are formulated. With the desired altitude Z dL , the sliding manifold is written as follows: In Equation (31), k 7 and k 8 represent the design constant. The error dynamics are defined as follows: i.e., e Z L = Z L − Z d L ,ė Z L =Ż L −Ż Ld L . Taking the first time derivative of Equation (32) yields the following expression: Equation (32) and Equation (3) are combined and expressed as follows: Using super twisting sliding mode theory, the altitude controller is derived as follows: Here represents the virtual control law. The stability proof is derived based on the same concepts presented for φ loop. The robust terms of Equation (34) are modified as follows: U 1Lsw = −k d7 k 8 |S Z L | 0.5 sgn(S Z L ) + v Z L ; where the term v Z L is calculated from the following expression:v Z L = − k d8 k 8 sgn(S Z L ). By combining the above terms with Equation (33) and robust term of (34),Ṡ ZL is expressed as follows: Let the Lyapunov function for Z loop is chosen as follows: and then the Lyapunov function is expressed as follows: V Z L = η T ZL P ZL η ZL The time derivative of the Lyapunov function along (35) yields the following relation [45]: where the new matrices are represented as follows: Q ZL = τ 7 2 2τ 8 + τ 2 7 −τ 7 −τ 7 1 and Applying the uncertainty bounds given in Assumption 4, expression (36) is simplified as follows [45]: where matrixQ ZL = 2τ 8 + τ 2 7 − 4τ 8 τ 7 + τ 7 ∆ 3L −τ 7 + 2∆ 3L −τ 7 + 2∆ 3L 1 . Equation (37) is negative definite only ifQ ZL > 0. If the gains satisfy the following criteria τ 7 > 2∆ 3L , τ 8 > τ 7 5∆ 3L τ 7 +4∆ 2 3L 2(τ 7 −2∆ 3L ) , thenQ ZL > 0 andV Z L < 0. Now to derive the XY controllers , we assume the following: With these expressions, Equations (1) and (2) are re-written for leader UAV in the following form:Ẍ Let the sliding manifolds for the position loops of leader UAV are expressed as follows: In Equation (40) k 9 ,k 10 ,k 11 ,k 12 are the design constants and the error dynamics are expressed as follows: By taking the first time derivative of Equation 40, and combining it with Equations (38) and (39) one acquires the following expressions: From Equation (41), the virtual controllers U XL and U YL are expressed as follows: The stability proof is derived based on the same concepts presented for Z loop. The robust terms of Equation (42) are modified as follows: U XLsw = −k d9 k 10 |S X L | 0.5 sgn(S X L ) + v X L ; where the term v X L is calculated from the following expression:v X L = − k d10 k 10 sgn(S X L ) and sgn(S Y L ). By combining the above terms with Equations (41) and (42),Ṡ XL andṠ YL are expressed as follows: Let the Lyapunov function for X loop dynamics is chosen as follows: V X L = 2τ 10 |S X L | + 0.5v 2 X L + 0.5(τ 9 |S X L | 0.5 sgn(S X L ) − v X L ) 2 , where as for Y loop dynamics is the Lyapunov function is given as follows: Define new matrices as follows: P XL = 4τ 10 + τ 2 9 −τ 9 −τ 9 2 ; P YL = 4τ 12 + τ 2 11 −τ 11 −τ 11 2 and then the Lyapunov functions are expressed as follows: V X L = ζ 1 η T XL P XL η XL + D T XL D XL ; V Y L = ζ 2 η T YL P YL η YL + D T YL D YL The time derivative of the Lyapunov functions along (43) yields the following relation [45]: where: Q XL = τ 9 2 2τ 10 + τ 2 9 −τ 9 −τ 9 1 ; Q YL = τ 11 2 2τ 12 + τ 2 11 −τ 11 −τ 11 1 and q T XL = 2τ 10 + 1 2 τ 2 9 − 1 2 τ 9 ;q T YL = 2τ 12 + 1 2 τ 2 11 − 1 2 τ 11 . From Equation (44), since D XL and D YL are scalar quantities so D T XL = D XL and D T YL = D YL , then adaptive laws are derived as follows: Applying the uncertainty bounds given in Assumption 4, and by combining Equation (44) with Equation (45), the simplified expressions of (44) are given as follows [45]: Here the matrices are defined as follows: Q XL = 2τ 10 + τ 2 9 − 4τ 10 τ 9 + τ 9 E∆ 1L −τ 9 + 2E∆ 1L −τ 9 + 2E∆ 1L 1 and The expressionsV XL andV YL are negative definite only ifQ XL > 0 andQ yL > 0.
If the gains satisfy the following criteria τ 9 > 2E∆ 1L , τ 10 > τ 9 Here the terms E∆ 1L = D XL estimated − D XL ; E∆ 2L = D YL estimated − D YL represent estimation error of the adaptive loops.

Remark 2.
Discontinuous projection operator is used to implement the adaptive lawsḊ XL ,Ḋ YL . The projection operator is defined as follows: In Equation (45), ζ 1 and ζ 2 represent the adaptation gains. To generate reference trajectories for θ dL and φ dL , the virtual controllers U XL and U YL are expressed as follows: Multiplying Equation (48) by sin ψ and Equation (49) by cos ψ and the deference of the resultant equations yields the following expression: Equation (50) is simplified to get the reference command for φ loop of the leader UAV as follows: Multiplying Equation (48) by cos ψ and Equation (49) by sin ψ and the summation of the resultant equations yields the following expression: Squaring Equation (50) on both hand sides and equating sin 2 φ L = 1 − cos 2 φ L , the expression is written in terms of cos φ L and given as follows: Now form Equations (51) and (52), reference command for θ dL is expressed as follows:

Leader Followers Formation Control
Before discussing the trajectory controllers for the followers UAV, it is necessary to derive the formation controller which will generate the desired trajectory for the followers UAV. Let the following sliding surface is defined for formation controller: where τ represents gain matrix of the sliding surface. By taking time derivative of Equation (55) and combining it with Equation (14), one obtains the following expression: For formation control, the desired longitudinal and lateral velocities of the followers UAV are calculated as follows: For the stability proof, the same procedures as adopted for XYZ loops are applied here.

Followers UAV Control Formulation
In this section, we briefly explain the trajectory control of followers UAV.As mentioned above, the reference position trajectories are generated using Equation (18), and governed by the formation controller of Equations (57) and (58). Thus, by defining the attitude, altitude and position errors for the followers UAV, the rest of the analysis used for the derivation of the subject controllers is the same as leader UAV. For simplicity, here the final control laws are included: let the attitude sliding manifolds are defined as follows: where j = [F 1 , F 2 ] and F 1 and F 2 represent follower1 and follower2 UAVs respectively. Also k 1j ,k 2j ,k 3j ,k 4j ,k 5j ,k 6j represent the constants of sliding surfaces for followers UAV. The respective errors are defined as follows: e φ j = φ j − φ dj ; e θ j = θ j − θ dj ; e ψ j = ψ j − ψ dj . Similarly position and altitude sliding manifolds for follower1 and 2 are given as follows: where k 7j ,k 8j ,k 9j ,k 10j ,k 11j ,k 12j represent the constants of sliding surfaces for followers UAV. The respective errors are defined as follows: e Z j = Z j − Z dj ; e X j = X j − X dj ; e Y j = Y j − Y dj . Now following the same procedures, the attitude, altitude and position controllers for followers UAV are formulated as follows: The reference commands for θ dj and φ dj are derived using the same procedures given in Equations (48)-(54).

Results and Discussion
In this section, the proposed ASTSMC controller is tested numerically for the system of multiple quad-rotors shown in Figure 2. The parameters of the leader and followers UAVs are identical and given in Table 1. The control parameters of leader and follower UAVs are given in Tables 2 and 3. Since the UAVs are identical, leader and followers UAV use the same control parameters.The parameters for formation control loops are chosen as follows: For the leader UAV, the reference position and altitude commands are set as follows: X L = sin t, Y L = cos t and Z L = t. Figure 3 shows the applied acceleration disturbance on X and Y dynamics of the leader and followers UAV. it is assumed that same type of disturbance acceleration is applied for all UAVs. Furthermore the disturbance acceleration has no effect on the Z dynamics of UAV. Moreover, the following parametric uncertainties are applied: a 1L = 2.5a 1L ;a 1j = 2.5a 1j ;a 2L = 2.5a 2L ;a 2j = 2.5a 2j ;a 3L = 2.5a 3L ;a 3j = 2.5a 3j ;a 4L = 2.5a 4L ;a 4j = 2.5a 4j ;a 5L = 2.5a 5L ;a 5j = 2.5a 5j ;b 1L = 1.75b 1L ;b 1j = 1.75b 1j ;b 2L = 1.75b 2L ;b 2j = 1.75b 2j ;b 3L = 1.75b 3L ;b 3j = 1.75b 3j .

Symbol
Value Unit 6.5 × 10 −5 kg·m 2 I xL = I xF1 = I xF2 7.5 × 10 −3 Ns 2 rad −1 I yL = I yF1 = I yF2 7.5 × 10 −3 Ns 2 rad −1 I zL = I zF1 = I zF2 1.3 × 10 −2 Ns 2 rad −1   Figure 4 shows the trajectory tracking simulations of leader followers UAV in presence of applied disturbance of Figure 3. From the presented results, it is concluded that in presence of disturbances, ASTSMC controllers ensure robust behaviour, while the fixed gain STSMC controllers exhibit steady state errors in the X and Y tracking responses of the leader and followers UAVs. To have a clear picture of the trajectory deviations under wind disturbance, Figures 5-8 show the trajectory tracking comparisons in XY plane for leader, follower1 and follower2 UAVs respectively. From the presented results, it is concluded that minimum deviations are observed in the trajectory tracking for all UAVs with ASTSMC controllers, while with fixed gain STSMC controllers all UAVs show significant drift from the reference trajectories in XY plane. Figure 8 shows the combined trajectories of leader followers UAVs with ASTSMC and fixed gain STSMC controllers in XY plane. From the presented results, it is obvious that the proposed ASTSMC controllers ensure robust formation control between the leader and followers UAVs, while with fixed gain STSMC controllers, all UAVs show drift in their trajectories.    To compare the trajectory tracking performance of the leader followers UAV quantitatively, the X and Y trajectories are individually plotted against time and the results are presented in Figures 9 and 10 for the leader UAV. From Figure 9, and at time t = 45 s, the e XL tracking error is 0.1 m with fixed gain STSMC controller while with ASTSMC controller, the measured error e XL is 0.05. ASTSMC ensures lowest error due the adaptive disturbance compensation term D XL and from the presented results of Figure 9, it is obvious that at time t = 45 s, the adaptive term D XL adds appropriate compensation to cancel the disturbance and it switches from 50 to −100. Similarly from Figure 10, e YL is measured 0.28 m and 0.05 m with fixed gain STSMC and ASTSMC controllers respectively. ASTSMC controller offers lowest error due to the adaptive estimator term D YL . From the presented results of Figure 10, it is obvious that at time t = 45 s, the adaptive term D YL adds appropriate compensation to cancel the disturbance and it switches from 0 to −150. Similarly for followers UAV, the X and Y tracking responses are plotted against simulation time and the results are shown in Figures 11-14. From the presented results and at time t = 45 s, the measured error signals with fixed gain STSMC controller are as follows: e X F1 = 0.2 m, e Y F1 = 0.3 m, e X F2 = 0.2 m, e Y F1 = 0.5 m , while with ASTSMC controllers , the errors are measured as follows: e X F1 = 0.1 m, e Y F1 = 0.1 m, e X F2 = 0.05 m, e Y F1 = 0.05 m. From the presented results of Figures 11-14, it is obvious that at time t = 45 s, the adaptive terms D X F1 ,D Y F1 , D X F2 ,D Y F2 add appropriate compensation to cancel the disturbances.   Figures 15 and 16 show θ and φ tracking responses for leader and follower UAVs with both fixed gain STSMC and ASTSMC controllers respectively. From the presented results and at time t = 45 s, it is evident that the proposed ASTSMC controllers generate appropriate reference commands for both θ and φ loops of leader and follower UAVs. To have a better understanding of the above claim, Figure 17 shows the difference of the generated reference θ and φ commands with ASTSMC and fixed gain STSMC controllers. From the presented results , it is obvious that at time t = 45 s, the proposed ASTSMC controllers generate appropriate reference commands for both θ and φ loops of leader and follower UAVs.       Figures 18 and 19 show Z and ψ loops tracking responses for leader and follower UAVs. Since no disturbances are applied on both these loops, the tracking responses under fixed gain STSMC and ASTSMC controllers are comparable. Finally, Figure 20 shows the robustness of the formation controllers for tracking the respective reference commands i.e., the distance between the leader and the followers in X , Y plane. From the presented results it is obvious that apart from the transient error, the formation controllers accurately maintain the desired distance between the leader-follower1 and leader-follower2 UAVs.     Figure 21 shows the simulation results of the control inputs using the proposed control schemes. Since the attitude loops are not adaptive and gains are fixed so the control signals chatters, however the control inputs are feasible for practical implementations and well bounded. Moreover, the virtual XY control outputs of the proposed control schemes offer vert low chattering (Figures 15 and 16) In future work, gains of the proposed control schemes will be tuned online to overcome the chattering phenomena.

Conclusions
This paper proposes adaptive super twisting sliding mode trajectory and formation controllers for multiple UAVs flying the leader follower configuration. Acceleration type disturbances and parametric uncertainties are applied to = X and Y dynamics of leader and follower UAVs. The formation control of UAVs is tested with the proposed ASTSMC and fixed gain STSMC controllers. The robust performance of the proposed control is verified from the following measured errors of the leader and follower UAVs. For leader UAV, e XL = 0.05 m, e YL = 0.05 m with ASTSMC control while with the fixed gain STSMC controller, the measured errors are as follows: e XL = 0.1 m, e YL = 0.25 m. Similarly for follower UAV,e X F1 = 0.09 m, e Y F1 = 0.05 m, e X F2 = 0.09 m, e Y F1 = 0.04 m, with ASTSMC control while with the fixed gain STSMC controller, the measured errors are as follows: e X F1 = 0.18 m, e Y F1 = 0.25 m, e X F2 = 0.2 m, e Y F1 = 0.14 m. Moreover, the settling time of XY states after the occurrence of disturbances is faster as compared to fixed gain STSMC control methods. From the quantitative comparison here; it is concluded that the proposed ASTSMC controllers show enhanced robust behaviour to the acceleration type disturbances and parametric uncertainties of the system.