Abstract

In this paper, a fuzzy adaptive output feedback dynamic surface sliding-mode control scheme is presented for a class of quadrotor unmanned aerial vehicles (UAVs). The framework of the controller design process is divided into two stages: the attitude control process and the position control process. The main features of this work are (1) a nonlinear observer is employed to predict the motion velocities of the quadrotor UAV; therefore, only the position signals are needed for the position tracking controller design; (2) by using the minimum learning technology, there is only one parameter which needs to be updated online at each design step and the computational burden can be greatly reduced; (3) a performance function is introduced to transform the tracking error into a new variable which can make the tracking error of the system satisfy the prescribed performance indicators; (4) the sliding-mode surface is introduced in the process of the controller design, and the robustness of the system is improved. Stability analysis proved that all signals of the closed-loop system are uniformly ultimately bounded. The results of the hardware-in-the-loop simulation validate the effectiveness of the proposed control scheme.

1. Introduction

Quadrotor unmanned aerial vehicles (UAVs), as a new product in the field of small UAVs, have become a research hotspot among research and scholars all over the word [15]. The main advantages of quadrotor UAVs, such as flying in any direction, take off and land vertically, and hover at an ideal attitude, make the quadrotor UAVs widely used in more important fields, such as providing medical assistance, transporting special resources, disaster monitoring, and agricultural mapping. However, quadrotor UAVs are a complex physical system with the following characteristics, such as multivariate, nonlinearity, underactuating, and strong coupling, which make it very difficult to design an effective adaptive robust flight controller.

In the recent years, various effective control techniques have been developed for quadrotor UAVs to achieve stabilized or automatic flight, such as adaptive PID linear quadratic regulator (LQR) control [2, 6, 7], LMI-based linear control [8, 9], and control [10, 11]. With the development of intelligence control theory, different kinds of advanced nonlinear control methods, which combine the linear control methods with intelligence control theory, such as feedback linearization control [12], model predictive control [13, 14], adaptive backstepping control [15, 16], adaptive sliding-mode control (SMC) [1720], fault-tolerant control [21], dynamic surface control (DSC) [2225], and adaptive fuzzy control [26, 27], have been proposed to achieve attitude and position trajectory tracking performance of quadrotor UAVs. In [28], a novel neural network-based output feedback controller is developed for a group of quadrotor UAVs. In [29], the prescribed performance backstepping dynamic surface control (DSC) scheme is proposed to solve the problem of trajectory tracking control for a quadrotor UAV with control input saturation. In [5], a fuzzy-based compound quantized control strategy is applied to the Quanser Qball-X4 quadrotor experimental platform, which achieved precise position control and tracking performance.

Among the above control schemes, the backstepping strategy has been widely used in the controller design for nonlinear systems. For instance, in [30], the trajectory tracking controller based on the backstepping approach was developed for the quadrotor model, while the PD control was used to attenuate the effects caused by system uncertainties. In [16], a nonlinear disturbance observer-based backstepping control method has been proposed to address the problem of loss of actuators’ effectiveness. However, one drawback of backstepping is the “explosion of complexity” caused by the recurrent derivation of the virtual control law in each design step. To deal with this problem, the DSC control method has been proposed for a class of nonlinear systems, by introducing a first-order low-pass filter in each design step, and the shortcoming was overcome [22, 3133].

An effective way to deal with the uncertainty of system parameters and unmodeled dynamics is to design an adaptive controller using the general approximation ability of the fuzzy logic system (FLS) and neural networks (NNs) [3436]. The number of adaptive laws which depends on the fuzzy rules or the NN weights will be significantly increased as the number of fuzzy rules or the NN weights increase. To overcome this problem, a new method by estimating the norm rather than each item of the weight vector was proposed in [37, 38]. Thus, the number of adaptive laws is reduced significantly. Actually, the quadrotor UAVs are not only time-varying coupled and nonlinear systems, but also suffer from perturbations such as payload variations and nonlinear friction. Therefore, it is necessary to design a controller with adaptive capability, fast convergence, and robust performance for the quadrotor UAVs.

As a widely used nonlinear control algorithm, the sliding-mode control (SMC) is known for its excellent performance properties for complex high-order nonlinear systems in the presence of uncertain conditions [39, 40]. In [18], the SMC trajectory tracking controller was proposed for quadrotor UAVs by considering the wind perturbations and external disturbance components. In [19], a hierarchical control strategy based on the double-loop integral sliding-mode controller was designed for the position and attitude tracking of quadrotor UAVs with sustained disturbances and parameter uncertainties. Most of the existing literature focuses on using the SMC method to solve the attitude control of quadrotor UAVs instead of the position trajectory tracking control design because the transformed dynamic equation has a preferred form for the attitude control.

However, a major constraint in the controller design of quadrotor UAVs is that all the state variables of the system are required to be measurable. But in practical application, under some unpredictable factors, it will cause the measurement sensor to fail. [4145]. In [41], an adaptive output feedback control scheme has been proposed for a class of uncertain SISO nonlinear systems under the constraint that only the system output can be obtained. In [43], a fuzzy state observer-based control method is designed for an uncertain MIMO nonlinear system, and by using the state observer, the problem of state immeasurability has been solved. Traditionally, the tracking performance in adaptive control schemes has been confined to ensure that the tracking error converges to a residual set, the size of which is determined by the explicit design parameters and some unknown bounded terms. The upper bounds of the tracking error are difficult to calculate, so it is a very practical work to make the prior selection of the tracking performance satisfy certain steady state behavior. In [46, 47], a prescribed performance control scheme has been proposed for a class of nonlinear systems, and by constructing a prescribed performance function, the tracking error of the system was transformed into a new variable to ensure that the convergence rate was no less than a prespecified value, and the steady-state error remains within the prescribed range. However, limited attention has been paid to this issue for the controller design of quadrotor UAVs.

Inspired by the aforementioned discussions, an adaptive output feedback dynamic surface sliding-mode control for a class of quadrotor UAV system is presented where the fuzzy approximators are used to approximate the unknown items of the system. The main contributions of the proposed control scheme are as follows:

Firstly, to our best knowledge, this is the first time to use the dynamic surface control techniques with the sliding-mode method to design and test the robust controller of quadrotor UAVs in the platform of hardware-in-the-loop simulation, leading to a greatly simplified structure of the controller and improved robustness of the system.

Secondly, by introducing performance and error transformed functions in the controller, the tracking error of the quadrotor UAVs is transformed into a new error constraint variable which can ensure the prescribed transient performance of the system.

Thirdly, by estimating the norm of the FLS weights instead of estimating each variables of the weight vector, there is only one parameter needed to be updated at each step. Thus, the computing time is reduced.

Finally, the nonlinear state observer is introduced to predict the unmeasurable state of the quadrotor such as the angular velocity state of the quadrotor. Then, only the measurable attitude and position information are required in the implementation of the controller of the quadrotor.

The rest of this paper is organized as follows. In Section 2, problem statement and preliminaries including the mathematical model of the quadrotor, fuzzy logic systems (FLSs), and performance function are introduced. The process of the controller design and analysis of stability are given in Sections 3 and 4, respectively. Section 5 shows the results of the hardware-in-the-loop simulation to validate the effectiveness of the proposed control scheme.

2. Problem Statement and Preliminaries

2.1. Dynamic Model of Quadrotor UAVs

The schematic configuration of the quadrotor in this paper is shown in Figure 1. The basic movements are vertical movement, front and back movement, left and right movement, pitch rotation, roll rotation, and yaw rotation. On changing the rotor speed altogether with the same quantity, the lift forces will change, in this case, affecting the attitude of the vehicle. The complicated motions of a quadrotor can be divided into two typical parts, and each part represents a subsystem with coupled terms. The first part is associated with the translational positions, and the second part is associated with the rotational angles. And in this section, we will deduce the mathematical model of a quadrotor UAV, including navigation equations and moment equations.

Define and , where , , and represent the roll angle, pitch angle, and yaw angle with respect to the inertia frame and , , and denote the angular velocity of the roll, pitch, and yaw with respect to the body-fixed frame. Let denote the transformation matrix between the inertia frame and the body-fixed frame using Euler–Lagrange formulation, which can be expressed aswhere and denote the and , respectively.

Let denote the position with respect to the inertia frame. According to Newton’s second law, the relationship between combined force and acceleration in the ground coordinate isand we can get the translational dynamic equations of the quadrotorwhere in which are the air drag coefficients; is the lift force generated by the rotors with respect to the body coordinate system, in which denote the rotary speed of the front, right, rear, and left rotors and is the lift coefficient of the rotor.

According to the kinetic equation, the relationship between and can be described aswhere denotes and the transformation matrix is bounded according to for a known constant provided and [23]. is defined as the torque provided by the rotors with respect to the body-fixed frame and is described as follows:where , is a constant, is the distance between a rotor and the center of mass of the quadrotor, denotes the lift provided by the rotor, and we get represents the antitorque coefficient. Using the Newton–Euler equation, we can get the rotational dynamic equation of the quadrotor:where is a symmetric positive definite constant matrix with , , and being the rotary inertia with respect to the , , and axes, the signal represents the cross multiplication, and and are the resultant torques due to the gyroscopic effects and the resultant of the aerodynamic frictions torque. They are given aswhere represents the moment of inertia of each rotor and , and are the corresponding aerodynamic drag coefficients.

From (6), the following equation can be obtained:where

With the help of (4), the following equations can be obtained:where .

Remark 1. It is worth noting to point out that the roll, pitch, and yaw angles are limited to which is physically meaningful.
By combing (3) and (10), a nonlinear equation of the quadrotor UAV is given as follows:where is the state vector and and are the smooth functions. The dynamic of quadrotor UAVs can be described as follows:where is the gravity acceleration and are the control inputs which can be expressed asDividing the unknown parameters into two parts, known part and unknown part , it is expressed as follows:

2.2. Fuzzy Logic Systems (FLSs)

The FLS is composed of three main components: fuzzy rule base, fuzzification, and defuzzification operators. The fuzzy rule base comprises a collection of fuzzy “IF-THEN” rules of the following form:: if is and is , and is , then is , ,where and are the FLS input and output, respectively, is the number of rules, and fuzzy sets and are associated with the fuzzy membership functions and . Through the singleton function, center average defuzzification, and product inference, the FLS can be expressed aswhere . The fuzzy basis function is defined as

Denoting and , FLS (15) can be rewritten as

Lemma 1. According to [34], FLSs can effectively approximate any continuous nonlinear function with any small approximated error in a compact set. It can be expressed aswhere the continuous nonlinear function with is a compact set, is an FLS (17), and is the approximated error. Therefore, can be described aswhere the minimum approximated error and is an ideal weight vector and can be defined as

2.3. Performance and the Error Transformation Functions

Similar to [46], the mathematical expression of the prescribed tracking performance is given bywhere , are the tracking errors, the performance function is defined as a smooth and decreasing positive function, and and are the given positive constants. Moreover, and represent the lower and upper bounds of the undershoot of and and are the maximum allowable size of .

The error transformation function is chosen aswhere is the transformed error variable and is a smooth strictly increasing function with the following properties:

Note that if is kept bounded, we have , and thus (21) holds. The inverse transformation of can be expressed aswhere , , and , are the transformed errors, the output tracking errors, and their corresponding performance functions.

In this paper, we chooseand differentiating (25) yieldswhere are the reference signals and are the output tracking errors. From the properties of the transformation, it is clear that and are bounded and .

Remark 2. It can be seen that a new variable is introduced through the above transformation process (21)–(25). If the designed control system can guarantee that is bounded, then the tracking error is bounded and meets formula (21). This means that the tracking error is always kept within the range or . The control objective is to design an adaptive controller so that is bounded.

2.4. Nonlinear Observer

For a class of nonlinear systems with in the observer canonical form is given bywithwith , , , and is the unknown smooth function. The vector is general and not in a restricted form. Only the output is assumed to be measurable [48]. For the uncertain system (27), the nonlinear state observer is established aswithwhere is the estimation of the state and is the observer gain vector; is chosen so that the characteristic polynomial of is strictly Hurwitz. Thus, for a given a matrix 0, there exists a positive definite matrix such that

The function is the estimation of . In the next section, we choose FLSs to approximate . According to (27) and (29), the observer error can be expressed aswhere

3. The Process of the Controller Design

In this section, an adaptive FLS dynamic surface sliding-mode control scheme is proposed for position and attitude trajectory tracking control. The structure of the proposed control scheme is shown in Figure 2. The recursive design procedure contains two parts. Part 1 is the position trajectory tracking control and part 2 is the attitude trajectory tracking control. Each part contains three design steps, which are shown in Tables 1 and 2. The details of the controller design process are shown in Appendix A.

In Table 1, are the surface errors and are the virtual control signals in Step 1, Step 3, and Step 5, respectively; (T1.3), (T1.9), and (T1.15) represent the virtual control signal pass through a first-order filter to obtain a new variable with the time constant ; (T1.6), (T1.12), and (T1.18) are the adaptive laws, and are the design positive parameters. It should be noted that are the virtual control given by , , and .

In Table 2, are the surface errors and are the virtual control signals in Step 7, Step 9, and Step 11, respectively; (T2.4), (T2.11), and (T2.18) represent the virtual control signal pass through a first-order filter to obtain a new variable with the time constant ; (T2.7), (T2.14), and (T2.21) are the adaptive laws, and are the design positive parameters; and are the observer gain.

It should be note that are the estimations of with , and and are the ideal weight vector and fuzzy basis function vector of FLSs which are used to approximate the unknown continuous nonlinear function at each design step.

Remark 3. For the attitude trajectory tracking control sysytem, , and can be regarded as known and the input can be solved. The denominator of will not cause singularity since the yaw angle is limited to .

Remark 4. In the traditional sliding-mode control method, the existence of the signum function will cause chattering in the control system. In practical applications, the saturation function [49] or the hyperbolic tangent function [50] are generally used to eliminate the chatting phenomenon.

4. Stability and Prescribed Tracking Performance Analysis

First of all, define the filter errorfrom (A.3), (A.13), (A.23), and (A.37). We have

Differentiating the boundary errors yieldsfrom which we have

By the same token, one haswhere , are the continuous functions. From (37), the following inequality holds:

Consider the Lyapunov Function candidate aswhere and are defined in (A.6)–(A.45).

Theorem 1. For the quadrotor system (12) and the state observer (29), the adaptive laws (T1.6), (T1.12), (T1.18), (T2.7), (T2.14), and (T2.21) and the control input (T2.6), (T2.13), and (T2.20) are given in Table 1 and Table 2; if satisfies , then, by properly selecting the design parameters and the time constant of the low-pass filter appropriately, all the signals in the closed-loop system are semiglobal uniformly bounded. In addition, the tracking error of position and attitude angle can converge to an arbitrarily residual set and is always kept in the prespecified cures.

Proof. Please see Appendix B for details.

5. Hardware in the Loop Simulation Results

In this paper, the hardware-in-loop testing platform is used to verify the effectiveness of the proposed control scheme. The experiment environment and the experimental system architecture are shown in Figures 3 and 4, where the following components are included:(i)NI PXIe-1082, the MT real-time simulator (RTS) with Kintex-7 325T FPGA chip and 16bits synchronized analog I/O with a data transfer rate of 1 MS/s. The simulator supports FPGA simulation for the quadrotor UAV system. The device accepts the control signal and calculates the response of the system in real time and outputs to the controller box.(ii)NI PXIe-1071, the MT Rapid Control Prototype (RCP), with Kintex-7 325T FPGA @Xilinx and 16 analog I/O channels with a transmission rate of 1 MS/s. This device is used to realize the real-time running of the control code and send the control signals to the quadrotor UAV simulation model which is running on the MT real-time simulator.(iii)Adapter plate: it is used to realize the signal connection between the RTS and RCP. The RTS, RCP, and the signal adapter board comprise a closed loop experimental system.(iv)Host computer: uses the StarSim RCP software to download the Matlab/simulink quadrotor UAV system model and the control algorithm into the RTS and RCP, respectively.

In this section, the effectiveness and the performance of the proposed adaptive dynamic surface sliding-mode output feedback control method are showed by the following experiments. Different scenarios are considered, including normal case and model uncertainty cases to demonstrate the robustness of the proposed controller\enleadertwodots The parameters of the quadrotor UAV adopted in this study are described in Table 3 [17]. In the experiment, the desired trajectory of the position and yaw angle desired trajectory is chosen as . The performance functions are selected as , with the parameters . . The controller parameters chosen for simulation are , and . To demonstrate the effectiveness of the proposed controller, the following different cases are considered and comparisons are conducted:Case 1: normal case: we assume that there are no uncertainties in the model, and all parameters of the quadrotor are normal. The initial state vector is set to be . And, the fuzzy membership function is adopted as , where and .Cases 2, 3, and 4: uncertainty ( added) in rotary inertia: in these cases, we consider three different model uncertainties and separately added in the yaw axis.

The experimental results are shown in Figures 515. Figures 59 illustrate the comparison experimental results of the tracking trajectory in Case1 between the proposed control method and the traditional PID control method. From Figures 69, it can be seen clearly that all the tracking errors of the position and yaw angles of the proposed control scheme are always kept within the performance function curves. That is, the control method proposed in this paper obtains much better control performance by comparing with the traditional PID control scheme. Figure 10 shows the control signals. Figure 11 shows the response curves of roll and pitch angles. Figures 12 and 13 show the change of six adaptive parameters. Figures 14 and 15 show the trajectories of and , . We can see that the proposed state observer can quickly approximate the output of the system. Figure 16 shows the 3D tracking trajectory with uncertainties and added in the yaw rotating axes. Figure 17 illustrates the results of the tracking error under cases 2, 3, and 4. Also, the maximum value of the tracking error (MVTE) and the root mean square value of the tracking error (RMSVTE) in the steady state () are shown in Table 4. From Table 4, we can see that the proposed control scheme has strong robustness, and even the uncertainty in yaw rotary inertia is up to .

6. Conclusion

In this paper, an adaptive dynamic surface sliding-mode output feedback controller has been proposed for attitude and position control of a class of quadrotor UAVs with consideration of parametric uncertainties and disturbances. By using the norm estimation approach, there is only one parameter which needs to be updated online at each design step regardless of the plant order and input-output dimension. Also, by introducing an error transformed function, the tracking performance of the quadrotor UAV has been achieved. The proposed control scheme can not only eliminate the problem of “explosion of complexity” existing in the backstepping control scheme but also improve the robustness of the system. The results of the hardware-in-loop simulation validate the effectiveness of the proposed control algorithm. Furthermore, our future research work will focus on applying the control algorithm of this paper to the quadrotor UAV test platform for experimental verification.

Appendix

A. The Procedures of Controller Design

The controller design process.

Step 1. After error transformation, let given by (25) be the first error variable. Then, the derivative of can be expressed asAccording to adaptive laws (T1.2), (A.1) can be rewritten aswhich suggests that we choose the virtual control signal as where is a positive constant. Introduce a new state variable , which can be obtained by the following first-order filter:Define the error surface (T1.4), and the time derivative of iswhere variables and are introduced. Since is unknown, we use FLSs to approximate the function :with respect to the unknown optimal weight vector in (A.5), define , and since is unknown, let be the estimation of and . Choosing the following proper sliding surface , consider the first Lyapunov function:where the differential of Lyapunov Function can be found as follows:Using Young’s inequality, it can be verified thatThen, (A.7) can be rewritten asThe stabilization of can be obtained by designing the virtual control (T1.5) and the adaptation law (T1.6), where are the positive constants and is an arbitrarily small positive constant. For the external disturbance encountered in the quadrotor flight process, the sliding surface is added to maintain system stability with . Substituting (T1.5) and (T1.6) into (A.9), we getSimilar design procedures can be used to design adaptive DSC sliding-mode laws for trajectory tracking control of axis position and axis position . Introduce two variables and . The corresponding control laws and adaptive laws are designed as follows.

Step 2. Let given by (26) be the second error variable. Then, the derivative of can be expressed asAccording to (29) and (T1.8), (A.11) can be rewritten aschoosing the virtual control signal where is a positive constant. Introduce a new state variable , which can be obtained by the following first-order filter:Define the error surface (T1.10), and the time derivative of iswhere . FLSs are used to approximate the unknown function :with respect to the unknown optimal weight vector in (A.15), define , and since is unknown, let be the estimation of and . Choosing the following proper sliding surface , consider the second Lyapunov function:Differentiating , we obtainUsing Young’s inequality, it can be verified thatThen, (A.17) can be rewritten asThe stabilization of can be obtained by designing the virtual control (T1.11) and the adaptation law (T1.12), where are the positive constants and is an arbitrarily small positive constant. For the external disturbance encountered in the quadrotor flight process, the sliding surface is added to maintain system stability with . Substituting (T1.11) and (T1.12) into (A.20), we get

Step 3. Let given by (26) be the third error variable. Then, the derivative of can be expressed asAccording to (29) and adaptive laws (T1.14), (A.21) can be rewritten asand the virtual control signal can be chosen as as with being a positive constant. Introduce a new state variable , which can be obtained by the following first-order filter:Define the error surface (T1.16), and the time derivative of iswhere . Utilizing FLSs to approximate the unknown function , we obtainDefine , and let be the estimation of and . Choosing the following proper sliding surface , consider the third Lyapunov function:The differentiation of is as follows:Using Young’s inequality, it can be verified thatThen, (A.27) can be rewritten asThe stabilization of can be obtained by designing the virtual control (T1.17) and the adaptive law (T1.18), where are the positive constants and is an arbitrarily small positive constant. For the external disturbance , the sliding surface is added to maintain system stability with . Substituting (T1.17) and (T1.18) into (A.29), we getBy associating , and , the virtual controllers are obtained asNotably, (A.31) has four degrees of freedom, namely, , , , and . We consider the reference trajectory for yaw angle , which is usually given in advance, and the corresponding DSC sliding-mode law is directly designed in the next section to ensure the rapid convergence of to . Thus, we regarded as known that can be replaced by in the controller, and the degrees of freedom in (A.31) is reduced so that , , and can be solved. The programs are as follows:where is one of the ultimate control laws, and in addition, . Introduce two new state variables , which can be obtained by the following first-order filters:For attitude trajectory tracking control, by taking as the desired attitude trajectory, the design of control laws contains three steps. The attitude dynamic system can be extracted as follows:

Step 4. Let given by (26) be the fourth error variable. Then, the derivative of can be expressed asAccording to (29) and adaptive laws (T2.2) and (T2.3), where , (A.35) can be rewritten aswhich suggests that we choose the virtual control signal as where is a positive constant. Introduce a new state variable , which can be obtained by the following first-order filter:Define the error surface (T2.5), and the time derivative of iswhere . Utilizing FLSs to approximate the unknown function , we obtainDefine , and let be the estimation of and . Choosing the following proper sliding surface , consider the following Lyapunov function:where the differential of Lyapunov function can be found as follows:Using Young’s inequality, it can be verified thatThen, (A.41) can be rewritten asThe stabilization of can be obtained by designing the virtual control (T2.6) and the adaptation law (T2.7), where are positive constants and is an arbitrarily small positive constant. For the external disturbance encountered in the quadrotor flight process, the sliding surface is added to maintain system stability with . Substituting (T2.6) and (T2.7) into (A.43), we getSimilarly, the adaptive DSC sliding-mode laws for trajectory tracking control of pitch angle and yaw angle can be designed as (T2.13), (T2.20), (T2.14), and (T2.21), where , are the output of each first-order filter, and , and are the positive constants.
The derivative of Lyapunov candidate for pitch angle and yaw angle are designed as follows:

B. Proof of Theorem 1

Taking the time derivative of and combing (38) yields

Consider the setswhere is also in a compact set. Then, the continuous functions have maximums on , say, . Thus,and using the following inequalities,which together with (B.1) implies thatwhere and is the smallest eigenvalue of matrix . Applying Young’s inequality and the fact the following inequation can be obtained:where and is the largest eigenvalue of matrix .

By using Young’s inequality, we obtain

Substituting (B.6) and (B.7) into (B.5), we havewhere

Choose the suitable design parameters to make and satisfy . The following conditions hold:where is the positive constant. Thus, we have . Hence, all the signals of the closed-loop system are semiglobal bounded. Particularly, the tracking errors of position and attitude angle can converge to an arbitrarily residual set and are always kept in the prespecified cures. This completes the proof.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

This work was supported in part by the NSF of China under Grant 61673101 and 61304015, the Natural Science Foundation of Jilin Province under Grant 20180201009SF and 20170414011GH, the Thirteenth Five-Year Science Research Plan of Jilin Province under Grant JJKH20200119KJ, and the Jilin Technological Innovation Development Plan under Grant 201831719.