Next Article in Journal
An Implicit Difference Scheme for the Fourth-Order Nonlinear Evolution Equation with Multi-Term Riemann–Liouvile Fractional Integral Kernels
Next Article in Special Issue
Split-Plot Designs with Few Whole Plot Factors Containing Clear Effects
Previous Article in Journal
Editorial for Special Issue “Fractional Calculus Operators and the Mittag–Leffler Function”
Previous Article in Special Issue
Stability Analysis for a Fractional-Order Coupled FitzHugh–Nagumo-Type Neuronal Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Synchronization of Fractional-Order Complex Dynamical Networks with Different Structures and Parameter Uncertainties

School of Mathematics and Statistics, Sichuan University of Science and Engineering, Zigong 643000, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(8), 441; https://doi.org/10.3390/fractalfract6080441
Submission received: 6 July 2022 / Revised: 3 August 2022 / Accepted: 8 August 2022 / Published: 14 August 2022
(This article belongs to the Special Issue Advances in Fractional-Order Neural Networks, Volume II)

Abstract

:
This paper deals with the asymptotic synchronization of fractional-order complex dynamical networks with different structures and parameter uncertainties (FCDNDP). Firstly, the FCDNDP model is proposed by the Riemann–Liouville (R-L) fractional derivative. According to the property of fractional calculus and the Lyapunov direct method, an original controller is proposed to achieve the asymptotic synchronization of FCDNDP. Our controller is more adaptable and effective than those in other literature. Secondly, a sufficient condition is given for the asymptotic synchronization of FCDNDP based on the asymptotic stability theorem and the matrix inequality technique. Finally, the numerical simulations verify the effectiveness of the proposed method.

1. Introduction

Fractional calculus is a popularization of integrals and differentials of integer orders. Although the history of fractional calculus and integer calculus is not very discrepant, the results show that fractional calculus is very important for the expression of the model. In addition to the memory function of fractional calculus, the characterization of a complex system has the advantages of simple modeling, clear physical meaning of parameters and exact description [1,2,3]. Because of these advantages, it can be applied in many fields, such as physics [4], engineering [2], chemistry [5], information processing [6], secure communication [7], thermal systems [8], and robot control problems [9], etc. Since fractional-order calculus makes a tool more precise with regard to the description of memory and genetic properties of multifarious materials and processes, it is extremely important to introduce fractional calculus into complex dynamic network models.
Complex networks have been widely used in different disciplines such as physics, biology and sociology because of complex network topology and spatiotemporal evolution of systems [10]. In real life, they might be applied in neural networks [11], social networks [12], ecological networks [13], electric networks [14] and other networks. These networks share or exchange opinions with each other in a certain link based on multi-nodes, so that the complex network can be used in nature, as well as various systems in society. This has made many scholars study complex networks.
Since chaotic synchronization is proposed, the synchronization of complex networks has become a popular topic in the field of complex networks. Because the fractional-order complex network with synchronous behavior is a common and vital nonlinear phenomenon, it exists widely in nature and human society. Therefore, understanding and controlling the fractional-order synchronization of complex networks has theoretical and practical significance. In addition, the controllability of fractional complex networks has turned into a popular topic, and our ability to control natural or technological systems is reflected in our understanding and ultimate proof of them [15]. Hence, the study the synchronization problem of fractional complex networks has important implications. So far, many synchronization methods are used to research complex networks, such as global synchronization [16], adaptive synchronization [17], projection synchronization [18] and finite time synchronization [19].
However, most papers study the synchronization problem of fractional-order systems with the same structure. In practical application, due to the inevitable mismatch of parameters and functions in genuine execution, the drive system and the response system are not quite the same [20]. From an engineering viewpoint, it is truly challenging to keep the two systems consistent. Therefore, it is of extraordinary importance to research the synchronization of complex networks with different structures.
The uncertainty of parameters should be considered due to the existence of environmental disturbances when we analyze the stability of complex networks. We know that the modified model parameters are usually unable to get an accurate value. The main reason is that when the model is applied to actual engineering, it will be affected by interference from the environment, leading to the uncertainty of parameters. In the dynamic behavior analysis of a nonlinear system, the influence of such parameter uncertainties cannot be ignored, because they may damage the stability, synchronization or other characteristics of the system [21]. Thus, it is indispensable to study the system problems with uncertain parameters.
Although the system with parameter uncertainty will make the system more complex, the complex network with parameter uncertainty is closer to reality. Therefore, it has become a hot area of research for scholars. In 2008, based on the linear matrix inequalities (LMI) technique, Shen gave a sufficient condition for the boundedness ensure of neural networks with uncertain parameters [22]. In 2012, Wong derived the robust synchronization of fractional-order complex networks with uncertain parameters based on the nature of a Kronecker product and the stability of fractional-order systems by applying nonlinear control [23]. In 2013, Li designed a state estimator for the problem of fractal-order complex networks with uncertain parameters by using LMI technology and matrix singular value decomposition [24]. In 2015, combining the Lyapunov stability theorem and homeomorphic mapping theorem, Samli obtained some prime sufficient conditions for the existence uniqueness and asymptotic stability [25]. In 2016, based on the fractional Lyapunov direct method, Ding obtained a new standard for global projection synchronization of the nonuniform fractional neural networks [26]. In 2018, using the property of fractional calculus and the Lyapunov direct method, Hu revealed some new properties of fractional calculus and asymptotic stability theorems for nonautonomous fractional calculus systems with Riemann–Liouville derivatives [27]. In 2020, Udpin derived a new exponential stability standard for uncertain discrete-time neural networks based on the discrete Halanay inequality and some other inequalities techniques [28]. In 2022, based on the property of fractional calculus and the fractional direct Lyapunov method, Aadhithiyan obtained a method to realize the asymptotic synchronization of the drive and response systems in nonidentical complex networks using neoteric control [29]. Because the system contains the parameter uncertainty, the interference of the system is very large still, so how to reduce or eliminate the influence of the uncertainty is an imminent issue. Coupled with the need for synchronization of nonidentical complex networks in the industry, the study of synchronization of fractional-order complex networks with different structures is still necessary.
The primary contribution of this paper can be briefly summed up as follows: (i) we propose a new fractal-order complex dynamical network model with different structures and parameter uncertainties; (ii) we conduct synchronization analysis of FCDNDP. Based on fractional differential theory and the Lyapunov direct method, we establish the criterion to guarantee the synchronization of the drive response model; (iii) contrasted with the existing results in the literature, the method is more general and less conservative.

2. Preliminaries

Definition 1
([1]). The Caputo derivative with ξ-order of function is defined as
t 0 C D t ξ f ( t ) = 1 Γ ( m ξ ) t 0 t ( t υ ) m ξ 1 f ( m ) ( υ ) d υ ,
where Γ ( · ) is Gamma function, t 0 , m 1 < ξ < m , m Z + . Particularly, when 0 < ξ < 1 ,
t 0 C D t ξ f ( t ) = 1 Γ ( 1 ξ ) t 0 t ( t υ ) ξ f ( m ) ( υ ) d υ .
Definition 2
([1]). The Riemann–Liouville derivative with ξ-order of function is defined as
t 0 R D t ξ f ( t ) = 1 Γ ( ξ ) t 0 t f ( υ ) ( t υ ) 1 ξ d υ ,
where Γ ( · ) is Gamma function, and 0 < ξ < 1 .
Lemma 1
([30]). When f ( t ) has a continuous first derivative, then
t 0 C D t ξ ( 1 2 f T ( t ) Q f ( t ) ) f T Q t 0 C D t ξ f ( t ) ,
where ξ ( 0 , 1 ) and Q is an arbitrary n-order positive definite matrix.
Lemma 2
([31]). Let ε > 0 , for W R n , H R n and the n × n matrix A, then the following formula holds
2 W T A H ε 1 W T A A T W + ε H T H .
Property 1
([1]). When 0 < ξ < 1 and f ( t ) is a continuous function, we have
t 0 R D t ξ t 0 C D t ξ f ( t ) = f ( t ) f t 0 .
Property 2
([1]). Let C be the constant; we can obtain
t o R D t ξ C = C t t o ξ Γ ( 1 ξ ) ,
where 0 < ξ < 1 and Γ ( · ) is Gamma function.
Lemma 3
([27]). Let f ( t ) be a continuous function, then the following formula holds
t 0 R D t ξ f ( t ) f t 0 = t 0 C D t ξ f ( t ) ,
where 0 < ξ < 1 .
Lemma 4
([27]). If ξ ( 0 , 1 ) , f ( t ) is the continuously differentiable function, then
t 0 R D t ξ f T ( t ) P f ( t ) ( t t 0 ) ξ Γ ( 1 ξ ) f T ( t 0 ) P f ( t 0 ) + 2 f T ( t ) P t 0 R D t ξ ( f ( t ) f ( t 0 ) ) ,
where P R n × n is a positive definite matrix.
Lemma 5
([27]). Let y = 0 be an equilibrium point of the nonautonomous fractional differential equation t 0 R D t ξ y ( t ) = g ( t , y ) which satisfies the Lipschitz constant c > 0 , and ξ ( 0 , 1 ) . If we suppose that there is a Lyapunov function v ( t ) = V ( t , y ) and some class-K functions λ p ( p = 1 , 2 , 3 ) , which satisfy the following equations
λ 1 ( y ) v ( t ) λ 2 ( y ) ,
t o R D t ξ v t t o ξ Γ ( 1 ξ ) v t o λ 3 ( y ( t ) ) ,
where for all t t 0 , then the equilibrium point of the fractional-order equation is asymptotically stable.

Model Description

In this paper, we can consider the following drive system and response system of fractional-order complex dynamical networks, as follows [29]:
t 0 C D t ξ ϕ p ( t ) = A ϕ p ( t ) + B γ ϕ p ( t ) + c q = 1 N d p q Λ 1 ϕ q ( t ) ,
and
t 0 C D t ξ φ p ( t ) = ( E 0 + Δ E ) φ p ( t ) + ( G 0 + Δ G ) η φ p ( t ) + k q = 1 N L p q Λ 2 φ q ( t ) + ω p ( t ) ,
where 0 < ξ < 1 , ϕ p ( t ) and φ p ( t ) are the state of the p -node of the above systems, A and E 0 are the real constant matrices, and B and G 0 indicate the weight association matrices. γ ϕ p ( t ) and η φ p ( t ) express the vector-valued nonlinear functions. D = ( d p q ) N × N and L = ( l p q ) N × N are the external coupling matrices which denote the link between node p and q. If they have a direct link between node p to q, d p q > 0 and l p q > 0 ; otherwise, d p q = 0 and l p q = 0 . Λ refers to the internal coupling matrix, a link between two subsystems. ω p ( t ) expresses the control input. Δ E and Δ G are the parameter uncertainties, such that,
E = ( E 0 + Δ E ) , G = ( G 0 + Δ G ) .
The parameter uncertainties Δ E and Δ G can be given by
E 0 ( t ) = H 1 Ξ 1 K 1 , G 0 ( t ) = H 2 Ξ 2 K 2 ,
where H i and K i ( i = 1 , 2 ) are positive constant matrices, and uncertain Ξ i ( t ) ( i = 1 , 2 ) are the matrices which are unknown and satisfy the following condition
Ξ i ( t ) ( Ξ i ( t ) ) T I ,
where I is an identity matrix.
We consider that the drive system is not quite the same as the response system; then, the error vector is defined by
υ p ( t ) = φ p ( t ) ϕ p ( t ) .
The error system can be expressed by
t 0 C D t ξ υ p ( t ) = t 0 C D t ξ φ p ( t ) t 0 C D t ξ ϕ p ( t ) .
By substituting (9) and (10) into (15), we can write
t 0 C D t ξ υ p ( t ) = ( E 0 + Δ E ) φ p ( t ) + ( G 0 + Δ G ) η φ p ( t ) + k q = 1 N L p q Λ 2 φ q ( t ) A ϕ p ( t ) B γ ϕ p ( t ) c q = 1 N d p q Λ 1 ϕ q ( t ) .
One can obtain
t 0 C D t ξ υ p ( t ) = ( E 0 + Δ E ) φ p ( t ) A ϕ p ( t ) ( E 0 + Δ E ) ϕ p ( t ) + ( E 0 + Δ E ) ϕ p ( t ) + ( G 0 + Δ G ) η φ p ( t ) B γ ϕ p ( t ) ( G 0 + Δ G ) η ϕ p ( t ) + ( G 0 + Δ G ) η ϕ p ( t ) + k q = 1 N L p q Λ 2 φ q ( t ) c q = 1 N d p q Λ 1 ϕ q ( t ) k q = 1 N L p q Λ 2 ϕ q ( t ) + k q = 1 N L p q Λ 2 ϕ q ( t ) + ω p ( t ) .
Assumption A1.
For any p = 1 , 2 , · · · , N , the nonlinear function f p ( · ) satisfies the Lipschitz condition if a constant c > 0 , such that
| f p χ 1 f p χ 2 | c | χ 1 χ 2 | , χ 1 , χ 2 R .

3. Main Results

In this section, studying the asymptotic synchronization problem of FCDNDP is similar to the stability analysis of the equilibrium point error model. By applying the R L fractional derivative, we initially propose an original controller and then present the vital criteria for asymptotic stability.
Theorem 1.
Let the Assumption 1 hold and scalar 0 < ξ < 1 ; the FCDNDP is asymptotically stable if the controller of form (ii) satisfies the criteria (i):
(i) Ω < 0 ,
(ii) ω p ( t ) = ( A E ) ϕ p ( t ) + B γ φ p ( t ) G η φ p ( t ) + c q = 1 N d p q Λ 1 φ q ( t ) k q = 1 N L p q Λ 2 ϕ q ( t ) , where
Ω = P E 0 + E 0 T ( t ) P + ϵ 1 1 P H 1 H 1 T P + ϵ 1 T 1 T 1 T + ϵ 2 1 P H 2 H 2 T P + ϵ 2 T 2 T 2 T + ϵ 3 θ 1 + ϵ 3 1 P G G T P + ϵ 4 1 P B B T X + ϵ 4 θ 2 + k q = 1 N L p q Λ 2 + c q = 1 N d p q Λ 1 .
Proof. 
It follows from the error system (17) that the Lyapunov function can be designed as
V ( t ) = υ p T ( t ) P υ p ( t ) .
Computing the derivative of the (19) along the error system, we can obtain
t 0 R D t ξ V p ( t ) = t 0 R D t ξ ( υ p T ( t ) P υ p ( t ) ) .
According to Lemma 4, (20) can be composed as
t 0 R D t ξ V p ( t ) = ( t t 0 ) ξ Γ ( 1 ξ ) ( υ p T ( t 0 ) P υ p ( t 0 ) ) + 2 υ p T ( t ) P t 0 R D t ξ ( υ p ( t ) υ p ( t 0 ) ) .
Based on the Lemma 3, (21) can be composed as
t 0 R D t ξ V p ( t ) = ( t t 0 ) ξ Γ ( 1 ξ ) ( υ p T ( t 0 ) P υ p ( t 0 ) ) + 2 υ p T ( t ) P t 0 C D t ξ υ p ( t ) .
Substituting the error system (17) into the (22), we can obtain
t 0 R D t ξ V p ( t ) = ( t t 0 ) ξ Γ ( 1 ξ ) ( υ p T ( t 0 ) P υ p ( t 0 ) ) + 2 υ p T ( t ) P { ( E 0 + Δ E ) φ p ( t ) A ϕ p ( t ) ( E 0 + Δ E ) ϕ p ( t ) + ( E 0 + Δ E ) ϕ p ( t ) + ( G 0 + Δ G ) η φ p ( t ) B γ ϕ p ( t ) ( G 0 + Δ G ) η ϕ p ( t ) + ( G 0 + Δ G ) η ϕ p ( t ) + k q = 1 N L p q Λ 2 φ q ( t ) c q = 1 N d p q Λ 1 ϕ q ( t ) k q = 1 N L p q Λ 2 ϕ q ( t ) + k q = 1 N L p q Λ 2 ϕ q ( t ) + ω p ( t ) } .
Applying the controller ω p ( t ) into (23), we have the following formula
t 0 R D t ξ V p ( t ) = ( t t 0 ) ξ Γ ( 1 ξ ) ( υ p T ( t 0 ) P υ p ( t 0 ) ) + 2 υ p T ( t ) P { ( E 0 + Δ E ) φ p ( t ) A ϕ p ( t ) ( E 0 + Δ E ) ϕ p ( t ) + ( E 0 + Δ E ) ϕ p ( t ) + ( G 0 + Δ G ) η φ p ( t ) B γ ϕ p ( t ) ( G 0 + Δ G ) η ϕ p ( t ) + ( G 0 + Δ G ) η ϕ p ( t ) + k q = 1 N L p q Λ 2 φ q ( t ) c q = 1 N d p q Λ 1 ϕ q ( t ) + B γ φ p ( t ) k q = 1 N L p q Λ 2 ϕ q ( t ) + k q = 1 N L p q Λ 2 ϕ q ( t ) + ( A E ) ϕ p ( t ) G η φ p ( t ) + c q = 1 N d p q Λ 1 φ q ( t ) k q = 1 N L p q Λ 2 ϕ q ( t ) } .
The above formula can be simplified as
t 0 R D t ξ V p ( t ) = ( t t 0 ) ξ Γ ( 1 ξ ) ( υ p T ( t 0 ) P υ p ( t 0 ) ) + 2 υ p T ( t ) P { ( E 0 + Δ E ) υ p ( t ) + ( G 0 + Δ G ) η υ p ( t ) + k q = 1 N L p q Λ 2 υ q ( t ) + B γ υ p ( t ) + c q = 1 N d p q Λ 1 υ q ( t ) } .
Utilizing the Lemma 2 into (25), we have
2 υ p T ( t ) P ( E 0 + Δ E ) υ p ( t ) 2 υ p T ( t ) P E 0 υ p ( t ) + 2 υ p T ( t ) P Δ E υ p ( t ) υ p T ( t ) ( P E 0 + E 0 T P ) υ p ( t ) + 2 υ p T ( t ) P H 1 Ξ 1 T 1 υ p ( t ) υ p T ( t ) ( P E 0 + E 0 T P ) υ p ( t ) + υ p T ( t ) ϵ 1 1 P H 1 H 1 T X υ p ( t ) + υ p T ( t ) ϵ 1 T 1 T 1 T υ p ( t ) υ p T ( t ) { P E 0 + E 0 T ( t ) P + ϵ 1 1 P H 1 H 1 T P + ϵ 1 T 1 T 1 T } υ p ( t ) ,
and
2 υ p T ( t ) P ( G 0 + Δ G ) η υ p ( t ) υ p T ( t ) { ϵ 2 1 P H 2 H 2 T P + ϵ 2 T 2 T 2 T + ϵ 3 1 P G G T P + ϵ 3 θ 1 } 2 υ p T ( t ) P B γ υ p ( t ) υ p T ( t ) { ϵ 4 1 P B B T P + ϵ 4 θ 2 } .
By substituting (26) and (27) into (25), we get
t 0 R D t ξ V p ( t ) ( t t 0 ) ξ Γ ( 1 ξ ) ( υ p T ( t 0 ) P υ p ( t 0 ) ) + 2 υ p T ( t ) { P E 0 + E 0 T ( t ) P + ϵ 1 1 P H 1 H 1 T P + ϵ 1 T 1 T 1 T + ϵ 2 1 P H 2 H 2 T P + ϵ 2 T 2 T 2 T + ϵ 3 1 P G G T P + ϵ 3 θ 1 + ϵ 4 1 P B B T P + ϵ 4 θ 2 + k q = 1 N L p q Λ 2 + c q = 1 N d p q Λ 1 } υ p ( t ) .
The above formula can be further simplified as
t 0 R D t ξ V p ( t ) ( t t 0 ) ξ Γ ( 1 ξ ) ( υ p T ( t 0 ) P υ p ( t 0 ) ) λ m i n ( Ω ) υ p ( t ) T υ p ( t ) ( t t 0 ) ξ Γ ( 1 ξ ) ( υ p T ( t 0 ) P υ p ( t 0 ) ) λ m i n ( Ω ) λ m a x ( P ) υ p ( t ) T P υ p ( t ) ( t t 0 ) ξ Γ ( 1 ξ ) ( υ p T ( t 0 ) P υ p ( t 0 ) ) χ V p ( t ) ,
where χ = λ m i n ( Ω ) λ m a x ( P ) , utilizing Lemma 5, and if Ω < 0 is satisfied, we infer that the error system (17) is asymptotically stable. Then, we find that the FCDNDP is asymptotically synchronized under the original feedback controller. □
In the light of the above Theorem 1, we can obtain the following corollary.
Corollary 1.
If the fractional-order complex network is composed of N nodes, which can be depicted as
t 0 C D t ξ ϕ p ( t ) = A ϕ p ( t ) + Q ϕ p ( t ) + B γ ϕ p ( t ) + c q = 1 N d p q Λ 1 ϕ q ( t ) ,
and
t 0 C D t ξ φ p ( t ) = E φ p ( t ) + Y φ p ( t ) + G η φ p ( t ) + k q = 1 N L p q Λ 2 φ q ( t ) + ω p ( t ) ,
which is similar to system (9),
Y = ( Y 0 + Δ Y ) , Δ Y = H 3 Ξ 3 K 3 ,
where Q and Y 0 are the real constant matrices, and other parameters are the same as system (9) and system (10). We can obtain the fractional-order complex network which is composed of the system (30) and system (31). They are asymptotically synchronized. The proof for Corollary 1 is similar to Theorem 1.

4. Numerical Simulation

In this part, two numerical instances are presented to illustrate the effectiveness of the proposed method.
Example 1.
Suppose that the following FCDNDP is composed of N nodes, and every node is a n-dimensional system which is given as follows
t 0 C D t ξ ϕ p ( t ) = A ϕ p ( t ) + B γ ϕ p ( t ) + c q = 1 N d p q Λ 1 ϕ q ( t ) ,
and
t 0 C D t ξ φ p ( t ) = ( E 0 + Δ E ) φ p ( t ) + ( G 0 + Δ G ) η φ p ( t ) + k q = 1 N L p q Λ 2 φ q ( t ) + ω p ( t ) ,
where p = 1 , 2 , 3 , 4 . c = 0.01 ; k = 0.02 ; ς ϕ p ( t ) = ( 0.25 tanh ( x 1 ( t ) ) , 0.25 tanh ( x 2 ( t ) ) , 0.25 tanh ( x 3 ( t ) ) , 0.25 tanh ( x 4 ( t ) ) ) ; η φ p ( t ) = ( 8.10 tanh ( x 1 ( t ) ) , 8.10 tanh ( x 2 ( t ) ) , 8.10 tanh ( x 3 ( t ) ) , 8.10 tanh ( x 4 ( t ) ) ) .
The parameter matrices and weight connection matrices can be expressed as
A = 20.15 0 0 0 0 5.98 0 0 0 0 12.47 0 0 0 0 16.56 , B = 0.2 0 0.5 0.78 0 1.15 1 0.35 0.3 1 1 0 0 0 0 2.15 ,
E = 2.14 0 0 0 0 0.16 0 0 0 0 5.57 0 0 0 0 0.23 , D = 1 0.025 0.1 0.6 2.1 1 0 1 0 1 2.25 1 1 0 1 1.35 .
Meanwhile, the inner coupling matrices can be expressed as
Λ 1 = 0.25 0 0 0 0 0.25 0 0 0 0 0.25 0 0 0 0 0.25 , Λ 2 = 0.35 0 0 0 0 0.35 0 0 0 0 0.35 0 0 0 0 0.35 .
The matrices of parameter uncertainties can be expressed as
H 1 = 0.1 0 0 0 0 0.5 0 0 0 0 0.2 0 0 0 0 0.3 , K 1 = 0.2 0 0 0 0 0.1 0 0 0 0 0.2 0 0 0 0 0.5 ,
Ξ 1 = cos ( x 1 ( t ) ) 0 0 0 0 cos ( x 2 ( t ) ) 0 0 0 0 cos ( x 3 ( t ) ) 0 0 0 0 cos ( x 4 ( t ) ) .
The external coupling matrices can be expressed as
d p q = 1 0 1.5 0 2.5 2.1 1 1 0.2 2.25 2 1.5 0 0.2 0.1 0.5 0 1.5 0.1 3 0 1 3 0 1 1.5 0 1.25 1.15 1 2.5 0 1.5 2 1.5 0 1.5 0 2 0.5 0 1 0 1 2.25 1.45 1.5 3 0.5 1.1 1.01 0 2 2 1 1.5 0 1 1 0 0.5 0 1 2.1 0 1.5 1.5 1.5 1.01 1.45 2 0.1 1.01 1 1.5 0 1 1.45 2.25 0.2 1 2.5 0 1.5 1 0.2 0.5 1 1.15 1.25 1.5 3 0.5 0.5 0 1.15 3 0 0.2 1 ,
L p q = 2 0 0 1.5 1.5 3 1 0 1 0.5 1 1 0 1 3.5 0 0.25 2 1.5 0 0 1 3.5 0.5 0.25 0 0 1.5 2.25 0.5 1 1 1.25 0.15 0 1 2 2.5 0 1.5 0.15 2.25 1 0 1 2 0 0.15 0.25 3 0 0.26 3 0 0 1 2 1 1 1 1.25 0 0.25 1.25 2 0 2.5 0.15 0 0.25 0 1.5 1 0 1 2.25 1.5 1 0.15 0 2.25 1.25 0 1.5 0.5 1 1 0 0.25 2 0.5 0 1 2 0 0.5 3 1 0 1.5 .
The simulation results are shown in Figure 1 and Figure 2, which show the time waveforms of the system errors e i 1 . Figure 1 shows the variation of the synchronization error system without a controller. Figure 2 shows the variation of the synchronization error system with a controller. Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 show the variation in the synchronization error systems e i 2 , e i 3 , e i 4 , e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) . According to the simulation results and figures, one can obtain that the error system is driven to initial point, i.e., the FCDNDP is synchronized, which verifies the validity of the proposed controller.
Example 2.
Suppose that the following FCDNDP is composed of N nodes and every node is a n-dimensional system, which is given as follows
t 0 C D t ξ ϕ p ( t ) = A ϕ p ( t ) + B γ ϕ p ( t ) + c q = 1 N d p q Λ 1 ϕ q ( t ) ,
and
t 0 C D t ξ φ p ( t ) = ( E 0 + Δ E ) φ p ( t ) + ( G 0 + Δ G ) η φ p ( t ) + k q = 1 N L p q Λ 2 φ q ( t ) + ω p ( t ) .
The parameter matrices and weight connection matrices can be expressed as
A = 30.15 0 0 0 0 7.99 0 0 0 0 15.73 0 0 0 0 24.89 , B = 1 0 0.5 0.35 0 2.15 1 0.35 0 1 1 0 0 0 0 0.15 ,
E = 3.15 0 0 0 0 0.52 0 0 0 0 8.51 0 0 0 0 0.45 , D = 1.048 0.025 0.1 0.6 2.1 0.85 0 1.87 0 1 3.25 1 1.86 0 1 1.35 .
Meanwhile, the inner coupling matrices can be expressed as
Λ 1 = 0.25 0 0 0 0 0.25 0 0 0 0 0.25 0 0 0 0 0.25 , Λ 2 = 0.35 0 0 0 0 0.35 0 0 0 0 0.35 0 0 0 0 0.35 .
The matrices of parameter uncertainties can be expressed as
H 2 = 1 0 0 0 0 0.9 0 0 0 0 1 0 0 0 0 0.3 , K 1 = 1 0 0 0 0 0.9 0 0 0 0 0.8 0 0 0 0 1 ,
Ξ 2 = 0.41 cos ( x 1 ( t ) ) 0 0 0 0 cos ( x 2 ( t ) ) 0 0 0 0 0.3 cos ( x 3 ( t ) ) 0 0 0 0 0.13 cos ( x 4 ( t ) ) .
The external coupling matrices can be expressed as
d p q = 1 0 1.5 0 2.5 2.1 1 1 0.2 2.25 2 1.5 0 0.2 0.1 0.5 0 1.5 0.1 3 0 1 3 0 1 1.5 0 1.25 1.15 1 2.5 0 1.5 2 1.5 0 1.5 0 2 0.5 0 1 0 1 2.25 1.45 1.5 3 0.5 1.1 1.01 0 2 2 1 1.5 0 1 1 0 0.5 0 1 2.1 0 1.5 1.5 1.5 1.01 1.45 2 0.1 1.01 1 1.5 0 1 1.45 2.25 0.2 1 2.5 0 1.5 1 0.2 0.5 1 1.15 1.25 1.5 3 0.5 0.5 0 1.15 3 0 0.2 1 ,
L p q = 2 0 0 1.5 1.5 3 1 0 1 0.5 1 1 0 1 3.5 0 0.25 2 1.5 0 0 1 3.5 0.5 0.25 0 0 1.5 2.25 0.5 1 1 1.25 0.15 0 1 2 2.5 0 1.5 0.15 2.25 1 0 1 2 0 0.15 0.25 3 0 0.26 3 0 0 1 2 1 1 1 1.25 0 0.25 1.25 2 0 2.5 0.15 0 0.25 0 1.5 1 0 1 2.25 1.5 1 0.15 0 2.25 1.25 0 1.5 0.5 1 1 0 0.25 2 0.5 0 1 2 0 0.5 3 1 0 1.5 .
Meanwhile, the other parameters are the same as in Example 1.
The simulation results are shown in Figure 11 and Figure 12, which show the time waveforms of the system errors e i 1 . Figure 11 shows the variation in the synchronization error system without a controller. Figure 12 shows the variation in the synchronization error system with a controller. Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19 and Figure 20 shows the variation in synchronization error systems e i 2 , e i 3 , e i 4 , e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) . According to the simulation results and figures, one can obtain that the error system is driven to the initial point, i.e., the FCDNDP is synchronized, which verifies the validity of the proposed controller.

5. Conclusions

In this paper, by applying the fractional differential equation theory and differential inclusion theory, a driving response system with different structures is established. The driving system is a general driving model, and the response system is a response model with parameter uncertainties. Based on the Lyapunov direct method, a sufficient condition for the asymptotic synchronization of FCDNDP is established. When we deal with a similar problem, this method can ensure that the obtained results have better application and less conservatism. Finally, the validity and feasibility of the theoretical outcomes are confirmed by two simulation instances. In the future, we hope to extend the technique proposed in this paper to the asymptotic synchronization of FCDNDP with time delay.

Author Contributions

X.H. proposed the main the idea and prepared the manuscript initially. T.L. gave the numerical simulation of this paper. D.L. revised the English grammar of this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work is partly supported by the Project of the Science and Technology Department in Sichuan Province (Grant No.2021ZYD0004), Fund of Sichuan University of Science and Engineering (Grant No. 2020RC26, 2020RC42), The Postgraduate Innovation Fund Project of Sichuan University of Science and Engineering (Grant No. y2022189).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no competing interest.

References

  1. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  2. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Hackensack, NJ, USA, 2000. [Google Scholar]
  3. Bagley, R.L. On the Fractional Calculus Model of Viscoelastic Behavior. J. Rheol. 1998, 30, 133–155. [Google Scholar] [CrossRef]
  4. Ge, Z.M.; Jhuang, W.R. Chaos, control and synchronization of a fractional order rotational mechanical system with a centrifugal governor. Chaos Solitons Fractals 2007, 33, 270–289. [Google Scholar] [CrossRef]
  5. Huang, L.L.; He, S.J. Stability of fractional state space system and its application to fractional order chaotic system. Acta Phys. Sin. Chin. Ed. 2011, 60, 119419573. [Google Scholar] [CrossRef]
  6. Hudson, J.L.; Kube, M.; Adomaitis, R.A.; Kevrekidis, I.G.; Lapedes, A.S.; Farber, R.M. Nonlinear signal processing and system identification: Applications to time series from electrochemical reactions. Chem. Eng. Sci. 1990, 45, 2075–2081. [Google Scholar] [CrossRef]
  7. Argenti, F.; De Angeli, A.; Del Re, E.; Genesio, R.; Pagni, P.; Tesi, A. Secure communications based on discrete time chaotic systems. Kybernetika 1997, 1, 41–50. [Google Scholar]
  8. Gerhards, M.; Schlerf, M.; Rascher, U.; Udelhoven, T.; Juszczak, R.; Alberti, G.; Miglietta, F.; Inoue, Y. Analysis of Airborne Optical and Thermal Imagery for Detection of Water Stress Symptoms. Remote Sens. 2018, 10, 1139. [Google Scholar] [CrossRef]
  9. Lei, Y.; Zheng, W. Research on robot automation and control problems. World Inverters 2015, 3, 86–88. [Google Scholar]
  10. Boccaletti, S.; Latora, V.; Moreno, Y.; Chavez, M.; Hwang, D.U. Complex networks: Structure and dynamics. Phys. Rep. 2006, 424, 175–308. [Google Scholar] [CrossRef]
  11. Wang, Z.; Shu, H.; Liu, Y.; Ho, D.W.; Liu, X. Robust stability analysis of generalized neural networks with discrete and distributed time delays. Chaos Solitons Fractals 2006, 30, 886–896. [Google Scholar] [CrossRef]
  12. Barabâsi, A.L.; Jeong, H.; Néda, Z.; Ravasz, E.; Schubert, A.; Vicsek, T. Evolution of the social network of scientific collaborations. Phys. Astatistical Mech. Its Appl. 2002, 311, 590–614. [Google Scholar] [CrossRef]
  13. Melián, C.J.; Bascompte, J.; Jordano, P.; Krivan, V. Diversity in a complex ecological network with two interaction types. Oikos 2010, 118, 122–130. [Google Scholar] [CrossRef]
  14. Zio, E.; Golea, L.R.; Mrs, C. Identifying groups of critical edges in a realistic electrical network by multi-objective genetic algorithms. Reliab. Eng. Syst. Saf. 2012, 99, 172–177. [Google Scholar] [CrossRef]
  15. Ma, T.; Zhang, J. Hybrid synchronization of coupled fractional-order complex networks. Neurocomputing 2015, 157, 166–172. [Google Scholar] [CrossRef]
  16. Li, H.L.; Cao, J.; Hu, C.; Zhang, L.; Wang, Z. Global synchronization between two fractional-order complex networks with non-delayed and delayed coupling via hybrid impulsive control. Neurocomputing 2019, 356, 31–39. [Google Scholar] [CrossRef]
  17. Yang, L.X.; Jiang, J. Adaptive synchronization of drive-response fractional-order complex dynamical networks with uncertain parameters. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 1496–1506. [Google Scholar] [CrossRef]
  18. Zhu, D.; Wang, R.; Liu, C.; Duan, J. Projective synchronization via adaptive pinning control for fractional-order complex network with time-varying coupling strength. Int. J. Mod. Phys. C 2019, 30, 268–968. [Google Scholar] [CrossRef]
  19. Li, H.L.; Cao, J.; Jiang, H.; Alsaedi, A. Graph theory based finite-time synchronization of fractional-order complex dynamical networks. J. Frankl. Inst. 2018, 355, 5771–5789. [Google Scholar] [CrossRef]
  20. Wu, H.; Wang, L.; Niu, P.; Wang, Y. Global projective synchronization in finite time of nonidentical fractional-order neural networks based on sliding mode control strategy. Neurocomputing 2017, 235, 264–273. [Google Scholar] [CrossRef]
  21. Lin, D.U.; Yong, Y.; Lei, Y. Synchronization in a fractional-order dynamic network with uncertain parameters using an adaptive control strategy. Appl. Math. Mech. 2018, 39, 353–364. [Google Scholar]
  22. Shen, Y.; Li, C. LMI-based finite-time boundedness analysis of neural networks with parametric uncertainties. Neurocomputing 2008, 71, 502–507. [Google Scholar] [CrossRef]
  23. Wong, W.K.; Li, H.; Leung, S. Robust synchronization of fractional-order complex dynamical networks with parametric uncertainties. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 4877–4890. [Google Scholar] [CrossRef]
  24. Li, H.J. State Estimation for Fractional-Order Complex Dynamical Networks with Linear Fractional Parametric Uncertainty. Abstr. Appl. Anal. 2013, 2013, 151–164. [Google Scholar] [CrossRef]
  25. Samli, R.; Yucel, E. Global robust stability analysis of uncertain neural networks with time varying delays. Neurocomputing 2015, 167, 371–377. [Google Scholar] [CrossRef]
  26. Ding, Z.; Shen, Y. Projective synchronization of nonidentical fractional-order neural networks based on sliding mode controller. Neural Netw. 2016, 76, 97–105. [Google Scholar] [CrossRef]
  27. Hu, T.; Zhang, X.; Zhong, S. Global asymptotic synchronization of nonidentical fractional-order neural networks. Neurocomputing 2018, 313, 39–46. [Google Scholar] [CrossRef]
  28. Suntonsinsoungvon, E.; Udpin, S. Exponential stability of discrete-time uncertain neural networks with multiple time-varying leakage delays. Math. Comput. Simul. 2020, 171, 233–245. [Google Scholar] [CrossRef]
  29. Aadhithiyan, S.; Raja, R.; Kou, B.; Selvam, G.; Niezabitowski, M.; Lim, C.P.; Cao, J. Asymptotic synchronization of fractional order non-identical complex dynamical networks with Parameter Uncertainties. Math. Methods Appl. Sci. 2022, 2022, 1–19. [Google Scholar] [CrossRef]
  30. Wei, J.; Pan, G. Design of a Sliding Mode Controller for Synchronization of Fractional-Order Chaotic Systems with Different Structures. J. Shanghai Jiaotong Univ. 2016, 50, 849–853, 860. [Google Scholar]
  31. Zeng, H.B.; He, Y.; Wu, M.; Xiao, H.Q. Improved Conditions for Passivity of Neural Networks With a Time-Varying Delay. IEEE Trans. Cybern. 2014, 44, 785–792. [Google Scholar] [CrossRef]
Figure 1. Synchronization errors e i 1 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) without the controllers.
Figure 1. Synchronization errors e i 1 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) without the controllers.
Fractalfract 06 00441 g001
Figure 2. Synchronization errors e i 1 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) with the controllers.
Figure 2. Synchronization errors e i 1 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) with the controllers.
Fractalfract 06 00441 g002
Figure 3. Synchronization errors e i 2 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) without the controllers.
Figure 3. Synchronization errors e i 2 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) without the controllers.
Fractalfract 06 00441 g003
Figure 4. Synchronization errors e i 2 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) with the controllers.
Figure 4. Synchronization errors e i 2 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) with the controllers.
Fractalfract 06 00441 g004
Figure 5. Synchronization errors e i 3 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) without the controllers.
Figure 5. Synchronization errors e i 3 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) without the controllers.
Fractalfract 06 00441 g005
Figure 6. Synchronization errors e i 3 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) with the controllers.
Figure 6. Synchronization errors e i 3 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) with the controllers.
Fractalfract 06 00441 g006
Figure 7. Synchronization errors e i 4 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) without the controllers.
Figure 7. Synchronization errors e i 4 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) without the controllers.
Fractalfract 06 00441 g007
Figure 8. Synchronization errors e i 4 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) with the controllers.
Figure 8. Synchronization errors e i 4 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (32) and (33) with the controllers.
Fractalfract 06 00441 g008
Figure 9. Synchronization errors e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) of the fractional-order system (32) and (33) without the controllers.
Figure 9. Synchronization errors e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) of the fractional-order system (32) and (33) without the controllers.
Fractalfract 06 00441 g009
Figure 10. Synchronization errors e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) of the fractional-order system (32) and (33) with the controllers.
Figure 10. Synchronization errors e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) of the fractional-order system (32) and (33) with the controllers.
Fractalfract 06 00441 g010
Figure 11. Synchronization errors e i 1 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) without the controllers.
Figure 11. Synchronization errors e i 1 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) without the controllers.
Fractalfract 06 00441 g011
Figure 12. Synchronization errors e i 1 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) with the controllers.
Figure 12. Synchronization errors e i 1 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) with the controllers.
Fractalfract 06 00441 g012
Figure 13. Synchronization errors e i 2 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) without the controllers.
Figure 13. Synchronization errors e i 2 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) without the controllers.
Fractalfract 06 00441 g013
Figure 14. Synchronization errors e i 2 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) with the controllers.
Figure 14. Synchronization errors e i 2 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) with the controllers.
Fractalfract 06 00441 g014
Figure 15. Synchronization errors e i 3 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) without the controllers.
Figure 15. Synchronization errors e i 3 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) without the controllers.
Fractalfract 06 00441 g015
Figure 16. Synchronization errors e i 3 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) with the controllers.
Figure 16. Synchronization errors e i 3 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) with the controllers.
Fractalfract 06 00441 g016
Figure 17. Synchronization errors e i 4 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) without the controllers.
Figure 17. Synchronization errors e i 4 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) without the controllers.
Fractalfract 06 00441 g017
Figure 18. Synchronization errors e i 4 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) with the controllers.
Figure 18. Synchronization errors e i 4 ( i = 1 , 2 , · · · , 10 ) of the fractional-order system (34) and (35) with the controllers.
Fractalfract 06 00441 g018
Figure 19. Synchronization errors e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) of the fractional-order system (34) and (35) without the controllers.
Figure 19. Synchronization errors e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) of the fractional-order system (34) and (35) without the controllers.
Fractalfract 06 00441 g019
Figure 20. Synchronization errors e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) of the fractional-order system (34) and (35) with the controllers.
Figure 20. Synchronization errors e i j ( i = 1 , 2 , · · · , 10 ; j = 1 , 2 , 3 , 4 ) of the fractional-order system (34) and (35) with the controllers.
Fractalfract 06 00441 g020
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

He, X.; Li, T.; Liu, D. Asymptotic Synchronization of Fractional-Order Complex Dynamical Networks with Different Structures and Parameter Uncertainties. Fractal Fract. 2022, 6, 441. https://doi.org/10.3390/fractalfract6080441

AMA Style

He X, Li T, Liu D. Asymptotic Synchronization of Fractional-Order Complex Dynamical Networks with Different Structures and Parameter Uncertainties. Fractal and Fractional. 2022; 6(8):441. https://doi.org/10.3390/fractalfract6080441

Chicago/Turabian Style

He, Xiliang, Tianzeng Li, and Dehui Liu. 2022. "Asymptotic Synchronization of Fractional-Order Complex Dynamical Networks with Different Structures and Parameter Uncertainties" Fractal and Fractional 6, no. 8: 441. https://doi.org/10.3390/fractalfract6080441

APA Style

He, X., Li, T., & Liu, D. (2022). Asymptotic Synchronization of Fractional-Order Complex Dynamical Networks with Different Structures and Parameter Uncertainties. Fractal and Fractional, 6(8), 441. https://doi.org/10.3390/fractalfract6080441

Article Metrics

Back to TopTop