Gu 2003
Gu 2003
Series Editor
William S. Levine
Department of Electrical and Computer Engineering
University of Maryland
College Park, MD 20742-3285
USA
Manfred Morari
ETH
ZUrich, Switzerland
Keqin Gu
Vladimir L. Kharitonov
Jie ehen
Jie Chen
University of California at Riverside
Department of Electrical Engineering
Riverside, CA 92521
USA
TJ213.G83 2003
629.~c21
2003041932
CIP
AII rights reserved. This work may not be translated or copied in whole or in part without the written
permission of the publisher (Springer Science+Business Media, LLC), except for brief excerpts in connec-
tion with reviews or scholarly analysis.Use in connection with any form ofinformation storage and retrieval,
electronic adaptation, computersoftware, or by similar or dissimilar methodology now known or hereafter
developed is forbidden .
The use of general descriptive names, trade names, trademarks , etc., in this publication, even if the former
are not especially identified, is not to be taken as a sign that such names, as understood by the Trade Marks
and Merchandise Marks Act, may accordingly be used freely by anyone.
987 6 5 4 3 2 1
KG dedicates the book to his parents Lijan Gu and Jieping Jiang for
their love and support.
JC dedicates the book to his wife Yan and sons Christopher and
Stephen for their love and trust.
Contents
Preface xiii
Notation xix
Bibliography 327
Index 351
Preface
This book is devoted to the study of stability of time-delay systems, an
area long pursued by mathematicians, physical and life scientists, engineers,
and economists. Given the history and maturity of the discipline, one may
ask what new material could be presented beyond the extensive existing
body of work. The answer lies in the enduring and surprising vitality of the
subject: in particular, over the last decade, research has yielded a number of
notable results. The study of time-delay systems has therefore undergone
a significant conceptual and practical leap. Our aim in this book is to
present some of the highlights of these advances and to further develop the
techniques and tools along recent trends.
By definition, the future evolution of a time-delay system depends not
only on its present state but also on its history. This particular cause and
effect relationship can be most succinctly captured and has been tradi-
tionally modeled by differential-difference equations, or more generally, by
functional differential equations. While in practice many dynamical systems
may be satisfactorily described by ordinary differential equations alone, for
which the system's future evolution depends solely on its current state,
there are times when the delay effect cannot be neglected. In short, time
delay is hardly a matter of rarity; numerous examples presented in this
book and elsewhere attest to its prevalence. And, because of their ubiqui-
tous presence, time-delay systems have been an active area of scientific re-
search in a wide range of natural and social sciences: they are often known
as hereditary systems, systems with after-effect, systems with time-lag,
and more generally, as a subclass of functional differential equations and
infinite-dimensional systems.
The field of time-delay systems had its origins in the 18th century, and
it received substantial attention in the early 20th century in works devoted
to the modeling of biological, ecological, and engineering systems. Stability
of time-delay systems became a formal subject of study in the 1940s, with
contributions from such towering figures as Pontryagin and Bellman. Over
the years its interest and popularity have grown steadily. In particular, in
the last 10 to 15 years there has been a surge in research and a prolifera-
tion of new techniques and results. A cursory glance at the large numbers
of articles published in international conferences, organized workshops, and
archive journals during this period indicates the scale and magnitude of this
progress. Advances in numerical methods and control theory, especially the
theory of robust stability and control, have had substantial impact on the
XIV Preface
field. We observe that the small gain theorem made it possible to formulate
the stability problem of time-delay systems as one of robust stability and
of the structured singular value. Techniques in robust stability analysis of
uncertain polynomials find generalizations to uncertain quasipolynomials,
which serve as models for uncertain time-delay systems. Efficient numeri-
cal algorithms for solving linear matrix inequalities (LMI) have generated
substantial interest to pose stability problems as LMI conditions.
Recognizing these advances, we concluded that a book project would be
the natural outcome of more than a decade's creative work which we hope
will contribute to future studies of the subject. Our primary objective is to
present a sufficiently thorough yet focused treatment of the key methods
and results of the past 10 years or so, and to systematically organize and
discuss in sufficient technical depth many of the results that are otherwise
scattered in various journals and conference proceedings. Since the volume
of such publications can be overwhelming, we hope this book will be a
convenient reference and will provide quick access to important issues and
developments in time-delay systems.
Our particular viewpoint is centered on computability, robust stability
and robust control. We emphasize conceptual significance over technical
details, without sacrificing mathematical rigor. We stress fundamental, in-
tuitive observations behind methods and results, and we try to build links
and relationships among concepts that, at first glance, may seem unre-
lated. In other words, we try to provide an "insider's story" based on our
experience. This perspective has guided us in both the selection and presen-
tation of material: while the facts and results are certainly worthwhile, we
also hope that the ideas and intuition will be helpful in understanding the
nitty-gritty technical details, especially to new readers. On the other hand,
we face the unpleasant yet inevitable task of limiting ourselves to some
selected topics, excluding materials that to some may be equally or even
more important. Examples include classical results, recent developments
on stabilization, control, and filtering of time-delay systems, and those on
stochastic time-delay systems, to name just several. Of course, this selec-
tion merely reflects the authors' viewpoints and preferences rather than
the value of the particular issues. Fortunately, we have been able to refer
readers to several classic and recently-published works on these topics.
Despite widespread interest, it is still uncommon to find time-delay sys-
tems adopted as a regular course at most institutions. For many, the first
and perhaps only encounter with the subject is likely to take place in an
undergraduate course on systems modeling or automatic control. For this
reason, our book has been conceived and written as a monograph aimed
at researchers, practicing engineers, and graduate students. We do not pre-
sume, although we do not preclude the possibility, that it can serve as a
textbook for an advanced graduate course. It is nevertheless appropriate as
a self-study text. Most of the book is easily accessible to a typical second-
Preface xv
Keqin Gu
Southern Illinois University at Edwardsville
Vladimir L. Kharitonov
CINVESTAV-IPN, Mexico
Jie Chen
University of California at Riverside
October 2002
Notation
Sets and Operators
Matrix-related notation
Thanspose, component-wise complex conjugate,
and Hermitian (conjugate) transpose of matrix A.
IAI Matrix formed by taking absolute value of each
entry of A
°
A> 0, A 2:': 0 The matrix A is positive (semi)definite « or ::; similar)
A>B,A2:':B A - B > 0, A - B 2:':
A(A), Ai(A) An eigenvalue and the ith eigenvalue of matrix A
xx Notation
Introduction to Time-Delay
Systems
1.1 Introduction
We begin this chapter by presenting a number of practical examples of time-
delay systems. This should give the reader a glimpse into how widely time
delays may occur in practice. To make the presentation more accessible,
we then analyze a simple time-delay system, which illustrates some of the
important concepts in a concrete manner. This analysis sets the stage for a
formal introduction to functional differential equations as a representation
of time-delay systems, and to such stability analysis tools as the Lyapunov-
Krasovskii Theorem and Razumikhin Theorem. Linear delay systems are
discussed in greater length with frequency domain descriptions. The chap-
ter is concluded with a brief outline of subsequent content.
Ordinary differential equations in the form of
It is known that with a proper design, the delayed resonator can signif-
icantly enhance the system performance, in reducing sensitivity to the
excitation frequency, eliminating the main structure's vibration with the
presence of damping, and simultaneously suppressing two excitation fre-
quencies.D
In this example, the time delay also occurs by way of feedback control
action. However, in contrast to the delay effect in Example 1.2, the delay
at present is intentionally introduced to enhance the system performance.
Time-delay systems may also arise from simplification of some partial
differential equations. Furthermore, time-delay systems are used to model
dynamical systems in many other scientific disciplines including engineer-
ing, biology, ecology, and economics. The Notes section at the end of this
chapter will give a brief historic overview to discuss some related literature.
1.2 A simple time-delay system 5
t = ~, °: ;
values of x(O) and x( -r). Similarly, in order to calculate x at the instant
~ < r, we need x(~) and x(~ - r). For such a ~, since -r ::;
~ -r < 0, x( ~ -r) cannot be generated in the solution process. It is therefore
clear that in order for the solution to be uniquely defined, it is necessary to
specify the value of x(t) for the entire interval -r ::; t ::; 0; in other words,
we need to specify the initial condition
x(t) = ¢(t), t E [-r, 0], (1.8)
with some given continuous function ¢ : [-r, 0] ~ R
Once the initial condition (1.8) is given, the solution is well defined.
Indeed, for t E [0, r], since x(t-r) is already known in this interval, we can
treat (1.7) as an ordinary differential equation, so as to obtain
Continuing this process will allow us to obtain the solution x(t), -r ::; t <
00. Such a process is known as the method of steps.
H(s) = C[h(t)]
exists. It can be shown that the solution of (1.7) satisfies the exponential
lt
bound
Ix(t)1 ::; ae bt (11¢llc + Ih(u)ldU)
is the continuous norm of ¢. This implies that the Laplace transform of x(t)
exists as well, and as such it can be used to solve the equation. Indeed, if
we take the Laplace transform of (1.7), together with the initial condition
(1.8), we obtain
where X(s) is the Laplace transform of x(t). Solving for X(s), we obtain
where
~(s) = s - ao - ale- rs (1.10)
is the characteristic quasipolynomial of the system. Let <I>(t) denote the
inverse Laplace transform of l/~(s),
h(t) = 0
and the initial condition
x(O) = { ~ 0=0,
-r ::; 0 < O.
1.2 A simple time-delay system 7
Given the fundamental solution (1.11), it is easy to see, using the convolu-
tion theorem, that the solution of (1.7) under an arbitrary initial condition
¢(t) and forcing function h(t) can be expressed as
Equation (1.12) is called the Cauchy formula (also known as the variation-
of-constant formula) for the system, which explicitly expresses the solution
in terms of the initial condition ¢ and the forcing function h.
From (1.12) it is clear that the growth of the general solution for any
bounded h is intimately related to the rate of exponential growth of the fun-
damental solution q>, which is determined by the poles (also known as the
characteristic roots) of the system, i.e., the solutions of the characteristic
equation
~(s) = O. (1.13)
Unlike for systems without delay, generally (1.13) has an infinite number
of solutions. According to complex variables theory, however, since ~(s) is
an entire function, it cannot have an infinite number of zeros within any
compact set lsi::; R, for any finite R > O. Therefore, "most" of the system's
poles go to infinity.
To understand the distribution of the poles, notice that (1.10) and (1.13)
imply that
(1.14)
When s ----> 00, the left hand side of (1.14) approaches infinity, and therefore
the right hand side must also approach infinity. But this means that
lim Re(s) = -00.
s---oo
Therefore, for any given real scalar a, there are only a finite number of
poles of the system with real parts greater than a. With this in mind, let
Si, i = 1, 2, ... be the poles of the system, and any a such that
i = 1, 2, ....
Using the residual theorem from complex variables theory, we may write
(1.11) as
2:
['JI.+jOO
q>(t) = Jo . e- st ~ -l(s)ds + Res[e sit ~ -l(Si)]. (1.15)
a.-JOO Re(si»a.
for t > O.
For the particular time-delay system (1. 7), we may now summarize the
above discussion in the following proposition.
8 1. Introduction to Time-Delay Systems
Proposition 1.1 For any a E JR, there are only a finite number of poles
with real parts greater than a. Let h(t) == 0 and Si, i = 1, 2, "', be the
poles of the system (1. 7), and let
ao = m?J(Re(si) .
•
Then for any a > ao, there exists a L > 0 such that the solution of (1.7)
with the initial condition (1.8) satisfies the inequality
Thus, in order for the solution to approach zero under an arbitrary initial
condition, it is sufficient that all the poles of the system have negative real
parts; in fact, this is also neccesary. This conclusion is the same as that for
systems without delay.
where x(t) E JRn and f : JRxC --+ JRn. Equation (1.16) indicates that the
derivative of the state variables x at time t depends on t and x(~) for
t - r ::; ~ ::; t. As such, to determine the future evolution of the state, it
is necessary to specify the initial state variables x(t) in a time interval of
length r, say, from to - r to to, i.e.,
Xto = ¢, (1.17)
Xl (t) y(t),
X2(t) y(t).
It is important to point out that in an RFDE, the highest order derivative
does not contain any delayed variables. When such a term does appear, then
we have a functional differential equation of neutral type. For example,
I. If the solution also satisfies the initial condition (1.17), we say that it
is a solution of the equation with the initial condition (1.17), or simply a
solution through (to, ¢). We denote it as x(to, ¢, I) when it is important
to specify the particular RFDE with the given initial condition. The value
of x(to, ¢, I) at t is denoted as x(t; to, ¢, I). We will omit I to write x(to, ¢)
or x(t;to,¢) when I is obvious from the context.
A fundamental issue in the study of both ordinary differential equations
and functional differential equations is the existence and uniqueness of so-
lution. It is not our aim in this book, however, to delve into this issue for
RFDEs. Instead, we state the following theorem without proof.
In nearly all cases the systems treated in this book satisfy the conditions
in this theorem, and hence the existence and uniqueness of the solutions
are ensured. Exceptions will be pointed out explicitly.
In this definition, the vector norm 11·11 represents the 2-norm 11.11 2 , It should
be pointed out that other norms such as I-norm and oo-norm can also be
1.4 Stability of time-delay systems 11
x(t) = °
Definition 1.1 For the system described by (1.16), the trivial solution
°
is said to be stable if for any to E JR. and any E > 0, there
exists a 0 = o(to, E) > such that Ilxto lie < 0 implies Ilx(t)11 < E for t ::::: to·
°
It is said to be asymptotically stable if it is stable, and for any to E JR. and
any E > 0, there exists a oa = Oa(to, E) > such that Ilxtolle < Oa implies
lim x( t) = 0. It is said to be uniformly stable if it is stable and 0(to, E)
t->oo
°
can be chosen independently of to. It is uniformly asymptotically stable if
it is uniformly stable and there exists a Oa > such that for any TJ > 0,
there exists aT = T(Oa,TJ), such that Ilxtolle < 0 implies Ilx(t)11 < TJ for
t ::::: to + T and to E JR.. It is globally (uniformly) asymptotically stable if
it is (uniformly) asymptotically stable and Oa can be an arbitrarily large,
finite number.
One should note that the stability notions herein are not at all different
from their counterparts for systems without delay, modulo to the different
assumptions on the initial conditions. Note also that in this book, we shall
only be concerned with asymptotic stability. Thus, while the Lyapunov-
Krasovskii theorem and Razumikhin theorem will be stated in complete-
ness in the next two sections, inclusive of both stability and asymptotic
stability, elsewhere an asymptotically stable system is often simply said to
be "stable." Correspondingly, we will sometimes refer to a system that is
not asymptotically stable as "unstable."
gives rise to
Proof. For any c > 0, we can find a sufficiently small 8 = 8(c) > 0 such
that v(8) < u(c). Hence, for any initial time to and any initial condition
Xto = ¢ with 11¢lle < 8, we have V(t,Xt) :::; 0, and therefore V(t,Xt) :::;
V(t, ¢), for any t ~ to. This implies that
u(llx(t)11) :::; V(t, Xt) :::; V(to, ¢) :::; v(II¢lie) :::; v(8) < u(c),
which in turn implies that Ilx(t)11 < c for t ~ to. This proves the uniform
stability.
To prove uniform asymptotic stability, let c > 0 and 8a > 0 such that
v(8 a ) < u(c). Then IIxtolle :::; 8a implies that IIx(t)1I < c for t ~ to. For this
8a and any TJ > 0, we need to show that there exists aT = T(8 a , TJ), such
that IIx(t)1I < TJ for t ~ to + T. Let 8b > 0 satisfy v(8 b ) < u(TJ). Then, it
suffices to show that IIXto+Tlle < 8b , which implies u(lIx(t)11) :::; V(t,Xt) <
v(8 b ) < u(TJ) for t ~ to+T and therefore IIx(t)1I < TJ. We establish this fact
by contradiction. Suppose that this is not true; in other words, no such T
exists. It follows that IIxtlle ~ 8b for all t ~ to. This means that a sequence
{tk} exists such that
to + (2k - l)r :::; tk :::; to + 2kr, k = 1,2"" ,
1.4 Stability of time-delay systems 13
and
Ilx(tk)112 Db.
For a sufficiently large L, due to the assumption on f, we have 1I±(t)1I =
Ilf(t,xt)11 < L for all t 2 to. Therefore, for t E h = [tk - ft,tk + due ft],
to the Mean Value Theorem, for some B E [0,1],
But this implies that V(tk, Xtk) < 0 for a sufficiently large k, which is a
contradiction. This proves the uniform asymptotic stability.
Finally, if lims-->CXl u(s) = 00, then Da above may be arbitrarily large,
and E can be chosen after oa is given to satisfy v(oa) < U(E), and therefore
global uniform asymptotic stability can be concluded. •
It is obvious from the above proof that u, V, W, and V (t, .) need only be
defined in a neighborhood of zero except for the case of global stability.
Notice also that the lower bound of V need only to be a positive function
of q';(0).
serves to measure the size of Xt. IfV(x(t)) < V(Xt), then Vex) > 0 does not
make V(Xt) grow. Indeed, for V(Xt) to not grow, it is only necessary that
V(x(t)) is not positive whenever V(x(t)) = V(Xt). The precise statement
is as follows.
Vet, x(t)) ::; -w(llx(t)ID whenever Vet + e, x(t + e)) ::; Vet, x(t)) (1.24)
eE
° °
for [-r, OJ, then the system (1.16) is uniformly stable.
If, in addition, w(s) > for s > 0, and there exists a continuous nonde-
creasing function pes) > s for s > such that condition (1.24) is strength-
ened to
Vet, x(t)) ::; -w(llx(t)ID if Vet + e, x(t + e)) ::; p(V(t, x(t))) (1.25)
e
for E [-r, 0], then the system (1.16) is uniformly asymptotically stable.
If in addition lim s _ co u(s) = 00, then the system (1.16) is globally
uniformly asymptotically stable.
°
(1.25), we need to show that for any sufficiently small 8a (such that there
exists an c > to satisfy v(8 a ) = uk)), and any arbitrarily given ""
", < c, there exists aT = T(8 a ,,,,), such that Ilxtll < ", for t ~ to + T,
°
<
that
implies
V(t, x(t)) ~ Vi + (N - k - l)a for t ~ tk + ah, (1.27)
especially, it implies (1.27) for tk+l - r ~ t ~ tk+!' An induction then
shows that V(t,x(t)) ~ Vi < u(1]) for t ~ tN = to + N(ah + r). This
implies jjx(t)jj ~ 1] for t ~ to +T, T = N(ah+r), which will complete the
proof.
Assume that V(tk' X(tk)) > Vi+(N -k-l)a. Then, according to condition
(1.25), there exists ad> 0 such that Vet, x(t)) ~ -'Y and
F(t,O) = 0
16 1. Introduction to Time-Delay Systems
and
L(t)¢ = lOr do [F(t, O)]¢(O). (1.29)
Here in (1.29), the Stieltjes integral is required in general, and the subscript
o is used to indicate that () (rather than t) is the integration variable. As
such, a general linear RFDE can be represented as
where Ak(t) and A(t,O) are given n x n real continuous matrix functions,
and rk(t) are given continuous functions representing time-varying delays,
which can be ordered with no loss of generality, so that
The fundamental solution <p(t, to) of the RFDE (1.30) is the n x n matrix
function satisfying the homogeneous equation
Note that although this initial condition does not satisfy the existence
condition in Theorem 1.2, the existence and uniqueness of the solution can
nevertheless be established. In fact, with the fundamental solution so given,
the solution of the linear RFDE (1.30) with the initial condition Xto = ¢
can be expressed as
x(t;to,¢) = <p(t,to)¢(O) + it to
<p(t,a)da[G(a;to,¢, h)], (1.33)
1.6 Linear time-invariant systems 17
where
Thus, for a linear RFDE, the stability of every solution is equivalent to the
stability of its trivial solution. For this reason, we will simply say that a
system is stable, instead that any particular solution is stable. It should
also be clear that a linear system is stable if and only if it is globally stable.
In this case, the fundamental solution <l>(t, to) only depends on t - to. As a
result, without loss of generality, we only need to consider to = O. We write
<l>(t) for <l>(t, 0).
Like their counterparts without delay, LTI time-delay systems can be
studied effectively using frequency domain methods. This distinguishing
advantage is made available by such frequency domain analysis tools as
Laplace transforms. The case study given in Section 1.2 has shed light on
how this may proceed, and is now extended below to general LTI systems.
Let
xo = ¢ (1.36)
be the initial condition. Taking the Laplace transform of (1.35) with the
initial condition (1.36), we obtain
18 1. Introduction to Time-Delay Systems
where X(s) and H(s) are the Laplace transform of x(t) and h(t), respec-
tively,
where
~(s) = sf - lOr eOSdF(B) (1.38)
det[~(s)] =0 (1.39)
is called the characteristic equation, and the expression det[~(s)] the char-
acteristic function, or alternatively, the characteristic quasipolynomial. The
solutions to (1.39) are called the characteristic roots or poles of the system.
Let h(t) = 0 and choose the initial condition (1.32). Then the right hand
side of (1.37) becomes ~ -1 (s). From this and the definition offundamental
solution, it follows that
C[<I>] = ~-l(s),
i.e., ~ -1 (s) is the Laplace transform of the fundamental solution. This
shows that the fundamental solution also satisfies the equation
Theorem 1.5 For any real scalar ,,(, the number of the solutions, count-
ing their multiplicities, to the characteristic equation (1.39) with real parts
greater than "( is finite. Define
(i) The LTI delay system (1.35) is stable if and only if ao < o.
(ii) FOT any a > ao, there exists an L > 1 such that any solution x(t) of
(1.35) with h(t) = 0 and the initial condition (1.36) is bounded by
(1.41)
m
= Po(s) + LPk(s)e- hkS ,
k=l
+ L LPkiy(i)(t -
n-l m
y(n)(t) h k ) = 0,
i=O k=O
20 1. Introduction to Time-Delay Systems
a(s, e- T) = Lak(s)e-kTS.
k=O
It is worth noting that for systems with incommensurate delays the char-
acteristic quasipolynomial is in effect a multivariate polynomial of several
variables, while for those with commensurate delays, it may be viewed
as a bivariate polynomial, or even further, a polynomial with coefficients
themselves as polynomials. This will prove to be a fundamental distinction.
Indeed, it will be shown in the subsequent chapters that for systems with
commensurate delays, the characteristic quasipolynomial is much easier to
analyze.
xo = ¢. (1.43)
D.(s) = 0 (1.44)
where
me/> = -r::;t::;O
max (11¢(t)11 + 11¢(t)ll)·
We will not give proof here. Interested readers are referred to the litera-
ture discussed in the Notes section.
It is interesting to compare the above theorem for neutral delay systems
to Theorem 1.5 for retarded delay systems. First, while there can be only a
finite number of characteristic roots on the right hand side of any vertical
line Re( s) = 'Y in a retarded delay system, this is not necessarily the case for
a neutral delay system. Second, in the definition of exponential exponent,
"max" in (1.40) is replaced by "sup" for a sufficient condition for stability
of neutral delay systems. In other words, in order to satisfy condition (i)
in the above theorem, all the characteristic roots need to be on the left
hand side of some vertical line Re( s) = -a for some a > O. Having all the
characteristic roots on the left hand side of the imaginary axis is a necessary
but not sufficient condition for the stability of neutral delay systems. There
exist neutral delay systems with all the roots on the left hand side of (but
may have roots arbitrarily close to) the imaginary axis, and the systems are
not asymptotically stable as defined in Definition 1.1. For retarded time-
delay systems, this complication does not occur since there are only a finite
number of characteristic roots on the right hand side of any given vertical
line. Third, for neutral time-delay systems, the stability exponent ao is not
necessarily continuous with respect to delays as they approach zero.
22 1. Introduction to Time-Delay Systems
where
J.Ll > J.L2 > ... > J.Lp ;::: O.
A strip with J.Lk > 0 is of retarded type, and a strip with J.Lk = 0 is of
neutml type. Clearly, within a retarded root strip, lim Re(s) = -00. For
8-+00
a neutral root strip, as s - 00, it may approach a vertical line, including
the imaginary axis. For a retarded delay system, all the root strips are of
retarded type. A neutral delay system may (but does not always) possess
a neutral root strip. The proof of this property can be carried out by a
rather sophisticated study of the relationship of .6.( s) and its comparison
function, and will not be discussed here.
sizes the critical role that the continuity of the stability exponent plays in
stability analysis. Second, it helps reveal the weakness of classical stability
results. These two aspects thus provide the guide for the ensuing devel-
opment, leading to the frequency-sweeping criteria of Section 2.3 and the
constant matrix tests of Section 2.4, both of which are necessary and suffi-
cient conditions for delay-dependent and delay-independent stability alike,
and both require only the computations of matrix pencils. The frequency-
sweeping tests are readily obtained based upon simple algebraic manipula-
tions, and are reminiscent of small gain conditions. On the other hand, the
derivation of the constant matrix tests draws heavily upon the theory of
polynomials, notably the results concerning the Schur-Cohn stability test.
While frequency-sweeping criteria are to be performed on the basis of fre-
quency gridding, the constant matrix tests require only the computation of
constant matrices.
Chapter 3 studies systems with incommensurate delays, which is both
motivated by and built upon a small gain approach. This chapter begins
with a concise introduction of an operator theoretic system description as
well as a statement of the small gain theorem, the rudiments needed in
a small gain approach. It continues with a brief narrative of the struc-
tured singular value, which is a key measure in robust control theory and
in the characterization of stability for LTI systems with incommensurate
delays. Indeed, it is shown that the necessary and sufficient condition for
delay-independent stability requires computing a structured singular value;
similarly, a sufficient condition is also obtained in the same manner for
delay-dependent stability. While these provide useful conceptual results,
the structured singular value is known to be difficult to compute, and in
fact poses an intractable computational problem in general. Consequently,
one is led to the contemplation on the inherent computational complexity
of the stability problem, for which we give a formal analysis drawing upon
concepts and techniques typically found in computing theory and operation
research, with such notions as NP-hard and NP-complete problems. The
analysis is preceded by a brief summary of necessary materials drawn from
complexity theory, and it led to the classification of the stability problem
as NP-hard. In light of this difficulty, we then develop a number of readily
computable sufficient conditions. The chapter ends with an extension to
neutral LTI delay systems, with both commensurate and incommensurate
delays.
Chapter 4 examines the robust stability of LTI delay systems described
by uncertain quasipolynomials with incommensurate delays. By robust sta-
bility we mean that the stability is maintained for the entire family of qua-
sipolynomials as the coefficients vary within a prescribed set. Fundamental
robust stability concepts and results such as the zero exclusion principle,
the edge theorem, and the phase growth condition are extended to and fur-
ther developed for uncertain quasipolynomials. More specifically, the chap-
24 1. Introduction to Time-Delay Systems
ter first presents the zero exclusion principle, which provides a general,
geometrical characterization on the zeros of uncertain quasipolynomials.
This leads to the first main result of the chapter, namely, the edge theorem.
The result is analogous to its counterpart in robust stability theory, and it
states that for a polytopic family of quasipolynomials, robust stability is
achieved if the edge members of the family are all stable. It is worth noting
that in the robust stability literature, the polytopic family remains the most
general uncertain polynomial for which a tractable necessary and sufficient
stability condition is available. As such, the edge theorem for the polytopic
quasipolynomial presents an equally strong result. The chapter continues
to identify some classes of poly topic quasipolynomial families whose sta-
bility is implied by the stability of a finite number of quasipolynomials in
this family. Furthermore, it presents a multivariate polynomial approach
for the interval and diamond multivariate polynomials, which in turn pro-
vide robust stability conditions for the interval and diamond families of
quasipolynomials.
Beginning with Chapter 5, we present time domain stability and robust
stability conditions. This is where the Lyapunov-Krasovskii Theorem and
the Razumikhin Theorem come to play their roles. Chapter 5 treats sys-
tems with a single delay, but without uncertainty. The chapter discusses
various Lyapunov-Krasovskii functionals as well as Razumikhin-type sta-
bility criteria, resulting in various delay-independent and delay-dependent
stability conditions, mostly posed as solutions to LMIs. The relationships
between various results are also explored. It is shown that the use of model
transformation in deriving some of the delay-dependent stability results
will inherently lead to conservatism, because of the resultant additional
dynamics. It is also shown that for necessary and sufficient stability con-
ditions, a complete quadratic Lyapunov-Krasovskii functional needs to be
used. A discretized Lyapunov functional method is used to reduce such
a Lyapunov-Krasovskii functional to a finite-dimensional LMI. Numerical
results show that conservatism is usually very small even for rather coarse
discretization.
Chapter 6 presents extensions to uncertain delay systems. Unlike Chap-
ter 4, the system uncertainty description herein is given on the system's
state matrices, which are either poly topic or norm-bounded. Both these
uncertainty descriptions are widely found in state-space models. Robust
delay-dependent and delay-independent stability conditions that also ad-
mit the form of LMIs are developed in parallel to Chapter 5. In particular,
for polytopic uncertainty, these LMIs are defined at the vertices of the
uncertain family.
Chapter 7 extends the results of Chapters 5 and 6 to systems with multi-
ple delays and distributed delays with piecewise constant coefficients, which
mayor may not have uncertainty. The ideas are very similar to the previous
two chapters but are technically more delicate.
1.9 Notes 25
1.9 Notes
1.9.1 A brief historic note
Time-delay systems are also known as hereditary systems, systems with
aftereffects, or systems with time lags. The first functional differential
equations were considered by, among others, such great mathematicians
as Euler, Bernoulli, Lagrange, Laplace, and Poisson in the eighteenth cen-
tury, arising from various geometry problems. In the early twentieth cen-
tury, a number of practical problems were also modelled using FDEs. These
include viscoelasticity problems studied in 1909, the predator-prey model
used in population dynamics in 1928-1931, both by Volterra, mathemat-
ical biology problems studied in 1934 by Kostyzin, and ship stabilization
problems in 1942 by Minorsky.
On the stability of time-delay systems, Pontryagin obtained some fun-
damental results concerning the zeros of quasipolynomials in 1942, and
Chebotarev published a number of papers devoted to the Routh-Hurwitz
criterion for quasipolynomials in early 1940s. A paper of Myshkis in 1949
26 1. Introduction to Time-Delay Systems
formulated for the first time the initial value problem. Those attempting to
extend Lyapunov's theory to time-delay systems initially encountered some
difficulties, which were not resolved until the work of Krasovskii. Krasovskii
was the first to emphasize the importance of adopting Xt rather than x(t)
as the state, and to develop further the idea to fruitation. This took place in
1956. The idea of bypassing functionals was due to Razumikhin, which led
directly to the Razumikhin Theorem. Further extensions along these lines
were made by Kolmanovskii and Nosov to neutral functional differential
equations. For a more comprehensive account of these early developments,
we refer to the book by Kolmanovskii and Nosov [152]. Other treatises of
time-delay systems include the books by Bellman and Cooke [9], G6recki
et al. [78], and Hale and Verduyn Lunel [105], who, like Kolmanovskii and
Nosov, detail the developments of their times and additionally cover other
aspects of time-delay systems. The recent monograph by Niculescu [204]
also documents briefly the historical developments.
Since 1990s there has been a substantial increase of research activities
on time-delay systems in the systems and control community. The main
achievements in this period are characterized by computational improve-
ment and in the study of robust stability. These advances form the core of
the later chapters, where brief reviews and references to the latest progress
will be provided accordingly.
Frequency Domain
Approach
2
where A o, Ak E IR nxn are given system matrices, and rk are delay times.
As discussed in Chapter 1, the stability of this system is fully determined
by its characteristic function
Specifically, the system is stable if and only if p (Sj e- r1S , " ' , e- rKS ) has
no zero, or root, in the closed right half plane C+. We state this formally
below, as a definition.
Hence, the delay system (2.1) is stable independent of delay if the stability
persists with respect to all possible nonnegative delays. On the other hand,
if it is stable only for a subset of nonnegative delays, then we say that the
stability is delay-dependent.
It is important to recognize that the characteristic function of a time-
delay system defines a real quasipolynomial of s. This feature makes it
possible to extend a number of key properties of polynomials to the study
of time-delay systems. Indeed, a critical ingredient in our subsequent de-
velopment is based on the continuity of the stability exponent with respect
to delays, a fact we have stated in Chapter 1. For a glimpse of its impor-
tance at the outset, let the system (2.1) be stable at the delay values r k,
k = 1, 2, "', K. Then by the continuity property of the stability ex-
ponent, the system will remain stable within a neighborhood of r k. Much
of our effort in determining delay-independent and delay-dependent sta-
bility then amounts to expanding that neighborhood, so that the stability
exponent remains negative. In turn, this amounts to finding the critical
delay values at which the characteristic roots intersect the stability bound-
ary, i.e., the imaginary axis, thus rendering the system unstable. In other
words, when the delay values deviate from r k, we want to determine the
smallest deviation of rk from r k, such that
(2.4)
our frequency domain stability analysis of LTI delay systems, in both this
chapter and Chapter 3. We shall further elaborate this point in the next
section, using systems with commensurate delays.
For the rest of this chapter, we shall consider LTI delay systems with
commensurate delays. Systems in this class can be represented by
m
(2.11)
This gives a general formula for computing the delay margin. More gener-
ally, we may assume, with no loss of generality, that
"11 < 'TJ2 < ... < 'TJN'
It follows that the system is stable for all T E ('TJi' 'TJi+l) whenever it is
stable for some T* E ('TJi' 'TJi+l)' This then allows us to ascertain a system's
°
stability in the full range of delay values, beyond the interval determined
by T. It also indicates that the stability at T = can be made without loss
of generality.
With the bivariate polynomial representation (2.9), it is evident that T
can be found by solving the imaginary roots s E ac+ and the unitary roots
z E allJ) of a(s, z), giving rise to a stability criterion commonly referred
to as a two-variable criterion. The two-variable criterion appears to be
the origin of many classical stability tests for systems with commensurate
delays, which attempt to solve the bivariate polynomial (2.9) in one way or
another. In particular, most of the classical tests attempt to accomplish this
by eliminating one variable, thus converting the stability problem to one
free of delay, and seeking the solution of polynomials of one single variable.
In the next section we discuss a number of tests in this spirit. These sample
tests are representative of classical results and should give the reader the
essential flavor of the two-variable criterion.
which maps 8 from the open right half complex plane C+ to A in the open
unit disk D. Construct the 2-D polynomial
a(8, z) 0, (2.13)
0,(8, z) O. (2.14)
When no solution exists, and when the system is stable in the delay-free
case, it must also be stable independent of delay. Otherwise, when the
two equations do admit a common solution, it is possible to eliminate one
variable, resulting in a polynomial of one single variable. For example, we
may eliminate 8 and obtain a polynomial in z, say b(z). If b(z) has no
unitary root, we may again conclude that the system is stable independent
of delay. Otherwise, we may proceed to find all the unitary roots Zi of b( z).
There are only a finite number of such roots since b( z) is a polynomial.
For each Zi, a(8, Zi) is a polynomial of the variable 8, which too admits
only a finite number of possible roots 8i E ac+. Thus, we find all the roots
(8i' Zi) E ac+ x aD such that a(8i' Zi) = O. Since the bivariate polynomial
satisfies the conjugate symmetry property, only those 8i on the positive
imaginary axis need to be considered. Thus, consider 8i = jWi, Zi = e-()i,
36 2. Systems with Commensurate Delays
where Wi > 0, and Bi E [0, 211'J. The delay margin can then be determined
using (2.11).
a(s, z)
o'(s, z)
Eliminating s yields the polynomial equation
z4 + z3 + z + 1 = (z + 1)2(z2 - Z + 1) = 0,
which has the solutions Zl = -1, Z2 = (1 + j.;3)/2, and Z3 = (1- j.;3)/2.
Correspondingly, Sl = 0, S2 = -j.;3, and S3 = j.;3. As a result, T is found
from S3 and Z3 = e- j 'lr/3, as T = 11'/(3.;3).0
c(T, s) := (1
1-
+ sT)qa ( s, 1 + sT .
ST)
Since c(O, s) = a(s, 1), the polynomial c(O, s) is stable by assumption.
Furthermore, for any 0 ::; T < 00, the factor (1 + sT)q has no effect on the
imaginary roots of a(s, ~+~~). Thus, for all 0 ::; T ::; 00, the imaginary
roots of a (s, ~+~~) consist of those of c(T, s), 0 < T < 00, and those
2.2 Some classical stability tests 37
of a(s, - 1); the latter polynomial accounts for the case T = 00. It then
suffices to find the imaginary roots of a( s, -1), and for finite values of T to
calculate the roots of c(T, s), which is a polynomial in s whose coefficients
are parameterized by T. Clearly, if neither of the two polynomials has any
imaginary root, then a(s, ~+:~) has no imaginary zero for all ~ T ~ °
00. As a result, a(s, z) has no root in 8C+ x 8][J), and hence the delay
system is stable independent of delay. However, if for some T = Ti > 0,
the polynomial c(1i, s) has an imaginary root Si = jWi, Wi > 0, then there
exists a unitary Zi = e- jOi E 8][J), ()i E [0, 27rJ, where
T ?: 0. (2.17)
a (s, e- TS ) = s2 + 4s + 4 - 0.25e- TS •
38 2. Systems with Commensurate Delays
Hence for any 0 < T < 00, c(T, s) is stable. Consider next a(s, - 1) =
s2+ 1, which has a pair of imaginary roots s = ±j. This leads to "7 = 11". 0
(2.19)
2.2 Some classical stability tests 39
(2.20)
the delay margin can be duly computed according to (2.11). When the
quasipolynomial is stable for T = 0 and (2.20) admits no solution, the
system is stable independent of delay, corresponding to the conditions
. ) I< 1,
Iaoal(jw)
(JW Vw>O, . (2.21)
I>
or
I:~g:~ 1,
Vw>O.
The latter, however, is not possible since ao(s) is a higher order polynomial
than al(s). Thus, we conclude that the necessary and sufficient condition
for stability independent of delay is that (2.21) holds, and that a(s, 1) =
ao(s) + ales) is stable. In this case, ao(s) is also stable, corresponding to
T= 00.
Tsypkin's Theorem states that under the condition that G(s), or equiv-
alently ao(s), is stable, then the closed-loop system will be stable for all
T ~ 0 whenever IG(jw)1 < 1 for all w E JR., which coincides with the
above sufficient condition. Finally, we note that the condition provides a
frequency-sweeping test in the form of a small gain theorem, which can
be checked either numerically by inspecting the equation (2.20), or graph-
ically by plotting the frequency-dependent measure lal(jw)/ao(jw)l. The
small gain theorem will play a predominant role in our development of
frequency-sweeping tests and will be discussed in Chapter 3.
where a and b are real constants. The quasipolynomial for this system is
Under the assumption that a+b > 0, the system is stable at T = O. Suppose
that this is the case. Then, the equation (2.20) is given by
w 2 + a 2 - b2 = 0,
which may have a nontrivial solution only when lal < Ibl, yielding the zero-
crossing frequency WI = )b2 - a 2. Clearly, this is possible only when either
a > 0, b > 0, or a < 0, b > 0. In the former case,
(}1 . b
= L. JWl +a - 11" = tan -1 (WI)
-I
a
I '
and thus
cos- 1 (¥)
T = -v'---r:b2~_=a=;;2:-
On the other hand, the system is stable independent of delay if and only if
a + b > 0, and a ~ Ibl. Indeed, the condition a ~ Ibl guarantees that
Note that when a = b > 0, lal(jW)/ao(jw)1 < 1 for all w > 0, but
lal(jw)/ao(jw)1 =
1 at w = 0. In this case, the system remains to be
stable independent of delay.
To further illustrate, consider the system given in Example 2.3, whose
quasipolynomial is in the form of (2.19), with
ao(s) = s2 + s + 1, al(s) = s.
The equation (2.20) is simplified to
w4 - 2w 2 + 1 = 0,
which yields the zero-crossing frequency WI = 1. In addition,
L. a1 (j) = 0.
ao(j)
Hence, 01 = 1r. This leads to the same T = 71', as obtained in Example 2.3.
o
In the case of multiple commensurate delays, the following iterative calcu-
lation can be employed to extend the above test. With the quasipolynomial
(2.6), define the new quasipolynomial
a(l) (s, e- rs ) :=ao(s)a(s, e- rs ) -aq(s)e-qTSa(-s, eTS ).
A simple calculation yields
q-l
L a~I)(s)e-kTS
k=O
q-l
L [ao( -s)ak(s) - aq(s)aq_k( -s)] e- kTS .
k=O
Furthermore,
q-l
a(l) (-s, eTS ) = L [ao(s)ak(-s) - aq(-s)aq_k(s)] ekTS .
k=O
Note that a(1) (s, e- rs ) and a(1) (-s, eTS ) may be considered as a frequency-
dependent linear transform of a (s, e- rs ) and a (-s, eTS ),
a (jw, e- jTW ) 0,
a (-jw, e jTW ) = 0,
42 2. Systems with Commensurate Delays
laq(jw)
ao(jw)
1< 1, Vw E JR.,
the two sets of roots coincide. Nevertheless, the order of a(1) (s, e- TS ) in
e- TS has been reduced to q-1. We may continue this procedure, by defining
That is,
q-i
a(i+l) (s, e- TS ) = :~:::>~i)(s)e-kTS,
k=O
with
ak(s),
a(i)
o
(-s)a(i)
k (s) _ a(i)
q-t.(s)a(i).
q-t-k (s).
For i =q- 1, the procedure generates the function
This effectively reduces the problem to that with one single delay, from
which all zero-crossing frequencies of a(q) (s, e- TS ), and hence those of
a (s, e- TS ), can be computed.
Example 2.5 Let us invoke the iterative procedure to test the stability of
the system (2.15). Note that ao(s) = s, al(s) = a2(s) = 1. Construct
1. The 2-D tests require that one variable be eliminated and a resultant
polynomial generated from two bivariate polynomials. This can be
difficult for high order systems with many delays, and there exists no
systematic procedure for performing this elimination.
Routh array. With the free parameter T, generating the Routh array
can be a cumbersome task.
3. Likewise, the direct method requires a complicated iterative proce-
dure which involves algebraic manipulations of polynomials.
Most notably, the elimination step in each of the three tests requires sym-
bolic computation, or calculation "by hand," which, from a computational
perspective, is highly undesirable. Thus, for high order systems with many
delays, their implementation may encounter a varying degree of difficulty,
rendering the tests inefficient and less applicable.
(2.24)
Theorem 2.1 The system (2.23) is stable independent of delay if and only
if
(i) Ao is stable,
(ii) Ao + Al is stable, and
(iii)
't/w > 0, (2.25)
where p(.) denotes the spectral radius of a matrix.
and (ii), respectively. Assume then that Ao is stable, so that (jw1 -AO)-I Al
is well defined for all w E R Suppose that (2.25) holds. This means that
for any eigenvalue Ai ((jw1 - Ao)-IA I ) of the matrix (jw1 - AO)-IA I ,
i = 1, ... , n,
Vw > 0,
Vw > 0
it follows that
Vw>O,
or equivalently,
Vw > O.
Moreover, note that the condition (ii) precludes the possibility of det(Ao +
AI) = O. Therefore, we have shown that for all wE JR,
that is, the quasipolynomial of the system does not intersect the imaginary
axis. Hence, the system is stable independent of delay. This proves the
sufficiency.
To establish the necessity, it suffices to show that (iii) is necessary. To-
ward this end, assume first that p ((jwo1 - AO)-I AI) = 1 for some Wo > O.
This implies that the matrix (jwoI - AO)-I Al has an eigenvalue ejOo , for
some eo E [0, 21l'J. Let TO = eo/woo It is clear that
or
det (jwo1 - Ao - Ale-jrowo) = O.
As such, the system becomes unstable at T = TO, and hence cannot be
stable independent of delay. Next, suppose that p ((jw1 - AO)-I AI) > 1
for some w > O. Since p ((jw1 - AO)-I AI) is a continuous function of w,
and since
46 2. Systems with Commensurate Delays
there exists some Wo E (w, 00), such that p ((jwof - AO)-1 AI) = 1. This,
as we have shown above, implies that the system cannot be stable inde-
pendent of delay. Consequently, we conclude that the conditions (i-iii) are
necessary for the system to be stable independent of delay, thus completing
the proof. _
The frequency-sweeping test given in Theorem 2.1 may be considered an
extended small gain condition, to be discussed later in Chapter 3. Stability
conditions of this type are routinely found in robust control theory and
are generally held as efficient measures in stability analysis; a subsequent
example will attest to this point. Evidently, our derivation of this test still
draws fundamentally upon the continuity property of the stability exponent
of quasipolynomials. Yet unlike the classical results, our test gets rid of the
variable elimination procedure and lends a readily implementable criterion.
The test can be easily checked by computing essentially the frequency-
dependent measure p ((jwf - AO)-1 AI)' which is amenable to computation
because the ease of computing the spectral radius.
While this frequency-sweeping test is more numerically inclined and in-
deed has a justifiable advantage in computation, as a small gain type con-
dition it may also be used for analysis purposes. This will be shown in
forthcoming examples. Moreover, it is of interest to point out that since
(sf - AO)-1 Al defines an analytic function in ~\, the spectral radius
p ( (sf - Ao) -1 AI) enjoys stronger properties than merely the continuity
required in establishing stability. In fact, it will be shown in Chapter 3
that, whenever Ao is stable, p ((sf - A O)-IA 1) is a sub harmonic function
in C+, and so is the function p ((sf - AO)-1 A 1e- rs ). One direct conse-
quence of this fact is that not only is p ((sf - AO)-1 A 1e- rs ) continuous in
C+, for any T :2: 0, but also it satisfies the well-known maximum modulus
principle; that is, it achieves its maximum on the boundary 8C+ of C+. In
other words,
Example 2.6 Let A E Rnxn be a stable matrix and consider the delay
system
x(t) = A x(t) +,BA x(t - T), T :2: 0, ,B E R (2.26)
We want to determine the set of ,B and A such that the system is stable
independent of delay. Note first that the condition (ii) dictates that,B > -1.
2.3 Frequency-sweeping tests 47
Hence, if Ai E R for all i, we have p ((3(jw1 - A)-1 A) < 1 for all w > 0 and
all (3 E (-1, 1]. In other words, the system is stable independent of delay
for all (3 E (-1, 1]. On the other hand, if A has complex eigenvalues, then
It is worth noting that this example has been studied elsewhere using clas-
sical stability tests, but the analysis was, understandably, restricted to 2 x 2
matrices. In general, for a matrix A of an arbitrary dimension, the classical
analysis becomes highly nontrivial. 0
Ao = ( -2 0)
0.5 -2 '
As a result,
Evidently, p ((jwI - Ao)-1 A 1) < 1 for all w > O. This leads to the same
conclusion that the system (2.17) is stable independent of delay. 0
that is, when certain eigenvalues of (jwI - AO)-I Al have a unit modulus.
In general, however, if only delay-dependent stability can be ensured, Ao
may have eigenvalues on ac+. This leads us to the consideration of the
frequency-dependent matrix pencil (jwI - Ao) - AA I . For a matrix pair
(A, B), denote its ith generalized eigenvalue by Ai(A, B). Furthermore,
define
e(A, B) := min {IAII det(A - AB) = O}.
It is a well-known fact that the number of finite generalized eigenvalues for
(A, B) is at most equal to the rank of B. Also, if the rank of B is constant,
then Ai(A, B) is continuous with respect to the elements of A and B.
The idea is then to use the frequency-dependent generalized eigenvalues
Ai(jwI - A o, AI) to replace the eigenvalues Ai ((jwI - AO)-I AI)' which
may not exist at certain frequencies.
Furthermore, define
That is, the system (2.23) is stable for all T E [0, T), but becomes unstable
at T = T.
Proof. Consider first the case T = 00. This corresponds to the condition
for all T ~ 0 and all w > O. Hence, it implies that for any T ~ 0,
Vw E (0, 00).
Ai:, a result, the system is stable for all T E [0, 00); in other words, the
system is stable independent of delay. Suppose now that T < 00. For any
T E [0, T), we claim that
Vw E [0, 00).
This is clearly the case for w =/:: wL for at these frequencies
for all i = 1, "', n. Additionally, at w = w~, since for any T E [0, T),
TW~ =/:: (h by definition, we have det (jwV - Ao - Ale-j7'W~) =/:: 0. This
proves the claim. Again, by an appeal to the continuity argument, we con-
clude that the system is stable for all T E [0, T). On the other hand, if
T = T, then there exists a pair (w~, O~) such that T = 01/wL and
Corollary 2.3 The system (2.23) is stable independent of delay if and only
if
(i) Ao is stable,
(ii) Ao + Al is stable, and
(iii)
Vw>o. (2.27)
°
gin. This can be seen by noting that whenever det (jwV - Ao - Al) =/:: 0,
the inequality < 01 < 21l' holds, and hence Ti > 0. One should note,
unsurprisingly, that the computation of the delay margin is more demand-
ing than the determination of the delay-independent stability. Note also
that the computation needs to be performed only over a finite frequency
interval, since the moduli of the generalized eigenvalues will inevitably ex-
ceed one beyond a certain frequency. This, together with the fact that the
generalized eigenvalues are continuous functions of w, is reassuring from a
numerical standpoint.
50 2. Systems with Commensurate Delays
det (~PkZk)
C pJ-( )}
0 1 0
dcl{Z 1 0 0 1
-Po -PI -Pm - I
Theorem 2.5 The system (2.28) is stable independent of delay if and only
if
(i) Ao is stable,
m
(ii) L:: Ak is stable, and
k=O
(iii)
p(Mm(jw)) < 1, 'v'w > 0, (2.29)
where
(s1 - ~O)-I Al
Mm(s):= ( .
o 1
Theorem 2.6 Suppose that the system (2.28) is stable at 7 = 0, and let
q = rank(Am). Furthermore, define
where
0 1 0
G(s),~ (
0
-(s1 - Ao)
0
Al
1
Am- I
).
H(s) := diag ( 1, ... , 1, -Am ) .
Then,
T:= min ri.
1:::;i:5q+n(m-l)
The system (2.28) is stable for all T E [0, r), but becomes unstable at T = r.
r
0 1 0
( 0
-Po
0
-PI
1
-Pm- 1
-r}
_PO-l Pm-l
~CTP' 0
1
(2.30)
Example 2.8 Consider the system (2.28), with four states and three com-
mensurate delays, given by
0, ~ ).
1 0
0 1
Ao = 0 0
-3 -5 -2
( -0.00 0.005 0.25
Al = 0.~05 0.005
0
0
0
~).
(T
-1 0 -0.5
0.0025 0
A2 =
-1
0
0
0.05
-0.5 -0.5
0 O~ ).
0 0.075
( Oog75 0.05 0.05 0~25 )
A3 0.05 0.05 0 o .
0 -2.5 0 -1
-0.6887 ± 1. 7636j
-0.3113 ± 0.6790j.
n,
1, 2, 3 modified to
(I ~} A2~P
0 0 0 0
0 0 0 0
Al 0 0 0 0
q
0 -1 -2 -1 -1
0 0
U
0 0
A3 =
0 0
-5 0 -2
1.4
1.3
1.2
1.1
~ 1
;3
'2
Ii'
::00.9
0.8
0.7
0.6
0.5
0 0.5 1.5 2 2.5 3 3.5 4 4.5 5
Frequency
LL
n-l q
y(n)(t) + akiy(i)(t - kr) = 0, r 2: O. (2.31)
i=O k=O
(2.31) by defining
1 o
Ao
o 1
-aO,n-l
).
)
o o
o o
-ak,n-l
Corollary 2.7 The system (2.31) is stable independent of delay if and only
if
)
-~ ... _aq_l{s) _ aq{s)
M.(s) ,~ ( ao{s)
~ ...
ao{s)
0
ao{s)
0
1 0
where
o 1
G.(s) ,~ (
Then,
The system (2.31) is stable for all T E [0, r), but becomes unstable at T = r.
Note that the dimensions of these matrices are equal to q, which may be
considerably lower than that of the matrices in Theorems 2.5 and 2.6. As
such, it is possible to check both conditions more efficiently.
p(z) := zq I)kZ-k.
k=O
W := pH (S)p(S) - pH (S)p(S),
It is known that the polynomial p( z) has all of its roots in lID if and only
if W is positive definite. In particular, the roots of p(z), Zi, i = 1, ... , q,
are related to W via the following Orlando formula.
Note that
(I P:,
pq Pq-I ... PI)
:).
o Pq ... P2
(
··· ..
. ' .
o 0 pq
Note also that C is the test matrix in the Schur-Cohn criterion. The fol-
lowing fact establishes a relationship between C and W.
Lemma 2.10
det(W) = (-I)q det(C).
Based on this fact and the Schur determinant formula, we find that
= ( )q (P(S) p(S))T
-1 det pH(S) pH(S)
)
0 0
( a,(,)
aleS) ao(s) 0
~l(S) :=
).
aq_l(s) aleS)
r
( a,(')
aq(s) a2(s)
E,(s) ,~ (2.33)
0 aq(s)
and
~(s) := (
~l(S)
~!j (s)
~2(S)
~{l (s) ).
In view of Lemmas 2.9 and 2.10, we assert that that for any s,
n
det(~(s)) = (-1)q laq(s)1 2q II (1 - ZiZj) , (2.34)
i,j=l
where Zi, i = 1, "', q, are the roots of a(z) for fixed s. Therefore, for
any w > 0, there is a z E aD such that a(jw, z) = a(z) = whenever °
°
det (~(jw)) = 0. The implication is then fairly clear: By solving the equa-
tion det(~(jw)) = 0, we find all such w > that a(z) = a(jw, z) has a
root on aD. Since det(~(jw)) defines a polynomial in w, the solutions can
be found by solving the eigenvalues of a constant matrix. Clearly, there are
only a finite number of solutions, and among them the positive solutions
constitute the zero-crossing frequencies. Once the zero-crossing frequencies
are found, the unitary roots of a(jw, z) can also be determined at each zero-
crossing frequency by solving the complex polynomial a(z) = a(jw, z). The
2.4 Constant matrix tests 59
aq-l,i
aqi
i = 0,1, ... , n - 1,
~ ),
ao, o
ah aOi
(
i = 0, 1, n-l,
(j)iTi (j)i Hi )
(_j)i HT (-j)iTT '
i = 0,1, n.
Furthermore, define
)
o I o
p~ (
).
1 0
F(s) ,~( r 0 1
-ao(s) -al(s) -aq_l(s)
where 0 -=I Wk E o-(P) n IR+ and (h E [0, 27rJ satisfy the relation e- jOk E
o-(F(jWk)' G(jWk))' The quasipolynomial (2.6) is stable for all T E [0, T),
but is unstable at T = T.
60 2. Systems with Commensurate Delays
Proof. Suppose first that a(P) n R+ = 0. This implies that for any
wE R+, det(wI - P) =f:. 0. In view of Lemma 2.4,
det(wI - P) =
L;2(jW) )
L;¥(jw) .
°
then a (jw, e-irw ) =f:. °
at T = 0, it is stable for all T E [0, (0). Note also that if a(P) nR+ = {O},
for all w E R+ and T E [0, (0), despite that
a(O, z) = may hold for some z on 8D. This follows from the fact that
a(jw, z) has no root on 8D for any w =f:. 0, and that a (jw, e-irw ) =f:. °
°
at w = 0, by the assumption that a (s, e- rs ) is stable at T = 0; the latter
implies that a(jw, 1) =f:. for all w E R+. Hence, again, we conclude that
a(s, e- TS ) is stable for all T E [0, (0) if a(P) nR+ = {O}.
°
Suppose now that a(P) nR+ =f:. 0 and a(P) nR+ =f:. {O}. In other words,
P does have positive real roots. Let =f:. Wk E a(P) n R+. Since according
to Lemma 2.4,
q
det(zG(s) - F(s)) = Lai(S)zi = a(s, z),
i=O
°°
As a result, we may conclude that T = 00 if a(F(jwk), G(jWk)) n 8D = 0
for all =f:. Wk E a(P) n R+. If, however, a(F(jwk), G(jWk)) n 8D =f:. 0
for some =f:. Wk E a(P) n R+, then there exist some Zk E 8D such that
Zk E a(F(jwk), G(jWk))' Consequently, there exist some (h E [0, 21r]
°
such that Zk= e- ifh and a (jWk' e- ifh ) = 0, or equivalently, e- j1h E
a(F(jwk), G(jwk))' It follows that a (jWk' e- jrwk ) = for T = (h/Wk,
and hence the system is unstable at T = T. However, for any T E [0, T),
°
since TWk =f:. Ok for all (Wk' Ok) such that e- iIJk E a(F(jwk), G(jWk)),
we have a (jWk' e- jrwk ) =f:. 0. This implies that a (jw, e- jrw ) =f:. for all
wE R+. Therefore, a (s, e- rs ) is stable for all T E [0, T). The proof is now
completed . •
Theorem 2.11 suggests that a two-step procedure can be employed to test
the stability of the quasipolynomial (2.6). First compute the eigenvalues of
the 2nq x 2nq matrix P. If P has no real eigenvalue or only one real eigen-
value at zero, we conclude that the quasipolynomial is stable independent
of delay. If this is not the case, compute next the generalized eigenvalues of
the q x q matrix pair (F(jwk), G(jWk)), with respect to each positive real
eigenvalue Wk of P. If for all such eigenvalues the pair (F(jwk), G(jwk)) has
2.4 Constant matrix tests 61
no generalized eigenvalue on the unit circle, we again conclude that the qua-
sipolynomial is stable independent of delay. Otherwise, we obtain the delay
margin T. Note that if aq(jw) -# 0 for all wE lR+, the pair (F(jw), G(jw))
will always have generalized eigenvalues on oJ[]) whenever (T(P) n lR+ -# 0.
Indeed, this implies that det(E(jwk)) = 0 for some Wk E (T(p)nlR+. Accord-
ing to (2.34), we may conclude that a(jwk' z) will always have roots on or
symmetric about oJ[]). Since the roots of a(jw, z) are continuous functions
of w, there must exist also some WI such that a(jwI, z) has roots on oJ[]).
Thus, under this circumstance, we may assert without further computing
the generalized eigenvalues of F(jw) and G(jw), that the quasipolynomial
cannot be stable independent of delay.
The following result gives an alternative test in the same spirit. It may
be executed by computing sequentially the eigenvalues of two matrices. A
minor restriction, however, is that the assumption ao (jw) -# 0 needs to be
imposed for all wE lR+.
Theorem 2.12 Suppose that ao(jw) -# 0 for all wE lR+. Suppose also that
the quasipolynomial (2.6) is stable at T = O. Let (A, B, C) be a minimal
realization ofEl1(s)En-s), where El(S) and E2(S) are given by (2. 33}.
Furthermore, define the Hamiltonian matrix
._ (
H .- A BBT) lD>2nqx2nq
-CTC _AT E~ .
_aq_l(s)
Ma(s),~ (
-~)
-~
ao(s) ao(s) ao(s)
o o
~
1 o
Then, if (T (Ma(jWk)) n oJ[]) = 0 for all 0 -# Wk E lR+ such that jWk E (T(H),
T = 00. In these cases the quasipolynomial (2.6) is stable independent of
delay. Otherwise,
where jWk E (T(H), Wk E lR+, Wk -# 0, and Ok E [0, 27r] satisfy the relation
e- j1h E (T(Ma(jWk)).
det(jwI - A) ao(jw) -# 0,
det(jwI + AT) (-ltao(-jw) -# 0,
62 2. Systems with Commensurate Delays
det(jwI - H)
jWI-A -BBT)
det ( eTc jw + AT
det(jwI - A)det (I - BT(_jwI - A)-TeTe(jwI - A)-IB)
det(jwI + AT)
det(jwI - A)det (I - E2(jw)EIH(jw)Ell(jw)Ef(jw))
det(jwI + AT)
= (_1)n lao(jw)1 2 det (I - E2(jw)EIH(jw)Ell(jw)Ef(jw))
= (-It lao(jw)12 det (EIH (jw)) det (E{f - Ef (jw)Ell(jw)E2(jw))
1-
( _l)nl ao (.JW )1 2 ao 1 (.JW )1 2 dt(
e
El(jW)
Ef (jw)
E2(jW))
E{f (jw)
= (-It det(E(jw)).
(8 4 + 28 3 + 58 2 + 38 + 2) + (8 2 + 2) e- rs
+ (8 2 + 8 + 2) e- 2rs + (28 3 + 58) e- 3rs .
We may compute T based on Theorems 2.11 and 2.12. Construct first the
matrices P and H, both of which are 24 x 24 matrices. Using Theorem
2.11, we found the positive real eigenvalues of P to be
0.27440341977960 - O.OOOOOOOOOOOOOOj,
0.91548178048665 - O.OOOOOOOOOOOOOOj,
1.30259246725682 + O.OOOOOOOOOOOOOOj,
and the corresponding unitary roots are
-0.89898090347872 + 0.43798782537941j,
-0.97120708373425 - 0.23823685798887j,
0.88081464844625 - 0.47346124982147j,
respectively. The reader may want to verify that the modulus of all the roots
is 1.00000000000000. From these results, we found T = 0.37864202327266.
The computation results based on Theorem 2.12 turn out to be identical
up to the indicated precision (16 digits), and thus the two sets of results
demonstrate for one another the numerical accuracy of the tests. 0
Assume that the system is stable at T = O. Then by fixing z E allJJ, the sta-
bility problem amounts to determining whether the matrix (E:;;'=o AkZk)
64 2. Systems with Commensurate Delays
In this case we may compute the delay margin once finding the pairs
(Wi, ()i)' Consequently, we are led to the problem of finding all such criti-
cal Z E 81ID that u (2:;;'=0 AkZk) has an imaginary eigenvalue, and that of
solving the imaginary eigenvalues for all such Z E 81ID. The Kronecker sum
of matrices discussed in Appendix A comes to our aid in this task.
Recall that for square matrices A and B,
(2.35)
Bm - k := I ® Af, k= 1, "', m.
Furthermore, define
)
o I o
o I
-Bo -B2m-l
det(zU - V) = det(~Bkzk)
= zn 2 m det (f
k=l
Bm+kZk + Bm + f
k=l
Bm_k z - k ) .
det(zU - V)
This implies that for any z E 8lIJJ and any W E lR.+, det(jwI - L:;;'=o AkZk) =1=
O. Thus, for any T E [0, (0) and any w E lR.+,
det (jWI - f
k=O
Ake-jkTW) =1= O. (2.36)
~ ), V-(~
1 0 0
U= ( 0 1 0 1 00)
010
001 - 0 o 0 1 '
000 -1 1 101
which in turn yields IT(V, U) = {-I, O.5±O.866j}. Clearly, IT(V, U) CaD.
For Z = -1, IT (2:~=0 AkZk) = {O}. However, for Z = 0.5±0.866j = e±i'Tr/3,
IT (2::~=oAkZk) = {:rJ3 j}. Hence, (7r/3)/J3 = 7r/(3J3). This value
T=
is identical to that obtained in Example 2.1. Similarly, if we are to invoke
Theorem 2.11, we need to construct
0 0
P=' ( 1 0
J -1 0
-1 -1
2.5 Notes
T ~ 0, (3.2)
that is, the system contains only one delay. We shall first develop small gain
type, frequency-sweeping stability conditions for general LTI systems in the
form of (3.1). The results employ a matrix measure known as the struc-
tured singular value, a measure used widely in robust control analyis and
design. The structured singular value approach constitutes a step forward
in the study of systems with incommensurate delays, in that it provides a
necessary and sufficient condition by means of a well-developed theoretical
notion as well as a computational tool, for which standard numerical al-
gorithms and commercial software programs are available. It thus enables
a systematic analysis of the stability problem from a computational per-
spective, which can be accomplished at least approximately. On the other
hand, the exact computation of the structured singular value is also known
to be difficult, and in fact is known to be intractable in general. This raises
the speculation that the stability analysis for systems of incommensurate
70 3. Systems with Incommensurate Delays
(1: IIf(t)II~dt)
where
IIfll£2 := 1/2,
n ) 1/2
IIf(t)II2:= ( ~ Ih(t)12
3.2 Small gain/ J.£ theorem 71
i:
Define the Fourier transform of ! (t) by
I(jw):= !(t)e-jwtdt,
(3.3)
In turn, the function spaces L2 and £2 form an isometrical isomorphism.
The Fourier transform thus provides an equivalent characterization of a
signal via its frequency domain representation.
By so identifying the input and output signals with elements of function
spaces, a system can be viewed as an operator mapping the input space to
the output space, as shown in Figure 3.1. In particular, a linear system can
be viewed as a linear operator. Let H: L2 ---+ L2 be a linear system. Then
its induced norm can be defined as
In other words, the induced norm coincides with the peak magnitude of
the system's frequency response, whereas the magnitude is measured by
the largest singular value and has the interpretation as the system's gain.
It follows that an LTI system is BIBO L 2 -stable if and only if its frequency
response is essentially bounded on OC+. If in addition an LTI BIBO L 2 -
stable system is also causal, then its transfer function matrix H(s) will
be analytic in C+ and essentially bounded on 8C+. It is known that the
collection of all such functions forms another function space, known as the
Hardy space of Hoc functions, endowed with the Hoc norm and defined as
Lemma 3.1 (Small Gain Theorem) Suppose that M and tl are both
causal, linear, and BlBO L 2 -stable. Then the M-tlloop is BlBO L 2 -stable
tl
if
IIM~1I2,2 < 1. (3.4)
Hence, the small gain theorem provides a sufficient condition that governs
the stability of feedback interconnections. One should note that in its full
generality, modulo appropriate modifications, the small gain theorem in
fact is applicable in a much broader sense. Nevertheless, the statement given
in Lemma 3.1 is sufficiently general for our present purpose. In particular,
with a more far-reaching implication, it allows the study of time-varying
systems and can therefore be applied to systems containing time-varying
delays. This will be illustrated shortly via examples in the present chapter
and be investigated in greater depth in Chapter 8.
The small gain theorem has played an especially prominent role in the ro-
bust control of uncertain LTI systems. In robust stability analysis, one typi-
cally models a system's uncertainty by deterministic, unknown but bounded
perturbations to its parameters or transfer function matrix. A specific un-
certainty model, known as unstructured uncertainty, describes unknown
perturbations as stable transfer function matrices with a prescribed 11.00
norm bound, by means of the set description
(3.5)
l)
In this case, M is constructed so that
and the delay operators are collectively represented by the block diagonal
uncertainty
Each of the nonzero, diagonal blocks in .6. is a delay operator: .6. kx(t) =
x(t - rk). Since II.6. k I1 2 ,2 = 1, it follows that 11.6.112,2 = 1. Consequently, a
sufficient, delay-independent condition for the system to be stable is that
Ao is stable and 11£11100 < 1.
More generally, the idea of representing a delay by a fictitious uncertainty
can be generalized to systems with time-varying delays. We illustrate this
point briefly below; the idea will be developed in details in Chapter 8.
Consider the linear delay system
M
dx/dt I (1/s)1 I X
,J I I~--~--~I~
I
I
I
'-------I
I A0 I I
I
I
I I
I
J
L _________________ - - __ I
I I
I I
1~ x 2(t - r(t))dt ~ 1: 00
x2(u)du,
and the linear operator ~x(t) = x(t-r(t)) satisfies the inequality 11~112.2 ~
1. As a result, a sufficient condition for the system to be stable is that Ao
is stable and lI(sf - Ao)-t At 1100 < 1. 0
r- - - - - - - - -- - --- - - - - - - - - - - - ----I
1 M 1
1 1
1 X 1
1 1
1 1
1 1
1 1
1 1
1 _______________________ I
1 1
:___________
A 1 1
,/ P(s)
-
e- rs
The examples given above underscore the role that the small gain type
results may play in stability analysis of time-delay systems. They demon-
strate that the small gain theorem can be applied effectively to yield a
sufficient condition for delay-independent stability; sometimes, this is true
even for systems with time-varying delays. In particular, Example 3.2 shows
that in certain situations it is possible to obtain a necessary and sufficient
condition, specifically when the delay can be modeled as an unstructured
uncertainty modulo a side constraint at s = 0, as in the case of the Tsypkin
criterion. This, however, ceases to be true when the system contains multi-
ple delays. Example 3.1 exhibits that in the latter case, one unstructured
uncertainty can no longer represent fully several delays without undue con-
servatism. It thus calls for a more refined description that can incorporate
the structural information of the uncertainty that results when several such
3.2 Small gain/ Ji, theorem 77
delays are lumped together. For this purpose, the so-called structured un-
certainty provides a remedy.
Then the M -.60 loop is EIBO L 2 -stable for all Li E ~oo if and only if
(3.6)
where
D-MHDM2:0. (3.10)
(i) Ao is stable,
m
(ii) L Ak is stable, and
k=O
(iii)
f..Lx", (M(jw)) < 1, Vw > 0, (3.12)
where
Proof. The conditions (i) and (ii) are required to ensure that the system
be stable for rk = 00 and rk = 0, for all k = 1, "', m, and hence are
necessary. Under the condition (i), the system is stable independent of delay
if and only if for all s E C+
(3.14)
Let
~(s) := diag (e- r1S I, "', e-r",s 1) .
Since M(s) and ~(s) are analytic in C+, the spectral radius p(M(s)~(s)),
either in its own right or as a special case of f..Lx", (M(s)~(s)), defines a
subharmonic function in C+. As a result, it achieves its maximum on ac+.
That is,
sup p(M(s)~(s)) = supp(M(jw)~(jw)).
SEC+ wER
3.3 Frequency-sweeping conditions 81
Therefore, whenever (3.12) is true, the condition (3.14) holds for all s E «:\
except at s = O. However, under the condition (ii), it is necessary that
This implies that (3.14) also holds at s = 0, and hence it does for all s E C+.
In other words, the system is stable independent of delay. The proof for
the sufficiency part is completed. To show the necessity, assume first that
J-tx", (M(jwo)) = 1 for some wo > O. This means that a matrix U E Urn
exists such that
(i) Ao is stable;
(ii) J.Lxm (M(O)) < 1, or J.Lxm (M(O)) = 1 but det (I - M(O)) =J 0;
(iii) J.Lxm (M(jw)) < 1, Vw> O.
Consequently, while the small J.L theorem requires that J.Lxm (M(O)) < 1,
it is possible that J.Lxm (M(O)) ::; 1 while the system may still be stable
independent of delay.
It is then clear that in general the condition for delay-independent stabil-
ity does not match exactly that for robust stability; the latter is stronger
and in fact provides only a sufficient condition for stability independent
of delay. The reason, evidently, points to the fact that the delays cannot
be taken exactly as uncertainties, a phenomenon we have already observed
in Example 3.2. This difference helps explain a discrepancy between any
necessary and sufficient condition and the classical two-variable criterion,
which was perceived for some time as a necessary and sufficient condition
for stability independent of delay. To further illustrate, consider the case
of one delay. The following corollary is a counterpart to Theorem 2.1 and
follows immediately from Theorem 3.6.
(i) Ao is stable;
(ii) p (Ai) I AI) < 1, or p (Ai}l AI) = 1 but det(Ao + AI) =J 0;
(iii) P (jwI - AO)-I AI) < 1, Vw > O.
Suppose that Ao is stable. A sufficient condition for stability independent
of delay is then given by
Vw~O. (3.15)
Vw ~ 0,
for .6. = {811 8 E C}. It thus follows that the condition (3.15) is equivalent
to
det(1 - (sl - Ao)-I AIZ) =J 0,
which is equivalent to the standard two-variable criterion
In other words, the two-variable criterion, which coincides with a small J.L
condition, is only a sufficient condition. One notes that for such a condition
to be also necessary, there must exist an "uncertainty" ~(s) E .6. 00 with
which the matrix I - M(s)~(s) loses rank. Since T ~ 0 can be selected
3.3 Frequency-sweeping conditions 83
(i) Ao is stable;
m
(ii) L Ak is stable;
k=O
m
(2) Let ak(s), k = 0, 1, ... , m, be defined in (3.18) and (3.19). Then the
system (3.17) is stable independent of delay if and only if
"iw > O.
(i) a > 0;
q q q
Proof. The statement (1) follows from Lemma 3.4 (viii), which leads to
m
fLx,,,(M(jw)) =L IVk(jwI - AO)-lul·
k=l
To prove the statement (2), note that it constitutes a special case of (1),
with
o
Ao
1
-aO,n-l
)
,
Finally, the statement (3) is easily recognized as a special case of (2), with
q
L lakl
/1. x..{M(jw)) = ;:~ + a2
The proof is then completed. •
We conclude this section by presenting a small gain type sufficient condi-
tion for delay-dependent stability, using also the structured singular value.
Before proceeding, we note that in principle a delay-dependent stability
condition can be obtained in the process of computing /1.xm (M(jw)) , in
verifying (3.12). Alternatively, it may also be obtained by solving a maxi-
mization problem of the form described by (3.8). In either case, the solution,
were it attainable, would yield certain critical, "destabilizing" .6. E Urn, at
the frequencies w where /1.x,jM(jw)) = 1. The critical delay values may
then be estimated from .6. and w. Evidently, while the estimates would
provide an exact range of delay values for which the system is stable, the
computation of such critical .6. is, unfortunately, as difficult as that of
/1.xm (M(jw)) , if not more, and hence from a computational standpoint is
unlikely to be tractable in general. The condition to be developed below
inherits the same complexity in computing /1., but is more amenable to
analysis.
Consider the system (3.1) with incommensurate delays rk. Using
x(t - ri) = x(t) -l~ri x(u)du = x(t) -l~ri (~AkX(U - rk)) du,
we may rewrite the state-space equation (3.1) as
(3.20)
The process of transforming the original equation (3.1) to the new repre-
sentation (3.20) is known as model transformation, which is widely used in
both frequency and time domain stability analysis of time-delay systems.
It should be pointed out that although the stability of (3.20) implies that
of (3.1), the reverse is not necessarily true due to the presence of additional
dynamics. This subtle point will be illustrated in Chapter 5. The Laplace
transform of the integral term is given by
#0,
86 3. Systems with Incommensurate Delays
(3.21)
This consequently leads to a delay-dependent stability condition character-
ized by the structured singular value.
Proof. For any Tk > 0, k = 1, ... , m, we may first rewrite the condition
(3.21) as
or equivalently
where
that is, asymptotically, the computation cost will not grow faster than
at a polynomial rate, in the worst case and hence in all instances of the
problem. It is customary to refer to this class of problems as the class P.
For a problem not in the class P, there exists at least one instance for
which the required computation time cannot be bounded by a polynomial
function. Problems in the latter category are generally held intractable.
One should be cautious in interpreting the practical difficulty in computing
problems in P. The formalism gives no indication and indeed de-emphasizes
the importance of nand q. It could well be possible that for all practical
3.4 Computational complexity analysis 89
PcNP.
An analogy can be made here between finding a difficult proof for a mathe-
matical theorem versus checking an existing proof. A longstanding problem
in the complexity theory, one that remains unresolved today and seemingly
will be open in the foreseeable future, questions whether
P=lNP.
While no conclusion has been reached either way, it has been generally
believed that this is likely to be true, a de facto proposition supported by
ample empirical evidence. The fact, under the premise that it does hold,
commands fundamental importance in complexity studies, due to a strong
property of problems in NP. It turns out that NP contains a subclass in
which every problem can be transformed via a polynomial-time algorithm
to another problem in the same subclass. Such a problem is said to be N P-
complete. More specifically, we say that a problem P is NP-hard if every
problem in NP is polynomially reducible to P, and NP-complete if it is
NP-hard and P E NP. Here by the polynomial reducibility of a problem
PI to another problem P2, we mean that every instance PI of PI can be
transformed by a polynomial-time algorithm to an instance P2 of P2, so that
PI admits a solution if and only if P2 does, and PI does not admit a solution
if and only if P2 does not. The implication then is that it is unlikely to find a
polynomial-time algorithm for an NP-hard problem P, for otherwise every
problem in NP can be polynomially reduced to P and consequently solved
via a polynomial-time algorithm, but this would necessarily contradict the
fact that P =I NP. Similarly, since an NP-complete problem is also NP-
hard, NP-complete problems are unlikely to be tractable as well.
That every problem in NP can be polynomially reduced to an NP-hard
problem suggests a proof for NP-hardness. It follows that if a problem
PI is NP-hard and polynomially reducible to some problem P2, then P2
must also be NP-hard. In other words, if one is to prove that P2 is NP-
hard, one may attempt to reduce polynomially a known NP-hard problem
PI to P2· In particular, it suffices to polynomially reduce PI to just one
instance of P2, since that particular instance of P2 is as difficult as any
90 3. Systems with Incommensurate Delays
ai (1 0 0 1 1 0 o )T, i = 1, n,
i = 1, ... , n + 1. (3.26)
Suppose that the Knapsack problem admits a solution X E {-I, l}n for
some nonzero C E zn. Let Zo = 1, and
if Xi = 1
i = 1, ... , n.
if Xi =-1
Lemma 3.11 For a given positive definite matrix A E lR nxn , the following
problems are NP-hard.
Hence,
sup IZr AZ21 S; sup IIQZII~ = sup IZ H AZI.
Zl,Z2E1iii n ZE1iii n zE1iii n
Similarly,
sup IZr AZ21 = sup IZ H AZI.
Zl,Z2E1['n ZE1['n
We claim that
1 ).
Then in light of Lemma 3.4 (v), it is easy to observe that with kl = ... =
k2n = 1,
sup p2 (M~)
sup p2 (MU)
UEU2n
sup IZ HAZI.
ZETn
1
(3.28)
Hence, if the 'lI'-Knapsack problem has a solution, then sup IZH AZI
zEliiin
1. Otherwise, sup IZH AZI < 1. Consequently, the 'lI'-Knapsack problem
ZEliii n
has been polynomially reduced to a special instance of the sesquilinear
programming problem, and similarly the bilinear programming problem.
We then conclude that both problems are NP-hard . •
The above development finally culminates at the main result of this sec-
tion. We are now ready to construct one specific system with incommen-
surate delays, for which we show that it is stable independent of delay if
and only if the 'lI'-Knapsack problem admits a solution. The implication
then is that the 'lI'-Knapsack problem is polynomially reducible to a par-
ticular instance of the stability problem, thereby enabling us to conclude
the NP-hardness of the latter. Consequently, we conclude that the compu-
tational complexity in checking the stability will grow exponentially with
the number of incommensurate delays.
(ii) Is the system (3.1) stable for all rk E [0, Tk), where Tk > 0, k =
1, "', m?
94 3. Systems with Incommensurate Delays
Ao:=
CA-I 0f
-Wo
0
-1
0
wof
0
_A-l
W~Io ) ' (3.29)
0 -wof 0 -1
and
eke~+l k = 1, ... , n
e n+1 eI_n k = n + 1, ... , 2n
(3.30)
ek-n+lern+2 k = 2n + 1, ... , 3n
e2n+2eI_2n+l k = 3n + 1, ... , 4n
where ek E 1R2n+2 is the kth Euclidean coordinate. Denote m = 2n + 2,
m
and 4> .. := Ao + L: Akzk . It is clear that for any Zk E iiii, k = 1, ... , m,
k=l
4> .. =
(
_A-l
ZT
2 f
wof
o 0)
wof
-Wo _A-l Z3 '
o zI -1
for some Zi E F, i = 1, 2, 3, 4. Suppose that sup IZH AZI < 1. This
ZEliiR
implies that
where
Za = Zl +Z2
2
-n
ElI),
~ _ ( Wz woI)
z- -woI Wz .
det (jWOI - Ao - f
k=l
AkZk) = O.
By equating
Ok/WO k= 1, "', n
(211' - Ok)/WO k= n+ 1, "', 2n
Ok/WO k= 2n + 1, "', 3n
(211' - Ok)/WO k= 3n + 1, "', 4n
96 3. Systems with Incommensurate Delays
we have
det (jWOI - Ao - f
k=l
Ake-jrkWo) = O.
That is, with the matrices Ak in (3.29) and (3.30), the system (3.1) can-
not be stable independent of delay. To summarize, we have shown that
if the 'lI'-Knapsack problem admits a solution, then we may construct via
a polynomial-time algorithm the matrix A in (3.28), and in turn a delay
system (3.1) with Ak given by (3.29) and (3.30), so that it is stable inde-
pendent of delay. Conversely, if the 'lI'-Knapsack problem has no solution,
then the same delay system cannot be stable independent of delay. In other
words, we have reduced polynomially the NP-hard 'lI'-Knapsack problem to
a specific instance of the problem of determining whether the system (3.1)
is stable independent of delay. Consequently, we have proven Theorem 3.12
(i); that is, for systems with independent, incommensurate delays, the prob-
lem of checking delay-independent stability is NP-hard. Since the case (i) is
a special case of (ii), it follows that Theorem 3.12 (ii) is also true: the prob-
lem of checking delay-dependent stability, for any delay region bounded by
rk, is NP-hard as well . •
Theorem 3.13 Suppose that Ao is stable. Then the system (3.2) is stable
independent of delay if one of the following conditions holds.
(iv) For any induced matrix norm 11·11 and the corresponding matrix mea-
sure v(·), IIAIII < -v(Ao).
(v) For any induced matrix norm 11·11 such that IleAotll :::; (e-T/t, Vt ~ 0
for some ( ~ 1 and 7] > 0, IIAdl < 7]/(.
and
-v(jwI) + v(Ao) :::; v( -jwI + Ao) :::; v( -jwI) + v(Ao).
We claim that v(jwI) = v( -jwI) = 0, and hence v( -jwI + Ao) = I/(Ao).
Indeed, by definition, we have
. I)
v (JW = 1m
1· III + a"jwIIi - 1
= 1·1m 11 + j(Jwl- 1 =
0
.
<1-40+ (J <1-40+ (J
Similarly, 1/( -jwI) = O. We are thus led to the inequality
1
IIA111 < II(jwI - Ao)-lll'
the condition (iv) follows. Finally, to prove (v), we use the Laplace trans-
form
98 3. Systems with Incommensurate Delays
has an eigenvalue on the imaginary axis, including the origin. The relation
between the Hoo norm and the Hamiltonian matrix H is well known, and
the bisection method is well documented in the Hoo control literature. The
connection between the Hoo norm and other time-domain measures will be
further explored in Chapter 8.
The next result also draws upon frequency-sweeping tests, which give
conditions based on the solution of the Lyapunov equation
Theorem 3.14 Suppose that Ao is stable, and let P be the unique positive
definite solution of {3.31}. Then the system {3.2} is stable independent of
delay if one of the following conditions holds.
a (<I>(jw) Ad ~ a(PAd·
The condition (i) then follows from Theorem 3.13 (iii). To establish (ii),
let x E en be the eigenvector of <I>(jW)AI associated with the eigenvalue >.
such that 1>'1 = p (<I>(jw)A I ). Then, pre- and post-multiplication of (3.33)
by x H and x leads to
X
H ('XPAI + 2('XPAI)H) x = 1\12 H
A X X.
Since
X
H ('XP Al + ('XP At) H) X <
2 -
we have
However,
This gives
p(<I>(jw)At} = 1>'1 ~ p (,PAti T 2+ IPAII).
(3.34)
indicating that Theorem 3.14 (i) is less conservative than Theorem 3.13 (iv),
when the norm in question is II . 112. Indeed, let x E ]Rn be the eigenvector
100 3. Systems with Incommensurate Delays
_ T (Ai' + Ao) _ xT X
x 2 x - u(P)"
Since
x T (Ai' +
2 Ao) x <
_ \
Amax
(Ai' +
2 Ao) x T x -_ V2 (A) T
0 xx,
the inequality (3.34) follows. Thus, whenever Theorem 3.13 (iv) is true,
so is Theorem 3.14 (i). The example below demonstrates the merits and
weaknesses of these sufficient conditions.
Example 3.3 The following delay system is modified after Example 2.8:
where
AO=(~ ~ ~ ~ ),
-0.05 0.005 0.25 0)
A _ ( 0.005 0.005 0 0
1- 0 0 0 0 '
-2 -3 -5 -2 -1 0 -0.5 0
Figure 3.7 plots the bound lI(jwI - AO)-l AlII, for Holder f1 f2, and foo
induced matrix norms, respectively. The reader may compare these plots
to assess the conservatism of Theorem 3.13 (ii-iii), with the three different
norms in use. A simple calculation shows that in all three cases, v(Ao) > O.
Thus, Theorem 3.13 (iv) ceases to be useful. The estimates of/3 based upon
other sufficient conditions are summarized in the following table. 0
0.9
0.8
0.7
0.6
~0.5
.a
'6,
~04
0.3
02
O.t
0
0 0.5 t.5 2.5 3.5 4.5
Frequency
_H,I
__ H,
t8 .... H.
t.6
II>
'Ot.2
.a
'2
Cl
~t
0.4
o
-bn - I
)
The conditions in Theorem 3.13 and Theorem 3.14 can be further simplified
for this class of systems. The following table lists the relevant quantities,
Al cb'l ' B
p(<P(S)AI) IbT<p(s)CI 1 bn_ls n +"+bo 'I
Sn+an_IS n I+".+ao
a(<p(s)At} 11<p( s)c11211b112 Isn+a n ISn
')'{s)
I+"'+aol
IIAIIIt m~ Ibil E ICil m~lbil
• i •
a(PA I ) IIPcI1211bl12 {3
p (IPAII'l~+IPAII) IW'IPcl+IIPcIl21IbI12 a+{3
2 2
n-I
where <p(s) := (s1 - A)-I, P = [Pij], a:= E IbiPinl, and
i=O
{3 .-
.- (~lbiI2) (~IPinI2) ,
b = (-b o -b i
This can be done in the same manner as in the proof for Theorem 3.14 (ii).
Accordingly, it follows that for any U E Um,
X
H (>.RU + 2(5:. RU t) x= I>. 12 x Hx
and
1.11',; p CRU + fRU)H) .
Since lUI = [IUij 11 = I and
The following conditions are then clear. Note that by setting dk = 1/ y'm,
k = 1, ... , m, these conditions replicate some of those given in Theorem
3.15.
Vw?: o.
3.5 Sufficient stability conditions 105
{i} For any absolute or unitarily invariant induced matrix norm II . II,
/lM(jw; T)II < 1, Yw :2: O.
{iii} 0' ([ TIAIG ... TmAmG]) < -V2 (~Ak) /vm(m + 1).
(iv) If for ,orne ( 20 1 and ~ > 0, 0' [exp (t. A,t)1 (e~"', Vt
<; 20 0,
106 3. Systems with Incommensurate Delays
(vi)
and
The proof of this result resembles that for Theorem 3.15 and hence is
omitted.
rk ::::: O. (3.35)
k=l k=l
We note that the two sides of (3.35) need not contain an equal number of de-
lays, but this can nevertheless be assumed with no loss of generality. More-
over, as in retarded systems, the delays can be commensurate or incom-
mensurate. For neutral systems of commensurate delays, the differential-
difference equation is specified to
m m
We show in this section how the preceding results may be extended to neu-
tral delay systems, with both commensurate and incommensurate delays.
Stability of neutral delay systems proves to be a more complex issue
because the system involves the derivative of the delayed state. In fact, the
stability notion itself needs to be modified in the first place. Define the
3.6 Neutral delay systems 107
characteristic function by
(3.37)
Based on Theorem 1.6 in Chapter 1, we give the following stability defini-
tion.
For the system (3.35) to be stable, it is necessary that the difference equa-
tion (3.39) be stable, whereas the stability is defined in the same sense as
in Definition 3.2.
(3.40)
108 3. Systems with Incommensurate Delays
°
Lemma 3.18 Let rk = kT, k = 1, "', m. Then the system (3.39) is
stable for all T ~ if and only if p(B) < 1, where
Bm - 1
°
I
°
k=l
Hence, if p(B) < 1, we may find an 0: < such that (3.40) holds. Otherwise,
assume that p(B) ~ 1. Then, B has an eigenvalue >..(B) = p(B)e j8 for
some B E [0, 211"]. Consequently, for any T > 0, one may find s = 0: +
jw E C+, where 0: = (logp(B)IT) ~
det (l - Be- TS ) = 0, and hence
°
and w = (BIT) ~ 0, such that
det (I -f
k=l
Bke- hS ) = 0.
Therefore, the system (3.39) is not stable. This completes the proof. _
In much the same manner, we may then obtain similar delay-dependent
and delay-independent stability conditions. With no loss of generality, we
shall assume that the system (3.36) is stable at T = 0. This enforces the
condition that
det(I - ~B}' 0
The following results resemble the frequency-sweeping criteria developed
for retarded systems, presented in Theorems 2.5 and 2.6, respectively. The
proofs are also similar and are left to the reader. We remind the reader of
the definition of delay margin, which can be generalized to neutral delay
systems in the same manner.
3.6 Neutral delay systems 109
Theorem 3.19 Suppose that the system (3.36) is stable at T = 0, and that
p(B) < 1. Let rank (Am) = q. Furthermore, define
where
1
G(s) :=
H(s) :=
( -(SI~~) °
Fl(S)
diag (1, ... , 1, - Fm(s» ,
Fk(S) := k=l, ... , m.
Then,
T:= min Ti.
l~i~q+n(m-l)
System (3.36) is stable for all T E [0, T), but becomes unstable at T = T.
(iii) (1 - f k=l
Bk) -1 f
k=O
Ak is stable, and
(iv)
p(Mm(jw)) < 1, Vw > 0, (3.41)
where
~(s).rl(S)
Mm(s) := ( .
Ai> with retarded systems, stability conditions in this case can be simplified.
Define the polynomials
q
b(z) := Lbkz k
k=l
bks n + ak(s), k = 1, ... , q,
where ao(s) and ak(s), k = 1, ... , q, are defined as in (3.18) and (3.19).
Based on Theorems 3.19 and 3.20, it is not difficult to extend Corollaries
2.7 and 2.8 to the neutral system (3.42).
(i) b(z) is Schur stable; that is, b(z) has all its zeros in the open unit
disc;
where
_E.!..W. ...
(
ao(s)
1 ... -~
ao(s)
o -~)
ao(s)
0
Ma(s) := b ·· ..
· .
1 0
where
1 0
Ga(s) :=
(-aL 0
-Pl(S)
1
-Pq-l(S)
).
Ha(s) := diag (1, ... , 1, Pq(s)).
3.6 Neutral delay systems 111
Then,
System {3.42} is stable for all 7 E [0, r), but becomes unstable at 7 = r.
B = (-~.8 g)
and
k)
o .
s+2
For (3 = 0.2625, p(B) = 1.05. Under this circumstance the system is unsta-
ble for any 7 2': O. For (3 = 0.09, p(B) = 0.9. Figure 3.8 plots the spectral
radius p(Ma(jw)), which enables us to conclude that the system is stable
independent of delay. 0
0.75
0.65
0.6
0.55
0.5'--~-~-~----'-~---'--~~-~---.J
o w g • ~ ~ ~ ~ • w
Frequency
Finally, extensions may also be pursued for neutral systems with in-
commensurate delays. In particular, by mimicking the proof of Theorem
3.5, we obtain the following necessary and sufficient condition for delay-
independent stability.
112 3. Systems with Incommensurate Delays
{iii} (1 - ~ k=l
Bk) -1 ~ Ak is stable,
k=O
and
where
B .-
~·-
M(s) :=
Vw > 0.
3.7 Summary 113
3.7 Summary
Our development of the frequency domain tests in this and the preceding
chapter was heavily influenced by robust control theory. Whether for sys-
tems with commensurate delays or those with incommensurate delays, the
stability problem has been interpreted as one of robust stability. Indeed,
from a conceptual standpoint, delay-independent and delay-dependent sta-
bility can both be held as notions of robust stability, whereas the delay
values are uncertain within the specified intervals. The application of the
continuity property of the stability exponent is reminiscent of the well-
known zero exclusion principle, which has played a critical role in robust
stability analysis and will continue to be an enabling tool in Chapter 4 in
our study of robust stability of time-delay systems. The use of the struc-
tured singular values as stability criteria, as well as the use of the spectral
radius for systems of commensurate delays, thus come more as a necessary
outcome than a sheer coincidence. More broadly, both can be interpreted
as generalized notions of gain, and the stability criteria as small gain con-
ditions. Small gain type results can be suitably applied to time-varying
systems; this has been demonstrated in the present chapter using examples
involving time-varying delays and will be further discussed in Chapter 8.
They can also be extended to address uncertain time-delay systems and
used as a synthesis tool in control design for time-delay systems.
While it appears justified to contend that the stability problem is largely
resolved in the case of commensurate delays, the problem remains open
for systems with incommensurate delays. The proof of NP-hardness in this
case leads to, on a firm theoretical basis, the somewhat pessimistic con-
clusion that it is unlikely to be tractable in general. One should therefore
resort to computable, albeit approximate, stability tests. The structured
singular values provide yet again relatively less conservative sufficient con-
ditions, by means of its bounds. Various relaxations and derivatives, with
varying degrees of conservatism versus computational ease, can be obtained
from the bounds, among which the least conservative is the D-scaled up-
per bound. This bound can be computed efficiently using commercial soft-
ware programs and yet still provide an acceptable answer. Amid the gener-
ally pessimistic prospect, on the other hand, it should be understood that
114 3. Systems with Incommensurate Delays
3.8 Notes
3.8.1 Small gain theorem and J-l
Robust control is an immensely rich subject of deep mathematical sophisti-
cation. There have been many well-written texts on the subject. The recent
survey article by Chen and Tits [36] provides a concise summary of key
progress in robust stability analysis. The structured singular value is cred-
ited to Doyle [57] and Safonov [238] and has found widespread applications
in robust control design. An in-depth study of this notion can be found in
Packard and Doyle [219], and Zhou, Doyle, and Glover [294]. A generalized
J.L with both complex- and real-valued uncertainties, suggested by Fan, Tits,
and Doyle [62], can be useful in studying robust stability of time-delay sys-
tems. The small gain theorem is credited to Zames [290, 291], whose more
recent treatise can be found in the texts by Desoer and Vidyasagar [54],
Zhou, Doyle, and Glover [294], and Chen and Gu [33]. The same references
contain detailed discussions on norms and spaces relevant to the control
context.
which all studied retarded as well as neutral delay systems. In fact, it was
due to [49] and [42] that the continuity property was firmly established.
Hale, Infante, and Tsen [102] investigated in meticulous detail the stability
problem with the general state-space formulation and gave the sweeping
conditions for scalar differential-difference equations. The approach via the
structured singular value, and the suggestion that the stability problem
be treated as one of robust stability, were initiated by Chen and Latch-
man [35], which led to the f-L characterization of conditions for stability
independent of delay. The idea was later extended in development of delay-
dependent conditions by Niculescu and Chen [208] and Huang and Zhou
[119]. The results in Section 3.3 are mainly based on [35, 208] and [208].
Chen and Latchman [35] and Huang and Zhou [118, 119] also addressed
robust stability and stabilization of delay systems in this framework.
(4.2)
of the system is less than -E. In this case f(s) is called Hurwitz stable. In
(4.2) coefficients aki are real or complex numbers and exponent coefficients
118 4. Robust Stability Analysis
are real and ordered as follows 0 = rO < rl < ... < rm. Function /(8) is
known as a qua8ipolynomial and may be written in different forms
m n
/(8) = LPi(8)e- riS
=L 1/Jk(8)8 n -k, (4.3)
i=O
where polynomials
and quasipolynomial
k = 0, 1, ... ,n.
Example 4.1 Let us consider the scalar equation
/(8) = 8 + 8e- rs + 1.
All zeros of the quasipolynomial have negative real part. To check this we
assume first that it has a zero 80 = a + j/3, with positive real part, a > O.
Then the real and the imaginary parts of /(80) are equal to zero
/3sin(r/3) + 1 = 0, (4.5)
/3[1 + cos(r/3)] O. (4.6)
cos(r/3) = -1.
This implies that sin(r/3) = 0, which contradicts (4.5). In other words /(8)
has no zeros in the closed right half-plane of the complex plane. Therefore
4.3 Zeros of a quasipolynomial 119
all zeros of the quasipolynomial have negative real part. On the other hand
the equation (4.4) is not exponentially stable because
sup {Res I f(s) = O} = O.
This statement can be easily checked if we observe that if So is a zero of
f(s) then
Remark 4.1 There are also examples of neutml systems where all the ze-
ros are on the open left half complex plane, but the system is not stable in
the sense of Lyapunov in addition to not exponentially stable.
where j is the imaginary unit and aki, bki , ai, /3i are real numbers. This
function may be written in two forms
m n
f(s) = LPi(s)e(Oi+if3i)8 = L 1Pk(S)sn-k,
i=O k=O
where
n
Pi(S) = L(aki + jbki)Sn-k, i = 0,1, ... , m,
k=O
and
m
1Pk(s) = L(aki + jbkde(Oi+if3i)S, k = 0, 1, ... , n.
i=O
In the following we assume that the
120 4. Robust Stability Analysis
• exponent coefficients (ai + j{3i), i = 0,1, ... ,m, are distinct complex
numbers;
• polynomials Pi (8), i = 0, 1, ... , m, are not trivial.
Under these conditions f( 8) may have a finite number of zeros only in the
case when m = 0. In the following we assume that m > 0, if not explicitly
stated otherwise.
Then a sufficiently small E > 0 may be chosen such that the E-sectors
ri(E) = {s=PeiI"IPE[O,OO)''PE['Pi-~E''Pi+~E]}'
i=1,2, ... ,N
Theorem 4.1 For every E > 0 there is R(E) > 0, such that all zeros of
f(s) with magnitudes greater than R(E) lie inside one of the E-sectors ri(E),
i=1,2, ... ,N.
Re {[(ai - ao) + j({3i - {30)] s} ::::: -p lsi sin G) , for i = 1,2, ... , m. (4.8)
Now we select R > 0 such that all zeros of all polynomials Pi(S), i =
0,1, ... , m, lie inside the disc
DR = { s I lsi < R } .
Hence there exists R(c) ~ R such that for all s E <I>(c) with magnitude
greater than R(c) the expression in square brackets in (4.9) is not zero.
The first two factors in (4.9) are also nonzero for such values of s. It means
that f(s) has no zeros in <I>(c) with magnitude greater than R(c). Repeating
these arguments for all other vertex points of the exponential diagram we
arrive to the theorem statement. _
From the above theorem one can easily derive the following consequence
about zeros of f(s).
Proof. It follows from geometric consideration that in this case one of the
rays points toward the open right half complex plane. So the corresponding
c-sector lies entirely in this half plane. According to Remark 4.2, inside the
sector f(s) has an infinite number of zeros. _
This corollary implies that the only case when f(s) may have all zeros
in the open left half complex plane is when all exponent coefficients of the
quasipolynomial have the same imaginary part. In this case, without any
loss of generality one may assume that all exponent coefficients are real
numbers.
= L cv s kv el'v s ,
N
f(s) (4.11)
v=O
where all complex coefficients Cv are supposed to be nonzero, and addi-
tionally we assume that there are no two terms in the sum for which both
kv and Iv are the same; that is, the terms have either different potential
factors or different exponential factors.
4.3 Zeros of a quasipolynomial 123
Let us superpose the b, k)-plane with the complex plane in such a way
that the ,-axis coincides with the real axis of the complex plane and in
turn the k-axis coincides with the imaginary one. On this joint plane we
mark points bv, kv) corresponding to all terms in (4.11).
Then we construct the upper part of the envelope of these points (see
Figure 4.2).
Im(k)
Re(y)
This upper part is known as the potential diagram of f(s), and it consists
of a finite number of segments. Let M be the number of these segments.
With each segment of the diagram we associate the logarithmic curve
AK = {s=x+iylx=ILKln(y), YE[l,oo)},
K=1,2, ... ,M.
Here ILK is a real number such that vector (ILK' 1) is along the direction of
the outer normal to the corresponding segment. If ILK > 0, then AK lies
in the right half complex plane, while the curve corresponding to ILK < 0
belongs to the left half complex plane. For ILk = 0 the corresponding AK
coincides with the imaginary axis.
For sufficiently small E > 0 the logarithmic E-sectors
_ +' I x E [(ILK -
{ s-x !E) In(y), (ILK + !E) In(y)], }
zy YE[l,oo) ,
K=1,2, ... ,M
have no common points except (0,1).
The following theorem describes distribution of zeros of quasipolynomial
(4.11).
Theorem 4.3 For every E > 0 there exists R1(E) > 0, such that all zeros
of f (s) in the upper half complex plane with magnitudes greater than Rl (E)
lie in the union of logarithmic E-sectors AK(E), K = 1,2, ... , M.
124 4. Robust Stability Analysis
s = ILln(y) + jy (4.12)
for some IL E (ILk' ILk+1)' and y > 1. Here ILk and ILk+l are coefficients ofthe
corresponding logarithmic curves Ak and Ak +1 (ILk < ILk+l)' If c: is positive
and
< ILk+l - ILk
c: 2'
then logarithmic sectors Ak (c:) Ak+ 1 (c:) have no common points except
(0,1). In this case we denote by w(c:) the part of the upper half com-
plex plane between the sectors. In the set w(c:) the quasipolynomial (4.11)
can be factorized as follows
Then simple geometric manipulations show that for all IL E [ILk + c:, ILk - c:]
the following inequality holds
(4.14)
Now we can see that along the logarithmic line s = ILln(y) + jy, y ~ 1,
where IL E [ILk + c:, ILk - c:J,
The right hand side of the last inequality approaches zero as y --4 00.
Therefore the sum in the square brackets in (4.13) also approaches zero as
4.4 Uncertain quasipolynomial 125
y -+ 00. In other words, for sufficiently large values of y, function I(s) has
no zeros along any logarithmic curve (4.12) where IL E [ILk -c:, ILk +c:J. This
observation proves the fact that there exists R > 0 such that no zero of I(s)
with magnitude greater than R lies in the domain 'If (c:). These arguments
can be applied for all other interior vertices of the potential diagram.
A slightly modified argument can also be applied for the two extreme
vertices of the diagram. _
Remark 4.3 The zeros oll(s) with large magnitudes in the lower hall
complex plane lie in the union 01 logarithmic c:-sectors obtained by the mir-
ror image 01 AI\: (c:) , K. = 1,2, ... ,M, with respect to the real axis.
Remark 4.4 Applying the principle 01 argument it may be shown that I(s)
has infinite (countable) number 01 zeros in every logarithmic c:-sector AI\:(C:)'
It follows from the theorem that I(s) has zeros with arbitrarily large
positive real parts when at least one of the values ILl\: is positive. Using the
potential diagram of I(s), one can conclude that such positive ILl\: exists
only if the outer normal of one of the segments forming the diagram points
toward the open right half complex plane. In order to guarantee the absence
of zeros of I(s) with arbitrary large positive real parts, one has to assume
that one of the terms in (4.11), for example, Cos ko e"YoB , satisfies the following
two conditions
Corollary 4.4 Quasipolynomial (4.11) may have all zeros in the open left
hall complex plane only il it has a principal term.
Without any loss of generality one may assume that aoo =I 0, then aoosne- rOB
With quasi polynomial (4.2) one may associate two vectors, the coefficient
vector
QF = { (a,r) I !(s,a,r) E F } C L
F = { !(s,a,r) I (a,r) E QF }
Remark 4.5 The assumptions Al-A3 are necessary conditions for the
robust stability of F.
Proof. The necessity part is nearly obvious. In fact, the first condition
of the theorem holds because all quasipolynomials in F are Hurwitz stable.
To check the second condition let us assume that VF(jwo) contains the
origin. There exists a quasipolynomial in the family that has jwo as a zero
and therefore the quasipolynomial is not Hurwitz stable. This contradicts
the Hurwitz stability of all members of F.
Now let us address the sufficiency part. It follows directly from assump-
tion A3 that there exists a positive value Ro such that no one member of
F has zeros in the complex half-plane
I
{ s Re( s) ~ - ~c }
128 4. Robust Stability Analysis
with magnitude greater or equal than Ro. According to the second condition
of the theorem there are no quasipolynomials in F that have zeros on the
imaginary axis.
Assume by contradiction that there exists f(s, 80, ro) in F that is not
Hurwitz stable. It means that the quasipolynomial has at least one zero
So with positive real part. By the first condition of theorem, there is
f(s, a1. rl) E F which is Hurwitz stable, that is, all zeros of f(s, a1. rl)
have negative real parts. According to assumption A4 there is a path in
QF:
a = aI" r = r JL , 11- E [0,1) ,
connecting points (ao, ro) and (a 1. rJ). Let us define
11-0 = inf {11- > 0 Isuch that f(s, aI" rJL) is Hurwitz stable}.
It follows directly from the definition of the value that 11-0 E (0,1). Poly-
nomial f(s, a JLo ' r JLo ) cannot be Hurwitz stable; otherwise by a continuity
argument, assumptions A1-A3 imply that there should be a positive II
such that for every 11- E (11-0 -II, 11-0 +11), quasipolynomial f(s, aI" rJL) is also
Hurwitz stable, which contradicts the definition of 11-0' On the other hand,
f(s,aJLo,r JLo ) cannot have a zero with positive real part because otherwise
for 11- a little bit smaller than 11-0, quasipolynomial f(s, aI" rJL) will also have
a zero in the open right half complex plane which contradicts again to the
choice of 11-0' It also cannot have zeros on the imaginary axis by the sec-
ond condition of the theorem. So f (s, a!Joo' r!Joo) cannot be unstable, either.
This contradiction has its origin in our assumption that f(s, 80, ro) is not
a Hurwitz stable quasipolynomial. This conclusion completes the proof of
the sufficiency part of the theorem. _
The second condition of the theorem admits certain relaxation.
• at least for one point, jwo, of the imaginary axis 0 fj. VF(jWo)i
• 0 fj. aVF(jw) for all other points of the imaginary axis. Here aVF(jw)
stands for the boundary of the value set VF (jw).
Proof. The statement follows directly from the property that the value
set changes continuously with respect to w. _
It is a difficult task in some applications to check the second condi-
tion of the theorem because there are no effective algorithms for the exact
construction of the value set. Also in many cases one cannot check this con-
dition for all points on the imaginary axis but only for a finite set of such
points. Having in mind this observation it seems quite natural to look for
4.5 Edge Theorem 129
L L a~~} 8n-ke-ris,
n m
f(v}(8) = V = 1,2, ... , N,
k=O i=O
of quasipolynomials, where as before, 0 = TO < T1 < ... < Tm. The poly-
topic family is defined as the convex hull of the quasi polynomials
A= {a t =
v=1
I ~ 0, v = 1,2, ... , Nj
JLva(v} JL v t
v=l
JLv = 1} ,
N
aOO = "'"
~ (v)
J.Lvaoo > 0.
v=l
Condition 4.2 There exist c > 0 and R > 0 such that members of 1}! have
no zeros in the half-plane
{s IRe( s) 2:: -c }
Remark 4.8 Condition 4.2 holds automatically when for all generators
deg(p~v») = n > deg(p~v»), i = 1,2, ... , m.
Remark 4.9 Let rI, r2, ... , rm be such numbers that ri = kir, where all ki'
i = 1,2, ... , m, are natural numbers. Then Condition 4.2 holds if and only
if all polynomials from
4.5 Edge Theorem 131
are anti-Schur, that is, all zeros of every polynomial from II lie outside the
closed unit disc of the complex plane. Here
m
1f (v) ( Z ) -_ '~
"' (v) Z ki , v -- 1 , 2 , ... , N .
a Oi
i=O
For the rest of this section a key role is played by the edges of A. An
edge is a segment connecting two vertices of A such that the set A without
the segment remains convex. In other words, an edge is a one-dimensional
face of polytope A. Every edge can be written as the convex hull of two
vertices
JLa(a) + (1 - JL)a(f3), where JL E [0, 1J .
Theorem 4.6 Let family P satisfy Conditions 4.1 and 4.2. All members
of the family are Hurwitz stable if and only if all quasipolynomials of the
edge subfamilies from Ep are Hurwitz stable.
°
unstable quasipolynomial II (s) E P. Then it should have a zero, say Sl,
with Re(sl) 2:: 0. It means that E Vp(Sl)' Observe that
for every f(s) E P. Hence, there exists an So with sufficiently large real part
such that no member of P has So as a zero; in other words, 0 ~ Vp(so).
Define now s-y = (1 - ,)so + ,SI and compute the value set Vp(s-y). For
, = 0 this set does not contain the origin, while for , = 1 it does. By
continuity arguments there exists '0
E (0,1], such that 0 E aVp(s-yo)'
It is clear that Re(s-yo) 2: O. According to our previous observation there
should be an edge family (4.16) such that 0 = ILof(e» (s-y) + (1- ILo)f({3) (s-y)
for some ILo E [0,1]. In other words, the edge quasipolynomial fo(s) =
ILof(e»(s) + (1 - ILo)f({3)(s) is not stable. This conclusion contradicts the
second condition of the theorem. The contradiction is a direct consequence
of our assumption that family P has at least one unstable quasipolynomial.
Theorem 4.6 reduces the robust stability analysis of a poly topic family to
•
a simpler problem of Hurwitz stability analysis of a finite set of one parame-
ter subfamilies of the type (4.16). Subfamilies (4.16) serve as a testing set
for the robust stability analysis of P. The number of such subfamilies may
in general be very high; the upper bound for this number is C'j. = N(~-I).
Analyzing the proof one can observe that not all edge subfamilies need to
be taken into account in the analysis. We only need to check those edges
whose images appear on the boundary of the value set for at least one point
jw. The number of such subfamilies may be significantly smaller than that
of all edge subfamilies. This observation sometimes allows one to obtain a
simpler testing set. But one question remains open: How can one check the
stability of an edge subfamily?
Claim 4.1 Let generator quasipolynomials fo (s) and h (s) be Hurwitz sta-
ble. Then all quasipolynomials f/1-(s), IL E [0,1] are Hurwitz stable if and
only if the complex curve
fo(jw)
z = -(--)'
h jw wE (-00, +00),
does not touch the negative real semi-axis of the complex plane.
4.5 Edge Theorem 133
Proof. The statement follows directly from the fact that the instability
of one of the quasipolynomials from (4.17) means that in (4.17) there is a
quasipolynomial with at least one zero on the imaginary axis. •
One interesting question is the following: Under what conditions does the
Hurwitz stability of vertex quasipolynomials fo(s) and h(s) imply that of
f/.L(s) for all J.L E [0, I]? This question leads us to the concept of convex
directions. First we rewrite (4.17) as
The next theorem gives necessary and sufficient conditions for g(s) to be
a convex direction.
1. 2r:S rm;
2. for all w > 0 such that q(jw) -:J. 0, the following inequality holds
o
134 4. Robust Stability Analysis
where
p~l)(S)
p~2)(s)
pP)(s)
p~4)(s)
These polynomials are such that the value set VIIj (jw) is a rectangle whose
corner points are p~II)(jw), lJ = 1,2,3,4:
The value set of IIie- ris at the point jw can be obtained by a rotation of
VIIi (jw) about the origin through the angle -riW. This observation allows
us to conclude that the value set V1(jw) of the interval quasipolynomial I
is a sum of such rotated rectangles:
• there are E > 0 and R > 0 such that in the interval family
Proof. This follows directly from the Edge Theorem and convex direc-
tion concept. •
(4.20)
136 4. Robust Stability Analysis
Proof. In fact, if f(s) has a zero, so, for which Re(so) ~ 0, then p(s, z)
has zero (so, zo) = (so, e- T1SO ) E ~ •
It is convenient to introduce the following change of variables to achieve
a more uniform formulation:
1-z
Sl = 8, and 82 = --.
l+z
Let n2 be the partial degree of p(8, z) with respect to z, then polynomial
Here 81,82, ... ,8 m are independent variables, and coefficients ail i2 ... i m are
real or complex constants.
Two polynomials are called relatively prime if they have no common
divisors except constant polynomials. Two primes are either coinciding up
to a nonzero constant multiplier, or are relatively prime.
4.6 Multivariate polynomial approach 137
A vector (S~o), S~o), ... , 8~)) with complex elements is called a zero of
(4. 21) , 1'f P((0) (0)) -- 0 •
8 1 ,8 2(0) , ... , 8 m
(4.22)
The main coefficient a~':) (S1, S2, ... , sm-d is a nonzero polynomial with
respect to the rest of the variables. Therefore, one can always fix these
variables such that the main coefficient is nonzero to arrive at an nm th
order polynomial of single variable Sm. This polynomial has exactly nm
roots. Each such root, together with the fixed variables, defines a zero of
p(S1, S2, ... , sm). Roughly speaking, there are a finite number of (m - 1)-
dimensional zero manifolds in the m-dimensional complex space.
One more essential difference between univariate and multivariate poly-
nomials is the following: Two multivariate polynomials may be relatively
prime and still possess common zeros.
(4.23)
With every polynomial P(Sl, 82, ... , Sm) from the set, one may associate
the real coefficient vector
from the N-dimensional coefficient space, where N = (nl + 1)(n2+ 1) ... (nm
+ 1).
In the following we will often denote the vector (s 1 , 82, ... , 8 m ) as s. The
following lemma is the standard extension of the well-known continuity
property of zeros of univariate polynomials with respect to small coefficient
variations.
Lemma 4.10 Let polynomial Po(s) E Pnl,n2, ... ,n", have a zero s(O) =
(s~O), s~O), ... , s~)). Denote by a(O) the coefficient vector of thi8 polynomial.
138 4. Robust Stability Analysis
For arbitrary c > 0 there exists 8 > 0 such that every polynomial p(s) with
the coefficient vector a in the 8-neighborhood of a(O):
(4.24)
Definition 4.2 A polynomial p(s) E Pnl,n2, ... ,n", is said to be strict sense
stable (SSS) if
p(s) =1=0, "is E r~).
Definition 4.3
2. for m > 1: given vector (n1, n2, ... , n m ), polynomial p( S1 , S2, ... , sm)
of degree (n1, n2, ... , n m ) is called stable if it satisfies the following
conditions:
The principal advantage of the stability concept is that the class of stable
polynomials is the biggest one for which polynomials preserve the stability
property under small coefficient variations.
Lemma 4.11 Let a polynomial p(s) E Pnl,n2, ... ,n", be stable, and assume
that nm > O. Define the transformed polynomial
Theorem 4.13 Let P(S1, S2, ... , sm) E P nl ,n2, ... ,n m be a stable polynomial.
There exists a E: > 0 such that every polynomial with a coefficient vector
from the E:-neighborhood of the coefficient vector ofp(s}, S2, ... , sm) is stable.
Theorem 4.14 Let p(s}, S2, ... , sm) E P nl ,n2, ... ,n m be a stable polynomial.
If nm > 0, then derivative
(1)( )
8 p( Sl , S2, ... , sm) nl g
'""' v a k S2, ... , sm k_
8s m L.t 8s Sl - ...
k=O m
nm-l (m-1)(
g )
'""' v a i Sl, ... , sm i
L.t g sm-1
i=O vSm
(4.26)
Corollary 4.15 Let p(s}, S2, ... , sm) E P nl ,n2, ... ,n m be a stable polynomial.
Then all coefficients in decomposition (4.22) are stable (m - I)-variate
polynomials from P nl ,n2 , ... ,nm-l .
140 4. Robust Stability Analysis
As it follows from the example below, for SSS polynomials this statement
is not true in general.
It was already mentioned that the main reason for the enormous sensitiv-
ity of the strict sense stability to small coefficient variations is the existence
of zeros close to the essential boundary (4.25) of the set r~). The following
theorem shows that stable polynomials have no such zeros.
Theorem 4.17 For every stable polynomial p(s) E Pnl,n2, ... ,n", there ex-
ists c > 0, such that it has no zeros in the c-neighborhood of the essential
boundary n.
Of course we assume that I c P n1 ,n2, ... ,n"" that is, all members of the
family have the same degree n = (nl, n2, ... , n m ). Moreover we assume that
all lower bounds are positive, that is Q:il ... i", > 0.
This set consists of 2m elements. With every vector from the set, one can
associate four polynomials as follows:
4.6 Multivariate polynomial approach 141
Given a sign vector (vt, V2, ... , v m ), the coefficients of the first polyno-
mial p~Vl.V2'''''V''') (SI, S2, ... , sm) are defined by the following rule: When the
index sum is even (il + i2 + ... + im = 2l) then
are defined by the rule: When the index sum is even (il + i2 + ... + im = 2l),
then
are defined by the rule: When the index sum is even (il +i2 + ... +i m = 2l),
then
are defined by the rule: When the index sum is even (il + i2 + ... + im = 2l)
then
142 4. Robust Stability Analysis
Proceeding in a similar way for all the elements from 3, one finally ob-
tains 4 . 2m polynomials. Taking into account that all polynomials in I
have only real coefficients, we may reduce the set by half and consider only
those of them for which Vm = 1. After this reduction, we obtain the desired
testing set of 4· 2m - 1 polynomials which is referred to in Theorem 4.18.
where a(O) denotes the center point and p(O) (s) = p(s, a(O») is the center
polynomial; r > 0 is the radius of the diamond. From here on we denote a
diamond family of multivariate polynomials by 'Dm = { p (', a) I a ED} .
According to the coefficient positivity assumption, all coefficients of sta-
ble polynomials have the same sign. Therefore, without loss of generality,
..
we may assume a~~12 i", > 0 and
(4.29)
It is clear that the last condition implies that 'Dm C Pnl,n2, ... ,n",.
(4.30)
{ p(O) (8) ± r [J.L ± (1 - J.L) SkJ s~a I J.L E [0, 1J}, k = 1,00', m. (4.31)
Here
l", = { n", - 1, if a = k, (4.32)
n"" if a # k,
where a E {I, 00', m}. There are exactly C;' = m subregions of this type.
And for everyone of them there are defined exactly 4m one-parametric
families.
Now fixing two of the indices, say a and f3 (a # (3), we define subregion
{p(O) (s) ± r [J.L ± (1 - J.L) SkJ s~as~ IJ.L E [0, 1J}, (4.33)
k=l,oo.,m; k#a, k#j3,
where l", and lf3 are defined in (4.32), and a,j3 E {l,oo.,m}, a # f3. The
number of such subregions is C~ = m( m - 1) /2.
Now fixing v < m of the indices, say a, f3, 00., 'Y (they are all distinct), we
can define subregion
where la, l{3, ... , l"( are defined as above. In this way we find C:;" subregions
of this type.
Continuing this process, we define the last subregion
where the boundary of the value set is defined by the following 4m one-
parametric families:
4.7 Notes
For the case of time-invariant systems, there is a strong connection between
the exponential stability of a system and the location on the complex plane
of the zeros of its characteristic function. For the case of time-delay sys-
tems of concentrated delays, the characteristic function is a finite sum of
potential and exponential factors.
The book [9] remains one of the best references for those who are in-
terested in the study of zeros of such functions. It contains an exhaustive
analysis of the location on the complex plane of zeros of quasipolynomials
from several interesting classes. A much more general and profound study
of zeros of quasipolynomial functions may be found in the Ph.D. disser-
tation of E. Schwengeler, see [241]. In our presentation of the exponential
and potential diagrams, we followed this dissertation. Example 4.1 is taken
from [152].
The Edge Theorem has been proved first for the case of retarded type
quasi polynomials in [70], and then has been extended to the case of neutral
type quasipolynornials in [71].
Stability analysis of an edge subfamily of quasipolynomials, including
conditions for convex directions, Theorem 4.7, may be found in [138]. Sta-
bility conditions for an interval quasipolynomial are also given in [138].
Basic results on applications of multivariate polynomials for stability
analysis of quasipolynomials are given in [131]-[133].
In presenting of stability conditions for multivariate polynomials, we fol-
lowed the publications [141]-[143].
Part II
where Ao and Al are given n x n real matrices. The usual initial condition
is in the form of
Xo = rjJ. (5.2)
We will defer the discussions on the uncertainties and systems with multiple
delays as well as distributed delays to later chapters.
We will use the Razumikhin Theorem (Theorem 1.4 in Chapter 1) and
the Lyapunov-Krasovskii Stability Theorem (Theorem 1.3 in Chapter 1) to
discuss the stability of the system. We will restrict ourselves to using the
bounded quadratic Lyapunov function or Lyapunov-Krasovskii functional,
and aim at arriving at stability criteria that can be written in the form
of Linear Matrix Inequalities (LMI) or a closely related form. Efficient
numerical methods are available to solve LMIs. The readers are encouraged
to read Appendix B to become familiar with the basic concepts of LMI
and a number of useful topics especially relevant to time-delay systems,
notably the variable elimination technique in LMI and quadratic integral
inequalities (Jensen Inequality). We will frequently refer to them in this
chapter.
The first part of this chapter will mainly use the Razumikhin Theorem.
Section 5.2 discusses delay-independent stability criteria. The first part dis-
cusses systems with a single delay expressed by (5.1), which can be derived
using the simplest form of Lyapunov function. The delay-independent sta-
bility criteria test if a system is asymptotically stable for arbitrary delay. It
is intuitively obvious that such a criterion would be very conservative if the
delay r is already known and small. It turns out that delay-independent
stability criteria effectively consider the delayed term, i.e., second term on
the right hand side of (5.1), to be always detrimental to the stability of
the system. Therefore, such criteria are useful when the system behavior is
148 5. Systems with Single Delay
dominated by the first term on the right hand side of (5.1) (i.e., Al is rather
small) and when we want to obtain a quick assurance that the delayed term
will not destabilize the system. Section 5.2 continues with a discussions on
the delay-independent stability of systems with distributed delays, which
is useful for deriving the delay-dependent stability criteria in Section 5.3.
Section 5.3 derives simple delay-dependent stability criteria using the
Razumikhin Theorem. In practice, as indicated in Example 1.2 in Chapter
1, there are systems that require feedback to improve the system perfor-
mance, and the feedback channel involves delays. In these cases, the delay-
independent stability criteria are clearly insufficient. In many cases, it is
more desirable that the system does not have a delay, and the presence
of delay is indeed detrimental to system stability and performance. It is
therefore natural to transform a delayed system into a system without de-
lay plus distributed delays, and to treat the term with distributed delay as
disturbance. This process is known as model transformation. We can use
the delay-independent stability criterion for distributed delay discussed in
Section 5.2 to obtain a stability criterion, which turns out to depend on
the size of delay. This is because the magnitude of the disturbance term is
closely related to the size of time delay. An important observation for such
a criterion is that the delay is always considered detrimental to stability.
Indeed, a necessary condition for such a criterion to be satisfied is that it
is satisfied for any smaller delay.
The delay-dependent stability criteria so arrived are still conservative.
In addition to the potential conservatism caused by applying the delay-
independent stability criterion to the system with distributed delays result-
ing from the model transformation, we will show that the model transfor-
mation process introduces spurious poles not present in the original system.
The dynamics represented by such spurious poles are known as additional
dynamics. As delay increases from zero, it is possible for the additional
dynamics to become unstable before the original system does. Therefore,
the model transformation process itself may introduce conservatism.
We will also show that it is possible to derive a simple delay-dependent
stability criterion without explicit model transformation. It turns out that
the resulting stability criteria include as special cases the delay-independent
and delay-dependent stability criteria derived in the earlier part of Section
5.3.
Sections 5.4 and 5.5 uses the Lyapunov-Krasovskii method to derive
delay-independent and simple delay-dependent stability criteria, parallel
to the results obtained by using Razumikhin Theorem. It is interesting to
notice that some resulting stability criteria take a very similar form as the
ones obtained using Razumikhin Theorem. The Razumikhin results can be
obtained from the Lyapunov-Krasovskii results by imposing additional con-
straints. While this may give the impression that the Razumikhin results
are more conservative, it will be shown in Chapter 6 that the Razumikhin
5.1 Introduction 149
(5.3)
and its derivative along the system trajectory V(x( t)) satisfies
(5.4)
whenever
V(x(t +~)) ~ pV(x(t)), - r ~ ~ ~ ° (5.5)
V(x) :S v(x).
satisfies
(5.7)
°
when K > is sufficiently large. _
We will refer to (5.6) as the Lyapunov-Krasovskii functional condition,
and (5.7) as the Lyapunov-Krasovskii derivative condition.
5.2 Delay-independent-Razumikhin Theorem 151
( P Ao +
Aip
Air P + oP PAl)
-oP <.
° (5.9)
(5.10)
is satisfied. Now consider the derivative of V along the trajectory of the
system (5.1),
. d
V(x(t)) = dt V(x(t)) = 2xT (t)P[A ox(t) + AIX(t - r)]. (5.11)
Whenever Xt satisfies
where
rPOr = (xT(t) xT(t - r) )T.
The inequality (5.9) implies that for some sufficiently small 8 > 0, p = 1+8,
PAo+AirP+opP
( Aip PAl)
-oP <,
°
152 5. Systems with Single Delay
subject to
Then, since the (1,1) entry of the matrix in (5.9) has to be negative definite,
any a satisfying (5.9) must satisfy
0< a < a.
For the least conservative result using this formulation, it is necessary to
solve (5.9) for all a in the interval (0, a). The sweeping of a may be compu-
tationally costly. Alternatively, we may choose to fix a = a/2 and accept
5.2 Delay-independent-Razumikhin Theorem 153
where (3 ~ 0 and r > O. For each given (3, we choose a = aj2 and check the
feasibility of (5.9). A bisection process is used to find (3max, the greatest (3
for (5.9) to be feasible. It was found that (3max = 0.9. 0
Based on the above dual form of stability criterion, state feedback sta-
•
bilization design is more convenient. For example, given the system
such that the closed loop system is stable independent of delay. Since the
closed loop system is
( AoQ+BVo+QA6+VlBT+aQ AIQ+BVI) 0
(AIQ + Bvif -aQ < ,
which is an LMI if a is fixed. Such a process is not directly possible based
on Proposition 5.3.
Most of the stability criteria discussed in this chapter have a dual form,
which can be obtained by a similar procedure. We will not discuss these
dual forms individually to avoid repetition.
( pa(O)P
AT(O)p
- R(O) PA(O) ) ( x(t)
-a(O)P x(t + 0)
) dO
.
With the above expression, and the fact that p > 1 can be arbitrarily close
to 1, (5.21) and (5.22) imply
Proposition 5.6 The conditions in Proposition 5.5 are satisfied for the
l:
system (5.19) if and only if they are also satisfied for the following system
(5.24)
The stability conditions in Proposition 5.5 for this new system are therefore
Rnew(B) _r R(_r B)
ruew rnew
This completes the proof. •
With arbitrary time delay, as the coefficient matrix Al of the delayed term
grows from zero, the system performance deteriorates and eventually the
system is in danger of losing stability.
However, in practice, one often encounters nominal open-loop systems of
the form (5.25) which are unstable. Feedback with small delay (due either
to measurement or control delay) is needed to stabilize the system. It is
therefore useful to take another point of view. Rewrite the system (5.1) in
the following form
x(t) = [Ao
+ i:+ Al)X(t)
[-AIAox(t + 0) - AIAIX(t - r + O))dO, (5.26)
where
</J(O) , -r 5: 0 5: 0,
1/;( 0) = {
solution of (5.1) with initial condition (5.2), 0< 0 5: r.
(5.28)
Therefore, the system described by (5.1) and (5.2) is embedded in the sys-
tem described by (5.26) and (5.21) without the initial condition constraint
158 5. Systems with Single Delay
where
(5.30)
pAo + A5' P + i 0
2r R(O)dO < 0, (5.34)
( ~(O)P
AT(O)P
- R{O) pA{o)
-a{O)P
)
< 0, - 2r ~ 0 < 0. (5.35)
-r ~ 0 < 0,
a{O)
-2r ~ 0 < -r,
R(O)
-r ~ 0 < 0,
-2r ~ 0 <-r
completes the proof. _
The stability criterion in Proposition 5.7 depends on the time delay r
and is therefore delay-dependent. We can eliminate the arbitrary matrices
Ro and Rl among the three matrix inequalities (5.32) to (5.33) to arrive
at the following equivalent form.
°
Corollary 5.8 The system described by (5.1) is asymptotically stable if
there exist real symmetric matrix P > and real scalars ao > 0, al >
such that
°
(5.36)
where
P > 0, (5.39)
!r [P{Ao + Ad + {Ao + Alf P] + aP < 0. (5.40)
160 5. Systems with Single Delay
ao > 0,
al > 0,
ao +al < 0,
for the existence of fixed ao and al such that the LMI (5.36) is feasible.
We can, of course, choose ao and al a priori within this region and accept
the conservatism. A reasonable choice seems to be
ao = (5.41)
2(IIAIAoll + II AlIl2)'
ollAlll2
(5.42)
Theorem 5.9 The transformed system described by {5.29} and {5.30} with
initial condition {5.31} is equivalent to the following system
I:
where
At I: z(t + B)dB.
det ( S1 - A o 0- e- Al
rs
.6.t (8) =
.6.a (8).6. 0 (8),
where
.6. 0 (8) = det(81 - A o - e- rs Ad (5.44)
is the characteristic quasipolynomial of the original system (5.1), and
1 -rs )
.6.a ( 8) = det ( 1 - - : AI.
Therefore, the poles of the transformed system consist of the poles of the
original system satisfying .6. 0 (8) = 0, and additional poles, which are the
solutions of
1 -rs )
.6. a (8) = det ( 1 - -: Al = O. (5.45)
It is clear that although the stability of the transformed system implies the
stability of the original system, the reverse is not necessarily true. As the
delay r increases, it is possible that one of the additional poles may cross
the imaginary axis before any of the poles of the original system do. It is
therefore possible that the original system is stable but the transformed
system is unstable. This is the stability implication of the additional dy-
namics.
To further characterize the additional dynamics, we can conclude from
(5.45) that
where Ai, is the ith eigenvalue of matrix AI. Let 8 = 8ik, k = 1,2,3, ... be
all the solutions of the equation
5.3 Simple delay-dependent-Razumikhin Theorem 163
Then Sik, i = 1,2, ... , n; k = 1,2,3, ... are all the additional poles of the
system described by (5.29) and (5.30). We will refer to Sik as the kth
additional pole corresponding to Ai (or corresponding to the ith eigenvalue
of AI)'
It is of interest to consider the trend of the additional poles Sik as the
delay r ~ 0+. We will state the following result without proof.
Proposition 5.10 For any given AI, all the additional poles satisfy
lim
r-O+
Re(sik) = -00.
As a result of this proposition, all the additional poles have negative real
parts for sufficiently small r. As r increases, some of the additional poles Sik
may cross the imaginary axis. It turns out that the exact value when that
happens can be analytically calculated. We will state the follow theorem,
again without proof.
k7r+ L(Ai)
r = rik = Im(Ai) > 0, k = 0, ±1, ±2, ....
Corresponding to a positive real eigenvalue Ai of AI, there is an additional
pole on the imaginary axis if and only if
1
r = Ai'
No additional poles corresponding to a negative real eigenvalue Ai of Al
will reach the imaginary axis for any finite delay.
Given a time delay r, we can use the above theorem to find the region
r r, such that all the additional poles (zeros of ~a (s))have strictly negative
real parts if and only if all the eigenvalues of Al lie in r r:
'm
~; f··~····~····~··~····~····~·,~·
····~····~·~~···~····.~
.. ~~~
. , ! 'r'
Ao:z
., E------~----~-c~-"----_+----_+1rh_ "
-~:~~~~~~~~
-f. ~· · · · · ·~·'· · · · · ~· · · · ~l=
•
Since all the eigenvalues of Al are real and negative, all the additional poles
of the transformed system are on the left half plane for any r. Therefore
the stability of this system and the transformed system are equivalent.
The conservatism shown in Example 5.2 is due to the application of the
Razumikhin Theorem, not due to the model transformation. 0
Example 5 .4 Consider
.
x(t) =
( -6
0.2 -~.8 ) x(t) + (~8 !8) x(t - r).
We can write
Vex) 2x T (t)P[A ox(t) + AIX(t - r)]
< 2x T (t)P[A ox(t) + AIX(t - r)]
+ lOr ( xT(t) ~xT(t+~) )
({r ~)(
x(t)
t~x(t + 0
) d~. (5.49)
166 5. Systems with Single Delay
(5.50)
where
x(t+e) )
Xtr(e)= ( x(t+e-r) . (5.51)
we have
Using (5.50) and (5.52) in (5.49), carrying out integration for terms in
(5.49) whenever possible, we arrive at
V(x(t))
<_ T (0) ( ~ + (0: + O:Or + 0:1 r )pP PAl - Y) (0)
x tr (PAl _ y)T -o:P Xtr
( ~+(0:+0:0r+0:1r)P PA1-y)
{PAl - y)T -o:P < 0, (5.53)
are satisfied.
We can eliminate the matrix variable Z in (5.48) and (5.54), then elimi-
nate X in the resulting matrix inequality and (5.53) to obtain the following
equivalent form.
-YA o
-aoP
(5.55)
°
°
where
To discuss the region of the scaling factors a, ao, and aI, first notice
that the (1,1) entry of (5.56) defined in (5.57) implies that ao and a1 must
lie in the triangular region ao > 0, a1 > 0, ao + a1 < ct, where ct is defined
in (5.38) to (5.40).
To find a bound for a, assume the matrix w = (Ao AI) has full row
rank, let w= (Air AT) T be a right inverse of w. We can left-multiply
168 5. Systems with Single Delay
M -PA1w- -W 1 TR )
'T < 0,
A
(
TTl
-w Al P - ;:IW -R - o:rw Pw
A
where
R = diag (o:oP O:lP),
and
1
M = -(PAo + Ao P)
A
r
T
+ (0:-r + 0:0 + 0:1 ) r
1
P - 2 (o:oAo PA o + 0:1A1 PAt}.
AT A AT A
(5.58)
Therefore, given 0:0,0:1, let &(0:0,0:1) be the solution of the following GEVP
all(o:l - A1)Aoil
(5.61)
3(11(0:1 - A1)Aoil + 11(0:1 - A 1)A1ID'
all(o:l - A 1 )A1 11
(5.62)
(5.63)
We can first set 0: = 0 in calculating 0:0 and 0:1 using (5.61) and (5.62),
then calculate 0: using 0:0 and 0:1 just calculated in (5.63), and update 0:0
and 0:1 using the newly calculated 0: (we can, of course, also update 0: using
the newly calculated 0:0 and 0:1).
Example 5.6 Consider the same system discussed in Example 5.2. Using
the method discussed here with 0:,0:0, and 0:1 calculated by (5.61) to (5.63),
with 0:0, 0:1, and 0: all updated, the calculation indicates that the the
system is asymptotically stable for r < rmax = 0.9476. For r = r max , the
scaling factors are 0: = 0.4813, 0:0 = 0.7160, and 0:1 = 0.6208. This is an
improvement over the method using explicit model transformation. 0
5.4 Delay-independent-Lyapunov-Krasovskii Theorem 169
where the matrices P and S are symmetric and positive definite. We will
write the above as
It can be easily calculated that the derivative of V along the system tra-
jectory is
- lOr ¢T (0)8(0)¢(0)dO.
The derivative condition will be satisfied if the two matrices in the above
are negative definite.
P > ° (5.73)
and real symmetric matrix functions R(O) and 8(0), such that
and
8(0) - R(O)
( AT(O)p PA(O))
-8(0) <,
°o E [-r, OJ. (5.75)
172 5. Systems with Single Delay
S( B) = a( B)P.
P > 0, (5.76)
-PA 1 A o -PA~)
-So 0 < 0, (5.77)
o -SI
where
1 T
M = -[P(A o + AI) + (Ao + AI) PJ + So + SI. (5.78)
r
5.5 Delay-dependent-simple Lyapunov-Krasovskii functional 173
R(O) { Ro,
Rb
-r ~ 0 < 0,
-2r ~ 0 < -r,
(5.79)
If one is to consider the stability of only the system (5.1), the above cri-
terion is clearly better than Proposition 5.8 since it is less conservative
and computationally more convenient because of the linearity of parame-
ters. The main value of Corollary 5.8 is again because it is also valid for
time-varying delay, as will be discussed in Chapter 6.
Example 5.8 Consider the same system discussed in Example 5.2. Using
Proposition 5.16, we can find that r max , the maximum delay for the system
to satisfy the condition, is r max = 1.000. This is less conservative than the
result in Example 5.2 obtained using the Razumikhin Theorem. However,
it is still very far from the true delay limit for stability r~~ytic = 6.1725.
o
(5.81)
174 5. Systems with Single Delay
where
cf>(~ + 0), - r ~ 0 ~ 0, - r ~ ~ ~ 0,
Aocf>(O) + A1cf>(-r).
Realizing
we can calculate
where
M = PAo+A~P+rA~ZAo +8,
cf>~r = (cf>T(O) cf>T(-r»).
Letting
T
T ·T
cf>oiJ = ( cf> (0) cf> (0) ) ,
we have
o ~ 1 0
-r
T(X Y)
cf>oiJ yT Z cf>oiJ dO
= rcf>T(O)Xcf>(O) + 2cf>T(0)y(cf>(0) - cf>(-r»
+ [Or ;l (O)Z¢(O)dO. (5.83)
(5.84)
we obtain
. T ( N PAl+rA~ZAl-y)A.
V(cf»~cf>or symmetric -8 + rA[ZAl 'l'Or'
where
N = P Ao + A~ P + r A~ ZAo + 8 + r X + y + yT,
5.6 Complete quadratic Lyapunov-Krasovskii functional 175
( N -Y)
(5.81) and
PAl +rA6'ZAI 0 (5.85)
symmetric -S+rAiZAl <.
To illustrate the connection between this result and the previous results
based on Lyapunov-Krasovskii functional methods, we eliminate the matrix
variable Z in (5.85) and (5.81) using Proposition B.6 in Appendix B to
obtain the following.
Example 5.9 Consider the system discussed in Example 5.8. Using the
above stability criterion, it can be concluded that the maximum delay for
stability is r max = 4.3588. This is much less conservative compared to the
corresponding result obtained using explicit model transformation. 0
5.6.1 Introduction
From the numerical examples presented so far in this chapter, it is obvious
that there is still substantial conservatism for the simple delay-dependent
176 5. Systems with Single Delay
If the control u(t) is set to zero, the system is clearly unstable because of
"negative damping." Classic control theory would suggest using a deriva-
tive feedback u(t) = -kx(t). Indeed, a control gain k > 0.1 will make the
system stable by achieving positive damping. It is well known that deriva-
tive feedback is usually not easy to achieve. It seems reasonable to use the
following finite difference
x(t) - x(t - r)
u(t) = x(t - r) - x(t) = -r--'--'--'----'- (5.88)
r
(5.89)
Recall that (5.89) has a positive definite solution for any given positive
definite matrix W if and only if the system (5.90) is asymptotically stable.
Indeed, if both P and Ware positive definite, we can choose Lyapunov
function
V(x) = xTpx,
5.6 Complete quadratic Lyapunov-Krasovskii functional 177
and verify that its derivative along the trajectory of the system (5.90) is
V(x) = xT(PA+ATp)x
= -xTWx.
Furthermore, the solution in this case can be explicitly expressed as
P= 1 00
cpT (t)WCP(t)dt,
where
cp(t) = etA
is the fundamental solution of the system (5.90), which satisfies the equa-
tion
ci>(t) = Acp(t),
and the initial condition
cp(O) = I.
Indeed, it will be shown that this idea can be extended to time-delay sys-
tems in that it is still possible to prescribe the derivative expression, and
"solve for" the Lyapunov-Krasovskii functional.
From (5.96), it can be easily concluded that q,{t) also satisfies the equation
d
dt q,(t) = q,(t)Ao + q,(t - r)A1' t ~ 0 (5.97)
and initial condition (5.94). Since the system (5.91) is asymptotically stable,
the stability exponent ao is negative. Therefore, we can find an a, ao <
a < 0, such that
11q,(t)11 ~ KeO'.t, for all t ~ 0 (5.98)
for some K > 1. In other words, the fundamental solution approaches zero
exponentially.
With the fundamental solution q,{t), for a symmetric matrix Wand real
scalar r, we can define the matrix Uw{r) as
Uw{r) = 1 00
q,T{t)Wq,{t+r)dt. (5.99)
This integration is well defined because q,(t) vanishes for t < 0 and ap-
proaches zero exponentially as t -+ +00. We will omit the subscript W
when no confusion may arise. For r ~ 0, with a transformation of integra-
tion variable ~ = t + r, and considering the symmetry of W, we write
UT(r) 10 00
q,T(t + r)Wq,(t)dt
1 q,T(~)Wq,{~
00
- r)~
= 10 00
q,T(~)Wq,(~ - r)d~.
In the last step, we have used the fact that q,(~ -r) vanishes when 0 < ~ ~
r. Thus we have arrived at the useful fact
U(-r) = UT(r). (5.100)
:t (q,T(t)Wq,(t))
= q,T(t)Wq,(t)Ao + q,T(t)Wq,(t - r)A 1
+A6q,T(t)Wq,(t) + Afq,T(t - r)Wq,(t).
5.6 Complete quadratic Lyapunov-Krasovskii functional 179
Integrating both sides of the above from 0 to 00, we arrive at another useful
fact
W + U(O)Ao + A5U(0) + UT(r)Al + AfU(r) = O. (5.101)
Differentiating (5.99) with respect to r using (5.97) and (5.100), we can
also write
d~ U(r) = U(r)Ao + UT(r - r)Al' r E [0, r]. (5.102)
With the above discussions, we can easily show that U(r), r E [0, r] can be
written as the solution of a two point boundary value problem as stated in
the following theorem.
Theorem 5.18 Matrix U(r) is the solution of the second order ordinary
differential equation
where
(5.107)
Let x(t,cjJ) be the solution of (5.91) with initial condition Xo = cjJ, then
it can be easily confirmed that
(5.108)
vw(cjJ)
cjJT (O)Uw (O)cjJ(O)
(5.111)
by
v(Xt)
Or more explicitly,
v( ¢) = ¢T (O)U(O)¢(O)
+2¢T(0) lOr U(-r - O)AI¢(O)dO
where
U(O) = UWl+rW2+W3(O).
Indeed we can prove the following.
< -8I1x(t)1I2
182 5. Systems with Single Delay
for 0 < 8 < Amin(H). Integrate from 0 to 00, noticing lim x(t) = 0, we
t->oo
obtain
Therefore
where
p = pT E IRnxn , (5.118)
and for all -r ::::: ~ ::::: 0 and -r ::::: 7] ::::: 0,
+2 lOr ¢T(OS(~)¢(~)d~.
- lOr ¢T(~)S(~)¢(~)d~
- lOr d~ lOr ¢T(~) [:~R(~,,,,) + :",R(~,,,,)] ¢(",)d",
+2¢T(0)[PAI - Q(-r)]¢(-r)
Recall that the system is asymptotically stable if and only if there exists a
Lyapunov-Krasovskii functional (5.117) such that the Lyapunov-Krasovskii
functional condition
(5.123)
and its derivative (5.122) satisfies the Lyapunov-Krasovskii derivative con-
dition
(5.124)
for some E: > O.
5. 7.2 Discretization
The basic idea of the DLF method is to divide the domain of definition of
matrix functions Q, R, and S into smaller regions, and choose these ma-
trix functions to be continuous piecewise linear, thus reducing the choice of
184 5. Systems with Single Delay
h = rlN.
Then
p=0,1,2, ... ,N.
This also divides the square S = [-r, OJ x [-r, OJ into N x N small squares
Spq = [Op, Op-IJ x [Oq, Oq-IJ. Each square is further divided into two triangles
{(Op+ah,Oq+,8h)
o S; ,8 S; 1,
I OS;aS;,8 }
,
~q {(Op+ah,Oq+,8h) I }.
The continuous matrix functions Q(~) and S(~) are chosen to be linear
within each segment Lp , and the continuous matrix function R( ~, ".,) is cho-
sen to be linear within each triangular region ~~ or T),q. Let
Qp Q(Op),
Sp S(Op),
Rpq = R(Op,Oq).
Then, since these functions are piecewise linear, they can be expressed in
terms of their values at the dividing points using a linear interpolation
formula, i.e., for 0 S; a S; 1, P = 1,2, ... , N,
N N N N
V( ¢) = ¢T (O)P¢(O) + 2¢T (0) L VQP +L L VRpq + L Vsp, (5.128)
p=1 p=lq=1 p=1
where
and
¢(P)(a) = ¢(Op + ah). (5.132)
(5.133)
Furthermore, it also satisfies
if
Sp > 0, p= 0,1, ... ,N. (5.135)
( ~T
Q
q
R+S
~) >0 (5.136)
186 5. Systems with Single Delay
(5.137)
(5.138)
(5.139)
and
W(N-l)(a)
W(N)(a)
= Qp-l"p(p)(l) -1 1
(Qp-l - Qp)"p(p) (a)da
VRpq 11 [1 1
¢(p)T(a)R(pq) (a, (3)¢(q) ({3)hd{3] hda
-1 1 (m(p;~a,{3)1/J(q)({3)d{3]hda
-1 1
1/J(p)T(a)[(Rp_1,q_1 - R p,q-1)1/J(q) (1)
p=l
N
~ 1 1
¢(p)T(a)S(p)(a)¢(p)(a)hda
~1
N 1
= ¢(p)T(a)[(l_ a)Sp + aSp-1]¢(P) (a)hda
+~
N-1
1 1
[a¢(p+1)T(a)Sp¢(p+1) (a) + (1- a)¢(p)T(a)Sp¢(p) (a)]hda.
188 5. Systems with Single Delay
N
~ l' [([
r
LVsP > ¢(1)({3)hd{3) T 80 ( [ ¢(I)({3)hd{3)
p=1
1
S(~) = 'h(Sp-1 - Sp),
1
Q(~) = 'h{Qp-1 - Qp),
8R{~,,,,) 8R(~, 7]) 1
8~ + 87] 'h{Rp-1,q-1 - Rpq).
-LL 11 </>(p)T(a)Rdpq</>(q)((3)dad(3
N N 1 1
p=1q=1 0 0
11 +
+2</>5r [D S (1 - 2a)Da]¢(a)hda, (5.146)
where
</>Or = (</>(0)
</>(-r)
)
, (5.147)
_
</>(a) =.
C'(a) )</>2 (a)
, (5.148)
</>N (a)
~ll
Rd12
( Rd21 Rd22 ~'N
Rd2N )
Rd = . , (5.155)
DS Div ), (5.157)
DpS (5.158)
Dop =
(5.159)
DIp = (5.160)
Da (5.161)
vap (5.162)
Vop (5.163)
Dfp = (5.164)
(5.166)
(5.167)
5.7 Discretized Lyapunov functional 191
V(¢) = -1 1
( ¢6r[D S + (1- 2a)Da] Jt (a) )
for arbitrary U. If (5.166) is satisfied, then use the Jensen Inequality (Propo-
sition B.8 in Appendix B) to the first term, we have
( U
-1
-1) ( ¢(a)
Sd
lDs + (1- 2a)DajT ¢Or ) da
( U
-1
-1) Jo{1 ( ¢(a)
Sd
lDs + (1 - 2a)Da]T ¢Or ) da
(
~_lDaUDaT
_D s :;'
_Ds
R+Sd
)
> °
are satisfied, the Lyapunov-Krasovskii derivative condition (5.124) is sat-
isfied. The arbitrary matrix function U can be eliminated from the above
inequality and (5.166) using Proposition B.6 in Appendix B to yield (5.167) .
•
5.7.5 Stability criterion and examples
We can summarize from the above discussion the following.
Proposition 5.22 The system with single time delay described by (5.116)
is asymptotically stable if there exist n x n matrices P = pT; Qp, Sp = SJ,
p = 0,1, ... ,N; Rpq = Rfp, p = 0,1, ... , N, q = 0,1, ... , N such that (5.136)
and (5.167) are satisfied, with notations defined in (5.137) to (5.139) and
(5.149) to (5. 164}.
192 5. Systems with Single Delay
SN > 0,
Sp-1 - Sp > 0, p = 1,2, ... ,N.
This implies
So> Sl > ... > SN > 0,
which implies (5.135).
Since both the Lyapunov-Krasovskii functional condition and Lyapunov-
Krasovskii derivative condition are satisfied, the system is asymptotically
stable. _
Clearly the criterion is in the form of LMI since all the system parameters
appear linearly. The following examples are provided to show the effective-
ness of the DLF method. The first example illustrates the convergence to
the analytical limit as N increases.
We want to calculate the maximum time delay rmax the system can tolerate
and still retain stability. Proposition 5.22 can be used to check the stability
of the system for any given r. A bisection process is used to estimate r max
for different N. The analytical stability limit is also listed for comparison.
The results are very close to the analytical limit. Even the result for N = 1
is sufficient for most practical applications, and is much better than the
ones obtained in earlier sections. The convergent trend to the analytical
limit is obvious. 0
Analytical
6.17258
This is obtained from the system (5.87) by applying feedback control (5.88).
This system is clearly unstable for r = O. To apply Proposition 5.22, the
system is written in the standard form
d ( x(t) ) ( 0 1 ) ( x(t) )
dt x(t) -2 0.1 x(t)
+ ( 01 0)
0
( x(tx(t -- r)r) ) .
The system is stable for r E (rmin' rmax). Proposition 5.22 with different N
is used to estimate rmin and rmax. The computation consists of sweeping
through a range of r with relatively large step size, and a bisection process
near the lower limit rmin and the upper limit rmax. The result is listed in
the following table along with the analytical limit. 0
N 1 2 3 Analytical
rmin 0.1006 0.1003 0.1003 0.1002
rmax 1.4272 1.6921 1.7161 1.7178
5.8 Notes
5.B.l Results based on the Razumikhin Theorem
An example of the results prior to the modern LMI approach is by Thowsen
[260]. There is a large amount of work done on this topic. The delay-
independent formulation in Proposition 5.3 and the dual form is equivalent
to Proposition 18 in Niculescu et al. [214]. The delay-dependent formulation
in Proposition 5.7 is equivalent to Proposition 19 in Niculescu et al. [214].
Some studies try to avoid the scaling factor and write the stability criteria
in LMI form. For example, it can be shown, using variable elimination
of LMI, that Li and de Souza [173] effectively choose the scaling factors
0:1 = 0:2 = 1 when specialized to the case of Proposition 5.7. It seems
clear that fixing the scaling factors are essential if one is to restrict the
formulation to LMI form. However, it is also essential that these scaling
factors should be proportional to the magnitude of the system matrices Ao
and AI. For example, if Ao and Al are multiplied by 1000 and the delay
r is divided by 1000 (such as choosing a different unit, say second versus
millisecond), this accelerates the evolution of the system without affecting
the stability of the system. A stability criterion with the scaling factor not
proportional to Ao and Al would conclude that the stability of these two
system may be different. The scaling scheme proposed here is from Gu and
Han [94].
The Razumikhin Theorem-based formulation with implicit model trans-
formation described in (5.53) and (5.54) was first proposed in Park [221] for
194 5. Systems with Single Delay
the special case of a = ao = a1. The forms presented here with connections
to other criteria is from Gu [92].
( yTX y)
Z =
(P~ArV
+ WT) rV -1 (
P+ W
1
TV
A)
1 .
The formulation without such a constraint was given in Park [221] with
extensions to considering uncertainty and cost function by Lee, Moon, and
Kwon [165]. The elimination of variable with connections to other delay-
independent and delay-dependent stability criteria was given by Gu [92].
The idea of a complete quadratic Lyapunov-Krasovskii functional was
proposed by Repin [236]. The special case of x being scalar was solved by
Datko [48]. The case of W 2 = W3 = 0 was solved by Infante and Castelan
[122]. Kharitonov and Zhabko [145] proposed the case W2 =I 0 and W3 =I 0
in studying the robust stability.
5.8 Notes 195
The uncertainty set characterizes the uncertainties and serves as basic in-
formation needed to carry out robust stability analysis. Notice also that the
coefficients may depend on the time t as well as the current and previous
state variable x(t+~), -r ::; ~ ::; O. For the sake of convenience, we will not
explicitly show these dependences or only show the dependence on time t
when no confusion may arise.
A good choice of uncertainty set is a compromise between minimizing
conservatism (and therefore, it is desirable to make the uncertainty set
"small") and the mathematical tractibility (and therefore, it is desirable
to make the uncertainty set structurally simple). In this section, we will
discuss some common uncertainty structures.
In addition to uncertainties in the coefficient matrices, the time delay r
may also involve uncertainties. In this chapter we will mainly discuss the
case where r = r(t) may be time-varying and is known to be within the
interval r(t) E (0, rmax]. However, we will reserve the term "uncertainty"
to coefficient uncertainty.
W = Wn + aWa + f3w,e.
Then, the uncertainty set possesses four vertices
W(I)
= Wn + aminWa + f3 mi nw ,e,
w(2)
= Wn + amaxWa + f3 minW,e,
w(3)
= Wn + aminWa + f3 max w,e,
w(4) Wn + amaxWa + f3 maxW,e.
He
r\ co {W( i ),2, -- 1, 2 , ... , nv }.
nv
f3i ~ 0, 'Lf3i = 1.
i=l
Proposition 6.1 The set 0 formed by np parameters, ai, i = 1,2, ... , np,
appearing multilinearly, with each parameter varying between a lower limit
and an upper limit
o -- {W
m,n
+ a m+l Wm,v I a(wm +l,min:$ a m+l :$ a m+l,maJ"
w ) E 0m
m,n, m,v
for some
2m
f3i ~ 0, :L:f3i = 1.
i=l
Therefore,
i=l
which can be obtained by letting each parameter assume its lower and
upper limit. •
As will be seen later, subpolytopic uncertainties are easy to treat in
time-domain stability analysis, since most stability criteria are in the form
of matrix inequalities with system matrices appearing linearly. Such ma-
trix inequalities are satisfied for all the possible system matrices in the
uncertainty set if and only if they are satisfied at the vertices.
There are, however, still practical difficulties. For example, for the case of
multilinear uncertain parameters, there are nv = 2np vertices of uncertainty
set if there are np uncertain parameters. The number nv can be very large
for a rather modest np. For example, np = 10 uncertain parameters will
result in nv = 1024. This phenomenon, known as combinatorial explosion,
limits the practical applications of the subpolytopic uncertainty model to
only a few uncertain parameters. There are also other uncertainties that
cannot be reasonably modeled by a subpolytopic uncertainty set with a
reasonable number of vertices. In such cases, norm-bounded uncertainties
are often used.
W = Wn + t::.w,
202 6. Robust Stability Analysis
i.e.,
Ao Aon + .6.Ao,
Al = Aln + .6.AI.
The uncertain part is written as
(6.1)
We will say that the uncertainty description (6.3) is well posed if (6.4)
is satisfied. The LF norm-bounded uncertainty can be motivated by the
following system
Ilu(t)112 ~ Ily(t)112.
6.2 Uncertainty characterization 203
From this setting, it may appear at first sight that the bounding (6.2) is
rather restrictive. For example, it is often desirable to model the uncertain
matrix F in (6.3) to satisfy
VVT P,
UTU Q-l,
and let
F* UFV, (6.9)
E* EU- l , (6.10)
Go = V-lGO, (6.11)
Gi V-lGl, (6.12)
D* V-lDU- l , (6.13)
Lemma 6.2 For any real matrices E, G and real symmetric positive defi-
nite matrix P, with compatible dimensions,
204 6. Robust Stability Analysis
Proof. Since
(EPE T + GT P-1G) - (EG + GT ET)
= (EP - GT)p-l(EP _ GT)T
> 0,
the conclusion is easy to see. _
IIXII ~ 1 (6.20)
if and only if there exists a real scalar). > 0 such that
P + )'GT G
( -ET E + )'GT D )
+ )'DTG -).(I - DT D) < 0, (6.21)
or equivalently,
Proof. In view of (6.18), the inequality (6.19) is well posed and is equiv-
alent to
-ep~ +2eE(I _XD)-lXG~ < 0 (6.23)
for all ~ =f O. The inequality (6.23) can be written as
or equivalently
(6.31)
We can properly scale the uncertain matrices Fi, i = 1,2, ... , m similar to
(6.9) to (6.13) such that they satisfy
Indeed it is often possible to extract all the uncertain elements from the
system and consider them as feedback. Then the system with all the un-
certain elements removed can be regarded as a multi-input-multi-output
206 6. Robust Stability Analysis
system. This process allows us to write the uncertain system in the stan-
dard block-diagonal form (6.32) to (6.35) and is known as "pulling out
uncertainties."
Equations (6.32) to (6.34) can be written as (6.5) to (6.7), with
(6.36)
which justifies the name "block-diagonal uncertainty." Notice also that the
block-diagonal uncertainty includes as a special case the independent un-
certainty case described by (6.14) to (6.17), since (6.14) and (6.15) can be
written as
where Iii and Iir are the identity matrices with the dimensions the same
as the number of rows and columns of F i , respectively. Then an important
observation is that
KIFK;l = F
for arbitrary ki' i = 1,2, ... , m. Therefore, we can rewrite (6.7) as
U = KIFK;ly,
or
(6.39)
where
If the stability criterion is satisfied for some ki' i = 1,2, ... , m, then the
system is stable. Many stability criteria can be written as LMI with either
k? or 11k? appearing linearly. The possibility of searching through all the
possible scaling factors significantly reduces the conservatism. It has been
argued that this scaling process allows one to reach "very close to the
optimum" in terms of reducing conservatism.
whenever
V(x(t + 0)) < pV(x(t)) for all - rmax :::; 0 :::; O. (6.47)
where
¢Or = (xT(t) xT(t - r) )T.
We can therefore conclude the following.
a> 0,
for all
±(t) = Ao(t)x(t) + 10
-r(t)
A(t, B)x(t + r(t,B))dB, (6.50)
with
(Ao(t),A(t,·)) E Or for all t 2: 0, (6.51)
°< r(t, B) :::; rmax for - r :::; £I < 0, (6.52)
where Or is compact. For the sake of convenience, we will often suppress
the notation of the explicit dependence of time t when no confusion may
arise. It is possible to write the system in a simpler form without sacrificing
the generality. For example, by reparameterizing £I, one may easily change
the integration limit in (6.50) to [-1,0]. However, since the main purpose
of studying such a system with distributed delay here is to derive the delay-
dependent stability criteria for systems with a single delay, we choose to
keep the formulation as close to the original form as possible. The study of
distributed delay in its own right will be conducted in Chapter 7.
We will study this system again using Lyapunov function V(x) as in
(6.45). The system is asymptotically stable if (6.46) is satisfied whenever
(6.47) is satisfied. For any a(t, B) 2: 0, -r(t) :::; £I :::; 0, when (6.47) is
satisfied,
+10 a(t,B)[pxT(t)Px(t)
-r(t)
_x T (t + r(t, B))Px(t + r(t, £I))] dB
0 T ( M(t,B) PA(t,B))
= 1-r(t) ¢O(J AT(t, B)P -a(t, B)P ¢O(JdB, (6.53)
where
M(t, B) = ~[PAo(t)
r
+ A6(t)P] + a(t, B)pP + R(B, Ao(t),A(t,')) (6.54)
and
_ ( x(t) )
¢0(J- x(t+r(t,B)) .
210 6. Robust Stability Analysis
The matrix function R((), Ao(t), A(t,·)) (which we will write R(()) for the
sake of convenience) is constrained to satisfy
/
0 R( ())d() = O. (6.55)
-r(t)
Proof. If (6.56) is satisfied, the (2,2) entry ofthe matrix in (6.56) implies
P>O, (6.57)
N(O) < O.
Using (6.55), the above implies
/
0 N(O)dO
-r(t)
1 0
-r(t)
R( O)dO "? 0 (6.60)
with x(t + 0) replaced by the right hand side of the system equation (6.42)
with appropriate time shift, we can obtain the transformed system
-r(t)
[Ao(t + O)x(t + 0)
+AI(t + O)x(t + 0 - r(t + O))]dO,
which can be written as
x(t) = Ao(t)x(t) + 1 0
-2r(t)
A(t, O)x(t + T(t, O))dO, (6.61)
where
Ao = Ao +Al }
A(O) = -AlAoo, A( -r + 0) = -AlAlO, -r ~ 0 < 0 ,
for (Ao,Ad E n and (Aoo,AlO) En, 0 E [-r,O]
(6.62)
and
0, -r(t) ~ 0 < 0,
r(t,O) = { 0 _ r(t + B), -2r(t) ~ 0 < r(t).
This model transformation, of course, also introduces additional dynamics.
These additional dynamics are more complicated and will not be discussed
here. Applying Proposition 6.6, we can obtain the following.
( Nk
-(PAlAkO)T
-PAlAkO )
-ak P < 0, k = 0,1 (6.63)
No = 1 [P(Ao + A 1) + (Ao
-2 + Ad T P] + aoP + R,
rmax
1 T
N1 -2-[P(Ao + A 1) + (Ao + Ad P] + alP - R.
rmax
Proof. Apply Proposition 6.6 to the transformed system (6.61) and un-
certainty set (6.62) with the selection
0.0 -r(t) < B ~ 0,
a(t, B) ={ 0.1 -2r(t) < B ~ -r(t),
and for -r < 0 ~ 0
R(B) R(Ao, A 1, Aoo, A 10 ),
R(O - r) -R(Ao, A 1, Aoo, A 10 ),
where Ako = Ak(t + 0), k = 0,1. •
I:
Krasovskii functional
II - ( PAo +A5 P
1- Aip
+S P_ASl) <0 (6.70)
W
(i) _
-
(A 0
(i) A(i»
' 1
. - 1, 2 , ... , n v ,
, 2-
•
Example 6.1 Consider the following uncertain system,
where
Ip(t)1 ::; 0.1.
This is a system with subpolytopic uncertainty. The uncertainty set n has
.
two vertIces (A(i)
0 ' A(i»).
1 ,z = 1 , 2 ,were
h
( -2.1
-1.1
-0.1
-0.1 )
-1.0 '
AP) = 0: (
-1 ~0.9 ),
A(2)
o
= (
-1.9 0.1 ) A~2) = 0: ( -0.9
~1.1 ) .
0.1 -0.8 ' -1.0
The calculation using Proposition 6.10 indicates that the system is stable
independent of delays for 0 ::; 0: < O:max = 0.6730. 0
P>O
and 8 such that
PAon + AInP + 8 + GIGo PAIn + GIG I PE + GID )
( AinP+GiGo -8 + GiG I GiD < O.
ETp+DTG o DTG I -I +DTD
(6.76)
Proof. Use Proposition 6.9; the system is stable if there exist a P' > 0
and a 8' such that
where
p = -
( P' A On
AT
In
P'
+ ATOn P' + 8' P'A In
-8' ), (6.78)
E ( ~'E ) , (6.79)
G = ( Go GI ). (6.80)
Use Lemma 6.3, (6.77) is satisfied for allllFll :::; 1 if and only ifthere exists
a >. > 0 such that
P = !p'
>. '
8 !8'
>. .
Corollary 6.12 The system described by (6.81) with the coefficient uncer-
tainty described by (6.72) to (6.75) is asymptotically stable if there exists
a symmetric matrix P and scalar a > 0 such that (6.76) is satisfied with
8=aP.
Proof. Since Proposition 6.11 is derived from Proposition 6.9, the result
is immediately clear with the observation of the connection between Propo-
sition 6.5 and Proposition 6.9, and the fact that P and 8 in Proposition
6.11 is proportional to P and 8 in Proposition 6.9 . •
There are other cases in this chapter that the Razumikhin Theorem-
based results can be easily read out from the corresponding Lyapunov-
Krasovskii functional-based results by a simple substitution. The readers
should not have difficulty identifying these cases.
6.4 Delay-independent-Lyapunov-Krasovskii functional 217
(6.82)
To be specific, let
m m
Fi E jRPiXqi, P = LPi' q = L qi· (6.83)
i=l i=l
where
AOn ( -2 0 )
o -0.9 '
( -0.5 0 )
-0.5 -0.5 '
218 6. Robust Stability Analysis
E = /31, D = 0.51, Go = GI = 1.
Using Proposition 6.11, we can conclude that the system is asymptotically
stable independent of delays for /3 < /3 max = 0.3378.
If D = 0 and all the other parameters remain the same, then the uncer-
tainty is norm-bounded. We can calculate /3 max = 0.6073. 0
with
(Ao(t),A(t,.)) En for all t 2: 0, (6.85)
where n is compact, and
V(¢) = ¢T(O)[PAo+A~P+lorS(O)dO]¢(O)
+2¢T(0) l : PA(O)¢(O)dO
- lOr ¢T(O)S(O)¢(O)dO,
(6.86)
we have
6.5 Delay-dependent-simple Lyapunov-Krasovskii functional 219
where
M(O) = !(PAo + A5' P) + S(O) + R(O, Ao, A(·». (6.88)
r
We can therefore conclude the following.
M(O)
( AT(O)p PA(O) )
-S(O) < 0, - r ~ 0 ~ 0
for all
(Ao,A(·» En,
where M(O) is defined in (6. 88}.
Notice in the above that M(O) also depends on the uncertain system
matrices Ao and A( 0) due to R. In practice, we may choose R to be inde-
pendent of the uncertain matrices and accept the additional conservatism.
As in the case of the Razumikhin Theorem-based formulation, an impor-
tant application of this delay-independent stability criterion is to derive a
delay-dependent stability criteria for systems with discrete delays. In such
case, it is essential to allow R to be dependent on the uncertainty such that
it can be analytically eliminated to result in a simple formulation.
where
k=~+~ }
A(O) = -A1AolJ,A(-r + 0) = -A1A11J, -r ~ 0 < 0 .
for (Ao,A1) E 0 and (AolJ,A11J) EO, 0 E [-r,OJ
(6.92)
As discussed earlier, the stability of the transformed system implies the
stability of the original system. However, the reverse is not necessarily true
due to the presence of additional dynamics.
Applying Proposition 6.14 to the transformed system, we can conclude
the following.
( Nk
-(PA1AklJ)T
-PAlAklJ )
-Sk < 0, k = 0,1 (6.93)
-r < 0 ~ 0,
-2r < 0 ~ -r,
and for -r < 0 ~ 0
6.5 Delay-dependent-simple Lyapunov-Krasovskii functional 221
P>O,
8k {=: ok P ,ok>Ofork=O,l,
r {=: Tmax.
W
(i)
= (A(i) A(i». 1, 2, ... , nv
0' 1 , 1 =
where
( ~Ao(t) ~AI(t) ) = EF(t) ( Go G1 ), (6.97)
and
IIF(t)11 ~ l. (6.98)
We can state the following.
Mn -PA1nAon -PA1nA ln
-A6'n A fn P -So + IlG6'Go IlG6' Gl
-AlnAtn P IlGiGo -St + IlGiGt
ETp 0 0
~(Go + G 1) -G1AOn -GlAln
-(PAInE)T 0 0
PE ~(Go + GI)T -PAInE
0 -(GtAOn)T 0
0 -(GtAIn)T 0
-I 0 0 < 0, (6.99)
0 -I -GtE
0 -(GIE)T -Ill
where
Mn = ![P(Aon
r
+ Aln ) + (Aon + Alnf P] + So + Sl'
where
and
IlFoll ~ 1. (6.103)
Using (6.95) to (6.97) in (6.94), we obtain
where
P -
( M.
-A~oATnP
-AlOAln P
-PAlnAoo
-So
0
-PAlnAlo
0
-Sl
),
E (fE),
G = ( ~(Go + Gd -GlAOO -GlAlO ).
for some A > O. Dividing A and using the Schur complement, the above is
equivalent to
_!P !E GT )
( t~T ~J 0 < O.
G 0 -J
Redefine
1
~P as P,
1
~So as So,
1
~Sl as Sl'
Then
-P + EFo(t)G < 0,
224 6. Robust Stability Analysis
where
-PA1nA1n
0
PE
0
~(Go +
-(G1AOn)T
Gt)']
-81 0 -(G1 Alnf ,
0 -f 0
-G 1Aln 0 -f
Using Lemma 6.4, the above is equivalent to the existence of a real scalar
> 0 such that
f.1.
ET ( i!.
+ f.1.GT G E ) 0
-f.1. f <,
which is (6.99). Notice, f.1. > 0 is implied by (6.99) . •
Example 6.3 Consider the system (6.89) with norm-bounded uncertainty
described by (6.95) to (6.98) with
where
W(O) = W3 + (r + 0)W2'
U(O) = UWl (0) + lOr UW2(0)dr + UW3(0)
= 1 00
<pT(t)[WI + rW2 + W3]<P(t + O)dt,
and <p(t) is the fundamental solution of the system (6.105). When WI > 0,
W2 ~ 0, W3 ~ 0, we may easily conclude that V(¢» ~ -c:l1¢>(0)11 2 and
V(¢» ~ O. A slight modification of V(¢» allows us to satisfy the condition
V(¢» ~ £11¢>(0)11 2 .
In this section we will explore applying the Lyapunov-Krasovskii func-
tional (6.107) to the uncertain system
i:(t) = [Aon + ~Ao(t)]x(t) + [Aln + ~AI(t)]X(t - r). (6.108)
We will show that although we may set W2 = W3 = 0 in systems without
uncertainty, these two matrices play important roles in robust analysis. The
purpose is to illustrate the idea rather than obtain the least conservative
condition. We will also assume the uncertainty satisfies
(6.109)
where Po > 0 and PI > 0 are real scalars, and R.n, R I , and S are positive
definite matrices. As indicated earlier, this can also be reformulated in the
standard norm-bound uncertainty. The uncertainties ~Ao and ~AI are
assumed to be independent.
We still use <p(t) to represent the fundamental solution of the nominal
system without uncertainty (6.105). Taking the derivative of the Lyapunov-
Krasovskii functional along the trajectory of the uncertain system (6.108)
yields
l:
. T
V(¢»I(6.108) = V(¢>)I(6.105) + 2[~Ao¢>(0) + ~A1¢>( -r)] .
Theorem 6.18 If the nominal system (6.105) is stable, then the uncertain
system (6.108) to (6.109) is asymptotically stable if we can find a scalar
JL > 0, matrices WI > 0, W 2 2': 0, and W3 2': 0 such that
WI > (1 + r) Po S + JLU T (O)(RQl + Rll )U(O),
JL
W2 > JLAinUT(-r-O)(RQl +Rl1)U(-r-O)Aln,
W3 > (1 + r) PI S.
JL
6.7 Discretized Lyapunov functional 227
where the system matrices are bounded by a known compact set n in the
following manner,
Use the process almost identical to the case with systems without uncer-
tainty described in Section 5.7 of Chapter 5: using Lyapunov functional
( ~T
Q
q -) > 0
R+S
(6.114)
and
(6.115)
228 6. Robust Stability Analysis
~o ~N
~l
( Rll
R ~10 ~IN
)
, (6.117)
~1I
Rdl2
Rd
( Rd21 Rd22 ~>N
Rd2N
.
)
, (6.125)
(6.127)
(6.128)
(6.129)
(6.130)
(6.131)
(6.132)
(6.133)
(6.134)
6.7 Discretized Lyapunov functional 229
(i) A(i)) .
(A 0 ' I ,Z = 1, 2 , ... ,nv ·
where
and
IIFII ~ 1. (6.138)
In spite of the apparent complexity, it can be seen that (6.114) does not
involve system parameters, and (6.115) has a structure very similar to
the delay-independent stability case, and much simpler than the delay-
dependent case using model transformation. Taking advantage of this struc-
ture, we can write (6.115) as
where
(6.140)
(6.141)
(6.142)
I!.T
(Q q -) > 0
R+S
(6.144)
and
Lln - LlG -D~
( _D~T Rd+Sd
_DaT 0 (6.145)
n
EJ, E;
are satisfied. In the above,
A (G'{;Go G'{;Gt)
~G = GiGo erG t '
6.8 Notes 231
and other notation are defined in (6.116) to (6.134) and (6.140) to (6.142),
with the subscript "n" indicating that Ao and Al should be replaced by the
nominal values Aon and A In , respectively, in the corresponding expressions.
Proof. It has already been shown above that (6.143) and (6.144) are
sufficient for stability. Divide both by .x, and redefine ±P, ±Qp, ±Sp, and
± Rpq as P, Qp, Sp, and Rpq , respectively. •
-1
Aon = ( -2
0
0 )
-0.9 ' A ln = ( -1 o)
-1
and
E = 0.21, Go = GI = I.
Using Proposition 6.20, we estimate the maximum delay r max such that
the system remains asymptotically stable for r < rmax. The computational
results using different discretization N are listed in the following table.
Similar to the case without uncertainty, even for the coarsest discretization
of N = 1, the result is much less conservative than the corresponding
result obtained using the simple Lyapunov-Krasovskii functional method
in Example 6.3 obtained using Proposition 6.17. 0
6.8 Notes
6.8.1 Uncertainty characterization
The uncertainty characterization has been extensively studied in the liter-
ature on robust stability and control. Earlier literature include the study
by Yakubovich [287)-[289], which is equivalent to the norm-bounded un-
certainty discussed here. Indeed, an equivalent of Lemma 6.3 was obtained
in these papers. The special case Lemma 6.4 was obtained by Petersen and
Hollot [224). The generality of block-diagonal uncertainty was shown in the
"pulling out the uncertainty" process by Doyle, Wall, and Stein [58). The
poly topic uncertainty and subpolytopic uncertainty were discussed in Ho-
risberger and Belanger [115]' Boyd and Yang [22), and Gu [82). Most of the
results on manipulating these uncertainties in the context of LMI can be
found in the book by Boyd et al. [21 J.
232 6. Robust Stability Analysis
L Ai(t)X(t - ri),
K
x(t) = (7.1)
i=O
where
o= ro < r1 < ... < r K = r, (7.2)
and
w(t) En, for all t 2: 0, (7.3)
and w represents all the coefficient matrices:
(7.4)
and n is a compact set. Using the Lyapunov-Krasovskii functional
V(¢) = ¢T(O)p¢(O) + ~
K
iTi
0
¢T(O)Si¢(O)dO,
where
pT >0, (7.5)
sf> 0, i = 1,2, ... ,K, (7.6)
we can easily calculate
K
+2¢T (O)P L Ai¢( -ri)
i=1
- L ¢T (-ri)Si¢( -ri).
K
i=1
Let
and
K
PAo +A5'P+ LSi PAl PA 2 PA K
i=l
Aip -S1 0 0
IlK = AIp 0 -S2 0 (7.7)
AT
KP 0 0 -SK
7.3 Simple delay-dependent-multiple delays 235
Then . -T-
V(¢» = ¢>KilK¢>K'
We can therefore conclude the following.
( -Ri
At PAi)
-8i < 0, i = 1,2, ... ,K. (7.11)
x(t) =
KKK
2:= Ai(t)X(t) - 2:= 2:=
i=O i=1 j=O
10
-ri
A(t)Aj(t + O)x(t + 0 - rj)dB.
236 7. Systems with Multiple and Distributed Delays
where P and Sij, i = 1,2, ... , K; j = 0,1, ... , K are symmetric positive
definite matrices. We can easily obtain
-247 (0)PLL
K
i=l j=o
K
1 0
-ri
Ai {t)Aj {t + €)¢{€ - rj)~
K K 0
- LLI ¢T(€-rj)Sij¢{€-rj)~,
i=1 j=O -ri
where
K 1 K
- L R 1j - ;- LUi + M, (7.14)
j=1 1 i=2
and
l(K K KK)
M = - rl Pt;Ai + t;ATP+ t;{;riSij . (7.15)
From the above, we can obtain the following simple delay-dependent sta-
bility criterion.
7.3 Simple delay-dependent-multiple delays 237
Note that Sij > 0 has already been implied by (7.17). We can restrict
Rij and Ui to be constant matrices and accept the additional conservatism.
We can also eliminate these matrix functions to arrive at the following
equivalent stability condition.
Corollary 7.3 The system with multiple delays described by (7.1) and
(7.4) is asymptotically stable if there exist symmetric matrices P, Sij,
i = 1,2, ... , K; j = 0, 1,2, ... , K, such that (7.16) is satisfied, and
M PA1We PA2We PAKWe
(PA1Wef *Sl 0 0
(PA2Wef 0 /2S2 0 >0
o o
is satisfied for all wEn and we E n, where
Si = diag (SiQ Si1 SiK ),
wand We are defined in (7.18) and (7.19), respectively.
where
1 .
-Ui, 2> 1,
ri
238 7. Systems with Multiple and Distributed Delays
This is equivalent to
This system includes, as special cases, the systems with multiple discrete de-
lays and distributed delays. The complete quadratic Lyapunov-Krasovskii
functional is very similar to the case with single delay. However, the deriva-
tion is far more sophisticated. In this section we will present the main idea.
Interested readers are referred to the literature discussed in the Notes sec-
tion for technical details.
Let <I>(t) denote the fundamental solution of the system (7.20), i.e., it
satisfies
<i>(t) = lOr d[F(B)]<I>(t + B),
with initial condition
<I>(t) = {I,0, t = 0,
t < 0.
For a real positive definite symmetric matrix W, let
7.4 Complete quadratic functional for general linear systems 239
where
Obviously
U&(r) = Uw(-r).
Based on Uw(r), for a asymptotically stable system, we can again prescribe
(7.21)
using the Lyapunov-Krasovskii functional
vw(¢) = 1 00
xT(t, ¢)Wx(t, ¢)dt,
Xo = ¢.
As shown by Huang [117], vw(¢) can be explicitly expressed in terms of ¢,
+ lOr 1° dv¢T(v)d7J[F(7])f
x(t) =
K
~AiX(t - ri) + lr°A(B)x(t + B)dB, (7.23)
lr°{3((}) IdF((}) I ~ ~~
K
~{3((}i)IIF((}i) - F((}i-l)II·
-r=8o<82< . ·<8K=0 -
Since F has bounded variation, the above is well defined. Especially for
l:
{3((}) == 1, the above becomes the total variation, which we denote as J.t
p. = IdF((})I·
Since
lOr [II¢(O)112 -
l:
V(¢) = -(1 + 8p.)II¢(0)W + 8 II¢(r)II 2)1dF(r)1
< -II¢(0)II2,
7.5 Discretized Lyapunov functional-multiple delays 241
inequality (7.26) is satisfied. To show (7.25), consider, for 0 < c < y'8/j.l,
v(¢) = V(¢) - cll¢(O)112.
Direct calculation using (7.21) yields
2. Many quantities involve two indices: one related to the delays, which
we usually use i or j; the other related to the discretization, which
we usually use p or q. We often need to manipulate these two indices
such as changing the order of summation.
3. The appropriate application of the quadratic integral inequality for
the S terms to render as small as possible the conservatism due to
overlapping of integration interval.
+ [Or ¢T (~)S(~)¢(~)d~
V(¢) = ¢T(O)P¢(O)+2L¢T(O)
K
i=l
1° Qi(~)¢(~)d~
iri°¢T(~)Si(~)¢(~)d~
-Ti
K
+~
where Qi, Si, and Rij are continuous matrix functions, and
P pT E jRnxn, (7.31)
Qi(~) E jRnxn, (7.32)
SiT(~) Si(~) E jRnxn, (7.33)
RijT(~,'fJ) Rji('fJ,~) E jRnxn, (7.34)
fori=1,2, ... ,K; j=1,2, ... ,K.
Since the additional functions Qi(~), Si(~), Rij (~, 'fJ), and RiK (~, 'fJ), i, j =
1,2, ... , K - 1 are introduced only to account for discontinuities, we will
constrain these functions to the following special forms without loss of
generality:
Taking the derivative of V in (7.30) with respect to time along the tra-
jectory of the system (7.28), carrying out integration by parts similar to
the single delay case, we can again write V(cp) as a quadratic expression of
cp:
K K
V(cp) - LLcpT(-ri)~ijcp(-rj)
i=O j=O
K
+2LLcpT(-ri)
i=O j=l
K
1° IIij(~)cp(~)d~
-ri
where
and
K-1
rrOj(~) = A~Qj +L Rij + RKj(O), (7.42)
i=1
K-1
rr OK (~) = A~ QK (~) +LRiK (~) + RKK (O,~) _ QK (~),
i=1
rr ij (~) AT Qj - Rij, 1:::; i :::; K - 1, 1:::; j :::; K - 1,
rr iK (~) = ATQK(~) - RiK(~), 1:::; i:::; K -1,
rrKj(O AIQj - RKj(-r), 1:::; j:::; K -1,
l1KK (~) AkQK(~) - RKK(_r,~). (7.43)
Recall that the system is asymptotically stable if there exist an c: > 0 and a
Lyapunov-Krasovskii functional (7.30) such that the Lyapunov-Krasovskii
functional condition
V(¢» ~ c:11¢>(O)W (7.44)
and the Lyapunov-Krasovskii derivative condition
(7.45)
are satisfied. The conditions are necessary and sufficient if there is no un-
certainty.
7.5.2 Discretization
For systems with multiple delays, it is essential that the division of the delay
interval [-r, 0] is compatible with the delays such that -ri, i = 1,2, ... , K-1
are among the division points. In other words, let ()p, P = 0,1, ... , N be the
dividing points
7.5 Discretized Lyapunov functional-multiple delays 245
then
-ri = {}N" i = 1,2, ... , K.
Thus each interval [-ri,Oj is divided into Ni intervals. Let the length of
pth segment be hp
hp = {}p-l - (}p.
For the sake of convenience, define
No 0, (7.46)
ho = 0, (7.47)
hN+1 0. (7.48)
We have
°= No < Nl < ... < N K = N,
and
Ni
ri = Lhp, i = 1,2, ... ,K.
p=l
°
The matrix functions QK (~), SK (~), RiK (~), and RKK (~,,,,) are again
chosen to be piecewise linear as follows: for ~ a ~ 1, p = 1,2, ... , N,
and for °
~ a ~ 1, ° /3 :::;
~ 1, p = 1,2, ... , N, q = 1,2, ... , N,
RKK ({}p + ahp, (}q + /3hq)
!
= R KK (pq)(a,/3)
(1 - a)R{;/ + /3R{;5.,q-l + (a - /3)R{;_~,q,
a'?:. /3;
(7.52)
(1- /3)R!// + aR!/_~,q_l + (/3 - a)R:':-l'
a < /3.
Thus, the Lyapunov-Krasovskii functional V is completely determined by
P , Q i , QK
P' Si 'P'
SK Rij ' P
RiK' RKK . ,')-- 1
pq , Z "
2.. "K - l',p,q-
- "1 .. " N ' °
As will be seen later, we often need to manipulate terms involving two
types of indices: one related to the delay, which we usually use i or j, the
other related to the discretization, which we usually use p or q. For this
purpose, it is useful to introduce the notation
(7,53)
246 7. Systems with Multiple and Distributed Delays
i if N i - 1 < p :::; N i ,
o if p = O.
MN; i, (7.54)
N Mp > P > N Mp - 1' (7.55)
Then, for any indexed expressions U;, i = 1,2, ... , Kj p = 1,2, ... , N, it is
often useful to change the order of summation as follows
K N; N K
LLU;
i=l p=l
LLU;,
p=l i=Mp
(7.57)
LLU;
i=l p=l
L L U;=L L U;, p=l i=Mp
(7.58)
L LU;
i=l p=l
L L U;=L L U;.
p=l i=Mp+l p=l i=Mp+ 1
(7.59)
In the last step of (7.58) and (7.59), we have used the convention (7.56)
and the fact that Mp = K > K -1 for p > N K - 1 .
V(¢)
(
~T ~ j;) (:(0) ) da
QKT il KT RKK w(a)
+L
K-1
i=l
1 0
-ri
¢T (~)Si¢(~)d~
+L
N
p=1 0
1 1
¢(p)T(a)SK(p) (a)¢(p) (a)hpda, (7.60)
where
e"
R12 R1,K-1
).
R21 R22 R2,K-1
R
~K-1,1 R K - 1,2 RK-I,K-I
RfK N
ilK (m
K
R~K R~K R'K
R2K
.
N ),
~fJ-I'K RK-I,K
I
RZ-I,K
~K
( RKK
10
RKK
01
RKK
11
RKK)
ON
RKK
IN
RKK . , (7.62)
RKK RKK RKK
NO NI NN
and
( J~", ¢(,)d<
f~r2 ¢(~)d~
)
<1>= (7.63)
~~rK-l ¢(~)d~ ,
248 7. Systems with Multiple and Distributed Delays
w(O)(a) ""(I)(a)
W(I)(a) ""(2) (a) + ",,(I)(a)
W(2) (a) ""(3) (a) + ",,(2) (a)
W(a) , (7.64)
""(p) (a) = hp i 1
¢(p) (f3)df3, p= 1,2, ... ,N, (7.66)
Proof. Similar to the single delay case, divide the integration interval
[-r,O] to the segments [Bp ,Bp _ 1 ], we have
K-l ° N
V(¢) = ¢T(0)P¢(0)+2L¢T(0)Qij ¢(~)d~+2L¢T(0)VQK(P)
~1 -~ p=1
K-l N
+ L Vs; +L VSK(p)
i=1 p=1
K-l N O N N
+2 L L j ¢T(~)VR;K(P)d~
+ LLVRKK(Pq),
i=1 p=1 -ri p=1 q=1
where
[Or; ¢T(~)Si¢(~)d~,
11 ¢(p)T(a)SK(P) (a)¢(p) (a)hpda,
Integrating by parts for VQK(P)' VRiK(p), VRKK(pq) similar to the single delay
case, and collecting terms, we may reach (7.60) . •
7.5 Discretized Lyapunov functional-multiple delays 249
If the last two terms involving Si and SK(p) in (7.30) do not exist, the
positive definiteness of the matrix in the first term of (7.60) is clearly a
necessary and sufficient condition for the Lyapunov-Krasovksii functional
condition (7.44). The extent of conservatism depends on how the terms
involving Si and SK(p) are incorporated to the rest of the quadratic ex-
pression.
fJ fl
F = (
f6 R
~:-1 ff-1
!pi = {I, p:S; Ni - 1 (or equivalently i 2 M p +1 ),
0, otherwise,
K-1
Sf
P S: + LSi, p=O,I, ... ,N. (7.72)
i=Mp +1
Proof. Let
250 7. Systems with Multiple and Distributed Delays
LOri ¢T(~)Si¢(~)d~
= E1
p=1 Jo
[1 ¢(p)T(a) [aS i + (1- a)Si]¢(P) (a)hpda
where the first two terms are of similar form to the expression of SK(p), we
can combine the first two terms of the above with the last term in (7.73)
to arrive at
Vs = Vs' + V,§, (7.74)
where
In view of (7.68) and (7.69), we can use the quadratic integral inequality
(Corollary B.9 of Appendix B) to obtain
N 1
Vs' = ~1 [a¢(p+l)T(a)S~</>(P+1)(a)hp+l
+(1- a)¢(p)T(a)S~¢(p)(a)hp]da
~ t(
p=oJo
[a[h p+ 1¢(P+1) (a)]T .!-S~[hp+l¢(P+1)(a)]
hp
11 wT(a)8w(a)da. (7.75)
7.5 Discretized Lyapunov functional-multiple delays 251
Similarly,
L
K-l 1
Vs > f t/J(N;)T(a)_I-Sit/J(N;)da
i=l 10 hN ;
11 E1 E1
T
11
t=l 0 p=O ' FO
V(¢)
( ~T ~+S
QKT
-SF ~;
(ilK _ SF)T R,KK + 8' + FTSF
) (:(0) )
w(a)
da.
-1 [11 a
( ¢T(o:) ¢T ((3) )
( ~KK R:faK ) (
o
~(o:)
¢((3)
) d{3] do:
,
(7.77)
da
where
ll.O1
ll.= COO ll.1O
.
ll.11 aOK
ll.1K
.
)
, (7.78)
[ K-1
ll.oo = - P Ao + A5 P + ~ (Qi + QiT + Si)
~I D~N
D82
( DB Df2 DfN )
DB = : 11 . ,
Dh Db D1<N
K-l
D ip
S
L..J hP (A!Qj
'"
t
_ Rij) + hp t p + QK
2 [A!(QK p + RiK
p-l ) _ (RiK )]
p-l'
j=Mp
K-l
D Kp L hp(AkQj - RWT)
DON)
D'lN
. ,
D'h . 1YkN
RK
ds12 dslN
Rt"
( Rlf21 RK RK )
RK
ds
ds22 ~lf2N ,
RK
dsNl RK
dsN2 RK
dsNN
K-1 K-1
RK
dspq = L hp(R iK q +
q-1 - RiK) L hq(RiKT
p-l _ RiKT)
p ,
i=Mp i=Mq
Proof. Dividing all integrations on the delay intervals [-ri,O] into inte-
gration on segments [Op, Op-l], with a change of variable such as
where
AT A
V~ ¢ !:l¢, (7.83)
Vn
K
L¢T(-ri)
i=O
K Nj
j=lp=1 0
1
LL
IIij(Op 1 + ahp)¢(p) (a)hpda, (7.84)
VSK L
N
p=1 0
1
1
¢(p)T(a)S~¢(p)(a)da, (7.85)
VnK
K-l N,
L L
i=1 p=1
1 0
1
¢(p)T(a)
VRKK LL
N
p=lq=1
N
1 1
¢(p)T(a)
[1 (1 8+ 18(38)
0
1
o hp 8a hq R KK (Op + ahp, Oq + (3h q)
¢(q) ((3)h qd(3] hpda. (7.87)
705 Discretized Lyapunov functional-multiple delays 255
YRK ~ t, i};, /.' ~(P)T [t. /.' (,,) (E,":, - E,K)~(q) (fJ)dfJ 1 hpda
where
KO
Rifl
RK
Rd12
RK
RK
dIN
RK
o
(
o o o
RK o
RK o
RK o
R dpq
o K
=
VRK =
o r </>(a)da
21 [10
1
A ] T
Ri; [ 10r </>(a)da
1
A ]
0
(7.89)
to write
VR~K = 210
1
[l~Jj ¢T({3)R:fsK ¢(a)da] d{3
210
1
[l~a ¢T(a)R:fsK¢({3)d{3] da. (7.92)
VR~K = 1 1
10 [10 ¢T(a)R:fsK ¢({3)d{3] da.
(7.93)
7.5 Discretized Lyapunov functional-multiple delays 257
> 0 (7.94)
E 0,
and
(7.95)
V(¢)
-1 1
( ¢T[DS + (1- 2a)Da] ¢t (a) )
-1 [1
1 Q
(¢T(a) ¢T((3)) (irffaK~K) ( ~~~~ )d(3] da.
Given
( U
-/
-/
S{f - W
)
> 0, (7.96)
258 7. Systems with Multiple and Distributed Delays
V(¢)
< -1 1
( Jt[D + (1 -
S 2a)Da] Jt (a) ) da
( U
-J
-~
Sd - W
) r
10
1
( lDs
¢(a)
+ (1- 2a)Da]T;P ) da
= - (;pT J; ;pT(a)da )
(~;;JDaUDaT SF~W+R:f.K +Rd~ ) (~;;P(a)da)
_11 [1° ( ;pT (a) J? (13) ) (~:aK WaK ) ( %~~~ )d13J da.
Therefore, the Lyapunov-Krasovskii derivative condition is satisfied if (7.96),
(7.95), and
D. - IDaUD aT _DB )
( (7.97)
-D.1- SK-W+RKK+R'K
d d. ds
>0
are satisfied. We may eliminate the variable U from (7.95) and (7.97) using
Proposition B.6 of Appendix B to obtain (7.94) . •
For the case of uniform mesh, since R:aK = 0, one can choose W = 0,
and (7.95) can be omitted, in which case (7.94) reduces to (5.167) if applied
to the single delay case.
Proof. According to Propositions 7.6 and 7.8, the system satisfies both
the Lyapunov-Krasovskii functional condition and Lyapunov-Krasovskii
derivative condition, and therefore is asymptotically stable, if (7.68), (7.69),
(7.70), (7.94), and (7.95) are satisfied. It remains to be shown that (7.68)
and (7.69) are not needed.
From (7.95), one concludes that W > 0, which together with the (3,3)
entries ofthe matrix in (7.94) implies (7.69). Examining the diagonal entries
of ..6. of the matrix in (7.94), one may conclude that (7.68) is implied by
(7.94) . •
To illustrate the effectiveness of the method, the following examples are
presented. The first example gives a sense of how close to the analytical
limit can be reached with a rather small N.
. = (-2 0 )
x(t) 0 -0.9 x(t)+ (-1 0 )
-1 -1 [0.05x(t - 0.5r) + 0.95x(t - r)].
The next example illustrates the case of polytopic uncertainty and in-
commensurate delays.
The next example considers the case where the system is unstable when
the delay vanishes.
1 )(Xl(t))+(O
( 0-1 0.1 0)(x 1(t-r/2))
X2(t) -1 0 X2(t - r/2)
00)(x1(t-r))
+( 1 0 X2(t - r) .
The stability criterion in Proposition 7.9 is used to find an interval of
r, [rmin' rmax], for which the system is guaranteed to be stable. A simple
search with uniform increment of 0.1 is conducted. Then a bisection near
the upper and lower bound is conducted to find rmin and rmax to within
0.001 for a given Ndl and N d2 . Using Ndl = Nd2 = 1, it was estimated
that [rmin, rmax] = [0.204,1.350]. For Ndl = Nd2 = 2, the interval can
be enlarged to [0.203,1.372]. It can be verified that for r = 0.2025, the
system has an imaginary pole at 1.0077i, and for r = 1.3723, the system
has an imaginary pole at 1. 3786i. Therefore, the exact stability interval is
(0.2025,1.3723) and these estimates are very accurate. 0
In other words,
The system matrices are uncertain and are bounded by a known compact
set
(7.102)
For the sake of simplicity, we will often suppress the time dependence of
the system matrices and write Ao, Ai, and A(O) instead of Ao(t), Ai(t) and
A(t,O).
A general quadratic Lyapunov functional will again be used:
where
P = pT E jRnxn , (7.104)
and for all -r S ~ S 0 and -r S 1] S 0,
Taking the derivative of V with respect to time along the trajectory of the
system (7.98), integration by parts when necessary similar to the single
delay case, one can write Vasa quadratic functional of ¢ as follows:
V(¢)
where
for some E > O. The conditions are necessary and sufficient if there is no
uncertainty in the system.
7.6.2 Discretization
Similar to the case of multiple delays, it is again important to allow nonuni-
form mesh in general, such that A(~) is independent of ~ in each segment.
Specifically, let 0 = 00 > 01 > ... > ON = -r be chosen such that
-ri = ON;, i = 0, 1,2, ... , K,
thus the intervals [-ri' OJ is divided into Ni segments [Op, Op-d of length
hp = Op-l - Op,
7.6 Discretized Lyapunov functional-distributed delays 263
(7.114)
We have the following.
V(¢»
11 (¢>T(O) wT(a)) ( ~T q
R
) (¢>(O)
w(a)
) da
+ ];
N
1 1
¢>(p)T(a)S(p)(a)¢>(p)(a)hpda (7.115)
if
Sp > 0, p = 0,1, ... ,N. (7.117)
264 7. Systems with Multiple and Distributed Delays
( ~T q
Q
-)
R+8
>0 (7.118)
Q ( Qo Q1 ... QN ), (7.119)
R01
R (~
~10 Rll ~N
~1N
)
, (7.120)
and
1l1(N-1)(a)
1l1(N) (a)
'l/J(p) (a) hp i 1
¢(p) ((3)d(3, (7.127)
of (7.115) as follows:
N
~VsP g ~
N
11
¢(p)T(a)S(p)(a)¢(p)(a)hpda
~
N
1 1
¢(p)T(a)[(l- a)Sp + aSp-1]¢(P) (a)hpda
11 ah 1¢(1)T(a) ~1 SOh1¢(1)(a)da
+ 11 (1- a)hN¢(N)T(a) h~ SNhN¢(N) (a)da
+ L r
N-1 1 [
ahp+1¢(P+1)T(a)~Sphp+1¢(P+1)(a)
P=l k ~1
+(1- a)hp¢(p)T(a) ~p Sphp¢(p) (a)] da
> r1 ah1¢(1)T(a)lSoh1¢(1)(a)da
Jo ho
We have
We will write
(7.129)
for the sake of convenience. Then, after the discretization, dividing the
integration interval [-r,O] into N segments [()p,()p-1], P = 1,2, ... ,N, the
derivative V of the Lyapunov-Krasovskii functional V along the trajectory
of the system described by (7.98) to (7.102) may be written as
-1 1
¢T(a)Sd¢(a)da
-1 [1" (
1
¢T (a) ¢T ((3) )
(°
R~a °
Rda ) ( ~(a) ) dfJ]]da
¢((3) ,
(7.130)
where
(7.131)
or more explicitly,
ES =
CII
E~l
Ef2
E~2 ElN
E~N )
· , (7.140)
Ea Cfl
E~l
Ef2
E~2 ~N
E~N )
· , (7.141)
hphq T )
E;q = -2- Ap (Qq + Qq-l , (7.142)
hphq T )
(7.143)
E;q -2- Ap (Qq - Qq-l ,
Rdsl2
Rds
C~U Rds2l Rds22 ~"N
R ds2N
·
)
, (7.144)
R dsNl RdsN2 R dsNN
1
Rdspq 2[(hp + hq)(Rp-l,q-l - Rpq),
+(hp - hq)(Rp,q-l - Rp-l,q)], (7.145)
R dal2
(~U
Rda21 Rda22 ~'N
Rda2N )
Rda = . , (7.146)
R daNl R daN2RdaNN
1
Rdapq = 2(hp - hq)(Rp-l,q-l - Rp-l,q - Rp,q-l + Rpq), (7.147)
268 7. Systems with Multiple and Distributed Delays
and
( ¢T (0) ¢T ( -r) ) ,
(¢(l)T(O!) ¢(2)T(0!) ... ¢(N)T(O!)).
V(¢)
< - (¢~r Io1 ¢T (O!)dO! )
_~DaUDaT + ~ -Ds - ~DaUEaT ) ( ¢Or )
(
_DsT - ~EaUDaT Sd - Esa + Rds - W I01 ¢(O!)dO!
-1 [l
1 CX
(¢T(O!) 4t(~)) (~ra ~a) ( :~~~ )d~] dO!,
(7.148)
provided
( U
-J Sd
-J) > O. (7.149)
> 0, (7.150)
> 0 (7.151)
Proof. The idea is very similar to the single delay and multiple delay
case. From (7.130), it can be verified that
V(¢)
= -1 (¢hE1
¢T(a)) (~I S:) ( ~f~ ) da
-¢6r fl.¢or + ¢6r [DSUD ST + ~DaUDaT] ¢Or
+2¢6r [DSUE ST + ~DaUEaT] 11 ¢((3)df3
for arbitrary matrix function U(Ao, AI, A2, ... , AK), where
-1 [l1 a
(¢T(a) ¢T(f3)) (irra ~a) ( ~~~~ ) df3 ] da
for arbitrary matrix W, where
E sa = E S + EsT + ~EaUEaT.
Therefore, the Lyapunov-Krasovskii derivative condition is satisfied if there
exist a matrix Wand a matrix function U such that (7.149), (7.151), and
(7.152)
Proof. This follows from Propositions 7.10 and 7.11. Similar to the sin-
gle delay and multiple delay case, condition (7.117) is already implied by
(7.150) . •
To illustrate the effectiveness of the discretized Lyapunov functional
method for systems with distributed delay, the following example is con-
sidered. The example considers a system resulting from a model transfor-
mation from a system with single pointwise delay.
x(t) = Ao(t)x(t) + 1 0
-r/2 AIX(t + ~)d~ + l-
-r
r2 2
/ A x(t + ~)d~, (7.153)
where
Ao ( 0.5-1.5 0 ),
-1
(7.154)
Al
( 02 2.5-0.5 ), (7.155)
A2
( 0 -1 ).
-1 0
(7.156)
This is the result of model transformation from the system with single delay
x. ()
t = (-3
1 -2.5)
0.5 x ()
t + (1.5
-0.5 2.5) x (t - r /2) ,
-1.5 (7.157)
7.7 Notes
Similar to the single delay case, there are numerous results available for
both delay-independent and delay-dependent stability results. See [214] for
an overview. For delay-indepent stability using Lyapunov-Krasovskii func-
tional method, see Ge, Frank, and Lin[74].
For a similar result to Corollary 7.3 using the Razumikhin Theorem, see
Cao, Sun, and Cheng [28] (similar to the single delay case, scaling factors
can be introduced). The derivation used here is original.
For systems with multiple delays, in addition to delay-independent and
delay-dependent stability, it is also possible to consider the concept of
mixed delay-independent and delay-dependent stability; see Kolmanovskii,
Niculescu, and Richard [159].
The extension to the general linear time-delay system of the complete
quadratic Lyapunov-Krasovskii functional (7.22) was reached by Huang
[117]. The modification (7.27) is along the same line as [145].
The discretized Lyapunov functional formulation for systems with mul-
tiple delays was first proposed in Gu [86]. The simplification to the current
form discussed here is based on Gu [90]. The expression of V(¢) in the
proof of Proposition 7.8 can be generalized to
V(¢)
-1 1
( 4/[DS + (1- 2a)Da] ;;/(a))
( U -J ) ([Ds+(1-2a)D a jT¢)d
-J S,f-aZ-(I-a)W ¢(a) a
-[1 1
¢(a)da]T (R:fsK + Rlf) [1 1
¢(a)da]
-1 [11 0
(¢T(a) ¢T({3)) (~:aK ~:aK) ( %~;j )df3] da.
The one used in the proof of Proposition 7.8 is the special case of the above
by setting Z = W. For nonuniform mesh with large difference of mesh size,
a stability criterion based on the above expression can be less conservative
but much more complicated.
The discretized Lyapunov functional method for systems with distributed
delay was proposed in Gu et al. [95]. The formulation presented here is
based on Gu [91].
Part III
Input-Output Approach
8
Input-Output Stability
8.1 Introduction
In this chapter we will discuss the input-output stability of time-delay
systems. The main purpose of presenting this formulation is to consider
the stability problem of systems under dynamical uncertainty. The input-
output stability framework discussed in this chapter is similar in flavor to
Chapter 3. A system is often written in the form of a nominal system and
an uncertain feedback, and stability is concluded based on a small gain
argument. However, the formulation here is more general. For example,
the nominal system may be a time-delay system, and the feedback may
be a more general uncertainty which is not necessarily due to time-delay
elements.
Based on the input-output framework, it is often possible to write the
stability criteria using the Lyapunov-Krasovskii functional formulation or,
in the case of linear time-invariant nominal systems, frequency sweeping
form. We will emphasize the Lyapunov-Krasovskii formulation in this chap-
ter since the frequency-sweeping form is very similar to Chapter 3.
Section 8.2 will discuss the basic concept of input-output stability and
the small gain condition for the preservation of input-output stability in a
feedback configuration. As a simple illustration, Section 8.3 discusses the
method of comparison systems. In this approach, the time-delay system
is embedded in a family of systems with a forward system without time
delay and a dynamical feedback uncertainty. The same LMI is obtained by
either a frequency domain approach using the Bounded Real Lemma, or
directly by a time domain approach. Then we begin a more comprehensive
discussion of the methodology of the input-output framework, with known
constant time delays considered as part of the nominal system rather than
uncertainty. In Section 8.4 we set up the scaled small gain problem. In
Section 8.5 we explain how the robust stability problem with dynamical
feedback uncertainty can be formulated as a scaled small gain problem and
discuss the uncertainty characterization. It will also be shown that some
robust stability conditions discussed in Chapter 6 are also sufficient for
dynamical uncertainty.
As an application of the above framework, in Section 8.6 we discuss the
approximation of time-delay elements in stability analysis. Specifically, we
discuss the approximation of time-varying delays by time-invariant delays,
276 8. Input-Output Stability
111112 = 100
111(t)1I 2 dt.
For any T > 0 and a function 1 : i+ -t IRn , the truncation operator (or
projection operator) P r is defined as
(FG)I = F(GJ).
A system H is linear if H(al + f3g) = aHI + f3Hg for all a,f3 E IR and
I, g E L 2e+. It is causal if it satisfies P rH = P rHP r for all T > O. In other
8.2 Input-output stability 277
satisfying
I(t, 0, 0) = 0,
h(t,O,O) = 0,
describe a system H. The mapping from u to y for zero initial condition
X(t) = 0, u(t) = 0, t S 0,
1'0(H) = IiH(s)lloo.
y = Gu, (8.2)
u = Fy+v. (8.3)
(8.4)
Condition (8.4) is known as the small gain condition. If the small gain
condition is satisfied, then we can easily find an upper bound of the gain
of the feedback system 1'0 (feedback(G, F)): from (8.2) and (8.3), we can
write
This indicates
1'0 (G)
1'0 (feedback(G, F)) $; (1 -l'o(Gho(F))"
In practice, the input v in (8.2) and (8.3) is often due to nonzero initial
conditions, in which case, we sometimes omit v and write
y = Gu,
u Fy.
We still say the above system is input-ouput stable to mean that the system
described by (8.2) and (8.3) is input-output stable.
1:
z(t) (Ao + C)z(t) + (AI - C)z(t - r) - C[z(t) - z(t - r)]
with feedback
Let uT = (uf uf) and yT = (yf yf). Then we can write the
forward system described by (8.7) to (8.9) as
y=Gu, (8.12)
u=ay. (8.13)
This casts the problem in an input-output form. From the small gain the-
orem, and the equivalence of stability for feedback (G x, F x) and that of
feedback(G, F), we can conclude that the system is input-output stable if
we can find a real nonsingular matrix X such that
< 1, (8.14)
< 1. (8.15)
6 = (~:), (8.18)
b (1 1 ~). (8.19)
Ul = Ll 1 Yb (8.20)
U2 = Ll2Y2, (8.21)
where Lll and Ll2 are described by (8.10) and (8.11), respectively. In fact,
both Lll and Ll2 are "linear scalar systems" with Hoo norms bounded by
one. In other words, the transfer function matrices of Llk are
(8.22)
is /'0 (G) =IIG(s)lloo < 1. We can obtain a stronger result by scaling. Let X
consist of nonsingular block-diagonal matrices of appropriate dimensions,
(8.25)
(8.26)
We can use the following Bounded Real Lemma to write (8.26) in the form
of a frequency-independent LM!.
pT =P >0 (8.28)
and
are satisfied.
8.3 Method of comparison systems 283
is satisfied if and only if there exists aPE jRnxn such that (8.28) and
are satisfied. From the above discussions, we can conclude the following.
P>O,
P(A 1 - C) + Air Z2 A I
-Zl + Aj Z2A I
rpc)
o < 0,
o -Z2
(8.31)
where
Proof. First, notice that (8.31) implies Zl > 0 and Z2 > 0 in view of
the (3,3) and (2,2) entries of the matrix. From the above discussion, the
system is stable if (8.30) is satisfied for some nonsingular X E X described
in (8.24). Left-multiply the second row of (8.30) by XT and right-multiply
the second column by X; then define
Clearly,
Z = {XTX I X EX}.
(8.32)
since
where
8.3 Method of comparison systems 285
lot uf(T)ZlUl(T)dT
= I:r yf(T)ZlYl(T)dT
Io yf(T)ZlYl(T)dT
t r
-
lot UnT)Z2U2(T)dT
= :2 lot [lOr Y2 (T + O)dO] T Z2 [lOr Y2 (T + O)dO] dT
< ~ lot [lOr ynT+O)Z2Y2(T+O)dO] dT
= ~ [Or [lot ynT + O)Z2Y2(T + O)dT] dO
~ ;1: [t' yi(T)z"l/2(T)dT] dO
where y = (yi yf) T represents the linear expression (8.8) and (8.9) of
x and u. If
f(t,O,O) 0, (8.46)
h(t, 0, 0) = 0, (8.47)
and x(t) E lR.n, u(t) E lR.m, y(t) E lR. m. Notice that we have constrained
u(t) and y(t) to have the same dimension. This constraint is imposed for
notational convenience and can be relaxed. Let X c lR. mxm be a given
set of real nonsingular matrices. The scaled small gain problem is to check
whether G is internally stable and there exists an X E X such that
The special case of X = {I}, i.e., check whether G is internally stable and
satisfies
l'o(G) < 1, (8.49)
is known as the small gain problem or 1£00 problem. When G is linear
and time invariant, a frequency domain technique can be used. Indeed, a
frequency-sweeping method is often possible to check (8.48). In the follow-
ing, we will discuss the time domain approach.
Let
z = {XT X I X EX}. (8.50)
Then, the scaled small gain problem is equivalent to finding whether there
exists a Z E Z such that
(8.51)
is satisfied for any t, u and y = Gu. To achieve (8.51), we have the following.
288 8. Input-Output Stability
°
proven that (8.48) is satisfied, or equivalently, (8.51) is satisfied. Integrating
(8.54) from to t, we obtain
where
for sufficiently small E > 0, and conditions (8.52) and (8.54) become
V(t, Xt) > Ellx(t)112, (8.55)
W(t,Xt,u(t)) < -Ellx(t)112 -Ellu(t)W. (8.56)
In this case it is often possible to reduce the conditions into an LMI prob-
lem.
As a simple illustration, consider a system G described by the following
equations
and
W(t,¢,u)
= -¢T(O)[-PAo - A5 P - Q(O) - QT(O) - S(O)]¢(O)
_¢T (-r)S( -r)¢( -r)
_[Or ¢T (~)S(~)¢(~)d~
- [Or d~ [Or ¢T(~) [:~R.(~,1]) + :1]R(~,1])] ¢(1])d1]
+2¢T(0)[PAl - Q(-r)]¢(-r)
u and the last term) in (8.62). Therefore, we can easily modify (5.146) in
Chapter 5 to obtain
DB = (~; ),
Da = (~; ),
and
D ua = h T( QI-QO Q2-Ql
-E QN -QN-l ).
2
Then the same manipulation as in Proposition 5.21 of Chapter 5 allows us
to conclude that (8.56) is satisfied if
(8.65)
Proof. The only thing that needs to be shown is that Bp > O. But it is
already implied by (8.65) in view of the (3,3) entry of the matrix in the
left-hand side. •
292 8. Input-Output Stability
The above proposition concerns the case with single delays. However, the
same ideas apply to the case where G is a system with multiple pointwise
delays or distributed delays with piecewise constant coefficients. We can
make similar modifications to the results in Chapter 7 to derive sufficient
conditions for 'Yo(Gx) < 1 in an LMI form.
Proof. This is a direct consequence of the small gain theorem and the
fact that the input-output stability of feedback(G, 4) is equivalent to the
input-output stability of feedback( G x, 4 x). •
Thus the robust stability problem under dynamical uncertainty is trans-
formed to the scaled small gain problem discussed in the last section.
(U[ ur ui ),
(y[ yf yJ ),
with Uk(t) E lR.m ", Yk(t) E lR. mk , k = 1,2, ... ,£, and A is such that Uk
depends only on Yk:
Uk = AkYk, Ak E'Dk .
and
'Yo(Ak) ~ 1 for all Ak E'Dk.
The scaling matrix set X also has a block-diagonal structure:
294 8. Input-Output Stability
where Uki(t) E JRmk, mki(t) E JRm\ mk = mk/£k, and each Uki depends
only on Yki with an identical linear mapping
Uki = AkYki, i = 1,2, ... , £k, Ak E ih,
where Ak is linear. The most common such linear blocks are (a) lin-
ear time-invariant system with transfer function matrix ~k(S): Uki(S) =
~k(S)Yki(S), and (b) time-varying memoryless gain matrix ~k(t): Uki(t) =
~k(t)Yki(t). Of course, linear time-varying dynamical systems also fall into
this category. Due to linearity, we can write
where the Jacobian matrix 8Uk/8Yk is symmetric and satisfies 118uk/8Ykll :::;
1. The symmetry of Jacobian matrix is satisfied by many system compo-
nents due to physical laws. The conservatism can be further reduced by
using multipliers. We will not discuss multipliers here.
When the systems discussed are linear and time invariant, it is also pos-
sible to have the scaling matrix X be frequency dependent, which allows
further reduction of conservatism. The theory of structured singular value
indicates that even with frequency-dependent scaling, there may still be
conservatism except when the number of blocks are sufficiently small.
8.5 Robust stability under dynamical uncertainty 295
x = {diag( 1m X) I XE"-"2},
Z {diag( 1m X T X ) I X E "-"2} .
If system G is described by the functional differential equation (8.44) and
(8.45) satisfying (8.46) and (8.47), then (8.70) is satisfied if the conditions
in Proposition 8.3 are satisfied for the above defined "-" and Z. To see this,
integrate (8.53) from 0 to t, considering zero initial conditions, to obtain
Io [yr(r)ZY2(r) - ur(r)Zu2(r)]dr ~ o.
t
296 8. Input-Output Stability
It It
Therefore,
IIY1{r)11 2dr ~ (I - c) lIul{r)1I2dr,
_I
t-r(t)
t ra
z{t - ra) - [Ao{r)z{r) + A 1{r)z{r - r{r))Jdr.
t-r(t)
8.6 Approximation of delay elements 297
-At(t) I t r
- ,, [Ao(r)x(r)dr
t-r(t)
+ At (r)x(r - r(r))]dr. (8.75)
Notice that we have changed notation for the state variable in (8.75). Sim-
ilar to the case discussed in Chapter 5, (8.75) is equivalent to (8.72) only
if a constraint is applied to the initial condition. Without the constraint,
the stability of (8.75) implies that of (8.72) but not vice versa because of
additional dynamics, although when rM - rm is sufficiently small, the ad-
ditional dynamics will not be unstable. The integration term in (8.75) is
considered to be the uncertainty.
System (8.75) can be written as
(8.76)
(8.77)
(8.83)
To obtain the smallest coefficient in the uncertainty term (the last term in
(8.78)), we can choose ra = (rM + r m )/2. We will show that
'Yo(A kXk ) :$ 1, for all nonsingular X k E lR nxn , k = 1,2. (8.84)
Once (8.84) is established, we can conclude the following.
We can use the method discribed in Section 8.4 to obtain an LMI con-
dition. If the coefficient matrices Ao and Al are constant matrices, we can
also use a frequency domain method.
Proof. From the above discussion, we only need to show (8.84), which
is equivalent to
Change integration variable 0 = p(T), considering (8.87) and the zero initial
condition, to obtain
(1- p) I t - r (t)
-r(O)
yf(O)ZIYI(O)q'(O)dO
rt-r(t)
(1- p) Jo yf(O)ZIYI(O)q'(B)dB
< lt yf(O)ZIYI(B)dO,
I2 = lt Uf(T)Z2U2(T)dT
12
rd
t
Jo
[I T
T-r(T)
-
ra
Y2(B)dB] T Z2 [I T
-
T-r(T)
ra
Y2(O)dB] dT.
lt l
Using Jensen's Inequality, we obtain
1 max {T-r(T),T-r a }
I2 ~ 2' Ir(T) - ral ynB)Z2Y2(O)dBdT.
rd 0 min{ T-r(T),T-r a }
[l
(Y2(O) = 0 for all B ~ 0), to obtain
1 It-min{ra,r(t)} min {Q(o),t} ]
I2 ~ 2' Ir(T) - ral dT ynO)Z2Y2(B)dO,
rd 0 min{,2(O),t}
8.6 Approximation of delay elements 299
where
l
we have
min{ma.x{6+ r",q(II)}.t}
o~ Ir(r) - raldr ~ r~.
min{min{6+r",q(II)},t}
Therefore,
x(t) = L In
N rllp - 1
where
Ap(t) = h1
p
l 1JP
9p
-
1
A(t, fJ)dO, (8.90)
lt uT(r)u(r)dr
< [r yT(r+fJ)UT(r,O)U(r,fJ)y(r+O)dO] dr
•
From the above lemma, we can immediately conclude the following.
8.6 Approximation of delay elements 301
(8.94)
1 N 18P-1
u(t) = - L (Ap(t) - A(t, B))y(t + B)dB. (8.95)
a p=l 8p
It is easy to conclude from Lemma 8.8 that 'Yo(/l.) :::; 1. Therefore, (8.93)
is sufficient for stability. •
In applying the above proposition, we may choose, for example,
N
a = r Lhpe~, (8.96)
p=l
where
ep = ~~ IIAp(t) - A(t, B)II. (8.97)
8E(8~8p-IJ
As pointed out in Section 8.4, a numerical method can be devised to check
(8.93).
When hp, p = 1,2, ... , N are small, a smaller error bound can often be
obtained by carrying out model transformation: for B E [Bp, Bp - 1 ],
t+8
x(T)dT
t+1I
+
8ap 0
-r
+ ~)d~] dT,(8.98)
where
Bap = (Bp + Bp - 1 )/2. (8.99)
Using the above in (8.91), we may conclude the following.
302 8. Input-Output Stability
Proof. Use (8.98) in (8.91), taking into account (8.90) and (8.100), to
obtain
1 N
u(t)=-73~lop t Jp
-
1
(Ap(t)-A(t,O)) [I +Oap 1_rA(T,~)Y(T+~)d~dT
t
t+9
0
] dO.
(8.103)
Therefore, the system (8.89) can be written as
y = Cu,
u = li y ,
where Ii is defined in (8.103). Define
(8.105)
- rOap
lop lop
r
(Ap(t) - A(t, O))dOz(t + T)dT] ,
8.6 Approximation of delay elements 303
1
or
0
u(t) = A2Z(t) = (31 -r H(t, r)z(t + r)dr,
where
Since
We have
1'0(.~) ~ 1'0(A2 )')'o(A 1 ) ~ l.
Therefore, (8.101) is sufficient for stability. •
To compare the two criteria, assume uniform division hp = h = r/N,
and use (8.96) to calculate 0:, then we have
When h is sufficiently small, (3 < 0:, and therefore, 1'o(G) < 1'o(G). There-
fore, Proposition 8.10 is less conservative than Proposition 8.9 when h is
sufficiently small.
where F is of bounded variation, continuous from the left on (-r, 0), and
(8.106)
h = r/N,
rk kh,
rmk (k - 1/2)h,
rMk (k + 1/2)h.
304 8. Input-Output Stability
We can write
Using
= x(t - rk) - i- rk
[[Or dF((})x(t + r + (})] dr
in (8.107), we obtain
N
:i:(t) = I: AkX(t - rk) + U(t), (8.108)
k=O
where
u(t) = - t l- r
k=O -rMk
"'k dF((})
JIJ
r-r k
Z(t + r)dr, (8.109)
and
z(t) = lOr dF((})x(t + (}).
Let
y(t) = ax(t), (8.110)
where
a = J.L2 yfrh,
and J.L is the total variation of F in [-r, 0]:
K
J.L = sup
K>O
I:
i=l
IIF((}i) - F((}i-l)1I = Var F((}).
[-r,O]
-r=IJ O <1J 1 < ... <IJK=O
y - Gu,
u 1l.y,
1
z(t) = -;; JO dF(B)y(t + B).
-r
With this framework, we will show in the following proposition that 'YO (.a) ::;
°
1, and therefore the stability is assured if 'Yo(G) < 1. An important obser-
vation is IIul12 ::; allxlb, and a -+ as h -+ 0, representing asymptotically
close approximation, which is intuitively clear.
Proof. We only need to show 'Yo(.a) ::; 1. For arbitrary interval [a, b]
and scalar function (3( B) 2:: 0, B E [a, b], define
1 b
(3(B)ldF(B)1 = !~
a=8o<8 1 <···<8K=b
8K
(3(Bi)IIF(Bi) - F(Bi-1)11·
The special case of (3(B) = 1 is the total variation of F(B) in [a, b]:
I a
b
IdF(B)1 = Var F(B),
[a,~
especially
I: IdF(B) I = Jj.
With this definition, we have another form of Jensen's Inequality (see [245]):
306 8. Input-Output Stability
We can write
Therefore
or
'Y~(al)
r
= ~i~~~k IIT(r)11 2dr
N h/2
= ~ih/2"T(r - kh)1I 2dr
r
h/2 N
= ih/2 ~ IIT(r - kh)1I 2dr
< ih/2
h/2 C
~ IIT(r - kh)1I dr
jh/2
< J.L 2dr
-h/2
hJ.L 2.
Therefore
'Yo(a 1) S JrJ.L2h.
We obtain
•
8.7 Notes
The input-output approach to stability and the small gain theorem are
discussed in more details in Desoer and Vidyasagar [54J and Vidyasagar
[277J, where more general settings are used.
The comparison system approach were explored by, for example, Halanay
[101], Lakshmikantham and Leela [164J, Driver [59J, and Huang and Zhou
[118J. The material presented here is a generalization of Zhang, Knopse,
and Tsiotras [292J. In [292]' C is constrained to AIM and the frequency
domain method was used. Equivalence was established with the stability
criterion proposed in Park [220J (see also the Notes in Chapter 5).
The scaled small gain problem and related passivity problem for finite-
dimensional linear time-invariant systems have been discussed in detail
in, for example, Boyd and Yang [22J and Packard and Doyle [219], with
308 8. Input-Output Stability
Matrix Facts
A.I Notation
In this appendix, we will discuss some useful facts about matrices. We
assume the readers are familiar with the basic matrix algebra. The material
presented here is not intended to be self-contained. Rather, it is intended
to fix the notation and point out some facts that may not be familiar to
some readers. For a systematic presentation of these materials, the readers
are referred to the references discussed in the Notes section at the end of
this appendix.
We use Rn to represent the set of real column vectors with n components,
and en of complex column vectors. Similarly, Rmxn and e mxn represent
the sets of real and complex m by n matrices, respectively. We will use Xi
to denote the ith component of column vector X, and aij to represent the
element of matrix A at the ith row and jth column. AT is the transpose
of A. For complex matrices, the Hermitian transpose defined as A H ~ AT
is often used, where A is the complex conjugate of A. In the following, we
will mainly discuss complex matrices, but with obvious specialization to
real matrices. I denotes an identity matrix. If we need to emphasize the
dimention, we may write, for example, In for the n-dimensional identity
matrix. The rank of A, the number of linearly independent columns, is
denoted as rank(A). If rank(A) = mini m, n}, then we say A is full rank.
We say A is of full row rank if rank (A) = m, and full column rank if
rank(A) = n. For a square matrix A E e nxn , we also use tr(A) to represent
the trace of A, i.e.,
n
tr(A) = L aii·
i=t
A.2 Determinant
For a square matrix A E cnxn , its determinant is denoted det (A). If A
and B are both square matrices, then it can be shown that
to obtain
det( ~ ~) = det(A) det(D - CA- l B), (A.3)
often use Ai(A) to denote the ith eigenvalue of A. The set of all the eigen-
values of matrix A is known as the spectrum of A and is denoted O'(A).
The spectrum radius is defined as
p(A) = m!1X IAi(A)I.
l~.~n
We sometimes write 0" max and 0" min as (j and Q, respectively. Singular values
have the following relations with the eigenvalues
A.5 Norms
A norm II . II is a function defined on en satisfying the following three
conditions
1. °
Ilxll = if and only if x = 0,
2. Ilaxll = lal '11xll for any scalar a,
3. Ilx + yll s; Ilxll + Ilyll·
For x E en, let
x = (Xl X2 ... xn )T .
1. I-norm:
n
Ilxll! = L IXil.
i=l
A.5 Norms 313
2. 2-norm:
n
Ilxlb = L
i=l
I iI 2, X
3. oo-norm
The 2-norm is geometrically the most appealing, but I-norm and 00-
norms are sometimes more convenient to use. Obviously, an m by n matrix
can be regarded as a vector with m x n components, and the norms can
be defined accordingly. For example, corresponding to the 2-norm, we can
define the following norm of a matrix A E mxn c
m n
IIAIIF = LL laijl2.
i=l j=l
However, the most useful norms for matrices are the induced norms defined
in terms of the corresponding vector norm
IIAxlik
IIAllk = IIxllk9
max IIAxllk = max - -
IIxllk;6o IIxlik .
It can be shown that the three most common induced norms are
1. I-norm
m
2. 2-norm
3. OQ-norm
allB
a a21 B
A®B= ( .
am1B
For square matrices A E c mxm , B E cnxn , we can also define their Kro-
necker sum as
If Ai, i = 1,2, ... , m are the eigenvalues of A, and Itj' j = 1,2, ... , n are
the eigenvalues of B, then the mn eigenvalues of A®B and AE9B are Ailtj,
and Ai + Itj' i = 1,2, ... , mj j = 1,2, ... , n, respectively.
A.8 Notes
Most of the materials presented here can be found in Golub and Van Loan
[76] and Horn and Johnson [114]. The Kronecker product and sum are
covered in Graham [79].
Appendix B
where
x = (Xl X2 ... Xm) TERm
F(1)(x)
( O °
F(2)(x)
316 Appendix B. LMI and Quadratic Integral Inequalities
For the sake of convenience, the variables often appear in matrix form.
For example,
(B.2)
(B.3)
implies
a > 0, d > 0, f > 0,
and
B.2 GEVP
The GEVP is to find
inf ,x, (B.4)
subject to the constraint
where A(x), B(x), and C(x) are matrices depending on the variable x E lRm
linearly. As a special case, the constraint (B.7) may be absent. The termi-
nology is motivated by the generalized eigenvalue problem in the classical
eigenvalue theory of matrices, which is also very useful in this book. To see
the connection, consider the special case of the above problem. When A
and B are constant matrices, and the additional constraint (B. 7) is absent,
B.2 GEVP 317
then the solution of the problem is the greatest >. to satisfy the following
generalized eigenvalue problem
A~ = >'B~, ~ # 0.
We will refer to the problem described by (B.4) to (B.7) by the abbreviation
GEVP in order to distinguish it from the classical generalized eigenvalue
problem.
A number of problems encountered in this book may be conveniently
transformed into the standard GEVP. For example, consider the problem
a = sup a, (B.8)
subject to
P(x) + aQ(x) > 0, (B.9)
P(x) > 0, (B.IO)
R(x) < 0. (B.ll)
This problem can be transformed to the standard GEVP
~a = inf,6,
subject to
,6P(x) + Q(x) > 0,
P(x) > 0,
-R(x) > 0.
Another problem that can be easily transformed to the standard GEVP is
a = sup a,
subject to
P(x) + aQ(x) < 0,
Q(x) > 0.
This is equivalent to the standard GEVP
~a = inf,6,
subject to
-,6P(x) - Q(x) > 0,
-P(x) > 0,
Q(x) > 0.
We will simply refer to these problems also as GEVP without further com-
ments.
318 Appendix B. LMI and Quadratic Integral Inequalities
(B.14)
A > 0, (B.15)
C - BT A-I B > O. (B.16)
B.4 S-procedure
S-procedure plays important role in robust stability theory. It can be stated
as follows.
Proposition B.3 Let Fi E IR nxn , i = 0,1,2, ... ,p. Then the following
statement is true
(B.19)
if and only if
(t T ~) > 0, (B.20)
Notice that (B.20) and (B.2I) are obtained by deleting the row and
column containing X in two different ways.
Proof. Necessity is obvious since the left hand side of both (B.20)
and (B.2I) are principal minors of the left hand side of (B.I9). For suf-
ficiency, left-multiply second row of (B.I9) by - V T R- 1 and add to the
320 Appendix B. LMI and Quadratic Integral Inequalities
third row, with symmetric operation for the columns. Equivalently, left-
multiply (B.19) by
( o~ ~_VTR- ~),
I I
n
Corollary B.5 There exists a matrix X such that
( P+XEo+ElxT Q+XE
(Q+xEf R > 0, (B.23)
XT VT
if and only if
( P+E,{;SEo
MT ~) > 0, (B.24)
(~ ~) > 0, (B.25)
where
M Q - Eifv T + Eif SEl,
N R- VE-ETVT +ETSE.
Proof. Left-multiply the third row of the left hand side of (B.23) by
- (Eo E) T and add to the second row, with symmetric operation on
the columns, then use Proposition BA, to realize that (B.23) is equivalent
to (B.24) and
R- VE _ETVT +ETSE V - ETS )
( (V - ET Sf S > O. (B.26)
( P2 +X Q2
QrR2 ) >
° (B.28)
if and only if
( PI + LP,LT Ql LQ2
°°
Qr Rl ) >0 (B.29)
QrLT R2
and
.bo = PI + LP2L T - QIRIIQf - LQ2R2 1 Qf LT > 0. (B.32)
Also, (B.27) and (B.28) are equivalent to (B.30), (B.31), and
.bo 1 PI - LXLT - QIRIIQf > 0, (B.33)
.bo2 = P2 + X - Q2R;;lQf > 0. (B.34)
Therefore, we only need to show that, given (B.30) and (B.31), there exists
an X satisfying (B.33) and (B.34) if and only if (B.32) is satisfied.
The existence of X to satisfy (B.33) and (B.34) implies (B.32) since
.bo = .bo 1 + L.bo 2LT. On the other hand, if (B.32) is satisfied, let
X = Q2R21Qr - P2 +d
for some sufficiently small E °
> results in
.bo2 = cI > 0,
and
.bo 1 = .bo - L.bo 2L T = .bo - ELLT > 0,
thus (B.33) and (B.34) are satisfied. •
The above results (Proposition BA, Corollary B.5, and Proposition B.6)
can be extended to the case where the constant matrices are replaced by
continuous matrix functions depending on variables within a compact set.
For example, Proposition BA can be extended to the following.
322 Appendix B. LMI and Quadratic Integral Inequalities
Corollary B.7 For given continuous matrix functions P( 0), Q( 0), R( 0),
V(O), and S(O), with variable 0 varying with a compact set 0 C lR.k, there
exists a continuous matrix function X(O) such that
R(O) V(O) )
(
VT(O) S(O)
> 0, (B.37)
for all 0 E o.
Proof. Based on Proposition B.4 and its proof, it is sufficient to prove
that X as expressed in (B.22) is continuous. From (B.36), R > 0 for all
o E O. This implies that R > cI for some constant c: > 0 since R is
continuous and 0 is compact. Therefore, R- 1 exists and is continuous and
bounded. Since Q and V are also continuous, X is continuous. •
> 11 (1 0
WI (f3)df3 + i1 W2(f3)df3 ) T
M (1 0
WI (f3)df3 +i 1 w2(f3)df3 ) da, (B.39)
11 (1 - a)wT(a)Mw(a)da
> 11 (1 0
W(f3)d f3 ) T ±M (1 0
W(f3)d f3 ) da (BAD)
> 11 (1 0
W(f3)df3 ) T M (1 0
w(f3)df3 ) da, (BA1)
11 awT(a)Mw(a)da
> (i
11 1w(f3)df3) T 1 ~ aM (i1 w(f3)df3) da (BA2)
1 0
wf(f3)Mwl(f3)df3 + i 1
wf(f3)Mw2(f3)df3
11 (1- a)wT(a)Mw(a)da
11 11 d(3wT(a)Mw(a)da
11 1f3 wT(a)Mw(a)dad(3
> l (J.' W(a)da) T ~M (J.' W(a)da) d~,
which is (BAO). In the above, the equality in the second step is the result
of exchanging the order of integration, and the inequality in the last step
results from the Jensen Inequality (B.38). Considering the fact that l/a 2': 1
in the integration interval, (BA1) is obvious. The proof of (BA2) and (BA3)
are similar. •
B.7 Notes
The importance of matrix inequalities, especially LMIs, in systems theory
and control has long been recognized since Lyapunov theory was published
around 1890 [185]. Indeed, solutions of some matrix inequalities have ap-
peared as early as early 1960s by Yakubovich [287]-[289], see also the work
by Horisberger and Belanger [115]. Largely because of the unavailability of
efficient numerical algorithms for the general form of LMI, most of the ear-
lier works on problems related to LMIs were reformulated in forms such as
Lyapunov equations and Algebraic Riccati Equations (AREs). The realiza-
tion that LMI is a convex optimization problem [20], and the development
of the efficient interior point method [203] have spurred tremendous inter-
est among researchers in control systems theory to formulate many control
problems in LMI form. See the book by Boyd et al. [21] for more concepts
on LMI, GEVP, and applications in system and control theory. For more
recent progress, and generalization to semidefinite programming, see [266]
and [75]. Commercial software of an LMI solver in MAT LAB is available
[72].
For a more comprehensive and systematic coverage of the materials in
Sections B.l to B.3, see the book by Boyd et al. [21].
The elimination of variables in LMI in Section B.5 was discussed in [87]
and [89]. Some other cases of variable elimination are discussed in [21].
B.7 Notes 325
[26J M. Briggs and R. Vinter. Linear filtering for time-delay systems. IMA
J. Math. Control Info., 6:167-178, 1989.
[28J Y.-Y. Cao, Y.-X. Sun, and C. Cheng. Delay-dependent robust stabi-
lization of uncertain systems with multiple state delays. IEEE 7rans.
Auto. Control, 43(11):1608-1612, 1998.
[29J Y.-Y. Cao and Y.-X. Sun. Robust stabilization of uncertain sys-
tems with time-varying multistate delays. IEEE 7rans. Auto. Con-
trol, 43(10):1484-1488, 1998.
[31J V.-S. Chellaboina. Robust stability and performance for linear and
nonlinear uncertain systems with structured uncertainty. Ph.D. Dis-
sertation, School of Aerospace Engineering, Georgia Institute of Tech-
nology, 1996.
[34J J. Chen, G. Gu, and C. N. Nett. A new method for computing delay
margins for stability of linear delay systems. Systems and Control
Letters, 26(2):107-117, Sept. 1995.
[51J R. Datko. Not all feedback stabilized hyperbolic systems are robust
with respect to small time delays in their feedbacks. SIAM J. Control
Optim., 26:697-713, 1988.
[68] M. Fu, H. Li, and S.-I. Niculescu. Robust stability and stabilization
of time-delay systems via integral quadratic constraint approach. In
L. Dugard and E. I. Verriest (Eds.), Stability and Control of Time-
Delay Systems. Lecture Notes in Control and Information Sciences,
228, Springer-Verlag, New York, 101-116, 1997.
[69] M. Fu. The real J.1, is hardly approximable. IEEE Trans. Auto. Con-
trol, AC-42:1286-1288, 1997.
[74] J.-H. Ge, P. M. Frank, and C.-F. Lin. Hoc control via output feedback
for state delayed systems. Int. J. Control, 64(1):1-8, 1996.
[116) C. Hou and J. Qian. On an estimate of the decay rate for applications
of Razumikhin-type theorem. IEEE Trans. Auto. Control, 43(7):958-
960, 1998.
[146] J.-H. Kim. Delay and its time-derivative dependent robust stability
of time-delayed linear systems with uncertainty. IEEE Trans. Auto.
Control, 46(5):789-792, 2001.
[147] J.-H. Kim. Robust stability of linear systems with delayed perturba-
tions. IEEE Trans. Auto. Control, 41(12):1820-1822, 1996.
[167] E. B. Lee and S. Zak. On spectrum placement for linear time invariant
delay systems. IEEE Trans. Auto. Control, AC-27(2):446-449, 1982.
340 Bibliography
[184] Z.-H. Luo, B.-Z. Guo, and O. Morgul. Stability and Stabilization of In-
finite Dimensional Systems with Applications. Springer-Verlag, Lon-
don, 1999.
[210] S.-I. Niculescu, J.-M. Dion, L. Dugard, and H. Li. Stability of linear
systems with delayed state: an L.M.I. approach. European Journal
of Automation, special issue on analysis and control of time-delay
systems, 31:955-970, 1997.
[244J J.-C. Shen, B.-S. Chen, and F.-C. Kung. Memoryless stabilization of
uncertain dynamic delay systems: Riccati equation approach. IEEE
Trans. Auto. Control, 36(5):638-640, 1991.
[262] O. Toker and H. Ozbay. On the complexity of purely complex J.L com-
putation and related problems in multidimensional systems. IEEE
Trans. Auto. Control, AC-43:409-414, 1998.
[270] E.1. Verriest. Stability of systems with distributed delays. IFAC Sys-
tem Structure and Control, Nantes, France, 294-299, 1995.
[284] M. Wirth and U. Brunner. Observer models for real-time traffic con-
trol. Proceedings of the American Control Conference, Baltimore,
MD, 923-926, 1994.
[285] H. Wu and K. Mizukami. Linear and nonlinear stabilizing continuous
controllers of uncertain dynamical systems including state delays.
IEEE Trans. Auto. Control, 41(1):116--121, 1996.
[286] H. Wu and K. Mizukami. Robust stability criteria for dynamical sys-
tems including delayed perturbation. IEEE Trans. Auto. Control,
40(3):487-490, 1995.
Bibliography 349
lI'-Knapsack problem, 90
time-delay, see delay
time-delay system, 9
linear, 15
linear time-invariant, 17
LTI, see linear
time-invariant
neutral, 20
triangle inequality, 99
trivial solution, 10
Tsypkin theorem, 40
Tsypkin's test, 75
two-variable criterion, 34, 82
Series Editor
William S. Levine
Department of Electrical and Computer Engineering
University of Maryland
College Park, MD 20742-3285
USA
The Birkhauser series Systems and Control: Foundations and Applications examines
the abstract and theoretical mathematical aspects of control. Control Engineering
complements this effort through a study of the industrial and applied implementation of
control- from techniques for analysis and design to hardware implementation, test, and
evaluation. While recognizing the harmony between abstract theory and physical applica-
tion, these publications emphasize real-world results and concerns. Problems and ex-
amples use the least amount of abstraction required, remaining committed to issues of
consequence, such as cost, tradeoffs, reliability, and power consumption.
Readership
Preparation of manuscripts