Next Article in Journal
Microgrid Frequency Regulation Based on a Fractional Order Cascade Controller
Next Article in Special Issue
Control Design for Fractional Order Leader and Follower Systems with Mixed Time Delays: A Resilience-Based Approach
Previous Article in Journal
Input-Output Finite-Time Stability of Fractional-Order Switched Singular Systems with D-Perturbation
Previous Article in Special Issue
Finite-Time Synchronization for Fractional Order Fuzzy Inertial Cellular Neural Networks with Piecewise Activations and Mixed Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pitch Angle Control of an Airplane Using Fractional Order Direct Model Reference Adaptive Controllers

by
Gustavo E. Ceballos Benavides
1,2,3,*,
Manuel A. Duarte-Mermoud
3,4,
Marcos E. Orchard
1 and
Juan Carlos Travieso-Torres
5
1
Department of Electrical Engineering, Faculty of Physical and Mathematical Sciences, University of Chile, Av. Tupper 2007, Santiago 8370451, Chile
2
Escuela de Negocios y Tecnología, Universidad Gabriela Mistral, Av. Andrés Bello 1337, Santiago 7500533, Chile
3
Advanced Mining Technology Center, University of Chile, Av. Tupper 2007, Casilla 412-3, Santiago 8370451, Chile
4
Facultad de Ingeniería y Arquitectura, Universidad Central de Chile, Av. Santa Isabel 1186, Santiago 8370292, Chile
5
Department of Industrial Technologies, Universidad de Santiago de Chile, El Belloto 3735, Santiago 9170022, Chile
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(4), 342; https://doi.org/10.3390/fractalfract7040342
Submission received: 4 March 2023 / Revised: 13 April 2023 / Accepted: 17 April 2023 / Published: 20 April 2023
(This article belongs to the Special Issue Robust and Adaptive Control of Fractional-Order Systems)

Abstract

:
This paper deals with the longitudinal movement control of an airplane (pitch angle) using fractional order adaptive controllers (FOACs). It shows the improvements achieved in the plane’s behavior, in terms of the minimization of a given performance index. At the same time, less control effort is needed to accomplish the control objectives compared with the classic integer order adaptive controllers (IOACs). In this study, fractional order direct model reference adaptive control (FO-DMRAC) is implemented at the simulation level, and exhibits a better performance compared with the classic integer order (IO) version of the DMRAC (IO-DMRAC). It is also shown that the proposed control strategy for FO-DMRAC reduces the resultant system control structure down to a relative degree 2 system, for which the control implementation is simpler and avoids oscillations during the transient period. Moreover, it is interesting to note that this is the first time that an FOAC with fractional adaptive laws is applied to the longitudinal control of an airplane. A suitable model for the longitudinal movement of the airplane related to the pitch angle θ as the output variable with the lifter angle (δe) as the control variable, is first analyzed and discussed to obtain a reliable mathematical model of the plane. All of the other input variables acting on the plane are considered as perturbations. For certain operating conditions defined by the flight conditions, an FO-DMRAC is designed, simulated, and analyzed. Furthermore, a comparison with the implementation of the classical adaptive general direct control (relative degree ≥ 2) is presented. The boundedness and convergence of all of the signals are theoretically proven based on the new Lemma 3, assuring the boundedness of all internal signals ω(t).

1. Introduction

Controlling the pitch angle θ ( t ) (attitude) of an aircraft (see Figure 1) is relevant in the aviation industry, since an important part of the pilots’ tasks is to maintain a specific attitude, that is to say, to achieve a straight and level flight, as well as to ascend or to descend with a certain angle θ (attitude) with respect to the horizon. Since these tasks require the pilot to be diligent, most sophisticated aircraft have an autopilot attitude to complete the job and achieve the following goals:
-
To relieve the pilot from manipulating the controls, reducing the loads on the plane, and improving navigation accuracy;
-
To fly the airplane without direct control over the longitudinal control surfaces (elevator or tail).
Small planes widely use classical controllers, such as proportional-integral-derivative (PID) controllers because of their simplicity to maintain the airplane under control [1,2,3,4,5]. The main disadvantage is that they do not have adaptation capabilities to face large variations in the operating conditions, external perturbations acting on the plane, or eventual plane parameter variations, since PID parameters are fixed.
In controlling plane trajectories, rather simple controllers, such as the integer order (IO) PID have been reported in the control literature. Only a few attempts have been made to use adaptive or more complex control strategies [6,7,8,9]. Moreover, techniques, such as μ-synthesis [10] and nonlinear control [11] have also been used to control airplanes with good results. In [11], a tracking controller consisting of feedforward and static state feedback was designed to guarantee uniform asymptotic trajectory tracking. The landing control of an F-18 fighter aircraft based on a PID controller with an active disturbance rejection (ADR) system was studied in [12]. In [13], an LQR controller based on genetic algorithms was designed for the longitudinal control of an F-16 fighter aircraft under different flight regimes. The longitudinal control of a B-1 bomber was analyzed in [14] using an optimal multivariable control approach [15,16,17] of the LQG/LTR type based on an adaptive observer [18,19].
Furthermore, in [9,12], the longitudinal control of an F-16 Falcon fighter aircraft was studied by applying advanced control techniques, such as MRAC in its combined version (CMRAC), which could serve as a basis for comparison for an extension of the CMRAC techniques to the fractional order case (FO-CMRAC).
Additionally, in [20], the longitudinal control using backstepping control was applied to a model X-plane similar to the NexStar plane with rectangular wings. In [21], the longitudinal and lateral control of unmanned vehicles, such as helicopters, were studied. The control approach used is of the adaptive type based on a continuous linearization for different operating points to then apply MRAC.
In [22], a typical MRAC approach was used (integer order plant and integer adaptive control laws) for the longitudinal control of an F-15 aircraft and the fractional part of this work was the fractional filter that approximates the plant that is newly approximate to an integer transfer function of the plant. Moreover, in this work, two new blocks were used (dynamic inversion and PI compensator blocks) achieving a more complex implementation.
It is important to mention that no attempt has been reported in the technical literature about the use of FO adaptive controllers for the longitudinal control of an airplane in which fractional adaptive laws are applied to an integer plant (airplane model) with no other additional special blocks.
The paper is organized as follows. In Section 2, the airplane model used in this study is presented for certain flight conditions and controlled by a standard DMRAC. Section 3 is devoted to the description of the DMRAC strategy. Section 4 presents some basic concepts of fractional order calculus that will be used in this study, and a new lemma (Lemma 3) that relaxes the stability condition is proposed. In Section 5, several simulation results are presented with numerical values of the parameters corresponding to the Cessna 182 plane used in this study [2]. Finally, in Section 6, some conclusions are drawn.

2. General Concepts on Plane Dynamics and Flight Control

2.1. Mathematical Model of the Longitudinal Movement of a Plane

In this Section, the basic concepts of a plane’s dynamics and its control are presented. In many cases, plane dynamic motion can be modeled assuming small disturbances concerning the airplane’s static stable path (operating point). The differential equations used in this study represent only the plane’s longitudinal movement, although the complete movement is in three dimensions. These equations are considered linear (or linearized) with constant coefficients corresponding to small deviations in the airplane’s operating point.
We consider the static stable path condition of a straight and level flight. The movement equations must consider the aerodynamic force and moment disturbances. By choosing ω 1 = 0 (angular velocity along the z-axis in rad/s), we obtain the plane small disturbances longitudinal equations separated from the small disturbances directional-lateral equations [2].
Next, for a better understanding of the variables in the equations that follow, Table 1 shows the glossary of terms with their respective units used in this paper.
To simplify the analysis, Figure 1 shows the geometry and the longitudinal angles of interest for the aircraft used in this study.
Figure 1. Fundamental angles of the longitudinal movement of an aircraft (figure uploaded by Baron Johnson [23]).
Figure 1. Fundamental angles of the longitudinal movement of an aircraft (figure uploaded by Baron Johnson [23]).
Fractalfract 07 00342 g001
V is the wind speed relative to the aircraft, θ is the pitch angle, α is the angle of attack, and (it is important to note that this symbol is the same as the order of fractional derivatives, nevertheless this should not be confusing because of the context) γ is the flight path angle. Note that γ = θ α .
The airplane used in this study for the simulations of different adaptive control strategies, in its integer and fractional order versions, is the Cessna-182 utility aircraft, which is very popular due to its low cost and high performance. In addition, this aircraft is widely used in the training process for civil pilots.
Using the plane general dynamic equations and the aerodynamic and trust force disturbances, we obtain the following equations describing the plane’s disturbed longitudinal movement [2,8].
u ˙ t = g θ t cos θ 1 + X u u t + X T u u t + X α α t + X δ e δ e t , U 1 α ˙ t U 1 θ ˙ t = g θ t s i n θ 1 + Z u u t + Z α α t + Z α ˙ α ˙ t + Z q θ ˙ t + Z δ e δ e t , θ ¨ t = M u u t + M T u u t + M α α t + M T α α t + M α ˙ α ˙ t + M q θ ˙ t + M δ e δ e t ,
where u t is the perturbed linear velocity along the longitudinal axis of the aircraft, α t is the angle of attack, θ t is the perturbed pitch angle, δ e t is the tail elevator angle and X i t , Z i t , and M i t are, respectively, the derivatives of cinematics forward, vertical, and moments variables concerning the variable of interest i.
Then, a linearized model around a specified operating point is defined allowing for the introduction of basic and advanced control concepts for the resulting dynamical system. First, it is convenient to perform the following change of notation:
x 1 ( t ) = θ t , x 2 ( t ) = θ ˙ t , x 3 ( t ) = u t , x 4 ( t ) = α t ,
to express the set of Equation (1) in the state space matrix form, represented as:
x ˙ ( t ) = A x ( t ) + B u ( t ) , y t = C x t ,
where A , B , and C are the constant matrices of proper dimensions. Following some algebraic manipulation, we arrive at the following system of equations expressed in state variables:
x ˙ 1 t = x 2 t , x ˙ 2 t = g s i n θ 1 M α ˙ Z α ˙ + U 1 x 1 t + M q + M α ˙ Z q + U 1 Z α ˙ + U 1 x 2 t + M u + M T u M α ˙ Z u Z α ˙ + U 1 x 3 t , x ˙ 3 t = g cos θ 1 x 1 t + X u + X T u x 3 t + X α x 4 t + X δ e δ e t , x ˙ 4 t = g s i n θ 1 Z α ˙ + U 1 x 1 t + Z q + U 1 Z α ˙ + U 1 x 2 t + Z u Z α ˙ + U 1 x 3 t + Z α Z α ˙ + U 1 x 4 t + Z δ e Z α ˙ + U 1 δ e t .
The matrix representation of System (2) becomes
x ˙ 1 x ˙ 2 x ˙ 3 x ˙ 4 = 0 1 0 0 g s i n θ 1 M α ˙ Z α ˙ + U 1 M q + M α ˙ Z q + U 1 Z α ˙ + U 1 M u + M T u M α ˙ Z u Z α ˙ + U 1 M α + M T α M α ˙ Z α Z α ˙ + U 1 g cos θ 1 0 X u + X T u X α g s i n θ 1 Z α ˙ + U 1 Z q + U 1 Z α ˙ + U 1 Z u Z α ˙ + U 1 Z α Z α ˙ + U 1 x 1 x 2 x 3 x 4 + 0 M δ e + M α ˙ Z δ e Z α ˙ + U 1 0 X δ e Z δ e Z α ˙ + U 1 δ e t , y t = 1 0 0 0 x 1 x 2 x 3 x 4 ,
where δ e t is the system input and y t = θ ( t ) is the system output.
The fractional order direct model reference adaptive control (FO-DMRAC) designed for this application will be later compared with its integer order counterpart (IO-DMRAC). The flight conditions will be those corresponding to a straight and level flight, considering the operating conditions shown in Table 2.
For more information on the operational and technical characteristics of this aircraft, the reader is referred to [24].
Additionally, Table 3 shows the values of the derivative coefficients of Equation (1) or Equation (2) for the operating conditions given in Table 2.
Performing a state space analysis, it can be shown that the dynamic model of the longitudinal movement of the aircraft for small variations has the input variable δ e t and the output variable y t for dimension 1.
Next, in Equation (4), the Cessna-182 model is represented in state variables for the operating conditions indicated in Table 1 and Table 2, resulting in
x ˙ 1 x ˙ 2 x ˙ 3 x ˙ 4 = 0 1 0 0 0 6.8485 0 19.2591 32.17 0 0.0456 19.4588 0 0.9877 0.0014 2.2329 x 1 x 2 x 3 x 4 + 0 34.7012 0 0.2162 δ e t , y t = 1 0 0 0 x 1 ( t ) x 2 ( t ) x 3 ( t ) x 4 ( t ) = x 1 t .
Furthermore, since the analysis interest is on the longitudinal movement, and in particular it is desired to control the pitch angle (rotation about the Y axis) of the aircraft, then the transfer function of interest is
W p s = x 1 s δ e s = θ s δ e s .
Therefore, according to Equation (3) or Equation (4), the output will be y t = x 1 t = θ t . Then, applying the Laplace transform to Equation (4) and imposing null initial conditions, the transfer function W p s has the form
W p ( s ) = θ ( s ) δ e ( s ) = C ( s I A ) 1 B + D
where A is a 4 × 4 matrix, B is a 4 × 1 matrix, C is a 1 × 4 matrix, and D = 0.
Then, under these operating conditions, the transfer function of the plant, that is, between the pitch angle θ t and the elevator angle δ e t , is given by
W p s = θ s δ e s = 34.7012 s + 0.0589 s + 2.0996 s 2 + 0.0444 s + 0.02533 s 2 + 9.0826 s + 34.3275 ,
or equivalently,
W p s = θ s δ e s = 34.7012 s 2 + 74.9025 s + 4.2914 s 4 + 9.127 s 3 + 34.7561 s 2 + 1.7542 s + 0.8695 .
Thus, the resulting airplane transfer function is of order 4 ( n = 4 ) and relative degree 2 ( n * = 2 ).
Controlling the pitch angle or attitude of the aircraft θ t is of central importance in aeronautics. Thus, an important pilot’s task is to maintain a specific attitude allowing to achieve a straight and level flight, as well as to ascend or descend with a certain degree θ of attitude concerning the artificial horizon. Since this task is quite demanding for a pilot, most sophisticated aircraft have an autopilot attitude to complete this job. This takes care of two important issues:
-
Relieving the pilot from the excessive manipulation of the controls, reducing the loads developed on the plane, and thus improving the navigation accuracy;
-
Flying the airplane without direct control over the longitudinal control surfaces (elevator).
It should be noted that in the state-of-the-art literature review, no implementations of attitude adaptive controllers (pitch angle control) using fractional order controllers have been reported.

2.2. Control Process Description

Figure 2 shows a simplified schematic block diagram for the DMRAC, in which parameters k and θ are adjusted over time for using their corresponding adaptive laws to maintain the error to be small or zero as t approaches infinity.
In our case, Figure 3 shows in detail the simplified block diagram of the DMRAC which will be used in this study, from analytical and simulation viewpoints. This Figure shows the plant, the reference model, the actuator, and the controller for the relative degree 2 ( n * = 2 ). A complete block diagram implementation of the controller can be seen in [25].
For simplicity, we will consider that the actuator dynamic is fast enough for considering this transfer function as a unity (the same assumption is made for the sensor dynamic which is not shown in the block diagram).
It is also interesting to show the behavior of the system output y p t = θ ( t ) , when the model reference output, denoted as y m t = θ m t , changes over time. Figure 4 shows the evolution of the system output y p t in open-loop when the input δ e ( t ) is a unit step at t = 0.
From Figure 4 we can deduce that for any step change in the reference, the pilot should constantly manipulate the plane for almost 200 s to maintain a straight and level flight, which is a demanding task besides the navigation activities, and this motivates the development of automatic control strategies to help pilot duties.

3. Control Strategies

For the sake of completeness, a brief description of the DMRAC algorithm for plants of relative degree 2 n * = 2 (that will be used to control the plane as an alternative), is given in what follows. For a more detailed explanation of the algorithm, the reader is referred to [25].

DMRAC Algorithm

Generally speaking, the main objective of the MRAC is to minimize the error between the plant output and the model reference output. In the DMRAC, the controller parameter adjustment is based on the control error (difference between the actual and the desired outputs) [25].
For the particular case of the adaptive control of a linear (or linearized) nth order plant with a relative degree 1 n * = 1 , it is used the DMRAC scheme shown in Figure 5, that gives rise to the well-known Error Model 1.
Under very well-known hypotheses [25], the control law for the DMRAC shown in Figure 5 has the form
u t = θ T ( t ) · ω ( t ) ,
where θ t = k t , θ 1 T t , θ 0 t , θ 2 T t T R 2 n and ω t = r t , ω 1 T t , y p t , ω 2 T t T R 2 n are the controller parameters vector and the auxiliary signals vector, respectively, and n is the order of the plant. The auxiliary signals ω 1 R n 1 , ω 2 R n 1 are defined by
ω ˙ 1 t = Λ ω 1 + l u t , ω ˙ 2 t = Λ ω 2 + l y t ,
with k t , θ 0 t , r t , y p t R and θ 1 ( t ) , θ 2 ( t ) , ω 1 ( t ) , ω 2 ( t ) n 1 . ( Λ , l ) is any arbitrary stable and controllable pair, with Λ n 1 × n 1 as an asymptotically stable matrix. For simplicity, in our study, we choose ( Λ , l ) in the controllable canonical form. In this particular case ( n * = 1 ) , the parameters adaptive laws are chosen as:
k ˙ t = s g n k p e 1 t r t , θ ˙ 0 t = s g n k p e 1 t y p t , θ ˙ 1 t = s g n k p e 1 t ω 1 t , θ ˙ 2 t = s g n k p e 1 t ω 2 t ,
and the output error e 1 t = y p t y m t can be expressed as e 1 ( t ) = k p k m W m ( s ) ϕ T ( t ) ω ( t ) , where W m ( s ) is a strictly positive real (SPR) transfer function and ϕ t = θ t θ * is the parametric error, with θ * as the ideal (but unknown) control parameters.
ε t = e 1 t + k t e 2 t e 2 t = θ T t ω ¯ t W m ( s ) θ T t ω ( t ) ω ¯ t = W m ( s ) I 2 n ω ( t )
where e 1 t = y p t y m ( t ) is the control error, ε ( t ) is the augmented control error, e 2 ( t ) is the auxiliary error, ω ¯ ( t ) is a filtered version of ω ( t ) vector, and k 1 ( t ) is an additional adjustable parameter. I 2 n denotes the identity matrix of dimension 2 n .
In our specific case study, since the plant’s relative degree exactly equals 2 ( n * = 2 ) , we can make a small modification to the diagram, as shown in Figure 5, (which is only valid for n * = 1 ), avoiding the introduction of the concept of augmented control error ε ( t ) and auxiliary error e 2 t , simplifying the adaptive controller when n * is exactly 2. Figure 6 shows the modification introduced to the block diagram of Figure 5, with the new block in red.
For the particular case of the longitudinal movement of the airplane under investigation, we will use a fourth order model with a relative degree 2, but with making use of the simpler relative degree equal to 1 implementation, as shown in Figure 5, with the modification included in Figure 6.
Additionally, the corresponding simplified block diagram of Figure 6, for this specific case, the control error e 1 ( t ) can be expressed, as shown in Figure 7.
From Figure 7, the output error can be expressed as e 1 ( t ) = k p k m W m ( s ) L ( s ) ϕ T ( t ) ω ( t ) L ( s ) with L s = s + a , such that W m ( s ) L ( s ) is now an SPR Hurwitz transfer function and ω ¯ t = ω ( t ) L ( s ) = ω ( t ) s + a is a new Hurwitz filtered version of ω ( t ) . In this case, we can use adaptive laws similar to Equation (9) to change ω ( t ) by the filtered version ω ¯ t , that is to say
k ˙ t = s g n k p e 1 t r ¯ t , θ ˙ 0 t = s g n k p e 1 t y ¯ p t , θ ˙ 1 t = s g n k p e 1 t ω ¯ 1 t , θ ˙ 2 t = s g n k p e 1 t ω ¯ 2 t ,
or written in the compact form
θ ˙ t = s g n k p e 1 t ω ¯ t ,
where ω ¯ t = 1 s + a r t y p t ω 1 t ω 2 t and e 1 ( t ) = k p k m W m ( s ) L ( s ) ϕ T ( t ) ω ( t ) L ( s ) .
Furthermore, the control law u t = θ T t ω t + θ ˙ T t ω ¯ t in view of (11) can be expressed as
u t = θ T t ω t s g n k p e 1 t ω ¯ T t ω ¯ t ,
that does not contain the parameter derivative.
Then, choosing L s = ( s + a ) , the realization of P L θ = L s θ t L 1 s [26] is free of drifts, and the subsequent stability test of the modified system under control becomes the essence of the solution to the control problem when the relative degree of the plant is 2 ( n * = 2 ). Prior to ending this section, it is interesting to show some properties of the linear operator P L θ used in the new adjustment control law. It satisfies the following properties:
(i) P L θ T = L s θ T t L 1 s (ii) P L θ * = θ *   if   θ * is constant; (iii) P L θ θ = P L ϕ ϕ ,   where   θ t = θ * + ϕ t

4. Fractional Calculus Overview

In this section, we present some definitions and concepts on the derivative and integral operators of the fractional order [26,27] that will be used in the stability analysis of the implementation of the FO-DMRAC for flight control.
Definition 1 ([27]).
The Riemann–Liouville fractional integral of order α > 0 of a function f t R is defined by
I t 0 α f t = 1 Γ α t 0 t f τ ( t τ ) 1 α d τ ,   t > t 0   and   R e α > 0 ,
where Γ α is the Gamma function defined as
Γ α = 0 t α 1 e t d t .
Definition 2 ([27]).
Let α 0 and n = α . The Caputo fractional derivative of order α of a function f t R is defined as
D t α t 0 C f t = 1 Γ n α t 0 t f n τ ( t τ ) α n + 1 d τ ;   a s   l o n g   a s   f n L 1 t 0 , t .
Some additional tools (lemmas and theorems) that are useful in the stability analysis of the fractional order adaptive control systems are presented in what follows.
Lemma 1 (Principle of fractional comparison).
Let e t R n be a vector of differentiable functions. Then, t t 0 , the following inequality holds [28,29,30,31]
D t α t 0 C e T t P e t 2 e T t P D t α t 0 C e t ,   α 0,1 ,
where P R n × n is a symmetric square matrix of constant coefficients and positive definite. A proof of Lemma 1 can be found in [28].
Another lemma that is useful in the study of the evolution of the output error in FO models is the following:
Lemma 2.
Let e · : R + R be a uniformly continuous and bounded function. If there exists an α 0,1 , such that
1 Γ α t 0 t e τ 2 t τ 1 α d τ < M ,   t t 0 ,   with   M 0 , ,
then
lim t t α ε t 0 t e τ 2 d τ t = 0 ,   ε > 0 .
A proof of this lemma is given in [30].
For completeness, we state Theorem 1 related to the boundedness and convergence of FO dynamical systems.
Theorem 1.
Let the state error e t and the output error e 1 t be represented by equations
D t β t 0 C e t = A m n e t + b m n ϕ T t ω t , e t 0 = e 0 , e 1 t = k p h m n T e t , e 1 t 0 = e 10 ,
where A m n R n × n is a Hurwitz matrix and such that given matrix Q = Q T > 0 R n × n . Then, there exists a matrix P = P T > 0 R n × n , such that
A m n T P + P A m n = Q , P b m n = h m n k p
This implies that the triplet A m n , b m n , h m n satisfies the conditions of the Kalman–Yakubovich–Popov Lemma [25]. Furthermore, k p is an unknown constant, but with a known sign, e t : R + R n corresponds to the vector of the non-accessible states error, e 1 t : R + R is the output error (accessible), b m n , h m n R n with ϕ t : R + R m is the parameter error vector defined as ϕ t = θ t θ * with θ t : R + R m , the estimated parameters (of the controller) and θ * : R + R m , the unknown ideal parameters (of the controller). ω t : R + R m is a vector of available auxiliary signals and β 0,1 is the fractional order of the plant, whose adaptive adjustment laws estimate the unknown controller parameters, which are given by
D t α t 0 C ϕ t = D t α t 0 C θ t = γ s g n k e 1 t ω t , ϕ t 0 = 0
with α < β and α 0,1 . Then, assuming that e t and ϕ t are differentiable and uniformly continuous functions, it holds that
  • The parametric error ϕ t , the state error e t , and the output error e 1 t remain bounded for all time;
  • Furthermore, if the auxiliary signal ω t is bounded, then D t α t 0 C ϕ t and D t β t 0 C e t also remain bounded;
  • The mean value of the squared norm of the state error e t 2 ¯ is t ε α   ε > 0 .
    The proof of this theorem can be found in [32].
Corollary 1.
From Theorem 1, it is evident that if (iii) holds, it must also hold that the mean value of the square norm of the output error e 1 t is t ε α ε > 0 , since e 1 t = h m n T e t with h m n T a vector has components that are constants.
Finally, a new lemma (Lemma 3) for the case when the relative degree is 2 ( n * = 2 ) is stated, that relaxes the hypothesis made on the auxiliary signals ω t in point (ii) of Theorem 1.
Lemma 3.
Let us consider a system of order n and relative degree n * = 2 represented by the state equations
x ˙ p t = A p x p t + b p u t , y p t = h p T x p t ,
whose FO adaptive laws to adjust the controller parameters are given by
D t α t 0 C ϕ t = D t α t 0 C θ t = γ e 1 t ω ¯ t , ϕ t 0 = 0 ,
and whose error equation is given by (in this case, the error equation is of integer order, that is to say, β = 1 with β > α . See Equation (15),
e ˙ t = A m n e t + b r p [ ϕ t ω ¯ t ] , e t 0 = e 0 , e 1 t = h m n T e t , e 1 t 0 = e 10 ,
A p R n × n is the state matrix of the plant, x p R n is the state vector, b p and h p are vectors R n whose input u t and output are y p t R , such that Z p s and R p s are monic polynomials of order n 2 and n, respectively. k p is s constant parameter called the high frequency gain, and its sign is known. Without a loss of generality, we will consider that it is positive throughout the entire analysis, and ( k p > 0 ) and Z p s is Hurwitz. In addition, it is assumed that the parameters of the plant (gain k p and coefficients of the polynomials Z p s and R p s ) are unknown and their reference model is also of the relative degree 2 ( n m * = 2 ) and Hurwitz, where A m n is also Hurwitz. Then, since e t and ϕ t are bounded by part (i) of Theorem 1, the auxiliary signals ω t and ω ¯ t will also be bounded. This is an important conclusion because we do not need to assume that the auxiliary signals ω t or ω ¯ t should be bounded, as is imposed in part (ii) of Theorem 1.
Proof. 
Assuming that e t and ϕ t are uniformly continuous and differentiable, and based on Theorem 1, in which it has been shown that both e t and ϕ t are bounded, then, e ˙ t will also be bounded, since e t is bounded, differentiable, and uniformly continuous. Then, since the auxiliary signals ω t and ω ¯ t are part of the equation for e ˙ t (Equation (18)), then ω t must be bounded and ω ¯ t will also be bounded, since ω ¯ t = 1 s + a ω t is with 1 s + a Hurwitz. This concludes the proof. □
In our case, the control law is u t = θ T t ω t s g n k p e 1 t ω ¯ T t ω ¯ t , but it can also be the classical adaptive control law u t = θ T t ω t , as long as the structure of the error equation is of Equation (18) type, in which case, ω ¯ t is replaced by ω t .
Remark 1.
Lemma 3 can be easily extended to the general case when n * 2 , as long as β of Equation (15) is equal to 1 β = 1 , as is the case of Equation (18).
This lemma allows us to relax the hypothesis that must be considered in Theorem 1, in order to prove that the FO derivatives are bounded, the auxiliary signal ω t should also be bounded, which is not necessary to impose in the case of β = 1 . In other cases, with β 0,1 , only Theorem 1 should be used by now.
Furthermore, if the auxiliary signal ω t is bounded (as shown in Lemma 3), then Theorem 1 guarantees that the squared norm of the state error e t 2 and the output e 1 t 2 tend to 0, as t tends to infinity. Thus, the stability of the proposed FO adaptive control system has been proved.

FO-DMRAC Algorithm

Figure 6 shows the block diagram of the FO-DMRAC for the specific case when n * = 2 . Since the equation of the plant is of integer order, the only equations that change are the adaptive laws. For comparison purposes, it is interesting to note that for the IO case and n * = 2 , these adaptive laws are given by
θ ˙ t = ϕ ˙ t = s g n k p e 1 t ω ¯ t
and for the FO case, these adaptive laws become
D t α t 0 C θ t = D t α t 0 C ϕ t = s g n ( k p ) e 1 t ω ¯ t
The control law in both cases is given by
u t = θ T t ω t s g n k p e 1 t ω ¯ T t ω ¯ t
where ω ¯ t = 1 s + a ω ( t ) with a , an arbitrary scalar is greater than 0. For simplicity, we will choose a = 1 . Furthermore, the high frequency gain of the plant is k p = 34.7012 (see Equation (7)), which is supposed to be unknown, but its sign is assumed to be known ( s g n ( k p ) > 0 ) .

5. Simulation Results and Comparisons

Computer simulations were performed in MATLAB-Simulink [33,34]. For all simulations, zero initial conditions were considered for the airplane. Table 4 shows the implementation details of the IO-DMRACs and the FO-DMRACs.
Fractional adaptive laws were implemented using the Ninteger Toolbox for MATLAB [35]. Specifically, the NID block was used, which is based on the Oustaloup approximation method [36]. In this study, five poles and five zeros were selected in the implementation, and the frequency interval chosen was ω [ 0.01 ; 100] [rad/s]. It is important to mention that better approximations can be achieved if the bandwidth is increased (e.g., [ 0.010 ;   1000 ] at the expense of increasing the simulation time.

5.1. IO-DMRAC v/s FO-DMRAC Using the Particle Swarm Optimization (PSO)

In this Section, we compare the control results using IO-DMRAC v/s FO-DMRAC, when the controller parameters of both implementations are optimized using the particle swarm optimization (PSO) technique [22,37], using 50 particles and 50 iterations. It should be noted that any other optimization technique could also be used to determine the controller’s parameters, and it is at the discretion of the designer. First, we analyze the open-loop response of the plant to a unit step input and zero initial conditions to observe the behavior of the oscillations during the transient period before reaching the steady-state regime. The result is shown in Figure 8.
From the open-loop response shown in Figure 8, we note that the time to reach the steady state regime is almost 300 s, which is considered too long. Therefore, an automatic control system has to be designed to reduce the time of the transient response while tracking the reference accurately.

5.2. Simulation Results Using the Exact Knowledge of n * = 2

In this section, we analyze the behavior of the controller given in Figure 6, i.e., assuming that the relative degree of the plant is exactly equal to 2.
The objective function used to optimize the process behavior is defined as
J n o r m = ω e s s e s s + ω d u t 0 = 0 t d u d t d t + t 0 = 0 t ω e · e 1 2 t σ e t + ω u · u 2 t σ u t d t
where e s s is the control steady-state error, d u d t is the rate of change of the input to the plant (or the controller output), e 1 ( t ) is the control error, ( t ) is the output of the controller, ω e and ω u are the weighted parameters, and σ e t and σ u t are the standard deviations of the control error and controller output as a time function, respectively. The reference is the unit step.
Remark 2.
In this study, a rather general objective function (Equation (22)) was chosen. Nevertheless, any other properly defined objective function could also be used by the designer.
Using the values and equations shown in Table 4 and considering the weights ω e s s = 1 , ω d u = 1 , ω e = 2 , and ω u = 1 , the design parameter values for each case (IO-DMRAC and FO-DMRAC) were obtained and used in the simulations.
Figure 9 shows the evolution of the fitness function (objective function) J n o r m as a function of the number of iterations used in the PSO algorithm.
The optimal parameter values obtained for the case of IO-DMRAC were:
γ k = 5.3633 ,   γ 11 = 0.1 ,   γ 12 = 0.1 ,   γ 13 = 49.58 , γ 0 = 0.1 ,   γ 21 = 0.1 ,   γ 22 = 7.6837   and   γ 23 = 0.1 .
As in the previous case, Figure 10 shows the evolution of the objective or fitness function while increasing the number of iterations of the PSO method in the case of the FO-DMRAC.
The best parameter values obtained for the FO-DMRAC case were
α k = 0.3 ,       α 11 = 0.5315 ,       α 12 = 1 ,       α 13 = 1 , α 0 = 0.3 ,       α 21 = 0.3 ,       α 22 = 1 ,       α 23 = 0.8608 , γ k = 5.6173 ,     γ 11 = 0.03 ,     γ 12 = 0.03 ,     γ 13 = 33.5375 , γ 0 = 0.03 ,       γ 21 = 10.0493 ,       γ 22 = 0.03 ,       γ 23 = 0.11957 .
From Figure 9 and Figure 10, it is interesting to observe a faster rate of convergence to a stable value of the objective function in the FO case (Figure 10) compared to the integer case (Figure 9), even though the number of optimization parameters is twice the one used in the IO case. This seems to be an indication that, since in the FO case there are more degrees of freedom in the controller, the behavior attained is better than the one obtained for the IO case.
Figure 11 and Figure 12 show the responses of both adaptive controllers regarding the tracking of the reference signal (in green). Both responses are similar to FO-DMRAC and IO-DMRAC and approximate the reference very well. Nevertheless, in the case of FO-DMRAC, the response is smoother than in the IO-DMRAC case.
Moreover, Figure 13 and Figure 14 show that the control effort in the IO case is greater than that required in the FO case. Particularly, this effort for the IO case is greater in the transient period.
Table 5 shows the optimal values of the cost function J n o r m when all controller parameters are varied.
From Table 5, it can be seen that a better performance is achieved when using FO-DMRAC, since all of the indices are smaller than the IO case. Furthermore, the control effort is smaller in the transient period, as is shown in Figure 14.
Remember that the weights used in the optimization processes were chosen as ω e s s = 1 , ω d u = 1 , ω e = 2 , and ω u = 1 .

5.3. Simulation Results Using the Generalized Controller for n * 2

In this section, we analyze the problem from a viewpoint different from the one used in Section 5.1, i.e., using the general implementation of the FO-DMRAC when the relative degree is greater than or equal to 2 ( n * 2 ). The idea is to compare the behavior with the one implemented in Section 5.1, which, as mentioned, only works for the case when the relative degree is exactly equal to 2 ( n * = 2 ).
The block diagram of the generalized DMRAC when n * 2 is shown in Figure 15 [25]. In this case, additional errors must be considered, such as the auxiliary error e 2 ( t ) and augmented error ε ( t ) , as well as an additional gain k 1 ( t ) . In this case, W m ( s ) and W p ( s ) are the transfer functions of the reference model and the plant, respectively.
It is important to mention that the structure of Figure 15 is the same for both the integer and fractional adaptive control cases. The only difference lies in the laws for the adjustment of the parameters (IO derivative in one case and FO derivative in the other one).
Then, to make the analysis comparable (when n * = 2 and FO-DMRAC) with the generalized adaptive implementation of the integer order DMRAC (or IO-DMRAC), we have determined the values of the γ’s and α’s optimized by PSO. The resulting values are:
γ k = 5.3633 ,       γ 11 = 0.1 ,       γ 12 = 0.1 ,       γ 13 = 49.58 , γ 0 = 0.1 ,       γ 21 = 0.1 ,       γ 22 = 7.6837 ,       γ 23 = 0.1 .
Table 6 shows the design parameters for the general case IO-DMRAC when the relative degree of the plant is greater than or equal to 2 ( n * 2 ).
A plot of the pitch angle θ ( t ) for the generalized IO-DMRAC case ( n * 2 ; b l a c k   c o l o r ) and FO-DMARC case ( n * = 2 ; r e d   c o l o r ) are shown in Figure 16 together with the output signal of the reference model θ m ( t ) in green.
θ m ( t ) is the reference signal or the desired evolution of the pitch angle (green color). The signal in black is the plant output when the generalized IO-DMRAC implementation is used (i.e., n * 2 ), and finally in red, the plant output is shown when FO-DMRAC is considered using the simple structure ( n * = 2 )
Without being exhaustive, it can be seen from Figure 16 and Figure 17 that the best response and least control effort is achieved with the FO-DMRAC for n * = 2 . Furthermore, the generalized IO-DMRAC case presents oscillations in the transient period, something that does not occur in the cases of the IO and FO controls when n * = 2 .
It is important to point out that the generalized case corresponds to the structure shown in Figure 15 (i.e., n * 2 ) . Furthermore, as discussed in the previous paragraph, the results for the generalized IO-DMRAC-PSO case are supported by the results shown in Table 7, where the best result is achieved with the FO-DMRAC with a relative degree exactly equal to 2 ( n * = 2 ), and with text highlighted in blue.
Finally, Figure 18 shows the boundedness of all auxiliary signals ω ¯ of the FO-DMRAC when n * = 2 . This is also true for the case when n * 2 is in the state of Lemma 3.

6. Conclusions

In this paper, we have analyzed and compared two types of adaptive controllers of the direct type that were applied to an integer order plant using integer and fractional order adjustment laws, namely, the classic or integer DMRAC (IO-DMRAC) and the FO-DMRAC. The simulation results allow us to verify the boundedness of all of the signals in the FO case, as assured by using Lemma 3. Moreover, through simulations, it was possible to show the advantages of the FO implementation, in particular when the relative degree of the plant is exactly equal to 2 ( n * = 2 ). It was also possible to prove that the auxiliary signals ω ( t ) remain bounded without the need of using hypothesis (ii) of Theorem 1, thus relaxing the stability analysis.
It is promising to say that, when possible, it is useful to decrease the relative degree of the plant to 2, since the performance is better than in the cases of larger relative degrees, in which case one is obliged to use the generalized scheme ( n * 2 ) that is more complex to implement due to errors and auxiliary blocks. In addition, the fact that the output signal to be controlled is subject to oscillations and overshoots in the transient period, this phenomenon does not occur when case of n * = 2 is considered. Finally, it is important to note, as said before, this is the first time that an FOAC is used in the longitudinal control of an airplane using fractional-order adaptive laws applied to an integer-order system (airplane model) without changing the classic structure of MRAC. Finally, the analysis of asymptotic stability is a subject still pending for future research work in the field of fractional-order MRAC control systems.

Author Contributions

Investigation, G.E.C.B.; Validation, M.A.D.-M., G.E.C.B.; Resources, M.A.D.-M.; Writing—original draft, G.E.C.B.; Writing—review & editing, M.A.D.-M. and J.C.T.-T.; Supervision, M.A.D.-M.; Funding acquisition, M.E.O. All authors have read and agreed to the published version of the manuscript.

Funding

The research presented in this paper has been funded by CONICYT-Chile, under grant FONDECYT 1190959, FONDECYT 1210031 and Advanced Center for Electrical and Electronic Engineering, AC3E, Basal Project FB0008, ANID.

Data Availability Statement

The data is generated automatically at the simulations level using the model and control laws of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Seckel, E. Stability and Control of Airplane and Helicopters; Department of Aeronautical Engineering, The James Forrestal Research Center School of Engineering and Applied Science, Princeton University: Princeton, NJ, USA, 1964. [Google Scholar]
  2. Roskam, J. Airplane Flight Dynamics and Automatic Flight Controls. Part 2; Roskam Aviation and Engineering Corporation: Lawrence, KS, USA, 1995. [Google Scholar]
  3. Singh, B.; Urooj, S.; Singh, S. Analysis of autopilot system, integrated with modelling and comparison of different controllers with the system. J. Discret. Math. Sci. Cryptogr. 2020, 23, 1059–1068. [Google Scholar] [CrossRef]
  4. Ponrani, M.A.; Godweena, A.K. Aircraft Pitch Control using PID Controller. In Proceedings of the International Conference on System, Computation, Automation and Networking (ICSCAN), Puducherry, India, 30–31 July 2021. [Google Scholar]
  5. Alves, L.; Brito, V.; Palma, L.B.; Gil, B. Position Control in Simulated Airplanes. In Proceedings of the 2019 International Young Engineers Forum, Costa da Caparica, Portugal, 10 May 2019. [Google Scholar]
  6. Blakelock, J. Automatic Control of Aircraft and Missiles; John Wiley & Sons, Inc.: New York, NY, USA, 1965. [Google Scholar]
  7. Modi, V.J.; Mokhtarian, F.; Yokomizo, T. Effect of moving surface on the airfoil boundary-layer control. AIAA-J. Aircr. 1990, 27, 42–50. [Google Scholar] [CrossRef]
  8. Duarte-Mermoud, M.A.; Rioseco, J.S.; González, R.I. Control of longitudinal movement of a plane using combined model reference adaptive control. Int. J. Aircr. Eng. Aerosp. Technol. 2005, 77, 199–213. [Google Scholar] [CrossRef]
  9. Lavretsky, E. Combined/composite model reference adaptive control. IEEE Trans. Autom. Control. 2009, 54, 2692–2697. [Google Scholar] [CrossRef]
  10. Mangiacasale, L. Airplane Control Systems: μ-Synthesis with MATLAB; Edizioni Levroto & Bella: Torino, Italy, 1996. [Google Scholar]
  11. Al-Hiddabi, S.A.; McClamroch, N.H. Tracking and maneuver regulation control for nonlinear nonminimum phase system: Applications to flight control. IEEE Trans. Control Syst. Technol. 2002, 10, 780–792. [Google Scholar] [CrossRef]
  12. Liang, Y.; Chen, X.; Xu, R. Research on Longitudinal Landing Track Control Technology of Carrier-based Aircraft. In Proceedings of the Chinese Control and Decision Conference (CCDC 2020), Hefei, China, 22–24 August 2020. [Google Scholar]
  13. Yuxing, Z.; Jie, X.; Zhikuo, C. LQR Controller for Longitudinal Control Augmentation System of Aircraft Based on Hybrid Genetic Algorithm. In Proceedings of the IEEE CSAA Guidance, Navigation and Control Conference (CGNCC), Xiamen, China, 10–12 August 2018. [Google Scholar]
  14. Lavretsky, E. Adaptive Output Feedback Design Using Asymptotic Properties of LQG/LTR Controllers. IEEE Trans. Autom. Control 2012, 57, 1587–1591. [Google Scholar] [CrossRef]
  15. Doyle, J.C.; Stein, G. Multivariable feedback design: Concepts for a classical/modern synthesis. IEEE Trans. Autom. Control 1981, 26, 4–16. [Google Scholar] [CrossRef]
  16. Safonov, M.; Laub, A.; Hartmann, G. Feedback properties of multivariable systems: The role and use of the return difference matrix. IEEE Trans. Autom. Control 1981, 26, 47–65. [Google Scholar] [CrossRef]
  17. Stein, G.; Athans, M. The LQG/LTR procedure for multivariable feedback control design. IEEE Trans. Autom. Control 1987, 32, 105–114. [Google Scholar] [CrossRef]
  18. Rimon, E.; Narendra, K. A new adaptive estimator for linear systems. IEEE Trans. Autom. Control 1992, 37, 410–412. [Google Scholar] [CrossRef]
  19. Kreisselmeier, G. Adaptive Observers with Exponential Rate of Convergence. IEEE Trans. Autom. Control 1977, 22, 2–8. [Google Scholar] [CrossRef]
  20. Rosario-Gabriel, I.; Cortes, H.R. Aircraft Longitudinal Control based on the Lanchester’s Phugoid Dynamics Model. In Proceedings of the International Conference on Unmanned Aircraft Systems (ICUAS), Dallas, TX, USA, 12–15 June 2018. [Google Scholar]
  21. Sheng, S.; Sun, C.; Duan, H.; Jiang, X.; Zhu, Y.; Shouzhao, S.; Chenwu, S.; Haibin, D.; Xiaoliang, J.; Yansong, Z. Longitudinal and Lateral Adaptive Flight Control Design for an Unmanned Helicopter with Coaxial Rotor and Ducted Fan. In Proceedings of the 33rd Chinese Control Conference, Nanjing, China, 28–30 July 2014. [Google Scholar]
  22. Ynineb, A.R.; Ladaci, S. MRAC Adaptive Control Design for an F15 Aircraft Pitch Angular Motion Using Dynamics Inversion and Fractional-Order Filtering. Int. J. Robot. Control Syst. 2022, 2, 240–252, ISSN 2775-2658. [Google Scholar] [CrossRef]
  23. Johnson, B.J. High Angle of Attack Flight Characteristics of a Small UAV with a Variable-Size Vertical Tail. Ph. D. Thesis, University of Florida, Gainesville, FL, USA, 2009. [Google Scholar]
  24. FAR 23, “Airworthiness Standards: Normal, Utility, Acrobatic and Commuter Category Airplanes”, Part 23, Chapter I, Federal Aviation Regulations, Title 14 Code of Federal Regulations (CFR), 1991. Available online: http://www.faa.gov (accessed on 10 February 2023).
  25. Narendra, K.S.; Annaswamy, A.M. Stable Adaptive Systems; Dover Publications Inc.: Mineola, NY, USA, 2005. [Google Scholar]
  26. Narendra, K.S.; Valavani, L.S. Stable Adaptive Controller Design-Direct Control. IEEE Trans. Autom. Control 1978, 23, 570–583. [Google Scholar] [CrossRef]
  27. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  28. Diethelem, K. The Analysis of Fractional Differential Equations; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  29. Duarte-Mermoud, M.A.; Aguila-Camacho, N.; Gallegos, J.A.; Castro-Linares, R. Using General Quadratic Lyapunov Functions to Prove Lyapunov Uniform Stability for Fractional Order Systems. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 650–659. [Google Scholar] [CrossRef]
  30. Aguila-Camacho, N.; Duarte-Mermoud, M.A.; Gallegos, J.A. Lyapunov Functions for Fractional Order Systems. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 2951–2957. [Google Scholar] [CrossRef]
  31. Aguila-Camacho, N.; Duarte-Mermoud, M.A. Boundedness of the solutions for certain classes of fractional differential equations with application to adaptive systems. ISA Trans. 2016, 60, 82–88. [Google Scholar] [CrossRef]
  32. NAguila-Camacho, N.; Gallegos, J.; Duarte-Mermoud, M.A. Analysis of fractional order error models in adaptive systems: Mixed order cases. Fract. Calc. Appl. Anal. 2019, 22, 1113–1132. [Google Scholar] [CrossRef]
  33. The Math Works Inc. Control System Toolbox User’s Guide. 1998. Available online: http://www.mathworks.com (accessed on 8 March 2022).
  34. Rauw, M.O. FDC 1.2-A Simulating Toolbox for Flight Dynamics and Control Analysis; The Math Works Inc.: Natick, MA, USA, 1998. [Google Scholar]
  35. Valerio, D.; Da Costa, J.S. Ninteger: A non-integer control toolbox for matlab. In Fractional Derivatives and Applications; OIFAC: Bordeaux, France, 2004. [Google Scholar]
  36. Sabatier, J.; Aoun, M.; Oustaloup, A.; Gregoire, G.; Ragot, F.; Roy, P. Fractional system identification for lead acid battery state of charge estimation. Signal Process. 2006, 86, 2645–2657. [Google Scholar] [CrossRef]
  37. Clerc, M. Particle Swarm Optimization; ISTE Ltd.: Hollywood, FL, USA, 2006. [Google Scholar]
Figure 2. Block diagram of the DMRAC scheme.
Figure 2. Block diagram of the DMRAC scheme.
Fractalfract 07 00342 g002
Figure 3. Detailed block diagram of the DMRAC applied to the plane.
Figure 3. Detailed block diagram of the DMRAC applied to the plane.
Fractalfract 07 00342 g003
Figure 4. Plant open−loop response under a unit step input.
Figure 4. Plant open−loop response under a unit step input.
Fractalfract 07 00342 g004
Figure 5. Structure of the DMRAC when the relative degree of the plant is n * = 1 .
Figure 5. Structure of the DMRAC when the relative degree of the plant is n * = 1 .
Fractalfract 07 00342 g005
Figure 6. Block diagram associated to the output error e 1 ( t ) valid only for the case of an nth order plant with a relative degree exactly equal to 2.
Figure 6. Block diagram associated to the output error e 1 ( t ) valid only for the case of an nth order plant with a relative degree exactly equal to 2.
Fractalfract 07 00342 g006
Figure 7. Simplified block diagram of the output error in the case of DMRAC for a plant of a relative degree n * = 2 .
Figure 7. Simplified block diagram of the output error in the case of DMRAC for a plant of a relative degree n * = 2 .
Fractalfract 07 00342 g007
Figure 8. Open−loop system response to a unit step reference change at t = 0.
Figure 8. Open−loop system response to a unit step reference change at t = 0.
Fractalfract 07 00342 g008
Figure 9. Evolution of the objective function J n o r m as the number of iterations increases using PSO for the IO-DMRAC.
Figure 9. Evolution of the objective function J n o r m as the number of iterations increases using PSO for the IO-DMRAC.
Fractalfract 07 00342 g009
Figure 10. Evolution of the objective function or cost function J n o r m as the number of iterations increases using PSO and FO-DMRAC.
Figure 10. Evolution of the objective function or cost function J n o r m as the number of iterations increases using PSO and FO-DMRAC.
Fractalfract 07 00342 g010
Figure 11. Pitch angle response θ t under IO-DMRAC and FO-DMRAC.
Figure 11. Pitch angle response θ t under IO-DMRAC and FO-DMRAC.
Fractalfract 07 00342 g011
Figure 12. Zoom of the variations of the pitch angle θ t shown in Figure 11 by using IO-DMRAC and FO-DMRAC.
Figure 12. Zoom of the variations of the pitch angle θ t shown in Figure 11 by using IO-DMRAC and FO-DMRAC.
Fractalfract 07 00342 g012
Figure 13. Control effort δ e t with FO-DMRAC and IO-DMRAC.
Figure 13. Control effort δ e t with FO-DMRAC and IO-DMRAC.
Fractalfract 07 00342 g013
Figure 14. Zoom of the control effort δ e t with FO-DMRAC and IO-DMRAC.
Figure 14. Zoom of the control effort δ e t with FO-DMRAC and IO-DMRAC.
Fractalfract 07 00342 g014
Figure 15. General structure of the DMRAC for the case when the relative degree of the plant is n * 2 .
Figure 15. General structure of the DMRAC for the case when the relative degree of the plant is n * 2 .
Fractalfract 07 00342 g015
Figure 16. Evolution of the pitch angle θ(t) with the generalized IO-DMRAC (black color) and FO-DMRAC (red color). The reference is shown in green (unit step).
Figure 16. Evolution of the pitch angle θ(t) with the generalized IO-DMRAC (black color) and FO-DMRAC (red color). The reference is shown in green (unit step).
Fractalfract 07 00342 g016
Figure 17. Zoom of the transient evolution of the control effort δ e t during the 10 first seconds. Black signal corresponds to the generalized IO-DMRAC structure case, and IO-DMRAC (blue color) and FO-DMRAC (red color) are the cases when the simple structure is used ( n * = 2 ).
Figure 17. Zoom of the transient evolution of the control effort δ e t during the 10 first seconds. Black signal corresponds to the generalized IO-DMRAC structure case, and IO-DMRAC (blue color) and FO-DMRAC (red color) are the cases when the simple structure is used ( n * = 2 ).
Fractalfract 07 00342 g017
Figure 18. Auxiliary signals ω ¯ (all bounded).
Figure 18. Auxiliary signals ω ¯ (all bounded).
Fractalfract 07 00342 g018
Table 1. Glossary of terms used in the control of the Cessna-182.
Table 1. Glossary of terms used in the control of the Cessna-182.
VariableMeaningVariableMeaning
θ 1   [ r a d ] Stationary pitch angle Z α   [ f t · r a d / s 2 ] Vertical acceleration per unit angle of attack
θ   [ r a d ] Pitch angle or perturbed pitch angle Z α ˙   [ f t · r a d / s ] Vertical acceleration per unit rate change of angle of attack
U 1   [ m p h ] Stationary longitudinal velocity Z q   [ f t · r a d / s ] Vertical acceleration per unit pitch rate
u   [ m p h ] Perturbed longitudinal velocity Z δ e   [ f t · r a d / s 2 ] Vertical acceleration per unit elevator angle
α   [ r a d ] Angle of attack or perturbed angle of attack M u   [ f t · r a d / s ] Pitch angular acceleration per speed unit change
γ   [ r a d ] Path flight angle M T u   [ f t · r a d / s ] Pitch angular acceleration per unit change in speed (due to thrust)
g   [ f t / s 2 ] Acceleration of gravity M α   [ 1 / s 2 ] Pitch angular acceleration per unit angle of attack
X u   [ 1 / s ] Forward acceleration per speed unit change M T α   [ 1 / s 2 ] Pitch angular acceleration per unit angle of attack (due to thrust)
X T u   [ 1 / s ] Forward acceleration per speed unit change (due to thrust) M α ˙   [ 1 / s ] Pitch angular acceleration per unit rate of change of attack angle
X α   [ f t · r a d / s 2 ] Forward acceleration per unit angle of attack M q   [ 1 / s ] Pitch angular acceleration per unit pitch rate
X δ e   [ f t · r a d / s 2 ] Forward acceleration per unit elevator angle M δ e   [ 1 / s 2 ] Pitch angular acceleration per unit elevator angle
Z u   [ 1 / s ] Vertical acceleration per speed unit change
Table 2. Cessna-182 operating conditions in straight and level flight.
Table 2. Cessna-182 operating conditions in straight and level flight.
Operating Conditions
Altitude [feet]5.000
Velocity [m.p.h.]130
Weight [pounds]2.650
Dynamic Pressure [p.s.i]49.6
Center of gravity in percent [%]26.4
Table 3. Derivative coefficients of the Cessna-182 airplane.
Table 3. Derivative coefficients of the Cessna-182 airplane.
Derivative Coefficients and Values
Derivative CoefficientUnits Derivative CoefficientUnits
X u = 0.95 f t / s 2 f t / s Z δ e = 0 f t / s 2 r a d
X T u = 0.47 f t / s 2 f t / s M u = 0 r a d / s 2 f t / s
X α = 605.76 f t / s 2 r a d M T u = 0 r a d / s 2 f t / s
X δ e = 0 f t / s 2 r a d M α = 19259.390 r a d / s 2 r a d
Z u = 9.09 f t / s 2 f t / s M T α = 0 r a d / s 2 r a d / s
Z α = 1425807 f t / s 2 r a d M α ˙ = 2542.51 r a d / s 2 r a d / s
Z α ˙ = 61.63 f t / s 2 r a d / s M q = 4336.60 r a d / s 2 r a d / s
Z q = 181.38 f t / s 2 r a d / s M δ e = 35251.27 r a d / s 2 r a d
Table 4. Implementation of the IO-DMRACs and FO-DMRACs.
Table 4. Implementation of the IO-DMRACs and FO-DMRACs.
Reference model W m s = 2 s 2 + 3 s + 2
Plant W p s = 34.7012 s 2 + 74.9025 s + 4.2914 s 4 + 9.127 s 3 + 34.7561 s 2 + 1.7542 s + 0.8695
Control laws θ t = k t θ 1 T t θ 0 t θ 2 T t T R 8
ω t = r t ω 1 T t y p t ω 2 T t T R 8
u t = θ T t ω t e 1 t ω ¯ T t ω ¯ t
Note: θ 1 T t and θ 2 T t   R 3
Auxiliary signals ω ˙ 1 t = Λ ω 1 t + l u t
ω ˙ 2 t = Λ ω 2 t + l y p t
Λ = 0 1 0 0 0 1 6 11 6
l = 0     0     1 T
First order filter L(s) with a = 1. L s = 1 s + a = 1 s + 1
Output error e 1 t = k p k m W ¯ m s ϕ T t ω ¯ t
W ¯ m s = s + 1 W m s y ω ¯ t = 1 s + 1 ω t
Integer order adaptive laws θ ˙ t = γ e 1 t ω ¯ t
Fractional order adaptive laws D t α t 0 C θ t = γ e 1 t ω ¯ t
Table 5. J n o r m cost function using optimization by PSO for IO-DMRAC and FO-DMRAC.
Table 5. J n o r m cost function using optimization by PSO for IO-DMRAC and FO-DMRAC.
ISEISI I S E + I S I J n o r m
IO-DMRAC-PSO0.02073.0563.076710.31
FO-DMRAC-PSO0.00903.0053.0148.989
Table 6. Implementation details of the general IO-DMRAC controller for the case n * 2 .
Table 6. Implementation details of the general IO-DMRAC controller for the case n * 2 .
Reference Model W m s = 2 s 2 + 3 s + 2
Plant W p s = 34.7012 s 2 + 74.9025 s + 4.2914 s 4 + 9.127 s 3 + 34.7561 s 2 + 1.7542 s + 0.8695
Generalized control law θ t = k t θ 1 T t θ 0 t θ 2 T t T R 8
ω t = r t ω 1 T t y p t ω 2 T t T R 8
u t = θ T t ω t
Auxiliary signals ω ˙ 1 t = Λ ω 1 t + l u t
ω ˙ 2 t = Λ ω 2 t + l y p t
Λ = 0 1 0 0 0 1 6 11 6
l = 0     0     1 T
u ¯ t = W m s u ( t )
ω ¯ t = W m s ω ( t )
Output error and augmented error e 1 t = y P t y m ( t )
e 2 t = θ T t ω ¯ t u ¯ ( t )
ε ( t ) = e 1 t + k t e 2 t
Integer order adaptive laws k ˙ 1 t = γ ε ( t ) e 2 t 1 + ω ¯ t ω ¯ T t
θ ˙ t = γ ε ( t ) ω ¯ t 1 + ω ¯ t ω ¯ T t
Table 7. J n o r m cost function using PSO for IO-DMRAC and FO-DMRAC.
Table 7. J n o r m cost function using PSO for IO-DMRAC and FO-DMRAC.
I S E + I S I J n o r m
IO-DMRAC-PSO3.076710.31
FO-DMRAC-PSO3.0148.989
Generalized IO-DMRAC-PSO3.2904613.5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ceballos Benavides, G.E.; Duarte-Mermoud, M.A.; Orchard, M.E.; Travieso-Torres, J.C. Pitch Angle Control of an Airplane Using Fractional Order Direct Model Reference Adaptive Controllers. Fractal Fract. 2023, 7, 342. https://doi.org/10.3390/fractalfract7040342

AMA Style

Ceballos Benavides GE, Duarte-Mermoud MA, Orchard ME, Travieso-Torres JC. Pitch Angle Control of an Airplane Using Fractional Order Direct Model Reference Adaptive Controllers. Fractal and Fractional. 2023; 7(4):342. https://doi.org/10.3390/fractalfract7040342

Chicago/Turabian Style

Ceballos Benavides, Gustavo E., Manuel A. Duarte-Mermoud, Marcos E. Orchard, and Juan Carlos Travieso-Torres. 2023. "Pitch Angle Control of an Airplane Using Fractional Order Direct Model Reference Adaptive Controllers" Fractal and Fractional 7, no. 4: 342. https://doi.org/10.3390/fractalfract7040342

APA Style

Ceballos Benavides, G. E., Duarte-Mermoud, M. A., Orchard, M. E., & Travieso-Torres, J. C. (2023). Pitch Angle Control of an Airplane Using Fractional Order Direct Model Reference Adaptive Controllers. Fractal and Fractional, 7(4), 342. https://doi.org/10.3390/fractalfract7040342

Article Metrics

Back to TopTop