Journal of Applied Mathematics

Journal of Applied Mathematics / 2012 / Article
Special Issue

Modeling and Control of Complex Dynamic Systems: Applied Mathematical Aspects

View this Special Issue

Research Article | Open Access

Volume 2012 |Article ID 749847 |

Jian Ma, Baodong Zheng, Chunrui Zhang, "A Matrix Method for Determining Eigenvalues and Stability of Singular Neutral Delay-Differential Systems", Journal of Applied Mathematics, vol. 2012, Article ID 749847, 11 pages, 2012.

A Matrix Method for Determining Eigenvalues and Stability of Singular Neutral Delay-Differential Systems

Academic Editor: Zhiwei Gao
Received29 Nov 2011
Revised06 Apr 2012
Accepted08 Apr 2012
Published04 Jun 2012


The eigenvalues and the stability of a singular neutral differential system with single delay are considered. Firstly, by applying the matrix pencil and the linear operator methods, new algebraic criteria for the imaginary axis eigenvalue are derived. Second, practical checkable criteria for the asymptotic stability are introduced.

1. Introduction

Nowadays, the time-delay systems have became an important natural models in physics, engineering, multibody mechanics, computer-aided design, and economic systems. The theory on ordinary differential equations with delays have been discussed for decades in a wide range, so there are very many results for them. Especially, the eigenvalues and the stability analysis of time-delay systems have received much attention of researchers and many excellent results have been obtained, see [1ā€“6]. Certainly most of them had been focused on the analytical methods or numerical methods, such as V-functional methods, Laplace transformation, Runge-Kutta methods, and linear multistep methods. In [7ā€“9], the numerical techniques for the computation of the eigenvalues were discussed. In [10], Zhu and Petzold researched the asymptotic stability of delay-differential-algebraic equations by applying the šœƒ-methods, Runge-Kutta methods, and linear multistep methods. These methods play the key roles at last. But in recent years, algebraic methods are developing fast, especially for the research on the more complex systems, such as the n-dimensional systems. Though the algebraic methods as a new and effective tool is also applied to analyze the time-delay systems [2, 3], the results are very few.

In this paper, we will discuss the differential-algebraic equations by the algebraic methods. Their dynamics have not been well understood yet.

Example 1.1. Consider the simple differential-algebraic system: Ģ‡š‘„1(š‘”)=š‘“1š‘„(š‘”),1(š‘”)āˆ’š‘„2(š‘”āˆ’šœ)=š‘“2(š‘”),(1.1) where š‘”ā‰„0, šœā‰„0, and š‘„1 and š‘„2 are given by continuous functions on the initial interval (āˆ’šœ,0]. So we have the solution: š‘„1(ī€œš‘”)=š‘”0š‘“1(š‘„š‘ )š‘‘š‘ +š‘,2(š‘”)=āˆ’š‘“2ī€œ(š‘”+šœ)+0š‘”+šœš‘“1(š‘ )š‘‘š‘ +š‘,(1.2) where š‘ is a constant. From the solution, we find that the solution depends on future integrals of the input š‘“(š‘”). This interesting phenomenon arrested many scholars to research. For a general š‘›-dimensional differential equations with delay, we can note by šµĢ‡š‘‹(š‘”)=š“0š‘‹(š‘”)+š“1š‘‹(š‘”āˆ’šœ),(1.3) where šµ,š“0,š“1āˆˆš‘…š‘›Ć—š‘›, Rankšµā‰¤š‘›, š‘”ā‰„0, šœā‰„0, and š‘‹(š‘”)āˆˆš‘…š‘› is given by continuous functions on the initial interval (āˆ’šœ,0]. When Rankšµ=š‘›, we called it the retarded differential equations. It can be improved as Ģ‡š‘‹(š‘”)=š“0š‘‹(š‘”)+š“1š‘‹(š‘”āˆ’šœ).(1.4) Many scholars have widely researched the delay-independent or delay-dependent stability and asymptotic stability by analytic methods or numerical methods. When Rankšµ<š‘›, it is called a singular (or degenerated) delay-differential equations. The imaginary axis eigenvalues are discussed by using matrix pencil, see [2]. But because of the complex nature of the singular differential systems with delay, the research is very difficult by using the analytical treatment. So few studies on the stability and the bifurcations have been conducted so far. Particulary, for the singular neutral differential systems with delays, there are hardly flexible and efficient verdicts.
In this paper, we will apply the algebraic methods to discuss the stability of a singular neutral differential system with a single delay, as follows: šµ0Ģ‡š‘„(š‘”)+šµ1Ģ‡š‘„(š‘”āˆ’šœ)=š“0š‘„(š‘”)+š“1š‘„(š‘”āˆ’šœ),(1.5) where šµ0,šµ1,š“0,š“1āˆˆā„š‘›Ć—š‘›, š‘„āˆˆā„š‘›, šœā‰„0. For the system (1.5), if detšµ0ā‰ 0, we can improve it as the form Ģ‡š‘„(š‘”)+šµ1Ģ‡š‘„(š‘”āˆ’šœ)=š“0š‘„(š‘”)+š“1š‘„(š‘”āˆ’šœ).(1.6) It is a neutral differential equation with a single delay. The problem of computing imaginary axis eigenvalues on the system (1.6) has been previously studied in [11]. Here, we consider the state rankšµ0ā‰¤š‘›. The solvability of the system (1.5), which is essentially the existence and the uniqueness of the solution, is determined by the regularity. The matrix pencil (šµ0,š“0) is said to be regular if š‘ šµ0+š“0 is not identically singular for any complex š‘ . If (šµ0,š“0) is regular, the zero š‘  of det(š‘ šµ0+š“0) is called the eigenvalue of the matrix pencil (šµ0,š“0). From [12], we know that the system (1.5) is solvable if and only if (šµ0,A0) is regular. So, in this paper, we suppose that (šµ0,š“0) is regular. In the following, we will analyze the eigenvalues and the stability of the system (1.5).

2. The Algebraic Criteria for Determining Imaginary Axis Eigenvalues

Firstly, we research an ordinary differential equation, which will motivate our analysis. Consideršµ0Ģ‡š‘‹(š‘”)+šµ1Ģ‡š‘Œ(š‘”)=š“0š‘‹(š‘”)+š“1Ģ‡š‘Œ(š‘”),š‘‹(š‘”)šµš‘‡1+Ģ‡š‘Œ(š‘”)šµš‘‡0=āˆ’š‘‹(š‘”)š“š‘‡1āˆ’š‘Œ(š‘”)š“š‘‡0,(2.1) where šµ0,šµ1,š“0,š“1āˆˆā„š‘›Ć—š‘›, š‘‹,š‘Œāˆˆā„‚š‘›Ć—š‘›. Let š‘‰ denote the vector space š‘‰=ā„‚š‘›Ć—š‘›Ć—ā„‚š‘›Ć—š‘› and š„,š… denote the operators on š‘‰, given byš„āŽ›āŽœāŽœāŽš‘‹š‘ŒāŽžāŽŸāŽŸāŽ =āŽ›āŽœāŽœāŽšµ0š‘‹+šµ1š‘Œš‘‹šµš‘‡1+š‘Œšµš‘‡0āŽžāŽŸāŽŸāŽ āŽ›āŽœāŽœāŽš‘‹š‘ŒāŽžāŽŸāŽŸāŽ =āŽ›āŽœāŽœāŽš“,š…0š‘‹+š“1š‘Œāˆ’š‘‹š“š‘‡1āˆ’š‘Œš“š‘‡0āŽžāŽŸāŽŸāŽ ,āˆ€š‘‹,š‘Œāˆˆā„‚š‘›Ć—š‘›.(2.2) With ī€·š‘(t)=š‘‹(š‘”)š‘Œ(š‘”)ī€ø, the system (2.1) can be written asš„Ģ‡š‘(š‘”)=š…š‘(š‘”).(2.3) Supposing ī‚ī‚€ī‚š‘‹ī‚š‘Œī‚=ī‚€š‘(š‘”)=š‘‹0š‘’š‘ š‘”š‘Œ0š‘’š‘ š‘”ī‚ is a matrix solution of the system (2.1), we haveī€·š‘ šµ0āˆ’š“0ī€øī‚ī€·š‘‹+š‘ šµ1āˆ’š“1ī€øī‚ī‚š‘‹ī€·š‘Œ=0,š‘ šµš‘‡1+š“š‘‡1ī€ø+ī‚š‘Œī€·š‘ šµš‘‡0+š“š‘‡0ī€ø=0.(2.4) For any complex š‘ , let š“=š“(š‘ ) be the operator š“=š‘ š„āˆ’š…; thenš“āŽ›āŽœāŽœāŽš‘‹š‘ŒāŽžāŽŸāŽŸāŽ =āŽ›āŽœāŽœāŽī€·š‘ šµ0āˆ’š“0ī€øī€·š‘‹+š‘ šµ1āˆ’š“1ī€øš‘Œš‘‹ī€·š‘ šµš‘‡1+š“š‘‡1ī€øī€·+š‘Œš‘ šµš‘‡0+š“š‘‡0ī€øāŽžāŽŸāŽŸāŽ ,āˆ€š‘‹,š‘Œāˆˆā„‚š‘›Ć—š‘›.(2.5) For any complex š‘ , let š“+=š“+(š‘ )āˆ¶š‘‰ā†’š‘‰, satisfyingš“+āŽ›āŽœāŽœāŽš‘‹š‘ŒāŽžāŽŸāŽŸāŽ =āŽ›āŽœāŽœāŽš‘‹ī€·š‘ šµš‘‡0+š“š‘‡0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øš‘Œī€·āˆ’š‘‹š‘ šµš‘‡1+š“š‘‡1ī€ø+ī€·š‘ šµ0āˆ’š“0ī€øš‘ŒāŽžāŽŸāŽŸāŽ ,āˆ€š‘‹,š‘Œāˆˆā„‚š‘›Ć—š‘›.(2.6) For any complex š‘ , let Ī›=Ī›(š‘ ), satisfyingī€·šš²š‘‹=š‘ šµ0āˆ’š“0ī€øš‘‹ī€·š‘ šµš‘‡0+š“š‘‡0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øš‘‹ī€·š‘ šµš‘‡1+š“š‘‡1ī€ø,āˆ€š‘‹āˆˆā„‚š‘›Ć—š‘›.(2.7) By simple computations, we can getī€·š‘ šµ0āˆ’š“0ī€øš‘‹ī€·š‘ šµš‘‡0+š“š‘‡0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øš‘‹ī€·š‘ šµš‘‡1+š“š‘‡1ī€ø=š‘‹0ī€·š‘ šµš‘‡0+š“š‘‡0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øš‘Œ0,ī€·š‘ šµ0āˆ’š“0ī€øš‘Œī€·š‘ šµš‘‡0+š“š‘‡0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øš‘Œī€·š‘ šµš‘‡1+š“š‘‡1ī€ø=āˆ’š‘‹0ī€·š‘ šµš‘‡1+š“š‘‡1ī€ø+ī€·š‘ šµ0āˆ’š“0ī€øš‘Œ0.(2.8) Expressing by the operator language, that is,š“+š“āŽ›āŽœāŽœāŽš‘‹š‘ŒāŽžāŽŸāŽŸāŽ =āŽ›āŽœāŽœāŽāŽžāŽŸāŽŸāŽ šš²š‘‹šš²š‘Œ.(2.9)

In the following, we will convert matrix ordinary differential equation (2.3) to vector form. Let šœ‰ be the elementary transform, šœ‰āˆ¶ā„‚š‘›Ć—š‘›ā†’ā„‚š‘›2, that is,āŽ›āŽœāŽœāŽœāŽœāŽš‘Žšœ‰š“=š‘‡1ā‹®š‘Žš‘‡š‘›āŽžāŽŸāŽŸāŽŸāŽŸāŽ āŽ›āŽœāŽœāŽœāŽœāŽš‘Ž,āˆ€š“=1ā‹®š‘Žš‘›āŽžāŽŸāŽŸāŽŸāŽŸāŽ ,š‘Žš‘‡š‘–āˆˆā„‚š‘›,š‘–=1,2,ā€¦,š‘›.(2.10) Let ī€·š‘„=šœ‰š‘‹,š‘¦=šœ‰š‘Œ,š‘§=š‘„š‘¦ī€ø, andšø0=āŽ›āŽœāŽœāŽšµ0āŠ—š¼šµ1āŠ—š¼š¼āŠ—šµ1š¼āŠ—šµ0āŽžāŽŸāŽŸāŽ ,š¹0=āŽ›āŽœāŽœāŽš“0āŠ—š¼š“1āŠ—š¼āˆ’š¼āŠ—š“1āˆ’š¼āŠ—š“0āŽžāŽŸāŽŸāŽ .(2.11) Using the property of the Kronecker product, we havešœ‰(š“š‘‹šµ)=š“āŠ—šµš‘‡šœ‰š‘‹,(2.12) where š“,šµ,š‘‹āˆˆā„‚š‘›Ć—š‘›. So (2.3) can be written asšø0Ģ‡š‘§(š‘”)=š¹0š‘§(š‘”).(2.13) Similarly, by denoting š‘‡0=š‘‡0(š‘ ),š‘‡+0=š‘‡+0(š‘ ),Ī›0=Ī›0(š‘ ) as followsš‘‡0=š‘‡0āŽ›āŽœāŽœāŽī€·(š‘ )=š‘ šµ0āˆ’š“0ī€øī€·āŠ—š¼š‘ šµ1āˆ’š“1ī€øī€·āŠ—š¼š¼āŠ—š‘ šµ1+š“1ī€øī€·š¼āŠ—š‘ šµ0+š“0ī€øāŽžāŽŸāŽŸāŽ =š‘ šø0āˆ’š¹0,š‘‡+0=š‘‡+0(āŽ›āŽœāŽœāŽī€·š‘ )=š¼āŠ—š‘ šµ0+š“0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øī€·āŠ—š¼āˆ’š¼āŠ—š‘ šµ1+š“1ī€øī€·š‘ šµ0āˆ’š“0ī€øāŽžāŽŸāŽŸāŽ ,Ī›āŠ—š¼0=Ī›0ī€·(š‘ )=š‘ šµ0āˆ’š“0ī€øāŠ—ī€·š‘ šµ0+š“0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øāŠ—ī€·š‘ šµ1+š“1ī€ø,(2.14) we haveš‘‡+0š‘‡0āŽ›āŽœāŽœāŽš‘„š‘¦āŽžāŽŸāŽŸāŽ =āŽ›āŽœāŽœāŽĪ›0š‘„Ī›0š‘¦āŽžāŽŸāŽŸāŽ ,thatis,š‘‡+0š‘‡0=āŽ›āŽœāŽœāŽĪ›000Ī›0āŽžāŽŸāŽŸāŽ .(2.15)

Lemma 2.1. For all complex š‘ , detš‘‡+0(š‘ )=detš‘‡0(š‘ ), and so detš‘‡0(š‘ )=Ā±detĪ›0(š‘ ).

Proof. Let š‘‡0=ī€·š‘Žš‘š‘š‘‘ī€ø, where š‘Ž=(š‘ šµ0āˆ’š“0)āŠ—š¼, š‘=(š‘ šµ1āˆ’š“1)āŠ—š¼, š‘=š¼āŠ—(š‘ šµ1+š“1), and š‘‘=š¼āŠ—(š‘ šµ0+š“0). Noting the Kronecker product identities (š¼āŠ—š‘€)(š‘āŠ—š¼)=(š‘āŠ—š¼)(š¼āŠ—š‘€)=š‘āŠ—š‘€, we have š‘‘š‘=š‘š‘‘, š‘š‘=š‘š‘. By regularity of (šµ0,š“0), we know š‘ šµ0+š“0 is nonsingular for enough complex š‘ . So by the property of the polynomial and easy computations, we can get that detš‘‡+0(š‘ )=detš‘‡0(š‘ ) for any complex š‘ . Elsewhere, š‘‡+0š‘‡0=āŽ›āŽœāŽœāŽĪ›000Ī›0āŽžāŽŸāŽŸāŽ ,(2.16) so ī€·detš‘‡0ī€ø2=detš‘‡+0detš‘‡0=ī€·detĪ›0ī€ø2.(2.17) Thus detš‘‡0(š‘ )=Ā±detĪ›0(š‘ ).(2.18)

Theorem 2.2. Any imaginary axis eigenvalue of the systems (1.5) is a zero point of detĪ›0(š‘ ) and thus also one of the eigenvalues of the matrix pencil (šø0,š¹0).

Proof. The matrix polynomial for the system (1.5) is š‘(š‘ ,š‘’āˆ’š‘ šœ)=š‘ šµ0āˆ’š“0+ī€·š‘ šµ1āˆ’š“1ī€øš‘’āˆ’š‘ šœ.(2.19) Let š‘ =š‘–š‘¤ be an imaginary axis eigenvalue of the system (1.5) and š‘£ is associated eigenvector, ā€–š‘£ā€–=1. We have š‘(š‘ ,š‘’āˆ’š‘ šœ)š‘£=0. By conjugating and transforming, we can get ī€·š‘ šµ0āˆ’š“0ī€øš‘£š‘£āˆ—ī€·š‘ šµš‘‡0+š“š‘‡0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øš‘£š‘£āˆ—ī€·š‘ šµš‘‡1+š“š‘‡1ī€ø=0.(2.20) Via the elementary transform šœ‰, we get ī€ŗī€·š‘ šµ0āˆ’š“0ī€øāŠ—ī€·š‘ šµ0+š“0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øāŠ—ī€·š‘ šµ1+š“1šœ‰ī€·ī€øī€»š‘£š‘£āˆ—ī€ø=0,(2.21) that is, Ī›0(š‘ )š‘¢=0,š‘¢=šœ‰(š‘£š‘£āˆ—). We know that detĪ›0(š‘ )=0, and so ī€·detš‘ šø0āˆ’š¹0ī€ø=Ā±detĪ›0(š‘ )=0.(2.22)
From Theorem 2.2, we know that all of the imaginary axis eigenvalues of the system (1.5) are zero points of the algebraic equation detī€ŗī€·š‘ šµ0āˆ’š“0ī€øāŠ—ī€·š‘ šµ0+š“0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øāŠ—ī€·š‘ šµ1+š“1ī€øī€»=0.(2.23)

Corollary 2.3. If det(šµ0āŠ—šµ0āˆ’šµ1āŠ—šµ1)ā‰ 0, then šø0 is invertible, and any imaginary axis eigenvalue of the systems (1.5) is the eigenvalue of š¹0šø0āˆ’1.

Proof. By proof of Lemma 2.1, we have detšø0=||||||šµ0āŠ—š¼šµ1āŠ—š¼š¼āŠ—šµ1š¼āŠ—šµ0||||||ī€·šµ=det0āŠ—šµ0āˆ’šµ1āŠ—šµ1ī€ø.(2.24) The corollary follows immediately from Theorem 2.2.

Corollary 2.4. Any imaginary axis eigenvalue of the system with single delay Ģ‡š‘„(š‘”)=š“0š‘„(š‘”)+š“1š‘„(š‘”āˆ’šœ)(2.25) is an eigenvalue of š¹0=ī‚€š“0āŠ—š¼š“1āŠ—š¼āˆ’š¼āŠ—š“1āˆ’š¼āŠ—š“0ī‚.

Proof. It follows immediately from Theorem 2.2.

In fact, the above result contains the system (1.6). For the system (1.6), we also have the following corollary.

Corollary 2.5 (see Jarlebring and Hochstenbach [2, Theoremā€‰1]). For the system (1.6), the imaginary axis eigenvalues are the roots of the equation detī€ŗī€·š‘ š¼āˆ’š“0ī€øāŠ—ī€·š‘ š¼+š“0ī€øāˆ’ī€·š‘ šµ1āˆ’š“1ī€øāŠ—ī€·š‘ šµ1+š“1ī€øī€»=0.(2.26)

Remark 2.6. In fact, (2.23) or (2.26) is usually called a polynomial eigenvalue problem. The classical and most widely used approach to research the polynomial eigenvalue problems is linearization, where the polynomial is converted into a larger matrix pencil with the same eigenvalues. There are many forms for the linearization: the companion form is most typically commission. The linearization method is also an important tool to research the characteristic equations in algebraic methods, see [2, 13].

From the above results, we find that the imaginary axis eigenvalues of the system (1.5) or (1.6) can be computed via the algebraic equation (2.23) and (2.26). The imaginary eigenvalues play an important role in the stability. Next, we will use the results to discuss the stability of the systems. At first, we will give the condition of the delay-independent stability on the system (1.5). Secondly, we will address the problem of finding the critical delays of the system (1.6), that is, the delay such that the system (1.6) has purely imaginary eigenvalues.

3. The Algebraic Criteria of the Asymptotic Stability

The stability of the delay ordinary differential equations has been widely discussed [14, 15]. It is well known that the Lyapunov-Krasovskii functional approach is the important analytic method to find the delay-independent stability criteria, which do not include any information on the size of delay. The main ideas for developing algebraic criteria of the stability analysis on the systems can be found in many works, such as ā€œThe Degenerate Differential Systems with Delayā€ (W. Jiang, 1998, [12]). The results for singular neutral differential equations are still very few, especially by algebraic methods.

Next we first research the delay-independent stability of the system (1.5). The characteristic equation for the system(1.5) is denoted again byš‘ƒī€ŗš‘ ī€·šµ(š‘ ,š‘§)=det0+šµ1š‘§ī€øāˆ’ī€·š“0+š“1š‘§ī€øī€»,āˆ€š‘ āˆˆā„‚,š‘§=š‘’āˆ’š‘ šœ.(3.1)

From [12], we known that the solution of the neutral time-delay systems is asymptotically stable if all roots of (3.1) have negative real part bounded away from 0, that is, there exists a number š›æ>0, such that Re(š‘ )ā‰¤āˆ’š›æ<0 for any root š‘  of (3.1). Especially for the system (1.5), Zhu and Petzold had found that there must exist the š›æ, if the condition |š‘¢š‘‡šµ0š‘¢|ā‰„|š‘¢š‘‡šµ1š‘¢| for all š‘¢āˆˆā„š‘› holds, see [10]. So from [3, 10], we can get the following theorem.

Theorem 3.1. Let š‘  be the zeros of (2.23). If the coefficient matrices of the system (1.5) satisfies the following conditions:(i)Rešœ†<0,šœ†āˆˆšœŽ(šµ0,āˆ’š“0), (ii)maxReš‘ =0šœŒ((š‘ šµ0āˆ’š“0)āˆ’1(š‘ šµ1āˆ’š“1))<1,then the system (1.5) is asymptotically stable for all šœā‰„0, that is, the stability of the systems (1.5) is delay independent.

Proof. It follows immediately from [10].
In the following, we consider the neutral system (1.6), whose characteristic equation is denoted by š‘ƒī€ŗš‘ ī€·(š‘ ,š‘§)=detš¼+šµ1š‘§ī€øāˆ’ī€·š“0+š“1š‘§ī€øī€»,āˆ€š‘ āˆˆā„‚,š‘§=š‘’āˆ’š‘ šœ.(3.2) It is well known that the spectrum of the neutral delay systems exhibits some discontinuity properties, that is to say, an infinitesimal change of the delay parameter may cause the stability of the system to shift. These discontinuity properties are closely related to the essential spectrum of the system. The critical condition for a stability switch of a neutral delay system is that the rightmost eigenvalue goes from the left complex half-plane into the right complex half-plane by passing the imaginary axis. So the appearance of the imaginary axis eigenvalue is the critical condition. In Section 2, we find all of the imaginary axis eigenvalues. Let the delay be a parameter. In the following, we will find the critical value of the delay parameter such that the stability switch occurs. It is known that, if a neutral delay system is stable, it is necessary that its neutral part must be stable. For the system (1.6), this requirement concerns the stability of the difference equation: š‘„(š‘”)+šµ1š‘„(š‘”āˆ’šœ)=0.(3.3) The eigenvalues of (3.3) are called the essential spectrum of the system (1.6). We know that (3.3) is stable if šœŒ(šµ1)<1. It is important to point out that, under this assumption, the condition Re(š‘ )ā‰¤āˆ’š›æ<0 can be improved to Re(š‘ )<0. Our task is to find the critical delay where the system (1.6) becomes unstable. So we have the following theorem.

Theorem 3.2. Supposing all of the eigenvalues of matrix š“0 have negative real part and šœŒ(šµ1)<1, one has the following.(i)If for any root š‘  of (3.2), Re(š‘ )<0, else if for any root š‘ =š‘–š‘¤,š‘¤>0, š‘§=š‘’āˆ’š‘ šœ,|š‘§|<1, then system (1.6) is asymptotic stability for any šœā‰„0, that is, stability is delay independent. (ii) Otherwise, for each root š‘ =š‘–š‘¤,š‘¤>0, |š‘§|=1, one can get the minimal critical value of delay parameter šœāˆ—, such that if šœāˆˆ[0,šœāˆ—), then system (1.6) is asymptotic stability, and when šœā‰„šœāˆ—, the stability of system (1.6) changes, that is, the system (1.6) is delay-dependent stable and generates bifurcation in šœ=šœāˆ—.

Example 3.3. Consider a neutral neural networks with a single delay, see [16]:Ģ‡š‘„1(š‘”)=āˆ’š‘„1(ī€·š‘„š‘”)+š‘Žš‘“1(ī€øī€·š‘„š‘”āˆ’šœ)+š‘š‘”2(ī€øī€·š‘„š‘”āˆ’šœ)+š‘š‘”3(ī€øš‘”āˆ’šœ)+š‘Ž2š‘“2ī€·Ģ‡š‘„1ī€ø(š‘”āˆ’šœ)+š‘2š‘”2ī€·Ģ‡š‘„2ī€ø(š‘”āˆ’šœ)+š‘2š‘”2ī€·Ģ‡š‘„3ī€ø,(š‘”āˆ’šœ)Ģ‡š‘„2(š‘”)=āˆ’š‘„2ī€·š‘„(š‘”)+š‘š‘”1ī€øī€·š‘„(š‘”āˆ’šœ)+š‘Žš‘“2ī€øī€·š‘„(š‘”āˆ’šœ)+š‘š‘”3ī€ø(š‘”āˆ’šœ)+š‘2š‘”2ī€·Ģ‡š‘„1ī€ø(š‘”āˆ’šœ)+š‘Ž2š‘“2ī€·Ģ‡š‘„2ī€ø(š‘”āˆ’šœ)+š‘2š‘”2ī€·Ģ‡š‘„3ī€ø,(š‘”āˆ’šœ)Ģ‡š‘„3(š‘”)=āˆ’š‘„3ī€·š‘„(š‘”)+š‘š‘”2ī€øī€·š‘„(š‘”āˆ’šœ)+š‘Žš‘“3ī€ø(š‘”āˆ’šœ)+š‘2š‘”2ī€·Ģ‡š‘„2ī€ø(š‘”āˆ’šœ)+š‘Ž2š‘“2ī€·Ģ‡š‘„3ī€ø.(š‘”āˆ’šœ)(3.4) Now we rewrite system (3.4) as the matrix equations: Ģ‡Ģ‡š‘‹(š‘”)+šµš‘‹(š‘”āˆ’šœ)=š“š‘‹(š‘”)+š¶1š‘“(š‘‹(š‘”āˆ’šœ))+š¶2š‘”(š‘‹(š‘”āˆ’šœ)).(3.5) The linearization of the system (3.4) around the origin is given by Ģ‡š‘‹(š‘”)+šµ1Ģ‡š‘‹(š‘”āˆ’šœ)=š“0š‘‹(š‘”)+š“1š‘‹(š‘”āˆ’šœ),(3.6) where š“0=āŽ›āŽœāŽœāŽœāŽœāŽāŽžāŽŸāŽŸāŽŸāŽŸāŽ āˆ’1000āˆ’1000āˆ’1,š“1=āŽ›āŽœāŽœāŽœāŽœāŽāŽžāŽŸāŽŸāŽŸāŽŸāŽ š‘Žš‘š‘š‘š‘Žš‘0š‘š‘Ž,šµ1=āŽ›āŽœāŽœāŽœāŽœāŽš‘Ž2š‘2š‘2š‘2š‘Ž2š‘20š‘2š‘Ž2āŽžāŽŸāŽŸāŽŸāŽŸāŽ .(3.7) Assume that š‘Ž=āˆ’2, š‘=āˆ’1.5, š‘Ž2=0.3, š‘2=āˆ’0.3, and š‘“(š‘„)=š‘”(š‘„)=tanh(š‘„). We carry out the numerical simulations for system (3.4). From (2.26), by MATLAB computation, we can get the imaginary eigenvalue š‘ ā‰Ā±0.445š‘–,šœāˆ—ā‰0.45. From Theorem (3.2), we know that the zero solution of the system (3.4) is delay-dependent stable. The direction of the Hopf bifurcation at šœ=šœāˆ— is supercritical and the bifurcating periodic solutions are asymptotically stable. The simulation results are shown in Figures 1 and 2.

4. Conclusion

In this paper, we consider a singular neutral differential system with a single delay. Via applying the algebraic method, that is, the matrix pencil and the linear operators, we discussed the eigenvalues and the stability of the time-delay systems (1.5) and (1.6). By using MATLAB, we could easily compute imaginary eigenvalues from the algebraic equation (2.23) or (2.26). In fact, we only find the imaginary axis eigenvalues, which are the small part of the infinite eigenvalues. So compared with the analytic methods and the numerical methods, the algebraic methods are more simple and more explicit for some time-delay system. Certainly, applying the algebraic methods to analyze the dynamical properties of the singular neutral differential systems with delays is a new and immature field. So we believe that the algebraic methods used to research the stability of the dynamical systems would be more interesting in the future.


The National Natural Science Foundation of China (10871056) and the Fundamental Research Funds for the Central Universities (DL12BB24) are greatly acknowledged.


  1. J. K. Hale, E. F. Infante, and F. S. P. Tsen, ā€œStability in linear delay equations,ā€ Journal of Mathematical Analysis and Applications, vol. 105, no. 2, pp. 533ā€“555, 1985. View at: Publisher Site | Google Scholar
  2. E. Jarlebring and M. E. Hochstenbach, ā€œPolynomial two-parameter eigenvalue problems and matrix pencil methods for stability of delay-differential equations,ā€ Linear Algebra and its Applications, vol. 431, no. 3-4, pp. 369ā€“380, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  3. S. L. Campbell and V. H. Linh, ā€œStability criteria for differential-algebraic equations with multiple delays and their numerical solutions,ā€ Applied Mathematics and Computation, vol. 208, no. 2, pp. 397ā€“415, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  4. Z. Gao, ā€œPD observer parametrization design for descriptor systems,ā€ Journal of the Franklin Institute, vol. 342, no. 5, pp. 551ā€“564, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  5. Z. Gao and S. X. Ding, ā€œActuator fault robust estimation and fault-tolerant control for a class of nonlinear descriptor systems,ā€ Automatica, vol. 43, no. 5, pp. 912ā€“920, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  6. Z. Gao, T. Breikin, and H. Wang, ā€œReliable observer-based control against sensor failures for systems with time delays in both state and input,ā€ IEEE Transactions on Systems, Man, and Cybernetics A, vol. 38, no. 5, pp. 1018ā€“1029, 2008. View at: Publisher Site | Google Scholar
  7. T. E. Simos, ā€œClosed Newton-Cotes trigonometrically-fitted formulae of high order for long-time integration of orbital problems,ā€ Applied Mathematics Letters, vol. 22, no. 10, pp. 1616ā€“1621, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  8. S. Stavroyiannis and T. E. Simos, ā€œOptimization as a function of the phase-lag order of nonlinear explicit two-step P-stable method for linear periodic IVPs,ā€ Applied Numerical Mathematics, vol. 59, no. 10, pp. 2467ā€“2474, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  9. T. E. Simos, ā€œExponentially and trigonometrically fitted methods for the solution of the Schrƶdinger equation,ā€ Acta Applicandae Mathematicae, vol. 110, no. 3, pp. 1331ā€“1352, 2010. View at: Publisher Site | Google Scholar
  10. W. Zhu and L. R. Petzold, ā€œAsymptotic stability of linear delay differential-algebraic equations and numerical methods,ā€ Applied Numerical Mathematics, vol. 24, no. 2-3, pp. 247ā€“264, 1997. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  11. J. Louisell, ā€œA matrix method for determining the imaginary axis eigenvalues of a delay system,ā€ Institute of Electrical and Electronics Engineers, vol. 46, no. 12, pp. 2008ā€“2012, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  12. W. Jiang, The Degenerate Differential Systems with Delay, Anhui University, Hefei, China, 1998.
  13. D. S. Mackey, N. Mackey, C. Mehl, and V. Mehrmann, ā€œVector spaces of linearizations for matrix polynomials,ā€ SIAM Journal on Matrix Analysis and Applications, vol. 28, no. 4, pp. 971ā€“1004, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  14. W. Michiels and T. VyhlĆ­dal, ā€œAn eigenvalue based approach for the stabilization of linear time-delay systems of neutral type,ā€ Automatica, vol. 41, no. 6, pp. 991ā€“998, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  15. E. Jarlebring, ā€œOn critical delays for linear neutral delay systems,ā€ in Proceedings of the European Control Conference (ECC '07), Kos, Greece, 2007. View at: Google Scholar
  16. L. Li and Y. Yuan, ā€œDynamics in three cells with multiple time delays,ā€ Nonlinear Analysis. Real World Applications, vol. 9, no. 3, pp. 725ā€“746, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH

Copyright Ā© 2012 Jian Ma et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.