1. Introduction
Practically, it has been shown that many real-world phenomena in various fields of science can be represented more accurately through mathematical models, including fractional differential equations [
1,
2,
3]. More detailed information on fractional calculus theory and fractional differential equations can be seen in the monographs of Kilbas et al. [
4] and Podlubny [
5]. Compared to fractional equations with Caputo-type derivatives, the fractional equations with Riemann–Liouville-type derivatives have been studied significantly less. The first obstacle was the lack of meaningful geometric and physical interpretations of the Riemann–Liouville type of integration and differentiation, which was overcome with the appearance of Ref. [
6]. The other problem was to derive an appropriate formulation of the initial conditions for the initial problem. Some considerations about them are given in [
7,
8], but as was mentioned in [
9] “The initial value problem is a subject that remains quite up-to-date”.
It is well known that the existence of explicit solutions or an integral representation (variation of constants formula) of the solutions of linear fractional differential equations and/or systems (ordinary or delayed) is a main tool in executing their qualitative analysis. That is why establishing either explicit solutions (see [
10]) or integral representation, for which existence of a fundamental matrix is needed (see [
11]), are important tasks for stability analysis, especially in the case of equations with Riemann–Liouville derivatives. But it is surprising that there are not many articles devoted to this problem. As far as we know, a survey concerning stability results for retarded and neutral fractional differential equations with Riemann–Liouville-type derivatives does not exist. In general, partially this gap can be fulfilled by reading the overview [
12] and the references therein. From the recent works devoted to the discussed themes concerning fractional differential equations with Riemann–Liouville-type derivatives we refer to [
13] for equations, to [
14,
15,
16] where retarded fractional differential equations are considered, and to [
17,
18] where the neutral case is considered. We suppose that the mentioned works with their references give a good enough picture of the studies in this area.
In the present article, we study a class of nonlinear disturbed neutral linear fractional systems with derivatives in the the Riemann–Liouville sense and distributed delay. We first study the important problem of the existence and uniqueness of the solutions of an initial problem (IP) for these systems in the case of discontinuous initial functions. As far as we know, there are no results devoted to the initial problem with discontinuous initial functions for neutral differential equations with derivatives in the Riemann–Liouville sense. Since the classically stability concepts are not directly applicable to systems with derivatives with the property that the derivative of a constant is not equal to zero (like the Riemann–Liouville fractional derivative, for example), we introduce a new concept for weighted stabilities in the Lyapunov, Ulam–Hyers, and Ulam–Hyers–Rassias senses, which coincides with the classical stability concepts for the cases of integer-order or Caputo-type derivatives.
The following abbreviations will be used in this manuscript: BV—bounded variation; GAS—globally asymptotically stable; UH—Ulam–Hyers; LAS—locally asymptotically stable; LT—Laplace transform; UHR—Ulam–Hyers–Rassias; IP—initial problem; PC—piecewise continuous; RL—Riemann–Liouville; WML—weighted Mittag–Leffler.
The paper is organized as follows:
Section 2 presents the necessary definitions and properties concerning the RL and Caputo fractional derivatives, the problem statement, and the needed auxiliary definitions and facts for our exposition.
Section 3 is devoted to the existence and the uniqueness of the solutions of the initial problem (IP) for the studied class of nonlinear neutral systems in the case when the initial function is discontinuous. In
Section 4, as a consequence of the derived result we prove the existence and uniqueness of a fundamental matrix and a generalized fundamental matrix for the linear homogeneous system, as well as establishing an integral representation of the solutions of the IP for the homogeneous system and the corresponding inhomogeneous system.
Section 5 is devoted to a new concept for weighted stabilities in the Lyapunov, UH, and UHR senses, the definition of which coincides with the classical stability concepts for the cases of integer-order or Caputo-type derivatives. It is proved that the zero solution of the homogeneous system is weighted stable if and only if all its solutions are weighted bounded. In addition, it is established that the weighted stabilities in the Lyapunov and UH senses are equivalent if and only if the inequality appearing in the UH definition possesses only bounded solutions. In
Section 6, natural sufficient conditions are obtained under which the property of weighted global asymptotic stability of the zero solution of the homogeneous system is preserved under nonlinear disturbances. Finally, in
Section 7 an illustrative example is presented. In the last
Section 8, we summarize some conclusions concerning the obtained results and propose some open problems.
2. Preliminaries and Problem Statement
As is standard to avoid possible misunderstandings, we recall some properties concerning the RL-type derivatives, as well as the needed definitions, conditions, and auxiliary results necessary for the exposition below. For more comprehensive information on fractional calculus we refer to Refs. [
4,
5].
For a function , we will say that some property holds locally if it is fulfilled on every compact subinterval . We will use the following notation for the real linear spaces: consists of all locally Lebesgue integrable functions, and are, respectively, the subspaces of all functions which are locally bounded and of all functions which have locally bounded variation.
We define the left-sided fractional RL integral operator of arbitrary order
with lower limit
for any
, via the relation
The corresponding left-hand side RL and Caputo fractional derivatives of arbitrary order
with lower limit
a are defined for any
by
For shortness, we will use the notation
.
Consider the neutral fractional nonlinear system with RL-type derivatives of order
with lower limit
and distributed delays in the following general form:
where
,
,
,
,
,
,
,
,
, (the notation
means column and superscript
T denotes the transposed vector),
,
,
, and
(Krasovskii-type functional) for
and
. In addition,
denote the identity and the zero matrices,
denotes the zero vector, and
,
.
For
,
, we use the norm
and
denotes the set of all jump points for any
. We will use the following Banach spaces of initial functions:
endowed with the norm
. In addition, we assume that in all spaces the functions
are right continuous at
(for
this assumption ultimately holds).
For any
we define the initial condition for the system (
1) as follows:
More details concerning other types of initial conditions are given in [
8] for the case when the initial functions are continuous.
Let us introduce the auxiliary integral system (whose system is considered in detail in Lemma 4 below) for any
and
:
Let and be arbitrary.
Definition 1 ([
4])
. The function , will be called right γ-continuous at a if the function i.e., the function is right continuous at a. By we denote the real linear space of all right γ-continuous-at-a functions and for any by the real Banach spacewith norm . Definition 2 ([
11])
. The vector function is a solution of IP (1), (2) in if satisfies the system (1) for all and the initial condition (2) too. Definition 3 ([
11])
. The vector function is a solution of IP (3), (2) in if satisfies system (3) for all and initial condition (2) too. The hypotheses
stated below in Definition 4, as in the cases of systems with derivatives of integer order or Caputo-type fractional order, will play a major role in the solvability of IP (
1), (
2) (IP (
3), (
2)) (see [
19,
20,
21]).
Definition 4 ([
19,
20,
21])
. We say that for the kernels the hypotheses hold if the following conditions are fulfilled for any and , :- (1)
The functions , are measurable in and normalized so that , for , for ,
for , , for any .
- (2)
For any fixed the kernels and are left continuous in θ on and , and , is uniformly bounded in and is uniformly nonatomic at zero [21] ( i.e., for every , there exists such that for any , we have that ). - (3)
For any fixed , the Lebesgue decomposition of the kernels and has the form: , , , and is the Heaviside function.
- (4)
The sets , do not have limit points and the relations hold for any .
Definition 5 ([
20,
21])
. We say that the vector-valued functional satisfies the modified Caratheodory conditions if the following conditions hold:- (1)
For almost all fixed the functional is continuous in arbitrary and for each fixed function the function .
- (2)
(Local Lipschitz-type condition) For any and for some its vicinity , there exists a function such that the inequalities hold for every .
In our investigations below, we will use the following auxiliary result:
Lemma 1 ([
4])
. Let and let be a Lebesgue measurable function on J.- (a)
If there exists a.e. (almost everywhere) the limit , then there also exists a.e. the limit .
- (b)
If there exist a.e. the limit and , then we have that .
Definition 6 ([
22])
. The low terminal a will be called a noncritical point (noncritical jump point) for some initial function relative to the delay , if the equality implies that there exists a constant (eventually depending on ) such that for . Definition 7 ([
22])
. The low terminal a for a function with will be called a critical point (critical jump point) relative to some delay if the equality implies that there exists a constant (eventually depending on ) such that for . Lemma 2 ([
22])
. Let the hypotheses hold and . Then, there exists (eventually depending from ), such that is continuous for . Lemma 3 ([
22])
. Let the hypotheses hold and .Then, for any initial function with , one of the following statements holds:
- (i)
The statement of Lemma 2 holds.
- (ii)
The low terminal a is a critical jump point for Φ relative to the kernel for some .
Theorem 1 ([
23], Krasnosel’skii’s fixed point theorem)
. Let be a Banach space, be a nonempty, closed, and convex subset of and the maps satisfy the following conditions:- (i)
The operator is a contraction with constant ;
- (ii)
The operator is continuous and the set is contained in a compact set;
- (iii)
For any , we have that .
Then, there exists a with .
Theorem 2 ([
24], Corollary 2)
. Suppose that and the following conditions hold:- 1.
The functions .
- 2.
The functions are nondecreasing, is bounded on J, and .
- 3.
The inequality holds for .
Then, for any we have that .
Remark 1. The main difference between the systems (even in the linear case) with different types of fractional derivatives from the point of view of their applicability as model tools of real-world phenomena is the answer for a concrete type of derivative, whether the derivative of a constant is identically equal to zero or not. Typical representatives of these two classes are the fractional derivatives of Caputo and RL types. So, for the study of retarded and neutral systems with Caputo derivatives a lot of the ideas and techniques known from the systems with integer-order derivatives can be used, because of the continuity of the corresponding Krasovskii functional. Note that in some important special cases it is possible that the Krasovskii functional can have a jump only of the first kind at the initial point, but this obstacle can be overcome with the help of some variants of Banach’s contraction principle. In contrast, for retarded and neutral systems with RL derivatives we have at least in the technical aspect several complications, mainly based on the availability of a discontinuity (jump of second kind) of the solutions at the low terminal, a fact which greatly complicates the use of the Riesz theorem for the representation of linear continuous functionals (Krasovskii functionals) on via a Lebesgue–Stieltjes integral. Even for the retarded systems with RL derivatives, more sophisticated techniques such as Weissinger theorem combined with application of the Mittag–Lafleur function must be used. We emphasize that the important neutral case with RL derivatives, the study of which is the aim of this article, is essentially more complicated even in comparison with the retarded systems with the same kind of derivatives, not only in the technical aspects, but in their applications as models.
3. The Initial Problem (1), (2) with Discontinuous Initial Functions
The next Lemma is well known for the case of systems with Caputo-type derivatives and for delayed systems with RL-type derivatives, and therefore, we will sketch only the differences appearing in the neutral case with RL-type derivatives and distributed delays.
Lemma 4. Let the following conditions be fulfilled:
- 1.
The hypotheses hold.
- 2.
For any , the vector-valued functional .
Then, every solution of IP (1), (2) is a solution of IP (3), (2) and vice versa. Proof. Let
be arbitrary and
be the corresponding solution of IP (
1), (
2). Then, according to formula 2.1.40 in [
5], we have
Then, obviously
will be a solution of IP (
3), (
2) if we prove
Condition (
) implies that there exists a constant
such that
for
, and since
is uniformly nonatomic at zero in
t, then for any
there exists
such that
uniformly for
. Then, for
we have
and, hence, for the first addend on the right-hand side of (
4) the relation
holds. Since
, then there exists
such that for
and any
we have that
Then, for the second addend on the right-hand side of (
4) for
we obtain that
which implies that
Thus,
is a solution of IP (
3), (
2).
Let
is a solution of IP (
3), (
2), and then, applying the operator
to both sides of (
3) we obtain immediately the opposite statement. □
Modifying the approach in [
11], for an arbitrary fixed number
we introduce the following real linear space:
and for any
introduce the linear subspaces:
endowed with the norm
.
Let
with
be arbitrary fixed and define the nonempty, closed, and convex subset
as follows:
It is simple to see that for any
we have that the norms
and
are equivalent.
For any
and
define the operator
as follows:
with the additional conditions
In addition, for arbitrary
,
and
we define the operators
as follows:
Then, the operator
, defined via (
5), (
6), for any
has the form
The next technical lemma is an immediate generalization of Theorem 4 in [
25]. The lemma below is more appropriate for applications and can be used in the cases of RL derivatives too. For simplicity we will consider only the case when the initial function
has only one jump point, but it is simple to see that the proof provided below can be used in the case of an arbitrary but finite number of jumps too.
Lemma 5. Let the following conditions be fulfilled:
- 1.
The hypotheses hold and with are arbitrary.
- 2.
The low terminal a is a jump point for Φ relative to the delay for some .
Then,
Proof. We will consider the more complicated case when
and
a is a critical jump point. Then, since
for all
we have that
. Let assume that there exists a sequence
such that
, and hence, the set
has a limit point which contradicts condition
. Then, there exists
such that
for any
. Let us assume that
,
and
is the number existing according to condition
such that
for any
. Then, for any
with
which is impossible. Thus, we have that
, and hence,
.
The proofs of the cases when a is a noncritical jump point and when are simpler versions of the proof in the case when (a is a critical jump point), and therefore, will be omitted. □
It is clear that without loss of generality it is possible to renumber all concentrated delays in the jump part and give the first numbers, of these delays, relative to which a is a critical jump point. Everywhere below we will assume that this renumbering is made.
Theorem 3. Let the following conditions be fulfilled:
- 1.
The conditions of Lemma 4 hold.
- 2.
where denote the partial derivatives with respect to θ.
- 3.
The condition holds.
Then, for any with there exists such that IP (3), (2) has at least one solution with an interval of existence . Proof. As above, we will consider the more complicated case when with . Let is arbitrary. As a first step we will prove that .
For any
from (
6) and Lemmas 2 and 3 it follows that
where
, and
.
The first addend on the right-hand side of (
10) is continuous for any
, and hence, belongs to
. Since
we can prolong
as a function from
. According to Lemma 1 in [
11] and Lemma 4 in [
22], since the initial function
is a continuous function for
, the second and the third addends on the right-hand side of (
10) are continuous functions at
. Taking into account that
a is a critical jump point relative to the delays with numbers
, we conclude that the fourth addend is a continuous function at
too. For the fifth addend, via the substitution
we obtain
Because
and
for any
are continuous at
, then the fifth addend is a continuous function at
by virtue of lemma 1 in [
19]. Since the right-hand side of (
10) is a continuous function at
, then the second additional relation in (
7) holds too, and hence, from (
7) it follows that
. Thus,
.
As a second step we will prove that the operator
is a contraction. Denote
,
and
Then, for arbitrary
from (
7) and Lemmas 2 and 3 it follows that
For the first addend on the right-hand side of (
12), since
and the low terminal
a is a critical jump point for
relative to the delays
, by virtue of Lemma 5 there exists
such that for any
we have the estimation
Hypothesis
implies that for
there exists
such that
for any
, where
, we obtain
Then, the second addend in (
12) for any
has the estimation
and, hence, for
, where
we obtain that
Multiplying both sides of (
12) by
and then from (
12), (
13), and (
14) it follows that for any
we obtain that
and, hence,
is a contraction operator in
.
For arbitrary
, from (
7) we obtain that
For the first addend on the right-hand side of (
15), for any
the following estimation holds:
and analogously, for the second one we obtain
Multiplying both sides of (
15) by
, and then, from (
15), (
16), and (
17), it follows that for any
we obtain that
and, hence,
is a continuous operator.
Let
and
be arbitrary fixed and consider the ball
For any
, from (
18) it follows that
and, hence, the set
is uniformly bounded.
To apply Theorem 1 we must prove that the set
is equicontinuous, and hence, it is at least relative compact. For any
and arbitrary
we have the estimation
Since for
from conditions
and (
19) we have that
then from (
20) and (
21) for any
we obtain
Then, since
is uniformly continuous at
, for any
there exists
such that if
we have that
and, hence, from (
22) it follows that the set
is equicontinuous. Thus, by virtue of the Arzella–Ascoli theorem,
is at least relative compact. For arbitrary
we consider the sum
. Taking into account that
and can be prolonged as
for any
, from (
7) it follows that
. Since
, then
too. Considering the fact that
then, we obtain that
and, hence,
. Applying Theorem 1, we obtain that there exists at least one fixed point
, i.e.,
,
. □
Theorem 4. Let the conditions of Theorem 3 be fulfilled.
Then, for any with the solution of IP (3), (2) is unique in the interval of its existence . Proof. Assume the contrary, that IP (
3), (
2) possess two different solutions
. Then, from (
3) for any
we obtain that
Since (
2) implies that
, then from (
23) it follows that
Denoting
and applying Theorem 2, we obtain for
the estimation
which contradicts our assumption. The proof of the case when
is almost the same as the proof in the case when
, and therefore, will be omitted. □
Definition 8. We say that the solution of IP (3), (2) is a continuation of a solution of IP (3), (2) if and on . The solution of IP (3), (2) which coincides with all of its continuations will be called the maximal solution. Theorem 5. Let the conditions of Theorem 3 be fulfilled.
Then, for any with the unique solution with interval of existence of IP (3), (2) can be continued as a unique solution of IP (3), (2) with interval of existence . Proof. Assume the contrary, that there exists a maximal solution
of IP (
3), (
2) with the interval of its existence
(which is closed from the right),
, and define the operator
with
via (
5) with the additional conditions
;
. Then, we can prove, fully analogous to the proof of Theorem 3, that for any
with
there exists
such that IP (
3), (
2) has at least one solution
with interval of existence
, and obviously,
for any
. Thus, we obtain that
is a continuation of
, which is impossible.
Let assume that that there exists a maximal solution
of IP (
3), (
2) and its interval of existence is open and finite, i.e.,
. Then, since
satisfies (
3) for any
and the additional condition
too, we conclude that
satisfies (
3) for any
, which contradicts our assumption. Thus, we have that
, which completes the proof. □
Remark 2. We note that from Theorems 3–5 it follows that the requirement in Theorem 2 [11], that in the Lebesgue decomposition of a singular part does not exist, is unnecessary. So, the results proved in Theorems 3–5 generalize the statement of Theorem 2 [11] even in the retarded case. 4. The Linear Case – Fundamental Matrices and Integral Representation of the Solutions
In this section, we establish the existence of two different types of fundamental matrices and of their bases; we obtain two types integral representations for the solutions of IP (
24), (
2) for different kinds of initial functions. In the case when
, system (
1) becomes an inhomogeneous linear system and has the following form:
where
and
. Consider the corresponding homogeneous linear system of (
1) and (
24):
and following [
11] introduce for any
and
the initial function
where
denotes the
j-th column of the identity matrix
.
From Theorem 5 it follows that for any
and
IP (
25), (
2) has a unique solution
, which satisfies (
25) for any
and the initial condition (
2) with the initial function (
26) too. Then, the matrix
will be called fundamental matrix of the system (
25).
Additionally, we introduce the following initial function for any
and
:
As above, for any
and
we can conclude that IP (
24), (
2) has a unique solution
which satisfies (
24) for any
and the initial condition (
2) with the initial function (
27) too. The matrix
we will call the generalized fundamental matrix.
Both matrices play crucial roles in the construction of different kinds of integral representations of the solutions of IP (
24), (
2) and IP (
25), (
2), having a lot of applications in the qualitative theory of the fractional systems.
It is clear that
if and only if
. First, we will establish some analytical properties of the matrices
and
, needed to obtain appropriate integral representations of the solutions of (
24) and (
25).
As in the case of Caputo-type derivatives [
26], for any
we define the function
where
, and
. It is necessary to establish for any
that the solutions
to IP (
25), (
2) are Lebesgue integrable in
s on
for any fixed
.
Remark 3. Note that to prove the statement that are Lebesgue integrable in s on for any fixed is not trivial, in contrast with the Caputo case where this fact can be proved via elementary application of a generalized Bellman–Gronwall inequality (see Theorem 2). The main tool in all these proofs of theorems devoted to integral representations of the solutions of fractional systems with Caputo-type derivatives is the Fubini theorem. Thus, from our point of view it is important to justify its correct application. The greatest obstacle in the neutral case with RL derivatives is the problem with the integrability of the integrands relative to the product measures, which are defined as the products of two Lebesgue–Stieltjes measures. Note that this problem arises in our case, since the integrand is with a singularity of order at the point a in contrast to the case of neutral systems with Caputo-type fractional derivatives, in which case the integrand is a Lebesgue measurable and locally bounded function. This difference leads at least to some additional technical complications, but not only that.
To overcome this obstacle, we need the statements of the next two lemmas.
Lemma 6. Let the following conditions be fulfilled:
- 1.
Conditions 1 and 2 of Theorem 3 hold.
- 2.
Then, for any compact interval , and we have that Proof. Let
and
be arbitrary fixed. Then, from (
25), since
we have
and, hence, by virtue of Theorem 2 we obtain
where in (
29)
is the one-parameter Mittag–Laffler function. This completes the proof. □
By virtue of Lemma 6 we can state that the integral on the right-hand side of (
28) is correctly defined.
Corollary 1. Let the conditions of Lemma 6 hold.
Then, for arbitrary fixed , is continuous in t for or when for . At , has a first-kind jump.
Lemma 7. Let the conditions of Lemma 6 hold.
Then, for any we have that .
Proof. Let
with
,
. Then, integrating by parts formally the integral in (
28) we have that
and since the integral on the left-hand side of (
30) exists, then the integral on the right-hand side exists too (see [
27], point 5, page 229). Let
be arbitrary,
and
be an arbitrary partition. Then, there exists a number
such that for any
we have that
and hence,
Thus,
. □
Theorem 6. Let the conditions of Lemma 6 hold.
Then, for arbitrary initial function the function , defined via (28), is the unique solution of IP (25), (2) with interval of existence . Proof. The proof is based on ideas used in [
22]. Thereupon, we will emphasize in detail those differences that arise from the influence on the neutral system of the Riemann–Liouville-type derivatives.
Theorem 5 implies that
for any
. Then, for the function
, defined via (
28), by virtue of Lemma 1 in [
22] we have that
too.
More concretely, for arbitrary fixed
the kernels
of
define two Lebesgue–Stieltjes measures
and
for any
, as well as the function
defining a measure
for any
. We introduce the product measures
and
of the rectangles in
. Then, to use the proposition 5.15 in [
28] we need to prove that the relations
and
hold. Since according to Lemmas 5 and 6 we have
and
for any
and
concerning the first argument, then taking into account that
and
we obtain
Then, applying proposition 5.15 in [
28] we have for
that
and in the same way we obtain
From (
25) and (
31) via the Fubini theorem it follows that
Hence, from (
25), (
32), and (
33) it follows that
, defined via (
28), satisfies (
25) for any
.
Let
be an arbitrary number, and then, from (
28) when
, it follows that
i.e.,
satisfies the initial condition (
2), and hence, it is the unique solution of the IP (
25), (
2) with interval of existence
. □
Following the idea in [
22], we introduce
where
, and
.
Theorem 7. Let the conditions of Lemma 6 hold.
Then, the function introduced via (34) is the unique solution of IP (24), (2) for the initial function and interval of existence . Proof. As in the theorem above, the proof is based on some ideas used in [
15]. Therefore, we will only sketch the similar parts and we will emphasize in detail those differences that arise from the influence on the neutral system of the Riemann–Liouville-type derivatives.
For any
and
, since
when
, we have that
Then, substituting
on the left-hand side of (
24) for any
via the Fubini theorem and using formula (2.211) in [
5] and (
35) we obtain that
For the third addend on the right-hand side of (36), we have that
Since
and
, for any
and
we obtain
and hence, from (36), (
37), and (
38) it follows that
Similarly as above, for the right-hand side of (
24) we obtain
and since
is a fundamental matrix, then the statement of the theorem follows from (
39) and (
40). □
Corollary 2. Let the following conditions be fulfilled:
- 1.
The conditions of Lemma 6 hold.
- 2.
For any initial function and , the functions and are defined via (28) and (34), respectively.
Then, the function is the unique solution of IP (24), (2) and has the following representation for any :where for and . Proof. It follows from the superposition principle and Theorems 5 and 6. □
5. Weighted Stabilities
It is worth noting that the standard definitions for stability in the Lyapunov or Ulam–Hyers senses introduced for the systems with integer-order derivatives (without or with delays) are applicable directly for systems with fractional derivatives only when the fractional derivative of a constant is equal to zero (as Caputo-type derivatives, etc.). As we have mentioned previously [
11], in general when the fractional derivative of a constant is not equal to zero (as for RL-type derivatives), then the standard definitions are not applicable since the solutions of the systems with RL-type derivatives have singularity at
(the low terminal) of power order
. That is why new types of definitions for the different kinds of stabilities applicable to systems with RL-type derivatives are needed. The aim of this section is to introduce definitions of weighted stability in the Lyapunov sense, as well as UH and UHR weighted stabilities for fractional systems (equations) with RL-type derivatives and to study the relations with the classical definitions and with the concept “stability in time” in the Lyapunov sense, introduced in [
9].
Let assume that for any . Below, we recall the classical definitions for Lyapunov stability, applicable only for the case when the fractional derivatives are, in the Caputo sense, of order (and not applicable for RL derivatives):
Definition 9 ([
29])
. The zero solution of IP (1), (2) (with derivatives in the Caputo sense!) is said to be:- (a)
Stable if for a given low terminal and any there is a such that for each initial function with the corresponding solution satisfies for each the inequality .
- (b)
Stable (uniformly) if for a given low terminal and any there is a such that for every initial function , with the corresponding solution satisfying for each the inequality .
- (c)
Locally asymptotically stable (LAS) if for a given low terminal there is a such that for every initial function with the corresponding solution is stable and .
- (d)
Globally asymptotically stable (GAS) if for every for every initial function the corresponding solution is stable and .
Both definitions below are applicable not only for initial problems with Caputo-type derivatives, but also for initial problems with fractional derivatives in the RL sense.
Define for any
the weighted function
as follows:
Definition 10. The zero solution of IP (1), (2) is said to be weighted stable (W-stable), weighted uniformly stable (W-uniformly stable), weighted locally asymptotically stable (W-LAS stable), or weighted globally asymptotically stable (W-GAS stable) if for the corresponding solution the product satisfies the conditions (a), (b), (c), or (d) from Definition 9. Definition 11 ([
9])
. The zero solution of IP (1), (2) is said to be:- (i)
Stable in time (Lyapunov in time stable) if for an arbitrary there exists a point and number such that for any initial function with , the corresponding solution satisfies for each the inequality .
- (ii)
Asymptotically stable in time if for any initial function the corresponding solution is stable in time and additionally .
Remark 4. It is simple to see that if (the order of differentiation) then Definition 10 for the weighted stability coincides with the classical one. In addition, when a solution is weighted stable then we can take for any in Definition 11, and hence, the weighted stability implies stability in time, as well as the weighted global asymptotically stability implying asymptotic stability in time.
Let
be an arbitrary number and consider the inequality
Definition 12. The function is a solution of IP (42), (2) in if satisfies the inequality (42) for and the condition (2) for , with initial function . Definition 13. The system (1) is said to be Ulam–Hyers weighted stable in if for any there exists a number and constant such that for any solution of (42) there exists an initial function , with for and for the corresponding unique solution of IP (1), (2), with initial function , the inequalityholds for any . The first result in this section clarifies the relations between the boundedness of all solutions and the weighted stabilities in the Lyapunov sense.
Theorem 8. Let the conditions of Lemma 6 hold.
Then, the zero solution of the system (25) is W-stable in the Lyapunov sense if and only if for all solutions of IP (25), (2) with initial function , the product is bounded in . Proof. Necessity: Let us assume that the zero solution is W-stable in the Lyapunov sense. Assume the contrary, that there exists
such that for the existing, according to Theorem 3, unique solution
of IP (
25), (
2) we have that
is unbounded in
. Since the zero solution is W-stable, then for any
there exists a number
such that for any
with
we have that the corresponding solution
satisfies the inequality
for any
. Consider the following initial function
with norm
. Then, since the system (
25) is a linear system the function
is the corresponding solution of IP (
25) (
2) for the initial function
, and hence, for
we have that
, which implies that
for
, which contradicts our assumption.
Sufficiency: Let for any
the product
be bounded in
, where
is the corresponding solution of IP (
25), (
2) with this initial function, and hence, according to Lemma 4, we have that
satisfies IP (
3), (
2) with
for
and vice versa. Let
be arbitrary and let
be fixed numbers with
. Then, for arbitrary initial functions
, their corresponding unique solutions of IP (
3), (
2) are
,
, respectively, and for any
we have
From (
44), for any
it follows that
Then, applying corollary 2 in [
24] to (
45), we obtain for any
the estimation
Let us define for any
a family of functionals
via the equality
, where
denotes the corresponding solution of IP (
25),(
2) for initial function
. From (
46), it follows that
is a continuous functional at any
for any
i.e., for every
we have that
, where
denotes the norm of the functional. Then applying the Banach–Steinhaus theorem we obtain that the norms of all functionals are uniformly bounded, i.e., there exists a constant
such that
for any
. Let
be arbitrary and choose
. Then, for any
with
for the corresponding solution
we have that
for any
, which implies that the zero solution is W-stable in the Lyapunov sense. □
Remark 5. It is worth emphasizing that the result from the application of any kind of stability definitions concerning the differential equation (system) with all types of derivatives (integer or fractional order) essentially depends on the functional type of the set of all solutions of the studied object. As an example, for equations (systems) with first-order derivatives the solutions are usually either continuous differentiable or absolutely continuous functions. So, the conclusions could be true for all solutions which are absolutely continuous or only for these solutions which are continuous differentiable. For retarded or neutral equations (systems) the situation is more complicated. The space of the initial functions plays a leading role in its functional type, because the analytical type of the solutions essentially depends on them. This is very important for the stability in the Ulam–Hyers sense for the following reason: the functional type of the initial function and the solutions of system (1) must be the same for the corresponding inequality (42), since all solutions of (1) are solutions of inequality (42) too. As a consequence, from the proved results, we obtain two necessary conditions for W-stability of the zero solution of the system (
25).
Theorem 9. Let the following conditions be fulfilled:
- 1.
The conditions of Lemma 6 hold.
- 2.
The zero solution of the system (25) is W-stable in the Lyapunov sense.
Then, the following relations hold:where Proof. Let us, as above for
and any
, define a family of functionals
via the equality
, where
denotes the corresponding solution of IP (
25), (
2) for this initial function. Then, according Theorem 6 and Lemma 6, for any
we have that
and hence,
is a continuous functional at any
for any
with norm
. Then, by virtue of the Banach-Steinhaus theorem we have that there exists a constant
such that
for any
. Thus, we have that the W-boundedness (
) of the matrix
is a necessary condition for the W-stability of the zero solution of system (
25).
By virtue of Theorem 8 we have that for any and , the relation holds, and taking into account that when , we can conclude that . □
Corollary 3. Let the following conditions be fulfilled:
- 1.
The conditions of Lemma 6 hold.
- 2.
The zero solution of the system (25) is W-GAS stable in the Lyapunov sense.
Then, the following relationshold, where . Proof. The statement follows immediately from Theorem 9, since according to condition 2 of Corollary 3, for any we have for and for any . □
Remark 6. Note that the differences in the proofs of Theorems 8 and 9 are caused from the differences in the spaces of the initial functions and , respectively.
The next theorem makes clear the relation between the W-stability in the Lyapunov sense and Ulam–Hyers W-stability.
Theorem 10. Let the following conditions be fulfilled:
- 1.
The conditions of Lemma 6 hold.
- 2.
The system (25) is Ulam–Hyers W-stable for any initial functions .
Then, the zero solution of the system (25) is stable in the Lyapunov sense if and only if for any and arbitrary IP (42), (2) has only W-bounded solutions. Proof. Necessity: Let the zero solution of the system (
25) be stable in the Lyapunov sense and assume that IP (
42), (
2) has a W-unbounded solution
for some
and
. Then, since (
25) is Ulam–Hyers W-stable for any
, there exists
and
with
for
such that for
inequality (
43) holds for
and
, where
is the corresponding solution of IP (
25), (
3) with initial function
. On the other hand, since according to Theorem
8 the function
is a W-bounded solution of IP (
25), (
3), and since
is a W-unbounded solution of (
42) according to our assumption, we conclude that inequality (
43) is impossible to be fulfilled for
and
, which contradicts our assumption. Thus, all solutions of the inequality (
43) are W-bounded.
Sufficiency: Let IP (
42), (
2) have only bounded solutions. Then, IP (
25), (
2) has, for any
, only W-bounded solutions. Then, by virtue of Theorem 8, the zero solution of system (
25) is W-stable in the Lyapunov sense. □
6. On the Preservation of the Stability Properties
The aim of this section is, based on the obtained results in the previous sections, to make clear the preservation of the W-stability properties of the system (
25) under nonlinear disturbances.
Definition 14. The function will be called weighted Mittag–Leffler bounded (WML-bounded) of order if there exist , a constant , and a function with such that for . By , we denote the subset of all WML-bounded functions of order ω in .
Definition 15. We say that the vector-valued functional is a damper of order ω for system (1) if for any for the function there exist a point and a constant such that the following estimation holds for . The next theorem establishes an a priori estimate of the solutions of IP (
1), (
2) with initial functions
.
Theorem 11. Let the following conditions be fulfilled.
- 1.
The conditions of Lemma 6 hold.
- 2.
The conditions hold.
- 3.
and .
Then, any solution of IP (3), (2) with initial function is WML-bounded of order α. Proof. Let
be an arbitrary initial function and
be the corresponding unique solution of IP (
3), (
2). Then, for
from (
3) we obtain that
Without loss of generality, we can assume that
when
. Then, from (47) we have
Since
is monotone decreasing for
, then
, and hence,
Applying to (
48) Theorem 2, we obtain that
which completes the proof. □
Remark 7. We note that Theorem 11 establishes that for Lipschitz-type nonlinear disturbances, all solutions of IP (3), (2) with initial functions are WML-bounded of order α. The proof of the next theorem is based on the following practical variant of the well-known Barbalat’s lemma.
Lemma 8 ([
30])
. Let be a bounded function and . Then, the relation holds. Theorem 12. Let the following conditions be fulfilled:
- 1.
The conditions of Theorem 11 hold.
- 2.
The zero solution of system (25) is W-GAS. - 3.
There exists such that the vector-valued functional is a damper of order ω.
Then, the zero solution of IP (1), (2) (or IP (3), (2)) with initial function is W-GAS. Proof. Let
be arbitrary and
be the corresponding unique solution of IP (
1), (
2) (or (
3), (
2)). By virtue of Theorem 11 we have that
is WML-bounded of order
. We introduce the function
,
, and using (
28) and (
34) define the functions
and
as follows:
Then, from Theorem 6 it follows that
is the unique solution of IP (
25), (
2) with initial function
, and by virtue of Theorem 7 the function
is the unique solution of IP (
24), (
2) with
and initial function
,
. Thus,
is a solution of IP (
2), (
1) (IP (
3), (
2)), and hence, we have that
From (
49) and Corollary 3 it follows that
Consideringthe function
obviously the integrand
is a continuous function at any
, and hence,
for
too. From (
53) it follows that
Since
and the zero solution of the system (
25) is W-GAS, then by virtue of Theorem 9 we have that
, and hence,
is bounded on
. Then, from (
53) and (
54) it follows that
and, hence,
. Then, by virtue of Lemma 8 we have that
, and thus, from (
51) and (
52) it follows that
, which completes the proof. □
7. An Example
Example 1. Let , the lower terminal , , , , , , and . Then, IP (25), (2) obtains the following form: According to the results obtained by us in Section 3, IP (55), (56) has a unique solution, which we will find an explicit representation of in this example. To obtain an explicit representation of the solution we will use the approach based on the Laplace transform (LT) introduced in [31] for the case of Caputo-type derivatives. Below, we assume that , and then, by virtue of theorem 3 in [17] the LT can be applied correctly to Equation (55). Denote the LT of by Applying the LT to both sides of (55) and taking into account (56) we haveand hence, For any with , we have the series expansionand then, from (57) it follows that Substituting into (58) the expression we obtain that the LT of has the following form: For simplicity we will assume that , , and then, we have that Consider the three-parameter Mittag–Leffler (Prabhakar) functionintroduced in [32], where . For , , with the LT of the function is and when (see theorem 1.31 in [33]) we have that the following relation holds: Denoting as usual by the greatest integer number less or equal to t for any , from (59), (60), and (61), it follows that the solution of IP (55), (56) possess the representation Remark 8. Note that the approach used by the presentation of the solution of IP (55), (56) can be extended to obtain a presentation of the solutions with step functions (with one jump of the first kind in the interval ) as initial functions, and hence, as a consequence, a presentation of the generalized fundamental matrix can be obtained. But this idea, as well as a numerical simulation via the methods developed in [34,35,36,37], is beyond the area of this article and can be a theme for future research. 8. Conclusions and Comments
In this paper, we have studied a general class of nonlinear disturbed neutral linear fractional systems with derivatives in the Riemann–Liouville sense and distributed delays. As motivation for this study we can refer to the meaningful physical interpretations of models with these Riemann–Liouville fractional derivatives presented in [
7]. Mainly, we have studied the most important, from the point of view of applications, and the technically more complicated case, when the lower terminal of the Riemann–Liouville derivatives coincides with the end point of the initial interval. First, it was proved that the initial problem for these systems with discontinuous initial functions possesses a unique solution under some natural assumptions. It is worth mentioning that the assumptions used to derive this result are similar to those used in the case of systems with first-order derivatives, and as far as we know this is the first result on this theme. Then, as a consequence of the obtained result, we have proved the existence and uniqueness of a fundamental matrix and a generalized fundamental matrix for the studied neutral linear homogeneous system. The existence of the fundamental matrices have allowed us to establish an integral representation for the solutions of the initial problem for the homogeneous system and also for the corresponding inhomogeneous system. Furthermore, for fractional systems with Riemann–Liouville derivatives we have introduced a new concept for weighted stabilities in the Lyapunov, Ulam–Hyers and Ulam–Hyers–Rassias senses. Note that the introduced concept coincides with the classical stability concept for the cases of integer-order or Caputo-type derivatives. To prove the applicability of the introduced concept it was proved that the zero solution of the homogeneous system is weighted stable if and only if all its solutions are weighted bounded; this result is well known for systems with first-order derivatives and is also established for the case of systems with Caputo-type derivatives in our former works. In addition, for the homogeneous system it was established that the weighted stability in the Lyapunov sense and weighted stability in the Ulam–Hyers sense are equivalent if and only if the inequality appearing in the Ulam–Hyers definition possesses only bounded solutions. Finally, we have derived natural sufficient conditions under which the property of weighted global asymptotic stability of the zero solution of the homogeneous system is preserved under appropriate nonlinear disturbances.
In our point of view the main contributions in the article can be highlighted as follows:
A new concept was introduced for weighted stabilities in the Lyapunov, Ulam–Hyers and Ulam–Hyers–Rassias senses; it coincides with the classical stability concept for the cases of integer-order or Caputo-type fractional derivatives.
Sufficient conditions have been obtained which guarantee that the weighted stability in the Lyapunov sense and the weighted stability in the Ulam–Hyers sense are equivalent.
Sufficient conditions have been given under which the property of weighted global asymptotic stability of the zero solution of the homogeneous system is preserved under appropriate nonlinear disturbances.
As a future perspective for research we think that it will be important from the point of view of applications to establish explicit-type sufficient conditions which guarantee weighted stability in the Lyapunov or Ulam–Hyers senses.
Author Contributions
Conceptualization, H.K., M.M., S.I.C. and A.Z.; methodology, H.K., M.M., S.I.C. and A.Z.; validation, H.K., M.M., S.I.C. and A.Z.; formal analysis, H.K., M.M., S.I.C. and A.Z.; investigation, H.K., M.M., S.I.C. and A.Z.; resources, H.K., M.M., S.I.C. and A.Z.; data curation, H.K., M.M., S.I.C. and A.Z.; writing—original draft preparation, H.K., M.M., S.I.C. and A.Z.; writing—review and editing; supervision and project administration, H.K., M.M., S.I.C. and A.Z. All authors have read and agreed to the published version of the manuscript.
Funding
The authors of this research have been partially supported as follows: Hristo Kiskinov by the Bulgarian National Science Fund, Grant KP-06-N52/9, Mariyan Milev by EU-NextGeneration, through the NRRPlan of Bulgaria, SUMMIT BG-RRP-2.004-0008-C01 and Andrey Zahariev by the Bulgarian National Science Fund under Grant KP-06-N52/4, 2021.
Data Availability Statement
No new data were created or analyzed in this study. Data sharing is not applicable to this article.
Conflicts of Interest
The authors declare no conflicts of interest.
References
- Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
- Saxena, R.K.; Kalla, S.L. On the solutions of certain fractional kinetic equations. Appl. Math. Comput. 2008, 199, 504–511. [Google Scholar] [CrossRef]
- Li, C.; Qian, D.; Chen, Y.Q. On Riemann-Liouville and Caputo derivatives. Discret. Dyn. Nat. Soc. 2011, 2011, 562494. [Google Scholar] [CrossRef]
- 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]
- Podlubny, I. Fractional Differential Equation; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
- Podlubny, I. Geometric and physical interpretation of fractional integration and fractional differentiation. Fract. Calc. Appl. Anal. 2002, 5, 367–386. [Google Scholar]
- Podlubny, I. Physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives. Rheol. Acta 2006, 46, 765–771. [Google Scholar]
- Kiskinov, H.; Zahariev, A. On fractional systems with Riemann-Liouville derivatives and distributed delays-Choice of initial conditions, existence and uniqueness of the solutions—Choice of initial conditions, existence and uniqueness of the solutions. Eur. Phys. J. Spec. Top. 2017, 9, 3473–3487. [Google Scholar] [CrossRef]
- Agarwal, R.; Hristova, S.; O’Regan, D. Stability concepts of Riemann-Liouville fractional-order delay nonlinear systems. Mathematics 2021, 9, 435. [Google Scholar] [CrossRef]
- Agarwal, R.; Hristova, S.; O’Regan, D. Explicit solutions of initial value problems for linear scalar Riemann-Liouville fractional differential equations with a constant delay. Mathematics 2020, 8, 32. [Google Scholar] [CrossRef]
- Kiskinov, H.; Madamlieva, E.; Zahariev, A. Hyers–Ulam and Hyers–Ulam–Rassias stability for linear fractional systems with Riemann–Liouville derivatives and distributed delays. Axioms 2023, 12, 637. [Google Scholar] [CrossRef]
- Li, C.P.; Zhang, F.R. A survey on the stability of fractional differential equations. Eur. Phys. J. Spec. Top. 2011, 193, 27–47. [Google Scholar] [CrossRef]
- Liu, S.; Wu, X.; Zhou, X.F.; Jiang, W. Asymptotical stability of Riemann-Liouville fractional nonlinear systems. Nonlinear Dyn. 2016, 86, 65–71. [Google Scholar] [CrossRef]
- Li, M.; Wang, J. Representation of solution of a Riemann–Liouville fractional differential equation with pure delay. Appl. Math. Lett. 2018, 85, 118–124. [Google Scholar] [CrossRef]
- Liang, C.; Wang, J.; O’Regan, D. Representation of a solution for a fractional linear system with pure delay. Appl. Math. Lett. 2018, 77, 72–78. [Google Scholar] [CrossRef]
- Agarwal, R.; Hristova, S.; O’Regan, D. Existence and Ulam type stability for nonlinear Riemann–Liouville fractional differential equations with constant delay. Electron. J. Qual. Theory Differ. Equ. 2020, 67, 1–18. [Google Scholar] [CrossRef]
- 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] [CrossRef]
- Chartbupapan, W.; Bagdasar, O.; Mukdasai, K. A novel delay-dependent asymptotic stability conditions for differential and Riemann-Liouville fractional differential neutral systems with constant delays and nonlinear perturbation. Mathematics 2020, 8, 82. [Google Scholar] [CrossRef]
- Myshkis, A. Linear Differential Equations with Retarded Argument; Nauka: Moscow, Russia, 1972. (In Russian) [Google Scholar]
- Hale, J.; Lunel, S. Introduction to Functional Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
- Kolmanovskii, V.; Myshkis, A. Introduction to the Theory and Applications of Functional Differential Equations; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1999. [Google Scholar]
- Kiskinov, H.; Madamlieva, E.; Veselinova, M.; Zahariev, A. Integral representation of the solutions for neutral linear fractional system with distributed delays. Fractal Fract. 2021, 5, 222. [Google Scholar] [CrossRef]
- Krasnosel’skii, M.A. Some problems of nonlinear analysis. Am. Math. Soc. Transl. Ser. 1958, 2, 10. [Google Scholar]
- Ye, H.; Gao, J.; Ding, Y. A generalized Gronwall inequality and its application to a fractional differential equations. J. Math. Anal. Appl. 2007, 328, 1075–1081. [Google Scholar] [CrossRef]
- Madamlieva, E.; Milev, M.; Stoyanova, T. On stability criteria induced by the resolvent kernel for a fractional neutral linear system with distributed delays. Mathematics 2023, 11, 626. [Google Scholar] [CrossRef]
- Kiskinov, H.; Milev, M.; Veselinova, M.; Zahariev, A. Continuous dependence on the initial functions and stability properties in Hyers–Ulam–Rassias sense for neutral fractional systems with distributed delays. Fractal Fract. 2023, 7, 742. [Google Scholar] [CrossRef]
- Natanson, I.P. Theory of Functions of a Real Variable, 5th ed.; Frederick Ungar Publishing Co.: New York, NY, USA, 1983. [Google Scholar]
- Reitano, R.R. Foundations of Quantitative Finance: 5. General Measure and Integration Theory; Brandeis International Business School: Waltham, MA, USA, 2018. [Google Scholar]
- Madamlieva, E.; Kiskinov, H.; Petkova, M.; Zahariev, A. On the preservation with respect to nonlinear perturbations of the stability property for nonautonomous linear neutral fractional systems with distributed delays. Mathematics 2022, 10, 2642. [Google Scholar] [CrossRef]
- Sun, Z. A gathering of Barbalat’s lemmas and their (unsung) cousins. arXiv 2023. [Google Scholar] [CrossRef]
- Tuan, H.T.; Thai, H.D.; Garrappa, R. An analysis of solutions to fractional neutral differential equations with delay. Commun. Nonlinear Sci. Numer. Simul. 2021, 100, 105854. [Google Scholar] [CrossRef]
- Prabhakar, T.R. A singular integral equation with a generalized Mittag Leffler function in the kernel. Yokohama Math. J. 1971, 19, 7–15. [Google Scholar]
- Schiff, J.L. The Laplace Transform: Theory and Applications; Springer: New York, NY, USA, 1999. [Google Scholar]
- Garrappa, R.; Popolizio, M. Computing the matrix Mittag-Leffler function with applications to fractional calculus. J. Sci. Comput. 2018, 77, 129–153. [Google Scholar] [CrossRef]
- Micula, S. An iterative numerical method for fractional integral equations of the second kind. J. Comput. Appl. Math. 2018, 339, 124–133. [Google Scholar] [CrossRef]
- Jafari, H.; Tuan, N.A.; Ganji, R.M. A new numerical scheme for solving pantograph type nonlinear fractional integro-differential equations. J. King Saud Univ.-Sci. 2021, 33, 101185. [Google Scholar] [CrossRef]
- Moumen, A.; Mennouni, A.; Bouye, M. A novel Vieta–Fibonacci projection method for solving a system of fractional integrodifferential equations. Mathematics 2023, 11, 3985. [Google Scholar] [CrossRef]
| Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).