Next Article in Journal
Multi-Monostatic Interferometric Radar with Radar Link for Bridge Monitoring
Next Article in Special Issue
YOLO MDE: Object Detection with Monocular Depth Estimation
Previous Article in Journal
Development of a 0.15 μm GaAs pHEMT Process Design Kit for Low-Noise Applications
Previous Article in Special Issue
Autonomous Trajectory Planning and Control of Anti-Radiation Loitering Munitions under Uncertain Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Angular Rate Constrained Sliding Mode Control of UAVs for Path Following

Department of Aerospace Engineering, Chosun University, Gwangju 61452, Korea
*
Author to whom correspondence should be addressed.
Electronics 2021, 10(22), 2776; https://doi.org/10.3390/electronics10222776
Submission received: 4 October 2021 / Revised: 8 November 2021 / Accepted: 11 November 2021 / Published: 12 November 2021

Abstract

:
In this work, a sliding-mode-based attitude controller constrained with the angular rate for unmanned aerial vehicles (UAVs) is addressed to withstand conditions below the allowable maximum angular velocity of UAVs in order to avoid the possibility of structural failure and to operate UAVs safely. The sliding mode controller suggested in this work defines a new sliding surface, inherently having two equilibrium points. These equilibrium points are carefully inspected, and the stability of the system controlled by means of the proposed approach is also analyzed using Lyapunov stability theory. To highlight the angular-rate constrained attitude control technique, a three-dimensional path is constructed using the Dubins path technique, and three-axis attitude commands for UAV are also generated by augmenting the line-of-sight algorithm. Compared with conventional sliding mode control measures, the excellent performance of the suggested control algorithm has been demonstrated by conducting numerical simulations.

1. Introduction

Various attitude control techniques designed to improve the maneuvering capability of emerging high agile unmanned aerial vehicles (UAVs) have been studied in the past several decades [1]. To meet the requirement for accomplishing complex missions that needed to be performed in a short time, modern UAVs demand increasingly high levels of maneuverability [2]. On the other hands, it is also known that the excessive maneuvering capability of UAVs increases the risk of structural damage or cracks due to induced gravitational forces and moments. Thus, resolving these opposite characteristics properly is very important in order to operate UAVs effectively and safely, which was the major motivation for this study.
Let us review the attitude control techniques applied to UAVs. The most commonly used control technique is the classical proportional-integral-derivative (PID) algorithm, which provides very reliable control performance and stability, despite its simple structure [3]. Nevertheless, due to its inherent limitations as a linear control method, this control law is not appropriate to apply to highly nonlinear systems, and it is quite vulnerable in its control of the system to external disturbances or internal uncertainty [4]. As an alternative, various nonlinear control theories, such as sliding mode control (SMC) [5], backstepping control [6], nonlinear model predictive control [7], and adaptive control [8], have been studied in order to improve the maneuverability of UAVs.
One of the options to increase the maneuverability of modern UAVs is the backstepping algorithm, which is a nonlinear technique that adjusts other state variables by sequentially setting each state as a virtual control input. The drawback of this approach is that it largely relies on having a quite accurate system model. To solve this uncertainty problem in the nonlinear platform, the adaptive backstepping control was augmented by approximating the uncertainty in the system model using a fuzzy-based approach [9]. Meanwhile, to control unsteady and highly nonlinear UAVs, the nonlinear model predictive control technique was introduced [10]. In this approach, the control input is obtained by selecting the current output among the several discretized control outputs, evaluated from the present time to a few seconds by means of the optimization technique. The advantage of this method is that it allows us to evaluate the optimal control input for a given cost function by taking into account the various constraints. However, this powerful merit ironically leads to a disadvantage in that it requires very high computational resources to evaluate the optimal value. That is, the larger the size of the constraint to be considered in the optimization process, the greater the computation burden required [11]. In an effort to reduce the amount of computational resources, some studies have been conducted. As an effective method of support for model predictive control, the neural network technique was introduced [12]. Note that this adaptive-like control involves a technique in which the gain value of the controller changes according to the state of the system, and has a disadvantage in that the control performance depends on the initial value of the gain. Therefore, learning-based control algorithms including the adaptive control technique have an implicit burden—the question of how to select the stable initial values or parameters when this technique is chosen as the primary option to control the actual platform.
The sliding mode control (SMC) approach for UAV attitude control applications has useful characteristics, including a guaranteed control response time. This approach is a control method that creates a sliding surface with state variables to control a nonlinear system and directs the state variables to the control target by sliding on the surface using a reaching law, which guarantees robustness [13]. This strong property leads to the biggest advantage of SMC. On the other hand, it is also known that there is the drawback of a high-frequency oscillation phenomenon, that is, a chattering problem, in the process of controlling and compensating for disturbance and uncertainty. To solve this problem, several studies have been carried out to find an effective reaching law [14,15].
In spite of the efforts to improve the attitude control capability of UAVs in the presence of external disturbances and internal uncertainty, less attention has been paid to structural constraints or dynamic limitations for safe operations in controller design. Merely increasing the maneuvering capability of UAVs does not cover the safe operation of UAVs to conduct complex missions. In this work, a nonlinear control law for driving below the allowable maximum angular velocity is considered in order to realizethe possibility of structural failures and safety uncertainties. To ensure the stable dynamic ranges and the capability of attitude control, the sliding mode control (SMC) technique is taken in account. Since this not only requires fewer computational resources, but also has strong properties against model uncertainty and perturbation, SMC could be one of the reliable options used in the implementation of actual UAVs. The newly designed control law in this work is called constrained SMC (CSMC). To guarantee the stable dynamic range by limiting the angular rate, a sliding surface is proposed. It is revealed that the suggested sliding surface has two equilibrium points in order to meet the exceptional requirement examined in this work. Thus, these equilibrium points are carefully observed, and the stability at each point is also investigated. Finally, the attitude control law, using the CSMC approach, is applied to UAVs and the stability is proven using Lyapunov stability theory. Next, a 3-dimensional guidance algorithm for path following of UAVs is additionally addressed to highlight the control law. To provide guidelines for the operation of UAVs under the maximum allowable angular rate, the concept of Dubins curves has been considered in order to generate 3-dimensional smooth guidance curves. Then, the path-following algorithm is augmented to CSMC, as proposed in this paper. The excellent performance of the proposed approach, in cooperation with the path-following technique, for safe UAV operations is demonstrated by integrated simulation studies involving external disturbances and internal uncertainties.
This paper is organized as follows. First, the equations of motion of a fixed-wing UAV are introduced in Section 2. Then, the conventional SMC is briefly reviewed, and the SMC-based attitude control method limiting the angular rate of UAVs is addressed in Section 3. The uniquely designed sliding surface is also described. Next, the stability analysis of the closed-loop system with the proposed control law is given for two equilibrium points. Then, the 3-dimensional path generation algorithm for UAV attitude command generation is followed in Section 4 using the Dubins curve concept. Finally, the proposed control scheme is demonstrated using numerical simulation studies. In this simulation, the CSMC technique is compared with the conventional SMC. Furthermore, to highlight the proposed control method, several harsh scenarios are provided to ensure that the UAV with the generated attitude command follows a given guidance route smoothly and safely.

2. Mathematical Model of UAVs

Various types of UAVs, for example, fixed-wing, rotary-wing, and multi-copter, have been designed and studied during the past several decades. Note that a fixed-wing UAV is employed in this paper to design the control law. Now, let us start by reviewing the coordinate system presented in Figure 1. The position vector of the illustrated UAV is defined by ξ = [ X Y Z ] T , and the attitude is defined by the Euler angle η = [ ϕ θ ψ ] T R 3 , which can also be converted to quaternion parameters [ q q 4 ] T R 4 by the attitude transformation. Note that most of the literature for the attitude representation of fixed wing UAVs has been focused on Euler angles rather than quaternions, since the Euler angle is much more intuitive for controller design. The attitude control law considered in this work, however, employs the quaternion representation due to its simple mathematical structure. Thus, the governing equations of motion of the fixed-wing UAV are given by [16].
(1) ξ ˙ = R ν (2) ν ˙ = R T 0 0 g ω × ν + 1 m F + 1 m T 0 0 (3) q ˙ = 1 2 ( q × + q 4 I 3 × 3 ) ω (4) q 4 ˙ = 1 2 q T ω (5) J ω ˙ = ω × J ω + M
where ν = [ u v w ] T represents the linear velocity vector, and ω is the angular velocity vector of the UAV in the body frame. F R 3 and M R 3 denote the aerodynamic force and moment in the body frame, respectively. m , g , and T represent the mass of the UAV, the gravitational constant, and the thrust of the UAV, respectively. Furthermore, I 3 × 3 R 3 × 3 represents the 3 × 3 identity matrix. The inertia matrix J is the symmetric positive matrix and is written in the body frame as
J = J x x 0 J x z 0 J y y 0 J z x 0 J z z
The attitude transformation matrix from the body frame to the inertia frame can be expressed by the quaternion as
R = ( 2 q 4 2 1 ) I 3 × 3 + 2 q q T 2 q 4 q ×
Note the operator × on any vector x = [ x 1 , x 2 , x 3 ] T represents a skew-symmetric matrix defined as
x × = 0 x 3 x 2 x 3 0 x 1 x 2 x 1 0
The aerodynamic force and moment F and M are written as [16]
(9) F = C x q ¯ S C z q ¯ S C z q ¯ S T (10) M = q ¯ S b C l q ¯ S c ¯ C m q ¯ S b C n T
where the aerodynamic coefficients can be given by
(11) C L = C L 0 + C L α + C L δ e δ e + C L α ˙ α ˙ + C L q q c ¯ / ( 2 V a ) (12) C D = C D 0 + C D δ e δ e + C D δ r δ r + C L + C L m i n / ( π e A R ) (13) C x = C D c o s α + C L s i n α (14) C y = C y β β + C y δ r δ r + C y p p + C y r r (15) C z = C D s i n α C L c o s α (16) C l = C l β β + C l δ a δ a + C l δ r δ r + C l p p + C l r r b / ( 2 V a ) (17) C m = C m 0 + C m α α + C m δ e δ e + C m α ˙ α ˙ + C m q q c ¯ / ( 2 V a ) (18) C n = C n β β + C n δ a δ a + C n δ r δ r + C n p p + C n r
Furthermore, b, S and c ¯ denote the wingspan, the wing area, and the mean aerodynamic chord (MAC) of the UAV, respectively. Furthermore, q ¯ represents dynamic pressure, e is the Oswald efficiency number, and A R is the wing aspect ratio. Lastly, δ a , δ e and δ r are the control input for the aileron, elevator, and rudder, respectively. The angle of attack and the sideslip angle are defined respectively as
α = arctan w u
β = arcsin v V a
Note that V a = ν T ν is the airspeed of the UAV.
To facilitate the derivation of an attitude controller design, the aerodynamic moment in Equation (10) can be linearly decoupled with two parts. Thus, the aerodynamic moment can be rewritten in a simple vector form with Equations (16)–(18) as
M = f + Λ u
where
(22) f = q ¯ S b 0 0 0 q ¯ S c ¯ 0 0 0 q ¯ S b C l β β + C l p p + C l r r b / ( 2 V a ) C m 0 + C m α α + C m α ˙ α ˙ + C m q q c ¯ / ( 2 V a ) C n β β + C n p p + C n r (23) Λ = q ¯ S b 0 0 0 q ¯ S c ¯ 0 0 0 q ¯ S b C l δ a 0 C l δ r 0 C m δ e 0 C n δ a 0 C n δ r (24) u = δ a δ e δ r T

3. Angular-Rate-Constrained Nonlinear Attitude Control

In this section, a new nonlinear attitude control law considering the maximum angular rate for fixed-wing UAVs is presented. To employ this unique approach, the sliding mode control (SMC) technology concept is utilized.

3.1. Sliding Mode Control

Let us first consider the baseline SMC. The first step in this technique is to reach the designed sliding surface. Then, the second step is the sliding stage, in which the sliding variable reaching the sliding surface goes to zero, which means that state errors converge to zero. For attitude control using the sliding control approach, the representative sliding surface s = [ s 1 s 2 s 3 ] T is chosen as
s = ω + a q e
and a simple reaching law s ˙ is written as
s ˙ = k s g n ( s )
where q e is the vector part of quaternion error defined as q e = q q c 1 and q c is the quaternion command, and the operator ⊗ refers to the quaternion multiplication. Furthermore, a and k are positive design parameters. Note that s g n ( ) is the signum function, defined as
s g n ( s ) = s i g n ( s 1 ) s i g n ( s 2 ) s i g n ( s 3 ) T
To compute the control output from the sliding surface, the time derivative of sliding surface in Equation (25) is expressed as
s ˙ = ω ˙ + a 1 2 ( q e × + q e , 4 I 3 × 3 ) ω
Inserting the angular acceleration, ω ˙ , in Equation (5) into the above equation leads to
s ˙ = J 1 ( ω × J ω + f + Λ u ) + a 1 2 ( q e × + q e , 4 I 3 × 3 ) ω
In terms of u from the above equation, the control input is expressed as
u = Λ 1 ω × J ω + f + a J 1 2 ( q e × + q e , 4 I 3 × 3 ) ω + k J s g n ( s )
It is known that the discontinuity in the reaching law introduces the chattering problem. To release the burden of chattering, the alternative reaching law is given by
s ˙ = k 1 s k 2 s ϵ s g n ( s )
where ϵ is a design parameter ranging from 0 to 1, and | s | ϵ R 3 × 3 is a matrix function defined as
s ϵ = d i a g ( | s 1 | ϵ , | s 2 | ϵ , | s 3 | ϵ )
and d i a g ( ) is the 3 × 3 diagonal matrix in this case. By inserting the above reaching law into Equation (30) to mitigate the chattering problem, the control input is rewritten as
u = Λ 1 ω × J ω + f + a J 1 2 ( q e × + q e , 4 I 3 × 3 ) ω + J k 1 s + k 2 s ϵ s g n ( s )
Note that the final form of the control input is the attitude sliding mode control law for UAVs, overcoming the inherently introduced chattering problem.
Lemma 1.
Once the sliding manifold s ( t ) = 0 is satisfied with properly chosen parameters, then the desired attitude maneuver can be achieved, i.e., the variable q e and ω will converge to zero. That is,
lim t q e ( t ) = 0
lim t ω ( t ) = 0
Proof. 
Assume that the sliding surface in Equation (25) is zero, and s = 0 , then it can be expressed as
ω = a q e
Substituting ω into Equation (4) and setting q q e introduces
q ˙ e , 4 = a 1 2 q e T q e
Due to the norm constraint of the quaternion given by q e T q e = 1 q e , 4 2 , the right-hand side of the equation is rewritten as
q ˙ e , 4 = a 1 2 ( 1 q e , 4 2 )
The closed-form solution of the differential equation for a given time is q e , 4 ( t ) = tanh a 2 t . As a sufficient time has elapsed, it can be seen that q e , 4 converges to 1:
lim t q e , 4 ( t ) = 1
With q 4 converging to one, this means that q e converges to the zero vector due to the norm constraint of the quaternion after a sufficient time has elapsed. Consequently, the sliding surface s approaches zero, which means that q e converges to zero independently. Furthermore, ω also converges to zero according to Equation (36). Thus, Lemma 1 is proven. □

3.2. Angular-Rate-Constrained Sliding Mode Control

In this subsection, a modified control law based on SMC is introduced by defining a sliding surface proposed in this work. Let us first assume that the fixed-wing UAV has limited maneuverability to prevent structural failure or cracks or to operate various missions safely. Without the loss of generality, the angular rate is directly linked with the magnitude of the centrifugal force according to the given airspeed of the UAV. Thus, it is natural that the maneuverability constraint can be converted to the angular rate limitation of the UAV. That is,
ω i ω m
where ω i is the angular rate of UAV for each body axis, and ω m is the allowable maximum angular speed of UAV to meet the maneuverability constraints. Now, to satisfy the angular rate constraint, the new sliding surface suggested in this work is written as
s = ω + a s a t L ( q e )
where the saturation vector function is defined as
s a t L ( q e ) = s a t L ( q e , 1 ) s a t L ( q e , 2 ) s a t L ( q e , 3 ) T
and the saturation scalar function can be given by
(43) s a t L ( q e , i ) = m i n ( L , q e , i ) s i g n ( q e , i ) (44) L = ω m a
Furthermore, L is the limiting parameter, obtained by dividing the maximum angular rate ω m by the design parameter a. The function m i n ( ) compares two components and selects a smaller value so that the saturation function selects a smaller value via the comparison between the component of the quaternion absolute error q e , i and the limiter L. From the inherent representation of the new sliding surface, there are two equilibrium points, setting s = 0 . That is, ω = q e = 0 and ω = a s a t L ( q e ) . The first equilibrium point is related to the attitude control purpose, and the second equilibrium point is deeply related to the angular rate limitation by forcing the UAV not to exceed the given limitation. Let us consider for the second equilibrium point that the quaternion absolute error q e , i is larger than the limiter L. Then, s a t L ( q e ) will be L according to Equation (43). Since L = ω m a , the saturation function can be the allowable maximum angular rate, that is, a s a t L ( q e , i ) = ω m s i g n ( q e , i ) . Thus, the second equilibrium point is related with the allowable maximum angular rate such that ω i = ω m s i g n ( q e , i ) .
For the angular rate constrained control law design, the time derivative of the sliding surface in Equation (41) is given by
s ˙ = ω ˙ + a D 1 2 ( q e × + q e , 4 I 3 × 3 ) ω
where D is the diagonal matrix with the element D i defined as
(46) D = d i a g ( D 1 , D 2 , D 3 ) (47) D i = 1 , if L q e , i L 0 , otherwise
Note that D i is differentiation of the scalar s a t L ( ) function, and it becomes 0 or 1 according to the algebraic comparison of L and q e , i . Thus, by substituting Equations (5) and (31) into (45), the constrained sliding mode control (CSMC) input can be expressed as
u = Λ 1 ω × J ω + f + a J D 1 2 ( q e × + q e , 4 I 3 × 3 ) ω + J k 1 s + k 2 s ϵ s g n ( s )
It is noted that the control input induced from the suggested sliding surface is the angular-rate constrained attitude control law for fixed-wing UAVs based on the sliding mode control. As seen in Equation (48), the control law is dedicated to the magnitude of attitude errors. If the attitude error q e , i for each axis is larger than the reference value of L, the term of a J 1 2 q × + q 4 I 3 × 3 ω is eliminated to improve the maneuverability. Otherwise, the term is activated, and the control law in Equation (48) supports both ω and q e to approach zero. In other words, it can be interpreted that the control law plays two important roles, since there are two equilibrium points. The first equilibrium point of the sliding surface, ω = a s a t L ( q e ) , is related to the case of D i = 0 . In this case, the control law allows us to reach the allowable maximum angular speed of the UAV in order to improve the maneuverability. This strategy is made possible by approaching the first equilibrium point. Next, for the second equilibrium point, ω = q e = 0 , connected with the case of D i = 1 , the control input causes the attitude error and angular velocity to converge to zero by controlling the sliding surface to reach the second equilibrium point. This property can be regarded as the unique characteristic of the constrained sliding mode control approach suggested in this paper.

3.3. Stability Analysis

Before verifying the stability of the control law, let us first inspect the overall process of the control law to ensure that the sliding surface works appropriately to achieve the control objective. As mentioned previously, there are two equilibrium points in the sliding surface. This means that the structure of the proposed control law can vary depending on the amount of attitude errors and the allowable maximum angular rate.
Suppose that a large attitude command is generated for the UAV in a steady state so that a large amount of attitude errors are applied to CSMC. It is matched with the case of D i = 0 . That is, the first step to achieve the goal is to adjust the angular rate ω i to converge to the allowable maximum angular rate ω m . This can done by sliding at the first equilibrium point. It is noted that approaching the first equilibrium point is deeply related with the maneuverability, since the angular rate increases quickly to the allowable maximum angular rate of the UAV. Consequently, due to the consistent angular rate applied to the UAV at the first equilibrium point, the magnitude of the attitude errors gradually decreases below the given reference L. This is the reason the control law can naturally move on to the second step, D i = 1 . As seen previously, the second step is to approach the second equilibrium point, ω = q e = 0 which is a step for controlling both the quaternion error q e , i and the angular rate ω i , causing them to converge to zero.
Let us carefully investigate the sliding surface in detail. According to the attitude error, the sliding surface defined in Equation (41) can be divided by two forms given by
(49) s i = ω i + ω m s i g n ( q e , i ) (50) s i = ω i + a q e , i
Equation (49) is the sliding surface in the case of large attitude errors | q e , i | enough to be larger than L, so that the point satisfying ω i = ω m s i g n ( q e , i ) is the equilibrium point. Furthermore, it is obvious that the variable of D i in Equation (47) becomes 0. Let us substitute the control law in Equation (48) into the governing equations of motion of UAVs in Equation (5). Then, it is found that the relationship between ω ˙ and s is given by
ω ˙ = k 1 s k 2 s ϵ s g n ( s )
From the above equation, it can be seen that the angular acceleration ω i ˙ is opposite to the sign of the sliding surface s i in Equation (49).
Assuming that the UAV is initially in a steady state, the angular rate of the UAV is zero. It holds s i = s g n ( q e , i ) ω m from Equation (49) so that the sign of s i is identical with the quaternion error, q e , i . This property is still valid, whereas s i is approaching to zero. In other words, It is valid before the angular rate of UAV, ω i , is identical with the allowable maximum angular speed, ω m . From Equation (51), it is obvious due to the fact that the signs of the quaternion error and the angular acceleration are opposite. That is, the angular rate increases while the quaternion error decreases. Consequently, ω i becomes ω m s i g n ( q e , i ) so that the sliding surface goes to zero, s i = 0 . Then the angular acceleration in Equation (51) becomes zero. Finally, it is noted that the angular rate reaching the allowable maximum angular rate remains unchanged, whereas the quaternion error decreases continuously to be smaller than L. Once the quaternion error becomes small after a sufficient time has elapsed, then the variable D i is switched to one. This is when the sliding surface of Equation (49) turns into Equation (50). It can be also seen that the sliding surface is the same as that of the conventional SMC in Equation (25). Therefore, the point where q e , i and ω i become zero is the equilibrium point.
Now, let us investigate the stability of the closed-loop attitude control system utilizing CSMC to ensure that the motion of the sliding surfaces work correctly. Stability analysis is required for each sliding surface. For the stability of the closed-loop system, the representative Lyapunov candidate by the first sliding surface in Equation (49), is defined as
V L = 1 2 s T s
Inserting Equation (45) into the time derivative of the Lyapunov candidate leads to
V ˙ L = s T s ˙ = s T ω ˙ + a D 1 2 ( q e × + q e , 4 I 3 × 3 ) ω
Then, let us substitute Equation (5) into the above equation, and replace the control input with Equation (48). Then, the time derivative of the Lyapunov candidate is rewritten as
V ˙ L = s T J 1 ( ω × J ω + f + Λ u ) = s T k 1 s k 2 s ϵ s g n ( s )
Note that D is zero in this case.
Moreover, the second term of the right-hand side of the above equation is always positive. That is,
k 2 s T s ϵ s g n ( s ) = k 2 i = 1 3 | s i | | s i | ϵ > 0
Thus, the time derivative of the Lyapunov candidate is given by
V ˙ L = k 1 s 2 k 2 i = 1 3 | s i | | s i | ϵ < 0
where s R denotes the two-norm of s . Since the time derivative of the Lyapunov candidate is always negative, the closed-loop system is asymptotically stable. This means that for a given initial condition of ω and q e , the sliding surface, s i , in Equation (49) will converge to the first equilibrium point, ω i = ω m s i g n ( q e , i ) .
Once again, for the closed-loop system stability by the second equilibrium point, the identical Lyapunov candidate by the sliding surface in Equation (50) is also defined as
V L = 1 2 s T s
By proceeding identically with the previous case, the time derivative of the Lyapunov candidate is also written as
V ˙ L = s T J 1 ( ω × J ω + f + Λ u ) + a D 1 2 ( q e × + q e , 4 I 3 × 3 ) ω = s T k 1 s k 2 s ϵ s g n ( s )
Note that the variable D does not disappear in this case. However, applying the control input in Equation (48), the remaining procedure is identical with that of the previous case. Since the closed-loop system is asymptotically stable for the given condition of L q e , i L , the sliding surface, s i , in Equation (50) will converge to the second equilibrium point, that is, ω i = q e , i = 0 , which is proven by Lemma 1.

3.4. Summary

For the attitude control of fixed-wing UAVs that are able to be operated within limited angular rates, the sliding mode control investigated in this section, similar to variable structure control technologies, is summarized as follows. This technique consists of two control laws separated by the amount of the attitude errors induced by the attitude commands and the allowable maximum angular rate of the UAV. If the attitude errors are larger than the limiter, for example, | q e , i | > L , then the related sliding surface and control law are given respectively by
(59) s = ω + ω m s g n ( q e ) (60) u = Λ 1 ω × J ω + f + J k 1 s + k 2 s ϵ s g n ( s )
otherwise, the relevant sliding surface and the control law are expressed respectively as
(61) s = ω + a q e (62) u = Λ 1 ω × J ω + f + a J 1 2 ( q e × + q e , 4 I 3 × 3 ) ω + J k 1 s + k 2 s ϵ s g n ( s )

4. 3D Path-Following Technique

In this section, a three-dimensional guidance algorithm for the path following of waypoints is additionally employed to ensure that the control law in Equation (48) works effectively. To provide the guidelines of the angular rate for a given UAV to be operated safely within the allowable forces and moment, the concept of the Dubins curve is introduced for the generation of three-dimensional smooth guidance curves. Thus, let us describe in this section the entire process used to generate the 3-dimensional Dubins path and the tracking algorithm with the turning radius of the UAV, determined for safe operations.

4.1. Dubins Path Generation

It is known that the Dubins path consists of only a few circles and lines so that it can be mathematically simplified for the generation of waypoint guidance for UAVs. The illustrative Dubins path indicating these two points is illustrated in Figure 2. The process of geometrically generating the Dubins path, for example, a CSC (Circular, Straight, Circular line) type in the 3D space is reviewed briefly in [17]. The minimum turning radius of the UAV, r m , can be first determined by the maximum allowable angular rate, ω m , and the given airspeed, V a . Assume that there is a UAV at a position ξ 1 heading to the unit vector, v ¯ ς , 1 . There is an another waypoint ξ 2 heading to the unit vector v ¯ ς , 2 that the UAV must pass. In Figure 2 the straight line, x ς , connecting two points of p 1 and p 2 is the Dubins line to generate the guidance route. The unit heading vector at a position ξ , which can be evaluated by means of the velocity vector of the UAV, is given by
ν ¯ = ξ ˙ ξ ˙
To define the three orthogonal unit vectors at the point, the vector perpendicular to the plane, consisting of the Dubins line and velocity vector, can be given by
u ς = x ς × ξ ˙
and the orthogonal unit vector is also written as
u ¯ ς = u ς u ς
Then, the third orthogonal unit vector, namely, the radius unit vector which is oriented toward the center of the circle from the position is evaluated as
w ¯ ς = u ¯ ς × ν ¯ ς
Then, using the unit vector and a given way point vector, the center point of the circle can be computed as
o ς = ξ + r m w ¯ ς
Next, the unit radius vector, y ¯ , leading from the point p to the center of the circle is given by
y ¯ ς = x ¯ ς × u ¯ ς
where x ¯ ς is the unit vector given by x ς x ς . The point p at which the tangent meets in the circle is also computed with respect to the center of the circle as
p = o ς r m y ¯ ς = ξ + r m w ¯ ς r m y ¯ ς
Consequently, since there are two waypoints illustrated in Figure 2, the straight line x ς connecting two points p 1 and p 2 can be evaluated as
x ς = p 2 p 1 = ( ξ 2 + r m w ¯ ς , 2 r m y ¯ ς , 2 ) ( ξ 1 + r m w ¯ ς , 1 r m y ¯ ς , 1 )
Furthermore, the angles denoted θ ς , 1 and θ ς , 2 shown in Figure 2 are computed using the inner product of the two unit vectors as, respectively,
cos θ ς , 1 = v ¯ ς , 1 · x ¯ ς
cos θ ς , 2 = v ¯ ς , 2 · x ¯ ς
From the geometrical relationships, the straight line vector is also given by
x ς = ( ξ 2 ξ 1 ) r m ( x ¯ ς + v ¯ ς , 1 ) tan θ ς , 1 2 r m ( x ¯ ς + v ¯ ς , 2 ) tan θ ς , 2 2
Now, for the given waypoints, the straight line can be obtained by resolving Equation (73) so that the guidance path from point p 1 to point p 2 can be designed. However, Equation (73) is the nonlinear multivariate system so that it is unable to obtain the closed-form solution. It is favorable to employ a numerical solver.

4.2. Path Tracking

In this subsection, a nonlinear UAV path tracking algorithm similar to the line-of-sight guidance algorithm is addressed to follow the Dubins path explained previously for the waypoint guidance. The kinematic equations of motion for UAVs with respect to time are simply given by [18]
(74) X ˙ = V a cos ψ cos γ (75) Y ˙ = V a sin ψ cos γ (76) Z ˙ = V a sin γ
where the airspeed V a is assumed to be constant, and γ = θ α is the flight path angle of the UAV.
Now, let us assume that there is a Dubins path, Γ , to be tracked. With respect to the distance variable, ς , the kinematic equations of the Dubins path can be also expressed as
(77) d X r ( ς ) d ς = cos ψ r ( ς ) cos γ r ( ς ) (78) d Y r ( ς ) d ς = sin ψ r ( ς ) cos γ r ( ς ) (79) d Z r ( ς ) d ς = sin γ r ( ς )
where ψ r and γ r are the reference angles for the UAV to be tracked, and X r , Y r and Z r are the reference locations to follow the Dubins path.
Also, let us define a straight line Ω tangent to the path Γ at a given point X r ( ς * ) , Y r ( ς * ) , Z r ( ς * ) . A reference point on the line Ω can be defined as
(80) X l ( ς * ) = X r ( ς * ) + 1 / K p cos ψ r * cos γ r * (81) Y l ( ς * ) = Y r ( ς * ) + 1 / K p sin ψ r * cos γ r * (82) Z l ( ς * ) = Z r ( ς * ) + 1 / K p sin γ r *
where ς * indicates the location that is closest to the UAV, and K p is the positive parameter.
When the positive parameter K p is chosen close to zero, the reference point ( X l , Y l , Z 1 ) on the line Ω will be located very far from the point X r ( ς * ) , Y r ( ς * ) , Z r ( ς * ) . Meanwhile, when K p is selected as almost K p , then the reference point is identical with the point X r ( ς * ) , Y r ( ς * ) , Z r ( ς * ) . In this work, the reference point is the actual point to be employed for the path tracking algorithm by adjusting K p .
If the UAV is located far from the guidance route, approaching the closest Dubins path could be the first priority for tracking the path. To accomplish this option, the parameter K p needs to be chosen as a large positive value so that the reference point can be identical with the point closest to the UAV. On the other hand, if the UAV is located close to the guidance route, heading toward the tangent direction of the Dubins path could be the priority to track the path appropriately. For this scenario, it can be achieved by selecting K p close to zero. Thus, adjusting K p with respect to the errors between the UAV position and the point X r ( ς * ) , Y r ( ς * ) , Z r ( ς * ) is very important to track the given route.
Therefore, the attitude command angles of ψ d , γ d are chosen with respect to the distance error between the reference point, X l , Y l and Z 1 and the UAV’s actual location. Thus, the relevant relationships are given geometrically as
(A9) cos ψ d = e ˜ x e ˜ x 2 + e ˜ y 2 , sin ψ d = e ˜ y e ˜ x 2 + e ˜ y 2 (84) cos γ d = e ˜ x 2 + e ˜ y 2 e ˜ , sin γ d = e ˜ z e ˜
where, e ˜ x = X l X , e ˜ y = Y l Y , e ˜ z = Z l Z , e ˜ = e ˜ x 2 + e ˜ y 2 + e ˜ z 2 .
As is known, most fixed-wing UAVs allow roll commands to follow the heading commands. Thus, the path in the lateral dimension for fixed-wing UAVs is usually controlled by a coordinated turn maneuver by controlling the ailerons to achieve the desired roll command ϕ d . When the UAV is banked at an angle ϕ , the lift F L by the wings can be divided into two components. For the successive coordinated turn maneuver, the requirement is given by
(85) F L sin ϕ = m a c m d (86) F L cos ϕ = m g
where a c m d is the lateral acceleration command. Furthermore, manipulating the two requirements leads simply to
tan ϕ = a c m d g
In order to obtain the roll command from the yaw command, let us assume a reference point ξ r on the desired trajectory shown in Figure 3. The distance between the reference point and the UAV is denoted as l, and the angle between the velocity vector ξ ˙ and the line l is the yaw error defined as ψ e = ψ ψ d . Then, the lateral acceleration command is given by [19]
a c m d = 2 V a 2 l cos ψ e
By inserting the acceleration command into Equation (87), the roll command ϕ d is expressed as
ϕ d = tan 1 2 V a 2 g l cos ψ e
Note that the pitch command is also evaluated by θ d = α + γ d . From the evaluated three-axis attitude command represented by Euler angles ( ϕ d , θ d , ψ d ) , the quaternion command, q c , is generated to employ the suggested control law to follow the given guidance route.

5. Numerical Simulations, Results, and Discussions

In this section, the proposed attitude control law and the path-following algorithm are examined using numerical simulations. Let us first explain the simulation conditions. The US 25e model [16,20] is employed for the numerical model of the fixed wing UAV in this work. The specifications of the model are listed in Table 1. The aerodynamic coefficients in Equations (11)–(18) for a given airspeed V a = 20 m/s are also listed in Table 2.
The proposed constrained sliding mode control (CSMC) approach is compared with the conventional SMC in this simulation. The gains for each control law are listed in Table 3 and Table 4, respectively. Note that by employing the angular rate constraint, ω m = 10.0 deg/s chosen for this simulation scenario, there is a slight modification of design parameters between the two control laws. Furthermore, to highlight the excellent performance of the proposed CSMC, wind dust is added to the simulation environment for the time from t = 25 s to t = 40 s. To simulate the disturbance varying sinusoidally, the external moment ranging from −0.2 Nm to 0.2 Nm shown in Figure 4 is applied to the UAV.
Now, let us select some reference waypoints to design the flight path for the UAV. To examine properly the control capability for CSMC, the randomly chosen waypoints and heading unit vectors are listed in Table 5. Then, the Dubins path generation algorithm is applied to the waypoints to design the guidance trajectory. Since the airspeed ( V a ) and the angular rate limit ( ω m ) of the UAV are chosen as 20 m/s and 10 deg/s, respectively, the turning radius for the Dubins circle can be also evaluated as r m = V a / ω m for this scenario. By applying these parameters and the given waypoints above, the guidance trajectory route from point 1 to point 5 can be calculated simply using Equation (73). The 3-D trajectory generated using only Dubins circles and lines results in Figure 5.
In order to verify the tracking capability of the proposed control law, the initial position of the UAV is chosen to be relatively far from waypoint 1. The resultant tracking trajectories by CSMC and SMC laws are indicated in Figure 6. It can be seen that both control laws are able to trace the reference trajectory without difficulty, even though the UAV is initially located far away from waypoint 1. The slight difference between the two tracking paths controlled by the two techniques is evident at the very beginning. This difference originates from the angular rate limitation. The detailed reasons for this are given below in the inspection of the steady state variables. Furthermore, it can be seen that the tracking path is slightly deviated on the first straight path, which shows that the UAV is pushed back due to wind gusts. Except for the first straight path, it can be noted that the UAVs controlled by means of SMC and CSMC follow the curved path and other straight paths without significant differences.
The detailed plots of the steady state variables for the given scenario are displayed in Figure 7, Figure 8 and Figure 9. Extracted from the attitude commands by Equations (83), (84) and (89) to track the designed route, the relevant UAV attitude commands are plotted in Figure 7. The attitude tracking histories using the control laws are also displayed in the Figure as well. Let us focus on the very beginning of the plots. It is obvious that in the case of SMC the attitude command in Figure 7a is quickly followed without errors, whereas in the case of CSMC the attitude tracking time illustrated in Figure 7b is slower than that of SMC due to the limitation of the angular rate. Then, let us inspect the angular rate of the UAV. Figure 8 is the graph of the angular rate. As expected, unrealistic large angular rates take place in Figure 8a to compensate for the generated large attitude errors in the beginning for the UAV controlled by means of SMC. Meanwhile, the angular rate of the UAV controlled by means of CSMC over the simulation time is bounded by the given maximum angular rate, which is plotted in Figure 8b. It is believed that several control transitions have occurred to ensure that the angular velocity does not exceed the angular velocity limit.
Lastly, by examining the control input histories applied to the UAV, it can be highlighted that the proposed CSMC has clear advantages for actual applications. When SMC is utilized, very unrealistic control inputs are applied to the UAV to compensate for the attitude errors introduced by the relatively large attitude command at the beginning of the scenario, which is plotted in Figure 9a. The CSMC approach proposed in this work does not generate the excessive deflection of the control inputs shown in Figure 9b. This characteristic originates from the control law restricting the angular rate, even though very large attitude commands are applied. Therefore, it is believed that the proposed technique providing the permissible control input can be applicable to actual applications.
Of course, the maneuverability of UAVs can deteriorate, when the proposed CSMC is chosen to be utilized primarily for the attitude control of UAVs. However, depending upon whether the emphasis is placed on significant emergencies for structural damage, cracks or safe operation as the first priority to accomplish complex missions, there are several options available for control.A numerical simulation according to the angular velocity limit is performed using CSMC using an identical scenario to that which has been conducted previously. Without the loss of generality, only the tracking histories are examined in this simulation study. Consequently, it can be noted in Figure 10 that when the proposed CSMC is applied, it is evident that the attitude maneuverability can even be identical to that of SMC. Therefore, the control law proposed can be applicable to actual applications by adjusting the angular rate limit, properly satisfying the attitude control capability requirements.

6. Conclusions

A sliding mode-based attitude controller was proposed in the consideration of angular rate constraints for UAVs in this work. The sliding surfaces suggested have two different equilibrium points. These equilibrium points were evaluated through stability analysis.One of these points is used for maneuverability improvements within the allowable maximum angular rate for the safe operation of UAVs; the other is used to determine the attitude control capability. Switching between the two equilibrium points defined by the sliding surfaces is conducted very smoothly according to the magnitude of the attitude errors. Numerical simulations were performed on the path generated by means of the Dubins algorithm in the presence of external disturbances. It was demonstrated clearly that the guidance path generated was followed by the UAVs using the proposed technique, while satisfying the angular rate constraints. Consequently, it is also believed that the proposed technique is appropriate for actual applications to accomplish complex missions safely by operating the UAV within the allowable structural force and moment conditions.

Author Contributions

Conceptualization, H.L.; methodology, H.L.; software, S.-h.J.; validation, S.-h.J., Y.Y. and H.L; formal analysis, S.-h.J., Y.Y. and H.L.; investigation, Y.Y.; resources, H.L.; data curation, S.-h.J. and H.L.; writing—original draft preparation, S.-h.J. and Y.Y.; writing—review and editing, H.L.; visualization, S.-h.J. and Y.Y.; supervision, H.L.; project administration, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study is a part of the “Leaders in INdustry-university Cooperation +” project, supported by the Ministry of Education and National Research Foundation of Korea. Also, this work (Grants No. S2905461 ) was supported partially by project for Cooperative R&D between Industry, Academy, and Research Institute funded Korea Ministry of SMEs and Startups in 2020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zulu, A.; John, S. A review of control algorithms for autonomous quadrotors. Open J. Appl. Sci. 2014, 4, 547–556. [Google Scholar] [CrossRef] [Green Version]
  2. Zhang, R.; Cao, S.; Zhao, K.; Yu, H.; Hu, Y. A Hybrid-Driven Optimization Framework for Fixed-Wing UAV Maneuvering Flight Planning. Electronics 2021, 10, 2330. [Google Scholar] [CrossRef]
  3. Warsi, F.A.; Hazry, D.; Ahmed, S.F.; Joyo, M.K.; Tanveer, M.H.; Kamarudin, H.; Razlan, Z.M. Yaw, Pitch and Roll controller design for fixed-wing UAV under uncertainty and perturbed condition. In Proceedings of the 2014 IEEE 10th International Colloquium on Signal Processing and its Applications, Kuala Lumpur, Malaysia, 7–9 March 2014; pp. 151–156. [Google Scholar]
  4. Sun, C.; Liu, M.; Liu, C.A.; Feng, X.; Wu, H. An Industrial Quadrotor UAV Control Method Based on Fuzzy Adaptive Linear Active Disturbance Rejection Control. Electronics 2021, 10, 376. [Google Scholar] [CrossRef]
  5. Shtessel, Y.; Edwards, C.; Fridman, L.; Levant, A. Sliding Mode Control and Observation; Springer: New York, NY, USA, 2014; Volume 10. [Google Scholar]
  6. Lungu, M. Auto-landing of fixed wing unmanned aerial vehicles using the backstepping control. ISA Trans. 2019, 95, 194–210. [Google Scholar] [CrossRef] [PubMed]
  7. Kang, Y.; Hedrick, J.K. Linear tracking for a fixed-wing UAV using nonlinear model predictive control. IEEE Trans. Control Syst. Technol. 2009, 17, 1202–1210. [Google Scholar] [CrossRef]
  8. Kahn, A. Adaptive control for small fixed-wing unmanned air vehicles. In Proceedings of the AIAA Guidance, Navigation, and Control Conference, Toronto, ON, Canada, 2–5 August 2010; p. 8413. [Google Scholar]
  9. Yang, Y.; Wu, J.; Zheng, W. Station-keeping control for a stratospheric airship platform via fuzzy adaptive backstepping approach. Adv. Space Res. 2013, 51, 1157–1167. [Google Scholar] [CrossRef]
  10. Ru, P.; Subbarao, K. Nonlinear model predictive control for unmanned aerial vehicles. Aerospace 2017, 4, 31. [Google Scholar] [CrossRef]
  11. Charitopoulos, V.M.; Papageorgiou, L.G.; Dua, V. Multi Set-Point Explicit Model Predictive Control for Nonlinear Process Systems. Processes 2021, 9, 1156. [Google Scholar] [CrossRef]
  12. Garcia, G.A.; Keshmiri, S.S.; Stastny, T. Robust and adaptive nonlinear model predictive controller for unsteady and highly nonlinear unmanned aircraft. IEEE Trans. Control Syst. Technol. 2014, 23, 1620–1627. [Google Scholar] [CrossRef]
  13. Al-Darraji, I.; Piromalis, D.; Kakei, A.A.; Khan, F.Q.; Stojemnovic, M.; Tsaramirsis, G.; Papageorgas, P.G. Adaptive Robust Controller Design-Based RBF Neural Network for Aerial Robot Arm Model. Electronics 2021, 10, 831. [Google Scholar] [CrossRef]
  14. Tseng, M.L.; Chen, M.S. Chattering reduction of sliding mode control by low-pass filtering the control signal. Asian J. Control 2010, 12, 392–398. [Google Scholar] [CrossRef]
  15. Hernandez, J.L.; González-Hernández, I.; Lozano, R. Super-twisting control in a solar unmanned aerial vehicle: Application to solar tracking. In Proceedings of the 2018 International Conference on Unmanned Aircraft Systems (ICUAS), Piscataway, NJ, USA, 12–15 June 2018; pp. 379–384. [Google Scholar]
  16. Paw, Y.C. Synthesis and Validation of Flight Control for UAV; University of Minnesota: Minneapolis, MN, USA, 2009. [Google Scholar]
  17. Hota, S.; Ghose, D. Optimal geometrical path in 3D with curvature constraint. In Proceedings of the 2010 IEEE/RSJ International Conference on Intelligent Robots and Systems, Taipei, Taiwan, 18–22 October 2010; pp. 113–118. [Google Scholar]
  18. Ambrosino, G.; Ariola, M.; Ciniglio, U.; Corraro, F.; De Lellis, E.; Pironti, A. Path Generation and Tracking in 3-D for UAVs. IEEE Trans. Control Syst. Technol. 2009, 17, 980–988. [Google Scholar] [CrossRef]
  19. Zhai, R.; Zhou, Z.; Zhang, W.; Sang, S.; Li, P. Control and navigation system for a fixed-wing unmanned aerial vehicle. Aip Adv. 2014, 4, 031306. [Google Scholar] [CrossRef]
  20. Ren, Z.; Fu, W.; Li, Y.; Yan, B.; Zhu, S.; Yan, J. Enhanced attitude control of unmanned aerial vehicles based on virtual angular accelerometer. IEEE Access 2019, 7, 104330–104343. [Google Scholar] [CrossRef]
Figure 1. The coordinate system of a UAV.
Figure 1. The coordinate system of a UAV.
Electronics 10 02776 g001
Figure 2. Geometry of the Dubins path connecting two waypoints of ξ 1 and ξ 2 .
Figure 2. Geometry of the Dubins path connecting two waypoints of ξ 1 and ξ 2 .
Electronics 10 02776 g002
Figure 3. Nonlinear lateral control for the UAV.
Figure 3. Nonlinear lateral control for the UAV.
Electronics 10 02776 g003
Figure 4. Scenario of external disturbances applied to the UAV.
Figure 4. Scenario of external disturbances applied to the UAV.
Electronics 10 02776 g004
Figure 5. A 3-dimensional Dubins path constructed by circles and lines for the five waypoints.
Figure 5. A 3-dimensional Dubins path constructed by circles and lines for the five waypoints.
Electronics 10 02776 g005
Figure 6. Tracking lines of the UAVs controlled by means of SMC and CSMC, respectively.
Figure 6. Tracking lines of the UAVs controlled by means of SMC and CSMC, respectively.
Electronics 10 02776 g006
Figure 7. Time histories of Euler angles.
Figure 7. Time histories of Euler angles.
Electronics 10 02776 g007
Figure 8. Angular rate responses of the UAV.
Figure 8. Angular rate responses of the UAV.
Electronics 10 02776 g008
Figure 9. The deflection angles of control surfaces generated by the control laws.
Figure 9. The deflection angles of control surfaces generated by the control laws.
Electronics 10 02776 g009
Figure 10. Tracking histories conducted by CSMC according to the angular rate limit variations.
Figure 10. Tracking histories conducted by CSMC according to the angular rate limit variations.
Electronics 10 02776 g010
Table 1. Characteristics for US 34e.
Table 1. Characteristics for US 34e.
ParameterValueParameterValue
m1.9 kgb1.27 m
S0.31 m 2 c ¯ 0.25 m
J x x 0.089 kg m 2 J y y 0.14 kg m 2
J z z 0.16 kg m 2 J x z 0.014 kg m 2
Table 2. Parameters of aerodynamic coefficients for US 25e.
Table 2. Parameters of aerodynamic coefficients for US 25e.
ParameterValueParameterValueParameterValue
C L 0 0.23 C L α 4.58 C L δ e 0.13
C L α ˙ 1.97 C L q 7.95 C D 0 0.043
C D δ e 0.014 C D δ r 0.03 C y β −0.83
C y δ r 0.191 C y p 0 C y r 0
C l β −0.04 C l δ a 0.068 C l δ r 0.017
C l p −0.41 C l r 0.4 C m 0 0.135
C m α −1.5 C m δ e −1.13 C m α ˙ −10.4
C m q −50.8 C n β 0.034 C n δ a −0.012
C n δ r −0.035 C n p −0.075 C n r −0.41
Table 3. Design parameters for SMC.
Table 3. Design parameters for SMC.
ParameterValueParameterValue
a12 ϵ 0.95
k 1 2.5 k 2 4.5
Table 4. Design parameters for CSMC.
Table 4. Design parameters for CSMC.
ParameterValueParameterValue
a8 ϵ 0.95
k 1 2 k 2 5.5
Table 5. The selected waypoints and heading unit vectors.
Table 5. The selected waypoints and heading unit vectors.
Waypont Index ξ i (meter) ν ¯ i
1 [ 0.000 0.000 100.00 ] T [ 0.8192 0.5736 0.0000 ] T
2 [ 1000.0 400.0 80.0 ] T [ 0.9848 0.0000 0.1736 ] T
3 [ 700.0 500.0 95.0 ] T [ 0.8627 0.4981 0.0872 ] T
4 [ 500.0 0.000 110.0 ] T [ 0.4924 0.8529 0.1736 ] T
5 [ 100.0 600.0 100.0 ] T [ 0.8192 0.5736 0.0000 ] T
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jang, S.-h.; Yang, Y.; Leeghim, H. Angular Rate Constrained Sliding Mode Control of UAVs for Path Following. Electronics 2021, 10, 2776. https://doi.org/10.3390/electronics10222776

AMA Style

Jang S-h, Yang Y, Leeghim H. Angular Rate Constrained Sliding Mode Control of UAVs for Path Following. Electronics. 2021; 10(22):2776. https://doi.org/10.3390/electronics10222776

Chicago/Turabian Style

Jang, Seok-ho, Youyoung Yang, and Henzeh Leeghim. 2021. "Angular Rate Constrained Sliding Mode Control of UAVs for Path Following" Electronics 10, no. 22: 2776. https://doi.org/10.3390/electronics10222776

APA Style

Jang, S. -h., Yang, Y., & Leeghim, H. (2021). Angular Rate Constrained Sliding Mode Control of UAVs for Path Following. Electronics, 10(22), 2776. https://doi.org/10.3390/electronics10222776

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