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
, which is absolutely continuous in
t on every compact subinterval of
. 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
be an arbitrary number and denote by
the linear space of all locally Lebesgue integrable functions
. Then for each
, each
and
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]):
If
then the next formula gives a direct definition of the Caputo left side derivative:
Everywhere below, the following notations will be used: , , denote the identity and zero matrix, respectively, denotes the k-th column of the identity matrix and is the zero element.
For we will use the notations , for we denote and for simplicity we will use the notation for the left side Caputo fractional derivative.
Consider the homogeneous linear delayed system of incommensurate type and distributed delay in the following general form
or described in more detailed form
and the corresponding nonhomogeneous one
where
.
Definition 1. A function is called piecewise absolutely continuous on some interval (denoted ), if it is with bounded variation (BV) in t on J, there is no singular term in the Lebesgue decomposition of and the set of discontinuity points of the function has no limit points in J.
Definition 2. With we denote the Banach space of right continuous vector functions with norm and the subspace of all absolutely continuous functions by , i.e., .
The set of jump points of every initial vector function we denote by . We emphasize that the set is finite for every compact interval and the case is not excluded.
For the system (
1) or (
2) introduce the following initial conditions:
We say that for the kernel the conditions (S) hold, if the following conditions are fulfilled:
- (S1)
The functions are measurable in and normalized so that for when and for all , and for .
- (S2)
The Lebesgue decomposition of the kernel
for
and
has the form:
where
,
are locally bounded on
,
is the Heaviside function, the delays
are bounded with
,
and
are continuous from left in
on
,
.
are locally bounded on and .
- (S3)
There exists a locally bounded function such that for and for every the following relations hold: .
- (S4)
The sets for every do not have limit points.
Consider the following auxiliary system in matrix form
where
, or for
in more detailed form
with the initial condition (
3).
Definition 3. The vector function is a solution of the IP (
2)
and (
3)
or IP (
3)
and (
4)
in if satisfies the system (
2)
, respectively, (
4)
for all and the initial condition (
3)
for each .
In virtue of Lemma 3.3 in [
15] every solution
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 and let the following conditions hold:
- 1.
There exists a sequence , with .
- 2.
The operator satisfies for each and for arbitrary the inequality
Then T has a uniquely fixed point and for every we have that .
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 be an arbitrary real Banach space.
Definition 4 ([25]).
The function is called a regulated function if it has one-sided (left and right) limits at every point .
Remark 2. If , where is an arbitrary compact interval and is a regulated function, then it is assumed that in the left (right) end on the interval the function has only a right (left) limit.
Theorem 2 ([25]). Let be an arbitrary function.
Then a necessary and sufficient condition for to be a regulated function is that must in every compact interval be a limit of n uniformly convergent sequence of step-functions (i.e., with respect to the supremum norm ).
Theorem 3 ([26]). Let and be arbitrary fixed numbers and the following conditions hold:
- 1.
The functions for some .
- 2.
The function for some and is nondecreasing.
- 3.
For every is fulfilled .
Then for the following inequality holds Remark 3. Note that the statement of Theorem 3 is proved in the partial case , but with small modifications the proof will be correct for arbitrary .
3. Existence and Uniqueness of the Solutions
Let it be that for every
consider the corresponding linear space
For each
define the set
where
is arbitrary and define a metric function
with
for each
.
Since
then for every
, according to a well-known result we conclude that
is a complete metric space concerning the metric
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 the set is a complete metric space concerning the metric too.
Proof. Let be an arbitrary fixed number and consider an arbitrary Cauchy sequence , i.e., . It is clear that there exists a vector valued function such that .
Let
and
be arbitrary numbers. There exist
and
such that for every
we have that
for each
. Since
then there exists
, such that for every finite sequence of pairwise disjoint subintervals
,
with
the inequality
holds. Then whenever when
we have that
Thus
for every
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 is arbitrary.
Then the IP (
1)
and (
3)
has a unique absolutely continuous solution in the interval for every . Proof. In the proof of this theorem we will use the approach introduced in [
23].
Let
and
be arbitrary. From (
5) in the case
for
we obtain the system
For every define the operator
via the operators
for
by
and the following additional condition:
Since
is PAC, the conditions (S) imply that for each
, the functions
are at least piecewise continuous on the interval
(see [
27], Lemma 1). Then the right side of (
7) is absolutely continuous on the interval
, which implies that the function
is absolutely continuous on the same interval and since
for
, then
. Thus the operator
ℜ maps
into
.
Since according to Lemma 1 the space is a full metric space concerning the metric , 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 .
We recall that the , has a local minimum at , where it attains the value . There exists such that and for every we denote with that number for which .
Let denote
and let
be arbitrary. Then from (
7) for
and every
we obtain
Let assume that for some
and for every
the inequality
holds. Obviously from (
9) it follows that the inequality (
10) holds for each
and every
at least for
. Denoting for simplicity
we obtain
Let us assume that (
10) holds for each
, every
and for some
. Then a similar way as in (
9) from (
10) and (
11) we obtain
Substituting
in the integral in the right side of (
12) and using the well-known relation between the gamma and beta functions we obtain
and hence (
10) holds for every
, for each
and every
.
For
from (
10) and (
13) it follows that
and denote
for
.
For all
from (
10) and (
13) we obtain that
and denote
for
.
Consider the one parameter Mittag–Leffler function
. It is simply to see that the series
is convergent because it is the considered Mittag–Leffler function evaluated at
. Then we have that
and then from Theorem 1 it follows that the IP (
1) and (
3) has a unique solution in
. □
Corollary 1. Let the conditions of Theorem 4 hold.
Then the IP (
1)
and (
3)
has a unique absolutely continuous solution in the interval . Proof. Let
be arbitrary and assume the contrary, that there exists
such that the solution
in the interval
does not possess a prolongation. Let
be an arbitrary number. In virtue of Theorem 4 the IP (
1) and (
3) has a unique absolutely continuous solution
in the interval
. The solution
obviously is a prolongation of the solution
, which contradicts of our assumption that
. □
For arbitrary fixed
consider the following auxiliary system
with the following condition
Corollary 2. Let the following conditions hold:
- 1.
Conditions (S) hold.
- 2.
The function has the form
Then for each and arbitrary the problem (
16)
and (
17)
has a unique solution, which satisfies Equation (
16)
for , the condition (
17)
for and is absolutely continuous in with a first kind jump at . Proof. Let
be an arbitrary fixed number, introduce the system
and consider the IP (
17) and (
18). Since Theorem 4 is proved for arbitrary
and
for arbitrary
is PAC on the interval
, then from Theorem 4 it follows that the IP (
17) and (
18) possess a unique solution
and moreover, from (
17) it follows also that
. Then for every
we have
and hence
satisfies the Equation (
16) for
and the condition (
17) for
.
Let consider an IP with Equation (
16) for
and initial condition
for
. Then obviously
is its unique solution in
. This completes the proof. □
4. Fundamental Matrix
Let
be an arbitrary fixed number and define the following matrix valued function
with
and denote
.
For arbitrary fixed number
consider the following matrix IP
Definition 5. The matrix valued function , is called a solution of the IP (
19)
and (
20)
in if is continuous for and satisfies the matrix Equation (
19)
on as well as the initial condition (
20)
too. Remark 5. Practically in condition (
20)
we need only the values of for , but for convenience we define also for . Then is prolonged as continuous in t function on . Theorem 5. Let the conditions (S) hold.
Then for every initial point , the matrix IP (
19)
and (
20)
has a unique absolutely continuous solution in the interval . Proof. The statement of the Theorem follows immediately from Corollary 2. □
Definition 6. The matrix , 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 is the fundamental matrix of the system (
1)
. Then for every the matrix function is locally bounded in s for and .
Proof. Let be an arbitrary fixed number, be arbitrary and consider the fundamental matrix . According to Remark 5 for we have that for and .
Taking into account Theorem 5 and Corollary 2 it is easy to be seen, that the unique solution
of IP (
19) and (
20) is a solution of the equation
Introduce for every
the notations:
(since
for
),
and then from (
21) it follows
and hence for each
From Theorem 3 and (
22) for
we obtain the estimation
and thus
is locally bounded in
s for
and
. □
Theorem 6. Let the conditions (S) hold and the matrix valued function be the fundamental matrix of system (
1)
. Then for every fixed the matrix function is continuous for when , for possess first kind jumps and hence is Lebesgue integrable in s on for each .
Proof. Let be an arbitrary fixed number, be arbitrary and consider the fundamental matrix . According to Remark 5 for we have that for and . In virtue of Theorem 2 and Lemma 2 it is enough to prove that has left and right limits for each and there exist the limits and .
- (i)
Let be arbitrary and let . Then it is simply to see that for every . Note that for we have that holds.
- (ii)
Let be arbitrary and let . Then and since , then we can conclude that exists and hence has jumps of first kind on the line for each .
- (iii)
Let . Then we have to consider two cases: either (left limit in ) or with (right limit in ).
- (iii.a)
For purposes of clarity we assume that
. According to Corollary 2 for each
, the IP (
16) and (
17) has unique solutions
and
for the initial functions
and
, respectively.
Taking into account that
and subtracting both sides of (
24) from the corresponding sides of (
23) we obtain
Let us denote
and
(
exists according to Lemma 2) and
. Then for each
and
from (
25) we obtain that
Define for the function .
It is simple to see that the function
p is monotonically decreasing in
for
when
. Furthermore, for arbitrary fixed
we have that
Then from (
26) and Theorem 3 it follows that for each
and
we obtain the following estimation:
From (
27) and (
28) it follows that
and hence
exists and
for
.
- (iii.b)
The case when and can be treated fully analogically to obtain that exists and .
□
Let
be an arbitrary fixed number and define the following matrix valued function
with
and consider following IP:
Theorem 7. Let the conditions (S) hold and be arbitrary.
Then the following statements hold:
- 1.
The matrix IP (
29)
and (
30)
has a unique absolutely continuous solution in t for for every . - 2.
The matrix function is continuous in s for each with .
- 3.
When with , possess first kind jumps and hence is Lebesgue integrable in s on .
Proof. - 1.
Let be an arbitrary fixed number. Then since is PAC for each , then the statement of point 1 follows from Theorem 4.
- 2.
Let with . Then if we have that . Consider the case . Then the same way as in the proof of point (iii) of Theorem 6 we obtain that is continuous in s when .
- 3.
Let with . Then obviously and 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
for the case when the function
is locally bounded.
Let
where
.
Theorem 8. Let the following conditions hold.
- 1.
The conditions (S) hold.
- 2.
The function is locally bounded and .
Then the vector function defined by equality (
31)
is a solution of IP (
2)
and (
3)
with initial condition . Proof. From (
31) and Theorem 5 it follows that
is an absolutely continuous function in
. Then we have that
Taking into account (
31) and (
32), Lemma 2.5 in [
1], Condition 2 of the Theorem and applying Fubini’s theorem we obtain
For the right side of (
33) we have
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 is locally bounded and .
Then the vector function defined by equality (
31)
is a solution of IP (
2)
and (
3)
with initial condition . Proof. For
with
is fulfilled
Then the proof is the same as in Theorem 8. □
Let
be a solution of the IP (
29) and (
30) for
. For arbitrary function
define the following function:
Theorem 9. Let the following conditions hold.
- 1.
The conditions (S) hold.
- 2.
The initial function is not a constant and has finitely many jumps.
Then the function defined with (
35)
is the unique solution of IP (
1)
and (
3).
Proof. Let
be an arbitrary fixed number and
. Then in virtue of Theorem 6 for every
the matrix function
is continuous in
s on
for
and when
possess first kind jumps. Thus we have that for
from (
30) it follows that
and for
we have
. Then
and hence
satisfies the initial condition (
3).
Theorem 7 implies that the matrix valued function
is an absolutely continuous function for
and then the vector valued
is defined by equality (
35) too. Then we have
and hence applying the Fubini’s theorem we have that
In the right side of (
1) for each
applying the unsymmetric Fubini theorem [
28] we obtain that
From (
36) and (
37) it follows that
and hence
satisfies (
1) for
. □
Theorem 10. Let the following conditions hold.
- 1.
The conditions of Theorem 9 hold.
- 2.
The function is locally bounded.
- 3.
Either or holds.
Then the functionwhere and 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 does not include a singular term.
Then the function defined by (
38)
possesses the following integral representation:where and are the jump term and the absolutely continuous term, respectively, in its Lebesgue decomposition and the summation is taken over all jump points . 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 is the fundamental matrix of the system (
1)
. Then for every fixed the matrix function for every fixed is absolutely continuous in s for every compact subinterval with .
Proof. Let be an arbitrary fixed number, be arbitrary and .
- (i)
When then for the fundamental matrix according to Remark 5 for arbitrary we have that and . Thus has jumps of first kind for and is absolutely continuous in s for .
- (ii)
Let
and
be an arbitrary fixed number. For purposes of clarity we assume that
with
. Then as in the proof of Theorem 6 we obtain that (
25) holds for
and introduce the same notations used there.
Since the function
is monotonically decreasing in
for
then in a similar way as in the proof of Theorem 6 for each
from (
25) we obtain that
when
. Then we can conclude that
for each
is
-Hoelder continuous in
s in every closed subinterval
with
. Thus
is absolutely continuous in
s for
.
- (iii)
Let
and
. Then as in the former case we conclude that (
40) holds and then
for each
is
-Hoelder continuous in
s for every closed subinterval
for which
. Thus
is absolutely continuous in
s for every compact subinterval
.
□
Remark 6. It is not difficult to see that under the conditions of Theorem 11 the statement of the theorem holds for too.
Corollary 5. Let the following conditions hold.
- 1.
The conditions of Theorem 9 hold.
- 2.
The Lebesgue decomposition of the function does not include a singular term.
- 3.
The delays for and .
Then the unique solution of IP (
1)
and (
3)
defined with (
35)
has the following representation: Proof. From conditions 1 and 2 of the theorem it follows that in virtue of Corollary 4 the unique solution
of IP (
1) and (
3) defined with (
35) has the representation:
Let
be an arbitrary fixed number and
be an arbitrary compact subinterval. Condition 3 of the theorem implies that for arbitrary fixed number
the set
is finite. Then according Theorem 11 for arbitrary fixed number
,
is absolutely continuous in
s for every compact subinterval
. Then integrating by parts the integral in (
42) we obtain:
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 , which has the following properties:
is absolutely continuous in t for PAC initial functions;
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
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 , 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.