LMI-Based Robust Stabilization of a Class of Input-Constrained Uncertain Nonlinear Systems with Application to a Helicopter Model

*is paper is concerned with the robust stabilization of a class of continuous-time nonlinear systems, with an application to the pitch dynamics of a simple helicopter model, via an affine state-feedback control law using the linear matrix inequality (LMI) approach. *e nonlinear dynamics is subject to norm-bounded parametric uncertainties and disturbances. In addition, the problem of actuator nonlinearity is addressed by considering the saturation effect of the control law. We demonstrate first that the synthesis problem of the saturated controller is expressed in terms of bilinear matrix inequalities (BMIs). *anks to the Schur complement lemma and the matrix inversion lemma, we convert these BMIs into LMIs allowing the simultaneous computation of the two gains of the affine controller. Furthermore, we address in this work the estimation problem of the domain of attraction using the invariant set concept. *is is solved by computing the largest attractive invariant ellipsoid. Compared with previous works, the research procedure of such ellipsoidal set is achieved in a single step with a reduced number of LMI constraints and then with less conservative conditions. A portfolio of numerical results is presented. *e effectiveness and robustness of the proposed saturated controller in the stabilization of the adopted helicopter pitch model toward parametric uncertainties and disturbances are illustrated through simulation results.


Background and Literature
Review. Mechatronics, such as aircraft, spacecrafts, launch vehicles, unmanned autonomous vehicles, missiles, walking robots, robot manipulators, electronic vehicles, and unmanned aerial vehicles, has played a very important role in modern industry-related applications. Nowadays, there is an ever-increasing demand of advanced control strategies for mechatronic systems with enhanced performances.
It is known, on the one hand, that almost all existing physical and mechatronic systems unavoidably include uncertainties and disturbances due to inaccurate modeling, measurement errors, exterior conditions, or parameter variations. e presence of uncertainties may cause instability and bad performances on a controlled system. us, considerable efforts have been assigned to the robust stability and stabilization of linear and nonlinear systems with parametric uncertainties. For a recent literature, we refer the readers to [1][2][3][4][5][6][7][8]. Two types of parametric uncertainties are very often considered in systems: norm-bounded uncertainty and polytopic uncertainty. In recent years, the linear matrix inequality (LMI) technique [9] has been widely used to solve the robust control for uncertain linear and nonlinear systems with polytopic uncertain parameters and norm-bounded uncertain parameters. However, most control synthesis problems cannot be written in a LMI form. However, they are written in terms of a more general form known as a bilinear matrix inequality (BMI), which is usually not exploitable numerically to solve. For some BMI problems with simple bilinearities, the YALMIP toolbox [10] was used to solve such problems.
ere are several approaches and different relaxed synthesis conditions proposed for conducting the BMIs into LMIs, such as in [5-7, 9, 11-16], just to mention a few.
us, some research works focused on writing the saturation function of the control input as a linear convex combination of some constrained variables [18,20,27,36,47,48,[50][51][52][53]. is approach leads to sufficient conditions in terms of a very large number of LMIs that are difficult to solve numerically [27][28][29]. e common goal of almost all of these previous research works is to find the largest attractive invariant ellipsoidal set in order to provide an estimation of the domain of attraction for initial conditions.
In addition, it is well recognized that most practical control systems are inherently nonlinear. e number of available results by taking the actuator saturation nonlinearity into account in the design and analysis of nonlinear control systems is still limited. Because of the difficulty of the problem itself, most research works have been focused on particular classes of nonlinear systems to design saturated controllers [5,25,26]. Furthermore, the simultaneous presence of parametric uncertainties and actuator saturation nonlinearity in physical systems has led many authors to combine the techniques of the robust control theory and those of constrained control [5,54].

1.2.
Objective of the Paper. In this paper, we are interested in the design of a robust controller to stabilize a class of nonlinear systems, for which the zero state is not the equilibrium point. As a motivation application, we consider the pitch dynamics of a simple helicopter model [55][56][57]. In fact, such helicopter pitch model has been considered as a testbed used in order to develop new advanced control strategies. Authors in [56] approximated the pitch dynamics of the helicopter with a set of piecewise affine stochastic systems. e proposed stochastic affine feedback controller is designed for the approximated model and implemented on the main nonlinear system. Authors in [57] designed an observer-based minimum variance control for the same objective by taking into consideration norm-bounded uncertainties. Such uncertainties arise in the form of the difference between the actual nonlinear dynamics and its piecewise-affine approximation. Furthermore, authors in [55] approximated also the simplified pitch model of the helicopter via a piecewise affine system. In these works [55][56][57], the saturation effect in the actuators was not taken into consideration in the control design.
In the present work, we propose an affine state-feedback control law for the adopted class of nonlinear systems, with a particular interest to the pitch dynamics of the helicopter model, where all the states are assumed to be available for direct measurement. Our main objective is to control, via such affine feedback control law, the adopted nonlinear systems to the zero state. Moreover, we consider the problem of norm-bounded parametric uncertainties and also the problem of external disturbances. Furthermore, we take into consideration the problem of actuator saturation in the design of the robust affine state-feedback control law. us, the main role of such control law is to asymptotically bring the trajectory of the nonlinear system to the zero state even if it is subject to the parametric uncertainties, the disturbances, and the actuator saturation. In addition, in this work, the problem of searching for appropriate feedback gain matrices of the affine state-feedback control law is realized also by solving the problem of estimating the largest attractive invariant set. Our design methodology of the saturated stabilizing affine state-feedback control law is based on the framework of LMIs. We show at first that the affine statefeedback controller is designed by solving an optimization problem subject to three BMI constraints. en, by using some judicious congruence transformations and some technical lemmas, we convert these BMIs into three LMI constraints. In the end of this work, we demonstrate through simulations the effectiveness of the proposed affine statefeedback control law in the robust stabilization of the pitch dynamics of the helicopter model while the largest attractive invariant set is guaranteed.

Contributions and Innovations.
e main contributions and innovations of the work in this paper can be summarized as follows: (1) e problem of robust stabilization of a class of nonlinear systems, with an application to the pitch dynamics of a simple helicopter model, under normbounded parametric uncertainties, external disturbance, and input saturation nonlinearity, is considered using the concept of the invariant set. e main objective is to control the disturbed uncertain nonlinear system to the zero state, which is not the equilibrium point in the open loop. related to the input saturation problem, into an LMI condition is achieved using the matrix inversion lemma and the Schur complement lemma. (4) An estimation of the ellipsoidal region of attraction for the nonlinear system under study is also realized. Compared with previous works [20, 47-50, 52, 58], our research procedure of the largest invariant attractive ellipsoid is achieved in only one step with a reduced number of LMI constraints.
1.4. Structure of the Paper. e rest of this paper is organized as follows. In Section 2, we present first some technical lemmas that will be used in this work, the adopted class of the nonlinear systems, and the problem formulation. e simple helicopter model, as a motivation application, and its pitch dynamics are also described in this section. e design of the robust affine state-feedback controller under saturation, based on the LMI approach, is discussed in Section 3. Transformation of the BMIs into LMIs is also realized in this section. e problem of computation/estimation of the domain of attraction (the largest attractive ellipsoid) is addressed in Section 4. Simulation results and some comparisons are presented in Section 5. Finally, concluding remarks and future works are drawn in Section 6.

Notations.
roughout this paper, A T represents the transpose of A, the symbol (⋆) in matrix inequality denotes the symmetric term of the matrix, for example, X Y (⋆) Z � X Y Y T Z and X + (⋆) � X + X T , X > 0 ( < 0) means X is a symmetric positive (negative) definite matrix, and diag(A, B, . . . , Z) represents a diagonal matrix. Moreover, O and I denote the zero matrix and the identity matrix, respectively, with appropriate dimensions.

Preliminaries and Problem Statement
In this section, we define the class of continuous-time nonlinear systems that will be investigated in this paper. e pitch dynamics of a simple helicopter model is given as an illustrative example. First of all, we present some technical lemmas that will be used subsequently.

Some Technical Lemmas
Lemma 1 (see [41]). e Lyapunov candidate function V(t) is bounded given that the initial condition V(0) is bounded, V(t) ≥ 0 is continuous, and if the following equation is true: where μ > 0 and η > 0.
Lemma 2 (the Young relation [12]). Given constant matrices X and Y with appropriate dimensions, the following inequality holds: for all positive scalar ε.
Lemma 3 (the matrix inversion lemma [9,12]). Given invertible matrices A and B such that A ∈ R n×n and B ∈ R m×m . Moreover, given matrices C and D with appropriate dimensions: C ∈ R n×m and D ∈ R m×n . en, Lemma 4 (the Schur complement lemma [9,12]). Given matrices Q � Q T , R � R T , and S with appropriate dimensions, the following propositions are equivalent: Lemma 5 (the S-procedure lemma [9,12]). Let F 0 , . . . , F p ∈ R n×n be symmetric matrices. We consider the following conditions on F 0 , . . ., F p : If there exist scalar variables τ 1 ≥ 0, . . . , τ p ≥ 0, such that then (5) holds.

Class of Nonlinear Systems.
A general class of continuous-time nonlinear systems with a control input u and under disturbance signal w is defined by the following form: where x ∈ R n a , u ∈ R n u , and w ∈ R n w . In this work, we will consider a particular class of these nonlinear systems, for which the mathematical model is given by the following differential equation: where A ∈ R n x ×n x , B ∈ R n x ×n u , C i ∈ R n x ×1 , E ∈ R n x ×n w , and the nonlinear functions f i (x) are scalars, i.e., f i (x) ∈ R for all i � 1, . . . , n f . Moreover, we will consider that, in the nonlinear system (7), we have F(0, 0, 0) ≠ 0. is means that the state x � 0 is not the equilibrium point. Hence, in the dynamics (8), we have In addition, we will consider, without loss of generality, that Remark 1. In this work, we consider a class of nonlinear systems (8), in which the nonlinear terms f i (x), for all i � 1, . . . , n f , satisfy the constraints in (10). Our LMI-based approach for the design of a stabilizing control law u, which will be designed in the sequel, can be applied for a general class of fully nonlinear systems such as (hyper) chaotic systems as in [59,60], robotic/mechatronic systems as in [43,46], and underactuated mechanical systems [61,62], just to mention a few. us, the nonlinearity term, saying , can be assumed to satisfy the following quadratic inequality [12]: where λ is a positive scalar and H is a constant matrix having appropriate dimension. Some other research papers have considered ellipsoidal condition or one-sided Lipschitz condition on the nonlinearity h(x). For this subject, we refer the readers to, for example, [3,8,60,[63][64][65] and references therein.

Motivation
Application: e Simplified Pitch Dynamics of the Helicopter Model. As an illustrative motivation application, we will consider in this work a simplified helicopter pitch dynamics model. e schematic model of the helicopter is shown in Figure 1. We refer our readers to [55][56][57] for further details on this simplified model. e pitch dynamics of the simplified helicopter model has two degrees of freedom with only one actuator. Such model was employed to develop new control strategies for mechatronic systems [55][56][57]. e nonlinear model of the pitch dynamics of the simplified helicopter is described by the following differential equations [55][56][57]: where x 1 and x 2 represent the pitch angle (θ in Figure 1) and pitch rate, respectively; I yy is the second moment around the y-axis; m heli is the mass of the helicopter; l cgx and l cgz are displacements from the center of mass (GC in Figure 1) relative to the rotation joint B shown in Figure 1; F vM is the pitch damping; s is the variance of the moment disturbance; u is the control torque exerted by the main blade of the helicopter around the y-axis; and w is the external disturbance, which is modeled as an additive white noise representing the turbulent moments on the helicopter [55][56][57]. For simplicity, let us introduce the following change of variables: a � − (F vM /I yy ), b � 1/I yy , c � − (m heli l cgx /I yy ), d � − (m heli l cgz g/I yy ), and e � s/I yy .
en, the nonlinear differential equations in (12) can be rewritten under the particular class (8) Notice that the two nonlinear terms f 1 (x) and f 2 (x) satisfy condition (10) and also constraint (9). Hence, the zero state, x � 0, is not the equilibrium (singular) point of the pitch dynamics of the helicopter model and then of the undisturbed uncontrolled nonlinear system (8), i.e., for w � 0 and u � 0.

Remark 2.
It is worth to note that the model of the simple pitch dynamics (8) of the helicopter in Figure 1 was adopted from [55][56][57]. It is a simple model where the propeller dynamics (including motor, propeller interaction, air stream interaction, gyroscopic moments, and propeller flexing), the perturbation modelling (including wind gust modelling and payload modelling), and the cyclic/collective mechanism modelling are not considered in this study. Only the turbulent moments are taken into account as an external disturbance input, w. We emphasize that all these previous unmodelled dynamics can be lumped into a single term, saying that δ(x) can be majored according to condition (11).
For more general and basic pitch models for helicopters that reproduce basic dynamic properties of these vehicles, we can refer to [66] for typical helicopter models or [67] for more modern applications.

Problem Statement.
In the sequel, and for simplicity in developing and designing LMI stability conditions, we will consider that, in (8), n f � 2, and then the class of nonlinear systems to consider is as follows: where f(x) � f 1 (x) and g(x) � f 2 (x) are the two scalar nonlinearities satisfying conditions (9) and (10). As noted before, the zero state x � 0 is not the equilibrium point of the undisturbed uncontrolled nonlinear system (8) and then the simplified nonlinear system (13).
us, we obtain Cf(0) + Dg(0) ≠ 0. Our objective in this work is to control dynamics (13) to the zero state by designing a feedback controller, even if the dynamics is under parametric uncertainties or subject to the disturbance vector w. us, we assume that the disturbance signal w is bounded such that for some positive scalar ρ.

Complexity
Furthermore, we assume that all parameters in the nonlinear system (13) are uncertain and are with a bounded norm. e values of these parameters presented in the nonlinear dynamics (13) are considered to be the nominal ones. us, the nonlinear system (13) will be considered as the nominal system, for which (A, B) is controllable.
Hence, under parametric uncertainties, the nominal system (13) will be rewritten as follows: where ΔA, ΔB, ΔC, Δ D, and ΔE are matrices containing parametric uncertainties and satisfying the following expressions: where the matrices A . . , q 4 , and i 5 � 1, . . . , q 5 , are with a bounded norm, i.e., for all i ∈ a, b, c, d, e { } and for all j ∈ i 1 , i 2 , i 3 , i 4 , i 5 . In addition, we will consider that the uncertain nonlinear dynamics (15) is subject to an input saturation of the controller u: where u i max , for all i � 1, . . . , n u , are prescribed positive scalars and u i is the ith element of the control vector u.
e problem we are addressing in this work is to find, for the uncertain disturbed nonlinear system (15), an affine state-feedback control law where K ∈ R n u ×n x and m ∈ R n u are the feedback gains to design, for which the control law u is constrained according to condition (18). e proposed control law (19) is actually a linear state-feedback law augmented with the constant term, m, to keep the system state around an operating point, that is, the zero-equilibrium point x � 0. Let us define the following set: as the region in the state space where the feedback control law m in (19) is linear (affine) in terms of the state vector x.

Complexity
In addition, let P ∈ R n x ×n x be a positive-definite symmetric matrix and denote the Lyapunov level set as where c is a prescribed positive scalar. e level set ε(P, c) defined by (21) is an invariant set [5,52]. Moreover, it is associated with the Lyapunov function V(x) � x T Px for the closed-loop system. Consequently, ε(P, c) is an invariant set of the nonlinear system (15) with the control law (19) in the sense of the following definition.

Definition 1.
A set S ⊂ R n x is said to be invariant with respect to motion of system (14), if for all initial conditions x(0) ∈ S, the motion of the system x(t) remains in S for all t > 0.
According to Lemma 1 and by posing η � μc, we will have us, on the boundary of the set ε(P, c), we have [5,47].
Condition (22) will be used in the sequel in order to design the robust stabilization conditions of the closed-loop nonlinear system (15) with the norm-bounded parametric uncertainties (17) under the bounded external disturbance (14) and subject to the saturation (18) of the affine statefeedback control law (19) under condition (35).
As we look for stability conditions for the uncertain disturbed nonlinear system (15) under the saturated feedback control law u inside the invariant ellipsoid (21), we should have the following constraint: erefore, if for a certain Lyapunov matrix P � P T > 0 satisfying conditions (22) and (23), then taking any initial condition in the invariant ellipsoid ε(P, c), we will have u � Kx + m belongs to the set L(K, m, u max ).
In the sequel, we will develop conditions under which the ellipsoid ε(P, c) is contractively invariant and the disturbed uncertain nonlinear system (15) is robustly stable by means of the saturated affine state-feedback control law (19). Moreover, we will show how to compute the gains K and m and the Lyapunov matrix P, which determine the largest ellipsoid ε(P, c) and hence obtain an estimate of the domain of attraction.
It is worth noting that, for the case of the pitch dynamics of the helicopter model, the uncertainty matrices are Moreover, it is easy to show that these uncertainty matrices ΔA, ΔB, ΔC, ΔD, and ΔE can be rewritten like so: where and δ e should satisfy condition (17) and In the sequel of this work, without loss of generality, we will take in (16a)-(16e), q 1 � 1, q 2 � 1, q 3 � 1, q 4 � 1, and q 5 � 1. Moreover, we will consider a single control input u and also a single disturbance signal w, i.e., n u � 1 and n w � 1. We will also consider the expression of the uncertainty matrices in (24).

Design of the Robust Saturated Affine State-Feedback Controller
In this section, we develop conditions satisfying robust stabilization of the uncertain disturbed nonlinear dynamics (15) under the affine state-feedback control law (19) subject to the saturation constraint (18). us, these stability conditions must verify constraints (22) and (23). We show first that the stability conditions are expressed in terms of BMIs. Next, by means of previously provided technical lemmas, we transform these BMIs into LMIs. Finally, we present an optimization problem providing the allowable maximum values of the parametric uncertainties δ a , δ b , δ c , δ d , and δ e .

BMI-Based Stability Conditions.
ese stability conditions of the disturbed uncertain nonlinear system (15) under the saturated affine state-feedback control law defined by expression (19), constraint (18), and condition (35) are presented in the following theorem. Theorem 1.
We now consider condition (22) on the Lyapunov function V(x) and also conditions (10), (14), and (17). en, using expressions of the uncertainties matrices in (24) and relying on the Young relation (Lemma 2), we obtain the following condition: is last inequality constraint can be rewritten as follows: where Hence, we obtain the matrix inequality (25).
Furthermore, from expressions (20) and (21), rewrite x ∈ ε(P, c) as and x ∈ L(K, m, u max ) as It is worth to note that the condition ε(P, c) ⊂ L(K, m, u max ) (constraint (23)) is nothing than the implication (32) ⟹ (33). erefore, by using the S-procedure lemma, this implication condition is equivalent to the existence of a positive scalar η such that the matrix inequality (26) holds.
As noted in the beginning of the previous section, the zero state x � 0 is not the equilibrium point. en, under the affine state-feedback control law (19) and for a certain and undisturbed dynamics (that is, for the nominal dynamics (13)), the zero-equilibrium point should verify the following condition: Posing C � Cf(0) + Dg(0), then relation (34) is recast as follows: It is worth to note that, in condition (35), the term Bm + C is a column vector. en, such relation in (35) can be written in a (scalar) matrix form like so: For a small enough positive scalar ε m , equality (36) can be transformed into condition (27). is completes the proof of eorem 1.
It is worth mentioning that the three matrix inequalities (25)- (27) in eorem 1 are BMIs, which are hardly tractable numerically. us, our objective is to convert these BMIs (25)- (27) into LMIs, which is the objective of the next section.
□ Remark 3. Recall that, in the previous development of the BMI stability conditions, we have considered a single control input u and then n u � 1. Hence, we used one row matrix gain K and one scalar gain m. is choice has led to the design of only one condition (26) satisfying the saturation condition of the controller u. For a general case and then for n u > 1, K ∈ R n u ×n x , m ∈ R n u , and u max ∈ R n u , the stability condition (26) becomes where e i � (0, . . . , 0, 1

Linearization of the BMI
where Based on the Schur complement lemma, matrix inequality (38) is equivalent to where (26). We premultiply and postmultiply the BMI (26) by the matrix diag(S, 1). As a result, we obtain

Linearization of the BMI
Let us introduce a new change of variables by posing en, by applying the Schur lemma to inequality (40), we obtain Relying on the matrix inversion lemma (see Lemma 3), its follows that with en, by substituting expression (43) in inequality (41), this last inequality (41) can be rewritten like so: rough expression (44), it is straightforward to show that erefore, based on this relation, we can note the following: us, as GM � G, substituting relations (47) and (48) in inequality (45) yields Accordingly, inequality (41) is transformed into inequality (49).
Our concern now is the auxiliary condition (42). As M > 0, then through the Schur lemma, inequality (42) is equivalent to By applying again the Schur lemma on inequality (50), we obtain Hence, it follows that H > 0. us, by multiplying inequality (49) by (η − 1 ) and relying on the Schur lemma, we obtain the following matrix inequality: Since α � − u 2 max , then we can write the following: with ηM � η 0 0 1/c . erefore, by taking into account relation (53), the Schur complement lemma states that inequality (52) Hence, the BMI (26) is converted into the LMI (54), where the two parameters c and η should be fixed a priori. 8 Complexity (27). We now consider the BMI (27). It can be rewritten like so:

Linearization of the BMI
By applying the Schur lemma, this matrix inequality (55) is equivalent to the following LMI:

LMI Conditions.
According to the previous linearization results, we state the following new theorem.
Theorem 2. Assume that, for some positive parameters c, μ, η, ε m , ρ, δ a , δ b , δ c , δ d , and δ e fixed a priori, there exist a symmetric matrix S, a matrix R, a scalar m, and some scalars ε 1 , ε 2 , ε 3 , ε 4 , ε 5 , ε 6 , ε 7 , ε 8 , and ε 9 so that the following set of LMIs is feasible: where en, the nonlinear system (15) with norm-bounded parametric uncertainties (17) and bounded disturbance (14) is robustly stabilizable by means of the saturated affine statefeedback control law (19) subject to constraint (18). Moreover, the matrix gain K of the control law (19) is Proof. e proof of this theorem was already achieved in the previous sections.
As noted in eorem 2, the parameters c, μ, η, ε m , ρ, δ a , δ b , δ c , δ d , and δ e should be fixed a priori to obtain a solution of LMIs (57)- (59). us, the computation complexity for solving these three LMIs is reduced. However, in this way, the design conservatism will increase and obtaining a solution becomes more difficult. Indeed, we should test several times the feasibility of these LMIs (57)- (59) with different values of all these fixed parameters. erefore, in order to reduce the conservatism of the designed LMIs, we should consider these parameters as decision variables, which should be optimized numerically. To this end, we next present improved LMI conditions. We note that the parameter ε m can be selected as a decision variable to be minimized. Nevertheless, we will consider it as a predefined constant scalar, which will be fixed to be very small. □ Remark 4. As noted in Remark 3, for a general case where n u > 1, we obtain the BMI (37) instead of the BMI (26). e linearization of the BMI stability condition (37) leads to the following LMI:

Enhanced LMI Stability Conditions.
Our immediate concern is to convert the constant parameters c, μ, η, ρ, δ a , δ b , δ c , δ d , and δ e into decision variables. First of all, it is evident that the two parameters δ c and δ d presented in the LMI (57) are independent of other unknown variables. en, they can be considered as decision variables. However, the two main difficulties lie on the one hand in the parameters c, μ, and η and on the other hand in the three parameters δ a , δ b , δ e , and ρ. Transformation of δ a , δ e , and ρ is simple [12]. Furthermore, transformation of the parameters c and δ b requires a certain judicious mathematical manipulation. Note that the parameter c appears in both LMI (57) and LMI (58). However, transformation of the two remaining parameters μ and η is quite difficult and perhaps impossible. To simplify computation, we can pose μ � η. Hence, we state the following enhanced version of eorem 2.
It is worth to mention that eorem 3 offers a wider choice of invariant ellipsoids ε(P, c) (as c varies) for optimization and will lead to less conservative estimation of the domain of attraction.
In addition, we stress that, in eorem 3, we imposed the two free parameters μ and η in eorem 2 to be equal, i.e., μ � η � λ, in order to obtain an LMI with a single parameter λ. Nevertheless, this choice, which is the only solution to simplify computation, injects some degree of conservatism in the LMI conditions (63)- (65). Furthermore, the parameter λ is the only one that should be fixed a priori to obtain a solution of the optimization problem in eorem 3.

Estimation of the Domain of Attraction
With all the ellipsoids satisfying the set invariance condition (21), we would like to choose from among them the largest one to get at least an estimation of the domain of attraction. Certain criteria are used in the literature to maximize the size of the attractive invariant ellipsoid ε(P, c), such as maximization of its volume and maximization of the sum of its semiaxis [5,9,18,21,27,47,48,52]. Actually, maximization of the volume of an ellipsoid is equivalent to minimization of the determinant of the Lyapunov matrix P, and maximization of the sum of its semiaxis is equivalent to minimization of the trace of the matrix P [9].
In this part, our main goal is to find a robust saturated affine state-feedback control law (18) associated to the largest invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ). Such largest invariant ellipsoid will be found for desired values of the maximum bounds of the parametric uncertainties δ a , δ b , δ c , δ d , and δ e and also the maximum bound of the disturbance ρ. us, we will use LMI (63) (or systematically LMI (57)) for predefined values of these maximum bounds. However, we should also use the two LMIs (64) and (65) to look for a maximum value of the parameter c and hence the associated gains K and m.
Next, the size of the attractive invariant ellipsoid ε(P, c) is measured through the sum of its semiaxis. For this subject, we introduce the following theorem.

Theorem 4.
e largest attractive invariant ellipsoid of the nonlinear system (15) with bounded disturbance (14) and norm-bounded parametric uncertainties (17), subject to the saturated affine state-feedback control law defined by expression (19) and constraint (18), is the ellipsoid ε(P, c), for a fixed positive scalar ε m ≪ 1 and some positive parameter λ known a priori, where P � S − 1 and c, together with ε 1 , ε 2 , ε 3 , ε 4 , ε 5 , ε 6 , ε 7 , ε 8 , and ε 9 , are a solution of the following LMIbased optimization problem: Proof. In expression (21), the attractive invariant ellipsoid E(P, c) can be recast as We emphasize that the trace of the matrix P/c is equal to the sum of the semiaxis of the invariant ellipsoid ε(P/c, 1) for some fixed c. Hence, the size (the length of the semiaxis) of ε(P/c, 1) increases as the trace of the matrix P/c decreases. As P � S − 1 , P/c � (cS) − 1 . erefore, maximization of the size of ε((cS) − 1 , 1) is equivalent to the maximization of the trace of the matrix cS. As a result, to guarantee the largest attractive invariant ellipsoid ε((cS) − 1 , 1), we should maximize the parameter c and also the trace of the matrix S. Notice that maximizing c and trace(S) is equivalent to minimizing (− trace(S) − c). is completes the proof of eorem 4.
It is worth to note that, for each fixed value of the parameter λ, the optimization problem in eorem 4 provides a solution associated to an attractive invariant ellipsoid ε(P, c) with a (local) largest size, which is obtained for such λ. erefore, in order to obtain the global maximum, we should vary the parameter λ, 0 < λ < ∞. To this end, the best manner is the use of the gridding method [68]. is method consists in making a new change of variable by defining κ � λ/(1 + λ) and then λ � κ/(1 − κ). We know that λ > 0 if and only if κ ∈ ]0, 1[. en, we assign a uniform subdivision of the interval ]0, 1[ and we solve the optimization problem in eorem 4 for each value of this subdivision. Once this interval is covered, we look for the maximum value between all the obtained results. is maximum value corresponds to the largest invariant ellipsoid, which was obtained for some κ ⋆ and then for some optimal values of the parameter λ, saying λ ⋆ .
In fact, finding the optimal value λ ⋆ and then the largest invariant attractive ellipsoid ε ⋆ takes a lot of time. Indeed, suppose that we uniformly subdivide the interval ]0, 1[, in which the parameter κ varies, and then we take N κ values.
us, for each value of κ, we run the optimization algorithm in eorem 4. Let us assume that a single run takes t κ . Hence, the total simulation time in order to obtain the largest ellipsoid ε ⋆ is about t c ≈ N κ × t κ . □ Remark 5. In several research works presented in the literature, the computation of the largest attractive invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ) is achieved mainly in two steps: (1) Hunt for an invariant ellipsoid ε(P, 1) by obtaining the matrix P (2) With such P, maximize the value of the parameter c Actually, these two steps are realized via two optimization algorithms. Compared with this hunting procedure, our proposed design of the saturated feedback control law guarantees the determination of the largest invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ) in only one step and then with only one optimization algorithm.

Simulation Results
In this section, we provide some simulation results to demonstrate the effectiveness of our developed method for the synthesis of the saturated affine state-feedback control law (19) for the robust stabilization of the disturbed uncertain nonlinear dynamics (15) of the pitch model of the helicopter. Moreover, we aim at identifying and analyzing the largest attractive invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ) in order to obtain an estimation of the domain of attraction. us, in the sequel, we take two values of the saturation level in (18) of the control law u: u max � 1 and u max � 10. Moreover, we fix the parameter ε m � 10 − 5 . e values of the parameters figured in the nonlinear dynamics (12) or of the pitch model of the helicopter are as follows [55][56][57]: I yy � 0.0283(kg/m 2 ), m heli � 0.9941(Kg), l cgx � 0.0134(m), l cgz � 0.0289(m), F vM � 0.0041 (Nm/rad/s), g � 9.81(m/s 2 ), and s � 0.0057(m/s 2 ).

Remark 6.
e choice of the parameters of the maximum value of the control input u, u max , is based on two facts: (1) in the literature, the common choice of the saturation limit is u max � 1 according to the invariance sets' concept and (2) in the present work, we have considered a general case of the actuator saturation limit in the development of LMI constraints and hence in the design of the state-feedback controller in order to show the effect of the saturation limit u max on the largeness of the invariant attractive ellipsoid.
us, in our investigation, we selected two values of the saturation limit: u max � 1 (which is the classical choice) and u max � 10. We can also choose another value of such saturation limit such as u max � 24, and we will obtain similar results.
Actually, the saturation limits are the intrinsic requirement of the physical actuators.
us, the adopted different values of the u max will be used in the sequel just to show the feasibility and the conservativeness of the proposed design method of the two gains of the affine state-feedback controller and also for the computation of the largest domain of attraction in different situations.

Remark 7.
It is worth to mention that for the computation of the solutions of the LMIs in eorem 2 and the LMI-based optimization problem in eorems 3 and 4, we use the LMI toolbox of MATLAB. We can also use the toolbox YALMIP [10] in order to solve these LMI problems. Moreover, we can use it to (locally) solve some BMI problems using the bisection method.

Numerical Results.
e optimization problem subject to LMI constraints in eorem 3 provides numerical results illustrated in Table 1 for the saturation level u max � 1 and for different values of the fixed parameter λ. In this table, we provide the two gains, K and m, of the saturated affine statefeedback control law u, the maximum bounds of the parametric uncertainties, i.e., δ a , δ b , δ c , δ d , and δ e , the maximum bound of the disturbance w, i.e., ρ, and the parameter c. It is obvious that the three uncertainty bounds, δ a , δ b , and δ e , are very small, around 10 − 4 . e two parameters, δ c and δ d , are found to be equal and decrease slightly as λ increases. However, the value of the allowable disturbance bound ρ is found to be very high, about 6 × 10 3 .
We emphasize that the obtained results for the allowable maximum bounds of the parametric uncertainties, δ a , δ b , and δ e , are dissatisfying (they are found to be very small). Moreover, the maximum disturbance bound is found to be very high. It is worth noting that the optimization algorithm in eorem 3 is characterized by an objective function that depends on 10 parameters ε 4 , ε 5 , ε 6 , ε 9 , α, β, φ, δ c , δ d , and ρ, which are optimized together. en, in order to improve the obtained results, we will fix in the next the parameter ρ and the remaining ones will be optimized. As discussed in the beginning of this paper, in the second section, the external disturbance w can represent the turbulent moments on the helicopter. Without loss of generality, we will fix ρ � 10.
us, for ρ � 10, the optimization algorithm subject to LMI constraints in eorem 3 gives numerical results illustrated in Table 2 for the case u max � 1 and in Table 3 for u max � 10. As in Table 1, we provide, for different values of the parameter λ, the gain matrix K, the gain constant m, the maximum bounds of the parametric uncertainties δ a , δ b , δ c , δ d , and δ e , and the parameter c.
We note from Tables 2 and 3 that the optimization problem provides identical values for the two uncertainty bounds δ c and δ d . Moreover, for all values of λ, the value of the gain m is about 0.0133 as in Table 1. In addition, it is evident that, for the same value of parameter λ, we have almost the same gain matrix K. Furthermore, as λ increases, the size of the gain matrix K and the uncertainty bounds, δ a , δ b , and δ e , increase. For u max � 1, δ c is found to be around 10 − 3 , whereas for u max � 10, δ c is about 14. In addition, it is evident that the values of the maximum bounds of the parametric uncertainties δ a , δ b , and δ e are bigger than those obtained in Table 1 since the parameter ρ is fixed here.
We can also observe from Tables 2 and 3 that, as the parameter λ increases, the value of the parameter c decreases. Actually, we can note that c ≈ u 2 max /λ. is shows that the size of the attractive invariant ellipsoid, which depends on c, varies with respect to the value of the saturation level u max .

Computation of the Largest Attractive Invariant Ellipsoid.
To obtain the largest attractive invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ), we solve the optimization problem in eorem 4. us, in order to obtain the global largest ellipsoid ε ⋆ (P ⋆ , c ⋆ ), we should vary the free parameter λ according to the gridding method as noted before. Moreover, in this optimization problem, we fix the bounds of the parametric uncertainties and we take a common value, δ, for all these bounds like so: for u max � 1, we take δ a � · · · � δ e � δ � 1, whereas for u max � 10, we fix δ a � · · · � δ e � δ � 10. Actually, with δ � 10 and u max � 1, the optimization problem in eorem 4 is unfeasible. It was found unfeasible also for all δ ≥ 2 and u max � 1.
en, we minimized the uncertainty bound to obtain a feasible solution. is will be demonstrated and discussed next. Figure 2 reveals evolution of the trace of the matrix cS, or systematically the sum of the semiaxis of the invariant ellipsoid ε (S − 1 , c), as the parameter λ varies. It is obvious from Figure 2 that, by sweeping through λ, we obtain the maximum trace of the matrix cS for the two cases u max � 1 and u max � 10. Such maximum trace corresponds to the largest invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ). Obtained results are summarized in Table 4. In this table, we provide the parameters, P ⋆ and c ⋆ , of the largest attractive invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ), the associated feedback gains, K ⋆ and m ⋆ , of the saturated affine state-feedback control law u, the nominal value, λ ⋆ , of the free parameter λ, and the size of the attractive invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ). According to Figure 2 and Table 4, in the first case u max � 1, the size of the largest invariant ellipsoid is about 1.7906. However, in the second case u max � 10, such ellipsoid has a size about 46.5870. We recall that, for u max � 1, we have chosen δ � 1 and for u max � 10, we have adopted δ � 10.
Our objective now is to see the effect of the maximum bound of the different parametric uncertainties δ a , δ b , δ c , δ d , and δ e and also the saturation limit u max on the size of the attractive invariant ellipsoid ε(P, c). To do this, we used the optimization algorithm in eorem 4 and we varied the parameter λ as mentioned previously. is task was achieved for a single value of the maximum uncertainty bound, saying δ max � δ a � · · · � δ e . us, we obtained the largest ellipsoid ε ⋆ (P ⋆ , c ⋆ ). Now, we will repeat the same procedure by also varying the parameter δ max and fixing the saturation level u max . Figure 3 reveals obtained results for u max � 1 and u max � 10. We stress first that, for u max � 1, the allowable maximum bound of the parametric uncertainty allowing to obtain an invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ) is δ max � 1.06. However, for u max � 10, this value of δ max increases considerably and reaches the limit 12.05.  Table 2: Numerical simulation results using eorem 3 for different values of the parameter λ and for ρ � 10 and u max � 1.  Moreover, it is obvious that, for each value of u max , the size of the largest ellipsoid increases rapidly as δ max decreases. A tiny decrease of the uncertainty bound δ max leads to a very important increase in the size of the ellipsoid ε ⋆ (P ⋆ , c ⋆ ). For example, for u max � 1 and for δ max � 0.5, the size is about 10 3 and for δ max � 0.1, the size is around 10 9 . Similarly, for the saturation limit u max � 10 and for δ max � 0.5, the size is around 10 8 and for δ max � 1.5, the size is around 10 5 .
Furthermore, it is clear that the size of the largest invariant ellipsoid decreases as δ max increases. e ellipsoid becomes smaller as the uncertainty bound δ max approaches its feasible limit (i.e., the value 1.06 for u max � 1 and 12.05 for u max � 10).
In addition, we remark that, for the same value of uncertainty bound δ max , the size of the largest ellipsoid ε ⋆ (P ⋆ , c ⋆ ) increases as the saturation limit u max increases (as noted previously for δ max � 0.5).
Using numerical results illustrated in Table 4, we plotted in Figure 4 the corresponding largest attractive invariant ellipsoids ε ⋆ (P ⋆ , c ⋆ ) for the two cases u max � 1 and u max � 10. In Figure 4, the two inclined lines (colored in pink and green) correspond to the bounds of the actuator saturation (20), which are defined as Kx + m � ±u max . e lower line (colored in pink) is depicted for Kx + m � u max .
It is obvious from Figure 4 that, in the two cases, the largest invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ) is contained in the region L(K, m, u max ), where |Kx + m| ≤ u max . erefore, the constraints on the affine state-feedback control law given in (18) are well respected.
In the end of this part, we can note the two following statements: (1) For a constant saturation limit u max , an increase of the maximum bound of the parametric uncertainties leads to the decrease of the largeness of the attractive invariant ellipsoid. (2) For a constant maximum bound of the uncertainty, an increase of the saturation limit u max causes the attractive invariant ellipsoid to become larger.

Remark 8.
In some research papers, the definition of the attractive invariant ellipsoid was chosen to be as follows: is set � ε(P) is equivalent to that adopted in the present work in (21), that is, ε(P, 1). Hence, we have c � 1. Our immediate concern is to show that the adopted ellipsoid (21) with a free parameter c gives less conservative results Here, for u max � 1 (resp. u max � 10), we have adopted δ � 1 (resp. δ � 10). Moreover, in the two cases, we have fixed ρ � 10.   Figure 3: Variation of the size of the largest attractive invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ) with respect to the variation of the allowable maximum bound of the parametric uncertainties for the two cases: u max � 1 and u max � 10. Here, we have fixed δ a � · · · � δ e � δ max , with always ρ � 10.
14 Complexity compared with the ellipsoid (72). We use the same optimization algorithm in eorem 4 where in the present case, we fix the parameter c as a constant one and then we take c � 1. Table 5 shows the numerical results obtained for a fixed parameter c � 1. Notice that the symbol (!) in Table 5 means that the LMI-based optimization problem in eorem 4 is unfeasible. To realize a comparison from conservatism point of view, we adopt the same parameters in Table 4. As noted previously, the maximum bound of the disturbance is fixed to ρ � 10. It is clear from Table 5 that, for the first case u max � 1 and δ � 1, the optimization problem in eorem 4 was found to be unfeasible, with c � 1. It was found to be feasible for all δ ≤ 0.648. Recall that, for a free parameter c, the optimization problem was found to be feasible for all δ < 2. is fact shows that, using a free parameter c in the definition of the attractive invariant ellipsoid ε(P, c), that is, the set (21), instead of fixing c � 1, gives less conservative stability conditions. Moreover, from Table 5 and for the case u max � 10 and δ � 10, the optimization problem with c � 1 is found to be feasible, as in Table 4. Nevertheless, the largest ellipsoid ε ⋆ ′ (P ⋆ ) obtained for c � 1, is too small compared with that obtained using a free parameter c, that is, the ellipsoid ε ⋆ (P ⋆ , c ⋆ ). Recall that the largeness of the attractive invariant ellipsoid is measured by the "trace" function as noted in Tables 4 and 5. As noted in Table 4, for a free parameter c, the volume of the largest ellipsoid was found to be about 46.5870. However, for a fixed parameter c � 1, the volume of the largest ellipsoid ε ⋆ ′ (P ⋆ ) is about 7.3840. e difference between the two sizes of the ellipsoid is evident. Another attractive result that can be observed from Table 5 is that the obtained matrix gain K ⋆ is too large compared with that obtained in Table 4.
Accordingly, we emphasize that the parameter c in the definition of the attractive ellipsoid (21) leads to less conservative results and then less restrictive LMI stability conditions. In addition, it contributes in obtaining a largest ellipsoid with a controller gain having a small size. Hence, the choice of c ≠ 0 reduces the conservatism of the LMIbased optimization problem in eorem 4.

Remark 9.
We noted in the end of Section 4, just before Remark 5, that the computation time for solving the LMIbased optimization problem in eorem 4 and then for the identification of the largest invariant attractive ellipsoid ε ⋆ (P ⋆ , c ⋆ ) depends on the number N κ of uniformly distributed points κ in the interval ]0, 1[ and the simulation time t κ for each value of the parameter κ. e total simulation time t c can be computed to be about t c ≈ N κ × t κ . Tables 6 and 7 present the simulation results for different values of the parameters N κ for u max � 1 ( Table 6) and for u max � 10 ( Table 7). Recall that, for u max � 1, we have δ � 1, whereas for u max � 10, we fixed δ � 10. Moreover, the maximum bound ρ of the disturbance is fixed to be ρ � 10 in the two cases. Actually, in order to obtain these results in Tables 6 and 7, we have fixed the interval in which the parameter κ varies as follows: [0.01: δ κ : 0.99], where we have selected four cases: δ κ � 0.01, δ κ � 0.001, δ κ � 0.0005, and δ κ � 0.0002. For these values of δ κ , it corresponds, respectively, the number N κ � 99, N κ � 981, N κ � 1961, and N κ � 4901. us, using these intervals and these values of the parameters N κ , we simulated the optimization problem under LMI constraints in eorem 4 and then we obtained the results in Tables 6 and 7. In these Tables, we give the two optimal gains, K ⋆ and m ⋆ , of the controller, λ ⋆ (notice that λ � κ/ (1 − κ)), c ⋆ , the size of the ellipsoid (trace(c ⋆ S ⋆ )), the simulation/computation time t c , and the number of feasible solutions N fs . It is important to note that not all the N κ -cases are feasible. From the results in Tables 6  and 7, it reveals that the number of feasible solutions is only about 4.8% for u max � 1 and about 6% for u max � 10. e computation time t c for finding the largest ellipsoid increases significantly as N κ increases too. We also stress that as N κ increases, the optimal solution converges to that already obtained in Table 4. We found the same values for N κ � 4901. Nevertheless, in this case, the computation time is very high, it is computed to be about 4.5 minutes for u max � 1 and it is about 2.2 minutes for the second case u max � 10.

Simulations for the Robust Stabilization of the Pitch
Dynamics. In this part, we would show the robustness and effectiveness of the designed saturated affine state-feedback control law (19) in the robust stabilization of the uncertain disturbed nonlinear dynamics (15) of the pitch model of the helicopter. e control gains K and m are provided in Table 4. We recall that, for u max � 1, we have Moreover, we recall that the external disturbance w in the nonlinear dynamics (15) is bounded according to condition (14) with ρ � 10. In addition, the uncertain parameters δ a , δ b , δ c , δ d , and δ e in (15) vary randomly over time so that the boundedness condition (17) is satisfied. Furthermore, for the simulation of the controlled nonlinear system (15), we will take first an initial condition, x 0 , located outside the largest attractive ellipsoid ε ⋆ (P ⋆ , c ⋆ ): for the case u max � 1, we take x 0 � 0.1 − 1 T , whereas for the case u max � 10, we take Each initial condition is located near the boundary of the corresponding ellipsoid (see Figure 4). Furthermore, we will take an initial condition located outside the invariant attractive ellipsoid: x 0 � 5 − 10 T , for the two saturation limits.
First of all, we will study the dynamics of the helicopter model under the saturated control law u in the nominal case, that is, without parametric uncertainties and disturbance. Figure 5(a) shows temporal evolution of the states x 1 and x 2 of the certain undisturbed controlled nonlinear system (15), and Figure 5(b) reveals the saturated affine state-feedback control law u for the saturation limit u max � 1. It is obvious that the state of the pitch dynamics converges to zero, which is our main objective in this paper. Moreover, the control law u converges to a constant value, which is found to be equal to the gain m. Actually, such constant value of the control input will ensure that the pitch dynamics of the helicopter stays at the zero state. In addition, we note that the reached maximum value of the control law u is about 0.23, which is less than the desired saturation limit u max � 1.
We have also analyzed the nominal pitch dynamics of the helicopter model for the case u max � 10. We have observed almost the same behavior in Figure 5. e reached   16 Complexity maximum value of the controller u is around 5. e controller converges to the constant value m ≈ 0.01336. Now, we will show simulation results for the controlled nonlinear system subject to both parametric uncertainties and external disturbance. Simulation results are shown in Figure 6 for u max � 1 and Figure 7 for u max � 10. Figure 6(a) (resp. Figure 7(a)) presents temporal evolution of the two states, x 1 and x 2 , of the nonlinear system (15) under uncertainties and disturbance. Figure 6(b) (resp. Figure 7(b)) reveals evolution of the saturated control law u for u max � 1 (resp. u max � 10). As noted previously, the initial condition is located inside the largest attractive invariant ellipsoid. Moreover, the disturbance w is injected into the system at t � 2[s] and during 3[s]. We recall that the external disturbance w satisfies condition (14) with ρ � 10.
It is obvious from Figure 6(a) that the motion of the pitch dynamics of the helicopter model experiences some fluctuations around the desired zero state when the disturbance is applied. Once the effect of the external disturbing torque w vanishes, the controlled system stabilizes again around its desired position with some very weak perturbations provoked by the injected parametric uncertainties. Moreover, in Figure 6(b), the control signal u varies around the value m � 0.0133. e control signal varies between two very small values: ±0.1. We stress that the effects of the parametric uncertainties and the external disturbance are compensated.
is shows the robustness of the control law u towards the parametric uncertainties and the external disturbance with high amplitude.
However, for u max � 10, the effect of the parametric uncertainties and the external disturbance w is clear in Figure 7. Figure 7(a) is almost identical to Figure 6(a). Moreover, when the disturbance is injected, the control signal u undergoes some fluctuations, which vary between ±1.5.
is happens because the maximum bound of the uncertainties is very large (δ a,...,e � 10) in this case compared with the first case, i.e., for u max � 1.
It is worth to note that, in the previous three cases, the saturation level ±u max of the controller was well respected.
We have tested another case by taking an initial condition located outside the largest attractive invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ). We take here only the case u max � 1, and then the corresponding largest ellipsoid is given by Figure 4(a). As noted previously, the initial condition is chosen to be x 0 � 5 − 10 T . Obtained results are illustrated in Figure 8. We observe first from Figure 8(a) that the state x 1 of the helicopter model experiences some smooth oscillations before its convergence and stabilization at the zero position despite the presence of the external disturbance and the parametric uncertainties. In fact, only weak fluctuations are observed as in Figure 7(a). However, the interesting phenomenon observed here is the saturation of the control law u depicted in Figure 8(b).

Remark 10.
It is worth to note that, for a predefined set of the system parameters, the previous established optimization problems under LMI constraints are (should be) simulated offline in order to compute the largest attractive ellipsoid and then the associated feedback gains K and m of the input-saturated affine state-feedback controller u before its application into a real helicopter system in practice.
us, for a prescribed saturation level u max , the possible maximum allowable bounds of the parametric uncertainties δ a , δ b , δ c , δ d , and δ e and the possible maximum bound of the external disturbance signal w, i.e., ρ, we look, as described previously, for the largest attractive invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ). us, once this set is found, the corresponding feedback gains K and m will be used in order to stabilize the pitch dynamics of the simple helicopter model in a real-world application.  Figure 5: Evolution of the two states, x 1 and x 2 , of the nonlinear system (15) and the saturated affine state-feedback control law u for u max � 1 in the nominal case. Here, the system is without uncertainties (δ a,b,c,d,e � 0) and disturbance (w � 0).

Conclusion and Future Works
In this paper, an LMI-based approach for designing a robust affine state-feedback control law to stabilize the pitch dynamics of a helicopter model was proposed. e nonlinear dynamics of the helicopter was subject to an external disturbance and norm-bounded parametric uncertainties. Moreover, the problem of the actuator saturation in the design of the control law was as well addressed. We showed that the stabilization problem is represented as a solving problem of BMI constraints. Furthermore, with a judicious utilization of the Schur lemma and the matrix inversion lemma, these BMIs were transformed into LMIs. We have also developed an optimization problem with enhanced LMI constraints permitting to compute the maximum bounds of the parametric uncertainties.
In addition, we have proposed an LMI-based approach for the maximization of the attractive invariant ellipsoid for the uncertain disturbed nonlinear dynamics of the helicopter model under the saturated affine state-feedback control law.  Figure 6: Temporal evolution of the states x 1 and x 2 of the nonlinear system (15) and the saturated affine state-feedback control law u for the saturation level u max � 1. Here, the system is subject to randomly time-varying uncertainties and also a randomly time-varying external disturbing torque w. Moreover, δ a,b,c,d,e � 1.  e largeness of the ellipsoid was measured by means of the length of its semiaxis. We have showed that the optimization problem computes efficiently the largest ellipsoid in only one step.
Finally, we have showed through numerical simulations the performance of the synthesized saturated controller in the robust stabilization of the nonlinear dynamics under the norm-bounded parametric uncertainties and the external disturbing torque.
In the literature, the estimation of the domain of attraction of the closed-loop system subject to actuation saturations is one of important issues and several methods have been achieved for both linear and nonlinear systems. e difference between these methods lies in the approach by which the saturation nonlinearity was handled [51]. Its treatment can be classified into three main approaches: (1) the first one is to treat the actuator saturation problem as achieved in the present work, (2) the second approach is to treat it as a locally/generalized sector bounded nonlinearity [17,69,70], (3) while the third method is to represent the saturation nonlinearity as a (polytopic or linear) differential inclusion (see [36,51] and references insides). As noted in Remark 5, the solution presented in this work leads to the computation of the largest attractive invariant ellipsoid in only one step, compared to some related approaches that used a two-step method. us, the approach presented in this paper is more simpler and is less restrictive.
In addition, it worth to note that, in this paper, a saturated affine state-feedback controller was designed. To the best of our knowledge, such problem has not been considered in the literature. Generally, a saturated linear feedback controller has been considered to stabilize linear and nonlinear systems. en, in order to realize a comparison and hence to show the possible efficiency of our design approach, a future work of this paper is the design of a saturated affine feedback controller using the sector-bounded nonlinearity method and the differential inclusion approach.
Furthermore, we aim at extending the present methodology of the saturation affine feedback controller for more complex nonlinear systems with different Lipschitzian conditions [65], and with measurement delays [71,72] and also for the design of observer-based feedback controllers for Lipschitz nonlinear systems [3,8,64]. Moreover, we hope to extend this work for impulsive hybrid nonlinear systems, such as the biped robots [73] and the impact mechanical oscillators subject to multiple rigid constraints [74][75][76].
In the present work, we considered the simple pitch dynamics of a helicopter model as an application. Moreover, the zero state is not the equilibrium point of such model. us, we have designed an affine feedback controller under saturation to achieve the robust stabilization at the zero state. Another important application that has such feature is the robot manipulators [77]. Indeed, because of the gravitational matrix, the dynamics of robot manipulators has an equilibrium point different to the zero state. us, several control approaches have been adopted for this subject [77]. Our objective is then to extend the design method of the saturated affine state-feedback controller for the case of manipulator robots.

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

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.  Figure 8: Temporal evolution of the two states x 1 and x 2 of the uncertain disturbed nonlinear system (15) and the saturated affine statefeedback control law u for the saturation level u max � 1. Here, the initial condition is selected to be outside the largest attractive invariant ellipsoid ε ⋆ (P ⋆ , c ⋆ ) in Figure 4(a).