`Journal of Applied MathematicsVolume 2012, Article ID 749847, 11 pageshttp://dx.doi.org/10.1155/2012/749847`
Research Article

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

1Department of Mathematics, Harbin Institute of Technology, Harbin 150001, China
2Department of Mathematics, Northeast Forestry University, Harbin 150040, China

Received 29 November 2011; Revised 6 April 2012; Accepted 8 April 2012

Academic Editor: Zhiwei Gao

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.

#### Abstract

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 . 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 , the numerical techniques for the computation of the eigenvalues were discussed. In , 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: where , , and and are given by continuous functions on the initial interval . So we have the solution: 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 where , , , , and is given by continuous functions on the initial interval . When , we called it the retarded differential equations. It can be improved as Many scholars have widely researched the delay-independent or delay-dependent stability and asymptotic stability by analytic methods or numerical methods. When , it is called a singular (or degenerated) delay-differential equations. The imaginary axis eigenvalues are discussed by using matrix pencil, see . 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: where , , . For the system (1.5), if , we can improve it as the form 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 . Here, we consider the state . 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 is said to be regular if is not identically singular for any complex . If is regular, the zero of is called the eigenvalue of the matrix pencil . From , we know that the system (1.5) is solvable if and only if is regular. So, in this paper, we suppose that 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 where , . Let denote the vector space and denote the operators on , given by With , the system (2.1) can be written as Supposing is a matrix solution of the system (2.1), we have For any complex , let be the operator ; then For any complex , let , satisfying For any complex , let , satisfying By simple computations, we can get Expressing by the operator language, that is,

In the following, we will convert matrix ordinary differential equation (2.3) to vector form. Let be the elementary transform, , that is, Let , and Using the property of the Kronecker product, we have where . So (2.3) can be written as Similarly, by denoting as follows we have

Lemma 2.1. For all complex , , and so .

Proof. Let , where , , , and . Noting the Kronecker product identities , we have , . By regularity of , we know is nonsingular for enough complex . So by the property of the polynomial and easy computations, we can get that for any complex . Elsewhere, so Thus

Theorem 2.2. Any imaginary axis eigenvalue of the systems (1.5) is a zero point of and thus also one of the eigenvalues of the matrix pencil .

Proof. The matrix polynomial for the system (1.5) is Let be an imaginary axis eigenvalue of the system (1.5) and is associated eigenvector, . We have . By conjugating and transforming, we can get Via the elementary transform , we get that is, . We know that , and so
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

Corollary 2.3. If , then is invertible, and any imaginary axis eigenvalue of the systems (1.5) is the eigenvalue of .

Proof. By proof of Lemma 2.1, we have The corollary follows immediately from Theorem 2.2.

Corollary 2.4. Any imaginary axis eigenvalue of the system with single delay is an eigenvalue of .

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

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, ). 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

From , 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 , such that 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 for all holds, see . 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), (ii),then the system (1.5) is asymptotically stable for all , that is, the stability of the systems (1.5) is delay independent.

Proof. It follows immediately from .
In the following, we consider the neutral system (1.6), whose characteristic equation is denoted by 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: The eigenvalues of (3.3) are called the essential spectrum of the system (1.6). We know that (3.3) is stable if . It is important to point out that, under this assumption, the condition can be improved to . 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 have negative real part and , one has the following.(i)If for any root of (3.2), , else if for any root , , then system (1.6) is asymptotic stability for any , that is, stability is delay independent. (ii) Otherwise, for each root , , one can get the minimal critical value of delay parameter , such that if , 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 : Now we rewrite system (3.4) as the matrix equations: The linearization of the system (3.4) around the origin is given by where Assume that , , , , and . We carry out the numerical simulations for system (3.4). From (2.26), by MATLAB computation, we can get the imaginary eigenvalue . 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.

Figure 1: For system (3.4), when , the equilibrium is asymptotically stable.
Figure 2: For system (3.4), when , the periodic solution bifurcates from the equilibrium.

#### 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.

#### Acknowledgments

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

#### References

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.
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.
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.
4. Z. Gao, “PD observer parametrization design for descriptor systems,” Journal of the Franklin Institute, vol. 342, no. 5, pp. 551–564, 2005.
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.
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.
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.
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.
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.
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.
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.
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.
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.
15. E. Jarlebring, “On critical delays for linear neutral delay systems,” in Proceedings of the European Control Conference (ECC '07), Kos, Greece, 2007.
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.