# Modified Projective Outer Synchronization between Two Fractional Order Complex Networks via Adaptive Control

Authors: Hong-juan Liu; Guo-qi Liu; Hai Yu; Zhi-liang Zhu
Journal: Discrete Dynamics in Nature and Society
 (2014)

---

## Abstract

The modified projective outer synchronization between two different fractional order complex dynamical networks with different node dynamics and different topological structures is investigated in this paper. Based on the stability theorem of fractional order systems, the adaptive controllers and sufficient criteria for achieving the synchronization are given by applying the nonlinear control. The proposed method is also applicable to modified projective outer synchronization between two fractional order complex networks with identical topological structures. In addition, there are fewer restrictions on the coupling configuration matrix and inner-coupling matrix. The corresponding numerical results coincide with theoretical analysis, which also show that feedback strength and fractional order can effectively affect the synchronization performance.


---

## Body

## 1. Introduction

Since the pioneering work of Pecora and Carroll synchronizing identical chaotic systems with different initial conditions [1], various synchronization scenarios have been studied for chaotic systems [2–4]. In recent years, research on complex networks has attracted increasing attention due to its wide and potential applications in various fields. A complex dynamical network consists of a set of interconnected nodes, which are dynamic systems with specific contents. The network model widely exists in nature and society, for example, the World Wide Web, Internet, social networks, biological neural networks, food webs, proteins, the transportation networks, and so forth [5–8]. With the development of the property and dynamics of complex networks, the chaos synchronization of complex networks has become a research focus owing to its engineering applications in secure communications, neural networks, biological systems, parallel image processing, and so forth. Therefore, much work has been done for the synchronization of complex networks in the literature [9–12].The above research results have laid a good foundation for further theoretical research and practical application of synchronization dynamics in complex network. However, existing research efforts mostly focused on synchronization of integer order complex networks (ICNs), whose dynamics are described by integer order differential equations. In fact, fractional order systems provide an excellent instrument for the description of memory and hereditary properties of various materials and processes. Moreover, they include traditional integer order derivatives as a special case. It would be much better if many practical problems are described by fractional order dynamical systems rather than integer order ones. On the other hand, due to the limited theories for analyzing the dynamics of fractional order systems, synchronization of fractional order complex systems (FCNs) is still a challenging research topic and studies about it are still few [13–16]. In [13], authors studied the synchronized motions in a star network of coupled fractional order systems in which the major element is coupled to each of the noninteracting individual elements. In [14], Chen et al. discussed cluster synchronization of complex dynamical networks with fractional order dynamical nodes. Only the nodes in one community which have direct connections to the nodes in other communities are needed to be controlled, resulting in reduced control cost. In [15], projective synchronization of a complex network is studied and the coupling functions of the connected nodes are identified. In [16], Wong et al. investigated robust synchronization of fractional order complex dynamical networks with parametric uncertainties. Based on the properties of the Kronecker product and the stability of the fractional order system, the robust synchronization criteria are derived by applying the nonlinear control.The synchronization among the nodes of a single network has been termed “inner synchronization” in the above literatures. In this case, it has been shown that the stability of synchronized states depends on the details of the underlying network topology. In contrast, the synchronization between the corresponding nodes of two coupled networks is called “outer synchronization.” In the real world, there are a great many examples about relationships between different networks, such as the original spreading of infectious diseases between two communities, the balance of beneficial bacteria and pathogenic bacteria in our digestive systems, and predator-prey interactions in ecological systems. These examples show that it is necessary and significant to investigate the dynamics between two coupled networks. At present, the outer synchronization between coupled ICNs has been studied extensively [17–19]. However, to the best of our knowledge, few papers investigate the synchronization between two complex networks with fractional order nodes and the only four attempts are the works in [20–23]. By designing nonlinear controllers, Wu and Lu first reported outer synchronization phenomena between two coupled FCNs [20]. Also, Asheghan et al. considered robust outer synchronization between two FCNs coupled by an open-plus-closed-loop scheme [21]. More recently, authors in [22] address the global outer synchronization problem between two FCNs coupled in a drive-response configuration. In particular, for a given FCN composed of Lur’e systems, an observer-type response network with nonfragile output feedback controllers is constructed. Authors in [23] propose a novel projective outer synchronization between unidirectionally coupled uncertain FCNs through scalar transmitted signals. However, among these references, only the work in [23] considers projective synchronization between two FCNs with identical nodes and identical topological structures. Thus, it is necessary to construct a more general projective synchronization scheme between two coupled FCNs.In this paper, we discussed the modified projective outer synchronization (MPOS) problem between two FCNs with different node dynamics and different topological structures. By using adaptive control scheme, a nonlinear controller is designed to achieve the MPOS. The controller is also applicable to the MPOS of two FCNs with identical topological structures. Based on the stability theory of fractional order systems, basic criteria in the format of linear matrix inequalities (LMIs) for MPOS are proposed. Finally, the effectiveness and feasibility of the designed control strategy for MPOS are demonstrated by numerical simulations on chaotic drive-response FCNs. These simulations indicate that increasing the fractional order and the feedback strength will accelerate the speed of outer synchronization.


## 2. Model Description and Preliminaries

There are several definitions of fractional derivatives. The Caputo derivative is more popular in the real applications [24], which is defined as
(1)
d
q
f
(
t
)
d
t
q
≡
D
q
f
(
t
)
=
1
Γ
(
m
-
q
)
∫
0
t
(
t
-
τ
)
m
-
q
-
1
f
(
m
)
(
τ
)
d
τ
,
where m is the smallest integer larger than q, D
q denotes the Caputo definition of the fractional derivative, f
(
m
)
(
t
) is the m-order derivative in the usual sense, and Γ stands for gamma function.Now, consider a class of complex networks consisting ofN identical coupled nodes, in which each node is an n-dimensional fractional order chaotic system. The standard derivative is replaced by a fractional derivative as the following:
(2)
D
q
x
i
(
t
)
=
f
(
x
i
(
t
)
)
+
a
∑
j
=
1
N
c
i
j
A
x
j
(
t
)
,
i
=
1,2
,
…
,
N
,
where x
i
∈
R
n is the state vector of the ith node, f
∈
R
n is a continuous differentiable vector function that describes the dynamics of the individual node, a
>
0 is the coupling strength, A is an inner-coupling matrix satisfying A
>
0, and C
=
(
c
i
j
)
N
×
N is the coupling configuration matrix representing the coupling strength and the topological structure of the network, in which c
i
j is defined as follows: if there is a connection from node i to node j
(
j
≠
i
), c
i
j
>
0; otherwise, c
i
j
=
0
(
j
≠
i
). The diagonal elements of matrix C are defined as
(3)
c
i
i
=
-
∑
j
=
1
,
j
≠
i
N
c
i
j
,
i
=
1,2
,
…
,
N
.In order to investigate outer synchronization between two different networks, we consider (2) as the master network and assume the response complex dynamical network containing N dynamical nodes as follows:
(4)
D
q
y
i
(
t
)
=
g
(
y
i
(
t
)
)
+
a
∑
j
=
1
N
d
i
j
A
y
j
(
t
)
+
u
i
(
t
)
,
i
=
1,2
,
…
,
N
,
where y
i
∈
R
n is the state vector of the node i, g
∈
R
n is a continuous differentiable vector function that describes the dynamics of the individual node, a
>
0 is the coupling strength, A has the same meaning as that in (2), D
=
(
d
i
j
)
N
×
N is the coupling configuration matrix which has the same meaning as matrix C, and u
i
(
t
) is the controller for node i to be designed according to the specific node dynamics and topological structures of the drive and response networks. Here, the coupling configuration matrices C and D are not necessarily symmetric or irreducible. We also do not impose any restriction on the inner-coupling matrix A.Next, according to the definition of modified projective synchronization between two chaotic systems, we will give the definition of MPOS between the drive FCN (2) and the response FCN (4).Definition 1 (MPOS).
For the drive network (2) and response network (4), it is said to be modified projective outer synchronization (MPOS) respect to the scaling factor matrix α
=
diag
⁡
(
α
1
,
α
2
,
…
,
α
N
), if there exists a controller u
i
(
t
), such that
(5)
lim
⁡
t
→
+
∞
∥
e
i
(
t
)
∥
=
lim
⁡
t
→
+
∞
∥
y
i
(
t
)
-
α
i
x
i
(
t
)
∥
=
0
,
i
=
1,2
,
…
,
N
.
always holds.Remark 2.
For unidirectional coupled FCN (2) and (4), we can see from Definition 1 that MPOS means that the corresponding node state vectors x
i
(
t
) and y
i
(
t
) could achieve synchronization with respect to the scaling factor α
i regardless of synchronization of inner networks. Therefore, the two networks may contain isolated clusters or nodes. The coupling configuration matrix is not necessary to be irreducible.Remark 3.
It is easy to see that arbitrary constant matrixα can be tuned to realize various outer synchronization behaviors between states x
i
(
t
) and y
i
(
t
). For example, if α
1
=
α
2
=
⋯
=
α
N
=
1, we can observe complete outer synchronization (COS); when α
1
=
α
2
=
⋯
=
α
N
=
-
1, antiphase out synchronization (AOS) is obtained. Moreover, we can also realize amplification or attenuation of the oscillation in the drive FCN (2) relative to that in the response FCN (4) by choosing |
α
i
|
>
1
,
i
=
1,2
,
…
,
N or |
α
i
|
<
1
,
i
=
1,2
,
…
,
N, respectively.Lemma 4 (see [25]).
For a given autonomous fractional order linear system(6)
D
q
x
(
t
)
=
f
(
x
)
=
W
x
,
where 0
<
q
<
1, x
∈
R
n is a state vector, and W
∈
R
n
×
n is a constant matrix, the system is asymptotically stable if and only if |
arg
(
λ
i
)
|
>
q
π
/
2
,
i
=
1,2
,
…
,
n, where arg
(
λ
i
) denotes the argument of the eigenvalue λ
i of  W. Also, system (6) is stable if and only if either it is asymptotically stable or those critical eigenvalues which satisfy |
arg
(
λ
i
)
|
=
q
π
/
2 have a geometric multiplicity one.Theorem 5 (see [26]).
For a class of autonomous fractional order nonlinear system(7)
D
q
x
(
t
)
=
f
(
x
)
=
W
(
x
)
x
,
where 0
<
q
<
1, x
∈
R
n is a state vector, and W
(
x
)
∈
R
n
×
n is a continuous coefficient matrix, which contains the state variables. If there exists a real symmetric positive definite matrix P
>
0, such that the equation
(8)
J
=
x
T
P
D
q
x
≤
0
always holds for any state x
(
t
), then the trivial solution x
=
0 of system (7) is asymptotically stable.Proof.
Obviously,x
=
0 is an equilibrium point of system (7) from f
(
x
)
=
W
(
x
)
x. Because J
=
x
T
P
D
q
x
≤
0, we have J
+
J
T
=
x
T
(
P
W
(
x
)
+
W
H
(
x
)
P
)
x
≤
0. Denoting P
W
(
x
)
+
W
H
(
x
)
P
=
Q
(
x
), we can conclude that P
W
(
0
)
+
W
H
(
0
)
P
=
Q
(
0
)
≤
0 always holds for x
=
0. Otherwise, if Q
(
0
)
>
0, by lim
⁡
x
→
0
Q
(
x
)
=
Q
(
0
), there must exist state x
*
≠
0, such that matrix Q
(
x
*
)
>
0; thus, we have x
*
T
(
P
W
(
x
*
)
+
W
H
(
x
*
)
P
)
x
*
=
x
*
T
Q
(
x
*
)
x
*
>
0, a contradiction. Therefore, the inequality P
W
(
0
)
+
W
H
(
0
)
P
=
Q
(
0
)
≤
0 always holds for x
=
0.
Assume thatλ is one of the eigenvalues of matrix W
(
0
) and the corresponding nonzero eigenvector is ξ; that is, W
(
0
)
ξ
=
λ
ξ. Multiply the left of P
W
(
0
)
+
W
H
(
0
)
P
=
Q
(
0
) by ξ
T and its right by ξ. We derive that
(9)
ξ
T
(
P
W
(
0
)
+
W
H
(
0
)
P
)
ξ
=
ξ
T
Q
(
0
)
ξ
≤
0
.
BecauseP
>
0, we can obtain λ
+
λ
¯
≤
0. Thus, we have |
arg
(
λ
)
|
≥
π
/
2
>
q
π
/
2 for 0
<
q
<
1. According to the stability theory of fractional-order nonlinear systems in Lemma 4, the trivial solution x
=
0 of system (7) is asymptotically stable.
The proof is completed.Assumption 6 (see [11]).
For the positive definite matrixA in network (2), there exists a constant matrix T such that
(10)
(
x
-
y
)
T
(
f
(
x
)
-
f
(
y
)
)
≤
(
x
-
y
)
T
T
A
(
x
-
y
)
,
∀
x
,
y
∈
R
n
.
always holds.
According to Definition1, the error vectors can be expressed as e
i
(
t
)
=
y
i
(
t
)
-
α
i
x
i
(
t
)
,
i
=
1,2
,
…
,
N. Subtracting network (4) from (2) gives the error dynamical network
(11)
D
q
e
i
(
t
)
=
g
(
y
i
(
t
)
)
+
a
∑
j
=
1
N
d
i
j
A
y
j
(
t
)
-
α
i
f
(
x
i
(
t
)
)
D
q
e
i
(
t
)
=
-
α
i
a
∑
j
=
1
N
c
i
j
A
x
j
(
t
)
+
u
i
(
t
)
,
i
=
1,2
,
…
,
N
.
Then the MPOS problem of drive network (2) and response network (4) is converted into the stability problem of the error dynamical network (11). Our target is to design some appropriate controllers to achieve lim
⁡
t
→
+
∞
∥
e
i
(
t
)
∥
=
lim
⁡
t
→
+
∞
∥
y
i
(
t
)
-
α
i
x
i
(
t
)
∥
=
0
,
i
=
1,2
,
…
,
N.


## 3. The Adaptive Controller for MPOS

In this section, we will study the MPOS between two FCNs with different topologies or identical topologies by using adaptive control method. We first provide the following main theorem about the MPOS between two FCNs with different topological structure.Theorem 7.
For the positive definite matrixA in network (2), suppose T satisfies Assumption 6 and let θ
=
∥
T
∥. If the following condition is satisfied,
(12)
θ
I
N
+
a
C
-
K
≤
0
,
then the solution e
=
0 of error dynamical network (11) is asymptotically stable under the following controllers:
(13)
u
i
(
t
)
=
α
i
f
(
x
i
(
t
)
)
-
g
(
α
i
x
i
(
t
)
)
u
i
(
t
)
=
-
a
∑
j
=
1
N
b
i
j
A
y
j
(
t
)
-
k
i
A
e
i
(
t
)
,
i
=
1,2
,
…
,
N
,
where I
N is N-dimensional identity matrix, difference matrix B
=
(
b
i
j
)
N
×
N
=
D
-
C
=
(
d
i
j
-
c
i
j
)
N
×
N, feedback strengths k
i
(
i
=
1,2
,
…
,
N
) are positive constants, and constant matrix K
=
diag
⁡
(
k
1
,
k
2
,
…
,
k
N
).Proof.
Consider the followingJ function:
(14)
J
(
t
)
=
∑
i
=
1
N
e
i
T
(
t
)
D
q
e
i
(
t
)
.
After substituting (11) and (13) into J
(
t
), we have
(15)
J
(
t
)
=
∑
i
=
1
N
e
i
T
(
t
)
[
a
∑
j
=
1
N
c
i
j
A
e
i
(
t
)
-
k
i
A
e
i
(
t
)
g
(
y
i
(
t
)
)
-
g
(
α
i
x
i
(
t
)
)
XXXXXXXXXX
+
a
∑
j
=
1
N
c
i
j
A
e
i
(
t
)
-
k
i
A
e
i
(
t
)
]
J
(
t
)
=
∑
i
=
1
N
e
i
T
(
t
)
[
g
(
y
i
(
t
)
)
-
g
(
α
i
x
i
(
t
)
)
]
J
(
t
)
=
+
a
∑
i
=
1
N
e
i
T
(
t
)
∑
j
=
1
N
c
i
j
A
e
i
(
t
)
-
∑
i
=
1
N
e
i
T
(
t
)
k
i
A
e
i
(
t
)
J
(
t
)
≤
∑
i
=
1
N
e
i
T
(
t
)
T
A
e
i
(
t
)
+
a
∑
i
=
1
N
∑
j
=
1
N
e
i
T
(
t
)
c
i
j
A
e
i
(
t
)
J
(
t
)
=
-
∑
i
=
1
N
e
i
T
(
t
)
k
i
A
e
i
(
t
)
J
(
t
)
≤
e
T
(
t
)
[
(
θ
I
N
+
a
C
-
K
)
⊗
A
]
e
(
t
)
,
where e
(
t
)
=
[
e
1
T
(
t
)
,
e
2
T
(
t
)
,
…
,
e
N
T
(
t
)
]
T and K
=
diag
⁡
(
k
1
,
k
2
,
…
,
k
N
). Because A
>
0, it is easy to conclude that, if (12) is satisfied, the solution e
=
0 of error dynamical network (11) is asymptotically stable.
The proof is completed.From Theorem7, we can conclude that as long as feedback strengths k
i in feedback matrix K are big enough, the MPOS between network (2) and network (4) can be achieved.As a special case of different topological structure, MOPS between two FCNs with identical topology can be achieved under the similar adaptive control. In this case, we haveC
=
D and the error dynamical systems become
(16)
D
q
e
i
(
t
)
=
g
(
y
i
(
t
)
)
-
α
i
f
(
x
i
(
t
)
)
D
q
e
i
(
t
)
=
+
a
∑
j
=
1
N
c
i
j
A
(
y
j
(
t
)
-
α
i
x
j
(
t
)
)
D
q
e
i
(
t
)
=
+
u
i
(
t
)
,
i
=
1,2
,
…
,
N
.Based on Theorem7, we can easily obtain the following corollary.Corollary 8.
For two FCNs with identical topology, supposeT satisfies Assumption 6 and let θ
=
∥
T
∥. If condition (12) is satisfied, then the solution e
=
0 of error dynamical network (16) is asymptotically stable under the following controllers:
(17)
u
i
(
t
)
=
α
i
f
(
x
i
(
t
)
)
-
g
(
α
i
x
i
(
t
)
)
-
k
i
A
e
i
(
t
)
,
i
=
1,2
,
…
,
N
.Proof.
When the topological structures of FCN (2) and (4) are the same, that is, C
=
D, difference matrix B
=
(
b
i
j
)
N
×
N becomes a zero matrix. Therefore, the controllers in (13) turn into the controllers described by (17). While substituting (16) and (17) into J function described by (14), the rest of the proof process is the same as that of Theorem 7.


## 4. Numerical Simulations

In this section, illustrative examples are given to verify the effectiveness of the theoretical criteria obtained in the preceding section. Some further analyses about the influence of fractional orderq and feedback strength k
i on MPOS are provided. For this purpose, we consider fractional order Chen system as the node of drive network and consider fractional order financial system as the node of response network.


### 4.1. The Verification of Theorem7

The well-known fractional order Chen chaotic system [27] is described by
(18)
D
q
x
1
=
a
(
x
2
-
x
1
)
,
D
q
x
2
=
(
c
-
a
)
x
1
-
x
1
x
3
+
c
x
2
,
D
q
x
3
=
x
1
x
2
-
b
x
3
.
When fractional order q holds q
∈
[
0.7,1
] and parameters are selected as a
=
35
,
b
=
3
,
and
c
=
28, the fractional order Chen system behaves chaotically. Figure 1 displays the attractor of the fractional order Chen chaotic system with q
=
0.96.Figure 1
The chaotic attractor of fractional order Chen system withq
=
0.96.Recently, a dynamic model of finance was reported in [28]. The time variation of interest rate y
1, investment demand y
2, and price index y
3 is described in the model. The fractional order financial chaotic system is given by
(19)
D
q
y
1
=
y
3
+
(
y
2
-
a
)
y
1
,
D
q
y
2
=
1
-
b
y
2
-
y
1
2
,
D
q
y
3
=
-
y
1
-
c
y
3
.The fractional order financial system is chaotic when fractional orderq changes in [
0.85,1
] and parameters are selected as a
=
3
,
b
=
0.1
,
and
c
=
1. The attractor of the fractional order financial chaotic system with q
=
0.96 is shown in Figure 2.Figure 2
The chaotic attractor of fractional order financial system withq
=
0.96.In the following, we will take the fractional order Chen chaotic system and the fractional order financial chaotic system as node dynamics to illustrate the proposed method for the MPOS. In all the simulations, let the fractional orders of the drive and the response networks beq
=
0.96. And a predictor-corrector scheme [29] is chosen to solve the fractional order differential equations. The network size is taken as N
=
5. Choose the coupling configuration matrices of network (2) and (4) as follows:
(20)
C
=
(
-
3
0
1
1
1
0
-
2
0
1
1
1
0
-
2
1
0
1
1
1
-
3
0
1
1
0
0
-
2
)
,
D
=
(
-
4
1
1
1
1
1
-
3
1
1
0
1
1
-
2
0
0
1
1
0
-
2
0
1
0
0
0
-
1
)
.For simplicity, the inner-coupling matrix is chosen asA
=
I
3. The initial values of the state vectors x
i
(
0
) and y
i
(
0
)
(
i
=
1,2
,
…
,
N
) are chosen arbitrarily in (
-
10,10
). The scaling factor matrix is chosen as α
=
diag
⁡
(
1
,
-
2,3
,
2
,
-
1
). The MPOS errors between two networks are defined as e
i
n
(
t
)
=
y
i
n
(
t
)
-
α
i
x
i
n
(
t
)
(
i
=
1,2
,
…
,
N
;
n
=
1,2
,
3
).According to Theorem7, condition (12) is only a sufficient condition for the MPOS between complex networks (2) and (4). Here, we give a choice of feedback matrix K
=
diag
⁡
(
20,15,5
,
16,10
) for realizing the MPOS, whose simulation results are shown in Figure 3. As seen in Figure 3, the trajectories of the synchronization errors e
i
n
(
i
=
1,2
,
…
,
10
;
n
=
1,2
,
3
) approach zero, which imply that the MPOS between networks (2) and (4) has been achieved with the adaptive controllers (13).Time evolutions of MPOS synchronization errore
i
n with time.
(a)
(b)
(c)


### 4.2. Further Analyses

Research results suggest that the values of fractional orderq and feedback strength k
i have a great influence on the synchronization effect. Some further analyses about synchronization effect are presented in this section.Here, in order to know the variations of all the MPOS errors, the norm of the total MPOS errors is defined by(21)
E
(
t
)
=
∑
i
=
1
N
[
e
i
1
2
(
t
)
+
e
i
2
2
(
t
)
+
e
i
3
2
(
t
)
]
N
.It is obvious that, whenE
(
t
) no longer increases, two fractional-order networks achieve the MPOS globally.First, we consider that feedback matrixK
=
diag
⁡
(
20,15,5
,
16,10
) is invariant and only fractional order q is changed. The values of synchronization error E
(
t
) with the increase of the fractional order q are described in Figure 4. From Figure 4, one can easily find that the synchronization effect becomes better with the increase of fractional order q.Time evolutions of synchronization errorE
(
t
) with different values of fractional order q.
(a)
(b)
(c)Second, let control gainK be changed while fractional order q maintains the same. In order to simplify the problem, let k
1
=
k
2
=
k
3
=
k
4
=
k
5 in current simulations. Figure 5 describes the time evolutions of MPOS error E
(
t
) with the fractional order q
=
0.96 while the feedback strength k
i is increasing from 5 to 35.Time evolutions of synchronization errorE
(
t
) with different values of k
i.
(a)
(b)It is seen from Figure5 that the synchronization effect turns better with the increase of the feedback strength k
i. Moreover, even if some of the k
i decrease to 0, the MPOS can still be achieved in the case of reducing synchronization speed. When the feedback strengths are selected as k
1
=
k
2
=
k
3
=
10
,
k
4
=
k
5
=
0 and k
1
=
10
,
k
2
=
k
3
=
k
4
=
k
5
=
0, the time evolutions of MPOS error E
(
t
) with the fractional order q
=
0.96 are illustrated in Figure 6.Time evolutions of synchronization errorE
(
t
) with some of k
i equal to 0.
(a)
(b)As mentioned below, different fractional orders and feedback strengths will lead to different rates towards synchronization. The appropriately large fractional order and feedback strength will lead to a fast outer synchronization. Furthermore, the MPOS can also be realized with some of feedback strengths equal to 0, although the speed of the synchronization will be much slower.


## 4.1. The Verification of Theorem7

The well-known fractional order Chen chaotic system [27] is described by
(18)
D
q
x
1
=
a
(
x
2
-
x
1
)
,
D
q
x
2
=
(
c
-
a
)
x
1
-
x
1
x
3
+
c
x
2
,
D
q
x
3
=
x
1
x
2
-
b
x
3
.
When fractional order q holds q
∈
[
0.7,1
] and parameters are selected as a
=
35
,
b
=
3
,
and
c
=
28, the fractional order Chen system behaves chaotically. Figure 1 displays the attractor of the fractional order Chen chaotic system with q
=
0.96.Figure 1
The chaotic attractor of fractional order Chen system withq
=
0.96.Recently, a dynamic model of finance was reported in [28]. The time variation of interest rate y
1, investment demand y
2, and price index y
3 is described in the model. The fractional order financial chaotic system is given by
(19)
D
q
y
1
=
y
3
+
(
y
2
-
a
)
y
1
,
D
q
y
2
=
1
-
b
y
2
-
y
1
2
,
D
q
y
3
=
-
y
1
-
c
y
3
.The fractional order financial system is chaotic when fractional orderq changes in [
0.85,1
] and parameters are selected as a
=
3
,
b
=
0.1
,
and
c
=
1. The attractor of the fractional order financial chaotic system with q
=
0.96 is shown in Figure 2.Figure 2
The chaotic attractor of fractional order financial system withq
=
0.96.In the following, we will take the fractional order Chen chaotic system and the fractional order financial chaotic system as node dynamics to illustrate the proposed method for the MPOS. In all the simulations, let the fractional orders of the drive and the response networks beq
=
0.96. And a predictor-corrector scheme [29] is chosen to solve the fractional order differential equations. The network size is taken as N
=
5. Choose the coupling configuration matrices of network (2) and (4) as follows:
(20)
C
=
(
-
3
0
1
1
1
0
-
2
0
1
1
1
0
-
2
1
0
1
1
1
-
3
0
1
1
0
0
-
2
)
,
D
=
(
-
4
1
1
1
1
1
-
3
1
1
0
1
1
-
2
0
0
1
1
0
-
2
0
1
0
0
0
-
1
)
.For simplicity, the inner-coupling matrix is chosen asA
=
I
3. The initial values of the state vectors x
i
(
0
) and y
i
(
0
)
(
i
=
1,2
,
…
,
N
) are chosen arbitrarily in (
-
10,10
). The scaling factor matrix is chosen as α
=
diag
⁡
(
1
,
-
2,3
,
2
,
-
1
). The MPOS errors between two networks are defined as e
i
n
(
t
)
=
y
i
n
(
t
)
-
α
i
x
i
n
(
t
)
(
i
=
1,2
,
…
,
N
;
n
=
1,2
,
3
).According to Theorem7, condition (12) is only a sufficient condition for the MPOS between complex networks (2) and (4). Here, we give a choice of feedback matrix K
=
diag
⁡
(
20,15,5
,
16,10
) for realizing the MPOS, whose simulation results are shown in Figure 3. As seen in Figure 3, the trajectories of the synchronization errors e
i
n
(
i
=
1,2
,
…
,
10
;
n
=
1,2
,
3
) approach zero, which imply that the MPOS between networks (2) and (4) has been achieved with the adaptive controllers (13).Time evolutions of MPOS synchronization errore
i
n with time.
(a)
(b)
(c)


## 4.2. Further Analyses

Research results suggest that the values of fractional orderq and feedback strength k
i have a great influence on the synchronization effect. Some further analyses about synchronization effect are presented in this section.Here, in order to know the variations of all the MPOS errors, the norm of the total MPOS errors is defined by(21)
E
(
t
)
=
∑
i
=
1
N
[
e
i
1
2
(
t
)
+
e
i
2
2
(
t
)
+
e
i
3
2
(
t
)
]
N
.It is obvious that, whenE
(
t
) no longer increases, two fractional-order networks achieve the MPOS globally.First, we consider that feedback matrixK
=
diag
⁡
(
20,15,5
,
16,10
) is invariant and only fractional order q is changed. The values of synchronization error E
(
t
) with the increase of the fractional order q are described in Figure 4. From Figure 4, one can easily find that the synchronization effect becomes better with the increase of fractional order q.Time evolutions of synchronization errorE
(
t
) with different values of fractional order q.
(a)
(b)
(c)Second, let control gainK be changed while fractional order q maintains the same. In order to simplify the problem, let k
1
=
k
2
=
k
3
=
k
4
=
k
5 in current simulations. Figure 5 describes the time evolutions of MPOS error E
(
t
) with the fractional order q
=
0.96 while the feedback strength k
i is increasing from 5 to 35.Time evolutions of synchronization errorE
(
t
) with different values of k
i.
(a)
(b)It is seen from Figure5 that the synchronization effect turns better with the increase of the feedback strength k
i. Moreover, even if some of the k
i decrease to 0, the MPOS can still be achieved in the case of reducing synchronization speed. When the feedback strengths are selected as k
1
=
k
2
=
k
3
=
10
,
k
4
=
k
5
=
0 and k
1
=
10
,
k
2
=
k
3
=
k
4
=
k
5
=
0, the time evolutions of MPOS error E
(
t
) with the fractional order q
=
0.96 are illustrated in Figure 6.Time evolutions of synchronization errorE
(
t
) with some of k
i equal to 0.
(a)
(b)As mentioned below, different fractional orders and feedback strengths will lead to different rates towards synchronization. The appropriately large fractional order and feedback strength will lead to a fast outer synchronization. Furthermore, the MPOS can also be realized with some of feedback strengths equal to 0, although the speed of the synchronization will be much slower.


## 5. Conclusions

In summary, the MPOS between two FCNs with different node dynamics and different topological structures is discussed. In the MPOS regime, we can observe different outer synchronization patterns by choosing the scaling factor, for example, complete outer synchronization, antiouter synchronization, and amplitude death. Based on the stability theory of the fractional order system and using the nonlinear control coupling schemes, some sufficient criteria are obtained and the suitable controllers are designed. The sufficient criteria and controllers are also applicable to the case with the same topological structures. This technique has been successfully applied to two fractional order complex networks, which consist of five fractional order Chen chaotic systems and five fractional order financial chaotic systems. Numerical evidence shows that both the feedback strength and the fractional order can be chosen appropriately to adjust the synchronization effect effectively. The larger the fractional order and the feedback strength are, the faster the outer synchronization will be achieved. The method can be extended to the condition that the fractional orders of chaotic system in drive networks and response networks are nonidentical or chaotic systems of complex network are with incommensurate orders.



---

*Source: 580378-2014-04-23.xml*