Abstract

The problem of stability analysis for a class of neutral systems with mixed time-varying neutral, discrete and distributed delays and nonlinear parameter perturbations is addressed. By introducing a novel Lyapunov-Krasovskii functional and combining the descriptor model transformation, the Leibniz-Newton formula, some free-weighting matrices, and a suitable change of variables, new sufficient conditions are established for the stability of the considered system, which are neutral-delay-dependent, discrete-delay-range-dependent, and distributed-delay-dependent. The conditions are presented in terms of linear matrix inequalities (LMIs) and can be efficiently solved using convex programming techniques. Two numerical examples are given to illustrate the efficiency of the proposed method.

1. Introduction

Delay (or memory) systems represent a class of infinite-dimensional systems largely used to describe propagation and transport phenomena or population dynamics [13]. Delay differential systems are assuming an increasingly important role in many disciplines like economic, mathematics, science, and engineering. For instance, in economic systems, delays appear in a natural way since decisions and effects are separated by some time interval. The presence of a delay in a system may be the result of some essential simplification of the corresponding process model. The problem of delay effects on the stability of systems including delays in the state, and/or input is a problem of recurring interest since the delay presence may induce complex behaviors (oscillation, instability, bad performances) for the schemes [1, 2]. Some improved methods pertaining to the problems of determining robust stability criteria and robust control design of uncertain time-delay systems have been reported; see, for example, [4, 5] and the references cited therein. When dealing with time-varying delays and the reduction of the level of design conservatism, one has to select appropriate Lyapunov-Krasovskii functional (LKF) with moderate number of terms [6].

Neutral delay systems constitute a more general class than those of the retarded type. Stability of these systems proves to be a more complex issue because the system involves the derivative of the delayed state. Especially in the past few decades, increased attention has been devoted to the problem of robust delay-independent stability or delay-dependent stability and stabilization via different approaches (e.g., model transformation techniques [2, 79], the improved bounding techniques [10, 11], and the properly chosen Lyapunov-Krasovskii functionals [12, 13]) for a number of different neutral systems with delayed state and/or input, parameter uncertainties, and nonlinear perturbations (see, e.g., [1425] and the references therein).

Among the existing results on neutral delay systems, the linear matrix inequality (LMI) approach is an efficient method to solve many problems such as stability analysis, stabilization [9, 15, 26, 27], ??8 control problems [2830], filter designs [31, 32], and guaranteed-cost (observer-based) control [3339]. Besides, for neutral systems with mixed neutral and discrete delays, most of the aforementioned methods can only provide neutral-delay-independent and discrete-delay-dependent results. Furthermore, the subject of the robust stability and feedback stabilization of continuous- and discrete-time systems (within the framework LMI) under additive perturbations which are nonlinear functions in time and state of the systems are investigated in [40, 41], respectively.

In the recent literature on neutral systems, He et al. in [42] proposed a new approach to analyze the stability of the systems with mixed delays by incorporating some free-weighting matrices, and the less conservative criteria, which were both discrete-delay-dependent and neutral-delay-dependent, were obtained without considering the model transformations. However, some of the free matrices did not serve to reduce the conservatism of the results that were obtained. Moreover, in [9, 20], the authors studied the problem of the robust stability of neutral systems with nonlinear parameter perturbations and mixed time-varying neutral and discrete delays and presented neutral-delay-independent stability criteria, that cannot be directly applied to the systems with different time-varying neutral, discrete, and distributed delays. Furthermore, from the published results, it appears that general results pertaining to neutral systems with mixed time-varying neutral, discrete, and distributed delays and nonlinear parameter perturbations are few and restricted; see [9, 10, 18, 20, 42] where most of the efforts were virtually neutral-delay-range-independent or were not centered on distributed delays.

In this paper, we develop new stability criteria for the stability analysis of the neutral systems with nonlinear parameter perturbations based on a descriptor model transformation. The dynamical system under consideration consists of time-varying neutral, discrete, and distributed delays without any restriction on upper bounds of derivatives of time-varying delays. By introducing a novel Lyapunov-Krasovskii functional and combining the descriptor model transformation, the Leibniz-Newton formula, some free-weighting matrices, and a suitable change of variables, new sufficient conditions are established for the stability of the considered system, which are neutral-delay-dependent, discrete-delay-range-dependent, and distributed-delay-dependent. The conditions are presented in terms of LMIs and can be easily solved by existing convex optimization techniques. Two numerical examples are given to demonstrate the less conservatism of the proposed results over some existence results in the literature.

Notations. The superscript ???? stands for matrix transposition; R?? denotes the n-dimensional Euclidean space; R??×?? is the set of all real m by n matrices. ?·? refers to the Euclidean vector norm or the induced matrix 2-norm. col{?} and diag{?} represent, respectively, a column vector and a block diagonal matrix, and the operator sym(??) represents ??+????. ??min(??) and ??max(??) denote, respectively, the smallest and largest eigenvalue of the square matrix ??. The notation ??>0 means that ?? is real symmetric and positive definite; the symbol * denotes the elements below the main diagonal of a symmetric block matrix.

2. Problem Description

Consider a class of linear neutral systems with different time-varying neutral, discrete, and distributed delays and nonlinear parameter perturbations represented by ??(??)-????(??-??(??))=????(??)+????(??-h(??))+??1??1(??,??(??))+??2??2(??,??(??-h(??)))+??3?????-??(??)??3(??,??(??))????+??4??4???,??(??-??(??))?,??(??)=??(??),???[-??,0],(2.1) where ??:=max{h2,??1,??1}, and ??(??)?R?? is the state vector. The time-varying vector valued initial function ??(??) is a continuously differentiable functional, and the time-varying delays h(??), ??(??), and ??(??) are functions satisfying, respectively, 0<h1=h(??)=h2,||?h(??)||=h3<8,(2.2a)0<??(??)=??1,||??(??)||=??2<8,(2.2b)0<??(??)=??1,||??(??)||=??2<8.(2.2c) The time-varying vector-valued functions ????:R+×R???R???? (??=1,,4) are continuous and satisfy ????(??,0)=0, and the Lipschitz conditions, that is, ?????(??,??0)-????(??,??0)?=?????(??0-??0)? for all ?? and for all ??0,??0?R?? such that ???? are some known matrices.

Remark 2.1. In this case, h(??) is called an interval-like or range-like time-varying delay [14]. It is also noted that this kind of time-delay describes the real situation in many practical engineering systems. For example, in the field of networked control systems, the network transmission induced delays (either from the sensor to the controller or from the controller to the plant) can be assumed to satisfy (2.2a) without loss of generality [43, 44].

Throughout the paper, the following assumptions are needed to enable the application of Lyapunov’s method for the stability of neutral systems [1]:

(A1)let the difference operator ??:??([-??,0],R??)?R?? given by ??????=??(??)-????(??-??(??)) be delay-independently stable with respect to all delays. A sufficient condition for (A1) is that(A2)all the eigenvalues of the matrix ?? are inside the unit circle.

Before ending this section, we recall the following lemmas, which will be used in the proof of our main results.

Lemma 2.2 (see [9]). For any arbitrary column vectors ??(??),??(??)?R??, any matrix ???R??×??, and positive-definite matrix ???R??×??, the following inequality holds: -2?????-??(??)??(??)????(??)????=?????-??(??)???(??)??(??)??????????*(????+??)????-1(????+??)????(??)??(??)?????.(2.3)

Lemma 2.3 (see [45]). Given matrices ??=????,??, ??, and ?? of appropriate dimensions with ??????=??, then the following matrix inequality holds: ??+sym(??????)<0,(2.4) for all ?? if and only if there exists a scalar ??>0 such that ??+????????+??-1??????<0.(2.5)

3. Main Results

In this section, new delay-range-dependent sufficient conditions for the asymptotic stability of the neutral system (2.1) are presented. By utilizing the Leibniz-Newton formula, the following two zero equations hold: ??1??(??)-??1??(??-h(??))-??1?????-h(??)??(??)????=0,(3.1a)??2??(??)-??2??(??-??(??))-??2?????-??(??)??(??)????=0,(3.1b) then, we can represent the system (2.1) as ??(??)-????(??-??(??))=?????(??)+?????(??-h(??))+??1??1(??,??(??))+??2??2(??,??(??-h(??)))-??2??(??-??(??))-??1?????-h(??)??(??)????-??2?????-??(??)??(??)????+??3?????-??(??)??3(??,??(??))????+??4??4(??,??(??-??(??))),(3.2) with ???:=??+??1+??2 and ???:=??-??1 where the matrices ??1,??2?R??×?? will be chosen in the following theorem.

Theorem 3.1. Under (A1), for given scalars ??,h1,h2,??1,??1>0,h3,??2,??2, the neutral system (2.1) is asymptotically stable, if there exist some scalars ??,??1,??2, matrices ??2,{????}20??=1,??1,??2, and positive-definite matrices ??1,{????}4??=1,{????}4??=1,???1,???2, such that the following LMI is feasible: ??=??????????????????????11??12??13??14??15??16??17??18??19??1,10??1,11*??22??230-????????170??27-????????19??29??2,100**??33000??370??3900***??4400-????150-????16??4,100****??55??56-????????18-??17??3-????????2000*****??66??67??68??6900******??77-??18??3-??18??4??7,100*******??88-????3????2000********-sym???20??4???9,100*********??10,100**********??11,11????????????????????<0(3.3) with ??11=sym??????????2??+?1+??1???1+?1+??2???2??1-????2???????2??+?1+??1???1+?1+??2???2?-??????2??????+O,??12=???????2??-?1+??1???1+????2-??1-??5+??9??????2??-???1+??1???1???,??13=???5??900?,??14=???-?1+??2???2+????14-??13-???1+??2???2???,??15=???????2??-????????17??????2??+????17???,??16=???????2??1????2??2??????2??1??????2??2???,??17=???????3+????15-????????18????18???,??18=???????2??3-????????19??????2??3+????19???,??19=???????2??4+????4+????16-????????20??????2??4+????20???,??1,10=?h12??1h12??5h12??9??1??130000?,??1,11=?h2???1+1????1??1???1+1????2h2?0??????1????1?????1?0??????2????2????,??22=????2??2-?1-h3???3-sym???2+??6-??10?,??23=???6??10?,??27=-????3-????7+????11-????????18,??29=-????4-????8+????12-????????20,??2,10=?h12??2h12??6h12??100?,??33=diag?-??1,-??2?,??37=???????7-????11???,??39=???????8-????12???,??44=-?1-??2???1-sym???14?,??4,10=?000??1??14?,??55=-?1-??2???2-sym???17???+????4??4,??56=?-??17??1-??17??2?,??66=diag{-??,-??},??67=???-????1????18-????2????18???,??68=???19??1??19??2???,??69=???-????1????20-????2????20???,??77=-??+??21??4,??88=-?1-??2???4-sym???19??3?,??4,10=?000??1??14?,??7,10=?h12??3h12??7h12??11??1??15?,??9,10=?h12??4h12??8h12??12??1??16?,??10,10=diag?-h12??4,-h12??5,-h12??4,-??1??3?,??1,11=diag?-h2???1,-??1???2,-h2???1,-??1???2?,(3.4) where (3.4)O=diag{??1+?3??=1????+????1??1+????3??3+sym(??1+??13),??2+??1??3+h12??4+h2??5}.

Proof Firstly, we represent () in an equivalent descriptor model form as ??(??)=??(??),0=-??(??)+?????(??)+????(??-??(??))+?????(??-h(??))-??2??(??-??(??))+??1??1(??,??(??))+??2??2(??,??(??-h(??)))-??1?????-h(??)??(??)????-??2?????-??(??)??(??)????+??3?????-??(??)??3(??,??(??))????+??4??4(??,??(??-??(??))).(3.4) Define the Lyapunov-Krasovskii functional ??(??)=6???=1????(??),(3.6) where ??1(??)=??(??)????1??(??):=??(??)????????(??),??2(??)=?????-??(??)???(??)????1??(??)+??(??)????2??(??)?????+?????-h(??)??(??)????3??(??)????+2???=1?????-h????(??)????????(??)????,??3(??)=?-h1-h2?????+????(??)????4??(??)????????+?0-h2?????+????(??)????5??(??)????????,??4(??)=?0-??1?????+????(??)????3??(??)????????,??5(??)=?0-h2?????+????(??)???0??1?????1?0??1???(??)????????+?0-??1?????+????(??)???0??2?????2?0??2???(??)????????,??6(??)=?????-??(??)????????3(??,??(??))?????????4????????3(??,??(??))???????????+???10?????-??(??-??+??)??3(??,??(??))????4??3(??,??(??))????????,(3.7) with ??(??):=col{??(??),??(??)}, ??=diag{??,0}, and ??=[??10??2??3], where ??1=????1>0.
On the other hand, noting that ??(??(??),??)=??min(??1)???(0)?2. According to [34], using the Cauchy-Schwarz inequality and after some manipulations, we obtain ??(??(??),??)=??(??(0),0)=??????(0)?2+?0-???????(??)??2?????,(3.8) where ??:=max(??1,??2) with ??1:=??max???1?+2??1??max???1?+2h1??max???1?+2h2??max???2?+2h2??max???3?+3??21??max?????3??4??3?,??2:=2??21??max???1?+??max???2?+2h22??max???3?+2h21??max???1?+2h22??max???2?+h2??max???5?+?h1+h2???max???4?+??1??max???3?+h2??max??0??1?????1?0??1??+??1??max??0??2?????2?0??2??+113??31??max?????3??4??3?.(3.9) Differentiating ??1(??) along the system trajectory becomes ???1(??)=2??(??)????1???(??)=2??(??)?????????(??)0?=2??(??)???????????(??)+?0?????(??-??(??))+?0??????(??-h(??))-?0??2???(??-??(??))+?0??1???1(??,??(??))+?0??2???2(??,??(??-h(??)))+??3?????-??(??)??3(??,??(??))????+??4??4(??,??(??-??(??)))?+??1(??)+??2(??),(3.10) where ??:=?0?????-???,??1(??)=-2?????-h(??)??(??)???????0??1???(??)????,??2(??)=-2?????-??(??)??(??)???????0??2???(??)????.(3.11)
Using Lemma 2.2 for ??(??)=col{0,????}??(??) and ??=??col{??(??),??(??)}, we obtain ??1(??)=h2??(??)???????????1??1+???????-11?????1??1+???????(??)+2??(??)??????????1??1?0??1?(??(??)-??(??-h(??)))+?????-h2??(??)???0??1?????1?0??1???(??)????,??2(??)=??1??(??)???????????2??2+???????-12?????2??2+???????(??)+2??(??)??????????2??2?0??2?(??(??)-??(??-??(??)))+?????-??1??(??)???0??2?????2?0??2???(??)????.(3.12) Differentiating the second Lyapunov term in (3.6) gives ???2(??)=??(??)?????1+3???=1???????(??)+??(??)????2??(??)-?1-?h(??)?????(??-h(??))??3??(??-h(??))-?1-??(??)?????(??-??(??))??1??(??-??(??))-?1-??(??)?????(??-??(??))??2??(??-??(??))-2???=1?????-h??????????????-h???=??(??)?????1+3???=1???????(??)+??(??)????2??(??)-?1-h3?????(??-h(??))??3??(??-h(??))-?1-??2?????(??-??(??))??1??(??-??(??))-?1-??2?????(??-??(??))??2??(??-??(??))-2???=1?????-h??????????????-h???,(3.13) and the time derivative of the third term of ??(??) in (3.6) is ???3(??)=??(??)???h12??4+h2??5???(??)-???-h1??-h2??(??)????4??(??)????-?????-h2??(??)????5??(??)????=??(??)???h12??4+h2??5???(??)-???-h(??)??-h2??(??)????4??(??)????-???-h1??-h(??)??(??)?????4+??5???(??)????,(3.14) and, similarly, ???4(??)=??1??(??)????3??(??)-?????-??1??(??)????3??(??)????=??1??(??)????3??(??)-?????-??(??)??(??)????3??(??)????,(3.15) and also the time derivative of the fifth and sixth terms of ??(??) in (3.6) are, respectively, ???5(??)=??(??)???h2?0??1?????1?0??1?+??1?0??2?????2?0??2????(??)-?????-h2??(??)???0??1?????1?0??1???(??)????-?????-??1??(??)???0??2?????2?0??2???(??)????,(3.16)???6(??)=-?1-??(??)???????-??(??)??3(??,??(??))?????????4??????-??(??)??3(??,??(??))?????+2?????-??(??)??3(??,??(??))????4????????3(??,??(??))?????????+???10????3(??,??(??))????4??3(??,??(??))????-???10?????-????3(??,??(??))????4??3(??,??(??))????????=?????-??(??)(??-??+??(??))???3(??,??(??))????4??3(??,??(??))+??3(??,??(??))????4??3(??,??(??))?????+???10????3(??,??(??))????4??3(??,??(??))????-?1-??2???????-??(??)??3(??,??(??))?????????4??????-??(??)??3(??,??(??))?????-?????-??1???-??+??1???3(??,??(??))????4??3(??,??(??))????=??21??3(??,??(??))????4??3(??,??(??))-?1-??2???????-??(??)??3(??,??(??))?????????4??????-??(??)??3(??,??(??))?????.(3.17) For nonlinear functions ????(·), we have 0=-??1(??,??(??))????1(??,??(??))+??(??)??????1??1??(??),0=-??2(??,??(??-h(??)))????2(??,??(??-h(??)))+??(??-h(??))??????2??2??(??-h(??)),0=-??3(??,??(??))????3(??,??(??))+??(??)??????3??3??(??),0=-??4(??,??(??-??(??)))????4(??,??(??-??(??)))+??(??-??(??))??????4??4??(??-??(??)).(3.18) Moreover, from the Leibniz-Newton formula and the system (2.1), the following equations hold for any matrices {????}10??=1 with appropriate dimensions: 2????(??)??1???(??)-??(??-h(??))-?????-h(??)??(??)?????=0,2????(??)??2??????-h1?-??(??-h(??))-???-h1??-h(??)??(??)?????=0,2????(??)??3???(??-h(??))-?????-h2?-???-h(??)??-h2??(??)?????=0,2????(??)??4???(??)-??(??-??(??))-?????-??(??)??(??)?????=0,2????(??)??5???(??)-????(??-??(??))-????(??)-????(??-h(??))-??1??1(??,??(??))-??2??2(??,??(??-h(??)))-??3?????-??(??)??3(??,??(??))????-??4??4(??,??(??-??(??)))?=0,(3.19) where ??1:=?????1,0,????2,0,,0???????6elements,????3,0,????4???,??2:=?????5,0,????6,0,,0???????6elements,????7,0,????8?????3:=?????9,0,????10,0,,0???????6elements,????11,0,????12???,??4:=?????13,0,,0???????4elements,????14,0,,0???????3elements,????15,0,????16???,??5:=?0,,0???????6elements,????17,0,0,????18,????19,????20???,??(??):=col???(??),??(??),??(??-h(??)),?????-h1?,?????-h2?,??(??-??(??)),??(??-??(??)),??1(??,??(??)),??2(??,??(??-h(??))),??3(??,??(??)),?????-??(??)??3(??,??(??))????,??4(??,??(??-??(??)))?.(3.20) Using the obtained derivative terms (3.10)–(3.17) and adding the right- and the left-hand sides of (3.18) and (3.19) into ???(??), the following result is obtained: ???(??)=6???=1?????(??)=??(??)??S??(??)-???-h1??-h(??)?????(??)??1+????(??)??4???-14?????(??)??1+????(??)??4???????-???-h1??-h(??)?????(??)??2+????(??)??5???-15?????(??)??2+????(??)??5???????-???-h1??-h(??)?????(??)??3+????(??)??4???-14?????(??)??3+????(??)??5???????-?????-??(??)?????(??)??4+????(??)??3???-13?????(??)??4+????(??)??3???????,(3.21) where S:=???+h12??1??-14????1+h12??2??-15????2+h12??3??-14????3+??1??4??-13????4, and the matrix ??? is given by ???=????????????????????11???12??13???14???15???16??17???18???19*??22??230-????????170??27-????????19??29**??33000??370??39***??4400-????150-????16****??55??56-????????18-??17??3-????????20*****??66??67??68??69******??77-??18??3-??18??4*******??88-????3????20********-sym???20??4??????????????????(3.22) with ???11=0x000a0sym0x000a0????????+?????1??1?0??1?+????2??2?0??2??[??0]??+?????h2?????1??1+???????-11?????1??1+???+??1?????2??2+???????-12?????2??2+??????+?0????h2?0??1?????1?0??1?+??1?0??2?????2?0??2???0?????+O,???12=?????0????-????????1??1?0??1?+?????2-??1-??5+??90?,???14=-?????0??2?-????????2??2?0??2?+?????14-??130?,???15=?????0???+???-????????17????17???,???16=?????00??1??2?,???18=?????0??3?+???-????????19????19???,???19=?????0??4?+???????4+????16-????????20????20???.(3.23) Now, if S<0 holds, then ???(??)<0 which means that the neutral system (2.1) is asymptotically stable. By applying the Schur complement, the matrix inequity S<0 results in ????????????????????????11???12??13???14???15???16??17???18???19??1,10???1,11*??22??230-????????170??27-????????19??29??2,100**??33000??370??3900***??4400-????150-????16??4,100****??55??56-????????18-??17??3-????????2000*****??66??67??68???6900******??77-??18??3-??18??4??7,100*******??88-????3????2000********-sym???20??4???9,100*********??10,100**********???11,11?????????????????????<0(3.24) with ???11=sym????????+?????1??1?0??1?+????2??2?0??2?????0???+O,???1,11=?h2?????????1??1+???????1??1?????????2??2+???????2h2?0????0??1?????1?0????0??2????,???11,11=diag?-h2??1,-??1??2,-h2??1,-??1??2?,(3.25) where ????=??-1??(??=1,2).
Following [34, 35], we choose ??3=????2, ?????, where ?? is a tuning scalar parameter (which may be restrictive). Let ??=diag???????,,?????????16elements,????,,???????????????????4elements?????.(3.26)Premultiplying ?? and postmultiplying ???? to the matrix inequality (3.24) and considering ????:=????2????, ?????:=??????????, and ????????=?????? (??=1,2) to eliminate the nonlinearities in the matrix inequality, we obtain the LMI (3.3). Moreover, from (2.1) and the fact that ??(??) is square integrable on [0,8), it follows that ???????????2[0,8). Under (A1), the later implies that ??(??-??(??))?????2[0,8). Therefore, by [1, Theorem 1.6], we conclude that the neutral system (2.1) is asymptotically stable.

Remark 3.2. The results given in Theorem 3.1 are derived for system (2.1) with time-varying delays h(??), ??(??), and ??(??) satisfying (2.2a), (2.2b), and (2.2c), where the derivatives of the time-varying delays are available. However, in many situations, the information on the derivative of time-varying delays is unknown a prior. In such circumstances, the corresponding delay-rate-independent stability analysis results for time-delays only satisfying 0<h1=h(??)=h2<8,0<??(??)=??1<8,0<??(??)=??1<8,(3.27) can be easily obtained by setting ??1=??2=??4=??3=0 in Theorem 3.1.

Remark 3.3. The reduced conservatism of Theorem 3.1 benefits from the construction of the new Lyapunov-Krasovskii functional in (3.6), utilizing Leibniz-Newton formula, using a free-weighting matrix technique, and no bounding technique is needed to estimate the inner product of the involved crossing terms (see, e.g., [12, 20]). It can be easily seen that results of this paper are quite different from most existing results in the recent literature in the following perspectives. (a) Theoretically stability analysis of neutral systems with different time-varying neutral, discrete, and distributed delays is much more complicated, especially, for the case where the delays are time-varying and different. (b) In this paper, the derived sufficient conditions are convex, neutral-delay-dependent, discrete-delay-range-dependent, and distributed-delay-dependent, which make the treatment in the present paper more general with less conservative in compare to most existing results in the literature which are independent of the neutral or distributed delays; see for instance [21, 22, 38].

4. Uncertainty Characterization

In this section, we will discuss the uncertainty characterization for the linear neutral system (2.1) with different time-varying neutral, discrete, and distributed delays and nonlinear parameter perturbations.

4.1. Polytopic Uncertainty

The first class of uncertainty frequently encountered in practice is the polytopic uncertainty [2]. In this case, the matrices of the system (2.1) are not exactly known, except that they are within a compact set O denoting O=?????????1??2??3?.(4.1) We assume that O=?????=1????O??(4.2) for some scalars ???? satisfying 0=????=1,?????=1????=1,(4.3) where the ?? vertices of the polytope are described by O??=???(??)??(??)??(??)??(??)1??(??)2??(??)3?.(4.4) In order to take into account the polytopic uncertainty in the system (2.1), we derive the following result from applying the same transformation that was used in deriving Theorem 3.1.

Theorem 4.1. Under (A1), for given scalars ??,h1,h2,??1,??1>0,h3,??2,??2, if the uncertainty set O is polytopic with vertices O??, ??=1,2,,??, then the system described by (2.1)), (2.2a), (2.2b), (2.2c), and (4.2)–(4.4) is asymptotically stable if there exist some scalars ??,??1,??2, matrices ??2,{????}20??=1,??1,??2, and positive-definite matrices ??1,{????}4??=1,{????}4??=1,???1,???2 such that LMI (3.3) is satisfied for all ?????????1??2??3?=???(??)??(??)??(??)??(??)1??(??)2??(??)3?,??=1,2,,??.(4.5)

Proof. It follows directly from the proof of Theorem 3.1 and using properties of (4.2)–(4.4).

4.2. Norm-Bounded Uncertainty

There are also other uncertainties that cannot be reasonably modeled by a polytopic uncertainty set with a number of vertices. In such a case, it is assumed that the deviation of the system parameters of an uncertain system from their nominal values is norm bounded [2]. In our case, consider the time-varying structured uncertain neutral system ??(??)-(??+???(??))??(??-??(??))=(??+???(??))??(??)+(??+???(??))??(??-h(??))+???1+???1(??)???1(??,??(??))+???2+???2(??)???2(??,??(??-h(??)))+???3+???3(??)??????-??(??)??3(??,??(??))????+???4+???4(??)???4???,??(??-??(??))?,??(??)=??(??),???[-??,0],(4.6) where the time-varying structured uncertainties are said to be admissible if the following form holds: ????(??)???(??)???(??)???1(??)???2(??)???3(??)???4(??)?=??1?(??)?????????????????1????2????3????4?,(4.7) where ????,????,????,????1,????2,????3,????4 are constant matrices with appropriate dimensions, ?(??) is an unknown, real, and possibly time-varying matrix with Lebesgue measurable elements, and its Euclidean norm satisfies ??(??)?=1,???.(4.8) In this section, we modify (A1)-(A2) in order to enable the application of Lyapunov’s method for the stability of the time-varying structured uncertain neutral system (4.6) as follows:

(A'1)let the difference operator ??:??([-??,0],R??)?R?? given by ??????=??(??)-(??+???(??))??(??-??(??)) be delay-independently stable with respect to all delays. A sufficient condition for (A'1) is that (A'2)all the eigenvalues of the matrix ??+???(??) are inside the unit circle, that is, ???+???(??)?<1.

Theorem 4.2. Under (A'1), for given scalars ??,h1,h2,??1,??1>0,h3,??2,??2, the neutral system (4.6) with admissible uncertainties (4.7)) and (4.8) is robustly asymptotically stable if there exist some scalars {????}4??=1>0, ??,??1,??2, matrices ??2,{????}20??=1,??1,??2, and positive-definite matrices ??1,{????}4??=1,{????}4??=1,???1,???2, such that the following LMI is feasible: ???????1?2*-?30**-?3????<0,(4.9) where ?1=[G??1G??2G??3G??4], ?2=[??1G????1??2G????2??3G????3??4G????4], and ?3=diag{??1??,??2??,??3??,??4??} with G??1=?????1??2??????1??20?0?????????18elements???,G??1=?????0????0?0?????????3elements????????1????20????3????40?0?????????8elements?,G??2=?0?0?????????6elements????1????170?0?????????13elements???,G??2=?-????0-????0?0?????????3elements-????0?0?????????13elements?,G??3=?0?0?????????7elements-????1????170?0?????????12elements???,G??3=?0?0?????????7elements????1????20????3????40?0?????????8elements?,G??4=?0?0?????????9elements????1????190?0?????????10elements???,G??4=?00-????0?0?????????4elements-????1-????20-????3-????40?0?????????8element???.(4.2)

Proof. If the matrices ??,??,??,??1,??2,??3,??4 in (3.3) are replaced with ??+??1?(??)????, ??+??1?(??)????, ??+??1?(??)????, ??1+??1?(??)????1, ??2+??1?(??)????2, ??3+??1?(??)????3, and ??4+??1?(??)????4, respectively, then (3.3) with the admissible uncertainties (4.7) is equivalent to the following condition: ??+4???=1sym?G???????(??)G?????<0.(4.11) By Lemma 2.3, a necessary and sufficient condition for (4.11) is that there exist some {????}4??=1>0 such that ??+4???=1???-1??G??????G????+????G??????G?????<0.(4.12) Applying Schur complements, we find that (4.12) is equivalent to (4.9).

Remark 4.3. It is noted that our approach is different from that in the reference [20] in several perspectives. (a) The system structure in [20] considers a system with different time-varying neutral and discrete delays and in compare to our case do not center on time-varying distributed delays, that is, the results in [20] cannot be directly applied to the systems with different time-varying neutral, discrete, and distributed delays. (b) Their system only considers the case that the range of the time-varying delay h(??) is from 0 to an upper bound in compare to our case that the time-varying discrete delay (h(??)) lies in a range, in which the lower bound is not 0. (c) The derived neutral-delay-range-independent conditions and using the inequality bounding technique [11, Lemma 1] employed for some cross terms encountered in their analysis conditions may produce conservative results in comparison with the present paper.

5. Numerical Results

In this section, two examples are provided to illustrate the effectiveness of the results obtained in the previous sections.

Example 5.1. Consider the neutral system (2.1) with the following matrices adopted from [9]: ??=?-1.20.1-0.1-1?,??=?-0.60.7-1-0.8?,??1=??2=??,??3=0,(5.1) and the nonlinear parameter perturbations ???1(??,??(??))?=??1???(??)? and ???2(??,??(??-h(??)))?=??2???(??-h(??))?, where ??1,??2=0.

Case 1. Assume that ??=0, h1=0, and ??4=0. Applying the criteria in [10, 18, 21, 22] and Theorem 3.1 in this paper, the maximum value of h2 for stability of system under consideration is listed in Table 1. Furthermore, in the case of h3=0.5 and ??1=??2=0.1 by the criteria in [20], the nominal system is asymptotically stable for any h(??) satisfying h2=1.0097. It is easy to see that the stability criterion in this paper gives a much less conservative result than that in [10, 2022], and also the maximum value of h2 decreases as ???? increases. Moreover, unlike the results obtained in [10, 18, 21, 22] our approach can also consider the case h3=1 and handle fast time-varying delays completely.

Case 2. Assume that ??=diag{0.1,0.1}, ??4=??, h1=0,h3=??2=0.5, and ???4(??,??(??-??(??)))?=??3???(??-??(??))?. We also calculate the maximum delay bounds h2 that guarantee the asymptotic stability of the system under consideration by Theorem 3.1 in this paper and developed methods in [9, 20] for different values of the parameter ??3, as listed in Table 2. It can be seen from the table that h2 decreases as ??3 increases, and the stability condition in Theorem 3.1 of this paper is less conservative than that in [9, 20].

Example 5.2 (see [46]). Consider the system (4.6) with the following matrices: ??=?-200-3?,??1=?0.5-0.1-0.2-0.3?,??2=?0.2-0.40.10.2?,??3=?0.40.30.10.2?,????=?000.20.2?,????1=?000.10.1?,????2=?000.30.3?,????3=?000.20.2?,??=??=??4=????=????=????4=0,??1=??,??1=??2=??3=??.(5.2) Recently, Kwon and Park in [46] derived the stability bound of h(??)=??(??) with h3=??2=1 as 0<h(??)=1.97 with the parameters above for the system (4.6). However, by applying Theorem 4.2 to the system under consideration, one can see that our criterion is feasible for 0<h(??)=2.0105. And for the condition h3=1, for fixed h2 or ??2, the upper bounds of delays h2 and ??2 are shown in Tables 3 and 4, respectively. From Tables 3 and 4, it can be seen that Theorem 4.2 provides a condition for guaranteeing stability with respect to given delay bounds h(??) and ??(??).

6. Conclusion

The problem of stability analysis has been presented in this paper for a class of neutral systems with different time-varying neutral, discrete, and distributed delays and nonlinear parameter perturbations using an appropriate Lyapunov-Krasovskii functional. By combining the descriptor model transformation, the Leibniz-Newton formula, some free-weighting matrices, and a suitable change of variables, new feasibility conditions, which are neutral-delay-dependent, discrete-delay-range-dependent, and distributed-delay-dependent, have been developed to ensure that the considered system is asymptotically stable. The conditions were presented in terms of linear matrix inequalities (LMIs) and solved by existing convex optimization techniques. Two numerical examples were given to demonstrate the less conservatism of the proposed results over some existence results in the literature.

Acknowledgments

This work has been partially funded by the European Union (European Regional Development Fund) and the Ministry of Science and Innovation (Spain) through the coordinated research project DPI2008-06699-C02-01. H. R. Karimi is grateful to the grant of Juan de la Cierva program of the Ministry of Science and Innovation (Spain), and M. Zapateiro is grateful to the FI grant of the Department for Innovation, University and Enterprise of the Government of Catalonia (Spain). The authors would also like to thank the associate editor and anonymous reviewers for their constructive comments and suggestions to improve the quality of the paper.