JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 929186 10.1155/2013/929186 929186 Research Article Stability Switches in a First-Order Complex Neutral Delay Equation Roales M. Rodríguez F. Saker Samir H. Departamento de Matemática Aplicada Universidad de Alicante Apartado Postal 99 03080 Alicante Spain ua.es 2013 27 2 2013 2013 23 07 2012 12 11 2012 2013 Copyright © 2013 M. Roales and F. Rodríguez. 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.

Stability of the first-order neutral delay equation x(t)+ax(tτ)=bx(t)+cx(tτ) with complex coefficients is studied, by analyzing the existence of stability switches.

1. Introduction

Delay differential equations (DDEs), and specifically DDEs of neutral type, appear in different scientific and technical problems, as in population dynamics, the modeling of networks containing lossless transmission lines, or the study of oscillations in elastic bars (see ).

The aim of this paper is to characterize the stability of the equation (1)x(t)+ax(t-τ)=bx(t)+cx(t-τ), where τ>0 is a constant delay and a, b, and c are complex parameters.

In , Cahlon and Schmidt characterized the asymptotic behavior of the zero solution of the retarded equation with complex coefficients (2)x(t)=px(t)+qx(t-τ), by transforming the complex equation into two coupled real DDEs. In , Wei and Zhang considered the same equation and, by studying the distribution of the roots of the characteristic equation for the associated real differential system with delay, analyzed the existence of stability switches .

Transforming a complex DDE into two coupled real DDEs to analyze its stability has some drawbacks, as, in general, the orders of the characteristic quasipolynomials to be analyzed double, and, since the study of the distribution of their roots is much more complicated as their degrees increase, it becomes very difficult to obtain necessary and sufficient conditions of stability [1, 9, 10].

To avoid this problem, recently Li et al.  presented a method for directly analyzing the stability of complex DDEs on the basis of stability switches. Their results generalize those for real DDEs, thus greatly reducing the complexity of the analysis. In this paper, the results developed in  will be used to study the stability switches of the zero solution of the neutral equation (1).

2. Previous Results

For the sake of clarity in the exposition, we recall some results that will be used later, that may be found in .

Consider time-delay systems with a characteristic equation in the following form: (3)Δ(λ,τ)=P(λ)+Q(λ)e-λτ, where τ0 and P(λ) and Q(λ) are complex polynomial of order n and m respectively, with either n>m or n=m and |α|>|β|, where α,β are the highest order coefficients of P(λ) and Q(λ), respectively. We will also assume that (4)P(0)+Q(0)0, which states that λ=0 is not a root of (3), and that P(λ) and Q(λ) have no roots on the imaginary axis simultaneously.

The assumed conditions are those required to use the result of Li et al.  stated below, but in other cases the stability of the system is mostly understood. Thus, if n<m, or n=m, and |α||β|, there are clusterings of roots on a vertical line in the right half plane or the imaginary axis, and in the first of these cases the zero solution is unstable for all delays. Also, if (4) fails or P(λ) and Q(λ) have common roots on the imaginary axis, these common roots, or λ=0, are roots of Δ(λ,τ) for all τ, and for all delays the zero solution is not asymptotically stable.

Consider the function (5)F(ω)=|P(iω)|2-|Q(iω)|2. If ω*0 is a zero of F(ω)=0, then there are an infinite number of delays τj corresponding to ω* satisfying (6)Δ(iω*,τj)=0. For these critical values, the following theorem [11, Theorem 1] characterizes the variation of the number of zeros with nonnegative real parts of Δ(λ,τ), in terms of the order and sign of the first nonzero derivative of F(ω) evaluated at ω=ω*, extending to the complex coefficients setting previous results valid only for real DDEs [8, 12, 13].

Theorem 1.

Assume that Δ(iω*,τj)=0,  j=0,1,2, Let N(τ) be the number of zeros with nonnegative real parts of Δ(λ,τ), and let M be an integer such that F(M)(ω*)0, and F(m)(ω*)=0 for all m<M. Then,

N(τ) keeps unchanged as τ increases along τj if M is even,

when M is odd, N(τ) increases by one if ω*F(M)(ω*)>0 and decreases by one if ω*F(M)(ω*)<0, as τ increases along τj.

It should be noted that it is possible for a particular delay τ to produce more than one pure imaginary root, so that, in Theorem 1, the change in N(τ) refers to the change resulting from the specific critical value ω*.

3. Stability Analysis of the First-Order Neutral DDE

Consider the complex DDE (1), where (7)a=a1+ia2,b=b1+ib2,c=c1+ic2. The characteristic equation associated with (1) is (8)λ+aλe-λτ-b-ce-λτ=0, so that (9)Δ(λ,τ)=P(λ)+Q(λ)e-λτ, where (10)P(λ)=λ-b,Q(λ)=aλ-c. Since the order of both polynomials is 1, we must demand that (11)|a|<1. As pointed out before, if |a|1 there are clusterings of roots on a vertical line in the right half plane, when |a|>1, in which case the system is unstable for all delays, or on the imaginary axis, when |a|=1. In this last critical case, the stability of the systems is not so clear, and a detailed, different type of analysis is required, as carried out in [8, pages 70–72] in the much simpler setting of real coefficients.

The following lemma gives N(0), the number of zeros with nonnegative real parts of Δ(λ,τ) when the delay is zero.

Lemma 2.

If (12)(b1+c1)(1+a1)+(b2+c2)a2<0, then N(0)=0. Otherwise, N(0)=1.

Proof.

Consider the equation (13)Δ(λ,0)=P(λ)+Q(λ)=λ-b+aλ-c=0. Then (14)λ=b+c1+a=(b1+c1)+i(b2+c2)1+a1+ia2=([(b1+c1)(1+a1)+(b2+c2)a2]+  i[(b2+c2)(1+a1)-(b1+c1)a2])×((1+a1)2+a22)-1, and therefore (15)(λ)<0(b1+c1)(1+a1)+(b2+c2)a2<0, and thus N(0)=0. If (λ)0, there is only one root with real part nonnegative, and hence N(0)=1.

Now consider the function (16)F(ω)=|P(iω)|2-|Q(iω)|2=(ω-b2)2+b12-(a1ω-c2)2-(a2ω+c1)2=ω2(1-a12-a22)+ω(2a1c2-2b2-2a2c1)+b12+b22-c12-c22=ω2(1-|a|2)-2(b+ac-)ω+|b|2-|c|2, and calculate its zeros. One gets (17)ω=(b+ac-)±B1-|a|2, where (18)B=((b+ac-))2-(1-|a|2)(|b|2-|c|2). First, we assume that (19)(b+ac-)>0, and consider several subcases.

If ((b+ac-))2-(1-|a|2)(|b|2-|c|2)<0, then F(ω) has no real root, and therefore the stability of the zero solution of (1) does not change for any τ>0.

If ((b+ac-))2-(1-|a|2)(|b|2-|c|2)>0, then F(ω) has two real roots, ω+, and ω-, such that (20)|b|2-|c|2>0ω+>ω->0,|b|2-|c|2<0ω+>0>ω-,|b|2-|c|2=0ω+>ω-=0, this last possibility being excluded from the analysis, since λ=0 is a root of (3), a contradiction to (4).

Consider the case where ω+>ω->0. Substituting λ=iω into (8), and separating the real and imaginary parts, one gets (21)ω-b2+(a2ω+c1)sinωτ+(a1ω-c2)cosωτ=0,-b1+(a1ω-c2)sinωτ-(a2ω+c1)cosωτ=0, obtaining the following two sets of values of τ for which there are roots, (22)τn,1=θ1ω++2nπω+,n=0,1,2,, where 0θ1<2π and (23)cosθ1=(b2-ω+)(a1ω+-c2)-b1(a2ω++c1)(a1ω+-c2)2+(a2ω++c1)2,sinθ1=(b2-ω+)(a2ω++c1)+b1(a1ω+-c2)(a1ω+-c2)2+(a2ω++c1)2,(24)τn,2=θ2ω-+2nπω-,n=0,1,2,, where 0θ2<2π and (25)cosθ2=(b2-ω-)(a1ω+-c2)-b1(a2ω-+c1)(a1ω--c2)2+(a2ω-+c1)2,sinθ2=(b2-ω-)(a2ω-+c1)+b1(a1ω--c2)(a1ω--c2)2+(a2ω-+c1)2.

Since (26)F(ω+)=-F(ω-)=B>0, one has (27)ω+F(ω+)>0,ω-F(ω-)<0. Therefore, according to Theorem 1, the number of the characteristic roots with nonnegative real parts increases by one as τ passes through τn,1 and decreases by one as τ passes through τn,2.

If N(0)=0, that is, if the zero solution of (1) is stable for τ=0, then it must follow that τ0,1<τ0,2, since N(τ) cannot become negative. There are stability switches when the delays are such that (28)τ0,1<τ0,2<τ1,1<τ1,2<. Since (29)τn+1,1-τn,1=2πω+<2πω-=τn+1,2-τn,2, the intervals become smaller with increasing n, so that eventually, for a certain k1, (30)τk-1,1<τk,1τk-1,2. Thus, the distribution of delays is (31)τ0,1<τ0,2<<τk-1,1<τk,1τk-1,2<τk+1,1<, and there is only a finite number of stability switches, with the system becoming unstable for τ>τk-1,1. Under these conditions, if N(0)=1, that is, if the zero solution of (1) is unstable for τ=0, the system cannot be stabilized.

If N(0)=1, stability switches occur when the distribution of delays is (32)τ0,2<τ0,1<<τk-1,2<τk-1,1<τk,1τk,2<τk+1,1<, in which case there are k switches from instability to stability to instability. Once τk,1τk,2, stability switch stops, so that the system becomes unstable for τ>τk-1,1.

The conditions on the delays for the above orderings to be valid can be formulated directly from (22) to (25).

Now we study the case when ω+>0>ω-. Proceeding as before, there are two sets of critical values of delays, τn,1 and τn,2, corresponding to ω+ and ω-, respectively, but now it holds that ω+F(ω+) and ω-F(ω-) are both positive. Hence, the number of the characteristic roots with nonnegative real parts increases by one as τ passes through τn,1 or τn,2. If N(0)=0, then there is τ0,1 such that the zero solution of (1) is stable for τ<τ0,1 and it becomes unstable for τ>τ0,1. If N(0)=1, the zero solution of (1) never becomes stable for any τ>0.

To finish the analysis when (b+ac-)>0, let (33)((b+ac-))2-(1-|a|2)(|b|2-|c|2)=0, so that F(ω) has a repeated real root, (34)ω*=2(b+ac-)1-|a|2>0. For this root, as in the previous cases, there is a set of critical delays τn. Since F(ω*)=0, we have to consider the second derivative, (35)F′′(ω*)=2(1-|a|2)0. According to Theorem 1, since M=2 is even, N(τ) keeps unchanged as τ increases along τn. We conclude that the stability of the zero solution of (1) does not change for any τ>0.

Finally, consider the case when (36)(b+ac-)<0. The analysis and results of this case are much similar to the previous one, and we will omit the details. There are only minor differences when (37)((b+ac-))2-(1-|a|2)(|b|2-|c|2)>0, and |b|2-|c|2>0.

In this subcase, the function F(ω) has two roots, 0>ω+>ω-, with corresponding sets of critical delays, (38)τn,1=θ1ω+-2nπω+,τn,2=θ2ω--2nπω-,n=0,1,2,, with -2πθ1<0 and -2πθ2<0 satisfying the expressions (23) and (25), respectively. It holds that ω+F(ω+)<0 and ω-F(ω-)>0, and hence, from Theorem 1, the number of the characteristic roots with nonnegative real parts decreases by one as τ passes through τn,1 and increases by one as τ passes through τn,2. Also, since (39)τn+1,2-τn,2=-2πω-<-2πω+=τn+1,1-τn,1, only a finite number of stability switches may exist.

If N(0)=0, then it must follow that τ0,2<τ0,1, and there are k switches from stability to instability when the distribution of the critical delays is (40)τ0,2<τ0,1<<τk-1,2<τk,2τk-1,1<τk+1,2<. Once τk,2τk-1,1, stability switch stops and the system becomes unstable for τ>τk-1,2.

Similarly, if N(0)=1, stability switches may occur for the distribution of critical delays (41)τ0,1<τ0,2<<τk-1,1<τk-1,2<τk,2τk,1<τk+1,2. Once τk,2τk,1, stability switch stops and the system remains unstable.

In summary, the following theorem has been established.

Theorem 3.

Under the hypotheses of Theorem 1, consider the first-order complex neutral delay equation (1), where a satisfies condition (11).

If (b+ac-)0, and ((b+ac-))2-(1-|a|2)(|b|2-|c|2)0,

if N(0)=0, then the zero solution of (1) is stable for all τ0;

if N(0)=1, then the zero solution of (1) is unstable for all τ0.

If (b+ac-)0,  ((b+ac-))2-(1-|a|2)(|b|2-|c|2)>0, and |b|2-|c|2<0,

if N(0)=0, then there is τ0,1 such that the zero solution of (1) is stable for τ<τ0,1 and it becomes unstable for τ>τ0,1;

if N(0)=1, the zero solution of (1) is always unstable.

If (b+ac-)>0, ((b+ac-))2-(1-|a|2)(|b|2-|c|2)>0, and |b|2-|c|2>0,

if N(0)=0 and the distribution of critical delays is (42)τ0,1<τ0,2<<τk-1,1<τk,1τk-1,2<τk+1,1<, then the zero solution of (1) is asymptotically stable for τ<τ0,1 and τj=0k-2(τj,2,τj+1,1) and unstable for τj=0k-2(τj,1,τj,2) and τ>τk-1,1;

if N(0)=1 and the distribution of critical delays is (43)τ0,2<τ0,1<<τk-1,2<τk-1,1<τk,1τk,2<τk+1,1<, then the zero solution of (1) is asymptotically stable for τj=0k-1(τj,2,τj,1) and unstable for τ<τ0,2,τj=0k-2(τj,1,τj+1,2), and τ>τk-1,1.

If (b+ac-)<0, ((b+ac-))2-(1-|a|2)(|b|2-|c|2)>0, and |b|2-|c|2>0,

if N(0)=0 and the distribution of critical delays is (44)τ0,2<τ0,1<<τk-1,2<τk,2τk-1,1<τk+1,2<, then the zero solution of (1) is asymptotically stable for τ<τ0,2 and τj=0k-2(τj,1,τj+1,2), and unstable for τj=0k-2(τj,2,τj,1) and τ>τk-1,2;

if N(0)=1 and the distribution of critical delays is (45)τ0,1<τ0,2<<τk-1,1<τk-1,2<τk,2τk,1<τk+1,2, then the zero solution of (1) is asymptotically stable for τj=0k-1(τj,1,τj,2) and unstable for τ<τ0,1, τj=0k-2(τj,2,τj+1,1), and τ>τk-1,2.

Finally, we will consider the case (46)(b+ac-)=0, for which (47)F(ω)=ω2(1-|a|2)+|b|2-|c|2.

If |b|2-|c|2>0, then F(ω) has no real root, and therefore the stability of the zero solution of (1) does not change for any τ0.

If |b|2-|c|2<0, then F(ω) has two real roots, (48)ω+=|c|2-|b|21-|a|2,ω-=-|c|2-|b|21-|a|2, with corresponding sets of critical delays (49)τn,1=θ1ω++2nπω+,n=0,1,2,, where 0θ1<2π, and (50)cosθ1=((b2-(|c|2-|b|2)/(1-|a|2))×(a1(|c|2-|b|2)/(1-|a|2)-c2))×(b12+((|c|2-|b|2)/(1-|a|2)-b2)2)-1-(b1(a2(|c|2-|b|2)/(1-|a|2)+c1))×(b12+((|c|2-|b|2)/(1-|a|2)-b2)2)-1,sinθ1=((b2-(|c|2-|b|2)/(1-|a|2))×(a2(|c|2-|b|2)/(1-|a|2)+c1))×(b12+((|c|2-|b|2)/(1-|a|2)-b2)2)-1+(b1(a1(|c|2-|b|2)/(1-|a|2)-c2))×(b12+((|c|2-|b|2)/(1-|a|2)-b2)2)-1,τn,2=θ2ω--2nπω-,n=0,1,2,, where -2πθ2<0, and (51)cosθ2=((b2+(|c|2-|b|2)/(1-|a|2))×(-a1(|c|2-|b|2)/(1-|a|2)-c2))×(b12+((|c|2-|b|2)/(1-|a|2)+b2)2)-1-(b1(-a2(|c|2-|b|2)/(1-|a|2)+c1))×(b12+((|c|2-|b|2)/(1-|a|2)+b2)2)-1,sinθ2=((b2+(|c|2-|b|2)/(1-|a|2))×(-a2(|c|2-|b|2)/(1-|a|2)+c1))×(b12+((|c|2-|b|2)/(1-|a|2)+b2)2)-1+(b1(-a1(|c|2-|b|2)/(1-|a|2)-c2))×(b12+((|c|2-|b|2)/(1-|a|2)+b2)2)-1.

In this case, it holds that ω+F(ω+) and ω-F(ω-) are both positive, and hence, from Theorem 1, the number of the characteristic roots with nonnegative real parts increases by one as τ passes either through τn,1 or τn,2. Thus, the following theorem has been established.

Theorem 4.

Under the hypotheses of Theorem 1, consider (1) satisfying condition (11). If (b+ac-)=0,

if |b|2-|c|2>0, then the zero solution of (1) is always stable if N(0)=0 and unstable if N(0)=1;

if |b|2-|c|2<0, if N(0)=0 then there is τ0,1 such that the zero solution of (1) is stable for τ<τ0,1 and it becomes unstable for τ>τ0,1. If N(0)=1, the zero solution of (1) never becomes stable for any τ0.

Acknowledgments

The authors wish to thank the reviewers for useful comments and suggestions, which improved the paper. F. Rodríguez is grateful for financial support from project FEEDBACK (CGL2011-30515-C02-01), funded by the Spanish Ministry of Science and Innovation.

Gopalsamy K. Stability and Oscillations in Delay Differential Equations of Population Dynamics 1992 Dordrecht, The Netherlands Kluwer 1163190 ZBL0833.35144 Driver R. D. A mixed neutral system Nonlinear Analysis: Theory, Methods and Applications A 1984 8 2 155 158 10.1016/0362-546X(84)90066-X 734448 ZBL0553.34042 Brayton R. K. Willoughby R. A. On the numerical integration of a symmetric system of difference-differential equations of neutral type Journal of Mathematical Analysis and Applications 1967 18 182 189 0213061 10.1016/0022-247X(67)90191-6 ZBL0155.47302 Cahlon B. Schmidt D. On stability of a first-order complex delay differential equation Nonlinear Analysis: Real World Applications 2002 3 3 413 429 10.1016/S1468-1218(01)00039-6 1905459 ZBL1010.34073 Wei J. Zhang C. Stability analysis in a first-order complex differential equations with delay Nonlinear Analysis: Theory, Methods and Applications A 2004 59 5 657 671 10.1016/j.na.2004.07.027 2096322 ZBL1067.34082 Cooke K. L. Grossman Z. Discrete delay, distributed delay and stability switches Journal of Mathematical Analysis and Applications 1982 86 2 592 627 10.1016/0022-247X(82)90243-8 652197 ZBL0492.34064 Beretta E. Kuang Y. Geometric stability switch criteria in delay differential systems with delay dependent parameters SIAM Journal on Mathematical Analysis 2002 33 5 1144 1165 10.1137/S0036141000376086 1897706 ZBL1013.92034 Kuang Y. Delay Differential Equations with Applications in Population Dynamics 1993 Arizona, Ariz, USA Academic Press 1218880 Ruan S. Wei J. On the zeros of a third degree exponential polynomial with applications to a delayed model for the control of testosterone secretion IMA Journal of Mathematics Applied in Medicine and Biology 2001 18 41 52 Li X. Wei J. On the zeros of a fourth degree exponential polynomial with applications to a neural network model with delays Chaos, Solitons and Fractals 2005 26 2 519 526 10.1016/j.chaos.2005.01.019 2143965 ZBL1098.37070 Li J. Zhang L. Wang Z. Two effective stability criteria for linear time-delay systems with complex coefficients Journal of Systems Science and Complexity 2011 24 5 835 849 10.1007/s11424-011-9252-4 2853684 Lee M. S. Hsu C. S. On the τ-decomposition method of stability analysis for retarded dynamical systems SIAM Journal on Control and Optimization 1969 7 215 233 Cooke K. L. van den Driessche P. On zeroes of some transcendental equations Funkcialaj Ekvacioj 1986 29 1 77 90 865215 ZBL0603.34069