Next Article in Journal
Fractional Optimal Control Model and Bifurcation Analysis of Human Syncytial Respiratory Virus Transmission Dynamics
Next Article in Special Issue
Fractional-Order Load Frequency Control of an Interconnected Power System with a Hydrogen Energy-Storage Unit
Previous Article in Journal
Pore Structure Quantification and Fractal Characterization of MSA Mortar Based on 1H Low-Field NMR
Previous Article in Special Issue
Optimal Design of Fractional-Order PID Controllers for a Nonlinear AWS Wave Energy Converter Using Hybrid Jellyfish Search and Particle Swarm Optimization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Controllability of Fractional Complex Networks

School of Mathematics and Computer Science, Northwest Minzu University, Lanzhou 730030, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2024, 8(1), 43; https://doi.org/10.3390/fractalfract8010043
Submission received: 5 December 2023 / Revised: 5 January 2024 / Accepted: 8 January 2024 / Published: 11 January 2024
(This article belongs to the Special Issue Fractional Order Controllers for Non-linear Systems)

Abstract

:
Controllability is a fundamental issue in the field of fractional complex network control, yet it has not received adequate attention in the past. This paper is dedicated to exploring the controllability of complex networks involving the Caputo fractional derivative. By utilizing the Cayley–Hamilton theorem and Laplace transformation, a concise proof is given to determine the controllability of linear fractional complex networks. Subsequently, leveraging the Schauder Fixed-Point theorem, controllability Gramian matrix, and fractional calculus theory, we derive controllability conditions for nonlinear fractional complex networks with a weighted adjacency matrix and Laplacian matrix, respectively. Finally, a numerical method for the controllability of fractional complex networks is obtained using Matlab (2021a)/Simulink (2021a). Three examples are provided to illustrate the theoretical results.

1. Introduction

In recent decades, due to the widespread existence of complex networks in nature and society, complex networks have been undergoing a period of rapid development in various interdisciplinary research fields, particularly in engineering [1], mathematics [2], physics [3], biology [4] and other related areas. The ultimate goal of researching complex networks is to control them to facilitate our lives. Despite significant efforts being dedicated to comprehending the interactions between complex networks and their dynamic behaviors, controlling these networks remains a prominent challenge [5].
In order to establish a framework for controlling complex networks, a fundamental and essential step is to study their controllability. Controllability ensures that by manipulating a subset of nodes with appropriate control inputs, every node within the network can be steered to the desired state [6,7]. Many scholars from the fields of control, physics, and mathematics [8,9,10,11] have been inspired by the pioneering paper on the controllability of complex networks [12]. Liu et al. [8] analyzed the controllability of discrete-time dynamic networks with both switching and fixed interaction topologies. Cai [9] utilized the condition number of a matrix as a metric to quantitatively assess controllability. Chen [10] incorporated pinning control strategies into the study of controllability in directed networks. Meanwhile, Whalen et al. [11] developed a group representational framework to tackle the controllability of nonlinear networks that exhibit explicit symmetries.
Due to its memory and hereditary properties, fractional calculus has garnered considerable attention across various scientific fields and is widely adopted in science and engineering [13,14,15,16,17,18,19,20]. Fractional operators have been integrated into traditional complex networks, significantly improving model accuracy. It is worth noting that while the controllability of integer-order complex networks is relatively mature, the controllability results of fractional networks are still in their infancy. Some papers have addressed the controllability of fractional systems [21,22,23,24,25], but only a few papers have delved into the controllability of fractional complex networks because of their inherent complexity and long memory characteristics. Zhang et al. [26] examined the controllability of linear fractional directed complex networks, and subsequently, they [27] extended this work to explore the controllability of linear fractional dynamical networks with specific topological structures. However, given the ubiquitous nature of nonlinearity, it is imperative to consider the control of fractional complex networks with nonlinear dynamics. As far as we know, no research has yet been explored on the controllability of nonlinear fractional complex networks.
Motivated by the above analysis, this manuscript aims to address the controllability issues for fractional complex networks through rigorous mathematical theory. The distinctive contributions are summarized as follows: (i) A concise technique is established to prove the controllability of fractional complex networks with linear dynamic behavior. This allows the controllability of similar networks to be easily determined. (ii) For the first time, nonlinear fractional complex networks with weighted adjacency matrices or Laplacian matrices are shown to be controllable, effectively making up for the deficiency left by existing methods. (iii) Matlab/Simulink is used to obtain a numerical implementation for the controllability of fractional complex networks. Three examples are presented to validate theoretical results.
The remainder of this paper is outlined as follows: Section 2 recalls necessary definitions and Lemmas, Section 3 provides a more concise proof, Section 4 contains the main results, Section 5 presents three corresponding examples, and finally, a brief conclusion is given.

2. Preliminaries

This section presents algebraic graph theory and some necessary theories since we will use them later.
A graph is defined by G = ( V , E ) , where V = { v 1 , v 2 , , v N } and E = { ( v j , v i ) , v j , v i V } V × V denote node sets and link sets, respectively. Additionally, ( v j , v i ) E represents a link from node v j to node v i . The concepts of a directed graph, undirected graph, and adjacency matrix can be found in [28,29,30].
Definition 1
([31]). The Caputo fractional derivative of order p of a function f is defined by
C D 0 , t p f ( t ) = 1 Γ ( n p ) 0 t ( t ς ) n p 1 f ( n ) ( ς ) d ς ,
where n 1 < p n N and Γ ( · ) is the Gamma function.
Definition 2
([6]). The two-parameter Mittag–Leffler matrix function for a matrix A R N × N is defined as
E p , q ( A ) = k = 0 A k Γ ( p k + q ) , p , q > 0 .
Note that when q = 1 , E p , 1 ( A ) = E p ( A ) becomes a one-parameter Mittag–Leffler matrix function. Additionally, the Laplace transform of t q 1 E p , q ( ± A t p ) is given by
L t q 1 E p , q ( ± A t p ) ; s = s p q s p I A .
Lemma 1
([32]). For n 1 < p n N , the Laplace transform of Caputo fractional derivative is given by
L C D 0 , t p f ( t ) ; s = s p F ( s ) k = 0 n 1 s p k 1 f ( k ) ( 0 ) .
In particular, if 0 < p 1 , then
L C D 0 , t p f ( t ) ; s = s p F ( s ) s p 1 f ( 0 ) ,
where F ( s ) = L f ( t ) ; s .
Lemma 2
(Schauder Fixed-Point theorem [33]). Let M be a non-empty, closed, bounded, convex subset of a Banach space X. Suppose G : M M is a compact operator, then G has a fixed point.

3. Controllability Analysis of Linear Fractional Complex Networks

In this part, a fractional complex network with linear dynamic behavior is described as follows:
C D 0 , t p x i ( t ) = c i x i ( t ) + j = 1 N a i j x j ( t ) + k = 1 m b i k u k ( t ) , i = 1 , 2 , , N ,
in which 0 < p 1 , x i ( t ) R represents the state of the node i, c i x i ( t ) denotes the intrinsic dynamics, and a i j 0 denotes the weight of the network. Here, m represents the number of controllers and u k ( t ) R is the outer controller. For example, the u k ( t ) can stand for the command from a leader in a social network or the signal from equipment in a sensor network. b i k represents the signal strength of the outer controller.
Let x ( t ) = [ x 1 ( t ) , x 2 ( t ) , , x N ( t ) ] T R N be the whole state of the network (4), and u ( t ) = [ u 1 ( t ) , u 2 ( t ) , , u m ( t ) ] T R m be the total outer controller. The linear network (4) can be reformulated as
C D 0 , t p x ( t ) = ( A + C ) x ( t ) + B u ( t ) ,
where A represents the weighted (unweighted) adjacency matrix, C = d i a g { c 1 , c 2 , , c N } , and B R N × m stands for the control matrix.
Definition 3
([7]). For a finite time t f , any initial state x ( 0 ) = x 0 R N and any final state x ( t f ) = x t f R N , if there exists a controller u ( t ) , t I = [ 0 , t f ] that satisfies x ( t f , x 0 , u ( t ) ) = x t f , then the linear fractional complex network (4) is called controllable on I.
Lemma 3
([34]). The linear fractional complex network (4) is controllable if the controllability Gramian matrix
W F = 0 t f E p , p ( ( A + C ) ( t f ς ) p ) B B * E p , p ( ( A + C ) * ( t f ς ) p ) d ς
is invertible.
For the linear fractional complex network (4), a controllability result is proposed in [26], but we give a more concise proof below.
Theorem 1.
The fractional complex network (4) is controllable if N × ( N m ) controllability matrix
Q F = [ B , ( A + C ) B , , ( A + C ) N 1 B ]
is of full rank, which is r a n k ( Q F ) = N .
Proof. 
Using the Cayley–Hamilton theorem [35], one gets
( A + C ) s = k = 0 N 1 a k s ( A + C ) k , s = N , N + 1 , , .
It immediately follows from (8) that
t p 1 E p , p ( A + C ) t p = k = 0 N 1 t p k + p 1 Γ ( p k + p ) ( A + C ) k + k = N t p k + p 1 Γ ( p k + p ) ( A + C ) k = k = 0 N 1 φ k ( t ) ( A + C ) k ,
where φ k ( t ) = t p k + p 1 Γ ( p k + p ) + s = N t p s + p 1 Γ ( p s + p ) a k s , k = 0 , 1 , , N 1 .
By virtue of Laplace transformation and inverse Laplace transformation, the solution of (5) can be derived as follows:
x ( t ) = E p ( A + C ) t p x 0 + 0 t ( t ς ) p 1 E p , p ( ( A + C ) ( t ς ) p ) B u ( ς ) d ς .
When t = t f , from (9) and (10), one has
x ( t f ) = E p ( ( A + C ) t f p ) x 0 + 0 t f ( t f ς ) p 1 E p , p ( ( A + C ) ( t f ς ) p ) B u ( ς ) d ς = E p ( ( A + C ) t f p ) x 0 + k = 0 N 1 0 t f φ k ( t f ς ) ( A + C ) k B u ( ς ) d ς .
Let
Ψ = E p ( ( A + C ) t f p ) x 0 .
Then, applying (11) and (12) leads to
x ( t f ) Ψ = k = 0 N 1 ( A + C ) k B 0 t f φ k ( t f ς ) u ( ς ) d ς = B , ( A + C ) B , , ( A + C ) N 1 B ξ 0 ξ 1 ξ N 1 ,
in which ξ k = 0 t f φ k ( t f ς ) u ( ς ) d ς , k = 0 , 1 , , N 1 . Note that for any x 0 and x t f in R N , the adequacy and necessity condition with an external controller vector u ( t ) satisfying (13) is that
rank B , ( A + C ) B , , ( A + C ) N 1 B = N .
The theorem is, thus, proved. □
Remark 1.
Noticeably, the linear fractional complex network (4) denoted by a pair of matrices ( A + C , B ) is controllable if the rank of Q F equals N. In particular, Theorem 1 is also true for network (5) with matrix C = 0 , i.e., when it (5) degenerates into
C D 0 , t p x ( t ) = A x ( t ) + B u ( t ) ,
a pair of ( A , B ) is controllable as well.

4. Controllability Analysis of Nonlinear Fractional Complex Networks

Complex networks with different topology structures can describe different connection relationships in complex worlds. This section mainly focuses on nonlinear fractional complex networks with a weighted adjacency matrix and Laplacian matrix.

4.1. Nonlinear Networks Represented by a Weighted Adjacency Matrix

Consider a nonlinear fractional complex network with a weighted adjacency matrix depicted by
C D 0 , t p x i ( t ) = f i ( x ( t ) ) + j = 1 N a i j x j ( t ) + k = 1 m b i k u k ( t ) , i = 1 , 2 , , N ,
where 0 < p 1 , x i ( t ) R expresses the state of node i. a i j 0 denotes the weight of a link ( j , i ) , m represents the number of controllers, u k ( t ) R is the outer controller, and f i ( x ( t ) ) : R N R is the nonlinear continuous function.
The nonlinear fractional complex network (16) can be recast into the following vector form
C D 0 , t p x ( t ) = f ( x ( t ) ) + A x ( t ) + B u ( t ) .
Here, x ( t ) = [ x 1 ( t ) , x 2 ( t ) , , x N ( t ) ] T R N expresses the state vector; u ( t ) = [ u 1 ( t ) , u 2 ( t ) , , u m ( t ) ] T R m represents the total outer controller; A = ( a i j ) N × N stands for the weighted adjacency matrix; B = ( b i k ) N × m represents the control matrix; f ( x ( t ) ) = [ f 1 ( x ( t ) ) , f 2 ( x ( t ) ) , , f N ( x ( t ) ) ] T : R N R N is the nonlinear vector function.
Similar to Definition 3, the definition of controllability for network (16) can also be given.
Henceforth, let us make the following assumptions:
( A 1 ) There is a positive constant M that satisfies
f ( x ( t ) M , for t I .
( A 2 ) The vector function f ( x ( t ) ) that is nonlinear meets the Lipschitz condition in vector form as follows:
f ( x ( t ) ) f ( y ( t ) ) k x ( t ) y ( t ) ,
where k > 0 , t I , and x ( t ) , y ( t ) R N are any vectors.
( A 3 ) Let
a 1 = sup t , ς I { E p , p ( A ( t ς ) p ) } , a 2 = sup t , ς I { | ( t ς ) p 2 | E p , p 1 ( A ( t ς ) p ) } , a 3 = sup t I { | t p 1 | A E p , p ( A t p ) x 0 } , b 1 = B B * W F 1 , b 2 = max t I { E p ( A t p ) x 0 } , b 3 = y 1 + b 2 , c 1 = a 1 3 b 1 p 1 t f p + 1 k , c 2 = a 1 2 a 2 p ( 2 p ) b 1 t f 2 k , r = m a x { r 1 , r 2 , 1 } , r 1 = b 2 + a 1 p 1 t f p M + a 1 2 t f b 1 ( b 3 + a 1 p 1 t f p M ) , r 2 = t f 1 p Γ ( 2 p ) ( a 3 + a 2 t f M + a 1 a 2 2 p b 1 t f 2 p ( b 3 + a 1 p 1 t f p M ) ) .
Theorem 2.
Assume that conditions ( A 1 ) ( A 3 ) are satisfied and system (15) is controllable, then the nonlinear fractional complex network (16) with a weighted adjacency matrix is controllable on I.
Proof. 
Give a Banach space
X = x ( t ) : x ( t ) C ( I , R N ) and C D 0 , t p x ( t ) C ( I , R N ) ,
which is endowed with the norm
x ( t ) X = max { x ( t ) , C D 0 , t p x ( t ) } .
By utilizing (15), it is possible to construct a controller u ( t ) that can regulate any x ( t ) . In this case, the controller u ( t ) can be constructed as
u ( t ) = ( t f t ) 1 p B * E p , p ( A * ( t f t ) p ) W F 1 y 1 E p ( A t f p ) x 0 0 t f ( t f ς ) p 1 E p , p ( A ( t f ς ) p ) f ( x ( ς ) ) d ς ] .
Now, we will prove that the nonlinear operator G : X X has a fixed point, which is also a solution of (17). The G is designed as
( G x ) ( t ) = E p ( A t p ) x 0 + 0 t ( t ς ) p 1 E p , p ( A ( t ς ) p ) f ( x ( ς ) ) d ς + 0 t ( t ς ) p 1 E p , p ( A ( t ς ) p ) B u ( ς ) d ς .
Substituting (18) into (19) yields
( G x ) ( t ) = E p ( A t p ) x 0 + 0 t ( t ς ) p 1 E p , p ( A ( t ς ) p ) f ( x ( ς ) ) d ς + 0 t ( t ς ) p 1 E p , p ( A ( t ς ) p ) B B * ( t f ς ) 1 p × E p , p ( A * ( t f ς ) p ) W F 1 y 1 E p ( A t f p ) x 0 0 t f ( t f ω ) p 1 E p , p ( A ( t f ω ) p ) f ( x ( ω ) ) d ω ] d ς .
If the nonlinear operator G satisfies ( G x ) ( t ) = x ( t ) , then ( G x ) ( t f ) = x ( t f ) = y 1 in view of Lemma 3; that is, the nonlinear fractional complex network (16) can be steered from initial state x 0 to desired final state y 1 by the controller u ( t ) within a finite time t f . By assuming ( A 3 ) and using Equation (20), it is true that
( G x ) ( t ) r 1 .
Taking into account derivatives of ( G x ) ( t ) , one gets
( G x ) ( t ) r 2 t f p 1 Γ ( 2 p ) .
From (22), it follows that
C D 0 , t p ( G x ) ( t ) 1 Γ ( 1 p ) 0 t ( t ς ) p ( G x ) ( ς ) d ς = r 2 .
Define the subspace of X as B r = { x ( t ) X : x ( t ) X r } . It is easy to see that the set B r is closed, bounded and convex. Let x ( t ) B r , from (21) and (23), one has
( G x ) ( t ) X = max { ( G x ) ( t ) , C D 0 , t p ( G x ) ( t ) } r .
Thus, we obtain that the G maps B r into itself.
In what follows, we will prove that the G is continuous on B r . Let { x n ( t ) } be a sequence of functions in B r with
x n ( t ) x ( t ) 0 as n , t I .
It can be deduced directly that
( G x n ) ( t ) ( G x ) ( t ) 0 t ( t ς ) p 1 E p , p ( A ( t ς ) p ) ( f ( x n ( ς ) ) f ( x ( ς ) ) ) d ς + 0 t ( t ς ) p 1 E p , p ( A ( t ς ) p ) B B * ( t f ς ) 1 p E p , p ( A * ( t f ς ) p ) × W F 1 0 t f ( t f ω ) p 1 E p , p ( A ( t f ω ) p ) ( f ( x n ( ω ) ) f ( x ( ω ) ) ) d ω ] d ς a 1 p 1 t f p k max t I { x n ( t ) x ( t ) } + c 1 max t I { x n ( t ) x ( t ) } = ( a 1 p 1 t f p k + c 1 ) max t I { x n ( t ) x ( t ) }
and
( G x n ) ( t ) ( G x ) ( t ) 0 t ( t ς ) p 2 E p , p 1 ( A ( t ς ) p ) ( f ( x n ( ς ) ) f ( x ( ς ) ) d ς + 0 t ( t ς ) p 2 E p , p 1 ( A ( t ς ) p ) B B * ( t f ς ) 1 p E p , p ( A * ( t f ς ) p ) × W F 1 0 t f ( t f ω ) p 1 E p , p ( A ( t f ω ) p ) ( f ( x n ( ω ) ) f ( x ( ω ) ) ) d ω ] d ς a 2 t f k max t I { x n ( t ) x ( t ) } + c 2 max t I { x n ( t ) x ( t ) } = ( a 2 t f k + c 2 ) max t I { x n ( t ) x ( t ) } .
From (25), one obtains
C D 0 , t p ( G x n ) ( t ) C D 0 , t p ( G x ) ( t ) = 1 Γ ( 1 p ) 0 t ( t ς ) p ( ( G x n ) ( ς ) ( G x ) ( ς ) ) d ς t f 1 p Γ ( 2 p ) ( a 2 t f k + c 2 ) max t I { x n ( t ) x ( t ) } .
Combining (24) and (26) yields
( G x n ) ( t ) ( G x ) ( t ) X = max { ( G x n ) ( t ) ( G x ) ( t ) , C D 0 , t p ( G x n ) ( t ) C D 0 , t p ( G x ) ( t ) } 0 as n .
Thus, the nonlinear operator G is continuous on B r .
Furthermore, we will demonstrate that the set G ( B r ) is relatively compact. Given that x ( t ) B r , t , τ I , 0 < t < τ < t f , it is shown that
( G x ) ( t ) ( G x ) ( τ ) E p ( A t p ) E p ( A τ p ) x 0 + τ t ( t ς ) p 1 E p , p ( A ( t ς ) p ) f ( x ( ς ) ) d ς + 0 τ [ ( t ς ) p 1 E p , p ( A ( t ς ) p ) ( τ ς ) p 1 E p , p ( A ( τ ς ) p ) ] f ( x ( ς ) ) d ς + τ t ( t ς ) p 1 E p , p ( A ( t ς ) p ) B B * ( t f ς ) 1 p E p , p ( A * ( t f ς ) p ) × W F 1 y 1 E p ( A t f p ) x 0 0 t f ( t f ω ) p 1 E p , p ( A ( t f ω ) p ) f ( x ( ω ) ) d ω ] d ς + 0 τ [ ( t ς ) p 1 E p , p ( A ( t ς ) p ) ( τ ς ) p 1 E p , p ( A ( τ ς ) p ) ] B B * × ( t f ς ) 1 p E p , p ( A * ( t f ς ) p ) W F 1 y 1 E p ( A t f p ) x 0 0 t f ( t f ω ) p 1 E p , p ( A ( t f ω ) p ) f ( x ( ω ) ) d ω ] d ς
and
C D 0 , t p ( G x ) ( t ) C D 0 , t p ( G x ) ( τ ) 1 Γ ( 1 p ) τ t ( t ς ) p ( G x ) ( ς ) d ς + 0 τ [ ( t ς ) p ( τ ς ) p ] ( G x ) ( ς ) d ς ) .
Clearly, one has
( G x ) ( t ) ( G x ) ( τ ) X = max { ( G x ) ( t ) ( G x ) ( τ ) , C D 0 , t p ( G x ) ( t ) C D 0 , t p ( G x ) ( τ ) } 0 as t τ .
As a result, G ( B r ) is equicontinuous, and it is not difficult to obtain that G ( B r ) is uniformly bounded. By applying the well-known Arzelà–Ascoli theorem [33], the set G ( B r ) is relatively compact. Moreover, since the nonlinear operator G is continuous on the bounded set B r in X and maps it into the relatively compact set G ( B r ) in X, it can be conluded that G is compact. According to Lemma 2, the nonlinear operator G exists at a fixed point, which means the nonlinear fractional complex network (16) is controllable. This concludes the demonstration of Theorem 2. □
Remark 2.
The nonlinear fractional complex network (16) represented by a pair ( A , B , f ) is controllable if the system ( A , B ) is controllable and conditions ( A 1 ) ( A 3 ) hold.

4.2. Nonlinear Networks Represented by a Laplacian Matrix

A nonlinear fractional complex network with a Laplacian matrix can be described as follows:
C D 0 , t p x i ( t ) = f i ( x ( t ) ) + j = 1 N a i j ( x j ( t ) x i ( t ) ) + k = 1 m b i k u k ( t ) , i = 1 , 2 , , N ,
in which all parameters represent the same meaning as (16).
The nonlinear network (27) can be recast as follows:
C D 0 , t p x ( t ) = f ( x ( t ) ) L x ( t ) + B u ( t ) ,
in which L = D A is the Laplacian matrix of network (27), D = d i a g { d 1 , d 2 , , d N } , and d i = j = 1 , j i N a i j is the out-degree of node i.
Theorem 3.
If A is replaced by L , conditions ( A 1 ) ( A 3 ) are true and system (15) is controllable, then network (27) with the Laplacian matrix can be controlled over interval I.
Proof of Theorem 3.
Similar to Theorem 2, Theorem 3 can be proved. Therefore, the proof is omitted here. □

5. Numerical Implementation

In this part, a numerical algorithm is given to demonstrate the controllability of the networks mentioned above with the help of the FOTF Toolbox [36]. Based on the maximum matching method [12], the minimum of the driven nodes of fractional complex networks can be obtained. For convenience, we take fractional order p = 0.5 and the final time t f = 1 .
Before the simulation, we will state two facts that are used to build the Matlab/Simulink simulation model.
Firstly, the control function u in (18) can be written as u ( t ) = g ( t ) · h ( t ) , where
g ( t ) = ( t f t ) 1 p B * E p , p ( A * ( t f t ) p ) W F 1
and
h ( t ) = y 1 E p ( A t f p ) x 0 0 t f ( t f ς ) p 1 E p , p ( A ( t f ς ) p ) f ( x ( ς ) ) d ς .
The Simulink simulation model corresponding to control u is shown in Figure 1. Here, the Interpreted Matlab Function 1 and Interpreted Matlab Function 2 blocks are Matlab functions of g ( t ) and h ( t ) , respectively. The integral in h ( t ) is computed using the trapezoidal formula.
Secondly, the Caputo derivative Simulink block is constructed using the Riemann–Liouville block created by Professor Oustaloup with the combination of the classical integrator block. As we can see in Figure 2, by connecting a Riemann–Liouville differentiator of order 0.7 to y ( t ) , the C D 0 , t 0.3 y ( t ) signal can be defined.
Example 1.
Consider a linear fractional complex network (4) comprising two nodes. The weighted adjacency matrix, control matrix, and diagonal matrix are presented as
A = 0 1 2 0 0 , B = 0 2 , C = 0 0 0 1 2 ,
respectively. In the light of Theorem 1, the fractional complex network (4) with linear dynamic behavior can be controlled. Specifically, we used the Matlab/Simulink method. The block diagram of the simulation is depicted in Figure 3.
Here, the initial and final values of nodes v 1 and v 2 are taken as [ x 1 ( 0 ) , x 2 ( 0 ) ] T = [ 0 , 1 ] T and y 1 = [ x 1 ( 1 ) , x 2 ( 1 ) ] T = [ 2 + 0.3 k , 2 + 0.3 k ] T , k = 1 , 2 , , 5 , respectively. It is known that the Mittag–Leffler matrix function for matrix A + C can be written as
E 0.5 , 0.5 ( ( A + C ) t 0.5 ) = M 1 2 ( t ) M 1 ( t ) M 2 ( t ) M 1 ( t ) M 2 ( t ) M 2 2 ( t ) ,
where M 1 ( t ) = E 0.5 , 0.5 ( 0.5 t 0.5 ) 1 π and M 2 ( t ) = E 0.5 , 0.5 ( 0.5 t 0.5 ) . The Gramian matrix of the network (4) is
W F = 0 1 E 0.5 , 0.5 ( ( A + C ) ( 1 ς ) 0.5 ) B B * E 0.5 , 0.5 ( ( A + C ) * ( 1 ς ) 0.5 ) d ς = 4 0 1 M 1 2 ( 1 ς ) M 1 ( 1 ς ) M 2 ( 1 ς ) M 1 ( 1 ς ) M 2 ( 1 ς ) M 2 2 ( 1 ς ) d ς = 1.4507 2.6830 2.6830 5.1886 .
Using Matlab/Simulink, the controlled trajectories x 1 ( t ) , x 2 ( t ) and steering control
u ( t ) = ( 1 t ) 0.5 B * E 0.5 , 0.5 ( A * ( 1 t ) 0.5 ) W F 1 [ y 1 E 0.5 ( A ) x 0 ]
are computed and are depicted in Figure 4 and Figure 5. Then, we can see that a linear fractional complex network with two nodes can be steered from the initial value [ 0 , 1 ] T to the desired value y 1 .
Example 2.
Consider a nonlinear fractional complex network (16) comprising two nodes, where weighted adjacency matrix, control matrix, and nonlinear function f ( x ( t ) ) satisfying conditions ( A 1 ) ( A 3 ) are given as follows:
A = 0 1 2 1 2 0 , B = 0 2 , f ( x ( t ) ) = 0 1 10 c o s ( x 2 ( t ) ) ,
respectively. According to Remark 1, the system ( A , B ) is controllable, which means that a nonlinear network (16) with a weighted adjacency matrix is controllable. To be specific, we used the Matlab/Simulink method. The block diagram of the simulation is depicted in Figure 6.
Here, the initial and final values of nodes v 1 and v 2 are taken as [ x 1 ( 0 ) , x 2 ( 0 ) ] T = [ 1 , 1 ] T and z 1 = [ x 1 ( 1 ) , x 2 ( 1 ) ] T = [ 2 + 0.5 k , 1 + 0.5 k ] T , k = 1 , 2 , , 5 , respectively. It is known that the Mittag–Leffler matrix function for a given matrix A can be written as
E 0.5 , 0.5 ( A t 0.5 ) = M 2 2 ( t ) M 1 ( t ) M 2 ( t ) M 1 ( t ) M 2 ( t ) M 1 2 ( t ) ,
where
M 1 ( t ) = 1 2 [ E 0.5 , 0.5 ( 0.5 t 0.5 ) + E 0.5 , 0.5 ( 0.5 t 0.5 ) ]
and
M 2 ( t ) = 1 2 [ E 0.5 , 0.5 ( 0.5 t 0.5 ) E 0.5 , 0.5 ( 0.5 t 0.5 ) ] .
The corresponding Gramian matrix of network (16) is
W F = 0 1 E 0.5 , 0.5 ( A ( 1 ς ) 0.5 ) B B * E 0.5 , 0.5 ( A * ( 1 ς ) 0.5 ) d ς = 4 0 1 M 2 2 ( 1 ς ) M 1 ( 1 ς ) M 2 ( 1 ς ) M 1 ( 1 ς ) M 2 ( 1 ς ) M 1 2 ( 1 ς ) d ς = 0.7026 1.1814 1.1814 2.1232 .
Using Matlab/Simulink, the controlled trajectories x 1 ( t ) , x 2 ( t ) and steering control
u ( t ) = ( 1 t ) 0.5 B * E 0.5 , 0.5 ( A * ( 1 t ) 0.5 ) W F 1 z 1 E 0.5 ( A ) x 0 0 1 ( 1 ς ) 0.5 E 0.5 , 0.5 ( A ( 1 ς ) 0.5 ) f ( x ( ς ) ) d ς ]
are computed and are depicted in Figure 7 and Figure 8. Then, we can see that the nonlinear network with two nodes can be steered from the initial value [ 1 , 1 ] T to the desired value z 1 .
Example 3.
Consider a nonlinear fractional directed network (27) with a Laplacian matrix. The Laplacian matrix, control matrix, and nonlinear function f ( x ( t ) ) satisfying conditions ( A 1 ) ( A 3 ) are given as
L = 0.5 0 0.5 0 , B = 1 0 , f ( x ( t ) ) = 1 100 c o s ( x 2 ( t ) ) 0 ,
respectively. On the basis of Remark 1, the system ( L , B ) is controllable, which indicates that a nonlinear directed network (27) with a Laplacian matrix is controllable. Concretely, we used the Matlab/Simulink method. The block diagram of the simulation is depicted in Figure 9.
Here, the initial and final values of nodes v 1 and v 2 are taken as [ x 1 ( 0 ) , x 2 ( 0 ) ] T = [ 1 , 0 ] T and w 1 = [ x 1 ( 1 ) , x 2 ( 1 ) ] T = [ 1 + 0.5 k , 3 + 0.5 k ] T , k = 1 , 2 , , 5 , respectively. It is known that the Mittag–Leffler matrix function for a given matrix L can be written as
E 0.5 , 0.5 ( L t 0.5 ) = M 1 2 ( t ) M 1 ( t ) M 2 ( t ) M 1 ( t ) M 2 ( t ) M 2 2 ( t ) ,
where M 1 ( t ) = E 0.5 , 0.5 ( 0.5 t 0.5 ) and M 2 ( t ) = E 0.5 , 0.5 ( 0.5 t 0.5 ) + 1 π . The corresponding Gramian matrix of the network (27) is
W F = 0 1 E 0.5 , 0.5 ( L ( 1 ς ) 0.5 ) B B * E 0.5 , 0.5 ( L * ( 1 ς ) 0.5 ) d ς = 0 1 M 1 2 ( 1 ς ) M 1 ( 1 ς ) M 2 ( 1 ς ) M 1 ( 1 ς ) M 2 ( 1 ς ) M 2 2 ( 1 ς ) d ς = 0.1157 0.0727 0.0727 0.0571 .
Using Matlab/Simulink, the controlled trajectories x 1 ( t ) , x 2 ( t ) and steering control
u ( t ) = ( 1 t ) 0.5 B * E 0.5 , 0.5 ( L * ( 1 t ) 0.5 ) W F 1 w 1 E 0.5 ( L ) x 0 0 1 ( 1 ς ) 0.5 E 0.5 , 0.5 ( L ( 1 ς ) 0.5 ) f ( x ( ς ) ) d ς ]
are computed and are depicted in Figure 10 and Figure 11. Then, we can see that a nonlinear fractional complex network with two nodes can be steered from the initial value [ 1 , 0 ] T to the desired value w 1 .

6. Conclusions

Considering the fact that controllability is one of the most basic problems in the field of fractional complex network control and complex networks are typically nonlinear in real-world applications, this paper investigates the controllability problems of nonlinear fractional complex networks with a weighted adjacency matrix or Laplacian matrix. What is worth mentioning is that all controllability results obtained are proved by means of rigorous mathematical theory. Our findings indicate that the fractional complex networks with linear dynamic behavior can be controlled if the controllability matrix is of full rank. More importantly, the controllability criteria of nonlinear fractional complex networks are deduced, which provides a theoretical framework for the controllability analysis. Finally, three corresponding examples clearly show that it is convenient and efficient to obtain the controllability of given networks by using a newly developed technique. It is believed that this manuscript will play a crucial role in the controllability analysis and controller design of nonlinear fractional complex networks. In the future, we will continue to focus on the following interesting topics:
Solve the other control problems like optimal control, approximate controllability, etc.;
Develop the controllability and observability on complex fractional time-varying systems;
Implement controllability for complex time-varying systems numerically using the Matlab/Simulink method.

Author Contributions

Conceptualization, W.M. and X.B.; methodology, X.B; writing—original draft preparation, X.B.; writing—review and editing, W.M. and X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Natural Science Foundation of Gansu Province (No. 22JR5RA184) and the Fundamental Research Funds for the Central Universities (Grant No. 31920230174).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cui, L.Y.; Kumara, S.; Albert, R. Complex networks: An engineering view. IEEE Circuits Syst. Mag. 2010, 10, 10–25. [Google Scholar] [CrossRef]
  2. MacArthur, B.D.; Sánchez-García, R.J.; Anderson, J.W. Symmetry in complex networks. Discret. Appl. Math. 2008, 156, 3525–3531. [Google Scholar] [CrossRef]
  3. Strogatz, S.H. Exploring complex networks. Nature 2001, 410, 268–276. [Google Scholar] [CrossRef]
  4. Barabási, A.L.; Oltvai, Z.N. Network biology: Understanding the cell’s functional organization. Nat. Rev. Genet. 2004, 5, 101–113. [Google Scholar] [CrossRef]
  5. Yuan, Z.; Zhao, C.; Di, Z.; Wang, W.-X.; Lai, Y.-C. Exact controllability of complex networks. Nat. Commun. 2013, 4, 2447. [Google Scholar] [CrossRef]
  6. Balachandran, K.; Govindaraj, V.; Rivero, M.; Trujillo, J.J. Controllability of fractional damped dynamical systems. Appl. Math. Comput. 2015, 257, 66–73. [Google Scholar] [CrossRef]
  7. Xiang, L.; Chen, F.; Ren, W.; Chen, G. Advances in network controllability. IEEE Circuits Syst. Mag. 2019, 19, 8–32. [Google Scholar] [CrossRef]
  8. Liu, B.; Chu, T.; Wang, L.; Xie, G. Controllability of a leader-follower dynamic network with switching topology. IEEE Trans. Autom. Control 2008, 53, 1009–1013. [Google Scholar] [CrossRef]
  9. Cai, N. On quantitatively measuring controllability of complex networks. Phys. A 2017, 474, 282–292. [Google Scholar] [CrossRef]
  10. Chen, G. Pinning control and controllability of complex dynamical networks. Int. J. Autom. Comput. 2017, 14, 1–9. [Google Scholar] [CrossRef]
  11. Whalen, A.J.; Brennan, S.N.; Sauer, T.D.; Schiff, S.J. Observability and controllability of nonlinear networks: The role of symmetry. Phys. Rev. X 2015, 5, 011005. [Google Scholar] [CrossRef]
  12. Liu, Y.Y.; Slotine, J.J.; Barabási, A.L. Controllability of complex networks. Nature 2011, 473, 167–173. [Google Scholar] [CrossRef]
  13. Ma, W.; Ma, N.; Dai, C.; Chen, Y.Q.; Wang, X. Fractional modeling and optimal control strategies for mutated COVID-19 pandemic. Math. Methods Appl. Sci. 2022, 1–25. [Google Scholar] [CrossRef]
  14. Sun, H.G.; Zhang, Y.; Baleanu, D.; Chen, W.; Chen, Y.Q. A new collection of real world applications of fractional calculus in science and engineering. Commun. Nonlinear Sci. Numer. Simul. 2018, 64, 213–231. [Google Scholar] [CrossRef]
  15. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equation; John Wiley and Sons: New York, NY, USA, 1993. [Google Scholar]
  16. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  17. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  18. Li, C.P.; Deng, W.H.; Xu, D. Chaos synchronization of the Chua system with a fractional order. Phys. A 2006, 360, 171–185. [Google Scholar] [CrossRef]
  19. Monje, C.A.; Chen, Y.Q.; Vinagre, B.M.; Xue, D.Y.; Feliu, V. Fractional-Order Systems and Controls; Springer-Verlag: London, UK, 2010. [Google Scholar]
  20. Ma, W.; Li, Z.; Ma, N. Synchronization of discrete fractional-order complex networks with and without unknown topology. Chaos 2022, 32, 013112. [Google Scholar] [CrossRef]
  21. Chen, Y.Q.; Ahn, H.S.; Xue, D.Y. Robust controllability of interval fractional order linear time invariant systems. Signal Process. 2006, 86, 2794–2802. [Google Scholar] [CrossRef]
  22. Balachandran, K.; Govindaraj, V.; Rodríguez-Germá, L.; Trujillo, J.J. Controllability results for nonlinear fractional order dynamical systems. J. Optim. Theory Appl. 2013, 156, 33–44. [Google Scholar] [CrossRef]
  23. Liu, J.; Liu, S.; Li, H. Controllability result of nonlinear higher order fractional damped dynamical system. J. Nonlinear Sci. Appl. 2017, 10, 325–337. [Google Scholar] [CrossRef]
  24. Ahmad, B.; Nieto, J.J. Existence results for a coupled system of nonlinear fractional differential equations with three-point boundary conditions. Comput. Math. Appl. 2009, 58, 1838–1843. [Google Scholar] [CrossRef]
  25. Govindaraj, V.; Balachandran, K.; George, R.K. Numerical Approach for the Controllability of Composite Fractional Dynamical Systems. J. Appl. Nonlinear Dyn. 2018, 7, 59–72. [Google Scholar] [CrossRef]
  26. Zhang, H.; Chen, D.Y.; Xu, B.B.; Zhang, R.F. Controllability of fractional-order directed complex networks. Mod. Phys. Lett. B 2014, 28, 1450211. [Google Scholar] [CrossRef]
  27. Zhang, H.; Chen, D.; Xu, B.B.; Zhou, R. Controllability of fractional-order directed complex networks with self loop and double edge structure. J. Circuits Syst. Comput. 2015, 24, 1550087. [Google Scholar] [CrossRef]
  28. Wang, L.; Chen, G.; Wang, X.; Tang, W.K.S. Controllability of networked MIMO systems. Automatica 2016, 69, 405–409. [Google Scholar] [CrossRef]
  29. García, P.M.I. Analyzing controllability of neural networks. Wseas Trans. Circuits Syst. 2019, 18, 1–6. [Google Scholar]
  30. Biggs, N.L. Algebraic Graph Theory; Cambridge University Press: Cambridge, UK, 1993. [Google Scholar]
  31. Li, C.P.; Deng, W.H. Remarks on fractional derivatives. Appl. Math. Comput. 2007, 187, 777–784. [Google Scholar] [CrossRef]
  32. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: San Diego, CA, USA, 1998. [Google Scholar]
  33. Zeidler, E. Nonlinear Functional Analysis and Its Applications I: Fixed-Point Theorems; Springer: New York, NY, USA, 1986. [Google Scholar]
  34. Govindaraj, V.; George, R.K. Controllability of fractional dynamical systems: A functional analytic approach. Math. Control Relat. Fields 2017, 7, 537. [Google Scholar] [CrossRef]
  35. Eǧecioǧlu, Ö.; Garsia, A.M. Lessons in Enumerative Combinatorics; Springer: Cham, Switzerland, 2021. [Google Scholar]
  36. Xue, D.Y. FOTF Toolbox. MATLAB Central File Exchange. 2023. Available online: https://www.mathworks.com/matlabcentral/fileexchange/60874-fotf-toolbox (accessed on 2 October 2023).
Figure 1. Simulink block diagram of the control u.
Figure 1. Simulink block diagram of the control u.
Fractalfract 08 00043 g001
Figure 2. The output y ( t ) signal and its Caputo derivative of order 0.3 .
Figure 2. The output y ( t ) signal and its Caputo derivative of order 0.3 .
Fractalfract 08 00043 g002
Figure 3. Matlab/Simulink block diagram for the network (4) in Example 1.
Figure 3. Matlab/Simulink block diagram for the network (4) in Example 1.
Fractalfract 08 00043 g003
Figure 4. The state trajectories of the linear fractional complex network in Example 1.
Figure 4. The state trajectories of the linear fractional complex network in Example 1.
Fractalfract 08 00043 g004
Figure 5. The trajectory of outer controller u ( t ) imposed on node v 2 in Example 1.
Figure 5. The trajectory of outer controller u ( t ) imposed on node v 2 in Example 1.
Fractalfract 08 00043 g005
Figure 6. Matlab/Simulink block diagram for the network (16) in Example 2.
Figure 6. Matlab/Simulink block diagram for the network (16) in Example 2.
Fractalfract 08 00043 g006
Figure 7. The state trajectories of a nonlinear fractional complex network with a weighted adjacency matrix in Example 2.
Figure 7. The state trajectories of a nonlinear fractional complex network with a weighted adjacency matrix in Example 2.
Fractalfract 08 00043 g007
Figure 8. The trajectories of outer controller u ( t ) imposed on node v 2 in Example 2.
Figure 8. The trajectories of outer controller u ( t ) imposed on node v 2 in Example 2.
Fractalfract 08 00043 g008
Figure 9. Matlab/Simulink block diagram for the network (27) in Example 3.
Figure 9. Matlab/Simulink block diagram for the network (27) in Example 3.
Fractalfract 08 00043 g009
Figure 10. The state trajectories of a nonlinear fractional complex network with Laplacian matrix in Example 3.
Figure 10. The state trajectories of a nonlinear fractional complex network with Laplacian matrix in Example 3.
Fractalfract 08 00043 g010
Figure 11. The trajectories of outer controller u ( t ) imposed on node v 1 in Example 3.
Figure 11. The trajectories of outer controller u ( t ) imposed on node v 1 in Example 3.
Fractalfract 08 00043 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bao, X.; Ma, W.; Li, X. Controllability of Fractional Complex Networks. Fractal Fract. 2024, 8, 43. https://doi.org/10.3390/fractalfract8010043

AMA Style

Bao X, Ma W, Li X. Controllability of Fractional Complex Networks. Fractal and Fractional. 2024; 8(1):43. https://doi.org/10.3390/fractalfract8010043

Chicago/Turabian Style

Bao, Xionggai, Weiyuan Ma, and Xin Li. 2024. "Controllability of Fractional Complex Networks" Fractal and Fractional 8, no. 1: 43. https://doi.org/10.3390/fractalfract8010043

APA Style

Bao, X., Ma, W., & Li, X. (2024). Controllability of Fractional Complex Networks. Fractal and Fractional, 8(1), 43. https://doi.org/10.3390/fractalfract8010043

Article Metrics

Back to TopTop