Chaos Control and Synchronization of a Complex Rikitake Dynamo Model

Abstract

A novel chaotic system called complex Rikitake system is proposed. Dynamical properties, including symmetry, dissipation, stability of equilibria, Lyapunov exponents and bifurcation, are analyzed on the basis of theoretical analysis and numerical simulation. Further, based on feedback control method, the complex Rikitake system can be controlled to any equilibrium points. Additionally, this paper not only proves the existence of two types of synchronization schemes in the complex Rikitake system but also designs adaptive controllers to realize them. The proposed results are verified by numerical simulations.

1. Introduction

Since the pioneer research work of Ott et al. [1], Pecora and Carroll [2], the topic of chaos control and synchronization has attracted a lot of researchers in diverse areas including mathematics, physics, biology, medicine, engineering, and so on. Lots of research has been paid to study chaos control for real systems, and plenty of control methods have been put forward, such as feedback control [3,4], sliding mode control [5,6], backstepping method [7], and so on. These control strategies can also be employed to realize various kinds of synchronization of real chaos. Further developments in this direction can be found in [8,9,10,11,12,13,14].

The quoted literature above are only related to real chaotic systems and do not consider the chaotic systems which consist of complex variables. As is known to all, in the real world, many cases exist in the form of complex variables. For instance, Fowler et al. [15] discovered the complex Lorenz system when they studied laser physics and baroclinic instability of the geophysical flows in 1982. Since then, the study on complex nonlinear systems has been paid a substantial amount of attentions and has become a hot topic due to its wide applications in chemical systems, optics and especially in secure communications [16,17,18]. A considerable amount of complex dynamical systems exhibit chaotic motion, such as the complex Chen system [19], the time-delay complex Lorenz system [20], the complex generalised Lorenz hyperchaotic system [21], just to name a few examples. Compared with real chaos, complex chaos has the diversity of synchronization types and results. On the one hand, a lot of authors extend some synchronization schemes of real chaos into complex space, for example, complete synchronization (CS) [22], anti-synchronization (AS) [23], lag synchronization (LS) [24], combination synchronization [25], etc. On the other hand, some new synchronization schemes have been proposed on the basis of the characteristics of complex systems, such as complex complete synchronization (CCS) [26], complex lag synchronization (CLS) [27], complex anti lag synchronization (CALS) [28], combination complex synchronization [29,30], and so forth. However, the existing results on complex chaos have three disadvantages: Firstly, chaos control of the complex dynamical systems has gained little attention. Secondly, the existence of the synchronization problem, which is fundamental theoretical base, has not been considered so far. Finally, most of the current designed controllers eliminate the nonlinear term of the system, which are not only complicated but also difficult to realize in engineering. Therefore, control and synchronization in complex chaotic systems needs to be further and extensively studied.

Motivated by the aforementioned discussion, the current investigation concentrates on chaos control and synchronization of a novel complex dynamical system named as complex Rikitake system, which is proposed based on the Rikitake system. Following the idea of studying dynamics in chaotic systems, this paper investigates symmetry, dissipation, stability of equilibria, Lyapunov exponents, Poincaré-sections and bifurcation of the complex Rikitake system. Thus, along with the deeper understanding of feedback control method presented in [9], we construct simple adaptive controllers to realize control and synchronization of the complex Rikitake system. Furthermore, we obtain a criterion to detect the existence of synchronization in the complex Rikitake system and further prove that there exist CS and the coexistence of CS and AS.

The main construct of the article is arranged as follows. We present the complex Rikitake system and analyze some basic dynamics in Section 2. In Section 3, adaptive controllers are designed to control the complex Rikitake system to any equilibrium points. Section 4 gives the main results on chaos synchronization of the complex Rikitake system. The conclusions are provided in Section 5.

2. A Complex Chaotic Rikitake Dynamo System

In 1958, Rikitake discovered the 3-D Rikitake dynamo system [31] whose equations are

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}\dot{x}=-\beta x+yz,\hfill \\ \dot{y}=-\beta y-\alpha x+xz,\hfill \\ \dot{z}=1-xy,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(1)

where $x,y,z\in \mathbb{R}$ are state variables, $\alpha ,\beta >0$ are parameters. As mentioned in [32], the Rikitake system (1) behaves chaotically for $\alpha =5$ and $\beta =2$ with $(x_0,y_0,z_0)=(3,1,6)$, which are shown in Figure 1.

A new system can be generated by assuming that x and y are complex states and changing cross coupled terms x and y to complex conjugate form. Thus, we call it complex Rikitake system, which can be described as

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}\dot{x}=-\beta x+yz,\hfill \\ \dot{y}=-\beta y-\alpha x+xz,\hfill \\ \dot{z}=1-\frac{1}{2}(\overline{x}y+x\overline{y}),\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(2)

where $x=x_1+j{x}_2$, $y=x_3+j{x}_4$, $z=x_5$, $j=-1$, $\overline{x}$ and $\overline{y}$ denote the complex conjugates of x and y. Replacing $x,y$ in system (2) with real and imaginary variables can lead to the following equivalent system

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{x}}_1=-\beta x_1+x_3{x}_5,\hfill \\ {\dot{x}}_2=-\beta x_2+x_4{x}_5,\hfill \\ {\dot{x}}_3=-\beta x_3-\alpha x_1+x_1{x}_5,\hfill \\ {\dot{x}}_4=-\beta x_4-\alpha x_2+x_2{x}_5,\hfill \\ {\dot{x}}_5=1-x_1{x}_3-x_2{x}_4.\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(3)

In the next subsection, we study some dynamical properties of this new system (3).

2.1. Symmetry

Given a coordinate transformation T as follows

$$T(x_1,x_2,x_3,x_4,x_5)⟶(-x_1,-x_2,-x_3,-x_4,x_5).$$

It is clear that each trajectory is symmetrical with respect to the $x_5$-axis. That means system (3) is invariant for the given transformation T.

2.2. Dissipation

The divergence of system (3) can be calculated as

$$\begin{array}{c}\hfill \nabla V=\sum _{l=1}^5\frac{\partial {\dot{x}}_l}{\partial x_l}=-4\beta .\end{array}$$

As a result, it follows from the condition $\beta >0$ that system (3) is dissipative.

2.3. Equilibria and Stability

In order to find the equilibria of system (3), we consider equations in the form

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}-\beta x_1+x_3{x}_5=0,\hfill \\ -\beta x_2+x_4{x}_5=0,\hfill \\ -\beta x_3-\alpha x_1+x_1{x}_5=0,\hfill \\ -\beta x_4-\alpha x_2+x_2{x}_5=0,\hfill \\ 1-x_1{x}_3-x_2{x}_4=0.\hfill \end{array}\right.\hfill \end{array}\end{array}$$

After computation, we obtain the following equilibrium points:

$$S=(\frac{cos\theta }{r},\frac{sin\theta }{r},rcos\theta ,rsin\theta ,\frac{\alpha +\sqrt{{\alpha }^2+4{\beta }^2}}{2}),$$

where $r=\sqrt{\frac{2\beta }{\alpha +\sqrt{{\alpha }^2+4{\beta }^2}}}$ and $\theta \in [0,2\pi ]$. Now, we consider the stability of S. The Jacobian of system (3) at point S is deduced as:

$$\begin{array}{c}\hfill J_S=\left(\begin{array}{ccccc}-\beta & 0& \frac{\alpha +\sqrt{{\alpha }^2+4{\beta }^2}}{2}& 0& rcos\theta \\ 0& -\beta & 0& \frac{\alpha +\sqrt{{\alpha }^2+4{\beta }^2}}{2}& rsin\theta \\ \frac{\sqrt{{\alpha }^2+4{\beta }^2}-\alpha }{2}& 0& -\beta & 0& \frac{cos\theta }{r}\\ 0& \frac{\sqrt{{\alpha }^2+4{\beta }^2}-\alpha }{2}& 0& -\beta & \frac{sin\theta }{r}\\ -rcos\theta & -rsin\theta & -\frac{cos\theta }{r}& -\frac{sin\theta }{r}& 0\end{array}\right).\end{array}$$

Furthermore, one can get the characteristic polynomial of $J_S$,

$$\begin{array}{c}\hfill \lambda (\lambda +2\beta )({\lambda }^3+2\beta {\lambda }^2+(\frac{4\beta }{\alpha +\sqrt{{\alpha }^2+4{\beta }^2}}+\frac{\alpha }{\beta })\lambda +\alpha +\sqrt{{\alpha }^2+4{\beta }^2}+\frac{4{\beta }^2}{\alpha +\sqrt{{\alpha }^2+4{\beta }^2}})=0.\end{array}$$

According to Routh–Hurwitz criterion, it is unstable for any given $\alpha >0$ and $\beta >0$.

2.4. Chaotic Behavior and Attractors

Assuming that $\alpha =5$, $\beta =2$ and $x\left(0\right)={(5-3j,1-4j,5.5)}^T$, the methods numerical analysis are used to obtain chaotic attractor, Poincaré map and bifurcation diagrams, see Figure 2, Figure 3 and Figure 4. Figure 2 shows chaotic attractors of the complex Rikitake system in different planes. The Poincaré diagrams of system (3) are depicted in Figure 3. As described in Figure 4a, basic bifurcation versus parameter $\beta \in (0,3)$ with $\alpha =5$. Figure 4b demonstrates system (3) is sensitivity to initial value. Furthermore, we apply numerical computation to obtain the corresponding Lyapunov exponents of system (3),

$$L{E}_1=0.117534,L{E}_2=-0.043679,L{E}_3=-0.051743,L{E}_4=-3.957388,\mathrm{and}L{E}_5=-4.064725.$$

Thus, using the formula of fractal dimension [33], we easily deduce that

$$D=j+\frac{1}{|L{E}_{j+1}|}\sum _{i=1}^{j}L{E}_i=3+\frac{L{E}_1+L{E}_2+L{E}_3}{|L{E}_4|}=3.0055875.$$

This implies that the Lyapunov dimension of system (3) is fractional. Consequently, it is noticeable that system (3) behaves chaotically for this choice of $\alpha =5$ and $\beta =2$.

3. Chaos Control

Adaptive technique is adopted to investigate the control problem of the complex Rikitake system. Before giving the conclusion of this section, we first introduce a lemma.

Lemma 1

([9]). Consider the nonlinear system

$$\begin{array}{c}\hfill \dot{\vartheta }=\phi \left(\vartheta \right),\end{array}$$

(4)

where $\vartheta \in {\mathbb{R}}^n$ is the state, $\phi \left(\vartheta \right)\in {\mathbb{R}}^n$ is continuous function with $\phi \left(0\right)=0$. Suppose that there exists a nonsingular coordinate transformation $\upsilon =T\vartheta$, which can convert system (4) into two subsystems

$$\begin{array}{c}\hfill \begin{array}{c}{\dot{\mu }}_1=G_1({\mu }_1,{\mu }_2),\hfill \\ {\dot{\mu }}_2=G_2({\mu }_1,{\mu }_2),\hfill \end{array}\end{array}$$

where ${\mu }_1={({\upsilon }_1,{\upsilon }_2,\dots ,{\upsilon }_r)}^T,r\ge 1$, ${\mu }_2={({\upsilon }_{r+1},{\upsilon }_{r+2},\dots ,{\upsilon }_n)}^T$, $G_1({\mu }_1,{\mu }_2)\in {\mathbb{R}}^r$, $G_2({\mu }_1,{\mu }_2)\in {\mathbb{R}}^{n-r}$, and the subsystem

$$\begin{array}{c}\hfill {\dot{\mu }}_2=G_2(0,{\mu }_2)\end{array}$$

is globally asymptotically stable (GAS). Then the controller is designed as

$$\begin{array}{c}\hfill u={(k{\mu }_1,0)}^T\end{array}$$

and the adaptation law is in the form of

$$\begin{array}{c}\hfill \dot{k}=-\sigma {\mu }_1^T{\mu }_1,\end{array}$$

where $\sigma >0$ is an arbitrary real number. That is to say, the controlled system

$$\begin{array}{c}\hfill \begin{array}{c}{\dot{\mu }}_1=G_1({\mu }_1,{\mu }_2)+k{\mu }_1,\hfill \\ {\dot{\mu }}_2=G_2({\mu }_1,{\mu }_2)\hfill \end{array}\end{array}$$

is asymptotically stable.

As discussed in Section 2, the complex Rikitake system has no stable equilibrium point. Next, we design a feedback controller to stabilize the complex Rikitake system to any fixed points. The equilibrium point of system (3) is recorded as $S=(s_1,s_2,s_3,s_4,s_5)$. Making the following coordinate transformation:

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{x}_1={\tilde{x}}_1+s_1,\hfill \\ x_2={\tilde{x}}_2+s_2,\hfill \\ x_3={\tilde{x}}_3+s_3,\hfill \\ x_4={\tilde{x}}_4+s_4,\hfill \\ x_5={\tilde{x}}_5+s_5,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

we further have the controlled system

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{\tilde{x}}}_1=-\beta {\tilde{x}}_1-{\tilde{x}}_3{\tilde{x}}_5+{\tilde{x}}_3{s}_5+{\tilde{x}}_5{s}_3+u_1,\hfill \\ {\dot{\tilde{x}}}_2=-\beta {\tilde{x}}_2-{\tilde{x}}_4{\tilde{x}}_5+{\tilde{x}}_4{s}_5+{\tilde{x}}_5{s}_4+u_2,\hfill \\ {\dot{\tilde{x}}}_3=-\beta {\tilde{x}}_3+{\tilde{x}}_1{\tilde{x}}_5+{\tilde{x}}_1{s}_5+{\tilde{x}}_5{s}_1-\alpha {\tilde{x}}_1+u_3,\hfill \\ {\dot{\tilde{x}}}_4=-\beta {\tilde{x}}_4+{\tilde{x}}_2{\tilde{x}}_5+{\tilde{x}}_2{s}_5+{\tilde{x}}_5{s}_2-\alpha {\tilde{x}}_2+u_4,\hfill \\ {\dot{\tilde{x}}}_5=-{\tilde{x}}_1{\tilde{x}}_3-{\tilde{x}}_1{s}_3-{\tilde{x}}_3{s}_1-{\tilde{x}}_2{\tilde{x}}_4-{\tilde{x}}_2{s}_4-{\tilde{x}}_4{s}_2+u_5,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(5)

where $u={(u_1,u_2,u_3,u_4,u_5)}^T$ is the controller to be designed. Thus, the problem of stabilizing system (3) to the equilibrium point S is converted to that of stabilizing system (5) at the origin. By Lemma 1, we have the following result.

Theorem 1.

System (3) can be controlled to the equilibrium point $S=(s_1,s_2,s_3,s_4,s_5)$ by constructing the following adaptive feedback controller

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}u={(k{\tilde{x}}_1,k{\tilde{x}}_2,0,0,0)}^T,\hfill \\ \dot{k}=-\sigma ({\tilde{x}}_1^2+{\tilde{x}}_2^2),\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(6)

where $\sigma >0$ is a chosen positive real number.

Proof. 

It is noticeable that when ${\tilde{x}}_1={\tilde{x}}_2=0$, the remainder subsystem of system (5) without a controller becomes

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{\tilde{x}}}_3=-\beta {\tilde{x}}_3+{\tilde{x}}_5{s}_1,\hfill \\ {\dot{\tilde{x}}}_4=-\beta {\tilde{x}}_4+{\tilde{x}}_5{s}_2,\hfill \\ {\dot{\tilde{x}}}_5=-{\tilde{x}}_3{s}_1-{\tilde{x}}_4{s}_2.\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(7)

The coefficient matrix of system (7) is

$$\begin{array}{c}\hfill J=\left(\begin{array}{ccc}-\beta & 0& s_1\\ 0& -\beta & s_2\\ -s_1& -s_2& 0\end{array}\right),\end{array}$$

and its corresponding characteristic equation is described by

$$\begin{array}{c}\hfill (\lambda +\beta )({\lambda }^2+\beta \lambda +s_1^2+s_2^2)=0.\end{array}$$

(8)

Through the analysis of Equation (8), we conclude all roots have negative real part. Thus, according to Routh–Hurwitz criterion, system (7) is GAS. By Lemma 1, one deduces that system (5) can be regulated at the origin by controller (6), that is, system (3) tends to the equilibrium point S. □

By the same argument, when ${\tilde{x}}_3={\tilde{x}}_4=0$, the subsystem of system (5) without controller is of the form

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{\tilde{x}}}_1=-\beta {\tilde{x}}_1+{\tilde{x}}_5{s}_3,\hfill \\ {\dot{\tilde{x}}}_2=-\beta {\tilde{x}}_2+{\tilde{x}}_5{s}_4,\hfill \\ {\dot{\tilde{x}}}_5=-{\tilde{x}}_1{s}_3-{\tilde{x}}_2{s}_4,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

which is GAS. We derive another result on stabilization of the complex Rikitake system from Lemma 1.

Theorem 2.

System (3) can be regulated to the equilibrium point $S=(s_1,s_2,s_3,s_4,s_5)$ by constructing the following adaptive feedback controller

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}u={(0,0,k{\tilde{x}}_3,k{\tilde{x}}_4,0)}^T,\hfill \\ \dot{k}=-\sigma ({\tilde{x}}_3^2+{\tilde{x}}_4^2),\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(9)

where $\sigma >0$ is an arbitrary real number.

Based on the proposed results, we will now give a numerical description on controlling the complex Rikitake system. In the following two cases, choose the parameters as $\alpha =5$, $\beta =2$, and fix the initial values as $x\left(0\right)={(5-3j,1-4j,5.5)}^T$.

For the choice of $\theta =\frac{\pi }{2}$, an unstable equilibrium point can be obtained as $S_1=(0,1.689,0,0.5923,5.702)$. From Theorem 1, we derive the controller (6) with $\sigma \left(0\right)=1$ and $k\left(0\right)=-1$. From Figure 5, one can see that the complex Rikitake system (3) can be regulated to its equilibrium point $S_1$ via the controller (6).

For the choice of $\theta =0$, an unstable equilibrium point can be obtained as $S_2=(-1.689,0,-0.5923,0,5.702)$. From Theorem 2, we derive the controller (9) with $\sigma \left(0\right)=3$ and $k\left(0\right)=-2$. As shown in Figure 6, the complex Rikitake system (3) converges to its equilibrium point $S_2$.

4. Synchronization Scheme

This section proves the existence of synchronization of the complex Rikitake system, and then realizes CS and the coexistence of CS and AS by feedback control method.

Let us consider two identical complex Rikitake systems with different initial conditions. The drive system is described by

$$\begin{array}{c}\hfill \begin{array}{c}\dot{Y}=h\left(Y\right),\hfill \end{array}\end{array}$$

(10)

where $Y={(y_1,y_2,y_3)}^T\in {\mathbb{C}}^3$, $y_1=y_1^r+j{y}_1^i$, $y_2=y_2^r+j{y}_2^i$, $h\left(Y\right)={(h_1\left(Y\right),h_2\left(Y\right),h_3\left(Y\right))}^T\in {\mathbb{C}}^3$, $y_3\in \mathbb{R}$, and

$$\begin{array}{c}\hfill \begin{array}{c}{h}_1\left(Y\right)=-\beta y_1+y_2{y}_3,\hfill \\ h_2\left(Y\right)=-\beta y_2-\alpha y_1+y_1{y}_3,\hfill \\ h_3\left(Y\right)=1-\frac{1}{2}({\overline{y}}_1{y}_2+y_1{\overline{y}}_2).\hfill \end{array}\end{array}$$

In the same way, the response system with controllers can be expressed as

$$\begin{array}{c}\hfill \begin{array}{c}\dot{Z}=h\left(Z\right)+u,\hfill \end{array}\end{array}$$

(11)

where $Z={(z_1,z_2,z_3)}^T\in {\mathbb{C}}^3$, $z_1=z_1^r+j{z}_1^i$, $z_2=z_2^r+j{z}_2^i$, $u={(u_1,u_2,u_3)}^T\in {\mathbb{C}}^3$ is the error feedback controller to be designed, $z_3\in \mathbb{R}$, and

$$\begin{array}{c}\hfill \begin{array}{c}{h}_1\left(Z\right)=-\beta z_1+z_2{z}_3,\hfill \\ h_2\left(Z\right)=-\beta z_2-\alpha z_1+z_1{z}_3,\hfill \\ h_3\left(Z\right)=1-\frac{1}{2}({\overline{z}}_1{z}_2+z_1{\overline{z}}_2).\hfill \end{array}\end{array}$$

The synchronization error is denoted as

$$\begin{array}{c}\hfill \begin{array}{c}e=Z-\delta Y,\hfill \end{array}\end{array}$$

where $\delta =diag\{{\delta }_1,{\delta }_2,{\delta }_3\}$ and ${\delta }_i\ne 0$ are real constants ($i=1,2,3$).

Following the results in [12], we introduce the relevant definition.

Definition 1.

For the drive system (10) and the response system (11),

1. 

Systems (10) and (11) are said to be CS if there exists a diagonal matrix $\delta =I_3,$ i.e., ${\delta }_i=1$ ($i=1,2,3$), such that $\underset{t\to \infty }{lim}\left|\right|e\left(t\right)\left|\right|=0$;

2. 

Systems (10) and (11) are said to be AS if there exists a diagonal matrix $\delta =-I_3,$ i.e., ${\delta }_i=-1$ ($i=1,2,3$), such that $\underset{t\to \infty }{lim}\left|\right|e\left(t\right)\left|\right|=0$;

3. 

Systems (10) and (11) are said to achieve the coexistence of CS and AS if there exist some ${\delta }_i=1$ while the remaining ${\delta }_j=-1$ ($i\ne j$, $i,j=1,2,3$), such that $\underset{t\to \infty }{lim}\left|\right|e\left(t\right)\left|\right|=0$.

4.1. The Existence of Synchronization in the Complex Rikitake System

Taking the derivative of $e\left(t\right)$ and using Equations (10) and (11), one obtains

$$\begin{array}{c}\hfill \begin{array}{c}\dot{e}=h\left(Z\right)-\delta h\left(Y\right)+u,\hfill \end{array}\end{array}$$

which is equivalent to the following equations

$$\begin{array}{c}\hfill \begin{array}{c}{\dot{e}}^r=h^r\left(Z\right)-\delta h^r\left(Y\right)+u^r,\hfill \end{array}\end{array}$$

and

$$\begin{array}{c}\hfill \begin{array}{c}{\dot{e}}^i=h^i\left(Z\right)-\delta h^i\left(Y\right)+u^i.\hfill \end{array}\end{array}$$

It is clear that $e=0$ implies $e^r=0$ and $e^i=0$. In order to implement a suitable controller, $e^r=0$ should be a fixed point of the error system without controllers (i.e., $u^r=0$)

$$\begin{array}{c}\hfill \begin{array}{c}{\dot{e}}^r=h^r\left(Z\right)-\delta h^r\left(Y\right)\hfill \end{array}\end{array}$$

(12)

and $e^i=0$ should be a fixed point of the error system in absence of controllers (i.e., $u^i=0$)

$$\begin{array}{c}\hfill \begin{array}{c}{\dot{e}}^i=h^i\left(Z\right)-\delta h^i\left(Y\right).\hfill \end{array}\end{array}$$

(13)

Thus, one has

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{h}^r\left(\delta Y\right)=\delta h^r\left(Y\right),\hfill \\ h^i\left(\delta Y\right)=\delta h^i\left(Y\right).\hfill \end{array}\right.\hfill \end{array}\end{array}$$

Furthermore, the following equality holds

$$\begin{array}{c}\hfill \begin{array}{c}h\left(\delta Y\right)=\delta h\left(Y\right).\hfill \end{array}\end{array}$$

Thus, we obtain the conclusion about the existence of the synchronization problem.

Theorem 3.

The existence of synchronization in complex chaotic system (10) iff $h\left(\delta Y\right)=\delta h\left(Y\right)$ has solutions for δ.

Proof. 

The proof is easily obtained by Theorem 1 in [12], so it is omitted here. □

Using the result of Theorem 3, one gets that the existence of synchronization in the complex Rikitake system (10) is converted to the following equations having solutions for $\delta$,

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{h}_1\left(\delta Y\right)-{\delta }_1{h}_1\left(Y\right)=({\delta }_2{\delta }_3-{\delta }_1)y_2{y}_3=0,\hfill \\ h_2\left(\delta Y\right)-{\delta }_2{h}_2\left(Y\right)=({\delta }_1{\delta }_3-{\delta }_2)y_1{y}_3-\alpha ({\delta }_1-{\delta }_2)y_1=0,\hfill \\ h_3\left(\delta Y\right)-{\delta }_3{h}_3\left(Y\right)=1-{\delta }_3-\frac{1}{2}({\delta }_1{\delta }_2-{\delta }_3)({\overline{y}}_1{y}_2+y_1{\overline{y}}_2)=0,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

which leads to

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\delta }_1={\delta }_2,\hfill \\ {\delta }_3=1,\hfill \\ |{\delta }_i|=1,i=1,2.\hfill \end{array}\right.\hfill \end{array}\end{array}$$

Furthermore, we have the following results:

I

${\delta }_1={\delta }_2={\delta }_3=1$, which implies CS in the complex Rikitake system (10) occurs.

II

${\delta }_1={\delta }_2=-1,{\delta }_3=1$, which implies the coexistence of CS and AS in the complex Rikitake system (10) exists.

4.2. CS of the Complex Rikitake System

Now, we consider CS of two complex Rikitake systems (10) and (11). When ${\delta }_1={\delta }_2={\delta }_3=1$, the CS error is defined as $e=Z-Y$. The error system is calculated as

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{e}}_1={\dot{z}}_1-{\dot{y}}_1=-\beta e_1+e_2{e}_3+y_3{e}_2+y_2{e}_3+u_1,\hfill \\ {\dot{e}}_2={\dot{z}}_2-{\dot{y}}_2=-\beta e_2-\alpha e_1+e_1{e}_3+y_3{e}_1+y_1{e}_3+u_2,\hfill \\ {\dot{e}}_3={\dot{z}}_3-{\dot{y}}_3=-\frac{1}{2}({\overline{e}}_1{e}_2+e_1{\overline{e}}_2+{\overline{e}}_1{y}_2+e_1{\overline{y}}_2+{\overline{y}}_1{e}_2+y_1{\overline{e}}_2)+u_3,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

which can be rewritten as

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{e}}_1^r=-\beta e_1^r+e_2^r{e}_3+y_3{e}_2^r+y_2^r{e}_3+u_1^r,\hfill \\ {\dot{e}}_1^i=-\beta e_1^i+e_2^i{e}_3+y_3{e}_2^i+y_2^i{e}_3+u_1^i,\hfill \\ {\dot{e}}_2^r=-\beta e_2^r-\alpha e_1^r+e_1^r{e}_3+y_3{e}_1^r+y_1^r{e}_3+u_2^r,\hfill \\ {\dot{e}}_2^i=-\beta e_2^i-\alpha e_1^i+e_1^i{e}_3+y_3{e}_1^i+y_1^i{e}_3+u_2^i,\hfill \\ {\dot{e}}_3=-(e_1^r{e}_2^r+e_1^i{e}_2^i+e_1^r{y}_2^r+e_1^i{y}_2^i+y_1^r{e}_2^r+y_1^i{e}_2^i)+u_3,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(14)

where $\tilde{u}={(u_1^r,u_1^i,u_2^r,u_2^i,u_3)}^T$ is a real controller to be designed. Thus, on the basis of Lemma 1, one has the following result.

Theorem 4.

Two identical complex Rikitake systems (10) and (11) can realize CS via the following adaptive controller

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}\tilde{u}={(0,0,k_1{e}_2^r,k_1{e}_2^i,0)}^T,\hfill \\ \dot{k}=-\sigma ({\left(e_2^r\right)}^2+{\left(e_2^i\right)}^2),\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(15)

where $\sigma >0$ is an arbitrary real number.

Proof. 

Let us consider the uncontrolled error dynamical system (14). It is clear that if $e_2^r=e_2^i=0$, then the subsystem of uncontrolled system (14) reads as

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{e}}_1^r=-\beta e_1^r+y_2^r{e}_3,\hfill \\ {\dot{e}}_1^i=-\beta e_1^i+y_2^i{e}_3,\hfill \\ {\dot{e}}_3=-y_2^r{e}_1^r-y_2^i{e}_1^i,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

which is GAS. From Lemma 1, system (14) with controller (15) approaches to the zero equilibrium point, i.e., CS of two identical complex Rikitake systems (10) and (11) can be realized by the designed controller (15). □

In the same argument, when $e_1^r=e_1^i=0$, the subsystem of system (14) in absence of controller is presented as

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{e}}_2^r=-\beta e_2^r+y_1^r{e}_3,\hfill \\ {\dot{e}}_2^i=-\beta e_2^i+y_1^i{e}_3,\hfill \\ {\dot{e}}_3=-y_1^r{e}_2^r-y_1^i{e}_2^i,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

which is GAS. Thus, the following result is deduced.

Theorem 5.

Two identical complex Rikitake systems (10) and (11) can realize CS by designing adaptive controller

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}\tilde{u}={(k{e}_1^r,k{e}_1^i,0,0,0)}^T,\hfill \\ \dot{k}=-\sigma ({\left(e_1^r\right)}^2+{\left(e_1^i\right)}^2),\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(16)

where $\sigma >0$ is an arbitrary real number.

In the next part, by giving the initial conditions as $y\left(0\right)={(15-3j,1-4j,5.5)}^T$, $z\left(0\right)={(4-j,2-3j,-0.3)}^T$, $k\left(0\right)=-1$, $\sigma =2$, and constructing controller (15), we have simulation results which are shown by the following Figure 7 and Figure 8. Figure 7a displays that the errors $e_1^r$, $e_1^i$, $e_2^r$, $e_2^i$ and $e_3$ can been regulated to the zero equilibrium point. Figure 8 depicts that state variables of system (11) are completely synchronized with state variables of system (10). That is, two identical complex Rikitake systems realize CS.

4.3. The Coexistence of CS and AS in the Complex Rikitake System

When ${\delta }_1={\delta }_2=-1$, ${\delta }_3=1$, AS error is denoted as $E_1=z_1+y_1$ and $E_2=z_2+y_2$, CS error is denoted as $e_3=z_3-y_3$. It is easy to obtain the following error dynamical system

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{E}}_1={\dot{z}}_1+{\dot{y}}_1=-\beta E_1+E_2{e}_3+y_3{E}_2-y_2{e}_3+u_1,\hfill \\ {\dot{E}}_2={\dot{z}}_2+{\dot{y}}_2=-\beta E_2-\alpha E_1+E_1{e}_3+y_3{E}_1-y_1{e}_3+u_2,\hfill \\ {\dot{e}}_3={\dot{z}}_3-{\dot{y}}_3=-\frac{1}{2}({\overline{E}}_1{E}_2+E_1{\overline{E}}_2-{\overline{E}}_1{y}_2-E_1{\overline{y}}_2-{\overline{y}}_1{E}_2-y_1{\overline{E}}_2)+u_3,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

which turns into

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{E}}_1^r=-\beta E_1^r+E_2^r{e}_3+y_3{E}_2^r-y_2^r{e}_3+u_1^r,\hfill \\ {\dot{E}}_1^i=-\beta E_1^i+E_2^i{e}_3+y_3{E}_2^i-y_2^i{e}_3+u_1^i,\hfill \\ {\dot{E}}_2^r=-\beta E_2^r-\alpha E_1^r+E_1^r{e}_3+y_3{E}_1^r-y_1^r{e}_3+u_2^r,\hfill \\ {\dot{E}}_2^i=-\beta E_2^i-\alpha E_1^i+E_1^i{e}_3+y_3{E}_1^i-y_1^i{e}_3+u_2^i,\hfill \\ {\dot{e}}_3=-E_1^r{E}_2^r-E_1^i{E}_2^i+E_1^r{y}_2^r+E_1^i{y}_2^i+y_1^r{E}_2^r+y_1^i{E}_2^i+u_3.\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(17)

Theorem 6.

Two identical complex Rikitake systems (10) and (11) can achieve the coexistence of CS and AS by virtue of the following adaptive controller

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}\tilde{u}={(0,0,k{E}_2^r,k{E}_2^i,0)}^T,\hfill \\ \dot{k}=-\sigma ({\left(E_2^r\right)}^2+{\left(E_2^i\right)}^2),\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(18)

where $\sigma >0$ is an arbitrary real number.

Proof. 

Let us consider the system (17) in absence of controller. Obviously, when $E_2^r=E_2^i=0$, the subsystem of system (17) without controller can be converted to

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{E}}_1^r=-\beta E_1^r-y_2^r{e}_3,\hfill \\ {\dot{E}}_1^i=-\beta E_1^i-y_2^i{e}_3,\hfill \\ {\dot{e}}_3=y_2^r{E}_1^r+y_2^i{E}_1^i,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(19)

which is GAS. From Lemma 1, system (17) can be governed at the origin by controller (18). That is to say, the coexistence of CS and AS in two identical complex Rikitake systems (10) and (11) can be realized by adaptive controller (18). □

Similarly, when $E_1^r=E_1^i=0$, the subsystem of system (17) without controller is described by

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}{\dot{E}}_2^r=-\beta E_2^r-y_1^r{e}_3,\hfill \\ {\dot{E}}_2^i=-\beta E_2^i-y_1^i{e}_3,\hfill \\ {\dot{e}}_3=y_1^r{E}_2^r+y_1^i{E}_2^i,\hfill \end{array}\right.\hfill \end{array}\end{array}$$

which is GAS. Thus, by means of Lemma 1, we obtain another result.

Theorem 7.

Two identical complex Rikitake systems (10) and (11) can realize the coexistence of CS and AS by designing the following controller

$$\begin{array}{c}\hfill \begin{array}{c}\left\{\begin{array}{c}\tilde{u}={(k{E}_1^r,k{E}_1^i,0,0,0)}^T,\hfill \\ \dot{k}=-\sigma ({\left(E_1^r\right)}^2+{\left(E_1^i\right)}^2),\hfill \end{array}\right.\hfill \end{array}\end{array}$$

(20)

where $\sigma >0$ is an arbitrary real number.

For numerical simulations, fix the initial values as $y\left(0\right)={(18+2j,1+2j,3)}^T$ and $z\left(0\right)={(4-j,1+2j,-0.3)}^T$. By constructing controller (20) with $k\left(0\right)=-4$ and $\sigma =5$, we can obtain the simulation results, see Figure 9 and Figure 10. As one can see from Figure 9 the errors $E_1^r$, $E_1^i$, $E_2^r$, $E_2^i$ and $e_3$ can be regulated to the zero equilibrium point. Figure 10 describes the change of state variables of systems (10) and (11). It is easy to see that $z_1^r$, $z_1^i$, $z_2^r$ and $z_2^i$ of system (11) anti-synchronize $y_1^r$, $y_1^i$, $y_2^r$ and $y_2^i$ of system (10) respectively, while $z_3$ of system (11) synchronizes completely with $y_3$ of system (10). Therefore, the coexistence of CS and AS in two identical complex Rikitake systems can be realized.

5. Conclusions

This paper centers on control and synchronization of a new complex chaotic system. Firstly, we propose a complex Rikitake system and investigate its dynamical behavior. Then, by means of feedback control, we design controllers to regulate the complex Rikitake system to any equilibrium points. Thus, we not only prove the existence of synchronization in the complex Rikitake system but also construct adaptive controllers to realize two types of synchronization schemes, such as CS and the coexistence of CS and AS. It is notable that the presented scheme is a single and linear feedback controller and it is easy to implement in engineering. Therefore, the control method will be widely applied in practice in the future.