Next Article in Journal
Evolution Equations in Hilbert Spaces via the Lacunae Method
Next Article in Special Issue
A New Adaptive Robust Sliding Mode Control Approach for Nonlinear Singular Fractional-Order Systems
Previous Article in Journal
The New Wave Structures to the Fractional Ion Sound and Langmuir Waves Equation in Plasma Physics
Previous Article in Special Issue
Image Enhancement Based on Rough Set and Fractional Order Differentiator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis of the Nabla Distributed-Order Nonlinear Systems

1
Department of Mathematics and Computer Science, Shanxi Normal University, Taiyuan 030006, China
2
Mechatronics, Embedded Systems and Automation (MESA) Lab, School of Engineering, University of California, Merced, 5200 North Lake Road, Merced, CA 95343, USA
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(5), 228; https://doi.org/10.3390/fractalfract6050228
Submission received: 1 February 2022 / Revised: 14 March 2022 / Accepted: 21 March 2022 / Published: 20 April 2022

Abstract

:
The stability of the nabla discrete distributed-order nonlinear dynamic systems is investigated in this paper. Firstly, a sufficient condition for the asymptotic stability of the nabla discrete distributed-order nonlinear systems is proposed based on Lyapunov direct method. In addition, some properties of the nabla distributed-order operators are derived. Based on these properties, a simpler criterion is provided to determine the stability of such systems. Finally, two examples are given to illustrate the validity of these results.

1. Introduction

Fractional calculus [1] is a generalization of integer-order calculus, which includes integration and differentiation of non-integer order. It is well known that the dynamic properties of more and more real systems can be characterized by fractional differential equations, which is due to the memory and hereditary properties of fractional calculus. Therefore, fractional calculus has been widely used in biology [2], economy [3], physics [4,5,6,7] and other fields [8,9,10,11,12].
At first, the study of fractional calculus mainly focused on the fractional operator with a constant and single-valued order. Although a constant value fractional order system can better and more accurately describe the actual physical system than an integer-order system, the constant and single-valued nature of the order still limits its ability to accurately capture certain complex phenomena whose underlying physics could either evolve in time or emerge as the result of the interplay of multiple orders. In relatively recent years, many scholars began to pay attention to the variable-order fractional operators [13,14] and the distribution-order operators [15,16]. The variable-order means that the derivative order can be a function of either dependent (e.g., state variables of the system) or independent (e.g., space or time) variables and can change value following the evolution of the system. The distributed-order definition of the operator allows considering an accumulation of orders and accounting for, as an example, physical phenomena such as memory effects in composite materials [17] or multi-scale effects [18]. A typical example that illustrates the capabilities of this class of operators is the mechanical behavior of viscoelastic materials having spatially varying properties [19,20]. Therefore distributed-order fractional calculus is a natural generalization of fractional calculus, which provides a tool to model more complex systems since it both retains the non-local and memory properties of the constant and single-valued order fractional calculus, and allows the case with multiple coexisting orders. Therefore, distributed-order systems can be applied in viscoelasticity, transportation process, robots and so on [21,22,23,24,25,26,27,28,29].
Discrete fractional calculus has also attracted great interest of researchers since most practical work in terms of discrete time series. There are many kinds of definitions for discrete fractional calculus, such as delta operator [30], nabla operator [31], etc. Among them, nabla discrete fractional calculus is especially noteworthy. The time–domain response and the infinite-dimensional nature of nabla discrete fractional-order systems are investigated in [32] and in [33,34], respectively. Some Lyapunov inequalities about nabla discrete fractional difference are presented in [35]. However, to our best knowledge, there are few results related to discrete distributed-order systems.
Stability analysis of dynamic systems is the most important and basic problem. The Lyapunov direct method is a powerful tool for analyzing the stability of a system since it does not need to solve explicit solutions of differential systems. Since then, many scholars have continuously improved and promoted the Lyapunov direct method. Li et al. [36] propose the fractional Lyapunov direct method to study the Mittag–Leffler stability of fractional-order nonlinear systems. Analogously, the definition of discrete Mittag–Leffler stability is given and the Lyapunov direct method is extended to study nabla discrete fractional-order systems in [37]. Recently, Fernández-Anaya et al. [38] generalize some properties of the Caputo fractional derivative to distributed-order cases and establish some stability conditions for such systems based on the distributed-order Lyapunov direct method. It may be noted that there are few references about the stability analysis of nabla discrete distributed-order nonlinear systems.
Motivated by the previous discussions, this paper mainly investigates the stability of the nabla discrete distributed-order nonlinear dynamic systems. We will establish a stability condition for the nabla discrete distributed-order nonlinear systems based on Lyapunov direct method and comparison principle. Additionally, in virtue of some important inequalities, we will propose a simpler stability criterion.
The rest of the paper is organized as follows. In Section 2, some basic definitions and preliminary knowledge are introduced. In Section 3, we propose some stability criteria for the nabla discrete distributed-order nonlinear systems. Moreover, some properties of the nabla distributed-order operators are derived. In Section 4, we illustrate the effectiveness of the results proposed in this paper via two examples. Finally, the main conclusions are drawn in Section 5.
Notations: R denotes the set of real numbers, R n denotes the n-dimensional Euclidean space and Z + denotes the set of positive integers. Given a real number a, we define set N a + 1 as N a + 1 : = { a + 1 , a + 2 , a + 3 , } . Function Γ ( z ) = 0 + e t t z 1 d t is the Gamma function.

2. Preliminaries

In this section, we firstly introduce some definitions, corresponding properties and lemmas on nabla fractional calculus.
Definition 1
([31]). The n-th order nabla backward difference of a function f : N a + 1 n R is defined by
n f ( k ) = j = 0 n ( 1 ) j n j f ( k j ) ,
where n Z + , k N a + 1 , a R and n j = Γ ( n + 1 ) Γ ( j + 1 ) Γ ( n j + 1 ) .
Definition 2
([31]). The α-th order nabla fractional sum of a function f : N a + 1 R is defined by
a k α f ( k ) = j = 0 k a 1 ( 1 ) j α j f ( k j ) ,
where k N a + 1 , a R and α ( n 1 , n ) , n Z + , and α j = ( 1 ) j Γ ( α + j ) Γ ( j + 1 ) Γ ( α ) .
Definition 3
([31]). The nabla Caputo and Riemann–Liouville fractional differences of a function f : N a + 1 n R are defined, respectively, by
a C k α f ( k ) = a k ( n α ) n f ( k ) = a k α n n f ( k ) ,
and
a R L k α f ( k ) = n a k ( n α ) f ( k ) = n a k α n f ( k ) ,
where k N a + 1 , a R and α ( n 1 , n ) , n Z + .
In this paper, we mainly adopt the Caputo definition method.
Definition 4.
The nabla distributed-order difference of a function f : N a + 1 n R is defined by
a C k c ( α ) f ( k ) = n 1 n c ( α ) a C k α f ( k ) d α ,
where c ( α ) denotes the weight function, c ( α ) 0 and α ( n 1 , n ) , n Z + .
Based on Definition 4, the nabla distributed-order difference of a constant function can be easily obtained as follows.
Property 1.
For 0 < α < 1 , the nabla distributed-order difference of a constant function f ( k ) = d ( k N a + 1 n ) is
a C k c ( α ) d = 0 .
Proof. 
Using the definition of the nabla Caputo fractional difference and since d is a constant, one has
a C k α d = a k ( 1 α ) 1 d = a k ( 1 α ) 0 .
Then, according to Definition 4, it is easily known that the nabla distributed-order difference of a constant function is also 0. □
Like Laplace transform in the standard calculus, the N -transform is a powerful tool to analyze properties of the nabla fractional calculus, so the definition and key properties on N -transform of a function are given as follows.
Definition 5
([31]). The N -transform of a function f : N a + 1 R is defined by
N a { f } ( s ) = k = 1 + ( 1 s ) k 1 f ( k + a ) .
Lemma 1
([33]). Let α ( n 1 , n ) and n Z + . If the N -transform of a function f : N a + 1 n R converges for | s 1 | < ρ for some ρ > 0 , then
N a a C k α f ( s ) = s α N a f ( s ) i = 0 n 1 s α i 1 i f ( k ) | k = a ,
for | s 1 | < ρ .
Based on the definition of the nabla distributed-order difference and Lemma 1, we can deduce that
N a a C k c ( α ) f ( s ) = n 1 n c ( α ) [ s α N a f ( s ) i = 0 n 1 s α i 1 i f ( k ) | k = a ] d α = n 1 n c ( α ) s α N a f ( s ) d α n 1 n c ( α ) [ i = 0 n 1 s α i 1 i f ( k ) | k = a ] d α = C ( s ) F ( s ) i = 0 n 1 1 s i + 1 C ( s ) i f ( k ) | k = a ] ,
where C ( s ) = n 1 n c ( α ) s α d α and F ( s ) = N a f ( s ) .
Lemma 2
([32]). For the discrete Mittag–Leffler function defined by F α , β ( λ , k , a ) = j = 0 + λ j Γ ( j α + β ) ( k a ) j α + β 1 ¯ , its N -transform is
N a { F α , β ( λ , k , a ) } ( s ) = s α β s α λ ,
for | s 1 | < 1 and | s α | > | λ | , where ( k a ) j α + β 1 ¯ = Γ ( k a + j α + β 1 ) Γ ( k a ) .
The Convolution Theorem builds the bridge between the time domain and the frequency domain, In the following, we will give the Convolution Theorem on nabla fractional calculus.
Lemma 3
([31]). (Convolution Theorem) Let functions f , g : N a + 1 R . Then N a { f g } ( s ) = N a { f } ( s ) N a { g } ( s ) , where * denotes the convolution operation, i.e.,
f g = j = a + 1 k f ( k j + a + 1 ) g ( j ) .
Next, we provide the Final Value Theorem on the N -transform, which plays a key role in the stability analysis.
Lemma 4
([32]). (Final Value Theorem) Let a function f : N a + 1 R and a R . If F ( s ) = N a f ( s ) and the poles of s F ( s ) satisfy | s 1 | > 1 , then
lim k + f ( k ) = lim s 0 s F ( s ) .

3. Stability Analysis

Consider the following nabla discrete distributed-order nonlinear systems
a C k c ( α ) x ( k ) = f ( x ( k ) , k ) ,
where α ( 0 , 1 ) , k N a + 1 , x ( k ) Ω , where Ω R n is a domain containing the origin, f : Ω × N a + 1 R n and x ( a ) is initial condition.
Throughout this paper, we assumed that 0 1 c ( α ) s α d α 0 and the solution of system (1) is existent and unique. Accordingly, the form of the unique solution of system (1) is given as follows.
Theorem 1.
Consider the nabla discrete distributed-order nonlinear system (1), its unique solution has the following form:
x ( k ) = x ( a ) + f N a 1 1 0 1 c ( α ) s α d α .
Proof. 
Taking the N -transform on both sides of Equation (1) yields
C ( s ) X ( s ) 1 s C ( s ) x ( a ) = F ( s ) ,
then we have
X ( s ) = 1 s x ( a ) + F ( s ) C ( s ) .
Taking the inverse N -transform to Equation (2) yields
x ( k ) = u ( k a 1 ) x ( a ) + f N a 1 1 C ( s ) .
Due to the discrete-time unit, step function is
u ( n ) = 1 n 0 , 0 n < 0 ,
and k N a + 1 , then u ( k a 1 ) = 1 .
Therefore,
x ( k ) = x ( a ) + f N a 1 1 0 1 c ( α ) s α d α .
Before discussing the stability problem, we need to present the definition of the equilibrium point of the nabla discrete distributed-order nonlinear system.
Definition 6.
For the nabla discrete distributed-order nonlinear system (1), the constant vector x ¯ is called its equilibrium point if f ( x ¯ , k ) = 0 .
Remark 1.
For convenience, we often assume that the equilibrium point is at the origin, i.e., x ¯ = 0 . This is no loss of generality, since if the equilibrium point is x ¯ 0 , we can take the change of variable y ( k ) = x ( k ) x ¯ , then
a C k c ( α ) y ( k ) = a C k c ( α ) ( x ( k ) x ¯ ) = f ( x ( k ) , k ) = f ( y ( k ) + x ¯ , k ) = g ( y ( k ) , k ) ,
where g ( 0 , k ) = 0 and the equilibrium point of the system about the new variable y ( k ) is at the origin.
In the following, we will investigate the stability of the nabla discrete distributed-order nonlinear system (1), the relevant stability definitions are firstly introduced.
Definition 7.
Let x ¯ = 0 is an equilibrium point of the nabla discrete distributed-order nonlinear system (1) and k = a is the initial time, the equilibrium point x ¯ is stable if for every ε > 0 , there exists δ > 0 such that x ( k ) < ε holds for all k a ( k N a + 1 ) , when x ( a ) < δ .
Definition 8.
The equilibrium point x ¯ = 0 of the nabla discrete distributed-order nonlinear system (1) is asymptotically stable if it is stable and δ can be chosen such that if x ( a ) < δ , then lim k + x ( k ) = 0 .
Remark 2.
The system stability we stated in the following refers to the stability of system at the equilibrium point x ¯ = 0 .
The following Lemma extends the comparison principle for the nabla fractional difference in [37] to the nabla distributed-order case.
Lemma 5.
(Nabla Distributed-order Comparison Principle) Let two functions x , y : N a R . If a C k c ( α ) x ( k ) a C k c ( α ) y ( k ) , where 0 < α < 1 , N a 1 1 C ( s ) 0 ( k N a + 1 ) , C ( s ) = 0 1 c ( α ) s α d α and x ( a ) = y ( a ) , then for any k N a + 1 , one has x ( k ) y ( k ) .
Proof. 
Since a C k c ( α ) x ( k ) a C k c ( α ) y ( k ) ( k N a + 1 ) , then there is a non-negative function z ( k ) satisfying
a C k c ( α ) x ( k ) = a C k c ( α ) y ( k ) + z ( k ) .
Taking the N -transform on the both sides of Equation (3), then
C ( s ) X ( s ) 1 s C ( s ) x ( a ) = C ( s ) Y ( s ) 1 s C ( s ) y ( a ) + Z ( s ) .
Since x ( a ) = y ( a ) , then
X ( s ) = Y ( s ) + Z ( s ) C ( s ) .
Applying the inverse N -transform on both side of Equation (4) and using Lemma 3, then we have
x ( k ) = y ( k ) + z N a 1 1 C ( s ) = y ( k ) + z g = y ( k ) + j = a + 1 k z ( k j + a + 1 ) g ( j ) ,
where g N a 1 1 C ( s ) .
For all k N a + 1 , since g ( k ) 0 and z ( k ) 0 , then x ( k ) y ( k ) . □
Remark 3.
It is worth nothing that if we choose a Dirac function as the distribution function, the comparison principle on the nabla fractional difference can be obtained (see Lemma 3 in [37]).
It is well known that Lyapunov direct method is the most effective tool for analyzing the stability of nonlinear systems. Therefore, in the following we will generalize Lyapunov direct method to establish the stability conditions for the nabla discrete distributed-order nonlinear system.
Theorem 2.
Consider the nabla discrete distributed-order nonlinear system (1), suppose there exists a Lyapunov function V ( x ( k ) , k ) : Ω × N a + 1 R such that
γ 1 x ( k ) b V ( x ( k ) , k ) γ 2 x ( k ) b c ,
and
a C k c ( α ) V ( x ( k ) , k ) γ 3 x ( k ) b c ,
where α ( 0 , 1 ) , k N a + 1 and b, c, γ i > 0   ( i = 1 , 2 , 3 ) . If the roots of C ( s ) + γ 3 γ 2 = 0 satisfy | s 1 | > 1 , N a 1 1 ( C ( s ) + γ 3 γ 2 ) 0 , C ( s ) = 0 1 c ( α ) s α d α , then this system is asymptotically stable.
Proof. 
Given ε > 0 , there must exist δ > 0 such that γ 2 δ b c < γ 1 ε b .
Based on the inequality (6), one has
a C k c ( α ) V ( x ( k ) , k ) 0 .
Applying Property 1 and Lemma 5, we can obtain
V ( x ( k ) , k ) V ( x ( a ) , a ) .
Therefore, when x ( a ) < δ , we get
γ 1 x ( k ) b V ( x ( k ) , k ) V ( x ( a ) , a ) γ 2 x ( a ) b c γ 2 δ b c γ 1 ε b ,
and thus x ( k ) < ε , which implies the stability of system (1).
Based on inequalities (5) and (6), we have
a C k c ( α ) V ( x ( k ) , k ) γ 3 γ 2 V ( x ( k ) , k ) ,
then there exists a non-negative function z ( k ) such that
a C k c ( α ) V ( x ( k ) , k ) + z ( k ) = γ 3 γ 2 V ( x ( k ) , k ) .
Taking the N -transform on both sides of Equation (7) results in
C ( s ) V ( s ) C ( s ) s V ( a ) + Z ( s ) = γ 3 γ 2 V ( s ) ,
then, we can find that
V ( s ) = C ( s ) V ( a ) s ( C ( s ) + γ 3 γ 2 ) Z ( s ) ( C ( s ) + γ 3 γ 2 ) .
Applying the inverse N -transform to (8) gives
V ( x ( k ) , k ) = N a 1 C ( s ) V ( a ) s ( C ( s ) + γ 3 γ 2 ) N a 1 Z ( s ) ( C ( s ) + γ 3 γ 2 ) .
The last term of Equation (9) can be rewritten as
N a 1 Z ( s ) ( C ( s ) + γ 3 γ 2 ) = z N a 1 1 ( C ( s ) + γ 3 γ 2 ) = z g = j = a + 1 k z ( k j + a + 1 ) g ( j ) ,
where g N a 1 1 ( C ( s ) + γ 3 γ 2 ) .
Considering that g 0 and z ( k ) 0 , then
N a 1 Z ( s ) ( C ( s ) + γ 3 γ 2 ) 0 .
Therefore, from Equation (9), we have
V ( x ( k ) , k ) N a 1 C ( s ) V ( a ) s ( C ( s ) + γ 3 γ 2 ) .
Since we assume that the roots of C ( s ) + γ 3 γ 2 = 0 satisfy | s 1 | > 1 , and according to Lemma 4, we can deduce that
lim k + N a 1 C ( s ) V ( a ) s ( C ( s ) + γ 3 γ 2 ) = lim s 0 C ( s ) V ( a ) ( C ( s ) + γ 3 γ 2 ) = 0 .
Combining (10), we get
lim k + V ( x ( k ) , k ) lim k + N a 1 C ( s ) V ( a ) s ( C ( s ) + γ 3 γ 2 ) = 0 .
It follows from (5) that
lim k + γ 1 x ( k ) b 0
Since b > 0 and γ 1 > 0 , then lim k + x ( k ) = 0 , which indicates that system (1) is asymptotically stable. □
Remark 4.
Note that if we choose an appropriate Dirac function as the distribution function, then we can get the asymptotic stability conditions of the nabla discrete fractional-order nonlinear systems (see Theorem 2 in [37]).
An important inequality for stability analysis based on the Lyapunov method is stated in following lemma, which will be generalized to the nabla distributed-order case.
Lemma 6
([35]). The following inequality holds
a C k α x 2 ( k ) 2 x ( k ) a C k α x ( k ) ,
for α ( 0 , 1 ) , k N a + 1 , x ( k ) R and a R .
Lemma 7.
Let x : N a + 1 R . Then, the following inequality holds
a C k c ( α ) x 2 ( k ) 2 x ( k ) a C k c ( α ) x ( k ) ,
where α ( 0 , 1 ) , k N a + 1 and a R .
Proof. 
Based on Lemma 6, one obtains
a C k α x 2 ( k ) 2 x ( k ) a C k α x ( k ) .
Due to c ( α ) 0 , then we have
c ( α ) a C k α x 2 ( k ) 2 c ( α ) x ( k ) a C k α x ( k ) .
Integrating both sides of the inequality (11) with respect to α from 0 to 1, we have
0 1 c ( α ) a C k α x 2 ( k ) d α 2 0 1 c ( α ) x ( k ) a C k α x ( k ) d α .
This completes the proof. □
Lemma 8.
For x : N a + 1 R n , Lemma 7 still holds, that is,
a C k c ( α ) x T ( k ) x ( k ) 2 x T ( k ) a C k c ( α ) x ( k ) ,
where α ( 0 , 1 ) .
Proof. 
By decomposing the inequality (12) into the sum of scalar products and applying Lemma 7, the result is obvious. □
Based on Lemmas 5 and 8 and Property 1, we will provide a simpler method to analyze the stability of the nabla discrete distributed-order nonlinear system (1).
Theorem 3.
Consider the nabla discrete distributed-order nonlinear system (1).
(i)
If x T ( k ) f ( x ( k ) , k ) 0 and N a 1 1 C ( s ) 0 , where C ( s ) = 0 1 c ( α ) s α d α , then system (1) is stable.
(ii)
Let ξ > 0 , if x T ( k ) f ( x ( k ) , k ) ξ x ( k ) 2 , the roots of C ( s ) + ξ = 0 satisfy | s 1 | > 1 , N a 1 1 ( C ( s ) + ξ 0 , then system (1) is asymptotically stable.
Proof. 
Choose the Lyapunov function
V ( x ( k ) , k ) = 1 2 x T ( k ) x ( k ) .
( i ) Using Lemma 8 to Equation (13) yields that
a C k c ( α ) V ( x ( k ) , k ) x T ( k ) a C k c ( α ) x ( k ) .
Since x T ( k ) f ( x ( k ) , k ) 0 and note that a C k c ( α ) x ( k ) = f ( x ( k ) , k ) , then
a C k c ( α ) V ( x ( k ) , k ) x T ( k ) f ( x ( k ) , k ) 0 .
Based on Property 1 and Lemma 5, one has
V ( x ( k ) , k ) V ( x ( a ) , a ) .
In terms of the definition of the function V ( x ( k ) , k ) , we obtain
x T ( k ) x ( k ) x T ( a ) x ( a ) .
Given ε > 0 and choose δ = ε , then when x ( a ) < δ , we can deduce that x ( k ) < δ for all k N a + 1 , which indicates x ( k ) < ε for all k N a + 1 , this means that system (1) is stable.
( ii ) If x T ( k ) f ( x ( k ) , k ) ξ x ( k ) 2 , then from the inequality (14), we have
a C k c ( α ) V ( x ( k ) , k ) ξ x ( k ) 2 .
Since the Lyapunov function
V ( x ( k ) , k ) = 1 2 x T ( k ) x ( k ) = 1 2 x ( k ) 2 ,
then
1 2 x ( k ) 2 V ( x ( k ) , k ) x ( k ) 2 .
According to Theorem 2, the asymptotic stability of system (1) can be obtained. □

4. Numerical Examples

Example 1.
Consider the following nabla discrete distributed-order nonlinear system
a C k c ( α ) x 1 ( k ) = x 1 ( k ) x 1 ( k ) x 2 2 ( k ) , a C k c ( α ) x 2 ( k ) = x 2 ( k ) + x 1 2 ( k ) x 2 ( k ) ,
where α ( 0 , 1 ) and c ( α ) = δ ( α 1 3 ) .
Choosing the Lyapunov function V ( x ( k ) , k ) = 1 2 ( x 1 2 ( k ) + x 2 2 ( k ) ) , then we have
a C k c ( α ) V ( x ( k ) , k ) = 1 2 a C k c ( α ) ( x 1 2 ( k ) + x 2 2 ( k ) ) x 1 ( k ) a C k c ( α ) x 1 ( k ) + x 2 ( k ) a C k c ( α ) x 2 ( k ) = ( x 1 2 ( k ) + x 2 2 ( k ) ) = x ( k ) 2 .
Due to 1 2 x ( k ) 2 V ( x ( k ) , k ) x ( k ) 2 , then one has γ 1 = 1 2 , γ 2 = 1 , γ 3 = 1 . The root of
C ( s ) + γ 3 γ 2 = s 1 3 + 1 = 0
is s = 1 , which satisfies | s 1 | > 1 . It follows from Lemma 2 that
N a 1 1 ( C ( s ) + γ 3 γ 2 ) = N a 1 1 s 1 3 + 1 = F 1 3 , 1 3 ( 1 , k , a ) .
Since F 1 3 , 1 3 ( 1 , k , a ) 0 , then N a 1 1 ( C ( s ) + γ 3 γ 2 ) 0 . From Theorem 2, we conclude that system (15) is asymptotically stable.
Example 2.
Consider the following nabla discrete distributed-order nonlinear system
a C k c ( α ) x 1 ( k ) = 4 ( x 1 ( k ) + x 2 ( k ) h ( x 1 ( k ) , x 2 ( k ) , k ) ) , a C k c ( α ) x 2 ( k ) = 4 ( x 2 ( k ) x 1 ( k ) h ( x 1 ( k ) , x 2 ( k ) , k ) ) ,
where α ( 0 , 1 ) , c ( α ) = δ ( α 2 3 ) + 4 δ ( α 1 3 ) , and h ( x 1 ( k ) , x 2 ( k ) , k ) is a differentiable function.
Since
x T ( k ) f ( x ( k ) , k ) = x 1 ( k ) x 2 ( k ) 4 ( x 1 ( k ) + x 2 ( k ) h ( x 1 ( k ) , x 2 ( k ) , k ) ) 4 ( x 2 ( k ) x 1 ( k ) h ( x 1 ( k ) , x 2 ( k ) , k ) ) = 4 ( x 1 2 ( k ) + x 2 2 ( k ) ) = 4 x ( k ) 2 ,
we can let ξ = 4 in Theorem 3.
Now, the distribution function c ( α ) = δ ( α 2 3 ) + 4 δ ( α 1 3 ) , then C ( s ) = 0 1 c ( α ) s α d α = s 2 3 + 4 s 1 3 , the roots of
C ( s ) + ξ = C ( s ) + 4 = s 2 3 + 4 s 1 3 + 4 = ( s 1 3 + 2 ) 2 = 0
are s 1 = 8 , s 2 = 8 , which satisfies that | s 1 | > 1 .
It follows from Lemma 2 that
N a 1 1 s 1 3 + 2 = F 1 3 , 1 3 ( 2 , k , a ) 0 .
According to Lemma 3, we can obtain that
N a 1 1 C ( s ) + ξ = N a 1 1 C ( s ) + 4 = N a 1 1 s 1 3 + 2 N a 1 1 s 1 3 + 2 0 .
Therefore, system (16) is asymptotically stable based on Theorem 3.

5. Conclusions

In this paper, the stability of the nabla discrete distributed-order nonlinear systems have been studied. The nabla distributed-order comparison principle is introduced. We generalize the Lyapunov direct method to establish the stability condition for the nabla discrete distributed-order systems. In addition, combined with some important inequalities, a simpler stability analysis method is provided. Finally, two examples are given to illustrate the validity of the obtained results. Based on the stability results established in this paper, one can investigate the controller design problem or the performance analysis problem of the the nabla distributed-order nonlinear systems.

Author Contributions

Conceptualization, C.W. and Y.C.; methodology, C.W.; software, C.W.; validation, C.W., T.Z. and Y.C.; formal analysis, C.W.; investigation, C.W.; resources, Y.C.; writing—original draft preparation, T.Z.; writing—review and editing, C.W., T.Z. and Y.C.; supervision, C.W.; funding acquisition, C.W. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by Foundation of National Natural Science of China under Grant No. 61907027 and Natural Science Foundation of Shanxi Normal University under Grant ZR1601.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.

Acknowledgments

We would like to express our great appreciation to the editors and reviewers.

Conflicts of Interest

The authors declare that they have no conflict of interest regarding the publication of this paper.

References

  1. Podlubny, I. Fractional Differential Equations: An introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netherlands, 1998. [Google Scholar]
  2. Toledo-Hernandez, R.; Rico-Ramirez, V.; Iglesias-Silva, G.A.; Diwekar, U.M. A fractional calculus approach to the dynamic optimization of biological reactive systems. Part I: Fractional models for biological reactions. Chem. Eng. Sci. 2014, 117, 217–228. [Google Scholar] [CrossRef]
  3. Wang, Z.; Huang, X.; Shi, G.D. Analysis of nonlinear dynamics and chaos in a fractional order financial system with time delay. Comput. Math. Appl. 2011, 62, 1531–1539. [Google Scholar] [CrossRef] [Green Version]
  4. Gritsenko, D.; Paoli, R. Theoretical analysis of fractional viscoelastic flow in circular pipes: General solutions. Appl. Sci. 2020, 10, 9093. [Google Scholar] [CrossRef]
  5. Drapaca, C.S.; Sivaloganathan, S. A fractional model of continuum mechanics. J. Elast. 2012, 107, 105–123. [Google Scholar] [CrossRef]
  6. Sumelka, W. Thermoelasticity in the framework of the fractional continuum mechanics. J. Therm. Stresses. 2014, 37, 678–706. [Google Scholar] [CrossRef]
  7. Saad, K.M.; Gómez-Aguilar, J.F. Analysis of reaction–diffusion system via a new fractional derivative with non-singular kernel. Physica A. 2018, 509, 703–716. [Google Scholar] [CrossRef]
  8. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  9. Zhang, X.F.; Chen, Y.Q. Admissibility and robust stabilization of continuous linear singular fractional order systems with the fractional order α: The 0<α<1 case. ISA Trans. 2018, 82, 42–50. [Google Scholar]
  10. Li, R.C.; Zhang, X.F. Adaptive sliding mode observer design for a class of T-S fuzzy descriptor fractional order systems. IEEE Trans. Fuzzy Syst. 2020, 28, 1951–1959. [Google Scholar] [CrossRef]
  11. Zhang, X.F.; Lin, C.; Chen, Y.Q.; Boutat, D. A unified framework of stability theorems for LTI fractional order systems with 0<α<2. IEEE Trans. Circuit Syst. II-Express 2020, 67, 3237–3241. [Google Scholar]
  12. Zhang, X.F.; Huang, W.K. Adaptive neural network sliding mode control for nonlinear singular fractional order systems with mismatched uncertainties. Fractal Fract. 2020, 4, 50. [Google Scholar] [CrossRef]
  13. Wu, G.C.; Deng, Z.G.; Baleanu, D.; Zeng, D.Q. New variable-order fractional chaotic systems for fast image encryption. Chaos 2019, 29, 083103. [Google Scholar] [CrossRef]
  14. Padron, J.P.; Perez, J.P.; Pérez Díaz, J.J.; Martinez Huerta, A. Time-delay synchronization and anti-synchronization of variable-order fractional discrete-time Chen–Rossler chaotic systems using variable-order fractional discrete-time PID control. Mathematics 2021, 9, 2149. [Google Scholar] [CrossRef]
  15. Caputo, M. Elasticitá e Dissipazione; Zanichelli: Bologna, Italy, 1969. [Google Scholar]
  16. Caputo, M. Mean fractional-order-derivatives differential equations and filters. Ann. Univ. Ferrara 1995, 41, 73–84. [Google Scholar] [CrossRef]
  17. Caputo, M.; Fabrizio, M. The kernel of the distributed order fractional derivatives with an application to complex materials. Fractal Fract. 2017, 1, 13. [Google Scholar] [CrossRef] [Green Version]
  18. Calcagni, G. Towards multifractional calculus. Front. Phys. 2018, 6, 58. [Google Scholar] [CrossRef]
  19. Lorenzo, C.F.; Hartley, T.T. Variable order and distributed order fractional operators. Nonlinear Dyn. 2002, 29, 57–98. [Google Scholar] [CrossRef]
  20. Ding, W.; Patnaik, S.; Sidhardh, S.; Semperlotti, F. Applications of distributed-order fractional operators: A review. Entropy 2021, 23, 110. [Google Scholar] [CrossRef]
  21. Hartley, T.T.; Lorenzo, C.F. Fractional-order system identification based on continuous order-distributions. Signal Process. 2003, 83, 2287–2300. [Google Scholar] [CrossRef]
  22. Atanackovic, T.M.; Budincevic, M.; Pilipovic, S. On a fractional distributed-order oscillator. J. Phys. A Math. Gen. 2005, 38, 6703. [Google Scholar] [CrossRef]
  23. Chechkin, A.V.; Gonchar, V.Y.; Gorenflo, R.; Korabel, N.; Sokolov, I.M. Generalized fractional diffusion equations for accelerating subdiffusion and truncated Lévy flights. Phys. Rev. E. 2008, 78, 021111. [Google Scholar] [CrossRef]
  24. Chechkin, A.; Sokolov, I.M.; Klafter, J. Natural and Modified Forms of Distributed-Order Fractional Diffusion Equations. In Fractional Dynamics: Recent Advances; World Scientific Publishing Co., Ltd.: Singapore, 2012; pp. 107–127. [Google Scholar]
  25. Al-Refai, M.; Luchko, Y. Analysis of fractional diffusion equations of distributed order: Maximum principles and their applications. Analysis 2016, 36, 123–133. [Google Scholar] [CrossRef]
  26. Fernández-Anaya, G.; Flores-Godoy, J.J.; Lugo-Peñaloza, A.F.; Muñoz-Vega, R. Stabilization and passification of distributed-order fractional linear systems using methods of preservation. J. Frankl. Inst. 2013, 350, 2881–2900. [Google Scholar] [CrossRef]
  27. Mahmoud, G.M.; Farghaly, A.A.; Abed-Elhameed, T.M.; Aly, S.A.; Arafa, A.A. Dynamics of distributed-order hyperchaotic complex van der Pol oscillators and their synchronization and control. Eur. Phys. J. Plus. 2020, 135, 1–16. [Google Scholar] [CrossRef]
  28. Al Themairi, A.; Farghaly, A. The dynamics behavior of coupled generalized van der pol oscillator with distributed order. Math. Probl. Eng. 2020, 2020, 1–10. [Google Scholar] [CrossRef]
  29. Abdelkawy, M.A.; Lopes, A.M.; Zaky, M. Shifted fractional Jacobi spectral algorithm for solving distributed order time-fractional reaction–diffusion equations. Comput. Appl. Math. 2019, 38, 1–21. [Google Scholar] [CrossRef]
  30. Abdeljawad, T. On delta and nabla Caputo fractional differences and dual identities. Discrete Dyn. Nat. Soc. 2013, 2013, 1–12. [Google Scholar] [CrossRef]
  31. Goodrich, C.; Peterson, A.C. Discrete Fractional Calculus; Springer: Cham, Switzerland, 2015. [Google Scholar]
  32. Wei, Y.H.; Gao, Q.; Cheng, S.S.; Wong, Y. Time-domain response of nabla discrete fractional order systems. arXiv 2018, arXiv:1812.11370. [Google Scholar]
  33. Wei, Y.H.; Chen, Y.Q.; Wang, J.C.; Wang, Y. Analysis and description of the infinite-dimensional nature for nabla discrete fractional order systems. Commun. Nonlinear Sci. Numer. Simul. 2019, 72, 472–492. [Google Scholar] [CrossRef]
  34. Cabada, A.; Dimitrov, N.D.; Jonnalagadda, J.M. Non-trivial solutions of non-autonomous nabla fractional difference boundary value problems. Symmetry 2021, 13, 1101. [Google Scholar] [CrossRef]
  35. Wei, Y.H.; Chen, Y.Q.; Liu, T.Y.; Wang, Y. Lyapunov functions for nabla discrete fractional order systems. ISA Trans. 2019, 88, 82–90. [Google Scholar] [CrossRef]
  36. Li, Y.; Chen, Y.Q.; Podlubny, I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag–Leffler stability. Comput. Math. Appl. 2010, 59, 1810–1821. [Google Scholar] [CrossRef] [Green Version]
  37. Wei, Y.D.; Wei, Y.H.; Chen, Y.Q.; Wang, Y. Mittag–Leffler stability of nabla discrete fractional-order dynamic systems. Nonlinear Dyn. 2020, 101, 407–417. [Google Scholar] [CrossRef]
  38. Fernández-Anaya, G.; Nava-Antonio, G.; Jamous-Galante, J.; Muñoz-Vega, R.; Hernández-Martínez, E.G. Asymptotic stability of distributed order nonlinear dynamical systems. Commun. Nonlinear Sci. Numer. Simul. 2017, 48, 541–549. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, C.; Zhu, T.; Chen, Y. Stability Analysis of the Nabla Distributed-Order Nonlinear Systems. Fractal Fract. 2022, 6, 228. https://doi.org/10.3390/fractalfract6050228

AMA Style

Wang C, Zhu T, Chen Y. Stability Analysis of the Nabla Distributed-Order Nonlinear Systems. Fractal and Fractional. 2022; 6(5):228. https://doi.org/10.3390/fractalfract6050228

Chicago/Turabian Style

Wang, Cuihong, Tianfen Zhu, and Yangquan Chen. 2022. "Stability Analysis of the Nabla Distributed-Order Nonlinear Systems" Fractal and Fractional 6, no. 5: 228. https://doi.org/10.3390/fractalfract6050228

APA Style

Wang, C., Zhu, T., & Chen, Y. (2022). Stability Analysis of the Nabla Distributed-Order Nonlinear Systems. Fractal and Fractional, 6(5), 228. https://doi.org/10.3390/fractalfract6050228

Article Metrics

Back to TopTop