Exponential stability of a PI plus reset integrator controller
by a sampled-data system approach
M. A. Davóa,∗, F. Gouaisbautb,c , A. Bañosd , S. Tarbouriechb,e , A. Seuretb,e

arXiv:1603.02458v1 [] 8 Mar 2016

a Gipsa-lab,

Grenoble Campus, 11 rue des mathématiques, BP 46, 38402 Saint Martin d’Héres Cedex,
France.
b CNRS, LAAS, 7 avenue du Colonel Roche, F-31400 Toulouse, France.
c Univ de Toulouse, UPS, LAAS, F-31400, Toulouse, France.
d Dpto. Informática y Sistemas, University of Murcia, 30071 Murcia, Spain.
e Univ de Toulouse, LAAS, F-31400, Toulouse, France.

Abstract
The paper deals with the stability analysis of time-delay reset control systems, for which
the resetting law is assumed to satisfy a time-dependent condition. A stability analysis
of the closed-loop system is performed based on an appropriate sampled-data system.
New linear matrix inequality (LMI) conditions are proposed to ensure the exponential
stability of the closed-loop resulting of the connection of a plant with a proportional and
integral controller together with a reset integrator (PI+RI).
Keywords: Time-delay; reset control systems; sampled-data systems; stability analysis;
LMI.

1. Introduction
Reset control systems represent a particular class of hybrid dynamical systems [25],
in which the system state (or part of it) is reset at the instants it intersects some reset
surface. When the reset is defined as a function of time, one can consider the reset control
system as an impulsive system, in which some states of the controller are subjected to
impulsive actions (they are usually set to zero) at some instants. The research field
of reset control systems is rather focused on feedback control with the main aim of
overcoming the performance of linear control systems while keeping low the order of
the system and its complexity. Since the seminal work of J. C. Clegg in 1958 ([12]),
a multitude of works have brought to light the benefits of the reset compensation over
linear feedback control schemes (e.g., [33, 30, 10, 40, 20, 3, 27, 46, 18, 19, 29]). Despite
the significant progress on reset control systems, only few works have considered reset
I c 2016.
This manuscript version is made available under the CC-BY-NC-ND 4.0 license
http://creativecommons.org/licenses/by-nc-nd/4.0/
∗ Corresponding author
Email addresses: Miguel.Davo-Navarro@gipsa-lab.fr (M. A. Davó), fgouaisb@laas.fr (F.
Gouaisbaut), abanos@um.es (A. Baños), sophie.tarbouriech@laas.fr (S. Tarbouriech),
aseuret@laas.fr (A. Seuret)

control strategies for processes in presence of time-delay ([17, 8, 38, 6, 7]). In particular,
it has been shown that a proportional and integral (PI) controller may be improved
without increasing the cost of feedback (sensitivity to sensor noise) by simply adding in
parallel a Clegg integrator (CI) [33] resulting in a reset controller known as PI+CI [8].
The benefits of this controller highlight specially in integrating processes [37]. Although
the PI+CI was originally designed to perform reset actions when the error signal is zero,
several modification of the resetting law have been propounded in the literature with the
aim of improving the performance of the closed-loop system (see, e.g., [3]): fixed reset
band, variable reset band, etc. In order to encompass the most general resetting law,
this article concerns with a generalization of the PI+CI controller in which the Clegg
integrator is replaced with a general reset integrator. This controller will be referred to
as PI+RI (RI stands for reset integrator), and its reset actions can be specified by an
arbitrary time/state-dependent resetting law, instead of only the zero crossing resetting
law (reset actions are performed when the input is zero).
Although, several real-world applications (see for instance [45, 14, 16, 5]) have shown
the practical advantages of the PI+CI controller (particular case of the PI+RI for a kind
of resetting laws), there does not exist results studying the stability of time-delay reset
control systems comprising this controller. Recently, different results on the stability of
time-delay reset control systems have been developed: delay-independent criteria [2, 7],
delay-dependent criteria [9, 41, 13], uncertain reset surface [28], input-to-output stability
[35, 36]. In general, these stability results are based on the Lyapunov-Krasovskii (LK)
theorem and functionals, and involve two conditions; the first one (flow condition) imposes constraints over the derivative of a functional between the reset instants and the
other condition (jump condition) deals with the instantaneous changes of the functional
due to the reset actions. The main limitation of these results comes from the flow condition, since it requires the stability of the base system (system without reset actions).
The challenge in the stability analysis of the PI+RI controller arises from its structure of
two integrators in parallel, that leads into a non-minimal realization of its base system,
and makes more difficult to fulfill the flow condition.
This work focuses on providing exponential stability criteria of a time-delay reset
control system whose feedback controller is a PI+RI. The main idea of this work is the
formulation of the reset integrator into the framework of sampled-data systems ([22]),
which allows representing the time-delay reset control system as a sampled-data system.
In this way, exponential stability of the reset control system can be concluded from
the analysis of the resulting sampled-data system. The advantages of this approach are
twofold: a minimal realization of the PI+RI controller is obtained, and the solutions of
the sampled-data system are continuous which leads into a continuous LK functional, i.e.,
no jump condition is required. The stability analysis of sampled-data systems has been
performed in many works using approaches related to robust analysis [24], a time-delay
modelling [23], an impulsive system interpretation [21, 39], the use of looped-functionals
[43] or a clock-dependent Lyapunov approach [11]. However, these methods cannot be
directly applied into the sampled-data system under study. This is due to its particular
structure in which the plant is also a closed-loop system with an internal time-delay.
The stability criterion provided in this work is based on the extension of the results of
[42, 43], and its combination with a new LK functional provided in [44]. Furthermore,
this work can be viewed as a comprehensive version of the conference work [15], where
only the asymptotic stability of the closed-loop system has been ensured. The paper can
2

be also complementary to [34] due to the type of LK functionals to be considered and
the temporal nature of the proposed conditions.
This article is organized as follows. Section 2 formulates the problem and establishes
an equivalence between reset control and sampled-data systems. Section 3 presents several propositions on exponential stability of time-delay reset control systems by means
of a sampled-data system formulation. Some examples of applications of the proposed
criterion are given in Section 4. Finally, Section 5 concludes the article.
Notation: Throughout the article, the sets N, R, Rn , Rn×n and Sn denote the sets of
natural numbers, real numbers, the n-dimensional Euclidean space and the set of n × n
matrices and symmetric matrices of dimension n × n, respectively. A column vector is
denoted by x ∈ Rn and x> , its transpose. Given two vectors x1 and x2 , we write (x1 , x2 )
> >
n
n
to denote [x>
1 x2 ] . The notation kxk is the Euclidean norm for x ∈ R . C([a, b], R )
n
stands for the set of continuous functions mapping [a, b] to R . The identity matrix
and the zero matrix of adequate dimensions are denoted by I and 0, respectively. The
notation P > 0 for P ∈ Sn means that P is positive definite (P < 0 means negative
definite). The set of symmetric positive definite matrices is denoted by Sn+ . For a matrix
A ∈ Rn×n , the notation He(A) refers to A + A> . For any matrices A, B, the notation
diag(A, B) stands for the block diagonal matrix [ A0 B0 ].
2. Problem statement
The system under study is the standard feedback control system, which is formed by
a plant P and a PI+RI reset controller. The plant is considered to be a linear and time
invariant (LTI) system without direct feedthrough from the input to the output. The
minimal state-space description of the plant is given by

ẋp (t) = Ap xp (t) + Bp up (t),
P :
(1)
yp (t) = Cp xp (t),
where xp ∈ Rnp and Ap ∈ Rnp ×np , Bp ∈ Rnp ×1 , Cp ∈ R1×np .
Originally, the Clegg integrator is a nonlinear integrator whose state is reset to zero
when its input is zero [12]. However, the term Clegg integrator has been also used to
refer to nonlinear integrators with more general state-dependent resetting laws. In this
work, we consider the term reset integrator (RI) to refer to a nonlinear integrator which
is reset to zero at those instants specified by an arbitrary time/state-dependent resetting
law. The parallel connection of a PI controller and a reset integrator (RI) leads into
the PI+RI controller (similarly to the PI+CI controller). In addition, it is assumed that
for a particular resetting law there exists for every execution of the system an strictly
increasing sequence of nonzero instants defined by
T = {tk }k∈N ,
for which the following assumptions stand
t0 = 0,

Tk := tk+1 − tk ∈ [Tm TM ],

lim tk = +∞,

k→+∞

(2)

where 0 ≤ Tm ≤ TM are positive numbers. The last equation means that T has no
accumulation points and that Zeno phenomenon is avoided. A state-space realization
3

of the PI+RI is obtained by using a state xr = (xi , xri ) ∈ R2 , where xi is the integral
term state and xri corresponds to the reset integrator state. As a result, the PI+RI is
described by the following impulsive differential equation

(t, xr (t), ur (t)) ∈
/ R,
 ẋr (t) = Br ur (t),
xr (t+ ) = Aρ xr (t),
(t, xr (t), ur (t)) ∈ R,
P I + RI :
(3)

yr (t) = Cr xr (t) + kp ur (t),
for some resetting set R and where Br = ki [ 11 ], Cr = [ 1 − pr pr ], Aρ = [ 10 00 ], and
xr (t+ ) = lim xr (t + )
→0
>0

is referred to as after-reset state. Besides the tuning parameters of the PI controller (kp
and ki ), the PI+RI possesses an extra dimensionless parameter pr ∈ [0, 1], referred to
as the reset ratio (as the PI+CI controller). This parameter represents the part of the
whole integral term over which the reset action is applied.
Note that if the resetting set R does not guaranteed the assumption on T, then a
time regularization technique (see, e.g., [32]) can be applied by augmenting the PI+RI
with an extra state.
The feedback interconnection of the plant P and the PI+RI controller is considered
to be affected by a time-delay h > 0. The connection between the plant and the PI+RI
controller is defined by up (t) = yr (t) and ur (t) = −yp (t − h). The selection of the system
state x(t) = (xp (t), xr (t)) ∈ Rn , with n = np + 2, leads into the following description of
the closed-loop system

/ S,
 ẋ(t) = Ax(t) + Ad x(t − h), (t, x(t)) ∈
x(t+ ) = AR x(t),
(t, x(t)) ∈ S,
(4)

x(t) = φ(t),
t ∈ [−h, 0],
where φ ∈ C([−h, 0], Rn ) is the initial condition function, and the matrices A, Ad , and
AR are given by






Ap Bp Cr
−kp Bp Cp 0
I 0
A=
, Ad =
, AR =
.
(5)
0
0
−Br Cp
0
0 Aρ
According to the resetting set of the reset controller PI+RI, the resetting set of the
system (4) is defined by
S = {(t, x) ∈ R × Rn : (t, xr ) ∈ R}.

(6)

Due to the assumption on the resetting sequence, it is possible to conclude from the
results in [1] that there exists a unique solution x(t, φ), or simply x(t), for t ∈ [−h, ∞).
2.1. Reset integrator as a sampled-data system
In this section, we show how the reset integrator can be reformulated into the framework of sampled-data systems. Consider a RI represented by the following equation

/ Rci
 ẋri (t) = uri (t), (t, xri (t), uri (t)) ∈
xri (t+ ) = 0,
(t, xri (t), uri (t)) ∈ Rci ,
(7)

yri (t) = xri (t),
4

whose reset instants are defined by the resetting set Rci . Again, it is assumed that
for every execution of the system there exists an strictly increasing sequence of nonzero
instants defined by T = {tk }k∈N for which (2) holds (note that a more general resetting
set can be considered as long as the assumption holds). Then, from the input/output
point of view, the above system is equivalent to the following sampled-data system

ẋs (t) = us (t),
(8)
ys (t) = xs (t) − xs (tk ), t ∈ (tk , tk+1 ],
where tk ∈ T. The output of both systems are the same for any initial condition xri (0) =
xs (0) and all input signal e. Therefore, a reset control system whose controller comprises
a RI (7) can be transformed into a sampled-data system by exchanging the RI with (8).
In this way, the evolution of both systems are equivalent except for the states xi of (7)
and xs of (8). However, it should be pointed out that yri (t) = 0 implies xri (t) = 0, but
the counterpart implication for (8) is false.
As a consequence, the internal stability of the reset control system does not imply
the internal stability of the sampled-data system. However, we will use the reserve
assumption, which is correct, to guarantee the asymptotic stability of the reset control
system as a consequence of the asymptotic stability of the sampled-data system.
2.2. PI+RI as a sampled-data system
Let first us focus on the sampled-data formulation of the PI+RI controller. This can
be done by simply considering the formulation (8) into the PI+RI, which results into


 ẋi (t)
= Br ur (t),
ẋs (t)
P I + RI :
(9)

yr (t) = (1−pr )xi (t) + pr (xs (t)−xs (tk )) + kp ur (t), t ∈ (tk , tk+1 ].
A simple inspection of the above equation shows that the evolution of both xi and xs
are equivalent, and thus, the realization of the system is non-minimal. However, a minimal realization of the PI+RI, equivalent (regarding input/output behavior) to (3), can
be obtained by removing xi and carefully considering the initial condition (xi (0), xri (0))
for (3). As a result, the PI+RI controller can be represented as follows:

 ẋs (t) = ki ur (t),
yr (t) = xs (t) − pr (xi (0) − xri (0)) + kp ur (t), t ∈ [0, t1 ],
(10)
P I + RI :

yr (t) = xs (t) − pr xs (tk ) + kp ur (t),
t ∈ (tk , tk+1 ].
Note that the equation yr (t) = xs (t) − pr (xi (0) − xri (0)) + kp ur (t) guarantees the
equivalence between the two realizations, (3) and (10), of the PI+RI controller.
Finally, considering the realization (10), the reset control system (4) is transformed
into the following sampled-data system

 Ẋ(t) = ΛX(t) + Λd X(t − h) + Λu(t),
(11)
u(t) = KX(tk ),
t ∈ (tk , tk+1 ],

X(t) = φX (t),
t ∈ [−h, 0],
5

where X(t) = (xp (t), xs (t)), φX the initial condition and the matrices Λ, Λd , and K
given by






Ap B p
−kp Bp Cp 0
0
0
Λ=
, Λd =
, K=
.
(12)
0
0
−ki Cp
0
0 −pr
Note that, if φX (t) = [ Inp +1 0 ]φ(t), t ∈ [−h, 0] and xi (0) − xri (0) = [ 0np 1 ]u(0),
then the evolution of the plant state, xp , is the same in both formulations and the
states xi and xri of (3) are directly retrieved from X(t) since xi (t) = xs (t) and xri (t) =
xs (t) + [ 0np 1 ]u(t).
3. Exponential stability analysis
3.1. Preliminaries
This section provides an exponential stability criterion for the time-delay reset control system (4). First, it is proved the equivalence between the systems (4) and (11),
regarding the exponential stability. Second, we propose sufficient conditions for the exponential stability of the sample-data system (11), allowing to conclude on the stability
of the system (4). As mentioned before, due to the particular structure of (11), the
results using the framework of sampled-data systems are not applicable.
The stability criterion developed in this section is based on a LK functional that
depends on both the solution X(t) and its derivative Ẋ(t). Therefore, it is assumed
that the initial conditions are restricted to the space of absolutely continuous functions
φX : [−h, 0] → Rn with square integrable derivative (see, e.g., [22]), which is denoted by
W [−h, 0] and its associated norm is
Z 0
 21
2
.
(13)
kφX kW = max kφX (θ)k +
kφ̇X (θ)k dθ
θ∈[−h,0]

−h

As discussed in [22], the stability results corresponding to continuous initial conditions
are equivalent to the case of absolutely continuous initial conditions.
Definition 1. The trivial solution of system (11) is said to be exponentially stable with
decay rate α > 0, if for every  > 0, there exists δ = δ() > 0 such that
kφX kW < δ

⇒

kX(t)k < e−αt ,

∀t ≥ 0.

(14)

The exponential stability of the time-delay reset control system (4) is defined equivalently with kφkW < δ and kx(t)k < e−αt , t ≥ 0.
Proposition 1. If the sampled-data system (11) is exponentially stable then the reset
control system (4) is exponentially stable.
Proof. The proof follows straightforwardly from the following relations:
φX (t)

=

[ Inp +1 0 ]φ(t),

xi (0) − xri (0)

=

[ 0np 1 ]u(0),

X(t)

=

[ Inp +1 0 ]x(t),

xri (t)

=

xs (t) − [ 0np 1 ]u(t),
6

t ∈ [−h, 0],
t ≥ 0,

t ≥ 0.

(15)

The next result provides sufficient conditions for the exponential stability of the
sampled-data system (11), which are based on the existence of a LK functional and a
family of functionals Vk . In these paper, the analysis is acheived thanks to the definition
of a lifted function, χk ∈ C([0, Tk ] × [−h, 0], Rn ), which is defined, for any reset instant
tk ∈ T as follows:
χk (τ, θ) = χ(tk + τ + θ).
(16)
Using this notation, the dynamics of system (11) is rewritten as
χ̇k (τ, 0) = Λχk (τ, 0) + Λd χk (τ, −h) + ΛKχk (0, 0),

∀τ ∈ [0, Tk ] ∀k ∈ N.

(17)

The following proposition is an adaptation of the technique of looped functionals used
in [42, 43].
Proposition 2. For given scalars 0 < α and 0 ≤ Tm ≤ TM , consider a functional
V : W [−h, 0] → R for which there exists some positive scalar µ1 , µ2 such that, for all
φ ∈ W [−h, 0], the following relations hold
µ1 kφ(0)k2 ≤ V (φ) ≤ µ2 kφk2W .

(18)

Moreover, if there exist a functional Vk : [0, Tk ] × C([0, Tk ] × [−h, 0], Rn ) → R, for
any k in N verifying, first, the looping condition
e2αTk Vk (Tk , χk ) = Vk (0, χk ),

(19)

Vk (0, χk ) ≤ ηV (χk (0, ·)),

(20)

−ηV (χk (0, ·)) ≤ e2ατ Vk (τ, χk )

(21)

and

for some η ≥ 0 and τ ∈ [0, Tk ], and, second,
Ẇ (σ, χk ) =


d
e2ασ (V (χk (σ, ·)) + Vk (σ, χk )) ≤ 0
dσ

(22)

for all σ ∈ [0, Tk ], then system (11) is exponentially stable with decay rate α > 0 for any
reset instants satisfying (2).
Proof. Consider the k-ith reset instant, then the integration of (22) with respect to σ
over the interval [0, τ ] lead into
e2ατ (V (χk (τ, ·)) + Vk (τ, χk )) − V (χk (0, ·)) − Vk (0, χk ) ≤ 0.

(23)

Applying inequalities (19), (20), and (21), it follows

and

V (χk (τ, ·)) ≤ (1 + 2η)e−2ατ V (χk (0, ·)),

(24)

V (χk (Tk , ·)) ≤ e−2αTk V (χk (0, ·)).

(25)

7

Considering the previous reset instants and the Lipschitz continuity of the timederivative of the state, there exists K > (1 + 2η) such that
V (χk (τ, ·)) ≤ Ke−2αt V (φχ ),

(26)

where t = tk + τ . Inequality (18) yield to
µ1 kχk (τ, 0)k2 ≤ V (χk (τ, ·)) ≤ Ke−2αt V (φχ ) ≤ Ke−2αt µ2 kφχ k2W .
(27)
q
µ1
Finally, for a given  > 0, set δ <
Kµ2 , and thus, system (11) is exponentially
stable with decay rate α.
3.2. Main result
The next result provides a sufficient criterion for the exponential stability of the
sampled-data system (11) expressed in terms of LMI. The following theorem is stated.
Proposition 3. For given scalars 0 < α and 0 ≤ Tm ≤ TM , the system (11) is exponentially stable with decay rate α if there exist a matrix PN ∈ S(N +1)n , matrices S, R, Qi ,
Xi , Ui ∈ Sn+ , matrices Zi ∈ Rn and Yi ∈ R(N +3)n×n such that the LMIs
PN +



e−2αh
diag(0n , S, 3S, . . . , (2N − 1)S) > 0,
h

(28)

h1 (0)Π1,i + h2 (0, Ti−1 )Π2,i + h4 (0, Ti )N2> Xi N2 < 0,

(29)

h1 (0)Π1,i + h2 (0, Ti )Π2,i + h4 (0, Ti )N2> Xi N2 < 0,

(30)


h1 (Ti−1 )Π1,i + h4 (Ti−1 , Ti )N2> Xi N2 h3 (Ti−1 , Ti )Yi
< 0,
?
−h3 (Ti−1 , Ti )Ui


h1 (Ti )Π1,i + h4 (Ti , Ti )N2> Xi N2 h3 (Ti , Ti )Yi
<0
?
−h3 (Ti , Ti )Ui

(31)
(32)

hold for all i ∈ {1, 2, . . . , M }, where the functions hi , i = 1, 2, 3, 4 are given by
h1 (τ )
h3 (τ, T )

= e2ατ ,

= e2αT

e2ατ −1
2α



,

h2 (τ, T )

=

h4 (τ, T )

=

e2αT −e2ατ
,
2α
 αT
2

e −4 e


+ 1 e2ατ + 2e2αT ,
(33)

and where the matrices Πj,i , i = 1, . . . m, j = 1, 2 are given by

Π1,i = He G> P H + ΣN + h2 F > RF − Γ>
N RN ΓN

>
>
− N12 Qi N12 − He N2 Zi N2 + He (Yi N12 )
>

>



>

Π2,i = F Ui F + He F Qi N12 + He F Zi N2

8



(34)





with N1 = In 0n,n(N+1) , N2 = 0n,n(N+1) In , N12 = N1 −N2 and where
e−2αh diag(R, 3R, . . . , (2N +1)R),


In
0n
0n,nN 0n
ΣN = diag(S, −e−2αh S, 0(N +1)n ), G =
,
0nN,n 0nN,n hInN 0n


F = [Λ Λd 0n,nN ΛK] ,
Γ(k) = I (−I)k+1 γk0 I . . . γkN −1 I 0n ,




F
Γ(0)



 Γ(1) 
 Γ(0) 
−(2i + 1)(1 − (−1)k+i ), i ≤ k,


i


, γk =
, ΓN  .
H=


..
0,
i > k.
 ..



.
Γ(N )
Γ(N−1)

RN =

(35)

Proof. The proof of this proposition is based on the application of Proposition 2. Therefore, for the sake of clarity, the proof is divided into three steps to show how each condition
of 2 can be verified. These steps are described as follows:
1) Construction of the functional V and an LMI condition ensuring the satisfaction of
inequality (18);
2) Construction of the looped-functional Vk , satisfying the looping condition (19) and
conditions (24) and (25);
3) Design of LMI conditions guaranteeing (22).
1) Contruction of V : In the literature, several class of functionals have been studied, leading to a wide class of stability conditions. In this paper, we focus on a recent
development on the stability analysis of time-delay systems based on Bessel inequality
and Legendre’s polynomials provided in [44]. Therein, the following LK functionals has
been introduced
Z 0
V (χk (τ, ·)) = χ̃>
(τ
)P
χ̃
(τ
)
+
e2αs χ>
N k
k
k (τ, s)Sχk (τ, s)ds
−h
(36)
Z 0 Z 0
2αs >
+h
e χ̇k (τ, s)Rχ̇k (τ, s)dsdβ,
−h

β

where the augmented state χ̃k (τ ) is given by

χk (τ, 0)
 R 0 L (s)χ (τ, s)ds
k

−h 0
χ̃k (τ ) = 
..

.
 R
0
L
(s)χ
k (τ, s)ds
−h N −1




,



(37)

where we use the Legendre polynomials considered over the interval [−h, 0] and given by
k
X
Lk = (−1)
(−1)l
k

l=0

k
l



k+l
l

9





u+h
h

l
,

∀k ∈ N,

(38)

where the notation ( ab ) with a, b in N denotes the binomial coefficients. The potential of
using of these polynomials has been demonstrated in [44] and arisen from their orthogonal
properties, the simple evaluation of the boundary values and the recursive relation to
compute their derivative using the same polynomial basis of lower degree. All these
details are already provided in [44] and are therefore omitted in this paper.
In order to prove that condition (18) holds for this Lyapunov functionals is inspired
from [44]. Firstly, considering the bounds e−2αh ≤ e2αs ≤ 1 for s ∈ [−h, 0], then, since
R > 0, we have
R0
V (χk (τ, ·)) ≤ χ̃>
+ e2αh −h χ>
k (τ )P
k (τ, s)Sχk (τ, s)ds

N χ̃k (τ ) −2αh
(39)
e
>
≤ χ̃k (τ ) PN + h diag(0n , S, 3S, . . . , (2N − 1)S) χ̃k (τ ).
The last inequality has been derived accordingly to [44] with the use of the BesselLegendre inequality. Then, if the condition (28) holds, it directly follows that the LK
functional (36) satisfies condition (18).
2) Contruction of Vk : Consider now the following functionals

Vk (τ, χk ) = f1 (τ, Tk ) ζk> (τ )Q(Tk )ζk (τ ) + He(ζk> (τ )Z(Tk )χk (0, 0))
Z τ
+ f1 (τ, Tk )
χ̇>
k (s, 0)U (Tk )χ̇k (s, 0)ds
+

0
f2 (τ, Tk )χ>
k (0, 0)X(Tk )χk (0, 0)

(40)

+ γe−2ατ kχk (0, 0)k2

for some γ > 0, and where ζk (τ ) = χk (τ, 0) − χk (0, 0), i ∈ {1, . . . , M } satisfies Tk ∈
[Ti−1 , Ti ] and
e2α(Tk −τ ) − 1
f1 (τ, Tk ) =
,
(41)
2α
f2 (τ, Tk ) =
and where



2 
e−2ατ  2αTk
e
− 1 e2ατ − 1 − e2ατ − 1
,
2α

Q(Tk ) =
Z(Tk ) =
U (Tk ) =
X(Tk ) =

Qi ,
Zi ,
Ui ,
Xi ,

if
if
if
if

Tk
Tk
Tk
Tk

∈ [Ti−1 ,
∈ [Ti−1 ,
∈ [Ti−1 ,
∈ [Ti−1 ,

Ti ],
Ti ],
Ti ],
Ti ].

(42)

(43)

This selection of functionals represents a new discretization method with respect to the
literature. The parameters of the looped-functional are piecewise constants and have
jumps only at the sampling instants. However, the due to the looping conditions, the
effects of these jumps vanishes and the functional Vk is still continuous.
The functionals Vk satisfy condition (19) for all k > 0 since e2αTk Vk (Tk , χk ) =
Vk (0, χk ) = γkχk (0, 0)k2 . Now consider the following bound of the LK functional
µ1 kχk (τ, 0)k2 ≤ V (χk (τ, ·)),

(44)

then by taking τ = 0, it can be obtained the following inequality
Vk (0, χk ) ≤

γ
V (χk (0, ·)).
µ1
10

(45)

Therefore, the functionals Vk and V satisfy condition (20) for η ≥ µγ1 .
Consider now the following term of the functionals Vk

f1 (τ, Tk ) ζk> (τ )Qi ζk (τ ) + He(ζk> (τ )Zi χk (0, 0)) + γe−2ατ kχk (0, 0)k2 .

(46)

The above term is nonnegative if the following LMI holds for τ ∈ [0, Tk ], all k > 0 and a
sufficiently large γ


f1 (τ, Tk )Qi f1 (τ, Tk )Zi
≥ 0.
(47)
?
γe−2ατ In
Note that Qi is positive definite, then applying the Schur complement, the above LMI
holds if
γIn − e2ατ f1 (τ, Tk )Zi> Q−1
(48)
i Zi ≥ 0.

Since e2ατ f1 (τ, Tk ) is bounded for τ ∈ [0, Tk ] and Tk ∈ [Tm , TM ], there exists a sufficiently
large γ such that the LMI (47) holds for all i ∈ {1, . . . , M }. In addition, the nonnegativeness of the term (46) and Ui > 0 imply that γ might be chosen sufficiently large
such that
2
2ατ
2ατ
−γkχ>
f2 (τ, Tk )χ>
Vk (τ, χk ).
k (0, 0)k ≤ e
k (0, 0)Xi χk (0, 0) ≤ e

(49)

Hence, considering equation (44), it follows that condition (21) is satisfied for η ≥ µγ1 .
Note that γ does not play any role in the conditions of Proposition 3, and thus, it
is not necessary to set a value for γ. The previous reasoning proves that there exists a
suitable value for γ, which is enough for the functionals Vk to fulfill the conditions in
Proposition 2.
3) Design of LMI conditions guaranteeing (22): Following [44], the computation
of the time-derivative of e2ατ V (χk (τ, ·)) along the trajectories of the system (11) is given
by


d 2ατ
>
˙
e V (χk (τ, ·)) = e2ατ He(χ̃>
k (τ )PN χ̃k (τ )) + 2αχ̃k (τ )PN χ̃k (τ )
dτ
−2αh >
+ χ>
χk (τ, −h)Sχk (τ, −h)
(50)
k (τ, 0)Sχk (τ, 0) − e
Z 0

+ h2 χ̇>
e2αs χ̇>
k (τ, 0)Rχ̇k (τ, 0) − h
k (τ, s)Rχ̇k (τ, s)ds .
−h

The first step consists in computing the time-derivative of the augmented state χ̃k .
Following [44] and the differentiation properties of the Legendre polynomials, we have
χ̃˙ k = Hξk (τ ),
where the matrix H is given in (35) and where the vector ξk is another augmented vector,
which gathers all the relevant information for the analysis and which is given by


χk (τ, 0)


χk (τ, −h)


 1 R 0 L (s)χ (τ, s)ds 
k

 h −h 0
.
(51)
ξk (τ ) = 
..


.



 1 R0
 h −h LN −1 (s)χk (τ, s)ds 
χk (0, 0)
11

Moreover noting that χ̃k (τ ) = Gξk (τ ), the following expression of the time-derivative
of e2ατ V (χk (τ, ·)) is rewritten as follows:


d 2ατ
e V (χk (τ, ·)) = e2ατ ξk> (τ )Φξk (τ ) − hI(χ̇k (τ, ·))
dτ
with
I(χ̇k (τ, ·)) =

Z

0

−h

e2αs χ̇>
k (τ, s)Rχ̇k (τ, s)ds.

(52)

(53)

The next step consists in expressing the time-derivative of W = e2ατ (V + Vk ). The
calculations lead to

Ẇ(τ, χk ) = h1 (τ ) ξk> (τ )Φξk (τ ) − hI(χ̇k (τ, ·))

− h1 (τ ) ζk> (τ )Qi ζk (τ ) + He(ζk> (τ )Zi χk (0, 0))
Z τ
χ̇>
− h1 (τ )
k (s, 0)Ui χ̇k (s, 0)ds
(54)
0

>
>
+ h2 (τ, Tk ) He(χ̇k (τ, 0)Qi ζk (τ )) + He(χ̇k (τ, 0)Zi χk (0, 0))

>
+ h2 (τ, Tk ) χ̇>
k (τ, 0)Ui χ̇k (τ, 0) + χ̇k (0, 0)Ui χ̇k (0, 0)

+ (e2αTk + 1)e2ατ − 2e4ατ χ>
k (0, 0)Xi χk (0, 0).
In the previous expression, there are two integral quadratic terms that cannot be
included directly into an LMI. Therefore, the integral terms in the right-hand side of the
above equation are bounded by the following two inequalities


N
X
1

(55)
I(χ̇k (τ, ·)) ≥ ξk (τ )>  (2j + 1)Γ>
N (j)RΓN (j) ξk (τ ),
h
j=0
−

Z
0

τ

>
>
−1 >
χ̇>
Y ξk (τ ).
k (s, 0)U χ̇k (s, 0)ds ≤ 2ξk (τ )Y ζk (τ ) + τ ξk (τ )Y U

(56)

The first inequality was proposed in [44] based on the properties of Legendre polynomials (38) and Bessel’s inequality. The second inequality can be found for instance in
[42, 43]. In addition, we can retrieve the following inequalities
τ h1 (τ ) ≤ h3 (τ, T ),
(e

2αTk

+ 1)e

2ατ

− 2e

4ατ

τ ∈ [0, T ],


≤ h4 (τ, Ti ),

Tk ∈ [Ti−1 , Ti ],
τ ∈ [0, Tk ].

(57)

(58)

Using (55), (56), (57), and (58), the time-derivative of the functional W = e2ατ (V +
Vk ) is bounded by Ẇ(τ, χk ) ≤ ξk> (τ )Ψi (τ, Tk )ξk (τ ) with
Ψi (τ, Tk ) = h1 (τ )Π1,i + h2 (τ, Tk )Π2,i + h3 (τ, Ti )Yi Ui−1 Yi> + h4 (τ, σ, Ti )N2> Xi N2 (59)
for some i ∈ {1, . . . , M }. Note that Ψi (τ, Tk ) is linear with respect to the term e2ατ ,
and thus it is necessary and sufficient to ensure the negativity at the edges. This leads
12

to Ψi (0, Tk ) < 0 and Ψi (Tk , Tk ) < 0. In addition, Ψi (τ, Tk ) is also linear with respect to
e2αTk , hence applying the same convexity arguments on Tk , it is obtained
Ψi (0, Ti−1 ) < 0,
Ψi (0, Ti ) < 0,

Ψi (Ti−1 , Ti−1 ) < 0,

(60)

Ψi (Ti , Ti ) < 0,

which are rewritten (using the Schur complement) as (29), (30), (31), and (32), respectively, and the proof is complete.
4. Examples
4.1. Destabilization of a stable base system
A multitude of works have shown the advantages of the reset control to stabilize
unstable systems (only non delayed systems). Consequently, one may wrongly think
that the reset control system must be stable at least whenever its base system is stable.
Several examples in [31, 4] have been proved (for non delayed systems) that intuition
fails, and reset actions can destabilize a stable system. This example shows that reset
actions can destabilize a stable base system, and how the addition of the bounds Tm
and TM lead into an exponentially stable reset control system. Consider a reset control
system that consists of a pure integrator plant, P (s) = 1/s, controlled by the PI+RI (3)
with kp = 1.4, ki = 0.3, pr = 0.5 and a zero crossing resetting law (a PI+CI controller).
The closed-loop matrices defined in (5) are given by






0 0.5 0.5
−1.4 0 0
1 0 0
0  , Ad =  −0.3 0 0  , AR =  0 1 0  .
A= 0 0
(61)
0 0
0
−0.3 0 0
0 0 0
The base system with a minimal realization has been proved to be asymptotically stable
for a time-delay h = 1 by using the results in [26]. Fig. 1 shows how the trajectory
of the base system with initial condition φ(t) = (1, 0), t ∈ [−1, 0] converges to the zero
solution. However, Fig. 2 shows the solution of the reset control system for the same
initial condition, and it can be seen that the trajectory is divergent. Now consider that
the reset instants are forced to satisfy Tm ≤ tk+1 − tk ≤ 1 = TM , then Table 1 shows
the minimum Tm for several values of M , such that Propositions 1 and 3 guarantee the
exponential stability of the reset control system with decay rate α = 1e−6 . It is important
to note that values of N greater than two do not provide further improvements for this
particular example.
4.2. Stabilization of an unstable base system
Consider a reset control system composed by a plant P with matrices


 


0 0
1
Ap =
, Bp =
, Cp = 0 1 ,
1 0.5
1

(62)

and a PI+RI compensator with kp = 1 and ki = 1. The base closed-loop system is not
stable independently of the time-delay. However, the exponential stability of the reset
13

x

1

0

−1
0

5

10

15

20

25
Time

30

35

40

45

50

Figure 1: Trajectory of the base system of the Example 4.1.

10

x

5

0

−5

−10

0

5

10

15

20

25
Time

30

35

40

45

50

Figure 2: Trajectory of the reset control system of the Example 4.1.

14

M

1

3

5

10

50

Tm

0.94

0.80

0.67

0.5

0.41

Table 1: Effect of the parameter M on the lower bound Tm for the Example 4.1. (TM = 1,
N = 2, and α = 10−6 ).
pr = 0.1
pr = 0.5
pr = 0.9

Period (T )

1.5

1

0.5

0

0

5 · 10−2

0.1

0.15

0.2

0.25

0.3

0.35

0.4

Time-delay (h)

Figure 3: Allowable reset period (T ) vs. time-delay (h) for several values of the reset
ratio (pr ) for the Example 4.2. (N = 2, M = 1, and α = 10−6 ).

control system can be guaranteed by the proposed results with a proper bounds Tm and
TM . First, we consider that the reset actions are periodic with period Tm = Tk = TM = T ,
for all k > 0. In Fig. 3, it can be seen the effect of the reset ratio and the time-delay
on the reset period. The results show that the time-delay can increase by increasing the
reset ratio, while the exponential stability is guaranteed, In addition, Fig. 4 shows how
the reset period affects the decay rate of the system. Finally, we consider asynchronous
reset instants with bounds Tm = 0.01TM (> 0). The effect of TM on the decay rate is
shown in Fig. 5.
5. Conclusions
New theoretical conditions through LMI conditions have been proved to ensure the
exponential stability of a class of time-delay reset control systems. The time-delay reset control system is constituted by a PI+RI controller (generalization of the PI+CI
controller) and a time-delay plant. The proposed approach is based on the extension
to results issues from the sampled-data system framework. This work lets some ques15

·10−2

10

Decay rate (α)

8
6
4
2
0

0

0.25

0.5

0.75
Period (T )

1

1.25

Figure 4: Bound of the decay rate (α) vs. reset period (T ) for the Example 4.2. (N = 2
and M = 1).

·10−2

10

Decay rate (α)

8
6
4
2
0

0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

TM
Figure 5: Bound of the decay rate (α) vs. asynchronous reset period (Tm = 0.01TM ) for
the Example 4.2. (N = 2 and M = 1).
16

tions open. In particular it could be interesting to propose conditions in order to design
the controller, by adding for example some performance condition (such as rejection of
perturbations).
Acknowledgements
This work has been supported by Ministerio de Economı́a e Innovación of Spain under project DPI2013-47100-C2-1-P (including FEDER co-funding) and the ANR project
LimICoS contract number 12-BS03-005-01.
References
[1] G. Ballinger and X. Liu. Existence, uniqueness and boundedness results for impulsive delay differential equations. Applicable Analysis: An international Journal, 74(1-2):71–93, 2000.
[2] A. Baños and A. Barreiro. Delay-independent stability of reset systems. IEEE Transactions on
Automatic Control, 54(2):341–346, 2009.
[3] A. Baños and A. Barreiro. Reset control systems. Advances in Industrial Control. Springer, 2012.
[4] A. Baños, J. Carrasco, and A. Barreiro. Reset times-dependent stability of reset control systems.
IEEE Transactions on Automatic Control, 56(1):217–223, 2011.
[5] A. Baños and M. A. Davó. Tuning of reset proportional integral compensation with a variable reset
ratio and reset band. IET Control Theory and Applications, 8(17):1949–1962, 2014.
[6] A. Baños, F. Perez, and J. Carrasco. Network-based reset control systems with time-varying delays.
IEEE Transactions on Industrial Informatics, 10(3):514 – 522, 2014.
[7] A. Baños, F. Perez, S. Tarbouriech, and L. Zaccarian. Low-Complexity Controllers for Time-Delay
Systems: Delay-independent stability via reset loops, chapter 8, pages 111–125. Advances in Delays
and Dynamics. Springer International Publishing, 2014.
[8] A. Baños and A. Vidal. Design of reset control systems: the PI+CI compensator. Journal of
Dynamic Systems, Measurement, and Control, 134(5), 2012.
[9] A. Barreiro and A. Baños. Delay-dependent stability of reset systems. Automatica, 46(1):216–221,
2010.
[10] O. Beker, C. V. Hollot, Y. Chait, and H. Han. Plant with integrator: an example of reset control
overcoming limitations of linear feedback. IEEE Transactions on Automatic Control, 46(11):1797–
1799, 2001.
[11] Corentin Briat. Theoretical and numerical comparisons of looped functionals and clock-dependent
lyapunov functions—the case of periodic and pseudo-periodic systems with impulses. International
Journal of Robust and Nonlinear Control, 2015.
[12] J. C. Clegg. A nonlinear integrator for servomechanisms. Transactions of the American Institute
of Electrical Engineers, Part II: Applications and Industry, 77(1):41–42, 1958.
[13] M. A. Davó and A. Baños. Delay-dependent stability of reset control systems with input/output
delays. In Proceedings of the 52nd Conference on Decision and Control, pages 2018–2023, Florence,
Italy, December 2013.
[14] M. A. Davó and A. Baños. Reset control of a liquid level process. In Proceedings of the 18th
Conference on Emerging Technologies & Factory Automation, pages 1–4, Cagliari, Italy, September
2013.
[15] M. A. Davó, F. Gouaisbaut, A. Baños, S. Tarbouriech, and A. Seuret. Stability of time-delay
reset control systems with time-dependent resetting law. In 5th IFAC Conference on Analysis and
Design of Hybrid Systems (ADHS), Atlanta (USA), October 2015.
[16] E. Delgado, M. D. Cacho, A. Baños, and A. Barreiro. Reset control of synchronous motors with
permanent magnet excitation. In 122nd Mediterranean Conference on Control and Automation,
pages 1347–1352, Palermo, Italy, June 2014.
[17] A. Fernández, A. B. Blas, J. Carrasco, and A. Baños. Reset control for passive bilateral teleoperation. IEEE Transactions on Industrial Electronics, 58(7):3037–3045, 2011.
[18] F. Fichera, C. Prieur, S. Tarbouriech, and L. Zaccarian. Improving the Performance of Linear
Systems by Adding a Hybrid Loop: the Output Feedback Case. In Proceedings of the 2012 American
Control Conference, pages 3192–3197, Montreal, Canada, 2012.

17

[19] F. Fichera, C. Prieur, S. Tarbouriech, and L. Zaccarian. Using Luenberger Observers and Dwelltime Logic for Feedback Hybrid Loops in Continuous-time Control Systems. International Journal
of Robust and Nonlinear Control, 23:1065–1086, 2013.
[20] F. Forni, D. Nešić, and L. Zaccarian. Reset passivation of nolinear controllers via suitable timeregular reset map. Automatica, 47(9):2099–2106, 2011.
[21] E. Fridman. A refined input delay approach to sampled-data control. Automatica, 46(2):421–427,
2010.
[22] E. Fridman. Introduction to Time-Delay Systems. Systems & Control: Foundations & Applications.
Birkhaäuser, 2014.
[23] E. Fridman, A. Seuret, and J.-P. Richard. Robust sampled-data stabilization of linear systems: an
input delay approach. Automatica, 40(8):1441–1446, 2004.
[24] H. Fujioka. Stability analysis of systems with aperiodic sample-and-hold devices. Automatica,
45(3):771–775, 2009.
[25] R. Goebel, R.G. Sanfelice, and A.R. Teel. Hybrid Dynamical Systems: Modeling, Stability and
Robustness. Princeton University Press, 2012.
[26] F. Gouaisbaut and D. Peaucelle. Delay-dependent stability of time-delay systems. In Proceedings of
the 5th IFAC Symposium on Robust Control Design, pages 453–458, Toulouse, France, July 2006.
[27] Y. Guo, Y. Wang, L. Xie, H. Li, and W. Gui. Optimal reset law design and its application to transient response improvement of HDD systems. IEEE Transactions on Control Systems Technology,
19(5):1160–1167, 2011.
[28] Y. Guo and L. Xie. Quadratic stability of reset control systems with delays. In Proceedings of the
10th World Conference on Intelligent Control and Automation, pages 2268–2273, Beijing, China,
July 2012.
[29] L. Hetel, J. Daafouz, S. Tarbouriech, and C. Prieur. Stabilization of linear impulsive systems
through a nearly-periodic reset. Nonlinear Analysis: Hybrid Systems, 7(1):4–15, 2013.
[30] I. M. Horowitz and P. Rosenbaum. Non-linear design for cost of feedback reduction in systems with
large parameter uncertainty. International Journal of Control, 21(6):977–1001, 1975.
[31] H. Hu, Y. Zheng, Y. Chait, and C. V. Hollot. On the zero-input stability of control systems
with Clegg integrators. In Proceedings of the 1997 American Control Conference, volume 1, pages
408–410, June 1997.
[32] K. H. Johansson, M. Egerstedt, J. Lygeros, and S. Sastry. On the regularization of Zeno hybrid
automata. Systems & Control Letters, 38(3):141–150, 1999.
[33] K. R. Kishnan and I. M. Horowitz. Synthesis of a non-linear feedback system with significant
plant-ignorance for prescribed system tolerances. International Journal of Control, 19(4):689–706,
1974.
[34] J. Liu and A.R. Teel. Hybrid dynamical systems with finite memory. In Recent Results on Nonlinear
Delay Control Systems, pages 261–273. Springer International Publishing, 2016.
[35] P. Mercader, J. Carrasco, and A. Baños. IQC analysis for time-delay reset control systems with
first order reset elements. In Proceedings of the 52nd Conference on Decision and Control, pages
2251–2256, Florence, Italy, December 2013.
[36] P. Mercader, M. A. Davó, and A. Baños. H∞ /H2 analysis for time-delay reset control systems. In
Proceedings of the 3rd International Conference on Systems and Control, pages 518–523, Algiers,
Algeria, October 2013.
[37] P. Mercader, M. A. Davó, and A. Baños. Performance analysis of PI and PI+CI compensation for an
IPDT process. In Proceedings of the 23rd Mediterranean Conference on Control and Automation,
Torremolinos, Spain, June 2015.
[38] J. C. Moreno, J. L. Guzmán, J. E. Normey-Rico, A. Baños, and M. Berenguel. A combined FSP and
reset control approach to improve the set-point tracking of dead-time processes. Control Engineering
Practice, 21(4):351–359, 2013.
[39] P. Naghshtabrizi, J.P. Hespanha, and A.R. Teel. Exponential stability of impulsive systems with
application to uncertain sampled-data systems. Systems & Control Letters, 57(5):378–385, 2008.
[40] D. Nešić, L. Zaccarian, and A. R. Teel. Stability properties of reset systems. Automatica, 44(8):2019–
2026, 2008.
[41] J. A. Prieto, A. Barreiro, S. Dormido, and S. Tarbouriech. Delay-dependent stability of reset control
systems with anticipative reset conditions. In Proceedings of the 7th IFAC Symposium on Robust
Control Design, volume 7, pages 219–224, Aalborg, Denmark, June 2012.
[42] A. Seuret. Stability analysis of networked control systems witn asynchronous sampling and input
delay. In Proceedings of the 2011 American Control Conference, pages 533–538, San Francisco,
California, USA, June 2011.

18

[43] A. Seuret. A novel stability analysis of linear systems under asynchronous samplings. Automatica,
48(1):177–182, 2012.
[44] Alexandre Seuret and Frédéric Gouaisbaut. Hierarchy of LMI conditions for the stability analysis
of time delay systems. Systems and Control Letters, 81:1–7, 2015.
[45] A. Vidal and A. Baños. Reset compensation for temperature control: Experimental application on
heat exchangers. Chemical Engineering Journal, 159(1):170–181, 2010.
[46] G. Zhao, D. Nešić, Y. Tan, and J. Wang. Improving L2 gain performance of linear systems by reset
control. In Proceedings of the 19th IFAC World Congress, Cape Town, South Africa, August 2014.

19

