Next Article in Journal
On Ball Numbers
Previous Article in Journal
Differential Equations Arising from the Generating Function of the (r, β)-Bell Polynomials and Distribution of Zeros of Equations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Approach for Exponential Stability Criteria of New Certain Nonlinear Neutral Differential Equations with Mixed Time-Varying Delays

1
Department of Mathematics, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
2
Department of Mathematics, Faculty of Science, University of Phayao, Phayao 56000, Thailand
3
Department of Mathematics, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(8), 737; https://doi.org/10.3390/math7080737
Submission received: 15 July 2019 / Revised: 4 August 2019 / Accepted: 9 August 2019 / Published: 12 August 2019
(This article belongs to the Section Mathematical Biology)

Abstract

:
This work is concerned with the delay-dependent criteria for exponential stability analysis of neutral differential equation with a more generally interval-distributed and discrete time-varying delays. By using a novel Lyapunov–Krasovkii functional, descriptor model transformation, utilization of the Newton–Leibniz formula, and the zero equation, the criteria for exponential stability are in the form of linear matrix inequalities (LMIs). Finally, we present the effectiveness of the theoretical results in numerical examples to show less conservative conditions than the others in the literature.

1. Introduction

Many scientists have worked on the stability criteria for the problem of neutral delayed differential equation (NDE) with constant delays [1,2,3,4,5,6,7,8,9] as follows:
d d t [ x ( t ) + p x ( t τ ) ] + a x ( t ) b tanh x ( t σ ) = 0 , t 0 ,
where | p | < 1 and a , τ , b, σ are real positive constants. With each solution x ( t ) of (1), the initial condition is satisfied:
x ( s ) = ϕ ( s ) , s [ r , 0 ] , where ϕ C ( [ r , 0 ] , R ) , r = max { τ , σ } .
The NDE (1) as a special case has been considered by many researchers for the study of the dynamical characteristics of the neural networks of the Hopfield type; see [10,11,12,13] and the references cited therein. As we all know, a small change in delay may destabilize a system [3,5,6,7,9]. Hence, many researchers increased their study of the stabilization of the system by proposing the criteria in several forms.
In [6], the delay-dependent stability criterion of certain neutral differential equations was considered by the integral inequalities approaches. However, the sufficient condition only depends on the discrete delay σ and does not depend on the neutral delay τ . An improved stability criterion for a class of neutral differential equations was studied in [8] by introducing a new Lyapunov functional and avoiding the use of the stability assumption on the main operators. In [1], the authors considered using a Lyapunov–Krasovskii functional technique and a descriptor model transformation approach to transform the system with distributed and discrete delays. The issue of exponential stability analysis of the neutral differential equation as in (1) has been presented [4,5] by utilizing the descriptor form, Lyapunov–Krasovskii functionals, and model transformation.
In the past few years, many researches have studied several methods in time-varying delay, s which can be considered as closer way to reflect the realistic situation because the time delays are usually not constant. In [1,14,15], the authors considered the problem of NDE for exponential stability with time-varying delays as follows:
d d t [ x ( t ) + p x ( t τ ( t ) ) ] + a x ( t ) b tanh x ( t σ ( t ) ) = 0 , t 0 ,
where | p | < 1 and a , b are real positive constants and σ ( t ) and τ ( t ) are discrete and neutral time-varying delays, respectively,
0 τ ( t ) τ , τ ˙ ( t ) < τ d , 0 σ ( t ) σ , σ ˙ ( t ) < σ d ,
where τ and σ are given real positive constants. For each solution x ( t ) of (2), the initial condition is satisfied:
x 0 ( t ) = ϕ ( t ) , where ϕ C ( [ r , 0 ] , R ) , r = max { τ , σ } .
In [14], the sufficient conditions were presented with the use of descriptor form, Lyapunov functions, the Leibniz–Newton formula, and radially-unbounded function.
On the other hand, distributed delays, as one of the most common types, have received considerable attention in recent years. In a practical system, the distributed delay is often useful in modeling processes that are (i) described by flow rates (aggregates) of entities that move at different rates through the given irreversible process (ii). Recently, a certain neutral differential equation with distributed and discrete time-varying delays was considered in [16] as follows:
d d t [ x ( t ) + p x ( t τ ( t ) ) ] + a x ( t ) b tanh x ( t σ ( t ) ) c t ρ ( t ) t x ( s ) d s = 0 , t 0 ,
where | p | < 1 and a , b , c are real positive constants. σ ( t ) , ρ ( t ) , and τ ( t ) are discrete, distributed, and neutral time-varying delays, respectively.
In practice, distributed time delays arise in models for a variety of applications including population dynamic [17], the predator-prey model [18,19], the neuronal model [20,21], neural networks [22,23,24,25], and a nonlinear system [26], which considered the distributed delay of the activation function [24]. From the discussions above, there still exists room for further improvement in NDE, and it might occur that there are connection between the distributed time-varying delay and tanh x ( s ) . As a result, it is more natural to consider NDE with the distributed delay of tanh x ( s ) . However, in recent works, there were no models [1,2,3,4,5,6,7,8,9,14,15,27,28,29,30] that include distributed time-varying delays.
Our main contributions can be summarized as the following aspects. First, our work is the first attempt to study for the problem of new certain nonlinear neutral differential equation for exponential stability with a more generally interval-distributed and discrete time-varying delays. The Lyapunov–Krasovkii functional, descriptor model transformation, the Newton–Leibniz formula, and utilization of the zero equation are applied. The exponential stability criteria are in the form of linear matrix inequalities (LMIs). Second, note that most of the existing NDE results [1,9,14,15,16] cannot be applied to the situation when there are two distributed time-varying delays. Finally, a comparison of the maximum delay bounds obtained by our proposed methods with those in many existing literature works [1,9,14,15] is conducted, which supports that the exponential stability criteria provided in this paper can enhance the feasible region effectively and significantly.
The rest of this paper is arranged as follows. Section 2 shows the definitions, notations, and lemmas for applying the main results of the paper. The exponential stability criteria for delay-dependent sufficient conditions are shown in Section 3. Numerical examples presented for effectiveness obtained consequences that are manifested in Section 4. This work ends with the conclusion in Section 5.
Notation: The following notations will be used throughout this manuscript. B T denotes the transpose of matrix B; B 0 ( B > 0 ) means that B is semi-positive definite (positive definite); · denotes the Euclidean norm of the vector or matrix, x s = sup r s 0 x ( s ) and x ˙ l = sup r s 0 x ˙ ( s ) l . The symbol * stands for the transposed elements in the symmetric positions.

2. Problem Formulation and Preliminary

Consider the neutral differential equation with time-varying delays as follows:
d d t [ x ( t ) + p x ( t τ ( t ) ) ] + a x ( t ) b tanh x ( t σ ( t ) ) e t h ( t ) t x ( s ) d s f t l ( t ) t tanh x ( s ) d s = 0 ,
where | p | < 1 and a , b , e, f are real positive constants. h ( t ) , l ( t ) , and σ ( t ) are distributed and discrete interval time-varying delays, and τ ( t ) is a neutral interval time-varying delay, satisfying,
0 τ 1 τ ( t ) τ 2 , τ ˙ ( t ) τ d ,
0 σ 1 σ ( t ) σ 2 , σ ˙ ( t ) σ d ,
0 h 1 h ( t ) h 2 ,
0 l 1 l ( t ) l 2 ,
where τ 1 , τ 2 , σ 1 , σ 2 , h 1 , h 2 , l 1 , l 2 , τ d , and σ d are given positive real constants. For each solution x ( t ) of (4), the initial condition is satisfied:
x 0 ( t ) = ϕ ( t ) , t [ r , 0 ] ,
where ϕ C ( [ r , 0 ] ; R ) and r = max { τ 2 , σ 2 , h 2 , l 2 } .
Definition 1
([14]). If there exist real positive constants λ, K such that:
x ( t ) K e λ t sup r s 0 x ( s ) = K e λ t x s ,
then the equilibrium point of Equation ( 2 ) is exponentially stable.
Lemma 1
(Jensen’s inequality [31]).For positive real number h, any symmetric positive definite matrix Q and vector function x ˙ ( t ) : [ h , 0 ] R n such that the following inequality holds:
h h 0 x ˙ T ( s + t ) Q x ˙ ( s + t ) d s + h 0 x ˙ ( s + t ) d s T Q h 0 x ˙ ( s + t ) d s 0 .
Lemma 2
(Cauchy’s inequality [31]).If there exists constant symmetric positive definite matrix P R n × n and a , b R n , we have:
± 2 a T b a T P a + b T P 1 b .
Lemma 3
([16]). For any constant Q R n × n is a symmetric positive definite matrix, h ( t ) is a discrete time-varying delay with (7), and vector function ω : [ h 2 , 0 ] R n such that the following inequality holds:
[ h 2 h 1 ] h 2 h 1 ω T ( s ) Q ω ( s ) d s h ( t ) h 1 ω T ( s ) d s Q h ( t ) h 1 ω ( s ) d s h 2 h ( t ) ω T ( s ) d s Q h 2 h ( t ) ω ( s ) d s .
Lemma 4
(Peng–Park’s integral inequality [32]).For positive scalars τ and τ ( t ) satisfying 0 < τ ( t ) < τ , vector function x ˙ ( s ) : [ τ , 0 ] R n , any matrix Z R n × n , Z = Z T , matrix S R n × n , and any matrix Z S * Z 0 such that the following concerned inequality holds:
τ t τ t x ˙ T ( s ) Z x ˙ ( s ) d s ω T ( t ) ω ( t )
where ω = [ x T ( t ) , x T ( t τ ( t ) ) , x T ( t τ ) ] and = Z Z S S * 2 Z + S + S T Z S * * Z .
Remark 1.
It should be pointed out that the problem of the delay differential equation with distributed delay function of tanh x ( s ) is more general than [16], where the delay is not considered. Generally applicable nonlinear system contain with mixed time-varying delays, for instance the population dynamics and predator-prey model that were studied in [17,18,19]. It should be noted that this is the first encountered when the distributed time-varying delay function of tanh x ( s ) is studied, which includes both distributed and discrete delays for the system.

3. Stability Analysis

For this part, we analyze the issue of the exponential stability criteria for the neutral differential Equation (4) with time-varying delays. From the model transformation technique, we get the Newton–Leibniz formula as follows:
t τ ( t ) t x ˙ ( s ) d s = x ( t ) x ( t τ ( t ) ) .
By utilizing the zero equation, we obtain:
d 1 t τ ( t ) t x ˙ ( s ) d s = d 1 x ( t ) d 1 x ( t τ ( t ) ) ,
( 1 d 2 ) t τ ( t ) t x ˙ ( s ) d s = ( 1 d 2 ) x ( t ) ( 1 d 2 ) x ( t τ ( t ) ) ,
where d 1 , d 2 R will be selected to guarantee the exponential stability of Equation (4). By (9)–(11), Equation (4) can be formulated as follows:
x ˙ ( t ) = ( d 1 a ) x ( t ) + b tanh x ( t σ ( t ) ) ( p + τ ˙ ( t ) d 2 τ ˙ ( t ) ) x ˙ ( t τ ( t ) ) d 1 x ( t τ ( t ) ) d 1 t τ ( t ) t x ˙ ( s ) d s + e t h ( t ) t x ( s ) d s + f t l ( t ) t tanh x ( s ) d s .
The following theorem presents the exponential stability for the neutral differential equation with time-varying delays in Equation (4).
Theorem 1.
Consider d 1 , d 2 R and positive real constants σ 1 , σ 2 , σ d , τ 1 , τ 2 , τ d , h 1 , h 2 , l 1 , l 2 . Equation (4) is exponentially stable with a delay rate α > 0 . For any positive real constants β i , η i , i = 1 , 2 , , 10 and real constants w k , y k , k = 1 , 2 , , 10 such that:
Ξ < 0 ,
where Ξ = [ Φ ( i , j ) ] and Φ ( j , i ) = Φ ( i , j ) ,
Φ ( 1 , 1 ) = 2 α β 1 + 2 β 1 d 1 2 β 1 a + 2 w 1 d 1 2 w 1 a + β 2 + β 3 + η 1 + η 2 + β 6 τ 2 2 + β 7 σ 2 2 + η 3 τ 1 2 + η 4 σ 1 2 β 8 + β 9 + β 10 σ 2 2 + η 6 σ 1 2 + 2 y 2 + η 7 h 2 2 + η 8 l 2 2 + η 9 h 1 2 + η 10 l 1 2
Φ ( 1 , 2 ) = w 1 + w 2 d 1 w 2 a + y 1 ,
Φ ( 1 , 3 ) = β 1 b + w 1 b + w 3 d 1 w 3 a + y 5 ,
Φ ( 1 , 4 ) = β 1 d 1 w 1 d 1 + β 8 s y 2 + y 3 ,
Φ ( 1 , 5 ) = β 1 p β 1 τ d + β 1 d 2 τ d w 1 p w 1 τ d + w 1 d 2 τ d + w 4 d 1 w 4 a ,
Φ ( 1 , 6 ) = β 1 d 1 w 1 d 1 + w 5 d 1 w 5 a y 2 + y 4 ,
Φ ( 1 , 7 ) = β 1 e + w 1 e ,
Φ ( 1 , 8 ) = β 1 f + w 1 f ,
Φ ( 1 , 9 ) = s ,
Φ ( 1 , 19 ) = y 6 ,
Φ ( 1 , 20 ) = y 7 ,
Φ ( 2 , 2 ) = 2 w 2 + β 4 + β 5 + β 8 τ 2 2 + η 5 τ 1 2 ,
Φ ( 2 , 3 ) = w 2 b w 3 ,
Φ ( 2 , 4 ) = w 2 d 1 y 1 ,
Φ ( 2 , 5 ) = w 2 p w 2 τ d + w 2 d 2 τ d w 4 ,
Φ ( 2 , 6 ) = w 2 d 1 w 5 y 1 ,
Φ ( 2 , 7 ) = w 2 e
Φ ( 2 , 8 ) = w 3 f ,
Φ ( 3 , 3 ) = 2 w 3 b β 9 e 2 α σ 2 + β 9 σ d ,
Φ ( 3 , 4 ) = w 3 d 1 y 5 ,
Φ ( 3 , 5 ) = w 3 p w 3 τ d + w 3 d 2 τ d + w 4 b ,
Φ ( 3 , 6 ) = w 3 d 1 + w 5 b y 5 ,
Φ ( 3 , 7 ) = w 3 e ,
Φ ( 3 , 8 ) = w 3 f ,
Φ ( 4 , 4 ) = 2 β 8 + 2 s 2 y 3 ,
Φ ( 4 , 5 ) = w 4 d 1 ,
Φ ( 4 , 6 ) = w 5 d 1 y 3 y 4 ,
Φ ( 4 , 9 ) = β 8 s
Φ ( 4 , 19 ) = y 6 ,
Φ ( 4 , 20 ) = y 7 ,
Φ ( 5 , 5 ) = 2 w 4 p 2 w 4 τ d + 2 w 4 d 2 τ d β 4 e 2 α τ 2 + β 4 τ d ,
Φ ( 5 , 6 ) = w 4 d 1 w 5 p w 5 τ d w 5 d 2 τ d ,
Φ ( 5 , 7 ) = w 4 e ,
Φ ( 5 , 8 ) = w 4 f ,
Φ ( 6 , 6 ) = 2 w 5 d 1 2 y 4 ,
Φ ( 6 , 7 ) = w 5 e ,
Φ ( 6 , 8 ) = w 5 f ,
Φ ( 6 , 19 ) = y 6 ,
Φ ( 6 , 20 ) = y 7 ,
Φ ( 7 , 7 ) = η 7 e 2 α h 2 ,
Φ ( 8 , 8 ) = η 8 e 2 α l 2 ,
Φ ( 9 , 9 ) = β 2 e 2 α τ 2 β 8 ,
Φ ( 10 , 10 ) = β 3 e 2 α σ 2 ,
Φ ( 11 , 11 ) = η 1 e 2 α τ 1 ,
Φ ( 12 , 12 ) = η 2 e 2 α σ 1 ,
Φ ( 13 , 13 ) = β 5 e 2 α σ 2 + β 5 σ d ,
Φ ( 14 , 14 ) = β 6 e 2 α τ 2 ,
Φ ( 15 , 15 ) = β 7 e 2 α σ 2 ,
Φ ( 16 , 16 ) = η 3 e 2 α τ 1 ,
Φ ( 17 , 17 ) = η 4 e 2 α σ 1 ,
Φ ( 18 , 18 ) = η 5 e 2 α τ 1 ,
Φ ( 19 , 19 ) = β 10 e 2 α σ 2 ,
Φ ( 20 , 20 ) = η 6 σ 1 e 2 α σ 1 ,
Φ ( 21 , 21 ) = η 7 e 2 α h 2 ,
Φ ( 22 , 22 ) = η 8 e 2 α l 2 ,
Φ ( 23 , 23 ) = η 9 e 2 α h 1 ,
Φ ( 24 , 24 ) = η 10 e 2 α l 1 ,
other terms are zero.
Proof. 
We study the following Lyapunov–Krasovskii functional for System (12):
V ( t , x t ) = i = 1 7 V i ( t , x t ) ,
where:
V 1 ( t , x t ) = β 1 e 2 α t x 2 ( t ) , V 2 ( t , x t ) = β 2 t τ 2 t e 2 α s x 2 ( s ) d s + β 3 t σ 2 t e 2 α s x 2 ( s ) d s + β 4 t τ ( t ) t e 2 α s x ˙ 2 ( s ) d s + β 5 t σ ( t ) t e 2 α s x ˙ 2 ( s ) d s + η 1 t τ 1 t e 2 α s x 2 ( s ) d s + η 2 t σ 1 t e 2 α s x 2 ( s ) d s , V 3 ( t , x t ) = β 6 τ 2 τ 2 0 t + θ t e 2 α s x 2 ( s ) d s d θ + β 7 σ 2 σ 2 0 t + θ t e 2 α s x 2 ( s ) d s d θ + η 3 τ 1 τ 1 0 t + θ t e 2 α s x 2 ( s ) d s d θ + η 4 σ 1 σ 1 0 t + θ t e 2 α s x 2 ( s ) d s d θ , V 4 ( t , x t ) = β 8 τ 2 τ 2 0 t + θ t e 2 α s x ˙ 2 ( s ) d s d θ + η 5 τ 1 τ 1 0 t + θ t e 2 α s x ˙ 2 ( s ) d s d θ , V 5 ( t , x t ) = β 9 t σ ( t ) t e 2 α s tanh 2 x ( s ) d s + β 10 σ 2 σ 2 0 t + θ t e 2 α s tanh 2 x ( s ) d s d θ + η 6 σ 1 σ 1 0 t + θ t e 2 α s tanh 2 x ( s ) d s d θ , V 6 ( t , x t ) = η 7 h 2 h 2 0 t + θ t e 2 α s x 2 ( s ) d s d θ + η 8 l 2 l 2 0 t + θ t e 2 α s tanh 2 x ( s ) d s d θ + η 9 h 1 h 1 0 t + θ t e 2 α s x 2 ( s ) d s d θ + η 10 l 1 l 1 0 t + θ t e 2 α s tanh 2 x ( s ) d s d θ , V 7 ( t , x t ) = γ e 2 α t x 2 ( t ) ,
where γ is a positive number that will be found later. Taking the derivative of V ( t , x t ) along any trajectory of the solution of System (12), we get:
V ˙ ( t , x t ) = i = 1 7 V i ˙ ( t , x t ) ,
where:
V ˙ 1 = 2 α β 1 e 2 α t x 2 ( t ) + 2 β 1 e 2 α t x ( t ) x ˙ ( t ) = 2 α β 1 e 2 α t x 2 ( t ) + 2 β 1 e 2 α t x ( t ) d 1 x ( t ) a x ( t ) d 1 x ( t τ ( t ) ) d 1 t τ ( t ) t x ˙ ( s ) d s + b tanh x ( t σ ( t ) ) p x ˙ ( t τ ( t ) ) + d 2 τ ˙ ( t ) x ˙ ( t τ ( t ) ) τ ˙ ( t ) x ˙ ( t τ ( t ) ) + e t h ( t ) t x ( s ) d s + f t l ( t ) t tanh x ( s ) d s + 2 w 1 e 2 α t x ( t ) d 1 x ( t ) a x ( t ) x ˙ ( t ) d 1 x ( t τ ( t ) ) d 1 t τ ( t ) t x ˙ ( s ) d s + b tanh x ( t σ ( t ) ) p x ˙ ( t τ ( t ) ) + d 2 τ ˙ ( t ) x ˙ ( t τ ( t ) ) τ ˙ ( t ) x ˙ ( t τ ( t ) ) + e t h ( t ) t x ( s ) d s + f t l ( t ) t tanh x ( s ) d s + 2 w 2 e 2 α t x ˙ ( t ) d 1 x ( t ) a x ( t ) x ˙ ( t ) d 1 x ( t τ ( t ) ) d 1 t τ ( t ) t x ˙ ( s ) d s + b tanh x ( t σ ( t ) ) p x ˙ ( t τ ( t ) ) + d 2 τ ˙ ( t ) x ˙ ( t τ ( t ) ) τ ˙ ( t ) x ˙ ( t τ ( t ) ) + e t h ( t ) t x ( s ) d s + f t l ( t ) t tanh x ( s ) d s + 2 w 3 e 2 α t tanh x ( t σ ( t ) ) d 1 x ( t ) a x ( t ) x ˙ ( t ) d 1 x ( t τ ( t ) ) d 1 t τ ( t ) t x ˙ ( s ) d s + b tanh x ( t σ ( t ) ) p x ˙ ( t τ ( t ) ) + d 2 τ ˙ ( t ) x ˙ ( t τ ( t ) ) τ ˙ ( t ) x ˙ ( t τ ( t ) ) + e t h ( t ) t x ( s ) d s + f t l ( t ) t tanh x ( s ) d s + 2 w 4 e 2 α t x ˙ ( t τ ( t ) ) d 1 x ( t ) a x ( t ) x ˙ ( t ) d 1 x ( t τ ( t ) ) d 1 t τ ( t ) t x ˙ ( s ) d s + b tanh x ( t σ ( t ) ) p x ˙ ( t τ ( t ) ) + d 2 τ ˙ ( t ) x ˙ ( t τ ( t ) ) τ ˙ ( t ) x ˙ ( t τ ( t ) ) + e t h ( t ) t x ( s ) d s + f t l ( t ) t tanh x ( s ) d s + 2 w 5 e 2 α t t τ ( t ) t x ˙ ( s ) d s d 1 x ( t ) a x ( t ) x ˙ ( t ) d 1 x ( t τ ( t ) ) d 1 t τ ( t ) t x ˙ ( s ) d s + b tanh x ( t σ ( t ) ) p x ˙ ( t τ ( t ) ) + d 2 τ ˙ ( t ) x ˙ ( t τ ( t ) ) τ ˙ ( t ) x ˙ ( t τ ( t ) ) + e t h ( t ) t x ( s ) d s + f t l ( t ) t tanh x ( s ) d s . V ˙ 2 = β 2 e 2 α t x 2 ( t ) e 2 α ( t τ 2 ) x 2 ( t τ 2 ) + β 3 e 2 α t x 2 ( t ) e 2 α ( t σ 2 ) x 2 ( t σ 2 ) + β 4 e 2 α t x ˙ 2 ( t ) e 2 α ( t τ ( t ) ) x ˙ 2 ( t τ ( t ) ) ( 1 τ ˙ ( t ) ) + β 5 e 2 α t x ˙ 2 ( t ) e 2 α ( t σ ( t ) ) x ˙ 2 ( t σ ( t ) ) ( 1 σ ˙ ( t ) ) + η 1 e 2 α t x 2 ( t ) e 2 α ( t τ 1 ) x 2 ( t τ 1 ) + η 2 e 2 α t x 2 ( t ) e 2 α ( t σ 1 ) x 2 ( t σ 1 ) .
Since τ 1 τ ( t ) τ 2 , σ 1 σ ( t ) σ 2 , τ ˙ ( t ) < τ d and σ ˙ ( t ) < σ d , we have:
V ˙ 2 e 2 α t [ β 2 x 2 ( t ) β 2 e 2 α τ 2 x 2 ( t τ 2 ) + β 3 x 2 ( t ) β 3 e 2 α σ 2 x 2 ( t σ 2 ) + β 4 x ˙ 2 ( t ) β 4 e 2 α τ 2 x 2 ( t τ ( t ) ) + τ d β 4 e 2 α τ 1 x ˙ 2 ( t τ ( t ) ) + β 5 x ˙ 2 ( t ) β 5 e 2 α σ 2 x ˙ 2 ( t σ ( t ) ) + σ d β 5 e 2 α σ 1 x ˙ 2 ( t σ ( t ) ) + η 1 x 2 ( t ) η 1 e 2 α τ 1 x 2 ( t τ 1 ) + η 2 x 2 ( t ) η 2 e 2 α σ 1 x 2 ( t σ 1 ) ] . V ˙ 3 = β 6 τ 2 τ 2 0 e 2 α t x 2 ( t ) d θ τ 2 0 e 2 α ( t + θ ) x 2 ( t + θ ) d θ + β 7 σ 2 σ 2 0 e 2 α t x 2 ( t ) d θ σ 2 0 e 2 α ( t + θ ) x 2 ( t + θ ) d θ + η 3 τ 1 τ 1 0 e 2 α t x 2 ( t ) d θ τ 1 0 e 2 α ( t + θ ) x 2 ( t + θ ) d θ + η 4 σ 1 σ 1 0 e 2 α t x 2 ( t ) d θ σ 1 0 e 2 α ( t + θ ) x 2 ( t + θ ) d θ .
By Lemma 1, we get:
V ˙ 3 e 2 α t [ β 6 τ 2 2 x 2 ( t ) + β 7 σ 2 2 x 2 ( t ) + η 3 τ 1 2 x 2 ( t ) + η 4 σ 1 2 x 2 ( t ) β 6 τ 2 e 2 α τ 2 t τ 2 t x ( s ) d s 2 β 7 σ 2 e 2 α σ 2 t σ 2 t x ( s ) d s 2 η 3 τ 1 e 2 α τ 1 t τ 1 t x ( s ) d s 2 η 4 σ 1 e 2 α σ 1 t σ 1 t x ( s ) d s 2 ] . V ˙ 4 = β 8 τ 2 τ 2 0 e 2 α t x ˙ 2 ( t ) d θ τ 2 0 e 2 α ( t + θ ) x ˙ 2 ( t + θ ) d θ + η 5 τ 1 τ 1 0 e 2 α t x ˙ 2 ( t ) d θ τ 1 0 e 2 α ( t + θ ) x ˙ 2 ( t + θ ) d θ .
By Lemma 4, we get:
V ˙ 4 e 2 α t β 8 τ 2 2 e 2 α τ 2 x ˙ 2 ( t ) + η 5 τ 1 2 x ˙ ( t ) + η 5 e 2 α τ 1 t τ 1 t x ˙ ( s ) d s 2 + e 2 α t x ( t ) x ( t τ ( t ) ) x ( t τ 2 ) β 8 β 8 s s * 2 β 8 + 2 s β 8 s * * β 8 x ( t ) x ( t τ ( t ) ) x ( t τ 2 ) . V ˙ 5 = β 9 e 2 α t tanh 2 x ( t ) β 9 e 2 α ( t σ ( t ) tanh 2 x ( t σ ( t ) ) + β 9 σ ˙ ( t ) e 2 α ( t σ ( t ) ) tanh 2 x ( t σ ( t ) ) + β 10 σ 2 σ 2 0 e 2 α t tanh 2 x ( t ) d θ β 10 σ 2 σ 2 0 e 2 α ( t + θ ) tanh 2 x ( t + θ ) d θ + η 6 σ 1 σ 1 0 e 2 α t tanh 2 x ( t ) d θ η 6 σ 1 σ 1 0 e 2 α ( t + θ ) tanh 2 x ( t + θ ) d θ .
Since tanh 2 x ( t ) x 2 ( t ) and using Lemma 1, we get:
V ˙ 5 e 2 α t [ β 9 x 2 ( t ) β 9 e 2 α σ 2 tanh 2 x ( t σ ( t ) ) + β 9 σ d e 2 α σ 1 tanh 2 x ( t σ ( t ) ) + β 10 σ 2 2 x 2 ( t ) β 10 σ 2 e 2 α σ 2 t σ 2 t tanh x ( s ) d s 2 + η 6 σ 1 2 x 2 ( t ) η 6 σ 1 e 2 α σ 1 t σ 1 t tanh x ( s ) d s 2 ] . V ˙ 6 = η 7 h 2 h 2 0 e 2 α t x 2 ( t ) e 2 α ( t + θ ) x 2 ( t + θ ) d θ + η 8 l 2 l 2 0 e 2 α t tanh 2 x ( t ) e 2 α ( t + θ ) tanh 2 x ( t + θ ) d θ + η 9 h 1 h 1 0 e 2 α t x 2 ( t ) e 2 α ( t + θ ) x 2 ( t + θ ) d θ + η 10 l 1 l 1 0 e 2 α t tanh 2 x ( t ) e 2 α ( t + θ ) tanh 2 x ( t + θ ) d θ .
Since tanh 2 x ( t ) x 2 ( t ) and using Lemma 3, we get:
V ˙ 6 e 2 α t η 7 h 2 2 x 2 ( t ) η 7 e 2 α h 2 t h ( t ) t x ( s ) d s 2 η 7 e 2 α h 2 t h 2 t h ( t ) x ( s ) d s 2 + e 2 α t [ [ η 8 l 2 2 x 2 ( t ) η 8 e 2 α l 2 t l ( t ) t tanh x ( s ) d s 2 η 8 e 2 α l 2 t l 2 t l ( t ) tanh x ( s ) d s 2 ] + e 2 α t [ η 9 h 1 2 x 2 ( t ) η 9 e 2 α h 1 t h 1 t x ( s ) d s 2 ] + e 2 α t η 10 l 1 2 x 2 ( t ) η 10 e 2 α l 1 t l 1 t tanh x ( s ) d s 2 .
From the Newton–Leibniz formula, we have:
2 y 1 x ˙ ( t ) + y 2 x ( t ) + y 3 x ( t τ ( t ) ) + y 4 t τ ( t ) t x ˙ ( s ) d s + y 5 tanh x ( t σ ( t ) ) + y 6 t σ 2 t tanh x ( s ) d s + y 7 t σ 1 t tanh x ( s ) d s x ( t ) x ( t τ ( t ) ) t τ ( t ) t x ˙ ( s ) d s .
Combining Equations (16)–(22), we get:
V ˙ * = i = 1 6 V i ( t , x t ) e 2 α t ω T ( t ) Ξ ω ( t ) ,
where ω ( t ) = x ( t ) , x ˙ ( t ) , tanh x ( t σ ( t ) ) , x ( t τ ( t ) ) , x ˙ ( t τ ( t ) ) , t τ ( t ) t x ˙ ( s ) d s , t h ( t ) t x ( s ) d s , t l ( t ) t tanh x ( s ) d s , x ( t τ 2 ) , x ( t σ 2 ) , x ( t τ 1 ) , x ( t σ 1 ) , x ˙ ( t σ ( t ) ) , t τ 2 t x ( s ) d s , t σ 2 t x ( s ) d s , t τ 1 t x ( s ) d s , t σ 1 t x ( s ) d s , t σ 1 t x ˙ ( s ) d s , t σ 2 t tanh x ( s ) d s , t σ 1 t tanh x ( s ) d s , t h 2 t h ( t ) x ( s ) d s , t l 2 t l ( t ) tanh x ( s ) d s , t h 1 t x ( s ) d s , t l 1 t tanh x ( s ) d s T ,
and Ξ is identified in Equation (13). Due to Ξ < 0 , we get V ˙ * e 2 α t ω T ( t ) Ξ ω ( t ) < 0 . Hence, there is a positive number λ such that:
V ˙ * λ e 2 α t [ | | x ˙ ( t ) | | 2 + | | x ( t ) | | 2 + | | tanh x ( t σ ( t ) ) | | 2 + | | x ( t τ ( t ) ) | | 2 + | | x ˙ ( t τ ( t ) ) | | 2 + | | t τ ( t ) t x ˙ ( s ) d s | | 2 + | | t h ( t ) t x ( s ) d s | | 2 + | | x ˙ ( t σ ( t ) ) | | 2 + | | t τ 2 t x ( s ) d s | | 2 + | | t σ 2 t x ( s ) d s | | 2 + | | t l ( t ) t tanh x ( s ) d s | | 2 + | | x ( t τ 2 ) | | 2 + | | x ( t σ 2 ) | | 2 + | | x ( t τ 1 ) | | 2 + | | x ( t σ 1 ) | | 2 + | | t τ 1 t x ( s ) d s | | 2 + | | t σ 1 t x ( s ) d s | | 2 + | | t σ 1 t x ˙ ( s ) d s | | 2 + | | t σ 2 t tanh x ( s ) d s | | 2 + | | t σ 1 t tanh x ( s ) d s | | 2 + | | t h 2 t h ( t ) x ( s ) d s | | 2 + | | t l 2 t l ( t ) tanh x ( s ) d s | | 2 + | | t h 1 t x ( s ) d s | | 2 + | | t l 1 t tanh x ( s ) d s | | 2 ] λ e 2 α t | | x ( t ) | | 2 .
Calculating the derivative of V 7 along any trajectory of Equation (4) and applying Lemma 2, then we obtain:
V ˙ 7 = 2 γ e 2 α t [ x ( t ) x ˙ ( t ) + α x 2 ( t ) ] = 2 γ e 2 α t x ( t ) a x ( t ) + b tanh x ( t σ ( t ) ) p x ˙ ( t τ ( t ) ) + e t h ( t ) t x ( s ) d s + f t l ( t ) t tanh x ( s ) d s + α x 2 ( t ) = 2 γ e 2 α t a x 2 ( t ) + b x ( t ) tanh x ( t σ ( t ) ) p x ( t ) x ˙ ( t τ ( t ) ) + e x ( t ) t h ( t ) t x ( s ) d s + f x ( t ) t l ( t ) t tanh x ( s ) d s + α x 2 ( t ) γ e 2 α t ( 2 a + 2 α ) x 2 ( t ) + b 2 tanh 2 x ( t σ ( t ) ) + p 2 x ˙ 2 ( t τ ( t ) ) + e 2 t h ( t ) t x ( s ) d s 2 + f 2 t l ( t ) t tanh x ( s ) d s 2 .
We select:
γ = λ 2 min { 1 b 2 , 1 p 2 , 1 e 2 , 1 f 2 } , if 2 a + 2 α 0 ; λ 2 min { 1 2 a + 2 α , 1 b 2 , 1 p 2 1 e 2 , 1 f 2 } , if otherwise .
We have V ˙ ( t , x t ) = i = 1 7 V i ( t , x t ) λ 2 e 2 α t | | x ( t ) | | 2 < 0 . From the conditions that 0 τ ( t ) τ 2 , 0 σ ( t ) σ 2 and V ˙ ( t , x t ) is negative definite, we get V ( t , x t ) V ( 0 , x 0 ) , t 0 , with:
V ( 0 , x 0 ) = i = 1 7 V i ( 0 , x 0 ) = β 1 x 2 ( 0 ) + β 2 τ 2 t e 2 α s x 2 ( s ) d s + β 3 σ 2 0 e 2 α s x 2 ( s ) d s + β 4 τ ( 0 ) 0 e 2 α s x ˙ 2 ( s ) d s + β 5 σ ( 0 ) t e 2 α s x ˙ 2 ( s ) d s + η 1 τ ( 0 ) 0 e 2 α s x 2 ( s ) d s + η 2 σ ( 0 ) t e 2 α s x 2 ( s ) d s + β 6 τ 2 τ 2 0 θ 0 e 2 α s x 2 ( s ) d s d θ + β 7 σ 2 σ 2 0 θ 0 e 2 α s x 2 ( s ) d s d θ + η 3 τ 1 τ 1 0 θ 0 e 2 α s x 2 ( s ) d s d θ + η 4 σ 1 σ 1 0 θ 0 e 2 α s x 2 ( s ) d s d θ + β 8 τ 2 τ 2 0 θ 0 e 2 α s x ˙ 2 ( s ) d s d θ + η 5 τ 1 τ 1 0 θ 0 e 2 α s x ˙ 2 ( s ) d s d θ + β 9 σ ( 0 ) 0 e 2 α s tanh 2 x ( s ) d s + β 10 σ 2 σ 2 0 θ 0 e 2 α s tanh 2 x ( s ) d s d θ + η 6 σ 1 σ 1 0 θ 0 e 2 α s tanh 2 x ( s ) d s d θ + η 7 h 2 h 2 0 θ 0 e 2 α s x 2 ( s ) d s d θ + η 8 l 2 l 2 0 θ 0 e 2 α s tanh 2 x ( s ) d s d θ + η 9 h 1 h 1 0 θ 0 e 2 α s x 2 ( s ) d s d θ + η 10 l 1 l 1 0 θ 0 e 2 α s tanh 2 x ( s ) d s d θ + γ x 2 ( 0 ) [ β 1 + β 2 τ 2 + β 3 σ 2 + β 4 τ 2 + β 5 σ 2 + η 1 τ 1 + η 1 σ 1 + β 6 τ 2 3 2 + β 7 σ 2 3 2 + η 3 τ 1 3 2 + η 4 σ 1 3 2 + β 8 τ 2 3 2 + η 5 σ 1 3 2 + β 9 σ 2 + β 10 σ 2 3 2 + η 6 σ 1 3 2 + η 7 h 2 3 2 + η 8 l 2 3 2 + η 9 h 1 3 2 + η 10 l 1 3 2 + γ ] max { | | x | | s 2 , | | x ˙ | | l 2 } = Δ max { | | x | | s 2 , | | x ˙ | | l 2 } = Δ N ,
where Δ = β 1 + β 2 τ 2 + β 3 σ 2 + β 4 τ 2 + β 5 σ 2 + η 1 τ 1 + η 1 σ 1 + β 6 τ 2 3 2 + β 7 σ 2 3 2 + η 3 τ 1 3 2 + η 4 σ 1 3 2 + β 8 τ 2 3 2 + η 5 σ 1 3 2 + β 9 σ 2 + β 10 σ 2 3 2 + η 6 σ 1 3 2 + η 7 h 2 3 2 + η 8 l 2 3 2 + η 9 h 1 3 2 + η 10 l 1 3 2 + γ , N = max { | | x | | s 2 , | | x ˙ | | l 2 } .
From γ e 2 α t x 2 ( t ) V ( x ( t ) ) Δ N , we obtain:
| | x ( t ) | | M e α t , M = Δ N γ .
This completes the poof. □
We show the delay-dependent criteria for the problem of the exponential stability of Equation (4) where τ 1 = σ 1 = e = f = 0 .
Corollary 1.
Consider real constants d 1 R and positive real constants τ 2 , τ d , σ 2 and σ d . Equation (4) is exponentially stable with a delay rate α > 0 , for any positive real constants β i , i = 1 , 2 , , 10 and real constants w k , m k , k = 1 , 2 , , 6 such that:
Ξ * < 0 ,
where Ξ * = [ Φ ( i , j ) * ] and [ Φ ( j , i ) * ] = [ Φ ( i , j ) * ] ,
Φ ( 1 , 1 ) * = 2 w 2 + β 4 + β 5 + β 8 τ 2 τ 2 2 w 6 ,
Φ ( 1 , 2 ) * = w 1 + w 2 d 1 w 2 a + y 1 w 6 a ,
Φ ( 1 , 3 ) * = w 2 d 1 y 1 ,
Φ ( 1 , 5 ) * = w 2 d 1 w 5 y 1 ,
Φ ( 1 , 6 ) * = w 2 b w 3 + w 6 b ,
Φ ( 1 , 7 ) * = w 2 p w 4 w 6 p ,
Φ ( 2 , 2 ) * = 2 α β 1 + 2 β 1 d 1 2 β 1 a + 2 w 1 d 1 2 w 1 a + β 2 + β 3 + β 6 τ 2 τ 2 + β 7 σ 2 σ 2 β 8 e 2 α τ 2 + β 9
+ β 10 σ 2 σ 2 + 2 y 2 ,
Φ ( 2 , 3 ) * = β 1 d 1 w 1 d 1 + β 8 e 2 α τ 2 s 1 y 2 + y 3 ,
Φ ( 2 , 4 ) * = s 1 ,
Φ ( 2 , 5 ) * = β 1 d 1 w 1 d 1 + w 5 d 1 a w 5 y 2 + y 4 ,
Φ ( 2 , 6 ) * = β 1 b + w 1 b + w 3 d 1 w 3 a + y 5 ,
Φ ( 2 , 7 ) * = β 1 p w 1 p + w 4 d 1 w 4 a ,
Φ ( 2 , 12 ) * = y 6 ,
Φ ( 3 , 3 ) * = 2 β 8 e 2 α τ 2 + 2 s 1 2 y 3 ,
Φ ( 3 , 4 ) * = β 8 e 2 α τ 2 s 1 ,
Φ ( 3 , 5 ) * = w 5 d 1 y 3 y 4 ,
Φ ( 3 , 6 ) * = w 3 d 1 y 5 ,
Φ ( 3 , 7 ) * = w 4 d 1 ,
Φ ( 3 , 12 ) * = y 6 ,
Φ ( 4 , 4 ) * = β 2 e 2 α τ 2 β 8 e 2 α τ 2 ,
Φ ( 5 , 5 ) * = 2 w 5 d 1 2 y 4 ,
Φ ( 5 , 6 ) * = w 3 d 1 + w 5 b y 5 ,
Φ ( 5 , 7 ) * = w 4 d 1 w 5 p ,
Φ ( 5 , 12 ) * = y 6 ,
Φ ( 6 , 6 ) * = 2 w 3 b β 9 e 2 α σ 2 + β 9 σ d ,
Φ ( 6 , 7 ) * = w 3 p + w 4 b ,
Φ ( 7 , 7 ) * = 2 w 4 p + τ d β 4 β 4 e 2 α τ 2 ,
Φ ( 8 , 8 ) * = β 3 e 2 α σ 2
Φ ( 9 , 9 ) * = β 5 e 2 α σ 2 + σ d β 5 ,
Φ ( 10 , 10 ) * = β 6 e 2 α τ 2 ,
Φ ( 11 , 11 ) * = β 7 e 2 α σ 2 ,
Φ ( 12 , 12 ) * = β 10 e 2 α σ 2
other terms are zero.
Proof. 
Consider the Lyapunov–Krasovskii functional candidates:
V ( t , x t ) = i = 1 6 V i ( t , x t ) ,
where:
V 1 ( t , x t ) = β 1 e 2 α t x 2 ( t ) , V 2 ( t , x t ) = β 2 t τ 2 t e 2 α s x 2 ( s ) d s + β 3 t σ 2 t e 2 α s x 2 ( s ) d s + β 4 t τ ( t ) t e 2 α s x ˙ 2 ( s ) d s + β 5 t σ ( t ) t e 2 α s x ˙ 2 ( s ) d s , V 3 ( t , x t ) = β 6 τ 2 τ 2 0 t + θ t e 2 α ( s θ ) x 2 ( s ) d s d θ + β 7 σ 2 σ 2 o t + θ t e 2 α ( s θ ) x 2 ( s ) d s d θ , V 4 ( t , x t ) = β 8 τ 2 τ 2 0 t + θ t e 2 α ( s θ ) x ˙ 2 ( s ) d s d θ , V 5 ( t , x t ) = β 9 t σ t t e 2 α s tanh 2 x ( s ) d s + β 10 σ 2 σ 2 0 t + θ t e 2 α ( s θ ) tanh 2 x ( s ) d s d θ , V 6 ( t , x t ) = γ e 2 α t x 2 ( t ) .
As stated by Theorem 1, we get the delay-dependent criteria for the exponential stability of Equation (4) where τ 1 = σ 1 = e = f = 0 . □
Remark 2.
We know that the delay-dependent condition is generally considered as less conservative than the delay-independent when the delay is small. Therefore, our criteria in Theorem 1 and Corollary 1 are delay-dependent on τ 2 , σ 2 , h 2 , and l 2 , which are less conservative than the delay-independent criteria considered in [2,6,8,27].
Remark 3.
We know that unbounded delays are more general than bounded delays, often including bounded delays as their special cases. In practical systems, there are many works such as the HIV-spread model, the modeling of oscillators, and the analysis of neural networks where infinite delays exist or modeling of one is need. Therefore, it is worth pointing out that we can extend this method to a more general system with infinite delays in future works.

4. Numerical Examples

In this part, numerical examples are presented to show the effectiveness of our main consequences, which compare the maximum allowable upper bounds σ and the parameter b in addition to investigating the rate of convergence. The feasibility of all criteria is obtained by using the LMI control toolbox in MATLAB.
Example 1.
Consider the neutral differential equation with mixed interval time-varying delays as follows:
d d t [ x ( t ) + 0.2 x ( t τ ( t ) ) ] + 0.6 x ( t ) b tanh x ( t σ ( t ) ) 0.1 t h ( t ) t x ( s ) d s 0.2 t l ( t ) t tanh x ( s ) d s = 0 .
Solving the LMI (13) when b = 0.15 , α = 0.0038 , τ ( t ) = 0.1 + cos 2 ( t ) 5 , σ ( t ) = 0.1 + sin 2 ( t ) 5 , h ( t ) = 0.1 + cos ( 2 t ) 3 and l ( t ) = 0.1 + sin ( 2 t ) 3 , we get a set of parameters to guarantee the exponential stability as follows:
η 1 = 108.2939 , η 2 = 108.2939 , η 3 = 264.6443 , η 4 = 264.6443 ,
η 5 = 268.6166 , η 6 = 505.1453 , η 7 = 290.0963 , η 8 = 427.4140 ,
η 9 = 264.6443 , η 10 = 264.6443 , β 1 = 1.5575 × 10 3 , β 2 = 104.3766 ,
β 3 = 108.3823 , β 4 = 244.4721 , β 5 = 338.3978 , β 6 = 173.6726 ,
β 7 = 173.6726 , β 8 = 248.9316 , β 9 = 336.7202 , β 10 = 173.6726 ,
w 1 = 229.6065 , w 2 = 487.9487 , w 3 = 124.9935 , w 4 = 1.5458 × 10 10 ,
w 5 = 14.5680 ,
Moreover, Table 1 lists the maximum upper bounds b for different values of σ , τ d , σ d . Table 2 lists the maximum upper bounds for the rate of convergent σ with different values of b , τ d , σ d . The maximum allowable bounds on b for different values of α , τ d , σ d can be found in Table 3. We investigate the maximum allowable bounds on the rate of convergent α in Table 4 with different values of b , τ d , σ d . In Table 5, we show the maximum allowable bounds on b for the exponential and asymptotic stabilities of Example 1 for different values of α , e . Table 6 provides the maximum allowable bounds on α for the exponential and asymptotic stabilities of Example 1 for different values of e , b . Table 7 provides the maximum allowable bounds on b for the asymptotic and exponential stabilities of Example 1 for different values of α , f . Furthermore, Table 8 presents the maximum allowable bounds on α for the asymptotic and exponential stabilities of Example 1 for different values of f , b .
Example 2.
Consider the equation focused on in [14,15] as follows:
d d t [ x ( t ) + 0.2 x ( t τ ( t ) ) ] + 0.6 x ( t ) 0.5 tanh x ( t σ ( t ) ) = 0 , t 0 ,
when σ 2 = 0.2 . and τ ( t ) = cos 2 ( t ) 10 .
Solving our criterion (23) for guaranteeing uniformly-exponential stability, α = 0.0038 is given, then we get the maximum allowable upper bound of the delay σ = 9.00 . Table 9 presents the comparison of the maximum allowable value of σ obtained in Corollary 1. We can see that our results are much better than those obtained in [14,15].
Example 3.
Consider the neutral differential equation focused on in [1,9,14,15] as follows:
d d t [ x ( t ) + 0.2 x ( t 0.1 ) ] + 0.6 x ( t ) 0.3 tanh x ( t σ ) = 0 , t 0 .
Solving our criterion (23) to guarantee uniformly-exponential stability, α = 0.0038 is given, then we get the maximum allowable upper bound of the delay σ = 175.3543 . Table 10 presents the comparison of the maximum allowable value of σ obtained in Corollary 1. We can see that our results are much better than the other works.

5. Conclusions

We improved the criteria for the exponential stability of a new certain nonlinear neutral differential equation with a more generally interval-distributed and discrete time-varying delays. By using the Lyapunov–Krasovkii functional, descriptor model transformation, the Leibniz–Newton formula, and utilization of the zero equation, the exponential stability criteria were in the form of linear matrix inequalities (LMIs). Finally, we presented numerical examples to show the effectiveness of the theoretical results and to illustrate less conservative conditions than many other works.

Author Contributions

The authors claim to have contributed significantly and equally to this work. All authors read and approved the final manuscript.

Funding

The second author was financially supported by the National Research Council of Thailand, Khon Kaen University 2019. The third author was supported by the Thailand Research Fund (TRF), the Office of the Higher Education Commission (OHEC) (Grant Number MRG6080061) and the University of Phayao, Phayao, Thailand 2019.

Acknowledgments

The authors thank the reviewers for their valuable comments and suggestions, which led to the improvement of the content of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, H. Some Improved Criteria on Exponential Stability of Neutral Differential Equation. Adv. Differ. Equ. 2012, 2012, 170. [Google Scholar] [CrossRef]
  2. El-Morshedy, H.A.; Gopalsamy, K. Nonoscillation, oscillation and convergence of a class of neutral equations. Nonlinear Anal. 2000, 40, 173–183. [Google Scholar] [CrossRef]
  3. Kwon, O.M.; Park, J.H. On improved delay-dependent stability criterion of certain neutral differential equations. Appl. Math. Comput. 2008, 199, 385–391. [Google Scholar] [CrossRef]
  4. Li, X. Global exponential stability for a class of neural networks. Appl. Math. Lett. 2009, 22, 1235–1239. [Google Scholar] [CrossRef] [Green Version]
  5. Nam, P.T.; Phat, V.N. An improved stability criterion for a class of neutral differential equations. Appl. Math. Lett. 2009, 22, 31–35. [Google Scholar] [CrossRef] [Green Version]
  6. Park, J.H. Delay-dependent criterion for asymptotic stability of a class of neutral equations. Appl. Math. Lett. 2004, 17, 1203–1206. [Google Scholar] [CrossRef]
  7. Park, J.H.; Kwon, O.M. Stability analysis of certain nonlinear differential equation. Chaos Solitons Fractals 2008, 37, 450–453. [Google Scholar] [CrossRef]
  8. Sun, Y.G.; Wang, L. Note on asymptotic stability of a class of neutral differential equations. Appl. Math. Lett. 2006, 19, 949–953. [Google Scholar] [CrossRef] [Green Version]
  9. Rojsiraphisal, T.; Niamsup, P. Exponential stability of certain neutral differential equations. Appl. Math. Comput. 2010, 217, 3875–3880. [Google Scholar] [CrossRef]
  10. Gopalsamy, K.; He, X.Z. Stability in asymmetric Hopfield nets with transmission delays. Phys. D Nonlinear Phenom. 1994, 76, 344–358. [Google Scholar] [CrossRef]
  11. Gopalsamy, K.; Leung, I. Delay induced periodicity in a neural netlet of excitation inhibition. Phys. D Nonlinear Phenom. 1996, 89, 395–426. [Google Scholar] [CrossRef]
  12. Gopalsamy, K.; Leung, I. Convergence under dynamical thresholds with delays. IEEE Trans. Neural Netw. 1997, 8, 341–348. [Google Scholar] [CrossRef] [PubMed]
  13. Gopalsamy, K.; Leung, I.; Liu, P. Global Hopf-bifurcation in a neural netlet. Appl. Math. Comput. 1998, 94, 171–192. [Google Scholar] [CrossRef]
  14. Keadnarmol, P.; Rojsiraphisal, T. Globally exponential stability of a certain neutral differential equation with time-varying delays. Adv. Differ. Equ. 2014, 2014, 32. [Google Scholar] [CrossRef] [Green Version]
  15. Chen, H.; Meng, X. An Improved Exponential Stability Criterion for a class of neutral delayed differential equations. Appl. Math. Lett. 2011, 24, 1763–1767. [Google Scholar] [CrossRef]
  16. Chatbupapan, W.; Mukdasai, K. New delay-range-dependent exponential stability criteria for certain neutral differential equations with interval discrete and distributed time-varying delay. Adv. Differ. Equ. 2016, 2016, 324. [Google Scholar] [CrossRef]
  17. Beretta, E.; Breda, D. Discrete or distributed delay? Effects on stability of population growth. Math. Biosci. Eng. 2016, 13, 19–41. [Google Scholar] [CrossRef] [PubMed]
  18. Liu, Q.; Jiang, D. Stationary distribution and extinction of a stochastic predator–prey model with distributed delay. Appl. Math. Lett. 2018, 78, 79–87. [Google Scholar] [CrossRef]
  19. Xu, C.; Shao, Y. Bifurcations in a predator-prey model with discrete and distributed time delay. Nonlinear Dyn. 2012, 63, 2207–2223. [Google Scholar] [CrossRef]
  20. Karmeshu, V.; Gupta, K.; Kadambari, V. Neuronal model with distributed delay: Analysis and simulation study for gamma distribution memory kernel. Biol. Cybern. 2011, 104, 369–383. [Google Scholar] [CrossRef]
  21. Sharma, S.K. Neuronal model with distributed delay: Emergence of unimodal and bimodal ISI distributions. IEEE Trans. Nanobiosci. 2013, 12, 1–12. [Google Scholar] [CrossRef]
  22. Maharajan, C.; Raja, R.; Cao, J.; Rajchakit, G. Novel global robust exponential stability criterion for uncertain inertial-type BAM neural networks with discrete and distributed time-varying delays via Lagrange sense. J. Frankl. Inst. 2018, 355, 4727–4754. [Google Scholar] [CrossRef]
  23. Maharajan, C.; Raja, R.; Cao, J.; Rajchakit, G.; Alsaedi, A. Novel results on passivity and exponential passivity for multiple discrete delayed neutral-type neural networks with leakage and distributed time-delays. Chaos Solitons Fractals 2018, 115, 268–282. [Google Scholar] [CrossRef]
  24. Samidurai, R.; Rajavel, S.; Zhu, Q.; Raja, R.; Zhou, H. Robust passivity analysis for neutral-type neural networks with mixed and leakage delays. Neurocomputing 2016, 175, 635–643. [Google Scholar] [CrossRef]
  25. Zhu, Q.; Senthilraj, S.; Raja, R.; Samidurai, R. Stability analysis of uncertain neutral systems with discrete and distributed delays via the delay partition approach. Int. J. Control Autom. Syst. 2017, 15, 2149–2160. [Google Scholar] [CrossRef]
  26. Khongja, N.; Botmart, T.; Niamsup, P.; Weera, W. Guaranteed cost control for exponential stability of a nonlinear system with mixed time-varying delays in states and controls. Adv. Differ. Equ. 2018, 2018, 435. [Google Scholar] [CrossRef]
  27. Agarwal, R.P.; Grace, S.R. Asymptotic Stability of Certain Neutral Differential Equations. Math. Comput. Model. Dyn. Syst. 2000, 31, 9–15. [Google Scholar] [CrossRef]
  28. El-Metwally, H.; Kulenovi, M.R.S.; Hadziomerspahi, S. Nonoscillatory solutions for system of neutral delay equation. Nonlinear Anal. 2003, 54, 63–81. [Google Scholar] [CrossRef]
  29. Kulenovi, M.R.S.; Hadziomerspahi, S. Existence of nonoscillatory solution of second order linear neutral delay equation. J. Math. Anal. Appl. 1998, 228, 436–448. [Google Scholar] [CrossRef]
  30. Kulenovic, M.R.S.; Ladas, G.; Meimaridou, A. Necessary and sufficient conditions for oscillations of neutral differential equations. J. Aust. Math. Soc. Ser. 1987, 28, 362–375. [Google Scholar] [CrossRef]
  31. Gu, K.; Kharitonov, V.L.; Chen, J. Stability of Time-Delay Systems; Birkhäuser: Berlin, Germany, 2003. [Google Scholar]
  32. Peng, C.; Fei, M.R. An Improved Result on the Stability of Uncertain T-S Fuzzy Systems with Interval Time-Varying Delay. Fuzzy Sets Syst. 2013, 212, 97–109. [Google Scholar] [CrossRef]
Table 1. Maximum allowable upper bounds b of Example 1 for α = 0.0038 .
Table 1. Maximum allowable upper bounds b of Example 1 for α = 0.0038 .
τ d = σ d σ = 0.2 σ = 0.3 σ = 0.4 σ = 0.5 σ = 0.6
0.10.42260.42240.42230.42210.4220
0.20.39840.39900.39780.39790.3970
0.30.37260.37240.37220.37200.3719
0.40.34600.34470.34450.34430.3440
0.50.31480.31450.31410.31420.3138
Table 2. Maximum allowable upper bounds σ of Example 1 for α = 0.0038 .
Table 2. Maximum allowable upper bounds σ of Example 1 for α = 0.0038 .
τ d = σ d b = 0.1 b = 0.2 b = 0.3 b = 0.4
0.05302.5884181.794490.514520.6635
0.1249.3654157.906778.036113.1826
0.15211.5763137.696566.63946.1042
0.2182.2537120.181356.15304.0409
Table 3. Maximum allowable upper bounds b of Example 1.
Table 3. Maximum allowable upper bounds b of Example 1.
τ d = σ d α = 0 α = 0.1 α = 0.2 α = 0.3 α = 0.4
0.10.42690.30360.19500.09710.0130
0.20.40240.28380.17930.08860.0108
0.30.37640.26220.16330.07920.0094
0.40.34850.24100.14560.06850.0078
0.50.31810.21290.12530.05580.0057
Table 4. Maximum allowable upper bounds α of Example 1.
Table 4. Maximum allowable upper bounds α of Example 1.
τ d = σ d b = 0.05 b = 0.1 b = 0.2 b = 0.3 b = 0.4
0.050.35630.30120.20090.11110.0296
0.10.35380.29690.19420.10310.0208
0.150.35090.29300.18800.09460.0116
0.20.34770.28670.17910.08550.0019
0.250.34390.28060.17050.07570.0110
Table 5. Maximum allowable upper bounds b of Example 1.
Table 5. Maximum allowable upper bounds b of Example 1.
e α = 0 α = 0.1 α = 0.15 α = 0.2 α = 0.25
0.050.33580.22920.18250.14010.1022
0.10.31810.21290.16700.12530.0883
0.150.30050.19650.15130.11060.0744
0.20.28280.18020.13570.10190.0605
0.250.26510.16380.12010.08200.0466
Table 6. Maximum allowable upper bounds α of Example 1.
Table 6. Maximum allowable upper bounds α of Example 1.
e b = 0.05 b = 0.1 b = 0.15 b = 0.2 b = 0.25
0.050.33170.25310.18790.13070.0792
0.10.30990.23350.16970.11350.0629
0.150.28810.21500.15150.09640.0464
0.20.26640.20100.13340.07920.0301
0.250.24470.17490.11540.06220.0138
Table 7. Maximum allowable upper bounds b of Example 1.
Table 7. Maximum allowable upper bounds b of Example 1.
f α = 0 α = 0.1 α = 0.15 α = 0.2 α = 0.25
0.050.37120.26190.21370.16970.1310
0.10.35360.24560.19810.15490.1161
0.150.33580.22920.18250.14010.1022
0.20.31810.21290.16690.12530.0883
0.250.30050.19650.15130.11060.0745
Table 8. Maximum allowable upper bounds α of Example 1.
Table 8. Maximum allowable upper bounds α of Example 1.
f b = 0.05 b = 0.1 b = 0.15 b = 0.2 b = 0.25
0.050.37480.29300.22420.16510.1130
0.10.35340.27260.20610.14790.1010
0.150.33170.25310.18790.13070.0792
0.20.30990.23350.16970.11350.0629
0.250.28810.21500.15150.09640.0464
Table 9. Maximum allowable upper bounds σ for Example 2.
Table 9. Maximum allowable upper bounds σ for Example 2.
Methods α = 0.0038 α = 0.028
H.Chen and X. Meng [15]infeasibleinfeasible
P. Keadnarmol and T. Rojsiraphisal [14]7.52310.0321
Our result (23) 9.00 0.1850
Table 10. Maximum allowable upper bounds σ of Example 3.
Table 10. Maximum allowable upper bounds σ of Example 3.
Methods α = 0.0038
H. Chen [1]175.2890
H. Chen and X. Meng [15] 10 21 (No α )
T. Rojsiraphisal and P. Niamsup [9]1.9470
P. Keadnarmol and T. Rojsiraphisal [14]175.3540
Our result (23)175.3543

Share and Cite

MDPI and ACS Style

Tranthi, J.; Botmart, T.; Weera, W.; Niamsup, P. A New Approach for Exponential Stability Criteria of New Certain Nonlinear Neutral Differential Equations with Mixed Time-Varying Delays. Mathematics 2019, 7, 737. https://doi.org/10.3390/math7080737

AMA Style

Tranthi J, Botmart T, Weera W, Niamsup P. A New Approach for Exponential Stability Criteria of New Certain Nonlinear Neutral Differential Equations with Mixed Time-Varying Delays. Mathematics. 2019; 7(8):737. https://doi.org/10.3390/math7080737

Chicago/Turabian Style

Tranthi, Janejira, Thongchai Botmart, Wajaree Weera, and Piyapong Niamsup. 2019. "A New Approach for Exponential Stability Criteria of New Certain Nonlinear Neutral Differential Equations with Mixed Time-Varying Delays" Mathematics 7, no. 8: 737. https://doi.org/10.3390/math7080737

APA Style

Tranthi, J., Botmart, T., Weera, W., & Niamsup, P. (2019). A New Approach for Exponential Stability Criteria of New Certain Nonlinear Neutral Differential Equations with Mixed Time-Varying Delays. Mathematics, 7(8), 737. https://doi.org/10.3390/math7080737

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop