Next Article in Journal
Collaborative Learning Based on Harry Potter for Learning Geometric Figures in the Subject of Mathematics
Next Article in Special Issue
Some Fractional Dynamic Inequalities of Hardy’s Type via Conformable Calculus
Previous Article in Journal
Secondary Education Students’ Beliefs about Mathematics and Their Repercussions on Motivation
Previous Article in Special Issue
Nonlinear Integro-Differential Equations Involving Mixed Right and Left Fractional Derivatives and Integrals with Nonlocal Boundary Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integral Representation for the Solutions of Autonomous Linear Neutral Fractional Systems with Distributed Delay

1
Faculty of Mathematics and Informatics, University of Plovdiv, 4000 Plovdiv, Bulgaria
2
Department of Mathematics, University of Architecture, Civil Engineering and Geodesy, 1046 Sofia, Bulgaria
3
Department of Matematics and Physics, University Of Food Technology, 4000 Plovdiv, Bulgaria
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(3), 364; https://doi.org/10.3390/math8030364
Submission received: 1 February 2020 / Revised: 2 March 2020 / Accepted: 2 March 2020 / Published: 6 March 2020

Abstract

:
The aim of this work is to obtain an integral representation formula for the solutions of initial value problems for autonomous linear fractional neutral systems with Caputo type derivatives and distributed delays. The results obtained improve and extend the corresponding results in the particular case of fractional systems with constant delays and will be a useful tool for studying different kinds of stability properties. The proposed results coincide with the corresponding ones for first order neutral linear differential systems with integer order derivatives.

1. Introduction and Notations

Fractional Calculus has a long history, but it has attracted considerable attention recently as an important tool for modeling of various real problems, such as viscoelastic systems, diffusion processes, signal and control processing, and seismic processes. Detailed information about the fractional calculus theory and its applications can be found in the monographs [1,2,3,4]. Some results for fractional linear systems with delays are in given in the book [5]. The monograph [6] is devoted to the impulsive differential and functional differential equations with fractional derivatives, as well as to some of their applications.
It is well known that the study of linear fractional equations (integral representation, several types of stability, etc.) is an evergreen theme for research. Concerning these fields of fundamental and qualitative investigations for linear fractional ordinary differential equations and systems we refer to [2,4,7] and the references therein. Using the Laplace transform method, several interesting results in this direction are obtained in [8,9] as well. Regarding works concerning fractional differential systems with constant delays, we point out [10,11,12,13]. Concerning the retarded differential systems with variable or distributed delays—fundamental theory and application (stability properties)—we refer to [11,14,15,16,17,18]. Neutral fractional systems with distributed delays are essentially studied less (see [19,20,21]). Stability properties of retarded fractional systems with derivatives of distributed order are studied in [22]. One of the existing best applications of fractional order equations with delays is modeling human manual control, in which perceptual and neuromuscular delays introduce a delay term. As interesting studies, we refer to [23,24].
The problem of establishing an integral representation for the solutions for neutral or delayed linear fractional differential equations and/or systems needs a theorem for the existence of a fundamental matrix, i.e., theorem for existence and uniqueness of the solution to the initial value problem (IVP) in the case of discontinuous initial functions. As far as we know, there are only a few results concerning the IVP for delayed and neutral systems with discontinuous initial function, for the delayed case [14,15,25,26,27] and for the neutral case [28].
The aim of the work is to prove an integral representation formula for the general solution of an autonomous linear fractional neutral system with Caputo type derivatives and distributed delays. Note that our results extend and improve the results obtained in [10,12,15]. The proposed results coincide with the corresponding ones for a first order neutral linear differential system with integer order derivatives.
The paper is organized as follows. In Section 2, we recall some necessary definitions of Riemann-Liouville and Caputo fractional derivatives, as well as part of their properties. In this section, we also present the linear neutral fractional system under consideration together with some conditions. In Section 3, as a main result, integral representations of the solutions of the IVP for autonomous linear fractional neutral system with Caputo type derivatives and distributed delays are obtained for the homogeneous and inhomogeneous case. In Section 4, we present an illustrative example. In Section 5 we explain the practical benefits and application options of the obtained theoretical results.
In what follows, we use the notations: N , R and C – the sets of natural, real and complex numbers, respectively; m , n – the set of integers m , m + 1 , , n ( m n ) ; R n × n – the space of real n × n matrices A with elements A p q ; R n = R n × 1 ; A – the transposed matrix A with elements ( A ) p q = A q p . The elements of R n are the real column n-vectors x = [ x 1 ; x 2 ; ; x n ] with elements x k . The row n-vectors are denoted as ξ = [ ξ 1 , ξ 2 , , ξ n ] (note that the elements of a vector column and a vector row are separated by “;” and “,”, respectively). The identity and the zero matrices are denoted by E and Θ , respectively.
We also denote C + = { p C | Re ( p ) > 0 } , C ¯ + = { p C | Re p 0 } , C = C C ¯ + , R + = ( 0 , ) , and J s = [ s , ) . For p C , y = [ y 1 ; y 2 ; ; y n ] C n and β = ( β 1 , β 2 , , β n ) , β k [ 1 , 1 ] we set I β ( p ) = diag p β 1 , p β 2 , , p β n and I β ( y ) = diag y 1 β 1 , y 2 β 2 , , y n β n . The linear space of locally Lebesgue integrable functions f : R R is denoted by L 1 loc ( R , R ) .

2. Preliminaries and Problem Statement

Below, the definitions of Riemann–Liouville and Caputo fractional derivatives and some of their properties necessary for our exposition are described in order to avoid possible misunderstandings. For more details and other properties, we refer to [2,3,4].
Let α ( 0 , 1 ) be an arbitrary number. Then for a R , each t > a and f L 1 loc ( R , R ) the left-sided fractional integral operator, the left side Riemann-Liouville and Caputo fractional derivatives of order α are defined by
( D a + α f ) ( t ) = 1 Γ ( α ) a t ( t s ) α 1 f ( s ) d s , R L D a + α f ( t ) = d d t ( D a + ( 1 α ) f ( t ) ) , C D a + α f ( t ) = R L D a + α [ f ( s ) f ( a ) ] ( t ) = R L D a + α f ( t ) f ( a ) Γ ( 1 α ) ( t a ) α ,
respectively. The following relations [4] involving fractional derivatives will be used
( D a + 0 f ) ( t ) = f ( t ) , C D a + α D a + α f ( t ) = f ( t ) , D a + α C D a + α f ( t ) = f ( t ) f ( a ) .
Concerning the Laplace transform L ,
L f ( p ) = 0 exp ( p t ) f ( t ) d t , p C ,
we shall need the relations
L D 0 + α f ( p ) = p α ( L f ) ( p ) , L R L D 0 + α f ( p ) = p α ( L f ) ( p ) [ R L D 0 + α 1 f ( t ) ] t = 0 , L C D 0 + α f ( p ) = p α ( L f ) ( p ) p α 1 f ( 0 ) .
In what follows, we consider the autonomous linear neutral fractional system with distributed delay
D α ( X ( t ) l = 1 r τ 0 d θ V l ( θ ) X ( t + θ ) ) = i = 0 m σ 0 d θ U i ( θ ) X ( t + θ ) + F ( t ) ,
as well as the corresponding homogeneous system
D α ( X ( t ) l = 1 r τ 0 d θ V l ( θ ) X ( t + θ ) ) = i = 0 m σ 0 d θ U i ( θ ) X ( t + θ ) ,
where
X , F : J 0 R n , U i , V l : R R n × n , U i ( θ ) = u k j i ( θ ) , V l ( θ ) = v k j l ( θ ) , τ , σ > 0 , τ r ( 0 , τ ] , l 1 , r , σ i ( 0 , σ ] , i 1 , m , h = max ( σ , τ ) , σ 0 = 0 , α = ( α 1 , α 2 , , α n ) , α k ( 0 , 1 ) , k 1 , n , J s = [ s , ) .
For simplicity, D α k denotes the left side Caputo fractional derivative C D 0 + α k in (1) and (2), and we use the notations
D α X ( t ) = D α 1 x 1 ( t ) ; D α 2 x 2 ( t ) ; ; D α n x n ( t ) , D α = diag ( D α 1 , D α 2 , , D α n ) , X ( t ) = [ x 1 ( t ) ; x 2 ( t ) ; ; x n ( t ) ] , F ( t ) = [ f 1 ( t ) ; f 2 ( t ) ; ; f n ( t ) ] .
Denote by B V [ h , 0 ] the linear space of matrix valued functions
W : R R n × n , W ( θ ) = ω k j ( θ )
with bounded variation in θ on [ h , 0 ] ,
Var [ h , 0 ] W ( . ) = k , j = 1 n Var [ h , 0 ] w k j ( . ) , | W ( θ ) | = k , j = 1 n | w k j ( θ ) | .
As a space of initial functions, we use the Banach space C ˜ = P C ( [ h , 0 ] , R n ) of the piecewise continuous on [ h , 0 ] vector functions Φ = [ ϕ 1 ; ϕ 2 ; ; ϕ n ] : [ h , 0 ] R n with norm
Φ = k = 1 n sup s [ h , 0 ] | ϕ k ( s ) | < .
The initial condition for the system (1) or (2) is
X ( t ) = Φ ( t ) , t [ h , 0 ] .
Definition 1.
The vector function X is a solution of the IVP (1), (3) in the interval J h if X | J 0 C ( J 0 , R n ) and if it satisfies the system (1) for t R + and the initial condition (3) for t [ h , 0 ] .
We say that for the kernels U i : R R n × n , V l : R R n × n the assumptions (SA) are fulfilled, if for each i 0 , m and l 1 , r the following conditions hold.
(SA1) The matrix valued functions θ U i ( θ ) and θ V l ( θ ) are measurable in θ R and normalized so that U i ( θ ) = 0 and V l ( θ ) = 0 for θ 0 , U i ( θ ) = U i ( σ i ) for θ σ i and V l ( θ ) = V l ( τ l ) for θ τ l .
(SA2) The kernels U i ( θ ) and V l ( θ ) are left continuous for θ ( σ , 0 ) and θ ( τ , 0 ] and U i ( · ) , V l ( · ) B V [ h , 0 ] .
(SA3) The Lebesgue decomposition of the kernels U i ( θ ) and V l ( θ ) for θ [ h , 0 ] is
U i ( θ ) = i ( θ ) + h θ B i ( s ) d s + Υ i ( θ ) , V l ( θ ) = ˜ l ( θ ) + h θ B ˜ l ( s ) d s + Υ ˜ l ( θ ) ,
where A i = a k j i , A ˜ l = a ˜ k j l R n × n and
i ( θ ) = a k j i H ( θ + σ i ) , ˜ l ( θ ) = a ˜ k j l H ( θ + τ l ) , Υ i ( θ ) = g k j i ( θ ) , Υ ˜ l ( θ ) = g ˜ k j l ( θ ) C ( R , R n × n ) , B i ( θ ) = b k j i ( θ ) , B ˜ l ( θ ) = b ˜ k j l ( θ ) L 1 loc ( R n × n , R n × n ) .
Remark 1.
The conditions (SA) are used essentially in the work [21] to establish an apriory estimate of all solutions of the IP (1), (3), which estimate guaranties that the Laplace transform can be correct applied to System (2) and to System (1) too, when the function F is exponentially bounded.
Let s 0 be an arbitrary number, J s = [ s , ) and consider the matrix IVP
D α Q ( t , s ) l = 1 r τ 0 d θ V l ( θ ) Q ( t + θ , s ) = i = 0 m σ 0 d θ U i ( t , θ ) Q ( t + θ , s )
with initial condition
Q ( t , t ) = I ; Q ( t , s ) = 0 , t < s .
Definition 2.
For each s 0 the matrix valued function
t Q ( t , s ) = γ k j ( t , s ) , Q ( · , s ) : J s R n × n ,
is called a solution of the IVP (4), (5) for t J s , if Q ( · , s ) is continuous in t on J s and satisfies the matrix Equation (4) for t ( s , ) and the initial condition (5).
It is well known that the problem of existence of a fundamental matrix for a linear homogeneous fractional system (delayed or neutral) leads to establishing that the corresponding IVP (4), (5) with discontinuous initial function has a unique solution. In the case when s = 0 , the matrix Q ( t ) = Q ( t , 0 ) will be called fundamental (or Cauchy) matrix of system (2).
Following [20,21], we introduce the characteristic matrix of System (2)
G ( p ) = I α ( p ) W ( p ) ,
where
W ( p ) = i = 0 m U i ( p ) + I α ( p ) l = 1 r V l ( p ) , i 0 , m , l 1 , r , U i ( p ) = h 0 exp ( p θ ) d u k j i ( θ ) , V l ( p ) = h 0 exp ( p θ ) d v k j l ( θ ) .

3. Main Results

The results in this section are a generalization of the results concerning the autonomous case obtained in [10,15,16,25].
Theorem 1.
Let us assume the conditions (SA) are satisfied. Then the IVP (4), (5) has a unique solution Q ( t , s ) in J s for every s 0 and the fundamental matrix Q ( t , 0 ) = Q ( t ) of Equation (2) is
Q ( t ) = L 1 I α 1 ( p ) G 1 ( p ) ( t ) .
Proof. 
Using the results from [28], we obtain that the IVP (4), (5) has a unique solution Q ( t , s ) in J s for every s 0 , and hence, a fundamental matrix Q ( t , 0 ) = Q ( t ) . In virtue of Theorem 3 [21], we can conclude that the Laplace transform can be applied to both sides of Equation (4). Substituting t + θ = η we obtain
0 exp ( p t ) i = 0 m h 0 d θ U i ( θ ) Q ( t + θ ) d t = i = 0 m h 0 d θ U i ( θ ) exp ( p θ ) θ 0 exp ( p η ) Q ( η ) d η + 0 exp ( p η ) Q ( η ) d η i = 0 m h 0 exp ( p θ ) d θ U i ( θ ) = L Q ( p ) i = 0 m h 0 exp ( p θ ) d θ U i ( θ ) .
In a similar way for the left-hand side of Equation (4), we have that
L D α Q ( t ) l = 1 r τ 0 d θ V l ( θ ) Q ( t + θ ) ( p ) = I α ( p ) L Q ( p ) E l = 1 r h 0 exp ( p θ ) d θ V l ( θ ) I α 1 ( p ) E l = 1 r τ 0 d θ V l ( θ ) Q ( θ ) .
From Equation (8), it follows that
L Q ( p ) I α I α ( p ) l = 1 r h 0 exp ( p θ ) d θ V l ( θ ) l = 1 m h 0 exp ( p θ ) d θ U i ( θ ) = I α 1 ( p )
and hence, L Q ( p ) = I α 1 ( p ) G 1 ( p ) , which completes the proof. □
Let us introduce the following functions:
Φ l ( t ) = Φ ( t ) , t [ τ l , 0 ] , Φ l ( t ) = 0 , t R [ τ l , 0 ] , l 1 , r ,
and
Φ i ( t ) = Φ ( t ) , t [ σ i , 0 ] , Φ i ( t ) = 0 , t R [ σ i , 0 ] , i 1 , m .
Then, applying the Laplace transform second shifting theorem, we obtain
θ 0 exp ( p ( η θ ) ) Φ l ( η ) d η = exp ( p θ ) L Φ l ( t ) ( p ) , θ 0 exp ( p ( η θ ) ) Φ i ( η ) d η = exp ( p θ ) L Φ i ( t ) ( p ) .
Now we are in position to prove the following theorem.
Theorem 2.
Let us assume the conditions (SA) are satisfied. Then for each Φ C ˜ the IVP (2), (3) has a unique solution X Φ ( t ) with the integral representation:
X Φ ( t ) = Q ( t ) Φ ( 0 ) l = 1 r h 0 d θ V l ( θ ) Φ l ( θ ) + l = 1 r h 0 d θ V l ( θ ) Φ l ( t + θ ) + l = 1 r h 0 d θ V l ( θ ) D 1 2 Q ( t ) D 1 2 Φ l ( t + θ ) + i = 0 m h 0 [ d θ U i ( θ ) ] D 1 α Q ( t ) Φ i ( t + θ ) + i = 0 m h 0 [ d θ U i ( θ ) ] D α Φ i ( t + θ ) .
Proof. 
Let Φ C ˜ . Then using the results from [26], we can conclude that the IVP (2), (3) has a unique solution X Φ ( t ) . In virtue of Theorem 3 from [21], we can conclude that the Laplace transform can be applied to both sides of Equation (2). Then, substituting X Φ ( t ) in Equation (2), applying the Laplace transform to Equation (2) and substituting t + θ = η , we obtain for the right-hand side of Equation (2)
L i = 0 m h 0 [ d θ U i ( θ ) ] X Φ ( t + θ ) ( p ) = 0 exp ( p η ) X Φ ( η ) d η i = 0 m h 0 exp ( p θ ) d θ U i ( θ ) + i = 0 m h 0 [ d θ U i ( θ ) ] ( θ 0 exp ( p ( θ η ) ) X Φ ( η ) d η ) = L X Φ ( t ) ( p ) i = 0 m h 0 exp ( p θ ) d θ U i ( θ ) + i = 0 m h 0 [ d θ U i ( θ ) ] θ 0 exp ( p ( θ η ) ) Φ i ( η ) d η .
Similarly, for the left-hand side of Equation (2), one obtains that
L D α X Φ ( t ) l = 1 r h 0 d θ V l ( θ ) X Φ ( t + θ ) ( p ) = I α 1 ( p ) Φ ( 0 ) l = 1 r h 0 d θ V l ( θ ) Φ l ( θ ) + L X Φ ( t ) ( p ) I α ( p ) I α ( p ) l = 1 r h 0 exp ( p θ ) d θ V l ( θ ) I α ( p ) l = 1 r h 0 [ d θ V l ( θ ) ] θ 0 exp ( p ( θ η ) ) Φ l ( η ) d η .
From Equations (11) and (12), it follows
L X Φ ( p ) I α ( p ) I α ( p ) l = 1 r h 0 exp ( p θ ) d θ V l ( θ ) i = 0 m h 0 exp ( p θ ) d θ U i ( θ ) = I α 1 ( p ) Φ ( 0 ) l = 1 r h 0 d θ V l ( θ ) Φ l ( θ ) + I α ( p ) l = 1 r h 0 d θ V l ( θ ) θ 0 exp ( p ( θ η ) ) Φ l ( η ) d η + i = 0 m h 0 d θ U i ( θ ) θ 0 exp ( p ( θ η ) ) Φ i ( η ) d η
and hence,
L X Φ ( p ) = G 1 ( p ) I α 1 ( p ) Φ ( 0 ) l = 1 r h 0 d θ V l ( θ ) Φ l ( θ ) + G 1 ( p ) I α ( p ) + l = 1 r h 0 d θ V l ( θ ) θ 0 exp ( p ( θ η ) ) Φ l ( η ) d η + G 1 ( p ) i = 0 m h 0 d θ U i ( θ ) θ 0 exp ( p ( θ η ) ) Φ i ( η ) d η .
The representations of Equations (7) and (13) imply that
L X Φ ( p ) = L Q ( t ) ( p ) Φ ( 0 ) l = 1 r h 0 [ d θ V l ( θ ) ] Φ l ( θ ) + I 1 ( p ) L Q ( t ) ( p ) l = 1 r h 0 d θ V l ( θ ) θ 0 exp ( p ( θ η ) ) Φ l ( η ) d η + I 1 α ( p ) L Q ( t ) ( p ) i = 0 m h 0 d θ U i ( θ ) θ 0 exp ( p ( θ η ) ) Φ i ( η ) d η .
In view of Equation (9), we obtain for the second term in the right-hand side of Equation (14) that
I 1 ( p ) L Q ( t ) ( p ) l = 1 r h 0 d θ V l ( θ ) θ 0 exp ( p ( θ η ) ) Φ l ( η ) d η = I 1 2 ( p ) L Q ( t ) ( p ) I 1 2 ( p ) l = 1 r h 0 d θ V l ( θ ) L Φ l ( t + θ ) ( p ) = L D 1 2 Q ( t ) ( p ) + I 1 2 ( p ) I 1 2 ( p ) l = 1 r h 0 d θ V l ( θ ) L Φ l ( t + θ ) ( p ) = L D 1 2 Q ( t ) ( p ) l = 1 r h 0 d θ V l ( θ ) I 1 2 ( p ) L Φ l ( t + θ ) ( p ) + l = 1 r h 0 d θ V l ( θ ) L Φ l ( t + θ ) ( p ) .
For the first term in the right-hand side of Equation (15) we have
L D 1 2 Q ( t ) ( p ) l = 1 r h 0 d θ V l ( θ ) I 1 2 ( p ) L Φ l ( t + θ ) ( p ) = l = 1 r h 0 d θ V l ( θ ) L D 1 2 Q ( t ) ( p ) L D 1 2 Φ l ( t + θ ) ( p )
and hence, from Equation (15) it follows
I 1 ( p ) L Q ( t ) ( p ) l = 1 r h 0 d θ V l ( θ ) θ 0 exp ( p ( θ η ) ) Φ l ( η ) d η = l = 1 r h 0 d θ V l ( θ ) ( L D 1 2 Q ( t ) ( p ) L D 1 2 Φ l ( t + θ ) ( p ) + l = 1 r h 0 [ d θ V l ( θ ) ] L Φ l ( t + θ ) ( p ) .
Analogously, for the third term in the right-hand side of Equation (14), we have
I 1 α ( p ) L Q ( t ) ( p ) i = 0 m h 0 d θ U i ( θ ) θ 0 exp ( p ( θ η ) ) Φ i ( η ) d η = I 1 α ( p ) L Q ( t ) ( p ) i = 0 m h 0 d θ U i ( θ ) L Φ i ( t + θ ) ( p ) = i = 0 m h 0 d θ U i ( θ ) L D 1 α Q ( t ) ( p ) L Φ i ( t + θ ) ( p ) + i = 0 m h 0 [ d θ U i ( θ ) ] L D α Φ i ( t + θ ) ( p )
From Equations (14), (16) and (17), it follows that
L X Φ ( p ) = L Q ( t ) ( p ) Φ ( 0 ) l = 1 r h 0 d θ V l ( θ ) Φ ( θ ) + l = 1 r h 0 d θ V l ( θ ) L Φ l ( t + θ ) ( p ) + l = 1 r h 0 d θ V l ( θ ) L D 1 2 Q ( t ) ( p ) L D 1 2 Φ l ( t + θ ) ( p ) + i = 0 m h 0 d θ U i ( θ ) L D 1 α Q ( t ) ( p ) L Φ i ( t + θ ) ( p ) + i = 0 m h 0 d θ U i ( θ ) L D α Φ i ( t + θ ) ( p )
Applying the inverse Laplace transform to both sides of Equation (18), we obtain Equation (10). □
Theorem 3.
Let the following conditions be satisfied:
(i)
The conditions (SA) hold.
(ii)
The function F L 1 loc ( R ¯ + , R n ) is exponentially bounded.
Then the solution X F ( t ) of the IVP (1), (3) with initial function Φ ( t ) 0 , t [ h , 0 ] has the following representation:
X F ( t ) = 0 t D 1 α Q ( t s ) F ( s ) d s + D α F ( t ) ,
where Q ( t ) is the fundamental matrix of System (2).
Proof. 
First we substitute X Φ ( t ) in Equation (1) and use the fact that X F ( t ) = 0 , t [ h , 0 ] . Since the function F is exponentially bounded, then we can apply to both sides the Laplace transform in order to get
L X F ( t ) ( p ) I α ( p ) I α ( p ) l = 1 r h 0 exp ( p θ ) d θ V l ( θ ) i = 0 m h 0 exp ( p θ ) d θ U i ( θ ) = L X F ( t ) ( p ) G ( p ) = L F ( t ) ( p ) .
Now it follows from the equality G 1 ( p ) = I 1 α ( p ) L Q ( t ) ( p ) that
L X F ( t ) ( p ) = I 1 α ( p ) I α 1 ( p ) G 1 ( p ) L F ( t ) ( p ) = I 1 α ( p ) L Q ( t ) ( p ) L F ( t ) ( p ) = L D 1 α Q ( t ) ( p ) + I α ( p ) L F ( t ) ( p ) = L D 1 α Q ( t ) ( p ) L F ( t ) ( p ) + L D α F ( t ) ( p ) .
Finally, we apply the inverse Laplace transform to Equation (21) and the representation Equation (19) follows. □
Corollary 1.
Let the conditions of Theorem 3 hold. Then for every initial function Φ C ˜ , the corresponding unique solution X Φ F ( t ) of the IVP (1), (3) has the integral representation
X ϕ F ( t ) = 0 t D 1 α Q ( t s ) F ( s ) d s + D α F ( t ) + Q ( t ) Φ ( 0 ) l = 1 r h 0 d θ V l ( θ ) Φ ( θ ) + l = 1 r h 0 d θ V l ( θ ) Φ ( t + θ ) + l = 1 r h 0 d θ V l ( θ ) D 1 2 Q ( t ) D 1 2 Φ l ( t + θ ) + i = 0 m h 0 d θ U i ( θ ) D 1 α Q ( t ) Φ i ( t + θ ) + i = 0 m h 0 d θ U i ( θ ) D α Φ i ( t + θ ) ,
where Q ( t ) is the fundamental matrix of System (2).
Proof. 
Let Φ C ˜ be an arbitrary initial function and let the functions X Φ ( t ) and X F ( t ) be defined by the Equalities (9) and (19), respectively. Then, according to the superposition principle, the function X Φ ( t ) + X F ( t ) is the unique solution of the IVP (1), (3). Now the statement of Corollary 1 follows immediately from Theorems 2 and 3. □

4. Example

First, we give some results needed for the illustrative example presented below:
The delayed Mittag-Leffler type matrix function E α , 1 B , τ : R R n × n for every matrix B R n × n and for τ R + is defined by
E τ B t α ( t ) : = I + k = 1 B k ( t ( k 1 ) τ ) α k Γ ( α k + 1 ) H ( k τ t ) , t 0
with E τ B t α ( 0 ) : = I , E τ B t α ( t ) : = Θ for t < 0 and H ( t ) is the Heavyside function with H ( 0 ) = 1 . This is a slight modification of the original definition in [29], and note that for each t 0 , the sum in Equation (22) is finite and for τ = 0 we have
E 0 B t α ( t ) : = E α ( B t α ) = k = 0 B k t α k Γ ( α k + 1 ) , t 0 ,
where the right side is the standard Mittag-Leffler type matrix function.
Example 1.
Consider the nonhomogeneous system for t > 0 :
D 0 + 0.5 x 1 ( t ) = x 1 ( t 1 ) + 1 D 0 + 0.5 ( x 2 ( t ) + x 1 ( t 1 ) + x 2 ( t 1 ) ) = x 2 ( t ) + x 2 ( t 1 ) + x 1 ( t 2 )
with the initial conditions
Φ ( t ) = ( 0 , 2 ) T , t [ 2 , 0 ] i . e . x 1 ( t ) = 0 , x 2 ( t ) = 2 f o r t [ 2 , 0 ] .
The homogenious system has the form
D 0 + 0.5 x ¯ 1 ( t ) = x ¯ 1 ( t 1 ) D 0 + 0.5 ( x ¯ 2 ( t ) + x ¯ 1 ( t 1 ) + x ¯ 2 ( t 1 ) ) = x ¯ 2 ( t ) + x ¯ 2 ( t 1 ) + x ¯ 1 ( t 2 )
and introduce the following initial conditions necessary for the calculating the fundamental matrix Q ( t ) :
1 . x 1 ( 0 ) = 1 , x 2 ( 0 ) = 0 a n d x 1 ( t ) = 0 , x 2 ( t ) = 0 f o r t [ 2 , 0 ) ;
2 . x 1 ( 0 ) = 0 , x 2 ( 0 ) = 1 a n d x 1 ( t ) = 0 , x 2 ( t ) = 0 f o r t [ 2 , 0 ) .
Let consider the IP (26), (27). Then the first Equation of (26) in virtue of Theorem 3.1 in [29] has the solution x ¯ 1 τ ( t ) = E 1 t 0 . 5 ( τ = 1 , α = 0.5 ). Taking into account Equation (27), it is simple to check that ( D 0 + 0 . 5 x ¯ 1 1 ( s 1 ) ) ( t ) = ( D 0 + 0 . 5 x ¯ 1 1 ( s ) ) ( t 1 ) , and then in virtue of Theorem 3.1 in [29] we have that ( D 0 + 0 . 5 x ¯ 1 1 ) ( t 1 ) = x ¯ 1 1 ( t 2 ) , and hense, from the second equation and Equation (27), we obtain that x ¯ 2 1 ( t ) 0 for t [ 2 , ) . Thus the IP (26), (27) have the following solution x ¯ 1 1 ( t ) = E 1 t 0 . 5 , x ¯ 2 1 ( t ) 0 for t [ 1 , ) .
Consider the IP (26), (28). Then obviously x ¯ 1 2 ( t ) 0 for t [ 2 , ) and the second equation become the form: D 0 + 0.5 ( x ¯ 2 ( t ) + x ¯ 2 ( t 1 ) ) = x ¯ 2 ( t ) + x ¯ 2 ( t 1 ) and by making the substitutuon y ( t ) = x ¯ 2 ( t ) + x ¯ 2 ( t 1 ) we obtain the equations D 0 + 0.5 y ( t ) = y ( t ) with initial codition y ( 0 ) = 1 , i.e., the following IP
D 0 + 0.5 y ( t ) = y ( t ) , t > 0 ; y ( 0 ) = 1 .
Applying Lemma 2.23 in [2] for the case when λ = 1 , τ = 1 , α = 0.5 we obtain that the solution of the IP (29) is the fuction y ( t ) = E 1 t 0.5 ( t ) = k = 0 t α k Γ ( α k + 1 ) . Then, using the step method, we obtain for each k N and t [ k 1 , k ) ] that x ¯ 2 2 ( t ) = k = 1 ( 1 ) k 1 E 1 t 0.5 ( t ( k 1 ) ) H ( k t ) for t > 0 . Thus, we obtain that the fundamental matrix have the form:
Q ( t ) = E 1 t 0 . 5 ( t ) 0 0 k = 1 ( 1 ) k 1 E 1 t 0 . 5 ( t ( k 1 ) ) H ( k t ) .
In the IP (24), (25) we have that: Φ ( t ) = ( 0 , 2 ) T , t [ 2 , 0 ] ; F ( t ) = ( 1 , 0 ) T . Then from Equation (19), we have
x 1 F ( t ) = 1 Γ ( 0.5 ) 0 t 0 t s ( t s η ) 0.5 E 1 t 0.5 ( η ) d η d s + t Γ ( 1.5 ) , x 2 F ( t ) = 0 .
From Equation (10), it follows
x 1 Φ ( t ) = 0 , x 2 Φ ( t ) = 2 + 2 0 t 0 s ( s η ) 0.5 k = 1 ( 1 ) k 1 E 1 t 0.5 ( η ( k 1 ) ) H ( k η ) d η d s + 2 t Γ ( 1.5 ) .
Then, the solution of the IP (24), (25), according Corollary 1, is
x 1 ( t ) = x 1 Φ ( t ) + x 1 F ( t ) = 1 Γ ( 0.5 ) 0 t 0 t s ( t s η ) 0.5 E 1 t 0.5 ( η ) d η d s + t Γ ( 1.5 ) ,
x 2 ( t ) = x 2 Φ ( t ) + x 2 F ( t ) = 2 + 2 0 t 0 s ( s η ) 0.5 k = 1 ( 1 ) k 1 E 1 t 0.5 ( η ( k 1 ) ) H ( k η ) d η d s + 2 t Γ ( 1.5 ) .

5. Conclusions

Following the investigations way in the case of functional differential systems with integer order derivatives, we proved a formula for integral representation of the solutions of Cauchy problem for fractional neutral systems, which improves and extends the corresponding former results obtained in the particular case of fractional systems with constant delays. However, the main idea is not only to make a standard generalization of existing results, but as in the case of systems with integer derivatives, the proved formula to be an useful tool for further study of different kinds stability properties of linear neutral fractional systems, which have a lot of practical applications.
As examples in this direction, we refer to the works [29,30], where finite time stability is studied by this approach, i.e., in the partial case of one constant delay. In the mentioned articles, first a formula for integral representation of the solutions of Cauchy problem is proved, and then, using the obtained result, sufficient conditions for finite time stability of the considered fractional delayed system are established. Furthermore, applying the same approach, in [16], the asymptotic stability properties of nonlinear perturbed linear fractional delayed systems are studied.

Author Contributions

The authors contribution in the article are equal. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors are grateful to the anonymous reviewers for their very helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Diethelm, K. The Analysis of Fractional Differential Equations, an Application–Oriented Exposition Using Differential Operators of Caputo Type; Lecture Notes in Mathematics; Springer: Berlin, Germany, 2010; Volume 2004. [Google Scholar]
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science B.V: Amsterdam, The Netherlands, 2006. [Google Scholar]
  3. Kiryakova, V. Generalized Fractional Calculus and Applications; Longman Scientific & Technical: Harlow, UK; John Wiley & Sons, Inc.: New York, NY, USA, 1994. [Google Scholar]
  4. Podlubny, I. Fractional Differential Equation; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  5. Kaczorek, T. Selected Problems of Fractional Systems Theory; Lecture Notes in Control and Information Sciences; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  6. Stamova, I.; Stamov, G. Functional and Impulsive Differential Equations of Fractional Order; Qualitative analysis and applications; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar]
  7. Bonilla, B.; Rivero, M.; Trujillo, J. On systems of linear fractional differential equations with constant coefficients. Appl. Comput. Math. 2007, 187, 68–78. [Google Scholar] [CrossRef]
  8. Li, K.; Peng, J. Laplace transform and fractional differential equations. Appl. Math. Let. 2011, 24, 2019–2023. [Google Scholar] [CrossRef]
  9. Lin, S.D.; Lu, C.H. Laplace transform for solving some families of fractional differential equations and its applications. Adv. Differ. Equ. 2013, 137, 1–9. [Google Scholar] [CrossRef] [Green Version]
  10. Golev, A.; Milev, M. Integral representation of the solution of the Cauchy problem for autonomous linear neutral fractional system. Int. J. Pure Appl. Math. 2018, 119, 235–247. [Google Scholar]
  11. Kiskinov, H.; Zahariev, A. Asymptotic stability of delayed fractional systems with nonlinear perturbation. AIP Conf. Proc. 2018, 2048, 050014. [Google Scholar]
  12. Zhang, H.; Cao, J.; Jiang, W. General solution of linear fractional neutral differential difference equations. Discret. Dyn. Nat. Soc. 2013. [Google Scholar] [CrossRef]
  13. Zhang, H.; Wu, D. Variation of constant formulae for time invariant and time varying Caputo fractional delay differential systems. J. Math. Res. Appl. 2014, 34, 549–560. [Google Scholar]
  14. Kiskinov, H.; Zahariev, A. On fractional systems with Riemann-Liouville derivatives and distributed delays—Choice of initial conditions, existence and uniqueness of the solutions. Eur. Phys. J. Spec. Top. 2017, 26, 3473–3487. [Google Scholar] [CrossRef]
  15. Krol, K. Asymptotic properties of fractional delay differential equations. Appl. Math. Comput. 2011, 218, 1515–1532. [Google Scholar]
  16. Milev, M.; Zlatev, S. A note about stability of fractional retarded linear systems with distributed delays. Int. J. Pure Appl. Math. 2017, 115, 873–881. [Google Scholar] [CrossRef] [Green Version]
  17. Veselinova, M.; Kiskinov, H.; Zahariev, A. Stability analysis of linear fractional differential system with distributed delays. AIP Conf. Proc. 2015, 1690, 040013. [Google Scholar]
  18. Veselinova, M.; Kiskinov, H.; Zahariev, A. About stability conditions for retarded fractional differential systems with distributed delays. Commun. Appl. Anal. 2016, 20, 325–334. [Google Scholar]
  19. Kiskinov, H.; Milev, N.; Zahariev, A. A comparison type theorem for linear neutral fractional systems with distributed delays. AIP Conf. Proc. 2017, 1910, 050009. [Google Scholar]
  20. Veselinova, M.; Kiskinov, H.; Zahariev, A. Stability analysis of neutral linear fractional system with distributed delays. Filomat 2016, 30, 841–851. [Google Scholar] [CrossRef] [Green Version]
  21. Veselinova, M.; Kiskinov, H.; Zahariev, A. Explicit conditions for stability of neutral linear fractional system with distributed delays. AIP Conf. Proc. 2016, 1789, 040005. [Google Scholar]
  22. Boyadzhiev, D.; Kiskinov, H.; Veselinova, M.; Zahariev, A. Stability analysis of linear distributed order fractional systems with distributed delays. Fract. Calc. Appl. Anal. 2017, 20, 914–935. [Google Scholar] [CrossRef]
  23. Martinez-Garcia, M.; Gordon, T.; Shu, L. Extended crossover model for human-control of fractional order plants. IEEE Access 2017, 5, 27622–27635. [Google Scholar] [CrossRef]
  24. Martinez-Garcia, M.; Zhang, Y.; Gordon, T. Memory pattern identification for feedback tracking control in human-machine systems. Hum. Factors 2019. [Google Scholar] [CrossRef] [Green Version]
  25. Boyadzhiev, D.; Kiskinov, H.; Zahariev, A. Integral representation of solutions of fractional system with distributed delays. Integral Transforms Spec. Funct. 2018, 29, 725–744. [Google Scholar] [CrossRef]
  26. Zahariev, A.; Kiskinov, H.; Angelova, E. Linear fractional system of incommensurate type with distributed delay and bounded Lebesgue measurable initial conditions. Dyn. Syst. Appl. 2019, 28, 491–506. [Google Scholar]
  27. Zahariev, A.; Kiskinov, H.; Angelova, E. Smoothness of the fundamental matrix of linear fractional system with variable delays. Neural Parallel Sci. Comput. 2019, 27, 71–83. [Google Scholar]
  28. Zahariev, A.; Kiskinov, H. Existence of fundamental matrix for neutral linear fractional system with distributed delays. Int. J. Pure and Appl. Math. 2018, 119, 31–51. [Google Scholar]
  29. Li, M.; Wang, J. Finite time stability of fractional delay differential equations. Appl. Math. Lett. 2017, 64, 170–176. [Google Scholar] [CrossRef]
  30. Li, M.; Wang, J. Exploring delayed Mittag-Leffler type matrix functions to study finite time stability of fractional delay differential equations. Appl. Mathe. Comput. 2018, 324, 254–265. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Madamlieva, E.; Konstantinov, M.; Milev, M.; Petkova, M. Integral Representation for the Solutions of Autonomous Linear Neutral Fractional Systems with Distributed Delay. Mathematics 2020, 8, 364. https://doi.org/10.3390/math8030364

AMA Style

Madamlieva E, Konstantinov M, Milev M, Petkova M. Integral Representation for the Solutions of Autonomous Linear Neutral Fractional Systems with Distributed Delay. Mathematics. 2020; 8(3):364. https://doi.org/10.3390/math8030364

Chicago/Turabian Style

Madamlieva, Ekaterina, Mihail Konstantinov, Marian Milev, and Milena Petkova. 2020. "Integral Representation for the Solutions of Autonomous Linear Neutral Fractional Systems with Distributed Delay" Mathematics 8, no. 3: 364. https://doi.org/10.3390/math8030364

APA Style

Madamlieva, E., Konstantinov, M., Milev, M., & Petkova, M. (2020). Integral Representation for the Solutions of Autonomous Linear Neutral Fractional Systems with Distributed Delay. Mathematics, 8(3), 364. https://doi.org/10.3390/math8030364

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