Asymptotic behavior of fractional-order nonlinear systems with two different derivatives

This paper addresses the asymptotic behavior of systems described by nonlinear differential equations with two fractional derivatives. Using the Mittag–Leffler function, the Laplace transform, and the generalized Gronwall inequality, a sufficient asymptotic stability condition is derived for such systems. Numerical examples illustrate the theoretical results.


Introduction
Fractional calculus is a generalization of the standard integer calculus, extending the derivative and integral operators to non-integer orders [1]. In 1 calculus experienced a rapid development, and its applications became common in mechanics, physics, chemistry, and other disciplines. Fractional-order (FO) models can describe many real-world systems better than their integer-order counterparts [2], namely when phenomena like long-range correlations, memory, and heredity effects are present.
State-space models of FO systems can be divided into two categories, namely commensurate and incommensurate. In general, if the state variables are of the same order, the system is called commensurate, otherwise it is denoted incommensurate. Multi-order systems are usually incommensurate, while commensurate systems are just special cases of the incommensurate ones. Studying the stability problem of incommensurate FO systems is crucial, but also demanding.
Stability analysis is one of the most important issues when studying dynamical systems [3][4][5]. The stability analysis of FO systems, meaning those described by FO differential equations, is often more difficult than that of integer-order systems, since fractional derivatives are non-local and have weakly singular kernels [6,7]. In the past decade, many results regarding the stability of commensurate FO systems were derived, including for linear [8], nonlinear [9], and delayed systems [10]. However, for incommensurate FO systems, fewer studies were reported in the literature, since their stability analysis is very demanding. The stability analysis of two-term FO differential equations was investigated in various papers [11][12][13][14][15]. The necessary and sufficient stability and instability conditions for linear FO differential equations with three Caputo derivatives were considered in the work [16]. Robust stability of FO systems described in pseudo-state space with incommensurate FO was addressed in the references [17,18].
We should note that all references mentioned previously, dealing with the stability of FO incommensurate systems, focus mainly on linear system [11][12][13][14][15][16][17][18]. Indeed, to the best of the authors knowledge, the stability of multidimensional FO nonlinear systems was not explored so far. In this paper, by using the Mittag-Leffler function, the Laplace transform, and the generalized Gronwall inequality, an asymptotic condition for FO nonlinear systems with two different derivatives is proposed.
The structure of this paper is as follows. Section 2 gives some preliminaries on fractional calculus and basic definitions. Section 3 is dedicated to proving the main results of the paper. Section 4 presents numerical examples for illustrating the theoretical findings. Finally, Sect. 5 lists the conclusions of the study.

Preliminaries and basic tools
Definition 2 [19] The Caputo fractional integral of order α ∈ R + of function x(t) is Definition 3 [19] The Caputo fractional derivative of order α ∈ R + of function x(t) is where n = α , with α being the smallest integer greater than or equal to α.
Proposition 1 [19] The Laplace transform of the Caputo fractional derivative is where X (s) represents the Laplace transform of function x(s).

Main results
The FO nonlinear system with two different derivatives is considered: where x(t) = (x 1 (t), x 2 (t), . . . , x n (t)) T ∈ R n denotes the state vector, A, B ∈ R n×n are constant matrices, f (x(t)) ∈ R n represents a nonlinear vector, and satisfies f (0) = 0, and the Lipschitz condition with Lipschitz constant L, that is, , and the FO p, q belong to 0 < p < 1 < 1+ p ≤ q < 2.
The system (1) can be used to describe a standard heat diffusion process, the standard voltage equations for impedance at the drive of lossy transmission lines, as well as other phenomena useful in practical engineering. where F(x(s)) = L( f (x(t))). Thus, we can obtain X (s) = As p−1 x(0) + Bs q−1 x(0) + Bs q−2 x (0) + F(x(s)) As p + Bs q = As p−1 x(0) As p + Bs q + Bs q−1 x(0) As p + Bs q + Bs q−2 x (0) As p + Bs q + F(x(s)) As p + Bs q . ( Applying the inverse Laplace transform on (2), it yields In view of (3), the solution of (2) can be obtained as follows: It follows from Lemma 1 that Combining (4) and (5) and denoting C = B −1 A, one obtain Since C is a stability matrix, one has e Ct q− p ≤ e −wt q− p .
Substituting (6) into (5), it yields Multiplying both sides of (7) by e wt q− p , one has Let us denote u(t) = e wt q− p x(t) , g(t) = L B −1 , and a(t) = t q− p x(0) + x(0) + t x (0) . Then, (8) can be rewritten as follows: From Lemma 2, it yields It follows from Lemma 3 that there exist two real constants C 1 , C 2 > 0 such that Multiplying both sides of (11) by e −wt q− p , one has According to the condition in Theorem 1, when t → ∞, x(t) → 0, which implies that system (1) has asymptotically stable behavior.
Remark 2 Theorem 1 is also valid for FO systems with irrational order.

Conclusion
We investigated the asymptotic behavior of systems described by nonlinear differential equations with two fractional derivatives, p, q, verifying 0 < p < 1 < 1+ p ≤ q < 2. The Mittag-Leffler function, the Laplace transform and properties of fractional calculus were adopted to derive a sufficient asymptotic stability condition. Numerical examples illustrated the effectiveness of the theoretical condition. Possible generalizations to the case of multi-term FO delayed differential equations will be addressed in future work.