Open AccessArticle

Adaptive Synchronization of Fractional-Order Complex Chaotic system with Unknown Complex Parameters

Ruoxun Zhang

^1,2,*,

Yongli Liu

¹ and

Shiping Yang

College of Primary Education, Xingtai University, Xingtai 054001, China

College of Physics Science and Information Engineering, Hebei Normal University, Shijiazhuang 050016, China

Author to whom correspondence should be addressed.

Entropy 2019, 21(2), 207; https://doi.org/10.3390/e21020207

Submission received: 27 December 2018 / Revised: 11 February 2019 / Accepted: 19 February 2019 / Published: 21 February 2019

(This article belongs to the Special Issue Nonlinear Dynamics and Entropy of Complex Systems with Hidden and Self-excited Attractors)

Download

Browse Figures

Versions Notes

Abstract

This paper investigates the problem of synchronization of fractional-order complex-variable chaotic systems (FOCCS) with unknown complex parameters. Based on the complex-variable inequality and stability theory for fractional-order complex-valued system, a new scheme is presented for adaptive synchronization of FOCCS with unknown complex parameters. The proposed scheme not only provides a new method to analyze fractional-order complex-valued system but also significantly reduces the complexity of computation and analysis. Theoretical proof and simulation results substantiate the effectiveness of the presented synchronization scheme.

Keywords:

synchronization; fractional-order; complex-variable chaotic system; unknown complex parameters

1. Introduction

In the past 20 years, fractional-order chaotic systems have been extensively studied due to their wide applications in the fields of secure communication, control engineering, finance, physical and mathematical science, entropy, encryption and signal processing [1,2,3,4]. Meanwhile, synchronization of such systems has aroused tremendous attention of many researchers. Lots of excellent results were obtained and some methods of synchronization have been presented [5,6,7,8,9,10,11,12,13,14,15,16,17]. In various synchronization methods, the adaptive control approach is an effective method to realize the synchronization of uncertain systems.

The aforementioned works mainly investigated the fractional-order systems with real variables, not involving complex variables. Because complex variables that double the number of variables can generate complicated dynamical behaviors, enhance anti-attack ability and achieve higher transmission efficiency [18,19,20], many researchers have taken complex variables into the fractional-order systems and investigated dynamics behavior, stability, stabilization and synchronization of FOCCS in recent years. In [21,22,23], fractional-order complex-variable Chen system, T system and Lorenz system have been investigated, respectively. Recently, Zhang et al. [24] have investigated synchronization of fractional-order complex-valued delayed neural networks. Li et al. [25] presented adaptive synchronization scheme for fractional-order complex-valued neural networks with discrete and distributed delays. Sun et al. [26] proposed real combination synchronization of three fractional-order complex-variable chaotic systems, Yadav et al. [27] studied Dual function projective synchronization of fractional order complex chaotic systems, Nian et al. [28] realized synchronization of fractional-order complex chaotic system with parametric and external disturbances via sliding mode control method and Jiang et al. have studied complex modified projective synchronization (CMPS) for FOCCS in [29]. However, in these papers [24,25,26,27,28,29], the parameters of the FOCCS are exactly known in priori. In fact, in many practical engineering situations, most of system parameters cannot be accurately determined in advance and chaos synchronization will be destroyed with these uncertain factors. Hence, it is an important problem to realize synchronization of FOCCS with unknown complex parameters.

Inspired by the above discussions, the synchronization problem of FOCCS with unknown complex parameters was investigated in this paper. Using the inequality of the fractional derivative containing complex variable and the stability theory for fractional-order complex-valued system, we realized synchronization of such systems by constructing a suitable response system. It should be noted that we deal with the synchronization problem of fractional-order uncertain complex-variable system in complex-valued domain. That is to say, it is not necessary to separate the complex-variable system into its real and imaginary parts. This greatly reduces the complexity of computation and the difficulty of theoretical analysis.

Notation:

ℂ^{n}

denotes complex n-dimensional space. For

z \in ℂ^{n}

z^{r}

z^{i}

\bar{z}

z^{T}

z^{H}

and

| | z | |

are the real part, imaginary part, conjugate, transpose, conjugate transpose and l₂-norm of

z

, respectively. For a matrix

A \in ℂ ​^{n \times n}

A^{H}

denotes its conjugate transpose.

2. Preliminaries

Definition 1

[30]. The fractional integral of order α for a function f is defined as:

I^{α} f (t) = {}_{t_{0}}D_{t}^{- α} f (t) = \frac{1}{Γ (α)} \int_{t_{0}}^{t} {(t - τ)}^{α - 1} f (τ) d τ

(1)

where t ≥ t₀ and α > 0.

Definition 2

[30]. Caputo’s fractional derivative of order

α

for a function

f \in ℝ^{n}

is defined by:

{}_{t_{0}}^{C}D_{t}^{α} f (t) = \frac{1}{Γ (n - α)} \int_{t_{0}}^{t} \frac{f^{(n)} (τ)}{{(t - τ)}^{α - n + 1}} d τ

(2)

where t ≥ t₀ and n is a positive integer such that n − < α < 1.

Lemma 1

[31]. Let

x (t) \in ℝ^{n}

be a continuous and derivable vector function. Then, for any time instant t ≥ t₀ and

α \in (0, 1)

\frac{1}{2} {}_{t_{0}}^{C}D_{t}^{α} [x^{T} (t) x (t)] \leq x^{T} (t) {}_{t_{0}}^{C}D_{t}^{α} x (t)

(3)

Corollary 1.

For a scalar derivable function

φ (t)

and a constants C, we have:

\frac{1}{2} {}_{t_{0}}^{C}D_{t}^{α} {(φ (t) - C)}^{2} \leq (φ (t) - C) {}_{t_{0}}^{C}D_{t}^{α} φ (t)

(4)

Lemma 2

[32] . Let

z \in ℂ^{n}

be a differentiable complex-valued vector. Then,

\forall t \geq t_{0}

and

α \in (0, 1]

, the following inequality holds:

{}_{t_{0}}^{C}D_{t}^{α} z^{H} (t) P z (t) \leq z^{H} (t) P {}_{t_{0}}^{C}D_{t}^{α} z (t) + {({}_{t_{0}}^{C}D_{t}^{α} z (t))}^{H} P z (t)

(5)

where

P \in ℂ^{n \times n}

is a constant positive definite Hermitian matrix.

Lemma 3.

For the fractional-order complex-variable systems:

{}_{0}^{C}D_{t}^{α} z (t) = h (z (t))

(6)

where

0 < α < 1

z (t) = {(z_{1}, z_{2}, \dots, z_{n})}^{T} \in ℂ^{n}

is the system complex state vector,

h \in ℂ^{n}

is a continuous nonlinear function complex vector, which satisfies the globally Lipschitz continuity condition in the complex domain. Let

z (t) = 0

be an equilibrium point of system (1) and let

V_{1} (t) = z^{H} (t) z (t)

and

V_{2} (z (t)) \geq 0

are continuously differentiable functions. If:

V (t) = V_{1} (t) + V_{2} (z (t))

(7)

and:

{}_{0}^{t}D_{t}^{α} V (t) \leq - θ V_{1} (t)

(8)

where

θ

is a positive constant. Then z(t) = 0 is asymptotically stable.

Proof.

See the Appendix A. It was pointed out [33,34,35] that a similar theorem is obtained for the real systems.

Remark 1.

Using Lemmas 2.2–2.3, one can directly analyze fractional order complex-variable system in the complex space.

3. Main Results

We considered a kind of FOCCS described by:

{}_{0}^{C}D_{t}^{α} z (t) = A z (t) + f (z (t))

(9)

where

0 < α < 1,

z (t) = {(z_{1}, z_{2}, \dots, z_{n})}^{T} \in ℂ^{n}

is the system complex state vector,

f \in ℂ^{n}

is a continuous nonlinear function vector, which satisfies the globally Lipschitz continuity condition in the complex domain and

A \in ℂ^{n \times n}

is unknown (complex or real) parameter matrix. Furthermore, Equation (10) can be rewritten as:

{}_{0}^{C}D_{t}^{α} z (t) = g (z (t)) θ + f (z (t))

(10)

where

g : ℂ^{n} \to ℂ^{n \times m}

is a complex function matrix, and

θ = {(θ_{1}, θ_{2}, \dots, θ_{m})}^{T}

is the system unknown complex parameter vector. For system (9) or (10), there are two propositions as follows.

Proposition 1.

There exists a positive constant l₁ such that the following inequality holds:

{(z - w)}^{H} [g (z) - g (w)] θ + {[g (z) - g (w)] θ}^{H} (z - w) \leq l_{1} {(z - w)}^{H} (z - w)

(11)

Proof.

Given that

g (z (t)) θ = A z (t)

results in:

\begin{array}{l} {(z - w)}^{H} [g (z) - g (w)] θ + {[g (z) - g (w)] θ}^{H} (z - w) \\ = {(z - w)}^{H} A (z - w) + {(z - w)}^{H} A^{H} (z - w) \\ = {(z - w)}^{H} (A + A^{H}) (z - w) \end{array}

Since

A + A^{H}

is Hermitian Matrix:

λ_{m} {(z - w)}^{H} (z - w) \leq {(z - w)}^{H} (A + A^{H}) (z - w) \leq λ_{M} {(z - w)}^{H} (z - w)

where

λ_{m}

and

λ_{M}

are the minimum and maximum eigenvalue of

A + A^{H}

, respectively [36,37].

Let

l_{1} = \max (| λ_{m} |, | λ_{M} |)

, then, one has:

{(z - w)}^{H} [g (z) - g (w)] θ + {[g (z) - g (w)] θ}^{H} (z - w) = {(z - w)}^{H} (A + A^{H}) (z - w) \leq l_{1} {(z - w)}^{H} (z - w)

Proposition 2

[38] . For the Lipschitz continuous function

f \in ℂ^{n}

, there exists a positive constant l₂ such that the following inequality holds:

{(z - w)}^{H} [f (z) - f (w)] + {[f (z) - f (w)]}^{H} (z - w) \leq l_{2} {(z - w)}^{H} (z - w)

(12)

Proof.

For

f \in ℂ^{n}

, Lipschitz is continuous, then

| | f (z) - f (w) | | \leq L | | z - w | |

, where

L \geq 0

is a constant. It follows:

\begin{array}{l} {(z - w)}^{H} [f (z) - f (w)] + {[f (z) - f (w)]}^{H} (z - w) = 2 Re {{(z - w)}^{H} [f (z) - f (w)]} \\ \leq 2 | z - w |^{T} | f (z) - f (w) | \leq | z - w |^{T} | z - w | + | f (z) - f (w) |^{T} | f (z) - f (w) | \\ = {(z - w)}^{H} (z - w) + | | f (z) - f (w) | |^{2} \leq {(z - w)}^{H} (z - w) + L^{2} | | z - w | |^{2} \\ = (1 + L^{2}) {(z - w)}^{H} (z - w) = l_{2} {(z - w)}^{H} (z - w) \end{array}

where

l_{2} = L^{2} + 1

| f (z) - f (w) | = {(| f_{1} (z) - f_{1} (w) |, | f_{2} (z) - f_{2} (w) |, \dots, | f_{n} (z) - f_{n} (w) |)}^{T}

and

| z - w | = {(| z_{1} - w_{1} |, | z_{2} - w_{2} |, \dots, | z_{n} - w_{n} |)}^{T}

Remark 2.

It is easy to check that many typical FOCCSs, such as the fractional-order complex-variable Chen system [21], T system [22] and Lorenz system [23] all satisfy Propositions 1 and 2.

Choose system (11) as the master system, then the controlled response system is given by:

{}_{0}^{C}D_{t}^{α} w (t) = g (w (t)) \hat{θ} + f (w (t)) + u (t)

(13)

where

w (t) = {(w_{1}, w_{2}, \dots, w_{n})}^{T}

is the complex state vector,

\hat{θ} \in ℂ^{m}

represents the estimate vector of unknown vector

θ

, and

u (t) = {(u_{1} (t), u_{2} (t), \dots, u_{n} (t))}^{T}

is controller to be determined.

Theorem 1.

Asymptotically synchronization and parameter identification of systems (13) and (11) can be achieved under adaptive controller:

u (t) = - k e (t)

(14)

and the complex update laws:

{}_{0}^{C}D_{t}^{α} k = σ e^{H} (t) e (t)

(15)

{}_{0}^{C}D_{t}^{α} e_{θ} = {}_{0}^{C}D_{t}^{α} \hat{θ} = - η g^{H} (w (t)) e (t)

(16)

where

e (t) = w (t) - z (t)

is the error vector,

e_{θ} = \hat{θ} - θ

is the parameter error, σ, η are two arbitrary positive constants.

Proof.

From the error vector and systems (11) and (13), it yields:

\begin{array}{l} {}_{0}^{C}D_{t}^{α} e (t) & = g (w (t)) \hat{θ} + f (w (t)) - g (z (t)) θ + f (z (t)) + u (t) \\ = g (w (t)) e_{θ} + [g (w (t)) - g (z (t))] θ + f (w (t)) - f (z (t)) + u (t) \end{array}

Let us present the Lyapunov function:

V (t, e (t)) = e^{H} (t) e (t) + \frac{1}{σ} {(k - k^{*})}^{2} + \frac{1}{η} e_{θ}^{H} (t) e_{θ} (t)

(17)

where

k^{*}

is to be determined.

Using Lemma 2.1, Corollary 2.1 and Lemma 2.2, we have:

\begin{array}{l} {}_{0}^{C}D_{t}^{α} V (t, e (t)) = {}_{0}^{C}D_{t}^{α} [e^{H} (t) e (t) + \frac{1}{σ} {(k - k^{*})}^{2} + \frac{1}{η} e_{θ}^{H} (t) e_{θ} (t)] \\ \leq e^{H} (t) {}_{0}^{C}D_{t}^{α} e (t) + {[{}_{0}^{C}D_{t}^{α} e (t)]}^{H} e (t) + \frac{2}{σ} (k - k^{*}) {}_{0}^{C}D_{t}^{α} k + \frac{1}{η} e_{θ}^{H} (t) {}_{0}^{C}D_{t}^{α} e_{θ} (t) + \frac{1}{η} {[{}_{0}^{C}D_{t}^{α} e_{θ} (t)]}^{H} e_{θ} (t) \\ \leq e^{H} (t) {g (w (t)) e_{θ} + [g (w (t)) - g (z (t))] θ + f (w (t)) - f (z (t)) - k e (t)} + {g (w (t)) e_{θ} \\ + [g (w (t)) - g (z (t))] θ + f (w (t)) - f (z (t)) - k e (t)}^{H} e + \frac{2}{σ} (k - k^{*}) {}_{0}^{C}D_{t}^{α} k \\ + \frac{1}{η} e_{θ}^{H} (t) {}_{0}^{C}D_{t}^{α} e_{θ} (t) + \frac{1}{η} {[{}_{0}^{C}D_{t}^{α} e_{θ} (t)]}^{H} e_{θ} (t) \\ \leq e^{H} (t) {g (w (t)) e_{θ} + [g (w (t)) - g (z (t))] θ + f (w (t)) - f (z (t)) - k e (t)} + {g (w (t)) e_{θ} \\ + [g (w (t)) - g (z (t))] θ + f (w (t)) - f (z (t)) - k e (t)}^{H} e + \frac{2}{σ} (k - k^{*}) {}_{0}^{C}D_{t}^{α} k \\ + \frac{1}{η} e_{θ}^{H} (t) {}_{0}^{C}D_{t}^{α} e_{θ} (t) + \frac{1}{η} {[{}_{0}^{C}D_{t}^{α} e_{θ} (t)]}^{H} e_{θ} (t) \\ \leq e^{H} (t) {g (w (t)) e_{θ} + [g (w (t)) - g (z (t))] θ + f (w (t)) - f (z (t)) - k e (t)} + {g (w (t)) e_{θ} \\ + [g (w (t)) - g (z (t))] θ + f (w (t)) - f (z (t)) - k e (t)}^{H} e + \frac{2}{σ} (k - k^{*}) {}_{0}^{C}D_{t}^{α} k \\ + \frac{1}{η} e_{θ}^{H} (t) {}_{0}^{C}D_{t}^{α} e_{θ} (t) + \frac{1}{η} {[{}_{0}^{C}D_{t}^{α} e_{θ} (t)]}^{H} e_{θ} (t) \end{array}

Substitute Equations (15) and (16) into the inequality above, we further have:

\begin{array}{l} {}_{0}^{C}D_{t}^{α} V (t, e (t)) \leq e^{H} (t) {g (w (t)) e_{θ} + [g (w (t)) - g (z (t))] θ + f (w (t)) - f (z (t)) - k e (t)} \\ + {g (w (t)) e_{θ} + [g (w (t)) - g (z (t))] θ + f (w (t)) - f (z (t)) - k e (t)}^{H} e \\ + 2 (k - k^{*}) e^{H} (t) e (t) - e_{θ}^{H} (t) g^{H} (w (t)) e (t) - {[g^{H} (w (t)) e (t)]}^{H} e_{θ} (t) \\ \leq e^{H} (t) [g (w (t) - g (z (t))] θ + {{[g (w (t) - g (z (t))] θ}}^{H} e (t) \\ + e^{H} (t) [f (w (t) - f (z (t))] + [f {(w (t) - f (z (t))]}^{H} e (t) - 2 k^{*} e^{H} (t) e (t) \end{array}

From Propositions 1 and 2, we can obtain:

e^{H} (t) [g (w (t)) - g (z (t)] θ + {[g (w (t)) {- g (z (t)] θ}}^{H} e (t) \leq l_{1} e^{H} (t) e (t)

and:

e^{H} (t) [f (w (t)) - f (z (t)] + [f (w (t)) - f {(z (t)]}^{H} e (t) \leq l_{2} e^{H} (t) e (t)

then, one has:

{}_{0}^{C}D_{t}^{α} V (t, e (t)) \leq e^{H} (t) [l_{1} + (l_{2} - 2 k^{*}) I] e (t)

Let

2 k^{*} = l_{1} + l_{2} + 1

, then:

{}_{0}^{C}D_{t}^{α} V (t, e (t)) \leq - e^{H} (t) e (t)

(18)

According to Lemma 2.3, one has

\lim_{t \to \infty} e^{H} (t) e (t) = 0,

which implies

\lim_{t \to \infty} e (t) = 0

, which shows that the systems (11) and (13) can obtain asymptotically synchronization. Meanwhile, according to Remark 1 of Theorem 1 in [39], parameter identification is achieved.

Remark 3.

In previous work [24,25,26,27,28,29], the common method to analyze fractional complex-valued systems is to separate into two real-valued systems according to their real and imaginary parts, and then the criteria on synchronization were obtained by investigating these real-valued systems. However, there are two problems with this approach. One is that the dimension of the real-valued system is twice that of the original complex-valued system, which adds the complexity of computation and analysis. The other is that this method requires that complex-valued functions be explicitly separated into real and imaginary parts. However, this separation is not always expressible in an analytical form. Unlike from previous works, in our proposed method, the entire analysis process is performed in the complex-valued domain, and the complex function theory is used to derive synchronization conditions without separating the original complex-valued chaotic system into two real-valued systems, which reduces the complexity of analysis and computation. Moreover, the proposed method can be applied to other complex-valued systems, such as complex networks with fractional-order complex-variable dynamics and fractional-order complex-valued neural network systems.

Remark 4.

If the system parameters are known, the update law will be reduced to (15) only.

4. Numerical Simulations

In this section, in order to show the effectiveness of the proposed scheme in preceding section, numerical example on fractional-order complex chaotic system will be provided. When numerically solving such systems, we first adopt the predictor–corrector method [40] by MATLAB. Lyapunov exponents of the systems are calculated by adopting the Wolf et al. algorithm [41] with some changes.

Consider the Lorenz-like fractional-order complex chaotic system with commensurate order:

(\begin{matrix} {}_{0}^{C}D_{t}^{α} z_{1} \\ {}_{0}^{C}D_{t}^{α} z_{2} \\ {}_{0}^{C}D_{t}^{α} z_{3} \end{matrix}) = (\begin{matrix} a (z_{2} + z_{1}) \\ - c z_{2} - z_{1} z_{3} \\ {\bar{z}}_{1} z_{1} - b z_{3} \end{matrix}) = (\begin{matrix} z_{2} + z_{1} & 0 & 0 \\ 0 & - z_{2} & 0 \\ 0 & 0 & - z_{3} \end{matrix}) (\begin{matrix} a \\ c \\ b \end{matrix}) + (\begin{matrix} 0 \\ - z_{1} z_{3} \\ {\bar{z}}_{1} z_{1} \end{matrix})

(19)

where

z_{1}, z_{2}, z_{3}

are the complex state variables and a, b, c are system parameters; let

a = 10 + i,

b = 3, c = 16 + 0.3 i

. The maximum Lyapunov exponent (MLE) spectrum is depicted in Figure 1a, and the bifurcation diagram is presented in Figure 1b. Figure 1a,b shows that system (19) is chaotic with fractional order

α \in [0.91, 0.98] \cup [0.985, 1]

and parameters

a = 10 + i, b = 3, c = 16 + 0.3 i

. When the fractional-order

α = 0.95

, the attractor trajectories are illustrated in Figure 2.

Recently, ref. [42] described how to perform a successful simulation and optimization, and how to synthesize the mathematical models using CMOS technology. The application of metaheuristics to optimize MLE by varying the parameters of the oscillators was discussed.

Field-programmable gate array (FPGA)-based implementation of chaotic oscillators has demonstrated its usefulness in the development of engineering applications in a wide variety of fields, such as: random number generators, robotics and chaotic secure communication systems, signal processing. Very recently, Pano-Azucena et al. [43] implemented the chaotic system using a field-programmable gate array (FPGA) based on trigonometric polynomials. Reference [44] detailed the FPGA-based implementation of all the fractional order chaotic oscillators applying Grünwald-Letnikov(G-L) method. Their work proved experimentally that applying G-L method with 256 elements of memory; it can observe different families of Fractional-order chaotic attractors having working frequencies between 77.59 MHz and 84.9 MHz. This is very beneficial to the development of fractional-order chaos in engineering applications. For the FPGA-based implementation of FOCCS, the FOCCS was first separated into two real-variable systems according to their real and imaginary parts, and then these real-variable systems can be implemented using FPGA by the method proposed by. In order to find much better behavior and characteristics of the FOCCS in the complex domain, we used the G-L method to numerically solve the system (19) again. The MLE spectrum with varying parameter aⁱ (the imaginary of a) is depicted in Figure 3a, the bifurcation diagram is presented in Figure 3b, the state trajectories are illustrated in Figure 4.

Taking the system (19) as master system, and assuming the parameters a, b and c are unknown, the response system is given as follows:

(\begin{matrix} {}_{0}^{C}D_{t}^{α} w_{1} \\ {}_{0}^{C}D_{t}^{α} w_{2} \\ {}_{0}^{C}D_{t}^{α} w_{3} \end{matrix}) = (\begin{matrix} w_{2} + w_{1} & 0 & 0 \\ 0 & - w_{2} & 0 \\ 0 & 0 & - w_{3} \end{matrix}) (\begin{matrix} \hat{a} \\ \hat{c} \\ \hat{b} \end{matrix}) + (\begin{matrix} 0 \\ - w_{1} w_{3} \\ {\bar{w}}_{1} w_{1} \end{matrix}) + (\begin{matrix} u_{1} \\ u_{2} \\ u_{3} \end{matrix})

(20)

where

\hat{a}, \hat{b}, \hat{c}

are parameter estimations.

u_{1}, u_{2}, u_{3}

are the controller.

According to Theorem 1, the controllers and the update rules are selected as:

u_{1} = - k e_{1}, u_{2} = - k e_{2}, u_{3} = - k e_{3} {}_{0}^{C}D_{t}^{α} k = σ e^{H} e = σ ({\bar{e}}_{1} e_{1} + {\bar{e}}_{2} e_{2} + {\bar{e}}_{3} e_{3}), (σ > 0)

(21)

(\begin{matrix} {}_{0}^{C}D_{t}^{α} \hat{a} \\ {}_{0}^{C}D_{t}^{α} \hat{c} \\ {}_{0}^{C}D_{t}^{α} \hat{b} \end{matrix}) = {(\begin{matrix} w_{1} + w_{2} & 0 & 0 \\ 0 & - w_{2} & 0 \\ 0 & 0 & - w_{3} \end{matrix})}^{H} (\begin{matrix} e_{1} \\ e_{2} \\ e_{3} \end{matrix}) = (\begin{matrix} - (w_{1} + w_{2}) e_{1} \\ {\bar{w}}_{2} e_{2} \\ {\bar{w}}_{3} e_{3} \end{matrix})

(22)

In the simulation, let

α = 0.95

, (a, b, c) = (30 + i, 3, 26 + 0.6i), the initial conditions

z (0) =

{(1 + i, - 2 - i, 6)}^{T}, w (0) = {(- 1 + i, - 3 + i, 10)}^{T},

(\hat{a} (0), \hat{b} (0), \hat{c} (0)) = (20, 2, 20)

k (0) = 0

and

σ = 6

. Two systems can achieve synchronization and the parameters are identified, as shown in Figure 5 and Figure 6.

5. Conclusions

We studied the adaptive synchronization of FOCCS with unknown complex parameters, and proposed a method for analyzing FOCCS without separating system into real and imaginary parts. By this method, the constructed response system can be asymptotically synchronized to an uncertain drive system with a desired complex scaling diagonal matrix. The proposed synchronization scheme retains the complex nature of fractional-order complex chaotic system. It not only provides a new method of analyzing FOCCS, but also significantly decreases the complexity of computation and analysis. We hope that the work performed will be helpful to further research of nonlinear fractional order complex-variable systems.

Author Contributions

Conceptualization, R.Z.; software, Y.L.; supervision, S.Y.

Funding

The author R.Z was supported by the Natural Science Foundation of Hebei province, China (No.A2015108010), and science and technology support program of Xingtai, China (No. 2016ZC191).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Lemma 3.

α

-integrating (8), we have:

V (t) - V (0) \leq - θ I^{α} V_{1} (t) = - θ I^{α} z^{H} (t) z (t) \leq 0

(A1)

Thus,

V (t) \leq V (0), t \geq 0

. From (6), we can obtain that

V_{1} (t) = z^{H} (t) z (t)

is bounded. Furthermore, since h(z(t)) satisfies the globally Lipschitz continuity condition, from (6), one has,

| | {}_{0}^{C}D_{t}^{α} z (t) | |

= | | h (z (t)) | |

\leq l | | z^{H} (t) z (t)) | |

, i.e.,

{[\sum_{i = 1}^{n} (\bar{{}_{0}^{C}D_{t}^{α} z_{i} (t)} {}_{0}^{C}D_{t}^{α} z_{i} (t))]^{1 / 2} \leq l [\sum_{i = 1}^{n} (\bar{z_{i} (t)} z_{i} (t))]}^{1 / 2}

, where l is positive constant. Given that

V_{1} (t) = z^{H} (t) z (t)

is bounded, we have

| z_{i} |

is bounded, and then

| {}_{0}^{C}D_{t}^{α} z_{i} (t) |

is bounded. Thus, there exists a constant M > 0, such that

| {}_{0}^{C}D_{t}^{α} z_{i} (t) | \leq M

. For

0 \leq t_{1} < t_{2}

and any

ε > 0,

| t_{2} - t_{1} | < δ (ε)

= {[\frac{ε ​ Γ (α + 1)}{2 M}]}^{\frac{1}{α}}

, one can get:

\begin{array}{l} | z_{i} (t_{1}) - z_{i} (t_{2}) | = | {}_{0}D_{t}^{- α} {}_{0}^{C}D_{t}^{α} z_{i} (t_{1}) + z_{i} (0) - [{}_{0}D_{t}^{- α} {}_{0}^{C}D_{t}^{α} z_{i} (t_{2}) + z_{i} (0)] | \\ = | {}_{0}D_{t}^{- α} {}_{0}^{C}D_{t}^{α} z_{i} (t_{1}) - {}_{0}D_{t}^{- α} {}_{0}^{C}D_{t}^{α} z_{i} (t_{2}) | \\ = \frac{1}{Γ (α)} | \int_{0}^{t_{1}} [{(t_{1} - τ)}^{α - 1} - {(t_{2} - τ)}^{α - 1}] {}_{0}^{C}D_{τ}^{α} z_{i} (τ) d τ - \int_{t_{1}}^{t_{2}} {(t_{2} - τ)}^{α - 1} {}_{0}^{C}D_{τ}^{α} z_{i} (τ) d τ | \\ \leq \frac{M}{Γ (α)} {\int_{0}^{t_{1}} [{(t_{1} - τ)}^{α - 1} - {(t_{2} - τ)}^{α - 1}] d τ + \int_{t_{1}}^{t_{2}} {(t_{2} - τ)}^{α - 1} d τ} \\ = \frac{M}{Γ (α + 1)} [(t_{1}^{α} - t_{2}^{α}) + 2 {(t_{2} - t_{1})}^{α}] \\ \leq \frac{2 M}{Γ (α + 1)} {(t_{2} - t_{1})}^{α} \leq ε \end{array}

Hence,

z_{i} (t)

is uniformly continuous. Therefore,

V_{1} (t) = z^{H} (t) z (t)

is uniformly continuous.

Next, we adopt contradiction to prove

\lim_{t \to \infty} V_{1} (t) = \lim_{t \to \infty} z^{H} (t) z (t) = 0

using the idea of Theorem 1 in Reference [45].

Suppose that

V_{1} (t) \neq 0

t \to \infty

. Then there exists a monotone increasing sequence

{(t_{k})}_{k \in N^{+}}

(

t_{k} \to \infty

k \to \infty

) and a positive constant

ε

> 0 such that

V_{1} (t_{k}) > ε

. Since the uniform continuity of

V_{1} (t)

, for the given

ε

\exists δ > 0 (δ \leq \inf_{j \in N^{+}} {t_{j + 1} - t_{j}},

which implies that the intervals

[t_{k}, t_{k} + δ]

are nonoverlapping) such that

| V_{1} (t_{k}) - V_{1} (t) | \leq ε / 2

. Then, for any

t \in [t_{k}, t_{k} + δ]

it follows that

V_{1} (t) = V_{1} (t_{k}) - V_{1} (t_{k}) + V_{1} (t)

\geq | V_{1} (t_{k}) | - | V_{1} (t_{k}) - V_{1} (t) | > ε / 2

. Thus, for any k = 1, 2, 3, …, from (A1), we have:

\begin{matrix} V (t_{k} + δ) - V (0) & \leq \frac{- θ}{Γ (α)} \int_{0}^{t_{k} + δ} \frac{V_{1} (τ)}{{(t_{k} + δ - τ)}^{1 - α}} d τ \\ = \frac{- θ}{Γ (α)} [\int_{0}^{t_{1}} + \int_{t_{1}}^{t_{1} + δ} + \int_{t_{1} + δ}^{t_{2}} + \int_{t_{2}}^{t_{2} + δ} + \int_{t_{2} + δ}^{t_{3}} + \int_{t_{3}}^{t_{3} + δ} + \dots + \int_{t_{k}}^{t_{k} + δ}] \frac{V_{1} (τ)}{{(t_{k} + δ - τ)}^{1 - α}} d τ \\ \leq \frac{- θ}{Γ (α)} [\int_{t_{1}}^{t_{1} + δ} + \int_{t_{2}}^{t_{2} + δ} + \int_{t_{3}}^{t_{3} + δ} + \dots + \int_{t_{k}}^{t_{k} + δ}] \frac{V_{1} (τ)}{{(t_{k} + δ - τ)}^{1 - α}} d τ \\ \leq \frac{- θ ε}{2 Γ (α)} [\int_{t_{1}}^{t_{1} + δ} + \int_{t_{2}}^{t_{2} + δ} + \int_{t_{3}}^{t_{3} + δ} + \dots + \int_{t_{k}}^{t_{k} + δ}] {(t_{k} + δ - τ)}^{α - 1} d τ \\ = \frac{- θ ε}{2 Γ (α)} \sum_{j = 1}^{k} \int_{t_{j}}^{t_{j} + δ} {(t_{k} + δ - τ)}^{α - 1} d τ \end{matrix}

Given that

{(t_{k} + δ - τ)}^{α - 1} \geq {(t_{k} + δ - t_{j})}^{α - 1}

for all

τ \in [t_{j}, t_{j} + δ]

results in:

\begin{matrix} V (t_{k} + δ) - V (0) & \leq \frac{- θ ε}{2 Γ (α)} \sum_{j = 1}^{k} \int_{t_{j}}^{t_{j} + δ} {(t_{k} + δ - τ)}^{α - 1} d τ \\ \leq \frac{- θ ε}{2 Γ (α)} \sum_{j = 1}^{k} \int_{t_{j}}^{t_{j} + δ} {(t_{k} + δ - t_{j})}^{α - 1} d τ \\ = \frac{- θ ε δ}{2 Γ (α)} \sum_{j = 1}^{k} {(t_{k} + δ - t_{j})}^{α - 1} \\ \leq \frac{- θ ε δ}{2 Γ (α)} \sum_{j = 1}^{k} {(t_{k} + δ - t_{1})}^{α - 1} \\ = \frac{- θ ε δ}{2 Γ (α)} \frac{k}{{(t_{k} + δ - t_{1})}^{1 - α}} \\ \leq \frac{- θ ε δ}{2 Γ (α)} \frac{k}{{(k d)}^{1 - α}} = \frac{- θ ε δ}{2 Γ (α)} \frac{k^{α}}{{(d)}^{1 - α}} \end{matrix}

where

d = \sup_{j \in N, 2 \leq j \leq k} {t_{j} - t_{j - 1}}

(since

V_{1} (t)

is a uniformly continuous function and assumed

V_{1} (t) \neq 0

, as

t \to \infty

, d is bounded). Obviously,

V (t_{k} + δ) - V (0) \leq \frac{- θ ε δ}{2 Γ (α)} \frac{k^{α}}{{(d)}^{1 - α}} \to - \infty

k \to \infty

, which contradict with

V (t) \geq 0

. Therefore,

\lim_{t \to \infty} V_{1} (t) = \lim_{t \to \infty} z^{H} (t) z (t) = 0

. i.e., z(t) = 0 is asymptotically stable. This completes the proof for Lemma 3.

References

Gorenlo, R.; Mainardi, F. Fractional Calculus; Springer: Berlin, Germany, 1997. [Google Scholar]
Caponetto, R. Fractional Order Systems: Modeling and Control Applications; World Scientiic: Singapore, 2010. [Google Scholar]
Magin, R.L. Fractional Calculus in Bioengineering; Begell House Redding: New York, NY, USA, 2006. [Google Scholar]
Radwan, A.G.; Shamim, A.; Salama, K.N. Theory of fractional order elements based impedance matching networks. IEEE Microwave Wirel. Comp. Lett. 2011, 21, 120–122. [Google Scholar] [CrossRef]
Radwan, A.G.; Abd-El-Haiz, S.K.; Abdelhaleem, S.H. Image encryption in the fractional-order domain. In Proceedings of the 1st International Conference on Engineering and Technology, Cairo, Egypt, 10–11 October 2012. [Google Scholar]
Kiani-B, A.; Fallah, K.; Pariz, N.; Leung, H. A chaotic secure communication scheme using fractional chaotic systems based on an extended fractional Kalman filter. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 863–879. [Google Scholar] [CrossRef]
Li, G.J.; Liu, H. Stability analysis and synchronization for a class of fractional-order neural networks. Entropy 2016, 18, 55. [Google Scholar] [CrossRef]
Yin, C.; Cheng, Y.; Zhong, S.M.; Bai, Z. Fractional-order switching type control law design for adaptive sliding mode technique of 3D fractional-order nonlinear systems. Complexity 2016, 21, 363–373. [Google Scholar] [CrossRef]
Liu, H.; Li, S.G.; Wang, H.X.; Huo, Y.H.; Luo, J.H. Adaptive synchronization for a class of uncertain fractional-order neural networks. Entropy 2015, 17, 7185–7200. [Google Scholar] [CrossRef]
Lu, J.G. Chaotic dynamics of the fractional-order Lü system and its synchronization. Phys. Lett. A 2006, 354, 305–311. [Google Scholar] [CrossRef]
Li, C.G.; Chen, G.R. Chaos in the fractional-order Chen system and it’s control. Chaos Solitons Fractals 2004, 22, 549–554. [Google Scholar] [CrossRef]
Agrawal, S.K.; Srivastava, M.; Das, S. Synchronization of fractional order chaotic systems using active control method. Chaos Solitons Fractals 2012, 45, 737–752. [Google Scholar] [CrossRef]
Chen, L.P.; Chai, Y.; Wu, R.C. Lag projective synchronization in fractional-order chaotic (hyperchaotic) systems. Phys. Lett. A 2011, 375, 2099–2110. [Google Scholar] [CrossRef]
Zhang, Y.; Sun, J. Chaotic synchronization and anti-synchronization based on suitable separation. Phys. Lett. A 2004, 330, 442–447. [Google Scholar] [CrossRef]
Zhang, R.X.; Yang, S.P. Robust synchronization of two different fractional-order chaotic systems with unknown parameters using adaptive sliding mode approach. Nonlinear Dyn. 2013, 71, 269–278. [Google Scholar] [CrossRef]
Zhang, R.X.; Yang, S.P. Adaptive synchronization of fractional-order chaotic systems via a single driving variable. Nonlinear Dyn. 2011, 66, 831–837. [Google Scholar] [CrossRef]
Zhang, R.X.; Yang, S.P. Robust chaos synchronization of fractional-order chaotic systems with unknown parameters and uncertain perturbations. Nonlinear Dyn. 2012, 69, 983–992. [Google Scholar] [CrossRef]
Wang, S.; Wang, X.; Han, B. Complex generalized synchronization and parameter identification of nonidentical nonlinear complex systems. PLoS ONE 2016, 11, e0152099. [Google Scholar] [CrossRef] [PubMed]
Jiang, C.; Zhang, F.; Li, T. Synchronization and anti-synchronization of N-coupled fractional-order complex chaotic systems with ring connection. Math. Meth. Appl. Sci. 2018, 41, 2625–2638. [Google Scholar] [CrossRef]
Zhang, F.F.; Liu, S. Self-time-delay synchronization of time-delay coupled complex chaotic system and its applications to communication. Int. J. Mod. Phys. C 2014, 25, 559–583. [Google Scholar] [CrossRef]
Luo, C.; Wang, X.Y. Chaos in the fractional-order complex Lorenz system and its synchronization. Nonlinear Dyn. 2013, 71, 241–257. [Google Scholar] [CrossRef]
Luo, C.; Wang, X. Chaos generated from the fractional-order complex Chen system and its application to digital secure communication. Int. J. Mod. Phys. C 2013, 24, 72–77. [Google Scholar] [CrossRef]
Liu, X.J.; Hong, L.; Yang, L.X. Fractional-order complex T system: Bifurcations, chaos control, and synchronization. Nonlinear Dyn. 2014, 75, 589–602. [Google Scholar] [CrossRef]
Zhang, W.W.; Cao, J.D.; Chen, D.Y.; Alsaadi, F.E. Synchronization in fractional-order complex-valued delayed neural networks. Entropy 2018, 20, 54. [Google Scholar] [CrossRef]
Li, L.; Wang, Z.; Lu, J.W.; Li, Y.X. Adaptive synchronization of fractional-order complex-valued neural networks with discrete and distributed delays. Entropy 2018, 20, 124. [Google Scholar] [CrossRef]
Sun, J.; Deng, W.; Cui, G.Z.; Wang, Y.F. Real combination synchronization of three fractional-order complex-variable chaotic systems. Optik 2016, 127, 11460–11468. [Google Scholar] [CrossRef]
Yadav, V.K.; Srikanth, N.; Das, S. Dual function projective synchronization of fractional order complex chaotic systems. Optik 2016, 127, 10527–10538. [Google Scholar] [CrossRef]
Nian, F.Z.; Liu, X.M.; Zhang, Y. Sliding mode synchronization of fractional-order complex chaotic system with parametric and external disturbances. Chaos Solitons Fractals 2018, 116, 22–28. [Google Scholar] [CrossRef]
Jiang, C.M.; Liu, S.T.; Zhang, F.F. Complex modified projective synchronization for fractional-order chaotic complex systems. Int. J. Auto Comput. 2017, 15, 603–615. [Google Scholar] [CrossRef]
Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999. [Google Scholar]
Aguila-Camacho, N.; Duarte-Mermoud, M.A.; Gallegos, J.A. Lyapunov functions for fractional order systems. Commun. Nonlinear Sci. Numer. Simulat. 2014, 19, 2951–2957. [Google Scholar] [CrossRef]
Quan, X.; Zhuang, S.; Liu, S.; Xiao, J. Decentralized adaptive coupling synchronization of fractional-order complex-variable dynamical networks. Neurocomputing 2016, 186, 119–126. [Google Scholar]
Zhang, R.X.; Liu, Y. A new Barbalat’s Lemma and Lyapunov stability theorem for fractional order systems. In Proceedings of the 29th Chinese Control and Decision Conference, Chongqing, China, 28–30 May 2017; pp. 3676–3681. [Google Scholar]
Wang, F.; Yang, Y.Q. Fractional order Barbalat’s Lemma and its applications in the stability of fractional order nonlinear systems. Math. Model. Anal. 2017, 22, 503–513. [Google Scholar] [CrossRef]
Wang, C. Adaptive fuzzy control for uncertain fractional-order financial chaotic systems subjected to input saturation. PLoS ONE 2016, 11, e0164791. [Google Scholar] [CrossRef]
Xu, Q.; Xu, X.H.; Zhuang, S.X.; Xiao, J.X.; Song, C.H. New complex projective synchronization strategies for drive-response networks with fractional complex-variable dynamics. Appl. Math. Comput. 2018, 338, 552–566. [Google Scholar] [CrossRef]
Tao, F.; Sun, J.T. Stability of complex-valued impulsive system with delay. Appl. Math. Comput. 2014, 240, 102–108. [Google Scholar]
Wu, Z.Y.; Chen, G.R.; Fu, X.C. Synchronization of a network coupled with complex-variable chaotic systems. Chaos 2012, 22, 023127. [Google Scholar] [CrossRef]
Gu, Y.; Yu, Y.; Wang, H. Synchronization-based parameter estimation of fractional-order neural networks. Physica A Stat. Mech. Appl. 2017, 483, 351–361. [Google Scholar] [CrossRef]
Diethelm, K.; Neville, J.F.; Freed, A.D. A predictor-corrector approach for the numerical solution of fractional differential equations. Nonlinear Dyn. 2002, 29, 1–4. [Google Scholar] [CrossRef]
Wolf, A.; Swift, J.B.; Swinney, H.L.; Vastano, J.A. Determining Lyapunov exponents from a time series. Physica D 1985, 16, 285–317. [Google Scholar] [CrossRef]
Carbajal-Gomez, V.H.; Tlelo-Cuautle, E.; Muñoz-Pacheco, J.M.; de la Fraga, L.G.; Sanchez-Lopez, C.; Fernandez-Fernandez, F.V. Optimization and CMOS design of chaotic oscillators robust to PVT variations: INVITED (In Press). Integration 2018. [Google Scholar] [CrossRef]
Pano-Azucena, A.D.; Tlelo-Cuautle, E.; Rodriguez-Gomez, G. de la Fraga, L.G. FPGA-based implementation of chaotic oscillators by applying the numerical method based on trigonometric polynomials. AIP Adv. 2018, 8, 075217. [Google Scholar] [CrossRef]
Pano-Azucena, A.D.; Ovilla-Martinez, B.; Tlelo-Cuautle, E.; Muñoz-Pacheco, J.M.; de la Fraga, L.G. FPGA-based implementation of different families of fractional-order chaotic oscillators applying Grounwald-Letnikov method. Commun. Nonlinear Sci. Numer. Simulat. 2019, 72, 516–527. [Google Scholar] [CrossRef]
Bao, H.; Park, J.H.; Cao, J. Adaptive synchronization of fractional-order memristor-based neural networks with time delay. Nonlinear Dyn. 2015, 82, 1343–1354. [Google Scholar] [CrossRef]

Figure 1. Dynamic behaviors of the fractional-order complex Lorenz-like System with commensurate order

α

(

a = 10 + i, b = 3, c = 16 + 0.3 i

). (a) maximal Lyapunov exponent; (b) bifurcation diagram.

Figure 1. Dynamic behaviors of the fractional-order complex Lorenz-like System with commensurate order

α

(

a = 10 + i, b = 3, c = 16 + 0.3 i

). (a) maximal Lyapunov exponent; (b) bifurcation diagram.

Figure 2. Chaotic attractors of fractional-order complex Lorenz-like system with

a = 10 + i, b = 3, c = 16 + 0.3 i

and commensurate order

α = 0.95

Figure 2. Chaotic attractors of fractional-order complex Lorenz-like system with

a = 10 + i, b = 3, c = 16 + 0.3 i

and commensurate order

α = 0.95

Figure 3. Dynamic behaviors of the fractional-order complex Lorenz-like System with commensurate order 0.95 (

a^{r} = 10, b = 3, c = 16 + 0.3 i

). (a) maximal Lyapunov exponent; (b) bifurcation diagram.

Figure 3. Dynamic behaviors of the fractional-order complex Lorenz-like System with commensurate order 0.95 (

a^{r} = 10, b = 3, c = 16 + 0.3 i

). (a) maximal Lyapunov exponent; (b) bifurcation diagram.

Figure 4. The state trajectories of fractional-order complex Lorenz-like system with

a = 10 + i, b = 3, c = 16 + 0.3 i

and commensurate order

α = 0.95

Figure 4. The state trajectories of fractional-order complex Lorenz-like system with

a = 10 + i, b = 3, c = 16 + 0.3 i

and commensurate order

α = 0.95

Figure 5. Synchronization errors e₁, e₂, e₃ of fractional-order complex Lorenz-like chaotic system.

Figure 6. Estimated complex parameters of fractional-order complex Lorenz-like chaotic system.

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Zhang, R.; Liu, Y.; Yang, S. Adaptive Synchronization of Fractional-Order Complex Chaotic system with Unknown Complex Parameters. Entropy 2019, 21, 207. https://doi.org/10.3390/e21020207

AMA Style

Zhang R, Liu Y, Yang S. Adaptive Synchronization of Fractional-Order Complex Chaotic system with Unknown Complex Parameters. Entropy. 2019; 21(2):207. https://doi.org/10.3390/e21020207

Chicago/Turabian Style

Zhang, Ruoxun, Yongli Liu, and Shiping Yang. 2019. "Adaptive Synchronization of Fractional-Order Complex Chaotic system with Unknown Complex Parameters" Entropy 21, no. 2: 207. https://doi.org/10.3390/e21020207

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

Article Menu

Adaptive Synchronization of Fractional-Order Complex Chaotic system with Unknown Complex Parameters

Abstract

1. Introduction

2. Preliminaries

3. Main Results

4. Numerical Simulations

5. Conclusions

Author Contributions

Funding

Conflicts of Interest

Appendix A

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI