Lecture 4 Chapter 4 Lyapunov Stability
Lecture 4 Chapter 4 Lyapunov Stability
Lecture 4 – p. 1/86
Autonomous Systems
Consider the autonomous system
ẋ = f (x) (1)
Lecture 4 – p. 2/86
Stability
Definition 4.1: The equilibrium point x = 0 of (1) is
stable if, for each ǫ > 0, there is δ = δ(ǫ) > 0 such that
Lecture 4 – p. 3/86
Stability
“Stability is a property of the equilibrium, not of the system"
Stability of the equilibrium is equivalent to stability of the
system only when there exists only one equilibrium (e.g.,
linear systems). In this case stability ≡ global stability.
The equilibrium point x = 0 of (1) is
attractive if there is δ > 0 such that
Lecture 4 – p. 4/86
Derivative along the trajectory
Let V : D → R be a continuously differentiable
Definition:
function defined in a domain D ∈ Rn that contains the
origin. The derivative of V along the trajectory (solution) of
(1), denoted by V̇ (x) is given by
n
X ∂V
V̇ (x) = ẋi
∂xi
i=1
n
X ∂V
= fi (x)
∂xi
i=1
∂V ∂V ∂V
= [ , ,..., ][f1 (x), f2 (x), . . . , fn (x)]T
∂x1 ∂x2 ∂xn
∂V
= f (x)
∂x
Lecture 4 – p. 5/86
Lyapunov Stability Theorem
Theorem 4.1:Let x = 0 be an equilibrium for (1) and D ∈ Rn
be a domain containing x = 0. Let V : D → R be a
continuously differentiable function, such that
V̇ (x) ≤ 0 in D (3)
Lecture 4 – p. 6/86
Lyapunov Stability Theorem
Proof:
Lecture 4 – p. 7/86
Lyapunov Stability Theorem
Lyapunov function candidate
Lecture 4 – p. 8/86
Lyapunov Stability Theorem
Lyapunov surface (level surface, level set)
{x|V (x) = c}
Lecture 4 – p. 9/86
Lyapunov Stability Theorem
Positive definite
V (0) = 0, V (x) ≥ 0, ∀x 6= 0
V (x) is negative (semi)definite if −V (x) is positive
(semi)definite
Lyapunov Theorem
V pdf + V̇ nsdf → stable
V pdf + V̇ ndf → asymptotically stable
Lecture 4 – p. 10/86
Lyapunov Stability Theorem
Example 4.4: Consider the pendulum equation with friction
ẋ1 = x2
g
k
ẋ2 = − sin x1 − x2
l m
Lecture 4 – p. 11/86
Region of Attraction
When the origin x = 0 is asymptotically stable, we are often
interested in determining how far from the origin the
trajectory can be and still converge to the origin as t → ∞.
This gives rise to the definition of region of attraction (also
called region of asymptotically stability, domain of
attraction, or basin).
Definition: Let φ(t, x) be the solution of (1) that starts at initial
state x at time t = 0. The, the region of attraction is defined
as the set of all points x such that limt→∞ φ(t, x) = 0
Question: Under what conditions will the region of attraction
be the whole space Rn ? In other words, for any initial state
x, under what conditions the trajectory φ(t, x) approaches
the origin as t → ∞, no matter how large kxk is. If an a.s.
equilibrium point at the origin has this property, it is said to
be globally asymptotically stable (g.a.s.).
Lecture 4 – p. 12/86
Global Lyapunov Stability Theorem
Let x = 0 be an equilibrium for (1). Let
Theorem 4.2:
V : Rn → R be a continuously differentiable function, such
that
V (0) = 0 and V (x) > 0, ∀x 6= 0 (5)
Lecture 4 – p. 13/86
Chetaev’s Instability Theorem
Theorem 4.3: Let x = 0 be an equilibrium for (1). Let
V : D → R be a continuously differentiable function, such
that V (0) = 0 and V (x0 ) > 0 for some x0 with arbitrarily
small kx0 k. Define a set
U = {x ∈ Br |V (x) > 0}
where
Br = {x ∈ Rn |kxk ≤ r}.
Suppose that V̇ (x) > 0 in U . Then x = 0 is unstable.
Lecture 4 – p. 14/86
Chetaev’s Instability Theorem
Proof:
ẋ1 = x1 + g1 (x)
ẋ2 = −x2 + g2 (x)
Lecture 4 – p. 16/86
Invariance Principle
Example 4.4: Consider the pendulum equation with friction
ẋ1 = x2
g
k
ẋ2 = − sin x1 − x2
l m
x22
g
k
V (x) = − (1 − cos x1 ) + ⇒ V̇ (x) = − x22
l 2 m
Lecture 4 – p. 17/86
Invariance Principle
Idea: If we can find a Lypunov function in a domain
containing the origin whose derivative along the trajectories
of the system is negative semidefinite, and if we can
establish that no trajectory can stay identicaly at points
where V̇ (x) = 0 except at the origin, then the origin is
asymptotically stable (LaSalle’s Invariance Principle).
Lecture 4 – p. 18/86
Invariance Principle
Let x(t) be a solution of the autonomous system ẋ = f (x).
Definition: The point p is a positive limit point of x(t) if exists
a sequence {tn } with tn → ∞ as n → ∞ such that x(t) → p
as n → ∞.
Definition: The set of all positive limit points of x(t) is called
the positive limit set of x(t).
Definition: A set M is (positively) invariant w.r.t. ẋ = f (x) if
x(0) ∈ M ⇒ x(t) ∈ M for all t ∈ R (t ∈ R+ ).
Definition: x(t) approaches M as t → ∞ if for each ǫ > 0
exists T > 0 such that dist(x(t), M )(= inf x∈M kp − xk) < ǫ for
all t > T .
x(t) → M as t → ∞ ⇒ limt→∞ dist(x(t), M ) = 0
x(t) → M as t → ∞ NOT ⇒ ∃ limt→∞ x(t)
Lecture 4 – p. 19/86
Invariance Principle
Lemma 4.1: If x(t), solution of ẋ = f (x), is bounded for all
t ≥ 0, then it has a nonempty positive limite set L+ which is
compact and invariant. Moreover,
x(t) → L+ as t → ∞
Lecture 4 – p. 20/86
Invariance Principle
Theorem 4.4 (LaSalle’s Theorem): Let Ω ⊂ D be a compact set
that is positively invariant w.r.t. ẋ = f (x). Let V : D → R be
a continuously differentiable function such that V̇ (x) ≤ 0 in
Ω. Let E be the set of all points in Ω where V̇ (x) = 0. Let M
be the largest invariant set in E . Then, every solution
starting in Ω approaches M as t → ∞.
Lecture 4 – p. 21/86
Invariance Principle
Proof:
Lecture 4 – p. 22/86
Invariance Principle
Corollary 4.1 (4.2):Let x = 0 be an equilibrium point of
ẋ = f (x). Let V : D(Rn ) → R be a continuously
differentiable (radially unbounded, positive definite) function
on a domain D containing the origin (on Rn ), such that
V̇ (x) ≤ 0 in D (in Rn ). Let S = {x ∈ D(Rn )|V̇ (x) = 0} and
suppose that no solution can stay identically in S , other
than the trivial solution. Then, the origin is (globally)
asymptotically stable.
Lecture 4 – p. 23/86
Invariance Principle
LaSalle’s theorem
relaxes the negative definiteness requirement for V̇ (x)
of Lyapunov’s theorem
does not require V (x) to be positive definite
gives an estimate of the region of attraction Ω which is
not necessarily a level set of V (x), i.e.,
Ωc = {x ∈ Rn |V (x) ≤ c}
applies not only to equilibrium points but also to
equilibrium sets
Lecture 4 – p. 24/86
Invariance Principle
Example:
ẋ = −|x|x + (1 − |x|)xy
1
ẏ = − (1 − |x|)x2
8
Lecture 4 – p. 25/86
Linear Systems and Linearization
The linear time-invariant system
ẋ = Ax
x = eAt x(0)
Lecture 4 – p. 26/86
Linear Systems and Linearization
For all A, there exists a nonsingular transformation P such
that
J1
−1
P AP = J =
. . .
Jr
where Ji is a Jordan block of order mi associated with
eigenvalue λi of A, i.e.,
λi 1 0 ... ... 0
0 λi 1 0 . . . 0
. ... ..
..
.
Ji = .
.. ...
0
..
...
. 1
Lecture 4 – p. 27/86
Linear Systems and Linearization
Then,
r X
X mi
eAt = P eJt P −1 = tk−1 e(λi t) Rik
i=1 k=1
Lecture 4 – p. 28/86
Linear Systems and Linearization
Theorem 4.5: The equilibrium point x = 0 of ẋ = Ax is
1. stable ⇔ Reλi ≤ 0 and for every eigenvalue with
Reλi = 0 and algebraic multiplicity qi ≥ 2,
rank(A − λi I) = n − qi , where n is the dimension of x
2. (globally) asymptotically stable ⇔ Reλi < 0 for all i
Proof:
If Reλi > 0, e(λi t) → ∞ as t → ∞. Moreover, an
eigenvalue on the imaginary axis (Reλi = 0) could give
rise to unbounded terms if order of associated Jordan
block is higher than one (rank(A − λi I) 6= n − qi ) due to
the term tk−1 .
If Reλi < 0, e(λi t) → 0 as t → ∞, i.e., e(At) → 0 as t → ∞
(a.s.). Since x(t) depends linearly on x(0), a.s. of origin
is global.
Lecture 4 – p. 29/86
Linear Systems and Linearization
When all eigenvalues of A satisfy Reλi < 0, A is called a
Hurwitz matrix or a stability matrix. The origin of ẋ = Ax is
a.s. if and only if A is Hurwitz.
Lyapunov function candidate
V (x) = xT P x, P > 0, P = PT
V̇ (x) = xT (P A + AT P )x = −xT Qx
Lyapunov equation
P A + AT P = −Q, Q > 0, Q = QT
Lecture 4 – p. 30/86
Linear Systems and Linearization
Theorem 4.6: A matrix A is a stability matrix or Hurwitz
matrix, that is, Reλi < 0 for all eigenvalues of A, if and only
if for any given positive definite symmetric matrix Q there
exists a positive definite symmetric matrix P that satisfies
the Lyapunov equation. Moreover, if A is stability matrix or
Hurwitz matrix, the P is the unique solution of the Lyapunov
equation.
Lecture 4 – p. 31/86
Linear Systems and Linearization
Proof:
Lecture 4 – p. 32/86
Linear Systems and Linearization
What if Q = C T C is positive semidefinite?
Lecture 4 – p. 33/86
Linear Systems and Linearization
Proof:
Lecture 4 – p. 34/86
Linear Systems and Linearization
Theorem 4.7 (Lyapunov’s First Method or Lyapunov’s Indirect
Method): Let x = 0 be an equilibrium point for the nonlinear
system
ẋ = f (x)
where f : D → Rn is continuously differentiable and D is a
neighborhood of the origin. Let
∂f
A= (x)
∂x x=0
Lecture 4 – p. 36/86
Nonautonomous Systems
Consider the nonautonomous system
ẋ = f (t, x) (8)
Lecture 4 – p. 37/86
Nonautonomous Systems
Suppose ȳ(τ ) is a solution of the system
dy
= g(τ, y)
dτ
defined for all τ ≥ a. The change of variables
x = y − ȳ(τ ), t=τ −a
˙ + a) , f (t, x)
˙ ) = g(t + a, x + ȳ(t + a)) − ȳ(t
ẋ = g(τ, y) − ȳ(τ
Since
˙ + a) = g(t + a, ȳ(t + a)),
ȳ(t ∀t ≥ 0
the origin x = 0 is an equilibrium point of the transformed
system at t = 0.
Lecture 4 – p. 38/86
Nonautonomous Systems
By examining the stability behavior of the origin as an
equilibrium point for the transformed system, we
determine the stability behavior of the solution ȳ(τ ) of
the original system.
Lecture 4 – p. 39/86
Stability
While the solutions of autonomous systems depend only on
t − t0 , the solutions of nonautonomous systems may
depend on both t and t0 . We then refine definitions.
Definition 4.4:The equilibrium point x = 0 of (8) is
stable if, for each ǫ > 0, there is δ = δ(ǫ, t0 ) > 0 such that
Lecture 4 – p. 40/86
Stability
asymptotically stable (a.s) if stable and attractive
uniformly asymptotically stable (u.a.s) if u.s. and u.a.
globally uniformly asymptotically stable (g.u.a.s) if u.a.s.
and c arbitrarily large
Lecture 4 – p. 41/86
Comparison Functions
The solutions of nonautonomous systems depend on both t
and t0 . We refine the stability definitions, so that they hold
uniformly in the initial time t0 , using special comparison
functions.
Definition 4.2: A continuous function α : [0, a) → [0, ∞) is said
to belong to class K if it is strictly increasing and α(0) = 0. It
is said to belong to class K∞ if a = ∞ and α(r) → ∞ as
r → ∞.
Definition 4.3: A continuous function
β : [0, a) × [0, ∞) → [0, ∞) is said to belong to class KL if for
each fixed s, the mapping β(r, s) belongs to class K with
respect to r and, for each fixed r, the mapping β(r, s) is
decreasing with respect to s and β(r, s) → 0 as s → ∞.
Lecture 4 – p. 42/86
Comparison Functions
Examples:
α(r) = rc ∈ K if c > 0
α(r) = rc ∈ K∞ if c > 0
β(r, s) = rc e−s ∈ KL if c > 0
Lecture 4 – p. 43/86
Comparison Functions
The next lemma states properties of class K and class KL
functions
Lemma 4.2: Let α1 and α2 be class K functions on [0, a), α3
and α4 be class K∞ functions, and β be a class KL
function. Denote the inverse of αi by αi−1 . Then,
α1−1 is defined on [0, α1 (a)) and belongs to class K
α3−1 is defined on [0, ∞) and belongs to class K∞
α1 ◦ α2 belongs to class K
α3 ◦ α4 belongs to class K∞
σ(r, s) = α1 (β(α2 (r), s)) belongs to class KL
Lecture 4 – p. 44/86
Comparison Functions
The next lemmas show how class K and class KL enter
into Lyapunov analysis
Lemma 4.3: Let V : D → R be a continuous positive definite
function defined on a domain D ⊂ Rn that contains the
origin. Let Br ⊂ D for some r > 0. Then, there exist class K
functions α1 and α2 , defined on [0, r), such that
Lecture 4 – p. 45/86
Comparison Functions
For a quadratic positive definite function V (x) = xT P x,
Lemma 4.3 follows from the inequality
y(t) = σ(y0 , t − t0 )
Lecture 4 – p. 46/86
Stability
Lemma 4.5: The equilibrium point x = 0 of (8) is
∂V ∂V
+ f (t, x) ≤ 0 (14)
∂t ∂x
for all t ≥ 0 and for all x ∈ D, where W1 (x) and W2 (x) are
continuous positive definite functions on D. Then, x = 0 is
uniformly stable (u.s.).
Lecture 4 – p. 48/86
Stability
Proof:
Lecture 4 – p. 49/86
Stability
Theorem 4.9: Let x = 0 be an equilibrium point for (8) and
D ⊂ Rn be a domain containing the origin. Let
V : [0, ∞) × D → R be a continuous differentiable function
such that
W1 (x) ≤ V (t, x) ≤ W2 (x) (15)
∂V ∂V
+ f (t, x) ≤ −W3 (x) (16)
∂t ∂x
for all t ≥ 0 and for all x ∈ D, where W1 (x), W2 (x) and W3 (x)
are continuous positive definite functions on D. Then, x = 0
is uniformly asymptotically stable (u.a.s.).
Lecture 4 – p. 50/86
Stability
Theorem 4.9 (cont’): Moreover, if r and c are chosen such
that Br = {kxk ≤ r} ⊂ D and c < minkxk=r W1 (x), then every
trajectory starting in {x ∈ Br |W2 (x) ≤ c} satisfies
Lecture 4 – p. 51/86
Stability
Proof:
Lecture 4 – p. 52/86
Stability
Theorem 4.10: Let x = 0 be an equilibrium point for (8) and
D ⊂ Rn be a domain containing the origin. Let
V : [0, ∞) × D → R be a continuous differentiable function
such that
k1 kxka ≤ V (t, x) ≤ k2 kxka (17)
∂V ∂V
+ f (t, x) ≤ −k3 kxka (18)
∂t ∂x
for all t ≥ 0 and for all x ∈ D, where k1 , k2 and k3 are
positive constants. Then, x = 0 is exponentially stable
(e.s.). Moreover, if the assumptions hold globally, then x = 0
is globally exponentially stable (e.s.).
Lecture 4 – p. 53/86
Stability
Proof:
Lecture 4 – p. 54/86
LTV Systems and Linearization
The stability behavior of the origin as an equilibrium point
for the linear time varying (LTV) system
ẋ = A(t)x (19)
Lecture 4 – p. 55/86
LTV Systems and Linearization
Theorem 4.11: The equilibrium point x = 0 of (19) is (globaly)
uniformly asymptotically stable if and only if the state
transition matrix satisfies the inequality
kΦ(t, t0 )k ≤ ke−λ(t−t0 ) , ∀t ≥ t0 ≥ 0
Lecture 4 – p. 56/86
LTV Systems and Linearization
For LTV systems, uniform asymptotic stability cannot be
characterized by the location of the eigenvalues of the
matrix A as it is done for LTI systems.
Example 4.22:
Lecture 4 – p. 57/86
LTV Systems and Linearization
Theorem 4.12: Let x = 0 be the exponentially stable
equilibrium point of (19). Suppose A(t) is continuous and
bounded. Let Q(t) be a continuous, bounded, positive
definite, symmetric matrix. Then, there is a continuously
differentiable, bounded, positive definite, symmetric matrix
P (t) that satisfies
Hence,
V (t, x) = xT P (t)x
is a Lyapunov function for the system that satisfies the
conditions of Theorem 4.10.
Note: It plays the role of Theorem 4.6 for LTV systems.
Lecture 4 – p. 58/86
LTV Systems and Linearization
Proof: Take
Z ∞
P (t) = ΦT (τ, t)Q(τ )Φ(τ, t)dτ
t
as gramian.
Lecture 4 – p. 59/86
LTV Systems and Linearization
Theorem 4.13:Let x = 0 be and equilibrium point for the
nonlinear system
ẋ = f (t, x)
where f : [0, ∞) × D → Rn is continuously differentiable,
D = {x ∈ Rn |kxk2 < r} and the Jacobian matrix [∂f /∂x] is
bounded and Lipschitz on D, uniformly in t. Let
∂f
A(t) = (t, x)
∂x x=0
ẋ = A(t)x
Lecture 4 – p. 60/86
LTV Systems and Linearization
Proof: Apply Theorem 4.12.
Lecture 4 – p. 61/86
Converse Theorems
Lyapunov Theorems:
Lecture 4 – p. 62/86
Converse Theorems
Theorem 4.14:Let x = 0 be an equilibrium point for the
nonlinear system
ẋ = f (t, x) (22)
Lecture 4 – p. 63/86
Converse Theorems
Then, there is a function V : [0, ∞) × D0 → R that satisfies
the inequalities
∂V ∂V
+ f (t, x) ≤ −c3 kxk2 (25)
∂t ∂x
∂V
≤ c4 kxk (26)
∂x
Lecture 4 – p. 65/86
Converse Theorems
Then, there is a function V : [0, ∞) × D0 → R that satisfies
the inequalities
∂V ∂V
+ f (t, x) ≤ −α3 (kxk) (30)
∂t ∂x
∂V
≤ α4 (kxk) (31)
∂x
Lecture 4 – p. 66/86
Converse Theorems
Theorem 4.15:Let x = 0 be an equilibrium point for the
nonlinear system
ẋ = f (t, x) (32)
Lecture 4 – p. 67/86
Converse Theorems
Let x = 0 be an equilibrium point for the
Corollary 4.3:
nonlinear system
ẋ = f (t, x) (35)
Lecture 4 – p. 68/86
Boundedness
Lyapunov analysis can be used to show boundedness of
the solution of the state equation, even when there is no
equilibrium point at the origin.
Example: ẋ = −x + δ sin(t), x(t0 ) = a, a > δ > 0
Lecture 4 – p. 69/86
Boundedness
Definition 4.6: The solutions of ẋ = f (t, x) are
uniformly bounded if there exists a positive constant c,
independent of t0 ≥ 0, and for every a ∈ (0, c), there is
β = β(a) > 0, independent of t0 such that
kx(t0 )k ≤ a ⇒ kx(t)k ≤ β, ∀t ≥ t0
kx(t0 )k ≤ a ⇒ kx(t)k ≤ b, ∀t ≥ t0 + T
∂V ∂V
+ f (t, x) ≤ −W3 (x), ∀ kxk ≥ µ > 0
∂t ∂x
for all t ≥ 0 and x ∈ D, where α1 , α2 are class K functions
and W3 (x) is a continuous positive definite function. Take
r > 0 such that Br ⊂ D and suppose that µ < α2−1 (α1 (r)).
Lecture 4 – p. 71/86
Boundedness
Then, there exists a class KL function β and for every initial
state x(t0 ), satisfying kx(t0 )k ≤ α2−1 (α1 (r)), there is T ≥ 0
(dependent on x(t0 ) and µ) such that the solution of
ẋ = f (t, x) satisfies
Lecture 4 – p. 72/86
Boundedness
Proof:
Lecture 4 – p. 73/86
Input-To-State Stability
Consider the system
ẋ = f (t, x, u) (39)
Lecture 4 – p. 74/86
Input-To-State Stability
Example: ẋ = Ax + Bu
Lecture 4 – p. 75/86
Input-To-State Stability
Suppose we have a Lyapunov function V (t, x) for the
unforced system (40) and let us calculate V̇ in presence of
u. Due to tthe boundedness of u, it is possible in some
cases to show that V̇ is negative outside a ball of radius µ,
where µ depends on sup kuk. This is expected when
f (t, x, u) is Lipschitz, i.e.,
Lecture 4 – p. 76/86
Input-To-State Stability
Definition 4.7:The system (39) is said to be input-to-state
stable (ISS) if there exist a class KL function β and a class
K function γ such that for any initial state x(t0 ) and any
bounded input u(t), the solution x(t) exists for all t ≥ t0 and
satisfies
kx(t)k ≤ β(kx(t0 )k, t − t0 ) + γ sup ku(τ )k (41)
t0 ≤τ ≤t
Lecture 4 – p. 77/86
Input-To-State Stability
The following Lyapunov-like theorem gives a sufficient
condition for ISS.
Theorem 4.19: Let V : [0, ∞) × Rn → R be a continuous
differentiable function such that
∂V ∂V
+ f (t, x) ≤ −W3 (x), ∀ kxk ≥ ρ(kuk) > 0
∂t ∂x
for all (t, x, u) ∈ [0, ∞) × Rn × Rm , where α1 , α2 are class K∞
functions, ρ is a class K function, and W3 (x) is a continuous
positive definite functions on Rn . Then, the system (39) is
input-to-state stable with γ = α1−1 ◦ α2 ◦ ρ.
The function V is called an ISS Lyapunov function.
Lecture 4 – p. 78/86
Input-To-State Stability
Proof: Direct application of Theorem 4.18 with µ = ρ(kuk).
Lecture 4 – p. 79/86
Input-To-State Stability
The following lemma is a direct consequence of the
Converse Lyapunov theorem for g.e.s. (Theorem 4.14)
Lemma 4.6: Suppose f (t, x, u is continuosly differentiable
and globally Lipschitz in (x, u), uniformly in t. If the unforced
system (40) has a globally exponentially stable equilibrium
point at the origin x = 0, then the system (39) is ISS.
Lecture 4 – p. 80/86
Input-To-State Stability
Converse also holds!
Theorem: For the system
ẋ = f (x, u)
Lecture 4 – p. 82/86
Input-To-State Stability
An interesting application of input-to-state stability arises in
the stability analysis of the cascade system
Lecture 4 – p. 83/86
Input-To-State Stability
Lemma 4.7: Under the stated assumptions, if the system
(43), with x2 as input, is input-to-state stable and the origin
of (44) is globally uniformly asymptotically stable, then the
origin of the cascade system (43)-(44) is globally uniformly
asymptotically stable.
Lecture 4 – p. 84/86
Input-To-State Stability
Proof:
Lecture 4 – p. 85/86
Input-To-State Stability
Suppose that in the system
Lecture 4 – p. 86/86