Lecture Notes On Linear System Theory
Lecture Notes On Linear System Theory
Lecture Notes On Linear System Theory
Automatic Control Laboratory, ETH Zurich
CH-8092, Zurich, Switzerland
[email protected]
1 Introduction 1
1.1 Objectives of the course . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Proof methods . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Functions and maps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2 Introduction to Algebra 11
2.1 Groups . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2 Rings and fields . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.3 Linear spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.4 Subspaces and bases . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.5 Linear maps . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.6 Linear maps generated by matrices . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
2.7 Matrix representation of linear maps . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.8 Change of basis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
3 Introduction to Analysis 33
3.1 Norms and continuity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.2 Infinite-dimensional normed spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.3 Completeness . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
3.4 The space of square-summable functions . . . . . . . . . . . . . . . . . . . . . . . . . 42
3.5 Induced norms and matrix norms . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
3.6 Ordinary differential equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
3.7 Existence and uniqueness of solutions . . . . . . . . . . . . . . . . . . . . . . . . . . 50
3.7.1 Background lemmas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
3.7.2 Proof of existence . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.7.3 Proof of uniqueness . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
i
4.2 Existence and structure of solutions . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
4.3 State transition matrix . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
6 Stability 83
6.1 Nonlinear systems: Basic definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
6.2 Linear time varying systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
6.3 Linear time invariant systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91
6.4 Systems with inputs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
6.5 Lyapunov equation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 95
ii
A Class composition and outline 149
A.1 Class composition . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 149
A.2 Class outline . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 149
iii
Chapter 1
Introduction
where x(t) Rn denotes the system state, u(t) Rm denotes the system inputs, y(t) Rp denotes
the system outputs, A(t) Rnn , B(t) Rnm , C(t) Rpn , and D(t) Rpm are matrices of
appropriate dimensions, and where, as usual, x(t) = dx
dt (t) denotes the derivative of x(t) with respect
to time.
Time varying linear systems are useful in many application areas. They frequently arise as models
of mechanical or electrical systems whose parameters (for example, the stiffness of a spring or the
inductance of a coil) change in time. As we will see, time varying linear systems also arise when
one linearizes a non-linear system around a trajectory. This is very common in practice. Faced with
a nonlinear system one often uses the full nonlinear dynamics to design an optimal trajectory to
guide the system from its initial state to a desired final state. However, ensuring that the system
will actually track this trajectory in the presence of disturbances is not an easy task. One solution
is to linearize the nonlinear system (i.e. approximate it by a linear system) around the optimal
trajectory; the approximation is accurate as long as the nonlinear system does not drift too far away
from the optimal trajectory. The result of the linearization is a time varying linear system, which
can be controlled using the methods developed in this course. If the control design is done well, the
state of the nonlinear system will always stay close to the optimal trajectory, hence ensuring that
the linear approximation remains valid.
A special class of linear time varying systems are linear time invariant systems, usually referred to
by the acronym LTI. LTI systems are described by state equations of the form
where the matrices A Rnn , B Rnm , C Rpn , and D Rpm are constant for all times
t R+ . LTI systems are somewhat easier to deal with and will be treated in the course as a special
case of the more general linear time varying systems.
1
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
2
The second and less obvious objective of the course is for students to experience something about
doing automatic control research, in particular developing mathematical proofs and formal logical
arguments. Linear systems are ideally suited for this task. There are two main reasons for this. The
first is that almost all the derivations given in the class can be carried out in complete detail, down to
the level of basic algebra. There are very few places where one has to invoke higher powers, such as
an obscure mathematical theorem whose proof is outside the scope of the course. One can generally
continue the calculations until he/she is convinced that the claim is true. The second reason is that
linear systems theory brings together two areas of mathematics, algebra and analysis. As we will
soon see, the state space, Rn , of the systems has both an algebraic structure (it is a vector space)
and a topological structure (it is a normed space). The algebraic structure allows us to perform
linear algebra operations, compute projections, eigenvalues, etc. The topological structure, on the
other hand, forms the basis of analysis, the definition of derivatives, etc. The main point of linear
systems theory is to exploit the algebraic structure to develop tractable algorithms that allow us
to answer analysis questions which appear intractable by themselves.
For example, consider the time invariant linear system
with x(t) Rn , u(t) Rm , A Rnn , and B Rnm . Given x0 Rn , T > 0 and a continuous
function u() : [0, T ] Rm (known as the input trajectory) one can show (Chapter 3) that there
exists a unique function x() : [0, T ] Rn such that
This function is called the state trajectory (or simply the solution) of system (1.1) with initial
condition x0 under the input u(). As we will see in Chapter 3, u() does not even need to be
continuous for (1.2) to be true, provided one appropriately qualifies the statement for all t [0, T ].
System (1.1) is called controllable (Chapter 8) if and only if for all x0 Rn , for all x Rn , and for
all T > 0, there exists u() : [0, T ] Rm such that the solution of system (1.1) with initial condition
x0 under the input u() is such that x(T ) = x. Controllability is clearly an interesting property for
a system to have. If the system is controllable then we can guide it from any initial state to any
final state by selecting an appropriate input. If not, there may be some desirable parts of the state
space that we cannot reach from some initial states. Unfortunately, determining whether a system is
controllable directly from the definition is impossible. This would require calculating all trajectories
that start at all initial conditions. Except for trivial cases (like the linear system x(t) = u(t)) this
calculation is intractable, since the initial states, x0 , the times T of interest, and the possible input
trajectories u() : [0, T ] Rm are all infinite. Fortunately, linear algebra can be used to answer the
question without even computing a single solution (Chapter 8).
Theorem 1.1 System (1.1) is controllable if and only if the matrix B AB . . . An1 B Rnnm
has rank n.
The theorem shows how the seemingly intractable analysis question is the system (1.1) control-
lable? can be answered by a simple algebraic calculation of the rank of a matrix.
either because they have been proven before, or because we postulate so (for example the axioms
of logic), or because we assume so in a certain context (for example, when we say Let V be a vector
space . . . we mean Assume that the set V verifies the axioms of a vector space . . . ).
Theorems of minor importance, or theorems whose main point is to establish an intermediate step in
the proof of another theorem, will be called Lemmas, Facts, or Propositions; An immediate
consequence of a theorem that deserves to be highlighted separately is usually called a Corollary.
And a logical statement that we think may be true but cannot prove so is called a Conjecture.
The logical statements we will most be interested in typically take the form
pq
Example (No smoke without fire) It is generally accepted that when there is smoke, there must
be some a fire somewhere. This knowledge can be encoded by the logical implication
This is a statement of the form p q with p the statement there is smoke and q the statement
there is a fire.
Hypotheses and consequences may typically depend on one or more free variables, that is, objects
that in the formulation of hypotheses and consequences are left free to change.
Example (Greeks) Despite recent economic turbulence, it is generally accepted that Greek citizens
are also Europeans. This knowledge can be encoded by the logical implication
A sentence like X is a . . . is the verbal way of saying something belongs to a set; for example the
above statement can also be written as
X Greeks X Europeans,
where Greeks and Europeans are supposed to be sets; the assertion that this implication is true
for arbitrary X (X, X Greeks X Europeans) is equivalent to the set-theoretic statement of
inclusion:
Greeks Europeans.
You can visualize the implication and its set-theoretic interpretation in Figure 1.1.
There are several ways of proving that logical statements are true. The most obvious one is a direct
proof: Start from p and establish a finite sequence of intermediate implications, p1 , p2 , . . . , pn
leading to q
p p1 p2 . . . pn q.
We illustrate this proof technique using a statement about the natural numbers.
Definition 1.1 A natural number n N is called odd if and only if there exists k N such that
n = 2k + 1. It is called even if and only if there exists k N such that n = 2k.
One can indeed show that all natural numbers are either even or odd, and no natural number is
both even and odd (Problem 1.1).
People
Europeans
Greeks
Proof:
n is odd k N : n = 2k + 1
k N : n2 = (2k + 1)(2k + 1)
k N : n2 = 4k 2 + 4k + 1
k N : n2 = 2(2k 2 + 2k) + 1
l N : n2 = 2l + 1 (namely, l = 2k 2 + 2k N)
n2 is odd
This proof principle can also be exploited to perform proof by induction. Proof by induction concerns
propositions, pk , indexed by the natural numbers, k N, and statements of the form
k N, pk is true.
One often proves such statements by showing that p0 is true and then establishing and infinite
sequence of implications
p0 p1 p2 . . . .
Clearly proving these implications one by one is impractical. It suffices, however, to establish that
pk pk+1 for all k N, or in other words
[p0 (pk pk+1 , k N)] [pk , k N] .
We demonstrate this proof style with another statement about the natural numbers.
Definition 1.2 The factorial, n!, of a natural number n N is the natural number n! = n (n 1)
. . . 2 1. By convention, if n = 0 we set n! = 1.
Proof: It is easy to check that the statement holds for the special cases m = k = 0, m = 0 and
k = 1, and m = 1 and k = 0. For the case m = k = 1, (m + k)! = 2! 1!1! = m!k!.
Assume now that for some m, k N, (m + k)! m!k! (we call this the induction hypothesis) and
show that for m, k + 1, (m + k + 1)! m!(k + 1)! (also (m + 1)!k! for the case m + 1, k, by symmetry).
(m + k + 1)! = (m + k)!(m + k + 1)
m!k!(m + k + 1) (by the induction hypothesis)
m!k!(k + 1) (since m N)
= m!(k + 1)!
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
5
Even though there is no direct way to illustrate proof by induction using statements about Greeks
and other nationalities, one could in principle use a similar line of reasoning to prove a statement
like once a Greek always a Greek by arguing that children with at least one Greek parent are
themselves Greek.
Sometimes direct proof p q is difficult. In this case we try to find other statements that are
logically equivalent to p q and prove these instead. An example of such a statement is q p,
or in logic notation
(q p) (p q).
Example (Greeks (cont.)) The statement that all the Greeks are Europeans is also equivalent to
non-Europeans non-Greeks.
If we stipulate a priori that all the possible X we may consider in our discourse belong to some big
set (i.e., People), in fact this is also equivalent to
where \ denotes the difference of two sets and the superscript c denotes the set-complement with
respect to People.
Exercise 1.1 Visualize this set theoretic interpretation by a picture similar to Figure 1.1.
Proof: Let p = n2 is odd, q = n is odd. Assume n is even (q) and show that n2 is even (p).
n is even k N : n = 2k
k N : n2 = 4k 2
l N : n = 2l (namely, l = 2k 2 N)
n2 is even.
Another common method that can be used to indirectly prove that p q is to suppose that p
is true, to suppose that q is false, and to apply other proof methods to derive a contradiction. A
contradiction is a proposition of the form r r (like There is smoke and there is no smoke, or
n is even and n is odd); all such statements are postulated to be false by virtue of their mere
structure, and irrespective of the proposition r. If, by assuming p is true and q is false we are able
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
6
to reach a false assertion, we must admit that if p is true the consequence q cannot be false, in other
words that p implies q. This method is known as proof by contradiction.
Example (Greeks and Chinese) Suppose the following implications: for all X,
X is a Greek X is a European
X is a Chinese X is an Asian
X is an Asian X is not a European
Indeed, suppose p(X) and the converse of q(X), that is, X is a Chinese. By direct deduction,
X is a Greek X is a Chinese
X is a European X is an Asian
X is a European X is not a European
Since the conclusion is a contradiction for all X, we must admit that p(X) q(X) is false or, which
is the same, that p(X) q(X). The set-theoretic interpretation is as follows: By postulate,
Europeans non-Europeans =
It follows that Greeks Chinese is also equal to the empty set. Therefore (here is the point of the
above proof), Greeks non-Chinese.
Exercise 1.2 Visualize this set theoretic interpretation by a picture similar to Figure 1.1.
We will illustrate this fundamental proof technique with another statement, about rational numbers.
Definition 1.3 The real number x R is called rational if and only if there exist integers n, m Z
with m 6= 0 such that x = n/m.
Theorem 1.5 (Pythagoras) 2 is not rational.
Proof: (Euclid) Assume,
for the sake of contradiction,
that 2 is rational. Then there exist n, m Z
with m 6= 0 such that 2 = n/m. Since 2 > 0, without loss of generality we can take n, m N; if
they happen to be both negative multiply both by 1 and replace them by the resulting numbers.
Without loss of generality, we can further assume that m and n have no common divisor; if they do,
divide both by their common divisors until there are no common divisors left and replace m and n
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
7
n2 is even
n is even (Theorem 1.4 and Problem 1.1)
k N : n = 2k
k N : 2m2 = n2 = 4k 2
k N : m2 = 2k 2
m2 is even
m is even (Theorem 1.4 and Problem 1.1).
Exercise 1.3 What is the statement p in Theorem 1.5? What is the statement q? What is the
statement r in the logical contradiction r r reached at the end of the proof?
Two statements are equivalent if one implies the other and vice versa,
(p q) is the same as (p q) (q p)
Usually showing that two statements are equivalent is done in two steps: Show that p q and then
show that q p. For example, consider the following statement about the natural numbers.
Proof: n is odd implies that n2 is odd (by Theorem 1.2) and n2 is odd implies that n is odd (by
Theorem 1.4). Therefore the two statements are equivalent.
This is argument is related to the canonical way of proving that two sets are equal, by proving two
set inclusions A B and B A. To prove these inclusions one proves two implications:
X AXB
XBXA
Example (Maximal natural number) Here is a proof that there is no number larger than 1.
Proof: Assume, for the sake of contradiction, that N > 1. Then N 2 is also a natural number and
N 2 > N . This contradicts the fact that N is the largest natural number. Therefore we must have
N = 1.
Obviously the theorem in this example is saying something quite silly. The problem, however,
is not that the proof is incorrect, but that the starting hypothesis let N be the largest natural
number is false, since there is no largest natural number.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
8
g
X Y
f g f
{y Y | x X : f (x) = y} Y
{(x, y) X Y | y = f (x)} X Y
Exercise 1.4 Show that composition is associative. In other words, for any three functions g : X
Y , f : Y Z and h : W X and for all w W , f (g h)(w) = (f g) h(w).
By virtue of this associativity property, we will simply use f gh : W Y to denote the composition
of three (or more) functions.
A special function that can always be defined on any set is the identity function, also called the
identity map, or simply the identity.
Definition 1.5 The identity map on X is the function 1X : X X defined by 1X (x) = x for all
x X.
Using the identity map one can also define various inverses of functions.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
9
f f
g L f = 1X X Y X Y f g R = 1Y
gL gR
f
g f = 1X X Y f g = 1Y
g
Figure 1.3: Commutative diagram of function inverses.
The commutative diagrams for the different types of inverses are shown in Figure 1.3.
Proof: Assume, for the sake of contradiction, that f is invertible but there exist two different
inverses, g1 : Y X and g2 : Y X. Since the inverses are different, there must exist y Y such
that g1 (y) 6= g2 (y). Let x1 = g1 (y) and x2 = g2 (y) and note that x1 6= x2 . Then
Problem 1.2 (Inverses of functions) Consider two sets X and Y and a function f : X Y .
Show that:
3. f is invertible if and only if it is bijective. In this case show also that all inverses (left-, right-
and two-sided) coincide.
4. If f is invertible then any two inverses (left-, right- or both) coincide.
Chapter 2
Introduction to Algebra
2.1 Groups
Definition 2.1 A group (G, ) is a set G equipped with a binary operation : G G G such
that:
(G, ) is called commutative (or Abelian) if and only if in addition to 1-3 above
4. is commutative: a, b G, a b = b a.
with the usual operation of matrix multiplication is a group1 . What is the identity? What is the
inverse?
Fact 2.1 For a group (G, ) the identity element, e, is unique. Moreover, for all a G the inverse
element, a1 , is unique.
1 For the time being, the reader is asked to excuse the use of matrices in the examples. Matrices will be formally
defined in the next section, but will be used in the meantime for informal illustrations.
11
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
12
Proof: To show the first statement, assume, for the sake of contradiction, that there exist two
identity elements e, e G with e 6= e . Then for all a G, e a = a e = a and e a = a e = a.
Then:
e = e e = e
which contradicts the assumption that e 6= e .
To show the second statement, assume, for the sake of contradiction, that there exists a G with
two inverse elements, say a1 and a2 with a1 6= a2 . Then
a1 = a1 e = a1 (a a2 ) = (a1 a) a2 = e a2 = a2 ,
am sm + am1 sm1 + . . . + a0
bn sn + bn1 sn1 + . . . + b0
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
13
an sn + an1 sn1 + . . . + a0
bn sn + bn1 sn1 + . . . + b0
for some n N with a0 , . . . , an , b0 , . . . , bn R with bn 6= 0 is a commutative ring. Note that an = 0
is allowed, i.e. it is possible for the degree of the numerator polynomial to be less than or equal to
that of the denominator polynomial.
Exercise 2.1 Show that for every ring (R, +, ) the identity elements 0 and 1 are unique. Moreover,
for all a R the inverse element (a) is unique.
1. For all a R, a 0 = 0 a = 0.
2. For all a, b R, (a) b = (a b) = a (b).
a + 0 = a a (a + 0) = a a a a + a 0 = a a
(a a) + a a + a 0 = (a a) + a a
0 + a 0 = 0 a 0 = 0.
The second equation is similar. For the second statement note that
(R[s], +, ) is not a field, since the multiplicative inverse of a polynomial is not a polynomial but a
rational function.
(R(s), +, ) is a field.
(Rp (s), +, ) is not a field, since the multiplicative inverse of a proper rational function is not neces-
sarily proper.
Exercise 2.2 If (F, +, ) is a field then for all a R the multiplicative inverse element a1 is unique.
Exercise 2.3 Show that if (F, +, ) is a field and a 6= 0 then a b = a c a = c. Is the same true
for a ring? Illustrate using an example from R22 .
with a11 , . . . , anm F . One can then define the usual matrix operations in the usual way: ma-
trix multiplication, matrix addition, determinant, identity matrix, inverse matrix, adjoint matrix,
. . . . We assume that the reader is familiar with these operations from basic linear algebra, a brief
summary is provided in Appendix C.
Matrices can also be formed from the elements of a commutative ring (R, +, ).
Example (Transfer function matrices) As we will see in Chapter 5, the set of transfer functions
of time invariant linear systems with m inputs and p outputs are p m matrices whose elements are
proper rational functions,
x(t) = Ax(t) + Bu(t)
G(s) = C(sI A)1 B + D Rp (s)pm .
y(t) = Cx(t) + Du(t)
Fact 2.3 Assume that (R, +, ) is a commutative ring, A Rnn , and Det(A) 6= 0. It is not always
the case that A1 exists.
Roughly speaking, what goes wrong is that even though DetA 6= 0, the inverse (DetA)1 may not
exist in R. Then the inverse of the matrix defined as A1 = (Det(A))1 Adj(A) is not defined. The
example of the transfer function matrices given above illustrates the point: Assuming m = p, the
elements of the inverse of a transfer function matrix are not necessarily proper rational functions,
and hence the inverse matrix cannot be the transfer function of a linear time invariant system.
Definition 2.4 A linear space (V, F, , ) is a set V (of vectors) and a field (F, +, ) (of scalars)
equipped with two binary operations, : V V V (called vector addition) and : F V V
(called scalar multiplication) such that:
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
15
Exercise 2.4 Let (F, +, ) be a field. Show that (F, F, +, ) is a linear space.
Exercise 2.5 For every linear space (V, F, , ) the identity element is unique. Moreover, for all
x V there exists a unique inverse element x.
Fact 2.4 If (V, F, , ) is a linear space and 0 is the addition identity element of F then for all
x V , 0 x = . Moreover, for all a F , x V , (a) x = (a x) = a (x).
Definition 2.5 If (V, F, V , V ) and (W, F, W , W ) are linear spaces over the same field, the
product space (V W, F, , ) is the linear space comprising all pairs (v, w) V W with defined
by (v1 , w1 ) (v2 , w2 ) = (v1 V v2 , w1 W w2 ), and defined by a (v, w) = (a V v, a W w).
Exercise 2.6 Show that (V W, F, , ) is a linear space. What is the identity element for addition?
What is the inverse element?
Two types of linear spaces will play a central role in these notes. The first is constructed by taking
repeatedly the product of a field with itself.
Example (Finite product of a field) For any field, (F, +, ), consider the product space F n . Let
x = (x1 , . . . , xn ) F n , y = (y1 , . . . , yn ) F n and a F and define : F n F n F n by
x y = (x1 + y1 , . . . , xn + yn )
and : F F n F n by
a x = (a x1 , . . . , a xn ).
Note that both operations are well defined since a, x1 , . . . , xn , y1 , . . . , yn all take values in the same
field, F .
Exercise 2.7 Show that (F n , F, , ) is a linear space. What is the identity element ? What is
the inverse element x of x F n ?
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
16
The most important instance of this type of linear space in these notes will be (Rn , R, +, ) with the
usual addition and scalar multiplication for vectors. The state, input, and output spaces of linear
systems will be linear spaces of this type.
The second class of linear spaces that will play a key role in linear system theory are function spaces.
Example (Function spaces) Let (V, F, V , V ) be a linear space and D be any set. Let F (D, V )
denote the set of functions of the form f : D V . Consider f, g F (D, V ) and a F and define
: F (D, V ) F (D, V ) F (D, V ) by
(f g) : D V such that (f g)(d) = f (d) V g(d) d D
and : F F (D, V ) F (D, V ) by
(a f ) : D V such that (a f )(d) = a V f (d) d D
Note that both operations are well defined since a F , f (d), g(d) V and (V, F, V , V ) is a linear
space.
Exercise 2.8 Show that (F (D, V ), F, , ) is a linear space. What is the identity element? What
is the inverse element?
The most important instance of this type of linear space in these notes will be (F ([t0 , t1 ], Rn ), R, +, )
for real numbers t0 < t1 . The trajectories of the state, input, and output of the dynamical systems
we consider will take values in linear spaces of this type. The state, input and output trajectories
will differ in terms of their smoothness as functions of time. We will use the following notation to
distinguish the level of smoothness of the function in question:
Exercise 2.9 Show that all of these sets are linear spaces. You only need to check that they are
closed under addition and scalar multiplication. E.g. if f and g are differentiable, then so is f g.
To simplify the notation, unless there is special reason to distinguish the operations and identity
element of a linear space from those of the field, from now on we will use the regular symbols +
and instead of and for the linear space operations of vector addition and scalar multiplication
respectively; in fact as for real numbers we will mostly ommit and simply write av instead of a v
for a F , v V . Likewise, unless explicitly needed we will also use 0 instead of to denote the
identity element of addition. Finally we will stop writing the operations explicitly when we define
the vector space and write (V, F ) or simply V whenever the field is clear from the context.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
17
Note that by definition a linear space and all its subspaces are linear spaces over the same field. The
equation provides a way of testing whether a given set W is a subspace or not: One needs to check
whether linear combinations of elements of W with coefficients in F are also elements of W .
Exercise 2.11 Show that if W is a subspace then for all n N, and ai F , wi W for i = 1, . . . , n
n
X
ai wi W.
i=1
1. R3
2. 2D planes through the origin.
3. 1D lines through the origin.
4. {0}.
For examples of subspace of function spaces consider (R[t], R) (polynomials of t R with real
coefficients). This is a linear subspace of C (R, R), which in turn is a linear subspace of C(R, R).
The set
{f : R R | t R, |f (t)| 1}
on the other hand is not a subspace of F (R, R).
It is easy to see that the family of subspaces of a given a linear space is closed under finite addition
and intersection.
Exercise 2.13 Let {(Wi , F )}ni=1 be a finite family of subspaces of a linear space (V, F ). Show that
(ni=1 Wi , F ) is also a subspace. Is (ni=1 Wi , F ) a subspace?
Exercise 2.14 Let (W1 , F ) and (W2 , F ) be subspaces of (V, F ) and define
W1 + W2 = {w1 + w2 | w1 W1 , w2 W2 }.
A subset of a linear space will of course not be a subspace in general. Each subset of a linear space
does, however, generate a subspace in a natural way.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
18
Definition 2.7 Let (V, F ) be a linear space and S V . The linear subspace of (V, F ) generated by S
is the smallest subspace of (V, F ) containing S.
Exercise 2.15 What is the subspace generated by {(x1 , x2 ) R2 | (x1 = 0) (x2 = 0)}? What is
the subspace generated by {(x1 , x2 ) R2 | x2 = x1 + 1}? What is the subspace of R2 generated by
{(1, 2)}?
Fact 2.5 Let (V, F ) be a linear space and S V . The subspace generated by S coincides with
Span(S).
Proof: The fact that Span(S) is a subspace and contains S is easy to check. To show that is is
the smallest subspace containing S, consider another subspace, W , that contains S and an arbitrary
v Span(S); we will show that v W and hence Span(S) W . Since v Span(S) it can be
written as
Xn
v= ai vi
i=1
The elements of Span(S) are known as linear combinations of elements of S. Notice that in general
the set S may contain an infinite number of elements; this was for example the case for the set
{(x1 , x2 ) R2 | (x1 = 0) (x2 = 0)} in Exercise 2.14. The span of S, however, is defined as the set
of all finite linear combinations of elements of S.
Definition 2.9 Let (V, F ) be a linear space. A set S V is called linearly independent if and only
if for all n N, vi S for i = 1, . . . , n
n
X
ai vi = 0 ai = 0, i = 1, . . . , n.
i=1
Note again that the set S may be infinite, but we only consider finite linear combinations to define
linear independence.
Exercise 2.16 Show that a set S V is linearly independent if and only if none of its elements
can be written as a finite linear combination of other elements in S. Show that every set S that
contains the identity element of vector addition is linearly dependent.
Example (Linearly independent set in Rn ) The following finite family of vectors are linearly
independent in (R3 , R):
{(1, 0, 0), (0, 1, 0), (0, 0, 1)}.
It is easy to show that any linear combination of these vectors is equal to (0, 0, 0) if and only if all
three coefficients are equal to 0.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
19
Example (Linearly independent functions) Consider C([1, 1], R) and let fk (t) = tk for all
t [1, 1] and k N. Clearly fk (t) C([1, 1], R).
Here C([1, 1], R) denotes the addition identity element for C([1, 1], R), i.e. the function
: [1, 1] R such that (t) = 0 for all t [1, 1].
One direction () is obvious. For the other () note that
n
X
f (t) = an fn (t) = a0 + a1 t + . . . + an tn = f (t) = 0 t [1, 1]
i=0
d dn
f (t) = f (t) = . . . = n f (t) = 0 t [1, 1]
dt dt
d dn
f (0) = f (0) = . . . = n f (0) = 0
dt dt
a0 = a1 = . . . = n!an = 0.
At this stage one may be tempted to ask how many linearly independent vectors can be found in
a given vector space. This number is a fundamental property of the vector space, and is closely
related to the notion of a basis.
Definition 2.10 Let (V, F ) be a linear space. A set of vectors S V is a basis of (V, F ) if and
only if they are linearly independent and Span(S) = V .
Fact 2.7 Let (V, F ) be a linear space. If a basis of (V, F ) with a finite number of elements exists,
then every other basis of (V, F ) has the same number of elements.
Proof: Consider two bases of (V, F ), S and S and assume that S has a finite number of elements
{vi }ni=1 . Assume, for the sake of contradiction that S does not have the same number of elements.
Without loss of generality assume further that S has more elements than S; if not, interchange S
and S since in this case S must also have a finite number of elements. Take n + 1 elements from S ,
{vi }n+1 n
i=1 and recall that since S is a basis they must be linearly independent. Since {vi }i=1 generate
V we can write v1 as a linear combination
v1 = a1 v1 + . . . + an vn .
Note that at least one of the ai must be non-zero; otherwise v1 = 0 and the set {vi }n+1
i=1 cannot be
linearly independent (Exercise 2.15). Assume, without loss of generality that a1 6= 0 and write
1 a2 an
v1 = v1 v2 . . . vn
a1 a1 a1
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
20
(where we make use of the identities in Problem 2.3). Since {vi }ni=1 is a basis any element v V
can be expressed as a linear combination of v1 , v2 , . . . , vn and hence, by the above equation, as a
linear combination of v1 , v2 , . . . , vn . In particular, we can write
v2 = b1 v1 + b2 v2 + . . . + bn vn .
Note again that the b2 , . . . , bn cannot all be zero; otherwise v2 = b1 v1 and {vi }n+1
i=1 cannot be linearly
independent. Assume, without loss of generality, that b2 6= 0 and write
1 b1 bn
v2 = v v . . . vn .
b2 2 b2 1 b2
This shows that every vector is V can be written as a linear combination of v1 , v2 , v3 , . . . , vn .
Repeat for v3 , etc. until finally
vn+1 = c1 v1 + c2 v2 + . . . + cn vn .
This, however, contradicts the assumption that the set {vi }ni=1 is linearly independent.
The fact provides conditions under which the number of elements of a basis of a linear space in well
defined, and independent of the choice of the basis. A similar statement can be made for linear
spaces which do not have a basis with a finite number of elements. In this case there are again
families of vectors whose span covers almost all elements of the vector space. The proof, however,
is considerably more involved. It also relies on concepts covered in Chapters 3 and 7 that go beyond
the purely algebraic structure considered in this chapter.
Definition 2.11 Let (V, F ) be a linear space. If a basis of (V, F ) with a finite number of elements
exists, the number of elements of this basis is called the dimension of (V, F ) and (V, F ) is called
finite dimensional. If not, (V, F ) is called infinite dimensional.
Exercise 2.17 If (V, F ) has dimension n then any set of n+ 1 or more vectors is linearly dependent.
Example (Bases) The vectors {(1, 0, 0), (0, 1, 0), (0, 0, 1)} form a basis for the linear space R3 over
the field R. This is called the canonical basis of R3 and is usually denoted by {e1 , e2 , e3 }. However
other choices of basis are possible, for example {(1, 1, 0), (0, 1, 0), (0, 1, 1)}. In fact any three linearly
independent vectors will form a basis for R3 . This shows that R3 is finite dimensional and of
dimension 3. In the same way, the vectors {(1, 0, 0), (0, 1, 0), (0, 0, 1)} form the canonical basis for
the linear space C3 over the field C; therefore C3 over C has dimension 3. On the other hand, C3
is a linear space also over the field R; in this case it has dimension 6, the following being a basis:
{(1, 0, 0), (i, 0, 0), (0, 1, 0), (0, i, 0), (0, 0, 1), (0, 0, i)}.
Exercise 2.18 Find a basis for Rn over the field R. Hence, show that Rn has dimension n over R.
Show that Cn has dimension n over C and 2n over R.
The linear space F ([1, 1], R) is infinite dimensional. We have already shown that the collection
{tk | k N} F ([1, 1], R) is linearly independent. The collection contains an infinite number of
elements and may or may not span F ([1, 1], R). Therefore any basis of F ([1, 1], R) (which must
by definition span the set) must contain at least as many elements.
Fact 2.8 The representation of a given x V with respect to a basis {b1 , . . . , bn } is unique.
Example (Representations) Let x = (x1 , x2 , x3 ) (R3 , R). The representation of x with respect
to the canonical basis is simply = (x1 , x2 , x3 ). The representation with respect to the basis
{(1, 1, 0), (0, 1, 0), (0, 1, 1)}, however, is = (x1 , x2 x1 x3 , x3 ) since
1 0 0 1 0 0
x = x1 0 + x2 1 + x3 0 = x1 1 + (x2 x1 x3 ) 1 + x3 1 .
0 0 1 0 0 1
Representations can also be defined for infinite dimensional spaces, but we will not get into the
details here. As an example, consider f (t) C ([1, 1], R). One can consider a representation of
f (t) with respect to the basis {tk | k N} defined through the Taylor series expansion. For example,
expansion about t = 0 gives
df 1 d2 f
f (t) = f (0) + (0)t + (0)t2 + . . . .
dt 2 dt2
2
which suggests that the representation of f (t) is = (f (0), df 1d f
dt (0), 2 dt2 (0), . . .). Making this state-
ment formal, however, is beyond the scope of these notes.
It turns out that all representations of an element of a linear space are related to one another:
Knowing one we can compute all others. To do this we need the concept of linear maps.
Note that both linear spaces have to be defined over the same field. For clarity we will sometimes
write A : (U, F ) (V, F ) if we need to specify the field over which the linear spaces are defined.
Example (Linear maps) Let (U, F ) = (Rn , R), (V, F ) = (Rm , R) and consider a matrix A Rmn .
Define
A : Rn Rm
u 7 A u.
It is easy to show that A is a linear map. Indeed:
A(a1 u1 + a2 u2 ) = A (a1 u1 + a2 u2 ) = a1 A u1 + a2 A u2 = a1 A(u1 ) + a2 A(u2 ).
Exercise 2.19 Show that the functions A and A are both linear.
It is easy to see that linear maps map the zero element of their domain to the zero element of their
co-domain.
Other elements of the domain may also be mapped to the zero element of the co-domain, however.
Null(A) = {u U | A(u) = V } U
Range(A) = {v V | u U : v = A(u)} V.
The word space in null space and range space is of course not accidental.
Fact 2.9 Show that Null(A) is a linear subspace of (U, F ) and Range(A) is a linear subspace of
(V, F ).
The proof is left as an exercise (Problem 2.5). It is easy to see that the properties of the null and
range spaces are closely related to the injectivity and surjectivity of the corresponding linear map,
and hence its invertibility (Problem ??).
3. If b Range(A) there exists a unique u U such that A(u) = b if and only if Null(A) =
{U }. In other words, A is injective if and only if Null(A) = {U }.
Finally, we generalise the concept of an eigenvalue to more general linear maps of a (potentially
infinite dimensional) linear space.
Definition 2.14 Let (V, F ) be a linear space and consider a linear map A : V V . An element
F is called an eigenvalue of A if and only if there exists v V such that v 6= V and A(v) = v.
In this case, v is called an eigenvector of A for the eigenvalue .
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
23
Example (Eigenvalues) For maps between finite dimensional spaces defined by matrices, the in-
terpretation of eigenvalues is the familiar one from linear algebra. Since eigenvalues and eigenvectors
are in general complex numbers/vectors we consider matrices as linear maps between complex finite
dimensional spaces (even if the entries of the matrix itself are real). For example, consider the linear
space (C2 , C) and the linear map A : C2 C2 defined by the matrix
0 1
A=
1 0
df
(A(f ))(t) = (t), t [t0 , t1 ].
dt
Exercise 2.21 Show that A is well defined, i.e. if f C ([t0 , t1 ], R) then A(f ) C ([t0 , t1 ], R).
Show further that A is linear.
One can see that in this case the linear map A has infinitely many eigenvalues. Indeed, any function
f : [t0 , t1 ] R of the form f (t) = et for R is an eigenvector with eigenvalue , since
d t
(A(f ))(t) = e = et = f (t), t [t0 , t1 ]
dt
which is equivalent to A(f ) = f .
Exercise 2.22 Let A : (V, F ) (V, F ) be a linear map and consider any F . Show that the set
{v V | A(v) = v} is a subspace of V .
Naturally, for linear maps between finite dimensional spaces that are generated by matrices a close
relation between the properties of the linear map and those of the matrix exist.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
24
The following facts can be used to link the nullity and rank of matrices. They will prove useful when
manipulating matrices later on.
1. Rank(A) + Nullity(A) = m.
2. 0 Rank(A) min{m, n}.
3. Rank(A) + Rank(B) m Rank(AB) min{Rank(A), Rank(B)}.
4. If P F mm and Q F nn are invertible then
Theorem 2.3 Let F be a field, A F nn be a matrix, and A : F n F n the linear map defined
by A(x) = Ax for all x F n . The following statements are equivalent:
1. A is invertible.
2. A is bijective.
3. A is injective.
4. A is surjective.
5. Rank(A) = n.
6. Nullity(A) = 0.
7. The columns aj = (a1j , . . . , anj ) F n form a linearly independent set {aj }nj=1 .
8. The rows ai = (ai1 , . . . , ain ) F n form a linearly independent set {ai }ni=1 .
Theorem 2.4 Let F be a field, A F nn be a matrix, and A : F n F n be the linear map defined
by A(x) = Ax for all x F n . The following statements are equivalent:
1. C is an eigenvalue of A.
2. Det(I A) = 0.
3. There exists v Cn such that v 6= 0 and Av = v. Such a v is called a right eigenvector of A
for the eigenvalue .
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
25
The proof is left as an exercise (Problem ??). The theorem shows (among other things) that the
eigenvalues of a linear map defined by a square matrix A F nn can be computed by finding the
roots of the characteristic polynomial of the matrix, defined by
A () = Det[I A] = sn + 1 sn1 + . . . + n1 s + n .
Notice that by definition the characteristic polynomial is a monic polynomial (the leading coefficient
is equal to 1) of degree n.
Definition 2.16 The spectrum of a matrix2 A F nn is the list of eigen values of A, denoted by
Spec[A] = {1 , . . . , n }.
Finally, combining the notion of eigenvalue with Theorem 2.3 leads to the following condition for
the invertibility of a matrix.
Theorem 2.5 A matrix A F nn is invertible if and only if none of its eigenvalues are equal to
the zero element 0 F .
yj = A(uj ) V for j = 1, . . . , n.
The vectors yj V can all be represented with respect to the basis {vi }m
i=1 of (V, F ). In other words
for all j = 1, . . . , n there exist unique aij F such that
m
X
yj = A(uj ) = aij vi .
i=1
Since representations are unique (Fact 2.8), the linear map A and the bases {uj }nj=1 and {vi }m
i=1
uniquely define the matrix A F mn .
2 Sometimes a distinction is made between the spectrum of the matrix (list of eigenvalues with repeated eigenvalues
included) and the spectra list of a matrix (list of eigenvalues without repetitions). We will however not make this
distinction here.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
26
Consider now an arbitrary x U . Again there exists unique representation F n of x with respect
to the basis {uj }nj=1 ,
Xn
x= j uj
j=1
Let us now see what happens to this representation if we apply the linear map A to x. Clearly
A(x) V , therefore there exists a unique representation F m of A(x) with respect to the basis
{vi }m
i=1 such that
Xm
A(x) = i vi .
i=1
What are the corresponding i F ? Recall that
Xn Xn n
X m
X m
X Xn
A(x) = A j uj = j A(uj ) = j aij vi = aij j vi .
j=1 j=1 j=1 i=1 i=1 j=1
By uniqueness of representation
n
X
i = aij j = A ,
j=1
where denotes the standard matrix multiplication. Therefore, when one looks at the representations
of vectors with respect to given bases, application of the linear map A to x U is equivalent to
multiplication of its representation (an element of F n ) with the matrix A F mn . To illustrate
this fact we will write things like
A
(U, F ) (V, F )
x 7 A(x)
mn
AF
{uj }nj=1 {vi }m
i=1
Fn 7 A F m
Theorem 2.6 The following relations between linear map operations and the corresponding matrix
representations hold:
1. Consider linear maps B : (U, F ) (V, F ) and A : (V, F ) (W, F ) where U , V and W
are finite dimensional linear spaces of dimensions n, m and p respectively. Then composition
C = A B : (U, F ) (W, F ) is also a linear map. Moreover, if we fix bases {uk }nk=1 , {vi }m
i=1
and {wj }pj=1 for the three spaces and
B A
(U, F ) (V, F ) (V, F ) (W, F )
BF mn and AF pm
{uk }nk=1 {vi }m
i=1 {vi }m
i=1 {wj }pj=1
then
C=AB
(U, F ) (W, F )
C=ABF pn
{uk }nk=1 {wj }pj=1
where denotes the standard matrix multiplication.
2. Consider an invertible linear map A : (V, F ) (V, F ) on an n-dimensional linear space V
and let A1 : (V, F ) (V, F ) denote its inverse. If A is the representation of A with respect
to a given basis of V , then A1 is the representation of A1 with respect to the same basis.
The proof is left as an exercise (Problem 2.8). Analogous statements can of course be made about
the representations of linear maps obtained by adding, or scaling other linear maps.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
27
In this section we investigate the relation between the two representations A and A.
Recall first that changing basis changes the representations of all vectors in the linear spaces. It is
therefore expected that the representation of a linear map will also change.
To derive the relation between A and A, consider first the identity map
1U
(U, F ) (U, F )
x 7 1U (x) = x
IF nn
{uj }nj=1 {uj }nj=1
QF nn
{uj }nj=1 {uj }nj=1
where 0 and 1 are the addition and multiplication identity of F . The argument used to derive the
representation of a linear map as a matrix in Section 2.7 suggests that the elements of Q F nn
are simply the representations of 1U (uj ) = uj (i.e. the elements of the basis {uj }nj=1 ) with respect
to the basis {uj }nj=1 . Likewise
1
V
(V, F ) (V, F )
x 7 1V (x) = x
IF mm
{vi }m
i=1 {vi }m
i=1
mm
P F
{vi }m
i=1 {vi }m
i=1
Exercise 2.23 Show that the matrices Q F nn and P F mm are invertible. (Recall that 1U
and 1V are bijective functions.)
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
28
and
1 2
R2 R
R2
x 7 x
1 1 R22
Q1 =
0 1
{(1, 0), (0, 1)} {(1, 0), (1, 1)}
(x1 , x2 ) 7 (x1 , x2 ) = (x1 x2 , x2 )
A = 1V A 1U : (U, F ) (V, F )
A = P A Q F mn
where A is the representation of A with respect to the bases {uj }nj=1 and {vi }m
i=1 and denotes
ordinary matrix multiplication. Since P F mm and Q F nn are invertible it is also true that
A = P 1 A Q1 .
Matrices that are related to each other in this way will play a central role in subsequent calculations.
We therefore give them a special name.
Definition 2.17 Two matrices A F mn and A F mn are equivalent if and only if there exist
Q F nn , P F mm both invertible such that A = P A Q.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
29
Theorem 2.7 Two matrices are equivalent if and only if they are representations of the same linear
map.
Proof: The if part follows from the discussion above. For the only if part use the matrices to
define linear maps from F n to F m .
1. Consider an arbitrary set S. Show that the set of bijective functions f : S S forms a group
under the operation of function composition. What does this group correspond to in the case
where S is finite?
2. Draw an equilateral triangle in R2 , say with vertices at (1, 0), ( 12 , 23 ) and ( 12 , 2 3 ). Con-
sider the set of symmetries of this triangle, i.e. the set of all functions f : R2 R2 that
map the triangle onto itself. Show that this set forms a group under the operation of function
composition. How many elements does this group contain? What is the identity element?
Repeat for the case of a square and a circle.
2. Let (R, +, ) be a ring. Show that there exists a unique element 0 R such that for all R,
+ 0 = 0 + = . Moreover, show that for all R there exists a unique () R such
that + () = () + = 0.
6= 0, = = .
Problem 2.3 (Identity and inverse properties) Let (V, F, , ) be a linear space and 0 be the
addition identity element of F . Show that for all x V , 0 x = . Moreover, show that for all
a F , x V , (a) x = (a x) = a (x).
x y = (1 + 1 , . . . , n + n )
x = ( 1 , . . . , n ).
Show that (F n , F, , ) is a linear space. What is the addition identity? What is the addition
inverse?
2. Let now (V, F, +, ) be a linear space and D an arbitrary set. Let F (D, V ) denote the set of all
functions f : D V . For f, g F (D, V ) and F define f g : D V and f : D V
by
for all d D. Show that (F (D, V ), F, , ) is a linear space. What is the addition identity?
What is the addition inverse.
Problem 2.5 (Subspaces) 1. Let U and V be linear spaces and let L(U, V ) denote the set of
linear functions A : U V . Show that L(U, V ) is a linear subspace of F (U, V ), the space
of all functions mapping U into V with the usual operations of function addition and scalar
multiplication.
2. Let U and V be linear spaces and A : U V be a linear function. Show that Range(A) is a
subspace of V and Null(A) is a subspace of U .
3. Let {Wi }ni=1 be a finite family of subspaces of V . Show that the the intersection and the direct
sum of these subspaces
n
\
Wi = {v V | i = 1, . . . , n, v Wi }
i=1
Mn
Wi = {v V | wi Wi , i = 1, . . . n, v = w1 + . . . wn }
i=1
Problem 2.6 (Basis and vector representation) Let V be a finite dimensional linear space.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
31
1. Let W be a subspace of V . Show that W is also finite dimensional and its dimension can be
no greater than that of V .
2. Show that the representation of a given x V with respect to a basis {b1 , . . . , bn } is unique.
Problem 2.7 (Rank and nullity) Let (F, +, ) be a field and consider the linear maps A : (F n , F )
(F m , F ) and B : (F p , F ) (F n , F ) represented by matrices A F mn and B F pn respectively.
Show that:
1. A is invertible.
2. None of the eigenvalues of A is equal to zero.
3. A is bijective.
4. A in injective.
5. A is surjective.
6. Rank(A) = n.
7. Nullity(A) = 0.
8. The columns aj = (a1j , . . . , anj ) F n form a linearly independent set {aj }nj=1 .
9. The rows ai = (ai1 , . . . , ain ) F n form a linearly independent set {ai }ni=1 .
1. Consider linear maps A : (V, F ) (W, F ) and B : (U, F ) (V, F ). Assume that U, V, W
have finite dimensions m, n, p respectively, and that A and B have representations A F pn
and B F nm with respect to given bases for the three spaces. Show that the composition
C = A B : (U, F ) (W, F ) has representation C = AB with respect to the same bases.
2. Consider a linear map A : (U, F ) (U, F ) where U has finite dimension n. Assume that A
has representation A F nn with respect to a given basis for U . Show that if A is invertible,
then A1 has representation A1 with respect to the same basis.
1. Consider a linear map A : (U, F ) (U, F ) where U has finite dimension n. Assume there
exists a vector b U such that the collection {b, A(b), A A(b), . . . , An1 (b)} forms a basis for
U . Derive the representation of A and b with respect to this basis.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
32
2. Consider a linear map A : (U, F ) (U, F ) where U has finite dimension n. Assume there exists
a basis bi , i = 1, . . . , n for U such that A(bn ) = bn and A(bi ) = bi + bi+1 , i = 1, . . . , n 1.
Derive the representation of A with respect to this basis.
3. Consider two matrices A, A Rnn related by a similarity transformation; i.e. there exists
Q Rnn invertible such that A = Q1 AQ. Show that Spec[A] = Spec[A].
1. C is an eigenvalue of A.
2. Det(I A) = 0.
3. There exists v Cn such that v 6= 0 and Av = v.
4. There exists Cn such that 6= 0 and T A = T .
Chapter 3
Introduction to Analysis
Consider a linear space (V, F ) and assume the F = R or F = C so that for a F the absolute value
(or modulus) |a| is well defined.
A linear space equipped with such a norm is called a normed linear space and is denoted by (V, F, kk).
A norm provides a notion of lenght for an element of a linear space. The norm can also be used
to define a notion of distance between elements of a linear space; one can think of kv1 v2 k as
the distance between two vectors v1 , v2 V .
Example (Normed spaces) In (Rn , R), the following are examples of norms:
n
X
kxk1 = |xi |, (1-norm)
i=1
v
u n
uX
kxk2 = t |xi |2 (Euclidean or 2-norm)
i=1
n
! p1
X
p
kxkp = |xi | for p 1, (p-norm)
i=1
kxk = max |xi | (infinity norm)
i=1,...,n
Exercise 3.1 Show that kxk1 , kxk2 , and kxk satisfy the axioms of a norm.
33
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
34
In fact the same definitions would also hold for the linear space (Cn , C). Different norms of course
lead to different notions of distance for the same linear space. One way to visualize these different
distance notions is through the so-called open balls.
Definition 3.2 If (V, F, kk) is a normed linear space, the (open) ball of radius r R+ centered at v V
is the set
B(v, r) = {v V | kv v k < r}.
It can be shown that B(v, r) is an open set in the sense of definition 3.3. The ball B(V , 1) is called
the unit ball of (V, F, k k).
Exercise 3.2 Draw the unit balls of the normed linear spaces (R2 , R, k k1 ), (R2 , R, k k2 ), and
(R2 , R, k k ).
For example, any ball B(v, r) is bounded, because B(v, r) B(V , kvk + r + 1).
To denote sequences, we will mostly use the notation {vi } i=0 , where i N plays the role of the
i
index variable. To indicate that a sequence converges to some vector v V we write vi v
or limi vi = v. Note that by the definition of the norm, the statement the sequence {vi } i=0
converges to v V is equivalent to limi kvi vk = 0 in the usual sense in R.
Using the definition one can also define open and closed sets.
Definition 3.5 Let (V, F, k k) be a normed space. A set K V is called closed if and only if it
contains all its limit points; i.e. if and only if for all sequences {vi }
i=1 K if vi v V then
v K. K is called open if and only if its complement V \ K is closed.
Some useful properties of open and closed sets are stated in ??.
Consider now functions between two normed spaces
f : (U, F, k kU ) (V, F, k kV )
Exercise 3.3 Show that the set of continuous functions f : (U, F, k kU ) (V, F, k kV ) is a linear
subspace of F (U, V ).
Fact 3.1 The norm k k, as a function between the normed linear spaces (V, F, k k) and (R, R, | |),
is continuous on V .
Continuity of functions has an intimate relationship with the convergence of sequences and hence
with open and closed sets. This relationship is highlighted in Problem ??.
Note that in the definitions of convergence and continuity we have left open the choice of the norm.
One may wonder if, by choosing another norm, a a convergent sequence may cease to converge, or a
continuous function can cease to be continuous. To investigate this question we need the notion of
equivallent norms.
Definition 3.7 Consider a linear space (V, F ) with two norms, k ka and k kb . The two norms are
equivalent if and only if
Exercise 3.4 Show that in this case it is also true that ml , mu R+ such that v V , ml kvkb
kvka mu kvkb .
Example (Equivalent norms) Consider x Rn and the kxk1 and kxk norms defined above.
Then
n
X
kxk = max |xi | |xi | = kxk1
i=1,...,n
i=1
n
X n
X
kxk1 = |xi | max |xj | = n max |xi | = nkxk .
j=1,...,n j=1,...,n
i=1 i=1
Therefore kxk kxk1 nkxk and the two norms are equivalent (take ml = 1 and mu = n).
Fact 3.2 Let kka and kkb be two equivallent norms on a linera space (V, F ) with F = R or F = C.
A sequence {vi }i=0 V converges to some v V in (V, F, k ka ) if and only if it converges to v in
(V, F, k kb ).
Proof: Suppose that the sequence {vi } i=0 converges to v V with respect to the norm k ka , that
is for all > 0 there exists N N such that kvm vka < for all m N . Due to equivalence of the
norms, there exists mu > 0 such that ml kvka kvkb mu kvka for all v. Fix an arbitrary > 0.
Then there exists N N such that, for all m N , kvm vka < mu ; But then, with the same N ,
for all m N
kvm vkb mu kvm vka < mu = .
mu
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
36
Since function continuity, open and closed sets were all defined in terms of convergence of sequences,
this fact further implies that open/closed sets defined using one norm remain open/closed for any
other equivallent norm. Likewise, continuous functions defined with respect to a pair of norms
remain continuous with respect to any other pair of respectively equivalent norms.
The fact that kxk1 and kxk are equivalent norms on Rn is not a coincidence. A remarkable result
states, indeed, that any two norms on a finite-dimensional space are equivalent. To show this, we
will use the following fact, which is indeed a corollary of two fundamental theorems in real analysis.
Fact 3.3 (A corollary to the theorems of Heine-Borel and of Weierstrass). A continuous function
f : S R defined on a closed and bounded subset S of a finite-dimensional linear space V attains
a maximum M R and a minimum m R on S. In other words, there exist xM , xm S such that
Proof: Assume that a certain norm k k has already been defined on V ; say that this is the norm
with respect to which the concepts of ball, convergence, and continuity have been defined. Let
kkA be another norm on V . The unit sphere S = {x V | kxk = 1} is bounded. Moreover, S is also
closed, as it is the inverse image of the closed set {1} under the continuous function k k : V R
(see Problem ?? and Fact ??). Therefore k kA , being continuous, admits a minimum m and a
maximum M on S. In particular M m > 0, because kxkA = 0 x = 0 kxk = 0 (it is not the
case for x S). Now, if x V and x 6= 0, x/kxk S. Hence,
x
kxkA = kxk
mkxk
kxk
A
x
kxkA = kxk
M kxk
kxk
A
The above inequalities are trivially true also when x = 0, and prove that k kA is equivalent to
k k. Since the choice of k kA was left open, it is immediate to see that any two norms, being both
equivalent to k k, are equivalent to each other.
kf k := sup |f (t)|
t[t0 ,t1 ]
is a norm on (C([t0 , t1 ], Rn ), R). Indeed, for all continuous functions f , kf k is a positive or null
quantity, and it is trivially equal to zero if and only if f (t) is the zero-function. Moreover,
Norms on function spaces can also be defined using integration. For example, for f C([t0 , t1 ], Rn )
we can define the following:
Z t1
kf k1 = kf (t)k2 dt
t
s0Z
t1
kf k2 = kf (t)k22 dt
t0
Z t1 p1
kf kp = kf (t)kp2 dt
t0
Here, the integral takes the role that the finite summation has in finite-dimensional spaces. Since
the 2-norm that appears in the integrands is a continuous function (Fact??), the integrands are also
continuous functions of the variable t, and all these quantities are well-defined. Moreover, all of them
are norms. Take for example the first one. Of course we have kf k1 0 for all f C([t0 , t1 ], Rn ); the
norm of the zero-function is of course equal to zero; on the other hand, given a continuous function
f , if f is nonzero at a point it must be nonzero in a whole subinterval of [t0 , t1 ], therefore kf k1 is
also nonzero. (Summing up, kf k1 = 0 if and only if f = 0.) Moreover,
Z t1 Z t1
kf k1 = |f (t)| dt = || |f (t)| dt = || kf k1
t0 t0
Z t1 Z t1
kf + gk1 = |f (t) + g(t)| dt (|f (t)| + |g(t)|) dt = kf k1 + kgk1
t0 t0
It is also easy to see that the use of the 2-norm in the integrands is arbitrary: Since norms in Rn
are all equivallent one could have used any other norm, leading to other norms on C([t0 , t1 ], Rn ).
Exercise 3.5 Show that the function norms obtained by replacing kf (t)k2 by another norm on Rn
in the integrands are equivalent.
Motivated by this fact and in analogy to finite-dimensional spaces, the reader may be tempted to
think that the all the norms introduced above are equivalent. Unfortunately, this is not the case.
Example (Non-equivalent norms) Consider the functions f C([0, 1], R) and the two norms
Z 1
kf k = sup |f (t)| and kf k1 = |f (t)|dt
t[0,1] 0
We will show that the two norms are not equivalent. Assume, for the sake of contradiction, that
there exist ml , mu R+ such that for all f C([0, 1], R) we have ml kf k1 kf k mu kf k1 . Since
Z 1 Z 1
kf k1 = |f (t)| dt kf k dt = kf k,
0 0
one can indeed choose ml = 1. On the other hand, it is impossible to find mu that satisfies the
other inequality for all f . Assume, for the sake of contradiction, that such an mu exists and consider
1
the family of functions fn (t) = tn C([0, 1], R) for n N. Clearly kfn k = 1 and kfn k1 = n+1 .
Therefore kfa k = (n + 1)kfa k1 . Selecting n + 1 > mu leads to a contradiction, showing that the
two norms are not equivalent.
3.3 Completeness
In our treatment of the solution of differential equations below we will be faced with the task of
establising the convergence of a series in a certain linear space. In principle, one should be able to
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
38
do so using Definition 3.2. In practice, however, this will not be possible as the limit is not known
beforehand. From the definition it is, however, apparent that if the sequence {vi } i=0 converges to
some v V the vi need to be getting closer and closer to v, and therefore closer and closer to each
other. Can we decide whether the sequence converges just by looking at kvi vj k?
Exercise 3.6 Show that a sequence is Cauchy if and only if for all > 0 there exists N N such
that kvm vn k < for all m, n N .
Definition 3.9 The normed space (V, F, k k) is called complete (or Banach) if and only if every
Cauchy sequence converges.
This observation is in fact is the reason why the real number system was introduced in the first
place, and is so fundamental that in many expositions of analysis it is taken as an axiom. This fact
is also the main tool used to establishing completeness of many other spaces.
Theorem 3.1 Let F = R or F = C and assume that (V, F ) is finite-dimensional. Then (V, F, k k)
is a Banach space for any norm k k.
Proof: Suppose that F = R. Let {b1 , , bd } be a basis of V , so that every vector v V can be
Pd
represented in an unique way as v = k=1 xk bk for some xk R. On Rd let the function k k be
defined by
d
X
kxk = k(x1 , , xd )k :=
xk bk
k=1
But then, we have the same property for each of the k-th components:
> 0 N N m, n N, |xnk xm
k | < .
In other words, the k-th components {xnk } n=0 form a Cauchy sequence in R. Since R is complete,
this sequence converges. Define x = (x1 , , xk , , xd ), where
xk := lim xnk
n
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
39
(we can take the limit within k k because the norm is continuous). Therefore, {xn } n=0 is conver-
gent, and the same happens for {vn } n=0 ; since the latter is arbitrary, (V, F ) is Banach.
To establish the same result for complex spaces (F = C) one can repeat the same proof after showing
that C itself is complete. This should now be easy, and is left as an exercise.
Corollary 3.1 Any finite-dimensional subspace W of a linear space (V, F, k k) is a closed subset
of V .
This fact will come in handy in Chapter 7 in the context of inner product spaces. The state, input
and output spaces of our linear systems will be of the form Rn , and will all be Banach spaces. One
may be tempted to think that this is more generally true. Unfortunately infinite-dimensional spaces
are less well-behaved.
Example (Non-Banach space) Consider the normed space (C([1, 1], R), R, k k1 ) of continuous
functions f : [1, 1] R with the 1norm
Z 1
kf k1 = |f (t)|dt.
1
It is easy to see that the sequence is Cauchy. Indeed, if we take N, m {1, 2, . . .} with m N ,
Z 1 Z 1/m Z 1/N
kfm fN k = |fm (t) fN (t)|dt = (mt N t)dt + (1 N t)dt
1 0 1/m
1/m 1/N
t2 t2 mN 1
= (m N ) + 1N = .
2 0 2 1/m 2mN 2N
Therefore, given any > 0 by selecting N > 1/(2) we can ensure that for all m N , kfm fN k < .
One can guess that the sequence {fi } i=1 converges to the function
0 if t < 0
f (t) =
1 if t 0.
Indeed Z 1/i
1
kfi f k1 = (1 it)dt = 0.
0 2i
The problem is that the limit function f F ([1, 1], R) is discontinuous and therefore f 6 C([1, 1], R).
Hence C([1, 1], R) is not a Banach space since the Cauchy sequence {fi }ni=1 does not converge to
an element of C([1, 1], R).
Mercifully, there are several infinite dimensional spaces that are Banach.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
40
kfm fN k1
fm (t)
1
kfm fN k
fN (t)
1 1
1 m N 1 t
Proof: We articulate the proof in three steps: First, given a sequence of continuous functions which
is Cauchy in the norm k k , we define a pointwise limit function; Second, we prove that the
sequence converges to this function; Third, we prove that this function is continuous.
Let therefore {fn } n=1 be a Cauchy sequence of functions defined and continuous on the closed
interval [t0 , t1 ]. For each t [t0 , t1 ], {fn (t)} n
n=1 is a Cauchy sequence of vectors in R (why?). Since
n
R is a Banach space (being finite-dimensional!), every such sequence has a limit. We define the
function f as follows:
f (t) := lim fn (t)
n
Now it happens that {fn }converges to f also with respect to the norm k k . Indeed, the fact
n=1
that {fn }
n=1 is Cauchy means that, for every > 0, there exists N N such that
Now fix an > 0. Since {fn }n=1 converges to f with respect to k k , there exists n such that
kf fn k < /3, and therefore kf (t) fn (t)k < /3 for all t and in particular kfn (t) f (t)k < /3.
On the other hand, since each function of the sequence is continuous, and in particular so is fn ,
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
41
there exists > 0 such that kfn (t) fn (t)k < /3 whenever |t t| < . Thus, for all > 0 there
exists (n and) > 0 such that, if |t t| < , then
kf (t) f (t)k /3 + /3 + /3 = ,
and f is continuous.
Summarizing, given an arbitrary Cauchy sequence {fn } n
n=1 in C([t0 , t1 ], R, R ), we can construct a
n
function f C([t0 , t1 ], R, R ) such that kfn f k 0.
Since every continuous function defined on the closed interval [t0 , t1 ] is bounded, its square-modulus
is also continuous, hence integrable (in the Riemann sense, that is, according to the definition of
integral given in calculus courses), so that k k2 is indeed a well-defined norm.
If we consider the space C(R+ , Rn ) of continuous functions over [0, +[ the discourse is different:
Even a continuous and bounded function does not need to be square-integrable (think of a nonzero
R +
constant function: 0 52 dt = +), in which case the norm is not well-defined. Nevertheless, the
space can be restricted. In principle, we could define
Z +
n 2
V = f C(R+ , R ) |f (t)| dt < + (3.3)
0
and it turns out that V would again be a well-behaved normed linear space.
Exercise 3.9 Show that the sum of square-integrable functions is also square-integrable (hence, V
is a subspace of C(R+ , Rn )).
Now, in the perspective of assuming C([t0 , t1 ], Rn ) and C(R+ , Rn ) as the environments for our
further developments, we could wishfully think of proving that C([t0 , t1 ], Rn ) and V are Banach.
Unfortunately, neither of them is complete (The previous example Non-Banach space shows this
fact for the space C([1, 1], R)).
If the only possible pathology was that of a single discontinuity of the limit function, or a finite
number of discontinuities, or even more, we could relax the continuity hypothesis and go for spaces
of piecewise-continuous functions or something like that. New issues would arise. Namely, if we allow
discontinuous functions to enter the picture, k k2 is not a well-defined norm anymore; Consider for
R +
example the function f that takes the value 1 at t = 0 and 0 everywhere else. Then 0 |f (t)|2 dt =
0, but f is not the zero function. This issue is easily solved by identifying those functions that are
equal except for a small set of points, such as finitely many points, or countably many points, or
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
42
anyway a set of points such that the integral of a function does not change if we alter the value of
the function on that set; f should therefore be identified with the zero function, because it differs
from it on only one point. Such functions are said to be equal almost everywhere.
But the situation is even worse. While the functions of interest are assumed to possess some sort of
continuity, there are always Cauchy sequences with respect to the norm k k2 that do not converge
to a function having that sort of continuity. Loosely speaking, what can happen is even that the
limit function can fail to be integrable in the sense of Riemann.
The canonical way in which this problem is solved is by extending the definition of integral. The
Lebesgue integral of a function is a technical construction that generalizes in many ways the notion of
integral that the reader remembers from calculus courses (the Riemann integral). When a function
f is integrable in the sense of Riemann (either over [t0 , t1 ] or over R+ ), the Lebesgue integral of
f always exists and coincides with its Riemann integral. But Lebesgue integration allows for a
much larger class of functions to be integrable, comprising functions way, way wilder than just
discontinuous at finitely many points. For example, the function d : [0, 1] R such that d(t) = 1
if t Q, and d(t) = 0 elsewhere is discontinuous at all points and is not integrable according to
R1
Riemann; but we mention, albeit without any justification, that the Lebesgue integral 0 d(t)dt
exists and is equal to 0. Lebesgue integration also provides more general convergence theorems, and
an easier way to reduce multiple integrals to iterated integrals of one variable (Fubinis theorem).
If the reader proceeds in the field of System Theory, he or she will have to seriously study Lebesgue
integration at some point. But since the long excursus into measure theory required by this study
leads far beyond the scope of this course, we will avoid this discussion altogether, leave to the reader
the bare statement that the Lebesgue integral is a more general notion of integral, and instead
rely, without providing a proof, on the following fundamental result of mathematical analysis:
Fact 3.5 Given a normed linear space (V, F, k k), there exist a complete normed linear space
(V , F, k k) and a norm-preserving linear function : V V (that is, k(x)k = kxk for all x V )
such that every point in V is a limit point of some sequence of points in (V ). V is called the
completion of V .
In other words, every normed space V can be mapped to a dense subspace of a complete space
V ; the map preserves the distances between points (it is an isometry), so that V and the dense
subspace can be identified to all practical purposes. To gain intuition, the reader can think at the
completion of V as the union of V and the set of all the missing limit points of the non-convergent
Cauchy sequences of V . For example, R is the completion of Q. On the other hand, if V is already
complete, then the completion of V is V itself (hence, the completion R is R). What is now of
greatest interest for us is the following
Fact 3.6 The completion of (C([t0 , t1 ], Rn ), kk2 ) is the space of Lebesgue square-integrable functions
Rt
over [t0 , t1 ], that is the set of all those functions f : [t0 , t1 ] Rn for which the integral t01 |f (t)|2 dt
exists (in the sense ofqLebesgue) and is finite; this space is a complete normed space with respect
R t1
to the norm kf k2 := 2
t0 |f (t)| dt, provided that functions which are equal almost everywhere
are identified. All the functions f whose square is integrable in the sense of Riemann belong to
this space, and their norm according to both the definitions of integral coincide. This completion is
denoted L2 ([t0 , t1 ], Rn ).
An analogous fact can be stated for the space (3.2), namely the completion of V is the Banach space
R +
L2 (R+ , Rn ) of those functions f : R+ Rn for which the Lebesgue integral 0 |f (t)|2 dt exists
qR
+
and is finite, equipped with the norm 0 |f (t)|2 dt.
L2 ([t0 , t1 ], Rn ) and L2 (R+ , Rn ) are now true Banach spaces. Let us restate it clearly: Every sequence
{fn }n=1 of Lebesgue square-integrable functions such that limm,n kfm fn k2 = 0 converges to
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
43
Exercise 3.10 Show that (Rmn , R) is a linear space with the usual operations of matrix addition
and scalar multiplication. (Hint: One way to do this is to identify Rmn with Rnm .)
More commonly used are the norms derived when one considers matrices as linear maps between
linear spaces. We start with a more general definition.
Definition 3.10 Consider the linear space of continuous functions f : U V between two normed
spaces (U, F, k kU ) and (V, F, k kV ). The induced norm of f is defined as
kf (u)kV
kf k = sup .
u6=U kukU
Fact 3.7 Consider two normed spaces (U, F, kkU ) and (V, F, kkV ) and a continuous linear function
A : U V . Then
kAk = sup kA(u)kV .
kukU =1
Notice that the induced norm depends not only on the function but also on the norms imposed on
the two spaces.
Example (Induced matrix norms) Consider A F mn and consider the linear map A : F n
F m defined by A(x) = A x for all x F n . By considering different norms on F m and F n different
induced norms for the linear map (and hence the matrix) can be defined:
kAxkp
kAkp = sup
xF n kxkp
Xm
kAk1 = max |aij | (maximum column sum)
j=1,...,n
i=1
1
kAk2 = max j (AT A) 2 (maximum singular value)
j=1,...,n
n
X
kAk = max |aij | (maximum row sum)
i=1,...,m
j=1
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
44
It turns out that the induced norm of linear maps is intimately related to their continuity.
Theorem 3.3 Consider two normed spaces (U, F, k kU ) and (V, F, k kV ) and a linear function
A : U V . The following statements are equivalent
1. A is continuous.
2. A is continuous at U .
3. The induced norm kAk is finite.
An easy corrolary is that linear functions between finite dimensional spaces are always continuous.
Corollary 3.2 All linear functions A : U V between two finite dimensional normed spaces
(U, F, k kU ) and (V, F, k kV ) are continuous.
Proof: Fix bases for the two spaces so that the linear function A is represented by multiplication
by a matrix A. Since both spaces are finite dimensional all norms are equivallent; for simplicity let
k kU and k kV both be the corresponding infinity norms.
Note that if A = 0 the function A is constant and hence continuous. Otherwise, consider arbitrary
an u U and note that kA(u) A(u )k kAk ku u k . Clearly kAk (the maximum row
sum) is a finite number. Therefore, for all > 0 there exists = /kAk such that ku u k <
implies that kA(u) A(u )k < . Hence the function is continuous.
Theorem 3.4 Consider linear functions A, A : (V, F, kkV ) (W, F, kkW ) and B : (U, F, kkU )
(V, F, k kV ) between normed linear spaces. Let k k denote the induced norms of the linear maps.
where
The difficult part conceptually is equation (3.3), a linear ordinary differential equation (ODE) with
time varying coefficients A(t) and B(t). Equation (3.3) is a special case of the (generally non-linear)
ODE
x(t) = f (x(t), u(t), t) (3.6)
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
45
with
t R, x(t) Rn , u(t) Rm
f : Rn Rm R Rn .
The only difference is that for linear systems the function f (x, u, t) is a linear function of x and u
for all t R.
In this section we are interested to find solutions (also known as trajectories, flows, . . . ) of
the ODE (3.5). In other words:
Given:
f : Rn Rm R Rn dynamics
n
(t0 , x0 ) R R initial condition
m
u() : R R input trajectory
Find:
x() : R Rn state trajectory
Such that:
x(t0 ) = x0
d
x(t) = f (x(t), u(t), t) t R
dt
While this definition is acceptable mathematically it tends to be too restrictive in practice. The
problem is that according to this definition x(t) should not only be continuous as a function of time,
but also differentiable; since the definition makes use of the derivative dx(t)/dt, it implicitly assumes
that the derivative is well defined. This will in general disallow input trajectories u(t) which are
discontinuous as a function of time. We could in principle only allow continuous inputs (in which
case the above definition of a solution should be OK) but unfortunately many interesting input
functions turn out to be discontinuous.
Example (Hasty driver) Consider a car moving on a road. Let y R denote the position of the
car with respect to a fixed point on the road, v R denote the velocity of the car and a R denote
its acceleration. We can then write a state space model for the car by defining
y
x= R2 , u = a R
v
Assume now that we would like to select the input u(t) to take the car from the initial state
0
x(0) =
0
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
46
for some Vmax > 0, amin < 0 < amax . It turns out that the optimal solution for u(t) is discontinuous.
Assuming that yF is large enough it involves three phases:
1. Accelerate as much as possible (u(t) = amax ) until the maximum speed is reached (x2 (t) =
Vmax ).
2. Coast (u(t) = 0) until just before reaching yF .
3. Decelerate as much as possible (u(t) = amin ) and stop exactly at yF .
Unfortunately this optimal solution is not allowed by our definition of solution given above.
To make the definition of solution more relevant in practice we would like to allow discontinuous
u(t), albeit ones that are not too wild. Measurable functions provide the best notion of how
wild input functions are allowed to be and still give rise to reasonable solutions for differential
equations. Unfortunately the proper definition of measurable functions requires some exposure to
measure theory and is beyond the scope of these notes; the interested reader is referred to [?] for a
discussion. For our purposes the following, somewhat simpler definition will suffice.
1. D left and right limits of u exist, i.e. limt + u(t) and limt u(t) exist and are finite.
Moreover, u( ) = limt + u(t).
2. t0 , t1 R with t0 < t1 , D [t0 , t1 ] contains a finite number of points.
We will use PC([t0 , t1 ], Rm ) to denote the linear space of piecewise continuous functions f : [t0 , t1 ]
Rm (and similarly for PC(R, Rm )).
Exercise 3.11 Show that (PC([t0 , t1 ], Rm ), R) is a linear space under the usual operations of func-
tion addition and scalar multiplication.
Definition 3.11 includes all continuous functions, the solution to the hasty driver example, square
waves, etc. Functions that are not included are things like 1/t or tan(t) (that go to infinity for some
t R), and wild things like
0
(t h1/2) (t 0)
1 1
u(t) = 1 t 2k+1 , 2k
h k = 1, 2, . . .
1 1 1
t ,
2(k+1) 2k+1
that have an infinite number of discontinuity points in the finite interval [0, 1/2].
Let us now return to our differential equation (3.5). Let us generalize the picture somewhat by
considering
x(t) = p(x(t), t) (3.7)
where we have obscured the presence of u(t) by defining p(x(t), t) = f (x(t), u(t), t). We will impose
the following assumption of p.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
47
We are now in a position to provide a formal definition of the solution of ODE (3.6).
1. (t0 ) = x0 , and
d
2. t R \ D, is differentiable at t, and dt (t) = p((t), t).
The solution only needs to be continuous and differentiable almost everywhere, i.e. everywhere
except at the set of discontinuity points, D, of p. This definition allows us to include input trajectories
that are discontinuous as functions of time, but still reasonably well behaved.
Is this a good definition? For example, is it certain that there exists a function satisfying the
conditions of the definition? And if such a function does exist, is it unique, or can there be other
functions that also satisfy the same conditions? These are clearly important questions in control
theory where the differential equations usually represent models of physical systems (airplanes, cars,
circuits etc.). In this case the physical system being modeled is bound to do something as time
progresses. If the differential equation used to model the system does not have solutions or has
many, this is an indication that the model we developed is inadequate, since it cannot tell us what
this something is going to be.
The answer to these questions very much depends on the function p. Notice that in all subsequent
examples the function p(x, t) is independent of t; therefore all potential problems have noting to do
with the fact that p is piecewise continuous as a function of time.
For the initial condition (t0 , x0 ) = (0, c) with c > 0, Definition 3.12 allows us to define the solution
(t) = c t for t c. There does not exist a function, however, that satisfies Definition 3.12 for
t > c.
In particular, for the initial condition (t0 , x0 ) = (0, 0) the solution is undefined for all t > 0.
The problem here is that the function p(x, t) is discontinuous in x. The example suggests that to
be able to use Definition 3.12 effectively we need to exclude p(x, t) that are discontinuous in x.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
48
Another alternative would be to relax the Definition 3.12 further to allow the so-called Filippov
solutions [?, ?], but this is beyond the scope of these notes.
Is the existence of a solution guaranteed if we restrict our attention to p(x, t) that are continuous in
x? Unfortunately not!
Notice that as t 1/c the solution escapes to infinity; no solution is defined for t 1/c. As in
the previous example, Definition 3.12 allows us to define a solution until a certain point in time,
but no further. There is a subtle difference between the two examples, however. In the latter case
the solution can always be extended for a little extra time; it is defined over the right-open interval
(, 1/c). In the former case, on the other hand, a dead end is reached in finite time; for initial
conditions (0, c) the solution is only defined over the closed interval (, c].
The problem here is that p(x, t) grows too fast as a function of x. To use Definition 3.12 we will
therefore need to exclude functions that grow too fast.
Even if these exclusions work and we manage to ensure that a solution according to Definition 3.12
exists, is this solution guaranteed to be unique, or can there be many functions (t) satisfying the
conditions of the definition for the same (t0 , x0 )? Unfortunately the answer is yes.
Notice that in this case the solution is not unique. In fact there are infinitely many solutions, one
for each a 0.
The problem here is that p(x, t) is too steep, since its slope goes to infinity as x tends to 0. To
use Definition 3.12 we will therefore need to also exclude functions that are too steep.
Functions that are discontinuous, grow too fast or are too steep can all be excluded at once by
the following definition.
Definition 3.13 The function p : Rn R Rn is globally Lipschitz in x if and only if there exists
a piecewise continuous function k : R R+ such that
x, x Rn , t R kp(x, t) p(x , t)k k(t)kx x k.
k(t) is called the Lipschitz constant of p at t R.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
49
Example (Lipschitz functions) One can easily verify that linear functions are Lipschitz; we will
do so when we introduce linear systems. All differentiable functions with bounded derivatives are
Lipschitz. However, not all Lipschitz functions are differentiable. For example, the absolute value
function | | : R R is Lipschitz (with Lipschitz constant 1) but not differentiable at x = 0.
All Lipschitz functions arepcontinuous, but not all continuous functions are Lipschitz. For exam-
ple, the functions x2 and |x| from R to R are both continuous, but neither is Lipschitz. x2 is
locally Lipschitz, i.e. x0 R, > 0, kx0 : R R+ piecewise continuous such that x R:
In the following section we will show that global Lipschitz continuity is sufficient to ensure the
existence and uniqueness of solutions for ordinary differential equations. Global Lipshitz continuity
is indeed a tight sufficient condition for existence and uniqueness of solutions, albeit not a necessary
one. The examples outlined earlier in this section demonstrate that there exist differential equation
defined by non-Lipschitz functions that do not posess unique solutions. On the other hand, one can
also find differential equation defined by non-Lipschitz functions that do posess unique solutions for
all intial conditions; to establish this fact, however, more work is needed on a case by case basis. The
existence and uniqueness results discussed here can be further fine tuned (assuming for example local
Lipschitz continuity to derive local existence of solutions). In the Chapter 4, however, we will show
that linear differential equations, the main topic of these notes, always satisfy the global Lipschitz
assumption. We will therefore not pursue such refinements here, instead we refer the interested
reader to [?, ?, ?].
1. (t0 ) = x0 .
d
2. t R \ D, dt (t) = p((t), t).
The proof of this theorem is rather sophisticated. We will build it up in three steps:
1. Background lemmas.
2. Proof of existence (construction of a solution).
3. Proof of uniqueness.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
50
Proof:(Sketch) Roughly speaking one can approximate the integral by a sum, use triangle inequality
on the sum and take a limit. Note that the absolute value on the right had side is necessary, since
the integral there will be negative if t1 < t0 .
1. m, k N, (m + k)! m! k!.
cm
2. c R, limm m! = 0.
Proof: Part 1 is easy to prove by induction (Theorem 1.2). For part 2, the case |c| < 1 is trivial,
since cm 0 while m! . For c > 1, take M > 2c integer. Then for m > M
cm cM cmM cM cmM cM 1
= mM
= mM 0.
m! M !(M + 1)(M + 2) . . . m M ! (2c) M! 2
d
f (t) = g(t).
dt
Notice that u(t) U (t), since the absolute value is redundant as u() and k() are non-negative
and t t0 . By the fundamental theorem of calculus U is continuous and wherever k and u are
continuous
d
U (t) = k(t)u(t).
dt
Then
Rt Rt
k( )d k( )d
u(t) U (t) u(t)k(t)e t0 U (t)k(t)e t0
R
d t k( )d t k( )d
R
U (t) e t0 U (t)k(t)e t0 0
dt
R
d t k( )d d Rtt k( )d
U (t) e t0 + U (t) e 0 0
dt dt
d t k( )d
R
U (t)e t0 0
dt R t
k( )d
U (t)e t0
decreases as t increases
Rt R t0
k( )d
U (t)e t0
U (t0 )e t0 k( )d
t t0
Rt
k( )d
U (t)e t0
U (t0 ) = c1 t t0
Rt
k( )d
u(t) U (t) c1 e t0
x(t) = p(x(t), t)
passing through (t0 , x0 ) R Rn using an iterative procedure. This second step of our existence-
uniqueness proof will itself involve three steps:
1. x0 (t) = x0 t R
Rt
2. xm+1 (t) = x0 + t0 p(xm ( ), )d, t R.
The generated sequence of functions is known at the Picard Iteration. Notice that all the functions
xm () generated in this way are continuous by construction.
Consider any t1 , t2 R such that t1 t0 t2 . Let
Notice that under the conditions of the theorem both k and T are non-negative and finite. Let k k
be the infinity norm on Rn . Then for all t [t1 , t2 ]
Z t Z t
kxm+1 (t) xm (t)k =
x
0 + p(x m ( ), )d x0 p(xm1 ( ), )d
0 t 0 t
Z t
=
[p(xm ( ), ) p(xm1 ( ), )]d
t0
Z t
kp(xm ( ), ) p(xm1 ( ), )k d (Fact 3.6)
t
Z 0t
k( ) kxm ( ) xm1 ( )k d (p is Lipschitz in x)
t0
Z t
k kxm ( ) xm1 ( )k d
t0
For m = 0 Z t Z
t2
kx1 (t) x0 (t)k kp(x0 , )k d kp(x0 , )k d = M,
t0 t1
for some non-negative, finite number M (which of course depends on t1 and t2 ). For m = 1
Z t
kx2 (t) x1 (t)k k kx1 ( ) x0 ( )k d
t0
Z t
k M d = kM |t t0 |.
t0
For m = 2,
Z t
kx3 (t) x2 (t)k k kx2 ( ) x1 ( )k d
t0
Z t 2
k (t t0 )2
k
kM | t0 |d = M .
t0 2
Step 2.2: We show that the sequence xm () is a Cauchy sequence in the Banach space (C([t1 , t2 ], Rn ), k
k ). Note that xm () is continuous by construction by the fundamental theorem of calculus. We
therefore need to show that xm () has bounded infinity norm for all m and that
Step 2.3: Finally, we show that the limit function () : [t1 , t2 ] Rn solves the differential equation.
We need to verify that:
1. (t0 ) = x0 ; and
d
2. t [t1 , t2 ] \ D, dt (t) = p((t), t).
For the first part, notice that by construction x0 (t0 ) = x0 and for all m > 0 since
Z t0
xm+1 (t0 ) = x0 + p(xm ( ), )d = x0 .
t0
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
54
Therefore Rt
xm+1 (t) = x0 + t0 p(xm ( ), )d
Rt
(t) = x0 + t0
p(( ), )d
By the fundamental theorem of calculus is continuous and
d
(t) = p((t), t) t [t1 , t2 ] \ D.
dt
Therefore, our iteration converges to a solution of the differential equation for all t [t1 , t2 ]. Since
[t1 , t2 ] are arbitrary, we can see that our iteration converges to a solution of the differential equation
for all t R, by selecting t1 small enough and t2 large enough to ensure that t [t1 , t2 ].
Can there be more solutions besides the that we constructed in Section 3.7.2? It turns out that
this is not the case. Assume, for the sake of contradiction, that there are two different solutions
(), () : R Rn . In other words
Therefore
d
(( ) ( )) = p(( ), ) p(( ), ), R \ D
d
Z t
(t) (t) = [p(( ), ) p(( ), )]d, t R
t0
Z t Z t
k(t) (t)k
kp(( ), ) p(( ), )kd k( )k( ) ( )kd
t0 t0
Z t
k(t) (t)k c1 + k( )k( ) ( )kd c1 0.
t0
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
55
Letting u(t) = k(t) (t)k and applying the Gronwall lemma leads to
Rt
k( )d
k(t) (t)k c1 e t0
, c1 0.
Letting c1 0 leads to
k(t) (t)k = 0 (t) = (t) t R
which contradicts the fact that () 6= ().
This concludes the proof of existence and uniqueness of solutions of ordinary differential equations.
From now on we can talk about the solution of the differential equation, as long as the conditions
of the theorem are satisfied. It turns out that the solution has several other nice properties, for
example it varies continuously as a function of the initial condition and parameters of the function
p.
Unfortunately, the solution cannot usually be computed explicitly as a function of time. In this case
we have to rely on simulation algorithms to approximate the solution using a computer. The nice
properties of the solution (its guarantees of existence, uniqueness and continuity) come very handy
in this case, since they allow one to design algorithms to numerically approximate the solutions and
rigorously evaluate their convergence properties. Unfortunately for more general classes of systems,
such as hybrid systems, one cannot rely on such properties and the task of simulation becomes much
more challenging.
One exception to the general rule that solutions cannot be explicitly computed is the case of linear
systems, where the extra structure afforded by the linearity of the function p allows us to study the
solution in greater detail. We pursue this direction in Chapter 4.
Problem 3.2 (Ball of a given norm) Let B(a; r), with a V and r R+ , denote a ball of a
normed vector space (V, F, || ||), and let be the zero vector. Show that:
(Hint: Show 1, 3 and 5, then show that 3 implies 2 and 5 implies 4.)
Problem 3.4 (Continuity) Let f : (U, F, k kU ) (V, F, k kV ) be a function between two normed
spaces. The following statements are equivalent:
1. f is continuous.
2. For all sequences {ui }
i=1 U
f 1 (K) = {u U | f (u) K}
is open.
4. For all K V closed, the set f 1 (K) is closed.
Problem 3.5 (Equivalent Norms) Let (V, F ) be a linear space. Let ||||a and ||||b be equivalent
norms on (V, F ). Let v V , X V and let {vi }iN be a sequence of elements of V . Show that:
1. {vi }iN is Cauchy w.r.t. || ||a {vi }iN is Cauchy w.r.t. || ||b ;
i i
2. vi v w.r.t. || ||a vi v w.r.t. || ||b ;
3. X is dense in V w.r.t. || ||a X is dense in V w.r.t. || ||b .
Problem 3.6 (Induced norms) Let (U, F, || ||U ) and (V, F, || ||V ) be normed spaces, and let U
be the zero vector of U .
is a well-defined norm for the space of all operators F : U V (you may assume that the
space of operators is a linear space over F ).
2. Let A : U V be a linear operator. For the induced norm ||A||, show that
where
t R, x(t) Rn , u(t) Rm , y(t) Rp
and
58
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
59
We assume that x(t) Rn and u(t) Rm are small (i.e. we start close to the optimal trajectory).
The plan is to use u(t) to ensure that x(t) remains small (i.e. we stay close to the optimal
trajectory). Note that
x(t) = x (t) + x(t) x(t) = x (t) + x(t)
= f (x (t), u (t)) + x(t).
Assume that
f1 (x, u)
..
f (x, u) = .
fn (x, u)
where the functions fi : Rn U R for i = 1, . . . , n are differentiable. For t [0, T ] define
f1 f1
x1 (x (t), u (t)) . . . xn (x (t), u (t))
f .. .. ..
(x (t), u (t)) = . . . = A(t) R
nn
x fn fn
x1 (x (t), u (t)) . . . xn (x (t), u (t))
f1 f1
u1 (x (t), u (t)) . . . um (x (t), u (t))
f .. .. ..
(x (t), u (t)) = . . . = B(t) R
nm
.
u fn f
u1 (x (t), u (t)) . . . um (x (t), u (t))
n
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
60
Then
A(t) : X X B(t) : U X
C(t) : X Y D(t) : U Y
p
parametrized by a real number t R. Fix bases, {ei }ni=1 for (X, R), {fi }m
i=1 for (U, R), and {gi }i=1
for (Y, R). Let A(t), B(t), C(t), and D(t) denote respectively the representation of the linear maps
A(t), B(t), C(t), and D(t) with respect to those bases,
A(t) B(t)
(X, R) (X, R) (U, R) (X, R)
nn
A(t)R B(t)Rnm
{ei }ni=1 {ei }ni=1 {fi }m
i=1 {ei }ni=1
C(t) D(t)
(X, R) (Y, R) (U, R) (Y, R)
C(t)Rpn B(t)Rnm
{ei }ni=1 {gi }pi=1 {fi }m
i=1 {gi }pi=1 .
For the rest of these notes we will be interested in dynamical systems of the form
where x(t) Rn , u(t) Rm , and y(t) Rp denote the representations of elements of (X, R), (U, R),
p
and (Y, R) with respect to the bases {ei }ni=1 , {fi }m
i=1 , and {gi }i=1 respectively.
To ensure that the system (4.4)-(4.5) is well-posed we will from now on impose following assumption:
Assumption 4.1 A(), B(), C(), D() and u() are piecewise continuous.
Fact 4.1 Given u() : R Rm and (t0 , x0 ) R Rn there exists a unique solution x() : R Rn
and y() : R Rp for the system (4.4)-(4.5).
1 The complex numbers, C, can also be used as the field, at the expense of some additional complication in
dimension counting. For simplicity we will think of linear spaces as defined over the field of real numbers, unless
otherwise specified (e.g. for eigenvalue calculations).
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
61
Proof:(Sketch) Define p : Rn R Rn by
Exercise 4.1 Show that under Assumption 4.1 p satisfies the conditions of Theorem 3.7. (Take D
to be the union of the discontinuity sets of A(), B(), and u()).
The existence and uniqueness of x() : R Rn then follows by Theorem 3.7. Defining y() : R Rp
by
y(t) = C(t)x(t) + D(t)u(t) t R
completes the proof.
mapping the input trajectory u() : R Rm and initial condition (t0 , x0 ) R Rn to the state and
output at time t R respectively.
Theorem 4.1 Let Dx be the union of the discontinuity sets of A(), B() and u() and Dy the union
of the discontinuity sets of C(), D() and u().
where Rn = (0, . . . , 0) denotes the zero element of Rn and PC () : R Rm denotes the zero
function PC (t) = (0, . . . , 0) for all t R.
Proof: Part 1 follows from the definition of the solution. Part 4 follows from Part 3, by setting
u1 = U , u2 = u, x01 = x0 , x02 = Rn , and a1 = a2 = 1. Part 2 follows from Part 4, by noting
that s(t, t0 , , u) = s(t, t0 , , U ) + s(t, t0 , Rn , u) and s(t, t0 , , U ) : Rn Rn is a linear function
between finite dimensional linear spaces (and hence continuous by Corollary 3.2); the argument for
is similar. So we only need to establish Part 3.
Let x1 (t) = s(t, t0 , x01 , u1 ), x2 (t) = s(t, t0 , x02 , u2 ), x(t) = s(t, t0 , a1 x01 + a2 x02 , a1 u1 + a2 u2 ), and
(t) = a1 x1 (t) + a2 x2 (t). We would like to show that x(t) = (t) for all t R. By definition
Therefore x(t) = (t) since the solution to the linear ODE is unique. Linearity of follows from the
fact that C(t)x + D(t)u is linear in x and u.
For simplicity, from now on we will just use 0 to denote zero vectors both for the finite dimensional
space Rn and for the infinite dimensional space PC (). The interpretation should be clear from the
context.
Moreover, by Part 3 of Theorem 4.1, the zero input components s(t, t0 , x0 , U ) and (t, t0 , x0 , U )
are linear in x0 Rn . Therefore, in the basis used for the representation of A(), s(t, t0 , x0 , U ) has
a matrix representation (which will in general depend on t and t0 ). This representation is called the
state transition matrix and is denoted by (t, t0 ); note that the representation will depend on t and
t0 . In other words,
s(t,t0 ,,U )
(X, R) (X, R)
nn
(t,t0 )R
{ei }ni=1 {ei }ni=1
and, assuming s(t, t0 , x0 , 0) refers to the representation of the solution with respect to the basis
{ei }ni=1 ,
s(t, t0 , x0 , 0) = (t, t0 )x0 .
Exercise 4.2 Show that the representation of (t, t0 , , 0) : X Y with respect to the bases {ei }ni=1
and {gi }pi=1 is given by C(t)(t, t0 ); in other words
Therefore the state transition matrix (t, t0 ) completely characterizes the zero input state transi-
tion and output response. We will soon see that, together with the input trajectory u(), it also
characterizes the complete state transition and output response.
1. (, t0 ) : R Rnn is the unique solution of the linear matrix ordinary differential equation
(t, t0 ) = A(t)(t, t0 ) with (t0 , t0 ) = I.
t
Hence it is continuous for all t R and differentiable everywhere except at the discontinuity
points of A(t).
2. For all t, t0 , t1 R, (t, t0 ) = (t, t1 )(t1 , t0 ).
1
3. For all t1 , t0 R, (t1 , t0 ) is invertible and its inverse is [(t1 , t0 )] = (t0 , t1 ).
Proof: Part 1. Recall that s(t, t0 , x0 , U ) = (t, t0 )x0 is the solution to the linear differential
equation
x(t) = A(t)x(t) with x(t0 ) = x0 .
For i = 1, . . . , n consider the solution xi () = s(, t0 , xi (t0 ), U ) to the linear differential equation
starting at the representation xi (t0 ) = (0, . . . , 1, . . . , 0) Rn of the basis vector ei
0 0
.. ..
. .
xi (t) = A(t)xi (t) with xi (t0 ) =
1 = xi (t) = (t, t0 ) 1 .
. .
.. ..
0 0
Note that xi (t) is equal to the ith column of (t, t0 ). Putting the columns next to each other
(t, t0 ) = x1 (t) x2 (t) . . . xn (t)
shows that
d 1 d 2 d n
(t, t0 ) = x (t) x (t) . . . x (t)
t dt dt dt
= A(t)x1 (t) A(t)x2 (t) . . . A(t)xn (t)
= A(t) x1 (t) x2 (t) . . . xn (t) = A(t)(t, t0 ).
Moreover,
1 0 ... 0
0 1 ... 0
(t0 , t0 ) = x1 (t0 ) x2 (t0 ) . . . xn (t0 ) = .. .. . . .. = I R
nn
.
. . . .
0 0 ... 1
Part 1 follows.
Part 2. Consider arbitrary t0 , t1 R and let L(t) = (t, t0 ) and R(t) = (t, t1 )(t1 , t0 ). We would
like to show that L(t) = R(t) for all t R. Note that (by Part 1)
L(t1 ) = (t1 , t0 )
R(t1 ) = (t1 , t1 )(t1 , t0 ) = I (t1 , t0 ) = (t1 , t0 ).
Therefore L(t) = R(t) by existence and uniqueness of solutions of linear differential equations.
Part 3. First we show that (t, t0 ) is nonsingular for all t, t0 R. Assume, for the sake of contra-
diction, that it is not, i.e. there exists t, t0 R such that (t, t0 ) is singular. Then the columns
of (t, t0 ) are linearly dependent (Theorem 2.3) and there exists x0 Rn with x0 6= 0 such that
(t, t0 )x0 = 0.
Let x( ) = (, t0 )x0 . Notice that x(t) = (t, t0 )x0 = 0 and
d d
x( ) = (, t0 )x0 = A( )(, t0 )x0 = A( )x( ).
d d
Therefore x( ) is the unique solution to the differential equation:
d
x( ) = A( )x( ) with x(t) = 0. (4.6)
d
The function x( ) = 0 for all R clearly satisfies (4.6); therefore it is the unique solution to (4.6).
Let now = t0 . Then
0 = x(t0 ) = (t0 , t0 )x0 = I x0 = x0
which contradicts the fact that x0 6= 0. Therefore (t, t0 ) cannot be non-singular.
To determine its inverse, recall that for all t, t0 , t1 R,
In addition to these, the state transition matrix also has several other interesting properties, some
of which can be found in the exercises. We can now show that the state transition matrix (t, t0 )
completely characterizes the solution of linear time varying differential equations.
Proof: Several methods exist for proving this fact. The simplest is to invoke the rule of Leibniz for
differentiating integrals.
"Z # Z
b(t) b(t)
d d d
f (t, )d = f (t, )d + f (t, b(t)) b(t) f (t, a(t)) a(t).
dt a(t) a(t) t dt dt
(for the sake of comparison, notice that the fundamental theorem of calculus is the special case
a(t) = t0 , b(t) = t and f (t, ) = f ( ) independent of t.)
We start by showing that
Z t
s(t, t0 , x0 , u) = (t, t0 )x0 + (t, )B( )u( )d
t0
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
65
Rt
If we let L(t) = s(t, t0 , x0 , u) and R(t) = (t, t0 )x0 + t0 (t, )B( )u( )d we would like to show
that for all t R, L(t) = R(t). Notice that by definition L(t0 ) = s(t0 , t0 , x0 , u) = x0 and
d d
L(t) = s(t, t0 , x0 , u) = A(t)s(t, t0 , x0 , u) + B(t)u(t) = A(t)L(t) + B(t)u(t).
dt dt
We will show that R(t) satisfies the same differential equation with the same initial condition; the
claim then follows by the existence and uniqueness theorem.
Note first that
Z t0
R(t0 ) = (t0 , t0 )x0 + (t, )B( )u( )d = I x0 + 0 = x0 = L(t0 ).
t0
Let us now analyze the zero state transition and response in greater detail. By Theorem 4.1, the zero
state transition and the zero state response are both linear functions s(t, t0 , 0, ) : PC(R, Rm ) Rn
and (t, t0 , 0, ) : PC(R, Rm ) Rp respectively.
s(t,t0 ,0,)
(PC(R, Rm ), R) (Rn , R)
Z t
u() 7 (t, )B( )u( )d
t0
and
(t,t0 ,0,)
(PC(R, Rm ), R) (Rp , R)
Z t
u() 7 C(t) (t, )B( )u( )d + D(t)u(t).
t0
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
66
For j = 1, . . . , m, consider the zero state transition2 s(t, t0 , 0, (,) (t)fj ) under input (,) (t)fj .
Since s(t0 , t0 , 0, (,) (t)fj ) = 0 and the input is zero until t = ,
Therefore
lim s(t, t0 , 0, (,)(t)fj ) = (t, )B()fj .
0
we obtain
s(t, t0 , 0, (t)fj ) = (t, )B()fj Rm .
The statement is formal since to pass the limit inside the function we first need to ensure that the
function is continuous.
Exercise 4.4 We already know that the function s(t, t0 , 0, ) : PC(R, Rm ) Rn is linear. What
more do we need to check to make sure that it is continuous?
Moreover, strictly speaking (t) is not an acceptable input function, since it is equal to infinity
at t = and hence not piecewise continuous. Indeed, (t) is not a real valued function at all,
it just serves as a mathematical abstraction for an input pulse of arbitrarily small length. This
mathematical abstraction is known as the impulse function or the Dirac pulse. Even though in
practice the response of a real system to such an impulsive input cannot be observed, by applying
2 Strictly speaking, we should use the Rm representation of the basis vector f U instead of f itself. The reader
i i
is asked to excuse this slight abuse of the notation.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
67
as input piecewise continuous functions (,) (t) for small enough one can approximate the state
transition s(t, t0 , 0, (t)fj ) arbitrarily closely. Notice also that for t
i.e. the zero state transition due to the impulse (t)fj is also a zero input transition starting with
state B()fj at time .
Repeating the process for all {fj }m
j=1 leads to m vectors (t, )B()fj for j = 1, . . . , m. Order-
ing these vectors according to their index j and putting them one next to the other leads to the
impulse transition matrix, K(t, ) Rnm , defined by
(t, )B() if t
K(t, ) =
0 if t < .
The (i, j) element of K(t, ) contains the trajectory of state xi when the impulse function (t) is
applied to input uj . Note that, even though these elements cannot be measured in practice, the
impulse transition matrix is still a well defined, matrix valued function for all t, R.
Substituting s(t, 0, 0, (t)fj ) into the equation y(t) = C(t)x(t) + D(t)u(t) leads to
Note that, unlike K(t, ), H(t, ) is in general not a well defined matrix valued function since it in
general contains an impulse in its definition (unless of course D(t) = 0 for all t R).
1. Let f : Rn Rn be Lipschitz with Lipschitz constant K [0, +). For t R, let x(t) be
the solution of x(t) = f (x(t)), with x(0) = x0 . Let x Rn be such that f (x) = 0. Show that
||x(t) x|| eKt ||x0 x||, t R+ . (Here || || is the Euclidean norm on (Rn , R) for which K
is defined. Hint: use the Gronwall Lemma).
2. Let A : R+ Rnn , B : R+ Rnm and u : R+ Rm be piecewise continuous functions.
Show that, for any x0 Rn , the linear ODE
(
x(t) = A(t)x(t) + B(t)u(t), t R+ ,
x(0) = x0
has a unique solution x() : R+ Rn . (Hint: you may assume that if A(t) is piecewise
continuous then so is its induced norm).
Problem 4.4 (Linear ODEs) Let A() : R Rnn be piecewise continuous. Consider the fol-
lowing linear ODE:
x(t) = A(t)x(t), (4.7)
and let (t, t0 ) be the state transition matrix.
1. Show that t (t0 , t) = (t0 , t)A(t).
(Hint: differentiate the identity (t0 , t)(t, t0 ) = I.)
2. Let X0 be a convex set (i.e. x, x X0 x + (1 )x X0 , [0, 1]). Let s(t, t0 , x0 ) be
the solution of (1) associated with the initial condition x(t0 ) = x0 . Show that the set
X(t) = {s(t, t0 , x0 ) : x0 X0 }
Let us now turn to the special case where all the matrices involved in the system dynamics are
constant, in other words
with t R+ , x(t) Rn , u(t) Rm , y(t) Rp and A Rnn , B Rnm , C Rpn and D Rpm
are constant matrices.
A2 t2 Ak tk
eAt = I + At + + ...+ + . . . Rnn . (5.3)
2! k!
Proof: Exercise. Show that eA(tt0 ) satisfies the conditions of Part 1 of Theorem 4.2 by taking the
derivative of the expansion in (5.3) with respect to t. The result follows by uniqueness of solutions
of ordinary differential equations.
The above theorem together with the different properties of (t, t0 ) established in Chapter 4 imme-
diately lead to the following corollary.
Corollary 5.1 The state transition matrix, solution, impulse transition, and impulse response of a
time invariant linear system satisfy the following properties:
1
1. For all t, t1 , t0 R, eAt1 eAt2 = eA(t1 +t2 ) and eAt = eAt .
2. For all t, t0 R, (t, t0 ) = (t t0 , 0).
69
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
70
From the above it becomes clear that for linear time invariant systems the solution is independent
of the initial time t0 ; all that matters is how much time has elapsed since then, i.e. t t0 . Without
loss of generality we will therefore take t0 = 0 and write
Z t
x(t) = s(t, 0, x0 , u) = eAt x0 + eA(t ) B( )u( )d
0
Z t
At
y(t) = (t, 0, x0 , u) = Ce x0 + C eA(t )B( )u( )d + Du(t)
0
At
e B if t 0
K(t) = K(t, 0) =
0 if t < 0.
At
Ce B + D0 (t) if t 0
H(t) = H(t, 0) =
0 if t < 0.
Notice that in this case the integral that appears in the state transition and output response is
simply the convolution of the input u() with the impulse transition and impulse response matrices
respectively,
Definition 5.1 A matrix A Rnn is semi-simple if and only if its right eigenvectors {vi }ni=1 Cn
are linearly independent in the linear space (Cn , C).
Theorem 5.2 Matrix A Rnn is semi-simple if and only if there exists a nonsingular matrix
T Cnn and a diagonal matrix Cnn such that A = T 1 T .
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
71
Proof: () Take
T = [v1 v2 . . . vn ]1 Cnn
(which is invertible since {vi }ni=1 are linearly independent). Then
A = T 1 T.
() Assume that there exists matrices T Cnn nonsingular and Cnn diagonal such that
A = T 1 T AT 1 = T 1 .
with i C. Then
Awi = i wi
and therefore wi is a right eigenvector of A with eigenvalue i . Since T 1 is invertible its columns
(and eigenvectors of A) {wi }ni=1 are linearly independent.
If A is semi-simple, then its eigenvectors are linearly independent and can be used as a basis. Let
us see what the representation of the matrix A and the state transition matrix eAt with respect to
this basis are.
A
(Cn , C) (Cn , C)
ARnn
{ei }ni=1 {ei }ni=1 (basis leading to representation of A by A)
A=T AT 1 =
{vi }ni=1 {vi }ni=1 (eigenvector basis)
Recall that if x = (x1 , . . . , xn ) Rn is the representation of x with respect to the basis {ei }ni=1 its
representation with respect to the complex basis {vi }ni=1 will be the complex vector x = T x Cn .
The above formula simply states that
x = T x = T Ax = T AT 1 x = x.
Therefore, if A is semi-simple, its representation with respect to the basis of its eigenvectors is the
diagonal matrix of its eigenvalues. Notice that even though A is a real matrix, its representation
is in general complex, since the basis {vi }ni=1 is also complex.
What about the state transition matrix?
Proof: Exercise. Show that Ak = T 1 k T and substitute into the expansion (5.3). Or differentiate
(??) and invoke existence-uniqueness.
In other words:
s(t,0,,U )
(Rn , R) (Rn , R)
x0 7 x(t)
eAt Rnn
{ei }ni=1 {ei }ni=1 (basis leading to representation of A by A)
t At 1
e =T e T
{vi }ni=1 {vi }ni=1 (eigenvector basis).
Once again, note that the matrices T , T 1 and et will general be complex. Fact 5.1, however,
shows that when taking their product the imaginary parts will miraculously all cancel and we will
be left with a real matrix.
Fact 5.1 shows that if a matrix is semi-simple the calculation of the matrix exponential is rather
straightforward. It would therefore be desirable to establish conditions under which a matrix is
semisimple.
Definition 5.2 A matrix A Rnn is simple if and only if its eigenvalues are distinct, i.e. i 6= j
for all i 6= j.
Proof: Assume, for the sake of contradiction, that A Rnn is simple, but not semi-simple. Then
i 6= j for all i 6= j but {vi }ni=1 are linearly dependent in (Cn , C). Hence, there exist a1 , . . . , an C
not all zero, such that
Xn
ai vi = 0.
i=1
Without loss of generality, assume that a1 6= 0 and multiply the above identity by (A 2 I)(A
3 I) . . . (A n I) on the left. Then
n
X
a1 (A 2 I)(A 3 I) . . . (A n I)v1 + ai (A 2 I)(A 3 I) . . . (A n I)vi = 0.
i=2
Concentrating on the first product and unraveling it from the right leads to
which leads to
n
X
a1 (1 2 )(1 3 ) . . . (1 n )v1 + ai (i 2 )(i 3 ) . . . (i n )vi = 0.
i=2
Each term of the sum on the right will contain a term of the form (i i ) = 0. Hence the sum on
the right is zero, leading to
a1 (1 2 )(1 3 ) . . . (1 n )v1 = 0.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
73
But, v1 6= 0 (since it is an eigenvector) and a1 6= 0 (by the assumption that {vi }ni=1 are linearly
dependent) and (1 i ) 6= 0 for i = 2, . . . , n (since the eigenvalues are distinct). This leads to a
contradiction.
In addition to simple matrices, other matrices are simple, for example diagonal matrices and orthog-
onal matrices. In fact one can show that semi-simple matrices are dense in Rnn , in the sense
that
A Rnn , > 0, A Rnn semi-simple : kA A k < .
(Here k k denotes any matrix norm, anyway they are all equivalent). The reason is that for a matrix
to be non semi-simple the matrix of its eigenvectors T Cnn must be singular, i.e. we must have
Det[T ] = 0. But this is a fragile condition, as arbitrary small perturbations of the matrix A will in
general lead to Det[T ] 6= 0. In fact it is not hard to convince ourselves that almost all matrices
in Rnn are semi-simple, since the condition Det[T ] = 0 imposes a single constraint on the n2
dimensional space of complex matrices1 .
In summary, almost all matrices are semi-simple, though not all.
is not semi-simple. Its eigenvalues are 1 = 2 = 3 = (hence it is not simple) but there is only
one eigenvector
1 0 x1 x1 x1 + x2 = x1 1
0 1 x2 = x2 x2 + x3 = x2 v1 = v2 = v3 = 0
0 0 x3 x3 x3 = x3 0
Note that A1 has the same eigenvalues as the matrices
1 0 0 0
A2 = 0 0 and A3 = 0 0
0 0 0 0
A2 is also not semi-simple, it has only two linearly independent eigenvectors v1 = v2 = (1, 0, 0) and
v3 = (0, 0, 1). A3 is semi-simple, with the canonical basis as eigenvectors. Notice that none of the
three matrices is simple.
Example (Non semi-simple matrices (cont.)) For the matrices considered above
Notice that in all three cases the algebraic multiplicity of the eigenvalue is 3. The problem with
A1 and A2 is that the geometric multiplicity of is smaller than its algebraic multiplicity. This
implies that there are not enough eigenvectors associated with eigenvalue (in other words linearly
independent vectors in the null space) to form a basis.
The above example suggests that the reason matrices fail to be semi-simple is a discrepancy between
their algebraic and geometric multiplicities. To alleviate this problem we need to complete the basis
with additional vectors. For this purpose we consider the so called generalized eigenvectors.
1. v 0 = 0,
2. {v j }j=1 are linearly independent, and
3. [A I]v j = v j1 for j = 1, 2, . . . , .
A Jordan chain {v j }j=0 is called maximal if it cannot be extended, i.e. there does not exist v Cn
linearly independent from {v j }j=1 such that [A I]v = v . The elements of all the maximal Jordan
chains at are the called generalized eigenvectors of .
Fact 5.2 Let {v j }j=0 Cn be a Jordan chain of length at eigenvalue C of the matrix
A Rnn :
Example (Non semi-simple matrices (cont.)) For the matrices considered above, A1 has one
maximal Jordan chain of length = 3 at , with
A2 has two maximal Jordan chains at , one of length 1 = 2 and the other of length 2 = 1, with
Theorem 5.4 Assume that the matrix A Rnn has k linearly independent eigenvectors v1 , . . . , vk
Cn with corresponding maximal Jordan chains {vij }j=0
i
Cn , i = 1, . . . , k. Then
Pk
1. The matrix v11 . . . v11 . . . vk1 . . . vkk . . . Cnn is invertible. In particular, i=1 i = n.
2. A can be written in the Jordan canonical form A = T 1 JT where J Cnn is block-diagonal
i 1 0 . . . 0 0
J1 0 . . . 0 0 i 1 . . . 0
0 J2 . . . 0 0
.. Ci i , i = 1, . . . , k
J = . . . . Cnn , Ji = ... .. .. . .
. .
.
. .. .
.. .. . . ..
0 0 0 . . . i 1
0 0 . . . Jk
0 0 0 ... 0 i
The proof of this fact is rather tedious and will be omitted, see [?]. Notice that there may be multiple
Jordan chains for the same eigenvalue , in fact their number will be the same as the number of
linearly independent eigenvectors associated with . If k = n (equivalently, all Jordan chains have
length 1) then the matrix is semi-simple, T 1 is the matrix of eigenvectors of A, and J = .
Example (Non semi-simple matrices (cont.)) In the above example, the three matrices A1 ,
A2 and A3 are already in Jordan canonical form. A1 comprises one Jordan block of size 3, A2 two
Jordan blocks of sizes 2 and 1 and A3 three Jordan blocks, each of size 1.
So the computation of the matrix exponential becomes easy once again. Notice that if k = n then we
are back to the semi-simple case. As in the semi-simple case, all matrices involved in the product will
in general be complex. However, the theorem ensures that when taking the product the imaginary
parts will cancel and the result will be a real matrix.
Notice that in all cases eAt consists of linear combinations of elements of the form
ei t , tei t , . . . , ti 1 ei t
for i Spec[A] and i the length of the longest Jordan chain at i . In other words
X
eAt = (t)et (5.5)
Spec[A]
where for each Spec[A], (t) C[t]nn is a matrix of polynomials of t with complex coefficients
and degree at most equal to the length of the longest Jordan chain at . In particular, if A is semi-
simple all Jordan chains have length equal to 1 and the matrix exponential reduces to
X
eAt = et
Spec[A]
where Cnn are constant complex matrices. Notice again that even though in general both
the eigenvalues and the coefficients of the corresponding polynomials (t) are complex numbers,
because the eigenvalues appear in complex conjugate pairs the imaginary parts for the sum cancel
out and the result is the real matrix eAt Rnn .
where the integral is interpreted element by element and we assume that it is well defined (for a
careful discussion of this point see, for example, [?]). The Laplace transform L{f (t)} transforms the
real matrix valued function f (t) Rnm of the real number t R+ to the complex matrix valued
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
77
function F (s) Cnm of the complex number s C. The inverse Laplace transform L1 {F (s)}
performs the inverse operation; it can also be expressed as an integral, even though in the calculations
considered here one mostly encounters functions F (s) that are Laplace recognizable transforms of
known functions f (t); in particular the functions F (s) will typically be proper rational functions of
s whose inverse Laplace transfrom can be computed by partial fraction expansion.
Fact 5.3 The Laplace transform (assuming that it is defined for all functions concerned) has the
following properties:
Fact 5.4 For all A Rnn and t R+ , L eAt = (sI A)1 .
The claim follows by multiplying on the left by (sI A); notice that the matrix is invertible for all
s C, except the eigenvalues of A.
Let us look somewhat more closely to the structure of the Laplace transform, (sI A)1 , of the
state transition matrix. This is an n n matrix of strictly proper rational functions of s, which as
we saw in Chapter 2 form a sub-ring of (Rp (s), +, ). To see this recall that by definition
Adj[sI A]
(sI A)1 =
Det[sI A]
The denominator is simply the characteristic polynomial,
M0 = I
Mi = Mi1 A + i I for i = 1, . . . , n 1
Mn1 A + n I = 0.
Since the last identity must hold for all s C the coefficients of the two polynomials on the left and
on the right must be equal. Equating the coefficients for sn leads to the formula for M0 , for sn1
leads to the formula for M1 , etc.
The theorem provides an easily implementable algorithm for computing the Laplace transform of
the state transition matrix without having to invert any matrices. The only thing that is needed is
the computation of the characteristic polynomial of A and some matrix multiplications. In addition,
the following useful fact about square matrices can be deduced as a corollary.
Theorem 5.7 (Cayley-Hamilton) Every square matrix A Rnn satisfies its characteristic poly-
nomial, i.e.
A (A) = An + 1 An1 + . . . + n I = 0 Rnn .
Proof: By the last equation of Theorem 5.6, Mn1 A + n I = 0. From the next to last equation,
Mn1 = Mn2 A + n1 I; substituting this into the last equation leads to
Mn2 A2 + n1 A + n I = 0.
Substituting Mn2 from the third to last equation, etc. leads to the claim.
Corollary 5.2 Let A Rnn be a square matrix. For any k N, Ak can be written as a linear
combination of {I, A, A2 , . . . , An1 }.
The proof is left as an exercise. Stated another way, the corollary shows that all powers of
an n-dimensional square matrix live in a low-dimansional subspace of the n2 dimensional linear
space of square matrices; this subspace has dimension at most n and is spanned by the matrices
{I, A, A2 , . . . , An1 }.
The Laplace transform can also be used to compute the response of the system. Notice that taking
Laplace transforms of both sides of the differential equation governing the evolution of x(t) leads to
L
x(t) = Ax(t) + Bu(t) = sX(s) x0 = AX(s) + BU (s)
leading to
X(s) = (sI A)1 x0 + (sI A)1 BU (s). (5.7)
How does this relate to the solution of the differential equation that we have already derived? We
have shown that Z t
At
x(t) = e x0 + eAt Bu( )d = eAt x0 + (K u)(t)
0
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
79
K(t) = eAt B
By (5.6)
Y (s) = C(sI A)1 x0 + (sI A)1 BU (s) + DU (s)
which for x0 = 0 (zero state response) reduces to
Comparing equation (5.7) with the zero state response that we computed earlier
Z t
y(t) = C eA(t ) Bu( )d + Du(t) = (H u)(t)
0
it is clear that the transfer function is the Laplace transform of the impulse response H(t) of the
system
G(s) = L{H(t)} = L{CeAt B + D0 (t)} = C(sI A)1 B + D.
Definition 5.6 The poles of the system are the values of s C are the roots of the denominator
polynomial of G(s).
From equation (5.9) it becomes apparent that all poles of the system are eigenvalues (i.e. are
contained in the spectrum) of the matrix A. Note, however, that not all eigenvalues of A are
necessarily poles of the system, since there may be cancellations of common factors in the numerator
and denominator when forming the fraction (5.9). It turns out that such cancellations are related
to the controllability and observability properties of the system. We will return to this point in
Chapter 8, after introducing these notions.
1. Derive the representation of the system w.r.t. bases {ui }m xi }ni=1 , {yi }pi=1 .
i=1 , {e
2. Compute the transition map (t,e t0 ) and the impulse response matrix H(t,
e ) w.r.t. the new
representation. How do they compare with the corresponding quantities (t, t0 ) and H(t, )
in the original representation?
Problem 5.2 (Time-invariant Systems) Consider the linear time-invariant system of Problem 5.1.
1. Show that k
X A(t t0 )
(t, t0 ) = exp A(t t0 ) = .
k!
k=0
2. Given two matrices A1 , A2 Rnn show that, if A1 A2 = A2 A1 , then A2 exp(A1 t) = exp(A1 t)A2
and exp (A1 + A2 )t = exp(A1 t) exp(A2 t). Also show that these properties may not hold if
A1 A2 6= A2 A1 .
3. Show that the impulse response matrix satisfies H(t, ) = H(t , 0) (i.e. it depends only on
the difference t ).
state and the output of system () sampled at times tk . Show that there exist matrices Ak ,
Bk , Ck and Dk such that
xk+1 = Ak xk + Bk uk ,
()
yk = Ck xk + Dk uk .
2. Now assume that () is time-invariant, i.e. A(t), B(t), C(t), D(t) = (A, B, C, D), t t0 ,
and that there exists a fixed T > 0, tk+1 tk = T , k. Provide simplified expressions for Ak ,
Bk , Ck and Dk and show that they are independent of k.
Problem 5.4 (Realization) Consider the following n-th order scalar differential equation with
constant coefficients:
y (n) (t) + a1 y (n1) (t) + . . . + an1 y (1) (t) + an y(t) = u(t), t R+ , (5.12)
where y (i) (t) denotes the i-th derivative of y at t, {ai } R and u() : R+ R is a piecewise
continuous input. Show that (5.11) can be put in the form (5.10) for an appropriate definition of
the state x(t) Rn and of matrices A, B, C, D.
where the elements not shown are zeroes. (Hint: in 3, consider the decomposition A = I + N , and
Pk k!
make use of (I + N )k = i=0 i!(ki)! (I)i N ki .)
Problem 5.6 (Jordan blocks and matrix exponential) For i = 1, . . . , m, let i C ni ni . De-
fine n = n1 + . . . + nm and the block diagonal matrix
1
2
= diag(1 , 2 , . . . , m ) = .. Rnn . ()
.
m
1. Show that exp() = diag exp(1 ), exp(2 ), . . . , exp(m ) .
2. Compute exp(i ) for the following definitions of i : (assume each entry is real)
1
1
.. ..
(a) i = .. ; (b) = . . ;
. i
1
ni
(c) i = R22 ,
where the elements not shown are zeroes. [Hint: in (b), consider the decomposition = I +N ,
P
and make use of (I + N )k = ki=0 i!(ki)!
k!
(I)i N ki .]
Problem 5.7 (Modal Analysis) For t R+ , consider the ODE x = Ax, x(0) = x0 . Compute
exp(At) for the cases listed below. In each case provide a rough plot of the parametric curves x2 (t)
vs. x1 (t) for some initial conditions x0 .
0
1. A= , for < 0 and > 0.
0
2. A= , for < 0 and > 0.
2 1
3. A= 1 , for 1 < 0 < 2 , 1 , 2 < 0 and 1 , 2 > 0.
0 2
1 1
4. A= , for < 0, > 0 and = 0.
1 +1
Problem 5.8 (Matrix powers) Consider A Rnn , C Rpn and B Rnm . Show that
Problem 5.9 (Nilpotent matrices) Show that a matrix A Rnn is nilpotent (i.e. there exists
k N such that Ak = 0) if and only if all of its eigenvalues are equal to zero. Show further that in
this case k n.
Problem 5.10 (Transfer function) Let A Rnn , B Rnm , C Rpn and D Rpm . Let
u() : R+ U be piecewise continuous. For t R+ , consider the linear time-invariant system:
(
x(t) = Ax(t) + Bu(t),
y(t) = Cx(t) + Du(t).
1. Consider a change of basis x = T x. Compute the transfer function with respect to the new
basis and compare it to the transfer function in the original basis.
2. Assume now that A is semisimple. Use your answer in part 1 to provide a simple formula for
G(s) in terms of the eigenvalues of A.
Chapter 6
Stability
Stability addresses the question of what happens to the solutions of this system as time, t, increases.
Do they remain bounded, will they get progressively smaller, or will they diverge to infinity. Stability
deals first and foremost with the properties of the differential equation. We will therefore ignore the
output equation for the time being. We will also start with the zero input case (u = U , or u(t) = 0
for all t R), i.e. by considering the solutions of
We will then return to inputs and outputs in Section 6.4. As we did for the definition of the solutions
of the differential equation, we will start by considering general nonlinear systems (Section 6.1, then
specialize to linear time varying systems (Section 6.2), then specialize further to linear time invariant
systems (Section 6.3).
Definition 6.1 A state x Rn is called an equilibrium of system (6.2) if and only if p(x, t) = 0 for
all t R.
83
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
84
d d
s(t, t0 , x) = x = 0 = p(x, t) = p(s(t, t0 , x), t).
dt dt
The conclusion follows by existence and uniqueness of solutions.
The fact shows that a solution which passes through an equilibrium, x, at some point in time is
doomed to stay on the equilibrium for all times. We call this constant solution the equilibrium
solution defined by the equilibrium x.
What if a solution passes close to the equilibrium, but not exactly through it? Clearly such a solution
will no longer be identically equal to the equilibrium, but will it move away from the equilibrium,
or will it remain close to it, even converge to it, and if so at what rate?
To answer these questions we first need to fix a norm on Rn to be able to measure distances. Any
norm will do since they are all equivalent, for simplicity we will use the Euclidean norm throughout.
Equiped with this norm, we can now formalize the above questions in the following definition.
1. Stable if and only if for all t0 R, and all > 0, there exists > 0 such that
4. Locally asymptotically stable if and only if it is stable and for all t0 R there exists M > 0
such that
kx0 xk M lim ks(t, t0 , x0 ) xk = 0.
t
5. Globally asymptotically stable if and only if it is stable and for all (t0 , x0 ) R Rn
lim ks(t, t0 , x0 ) xk = 0.
t
6. Locally exponentially stable if and only if for all t0 R there exist , m, M > 0 such that for
all all t t0
kx0 xk M ks(t, t0 , x0 ) xk mkx0 xke(tt0 ) .
7. Globally exponentially stable if and only if for all t0 R there exist , m > 0 such that for all
x0 Rn and all t t0
Special care is needed in the above definition: The order of the quantifiers is very important. Note
for example that the definition of stability implicitly allows to depend on t0 and ; one sometimes
writes (t0 , ) to highlight this dependence. On the other hand, in the definition of uniform stability
can depend on but not on t0 , i.e. the same must work for all t0 ; one sometimes uses the
notation () to highlight this fact. Likewise, the definition of global exponential stability requires
and m to be independent of x0 , i.e. the same and m must work for all x0 Rn ; a variant of
this definition where m and are allowed to depend on x0 is sometimes referred to as semi-global
exponential stability.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
85
The definition distinguishes stability concepts along three axes. The most fundamental distinction
deals with the convergence of nearby solutions to the equilibrium. The equilibrium is called unstable
if we cannot keep solutions close to it by starting sufficiently close, stable if we can keep solutions
as close as we want by starting them sufficiently close, asymptotically stable if in addition nearby
solutions converge to the equilibrium, and exponentially stable if they do so at an exponential
rate. The second distinction deals with how these properties depend on the starting time, t0 : For
uniform stability the starting time is irrelevant, the property holds the same way irrespective of
when we look at the system. The third distinction deals with how these properties depend on the
starting state, x0 : Local implies that the property holds provided we start close enough to the
equilibrium, whereas global requires that the property holds irrespective of where we start. Note
that this distinction is irrelevant for stability and uniform stability, since the conditions listed in
the definition are required to hold provided we start close enough. One can also combinatorially
mix these qualities to define other variants of stability notions: Uniform local asymptotic stability
(where the equilibrium is uniformly stable and the convergence rate is independent of the starting
time), uniform global exponential stability, etc. We will not pursue these variants of the definitions
here, since most of them turn out to be irrelevant when dealing with linear systems.
It is easy to see that the notions of stability introduced in Definition 6.2 are progressively stronger.
Fact 6.2 Consider an equilibrium of system 6.2. Then the following statements are true:
Proof: Most of the statements are obvious from the definition. Asymptotic stability requires
stability, global asymptotic stability implies that the conditions of local asymptotic stability hold
for any M > 0, etc. The only part that requires any work is showing that local/global exponential
stability implies local/global asymptotic stability.
Consider a globally1 exponentially stable equilibrium x, i.e. assume that for all t0 there exist
, m > 0 such that for all x0 , ks(t, t0 , x0 ) xk mkx0 xke(tt0 ) for all t t0 . For t0 R and
> 0 take = /m. Then for all x0 Rn such that kx0 xk < and all t t0
ks(t, t0 , x0 ) xk mkx0 xke(tt0 ) < me(tt0 ) = e(tt0 ) .
Hence the equilibrium is stable. Moreover, since by the properties of the norm ks(t, t0 , x0 ) xk 0,
0 lim ks(t, t0 , x0 ) x0 k lim mkx0 xke(tt0 ) = 0.
t t
It is also easy to see that the stability notions of Definition 6.2 are strictly stronger one from the
other; in other words the converse implications in the table of Fact 6.2 are in general not true. We
show this through a series of counter-examples.
Example (Stable, non-uniformly stable equilibrium) For x(t) R consider the linear, time
varying system
2t
x(t) = x(t) (6.3)
1 + t2
Exercise 6.1 Show that the system has a unique equilibrium at x = 0. Show further that
1 + t20
s(t, t0 , x0 ) = x0
1 + t2
by differentiating and invoking existence and uniqueness of solutions.
Typical trajectories of the system for x0 = 1 and different values of t0 are shown in Figure ??.
It is easy to see that x = 0 is a stable equilibrium. Indeed, given t0 R and > 0 let = /(1 + t20).
Then for all x0 Rn such that kx0 k < ,
1 + t20
ks(t, t0 , x0 )k =
x0
< .
1+t 2
1 + t2
However, the equilibrium is not uniformly stable: For a given we cannot find a that works for all
t0 R. To see this, notice that for t0 0, ks(t, t0 , x0 )k reaches a maximum of (1 + t20 )kx0 k at t = 0.
Hence to ensure that ks(t, t0 , x0 )k < we need to ensure that (1 + t20 )kx0 k < which is impossible
to do by restricting
p x0 alone; for any 0 < < and kx0 k < we can make ks(0, t0 , x0 )k > by
taking t0 < / 1.
Example (Stable, non asymptotically stable equilibrium) For x(t) R2 consider the linear,
time invariant system
0 x01
x(t) = x(t) with x(0) = x0 = .
0 x02
Since the system is linear time invariant we can take t0 = 0 without loss of generality.
Exercise 6.2 Show that the system has a unique equilibrium x = 0. Show further that
cos(t) sin(t)
(t, 0) =
sin(t) cos(t)
Typical trajectories of the system for different values of x0 are shown in Figure ??. Clearly
x01 cos(t) x02 sin(t)
x(t) = (t, 0)x0 = .
x01 sin(t) + x02 cos(t)
Exercise 6.3 Show that linear time invariant system x(t) = 0 with x(t) R also has a stable, but
not asymptotically stable equilibrium x = 0. Does this system have any other equilibria? Are they
stable? Asymptotically stable?
for all t0 , x0 . Hence the equilibrium is globally asymptotically stable. It is, however not exponentially
stable (not even locally). Assume, for the sake of contradiction, that for all t0 R there exist
, m, M > 0 such that for all x0 R with kx0 k M and all t t0 ,
Notice that so far all the counter-examples have involved linear systems (sometimes time varying).
As we will see in the next section, for linear systems local and global asymptotic stability are
equivalent concepts. To distinguish between them we therefore need to resort to non-linear systems.
x1 (t) = x2 (t)
x2 (t) = x1 (t) x1 (t)3 x2 (t).
It is easy to see that this system has three equilibria, at (0, 0), (1, 0) and (1, 0) respectively.
Exercise 6.4 Verify that these are indeed the only equilibria of this system. Is the system linear or
nonlinear? What is the function p in x(t) = p(x(t), t)? Is it globally Lipshitz in its first argument?
Clearly none of the three equilibria can be globally asymptotically stable: It is impossible for all
trajectories to converge to a particular equilibrium, since those starting at another equilibrium will
stay put. To study whether some of these equilibria are stable or locally asymptotically stable one
can compute the linerization of the system about each of the equilibria and study the stability of
the resulting linear system using the methods presented later in this chapter. It is reasonable to
assume that nearby the equilibrium, where the terms neglected in the linearization are small, the
behavior of the nonlinear system will be similar to that of its linearization. Hence if the linearization
is asymptotically stable one whould expect the equilibrium to be locally asymptotically stable for
the nonlinear system. This argument can in fact be formalized, leading to the so-called Lyapunov
Indirect Method for checking local asymptotic stability, or instability of nonlinear systems. The
interested reader is referred to [?] or [?] for more information on this topic.
For this system, linearization suggests that one would expect the equilibria (1, 0) and (1, 0) to be
locally asymptotically stable and the equilibrium (0, 0) to be unstable (Problem ??). This can be
visually confirmed by simulating the system for various initial conditions and plotting the trajectories
x(t). The most informative way of doing this is to generate a parametric plot of x1 (t) against x2 (t)
parametrized by t. This so-called phase plane plot for this system is shown in Figure ??.
Before restricting our attention to linear system we point out two more general facts about the
stability concepts introduced in Definition 6.2. The first is an intimate relation between stability
and continuity. To expose this link we need to think of the function mapping initial conditions to
state trajectories from the initial time t0 onwards. Since state trajectories are continuous functions
of time, for each t0 R one can think of this function as a map between the state space Rn and the
space of continuous functions C([t0 , ), Rn )
s(, t0 , ) : Rn C([t0 , ), Rn )
x0 7 {s(, t0 , x0 ) : [t0 , ) Rn }.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
88
The strange notation is meant to alert the reader to the fact that we consider s(, t0 , ) for fixed t0
as a function mapping a vector (denoted by the placeholder ) to a function of time (denoted by
the placeholder , left over after x0 is substituted for ).
Recall that for the stability definitions we have equipped Rn with a norm k k. We now equip
C([t0 , ), Rn ) with the corresponding infinity norm
where we include t0 in the notation to make the dependence on initial time explicit. Notice that the
first norm in Equation (6.4) is a norm on the infinite dimensional function space (i.e., s(, t0 , x0 )
C([t0 , ), Rn ) is though of as a function of time), whereas the second norm is a norm on the finite
dimensional state space (i.e., s(t, t0 , x0 ) Rn is the value of this function for the specific time
t [t0 , )).
Fact 6.3 An equilibrium, x, of system (6.2) is stable if and only if for all t0 R the function
s(, t0 , ) mapping the normed space (Rn , k k) into the normed space (C([t0 , ), Rn ), k kt0 , ) is
continuous at x.
Proof: The statement is effectively a tautology. Fix t0 R and recall that, according to Defini-
tion 3.6, s(, t0 , ) is continuous at x if and only for all > 0 there exists > 0 such that
A similar relation between uniform stability and uniform continuity (where the above is indepen-
dent of t0 ) can also be derived in the same way.
The second general fact relates to the possible rate of convergence. The strongest notion of stability
in Definition 6.2, namely exponential stability, requires that solutions converge to the equilibrium
exponentially (i.e. rather quickly) in time. Could they converge even faster? Could we, for example,
introduce another meaningful stability definition that requires solutions to coverge with a rate of
2
et for some > 0? And if not, can we at least increase in the exponential convergence? The
following fact reveals that Lipschitz continuity imposes a fundamental limit on how fast convergence
can be.
Fact 6.4 Let x be an equilibrium of system (6.2) and assume that there exists k > 0 such that for
all x, x Rn , kp(x, t) p(x , t)k kkx x k. Then for all t0 R and all t t0
Proof: If x0 = x the claim is trivially true, we therefore restrict attention to the case x0 6= x. Note
that in this case we must have s(t, t0 , x0 ) 6= x for all t; if s(t, t0 , x0 ) = x for some t then s(t, t0 , x0 )
must be the equilibrium solution and s(t, t0 , x0 ) = x for all t which, setting t = t0 contradicts the
fact that x0 6= x.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
89
Recall that for simplicity we are using the Euclidean norm. Hence ks(t, t0 , x0 ) xk2 = (s(t, t0 , x0 )
x)T (s(t, t0 , x0 ) x) and
d
ks(t, t0 , x0 ) xk2 = d s(t, t0 , x0 )T (s(t, t0 , x0 ) x) + (s(t, t0 , x0 ) x)T d s(t, t0 , x0 )
dt dt dt
= p(s(t, t0 , x0 ), t)T (s(t, t0 , x0 ) x) + (s(t, t0 , x0 ) x)T p(s(t, t0 , x0 ), t)
p(s(t, t0 , x0 ), t)T (s(t, t0 , x0 ) x) + (s(t, t0 , x0 ) x)T p(s(t, t0 , x0 ), t)
kp(s(t, t0 , x0 ), t)T k ks(t, t0 , x0 ) xk + k(s(t, t0 , x0 ) x)T k kp(s(t, t0 , x0 ), t)k
= 2ks(t, t0 , x0 ) xk kp(s(t, t0 , x0 ), t)k
= 2ks(t, t0 , x0 ) xk kp(s(t, t0 , x0 ), t) p(x, t)k
2kks(t, t0 , x0 ) xk ks(t, t0 , x0 ) xk.
or in other words
d
kks(t, t0 , x0 ) xk ks(t, t0 , x0 ) xk kks(t, t0 , x0 ) xk.
dt
Applying the Gronwall Lemma (Theorem 3.9) to the right inequality leads to
From the right inequality (adapting the steps of the proof of the Gronwall Lemma) we have
d d d
ks(t, t0 , x0 ) xkek(tt0 ) = (ks(t, t0 , x0 ) xk) ek(tt0 ) + ks(t, t0 , x0 ) xk ek(tt0 )
dt dt dt
kks(t, t0 , x0 ) xkek(tt0 ) + ks(t, t0 , x0 ) xkkek(tt0 ) = 0.
which leads to
ks(t, t0 , x0 ) xk kx0 xkek(tt0 ) .
The fact reveals that convergence to an equilibrium can be at most exponential and the rate can
be at most equal to the Lipschitz constant of p. On the bright side, the fact also shows that if
an equilibrium is unstable divergence cannot be any faster than exponential. Even though a fixed
Lipshitz constant is assumed to simplify the proof it is easy to see that the claim still holds if the
Lipschitz constant is time varying but bounded from above and below; one simply needs to replace
k by its lower bound in the left inequality and its upper bound in the right inequality. The lower
bound on the Lipschitz constant also provides a bound for the rate of exponential convergence.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
90
is the solution of (6.1). Another way to think of this observation is that if x(t0 ) = 0 for some
t0 R+ , then
x(t0 ) = A(t)x(t0 ) = 0
therefore the solution of the differential equation does not move from 0. Either way, a zero input
solution that passes through the state x(t0 ) = 0 at some time t0 R+ will be identically equal to
zero for all times.
Exercise 6.5 Can there be other x0 6= 0 such that s(t, t0 , x0 , 0) = x0 for all t R+ ?
1. Stable if and only if for all t0 R, there exists M > 0 such that
ks(t, t0 , x0 , U )k M kx0 k t t0 , x0 Rn .
Proof: Part 1: By Fact ??, the equilibrium x = 0 is stable if for all t0 R the function
is continuous at x0 = 0. Recall that for all t0 R+ this function is linear in x0 (by Theorem 4.1).
Therefore (by Theorem 3.5) the function is continuous at x0 = 0 if and only if it is continuous at all
x0 Rn , if and only if its induced norm is finite. In other words, the equilibrium solution is stable
if and only if for all t0 R+ there exists M > 0 such that
ks(, t0 , x0 , U )kt0 ,
sup M ks(, t0 , x0 , U )kt0 , M kx0 k2 x0 Rn
x0 6=0 kx0 k2
ks(t, t0 , x0 , U )k2 M kx0 k2 t t0 , x0 Rn .
Part 2: Assume first that limt (t, 0) = 0 and show that the equilibrium solution is asymptotically
stable. For all (x0 , t0 ) Rn R+
ks(t, t0 , x0 , 0)k2 = k(t, t0 )x0 k2 = k(t, 0)(0, t0 )x0 k2 k(t, 0)k k(0, t0 )x0 k2
(where we use the induced 2-norm for the matrix (t, 0)). Therefore, for all x0 Rn , t0 R,
lim ks(t, t0 , x0 , U )k = 0.
t
To establish that the equilibrium is asymptotically stable it therefore suffices to show that it is stable.
Recall that the function s(, t0 , x0 , U ) = (, t0 )x0 is continuous (by Theorem 4.1). Therefore, since
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
91
limt ks(t, t0 , x0 , 0)k = 0, we must have that s(, t0 , x0 , 0) is bounded, i.e. there exists M > 0 such
that for all x0 Rn
s(t, t0 , x0 , U ) M kx0 k t t0 .
Hence the equilibrium solution is stable by Part 1.
Conversely, assume that the equilibrium solution is locally asymptotically stable. Then there exists
M > 0 such that for all t R and all x0 Rn with kx0 k M ,
lim s(t, t0 , x0 , 0) = lim (t, t0 )x0 = lim ((t, 0)(0, t0 )x0 ) = lim ((t, 0)) (0, t0 )x0 = 0 (6.5)
t t t t
From this theorem one can further conclude that when it comes to linear systems (even linear time
varying systems) local stability notions are equivallent to global ones.
Fact 6.6 The equilibrium x = 0 of the linear time invariant system x(t) = Ax(t) is uniformly stable
if and only if it is stable.
Theorem 6.2 For linear time invariant systems the following statements are equivalent:
Lemma 6.1 For all > 0 there exists m > 0 such that for all t R+
keAt k me(+)t
where k k denotes any induced norm on Rnn and = max{Re[] | Spec[A]}.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
92
Proof: Recall that the existence of the Jordan canonical form implies that
X
eAt = (t)et
Spec[A]
where (t) C[t]nn are n n matrices of polynomials in t with complex coefficients. Consider
(for simplicity) the infinity induced norm k k for Cnn . Then
X
X
keAt k =
(t)e t
k (t)k et .
Spec[A]
Spec[A]
Recall that k (t)k is the maximum among the rows of (t) of the sum of the magnitudes of the
elements in the row. Since all entries are polynomial, then there exists a polynomial p (t) R[t]
such that p (t) k (t)k for all t R+ . If we define the polynomial p(t) R[t] by
X
p(t) = p (t),
Spec[A]
then
keAt k p(t)et . (6.6)
t
Since p(t) is a polynomial in t R+ , for any > 0 the function p(t)e is continuous and
limt p(t)et = 0. Therefore p(t)et is bounded for t R+ , i.e. there exists m > 0 such
that p(t)et m for all t R+ . Substituting this into equation (6.6) leads to
Proof: (Of Theorem 6.2) We have already seen that 2 1 (Fact 6.2).
3 2: If all eigenvalues have negative real part then
Consider (0, ) and set = ( + ) > 0. By Lemma 6.1 there exists m > 0 such that
keAt k met .
s(t, 0, v, 0) = et v Cn .
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
93
is bounded away from zero; it is either constant (if = 0) or diverges to infinity (if > 0). Hence
ks(t, 0, 2v1 , 0)k cannot converge to 0 as t tends to infinity and the equilibrium solution cannot be
asymptotically stable.
We have shown that 1 3 2 1. Hence the three statements are equivalent.
Theorem 6.2 does not hold is the system is time varying: We have already seen an example of
a system that is asymptotically but not exponentially stable, so the three statements cannot be
equivalent. How about 3 1? In other words if the eigenvalues of the matrix A(t) Rnn have
negative real parts for all t R is the equilibrium solution of x(t) = A(t)x(t) asymptotically stable?
This would provide an easy test of stability for time varying systems. Unfortunately, however, this
statement is not true.
Example (Unstable time varying system) For a (1, 2) consider the matrix
1 + a cos2 (t) 1 a cos(t) sin(t)
A(t) =
1 + a cos(t) sin(t) 1 + a sin2 (t)
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
94
As the example suggests, determining the stability of time varying linear systems is in general
difficult. Some results exist (for example, if A(t) has negative eigenvalues, is bounded as a function
of time and varies slowly enough) but unfortunately, unlike time invariant linear systems, there
are no simple, general purpose methods that can be used to investigate the stability of time varying
linear systems.
An argument similar to that of Theorem 6.2 can be used to derive stability conditions for time
invariant systems in terms of the eigenvalues of the matrix A.
Theorem 6.3 The equilibrium solution of a linear time invariant systems is stable if and only if
the following two conditions are met:
Similar definitions can be given for the matrix valued functions A(t), B(t), etc. using the corre-
sponding induced norms. Recall also that for existence and uniqueness of solutions we assume that
A(t), B(t), etc. are piecewise continuous functions of time.
1. The equilibrium solution is exponentially stable, i.e. there exists m, > 0 such that for all
(x0 , t0 ) Rn R+ and all t R+ , ks(t, t0 , x0 , U )k mkx0 ke(tt0 ) .
2. For all t0 R, kA()kt0 , , kB()kt0 , , kC()kt0 , , kD()kt0 , are bounded.
Then for all (x0 , t0 ) Rn R and all u() : [t0 , ) Rm with ku()kt0 , bounded,
m
ks(, t0 , x0 , u)kt0 , mkx0 ket0 + kB()kt0 , ku()kt0 ,
hm i
k(, t0 , x0 , u)kt0 , mkC()kt0 , kx0 ket0 + kC()kt0 , kB()kt0 , + kD()kt0 , ku()kt0 , .
If in addition limt u(t) = 0 then limt s(t, t0 , x0 , u) = 0 and limt (t, t0 , x0 , u) = 0.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
95
One can see that also in the case of systems with inputs there is an intimate relation between stability
and continuity of the solution functions in an appropriate function space. To see this consider the
functions
The strange notation is again meant to alert the reader to the fact that we consider s(, t0 , x0 , ) for
fixed t0 , x0 as a function mapping a piecewise continuous function of time (denoted by the placeholder
) to a continuous function of time (denoted by the placeholder , left over after u is substituted for
). Then Theorem ?? directly implies the following.
Corollary 6.1 If the linear time varying system is exponentially stable then the functions s(, t0 , x0 , )
and (, t0 , x0 , ) are continuous.
Example (Unstable with input, asymptotically stable without) For x(t), u(t) R and
t R+ consider the system
1
x(t) = x + u.
1+t
Let x(0) = 0 and apply the constant input u(t) = 1 for all t R+ . Then
Z t Z t
1+ t + t2 /2
s(t, 0, 0, 1) = (t, )u( )d = d =
0 0 1+t 1+t
even though the equilibrium solution is asymptotically stable and the input is bounded since
ku()k = 1 < .
AT P + P A = Q
(known as the Lyapunov equation). To formally state this fact recall that a matrix P Rnn is
called symmetric if and only if P T = P . A symmetric matrix, P = P T Rnn is called positive
definite if and only if for all x Rn with x 6= 0, xT P x > 0; we then write P = P T > 0.
The proof of this fact is deferred to the next section. For the time being we simply show what can
go wrong with the Lyapunov equation if the system is not asymptotically stable.
Example (Lyapunov equation for scalar systems) Let x(t) R and consider the linear system
x(t) = ax(t) for a R; recall that the solution is simply x(t) = eat x(0). For this system the
Lyapunov equation becomes
aT p + pa = 2pa = q
for p, q R. The matrix q is automatically symmetric and it is positive definite if and only if
q > 0. If a < 0 (i.e. the equilibrium solution is asymptotically stable) the Lyapunov equation has
the unique positive definite solution p = q/(2a) > 0. If a = 0 (i.e. the equilibrium solution is stable
but not asymptotically stable) the Lyapunov equation does not have a solution. Finally, if a > 0 (i.e.
the equilibrium solution is unstable) the Lyapunov equation has a unique solution p = q/(2a) < 0
but no positive definite solution.
1. Show that
e(a1)t cos(t) et sin(t)
(t, 0) =
e(a1)t sin(t) et cos(t).
2. Deduce how the stability properties of the system change with the value of a.
3. Find the values of a for which the eigenvalues of A(t) have negative real part for all t 0.
Compare the result with that of 2.
1. Let A Rnm , B Rmp and || || be any induced norm on the corresponding matrix spaces.
Show that ||A B|| ||A|| ||B||.
2. Let p(t) be any polynomial in t. Show that > 0 m > 0 : |p(t)| met , t 0.
1. Let x(0) = 0. Make use of the Laplace transform to derive an input-output (ARMA) model
of the system of the following form: for ai R and Bi Rpm , i = 0, 1, . . . , n,
a0 y (n) (t) + a1 y (n1) (t) + . . . + an y(t) = B0 u(n) (t) + B1 u(n1) (t) + . . . + Bn u(t),
where y (i) (t) (resp. u(i) (t)) denotes the i-th derivative of y (resp. of u) at t. You may assume
that y (i) (0) = 0 and u(i) (0) = 0 for i = 0, . . . n.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
97
A11 0 0
2. Let A = ,B= and C = C1 C2 , where, for 1 n2 n, A22 Rn2 n2 ,
A21 A22 B2
B2 Rn2 m , C2 Rpn2 and all other matrix blocks have consistent dimensions. Compute
the transfer function of the system and use the result to build a state-space representation of
order (i.e. size of the state evolution matrix) n2 having the same transfer function.
M11 0
(Hint: verify that the inverse of an invertible block triangular matrix M = , with
M21 M22
M11 0
M11 and M22 square, is also block triangular, i.e. M 1 = ; in particular, M22
M21 M22
must be equal to... .)
d
1. Consider a system x(t) = f (x(t)). Suppose that we have dt x(t)T P x(t) x(t)T Qx(t),
where P and Q are symmetric positive definite matrices (here the time-derivative is taken
along solutions of the system). Prove that under this condition the system is exponentially
stable, in the sense that its solutions satisfy kx(t)k cet kx(0)k for some c, > 0. Note
that this statement is true whether the system is linear or not. Further, if f (x) = Ax, then
show that the above condition is equivalent to AT P + P A Q.
2. Consider the system x(t) = Ax(t) + Bu(t) such that
eAt
cet for some c, > 0.
(a) Prove that if u is bounded over all time (in the sense that sup0t< ku(t)k M for some
M ), then x is also bounded, for any initial condition.
(b) Now restrict attention to the zero initial condition (x0 = 0). We can view the system above
as a linear operator from the normed linear space of bounded functions u : [0, [ Rm
with norm kuk := sup0t< ku(t)k, to the normed linear space of functions x : [0, [
Rn with the norm kxk := sup0t< kx(t)k. What can you say about the induced norm
of this operator, using the calculations you made in (a)? What can you say about its
continuity?
Problem 6.5 (Discrete-time Systems) Consider a discrete-time linear system x(k+1) = Ax(k)+
Bu(k), k = 0, 1, . . ..
1. Write the formula for the solution x(k) at time k starting from some initial state x(0) at time
0.
2. Assume A is semi-simple. Under what conditions on the eigenvalues of A is the discrete-time
system x(k + 1) = Ax(k) (no controls) asymptotically stable? Stable? Justify your answers.
(Stability definitions are the same as for continuous-time systems, just replace t by k.)
3. Lyapunovs second method for discrete-time system x(k + 1) = f (x(k)) involves the difference
V (x) := V (f (x)) V (x) instead of the derivative V (x); with this substitution, the statement
is the same as in the continuous-time case. Derive the counterpart of the Lyapunov equation
for the LTI discrete-time system x(k + 1) = Ax(k).
Problem 6.6 (Linearization example) Consider a two dimensional state vector x = (x1 , x2 )
R2 whose evolution is goverened by the following differential equations
x1 (t) = x2 (t)
x2 (t) = x1 (t) x1 (t)3 x2 (t).
2. Compute the linearization of the system about its equilibria. The calculation is the same as
the one in Section 4.1, replacing the optimal trajectory, x (t), with the equilibrium solution
x (t) = x for all t, where x is each of the equilibria of the system.
3. Compute the eigenvalues of the matrices of the resulting linearizations and speculate about
the stability of the equilibria of the nonlinear system by invoking Theorem 6.2. Simulate the
system to confirm your intuition.
Problem 6.7 (When local is global) The equilibrium solution x = 0 of the linear time varying
system x(t) = A(t)x(t) is
Is it possible for the linear time varying system to have more than one equilibria? If yes, is it possible
for any of them to be locally asymptotically stable in such a case?
Problem 6.8 (When all is uniform) Show that the equilibrium x = 0 of the linear time invariant
system x(t) = Ax(t) is uniformly stable if and only if it is stable. Show further that this is the case
if and only if the following two conditions are met:
Problem 6.9 (BIBO and BIBS stability) Consider the time varying linear system
1. The equilibrium solution is exponentially stable, i.e. there exists m, > 0 such that for all
(x0 , t0 ) Rn R+ and all t R+ , ks(t, t0 , x0 , U )k mkx0 ke(tt0 ) .
2. For all t0 R, kA()kt0 , , kB()kt0 , , kC()kt0 , , kD()kt0 , are bounded.
Show that:
1. For all (x0 , t0 ) Rn R and all u() : [t0 , ) Rm with ku()kt0 , bounded,
m
ks(, t0 , x0 , u)kt0 , mkx0 ket0 + kB()kt0 , ku()kt0 ,
hm i
k(, t0 , x0 , u)kt0 , mkC()kt0 , kx0 ket0 + kC()kt0 , kB()kt0 , + kD()kt0 , ku()kt0 , .
are continuous.
Chapter 7
We return briefly to abstract vector spaces for a while to introduce the notion of inner products,
that will form the basis of our discussion on controllability.
Definition 7.1 Let (H, F ) be a linear space. A function h, i : H H F is called an inner product
if and only if for all x, y, z H, F ,
Exercise 7.1 For all x, y, z H and all a F , hx + y, zi = hx, zi + hy, zi and hax, yi = ahx, yi.
Moreover, hx, 0i = h0, xi = 0 and hx, xi = 0 if and only if x = 0.
p
Fact 7.1 If (H, F, h, i) then the function k k : H F defined by kxk = hx, xi is a norm on
(H, F ).
Note that the function is well defined by property 3 of the inner product definition. The proof is
based on the following fact.
Theorem 7.1 (Schwarz inequality) With k k defined as in Fact 7.1, |hx, yi| kxk kyk.
Proof: If x = 0 the claim is obvious by Exercise 7.1. If x 6= 0, select F such that || = 1 and
hy,xi
hx, yi = |hx, yi|; i.e. if F = R take to be the sign of hx, yi and if F = C take = |hx,yi| . Then
99
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
100
for all R
kx + yk2 = hx + y, x + yi
= hx + y, xi + hx + y, yi
= hx + y, xi + hx + y, yi
= hx, x + yi + hy, x + yi
= hx, xi + hx, yi) + (hy, xi + hy, yi)
= 2 kxk2 + |hx, yi| + |hx, yi| + ||2 kyk2
= 2 kxk2 + 2|hx, yi| + kyk2 .
This is a quadratic in that must be non-negative for all R. This will be the case if and only if
it is non-negative at its minimum point. Differentiating with respect to shows that the minimum
occurs when
|hx, yi|
=
kxk2
Substituting this back into the quadratic we see that
Definition 7.2 Let (H, F, h, i) be an inner product space. The norm defined by kxk2 = hx, xi
for all x H is called the norm defined by the inner product. If the normed space (H, F, k k) is
complete (a Banach space) then (H, F, h, i) is called a Hilbert space.
To demonstrate the above definitions we consider two examples of inner products, one for finite
dimensional and one for infinite dimensional spaces.
Example (Finite dimensional inner product space) For F = R or F = C, consider the linear
space (F n , F ). Define the inner product h, i : F n F n F by
n
X
hx, yi = xi yi = xT y
i=1
Exercise 7.3 Show that this satisfies the axioms of the inner product.
is the Euclidean norm. We have already seen that (F n , F, k k) is a complete normed space, hence
(F n , F, h, i) is a Hilbert space.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
101
Exercise 7.4 Show that this satisfies the axioms of the inner product.
As before, to make this norm well defined we need to identify functions that are zero almost every-
where (i.e. are zero over a dense set) with the zero function.
Example (Orthogonal vectors) Consider the inner product space (Rn , R, h, i) with the inner
project hx, yi = xT y defined above. Given two non-zero vectors x, y Rn one can define the angle,
, between them by
hx, yi
= cos1
kxk kyk
x, y are orthogonal if and only if = (2k + 1)/2 for k Z.
Exercise 7.5 Show that is well defined by the Schwarz inequality). Show further that =
(2k + 1)/2 if and only if x and y are orthogonal.
This interpretation of the angle between two vectors motivates the following generalisation of a
well known theorem in geometry.
Theorem 7.2 (Pythagoras theorem) Let (H, F, h, i) be an inner product space. If x, y H are
orthogonal then kx + yk2 = kxk2 + kyk2 , where k k is the norm defined by the inner product.
Proof: Exercise.
Note that H does not need to be R3 (or even finite dimensional) for the Pythagoras theorem to hold.
Given a subspace M H of a linear space, consider now the set of all vectors y H which are
orthogonal to all vectors x M .
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
102
Definition 7.4 The orthogonal complement of a subspace, M , of an inner product space (H, F, h, i)
is the set
M = {y H | hx, yi = 0 x M }.
Consider now an inner product space (H, F, h, i) with the norm kyk2 = hy, yi defined by the
inner product. Recall that a sequence {yi } i=0 is said to converge to a point y H if an only if
limi ky yi k = 0. In this case y is called the limit point of the sequence {yi }
i=0 . Recall also that
a subset K H is called closed (Definition 3.3) if and only if it contains the limit points of all the
sequences {yi }
i=0 K.
Fact 7.2 Let M be a subspace of the inner product space (H, F, h, i). M is a closed subspace of
H and M M = {0}.
hx, yi = hx, yi + y yi i
= hx, yi i + hx, y yi i
= hx, y yi i.
Definition 7.5 Let M, N be subspaces of a linear space (H, F ). The sum of M and N is the set
If in addition M N = {0} then M + N is called the direct sum of M and N and is denoted by
M N.
Fact 7.3 V = M N if and only if for all x V there exists unique u M and v N such that
x = u + v.
Proof: (exercise) (). Let V = M N . Then M N = {0} and for all x V there exist u M
and v N such that x = u + v. It remains to show that the u and v are unique. Assume, for the
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
103
sake of contradiction, that they are not. Then there exist u M , v N with (u , v ) 6= (u, v) such
that x = u + v . Then
u + v = u + v u u = v v .
Moreover, u u M and v v N (since M, N are subspaces). Hence M N {u u, 0} =
6 {0}
which is a contradiction.
(). Assume that for all x V there exist unique u M and v N such that x = u + v.
Then V = M + N by definition. We need to show that M N = {0}. Assume, for the sake of
contradiction that this is not the case, i.e. there exists y 6= 0 such that y M N . Consider and
arbitrary x V and the unique u M and v N such that x = u + v. Define u = u + y and
v = v y. Note that u M and v N since M and N are subspaces and y M N . Moreover,
u + v = u + y + v y = u + v = x but u 6= u and v 6= v since y 6= 0. This contradicts the
uniqueness of u and v.
Theorem 7.3 Let M be a closed subspace of a Hilbert space (H, F, h, i). Then:
1. H = M M (denoted by M M ).
2. For all x H there exists a unique y M such that x y M . This y is called the
orthogonal projection of x onto M .
3. For all x H the orthogonal projection y M is the unique element of M that achieves the
minimum
kx yk = inf{kx uk | u M }.
for all u U , v V .
Therefore the adjoint of the linear map defined by the matrix A F mn is the linear map defined by
T
the matrix A = [aji ] F nm , the complex conjugate transpose (known as the Hermitian transpose)
of A. If in addition F = R, then A = AT = [aji ] F nm is simply the transpose of A.
As we have seen in Chapter 2, any linear map between finite dimensional vector spaces can be
represented by a matrix, by fixing bases for the domain and co-domain spaces. The above example
therefore demonstrates that when it comes to linear maps between finite dimensional spaces the
adjoint operation always involves taking the complex conjugate transpose of a matrix. Recall also
that linear maps between finite dimensional spaces are always continuous (Theorem 3.5). For infinite
dimensional spaces the situation is in general more complicated. In some cases, however, the adjoint
may be simple to represent.
Example (Infinite dimensional adjoint) Let U = (L2 ([t0 , t1 ], F m ), F, h, i2 denote the space of
square integrable functions u : [t0 , t1 ] F m and V = (F n , F, h, iF n . Consider G() L2 ([t0 , t1 ], F nm )
(think for example of G(t) = (t1 , t)B(t) for t [t0 , t1 ] for a linear system). Define a function
A : U V by Z t1
A(u()) = G(t)u(t)dt.
t0
It is easy to see that A is well defined and linear. Assume also that G(t) is well behaved enoug (e.g.
it is bounded) so that that A is also continuous.
For arbitrary x V , u() U
Z t1
T
hx, A(u())iF m = x G(t)u(t)dt
t0
Z t1
= xT G(t)u(t)dt
t0
Z t1 T
T
= G(t) x u(t)dt
t0
T
= hG x, ui2 .
Definition 7.7 Let (H, F, h, i) be a Hilbert space and A : H H be linear and continuous. A is
called self-adjoint if and only if A = A, in other words for all x, y H
Example (Finite dimensional self-adjoint map) Let H = F n with the standard inner product
and A = [aij ] F nn . The linear map A : F n F n defined by A(x) = Ax is self adjoint if and
T
only if A = A, or in other words aij = aji for all i, j = 1, . . . n. Such matrices are called Hermitian.
If in addition F = R, A is self-adjoint if and only if aij = aji , i.e. AT = A. Such matrices are called
symmetric.
Example (Infinite dimensional self-adjoint map) Let H = L2 ([t0 , t1 ], R) and K(, ) : [t0 , t1 ]
[t0 , t1 ] R such that
Z t1 Z t1
|K(t, )|2 dtd <
t0 t0
Let (H, F ) be a linear space and A : H H a linear map. Recall that (Definition 2.14) an element
F is called an eigenvalue of A if and only if there exists v H with v 6= 0 such that A(v) = v;
in this case v is called an eigenvector corresponding to .
Fact 7.4 Let (H, F, h, i) be a Hilbert space and A : H H be linear, continuous and self-adjoint.
Then
Range(A) = {v V | u U, v = A(u)}
Null(A) = {u U | A(u) = 0}
Theorem 7.5 (Finite Rank Lemma) Let F = R or F = C, let (H, F, h, i) be a Hilbert space and
recall that (F m , F, h, iF m ) is a finite dimensional Hilbert space. Let A : H F m be a continuous
linear map and A : F m H be its adjoint. Then:
A|Range(A ) : Range(A ) F m
A |Range(A) : Range(A) H
Proof: Part 1: A is linear and continuous and hence A A and A A are linear, continuous and
self-adjoint by Theorem 7.4.
Part 2: Recall that Range(A) F m and therefore is finite dimensional. Moreover, Dim[Range(A )]
Dim[F m ] = m therefore Range(A ) is also finite dimensional. Therefore both Range(A) and
Range(A ) are closed (by Fact ??), hence by Theorem 7.3
But
Therefore, Null(A ) = Range(A) therefore F m = Range(A) Null(A ). The proof for H is
similar.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
107
Part 3: Consider the restriction of the linear map A to the range space of A
A|Range(A ) : Range(A ) F m .
A x = 0 A A (x) = 0 x Null(A A )
Range(A) = {y F m | x H, y = A(x)}
= {y F m | x Range(A ), y = A(x)} (by Part 3)
= {y F m | u F m , y = A(A (u))}
= {y F m | u F m , y = (A A )(u)}
= Range(A A ).
The proofs of Part 5 and Part 6 are analogous to those of Part 3 and Part 4 respectively.
Notice that the assumption that the co-domain of A is finite dimensional is only used to establish
that the ranges of A and A are closed sets, to be able to invoke Theorem 7.3. In fact, Theorem 7.5
holds more generally if we replace the spaces by their closures.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
108
Prove that A is linear and continuous. Also prove that A is self-adjoint if K is symmetric.
Chapter 8
Can the input u be used to steer the state of the system to an arbitrary value?
Can we infer the value of the state by observing the input and output?
We will again start by a general discussion of these questions for non-linear systems. We will then
specialize to the case of time varying linear systems and then further to the case of time invariant
linear systems.
of system (8.1)(8.2). For simplicity, assume further that r is continuous in x and piecewise contin-
uous in t.
Definition 8.1 Consider t0 , t1 R with t0 t1 . The input trajectory u() PC([t0 , t1 ], Rm ) steers
(x0 , t0 ) Rn R to (x1 , t1 ) Rn R if and only if s(t1 , t0 , x0 , u) = x1 . The system (8.1)(8.2) is
controllable on [t0 , t1 ] if and only if for all x0 , x1 Rn there exists u() PC([t0 , t1 ], Rm ) that steers
(x0 , t0 ) to (x1 , t1 ).
109
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
110
FS
Notice that controllability has nothing to do with the output of the system, it is purely an input
to state relationship; we will therefore talk about controllability of the system (8.1) ignoring
Equation (8.2) for the time being. The following fact is an immediate consequence of the definition.
Fact 8.1 The system (8.1)(8.2) is controllable on [t0 , t1 ] if and only if for all x0 Rn the function
s(t1 , to , x0 , ) : PC([t0 , t1 ], Rm ) Rn is surjective.
It is easy to see that controllability is preserved under state and output feedback. Consider a state
feedback map FS : Rn R Rm and a new input variable v(t) Rm . Let u(t) = v(t) + FS (x(t), t)
and define the state feedback system
x(t) = f (x(t), w(t) + FO (y(t), t), t) = f (x(t), w(t) + FO (r(x(t), t), t), t) = fO (x(t), w(t), t) (8.4)
Proof: The proof requires only careful application of the definitions. For example, to show that if
the system of Equation (8.1) is controllable on [t0 , t1 ] then so is the system of Equation (8.3), consider
arbitrary x0 , x1 Rn and look for v() PC([t0 , t1 ], Rm ) that steers (x0 , t0 ) to (x1 , t1 ). Since (8.1) is
controllable on [t0 , t1 ] there exists u() PC([t0 , t1 ], Rm ) such that the unique solution, s(t, t0 , x0 , u),
of (8.1) satisfies s(t1 , t0 , x0 , u) = x1 . Define v(t) = u(t) FS (s(t, t0 , x0 , u), t); clearly v() : [t0 , t1 ]
Rm is piecewise continuous since u is piecewise continuous, FS is assumed to be continuous in its
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
111
replacemen
w(t) u(t)
+ x(t) R x(t) y(t)
f r
FO
first argument and piecewise continuous in the second, and, by definition, s(, t0 , x0 , u) is continuous.
Moreover, under v() to system (8.3) evolves according to
x(t) = f (x(t), v(t) + FS (x(t), t), t) = f (x(t), u(t) FS (s(t, t0 , x0 , u), t) + FS (x(t), t), t)
starting at x(t0 ) = x0 . Clearly x(t) = s(t, t0 , x0 , u) satisfies both the initial condition and the
differential equation and therefore is also the unique solution for system (8.3) under the input v().
Hence, the proposed input v() steers (x0 , t0 ) to (x1 , t1 ).
The remaining parts are similar and are left as an exercise.
Fact 8.2 The system (8.1)(8.2) is observable on [t0 , t1 ] if and only if for all u() PC([t0 , t1 ], Rm )
the function
(, t0 , , u) : Rn PC([t0 , t1 ], Rp )
x0 7 (, t0 , x0 , u) : [t0 , t1 ] Rp
is injective.
It is easy to see that if we establish the value of x0 Rn we can in fact reconstruct the value of x(t)
for all t [t0 , t1 ]; this is because, by uniqueness, x(t) = s(t, t0 , x0 , u) is uniquely determined once we
know t, t0 , x0 and u : [t0 , t] Rm .
It is easy to see that observability is preserved under output feedback and input feed-forward.
Consider a input feed-forward map FF : Rm R Rp and a new output variable z(t) Rp . Let
z(t) = y(t) + FF (u(t), t) and define the input feed-forward system with
FF
Proof: Exercise.
Notice that observability is not preserved by state feedback. If system (8.1)(8.2) is observable,
system (8.3)(8.2) may or may not be observable and vice versa.
Exercise 8.1 Provide a simple example of a linear system whose observability properties are altered
by state feedback.
In nonlinear control, feedback linearization provides a dramatic demonstration of the fact that
state feedback can alter observability. In feedback linearization state feedback and a coordinate
transformation is used to alter the dynamics of some of the states, so that system behaves like a
linear system. The price to pay for this is that the remaining states (the so-called zero dynamics)
are disconnected from the output, i.e. are rendered unobservable; the reader is referred to [?, ?, ?]
for more information on this topic.
2 m
for u, u L ([t0 , t1 ], R ). Note that by considering this extended space we can no longer assume
that the input trajectories will be piecewise continuous, and hence are no longer covered by the
existence uniqueness argument of Theorem 3.7. We will alleviate this difficulty, however, by showing
that the properties of interest here (e.g. steering the system from a given initial state to a desired
final state) can indeed be achieved by piecewise continuous trajectories.
Since controllability is an input to state property only the matrices A() and B() come into play.
To simplify the notation we give the following definition.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
113
Definition 8.3 The pair (A(), B()) is controllable on [t0 , t1 ] if and only if for all x0 , x1 Rn there
exists u : [t0 , t1 ] Rm that steers (x0 , t0 ) to (x1 , t1 ), i.e.
Z t1
x1 = (t1 , t0 )x0 + (t1 , )B( )u( )d.
t0
Exercise 8.2 The three statements of Theorem 8.3 are not equivalent for nonlinear systems. Where
does linearity come into play in the proof?
Theorem 8.3 states that for linear systems controllability is equivalent to controllability to zero and
to reachability from zero. In the subsequent discussion will use the last of the three equivalent
statements to simplify the analysis; Theorem 8.3 implies that this can be done without loss of
generality. In line with this, we introduce the following definition.
Definition 8.4 A state x1 Rn is reachable on [t0 , t1 ] by the pair (A(), B()) if and only if there
exists u() L2 ([t0 , t1 ], Rm ) that steers (0, t0 ) to (x1 , t1 ). The reachability map on [t0 , t1 ] of the pair
(A(), B()) is the function
Lr : L2 ([t0 , t1 ], Rm ) Rn
R t1
u 7 t0 (t1 , )B( )u( )d.
Fact 8.3 Lr is linear and continuous. The set of reachable states is the linear subspace Range(Lr ).
Proof: The fact that Lr is linear is obvious from the properties of the integral. The fact that
reachable states are equal to Range(Lr ) is immediate from the definition of reachable states. To
show that Lr is continuous let kk denote the 2-norm in Rn and note that for all u() L2 ([t0 , t1 ], Rm )
Z t1
2
2
kLr (u)k =
(t1 , )B( )u( )d
t0
Z t1
2
k(t1 , )B( )k ku( )k2 d.
t0
Since (t1 , )B() is a piecewise continuous function on [t0 , t1 ], there exists some M > 0 such that
k(t1 , )B( )k M for all [t0 , t1 ]. Hence
Z t1
kLr (u)k2 M 2 ku( )k2 d = M 2 ku()k22
t0
Definition 8.5 The controllability gramian of the pair (A(), B()) on [t0 , t1 ] is the matrix
Z t1
Wr (t0 , t1 ) = (t1 , )B( )B( )T (t1 , )T d Rnn .
t0
The controllability gramian will be our primary vehicle for studying the controllability of time
varying linear systems. We start by establishing some of its basic properties.
Fact 8.4 The controllability gramian, Wr (t0 , t1 ), has the following properties:
1. It is symmetric.
2. It is positive semi-definite.
3. For all t0 t0 , Wr (t0 , t1 ) Wr (t0 , t1 ) in the sense that xT [Wr (t0 , t1 ) Wr (t0 , t1 )]x 0 for
all x Rn .
The following theorem allows us to answer controllability questions by simply checking the rank of
the controllability gramian.
Proof: 1 2: By Theorem 8.3, (A(), B()) is controllable on [t0 , t1 ] if and only if all states are
reachable from 0 on [t0 , t1 ], i.e. if and only if Lr is surjective, i.e. if and only if Range(Lr ) = Rn .
2 3: By the finite rank lemma.
3 4: Notice that Lr Lr : Rn Rn is a linear map between two finite dimensional spaces.
Therefore it admits a matrix representation. We will show that Wr (t0 , t1 ) Rnn is the matrix
representation of this linear map. Then Lr Lr is surjective if and only if Wr (t0 , t1 ) is invertible i.e.
if and only if Det[Wr (t0 , t1 )] 6= 0.
Recall that Z t1
Lr : u 7 (t1 , )B( )u( )d
t0
and hLr x, ui = hx, Lr ui. We have already seen that (Section 7.3)
Therefore Z
t1
[Lr Lr ](x) = (t1 , )B( )B( )T (t1 , )T d x = Wr (t0 , t1 )x.
t0
Exercise 8.3 Show that (A(), B()) is controllable on [t0 , t1 ] if an only if Wr (t0 , t1 ) is positive
definite (i.e. xT Wr (t0 , t1 )x > 0 for all x Rn with x 6= 0).
Lr (u) = x1 .
Lr (u + u) = Lr (u) + Lr (u) = x1 + 0 = x1 .
Therefore for any input u that steers (0, t0 ) to (x1 , t1 ) any other input of the form u + u with
u Null(Lr ) will do the same. So in general there will be an infinite number of inputs that get
the job done. The question now becomes can we somehow find a canonical one? It turns out that
this can be done by a projection onto Range(Lr ).
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
116
L2 ([t0 , t1 ], Rp ) Rn
Range(Lr ) Lr Range(Lr )
u
x1 Lr Lr
u
Lr
Null(Lr ) = {0}
u
Null(Lr )
is a bijection. Therefore, for all x1 Range(Lr ) there exists a unique u Range(Lr ) such that
Lr (u) = x1 . We have seen that for any u Null(Lr ) if we let u = u + u then Lr (u) = x1 . However,
But u Range(Lr ), u Null(Lr ) and by the finite rank lemma Range(Lr ) is orthogonal to
Null(Lr ). Therefore hu, ui2 = hu, ui2 = 0 and
since kuk22 0. Therefore, among all the u L2 ([t0 , t1 ], Rm ) that steer (0, t0 ) to (x1 , t1 ), the unique
u Range(Lr ) is the one with the minimum 2norm. The relation of the different spaces involved
is depicted in Figure 8.4.
We now return to the general controllability case and formally state our findings.
Theorem 8.5 Assume (A(), B()) are controllable in [t0 , t1 ]. Given x0 , x1 Rn define u : [t0 , t1 ]
Rm by u = Lr (Lr Lr )1 [x1 (t1 , t0 )x0 ], i.e.
Then
Proof:
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
117
Part 1:
Z t1
x(t1 ) = (t1 , t0 )x0 + (t1 , t)B(t)u(t)dt
t0
Z t1
= (t1 , t0 )x0 + (t1 , t)B(t)B(t)T (t1 , t)T Wr (t0 , t1 )1 [x1 (t1 , t0 )x0 ] dt
t0
Z t1
= (t1 , t0 )x0 + (t1 , t)B(t)B(t)T (t1 , t)T dt Wr (t0 , t1 )1 [x1 (t1 , t0 )x0 ]
t0
= (t1 , t0 )x0 + Wr (t0 , t1 )Wr (t0 , t1 )1 [x1 (t1 , t0 )x0 ]
= x1 .
Part 2 is obvious from the formula, since B(t) is piecewise continuous, (t1 , t) is continuous and
all other quantities are constant. In fact the only discontinuity points are those of B(t). Hence if
B(t) is continuous (e.g. for time invariant linear systems) then the minimum energy controls are
continuous.
Part 3:
Z t1
kuk22 = u(t)T u(t)dt
t0
Z t1
= [x1 (t1 , t0 )x0 ]T (Wr (t0 , t1 )1 )T (t1 , t)B(t)B(t)T (t1 , t)T dt
t0
Wr (t0 , t1 )1 [x1 (t1 , t0 )x0 ]
T
= [x1 (t1 , t0 )x0 ] (Wr (t0 , t1 )T )1 Wr (t0 , t1 )Wr (t0 , t1 )1 [x1 (t1 , t0 )x0 ]
= [x1 (t1 , t0 )x0 ]T Wr (t0 , t1 )1 [x1 (t1 , t0 )x0 ]
The perceptive reader will have noticed that this calculation is effectively the same as the one carried
out in Section 7.5 to compute the pseudo-inverse of a matrix. The only difference is that here one
of the two spaces is infinite dimensional. Compare also Figure 8.4 with Figure ??.
The theorem provides several interesting insights into the structure of these so-called minimum
energy controls. For one, note that even though for mathematical reasons we were forced to allow
more general control inputs u L2 ([t0 , t1 ], Rm ) to apply the Finite Rank Lemma, the minimum
energy controls (the unique u Range(Lr ) that steers (x0 , t0 ) to (x1 , t1 )) happens to be piecewise
continuous. In other words, if there is a control in u L2 ([t0 , t1 ], Rm ) that will get the job done
there is also one in PC([t0 , t1 ], Rm ) (in fact, a better one!) that will do the same. This is good
to know, since we do not have to worry about existence of solutions which we only established for
piecewise continuous input functions in Chapter 3.
The structure of the minimum energy itself, kuk22 , is also revealing. kuk22 is a quadratic function of
x1 (t1 , t0 )x0 , a term which reflects the distance between where we want the system to be at time
t1 (namely x1 ) and where it would end up if we left it alone (namely (t1 , t0 )x0 ). In other words, the
further we want to push the system from its natural course the more energy we need to expend. The
term in the middle, Wr (t0 , t1 )1 , in a sense reflects how controllable the system is. By Theorem 8.4,
if the system is uncontrollable Wr (t0 , t1 ) will be singular and, loosely speaking, Wr (t0 , t1 )1 would
be infinite, suggesting that we may not be able to push the system where we want however hard
we try. If the system is weakly controllable Wr (t0 , t1 ) will be positive definite but small, hence
Wr (t0 , t1 )1 will be large. In other words, the less controllable a system is the more energy we
need to expend to get it where we want. The observation that Wr (t0 , t1 ) Wr (t0 , t1 ) for t0 t0 and
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
118
t1 t1 (Fact 8.4) also admits a natural interpretation: The more time we have to take the system
where we want the less energy we need to do so (i.e., the more relaxed we can be about it!)
The proof of the theorem reveals that an even more general statement is possible. Even if the system
is not controllable, if it happens that x1 (t1 , t0 )x0 Range(Lr ) then we can still find a minimum
energy input that steers (x0 , t0 ) to (x1 , t1 ). The construction will be somewhat different (in this case
Wr (t0 , t1 ) will not be invertible) but the general idea is the same: Consider the restricted map
Lr |Range(L ) : Range(Lr ) Range(Lr )
r
which is bijective even if Range(Lr ) 6= Rn and select the unique u Range(Lr for which Lr (u) =
x1 (t1 , t0 )x0 . Indeed, even if x1 (t1 , t0 )x0 6 Range(Lr ) we can still formulate a more
general problem where the objective is to minimize kx1 x(t1 )k using a minimum energy input.
The construction will be even more complicated in this case, but the idea is again similar to the
pseudo-inverse calculations: Project x1 (t1 , t0 )x0 onto Range(Lr ) and among the u that map
to the projection select the unique one that belongs to Range(Lr ).
Finally, we note that Theorem 8.5 provides a simple example of an optimal control problem. The u
of the theorem is the input that drives the system from (x0 , t0 ) to (x1 , t1 ) and minimizes
Z t1
2
kuk2 = u(t)T u(t)dt.
t0
This is the starting point for more general optimal control problems where costs of the form
Z t1
xT (t)Q(t)x(t) + u(t)T R(t)u(t) dt
t0
nn
for some Q(t) R and R(t) Rmm symmetric and positive definite are minimized, subject
to the constraints imposed by the dynamics of the system and the desired initial and final states.
An example are the so-called Linear Quadratic Regulator (LQR) problems. The interested reader
is referred to [?, ?, ?] for more information on this topic.
Definition 8.6 The pair of matrices (C(), A()) is called observable on [t0 , t1 ] if and only if for
all x0 Rm and all u : [t0 , t1 ] Rm , one can uniquely determine x0 from the information
{(u(t), y(t)) | t [t0 , t1 ]}.
Note that once we know x0 and u : [t0 , t] Rm we can reconstruct x(t) for all t [t0 , t1 ] by
Z t
x(t) = (t, t0 )x0 + (t, )B( )u( )d.
t0
the last two terms in the summation can be reconstructed if we know u : [t0 , t1 ] Rm . There-
fore the difficult part in establishing observability is to determine x0 from the zero input response
C(t)(t, t0 )x0 with t [t0 , t1 ]. Without loss of generality we will therefore restrict our attention to
the case u(t) = 0 for all t [t0 , t1 ].
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
119
Definition 8.7 A state x0 Rn is called unobservable on [t0 , t1 ] if and only if C(t)(t, t0 )x0 = 0
for all t [t0 , t1 ].
Clearly the state x0 = 0 is unobservable. The question then becomes whether there are any x0 6= 0
that are also unobservable. This would be bad, since the zero input response from these states would
be indistinguishable from that of the 0 state.
Definition 8.8 The observability map of the pair (C(), A()) is the function
Lo : Rn L2 ([t0 , t1 ], Rp )
x0 7 C(t)(t, t0 )x0 t [t0 , t1 ].
Fact 8.5 Lo is a linear, continuous function of x0 . Moreover, for all x0 Rn , Lo (x0 ) PC([t0 , t1 ], Rp .
The state x0 Rn is unobservable if and only if x0 Null(Lo ). The pair of matrices (C(), A())
is observable if and only if Null(Lo ) = {0}.
The proof is left as an exercise. Notice that to be able to apply the Finite Rank Lemma we are forced
to consider functions y() : [t0 , t1 ] Rp that belong to the space L2 ([t0 , t1 ], Rp ) in Definition 8.8.
Fact 8.5, however, clarifies that this is nothing to worry about: The only functions we will encounter
will live in Range(Lo ) and will be piecewise continuous, in fact the only discontinuity points will
be those of C(t).
Definition 8.9 The observability gramian of the pair (C(), A()) is the matrix
Z t1
Wo (t0 , t1 ) = (, t0 )T C( )T C( )(, t0 )d Rnn .
t0
The following theorem provides a complete characterization of observability of linear time varying
systems in terms of the properties of the matrix Wo (t0 , t1 ).
Proof: The proof is effectively the same as that of Theorem 8.4 with Lo in place of Lr . Indeed,
Wo (t0 , t1 ) turns out to be the representation of the linear map Lo Lo : Rn Rn in the basis used
for the representation of the matrix A().
The similarity between the controllability and observability theorems is further highlighted by the
decomposition into subspaces induced by Lo , as shown in Figure 8.5. Comparing Figure 8.5 to
Figure 8.4, it becomes apparent that controllability and observability are in fact dual concepts, they
are just the two faces of the Finite Rank Lemma. To formalize this statement, note that given the
system
Rn L2 ([t0 , t1 ], Rm )
Range(Lo ) Lo Range(Lo )
y
Lo Lo x0
y
Lo
Null(Lo ) = {0}
y
Null(Lo )
with x(t) Rn , u(t) Rm , y(t) Rp , one can define the dual system
x(t) = A(t)T x(t) C(t)T u(t) (8.8)
T T
y(t) = B(t) x(t) + D(t) u(t) (8.9)
with x(t) Rn , u(t) Rp and y(t) Rm . The two systems are closely related in terms of their
controllability and observability properties.
Theorem 8.7 (Duality theorem) Let (t, t0 ), (t, t0 ) Rnn denote the state transition matri-
ces of system (8.6)(8.7) and system (8.8)(8.9) respectively.
Proof: For Part 1, recall that by Theorem 4.2, (t0 , t0 ) = (t0 , t0 )T = I and t (t, t0 ) =
A(t)T (t, t0 ). Moreover
(t0 , t)(t, t0 ) = (t0 , t0 ) [(t0 , t)(t, t0 )] = [(t0 , t)] (t, t0 ) + (t0 , t) (t, t0 ) = 0
t t t
[(t0 , t)] (t, t0 ) = (t0 , t)A(t)(t, t0 )
t
(t0 , t)T = A(t)T (t0 , t)T A(t).
t
The conclusion follows by existence and uniqueness.
For Part 2, by Theorem 8.6 system (8.8)(8.9) is observable on [t0 , t1 ] if and only if the matrix
R t1 T
t0
(, t0 )T B( )T B( )T (, t0 )d is invertible. Substituting (t, t0 ) from Part 1 this matrix
is the same as
Z t1 Z t1
T T
(t0 , )B( )B( ) (t0 , ) d = (t0 , t1 )) (t1 , )B( )B( )T (t1 , )T d (t0 , t1 )T .
t0 t0
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
121
Rn L2 ([t0 , t1 ], Rp ) L2 ([t0 , t1 ], Rm ) Rn
Lo Lr
Lo Lo Range(L
o) Range(Lo ) Range(L
r) Range(Lr ) Lr Lr
Lo Lr
Null(L
o) Null(Lr )
Observability Controllability
Since (t0 , t1 ) is invertible by Theorem 4.2, the above is invertible if and only if
Z t1
(t1 , )B( )B( )T (t1 , )T d
t0
is invertible, i.e. if and only if system (8.6)(8.7) is controllable (by Theorem 8.4).
The proof of Part 3 is similar to that of Part 2 and is left as an exercise.
Theorem 8.8 Assume that the pair of matrices (C(), A()) is observable and consider an arbitrary
y L2 ([t0 , t1 ], Rp ). Then
Z t1
x0 = (Lo Lo )1 Lo (y) = [Wo (t0 , t1 )]1 (, t0 )T C( )T y( )d
t0
Proof: By the Finite Rank Lemma, y uniquely decomposes into a sum y = y+y with y Range(Lo )
and y Null(Lo ) (see Figure 8.5). For all x Rn , we clearly have Lo (x) Range(Lo ) and hence
y Lo (x) Range(Lo ) (recall that Range(Lo ) is a subspace of L2 ([t0 , t1 ], Rp )). Since Range(Lo )
and Null(Lo ) are orthogonal
Lo (x) : Rn Range(Lo )
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
122
is bijective . Therefore, there exists a unique x0 Rn such that Lo (x0 ) = y. Note that
Writing the above in the coordinates used in the representation of A() and recalling the infinite
dimensional adjoint calculation in Section 7.3 shows that
Z t1
1
x0 = [Wo (t0 , t1 )] (, t0 )T C( )T y( )d
t0
Notice that the theorem does not require the observed output, y(), to lie in Range(Lo ), i.e. to
correspond exactly to the zero input response of some initial state x0 Rn . In fact, the observed
y() does even have to be piecewise continuous, as one would expect from a zero input response. The
formula provided in the theorem generates the x0 whose zero input res Ponce most closely matches
the observed y() in the 2-norm sense. This is very useful, since in practice the measured system
response will be corrupted by noise and other disturbances, hence is unlikely to correspond exactly
to some initial state.
As in the case of controllability, the above calculation can also be extended to the case where the
system is not observable. The minimizer will not be unique in this case and the computation will
be more involved since W0 (t0 , t1 ) will not be invertible, but the idea is similar: Find the projection
of the measured output y onto Range(Lo ) and the minimum norm x0 that corresponds to it.
Proof: Consider any [t0 , t1 ] and without loss of generality assume that t0 = 0.
Part 1: x0 is unobservable if and only if x0 Null(Lo ), i.e. if and only if
Part 2:
Part 3: Consider x0 Null(O), i.e. CAk x0 = 0 for all k = 0, . . . , n 1. We would like to show
that x = Ax0 Null(O). Indeed Cx = CAx0 = 0, CAx = CA2 x0 = 0, . . . , until
CAn1 x = CAn x0 = C(n I + n1 A + . . . + 1 An1 )x0 = 0
s(t, 0, v, U ) = eAt e = et v
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
124
CeAt e = et Cv = 0, t 0
and {w1 , . . . , wnr , e1 , . . . , er } form a basis for Rn . The representation of x Rn with respect to
this basis decomposes into two parts, x = (x1 , x2 ), with x1 Rnr (observable) and x2 Rr
(unobservable). For the rest of this proof we assume that all vectors are represented with respect
to this basis.
Recall that, from Part 2, Null(O) is an Ainvariant subspace, i.e. for all x Null(O), Ax
Null(O). Note that
0
x Null(O) x = .
x2
Moreover,
A11 A12 0 A12 x2
Ax = =
A21 A22 x2 A22 x2
Therefore, since Ax Null(O) we must have that A12 x2 = 0 for all x2 Rr , which implies that
A12 = 0 R(nr)r (think of using the elements of the canonical basis of Rr as x2 one after the
other). Notice further that if x Null(O) then in particular x Null(C) and
0
Cx = C1 C2 = C2 x2 = 0 x2 Rr .
x2
Hence C2 = 0 Rpr .
In summary, in the new coordinates the system representation is
Take C an eigenvalue of A22 Rrr (and note that will also be an eigenvalue of A Rnn ).
Let v Cr with v 6= 0 be the corresponding eigenvector. Then
I A 0 0
= .
C v 0
The theorem provides two easy ways of checking the observability of a time invariant linear system,
by checking the rank of matrices. Notice that the conditions of the theorem are independent of
1 This argument is fine as long as is real. If is complex, strictly speaking, we have not constructed a real vector
x0 Rn which is unobservable. This can be done, however, by taking linear combinations of the complex eigenvectors,
as in the proof of Theorem 6.2. The construction is left as an exercise.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
125
[t0 , t1 ]. Therefore, (C, A) is observable on some [t0 , t1 ] with t0 < t1 if and only if it is observable
on all [t0 , t1 ]. To put it another way, observations {(u(t), y(t) | t [0, ]} over an arbitrary small
interval are sufficient to reconstruct the initial condition x0 and hence all future states of the system
(at least in theory).
The most commonly used test for observability involves testing the rank of the matrix O. Though
easy to implement for small systems, this test may face numerical difficulties for larger systems, since
it requires the computation of powers of A up to An1 , which may lead to numerical conditioning
problems. Testing the rank of
I A
C
tends to be more robust numerically on the other hand. Notice that even though the theorem
suggests that this condition needs to be tested for all C (clearly an impossible task) the proof
indicates that the condition need only be checked for Spec[A], which is feasible.
The decomposition of the state developed in the proof of Part 4 also shows that for any (C, A),
the state of the system decomposes into an observable and and unobservable part. An immediate
corollary of the proof of the theorem is that the observable part is indeed observable.
Corollary 8.1 For any matrices C Rpn and A Rnn there exists a change of basis T Rnn
with Det[T ] 6= 0, such that in the new coordinates the representation of the matrices decomposes
into
A11 0
A = T AT 1 = and C = CT 1 = C1 0
A21 A22
and the pair of matrices (C1 , A11 ) is observable.
Proof: The decomposition was already established in the proof of Theorem 8.9. To show that
(C1 , A11 ) are observable assume, for the sake of contradiction, that they are not. Then there exists
x1 R(nr)r with x1 6= 0 such that
C1
C1 A11
.. x1 = 0.
.
nr1
C1 A11
The structure of the matrices implies that
C1 0
C1 A11 0
x1 x1
O = .. .. = 0.
0 . . 0
nr1
C1 A11 0
Hence (x1 , 0) Null(O). By the choice of basis, however, (x1 , 0) Span{w1 , . . . , wnr } =
Null(O) . Since (x1 , 0) 6= 0 this is a contradiction (by Fact 7.2).
to infinity without any indication of this at the output. Indeed, if the system is unobservable and one
of the eigenvalues of the matrix A22 above has positive real part the state x2 (t) may go to infinity
as t increases. Since, however, (0, x2 ) is in the nullspace of O, the output y(t) will be unaffected by
this (either directly or through the state x1 ) and may remain bounded. For unobservable systems
one would at least hope that the unobservable modes (the eigenvalues of the matrix A22 ) are stable.
In this case, if the output of the system is bounded one can be sure that the state is also bounded.
This requirement, which is weaker than observability, is known as detectability.
and use it to test whether the systems is controllable. The main facts in this direction are summarized
in the following theorem.
Rank [I A B] = n.
The proof is effectively the same as that of Theorem 8.9 by duality. As before, the condition of
Part 4 need only be tested for Spec[A]. Moreover, for all (A, B) the state decomposes into a
controllable and an uncontrollable part. For an appropriate choice of basis the system representation
becomes:
x1 (t) = A11 x1 (t) +A12 x2 (t) +B1 u(t)
x2 (t) = +A22 x2 (t)
y(t) = C1 x1 (t) C2 x2 (t) +Du(t)
where the pair of matrices (A11 , B1 ) is controllable. Likewise, the spectrum of A decomposes into
where the former set contains eigenvalues whose eigenvectors are reachable (controllable modes) while
the latter contains eigenvalues whose eigenvectors are unreachable (uncontrollable modes).
Dually to observability, an immediate danger with uncontrollable systems is that the state may be
diverging to infinity without the input being able to prevent it. If the system is uncontrollable and
one of the eigenvalues of the matrix A22 above has positive real part the state x2 (t) may go to
infinity as t increases. Since the evolution of this state is not affected by the input (neither directly
nor through the state x1 ) there is nothing the input can do to prevent this. For uncontrollable
systems one would at least hope that the uncontrollable modes (the eigenvalues of the matrix A22 )
are stable. This requirement, which is weaker than controllability, is known as stabilizability.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
127
Theorem 8.11 (Kalman decomposition) For an appropriate change of basis for Rn the state
vector is partitioned into
x1 controllable, observable
x2 controllable, unobservable
x=
x3 .
uncontrollable, observable
x4 uncontrollable, unobservable
with
A11 0 B1 A11 A13
, controllable and C1 C3 , observable.
A21 A22 B2 0 A33
where the first set contains controllable and observable modes, the second set controllable and unob-
servable modes, the third set uncontrollable and observable modes and the fourth set uncontrollable
and unobservable modes.
Pictotially the interdependencies of the sub-systems revealed by the Kalman decomposition are
shown in Figure 8.7. Note that this figure is not a block diagram, the arrows represent dependencies
between subsystems rather than signals. The labels next to the arrows indicate which of the matrices
in Theorem 8.11 is responsible for the corresponding dependence.
u(t) B1 C1 y(t)
A11 +
A21 A13
B2 A23 C3
A22 A33
A33 A43
A44
Wr (t1 , t0 ) Wr (t1 , t0 ) 0 t0 t0 . ()
Problem 8.2 (Controllability) Consider a harmonic oscillator with control input u, satisfying
the equation
x1 0 1 x1 0
= + u
x2 1 0 x2 1
Suppose we want to drive the system from the state [1, 0]T to the state [0, 0]T in 2 units of time.
Do there exist constants u1 , u2 , u3 such that we can make the transfer from [1, 0]T at t = 0 to
[0, 0]T at time t = 2?
3. Suppose that the system is initially at rest, i.e. x(0)=0, and for x Rn we wish to find a
piecewise continuous control u : [0, +) Rm such that x(t) = x for all t 1. Such a control
can be found for all x Rn .
Problem 8.4 (Observability) For some matrices A Rnn , B Rn1 , C R1n and D R,
consider a single-input, single-output time-invariant system
x = A(t) x
e(t) e(t) C(t) u
e(t)
ye(t) = B(t) xe(t) + D(t) ue(t)
Knowing whether a system is controllable or observable is usually not an end in itself. One wants
to subsequently exploit this knowledge to design controllers that will make the system behave in
some desired way: Ensure stability, force the output to track some desired trajectory, force the state
to converge to an equilibrium at some desired rate, etc. We have already seen how for controllable
systems one can generate input trajectories to steer the system from its initial state to some desired
final state (Theorem 8.5). In practice, however, this may not be sufficient. Due to small deviations
in the initial state or in the matrices involved in the dynamics, these so-called open loop input
trajectories will invariably fail to steer the system as expected. Indeed, if the system is unstable
the deviation may get arbitrarily large. A much more robust way of steering the system is to
measure its state as we go along and if we find that it deviates from the desired trajectory adapt
the open loop inputs to correct the deviation. The simplest case of such a feedback control involves
stabilization, where the desired trajectory is simply an equilibrium solution and the objective is to
design a controller to make the system asymptotically stable. Similar (dual) comments can also be
made for the merits of closed loop state reconstruction by incorporating output measurements
on-line, as opposed to the batch processing approach of Theorem 8.8.
In this chapter we investigate how such controllers and state estimators can be designed. We first
consider the so-called state feedback controllers, where one assumes that the entire state of the
system can be measured and used when making decisions. We then extend the approach to output
feedback, where one assumes that only the outputs of the system are available for measurement.
In particular, we consider the case of observer based output feedback, where an observer is used to
reconstruct the value of the state from the output measurements; a state feedback controller then
uses the reconstructed state to steer the system.
The discussion is restricted to linear time invariant systems. In this case, one makes use of the
controllability and observability matrices derived in Chapter 8 to transform the system into special
forms for which controllers and observers can easily be designed. The transformations involve changes
of basis using invertible matrices derived from the controllability and observability matrices. For
clarity of exposition we first recall some facts about basis changes.
130
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
131
with A Rnn , B Rnm , C Rpn and D Rpm . Throughout these notes several facts
regarding the relation between these two descriptions of the system have been established. We
summarize the most important ones in the following theorem.
Theorem 9.1 Consider the linear time invariant systems (9.1)(9.2) and (9.3)(9.4) related through
the change of coordinates x(t) = T x(t) for all t R+ , where T Rnn is an invertible matrix. The
following hold:
4. The two systems have the same impulse response matrix, i.e.
Exercise 9.1 Are the state transition and impulse state transition matrices of the two systems the
same?
In this chapter we will repeatedly use the properties of coordinate transformations to bring the
system equations into special forms that make it easier to design controllers and observers. We will
then exploit the properties listed in Theorem 9.1 to transfer the design to the original coordinates.
The two forms we will especially be interested in are the so-called controllable canonical form and
observable canonical form. Both of these will first be presented for the easier case of single input,
single output systems and then generalized.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
132
where u(t) R and B Rn is a vector. The output equations will play no role in this section and
will be omitted for simplicity, we will return to them in Section 9.3.
Recall that the system (9.5) is controllable if and only if the matrix
P = B AB . . . An1 B Rnn (9.6)
has full rank, or, equivalently (since the matrix is square), if and only if it is invertible, or, equiva-
lently, if and only if its columns are linearly independent. We will use the last fact to construct an
invertible matrix T Rnn to use as a change of coordinates. We will start by defining a matrix
S Rnn that will be invertible whenever the system is controllable and will eventually become
T 1 .
Let the characteristic polynomial of the matrix A Rnn be
Det[I A] = n + 1 n1 + . . . + n1 + n .
sn = B
sn1 = Asn + 1 B = AB + 1 B
sn2 = Asn1 + 2 B = A2 B + 1 AB + 2 B
..
.
s1 = As2 + n1 B = An1 B + 1 An2 B + . . . + n1 B.
Lemma 9.1 As1 + n B = 0. Moreover, the matrix S = [s1 . . . sn ] Rnn is invertible if and only
if the system (9.5) is controllable.
Proof: By definition
As1 + n B = An B + 1 An1 B + . . . + n1 AB + n B
= An + 1 An1 + . . . + n1 A + n I B = 0
which is not equal to zero if and only if the system is controllable; notice that the determinant of
the matrix on the right is email to 1.
The lemma shows that whenever system (9.5) is controllable we can use S to define a change of
coordinates. It turns out that in the new coordinates the system matrices take a particularly useful
form.
Theorem 9.2 The system (9.5) is controllable if and only if there exists a change of coordinates
x(t) = T x(t) with T Rnn invertible, such that the matrices A = T AT 1 and B = T B satisfy:
0 1 0 ... 0 0
0 0 1 ... 0 0
.. .
.. .
.. . .. .
.. , B = .
A = . .. (9.7)
0 0 0 ... 1 0
n n1 n2 . . . 1 1
Proof: (): Assume that the system is controllable. Define S Rnn as above and let T = S 1 .
Note that
0 0
.. ..
B = S . B = T B = S 1 B = . .
0 0
1 1
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
134
Consider next the matrix AS = A [s1 . . . sn ]. Reading the columns from the right we have
0 0
.. ..
. .
Asn = sn1 1 B = S 0 , Asn1 = sn2 2 B = S 1 , etc. until
1 0
1 2
0
..
.
As1 = n B = S 0 (by Lemma 9.1).
0
n
Hence,
0 1 0 ... 0
0 0 1 ... 0
.. .. .. .. ..
AS = S . . . . .
0 0 0 ... 1
n n1 n2 . . . 1
which implies that
0 1 0 ... 0
0 0 1 ... 0
.. .. .. .. ..
. . . . . = S 1 AS = T AT 1 = A.
0 0 0 ... 1
n n1 n2 . . . 1
(): Assume that there exists a change of coordinates bringing the systems in the form (9.7). Note
that in the new coordinates
0 0 ... 0 1
... 1
h i 0 0
. . .. ..
P = B . . . An1 B = .. .. . .
0 1 ...
1 1 . . .
h i
where stands for some number which depends on 1 , . . . , n . Hence Det P = 1 and the system
in the new coordinates is controllable. Theorem 9.1 then implies that the system in the original
coordinates is also controllable.
The system representation of equation (9.7) is known as the controllable canonical form. The term
canonical refers to the fact that any controllable system can be brought to this form through a
change of coordinates, as Theorem 9.2 shows. The controllable canonical form is useful for designing
state feedback controllers to modify the behavior of system (9.5). A state feedback controller is a
method for selecting the inputs of the system as a function of the state to force the system to exhibit
some desired behavior. Here we will restrict our attention to linear time invariant state feedback
controllers of the form
u(t) = Kx(t) + r(t), (9.8)
where K Rmn is known as the gain matrix, and r(t) Rm is an external input vector; in the
single input case considered in this section K = [k1 . . . kn ] R1n is just a row vector.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
135
The system obtained by connecting a state feedback controller of the form (9.8) with a system of
the form (9.5) is known as the closed loop system (Figure 9.1). Its dynamics are described by the
linear state equations
x(t) = (A + BK)x(t) + Br(t). (9.9)
We would like to select the gain matrix K such that the closed loop systems exhibits some desired
behavior, for example is asymptotically stable and its state tends to zero sufficiently fast. To ensure
properties like these one typically needs to select K so that the eigenvalues of the closed loop system
matrix A + BK coincide with some desired set {1 , . . . , n } C. Clearly for this to be possible the
set must contain the complex conjugates of all its members, i.e. {1 , . . . , n } = {1 , . . . , n }; we
call such a set of complex numbers a complex conjugate set.
Theorem 9.3 System (9.5) is controllable if and only if for all complex conjugate sets {1 , . . . , n }
C there exists K Rmn such that Spec[A + BK] = {1 , . . . , n }.
Proof: (). Assume that (9.5) is controllable. Given a desired complex conjugate set {1 , . . . , n }
C we would like to select feedback gains K such that the closed loop matrix A + BK has this set as
eigenvalues. In other words, we would like
Det[I (A + BK)] = ( 1 )( 2 ) . . . ( n )
= n + d1 n1 + . . . + dn1 + dn (9.10)
Notice that the system is still in controllable canonical form. Therefore the terms in the last row
are the coefficients of the characteristic polynomial
If we set
k1 = 1 d1 , . . . , kn = n dn
this becomes identical to the desired characteristic polynomial (9.10). Hence with this selection for
the matrix K we have made the eigenvalues of A+ B K equal to the desired set {1 , . . . , n }. Finally,
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
136
since eigenvalues are not affected by coordinate transformations (Theorem 9.1) setting K = KT
ensures that with the state feedback u(t) = KT x(t) + r(t) the eigenvalues of the matrix A + BK of
the closed loop system are also equal to the desired list.
(). By contraposition. Assume that the system is not controllable. Then, by Theorem 8.10, there
exists a change of coordinates x(t) = T x(t) such that
A11 A12 B1
x(t) = x(t) + u(t).
0 A22 0
Notice that part of the spectrum (the eigenvalues of A22 ) is not affected by the feedback gains K.
Therefore given a list of desired eigenvalues {1 , . . . , n } C it is impossible, in general, to select K
such that Spec[A + B K] = {1 , . . . , n }; this is only possible if Spec[A22 ] {1 , . . . , n }. Setting
K = KT and noting that Spec[A + BK] = Spec[A + B K] (Theorem 9.1) completes the proof.
Theorem 9.3 provides a method for moving the eigenvalues of a controllable linear system to arbitrary,
complex conjugate values by selecting a feedback gain matrix K and applying linear state feedback
u(t) = Kx(t) + r(t). This method is known in the literature as pole placement, eigenvalue placement,
or eigenvalue assignment. Even though for the proof of the theorem we first have to bring the system
into controllable canonical form, in practice this is not necessary. Given matrices A Rnn and
B Rnm we can directly compute the characteristic polynomial Det[I (A + BK)], treating the
entries in the matrix K as variables. The coefficients of the characteristic polynomial will turn out
to be linear functions of the entries of the gain matrix K Rmn . Equating these coefficients with
those of the desired characteristic polynomial (9.10) leads to a system of linear equations where the
entries of K play the role of the unknowns. Theorem 9.3 then ensures that this system will have
a solution which, if used in the feedback function u(t) = Kx(t) + r(t), will lead to a closed loop
system with the desired eigenvalues. As we will see below, this procedure works more generally,
even if m > 1. The difference is that for single input systems (m = 1) the system of linear equations
has the same number of equations as unknowns (both n). Thus (assuming controllability) there
is a unique solution to this system of linear equations and hence a unique feedback gain matrix
K R1n that will get the job done. For multi-input systems (m > 1), on the other hand, there
are more unknowns (nm) than equations (n). Therefore there will be multiple choices of K that
will lead to the same eigenvalues for the close loop system; in fact a whole n(m 1) dimensional
subspace of them! Unless other considerations are present (minimizing some cost criterion, etc.) one
can eliminate the redundant degrees of freedom by setting some of the elements of K equal to zero.
This leads to a sparser feedback gain matrix, which both simplifies the calculations and makes the
feedback function easier to implement in practice.
Theorem 9.3 also shows that the procedure outlined above for designing feedback gain matrices will
in general fail if the system is not controllable. The same steps can again be executed, but in this case
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
137
System (A, B, C, D)
Output feedback
the system of linear equations in the entries of K will not have any solution in general. The reason
is the presence of the uncontrollable modes (the elements of Spec[A22 ] in the proof of Theorem 9.3)
which are not affected by our choice of K. If Spec[A22 ] contains some eigenvalues with positive real
part this is usually bad news: Whatever gain matrix K we select, it will not be possible to obtain a
stable closed loop system. If stability is what we are after, then drastic system re-design is needed.
We may, for example, try to add, or modify one of the system actuators to augment the matrix B
and make the offending modes controllable. If, on the other hand, Spec[A22 ] lies entirely on the left
half of the complex plane (i.e. the system is stabilizable) it is still possible to stabilize the system.
We just have to resign ourselves to the fact that the eigenvalues in Spec[A22 ] will stay put whatever
we do and select K to move the remaining eigenvalues (the controllable modes in Spec[A11 ]) to
desired locations. To do this, we first we have to determine the uncontrollable modes (e.g. using the
rank test in Theorem 8.10, condition 4), then construct a list of desired eigenvalues containing the
uncontrollable modes, compute the characteristic polynomial Det[I (A + BK)], and equate its
coefficients with those of the desired characteristic polynomial. The proof of Theorem 9.3 together
with the fact that the uncontrollable modes are included in the list of desired eigenvalues will ensure
that the resulting system of linear equations will have a solution for the entries of K. The solution
will not be unique, however, even for single input systems. Some of the degrees of freedom (those
corresponding to K2 in the proof of Theorem 9.3) will be redundant. Unless other considerations
are present, one can again select these redundant degrees of freedom to make the matrix K as sparse
as possible.
for the input to apply. An algorithm that generates an estimate of the state using the past and
present inputs and outputs is known as a state estimator or a state observer. The resulting feedback
arrangement will look like Figure 9.3.
It turns out that for linear systems it suffices to consider linear state observers. Analogously to
the case of state feedback design we will start by describing the design of such observers for single
output systems (p = 1) and return to the general case of multi-output systems (p > 1) toward the
end of this chapter. We will consider the usual linear system equations
x(t) = Ax(t) + Bu(t) (9.12)
y(t) = Cx(t) + Du(t) (9.13)
where C R1n and D R1m are just row vectors.
A linear state observer generates an estimate of the state, x(t), at time t R+ using the values
of the input and output for all times up to t, {(u( ), y( )) | [0, t]}. The observer can itself be
thought of as a linear system, with n states (x(t) Rn ), m + p inputs ((u(t), y(t)) Rm+p ), and
n outputs (the state estimate itself x(t) Rn ). The equations of the observer try to mimic those
of the linear system (9.12)(9.13), with a correction term to account for any mismatch between the
measured output, y(t), at time t and the output predicted by the observer (denoted by y(t)) based
on its own state estimate, x(t). More precisely,
x(t) = Ax(t) + Bu(t) + L(y(t) y(t)) (9.14)
y(t) = C x(t) + Du(t). (9.15)
Notice that at time t the observer only makes use of its own estimate of the state, x(t), the input
u(t) and the output y(t) to update its estimate of the state; in particular the unknown value of the
state, x(t), is not used anywhere. The linear correction term, L(y(t) y(t)), introduces an indirect
estimate of the difference between the estimated and the true value of the state which, as we will
see below, can be used to make the two converge to each other. The matrix L Rnp is known
as the observer gain matrix. Note also that only the present values of y(t) and u(t) are used in
Equations (9.14)(9.15). The state estimate x(t) at time t, however, will depend on all of the past
values of the inputs and outputs, {(u( ), y( )) | [0, t]} (as well as the initial condition x0 that we
select). In a sense, x(t) encapsulates all the relevant information contained in the input and output
trajectories up to time t.
Under what conditions will x(t) be a good estimate of x(t)? To study this question we can look at the
estimation error, e(t) = x(t) x(t) Rn . If this error (or more precisely its norm) is small then the
state estimate generated by the observer is accurate. Clearly since initially we know nothing about
the true value of x(0) our initial estimate, x(0), and hence the initial error, e(0), will be arbitrary.
The question is whether the estimate will get better as time goes on and more information is collected
through u(t) and y(t). To answer this question we can look at how the error, e(t), evolves over time:
e(t) = x(t) x(t)
= Ax(t) + Bu(t) Ax(t) Bu(t) L(y(t) y(t))
= A(x(t) x(t)) L(Cx(t) + Du(t) Cx(t) Du(t))
= (A LC)(x(t) x(t))
= (A LC)e(t).
Notice that the evolution of the estimation error is also governed by a linear system without any
inputs. The estimation error will therefore converge to zero if and only if the eigenvalues of the
matrix A LC have negative real part. Otherwise, as time goes on our estimate of the state will
remain as bad as our initial guess, or get even worse, despite the fact that we have collected more
information through the output of the system.
The question now becomes can we select the observer gain matrix L so that the eigenvalues of ALC
have negative real parts? It is easy to see that this question is the dual of the question addressed
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
139
x(t)
K Observer
Figure 9.3: Linear feedback with state estimation. x(t) denotes the estimate of the state
in Section 9.2. If the system is observable, one can in fact mimic the pole placement method step
by step to develop a method for placing the eigenvalues of the estimation error dynamics in some
desired, complex conjugate values. First, one can write the system in observable canonical form
using a change of coordinates.
Theorem 9.4 The system (9.12)(9.13) is observable if and only if there exists a change of coor-
dinates x(t) = T x(t) with T Rnn invertible, such that the matrices A = T AT 1 and C = BT 1
satisfy:
0 0 ... 0 n
1 0 . . . 0 n1
A = . . .. . . .. ,
.. .. . . .
0 0 ... 1 1
C = 0 0 ... 0 1 .
The proof is the dual to that of Theorem 9.2: One uses the rows of the observability matrix
C
CA
O= .. Rnn
.
CAn1
which under the observability assumption are linearly independent to generate a change of coordi-
nates and bring the system to the desired form. One can then use the observable canonical form to
show the following.
Theorem 9.5 System (9.12)(9.13) is observable if and only if for all complex conjugate sets
{1 , . . . , n } C there exists L Rmn such that Spec[A LC] = {1 , . . . , n }.
The procedure for selecting the observer gain matrix is also similar. One uses the desired eigenvalues
{1 , . . . , n } to generate a desired characteristic polynomial. Then one computes the polynomial
Det[I (A LC)] and equates its coefficients to those of the desired polynomial, giving rise to a
system of n linear equations in np unknowns (the entries of the observer gain matrix L Rnp ). If
the system (9.12)(9.13) is observable, then the system of linear equations has a unique solution (if
p = 1) or multiple solutions (if p > 1), for which the eigenvalues of A LC are equal to the desired
ones. If the system (9.12)(9.13) is not observable, on the other hand, the system of linear equation
will in general have no solution. If the system is not even detectable (i.e. some of the unobservable
modes have non-negative real parts) then there is no hope of designing a stable observer that will
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
140
generate a correct estimate of the state, even asymptotically. If the system is detectable one can
still design a stable observer by computing the unobservable modes, including them in the list of
desired eigenvalues, and placing the remaining eigenvalues in some desired locations as above.
The closed loop system is therefore itself a linear time invariant system, with 2n states (x(t), x(t))
R2n ), m inputs (r(t) Rm ) and p outputs (y(t) Rp ) and state equations
x(t) A BK x(t) B
= + r(t)
x(t) LC A + BK LC x(t) B
x(t)
y(t) = C DK + Dr(t).
x(t)
To determine whether this system is stable it is easier to start with a coordinate transformation,
making the state estimation error one of the states. Note that
x(t) x(t) I 0 x(t) x(t)
= = =T
e(t) x(t) x(t) I I x(t) x(t)
In summary, the 2n eigenvalues of the closed loop system coincide with the n eigenvalues of the
system with perfect state feedback and the n eigenvalues of the state estimation error. This fact
greatly simplifies the design of output feedback controllers for linear systems. One can design the
observer and the state feedback gain matrix separately, put them together, and the resulting closed
loop system will work as planned. This separation of the controller from the observer design, known
as the principle of separation, is a fundamental property of linear systems. Unfortunately, even
though related principles have been established for other classes of systems, this nice property does
not always hold.
From the above discussion it appears that, assuming that the system is controllable and observable,
the eigenvalues of the closed loop system can be arbitrarily placed. Therefore in principle one can
make the closed loop system respond arbitrarily fast, by making the real part of all its eigenvalues
sufficiently negative. There are several good practical reasons, however, why this temptation should
be resisted. The most important is modeling inaccuracies. Any model of a system used for design
or analysis is an approximate mathematical abstraction of a physical process; in fact a linear model
is usually a rather crude approximation of reality, valid only for certain values of the state and
input.The real physical process is bound to exhibit nonlinearities, additional dynamics, noise in the
sensors and actuators, delays, and other phenomena not adequately captured by the model. Faster
response, small eigenvalues and high gains tend to make the closed loop system more sensitive to all
these unmodeled factors.
In the end, the choice of the eigenvalues of the closed loop system (and hence the gains of the
controller and observer) usually comes down to a trade-off between several such considerations.
For linear systems, methods that in some cases allow one to establish optimal trade-offs have been
developed, in areas such as optimal filtering, linear quadratic Gaussian control, robust control, etc.
Very often, however, insight, intuition, and trial-and-error play a central role.
nn
Consider again system
AR
(9.5) with and B Rnm , possibly with m > 1. Consider the
columns of B = b1 . . . bm and assume that they are linearly independent. This assumption
can be made without loss of generality. Indeed, if the columns of B are linearly dependent there
exists a change of coordinates in Rm u(t) = T u(t) with T Rmn invertible such that for B = B T 1
and an appropriate partition of u(t):
u1 (t)
Bu(t) = B u = B1 0
u2 (t)
and the columns of B1 are linearly independent; 0 denotes a matrix of appropriate dimensions all
of whose elements are equal to zero. Hence in the new coordinates the inputs u2 (t) are redundant
and we can work with just the reduced system
has rank equal to n. As for the single input case, the idea is to use this matrix to generate a change of
coordinates, x(t) = T x(t) for some T Rnn invertible, so that the resulting matrices A = T AT 1
and B = T B are such that one can easily design feedback controllers for them. Since the matrix
has rank n, out of its nm columns,
we can select n linear independent ones. The way we do the selection is not particularly important.
Different choices will lead to somewhat different structures in the matrices A and B, hence the
resulting form we obtain in the end will be somewhat less canonical than the controllable canonical
form is for single input systems. The most common choice is to start selecting the columns from
the left, skipping any columns that are linearly dependent on the previously selected ones, until
n linearly independent columns have been selected. Since we have assumed that the columns of
B are linearly independent, b1 , . . ., bm will be among the selected columns. In addition, for each
j = 1, . . . , m there is a maximal integer kj = 1, . . . , n such that Akj 1 bj will be among the selected
columns. Or, equivalently, a minimal integer such that Akj bj can be written as a linear combination
of the previously selected columns.
Exercise 9.3 Show that if Akj 1 bj is among the selected columns for some kj then so is Akj i bj
for all i = 1, . . . , kj . Or, equivalently, that if Akj bj can be written as a linear combination of the
previously selected columns then so can Ai bj for all i kj .
b1 , . . . , Ak1 1 b1
b2 , . . . , Ak2 1 b2
..
.
bm , . . . , Akm 1 bm .
We will assume that the controllability indices come in increasing order k1 k2 . . . km ; if not,
reorder the columns of B and P below so that this is the case.
Arrange the selected columns in a new matrix
P = b1 . . . Ak1 1 b1 b2 . . . Ak2 1 b2 . . . bm . . . Akm 1 bm Rnn .
Note that P is invertible since, by construction, its columns are linearly independent. Compute the
inverse of P and consider its rows
p1
p2
P 1 = . Rnn , p1 , . . . , pn R1n .
..
pn
Among these rows select the ones corresponding to the controllability indices pk1 , pk1 +k2 , . . . ,
pk1 +...+km = pn and form the matrix
pk1
pk1 A
..
.
pk1 Ak1 1
..
T = . R(k1 +...+km )n = Rnn .
pn
p A
n
.
..
km 1
pn A
Lemma 9.2 pk1 Ak1 1 b1 = pk1 +k2 Ak2 1 b2 = . . . = pn Akm 1 bm = 1. Moreover, pk1 Ai b1 = 0 for
i = 0, . . . , k1 2 and pk1 Ai bj = 0 for j = 2, . . . , m and i = 0, . . . , kj 1. Likewise, pk1 +k2 Ai b2 = 0
for i = 0, . . . , k2 2 and pk1 +k2 Ai bj = 0 for j = 1, 3, . . . , m and i = 0, . . . , kj , etc. The matrix T is
invertible.
Proof: By definition
p1
P 1 P = ... b1 . . . Ak1 1 b1 . . . bm . . . Akm 1 bm
pn
p1 b1 . . . p1 Ak1 1 b1 . . . p1 Akm 1 bm
.. .. .. .. ..
. . . . .
= pk1 b1 . . . pk1 Ak1 1
b1 . . . pk1 A km 1
bm
= I.
. . . . .
.. .. .. .. ..
k1 1 km 1
pn b 1 . . . pn A b 1 . . . pn A bm
Therefore, the elements on the diagonal are all equal to 1; in particular,
pk1 Ak1 1 b1 = pk1 +k2 Ak2 1 b2 = . . . = pn Akm 1 bm = 1.
Likewise, the off-diagonal elements are all equal to 0; in particular
pk1 b1 = . . . = pk1 Ak1 2 b1 = pk1 b2 = . . . = pk1 Ak2 1 b2 = . . . = pk1 bm = . . . = pk1 Akm 1 bm = 0
..
.
pn b1 = . . . = pn Ak1 1 b1 = pn b2 = . . . = pn Ak2 1 b2 = . . . = pn bm = . . . = pn Akm 2 bm = 0.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
144
In the last row we have decomposed the n n matrix into blocks Pij Rki kj . By the first part of
the lemma, every block of the form Pii , for i = 1, . . . , m has 1 on its anti-diagonal (elements (1, ki ),
(2, ki 1), . . ., (ki , 1)) and 0 everywhere else. Every block of the form Pij for i 6= j, on the other hand
has all its elements equal to zero. Therefore the last matrix is invertible (indeed, its determinant is
equal to either +1 or to 1). Since Det[T P ] = Det[T ]Det[P ] and P is invertible by construction
(recall that its columns were selected to be linearly independent), we must have that Det[T ] 6= 0.
Hence the matrix T is invertible.
Consider now the coordinate transformation x(t) = T x(t). It is easier to think of the matrices
A = T AT 1 and B = T B as decomposed into blocks
A11 . . . A1m B1
.. and B = .. ,
A = ... ..
. . .
Am1 . . . Amm Bm
Scanning the elements of B and comparing to Lemma 9.2 it becomes apparent that all of them are
equal to 0, except one element in each column, which is equal to 1, namely
The first row of A11 will have its second entry equal to 1 (the coefficient of pk1 A), and all other
entries equal to 0 (the coefficients of pk1 , pk1 A2 etc.).
The second row of A11 will have its third entry equal to 1 (the coefficient of pk1 A2 ) and all
other entries equal to 0, etc. until,
row k1 1 of A11 , which will have its last entry equal to 1 and all other entries equal to 0.
Likewise,
the first row of A12 will have all its entries equal to 0 (coefficients of pk1 +k2 , . . ., pk1 +k2 Ak2 1 ),
etc. until,
Row k1 1 of A1m , which will have all its entries equal to 0.
Row k1 is equal to pk1 Ak1 and (since the rows of T are linearly independent and span R1n ) can be
written as a linear combination
pk1 Ak1 = ak1 1 pk1 + ak1 2 pk1 A + . . . ak1 n pn Akm 1
for some ak1 1 , . . ., ak1 n R; we will not bother too much with the exact form of these coefficients.
Repeating the process for the remaining blocks shows that in the new coordinates the system matrices
A and B become respectively
0 1 ... 0 0 0 ... 0 0 ... 0
.. .. .. .. .. .. .. .. .. . . .
. . . . ... . . . . .
. ..
0 0 ... 1 0 0 ... 0 0 ... 0
ak1 1 ak1 2 . . . ak1 k1 ak1 (nkm +1) ak1 (nkm +2) . . . ak1 n
1 ... 0
.. .. . . .
. . .. , .. . . . .. .
0 0 ... 0 0 1 ... 0 0 ... 0
.. .. .. .. .. .. .. .. . . ..
. . . . . . . . . . . .
... .
0 0 ... 0 0 0 ... 1 0 ... 0
an1 an2 . . . ank1 an(nkm +1) an(nkm +2) . . . ann 0 ... 1
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
146
This special form is known as the Brunowski normal form and it the analogue of the controllable
canonical form for multi-input systems. We note, however, that, unlike the controllable canonical
form, the Brunowski normal form is not unique, since it depends on the choice of columns of the
matrix P , the way these columns are rearranged, etc.
How does the Brunowski normal form help us to design controllers for the system? Assume that we
are given a complex conjugate set of desired eigenvalues {1 , . . . , n }. As in the single input case
we would like to select a linear state feedback gain matrix K Rmn so that under the feedback
the closed loop matrix A+ BK has the desired eigenvalues. To see that this is possible first bring the
system in Brunowski normal form through a coordinate transformation x(t) = T x(t) and determine
the controllability indices k1 , . . ., km . Write also the state feedback in the new coordinates
Split the set of desired eigenvalues into m subsets {1 , . . . , k1 }, . . ., {nkm +1 , . . . , n }. For sim-
plicity assume that this can be done so that each of the subsets is itself complex conjugate. Form
the m characteristic polynomials
( 1 ) . . . ( k1 ) = k1 + d11 k1 1 + . . . + d1k1
...
( nkm +1 ) . . . ( n ) = km + dm1 km 1 + . . . + dmkm .
The idea is to use each row Ki of the feedback matrix to ensure that the subsystem Aii has the
corresponding characteristic polynomial. To do this select
k11 = (ak1 1 + d1k1 ), . . . , k1k1 = (ak1 k1 + d11 ), . . . , k1(k1 +1) = ak1 (k1 +1) , . . . , k1n = ak1 n
...
km1 = an1 , . . . , km(nkm ) = an(nkm ) , . . . , kmn = (ak1 n + dm1 ).
The elements of the off-diagonal blocks Aij for i 6= j are then eliminated and the resulting closed
loop system matrix becomes
0 1 ... 0 0 0 ... 0
.. .. . . .. .. .. .. .
. . . . ... . . . ..
0 0 . . . 1 0 0 ... 0
d1k1 d1(k1 1) . . . d11 0 0 ... 0
.. .. ..
A + B K = . . . .
0 0 ... 0 0 1 ... 0
.
.. .. . . . ...
. ..
.
..
.
..
.
..
.
...
0 0 ... 0 0 0 ... 1
0 0 ... 0 dmkm dm(km 1) . . . dm1
Notice that the matrix is block-diagonal and that each of the diagonal blocks is in controllable
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
147
which is the desired characteristic polynomial. Setting K = KT completes the design (recall that
coordinate transformations do not affect the eigenvalues by Theorem 9.1).
As one can imagine from the above discussion, bringing a system into Brunowski normal form can be
rather tedious. Fortunately this is not necessary if the end goal is to design feedback controllers for
the system. As for single input systems, to do this it suffices to form the characteristic polynomial
of the matrix (A + BK), treating the elements of the feedback gain matrix K Rmn as unknowns.
Equating the coefficients of this polynomial to those of the desired characteristic polynomial (
1 ) . . . ( n ) leads to a system of n linear equations with the nm elements of K as unknowns.
The Brunowski normal form construction guarantees that if the system is controllable this system
of equations will have a solution; for multi-input systems it will in fact have an infinite number of
solutions, since the number of equations is smaller than the number of unknowns. If the system
is not controllable but merely stabilizable then one has to make sure that the uncontrollable (but
stable) modes are included in the set of desired eigenvalues and repeat the process.
In a similar way one can also construct observers and implement output feedback controllers for the
system. The construction is just dual to that of the Brunowski normal form and will not be given
in detail.
Compute the observer gain L such that the dynamics of the error e := x x have three poles
at s = 3.
Problem 9.2 (Controllable canonical form) Use the procedure of Section 9.5 to re-derive the
controllable canonical form for single input systems. Explain how the coefficients in the last row of
matrix A arise.
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
148
Problem 9.3 (Controllability, observability and state feedback) Consider the linear state
feedback arrangement of Figure 9.1.
1. Show that the closed loop system is controllable if and only if the open loop system is control-
lable.
2. Assume now that the open loop system is observable. Is it true that the closed loop system
will always be observable? Provide a proof, or a counter-example.
Problem 9.4 (Controllability, observability and observer feedback) Consider the observer
feedback arrangement of Figure 9.3. Assume that the open loop system is observable and control-
lable. Will the closed loop system always be observable? Will it be controllable? In both cases
provide either a proof, or a counter-example.
Appendix A
Coverage:
149
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
150
No nonlinear systems
No discrete event or hybrid systems
No infinite dimensional systems (no PDEs, occasional mention of delays)
Logistics:
B.1 Shorthands
Def. = Definition
Thm. = Theorem
Ex./ = Example or exercise
iff = if and only if
wrt = with respect to
wlog = without loss of generality
ftsoc = for the sake of contradiction
= therefore
= contradiction
: = such that
= Q.E.D. (quod erat demonstrandum)
B.2 Sets
, 6, , 6, , ,
For A, B X, Ac stands for the complement of a set and \ for set difference, i.e. Ac = X \ A
and A \ B = A B c .
R real numbers, R+ non-negative real numbers, [a, b], [a, b), (a, b], (a, b) intervals
Q rational numbers
Z integers, N non negative integers (natural numbers)
C complex numbers, C+ complex number with non-negative real part
Sets usually defined through their properties as in: For a < b real
151
c J. Lygeros & F. Ramponi, 2013
Lecture Notes on Linear System Theory,
152
Cartesian products of sets: X, Y two set, X Y is the set of ordered pairs (x, y) such that
x X and y Y .
B.3 Logic
, , , , .
! = exists unique.
= and
= or
= not
Exercise B.1 Is the statement (! x R : x2 = 1) true or false? What about the statement
( x R : x2 = 1)?
Appendix C
153
Bibliography
154