Journal of Control Science and Engineering

Volume 2015, Article ID 720203, 10 pages

http://dx.doi.org/10.1155/2015/720203

## Less Conservative Optimal Robust Control of a 3-DOF Helicopter

^{1}Universidade Tecnológica Federal do Paraná (UTFPR), 86812-460 Apucarana, PR, Brazil^{2}Research Laboratory in Control, Department of Electrical Engineering, Universidade Estadual Paulista (UNESP), Campus of Ilha Solteira, 15385-000 Ilha Solteira, SP, Brazil^{3}Universidade Tecnológica Federal do Paraná (UTFPR), 86300-000 Cornélio Procópio, PR, Brazil

Received 6 October 2014; Accepted 14 January 2015

Academic Editor: Yongji Wang

Copyright © 2015 L. F. S. Buzachero et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This work proposes an improved technique for design and optimization of robust controllers norm for uncertain linear systems, with state feedback, including the possibility of time-varying the uncertainty. The synthesis techniques used are based on LMIs (linear matrix inequalities) formulated on the basis of Lyapunov’s stability theory, using Finsler’s lemma. The design has used the addition of the decay rate restriction, including a parameter *γ* in the LMIs, responsible for decreasing the settling time of the feedback system. Qualitative and quantitative comparisons were made between methods of synthesis and optimization of the robust controllers norm, seeking alternatives with lower cost and better performance that meet the design restrictions. A practical application illustrates the efficiency of the proposed method with a failure purposely inserted in the system.

#### 1. Introduction

The history of linear matrix inequalities (LMIs) in the analysis of dynamical systems dates back to over 100 years. The story begins around 1890 when Lyapunov published his work introducing what is now called Lyapunov’s theory [1]. The researches and publications involving Lyapunov’s theory have grown up a lot in recent decades [2], opening a very wide range for various approaches such as robust stability analysis of linear systems [3].

In addition to the various current controllers design techniques, the design of optimal robust controllers (or controller design by quadratic stability) using LMI stands out for solving problems that previously had no known solution [4]. These designs use specialized computer packages [5].

Recent publications have found a certain conservatism inserted in the analysis of quadratic stability, which led to a search for solutions to eliminate this conservatism. Finsler’s lemma [6] has been widely used in control theory for the stability analysis by LMIs, with results similar to the LMI quadratic stability, but with extra matrices, allowing a certain relaxation in the stability analysis (referred as “extended” stability), by obtaining a larger region of feasibility. However, these methods only guarantee stability for uncertain time-invariant systems or with very small variation rate [7].

The main focus of this paper is to propose an improved method of design and optimization technique than the one presented in [8], including the possibility of time-varying the uncertainty, searching for lower gains of controllers that have the same dynamic performance when the system is forced to have a faster transient with the inclusion of decay rate restriction (a -stability performance index [9] which ensures that the eigenvalues of the uncertain system are on the left of scalar) in the LMIs formulation.

Comparisons will be made through a practical application in a Quanser’s 3-DOF helicopter [10], considering a failure during the landing treated as a time variant uncertainty, and a generic analysis involving 1000 randomly generated polytopic uncertain systems.

The notation used throughout the paper is standard. The symbol indicates transpose and indicates the inverse matrix; means that is symmetric positive definite. denotes the nonnegative integers and the real numbers.

#### 2. Robust Stability Using Decay Rate as a Performance Index

Consider the controllable uncertain linear time-invariant system described in the state space form:

This system can be described as convex combination of the polytope vertexes: with where is the number of the polytope vertexes [1].

Considering the uncertain system (2) and exiting Lyapunov theory for designing controllers, the following theorem is stated [1].

Theorem 1. *A sufficient condition for the uncertain system (2) stability guarantee subject to decay rate greater than or equal to is the existence of matrices and , such that
**
with .**When the LMIs (4) and (5) are feasible, a state feedback matrix that stabilizes the system can be found as
*

*Proof. *See [1].

*Thus, the feedback of the uncertain system presented in (1) can be done, with (4) and (5) being sufficient conditions for asymptotic stability of the polytope for a state feedback system with decay rate restriction. If the solution of LMIs is feasible, the uncertain system’s stability is guaranteed.*

*In many situations the norm of the state feedback matrix is high, precluding its practical application. In [11] an optimization method that minimizes the gains of the controller designed via LMIs (4) and (5) was proposed, but this does not eliminate the inherent conservatism.*

*Thus, in [8] a new way to optimize the norm of controller was presented eliminating the conservatism of the LMIs by Finsler’s lemma.*

*We can use the Finsler’s lemma to express the stability in terms of LMIs, with advantages over the existing Lyapunov’s theory [1], once it introduces new variables (, ) under conditions which involve only , , and .*

*Lemma 2 (Finsler). Consider , , and with and a basis for the null space of (i.e., ). Then the following conditions are equivalent: (1), (2)(3)(4).*

*2.1. Robust Stability of Systems Using Finsler’s Lemma with Decay Rate Restriction*

*2.1. Robust Stability of Systems Using Finsler’s Lemma with Decay Rate Restriction*

*Defining , , , and , note that corresponds to the feedback system with and represents the stability restriction with decay rate formulated from the quadratic Lyapunov function [1]. In this case the dimensions of Lemma’s 2 variables are and .*

*Thus, it is possible to characterize stability through the quadratic Lyapunov function (), generating new degrees of freedom for the synthesis of controllers.*

*From existing proof of Finsler’s lemma it can be concluded that the properties of one to four are equivalent. Thus, we can rewrite fourth property as follows:(4), such that
conveniently choosing the matrix variables , with nonsymmetric and a relaxation constant that has the function of flexible matrix in the LMI [13]. This constant can be defined by making a one-dimensional search. Applying the congruence transformation at the left and at the right, in the fourth property, and making , , and , the following LMIs were found:
with , , , and .*

*These LMIs meet the restrictions for the asymptotic stability of the system with state feedback. The stability resulting of the LMIs derived from Finsler’s lemma referred to as extended stability [14]. The advantage of using Finsler’s lemma formulation for robust stability analysis is the freedom of Lyapunov’s function, now defined as , , and . As depends on , the Lyapunov matrix use fits to time-invariant polytopic uncertainties, with permitted rate of variation being sufficiently small. Thus in [8] the following theorem was presented.*

*Theorem 3. In order to guarantee the stability of the uncertain system (2) subject to decay rate greater than or equal to , a sufficient condition is the existence of matrices , , and , such that
with .*

When LMIs (9) and (10) are feasible, a state feedback matrix that stabilizes the system can be given by

*Proof. *See [8].

*Thus, it can be feedback into the uncertain system with (9) and (10) sufficient conditions for asymptotic stability of the polytope.*

*2.2. Optimization of Norm Using Finsler’s Lemma*

*2.2. Optimization of Norm Using Finsler’s Lemma**In [8] there was a difficulty in applying the existing theory for the optimization of matrix norm [11] for the new structure of LMIs. This was due to the nonsymmetry of matrix for the controller synthesis, condition that was necessary to the LMI development when the controller synthesis matrix was . The solution found was using the idea of the optimization procedure for redesign presented in [15]. Thus, in [8] the new optimization method adequacy was proposed, with the minimization of a scalar obeying the relation with the Lyapunov function, to the new relaxed parameters through the following theorem.*

*Theorem 4. A constraint for the matrix norm of state feedback can be obtained, with and , being , , and , finding the minimum , , such that , . One can get the optimal value of solving the optimization problem with the LMIs:
where denotes the identity matrix of order.*

*Proof. *
See [8].

*This way of optimizing the norm of showed better results than the one presented in [11]. However, the optimization LMI, in order to be bound by the Lyapunov’s matrix , still does not have the minimal gains that would be found to meet the design requirements.*

*In order to improve the optimization performance, an alternative is presented in Section 2.4 to optimize the norm of controller reducing the conservatism of the LMIs design with a convenient manipulation.*

*2.3. Robust Stability of Systems Using Reciprocal Projection Lemma with the Decay Rate Restriction*

*2.3. Robust Stability of Systems Using Reciprocal Projection Lemma with the Decay Rate Restriction*

*In order to verify the advantages of the new formulation proposed in Section 2.4, one of the best examples of optimal design of was borrowed from [8] for comparison purposes. As in the extended stability case, the advantage of using the reciprocal projection lemma for robust stability analysis is Lyapunov’s function degree of freedom, now defined as , , and . As described before Theorem 3, the use of fits to time-invariant polytopic uncertainties, with permitted rate of variation being sufficiently small. To verify this, Theorem 5 follows.*

*Theorem 5. A sufficient condition which guarantees the stability of the uncertain system (2) is the existence of matrices , , and , such that LMIs (13) are met as follows:
with .*

When the LMIs (13) are feasible, a state feedback matrix which stabilizes the system can be given by

*Proof. *
See [8].

*Theorem 6 shows the optimization of for LMIs (13).*

*Theorem 6. A constraint for the matrix norm of state feedback is obtained, with , , and finding the minimum , , such that , being and therefore . We can get the optimal value of solving the optimization problem with the LMIs as follows:
where and denote the identity matrices of and order, respectively.*

*Proof. *See [8].

*2.4. New Formulation for Robust Stability of Systems Using Finsler’s Lemma with Decay Rate Restriction*

*2.4. New Formulation for Robust Stability of Systems Using Finsler’s Lemma with Decay Rate Restriction*

*Defining , and also consider . Check that the lemma variables (1) dimensions are and . Considering that is the matrix used to define the quadratic Lyapunov function, we will have the second propriety of Finsler’s lemma written as follows:(2) such that
which results in the necessary and sufficient condition for the system’s stabilizability, including decay rate:(2). *

*Thus, it is possible to characterize stability through the quadratic Lyapunov function (), generating further degrees of freedom for the controllers synthesis.*

*From existing Finsler’s lemma proof (Lemma 2), it can be concluded that the second and fourth properties are equivalent. Thus, we can rewrite the fourth property as follows:(4), such that
*

*Conveniently choose the matrix variables , with . Developing the fourth property, we have
*

*Thus the following LMIs subject to decay rate greater than or equal to were found:
with and , and , , and , .*

*These LMIs meet the restrictions for the asymptotic stability of the system with state feedback. It can be seen that the first principal minor of LMI (19) has the structure of the results found with the theorem of stability with decay rate. However, there is also, as stated in Finsler’s lemma, a high degree of freedom, due to the relaxation variable matrices and , without being symmetric, and for a robust stability approach, they may be polytopic: and , , and . Therefore the following theorem is proposed.*

*Theorem 7. In order to guarantee the stability of the uncertain system (2) subject to decay rate greater than or equal to a sufficient condition is the existence of matrices and , , and , such that*

When LMIs (21), (22), and (23) are feasible, a state feedback matrix that stabilizes the system can be given by

*Proof. *Assume that LMIs (21), (22), and (23) are feasible. Considering for , then we have

Knowing that generically
thus the following equivalences are true:

So, (25) can be rewritten as

and consequently

Thus (29) can be rewritten aswhere and , with , and .

*There is a great advantage in using the new formulations (21) and (22) using Finsler’s lemma compared with the old one (9) due to insertion of two polytopic matrices and , relaxing more the LMIs when compared to the old formulation where there was only the polytopic Lyapunov matrix ; moreover as described before Theorem 3, the use of fits to time-invariant polytopic uncertainties, with permitted rate of variation being sufficiently small, unlike this new formulation where the Lyapunov matrix is not polytopic, allowing variations in .*

*Below we have an alternative optimization of the norm of , which together with LMIs (21), (22), and (23) shows better results for the controller gains, as will be seen in Section 4.*

*Theorem 8. Given a constant , a constraint for the state feedback matrix norm is obtained, with , , , and by finding the minimum of , such that . We can get the minimum solving the follwoing optimization problem:
where and denote identity matrices of and order, respectively.*

*Proof. *Applying the Schur complement for the first inequality of (31) results in

Thus, from (35) we obtain

Replacing in (36) results in

From (32) we obtain

So from (37) and (38),
on which is the optimal controller associated with (19) and (20) or (21), (22), and (23).

*3. 3-DOF Helicopter*

*3. 3-DOF Helicopter*

*Consider the schematic model in Figure 2 of the 3-DOF helicopter [10] shown in Figure 1. This equipment is a patrimony of Control Research Laboratory at FEIS, UNESP. Two DC motors are mounted at the two ends of a rectangular frame and drive two propellers. A positive voltage applied to the front motor causes a positive pitch while a positive voltage applied to the back motor causes a negative pitch (pitch angle ()). A positive voltage to either motor also causes an elevation of the body (elevation angle () of the arm). If the body pitches, the thrust vectors result in a travel of the body (travel angle () of the arm) as well. The objective of this experiment is to design a control system to track and regulate the elevation and travel of the 3-DOF helicopter [10].*