Next Article in Journal
Including Jumps in the Stochastic Valuation of Freight Derivatives
Next Article in Special Issue
Some Existence and Dependence Criteria of Solutions to a Fractional Integro-Differential Boundary Value Problem via the Generalized Gronwall Inequality
Previous Article in Journal
Weighted Wiener Indices of Molecular Graphs with Application to Alkenes and Alkadienes
Previous Article in Special Issue
Stability Results for Two-Dimensional Systems of Fractional-Order Difference Equations
 
 
Correction published on 3 June 2021, see Mathematics 2021, 9(11), 1282.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence of Absolutely Continuous Fundamental Matrix of Linear Fractional System with Distributed Delays

by
Hristo Kiskinov
*,
Ekaterina Madamlieva
,
Magdalena Veselinova
and
Andrey Zahariev
Faculty of Mathematics and Informatics, University of Plovdiv, 4000 Plovdiv, Bulgaria
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(2), 150; https://doi.org/10.3390/math9020150
Submission received: 23 December 2020 / Revised: 7 January 2021 / Accepted: 10 January 2021 / Published: 12 January 2021
(This article belongs to the Special Issue Qualitative Theory of Fractional-Order Systems)

Abstract

:
The goal of the present paper is to obtain sufficient conditions that guaranty the existence and uniqueness of an absolutely continuous fundamental matrix for a retarded linear fractional differential system with Caputo type derivatives and distributed delays. Some applications of the obtained result concerning the integral representation of the solutions are given too.
MSC:
34A08; 34A30; 26A33; 34A12

1. Introduction

It is well known that many mathematical models of real world phenomena can be described more accurately through fractional derivative formulation. For more details on fractional calculus theory and fractional differential equations, we recommend the monographs of Kilbas et al. [1] and Podlubny [2]. For distributed order fractional differential equations see Jiao at all [3] and for an application-oriented exposition Diethelm [4]. Impulsive differential and functional differential equations with fractional derivative and some applications are studied by Stamova and Stamov [5].
The theme of the integral representation (variation of constants formula) of the solutions of linear fractional differential equations and/or systems (ordinary or with delay) is an “evergreen” theme for research. This explains why a lot of papers are devoted to different aspects of this problem. For linear fractional ordinary differential equations and systems, we refer the works [1,2,6,7,8,9,10] and the references therein. Relatively, as far as we know, there are not many works devoted to the variation of constants formula for delayed linear fractional systems [11,12,13,14,15]. The case of neutral fractional systems is studied in [16,17,18,19,20].
The establishment of a fundamental matrix with appropriate properties (for example in [21] smoothness is obtained) is the basis for obtaining any integral representation and is a key tool in the study of different types of stability of linear and nonlinear disturbed systems (see [20]).
In the present work, we consider linear fractional systems with distributed delays and incommensurate order derivatives in the Caputo sense. The first goal of the work is to establish sufficient conditions for existence and uniqueness of a fundamental matrix C ( t , s ) , which is absolutely continuous in t on every compact subinterval of R . The second one is to clarify the analytic properties in s, which are very similar to these in the integer case. As an application of the obtained results, some results concerning the integral representation of the solutions given in [15,18] are improved.
The paper is organized as follows. In Section 2 we recall the definitions of Riemann–Liouville and Caputo fractional derivatives with some of their properties. In the same section is the statement of the problem, as well as some necessary preliminary results used later. Section 3 is devoted to the existence and the uniqueness of the solutions of the Initial Problem (IP) for linear fractional systems with distributed delays and incommensurate order derivatives in the Caputo sense and special type discontinuous initial function. In Section 4 the existence and uniqueness of an absolutely continuous fundamental matrix is proved und its analytical properties are studied. Using the obtained results, in Section 5 we establish new integral representations for the solutions of the studied systems.

2. Preliminaries and Problem Statement

For readers convenience, below we recall the definitions of Riemann–Liouville and Caputo fractional derivatives as well as some needed properties. For details and other properties we refer to [1,2,3].
Let α ( 0 , 1 ) be an arbitrary number and denote by L 1 l o c ( R , R ) the linear space of all locally Lebesgue integrable functions f : R R . Then for each a R , each t > a and f L 1 l o c ( R , R ) the definitions of the left-sided fractional integral operator, the left side Riemann–Liouville and Caputo fractional derivatives of order α and some properties are given below (see [1]):
( 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 ) ;
C D a + α f ( t ) = R L D a + α f ( t ) f ( a ) Γ ( 1 α ) ( t a ) α ;
( a ) C D a + α D a + α f ( t ) = ( D a + 0 f ) ( t ) = f ( t ) ; ( b ) D a + α C D a + α f ( t ) = f ( t ) f ( a )
If f A C ( R , R ) then the next formula gives a direct definition of the Caputo left side derivative:
C D a + α f ( t ) = 1 Γ ( 1 α ) a t ( t s ) α f ( s ) d s ,
Everywhere below, the following notations will be used: R + = ( 0 , ) , R ¯ + = [ 0 , ) , J a = [ a , ) , a R , J s + M = [ s , s + M ] , s J a , M R + , n = { 1 , 2 , , n } , m 0 = m { 0 } , n , m N , I , Θ R n × n denote the identity and zero matrix, respectively, I k , k n denotes the k-th column of the identity matrix and 0 R n is the zero element.
For Y ( t ) = ( y 1 ( t ) , , y n ( t ) ) T : J a R n , β = ( β 1 , , β n ) , β k [ 1 , 1 ] , k n we will use the notations I β ( Y ( t ) ) = diag ( ( y 1 ( t ) ) β 1 , , ( y n ( t ) ) β n ) , for W ( t ) = { w k j ( t ) } k , j = 1 n : J a R n × n we denote | W ( t ) | = k , j = 1 n | w k , j ( t ) | , t J a and for simplicity we will use the notation D a + α = C D a + α for the left side Caputo fractional derivative.
Consider the homogeneous linear delayed system of incommensurate type and distributed delay in the following general form
D a + α X ( t ) = h 0 [ d θ U ( t , θ ) ] X ( t + θ ) , t > a
or described in more detailed form
D a + α k x k ( t ) = j = 1 n h 0 x j ( t + θ ) d θ u k j ( t , θ ) , k n , t > a
and the corresponding nonhomogeneous one
D a + α X ( t ) = h 0 [ d θ U ( t , θ ) ] X ( t + θ ) + F ( t ) , t > a
where h R + , X ( t ) = ( x 1 ( t ) , , x n ( t ) ) T , F ( t ) = ( f 1 ( t ) , , f n ( t ) ) T : J a R n , α = ( α 1 , , α n ) , α k ( 0 , 1 ) , α m = min k n α k , D a + α = diag ( D a + α 1 , , D a + α n ) T , U : J a × R R n × n , U ( t , θ ) = { u k j ( t , θ ) } k , j = 1 n .
Definition 1.
A function Z ( t ) = ( z 1 ( t ) , , z n ( t ) ) T : J R n is called piecewise absolutely continuous on some interval J R (denoted Z ( t ) P A C ( J , R n ) ), if it is with bounded variation (BV) in t on J, there is no singular term in the Lebesgue decomposition of Z ( t ) and the set of discontinuity points of the function Z ( t ) has no limit points in J.
Definition 2.
With C a , a R we denote the Banach space of right continuous vector functions Φ ( t ) P A C ( [ a h , a ] , R n ) with norm | | Φ | | = sup t [ a h , a ] | Φ ( t ) | = sup t [ a h , a ] k = 1 n | ϕ k ( t ) | < and the subspace of all absolutely continuous functions by C a = A C ( [ a h , a ] , R n ) , i.e., C a C a .
The set of jump points of every initial vector function Φ C a we denote by S Φ . We emphasize that the set S Φ K is finite for every compact interval K R and the case S Φ = is not excluded.
For the system (1) or (2) introduce the following initial conditions:
X ( t ) = Φ ( t ) ( x k ( t ) = ϕ k ( t ) , k n ) , t [ a h , a ] , Φ C a
We say that for the kernel U : J a × R R n × n the conditions (S) hold, if the following conditions are fulfilled:
(S1)
The functions ( t , θ ) U ( t , θ ) are measurable in ( t , θ ) J a × R and normalized so that for t J a , U ( t , θ ) = 0 when θ R ¯ + and U ( t , θ ) = U ( t , h ) for all θ ( , h ] , h R + and V a r θ [ h , 0 ] U ( t , · ) < for t J a .
(S2)
The Lebesgue decomposition of the kernel U ( t , θ ) for t J a and θ [ h , 0 ] has the form:
U ( t , θ ) = U J ( t , θ ) + U A C ( t , θ ) + U S ( t , θ )
where U J ( t , θ ) = i = 0 m A i ( t ) H ( θ + σ i ( t ) ) , m N , A i ( t ) = { a k j i ( t ) } k , j = 1 n L 1 l o c ( J a , R n × n ) are locally bounded on J a , H ( t ) is the Heaviside function, the delays σ i ( t ) C ( J a , R ¯ + ) are bounded with σ i = sup t J a σ i ( t ) , max i m 0 σ i h and A i ( t ) H ( θ + σ i ( t ) ) are continuous from left in θ on ( σ i , 0 ) , i m , σ 0 ( t ) 0 .
U A C ( t , θ ) = h 0 B ( t , θ ) d θ , B ( t , θ ) = { b k j ( t , θ ) } k , j = 1 n L 1 l o c ( J a × R , R n × n ) are locally bounded on J a and U S ( t , θ ) C ( J a × R , R n × n ) .
(S3)
There exists a locally bounded function z u L 1 l o c ( J a , R + ) such that V a r [ h , 0 ] U ( t , · ) z u ( t ) for t J a and for every t J a the following relations hold: lim t t h 0 | U ( t , θ ) U ( t , θ ) | d θ = 0 .
(S4)
The sets S Φ i = { t J a | t σ i ( t ) S Φ } for every i m do not have limit points.
Consider the following auxiliary system in matrix form
X ( t ) = Φ ( a ) + I 1 ( Γ ( α ) ) a t I α 1 ( t η ) F ( η ) d η + I 1 ( Γ ( α ) ) a t I α 1 ( t η ) h 0 [ d θ U ( η , θ ) ] X ( η + θ ) d η , t > a
where I 1 ( Γ ( α ) ) = diag ( Γ 1 ( α 1 ) , , Γ 1 ( α n ) ) , or for k n in more detailed form
x k ( t ) = ϕ k ( a ) + 1 Γ ( α k ) a t ( t η ) α k 1 f k ( η ) d η + 1 Γ ( α k ) a t ( t η ) α k 1 ( j = 1 n h 0 x j ( η + θ ) d θ u k j ( η , θ ) ) d η , t > a
with the initial condition (3).
Definition 3.
The vector function X ( t ) = ( x 1 ( t ) , , x n ( t ) ) T is a solution of the IP (2) and (3) or IP (3) and (4) in J a + M ( J a ) , M R + if X C ( [ a , a + M ] , R n ) ( X C ( J a , R n ) ) satisfies the system (2), respectively, (4) for all t ( a , a + M ] ( t ( a , ) ) and the initial condition (3) for each t [ a h , a ] .
In virtue of Lemma 3.3 in [15] every solution X ( t ) of IP (2) and (3) is a solution of the IP (3) and (4) and vice versa.
We will need a slightly modified version of the Weissinger generalization of the Banach’s fixed point theorem for complete metric spaces (see [22], Fixpunktsatz, p. 195).
Theorem 1.
Let Ω be a complete metric space with metric d Ω and let the following conditions hold:
1. 
There exists a sequence γ q 0 , q N , with q = 1 γ q < .
2. 
The operator T : Ω Ω satisfies for each q N and for arbitrary x , y Ω the inequality
d Ω ( T q x , T q y ) γ q d Ω ( x , y )
Then T has a uniquely fixed point x Ω ( x = T x ) and for every x Ω we have that lim q T q x = x .
Remark 1.
This modification of the Weissinger generalization of the Banachs fixed point is not new. It is used in [23] and in the case when Ω is a Banach space in [24]. It is clear that the original Weissinger proof is correct for the presented variant too, with elementary modifications.
Let B be an arbitrary real Banach space.
Definition 4 ([25]).
The function f ( t ) : R B is called a regulated function if it has one-sided (left and right) limits at every point t R .
Remark 2.
If f ( t ) : J K B , where J K R is an arbitrary compact interval and f ( t ) is a regulated function, then it is assumed that in the left (right) end on the interval J K the function f ( t ) has only a right (left) limit.
Theorem 2 ([25]).
Let f ( t ) : R B be an arbitrary function.
Then a necessary and sufficient condition for f ( t ) to be a regulated function is that f ( t ) must in every compact interval J K R be a limit of n uniformly convergent sequence of step-functions (i.e., with respect to the supremum norm · = sup t J K · B ).
Theorem 3 ([26]).
Let t 0 J a and α > 0 be arbitrary fixed numbers and the following conditions hold:
1. 
The functions p ( t ) , u ( t ) L l o c 1 ( [ t 0 , T ) , R ¯ + ) for some T .
2. 
The function g ( t ) C ( [ t 0 , T ) , [ 0 , M ] ) for some M R + and is nondecreasing.
3. 
For every t J a is fulfilled u ( t ) p ( t ) + g ( t ) t 0 t ( t η ) α 1 u ( η ) d η .
Then for t J a the following inequality holds
u ( t ) p ( t ) + t 0 t [ q = 1 ( g ( η ) ( Γ ( α ) ) q Γ ( α q ) ( t η ) α q 1 ] p ( η ) d η .
Remark 3.
Note that the statement of Theorem 3 is proved in the partial case t 0 = 0 , but with small modifications the proof will be correct for arbitrary t 0 J a .

3. Existence and Uniqueness of the Solutions

Let it be that for every Φ C a consider the corresponding linear space
E Φ = { G : J a R n | G | J a A C ( J a , R n ) , G ( t ) = Φ ( t ) , t [ a h , a ] }
For each M R + define the set
E M Φ = { G M : J a + M R n | G M = G | [ a , a + M ] , G E Φ ) ,
where Φ C a is arbitrary and define a metric function d M : E M Φ × E M Φ R ¯ + with
d M ( G M , G M ) = k = 1 n sup t J a + M | g k ( t ) g k ( t ) |
for each G M , G M E M Φ .
Since G M ( a ) = Φ ( a ) = G M ( a ) then for every M R + , according to a well-known result we conclude that E M Φ is a complete metric space concerning the metric
d M V a r ( G M , G M ) = V a r [ a , a + M ] ( G M ( t ) G M ( t ) ) = k = 1 n V a r [ a , a + M ] ( g k ( t ) g k ( t ) ) .
It was a very strange for us that we could not find a result from which the statement of the next lemma directly follows.
Lemma 1.
For every M R + the set E M Φ is a complete metric space concerning the metric d M too.
Proof. 
Let M R + be an arbitrary fixed number and consider an arbitrary Cauchy sequence { G M l ( t ) = ( g 1 l ( t ) , , g 1 l ( t ) ) T } l = 1 E M Φ , i.e., lim l , r d M ( G M l ( t ) , G M r ( t ) ) = 0 . It is clear that there exists a vector valued function G M 0 ( t ) = ( g 1 0 ( t ) , , g 1 0 ( t ) ) T C ( J a + M , R n ) , G M 0 ( t ) = Φ ( t ) , t [ a h , a ] such that lim l d M ( G M l ( t ) , G M 0 ( t ) ) = 0 .
Let ε > 0 and k N be arbitrary numbers. There exist n 0 = n 0 ( ε ) N and δ = δ ( ε , k ) > 0 such that for every n n 0 we have that | G n ( t ) G 0 ( t ) | < ε k for each t [ a , a + M ] . Since G n 0 A C ( [ a , a + M ] ) then there exists δ = δ ( ε , n 0 ) ( 0 , δ ) , such that for every finite sequence of pairwise disjoint subintervals { [ a j , b j ] } j k , j k [ a j , b j ] [ a , a + M ] with j = 1 k ( b j a j ) < δ the inequality j = 1 k | G n 0 ( a j ) G n 0 ( b j ) | < ε holds. Then whenever when j = 1 k ( b j a j ) < δ we have that
j = 1 k | G 0 ( a j ) G 0 ( b j ) | = j = 1 k | G 0 ( a j ) G 0 ( b j ) + G n 0 ( a j ) G n 0 ( a j ) + G n 0 ( b j ) G n 0 ( b j ) | j = 1 k | G 0 ( a j ) G n 0 ( a j ) | + j = 1 k | G n 0 ( b j ) G 0 ( b j ) | + j = 1 k | G n 0 ( a j ) G n 0 ( b j ) | 3 ε
Thus G 0 A C ( [ a , a + M ] ) for every M R + and the statement is proved. □
Remark 4.
It must be noted that the statement of the next theorem cannot be obtained directly as a corollary from analogical results for the considered initial problems in [15,18,23].
Theorem 4.
Let the following conditions hold:
1. 
Conditions (S) hold.
2. 
The function Φ C a is arbitrary.
Then the IP (1) and (3) has a unique absolutely continuous solution in the interval [ a , a + M ] for every M max { h , 1 } .
Proof. 
In the proof of this theorem we will use the approach introduced in [23].
Let Φ C a and M max { h , 1 } be arbitrary. From (5) in the case f k ( t ) 0 for k n we obtain the system
x k ( t ) = ϕ k ( a ) + 1 Γ ( α k ) a t ( t η ) α k 1 ( j = 1 n h 0 x j ( η + θ ) d θ u k j ( η , θ ) ) d η
For every G ( t ) = ( g 1 ( t ) , , g n ( t ) ) T E M Φ define the operator ( G ) ( t ) = ( 1 g 1 ( t ) ,
, n g n ( t ) ) T via the operators k g k ( t ) for t ( a , a + M ] , k n by
k g k ( t ) = ϕ k ( a ) + 1 Γ ( α k ) a t ( t s ) α k 1 ( j = 1 n h 0 g j ( s + θ ) d θ u k j ( s , θ ) ) d s ,
and the following additional condition:
k g k ( t ) = ϕ k ( t ) , t [ a h , a ]
Since Φ is PAC, the conditions (S) imply that for each k , j n , the functions t h 0 g j ( t + θ ) d θ u k , j ( t , θ ) are at least piecewise continuous on the interval [ a , a + M ] (see [27], Lemma 1). Then the right side of (7) is absolutely continuous on the interval [ a , a + M ] , which implies that the function k g k ( t ) is absolutely continuous on the same interval and since lim t a + k g k ( t ) = ϕ k ( a ) for k n , then ( G ) E M Φ . Thus the operator maps E M Φ into E M Φ .
Since according to Lemma 1 the space E M Φ is a full metric space concerning the metric d M , it is enough to check that for the operator the conditions of Theorem 1 hold and then in virtue of Theorem 1 we will obtain that the operator has a unique fixed point in E M Φ .
We recall that the Γ ( z ) , z R + , has a local minimum at z m i n 1.46163 , where it attains the value Γ ( z m i n ) 0.885603 . There exists q 0 N such that q 0 + 1 > α m 1 q 0 and for every q q 0 + 1 we denote with α q , α q { α 1 , , α n } that number for which Γ ( 1 + α q q ) = min k n Γ ( 1 + α k q ) .
Let denote P = max k n ( j = 1 n sup s [ a , a + M ] V a r θ [ h , 0 ] u k j ( s , θ ) ) and let G , G ¯ E M Φ be arbitrary. Then from (7) for k n and every t [ a , a + M ] we obtain
| k g k ( t ) k g ¯ k ( t ) | 1 Γ ( α k ) a t ( t s ) α k 1 ( j = 1 n | h 0 ( g j ( s + θ ) g ¯ j ( s + θ ) ) d u k j ( s , θ ) | ) d s ( t a ) α k Γ ( 1 + α k ) j = 1 n sup s [ a , a + M ] V a r θ [ h , 0 ] u k j ( s , θ ) sup t [ a , a + M ] | g j ( t ) g ¯ j ( t ) | ( t a ) α k P Γ ( 1 + α k ) j = 1 n sup t [ a , a + M ] | g j ( t ) g ¯ j ( t ) | ( t a ) α k P Γ ( 1 + α k ) d M ( G , G ¯ )
Let assume that for some q N , k n and for every t [ a , a + M ] the inequality
| k q g k ( t ) k q g ¯ k ( t ) | ( t a ) q α k P q Γ ( 1 + q α k ) d M ( G , G ¯ )
holds. Obviously from (9) it follows that the inequality (10) holds for each k n and every t [ a , a + M ] at least for q = 1 . Denoting for simplicity q G ( t ) = Y ( t ) = ( y ( t ) , , y n ( t ) ) T , q G ¯ ( t ) = Y ¯ ( t ) = ( y ¯ ( t ) , , y ¯ n ( t ) ) T , we obtain
| k q + 1 g k ( t ) k q + 1 g ¯ k ( t ) | = | k q g k ( t ) k q g ¯ k ( t ) | = | y k ( t ) y ¯ k ( t ) | .
Let us assume that (10) holds for each k n , every t [ a , a + M ] and for some q N . Then a similar way as in (9) from (10) and (11) we obtain
| k y k ( t ) y ¯ k ( t ) | 1 Γ ( α k ) a t ( t s ) α k 1 j = 1 n | h 0 ( y j ( s + θ ) y ¯ j ( s + θ ) ) d θ u k j ( s , θ ) | d s P Γ ( α k ) a t ( t s ) α k 1 j = 1 n sup η [ a , a + M ] | y j ( η ) y ¯ j ( η ) | d s P P q Γ ( α k ) Γ ( 1 + q α k ) d M ( G , G ¯ ) a t ( t s ) α k 1 ( s a ) q α k d s
Substituting s a = z ( t a ) in the integral in the right side of (12) and using the well-known relation between the gamma and beta functions we obtain
| k y k ( t ) y ¯ k ( t ) | P q + 1 ( t a ) α k ( q + 1 ) Γ ( α k ) Γ ( 1 + q α k ) d M ( G , G ¯ ) 0 1 ( 1 z ) α k 1 z q α k d z Γ ( α k ) Γ ( 1 + q α k ) ( P ( t a ) α k ) q + 1 Γ ( α k ) Γ ( 1 + q α k ) Γ ( 1 + ( q + 1 ) α k ) d M ( G , G ¯ ) ( P ( t a ) α k ) q + 1 Γ ( 1 + ( q + 1 ) α k ) d M ( G , G ¯ ) ( P M α k ) ( q + 1 ) Γ ( 1 + ( q + 1 ) α k ) d M ( G , G ¯ ) ( P M ) ( q + 1 ) Γ ( 1 + ( q + 1 ) α k ) d M ( G , G ¯ )
and hence (10) holds for every q N , for each k n and every t [ a , a + M ] .
For q q 0 + 1 from (10) and (13) it follows that
d M ( q + 1 G , q + 1 G ¯ ) n ( M P ) ( q + 1 ) Γ ( 1 + α q ( q + 1 ) ) d M ( G , G ¯ ) .
and denote γ q = n ( M P ) q Γ ( 1 + α q q ) for q q 0 + 1 .
For all q > q 0 + 1 from (10) and (13) we obtain that
d M ( q + 1 G , q + 1 G ¯ ) n ( M P ) ( q + 1 ) Γ ( 1 + α m ( q + 1 ) ) d M ( G , G ¯ ) .
and denote γ q = n ( M P ) q Γ ( 1 + α m q ) for q > q 0 + 1 .
Consider the one parameter Mittag–Leffler function E α m , 1 ( z ) = q = 1 z q Γ ( 1 + α m q ) , z R ¯ + . It is simply to see that the series q = 1 ( M P ) q Γ ( 1 + α m q ) is convergent because it is the considered Mittag–Leffler function evaluated at z = M P . Then we have that
q = 1 γ q = n ( q = 1 q 0 ( M P ) q Γ ( 1 + α q q ) + q = q 0 + 1 ( M P ) q Γ ( 1 + α m q ) ) <
and then from Theorem 1 it follows that the IP (1) and (3) has a unique solution in t [ a , a + M ] . □
Corollary 1.
Let the conditions of Theorem 4 hold.
Then the IP (1) and (3) has a unique absolutely continuous solution in the interval J a .
Proof. 
Let Φ C a be arbitrary and assume the contrary, that there exists M m a x < such that the solution X m a x ( t ) in the interval t [ a , a + M m a x ] does not possess a prolongation. Let M > M m a x be an arbitrary number. In virtue of Theorem 4 the IP (1) and (3) has a unique absolutely continuous solution X ( t ) in the interval t [ a , a + M ] . The solution X ( t ) obviously is a prolongation of the solution X m a x ( t ) , which contradicts of our assumption that M m a x < . □
For arbitrary fixed s > a consider the following auxiliary system
D a + α X ( t ) = h 0 [ d θ U ( t , θ ) ] X ( t + θ )
with the following condition
X ( t ) = Φ ˜ j s ( t ) , t ( , s ] .
Corollary 2.
Let the following conditions hold:
1. 
Conditions (S) hold.
2. 
The function Φ ˜ j s has the form
Φ ˜ j s ( t ) = I j , t = s , j n 0 , t < s
Then for each s > a and arbitrary j n the problem (16) and (17) has a unique solution, which satisfies Equation (16) for t > s , the condition (17) for t s and is absolutely continuous in ( , s ) ( s , ) with a first kind jump at t = s .
Proof. 
Let s > a be an arbitrary fixed number, introduce the system
D s + α X ( t ) = h 0 [ d θ U ( t , θ ) ] X ( t + θ ) , t > s
and consider the IP (17) and (18). Since Theorem 4 is proved for arbitrary a R and Φ ˜ j s ( t ) for arbitrary j n is PAC on the interval [ s h , s ] , then from Theorem 4 it follows that the IP (17) and (18) possess a unique solution X j ( · ) A C ( J s , R n ) and moreover, from (17) it follows also that X j ( · ) A C ( [ a , s ) J s , R n ) . Then for every k n we have
D a + α k x j k ( t ) = 1 Γ ( 1 α k ) a t ( t η ) α k x j k ( η ) d η = 1 Γ ( 1 α k ) a s ( t η ) α k x j k ( η ) d η + 1 Γ ( 1 α k ) s t ( t η ) α k x j k ( η ) d η = 1 Γ ( 1 α k ) s t ( t η ) α k x j k ( η ) d η = D s + α k x j k ( t )
and hence X j ( t ) = ( x j 1 , , x j n ) T satisfies the Equation (16) for t ( s , ) and the condition (17) for t ( , s ] .
Let consider an IP with Equation (16) for t > a and initial condition X ( t ) = X j ( t ) = 0 for t [ a h , a ] . Then obviously X j ( t ) = 0 is its unique solution in t [ a , s ) . This completes the proof. □

4. Fundamental Matrix

Let s J a be an arbitrary fixed number and define the following matrix valued function Φ ¯ ( t , s ) = ( φ ¯ k j ( t , s ) ) k , j = 1 n : R × J a R n × n with
Φ ¯ ( t , s ) = I , t = s Θ , t < s
and denote Φ ¯ j ( t , s ) = ( φ ¯ 1 j ( t , s ) , , φ ¯ n j ( t , s ) ) T , j n .
For arbitrary fixed number s J a consider the following matrix IP
D a + α C ( t , s ) = h 0 [ d θ U ( t , θ ) ] C ( t + θ , s ) , t ( s , )
C ( t , s ) = Φ ¯ ( t , s ) , t ( , s ] .
Definition 5.
The matrix valued function t C ( t , s ) = ( C 1 ( t , s ) , , C n ( t , s ) ) = { c k j ( t , s ) } k , j = 1 n , s J a , is called a solution of the IP (19) and (20) in J s if C ( · , s ) : [ s , ) R n × n is continuous for t [ s , ) and satisfies the matrix Equation (19) on t [ s , ) as well as the initial condition (20) too.
Remark 5.
Practically in condition (20) we need only the values of Φ ¯ ( · , s ) for t [ s h , s ] , but for convenience we define C ( t , s ) = Θ also for t ( , s h ) . Then C ( t , s ) is prolonged as continuous in t function on ( , s ) .
Theorem 5.
Let the conditions (S) hold.
Then for every initial point s J a , the matrix IP (19) and (20) has a unique absolutely continuous solution t C ( t , s ) in the interval J s .
Proof. 
The statement of the Theorem follows immediately from Corollary 2. □
Definition 6.
The matrix C ( t , s ) , which is a solution of the IP (19) and (20), will be called fundamental (or Cauchy) matrix for the homogeneous system (1).
Lemma 2.
Let the conditions (S) hold and the matrix valued function t C ( t , s ) is the fundamental matrix of the system (1).
Then for every t ¯ J a the matrix function C ( t , · ) : [ a , t ¯ ] R n × n is locally bounded in s for s [ a , t ¯ ] and t ( , t ¯ ] .
Proof. 
Let t ¯ J a be an arbitrary fixed number, s [ a , t ¯ ] be arbitrary and consider the fundamental matrix C ( t , · ) : [ a , t ¯ ] R n × n . According to Remark 5 for s [ a , t ¯ ] we have that C ( t , s ) = Θ for t ( , a ) and C ( a , a ) = I .
Taking into account Theorem 5 and Corollary 2 it is easy to be seen, that the unique solution C ( t , s ) of IP (19) and (20) is a solution of the equation
C ( t , s ) = Φ ¯ ( s , s ) + I 1 ( Γ ( α ) ) s t I α 1 ( t η ) ( h 0 [ d θ U ( η , θ ) ] C ( η + θ , s ) ) d η .
Introduce for every k , j n the notations: w k j ( t , s ) = sup ξ [ s , t ] | c k j ( ξ , s ) | = sup ξ [ a , t ] | c k j ( ξ , s ) | (since c k j ( ξ , s ) = 0 for ξ < s ), V U = sup t [ a , t ¯ ] V a r θ [ h , 0 ] U ( t , θ ) , b = | I 1 ( Γ ( α ) ) | V U and then from (21) it follows
| C ( t , s ) | | Φ ¯ ( s , s ) | + | I 1 ( Γ ( α ) ) | | s t I α 1 ( t η ) ( σ 0 [ d θ U ( η , θ ) ] C ( η + θ , s ) ) d η |
and hence for each k , j n
w k j ( t , s ) | φ ¯ k j ( s , s ) | + V U Γ ( α k ) s t ( t η ) α k 1 w k j ( η , s ) d η | φ ¯ k j ( s , s ) | + b a t ( t η ) α k 1 w k j ( η , s ) d η
From Theorem 3 and (22) for a t s t ¯ we obtain the estimation
w k j ( t , s ) 1 + a t ( q = 1 ( b Γ ( α k ) ) q Γ ( α k q ) ( t η ) α k q 1 ) d η 1 + a t ¯ ( q = 1 ( b Γ ( α k ) ) q Γ ( α k q ) ( t η ) α k q 1 ) d η
and thus C ( t , s ) is locally bounded in s for s [ a , t ¯ ] and t ( , t ¯ ] . □
Theorem 6.
Let the conditions (S) hold and the matrix valued function t C ( t , s ) be the fundamental matrix of system (1).
Then for every fixed t ¯ J a the matrix function C ( t , · ) : [ a , t ¯ ] R n × n is continuous for s [ a , t ¯ ] when s t , for s = t possess first kind jumps and hence is Lebesgue integrable in s on [ a , t ¯ ] for each t ( , t ¯ ] .
Proof. 
Let t ¯ J a be an arbitrary fixed number, s [ a , t ¯ ] be arbitrary and consider the fundamental matrix C ( t , · ) : [ a , t ¯ ] R n × n . According to Remark 5 for s ( a , t ¯ ] we have that C ( t , s ) = Θ for t ( , a ] and C ( a , a ) = I . In virtue of Theorem 2 and Lemma 2 it is enough to prove that C ( t , s ) has left and right limits for each s ( a , t ¯ ) and there exist the limits C ( t ¯ , t ¯ 0 ) and C ( a , a + 0 ) .
(i)
Let s [ a , t ¯ ] be arbitrary and let t < s . Then it is simply to see that for every s [ a , t ¯ ) lim s s + 0 C ( t , s ) = Θ = lim s s 0 C ( t , s ) . Note that for s = t ¯ we have that lim s t ¯ 0 C ( t , s ) = Θ holds.
(ii)
Let s [ a , t ¯ ] be arbitrary and let t = s . Then lim s s + 0 C ( s , s ) = Θ and since C ( s , s ) = I , then we can conclude that C ( s , s + 0 ) exists and hence C ( t , s ) has jumps of first kind on the line t = s for each s [ a , t ¯ ) .
(iii)
Let s [ a , t ¯ ] , t > s . Then we have to consider two cases: either s ( a , t ¯ ] , s ( a , s ) (left limit in s ) or s [ a , t ¯ ) with s ( s , t ) (right limit in s ).
(iii.a)
For purposes of clarity we assume that s ( a , t ¯ ] , s ( a , s ) . According to Corollary 2 for each j n , the IP (16) and (17) has unique solutions C j ( t , s ) and C j ( t , s ) for the initial functions Φ ¯ j ( t , s ) = ( φ ¯ 1 j ( t , s ) , , φ ¯ n j ( t , s ) ( t , s ) ) T and Φ ¯ j ( t , s ) = ( φ ¯ 1 j ( t , s ) , , φ ¯ n j ( t , s ) ) T , respectively.
Then we obtain
C j ( t , s ) = Φ ¯ j ( s , s ) + I 1 ( Γ ( α ) ) s t I α 1 ( t η ) ( h 0 [ d θ U ( η , θ ) ] C j ( η + θ , s ) ) d η
C j ( t , s ) = Φ ¯ j ( s , s ) + I 1 ( Γ ( α ) ) s t I α 1 ( t η ) ( h 0 [ d θ U ( η , θ ) ] C j ( η + θ , s ) ) d η .
Taking into account that Φ ¯ ( s , s ) = Φ ¯ ( s , s ) = I and subtracting both sides of (24) from the corresponding sides of (23) we obtain
C j ( t , s ) C j ( t , s ) = I 1 ( Γ ( α ) ) s s I α 1 ( t η ) ( h 0 [ d θ U ( η , θ ) ] C j ( η + θ , s ) ) d η + I 1 ( Γ ( α ) ) s t I α 1 ( t η ) ( h 0 [ d θ U ( η , θ ) ] ( C j ( η + θ , s ) C j ( η + θ , s ) ) ) d η .
Let us denote y k j ( t , s , s ) = sup ξ [ a , t ] | c k j ( ξ , s ) c k j ( ξ , s ) | and C ¯ j = sup ξ , s [ a , t ¯ ] C j ( ξ , s ) | ( C ¯ j exists according to Lemma 2) and b = | I 1 ( Γ ( α ) ) | V U . Then for each k , j n and s ( a , s ) from (25) we obtain that
y k j ( t , s , s ) b C ¯ j s s ( t η ) α k 1 d η + b V U s t ( t η ) α k 1 y k j ( η , s , s ) d η = b C ¯ j ( t s ) α k ( t s ) α k α k + b V U s t ( t η ) α k 1 y k j ( η , s , s ) d η
Define for η [ s , t ] the function p ( η , s , s ) = b C ¯ j ( η s ) α k ( η s ) α k α k .
It is simple to see that the function p is monotonically decreasing in η for η [ s , t ] when s ( a , s ) . Furthermore, for arbitrary fixed η [ s , t ] we have that
lim s s p ( η , s , s ) = 0 .
Then from (26) and Theorem 3 it follows that for each j n , k n and s ( a , s ) we obtain the following estimation:
y k j ( t , s , s ) p ( t , s , s ) + s t [ q = 1 ( b V U ( Γ ( α k ) ) q Γ ( α k q ) ( t η ) α k q 1 ] p ( η , s , s ) d η p ( t , s , s ) + p ( s , s , s ) s t [ q = 1 ( b V U ( Γ ( α k ) ) q Γ ( α k q ) ( t η ) α k q 1 ] d η
From (27) and (28) it follows that lim s s y k j ( η , s , s ) = lim s s sup ξ [ a , t ] | c k j ( ξ , s ) c k j ( ξ , s ) | = 0 and hence c k j ( t , s 0 ) exists and c k j ( t , s 0 ) = c k j ( t , s ) for t > s .
(iii.b)
The case when t > s and s ( s , t ) can be treated fully analogically to obtain that c k j ( t , s + 0 ) exists and c k j ( t , s + 0 ) = c k j ( t , s ) .
Let s [ a h , a ] be an arbitrary fixed number and define the following matrix valued function Φ s ( t , s ) = ( φ k j ( t , s ) ) k , j = 1 n : R × [ a h , a ] R n × n with
Φ s ( t , s ) = I , s s t a Θ , t < s o r s < s
and consider following IP:
D a + α T s ( t , s ) = h 0 [ d θ U ( t , θ ) ] T s ( t + θ , s ) , t > a
T s ( t , s ) = Φ s ( t , s ) , t ( , a ] .
Theorem 7.
Let the conditions (S) hold and t ¯ J a be arbitrary.
Then the following statements hold:
1. 
The matrix IP (29) and (30) has a unique absolutely continuous solution T s ( t , s ) in t for t J a for every s [ a h , a ] .
2. 
The matrix function T a h ( t , · ) : [ a , t ¯ ] R n × n is continuous in s for each s [ a , t ¯ ] with s t .
3. 
When s = t with t a , T a h possess first kind jumps and hence is Lebesgue integrable in s on [ a , t ¯ ] .
Proof. 
1.
Let s [ a h , a ] be an arbitrary fixed number. Then since Φ s is PAC for each s [ a h , a ] , then the statement of point 1 follows from Theorem 4.
2.
Let s [ a , t ¯ ] with s t . Then if t < s we have that T a h ( t , s ) = 0 . Consider the case t > s . Then the same way as in the proof of point (iii) of Theorem 6 we obtain that T a h ( t , s ) is continuous in s when s t .
3.
Let s = t with t a . Then obviously lim t s + 0 T a h ( t , s ) = I and lim t s 0 T a h ( t , s ) = 0 and this completes the proof.

5. Applications

We will demonstrate that the obtained results concerning the fundamental matrix allow to improve the integral representation of the solution of the IP (2) and (3) and simplify the proofs.
As usual according the superposition principle we will seek a solution of IP (2) and (3) with initial condition Φ ( t ) 0 , t [ a h , a ] for the case when the function F L 1 l o c ( J a , R n ) is locally bounded.
Let
X F ( t ) = a t K ( t , s ) d s ,
where K ( t , s ) = C ( t , s ) T ( s ) , T ( s ) = R L D a + 1 α F ( s ) .
Theorem 8.
Let the following conditions hold.
1. 
The conditions (S) hold.
2. 
The function F L 1 l o c ( J a , R n ) is locally bounded and D a + α F ( t ) L 1 l o c ( J a , R n ) .
Then the vector function X F ( t ) defined by equality (31) is a solution of IP (2) and (3) with initial condition X F ( t ) = Φ ( t ) 0 , t [ a h , a ] .
Proof. 
From (31) and Theorem 5 it follows that X F ( t ) is an absolutely continuous function in J a . Then we have that
d d s X F ( s ) = d d s a s K ( s , η ) d η = a s C s ( s , η ) T ( η ) d η + C ( t , t ) T ( t ) = a s C s ( s , η ) T ( η ) d η + T ( t ) .
Taking into account (31) and (32), Lemma 2.5 in [1], Condition 2 of the Theorem and applying Fubini’s theorem we obtain
D a + α X F ( t ) = 1 Γ ( 1 α ) a t ( t s ) α ( a s T ( η ) C s ( s , η ) ) d η ) d s + 1 Γ ( 1 α ) a t ( t s ) α T ( s ) d s = a t T ( η ) ( 1 Γ ( 1 α ) ( η t ( t s ) α C s ( s , η ) ) d s ) d η + D a + α 1 T ( t ) = a t T ( η ) ( 1 Γ ( 1 α ) ( a t ( t s ) α C s ( s , η ) ) d s ) d η + D a + α 1 R L D a + 1 α F ( t ) = a t T ( η ) C D a + α C ( t , η ) d η + F ( t )
For the right side of (33) we have
a t C D a + α C ( t , η ) T ( η ) d η = a t T ( η ) [ i = 0 m h 0 [ d θ U i ( t , θ ) ] C ( t + θ , η ) ] d η = i = 0 m h 0 [ d θ U i ( t , θ ) ] ( a t C ( t + θ , η ) T ( η ) d η ) = i = 1 m h 0 [ d θ U i ( t , θ ) ] X F ( t + θ )
and then the statement of the corollary follows from (33) and (34). □
Corollary 3.
Let the following conditions hold.
1. 
The conditions (S) hold.
2. 
The function F L 1 l o c ( J a , R n ) is locally bounded and F ( a ) = 0 .
Then the vector function X F ( t ) defined by equality (31) is a solution of IP (2) and (3) with initial condition X F ( t ) = Φ ( t ) 0 , t [ a h , a ] .
Proof. 
For F L 1 l o c ( J a , R n ) with F ( a ) = 0 is fulfilled
D a + α 1 T ( t ) = D a + α 1 R L D a + 1 α F ( t ) = D a + α 1 C D a + 1 α F ( t ) = F ( t ) F ( a ) = F ( t ) .
Then the proof is the same as in Theorem 8. □
Let T a h ( t , s ) be a solution of the IP (29) and (30) for s [ a h , a ] . For arbitrary function Φ ( t ) B V ( [ a h , a ] , R n ) define the following function:
X Φ ( t ) = a h a T a h ( t , s ) d s Φ ( s ) + T a h ( t , a h ) Φ ( a h ) , t J a
Theorem 9.
Let the following conditions hold.
1. 
The conditions (S) hold.
2. 
The initial function Φ ( t ) B V ( [ a h , a ] , R n ) is not a constant and has finitely many jumps.
Then the function X Φ ( t ) defined with (35) is the unique solution of IP (1) and (3).
Proof. 
Let t [ a h , a ] be an arbitrary fixed number and s [ a h , a ] . Then in virtue of Theorem 6 for every t [ a h , a ] the matrix function T a h ( t , · ) : [ a h , a ] R n × n is continuous in s on [ a h , a ] for s t and when t = s possess first kind jumps. Thus we have that for s [ t , a ] from (30) it follows that T a h ( t , s ) = Θ and for s [ a h , t ] we have T a h ( t , s ) = I . Then X Φ ( t ) = Φ ( t ) and hence X Φ ( t ) satisfies the initial condition (3).
Theorem 7 implies that the matrix valued function T a h ( t , s ) is an absolutely continuous function for t J a and then the vector valued X Φ ( t ) is defined by equality (35) too. Then we have
d d t X Φ ( t ) = d d t a h a T a h ( t , s ) d s Φ ( s ) = a h a T a h ( t , s ) t d s Φ ( s )
and hence applying the Fubini’s theorem we have that
C D a + α X Φ ( t ) = C D a + α a h a T a h ( t , s ) d s Φ ( s ) = a t ( t η ) α ( a h a T a h ( η , s ) η d s Φ ( s ) ) d η = a h a ( a t ( t η ) α T a h ( η , s ) η d s Φ ( s ) ) = a h a C D a + α T a h ( t , s ) d s Φ ( s )
In the right side of (1) for each i m applying the unsymmetric Fubini theorem [28] we obtain that
h 0 [ d θ U ( t , θ ) ] X Φ ( t + θ ) = h 0 [ d θ U ( t , θ ) ] ( a h a T a h ( t + θ , s ) d s Φ ( s ) ) = a h a ( h 0 [ d θ U ( t , θ ) ] ( T a h ( t + θ , s ) d s Φ ( s ) )
From (36) and (37) it follows that
a h a [ C D a + α T a h ( t , s ) h 0 [ d θ U i ( t , θ ) ] T a h ( t + θ , s ) ] d s Φ ( s ) = 0
and hence X Φ ( t ) satisfies (1) for t > a . □
Theorem 10.
Let the following conditions hold.
1. 
The conditions of Theorem 9 hold.
2. 
The function F L 1 l o c ( J a , R n ) is locally bounded.
3. 
Either F ( 0 ) = 0 or D a + α F ( t ) L 1 l o c ( J a , R n ) holds.
Then the function
X f Φ ( t ) = X Φ ( t ) + X f ( t ) = a t C ( t , s ) R L D a + 1 α F ( s ) d s + a h a T a h ( t , s ) d s Φ ( s ) + T a h ( t , a h ) Φ ( a h )
where X F ( t ) and X Φ ( t ) are defined by (31) and (35), respectively, is the unique solution of IP (2) and (3).
Proof. 
The statement of the theorem follows immediately from the superposition principle, Corollary 3 and Theorems 8 and 9. □
Corollary 4.
Let the following conditions hold.
1. 
The conditions of Theorem 10 hold.
2. 
The Lebesgue decomposition of the function Φ ( t ) B V ( [ a h , a ] , R n ) does not include a singular term.
Then the function X F Φ ( t ) defined by (38) possesses the following integral representation:
X f Φ ( t ) = T a h ( t , a ) ( Φ J ( a + ) Φ J ( a ) ) + i T a h ( t , s i ) ( ( Φ J ( s i + ) Φ J ( s i ) ) + a h a T a h ( t , s ) Φ A ( s ) d s + a t C ( t , s ) R L D a + 1 α F ( s ) d s + T a h ( t , a h ) Φ ( a h )
where Φ ( t ) = Φ J ( t ) + Φ A ( t ) and Φ J ( t ) , Φ A ( t ) are the jump term and the absolutely continuous term, respectively, in its Lebesgue decomposition and the summation is taken over all jump points s i S Φ .
Proof. 
The statement of the corollary follows immediately from Theorem 10. □
One of the important questions that arise when we use a fundamental matrix is what kind are its analytical properties concerning the variable s. It is well known that for the integer-order linear differential systems without delay, this problem is fully solved, i.e., both variables are symmetric. Generally speaking, this is not true even in the integer case for the delayed differential systems (the symmetry disappears for the non-autonomous systems and in some cases for the autonomous systems too). As far as we know there are no results in this direction for delayed non-autonomous fractional differential systems. The next result is a first attempt to establish some analytical properties of the fundamental matrix in the mentioned case concerning the variable s.
Theorem 11.
Let the conditions of Theorem 9 hold and the matrix valued function t C ( t , s ) is the fundamental matrix of the system (1).
Then for every fixed a ¯ J a the matrix function C ( t , · ) : [ a , a ¯ ] R n × n for every fixed t ( , s ) ( s , ) is absolutely continuous in s for every compact subinterval s [ a 1 , a 2 ] ( a , a ¯ ) with t [ a 1 , a 2 ] .
Proof. 
Let a ¯ J a be an arbitrary fixed number, s [ a , a ¯ ] be arbitrary and t ( , s ) ( s , ) .
(i)
When t ( , s ) then for the fundamental matrix C ( t , · ) : [ a , a ¯ ] R n × n according to Remark 5 for arbitrary s ( a , a ¯ ) we have that C ( t , s ) = Θ and C ( s , s ) = I . Thus C ( t , s ) has jumps of first kind for t = s and C ( t , s ) is absolutely continuous in s for s ( a , a ¯ ) .
(ii)
Let s ( a , a ¯ ) and t ¯ ( a ¯ , ) be an arbitrary fixed number. For purposes of clarity we assume that s , s ( a , a ¯ ) with s < s . Then as in the proof of Theorem 6 we obtain that (25) holds for t [ a , a ¯ ] and introduce the same notations used there.
Since the function p ( η , s , s ) = b C ¯ j ( η s ) α k ( η s ) α k α k is monotonically decreasing in η for η [ s , a ¯ ] then in a similar way as in the proof of Theorem 6 for each k , j n from (25) we obtain that
y k j ( t , s , s ) b C ¯ j s s ( t η ) α k 1 d η + b V U s t ( t η ) α k 1 y k j ( η , s , s ) d η = b C ¯ j ( t s ) α k ( t s ) α k α k + b V U s t ( t η ) α k 1 y k j ( η , s , s ) d η [ ( t s ) α k ( t s ) α k ] b C ¯ j + C ¯ j b V U α k ( s s ) α k b C ¯ j + C ¯ j b V U α k ( s s ) α m b C ¯ j + C ¯ j b V U α m
when s s 1 . Then we can conclude that c k j ( t , s ) for each k , j n is α m -Hoelder continuous in s in every closed subinterval [ a 1 , a 2 ] ( a , a ¯ ) with a 2 a 1 1 . Thus c k j ( t , s ) is absolutely continuous in s for s ( a , a ¯ ) .
(iii)
Let s ( a , a ¯ ) , t ¯ ( s , a ¯ ] and s < s < t ¯ . Then as in the former case we conclude that (40) holds and then c k j ( t , s ) for each k , j n is α m -Hoelder continuous in s for every closed subinterval [ a 1 , a 2 ] ( a , a ¯ ) for which a 2 a 1 1 . Thus c k j ( t , s ) is absolutely continuous in s for every compact subinterval [ a 1 , a 2 ] ( a , a ¯ ) .
Remark 6.
It is not difficult to see that under the conditions of Theorem 11 the statement of the theorem holds for T a h ( t , s ) too.
Corollary 5.
Let the following conditions hold.
1. 
The conditions of Theorem 9 hold.
2. 
The Lebesgue decomposition of the function Φ ( t ) B V ( [ a h , a ] , R n ) does not include a singular term.
3. 
The delays σ i ( t ) C 1 ( J a , R ¯ + ) , σ i ( t ) < 1 for t J a and i m .
Then the unique solution X Φ ( t ) of IP (1) and (3) defined with (35) has the following representation:
X Φ ( t ) = T a h ( t , a ) Φ J ( a + ) a h a T a h ( t , s ) s Φ A ( s ) d s + T a h ( t , a h ) Φ ( a h )
Proof. 
From conditions 1 and 2 of the theorem it follows that in virtue of Corollary 4 the unique solution X Φ ( t ) of IP (1) and (3) defined with (35) has the representation:
X Φ ( t ) = T a h ( t , a ) ( Φ J ( a + ) Φ J ( a ) ) + i T a h ( t , s i ) ( Φ J ( s i + ) Φ J ( s i ) ) + a h a T a h ( t , s ) Φ A ( s ) d s + T a h ( t , a h ) Φ ( a h )
Let a ¯ J a be an arbitrary fixed number and [ a 1 , a 2 ] ( a , a ¯ ) be an arbitrary compact subinterval. Condition 3 of the theorem implies that for arbitrary fixed number t R the set S t σ = { s [ a 1 , a 2 ] | s = t σ i ( t ) , i m 0 } is finite. Then according Theorem 11 for arbitrary fixed number t R , T a h ( t , s ) is absolutely continuous in s for every compact subinterval [ a 1 , a 2 ] ( a , a ¯ ) . Then integrating by parts the integral in (42) we obtain:
a h a T a h ( t , s ) Φ A ( s ) d s = a h a T a h ( t , s ) d Φ A ( s ) = T a h ( t , a ) Φ J ( a ) + i T a h ( t , s i ) ( Φ J ( s i ) Φ J ( s i + ) ) a h a T a h ( t , s ) s Φ A ( s ) d s
Then the statement of the corollary follows from (42) and (43). □

6. Conclusions

In this article, first the existence and uniqueness of the solutions of an initial problem for linear differential systems with incommensurate order Caputo fractional derivatives and with piecewise absolutely continuous (PAC) initial function is proved.
Then we prove the existence and uniqueness of an absolutely continuous fundamental matrix C ( t , s ) , which has the following properties:
  • C ( t , s ) is absolutely continuous in t for PAC initial functions;
  • C ( t , s ) is absolutely continuous in s (with appropriate additional assumptions).
It must be noted that when the fundamental matrix is absolutely continuous in t and in s, the fundamental matrix has integrable derivatives in t and in s and this allows simpler and more applicable formulas to be obtained in the integral representations, as well as simpler and shorter proofs.
As far we know there are no other articles where such properties of the fundamental matrix C ( t , s ) concerning the variable s for delayed non-autonomous fractional differential systems are obtained. A brief comparison with similar fundamental matrix studies shows that the same system was studied in [15], but there is proof of existence of a continuous fundamental matrix, which is only continuous in t for initial functions with bounded variation. Our result is more general than that obtained in [21] where the smoothness of the fundamental matrix is proven.
Finally, using the properties of the fundamental matrix thus obtained, integral representations are obtained in the paper for the particular solution of the inhomogeneous system with zero initial conditions and for the general solution of the homogeneous system.
A comparison, for example, with the integral representations obtained in [15] shows, that all the proofs are shorter and the obtained formula for the general solution of the homogeneous system is simpler and more applicable.
A general comparison with the analogous results for integer order derivatives shows that those obtained in the article results coincide with them at α = 1 , which means that they are a generalization of the classical ones.
We hope that the results obtained will be useful both for future research and generalizations from a mathematical point of view, as well as for modeling of real-world phenomena.

Author Contributions

Conceptualization, H.K., E.M., M.V. and A.Z. Writing—Review and Editing, H.K., E.M., M.V. and A.Z. All 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.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

In the paper are used the following abbreviations:
ACAbsolutely Continuous
PACPiecewise Absolutely Continuous
BVBounded Variation
IPInitial Problem

References

  1. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science BV: Amsterdam, The Netherlands, 2006. [Google Scholar]
  2. Podlubny, I. Fractional Differential Equation; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  3. Jiao, Z.; Chen, Y.Q.; Podlubny, I. Distributed-Order Dynamic Systems: Stability, Simulation, Applications and Perspectives; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  4. Diethelm, K. The Analysis of Fractional Differential Equations, an Application-Oriented Exposition Using Differential Operators of Caputo Type; Lecture Notes in Mathematics; Springer: Berlin/Heidelberg, Germany, 2010; Volume 2004. [Google Scholar]
  5. 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]
  6. Kaczorek, T.; Idczak, D. Cauchy formula for the time-varying linear systems with Caputo derivative. Fract. Calc. Appl. Anal. 2017, 20, 494–505. [Google Scholar] [CrossRef]
  7. Gomoyunov, M.I. On representation formulas for solutions of linear differential equations with Caputo fractional derivatives. Fract. Calc. Appl. Anal. 2020, 23, 1141–1160. [Google Scholar] [CrossRef]
  8. Idczak, D.; Kamocki, R. On the existence and uniqueness and formula for the solution or R–L fractional Cauchy problem in R n . Fract. Calc. Appl. Anal. 2011, 14, 538–553. [Google Scholar] [CrossRef]
  9. Bonilla, B.; Rivero, M.; Trujillo, J. On systems of linear fractional differential equations with constant coefficients. Appl. Math. Comput. 2007, 187, 68–78. [Google Scholar] [CrossRef]
  10. Odibat, Z.M. Analytic study on linear systems of fractional differential equations. Comput. Math. Appl. 2010, 59, 1171–1183. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, X. Some results of linear fractional order time-delay system. Appl. Math. Comput. 2008, 197, 407–411. [Google Scholar] [CrossRef]
  12. 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]
  13. Krol, K. Asymptotic properties of fractional delay differential equations. Appl. Math. Comput. 2011, 218, 1515–1532. [Google Scholar] [CrossRef] [Green Version]
  14. Wei, Y. The constant variation formulae for singular fractional differential systems with delay. Comput. Math. Appl. 2010, 59, 1184–1190. [Google Scholar] [CrossRef] [Green Version]
  15. Boyadzhiev, D.; Kiskinov, H.; Zahariev, A. Integral representation of solutions of fractional system with distributed delays. Integral Transf. Spec. Funct. 2018, 29. [Google Scholar] [CrossRef]
  16. Zhang, H.; Cao, J.; Jiang, W. General solution of linear fractional neutral differential difference equations. Discret. Dyn. Nat. Soc. 2013, 489521. [Google Scholar] [CrossRef]
  17. 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]
  18. Zahariev, A.; Kiskinov, H. Existence of fundamental matrix for neutral linear fractional system with distributed delays. Int. J. Pure Appl. Math. 2018, 1, 31–51. [Google Scholar]
  19. 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. [Google Scholar] [CrossRef] [Green Version]
  20. Zahariev, A.; Kiskinov, H. Asymptotic stability of the solutions of neutral linear fractional system with nonlinear perturbation. Mathematics 2020, 8, 390. [Google Scholar] [CrossRef] [Green Version]
  21. Zahariev, A.; Kiskinov, H.; Angelova, E. Smoothness of the fundamental matrix of linear fractional system with variable delays. Neural Parall. Sci. Comput. 2019, 27, 71–83. [Google Scholar]
  22. Weissinger, J. Zur Theorie und Anwendung des Iterationsverfahrens. Math. Nachr. 1952, 8, 193–212. [Google Scholar] [CrossRef]
  23. 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]
  24. Diethelm, K.; Ford, N.J. Analysis of fractional differential equations. J. Math. Anal. Appl. 2002, 265, 229–248. [Google Scholar] [CrossRef] [Green Version]
  25. Dieudonne, J. Foundations of Modern Analysis, (Treatise on Analysis 7.6); Academic Press: New York, NY, USA; London, UK, 1969. [Google Scholar]
  26. Ye, H.; Gao, J.; Ding, Y. A generalized Gronwall inequality and its application to a fractional differential equation. J. Math. Anal. Appl. 2007, 328, 1075–1081. [Google Scholar] [CrossRef] [Green Version]
  27. Myshkis, A. Linear Differential Equations with Retarded Argument; Nauka: Moscow, Russia, 1972. (In Russian) [Google Scholar]
  28. Johnson, G.W. An unsymmetric Fubini theorem. Am. Math. Mon. 1984, 91, 131–133. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kiskinov, H.; Madamlieva, E.; Veselinova, M.; Zahariev, A. Existence of Absolutely Continuous Fundamental Matrix of Linear Fractional System with Distributed Delays. Mathematics 2021, 9, 150. https://doi.org/10.3390/math9020150

AMA Style

Kiskinov H, Madamlieva E, Veselinova M, Zahariev A. Existence of Absolutely Continuous Fundamental Matrix of Linear Fractional System with Distributed Delays. Mathematics. 2021; 9(2):150. https://doi.org/10.3390/math9020150

Chicago/Turabian Style

Kiskinov, Hristo, Ekaterina Madamlieva, Magdalena Veselinova, and Andrey Zahariev. 2021. "Existence of Absolutely Continuous Fundamental Matrix of Linear Fractional System with Distributed Delays" Mathematics 9, no. 2: 150. https://doi.org/10.3390/math9020150

APA Style

Kiskinov, H., Madamlieva, E., Veselinova, M., & Zahariev, A. (2021). Existence of Absolutely Continuous Fundamental Matrix of Linear Fractional System with Distributed Delays. Mathematics, 9(2), 150. https://doi.org/10.3390/math9020150

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