Asymptotic Analysis of Linear and Interval Linear Fractional-Order Neutral Delay Differential Systems Described by the Caputo-Fabrizio Derivative

Asymptotic stability of linear and interval linear fractional-order neutral delay differential systems described by the Caputo-Fabrizio (CF) fractional derivatives is investigated. Using Laplace transform, a novel characteristic equation is derived. Stability criteria are established based on an algebraic approach and norm-based criteria are also presented. It is shown that asymptotic stability is ensured for linear fractional-order neutral delay differential systems provided that the underlying stability criterion holds for any delay parameter. In addition, sufficient conditions are derived to ensure the asymptotic stability of interval linear fractional order neutral delay differential systems. Examples are provided to illustrate the effectiveness and applicability of the theoretical results.


Introduction
Fractional calculus is attracting more and more researchers in applied sciences and engineering because of many advantages of fractional derivatives which provide important tools in modelling natural phenomena, see, e.g., [1] [2]. There are different types of fractional derivatives, those of Riemann-Liouville and Caputo are the most popular in the literature [3] [4] [5] [6].
In their recent work, Caputo and Fabrizio [7] introduced a new fractional order derivative with a nonsingular kernel, hereinafter called the fractional Caputo-Fabrizio (CF) derivative. This new fractional derivative is less affected by the past compared to the Caputo fractional derivative, which may exhibit slow stabilization [8]. The properties and numerical aspects of the CF derivative and their corresponding fractional integrals have been studied in [9]- [15]. In this paper, we are interested in linear and interval linear fractional-order neutral delay differential equations described by the CF derivative. Recently, a great attention has been paid to fractional delay differential systems. One of the most important research topics of the theory of such systems is the stability analysis.
However, in the literature cited above, there are only few results (see, e.g., [8] [16]) on this topic. In [17] the authors study the stability analysis of linear fractional-order ordinary differential equations described by the CF derivative, whereas the authors of [18] consider the stability analysis of linear fractional-order systems with time delay, establish a characteristic equation using the Laplace transform and provide some brief sufficient stability conditions. In [19] the authors consider the asymptotic stability for uncertain singular neutral delay systems and in [20] [21] [22] the authors study the stability analysis of interval linear fractional ordinary differential systems, interval linear fractional neutral differential systems described by the Caputo derivative and interval linear fractional neutral differential-algebraic systems described by the Caputo derivative, respectively. While being different, we extend the analysis carried out in the above cited references. We apply a spectrum based approach to establish asymptotic stability criteria for fractional-order neutral delay differential systems and the novelty of this work lies in the following aspects. Firstly, the CF definition of the fractional derivative is applied to analyze linear and interval linear fractional-order differential systems including neutral time delay. Secondly, by using the Laplace transform, we establish a novel characteristic equation. Thirdly, we apply an algebraic approach to establish sufficient asymptotic stability criteria ensuring that all the roots of characteristic equation lie in open left half of the complex plane.
Fourthly, we present sufficient norm criteria to ensure the asymptotic stability of such systems. A nice consequence is that these stability criteria avoid solving the characteristic (transcendental) equation.
The rest of the paper is organized as follows: In Section 2 we formulate the problem and introduce the notation that will be used throughout the paper. In Section 3 we establish the main results for the asymptotic stability criteria for linear and interval linear fractional neutral delay differential equations. Finally, some concluding remarks are given in Section 4.

Problem Formulation and Notation
We consider linear fractional-order neutral delay differential systems of the form where 0 1 α < < , ( ) y t is a real vector of size n (the state vector); A, B and C are real n n × matrices; τ is a positive parameter (the time delay) and ϕ is a consistent initial function. The notation ( ) CF D z t α stands for Caputo-Fabrizio fractional order derivative of order α of z(t) defined by (see [7]) We study two cases. In the first case the matrices A, B, and C are constant and in the second case they are interval matrices, that is are given. The equations associated with the constant and interval cases will be denoted respectively by (LFNDDE) and (ILFNDDE). When a property applies equally to both cases (as in Definition 1 for example), the corresponding system will simply be referred to as system (1).
The following notation will be used throughout the paper. For a square complex matrix X, * X denotes conjugate transpose of X, denotes the spectral radius of X, X denotes the spectral norm of X; [ ] X µ denotes the logarithmic norm of X; The following inequalities, which will be used in the sequel, can be found in the literature, see, e.g., [23]. If ,

Main Results
Applying the Laplace transform to (1), we obtain The characteristic equation of (1) is defined by As we will see, it plays an important role in the stability analysis of system (1). From Definition 1 we obtain the following theorem whose proof is similar to that given in [24]. Theorem 1. System (1) is asymptotically stable if the roots of the characteristic equation (4) lie in the open left half of the complex plane and are uniformly bounded away from the imaginary axis.
Next, several sufficient conditions of stability of system (1) are given.
Proof. To simplify the notation, let us denote by The condition (C 2 ) can then be written Remark 1. Note that the condition (C 3 ) is a transcendental inequality, which is difficult to solve. The sufficient conditions in Lemma 1 avoid this difficulty.
Since the spectral radius of a matrix is always less than or equal to any induced norm of the matrix, Lemma 1 remains valid under the conditions (C 1 ) and As in the proof of Lemma 1, we use the notation is a continuous function of s for Re s ≥ , similar to [25] and from the condition (C 2 ), there exists 0 For sufficiently large m there exists a positive constant * ε ( Then it follows from (5) and (6) then the following equality holds for .
Proof. Note first that there is a z satisfying (7). Indeed such a z is given by Combining the expressions of I zJ − and F zG + above leads to the desired result.
In view of Theorem 2 and Lemma 2, we obtain the following  In the following, we give other sufficient conditions ensuring asymptotic stability. For this purpose, we need the following lemma whose proof can be found, for example, in [26].
In particular, (10) and (11) show that ( ) T s is nonsingular and therefore the roots of the characteristic Equation (4) contradicting the assumption that ( m s ) is a sequence of roots of the characteristic equation. Now we consider the stability of system (1) in the case (ILFNDDE). From (2) we have Then it is easy to see that , , .
The following theorem gives sufficient conditions for the stability of (ILFNDDE).
Theorem 5. If the family of matrices { } , , A B C defined in (ILFNDDE) satisfies the following conditions Re s ≥ , It is straightforward to verify that The condition (P 2 ) says that ( ) 1 R s is nonsingular. From (12), (13) and (14) we  (4) belong to the open left half of the complex plane. The proof that the roots remain bounded away from the imaginary axis is a simple repetition of the proof of Theorem 2.    We deduce form Theorem 5 that the system is asymptotically stable.

Conclusion
This paper is concerned with the asymptotic stability of linear fractional-order neutral differential delay systems described by the Caputo-Fabrizio derivative. The systems matrices can be constant or interval matrices. Using the Laplace transform, we have derived a new characteristic equation for such systems. This characteristic equation, though interesting in theory, involves a transcendental term which makes it difficult to use in practice and in particular to study the asymptotic stability of such systems. To overcome this difficulty, sufficient algebraic criteria have been given to ensure the asymptotic stability of such systems in the case of constant matrices (Lemma 1 and Lemma 2). We have successfully shown that under these algebraic criteria, asymptotic stability holds (Theorem 2 and Theorem 3). Also, norm-based criteria have been given to ensure the asymptotic stability. In particular, in case of constant matrices, Theorem 4 involves easily verifiable criteria based only on matrix norm computations. Finally, sufficient conditions are derived to ensure the asymptotic stability of such systems in the case of interval matrices. We have shown that this system is asymptotically stable for any delay parameter (Theorem 5). The effectiveness of the theoretical results has been illustrated by numerical examples.