Open access

Robust H Fuzzy Control Design for Nonlinear Two-Time Scale System with Markovian Jumps based on LMI Approach

Written By

Wudhichai Assawinchaichote, Sing Kiong Nguang and Non-members

Published: 01 October 2009

DOI: 10.5772/8230

From the Edited Volume

Advanced Technologies

Edited by Kankesu Jayanthakumaran

Chapter metrics overview

1,929 Chapter Downloads

View Full Metrics

1. Introduction

Markovian jump systems consists of two components; i.e., the state (differential equation) and the mode (Markov process). The Markovian jump system changes abruptly from one mode to another mode caused by some phenomenon such as environmental disturbances, changing subsystem interconnections and fast variations in the operating point of the system plant, etc. The switching between modes is governed by a Markov process with the discrete and finite state space. Over the past few decades, the Markovian jump systems have been extensively studied by many researchers; see (Kushner, 1967; Dynkin, 1965; Wonham, 1968; Feng et al., 1992; Souza & Fragoso, 1993; Boukas & Liu, 2001 ; Boukas & Yang, 1999; Rami & Ghaoui, 1995; Shi & Boukas, 1997; Benjelloun et al., 1997; Boukas et al., 2001 ; Dragan et al., 1999). This is due to the fact that jumping systems have been a subject of the great practical importance. For the past three decades, singularly perturbed systems have been intensively studied by many researchers; see (Dragan et al., 1999; Pan & Basar, 1993; Pan & Basar, 1994; Fridman, 2001; Shi & Dragan, 1999; Kokotovic et al., 1986).

Multiple time-scale dynamic systems also known as singularly perturbed systems normally occur due to the presence of small “parasitic” parameters, typically small time constants, masses, etc. In state space, such systems are commonly modelled using the mathematical framework of singular perturbations, with a small parameter, say ε, determining the degree of separation between the “slow” and “fast” modes of the system. However, it is necessary to note that it is possible to solve the singularly perturbed systems without separating between slow and fast mode subsystems. But the requirement is that the “parasitic” parameters must be large enough. In the case of having very small “parasitic” parameters which normally occur in the description of various physical phenomena, a popular approach adopted to handle these systems is based on the so-called reduction technique. According to this technique the fast variables are replaced by their steady states obtained with “frozen” slow variables and controls, and the slow dynamics is approximated by the corresponding reduced order system. This time-scale is asymptotic, that is, exact in the limit, as the ratio of the speeds of the slow versus the fast dynamics tends to zero.

In the last few years, the research on singularly perturbed systems in the H sense has been highly recognized in control area due to the great practical importance. H optimal control of singularly perturbed linear systems under either perfect state measurements or imperfect state measurements has been investigated via differential game theoretic approach. Although many researchers have studied the H control design of linear singularly perturbed systems for many years, the H control design of nonlinear singularly perturbed systems remains as an open research area. This is due to, in general, nonlinear singularly perturbed systems can not be decomposed into slow and fast subsystems.

Recently, a great amount of effort has been made on the design of fuzzy H for a class of nonlinear systems which can be represented by a Takagi-Sugeno (TS) fuzzy model; see (Nguang & Shi, 2001; Han & Feng, 1998; Chen et al., 2000; Tanaka et al., 1996). Recent studies (Nguang & Shi, 2001; Han & Feng, 1998; Chen et al., 2000; Tanaka et al., 1996; Wang et al., 1996) show that a fuzzy model can be used to approximate global behaviors of a highly complex nonlinear system. In this fuzzy model, local dynamics in different state space regions are represented by local linear systems. The overall model of the system is obtained by “blending” of these linear models through nonlinear fuzzy membership functions. Unlike conventional modelling which uses a single model to describe the global behavior of a system, fuzzy modelling is essentially a multi-model approach in which simple submodels (linear models) are combined to describe the global behavior of the system. Employing the existing fuzzy results (Nguang & Shi, 2001; Han & Feng, 1998; Chen et al., 2000; Tanaka et al., 1996; Wang et al., 1996) on the singularly perturbed system, one ends up with a family of ill-conditioned linear matrix inequalities resulting from the interaction of slow and fast dynamic modes. In general, ill-conditioned linear matrix inequalities are very difficult to solve.

What we intend to do in this chapter is to design a robust H fuzzy controller for a class of uncertain nonlinear two time-scale dynamic systems with Markovian jumps. First, we approximate this class of uncertain nonlinear two time-scale dynamic systems with Markovian jumps by a Takagi-Sugeno fuzzy model with Markovian jumps. Then based on an LMI approach, we develop a technique for designing a robust H fuzzy controller such that the L 2-gain of the mapping from the exogenous input noise to the regulated output is less than a prescribed value. To alleviate the ill-conditioned linear matrix inequalities resulting from the interaction of slow and fast dynamic modes, these ill-conditioned LMIs are decomposed into ε-independent LMIs and ε-dependent LMIs. The ε-independent LMIs are not ill-conditioned and the ε-dependent LMIs tend to zero when ε approaches to zero. If ε is sufficiently small, the original ill-conditioned LMIs are solvable if and only if the ε-independent LMIs are solvable. The proposed approach does not involve the separation of states into slow and fast ones, and it can be applied not only to standard, but also to nonstandard singularly perturbed systems.

This chapter is organized as follows. In Section 2, system descriptions and definition are presented. In Section 3, based on an LMI approach, we develop a technique for designing a robust H fuzzy controller such that the L 2-gain of the mapping from the exogenous input noise to the regulated output is less than a prescribed value for the system described in Section 2. The validity of this approach is demonstrated by an example from a literature in Section 4. Finally, conclusions are given in Section 5.

Advertisement

2. System Descriptions and Definitions

The class of uncertain nonlinear two time-scale dynamic systems with Markovian jumps under consideration is described by the following TS fuzzy model with Markovian jumps:

E ε x ˙ ( t ) = i = 1 r μ i ( υ ( t ) ) [ [ A i ( η ( t ) ) + Δ A i ( η ( t ) ) ] x ( t ) + [ B 1 i ( η ( t ) ) + Δ B 1 i ( η ( t ) ) ] ω ( t ) + [ B 2 i ( η ( t ) ) + Δ B 2 i ( η ( t ) ) ] u ( t ) ] x ( 0 ) = 0 z ( t ) = i = 1 r μ i ( υ ( t ) ) [ [ C 1 i ( η ( t ) ) + Δ C 1 i ( η ( t ) ) ] x ( t ) + [ D 12 i ( η ( t ) ) + Δ D 12 i ( η ( t ) ) ] u ( t ) ] E1

where E ε = [ I 0 0 ε I ] ε 0 is the singular perturbation parameter, υ ( t ) = [ υ 1 ( t ) υ ϑ ( t ) ] is the premise variable that may depend on states in many cases, μ i ( υ ( t ) ) denote the normalized time varying fuzzy weighting functions for each rule, ϑ is the number of fuzzy sets, x ( t ) n is the state vector u ( t ) m is the input, w ( t ) p is the disturbance which belongs to L 2[0,∞), z ( t ) s is the controlled output, the matrix functions A i ( η ( t ) ) B 1 i ( η ( t ) ) B 2 i ( η ( t ) ) C 1 i ( η ( t ) ) D 12 i ( η ( t ) ) Δ A i ( η ( t ) ) Δ B 1 i ( η ( t ) ) Δ B 2 i ( η ( t ) ) Δ C 1 i ( η ( t ) ) and Δ D 12 i ( η ( t ) ) are of appropriate dimensions. { ( η ( t ) ) } is a continuous-time discrete-state Markov process taking values in a finite set S = {1,2, , s} with transition probability matrix P r Δ _ _ { P ι k ( t ) } given by

P ι k ( t ) = P r ( η ( t + Δ ) = k | η ( t ) = ι ) = { λ ι k Δ + O ( Δ ) i f ι k 1 + λ ι ι Δ + O ( Δ ) i f ι = k E2

where Δ > 0, and lim Δ 0 O ( Δ ) Δ = 0 . Here λ ι k 0 is the transition rate from mode ι (system operating mode) to mode k ( ι k ) , and

λ ι ι = k = 1, k ι s λ ι k E3

For the convenience of notations, we let μ i Δ _ _ μ i ( υ ( t ) ) η = η ( t ) and any matrix M ( μ ι ) Δ _ _ M ( μ η = ι ) . The matrix functions Δ A i ( η ) Δ B 1 i ( η ) Δ B 2 i ( η ) Δ C 1 i ( η ) and Δ D 12 i ( η ) re-present the time-varying uncertainties in the system and satisfy the following assumption.

Assumption 1:

Δ A i ( η ) = F ( x ( t ) η t ) H 1 i ( η ) Δ B 1 i ( η ) = F ( x ( t ) η t ) H 2 i ( η ) Δ B 2 i ( η ) = F ( x ( t ) η t ) H 3 i ( η ) Δ C 1 i ( η ) = F ( x ( t ) η t ) H 4 i ( η )

and Δ D 12 i ( η ) = F ( x ( t ) η t ) H 5 i ( η )

where H j i ( η ) j = 1,2, ,5 are known matrices which characterize the structure of the uncertainties. Furthermore, there exists a positive function ρ(η) such that the following inequality holds:

F ( x ( t ) η t ) ρ ( η ) E4

We recall the following definition.

Definition 1: Suppose γ is a given positive number. A system of the form (1) is said to have the L 2-gain less than or equal to γ if

Ε [ 0 T f { z T ( t ) z ( t ) γ 2 w T ( t ) w ( t ) } d t ] 0 x ( 0 ) = 0 E5

where Ε [ ] stands for the mathematical expectation, for all Tf and all w(t) ∈ L 2 [0, Tf].

Note that for the symmetric block matrices, we use (*) as an ellipsis for terms that are induced by symmetry.

Advertisement

3. Robust H Fuzzy Control Design

This section provides the LMI-based solutions to the problem of designing a robust H fuzzy controller that guarantees the L 2-gain of the mapping from the exogenous input noise to the regulated output to be less than some prescribed value.

First, we consider the following H∞ fuzzy state feedback which is inferred as the weighted average of the local models of the form:

u ( t ) = j = 1 r μ j K j ( l ) x ( t ) E6

Then, we describe the problem under our study as follows.

Problem Formulation: Given a prescribed H performance γ > 0, design a robust H fuzzy state-feedback controller of the form (6) such that the inequality (5) holds.

Before presenting our first main result, the following lemma is needed.

Lemma 1: Consider the system (1). Given a prescribed H performance γ > 0, for ι = 1, 2, , s, if there exist matrices P ε ( ι ) = P ε T ( ι ) , positive constants δ(ι) and matrices Yj (ι), j = 1, 2, , r such that the following ε-dependent linear matrix inequalities hold:

P ε ( ι ) 0 E7
ψ i i ( ι ε ) 0, i = 1,2, r E8
ψ i j ( ι ε ) + ψ j i ( ι ε ) 0, i j r E9

where

Ψ i j ( ι ε ) = ( Φ i j ( ι ε ) ( * ) T ( * ) T ( * ) T R ( ι ) B 1 i T ( ι ) γ R ( ι ) ( * ) T ( * ) T φ i j ( ι ε ) 0 γ R ( ι ) ( * ) T Z T ( ι ε ) 0 0 P ( ι , ε ) ) E10
Φ i j ( ι ε ) = A i ( ι ) E ε 1 P ε ( ι ) + E ε 1 P ε ( ι ) A i T ( ι ) + B 2 i ( ι ) Y j ( ι ) + Y j T ( ι ) B 2 i T ( ι ) + λ ι ι E ε 1 P ε ( ι ) E11
ϕ i j ( ι ε ) = C ˜ 1 i ( ι ) E ε 1 P ε ( ι ) + D ˜ 12 i ( ι ) Y j ( ι ) E12
R ( ι ) = d i a g { δ ( ι ) I I δ ( ι ) I I } E13
Z ( ι , ε ) = ( λ ι 1 E ε 1 P ε ( ι ) λ ι ( ι 1 ) E ε 1 P ε ( ι ) λ ι ( ι + 1 ) E ε 1 P ε ( ι ) λ ι s E ε 1 P ε ( ι ) ) E14
P ( ι , ε ) = d i a g { E ε 1 P ε ( 1 ) E ε 1 P ε ( ι 1 ) E ε 1 P ε ( ι + 1 ) E ε 1 P ε ( s ) } E15

with

B ˜ 1 i ( ι ) = [ I I I B 1 i ( ι ) ] C ˜ 1 i ( ι ) = [ γ ρ ( ι ) H 1 i T ( ι ) 2 ( ι ) ρ ( ι ) H 4 i T ( ι ) 0 2 ( ι ) C 1 i T ( ι ) ] T D ˜ 12 i ( ι ) = [ 0 2 ( ι ) ρ ( ι ) H 5 i T ( ι ) γ ρ ( ι ) H 3 i T ( ι ) 2 ( ι ) D 12 i T ( ι ) ] T ( ι ) = ( I + ρ 2 ( ι ) i = 1 r j = 1 r [ H 2 i T ( ι ) H 2 j ( ι ) ] ) 1 2

then the inequality (5) holds. Furthermore, a suitable choice of the fuzzy controller is

u ( t ) = j = 1 r μ j K ε j ( ι ) x ( t ) E16

where

K ε j ( ι ) = Y j ( ι ) ( P ε ( ι ) ) 1 E ε E17

Proof: The closed-loop state space form of the fuzzy system model (1) with the controller (16) is given by

E ε x ˙ ( t ) = i = 1 r j = 1 r μ i μ j ( [ A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ] x ( t ) + [ Δ A i ( ι ) + Δ B 2 i ( ι ) K ε j ( ι ) ] x ( t ) + [ B 1 i ( ι ) + Δ B 1 i ( ι ) ] w ( t ) ) x ( 0 ) = 0 E18

or in a more compact form

E ε x ˙ ( t ) = i = 1 r j = 1 r μ i μ j ( [ A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ] x ( t ) + B ˜ 1 i ( ι ) R ( ι ) w ˜ ( t ) ) x ( 0 ) = 0 E19

where

B ˜ 1 i ( ι ) = [ I I I B 1 i ( ι ) ] E20
w ˜ ( t ) = R - 1 ( ι ) [ F ( x ( t ) ι t ) H 1 i ( ι ) x ( t ) F ( x ( t ) ι t ) H 2 i ( ι ) w ( t ) F ( x ( t ) ι t ) H 3 i ( ι ) K ε j ( ι ) x ( t ) w ( t ) ] E21

Consider a Lyapunov functional candidate as follows:

V ( x ( t ) ι ε ) = γ x T ( t ) Q ε ( ι ) x ( t ) ι S E22

Note that Qε (ι) is constant for eachι. For this choice, we have V(0, ι 0 , ε) = 0 and V(x(t),ι,ε)→∞ only when ||x(t)|| → ∞.

Now let consider the weak infinitesimal operator Δ ˜ of the joint process {(x(t), ı, ε), t ≥ 0}, which is the stochastic analog of the deterministic derivative. {(x(t), ı, ε), t ≥ 0} is a Markov process with infinitesimal operator given by (Souza & Fragoso, 1993),

Δ ˜ V ( x ( t ) ι ε ) = γ x ˙ T ( t ) Q ε ( ι ) x ( t ) + γ x T ( t ) Q ε ( ι ) x ˙ ( t ) + γ x T ( t ) k = 1 s λ ι k Q ε ( k ) x ( t ) = i = 1 r j = 1 r μ i μ j ( γ x T ( t ) Q ε ( ι ) [ ( A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ) ] x ( t ) + γ x T ( t ) [ A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ] T Q ε ( ι ) x ( t ) + γ x T ( t ) Q ε ( ι ) B ˜ 1 i ( ι ) R ( ι ) w ˜ ( t ) + γ w ˜ T ( t ) R ( ι ) B ˜ 1 i T ( ι ) Q ε ( ι ) x ( t ) + γ x T ( t ) k = 1 s λ ι k Q ε ( k ) x ( t ) ) E23

Adding and subtracting 2 ( ι ) z T ( t ) z ( t ) + γ 2 i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n [ w ˜ T ( t ) R ( ι ) w ˜ ( t ) ] to and from (23), we get

Δ ˜ V ( x ( t ) ι ε ) = 2 ( ι ) z T ( t ) z ( t ) + γ 2 i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n [ w ˜ T ( t ) R ( ι ) w ˜ ( t ) ] + 2 ( ι ) z T ( t ) z ( t ) + γ 2 i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n [ x ( t ) w ˜ ( t ) ] T × ( ( [ A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ] T Q ε ( ι ) Q ε ( ι ) [ A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ] + k = 1 s λ ι k Q ε ( k ) ) ( * ) T R ( ι ) B ˜ 1 i T ( ι ) Q ε ( ι ) γ R ( ι ) ) [ x ( t ) w ˜ ( t ) ] E24

Now let us consider the following terms:

γ 2 i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n [ w ˜ T ( t ) R ( ι ) w ˜ ( t ) ] ρ 2 ( ι ) γ 2 δ ( ι ) i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n x T ( t ) × { H 1 i T ( ι ) H 1 m ( ι ) + K ε j T ( ι ) H 3 i T ( ι ) H 3 m ( ι ) K ε n ( ι ) } x ( t ) + 2 ( ι ) γ 2 w T ( t ) w ( t ) E25

and

2 ( ι ) z T ( t ) z ( t ) 2 2 ( ι ) i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n x T ( t ) { [ C 1 i ( ι ) + D 12 i ( ι ) K ε j ( ι ) ] T × [ C 1 m ( ι ) + D 12 m ( ι ) K ε n ( ι ) ] + ρ 2 ( ι ) [ H 4 i ( ι ) + H 5 i ( ι ) K ε j ( ι ) ] T × [ H 4 m ( ι ) + H 5 m ( ι ) K ε n ( ι ) ] } x ( t ) E26

where ( ι ) = ( 1 + ρ 2 ( ι ) i = 1 r j = 1 r [ H 2 i T ( ι ) H 2 j ( ι ) ] ) 1 2 Hence,

γ 2 i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n [ w ˜ T ( t ) R ( ι ) w ˜ ( t ) ] + 2 ( ι ) z T ( t ) z ( t ) i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n ( x T ( t ) [ C ˜ 1 i ( ι ) + D ˜ 12 i ( ι ) K ε j ( ι ) ] T R - 1 ( ι ) × [ C ˜ 1 m ( ι ) + D ˜ 12 m ( ι ) K ε n ( ι ) ] x ( t ) + 2 ( ι ) γ 2 w T ( t ) w ( t ) E27

where

C ˜ 1 i ( ι ) = [ γ ρ ( ι ) H 1 i T ( ι ) 2 ( ι ) ρ ( ι ) H 4 i T ( ι ) 0 2 ( ι ) C 1 i T ( ι ) ] T D ˜ 12 i ( ι ) = [ 0 2 ( ι ) ρ ( ι ) H 5 i T ( ι ) γ ρ ( ι ) H 3 i T ( ι ) 2 ( ι ) D 12 i T ( ι ) ] T

Substituting (27) into (24), we have

Δ ˜ V ( x ( t ) ι ε ) 2 ( ι ) z T ( t ) z ( t ) + γ 2 2 ( ι ) w T ( t ) w ( t ) + γ 2 i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n × [ x ( t ) w ˜ ( t ) ] T Φ i j m n ( ι ε ) [ x ( t ) w ˜ ( t ) ] E28

where

Φ i j m n ( ι ε ) = ( ( ( A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ) T Q ε ( ι ) + Q ε ( ι ) [ A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ] + 1 γ [ C ˜ 1 i ( ι ) + D ˜ 12 i ( ι ) K ε j ( ι ) ] T × R - 1 ( ι ) [ C ˜ 1 m ( ι ) + D ˜ 12 m ( ι ) K ε n ( ι ) ] + k = 1 s λ ι k Q ε ( k ) ) ( * ) T R ( ι ) B ˜ 1 i T ( ι ) Q ε ( ι ) γ R ( ι ) )

Using the fact

i = 1 r j = 1 r m = 1 r n = 1 r μ i μ j μ m μ n M i j T ( ι ) N m n ( ι ) 1 2 i = 1 r j = 1 r μ i μ j [ M i j T ( ι ) M i j ( ι ) + N i j ( ι ) N i j T ( ι ) ]

we can rewrite (28) as follows:

Δ ˜ V ( x ( t ) ι ε ) 2 ( ι ) z T ( t ) z ( t ) + γ 2 2 ( ι ) w T ( t ) w ( t ) + γ 2 i = 1 r j = 1 r μ i μ j [ x ( t ) w ˜ ( t ) ] T Φ i j ( ι ε ) [ x ( t ) w ˜ ( t ) ] = 2 ( ι ) z T ( t ) z ( t ) + γ 2 2 ( ι ) w T ( t ) w ( t ) + γ i = 1 r μ i 2 [ x ( t ) w ˜ ( t ) ] T Φ i i ( ι ε ) [ x ( t ) w ˜ ( t ) ] + γ i = 1 r i j r μ i μ j [ x ( t ) w ˜ ( t ) ] T ( Φ i j ( ι ε ) + Φ i j ( ι ε ) ) [ x ( t ) w ˜ ( t ) ] E29

where

Φ i j ( ι ε ) = ( ( ( A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ) T Q ε ( ι ) + Q ε ( ι ) ( A i ( ι ) + B 2 i ( ι ) K ε j ( ι ) ) + 1 γ [ C ˜ 1 i ( ι ) + D ˜ 12 i ( ι ) K ε j ( ι ) ] T × R - 1 ( ι ) [ C ˜ 1 m ( ι ) + D ˜ 12 m ( ι ) K ε n ( ι ) ] + k = 1 s λ ι k Q ε ( k ) ) ( * ) T R ( ι ) B ˜ 1 i T ( ι ) Q ε ( ι ) γ R ( ι ) ) E30

Pre and post multiplying (30) by

Ξ ( ι ) = ( P ε ( ι ) 0 0 I )

with P ε ( ι ) = Q ε 1 ( ι ) we obtain

Ξ ( ι ) Φ i j ( ι ε ) Ξ ( ι ) = ( ( P ε ( ι ) A i T ( ι ) + Y j T ( ι ) B 2 i T ( ι ) + A i ( ι ) P ε ( ι ) + B 2 i ( ι ) Y j ( ι ) + 1 γ [ C ˜ 1 i ( ι ) P ε ( ι ) + D ˜ 12 i ( ι ) Y j ( ι ) ] T × R - 1 ( ι ) [ C ˜ 1 m ( ι ) P ε ( ι ) + D ˜ 12 m ( ι ) Y j ( ι ) ] + k = 1 s λ ι k P ε ( ι ) P ε 1 ( k ) P ε ( ι ) ) ( * ) T R ( ι ) B ˜ 1 i T ( ι ) γ R ( ι ) ) E31

Note that (31) is the Schur complement of Ψ ij (ι, ε) defined in (10). Using (8)-(9), we learn that

Φ i i ( ι ε ) 0, i = 1,2, r E32
Φ i j ( ι ε ) + Φ j i ( ι ε ) 0, i j r E33

Following from (29), (32) and (33), we know that

Δ ˜ V ( x ( t ) ι ε ) 2 ( ι ) z T ( t ) z ( t ) + γ 2 2 ( ι ) w T ( t ) w ( t ) E34

Applying the operator E [ 0 T f ( ) d t ] on both sides of (34), we obtain

E [ 0 T f Δ ˜ V ( x ( t ) ι ε ) d t ] E [ 0 T f ( 2 ( ι ) z T ( t ) z ( t ) + γ 2 2 ( ι ) w T ( t ) w ( t ) ) d t ] E35

From the Dynkin’s formula [2], it follows that

E [ 0 T f Δ ˜ V ( x ( t ) ι ε ) d t ] = E [ V ( x ( T f ) ι ( T f ) ε ) ] E [ V ( x ( 0 ) ι ( 0 ) ε ) ] E36

Substitute (36) into (35) yields

0 E [ 0 T f ( 2 ( ι ) z T ( t ) z ( t ) + γ 2 2 ( ι ) w T ( t ) w ( t ) ) d t ] E [ V ( x ( T f ) ι ( T f ) ε ) ] + E [ V ( x ( 0 ) ι ( 0 ) ε ) ]

Using (34) and the fact that we have

E [ 0 T f { z T ( t ) z ( t ) γ 2 w T ( t ) w ( t ) } d t ] 0 E37

Hence, the inequality (5) holds. This completes the proof of Lemma 1.

Remark 1: The linear matrix inequalities given in Lemma 1 becomes ill-conditioned when ε is sufficiently small, which is always the case for the uncertain nonlinear two time-scale dynamic systems. In general, these ill-conditioned linear matrix inequalities are very difficult to solve. Thus, to alleviate these ill-conditioned linear matrix inequalities, we have the following theorem which does not depend on ε.

Now we are in the position to present our first result.

Theorem 1: Consider the system (1). Given a prescribed H performance γ > 0, for ι = 1, 2, , s, if there exist matrices P(ι), positive constants δ(ι) and matrices Yj (ι), j = 1, 2, , r such that the following ε-independent linear matrix inequalities hold:

E P ( ι ) + P ( ι ) D 0 E38
Ψ i i ( ι ) 0 i = 1,2, r E39
Ψ i j ( ι ) + Ψ j i ( ι ) 0 i j r E40

where E P ( ι ) = P T ( ι ) E P ( ι ) D = D P T ( ι ) E = ( I 0 0 0 ) D = ( 0 0 0 I )

Ψ i j ( ι ) = ( Φ i j ( ι ) ( * ) T ( * ) T ( * ) T R ( ι ) B 1 i T ( ι ) γ R ( ι ) ( * ) T ( * ) T φ i j ( ι ) 0 γ R ( ι ) ( * ) T Z T ( ι ) 0 0 P ( ι ) ) E41
Φ i j ( ι ) = A i ( ι ) P ( ι ) + P T ( ι ) A i T ( ι ) + B 2 i ( ι ) Y j ( ι ) + Y j T ( ι ) B 2 i T ( ι ) + λ ι ι P ¯ ˜ ( ι ) E42
ϕ i j ( ι ) = C ˜ 1 i ( ι ) P ( ι ) + D ˜ 12 i ( ι ) Y j ( ι ) E43
R ( ι ) = d i a g { δ ( ι ) I I δ ( ι ) I I } E44
Z ( ι ) = ( λ ι 1 P ¯ ˜ ( ι ) λ ι ( ι 1 ) P ¯ ˜ ( ι ) λ ι ( ι + 1 ) P ¯ ˜ ( ι ) λ ι 1 P ¯ ˜ ( ι ) ) E45
P ( ι ) = d i a g { P ¯ ˜ ( 1 ) P ¯ ˜ ( ι 1 ) P ¯ ˜ ( ι + 1 ) P ¯ ˜ ( s ) } E46
P ¯ ˜ ( ι ) = P ( ι ) + P T ( ι ) 2 E47

with

B ˜ 1 i ( ι ) = [ I I I B 1 i ( ι ) ] C ˜ 1 i ( ι ) = [ γ ρ ( ι ) H 1 i T ( ι ) 2 ( ι ) ρ ( ι ) H 4 i T ( ι ) 0 2 ( ι ) C 1 i T ( ι ) ] T D ˜ 12 i ( ι ) = [ 0 2 ( ι ) ρ ( ι ) H 5 i T ( ι ) γ ρ ( ι ) H 3 i T ( ι ) 2 ( ι ) D 12 i T ( ι ) ] T ( ι ) = ( I + ρ 2 ( ι ) i = 1 r j = 1 r [ H 2 i T ( ι ) H 2 j ( ι ) ] ) 1 2

then there exists a sufficiently small ε ^ 0 such that the inequality (5) holds for ε ( 0, ε ] Furthermore, a suitable choice of the fuzzy controller is

u ( t ) = j = 1 r μ j K j ( ι ) x ( t ) E48

where

K j ( ι ) = Y j ( ι ) ( P ( ι ) ) 1 E49

Proof: Suppose there exists a matrix P(ι) such that the inequality (38) holds, then P(ι) is of the following form:

P ( ι ) = ( P 1 ( ι ) 0 P 2 T ( ι ) P 3 ( ι ) ) E50

with P 1 ( ι ) = P 1 T ( ι ) 0 and P 3 ( ι ) = P 3 T ( ι ) 0 . Let

P ε ( ι ) = E ε ( P ( ι ) + ε P ˜ ( ι ) ) E51

with

P ~ ( ι ) = ( 0 P 2 ( l ) 0 0 ) E52

Substituting (50) and (52) into (51), we have

P ε ( ι ) = ( P 1 ( ι ) ε P 2 ( ι ) ε P 2 T ( ι ) ε P 3 ( ι ) ) E53

Clearly, P ε ( ι ) = P ε T ( ι ) , and there exists a sufficiently small ε ^ such that for ε ( 0, ε ^ ] P ε ( ι ) 0 .

Using the matrix inversion lemma, we learn that

P ε 1 ( ι ) = [ P 1 ( ι ) + ε M ε ( ι ) ] E ε 1 E54

where M ε ( ι ) = P 1 ( ι ) P ˜ ( ι ) ( I + ε P 1 ( ι ) P ˜ ( ι ) ) 1 P 1 ( ι ) . Substituting (51) and (54) into (10), we obtain

Ψ i j ( ι ) + ψ i j ( ι ) E55

where the ε-independent linear matrix Ψ ij (ι) is defined in (41) and the ε-dependent linear matrix is

Ψ i j ( ι ) = ε ( ( A i ( ι ) P ˜ ( ι ) + P ˜ T ( ι ) A i T ( ι ) + B 2 i ( ι ) Y ε j ( ι ) + Y ε j T ( ι ) B 2 i T ( ι ) + λ ι ι P ¯ ˜ ( ι ) ) ( * ) T ( * ) T ( * ) T 0 0 ( * ) T ( * ) T ( C ˜ 1 i ( ι ) P ˜ ( ι ) + D ˜ 12 i ( ι ) Y ε j ( ι ) ) T 0 0 ( * ) T Z ε T ( ι ) 0 0 P ε ( ι ) ) E56

with Z ε ( ι ) = ( λ ι 1 P ¯ ^ ( ι ) λ ι ( ι 1 ) P ¯ ^ ( ι ) λ ι ( ι + 1 ) P ¯ ^ ( ι ) λ ι 1 P ¯ ˜ ( ι ) ) P ε ( ι ) = d i a g { P ¯ ^ ( 1 ) P ¯ ^ ( ι 1 ) P ¯ ^ ( ι + 1 ) P ¯ ^ ( s ) } P ¯ ^ ( ι ) = P ˜ ( ι ) + P ˜ T ( ι ) 2 and Y ε j ( ι ) = K j ( ι ) M ε 1 ( ι ) . Note that the ε- dependent linear matrix tends to zero when ε approaches zero.

Employing (38)-(40) and knowing the fact that for any given negative definite matrix W, there exists an ε > 0 such that W + εI < 0, one can show that there exists a sufficiently small ε ^ 0 such that for ε ( 0, ε ^ ] (8) and (9) hold. Since (7)-(9) hold, using Lemma 1, the inequality (5) holds for ε ( 0, ε ^ ] .

Advertisement

4. Illustrative Example

Consider a modified series dc motor model based on (Mehta & Chiasson, 1998) as shown in Figure 1 which is governed by the following difference equations:

J d ω ˜ ( t ) d t = K m L f i ˜ 2 ( t ) ( D + Δ D ) ω ˜ ( t ) L d i ˜ ( t ) d t = R i ˜ ( t ) K m L f i ˜ ( t ) ω ˜ ( t ) + V ˜ ( t ) E57

where ω ˜ ( t ) = ω ( t ) ω r e f ( t ) is the deviation of the actual angular velocity from the desired angular velocity, i ˜ ( t ) = i ( t ) i r e f ( t ) is the deviation of the actual current from the desired current, V ˜ ( t ) = V ( t ) V r e f ( t ) is the deviation of the actual input voltage from the desired input voltage, J is the moment of inertia, Km is the torque/back emf constant, D is the viscous friction coefficient, and R a , R f , L a and L f are the armature resistance, the field winding resistance, the armature inductance and the field winding inductance, respectively, with R Δ _ _ R f + R a and L Δ _ _ L f + L a . Note that in a typical series-connected dc motor, the condition L f L a holds. When one obtains a series connected dc motor, i ( t ) = i a ( t ) = i f ( t ) we have L f L a . Now let us assume that |ΔJ| ≤ 0.1J.

Figure 1.

A modified series dc motor equivalent circuit.

Giving x 1 ( t ) = ω ˜ ( t ) , x 2 ( t ) = i ˜ ( t ) and u ( t ) = V ˜ ( t ) , (57) becomes

[ x ˙ 1 ( t ) ε x ˙ 2 ( t ) ] = [ D J + Δ J K m L f J + Δ J x 2 ( t ) K m L f x 2 ( t ) R ] [ x 1 ( t ) x 2 ( t ) ] + [ 0 1 ] u ( t ) E58

where ε = L represents a small parasitic parameter. Assume that, the system is aggregated into 3 modes as shown in Table 1:

Mode Moment of Inertia J ( ) ± ∆ J ( )(kg· m 2 )
1 Small 0.0005 ± 10%
2 Normal 0.005 ± 10%
3 Large 0.05 ± 10%

Table 1.

System Terminology.

The transition probability matrix that relates the three operation modes is given as follows:

P ι k = [ 0.67 0.17 0.16 0.30 0.47 0.23 0.26 0.10 0.64 ]

The parameters for the system are given as R = 10 Ω, Lf = 0.005 H, D = 0.05 N m/rad/s and Km = 1 N m/A. Substituting the parameters into (58), we get

[ x ˙ 1 ( t ) ε x ˙ 2 ( t ) ] = [ 0.05 J ( ι ) 0.005 J ( ι ) x 2 ( t ) 0.005 x 2 ( t ) 10 ] [ x 1 ( t ) x 2 ( t ) ] + [ 0 0 0.1 0 ] w ( t ) + [ 0 1 ] u ( t ) + [ 0.05 Δ J ( ι ) 0.005 Δ J ( ι ) x 2 ( t ) 0 0 ] [ x 1 ( t ) x 2 ( t ) ] z ( t ) = [ 1 0 0 1 ] [ x 1 ( t ) x 2 ( t ) ] + [ 0 1 ] u ( t ) E59

where x ( t ) = [ x 1 T ( t ) x 2 T ( t ) ] T is the state variables, w ( t ) = [ w 1 T ( t ) w 2 T ( t ) ] T is the disturbance input, u(t) is the controlled input and z(t) is the controlled output.

The control objective is to control the state variable x2(t) for the range x2(t) ∈ [N1 N2]. For the sake of simplicity, we will use as few rules as possible.

Note that Figure 2 shows the plot of the membership function represented by

M 1 ( x 2 ( t ) ) = x 2 ( t ) + N 2 N 2 N 1 and M 2 ( x 2 ( t ) ) = x 2 ( t ) + N 2 N 2 N 1

Knowing that x2(t) ∈ [N1 N2], the nonlinear system

Figure 2.

Membership functions for the two fuzzy set.

(59) can be approximated by the following TS fuzzy model

E ε x ˙ ( t ) = i = 1 r μ i [ [ A i ( ι ) + Δ A i ( ι ) ] x ( t ) + B 1 i ω ( t ) + B 2 i u ( t ) ] x ( 0 ) = 0, z ( t ) = i = 1 r μ i [ C 1 i x ( t ) + D 12 i u ( t ) ]

where μi is the normalized time-varying fuzzy weighting functions for each rule, i = 1, 2, x(t)

= [ x 1 ( t ) x 2 ( t ) ] E ε = [ 1 0 0 ε ] Δ A 1 ( ι ) = F ( x ( t ) ι t ) H 1 1 ( ι ) Δ A 2 ( ι ) = F ( x ( t ) ι t ) H 1 2 ( ι )

A 1 ( 1 ) = [ 100 10 N 1 0.005 N 1 10 ] A 2 ( 1 ) = [ 100 10 N 2 0.005 N 2 10 ] A 1 ( 2 ) = [ 10 N 1 0.005 N 1 10 ] A 2 ( 2 ) = [ 10 N 2 0.005 N 2 10 ] A 1 ( 3 ) = [ 1 0.1 N 1 0.005 N 1 10 ] A 2 ( 3 ) = [ 1 0.1 N 2 0.005 N 2 10 ]

B 1 1 ( ι ) = B 1 2 ( ι ) = [ 0 0 0.1 0 ] B 2 1 ( ι ) = B 2 2 ( ι ) = [ 0 1 ] C 1 1 ( ι ) = C 1 2 ( ι ) = [ 1 0 0 1 ]

and D 12 1 ( ι ) = D 12 2 ( ι ) = [ 0 1 ]

with ||F(x(t),ι,t)|| ≤ 1. Then we have

H 1 1 ( ι ) = [ 0.05 J ( ι ) 0.05 J ( ι ) N 1 0 0 ] and H 1 2 ( ι ) = [ 0.05 J ( ι ) 0.05 J ( ι ) N 2 0 0 ]

In this simulation, we select N1 = −3 and N2 = 3. Using the LMI optimization algorithm and Theorem1 with ε = 0.005, γ = 1 and δ(1) = δ(2) = δ(3) = 1, we obtain the results given in Figure 3, Figure 4 and Figure 5

Remark 2: Employing results given in (Nguang & Shi, 2001; Han & Feng, 1998; Chen et al., 2000; Tanaka et al., 1996; Wang et al., 1996) and Matlab LMI solver [28], it is easy to realize that when ε< 0.005 for the state-feedback control design, LMIs become ill-conditioned and Matlab LMI solver yields an error message, “Rank Deficient”. However, the state-feedback fuzzy controller proposed in this paper guarantee that the inequality (5) holds for the system (59). Figure 3 shows the result of the changing between modes during the simulation with the initial mode at mode 1 and ε = 0.005. The disturbance input signal, w(t), which was used during simulation is given in Figure 4. The ratio of the regulated output energy to the disturbance input noise energy obtained by using the H fuzzy controller is depicted in Figure 5. The ratio of the regulated output energy to the disturbance input noise energy tends to a constant value which is about 0.0094. So γ = 0.0094 = 0.0970 which is less than the pres-cribed value 1. Finally, Table 2 shows the performance index, γ, for different values of ε.

Figure 3.

The result of the changing between modes during the simulation with the initial mode at mode.

Figure 4.

The disturbance input, w(t).

Figure 5.

The ratio of the regulated output energy to the disturbance noise energy, ( 0 T f z T ( t ) z ( t ) d t 0 T f w T ( t ) w ( t ) d t ) with ε = 0.005.

The performance index γ
ε State-feedback control design
0.005 0.0970
0.10 0.4796
0.30 0.8660
0.40 0.9945
0.41 > 1

Table 2.

The performance index γ for different values of ε.

Advertisement

5. Conclusion

This chapter has investigated the problem of designing a robust H controller for a class of uncertain nonlinear two time-scale dynamic systems with Markovian jumps that guarantees the L2 - gain from an exogenous input to a regulated output to be less or equal to a prescribed value. First, we approximate this class of uncertain nonlinear two time-scale dynamic systems with Markovian jumps by a class of uncertain Takagi-Sugeno fuzzy models with Markovian jumps. Then, based on an LMI approach, sufficient conditions for the uncertain nonlinear two time-scale dynamic systems with Markovian jumps to have an H performance are derived. The proposed approach does not involve the separation of states into slow and fast ones and it can be applied not only to standard, but also to nonstandard uncertain nonlinear two time-scale dynamic systems. An illustrative example is used to illustrate the effectiveness of the proposed design techniques.

References

  1. 1. Assawinchaichote W. Nguang S. K. 2002 Fuzzy control design for singularly perturbed nonlinear systems: An LMI approach, ICAIET, 146 151 , Kota Kinabalu, Malaysia.
  2. 2. Assawinchaichote W. Nguang S. K. 2002 Fuzzy observer-based controller design for singularly perturbed nonlinear systems: An LMI approach, Proc. IEEE Conf. Decision and Contr., 2165 2170 , Las Vegas.
  3. 3. Benjelloun K. Boukas E. K. Costa O. L. V. 1997 H control for linear time delay with Markovian jumping parameters, J. of Opt. Theory and App., 105 page (73-95).
  4. 4. Boukas E. K. Liu Z. K. 2001 Suboptimal design of regulators for jump linear system with timemultiplied quadratic cost, IEEE Tran. Automat. Contr., 46 page (944-949).
  5. 5. Boukas E. K. Yang H. 1999 Exponential stabilizability of stochastic systems with Markovian jump parameters, Automatica, 35 page (1437-1441).
  6. 6. Boukas E. K. Liu Z. K. Liu G. X. 2001 Delay dependent robust stability and H control of jump linear systems with time-delay,” Int. J. of Contr.
  7. 7. Boyd S. El Ghaoui L. Feron E. Balakrishnan V. 1994 Linear Matrix Inequalities in Systems and Control Theory, 15 SIAM, Philadelphia.
  8. 8. Chen B. S. Tseng C. S. He Y. Y. 2000 Mixed H/H fuzzy output feedback control design for nonlinear dynamic systems: An LMI approach, IEEE Trans. Fuzzy Syst., 8 page (249-265).
  9. 9. Dragan V. Shi P. Boukas E. K. 1999 Control of singularly perturbed system with Markovian jump parameters: An H∞ approach, Automatica, 35 page (1369-1378),.
  10. 10. Dynkin E. B. 1965 Markov Process, Springer-Verlag, Berlin.
  11. 11. Farias D. P. Geromel J. C. Val J. B. R. Costa O. L. V. 2000 Output feedback control of Markov jump linear systems in continuous-time, IEEE Trans. Automat. Contr., 45 page (944-949).
  12. 12. Feng X. Loparo K. A. Ji Y. Chizeck H. J. 1992 Stochastic stability properties of jump linear system, IEEE Tran. Automat. Contr., 37 page (38-53).
  13. 13. Fridman E. 2001 State-feedback H control of nonlinear singularly perturbed systems, Int. J. Robust Nonlinear Contr., 6 page (25-45).
  14. 14. Gahinet P. Nemirovski A. Laub A. J. Chilali M. 1995 LMI Control Toolbox- For Use with MATLAB, The MathWorks, Inc., MA.
  15. 15. Han Z. X. Feng G. 1998 State-feedback H∞ controller design of fuzzy dynamic system using LMI techniques, Fuzzy-IEEE’98, page (538-544).
  16. 16. Ikeda K. T. Wang H. O. 1996 Robust stabilization of a class of uncertain nonlinear systems via fuzzy control: Quadratic stability, H control theory, and linear matrix inequality, IEEE Trans. Fuzzy. Syst., 4 page (1-13).
  17. 17. Kokotovic P. V. Khalil H. K. O’Reilly J. 1986 Singular Perturbation Methods in Control: Analysis and Design, Academic Press, London.
  18. 18. Kushner H. J. 1967 Stochastic Stability and Control, Academic Press, New York.
  19. 19. Mehta S. Chiasson J. 1998 Nonlinear control of a series dc motor: Theory and experiment, IEEE Trans. Ind. Electron., 45 page (134-141).
  20. 20. Nguang S. K. Shi P. 2001 H fuzzy output feedback control design for nonlinear systems: An LMI approach, Proc. IEEE Conf. Decision and Contr., pages (4352-4357).
  21. 21. Nguang S. K. Shi P. 2003 H fuzzy output feedback control design for nonlinear systems: An LMI approach, IEEE Trans. Fuzzy Syst., 11 page (331-340).
  22. 22. Nguang S. K. 1996 Robust nonlinear H output feedback control, IEEE Trans Automat. Contr.,41 page (1003-1008).
  23. 23. Pan Z. Basar T. 1993 H-optimal control for singularly perturbed systems Part I: Perfect state measurements, Automatica, 29 page (401-423).
  24. 24. Pan Z. Basar T. 1994 H-optimal control for singularly perturbed systems Part II: Imperfect state measurements, IEEE Trans. Automat. Contr., 39 page(280-299).
  25. 25. Rami M. A. El Ghaoui L. 1995 H statefeedback control of jump linear systems, Proc.Conf. Decision and Contr., page (951-952),
  26. 26. Shi P. Boukas E. K. 1997 H control for Markovian jumping linear system with parametric uncertainty, J. of Opt. Theory and App., 95 page numbers (75-99),.
  27. 27. Shi P. Dragan V. 1999 Asymptotic H control of singularly perturbed system with parametric uncertainties, IEEE Trans. Automat. Contr., 44 page (1738-1742).
  28. 28. Souza C. E. Fragoso M. D. 1993 H control for linear systems with Markovian jumping parameters, Control-Theory and Advanced Tech., 9 page (457-466).
  29. 29. Wang H. Tanaka O. K. Griffin M. F. 1996 An approach to fuzzy control of nonlinear systems: Stability and design issues, IEEE Trans. Fuzzy Syst., 4 1 page (14-23).
  30. 30. Wonham W. M. 1968 On a matrix Riccati equation of stochastic control, SIAM J. Contr., 6 page (681-697).

Written By

Wudhichai Assawinchaichote, Sing Kiong Nguang and Non-members

Published: 01 October 2009