Open access peer-reviewed chapter

Attitude Control of a Quadcopter Using Adaptive Control Technique

By Ramiro Ibarra Pérez, Gerardo Romero Galvan, Aldo Jonathan Muñoz Vázquez, Silvia Florida Melo and David Lara Alabazares

Submitted: August 9th 2017Reviewed: October 2nd 2017Published: December 20th 2017

DOI: 10.5772/intechopen.71382

Downloaded: 207

Abstract

This chapter presents an adaptive control technique to stabilize the attitude dynamics of unmanned aerial vehicle (UAV) type quadrotor in the presence of disturbances and/or uncertainties in the parameters due to changes in the payload, nonlinear actuators, and change in environmental conditions. To address the above problem, MRAC (model reference adaptive control) strategy is used. In this schema, a cost function is defined as a function of the error between the output of the system and a desired response from a reference model. Based on this error, the controller parameters are updated. To guarantee the global asymptotic stability of the system, Lyapunov’s theory is applied. Simulation results using MATLAB-Simulink platform are presented to demonstrate the convergence of the controller parameters.

Keywords

  • Quadrotor UAV
  • adaptive control
  • Lyapunov’s method
  • aircraft
  • attitude control

1. Introduction

Unmanned aerial vehicles (UAVs) are aircrafts that do not require a pilot on board to be controlled. In the beginning, they were solely used for military purposes. One of the first applications of these vehicles was aerial photography. In the 1883, an Englishman named Douglas Archibald provided one of the world’s first reconnaissance UAVs. However, it was not until the World War I that UAVs became recognized systems. Since then, they are being widely used in military missions such as surveillance of enemy activity, airfield base security, airfield damage assessment, elimination of unexploded bombs, etc.

In the last decades, advances in technology and costs reduction permitted to adopt this technology in civil applications such as aerial photography, video and mapping, pollution and land monitoring, powerline inspection, fire detection, agriculture, and among other applications [1].

Quadrotors are a kind of mini-UAV’s with vertical take-off and landing, controlled completed through four rotors mounted on each end of the crossed axes, which provide lift forces for the quadrotor move, this vehicle is very popular in the research community due to their special features like strong coupling subsystems, unknown physical parameters, and nonparametric uncertainties in inputs and external disturbances. Therefore, a suitable control system for stabilizing the closed loop control system is required; to do this, various control techniques, linear and nonlinear, have been implemented such as control PD [2, 3], control PID [4, 5], control of position and orientation by vision [6], sliding mode control [1, 7], fuzzy logic [8, 9], and adaptive control in [10].

The dynamic behavior of quadrotor has been published in varying complexity [11, 12]; particularly, the model we used is based on [13], where an extended mathematical description by the full consideration of nonlinear coupling between the axes is presented. We assume elastic deformations sufficient stiffness and realistic flight maneuvers were omitted, mass distributions of the quadrotor are symmetrical in the x-y plane, drag factor and thrust factor of the quadrotor are constant, and air density around the quadrotor is constant.

In this work, we develop an adaptive control strategy to stabilize the attitude dynamics of a quadrotor UAV. The adaptive control permit deals with modeling errors and disturbance uncertainty, variations of the mass, inertia, actuators faults, nonlinear aerodynamics, etc.

This chapter is organized as follows: Section 2 presents the mathematical model of the quadrotor obtained using Newton-Euler equations. Section 3 provides brief introduction about the adaptive control theory and the methodology used is based to obtain the adaptive law equations. In Section 4, simulations and analysis are presented to verify the performance of MRAC schema, and finally in Section 5, conclusions of this work are presented.

2. Mathematical model

2.1. ESC + motor + propeller

The rotors of quadrotor are modeled by

Ti=KωMs+ωMuiE1
vi=ωMs+ωMuiE2

For i = (1, 2…, 4), where ui represents the PWM input, the gain K > 0 and ωM are the rotor bandwidth, vi is the actuator dynamics. We assume K and ωM are the same for all rotors [5].

2.2. Description of Quadrotor system

The configurations of quadrotor are described in Figure 1 and consist of a body-fixed frame denoted by Fb = {Fbx , Fby  , Fbz }, the center mass that coincident with the body-frame origin denoted by Cm , the earth-fixed inertial frame is denoted by Fe = {Fex, Fey , Fez}, the position in earth-fixed inertial frame is given by a vector rI  = (x, y, z) T , four rotors denoted each one as Mi with i = (1, 2, …, 4).

Figure 1.

Quadrotor configuration.

To obtain the position of CM in the NED (North-East-Down) coordinate system with respect to the inertial frame Fe , we use the Euler angles (roll ϕ, pitch θ, yaw ψ) combined in a vector Ω = (ϕ, θ, ψ)T, where roll angle is generated by differential thrust betweet M 2 and M 4, pitch angle is generated by differential thrust between M 1 and M 3, and yaw angle is generated by diferential torque between clock wise and anticlock wise rotors, i.e., (M 1 − M 2 + M 3 − M 4). This kind of device is considered an underactuated system with six degree of freedom (DOF). Rotation matrix R is used to map from Fb to Fe, applying three consecutive rotations denoted by

Rz=CψSψ0SψCψ0001Ry=Cθ0Sθ010Sθ0CθRx=1000CϕSϕ0SϕCϕR=RzRyRx=CψCθCψSθSϕSψCϕCψSθCϕ+SψSϕSψCθSψSθSϕ+CψCϕSψSθCϕCψSϕSθCθSϕCθCϕE3

where S and C refer to sin and cos function, respectively.

2.3. Rigid body dynamics

The dynamics of a rigid body under external forces and moments in the f b can be formulated as follows [14]

mv̇=FiE4
JẆ+W×JW=MiE5

where i = (1, 2, …, 4), ∑Fi  ∈ R3represent the sum of thrust forces from each rotor, ∑Mi  ∈ R3represent the sum of moments acting on the vehicle, v is the velocity in the Fb , m is the quadrotor mass, W = (wx, wy, wz) T determines the angular velocity vector which has its origin in at the center of mass of the rigid body, and J is the inertial tensor of the symmetric rigid body around its center of mass.

To obtain the thrust force contribution of each rotor applied to Fb , we use the following formula

Fi=fxi+fyi+fzi=kωi2E6

with fxi force from rotor i applied in Fb respect x axes, fyi force from rotor i applied in Fb respect y axes and fzi force from rotor i applied in Fb respect z axes, ωi as rotor’s angle speed, k as proportional gain related with air density, the geometry of the rotor blade and its pitch angles. As fxi and fyi is zero, since they form an angle of 90 degrees respect to the rotor thrust, rewriting Eq. (6), we have

Fi=kωi2 001TE7

By transforming Fi to F e and making use of the principle of linear momentum, the following equations can be introduced for a quadrotor of mass m under gravity g:

r¨=00gT+R14FimE8

The rotatory moment of the body is described by Eq. (4) can be rewrite as

JẆ+W×JW=MBMGE9

where MB is the vector of external torques and is composed of the thrust differences and drag moments of the individual rotors and under considerations of the rotatory directions and can be calculated by

MB=Lbω22ω42Lbω12ω32dω12ω22+ω32ω42E10

where d is called drag factor and it is related with the air resistance, L is a parameter which represents the length of the lever between center of mass and the four rotors. MG is the vector of gyroscopic torques generated due to the propellers rotational movements and can be calculated by the follow formula

MG=JRW×001ω1ω2+ω3ω4E11

where JR represents the inertia of rotating rotors. The follow equations permit us to determine the rate of change of the Euler angles in the inertial frame Fe

Ω̇=HΦW=1tanθsinϕtanθcosϕ0cosϕsinϕ0sinϕcosθcosϕcosθWE12

The relationship that maps rotor angular velocities to forces and moments on the vehicle is

u=u1u2u3u4=bbbb0b0bb0b0ddddω12ω22ω32ω42E13

where the variable u 1 = F 1 + F 2 + F 3 + F 4. The variables u 2 and u 3 correspond the forces from rotors necessary to generate the pitch and roll moments, and u 4 represents the yaw moment.

The translation motion is obtained combining Eqs. (3), (7), (8), and (13), as a result, we have the follow equations

r¨=x¨y¨z¨=CψSθCϕ+SψSϕu1mSψSθCϕCψSϕu1mgCθCϕu1mE14

On the other hand, if we assume small perturbations in hover flight, Ω̇Wand combining Eqs. (9), (11), and (13), we obtain the follow equations

ϕ¨=θ̇ψ̇JyJzJxJrJxθ̇hu+LJxu2θ¨=ϕ̇ψ̇JzJxJy+JrJyϕ̇hu+LJyu3ψ¨=ϕ̇θ̇JxJyJz1Jzu4E15

2.4. Linearization of quadrotor dynamics

To simplify the earlier equations, we have linearized the rotatory system around hover state assuming small change in Euler angles, u 1 ≈ mg in x and y directions, cos(α) = 1, sin(α) = α and neglected the gyroscopic torques; therefore,

x¨=θθgy¨=ϕgz¨=u1mgϕ¨=LJxu2θ¨=LJyu3ψ¨=1Jzu4E16

This work only considers the angular momentum corresponding to the orientation of the vehicle, so that the displacement in (x, y, z) axes will not be used for mathematical analysis.

Therefore, the transfer function for dynamic correspondence of orientation is expressed as

ϕs=Ls2Jxu2θs=Ls2Jyu3ψs=1s2Jzu4E17

and adding rotor’s model, we have

Gseτs=Kp1s2s+aϕ000Kp2s2s+aθ000Kp3s2s+aψE18

where K p1 = Lkϕ /Jx , K p2 = Lkθ /Jy , K p3 = Lkψ /Jz . These parametres were computed for a X4-flyer and have the following values: K p1 = K p2 = 33.23, K p3 = 16.95,  =  =  = 4.1 and τ = 0.017. This mathematical model will be used to synthesize the control law to stabilize the closed-loop system [15].

The overall control schema is showed in Figure 2 , where position and attitude control are presented. This schema consist in two loops, first one is used to perform the quadrotor tracking of desired trajectory rd , while the second one is used to achieve the desired Euler (ϕd , θd , ψd ). Due that we are only interested in attitude control, the above schema is redrawn as is showed in Figure 3 .

Figure 2.

Overall control system.

Figure 3.

Attitude control system of a quadrotor.

3. Adaptive controller design

The adaptive control is an advance control technique which provides a systematic approach for automatic adjustment of controllers in real time, in order to achieve or to maintain a desired level of control system performance, when the parameters of the plant dynamic model are unknown and/or change in time [16]. Two different approaches can be distinguished: indirect and direct approaches. In the first approach, the plant parameters are estimate online and used to calculate the controller parameters. In the second, the plant model is parameterized in terms of the controller parameters that are estimated directly without intermediate calculations involving plant parameter estimates [17].

The model reference adaptive control or MRAC is a direct adaptive strategy which consists of some adjustable controller parameters and an adjusting mechanism to adjust them. The goal of the MRAC approach is adjusting the controller parameters so that the output of the plant tracks the output of the reference model having the same reference input. The MRAC schema is combine two loops: the inner or primary loop where controller and plant are feedback as in normal loop and outer loop or also called adjustment loop where some adaptive mechanisms and a model reference are used to obtain the some performance [18]. In Figure 4 , an overall MRAC schema is presented.

Figure 4.

MRAC control schema.

This section presents the design of an adaptive controller MRAC using the Lyapunov‘s stability theory. This allows us to ensure the tracking trajectory of an X4 to our reference model and makes the system insensitive to parameter variation and external disturbances, leading the state error to zero. Based on this, the process model can be represented by state space as follows:

ẋ=Apx+BpuE19
y=CxE20

where Ap and Bp represent the matrix and the vector of unknown constant parameters of the system, u is the output signal of the controller, and x is the state vector.

The reference model is defined as follows:

ẋm=Amxm+BmucE21
y=CmxmE22

Then, the control law is selected as:

u=LrucLx=θ4ucθ1x1θ2x2θ3x3E23

where Lr and L are the matrix containing the parameters of the controller, which can be freely selected, and uc is the reference signal.

Substituting (Eq. (22)) in (Eq. (18)), the closed loop system is expressed as follows:

ẋ=ApBpLx+BpLruc=Aθx+BθucE24

Now we introduce the error equation as follows

e=xxmE25

Differentiating the error with respect to time, we obtain:

dedt=dxdtdxmdt=Ax+BuAmxmBmucE26

and adding Amx and subtracting to the left side of the equation

dedt=Ame+AAmx+BBmucE27

The error goes to zero if Am is stable and

AθAm=0E28
BθBm=0E29

If we assume that the closed-loop system can be described by (Eq. (14)), where the matrices A and B depend on the parameter θ, and it is some combination of Lr and L, then we can define the following Lyapunov function for the parameter adaptation law:

V=eTPe+trAθAmTQaAθAm+trBθBmTQbBθBmE30

where PR3x3is a positive definite matrix and V is a positive definite function. If its first derivative to time V is not positive definite, then V is Lyapunov function. Now we can solve the derivative of V to time t.

dVdt=eTPAme+eTAmTPe+2trAAmTQaȦ+PexT+2trBBmTQbḂ+PeucTE31

where Q is a positive definite matrix such that meet the follow equation:

AmTP+PAm=QE32

Therefore, if we choose a stable Am matrix, we will get always a P and Q positive definite matrix.

Accordingly, the derivative to time of function V is

dVdt=eTQeE33

Where the function V is a Lyapunov function negative semi-definite ensures the output error between the real system and the reference model will tend to be zero, and the system is asymptotically stable.

Thus, we obtain the following parameter adaptation laws

dθ1dt=p11e1x1+p12e2x1+p13e3x1E34
dθ2dt=p11e1x2+p12e2x2+p13e3x2E35
dθ3dt=p11e1x3+p12e2x3+p13e3x3E36
dθ4dt=p11e1r+p12e2r+p13e3rE37

Then, the resulting control diagram is shown in Figure 5 .

Figure 5.

MRAC control of quadrotor: block diagram.

4. Simulation and analysis

This section presents several simulations test made to prove the performance of MRAC controller to stabilize a mini-UAV quadrotor. As mentioned before, only orientation dynamic (angle position, angular velocity and acceleration) are considered. Analyzing Eq. (18), it is easy to see that roll, pitch, and yaw dynamics are very similar; for this reason, only roll moment is used as example in simulation.

The test begins considering controller parameters are unknown and by using an online adaptive mechanism to determine the values that permit the convergence of plant response to reference model response. It is important to note that the MRAC approach seeks to keep the tracking error (x − xm ) equal to zero by an adjustment of the controller parameters and do not to seek to identify the real parameters of the plant.

Figures 6 8 shows a comparative between the states response of the plant and the state response of the model reference, where can be observer than all states of the plant converge asymptotical to the reference model states. To verify this, the Figures 9 11 are presented; these figures show the tracking error of states goes to zero. Additionally, in Figure 12 , the reference input is compared to the plant and model reference output.

Figure 6.

Comparison between x 1 and x m1.

Figure 7.

Comparison between x 2 and x m2.

Figure 8.

Comparison between x 3 and x m3.

Figure 9.

Error between x 1 to x m1.

Figure 10.

Error between x 2 to x m2.

Figure 11.

Error between x 3 to x m3.

Figure 12.

Roll response.

Figure 13.

Tuning of controller parameters.

Figure 14.

Simulink diagram.

Figure 13 shows the variations of the controller parameters during the adjustment process. This mechanism started many times is necessary to assure the perfect tracking of the plant to the desire response. Finally, in Figure 14 , the simulation diagram is presented.

5. Conclusions

This work presents an adaptive control technique to stabilize the attitude of a quadrotor UAV using MRAC schema, which requires no information of the plant model. The asymptotic stability was demonstrated using the well-known Lyapunov’s theory, obtaining in this way the adaptation law of the controller parameters. Simulations results demonstrate that the adaptive control approach proposed have a good performance to perform the asymptotic tracking of model reference output. It is important to note that the adaptive mechanism is started only when it is needed.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Ramiro Ibarra Pérez, Gerardo Romero Galvan, Aldo Jonathan Muñoz Vázquez, Silvia Florida Melo and David Lara Alabazares (December 20th 2017). Attitude Control of a Quadcopter Using Adaptive Control Technique, Adaptive Robust Control Systems, Le Anh Tuan, IntechOpen, DOI: 10.5772/intechopen.71382. Available from:

Embed this chapter on your site Copy to clipboard

<iframe src="http://www.intechopen.com/embed/adaptive-robust-control-systems/attitude-control-of-a-quadcopter-using-adaptive-control-technique" />

Embed this code snippet in the HTML of your website to show this chapter

chapter statistics

207total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Matlab-Simulink-Based Compound Model Reference Adaptive Control for DC Motor

By Marian Găiceanu

Related Book

Frontiers in Guided Wave Optics and Optoelectronics

Edited by Bishnu Pal

First chapter

Frontiers in Guided Wave Optics and Optoelectronics

By Bishnu Pal

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More about us