Next Article in Journal
Developing a New Collision-Resistant Hashing Algorithm
Next Article in Special Issue
Improved Nonfragile Sampled-Data Event-Triggered Control for the Exponential Synchronization of Delayed Complex Dynamical Networks
Previous Article in Journal
A Parallel Convolution and Decision Fusion-Based Flower Classification Method
Previous Article in Special Issue
Impulsive Stabilization on Hyper-Chaotic Financial System under Neumann Boundary
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Stability Criteria for Delayed Neural Networks via a Relaxed Delay-Product-Type Lapunov–Krasovskii Functional

1
School of Computer, Chengdu University, Chengdu 610106, China
2
School of Electronic Information and Electrical Engineering, Chengdu University, Chengdu 610106, China
3
School of Aeronautics and Astronautics, University of Electronic Science and Technology of China, Chengdu 611731, China
*
Author to whom correspondence should be addressed.
Submission received: 14 July 2022 / Revised: 31 July 2022 / Accepted: 2 August 2022 / Published: 4 August 2022

Abstract

:
In this paper, the asymptotic stability problem of neural networks with time-varying delays is investigated. First, a new sufficient and necessary condition on a general polynomial inequality was developed. Then, a novel augmented Lyapunov–Krasovskii functional (LKF) was constructed, which efficiently introduces some new terms related to the previous information of neuron activation function. Furthermore, based on the suitable LKF and the stated negative condition of the general polynomial, two criteria with less conservatism were derived in the form of linear matrix inequalities. Finally, two numerical examples were carried out to confirm the superiority of the proposed criteria, and a larger allowable upper bound of delays was achieved.

1. Introduction

Neural networks (NNs) are generally recognized as one of the mathematical models that simulates the intelligent activities of the human brain [1]. Because of its strong ability to model complex and nonlinear relationships between the input and output, NNs have attracted widespread attention over the past decades and have been applied in many areas, such as image processing [2], optimization problems [3], signal processing [4], parallel computation [5] and other scientific fields. However, due to the limited switching frequency of amplifiers and the inherent communication time among neurons, time delays are inevitable, which may affect the performance of NNs and lead to undesired dynamics, such as oscillation or even instability [6,7,8]. Therefore, it is essential to consider an allowable upper bound of delays such that the NNs with a delay less than this bound remain stable, and many excellent results have been achieved to reduce the conservatism of the stability conditions [9,10,11,12].
It is well recognized that the Lyapunov–Krasovskii functional (LKF) method is powerful in analyzing such delayed NNs. The key point is concerned with the choice of LKF [13,14]. It has been verified that an LKF with delay-dependent matrices can yield less conservatism. A delay-product-type LKF is proposed in [15], which contains a matrix-valued polynomial of any chosen degree in time-varying delay. It shows that the conservatism is affected by the maximal degree of time-varying delay included in the delay-product-type terms of the LKF. Additionally, by taking into consideration more information on state-related vectors and the activation function, some newly augmented LKFs are presented in [13,14,16,17]. As a result, less conservative conditions are predicted, as more information about the state, delay or neuron activation functions are beneficial for reducing conservatism. Therefore, one of the research contents of this work is to build a delay-product-type LKF that contains more information on state-related vectors and the activation function, and hence is able to lead to a less conservative stability criterion.
Besides the construction of a proper LKF, the estimation of an integral term originating in the derivative of LKF is strongly related to the conservatism of criteria [18]. A common term with the form t h t θ t z T ( s ) T z ( s ) d s d θ is always included in different LKFs. Then, its derivative contains the term as follows: t h t z T ( t ) T z ( t ) d s . Various integral inequalities techniques have been proposed to bound such integral terms. The Jensen inequality [19] is known as an efficient tool to estimate the integral terms, but the reciprocally convex combination terms are usually encountered with respect to the time-varying delays, and the reciprocally convex lemma must be used to handle them [20,21]. Based on this work, many researchers have focus on the investigation of novel integral inequalities to obtain less conservative results of stability for the systems. Wirtinger-based integral inequality is investigated in [18], which encompasses Jensen inequality and obtains better results than Jensen inequality. Free-matrix-based integral inequality with a tighter bound has been proposed, while the reciprocally convex lemma is no longer required, but this increased the computational complexity [22]. Several zero-value equalities that add cross terms to the time derivative LKF terms are proposed, and improved results are derived. A novel estimating method called the general free-weighting-matrix approach is created in [17]. It is shown that this approach can incorporate Wirtinger-based inequality and free-matrix-based integral inequality. It has to be pointed out that the construction of an LKF should take integral inequalities into account to ensure that the derivative of the LKF includes those vectors appearing in the integral inequality used [23]. Therefore, utilizing a valid integral inequality technique to bound the integral terms is strongly connected with the conservatism of stability criteria. This paper investigates the stability of NNs following this direction.
In this paper, we further investigate the stability problem for NNs with time-varying delays. The contributions of this paper are summarized as follows:
(i)
A new sufficient and necessary condition for the general polynomial inequalities was developed by introducing some slack variables. The negative definiteness determination method was applied to derive less conservative stablity criteria.
(ii)
A suitable delay-product-type LKF was constructed by efficiently introducing some new terms relating to the previous information of neuron activation function. As a result, the stability criterion based on the established LKF becomes more dependent on the technique for dealing with the single integrations of activation function.
(iii)
Based on the LKF and the stated negative condition of the general polynomial, two less conservative criteria were derived, which are demonstrated through two numerical examples, and a larger allowable upper bound of delays is achieved.
Notations: Throughout this paper, N + represents the set of positive integers. R n and R m × n denote the n-dimensional vector space and the set of all m × n real matrices, respectively. Z 1 and Z T stand for the inverse and the transpose of the matrix Z, respectively. P > 0 ( P 0 ) means that P is a symmetric and positive definite (semi-definite) matrix. I and 0 denote the identity matrix and zero matrix. diag ( x 1 , , x n ) represents a diagonal matrix, in which, its diagonal elements are x 1 , , x n . col ( x 1 , , x n ) = [ x 1 T , , x n T ] T , Sym { Z } = Z + Z T . ∗ stands for the symmetric terms in a symmetric matrix. n m represents n ! m ! ( n m ) ! .

2. Preliminaries

In this section, we consider the following NNs system:
x ˙ ( t ) = A x ( t ) + W 0 f ( W x ( t ) ) + W 1 f ( W x ( t d ( t ) ) ) ,
where x ( t ) R n denotes the neuron state vector and A R n × n is a diagonal matrix. W = col W ¯ 1 , W ¯ 2 , , W ¯ n , W 0 and W 1 are the interconnection weight matrices, and f ( W x ( t ) ) = col f 1 ( W ¯ 1 x 1 ( t ) ) , f 2 ( W ¯ 2 x 2 ( t ) ) , , f n ( W ¯ n x n ( t ) R n is the activation function with f i ( 0 ) = 0 . The time-varying delay d ( t ) satisfies
0 d ( t ) h , μ d ˙ ( t ) μ ,
where h and μ are known positive constants. The neuron activation function f ( · ) is bounded by
δ i f i u 1 f i u 2 u 1 u 2 δ i + , i = 1 , 2 , , n ,
where u 1 , u 2 R , u 1 u 2 , δ i and δ i + are real scalars. For convenience, define K p = diag ( δ 1 + , δ 2 + , , δ n + ) , K m = diag ( δ 1 , δ 2 , , δ n ) .
In order to derive the stability criteria for the system (1), the following lemmas were utilized.
Lemma 1
([24]). For an integer m 0 , two scalars a and b with a < b , a real positive matrix R and a differentiable function x : [ a , b ] R n such that the integration below is well defined, then
( b a ) a b x ˙ T ( s ) R x ˙ ( s ) d s ( Θ m Λ m ϖ m ) T R m Θ m Λ m ϖ m
where R m = diag ( R , 3 R , , ( 2 m + 1 ) R ) ,
Θ m = I 0 0 I ( 1 ) 1 1 1 2 1 I 0 I ( 1 ) 1 m 1 m + 1 1 I ( 1 ) m m m 2 m m I ,
Λ m = I I 0 0 0 0 I I 0 0 0 I 0 2 I 0 0 I 0 0 m I ,
ϖ m = col ( x ( b ) , x ( a ) , γ 1 ( a , b ) , , γ m ( a , b ) ) , γ k ( a , b ) = 1 ( b a ) k a b ( b s ) k 1 x ( s ) d s , ( k = 1 , 2 , , m ) .
Lemma 2
([21]). For a given scalar α ( 0 , 1 ) , vectors ω 1 , ω 2 R n and matrices R > 0 , X 1 > 0 , X 2 > 0 , Y 1 , Y 2 R n × n , such that
R X 1 Y 1 R 0 , R Y 2 R X 2 0 ,
the following inequality holds.
1 α ω 1 T R ω 1 + 1 1 α ω 2 T R ω 2 ω 1 R + ( 1 α ) X 1 ω 1 + ω 2 R + α X 2 ω 2 + 2 ω 1 α Y 1 + ( 1 α ) Y 2 ω 2 .
Lemma 3
([25]). Consider an integrable function x ( u ) in [ a , b ] , R > 0 and a matrix T with appropriate dimensions. The following inequality is true.
a b x T ( s ) R x ( s ) d s Sym β T T λ + ( b a ) β T T R 1 T T β ,
where β is an arbitrary vector, and λ = a b x ( s ) d s .
Lemma 4.
Consider the following polynomial described by f ( z ) = a 2 n z 2 n + a 2 n 1 z 2 n 1 + + a 1 z + a 0 , where n N + , z 0 , h with h is a constant; then,
(i)
f ( z ) < 0 for all z 0 , h if and only if there exists a scalar m 0 such that
Y 2 = a 0 q 1 0 0 q 2 0 0 q 2 n 2 q 2 n 1 q 2 n < 0 ,
where q 2 n 1 = 1 2 a 2 n 1 + h m , and q 2 n = a 2 n 2 m .
(ii)
f ( z ) > 0 for all z 0 , h if and only if there exists a scalar m 0 such that
Y 1 = a 0 p 1 0 0 p 2 0 0 p 2 n 2 p 2 n 1 p 2 n > 0 ,
where p 2 n 1 = 1 2 a 2 n 1 h m , and p 2 n = a 2 n + 2 m .
Proof. 
(i) First, for z 0 , h , we define g ( z ) = z z 2 n 1 2 h z 2 n 2 + + z 3 2 h z 2 + z 2 h 0 . f ( z ) and g ( z ) can be represented as f ^ ( z ^ ) and g ^ ( z ^ ) , respectively.
f ^ ( z ^ ) = z ^ T Λ 1 z ^ , g ^ ( z ^ ) = z ^ T Λ 2 z ^ , z ^ = col ( 1 , z , , z N 1 , z N ) ,
Λ 1 = a 0 1 2 a 1 0 0 1 2 a 1 a 2 0 0 0 0 a 2 n 2 1 2 a 2 n 1 0 0 1 2 a 2 n 1 a 2 n , Λ 2 = 0 h 0 0 h 2 0 0 0 0 2 h 0 0 h 2 .
(Sufficiency) If there exists a scalar m 0 satisfying f ( z ) m g ( z ) = z ^ T Y 1 z ^ < 0 , then f ( z ) < 0 for all z [ 0 , h ] .
(Necessity) H = ( f ^ ( z ^ ) , g ^ ( z ^ ) ) : z ^ R n is a convex set. Then, we set an open convex cone as T = { ( α 1 , α 2 ) : α 1 > 0 , α 2 < 0 } . If there exists a nonzero two-tuple β 1 , β 2 , we have
α 1 β 1 α 2 β 2 > 0 , ( α 1 , α 2 ) T
β 1 f ^ ( z ^ ) β 2 g ^ ( z ^ ) 0 , ( f ^ ( z ^ ) , g ^ ( z ^ ) ) H
For any given ε > 0 , we have ε , α 2 T , for all α 1 > 0 . Thus, (11) implies that β 2 > 0 . Since g ^ ( z ^ ) < 0 is regular, (12) shows that β 1 > 0 . Hence, the following inequality can be obtained by using β 1 > 0 :
f ^ ( z ^ ) β 2 β 1 g ^ ( z ^ ) 0
for m = β 2 β 1 0 . This means that there exists m 0 such that f ^ ( z ) m g ^ ( z ) = z ^ T Y 1 z ^ < 0 . This completes the proof of (i).
(ii) Set f ˜ ( z ) = f ( z ) > 0 ; then, g ( z ) can be rewritten as g ˜ ( z ^ ) = z ^ T Λ 3 z ^ , where
Λ 3 = 0 h 0 0 h 2 0 0 0 0 2 h 0 0 h 2 ,
Then, the proof is straightforward from the proof of (i). □
Remark 1.
Lemma 4 provides necessary and sufficient conditions of high-order polynomial inequality f ( s ) < 0 and f ( s ) > 0 for z 0 , h . When a 2 n = 0 ,   f ( z ) represents an odd polynomial on z if a 2 n 1 0 . Thus, Lemma 4 covers both even and odd polynomials. It is worth mentioning that Lemma 2 in [26] is a special case of Lemma 4 with n = 1 . A quadratic function f ( z ) = a 2 z 2 + a 1 z + a 0 < 0 for all z 0 , h if and only if there exists a scalar m 0 such that
a 0 1 2 a 1 + h m 1 2 a 1 + h m a 2 2 m < 0 .

3. Main Results

In this section, two improved stability criteria are developed. For simplification, the following notations are defined.
  • v 1 ( t ) = col t d ( t ) t x ( s ) d ( t ) d s , t d ( t ) t ( t s ) x ( s ) d 2 ( t ) d s ,
  • v 2 ( t ) = col t h t d ( t ) x ( s ) h d ( t ) d s , t h t d ( t ) ( t d ( t ) s ) x ( s ) ( h d ( t ) ) 2 d s ,
  • v 3 ( t ) = t d ( t ) t f ( W x ( s ) ) d s , v 4 ( t ) = t h t d ( t ) f ( W x ( s ) ) d s ,
  • η a ( t ) = col x ( t ) , x ( t d ( t ) ) , x ( t h ) , η b ( t ) = col v 3 ( t ) , v 4 ( t ) ,
  • η 1 ( t ) = col η a ( t ) , η b ( t ) , d ( t ) v 1 ( t ) , ( h d ( t ) ) v 2 ( t ) ,
  • η 2 ( s , t ) = col x ˙ ( s ) , x ( s ) , x ( t ) , t d ( t ) s x ( θ ) d θ , s t x ( θ ) d θ ,
  • η 3 ( s , t ) = col x ( s ) , x ( s ) , x ( t ) , t h s x ( θ ) d θ , s t d ( t ) x ( θ ) d θ , η 4 ( t ) = col x ˙ ( t ) , x ( t ) , f ( W x ( t ) ) ,
  • f ( t ) = f ( W x ( t ) ) , f d ( t ) = f ( W ( t d ( t ) ) ) , f h ( t ) = f ( W x ( t h ) ) ,
  • ξ ( t ) = col η a ( t ) , x ˙ ( t d ( t ) ) , x ˙ ( t h ) , v 1 ( t ) , v 2 ( t ) , f ( t ) , f d ( t ) , f h ( t ) , η b ( t ) ,
  • e i = 0 n × ( i 1 ) n I n 0 n × ( 14 i ) n ( i = 1 , 2 , , 14 ) , e s = A e 1 + W 0 e 10 + W 1 e 11 .
Now, we present the following stability conditions.
Theorem 1.
For given scalars h > 0 and u > 0 , the system (1) is asymptotically stable if there exist symmetric matrices P i > 0 , Q i > 0 , R i > 0 ( i = 1 , 2 ) , Z > 0 , T > 0 , diagonal matrices U j > 0 , V j > 0 ( j = 1 , 2 , 3 ) , Δ k > 0 ( k = 1 , , 6 ) and any matrices X l , Y l , M l , N l ( l = 1 , 2 ) , such that
Z ˜ X 1 Y 1 Z ˜ 0 , Z ˜ Y 2 Z ˜ X 2 0 ,
M i + M i T N 1 T N 2 T h T 0 h T > 0 ,
Φ 0 ( μ ) 0.5 Φ 1 ( μ ) + h M 1 h N 2 T 0 Φ 2 ( μ ) M 1 M 1 T N 2 T N 1 T h T 0 h T < 0 ,
Φ 0 ( μ ) 0.5 Φ 1 ( μ ) + h M 2 h N 2 T 0 Φ 2 ( μ ) M 2 M 2 T N 2 T N 1 T h T 0 h T < 0 ,
where Z ˜ = diag ( Z , 3 Z , 5 Z ) , Φ ( d ( t ) , d ˙ ( t ) ) = Φ 2 ( d ˙ ( t ) ) d 2 ( t ) + Φ 1 ( d ˙ ( t ) ) d ( t ) + Φ 0 ( d ˙ ( t ) ) with
Φ 0 ( d ˙ ( t ) ) = Sym h r 11 T P 2 r 2 + d ˙ ( t ) r 11 T ( P 1 P 2 ) r 11 + r 31 T Q 1 r 31 r 41 T Q 2 r 41 ( 1 d ˙ ( t ) ) r 51 T Q 1 r 51 r 61 T Q 2 r 61 + Sym r 9 T Q 1 r 71 + r 10 T Q 2 r 81 + Sym K p W e 1 e 10 T Δ 1 + e 10 K m W e 1 T Δ 2 W e s + Sym K p W e 2 e 11 T Δ 3 + e 11 K m W e 2 T Δ 4 W e 4 ( 1 d ˙ ( t ) ) + Sym K p W e 3 e 12 T Δ 5 + e 12 K m W e 3 T Δ 6 W e 5 + Sym c = 1 3 K p W e c e c + 9 T U c e c + 9 K m W e c + Sym c = 1 2 K p W e c e c + 1 e c + 9 e c + 10 T V c e c + 9 e c + 10 K m W e c e c + 1 + Sym K p W e 1 e 3 e 10 e 12 T V 3 e 10 e 12 K m W e 1 e 3 + ρ 1 T R 1 ρ 1 + ( 1 d ˙ ( t ) ) ρ 2 T R 2 R 1 ρ 2 ρ 3 T R 2 ρ 3 + h 2 e s T Z e s ρ 4 T Z ˜ + X 1 ρ 4 ρ 5 T Z ˜ ρ 5 Sym ρ 4 T Y 2 ρ 5 + h e 10 T T e 10 + Sym e 13 T N 1 + e 14 T N 2 , Φ 1 ( d ˙ ( t ) ) = Sym r 11 T P 1 P 2 r 2 + h r 12 T P 2 r 2 + d ˙ ( t ) r 11 T ( P 1 P 2 ) r 12 + r 31 T Q 1 r 32 r 41 T Q 2 r 42 ( 1 d ˙ ( t ) ) r 51 T Q 1 r 52 r 61 T Q 2 r 62 + r 9 T Q 1 r 72 + r 10 T Q 2 r 82 1 h Sym ρ 4 T Y 1 Y 2 ρ 5 + 1 h ρ 4 T X 1 ρ 4 ρ 5 T X 2 ρ 5 , Φ 2 ( d ˙ ( t ) ) = Sym r 12 T ( P 1 P 2 ) r 2 + d ˙ ( t ) r 12 T ( P 1 P 2 ) r 12 + r 32 T Q 1 r 32 r 42 T Q 2 r 42 ( 1 d ˙ ( t ) ) r 52 T Q 1 r 52 r 62 T Q 2 r 62 + Sym r 9 T Q 1 r 73 + r 10 T Q 2 r 83 ,
r 1 = r 11 + d ( t ) r 12 , r 11 = col ( e 1 , e 2 , e 3 , e 13 , e 14 , 0 , 0 , h e 8 , h e 9 ) , r 12 = col ( 0 , 0 , 0 , 0 , 0 , e 6 , e 2 , e 8 , e 9 ) , r 2 = col ( e s , ( 1 d ˙ ( t ) ) e 4 , e 5 , e 10 ( 1 d ˙ ( t ) ) e 11 , ( 1 d ˙ ( t ) ) e 11 e 12 , e 1 ( 1 d ˙ ( t ) ) e 2 , ( 1 d ˙ ( t ) ) e 2 + e 6 d ˙ ( t ) e 7 , ( 1 d ˙ ( t ) ) e 2 e 3 , e 3 + ( 1 d ˙ ( t ) ) e 8 + d ˙ ( t ) e 9 ) , r i = r i 1 + d ( t ) r i 2 ( i = 3 , 4 , 5 , 6 ) , r j = r j 1 + d ( t ) r j 2 + d 2 ( t ) r j 3 ( j = 7 , 8 ) , r 31 = col ( e s , e 1 , e 1 , 0 , 0 ) , r 32 = col ( 0 , 0 , 0 , e 6 , 0 ) , r 41 = col ( e 5 , e 3 , e 1 , 0 , h e 8 ) , r 42 = col ( 0 , 0 , 0 , 0 , e 8 ) , r 51 = col ( e 4 , e 2 , e 1 , 0 , 0 ) , r 52 = col ( 0 , 0 , 0 , 0 , e 6 ) , r 61 = col ( e 4 , e 2 , e 1 , h e 8 , 0 ) , r 62 = col ( 0 , 0 , 0 , e 8 , 0 ) , r 71 = col ( e 1 e 2 , 0 , 0 , 0 , 0 ) , r 72 = col ( 0 , e 6 , e 1 , 0 , 0 ) , r 73 = col ( 0 , 0 , 0 , e 7 , e 6 e 7 ) , r 83 = col ( 0 , 0 , 0 , e 9 , e 8 e 9 ) , r 81 = col ( e 2 e 3 , h e 8 , h e 1 , h 2 e 9 , h 2 ( e 8 e 9 ) , r 82 = col ( 0 , e 8 , e 1 , 2 h e 9 , 2 h ( e 8 e 9 ) ) , r 9 = col ( 0 , 0 , e s , ( 1 d ˙ ( t ) ) e 2 , e 1 ) , r 10 = col ( 0 , 0 , e s , e 3 , ( 1 d ˙ ( t ) ) e 2 ) , ρ 1 = col ( e s , e 1 , e 10 ) , ρ 2 = col ( e 4 , e 2 , e 11 ) , ρ 3 = col ( e 5 , e 3 , e 12 ) , ρ 4 = col ( e 1 e 2 , e 1 + e 2 2 e 6 , e 1 e 2 6 e 6 + 12 e 7 ) , ρ 5 = col ( e 2 e 3 , e 2 + e 3 2 e 8 , e 2 e 3 6 e 8 + 12 e 9 ) .
Proof. 
To prove the stability of the system (1), we constructed the following LKF:
V ( t ) = i = 1 6 V i ( t ) ,
where
V 1 ( t ) = η 1 T ( t ) P ( d ( t ) ) η 1 ( t ) , V 2 ( t ) = t d ( t ) t η 2 T ( s , t ) Q 1 η 2 ( s , t ) d s + t h t d ( t ) η 3 T ( s , t ) Q 2 η 3 ( s , t ) d s , V 3 ( t ) = 2 i = 1 n 0 W i x ( t ) λ 1 i δ i + s f i ( s ) + λ 2 i f i ( s ) δ i s d s , + 2 i = 1 n 0 W i x ( t d ( t ) ) λ 3 i δ i + s f i ( s ) + λ 4 i f i ( s ) δ i s d s , + 2 i = 1 n 0 W i x ( t h ) λ 5 i δ i + s f i ( s ) + λ 6 i f i ( s ) δ i s d s , V 4 ( t ) = t d ( t ) t η 4 T ( s ) R 1 η 4 ( s ) d s + t h t d ( t ) η 4 T ( s ) R 2 η 4 ( s ) d s , V 5 ( t ) = h t h t ( h t + s ) x ˙ T ( s ) Z x ˙ ( s ) d s , V 6 ( t ) = t h t θ t f T ( W x ( s ) ) T f ( W x ( s ) ) d s d θ .
where P ( d ( t ) ) = d ( t ) P 1 + ( h d ( t ) ) P 2 .
Then, taking the derivatives of V ( t ) along the trajectory of system (1) yields
V ˙ 1 ( t ) = 2 η 1 T ( t ) P ( d ( t ) ) η ˙ 1 ( t ) + η 1 T ( t ) P ˙ ( d ( t ) ) η 1 ( t ) , V ˙ 2 ( t ) = η 2 T ( t , t ) Q 1 η 2 ( t , t ) ( 1 d ˙ ( t ) ) η 2 T ( t d ( t ) , t ) Q 1 η 2 ( t d ( t ) , t ) η 3 T ( t h , t ) Q 2 η 3 ( t h , t ) + ( 1 d ˙ ( t ) ) η 3 T ( t d ( t ) , t ) Q 2 η 3 ( t d ( t ) , t )
+ 2 t η 2 T ( s , t ) Q 1 t d ( t ) t η 2 ( s , t ) d s + 2 t η 3 T ( s , t ) Q 2 t h t d ( t ) η 3 ( s , t ) d s , V ˙ 3 ( t ) = Sym K p W x ( t ) f ( t ) T Δ 1 + f ( t ) K m W x ( t ) T Δ 2 W x ˙ ( t ) + ( K p W x ( t d ( t ) ) f d ( t ) ) T Δ 3 + f d ( t ) K m W x ( t d ( t ) ) T Δ 4 W x ˙ ( t d ( t ) ) ( 1 d ˙ ( t ) )
+ K p W x ( t h ) f h ( t ) T Δ 5 + f h ( t ) K m W x ( t h ) T Δ 6 W x ˙ ( t h ) , V ˙ 4 ( t ) = η 4 T ( t ) R 1 η 4 ( t ) + ( 1 d ˙ ( t ) ) η 4 T ( t d ( t ) ) R 2 R 1 η 4 ( t d ( t ) )
η 4 T ( t h ) R 2 η 4 ( t h ) ,
V ˙ 5 ( t ) = h 2 x ˙ T ( t ) Z x ˙ ( t ) h t h t x ˙ T ( s ) Z x ˙ ( s ) d s ,
V ˙ 6 ( t ) = h f T ( W x ( t ) ) T f ( W x ( t ) ) t h t f T ( W x ( t ) ) T f ( W x ( t ) ) d s .
Based on Lemmas 1 and 2, it follows from (14) that
h t h t x ˙ T ( s ) Z x ˙ ( s ) d s = h t h t d ( t ) x ˙ T ( s ) Z x ˙ ( s ) d s h t d ( t ) t x ˙ T ( s ) Z x ˙ ( s ) d s h d ( t ) ξ T ( t ) ρ 4 T Z ρ 4 ξ ( t ) h h d ( t ) ξ T ( t ) ρ 5 T Z ρ 5 ξ ( t ) ξ T ( t ) ρ 4 T Z ˜ + X 1 ρ 4 ρ 5 T Z ˜ ρ 5 Sym ρ 4 T Y 2 ρ 5 + d ( t ) h ρ 4 T X 1 ρ 4 ρ 5 T X 2 ρ 5 Sym ρ 4 T Y 1 Y 2 ρ 5 ξ ( t ) ,
Using Lemma 3, one has the following inequality:
t h t f T ( W x ( t ) ) T f ( W x ( t ) ) d s ξ T ( t ) ( Sym e 13 T N 1 + e 14 T N 2 + d ( t ) N 1 T T 1 N 1 + ( h d ( t ) ) N 2 T T 1 N 2 ) ξ ( t ) .
To make use of the information of the activation function, we can obtain
ϕ j ( u ) = 2 K p W x ( u ) f ( W x ( u ) ) T U j f ( W x ( u ) ) K m W x ( u ) 0 , φ j u 1 , u 2 = 2 K p W ( x ( u 1 ) x ( u 2 ) ) ( f ( W x ( u 1 ) ) f ( W x ( u 2 ) ) ) T
× V j ( f ( W x ( u 1 ) ) f ( W x ( u 2 ) ) ) K m W ( x ( u 1 ) x ( u 2 ) ) 0 ,
where u 1 , u 2 R , U j = diag u j 1 , u j 2 , , u j n S + n , V j = diag v j 1 , v j 2 , , v j n S + n ( j = 1 , 2 , 3 ) ; then, it is clear that
ϕ 1 ( t ) + ϕ 2 ( t d ( t ) ) + ϕ 3 t h 0 ,
φ 1 ( t , t d ( t ) ) + φ 2 t d ( t ) , t h + φ 3 t , t h 0 .
From (19) to (30), the following inequality can be obtained:
V ˙ ( t ) ξ T ( t ) Ω ξ ( t ) ,
where
Ω = Φ 2 ( d ˙ ( t ) ) d 2 ( t ) + Φ 1 ( d ˙ ( t ) ) + N 1 T T 1 N 1 N 2 T T 1 N 2 d ( t ) + Φ 0 ( d ˙ ( t ) ) + h N 2 T T 1 N 2 .
To guarantee that V ˙ ( t ) < 0 , applying Lemma 2 and M ¯ i + M ¯ i T > 0 ( i = 1 , 2 ) , Ω < 0 is equivalent to
Φ 0 ( d ˙ ( t ) ) 0.5 Φ 1 ( d ˙ ( t ) ) + h M i Φ 2 ( d ˙ ( t ) ) M i M i T + h N 1 T T 1 N 1 N 1 T T 1 N 1 1 h N 1 T T 1 N 1 + 1 h N 2 T 1 N 2 T < 0 ,
where M ¯ i = 1 2 h N 1 T T 1 N 1 1 2 h N 2 T 1 N 2 T + M i ( i = 1 , 2 ) . Finally, if (16) and (17) are satisfied, then V ˙ ( t ) < 0 for M i satisfying (15) by a Schur complement [27], which means that the NNs (1) are asymptotically stable. The proof is completed. □
Remark 2.
To derive the stability criterion in terms of linear matrix inequality, the matrix M ¯ i ( i = 1 , 2 ) should be appropriately designed. In Theorem 1, M ¯ i ( i = 1 , 2 ) are set as M ¯ i = 1 2 h N 1 T T 1 N 1 1 2 h N 2 T 1 N 2 T + M i ( i = 1 , 2 ) to handle 1 2 N 1 T T 1 N 1 and 1 2 N 2 T 1 N 2 T of Ω. The slack variable M ¯ i ( i = 1 , 2 ) helps in reducing the conservativeness by giving the freedom for checking the feasibility of the linear matrix inequality conditions.
Remark 3.
It should be pointed out that most literature [11,28] only propose the information of activation functions, such as f ( x ( t ) ) , f ( x ( t d ( t ) ) ) and f ( x ( t d ) ) . In this paper, the information of single integrations of activation functions, i.e., t d ( t ) t f ( W x ( s ) ) d s and t h t d ( t ) f ( W x ( s ) ) d s , are considered in V 1 ( t ) . According to this construction, more activation function information on NNs can be adequately utilized, which leads to less conservative results. t h t θ t f T ( W x ( s ) ) T f ( W x ( s ) ) d s d θ is firstly proposed in [11]; different from its derivative handled method, this paper estimates the integral term by Lemma 3, which makes full use of cross-terms from the single integrations of activation function to others, and leads to a tighter bound.
Corollary 1.
For given scalars h > 0 and u > 0 , the system (1) is asymptotically stable if there exist symmetric matrices P > 0 , Z > 0 , T > 0 , Q i > 0 , R i > 0 ( i = 1 , 2 ) , diagonal matrices U j > 0 , V j > 0 ( j = 1 , 2 , 3 ) , Δ k > 0 ( k = 1 , , 6 ) and any matrices X l , Y l , M l , N l ( l = 1 , 2 ) , such that
Z ˜ X 1 Y 1 Z ˜ 0 , Z ˜ Y 2 Z ˜ X 2 0 ,
M i + M i T N 1 T N 2 T h T 0 h T > 0 ,
Φ 0 ( μ ) 0.5 Φ 1 ( μ ) + h M 1 h N 2 T 0 Φ 2 ( μ ) M 1 M 1 T N 2 T N 1 T h T 0 h T < 0 ,
Φ 0 ( μ ) 0.5 Φ 1 ( μ ) + h M 2 h N 2 T 0 Φ 2 ( μ ) M 2 M 2 T N 2 T N 1 T h T 0 h T < 0 ,
where
Φ 0 ( d ˙ ( t ) ) = Sym r 11 T P r 2 + r 31 T Q 1 r 31 r 41 T Q 2 r 41 ( 1 d ˙ ( t ) ) r 51 T Q 1 r 51 r 61 T Q 2 r 61 + Sym r 9 T Q 1 r 71 + r 10 T Q 2 r 81 + Sym K p W e 1 e 10 T Δ 1 + e 10 K m W e 1 T Δ 2 W e s + Sym K p W e 2 e 11 T Δ 3 + e 11 K m W e 2 T Δ 4 W e 4 ( 1 d ˙ ( t ) ) + Sym K p W e 3 e 12 T Δ 5 + e 12 K m W e 3 T Δ 6 W e 5 + Sym c = 1 3 K p W e c e c + 9 T U c e c + 9 K m W e c + Sym c = 1 2 K p W e c e c + 1 e c + 9 e c + 10 T V c e c + 9 e c + 10 K m W e c e c + 1 + Sym K p W e 1 e 3 e 10 e 12 T V 3 e 10 e 12 K m W e 1 e 3 + ρ 1 T R 1 ρ 1 + ( 1 d ˙ ( t ) ) ρ 2 T R 2 R 1 ρ 2 ρ 3 T R 2 ρ 3 + h 2 e s T Z e s ρ 4 T Z ˜ + X 1 ρ 4 ρ 5 T Z ˜ ρ 5 Sym ρ 4 T Y 2 ρ 5 + h e 10 T T e 10 + Sym e 13 T N 1 + e 14 T N 2 Φ 1 ( d ˙ ( t ) ) = Sym r 12 T p r 2 + r 31 T Q 1 r 32 r 41 T Q 2 r 42 ( 1 d ˙ ( t ) ) r 51 T Q 1 r 52 r 61 T Q 2 r 62 + r 9 T Q 1 r 72 + r 10 T Q 2 r 82 1 h Sym ρ 4 T Y 1 Y 2 ρ 5 + 1 h ρ 4 T X 1 ρ 4 ρ 5 T X 2 ρ 5 , Φ 2 ( d ˙ ( t ) ) = r 32 T Q 1 r 32 r 42 T Q 2 r 42 ( 1 d ˙ ( t ) ) r 52 T Q 1 r 52 r 62 T Q 2 r 62 + Sym r 9 T Q 1 r 73 + r 10 T Q 2 r 83 ,
where the other notions are defined well in Theorem 1.
Proof. 
From the proof of Theorem 1, we choose the LKF (19) with P ( d ( t ) ) = P in V 1 ( t ) . The rest of the process is similar to that in the proof of Theorem 1; hence, we omit it here. □
Remark 4.
Theorem 1 and Corollary 1 are related to the stability problem of NNs with time-varying delays. It is possible to apply the proposed methods to other control problems, e.g., H control [29] and sampled-data control [30,31].

4. Illustrative Examples

In this section, two numerical examples are presented to demonstrate the superiority and the effectivity of the derived criteria.
Example 1.
Consider the delayed NNs (1) with the following parameters:
A = diag ( 1.5 , 0.7 ) , W = diag ( 1 , 1 ) , K m = diag ( 0 , 0 ) , K p = diag ( 0.3 , 0.8 ) , W 0 = 0.0503 0.0454 0.0987 0.2075 , W 1 = 0.2381 0.9320 0.0388 0.5062 .
For various μ , the allowable maximum delay bounds by Theorem 1 and Corollary 1 in this paper and the related methods in [17,32,33,34] are listed in Table 1. The results show that Theorem 1 provides a bigger allowable maximum upper bound of delays than the related methods of [17,32,33,34] do, which verifies the lower conservatism of the proposed criteria. Theorem 1 and Corollary 1 are derived via a similar approach; it can be seen that Theorem 1 provides larger allowable upper delay bounds than Corollary 1, which means that the LKF that contains delay-dependent matrices can lead to a less conservative result. In addition, when h = 42.2617 , μ = 0.4 , f ( x ) = col ( 0.3 tanh ( x 1 ) , 0.8 tanh ( x 2 ) ) , the initial state x ( t ) = col ( 1 , 1 ) , t [ 42.2617 , 0 ] and d ( t ) = 42.2617 2 + 42.2617 2 sin ( 0.8 42.2617 t ) . The state trajectories of the x ( t ) given in Figure 1 show that the NNs (1) are stable at their equilibrium points, which verifies the effectiveness of our methods.
Example 2.
Consider the delayed NNs (1) with the following parameters:
A = diag ( 7.3458 , 6.9987 , 5.5949 ) , K m = diag ( 0 , 0 , 0 ) , K p = diag ( 0.3680 , 0.1795 , 0.2876 ) , W 0 = diag ( 0 , 0 , 0 ) , W 1 = diag ( 1 , 1 , 1 ) , W = 13.6014 2.9616 0.6936 7.4736 21.6810 3.2100 0.7290 2.6334 20.1300 .
It can be easily found that the proposed stability criteria can produce a larger allowable maximum upper bound for all cases than those given in the existing literature. This shows that the proposed criteria are indeed less conservative than the ones in the literature. Table 2 shows the allowable maximum upper bound of delays for various μ . When h = 2.4705 , μ = 0.5 , d ( t ) = 2.4705 2 + 2.4705 2 sin ( 1 2.4705 t ) , the initial state x ( t ) = col ( 0.2 , 0.1 , 0.3 ) , t [ 2.4705 , 0 ] and the activation function f ( x ) = col ( 0.3680 tanh ( x 1 ) , 0.1795 tanh ( x 2 ) , 0.2876 tanh ( x 3 ) ) ; the state trajectories of NNs in Example 2 are shown in Figure 2. The results given reveal that the LKF with a more integral term of activation function (i.e., t d ( t ) t f ( W x ( s ) ) d s and t h t d ( t ) f ( W x ( s ) ) d s ) has improved the results.

5. Conclusions

This paper has investigated the stability problem of NNs with time-varying delays. A positivity and negativity definiteness condition for a general polynomial with any chosen degree in time-varying delay has been established and a relaxed LKF that includes more previous information of the neuron activation function has been proposed. Benefiting from the relaxed LKF and the negative definiteness determination method, two stability criteria with less conservatism have been derived, which has been demonstrated by two numerical examples. Nevertheless, the definiteness condition is still proposed with conservatism, and future work on reducing its conservatism is meaningful.

Author Contributions

Conceptualization, S.W. and K.S.; formal analysis, S.W. and K.S.; funding acquisition, K.S.; methodology, S.W., K.S. and J.Y.; project administration, K.S.; software, S.W.; supervision, K.S.; validation, S.W. and J.Y.; visualization, S.W. and J.Y.; writing—original draft, S.W.; writing—review and editing, K.S. and J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China under Grant (No. 61703060 and 12061088), the Sichuan Science and Technology Program under Grant No. 21YYJC0469, the Project funded by China Postdoctoral Science Foundation under Grant Nos. 2020M683274 and 2021T140092, Central Government Funds of Guiding Local Scientific and Technological Development for Sichuan Province of China under grant 2021ZYD0012, the Opening Fund of Geomathematics Key Laboratory of Sichuan Province (scsxdz2020zd01).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, T.H.; Park, M.J.; Park, J.H. An improved stability criterion of neural networks with time-varying delays in the form of quadratic function using novel geometry-based conditions. Appl. Math. Comput. 2021, 404, 126–226. [Google Scholar] [CrossRef]
  2. Chua, L.; Yang, L. Cellular neural networks: Applications. IEEE Trans. Circuits Syst. 1988, 35, 1273–1290. [Google Scholar] [CrossRef]
  3. Mestari, M.; Benzirar, M.; Saber, N.; Khouil, M. Solving Nonlinear Equality Constrained Multiobjective Optimization Problems Using Neural Networks. IEEE Trans. Neural Netw. Learn. Syst. 2015, 26, 2500–2520. [Google Scholar] [CrossRef]
  4. Zhang, H.; Wang, Z.; Liu, D. A Comprehensive Review of Stability Analysis of Continuous-Time Recurrent Neural Networks. IEEE Trans. Neural Netw. Learn. Syst. 2014, 25, 1229–1262. [Google Scholar] [CrossRef]
  5. Chen, Y.H.; Fang, S.C. Neurocomputing with time delay analysis for solving convex quadratic programming problems. IEEE Trans. Neural Netw. Learn. Syst. 2000, 11, 230–240. [Google Scholar] [CrossRef] [PubMed]
  6. Melin, P.; Castillo, O. An intelligent hybrid approach for industrial quality control combining neural networks, fuzzy logic and fractal theory. Inf. Sci. 2007, 177, 1543–1557. [Google Scholar] [CrossRef]
  7. Melin, P.; Mendoza, O.; Castillo, O. Face Recognition With an Improved Interval Type-2 Fuzzy Logic Sugeno Integral and Modular Neural Networks. IEEE Trans. Syst. Man Cybern. A Syst. Hum. 2011, 41, 1001–1012. [Google Scholar] [CrossRef]
  8. Melin, P.; Sánchez, D.; Castillo, O. Genetic optimization of modular neural networks with fuzzy response integration for human recognition. Inf. Sci. 2012, 197, 1–19. [Google Scholar] [CrossRef]
  9. Zhang, X.M.; Han, Q.L.; Ge, X.; Ding, D. An overview of recent developments in Lyapunov–Krasovskii functionals and stability criteria for recurrent neural networks with time-varying delays. Neurocomputing 2018, 313, 392–401. [Google Scholar] [CrossRef]
  10. Zhang, X.M.; Han, Q.L.; Ge, X. An overview of neuronal state estimation of neural networks with time-varying delays. Inf. Sci. 2019, 478, 83–99. [Google Scholar] [CrossRef]
  11. Chen, J.; Park, J.H.; Xu, S. Stability Analysis for Delayed Neural Networks With an Improved General Free-Matrix-Based Integral Inequality. IEEE Trans. Neural Netw. Learn. Syst. 2020, 31, 675–684. [Google Scholar] [CrossRef]
  12. Xu, B.; Li, B. Event-Triggered State Estimation for Fractional-Order Neural Networks. Mathematics 2022, 10, 325. [Google Scholar] [CrossRef]
  13. Kwon, O.M.; Park, M.J.; Lee, S.M.; Park, J.H.; Cha, E.J. Stability for Neural Networks With Time-Varying Delays via Some New Approaches. IEEE Trans. Neural Netw. Learn. Syst. 2013, 24, 181–193. [Google Scholar] [CrossRef]
  14. Zhang, X.M.; Han, Q.L. Global Asymptotic Stability for a Class of Generalized Neural Networks with Interval Time-Varying Delays. IEEE Trans. Neural Netw. Learn. Syst. 2011, 22, 1180–1192. [Google Scholar] [CrossRef]
  15. Long, F.; Zhang, C.K.; He, Y.; Wang, Q.G.; Gao, Z.M.; Wu, M. Hierarchical Passivity Criterion for Delayed Neural Networks via A General Delay-Product-Type Lyapunov-Krasovskii Functional. IEEE Trans. Neural Netw. Learn. Syst. 2021, 1–12. [Google Scholar] [CrossRef]
  16. Kwon, O.M.; Park, J.H.; Lee, S.M.; Cha, E.J. New augmented Lyapunov–Krasovskii functional approach to stability analysis of neural networks with time-varying delays. Nonlinear Dyn. 2014, 76, 221–236. [Google Scholar] [CrossRef]
  17. Zhang, C.K.; He, Y.; Jiang, L.; Lin, W.J.; Wu, M. Delay-dependent stability analysis of neural networks with time-varying delay: A generalized free-weighting-matrix approach. Appl. Math. Comput. 2017, 294, 102–120. [Google Scholar] [CrossRef]
  18. Seuret, A.; Gouaisbaut, F. Wirtinger-based integral inequality: Application to time-delay systems. Automatica 2013, 49, 2860–2866. [Google Scholar] [CrossRef] [Green Version]
  19. Gu, K. An integral inequality in the stability problem of time-delay systems. In Proceedings of the Proceedings of the 39th IEEE Conference on Decision and Control (Cat. No.00CH37187), Sydney, Australia, 12–15 December 2000; Volume 3, pp. 2805–2810. [Google Scholar]
  20. Zhang, X.M.; Lin, W.J.; Han, Q.L.; He, Y.; Wu, M. Global Asymptotic Stability for Delayed Neural Networks Using an Integral Inequality Based on Nonorthogonal Polynomials. IEEE Trans. Neural Netw. Learn. Syst. 2018, 29, 4487–4493. [Google Scholar] [CrossRef]
  21. Park, P.; Ko, J.W.; Jeong, C. Reciprocally convex approach to stability of systems with time-varying delays. Automatica 2011, 47, 235–238. [Google Scholar] [CrossRef]
  22. Zeng, H.B.; He, Y.; Wu, M.; She, J. Free-Matrix-Based Integral Inequality for Stability Analysis of Systems with Time-Varying Delay. IEEE Trans. Autom. Control 2015, 60, 2768–2772. [Google Scholar] [CrossRef]
  23. Zhang, X.M.; Han, Q.L.; Seuret, A.; Gouaisbaut, F. An improved reciprocally convex inequality and an augmented Lyapunov–Krasovskii functional for stability of linear systems with time-varying delay. Automatica 2017, 84, 221–226. [Google Scholar] [CrossRef] [Green Version]
  24. Zhang, X.M.; Han, Q.L.; Ge, X. Novel stability criteria for linear time-delay systems using Lyapunov-Krasovskii functionals with a cubic polynomial on time-varying delay. IEEE/CAA J. Autom. Sin. 2021, 8, 77–85. [Google Scholar] [CrossRef]
  25. Park, M.; Kwon, O.; Ryu, J. Passivity and stability analysis of neural networks with time-varying delays via extended free-weighting matrices integral inequality. Neural Netw. 2018, 106, 67–78. [Google Scholar] [CrossRef] [PubMed]
  26. Park, J.; Park, P. Finite-interval quadratic polynomial inequalities and their application to time-delay systems. J. Frankl. Inst. 2020, 357, 4316–4327. [Google Scholar] [CrossRef]
  27. Ouellette, D.V. Schur complements and statistics. Linear Algebra Appl. 1981, 36, 187–295. [Google Scholar] [CrossRef] [Green Version]
  28. Zeng, H.B.; Park, J.H.; Zhang, C.F.; Wang, W. Stability and dissipativity analysis of static neural networks with interval time-varying delay. J. Frankl. Inst. 2015, 352, 1284–1295. [Google Scholar] [CrossRef]
  29. Chang, X.H.; Yang, G.H. Nonfragile H Filtering of Continuous-Time Fuzzy Systems. IEEE Trans. Signal Process. 2011, 59, 1528–1538. [Google Scholar] [CrossRef]
  30. Lee, T.H.; Wu, Z.G.; Park, J.H. Synchronization of a complex dynamical network with coupling time-varying delays via sampled-data control. Appl. Math. Comput. 2012, 219, 1354–1366. [Google Scholar] [CrossRef]
  31. Zeng, H.B.; Teo, K.L.; He, Y.; Wang, W. Sampled-data stabilization of chaotic systems based on a T-S fuzzy model. Inf. Sci. 2019, 483, 262–272. [Google Scholar] [CrossRef]
  32. Chen, J.; Park, J.H.; Xu, S. Stability Analysis for Neural Networks with Time-Varying Delay via Improved Techniques. IEEE Trans. Cybern. 2019, 49, 4495–4500. [Google Scholar] [CrossRef]
  33. Zhang, X.M.; Han, Q.L.; Zeng, Z. Hierarchical Type Stability Criteria for Delayed Neural Networks via Canonical Bessel–Legendre Inequalities. IEEE Trans. Cybern. 2018, 48, 1660–1671. [Google Scholar] [CrossRef]
  34. Chen, Y.; Chen, G. Stability analysis of delayed neural networks based on a relaxed delay-product-type Lyapunov functional. Neurocomputing 2021, 439, 340–347. [Google Scholar] [CrossRef]
  35. Li, Z.; Yan, H.; Zhang, H.; Zhan, X.; Huang, C. Stability Analysis for Delayed Neural Networks via Improved Auxiliary Polynomial-Based Functions. IEEE Trans. Neural Netw. Learn. Syst. 2019, 30, 2562–2568. [Google Scholar] [CrossRef] [PubMed]
Figure 1. State trajectories for μ = 0.4 and h = 42.2617 .
Figure 1. State trajectories for μ = 0.4 and h = 42.2617 .
Mathematics 10 02768 g001
Figure 2. State trajectories for μ = 0.5 and h = 2.4705 .
Figure 2. State trajectories for μ = 0.5 and h = 2.4705 .
Mathematics 10 02768 g002
Table 1. The allowable upper bound of delays for different μ (Example 1).
Table 1. The allowable upper bound of delays for different μ (Example 1).
Methods μ = 0.4 μ = 0.5
Theorem 1 (Case II) [32]8.56697.2809
Theorem 3 [17]13.867110.0050
Proposition 3 ( N = 2 )  [33]23.840914.8593
Corollary 1 [34]24.232215.0214
Theorem 1 [34]33.582818.4200
Corollary 127.459414.9889
Theorem 142.261722.4641
Improvement25.8433%21.9549%
Table 2. The allowable upper bound of delays for different μ (Example 2).
Table 2. The allowable upper bound of delays for different μ (Example 2).
Methods μ = 0.1 μ = 0.5
Theorem 1 ( N = 2 ) [11]1.15450.6010
Proposition 3 ( N = 3 ) [33]1.15540.5967
Corollary 1 [34]1.16120.6162
Theorem 1 [34]1.18450.6572
Theorem 1 [35]2.04970.9860
Theorem 1 [1]2.36791.0187
Corollary 17.91451.5412
Theorem 110.93342.4705
Improvement361.7340%142.5149%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, S.; Shi, K.; Yang, J. Improved Stability Criteria for Delayed Neural Networks via a Relaxed Delay-Product-Type Lapunov–Krasovskii Functional. Mathematics 2022, 10, 2768. https://0-doi-org.brum.beds.ac.uk/10.3390/math10152768

AMA Style

Wang S, Shi K, Yang J. Improved Stability Criteria for Delayed Neural Networks via a Relaxed Delay-Product-Type Lapunov–Krasovskii Functional. Mathematics. 2022; 10(15):2768. https://0-doi-org.brum.beds.ac.uk/10.3390/math10152768

Chicago/Turabian Style

Wang, Shuoting, Kaibo Shi, and Jin Yang. 2022. "Improved Stability Criteria for Delayed Neural Networks via a Relaxed Delay-Product-Type Lapunov–Krasovskii Functional" Mathematics 10, no. 15: 2768. https://0-doi-org.brum.beds.ac.uk/10.3390/math10152768

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