# Dynamics and Synchronization of a Memristor-Based Chaotic System with No Equilibrium

Authors: Hong-Min Li; Yan-Feng Yang; Yang Zhou; Chun-Lai Li; Kun Qian; Zhao-Yu Li; Jian-Rong Du
Journal: Complexity
 (2019)

---

## Abstract

The topics of memristive system and synchronization are two hot fields of research in nonlinear dynamics. In this paper, we introduce a memristor-based chaotic system with no equilibrium. It is found that the memristor-based system under investigation exhibits fruitful dynamic behaviors such as coexisting bifurcation, multistability, transient chaos, and transient quasiperiod. Thus, it is difficult to reproduce the accurate dynamics of the system, which is highly advantageous in encryption and communication. Then, a simple intermittent control scheme with adaptive mechanism is developed to achieve complete synchronization for the introduced system. Because the output signal is transmitted intermittently to the receiver system, more channel capacity can be saved and the security performance can be improved naturally in practical communication.


---

## Body

## 1. Introduction

As the fourth basic circuital element along with resistor, inductor, and capacitor, the memristor was postulated by Chua in 1971 [1], and it was then successfully fabricated by the HP laboratories in 2008 [2]. Since then, the memristor was recognized to perfect the symmetry of the four fundamental circuital variables and has aroused wide interest in academia [3–5]. The memristor is commonly defined as a two-terminal nonlinear component with controllable resistance called memristance that varies according to the amount of charge or flux flowing through it [6]. The fingerprint of a memristor is composed of a current-voltage characteristic curve, which shows a pinched hysteresis loop whose shape varies with frequency and converges to a straight line with the increase of frequency [7].The memristor is currently used to design flash memory, improve neural networks, and construct chaotic circuits, for the intrinsic characteristics of memory, nanoscale device, and inherent nonlinearity. Itoh and Chua constructed the memristive chaotic oscillator in 2008, by replacing Chua’s diodes in Chua’s circuit with the piecewise linear memristor [8]. Afterwards, many memristor-based chaotic oscillators were constructed. For example, by replacing the single diode with a memristor in the original circuit, Pelap postulated an emendatory Tamasevicius oscillator [9]. Bi-Rong designed a simple chaotic circuit consisting of an inductor, a capacitor, and a voltage-controlled memristor [10]. Zhao et al. proposed a memristor-based chaotic system by replacing the nonlinear diode in the Chua circuit with an active flux-controlled memristor [11]. In order to increase the complexity of memristor-based system, Teng et al. used a fourth-degree polynomial memristance to produce a multiscroll chaotic attractor [12]. By replacing Chua’s diode with a physical SBT memristor and a negative conductance in the canonical Chua’s circuit, a new memristor-based modified Chua’s circuit is constructed [13]. There usually emerges special dynamics in this kind of memristive systems, such as initial sensitivity, coexisting bifurcation, coexistence attractors, and transient dynamics. Therefore, the memristive chaotic system will provide more complex dynamics and facilitates the engineering applications of information encryption, secure communication, and signal processing [14–17].Meanwhile, close attention was paid to chaotic system without equilibrium [18–20]. From the computer-processing perspective, it is challenging to numerically localize the attractor in such system since there is no transient process leading from the vicinity of unstable equilibrium point. In other words, the attracting basin of such system does not intersect with any small neighborhood of its equilibrium point, or the attractor is “hidden” [21–24]. Up to now, little information is known about the dynamical behavior in such system, and what is worse is that the Shilnikov criteria cannot be employed to prove the chaos for the absence of heteroclinic or homoclinic orbit [25, 26].Because of its application in secure communication, digital signal, neural network, and other fields, the synchronization of chaotic system is a fashionable subject in nonlinear science. Since the first scheme was carried out by Pecora and Carroll for the synchronization of two identical chaotic systems [27], a great diversity of methods have been proposed to synchronize chaotic systems, such as active control, adaptive control, impulsive control, sliding mode control, intermittent control, pinning control, and hybrid control [28–31]. Generally, a chaotic communication system can be constructed based on master-slave synchronization, where the message is modulated by the transmitting system and is then sent to the receiving system. Also, in the receiver, the designed synchronization scheme is used to demodulate the received signal and extract the message [32]. The intermittent synchronization implies that the slave system receives the demodulated information from the master system intermittently. Therefore, the intermittent synchronization scheme will decrease the amount of conveyed information and the communication channel capacity will be reserved for more message transmission. Also, accordingly, the security of the chaotic communication system will be improved since the redundancy of the synchronization information in the channel is reduced. Therefore, the intermittent synchronization scheme is especially fit for the design of practical chaos-based communication system.In this paper, we introduce a memristor-based chaotic system with no equilibrium. The dynamical evolution of the memristive system is studied by using phase diagram, time-domain trajectory, bifurcation diagram, and Lyapunov exponent. It is found that by changing system parameters or initial condition, the reported system exhibits different topological structures of coexisting bifurcation, multistability, transient dynamics. The coexisting hidden attractors signify that the system has fruitful and complex dynamic behaviors, which is highly advantageous in encryption and communication for the difficulty of reproducing the accurate dynamics of the system. Then, a simple intermittent control scheme with adaptive mechanism is developed to achieve complete synchronization for the introduced memristive system. Since the output signal is transmitted intermittently to the receiver system, more channel capacity can be saved and the security performance of the communication system can be improved naturally in practical communication. Theoretical analysis and illustrative examples are executed to verify the effectiveness of the proposed synchronization scheme.


## 2. Memristor-Based Chaotic System with No Equilibrium

### 2.1. Model Description

Based on Sprott A system, the constructed memristive chaotic system can be described by the following differential equations:(1)x˙=y,y˙=−Wx+yz,z˙=c−y2,where x, y, z are state variables; the function W represents the model of a flux-controlled memristor, depicted as W=3ax2+b; and a, b, c are positive parameters.The dissipativity is decided by∇V=∂x˙/x+∂y˙/y+∂z˙/z=z; thus, system (1) is non-Hamiltonian conservative of phase volume [33]. It is palpable that there exists no equilibrium in system (1). Therefore, the strange attractor is “hidden” in the sense of classification method described by Leonov et al. [21], and the Shilnikov method cannot be employed to verify the emergence of chaos since there is no heteroclinic or homoclinic orbit in this system. It is easy to know that the system is symmetric with respect to the y-axis in the sense of coordinate transformation (x, y, z, t) ⟶ (−x, y, −z, −t).When choosing the parametersa = 3, b = −1, and c = 1 and initial condition (0.2, 0.15, 0), system (1) appears a chaotic state with the Lyapunov exponents 0.1062, 0, −0.1062, as illustrated by the phase portrait in Figure 1.Figure 1
Phase portrait projected onto the plane of (a)x–y; (b) y–z; (c) x–z.
(a)
(b)
(c)


### 2.2. Coexisting Bifurcation and Multiple Attractors

It is found that the memristive system under consideration can experience rich bifurcation structures when continuously monitoring the bifurcation parameter. Also, the memristive system has completely different bifurcation behaviors when the initial conditions are set to different values.We assign the parametersb and c of system (1) as b = −1 and c = 1, and select parameter a serving as the representative bifurcation parameter. The coexisting bifurcation diagrams, produced by the local maxima of the state variable z in terms of control parameter a, are depicted in Figure 2(a) when the system starts with the initial states (0.1, 0, 0) and (0.2, 0, 0). Also, the corresponding maximal Lyapunov exponents are depicted in Figure 2(b). This strategy represents a convenient and intuitive approach to identify the window in which the multiple coexisting attractors arise. We further draw in Figure 3 the enlarged bifurcation diagrams of Figure 2(a) to show the typical regions of multiple coexisting attractors. By taking different initial conditions (0.1, 0, 0) and (0.2, 0, 0), respectively, we plot the coexisting multiple attractors with different parameter a in Figure 4. It is found that there may emerge rich dynamical structures of coexisting chaos, quasiperiod or period with different shape for the same parameter a when starting from different initial states.Figure 2
(a) Bifurcation diagrams and (b) spectrums of maximal Lyapunov exponent versus parametera, the initial conditions are (0.1, 0, 0) for blue diagrams and (0.2, 0, 0) for red diagrams.
(a)
(b)Figure 3
Enlargement of the bifurcation diagrams of Figure2(a) in ranges of (a) [0.8, 1.2] and (b) [3, 4]; the initial conditions are (0.1, 0, 0) for blue diagrams and (0.2, 0, 0) for red diagrams.
(a)
(b)Figure 4
Coexistence of different attractors with different parametera. (a) Quasiperiods with different shape; (b) chaos and quasiperiod; (c) chaos and complicated period; (d) chaos and period-1; (e) chaos and quasiperiod; (f) period and chaos, for the initial conditions (0.1, 0, 0) and (0.2, 0, 0), respectively.
(a)
(b)
(c)
(d)
(e)
(f)We also study the dynamics evolution of system (1) by using the initial value x0 served as the bifurcation parameter. The system parameters are fixed as a = 3, b = −1, c = 1, and the rest initial conditions are assigned to be y0 = 0, z0 = 0. The bifurcation diagram and spectra of Lyapunov exponent are depicted in Figures 5(a) and 5(b), respectively. It can be found that when increasing initial value x0 from −2 to 2, there emerge periodic windows embedded in the chaotic region, and the dynamics is symmetrically distributed with respect to zero value. In fact, the dynamics is also symmetrically distributed with respect to y-axis and z-axis, as depicted by the dynamical map of x0 versus y0 and the dynamical map of y0 versus z0 in Figures 6(a) and 6(b), respectively. In the dynamical map, the system is chaotic in the cyan region and periodic in the pink region. Some representative coexistence attractors of symmetric distribution with respect to x0 are displayed in Figure 7.Figure 5
(a) Bifurcation diagram; (b) Lyapunov exponent spectrum versusx0.
(a)
(b)Figure 6
Dynamic distributions of system (1) with a = 3, b = −1, c = 1: (a) x0 versus y0 with z0 = 0; (b)y0 versus z0 with x0 = 0.1. The chaotic regions are shown in cyan and the periodic regions are shown in pink.
(a)
(b)Figure 7
Symmetric coexisting attractors with respect tox0. (a) Period with x0 = 1.3 and −1.3; (b) chaos with x0 = 1.6 and −1.6.
(a)
(b)


### 2.3. Transient Dynamics

It is surprising to see in Figure2(b) that a periodic window appears in the parameter region of 2.586 ≤ a ≤ 2.688, but it does show the chaotic behavior in Figure 2(a). The emergence of different dynamical modes is due to the transition from the long term transient period to steady chaos with the time evolutions of system. The transient dynamics can be represented by chaotic orbit before entering the final nonchaotic behavior, and the inverse process is also correct.Firstly, the case of system parametersa = 3, b = −1, and c = 1 and initial condition x0 = 0.1, y0 = 0.15, z0 = 0.1 is considered. The time trajectory in the region of [0 s, 800 s] and the phase diagrams in two different time intervals of [0 s, 400 s] and [450 s, 800 s] are depicted in Figure 8, which illustrates the dynamics transformation from transient period to steady chaos.Figure 8
(a) Time-domain waveform ofz in the region of [0 s, 800 s]; (b) the phase portrait in time interval of [0 s, 400 s]; (c) the phase portrait in time interval of [450 s, 800 s].
(a)
(b)
(c)Then we take the selection of system parametersa = 3, b = −1, and c = 1 and initial condition x0 = 0.3, y0 = 0.1, z0 = 0.6. The time trajectory in the region of [0 s, 2000 s] and the phase diagrams in two different time intervals of [0 s, 900 s] and [1100 s, 2000 s] are depicted in Figure 9. It is observed in Figure 9 that the trajectory of system (1) starts from a quasiperiodic orbit for a long time and then it transforms into a chaotic state at t = 1050 s. Similarly to transient chaos, we call this dynamic phenomenon as transient quasiperiod.Figure 9
(a) Time-domain waveform ofz in the region of [0 s, 2000 s]; (b) the phase portrait in time interval of [0 s, 900 s]; (c) the phase portrait in time interval of [1100 s, 2000 s].
(a)
(b)
(c)


## 2.1. Model Description

Based on Sprott A system, the constructed memristive chaotic system can be described by the following differential equations:(1)x˙=y,y˙=−Wx+yz,z˙=c−y2,where x, y, z are state variables; the function W represents the model of a flux-controlled memristor, depicted as W=3ax2+b; and a, b, c are positive parameters.The dissipativity is decided by∇V=∂x˙/x+∂y˙/y+∂z˙/z=z; thus, system (1) is non-Hamiltonian conservative of phase volume [33]. It is palpable that there exists no equilibrium in system (1). Therefore, the strange attractor is “hidden” in the sense of classification method described by Leonov et al. [21], and the Shilnikov method cannot be employed to verify the emergence of chaos since there is no heteroclinic or homoclinic orbit in this system. It is easy to know that the system is symmetric with respect to the y-axis in the sense of coordinate transformation (x, y, z, t) ⟶ (−x, y, −z, −t).When choosing the parametersa = 3, b = −1, and c = 1 and initial condition (0.2, 0.15, 0), system (1) appears a chaotic state with the Lyapunov exponents 0.1062, 0, −0.1062, as illustrated by the phase portrait in Figure 1.Figure 1
Phase portrait projected onto the plane of (a)x–y; (b) y–z; (c) x–z.
(a)
(b)
(c)


## 2.2. Coexisting Bifurcation and Multiple Attractors

It is found that the memristive system under consideration can experience rich bifurcation structures when continuously monitoring the bifurcation parameter. Also, the memristive system has completely different bifurcation behaviors when the initial conditions are set to different values.We assign the parametersb and c of system (1) as b = −1 and c = 1, and select parameter a serving as the representative bifurcation parameter. The coexisting bifurcation diagrams, produced by the local maxima of the state variable z in terms of control parameter a, are depicted in Figure 2(a) when the system starts with the initial states (0.1, 0, 0) and (0.2, 0, 0). Also, the corresponding maximal Lyapunov exponents are depicted in Figure 2(b). This strategy represents a convenient and intuitive approach to identify the window in which the multiple coexisting attractors arise. We further draw in Figure 3 the enlarged bifurcation diagrams of Figure 2(a) to show the typical regions of multiple coexisting attractors. By taking different initial conditions (0.1, 0, 0) and (0.2, 0, 0), respectively, we plot the coexisting multiple attractors with different parameter a in Figure 4. It is found that there may emerge rich dynamical structures of coexisting chaos, quasiperiod or period with different shape for the same parameter a when starting from different initial states.Figure 2
(a) Bifurcation diagrams and (b) spectrums of maximal Lyapunov exponent versus parametera, the initial conditions are (0.1, 0, 0) for blue diagrams and (0.2, 0, 0) for red diagrams.
(a)
(b)Figure 3
Enlargement of the bifurcation diagrams of Figure2(a) in ranges of (a) [0.8, 1.2] and (b) [3, 4]; the initial conditions are (0.1, 0, 0) for blue diagrams and (0.2, 0, 0) for red diagrams.
(a)
(b)Figure 4
Coexistence of different attractors with different parametera. (a) Quasiperiods with different shape; (b) chaos and quasiperiod; (c) chaos and complicated period; (d) chaos and period-1; (e) chaos and quasiperiod; (f) period and chaos, for the initial conditions (0.1, 0, 0) and (0.2, 0, 0), respectively.
(a)
(b)
(c)
(d)
(e)
(f)We also study the dynamics evolution of system (1) by using the initial value x0 served as the bifurcation parameter. The system parameters are fixed as a = 3, b = −1, c = 1, and the rest initial conditions are assigned to be y0 = 0, z0 = 0. The bifurcation diagram and spectra of Lyapunov exponent are depicted in Figures 5(a) and 5(b), respectively. It can be found that when increasing initial value x0 from −2 to 2, there emerge periodic windows embedded in the chaotic region, and the dynamics is symmetrically distributed with respect to zero value. In fact, the dynamics is also symmetrically distributed with respect to y-axis and z-axis, as depicted by the dynamical map of x0 versus y0 and the dynamical map of y0 versus z0 in Figures 6(a) and 6(b), respectively. In the dynamical map, the system is chaotic in the cyan region and periodic in the pink region. Some representative coexistence attractors of symmetric distribution with respect to x0 are displayed in Figure 7.Figure 5
(a) Bifurcation diagram; (b) Lyapunov exponent spectrum versusx0.
(a)
(b)Figure 6
Dynamic distributions of system (1) with a = 3, b = −1, c = 1: (a) x0 versus y0 with z0 = 0; (b)y0 versus z0 with x0 = 0.1. The chaotic regions are shown in cyan and the periodic regions are shown in pink.
(a)
(b)Figure 7
Symmetric coexisting attractors with respect tox0. (a) Period with x0 = 1.3 and −1.3; (b) chaos with x0 = 1.6 and −1.6.
(a)
(b)


## 2.3. Transient Dynamics

It is surprising to see in Figure2(b) that a periodic window appears in the parameter region of 2.586 ≤ a ≤ 2.688, but it does show the chaotic behavior in Figure 2(a). The emergence of different dynamical modes is due to the transition from the long term transient period to steady chaos with the time evolutions of system. The transient dynamics can be represented by chaotic orbit before entering the final nonchaotic behavior, and the inverse process is also correct.Firstly, the case of system parametersa = 3, b = −1, and c = 1 and initial condition x0 = 0.1, y0 = 0.15, z0 = 0.1 is considered. The time trajectory in the region of [0 s, 800 s] and the phase diagrams in two different time intervals of [0 s, 400 s] and [450 s, 800 s] are depicted in Figure 8, which illustrates the dynamics transformation from transient period to steady chaos.Figure 8
(a) Time-domain waveform ofz in the region of [0 s, 800 s]; (b) the phase portrait in time interval of [0 s, 400 s]; (c) the phase portrait in time interval of [450 s, 800 s].
(a)
(b)
(c)Then we take the selection of system parametersa = 3, b = −1, and c = 1 and initial condition x0 = 0.3, y0 = 0.1, z0 = 0.6. The time trajectory in the region of [0 s, 2000 s] and the phase diagrams in two different time intervals of [0 s, 900 s] and [1100 s, 2000 s] are depicted in Figure 9. It is observed in Figure 9 that the trajectory of system (1) starts from a quasiperiodic orbit for a long time and then it transforms into a chaotic state at t = 1050 s. Similarly to transient chaos, we call this dynamic phenomenon as transient quasiperiod.Figure 9
(a) Time-domain waveform ofz in the region of [0 s, 2000 s]; (b) the phase portrait in time interval of [0 s, 900 s]; (c) the phase portrait in time interval of [1100 s, 2000 s].
(a)
(b)
(c)


## 3. Synchronization Control of Memristor-Based Chaotic System

### 3.1. Synchronization Scheme

We consider the master-slave synchronization scheme for the introduced chaotic system, and the corresponding master-slave systems are described by the compact form:(2)master systemdxdt=fx,x=x1,x2,x3=x,y,z,(3)slave systemdydt=fy,y=y1,y2,y3=x′,y′,z′.In which,f·∈R3 is the smooth vector field satisfying the Lipschitz condition:(4)fix−fiy≤kxi−yi≤kx−y∞,x−y∞=maxixi−yi,i=1,2,3.To realize the synchronization of systems (2) and (3), we add a single linear controller to the i-th equation of the slave system, as depicted by(5)dydt=fy+ui,ui=μxi−yi.The synchronization error is defined byei=xi−yi,i=1,2,3. Also, we construct the candidate Lyapunov function as V=1/2∑i=13xi−yi2. The time derivative of V along the synchronization error is deduced by(6)V˙=∑i=13xi−yix˙i−y˙i=∑i=13xi−yifix−fy−ui=∑i=13xi−yifix−fy−μxi−yi∑i=13xi−yi≤3k−μx−y∞2.Whenμ≥3k, we have V˙≤0. Thus, the controlled slave system will asymptotically synchronize with the master system with the simple controller ui=μxi−yi when μ≥3k.In fact, to reduce the control consumption, we can optimize the controller as(7)ui=μxi−yi⋅Hxi−yi−ε,where ε is a small positive constant; the function H⋅ is described as Hz = 1 when z ≥ 0 and Hz = 0 when z < 0.The practical significance of the optimized control scheme is that one imposes the controllerμxi−yi to the slave system when xi−yi≥ε, but the controller does not work when xi−yi<ε. Thus, the controller can realize the system synchronization intermittently with the adaptive mechanism, according to the characteristics of motion trajectories. Therefore, compared with continuous synchronization schemes, intermittent synchronization will reduce the amount of conveyed information, which is of significance in the practical communication since the communication channel capacity will be reserved for more message transmission. In addition, the security of chaotic communication system will be improved due to the reduction of redundancy of synchronization information in the channel.


### 3.2. Numerical Simulation

We impose the controllerμx2−y2⋅Hx2−y2−ε to the second term of the slave system. The system parameters are set as a = 3, b = −1, and c = 1.We first choose the initial condition of system (2) as x0=0.1,0,0 with which system (2) is chaotic, and the initial states of system (4) are taken as y0=−0.1,0,0 for also displaying chaotic, as shown in Figure 10. The synchronization result is shown in Figure 11 when the controller gain μ equals 3 and ε is set to be 0.02.Figure 10
(a) Phase diagram and (b) state trajectories withx0=0.1,0,0 and y0=−0.1,0,0.
(a)
(b)Figure 11
(a) Time response of the states and (b) synchronization error withμ = 3, x0=0.1,0,0 and y0=−0.1,0,0.
(a)
(b)Then, we also setx0=0.1,0,0 with which system (2) is chaotic, but the initial states of system (4) are taken as y0=0.1,1,0 for displaying quasiperiodic, as shown in Figure 12. The synchronization result is shown in Figure 13 when controller parameters μ = 2 and ε = 0.02.Figure 12
(a) Phase diagram and (b) state trajectories withx0=0.1,0,0 and y0=0.1,1,0.
(a)
(b)Figure 13
(a) Time response of the states and (b) synchronization error withμ = 2, x0=0.1,0,0 and y0=0.1,1,0.
(a)
(b)We know that no matter what the dynamic state of the memristive system is, the synchronization control of the memory system can be easily realized by adopting the designed method.


## 3.1. Synchronization Scheme

We consider the master-slave synchronization scheme for the introduced chaotic system, and the corresponding master-slave systems are described by the compact form:(2)master systemdxdt=fx,x=x1,x2,x3=x,y,z,(3)slave systemdydt=fy,y=y1,y2,y3=x′,y′,z′.In which,f·∈R3 is the smooth vector field satisfying the Lipschitz condition:(4)fix−fiy≤kxi−yi≤kx−y∞,x−y∞=maxixi−yi,i=1,2,3.To realize the synchronization of systems (2) and (3), we add a single linear controller to the i-th equation of the slave system, as depicted by(5)dydt=fy+ui,ui=μxi−yi.The synchronization error is defined byei=xi−yi,i=1,2,3. Also, we construct the candidate Lyapunov function as V=1/2∑i=13xi−yi2. The time derivative of V along the synchronization error is deduced by(6)V˙=∑i=13xi−yix˙i−y˙i=∑i=13xi−yifix−fy−ui=∑i=13xi−yifix−fy−μxi−yi∑i=13xi−yi≤3k−μx−y∞2.Whenμ≥3k, we have V˙≤0. Thus, the controlled slave system will asymptotically synchronize with the master system with the simple controller ui=μxi−yi when μ≥3k.In fact, to reduce the control consumption, we can optimize the controller as(7)ui=μxi−yi⋅Hxi−yi−ε,where ε is a small positive constant; the function H⋅ is described as Hz = 1 when z ≥ 0 and Hz = 0 when z < 0.The practical significance of the optimized control scheme is that one imposes the controllerμxi−yi to the slave system when xi−yi≥ε, but the controller does not work when xi−yi<ε. Thus, the controller can realize the system synchronization intermittently with the adaptive mechanism, according to the characteristics of motion trajectories. Therefore, compared with continuous synchronization schemes, intermittent synchronization will reduce the amount of conveyed information, which is of significance in the practical communication since the communication channel capacity will be reserved for more message transmission. In addition, the security of chaotic communication system will be improved due to the reduction of redundancy of synchronization information in the channel.


## 3.2. Numerical Simulation

We impose the controllerμx2−y2⋅Hx2−y2−ε to the second term of the slave system. The system parameters are set as a = 3, b = −1, and c = 1.We first choose the initial condition of system (2) as x0=0.1,0,0 with which system (2) is chaotic, and the initial states of system (4) are taken as y0=−0.1,0,0 for also displaying chaotic, as shown in Figure 10. The synchronization result is shown in Figure 11 when the controller gain μ equals 3 and ε is set to be 0.02.Figure 10
(a) Phase diagram and (b) state trajectories withx0=0.1,0,0 and y0=−0.1,0,0.
(a)
(b)Figure 11
(a) Time response of the states and (b) synchronization error withμ = 3, x0=0.1,0,0 and y0=−0.1,0,0.
(a)
(b)Then, we also setx0=0.1,0,0 with which system (2) is chaotic, but the initial states of system (4) are taken as y0=0.1,1,0 for displaying quasiperiodic, as shown in Figure 12. The synchronization result is shown in Figure 13 when controller parameters μ = 2 and ε = 0.02.Figure 12
(a) Phase diagram and (b) state trajectories withx0=0.1,0,0 and y0=0.1,1,0.
(a)
(b)Figure 13
(a) Time response of the states and (b) synchronization error withμ = 2, x0=0.1,0,0 and y0=0.1,1,0.
(a)
(b)We know that no matter what the dynamic state of the memristive system is, the synchronization control of the memory system can be easily realized by adopting the designed method.


## 4. Conclusions

In this paper, we introduce a memristor-based chaotic system with no equilibrium. Various tools including phase diagram, time-domain trajectory, bifurcation diagram, and Lyapunov exponent are exploited to establish the connection between the system parameters and dynamical behaviors. It is found that the reported system exhibits complex dynamics such as coexisting bifurcation, multistability, symmetric coexisting attractors, and transient dynamics, which is helpful for the security improvement of encryption and communication due to the difficulty of reproducing the accurate dynamics. Then, a simple control scheme with single linear couple is developed to achieve complete synchronization for the memristive system. Since the output signal is transmitted intermittently to the receiver system with the adaptive mechanism, the communication channel capacity will be reserved for more message transmission. Also, the security of chaotic communication system will be improved for the reduction of redundancy of synchronization information in the channel.



---

*Source: 4647608-2019-10-28.xml*