Explicit solutions of initial value problems for systems of linear Riemann–Liouville fractional differential equations with constant delay

A system of linear Riemann–Liouville fractional differential equations with constant delay is studied. The initial condition is set up similar to the case of the ordinary derivative. Explicit formulas for the solutions are obtained for various initial functions.

studied the linear homogenous Caputo fractional delay differential equations and gave a representation of the solution. Also, in [15,17] representations of the solution of linear non-homogeneous Caputo fractional delay differential equations are provided.
However, little is known regarding Riemann-Liouville (RL) fractional differential equations with delays. For some related contributions about RL factional differential equations, one can refer to previous works [1,3]. Note that linear systems of RL fractional differential equations without any delay are studied in [18] and explicit formulas for the solutions are obtained. RL fractional differential equations with delays are not well studied. We mention the papers [14,16] where the lower bound of the RL fractional derivative coincides with the left side end of the initial interval, but we note that this does not correspond to the idea in the case of delay differential equations with ordinary derivatives and the idea of the initial value problem of RL fractional differential equations.
In this paper we study initial value problems of systems of linear RL fractional differential equations with constant delay of the type where A, B are constant matrices, τ > 0 is a constant delay, and T ≤ ∞. Similar to the case of the ordinary derivative, the differential equation is given to the right of the initial time interval. It requires the lower bound of the RL fractional derivative to coincide with the right side end of the initial interval (usually this point is zero). Note that in this case any solution of an initial problem (IVP) with RL fractional derivatives is not continuous at the initial point. That is why RL fractional delay differential equations are convenient for modeling a process with impulsive types of initial conditions. These types of processes can be found in physics, chemistry, engineering, biology, and economics. To determine the law of the initial impulsive reaction, we need to add to the usual initial condition (for example, x(t) = φ(t) on the initial interval [-τ , 0], τ > 0 is the delay) a fractional condition. This conclusion is based on the results obtained in [9] and [21] concerning the physical interpretation of RL fractional derivatives and initial conditions which include derivatives of the same kind. Based on the above, we set up appropriate IVPs for RL linear fractional differential equations with the lower limit of the RL derivative equal to the right side point of the initial interval, i.e., we study initial conditions of the type where C is a constant vector.
Explicit formulas for the solutions of initial value problems with both zero and nonzero initial functions are obtained. Also, the cases of homogeneous as well as non-homogeneous equations are studied. In the case A = 0 the explicit formulas are comparatively easy to be applied, and in the case A = 0 the q-matrix functions (or the matrix Mittag-Leffler functions) are used. The scalar case of the linear RL fractional differential equations with the above mentioned initial conditions is studied in [2] and explicit solutions are obtained. Note that in [2] the Mittag-Leffler function is applied in all formulas for the explicit solutions. In the case of systems and vector functions, the Mittag-Leffler function is not applicable and it is replaced with the q-matrix exponential function defined and used in [4]. The application of this function not only leads to more complicated calculations but also to new formulas for the exacts solutions.

Preliminary notes on fractional derivatives and equations
Let m ∈ L loc 1 ([t 0 , T], R) and t 0 , T ≥ 0 : t 0 < T ≤ ∞ (in the case T = ∞ the intervals (t 0 , T] and [t 0 , T] are (t 0 , T) and [t 0 , T), respectively). In this paper we use the following definitions for fractional derivatives and integrals: -Riemann-Liouville fractional integral of order q ∈ (0, 1) [7,20] where Γ (·) is the gamma function. Note that sometimes the We give fractional integrals and RL fractional derivatives of some elementary functions which will be used later.

Proposition 1
For t > t 0 and β > 0, the following equalities are true: Let M n×n (R) be the set of all matrices A = {a ij } n i,j=1 with a ij ∈ R. We will use the notation I for the unit matrix from M n×n (R). For any matrix A ∈ M n×n (R), we will use the notation A 0 = I and define the q-matrix exponential function [4] by , t ≥ t 0 , and the Mittag-Leffler matrix function by Therefore, The definitions of the initial condition for systems of fractional differential equations with RL-derivatives are based on the following result for the linear RL matrix fractional equation: with X ∈ L loc 1 ([t 0 , T], R n ) and A ∈ M n×n (R).
Proposition 2 (Theorem 3.1 [25]) Let q ∈ (0, 1). Then problem (1) with the initial condition is equivalent to problem (1) with the initial condition Remark 1 According to Proposition 2, it is enough to study one of the initial conditions (2) or (3). Following this result, we will study only the initial condition of type (2).
The explicit formula of the initial value problem (IVP) (1), (2) is given by (see Theorem 3.2 [25]) In the case of a system of non-homogeneous fractional differential equations with RLderivatives with F ∈ C([t 0 , T], R n ), the explicit formula of IVP (5), (2) is given by (see Theorem 3 and Remark 1 [18]) In the case of a scalar linear RL fractional differential equation, we have the following result.
In the case of a system of RL fractional differential equations with constant coefficients, we have the following result, which is a special case of Theorem 3 [18].

Proposition 4 ([18]) The solution to the initial value problem for the system
has the following form: or its equivalent form

Explicit formula for the solutions of scalar linear RL fractional equations with delays and zero initial values
Throughout the paper we assume l i=n ( * ) = 0 for the integers n, l, n > l.
Remark 2 Without loss of generality we assume that there exists a natural number N ≤ ∞ such that T = (N + 1)τ .
Continue this process and the claim is established.
Remark 4 In the scalar case n = 1 the system of RL delay fractional equations (10) is reduced to a scalar linear delay RL fractional equation, the vector 1 is reduced to the constant 1, and formula (13) coincides with equation (7) [2].
Proof Let t ∈ (0, τ ]. Use the variation of constants method and we will search for solutions in the form where K ∈ C([0, τ ], R n ), K = (K 1 , K 2 , . . . , K n ) T is the unknown function to be obtained. According to the initial condition (11), we have Then, applying From (27) and (28) we get K i (t) ≡ F i (t), i = 1, 2, . . . , n, i.e., the solution X(t) of IVP (23), (11), (24) for t ∈ (0, τ ] is given by Note that it is easy to check the validity of condition (24) for the solution X(t) defined by (29).
Remark 5 Note that the formula for the solution in the homogeneous case does not follow from the one in the non-homogeneous case because of fractional conditions (12), respectively, (24).
Remark 6 The explicit formula (25) of solution in Theorem 2 is a generalization of formula (20) [2] for the scalar case n = 1.

Explicit formula for the solutions of scalar linear RL fractional equations with delays and non-zero initial values
Consider the linear non-homogeneous RL fractional differential equation (23) with nonzero initial value: where Remark 7 Note that the function G(t) = t q-1 1 is not applicable in this case as an initial function.

Zero initial function
Consider the system of non-homogeneous linear Riemann-Liouville fractional differential equations with constant delay: with the initial conditions where X ∈ R n , F ∈ C(R + , R n ), A, B ∈ M n×n (R), τ > 0 is a real constant, C ∈ R n .
Remark 8 In the case τ = 0 (no delay), B = 0, and F ≡ 0, the system of RL fractional differential equations (47) is reduced to the one studied in [25] with initial condition (49), and the formula obtained in Theorem 4 reduces to the formula obtained in Theorem 3.2 [25].
Remark 9 In the case τ = 0 (no delay), B = 0, the system of RL fractional differential equations (47) is reduced to the studied system (34) in [4] with initial condition (49), and the formula obtained in Theorem 4 reduces to the formula obtained in Theorem 3 [4].
Remark 10 Note that in the case the initial time 0 is replaced with arbitrary t 0 , all the results in the paper are true with slight changes.

Conclusions
The formulas for the exact solutions are important tools in fractional models. Often it is quite complicated to find the exact solution for RL fractional differential equations even in the linear scalar case. In this paper we study various types of systems of linear RL fractional differential equations with constant delays. We set up an initial value problem in an appropriate way based on the physical meaning to initial conditions expressed in terms of Riemann-Liouville fractional derivatives or integrals [9]. Explicit formulas for the solutions of initial value problems with both zero and nonzero initial functions are obtained, and systems with homogeneous and non-homogeneous equations are studied. The q-matrix exponential function is successfully applied in explicit solutions [4].
The obtained formulas will be very helpful in the theoretical study of fractional models with RL fractional derivative, for linearization of multi-dimensional nonlinear models, for the monotone-iterative technique, and for systems of RL fractional differential equations with delays.