Z. Li et al. / Automatica 63 (2016) 287–291 289
3. Control design and analysis
Define the filtered tracking error s
i
(Slotine & Li, 1991)
s
i
= e
(n−1)
i
+ λ
i1
e
(n−2)
i
+ · · · + λ
i,n−1
e
i
, (15)
where e
i
= y
i
− y
di
, e
(n−1)
i
(i = 1, 2, . . . , m) is the (n −
1)th derivative of e
i
, λ
i1
, . . . , λ
i,n−1
are positive constants and are
appropriately chosen coefficient vectors such that e
i
→ 0 as s
i
→
0 (i.e. r
m−1
+ λ
i1
r
m−2
+ · · · + λ
i,n−1
is Hurwitz).
From (15), we have
˙
S = B
−1
d
(x)
F (x(t − τ (t))) + K
T
ΦV + ∆
P
+ ν, (16)
where S = [s
1
, . . . , s
i
, . . . , s
m
]
T
and ν = [ν
1
, . . . , ν
i
, . . . , ν
m
]
T
with
ν
i
= −y
(n)
di
+ λ
i1
e
(n−1)
i
+ · · · + λ
i,n−1
. (17)
We now construct a new high dimensional Lyapunov–
Krasovskii functional (see Eq. (31)). The first part of the Lyapunov–
Krasovskii functional is chosen as
V
1
= S
T
B
ϑ
S, (18)
where
B
ϑ
=
1
0
ϑB
α
dϑ = diag
1
0
ϑB
αii
(
¯
x
i
)dϑ
(19)
with B
α
= B
d
α = diag[b
dii
α
ii
]
m×m
and matrix α ∈ R
m×m
.
For easy analysis, we choose α
11
= · · · = α
mm
. By exchanging
x
ni
in x with ϑ s
i
+ ζ
i
(i = 1, 2, . . . , m), we define
¯
x
i
= [x
T
1
,
x
T
2
, . . . , x
T
n−1
, x
n1
, x
n2
, . . . , ϑs
i
+ ζ
i
, . . . , x
nm
]
T
∈ R
nm
where ζ
i
=
y
(n−1)
di
− ξ
i
with ξ
i
= λ
i1
e
(n−2)
i
+ · · · + λ
i,n−1
e
i
. ϑ is a scalar and
independent of
¯
x
i
. We can choose suitable B
d
(x) and α, such that
b
dii
α
ii
> 0.
Because B
ϑ
in (18) depends on time t, the time derivative of
V
1
includes the differentiation of matrix B
ϑ
with regard to time
t. To facilitate computation of its derivative, according to Gentle
(2007), we introduce a matrix operator for derivative operation
of matrix-value function with respect to time t, i.e., for a time-
dependent matrix A ∈ R
m×n
and a vector b(t) ∈ R
l
, a matrix
operator M
∂
(A, b) ∈ R
m×n
is defined with the entry of its ith row
and jth column being M
∂ij
(A, b) =
∂A
ij
∂b
T
˙
b with i = 1, 2, . . . , m and
j = 1, 2, . . . , n.
Differentiating (18) with respect to t gives
˙
V
1
= 2S
T
B
ϑ
˙
S + S
T
M
∂
(B
ϑ
, S)S
+ S
T
M
∂
(B
ϑ
, x)S + S
T
M
∂
(B
ϑ
, ζ )S, (20)
where ζ = [ζ
1
, ζ
2
, . . . , ζ
m
]
T
, M
∂
(B
ϑ
, S) ∈ R
m×m
, M
∂
(B
ϑ
, x) ∈
R
m×m
and M
∂
(B
ϑ
, ζ ) ∈ R
m×m
are given below
M
∂
(B
ϑ
, S) = diag
1
0
ϑ
∂B
αii
∂s
i
˙
s
i
dϑ
, (21)
M
∂
(B
ϑ
, x) = diag
1
0
ϑ
nm
j=1,j=i
∂B
αii
∂x
j
˙
x
j
dϑ
, (22)
M
∂
(B
ϑ
, ζ ) = diag
1
0
ϑ
∂B
αii
∂ζ
i
˙
ζ
i
dϑ
, i = 1, . . . , m. (23)
Let σ = ϑs
i
(i = 1, 2, . . . , m), we can obtain
∂B
αii
∂s
i
=
∂B
αii
∂σ
∂σ
∂s
i
= ϑ
∂B
αii
∂σ
, (24)
∂B
αii
∂ϑ
=
∂B
αii
∂σ
∂σ
∂ϑ
=
∂B
αii
∂σ
s
i
. (25)
Noting that ϑ is a scalar and independent of ζ
i
, and the fact
˙
ζ
i
= −ν
i
(i = 1, 2, . . . , m), we have
1
0
ϑ
∂B
αii
∂ζ
i
˙
ζ
i
s
i
dϑ = −ν
i
1
0
ϑ
∂B
αii
∂s
i
dϑ. (26)
It is easy to check that Ψ (z) is well-defined even if S approaches
zero. We design an adaptive control
V = Λ
sgn
u
1
+ u
2
, (27)
u
1
= −β
−1
0
(K
1
+
m
2
α)S
|·|
+
ˆ
W
T
|·|
Ψ
|·|
(z) + Υ
|·|
,
(28)
u
2
= −β
−1
0
S
∥S∥
ρ, (29)
ρ =
γ
1
+ γ
2
∥u
1
∥
1 − γ
2
∥β
−1
0
∥
, if S = 0,
0, if S = 0,
(30)
where Λ
sgn
= diag[sgn(s
i
)]; (∗)
|·|
denotes matrix or vector that
its every element is the absolute value of (∗)’s corresponding
element; ρ is positive when γ
2
< ∥β
−1
0
∥
−1
;
ˆ
W is the estimate
of W ; K
1
is a positive diagonal matrix; Υ will be defined later;
and β
0
= diag[β
10
, . . . , β
m0
] where β
i0
has been defined in
Assumption 2.1.
Consider the following Lyapunov function candidate as
V
2
= V
1
+ V
a
+
m
j=1
V
U
j
(t), (31)
V
a
=
m
i=1
1
2
˜
W
T
i
Ω
−1
i
˜
W
i
, (32)
where
˜
W
i
= W
i
−
ˆ
W
i
.
The adaption law is designed as
˙
ˆ
W
i
= Ω
i
Ψ
i
(z)α
ii
s
i
, (33)
where Ω
i
> 0(i = 1, 2, . . . , m) is a diagonal constant matrix to be
designed.
V
U
j
(t) is introduced to overcome unknown time-delays τ
1
(t),
τ
2
(t), . . . , τ
m
(t) and defined as
V
U
j
(t) =
1
2(1 − ¯τ
max
)
m
k=1
t
t−τ
k
(t)
ϱ
2
jk
(x
k
(τ ))dτ . (34)
From the definition of V
a
, we have
˙
V
a
=
m
i=1
˜
W
T
i
Ω
−1
i
˙
˜
W
i
. (35)
The time derivative of V
U
j
(t) is
˙
V
U
j
(t) =
1
2(1 − ¯τ
max
)
m
k=1
ϱ
2
jk
(x
k
(t))
− ϱ
2
jk
(x
k
(t − τ
k
(t)))(1 − ˙τ
k
(t))
. (36)
Thus, the time derivative of V
2
is
˙
V
2
=
˙
V
1
+
˙
V
a
+
m
j=1
˙
V
U
j
(t)
≤ S
T
α(−K
1
S − Υ ) +
1
2(1 − ¯τ
max
)
m
j=1
m
k=1
ϱ
2
jk
(x
k
(t)), (37)