Next Article in Journal
Uniformly Resolvable Decompositions of Kv-I into n-Cycles and n-Stars, for Even n
Previous Article in Journal
Unpacking Mathematics Pedagogical Content Knowledge for Elementary Number Theory: The Case of Arithmetic Word Problems
Previous Article in Special Issue
Existence and Stability Analysis for Fractional Impulsive Caputo Difference-Sum Equations with Periodic Boundary Condition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Stability of Linear Incommensurate Fractional-Order Difference Systems

1
Laboratory of Mathematics, Informatics and Systems (LAMIS), University of Laarbi Tebessi, Tebessa 12002, Algeria
2
Department of Mathematics, University of Jordan, Amman 11942, Jordan
3
Dipartimento Ingegneria Innovazione, Universita del Salento, 73100 Lecce, Italy
4
Nonlinear Systems and Applications, Faculty of Electrical and Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City 758307, Vietnam
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(10), 1754; https://doi.org/10.3390/math8101754
Submission received: 5 July 2020 / Revised: 26 July 2020 / Accepted: 29 July 2020 / Published: 13 October 2020
(This article belongs to the Special Issue Stability Analysis of Fractional Systems)

Abstract

:
To follow up on the progress made on exploring the stability investigation of linear commensurate Fractional-order Difference Systems (FoDSs), such topic of its extended version that appears with incommensurate orders is discussed and examined in this work. Some simple applicable conditions for judging the stability of these systems are reported as novel results. These results are formulated by converting the linear incommensurate FoDS into another equivalent system consists of fractional-order difference equations of Volterra convolution-type as well as by using some properties of the Z-transform method. All results of this work are verified numerically by illustrating some examples that deal with the stability of solutions of such systems.

1. Introduction

Undoubtedly, it has been demonstrated, over the past few decades, that the non-integer calculus is a forceful mathematical argument for providing many and more dynamics for lots of ancient as well as modern models. For instance, there are several dynamical models and many real-world applications that were recently handled via this tool such as diffusion modeling [1], robot manipulators [2], economics [3], and many more. As a matter of fact, such branch of calculus, in its two structures (discrete and continuous), requires the order of the functional operator of calculus to be in its fractional-order case for a number of core concepts such as derivatives, integrals, and differences. This indeed allows for creating appropriate mathematical models together with extremely rich and complex dynamics [4,5]. In the area of the discrete fractional calculus, the explanation of its main concept is referred to in Diaz and Olseronly in 1974 [6]. After that time, for at least a decade or more, some serious efforts in this field began to appear. Among those efforts, to name a few, is what Miller and Ross proposed in [7]. They actually established some basic definitions, primary schemes, and properties related to the fundamental theorem [8]. To date, a lot of researchers are competing in proposing impressive results and methods in regard to this field. Anyhow, for a complete comprehensive description about this branch of calculus, the reader may refer to [9,10,11].
Over the past few years, modeling several chemical and physical phenomena have broadly been carried out using the theory of Fractional-order Difference Systems (FoDSs) [12]. In fact, the FoDS is employed, with its particular digital data, for the purpose of approximating its corresponding fractional-order differential equations. This allows one to enter them into suitable computer programs and then simulate the obtained results [12]. However, it was reported in [12] that the Z-transform method can solve the linear FoDSs, as it can perform the same matter for the linear Fractional-order Difference Equation (FoDE). The Z-transform method can be, at the same time, employed as a powerful aid tool in discussing the stability analysis of such systems, as declared in [13,14,15,16]. However, even so, there are formidable challenges, generally, in proposing a proper tool for performing this task [13]. Therefore, it might be said that the stability analysis of the FoDSs is not yet developed [17]. This has motivated, in recent years, many researchers to deal with this dilemma, and as a consequence of this, several other works have preferred it to be the main target of their investigations (see [15,17,18,19]). In particular, Abu-Saris and Al-Mdallal developed in [15] a theory about the stability of linear FoDSs. However, unfortunately, it was shown that their theorem is difficult to apply. That matter remained until 2015. This year, Čermák et al. developed another simplified and applicable theorem for the same purpose (see [17]). As far as we know, discussing the stability analysis by providing novel simplified results for Linear FoDS with incommensurate orders remains up until now a recent and mostly unexamined topic. In light of this urgent need, this paper presents some simple applicable conditions for judging the stability of such system by first converting it into another equivalent form that includes FoDEs of the Volterra convolution-type as well as by using the properties of the Z-transform method. All results of this work are provided by examples, so that all plots of the stability of the solutions for systems are exhibited according to given incommensurate orders. However, this paper is organized in the following order. Section 2 introduces some primary preliminaries associated with discrete fractional calculus, while Section 3 discusses some recently established results that have handled the stability of linear commensurate FoDS for the purpose of getting new results in regard to incommensurate orders. Section 4 exhibits several examples to verify all findings, followed by the last section that summarizes the achievements of the whole work.

2. Preliminaries

This section briefly introduces some basic definitions and preliminaries associated with discrete fractional calculus. In all of the definitions below, the function f is defined on N a , where N a = { a , a + 1 , a + 2 , } , for a R .
Definition 1
([20]). The fractional sum operator Δ a δ of a function f : N a R is expressed as:
Δ a δ f ( t ) : = 1 Γ ( δ ) τ = a t δ ( t τ 1 ) ( δ 1 ) f ( τ ) ,
where δ > 0 is the order of the operator and t N a + δ .
Definition 2
([20]). The fractional Riemann–Liouville difference operator of a function f is defined by:
Δ a δ f ( t ) : = Δ n Δ a ( n δ ) f ( t ) = 1 Γ ( n δ ) Δ n τ = a t ( n δ ) ( t τ 1 ) ( n δ 1 ) f ( τ ) ,
where δ > 0 is the order of the operator, n = [ δ ] + 1 , and t N a + n δ .
Definition 3
([20]). The fractional Caputo difference operator of a function f is expressed as:
C Δ a δ f ( t ) : = Δ a ( n δ ) Δ n f ( t ) = 1 Γ ( n δ ) τ = a t ( n δ ) ( t τ 1 ) ( n δ 1 ) Δ n f ( τ ) ,
where δ > 0 is the order of the operator, n = [ δ ] + 1 and t N a + n δ .
Lemma 1
([17]). Let δ > 0 and δ N . Then,
( C Δ 0 δ x ) ( j + 1 δ ) = τ = 0 j + 1 ( 1 ) j τ + 1 δ j τ + 1 x ( τ ) ( 1 ) j + 1 δ 1 τ + 1 x ( 0 ) , j = 0 , 1 , ,
where δ β is the general binomial coefficient.
Proof. 
By Definition 2, one can obtain:
( C Δ 0 δ x ) ( j + 1 δ ) = 1 Γ ( 1 δ ) τ = 0 n ( j δ τ ) ( δ ) Δ x ( τ ) = 1 Γ ( 1 δ ) τ = 0 j Γ ( j δ τ + 1 ) Γ ( j τ + 1 ) Δ x ( τ ) = τ = 0 j j δ τ j τ Δ x ( τ ) = τ = 1 j j δ τ + 1 j τ + 1 j δ τ j τ x ( τ ) j δ j x ( 0 ) + δ 0 x ( j + 1 ) .
Using Pascal rule yields
( C Δ 0 δ x ) ( j + 1 δ ) = τ = 1 j j δ τ j τ + 1 x ( τ ) j δ j x ( 0 ) + x ( j + 1 ) = τ = 0 j j δ τ j τ + 1 x ( τ ) j δ j + j δ j + 1 x ( 0 ) + x ( j + 1 ) = τ = 0 j + 1 ( 1 ) j τ + 1 δ j τ + 1 x ( τ ) ( 1 ) j + 1 δ 1 j + 1 x ( 0 ) .
 □

3. Stability Analysis of the Linear Incommensurate FoDS

In this section, we state some additional results, reported in [15,17], in order to pave the way for introducing the main results of this work, which will be verified numerically, later on, in Section 4. Such results will clearly show the stability of the linear FoDS with incommensurate orders via some useful conditions formulated as theorems. First of all, consider the following linear incommensurate FoDS:
C Δ 0 δ 1 x 1 ( j + 1 δ 1 ) = a 11 x 1 ( j ) + a 12 x 2 ( j ) + + a 1 n x n ( j ) , C Δ 0 δ 2 x 2 ( j + 1 δ 2 ) = a 21 x 1 ( j ) + a 22 x 2 ( j ) + + a 2 n x n ( j ) , C Δ 0 δ n x n ( j + 1 δ n ) = a n 1 x 1 ( j ) + a n 2 x 2 ( j ) + + a n n x n ( j ) , j = 0 , 1 , ,
where C Δ 0 δ i is the fractional Caputo difference operator of order δ i , and 0 < δ i 1 , for i = 1 , 2 , , n .
The main objective of this work does not explore the separate investigation of the FoDSs, but it especially explores the generalization of some notions and results achieved in the framework of these discrete settings. One of the main features that is reflected from studying such systems is memory. This feature can be revealed by tracking the evolution of any state of these systems, which not only depends on the current solution, but it also depends on the whole history of its development. In general, many literature works seem to agree that these systems possess superior properties as compared to their standard counterpart. It has been shown that the core of their general dynamics solution is heavily dependent on the variations of the fractional-order values. However, most of those literature works are concerned only with studying the dynamic characteristics of linear commensurate FoDSs, which is of course a special case of systems with incommensurate fractional-orders. From the perspective of such orders having different effects on the linear incommensurate FoDSs, several physical phenomena are formulated and accurately described for the purpose of improving the complexity of their solutions over the original systems. For instance, if δ 1 = δ 2 = = δ n = δ , then system (5) might be expressed in the commensurate form as follows:
C Δ 0 δ X ( j + 1 δ ) = A X ( j ) , j = 0 , 1 , ,
where X ( t ) = ( x 1 ( k ) , x 2 ( k ) , , x n ( k ) ) T R n , A R n × n , and where C Δ 0 δ X ( j + 1 δ ) = ( C Δ 0 δ x 1 ( j + 1 δ ) , C Δ 0 δ x 2 ( j + 1 δ ) , , C Δ 0 δ x n ( j + 1 δ ) ) T . In [15], Abu-Saris and Al-Mudalal established some results associated with the stability of the system given in (6) by transforming it into a northern difference system of the Volterra convolution-type as well as using some properties of the Z-transform method. They actually deduced Theorem 1 given below.
Theorem 1
([15]). Suppose 0 < δ < 1 . The isolated zeroes, of the nonnegative real axis, of
det I n z 1 ( 1 z 1 ) δ A
lie inside the unit disk if and only if the zero solution of system (6) is asymptotically stable, where I n is the identity matrix.
One might notice that, despite their good work, applying such results is, indeed, an extremely difficult task. The justification of that belief is referred to difficulties in determining the zeros of (7). This matter, however, encouraged Čermák et al. to overcome its difficulties. For that purpose, they transformed system (6) into another difference system of Volterra convolution-type and then applied the deduced results of Elaydi et al. in [21]. In other words, they formulated the following result.
Theorem 2
([17]). Suppose that δ ( 0 , 1 ) and ξ is an eigenvalue of A. Then, the zero solution of system (6) is asymptotically stable if ξ S δ = z C : z < 2 cos arg z π 2 δ δ and arg z > δ π 2 , ξ . In this regard, every solution X of system (6) tends to zero algebraically, i.e.,
X ( j ) = O ( j δ ) as j .
In addition, if ξ C c l ( S δ ) , then the zero solution of system (6) is not stable, where c l ( S δ ) denotes to the closure of S δ .
Remark 1.
The affirmations of the two above results outline the same region of stability, although they have two analytical different descriptions. For instance, the condition declared in Theorem 2 appears to be more suitable for practical intents because of S δ , which has been established in an explicit form.
Thus, one can definitely observe that the two aforementioned theorems could be implemented in dealing with linear commensurate FoDSs. To follow the progress in this context, the linear incommensurate FoDSs will be handled next. Actually, discussing the stability of these systems, as far as we know, remains, up to now, a recent and generally unexamined topic. In this regard, we provide the following result:
Theorem 3.
Consider system (5) subject to the initial vector condition X ( 0 ) = X 0 R n . Then,
  • If all roots of the following characteristic equation:
    det d i a g z ( 1 z 1 ) δ 1 , z ( 1 z 1 ) δ 2 , , z ( 1 z 1 ) δ n A = 0 ,
    lie inside the unit disk, then the zero solution of system (5) is asymptotically stable.
  • If there exists a zero, say z * , of (9) such that z * > 1 , then the zero solution of system (5) is not stable.
Proof. 
Using Lemma 1, we can rewrite system (5) as follows:
x 1 ( j + 1 ) = τ = 0 j ( 1 ) j τ δ 1 j τ + 1 x 1 ( τ ) + ( 1 ) j + 1 δ 1 1 j + 1 x 1 ( 0 ) + a 11 x 1 ( j ) + a 12 x 2 ( j ) + + a 1 n x n ( j ) , x 2 ( j + 1 ) = τ = 0 j ( 1 ) j τ δ 2 j τ + 1 x 2 ( τ ) + ( 1 ) j + 1 δ 2 1 j + 1 x 2 ( 0 ) + a 21 x 1 ( j ) + a 22 x 2 ( j ) + + a 2 n x n ( j ) , x n ( j + 1 ) = τ = 0 j ( 1 ) j τ δ n j τ + 1 x n ( τ ) + ( 1 ) j + 1 δ n 1 j + 1 x n ( 0 ) + a n 1 x 1 ( j ) + a n 2 x 2 ( j ) + + a n n x n ( j ) , j = 0 , 1 , .
One might take the Z-transform to (10). This yields the following system:
z x ˜ 1 ( z ) z x 1 ( 0 ) = ( z z ( 1 1 z ) δ 1 ) x ˜ 1 ( z ) + ( z ( 1 1 z ) δ 1 1 z ) x 1 ( 0 ) + a 11 x ˜ 1 ( z ) + a 12 x ˜ 2 ( z ) + + a 1 n x ˜ n ( z ) , z x ˜ 2 ( z ) z x 2 ( 0 ) = ( z z ( 1 1 z ) δ 2 ) x ˜ 2 ( z ) + ( z ( 1 1 z ) δ 2 1 z ) x 2 ( 0 ) + a 21 x ˜ 1 ( z ) + a 22 x ˜ 2 ( z ) + + a 2 n x ˜ n ( z ) , z x ˜ n ( z ) z x n ( 0 ) = ( z z ( 1 1 z ) δ n ) x ˜ n ( z ) + ( z ( 1 1 z ) δ n 1 z ) x n ( 0 ) + a n 1 x ˜ 1 ( z ) + a n 2 x ˜ 2 ( z ) + + a n n x ˜ n ( z ) ,
where x ˜ i ( z ) indicates the Z -transform of x i ( j ) , (i.e., x ˜ i ( z ) = Z [ x i ( j ) ] , 1 i n ). Consequently, we can rewrite system (11) as follows:
M ( z ) . x ˜ 1 ( z ) x ˜ 2 ( z ) x ˜ n ( z ) = z ( 1 1 z ) δ 1 1 x 1 ( 0 ) z ( 1 1 z ) δ 2 1 x 2 ( 0 ) z ( 1 1 z ) δ n 1 x n ( 0 ) ,
in which
M ( z ) = z ( 1 1 z ) δ 1 a 11 a 12 a 1 n a 21 z ( 1 1 z ) δ 2 a 22 a 2 n a n 1 a n 2 z ( 1 1 z ) δ n a n n .
Multiplying both sides of (12) by ( z 1 ) gives:
M ( z ) . ( z 1 ) x ˜ 1 ( z ) ( z 1 ) x ˜ 2 ( z ) ( z 1 ) x ˜ n ( z ) = z 2 ( 1 1 z ) δ 1 x 1 ( 0 ) z 2 ( 1 1 z ) δ 2 x 2 ( 0 ) z 2 ( 1 1 z ) δ n x n ( 0 ) .
Now, we should note that, if all roots of det M ( z ) = 0 lie inside the unit disk, then system (14) will be considered such that z satisfies z R , with R 1 , (R is the radius of convergence of X ˜ ( z ) ) . Actually, system (14) has a unique solution in this limited area represented by ( ( z 1 ) x ˜ 1 ( z ) , ( z 1 ) x ˜ 2 ( z ) , , ( z 1 ) x ˜ n ( z ) ) . Accordingly, we have:
lim z 1 ( z 1 ) x ˜ i ( z ) = 0 , i = 1 , 2 , , n .
Based on the assumption stated in first part of this theorem, and based also on the Final-Value Theorem associated with Z-transform, we obtain:
lim j x i ( j ) = lim z 1 ( z 1 ) x ˜ i ( z ) = 0 , i = 1 , 2 , , n .
On the other hand, considering the second part of this theorem implies that the convergence radius R of the series:
j = 0 X ( j ) z j = X ˜ ( z )
is greater than 1 (i.e., R > 1 ). Therefore, there exists i 0 , where 1 i 0 n , which makes the convergence radius R i 0 of the series:
j = 0 x i 0 ( j ) z j = x ˜ i 0 ( z )
also be greater than 1 (i.e., R i 0 > 1 ). Thus, by using the Cauchy–Hadamard Theorem, we obtain:
R i 0 = lim j sup x i 0 ( j ) j > 1 .
Consequently, lim j sup x i 0 ( j ) = . This, however, implies that x will be never bounded and hence (5) is not stable.  □
Remark 2.
If δ 1 = δ 2 = = δ n = δ , then Theorem 3 becomes similar to Theorem 1 which is equivalent, as it is known, to Theorem 2.
Corollary 1.
Suppose that A is a triangular matrix with diagonal elements ξ i , i = 1 , , n . If 2 δ i < ξ i < 0 , i , then the zero solution of system (5) is asymptotically stable. Furthermore, such solution is not stable if either ξ i > 0 or ξ i < 2 δ i , for some i.
Proof. 
Consider the first part of Theorem 3 and A as assumed here. This will turn (9) into the form:
i = 1 n z ( 1 z 1 ) δ i ξ i = 0 .
It means that there exists i, where 0 i n , such that:
z ( 1 z 1 ) δ i ξ i = 0 .
Now, according to the assumption that supposes all roots of (9) lie inside the unit disk, we deduce that all ξ i ’s belong to the set z ( 1 z 1 ) δ i , z C and z < 1 , for 0 i n . Based on the proof of Theorem 1.4 described in [17], we also deduce:
z ( 1 z 1 ) δ i , z C and z < 1 = z C : z < 2 cos arg z π 2 δ δ i and arg z > δ i π 2 ,
where ξ i R , 0 i n . This means that 2 δ i < ξ i < 0 , as desired.  □
As a matter of fact, the above result, represented by Corollary 1, is deemed one of the main significant results in this work. It can be easily implemented for exploring the stability of some linear incommensurate FoDSs which involve just a triangular matrix A in their forms. On the contrary, one can find it extremely hard to verify condition (9) in the proposed result represented by Theorem 3. Actually, this condition is related to the full matrix that might be found in the linear incommensurate FoDSs. To deal with this problem, we present below another more practical result which is equivalent to Theorem 3.
Theorem 4.
Consider 0 < δ i < 1 , for i = 1 , 2 , , n , and M is the lowest common multiple of the denominators u i of δ i ’s, in which δ i = v i u i and ( u i , v i ) = 1 , where u i , v i Z + , i . Then, the zero solution of system (5), subject to the initial vector condition X ( 0 ) = X 0 R n , is:
  • Asymptotically stable if any zero solution of the polynomial:
    det d i a g ξ M δ 1 , ξ M δ 2 , , ξ M δ n ( 1 ξ M ) A
    lies inside the set
    C K γ ,
    where γ = 1 M and where
    K γ = z C : z 2 cos arg z γ γ and arg z γ π 2 .
  • Not stable, furthermore, if there is a zero, say ξ, of (22) such that ξ I n t K γ .
Proof. 
In accordance with Theorem 3, system (5) is asymptotically stable if all zeros of the following characteristic equation:
det d i a g z ( 1 z 1 ) δ 1 , z ( 1 z 1 ) δ 2 , , z ( 1 z 1 ) δ n A = 0
are located inside the unit circle. However, setting
1 1 z = ξ M z = 1 1 ξ M , ξ M 1 ,
will turn the above characteristic equation to be in the form:
det d i a g ξ M δ 1 1 ξ M , ξ M δ 2 1 ξ M , , ξ M δ n 1 ξ M A = 0 .
Multiplying both sides of (24) by ( 1 ξ M ) n yields:
det d i a g ξ M δ 1 , ξ M δ 2 , , ξ M δ n ( 1 ξ M ) A = 0 .
Now, one finds that it is necessary to prove the two assertions z z C : 0 < z < 1 ξ C K γ and z z C : z > 1 ξ I n t K γ . For achieving those goals, consider the following steps:
  • Step 1: (Defining the stability boundary). Consider the following curve:
    L γ = 1 1 z γ : z = 1 ,
    which defines the stability boundary for system (5) and also describes its structure. Suppose z = e i ϕ and 1 z 1 = r e i ω , for 0 ϕ < 2 π and r ( ϕ ) = r 0 , and also suppose ω ( ϕ ) = ω , where ω [ 0 , 2 π ) . Then,
    1 e i ϕ = r e i ω .
    This equation, after the imaginary and real parts are equated, will be turned into the following two components:
    sin ϕ = r sin ω , 1 cos ϕ = r cos ω .
    Observe that, when ϕ = 0 , then r = 0 . Otherwise, we have:
    tan ω = sin ϕ 1 cos ϕ .
    In view of the fact that:
    sin ϕ 1 cos ϕ = 2 sin ( ϕ / 2 ) cos ( ϕ / 2 ) 2 sin 2 ( ϕ / 2 ) = cot ϕ 2 = tan ( π 2 ϕ 2 ) ,
    we can write r and ω as r = 2 sin ϕ 2 and ω = π / 2 ϕ / 2 , respectively. From here, we obtain:
    L γ = 2 sin ϕ 2 γ exp ( i ( π ϕ ) γ 2 ) : 0 ϕ < 2 π .
    Observe that setting ψ = ω = ϕ / 2 π / 2 will turn L γ to be as follows:
    L γ = 2 cos ψ γ exp ( i γ ψ ) : π 2 ψ < π 2 .
    One can use the polar form represented by z = 2 cos ψ γ ), where arg z = γ ψ , to obtain:
    L γ = z C : z = 2 cos arg z γ γ and arg z γ π 2 .
  • Step 2: (Showing that w γ ( z ) = ( 1 1 z ) γ maps the set D E = z C : z > 1 onto I n t K γ , with noting that w γ ( 2 ) = 2 γ I n t K γ ). In view of the Open Mapping Theorem, and since w γ is nonconstant holomorphic on D E , then it maps D E to an open set. In other words, we have a neighborhood of w γ ( z * ) included in w γ ( D E ) , z * D E . This implies that the boundary of w γ ( D E ) can not be mapped by any point of D E . This means that w γ ( D I ) I n t K γ . Similarly, one can prove that w γ ( D I ) C K γ , where D I = z C : 0 < z < 1 . In view of w γ ( z C : z = 1 ) = L γ (see the previous step), and also in view of the continuity of w γ , the above arguments imply that w γ ( D E ) = K γ .
  • Step 3: For the purpose of showing the other part this theorem, we first assume that there is a solution ξ of (22) with ξ I n t K γ . This implies that z = 1 1 ξ M is a solution of (9) with z > 1 . Thus, we can deduce, in view of Theorem 3, that there is an instability of the zero solution of system (5). On the other hand, if each solution of (22) belongs to C K γ , then all solutions of (9) will belong to z C : z < 1 , which makes the zero solution of system (5), via Theorem 3, asymptotically stable, as required.
 □

4. Numerical Simulations

To highlight the primary outcomes of this paper in exploring the stability of the linear incommensurate FoDSs, Theorem 4 will be utilized to investigate two examples that explore such stability when these systems involve full matrices in their forms.
Example 1.
Consider the following linear incommensurate FoDS:
C Δ 0 1 2 x 1 ( j + 1 1 2 ) C Δ 0 1 4 x 2 ( j + 1 1 4 ) = 1 1 9 16 1 2 x 1 ( j ) x 2 ( j ) , j = 0 , 1 , 2 ,
One can obtain M to be equal 4. This, however, implies:
det ξ 2 0 0 ξ ( 1 ξ 4 ) 1 1 9 16 1 2 = 0
det ξ 2 + ( 1 ξ 4 ) ( 1 ξ 4 ) 9 16 ( 1 ξ 4 ) ξ 1 2 ( 1 ξ 4 ) = 0
1 16 ξ 8 + 1 2 ξ 6 ξ 5 1 8 ξ 4 + ξ 3 1 2 ξ 2 + ξ + 1 16 = 0
Accordingly, the solution of (29) will be in the following form:
ξ ¯ = 1.1634 6.0451 × 10 2 0.78732 + 3.1894 i 0.78732 3.1894 i 1.3269 0.4875 i 1.3269 + 0.4875 i 7.2415 × 10 2 + 0.80874 i 7.2415 × 10 2 0.80874 i
As per Theorem 4, and due to ξ i C K 1 4 , i = 1 , 2 , , 8 , then system (28) is asymptotically stable about its zero solution.
In order to demonstrate the validity of the obtained outcomes, one can observe, from Figure 1, that the two states of system (28) converge to zero, and hence it is stable.
Example 2.
Consider the following linear incommensurate FoDS:
C Δ 0 1 2 x 1 ( t + 1 1 2 ) C Δ 0 1 3 x 2 ( t + 1 1 3 ) C Δ 0 2 3 x 3 ( t + 1 2 3 ) = 1 0 0.2 3.4 1 0.2 0 0 1 x 1 ( t ) x 2 ( t ) x 3 ( t )
One can find M = 6 , which leads to:
det ξ 3 0 0 0 ξ 2 0 0 0 ξ 4 ( 1 ξ 6 ) 1 0 0.2 3.4 1 0.2 0 0 1 = 0
ξ 18 + ξ 16 + ξ 15 + ξ 14 ξ 13 + 2 ξ 12 ξ 11 2 ξ 10 ξ 9 2 ξ 8 + ξ 7 2 ξ 6 + ξ 5 + ξ 4 + ξ 3 + ξ 2 + 1 = 0
Consequently, the solution of (32) will be in the form:
ξ ¯ = 0.32131 0.87498 i 0.32131 + 0.87498 i 0.32131 0.87498 i 0.32131 + 0.87498 i 0.85180 0.54463 + 0.7276 i 0.54463 0.7276 i 0.54463 + 0.7276 i 0.54463 0.7276 i 0.42590 0.73768 i 0.42590 + 0.73768 i 1.1510 1.1510 1.2106 1.2106 0.58699 1.0167 i 0.58699 + 1.0167 i 1.1740
where
K 1 6 = z C : z 2 cos 6 arg z 1 6   a n d   arg z π 12
Hence, in view of Theorem 4, system (31) is also asymptotically stable about its zero solution.
To confirm the final inference of Example 2, Figure 2 illustrates such stability by exhibiting the convergence of all system’s states to zero, which shows the validity of the proposed results.

5. Conclusions and Future Works

In the present work, some simple applicable conditions for judging the stability of the linear incommensurate Fractional-order Difference Systems have been reported as novel results. These results have really been verified numerically by illustrating the stability of the solutions of such systems via several examples. All results of this work are applicable to be implemented in lots of difference systems, like e.g., the Duffing oscillator system which has been successfully employed to model a set of physical schemes such as beam buckling, ionization waves in plasmas, nonlinear electronic circuits, stiffening springs, and superconducting Josephson parametric amplifiers. Such investigation together with studying the dynamics of the linear incommensurate FoDSs will be some of several targets that are left for future consideration.

Author Contributions

Conceptualization, A.O.; Data curation, N.D.; Formal analysis, N.D.; Funding acquisition, G.G.; Investigation, N.D. and A.O.; Methodology, A.O.; Resources, I.M.B. and Viet-Thanh Pham; Software, I.M.B.; Validation, G.G.; Visualization, V.-T.P.; Writing—original draft, I.M.B. and V.-T.P.; Writing—review & editing, G.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sun, H.; Chen, W.; Chen, Y. Variable–order fractional differential operators in anomalous diffusion modeling. Phys. Stat. Mech. Its Appl. 2009, 388, 4586–4592. [Google Scholar] [CrossRef]
  2. Wang, Y.; Gu, L.; Xu, Y.; Cao, X. Practical tracking control of robot manipulators with continuous fractional-order nonsingular terminal sliding mode. IEEE Trans. Ind. Electron. 2016, 63, 6194–6204. [Google Scholar] [CrossRef]
  3. Machado, J.; Mata, M.E. A fractional perspective to the bond graph modeling of world economies. Nonlinear Dyn. 2015, 80, 1839–1852. [Google Scholar] [CrossRef] [Green Version]
  4. Mozyrska, D.; Wyrwas, M. Stability analysis of discrete-time fractional-order Caputo-type Lorenz system. In Proceedings of the 2019 42nd International Conference on Telecommunications and Signal Processing (TSP), Budapest, Hungary, 3–5 July 2019. [Google Scholar]
  5. May, R.M. Simple mathematical models with very complicated dynamics. Nature 1976, 261, 459–467. [Google Scholar] [CrossRef] [PubMed]
  6. Diaz, J.B.; Olser, T.J. Differences of fractional order. Math. Comput. 1974, 28, 185–202. [Google Scholar] [CrossRef]
  7. Srivastava, H.M.; Owa, S. Univalent Functions, Fractional Calculus, and Their Applications; Ellis Horwood Ltd.: Hertfordshire, UK, 1989. [Google Scholar]
  8. Jarad, F.; Taş, K. On Sumudu transform method in discrete fractional calculus. Abstr. Appl. Anal. 2012, 2012, 270106. [Google Scholar] [CrossRef] [Green Version]
  9. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific Publishing Company: Singapore, 2000. [Google Scholar]
  10. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies, 204; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  11. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  12. Mozyrska, D.; Wyrwas, M. Solutions of fractional linear difference systems with Riemann–Liouville–type operator via transform method. In Proceedings of the ICFDA’14 International Conference on Fractional Differentiation and Its Applications, Catania, Italy, 23–25 June 2014. [Google Scholar]
  13. Mozyrska, D.; Ostalczyk, P.; Wyrwas, M. Stability conditions for fractional-order linear equations with delays. Bull. Pol. Acad. Sci. Tech. Sci. 2018, 66, 449–454. [Google Scholar]
  14. Mozyrska, D.; Wyrwas, M. The Z-transform method and delta type fractional difference operators. Discret. Dyn. Nat. Soc. 2015, 2015, 852734. [Google Scholar] [CrossRef] [Green Version]
  15. Abu-Saris, R.; Al-Mdallal, Q. On the asymptotic stability of linear system of fractional-order difference equations. Int. J. Theory Appl. Fract. Calc. Appl. Anal. 2013, 16, 613–629. [Google Scholar] [CrossRef]
  16. Mozyrska, D.; Wyrwas, M. Explicit criteria for stability of two-dimensional fractional difference systems. Int. J. Dynam. Control 2017, 5, 4–9. [Google Scholar] [CrossRef] [Green Version]
  17. Čermák, J.; Gyori, I.; Nechvátal, L. On explicit stability conditions for a linear farctional difference system. Fract. Calc. Appl. Anal. 2015, 18, 651–672. [Google Scholar] [CrossRef]
  18. Atici, F.M.; Eloe, P.W. Linear systems of fractional nabla difference equations. Rocky Mt. J. Math. 2011, 41, 353–370. [Google Scholar] [CrossRef]
  19. Chen, F. Fixed points and asymptotic stability of nonlinear fractional difference equations. El. J. Qualit. Theory Differ. Equ. 2011, 39, 1–18. [Google Scholar] [CrossRef]
  20. Abdeljawad, T. On Riemann and Caputo fractional differences. Comput. Math. Appl. 2011, 62, 1602–1611. [Google Scholar] [CrossRef] [Green Version]
  21. Elaydi, S.; Murakami, S. Asymptotic stability versus exponential stability in linear volterra difference equations of convolution type. J. Differ. Equ. Appl. 1996, 940, 35–46. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The stability of the zero solution of system (28).
Figure 1. The stability of the zero solution of system (28).
Mathematics 08 01754 g001
Figure 2. The stability of the zero solution of system (31).
Figure 2. The stability of the zero solution of system (31).
Mathematics 08 01754 g002

Share and Cite

MDPI and ACS Style

Djenina, N.; Ouannas, A.; Batiha, I.M.; Grassi, G.; Pham, V.-T. On the Stability of Linear Incommensurate Fractional-Order Difference Systems. Mathematics 2020, 8, 1754. https://doi.org/10.3390/math8101754

AMA Style

Djenina N, Ouannas A, Batiha IM, Grassi G, Pham V-T. On the Stability of Linear Incommensurate Fractional-Order Difference Systems. Mathematics. 2020; 8(10):1754. https://doi.org/10.3390/math8101754

Chicago/Turabian Style

Djenina, Noureddine, Adel Ouannas, Iqbal M. Batiha, Giuseppe Grassi, and Viet-Thanh Pham. 2020. "On the Stability of Linear Incommensurate Fractional-Order Difference Systems" Mathematics 8, no. 10: 1754. https://doi.org/10.3390/math8101754

APA Style

Djenina, N., Ouannas, A., Batiha, I. M., Grassi, G., & Pham, V. -T. (2020). On the Stability of Linear Incommensurate Fractional-Order Difference Systems. Mathematics, 8(10), 1754. https://doi.org/10.3390/math8101754

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop