## Services on Demand

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Sba: Controle & Automação Sociedade Brasileira de Automatica

*Print version* ISSN 0103-1759

### Sba Controle & Automação vol.23 no.4 Campinas July/Aug. 2012

#### http://dx.doi.org/10.1590/S0103-17592012000400001

**CONTROLE ROBUSTO**

**A new approach for modeling and control of nonlinear systems via norm-bounded linear differential inclusions**

**Uma nova abordagem para modelagem e controle de sistemas não-lineares via inclusões diferenciais lineares limitadas por norma**

**Roman Kuiava ^{I}; Rodrigo A. Ramos^{II}; Hemanshu R. Pota^{III}**

^{I}Federal University of Parana (UFPR), Polytechnic Center, Department of Electrical Engineering, Rua Cel. Francisco Heraclito dos Santos, 100, Jardim das Americas, 81531-980, Curitiba, Brazil. kuiava@eletrica.ufpr.br

^{II}University of Sao Paulo (USP), Engineering School of Sao Carlos (EESC), Department of Electrical Engineering, Av. Trabalhador Saocarlense, 400, Centro, 13566-590, Sao Carlos, Brazil. ramos@sc.usp.br

^{III}University of New South Wales at Australian Defence Force Academy (UNSW@ADFA), School of Information Technology and Electrical Engineering, Northcott Drive, Canberra, Australia. h.pota@adfa.edu.au

**ABSTRACT**

A systematic approach to model nonlinear systems using norm-bounded linear differential inclusions (NLDIs) is proposed in this paper. The resulting NLDI model is suitable for the application of linear control design techniques and, therefore, it is possible to fulfill certain specifications for the underlying nonlinear system, within an operating region of interest in the state-space, using a linear controller designed for this NLDI model. Hence, a procedure to design a dynamic output feedback controller for the NLDI model is also proposed in this paper. One of the main contributions of the proposed modeling and control approach is the use of the mean-value theorem to represent the nonlinear system by a linear parameter-varying model, which is then mapped into a polytopic linear differential inclusion (PLDI) within the region of interest. To avoid the combinatorial problem that is inherent of polytopic models for medium- and large-sized systems, the PLDI is transformed into an NLDI, and the whole process is carried out ensuring that all trajectories of the underlying nonlinear system are also trajectories of the resulting NLDI within the operating region of interest. Furthermore, it is also possible to choose a particular structure for the NLDI parameters to reduce the conservatism in the representation of the nonlinear system by the NLDI model, and this feature is also one important contribution of this paper. Once the NLDI representation of the nonlinear system is obtained, the paper proposes the application of a linear control design method to this representation. The design is based on quadratic Lyapunov functions and formulated as search problem over a set of bilinear matrix inequalities (BMIs), which is solved using a two-step separation procedure that maps the BMIs into a set of corresponding linear matrix inequalities. Two numerical examples are given to demonstrate the effectiveness of the proposed approach.

**Keywords:** nonlinear systems; linear dynamic output feedback control; linear differential inclusions; robust control; linear matrix inequalities.

**RESUMO**

Este artigo propõe um método sistematico para modelagem de sistemas não-lineares na forma de inclusões diferenciais limitadas por norma (IDLNs). O modelo de IDLN resultante é adequado para aplicação de técnicas de projeto de controle linear, o que possibilita atender critérios específicos de desempenho dinâmico para o sistema não linear original em uma região de operação de interesse no espaço de estados, a partir de um controlador linear projetado para a sua representação na forma de IDLN. Então, um procedimento para projeto de um controlador por realimentação dinâmica de saída para um sistema descrito na forma de IDLN é também proposto neste artigo. Uma das principais contribuições da abordagem proposta de modelagem e controle é a aplicação do teorema do valor intermediário para representar sistemas não-lineares na forma de um modelo linear com parâmetros variantes no tempo, o qual é então mapeado em uma inclusão diferencial linear politópica (IDLP). Para evitar o problema combinatório inerente aos modelos politópicos para sistemas de médio e grande porte, a IDLP é transformada em uma IDLN, e este processo é feito de tal forma que todas as trajetórias do sistema não-linear original sejam também trajetórias do modelo resultante de IDLN. Além do mais, é também possível escolher uma estrutura particular para os parâmetros da IDLN de forma à reduzir o conservadorismo na representação do sistema não-linear pelo modelo de IDLN, e esta característica é também uma importante contribuição deste artigo. Quanto ao projeto do controlador, ele é formulado como um problema de busca por uma solução que satisfaça um conjunto de restrições escritas na forma de desigualdades matriciais bilineares (ou BMIs, do inglês bilinear matrix inequalities). Tal solução é então encontrada usando-se um procedimento de separação em duas etapas que transforma o conjunto original de BMIs em um conjunto correspondente de desigualdades matriciais lineares (ou LMIs, do inglês linear matrix inequalities). Dois exemplos numéricos são apresentados para demonstrar a eficiência da abordagem proposta.

**Palavras-chave:** sistemas não-lineares; controle por realimentação dinâmica de saída; inclusões diferenciais lineares; controle robusto; desigualdades matriciais lineares.

**1 INTRODUCTION**

Stability analysis and control synthesis for complex systems (such as, nonlinear systems involving saturation or time-varying uncertainties) can be simplified by the use of linear descriptions of these dynamical systems in the form of linear differential inclusions (LDIs). Generally speaking, it may be possible to ensure that the trajectories of such a system exhibit certain features by analyzing a corresponding model in the form of an LDI. Necessary conditions that guarantee, for example, the existence of a polytopic LDI (PLDI) containing all trajectories of a particular nonlinear system are discussed and demonstrated in (Boyd et al., 1994; Hu and Chen, 2007).

A wide variety of analysis and control problems have been formulated and solved for different classes of LDIs, most of them in the form of linear matrix inequalities (LMIs) (e.g., (Boyd et al., 1994); (Xie and de Souza, 1992; Bernard et al., 1997)). The LMI approach makes it easy to include stability and performance specifications (such as, a minimum decay rate and bounds on the output peak values) in the formulation of analysis and synthesis problems for LDIs. Once the problem is written in terms of LMIs, efficient convex search or optimization methods can be used to find a solution to it (Boyd et al., 1994). Based on this feature, it may be easier to design a controller for a nonlinear system using a linear description of it in the form of an LDI.

The modeling, analysis and control of nonlinear systems via linear models are not restricted to approaches based on LDIs. In this sense, Mamdani and Takagi Sugeno fuzzy systems play also an important role (Montagner et al., 2010; Mozelli et al., 2010; Tognetti and Oliveira, 2010). These type of fuzzy systems allow nonlinear systems to be approximated by means of an averaged sum of linear models. Then, the problems of analysis and synthesis can be written in terms of LMIs (Montagner et al., 2010; Mozelli et al., 2010). Also, based on feedback linearization technique, adaptative neural networks or fuzzy control schemes have been introduced to approximate nonlinear systems into linear models (Chen et al., 1996).

This paper addresses both modeling and control of nonlinear systems via a particular class of LDIs: the class of the norm-bounded LDIs (NLDIs). The main idea behind the modelling technique is to represent the nonlinear system by a linear parameter-varying (LPV) system using the mean value theorem (Vidyasagar, 1993; Zemouche et al., 2005; Hossain et al., 2009). This LPV system can be particularly represented by a PLDI, provided that a certain set of conditions (that will be presented later in the paper) are satisfied.

It is possible to use LMIs to check the quadratic stability of a PLDI, but this requires the solution of one LMI for each of the vertices of the polytopic domain. One major drawback of this approach is the fact that, in general, the number of vertices of this polytope is high (see the numerical examples in (Hu and Chen, 2007; Hu, 2007)), and this translates into significant (and sometimes intractable) computational burden.

For some applications - especially the ones involving high dimensional problems (e.g. (Ramos et al., 2004; de Oliveira et al., 2009; Hossain et al., 2009)) - it may be easier to use a system description in the form of an NLDI, because checking its quadratic stability involve only two LMI constraints (Boyd et al., 1994). In this case, an efficient outer approximation (or overbounding) of a PLDI can be given by an NLDI (Boyd et al., 1994), which means that every trajectory of the former is also a trajectory of the latter. The formulation proposed in (Boyd et al., 1994) is based on the solution of an LMI feasibility problem and considers only the case where a specific matrix parameter of the NLDI is square and invertible. This condition is substituted here by a weaker one in which this matrix parameter must only have full column rank (and this is another important contribution of this paper). Using this weaker condition, it is possible to choose a particular structure of the NLDI parameters in order to reduce the conservatism in the representation of the underlying nonlinear system by the NLDI model.

Based on this approach for modeling a nonlinear system via an NLDI, this paper proposes the design of a linear dynamic output feedback (LDOF) controller for the NLDI model, in such a way that stability and performance specifications can be satisfied for the closed-loop nonlinear system using a linear controller. The motivation for focusing on this type of control comes from problems in the area of power systems stability (see, for example, (Basler and Schaefer, 2008; Ramos et al., 2004; Hossain et al., 2009)), in which linear controllers are often applied to a system with highly nonlinear behaviors, and these controllers must satisfy stability and performance specifications over a wide range of operating conditions. The control problem is formulated in this paper using quadratic Lyapunov functions and constraints in the form of bilinear matrix inequalities (BMIs).

Computational methods for solving BMIs are still under development and, in general, the existing ones are restricted to small dimensional problems (Polanski, 1997; Yfoulis and Shorten, 2004). However, it is possible to transform the particular set of BMIs resulting from the approach proposed in this paper into a set of LMIs, using a separation procedure presented in (de Oliveira et al., 2000). With this technique, it becomes possible to design a linear controller for the NLDI model of the original nonlinear system and, therefore, ensure that the closed-loop nonlinear system fulfills the stability and performance requirements within the whole operating region of interest.

This paper is organized as follows. Section 2 provides the problem formulation. Section 3 deals with the modeling problem, and this section is divided into five parts. The first part gives the basis for rewriting a nonlinear system as an LPV system; the second one defines a PLDI using the previously obtained LPV system; the third part presents a result that enables the overbounding of the PLDI by the NLDI model; the fourth part discusses the estimation of regions of attraction of nonlinear systems using NLDIs models and, the last part presents a numerical example that illustrate the application of the modeling procedure. Section 4 describes the fundamentals of the proposed controller design method, based on the NLDI model that was previously obtained. Section 6 presents some tests of the proposed control design procedure and its corresponding results. Finally, section 7 contains the conclusions and some final remarks on the proposed approach.

**Notation:** the notation used throughout this paper is standard. ^{n} denotes the *n*-dimensional Euclidean space, ^{n×m} is the set of *n×m* real matrices. The closed convex polytope defined by a finite number of vertices, say *S*_{1}, *S*_{2}, ..., *S _{v}*, is defined as the convex hull of those elements and is denoted simply as

*Co*(

*S*

_{1},

*S*

_{2}, ...,

*S*). For two elements

_{v}*a*and

*b*in

^{n}, {

*a,b*} denotes the set constituted by only these two elements, while [

*a,b*] denotes the set containing all the points in the line segment between

*a*and

*b*. For matrices and vectors ()' means transposition. For a symmetric matrix

*P, P*0 (

*P*0) denotes positive (negative) definitess. Positive (negative) semi-definiteness is denoted by

*P*0 (

*P*0). An identity matrix with appropriate dimensions is denoted simply by

*I*. For a matrix, ||·|| denotes the largest singular value of the matrix. For singular matrices, (·)

^{+}denotes the pseudo-inverse of the matrix.

**2 PROBLEM FORMULATION**

Consider a continuous-time nonlinear system described in the state-space form by

where *x*(*t*) = [*x*_{1}(*t*) ... *x _{n}*(

*t*)]' ∈

*X*⊂

^{n}is the state vector,

*u*(

*t*) ∈

*U*⊂ is the control input,

*B*is a constant matrix of proper dimension and

*f:X*

^{n}is a nonlinear function of class

*C*

^{1}. Here, the subset

*X*of

^{n}represents a state-space region of interest of (1) given by (Rohr et al., 2009)

where *a _{k}* ∈

^{n}are given constant row vectors and

*n*is the number of edges of

_{e}*X*.

This form of system (1), in which the nonlinearities are present only in the dynamics of the state with respect to itself, is a peculiarity of some practical systems (see Ramos et al., 2004; de Oliveira et al., 2009)). For convinience, we assumed *x* = 0 (with *u* = 0) as being the equilibrium point of interest, so *f*(0) = 0, 0 ∈ *X* and 0 ∈ *U*. This assumption is quite standard and can be satisfied by a simple change of coordinates (Vidyasagar, 1993).

Initially, our goal is to model the nonlinear system (1) in the form of a linear differential inclusion defined as

where *D*(*x*(*t*),*u*(*t*)) := {*z*(*t*) : *z*(*t*) = Λ*x*(*t*) + *Bu*(*t*), ∀ Λ ∈ Ω}, being Ω a non-empty, closed set of real (*n×n*)-matrices (Boyd et al., 1994; Pyatnitskiy and Rapoport, 1996). Depending on the form of the set Ω, different types of LDIs can be obtained. In this paper, we are interested in the class of norm-bounded LDIs, in which the set Ω has the particular form given by

where *A*_{0} ∈ ^{n×n}*, F* ∈ and G ∈ are known fixed matrices, while E is any real (*n _{p}×n_{q}*)-matrix satisfying ||

*E*||

__<__1. As is already well-known (Boyd et al, 1994), this NLDI (3)-(4) is equivalent to the linear time-varying system

where Δ*A*(*t*) = *FE*(*t*)*G*, being *E*(*t*) an unknown matrix satisfying ||*E*|| __<__ 1 for all *t* > 0.

Consider again the nonlinear system (1). The mean value theorem guarantees the existence of a matrix *J*(*t*) such that *f*(*x*(*t*)) = *J*(*t*)*x*(*t*), for every *x*(*t*) ∈ *X*. Hence, by an adequate choice of matrices *A*_{0}, *F* and *G* of (5), it can be possible to guarantee that *J*(*t*) ∈ Ω_{NLDI} for all *t* > 0, which follows immediately that every trajectory of the nonlinear system (1) is also a trajectory of the NLDI (5).

In this paper, *A*_{0} is assumed to be the Jacobian matrix obtained by truncating the Taylor series expansion of (1) (with *u* = 0) at the first-order term, so the local properties of (1) are well described by (5). On the other hand, by a proper choice of matrices *F* and *G*, the nonlinear behavior of (1) within the set *X* is expected to be well described by the term Δ*A*(*t*). Section 3 discusses the main fundamentals of the proposed procedure to calculate these two matrices of the model.

Once we have a model of (1) in the form of the NLDI (5), the problem of interest is the determination of a linear dynamic output feedback controller in the state-space form

that stabilizes the system (5) and guarantee a desirable performance to the controlled system, where *x _{c}*(

*t*) ∈ and

*y*(

*t*) =

*Cx*(

*t*) is the measured output of (1). This will be discussed in section 4.

**3 MODELING A NONLINEAR SYSTEM VIA NLDI**

The procedure proposed in this paper to calculate matrices *F* and *G* is based on the description of (1) in the form of an LPV system, within the state-space region of interest specified by *X* in the previous section. This is done using the mean value theorem (Vidyasagar, 1993; Zemouche et al., 2005), as it is discussed in section 3.1. Using this reformulation of (1), section 3.2 shows that it is possible to define a politopic LDI whose set of trajectories contains all the solutions of the LPV system. Finally, in section 3.3 we calculate the matrices *F* and *G* by solving an optimization problem in the form of LMIs which guarantees that Ω_{NLDI} ⊇ Ω_{PLDI}, where Ω* _{PLDI}* is the particular form of Ω associated to the politopic LDI (Boyd et al., 1994; Pyatnitskiy and Rapoport, 1996).

**3.1 Rewriting the nonlinear system as an LPV system**

We present the version of the mean value theorem that is applicable to the general case where *f*(*x*) = [*f*_{1}(*x*) ... *f _{q}*(

*x*)]', with

*f*:

_{i }^{n },

*i*= 1,...,

*q*. Consider the canonical basis of the vectorial space

^{s}, for

*s*

__>__1, given by

*E*= {

_{s}*e*(

_{s}*i*) :

*e*(

_{s}*i*) = (0,...,0,1,0,...,0)',

*i*= 1,...,

*s*}. Using the canonical basis

*E*of the vectorial space

_{q}^{q}, it is possible to write

*f*(

*x*) as

*f*(

*x*) =

*e*(

_{q}*i*)

*f*(

_{i}*x*). Now, we can state the following proposition (Zemoucge et al., 2005).

**Proposition 1** Let *f*(*x*) : ^{n } ^{q}. Let *a* and *b* two elements in ^{n}, and assume that *f* is differentiable on *Co*(*a,b*). Then, there are constant vectors *c _{i}* ∈

*Co*(

*a,b*),

*c*≠

_{i}*a, c*≠

_{i}*b, i*= 1,...,

*q*, such that

For the proof, see (Zemouche et al., 2005). In this paper, we use this version of the mean value theorem to rewrite the nonlinear system (1) as an LPV system. We first assume that *f* is differentiable on the line segment between *x*(*t*) and the equilibrium point at the origin *x* = 0, *i.e.*, on the set *Co*(*x*(*t*),0) = { λ*x*(*t*) : λ ∈ [0,1]}, ∀*t* > 0. Hence, Proposition 1 guarantees the existence of *n* vectors (*t*) ∈ *Co*(*x*(*t*),0), *x _{si}*(

*t*) ≠

*x*(

*t*), (

*t*) ≠ 0,

*i*= 1,...,

*n*, such that

where *Q _{ij} = e_{n}*(

*i*)

*e*(

_{n}*j*)'. Now, let us define the functions

where *i,j* = 1,...,*n* and the second term in the right side of the equality (9) is the element *i×j* of the Jacobian matrix obtained by truncating the Taylor series expansion of (1) at the first-order term. This Jacobian matrix will be denoted simply by matrix *A*_{0} and it describes the linear properties of (1) in the neighbourhood of the origin. So, we can write (8) as being

Observe that the second term on right side of the equality (10) captures only the nonlinear behavior of (1). It is important to point out that the mean value theorem guarantees the existence of the vectors (*t*), *i* = 1,...,*n*, for which the equality (10) is satisfied, but it does not provide any means to calculate them. We can deal with this problem, however, by setting a bounded range for the functions *h _{ij}*(

*t*),

*i,j*= 1,...,

*n*. To do so, let us assume that,

This allows us to specify a lower bound *h** _{ij}* and an upper bound

*by*

_{ij}We can use the fact that the functions *h _{ij}* are bounded (and their bounds are given by (12)) to obtain a description of the nonlinear system (1) within the region

*X*via the following LPV system with bounded parameters

where θ* _{ij }*: ℜ

_{+}[

*h**,*

_{ij}_{ij}],

*i,j*= 1,...,

*n*and

*x*(0) =

*x*

_{0}.

Observe that the LPV system (13) captures the effects of the nonlinearities of (1) on the system dynamics by the time varying-parameters θ* _{ij}*(

*t*),

*i,j*= 1,...,

*n*. It is important to emphasize, however, that this description of (1) may be quite conservative, since the exact relationship between θ

_{ij}(

*t*) and

*h*(

_{ij}*t*) (given by (9)) is neglected in (13). As a consequence, many trajectories of the LPV system (13) may not be trajectories of (1). This is the price we pay for representing the nonlinear system (1) in a linear form described by (13).

In relation to the practical meaning of the assumption (11) and the set *X*, notice that most of the physical systems operate within a bounded range of their variables (such as voltages, frequencies and rotor angles, for an electric power system (Basler and Schaefer, 2008; Ramos et al., 2004; Hossain et al., 2009)), which means that it is generally possible to set lower and upper bounds to the states of (1). These bounds will define a certain state-space region of interest *X* ⊂ ^{n} around the equilibrium point of interest. Hence, this set *X* contains all the practical values of the state variables. Using these pre-specified bounds for the states, it is possible then, by an analysis of the mathematical expressions of (9), to define the upper and lower bounds to these functions *h _{ij}*(

*t*),

*i,j*= 1,...,

*n*. For some cases (when, for example, the functions

*h*(

_{ij}*t*),

*i,j*= 1,...,

*n*are written as combinations of sine and cosine functions), however, it may be possible to ensure the assumption (11) for the entire state-space region

^{n}. In these cases, global stabilization of the nonlinear system (1) can be evaluated via the NLDI model, but this possibility is not investigated in details in this paper.

Comparing the LPV system (13) to the NLDI in the form of (5), it becomes clear now that the proposal of this paper is to model the term *Q _{ij}* θ

*(*

_{ij}*t*) by the term Δ

*A*(

*t*) of (5), which amounts to choosing adequately the matrices

*F*and

*G*. For that, notice that the LPV system (13) can be particularly represented in a polytopic form when (

*A*

_{0}+

*Q*θ

_{ij}_{ij}(

*t*)) involves a polytopic domain of vertices. This particular form of (13) can be interpreted as a polytopic LDI, where all trajectories of the LPV system (13) will also be trajectories of this PLDI. This is discussed in the next section.

**3.2 Specifying a PLDI containing all the trajectories of the LPV system**

In this section, we are interested to model the LPV system (13) as a PLDI model defined in the form of the following linear time-varying system (Pyatnitskiy and Rapoport, 1996):

where the elements of the vector function α(*t*) = [ α_{1}(*t*) ... α_{v}(*t*)]' satisfy the conditions

and *S _{i}* ∈

^{n×n}

*, i*= 1,...,

*v*, is the

*i*vertex of the convex set Ω

^{th}*:=*

_{PLDI}*Co*(

*S*

_{1},

*S*

_{2},...,

*S*). So, we have α

_{v}*(*

_{i}*t*)

*S*∈ Ω

_{i}_{PLDI}for every α(

*t*) whose elements satisfy (15).

Considering the LPV system (13), it is possible to define a set Ω* _{PLDI}* such that (

*A*

_{0}+

*Q*θ

_{ij}*) ∈ Ω*

_{ij}*, for all θ*

_{PLDI}*∈ [*

_{ij}

*h**,*

_{ij}*],*

_{ij}*i,j*= 1,...,

*n*, from which it follows immediately that every trajectory of (13) is also a trajectory of (14). For that, the set of the vertices of Ω

_{PLDI}is defined by

The next step of the procedure consists in searching for a set Ω_{NLDI} that overbounds Ω_{PLDI}. This is discussed in the next section.

**3.3 Overbounding the PLDI by the proposed NLDI model**

Our goal now is to find matrices *F* and *G* for (5) such that Ω_{NLDI} ⊇ Ω_{PLDI}, with the size of the set Ω_{NLDI} as small as possible. This will give an efficient outer approximation (or overbounding) of the PLDI (14) by the NLDI (5), from which it follows immediately that every trajectory of (14) is also a trajectory of (5) (Boyd et al., 1994). Reference (Boyd et al., 1994) formulates this problem in the form of LMIs. The result proposed in (Boyd et al., 1994) considers only the case with matrix *F* square and nonsingular. Here, we relax this condition by a weaker one in which *F* must only be full column rank.

**Proposition 2** Let Ω_{NLDI} and Ω_{PLDI} be the sets defined for systems (5) and (14), respectively. We have Ω_{NLDI} ⊇ Ω* _{PLDI}* if there are a full column rank matrix

*F*and a matrix

*G*with proper dimension, such that

where *A*_{0} is a known matrix in Ω_{NLDI} and *S _{i}, i* = 1,...,

*v*, are the vertices of Ω

_{PLDI}.

**Proof** We have Ω_{NLDI} ⊇ Ω* _{PLDI}* if, for every

*x*(

*t*) and α

*(*

_{i}*t*),

*i*= 1,...,

*v*, such that 0

__<__α

_{i}(

*t*)

__<__1, α

*(*

_{i}*t*) = 1, there exists

*E*(

*t*) satisfying

where *S _{i}* ∈

*V*is the

_{PLDI}*i*vertex of Ω

^{th}_{PLDI}. Notice that (19) is equivalent to the norm condition ||

*E*||

__<__1. We complete the proof by showing that (17) provides an equivalent expression of (18)-(19).

For that, it follows from (17) that

for all α_{i}(*t*), *i* = 1,...,*v*, satisfying 0 __<__ α* _{i}*(

*t*)

__<__1 and α

*(*

_{i}*t*) = 1. Therefore, using Schur complements, the matrix inequality (20) is equivalent to

for all *x*(*t*) and α* _{i}*(

*t*),

*i*= 1,...,

*v*, satisfying 0

__<__α

*(*

_{i}*t*)

__<__1 and α

*(*

_{i}*t*) = 1. It can be concluded from this development that, if (17) is feasible for all

*i*= 1,...,

*v*, then the inequality (21) holds for all

*x*(

*t*) and for all admissibles values of α

*(*

_{i}*t*),

*i*= 1,...,

*v*.

Hence, to complete the proof we only need to show that (21) is equivalent to (18)-(19). This equivalence is easily verified by assuming that *F* is a full column rank matrix. This assumption allows us to rewrite the equation (18) as *E*(*t*)*Gx*(*t*) = *F*^{+}( α* _{i}*(

*t*)

*S*

_{i}- A_{0})

*x*(

*t*). Now, the substitution of this expression into the inequality (19) leads exactly to (21), which completes the desired equivalence and the proof.

To solve the set of matrix inequalities (17) in the form of LMIs it is necessary to introduce the new variables

where *V* ∈ ^{n×n} must be a symmetric and positive semidefinite matrix and *W* ∈ . The condition of positive semidefinitess of matrix *V* can be guaranteed by the additional constraint *V* 0. The condition on matrix *W* can be dealt using LMI solvers specialized in rank constraints (Orsi et al., 2006), but this alternative is also not investigated in this paper. In this paper, the condition *rank*(*W*) = *n _{p}* is reinforced by imposing a particular, desired structure for

*W*and checked after the solution of the matrix inequalities (17) is obtained.

With the change of variables suggested by (22), our problem is to find matrices *V* and *W* such that

where matrices *A*_{0} and *S _{i}, i* = 1,...,

*v*,, were already previously introduced.

Once we have calculated matrices *V* and *W* by solving the LMIs (23), matrix *F* can be easily recovered as *F = W*^{+}. On the other hand, applying the Cholesky-like covariance decomposition to the matrix *V* (and this is possible, once that *V* is a positive semidefinite matrix), we find a matrix *G* with *n _{q} = rank*(

*V*). In order to obtain Ω

_{NLDI}⊇ Ω

*, with the set Ω*

_{PLDI}_{NLDI}as small as possible, reference (Boyd et al., 1994) proposes to solve the LMIs (23) as an optimization problem by minimizing the trace of matrix

*V*. This alternative was adopted in this paper.

The weaker condition imposed on matrix *F* (in which it must only be a full column rank matrix) allow us to adopt a particular structure to the matrix *F* in such a way that the term Δ*A*(*t*) in (5) can better represent the nonlinear characteristics of (1). For example, if the *i ^{th}* equation of (1) is linear, it follows as a suggestion to set the elements of the

*i*row of

^{th}*F*to be equal to zero (which means that the

*i*row of Δ

^{th}*A*(

*t*) will always be equal to zero, independently of the values of

*E*(

*t*)). To better clarify this point and the whole modeling procedure, consider the numerical example of the section 3.5. Before that, however, in section 3.4 the estimation of stability regions for the nonlinear system (1) using NLDIs is discussed.

**3.4 Estimation of regions of attractions for nonlinear systems using NLDIs**

In the sequence, we present a basic result from the Lyapunov theory that provides an estimate *D* ⊂ *X* of the region of attraction of system (1) (Kiyama and Iwasaki, 2000; Coutinho and da Silva Jr., 2010).

**Lemma 3** Consider the nonlinear system (1). Suppose there exist positive scalars є_{1}, є_{2} and є_{3} and a continuously differentiable function *V : X* that satisfies the following conditions:

Then, *V*(*x*) is a Lyapunov Function in *X*. Moreover, for all *x*(0) ∈ *D*, the trajectory *x*(*t*) belongs to *D* and approaches the origin as *t* → ∞.

Next Lemma provides sufficient conditions to ensure that the region of attraction *D* as defined by *D *:= {*x : x'Px* __<__ 1, *P = P' * 0} is bounded by the state-space region *X* (Rohr et al., 2009).

**Lemma 4 **Consider the state-space region *X* defined as (2). The condition *x* ∈ *X* can be written as

Let the domain *D* be defined as *D* := {*x : x'Px* __<__ 1, *P = P'* 0}. Thus, if *x* ∈ *D*, then

hus, the condition x ∈ D ⊂ X is guaranteed if the following inequality holds

Now, from Lemmas 3 and 4 and the ideas presented in (Rohr et al., 2009), it is possible to derive some sufficient conditions to ensure that *D* with *V*(*x*) = *x'Px, P = P' * 0, is an estimate of the region of attraction of the nonlinear system (1) within the set *X*. It is important to emphasize that the calculation of this region *D* is done for the representation of (1) in the form of the NLDI (5). This is allowed to do, once that the NLDI model obtained from the proposed modeling procedure represents the behavior of the nonlinear system (1) within the subset *X*.

**Lemma 5 ** Consider the nonlinear system (1) and its representation in the form of the NLDI (5). Suppose there exist a matrix *P = P' * 0 and a scalar λ > 0 satisfying the following LMIs:

Then, *V*(*x*) = *x'Px* is a Lyapunov Function in *X*. Moreover, for all *x*(0) ∈ *D* the trajectory *x*(*t*) approaches the origin when *t* → ∞, where *D* := {*x : V*(*x*) __<__ 1} ⊂ *X*.

**Proof** Let є_{1} and є_{2} be, respectively, the smallest and largest eigenvalues of *P*. Then, the following inequalities hold for all *x* ∈ *X*

which leads to the condition (24) of Lemma 3. Now, from (Boyd et al., 1994), it follows that if the LMI (30) is satisfied for a matrix *P = P' * 0 and a scalar λ > 0, then we have (*x*) < 0 for all *x* ∈ *X*, where *V*(*x*) = *x'Px*. As *x* is bounded, there exists a sufficiently small positive scalar є_{3} such that

which leads to the condition (25) of Lemma 3. Now, we only have to prove that estimate *D* := {*x : V*(*x*) __<__ 1} is bounded by *X*. To do so, pre- and post-multiply LMI (31) by [1 - *x*]' and its transpose. It provides

Notice that (34) is the matrix form of (29). This completes the proof.

Lemma 5 provides an estimate *D* ⊂ *X* of the region of attraction of the nonlinear system (1), using the representation of this system in the form of the NLDI (5). In order to find the largest estimate *D* inside *X*, we may solve the LMIs in Lemma 5 as an optimization problem by minimizing the trace of matrix *P*, as suggested in (Rohr et al., 2009).

The application of the modeling procedure proposed in this section and the estimate of the system region of attraction are illustrated in the numerical example given in the next section.

**3.5 Numerical example 1**

Consider the following nonlinear system, where the state vector *x*(*t*) is given by *x*(*t*) = [*u*(*t*) *v*(*t*) *z*(*t*)]' (Topcu and Packard, 2009):

Notice that *f*_{1}(*x*(*t*)) and *f*_{2}(*x*(*t*)) are nonlinear functions of the state variables *u*(*t*), *v*(*t*) and *z*(*t*) and *f*_{3}(*x*(*t*)) is a linear function of the state *u*(*t*). Let us consider the equilibrium point at the origin *x _{e}* = [0 0 0]'. The Jacobian matrix

*A*

_{0}was calculated via Taylor series expansion of system (35)-(37) around the equilibrium point

*x*. The result of this process is given by

_{e}The first step of the modeling procedure is to rewrite the nonlinear system (35)-(37) in the form of an LPV system. For that, Proposition 1 of section 3.1 guarantees the existence of two vectors ∈ *Co*(*x*(*t*),0) and ∈ *Co*(*x*(*t*),0), ≠ *x*(*t*), ≠ 0, *i* = 1,2, such that

for all *t* > 0. Notice that, once the function *f*_{3}(*x*(*t*)) is linear, it is not altered by the application of the mean value theorem. As discussed in section 3.1, we can write this last equation as being

where,

Calculating the functions *h _{ij}*(

*t*) (for

*i*= 1,2 and

*j*= 1,2,3) from (38)-(43), we have

The bounds *h** _{ij}* and

_{ij}as defined by (12) were calculated by analysing the mathematical expressions of the functions (44)-(49) considering the following state-space region

This state-space region was specified by assuming it as being an operation region of the system or the region containing the practical values of the states. Once we have delimited the state-space region of interest, the points (*t*), *i* = 1,2, are also bounded, once that (*t*) ∈ *Co*(*x*(*t*),0), for all *t* > 0. So, by analysing the mathematical expressions of the functions (44)-(49) it was possible to define the following upper and lower bounds: *h*_{11} = -0.126, _{11} = 0.126, *h*_{12} = 0.150, _{12} = 3.106, *h*_{13} = -0.015, _{13} = 0.015, *h*_{22} = -0.329, _{22} = 0.329, *h*_{23} = -0.329, _{23} = 0.329.

To better understand the calculation of these bounds, let us take *h*_{11}(*t*) as an example. It is function of the point (*t*), whose value must be in the line segment between - and , for all *t *> 0. Hence, the maximum value that *h*_{11}(*t*) can achieve (*i.e*, the upper bound _{11}) is 0.126 and the smallest value (*i.e*, the lower bound *h*_{11}) is -0.126.

With the specified bounds to the functions (44)-(49), we have a description of the nonlinear system (35)-(37) in the form of the following LPV system

where θ* _{ij }*: ℜ

_{+}→ [

*h**,*

_{ij}*],*

_{ij}*i*= 1,2,

*j*= 1,2,3 (

*j*≠ 1 when

*i*= 2).

The second step of the modeling procedure consists of specifying a PLDI in the form of (14) containing all the trajectories of the LPV system (51). For that, we only need to construct a set Ω_{PLDI} as defined by (16). Once that the LPV system (51) has 5 non-zero functions θ_{ij}(*t*), so the corresponding set Ω_{PLDI} has 32 vertices.

To complete the modeling procedure, we have to overbound the PLDI by the proposed NLDI model. This is done by solving the LMI problem introduced in section 3.3. The two matrices of the model that we must calculate are *G* and *F*. The standard method proposed in (Boyd et al., 1994) requires that *F* be a square and nonsingular matrix so, in this case, we should have to calculate a full rank 3×3 matrix. As a result, the third line of the term Δ*A*(*t*) of the NLDI model can have non-zero elements. But, notice that *f*_{3}(*x*(*t*)) of our example is linear, so the third line of Δ*A*(*t*) is expected to be equal to zero (which means that all the elements of this line are equal to zero) to avoid conservatism.

On the other hand, by applying the relaxing condition on matrix *F*, as proposed in Proposition 2 of section 3.3, we can set a matrix *F* in the form *F* = [*f*_{11} *f*_{21} 0]'. In this case, the row with the element equal to zero (row 3) forces the respective row of the matrix Δ*A*(*t*) to also have all its elements equal to zero, independently of any value of *E*(*t*) and *G*. Hence, the obtained NLDI model will mantain the characteristic of the study system with respect to the fact that *f*_{3}(*x*(*t*)) is linear. Another possible advantage of Proposition 2 (in comparison to the standard method proposed in (Boyd et al., 1994)) is the fact that the structure of the matrix *F* can be choosen in order to reduce the number of elements to be determined by the LMI optimization problem. To better clarify this point, notice that matrix *F* has 9 elements to be calculated by applying the standard method, while the proposed method allows us to specify a matrix *F* with only 2 elements to be calculated.

For the matrix *G*, we set it as a 3×3 matrix, so the matrix variable *V* was choosen to be a 3×3 symmetric matrix. In addition, we imposed the following structure to the matrix variable *W* in order to obtain a matrix *F* in the form of *F* = [*f*_{11} *f*_{21} 0]':

From the solution of the optimization problem suggested in section 3.3 using the SeDuMi solver (Sturm, 1999) in conjunction with YALMIP (Lofberg, 2004) it was obtained the following matrices *F* and *G*

This complete the modeling procedure. Now, let us estimate the region of attraction of the nonlinear system (35)-(37) with respect to the origin by using its representation in the form of the NLDI model previously calculated. This was done using Lemma 5. In order to obtain a solution for the LMI (31), we have to describe the set *X* in the form given by (2). This was done by following the ideas presented in (Rohr et al., 2009). At first, notice that (50) is equivalent to a polytope *X* ⊂ ^{3} whose vertices are defined by Θ := {*c*_{1},*c*_{2},...,*c*_{8}}, where

A vertex representation of *X* is defined as the convex hull of *c*_{1},*c*_{2},...,*c*_{8}, *i.e., X = Co*(*c*_{1},*c*_{2},...,*c*_{8}). Equivalently, we can define this vertex form of *X* as (2):

In order to calculate the vectors *a _{k}* ∈

^{3}, notice that each inequality defines an hyperplane {

*x : a*

_{k}x__<__1} for which two vertices belongs to its edge, that is

*c*_{1},*c*_{2} ∈ {*x : a*_{1}*x* = 1}, *c*_{2},*c*_{3} ∈ {*x : a*_{2}*x* = 1},

*c*_{3},*c*_{4} ∈ {*x : a*_{3}*x* = 1}, *c*_{1},*c*_{4} ∈ {*x : a*_{4}*x* = 1},

*c*_{4},*c*_{5} ∈ {*x : a*_{5}*x* = 1}, *c*_{5},*c*_{6} ∈ {*x : a*_{6}*x* = 1},

*c*_{6},*c*_{7} ∈ {*x : a*_{7}*x* = 1}, *c*_{7},*c*_{8} ∈ {*x : a*_{8}*x* = 1},

*c*_{5},*c*_{8} ∈ {*x : a*_{9}*x* = 1}, *c*_{1},*c*_{8} ∈ {*x : a*_{10}*x* = 1},

*c*_{2},*c*_{7} ∈ {*x : a*_{11}*x* = 1}, *c*_{3},*c*_{6} ∈ {*x : a*_{12}*x* = 1}.

Figure 1 shows the set *X* with its vertices and illustration of some of its edges ({*x : a*_{1}*x* = 1}, {*x : a*_{7}*x* = 1} and {*x : a*_{11}*x* = 1}).

For a given set Θ, the row vectors *a _{k}* can be determined by solving the following set of linear systems:

yielding the following row vectors

*a*_{1} =[ 0.29 -0.45 0 ], *a*_{2} = [ 0 -0.32 0.32 ],

*a*_{3} = [ -0.29 -0.45 0 ], *a*_{4} = [ 0 -0.32 -0.32 ],

*a*_{5} = [ -0.29 0 -0.45 ], *a*_{6} = [ -0.29 0.45 0 ],

*a*_{7} = [ 0 0.32 0.32 ], *a*_{8} = [ 0.29 0.45 0 ],

*a*_{9} = [ 0 0.32 -0.32 ], *a*_{10} = [ 0.29 0 -0.45 ],

*a*_{11} = [ 0.29 0 0.45 ], *a*_{12} = [ -0.29 0 0.45 ].

The ellipsoid *D* := {*x* ∈ ^{3} | *x'P x* __<__ 1} was used for the estimation of the region of attraction of the nonlinear system. To calculate the largest estimate in region *X*, the LMIs of Lemma 5 were solved as an optimization problem by minimizing the trace of matrix *P*. The calculated matrix *P* is

Figure 2 shows the estimated region of attraction *D* and the state-space region *X* defined as (50). It is important to emphasize that *D* is an estimate of the region of attraction of the nonlinear system (35)-(37) with respect to the origin. It was calculated via Lemma 5 by using the representation of the nonlinear system in the form of the NLDI previously calculated.

Next section presents the theory towards the design technique of the dynamic output feedback controller.

**4 DYNAMIC OUTPUT FEEDBACK CONTROLLER DESIGN**

Let the controlled system formed by the interconnection between (5) and (6)-(7) (with *y*(*t*) = *Cx*(*t*)) be described as

where Δ*Ã*(*t*) = *E*(*t*), with* E*(*t*) satisfying ||*E*(*t*)|| __<__ 1, for all *t* > 0. Besides, (*t*) = [*x*(*t*) *x _{c}*(

*t*)]', and

In many applications, it is desirable that the system trajectories approach the origin as fast as possible (Folcher and Ghaoui, 1994; Silva and Junior, 2006). This practical requirement can be fulfilled by guaranteeing that lim_{t }_{→ ∞}*e* ^{α}^{t}||*x*(*t*)|| = 0 for all trajectories of (55), where α is defined as being the decay rate of (55). This can be done by satisfying the condition ((*t*)) __<__ -2 α*V*((*t*)) (Boyd et al., 1994), where *V*((*t*)) = (*t*)'(*t*) is the Lyapunov function adopted to the NLDI model. This condition is equivalent to the existence of a positive definite symmetric matrix ∈ and matrices *A _{c}, B_{c}* and

*C*such that

_{c}In addition, it may be desirable that all modes of response (eigenvalues) of the closed-loop matrix *Ã*_{0} have damping ratios larger than a minimum, pre-defined value (see the applications in (Ramos et al., 2004; de Oliveira et al., 2009)). To comply with this, we impose an additional restriction on the problem formulation using the regional pole placement (RPP) technique (Chiali et al., 1999). This technique consists in the definition of a region for pole placement in the complex plane where the design objective is fulfilled. This region is defined by all the complex numbers that have a damping ratio ξ higher (or equal) than ξ* _{min}*, and it can be viewed in Figure 3, where ξ

*is the desired minimum damping ratio for the eigenvalues of matrix*

_{min}*Ã*

_{0}and δ = cos

^{-1}ξ

_{min}. Note that the damping ratio is a local property, which only makes sense when the trajectory is close enough to the equilibrium point of interest.

All the eigenvalues of matrix *Ã*_{0} are located within the region specified in Figure 3 if there exists a positive definite symmetric matrix and matrices *A _{c}, B_{c}* and

*C*such that (Chiali et al., 1999)

_{c}**4.1 Robustness of the controller with respect to the variations in the operating conditions of the system**

To deal with the robustness of the controller with respect to the different operating points of the system, we obtain a description of the system (1) in a certain state-space region around each equilibrium point of interest. The idea is to design a fixed parameter controller in the form (6)-(7) that exhibits an effective performance in all of these regions. This controller, however, must not change the equilibrium points of the open-loop sytem. In fact, once that a certain initial condition is within the region of attraction of a particular equilibrium (say, *x _{e}*

_{1}), then the control objective is to force the system trajectories to approach the point

*x*

_{e}_{1}as fast as possible, in accordance to the practical requirements discussed in the previous section.

The set of resulting controlled systems are described in state space form by

where Δ*Ã _{i}*(

*t*) =

*(*

_{i}E_{i}*t*)

*, with*

_{i}*E*(

_{i}*t*) satisfying ||

*E*(

_{i}*t*)||

__<__1, for all

*t*> 0. We also have

being *x _{ei}* the

*i*equilibrium point of interest of (1).

^{th}**4.2 The complete control problem**

Grouping the constraints (56)-(57) with the robusteness requirement shown in previous subsection, it is possible to formulate the complete control problem as a search for positive definite symmetric matrices _{i}*, i* = 1,...,*n _{p}*, and matrices

*A*and

_{c}, B_{c}*C*of proper dimensions such that

_{c}Notice that (59)-(60) are bilinear matrix inequalities (BMIs), since there are cross-products among the controller variables (*i.e., A _{c}, B_{c}* and

*C*) and matrices

_{c}_{i}

*, i*= 1,...,

*n*. In this paper, we apply a two-step separation procedure that allows to transform the BMIs (59)-(60) into a set of LMIs (de Oliveira et al., 2000). Basically, this separation procedure consists on the parameterization of some matrix variables and the definition of some new variables.

_{p}

**5 SOLVING THE BMI PROBLEM: THE TWO-STEP SEPARATION PROCEDURE**

The concepts and procedures described in this section are derived from the ideas presented in (de Oliveira et al., 2000). Consider the set of BMIs (59)-(60). Let us assume that: (*i*) = _{i}*, i* = 1,...,*n _{p}*; (

*ii*) the dimension of the controller is equal to the dimension of the plant to be controlled,

*i.e., n*, and; (

_{c}= n*iii*) the controller output matrix

*C*is previously known. Let us partionate the matrices and its inverse

_{c}^{-1}and define a matrix as follows

where *X,U,X _{c},Y,Y_{c}* ∈ ℜ

^{n×n}and the dimensions of and its submatrices are implicitly determined by . Moreover, the following changes of variables are carried out

where the dimensions of matrices *V, P* and *S* are implicity determined by the transformations.

Now, the set of BMIs (59)-(60) can be transformed into a set of LMIs. For that, multiply , (59) and (60) on the right and the left by *T, diag*[*T,I*] and *diag*[*T,T*], respectively; introduce the new variables *V, P* and *S*, and simplify the expressions by algebraic manipulations, remembering that ^{-1} = *I*. The resulting set of LMIs are given by

where, *D*_{11} = *P* + *P*[ A]_{0i }+ *G _{i }*+2 α

*P*,

*D*

_{12}=

*PA*

_{0i }+

*X*+

*V*+

^{T }*S*+

*G*+ 2 α

_{i }*P, D*

_{13}=

*PF*

_{i}, D_{22}=

*X*+

*XA*

_{0i }+

*VC*+

_{i }*V*+

^{T }*G*+ 2 α

_{i }*X, D*

_{23}=

*XF*

_{i}, D_{33}= -

*I, N*

_{11}= (

*P*+

*P*[ A]

_{0i})sin( δ),

*N*

_{12}= (

*X*+

*P*[ A]

_{0i }+

*V*+

^{T }*S*sin( δ),

*N*

_{13}= (

*P*-

*P*[ A]

_{0i})cos( δ),

*N*

_{14}= (

*X*-

*PA*

_{0i }+

*V*)cos( δ),

^{T }+ S*N*

_{22}= (

*X + XA*

_{0i }+

*V*)sin( δ),

^{T }+ VC_{i}*N*

_{24}= (

*X - XA*

_{0i }+

*V*)cos( δ),

^{T }- VC_{i}*N*

_{23}=

*N'*

_{14},

*N*

_{33}=

*N*

_{11},

*N*

_{34}=

*N*

_{12},

*N*

_{44}=

*N*

_{22}, [ A]

_{0i}=

*A*

_{0i }

*+ B*= 1,...,

_{i}C_{c}, i*n*.

_{p} Solving this set of LMIs in the matrices variables *V, P, S* and *X*, the matrices *A _{c}* and

*B*of the controller (remembering that matrix

_{c}*C*must be pre-specified) can be calculated by

_{c}where *U = P - X*. The success in solving this LMI problem depends, obviously, in a proper choice of matrix *C _{c}*. In this paper, we calculate this matrix by setting up a state feedback gain

*K*, in which the control law

*u*(

*t*) =

*Kx*(

*t*) stabilizes the system (5) and fulfills the constraints discussed in section 4. This matrix

*K*can be found by solving the following LMIs in the matrix variables L and Y (Ramos et al., 2004):

where, *Q*_{11} = *Y* + *A*_{0i}*Y + L ^{T}* +

*B*+ 2 α

_{i}L*Y, Q*

_{12}=

*Y*,

*Q*

_{22}= -

*I, R*

_{11}= (

*Y*+

*A*

_{0i}

*Y + L*+

^{T}*B*)sin( δ),

_{i}L*R*

_{12}= (

*Y*-

*A*

_{0i}

*Y + L*-

^{T}*B*)cos( δ),

_{i}L*R*

_{22}=

*R*

_{11},

*i*= 1,...,

*n*.

_{p} Once we have solved this set of LMIs (65), it is setled *C _{c}: = K*, where

*K = LY*

^{-1}.

It is important to emphasize that given the non-convex nature of the undepinning BMI, we cannot guarantee that the LMI problem that is created by setting the matrix Cc obtained from the solution of the state feedback problem will always have a solution. Recent papers, however have shown that this heuristics for calculation of matrix *C _{c}* provides satisfactory results for the overall design procedure, as seen, for example, in (Ramos et al., 2004; Ramos et al., 2005; Kuiava et al., 2009).

**6 TESTS AND RESULTS**

Consider the following nonlinear system:

We consider three equilibrium points of interest: *x _{e}*

_{1}= [0 0 0]

^{T}

*, x*

_{e}_{2}= [0.3074 2.1277 -4.2718]

^{T}and

*x*

_{e}_{3}= [-0.034 -0.5368 1.4405]

^{T}. The Jacobian matrices

*A*

_{01},

*A*

_{02}and

*A*

_{03}were calculated via Taylor series expansion of system (66)-(68) around the equilibrium points

*x*

_{e}_{1},

*x*

_{e}_{2}and

*x*

_{e}_{3}, respectively. They are:

An eigenvalue analysis shows that these three equilibrium points are all locally asymptotically stable. Our goal is then to design a dynamic output feedback controller that improves the decay rate of the system trajectories, as well as, the damping ratio of these trajectories as they are approaching the equilibrium points. For that, we first modelled the nonlinear system (66)-(68) via three NLDIs in the form (5), each one describing a certain neighborhood of the points *x _{e}*

_{1},

*x*

_{e}_{2}and

*x*

_{e}_{3}

^{1}.

Let us discuss the construction of an NLDI model describing a certain state-space region of the studied nonlinear system containing the equilibrium point at the origin. Then, the same approach was applied to determine the other NLDIs associated to the points *x _{e}*

_{2}and

*x*

_{e}_{3}. Applying the mean value theorem to the nonlinear system (66)-(68) with respect to

*x*

_{e}_{1}, we have a reformulation of it in the form of the following LPV system

where θ* _{ij}* : ℜ

_{+}→ [

*h**,*

_{ij}*],*

_{ij}*i,j*= 1,...,3.

The bounds *h** _{ij}* and

*were specified by analysing the mathematical expressions of the functions*

_{ij}*h*(

_{ij}*t*) (which are calculated by (9)) considering the following state-space region:

The region *X*_{1} was defined by considering it as being the operating range with respect to the equilibrium point at the origin. Thus, we specified the following bounds for the non-zero functions *h _{ij}*(

*t*):

__h___{11}= -0.020,

_{11}= 0.024,

*h*_{12}= -0.236,

_{12}= 0.316,

*h*_{13}= -0.0025,

_{13}= 0.003,

*h*_{22}= -0.09,

_{22}= 0.086,

*h*_{23}= -0.05,

_{23}= 0.055,

*h*_{32}= -0.48,

_{32}= 0.56,

*h*_{33}= -0.025 and

_{33}= 0.025.

The LPV system (71) has 7 non-zero functions *h _{ij}*(

*t*), which means that the corresponding set Ω

_{PLDI}has 128 vertices (as defined by (16)). To complete the procedure, we have to calculate matrices

*F*

_{1}and

*G*

_{1}. Different from the numerical example of section 3.5, all the equations of the system (66)-(68) are nonlinear. Hence, none of the elements of

*F*

_{1}was forced to be equal to zero. The dimension of this matrix

*F*

_{1}was settled to 3×1. This dimension was chosen instead of 3×2 or 3×3 in order to reduce the number of elements to be calculated by the LMI optimization problem. The dimension of matrix

*G*

_{1}was chosen to be 3×3.

By especifying these dimensions to matrices *F*_{1} and *G*_{1}, matrices *W* and *V* were defined as being a 1×3 full row rank matrix and a 3×3 symmetric matrix, respectively. From the solution of the optimization problem suggested in section 3.3 using the SeDuMi solver (Sturm, 1999) in conjunction with YALMIP (lofberg, 2004) it was possible to obtain the following matrices

The NLDIs describing the system dynamics around the equilibrium points *x _{e}*

_{2}and

*x*

_{e}_{3}were calculated considering, respectively, the following regions:

The resulting set of NLDIs are given in the form

where *i* = 1,2,3 and *x*^{(i)}(*t*) = *x*(*t*) - *x _{ei}*. In addition,

*A*

_{01},

*A*

_{02}and

*A*

_{03}are given by (70), respectively;

*F*

_{1}and

*G*

_{1}by (73) and the other matrices are

Finally, matrices *B* and *C* are directly determined by looking the nonlinear equations (66)-(69).

Now, let us estimate the region of attraction *D*_{1} of the equilibrium point *x _{e}*

_{1}of the nonlinear system (66)-(69) using the corresponding NLDI model previously calculated. This was done by using Lemma 5. As already shown in section 3.5 for the numerical example 1, we have first to describe the set

*X*

_{1}in the form of (2). This description can be obtained using the same procedure adopted for that example of section 3.5. As a result, we have a set

*X*

_{1}described in the form given by (2), where

*a*_{1} = [ 2.00 -2.00 0 ], *a*_{2} = [ 0 -1.64 1.97 ],

* a*_{3} = [ -2.00 -2.00 0 ] , *a*_{4} = [ 0 -1.64 -1.97 ] ,

* a*_{5} = [ -1.64 0 -1.97 ] , *a*_{6} = [ -1.64 1.97 0 ] ,

* a*_{7} = [ 0 1.67 1.67 ] , *a*_{8} = [ 1.64 1.97 0 ] ,

* a*_{9} = [ 0 1.67 -1.67 ] , *a*_{10} = [ 1.64 0 -1.97 ] ,

* a*_{11} = [ 1.64 0 1.97 ] , *a*_{12} = [ -1.64 0 1.97 ] .

Figure 4 shows the ellipsoid *D*_{1} and the state-space region *X*_{1} defined as (72).

Now, our goal is to design a dynamic output feedback controller to the nonlinear system (66)-(69) using a description of it in the form of (76)-(77). A decay rate of 0.01 and a minimum damping ratio equal to 15% were imposed as design objectives to the controlled system. The SeDuMi solver, used in conjunction with YALMIP, was used to solve the set of LMIs related to this control problem. The calculated controller is given by

The performance of the controlled system was verified via nonlinear simulations. Figs. 5, 6 and 7 show the response of variables *x*_{1}(*t*), *x*_{2}(*t*) and *x*_{3}(*t*), respectively, with respect to different initial conditions. It is interesting to observe that the trajectory of the controlled system is less oscillatory and approaches the equilibrium faster than the trajectory of the open-loop system. This result shows the effectiveness of the designed controller.

The region of attraction of the equilibrium point _{e1} of the closed-loop nonlinear system, where _{e1} = [*x' _{e}*

_{1}0 0 0]', was also estimated using the corresponding NLDI model in the form of (58). For that, we defined the ellipsoid

_{1 }:= { ∈

^{6}| '

_{1}

__<__1}.

Figure 8 shows the ellipsoid _{1}, considering the state of the controller equal to zero. This plot corresponds to the cut of the actual estimate of ellipsoid _{1} in the hyperplane defined by the system states. Hence, from the result shown in Figure 8, the effectiveness of the designed controller for the nonlinear system (66)-(69) is guaranteeded within the region _{1}.

**7 CONCLUSION**

In the first part of this paper, a method to calculate the parameters of an NLDI model was presented. The objective is to obtain a suitable linear representation of a nonlinear system for control purposes, in such a way that a linear controller can be designed to guarantee some desired features for the nonlinear closed-loop system. In the second part, a robust control design method, written in terms of LMIs was presented, in order to design a linear dynamic output feedback controller for a nonlinear system using the NLDI model obtained by the approah proposed in the first part.

The numerical examples presented in the previous section have shown the effectiveness of the modeling and control approach proposed in this paper. However, it is important to emphasize that the application of the modeling procedure may be restricted to small and medium size systems, once that the number of vertices of the obtained PLDI may become excessively high as the number of system nonlinearities increases. With respect to the proposed control design, one of its difficulties is the requirement, due to the nature of the control problem formulation, that the order of the designed LDOF controller must be equal to the order of the system to be controlled. In order to achieve the goal of producing low order controllers (in the cases where the dimension of the system is sufficiently high), the control design may be combined with a final step of model order reduction.

The design objectives targeted in this paper were a minimum decay rate for all possible system trajectories within the operating region of interest and a minimum damping ratio to the oscillation modes of the state matrix of the closed-loop system. These design objectives are very important for the problem of robust damping control in power systems, for example, and an application of the approach proposed in this paper to such a problem is the next predicted step of this research. It is important to point out, however, that other practical requirements can be included in the control problem formulation, such as the minimization of an output energy function or restrictions on the controller bandwidth. These are also possible extensions of this research foreseen to the sequence of this research.

**ACKNOWLEDGE**

The authors would like to thank the reviewers for their helpful comments and suggestions.

The authors would like also to thank Prof. Luis Fernando Costa Alberto for providing the Matlab routines used for the estimates of the regions of attractions in examples 1 and 2.

**REFERENCES**

Basler, M. J. and Schaefer, R. C. (2008). Understanding power system stability, *IEEE Trans. Indus. Appl.* **44**(2): 463–474. [ Links ]

Bernard, F., Dufour, F. and Bertrand, P. (1997). On the jlq problem with uncertainty, *IEEE Trans. Automat. Contr.* **42**(6): 869–872. [ Links ]

Boyd, S., Ghaoui, L. E., Feron, E. and Balakrishnam, V. (1994). *Linear Matrix Inequalities in System and Control Theory*, Society for industrial and applied mathematics. [ Links ]

Chen, B. S., Lee, C. H. and Chang, Y. C. (1996). H_{∞} tracking design of linear systems: adaptative fuzzy approach, *IEEE Trans. Fuzzy Syst.* **4**(5): 32–43. [ Links ]

Chiali, M., Gahinet, P. and Apkarian, P. (1999). Robust pole placement in lmi regions, *IEEE Trans. Automat. Contr.* **44**(12): 2257–2270. [ Links ]

Coutinho, D. F. and da Silva Jr., J. M. G. (2010). Computing estimates of the region of attraction for rational control systems with saturating actuators, *IET Control Theory Applications* **4**(3): 315–325. [ Links ]

de Oliveira, M. C., Geromel, J. C. and Bernussou, J. (2000). Design of dynamic output feedback decentralized controllers via a separation procedure, *Int. J. Contr.* **73**(5): 371–381. [ Links ]

de Oliveira, R. V., Kuiava, R., Ramos, R. A. and Bretas, N. G. (2009). Automatic tuning method for the design of supplementary damping controllers for flexible alternating current transmission system devices, *IET Gen., Trans. & Distrib.* **3**(10): 919–929. [ Links ]

Folcher, J. P. and Ghaoui, E. L. (1994). State-feedback design via linear matrix inequalities: application to a benchmark problem, *Proc. of the IEEE Conf. on Contr. Appl. [ Links ]*

Hossain, J., Pota, H. R., Ugrinovskii, V. and Ramos, R. A. (2009). A novel STATCOM control to augment LVRT of fixed speed wind generators, *Proc. of the IEEE Conf. on Dec. and Contr. [ Links ]*

Hu, T. (2007). Nonlinear control design for linear differential inclusions via convex hull of quadratics, *Automatica* **43**(4): 685–692. [ Links ]

Hu, X. B. and Chen, W. H. (2007). Model predictive control of nonlinear systems: stability region and feasible initial control, *Int. J. Automation and Comput.* **4**(2): 195–202. [ Links ]

Kiyama, T. and Iwasaki, T. (2000). On the use of multi-loop circle criterion for saturating control synthesis, *System and Control Letters* **41**: 105–114. [ Links ]

Kuiava, R., Ramos, R. A. and Bretas, N. G. (2009). Robust control methodology for the design of supplementary damping controllers for FACTS devices, *Revista Controle & Automao* **20**(2): 192–205. [ Links ]

Lofberg, J. (2004). Yalmip: a toolbox for modeling and optimization in matlab, *Proc. of the CACSD Conf.*, Taipei, Taiwan. *http://control.ee.ethz.ch/joloef/yalmip.php [ Links ]

Montagner, V., Oliveira, R. C. L. F. and Peres, P. L. D. (2010). Relaxaes convexas de convergncia garantida para o projeto de controladores para sistemas nebulosos de takagi-sugeno, *Revista Controle & Automao ***22**(1): 82–95. [ Links ]

Mozelli, L. A., Palhares, R. M., de Avellar, G. S. C. and dos Santos, R. F. (2010). Condies lmis alternativas para sistemas takagi-sugeno via funo de lyapunov fuzzy, *Revista Controle & Automao* **21**(1): 96–107. [ Links ]

Orsi, R., Helmke, U. and Moore, J. B. (2006). A newtonlike method for solving rank constrained linear matrix inequalities, *Automatica* **42**(11): 1875–1882. [ Links ]

Polanski, A. (1997). Lyapunov function construction by linear programming, *IEEE Trans. Automat. Contr.* **42**(7): 1013–1016. [ Links ]

Pyatnitskiy, E. S. and Rapoport, L. B. (1996). Criteria of asymptotic stability of differential inclusions and periodic motions of time-varying nonlinear control systems, *IEEE Trans. Circ. Syst. I* **43**(3): 219–229. [ Links ]

Ramos, R. A., Alberto, L. F. C. and Bretas, N. G. (2004). A new methodology for the coordinated design of robust decentralized power system damping controllers, *IEEE Trans. Power Syst.* **19**(1): 444–454. [ Links ]

Ramos, R. A., Martins, A. C. P. and Bretas, N. G. (2005). Improved methodology for the design of power system damping controllers, *IEEE Trans. Power Syst.* **20**(4): 1938–1945. [ Links ]

Rohr, E. R., Pereia, L. F. A. and Coutinho, D. F. (2009). Robustness analysis of nonlinear systems subject to state feedback linearization, *SBA Controle e Automao* **20**(4): 482–489. [ Links ]

Silva, S. and Junior, V. L. (2006). Active flutter suppression in a 2-d airfoil using linear matrix inequalities techniques, *J. Braz. Soc. Mech. Sci. Eng.* **28**(1): 84–93. [ Links ]

Sturm, J. F. (1999). Using SeDuMi 1.02, a Matlab toolbox for optimization over symmetric cones, *Optimization Methods and Software* **11**(1): 625–653. [ Links ]

Tognetti, E. S. and Oliveira, V. A. (2010). Fuzzy pole placement based on piecewise lyapunov functions, *Int. J. Robust and Nonlin. Contr.* **20**(1): 571–578. [ Links ]

Topcu, U. and Packard, A. (2009). Local stability analysis for uncertain nonlinear systems, *IEEE Trans. Automat. Contr.* **54**(5): 1042–1047. [ Links ]

Vidyasagar, M. (1993). *Nonlinear systems analysis*, Englewood Cliffs, N.J: Prentice Hall. [ Links ]

Xie, L. and de Souza, C. E. (1992). Robust control for linear systems with norm-bounded time-varying uncertainty, *IEEE Trans. Automat. Contr.* **37**(8): 1188–1191. [ Links ]

Yfoulis, C. A. and Shorten, R. (2004). A numerical technique for the stability analysis of linear switched systems, *Int. J. Contr.* **77**(11): 1019–1039. [ Links ]

Zemouche, A., Boutayeb, M. and Bara, G. (2005). Observer design for nonlinear systems: an approach based on the differential mean value theorem, *Proc. of the IEEE Conf. on Dec. and Contr. [ Links ]*

Artigo submetido em 24/03/2011 (Id.: 1308)

Revisado em 05/05/2011, 24/08/2011, 06/09/2011

Aceito sob recomendação do Editor Associado Prof. Daniel Coutinho