Open access peer-reviewed chapter

Non-Fragile Guaranteed Cost Control of Nonlinear Systems with Different State and Input Delays Based on T-S Fuzzy Local Bilinear Models

Written By

Junmin Li, Jinsha Li and Ruirui Duan

Submitted: 23 January 2017 Reviewed: 16 May 2017 Published: 30 August 2017

DOI: 10.5772/intechopen.69777

From the Edited Volume

Modern Fuzzy Control Systems and Its Applications

Edited by S. Ramakrishnan

Chapter metrics overview

1,318 Chapter Downloads

View Full Metrics

Abstract

This paper focuses on the non-fragile guaranteed cost control problem for a class of Takagi-Sugeno (T-S) fuzzy time-varying delay systems with local bilinear models and different state and input delays. A non-fragile guaranteed cost state-feedback controller is designed such that the closed-loop T-S fuzzy local bilinear control system is delay-dependent asymptotically stable, and the closed-loop fuzzy system performance is constrained to a certain upper bound when the additive controller gain perturbations exist. By employing the linear matrix inequality (LMI) technique, sufficient conditions are established for the existence of desired non-fragile guaranteed cost controllers. The simulation examples show that the proposed approach is effective and feasible.

Keywords

  • fuzzy control
  • non-fragile guaranteed cost control
  • delay-dependent
  • linear matrix inequality (LMI)
  • T-S fuzzy bilinear model

1. Introduction

In recent years, T-S (Takagi-Sugeno) model-based fuzzy control has attracted wide attention, essentially because the fuzzy model is an effective and flexible tool for the control of nonlinear systems [18]. Through the application of sector nonlinearity approach, local approximation in fuzzy partition spaces or other different approximation methods, T-S fuzzy models will be used to approximate or exactly represent a nonlinear system in a compact set of state variables. The merit of the model is that the consequent part of a fuzzy rule is a linear dynamic subsystem, which makes it possible to apply the classical and mature linear systems theory to nonlinear systems. Further, by using the fuzzy inference method, the overall fuzzy model will be obtained. A fuzzy controller is designed via the method titled ‘parallel distributed compensation (PDC)’ [36], the main idea of which is that for each linear subsystem, the corresponding linear controller is carried out. Finally, the overall nonlinear controller is obtained via fuzzy blending of each individual linear controller. Based on the above content,, T-S fuzzy model has been widely studied, and many results have been obtained [18]. In practical applications, time delay often occurs in many dynamic systems such as biological systems, network systems, etc. It is shown that the existence of delays usually becomes the source of instability and deteriorating performance of systems [38]. In general, when delay-dependent results were calculated, the emergence of the inner product between two vectors often makes the process of calculation more complicated. In order to avoid it, some model transformations were utilized in many papers, unfortunately, which will arouse the generation of an inequality, resulting in possible conservatism. On the other hand, due to the influence of many factors such as finite word length, truncation errors in numerical computation and electronic component parameter change, the parameters of the controller in a certain degree will change, which lead to imprecision in controller implementation. In this case, some small perturbations of the controllers’ coefficients will make the designed controllers sensitive, even worse, destabilize the closed-loop control system [9]. So the problem of non-fragile control has been important issues. Recently, the research of non-fragile control has been paid much attention, and a series of productions have been obtained [1013].

As we know, bilinear models have been widely used in many physical systems, biotechnology, socioeconomics and dynamical processes in other engineering fields [14, 15]. Bilinear model is a special nonlinear model, the nonlinear part of which consists of the bilinear function of the state and input. Compared with a linear model, the bilinear models have two main advantages. One is that the bilinear model can better approximate a nonlinear system. Another is that because of nonlinearity of it, many real physical processes may be appropriately modeled as bilinear systems. A famous example of a bilinear system is the population of biological species, which can be showed by d θ d t = θ v . In this equation, v is the birth rate minus death rate, and θ denotes the population. Obviously, the equation cannot be approximated by a linear model [14].

Most of the existing results focus on the stability analysis and synthesis based on T-S fuzzy model with linear local model. However, when a nonlinear system has of complex nonlinearities, the constructed T-S model will consist of a number of fuzzy local models. This will lead to very heavy computational burden. According to the advantages of bilinear systems and T-S fuzzy control, so many researchers paid their attentions to the T-S fuzzy models with bilinear rule consequence [1618]. From these papers, it is evident that the T-S fuzzy bilinear model may be suitable for some classes of nonlinear plants. In Ref. [16], a nonlinear system was transformed into a bilinear model via Taylor’s series expansion, and the stability of T-S fuzzy bilinear model was studied. Moreover, the result was stretched into the complex fuzzy system with state time delay [17]. Ref. [18] presented robust stabilization for a class of discrete-time fuzzy bilinear system. Very recently, a class of nonlinear systems is described by T-S fuzzy models with nonlinear local models in Ref. [19], and in this paper, the scholars put forward a new fuzzy control scheme with local nonlinear feedbacks, the advantage of which over the existing methods is that a fewer fuzzy rules and less computational burden. The non-fragile guaranteed cost controller was designed for a class of T-S discrete-time fuzzy bilinear systems in Ref. [20]. However, in Refs. [19, 20], the time-delay effects on the system is not considered. Ref. [17] is only considered the fuzzy system with the delay in the state and the derivatives of time-delay, d ˙ ( t ) < 1 is required. Refs. [2123] dealt with the uncertain fuzzy systems with time-delay in different ways. It should be pointed out that all the aforementioned works did not take into account the effect of the control input delays on the systems. The results therein are not applicable to systems with input delay. Recently, some controller design approaches have been presented for systems with input delay, see [2, 3, 4, 18, 2432] for fuzzy T-S systems and [8, 15, 33, 34] for non-fuzzy systems and the references therein. All of these results are required to know the exact delay values in the implementation. T-S fuzzy stochastic systems with state time-vary or distributed delays were studied in Refs. [3539]. The researches of fractional order T-S fuzzy systems on robust stability, stability analysis about “0 < α < 1”, and decentralized stabilization in multiple time delays were presented in Refs. [4042], respectively. For different delay types, the corresponding adaptive fuzzy controls for nonlinear systems were proposed in Refs. [33, 43, 44]. In Refs. [45, 46], to achieve small control amplitude, a new T-S fuzzy hyperbolic model was developed, moreover, Ref. [46] considered the input delay of the novel model. In Ref. [25, 47], the problems of observer-based fuzzy control design for T-S fuzzy systems were concerned.

So far, the problem of non-fragile guaranteed cost control for fuzzy system with local bilinear model with different time-varying state and input delays has not been discussed.

In this paper, the problem of delay-dependent non-fragile guaranteed cost control is studied for the fuzzy time-varying delay systems with local bilinear model and different state and input delays. Based on the PDC scheme, new delay-dependent stabilization conditions for the closed-loop fuzzy systems are derived. No model transformation is involved in the derivation. The merit of the proposed conditions lies in its reduced conservatism, which is achieved by circumventing the utilization of some bounding inequalities for the cross-product between two vectors as in Ref. [17]. The three main contributions of this paper are the following: (1) a non-fragile guaranteed cost controller is presented for the fuzzy system with time-varying delay in both state and input; (2) some free-weighting matrices are introduced in the derivation process, where the constraint of the derivatives of time-delay, d ˙ ( t ) < 1 and h ˙ ( t ) < 1 , is eliminated; and (3) the delay-dependent stability conditions for the fuzzy system are described by LMIs. Finally, simulation examples are given to illustrate the effectiveness of the obtained results.

The paper is organized as follows. Section 2 introduces the fuzzy delay system with local bilinear model, and non-fragile controller law for such system is designed based on the parallel distributed compensation approach in Section 3. Results of non-fragile guaranteed cost control are given in Section 4. Two simulation examples are used to illustrate the effectiveness of the proposed method in Section 5, which is followed by conclusions in Section 6.

Notation: Throughout this paper, the notation P > 0(P ≥ 0) stands for P being real symmetric and positive definite (or positive semi-definite). In symmetric block matrices, the asterisk (*) refers to a term that is induced by symmetry, and diag{….} denotes a block-diagonal matrix. The superscript T means matrix transposition. The notion i , j = 1 s is an abbreviation of i = 1 s j = 1 s . Matrices, if the dimensions are not explicitly stated, are assumed to have compatible dimensions for algebraic operations.

Advertisement

2. System description and assumptions

In this section, we introduce the T-S fuzzy time-delay system with local bilinear model. The ith rule of the fuzzy system is represented by the following form:

P l a n t   R u l e   i : I F ϑ 1 ( t )   i s   F i 1   a n d     a n d   ϑ v ( t )   i s   F i v T H E N   x ˙ ( t ) = A i x ( t ) + A d i x ( t d ( t ) ) + B i u ( t ) + B h i u ( t h ( t ) ) + N i x ( t ) u ( t ) + N d i x ( t d ( t ) ) u ( t h ( t ) )              x ( t ) = φ ( t ) ,     t [ τ 1 , 0 ]  ,  i = 1 , 2 , , s E1

where Fij is the fuzzy set, s is the number of fuzzy rules, x(t) ∈ Rn is the state vector, and u(t) ∈ R is the control input, ϑ1(t), ϑ2(t),…,ϑv(t) are the premise variables. It is assumed that the premise variables do not depend on the input u(t). A i , A d i , N i , N d i R n × n , B i , B h i R n × 1 denote the system matrices with appropriate dimensions. d(t) is a time-varying differentiable function that satisfies 0 ≤ d(t) ≤ τ1, 0 ≤ h(t) ≤ τ2, where τ1, τ2 are real positive constants as the upper bound of the time-varying delay. It is also assumed that d ˙ ( t ) σ 1 , h ˙ ( t ) σ 2 , and σ1, σ2 are known constants. The initial conditions φ(t), ϕ(t) are continuous functions of t, t [ - τ ,    0 ] , τ = min ( τ 1 , τ 2 ) .

Remark 1: The fuzzy system with time-varying state and input delays will be investigated in this paper, which is different from the system in Ref. [17]. In Ref. [17], only state time-varying delay is considered. And also, here, we assume that the derivative of time-varying delay is less than or equal to a known constant that may be greater than 1; the assumption on time-varying delay in Ref. [17] is relaxed.

By using singleton fuzzifier, product inferred and weighted defuzzifier, the fuzzy system can be expressed by the following globe model:

x ˙ ( t ) = i = 1 s h i ( ϑ ( t ) ) [ A i x ( t ) + A d i x ( t d ( t ) ) + B i u ( t ) + B h i u ( t h ( t ) ) + N i x ( t ) u ( t ) + N d i x ( t d ( t ) ) u ( t h ( t ) ) ] E2

where

h i ( ϑ ( t ) ) = ω i ( ϑ ( t ) ) / i = 1 s ω i ( ϑ ( t ) ) , ω i ( ϑ ( t ) ) = j = 1 v μ i j ( ϑ ( t ) ) , μij (ϑ(t)) is the grade of membership of ϑi(t) in Fij. In this paper, it is assumed that ω i ( ϑ ( t ) ) 0 , i = 1 s ω i ( ϑ ( t ) ) > 0 for all t. Then, we have the following conditions hi(ϑ(t)) ≥ 0, i = 1 s h i ( ϑ ( t ) ) = 1 for all t. In the consequent, we use abbreviation hi, hhi, xd(t), ud(t), xh(t), uh(t), to replace hi(ϑ(t)), hi(ϑ(th(t))), x(td(t)), u(td(t)), x(th(t)), u(th(t)), respectively, for convenience.

The objective of this paper is to design a state-feedback non-fragile guaranteed cost control law for the fuzzy system (2).

Advertisement

3. Non-fragile guaranteed cost controller design

Extending the design concept in Ref. [17], we give the following non-fragile fuzzy control law:

I F   ϑ 1 ( t )   i s   F 1 i   a n d     a n d   ϑ v ( t )   i s   F v i T H E N   u ( t ) = ρ ( K i + Δ K i ) x ( t ) 1 + x T ( K i + Δ K i ) T ( K i + Δ K i ) x = ρ sin θ i = ρ cos θ i ( K i + Δ K i ) x ( t ) E3

where ρ > 0 is a scalar to be assigned, and KiRl×n is a local controller gain to be determined. ΔKi represents the additive controller gain perturbations of the form ΔKi = HiFi(t)Eki with Hi and Eki being known constant matrices, and Fi(t) the uncertain parameter matrix satisfying F i T ( t ) F i ( t ) I . sin θ i = K ¯ i x ( t ) 1 + x T K ¯ i T K ¯ i x , cos θ i = 1 1 + x T K ¯ i T K ¯ i x , θ i [ π 2 , π 2 ] , K ¯ i = K i + Δ K i ( t ) = K i + H i F i ( t ) E k i .

The overall fuzzy control law can be represented by

u ( t ) = i = 1 s h i ρ K ¯ i x ( t ) 1 + x T K ¯ i T K ¯ i x = i = 1 s h i ρ sin θ i = i = 1 s h i ρ cos θ i K ¯ i x ( t ) E4

When there exists an input delay h(t), we have that

u h ( t ) = l = 1 s h h l ρ sin ϕ l = l = 1 s h h l ρ cos ϕ l K ˜ l x h ( t ) E5

where sin ϕ l = K ˜ l x h ( t ) 1 + x h T K ˜ l T K ˜ l x h , cos ϕ l = 1 1 + x h T K ˜ l T K ˜ l x h , ϕ l [ π 2 , π 2 ] , K ˜ l = K l + Δ K l ( t h ( t ) ) = K l + H l F l ( t h ( t ) ) E k l .

So, it is natural and necessary to make an assumption that the functions hi are well defined all t ∈ [−τ2, 0], and satisfy the following properties:

h i ( ϑ ( t h ( t ) ) ) 0 , for i = 1,2,…,s, and i = 1 s h i ( ϑ ( t h ( t ) ) ) = 1 .

By substituting Eq. (5) into Eq. (2), the closed-loop system can be given by

x ˙ ( t ) = i , j , l = 1 s h i h j h h l ( Λ i j x ( t ) + Λ d i j x d ( t ) + Λ h i l x h ( t ) ) E6

where

Λ i j = A i + ρ sin θ j N i + ρ cos θ j B i K ¯ j , Λ d i l = A d i + ρ sin ϕ l N d i , Λ h i l = ρ cos ϕ l B h i K ˜ l .

Given positive-definite symmetric matrices SRn×n and WR, we take the cost function

J = 0 [ x T ( t ) S x ( t ) + u T ( t ) W u ( t ) ] d t E7

Definition 1. The fuzzy non-fragile control law u(t) is said to be non-fragile guaranteed cost if for the system (2), there exist control laws (4) and (5) and a scalar J0 such that the closed-loop system (6) is asymptotically stable and the closed-loop value of the cost function (7) satisfies JJ0.

Advertisement

4. Analysis of stability for the closed-loop system

Firstly, the following lemmas are presented which will be used in the paper.

Lemma 1 [20]: Given any matrices M and N with appropriate dimensions such that ε > 0, we have MT N + NT Mε MT M + ε−1 NT N.

Lemma 2 [21]: Given constant matrices G, E and a symmetric constant matrix S of appropriate dimensions. The inequality S + GFE + ETFTGT < 0 holds, where F(t) satisfies FT(t) F(t) ≤ I if and only if, for some ε > 0, S + εGGT + ε−1ET E < 0.

The following theorem gives the sufficient conditions for the existence of the non-fragile guaranteed cost controller for system (6) with additive controller gain perturbations.

Theorem 1. Consider system (6) associated with cost function (7). For given scalars ρ > 0, τ1 > 0, τ2 > 0, σ1 > 0, σ2 > 0, if there exist matrices P > 0, Q1 > 0, Q2 > 0, R1 > 0, R2 > 0, Ki, i = 1, 2,…, s, X1, X2, X3, X4, Y1, Y2, Y3, Y4, and scalar ε > 0 satisfying the inequalities (8), the system (6) is asymptotically stable and the control law (5) is a fuzzy non-fragile guaranteed cost control law, moreover,

J x T ( 0 ) P x ( 0 ) + d ( 0 ) 0 x T ( s ) Q 1 x ( s ) d s + τ 1 0 θ 0 x ˙ T ( s ) R 1 x ˙ ( s ) d s d θ + h ( 0 ) 0 x T ( s ) Q 2 x ( s ) d s + τ 2 0 θ 0 x ˙ T ( s ) R 2 x ˙ ( s ) d s d θ = J 0
[ T i j l * * τ 1 X T τ 1 R 1 * τ 2 Z T 0 τ 2 R 2 ] < 0 ,     i , j , l = 1 , 2 , , s    E8

where T i j l = [ T 11 , i j * * * T 21 , i T 22 , i * * T 31 , i T 32 , i j T 33 , i l * T 41 , i T 42 , i T 43 T 44 ] ,   

T 11 , i j = Q 1 + Q 2 + X 1 + X 1 T + Y 1 A i + A i T Y 1 T + S + 2 ε ρ 2 Y 1 Y 1 T + 4 ε 1 ( N i T N i + ( B i K ¯ j ) T ( B i K ¯ j ) )           +   Z 1 T + Z 1 + ρ 2 K ¯ i T W K ¯ i , T 21 , i = X 1 T + X 2 + Z 2 + Y 2 A i + A d i T Y 1 T ,       T 31 , i = Z 3 Z 1 + X 3 + Y 3 A i , T 22 , i j = ( 1 σ 1 ) Q 1 X 2 X 2 T + Y 2 A d i + A d i T Y 2 T + 2 ε ρ 2 Y 2 Y 2 T + 4 ε 1 N d i T N d i , T 32 , i = X 3 + Y 3 A d i Z 2 T ,       T 33 , i l = ( 1 σ 2 ) Q 2 Z 3 Z 3 T + 2 ε ρ 2 Y 3 Y 3 T + 4 ε 1 ( B h i K ˜ l ) T B h i K ˜ l T 41 , i = P + X 4 + Z 4 + Y 4 A i Y 1 T , T 42 , i = X 4 + Y 4 A i Y 2 T , T 43 = Z 4 Y 3 T , T 44 = τ 1 R 1 + τ 2 R 2 Y 4 Y 4 T + 2 ε ρ 2 Y 4 Y 4 T .

Proof: Take a Lyapunov function candidate as

V ( x ( t ) , t ) = x T ( t ) P x ( t ) + t d ( t ) t x T ( s ) Q 1 x ( s ) d s + τ 1 0 t + θ t x ˙ T ( s ) R 1 x ˙ ( s ) d s d θ + t h ( t ) t x T ( s ) Q 2 x ( s ) d s + τ 2 0 t + θ t x ˙ T ( s ) R 2 x ˙ ( s ) d s d θ E9

The time derivatives of V(x(t),t), along the trajectory of the system (6), are given by

V ˙ ( x ( t ) , t ) = 2 x T ( t ) P x ˙ ( t ) + x T ( t ) ( Q 1 + Q 2 ) x ( t ) ( 1 d ˙ ( t ) ) x d T ( t ) Q 1 x d ( t ) + x ˙ T ( t ) ( τ 1 R 1 + τ 2 R 2 ) x ˙ ( t ) t τ 1 t x ˙ T ( s ) R 1 x ˙ ( s ) d s ( 1 h ˙ ( t ) ) x h T ( t ) Q 2 x h ( t ) t τ 2 t x ˙ T ( s ) R 2 x ˙ ( s ) d s E10

Define the free-weighting matrices as X = [ X 1 T X 2 T X 3 T X 4 T ] T , Y = [ Y 1 T Y 2 T Y 3 T Y 4 T ] T , Z = [ Z 1 T Z 2 T Z 3 T Z 4 T ] T , where XkRn×n, YkRn×n, ZkRn×n, k = 1, 2, 3, 4 will be determined later.

Using the Leibniz-Newton formula and system equation (6), we have the following identical equations:

[ x T ( t ) X 1 + x d T ( t ) X 2 + x h T ( t ) X 3 + x ˙ T ( t ) X 4 ] [ x ( t ) x d ( t ) t d ( t ) t x ˙ ( s ) d s ] 0 , [ x T ( t ) Z 1 + x d T ( t ) Z 2 + x h T ( t ) Z 3 + x ˙ T ( t ) Z 4 ] [ x ( t ) x h ( t ) t h ( t ) t x ˙ ( s ) d s ] 0 , i , j = 1 s h i h j h l [ x T ( t ) Y 1 + x d T ( t ) Y 2 + x h T ( t ) Y 4 + x ˙ T ( t ) Y 4 ] [ Λ i j x ( t ) + Λ d i l x d ( t ) + Λ h i l x h ( t ) x ˙ ( t ) ] 0 E11

Then, substituting Eq. (12) into Eq. (11) yields

V ˙ ( x ( t ) , t ) = 2 x T ( t ) P x ˙ ( t ) + x T ( t ) ( Q 1 + Q 2 ) x ( t ) + x ˙ T ( t ) ( τ 1 R 1 + τ 2 R 2 ) x ˙ ( t ) ( 1 d ˙ ( t ) ) x d T ( t ) Q 1 x d ( t ) ( 1 h ˙ ( t ) ) x h T ( t ) Q 2 x h ( t ) t τ 1 t x ˙ T ( s ) R 1 x ˙ ( s ) d s + 2 η T ( t ) X [ x ( t ) x d ( t ) t d ( t ) t x ˙ ( s ) d s ] t τ 2 t x ˙ T ( s ) R 2 x ˙ ( s ) d s + 2 η T ( t ) Z [ x ( t ) x h ( t ) t h ( t ) t x ˙ ( s ) d s ] + 2 η T ( t ) Y i , j , l = 1 s h i h j h h l [ Λ i j x ( t ) + Λ d i l x d ( t ) + Λ h i l x h ( t ) x ˙ ( t ) ] 2 x T ( t ) P x ˙ ( t ) + x T ( t ) ( Q 1 + Q 2 ) x ( t ) + x ˙ T ( t ) ( τ 1 R 1 + τ 2 R 2 ) x ˙ ( t ) ( 1 σ 1 ) x d T ( t ) Q 1 x d ( t ) ( 1 σ 2 ) x h T ( t ) Q 2 x h ( t ) t d ( t ) t x ˙ T ( s ) R 1 ( s ) x ˙ ( s ) d s + 2 η T ( t ) X [ x ( t ) x d ( t ) t d ( t ) t x ˙ ( s ) d s ] t h ( t ) t x ˙ T ( s ) R 2 ( s ) x ˙ ( s ) d s + 2 η T ( t ) Z [ x ( t ) x h ( t ) t h ( t ) t x ˙ ( s ) d s ] + 2 η T ( t ) Y i , j , l = 1 s h i h j h h l [ Λ i j x ( t ) + Λ d i l x d ( t ) + Λ h i l x h ( t ) x ˙ ( t ) ] + x T ( t ) S x ( t ) + i , j = 1 s h i h j ρ 2 x T ( t ) K ¯ i T cos θ i W K ¯ j cos θ j x ( t ) [ x T ( t ) S x ( t ) + u T ( t ) W u ( t ) ] E12

where η ( t ) = [ x T ( t ) ,   x d T ( t ) ,   x h T ( t ) ,   x ˙ T ( t ) ] T .

Applying Lemma 1, we have the following inequalities:

2 x T ( t ) Y 1 Λ i j x ( t ) 2 x T ( t ) Y 1 A i x ( t ) + ε ρ 2 x T ( t ) Y 1 Y 1 T x ( t ) + ε 1 x T ( t ) ( N i T N i + ( B i K ¯ j ) T ( B i K ¯ j ) ) x ( t ) , 2 x T ( t ) Y 1 Λ d i l x d ( t ) 2 x T ( t ) Y 1 A d i x d ( t ) + ε ρ 2 sin 2 φ l x T ( t ) Y 1 Y 1 T x ( t ) + ε 1 x d T ( t ) N d i T N d i x d ( t ) , 2 x T ( t ) Y 1 Λ h i l x h ( t ) ε ρ 2 cos 2 φ l x T ( t ) Y 1 Y 1 T x ( t ) + ε 1 x h T ( t ) ( B h i K ˜ l ) T ( B h i K ˜ l ) x h ( t ) , 2 x d T ( t ) Y 2 Λ i j x ( t ) 2 x d T ( t ) Y 2 A i x ( t ) + ε ρ 2 x d T ( t ) Y 2 Y 2 T x d ( t ) + ε 1 x T ( t ) ( N i T N i + ( B i K ¯ j ) T ( B i K ¯ j ) ) x ( t ) , 2 x T ( t ) Y 1 Λ h i l x h ( t ) ε ρ 2 cos 2 φ l x T ( t ) Y 1 Y 1 T x ( t ) + ε 1 x h T ( t ) ( B h i K ˜ l ) T ( B h i K ˜ l ) x h ( t ) , 2 x d T ( t ) Y 2 Λ i j x ( t ) 2 x d T ( t ) Y 2 A i x ( t ) + ε ρ 2 x d T ( t ) Y 2 Y 2 T x d ( t ) + ε 1 x T ( t ) ( N i T N i + ( B i K ¯ j ) T ( B i K ¯ j ) ) x ( t ) , 2 x d T ( t ) Y 2 Λ d i l x d ( t ) 2 x d T ( t ) Y 2 A d i x d ( t ) + ε ρ 2 sin 2 φ l x d T ( t ) Y 2 Y 2 T x d ( t ) + ε 1 x d T ( t ) N d i T N d i x d ( t ) , 2 x d T ( t ) Y 2 Λ h i l x h ( t ) ε ρ 2 cos 2 φ l x d T ( t ) Y 2 Y 2 T x d ( t ) + ε 1 x h T ( t ) ( B h i K ˜ l ) T ( B h i K ˜ l ) x h ( t ) , 2 x h T ( t ) Y 3 Λ i j x ( t ) 2 x h T ( t ) Y 3 A i x ( t ) + ε ρ 2 x h T ( t ) Y 3 Y 3 T x h ( t ) + ε 1 x T ( t ) ( N i T N i + ( B i K ¯ j ) T ( B i K ¯ j ) ) x ( t ) , 2 x h T ( t ) Y 3 Λ d i l x d ( t ) 2 x d T ( t ) Y 3 A d i x d ( t ) + ε ρ 2 sin 2 φ l x h T ( t ) Y 3 Y 3 T x h ( t ) + ε 1 x d T ( t ) N d i T N d i x d ( t ) , 2 x h T ( t ) Y 3 Λ h i l x h ( t ) ε ρ 2 cos 2 φ l x h T ( t ) Y 3 Y 3 T x h ( t ) + ε 1 x h T ( t ) ( B h i K ˜ l ) T ( B h i K ˜ l ) x h ( t ) , 2 x ˙ T ( t ) Y 4 Λ i j x ( t ) 2 x ˙ T ( t ) Y 4 A i x ( t ) + ε ρ 2 x ˙ T ( t ) Y 4 Y 4 T x ˙ ( t ) + ε 1 x T ( t ) ( N i T N i + ( B i K ¯ j ) T ( B i K ¯ j ) ) x ( t ) , 2 x ˙ T ( t ) Y 4 Λ d i l x d ( t ) 2 x ˙ T ( t ) Y 4 A d i x d ( t ) + ε ρ 2 x ˙ T ( t ) Y 4 Y 4 T x ˙ ( t ) + ε 1 x d T ( t ) N d i T N d i x d ( t ) , 2 x ˙ T ( t ) Y 4 Λ h i l x h ( t ) ε ρ 2 cos 2 φ l x ˙ T ( t ) Y 4 Y 4 T x ˙ ( t ) + ε 1 x h T ( t ) ( B h i K ˜ l ) T ( B h i K ˜ l ) x h ( t ) E13

Substituting Eq. (13) into Eq. (12) results in

V ˙ ( x ( t ) , t ) i , j , l = 1 s h i h j h h l η T ( t ) T i j η ( t ) t d ( t ) t x ˙ T ( s ) R 1 x ˙ ( s ) d s t h ( t ) t x ˙ T ( s ) R 2 x ˙ ( s ) d s 2 η T ( t ) X t d ( t ) t x ˙ ( s ) d s 2 η T ( t ) Z t h ( t ) t x ˙ ( s ) d s [ x T ( t ) S x ( t ) + u T ( t ) W u ( t ) ] i , j , l = 1 s h i h j h h l η T ( t ) ( T i j l + τ 1 X R 1 1 X T + τ 2 Z R 2 1 Z T ) η ( t ) t d ( t ) t ( η T ( t ) X + x ˙ T ( s ) R 1 ) R 1 1 ( η T ( t ) X + x ˙ T ( s ) R 1 ) T d s t h ( t ) t ( η T ( t ) Z + x ˙ T ( s ) R 2 ) R 2 1 ( η T ( t ) X + x ˙ T ( s ) R 2 ) T d s [ x T ( t ) S x ( t ) + u T ( t ) W u ( t ) ]         i , j , l = 1 s h i h j h h l η T ( t ) ( T ˜ i j l + τ 1 X R 1 1 X T + τ 2 Z R 2 1 Z T ) η ( t ) [ x T ( t ) S x ( t ) + u T ( t ) W u ( t ) ] E14

where

   T ˜ i j l = [ T ˜ 11 , i j * * * T 21 , i T 22 , i * * T 31 , i T 32 , i j T 33 , i l * T 41 , i T 42 , i T 43 T 44 ] ,    T ˜ 11 , i j = T 11 , i j + ρ 2 K ¯ i T cos θ i W K ¯ j cos θ j ρ 2 K ¯ i T W K ¯ j . E107

In light of the inequality K ¯ i T W K ¯ j + K ¯ j T W K ¯ i K ¯ i T W K ¯ i + K ¯ j T W K ¯ j , we have

V ˙ ( x ( t ) , t ) i , j , l = 1 s h i h j h h l η T ( t ) ( Τ i j l + τ 1 X R 1 1 X T + τ 2 Z R 2 1 Z T ) η ( t ) [ x T ( t ) S x ( t ) + u T ( t ) W u ( t ) ] E15

Applying the Schur complement to Eq. (8) yields

Τ i i + + τ 1 X R 1 1 X T + τ 2 Z R 2 1 Z T < 0 , Τ i j + Τ j i + 2 τ 1 X R 1 1 X T + 2 τ 2 Z R 2 1 Z T < 0. E1259

Therefore, it follows from Eq. (15) that

V ˙ ( x ( t ) , t ) [ x T ( t ) S x ( t ) + u T ( t ) W u ( t ) ] < 0 E16

which implies that the system (6) is asymptotically stable.

Integrating Eq. (16) from 0 to T produces

0 T [ x T ( t ) S x ( t ) + u T ( t ) W u ( t ) ] d t V ( x ( T ) , T ) + V ( x ( 0 ) , 0 ) < V ( x ( 0 ) , 0 ) E22651

Because of V (x(t),t) ≥ 0 and V ˙ ( x ( t ) , t ) < 0 , thus lim T V ( x ( T ) , T ) = c , where c is a nonnegative constant. Therefore, the following inequality can be obtained:

J x T ( 0 ) P x ( 0 ) + d ( 0 ) 0 x T ( s ) Q 1 x ( s ) d s + τ 1 0 θ 0 x ˙ T ( s ) R 1 x ˙ ( s ) d s d θ + h ( 0 ) 0 x T ( s ) Q 2 x ( s ) d s + τ 2 0 θ 0 x ˙ T ( s ) R 2 x ˙ ( s ) d s d θ = J 0 E17

This completes the proof.

Remark 2: In the derivation of Theorem 1, the free-weighting matrices XkRn×n, YkRn×n, k = 1, 2, 3, 4 are introduced, the purpose of which is to reduce conservatism in the existing delay-dependent stabilization conditions, see Ref. [17].

In the following section, we shall turn the conditions given in Theorem 1 into linear matrix inequalities (LMIs). Under the assumptions that Y1, Y2, Y3, Y4 are non-singular, we can define the matrix Y i T = λ Z , i = 1, 2, 3, 4, Z = P−1,λ > 0.

Pre- and post-multiply (8) and (9) with Θ = d i a g { Y 1 1 , Y 2 1 , Y 3 1 , Y 4 1 , Y 4 1 , Y 4 1 } and Θ T = d i a g { Y 1 T , Y 2 T , Y 3 T , Y 4 T , Y 4 T , Y 4 T } , respectively, and letting Q ¯ 1 = Y 1 1 Q 1 Y 1 T , Q ¯ 2 = Y 1 1 Q 2 Y 1 T , R ¯ k = Y 4 1 R k Y 4 T , k = 1 , 2 , X ¯ i = Y i 1 X i Y i T , Z ¯ i = Y i 1 Z i Y i T , i = 1, 2, 3, 4, we obtain the following inequality (18), which is equivalent to (8):

[ T ¯ 11 , i j * * * * * T ¯ 21 , i T ¯ 22 , i * * * * T ¯ 31 , i T ¯ 32 , i T ¯ 33 , i l * * * T ¯ 41 , i T ¯ 42 , i T ¯ 43 T ¯ 44 * * τ 1 X 1 τ 1 X 2 τ 1 X 3 τ 1 X 4 τ 1 R ¯ 1 * τ 2 Z ¯ 1 τ 2 Z ¯ 2 τ 2 Z ¯ 3 τ 2 Z ¯ 4 0 τ 2 R ¯ 2 ] < 0 ,      i , j , l = 1 , 2 , , s E18

where

Τ ¯ 11 , i j = Q ¯ 1 + Q ¯ 2 + X ¯ 1 + X ¯ 1 T + λ A i Z + λ Z A i T + λ 2 Z S Z + 2 ε ρ 2 I + 4 ε 1 λ 2 Z N i T N i Z + Z 1 + Z 1 T + 4 ε 1 λ 2 ( B i K ¯ j Z ) T ( B i K ¯ j Z ) + ρ 2 λ 2 Z K ¯ i T W K ¯ i Z , Τ ¯ 21 , i = X ¯ 1 T + X ¯ 2 + Z ¯ 2 + λ A i Z + λ Z A d i T ,     Τ ¯ 31 , i = Z ¯ 3 Z ¯ 1 + X ¯ 3 + λ A i Z , Τ ¯ 41 , i = λ 2 Z + λ A i Z λ Z + X ¯ 4 + Z ¯ 4 , Τ ¯ 22 , i = ( 1 σ 1 ) Q ¯ 1 X ¯ 2 X ¯ 2 T + λ A d i Z + λ Z A d i T + 2 ε ρ 2 I + 4 ε 1 λ 2 Z N d i T N d i Z , Τ ¯ 32 , i = X ¯ 3 Z ¯ 2 + λ A d i Z λ Z A d i T ,   Τ ¯ 42 , i X ¯ 4 + λ A i Z λ Z ,   Τ ¯ 33 , i l = ( 1 σ 2 ) Q ¯ 2 Z ¯ 3 Z ¯ 3 T + 4 ε 1 λ 2 ( B h i K ˜ l Z ) T B h i K ˜ l Z + 2 ε ρ 2 I , Τ ¯ 43 = Z ¯ 4 λ Z ,   Τ ¯ 44 = τ 1 R 1 + τ 2 R 2 λ Z λ Z T + 2 ε ρ 2 I . E2256314

Applying the Schur complement to Eq. (18) results in

Γ i j l = [ Φ 11 , i * * Φ 21 , i j Φ 22 * Φ 31 , i l 0 Φ 33 ] < 0 ,      i , j , l = 1 , 2 , , s E19

where

Φ 11 , i = [ T ¯ ¯ 11 , i * * * * * T ¯ 21 , i T ¯ ¯ 22 , i * * * * T ¯ 31 , i T ¯ 32 , i T ¯ ¯ 33 , i l * * * T ¯ 41 , i T ¯ 42 , i T ¯ 43 T ¯ 44 * * τ 1 X 1 τ 1 X 2 τ 1 X 3 τ 1 X 4 τ 1 R ¯ 1 * τ 2 Z ¯ 1 τ 2 Z ¯ 2 τ 2 Z ¯ 3 τ 2 Z ¯ 4 0 τ 2 R ¯ 2 ] E2589631426

With

Τ ¯ ¯ 11 , i = Q ¯ 1 + Q ¯ 2 + X ¯ 1 + X ¯ 1 T + λ A i Z + λ Z A i T + Z 1 + Z 1 T + 2 ε ρ 2 I , Τ ¯ ¯ 22 , i = ( 1 σ 1 ) Q ¯ 1 X ¯ 2 X ¯ 2 T + λ A d i Z + λ Z A d i T + 2 ε ρ 2 I , Τ ¯ ¯ 33 = ( 1 σ 2 ) Q ¯ 2 Z ¯ 3 Z ¯ 3 T + 2 ε ρ 2 I . E25896327
Φ 21 , i j = [ λ Z 0 0 0 λ N i Z 0 0 0 λ B i K ¯ j Z 0 0 0 ρ λ K ¯ j Z 0 0 0 ] Φ 31 , i l = [ 0 λ B h i K ˜ l Z 0 0 0 0 λ N d i Z 0 ] Φ 22 = [ S 1 0 0 0 0 ε 4 I 0 0 0 0 ε 4 I 0 0 0 0 W 1 ] Φ 33 = [ ε 4 I 0 0 ε 4 I ] E22365418

Obviously, the closed-loop fuzzy system (6) is asymptotically stable, if for some scalars λ > 0, there exist matrices Z > 0 , Q ¯ > 0 , R ¯ > 0 and X ¯ 1 , X ¯ 2 , X ¯ 3 , K ¯ i , i = 1 , 2 , .. , s satisfying the inequality (19).

Theorem 2. Consider the system (6) associated with cost function (7). For given scalars ρ > 0, τ1 > 0, τ2 > 0, σ1 > 0, σ2 > 0 and λ > 0, δ > 0, if there exist matrices Z > 0 , Q ¯ 1 > 0 , R ¯ 1 > 0 , Q ¯ 2 > 0 , R ¯ 2 > 0 and X ¯ 1 , X ¯ 2 , X ¯ 3 , X ¯ 4 , Mi, i = 1,2,…,s and scalar ε > 0 satisfying the following LMI (20), the system (6) is asymptotically stable and the control law (5) is a fuzzy non-fragile guaranteed cost control law

[ Θ 1 , i j l Θ 2 , i j l Θ 3 ] < 0 ,      i , j , l = 1 , 2 , , s E20

Moreover, the feedback gains are given by

K i = M i Z 1 , i = 1 , 2 , , s E325630

and

J x T ( 0 ) P x ( 0 ) + d ( 0 ) 0 x T ( s ) Q 1 x ( s ) d s + τ 1 0 θ 0 x ˙ T ( s ) R 1 x ˙ ( s ) d s d θ + h ( 0 ) 0 x T ( s ) Q 2 x ( s ) d s + τ 2 0 θ 0 x ˙ T ( s ) R 2 x ˙ ( s ) d s d θ = J 0 E336251

where

Θ 2 , i j l = [ λ E k j Z 0 0 0 0 0 0 0 0 0 0 0 λ E k i Z 0 0 0 0 0 0 0 0 0 0 0 0 0 λ E k l Z 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ( B i H j ) T 0 0 0 0 0 0 0 0 0 0 0 0 ρ H i T 0 0 0 0 0 0 0 0 0 0 0 0 0 ( B h i H l ) T ] , Θ 3 = diag { δ I , δ I , δ I , δ 1 I , δ 1 I , δ 1 I } , E32653562
Θ 1 , i j l = [ Τ ¯ ¯ 11 , i * * * * * * * * * * * Τ ¯ 21 , i Τ ¯ ¯ 22 , i * * * * * * * * * * Τ ¯ 31 , i Τ ¯ 32 , i Τ ¯ 33 * * * * * * * * * Τ ¯ 41 , i Τ ¯ 42 , i Τ ¯ 43 , i Τ ¯ 44 * * * * * * * * τ 1 X ¯ 1 τ 1 X ¯ 2 τ 1 X ¯ 3 τ 1 X ¯ 4 τ 1 R ¯ 1 * * * * * * * τ 2 Z 1 τ 2 Z 2 τ 2 Z 3 τ 2 Z 4 0 τ 2 R ¯ 2 1 * * * * * * λ Z 0 0 0 0 0 S 1 * * * * * λ N i Z 0 0 0 0 0 0 ε 4 I * * * * λ B i M j 0 0 0 0 0 0 0 ε 4 I * * * ρ λ M i 0 0 0 0 0 0 0 0 W 1 * * 0 λ N d i Z 0 0 0 0 0 0 0 0 ε 4 I * 0 0 λ B h i M l 0 0 0 0 0 0 0 0 ε 4 I ] E32563253

Proof: At first, we prove that the inequality (20) implies the inequality (19). Applying the Schur complement to Eq. (20) results in

Φ 1 , i j l + δ [ 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 B i H j 0 0 0 ρ H i 0 0 0 0 0 0 B h i H l ] [ 0 0 0 0 0 0 0 0 ( B i H j ) T 0 0 0 0 0 0 0 0 0 0 0 0 ρ H i T 0 0 0 0 0 0 0 0 0 0 0 0 0 ( B h i H i ) T ] E3256324
+ δ 1 [ ( λ E k j Z ) T ( λ E k i Z ) T 0 0 0 0 0 0 ( λ E k l Z ) T 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] [ λ E k j Z 0 0 0 0 0 0 0 0 0 0 0 λ E k i Z 0 0 0 0 0 0 0 0 0 0 0 0 0 λ E k l Z 0 0 0 0 0 0 0 0 0 ] < 0 E21

Using Lemma 2 and noting Mi = KiZ, by the condition (21), the following inequality holds:

Φ 1 , i j l + [ 0 * * * * * * * * * * * 0 0 * * * * * * * * * * 0 0 0 * * * * * * * * * 0 0 0 0 * * * * * * * * 0 0 0 0 0 * * * * * * * 0 0 0 0 0 0 * * * * * * 0 0 0 0 0 0 0 * * * * * 0 0 0 0 0 0 0 0 * * * * λ B i Δ K ¯ j Z 0 0 0 0 0 0 0 0 * * * ρ λ Δ K ¯ i Z 0 0 0 0 0 0 0 0 0 * * 0 0 0 0 0 0 0 0 0 0 0 * 0 0 λ B h i Δ K ˜ l Z 0 0 0 0 0 0 0 0 0 ] < 0 E22

where Δ K ˜ i = Δ K i ( t d ( t ) ) .

Therefore, it follows from Theorem 1 that the system (6) is asymptotically stable and the control law (5) is a fuzzy non-fragile guaranteed cost control law. Thus, we complete the proof.

Now consider the cost bound of

J x T ( 0 ) P x ( 0 ) + d ( 0 ) 0 x T ( s ) Q 1 x ( s ) d s + τ 1 0 θ 0 x ˙ T ( s ) R 1 x ˙ ( s ) d s d θ + h ( 0 ) 0 x T ( s ) Q 2 x ( s ) d s + τ 2 0 θ 0 x ˙ T ( s ) R 2 x ˙ ( s ) d s d θ = J 0 E3256327

Similar to Ref. [23], we supposed that there exist positive scalars α1, α2, α3, α4, α5, such that Z 1 α 1 I , 1 λ 2 P Q ¯ 1 P α 2 I , 1 λ 2 P Q ¯ 2 P α 3 I , 1 λ 2 P R ¯ 1 P α 4 I , 1 λ 2 P R ¯ 2 P α 5 I .

Then, define S Q 1 = Q ¯ 1 1 , S Q 2 = Q ¯ 2 1 , S R 1 = R ¯ 1 1 , S R 2 = R ¯ 2 1 , by Schur complement lemma, we have the following inequalities:

[ α 1 I I I Z ] 0 , [ α 2 I 1 λ P 1 λ P S Q 1 ] 0 , [ α 3 I 1 λ P 1 λ P S Q 2 ] 0 , [ α 4 I 1 λ P 1 λ P S R 1 ] 0 , [ α 5 I 1 λ P 1 λ P S R 2 ] 0 , [ Z I I Z ] 0 , [ S Q 1 I I Q ¯ 1 ] 0 , [ S Q 2 I I Q ¯ 2 ] 0 , [ S R 1 I I R ¯ 1 ] 0 , [ S R 2 I I R ¯ 2 ] 0 , E23

Using the idea of the cone complement linear algorithm in Ref. [24], we can obtain the solution of the minimization problem of upper bound of the value of the cost function as follows:

minimize { t r a c e ( P Z + S Q 1 Q ¯ 1 + S Q 2 Q ¯ 2 + S R 1 R ¯ 1 + S R 2 R ¯ 2 + α 1 x T ( 0 ) x ( 0 ) + α 2 d ( 0 ) 0 x T ( s ) x ( s ) d s + α 4 τ 1 0 θ 0 x ˙ T ( s ) x ˙ ( s ) d s d θ + α 3 h ( 0 ) 0 x T ( s ) x ( s ) d s + α 5 τ 2 0 θ 0 x ˙ T ( s ) x ˙ ( s ) d s d θ } subject to (20),(23),   ε >0, Q ¯ 1 > 0 , Q ¯ 2 > 0 , R ¯ 1 > 0 , R ¯ 2 > 0 , Z > 0 , α i > 0 , i = 1 , , 5 E24

Using the following cone complement linearization (CCL) algorithm [24] can iteratively solve the minimization problem (24). □

Advertisement

5. Simulation examples

In this section, the proposed approach is applied to the Van de Vusse system to verify its effectiveness.

Example: Consider the dynamics of an isothermal continuous stirred tank reactor for the Van de Vusse

x ˙ 1 = 50 x 1 10 x 1 3 + u ( 10 x 1 ) + u ( t h ) + u ( t h ) ( 0.5 x 1 ( t d ) + 0.2 x 2 ( t d ) ) + 5 x 2 ( t d ) x ˙ 2 = 50 x 1 100 x 2 u ( t h ) + u ( t h ) ( 0.3 x 1 ( t d ) 0.2 x 2 ( t d ) ) + 10 x 2 ( t d ) 5 x 1 ( t d ) E25

From the system equation (25), some equilibrium points are tabulated in Table 1. According to these equilibrium points, [xe ue], which are also chosen as the desired operating points, [ x e   u e ] , we can use the similar modeling method that is described in Ref. [16].

x e T x d e T ue ude
[2.0422 1.2178] [2.0422 1.2178] 20.3077 20.3077
[3.6626 2.5443] [3.6626 2.5443] 77.7272 77.7272
[5.9543 5.5403] [5.9543 5.5403] 296.2414 296.2414

Table 1.

Data for equilibrium points.

Thus, the system (25) can be represented by

R 1 :   i f   x 1   i s   a b o u t   2.0422 then

x ˙ δ ( t ) = A 1 x δ ( t ) + A d 1 x d δ ( t ) + B 1 u δ ( t ) + B h 1 u d δ ( t ) + N 1 x δ ( t ) u δ ( t ) + N d 1 x d δ ( t ) u h δ ( t ) E41
R 2 :   i f   x 1   i s   a b o u t   3.6626 ,   t h e n            x ˙ δ ( t ) = A 2 x δ ( t ) + A d 2 x d δ ( t ) + B 2 u δ ( t ) + B d 2 u h δ ( t ) + N 2 x δ ( t ) u δ ( t ) + N d 2 x d δ ( t ) u h δ ( t ) R 3 :   i f   x 1   i s   a b o u t   5.9543 , t h e n       x ˙ δ ( t ) = A 3 x δ ( t ) + A d 3 x d δ ( t ) + B 3 u δ ( t ) + B d 3 u h δ ( t ) + N 3 x δ ( t ) u δ ( t ) + N d 3 x d δ ( t ) u h δ ( t ) E26

where

A 1 = [ 75.2383 7.7946 50 100 ] , A 2 = [ 98 .3005 11.7315 50 100 ] , A 3 = [ 122.1228 8.8577 50 100 ] , N 1 = N 2 = N 3 = [ 1 0 0 1 ] ; B 1 = B 2 = B 3 = [ 10 0 ] ; A d 1 = A d 2 = A d 3 = [ 0 5 10 5 ] , N d 1 = N d 2 = N d 3 = [ 0.5 0.2 0.3 0.2 ] , B h 1 = B h 2 = B h 3 = [ 1 0 ] , x δ = x ( t ) x e ' , u δ = u ( t ) u e ' , x d δ = x ( t d ) x d e ' , u h δ = u ( t d ) u h e ' . E43623

The cost function associated with this system is given with S = [ 1 0 0 1 ] , W = 1 . The controller gain perturbation ΔK of the additive form is give with H1 = H2 = H3 = 0.1, Ek1 = [0.05 −0.01], Ek2 = [0.02 0.01], Ek3 = [−0.01 0].

The membership functions of state x1 are shown in Figure 1.

Figure 1.

Membership functions.

Then, solving LMIs (23) and (24) for ρ = 0.45, λ = 1.02 and δ = 0.11 , τ 1 = τ 2 = 2 , σ 1 = 0 , σ 2 = 0 gives the following feasible solution:

P = [ 4.2727 1.3007 1.3007 6.4906 ] , Q 1   = [ 14.1872 - 1.9381 - 1.9381 13.0104 ] , Q 2   = [ 3.1029 1.2838 1.2838 2.0181 ] , R 1 = [ 8 .3691 - 1 .3053 - 1 .3053 7 .0523 ] , R 1 = [ 5 .2020 2 .2730 2 .2730 1 .0238 ] , ε = 1.8043 , K 1 = [ - 0.4233   - 0.5031 ] , K 2 = [ - 0.5961   - 0.7049 ] , K 1 = [ - 0.4593   - 0.3874 ] . E43232623

Figures 24 illustrate the simulation results of applying the non-fragile fuzzy controller to the system (25) with x e = [ 3 .6626 2 .5443 ] T and u e = 77 .7272 under initial condition ϕ(t) = [1.2 −1.8]T, t ∈ [−2 0]. It can be seen that with the fuzzy control law, the closed-loop system is asymptotically stable and an upper bound of the guaranteed cost is J0 = 292.0399. The simulation results show that the fuzzy non-fragile guaranteed controller proposed in this paper is effective.

Figure 2.

State responses of x1(t).

Figure 3.

State responses of x2(t).

Figure 4.

Control trajectory of system.

Advertisement

6. Conclusions

In this paper, the problem of non-fragile guaranteed cost control for a class of fuzzy time-varying delay systems with local bilinear models has been explored. By utilizing the Lyapunov stability theory and LMI technique, sufficient conditions for the delay-dependent asymptotically stability of the closed-loop T-S fuzzy local bilinear system have been obtained. Moreover, the designed fuzzy controller has guaranteed the cost function-bound constraint. Finally, the effectiveness of the developed approach has been demonstrated by the simulation example. The robust non-fragile guaranteed cost control and robust non-fragile H-infinite control based on fuzzy bilinear model will be further investigated in the future work.

Advertisement

Acknowledgments

This work is supported by NSFC Nos. 60974139 and 61573013.

References

  1. 1. Pang CT, Lur YY. On the stability of Takagi-Sugeno fuzzy systems with time-varying uncertainties. IEEE Transactions on Fuzzy Systems. 2008;16:162–170
  2. 2. Zhou SS, Lam J, Zheng WX. Control design for fuzzy systems based on relaxed non-quadratic stability and H∞ performance conditions. IEEE Transactions on Fuzzy Systems. 2007;15:188–198
  3. 3. Zhou SS, Li T. Robust stabilization for delayed discrete-time fuzzy systems via basis-dependent Lyapunov-Krasovskii function. Fuzzy Sets and Systems. 2005;151:139–153
  4. 4. Gao HJ, Liu X, Lam J. Stability analysis and stabilization for discrete-time fuzzy systems with time-varying delay. IEEE Transactions on Systems, Man, Cybernetics, Part B. 2009;39:306–316
  5. 5. Wu HN, Li HX. New approach to delay-dependent stability analysis and stabilization for continuous-time fuzzy systems with time-varying delay. IEEE Transactions on Fuzzy Systems. 2007;15:482–493
  6. 6. Chen B, Liu XP. Delay-dependent robust H∞ control for T-S fuzzy systems with time delay. IEEE Transactions on Fuzzy Systems. 2005;13:238–249
  7. 7. Chen M, Feng G, Ma H, Chen G. Delay-dependent H∞ filter design for discrete-time fuzzy systems with time-varying delays. IEEE Transactions on Fuzzy Systems. 2009;17:604–616
  8. 8. Zhang J, Xia Y, Tao R. New results on H∞ filtering for fuzzy time-delay systems. IEEE Transactions on Fuzzy Systems. 2009;17:128–137
  9. 9. Keel LH, Bhattacharryya SP. Robust, fragile, or optimal?. IEEE Transactions on Automatic Control. 1997;42:1098–1105
  10. 10. Yang GH, Wang JL, Lin C. H∞ control for linear systems with additive controller gain variations. International Journal of Control. 2000;73:1500–1506
  11. 11. Yang GH, Wang JL. Non-fragile H∞ control for linear systems with multiplicative controller gain variations. Automatica. 2001;37:727–737
  12. 12. Zhang BY, Zhou SS, Li T. A new approach to robust and non-fragile H∞ control for uncertain fuzzy systems. Information Sciences. 2007;177:5118–5133
  13. 13. Yee JS, Yang GH, Wang JL. Non-fragile guaranteed cost control for discrete-time uncertain linear systems. International Journal of Systems Science. 2001;32:845–853
  14. 14. Mohler RR. Bilinear Control Processes. New York, NY: Academic; 1973
  15. 15. Elliott DL. Bilinear Systems in Encyclopedia of Electrical Engineering. New York, NY: Wiley; 1999
  16. 16. Li THS, Tsai SH. T-S fuzzy bilinear model and fuzzy controller design for a class of nonlinear systems. IEEE Transactions on Fuzzy Systems. 2007;15:494–505
  17. 17. Tsai SH, Li THS. Robust fuzzy control of a class of fuzzy bilinear systems with time-delay. Chaos, Solitons and Fractals. 2009;39:2028–2040
  18. 18. Li THS, Tsai SH, et al. Robust H∞ fuzzy control for a class of uncertain discrete fuzzy bilinear systems. IEEE Transactions on Systems, Man, Cybernetics, Part B. 2008;38:510–526
  19. 19. Dong J, Wang Y, Yang G. Control synthesis of continuous-time T-S fuzzy systems with local nonlinear models. IEEE Transactions on Systems, Man, Cybernetics, Part B. 2009;39:1245–1258
  20. 20. Zhang G, Li JM. Non-fragile guaranteed cost control of discrete-time fuzzy bilinear system. Journal of Systems Engineering and Electronics. 2010;21:629–634
  21. 21. Ho DWC, Niu Y. Robust fuzzy design for nonlinear uncertain stochastic systems via sliding-mode control. IEEE Transactions on Fuzzy Systems. 2007;15:350–358
  22. 22. Yang DD, Cai KY. Reliable guaranteed cost sampling control for nonlinear time-delay systems. Mathematics and Computers in Simulation. 2010;80:2005–2018
  23. 23. Chen WH, Guan ZH, Lu XM. Delay-dependent output feedback guaranteed cost control for uncertain time-delay systems. Automatica. 2004;44:1263–1268
  24. 24. Chen B, Lin C, Liu XP, Tong SC. Guaranteed cost control of T–S fuzzy systems with input delay. International Journal of Robust Nonlinear Control. 2008;18:1230–1256
  25. 25. Chen B, Liu XP, Tong SC, Lin C. Observer-based stabilization of T–S fuzzy systems with input delay. IEEE Transactions on Fuzzy Systems. 2008;16:652–663
  26. 26. Chen B, Liu X, Tong S, Lin C. guaranteed cost control of T-S fuzzy systems with state and input delay. Fuzzy Sets and Systems. 2007;158:2251–2267
  27. 27. Du BZ, Lam J, Shu Z. Stabilization for state/input delay systems via static and integral output feedback. Automatica. 2010;46:2000–2007
  28. 28. Kim JH. Delay-dependent robust and non-fragile guaranteed cost control for uncertain singular systems with time-varying state and input delays. International Journal of Control, Automation, and Systems. 2009;7:357–364
  29. 29. Li L, Liu XD. New approach on robust stability for uncertain T–S fuzzy systems with state and input delays. Chaos, Solitons and Fractals. 2009;40:2329–2339
  30. 30. Yu KW, Lien CH. Robust H-infinite control for uncertain T–S fuzzy systems with state and input delays. Chaos, Solitons and Fractals. 2008;37:150–156
  31. 31. Yue D, Lam J. Non-fragile guaranteed cost control for uncertain descriptor systems with time-varying state and input delays. Optimal Control Applications and Methods. 2005;26:85–105
  32. 32. Zhang G, Li JM. Non-fragile guaranteed cost control of discrete-time fuzzy bilinear system with time-delay. Journal of Dynamic Systems, Measurement and Control. 2014;136:044502–044504
  33. 33. Yue HY, Li JM. Output-feedback adaptive fuzzy control for a class of nonlinear systems with input delay and unknown control directions. Journal of the Franklin Institute. 2013;350:129–154
  34. 34. Yue HY, Li JM. Adaptive fuzzy tracking control for a class of perturbed nonlinear time-varying delays systems with unknown control direction. International Journal of Uncertainty, Fuzziness and Knowledge-based Systems. 2013;21:497–531
  35. 35. Wang RJ, Lin WW, Wang WJ. Stabilizability of linear quadratic state feedback for uncertain fuzzy time-delay systems. IEEE Transactions on Systems, Man, Cybernetics, Part B. 2004;34:1288–1292
  36. 36. Xia ZL, Li JM, Li JR. Delay-dependent fuzzy static output feedback control for discrete-time fuzzy stochastic systems with distributed time-varying delays. ISA Transaction. 2012;51:702–712
  37. 37. Xia ZL, Li JM. Switching fuzzy filtering for nonlinear stochastic delay systems using piecewise Lyapunov-Krasovskii function. International Journal of Fuzzy Systems. 2012;14:530–539
  38. 38. Li JR, Li JM, Xia ZL. Delay-dependent generalized H2 fuzzy static-output-feedback control for discrete T-S fuzzy bilinear stochastic systems with mixed delays. Journal of Intelligent and Fuzzy Systems. 2013;25:863–880
  39. 39. Xia ZL, Li JM, Li JR. Passivity-based resilient adaptive control for fuzzy stochastic delay systems with Markovian switching. Journal of the Franklin Institute-Engineering and Applied Mathematics. 2014;351:3818–3836
  40. 40. Li JM, Li YT. Robust stability and stabilization of fractional order systems based on uncertain T-S fuzzy model with the fractional order. Journal of Computational and Nonlinear Dynamics. 2013;8:041005
  41. 41. Li YT, Li JM. Stability analysis of fractional order systems based on T-S fuzzy model with the fractional order α: 0<α<1. Nonlinear Dynamics. 2014;78:2909–2919
  42. 42. Li YT, Li JM. Decentralized stabilization of fractional order T-S fuzzy interconnected systems with multiple time delays. Journal of Intelligent and Fuzzy Systems. 2016;30:319–331
  43. 43. Li JM, Yue HY. Adaptive fuzzy tracking control for stochastic nonlinear systems with unknown time-varying delays. Applied Mathematics and Computation. 2015;256:514–528
  44. 44. Yue HY, Yu SQ. Adaptive fuzzy tracking control for a class of stochastic nonlinearly parameterized systems with distributed input delay. Journal of the Franklin Institute-Engineering and Applied Mathematics. 2016;353:713–734
  45. 45. Duan RR, Li JM, Zhang YN, Yang Y, Chen GP. Stability analysis and H-inf control of discrete T-S fuzzy hyperbolic systems. International Journal of Applied Mathematics and Computer Science. 2016;26:133–145
  46. 46. Wang JX, Li JM. Stability analysis and feedback control of T-S fuzzy hyperbolic delay model for a class of nonlinear systems with time-varying delay. Iranian Journal of Fuzzy Systems. 2016;13:111–134
  47. 47. Li JR, Li JM, Xia ZL. Observer-based fuzzy control design for discrete time T-S fuzzy bilinear systems. International Journal of Uncertainty, Fuzziness and Knowledge-Based Systems. 2013;21:435–454

Written By

Junmin Li, Jinsha Li and Ruirui Duan

Submitted: 23 January 2017 Reviewed: 16 May 2017 Published: 30 August 2017