Next Article in Journal
Biomechanical Effects of Medializing Calcaneal Osteotomy on Bones and the Tissues Related to Adult-Acquired Flatfoot Deformity: A Computational Study
Next Article in Special Issue
Asymptotic Behavior of Solutions in Nonlinear Neutral System with Two Volterra Terms
Previous Article in Journal
Optimal Private Health Insurance Contract towards the Joint Interests of a Policyholder and an Insurer
Previous Article in Special Issue
Exponential Stability of Nonlinear Time-Varying Delay Differential Equations via Lyapunov–Razumikhin Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Criteria on Exponential Incremental Stability of Dynamical Systems with Time Delay

1
School of Mathematical Sciences, Fudan University, Shanghai 200433, China
2
Shanghai Center for Mathematical Sciences, Fudan University, Shanghai 200433, China
3
Shanghai Key Laboratory for Contemporary Applied Mathematics, Shanghai 200433, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(10), 2242; https://doi.org/10.3390/math11102242
Submission received: 14 April 2023 / Revised: 7 May 2023 / Accepted: 8 May 2023 / Published: 10 May 2023

Abstract

:
Incremental stability analysis for time-delay systems has attracted more and more attention for its contemporary applications in transportation processes, population dynamics, economics, satellite positions, etc. This paper researches the criteria for exponential incremental stability for time-delay systems with continuous or discontinuous right-hand sides. Firstly, the sufficient conditions for exponential incremental stability for time-delay systems with continuous right-hand sides are studied, and several corollaries for specific cases are provided. As for time-delay systems with discontinuous right-hand sides, after expounding the relevant conditions for the existence and uniqueness of the Filippov solution, by using approximation methods, sufficient conditions for exponential incremental stability are obtained. The conclusions are applied to linear switched time-delay systems and Hopfield neural network systems with composite right-hand sides.

1. Introduction

Stability analysis in mathematics mainly refers to the relevant research on the long-term performance of the dynamical system’s steady state. With more and more applications of neural networks and complex systems, the stability analysis of differential equations, including time-delay dynamical systems, has attracted more and more attention from academia and industry. Incremental stability [1,2] has been presented to be a perfect instrument for stability analysis, which is able to address problems of synchronization of coupled systems.
Incremental stability means that as time approaches infinity, the solutions of the dynamical system in different initial states will approach each other, that is, the state variables of the system with different initial states will gradually converge to the same trajectory. This property has a very wide range of applications in different fields of academia and industry. In recent years, due to the increasing potential application value in many frontier fields, such as PI-controlled missiles, Ref. [3] as well as the synchronization problem of network dynamics [4,5,6], there is already a lot of literature available on incremental stability, e.g., Refs. [7,8] provided a systematic exposition and discussion of related issues, and Ref. [3] provided specific examples of incremental stability-related applications.
In dynamical system analysis, the theoretical research on differential equations with discontinuous right-hand sides has also been highly valued because of its wide application. In some fields, such as mechanical engineering, electronic engineering, and automatic control theory, many problems rely on relevant theories of these ‘discontinuous’ differential equations [9]. Among them, switched systems, as a type of differential equation with discontinuous right-hand functions, are particularly commonly used in the field of automatic control, thus driving the development of related theories [10,11,12,13].
In 1964, Filippov [14] studied the motion of Coulomb friction oscillators and proposed a differential equation with a discontinuous right-hand side. In order to study the trajectory of the solution, ‘differential inclusion’ and set-valued mapping were introduced, and the existence and uniqueness of the solution of the discontinuous differential equation were discussed. A detailed discussion on this type of discontinuous differential equation can be found in reference [15]. Before we study the compressibility of the system, we need to first ensure the existence and uniqueness of the Carathéodory solution of the system. Filippov’s theory mainly focuses on the existence and uniqueness of the solution of the non-smooth dynamical system. The relevant conclusions have been listed in [15].
Research on time-delay differential equations first began in the early 20th century by Volterra [16,17]. The ordinary differential equations with time delay are commonly used in contemporary applications, gradually pushing the relevant theories of time-delay systems to integrity and maturity and producing rich results. There are many important achievements emerging one after another. It is worth mentioning Hale and Verduyn Lunel’s comprehensive work [18], which discussed in detail the properties of the solutions for some time-delay differential equations, e.g., uniqueness, continuous dependence of parameters, continuity and compactness of solutions, stability and invariance, etc. At the same time, some concrete analyses and methods on the properties of these solutions for time-delay systems were proposed in the literature [19,20,21,22,23], among which, Ref. [20] includes an introductory chapter that provides detailed examples of time-delay differential equations used to control computer systems, transportation processes, population dynamics, economics, satellite positions, urban transportation, and so on.
Set-valued dynamical systems, also named Filippov systems, whose right-hand sides are set-valued mappings, are widely used in these applications mentioned above. These set-valued dynamical systems are perfect instruments to represent the time-delay differential equations with discontinuous right-hand sides [20,24,25,26,27,28] or control systems with time delay [20,29,30]. Therefore, naturally, a large amount of the literature has discussed the problems along these lines. One of the most important achievements in this field is [31], in which Haddad focused on upper semicontinuous dynamics, elaborated on the existence and compactness of the solution set, and also proved the upper semicontinuity of the solution. Haddad’s work [31] is given under functional differential inclusion, where the corresponding time-delay term acts on the infinite-dimensional space of continuous functions.
In recent years, there have been many related analytical studies and achievements on the incremental stability of Filippov systems. In the case that the local Lipschisz condition is satisfied, Ref. [32] provided a sufficient condition for the local stability of Filippov solutions. Ref. [33] used the concept of Filippov solutions to analyze a class of time-delay dynamical systems with discontinuous right-hand sides. In the sense of the Filippov solution, Ref. [34] proposed the conditions for global asymptotic stability of the error system of the time-delay neural network with a discontinuous activation function. While [35] put forward an approximation method and gave the specific sufficient conditions for the exponential incremental stability of the switched system.
In this paper, first, several preliminary definitions are given in Section 2. In Section 3, we research the criteria on the exponential incremental stability of the solutions for time-delay systems with continuous right-hand sides, involving several specific cases, and relevant corollaries are provided. Then, in Section 4, under the hypothesis that the system has a unique solution, we extend the sufficient conditions for exponential incremental stability to the time-delay system with a discontinuous right-hand side in the sense of the Filippov solution by using a sequence of ‘continuous systems’ to approach the corresponding Filippov system. In this section, we also provide the conditions for the existence and uniqueness of the solution for the time-delay dynamical system before stability analysis in Section 4.1. The applications on linear switched systems and Hopfield neural network systems with time delay are given in Section 5, respectively, and corresponding numerical examples are given in Section 6.

2. Preliminaries

Here, we first introduce some primary definitions, including matrix measure and multiple norms, incremental stability property, exponential incremental stability, and so on.
Definition 1 
(Definition 1 in [35]). For any real matrix A R n × n and a given norm · , we define the corresponding matrix measure ν ( A ) as
ν ( A , · ) = lim h 0 + I + h A 1 h .
The matrix measure above can be considered as the one-sided directional derivative of the induced matrix norm function · , evaluated at the point I, in the direction of A.
In the following parts, we will study the incremental stability property of time-delay systems under multiple norms. Here, we also list the definitions of multiple norms with subscript χ ( t ) and the corresponding measures. Note that the function χ ( t ) is a piecewise right-continuous function.
Definition 2. 
For real matrix A R n × n and the matrix norm · χ ( t ) , (the corresponding vector measure | · | χ ( t ) ), we here define the corresponding matrix measure μ χ ( t ) ( A ) as follows:
μ χ ( t ) ( A ) = lim h 0 + I + h A χ ( t ) 1 h .
If lim h 0 + | x | χ ( t + h ) ( lim h 0 | x | χ ( t + h ) ) exists, then we denote the right (left) limit of the norm | · | χ ( t ) at time point t by | · | χ ( t ± ) . We say the norm | · | χ ( t ) is continuous at t, if and only if | · | χ ( t + ) = | · | χ ( t ) = | · | χ ( t ) , that is, | · | χ ( t ) is right-continuous and left-continuous as well. If there exists D > 0 , such that | x | χ ( t ) < D | x | χ ( s ) holds for all t , s R + , we say | · | χ ( t ) is uniformly equivalent.
Then, we extend the definition for matrix measure in the sense of multiple norms. Consider the time-varying nature of χ ( t ) , as follows:
Definition 3 
(Definition 3 in [36]). If the following limit exists, the switched matrix measure with respect to vector norm | · | χ ( t ) is defined as follows:
ν χ ( t ) ( A ) = lim h 0 + ¯ 1 h sup | x | χ ( t ) = 1 | ( I n + h A ) x | χ ( t + h ) 1
where lim ¯ stands for the upper superior.
Remark 1. 
The definition of multiple matrix norms · χ ( t ) can be thought of as matrix norms induced by vector norm | · | χ ( t ) : for a matrix A,
A χ ( t ) = sup | x | χ ( t ) = 1 | A x | χ ( t )
This implies that, if χ ( t ) is constant over an internal [ a , a + δ ) , then it holds that ν χ ( t ) ( A ) = μ χ ( t ) ( A ) in [ a , a + δ ) .
According to the existence of the switched matrix measure, the definition is as follows,
Definition 4 
([36]). Define the partial differential of the switched norm | · | χ ( t ) as follows,
¯ t ( | · | χ ( t ) ) = lim h 0 + ¯ sup | x | χ ( t ) = 1 | x | χ ( t + h ) 1 h
If ¯ t ( | · | χ ( t ) ) exists at t, we say the multiple norm | · | χ ( t ) is right regular at time t.
According to the ‘right regular’ property, we have the following proposition [36]:
Proposition 1 
([36]). If the multiple norm | · | χ ( t ) is right regular, then
1. 
The multiple norm | · | χ ( t ) is right-continuous at time t;
2. 
ν χ ( t ) ( · ) exists at time t.
For a clearer statement in the following part, we define a transaction function between norms | · | χ and | · | χ ,
Definition 5. 
Function C ( χ , χ ) > 0 is the transaction function between norms | · | χ and | · | χ , satisfying that
| · | χ C ( χ , χ ) | · | χ .
Definition 6. 
If the function C ( χ ( t ) , χ ( t ) ) is well-defined for all t and t , we say the multiple norm | · | χ ( t ) is equivalent for all t.
Here, we consider the following dynamical time-delay system:
x ˙ = f ( x , x τ 1 ( t ) , , x τ m ( t ) , r ( t ) ) , t t 0 x ( s ) = ϕ ( s ) , s [ t 0 τ ¯ , t 0 ]
where x R n , x τ k ( t ) = x ( t τ k ( t ) ) ( k = 1 , , m ) is the time-delay term, each τ k is a bounded function, max k s u p t [ t 0 , ) τ k ( t ) = τ ¯ , function r ( t ) : [ t 0 , + ) R is an upper continuous staircase function, f = ( f 1 , , f n ) : R n × ( m + 1 ) × [ t 0 , + ) R n . The function ϕ ( · ) represents the initial value function, and ϕ ( · ) C 1 ( [ t 0 τ ¯ , t 0 ] , R n ) . Let x ( t ; ϕ , r t ) be the solution of system (1).
Then, we will research the sufficient conditions for incremental stability for time-delay dynamical system (1), which is defined as follows.
We have the following definition for several types of incremental stability (IS):
Definition 7. 
Let ϕ ( t ; t 0 , x 0 ) be the solution of system (1) with the initial time t 0 and initial value function x 0 ( · ) . If there exist a function β ( s , t ) of class KL , and some norm | · | [ t 0 τ ¯ , t 0 ] induced by vector norm | · | , defined as | φ | [ t 0 τ ¯ , t 0 ] = sup s [ t 0 τ ¯ , t 0 ] | φ ( s ) | , such that for any initial function x 0 ( · ) , y 0 ( · ) ,
| ϕ ( t + t 0 ; x 0 , t 0 ) ϕ ( t + t 0 ; y 0 , t 0 ) | β ( | x 0 y 0 | [ t 0 τ ¯ , t 0 ] , t ) ,
then, we say that system (1) is Incrementally Asymptotically Stable (δAS) in the region Σ R n . If β ( s , t ) is independent of initial time t, then we say system (1) is Incrementally Uniformly Asymptotically Stable (δUAS). If function β ( s , t ) is of class EKL , then we say system (1) is Incrementally Uniformly Exponentially Asymptotically Stable (δUEAS).
Moreover, if there exists a constant M and c > 0 such that β ( · , · ) of class KL satisfies (2) with some norm | · | , then system (1) is said to be exponentially incrementally stable.
β ( | x 0 y 0 | [ t 0 τ ¯ , t 0 ] , t ) = M e c t | x 0 y 0 | [ t 0 τ ¯ , t 0 ] ,

3. Contraction Theory for Time-Delay Systems

Here, we list the following hypothesis, denoted by Assumption 1, including the Carathéodory condition to guarantee the existence and uniqueness of the solution of the time-delay dynamical system (1) with a continuous right-hand side.
Assumption 1. 
Dynamical time-delay system (1) satisfies the following conditions:
1. 
f ( x , x τ 1 ( t ) , , x τ m ( t ) , r ( t ) ) is continuously differentiable with respect to x, and continuous with respect to ( x , x τ 1 ( t ) , , x τ m ( t ) , r ) except for the switching time points { t 1 , , t j , } .
2. 
τ k ( t ) is upper bounded and has a positive lower bound for each k, and max k sup t [ t 0 , ) τ k ( t ) = τ ¯ , min k inf t [ t 0 , ) τ k ( t ) = τ ̲ .
3. 
f i : R n × ( m + 1 ) × [ t 0 , + ) R ( i = 1 , , n ) is locally Lipschitz.
Thus, under Assumption 1, system (1) has a unique solution. (Refer to the Reference [37] for details.)
Here, we try to research the contraction property for time-delay system (1). First, enlightened by Ref. [38], we prove the following lemma:
Lemma 1. 
Assume that | · | χ ( t ) is right regular, χ ( t ) is a right-continuous staircase function with discontinuities in { t 1 , , t j , } and | · | χ ( t ) C 0 | · | χ ( t ) holds for any t , t t 0 . x ( t ) is continuous, t 0 is the initial time, and τ ¯ = max k sup s t 0 τ k ( s ) . For the following time-delay system:
x ˙ ( t ) = A ( t ) x ( t ) + k = 1 n B k ( t ) x ( t τ k ( t ) )
where x ( t ) R n , A ( t ) , B k ( t ) R n × n ( k = 1 , 2 , , m ) is piecewise continuous with respect to t, and the discontinuities belong to { t 1 , t 2 , , t i , } , which is a countable set. If there exists a piecewise right-continuous function m ( t ) > 0 whose discontinuities belong to { t 1 , t 2 , , t i , } , and matrix-valued functions B k ( 1 ) ( t ) , B k ( 2 ) ( t ) such that
ν χ ( t ) ( A ( t ) + k = 1 m B k ( 1 ) ( t ) ) + k = 1 m B k ( 2 ) ( t ) χ ( t ) m ( t ) m ( t τ k ( t ) ) + k = 1 m τ k ( t ) B k ( 1 ) ( t ) χ ( t ) ( A ˜ k ( t ) + B ˜ k ( t ) ) D + m ( t ) m ( t )
where D + represents the Dini derivative, and
A ˜ k ( t ) = sup t τ k ( t ) s t A ( s ) χ ( s ) m ( t ) m ( s ) , B ˜ k ( t ) = sup t τ k ( t ) s t B k ( s ) χ ( s ) m ( t ) m ( s τ k ( s ) ) .
Then, for each t [ t j , t j + 1 ) , we have
| x ( t ) | χ ( t j ) m ( t ) sup t 0 τ ¯ s t | x ( s ) | χ ( t j ) m ( s ) C 0 sup t 0 τ ¯ s t j | x ( s ) | χ ( t j ) m ( s )
Proof. 
Let V ( t ) = sup t 0 τ ¯ θ t | x ( θ ) | χ ( θ ) m ( θ ) . Assume that V ( t ) is strictly increasing at time point t * , it implies that V ( t * ) = | x ( t * ) | χ ( t * ) m ( t * ) . There exists j such that t * [ t j , t j + 1 ) , and here we calculate the Dini derivative of | x ( t ) | χ ( t j ) m ( t ) ,
D + [ | x ( t ) | χ ( t j ) m ( t ) ] = m ˙ ( t ) | x ( t ) | χ ( t j ) + m ( t ) lim h 0 | x ( t + h ) | χ ( t j ) | x ( t ) | χ ( t j ) h ,
in which x ˙ ( t ) can be rewritten as
x ˙ ( t ) = A ( t ) x ( t ) + k = 1 m B k ( t ) x ( t τ k ( t ) ) = [ A ( t ) + k = 1 m B k ( 1 ) ( t ) ] x ( t ) + k = 1 m B k ( 1 ) ( t ) t t τ k ( t ) x ˙ ( s ) d s + k = 1 m B k ( 2 ) ( t ) x ( t τ k ( t ) )
Thus, we have
lim ¯ h 0 | x ( t + h ) | χ ( t j ) | x ( t ) | χ ( t j ) h = lim ¯ h 0 1 h [ | x ( t ) + h ( A ( t ) + k = 1 m B k ( 1 ) ( t ) ) x ( t ) + h k = 1 m B k ( 1 ) ( t ) t t τ k ( t ) x ˙ ( s ) d s + h k = 1 m B k ( 2 ) ( t ) x ( t τ k ( t ) ) | χ ( t j ) | x ( t ) | χ ( t j ) ] = lim ¯ h 0 1 h | x ( t ) | χ ( t j ) I + h ( A ( t ) + k = 1 m B k ( 1 ) ( t ) ) χ ( t j ) 1 + k = 1 m B k ( 1 ) ( t ) χ ( t j ) t τ k ( t ) t | x ˙ ( s ) | χ ( t j ) d s + k = 1 m B k ( 2 ) ( t ) χ ( t j ) | x ( t τ k ( t ) ) | χ ( t j ) = ν χ ( t j ) ( A ( t ) + k = 1 m B k ( 1 ) ( t ) ) | x ( t ) | χ ( t j ) + k = 1 m B k ( 1 ) ( t ) χ ( t j ) t τ k ( t ) t | x ˙ ( s ) | χ ( t j ) d s + k = 1 m B k ( 2 ) ( t ) χ ( t j ) | x ( t τ ( t ) ) | χ ( t j )
that is,
D + [ | x ( t ) | χ ( t j ) m ( t ) ] = m ˙ ( t ) m ( t ) | x ( t ) | χ ( t j ) m ( t ) + ν χ ( t j ) ( A ( t ) + k = 1 m B k ( 1 ) ( t ) ) | x ( t ) | χ ( t j ) m ( t ) + k = 1 m B k ( 2 ) ( t ) χ ( t j ) | x ( t τ k ( t ) ) | χ ( t j ) + k = 1 m B k ( 1 ) ( t ) χ ( t j ) m ( t ) t τ k ( t ) t | A ( s ) x ( s ) + k = 1 m B k ( s ) x ( s τ k ( s ) ) | χ ( t j ) d s = m ˙ ( t ) m ( t ) | x ( t ) | χ ( t j ) m ( t ) + ν χ ( t j ) ( A ( t ) + k = 1 m B k ( 1 ) ( t ) ) | x ( t ) | χ ( t j ) m ( t ) + k = 1 m B k ( 2 ) ( t ) χ ( t j ) | x ( t τ k ( t ) ) | χ ( t j ) m ( t τ k ( t ) ) m ( t ) m ( t τ k ( t ) ) + k = 1 m B k ( 1 ) ( t ) χ ( t j ) t τ k ( t ) t A ( s ) χ ( t j ) | x ( s ) | χ ( t j ) m ( s ) m ( t ) m ( s ) + k = 1 m B k ( s ) χ ( t j ) | x ( s τ k ( s ) ) | χ ( t j ) m ( s τ k ( s ) ) m ( t ) m ( s τ k ( s ) ) d s
Let
A ˜ k ( t ) = sup t τ k ( t ) s t A ( s ) χ ( t j ) m ( t ) m ( s ) B ˜ k ( t ) = sup t τ k ( t ) s t B k ( s ) χ ( t j ) m ( t ) m ( s τ k ( s ) )
then, for t = t * , we have
D + V ( t ) [ m ˙ ( t ) m ( t ) + ν χ ( t j ) ( A ( t ) + k = 1 m B k ( 1 ) ( t ) ) + k = 1 m τ k ( t ) B k ( 1 ) ( t ) χ ( t j ) ( A ˜ k ( t ) + B ˜ k ( t ) ) + k = 1 m B k ( 2 ) ( t ) χ ( t j ) m ( t ) m ( t τ k ( t ) ) ] V ( t ) 0
which is contradictory to the hypothesis that V ( t ) is strictly increasing at t * . Therefore, V ( t ) is decreasing for t [ t j , t j + 1 ) , that is,
V ( t ) V ( t j )
that is,
| x ( t ) | χ ( t j ) m ( t ) sup t 0 τ ¯ s t j | x ( s ) | χ ( t j ) m ( s ) C 0 sup t 0 τ ¯ s t j | x ( s ) | χ ( t j ) m ( s )
The lemma is proved. □
With the conclusion in Lemma 1, we can select a proper function m ( t ) > 0 , and prove the following theorem for exponential incremental stability property of system (1).
Theorem 1. 
Suppose that Assumption 1 holds for time-delay system (1). χ ( t ) is a right-continuous staircase function with discontinuities in { t 1 , , t j , } . Let N ( t ) = # { j : t t j , j = 1 , 2 , } , and there exist that
1. 
A constant T 0 > 0 ,
2. 
Positive constants γ k ( k = 0 , 1 , ) , γ 0 γ 1 ,
3. 
C ( χ ( t ) , χ ( t ) ) C 0 for any t , t t 0 ,
such that the following conditions hold:
Let A ( t ) = ( f / x ) ( x , x τ 1 , , x τ m , r ( t ) ) , B k ( t ) = ( f / x τ k ) ( x , x τ 1 , , x τ m , r ( t ) ) , k = 1 , , m be piecewise continuous with respect to t, and the discontinuities belong to { t 1 , t 2 , , t i , } . There exist matrix-valued functions B k ( 1 ) ( t ) , B k ( 2 ) ( t ) , B k ( 1 ) ( t ) + B k ( 2 ) ( t ) = B k ( t ) such that
ν χ ( t ) ( A ( t ) + k = 1 m B k ( 1 ) ( t ) ) + k = 1 m τ k ( t ) B k ( 1 ) ( t ) χ ( t ) ( A ˜ k ( t ) + B ˜ k ( t ) ) + k = 1 m B k ( 2 ) ( t ) χ ( t ) exp ( γ N ( t ) γ N ( t τ k ( t ) ) ) t exp γ N ( t τ k ( t ) ) τ k ( t ) γ N ( t )
where
A ˜ t k = sup t τ k ( t ) s t A ( s ) χ ( t ) exp γ N ( t ) t γ N ( s ) s B ˜ t k = sup t τ k ( t ) s t B k ( s ) χ ( t ) exp γ N ( t ) t γ N ( s τ k ( s ) ) ( s τ k ( s ) ) .
and for any t > T 0 ,
1 t i = 0 N ( t ) 1 [ ( γ i + 1 γ i ) t i + 1 + log β i + 1 ] + γ 0 t 0 γ N ( t ) t c
where β k = C ( χ ( t k ) , χ ( t k ) ) , then system (1) is exponentially incrementally stable.
Proof. 
For any initial state x 0 , y 0 , denote the corresponding initial function by x 0 ( · ) , y 0 ( · ) , x 0 ( · ) , y 0 ( · ) C ( [ t 0 τ ¯ , t 0 ] , R m ) . Here, we define a function φ λ ( · ) as φ λ ( s ) = ( 1 λ ) x 0 ( s ) + λ y 0 ( s ) , λ [ 0 , 1 ] , s [ t 0 τ ¯ , t 0 ] , which is the initial value function of the initial state φ λ . Let φ 0 = x 0 , φ 1 = y 0 . f ( x , x τ 1 , , x τ m , r ( t ) ) is continuous with respect to ( x , x τ 1 , , x τ m , r ) except for the switching time points { t j } , and is continuously differentiable with respect to x, so we have the solution ψ ( t , λ ) = x ( t ; φ λ , r t ) , which is continuously differentiable with respect to φ λ . Let ω = ψ / λ , thus we can conclude that ω is well-defined and continuous. From chain rule, ω ( t , λ ) is the solution of the following system:
ω ˙ = 2 ψ λ t = f λ = f x ω + k = 1 m f x τ k ω τ k ω ( s , λ ) = φ λ λ ( s ) , s [ t 0 τ ¯ , t 0 ]
ψ / λ is well-defined and continuous, with condition (3), from the proof of Lemma 1, select a piecewise right-continuous function m ( t ) = exp ( γ N ( t ) t ) > 0 whose discontinuities belong to { t 1 , t 2 , , t i , } , we have
m ( t ) | ω ( t , λ ) | χ ( t ) V ( t ) V ( t j ) = sup t 0 τ ¯ s t j m ( s ) | ω ( s , λ ) | χ ( t j )
holds for any t [ t j , t j + 1 ) . Thus, under the condition (4), we have
m ( t ) | ω ( t , λ ) | χ ( t ) V ( t ) V ( t j ) = sup t 0 τ ¯ s t j m ( s ) | ω ( s , λ ) | χ ( t j ) exp ( γ j γ j 1 ) t j · β j sup t 0 τ ¯ s t j m ( s ) | ω ( s , λ ) | χ ( t j ) exp ( γ j γ j 1 ) t j · β j sup t 0 τ ¯ s t j 1 m ( s ) | ω ( s , λ ) | χ ( t j 1 ) j = 1 N ( t ) exp ( γ j γ j 1 ) t j · β j · sup t 0 τ ¯ s t 0 m ( s ) | ω ( s , λ ) | χ ( t 0 ) j = 1 N ( t ) exp ( γ j γ j 1 ) t j · β j · sup θ [ t 0 τ ¯ , t 0 ] m ( θ ) sup t 0 τ ¯ s t 0 | ω ( s , λ ) | χ ( t 0 )
then, we have
| ω ( t , λ ) | χ ( t ) exp ( j = 1 N ( t ) ( γ j γ j 1 ) t j + log β j ) exp γ 0 t 0 γ N ( t ) t sup t 0 τ ¯ s t 0 | ω ( t 0 , λ ) | χ ( t 0 ) e c t sup t 0 τ ¯ s t 0 | ω ( t 0 , λ ) | χ ( t 0 ) .
Therefore, together with the condition (4), we conclude that the time-delay system (1) is incrementally uniformly asymptotically stable:
| x ( t ; y 0 , r t ) x ( t ; x 0 , r t ) | χ ( t 0 ) C 0 | x ( t ; y 0 , r t ) x ( t ; x 0 , r t ) | χ ( t ) C 0 | 0 1 ψ ( t , λ ) λ d λ | χ ( t ) C 0 0 1 | ω ( t , λ ) | χ ( t ) d λ C 0 e c t 0 1 sup t 0 τ ¯ s t 0 | ω ( t 0 , λ ) | χ ( t 0 ) d λ C 0 e c t sup θ [ t 0 τ ¯ , t 0 ] | x 0 ( θ ) y 0 ( θ ) | χ ( t 0 )
Actually, the key thought of Theorem 1 is replacing the time-delay term B k ( t ) x τ k ( t ) with B k ( 1 ) ( t ) x τ k ( t ) + B k ( 2 ) ( t ) x ( t ) + B k ( 2 ) ( t ) t t τ k ( t ) x ˙ ( s ) d s . In some special cases, we can set B k ( 1 ) ( t ) = 0 or B k ( 2 ) ( t ) = 0 for k = 1 , 2 , , m , which infers the following corollaries.
Corollary 1. 
Suppose that Assumption 1 holds. χ ( t ) is a right-continuous staircase function with discontinuities in { t 1 , , t j , } . Let N ( t ) = # { j : t t j , j = 1 , 2 , } , and there exist
1. 
A constant T 0 > 0 ,
2. 
Positive constants γ k ( k = 0 , 1 , ) , γ 0 γ 1 ,
3. 
C ( χ ( t ) , χ ( t ) ) C 0 for any t , t t 0 ,
such that the following conditions hold:
Let A ( t ) = ( f / x ) ( x , x τ 1 , , x τ m , r ( t ) ) , B k ( t ) = ( f / x τ k ) ( x , x τ 1 , , x τ m , r ( t ) ) , k = 1 , , m be piecewise continuous with respect to t, and the discontinuities belong to { t 1 , t 2 , , t i , } ,
ν χ ( t ) ( A ( t ) ) + k = 1 m B k ( t ) χ ( t ) exp γ N ( t ) t exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) γ N ( t )
and for any t > T 0 ,
1 t i = 0 N ( t ) 1 [ ( γ i + 1 γ i ) t i + 1 + log β i + 1 ] + γ 0 t 0 γ N ( t ) t c ,
where β k = C ( χ ( t k ) , χ ( t k ) ) , then system (1) is exponentially incrementally stable.
Corollary 2. 
Suppose that Assumption 1 holds. χ ( t ) is a right-continuous staircase function with discontinuities in { t 1 , , t j , } . Let N ( t ) = # { j : t t j , j = 1 , 2 , } , and there exist
1. 
A constant T 0 > 0 ,
2. 
Positive constants α k > 0 ( k = 0 , 1 , ) ,
3. 
C ( χ ( t ) , χ ( t ) ) C 0 for any t , t t 0 ,
such that the following conditions hold:
Let A ( t ) = ( f / x ) ( x , x τ 1 , , x τ m , r ( t ) ) , B k ( t ) = ( f / x τ k ) ( x , x τ 1 , , x τ m , r ( t ) ) , k = 1 , , m , there exist matrix-valued functions B k ( 1 ) ( t ) , B k ( 2 ) ( t ) , B k ( 1 ) ( t ) + B k ( 2 ) ( t ) = B k ( t ) such that
ν χ ( t ) ( A ( t ) + k = 1 m B k ( t ) ) + k = 1 m τ k ( t ) B k ( t ) χ ( t ) ( A ˜ k ( t ) + B ˜ k ( t ) ) γ N ( t )
where
A ˜ t k = sup t τ k ( t ) s t A ( s ) χ ( t ) exp γ N ( t ) t γ N ( s ) s B ˜ t k = sup t τ k ( t ) s t B k ( s ) χ ( t ) exp γ N ( t ) t γ N ( s τ k ( s ) ) ( s τ k ( s ) )
and for any t > T 0 ,
1 t i = 0 N ( t ) 1 [ ( γ i + 1 γ i ) t i + 1 + log β i + 1 ] + γ 0 t 0 γ N ( t ) t c ,
where β k = C ( χ ( t k ) , χ ( t k ) ) , then system (1) is exponentially incrementally stable.

4. Incremental Stability for Time-Delay Dynamical Systems with Discontinuous Right-Hand Sides

Here, we consider the time-delay dynamical systems [39] with discontinuous right-hand sides using multiple norms, formulated as follows:
x ˙ = f ( x ( t ) , x τ 1 ( t ) , , x τ m ( t ) , r ( t ) )
where x , x τ 1 , , x τ m R n , x τ k ( t ) = x ( t τ k ( t ) ) , ( k = 1 , , m ) represents the time-delay term. The right-hand function f may be discontinuous with respect to ( x , x τ 1 , , x τ m ) . The solution of the system (5) can be defined as a solution of the following differential inclusion, which is named a (time-delay) Filippov system,
x ˙ F ( x ( t ) , x τ 1 ( t ) , , x τ m ( t ) , r ( t ) ) , t [ 0 , T ] x ( s ) = ϕ ( s ) , s [ t 0 τ ¯ , t 0 ]
where x τ k ( t ) = x ( t τ k ( t ) ) , τ k ( t ) is non-negative for t [ t 0 , T ] , T ( t 0 , + ] , τ ¯ = max k sup s t 0 τ k ( s ) , the initial function ϕ ( · ) is defined on [ t 0 τ ¯ , t 0 ] , and F is a set-valued mapping defined as follows,
F ( x ( t ) , x τ 1 ( t ) , , x τ m ( t ) , r ( t ) ) = K [ f ] ( x ( t ) , x τ 1 ( t ) , , x τ m ( t ) , r ( t ) ) = ϵ > 0 μ ( P ) = 0 c o ¯ { f ( B ( ( x , x τ 1 , , x τ m ) , ϵ ) P , r ( t ) ) } ,
where μ ( · ) stands for the Lebesgue measure, B ( ( x , x τ 1 , , x τ m ) , ϵ ) = { ( y , y τ 1 , , y τ m ) : | ( y , y τ 1 , , y τ m ) ( x , x τ 1 , , x τ m ) | ϵ } represents the ϵ -neighborhood of ( x , x τ 1 , , x τ m ) with the given vector norm | · | , and c o ¯ represents convex closure.

4.1. Existence and Uniqueness of the Solution

Before the main theorem, the existence and uniqueness of the Cauchy problem of the Filippov system (6) should be proved first. Herein, several of the existing results on dynamical systems without time-delay terms are presented as follows. Readers are referred to [15,37,40] for the details. For some dynamical systems without time-delay terms, formulated as follows,
x ˙ = f ( x , t )
in which x R n , we have the corresponding Filippov system,
x ˙ = F ( x , t ) = K [ f ] ( x , t ) = ϵ > 0 μ ( P ) = 0 c o ¯ { f ( B ( x , ϵ ) P , t ) }
Then, from Definition 4 and 5 in [40]. We have concluded that under the following Assumption 2, it can be guaranteed that system (7) has at least one solution.
Assumption 2 
([40]). The set-valued mapping F : R n × R + R n satisfies that for all ( x , t ) R n × R + , F ( x , t ) is non-empty, bounded, convex, and closed, and F is upper semicontinuous at ( x , t ) .
In Assumption 2, ‘upper semicontinuity’ for the set-valued mapping F is defined as follows.
Definition 8 
(Section 1, Chapter 2 in [37]). A set-valued mapping F : R n × R + Y is called upper semicontinuous at ( x , t ) R n × R + if and only if for any neighborhood U of F ( x , t ) , δ > 0 , such that ( x ˜ , t ˜ ) B ( ( x , t ) , δ ) , F ( x ˜ , t ˜ ) U .
Assumption 3. 
With respect to a given Euclid norm | · | defined on a n-dimensional space, we have the following hypothesis:
1. 
For any ( x , y 1 , , y m , r ) R n × ( m + 1 ) × R , F ( x , y 1 , , y m , r ( t ) ) is non-empty, convex, closed in R n , and set-valued mapping F is upper semicontinuous with respect to ( x , y 1 , , y m , r ) .
2. 
(Linearly increasing) There exists α > 0 such that
sup { | v | : v F ( x , y 1 , , y m , r ) } α ( | x | + | y 1 | + + | y m | + | r | + 1 )
holds for any ( x , y 1 , , y m , r ) R n × ( m + 1 ) × R . With Gronwall inequality [41], it can easily be seen as equivalent to: there exists Θ > 0 such that
sup { | v | : v F ( x , y 1 , , y m , r ) } Θ
holds for any ( x , y 1 , , y m , r ) R n × ( m + 1 ) × R .
3. 
Function τ k : [ t 0 , T ] [ 0 , ) ( k = 1 , , m ) is continuously differentiable and bounded, with its upper bound τ ¯ : = max { τ k ( t ) : t [ t 0 , T ] , k = 1 , , m } and lower bound τ ̲ = min { τ k ( t ) : t [ t 0 , T ] , k = 1 , , m } > 0 .
4. 
The initial function ϕ ( · ) L ( [ t 0 τ ¯ , t 0 ] , R n ) is measurable.
5. 
For any ( x , y 1 , , y m , r ) R n × ( m + 1 ) × R n × R + , there exists continuous function h ( · ) : R + R + , such that | F ( x 1 , y 1 , , y m , r ( t ) ) F ( x 2 , y 1 , , y m , r ( t ) ) | h ( | x 1 x 2 | ) holds.
First, we fix a continuous initial function ϕ ( · ) L ( [ t 0 τ ¯ , t 0 ] , R n ) , then select a measurable function ψ : [ t 0 τ ¯ , t 0 ] R n such that ψ ( s ) K [ ϕ ] ( s ) holds for s [ t 0 τ ¯ , t 0 ] almost everywhere. For t [ t 0 , t 0 + τ ̲ ] , consider the following differential inclusion [42]:
x ˙ ( t ) F ( x ( t ) , ψ ( t τ 1 ( t ) ) , , ψ ( t τ m ( t ) ) , r ( t ) ) x ( 0 ) = ϕ ( 0 ) ,
From Assumption 3, together with Assumption 2 and the conclusion for the existence of the solution of system (8), the inclusion (6) has at least one solution defined in [ t 0 , t 0 + τ ̲ ] .
Therefore, similarly, the solution can be extended to [ t 0 , + ) . Assume, as an inductive step, that the solution x is defined on [ t 0 τ ¯ , t 0 + N τ ̲ ] , for some N = 1 , 2 , . Then, one can consider the vector α = x ( t 0 + N τ ̲ ) as the initial state of the following differential inclusion,
x ˙ ( t ) F ( x ( t ) , ψ ( t τ 1 ( t ) ) , , ψ ( t τ m ( t ) ) , r ( t ) ) x ( 0 ) = ϕ ( 0 ) t [ t 0 + N τ ̲ , t 0 + ( N + 1 ) τ ̲ ] almost everywhere .
Then, one can extend x ( t ) , t [ t 0 τ ¯ , t 0 + N τ ̲ ] to a right neighborhood of t 0 + N τ ̲ , the interval [ t 0 τ ¯ , t 0 + ( N + 1 ) τ ̲ ] . That is, we have the following lemma,
Lemma 2. 
If the time-delay Filippov system (6) satisfies Assumption 3, then, the system (6) has at least on solution in [ t 0 , T ] , T ( 0 , + ] .
Then, when it comes to the problem of uniqueness, enlightened by Chapter 2, Section 10, Theorem 1 in [15], we similarly research the conditions for uniqueness of the solution of system (6).
Theorem 2. 
Suppose that function f ( x ( t ) , x τ 1 ( t ) , , x τ m ( t ) , r ( t ) ) defined on region D R n × R n × R + is discontinuous on zero measure set M, and there exists an integral function l ( t ) such that | f ( x , x τ 1 , , x τ m , r ) | l ( t ) holds for any ( x , x τ 1 , , x τ m , r ) D , and l ( t ) < almost everywhere. Let ϵ 0 > 0 , for any ( x 0 , x τ 1 , , x τ m , r ) , ( y 0 , y τ 1 , , y τ m , r ) satisfying that | x 0 y 0 | < ϵ 0 and | x τ k y τ k | | x 0 y 0 | , ( k = 1 , , m ) ,
( x 0 y 0 ) · ( f ( x 0 , x τ 1 , , x τ m , t ) f ( y 0 , y τ 1 , , y τ m , r ( t ) ) ) l ( t ) | x 0 y 0 | 2
Then, under the simplest convex definition (Page 50 in [15]), equation x ˙ = f ( x , x τ 1 , , x τ m , r ( t ) ) is right-unique on D.
Proof. 
If for all t R + , any ϵ 0 > 0 , x ( · ) and y ( · ) satisfying that | x ( · ) y ( · ) | [ t τ ¯ , t ] | x ( t ) y ( t ) | < ϵ 0 , on the basis of (9), we have
( x ( t ) y ( t ) ) · ( f ( x ( t ) , x τ 1 ( t ) , , x τ m ( t ) , r ( t ) ) f ( y ( t ) , y τ 1 ( t ) , , y τ m ( t ) , r ( t ) ) ) l ( t ) | x ( t ) y ( t ) | 2 ,
holds. Thus, we have
1 2 d d t | x ( t ) y ( t ) | 2 = ( x ( t ) y ( t ) ) ( x ˙ ( t ) y ˙ ( t ) ) = ( x ( t ) y ( t ) ) · ( f ( x ( t ) , x τ 1 ( t ) , , x τ m ( t ) , r ( t ) ) f ( y ( t ) , y τ 1 ( t ) , , y τ m ( t ) , r ( t ) ) ) l ( t ) | x ( t ) y ( t ) | 2 .
Therefore, similar to the proof of Theorem 1, in Section 10, Chapter 2 in [15], it can be seen that
d d t | x ( t ) y ( t ) | 2 e L ( t ) 0
where L ( t ) = t 0 t l ( s ) d s , that is, | x ( t ) y ( t ) | 2 e L ( t ) is decreasing with respect to time t. If | x ( · ) y ( · ) | [ t 0 τ ¯ , t 0 ] = sup { | x ( t ) y ( t ) | : t [ t 0 τ ¯ , t 0 ] } = 0 , then | x ( t ) y ( t ) | = 0 holds for t > t 0 , then we obtained the right uniqueness of the solution of system (6). □
If system (5) is a switched system, according to Theorem 2, one can prove the right uniqueness of the solution.
Here, we formulate a switched system with time delay, that is, the right-hand function is switched with respect to ( x , x τ 1 , , x τ m ) .
f ( x , t ) = f i ( x , x τ 1 , , x τ m , r ( t ) ) , ( x , x τ 1 , , x τ m , t ) R i
in which f i : R n × ( m + 1 ) × R + R n , regions R i R n × ( m + 1 ) × R + , i = 1 , , K . All of the regions R i have non-empty interiors. The discontinuities are composed of several smooth hypersurfaces of dimension d ( d < ( m + 1 ) n + 1 ) . Suppose that { S i } i = 1 N is a sequence ( ( m + 1 ) n ) -dimensional smooth hypersurfaces, S i = { ( x , x τ 1 , , x τ m , t ) R n × ( m + 1 ) × R + ϕ i ( x , x τ 1 , , x τ m , t ) = 0 } , in which ϕ i ( x , x τ 1 , , x τ m , t ) C 1 ( R n × ( m + 1 ) × R + , R ) , the continuous region of the function f is a sequence of connected regions, whose boundaries are the switching surfaces. Suppose that the switching surfaces never intersect each other. Denote one of the connected continuous region of f by G i + ( G i ) , then it satisfies that
1.
G i + ( G i ) k S k ;
2.
f ( x , x τ 1 , , x τ m , r ( t ) ) is continuous in G i + ( G i ) ;
3.
ϕ i ( x , x τ , , x τ m , t ) > 0 ( < 0 ) holds in G i + ( G i ).
in which G i + and G i are two different regions with their common boundary on S i .
Take a switched system defined as (10) with K = 2 , N = 1 as an example.In the domain G, consider one of the switching hypersurface S, f + ( x , x τ 1 , , x τ m , r ( t ) ) and f ( x , x τ 1 , , x τ m , r ( t ) ) represent the limiting values of the function f ( x , x τ 1 , , x τ m , r ( t ) ) at point ( x , x τ 1 , , x τ m , t ) from the regions G + and G , respectively. f N + ( x , x τ 1 , , x τ m , r ( t ) ) and f N ( x , x τ 1 , , x τ m , r ( t ) ) represent the projections of the vectors f + ( x , x τ 1 , , x τ m , r ( t ) ) and f ( x , x τ 1 , , x τ m , r ( t ) ) onto the normal vector to S directed from G + to G , respectively, at the point ( x , x τ 1 , , x τ m , t ) .
Together with Theorem 2, according to bimodal time-delay systems, we have the following conclusion:
Theorem 3. 
Under the notations defined above, for all t [ t 0 , + ) and point ( x , x τ 1 , , x τ m , t ) S , if the inequality f N ( x , x τ 1 , , x τ m , r ( t ) ) f N + ( x , x τ 1 , , x τ m , r ( t ) ) > 0 is fulfilled, then right uniqueness of Filippov solution for the bimodal system (10) ( K = 2 , N = 1 ) occurs. (Possibly different inequalities for different ( x , x τ 1 , , x τ m ) and t.)
Proof. 
We use the conclusion in Theorem 2 to prove Theorem 3 above.
For any point ( z , z τ 1 , , z τ m , t ) on the switching surface S, and ( x , x τ 1 , , x τ m , t ) R 1 , ( y , y τ 1 , , y τ m , t ) R 2 , which satisfy that | y τ k z τ k | | y z | , | x τ k z τ k | | x z | , | x τ k y τ k | | y x | . Since f 1 / x , f 2 / x and f 1 / x τ i , f 2 / x τ i ( i = 1 , , m ) are bounded, this implies that there exist l , k i > 0 ( i = 1 , , m ) such that
| f ( x , x τ 1 , , x τ m , r ( t ) ) f 1 ( z , z τ 1 , , z τ m , r ( t ) ) | l | x z | + i = 1 m k i | x τ i z τ i | , | f ( y , y τ 1 , , y τ m , r ( t ) ) f 2 ( z , z τ 1 , , z τ m , r ( t ) ) | l | y z | + i = 1 m k i | y τ i z τ i | .
Similar to the proof of Theorem 2 in Chapter 2 in [15], if f N ( x , x τ 1 , , x τ m , r ( t ) ) f N + ( x , x τ 1 , , x τ m , r ( t ) ) > 0 holds, then
( x z ) ( f 1 ( z , z τ 1 , , z τ m , r ( t ) ) f 2 ( z , z τ 1 , , z τ m , r ( t ) ) ) 0 .
The inequality (12) still holds if vector x y is substituted with vector x z , which is in the same direction. Together with (11) and (12), it infers that
( x y ) ( f ( x , x τ 1 , , x τ m , r ) + f 2 ( z , z τ 1 , , z τ m , r ) f 1 ( z , z τ 1 , , z τ m , r ) f ( y , y τ 1 , , y τ m , r ) ) | y x | · ( | f ( x , x τ 1 , , x τ m , r ) f 1 ( z , z τ 1 , , z τ m , r ) | + | f ( y , y τ 1 , , y τ m , r ) f 2 ( z , z τ 1 , , z τ m , r ) | ) < | y x | · ( l | x z | + i = 1 m k i | x τ i z τ i | + l | y z | + i = 1 m k i | y τ i z τ i | ) .
That is, if | y τ i z τ i | | y z | , | x τ i z τ i | | x z | , | x τ i y τ i | | y x | ( i = 1 , , m ) , it holds that
( x y ) ( f ( x , x τ 1 , , x τ m , r ( t ) ) f ( y , y τ 1 , , y τ m , r ( t ) ) ) < 2 ( l + i = 1 m k i ) | y x | 2 .
Thus, together with Theorem 2, at any point ( z , z τ 1 , , z τ m , t ) in the domain, right uniqueness of the Filippov solution for system (10) ( K = 2 , N = 1 ) occurs for t [ t 0 , + ) . □

4.2. Criteria for Incremental Stability for Filippov Systems with Time Delay

Here, with the conclusions above, we then research the conditions of incremental stability for time-delay systems with discontinuous right-hand sides.
x ˙ ( t ) = f ( x , x τ 1 , , x τ m , r ( t ) )
where x τ k ( t ) = x ( t τ k ( t ) ) , τ k ( t ) ( k = 1 , 2 , , m ) is a bounded function with respect to time t, the function f ( x , x τ 1 , , x τ m , r ( t ) ) = ( f 1 , f 2 , , f n ) : R n × ( m + 1 ) × R + R n is discontinuous with respect to ( x , x τ 1 , , x τ m , t ) on a zero measure set. r ( t ) is a piecewise right-continuous switched function, with its discontinuities in { t 1 , t 2 , , t j , } .
Then, according to the discontinuous right-hand side f ( x , x τ 1 , , x τ m , r ( t ) ) , with corresponding set-valued mapping F ( x , x τ 1 , , x τ m , r ( t ) ) = K [ f ] ( x , x τ 1 , , x τ m , r ( t ) ) , we construct a sequence of functions { f p ( x , x τ 1 , , x τ m , r ( t ) ) } p = 1 satisfying the following conditions, denoted by Condition C t i m e d e l a y ( Σ ) , where Σ R n × ( m + 1 ) :
1.
f p ( x , x τ 1 , x τ m , r ( t ) ) = ( f 1 p , f 2 p , , f m p ) : R n × ( m + 1 ) × R + R n is continuous and continuously differentiable with respect to ( x , x τ 1 , x τ m ) , and continuous with respect to ( x , x τ 1 , x τ m , r ) . Moreover, f i p , i = 1 , , m satisfies local Lipschitz conditions for ( x , x τ 1 , x τ m ) Σ .
2.
For each t t 0 and compact set Σ R n × ( m + 1 ) ,
lim m d H G r a p h ( f p ( Σ , r ( t ) ) ) , G r a p h ( F ( Σ , r ( t ) ) ) = 0
holds, where F ( x , x τ 1 , x τ m , r ( t ) ) = K [ f ] ( x , x τ 1 , x τ m , r ( t ) ) and d H represents the Hausdorff metric. G r a p h ( F ) and G r a p h ( f p ) are considered on R m × ( n + 1 ) × [ t 0 , + ) , where G r a p h ( F ( x , x τ 1 , x τ m , r ( t ) ) ) = { ( x , t , y ) : y F ( x , x τ 1 , x τ m , r ( t ) ) , ( x , x τ 1 , x τ m ) Σ , t t 0 } , and so it is with G r a p h ( f p ( x , x τ 1 , x τ m , r ( t ) ) ) .
3. 
For any compact set Σ R m × ( n + 1 ) , there exists measure w ( · ) , defined as w ( Σ ) = q ( λ ( Σ ) ) , in which λ represents the Lebesgue measure, q is a measurable function mapping R + to R + , such that | f p ( x , x τ 1 , x τ m , t ) | w ( Σ ) holds for each ( x , x τ 1 , x τ m ) Σ and t t 0 .
Thus, we have the following conclusion on incremental stability of time-delay systems with discontinuous right-hand sides:
Theorem 4. 
Suppose that system (13) has a unique solution, and there exists a sequence of functions { f p } p = 1 satisfying Condition C t i m e d e l a y ( Σ ) . χ ( t ) is a right-continuous staircase function with discontinuities in { t 1 , , t j , } . Let N ( t ) = # { j : t t j , j = 1 , 2 , } , and suppose that there exist
1. 
A constant T 0 > 0 ,
2. 
Positive constants γ k ( k = 0 , 1 , ) , γ 0 γ 1 ,
3. 
C ( χ ( t ) , χ ( t ) ) C 0 for any t , t t 0 ,
such that the following conditions hold:
Let A p ( t ) = ( f p / x ) ( x , x τ 1 , , x τ m , r ( t ) ) , B k p ( t ) = ( f p / x τ k ) ( x , x τ 1 , , x τ m , r ( t ) ) , k = 1 , , m , there exist matrix-valued functions B 1 k p ( t ) , B 2 k p ( t ) , B 1 k p ( t ) + B 2 k p ( t ) = B k p ( t ) such that
ν χ ( t ) ( A p ( t ) + k = 1 m B 1 k p ( t ) ) + k = 1 m τ k ( t ) B 1 k p ( t ) χ ( t ) ( A ˜ k p ( t ) + B ˜ k p ( t ) ) + k = 1 m B 2 k p ( t ) χ ( t ) exp γ N ( t ) t exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) γ N ( t )
where
A ˜ k p ( t ) = sup t τ k ( t ) s t A p ( s ) χ ( t ) exp γ N ( t ) t γ N ( s ) s B ˜ k p ( t ) = sup t τ k ( t ) s t B k p ( s ) χ ( t ) exp γ N ( t ) t γ N ( s τ k ( s ) ) ( s τ k ( s ) ) .
and for any t > T 0 ,
1 t i = 0 N ( t ) 1 [ ( γ i + 1 γ i ) t i + 1 + log β i + 1 ] + γ 0 t 0 γ N ( t ) t c ,
where β k = C ( χ ( t k ) , χ ( t k ) ) , then system (6) is exponentially incrementally stable in Σ for t [ t 0 , + ) .
Proof. 
First, we construct a sequence of functions satisfying Condition C t i m e d e l a y ( Σ ) , and the corresponding sequence of time-delay system:
x ˙ ( t ) = f p ( x , x τ 1 , x τ m , r ( t ) )
in which, according to Condition C t i m e d e l a y ( Σ ) , f i , i = 1 , , m satisfy local Lipschitz conditions, so it can be seen that (14) has a unique solution for t R + .
Therefore, from Theorem 1, for each p, system (14) is exponentially incrementally stable.
Denote any two of the solutions of (6) with different initial function x 0 ( · ) and y 0 ( · ) by x ( t ) and y ( t ) . Then, here, try to approximate x ( t ) and y ( t ) by two sequences of solutions of (14), denoted by x p and y p , respectively. The initial functions x 0 ( s ) and y 0 ( s ) here are defined for s [ t 0 τ ¯ , t 0 ] , where τ ¯ = max k sup s [ t 0 , + ) τ k ( s ) .
With f p satisfying condition C t i m e d e l a y ( Σ ) , for each given T, t 0 t T , we have x p ( t ) and x ˙ p ( t ) are bounded regarding ( x 0 ( · ) , T ) . It is the same with y p ( t ) . Thus, x p ( t ) ( y p ( t ) ) is uniformly bounded on [ t 0 , T ] , and x ˙ p ( t ) ( y ˙ p ( t ) ) is uniformly bounded on [ t 0 , T ] as well.
Because of Condition 1 in C t i m e d e l a y ( Σ ) , x p ( t ) ( y p ( t ) ) is continuous with respect to t. Together with Condition 3 in C t i m e d e l a y ( Σ ) , we conclude that x p ( t ) ( y p ( t ) ), p N + are equicontinuous for t [ t 0 , T ] .
Here, we present the Arzela–Ascoli lemma (similar to Theorem 2.2 in [43]):
Lemma 3. 
(Arzela–Ascoli Lemma) X is a compact set on R n . If a sequence { f n } 1 in C ( X ) is bounded and equicontinuous, then it has a uniformly convergent subsequence.
From Lemma 3, one can find a sub-sequence of { x p ( t ) } p N + and { y p ( t ) } p N + (still denoted by x m ( t ) and y m ( t ) ) satisfying that x p ( t ) ( y p (t)) uniformly converges to a continuous function x * ( t ) ( y * ( t ) ) on [ t 0 , T ] , T ( t 0 , + ] . For all k N + , one can find a subsequence { x k j ( t ) } j N + of { x p ( t ) } p N + such that | x k j ( t ) x * ( t ) | < 1 j holds on the interval [ t 0 , t 0 + k + 1 ] . Then, by diagonal selection principle, select a new subsequence { x k k ( t ) } k N + such that { x k k ( t ) } k N + uniformly converges to a continuous function x * ( t ) on [ t 0 , T ] , T ( t 0 , + ] . It is the same with x τ k * ( t ) , for k = 1 , , m .
The system x ˙ p ( t ) = f p ( x , x τ 1 , , x τ m , r ( t ) ) has a unique solution. x p ( t ) satisfies Lipschitz condition:
| x p ( t ) x p ( t ) | χ ( t 0 ) L | t t | χ ( t 0 )
where t , t [ t 0 , T ] , L > 0 . Because of norm equivalence, the Lipschitz condition above also holds with other norms defined in R n . Therefore, x * ( t ) ( y * ( t ) ) also satisfies the Lipschitz condition, that is, x ˙ * ( t ) ( y ˙ * ( t ) ) exists and is bounded and measurable for [ t 0 , T ) , T ( t 0 , + ] .
We then have the conclusion that x ˙ p ( t ) ( y ˙ p ( t ) ) weakly converges to x ˙ * ( t ) ( y ˙ * ( t ) ) on the space L 1 ( [ t 0 , T ] , R n ) , the demonstrations are as follows.
C 0 ( [ t 0 , T ] , R n ) is dense in the Banach space L ( [ t 0 , T ] , R n ) , which is the conjugate space L 1 ( [ t 0 , T ] , R n ) . Therefore, the following equation
t 0 T x ˙ p ( t ) x ˙ * , q ( t ) d t = t 0 T q ˙ ( t ) , x p ( t ) x * d t .
holds for each q ( t ) C 0 ( [ t 0 , T ] , R n ) . Since { x ˙ p ( t ) } is bounded for each p, from the Lebesgue-dominant convergence theorem we have
lim p t 0 T x ˙ p ( t ) x ˙ * , q ( t ) d t = t 0 T q ˙ ( t ) , lim p x p ( t ) x * ( t ) d t = 0 .
That is, { x ˙ p ( t ) } weakly converges to x ˙ * ( t ) on the space L 1 ( [ t 0 , T ] , R n ) .
With Mazur’s convexity theorem [43], one can find a l n ( b l n ) with l = 1 m a l m = 1 ( l = 1 m b l m = 1 ) such that x ˜ ˙ p converges to x ˙ * ( t ) almost everywhere on [ t 0 , T ] , where x ˜ p ( t ) = l = 1 p a l p x p . Notice that x ˜ p is in the convex closure of { x p } , x ˜ ˙ p converges to x ˙ * uniformly. So it is with y ˜ ˙ p ( t ) with y ˜ p ( t ) = l = 1 p b l p y p , and it is the same with { x ˙ τ k p ( t ) } and { x ˙ τ k * ( t ) } .
Recall Condition 3 in C t i m e d e l a y ( Σ ) . For Σ R n , it holds that
lim p d H G r a p h ( f p ( Σ , t ) ) , G r a p h ( F ( Σ , t ) ) = 0 , t t 0 ,
With x ˜ ˙ p ( t ) in the convex closure of { f p ( x p , x τ 1 p , , x τ m p , r ( t ) ) } , for any ϵ > 0 , there exists N > 0 such that x ˙ p ( t ) B ( F ( x * , x τ 1 * , , x τ m * , r ( t ) ) , ϵ ) for all p > N and ( x * , x τ 1 * , , x τ m * ) , ( x p , x τ 1 p , , x τ m p ) Σ , t [ t 0 , T ] .
Since ϵ can be arbitrarily small, it can be seen that x ˙ * ( t ) F ( x * , x τ 1 * , , x τ m * , r ( t ) ) with x * Σ , which infers that the solution of (6) equals to x * ( t ) in Σ almost everywhere on [ t 0 , T ] . For x ( t ) and x * , both are continuous because x ˜ p converges to x * uniformly on [ t 0 , T ] , it can be seen that x * is the solution of (6). So it is with y * ( t ) .
Because system (6) has a unique solution, x ( t ) = x * ( t ) and y ( t ) = y * ( t ) almost everywhere for t t 0 . That is, x p converges to x ( t ) uniformly in [ t 0 , T ] , T ( t 0 , + ] . A similar proof can be applied to y p ( t ) and y ( t ) .
For { f p ( x , x τ 1 , , x τ m , r ( t ) ) } is exponentially incrementally stable from Theorem 1, there exists some M > 0 and α > 0 ,
| x p ( t ) y p ( t ) | χ ( t 0 ) M e α ( t t 0 ) sup s [ t 0 τ ¯ , t 0 ] | x 0 ( s ) y 0 ( s ) | χ ( t 0 )
for all p N + and t t 0 . For each given t t 0 , let ϵ ( t ) = 3 M e α ( t t 0 ) sup s [ t 0 τ ¯ , t 0 ] | x 0 ( s ) y 0 ( s ) | χ ( t 0 ) , there exists some p 0 ( ϵ , t ) with which | x p ( t ) x ( t ) | χ ( t 0 ) ϵ / 3 and | y p ( t ) y ( t ) | χ ( t 0 ) ϵ / 3 hold for p p 0 ( ϵ , t ) , which implies that
| x ( t ) y ( t ) | χ ( t 0 ) | x ( t ) x p ( t ) | χ ( t 0 ) + | y ( t ) y p ( t ) | χ ( t 0 ) + | x p ( t ) y p ( t ) | χ ( t 0 ) ϵ = 3 M e α ( t t 0 ) sup s [ t 0 τ ¯ , t 0 ] | x 0 ( s ) y 0 ( s ) | χ ( t 0 ) .
This completes the proof. □
Remark 2. 
Similar to Theorem 4, Corollary 1 and Corollary 2 can also be extended to discontinuous cases.

5. Applications

In this section, with the conclusion in Theorem 4, the applications to switched time-delay systems and Hopfield time-delay systems are given.

5.1. Linear Switched Time-Delay System

Consider the following linear switched time-delay system:
x ˙ ( t ) = f ( x , x τ 1 , , x τ m , t ) = A 1 ( t ) x ( t ) + k = 1 m B 1 k ( t ) x ( t τ k ( t ) ) + J 1 ( t ) , ϕ ( x , x τ 1 , , x τ m ) > 0 A 2 ( t ) x ( t ) + k = 1 m B 2 k ( t ) x ( t τ k ( t ) ) + J 2 ( t ) , ϕ ( x , x τ 1 , , x τ m ) < 0 .
in which x ( t ) R n , A 1 ( t ) , A 2 ( t ) R n × n , B 1 k ( t ) , B 2 k ( t ) ( k = 1 , , m ) R n × n is piecewise continuous, bounded matrix-valued functions, whose discontinuities are in { t 1 , t 2 . t i , } . Take a simple bimodal system as an example, with the switching surface S = { ( x , x τ 1 , , x τ m ) : ϕ ( x , x τ 1 , , x τ m ) = 0 } .
In order to guarantee the uniqueness of the Filippov solution for system (15), the linear time-delay system satisfies the following hypothesis:
Assumption 4. 
The linear time-delay system (15) satisfies:
1. 
The right-hand function of system (15) satisfies Assumption 3.
2. 
For each point ( x , x τ 1 , , x τ m ) S , the time-delay system (15) satisfies
ϕ x ( x , x τ 1 , , x τ m ) · ( A 1 ( t ) A 2 ( t ) ) x + k = 1 m ( B 1 k ( t ) B 2 k ( t ) ) x τ k + J 1 ( t ) J 2 ( t ) < 0
Under the definitions and Assumption 4 above, from Lemma 2 and Theorem 3, it can be obtained that system (15) has a unique solution.
Assumption 5. 
Suppose that in the continuous regions, the right-hand side f ( x , x τ 1 , , x τ m , t ) satisfies the conditions for exponential incremental stability in Theorem 1, which are formulated as follows.
Let N ( t ) = # { j : t t j , j = 1 , 2 , } , χ ( t ) be a right-continuous staircase function with its discontinuous points belonging to { t 1 , , t j , } , and there exist
1. 
A constant T 0 > 0 ,
2. 
Positive constants γ k ( k = 0 , 1 , ) , γ 0 γ 1 ,
3. 
C ( χ ( t ) , χ ( t ) ) C 0 for any t , t t 0 ,
4. 
Matrix B 1 k ( 1 ) , B 1 k ( 2 ) , B 2 k ( 1 ) , B 2 k ( 2 ) R n × n satisfying that B 1 k ( 1 ) ( t ) + B 1 k ( 2 ) ( t ) = B 1 k ( t ) and B 2 k ( 1 ) ( t ) + B 2 k ( 2 ) ( t ) = B 2 k ( t ) ,
such that the following conditions are satisfied:
L 1 ( t ) = ν χ ( t ) ( A 1 ( t ) + k = 1 m B 1 k ( 1 ) ( t ) ) + k = 1 m τ k ( t ) ( A ˜ 1 k ( t ) + B ˜ 1 k ( t ) ) B 1 k ( 1 ) ( t ) χ ( t ) + k = 1 m B 1 k ( 2 ) ( t ) χ ( t ) exp γ N ( t ) t exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) γ N ( t ) L 2 ( t ) = ν χ ( t ) ( A 2 ( t ) + k = 1 m B 2 k ( 1 ) ( t ) ) + k = 1 m τ k ( t ) ( A ˜ 2 k ( t ) + B ˜ 2 k ( t ) ) B 2 k ( 1 ) ( t ) χ ( t ) + k = 1 m B 2 k ( 2 ) ( t ) χ ( t ) exp γ N ( t ) t exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) γ N ( t )
in which, for i = 1 , 2 ,
A ˜ i k ( t ) = sup t τ k ( t ) s t A i ( s ) χ ( t ) exp γ N ( t ) t γ N ( s ) s B ˜ i k ( t ) = sup t τ k ( t ) s t B i k ( s ) χ ( t ) exp γ N ( t ) t γ N ( s τ k ( s ) ) ( s τ k ( s ) )
and for all t > T 0 ,
1 t i = 0 N ( t ) 1 [ ( γ i + 1 γ i ) t i + 1 + log β i + 1 ] + γ 0 t 0 γ N ( t ) t c ,
where β k = C ( χ ( t k ) , χ ( t k ) ) .
Thus, together with Theorem 4, we have the following corollary on incremental stability for bimodal linear time-delay systems.
Corollary 3. 
Suppose that time-delay system (15) satisfies Assumptions 4 and 5. Meanwhile, in the neighborhood of switching surface S, that is, when δ / 2 < ϕ ( x , x τ 1 , , x τ m ) < δ / 2 , let w = ( A 1 ( t ) A 2 ( t ) ) x ( t ) + k = 1 m ( B 1 k ( t ) B 2 k ( t ) ) x ( t τ k ( t ) ) + J 1 ( t ) J 2 ( t ) , if it is satisfied that
p L 1 ( t ) + ( 1 p ) L 2 ( t ) + 1 δ ν χ ( t ) w · ϕ x + 1 δ k = 1 m w · ϕ x τ k χ ( t ) exp ( γ N ( t ) t ) exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) + p ( 1 p ) k = 1 m τ k ( t ) [ ( B 1 k ( 1 ) χ ( t ) B 2 k ( 1 ) χ ( t ) ) ( A ˜ 2 k ( t ) A ˜ 1 k ( t ) + B ˜ 2 k ( t ) B ˜ 1 k ( t ) ) + 1 δ p B 1 k ( 1 ) + ( 1 p ) B 2 k ( 1 ) χ ( t ) sup t τ k ( t ) s t w · ϕ x χ ( t ) exp γ N ( t ) t exp γ N ( s ) s + w · ϕ x τ k χ ( t ) exp γ N ( t ) t exp γ N ( s τ k ( s ) ) ( s τ k ( s ) ) ] γ N ( t )
holds for δ 2 < ϕ ( x , x τ 1 , , x τ m ) < δ 2 and p = σ ( ϕ ( x , x τ 1 , , x τ m ) / δ ) [ 0 , 1 ] , then the linear switched time-delay system (15) is exponentially incrementally stable.
Proof. 
With the conclusion in Theorem 4, by constructing a sequence of time-delay systems with continuous right-hand sides satisfying Condition C t i m e d e l a y ( Σ ) as follows, we can prove the exponential incremental stability property of system (15):
x ˙ = f δ = σ ( ϕ ( x , x τ 1 , , x τ m ) δ ) [ A 1 ( t ) x ( t ) + k = 1 m B 1 k ( t ) x ( t τ k ( t ) ) + J 1 ( t ) ] + ( 1 σ ( ϕ ( x , x τ 1 , , x τ m ) δ ) ) [ A 2 ( t ) x ( t ) + k = 1 m B 2 k ( t ) x ( t τ k ( t ) ) + J 2 ( t ) ]
in which σ ( · ) is defined as
σ ( ρ ) = 1 , ρ > 1 / 2 , ρ + 1 / 2 , ρ [ 1 / 2 , 1 / 2 ] , 0 , ρ < 1 / 2 .
let w = ( A 1 ( t ) A 2 ( t ) ) x ( t ) + k = 1 m ( B 1 k ( t ) B 2 k ( t ) ) x ( t τ k ( t ) ) + J 1 ( t ) J 2 ( t ) and p = σ ( ϕ ( x , x τ 1 , , x τ m ) / δ ) , the partial derivative of the right-hand function f δ with respect to x and x τ k are as follows,
f δ ( x , x τ 1 , , x τ m , r ( t ) ) x = σ ( ϕ ( x , x τ 1 , , x τ m ) δ ) A 1 ( t ) + [ 1 σ ( ϕ ( x , x τ 1 , , x τ m ) δ ) ] A 2 ( t ) + 1 δ w · ϕ x f δ ( x , x τ 1 , , x τ m , r ( t ) ) x τ k = σ ( ϕ ( x , x τ 1 , , x τ m ) δ ) B 1 k ( t ) + [ 1 σ ( ϕ ( x , x τ 1 , , x τ m ) δ ) ] B 2 k ( t ) + 1 δ w · ϕ x τ k
Thus, it needs to satisfy
L ( f δ , t ) = ν χ ( t ) p A 1 + ( 1 p ) A 2 + 1 δ w · ϕ x + k = 1 m [ p B 1 k ( 1 ) + ( 1 p ) B 2 k ( 1 ) ] + k = 1 m τ k ( t ) ( U ˜ k ( t ) + V ˜ k ( t ) ) p B 1 k ( 1 ) + ( 1 p ) B 2 k ( 1 ) χ ( t ) + k = 1 m p B 1 k ( 2 ) + ( 1 p ) B 2 k ( 2 ) + 1 δ w · ϕ x τ k χ ( t ) exp γ N ( t ) t exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) γ N ( t )
in which
U ˜ k ( t ) = sup t τ k ( t ) s t f δ x ( x ( s ) , x τ 1 ( s ) , , x τ m ( s ) , s ) χ ( t j ) exp γ N ( t ) t γ N ( s ) s V ˜ k ( t ) = sup t τ k ( t ) s t f δ x τ ( x ( s ) , x τ 1 ( s ) , , x τ m ( s ) , s ) χ ( t j ) exp ( γ N ( t ) t γ N ( s τ k ( s ) ) ( s τ k ( s ) ) )
where t [ t j , t j + 1 ) . That is,
L ( f δ , t ) p L 1 ( t ) + [ 1 p ] L 2 ( t ) + 1 δ ν χ ( t ) w · ϕ x + 1 δ k = 1 m w · ϕ x τ k χ ( t ) exp γ N ( t ) t exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) + p ( 1 p ) k = 1 m τ k ( t ) [ ( B 1 k ( 1 ) χ ( t ) B 2 k ( 1 ) χ ( t ) ) ( A ˜ 2 k ( t ) A ˜ 1 k ( t ) + B ˜ 2 k ( t ) B ˜ 1 k ( t ) ) + 1 δ p B 1 k ( 1 ) + ( 1 p ) B 2 k ( 1 ) χ ( t ) sup t τ k ( t ) s t w · ϕ x χ ( t ) exp γ N ( t ) t exp γ N ( s ) s + w · ϕ x τ k χ ( t ) exp γ N ( t ) t exp γ N ( s τ k ( s ) ) ( s τ k ( s ) ) ]
Together with (17), it infers that L ( f δ , t ) γ N ( t ) , that is, system (18) is exponentially incrementally stable for each δ = 1 m , m N + . From Theorem 4, it can be proved that the time-delay system is exponentially incrementally stable. □
Here, let B 1 k ( 1 ) = B 2 k ( 1 ) = 0 , ( k = 1 , , m ) , we have the following corollary, which is a variant of Corollary 3.
Corollary 4. 
Let N ( t ) = # { j : t t j , j = 1 , 2 , } , and suppose that system (15) satisfies Assumption 4, and there exist
1. 
A constant T 0 > 0 ,
2. 
Positive constants γ k ( k = 0 , 1 , ) , γ 0 γ 1 ,
3. 
C ( χ ( t ) , χ ( t ) ) C 0 for any t , t t 0 ,
such that the following conditions are satisfied:
L 1 ( t ) = ν χ ( t ) ( A 1 ( t ) ) + k = 1 m B 1 k ( t ) χ ( t ) exp γ N ( t ) t exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) γ N ( t ) L 2 ( t ) = ν χ ( t ) ( A 2 ( t ) ) + k = 1 m B 2 k ( t ) χ ( t ) exp γ N ( t ) t exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) γ N ( t )
and for all t > T 0 ,
1 t i = 0 N ( t ) 1 [ ( γ i + 1 γ i ) t i + 1 + log β i + 1 ] + γ 0 t 0 γ N ( t ) t c ,
where β k = C ( χ ( t k ) , χ ( t k ) ) . Moreover, in the neighborhood of switching surface S, that is, when δ / 2 < ϕ ( x , x τ 1 , , x τ m ) < δ / 2 , let w = ( A 1 ( t ) A 2 ( t ) ) x ( t ) + k = 1 m ( B 1 k ( t ) B 2 k ( t ) ) x ( t τ k ( t ) ) + J 1 ( t ) J 2 ( t ) , if it is satisfied that
p L 1 ( t ) + ( 1 p ) L 2 ( t ) + 1 δ ν χ ( t ) w · ϕ x + 1 δ k = 1 m w · ϕ x τ k χ ( t ) exp γ N ( t ) t exp γ N ( t τ k ( t ) ) ( t τ k ( t ) ) γ N ( t )
holds for δ / 2 < ϕ ( x , x τ 1 , , x τ m ) < δ / 2 and p = ϕ ( x , x τ 1 , , x τ m ) / δ [ 0 , 1 ] , then the linear switched time-delay system (15) is exponentially incrementally stable.

5.2. Hopfield Neural Network Systems with Time Delay

Consider the following Hopfield neural network system with time delay:
x ˙ ( t ) = D ( t ) x ( t ) + T ( t ) g ( x ( t ) ) + S ( t ) u ( x τ ( t ) ) ) + J ( t )
where x = ( x 1 , x 2 , , x n ) is the state variable, the time-delay term x τ ( t ) = x ( t τ ( t ) ) . For any t R + , D ( t ) = d i a g { d 1 ( t ) , , d n ( t ) } , T ( t ) = ( T i j ( t ) ) R n , n , S ( t ) = ( S i j ( t ) ) R n , n , J = ( J 1 , J 2 , , J n ) R n is the input vector, g ( x ) = ( g 1 ( x 1 ) , g 2 ( x 2 ) , , g n ( x n ) ) , u ( x ) = ( u 1 ( x 1 ) , u 2 ( x 2 ) , , u n ( x n ) ) .
Here, we list the following hypothesis, denoted by Condition C 2 :
1.
There exists D ˜ = d i a g { D 1 , D 2 , , D n } , D i > 0 , such that d i ( · ) is continuous and d i ( ζ 1 ) d i ( ζ 2 ) ζ 1 ζ 2 D i holds for i = 1 , 2 , , n and ζ 1 ζ 2 .
2.
g i ( · ) is non-decreasing and non-trivial in any compact set in R , and each g i ( · ) has only finite discontinuous points. Therefore, in any compact set in R , except a finite points ρ k , where there exist finite right and left limits g i ( ρ k + ) and g i ( ρ k ) with g i ( ρ k + ) > g i ( ρ k ) , g i ( · ) is continuous.
3.
u i ( · ) is non-decreasing and non-trivial in any compact set in R , and each u i ( · ) has only finite discontinuous points. Therefore, in any compact set in R , except a finite points η k , where there exist finite right and left limits u i ( η k + ) and u i ( η k ) with u i ( η k + ) > u i ( η k ) , u i ( · ) is continuous.
4.
Here, define a matrix measure ν ξ , 1 ( A ) = max j [ a j j + i j | ξ i ξ j 1 a i j | ] for matrix A = ( a i j ) , with respect to vector norm | · | ξ , 1 and matrix norm A ξ , 1 = ξ A ξ 1 1 , where ξ = d i a g { ξ 1 , , ξ n } . There exists a positive diagonal matrix ξ ( r ( t ) ) = d i a g { ξ 1 ( r ( t ) ) , ,
ξ n ( r ( t ) ) } such that
ν ξ ( r ( t ) ) , 1 ( T ) = max j { T j j ( t ) + i j | ξ i ( r ( t ) ) ξ j ( r ( t ) ) 1 T i j ( t ) | } 0
holds for t R + .
Therefore, we have the following corollary:
Corollary 5. 
Suppose the system (19) has a unique solution for t R + , and satisfies Condition C 2 above. Let N ( t ) = # { j : t t j , j = 1 , 2 , } , and there exists
1. 
Positive piecewise right-continuous function m ( t ) > 0 ,
2. 
A constant T 0 > 0 ,
3. 
Positive constants α k > 0 ( k = 0 , 1 , ) ,
4. 
C ( r ( t ) , r ( t ) ) C 0 for any t , t t 0 ,
such that in the continuous regions of g and u,
ν ξ ( r ( t ) ) , 1 ( D ( t ) + T ( t ) g x ( x , x τ ) ) + exp γ N ( t ) t exp γ N ( t τ ( t ) ) ( t τ ( t ) ) S ( t ) u x τ ( x , x τ ) ξ ( r ( t ) ) , 1 γ N ( t )
Let K i = min k [ g i ( ρ k + ) g i ( ρ k ) ] , M i = max k [ u i ( γ k + ) u i ( γ k ) ] , and K = min i K i , M = max i M i ,
K · ν ξ ( r ( t ) ) , 1 ( T ( t ) ) + M · exp γ N ( t ) t exp γ N ( t τ ( t ) ) ( t τ ( t ) ) S ( t ) ξ ( r ( t ) ) , 1 < 0
For any t > T 0 ,
1 t i = 0 N ( t ) 1 [ ( γ i + 1 γ i ) t i + 1 + log β i + 1 ] + γ 0 t 0 γ N ( t ) t c ,
where β k = C ( r ( t k ) , r ( t k ) ) , then system (19) is exponentially incrementally stable.
Proof. 
Suppose that system (19) has a unique solution for t R + , then we construct a sequence of ‘continuous systems’ as follows,
x ˙ ( t ) = f δ ( x , x τ , t ) = D ( t ) x ( t ) + T ( t ) g ˜ ( x ( t ) ) + S ( t ) u ˜ ( x τ ( t ) ) )
where g ˜ ( x ) = ( g ˜ 1 ( x 1 ) , g ˜ 2 ( x 2 ) , , g ˜ n ( x n ) ) , u ˜ ( x ) = ( u ˜ 1 ( x 1 ) , u ˜ 2 ( x 2 ) , , u ˜ n ( x n ) ) . For each i, denote one of the discontinuous points of g i ( x ) by ρ i , and one of the discontinuous point of u i ( x ) by η i .
Function g ˜ δ ( x ) is formulated as follows, if x i [ ρ i δ 2 , ρ i + δ 2 ] , g ˜ i δ ( x i ) = g i ( x i ) , and if x i [ ρ i δ 2 , ρ i + δ 2 ] ,
g i ˜ δ ( x i ) = g i ( ρ i + δ 2 ) g i ( ρ i δ 2 ) δ [ x i ρ i + δ 2 ] + g i ( ρ i δ 2 ) .
Function u ˜ ( x ) is similarly constructed. It can be seen that when δ 0 , the function sequence f δ ( x , x τ , t ) converges to the Filippov differential inclusion of the right-hand side of (19), that is, { f δ ( x , x τ , t ) } satisfies Condition C t i m e d e l a y ( R n × 2 ) . Let g ¯ i ( ρ i , δ ) = g i ( ρ i + δ 2 ) g i ( ρ i δ 2 ) δ , u ¯ i ( η i , δ ) = u i ( η i + δ 2 ) u i ( η i δ 2 ) δ . Here, denote the Jacobi matrix of functions g ˜ δ ( · ) and u ˜ δ ( · ) by matrix G δ and U δ . In the neighborhood of the discontinuous point,
G δ = g ˜ ( x ) x x = ( ρ 1 , ρ 2 , , ρ n ) = d i a g { g ¯ 1 ( ρ 1 , δ ) , g ¯ 2 ( ρ 2 , δ ) , , g ¯ n ( ρ n , δ ) } U δ = u ˜ ( x τ ) x τ x τ = ( η 1 , η 2 , , η n ) = d i a g { u ¯ 1 ( η 1 , δ ) , u ¯ 2 ( η 2 , δ ) , , u ¯ n ( η n , δ ) }
that is, we have
f δ ( x , x τ , t ) x = D ( t ) + T ( t ) G δ , f δ ( x , x τ , t ) x τ = S ( t ) U δ ,
In the continuous regions, the function sequence { f δ ( x , x τ , t ) } satisfies (20), and ν ξ ( r ( t ) ) , 1 ( T ( t ) ) < 0 . Meanwhile, in the neighborhood of the discontinuities, together with the condition (21), it holds that,
ν ξ ( r ( t ) ) , 1 ( D ( t ) + T ( t ) G δ ) + exp γ N ( t ) t exp γ N ( t τ ( t ) ) ( t τ ( t ) ) S ( t ) U δ ξ ( r ( t ) ) , 1 ν ξ ( r ( t ) ) , 1 ( D ( t ) ) + ν ξ ( r ( t ) ) , 1 ( T ( t ) G δ ) + exp γ N ( t ) t exp γ N ( t τ ( t ) ) ( t τ ( t ) ) S ( t ) U δ ξ ( r ( t ) ) , 1 ν ξ ( r ( t ) ) , 1 ( D ( t ) ) + min i g ¯ i ( ρ i , δ ) ν ξ ( r ( t ) ) , 1 ( T ( t ) ) + max i u ¯ i ( η i , δ ) exp γ N ( t ) t exp γ N ( t τ ( t ) ) ( t τ ( t ) ) S ( t ) ξ ( r ( t ) ) , 1 γ N ( t )
for δ 0 . Therefore, if there exists a large enough number m 0 > 0 , such that for each m > m 0 , m N + ,
ν ξ ( r ( t ) ) , 1 ( D ( t ) + T ( t ) G δ ) + exp γ N ( t ) t exp γ N ( t τ ( t ) ) ( t τ ( t ) ) S ( t ) U δ ξ ( r ( t ) ) , 1 γ N ( t )
holds for δ = 1 / m . According to Theorems 1 and 4, the time-delay system (19) is exponentially incrementally stable. □

6. Numerical Experiments

6.1. Linear Time-Delay System

Consider a linear switched time-delay system formulated as follows:
x ˙ ( t ) = A 1 x ( t ) + B 1 x ( t τ ( t ) ) + J 1 , z x > 0 A 2 x ( t ) + B 2 x ( t τ ( t ) ) + J 2 , z x < 0
where x , x τ R 2 . The switching surface of system (22) is { x : ϕ ( x ) = z x = 0 } where z = [ 1 , 3 ] .
Let τ ( t ) = 1 and
A 1 = 6 1 0 7 , A 2 = 7 2 1 10 , B 1 = B 2 = 1 0 0 1 , J 1 = 1 2 , J 2 = 0 1 .
Here, we first prove that system (22) has a unique solution for t [ 1 , + ) . On the switching surface { x : ϕ ( x ) = z x = 0 } , we have
d ϕ d x ( x ) · ( A 1 A 2 ) x + J 1 J 2 = [ 1 , 3 ] 1 3 1 3 x ( t ) + 1 3 = 4 ,
According to Assumption 4, the uniqueness of the Filippov solution for (22) occurs for t [ 0 , + ) . Let γ N ( t ) = 1 1 / 2 [ t ] and the number [ t ] represents the floor of time t. So, for t > 1 ,
L 1 ( t ) = ν 2 ( A 1 ) + B 1 2 e γ N ( t ) t e γ N ( t ) 1 ( t 1 ) = 5.75 + e ( 1 1 / 2 [ t ] ) t e ( 1 1 / 2 [ t 1 ] ) ( t 1 ) < 3.5 < 1 + 1 2 [ t ] L 2 ( t ) = ν 2 ( A 2 ) + B 2 2 e γ N ( t ) t e γ N ( t ) 1 ( t 1 ) = 6.35 + e ( 1 1 / 2 [ t ] ) t e ( 1 1 / 2 [ t 1 ] ) ( t 1 ) < 4 < 1 + 1 2 [ t ]
where · 2 stands for 2-norm. Then, for δ / 2 < ϕ ( x ) < δ / 2 and p = ϕ ( x ) / δ + 1 / 2 [ 0 , 1 ] ,
p L 1 ( t ) + ( 1 p ) L 2 ( t ) + 1 δ ν 2 1 3 1 3 x ( t ) + 1 3 · [ 1 , 3 ] < 1.5 < γ N ( t ) ,
and for all t > T 0 ,
1 t i = 0 N ( t ) 1 [ ( γ i + 1 γ i ) t i + 1 + log β i + 1 ] + γ 0 t 0 γ N ( t ) t 1 t i = 0 N ( t ) 1 [ 1 2 i + 1 ( i + 1 ) ] ( 1 1 2 [ t ] ) t 0.5
Thus, from Corollary 4, the switched time-delay system (22) is exponentially incrementally asymptotically stable. With the initial state x 0 ( s ) = [ 5 × ( 0.5 + s ) 2 , 6 × ( 1 s ) 1.5 ] and y 0 ( s ) = [ 3 × ( 1 + s ) 3 2 , 3 × ( 1 s ) ] for s [ 0 , 1 ] , the corresponding solution of system (22) are x ( t ) = ( x 1 ( t ) , x 2 ( t ) ) and y ( t ) = ( y 1 ( t ) , y 2 ( t ) ) , respectively.
Figure 1 shows the dynamical trajectories of two of the solutions with the initial function defined as x 0 ( · ) and y 0 ( · ) for system (22). Moreover, Figure 2 shows the errors between the two dynamical trajectories of their segments.

6.2. Hopfield Neural Network with Time Delay

Here, we take an example of the Hopfield system with time delay for illustration. The system is formulated as follows,
x ˙ ( t ) = D ( t ) x ( t ) + T ( t ) g ( x ( t ) ) + σ ( t ) S ( t ) u ( x τ ( t ) ) ) + J ( t )
where x = [ x 1 , x 2 ] is the state vector, σ ( t ) is the switched function with respect to time t, takes value between 0 and 1. The parameter matrix are
σ ( t ) = 0 t [ k T 0 , k T 0 + 1 2 T 0 ) 1 t [ k T 0 + 1 2 T 0 , ( k + 1 ) T 0 ) , τ ( t ) = t t 0.2 0.1 × | sin ( π t / 0.4 ) | + 0.1 t > 0.2
D ( t ) = 0.3 + 0.1 × sin ( t ) 0 0 0.3 + 0.1 × cos ( t ) , J ( t ) = 10 sin ( 2 t ) 10 sin ( 2 t ) ,
T ( t ) = 5 sin ( t ) 2.5 + cos ( t ) 2.5 cos ( t ) 5 sin ( t ) , S ( t ) = 0.1 × sin ( t ) 0.1 + 0.1 × cos ( t ) 0.1 0.1 × cos ( t ) 0.1 × sin ( t )
g i ( x ) = x + 2.5 x > 0 x 2.5 x < 0 , u i ( x ) = x + 1 x > 0 x 1 x < 0 , i = 1 , , n
where k N 0 , T 0 = 1 . According to the uniqueness conditions for the Filippov solution of the time-delay system, Lemma 2 and Theorem 3, it can be seen that (23) has a unique solution in compact set Σ R 2 .
Here, define the norm with subscript σ ( t ) as | x | 0 = 0.8 | x 1 | + | x 2 | , | x | 1 = | x 1 | + | x 2 | , β 01 = 1.25 , β 10 = 1 , then ν 0 ( T ) < 0.5 , ν 1 ( T ) = 2.5 + 2 < 1 , S 1 < 0.3 , S 0 < 0.25 , therefore, let N ( t ) = j , t [ t j 1 , t j ) , γ j = max { 1 / 4 , ( 1 / 2 ) ( 1 / 2 ) j } , t j = ( T 0 / 2 ) · j , j = 1 , 2 , , such that in the continuous region of the right-hand side of system (23), it holds that
α N ( t ) + ν 0 ( D ( t ) + T ( t ) ) < 0 α N ( t ) + ν 1 ( D ( t ) + T ( t ) ) + exp γ N ( t ) t exp γ N ( t τ ( t ) ) ( t τ ( t ) ) S ( t ) 1 < 0.5 1 + 0.25 < 0
As for condition (21), after calculation, we have K = 5 , M = 2 ; thus, condition (21) holds for σ ( t ) = 0 , 1 . Moreover, for t > 10 , we have
i = 0 N ( t ) 1 [ ( γ i + 1 γ i ) t i + 1 + log β i + 1 ] + γ 0 t 0 γ N ( t ) t 1.75 + log ( 1.25 ) t 0.49 t < 0.25 t
Together with Corollary 5, the switched time-delay system (23) is exponentially incrementally asymptotically stable.With the initial state x 0 ( s ) = [ 5 × ( 0.5 + s ) 2 , 6 × ( 1 s ) 1.5 ] and y 0 ( s ) = [ 2 ( 1 + s ) 2 2 , 3 × ( 1 s ) ] for s [ 0 , 0.5 ] , the corresponding solution of system (23) are x ( t ) = ( x 1 ( t ) , x 2 ( t ) ) and y ( t ) = ( y 1 ( t ) , y 2 ( t ) ) , respectively. When t > 0.5 , σ ( t ) = 1 and time lag occurs.
Figure 3 shows the dynamical trajectories of two of the solutions with the initial function defined as x 0 ( · ) and y 0 ( · ) for system (23). Moreover, Figure 4 shows the errors between the two dynamical trajectories of their segments.

7. Conclusions

This paper researches the criteria for incremental stability for time-delay dynamical systems, including systems with continuous right-hand sides and systems with discontinuous right-hand sides, respectively. In this paper, the corresponding sufficient conditions for the exponential incremental stability of solutions for time-delay dynamical systems with continuous right-hand sides are proposed and proved. Before studying sufficient conditions for incremental stability of the systems with discontinuous right-hand sides, we first provide the conditions for the existence and uniqueness of the Filippov solution. Then, by constructing a sequence of systems with continuous right-hand sides and using the approximation method, sufficient conditions for exponential incremental stability of the systems with discontinuous right-hand sides are obtained.
There still needs to be much further work. Our theorem can be helpful for applications in other more complex scenarios and we may propose more corollaries for some more complex systems in the future. Furthermore, we may seek some other approaches to construct “continuous systems” to approximate discontinuous systems.

Author Contributions

Conceptualization, W.L.; methodology, Y.L. and W.L.; software, Y.L.; validation, Y.L. and W.L.; formal analysis, Y.L.; investigation, Y.L.; resources, W.L.; data curation, W.L.; writing—original draft preparation, Y.L.; writing—review and editing, W.L.; visualization, Y.L.; supervision, W.L.; project administration, W.L.; funding acquisition, W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is jointly supported by the National Natural Sciences Foundation of China under Grant (No. 62072111), Shanghai Municipal Science and Technology Major Project under Grant 2018SHZDZX01. This work is conducted in joint accordance with the ZJLab and the Science and Technology Commission of Shanghai Municipality (No. 19JC1420101).

Institutional Review Board Statement

This work is conducted in joint accordance with the ZJLab and the Science and Technology Commission of Shanghai Municipality (No. 19JC1420101).

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created.

Acknowledgments

We here thank Hao Qiang for his assistance with some preliminary research.

Conflicts of Interest

The authors declare no conflict of interest.

Notations

| · | χ ( t ) Vector norm with subscript χ ( t )
· χ ( t ) Matrix norm induced by | · | χ ( t )
ν χ ( t ) Matrix measure induced by | · | χ ( t )
χ ( t ) A right-continuous staircase function with respect to t, with switching points belonging to { t j }
r ( t ) A piecewise right-continuous function with respect to t, with switching points { t j }
t 0 The initial time
τ ¯ The upper bound of τ k : max k s u p t [ t 0 , ) τ k ( t ) = τ ¯
τ ̲ The lower bound of τ k : min k i n f t [ t 0 , ) τ k ( t ) = τ ̲
N ( t ) N ( t ) = # { j : t t j , j = 1 , 2 , }

References

  1. Pavlov, A.; Pogromsky, A.; Wouw, N.V.D.; Nijmeijer, H. Convergent dynamics, a tribute to Boris Pavlovich Demidovich. Syst. Control Lett. 2004, 52, 257–261. [Google Scholar] [CrossRef]
  2. Rüffer, B.S. Convergent systems vs. incremental stability. Syst. Control Lett. 2013, 62, 277–285. [Google Scholar] [CrossRef]
  3. Fromion, V.; Scorletti, G.; Ferreres, G. Nonlinear performance of a PI controlled missile: An explanation. Int. J. Robust Nonlinear Control. IFAC-Affil. J. 1999, 9, 485–518. [Google Scholar] [CrossRef]
  4. Hamadeh, A.; Stan, G.B.; Sepulchre, R.; Goncalves, J. Global State Synchronization in Networks of Cyclic Feedback Systems. IEEE Trans. Autom. Control 2012, 57, 478–483. [Google Scholar] [CrossRef]
  5. Russo, G.; Di Bernardo, M. Correction to Section V-B of “Contraction Theory and the Master Stability Function: Linking Two Approaches to Study Synchronization in Complex Networks” [Feb 09 177-181]. Circuits Syst. Express Briefs IEEE Trans. 2014, 61, 915. [Google Scholar] [CrossRef]
  6. Kim, H.; Persis, C.D. Output synchronization of Lur’e-type nonlinear systems in the presence of input disturbances. IEEE Conf. Decis. Control 2015, 4145–4150. [Google Scholar]
  7. Promion, V. Some results on the behavior of Lipschitz continuous systems. In Proceedings of the 1997 European Control Conference (ECC), Brussels, Belgium, 1–7 July 1997; pp. 2011–2016. [Google Scholar]
  8. Lohmiller, W.; Slotine, J.J.E. On contraction analysis for non-linear systems. Automatica 1998, 34, 683–696. [Google Scholar] [CrossRef]
  9. Andronov, A.A.; Vitt, A.A.; Khaikin, S.E. Vibration Theory; Fizmatgiz: Moscow, Russia, 1959. (In Russian) [Google Scholar]
  10. Barbashin, E.A. Introduction to Stability Theory; Nauka: Moscow, Russia, 1967. (In Russian) [Google Scholar]
  11. Neimark, Y.I. Points Mapping Method in the Theory of Nonlinear Vibrations; Nauka: Moscow, Russia, 1972. (In Russian) [Google Scholar]
  12. Bautin, N.N.; Leontovich, E.A. Methods and Techniques of Qualitative Study of Dynamic Systems in a Plane; Nauka: Moscow, Russia, 1976. (In Russian) [Google Scholar]
  13. Gelig, A.K.; Leonov, G.A.; Yakubovich, V.A. Stability of Nonlinear Systems with Nonunique Equilibrium State; Nauka: Moscow, Russia, 1978. (In Russian) [Google Scholar]
  14. Filippov, A. Differential equations with discontinuous right-hand side. Am. Math. Soc. Transl. 1964, 42, 199–231. [Google Scholar]
  15. Filippov, A.F. Differential Equations with Discontinuous Right-Hand Sides; Kluwer Academic: Boston, MA, USA, 1988. [Google Scholar]
  16. Volterra, V. Sur la théorie mathématique des phénom enes héréditaires. J. Math. Pures Appl. 1928, 7, 249–298. [Google Scholar]
  17. Volterra, V. Theorie Mathematique de la Lutte Pour la vie; Gauthier-Villars: Paris, France, 1931. [Google Scholar]
  18. Hale, J.K.; Lunel, S.V. Introduction to Functional Differential Equations(Applied Mathematical Sciences, 99); Springer: Berlin, Germany, 1993. [Google Scholar]
  19. Cooke, K.L. Differential—Difference equations. In International Symposium on Nonlinear Differential Equations and Nonlinear Mechanics; Academic Press: Cambridge, MA, USA, 1963; pp. 155–171. [Google Scholar]
  20. Nihtilä, M.T. Analysis and synthesis of time delay systems: H. Gorechi, S. Fuksa, P. Grabowski and A. Korytowski. Automatica 1991, 27, 1065–1067. [Google Scholar] [CrossRef]
  21. Kaminogo, T. Kneser’s property and boundary value problems for some retarded functional differential equations. Tohoku Math. J. 1978, 30, 471–486. [Google Scholar] [CrossRef]
  22. Kaminogo, T. Topological properties of solution mappings in functional differential equations. Nonlinear Anal. Theory Methods Appl. 1997, 30, 871–876. [Google Scholar] [CrossRef]
  23. Tunç, C. Qualitative properties in nonlinear Volterra integro-differential equations with delay. J. Taibah Univ. Sci. 2017, 11, 309–314. [Google Scholar] [CrossRef]
  24. Deimling, K. Multivalued Differential Equations; De Gruyter: Berlin, Germany; New York, NY, USA, 1992. [Google Scholar]
  25. Haddad, T.; Haddad, T. Delay perturbed state-dependent sweeping process. Appl. Anal. 2016, 95, 270–282. [Google Scholar] [CrossRef]
  26. Marques, M.D.P.M.; Dzonou, R. Dynamics of a particle with friction and delay. Comptes Rendus Mec. 2017, 346, 237–246. [Google Scholar] [CrossRef]
  27. Pust, L.; Pešek, L. Influence of delayed excitation on vibrations of turbine blades couple. Appl. Comput. Mech. 2013, 7, 39–52. [Google Scholar]
  28. Smirnov, G.V. Introduction to the Theory of Differential Inclusions; University of Porto: Porto, Portugal, 2002. [Google Scholar]
  29. Ali, M.; Hou, Z.; Noori, M. Stability and performance of feedback control systems with time delays. Comput. Struct. 1998, 66, 241–248. [Google Scholar] [CrossRef]
  30. Callender, A.B.; Hartree, D.R.; Porter, A.R. Time-Lag in a Control System. Philos. Trans. R. Soc. A 1936, 235, 415–444. [Google Scholar]
  31. Haddad, G. Topological properties of the sets of solutions for functional differential inclusions. Nonlinear Anal. Theory Methods Appl. 1981, 5, 1349–1366. [Google Scholar] [CrossRef]
  32. Bacciotti, A.; Ceragioli, F. Optimal Regulation and Discontinuous Stabilization. IFAC Proc. Vol. 2001, 34, 39–44. [Google Scholar] [CrossRef]
  33. Lu, W.; Chen, T. Almost periodic dynamics of a class of delayed neural networks with discontinuous activations. Neural Comput. 2008, 20, 1065–1090. [Google Scholar] [CrossRef] [PubMed]
  34. Liu, X.; Cao, J. Synchronization control of discontinuous neural networks via approximation. In Proceedings of the 2010 Chinese Control and Decision Conference, Xuzhou, China, 26–28 May 2010; pp. 782–787. [Google Scholar]
  35. Fiore, D.; Hogan, S.J.; di Bernardo, M. Contraction analysis of switched systems via regularization. Automatica 2016, 73, 279–288. [Google Scholar] [CrossRef]
  36. Lu, W.; Bernardo, M. Contraction and incremental stability of switched Carathéodory systems using multiple norms. Automatica 2016, 70, 1–8. [Google Scholar] [CrossRef]
  37. Aubin, J.P.; Bayen, A.M.; Saint-Pierre, P. Viability Theory; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  38. Qiang, H.; Lu, W. The Coupling Mechanism and Dynamic Behavior of Complex Network Systems; Fudan University: Shanghai, China, 2021. [Google Scholar]
  39. Ortiz-Robinson, N.; Rios, V. Qualitative Properties of the Solution Set for Time-Delayed Discontinuous Dynamics. Set-Valued Var. Anal. 2020, 28, 327–344. [Google Scholar] [CrossRef]
  40. Liu, B.; Lu, W.; Chen, T. New conditions on synchronization of networks of linearly coupled dynamical systems with non-Lipschitz right-hand sides. Neural Netw. 2012, 25, 5–13. [Google Scholar] [CrossRef]
  41. Ortiz-Robinson, N.; Rios, V. Forward Euler Solutions and Weakly Invariant Time-Delayed Systems. Abstr. Appl. Anal. 2012, 3, 327–344. [Google Scholar] [CrossRef]
  42. Forti, M.; Nistri, P.; Papini, D. Global exponential stability and global convergence in finite time of delayed neural networks with infinite gain. IEEE Trans. Neural Netw. 2005, 16, 1449–1463. [Google Scholar] [CrossRef]
  43. Haddad, G. Monotone viable trajectories for functional differential inclusions. J. Differ. Equ. 1981, 42, 1–24. [Google Scholar] [CrossRef]
Figure 1. Dynamical trajectories of the solutions for time-delay system (22).
Figure 1. Dynamical trajectories of the solutions for time-delay system (22).
Mathematics 11 02242 g001
Figure 2. A diagram for exponential incremental uniform stability (the error of each segment) of time-delay system (22).
Figure 2. A diagram for exponential incremental uniform stability (the error of each segment) of time-delay system (22).
Mathematics 11 02242 g002
Figure 3. Dynamical trajectories of the solutions for time-delay system (23).
Figure 3. Dynamical trajectories of the solutions for time-delay system (23).
Mathematics 11 02242 g003
Figure 4. A diagram for exponential incremental uniform stability (the error of each segment) for system (23).
Figure 4. A diagram for exponential incremental uniform stability (the error of each segment) for system (23).
Mathematics 11 02242 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lang, Y.; Lu, W. Criteria on Exponential Incremental Stability of Dynamical Systems with Time Delay. Mathematics 2023, 11, 2242. https://doi.org/10.3390/math11102242

AMA Style

Lang Y, Lu W. Criteria on Exponential Incremental Stability of Dynamical Systems with Time Delay. Mathematics. 2023; 11(10):2242. https://doi.org/10.3390/math11102242

Chicago/Turabian Style

Lang, Yingying, and Wenlian Lu. 2023. "Criteria on Exponential Incremental Stability of Dynamical Systems with Time Delay" Mathematics 11, no. 10: 2242. https://doi.org/10.3390/math11102242

APA Style

Lang, Y., & Lu, W. (2023). Criteria on Exponential Incremental Stability of Dynamical Systems with Time Delay. Mathematics, 11(10), 2242. https://doi.org/10.3390/math11102242

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