A discussion of stability analysis for systems of differential equations with multiple and distributed delays

In this article, we consider a class of systems of multiple delay differential equations (MDDEs). We ﬁrst deﬁne a characteristic matrix equation that can be used to analyze the stability of the equilibrium of a system of MDDEs. Then we construct a matrix based on the coefﬁcients of the characteristic matrix equation and use the spectrum of this matrix to derive necessary and sufﬁcient conditions for the system to be stable. Next we discuss a comparison of the stability equivalency between a system of delay differential equations (DDEs) to the system of MDDEs and relate our results to distributed delay systems (DDSs). Numerical examples are given to justify our theory.


Introduction
Studying the stability of delay differential equations (DDEs) has become increasingly important in recent times; see [1][2][3][4][5][6] and the growing body of literature in the field.In the past few years, there has been a significant amount of research focused on MDDEs; see [3,[7][8][9][10][11][12], as well as their applications in various fields such as the dynamics of electrical power systems, macroeconomic models, electricity market models, and more, as outlined in [8,[13][14][15][16][17][18][19].In addition to these areas, there has been an increasing interest in studying other types of systems that exhibit memory effects, including systems of fractional differential and difference equations, as well as systems of fractional nabla difference equations, as highlighted in [20][21][22].Conventional approaches for the stability analysis of DDEs are based on Lyapunov functional method (LFM) and techniques that require the solution of a linear matrix inequality (LMI) problem; see [23][24][25].The complexity to construct the Lyapunov function and the heavy computational burden to solve the LMI problem limit the B Ioannis Dassios ioannis.dassios@ucd.ieMuyang Liu Muyang.liu@xju.edu.cnFederico Milano federico.milano@ucd.ie 1 University College Dublin, Dublin, Ireland 2 Xinjiang University, Ürümqi, China application of LFMs on engineering fields.Moreover, as LFMs provide only sufficient but not necessary conditions for system stability, they tend to be conservative.There also exist a variety of frequency-domain approaches to solve the stability of DDEs, [2,4,[26][27][28][29][30]. Most of these techniques are based on the solution of an eigenvalue problem.This consists in estimating the dominant modes of the DDEs through the solution of the characteristic equation of the system.In [6,10,16,17], a general eigenvalue analysis approach is developed to solve the stability of large system described by a set of delay differential algebraic equations (DDAEs).Compared to LFMs, eigenvalue-based approaches are less computationally intensive and provide a more accurate stability analysis.For this reason, we consider an eigenvalue-based approach also in this paper.
In this article, we will firstly provide stability criteria for a class of systems of DDEs with multiple delays based on eigenvalue analysis.Additionally, we explore the relationship between the stability of a DDE system with one delay and a system of MDDEs.This discussion offers a fresh perspective and novel insights that may inspire future research in this field.We are interested in the evaluation of the small signal stability of a nonlinear system of MDDEs in the following form: where τ i > 0 is constant time delay, Y ∈ R m are the state variables, f : R (n+1)×m → R m are the differential equa-tions and can be implicit with its partial derivatives at an equilibrium point to be singular matrices, i.e., det We consider only small disturbances, e.g., disturbances whose effects on the stability of a given equilibrium point can be studied through the linearized set of the equations that model the system.If we consider small disturbances, e.g., disturbances whose effects on the stability of a given equilibrium can be studied through the linearized set of the equations that model the system.The linearized systems of MDDEs have the form: where δY = Y − Y eq , Y eq is equilibrium and f Ẏ full rank at an equilibrium.The characteristic equation is then given by and its characteristic roots will provide the necessary information for small signal stability of the system of MDDE.To sum up the small signal stability of the nonlinear MMDE at a given equilibrium can be studied from the following linear system MDDE: ( with characteristic equation where A necessary and sufficient condition for the equilibrium solution to be asymptotically stable is that the roots of the characteristic equation all have negative real parts; see [3]. In the remainder of the paper I m denotes the identity matrix m × m, 0 i, j the zero matrix of i rows, j columns, and • a natural norm.The remainder of this paper is organized as follows.In Sect.2, we present a theorem that establishes a stability criterion for MDDE systems in the form of (1).In Sect.3, we explore the stability equivalence between systems of DDEs and systems of MDDEs and discuss the relevance of our findings for DDSs.Section 4 is devoted to numerical examples, which demonstrate the effectiveness of our approach.Finally, we summarize our key findings in the Conclusions section.

Stability analysis of multiple delay systems
In this section, we present our main results.We will use the following definition: Initially we consider τ k = kτ in (1) and provide the following Theorem: Theorem 2.1 We consider system (1) with τ k = kτ .Then the following conditions must hold for delay-independent stability of (1): 1.The matrix A 0 is stable; 2. The matrix A 0 + n k=1 A k is stable; 3. The spectral radius of A ω is less than 1, ∀ω > 0, where Hence, the matrix A 0 has to be stable in order to have stability for (1) at the equilibrium state.If τ = 0 then (2) takes the form which means that the matrix A 0 + n k=1 A k has to be stable in order to have stability for (1) at the equilibrium state.By applying the Fourier transform F(Y ) = (ω) into (1), we get: or, equivalently, Then det( jωI m − A 0 − n k=1 A k e − jωkτ ) = 0 is the characteristic equation of (1).We adopt the following notation: or, equivalently, or, equivalently, in matrix form where . . .
and F 1 , F 2 are given by It is worth noting that the system under consideration is a set of DDAEs.Therefore, since we have established the equivalence of the two systems, we can conclude that the characteristic equation of (1) is We have that: Then the characteristic equation of the DDAEs can be written as: Let ρ(•) be spectral radius of a matrix.Using Theorem 2.1 in [3] we get that ∀ω > 0, if ρ(F −1  ω [e − jωτ G]) < 1, then the set of DDAEs is stable independent of delay, where , and Note that ∀ω > 0: Since the set of DDAEs is equivalent to (1), we have that ∀ω > 0 if ρ(A ω ) < 1 holds then (1) is stable independent of delay.

Discussion on the equivalency of MDDEs to a system of DDEs and DDSs
We consider now the system of DDEs: where A, A i ∈ R m×m , i =, 1, . . ., n, Y : [0, +∞] → R m×1 , and A 0 ∈ R m×m is stable matrix.In the following discussion we will attempt to relate systems (1), (3) for small disturbances, i.e., , and investigate their stability such that if (3) is asymptotically stable, then (1) is also asymptotically stable.An implicit solution of system (1) is given by: or, equivalently, if we apply the Weierstrass canonical form; see [22], Similarly, an implicit solution of ( 3) is given by: whereby applying the Weierstrass canonical form we get: The matrices P, Q have as columns the left, and right respectively linear independent eigenvectors of A 0 , while J is the Jordan matrix of the eigenvalues of A 0 .By subtracting these two solutions we get: By applying a natural norm we get: whereby using the property of the norm we have that: If λ i is an eigenvalue of A 0 with algebraic multiplicity p i , the Jordan matrix has the form: where In addition: where By taking the norm • 1 of e J t we have Let By using ( 5), ( 6) into ( 4) we get: and hence for λ i < 0 we get: In the following remark we discuss the possibility of A = g(A i ) for certain cases.

Remark 3.1
We discussed the stability equivalency between a system of DDEs with one delay and a system of MDDEs.The significance of this aspect was to devise a new concept that can provide novel perspectives for researchers.By exploring this idea in our discussion, we aim to pave the way for future research to advance the concept of stability equivalence.We considered only small disturbances.Let in ( 1), (3), Hence in this special case by applying the Laplace transform into (1), ( 3) we can observe that the real parts of the rightmost eigenvalues of (3) should also converge to that of (1) if

Consequently, one of the practical options for
The proof of Theorem 2.5 in [3] provides the necessary and sufficient condition to ensure that for the special case that A = g(A k ) = n k=0 A k , there exists a τ ≤ τ such that the systems (1), (3) have the same stability assertion.Hence, for the case that A = g(A k ) = c n k=0 A k , particularly c = 1 is an appropriate selection to apply the idea described above, and one can obtain that the MDDE system (1) is asymptotically stable if the single delay system (3) with A = n k=0 A k and τ ≤ τ = n k=1 τ k n is asymptotically stable.The results discussed and obtained in Sect. 3 can also be used for the stability analysis of DDSs in the following form: where π(•) is the probability distribution of ξ that satisfies the following property: and π(ξ ) is a non-negative function.
From [28,31,32], we know that the DDS in the form of (7) has the same spectrum as the comparison system: and κ is a weighting parameter decided by the interpolation method used.With z m truncated at a finite value, one obtains the following system: which is a system of MDDEs in the form of (1), and its spectrum can be studied according to Theorem 2.1, and the discussions in this section.Hence, with a proper interpolation method and a fixed z m , the critical eigenvalues of the DDS (7) converge to those of system (8).

Numerical examples
In this section, we provide numerical examples to illustrate the theory presented in Sects. 2 and 3.
To utilize Theorem 2.1, the following steps should be followed.
Step 1 we compute the matrix jω − A 0 : Step 2 we compute the matrices ( jω  Step 3 we form the matrix A ω : Obviously the matrices A 0 , A 0 + 2 k=1 A k are both stable since all their eigenvalues real and negative; in addition, ρ(A ω ) = 0.0014 < 1 and hence from Theorem 2.1 the system of MDDEs is delay independent stable.

Example 4.2
We consider now the DDS (7) with: Then if z m = 2" for the comparison system (8) we have: Hence by setting κ = 2 τ the DDS (8) takes the form: or, equivalently, by setting τ = τ 2 : which is the system of MDDEs (1) for n = 2 and Then the matrices A 0 , A 0 + 2 k=1 A k are both stable since all their eigenvalues real and negative, and in addition ρ(A ω ) = 0.0014 < 1.Thus, by applying Theorem 2.1, it can be concluded that the system of MDDEs is delay independent stable.Since the DSS is interconnected with the system of MDDEs, it also inherits the same property of delay independence stability.

Conclusions
In this article, we derived simple and practical conditions for the small signal stability analysis of the MDDE system (1).Furthermore, we explored the connection between the stability analysis of MDDEs and that of DDEs.In addition to MDDEs, we discussed how the stability criterion can be extended to DDSs.We demonstrated the applicability of our approach with numerical examples, which showed that the proposed method is effective.In summary, this article provided a comprehensive approach for small signal stability analysis of MDDEs and DDSs.Our findings have significant practical implications for the design and analysis of these types of systems.
As a future direction, we plan to expand the scope of our research to encompass other types of systems that exhibit memory effects, such as systems of fractional differential and difference equations, see [21,22].Additionally, we intend to explore promising applications where delays are significant, such as in the dynamics of electrical power systems, macroeconomic models, and electricity market models, among others; see [8,13,15,18].Furthermore, we aim to extend our theoretical findings to systems that employ forward operators.Specifically, we plan to investigate systems where instead of a vector with delay in the form Y (t − τ ), a vector with delay in the form Y (t + τ ) is used, as discussed in [20,33].There is already some research in progress in these areas, and we hope to contribute to these efforts with our work.

Example 4 . 1
We consider the system of MDDEs (1) for n =

Table 1
The stability assertions of (9) according to Theorem 2.1

Table 2
The rightmost eigenvalues of the different scenarios, where λ k are the rightmost eigenvalues of (9); λ