## Networked control systems in the presence of scheduling protocols and communication delays ∗

Kun Liu , Emilia Fridman, Laurentiu Hetel † ‡ §

August 4, 2014

## Abstract

This paper develops the time-delay approach to Networked Control Systems (NCSs) in the presence of variable transmission delays, sampling intervals and communication constraints. The system sensor nodes are supposed to be distributed over a network. Due to communication constraints only one node output is transmitted through the communication channel at once. The scheduling of sensor information towards the controller is ruled by a weighted Try-Once-Discard (TOD) or by Round-Robin (RR) protocols. Differently from the existing results on NCSs in the presence of scheduling protocols (in the frameworks of hybrid and discrete-time systems), we allow the communication delays to be greater than the sampling intervals. A novel hybrid system model for the closed-loop system is presented that contains time-varying delays in the continuous dynamics and in the reset conditions . A new Lyapunov-Krasovskii method, which is based on discontinuous in time Lyapunov functionals is introduced for the stability analysis of the delayed hybrid systems. Polytopic type uncertainties in the system model can be easily included in the analysis. The efficiency of the time-delay approach is illustrated on the examples of uncertain cart-pendulum and of batch reactor.

Key words networked control systems, time-delay approach, scheduling protocols, hybrid systems, Lyapunov-Krasovskii method.

AMS subject classification. 93D15, 93D05.

## 1 Introduction

Networked Control Systems (NCSs) are systems with spatially distributed sensors, actuators and controller nodes which exchange data over a communication

∗ This work was partially supported by Israel Science Foundation (grant No 754/10), the Knut and Alice Wallenberg Foundation and the Swedish Research Council.

† ACCESS Linnaeus Centre and School of Electrical Engineering, KTH Royal Institute of Technology, SE-100 44 Stockholm, Sweden, kunliu@kth.se

‡ School of Electrical Engineering, Tel Aviv University, Tel Aviv, 69978 Israel, emilia@eng.tau.ac.il

§ University Lille Nord de France, LAGIS, FRE CNRS 3303, Ecole Centrale de Lille, Cite Scientifique, BP 48, 59651 Villeneuve d'Ascq cedex, France, laurentiu.hetel@ec-lille.fr

data channel [1]. In many NCSs, only one node is allowed to use the communication channel at once. The communication along the data channel is orchestrated by a scheduling rule called protocol. The introduction of communication network media offers several practical advantages: reduced costs, ease of installation and maintenance and increased flexibility.

Three main approaches have been used to model the sampled-data control and later to the NCSs: a discrete-time [6, 10], an impulsive/hybrid system [18, 20] and a time-delay [7, 9, 11, 12] approaches. The hybrid system approach, which was inspired by [23], has been applied to nonlinear NCSs under Try-OnceDiscard (TOD) and Round-Robin (RR) protocols in [15, 20]. In the framework of discrete-time approach, network-based stabilization of linear time-invariant systems with TOD/RR protocols and communication delays has been considered in [6]. Variable sampling intervals and/or small communication delays (that are smaller than the sampling intervals ) have been considered in the above works.

Note that in the absence of scheduling protocols, all the three approaches are applicable to non-small communication delays (see e.g., [3, 19]). The time-delay approach that was recently suggested in [16] allowed, for the first time, to treat NCSs under RR protocol in the presence of non-small communication delays. In [16] the closed-loop system was presented as a switched system with multiple ordered delays.

In the present paper, we consider linear (probably, uncertain) NCS with additive essentially bounded disturbances in the presence of scheduling protocols, variable sampling intervals and transmission delays. Our first goal is to extend the time-delay approach to NCSs under TOD protocol in the presence of communication delays that are allowed to be non-small. This leads to a novel hybrid system model for the closed-loop system, where time-varying delays appear in the dynamics and in the reset equations . Since a similar hybrid system model corresponds to RR protocol, we derive new conditions for Input-To-State (ISS) under RR protocol as well. These conditions are computationally simpler than the existing ones of [16] though may lead to more conservative results. A novel Lyapunov-Krasovskii method is introduced for hybrid delayed systems, which is based on discontinuous in time Lyapunov functionals.

Polytopic type uncertainties in the system model can be easily included in the analysis. The efficiency and advantages of the presented approach are illustrated by two examples. Some preliminary results were presented in [17].

Notation: Throughout the paper, the superscript ' T ' stands for matrix transposition, R n denotes the n dimensional Euclidean space with vector norm | · | , R n × m is the set of all n × m real matrices, and the notation P &gt; 0 , for P ∈ R n × n means that P is symmetric and positive definite. The symmetric elements of the symmetric matrix will be denoted by ∗ , λ min ( P ) denotes the smallest eigenvalue of matrix P . The space of functions φ : [ -τ M , 0] → R n , which are absolutely continuous on [ -τ M , 0] , and have square integrable first order derivatives is denoted by W [ -τ M , 0] with the norm ‖ φ ‖ W = max θ ∈ -[ τ M , 0] | φ θ ( ) | +

[ ∫ 0 -τ M | ˙ ( φ s ) | 2 ds ] 1 2 . Z + denotes the set { 0, 1, 2, . . . }, whereas N denotes the natural numbers. The symbol x t denotes x t ( θ ) = x t ( + ) θ , θ ∈ -[ τ M , 0] , whereas

Figure 1: System architecture with N sensors

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0147v1-figure-a553109d-86c6-4dc2-9761-c94d470f5372.jpg)

‖ w t , t [ 0 ] ‖ ∞ stands for the essential supremum of the Euclidean norm | w t [ 0 , t ] | , where w : [ t 0 , t ] → R n w . MATI and MAD denote Maximum Allowable Transmission Interval and Maximum Allowable Delay, respectively.

## 2 Problem formulation

## 2.1 The description of NCS and the hybrid model

Consider the system architecture in Figure 1 with plant

$$\dot { x } ( t ) = A x ( t ) + B u ( t ) + D \omega ( t ), \ \ t \geq 0,$$

where x t ( ) ∈ R n is the state vector, u t ( ) ∈ R m is the control input, ω t ( ) ∈ R q is the essentially bounded disturbance. Assume that there exists a real number ∆ &gt; 0 such that ‖ ω [0 , t ] ‖ ∞ ≤ ∆ for all t ≥ 0 . The system matrices A, B and D can be uncertain with polytopic type uncertainties.

The system has several nodes ( N distributed sensors, a controller node and an actuator node) which are connected via the network. The measurements are given by y i ( ) = t C x t i ( ) ∈ R n i , i = 1 , . . . , N, ∑ N i =1 n i = n y and we denote C = [ C T 1 · · · C T N ] T , y t ( ) = [ y T 1 ( ) t · · · y T N ( ) t ] T ∈ R n y . Let s k denote the unbounded monotonously increasing sequence of sampling instants

$$0 = s _ { 0 } < s _ { 1 } < \dots < s _ { k } < \dots, \ k \in \mathbb { Z } _ { + }, \lim _ { k \to \infty } s _ { k } = \infty.$$

At each sampling instant s k , one of the outputs y i ( s k ) ∈ R n i is transmitted via the sensor network. We suppose that data loss is not possible and that the transmission of the information over the network is subject to a variable delay η k . Then t k = s k + η k is the updating time instant of the Zero-Order Hold (ZOH).

Differently from [6, 15], we do not restrict the network delays to be small with t k = s k + η k &lt; s k +1 , i.e. η k &lt; s k +1 -s k . As in [19] we allow the delay

to be non-small provided that the old sample cannot get to the destination (to the controller or to the actuator) after the most recent one. Assume that the network-induced delay η k and the time span between the updating and the most recent sampling instants are bounded:

$$( 2 ) \quad \ t _ { k + 1 } - t _ { k } + \eta _ { k } \leq \tau _ { M }, \ 0 \leq \eta _ { m } \leq \eta _ { k } \leq M \text{AD}, \ k \in \mathbb { Z } _ { + },$$

where τ M denotes the maximum time span between the time

$$s _ { k } = t _ { k } - \eta _ { k }$$

at which the state is sampled and the time t k +1 at which the next update arrives at the destination. Here η m and MAD are known bounds and τ M = MATI+MAD . Since MATI = τ M -MAD ≤ τ M -η m , η m &gt; τ M 2 implies that the network delays are non-small due to η k ≥ η m &gt; τ M -η m . In the examples of Section 6, we will show that our method is applicable for η m &gt; τ M 2 .

Denote by ˆ( y s k ) = [ ˆ y T 1 ( s k ) · · · ˆ y T N ( s k ) ] T ∈ R n y the output information submitted to the scheduling protocol. At each sampling instant s k , one of the system nodes i ∈ { 1 , . . . , N } is active, that is only one of ˆ ( y i s k ) values is updated with the recent output y i ( s k ) . Let i ∗ k ∈ { 1 , . . . , N } denote the active output node at the sampling instant s k , which will be chosen due to scheduling protocols.

Consider the error between the system output y s ( k ) and the last available information ˆ( y s k -1 ) :

$$e ( t ) = \text{col} \{ e _ { 1 } ( t ), \cdots, e _ { N } ( t ) \} \equiv \hat { y } ( s _ { k - 1 } ) - y ( s _ { k } ), \\ t \in [ t _ { k }, t _ { k + 1 } ), \ k \in \mathbb { Z } _ { + }, \ \hat { y } ( s _ { - 1 } ) \stackrel { \Delta } { \equiv } 0, \ e ( t ) \in \mathbb { R } ^ { n _ { \nu } }. \\ \dashrightarrow \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots \quad \cdots.$$

We suppose that the controller and the actuator are event-driven (in the sense that the controller and the ZOH update their outputs as soon as they receive a new sample). The most recent output information at the controller level is denoted by ˆ( y s k ) .

## Static output feedback control

Assume that there exists a matrix K = [ K 1 · · · K N ] , K i ∈ R m n × i such that A + BKC is Hurwitz. Then, the static output feedback controller has a form

/negationslash where i ∗ k is the index of the active node at s k and η k is communication delay. We obtain thus the impulsive closed-loop model with the following continuous dynamics:

$$( 5 ) \quad u ( t ) = K _ { i _ { k } ^ { * }$$

/negationslash

$$\begin{array} { c c } \dot { x } ( t ) = A x ( t$$

where A 1 = BKC, B i = BK , i i = 1 , . . . , N.

Taking into account (4), we obtain

$$\begin{array} { r c l } e _ { i } ( t _ { k + 1$$

and

/negationslash

$$\begin{array} { r c l } \text{nd} \\ e _ { i } ($$

Thus, the delayed reset system is given by

$$\begin{array} {$$

/negationslash

Therefore, (6)-(7) is the hybrid model of the NCS. Since x t ( k -η k ) = x t ( -τ ( )) t for t ∈ [ t k , t k +1 ) with τ ( ) t = t -t k + η k ∈ [ η m , τ M ] (cf. (2)), the hybrid system model (6)-(7) contains the piecewise-continuous delay τ ( ) t in the continuous-time dynamics (6). Even for η k = 0 we have the delayed state x t ( k ) = x t ( -τ ( )) t with τ ( ) = t t -t k .

Note that the first updating time t 0 corresponds to the time instant when the first data is received by the actuator. Assume that initial conditions for (6)-(7) are given by x t 0 ∈ W [ -τ M , 0] and e t ( 0 ) = -Cx t ( 0 -η 0 ) = -Cx 0 .

## Dynamic output feedback control

Assuming that the controller is directly connected to the actuator, consider a dynamic output feedback controller of the form

$$\dot { x } _ { c } ( t ) = A _ { c } x _ { c } ( t ) + B _ { c } \hat { y } ( s _ { k } ), \\ u ( t ) = C _ { c } x _ { c } ( t ) + D _ { c } \hat { y } ( s _ { k } ), \ t \in [ t _ { k }, t _ { k + 1 } ), \ k \in \mathbb { Z } _ { + },$$

where x c ( ) t ∈ R n c , A , B , C c c c and D c are the matrices with appropriate dimensions. Let e i ( )( t i = 1 , . . . , N ) be defined by (4). The closed-loop system can be presented in the form of (6)-(7), where x , e i and the matrices are changed by the ones with the bars as follows:

$$\begin{array} { l l l } \intertext { the ones with the bars as follows: } \bar { x } = \left [ \begin{array} { c c } x \\ x _ { c } \end{array} \right ], \, \bar { A } = \left [ \begin{array} { c c } A & B C _ { c } \\ 0 _ { n _ { c } \times n } & A _ { c } \end{array} \right ], \, \bar { B } _ { i } = \left [ \begin{array} { c c } B D _ { c } \\ B _ { c } \end{array} \right ], \, \bar { D } = \left [ \begin{array} { c c } D \\ 0 _ { n _ { c } \times q } \end{array} \right ], \\ \bar { A } _ { 1 } = \left [ \begin{array} { c c } B D _ { c } C & 0 _ { n _ { c } \times n _ { c } } \\ B _ { c } C & 0 _ { n _ { c } \times n _ { c } } \end{array} \right ], \, \bar { C } _ { 1 } = \left [ \begin{array} { c c } C _ { 1 } & 0 \\ 0 & 0 \end{array} \right ], \, \bar { C } _ { i } \in \mathbb { R } ^ { n _ { \nu } \times ( n + n _ { c } ) }, \\ \bar { C } _ { 2 } = \left [ \begin{array} { c c } 0 _ { n \times n _ { 1 } } & C _ { 2 } ^ { T } & 0 \\ 0 _ { n \times n _ { 1 } } & 0 & 0 \end{array} \right ], \, \dots, \, \bar { C } _ { N } = \left [ \begin{array} { c c } 0 & C _ { N } ^ { T } \\ 0 & 0 \end{array} \right ] ^ { T }, \\ \bar { e } _ { 1 } ( t ) = [ e _ { 1 } ^ { T } ( t ) \ 0 ] ^ { T }, \, \bar { e } _ { 2 } ( t ) = [ 0 _ { 1 \times n _ { 1 } } & e _ { 2 } ^ { T } ( t ) \ 0 ] ^ { T }, \cdots \, \bar { e } _ { N } ( t ) = [ 0 \, \ e _ { N } ^ { T } ( t ) ] ^ { T }, \, \bar { e } _ { i } ( t ) \in \mathbb { R } ^ { n _ { \nu } }. \\ \text{Remark $2.1 $ In $ [15]$, $ a $ piecewise-continuous $error $e(t)$ = $\hat{y}(t_{k}) - y(t), t \, \in$}$$

Remark 2.1 In [15], a piecewise-continuous error e t ( ) = ˆ( y t k ) -y t , t ( ) ∈ [ s k , s k +1 ] is defined, which leads to the non-delayed continuous dynamics. The

derivation of reset equations is based on the assumption of small communication delays, that is avoided in our approach. In our approach e t ( ) is different: it is given by (4) and is piecewise-constant. As a result, our hybrid model is different with the delayed continuous dynamics. Moreover, in the absence of scheduling protocols, the closed-loop system is given by non-hybrid system (6), where e t ( ) ≡ 0 . The latter is consistent with the time-delay model considered e.g. in [11, 12].

## 2.2 Scheduling protocols

## TOD protocol

In TOD protocol, the output node i ∈ { 1 , . . . , N } with the greatest (weighted) error will be granted the access to the network.

Definition 2.2 (Weighted TOD protocol) Let Q i &gt; 0( i = 1 , . . . , N ) be some weighting matrices. At the sampling instant s k , the weighted TOD protocol is a protocol for which the active output node with the index i ∗ k is defined as any index that satisfies

$$( 8 ) \quad | \sqrt { Q } _ { i _ { k } ^ { * } } e _ { i _ { k } ^ { * } } ( t ) | ^ { 2 } \geq | \sqrt { Q } _ { i } e _ { i } ( t ) | ^ { 2 }, \ t \in [ t _ { k }, t _ { k + 1 } ), \ k \in \mathbb { Z } _ { + }, \ i = 1, \dots, N.$$

A possible choice of i ∗ k is given by

$$i _ { k } ^ { * } = \min \{ \arg \max _ { i \in \{ 1, \dots, N \} } | \sqrt { Q } _ { i } \left ( \hat { y } _ { i } ( s _ { k - 1 } ) - y _ { i } ( s _ { k } ) \right ) | ^ { 2 } \}. \\ \text{conditions for commutino the weiothino matrices } O. \quad O, \text{ will be } c$$

The conditions for computing the weighting matrices Q , . . . , Q 1 N will be given in Theorem 3.2 below.

Remark 2.3 For implementation of TOD protocol in wireless networks we refer to [4].

## RR protocol

The active output node is chosen periodically:

$$\begin{array} { c } i _ { k } ^ { * } = i _ { k + N } ^ { * }, \, f o r \ a l l \, k \in \mathbb { Z } _ { + }, \\ i _ { j } ^ { * } \neq i _ { l } ^ { * }, \, f o r \ 0 \leq j < l \leq N - 1. \end{array}$$

/negationslash

Remark 2.4 Note that another model for the closed-loop system under RR protocol was given in [16]. The model in [16] is a switched system with ordered delays τ 1 ( ) t &lt; · · · &lt; τ N ( ) t , where τ i ( ) t = t -t k -i +1 + η k -i +1 , i = 1 , . . . , N . A Lyapunov-Krasovskii analysis of the latter model is based on the standard time-independent Lyapunov functional for interval delay.

Definition 2.5 The hybrid system (6)-(7) with essentially bounded disturbance ω is said to be partially ISS with respect to x (or x -ISS) if there exist constants b &gt; 0 , δ &gt; 0 and c &gt; 0 such that the following holds for t ≥ t 0

$$| x ( t ) | ^ { 2 } \leq b e ^ { - \delta ($$

for the solutions of the hybrid system initialized with x t 0 = φ ∈ W [ -τ M , 0] and e t ( 0 ) ∈ R n y . The hybrid system (6)-(7) is ISS if additionally the following bound is valid for t ≥ t 0

$$| e ( t ) | ^ { 2 } \leq & \, b e ^ { - \$$

Our objective is to derive Linear Matrix Inequality (LMI) conditions for the partial ISS of the hybrid system (6)-(7) with respect to the variable of interest x . In [5], the notion of partial stability was also used. In Section 3 below, ISS of (6)-(7) under TOD protocol with N sensor nodes will be studied. For N = 2 , less restrictive conditions will be derived in Section 4, and it will be shown that the same conditions guarantee x -ISS of (6)-(7) under RR protocol. In Section 5, the latter conditions will be extended to RR protocol with N ≥ 2 .

## 3 ISS under TOD protocol: general N

/negationslash

Note that the differential equation for x given by (6) depends on e i ( ) t = e i ( t k ) , t ∈ [ t k , t k +1 ) with i = i ∗ k only. Consider the following Lyapunov functional:

$$\begin{array} {$$

where x t ( θ ) ∆ = x t ( + θ , ) θ ∈ [ -τ M , 0] and where we define (for simplicity) x t ( ) = x , 0 t &lt; 0 .

Here the terms

$$e _ { i } ^ { T } ( t ) Q _ { i } e _ { i }$$

are piecewise-constant, ˜ ( V t, x t , x ˙ t ) presents the standard Lyapunov functional for systems with interval delays τ ( ) t ∈ [ η m , τ M ] . The novel piecewise-continuous

in time term V G is inserted to cope with the delays in the reset conditions . It is continuous on [ t k , t k +1 ) and do not grow in the jumps (when t = t k +1 ), since

$$\text{comunous on } \lfloor \iota _ { k }, \i$$

where we applied Jensen's inequality (see e.g., [14]). The function V e ( ) t is thus continuous and differentiable over [ t k , t k +1 ) . The following lemma gives sufficient conditions for the x -ISS of (6)-(8):

Lemma 3.1 Let there exist positive constants α, b, 0 &lt; Q i ∈ R n i × n i , 0 &lt; U i ∈ R n i × n i , 0 &lt; G i ∈ R n i × n i , i = 1 , . . . , N, and V e ( ) t of (10) such that along (6) the following inequality holds

/negationslash

$$& \text{are your own any one} \\ & ( 1 2 ) & \dot { V } _ { e } ( t ) + 2 \alpha V _ { e } ( t ) - \frac { 1 } { \tau _ { M } - \eta _ { m } } \sum _ { i = 1, i \neq i _ { k } ^ { * } } ^ { N } | \sqrt { U } _ { i } e _ { i } ( t ) | ^ { 2 } \\ & - 2 \alpha | \sqrt { Q _ { i _ { k } ^ { * } } } e _ { i _ { k } ^ { * } } ( t ) | ^ { 2 } - b | \omega ( t ) | ^ { 2 } \leq 0, \ t \in [ t _ { k }, t _ { k + 1 } ). \\ & \text{assume additionally that}$$

$$( 1 3 ) \quad \Omega _ { i } \stackrel { \Delta } { = } \left [ \begin{smallmatrix} - \frac { 1 - 2 \alpha ( \tau _ { M } - \eta _ { m } ) } { N - 1 } Q _ { i } + U _ { i } } & Q _ { i } \\ * & Q _ { i } - G _ { i } e ^ { - 2 \alpha \tau _ { M } } \end{smallmatrix} \right ] < 0, \ i = 1, \dots, N. \\ \Gamma _ { m } \quad.. \quad.. \quad.. \quad.. \quad.. \quad. \quad. \quad. \quad. \quad. \quad.$$

Then V e ( ) t does not grow in the jumps along (6)-(8):

/negationslash

Moreover, the following bounds hold for a solution of (6)-(8) initialized by x t 0 ∈ W [ -τ M , 0] , e ( t 0 ) ∈ R n y :

$$\begin{array} { l } \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ ( 1 4 ) & \Theta \stackrel { \Delta } { = } } V _ { e } ( t _ { k + 1 } ) - V _ { e } ( t _ { k + 1 } ^ { - } ) + \sum _ { i = 1, i \neq i _ { k } ^ { * } } ^ { N } | \sqrt { U _ { i } } e _ { i } ( t _ { k } ) | ^ { 2 } \\ & + 2 \alpha ( \tau _ { M } - \eta _ { m } ) | \sqrt { Q _ { i _ { k } ^ { * } } } e _ { i _ { k } ^ { * } } ( t _ { k } ) | ^ { 2 } \leq 0. \end{array} \\ \begin{array} { l } \Theta \equiv \Delta & V _ { e } ( t _ { k + 1 } ) - V _ { e } ( t _ { k + 1 } ^ { - } ) + \sum _ { i = 1, i \neq i _ { k } ^ { * } } ^ { N } | \sqrt { U _ { i } } e _ { i } ( t _ { k } ) | ^ { 2 } \\ & + 2 \alpha ( \tau _ { M } - \eta _ { m } ) | \sqrt { Q _ { i _ { k } ^ { * } } } e _ { i _ { k } ^ { * } } ( t _ { k } ) | ^ { 2 } \leq 0. \end{array} \\ \text{Moreover, the following bounds holds for a solution of } ( 6 ) - ( 8 ) \text{ initialized by } x _ { t } \\ W | _ { - \tau }.. \ \partial | _ { \ e } ( t _ { \wedge } \searrow \mathbb { P } ^ { n _ { u } }. \end{array}$$

$$& V ( t, x _ { t }, \dot { x } _ { t } ) \leq e ^ { - 2 \alpha ( t - t _ { 0 } ) } V _ { e } ( t _ { 0 } ) + \frac { b } { 2 \alpha } \Delta ^ { 2 }, \ t \geq t _ { 0 }, \\ & V _ { e } ( t _ { 0 } ) = V ( t _ { 0 }, x _ { t _ { 0 } }, \dot { x } _ { t _ { 0 } } ) + \sum _ { i = 1 } ^ { N } | \sqrt { Q _ { i } } e _ { i } ( t _ { 0 } ) | ^ { 2 }, \\ & \text{and}$$

$$& \dots \\ & ( 1 6 ) \quad \sum _ { i = 1 } ^ { N } | \sqrt { Q _ { i } } e _ { i } ( t ) | ^ { 2 } \leq \tilde { c } e ^ { - 2 \alpha ( t - t _ { 0 } ) } V _ { e } ( t _ { 0 } ) + \frac { b } { 2 \alpha } \Delta ^ { 2 }, \\ & \text{where } \tilde { c } = e ^ { 2 \alpha ( \tau _ { M } - \eta _ { m } ) }, \text{ simplify ISS of } ( 6 ) \text{-} ( 8 ).$$

where ˜ = c e 2 α τ ( M -η m ) , implying ISS of (6)-(8).

Proof. Since ∫ t t k e -2 α t ( -s ) ds ≤ τ M -η m , t ∈ [ t k , t k +1 ) and | ω t ( ) | ≤ ∆ , by the comparison principle, (12) implies (17)

/negationslash

$$V _ { e } ( t ) & \leq e ^ { - 2 \alpha ( t -$$

Note that (13) yields 0 &lt; 2 α τ ( M -η m ) &lt; 1 and U i ≤ 1 -2 α τ ( M -η m ) N -1 Q i ≤ Q , i i = 1 , . . . , N . Hence,

$$& ( 1 8 ) \quad V ( t, x _ { t }, \dot { x }$$

Since ˜ V | t = t k +1 = ˜ V | t = t -k +1 and e t ( -k +1 ) = e t ( k ) , we obtain

/negationslash

Then taking into account (11) we find

$$\begin{array} {$$

/negationslash

/negationslash

Note that under TOD protocol

$$\begin{array} { l } \text{em using no account} \, ( 1 1 )$$

/negationslash

Denote ζ i = col { e i ( t k ) , C i [ x s ( k ) -x s ( k +1 )] } . Then, employing (7) and (3) we arrive at

$$\dots \dots \dots - \dots - r \dots \dots$$

/negationslash that yields (14).

$$\Theta & \leq - | \sqrt { G _ { i _ { k } ^ { * } } e ^ { - 2 \alpha \tau _ { M } } - Q _ { i _ { k } ^ { * } } } C _ { i _ { k } ^ { * } } [ x ( s _ { k } ) - x ( s _ { k + 1 } ) ] | ^ { 2 } + \sum _ { i = 1, i \neq i _ { k } ^ { * } } ^ { N } \zeta _ { i } ^ { T } \Omega _ { i } \zeta _ { i } & \leq 0, \\ \text{that yields (14).}$$

The inequalities (14) and (17) with t = t -k +1 imply

Then

$$V _ { e } ( t _ { k + 1 } ) \leq e ^ { - 2 \alpha ( t _ { k + 1 } - t _ { k } ) } V _ { e } ( t _ { k } ) + b \Delta ^ { 2 } \int _ { t _ { k } } ^ { t _ { k + 1 } } e ^ { - 2 \alpha ( t _ { k + 1 } - s ) } d s. \\ \text{Then}$$

$$V _ { e } ( t _ { k + 1 } ) & \leq e ^ { - 2 \alpha ( t _ { k + 1 } - t _ { k - 1 } ) } V _ { e } ( t _ { k - 1 } ) + b \Delta ^ { 2 } \int _ { t _ { k - 1 } } ^ { t _ { k + 1 } } e ^ { - 2 \alpha ( t _ { k + 1 } - s ) } d s \\ & \leq e ^ { - 2 \alpha ( t _ { k + 1 } - t _ { 0 } ) } V _ { e } ( t _ { 0 } ) + b \Delta ^ { 2 } \int _ { t _ { 0 } } ^ { t _ { k + 1 } } e ^ { - 2 \alpha ( t _ { k + 1 } - s ) } d s. \\ \underset { \sim \epsilon \, \neq \, \partial \, \left ( \partial \, \times \, \mathfrak { k } \, \underset { \sim \, \partial \, \dots \, \dots } { \dots } } { \dots } }$$

Replacing in (19) k +1 by k and using (18), we arrive at (15), which yields x -ISS of (6)-(8) because

$$\lambda _ { m i n } ( P ) | x ( t ) | ^ { 2 } \leq V ( t, x _ { t }, \dot { x } _ { t } ), \, V ( t _ { 0 }, x _ { t _ { 0 } }, \dot { x } _ { t _ { 0 } } ) \leq \delta \| x _ { t _ { 0 } } \| _ { W } ^ { 2 }$$

for some scalar δ &gt; 0 . Moreover, (19) with k + 1 replaced by k implies (16) since for t ∈ [ t k , t k +1 )

$$e ^ { - 2 \alpha ( t _ { k } - t _ { 0 } ) } = e ^ { - 2 \alpha ( t - t _ { 0 } ) } e ^ { - 2 \alpha ( t _ { k } - t ) } \leq \tilde { c } e ^ { - 2 \alpha ( t - t _ { 0 } ) }.$$

By using Lemma 3.1 and the standard arguments for the delay-dependent analysis, we derive LMI conditions for ISS of (6)-(8) (see Appendix A for the proof):

Theorem 3.2 Given 0 ≤ η m &lt; τ M , α &gt; 0 , assume that there exist positive scalar b , n × n matrices P &gt; 0 , S 0 &gt; 0 , R 0 &gt; 0 , S 1 &gt; 0 , R 1 &gt; 0 , S 12 , n i × n i matrices Q &gt; i 0 , U &gt; i 0 , G &gt; i 0 , i = 1 , . . . , N, such that (13) and the following inequalities are feasible:

$$\Phi = \left [ \begin{array} { c c } R _ { 1 } & S _ { 1 2 } \\ * & R _ { 1 } \end{array} \right ] \geq 0,$$

$$\begin{pmatrix} \Sigma _ { i } - ( F ^ { i } ) ^ { T } \Phi F ^ { i } e ^ { - 2 \alpha \tau _ { M } } & \Xi _ { i } ^ { T } H \\ * & - H \end{pmatrix} < 0, \ i = 1, \dots, N,$$

where (22)

/negationslash

/negationslash

$$\text{where} ( 2 2 ) \\ H = \eta _ { m } ^ { 2 } R _ { 0 } + ( \tau _ { M } - \eta _ { m } ) ^ { 2 } R _ { 1 } + ( \tau _ { M } - \eta _ { m } ) \sum _ { l = 1 } ^ { N } C _ { l } ^ { T } G _ { l } C _ { l }, \\ \Sigma _ { i } = ( F _ { 1 } ^ { i } ) ^ { T } P \Xi _ { i } + ( \Xi ^ { i } ) ^ { T } P F _ { 1 } ^ { i } + \Upsilon _ { i } - ( F _ { 2 } ^ { i } ) ^ { T } R _ { 0 } F _ { 2 } ^ { i } e ^ { - 2 \alpha \eta _ { m } }, \\ F _ { 1 } ^ { i } = [ I _ { n } \ 0 _ { n \times ( 3 n + n _ { y } - n _ { i } + q ) } ], \\ F _ { 2 } ^ { i } = [ I _ { n } \ - I _ { n } \ 0 _ { n \times ( 2 n + n _ { y } - n _ { i } + q ) } ], \\ F _ { 3 } ^ { i } = [ I _ { n } \ - I _ { n } \ 0 _ { n \times n } \ - I _ { n } \ - I _ { n } \ 0 _ { n \times ( n _ { y } - n _ { i } + q ) } \ ], \\ \Xi _ { i } = [ A \ 0 _ { n \times n } \ A _ { 1 } \ 0 _ { n \times n } \ B _ { 2 } \ \cdots B _ { N } \ D ], \ i = 1, \\ \Xi _ { i } = [ A \ 0 _ { n \times n } \ A _ { 1 } \ 0 _ { n \times n } \ B _ { 1 } \ \cdots \ B _ { j _ { j } \neq i } \cdots B _ { N } \ D ], \ i = 2, \dots N, \\ \Upsilon _ { i } = \text{diag} \{ S _ { 0 } + 2 \alpha P, - ( S _ { 0 } - S _ { 1 } ) e ^ { - 2 \alpha \eta _ { m } }, 0 _ { n \times n }, - S _ { 1 } e ^ { - 2 \alpha \tau _ { M } }, \psi _ { 2 }, \cdots, \psi _ { N }, - b I _ { q } \}, \ i = 1, \\ \Upsilon _ { i } = \text{diag} \{ S _ { 0 } + 2 \alpha P, - ( S _ { 0 } - S _ { 1 } ) e ^ { - 2 \alpha \eta _ { m } }, 0 _ { n \times n }, - S _ { 1 } e ^ { - 2 \alpha \tau _ { M } }, \psi _ { 1 } ; \cdots \psi _ { j _ { j } \neq i } ; \cdots \psi _ { N }, - b I _ { q } \}, \\ i = 2, \dots N, \ \psi _ { j } = - \frac { \tau _ { M } - \eta _ { m } } { \tau _ { M } - \eta _ { m } } U _ { j } + 2 \alpha Q _ { j }, \ j = 1, \dots, N. \\ \text{Then solutions of the hybrid system} ( 6 ) - ( 8 ) \text{ satisfy the bound} ( 15 ), \text{where} V ( t, x _ { t }, \dot { x } _ { t } ) \\ \text{is} \ q i v e n b y \ ( 1 0 ), \ i m p l y i n q \ ISS \text{ of} ( 6 ) - ( 8 ). \ \text{ If the above LMIs are feasible with}$$

Then solutions of the hybrid system (6)-(8) satisfy the bound (15), where V ( t, x t , x ˙ t ) is given by (10), implying ISS of (6)-(8). If the above LMIs are feasible with α = 0 , then the bound (15) holds with a small enough α 0 &gt; 0 .

## 4 ISS under TOD/RR protocol: N = 2

For N = 2 less restrictive conditions than those of Theorem 3.2 for the x -ISS of (6)-(7) will be derived via a different from (10) Lyapunov functional:

/negationslash

$$V _ { e } ( t ) & = V ( t, x _ { t }, \dot { x } _ { t } ) + \frac { t _ { k + 1 } - t } { \tau _ { M } - \eta _ { m } } \{ e _ { i } ^ { T } ( t ) Q _ { i } e _ { i } ( t ) \} _ { | i \neq i _ { k } ^ { * } }, \\ Q _ { 1 } & > 0, \ Q _ { 2 } > 0, \ \alpha > 0, \ t \in [ t _ { k }, t _ { k + 1 } ), \ k \in \mathbb { Z } _ { + }, \\. \quad. \quad. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

where i ∗ k ∈ { 1 2 , } and V ( t, x , x t ˙ t ) is given by (10) with G i = Q e i 2 ατ M . The term t k +1 -t τ M -η m { e T i ( t k ) Q e i i ( t k ) } is inspired by the similar construction of Lyapunov functionals for the sampled-data systems [8, 18, 22]. The following statement holds:

Lemma 4.1 Given N = 2 , if there exist positive constants α, b and V e ( ) t of (23) such that along (6)-(8) ((6), (7), (9)) the following inequality holds

$$\dot { V } _ { e } ( t ) + 2 \alpha V _ { e } ( t ) - b | \omega ( t ) | ^ { 2 } \leq 0, \ \ t \in [ t _ { k }, t _ { k + 1 } ).$$

Then V e ( ) t does not grow in the jumps along (6)-(8) ((6), (7), (9)), where

$$\Theta & \stackrel { \Delta } { = } V _ { e } ( t _ { k + 1 } ) - V _ { e } ( t _ { k + 1 } ^ { - } ) \leq 0.$$

The bound (15) is valid for a solution of (6)-(8) ((6), (7), (9)) with the initial condition x t 0 ∈ W [ -τ M , 0] , e t ( 0 ) ∈ R n y , implying the x -ISS of (6)-(8) ((6), (7), (9)).

Proof. Since | ω t ( ) | ≤ ∆ , (24) implies

$$& V _ { e } ( t ) \leq e ^ { - 2 \alpha ( t - t _ { k } ) } V _ { e } ( t _ { k } ) + b \Delta ^ { 2 } \int _ { t _ { k } } ^ { t } e ^ { - 2 \alpha ( t - s ) } d s, \ t \in [ t _ { k }, t _ { k + 1 } ). \\ & \text{Noting that}$$

Noting that

/negationslash we obtain employing (11)

$$V _ { e } ( t _ { k + 1 } ) & \leq \tilde { V } _ { | t = t _ { k + 1 } } + | \sqrt { Q } _ { i } e _ { i } ( t _ { k + 1 } ) | _ { | i \neq i _ { k + 1 } ^ { * } } \\ & \quad + \sum _ { i = 1 } ^ { 2 } ( \tau _ { M } - \eta _ { m } ) \int _ { t _ { k + 1 } - \eta _ { k + 1 } } ^ { t _ { k + 1 } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } | \sqrt { G _ { i } } C _ { i } \dot { x } ( s ) | ^ { 2 } d s, \\ \text{we obtain employing (11)}$$

/negationslash

/negationslash

We will prove that Θ ≤ 0 under TOD and RR protocols, respectively. Under TOD protocol we have

$$\Theta & \leq e _ { i } ^ { T } ( t _ { k + 1 } ) Q _ { i } e _ { i } ( t _ { k + 1 } ) _ { | i \neq i _ { k + 1 } ^ { * } } + V _ { G | t = t _ { k + 1 } } - V _ { G | t = t _ { k + 1 } ^ { - } } \\ & \leq e _ { i } ^ { T } ( t _ { k + 1 } ) Q _ { i } e _ { i } ( t _ { k + 1 } ) _ { | i \neq i _ { k + 1 } ^ { * } } - \sum _ { i = 1 } ^ { 2 } | \sqrt { Q _ { i } } C _ { i } [ x ( t _ { k } - \eta _ { k } ) - x ( t _ { k + 1 } - \eta _ { k + 1 } ) ] | ^ { 2 }. \\ \text{We will prove that } \Theta \leq 0 \text{ under TOD and RR protocols, respectively. } \text{ Under}$$

/negationslash

$$e _ { i } ^ { T } ( t _ { k + 1 } ) Q _ { i } e _ { i } ( t _ { k + 1 } ) _ { | i \neq i _ { k + 1 } ^ { * } } \leq e _ { i _ { k } ^ { * } } ^ { T } ( t _ { k + 1 } ) Q _ { i _ { k } ^ { * } } e _ { i _ { k } ^ { * } } ( t _ { k + 1 } )$$

for i ∗ k +1 = i ∗ k , whereas

/negationslash

$$e _ { i } ^ { T } ( t _ { k + 1 } ) Q _ { i } e _ { i } ( t _ { k + 1 } ) _ { | i \neq i _ { k + 1 } ^ { * } } = e _ { i _ { k } ^ { * } } ^ { T } ( t _ { k + 1 } ) Q _ { i _ { k } ^ { * } } e _ { i _ { k } ^ { * } } ( t _ { k + 1 } )$$

for i ∗ k +1 = i ∗ k . Then, taking into account (7) we obtain

$$\L k + 1 / \L k \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \\ \Theta & \leq | \sqrt { Q _ { i } } C _ { i } [ x ( t _ { k } - \eta _ { k } ) - x ( t _ { k + 1 } - \eta _ { k + 1 } ) ] | _ { | i = i _ { k } ^ { * } } ^ { 2 } \\ & - \sum _ { i = 1 } ^ { 2 } | \sqrt { Q _ { i } } C _ { i } [ x ( t _ { k } - \eta _ { k } ) - x ( t _ { k + 1 } - \eta _ { k + 1 } ) ] | ^ { 2 } \leq 0. \\ \intertext {ler RR protocol we have } i _ { k + 1 } ^ { * } \neq i _ { k } ^ { * } \text{ meaning that } ( 2 7 ) \text{ holds and that } \Theta \leq \\ \text{$w$ to $real$ for $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $geq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq$neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neeq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $peq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq  neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq e $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $ neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $ne$neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $pe$neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $ge$neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $ neeq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq$peq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $ne$geq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $ne$peq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $peq$neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq$ne$neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq $neq$$

/negationslash

/negationslash

Under RR protocol we have i ∗ k +1 = i ∗ k meaning that (27) holds and that Θ ≤ 0 . Then the result follows by the arguments of Lemma 3.1.

Remark 4.2 Differently from Lemma 3.1, Lemma 4.1 guarantees (19) that does not give a bound on e i ∗ k ( t k ) since V e ( ) t for t ∈ [ t k , t k +1 ) does not depend on e i ∗ k ( t k ) . That is why Lemma 4.1 guarantees only x -ISS. However, as explained in Remark 5.3 below, under RR protocol x -ISS implies boundedness of e .

In the next section, we will extend the result of Lemma 4.1 under RR protocol to the case of N ≥ 2 . Theorem 5.2 below (in the particular case of N = 2 ) will provide LMIs for the x -ISS of (6)-(8) ((6), (7), (9)).

## 5 ISS under RR protocol: N ≥ 2

Under RR protocol (9), the reset system (7) can be rewritten as

$$\begin{array} { c } x ( t _ { k + 1 } ) = x ( t _ { k + 1 } ^ { - } ), \\ e _ { i _ { k - j } ^ { * } } ( t _ { k + 1 } ) = C _ { i _ { k - j } ^ { * } } [ x ( s _ { k - j } ) - x ( s _ { k + 1 } ) ], \\ j = 0, \dots, N - 1 \text{ if } k \geq N - 1, \end{array}$$

where the index k -j corresponds to the last updated measurement in the node i ∗ k -j .

Consider the following Lyapunov functional:

$$\begin{array} { c } V _ { e } ( t ) = V ( t, x _ { t }, \dot { x } _ { t } ) + V _ { Q }, \quad t \geq t _ { N - 1 }, \\ V ( t, x _ { t }, \dot { x } _ { t } ) = \tilde { V } ( t, x _ { t }, \dot { x } _ { t } ) + V _ { G }, \end{array}$$

where ˜ ( V t, x t , x ˙ t ) is given by (10). The discontinuous in time terms V Q and V G are defined as follows:

(30)

$$& ( 3 0 ) \\ & V _ { Q } = \sum _ { j = 1 } ^ { N - 1 } \frac { t _ { k + 1 } - t } { j ( \tau _ { M } - \eta _ { m } ) } | \sqrt { Q _ { i _ { k - j } ^ { * } } } e _ { i _ { k - j } ^ { * } } ( t ) | ^ { 2 }, \ k \geq N - 1, \ t \in [ t _ { k }, t _ { k + 1 } ), \\ & V _ { G } = \begin{cases} \sum _ { i = 1 } ^ { N } ( \tau _ { M } - \eta _ { m } ) \int _ { s _ { k } } ^ { t } e ^ { 2 \alpha ( s - t ) } | \sqrt { G _ { i } } C _ { i } \dot { x } ( s ) | ^ { 2 } d s, \ k \geq N, \ t \in [ t _ { k }, t _ { k + 1 } ), \\ \sum _ { i = 1 } ^ { N } ( \tau _ { M } - \eta _ { m } ) \int _ { s _ { 0 } } ^ { t } e ^ { 2 \alpha ( s - t ) } | \sqrt { G _ { i } } C _ { i } \dot { x } ( s ) | ^ { 2 } d s, \ t \in [ t _ { N - 1 }, t _ { N } ), \end{cases} \\ & \text{where for } i = 1, \dots, N$$

where for i = 1 , . . . , N

$$G _ { i } = ( N - 1 ) Q _ { i } e ^ { 2 \alpha$$

Here V e does not depend on e i ∗ k ( t k ) . Note that given i = 1 , . . . , N , e i -term appears N -1 times in V Q for every N intervals [ t k + j , t k + +1 j ) , j = 0 , . . . , N -1 (except of the interval with i ∗ k + j = i ). This motivates N -1 in (31) because V G is supposed to compensate V Q term.

As in the previous sections, the term V G is inserted to cope with the delays in the reset conditions. It is continuous on [ t k , t k +1 ) and does not grow in the jumps (when t = t k +1 ), since for k &gt; N -1 (cf. (11)) (32)

and for k = N -1

$$V _ { G | t = t _ { k + 1 } } - V _ { G | t$$

$$\dots \dots \dots \dots \dots \dots \dots \$$

The term V Q grows in the jumps as follows:

$$\begin{array} { r c l } \kappa \, - & & \cdot &$$

where we have used Jensen's inequality and the bound

$$\begin{array} {$$

Since 1 ≤ e 2 α τ [ M +( N -2)( τ M -η m )] e 2 α s ( k -j -t k +1 ) for j = 0 , . . . , N -2 , we obtain

$$V _ { Q | _ { t = t _ { k + 1 } } } - V _ { Q | _ { t = t _ { k + 1 } } } & \leq \sum _ { j = 0 } ^ { N - 2 } ( \tau _ { M } - \eta _ { m } ) e ^ { 2 \alpha [ \tau _ { M } + ( N - 2 ) ( \tau _ { M } - \eta _ { m } ) ] } \\ & \quad \times \int _ { s _ { k - j } } ^ { s _ { k + 1 } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } | \sqrt { Q _ { i _ { k - j } ^ { * } } } C _ { i _ { k - j } ^ { * } } \dot { x } ( s ) | ^ { 2 } d s. \\ \text{The following lemma gives sufficient conditions for the $x$-ISS of $(6),(9),(28)$}.$$

The following lemma gives sufficient conditions for the x -ISS of (6), (9), (28) (see Appendix B for proof):

Lemma 5.1 If there exist positive constants α, b and V e ( ) t of (29) such that along (6), (9), (28) the inequality (24) is satisfied for k ≥ N -1 . Then the following bound holds along the solutions of (6), (9), (28):

$$V _ { e } ( t _ { k + 1 } ) & \leq e ^ { - 2 \alpha ( t _ { k + 1 } - t _ { N - 1 } ) } V _ { e } ( t _ { N - 1 } ) + \Psi _ { k + 1 } \\ & \quad + b \Delta ^ { 2 } \int _ { t _ { N - 1 } } ^ { t _ { k + 1 } } e ^ { - 2 \alpha ( t _ { k + 1 } - s ) } d s, \ k \geq N - 1, \\ \text{where}$$

Moreover, for all t ≥ t N -1

$$\vdash \vdash \ J _ { t _ { N - 1 } } \lambda \lambda \lambda \nu, \nu \varepsilon - \lambda \nu, \\ \text{where} \\ ( 3 7 ) \\ \Psi _ { k + 1 } = - ( \tau _ { M } - \eta _ { m } ) e ^ { 2 \alpha [ \tau _ { M } + ( N - 2 ) ( \tau _ { M } - \eta _ { m } ) ] } \\ \times & \left [ \sum _ { j = 0 } ^ { N - 3 } ( N - 2 - j ) \int _ { s _ { k - j - 1 } } ^ { s _ { k + 1 } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } | \sqrt { Q _ { i _ { k - j } } } C _ { i _ { k - j } ^ { * } } \dot { x } ( s ) | ^ { 2 } d s \\ + ( N - 1 ) \int _ { s _ { k } } ^ { s _ { k + 1 } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } | \sqrt { Q _ { i _ { k + 1 } } } C _ { i _ { k + 1 } ^ { * } } \dot { x } ( s ) | ^ { 2 } d s \right ] \leq 0. \\ \text{Moreover}, \text{ for all } t \geq t _ { N - 1 } \\ \cdot \quad \cdot \quad \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \Lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \lambda \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \phi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \chi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \xi \x \,.$$

$$\begin{array} { c } V ( t, x _ { t }, \dot { x } _ { t } ) \leq e ^ { - 2 \alpha ( t - t _ { N - 1 } ) } V _ { e } ( t _ { N - 1 } ) + \frac { b } { 2 } \Delta ^ { 2 }, \\ V _ { e } ( t _ { N - 1 } ) = V ( t _ { N - 1 }, x _ { t _ { N - 1 } }, \dot { x } _ { t _ { N - 1 } } ) + \sum _ { i = 1 } ^ { N } | \sqrt { Q _ { i } } e _ { i } ( t _ { N - 1 } ) | ^ { 2 }. \\ \intertext { The latter inequality guarantees the $x$-ISS of $(6), $(9)$, $(28)$ for $t\geq t_{N-1}$.} } \end{array}$$

The latter inequality guarantees the x -ISS of (6), (9), (28) for t ≥ t N -1 .

By using Lemma 5.1, arguments of Theorem 3.2 and the fact that for j = 1 , . . . N -1

$$\frac { d } { d t } \frac { t _ { k + 1 } - t } { j ( \tau _ { M } - \eta _ { m } ) } = - \frac { 1 } { j ( \tau _ { M } - \eta _ { m } ) } \leq - \frac { 1 } { ( N - 1 ) ( \tau _ { M } - \eta _ { m } ) }, \\ \omega = \omega \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ a \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ o n \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ u n \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e \ e e$$

we arrive at the the following result:

Theorem 5.2 Given 0 ≤ η m &lt; τ M and α &gt; 0 , assume that there exist positive scalar b , n × n matrices P &gt; 0 , S 0 &gt; 0 , R 0 &gt; 0 , S 1 &gt; 0 , R 1 &gt; 0 , S 12 and n i × n i matrices Q &gt; i 0 ( i = 1 , . . . , N ) such that (20) and (21) are feasible with U i = Q i N -1 , where G i is given by (31) . Then for N &gt; 2 solutions of the hybrid system (6), (9), (28) satisfy the bound (38) with V ( t, x t , x ˙ t ) given by (29), meaning x -ISS for t ≥ t N -1 . For N = 2 solutions of the hybrid system (6)-(8) ((6), (7), (9)) satisfy the bound (15) meaning x -ISS (for t ≥ t 0 ). Moreover, if the above LMIs are feasible with α = 0 , then the solution bounds hold with a small enough α 0 &gt; 0 .

Remark 5.3 For N = 2 and α = 0 , the LMIs of Theorem 3.2 are more restrictive than those of Theorem 5.2: (13) of Theorem 3.2 yields ( N -1) U &lt; Q &lt; G i i i , whereas in Theorem 5.2 we have ( N -1) 2 U i = ( N -1) Q i = G i that leads to larger U i for the same G i . The latter helps for the feasibility of (21) , where U i &gt; 0 appears on the main diagonal only (with minus). However, Theorem 3.2 achieves ISS with respect to the full state col { x, e } and provides the solution bound for t ≥ t 0 , while Theorem 5.2 guarantees only x -ISS.

Note that Theorem 5.2 under RR protocol guarantees boundedness of e as well. Indeed, since e t ( N ) in (28) depends on x (0) , . . . , x ( t N -η N ) and t N ≤ Nτ M (this can be verified similar to (34)), relations (28) yield

$$| e _ { i } ( t ) | ^ { 2 } \leq c ^ { \prime } s u p _ { \theta \in [ - N \tau _ { M }, 0 ] } | x ( t + \theta ) | ^ { 2 }, \ t \geq t _ { N }$$

with some c ′ &gt; 0 , which together with (38) imply

$$| e _ { i } ( t ) | ^ { 2 } \leq c ^ { \prime \prime } [ e ^ { - 2 \alpha ( t - t _ { N - 1 } ) } V _ { e } ( t _ { N - 1 } ) + \Delta ^ { 2 } ]$$

for some c ′′ &gt; 0 and all t ≥ t N + Nτ M .

Remark 5.4 The LMIs of Theorems 3.2 and 5.2 are affine in the system matrices. Therefore, in the case of system matrices from an uncertain time-varying polytope

$$\begin{array} { c } \Omega = \sum _ { j = 1 } ^ { M } g _ { j } ( t ) \Omega _ { j }, \quad 0 \leq g _ { j } ( t ) \leq 1, \\ \sum _ { j = 1 } ^ { M } g _ { j } ( t ) = 1, \quad \Omega _ { j } = \left [ \ A ^ { ( j ) } \quad B ^ { ( j ) } \quad D ^ { ( j ) } \ \right ], \\ e \, \text{ to solve these LMIs simulataneously for all the $M$ vertices $\Omega_{j}$,} \\ e \, \text{ decision matrices.} \end{array}$$

one have to solve these LMIs simultaneously for all the M vertices Ω j , applying the same decision matrices.

## 6 Examples

## 6.1 Example 1: uncertain inverted pendulum

Consider an inverted pendulum mounted on a small car. We focus on the stability analysis in the absence of disturbance. Following [13], we assume that the friction coefficient between the air and the car, f c , and the air and the bar, f b , are not exactly known and time-varying: f c ( ) t ∈ [0 15 . , 0 25] . and f b ( ) t ∈ [0 15 . , 0 25] . . The linearized model can be written as (1), where the matrices A = E -1 A f and B = E -1 B 0 are determined from

$$A = E ^ { - \ell } A _ { f } \text{ and } B = E ^ { - \ell } B _ { 0 } \text{ are determined from} \\ E = \left [ \begin{array} { c c c c } 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 3 / 2 & - 1 / 4 \\ 0 & 0 & - 1 / 4 & 1 / 6 \end{array} \right ], \\ A _ { f } = \left [ \begin{array} { c c c c } 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & - ( f _ { c } + f _ { b } ) & f _ { b } / 2 \\ 0 & 5 / 2 & f _ { b } / 2 & - f _ { b } / 3 \end{array} \right ] \quad \text{and} \quad B _ { 0 } = \left [ \begin{array} { c c c c } 0 \\ 0 \\ 1 \\ 0 \end{array} \right ].$$

Here A belongs to uncertain polytope, defined by four vertices corresponding to f /f c b = 0 15 . and f /f c b = 0 25 . . The pendulum can be stabilized by a state feedback u t ( ) = Kx t ( ) , where x = [ x , x 1 2 , x 3 , x 4 ] T , with the gain

$$K = [ 1 1. 2 0 6 2 \ - 1 2 8. 8 5 9 7 \ 1 0. 7 8 2 3 \ - 2 2. 2 6 2 9 ].$$

In this model, x 1 and x 2 represent cart position and velocity, whereas x , x 3 4 represent pendulum angle from vertical and its angular velocity respectively. In practice x , x 1 2 and x 3 , x 4 (presenting spatially distributed components of the state of the pendulum-cart system) are not accessible simultaneously. Suppose that the state variables are not accessible simultaneously. Consider first N = 2 and

$$C _ { 1 } = \left [ \begin{array} { c c c } 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \end{array} \right ], \, C _ { 2 } = \left [ \begin{array} { c c c } 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{array} \right ].$$

The applied controller gain K has the following blocks:

$$K _ { 1 } = \left [ \ 1 1. 2 0 6 2 \ - 1 2 8. 8 5 9 7 \ \right ], \ K _ { 2 } = \left [ \ 1 0. 7 8 2 3 \ - 2 2. 2 6 2 9 \ \right ]. \\ \text{For the values of } \eta _ { m } \text{ given in Table 1, we apply Theorems 3.2 and 5.2}$$

For the values of η m given in Table 1, we apply Theorems 3.2 and 5.2 with α = 0 , b = 0 via Remark 5.4 and find the maximum values of τ M = MATI+MAD that preserve the stability of the hybrid system (6)-(7) with ω t ( ) = 0 with respect to x . From Table 1, it is observed that under TOD or RR protocol the conditions of Theorem 5.2 possess less decision variables, and stabilize the system for larger τ M than the results in [16] under RR protocol. Moreover, when η m &gt; τ M 2 ( η m = 0 02 0 04) . , . , our method is still feasible (communication delays are larger than the sampling intervals). The computational time for solving the LMIs (in seconds) under the TOD protocol is essentially less than that under RR protocol in [16] (till 36% decrease).

Table 1: Example 1 (N=2): max. value of τ M = MATI+MAD

| τ M \ η m            |     0 |   0.005 |   0.01 |   0.02 |   0.04 |   Decision variables |
|----------------------|-------|---------|--------|--------|--------|----------------------|
| [16] (RR)            | 0.023 |   0.026 |  0.029 |  0.035 |  0.046 |                  146 |
| Theorem 3.2 (TOD)    | 0.014 |   0.018 |  0.021 |  0.029 |  0.044 |                   84 |
| Theorem 5.2 (TOD/RR) | 0.025 |   0.028 |  0.031 |  0.036 |  0.047 |                   72 |

Consider next N = 4 , where C , . . . , C 1 4 are the rows of I 4 and K ,. . . , K 1 4 are the entries of K given by (39). Here the maximum values of τ M that preserve the stability of (6)-(7) with ω t ( ) = 0 with respect to x are given in Table 2. Also here Theorem 5.2 leads to less conservative results than Theorem 3.2.

Table 2: Example 1 (N=4): max. value of τ M = MATI+MAD

| τ M \ η m         |     0 |   0.01 |
|-------------------|-------|--------|
| Theorem 3.2 (TOD) | 0.003 |  0.012 |
| Theorem 5.2 (RR)  | 0.006 |  0.015 |

## 6.2 Example 2: batch reactor

We illustrate the efficiency of the given conditions on the example of a batch reactor under the dynamic output feedback (see e.g., [15]), where N = 2 and

$$A = \left [ \begin{array} { c c c c } 1. 3 8 0 & - 0. 2 0 8 & 6. 7 1 5 & - 5. 6 7 6 \\ - 0. 5 8 1 & - 4. 2 9 0 2 & 0 & 0. 6 7 5 \\ 1. 0 6 7 & 4. 2 7 3 & - 6. 6 5 4 & 5. 8 9 3 \\ 0. 0 4 8 & 4. 2 7 3 & 1. 3 4 3 & - 2. 1 0 4 \end{array} \right ],$$

For the values of η m given in Table 3, we apply Theorems 3.2 and 5.2 with α = 0 , b = 0 and find the maximum values of τ M = MATI+MAD that preserve the stability of the hybrid system (6)-(7) with ω t ( ) = 0 with respect to x . From Table 3 it is seen that the results of our method essentially improve the results in [15], and are more conservative than those obtained via the discrete-time approach. Recently in [2] the same result τ M = 0 035 . as ours in Theorem 5.2 for η m = 0 , MAD = 0 01 . has been achieved. In [2] the sum of squares method is developed in the framework of hybrid system approach. We note that the sum of squares method has not been applied yet to ISS. Moreover, our conditions are simple LMIs with a fewer decision variables. When η m &gt; τ M 2 ( η m = 0 03 0 04) . , . , our method is still feasible (communication delays are larger than the sampling intervals). The computational time under the TOD protocol is essentially less than that under RR protocol in [16] (till 32% decrease).

$$A = \left [ \begin{array} { c c c c } 1. 3 8 0 & - 0. 2 8 & 6. 7 1 5 & - 5. 6 7 6 \\ - 0. 5 8 1 & - 4. 2 9 0 2 & 0 & 0. 6 7 5 \\ 1. 0 6 7 & 4. 2 7 3 & - 6. 6 5 4 & 5. 8 9 3 \\ 0. 0 4 8 & 4. 2 7 3 & 1. 3 4 3 & - 2. 1 0 4 \end{array} \right ], \\ B = \left [ \begin{array} { c c c c } 0 & 0 \\ 5. 6 7 9 & 0 \\ 1. 1 3 6 & - 3. 1 4 6 \\ 1. 1 3 6 & 0 \end{array} \right ], \, C = \left [ \frac { C _ { 1 } } { C _ { 2 } } \right ] = \left [ \frac { 1 } { 0 } & 1 & - 1 \\ 0 & 1 & 0 & 0 \end{array} \right ], \\ \left [ \frac { A _ { c } } { C _ { c } } \left | \ B _ { c } \\ \frac { C _ { c } } { C _ { c } } \right | \ D _ { c } \right ] = \left [ \frac { 0 } { - 2 } & 0 & 1 & 0 \\ 0 & 8 & 5 & 0 \end{array} \right ].$$

## 7 Conclusions

In this paper, a time-delay approach has been developed for the ISS of NCS with scheduling protocols, variable transmission delays and variable sampling intervals. A novel hybrid system model with time-varying delays in the continuous dynamics and in the reset equations is introduced and a new Lyapunov-

Table 3: Example 2: max. value of τ M = MATI+MAD for different η m

| τ M \ η m             |      0 |   0.004 | 0.01   | 0.02   | 0.03   | 0.04   |
|-----------------------|--------|---------|--------|--------|--------|--------|
| [15]( MAD = 0 . 004 ) | 0.0108 |  0.0133 | -      | -      | -      | -      |
| [6]( MAD = 0 . 03 )   | 0.069  |  0.069  | 0.069  | 0.069  | 0.069  | -      |
| Theorem 3.2 (TOD)     | 0.019  |  0.022  | 0.027  | 0.034  | 0.042  | 0.050  |
| Theorem 5.2 (TOD/RR)  | 0.035  |  0.037  | 0.041  | 0.047  | 0.053  | 0.059  |
| [16] (RR)             | 0.042  |  0.044  | 0.048  | 0.053  | 0.058  | 0.063  |

Krasovskii method is developed. The ISS conditions of the delayed hybrid system are derived in terms of LMIs. Differently from the existing (hybrid and discrete-time) methods on the stabilization of NCS with scheduling protocols, the time-delay approach allows non-small network-induced delay (which is not smaller than the sampling interval). Future work will involve consideration of more general NCS models, including packet dropouts, packet disordering, quantization and scheduling protocols for the actuator nodes.

## A Proof of Theorem 3.2

Proof. Consider t ∈ [ t k , t k +1 ) , k ∈ Z + and define ξ i ( ) t = col { x t , x t ( ) ( -η m ) , x ( t -τ ( )) t , x ( t -τ M ) , e 1 ( ) t , · · · , e j ( ) t , · · · , e N ( ) t , ω ( ) t } with i = i ∗ k ∈ N , j = i . Differentiating V e ( ) t along (6) and applying Jensen's inequality, we have

/negationslash

$$j \neq i. \ D i f e r e n t i a l i n g \, V _ { e } ( t ) \, \text{ along} \, ( 6 ) \, \text{ and applying Jensen's inequality, we have} \\ \eta _ { m } \int _ { t - \eta _ { m } } ^ { t } \dot { x } ^ { T } ( s ) R _ { 0 } \dot { x } ( s ) d s \geq \int _ { t - \eta _ { m } } ^ { t } \dot { x } ^ { T } ( s ) d s R _ { 0 } \int _ { t - \eta _ { m } } ^ { t } \dot { x } ( s ) d s \\ = \xi _ { i } ^ { T } ( t ) ( F _ { 2 } ^ { i } ) ^ { T } R _ { 0 } F _ { 2 } ^ { i } \xi _ { i } ( t ), \\ - ( \tau _ { M } - \eta _ { m } ) \int _ { t - \tau _ { M } } ^ { t - \eta _ { m } } \dot { x } ^ { T } ( s ) R _ { 1 } \dot { x } ( s ) d s \\ = - ( \tau _ { M } - \eta _ { m } ) \int _ { t - \tau ( t ) } ^ { t - \eta _ { m } } \dot { x } ^ { T } ( s ) R _ { 1 } \dot { x } ( s ) d s - ( \tau _ { M } - \eta _ { m } ) \int _ { t - \tau _ { M } } ^ { t - \tau ( t ) } \dot { x } ^ { T } ( s ) R _ { 1 } \dot { x } ( s ) d s \\ \leq - \frac { \tau _ { M } - \eta _ { m } } { \tau ( t ) - \eta _ { m } } \xi _ { i } ^ { T } ( t ) \left [ \left [ I _ { n } \, 0 _ { n \times n } \right ] F ^ { i } \right ] ^ { T } R _ { 1 } [ I _ { n } \, 0 _ { n \times n } ] F ^ { i } \xi _ { i } ( t ) \\ - \frac { \tau _ { M } - \eta _ { m } } { \tau _ { M } - \tau ( t ) } \xi _ { i } ^ { T } ( t ) \left [ [ 0 _ { n \times n } \, I _ { n } ] F ^ { i } \right ] ^ { T } R _ { 1 } [ 0 _ { n \times n } \, I _ { n } ] F ^ { i } \xi _ { i } ( t ) \\ \leq - \xi _ { i } ^ { T } ( t ) ( F ^ { i } ) ^ { T } \Phi F ^ { i } \xi _ { i } ( t ). \\ \text{The latter inequality holds if} \, ( 2 0 ) \, \text{ is feasible} \, [ 2 1 ]. \ \text{ Then}$$

The latter inequality holds if (20) is feasible [21]. Then

/negationslash

$$\begin{array} {$$

if Σ + Ξ i T i H Ξ i -( F i ) T Φ F e i -2 ατ M &lt; 0 , i.e., by Schur complement, if (21) is feasible. Thus due to Lemma 3.1, inequalities (13), (20) and (21) imply (15).

## B Proof of Lemma 5.1

Proof. Since | ω t ( ) | ≤ ∆ , (24) implies

$$( 4 0 ) \quad & V ( t, x _ { t }, \dot { x }$$

Note that V e ( t k +1 ) = ˜ V | t = t k +1 + V Q | t = t k +1 + V G | t = t k +1 . Taking into account (28) and the relations ˜ V t = t k +1 = ˜ V t = t -k +1 , e t ( -k +1 ) = e t ( k ) , we obtain due to (31), (32) and (35) for k &gt; N -1 (41)

$$( 4 1 ) \overset { \cdot } { \equiv } & V _ { e$$

whereas for k = N -1 due to (33) and (35)

$$\text{whereas for $k=1\nu-1$ use $\omega(\,\$$

We will prove (36) by induction. For k = N -1 we have

$$\cdots \cdots \cdots \cdots \cdots \cdots \cdots \$$

which implies (36).

Assume that (36) holds for k -1 ( k ≥ N -1 ):

$$V _ { e } ( t _ { k } ) \leq e ^ { - 2 \alpha ( t _ { k } - t _ { N - 1 } ) } V _ { e } ( t _ { N - 1 } ) + \Psi _ { k } + b \Delta ^ { 2 } \int _ { t _ { N - 1 } } ^ { t _ { k } } e ^ { - 2 \alpha ( t _ { k } - s ) } d s. \\ \text{Then due to } ( 4 0 ) \text{ for } t = t _ { k + 1 } ^ { - } \text{ we obtain}$$

Then due to (40) for t = t -k +1 we obtain

$$V _ { e } ( t _ { k + 1 } ) & \leq \Theta _ { k + 1 } + e ^ { - 2 \alpha ( t _ { k + 1 } - t _ { k } ) } \Psi _ { k } + e ^ { - 2 \alpha ( t _ { k + 1 } - t _ { N - 1 } ) } V _ { e } ( t _ { N - 1 } ) \\ & \quad + b \Delta ^ { 2 } \int _ { t _ { N - 1 } } ^ { t _ { k + 1 } } e ^ { - 2 \alpha ( t _ { k + 1 } - s ) } d s.$$

We have

$$\text{We have} \\ e ^ { - 2 \alpha ( t _ { k + 1 } - t _ { k } ) } \Psi _ { k } = - ( \tau _ { M } - \eta _ { m } ) e ^ { 2 \alpha [ \tau _ { M } + ( N - 2 ) ( \tau _ { M } - \eta _ { m } ) ] } \\ \times & \left [ \sum _ { l = 0 } ^ { N - 3 } ( N - 2 - l ) \int _ { s _ { k - l - 2 } } ^ { s _ { k } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } \times | \sqrt { Q _ { i _ { k - 1 } - l } } C _ { i _ { k - 1 } - l } ^ { * } \dot { x } ( s ) | ^ { 2 } d s \right ] \\ + & ( N - 1 ) \int _ { s _ { k - 1 } } ^ { s _ { k } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } | \sqrt { Q _ { i _ { k } } ^ { * } } C _ { i _ { k } ^ { * } } \dot { x } ( s ) | ^ { 2 } d s \right ] \\ = & - ( \tau _ { M } - \eta _ { m } ) e ^ { 2 \alpha [ \tau _ { M } + ( N - 2 ) ( \tau _ { M } - \eta _ { m } ) ] } \left [ \sum _ { j = 0 } ^ { N - 2 } ( N - 1 - j ) \\ \times & \int _ { s _ { k - j - 1 } } ^ { s _ { k } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } | \sqrt { Q _ { i _ { k - j } } ^ { * } } C _ { i _ { k - j } ^ { * } } \dot { x } ( s ) | ^ { 2 } d s \right ]. \\ \text{Then, taking into account } ( 4 1 ), we find$$

Then, taking into account (41), we find

$$& \text{Then, taking into account (41), we find} \\ & \Theta _ { k + 1 } + e ^ { - 2 \alpha ( t _ { k + 1 } - t _ { k } ) } \Psi _ { k } \leq ( \tau _ { M } - \eta _ { m } ) e ^ { 2 \alpha [ \tau _ { M } + ( N - 2 ) ( \tau _ { M } - \eta _ { m } ) ] } \\ & \quad \times \left [ \sum _ { j = 0 } ^ { N - 2 } \int _ { s _ { k } - j - 1 } ^ { s _ { k + 1 } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } | \sqrt { Q _ { i ^ { * } _ { k - j } } } C _ { i ^ { * } _ { k - j } } \dot { x } ( s ) | ^ { 2 } d s \\ & \quad \ \ - \sum _ { i = 1 } ^ { N } ( N - 1 ) \int _ { s _ { k } } ^ { s _ { k + 1 } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } | \sqrt { Q _ { i ^ { * } _ { k - j } } } C _ { i ^ { * } _ { k - j } } \dot { x } ( s ) | ^ { 2 } d s \\ & \quad \ - \sum _ { j = 0 } ^ { N - 2 } ( N - 1 - j ) \int _ { s _ { k - j - 1 } } ^ { s _ { k } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } \times | \sqrt { Q _ { i ^ { * } _ { k - j } } } C _ { i ^ { * } _ { k - j } } \dot { x } ( s ) | ^ { 2 } d s \right ] \\ & \quad \times \left [ \sum _ { j = 0 } ^ { N - 2 } ( N - 2 - j ) \int _ { s _ { k - j - 1 } } ^ { s _ { k + 1 } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } \times | \sqrt { Q _ { i ^ { * } _ { k - j } } } C _ { i ^ { * } _ { k - j } } \dot { x } ( s ) | ^ { 2 } d s \\ & \quad \ + ( N - 1 ) \int _ { s _ { k } } ^ { s _ { k + 1 } } e ^ { 2 \alpha ( s - t _ { k + 1 } ) } | \sqrt { Q _ { i ^ { * } _ { k + 1 } } } C _ { i ^ { * } _ { k + 1 } } \dot { x } ( s ) | ^ { 2 } d s \right ] \\ & = \Psi _ { k + 1 }, \\ & \text{which implies (36). Hence, (36) and (40) yield (38). } \ \text{m}$$

which implies (36). Hence, (36) and (40) yield (38).

## References

- [1] P. Antsaklis and J. Baillieul. Special issue on technology of networked control systems. Proceedings of the IEEE , 95(1), 5-8, 2007.
- [2] N.W. Bauer, P.J.H. Maas and W.P.M.H. Heemels. Stability analysis of networked control systems: A sum of squares approach. Automatica , 48, 1514-1524, 2012.
- [3] M.B.G. Cloosterman, L. Hetel, N. van de Wouw, W.P.M.H. Heemels, J. Daafouz and H. Nijmeijer. Controller synthesis for networked control systems. Automatica , 46(10), 1584-1594, 2010.
- [4] D. Christmann, R. Gotzhein, S. Siegmund and F. Wirth. Realization of Try-Once-Discard in Wireless Multi-hop Networks. In Proceedings of the 18th World IFAC Congress , Milano, Italy, 2011.
- [5] D. Dacic and D. Nesic. Quadratic stabilization of linear networked control systems via simultaneous protocol and controller design. Automatica , 43, 1145-1155, 2007.

| [6]   | M.C.F. Donkers, W.P.M.H. Heemels, N. van de Wouw and L. Hetel. Sta- bility analysis of networked control systems using a switched linear sys- tems approach. IEEE Transactions on Automatic Control , 56(9), 2101- 2115, 2011.                             |
|-------|------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [7]   | E. Fridman. Use of models with aftereffect in the problem of design of op- timal digital control. Automation and Remote Control , 53(10), 1523-1528, 1992.                                                                                                 |
| [8]   | E. Fridman. A refined input delay approach to sampled-data control. Au- tomatica , 46, 421-427, 2010.                                                                                                                                                      |
| [9]   | E. Fridman, A. Seuret and J.P. Richard. Robust sampled-data stabilization of linear systems: an input delay approach. Automatica , 40, 1441-1446, 2004.                                                                                                    |
| [10]  | H. Fujioka. A discrete-time approach to stability analysis of systems with aperiodic sample-and-hold devices. IEEE Transactions on Automatic Con- trol , 54(10), 2440-2445, 2009.                                                                          |
| [11]  | H. Gao, T. Chen and T. Chai. Passivity and passification for networked control systems. SIAM Journal on Control and Optimization , 46(4): 1299- 1322, 2007.                                                                                                |
| [12]  | H. Gao, X. Meng, T. Chen and J. Lam. Stabilization of networked control systems via dynamic output-feedback controllers. SIAM Journal on Control and Optimization , 48(5): 3643-3658, 2010.                                                                |
| [13]  | J.C. Geromel, R.H. Korogui and J. Bernussou. H 2 and H ∞ robust out- put feedback control for continuous time polytopic systems. IET Control Theory & Applications , 1(5), 1541-1549, 2007.                                                                |
| [14]  | K. Gu, V. Kharitonov and J. Chen. Stability of time-delay systems . Birkhauser: Boston, 2003.                                                                                                                                                              |
| [15]  | W.P.M.H. Heemels, A.R. Teel, N. van de Wouw and D. Nesic. Networked control systems with communication constraints: tradeoffs between trans- mission intervals, delays and performance. IEEE Transactions on Auto- matic Control , 55(8), 1781-1796, 2010. |
| [16]  | K. Liu, E. Fridman and L. Hetel. Stability and L 2 -gain analysis of net- worked control systems under Round-Robin scheduling: a time-delay ap- proach. Systems & Control Letters , 61(5): 666-675, 2012.                                                  |
| [17]  | K. Liu, E. Fridman and L. Hetel. Network-based control via a novel analysis of hybrid systems with time-varying delays. In Proceedings of the 51th IEEE Conference on Decision and Control , Maui, Hawaii, 2012.                                           |
| [18]  | P. Naghshtabrizi, J. Hespanha and A.R. Teel. Exponential stability of im- pulsive systems with application to uncertain sampled-data systems. Sys- tems & Control Letters , 57, 378-385, 2008.                                                             |

- [19] P. Naghshtabrizi, J. Hespanha and A.R. Teel. Stability of delay impulsive systems with application to networked control systems. Trans. of the Inst. of Measurement and Control, Special Issue on Hybrid and Switched Systems , 32(5), 511-528, 2010.
- [20] D. Nesic and A.R. Teel. Input-output stability properties of networked control systems. IEEE Transactions on Automatic Control , 49(10), 16501667, 2004.
- [21] P. Park, J. Ko and C. Jeong. Reciprocally convex approach to stability of systems with time-varying delays. Automatica , 47, 235-238, 2011.
- [22] A. Seuret. A novel stability analysis of linear systems under asynchronous samplings. Automatica , 48, 177-182, 2012.
- [23] G.C. Walsh, H. Ye and L.G. Bushnell. Stability analysis of networked control systems. IEEE Control Systems Technology , 10, 438-446, 2002.
Kun Liu Emilia Fridman Laurentiu Hetel
Networked control systems in the presence of scheduling protocols and communication delays
This paper develops the time-delay approach to Networked Control Systems
(NCSs) in the presence of variable transmission delays, sampling intervals and
communication constraints. The system sensor nodes are supposed to be
distributed over a network. Due to communication constraints only one node
output is transmitted through the communication channel at once. The scheduling
of sensor information towards the controller is ruled by a weighted
Try-Once-Discard (TOD) or by Round-Robin (RR) protocols. Differently from the
existing results on NCSs in the presence of scheduling protocols (in the
frameworks of hybrid and discrete-time systems), we allow the communication
delays to be greater than the sampling intervals. A novel hybrid system model
for the closed-loop system is presented that contains {\it time-varying delays
in the continuous dynamics and in the reset conditions}. A new
Lyapunov-Krasovskii method, which is based on discontinuous in time Lyapunov
functionals is introduced for the stability analysis of the delayed hybrid
systems. Polytopic type uncertainties in the system model can be easily
included in the analysis. The efficiency of the time-delay approach is
illustrated on the examples of uncertain cart-pendulum and of batch reactor.