Nonsingular Fixed Time Consensus Tracking For Second Order Multi Agent Networks
Nonsingular Fixed Time Consensus Tracking For Second Order Multi Agent Networks
net/publication/273480475
CITATIONS READS
678 1,685
1 author:
Zongyu Zuo
Beihang University (BUAA)
170 PUBLICATIONS 7,701 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Fixed-time stabilization control with applications to multi-agent systems and multi-UAV systems View project
All content following this page was uploaded by Zongyu Zuo on 08 April 2018.
Automatica
journal homepage: www.elsevier.com/locate/automatica
Technical communique
by this, this paper extends the idea in Zuo and Tie (2014a,b) where x0 ∈ R, v0 ∈ R and u0 ∈ R denote the position, velocity
and constructs global well-defined fixed-time consensus protocols and input of the reference system. We introduce a virtual leader
for second-order multi-agent networks with directed information π0 with the states ξ0 = [x0 , v0 ]T for the multi-agent system in
flow. The key contributions include: (i) a guaranteed settling time (6). Note, however, that ξ0 is available not to all agents but to only
independent of initial conditions is obtained, i.e., fixed-time con- a portion of agents. Here, define a nonnegative diagonal matrix
sensus; (ii) a continuous sinusoid function is introduced into the B = diag{b1 , b2 , . . . , bN } to indicate the accessibility of ξ0 by the
protocols to eliminate the singularity. To the author’s knowledge, agents, where bi = 1 if ξ0 is accessible by the ith agent, and bi = 0
no results on fixed-time consensus tracking with assignable set- otherwise. The directed graph incorporating π0 into G is denoted
tling time for second-order multi-agent networks are available till by Ge .
now.
Assumption 3. Ge has a spanning tree with π0 being its root
2. Preliminaries vertex, i.e., B ̸= 0.
2.2. Some lemmas To clarify the core idea, we first consider the fixed-time control
of a single system defined by
Lemma 1 (Zuo & Tie, 2014a). Let ξ1 , ξ2 , . . . , ξN ≥ 0. Then ż1 (t ) = z2 (t ), ż2 (t ) = u(t ) (9)
where z = [z1 , z2 ]T ∈ R2 denotes the system state vector, u ∈ R
p
N
N
ξip ≥ ξi if 0 < p ≤ 1 (1) the control input.
i =1 i=1 Let mk , nk , pk , qk be positive odd integers satisfying mk > nk ,
p p1 < q1 < 2p1 , p2 < q2 and m1 /n1 − p1 /q1 > 1, and αk , βk be
N N
p
positive constants, where k = 1, 2. Define a C 1 function µτ (·) :
ξ ≥N 1−p
ξi if 1 < p < ∞. (2)
i [0, +∞) → [0, 1] as
i =1 i =1
π x
sin · if x ≤ τ
µτ ( x ) = 2 τ (10)
Lemma 2 (Zuo & Tie, 2014b). Consider a scalar system 1 otherwise
m
ẏ = −α y n − β y q ,
p
y(0) = y0 (3) where τ is a positive constant.
To circumvent the singularity problem (Feng et al., 2013), a new
where α > 0, β > 0, m, n, p, q are positive odd integers satisfying sliding surface is proposed:
m > n and p < q. The equilibrium of (3) is globally fixed-time stable q1
with settling time T bounded by s = z1 + [κ(z1 ) · z2 ] p1 (11)
z2 = −α1 z1 1 − β1 z1 1 .
n q
Theorem 5. Consider the second-order system in (9) with the control 4. Nonsingular fixed-time consensus tracking
feedback law designed as (13). Then the state z = [z1 , z2 ]T is globally
Let ep = [e1 , e2 , . . . , eN ]T and ev = [ev1 , ev2 , . . . , evN ]T be the po-
p p p
fixed-time stable and the settling time estimate is derived as
sition and velocity disagreement vectors, respectively, with their
T < Tmax := T1 + T2 + ϵ(τ ) (14) p
elements defined by ei = j∈Ni aij (xi − xj ) + bi (xi − x0 ) and
where Tk , nk /[αk (mk − nk )] + qk /[βk (qk − pk )], (k = 1, 2), and evi = (v v ) (v − v0 ), taking time derivatives of
j∈Ni a ij i − j + b i i
ϵ(τ ) denotes a small time margin related to τ . which yields
Proof. Differentiating s in (11) against time yields
ėp (t ) = ev (t ), ėv (t ) = (L + B)u(t ) − B1u0 (t ) (17)
q1
−1
q1 (κ z2 )
m1 p1
where u = [u1 , u2 , . . . , uN ] .
p1
m1 p1 n − q −1
T
ṡ = z2 − α1 − z1 1 1 (κ z2 )2 − κ u .
p1 n1 q1 To solve the fixed-time consensus tracking problem, a new
nonsingular distributed protocol may be derived as
Substituting (13) into the preceding equation obtains
−1
q1 m p2
−1 2
ṡ = −µτ
p1
z2 · α2 s n2 + β2 s q2 . (15) ui = aij + bi udi + aij uj − b i uM (si )
0 sgn (18)
j∈Ni j∈Ni
Observe that µτ (·) > 0 if z2 ̸= 0. For the convenience of the proof, where si denotes a sliding surface, defined by
the state space z ∈ R2 is divided into two different areas, as shown
q1
in Fig. 1,
si = ei + [κi (ei ) · evi ] p1
p p
(19)
q1 q1
p −1 p −1
S1 = (z1 , z2 )|z2 1 ≥τ , S2 = (z1 , z2 )|z2 1 <τ with the virtual control signal udi defined by
(i) When the system states (z1 , z2 ) belong to S1 , the function 1 m1 p1 m1 p
− 1 −1
udi = α1 − (epi ) n1 q1 (κi evi )2
µτ (·) takes value one. Applying Lemma 2, the states (z1 , z2 ) will κi n1 q1
reach the sliding surface s = 0 or enter S2 within fixed-time. q
p1 1− p11 v 2− pq1
p1 − p1
q q1
− κi 1 µτi (evi ) p1
−1
(ii) In S2 , 0 < µτ < 1 when z2 ̸= 0. Applying Lemma 2 for − κi (ei ) 1
where p1 /(q1 − p1 ) > 1. This implies that ϵ(τ ) can be made very
small by choosing sufficiently small τ . For very small z1 including
q1
q1 q1
κi 1 (evi ) p1
−1
0 sgn(si ) + bi u0 .
p
b i uM
the case z1 = 0, z2 → 0 in S2 implies the system trajectory − (23)
p1
308 Z. Zuo / Automatica 54 (2015) 305–309
i=1
where the virtual signal ūdi is given by
N q1
q1 1 m1 p1 m1 p
− q 1 −1
v −1 ūdi α1 p
(ēi ) n1 (κi evi )2
κi (ei ) (
p1
− p1
b i uM − |u0 |)|si | = − 1
0 κi n1 q1
i=1
q q
p1 1− p1 1 q p1 − p11 q1
v 2− p1
µτi (evi ) p1
−1
N
κi 1
(ei ) κi
τ
m2 +n2 p2 +q2 − −
≤ −µm α ( )2
2 si
2n2
+β ( ) 2
2 si
2q2
q1 q1
i =1 q
m2 p2
1− p1
· (evi ) α2 s̄i + β2 s̄i
n2 q2
n2 −m2 m2 +n2 p2 +q2 1 (27)
≤ −µτm α2 N 2n2
(2V2 ) 2n2
+ β2 (2V2 ) 2q2
q /p
where κi 1 1 · (evi )q1 /p1 −1 ≥ 0, µτm , min{µτ1 , . . . , µτN }, Lemma 1 Theorem 9. Consider a multi-agent system (6) and (7) in networks
√
is used in the last inequality. If V2 ̸= 0, then let y2 = 2V2 be the with interaction topology satisfying Assumptions 3 and 4. The
solution to the differential equation nonsingular distributed protocol proposed in (26) achieves fixed-time
n2 −m2 m2 p2 consensus tracking with the settling time bounded by
ẏ2 (t ) = −µτm α2 N 2n2 y2 2 (t ) + β2 y2 2 (t )
n q
ẏ1 (t ) = −α1 N y1 1 (t ) − β1 y1 1 (t ).
2n1 n q
5. Simulation example
p
Similarly as (i), we have = 0 with settling time t1
limt →t1 ei
Consider a group of six agents in a network with the interaction
bounded by t1 < T1 for all i ∈ IN .
Thus, ep converges to zero within the period T < T1 + T2 +ϵ(τ ). graph shown in Fig. 2. The time-varying control input of the
Applying Theorem 4 in Khoo et al. (2009), (8) follows, i.e., the fixed- virtual leader (7) is designed as u0 = − sin(x0 )/(1 + exp(−t ))
time consensus tracking is achieved. verifying uM0 = 1. Consider the initial scenario that x(0) = [−5,
−3, 7, 9, 4, 5]T and v(0) = 0 are set for the group of agents (6) and
To incorporate the interaction topology into the sliding surface, x0 (0) = π /2 and v0 (0) = 0 for the virtual leader (7). The design
we may modify the sliding variables as parameters α1 = β1 = α2 = β2 = 2, m1 = 9, n1 = 5, p1 =
q1
7, q1 = 9, m2 = 11, n2 = 9, p2 = 5, q2 = 7 and τ = 0.1 are
s̄i = ēi + [κi (ēi ) · evi ] p1
p p
(25) set for protocols (18) and (26). It can be verified that the parameter
p
( p p
)+ p
with âij = (aij + aji )/2 for constraints m1 /n1 − p1 /q1 > 1 and p1 /q1 > 1/2 are satisfied. The
where ēi = j∈Ni âij ei − ej bi e i
bounds for settling time T in (21) and (28) are 8.025 s and 24.985 s,
all i, j ∈ IN . Let L̂ = (L + L )/2 be the Laplacian matrix of G(Â)
T
respectively. Since ε2 = 7/9 < 1, a less conservative bounds can
with  , [âij ]n×n . Thus, the corresponding protocol for fixed-time be calculated as 7.889 and 24.425 s.
consensus tracking is then obtained: The simulation results in Fig. 3 show that the consensus tracking
−1
objective is achieved within fixed time. Note that the settling times
ui = aij + bi ūdi + aij uj − bi umax
0 sgn(s̄i ) (26) under the two protocols are about respectively 2 s and 20 s, which
j∈Ni j∈Ni demonstrate the performance claimed in Theorems 8 and 9. It can
Z. Zuo / Automatica 54 (2015) 305–309 309
References