Next Article in Journal
Low-Complexity Channel Estimation in 5G Massive MIMO-OFDM Systems
Next Article in Special Issue
Framework for Onboard Bus Comfort Level Predictions Using the Markov Chain Concept
Previous Article in Journal
Pilot-Based Adaptive Channel Estimation for Underwater Spatial Modulation Technologies
Previous Article in Special Issue
A Dynamic Simulation of the Immune System Response to Inhibit and Eliminate Abnormal Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Cauchy’s Interlacing Theorem and the Stability of a Class of Linear Discrete Aggregation Models Under Eventual Linear Output Feedback Controls

Institute of Research and Development of Processes IIDP, University of the Basque Country, Campus of Leioa, Leioa, PO Box 48940, Bizkaia, Spain
Submission received: 8 May 2019 / Revised: 19 May 2019 / Accepted: 21 May 2019 / Published: 24 May 2019
(This article belongs to the Special Issue Symmetry in Mathematical Analysis and Applications)

Abstract

:
This paper links the celebrated Cauchy’s interlacing theorem of eigenvalues for partitioned updated sequences of Hermitian matrices with stability and convergence problems and results of related sequences of matrices. The results are also applied to sequences of factorizations of semidefinite matrices with their complex conjugates ones to obtain sufficiency-type stability results for the factors in those factorizations. Some extensions are given for parallel characterizations of convergent sequences of matrices. In both cases, the updated information has a Hermitian structure, in particular, a symmetric structure occurs if the involved vector and matrices are complex. These results rely on the relation of stable matrices and convergent matrices (those ones being intuitively stable in a discrete context). An epidemic model involving a clustering structure is discussed in light of the given results. Finally, an application is given for a discrete-time aggregation dynamic system where an aggregated subsystem is incorporated into the whole system at each iteration step. The whole aggregation system and the sequence of aggregated subsystems are assumed to be controlled via linear-output feedback. The characterization of the aggregation dynamic system linked to the updating dynamics through the iteration procedure implies that such a system is, generally, time-varying.

1. Introduction

Stability and convergence properties are very important topics when dealing with both continuous- and discrete-time controlled dynamic systems. In this context, one of the most important design tools is the closed-loop stabilization of control systems via the appropriate incorporation of stabilizing controllers; see, for instance, [1,2,3,4] and references therein. In particular, in [1], and in some references therein, the robust stable adaptive control of tandem of master-slave robotic manipulators using a multi-estimation scheme is discussed. There are several questions of interest in the analysis, such as the fact that the dynamics may be time-varying and imperfectly known, and the fact that a parallel multi-estimation with eventual switching through time is incorporated into the adaptive controller to improve the transient behavior. The speed estimation and stable control of an induction motor based on the use of artificial neural networks is analyzed in [2]. Strategies of decentralized control, including several applications and stabilization tools, are given in [3,4]. In particular, decentralized control is useful when the various subsystems which are integrated in a whole integrated system are located in separate areas, or when the amount of information needed presents difficulties with regards to obtaining completely optimal suitable performance. Thus, the individual controllers associated with the various subsystems get local information about the corresponding subsystems, and eventually some extra partial information about the remaining ones to achieve stabilization, provided that the neglected coupling dynamics are weak enough. Stabilizing decentralized control designs are described in [3] for networked composite systems. Some technical aspects and the results of non-negative matrices of usefulness to describe the properties and behavior of positive dynamic systems, the robustness of matrices against numerical parameterization perturbations of their entries, and the properties of linear dynamic systems are discussed in [5,6,7,8].
This paper focuses on the study of sequences of Hermitian matrices of increasing order which are built via block partition aggregation at each iteration, in such a way that both the current iteration and the next one are Hermitian matrices. The basic mathematical tool is the use of the interlacing Cauchy’s theorem of the matrix eigenvalues of the matrices of the sequence, which orders the sequences of the eigenvalues as the iteration progresses [8]. Our main objective is to adapt the interlacing theorem in order to use it to derive stability or convergence conditions of the sequence of matrices, and to use the results for the stability of a large-scale discrete aggregation-type dynamic system [9,10,11,12,13,14]. The paper is organized as follows. Section 2 is devoted to investigating the properties of boundedness and convergence of the sequences of the determinants and the sequences of eigenvalues as the iteration progresses by aggregation of the updated information while maintaining a Hermitian structure. In the particular case when the matrices describing the problem are real, the updated information has a symmetrical structure. The results are used, in particular, to give stability or anti-stability (in the sense that all the matrix eigenvalues of the matrices of the iterative sequence are unstable) conditions to the matrices used in the standard factorization of Hermitian positive definite matrices. Section 3 extends some of the above results to the convergence of sequences of partitioned Hermitian matrices constructed by aggregation of the updated information. Note that the concept of the convergence of matrices is a discrete counterpart of the matrix stability property in the continuous-time domain, since matrices are stability matrices if all their eigenvalues are in the open complex left-half plane. The basic idea that complex square matrices are convergent if their eigenvalues are within the open unit circle centered at zero is taken into account. An example is discussed concerning a SIR epidemic model with contagions between populations of adjacent clusters in Section 4. Section 5 is devoted to developing an application for the stability of an aggregation discrete-time dynamic controlled system whose order increases by successive incorporations of new subsystems as the iteration index progresses, and whose structure keeps a symmetry. Finally, some conclusions are presented at the end the paper. The relevant mathematical proofs are given in the appendix in order to facilitate a direct reading of the manuscript. The system is assumed to be parameterized by real parameters and controlled by linear output-feedback control laws; it is also assumed that the former whole aggregation system and each new aggregated subsystem at each iteration might eventually be coupled.

Notation and Mathematical Symbols

If M is a square Hermitian matrix, then M 0 denotes that it is positive definite and M _ 0 denotes it is positive semidefinite. Also, M 0 denotes, that it is negative definite.
Z 0 + = { z Z : z 0 } ;   Z + = { z Z : z > 0 } ,
R 0 + = { z R : z 0 } ;   R + = { z R : z > 0 } ,
C 0 + = { z C : Re   z 0 } ;   C + = { z C : Re z > 0 } ,
n ¯ = { 1 , 2 , , n }
I is an identity matrix specified by I n if it denotes the n -th identity matrix, A 0 denotes that the square matrix A is positive definite (positive semidefinite), A 0 denotes that the square matrix A is negative definite (respectively, negative semidefinite), A B , A _ B , A B , A _ B denote, respectively, that A B 0 , A B _ 0 , A B 0 and A B _ 0 , λ min ( M ) and λ max ( M ) denote, respectively, the minimum and maximum eigenvalue of a square real symmetric matrix M , r ( M ) is the spectral radius of any square complex matrix M , s p ( M ) is the set of eigenvalues of the Hermitian matrix M . If such a set is ordered with respect to the partial order relation then the ordered spectrum is denoted by s p ( M ) . The superscripts * and T stand, respectively, for complex conjugates or transposes of any vector or matrix, A B is the Kronecker product of the matrices A , if A C n × m then its vectorization is a vector v e c ( A ) C n × m whose components are all the rows of A written in column in its order and respecting the order of its respective entries, A is the Moore-Penrose pseudoinverse of the matrix A .

2. Technical Results on Partitioned Hermitian Matrices, Cauchy’s Interlacing Theorem and Stability

The subsequent result relies on the conditions for the non-singularity of a partitioned Hermitian matrix of order ( n + 1 ) which is built by aggregation from a principal Hemitian sub-matrix of order n . Mathematical proof is given in Appendix A.
Lemma 1.
Consider the partitioned matrix M = [ M m m d + d ˜ ] C ( n + 1 ) × ( n + 1 ) for any n Z + , where M C n × n is Hermitian, m C n and d , d ˜ R . Then,
(i) 
M is non-singular if and only if det [ M m m d ˜ ] d det M , equivalently, if and only if det [ M m m 0 ] det [ M 0 m d + d ˜ ] .
(ii) 
Assume that M 0 and d > 0 . Then, M 0 if and only if det [ M m m d ˜ ] > det [ M 0 m d ] . If M 0 and d 0 then M _ 0 if det [ M m m d ˜ ] 0 . If M 0 and d 0 then M _ 0 if det [ M m m d ˜ ] 0 .
The subsequent result relies on some conditions which guarantee the boundedness of the determinant and eigenvalues of a recursive sequence of Hermitian matrices which were obtained and supported by Lemma 1 and Cauchy’s interlacing theorem.
Lemma 2.
Consider the recursive sequence of Hermitian matrices { M ( n ) } n = n 0 for a given initial M ( n 0 ) C n 0 × n 0 for some given arbitrary n 0 Z + , where M ( n ) C n × n ; n ( n 0 ) Z + , defined by M ( n + 1 ) = [ M ( n ) m ( n ) m ( n ) d ( n ) + d ˜ ( n ) ] ; n ( n 0 ) Z + and assume that there is a real sequence { ε ( n ) } n = n 0 [ 0 , 1 ) such that 1 k M ( n ) k d ˜ ( n ) m ( n ) * m ( n ) ε ( n ) , equivalently k M ( n ) 1 k d ˜ ( n ) ε ( n ) m ( n ) * m ( n ) ; n ( n 0 ) Z + , where k M ( n ) = | λ min ( M ( n ) ) | K M ( n ) = | λ max ( M ( n ) ) | ; | d ( n ) | [ k d ( n ) , K d ( n ) ] ; | d ˜ ( n ) | [ k d ˜ ( n ) , K d ˜ ( n ) ] ; n ( n 0 ) Z + with K d ( n ) K d ˜ ( n ) and k d ( n ) k d ˜ ( n ) ; n ( n 0 ) Z + . Then, the following properties hold:
(i) 
lim sup n ( | det M ( n + 1 ) | | det M ( n ) | ) 0 , | det M ( n + 1 ) | | det M ( n ) | for any given n ( n 0 ) Z + if d ( n ) , d ˜ ( n ) and m ( n ) satisfy the constraint K d ˜ ( n ) 1 K d ( n ) 1 + ( 2 n + 1 1 ) ε ( n ) ; n ( n 0 ) Z + with K d ( n ) 1 , which becomes | d ˜ ( n ) | 1 | d ( n ) | 1 + ( 2 n + 1 1 ) ε ( n ) ; n ( n 0 ) Z + if | d ( n ) | = K d ( n ) = k d ( n ) 1 and | d ˜ ( n ) | = K d ˜ ( n ) = k d ˜ ( n ) ; n ( n 0 ) Z + .
(ii) 
Assume that s p ( M ( n ) ) = { λ 1 ( n ) , λ 2 ( n ) , , λ n ( n ) } and that K d ˜ ( n ) 1 K d ( n ) 1 + ( 2 n + 1 1 ) ε ( n ) with K d ( n ) 1 , n ( n 0 ) Z + , for some given n ( n 0 ) Z + . Then, the following relations hold:
| λ n + 1 ( n + 1 ) | | i = 1 n λ i ( n ) i = 1 n λ i ( n + 1 ) | = | det M ( n ) | | i = 1 n λ i ( n + 1 ) |
If, furthermore, M ( n ) _ 0 and M ( n + 1 ) _ 0 for some n ( n 0 ) Z + then
λ n ( n ) λ n + 1 ( n + 1 ) i = 1 n λ i ( n ) i = 1 n λ i ( n + 1 ) ;   1 λ n + 1 ( n + 1 ) λ n ( n ) i = 1 n 1 λ i ( n ) i = 1 n λ i ( n + 1 ) ;   i = 2 n + 1 λ i ( n + 1 ) i = 2 n λ i ( n ) λ 1 ( n ) λ 1 ( n + 1 ) 1
(iii) 
Assume that the constraints of Property (ii) hold with M ( n 0 ) _ 0 ; n ( n 0 ) Z + and, furthermore, lim sup n ( d ( n ) + d ˜ ( n ) ) 1 and m ( n ) = min ( o ( M ( n ) ) , o ( | d ˜ ( n ) | ) ) , which is guaranteed if m ( n ) = min ( o ( K M ( n ) ) , o ( K d ˜ ( n ) ) ) . Then, M ( n ) _ 0 and M ( n + 1 ) _ 0 ; n ( n 0 ) Z + , { det M ( n ) } n = n 0 is bounded and the sequence { s p ( M ( n ) ) } n = n 0 is bounded, if det M ( n 0 ) is finite, and then lim sup n ( det M ( n + 1 ) det M ( n ) ) 0 .
Remark 1.
Concerning Lemma 2 (i), we can focus on the following particular cases of interest A:
(a) 
m ( n ) = 0 and | det M ( n + 1 ) | | det M ( n ) | fails for all n ( n ¯ 0 ) Z + and some n ¯ 0 ( n 0 ) Z + . Then, | det M ( n + 1 ) | = | det M ( n ) | | d ( n ) + d ˜ ( n ) | > | det M ( n ) | so that K d ( n ) + K d ˜ ( n ) | d ( n ) + d ˜ ( n ) | > 1 and 1 K d ( n ) 1 + ( 2 n + 1 1 ) ε ( n ) K d ˜ ( n ) > 1 K d ( n ) so that 1 1 + ( 2 n + 1 1 ) ε ( n ) < 1 , a contradiction. Thus, one has | det M ( n + 1 ) | | det M ( n ) | for any n ( n ¯ 0 ) Z + such that m ( n ) = 0 and also lim sup n ( | det M ( n + 1 ) | | det M ( n ) | ) 0 .
(b) 
d ( n ) + d ˜ ( n ) = 0 and | det M ( n + 1 ) | | det M ( n ) | fails for all n ( n ¯ 0 ) Z + and some n ¯ 0 ( n 0 ) Z + . Note from the definition of the recursive sequence { M ( n ) } n = n 0 that
det M ( n + 1 ) = det [ M ( n ) 0 m ( n ) d ( n ) ] + det [ M ( n ) m ( n ) m ( n ) d ˜ ( n ) ] = d ( n ) det M ( n ) + det ( [ M ( n ) 0 0 d ˜ ( n ) ] ( I n + 1 + [ M ( n ) 1 0 0 1 / d ˜ ( n ) ] [ 0 m ( n ) m ( n ) * 0 ] ) ) = d ˜ ( n ) det M ( n ) ( 1 det ( I n + 1 + [ M ( n ) 1 0 0 1 / d ( n ) ] [ 0 m ( n ) m ( n ) * 0 ] ) ) ;   n ( n ¯ 0 ) Z +
Since 1 k M ( n ) k d ˜ ( n ) m ( n ) * m ( n ) ε ( n ) and { ε ( n ) } n = n 0 0 then
lim sup n ( | det M ( n + 1 ) / det M ( n ) | | d ˜ ( n ) | | 1 det ( I n + 1 + [ 0 M ( n ) 1 m ( n ) ( 1 / d ( n ) ) m ( n ) * 0 ] )    | ) = lim sup n | det M ( n + 1 ) / det M ( n ) | = 0
and lim sup n ( | det M ( n + 1 ) | | det M ( n ) | ) 0 since, otherwise, if lim sup n ( | det M ( n + 1 ) | | det M ( n ) | ) > 0 then lim sup n | det M ( n + 1 ) / det M ( n ) | ( , 1 ) ( 1 , + ) , a contradiction to lim sup n | det M ( n + 1 ) / det M ( n ) | = 0 . Note that if d ˜ ( n ) = 0 then | d ( n ) | 1 under the given constraints so that if { m ( n ) } n = n 0 0 ;
(c) 
n ( n ¯ 0 ) Z + and some n ¯ 0 n 0 then lim sup n ( | det M ( n + 1 ) | | det M ( n ) | ) 0 .
Now, one gets from Lemma 2 [(ii), (iii)] the subsequent dual result concerning the recursion obtained from the inverse of M ( n ) . The use of this result will make it possible to give sufficiency-type conditions regarding the non-singularity of the recursive calculation for any positive integer n , and also as n tends to infinity.
Lemma 3.
For some given arbitrary n 0 Z + and all n ( n 0 ) Z + , define:
M ¯ ( n ) = M ( n ) 1 ( I n + m ( n ) ( d ( n ) + d ˜ ( n ) m ( n ) M ( n ) 1 m ( n ) ) 1 m ( n ) M ( n ) 1 )
m ¯ ( n ) = M ( n ) 1 m ( n ) ( d ( n ) + d ˜ ( n ) m ( n ) M ( n ) 1 m ( n ) ) 1
d ¯ ( n ) + d ˜ ¯ ( n ) = ( d ( n ) + d ˜ ( n ) m ( n ) M ( n ) 1 m ( n ) ) 1
and assume that:
(1) 
there is a real sequence { ε ¯ ( n ) } n = n 0 [ 0 , 1 ) such that 1 k ¯ M ( n ) ( n ) k ¯ d ( n ) m ¯ ( n ) * m ¯ ( n ) ε ¯ ( n ) ; n ( n 0 ) Z + , where
k ¯ M ( n ) = | λ min ( M ¯ ( n ) ) | K ¯ M ( n ) = | λ max ( M ¯ ( n ) ) | ;   d ¯ ( n ) [ k ¯ d ( n ) , K ¯ d ( n ) ] ; d ˜ ¯ ( n ) [ k ¯ d ˜ ( n ) , K ¯ d ˜ ( n ) ] ;   n ( n 0 ) Z + ,
d ˜ ( n ) m ( n ) M ( n ) 1 m ( n ) d ( n )
(2) 
M ¯ ( n 0 ) _ 0 , lim sup n ( d ¯ ( n ) + d ˜ ¯ ( n ) ) 1 and m ¯ ( n ) = min ( o ( M ¯ ( n ) ) , o ( d ˜ ¯ ( n ) ) ) , which is guaranteed if m ¯ ( n ) = min ( o ( K ¯ M ( n ) ) , o ( K ¯ d ˜ ( n ) ) ) .
Then, M ¯ ( n ) _ 0 and M ¯ ( n + 1 ) _ 0 ; n ( n 0 ) Z + , { det M ¯ ( n ) } n = n 0 is bounded, the sequence { s p ( M ¯ ( n ) ) } n = n 0 is bounded, if det M ¯ ( n 0 ) is finite, and then lim sup n ( det | M ¯ ( n + 1 ) | det | M ¯ ( n ) | ) 0 .
One gets by combining Lemma 2 and Lemma 3 the two subsequent direct results:
Lemma 4.
Assume that M ( n 0 ) 0 for some given arbitrary n 0 Z + and assume also that the conditions of Lemma 2 (iii) and Lemma 3 hold. Then, M ( n + 1 ) 0 ; n ( n 0 ) Z + .
Lemma 5.
Assume that, for some finite n 0 Z + , A ( n 0 ) C n 0 × n 0 is a stability matrix and construct a sequence { M ( n ) } n = n 0 according to the recursive rule:
M ( n + 1 ) = A ( n + 1 ) * A ( n + 1 ) = [ M ( n ) m ( n ) m ( n ) d ( n ) + d ˜ ( n ) ] = [ A ( n ) * A ( n ) m ( n ) m ( n ) d ( n ) + d ˜ ( n ) ] ;   n n 0
with initial condition M ( n 0 ) = A ( n 0 ) * A ( n 0 ) 0 . Assume also that { M ( n ) } n = n 0 and the sequence of its inverses satisfy the constraints of Lemma 2 [(ii),(iii)] and Lemma 3.
Then, { A ( n ) } n = n 0 is a sequence of stability matrices.
The above result can be directly extended for the case when A ( n 0 ) C n 0 × n 0 is antistable, that is, when all its eigenvalues have positive real parts and M ( n 0 ) 0 . Then, by using similar arguments, as in the proof of Lemma 5 based on the continuity of the matrix eigenvalues with respect to its entries and supported by Lemmas 2,3, according to Cauchy’s interlacing theorem, one concludes that { A ( n ) } n = n 0 consists of antistable members.
Lemma 6.
Lemma 5 holds “mutatis-mutandis” if A ( n 0 ) C n 0 × n 0 is antistable.

3. Some Extended Results Related to Sequences of Convergent Matrices

In order to be able to adapt the above results to discrete dynamic systems, the well-known result that that the stability domain of a convergent matrix (i.e., a “stable” discrete matrix) is the open unit circle of the complex plane centered at zero has to be taken into account. Note that, in particular, A C n × n is convergent if s p ( A ) C 1 = { z C : | z | < 1 } so that { A m } 0 0 as m . It turns out that convergent matrices describe the stability property in the discrete sense. In other words, the solution of the discrete difference vector equation z m + 1 = A z m , where A C n × n , converges to 0 R n for any given z 0 C n if and only if A C n × n is convergent. The relevant results of Section 2 can be extended to this situation as follows, provided that A C n × n is also Hermitian. Consider the following cases:
  • Case a: s p ( A ) C 1 + = { z C : 0 < z < 1 } C 1 so that A ( C n × n ) 0 . Then, it is convergent if and only if ( I n A ) 0 . The proof is direct since if ( I n A ) 0 then x * A x < x * x for any x ( 0 ) C n . Thus, by taking any λ ( 0 ) s p ( A ) of eigenvector x C n , one determines that λ < 1 if λ 0 and λ = 0 directly fulfills the constraint. This proves the sufficiency part. The “only if part” follows, since if ( I n A ) 0 fails, there is λ ( 0 ) s p ( A ) of eigenvector x C n such that x * A x x * x then λ 1 and A is not convergent.
  • Case b: s p ( A ) C 1 = { z C : 1 < z < 0 } C 1 so that A ( C n × n ) 0 . Then, it is convergent if and only if ( I n + A ) 0 . The proof is direct, since if ( I n + A ) 0 , then x * A x > x * x for any x ( 0 ) C n . Thus, by taking any λ ( 0 ) s p ( A ) of eigenvector x C n , one determines that 0 > λ > 1 . The remainder of the proof follows Case a closely.
  • Case c: s p ( A ) C 1 so that A 2 ( C n × n ) 0 so that s p ( A 2 ) C 1 + . Then, it is convergent if and only if ( I n A 2 ) 0 according to Case a by replacing A A 2 . Note that Case c is included Case a and Case b.
Now, for Case a, replace M ( n + 1 ) , defined in Lemma 2, by I n + 1 M ( n + 1 ) = [ I n M ( n ) m ( n ) m ( n ) 1 ( d ( n ) + d ˜ ( n ) ) ] and it has to be guaranteed that if M ( n ) is Hermitian, then ( I n M ( n ) ) is also Hermitan, and ( I n M ( n 0 ) ) 0 for some n 0 Z + then ( I n + 1 M ( n + 1 ) ) 0 ; n n 0 .
For Case b, replace M ( n + 1 ) ( I n + 1 + M ( n + 1 ) ) = [ I n + M ( n ) m ( n ) m ( n ) 1 + d ( n ) + d ˜ ( n ) ] and it has to be guaranteed that if M ( n ) is Hermitian, then ( I n + M ( n ) ) is also Hermitan, and ( I n + M ( n 0 ) ) 0 for some n 0 Z + then ( I n + 1 + M ( n + 1 ) ) 0 ; n n 0 .
For Case c, note that ( I n M ( n ) 2 ) = ( I n + M ( n ) ) ( I n M ( n ) ) so that
( I n M ( n ) 2 ) = [ I n + M ( n ) m ( n ) m ( n ) 1 + d ( n ) + d ˜ ( n ) ] [ I n M ( n ) m ( n ) m ( n ) 1 ( d ( n ) + d ˜ ( n ) ) ]
then, replace
M ( n + 1 ) ( I n + 1 M ( n + 1 ) 2 ) = [ I n M ( n ) 2 ( ( d ( n ) + d ˜ ( n ) ) I n + M ( n ) ) m ( n ) m ( n ) ( ( d ( n ) + d ˜ ( n ) ) I n + M ( n ) * ) 1 ( d ( n ) + d ˜ ( n ) ) 2 m ( n ) m ( n ) ]
and it has to be guaranteed that if M ( n ) 2 is Hermitian and ( I n M ( n 0 ) 2 ) 0 for some n 0 Z + then ( I n + 1 M ( n + 1 ) 2 ) 0 ; n n 0 . Since A 0 is Hermitian, it is of the form A = E * E for some full rank n -matrix E . Then, A 2 = ( E * E ) 2 = E * E E * E = A * A . If
A 0 then ( A ) 0 and ( A ) 2 = ( F * F ) 2 = F * F F * F = ( A * ) ( A ) = A * A for some full rank n -matrix F . Then, Cases a and b can be dealt with using Case c by replacing A 2 A * A .
By taking advantage from the fact that a complex square matrix A is convergent (i.e., stable in the discrete sense) if and only if the Hermitian matrix A * A is convergent, we now build a sequence { M ( n ) } n = 0 of Hermitian matrices as follows, in order to discuss the convergence of its members, provided that M ( n 0 ) is convergent for some given n 0 Z 0 + or, with no loss in generality, provided that M ( 0 ) is convergent. Then,
M ( n + 1 ) = [ M ( n ) λ ( n ) m ( n ) λ ( n ) m ( n ) * δ ( n ) ] = [ M ( n ) 0 0 δ ( n ) ] + [ 0 λ ( n ) m ( n ) λ ( n ) m ( n ) * 0 ] ;   n Z 0 +
with { λ ( n ) } n = 0 [ 0 , 1 ) .Then,
| λ max [ M ( n ) λ ( n + 1 ) m ( n ) λ ( n ) m ( n ) * δ ( n ) ] | | λ max [ M ( n ) 0 0 δ ( n ) ] | + | λ max [ 0 λ ( n ) m ( n ) λ ( n ) m ( n ) * 0 ] | = max ( | λ max ( M ( n ) ) | , | δ ( n ) | ) + λ ( n ) m ( n ) * m ( n ) max ( 1 ε ( n ) , | δ ( n ) | ) + λ ( n ) m ( n ) * m ( n ) < 1 ε ( n ) ;   n Z 0 +
which holds if
m ( n ) * m ( n ) < 1 λ ( n ) 2 ( 1 ε ( n ) max ( 1 ε ( n ) , | δ ( n ) | ) ) 2 ;   n Z 0 +
or, m ( n ) 2 < 1 λ ( n ) 2 ( 1 ε ( n ) max ( 1 ε ( n ) , | δ ( n ) | ) ) ; n Z 0 + , provided that ε ( n + 1 ) < ε ( n ) ; n Z 0 + , that is { ε ( n ) } n = 0 [ 0 , 1 ) is strictly decreasing, so { ε ( n ) } n = 0 0 , and | δ ( n ) | < 1 ε ( n + 1 ) ; n Z 0 + .
Now, assume that the iterations to build { M ( n + 1 ) } n = 0 do not add a new row and column to obtain M ( n + 1 ) from M ( n ) via the contribution of the members of an updating sequence { M ¯ ( n ) } n = 0 ; n Z 0 + but a set of the, in general. Then, one may get that:
M ( n + 1 ) = [ M ( n ) λ ( n ) M ¯ ( n ) λ ( n ) M ¯ ( n ) * Δ ( n ) ] = [ M ( n ) 0 0 Δ ( n ) ] + [ 0 λ ( n ) M ¯ ( n ) λ ( n ) M ¯ ( n ) * 0 ] ;   n Z 0 +
so that
| λ max M ( n + 1 ) | = | λ max [ M ( n ) λ ( n ) M ¯ ( n ) λ ( n ) M ¯ ( n ) * Δ ( n ) ] | | λ max [ M ( n ) 0 0 Δ ( n ) ] | + | λ max [ 0 λ ( n ) M ¯ ( n ) λ ( n ) M ¯ ( n ) * 0 ] | = max ( | λ max ( M ( n ) ) | , | λ max ( Δ ( n ) ) | ) + λ ( n ) λ max 1 / 2 ( M ¯ ( n ) * M ¯ ( n ) ) max ( 1 ε ( n ) , | λ max ( Δ ( n ) ) | ) + λ ( n ) λ max 1 / 2 ( M ¯ ( n ) * M ¯ ( n ) ) < 1 ε ( n + 1 ) ; n Z 0 +
which holds by complete induction if M ( 0 ) is convergent and
M ¯ ( n ) 2 2 = λ max ( M ¯ ( n ) * M ¯ ( n ) ) < 1 λ ( n ) 2 ( 1 ε ( n + 1 ) max ( 1 ε ( n ) , | λ max ( Δ ( n ) ) | ) ) 2 ; n Z 0 +
or, M ¯ ( n ) 2 < 1 λ ( n ) ( 1 ε ( n + 1 ) max ( 1 ε ( n ) , | λ max ( Δ ( n ) ) | ) ) ; n Z 0 + , provided that ε ( n + 1 ) < ε ( n ) , that is { ε ( n ) } n = 0 [ 0 , 1 ) is strictly decreasing, so { ε ( n ) } n = 0 0 , and | λ max ( Δ ( n ) ) | < 1 ε ( n + 1 ) ; n Z 0 + . This implies that { M ( n ) } n = 0 is convergent.
In the particular case that for some λ [ 0 , 1 ) , λ ( n ) = λ n ; n Z 0 + , such a λ is a forgetting factor of the iteration.
We now consider the matrix factorization M ( n ) = A ( n ) * A ( n ) ; n Z 0 + . By construction, M ( n ) is Hemitian (then square), even if A ( n ) is not square; n Z 0 + . In the case when A ( n ) is not square, and since its order strictly increases as n increases, it is possible to consider the convergence of the sequence { A ( n ) } n = 0 (without invoking the values of its eigenvalues) as the following property { A ( n ) } n = 0 is asymptotically convergent if lim m { | A ( n + ξ ) A ( n ) | m } n = 0 = 0 for any given ξ Z + . { A ( n ) } n = 0 is convergent if lim m { A ( n ) m } n = 0 = 0 . The following related results are direct of simple proofs given in Appendix A:
Lemma 7.
If { A ( n ) } n = 0 is convergent then it is asymptotically convergent. The inverse is, in general, not true.
Lemma 8.
Assume that M ( n ) = A ( n ) * A ( n ) ; n Z 0 + is a complex square matrix of any arbitrary order. Then:
(i) 
If A ( n ) 2 < 1 ; n Z 0 + then { A ( n ) } n = 0 and { M ( n ) } n = 0 are convergent sequences.
(ii) 
If M ( n ) 2 < 1 ; n Z 0 + then { M ( n ) } n = 0 and { A ( n ) } n = 0 are convergent sequences.
(iii) 
For any n Z 0 + , A ( n ) 2 < 1 if and only if M ( n ) 2 < 1 . { A ( n ) } n = 0 is convergent if and only if { M ( n ) } n = 0 is convergent.

4. Example of SIR-Type Epidemic Models of Inter-Community Clusters

The stability of the equilibrium points of the epidemic models is an interesting topic which is of great relevance to healthcare management. See, for instance, [15,16,17,18,19]. Now, we discuss an epidemic based-model related to stabilization under the given framework of the Cauchy’s interlacing theorem.
Example 1.
Consider the subsequent continuous-time linearized epidemic model with Q community clusters:
S ˙ i + 1 = ν i ( 1 μ i ) I i β i + 1 S i + 1
I ˙ i + 1 = β i + 1 S i + 1 ν i + 1 I i + 1
R ˙ i + 1 = ν i + 1 μ i + 1 I i
for i = 0 , 1 , , Q 1 with S i ( 0 ) = S i 0 0 , I i ( 0 ) = I i 0 0 , R i ( 0 ) = R i 0 0 ; i Q ¯ are the initial conditions of the susceptible, infectious and recovered subpopulations, respectively, and I 0 ( 0 ) = 0 , ν 0 = 1 . In this model, the infectious subpopulation I i of a community i Q ¯ = { 1 , 2 , , Q } may infect the population of the neighboring community ( i + 1 ) . The parameterization is as follows: β ( . ) are the disease transmission rates, ν ( . ) are the removal rates and μ ( . ) are the separation constants which bifurcate the disease rate between the local community and the total community. Note that the assumption μ 0 = 1 implies that the first cluster is not affected by contagions from any other cluster, [15]. A simple analysis of the trajectory solution of the first cluster shows that
S 1 ( t ) = e β 1 t S 10 0   as   t
at an exponential rate, irrespective of the initial conditions, and is definitively bounded,
I 1 ( t ) = e ν 1 t I 10 + 0 t e ν 1 ( t τ ) S 1 ( τ ) d τ = e ν 1 t ( I 10 + S 10 0 t e ( ν 1 β 1 ) τ d τ ) = ( e ν 1 t I 10 + e β 1 t e ν 1 t ν 1 β 1 S 10 ) 0 as   t
if ν 1 β 1 . If ν 1 = β 1 then I 1 ( t ) = e ν 1 t ( I 10 + S 10 t ) 0 as t . In both cases, the convergence is of exponential order, irrespective of the initial conditions, and is definitively bounded, and
R 1 ( t ) = R 10 + ν 1 μ 1 0 t I i ( τ ) d τ = R 10 + ν 1 μ 1 0 t ( e ν 1 τ I 10 + e β 1 τ e ν 1 τ ν 1 β 1 S 10 ) d τ = R 10 + ν 1 μ 1 0 t ( e ν 1 τ I 10 + e β 1 τ e ν 1 τ ν 1 β 1 S 10 ) d τ = R 10 + ν 1 μ 1 I 10 1 e ν 1 t ν 1 + ν 1 μ 1 S 10 ν 1 β 1 ( e β 1 t e ν 1 t ) R 10 + μ 1 I 10 + ν 1 μ 1 S 10 ν 1 β 1 ( e β 1 t e ν 1 t ) as   t
if ν 1 β 1 with solution which is definitively bounded, and, if ν 1 = β 1 then
R 1 ( t ) = R 10 + ν 1 μ 1 0 t e ν 1 τ ( I 10 + S 10 τ ) d τ = R 10 + μ 1 [ ( 1 e ν 1 t ) I 10 + ( 1 ν 1 ( 1 e v 1 t ) t e ν 1 t ) S 10 ] R 10 + μ 1 [ I 10 + S 10 ν 1 ]   as t
with a solution which is definitively bounded. As a result, the total subpopulation at the first cluster is also definitively bounded, and it converges asymptotically to the limit value of the recovered subpopulation. The solution trajectory is also definitively nonnegative since the matrix of dynamics [ β 1 0 0 β 1 ν 1 0 0 ν 1 μ 1 0 ] is a Metzler matrix and the initial conditions are non-negative. Interpretation shows that the total equilibrium subpopulation is that of the disease-free equilibrium which only has a recovered subpopulation. It can be surprising at a first glance to see that the usual nonlinear term β 1 S 1 ( t ) I 1 ( t ) is the susceptible and infectious subpopulations evolutions of the corresponding SIR Kermack-Mcendrick model counterpart is replaced by a linear term. However, for stability purposes, there is no substantial distinct qualitative behavior between both models, since in this case, S 1 ( t ) = e β 1 0 t I 1 ( τ ) d τ S 10 is strictly decreasing for t 0 and I 1 ( t ) = e 0 t ( β 1 S 1 ( τ ) ν 1 ) d τ I 10 is also strictly decreasing for t 0 provided that S 10 < ν 1 β 1 . Now describe the whole model (14) of N clusters in a more compact way through the individual states x i = ( S i , I i , R i ) T ; i N ¯ and associated matrices of self-dynamics for each i N ¯ and coupled dynamics with the respective preceding cluster ( i 1 ) N ¯ :
A i = [ β i 0 0 β i ν i 0 0 ν i μ i 0 ] ,   i N ¯ ;   A i , i 1 = [ 0 ν i 1 ( 1 μ i 1 ) 0 0 0 0 0 0 0 ] , i ( 2 ) N ¯ ;   A 0 , 1 = 0
so that (13) is equivalently described as
x ˙ i ( t ) = A i x i ( t ) + A i , i 1 x i 1 ( t ) ,   x i ( 0 ) = x i 0 ( 0 ) ;   i N ¯
with x 0 ( t ) 0 for t 0 , and compactly, as follows:
x ˙ ( t ) = A x ( t ) ,   x ( 0 ) = x 0
where x ( t ) = ( x 1 T ( t ) , x 2 T ( t ) , , x N T ( t ) ) T and
A = [ A 1 0 0 0 A 21 A 2 0 0 0 A 32 A 3 0 0 0 0 0 A N , N 1 A N ]
Thus, system (15), like (16), (17), can be interpreted as an aggregation model given by starting with the first cluster and successively incorporating the dynamics of the remaining clusters. Now, define the symmetric matrix M = M ( N ) = A T A . Then, define:
M ( 1 ) = A 1 T A 1
M ( 2 ) = [ A 1 T A 21 T 0 A 2 T ] [ A 1 0 A 21 A 2 ] = [ A 1 T A 1 + A 21 T A 21 A 21 T A 2 A 2 T A 21 A 2 T A 2 ] = [ M ( 1 ) A 21 T A 2 A 2 T A 21 A 2 T A 2 ] + M ˜ ( 2 )
M ( 3 ) = [ A 1 T A 21 T 0 0 A 2 T A 32 T 0 0 A 3 T ] [ A 1 0 0 A 21 A 2 0 0 A 32 A 3 ] = [ M ( 2 ) 0 A 32 T A 3 0 A 3 T A 32 A 3 T A 3 ] + M ˜ ( 3 ) M ( N ) = [ M ( N 1 ) m ( N 1 ) m ( N 1 ) T A 3 T A 3 ] + M ˜ ( N )
M ˜ ( 2 ) = [ A 21 T A 21 0 0 0 ] ,   M ˜ ( 3 ) = [ A 21 T A 21 0 0 0 A 32 T A 32 0 0 0 0 ]
M ˜ ( N ) = [ A 21 T A 21 0 0 0 A 32 T A 32 0 0 0 0 A N , N 1 T A N , N 1 0 ]
if N 2 and M ˜ ( 1 ) = 0 . By inspection of (18), one concludes that M ˜ ( n ) = 0 ( n = 2 n sup 2 m ( n ) ) for i = 2 , 3 , , N which concludes that, if { A n 1 , n } n = 2 N 0 as N in such a way that { m ( j ) } j = 2 N 0 as N , for instance, if the convergence is at exponential rate, then lim N n = 2 N sup 2 m ( n ) K m < + . Furthermore, if A 1 is a stability matrix of absolute stability abscissa which are sufficiently larger than K m , then the dynamic system (16),(17) is globally asymptotically stable according to Lemma 5. In particular, note that if there is any pair of stable complex conjugate eigenvalues s 1 , 2 = μ ± i ν ( μ < 0 , ν > 0 ) for the first cluster, then there is a submatrix of M ( 1 ) ,
M s ( 1 ) = A 1 s T A 1 s = [ μ ν ν μ ] [ μ ν ν μ ] = [ μ 2 + ν 2 0 0 μ 2 + ν 2 ] = [ λ 1 ( 1 ) 0 0 λ 2 ( 1 ) ]
in the real canonical form. Since μ ( 1 ) = μ < 0 then λ 2 ( 1 ) = λ 2 ( 1 ) = μ 2 + ν 2 > 0 . The maximum and minimum corresponding eigenvalues of M s ( 2 ) are no less than λ 2 ( 1 ) and no larger than λ 1 ( 1 ) , respectively, from Cauchy’s interlacing theorem. Since the eigenvalues are continuous functions of the matrix entries, and since M s ( n ) is positive definite any critical eigenvalue of a member M s ( j ) of the sequence, { M s ( n ) } n = 2 implies a lot of stability of the corresponding A j s . This is avoided if lim N n = 2 N sup 2 m ( n ) K m < + , implying also that { m ( n ) } 2 N 0 as N , and the sequence of separation constants { μ n } n = 2 N 1 as N if K m is small enough related to μ . The physical interpretation relies on the fact that the contagion link from a cluster to the next one is weakened sufficiently quickly as the cluster index increases, due to the fact of the numbers of the infected subpopulations are rapidly decreasing as the cluster index increases at a sufficiently large rate.

5. Dynamic Linear Discrete Aggregation Model with Output Delay and Linear Feedback Control

In this section, the convergence results of Section 3 are applied to a dynamic discrete system which is built by the aggregation of discrete dynamic subsystems subject to linear output feedback control. Since we are dealing with a physical system, it turns out that the formalism of Section 2 can be developed by invoking conditions related to real symmetric systems, rather than to complex Hermitian ones, when necessary. It would suffice to describe the state by expressing the matrix of dynamics in the real canonical form and to transform the control and output matrices by the appropriate similarity matrix. The necessary mathematical proof is given in Appendix A.
Consider the aggregation linear discrete dynamic system subject to r point delays under linear output-feedback:
x 0 ( n + 1 ) = A 0 ( n ) x ( n ) + A ^ 0 ( n ) x ^ ( n ) + j = 1 r B j ( n ) y ( n j ) + B 0 ( n ) u ( n )
y ( n ) = C ( n ) x ( n )
u ( n ) = j = 0 r n K j ( n ) y ( n j ) = j = 0 r n K j ( n ) C x ( n j )
n Z 0 + , with initial conditions x 0 ( 0 ) = x 0 , where { n i } i = 0 is a sequence of positive integer numbers, x ( n ) R i = 0 n n i is the “a priori” vector state at the n-th iteration, x ^ ( n ) R n n is the aggregated “a priori” new substate at the n -th iteration (that is basically, the new information needed to update the state vector and its dimension) and x 0 ( n + 1 ) R i = 0 n n i is the “a priori” whole state at the ( n + 1 ) -th iteration. Also, x 0 ( n ) R i = 0 n 1 n i , u ( n ) R i = 0 n m i and y ( n ) R i = 0 n p i are, respectively, the “a priori” input and measurable output vectors at the n -th iteration and r Z + is a sequence of delays influencing the global dynamics. The sequences of matrices of dynamics { A 0 ( n ) } n = 0 and { A ^ 0 ( n ) } n = 0 , control { B 0 ( n ) } n = 0 , output-state coupling { B j ( n ) } n = 0 for j r ¯ are of members A 0 ( n ) R ( i = 0 n n i ) × ( i = 0 n n i ) , A ^ ( n ) R ( i = 0 n n i ) × n n , and B 0 ( n ) R ( i = 0 n n i ) × ( i = 0 n m i ) and B j ( n ) R ( i = 0 n n i ) × ( i = 0 n j p i ) for j r ¯ and the output matrix C ( n ) R ( i = 0 n p i ) × ( i = 0 n n i ) . The sequences of matrices { K j ( n ) } n = 0 , with K j ( n ) R ( i = 0 n m i ) × ( i = 0 n p i ) for j r ¯ { 0 } , are the output-feedback control gains which generate the control law sequence { u ( n ) } n = 0 .
The dynamics of the new dynamics at the ( n + 1 ) -th iteration aggregated to the former global aggregation system of state x ( n ) obtained at the n -th iteration, are assumed to be described by:
x ^ ( n + 1 ) = A ^ ( n + 1 ) x ( n ) + ( D ^ ( n ) + D ^ ˜ ( n ) ) x ^ ( n ) + j = 1 r ( B ^ j a ( n ) y ( n j ) + B ^ j ( n ) y ^ ( n j ) ) + B ^ 0 ( n ) u ^ ( n )
y ^ ( n ) = C ^ ( n ) x ^ ( n )
u ^ ( n ) = j = 0 r ( K ^ j ( n ) y ^ ( n j ) + K ^ j a ( n ) y ( n j ) ) = j = 0 r ( K ^ j ( n ) C ^ ( n ) x ^ ( n j ) + K ^ j a ( n ) C ( n ) x ( n j ) )
n Z 0 + , where x ^ ( n + 1 ) R n n + 1 is the “a posteriori” state of the aggregated subsystem at the ( n + 1 ) -th iteration whose “a priori” value is x ^ ( n ) R n n , y ^ ( n ) R p n , u ^ ( n ) R m n , A ^ ( n + 1 ) R n n + 1 × ( i = 0 n n i ) , B ^ 0 ( n ) R n n + 1 × m n , B ^ j ( n ) R n n + 1 × p n j , B ^ j a ( n ) R n n + 1 × ( i = 0 n j p i ) for j r ¯ ; C ^ ( n j ) R p n j × n n j , K ^ j ( n ) R m n + 1 × p n j , K ^ j a ( n ) R m n + 1 × ( i = 0 n j p i ) for j r ¯ { 0 } , and D ^ ( n ) , D ^ ˜ ( n ) R n n + 1 × n n and B ^ 0 ( n ) R n n + 1 × m n + 1 , B ^ j ( n ) R n n + 1 × p n j for j r ¯ ; n Z 0 + .
Note that the aggregated subsystem (22)–(24) is coupled to the former global state x ( n ) describing the total system’s dynamics prior to the aggregation action. It can be seen that the coupling terms do not necessary demonstrate infinite memory requirements as n tends to infinity, since the matrices A ( n + 1 ) , B ^ j a ( n ) and K ^ j a ( n ) can contain nonzero columns associated with the most recent state/output data related to the previous aggregation system; see, for instance, [12]. Note also that, due to the coupling between the a priori whole state at the n -th iterations with the a priori new aggregated substate, it can happen that the a posteriori vector after the new aggregated substate has a higher dimension than its a priori version. The various dynamics, control and output matrices have the appropriate orders.
After incorporating the control law, we can write this whole system of extended states x ( n ) = ( x 0 ( n ) T , x ^ ( n ) T ) T R i = 0 n n i ; n Z 0 + in a compact way:
x ( n + 1 ) = [ x 0 ( n + 1 ) x ^ ( n + 1 ) ] = [ A 0 ( n ) + B 0 ( n ) K 0 ( n ) C ( n ) A ^ 0 ( n ) A ^ ( n + 1 ) + B ^ 0 ( n ) K ^ 0 a ( n ) C ( n ) D ^ ( n ) + D ^ ˜ ( n ) + B ^ 0 ( n ) K ^ 0 ( n ) C ^ ( n ) ] [ x ( n ) x ^ ( n ) ] + j = 1 r [ ( B j ( n ) + B 0 ( n ) K j ( n ) ) C ( n ) 0 ( B ^ j a ( n ) + B ^ 0 ( n ) K ^ j a ( n ) ) C ( n j ) ( B ^ j ( n ) + B ^ 0 ( n ) K ^ j ( n ) ) C ^ ( n j ) ] x ( n j ) ;   n Z 0 +
so that x ( n ) , x 0 ( n + 1 ) R i = 0 n n i and x ^ ( n ) R n n and x ^ ( n + 1 ) R n n + 1 imply that x ( n + 1 ) R i = 0 n + 1 n i , x ( n + 1 ) = ( x 0 ( n + 1 ) T , x ^ ( n + 1 ) T ) T R i = 0 n + 1 n i .
In order to construct a state vector which includes delayed dynamics, we now define the modified extended state x ¯ ( n ) defined by x ¯ ( n ) = ( x ( n ) T , x ( n 1 ) T , , x ( n r ) T ) T R j = n r n i = 0 j n i ; n Z 0 + . Thus, one determines from (25) that:
x ¯ ( n + 1 ) = A ¯ ( n ) x ¯ ( n ) ;   n Z 0 +
where
A ¯ ( n ) = [ A 0 ( n ) + B 0 0 ( n ) K 0 ( n ) C ( n ) A ^ 0 ( n ) A ^ ( n + 1 ) + B ^ 0 ( n ) K ^ 0 a ( n ) C ( n ) D ( n ) + D ˜ ( n ) + B ^ 0 ( n ) K ^ 0 ( n ) C ^ ( n ) ( B 1 ( n ) + B 0 ( n ) K 1 ( n ) ) C ( n 1 ) 0 ( B ^ 1 a ( n ) + B ^ 0 ( n ) K ^ 1 a ( n ) ) C ( n 1 ) ( B ^ 1 ( n ) + B ^ 0 ( n ) K ^ 1 ( n ) ) C ^ ( n 1 ) ( B r ( n ) + B 0 ( n ) K r ( n ) ) C ( n r ) 0 ( B ^ r a ( n ) + B ^ 0 ( n ) K ^ r a ( n ) ) C ( n r ) ( B ^ r ( n ) + B ^ 0 ( n ) K ^ r ( n ) ) C ^ ( n r ) I i = 0 n 1 n i I i = 0 n r n i 0 0 ]
n Z 0 + , with A ¯ ( n ) R ( j = n r n + 1 i = 0 j n i ) × ( j = n r n i = 0 j n i ) . Now, consider the symmetric matrices:
M ¯ ( n ) = A ¯ ( n ) T A ¯ ( n ) = [ A ¯ ( n ) T A ¯ ( n ) A ¯ ( n ) T B ¯ ( n ) B ¯ ( n ) T A ¯ ( n ) B ¯ ( n ) T B ¯ ( n ) ]
M ( n ) = [ A 0 ( n + 1 ) T A 0 ( n + 1 ) A 0 ( n + 1 ) T B 0 ( n + 1 ) B 0 ( n + 1 ) T A 0 ( n + 1 ) B 0 ( n + 1 ) T B 0 ( n + 1 ) ] R ( i = 0 n + 1 n i + j = n r n ( i = 0 j n i + n j ) ) × ( i = 0 n + 1 n i + j = n r n ( i = 0 j n i + n j ) ) ;   n Z 0 +
where the relations between the a priori dynamics of the new iteration after the aggregation of a new substate to the whole dynamics with the a posteriori dynamics of the former iteration are given by:
A 0 ( n + 1 ) = [ A ( n )    0 ( i = 0 n + 1 n i ) × n n + 1 ] = [ A 0 ( n ) + B 0 ( n ) K 0 ( n ) C ( n ) A ^ 0 ( n ) A ^ ( n + 1 ) + B ^ 0 ( n ) K ^ 0 a ( n ) C ( n ) D ^ ( n ) + D ^ ˜ ( n ) + B ^ 0 ( n ) K ^ 0 ( n ) C ^ ( n ) 0 0 ] R ( i = 0 n + 1 n i ) × ( i = 0 n + 1 n i )
B 0 ( n + 1 ) = B ( n ) = [ ( B 1 ( n ) + B 0 ( n ) K 1 ( n ) ) C ( n 1 ) 0 ( B ^ 1 a ( n ) + B ^ 0 ( n ) K ^ 1 a ( n ) ) C ( n 1 ) ( B ^ 1 ( n ) + B ^ 0 ( n ) K ^ 1 ( n ) ) C ^ ( n 1 ) ( B j ( n ) + B 0 ( n ) K r ( n ) ) C ( n r ) 0 ( B ^ r a ( n ) + B ^ 0 ( n ) K ^ r a ( n ) ) C ( n r ) ( B ^ j ( n ) + B ^ 0 ( n ) K ^ r ( n ) ) C ^ ( n r ) ] R ( i = 0 n + 1 n i ) × ( j = n r n 1 ( i = 0 j n i + n j ) )
n Z 0 + . which are built in order to complete a square a priori matrix of dynamics of the ( n + 1 ) - the aggregated system which was obtained after the aggregation of the ( n + 1 ) -th subsystem.
The stability of the aggregation dynamic system (19) to (24) under discrete delays is now discussed via the modified extended system (26), subject to (27), which can be obtained via Lemmas 7,8 from the convergence of the symmetric matrix (28), subject to (29),(30). The following result holds:
Theorem 1.
The following properties hold:
(i) 
{ M ( n ) } n = 0 and { A ¯ ( n ) } n = 0 are convergent, and also asymptotically convergent, if and only if lim m A ¯ ( n ) m = 0 ; n Z 0 + .
(ii) 
{ M ( n ) } n = 0 and { A ¯ ( n ) } n = 0 are asymptotically convergent if and only if lim m { | A ¯ ( n + ξ ) A ¯ ( n ) | m } n = 0 = 0 for any given n Z 0 + , ξ Z + .
(iii) 
If { A ¯ ( n ) } n = 0 (and then { M ( n ) } n = 0 ) is convergent, then the state of the modified extended system, (26), converges asymptotically to zero, i.e. x ¯ ( n + m ) 0 and also x ( n + m ) 0 as m for any given initial condition x ( 0 ) and any n Z 0 + so that the aggregation system is globally asymptotically stable.
(iv) 
If { A ¯ ( n ) } n = 0 (and then { M ( n ) } n = 0 ) is asymptotically convergent then x ¯ ( n + m + ξ ) x ¯ ( n + m ) 0 as m for any given n Z 0 + , ξ Z + and any given initial condition x ( 0 ) and also x ( n + m + ξ ) x ( n + m ) 0 as m for any given initial condition x ( 0 ) and any given n Z 0 + , ξ Z + so that the incremental aggregation system is globally asymptotically stable.
(v) 
Assume that M ( 0 ) is convergent and that | λ max ( B ( n ) T B ( n ) ) | < min ( 1 ε ( n + 1 ) , ε ( n ) ε ( n + 1 ) ) ; n Z 0 + for some strictly decreasing real sequence { ε ( n ) } n = 0 [ 0 , 1 ) . Then, M ( n ) 2 < 1 ; n Z 0 + and { M ( n ) } n = 0 and { A ( n ) } n = 0 are convergent sequences.
It is of interest to now discuss how the stability properties of the aggregation system of Theorem 1 can be guaranteed or addressed by the synthesis of the basic controller (21) on the current aggregated system, and how its updated rule (24) can be applied to the new aggregated subsystem to generate the aggregated system for the next iteration step. This discussion invokes conditions to guarantee that the equation of dimensionally compatible real matrices
B K C = A m A
is solvable in K for a given quadruple ( A , B , C , A m ) with A and A m being square, A m being convergent (basically stable in the discrete context) and defining the closed-loop system dynamics after linear output-feedback control u = K y = K C x via the linear stabilizing controller of gain K ; A , B and C are the open-loop dynamics (i.e., the one being got for K = 0 ) and B and C are the control and output matrices. Equation (31) is written in equivalent vector form for the unknown K as follows:
( B C T ) v e c ( K ) = v e c ( A m A )
It turns out that (31) is solvable in K if and only if (32) is solvable in v e c ( K ) , that is, if r a n k ( B C T ) = r a n k [ B C T , v e c ( A m A ) ] according to the Rouché-Froebenius theorem for solvability of linear systems of algebraic equations. Note that if A m satisfies the constraint A m = E A , for some square matrix E of the same order as A , then A m 2 < 1 (so that A m is convergent) if E 2 < 1 / A 2 . In particular, if A m = ρ A with ρ R then A m is convergent if | ρ | < 1 / A 2 . A preliminary technical result concerning the solvability if the concerned algebraic system (31), or equivalently (32), is (either indeterminate or determinate) compatible to be then used follows:
Lemma 9.
Assume that A R n × n , B R m × n and C R p × n . Then, the following properties hold:
(i) 
linear output-feedback controller exists which stabilizes the closed-loop matrix of dynamics A m = E A for some E with E 2 < 1 / A 2 , [5,7], which satisfies the rank constraint:
r a n k ( B C T ) = r a n k [ B C T , ( I A T ) v e c ( E ) v e c ( A ) ]
If (33) holds, then the set of stabilizing linear-output feedback controllers of gains K which solve (32), equivalently (31), which is a compatible algebraic linear system, for A m = E A , are given by
v e c ( K ) = ( B C T ) [ ( I A T ) v e c ( E ) v e c ( A ) ] + [ I ( B C T ) ( B C T ) ] k w
with k w being any arbitrary real vector of the same dimension as v e c ( K ) . Assume that ( B C T ) = ( B C T ) 1 (a necessary condition being min ( m , p ) n ). Then (32) for A m = E A is a compatible determinate, and the unique solution to (33) is
v e c ( K ) = ( B C T ) [ ( I A T ) v e c ( E ) v e c ( A ) ]
If A m = ρ A with | ρ | < 1 / A 2 then (33), (34) and (35) become, in particular,
r a n k ( B C T ) = r a n k [ B C T , ( ρ 1 ) v e c ( A ) ]
v e c ( K ) = ( ρ 1 ) ( B C T ) v e c ( A ) + [ I ( B C T ) ( B C T ) ] k w
and
v e c ( K ) = ( B C T ) ( B C T ) v e c ( A )
(ii) 
Assume that A m = E A and
r a n k [ B C T , ( I A T ) v e c ( E ) v e c ( A ) ] = r a n k ( B C T ) + 1
Then, (32), equivalently (31), is an algebraically incompatible system of equations, and
v e c ( K ) = ( B C T ) [ ( I A T ) v e c ( E ) v e c ( A ) ] ,   i . e .
i.e., Equation (34) for k w = 0 , is the best least-squares approximated solution to (32) in the sense that the corresponding controller gain minimizes the norm error B K C + A A m 2 2 . If (39) holds for any A m of the form A m = E A then there is no solution to (31) in K ; only best approximation solutions exist.
Particular cases of interest which are well-known from basic Control Theory (see e.g., [13]) are:
(1)
p = n , C R n × n is non-singular and ( A , B ) is stabilizable, i.e., any unstable or critically unstable mode of the open-loop dynamics can be closed-loop stabilized under linear state feed-back control. Thus, r a n k [ λ I n A , B ] = n ; λ C with | λ | 1 (discrete form of Popov-Belevitch-Hautus stabilizability test [6,13,14]). Then (31) becomes B K = ( A m A ) C 1 , which is solvable in K , and there is always an output-feedback stabilizing linear controller generating a stabilizing controller of gain K , generating a control u = K y = K C x , such that the closed-loop dynamics is defined by a convergent matrix A m = A + B K C .
(2)
In Case 1, C = I n . Then, the control law is a linear state-feedback control, and a state-feedback stabilizing linear controller generating a control u = K x exists, leading to closed-loop dynamics defined by the convergent matrix A m = A + B K .
Lemma 9 is useful to guarantee the relevant results of Theorem 1 in terms of the controller gains choices under certain algebraic solvability conditions. This feature is addressed in the subsequent result:
Theorem 2.
Assume that:
(1) r a n k ( B 0 ( n ) C ( n ) T ) = r a n k ( B 0 ( n ) C ( n ) T , A f ( n ) A 0 ( n ) ) so that A 0 ( n ) + B 0 ( n ) K 0 ( n ) C ( n ) = A f ( n ) is solvable in K 0 ( n ) for some convergent matrix A f ( n ) of appropriate order; n Z 0 + , r a n k ( B ^ 0 ( n ) C ( n ) T ) = r a n k ( B ^ 0 ( n ) C ( n ) T , A ^ f a ( n + 1 ) A ^ ( n + 1 ) ) so that A ^ ( n + 1 ) + B ^ 0 ( n ) K ^ 0 ( n ) C ( n ) = A ^ f a ( n + 1 ) is solvable in K ^ 0 a ( n ) for some matrix A ^ f a ( n + 1 ) of appropriate order; n Z 0 + , r a n k ( B ^ 0 ( n ) C ^ ( n ) T ) = r a n k ( B ^ 0 ( n ) C ^ ( n ) T , A ^ f ( n + 1 ) D ^ ( n ) D ^ ˜ ( n ) ) so that A ^ f ( n ) = D ^ ( n ) + D ^ ˜ ( n ) + B ^ 0 ( n ) K ^ 0 ( n ) C ^ ( n ) is solvable in K ^ 0 ( n ) for some matrix A ^ f ( n ) of appropriate order; n Z 0 + ,
(2) and that subsequent rank conditions hold:
r a n k ( B 0 ( n ) C ( n i ) T ) = r a n k ( B 0 ( n ) C ( n i ) T , B i ( n ) C ( n i ) ) r a n k ( B ^ 0 ( n ) C ^ ( n i ) T ) = r a n k ( B ^ 0 ( n ) C ^ ( n i ) T , B ^ i ( n ) C ^ ( n i ) ) r a n k ( B ^ 0 ( n ) C ( n i ) T ) = r a n k ( B ^ 0 ( n ) C ( n i ) T , B ^ i a ( n ) C ( n i ) )
i r ¯ so that the following matrix equations are solvable in the delayed controller gains K i ( n ) , K ^ i ( n ) and K ^ i a ( n ) :
B 0 ( n ) K i ( n ) C ( n i ) = B i ( n ) C ( n i ) ; B ^ 0 ( n ) K ^ i ( n ) C ^ ( n i ) = B ^ i ( n ) C ^ ( n i ) ; B ^ 0 ( n ) K ^ i a ( n ) C ( n i ) = B ^ i ( n ) C ( n i ) ; i r ¯ ; n Z 0 + .
Then, the matrix equations
A 0 ( n ) + B 0 ( n ) K 0 ( n ) C ( n ) = A f ( n ) ,   A ^ ( n + 1 ) + B ^ 0 ( n ) K ^ 0 a ( n ) C ( n ) = A ^ f a ( n + 1 ) , A ^ f ( n ) = D ^ ( n ) + D ^ ˜ ( n ) + B ^ 0 ( n ) K ^ 0 ( n ) C ^ ( n ) ,   B 0 ( n ) K i ( n ) C ( n i ) = B i ( n ) C ( n i ) ; B ^ 0 ( n ) K ^ i ( n ) C ^ ( n i ) = B ^ i ( n ) C ^ ( n i ) ;   i r ¯ ;   n Z 0 + )
are solvable in the controller gains K 0 ( n ) , K i ( n ) and K ^ i ( n ) ; i r ¯ ; n Z 0 + leading to the solutions
v e c ( K 0 ( n ) ) = ( B 0 ( n ) C ( n ) T ) v e c ( A f ( n ) A 0 ( n ) ) + ( I ( B 0 ( n ) C ( n ) T ) ( B 0 ( n ) C ( n ) T ) ) v e c ( K 0 v ( n ) ) v e c ( K ^ 0 a ( n ) ) = ( B ^ 0 ( n ) C ( n ) T ) v e c ( A ^ f a ( n + 1 ) A ^ ( n + 1 ) ) + ( I ( B ^ 0 ( n ) C ( n ) T ) ( B ^ 0 ( n ) C ( n ) T ) ) v e c ( K ^ 0 a v ( n ) ) v e c ( K ^ 0 ( n ) ) = ( B ^ 0 ( n ) C ^ ( n ) T ) v e c ( A ^ f ( n + 1 ) D ^ ( n ) D ^ ˜ ( n ) ) + ( I ( B ^ 0 ( n ) C ^ ( n ) T ) ( B ^ 0 ( n ) C ^ ( n ) T ) ) v e c ( K ^ 0 v ( n ) ) v e c ( K i ( n ) ) = ( B 0 ( n ) C ( n i ) T ) v e c ( B i ( n ) C ( n i ) ) + ( I ( B 0 ( n ) C ( n i ) T ) ( B 0 ( n ) C ( n i ) T ) ) v e c ( K i v ( n ) ) v e c ( K ^ i ( n ) ) = ( B ^ 0 ( n ) C ^ ( n i ) T ) v e c ( B ^ i ( n ) C ^ ( n i ) ) + ( I ( B ^ 0 ( n ) C ^ ( n i ) T ) ( B ^ 0 ( n ) C ^ ( n i ) T ) ) v e c ( K ^ i v ( n ) )
v e c ( K ^ i a ( n ) ) = ( B ^ 0 ( n ) C ( n i ) T ) v e c ( B ^ i a ( n ) C ( n i ) ) + ( I ( B ^ 0 ( n ) C ( n i ) T ) ( B ^ 0 ( n ) C ( n i ) T ) ) v e c ( K ^ i a v ( n ) )
with K 0 v ( n ) , K 0 a v ( n ) , K ^ 0 v ( n ) , K i v ( n ) , K ^ i v ( n ) and K ^ i a v ( n ) ; i r ¯ ; n Z 0 + being arbitrary matrices of appropriate orders for the corresponding equation (above) in each case whose equivalent vector expressions are denoted by v e c ( . ) .
It should be pointed out that it can be of interest to apply the results on interlacing Cauchy’s theorem and some of its extensions (see e.g., [20,21,22]) to the stability of aggregation models based on dynamic systems formulated via differential, difference or hybrid differential/difference equations.

6. Conclusions

This paper relies on partitioned Hermitian matrices and Cauchy’s interlacing theorem and the associated stability results. Based on the fact that convergent matrices are a discrete counterpart of stability matrices, the results presented above are then extended to sequences of convergent matrices. Then, an example of a SIR-type epidemic model continuous-time consisting of intercommunity clusters is discussed relative to the previously given stability theoretic results under the proposed framework based on Cauchy’s interlacing theorem, and which may be of interest for healthcare management. Later, a dynamic linear discrete aggregation model is discussed, which involves output delay and linear output feedback, and which can be also be reformulated for linear-state feedback by identifying state and output, that is, by taking the output matrix equal to the identity, provided that the state is available for measurement. The studied aggregation model is built through the successive incorporation of discrete subsystems with particular coupled dynamics. Stabilizing decentralized controllers are proposed and discussed for this type of aggregation model.

Author Contributions

The author contributed by himself the whole manuscript.

Funding

This research has been jointly supported by the Spanish Government and the European Commission with Grants RTI2018-094336-B-I00 (MINECO/FEDER, UE) and DPI2015-64766-R (MINECO/FEDER, UE).

Acknowledgments

The author is grateful to the Spanish Government for Grants RTI2018-094336-B-I00 and DPI2015-64766-R (MINECO/FEDER, UE).

Conflicts of Interest

The author declares that he has no competing interests regarding the publication of this article.

Appendix A. Mathematical Proofs

Proof of Lemma 1.
From the given assumptions, M is Hermitian, and note that
det M = det [ M 0 m d ] + det [ M m m d ˜ ] = det [ M 0 m d ˜ ] + det [ M m m d ] = det [ M 0 m d + d ˜ ] + det [ M m m 0 ] .
Thus, Properties (i) to (iii) follow directly from the above relations by noting, furthermore, that det [ M 0 m c ] = c det M . □
Proof of Lemma 2.
First, note that
k d ( n ) k d ˜ ( n ) | | d ( n ) | | d ˜ ( n ) | | | d ( n ) + d ˜ ( n ) | K d ( n ) + K d ˜ ( n ) ;   n ( n 0 ) Z + .
Thus,
det M ( n + 1 ) = det [ M ( n ) 0 m ( n ) d ( n ) ] + det [ M ( n ) m ( n ) m ( n ) d ˜ ( n ) ] = d ( n ) det M ( n ) + det ( [ M ( n ) 0 0 d ˜ ( n ) ] ( I n + 1 + [ M ( n ) 1 0 0 1 / d ˜ ( n ) ] [ 0 m ( n ) m ( n ) * 0 ] ) ) = d ( n ) det M ( n ) + d ˜ ( n ) det M ( n ) det ( I n + 1 + [ M ( n ) 1 0 0 1 / d ˜ ( n ) ] [ 0 m ( n ) m ( n ) * 0 ] ) ; n ( n 0 ) Z +
if { M ( n ) 1 } n = n 0 exists so that det M ( n ) 0 , that is, if 0 < k M ( n ) = | λ min ( M ( n ) ) | K M ( n ) = | λ max ( M ( n ) ) | < + ; n ( n 0 ) Z + . Thus,
| det M ( n + 1 ) | | d ( n ) det M ( n ) | + | d ˜ ( n ) det M ( n ) | | det ( I n + 1 + [ M ( n ) 1 0 0 1 / d ˜ ( n ) ] [ 0 m ( n ) m ( n ) * 0 ] ) | n K d ( n ) K M ( n ) + n K d ˜ ( n ) K M ( n ) ( 1 + 1 k M ( n ) k d ˜ ( n ) m ( n ) * m ( n ) ) n + 1 n K d ( n ) K M ( n ) + n K d ˜ ( n ) K M ( n ) ( 1 + ε ( n ) ) n + 1 n K d ( n ) K M ( n ) + n K d ˜ ( n ) K M ( n ) ( 1 + ( 2 n + 1 1 ) ε ( n ) ) ; n ( n 0 ) Z +
since λ max ( m ( n ) * m ( n ) ) = m ( n ) * m ( n ) ; n ( n 0 ) Z + . Since it is assumed the existence of a real sequence { ε ( n ) } n = n 0 [ 0 , 1 ) such that 1 k M ( n ) k d ˜ ( n ) m ( n ) * m ( n ) ε ( n ) ; n ( n 0 ) Z + and since
( 1 + ε ( n ) ) n + 1 = i = 0 n + 1 ( n + 1 i ) ε ( n ) i = 1 + i = 1 n + 1 ( n + 1 i ) ε ( n ) i 1 + ( 2 n + 1 1 ) ε ( n ) ; n ( n 0 ) Z +
thus, it follows that | det M ( n + 1 ) | ( n + 1 ) K M ( n + 1 ) n K M ( n ) holds for any given n ( n 0 ) Z + if
n K d ( n ) K M ( n ) + n K d ˜ ( n ) K M ( n ) ( 1 + ( 2 n + 1 1 ) ε ( n ) ) n K M ( n )
for any given n ( n 0 ) Z + , that is, if
n K M ( n ) ( K d ( n ) 1 ) + n K d ˜ ( n ) K M ( n ) ( 1 + ( 2 n + 1 1 ) ε ( n ) ) 0
for such n Z + , equivalently, if n K d ˜ ( n ) ( 1 + ( 2 n + 1 1 ) ε ( n ) ) n ( 1 K d ( n ) ) , equivalently if K d ( n ) 1 , and furthermore, and accordingly to the restriction on the sequence { ε ( n ) } n = n 0 , if
K d ˜ ( n ) ( 1 K d ( n ) ) 1 + ( 2 n + 1 1 ) ε ( n ) ( 1 K d ( n ) ) k M ( n ) k d ( n ) k M ( n ) k d ( n ) + ( 2 n + 1 1 ) λ max ( m ( n ) * m ( n ) )
for such an n Z 0 + . In particular, if | d ( n ) | = K d ( n ) = k d ( n ) 1 and | d ˜ ( n ) | = K d ˜ ( n ) = k d ˜ ( n ) for some given n ( n 0 ) Z 0 + , then (A4) holds if
| d ˜ ( n ) | ( 1 | d ( n ) | ) 1 + ( 2 n + 1 1 ) ε ( n ) ( 1 | d ( n ) | ) k M ( n ) | d ( n ) | k M ( n ) | d ( n ) | + ( 2 n + 1 1 ) λ max ( m ( n ) * m ( n ) )
and Property (i) has been proved. On the other hand, since Property (ii) assumes that the constraints of Property (i) hold the constraint (1) is a direct consequence of Property (i). Also, if M ( n ) _ 0 and M ( n + 1 ) _ 0 then the constraints (2) follow directly from (1) and Cauchy’s interlacing theorem of the eigenvalues which are real non-negative. Property (ii) has been proved. Property (iii) follows since lim sup n ( det | M ( n + 1 ) | det | M ( n ) | ) 0 directly from (A1) and the given assumptions. □
Proof of Lemma 3.
By the inversion of a block partitioned matrix, one gets:
M ( n + 1 ) 1 = [ M ( n ) m ( n ) m ( n ) d ( n ) + d ˜ ( n ) ] 1 = [ M ¯ ( n ) m ¯ ( n ) m ¯ ( n ) * d ¯ ( n ) + d ˜ ¯ ( n ) ]
where
M ¯ ( n ) = M ( n ) 1 ( I n + m ( n ) ( d ( n ) + d ˜ ( n ) m ( n ) M ( n ) 1 m ( n ) ) 1 m ( n ) M ( n ) 1 )
m ¯ ( n ) = M ( n ) 1 m ( n ) ( d ( n ) + d ˜ ( n ) m ( n ) M ( n ) 1 m ( n ) ) 1
d ¯ ( n ) + d ˜ ¯ ( n ) = ( d ( n ) + d ˜ ( n ) m ( n ) M ( n ) 1 m ( n ) ) 1
since M ( n ) is non-singular and d ˜ ( n ) m ( n ) M ( n ) 1 m ( n ) d ( n ) (i.e., ( d ( n ) + d ˜ ( n ) m ( n ) M ( n ) 1 m ( n ) ) is non-singular). The proof follows directly from Lemma 2 [(ii)-(iii)] by replacing M ( n + 1 ) M ¯ ( n + 1 ) = M ¯ ( n + 1 ) 1 ; n ( n 0 ) Z + . □
Proof of Lemma 4.
Since M ( n 0 ) 0 then its inverse is also positive definite and then both of them fulfil the positive semi-definiteness constraint of Lemma 2 [(ii),(iii)] and Lemma 3. The proof follows directly since sup n n 0 max ( | det M ( n ) | , | det M ( n ) 1 | ) < + and lim sup n ( sup max ( | det M ( n ) | , | det M ( n ) 1 | ) ) < + then M ( n ) 0 for n n 0 if M ( n 0 ) 0 , lim sup n M ( n ) 0 and lim inf n M ( n ) 0 . □
Proof of Lemma 5.
Note that M ( n 0 ) = A ( n 0 ) * A ( n 0 ) 0 since A ( n 0 ) is a stability matrix. From Lemma 2 and Lemma 3, the sequence { M ( n ) } n = n 0 consists of positive definite members. Thus, the singular values of the elements of the sequence { A ( n ) } n = n 0 are positive and bounded. Since A ( n 0 ) is stable and since the eigenvalues of any square matrix are continuous functions of its entries there is no zero eigenvalue in any member of the sequence { A ( n ) } n = n 0 . Any member of this sequence has no eigenvalues at the imaginary complex axis other than zero (i.e., any nonzero critically stable eigenvalues) since then the corresponding M ( n ) = A ( n ) * A ( n ) is not positive definite contradicting the given assumption. As a result, no member of the sequence { A ( n ) } n = n 0 has a critically stable eigenvalue (i.e., located on the imaginary complex axis) or unstable eigenvalue (i.e., located on the complex open right half plane).
M ( n 0 ) 0 for some given arbitrary n 0 Z + and assume also that the conditions of Lemma 2 (iii) and Lemma 3 hold. Then, M ( n + 1 ) 0 ; n ( n 0 ) Z + . □
Proof of Lemma 7.
It is direct since lim m { A ( n ) m } n = 0 = 0 implies that lim m { | A ( n + ξ ) A ( n ) | m } n = 0 = 0 for any given ξ Z + . □
Proof of Lemma 8.
If for any given n Z 0 + , A ( n ) 2 < 1 then lim m A ( n ) m = 0 for any n Z 0 + , then lim m { A ( n ) m } n = 0 = 0 , M ( n ) 2 2 = A ( n ) * A ( n ) 2 2 A ( n ) 2 4 < 1 and then M ( n ) 2 < 1 , lim m M ( n ) m = 0 and lim m { M ( n ) m } n = 0 = 0 for any given n Z 0 + . Thus, if { A ( n ) } n = 0 is convergent then { M ( n ) } n = 0 is convergent, hence Property (i) holds. On the other hand, since M ( n ) is semidefinite positive Hermitian by construction; n Z 0 + then the condition M ( n ) 2 < 1 ; n Z 0 + leads to the convergence of { M ( n ) } n = 1 , and if M ( n ) 2 < 1 for some n Z 0 + then
r 2 ( M ( n ) ) = λ max ( M ( n ) * M ( n ) ) = λ max 2 ( A ( n ) * A ( n ) ) = λ max 2 ( M ( n ) ) = M ( n ) 2 2 < 1
Then, A ( n ) 2 = λ max 1 / 2 ( A ( n ) * A ( n ) ) = M ( n ) 2 1 / 2 < 1 . Thus, if the above holds for any n Z 0 + , one concludes that { A ( n ) } n = 0 is convergent if { M ( n ) } n = 0 is convergent. Hence, Property (ii) follows. Property (iii) is a combination of the other two properties since for any n Z 0 + , M ( n ) 2 < 1 , A ( n ) 2 < 1 implies that M ( n ) 2 < 1 and M ( n ) 2 < 1 implies that A ( n ) 2 < 1 . □
Proof of Theorem 1.
Properties (i),(ii) follows from Lemma 7and Lemma 8 by taking into account the factorization (28). Property (iii) follows from Property (i) and (20) since
x ¯ ( n + m ) ( i = 0 m 1 A ¯ ( n i ) ) x ¯ ( n ) max 0 i m 1 A ¯ ( n i ) m x ¯ ( n ) ;   n Z 0 + , m Z +
and then x ¯ ( n + m ) 0 as m for any n Z 0 + . Property (iv) is proved in the same way as Property (iii) via Property (ii). The proof of Property (v) is made by comparing (28) with (10)–(12) by replacing λ ( n ) M ¯ ( n ) A ( n ) T B ( n ) and Δ ( n ) B ( n ) T B ( n ) . One gets, via complete induction, that if M ( 0 ) is convergent and, furthermore,
λ max ( B ( n ) T A ( n ) A ( n ) T B ( n ) ) < ( 1 ε ( n + 1 ) max ( 1 ε ( n ) , | λ max ( B ( n ) T B ( n ) ) | ) ) 2 ; n Z 0 +
Then, { M ( n ) } n = 0 is convergent, since M ( 0 ) is convergent, provided that ε ( n + 1 ) < ε ( n ) , that is { ε ( n ) } n = 0 [ 0 , 1 ) is strictly decreasing, and
| λ max ( B ( n ) T B ( n ) ) | < 1 ε ( n + 1 ) ;   n Z 0 +
The constrains (A10), (A11) are jointly fulfilled if | λ max ( B ( n ) T B ( n ) ) | < min ( 1 ε ( n + 1 ) , ε ( n ) ε ( n + 1 ) ) ; n Z 0 + .□
Proof of Lemma 9.
Since A m = E A = I E A then v e c ( E A ) = v e c ( I E A ) = ( I A T ) v e c ( E ) , and (32) becomes for I being the identity matrix of the same order as E and A :
( B C T ) v e c ( K ) = v e c ( ( E I ) A ) = ( I A T ) v e c ( E ) v e c ( A )
whose solutions are given by (34) if (33) holds and the whole set of solutions reduces to (35) if the solution is unique. The corresponding Equations (36)–(38) are got for the particular case when A m = ρ A with | ρ | < 1 / A 2 . Property (i) has been proved. Property (ii) follows since the least-squares best approximation to the corresponding incompatible algebraic system (31), or (32), is (40), that is, (34) for k w = 0 , [13,14]. □
Proof of Theorem 2.
The solvability of (41) in the form (42) follows from the Rouché-Froebenius rank conditions from the algebraic compatibility under Assumptions 1,2. By defining
A 0 ( n + 1 ) = [ A ( n ) 0 ( i = 0 n + 1 n i ) × n n + 1 ] ;   B 0 ( n + 1 ) = B ( n )
in order to complete a square “a priori” matrix of dynamics of the ( n + 1 ) -the aggregated system obtained after the aggregation of the ( n + 1 ) -th subsystem, note that
A 0 ( n + 1 ) = [ A 0 ( n ) + B 0 ( n ) K 0 ( n ) C ( n ) A ^ 0 ( n ) A ^ ( n + 1 ) + B ^ 0 ( n ) K ^ 0 a ( n ) C ( n ) D ^ ( n ) + D ^ ˜ ( n ) + B ^ 0 ( n ) K ^ 0 ( n ) C ^ ( n )    0 0 ] R ( i = 0 n + 1 n i ) × ( i = 0 n + 1 n i )
B 0 ( n + 1 ) = [ ( B 1 ( n ) + B 0 ( n ) K 1 ( n ) ) C ( n 1 ) 0 ( B ^ 1 a ( n ) + B ^ 0 ( n ) K ^ 1 a ( n ) ) C ( n 1 ) ( B ^ 1 ( n ) + B ^ 0 ( n ) K ^ 1 ( n ) ) C ^ ( n 1 ) ( B j ( n ) + B 0 ( n ) K r ( n ) ) C ( n r ) 0 ( B ^ r a ( n ) + B ^ 0 ( n ) K ^ r a ( n ) ) C ( n r ) ( B ^ j ( n ) + B ^ 0 ( n ) K ^ r ( n ) ) C ^ ( n r ) ] R ( i = 0 n + 1 n i ) × ( j = n r n 1 ( i = 0 j + 1 n i ) )
n Z 0 + . From (41), with corresponding associated controller explicit solutions (42), one gets that the ( n + 1 ) -th aggregated delay-free dynamics is described by the matrix:
A 0 ( n + 1 ) = [ A ( n ) 0 ( i = 0 n + 1 n i ) × n n + 1 ] = [ A f ( n ) A ^ 0 ( n ) A ^ f a ( n + 1 ) A ^ f ( n ) 0 0 ] R ( i = 0 n + 1 n i ) × ( i = 0 n + 1 n i ) ;   n Z 0 +
Having in mind (27), construct
M ( n + 1 ) = A ¯ ( n ) A ¯ ( n ) T = [ A f ( n ) A ^ 0 ( n ) A ^ f a ( n + 1 ) A ^ f ( n )    0 0    B ¯ 1 0 ( n ) B ¯ 2 0 ( n ) i ( B ¯ 0 ( n ) ) I       0 ] [ A f ( n ) A ^ 0 ( n ) A ^ f a ( n + 1 ) A ^ f ( n ) 0 0 B ¯ 1 0 ( n ) B ¯ 2 0 ( n ) i ( B ¯ 0 ( n ) ) I 0 ] T = [ A f ( n ) A f ( n ) T + A ^ 0 ( n ) A ^ 0 ( n ) T + B ¯ 1 0 ( n ) B ¯ 1 0 ( n ) T A f ( n ) A ^ f a ( n + 1 ) T + A ^ 0 ( n ) A ^ f ( n ) T + B ¯ 1 0 ( n ) B ¯ 2 0 ( n ) T [ A f ( n ) A ^ 0 ( n ) A ^ f a ( n + 1 ) A ^ f ( n ) ] [ A f ( n ) A ^ 0 ( n ) A ^ f a ( n + 1 ) A ^ f ( n ) ] T i ( B ¯ 0 ( n ) ) I ] = [ A f ( n ) A f ( n ) T + A ^ 0 ( n ) A ^ 0 ( n ) T + B ¯ 1 0 ( n ) B ¯ 1 0 ( n ) T A f ( n ) A ^ f a ( n + 1 ) T + A ^ 0 ( n ) A ^ f ( n ) T + B ¯ 1 0 ( n ) B ¯ 2 0 ( n ) T [ A f ( n ) A ^ f a ( n + 1 ) ] [ A ^ 0 ( n ) A ^ f ( n ) ] [ A f ( n ) T A ^ f a ( n + 1 ) T ] [ A ^ f 0 ( n ) T A ^ f ( n ) T ]      i ( B ¯ 0 ( n ) ) δ ( n ) I ] + i ( B ¯ 0 ( n ) ) [ 0 0 0 ( 1 δ ( n ) ) I ]
where
B ¯ 0 ( n ) = [ B ¯ 1 0 ( n ) B ¯ 2 0 ( n ) ] = [ ( B 1 ( n ) + B 0 ( n ) K 1 ( n ) ) C ( n 1 ) 0 ( B ^ 1 a ( n ) + B ^ 0 ( n ) K ^ 1 a ( n ) ) C ( n 1 ) ( B ^ 1 ( n ) + B ^ 0 ( n ) K ^ 1 ( n ) ) C ^ ( n 1 ) ( B r ( n ) + B 0 ( n ) K r ( n ) ) C ( n r ) 0 ( B ^ r a ( n ) + B ^ 0 ( n ) K ^ r a ( n ) ) C ( n r ) ( B ^ r ( n ) + B ^ 0 ( n ) K ^ r ( n ) ) C ^ ( n r ) ]
and
(a)
i ( B ¯ 0 ( n ) ) is a binary indicator function defined as i ( B ¯ 0 ( n ) ) = 1 if ( B ¯ 0 ( n ) ) 0 and i ( B ¯ 0 ( n ) ) = 0 if ( B ¯ 0 ( n ) ) = 0 . The reason of the use of this indicator is that, in fact, if the delayed dynamics is zero then the dimension of the extended state, so that of A ¯ ( n ) , decreases since the resulting block identity matrices are removed,
(b)
{ δ ( n ) } n = 0 [ 0 , 1 ] is a design sequence which satisfies { δ ( n ) } n = 0 0 .
Note that the fact that I = i ( B ¯ 0 ( n ) ) I = i ( B ¯ 0 ( n ) ) δ ( n ) I + i ( B ¯ 0 ( n ) ) ( 1 δ ( n ) ) I justifies (A16). □

References

  1. Ibeas, A.; De La Sen, M. Robustly stable adaptive control of a tandem of master–slave robotic manipulators with force reflection by using a multiestimation scheme. IEEE Trans. Syst. Man Cybern. Part B (Cybernetics) 2006, 36, 1162–1179. [Google Scholar] [CrossRef]
  2. Barambones, O.; Garrido, A.J.; Garrido, I. Robust speed estimation and control of an induction motor drive based on artificial neural networks. Int. J. Adapt. Control Signal Process. 2008, 22, 440–464. [Google Scholar] [CrossRef]
  3. Bakule, L.; de la Sen, M. Decentralized stabilization of networked complex composite systems with nonlinear perturbations. In Proceedings of the 2009 IEEE International Conference on Control and Automation, Christchurch, New Zealand, 9–11 December 2009; pp. 2272–2277. [Google Scholar]
  4. Singh, M.G. Decentralised Control; North-Holland: Systems and Control Series; North Holland Publishing Company: New York, NY, USA, 1981; Volume 1. [Google Scholar]
  5. Berman, A.; Plemmons, R.J. Nonnegative Matrices in the Mathematical Sciences; Academic Press: New York, NY, USA, 1979. [Google Scholar]
  6. Kailath, T. Linear Systems; Prentice-Hall Inc.: Englewood Cliffs, NJ, USA, 1980. [Google Scholar]
  7. Ortega, J.M. Numerical Analysis; Academic Press: New York, NY, USA, 1972. [Google Scholar]
  8. Kouachi, S. The Cauchy interlace theorem for symmetrizable matrices. arXiv Preprint arXiv:1603.04151 2016, 2016.
  9. Almutairi, S.; Thenmozhi, V.; Watkins, J.; Sawan, M.E. Optimal Control Design of Large Scale Systems with Uncertainty. In Proceedings of the 2018 IEEE 61st International Midwest Symposium on Circuits and Systems (MWSCAS), Windsor, ON, Canada, 5–8 August 2018; pp. 372–375. [Google Scholar]
  10. Auger, P.; de La Parra, R.B.; Poggiale, J.C.; Sánchez, E.; Sanz, L. Aggregation methods in dynamical systems and applications in population and community dynamics. Phys. Life Rev. 2008, 5, 79–105. [Google Scholar] [CrossRef] [Green Version]
  11. Aoki, M. Control of large-scale dynamic systems by aggregation. IEEE Trans. Autom. Control 1968, 13, 246–253. [Google Scholar] [CrossRef]
  12. Darbha, S.; Rajagopal, K.R. Aggregation of a class of linear interconnected systems. In Proceedings of the American Control Conference, San Diego, CA, USA, 2–4 June 1999; pp. 1496–1501. [Google Scholar]
  13. Mackenroth, U. Robust Control Systems. Theory and Case Studies; Springer: Berling/Heidelberg, Germany, 2003. [Google Scholar]
  14. Barnett, S. Matrices in Control Theory with Applications to Linear Programming; Van Nostrand Reinhold Company: London, UK, 1971. [Google Scholar]
  15. Verma, R.; Sehgal, V.K. Computational stochastic modelling to handle the crisis occurred during community epidemic. Ann. Data Sci. 2016, 3, 119–133. [Google Scholar] [CrossRef]
  16. Iggidr, A.; Souza, M.O. State estimators for some epidemiological systems. J. Math. Biol. 2019, 78, 225–256. [Google Scholar] [CrossRef] [PubMed]
  17. Yang, H.M.; Freitas, A.R. Biological view of vaccination described by mathematical modellings: from rubella to dengue vaccines. Math. Biosci. Eng. 2019, 16, 3195–3214. [Google Scholar] [CrossRef]
  18. De la Sen, M.; Agarwal, R.P.; Ibeas, A.; Alonso-Quesada, S. On a generalized time-varying SEIR epidemic model with mixed point and distributed time-varying delays and combined regular and impulsive vaccination controls. Adv. Differ. Equ. 2010, 1, 1–42. [Google Scholar] [CrossRef]
  19. De la Sen, M.; Ibeas, A.; Alonso-Quesada, S.; Nistal, R. On a SIR model in a patchy environment under constant and feedback decentralized controls with asymmetric parameterizations. Symmetry 2019, 11, 430. [Google Scholar] [CrossRef]
  20. Hu, S.A.; Smith, R.L. The Schur complement interlacing theorem. SIAM J. Matrix Anal. Appl. 1995, 16, 1013–1023. [Google Scholar] [CrossRef]
  21. Mercer, A.M.; Mercer, P.R. Cauchy’s interlace theorem and lower bounds for the spectral radius. Int. J. Math. Math. Sci. 2000, 23, 563–566. [Google Scholar] [CrossRef]
  22. Hwang, S.G. Cauchy’s interlace theorem for eigenvalues of Hermitian matrices. Am. Math. Mon. 2004, 111, 157–159. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

De la Sen, M. On Cauchy’s Interlacing Theorem and the Stability of a Class of Linear Discrete Aggregation Models Under Eventual Linear Output Feedback Controls. Symmetry 2019, 11, 712. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11050712

AMA Style

De la Sen M. On Cauchy’s Interlacing Theorem and the Stability of a Class of Linear Discrete Aggregation Models Under Eventual Linear Output Feedback Controls. Symmetry. 2019; 11(5):712. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11050712

Chicago/Turabian Style

De la Sen, Manuel. 2019. "On Cauchy’s Interlacing Theorem and the Stability of a Class of Linear Discrete Aggregation Models Under Eventual Linear Output Feedback Controls" Symmetry 11, no. 5: 712. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11050712

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop