Modelling and control of nonlinear negative imaginary systems

This article investigates mathematical modelling and control of a nonlinear negative imaginary system using an illustrative benchmark physical example of a quadrotor dynamic model. A generalized methodology based on Euler–Lagrange equation is applied to obtain nonlinear negative imaginary dynamic model for the quadrotor. In this method, the Kronecker product is employed to formulate the Coriolis matrix, which is then used to construct a mathematical model of a quadrotor. This article further presents nonlinear negative imaginary systems theory-based analysis and synthesis framework to find the control solution for quadrotor’s attitude stability problem while hovering. The multi-loop control scheme is applied to the quadrotor that directly uses the Euler angles (angular positions) instead of angular velocity measurements. Numerical simulation results of this paper show that the investigated control strategy ensures the asymptotic stabilization of quadrotor attitude system model in the presence of external disturbance.

Keywords Nonlinear systems · Negative imaginary systems · Quadrotor · Modelling · Feedback stability · Multi-loop control system · Attitude control

List of symbols ⊗
Kronecker product operator R Fields of real number C Fields of Complex number R m×n Real matrix of m × n dimension C m×n Complex matrix of m × n dimension R(s) Real part of s ∈ C RH m×m ∞ Proper real stable transfer function matrix of m × m dimension A T Transpose of matrix A A * Complex-conjugate transpose of A C 1 Differentiable functions whose derivative is continuous (·) Signal constant value

Introduction
The negative imaginary (NI) systems are one type of dynamical systems where systems exhibit negative imaginary frequency response. The concept of NI systems theory was first proposed in [1,2] as a system modelling and controller synthesis framework for the dynamics of flexible structure, and subsequently, over the past decade, a lot of research works have been found in the literature towards analysis and synthesis for several dynamical systems (see the work of [3][4][5][6][7][8][9][10][11][12][13][14], etc.). The development of NI systems theory is motivated by inertial systems where second Newton's law of motion is used to describe everything. In such system, actuation is given by a torque or force input, while sensing is acquired by measuring the position or angle output. Therefore, by selecting the suitable input and output, a variety of dynamical systems may be modelled as NI system. For example, NI transfer function can be found in any flexible dynamic structure model with colocated force actuator (input) and position sensor (output) [15]. NI systems theory and positive real (PR) systems theory are closely related. However, the PR systems theory applies to the systems with relative degrees ranging from zero to one, whereas NI systems theory expanded this range to include the systems with relative degrees of zero to two.
For the analysis and control design of nonlinear systems, the classical dissipativity and passivity theories are most widely utilized. These theories characterize the dissipation of energy in the system with respect to a supply energy rate from the environment. However, many systems do not come under this classical theory framework. There are many physical systems where u is the force input (mechanical or electrical) and y is the corresponding displacement output, but the associated energy pair is (u,ẏ) instead of (u, y), and in this case, systems are dissipative (passive) from u toẏ [16]. For example, dynamical models of flexible structures having force input and position output are passive (or dissipative) from the input to the derivative of the output (instead of the output as in the classical passivity theory). More examples of such systems are discussed in [17]. For such passive (or dissipative) systems, NI systems theory is significant tool for the analysis and design of control systems. In [14], it has been shown that the negative imaginary control of flexible structure dynamic models through position sensors is effective and robust against plant uncertainty and unmodelled dynamics. This motivates us to find a nonlinear negative imaginary mathematical model of a benchmark physical system having force input and position output. In this work, a benchmark model of quadrotor system is investigated for modelling and control of nonlinear NI system.
Quadrotor is one type of unmanned aerial vehicles and there has been a lot of work published in the literature on modelling and control of quadrotor systems (see [18][19][20][21][22][23][24]) and many approaches for obtaining the dynamical model of quadrotor based on Euler-Lagrange equations are available in the literature. The skew symmetry of dynamic equations is an essential property that is frequently used to construct a variety of control methodologies. The typical existing approach for achieving the skew-symmetric property is to obtain the Coriolis matrix using Christoffel symbols of the first kind [25], which is a time-consuming and complex process. In [18,24], such Coriolis matrix based on Christoffel symbols is directly employed for the dynamic modelling of quadrotor. In this paper, a quite simple method based on Kronecker product is proposed to formulate the Coriolis matrix for quadrotor, which ensures the skew-symmetric property of quadrotor dynamical equations. In [18][19][20][21]24], the attitude of quadrotor is determined using measurements of angular velocity. However, in quadrotors, velocity measurements can be inaccurate or even unavailable. Moreover, adopting observer-based methodologies to reconstruct angular velocity results in faulty velocity measurement, which affects attitude control performance [23]. This recommends that rather than employing indirect estimations, direct measures of attitude should be used. The proposed control scheme of this paper is based on multi-loop control architecture that directly uses the angular positions (attitude), i.e. Euler angles φ (roll), θ (pitch) and ψ (yaw), of the quadrotor instead of angular velocity measurements.
In [26][27][28], nonlinear NI systems theory was recently introduced via the combination of Lyapunov and dissipativity theory techniques. The nonlinear NI property is similar to the dissipativity property with a supply rate uẏ = W (u,ẏ) for a nonlinear system with input u and output y, where y is the only function of state x. The concept of output strictly passive (OSP) systems developed a subclass of NI systems called output strictly negative imaginary (OSNI) systems (see [29,30]). The OSNI property of a system having minimum realization (A, B, C, D) with input u and output y is identical to the dissipativity property with a supply ratė [26,29]. Many existing results of NI linear systems theory were extended for nonlinear systems in [26][27][28]. In [26], an asymptotic stability theorem for closed-loop positive interconnection between a NI nonlinear system and a output strictly NI (OSNI) nonlinear system was also established under some assumptions. Along the same line of the research as discussed in [26][27][28], this work extends nonlinear NI control theory further and applies it to a benchmark unmanned aerial vehicle model. Moreover, in this work, nonlinear NI systems are defined for MIMO (multi-input multi-output) systems rather than SISO (single-input single-output) systems. In [27,28], a weak strict NI nonlinear system is investigated in positive feedback interconnection with NI nonlinear system, but in this work a OSNI nonlinear system is demonstrated as a controller, and then, asymptotic stability is proved between quadrotor (NI nonlinear) system and controller (OSNI nonlinear). In this article, the primary purpose is to apply NI nonlinear theory to provide control solution for the attitude stability problem of quadrotor while it is hovering. The numerical simulation results of this paper justify the proposed control synthesis method for the quadrotor dynamical model.
The major contribution of this article is to obtain dynamic model of quadrotor using Kronecker product-based Coriolis matrix and then to show that quadrotor dynamical model is nonlinear negative imaginary between force input and position output. The nonlinear NI theory-based multi-loop control scheme for quadrotor is proposed to obtain robust control solution for attitude stability of quadrotor while it is hovering. In particular, the contributions of this article are as follows: (a) formulation of Coriolis matrix using Kronecker product, (b) establishment of nonlinear NI dynamic model for unmanned aerial vehicle, i.e. quadrotor, and (c) nonlinear NI theory-based multi-loop control scheme for quadrotor and its asymptotic stabilization scheme. The remaining article is structured as follows: Sect. 2 illustrates some preliminaries related to vector calculus and basic notions of linear NI systems that enhance better comprehension of the article's key points. Section 3 discusses some fundamental concepts of nonlinear negative imaginary (NI) systems which are critical to this work. In Sect. 4, a nonlinear negative imaginary model of quadrotor systems is established using Euler-Lagrange equation. In Sect. 5, nonlinear state-space model of the quadrotor is obtained for control synthesis. Section 6 demonstrates the nonlinear NI approach for stabilizing the attitude of quadrotor. The numerical simulation and analysis of the results are shown in Sect. 7. Finally, Sect. 8 concludes the article with few comments.

Preliminaries
This section presents some mathematical preliminaries related to Kronecker product, time derivative of matrix and the partial derivatives of two matrices product with respect to vector. These mathematical preliminaries briefs few consistent definitions and theorems, which make article self-explanatory. Furthermore, some basic notions of linear negative imaginary systems are also introduced, which helps to streamline the concepts of this article.

Mathematical preliminaries
Definition 1 [31] If A ∈ R p×q and B ∈ R r ×s are two given matrices, then A ⊗ B ∈ R pr×qs is called Kronecker product of A with B, and it is given by Definition 2 [31] If matrix A( y) ∈ R p×q and vector y ∈ R r , the partial derivative of A( y) with respect to y is given by ∂ y 1 · · · ∂a 11 ∂ y r ∂a 12 ∂ y 1 · · · ∂a 12 ∂ y r · · · ∂a 1q ∂ y r · · · ∂a 1q ∂ y r ∂a 21 ∂ y 1 · · · ∂a 21 ∂ y r ∂a 22 The following theorems and lemma are used in Sect. 4 to establish the nonlinear dynamical model of quadrotor.
Theorem 1 [31] If A( y) ∈ R p×q and B( y) ∈ R q×s are two matrices and y ∈ R r be any vector, then the partial derivative of product of A( y) and B( y) with respect to y is defined as where I r ∈ R r ×r denotes the identity matrix.

Theorem 2
If matrix A( y) ∈ R p×q and vector y ∈ R r , the time derivative of A( y) with respect to y is obtained aṡ The proof of the Theorems 1 and 2 can be found in [31].
To verify the existence of skew-symmetric property for Coriolis matrix in quadrotor dynamic equation, the following lemma will be used (see Sect. 4).

Lemma 3
Suppose I m ∈ R m×m is the identity matrix and y ∈ R n is a vector. The two products (I m ⊗ y) y and ( y ⊗ I m ) y exist if n = m, and the following rule has been fulfilled.
Definition 3 [34] A proper real rational transfer matrix G :

Definition 4 A proper real rational function matrix
Definition 5 [26,27,30] A transfer function matrix G(s) ∈ RH m×m ∞ is called output SNI (OSNI) system, if there is a scalar value δ > 0 such that The scalar δ > 0 is the index that determines the output strictness level of OSNI system. Definition 6 [29] A dynamical system G defined by (7) is said to be dissipative with respect to an energy supply rate W (u, y) if there exists a differentiable storage function V : Theorem 4 [26,29,30] Assume that G be the LTI square, causal system with the following state equationṡ where A B C D is the minimal realization of the system with D T = D and matrix A is Hurwitz. Consider y − Du =ỹ and if G(s) be the associated the transfer function matrix, the followings are equivalent for a given scalar δ > 0: 1. G(s) is output strictly NI (OSNI) with the output strictness level δ; 2. There is a real matrix X = X T > 0 so that AX + X A T + δ (C AX) T (C AX) ≤ 0 and B = −AXC T ;

Fundamentals of nonlinear NI systems
This section illustrates some fundamental concepts of nonlinear NI systems along with significant findings, which are critical to this work.
Definition 7 [26,27] The dynamical systems given bẏ are called nonlinear NI if there exists a positive definite storage function V : R p → R of class C 1 function such that the dissipative inequality given bẏ Positive closed-loop interconnection between NI nonlinear and OSNI nonlinear system is valid for all t ≥ 0, where f : R p × R m → R p and h : R p → R m are continuous Lipschitz and class C 1 function, respectively.
Definition 8 [26] The system (8) is called nonlinear OSNI if there exists a positive definite storage function V : R p → R of class C 1 function and a scalar δ > 0 such thaṫ for all t ≥ 0.
The scalar value δ is known as output strictness level of the system (8).

Asymptotic feedback stability for nonlinear NI systems
Consider the positive feedback loop between a nonlinear NI system H 1 and a nonlinear output SNI (OSNI) system H 2 as illustrated in Fig. 1. This positive interconnection determines the asymptotic stability of feedback system under certain assumptions [26,28].
The following equations are assumed to describe the systems H i (i = 1, 2): Consider the following assumptions that satisfy Eqs. (11) and (12) for every x i (t), y i (t) and u i (t).
Assumption 1 [26] For any interval of time [t 1 , t 2 ] in which Assumption 2 [26] For any interval of time [t 1 , t 2 ] in which Assumption 3 [26] Consider the open-loop interconnection as displayed in Fig. 2. Let y 1 (t) be the output of H 1 for any is the input of H 2 and y 2 (t) ≡ȳ 2 is the corresponding constant output of H 2 , there exists a scalar 0 < γ < 1 such that for constant inputū 1 and its corresponding constant outputȳ 2 , we havē The next theorem provides the concept of feedback stability for positive interconnection as shown in Fig. 1. In other words, this theorem illustrates concept of asymptotic stability, i.e. states of the system converge to the origin.
Theorem 5 [26] Consider the positive feedback interconnection of the nonlinear NI system H 1 and nonlinear OSNI system H 2 as displayed in Fig. 1 such that Assumptions 1-3 are confirmed. If a storage function W (x 1 , x 2 ) given by is positive definite, where the functions V 1 (x 1 ) and V 2 (x 2 ) are both positive definite that satisfy (9) for H 1 and (10) for H 2 , respectively, and functions h 1 (x 1 ) and h 2 (x 2 ) are the output of H 1 and H 2 , respectively, then the positive feedback system depicted in Fig. 1 is asymptotically stable.
The concept of feedback stability for positive interconnection as shown in Fig. 1 will be further used in outer control loop of multi-loop control architecture (Fig. 4) for quadrotor system (see Sect. 6).

Dynamic model of nonlinear NI systems
Consider a quadrotor structure as shown in Fig. 3, where the angular velocities, torque and forces are generated by four rotors. Let {R} (O, x, y, z) represent the inertial reference frame fixed to the earth and The linear position and the angular position (i.e. attitude) with three Euler angles of the quadrotor are represented by ξ and η, respectively, and q is the vector that contains both ξ and η. The rotations of quadrotor around x, y and z axes are determined by the Roll angle φ, pitch angle θ and yaw angle ψ, respectively.
The Euler angles are also assumed to be constrained as follows [18]: It is possible to transform the vector in body frame to the vector in inertial frame. If the force F B is represented in body reference frame coordinates, the force F can be described as follows in inertial reference frame coordinates where R B→E is called transformation (rotation) matrix from body reference frame to the inertial reference frame and it is given by where C and S denote cos and sin terms, respectively. Since the transformation matrix R is orthogonal, the transformation matrix from inertial reference frame to the body reference frame is given by R −1 = R T . If W η is the transformation matrix for angular velocity from inertial reference frame to the body reference frame, then the relationship between the angular velocity vector [ p q r] T = in inertial reference frame and angular velocity vector [φθψ] T =η in the body reference frame can be expressed as: The structure of quadrotor having four arms aligned with x and y axes of the body is assumed to be symmetric, and therefore, corresponding diagonal inertia matrix I with I x x = I yy is given by Here I x x , I yy and I zz denote roll, pitch and yaw moment of inertia. The angular velocity, denoted by ω i of rotor i, generates f i force in rotor axis direction. The angular velocity ω i and accelerationω i of the rotor produce torque τ Mi around Lift and drag constants are denoted by k and b, respectively, whereas I M denotes the inertia moment of rotor. In most cases, the impact ofω is assumed to be minor, and hence, it is neglected. The total forces of the rotors produce T thrust in z-axis direction of the body. The torques τ φ , τ θ and τ ψ in the direction of respective body frame angles constitute the τ B torque.
The length between the quadrotor's centre of mass and rotor is given by l. The roll moment is acquired by reducing the second rotor's velocity while raising the fourth rotor's velocity. The pitch moment is achieved by reducing the first rotor's velocity and raising the third rotor's velocity. By increasing the two opposite rotor's angular velocities, and reducing the velocities of other two, yaw moment is achieved .

Euler-Lagrange equation-based nonlinear NI system
The Euler-Lagrange approach is followed to develop a dynamic model of nonlinear negative imaginary system that explains the quadrotor's translational and rotational motion. A Lagrangian can be determined by modelling the system's energy, and therefore, the kinetic and potential energy differences are considered while deriving the Lagrangian. The system's kinetic energy is modelled for both rotational and translational motion, whereas the system's potential energy is solely dependent on altitude of quadrotor. If E trans and E rot denote the kinetic energy for translational and rotational motion and the potential energy is denoted by E pot , then Lagrangian L for the quadrotor can be obtained as If J(η) denotes the positive definite Jacobian matrix from toη, Hence, the rotational energy is represented in the earth frame (inertial frame) as Thus, Lagrangian L is given by where Here N is positive definite matrix. The Euler-Lagrange equation with torques and external forces of the quadrotor are given by If we assume h = Nq, then Eq. (27) is written as Using Theorem 1 including Definition 1 and Definition 2, we can obtain ∂L ∂q from Eq. (30) as where I n ∈ R n×n denotes the identity matrix. If = h 1 e T 1 I n · · · h n e T n I n = h 1 · · · h n = h T =q T N, (32) where e i is the ith column vector of I n , anḋ then, using Eqs. (32) and (33) in Eq. (31, we obtain Similarly, we can find ∂L ∂q from Eq. (30) as where gravity and Coriolis/centrifugal terms are as follows: To extract the Coriolis matrix C from Eq. (42), we need to show that matrixṄ −2C is having skew-symmetric property. For this, we use Lemma 3 for the identity matrix I n ∈ R n×n and the vectorq ∈ R n , and obtain I n ⊗q q = q ⊗ I n q (43) By dividing the first term of right-hand side of Eq. (42) into two equal parts and substituting Eq. (43) into one of them, we get The proposed Coriolis matrix formulation is constructed as The skew-symmetric property of matrixṄ − 2C is strongly guaranteed by the Coriolis matrix C as given in (45). This can be proved, if we assume According to Theorem 2, we can writeṄ = ∂ N ∂q I n ⊗q = α Since β T − β is skew-symmetric with an arbitrary square matrix β,Ṅ − 2C(q,q) skew-symmetric matrix.

Remark 1
The matrixṄ − 2C(q,q) is skew-symmetric; this is equivalent to the fact that the fictitious internal forces do no work on the system. Now it remains to show that the dynamic model of quadrotor given by Eq. (40) is nonlinear negative imaginary system. The following lemma is used for this purpose.
Since G(q) is the scalar quantity, It is well established in [25]. This means that the so-called fictitious forces, defined by C(q,q)q, do not work on the system.
According to Definition 7, this demonstrates that the quadrotor dynamical model is nonlinear NI system from force input

Nonlinear state-space model of quadrotor
The angular and linear components are independent of one another, and they may be analysed individually. The external linear force is the total thrust of the rotors. As a result, we may formulate the following equation via Euler-Lagrange Eq. (29) Here the translational force is applied to the quadroter because of the main control input u 1 = k 4 i=1 ω 2 i in z axis direction. Equation (48) represents the translational subsystem model of quadroter. To obtain the rotational subsystem model, at first, we apply Eq. (25) to calculate the rotational inertia Jacobian matrix J(η) from toη and obtain The angular external force of the quadrotor is the torques of the rotors. Therefore, using Euler-Lagrange Eq. (29), we can get the following equation.
This can be expressed as According to Eq. (42), we can rewrite Eq. (50) as where represents the Coriolis matrix, containing the gyroscopic and centripetal terms. On evaluating Eq. (52) using MATLAB Symbolic Math Toolbox, we can obtain the Coriolis matrix C(η,η) in the following form: Consider torque τ B consists of the torques τ φ , τ θ and τ ψ is applied to the quadrotor due to control input u 2 , u 3 and u 4 in direction of the respective body frame angles. Hence, according to Eq. (23), we can write From Eqs. (51) and 54, the rotational subsystem model of quadrotor is obtained in the following form: The quadroter system described by Eqs. (48) 55 can easily be written in state-space forṁ and state vector x=

Control synthesis method
In this paper, we are interested to find the control solution for the problem of stabilizing the attitude of quadrotor subsystem (55) around specified reference angles by applying a nonlinear NI technique. For this purpose, we consider controller synthesis for the rotational subsystem of quadrotor whose states are given by pitch, roll and yaw angles and its time derivatives. The nonlinear state-space model of rotational systeṁ can be obtained from Eq. (55), where Consider the multi-loop control architecture as shown in Fig. 4. In this section, we investigate the inner-outer loop control method based on nonlinear NI systems theory to stabilize the quadrotor's attitude, i.e. Euler angles η = [x 7 , x 9 ,

Inner control loop
The inner control loop is constructed due to the free dynamic nature of quadrotor. Consider the following feedback control law is the diagonal positive matrix and it's elements k φ , k θ and k ψ are treated as tuning control parameters and ν = [ν 1 , ν 2 , ν 3 ] T is new updated input torque for the rotational quadrotor system in the inner control loop. Thus, the torque τ B is defined as follows: If we consider the desired reference attitude then the nonlinear quadrotor rotational subsystem can be written in the following form: and matrix C(η,η) can be obtained from (53). Equivalently, nonlinear quadrotor rotational system can be expressed as The rotational dynamical model of quadrotor described by Eq. (60) is a nonlinear NI system, according to the following lemma: Since it is well established in [25] thatη T J − 2C η = 0, According to Definition 7, this shows that the quadrotor rotational subsystem is nonlinear NI between force input ν and displacement output η.

Outer control loop
Consider the outer positive feedback control loop as shown in Fig. 4. The following integral resonant controller (IRC) [1,38] is selected as a outer control loop controller where and are positive definite matrix given by Here and are considered as the tuning parameters. Integral resonant controller employ position actuator, force sensor and integral feedback to generate integral force feedback. Adding a direct feedthrough to a resonant system with a colocated force actuator and position sensor yields integral resonant controller [1,38]. The use of integral feedback causes the resonant poles to be dampened. The form of controller (62) is produced by combining the direct feedthrough with the integral feedback. According to the following theorem, the integral resonant controller (62) is strictly negative imaginary.
Theorem 8 [1] The transfer function matrix (62) with positive definite matrices and is strictly negative imaginary system.
In order to apply Theorem 5 in the positive feedback outer loop as shown in Fig. 4 to achieve asymptotic attitude stability of the quadrotor, we need to verify the followings: 1. SNI controller given by (62) is equivalence of nonlinear OSNI system; 2. Validity of the Assumptions 1 to 3 for the positive feedback outer loop as shown in Fig. 4.
The following proposition from [30] is used to verify that SNI system (62) belongs to the OSNI class.

Proposition 1
The subset of SNI systems that satisfy the condition belong to the OSNI class.
If we evaluate the condition (63) for the SNI controller (62), we get Thus, outer loop controller C ν given by Eq. (62) is OSNI.
is dissipative with respect to the supply rate given byẏ T u − δ ẏ 2 = W (u,ẏ). Therefore, there exists a positive definite storage function V (x) and scalar δ > 0 which satisfieṡ Therefore, by considering the storage function Since y = x, so we geṫ For δ > 0, the storage function V = 1 2 x T x satisfies the dissipative inequality. Hence, outer loop controller (62) is equivalent to nonlinear OSNI system. Now it remains to show that the validity of Assumptions 1to 3 for the positive feedback outer loop as shown in Fig. 4. In steady-state case, we have f r (x r ,τ B ) = 0; this implies thatx 8 =x 10 =x 12 = 0 and So, for every constant value of input ν, we can get the corresponding constant value of x r . As a result, Assumption 2 is valid. Because controller (62) is a linear system and in steady-state caseȳ =x r , Assumption 1 is also trivially true.
To verify the Assumption 3, consider an open-loop connection as illustrated in Fig. 5. At steady state, DC gain of the controller (62) is given by Fig. 5, we can writē Using Eqs. (65), we can writē Therefore, to guarantee that Assumption 3 is true such that 0 < γ < 1, we select the controller tuning parameter so that According to inequality (66), we have obtained a lower bound on which is required to attain asymptotic attitude stability of the quadrotor by applying Theorem 5.

Remark 2
If the parameters of controller (62) are chosen in such a way that > 0, >0 and condition (66) are satisfied, then Theorem 5 assures the robust stability of quadrotor rotational system (57) in a multi-loop control architecture as demonstrated in Fig. 4.

Numerical simulations and results analysis
In order to show the effectiveness of investigated control synthesis method for the stability of nonlinear quadrotor rotational model, the numerical simulation results are presented. The aim of the simulation is to show the attitude stability of quadrotor while it is hovering. Therefore, the desired reference attitude η d is set to zero, i.e. φ d = θ d = ψ d = 0 • . The values of the quadrotor model parameters used in the simulation are given as follows [19]: The initial state of the attitude vector is selected as φ = 30 • , θ = −20 • , ψ = 10 • , and the control objective is to maintain stability of the quadrotor at a hovering position, i.e. φ d = θ d = ψ d = 0 • . The mathematical model of the quadrotor with multi-loop control architecture is implemented for simulation in MAT-LAB/Simulink using S-function programming, and the roll (φ), pitch (θ ) and yaw (ψ) angles responses are obtained. The angles responses are represented in radians rather than degrees. The simulations results in Fig. 6 demonstrates that the proposed synthesis method stabilizes the quadrotor's attitude to the desired hovering position. To test the robustness of attitude control, we added the step disturbance to the quadrotor's attitude system and executed simulations. The simulation results shown in Figs. 7, 8 and 9 demonstrate that the controller rejects disturbances in minimum time. There is an overshoot in responses that is due to fast response of the system.

Conclusion
In this article, modelling and control of nonlinear negative imaginary system are investigated using an illustrative physical example of quadrotor model. The article first describes the complete generalized method via Euler-Lagrange equation to obtain a nonlinear NI dynamic model for the quadrotor. Moreover, using the Kronecker product, the Coriolis matrix formulation for quadrotors is also proposed. This article further illustrates the nonlinear NI theory-based multi-loop control scheme for the attitude stability of quadrotor while it is hovering. This control scheme directly employs the Euler angles (angular positions) instead of angular velocity measurements. The numerical simulation results of this paper justify the control synthesis method for the quadrotor system. Since the nonlinear NI structure of the quadrotor model is illustrated by underlying analysis of this paper, we can accomplish a variety of control goals (e.g. trajectory tracking, altitude control and so on) for the quadrotor in future.
Funding Open access funding provided by Manipal Academy of Higher Education, Manipal.

Data Availability Not applicable.
Code Availability Not applicable.

Conflicts of interest
The authors have no conflicts of interest to declare. All co-authors have seen and agreed with the contents of the manuscript and there is no financial interest to report.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.