On asymptotic stability of a class of time–delay systems

where a is a real number, A is a 2 2 real constant matrix, and k, l are positive numbers such that k > l. Time–delay systems are a type of differential equations in which the derivative of the unknown function at a certain time is given in terms of the values of the function at previous times. Also, they are called delay differential equations, retarded differential equations or differential-difference equations. On the other hand, since asymptotic stability is an interdisciplinary material, the asymptotic stability of these systems has a wide range of applications as biology,


Introduction
In this paper, we study the asymptotic stability of the solutions of time-delay systems of the form where a is a real number, A is a 2 Â 2 real constant matrix, and k, l are positive numbers such that k > l. Time-delay systems are a type of differential equations in which the derivative of the unknown function at a certain time is given in terms of the values of the function at previous times. Also, they are called delay differential equations, retarded differential equations or differential-difference equations. On the other hand, since asymptotic stability is an interdisciplinary material, the asymptotic stability of these systems has a wide range of applications as biology,

ABOUT THE AUTHOR
My key research activities include Delay Differential Equations, Delay Difference Equations, Neutral Differential Equations, Neutral Difference Equations and stability of these equations. More generally, we can say Differential Equations, Difference Equations and stability of these equations. The research reported in this paper relates to the stability of the systems that can be modeled by delay differential equations such as biology, physics and medicine. Thus, it is useful to use qualitative approaches to investigate the asymptotic stability of these systems. Our research represents a generalized method for describing the asymptotic stability of systems within of science branches such as biology and physics.

PUBLIC INTEREST STATEMENT
In this paper, asymptotic stability of a class of time-delay systems was investigated. This system is a generalized version of modeling compound optical resonators by using a matrix instead of a scaler. In addition, since asymptotic stability is an interdisciplinary material, the asymptotic stability of these systems has a wide range of applications as biology, physics and medicine.
We know that for the time-delay equations, an equation is asymptotically stable if and only if all roots of the associated characteristic equation have negative real parts. Stability analysis, however, does not require the exact calculation of the characteristic roots. This analysis can be performed by D-subdivision method, which gives a necessary and sufficient condition for stability based on the coefficients or delay parameters of the characteristic equation. Thus, this method reveals a qualitative approach in order to prove the asymptotic stability of the systems. Consequently, we have created new necessary and sufficient conditions for the systems, which are asymptotically stable.
physics, and medicine. For instance, Naresh, Tripathi, Tchuenche, and Sharma (2009) formulated a nonlinear mathematical model to study the framework of disease Epidemiology. As another example Ikeda and Watanabe (2014) were investigated the stochastic differential equations and diffusion processes in their study about physic. Kruthika, Mahindrakar, and Pasumarthy (2017) were studied stability analysis of nonlinear time-delayed systems with application to biological models. They analyzed the local stability of a gene-regulatory network and immunotherapy for cancer modeled as nonlinear time-delay systems.

Preliminaries
It is known that for the time-delay equations, an equation is asymptotically stable if and only if all roots of the associated characteristic equation have negative real parts. Stability analysis, however, does not require the exact calculation of the characteristic roots; only the sign of the real part of the critical root must be determined. This analysis can be performed by D-subdivision method (see, e.g., Insperger & Stépán, 2011;Stepan, 1998), which gives a necessary and sufficient condition for stability based on the coefficients of the characteristic equation. The aim of this paper is to obtain new results for the asymptotic stability of zero solution of system (1.1), while the characteristic equation of system (1.1) has roots on the imaginary axis when A is a constant matrix. If we obtain x t ð Þ ¼ Py t ð Þ for a regular matrix P in (1.1), then we have the following system: Thus, matrix A can be given one of the following two matrices in Jordan form 7 ½ : ; q 1 ; q 2 and r are real constants; II ð Þ A ¼ q cos θ À sin θ sin θ cos θ ; q; θ are real constants and θ j j π 2 : Here we consider the case ðIIÞ; the other case should be considered similarly. The characteristic equation of system (1.1) is given as FðλÞ : ¼ det λI 2 þ 1 À a ð ÞI 2 þ A e Àλk þ e Àλl ¼ 0; (2:1) where I 2 is the 2 Â 2 identity matrix. Using 2:1 ð Þ; we obtain Þþq cos θ e Àλk þ e Àλl 2 þ q sin θ e Àλk þ e Àλl 2 ¼ λ þ 1 À a ð Þþq cos θ j j e Àλk þ e Àλl 2 À iq sin θ j j e Àλk þ e Àλl 2 ¼ λ þ 1 À a ð Þþqe i θ j j e Àλk þ e Àλl λ þ 1 À a ð Þþqe Ài θ j j e Àλk þ e Àλl : If we let (2:2) then we have where λ is the complex conjugate of any complex λ: Note that f λ ð Þ ¼ 0 implies f λ ð Þ ¼ 0:

Some auxiliary lemmas
In this section, we will investigate the distribution of the zeros of the characteristic equation of system (1.1). Thus, we state and prove some basic results on the roots of the characteristic equation of system (1.1).
Lemma 1 (Stepan, 1989) The zero solution of (1.1) is asymptotically stable if and only if all the roots of equation Lie in the left half of the complex plane.
Since f is an analytic function of λ; k and l for the fixed numbers a; q and θ; one can regard the root λ ¼ λ k; l ð Þ of (3.1) as a continuous function of k and l: The next lemma plays very important role for the main theorem.
Lemma 2 (Cooke & Grossman, 1982) As k and l vary, the sum of the multiplicities of the roots of (3.1) in the open right half-plane can change only if a root appears on or crosses the imaginary axis.
Consequently, we claim that (3.1) has only imaginary roots AE iω. We will determine how the value of k and l change as Equation (3.1) has roots on the imaginary axis. Now, we can write the characteristic Equation (3.1) as follows: Using the real part and the imaginary part of (3.3) which is equivalent to is obtained.
When q > 0; we have the following analogous result.
Lemma 5. Suppose that 0 < θ π 2 ; ω 2 Àπ kÀl ; π kÀl À Á À 0 f g and Also, the following conditions (3:11) are provided. Then all the roots of Equation (3.1) on the imaginary axis move in the right half-plane as k and l increase.

Main results
We will show that the stability analysis with a qualitative approach, as we have already mentioned in section 1.
Theorem 1. Suppose that 0 < θ π 2 and the conditions of Lemma 5 are satisfied. Let the matrix A of system (1.1) be in the form II ð Þ. Then system (1.1) is asymptotically stable if and only if either Proof. In the case of k ¼ 0 and l ¼ 0; the root of (3.1) is only λ 0; 0 ð Þ¼a À 1 À 2b cos θ À i2b sin θ: Thus, the root of the Equation (3.1) has a negative real part. By the continuity of the roots with respect to k and l; we can say that all the roots of (3.1) lie in the left half plane for k > 0 and l > 0 sufficiently small.
For the sufficiency, here our claim is: If either condition (4.1) or (4.2) holds, then (3.1) does not have a root on the imaginary axis. By condition (4.1) and Lemma 3, our claim is true for k > 0 and l > 0: Now, suppose that condition (4.2) holds: Since and Remark 1, we obtain k þ l ð Þ Þ k n þ l n ð Þ AE , k þ l ð ÞÞ α n þ β n ð Þ AE for n ¼ 0; 1; 2; . . . Thus, we obtain the contraposition with Lemma 1, our other claim is also true. By the above argument and Lemma 2, we can say that if either condition (4.1) or (4.2) holds, then all the roots of (3.1) lie in the left half plane.
For the necessity, we will show the following contraposition: either Thus, if (4.3) and (4.4) hold, then there exists roots λ j of (3.1) such that Re λ j À Á > 0 for j ¼ 1; 2: Assume that (4.3) holds and let λ 1 k; l ð Þ be the branch of the root of satisfying λ 1 0; 0 ð Þ¼a À 1 À 2b cos θ À i2b sin θ: Then, Lemma 5 or the continuity of λ 1 k; l ð Þ implies that Re λ 1 ð Þ > 0 for k > 0 and l > 0 sufficiently small. From here, we can say that λ 1 k; l ð Þ cannot move in the left half-plane crossing on the imaginary axis as k and l increase. Hence, we have Re λ 1 ð Þ> 0 for all k > 0 and l > 0: Assume that (4.4) holds and let λ 2 k; l ð Þ be the branch of the root of satisfying λ 2 k 0 ; l 0 ð Þ¼Àisgn q ð Þ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2q cos ω kÀl ð Þ 2 2 À a À 1 ð Þ 2 r : Then, Lemma 5 or the continuity of λ 2 k; l ð Þ implies that Re λ 2 ð Þ> 0 for k À k 0 > 0 and l À l 0 > 0 sufficiently small. From here, we can say that λ 2 k; l ð Þ cannot move in the left half-plane crossing on the imaginary axis as k and l increase. Hence, we have Re λ 2 ð Þ> 0 for all k > k 0 and l > l 0 : The proof is completed.
Remark 2. We consider the delay differential system (1.1) where matrix A is given as in case I ð Þ, i.e., Then, characteristic equation of (4.5) is as follows: λ þ 1 À a ð Þþq 1 e Àλk þ e Àλl λ þ 1 À a ð Þþq 2 e Àλk þ e Àλl ¼ 0: (4:6) It is obvious that for a ¼ 1 and i ¼ 1; 2, the equation λ þ q i e Àλk þ e Àλl À Á ¼ 0 is the characteristic equation of (4.6) with q ¼ q i , and so one can immediately obtain the following corollary from the previous result given by Kuang (1993).
Corollary 1. Suppose that a ¼ 1 for system (4.6) Let the matrix A of system (4.6) is written as the form I ð Þ: Then system (4.6) is asymptotically stable if and only if for i ¼ 1; 2 2q i k þ l ð Þcos k À l k þ l π 2 < π: (4:7) Theorem 2. Suppose that conditions of Lemma 5 are satisfied. Let the matrix A of system (1.1) be in the form I ð Þ: Proof. The proof is similar to Theorem 1.

An extension to a system of higher dimension
Finally, a higher dimensional linear delay differential system with two delays is considered where a is a real number, A is a d Â d real constant matrix, and k, l are positive numbers such that k > l.
Proof. Since q j e iθ j j ¼ 1; 2; . . . ; d ð Þbe the eigenvalues of matrix A′ the characteristic equation of system (5.1) is given by Þþq j e Àλkþi θ j j j þ e Àλlþi θ j j j ¼ 0: Thus, Theorem 3 can be seen as a result of Theorems 1 and 3.