NonlinearBook10pt November 20, 2007
ROBUST NONLINEAR CONTROL 661
Proof. The result is a direct consequence of Theorem 11.1 with f (x) =
(A + ∆A)x, f
0
(x) = Ax, L(x) = x
T
Rx, V (x) = x
T
P x, Γ(x) = x
T
Ω(P )x,
and D = R
n
. Specifically, conditions (11.6) and (11.7) are trivially satisfied.
Now, V
′
(x)f(x) = x
T
(A
T
P + P A)x + x
T
(∆A
T
P + P ∆A)x, and hence,
it follows fr om (11.52) that V
′
(x)f(x) ≤ V
′
(x)f
0
(x) + Γ(x) = x
T
(A
T
P +
P A + Ω(P ))x, for all ∆A ∈ ∆
A
. Furthermore, it follows from (11.53) that
L(x)+V
′
(x)f
0
(x)+Γ(x) = 0, and hence, V
′
(x)f
0
(x)+Γ(x) < 0 for all x 6= 0,
so that all the conditions of Theorem 11.1 are satisfied. Finally, since V (x),
x ∈ R
n
, is rad ially unbounded, (11.50) is globally asymptotically stable for
all ∆A ∈ ∆
A
.
Corollary 11.2 is the deterministic version of Theorem 4.1 of [50]
involving quadratic Lyapunov bounds for addr essin g robust stability and
performance analysis of linear uncertain systems. For convenience we shall
say that Ω(·) bounds ∆
A
if (11.52) is satisfied. To apply Corollary 11.2,
we fir st specify a function Ω(·) and an uncertainty set ∆
A
such that Ω(·)
bounds ∆
A
. If the existence of a positive-definite solution P to (11.53) can be
determined analytically or numerically, then robust stability is guaranteed
and the performance bound (11.54) can be computed. We can then enlarge
∆
A
, modify Ω(·), and again attempt to solve (11.53). If, however, a
positive-definite solution to (11.53) cannot be determined, then ∆
A
must be
decreased in size until (11.53) is solvable. For example, Ω(·) can be replaced
by εΩ(·) to bound ε∆
A
, where ε > 1 enlarges ∆
A
and ε < 1 shrinks ∆
A
.
Of course, the actual range of uncertainty that can be bounded depends on
the nominal matrix A, the function Ω(·), and the structure of ∆
A
.
Since the ord ering induced by the cone of nonnegative-definite matrices
is only a partial ordering, it should not be expected that there exists an
operator Ω(·) satisfying (11.52), which is a least upper bound. Next, we
illustrate two bounding functions for two different uncertainty characteriza-
tions. Specifically, we assign exp licit structure to the uncertainty set ∆
A
and the bound Ω (·) satisfying (11.52). First, we assume that the uncertainty
set ∆
A
to be of the form
∆
A
=
(
∆A ∈ R
n×n
: ∆A =
p
X
i=1
δ
i
A
i
,
p
X
i=1
δ
2
i
α
2
i
≤ 1
)
, (11.56)
where for i = 1, . . . , p, A
i
is a fixed matrix denoting the str ucture of the
parametric uncertainty; α
i
is a given positive number; and δ
i
is an uncertain
parameter. Note that the uncertain parameters δ
i
are assumed to lie in a
specified ellipsoidal region in R
p
. For the structure ∆
A
as specified by
(11.56), the bound Ω(·) satisfying (11.52) can now be given a concrete form.