Next Article in Journal
A New Method to Perform Direct Efficiency Measurement and Power Flow Analysis in Vibration Energy Harvesters
Next Article in Special Issue
Application of a Drone Magnetometer System to Military Mine Detection in the Demilitarized Zone
Previous Article in Journal
vis–NIR and XRF Data Fusion and Feature Selection to Estimate Potentially Toxic Elements in Soil
Previous Article in Special Issue
Matheuristics for Multi-UAV Routing and Recharge Station Location for Complete Area Coverage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Velocity-Free Attitude Control of Quadrotors: A Nonlinear Negative Imaginary Approach

Research School of Engineering, The Australian National University, Canberra, ACT 2601, Australia
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(7), 2387; https://doi.org/10.3390/s21072387
Submission received: 12 January 2021 / Revised: 22 March 2021 / Accepted: 22 March 2021 / Published: 30 March 2021

Abstract

:
In this paper, we propose a new approach to the attitude control of quadrotors, by which angular velocity measurements or a model-based observer reconstructing the angular velocity are not needed. The proposed approach is based on recent stability results obtained for nonlinear negative imaginary systems. In specific, through an inner-outer loop method, we establish the nonlinear negative imaginary property of the quadrotor rotational subsystem. Then, a strictly negative imaginary controller is synthesized using the nonlinear negative imaginary results. This guarantees the robust asymptotic stability of the attitude of the quadrotor in the face of modeling uncertainties and external disturbances. First simulation results underline the effectiveness of the proposed attitude control approach are presented.

1. Introduction

In recent decades, unmanned aerial vehicles (UAVs) have seen increasing interest within the research communities and industry due to their potential for numerous applications including, for instance, inspection, surveillance, data acquisition, and military applications. Their potential future applications include search and rescue, border patrol, surveillance of wildfires, surveillance of traffic and land surveys. An important type of UAVs are quadrotors, which have useful properties such as a simple structure, and low operation and manufacturing costs [1,2]. Numerous control methods have been proposed in order to tackle the quadrotor stability problem, see, e.g., [3,4,5,6].
In many of the aforementioned applications of quadrotors, designing an attitude controller with high level of performance and reliability is crucial. In the existing literature, a typical method of estimating the quadrotor’s attitude assumes that the angular velocity measurements are required. However, the velocity measurements in quadrotors can be noisy or even not available. Furthermore, using observer-based methods to reconstruct the angular velocity normally lead to inaccurate estimation of the velocity and, hence, degrade the attitude control performance [7,8,9,10]. This suggests using direct measurements of the attitude rather than indirect estimations.
Theoretically, different methods can be used to measure the attitude of the quadrotor such as: Euler angles, direction cosine matrices, and quaternions. However, the Euler angles have the advantage of being more physically sensible and easier to use, as each of the three Euler angles represents elementary rotations around the three principles axes of the quadrotor: roll, pitch, and yaw. Euler angles can be determined directly based on measurements earth’s magnetic field on the three body axes of the quadrotor. Sensors like AMR (anisotropic magneto resistive) magnetometers [11] are used to measure earth’s magnetic field.
Motivated by the forgoing, we aim in this paper to design an attitude control scheme for quadrotor systems where the angular velocity measurements or a model-based observer reconstructing the angular velocity are not required. The proposed approach rely on using the nonlinear negative imaginary systems framework, which is recently developed in [12] for nonlinear systems which are passive from the input to the derivative of the output (rather than the output as in the classical passivity theory). The nonlinear negative imaginary property of the quadrotor system will be established, and we shall employ the stability robustness results of [12] to design a velocity-free attitude controller by direct use of the Euler angles. For future work, the nonlinear negative imaginary approach can be employed or combined with other techniques as in [13,14,15] to provide the full flight control to quadrotor systems.
The rest of the paper is organized as follows: In Section 2, main related definitions and robust stability results from the linear/nonlinear negative imaginary literature are reviewed. In Section 3, we use the Euler–Lagrange dynamics of quadrotors systems to establish the nonlinear negative imaginary property of these systems. We use in Section 4, an inner-outer loop technique to design a velocity-free stabilizing controller for the attitude of quadrotors based on recently obtained nonlinear NI stability results. Finally, in Section 5, simulation results and some concluding remarks along with future directions are provided.

2. Preliminaries

Negative imaginary (NI) systems theory has been introduced in [16] for the control of flexible structures with colocated force actuators and position sensors. NI systems theory has seen significant progress in theory and application in the last decade, see for instance [17,18]. In this section, we review some of the related definitions and stability results from the NI literature in both the linear and nonlinear case.

2.1. Negative Imaginary Systems: Linear Case

We consider here the following linear time invariant (LTI) system:
x ˙ ( t ) = A x ( t ) + B u ( t ) ,
y ( t ) = C x ( t ) + D u ( t )
where the matrices A R n × n , B R n × m , C R m × n , and D R m × m . Assume that the system (1) and (2) has the m × m real-rational proper transfer function G ( s ) : = C ( s I A ) 1 B + D . The frequency domain characterization of the NI property of the above LTI system is given in the following definition.
Definition 1
([19]). A square transfer function matrix G ( s ) is called negative imaginary if the following conditions are satisfied:
1. 
G ( s ) has no pole at the origin and in [ s ] > 0 ;
2. 
For all ω > 0 , such that j ω is not a pole of G ( s ) , and j G ( j ω ) G ( j ω ) T 0 ;
3. 
If j ω 0 ; ω 0 ( 0 , ) , is a pole of G ( j ω ) , it is at most a simple pole and the residue matrix K 0 = lim s j ω 0 ( s j ω 0 ) s G ( s ) is positive semidefinite Hermitian.
A linear time invariant system of the form (1) and (2) is NI if its transfer function is NI. An equivalent time-domain definition of the NI property for the LTI system (1) and (2) is given in the following lemma.
Lemma 1.
Suppose that the system (1) and (2) (with D = 0 ) is controllable and observable. Then, G ( s ) is negative imaginary if and only if there exists matrix P as in LMI (4) such that along the trajectories of the system, the function V ( x ) = 1 2 x T P x satisfies
V ˙ ( x ( t ) ) y ˙ ( t ) u ( t ) , t 0 .
A strict notion of the negative imaginary property of the above LTI system is provided in the following definition.
Definition 2
([19]). A square transfer function matrix G ( s ) is strictly negative imaginary (SNI) if:
1. 
G ( s ) has no poles in [ s ] 0 ;
2. 
j [ G ( j ω ) G T ( j ω ) ] > 0 for ω ( 0 , ) .
The next two lemmas provide a state-space characterization of the NI and SNI properties for the LTI system (1) and (2), respectively.
Lemma 2
([20]). Let ( A , B , C , D ) be a minimal state-space realization of the transfer function matrix G ( s ) . Then G ( s ) is negative imaginary if and only if d e t ( A ) 0 , D = D T and there exist matrices P = P T > 0 , W R m × m , and L R m × n such that the following LMI is satisfied:
P A + A T P P B A T C T B T P C A ( C B + B T C T ) = L T L L T W W T L W T W 0 .
Lemma 3
([21]). Let ( A , B , C , D ) be a minimal state-space realization of the transfer function matrix G ( s ) . Then G ( s ) is strictly negative imaginary if and only if:
1. 
det ( A ) 0 , D = D T ;
2. 
there exists a matrix P = P T > 0 , P R n × n , such that
A P 1 + P 1 A T 0 , and B + A P 1 C T = 0 ;
3. 
the transfer function matrix M ( s ) A B L P A 1 0 has full column rank at s = j w for any ω ( 0 , ) where L T L = A P 1 P 1 A T . That is, rank M ( j ω ) = m for any ω ( 0 , ) .
The stability robustness of a positive feedback interconnection of NI system is established in the following theorem:
Theorem 1.
[16] Assume G ( s ) is a negative imaginary system with no poles at the origin and H ( s ) is a strictly negative imaginary system such that G ( ) H ( ) = 0 and H ( ) 0 . Then, the positive feedback interconnection of G ( s ) and H ( s ) , as in Figure 1, is internally stable if and only if
λ m a x ( G ( 0 ) H ( 0 ) ) < 1 ,
where λ m a x ( · ) denotes the maximum eigenvalue of a matrix with only real eigenvalues.

2.2. Nonlinear Negative Imaginary Systems

In [12], negative imaginary systems theory has been recently generalized to nonlinear systems. Here, we review the main definitions and stability results of nonlinear negative imaginary systems. We consider the following multi-input multi-output (MIMO) general nonlinear system of the form
x ˙ = f ( x , u )
y = h ( x )
where f : R n × R m R n is Lipschitz continuous function and h : R n R m is continuously differentiable function. The following definitions give a time-domain characterization of the nonlinear negative imaginary property of the nonlinear system (5) and (6).
Definition 3
([12]). The system (5) and (6) is said to be nonlinear negative imaginary if there exists a non-negative function V : R n R of a class C 1 such that the following dissipative inequality
V ˙ ( x ( t ) ) y ˙ T ( t ) u ( t ) ,
holds for all t 0 . Here, the function V is called a storage function.
Analogously to [22] for passive systems, we introduce a slightly stronger notions of the above definition for the purpose of stability analysis. We have the following two definitions.
Definition 4
([12]). The system (5) and (6) is said to be a marginally strictly nonlinear NI system if the dissipative inequality (7) is satisfied, and for all u and x such that
V ˙ ( x ) = y ˙ T ( t ) u ( t )
for all t > 0 , then lim t u ( t ) = c , and c is a constant vector.
Definition 5
([12]). The system (5) and (6) is said to be weakly strictly nonlinear NI system if it is marginally strictly nonlinear NI and globally asymptotically stable with u = 0 .
Remark 1.
For LTI systems of the form (1) and (2), weakly strictly nonlinear negative imaginary property becomes strictly negative imaginary property. This can be readily seen by considering a positive definite storage function as V ( x ) = 1 2 x T P x , where x is the state vector of the system and the matrix P is positive definite symmetric matrix which satisfies the LMI (4). Differentiating V with respect to time, we have V ˙ ( x ( t ) ) = y ˙ T u y ˜ T y ˜ , where y ˜ is the output of the auxiliary system given by H ( s ) = s M ( s ) = L P s I A 1 B L P A 1 B which has no zeros on the imaginary axis except at the origin (see Lemma 3). Since the system is stable, u ( t ) can consist only of exponentially decaying terms and sinusoids (including zero frequency). Note that Y ˜ ( s ) = H ( s ) U ( s ) . The condition that H ( j ω ) has no zeros on the imaginary axis implies it is not zero for all ω 0 , which guarantees the convergence of u ( t ) to a (possibly zero) constant.
In what follows, we highlight the main robust stability result introduced in [12] for the positive feedback interconnection of two nonlinear NI systems. This nonlinear stability result will be used later to robustly stabilize the attitude of the quadrotor system.
Now, consider the following general two MIMO nonlinear systems described by:
H 1 : x ˙ 1 = f 1 ( x 1 , u 1 ) y 1 = h 1 ( x 1 )
and
H 2 : x ˙ 2 = f 2 ( x 2 , u 2 ) y 2 = h 2 ( x 2 )
where h i : R n R n is a C 1 function with h i ( 0 ) = 0 , f i : R n × R n R n is continuous and locally Lipschitz in x i for bounded u i , and where f i ( 0 , 0 ) = 0 . We shall consider the open-loop interconnection of the systems H 1 and H 1 as shown in Figure 2. This interconnected system determines the stability properties of the closed-loop system, see [12]. We have the following assumptions for the open-loop interconnection of systems H 1 and H 2 .
Assumption 1.
For any constant u ¯ 1 , there exists a unique solution x ¯ 1 , y ¯ 1 to the equations
0 = f 1 ( x ¯ 1 , u ¯ 1 ) y ¯ 1 = h 1 ( x ¯ 1 )
such that u ¯ 1 0 implies x ¯ 1 0 and the mapping u ¯ 1 x ¯ 1 is continuous.
Assumption 2.
For any constant u ¯ 2 , there exists a unique solution ( x ¯ 2 , y ¯ 2 ) to the equations
0 = f 2 ( x ¯ 2 , u ¯ 2 ) y ¯ 2 = h 2 ( x ¯ 2 )
such that u ¯ 2 0 implies x ¯ 2 0 and the mapping u ¯ 2 x ¯ 2 is continuous.
Assumption 3.
h 1 T ( x ¯ 1 ) h 2 ( x ¯ 2 ) 0 , for any constant u ¯ 1 where u ¯ 2 = y ¯ 1 .
Assumption 4.
For any constant u ¯ 1 , let ( x ¯ 1 , y ¯ 1 ) be defined as in Assumption 1 and ( x ¯ 2 , y ¯ 2 ) be defined as in Assumption 2 where u ¯ 2 = y ¯ 1 . Then there exits a constant 0 < γ < 1 such that for any u ¯ 1 and with y ¯ 2 defined as in Assumption 2, the following sector bound condition:
y ¯ 2 T y ¯ 2 γ 2 u ¯ 1 T u ¯ 1 ,
holds.
The stability robustness of the positive feedback interconnections of systems H 1 and H 2 has been given in the next theorem using the Lyapunov theory and the LaSalle’s invariance principal.
Theorem 2
([12]). Consider a positive feedback interconnection of systems H 1 and H 2 where u 1 = y 2 , u 2 = y 1 . Suppose that the system H 1 is nonlinear NI and zero-state observable, and the system H 2 is weakly strictly nonlinear NI. Moreover, suppose that Assumptions 1–4 are satisfied. Then, the equilibrium point x 1 , x 2 = ( 0 , 0 ) of the closed-loop system of H 1 and H 2 is asymptotically stable.

3. Quadrotor System

In this section, we aim to reveal the nonlinear negative imaginary structure of the quadrotor system using the Euler–Lagrange equation of the quadrotor. First we recall the kinematics and dynamics model of the quadrotor with parameters as shown in Table 1.

3.1. Kinematics Model

Two reference frames are used to study the quadrotor system (see Figure 3): a reference frame fixed to the earth { R } ( O , x , y , z ) , and a body-fixed frame { R B } { O B , x B , y B , z B } , where O B is fixed to the center of mass of the quadrotor. { R B } is related to { R } by a position vector ξ = [ x , y , z ] T , describing the position of the center of gravity in { R B } relative to { R } and by a vector of three independent angles, known as Euler angles and denoted by η = [ ϕ , θ , ψ ] T , which represent roll, pitch, and yaw angles of the quadrotor. It is assumed that the Euler angles are bounded as follows:
ϕ ( π / 2 , π / 2 ) , θ ( π / 2 , π / 2 ) , ψ ( π , π ] .
A vector in the body reference frame can be transformed into vectors in the earth reference frame. For example, given a force F B , expressed using the coordinates of the body frame, the force F can be expressed in the coordinates of the earth frame as follows:
F = R B E F B
where R B E is the transformation (rotation) matrix given by
R : = C θ C ψ C ψ S θ S ϕ C ϕ S ψ C ϕ C ψ S θ + S ϕ S ψ C θ S ψ S θ S ϕ S ψ + C ϕ C ψ C ϕ S θ S ψ C ψ S ϕ S θ C θ S ϕ C θ C ϕ .
Here S ( · ) and C ( · ) represent the functions sin ( · ) and cos ( · ) , respectively. Likewise, the relation between the angular velocity vector in the inertial frame ω = [ p , q , r ] T , and the angular velocity in the body frame η = [ ϕ , θ , ψ ] T is given as:
ω = W η η ˙
where
W η : = 1 0 S θ 0 C ϕ S ϕ C θ 0 S ϕ C ϕ C θ .
Remark 2.
For small-angle approximation, we obtain an equality between the Euler rates η ˙ = [ ϕ ˙ , θ ˙ , ψ ˙ ] T and the angular velocity vector ω = [ p , q , r ] T , i.e., W η = I .

3.2. Euler–Lagrange Model: Nonlinear Negative Imaginary Structure

Here, an Euler–Lagrange approach is adopted in order to write the equations which describe the translational and rotational motion of the quadrotor. The Euler–Lagrange dynamics of the quadrotor are given by
d d t L q ˙ L q = F
where q = [ ξ T , η T ] T = [ x , y , z , ϕ , θ , ψ ] T is the generalized coordinates vector for the quadrotor, and F = [ F ξ T , τ T ] T , is the input (generalized forces) of the system where F ξ is the thrust force and τ is total torque. The Lagrangian L ( q , q ˙ ) of the quadrotor is the difference between the total kinetic energy T and the potential energy P; that is
L = T t r a n s + T r o t P .
Herein, T t r a n s represents the translational kinetic energy, and is given by
T t r a n s = 1 2 m ξ ˙ T ξ ˙ ,
where m denotes the whole mass of the quadrotor. The term T r o t represents the rotational kinetic energy, and is given by
T r o t = 1 2 η ˙ T J η ˙ ,
where J is positive-definite matrix denotes the rotational inertia matrix of the quadrotor and is defined in the body frame as
J = J x 0 0 0 J y 0 0 0 J z .
The rotational inertia matrix J is diagonal due to the quadrotor’s symmetry with the four arms aligned with the body x and y axis. The potential energy is P = m g z , where g is the acceleration due to gravity. Then, the Lagrangian of a quadrotor is given as follows:
L ( q , q ˙ ) = 1 2 m ξ ˙ T ξ ˙ + 1 2 η ˙ T J η ˙ m g z = 1 2 q ˙ T M q ˙ + G ( q ) ,
where
M = m I 3 × 3 0 3 × 3 0 3 × 3 J , G ( q ) = 0 0 m g 0 0 0 T .
To establish the nonlinear negative imaginary property of the quadrotor we express the equation (14) in terms of the generalized coordinates. Using the Lagrangian (19), Equation (14) can be written as follows
M ( q ) q ¨ + C ( q , q ˙ ) q ˙ + G ( q ) = F ,
where M ( q ) denotes inertia matrix (given in (20)) and is symmetric and positive definite. C ( q , q ˙ ) is the Coriolis and centrifugal matrix where C ( q , q ˙ ) = d d t M ( q ) 1 2 q q ˙ T M , the term G ( q ) is the gravitational vector where G ( q ) = P ( q ) q , and F is the input of the quadrotor system. We assume that the matrix C ( q , q ˙ ) is defined using the Christoffel symbols; then M ˙ 2 C ( q , q ˙ ) is skew-symmetric [23,24]. Furthermore, P ( q ) is assumed to have an absolute minimum at q = 0 .
In the following lemma, we show that the quadrotor system (21) with F as input and q as output is nonlinear negative imaginary.
Lemma 4.
Consider the system (21) with input F and output q . Then the system (21) is nonlinear negative imaginary with a positive-definite storage function given by
V ( q , q ˙ ) = 1 2 q ˙ T M ( q ) q ˙ + P ( q ) .
Proof. 
It can be easily shown that V is positive definite since where M ( q ) is positive-definite matrix, and P ( q ) is a nonnegative scalar quantity. Taking the first time derivative of the function V we obtain
d V d t ( q , q ˙ ) = q ˙ T M ( q ) q ¨ + 1 2 q ˙ T M ˙ ( q ) q ˙ + G ( q ) q ˙ = q ˙ T C ( q , q ˙ ) q ˙ G ( q ) + F + 1 2 q ˙ T M ˙ ( q ) q ˙ + G ( q ) q ˙ = q ˙ T F + 1 2 q ˙ T M ˙ ( q ) 2 C ( q , q ˙ ) q ˙ = q ˙ T F ,
which shows that the systems is nonlinear negative system from F to q . □
Remark 3.
By revealing the nonlinear negative imaginary structure of the quadrotor model, we can leverage powerful techniques from negative imaginary systems theory to be utilized in future work to achieve a broad range of control objectives for the quadrotor system.

3.3. State-Space Model

Since the stability results of [12] deal with state-space representation of nonlinear systems, we aim here to use the state-space representation of the quadrotor rotational subsystem to design quadrotor attitude system, which is the main scope of this paper. We see that the Lagrangian contains no cross-terms in the kinetic energy combining ξ ˙ and η ˙ , so the Euler–Lagrange Equation (14) partitions into two parts; that is, the translational and rotational components. The translational equation of the quadrotor is described by the following equation
m ξ ¨ + 0 0 m g = F ξ ,
with F ξ is the thrust force generated by the four rotors and is given by F ξ = i = 1 4 b ω i 2 , where ω i is ith-rotor’s speed and b is the thrust factor. The rotational subsystem describing the roll, pitch and yaw rotations of the quadrotor is described by the following equation
J ( η ) η ¨ + d d t { J ( η ) } η ˙ 1 2 η ( η ˙ T J ( η ) η ˙ ) = τ ,
or, by appropriate definition of variables,
J ( η ) η ¨ + C ( η , η ˙ ) η ˙ = τ ,
where the input vector τ = [ τ ϕ , τ θ , τ ψ ] T is the total torque in the pitch, roll, and yaw. The term J ( η ) 2 C ( η , η ˙ ) can be proved a skew-symmetric by using specific representation [24].
In matrix form, the vector τ is defined in terms of the four rotor speeds as follows,
τ = l b ω 2 2 ω 4 2 l b ω 1 2 ω 3 2 d ω 1 2 + ω 3 2 ω 2 2 ω 4 2 ,
where l is the arm length, the distance from the axis of rotation of the rotors to the center of the quadrotor, and d is the drag force. By defining the input of the quadrotor as follows
F ξ τ ϕ τ θ τ ψ = u 1 u 2 u 3 u 4 = b b b b 0 b 0 b b 0 b 0 d d d d ω 1 2 ω 2 2 ω 3 2 ω 4 2 ,
we obtain the overall quadrotor dynamic model in the following form
x ¨ = ( cos ϕ sin θ cos ψ + sin ϕ sin ψ ) · u 1 m y ¨ = ( cos ϕ sin θ sin ψ sin ϕ cos ψ ) · u 1 m z ¨ = g ( cos ϕ cos θ ) · u 1 m ϕ ¨ = θ ˙ ψ ˙ J y J z J x J r J x θ ˙ g ( u ) + l J x u 2 θ ¨ = ϕ ˙ ψ ˙ J z I x J y + J r J y ϕ ˙ g ( u ) + l J y u 3 ψ ¨ = ϕ ˙ θ ˙ J x J y J z + 1 J z u 4
Introducing the abbreviation g ( u ) = ω 1 ω 2 + ω 3 ω 4 , the system (28) can be represented in the form x ˙ = f ( x , u ) with the 12-dimensional state vector x = [ x 1 , , x 12 ] T = [ ϕ , ϕ ˙ , θ , θ ˙ , ψ , ψ ˙ , z , z ˙ , x , x ˙ , y , y ˙ ] T , and the input vector u = [ u 1 , u 2 , u 3 , u 4 ] T :
f ( x , u ) = x 2 x 4 x 6 a 1 x 4 a 2 g ( u ) + b 1 u 2 x 4 x 2 x 6 a 3 + x 2 a 4 g ( u ) + b 2 u 3 x 6 x 4 x 2 a 5 + b 3 u 4 x 8 g u 1 m cos x 1 cos x 3 x 10 u 1 m sin x 1 sin x 5 + cos x 1 sin x 3 cos x 5 x 12 u 1 m sin x 1 cos x 5 cos x 1 sin x 3 sin x 5
where
a 1 = J y J z / J x b 1 = l / J x a 2 = J r / J x b 2 = l / J y a 3 = J z J x / J y b 3 = 1 / J z a 4 = J r / J y a 5 = J x J y / J z

4. Attitude Control Design

As it can be seen from the previous section, the state-space model of the quadrotor can be divided into two subsystems, one of which, the rotational subsystem, describes the dynamics of the attitude (i.e., the angles) and the other describes the translation of the quadrotor. In this paper, we are interested in the problem of stabilizing the attitude of the quadrotor around a desired reference signal by using a nonlinear negative imaginary approach. For this purpose we confine ourselves to the rotational subsystem whose state is a restriction to the last 6 components of x representing the roll, pitch and yaw angles and their time derivatives. The rotational subsystem is then described by the differential equation
x ˙ = f α ( x , u ) = x 2 x 4 x 6 a 1 x 4 a 2 g ( u ) + b 1 u 2 x 4 x 2 x 6 a 3 + x 2 a 4 g ( u ) + b 2 u 3 x 6 x 4 x 2 a 5 + b 3 u 4 .
In recent studies, the use of a multi-loop control architecture has been proposed for a variety of quadrotor control problems; see for instance [25,26]. In this section, we propose an inner-outer loop architecture based on the nonlinear negative theory [12] to robustly stabilize the attitude, i.e., the Euler angles, around desired reference signal η d = [ ϕ d , θ d , ψ d ] T = [ x 1 d , x 3 d , x 5 d ] T , where the angular velocity measurements are not needed.

4.1. Inner-Control Loop

The inner-control loop is mainly designed due to the free motion behavior of the quadrotor. We define the following feedback control law:
τ = K p ( η η d ) + v
where v = [ v 1 , v 2 , v 3 ] T denotes the new input torque of the quadrotor, and K p = diag k p ϕ , k p θ , k p ψ is a positive diagonal matrix and the diagonal elements are used as tuning parameters. The architecture of the inner-loop can be interpreted as a proportional-only controller as seen in Figure 4 below.
The designed torque τ is then defined as follows,
τ = u 2 u 3 u 4 = k p ϕ x 1 x 1 d + v 1 k p θ x 3 x 3 d + v 2 k p ψ x 5 x 5 d + v 3
Using (32), and setting the desired reference signal η d = [ x 1 d , x 3 d , x 5 d ] T = 0 , the rotational subsystem can be put in the form x ˙ = f ˜ α ( x , u ) :
f ˜ α ( x , u ) = x 2 x 4 x 6 a 1 x 4 a 2 g ( u ) b 1 k p ϕ x 1 x 1 d + b 1 v 1 x 4 x 2 x 6 a 3 + x 2 a 4 g ( u ) b 2 k p θ x 2 x 3 d + b 2 v 2 x 6 x 4 x 2 a 5 b 3 k p ψ x 3 x 3 d + b 3 v 3 ,
or equivalently,
J ( η ) η ¨ + C ( η , η ˙ ) η ˙ = K p η + v .
The above rotational dynamical system (33) can be seen as a nonlinear negative imaginary system from the input v to the output η according to the following lemma.
Lemma 5.
Consider the quadrotor rotational subsystem (33) with v as input and η as output. Then, the system (33) is a nonlinear negative imaginary system with respect to the following positive-definite storage function,
V ( η , η ˙ ) = 1 2 η ˙ T J ( η ) η ˙ + 1 2 η T K p η .
Proof. 
It can be easily shown that V is a valid storage function; since the rotational inertia matrix J is positive definite and K p > 0 . Taking the time derivative of V we obtain
d V d t ( η , η ˙ ) = η ˙ T J ( η ) η ¨ + 1 2 η ˙ T J ˙ ( η ) η ˙ + η T K p η ˙ = η ˙ T [ C ( η , η ˙ ) η ˙ K p η + v ] + 1 2 η ˙ T J ˙ ( η ) η ˙ + η T K p η ˙ = η ˙ T v + 1 2 η ˙ T [ J ˙ ( η ) 2 C ( η , η ˙ ) ] η ˙ = η ˙ T v ,
which shows that the system (33) is nonlinear negative imaginary from the input v to the output η . □

4.2. Outer-Control Loop

We aim here to design the outer control loop of the rotational subsystem in order to get a positive-feedback closed-loop system which guarantees the asymptotic stability of the quadrotor attitude in view of Theorem 2. For simplicity, we will use the following linear MIMO integral resonant controller as the outer control loop controller,
C v ( s ) = [ s I + Γ Δ ] 1 Γ .
Here, Γ and Δ are positive-definite matrices given by Δ = diag δ , δ , δ , and Γ = diag Γ , Γ , Γ , where δ and Γ are tuning parameters. The transfer function matrix C v ( s ) is strictly negative imaginary [16]. The dc-gain (the gain of the system at steady-state) of the controller is C v ( 0 ) = Δ 1 .
In order to apply Theorem 2, we need to validate Assumptions 1–4 on the open-loop interconnection (as shown in Figure 5) of the quadrotor rotational subsystem (33) and the SNI controller (35) in the steady-state case. We have f ˜ α ( x ¯ , u ¯ ) = 0 , this implies x ¯ 2 = x ¯ 4 = x ¯ 6 = 0 and
x ¯ 1 = v ¯ 1 k p ϕ , x ¯ 3 = v ¯ 2 k p θ , x ¯ 5 = v ¯ 3 k p ψ ,
which shows for every constant value of v ¯ there is a corresponding unique value of x ¯ and hence Assumption 1 holds. Furthermore, Assumption 2 trivially holds since the controller (35) is a linear system.
We can easily see that Assumption 3 is valid since we have y ¯ c = Δ 1 u ¯ c = 1 δ y ¯ , where δ > 0 , it yields
y ¯ T y ¯ c = 1 δ y ¯ T y ¯ 0 .
Lastly, using (36) we see that
y ¯ c T y ¯ c = 1 δ 2 x ¯ 1 + x ¯ 3 + x ¯ 5 = 1 δ 2 v ¯ 1 2 k p ϕ 2 + v ¯ 2 2 k p θ 2 + v ¯ 3 2 k p ψ 2 1 δ 2 max i = ϕ , θ , ψ 1 k p i 2 v ¯ 1 2 + v ¯ 2 2 + v ¯ 3 2 = 1 δ 2 max i = ϕ , θ , ψ 1 k p i 2 v ¯ T v ¯ = γ 2 v ¯ T v ¯ ,
where
γ 2 = 1 δ 2 max i = ϕ , θ , ψ 1 k p i 2 .
Thus, to ensure that Assumption 4 is valid such that γ 2 < 1 , we choose the tuning parameter δ such that:
δ 2 > max i = ϕ , θ , ψ 1 k p i 2 .
Based upon the above arguments, we have found a lower bound on the DC-gain of the outer controller that is necessary to achieve asymptotic stability of the attitude vector around the reference signal by virtue of Theorem 2. We summarize the above result in the following theorem.
Theorem 3.
Consider the closed-loop system, as in Figure 6, of the quadrotor rotational subsystem (33) and the strictly negative imaginary controller (35). Assume that the system (33) is zero-state observable, and the condition (37) is satisfied. Then the closed-loop system is asymptotically stable.
Proof. 
The proof follows from the proof of Theorem 2 along with remark 1. □
Figure 6. Block diagram of the proposed attitude control system comprised of a positive feedback interconnection of the rotational subsystem (33) and SNI controller (35).
Figure 6. Block diagram of the proposed attitude control system comprised of a positive feedback interconnection of the rotational subsystem (33) and SNI controller (35).
Sensors 21 02387 g006
Remark 4.
The above stability result leads to a robust control system since the stability is guaranteed irrespective of the quadrotor and controller parameters so long as the condition (37) is satisfied.

5. Simulation Results

In order to verify the proposed attitude control method in this paper, we present a simulation results of the underlying quadrotor system tracking a desired reference attitude η d . The quadrotor parameters used in the simulation are given as follows: m = 0.5 kg , I x x = I y y = 4.85 × 10 3 kg · m 2 , I z z = 8.81 × 10 3 kg · m 2 , g = 9.81 m / s 2 , b = 2.92 × 10 6 Ns 2 , d = 1.12 × 10 7 Nms 2 . Using these parameters, the quadrotor rotational subsystem (33) is then implemented in MATLAB/Simulink for a simulation. The parameters of the inner-loop controller are chosen as k p ϕ = k p θ = k p ψ = 5 . Based on the result of Theorem 3, the tuning gains of the outer controller are set as Δ = diag 0.3 , 0.3 , 0.3 , and Γ = diag 160 , 160 , 160 . We assume that the initial state of the attitude vector is 30 , 20 , 10 T , where the control goal is to stabilize the quadrotor at a hovering position, i.e., η d = 0 . The obtained control result is shown in Figure 7 as a time plot of all Euler angles of the quadrotor system. The simulations show that the proposed control method asymptotically stabilize the attitude of the quadrotor to the desired reference signal.

6. Conclusions and Future Directions

The nonlinear negative imaginary systems theory has been employed to control the attitude of the quadrotor model based on the nonlinear NI property of the quadrotor rotational subsystem. The designed attitude control doesn’t rely on the angular velocity measurements through an appropriate design. In specific, an inner-outer loop architecture has been proposed to design a positive feedback control system that robustly stabilizes the quadrotor’s attitude in the face of model uncertainties and disturbances. Future direction of the proposed nonlinear negative imaginary approach is to address other quadrotor control problems (trajectory tracking, altitude control, etc.) and then complement the control applied to the quadrotor.

Author Contributions

Conceptualization, A.G.G. and I.R.P.; formal analysis, A.G.G. and I.R.P.; methodology, A.G.G. and I.R.P.; writing—original draft preparation, A.G.G.; writing—review and editing, A.G.G. and I.R.P.; supervision, I.R.P.; funding acquisition, I.R.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Australian Research Council under grant DP160101121.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Du, H.; Zhu, W.; Wen, G.; Duan, Z.; Lü, J. Distributed formation control of multiple quadrotor aircraft based on nonsmooth consensus algorithms. IEEE Trans. Cybern. 2017, 49, 342–353. [Google Scholar] [CrossRef]
  2. Özbek, N.S.; Önkol, M.; Efe, M.Ö. Feedback control strategies for quadrotor-type aerial robots: A survey. Trans. Inst. Meas. Control 2016, 38, 529–554. [Google Scholar] [CrossRef]
  3. Bouadi, H.; Bouchoucha, M.; Tadjine, M. Modelling and stabilizing control laws design based on backstepping for an UAV type-quadrotor. IFAC Proc. Vol. 2007, 40, 245–250. [Google Scholar] [CrossRef]
  4. Li, S.; Li, B.; Geng, Q. Adaptive sliding mode control for quadrotor helicopters. In Proceedings of the 33rd Chinese Control Conference, Kunming, China, 22–24 May 2014; pp. 71–76. [Google Scholar]
  5. Tran, V.P.; Santoso, F.; Garratt, M.; Anavatti, S. Distributed Artificial Neural Networks-Based Adaptive Strictly Negative Imaginary Formation Controller for Unmanned Aerial Vehicles in Time-Varying Environments. IEEE Trans. Ind. Inform. 2020, 17, 3910–3919. [Google Scholar] [CrossRef]
  6. Alexis, K.; Nikolakopoulos, G.; Tzes, A. Model predictive quadrotor control: Attitude, altitude and position experimental studies. IET Control Theory Appl. 2012, 6, 1812–1827. [Google Scholar] [CrossRef] [Green Version]
  7. Mahony, R.; Hamel, T.; Pflimlin, J.M. Nonlinear complementary filters on the special orthogonal group. IEEE Trans. Autom. Control 2008, 53, 1203–1218. [Google Scholar] [CrossRef] [Green Version]
  8. Crassidis, J.L.; Markley, F.L.; Cheng, Y. Survey of nonlinear attitude estimation methods. J. Guid. Control Dyn. 2007, 30, 12–28. [Google Scholar] [CrossRef]
  9. Hua, M.D.; Ducard, G.; Hamel, T.; Mahony, R.; Rudin, K. Implementation of a nonlinear attitude estimator for aerial robotic vehicles. IEEE Trans. Control Syst. Technol. 2013, 22, 201–213. [Google Scholar] [CrossRef] [Green Version]
  10. Chang, J.; Cieslak, J.; Dávila, J.; Zhou, J.; Zolghadri, A.; Guo, Z. A two-step approach for an enhanced quadrotor attitude estimation via IMU data. IEEE Trans. Control Syst. Technol. 2017, 26, 1140–1148. [Google Scholar] [CrossRef]
  11. Ilg, M.; Chang, B.C.; Hepner, D.; Thompson, A. A micro-controller solution for AMR magnetic sensing in flying munition systems. In Proceedings of the 2005 IEEE International Conference on Mechatronics, Taipei, Taiwan, 10–12 July 2005. [Google Scholar]
  12. Ghallab, A.G.; Mabrok, M.A.; Petersen, I.R. Extending Negative Imaginary Systems Theory to Nonlinear Systems. In Proceedings of the 2018 IEEE Conference on Decision and Control (CDC), Miami, FL, USA, 17–19 December 2018; pp. 2348–2353. [Google Scholar]
  13. Liu, Z.; Theilliol, D.; Yang, L.; He, Y.; Han, J. Interconnection and damping assignment passivity-based control design under loss of actuator effectiveness. J. Intell. Robot. Syst. 2020, 100, 29–45. [Google Scholar] [CrossRef]
  14. Guerrero-Sánchez, M.E.; Hernández-González, O.; Lozano, R.; García-Beltrán, C.D.; Valencia-Palomo, G.; López-Estrada, F.R. Energy-Based Control and LMI-Based Control for a Quadrotor Transporting a Payload. Mathematics 2019, 7, 1090. [Google Scholar] [CrossRef] [Green Version]
  15. Dong, J.; He, B. Novel fuzzy PID-type iterative learning control for quadrotor UAV. Sensors 2019, 19, 24. [Google Scholar] [CrossRef] [Green Version]
  16. Petersen, I.R.; Lanzon, A. Feedback Control of Negative Imaginary Systems. IEEE Control Syst. Mag. 2010, 30, 54–72. [Google Scholar]
  17. Petersen, I.R. Negative imaginary systems theory and applications. Annu. Rev. Control 2016, 42, 309–318. [Google Scholar] [CrossRef]
  18. Lanzon, A.; Chen, H.J. Feedback stability of negative imaginary systems. IEEE Trans. Autom. Control 2017, 62, 5620–5633. [Google Scholar] [CrossRef] [Green Version]
  19. Lanzon, A.; Petersen, I.R. Stability Robustness of a Feedback Interconnection of Systems With Negative Imaginary Frequency Response. IEEE Trans. Autom. Control 2008, 53, 1042–1046. [Google Scholar] [CrossRef] [Green Version]
  20. Mabrok, M.A.; Kallapur, A.G.; Petersen, I.R.; Lanzon, A. Stability analysis for a class of negative imaginary feedback systems including an integrator. In Proceedings of the 8th Asian Control Conference (ASCC), Kaohsiung, Taiwan, 15–18 May 2011; pp. 1481–1486. [Google Scholar]
  21. Xiong, J.; Petersen, I.R.; Lanzon, A. A negative imaginary lemma and the stability of interconnections of linear negative imaginary systems. IEEE Trans. Autom. Control 2010, 55, 2342–2347. [Google Scholar] [CrossRef]
  22. Isidori, A.; Joshi, S.; Kelkar, A. Asymptotic stability of interconnected passive non-linear systems. Int. J. Robust Nonlinear Control IFAC Affil. J. 1999, 9, 261–273. [Google Scholar] [CrossRef]
  23. Lewis, F.L.; Dawson, D.M.; Abdallah, C.T. Robot Manipulator Control: Theory and Practice; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  24. Ortega, R.; Spong, M.W. Adaptive motion control of rigid robots: A tutorial. Automatica 1989, 25, 877–888. [Google Scholar] [CrossRef]
  25. Bolandi, H.; Rezaei, M.; Mohsenipour, R.; Nemati, H.; Smailzadeh, S.M. Attitude Control of a Quadrotor with Optimized PID Controller. Intell. Control and Autom. 2013, 4, 335–342. [Google Scholar]
  26. Thanh, H.L.N.N.; Phi, N.N.; Hong, S.K. Simple nonlinear control of quadcopter for collision avoidance based on geometric approach in static environment. Int. J. Adv. Robot. Syst. 2018, 15. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Positive feedback control systems of NI systems.
Figure 1. Positive feedback control systems of NI systems.
Sensors 21 02387 g001
Figure 2. Open-loop interconnection of H 1 and H 2 .
Figure 2. Open-loop interconnection of H 1 and H 2 .
Sensors 21 02387 g002
Figure 3. Quadrotor configuration with body-fixed frame and inertial frame.
Figure 3. Quadrotor configuration with body-fixed frame and inertial frame.
Sensors 21 02387 g003
Figure 4. Inner-control loop.
Figure 4. Inner-control loop.
Sensors 21 02387 g004
Figure 5. Open-loop interconnection (in the steady-state case) of the quadrotor rotational subsystem (33) and the SNI controller (35) (where ‘c’ refers to the controller).
Figure 5. Open-loop interconnection (in the steady-state case) of the quadrotor rotational subsystem (33) and the SNI controller (35) (where ‘c’ refers to the controller).
Sensors 21 02387 g005
Figure 7. Euler angles during attitude control.
Figure 7. Euler angles during attitude control.
Sensors 21 02387 g007
Table 1. Quadrotor Parameters.
Table 1. Quadrotor Parameters.
DefinitionParameterUnit
Quadrotor massm kg
Gravitational accelerationg m / s 2
Arm lengthl m
Thrust coefficientb N · s 2 / rad 2
Drag coefficientd N · s 2 / rad 2
Roll inertia x-axis J x kg · m 2
Pitch inertia y-axis J y kg · m 2
Yaw inertia J z kg · m 2
Rotational Inertial J R kg · m 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ghallab, A.G.; Petersen, I.R. Velocity-Free Attitude Control of Quadrotors: A Nonlinear Negative Imaginary Approach. Sensors 2021, 21, 2387. https://doi.org/10.3390/s21072387

AMA Style

Ghallab AG, Petersen IR. Velocity-Free Attitude Control of Quadrotors: A Nonlinear Negative Imaginary Approach. Sensors. 2021; 21(7):2387. https://doi.org/10.3390/s21072387

Chicago/Turabian Style

Ghallab, Ahmed G., and Ian R. Petersen. 2021. "Velocity-Free Attitude Control of Quadrotors: A Nonlinear Negative Imaginary Approach" Sensors 21, no. 7: 2387. https://doi.org/10.3390/s21072387

APA Style

Ghallab, A. G., & Petersen, I. R. (2021). Velocity-Free Attitude Control of Quadrotors: A Nonlinear Negative Imaginary Approach. Sensors, 21(7), 2387. https://doi.org/10.3390/s21072387

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop