Next Article in Journal
Evaluation of Bank Innovation Efficiency with Data Envelopment Analysis: From the Perspective of Uncovering the Black Box between Input and Output
Next Article in Special Issue
A New Construction and Convergence Analysis of Non-Monotonic Iterative Methods for Solving ρ-Demicontractive Fixed Point Problems and Variational Inequalities Involving Pseudomonotone Mapping
Previous Article in Journal
Analytical Solution of the Three-Dimensional Laplace Equation in Terms of Linear Combinations of Hypergeometric Functions
Previous Article in Special Issue
A General Class of Differential Hemivariational Inequalities Systems in Reflexive Banach Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New and Improved Criteria on Fundamental Properties of Solutions of Integro—Delay Differential Equations with Constant Delay

1
Department of Mathematics, Faculty of Sciences, Van Yuzuncu Yil University, Van 65080, Turkey
2
Department of Mathematics, Zhejiang Normal University, Jinhua 321004, China
3
Department of Computer Programing, Baskale Vocational School, Van Yuzuncu Yil University, Van 65080, Turkey
4
Research Center for Interneural Computing, China Medical University Hospital, China Medical University, Taichung 404332, Taiwan
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 16 November 2021 / Revised: 10 December 2021 / Accepted: 14 December 2021 / Published: 20 December 2021
(This article belongs to the Special Issue New Advances in Functional Analysis)

Abstract

:
This paper is concerned with certain non-linear unperturbed and perturbed systems of integro-delay differential equations (IDDEs). We investigate fundamental properties of solutions such as uniformly stability (US), uniformly asymptotically stability (UAS), integrability and instability of the un-perturbed system of the IDDEs as well as the boundedness of the perturbed system of IDDEs. In this paper, five new and improved fundamental qualitative results, which have less conservative conditions, are obtained on the mentioned fundamental properties of solutions. The technique used in the proofs depends on Lyapunov-Krasovski functionals (LKFs). In particular cases, three examples and their numerical simulations are provided as numerical applications of this paper. This paper provides new, extensive and improved contributions to the theory of IDDEs.
MSC:
93D20; 26D15; 34K06; 24K20

1. Introduction

Integral equations (IEs) and integro-differential equations (IDEs) occur widely in diverse areas of applied mathematics and physics. IEs and IDEs offer powerful techniques to solve various practical problems in numerous different scientific fields, and these equations have also attracted attentions for these reasons by numerous researchers up to now. Indeed, IDEs appear when we convert initial value problems to IEs. During the last a few decades, there has been much work on fundamental properties of solutions of IDEs and IDDEs by using the LKF method, the Jensen inequality method, the fixed point method, the variation of parameters formula, Wirtinger-based integral inequality, the free-weighting matrix method and so on (see, for example, [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21]). Among these methods, the LKF method can be used very effectively during the investigation of the fundamental properties solutions of IDEs and IDDEs. The positive effects of the LKF method depend on construction or choice or definition of a suitable LKF. Indeed, the major advantage of the LKF method is that the stability in the large can be obtained without any prior knowledge of solutions. Without loss of generality, during the investigations, the LKF method requires the construction or definition of suitable functional(s) as a basic tool, known as the Lyapunov-Krasovski functional(s), to prove qualitative results on fundamental properties of solutions, stability, instability, boundedness, convergence, etc., of FDEs, IDDE, and so on. The major disadvantage of this method is that constructing or defining suitable LKFs remains an open problem in the literature up to now. However, when a suitable LKF is found, the LKF method is very effective to obtain qualitative results under weaker and effective conditions. In this paper, we will use the LKF method as a basic tool for problems to be considered here.
First, we would like to present some recent works related to this paper.
In 2020, Zeng et al., [20] studied the asymptotic stability (AS) of the delay linear system:
x ˙ ( t ) = A x ( t ) + A 1 x ( t h ) + A 2 t h t x ( s ) d s ,
where x ( t ) R n , t [ h , ) , x ( t ) = ϕ ( t ) on [ h , 0 ] , A , A 1 , A 2 R n × n and h is the constant time retardation. In [20], the new integral inequality was devised and used with the LKF to investigate the AS of the above delay linear system.
In 2020, Tian et al., [10] considered the delay linear system given by
x ˙ ( t ) = A x ( t ) + B x ( t ) + C t h t x ( s ) d s ,
where x ( t ) R n is the sate vector, t [ h , ) , x ( t ) = ϕ ( t ) on [ h , 0 ] and A , B , C R n × n . Tian et al., [10] proved the AS of the above system. In [10], the new integral inequality is obtained. Later, the delay-dependent stability criteria are established using this inequality and the LKF.
In 2020, Tian and Wang [11] investigated the AS of the following delay linear system:
x ˙ ( t ) = A x ( t ) + A d x ( t h ) + A D t h t x ( s ) d s ,
where x ( t ) R n is the sate vector, t [ h , ) , x ( t ) = ϕ ( t ) on [ h , 0 ] , A , A d , A D R n × n and h is the constant delay. In [11], firstly, the vectors-based multiple integral inequality is constructed. Then, based on this integral inequality and the LKF, new delay-dependent stability criteria are derived for the above system.
In two nice papers published recently, Liu and Zhao [7] and Zhao et al., [18] investigated (AS) of the linear system of IDDEs with state and distributed delays:
x ˙ ( t ) = A x ( t ) + A d x ( t h ) + A D t h t x ( s ) d s ,
where x ( t ) R n , A , A d , A D R n × n are constant matrices, x ( t ) = ϕ ( t ) on [ h , 0 ] , ϕ C ( [ h , 0 ] , R n ) is an initial function and h is the constant delay.
In Zhao et al., [18], Theorem 1, first the double integral inequality was defined by the authors. Then, the proposed double integral inequality was used to derive new delay-dependent stability criteria by using LMIs (linear matrix inequalities) to guarantee the AS of system (1). In [18], two numerical examples were also given to demonstrate the effectiveness and merits of the AS result.
When we look at the above delay systems all of them are linear and have simple forms. Indeed, the result of the paper [18] is very nice and interesting. To the best of knowledge, the AS result of this paper has conservative conditions because of the LKF used therein.
In this paper, particularly motivated from the result of Zhao et al., [18], Theorem 1 and the works mentioned above, we consider the following perturbed nonlinear system of IDDEs with state and distributed delays:
x ˙ ( t ) = A ( t ) x ( t ) + A d ( t ) F ( x ( t h ) ) + t h t A D ( s ) G ( x ( s ) ) d s + H ( t , x ( t ) , x ( t h ) ) ,
where x ( t ) R n , t [ h , ) , x ( t ) = ϕ ( t ) on [ h , 0 ] , h R , h > 0 is delay term, ϕ C ( [ h , 0 ] , R n ) is the initial function, A ( t ) = ( a i j ( t ) ) , A d = ( a d i j ( t ) ) , A ( t ) , A d ( t ) C [ R + , R n × n ] , R + = [ 0 , ) , A D = ( a D i j ( s ) ) , i , j = 1 , 2 , . . . , n , h s t < , A D C [ [ h , ) , R n × n ] , F , G C [ R n , R n ] , F ( 0 ) = G ( 0 ) = 0 and H C ( R + × R n × R n , R n ) .
Since the functions F, G and H are continuous, the continuity of these functions is a sufficient conditions for the existence of the solutions of the system of IDDEs (2). We also assume that the functions F, G and H satisfy the Lipschitz condition in their respective arguments x ( t h ) , x and x, x ( t h ) , respectively. This assumption is a sufficient condition for the uniqueness of the solutions of the system of IDDEs (2).
It is obvious that the system of IDDEs (2) includes all of the above systems and extends these systems from the linear case to a very general nonlinear delay system. Next, in the mentioned papers, it was focused on the maximum admissible upper bounds of the delay term. At the end of this paper, we also give some comments on the upper bounds of the delay term.
Here, in particular, we focus on the paper of Zhao et al., [18] with some details. The aim of the paper is to obtain the result of [18] under weaker conditions, to extend and improve it. We also added four new results with some numerical simulations. The information given above shortly outline the novelty, originality and contributions of this paper.
We should also mention that if we replace h by—h in the perturbed nonlinear system of IDDEs (2), then we will obtain an IDE with advanced argument. To the best of information, there is no work in the literature on the fundamental properties of solutions of IDEs with advanced arguments. This is an open problem and this work could be extended for the advanced case in the future.
The remainder of this paper is organized as follows: Section 2 includes the preliminaries such as background information of the paper and the result of Zhao et al., [18], Theorem 1. Section 3 includes the main results of this paper, the stability and integrability results, Theorems 2–4. In Section 4, an illustrative example is given to support our new results, Theorems 2–4. In Section 5 and Section 6, the instability theorem, Theorem 5 and the related illustrative example, Example 2, are provided, respectively. Similarly, in Section 7 and Section 8, the boundedness theorem, Theorem 6, and an illustrative example, Example 3, are given, respectively. Finally, Section 9 and Section 10 include the contributions and the conclusion of the paper, respectively

2. Preliminaries

We begin by couching our problem in a more general setting. Consider the system of functional differential equations (FDEs) (see, Burton [3]):
d x d t = G ( t , x t ) ,
with G : R + × C R n being continuous, where R + = [ 0 , ) , R = ( , ) , and C is the set of continuous functions ϕ : [ τ , 0 ] R n , τ > 0 . We suppose that G takes closed bounded sets of R + × C into bounded sets in R n .
As the next step, for any a 0 , t 0 0 , and x C ( [ t 0 τ , t 0 + a ] , R n ) , we assume that x t = x ( t + θ ) for θ [ τ , 0 ] .
For x R n , the norm . is denoted by
x = i = 1 n x i .
If A R n × n , then, A is defined by
A = max 1 j n i = 1 n a i j .
For any ϕ C , let
ϕ C = sup θ [ τ , 0 ] ϕ ( θ ) = ϕ ( θ ) [ τ , 0 ]
and
C H = { ϕ : ϕ C and ϕ C H < } .
Throughout this paper, we tacitly assume that the function G satisfies conditions which guarantee the uniqueness of solutions of system of FDEs (3). We should point out that IDDEs (1) and (2) are included by system of FDEs (3).
We suppose that x ( t ) = x ( t , t 0 , ϕ ) is a solution of system of FDEs (3) with x ( t ) = ϕ ( t ) on [ t 0 τ , t 0 ] .
Zhao et al., [18], Theorem 1 proved the following theorem.
Theorem 1.
(Zhao et al., [18], Theorem 1). Given h > 0 , system (1) is AS if there exist positive definite matrices P R 4 n × 4 n and Q , S , R R n × n such that the following LMI holds:
Ψ = s y m ( II 1 T P II 2 ) + e 1 T Q e 1 e 2 T Q e 2 + h 2 e 0 T S e 0 + h 2 2 e 0 T Re 0 II 3 T S II 3 3 II 4 T S II 4 5 II 5 T S II 5 2 II 6 T R II 6 4 II 7 T R II 7 6 II 8 T R II 8 < 0 ,
where
II 1 = [ e 1 T , e 3 T , e 4 T , e 5 T ] T ,
II 2 = [ e 3 T , e 1 T e 2 T , h e 1 T e 3 T , h 2 2 e 1 T e 4 T ] T ,
II 3 = e 1 e 2 ,
II 4 = e 1 + e 2 2 h e 3 ,
II 5 = e 1 e 2 + 6 h e 3 12 h 2 e 4 ,
II 6 = e 1 1 h e 3 ,
II 7 = e 1 + 2 h e 3 6 h 2 e 4 ,
II 8 = e 1 3 h e 3 + 24 h 2 e 4 60 h 3 e 5 ,
e 0 = A e 1 + A d e 2 + A D e 3 ,
and e i R n × 5 n is defined as e i = [ 0 n × ( i 1 ) n , I n , 0 n × ( 5 i ) n ] for i = 1 , 2 , . . . , 5 .
To prove this theorem, Zhao et al., [18], Theorem 1, the following LKF is defined:
V ( t ) = i = 1 4 V i ( t ) ,
where
V 1 ( t ) = η T ( t ) P η ( t ) , V 2 ( t ) = t h t x T ( s ) Q x ( s ) d s , V 3 ( t ) = h t h t u t x ˙ T ( s ) S x ˙ ( s ) d s d u , V 4 ( t ) = h t h t u t s t x ˙ T ( s ) R x ˙ ( r ) d r d s d u , η ( t ) = [ x T ( t ) , t h t x T ( s ) d s , t h t u t x T ( s ) d s d u , t h t u t s t x T ( r ) d r d s d u ] T .
Later, they obtained the derivative of this LKF along system of IDDEs (1) as the following:
V ˙ ( t ) = ξ T ( t ) Ψ ξ ( t ) ,
where
ξ ( t ) = [ x T ( t ) , x T ( t h ) , t h t x T ( s ) d s , t h t u t x T ( s ) d s d u , t h t u t s t x T ( r ) d r d s d u ] T
and the matrix Ψ is defined by (4).

3. US, UAS and Integrability of Solutions

Consider the un-perturbed system
x ˙ ( t ) = A ( t ) x ( t ) + A d ( t ) F ( x ( t h ) ) + t h t A D ( s ) G ( x ( s ) ) d s ,
which is obtained from system (2).
This section is devoted to improve and generalize the AS result of Zhao et al., [18], Theorem 1. Here, the AS result of Zhao et al., [18], Theorem 1 can be obtained under less conservative conditions and we give additional new results for system (5). The technique used in the proofs is based upon the LKF approach.
A. Assumptions
Throughout this and next sections, the following conditions are needed.
(A1)
The functions F , G C [ R n , R n ] satisfy
F ( x ) F 0 x for all x R n with F ( 0 ) = 0
and
G ( x ) G 0 x for all x R n with G ( 0 ) = 0 ,
where F 0 and G 0 are positive constants.
(A2)
There are constants F 0 , G 0 from (A1), A 0 > 0 and h > 0 such that
a i i ( t ) + j = 1 , j i n a j i ( t ) A 0 for all t R +
and
A 0 h G 0 A D ( t ) F 0 A d ( t + h ) 0 for all t [ h , ) .
(A3)
There are constants A 0 , G 0 , F 0 , h > 0 from (A2) and K 0 > 0 such that
a i i ( t ) + j = 1 , j i n a j i ( t ) A 0 for all t R +
and
A 0 h G 0 A D ( t ) F 0 A d ( t + h ) K 0 for all t [ h , ) .
The next theorem is the first main result of this paper.
Theorem 2.
We assume that conditions (A1) and (A2) are satisfied. Then, the solution x ( t ) 0 of system (5) is uniformly stable.
Proof. 
Define a new LKF II : = II ( t , x t ) by
II ( t , x t ) : = x ( t ) + t h t A d ( s + h ) F ( x ( s ) ) d s + h 0 t + η t A D ( s ) G ( x ( s ) ) d s d η ,
where h η 0 .
As the next step, it is notable that the KLF II satisfies
II ( t , 0 ) = 0 , II ( t , x t ) x .
Differentiating the LKF II in (6) along system (5), we find that
d d t II ( t , x t ) = i = 1 n x i ( t ) s g n x i ( t + 0 ) + h 0 A D ( t ) G ( x ( t ) ) A D ( t + η ) G ( x ( t + η ) ) d η + A d ( t + h ) F ( x ( t ) ) A d ( t ) F ( x ( t h ) ) = i = 1 n x i ( t ) s g n x i ( t + 0 ) + h A D ( t ) G ( x ( t ) ) t h t A D ( s ) G ( x ( s ) ) d s + A d ( t + h ) F ( x ( t ) ) A d ( t ) F ( x ( t h ) ) i = 1 n x i ( t ) s g n x i ( t + 0 ) + h G 0 A D ( t ) x ( t ) t h t A D ( s ) G ( x ( s ) ) d s + A d ( t + h ) F ( x ( t ) ) A d ( t ) F ( x ( t h ) ) .
Using conditions (A1) and (A2), we derive that
i = 1 n s g n x i ( t + 0 ) x i ( t ) i = 1 n a i i ( t ) x i ( t ) + i = 1 n j = 1 , j i n a j i ( t ) x i ( t ) + i = 1 n j = 1 n a d i j ( t ) F j ( x ( t h ) ) + t h t A D ( s ) G ( x ( s ) ) d s i = 1 n a i i ( t ) + j = 1 , j i n a j i ( t ) x i ( t ) + A d ( t ) F ( x ( t h ) ) + t h t A D ( s ) G ( x ( s ) ) d s A 0 x ( t ) + A d ( t ) F ( x ( t h ) ) + t h t A D ( s ) G ( x ( s ) ) d s .
A combination of this inequality, inequality (7) and the use of condition of (A1) imply that
d d t II ( t , x t ) A 0 x ( t ) + A d ( t ) F ( x ( t h ) ) + t h t A D ( s ) G ( x ( s ) ) d s + h G 0 A D ( t ) x ( t ) t h t A D ( s ) G ( x ( s ) ) d s + A d ( t + h ) F ( x ( t ) ) A d ( t ) F ( x ( t h ) ) A 0 x ( t ) + h G 0 A D ( t ) x ( t ) + F 0 A d ( t + h ) x ( t ) .
Hence, we obtain
d d t II ( t , x t ) A 0 h G 0 A D ( t ) F 0 A d ( t + h ) x ( t ) .
Using condition (A2), from (8) we obtain
d d t II ( t , x t ) 0 .
Thus, it is proved that the zero solution of system (5) is uniformly stable. □
The next theorem, Theorem 3, is our second main result.
Theorem 3.
We assume that conditions (A1) and (A3) are satisfied. Then, the solution x ( t ) 0 of system (5) is uniformly asymptotically stable.
Proof. 
Using condition (A1) and revising the inequality (8) according to condition (A3), we derive that
II ( t , x t ) : = x ( t ) + t h t A d ( s + h ) F ( x ( s ) ) d s + h 0 t + η t A D ( s ) G ( x ( s ) ) d s d η x ( t ) + F 0 t h t A d ( s + h ) x ( s ) d s + G 0 h 0 t + η t A D ( s ) x ( s ) d s d η x ( t ) + h F 0 sup t h s t A d ( s + h ) x ( s ) h η G 0 sup t + η s t A D ( s ) x ( s ) , h η 0 , II ( t , x t ) x ( t )
and
d d t II ( t , x t ) A 0 h G 0 A D ( t ) F 0 A d ( t + h ) x ( t ) K 0 x ( t ) < 0 , x ( t ) 0 .
From the above discussion, it follows that the LKF II ( t , x t ) has lower and upper bounds as well as its time derivative along the un-perturbed system (5) is negative definite. From the basic information of qualitative theory of IDDEs, these results are sufficient conditions for all solutions of system (5) uniformly converging to the origin (see, Burton [3], Theorem 4.4.3).
Thus, the solution x ( t ) 0 of (5) is uniformly asymptotically stable. This is the end of the proof. □
The integrability of solutions is given by Theorem 4 at the following.
Theorem 4.
We assume that conditions (A1) and (A3) are satisfied. Then, the norms of solutions of system (5) are integrable on R + in sense of Lebesgue.
Proof. 
Using the LKF II : = II ( t , x t ) , it follows that (see the proof of Theorem 2)
d d t II ( t , x t ) K 0 x ( t ) .
From this inequality, an integration gives that
II ( t , x t ) II ( t 0 , ϕ ( t 0 ) ) K 0 t 0 t x ( s ) d s .
Hence, as the next step, we derive that
t 0 x ( s ) d s K 0 1 II ( t 0 , ϕ ( t 0 ) < + .
The proof is completed. □

4. Numerical Application

Example 1.
We consider the nonlinear delay system
x 1 ( t ) x 2 ( t ) = 51 + cos t 1 1 51 + cos t x 1 ( t ) x 2 ( t ) + exp ( t ) 1 1 exp ( t ) x 1 ( t ) 1 + x 1 2 ( t 5 1 ) x 2 ( t ) 1 + x 2 2 ( t 5 1 ) + t 1 5 t 1 1 + exp ( s ) 1 2 1 1 + exp ( s ) x 1 ( t ) 1 + x 1 2 ( s ) x 2 ( t ) 1 + x 2 2 ( s ) d s ,
where t 5 1 .
Comparing systems (9) and (5), we obtain the following relations:
A ( t ) = 51 + cos t 1 1 51 + cos t = a 11 ( t ) a 12 ( t ) a 21 ( t ) a 22 ( t ) ,
A d ( t ) = exp ( t ) 1 1 exp ( t ) = a d 11 ( t ) a d 12 ( t ) a d 21 ( t ) a d 22 ( t ) ,
F ( x ( t h ) ) = F ( x 1 ( t 5 1 ) , x 2 ( t 5 1 ) ) = x 1 ( t ) 1 + x 1 2 ( t 5 1 ) x 2 ( t ) 1 + x 2 2 ( t 5 1 )
G ( x ) = G ( x 1 , x 2 ) = x 1 1 + x 1 2 x 2 1 + x 2 2 , where x = x 1 x 2 ,
A D ( t ) = 1 1 + exp ( t ) 1 2 1 1 + exp ( t ) .
Considering the above relations and proceeding some simple elementary calculations, we have the following estimates:
F ( 0 ) = 0 , F ( x ) = x 1 1 + x 1 2 ( t 5 1 ) x 2 1 + x 2 2 ( t 5 1 )
= x 1 1 + x 1 2 ( t 5 1 ) + x 2 1 + x 2 2 ( t 5 1 ) F 0 x , F 0 = 1 ,
G ( 0 ) = 0 , G ( x ) = x 1 1 + x 1 2 x 2 1 + x 2 2 = x 1 1 + x 1 2 + x 2 1 + x 2 2 G 0 x , G 0 = 1 ,
a i i ( t ) + j = 1 , j i n a j i ( t ) 49 = A 0
since
a 11 ( t ) + a 21 ( t ) = 51 + cos t + 1 49 = A 0
and
a 22 ( t ) + a 12 ( t ) = 51 + cos t + 1 49 = A 0 ,
A d ( t + h ) = exp ( t h ) 1 1 exp ( t h ) = max 1 j 2 i = 1 2 a d i j ( t + h )
= max exp ( t h ) + 1 , 1 + exp ( t h )
= 1 + exp ( t h ) , where h = 1 5 , t 5 1 ,
A D ( t ) = 1 1 + exp ( t ) 1 2 1 1 + exp ( t ) = max 1 j 2 i = 1 2 a D i j ( t )
= max 1 1 + exp ( t ) + 2 , 1 1 + exp ( t ) + 1
= 1 1 + exp ( t ) + 2 3 = D 0 .
As for the last step, from these estimates, we arrive that
A 0 h G 0 A D ( t ) F 0 A d ( t + h ) 49 5 1 × 3 1 exp ( t 5 1 ) > 46 = K 0 .
Thus, we have satisfaction of conditions (A1)–(A3), so the solution x ( t ) 0 of system (9) is uniformly stable, uniformly asymptotically stable and the norms of the solutions of the system are integrable.
In Figure 1 and Figure 2, the two dimensional system of non-linear DDEs (9) was solved by MATLAB software. The behaviors of stable solutions of system (9) are shown in Figure 1 and Figure 2.

5. Instability of Solutions

B. Assumptions
For our instability result, the following conditions are needed.
(A4)
Condition (A1) from Theorem 1 is satisfied with
A d ( t + h ) d 0 for all t R + , A D ( s ) D 0 for all s [ h , ) ,
where d 0 , D 0 , F 0 , G 0 and h are positive constants such that
h < 1 d 0 F 0 η D 0 G 0 , h η 0 .
(A5)
There are constants F 0 , G 0 , h from (A2) and M 0 > 0 , ε > 0 such that
a i i ( t ) j = 1 , j i n a j i ( t ) M 0 for all t R +
and
M 0 F 0 A d ( t + h ) h G 0 A D ( t ) ε .
The next theorem, Theorem 5, is our instability result.
Theorem 5.
We assume that conditions (A4) and (A5) are satisfied. Then, the solution x ( t ) 0 of system (5) is unstable.
Proof. 
We define a new LKF II 1 : = II 1 ( t , x t ) by
II 1 ( t , x t ) : = x ( t ) t h t A d ( s + h ) F ( x ( s ) ) d s h 0 t + η t A D ( s ) G ( x ( s ) ) d s d η ,
where h η 0 .
Using condition (A4), we obtain from (10) that
II 1 ( t , x t ) x ( t ) t h t d 0 F 0 x ( s ) d s h 0 t + η t D 0 G 0 x ( s ) d s d η x ( t ) h d 0 F 0 sup t h s t x ( s ) + h η D 0 G 0 sup t + η s t x ( s ) x ( t ) h ( d 0 F 0 η D 0 G 0 ) sup t h s t x ( s ) = 1 h ( d 0 F 0 η D 0 G 0 ) sup t h s t x ( s ) > 0
provided x ( t ) = sup t h s t x ( s ) , h < 1 d 0 F 0 η D 0 G 0 , h η 0 and x ( t ) 0 .
Next, the derivative of LKF II 1 in (10) along solutions of system (5) gives that
d d t II 1 ( t , x t ) = i = 1 n x i ( t ) s g n x i ( t + 0 ) A d ( t + h ) F ( x ( t ) ) + A d ( t ) F ( x ( t h ) ) h A D ( t ) G ( x ( t ) ) + t h t A D ( s ) G ( x ( s ) ) d s .
Using condition (A4), for the first term of (11), we obtain
i = 1 n s g n x i ( t + 0 ) x i ( t ) i = 1 n a i i ( t ) x i ( t ) i = 1 n j = 1 , j i n a j i ( t ) x i ( t ) i = 1 n j = 1 n a d i j ( t ) F j ( x ( t h ) ) t h t A D ( s ) G ( x ( s ) ) d s i = 1 n a i i ( t ) j = 1 , j i n a j i ( t ) x i ( t ) A d ( t ) F ( x ( t h ) ) t h t A D ( s ) G ( x ( s ) ) d s M 0 x ( t ) A d ( t ) F ( x ( t h ) ) t h t A D ( s ) G ( x ( s ) ) d s .
On gathering the relations (11), (12) and using conditions (A4) and (A5), we have
d d t II 1 ( t , x t ) M 0 x ( t ) A d ( t ) F ( x ( t h ) ) t h t A D ( s ) G ( x ( s ) ) d s A d ( t + h ) F ( x ( t ) ) + A d ( t ) F ( x ( t h ) ) h A D ( t ) G ( x ( t ) ) + t h t A D ( s ) G ( x ( s ) ) d s = M 0 x ( t ) A d ( t + h ) F ( x ( t ) ) h A D ( t ) G ( x ( t ) ) M 0 x ( t ) F 0 A d ( t + h ) x ( t ) h G 0 A D ( t ) x ( t ) = M 0 F 0 A d ( t + h ) h G 0 A D ( t ) x ( t ) ε x ( t ) > 0 , x ( t ) 0 .
As a result of the above study, we conclude that the solution x 1 ( t ) x 2 ( t ) = 0 0 of system (5) is unstable provided conditions (A4) and (A5) are satisfied. □

6. Numerical Application

Example 2.
Let us consider the following nonlinear delay system:
x 1 ( t ) x 2 ( t ) = 51 cos t 1 1 51 cos t x 1 ( t ) x 2 ( t ) + exp ( t ) 1 1 exp ( t ) x 1 ( t ) 1 + x 1 2 ( t 5 1 ) x 2 ( t ) 1 + x 2 2 ( t 5 1 ) + t 1 5 t 1 1 + s 2 + exp ( t ) 1 2 1 1 + s 2 + exp ( t ) x 1 ( s ) 1 + x 1 2 ( s ) x 2 ( s ) 1 + x 2 2 ( s ) d s ,
where t 5 1 .
Comparing systems (13) and (5), it follows that
A ( t ) = 51 cos t 1 1 51 cos t = a 11 ( t ) a 12 ( t ) a 21 ( t ) a 22 ( t ) .
Next, the matrices A d ( t ) , A D ( t ) and the functions G ( x 1 , x 2 ) , F ( x 1 ( t 5 1 ) , x 2 ( t 5 1 ) ) are the same as those given by Example 1. As for the following step, it is verified that
a i i ( t ) j = 1 , j i n a j i ( t ) = 51 cos t 1 > 48 = M 0
since
a 11 ( t ) a 21 ( t ) = 51 cos t 1 > 48 = M 0
and
a 22 ( t ) a 12 ( t ) = 51 cos t 1 > 48 = M 0 .
As for the next step, from the above estimates and those in Example 1, we obtain that
h = 1 5 < 1 2 + 3 5 1 d 0 F 0 η D 0 G 0 , with η = 1 5 ,
M 0 F 0 A d ( t + h ) h G 0 A D ( t )
48 2 5 1 × 3 > 45 = ε .
Therefore, the solution x 1 ( t ) x 2 ( t ) = 0 0 of system (13) is unstable.
In Figure 3 and Figure 4, the two dimensional system of non-linear DDEs (13) was solved by MATLAB software. The behaviors of unstable solutions of system (13) are shown in Figure 3 and Figure 4.

7. Boundedness

We consider the perturbed system (2) and give a new theorem on the boundedness of solutions.
For the boundedness result of this paper, we need the following additional assumption.
C. Assumption
(A6)
The function H C ( R + × R n × R n , R n ) satisfies
H ( t , x ( t ) , x ( t h ) ) H 0 ( t ) x ( t )
for all t R + and for all x , x ( t h ) R n , where H 0 C ( R + , R ) , the function H 0 is bounded and there are constants F 0 , G 0 , A 0 > 0 and h > 0 from (A2) such that
a i i ( t ) + j = 1 , j i n a j i ( t ) A 0 for all t R +
and
A 0 h G 0 A D ( t ) F 0 A d ( t + h ) H 0 ( t ) 0 for all t [ h , ) .
The following theorem, Theorem 6, is the boundedness result of this paper.
Theorem 6.
The solutions of perturbed system (2) are bounded at infinity if conditions (A1) and (A6) are satisfied
Proof. 
In this proof, we use again the LKF given by (6). Clearly, using conditions (A1) and (A6), one can find that
d d t II ( t , x t ) A 0 F 0 A d ( t + h ) h G 0 A D x ( t ) + H ( t , x ( t ) , x ( t h ) ) A 0 F 0 A d ( t + h ) h G 0 A D H 0 ( t ) x ( t ) 0 .
As a result of this inequality, we find
II ( t , x t ) II ( t 0 , ϕ ( t 0 ) ) .
Hence, it is notable from the definition of LKF (6) that
x ( t ) II ( t 0 , ϕ ( t 0 ) ) , t t 0 .
If t + , then
lim t + x ( t ) lim t + II ( t 0 , ϕ ( t 0 ) ) = II ( t 0 , ϕ ( t 0 ) ) a positive constant .
Thus, the solutions of system (2) are bounded as t + . The proof of Theorem 6 is completed. □

8. Numerical Application

Example 3.
We now take into consideration the following nonlinear delay system:
x 1 ( t ) x 2 ( t ) = 51 + cos t 1 1 51 + cos t x 1 ( t ) x 2 ( t ) + exp ( t ) 1 1 exp ( t ) x 1 ( t ) 1 + x 1 2 ( t 5 1 ) x 2 ( t ) 1 + x 2 2 ( t 5 1 ) + t 1 5 t 1 1 + exp ( s ) 1 2 1 1 + exp ( s ) x 1 ( s ) 1 + x 1 2 ( s ) x 2 ( s ) 1 + x 2 2 ( s ) d s + sin x 1 ( t ) 1 + exp ( t ) + exp ( x 1 ( t 5 1 ) ) sin x 2 ( t ) 1 + exp ( t ) + exp ( x 2 ( t 5 1 ) ) .
Comparing systems (14) and (2), it follows that
H ( t , x 1 , x 2 , x 1 ( t 5 1 ) , x 2 ( t 5 1 ) ) = sin x 1 1 + exp ( t ) + exp ( x 1 ( t 5 1 ) ) sin x 2 1 + exp ( t ) + exp ( x 2 ( t 5 1 ) ) , H ( t , x 1 , x 2 , x 1 ( t 5 1 ) , x 2 ( t 5 1 ) ) sin x 1 1 + exp ( t ) + exp ( x 1 ( t 5 1 ) ) + sin x 2 1 + exp ( t ) + exp ( x 2 ( t 5 1 ) ) sin x 1 + sin x 2 1 + exp ( t ) x 1 + x 2 1 + exp ( t ) = 1 1 + exp ( t ) x , where H 0 ( t ) = 1 1 + exp ( t ) .
Next, the matrices A ( t ) , A d ( t ) , A D ( t ) and the functions G ( x 1 , x 2 ) , F ( x 1 ( t 5 1 ) , x 2 ( t 5 1 ) ) are the same as those in Example 1. Hence, it is seen that
A 0 h G 0 A D ( t ) F 0 A d ( t + h ) H 0 ( t ) 49 5 1 × 3 2 1 1 + exp ( t ) > 0 .
Then, it is obvious that all solutions x 1 ( t ) x 2 ( t ) of system (14) are bounded.
In Figure 5 and Figure 6, the two dimensional system of non-linear DDEs (14) was solved by MATLAB software. The behaviors of boundedness solutions of system (14) are shown in Figure 5 and Figure 6.
Remark 1.
In the three previous sections, Section 4, Section 6 and Section 8, called “Numerical application”, to illustrate Theorems 2–5 and to show their applications in particular cases, Examples 1–3 were solved using the 4th order Runge-Kutta method via MATLAB Simulink software by taking the step size as 0.001.

9. Contributions

Here, we would like to summarize the contributions and novelty of this paper and to explain comparisons with some former literature.
(1)
The system of IDDEs considered in Zhao et al., [18], system (1), is linear. System (1) also has a simple form. However, perturbed and unperturbed systems of IDDEs (2) and (5), respectively, are non-linear and they have more general forms than system (1). Indeed, instead of the constant matrices A , A d , A D and linear term x ( t h ) of (1), we took the variable matrices A ( t ) , A d ( t ) , A D ( s ) and the non-linear functions F ( x ( t h ) ) , G ( x ) in both of systems (2) and (5). We also added the nonlinear perturbed term H ( t , x ( t ) , x ( t h ) ) to system (2). These mean that system (2) and (5) extend and improve system (1) from linear case to the nonlinear case.
(2)
Zhao et al., [18], Theorem 1 proved a nice result on the asymptotic stability of system (1) (see the above Theorem 1) by using a very interesting approach. Indeed, the authors used the new double integral inequality in terms of LMIs and the LKF (see Section 2) to prove [18], Theorem 1. However, the use of double integral inequality, LMIs and the LKF can cause stronger conditions for the asymptotic stability result of [18], Theorem 1. This case can be seen when we look at the conditions of [18], Theorem 1. In this paper, we did not use the double integral inequality, LMIs and the LKF of [18], Theorem 1. From this point of view, we defined two new, different, interesting and effective LKFs (see the LKFs given by (6) and (12)) and used these two LKFs. Here, the LKFs (6), (12) can lead very less conservative, more suitable and practicable conditions. This case can be checked clearly when we compare the LKF used by Zhao et al., [18], Theorem 1 (see, Section 2) and the LKF (6) used to prove the uniform asymptotic stability result of this paper, Theorem 3.
Indeed, to the best of information, up to now, the LKFs
II ( t , x t ) : = x ( t ) + t h t A d ( s + h ) F ( x ( s ) ) d s + h 0 t + η t A D ( s ) G ( x ( s ) ) d s d η .
and
II 1 ( t , x t ) : = x ( t ) t h t A d ( s + h ) F ( x ( s ) ) d s h 0 t + η t A D ( s ) G ( x ( s ) ) d s d η .
were not used in the literature to discuss the fundamental behaviors of solutions related to IDDEs and DDEs. Next, indeed, as it is well-known, the investigation of fundamental behaviors of trajectories of non-linear systems of IDDEs and DDEs is a difficult task. However, the effective roles of the LKFs (6) and (12) come from the integral terms
t h t A d ( s + h ) F ( x ( s ) ) d s , h 0 t + η t A D ( s ) G ( x ( s ) ) d s d η
and
t h t A d ( s + h ) F ( x ( s ) ) d s and h 0 t + η t A D ( s ) G ( x ( s ) ) d s d η .
Without loss of generality, if suitable LKFs could be found, then very less conservative qualitative criteria can be obtained.
(3)
The UAS result of this paper, Theorem 3, extends and improves the result of Zhao et al., [18], Theorem 1 under weaker and more suitable conditions. Next, indeed, UAS implies US, but its converse is not true. Additionally, to the best of information, there is no result on the instability of system of IDDEs (5). Theorem 5 is the first instability result in the literature. The uniform stability and integrability results of unperturbed system (5) and boundedness result of perturbed system (2), i.e., Theorems 2, 4 and 6, respectively, are new and original. These theorems have also very suitable and less conservative conditions for proper and related applications in various fields of sciences and engineering.
Next, for the cases of US and UAS, since
A 0 F 0 A d ( t + h ) > 0 ( see ( A 2 ) ) ,
the upper bound for the constant delay term h is determined as the following, respectively:
h A 0 F 0 A d ( t + h ) G 0 A D ( t ) and h < A 0 F 0 A d ( t + h ) G 0 A D ( t ) , A D ( t ) 0 .
As for the case of the instability, since
M 0 F 0 A d ( t + h ) h G 0 A D ( t ) > 0 ( see ( A 5 ) ) ,
h < 1 d 0 F 0 η D 0 G 0 , h η 0 , ( see ( A 4 ) ) ,
the upper bound for the constant delay term h is determined as the following:
h = min { M 0 F 0 A d ( t + h ) G 0 A D ( t ) , 1 d 0 F 0 η D 0 G 0 } , h η 0 .
Finally, for the case of boundedness, since
A 0 F 0 A d ( t + h ) H 0 ( t ) > 0 ,
it follows that
h A 0 F 0 A d ( t + h ) H 0 ( t ) G 0 A D ( t ) , A D ( t ) 0 ,
where the function H 0 ( t ) is bounded such that
H 0 ( t ) B 0 with A 0 F 0 A d ( t + h ) B 0 > 0 .
(4)
In this paper, three examples, Examples 1–3, are provided, which satisfy the conditions of the related theorems. These examples are solved by the 4th order Runge-Kutta method in MATLAB software.

10. Conclusions

In this article, we have discussed the US, UAS, integrability, instability and boundedness problems of un-perturbed and perturbed systems of non-linear IDDEs, respectively. Five new theorems have been obtained and the obtained theorems are new and original with scientific novelties. The proofs of these theorems are performed by defining two new LKFs. Three examples are given as numerical applications of our theorems. We believe that the results of this paper provide new and effective contributions to the literature and may be useful to researchers working on these concepts.

Author Contributions

Conceptualization, O.T., C.T. and J.-C.Y.; Data curation, O.T., C.T.; Formal analysis, O.T., C.T., Y.W. and J.-C.Y.; Funding acquisition, Y.W.; Methodology, C.T., O.T. and Y.W.; Project administration, J.-C.Y.; Supervision, Y.W. and J.-C.Y.; Validation, O.T., C.T.; Visualization, O.T.; Writing—original draft, O.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by financial grants from the National Natural Science Foundation of China (No.: 12171435).

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. Bohner, M.; Tunç, O.; Tunç, C. Qualitative analysis of Caputo fractional integro-differential equations with constant delays. Comp. Appl. Math. 2021, 40, 214. [Google Scholar] [CrossRef]
  2. Berezansky, L.; Diblík, J.; Svoboda, Z.; Šmarda, Z. Uniform exponential stability of linear delayed integro-differential vector equations. J. Differ. Equ. 2021, 270, 573–595. [Google Scholar] [CrossRef]
  3. Burton, T.A. Volterra Integral and Differential Equations, 2nd ed.; Mathematics in Science and Engineering, 202; Elsevier: Amsterdam, The Netherlands, 2005. [Google Scholar]
  4. Graef, J.R.; Tunç, C. Continuability and boundedness of multi-delay functional integro-differential equations of the second order. RACSAM 2015, 109, 169–173. [Google Scholar] [CrossRef]
  5. Hale, J.K.; Verduyn Lunel, S.M. Introduction to Functional-Differential Equations; Applied Mathematical Sciences, 99; Springer: New York, NY, USA, 1993. [Google Scholar] [CrossRef]
  6. Huyen, D.T.K.; Yao, J.-C. Solution stability of a linearly perturbed constraint system and applications. Set-Valued Var. Anal. 2019, 27, 169–189. [Google Scholar] [CrossRef]
  7. Liu, X.; Zhao, D. New stability criterion for time-delay systems via an augmented Lyapunov-Krasovskii functional. Appl. Math. Lett. 2021, 116, 7. [Google Scholar] [CrossRef]
  8. Nieto, J.J.; Tunç, O. An application of Lyapunov-Razumikhin method to behaviors of Volterra integro-differential equations. RACSAM 2021, 115, 197. [Google Scholar] [CrossRef]
  9. Raffoul, Y. Exponential stability and instability in finite delay nonlinear Volterra integro-differential equations. Dyn. Contin. Discrete Impuls. Syst. Ser. A Math. Anal. 2013, 1, 95–106. [Google Scholar]
  10. Tian, J.; Ren, Z.; Zhong, S. A new integral inequality and application to stability of time-delay systems. Appl. Math. Lett. 2020, 101, 7. [Google Scholar] [CrossRef]
  11. Tian, Y.; Wang, Z. A new multiple integral inequality and its application to stability analysis of time-delay systems. Appl. Math. Lett. 2020, 105, 8. [Google Scholar] [CrossRef]
  12. Tunç, C.; Tunç, O. New qualitative criteria for solutions of Volterra integro-differential equations. Arab. J. Basic Appl. Sci. 2018, 3, 158–165. [Google Scholar] [CrossRef]
  13. Tunç, C.; Tunç, O. On the stability, integrability and boundedness analyses of systems of integro-differential equations with time-delay retardation. RACSAM 2021, 115, 1–17. [Google Scholar] [CrossRef]
  14. Tunç, O. On the behaviors of solutions of systems of non-linear differential equations with multiple constant delays. RACSAM 2021, 115, 164. [Google Scholar] [CrossRef]
  15. Tunç, O.; Atan, Ö.; Tunç, C.; Yao, J.-C. Qualitative analyses of integro-fractional differential equations with Caputo derivatives and retardations via the Lyapunov-Razumikhin Method. Axioms 2021, 10, 58. [Google Scholar] [CrossRef]
  16. Xu, X.; Liu, L.; Feng, G. Stability and stabilization of infinite delay systems: A Lyapunov-based approach. IEEE Trans. Automat. Control 2020, 11, 4509–4524. [Google Scholar] [CrossRef]
  17. Wang, Q. The stability of a class of functional differential equations with infinite delays. Ann. Differ. Equ. 2000, 1, 89–97. [Google Scholar]
  18. Zhao, N.; Lin, C.; Chen, B.; Wang, Q.G. A new double integral inequality and application to stability test for time-delay systems. Appl. Math. Lett. 2017, 65, 26–31. [Google Scholar] [CrossRef]
  19. Zhao, J.; Meng, F. Stability analysis of solutions for a kind of integro-differential equations with a delay. Math. Probl. Eng. 2018, 2018, 9519020. [Google Scholar] [CrossRef] [Green Version]
  20. Zeng, H.B.; He, Y.; Wu, M.; Jinhua, S. New results on stability analysis for systems with discrete distributed delay. Autom. J. IFAC 2015, 60, 189–192. [Google Scholar] [CrossRef]
  21. Zhu, X.L.; Yang, G.H. New stability criteria for systems with time-varying delay. In Proceedings of the 2008 Chinese Control and Decision Conference, Yantai, China, 2–4 July 2008. [Google Scholar] [CrossRef]
Figure 1. The trajectories of x 1 ( t ) of system (9) and their behaviors.
Figure 1. The trajectories of x 1 ( t ) of system (9) and their behaviors.
Mathematics 09 03317 g001
Figure 2. The trajectories of x 2 ( t ) of system (9) and their behaviors.
Figure 2. The trajectories of x 2 ( t ) of system (9) and their behaviors.
Mathematics 09 03317 g002
Figure 3. The trajectories of x 1 ( t ) of system (13) and their behaviors.
Figure 3. The trajectories of x 1 ( t ) of system (13) and their behaviors.
Mathematics 09 03317 g003
Figure 4. The trajectories of x 2 ( t ) of system (13) and their behaviors.
Figure 4. The trajectories of x 2 ( t ) of system (13) and their behaviors.
Mathematics 09 03317 g004
Figure 5. The trajectories of x 1 ( t ) of system (14) and their behaviors.
Figure 5. The trajectories of x 1 ( t ) of system (14) and their behaviors.
Mathematics 09 03317 g005
Figure 6. The trajectories of x 2 ( t ) of system (14) and their behaviors.
Figure 6. The trajectories of x 2 ( t ) of system (14) and their behaviors.
Mathematics 09 03317 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tunç, C.; Wang, Y.; Tunç, O.; Yao, J.-C. New and Improved Criteria on Fundamental Properties of Solutions of Integro—Delay Differential Equations with Constant Delay. Mathematics 2021, 9, 3317. https://0-doi-org.brum.beds.ac.uk/10.3390/math9243317

AMA Style

Tunç C, Wang Y, Tunç O, Yao J-C. New and Improved Criteria on Fundamental Properties of Solutions of Integro—Delay Differential Equations with Constant Delay. Mathematics. 2021; 9(24):3317. https://0-doi-org.brum.beds.ac.uk/10.3390/math9243317

Chicago/Turabian Style

Tunç, Cemil, Yuanheng Wang, Osman Tunç, and Jen-Chih Yao. 2021. "New and Improved Criteria on Fundamental Properties of Solutions of Integro—Delay Differential Equations with Constant Delay" Mathematics 9, no. 24: 3317. https://0-doi-org.brum.beds.ac.uk/10.3390/math9243317

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