
(a) The states response. (b) The control signals response.
Figure 1: Simulation of protocol (2) perturbed by a step at 𝑡=𝑡𝑑.
by (2) are shown in Fig. 1 in the time interval [0, 𝑡𝑑]. Observe
that on this time interval the states converge exponentially to
the average of their initial conditions and the control signals all
asymptotically vanish.
This might no longer be the case if the agents are affected by
load disturbances 𝑑𝑖, viz.
¤𝑥𝑖(𝑡)=𝑢𝑖(𝑡) + 𝑑𝑖(𝑡).(3)
An example of what happens in such situations is also shown
in Fig. 1. At the time instance 𝑡=𝑡𝑑one agent is affected by a
unit step disturbance. As a result, all states cease to agree and
start to diverge when 𝑡 > 𝑡𝑑, whereas the control signals reach
non-zero steady-state values.
The apparent instability of the whole system,manifested in the
unboundedness of the states, can be explained by the well-known
fact that the consensus protocol has a closed-loop eigenvalue at
the origin (Olfati-Saber et al., 2007). Nevertheless, the bound-
edness of the control signals under such conditions is intriguing.
Situations wherein some signals in the closed-loop system are
bounded while some others are not normally indicate unstable
pole-zero cancellations in the feedback loop (Zhou et al., 1996,
Sec. 5.3). However, controller (2) is static and thus has no zeros.
1.2. Contribution
The example above suggests that a deeper inspection of the
internal stability property could offer insight into the behavior
of diffusively-coupled systems. The internal stability of any
feedback interconnection requires the stability of all possible
input / output relations in the system, see (Zhou et al., 1996;
Skogestad and Postlethwaite, 2005). However, to the best of
our knowledge, internal stability has not been explicitly studied
in the context of diffusively-coupled architectures of MASs yet.
In this paper we show that diffusively-coupled systems of
LTI (linear time-invariant) agents might not be internally sta-
bilizable. Loosely speaking, this happens if the agents share
common unstable dynamics, directions counting. This, for ex-
ample, is always the case in a group of homogeneous unstable
agents, like those discussed in §1.1.
When restricting the result to finite-dimensional agents, we
also explain the mechanism behind the shown internal instabil-
ity. It is caused by unstable cancellations in the cascade of the
aggregate plant and a diffusive controller. Important is that these
cancellations are caused not by controller zeros, but rather by an
intrinsic spatial deficiency of the diffusively-coupled configura-
tion. These cancellations are intrinsic to the diffusive structure
and cannot be affected by controller dynamics. Consequently,
the internal stability of feedback systems utilizing only relative
measurements depends solely on the agent dynamics.
In addition to providing a rigorous analysis of the internal
stability of diffusively-coupled systems, we show how the anal-
ysis is readily applied to common extensions found in the lit-
erature. In particular, we discuss more general symmetrically
coupled multi-agent systems (i.e. not restricted to only diffu-
sive coupling), asymmetric coupling (i.e. MASs over directed
graphs), unstable systems with no closed right-half plane poles,
and MASs over time-varying networks. Numerous examples
are also provided along the way to illustrate the main results.
The paper is organized as follows. The problem is set up in
Section 2 and the main result is presented in Section 3, with sev-
eral generalizations discussed in §3.2. Section 4 addresses the
case of finite-dimensional agents, reformulating the main result
in a more transparent form and revealing the underlying reason
for the reported behavior. Concluding remarks are provided in
Section 5. Two appendices collect definitions and technical re-
sults about coprime factorizations over 𝐻∞and poles and zero
directions of multivariable real-rational transfer functions.
Notation
The sets of integer, real, and complex numbers are ℤ,ℝ,
and ℂ, respectively, with subsets ℕ𝜈≔{𝑖∈ℤ|1≤𝑖≤𝜈},
ℂ0≔{𝑠∈ℂ|Re 𝑠 > 0}, and ¯
ℂ0≔{𝑠∈ℂ|Re 𝑠≥0}. By 𝐼𝜈
and 𝟙𝜈we denote the 𝜈×𝜈identity matrix and 𝜈-dimensional
vector of ones, respectively. When the dimension is immaterial
or clear from context, we use 𝐼and 𝟙. The complex-conjugate
transpose of a matrix 𝐴is denoted by 𝐴⊤, the set of all its
eigenvalues by spec(𝐴), and its minimal singular value by 𝜎(𝐴).
The notation diag{𝐴𝑖}stands for a block-diagonal matrix with
diagonal elements 𝐴𝑖. The image (range) and kernel (null)
spaces of a matrix 𝐴are notated Im 𝐴and ker 𝐴, respectively.
Given two matrices 𝐴and 𝐵,𝐴⊗𝐵denotes their Kronecker
product.
By the stability of a system 𝐺we understand its 𝐿2-stability.
It is known (Curtain and Zwart, 2020, §A.6.3) that a 𝑝×𝑚LTI
system is causal and stable iff its transfer function 𝐺(𝑠)be-
longs to 𝐻𝑝×𝑚
∞, which is the space of holomorphic and bounded
functions ℂ0→ℂ𝑝×𝑚(we write 𝐻∞when the dimensions are
clear). Given a real-rational transfer function 𝐺(𝑠), its McMil-
lan degree is denoted by deg(𝐺). By nrank 𝐺(𝑠)we understand
the normal rank of a function 𝐺(𝑠).
A digraph G=(V,E) consists of a vertex set Vand an
edge set E ⊂ V × V, see (Godsil and Royle, 2001) for more
details. The (oriented) incidence matrix of Gis denoted by 𝐸G
or simply 𝐸when the association with a concrete graph is clear.
It is a |V | × |E| matrix, whose (𝑖, 𝑗)entry is
[𝐸G]𝑖 𝑗 =
1 if vertex 𝑖is the head of edge 𝑗
−1 if vertex 𝑖is the tail of edge 𝑗
0 if vertex 𝑖does not belong to edge 𝑗
.
2