This paper is devoted to providing a novel method to global exponential stability of impulsive delayed differential equations. By utilizing relative nonlinear measure method, several global exponential stability criteria are presented for the impulsive delayed differential equations. Compared with the Razumikhin technique and Lyapunov function method, our method is less conservative and gives a convergence rate, and one of our stability criteria is more flexible by incorporating an adjustable matrix. An example and its simulation are provided to illustrate that our method is efficient and our results are new and correct.
1. Introduction
Impulsive delayed differential equations are a kind of important mathematical models because they can reflect many real processes and phenomena appearing in biology, mechanics, medicine, electronics, economics, and so forth and have been successfully applied to these fields [1]. Exponential stability is a favorable quality for some applications, such as, the constructions and applications of artificial neural networks and the predictions of future economy [2]. Hence, exponential stability analysis of these equations has attracted interest in theoretical and practical applied fields [1, 3–7], which have presented some significant stability criteria. Among these existing methods, Lyapunov-Razumikhin method is undoubtedly the most popular (see [1, 3–7] and the relevant references therein). By analyzing their proofs, we find that some assumptions may be caused by Lyapunov-Razumikhin method itself. For example, the time delays are assumed to be smaller than the length of each impulsive interval [6] in order to employ Lyapunov-Razumikhin method although this assumption is proved to be unnecessary for the stability [7]. Moreover, the impulsive magnitude functions and the length of each impulsive interval are assumed to satisfy some inequality relations by means of v0 class functions (see, conditions (iii) and (iv) of Theorem 2.1 in [7] and Theorem 3.1 in [6]), and the supremum or infimum of lengths of impulsive intervals are supposed to meet the specific relations (see, Lemma 1 and Theorem 1 in [3]). These mean that these assumptions may be unnecessary from the point of the equation itself. Moreover, these unnecessary assumptions may bring more difficulty to the construction of Lyapunov function and further reduce the application range of the derived stability criteria. Therefore, it is significant to explore a novel method to discuss exponential stability of the impulsive delayed differential equations from the point of the equations themselves.
Being a nonlinear generalization of matrix measure (or logarithmic norm) of a matrix, relative nonlinear measure is presented [8] and employed to effectively characterize exponential stability of continuous delayed differential equations [2, 8]. Compared with the Lyapunov-Razumikhin method, relative nonlinear measure of coefficient operators of a differential equation can better and more essentially reflect the stability information of the differential equation because it does not need any additional assumptions on itself. Moreover, the relative nonlinear measure method has been to successfully characterize exponential stability of delayed and impulsive cellular neural networks [9]. It implies that the relative nonlinear measure method is also effective for describing exponential stability of a special kind of impulsive delayed differential equations. Motivated by this, the main aim of this paper is to utilize relative nonlinear measure to present several stability criteria for a general impulsive delayed differential equations from the viewpoint of theory.
This paper is organized as follows. In Section 2, we introduce some basic definitions and related knowledge. In Section 3, we present several criteria for exponential stability of the impulsive delayed differential equations, an illustrative example, and its simulation. Finally, conclusions are given in Section 4.
2. Preliminaries
Throughout this paper, the n-dimensional real vector space ℝn is endowed with l1-norm ∥·∥1; that is, ∥x∥1=∑i=1n|xi| for every x=(x1,x2,…,xn)T∈ℝn, where the superscript “T” denotes the transpose. Let ℝ+=[0,∞), J⊂ℝ and define PC(J,ℝn) by
(1)PC(J,ℝn)={x:J⟶ℝn∣xiscontinuous=everywhereexceptatthepointst=tk∈J=wherex(tk-)andx(tk+)existandx(tk)=x(tk-)}.
For τ>0, we equip the space PC([-τ,0],ℝn) with the norm ∥·∥ defined by
(2)∥ϕ∥=sup-τ≤s≤0∥ϕ(s)∥1,∀ϕ∈PC([-τ,0],ℝn).
In what follows, the symbol PC is used to denote PC([-τ,0],ℝn) and let ℕ={1,2,3,…}.
Firstly, we consider the following impulsive delayed differential equation:
(3)dx(t)dt=f(t,xt),t≥t0,t≠tk,Δx(tk)=Ik(tk,xtk),k∈ℕ,xt0=ϕ,
where functions f,Ik:ℝ+×PC→ℝn and ϕ∈PC satisfy all necessary conditions for the global existence, uniqueness, and continuality of solutions for t≥t0 [1, 10]; 0≤t0<t1<t2<⋯ is a strictly increasing sequences such that limk→+∞tk=+∞; Δx(tk)=x(tk+)-x(tk), x(tk+) denotes the right-hand limit of x(t) at t=tk; xt,xt+∈PC are defined by xt(s)=x(t+s) and xt+(s)=x(t++s) for all s∈[-τ,0]. For simplicity, we assume that f(t,0)=Ik(t,0)=0 for all t∈ℝ+ and k∈ℕ such that (3) admits the trivial solution. Denote by x(t)=x(t,t0,ϕ) the solution of (3) with the initial function xt0=ϕ. We further assume that all solutions x(t) of (3) are continuous except at tk,k∈ℕ, at which x(t) is left continuous; that is, x(tk)=x(tk-), k∈ℕ.
Definition 1.
Let f:ℝ+×PC→ℝn be a nonlinear operator. The constant
(4)ℳ(f,0)=supt∈ℝ+,xt∈PC,xt≠0〈f(t,xt)-f(t,0),sign(x(t))〉∥xt∥
is called relative nonlinear measure of f at xt=0. Here 〈·,·〉 denotes the inner product in ℝn and sign(x)=(sign(x1),sign(x2),…,sign(xn))T the sign vector of x∈ℝn, where sign(r) represents the sign function of r∈ℝ.
Definition 2.
Let f:ℝ+×PC→ℝn be a nonlinear operator. The constant
(5)ℒp(f,0)=supt∈ℝ+,xt∈PC,xt≠0∥f(t,xt)-f(t,0)∥1∥xt∥
is called minimal partial Lipschitz constant of f with respect to xt=0. The operator f is called partially Lipschitz continuous with respect to xt=0 if ℒp(f,0)<+∞.
Definition 3.
The trivial solution of (3) is said to be globally exponentially stable if there exist constants M≥1 and γ>0 such that
(6)∥x(t)∥1≤Me-γ(t-t0)∥ϕ∥
holds for t≥t0, where x(t) is the unique solution of (3) initiated from the function ϕ∈PC.
Lemma 4 (see [11]).
If a>c≥0, for every nonnegative real number b,
(7)0=λ-a+ceλb
has a unique positive solution.
Lemma 5 (see [12, Lemma D, Pages 389-390]).
Let γ and p be constants with 0<p<γ. Let v(·) be a continuous nonnegative function satisfying v′(t)≤-γv(t)+p∥vt∥ for t≥t0. Then one has
(8)v(t)≤∥vt0∥e-λ(t-t0)
for all t≥t0, where ∥vt0∥=supt0-b≤s≤t0∥v(s)∥ and λ is the unique positive solution of
(9)λ=γ-peλb.
Remark 6.
It is well known that an impulsive delayed differential equation is considered to be a continuous differential equation coupled with a difference equation. At present, the impulsive delayed differential equation is usually studied by main means of the research methods of the continuous delayed differential equation. Consequently, we can also directly apply Lemma 5 to characterize the asymptotic behaviors of the continuous parts of the impulsive delayed differential equation (3).
In what follows, we will frequently utilize the following relations between inner product 〈·,·〉 and l1-norm ∥·∥1 in ℝn:
(10)∥x∥1=〈x,sign(x)〉,∥x∥1≥〈x,sign(y)〉,∀x,y∈ℝn.
3. The Relative Nonlinear Measure Method
In order to discuss the exponential stability of the impulsive delayed differential equation (3), we make the following assumption.
The impulsive magnitude operator Ik:[t0,+∞)×PC→ℝn satisfies
(11)∥Ik(tk,xtk)∥1≤N∥x(tk)∥1
for k∈ℕ and some constant N>0.
Theorem 7.
Let the assumption (H) hold. If the relative nonlinear measure of f at xt=0 is negative, that is,
(12)ℳ(f,0)<0,
then the trivial solution of (3) is globally exponentially stable. That is, there exists a constant M≥1 such that the solution x(t) of (3) initiated from the function ϕ∈PC satisfies
(13)∥x(t)∥1≤Me-σ(t-t0)·∥ϕ∥
for all t≥t0, where σ is the unique positive solution of
(14)σ+ℳ(f,0)eστ=0.
Proof.
Let x(t) be any solution of the impulsive delayed differential equation (3) with the initial function xt0=ϕ∈PC. Let t∈[t0,t1)∪(tk,tk+1), k∈ℕ, and from the relations (10) we derive that
(15)∥x(t)∥1-∥x(t-s)∥1s≤1s〈x(t)-x(t-s),sign(x(t))〉
holds for all sufficiently small s>0 such that x(t) and x(t-s) simultaneously belong to some interval [t0,t1) or (tk0,tk0+1). Consequently, the function t↦∥x(t)∥1 is almost everywhere absolutely continuous in [t0,+∞), which implies that derivatives of ∥x(t)∥1 exist almost everywhere in [t0,+∞). Furthermore, along the solution trajectories of (3) we conclude that derivatives of ∥x(t)1∥ satisfy
(16)d∥x(t)∥1dt=〈dx(t)dt,sign(x(t))〉=〈f(t,xt),sign(x(t))〉=〈f(t,xt)-f(t,0),sign(x(t))〉≤ℳ(f,0)∥xt∥.
The combination of condition (12) and Lemmas 4 and 5 implies that
(17)∥x(t)∥1≤e-σ(t-t0)∥ϕ∥
holds for all t∈[t0,t1)∪(tk,tk+1), k∈ℕ, where σ is the unique positive solution of
(18)σ+ℳ(f,0)eστ=0.
According to the assumption (H) and Ik(tk,xtk)=x(tk+)-x(tk), we enjoy
(19)∥x(tk+)∥1=∥x(tk)+Ik(tk,xtk)∥1≤∥x(tk)∥1+∥Ik(tk,xtk)∥1≤(1+N)∥x(tk)∥1,k∈ℕ.
From (17) and (19) we derive
(20)∥x(tk+)∥1≤(1+N)∥x(tk)∥1=(1+N)∥x(tk-)∥1≤(1+N)e-σ(t-t0)∥ϕ∥,k∈ℕ.
Taking M=1+N, we can easily derive that (13) holds for all t≥t0.
Remark 8.
From the point of impulsive stabilization of unstable delayed differential equations, papers [6, 7] have discussed (3) by Lyapunov-Razumikhin method and have derived some significant stabilization conditions. In order to utilize Lyapunov-Razumikhin method, paper [6] assumed the time delays to be smaller than the length of each impulsive interval. Moreover, both of papers [6, 7] required the impulsive magnitude operator Ik and the length of each impulsive interval to satisfy special inequality relations by means of v0 class functions (see, conditions (iii) and (iv) of Theorem 2.1 in [7] and Theorem 3.1 in [6]). However, our method does not need these assumptions. Compared with results in [6, 7], hence, our results are new and may even be improvements of their results from the point of verifying the exponential stability of impulsive delayed differential equations.
Remark 9.
Paper [3] has studied exponential stability of the impulsive delayed differential equation (3) by Lyapunov-Razumikhin method and has presented several significant stability criteria with the assumption on the supremum or infimum of lengths of impulsive intervals satisfying special relations. Our method is not dependent on the assumption. Compared with results in [3], our result is new.
Equation (3) enjoys a general form of impulsive functional differential equations. The functional f sometimes acts on two separate parts: undelayed state x(t) and delayed one xt; that is,
(21)dx(t)dt=F(t,x(t))+G(t,xt),t≠tk,t≥t0,Δx(tk)=Ik(tk,xtk),k∈ℕ,xt0=ϕ,
where the nonlinear coefficient operators F and G satisfy assumptions of [1, 10] on their respective domains. Other assumptions of (21) are the same as ones of (3). For this special form, we can present the following corresponding stability criterion by means of the definition of the relative nonlinear measure of F at x=0∈ℝn.
Theorem 10.
Let the assumption (H) hold. F:ℝ+×ℝn→ℝn is a nonlinear operator and ℳ(F,0) denotes the relative nonlinear measure of F at x=0; that is,
(22)ℳ(F,0)=supt∈ℝ+,x(t)∈ℝn,x(t)≠0〈F(t,x(t))-F(t,0),sign(x(t))〉∥x(t)∥1.G:ℝ+×PC→ℝn is partially Lipschitz continuous with respect to xt=0. If
(23)ℳ(F,0)+ℒp(G,0)<0
holds, then the trivial solution of (21) is globally exponentially stable; that is, there exists a constant M≥1 such that
(24)∥x(t)∥1≤Me-σ(t-t0)∥ϕ∥
holds for all t≥t0, where σ is the unique positive solution of
(25)0=σ+ℳ(F,0)+ℒp(G,0)·eστ.
Proof.
Proof of the theorem is similar to that of Theorem 7. In detail, let x(t) be any solution of the impulsive delayed differential equation (21) with the initial function xt0=ϕ. Let t∈[t0,t1)∪(tk,tk+1),k∈ℕ, and from the relations (10) we derive that
(26)∥x(t)∥1-∥x(t-s)∥1s≤1s〈x(t)-x(t-s),sign(x(t))〉
holds for all sufficiently small s>0 such that x(t) and x(t-s) simultaneously belong to some interval [t0,t1) or (tk,tk+1). Consequently, the function t↦∥x(t)∥1 is almost everywhere absolutely continuous in [t0,+∞), which implies that derivatives of ∥x(t)∥1 exist almost everywhere in [t0,+∞). Furthermore, let t∈[t0,t1)∪(tk,tk+1), k∈ℕ, and from (21) we conclude that derivatives of ∥x(t)∥1 satisfy
(27)d∥x(t)∥1dt=〈dx(t)dt,sign(x(t))〉=〈F(t,x(t))+G(t,xt),sign(x(t))〉=〈F(t,x(t))-F(t,0),sign(x(t))〉+〈G(t,xt)-G(t,0),sign(x(t))〉≤〈F(t,x(t))-F(t,0),sign(x(t))〉+∥G(t,xt)-G(t,0)∥≤〈F(t,x(t))-F(t,0),sign(x(t))〉+ℒp(G,0)∥xt∥≤ℳ(F,0)∥x(t)∥1+ℒp(G,0)∥xt∥.
The combination of condition (23) and Lemmas 4 and 5 implies that
(28)∥x(t)∥1≤e-σ(t-t0)∥ϕ∥
holds for all t∈[t0,t1)∪(tk,tk+1), k∈ℕ, where σ is the unique positive solution of
(29)0=σ+ℳ(F,0)+ℒp(G,0)·eστ.
According to the assumption (H), we enjoy
(30)∥x(tk+)∥1=∥x(tk)+Ik(tk,xtk)∥1≤(1+N)∥x(tk)∥1,k∈ℕ.
From (28) and (30) we derive
(31)∥x(tk+)∥1≤∥(1+N)x(tk)∥1=∥(1+N)x(tk-)∥1≤(1+N)e-σ(t-t0)∥ϕ∥,k∈ℕ.
Consequently, we conclude that (24) holds for all t≥t0 if we take M=1+N.
In conclusion, Theorems 7 and 10 can be applied to discuss the exponential stability of a general nonlinear impulsive delayed differential equations. Obviously, they can be used to characterize the exponential stability of the following linear impulsive delayed differential equation:
(32)dx(t)dt=Ax(t)+Bx(t-τ(t)),t≥0,t≠tk,Δx(t)=Ck(x(t)),t=tk,k∈ℕ,x0=ϕ,
where matrices A,B:ℝn→ℝn and Ck satisfies (H), 0≤τ(t)≤τ.
It is worth being noted that relative nonlinear measure of a matrix A is just logarithmic norm (matrix measure) μ(A) of the matrix, and minimal partial Lipschitz constant of a matrix is just the matrix norm [13]. Based on this, combined with the technique of the adjustable matrix in [9], we can immediately derive the following result.
Corollary 11.
If there exists some diagonal matrix D=diag(d1,d2,…,dn) with di>0 such that the inequality
(33)μ1(AD)+∥BD∥1<0
holds, the trivial solution of (32) is globally exponentially stable; that is, there exists a constant M≥1 such that
(34)∥x(t)∥1≤Me-σt∥ϕ∥
holds for all t≥0, where σ is the unique positive solution of
(35)0=σ·min1≤i≤ndi+μ1(AD)+∥BD∥1·eστ.
Here μ1(AD) and ∥BD∥1 represent matrix measure of the matrix AD and matrix norm of the matrix BD with respect to ∥·∥1, respectively.
Remark 12.
According to definitions of matrix measure, matrix norm, and ∥·∥1, we can easily derive
(36)μ1(AD)=max1≤j≤n{djajj+∑i=1,i≠jndj|aij|},(37)∥BD∥1=max1≤j≤n{∑i=1ndj|bij|}.
Hence, condition (33) can be changed into the following form conveniently verified:
(38)max1≤j≤n{djajj+∑i=1,i≠jndj|aij|}+max1≤j≤n{∑i=1ndj|bij|}<0.
Example 13.
Consider the following linear impulsive delayed differential equation:
(39)x′(t)=Ax(t)+Bx(t-32(2+sin(t))),t≠tk,t≥0,Δx(t)=Ckx(t),t=tk,k∈ℕ,x0=ϕ,
where Ck=diag(-2,-5,-7),
(40)A=(-23-1.50-1510.158-3),B=(0.11.20.11-0.130.150.20.12-1.80.037).
Obviously, τ(t)=(3/2)(2+sin(t))≤4.5=τ and t-τ(t)=t-(3/2)(2+sin(t)) are not strictly increasing. Moreover, none of any restriction is imposed on lengths of the impulsive intervals. Consequently, the stability principles in [3, 6, 7] are unapplicable to deal with (39). Taking d1=2, d2=0.1, and d3=1, we can easily verify
(41)max1≤j≤n{djajj+∑i=1,i≠jndj|aij|}+max1≤j≤n{∑i=1ndj|bij|}=-0.5+0.475=-0.025<0.
From condition (38) we conclude that the trivial solution of the delayed impulsive equation (39) is globally exponentially stable with the convergence rate which is the unique solution of
(42)0.1σ-0.4+0.7e4.5σ=0.
The numerical simulation of this delayed impulsive equation is given in Figure 1 with respect to the initial functions:
(43)ϕ1(s)={0,s∈[-1,0),-2.3,s=0;ϕ2(s)={0,s∈[-1,0),1.5,s=0;ϕ3(s)={0,s∈[-1,0),-3.4,s=0.
Simulation of exponential convergence of (39).
4. Conclusions
The global exponential stability of the impulsive delayed differential equations has been discussed by means of the relative nonlinear measure. Several global exponential stability criteria have been provided for different forms of impulsive delayed differential equations. Compared with the existing methods, our method is less conservative and gives the convergence rate, and one of our stability criteria is more flexible by incorporating an adjustable matrix. An example and its simulation have been given to illustrate that our method is efficient and our results are new and correct.
Acknowledgments
This work was supported by the Natural Science Foundation of China under the Contacts nos. 11201038 and 11131006, Natural Science Foundation of Shaanxi Province under the Contact no. 2011JQ1004, and the Special Fund for Basic Scientific Research of Central Colleges under the Contact nos. 2013G2121017 and CHD2012TD015.
StamovaI.2009New York, NY, USAWalter de Gruyter10.1515/9783110221824MR2604930ZBL1189.34001SongX.PengJ.Exponential stability of equilibria of differential equations with time-dependent delay and non-Lipschitz nonlinearity20101153628363810.1016/j.nonrwa.2010.01.009MR2683818ZBL1206.34096ChenW. H.ZhengW. X.Exponential stability of nonlinear time-delay systems with delayed impulse effects20114751075108310.1016/j.automatica.2011.02.031MR2878380ZBL1233.93080LiuB.LiuX.TeoK. L.WangQ.Razumikhin-type theorems on exponential stability of impulsive delay systems2006711476110.1093/imamat/hxh091MR2203043ZBL1128.34047WangQ.LiuX.Exponential stability for impulsive delay differential equations by Razumikhin method2005309246247310.1016/j.jmaa.2004.09.016MR2154129ZBL1084.34066WangQ.LiuX.Impulsive stabilization of delay differential systems via the Lyapunov-Razumikhin method200720883984510.1016/j.aml.2006.08.016MR2323119ZBL1159.34347WuQ.ZhouJ.XiangL.Global exponential stability of impulsive differential equations with any time delays201023214314710.1016/j.aml.2009.09.001MR2559457ZBL1210.34105QiaoH.PengJ.XuZ. B.Nonlinear measures: a new approach to exponential stability analysis for Hopfield-type neural networks20011223603702-s2.0-003527237210.1109/72.914530SongX.XinX.HuangW.Exponential stability of delayed and impulsive cellular neural networks with partially Lipschitz continuous activation functions201229-3080902-s2.0-8486034160110.1016/j.neunet.2012.01.006BallingerG.LiuX.Existence and uniqueness results for impulsive delay differential equations199951–4579591MR1678279ZBL0955.34068PengJ.QiaoH.XuZ. B.A new approach to stability of neural networks with time-varying delays2002151951032-s2.0-0036125925DriverR. D.1977New York, NY, USASpringerMR0477368MoriT.FukumaN.KuwaharaM.On an estimate of the decay rate for stable linear delay systems1982361959710.1080/00207178208932877MR672704ZBL0483.93057