TABLE OF CONTENTS
1 Introduction and Generalities 1
1.1 Basic concepts from linear functional analysis . . . . . . . . . . . . . . . . . 2
1.1.1 Linear operators . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2 Matrix calculus and basic Operator theory . . . . . . . . . . . . . . . . . . . 5
1.2.1 Matrices, eigenvalues and eigenvectors . . . . . . . . . . . . . . . . . 5
1.2.2 Limits of sequences of operators . . . . . . . . . . . . . . . . . . . . . 8
1.3 Review of calculus . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.3.1 The mean value theorem . . . . . . . . . . . . . . . . . . . . . . . . . 12
1.3.2 Integration in Banach space . . . . . . . . . . . . . . . . . . . . . . . 12
2 Basic notions of ordinary dierential equations 13
2.1 Generalities . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.2 Existence and Uniqueness of Solutions to a System . . . . . . . . . . . . . . 15
2.2.1 General theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.2.2 Linear systems of ODE . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.3 Stability theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.3.1 Phase space . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
2.3.2 General denition of stability . . . . . . . . . . . . . . . . . . . . . . 26
2.3.3 Stability of linear systems . . . . . . . . . . . . . . . . . . . . . . . . 28
3 Floquet theory: Presentation and stability of periodic solutions 32
3.1 Linear systems with periodic coecients: Floquet theory . . . . . . . . . . . 32
3.1.1 Nonhomogeneous linear systems . . . . . . . . . . . . . . . . . . . . . 36
3.2 Stability of periodic solutions . . . . . . . . . . . . . . . . . . . . . . . . . . 37
3.2.1 Autonomous systems . . . . . . . . . . . . . . . . . . . . . . . . . . . 38
3.2.2 Nonautonomous systems . . . . . . . . . . . . . . . . . . . . . . . . . 40
3.3 Application . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42
Introduction and Generalities
In physical sciences (e.g, elasticity, astronomy) and natural sciences (e.g, ecology) among others,
the consideration of periodic environmental factor in the dynamics of multi-phenomena
interractions (or multi-species interractions) leads to the study of dierential systems with
periodic data. Therefore, it is worth investigating, the fundamental questions inherent in systems
of periodic (ordinary) dierential equations such as: existence, uniqueness and stability.
This work is specially devoted to one of the celebrated tools which is crucial in the analysis
and control of time-periodic systems, namely Floquet Theory. This basic theory which is
about a systematic study of linear systems of ordinary dierential equations with periodic
coecients, has had a rich history since the pioneering work of Floquet (1883) followed by the
contribution of Lyapunov (1892). Futhermore it gives a change of variables that transforms a
linear dierential system with periodic coecients, into a system with constant coecients,
and also provides a representation formula of the solutions of (1).
The aim of this project is to present Floquet theory and to use it to assess stability of
periodic solution of periodic dierential systems (linear or nonlinear).
The work is organized as follows:
First, we review some basic notions from Analysis and Linear algebra that will be used in
the subsequent chapters.
In chapter 2, we present the basic notions and theorems of the general theory of Ordinary
Dierential Equations (ODE).
Finally in chapter 3, we state Floquet’s theorem, prove it and use it to address stability of
periodic solutions of linear as well as nonlinear periodic systems. Furthermore, we illustrate
our method by studying Hill’s equation which is a generalization of Mathieu equation.
The purpose of this chapter is to give some concepts and theorems of linear algebra and
Mathematical analysis that will be of importance in the subsequent chapters
1.1 Basic concepts from linear functional analysis
let’s recall some denitions and results from linear functional analysis
Denition 1.1.1 Let X be a linear space over a eld K; where K holds either for Rn or C.
A mapping k:k: X ! R is called a norm provided that the following conditions hold:
i) kxk 0; for all x 2 X and kxk= 0 , x = 0
ii) kxk= jjkxk for all 2 K; x 2 X
iii) kx + yk kxk+kyk for arbitrary x; y 2 X
If X is a linear space and k:k is a norm on X; then the pair (X; k:k) is called a normed
linear space over K.
Should no ambiguity arise about the norm, we simply abbreviate this pair by saying that X
is a normed linear space over K.
Example 1. Each of the following expressions denes on the vector space Rn a norm which
is in common use.
i) The absolute norm :
jxij ; for every x = (x1; : : : ; xn) 2 Rn:
ii) The euclidean norm :
; for every x = (x1; : : : ; xn) 2 Rn:
iii) The maximum norm :
kxk1 = max
jxij; for every x = (x1; : : : ; xn) 2 Rn:
Example 2. Let X = C([0; 1]) be the space of all real-valued continuous functions on [0; 1].
Each of the following expressions denes on the vector space C([0; 1]) a norm which is in
0 jf(t)jdt for every f 2 C([0; 1]).
2 for every f 2 C([0; 1]).
iii) kfk1= max
jf(t)j: t 2 [0; 1]
Denition 1.1.2 (Equivalent norms)
Two norms k:k1 and k:k2 dened on a normed linear space X are said to be equivalent if
there exists > 0 and > 0 constants such that
kxk1 kxk2 kxk1 8x 2 X:
Theorem 1.1.3 In a nite dimensional normed linear space, all the norms are equivalent.
Denition 1.1.4 Every normed linear space E is canonically endowed with a metric d de-
ned on E E by
d(x; y) = jjx yjj 8 x; y 2 E:
Denition 1.1.5 (Cauchy sequence)
A sequence (xn)n1 of elements of a normed vector space X is a Cauchy sequence if
kxn xmk= 0:
That is, for any > 0 there is an interger N = N() such that kxn xmk< whenever
n N and m N.
Remark. In a normed linear sapce, every Cauchy sequence (xn)n1 is bounded; i.e, there
exists a constant M 0 such that jjxnjj M ; 8n 1: (See also Denition 1.1.11 below)
Denition 1.1.6 (convergent sequence)
A sequence (xn)n1 of elements of a normed vector space X converges to an element x 2 X
kxn xk= 0:
In such a case, we say that (xn)n1 is a convergent sequence.
Remark. In a normed linear space, every convergent sequence is a Cauchy sequence.
Denition 1.1.7 A normed linear space is complete if every Cauchy sequence in X has the
limit in X. A complete normed linear space is called a Banach space.
Remark. The notion of completeness is also dened for metric spaces which need not have
any linear structure.
Example (Banach space). The normed linear space
C([0; 1]); k k1
is a Banach space.
Denition 1.1.8 (Open sets and closed sets)
Let X be a normed linear space. We dene open (respectively closed) ball with center at a
point x 2 X and radius r > 0 by
Br(x) = fx 2 X : kxk< rg (respectively Br(x) = fx 2 X : kxk rg ) :
A nonempty subset A of a normed linear space X is said to be open if for all x 2 A; there
exists r > 0 such that Br(x) A: And A is said to be closed if XnA is open.
Proposition 1.1.9 A subset A of a normed linear space is closed if and only if every con-
vergent sequence (an)n1 of elements of A has its limit in A:
Denition 1.1.10 (Closure and interior) Let A be a subset of a normed linear space X.
The interior of A denoted by intA is dened as the union of all open sets contained in A
and the closure of A denoted by cl(A) or A is dened as the intersection of all closed sets
Theorem 1.1.11 Let A be a subset of a normed linear space X and x 2 X then,
a) x 2 intA if and only if 9 r > 0 : B(x; r) A:
b) x 2 clA if and only if 8 r > 0; B(x; r) \ A 6= ?:
Remark. Given a subset A of a normed linear space X; we have :
x 2 A () 9 (an)n A such that lim
an = x:
Denition 1.1.12 Let X be a normed linear space, x 2 X and let V be a subset of X
containing x: We say that V is a neighbourhood of x if there exists an open set U of X
containing x and contained in V: We denote by N(x) the collection of all neighbourhoods of
Denition 1.1.13 A subset of a normed linear space X is said to be bounded if it can be
included in some ball.
Theorem 1.1.14 (Riesz/ Heine-Borel) A normed linear space is nite dimensional if
and only if its closed unit ball is compact, i.e., every bounded sequence of the closed unit ball,
has a convergent subsequence.
1.1.1 Linear operators
In this section X and Y are normed linear spaces over K.
Denition 1.1.15 A K-linear operator T from X into Y is a map T : X ! Y such that
T(x + y) = Tx + Ty
for all ; 2 K and all x; y 2 X:
When Y = K, such a map is called a linear functional or a linear form.
Proposition 1.1.16 The set of K-linear operators from X into Y has a natural structure
of linear space over K and is denoted by L(X; Y ). Note that L(X;X) is simply denoted by
Proposition 1.1.17 If Z is also a linear space, then
f 2 L(X; Y ) and g 2 L(Y;Z) =) gof 2 L(X;Z) :
Theorem 1.1.18 Let T 2 L(X; Y ). Then the following are equivalent
i) T is continuous at the origin (in the sense that if fxngn is a sequence in X such that
xn ! 0 as n ! 1, then T(xn) ! 0 in Y as n ! 1.
ii) T is Lipschitz, i.e., there exists a constant K 0 such that for every x 2 X,
jjT(x)jj Kjjxjj :
iii) The image of the closed unit ball, T
, is bounded.
Denition 1.1.19 A linear operator T : X ! Y is said to be bounded if there exist some
k 0 such that
for all x 2 X:
If T is bounded, then the norm of T is dened by
kTk= inffk : kT(x)k kxk; x 2 Xg:
The set of bounded linear operators from X into Y is denoted B(X; Y ): If X = Y; one simply
Proposition 1.1.20 Suppose X 6= f0g and T 2 B(X), then we have the following:
1.2 Matrix calculus and basic Operator theory
1.2.1 Matrices, eigenvalues and eigenvectors
Denition 1.2.1 An m n matrix A is a rectangular array of numbers, real or complex,
with m rows and n columns. We shall write aij for the number that appears in the ith row
and the jth column of A; this is called the (i; j) entry of A: We can either write A in the
extended form 0
a11 a12 a1n
a21 a22 a2n
. . .
am1 am2 amn
or in the more compact form
We will denoted A = (aij)m;n with 1 i m; 1 j n:
Associated with each matrix A a matrix At; known as the transpose of A, and obtained
from A by interchanging the rows and the columns of A. Thus, if A = (aij)m;n then At =
(aji)n;m: The trace denoted by tr of A is the sum of a diagonal elements of A. Two m n
matrices A and B are said to be equal if all corresponding elements are equal, that is , if
aij = dij for each i and j.
Let D = (dij)m;n, A = (aij)m;n and B = (bij)n;p be a matrices.
1) The sum of two m n matrices A and D is dened as matrix obtained by adding
D + A = (dij + aij)m;n:
Similarly, the dierence is
D A = (dij aij)m;n:
2) The multiplication of a matrix A by a scalar is dened as follows: A = (aij)m;n:
3) The product of a m n matrix A and a n p matrix B is
i) A row-vector or 1 m matrix: A = (a1; a2; :::; an) where the a0
is are scalars.
ii) A column-vector or n 1 matrix: A =
where the a0
is are scalar.
iii) A zero matrix is a matrix which all of whose entries are zero. The zero m n matrix
is denoted by 0m;n or simply 0:
iv) A square matrix is a matrix with the same number of rows and columns.
v) The identity matrix of order n has one on the principal diagonal, that is from top
left to bottom right, and zeros elsewhere; it is denoted by In = (ij) where ij is the
Kronecker’s symbol. From the denition of matrix multiplication we have
AI = IA = A
for any square matrix A.
iv) A square matrix A is regular or non-singular if its column vectors are linearly independant
or equivalently its determinant, det(A), is nonzero. Otherwise A is said to be
singular or degenerated when det(A) = 0. (See details/recalls below).
v) A square matrix in which all the non-zero elements lie on the principal diagonal is
called a diagonal matrix.
vi) An n n matrix N is said to be nilpotent, if there is a positive integer k such that
Nk = 0:
Denition 1.2.2 Let A and B be two n n matrices. We say that A and B are similar,
notation A B; if and only if there exists a nonsingular matrix T such that T1AT = B
Eigenvalues and Eigenvectors
Denition 1.2.3 Let A be an nn square matrix. Then a scalar is called an eigenvalue
of A, if there exists a nonzero vector v 2 Rn such that Av = v:
In this case, the is called an eigenvalue of A and the v; eigenvector associated with :
The eigenvalues of A are also the roots of the characteristic polynomial p() = det(AI)
with p() of degree n:
Theorem 1.2.4 Let A be a square matrix, if the eigenvalues of A are all distincts, A is
similar to a diagonal matrix (whose diagonal entries are the eigenvalues of A).
Now if the matrix A has repeated eigenvalues, then it is not possible to diagonalize it. In
this case, we introduce the concept of generalized eigenvector.
Denition 1.2.5 A nonzero vector v is called a generalized eigenvector of rank k of A;
associated with an eigenvalue if and only if
(A I)kv = 0 and (A I)k1v 6= 0:
Lemma 1.2.6 If v is a generalized eigenvector of rank k; then the vectors
v; (A I)v; :::; (A I)k1v are linearly independent.
Recall that a set of vectors v1; v2; : : : ; vk is linearly dependent if there exist scalars c1; c2; : : : ; ck
not all zero, such that
c1v1 + c2v2 + : : : + ckvk = 0:
A set of vectors v1; v2; : : : ; vk is linearly independent if it is not linearly dependent.
From the lemma above, we construct a new basis for Cn such that the matrix representation
of A with respect to this new basis is the one we call Jordan canonical form denoted by J:
Theorem 1.2.7 For every nn complex matrix A with eigenvalues 1; :::; s (not necessarily
distinct) of multiplicities n1; :::; ns respectively, there exists a nonsingular n n matrix P
P1AP = J = diag(J1; :::; Js);
where each of block matrices J1; :::; Js is of the form
k 1 0 0
0 k 1 0
0 k 1
0 0 k
; k = 1; :::; s;
k=1 nk = n:
For the proof of this theorem, refer to Coddington and Levinson [C/L] or Hirsch and Smale
The block matrices J1; :::; Js are called Jordan blocks, and J is called the Jordan canonical
form of A: Note, any Jordan block Jk() can be written as Jk = kI + Nk; where Nk is
nilpotent of order k and that A is similar to J:
1.2.2 Limits of sequences of operators
We introduce the concepts of limit in the norm and of strong limit of a sequence of operators.
We shall then introduce the derivative and integral of operators depending on a parameter
and shall discuss series of operators.
Let X be a Banach space and (An) a sequence of operators in L(X):
– We say that (An) converges in norm to the operator A 2 L(X); if
jjAn Ajj = 0: (1.2.1)
– If, for each element x 2 X;
jjAnx Axjj = 0; (1.2.2)
we shall say that An converges strongly to the operator A 2 L(X):
It is immediate that (1.2.1) implies (1.2.2). The converse is not in general true. However, in
the case of nite dimensional Banach spaces it is true, for if we put x = e(i); we have from
i1 ; :::; a(k)
ik )jj = 0 for i = 1; :::; k:
This ensures that all the components tend to 0; and since there are only a nite number of
them, the limit is uniform.
It is now possible to dene the meaning of the series
As is said to be convergent if the sequence made by the partial sums
converges in L(X).
In this case of matrices,
As corresponds to a matrix of which elements are
ij = aij :
– We say that the series of operators
As converges absolutely, if the series
In the case of matrices, this happens if and only if, for every component,
converges with As = (a(s)
Proposition 1.2.8 Let X be a normed linear space. If A 2 L(X); then the series
Proof. If suces to show that the sequence of partial sums fSNg1N
=1 for the series
n! is a Cauchy sequence. Let us dene
n! : Note that the partial sums of the convergent
series of real number
n! = ekAk form a cauchy sequence. Using this fact, it follows
that SN is a cauchy sequence in L(X).
Dene the exponential map exp : L(X) ! L(X) by exp(A) = eA =
The main properties of exponential map are summarized in the following proposition
Proposition 1.2.9 Suppose that A;B 2 L(Rn):
i) If A 2 L(Rn; ) then eA 2 L(Rn):
ii) If B is nonsingular, then B1eAB = eB1AB:
iii) eA = (eA)1:
iv) keAk ekAk:
– Let A be an operator depending on a real parameter t with a t b: Let t0 2 [a; b]; if
h is suciently small, the operator
A(t0 + h) A(t0)
can be dened. If its limit as h ! 0 exists, we say that the operator A is dierentiable at t
with respect to t: The limiting operator is denoted by d
dtA(t): We thus have
A(t0 + h) A(t0)
In the case of matrices, the limits exists if and only if each component aij is dierentiable,
and we have
In general, if A and B are two dierentiable operators in L(X); the following product rule
B + A
We say that the series of operators depending on a parameter converges uniformly if, for
every > 0; there exists a > 0 such that for every t 2 [a; b] and every >
If the operator As(t) are dierntiable and the series
converges uniformly, then also the operator A =
As(t) is dierntiable, and we have
We illustrate this notion with the following result: If A 2 L(X) then tA 2 L(X) for each
t 2 R so that the function t 7! etA is dierentiable and
(etA) = AetA:
If A is a nonsingular matrix, then logarithm of A denoted by ln(A) is well-dened matrix.
This important result is stated in the following theorem.
Theorem 1.2.10 Let A be a nonsingular n n matrix, then there exists an n n matrix
B (called logarithm of A) such that A = eB:
Let A(t) be a matrix depending on a parameter t; we suppose that the component aij are
all integrable functions over the interval [t0; t]: We shall say that the matrix
is the intergral of the matrix A(t) between t0 and t:
According to the following denition
we have that if the functions aij are continuous, then
A( )d = A(t):
If, nally A is an m n matrix, and L is a number greater than or equal to the absolute
values of aij ; then
1.3 Review of calculus
We now give the general denition of dierentiability. Let U be a nonempty open subset of
a Banach space X, let Y denote a Banach space, and let the symbol k:k denote the norm in
both Banach spaces.
Denition 1.3.1 A function f : U ! Y is dierentiable at x 2 U if there is a map
A 2 L(X; Y ) such that
kf(x + h) + f(x) Ahk
Remark. If such a linear map exists, then it is unique and we write it as A = f0(x) called
the derivative of f at x:
Other common notations for the derivative are Df and fx
The following is a list of standard facts about the derivative. For the statement in the list,
the symbols X; Y; Xi; and Yi denoted Banach spaces.
i) If f : X ! Y is dierentiable at a 2 X, then f is continuous at a.
ii) If f : X ! Y1 ::: Yn is given by f(x) = (f1(x); :::fn(x)); and if fi is dierentiable
for each i, then so is f and,
Df(x) = (Df1(x); :::;Dfn(x)):
iii) If the function f : X1 X2 ::: Xn ! Y is given by (x1; x2; :::; xn) 7!
f(x1; x2; :::; xn), then the ith partial derivative of f at (a1; a2:::; an) 2 X1X2:::Xn
is the derivative of the function g : Xi ! Y dened by g(xi) = f(a1; :::; ai1; xi; ai+1; :::; an).
This derivatives is denoted by Dif(a): If f is dierentiable, then all its partial derivatives
exist and, if we dene h = (h1; h2; :::; hn); we have
The converse is not true in general, but if all the partial derivatives of f exist and are
continuous in an open set U X1 X2 :::Xn then f is continuously dierentiable
1.3.1 The mean value theorem
Theorem 1.3.2 Suppose that [a; b] is a closed interval, and f : [a; b] ! Y is a continuous
function. If f is dierentiable on the open interval (a; b) and there is some number M > 0
such that kf0(t)k M for all t 2 (a; b); then
kf(b) f(a)k M(b a)
Theorem 1.3.3 (Mean value theorem) Suppose that f : X ! Y is dierentiable on an
open set U X with a; b 2 U and a + t(b a) 2 U for 0 t 1: If there is some M > 0
kDf(a + t(b a)k M;
kf(b) f(a)k Mkb ak:
1.3.2 Integration in Banach space
Let X be a Banach space and I = [a; b] R where a < b.
Denition 1.3.4 If f : I ! X is continuous on I, we dene its integral in the sense of
Rieman by the following formula:
f(t)dt = lim
f(a + k
It is easily seen that
Theorem 1.3.5 Suppose that U is a nonempty open subset of X: If f : X ! Y is a
dierentiable function, and x + ty 2 U for 0 t 1; then
f(x + y) f(x) =
Df(x + ty)ydt: