# Composite One- to Six-Scroll Hidden Attractors in a Memristor-Based Chaotic System and Their Circuit Implementation

Authors: Ying Li; Xiaozhu Xia; Yicheng Zeng; Qinghui Hong
Journal: Complexity
 (2020)

---

## Abstract

Chaotic systems with hidden multiscroll attractors have received much attention in recent years. However, most parts of hidden multiscroll attractors previously reported were repeated by the same type of attractor, and the composite of different types of attractors appeared rarely. In this paper, a memristor-based chaotic system, which can generate composite attractors with one up to six scrolls, is proposed. These composite attractors have different forms, similar to the Chua’s double scroll and jerk double scroll. Through theoretical analysis, we find that the new system has no fixed point; that is to say, all of the composite multiscroll attractors are hidden attractors. Additionally, some complicated dynamic behaviors including various hidden coexisting attractors, extreme multistability, and transient transition are explored. Moreover, hardware circuit using discrete components is implemented, and its experimental results supported the numerical simulations results.


---

## Body

## 1. Introduction

Recently, much effort has been devoted to the analysis of various chaotic systems owing to its potential applications [1–4]. To our knowledge, chaotic attractors are categorized as either self-excited attractors [5–7] or hidden attractors [8–11]. It should be particularly pointed out that hidden attractors have neither homoclinic nor heteroclinic orbits; thus, the Shil’nikov theorem [12] cannot be utilized to verify the existence of chaos, so that hidden chaotic systems attract great attention of scholars, and increasing number of researchers begin to study the rich dynamic behaviors of them.As is known to all, multiscroll chaotic systems have aroused extensive interests for their much complicated dynamical properties than single-scroll ones. So, it is a more desirable task to explore the hidden multiscroll chaotic system, and there have been emerging related literatures on this topic. In 2016, Jafari et al. presented a novel no-equilibrium chaotic system with multiscroll hidden chaotic sea by introducing a sine function into a 3D chaotic system [13]. In the same year, Hu et al. constructed two simple 3D chaotic systems without equilibrium based on an improved Sprott A system by adding a nonlinear function, from which hidden multiscroll attractors can be obtained [14]. Later, a new class of PWL dynamical system without equilibrium whose hidden chaotic attractors can display grid multiscroll has been introduced [15]. Deng and Wang proposed a multiscroll hidden chaotic system that has only two stable node-foci equilibrium points [16]. Hong et al. investigated a novel method for designing multidirection multibutterfly chaotic attractors (MDMBCA) without reconstructing nonlinear functions [17]. Furthermore, some memristive chaotic systems, which can generate hidden multiscroll or multiwing chaotic attractors, are reported [18, 19]. However, all of the hidden multiscroll attractors discussed above are of the same type attractor repeated in unidirectional or multidirection ones. Are there hidden multiscroll attractors, which are composed of different attractors? We are surprised that the answer is positive.If an attractor is compounded from two or more different single-scroll attractors, which is called composite attractor, its dynamic characteristics are more complicated. The composite attractor has been already discovered in some self-excited chaotic systems. For example, Cang et al. proposed a Lorenz-like system with composite structure by controlling system parameters [20]. Zhang et al. reported an autonomous-system-based approach for creating composite chaotic attractors from a class of generalized Lorenz systems via switching control [21]. And then, Xiong et al. constructed a chaotic system, which can generate composite four-scroll attractor by adding a symbolic function into a 3D jerk system [22]. However, although the complexity of dynamic behavior of these chaotic systems, which can produce composite attractors, has increased, the operation methods and state equations have also been correspondingly complicated, and the circuit experiment is difficult to realize. Besides, these composite multiscroll chaotic attractors mentioned above are limited to self-excited attractors, and the hidden composite multiscroll attractors have not been reported up to now. It is obvious that designing a simple chaotic system, which can generate hidden composite multiscroll attractors, is a meaningful task.Based on these considerations, we introduce a quadratic flux-controlled memristor into the 3D jerk chaotic system to construct a memristor-based hidden composite multiscroll chaotic system, which can produce composite one- to six-scroll hidden attractors through changing only one system parameter. It is pointed out that the composite six-scroll hidden attractor is composed of a Chua’s double-scroll attractor and two jerk double-scroll attractors. Compared with [22], this new memristor-based hyperchaotic system has more complex dynamic behaviors including hyperchaotic composite multiscroll hidden attractors and extreme multistability phenomenon. These composite multiscroll attractors break the previous pattern of the same single-scroll attractor that appeared repeatedly and has not been reported in previous literature.The rest of this paper is organized as follows. In Section2, the mathematical model of the proposed memristive hyperchaotic system and its typical attractor are presented. In Section 3, complex hidden dynamic behaviors of this new system including controllable multiscroll hidden chaotic or hyperchaotic attractors, coexistence of hidden attractors, extreme multistability, and transient transition behavior are described. In Section 4, the corresponding circuit implementation and experiment results are presented. The conclusion is presented in the last section.


## 2. Memristive Hyperchaotic System and Its Typical Attractors

In 2016, Kengne et al. performed a simple 3D autonomous jerk system with cubic nonlinearity [23], which is described as(1)x˙=y,y˙=az,z˙=x−by−z−x3,where x, y, z are state variables, and a, b are positive tunable parameters. In order to construct a memristor-based hyperchaotic system, a quadratic flux-controlled memristor model [24] described in equation (2) is added into system (1):(2)u˙=v,i=m+nuv,where u, v and i denote the internal state of the memristor, the input of the memristor, and the output of the memristor model, respectively. The m, n are two real parameters of the memristor model. So, we constructed a 4D memristor-based chaotic system, which can generate composite multiscroll hidden attractors, and it is described as(3)x˙=y+z,y˙=az−sgnz+p,z˙=2x−by−z−kWuy−x3,u˙=cy,in which x, y, z are state variables, a, b, c, k are real parameters, W (u) = (m + nu) is the memductance function, and p is positive control parameter. From a dissipative perspective, the divergence of the new system (3) can easily be obtained as follows:(4)∇V=∂x˙∂x+∂y˙∂y+∂z˙∂z+∂u˙∂u=−1,at any given point (x, y, z, u) of the state space. Consequently, the general condition of dissipativity related to the existence of attractive sets in our model is satisfied. The equilibrium points of system (3) can be obtained by solving the following equations: x˙=y˙=z˙=u˙=0. It is obvious that the equations have no real solutions when parameter p is positive. That is to say, system (3) has no equilibrium, and it belongs to hidden chaotic system.When fixing the control parametersa = 3.5, b = 0.8, m = 0.1, n = 0.3, p=0.01, k = 0.2, c = 0.01 and selecting the initial conditions as (0.1, 0.1, 0.1, 0.1), a composite six-scroll hidden chaotic attractor and its Poincaré map on y = 0 section as well as power spectrum are displayed in Figures 1(a)–1(c). As is depicted in Figure 1(a), the six-scroll attractor is composed of a Chua’s double-scroll attractor and two jerk double-scroll attractors, and the two small jerk double-scroll attractors are inside the large Chua’s double-scroll attractor, respectively. It is surprising that the composite six-scroll hidden chaotic attractors have never appeared before as far as we know.Figure 1
(a) The composite six-scroll hidden chaotic attractor. (b) The corresponding Poincaré mapping ony = 0 section. (c) The corresponding power spectrum.
(a)(b)(c)


## 3. Dynamic Analysis and Numerical Investigations

In this section, multifarious complicated dynamic behaviors including controllable multiscroll hidden attractors, coexisting of hidden chaotic or hyperchaotic attractors, extreme multistability, and transient transition behaviors of proposed system (3) are investigated in detail by means of Lyapunov exponent spectra, bifurcation diagrams, phase portraits, and time series. Note that all numerical simulations are carried out on MATLAB software, and the first 20000 points are discarded in the phase portraits.


### 3.1. Controllable Multiscroll Hidden Chaotic or Hyperchaotic Attractors

One interesting phenomenon is that the newly proposed system (3) can generate composite one- to six-scroll hidden attractors by changing only one system parameter. Fixing the control parameters b = 0.8, m = 0.1, n = 0.3, k = 0.2, c = 0.01, p=0.01 and changing the parameter a in the interval of [1.3, 5] with the initial conditions (0.1, 0.1, 0.1, 0.1), the corresponding bifurcation diagram is drawn in Figure 2(a). Additionally, using the famous Wolf method [25], the corresponding Lyapunov exponent spectrum proves the complexity of system (3) as depicted in Figure 2(b).Figure 2
Dynamic behavior with increasing the parameter a in the interval of [1.3, 5]. (a) Bifurcation diagrams ofx. (b) Lyapunov exponent spectrum.
(a)(b)In reference to Figure2(b), system (3) undergoes chaotic and hyperchaotic routes as follows: chaotic (1.3 ≤ a < 1.5) ⟶ hyperchaotic (1.5 ≤ a < 2.65) ⟶ chaotic (2.56 ≤ a < 3.06) ⟶ hyperchaotic (3.06 ≤ a < 3.42) ⟶ chaotic (3.42 ≤ a < 5). What surprised us is that system (3) generates composite one- to six-scrolls- hidden attractors when selecting appropriate parameters, and these attractors may exhibit either a symmetric or asymmetric structure. It is rare to see such attractors in an asymmetric system. Some typical hidden chaotic or hyperchaotic attractors with different values of parameter a are shown in Figures 3(a)–3(l). Note that Figures 3(a) and 3(b) display a pair of asymmetric single-scroll hidden chaotic attractors at a = 1.38 and 2.92. Figures 3(c) and 3(d) show a pair of symmetric double-scroll hidden hyperchaotic attractors at a = 2.3 and 2.55, and the corresponding Lyapunov exponents are (0.213, 0.179, 0, −1.381) and (0.255, 0.209, 0, −1.444), respectively. Figures 3(e) and 3(f) exhibit two different double-scroll hidden chaotic attractors at a = 2.946 and 5. And a pair of asymmetric composite 3-scroll hidden chaotic and hyperchaotic attractors are drawn in Figures 3(g) and 3(h) with a = 2.974 and 3.117. Furthermore, a pair of symmetric composite 4-scroll hidden chaotic and hyperchaotic attractors, a composite 5-scroll hidden chaotic attractor, and a composite 6-scroll hidden hyperchaotic attractor are presented in Figures 3(i) and 3(l) with a = 3.03, 3.076, 3.055, and 3.125, respectively. To explain more intuitively and clearly, the different types of attractors with different values of parameter a are listed in Table 1. Consequently, the unusual controllable multiscroll hidden chaotic or hyperchaotic attractors of system (3) are innovative.Figure 3
Some typical phase portraits with different parametera. (a, b) A pair of asymmetric single-scroll chaotic attractors. (b) A pair of asymmetric single-scroll chaotic attractors. (c, d) A pair of symmetric double-scroll chaotic attractors. (e, f) Two different double-scroll chaotic attractors. (g, h) Two different composite 3-scroll chaotic attractors. (i, j) A pair of symmetric composite 4-scroll chaotic and hyperchaotic attractors. (k) A composite 5-scroll chaotic attractor. (l) A composite 6-scroll hyperchaotic attractor.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)Table 1
Composite one- to six-scroll attractors of system (3) with variation of parameter a.
ParameteraLESChaotic or hyperchaoticThe number of scrolls of hidden attractorDiagrams1.38(+, 0, −, −)ChaoticSingle-scrollFigure3(a)2.92(+, 0, −, −)ChaoticSingle-scrollFigure3(b)2.3(+, +, 0, −)HyperchaoticDouble-scrollFigure3(c)2.55(+, +, 0, −)HyperchaoticDouble-scrollFigure3(d)2.946(+, 0, −, −)ChaoticDouble-scrollFigure3(e)5(+, 0, −, −)ChaoticDouble-scrollFigure3(f)2.974(+, 0, −, −)Chaotic3-scrollFigure3(g)3.117(+, +, 0, −)Hyperchaotic3-scrollFigure3(h)3.03(+, 0, −, −)Chaotic4-scrollFigure3(i)3.076(+, +, 0, −)Hyperchaotic4-scrollFigure3(j)3.055(+, 0, −, −)Chaotic5-scrollFigure3(k)3.125(+, +, 0, −)Hyperchaotic6-scrollFigure3(l)


### 3.2. Coexistence of Hidden Attractors

Nowadays, attractor coexisting is one of the common nonlinear dynamic phenomena, which is a chaotic system performing more than one concurrent attractor for a given set of system parameters under the different initial conditions [26–32]. In this section, the coexistence of hidden attractors with some specific parameters setting and initial conditions in system (3) is discussed. By selecting the control parameters as b = 0.8, m = 0.1, n = 0.3, k = 0.2, c = 0.01, p=0.01, the coexisting bifurcation diagram of the state variable x with respect to the parameter a is described in Figure 4(a). Figure 4(b) displays the Lyapunov exponent spectrum under the initial conditions (−0.1, −0.1, −0.1, −0.1), which is in accordance with the Lyapunov exponent spectrum under the initial conditions (0.1, 0.1, 0.1, 0.1) in Figure 2(b). Various composite multiscroll chaotic or hyperchaotic hidden coexisting attractors are shown in Figure 5; they are including symmetric attractors and asymmetric attractors coexisting, odd number of attractors, and even number of attractors coexisting. The more detailed information is displayed in Table 2. All in all, it is informative and rare that so many composite multiscroll attractors coexist, which injects new blood into the field of chaos coexisting.Figure 4
Coexisting bifurcation behaviors with increasing the parameter a in the interval of [1.3, 5]. (a) Coexisting bifurcation diagrams of the state variable (x). (b) Lyapunov exponent spectrum with initial conditions (−0.1, −0.1, −0.1, −0.1).
(a)(b)Figure 5
The coexisting chaotic or hyperchaotic hidden attractors with different values of parametera.
(a)(b)(c)(d)(e)(f)(g)(h)(i)Table 2
The coexisting multiscroll attractors of system (3) with different values of parameter a.
ParameteraChaotic or hyperchaoticCoexisting attractorsDiagrams1.38ChaoticA single scroll and a double-scrollFigure5(a)2.3HyperchaoticTwo symmetric double-scrollFigure5(b)2.92ChaoticTwo symmetric single-scrollFigure5(c)2.974ChaoticA 3-scroll and a single-scrollFigure5(d)3.055ChaoticA 5-scroll and a 3-scrollFigure5(e)3.095HyperchaoticA 3-scroll and a double-scrollFigure5(f)3.105HyperchaoticA 5-scroll and a double-scrollFigure5(g)3.117HyperchaoticA 3-scroll and a 4-scrollFigure5(h)3.125HyperchaoticA 6-scroll and a double-scrollFigure5(i)


### 3.3. Extreme Multistability Relying on Memristor Initial Condition

In this subsection, we further studied the extreme multistability phenomenon [33–38] in system (3). To reveal extreme multistability phenomenon relying on memristor initial condition, the typical parameters a = 3.05, b = 0.8, c = 0.01, m = 0.1, n = 0.3, k = 0.2, p=0.01 and initial conditions x (0) = 0.1, y (0) = 0.1, z (0) = 0.1 are fixed. When the memristor initial condition u (0) is increased gradually from −4.8 to 10, the bifurcation diagram of the state variable x and the Lyapunov exponent spectrum are plotted in Figures 6(a) and 6(b), respectively. As can be seen from Figure 6, the system has extremely rich dynamic characteristics.Figure 6
Bifurcation diagram and the Lyapunov exponent spectrum with respect to the initial conditionsu (0) in region of [−4.8, 10].
(a)(b)Figure 7
The coexistence of infinite number of attractors onx–z plane with different initial condition u (0).
(a)(b)(c)(d)(e)(f)To explain the coexistence of infinite number of attractors in system (2) more intuitively, several special values of u are selected as shown in Table 3 and the phase diagrams of the infinite coexisting attractors as shown in Figure 7.Table 3
Extreme multistability behaviors of system (3).
Initial conditionsTypes of coexisting multiple attractorsDiagrams(0.1, 0.1, 0.1, −1.76); (0.1, 0.1, 0.1, 0.8); (0.1, 0.1, 0.1, 0.14); (0.1, 0.1, 0.1, 0.59)Two single scroll and two double-scrollFigure7(a) red, green, cyan, magenta.(0.1, 0.1, 0.1, −4.3); (0.1, 0.1, 0.1, 0.4); (0.1, 0.1, 0.1, 7.2)Three double-scrollFigure7(b) magenta, green, cyan.(0.1, 0.1, 0.1, 0.23); (0.1, 0.1, 0.1, 0.37)Two 3-scrollFigure7(c) magenta, green.(0.1, 0.1, 0.1, −0.18); (0.1, 0.1, 0.1,−0.1)Two 4-scrollFigure7(d) magenta, green.(0.1, 0.1, 0.1, −1.29); (0.1, 0.1, 0.1, − 0.42)Two 5-scrollFigure7(e) magenta, green.(0.1, 0.1, 0.1, −3.4); (0.1, 0.1, 0.1, −1.8)Two 6-scrollFigure7(f) magenta, green.Specially, when settinga = 20, b = 0.8, m = 0.1, n = 0.3, k = 0.2, c = 0.01, p=0.01 and initial conditions x (0) = 0.1, y (0) = 0.1, z (0) = 0.1, and change the memristor initial condition u (0) in the region of [0, 8], the bifurcation diagram is plotted in Figure 8. In light of Figure 8, it is clear that the system goes through reverse period doubling bifurcation route from chaos to periodic. In order to illustrate this phenomenon of extreme multistability, there are seven types of hidden attractors, which are presented in Figure 9. In fact, more different types of hidden attractors can be obtained under the different initial conditions. In general, Figures 7–9 intuitively reveal the coexisting infinitely many hidden hyperchaotic attractors, which indicates that system (3) possesses hidden extreme multistability indeed.Figure 8
Bifurcation diagram with respect to the initial conditionsu (0) in region of [0, 8].Figure 9
The coexistence of infinite number of attractors onx-z plane with different initial condition u (0), the different attractors corresponding to the different values of u (0) are marked with different colors.


### 3.4. Transient Transition Behavior

In system exhibiting transient transition, the orbit is from one state (chaotic, periodic) to another state (chaotic, periodic) for a finite time interval before settling into a final state. Transient chaos with boundary crisis is often encountered in nonlinear dynamic systems [22, 39–43]. For the sake of verifying transient transition behavior in system (3), the system parameters are considered as a = 3.115, b = 0.8, m = 0.1, n = 0.3, k = 0.2, c = 0.01, p=0.01, while the initial conditions are selected as (0.1, 0.1, 0.1, 0.1) and the time is t = 1000 s, and a phase diagram of six scrolls and the corresponding time-domain waveform with state variables x and z are shown in Figures 10(a)–10(c). Figure 11 displays three different phase diagrams with different simulation time and the Lyapunov exponent spectrum with increased time. From the Lyapunov exponent spectrum in Figure 11(d), it is clearly seen that there are two Lyapunov exponents greater than zero at t = 80 s; that is to say, system (3) transitions from chaotic to hyperchaotic in the time of t = 80 s.Figure 10
The transient transition behaviors of system (3). (a) A phase diagram of six-scroll hidden chaotic attractors. (b) Time-domain waveform of the variable (x). (c) Time-domain waveform of the variable z.
(a)(b)(c)Figure 11
The transient transition behaviors with system (3). (a) A phase diagram of double-scroll hidden hyperchaotic attractors with t = 80 s. (b) A phase diagram of 3-scroll hidden chaotic attractors with t = 250 s. (c) A phase diagram of 5-scroll hidden hyperchaotic attractors with t = 750 s. (d) The Lyapunov exponent spectrum with increased time in the region of [0, 1000].
(a)(b)(c)(d)


## 3.1. Controllable Multiscroll Hidden Chaotic or Hyperchaotic Attractors

One interesting phenomenon is that the newly proposed system (3) can generate composite one- to six-scroll hidden attractors by changing only one system parameter. Fixing the control parameters b = 0.8, m = 0.1, n = 0.3, k = 0.2, c = 0.01, p=0.01 and changing the parameter a in the interval of [1.3, 5] with the initial conditions (0.1, 0.1, 0.1, 0.1), the corresponding bifurcation diagram is drawn in Figure 2(a). Additionally, using the famous Wolf method [25], the corresponding Lyapunov exponent spectrum proves the complexity of system (3) as depicted in Figure 2(b).Figure 2
Dynamic behavior with increasing the parameter a in the interval of [1.3, 5]. (a) Bifurcation diagrams ofx. (b) Lyapunov exponent spectrum.
(a)(b)In reference to Figure2(b), system (3) undergoes chaotic and hyperchaotic routes as follows: chaotic (1.3 ≤ a < 1.5) ⟶ hyperchaotic (1.5 ≤ a < 2.65) ⟶ chaotic (2.56 ≤ a < 3.06) ⟶ hyperchaotic (3.06 ≤ a < 3.42) ⟶ chaotic (3.42 ≤ a < 5). What surprised us is that system (3) generates composite one- to six-scrolls- hidden attractors when selecting appropriate parameters, and these attractors may exhibit either a symmetric or asymmetric structure. It is rare to see such attractors in an asymmetric system. Some typical hidden chaotic or hyperchaotic attractors with different values of parameter a are shown in Figures 3(a)–3(l). Note that Figures 3(a) and 3(b) display a pair of asymmetric single-scroll hidden chaotic attractors at a = 1.38 and 2.92. Figures 3(c) and 3(d) show a pair of symmetric double-scroll hidden hyperchaotic attractors at a = 2.3 and 2.55, and the corresponding Lyapunov exponents are (0.213, 0.179, 0, −1.381) and (0.255, 0.209, 0, −1.444), respectively. Figures 3(e) and 3(f) exhibit two different double-scroll hidden chaotic attractors at a = 2.946 and 5. And a pair of asymmetric composite 3-scroll hidden chaotic and hyperchaotic attractors are drawn in Figures 3(g) and 3(h) with a = 2.974 and 3.117. Furthermore, a pair of symmetric composite 4-scroll hidden chaotic and hyperchaotic attractors, a composite 5-scroll hidden chaotic attractor, and a composite 6-scroll hidden hyperchaotic attractor are presented in Figures 3(i) and 3(l) with a = 3.03, 3.076, 3.055, and 3.125, respectively. To explain more intuitively and clearly, the different types of attractors with different values of parameter a are listed in Table 1. Consequently, the unusual controllable multiscroll hidden chaotic or hyperchaotic attractors of system (3) are innovative.Figure 3
Some typical phase portraits with different parametera. (a, b) A pair of asymmetric single-scroll chaotic attractors. (b) A pair of asymmetric single-scroll chaotic attractors. (c, d) A pair of symmetric double-scroll chaotic attractors. (e, f) Two different double-scroll chaotic attractors. (g, h) Two different composite 3-scroll chaotic attractors. (i, j) A pair of symmetric composite 4-scroll chaotic and hyperchaotic attractors. (k) A composite 5-scroll chaotic attractor. (l) A composite 6-scroll hyperchaotic attractor.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)(k)(l)Table 1
Composite one- to six-scroll attractors of system (3) with variation of parameter a.
ParameteraLESChaotic or hyperchaoticThe number of scrolls of hidden attractorDiagrams1.38(+, 0, −, −)ChaoticSingle-scrollFigure3(a)2.92(+, 0, −, −)ChaoticSingle-scrollFigure3(b)2.3(+, +, 0, −)HyperchaoticDouble-scrollFigure3(c)2.55(+, +, 0, −)HyperchaoticDouble-scrollFigure3(d)2.946(+, 0, −, −)ChaoticDouble-scrollFigure3(e)5(+, 0, −, −)ChaoticDouble-scrollFigure3(f)2.974(+, 0, −, −)Chaotic3-scrollFigure3(g)3.117(+, +, 0, −)Hyperchaotic3-scrollFigure3(h)3.03(+, 0, −, −)Chaotic4-scrollFigure3(i)3.076(+, +, 0, −)Hyperchaotic4-scrollFigure3(j)3.055(+, 0, −, −)Chaotic5-scrollFigure3(k)3.125(+, +, 0, −)Hyperchaotic6-scrollFigure3(l)


## 3.2. Coexistence of Hidden Attractors

Nowadays, attractor coexisting is one of the common nonlinear dynamic phenomena, which is a chaotic system performing more than one concurrent attractor for a given set of system parameters under the different initial conditions [26–32]. In this section, the coexistence of hidden attractors with some specific parameters setting and initial conditions in system (3) is discussed. By selecting the control parameters as b = 0.8, m = 0.1, n = 0.3, k = 0.2, c = 0.01, p=0.01, the coexisting bifurcation diagram of the state variable x with respect to the parameter a is described in Figure 4(a). Figure 4(b) displays the Lyapunov exponent spectrum under the initial conditions (−0.1, −0.1, −0.1, −0.1), which is in accordance with the Lyapunov exponent spectrum under the initial conditions (0.1, 0.1, 0.1, 0.1) in Figure 2(b). Various composite multiscroll chaotic or hyperchaotic hidden coexisting attractors are shown in Figure 5; they are including symmetric attractors and asymmetric attractors coexisting, odd number of attractors, and even number of attractors coexisting. The more detailed information is displayed in Table 2. All in all, it is informative and rare that so many composite multiscroll attractors coexist, which injects new blood into the field of chaos coexisting.Figure 4
Coexisting bifurcation behaviors with increasing the parameter a in the interval of [1.3, 5]. (a) Coexisting bifurcation diagrams of the state variable (x). (b) Lyapunov exponent spectrum with initial conditions (−0.1, −0.1, −0.1, −0.1).
(a)(b)Figure 5
The coexisting chaotic or hyperchaotic hidden attractors with different values of parametera.
(a)(b)(c)(d)(e)(f)(g)(h)(i)Table 2
The coexisting multiscroll attractors of system (3) with different values of parameter a.
ParameteraChaotic or hyperchaoticCoexisting attractorsDiagrams1.38ChaoticA single scroll and a double-scrollFigure5(a)2.3HyperchaoticTwo symmetric double-scrollFigure5(b)2.92ChaoticTwo symmetric single-scrollFigure5(c)2.974ChaoticA 3-scroll and a single-scrollFigure5(d)3.055ChaoticA 5-scroll and a 3-scrollFigure5(e)3.095HyperchaoticA 3-scroll and a double-scrollFigure5(f)3.105HyperchaoticA 5-scroll and a double-scrollFigure5(g)3.117HyperchaoticA 3-scroll and a 4-scrollFigure5(h)3.125HyperchaoticA 6-scroll and a double-scrollFigure5(i)


## 3.3. Extreme Multistability Relying on Memristor Initial Condition

In this subsection, we further studied the extreme multistability phenomenon [33–38] in system (3). To reveal extreme multistability phenomenon relying on memristor initial condition, the typical parameters a = 3.05, b = 0.8, c = 0.01, m = 0.1, n = 0.3, k = 0.2, p=0.01 and initial conditions x (0) = 0.1, y (0) = 0.1, z (0) = 0.1 are fixed. When the memristor initial condition u (0) is increased gradually from −4.8 to 10, the bifurcation diagram of the state variable x and the Lyapunov exponent spectrum are plotted in Figures 6(a) and 6(b), respectively. As can be seen from Figure 6, the system has extremely rich dynamic characteristics.Figure 6
Bifurcation diagram and the Lyapunov exponent spectrum with respect to the initial conditionsu (0) in region of [−4.8, 10].
(a)(b)Figure 7
The coexistence of infinite number of attractors onx–z plane with different initial condition u (0).
(a)(b)(c)(d)(e)(f)To explain the coexistence of infinite number of attractors in system (2) more intuitively, several special values of u are selected as shown in Table 3 and the phase diagrams of the infinite coexisting attractors as shown in Figure 7.Table 3
Extreme multistability behaviors of system (3).
Initial conditionsTypes of coexisting multiple attractorsDiagrams(0.1, 0.1, 0.1, −1.76); (0.1, 0.1, 0.1, 0.8); (0.1, 0.1, 0.1, 0.14); (0.1, 0.1, 0.1, 0.59)Two single scroll and two double-scrollFigure7(a) red, green, cyan, magenta.(0.1, 0.1, 0.1, −4.3); (0.1, 0.1, 0.1, 0.4); (0.1, 0.1, 0.1, 7.2)Three double-scrollFigure7(b) magenta, green, cyan.(0.1, 0.1, 0.1, 0.23); (0.1, 0.1, 0.1, 0.37)Two 3-scrollFigure7(c) magenta, green.(0.1, 0.1, 0.1, −0.18); (0.1, 0.1, 0.1,−0.1)Two 4-scrollFigure7(d) magenta, green.(0.1, 0.1, 0.1, −1.29); (0.1, 0.1, 0.1, − 0.42)Two 5-scrollFigure7(e) magenta, green.(0.1, 0.1, 0.1, −3.4); (0.1, 0.1, 0.1, −1.8)Two 6-scrollFigure7(f) magenta, green.Specially, when settinga = 20, b = 0.8, m = 0.1, n = 0.3, k = 0.2, c = 0.01, p=0.01 and initial conditions x (0) = 0.1, y (0) = 0.1, z (0) = 0.1, and change the memristor initial condition u (0) in the region of [0, 8], the bifurcation diagram is plotted in Figure 8. In light of Figure 8, it is clear that the system goes through reverse period doubling bifurcation route from chaos to periodic. In order to illustrate this phenomenon of extreme multistability, there are seven types of hidden attractors, which are presented in Figure 9. In fact, more different types of hidden attractors can be obtained under the different initial conditions. In general, Figures 7–9 intuitively reveal the coexisting infinitely many hidden hyperchaotic attractors, which indicates that system (3) possesses hidden extreme multistability indeed.Figure 8
Bifurcation diagram with respect to the initial conditionsu (0) in region of [0, 8].Figure 9
The coexistence of infinite number of attractors onx-z plane with different initial condition u (0), the different attractors corresponding to the different values of u (0) are marked with different colors.


## 3.4. Transient Transition Behavior

In system exhibiting transient transition, the orbit is from one state (chaotic, periodic) to another state (chaotic, periodic) for a finite time interval before settling into a final state. Transient chaos with boundary crisis is often encountered in nonlinear dynamic systems [22, 39–43]. For the sake of verifying transient transition behavior in system (3), the system parameters are considered as a = 3.115, b = 0.8, m = 0.1, n = 0.3, k = 0.2, c = 0.01, p=0.01, while the initial conditions are selected as (0.1, 0.1, 0.1, 0.1) and the time is t = 1000 s, and a phase diagram of six scrolls and the corresponding time-domain waveform with state variables x and z are shown in Figures 10(a)–10(c). Figure 11 displays three different phase diagrams with different simulation time and the Lyapunov exponent spectrum with increased time. From the Lyapunov exponent spectrum in Figure 11(d), it is clearly seen that there are two Lyapunov exponents greater than zero at t = 80 s; that is to say, system (3) transitions from chaotic to hyperchaotic in the time of t = 80 s.Figure 10
The transient transition behaviors of system (3). (a) A phase diagram of six-scroll hidden chaotic attractors. (b) Time-domain waveform of the variable (x). (c) Time-domain waveform of the variable z.
(a)(b)(c)Figure 11
The transient transition behaviors with system (3). (a) A phase diagram of double-scroll hidden hyperchaotic attractors with t = 80 s. (b) A phase diagram of 3-scroll hidden chaotic attractors with t = 250 s. (c) A phase diagram of 5-scroll hidden hyperchaotic attractors with t = 750 s. (d) The Lyapunov exponent spectrum with increased time in the region of [0, 1000].
(a)(b)(c)(d)


## 4. Circuit Implementation and Experiment Results

Circuit implementation provides an alternative approach to explore the novel no-equilibrium memristive hyperchaotic system (3), and it is the key to the chaos application. In this part, a suitable electrical circuit is designed and a hardware circuit is implemented to verify the theoretical results obtained previously. The circuit diagram of the proposed system (3) is provided in Figure 12, and the circuits in the dashed box are constructed for the symbolic function circuit and memristor equivalent circuit.Figure 12
Circuit implementation of system (3).The operational amplifiers of Figure12 are LF353, whose supply voltages are ±15 V. And multipliers A1, A2, A3 are AD633 with output coefficient of 0.1. Based on system (3), select the system parameters as a = 3.5, b = 0.8, m = 0.1, n = 0.3, k = 0.2, c = 0.01, p=0.01 and let γ = γ0 t, where γ0 = 10000 is the time-scale transformation factor and system (3) can be rewritten by(5)x˙=10000y+10000z,y˙=35000z−10000signz+100,z˙=20000x−8000y−10000z−2000Wuy−10000x3,v˙=100y,where W(u) = (m + nu). Considering Kirchhoff’s law, the circuit equation of system (3) can be written as(6)dvxdt=1R1C1vy−1R2C1vz,dvydt=1R3C2vy−1R4C2signvz+1R5C2V1,dvzdt=1R6C3vx−1R7C3vy−1R8C3vz−1RmC3vy−110RnC3vuvy−1100R9C3vx3,dvudt=1R10Cuvy.Here, vx, vy, vz and vu are the voltages on capacitors C1, C2, C3 and Cu, respectively, and the capacitors C1 = C2 = C3 = Cu = 10 nF, V1 = −1 V. Comparing equation (5) with equation (6) and keeping the corresponding coefficients equal, one gets R1 = R2 = R4 = R8 = 10 kΩ, R3 = 2.86 kΩ, R5 = R10 = 10 MΩ, R6 = 5 kΩ, R7 = 12.5 kΩ, R9 = 100 Ω, Rm = 500 kΩ, Rn = 16.7 kΩ. Then, according to the circuit diagram in Figure 12, some off-the-shelf discrete components are used on the breadboard to build the hardware circuit as shown in Figure 13 and the parameters in the circuit are set as above. We can see a composite six-scroll chaotic attractor as shown in Figure 14(a) by using digital oscilloscope, which is well consistent with the phase diagram simulated with MATLAB in Figure 1(a). Lastly, keeping the values of R1, R2, R4, R5, R6, R7, R8, R9, Rm and Rn unchanged, and adjusting the values of R3 to appropriate resistance, other composite multiscroll attractors can be obtained from this hardware circuit as shown in Figures 14(b)–14(f), which also agree with the results of MATLAB simulation aforementioned. All of what is discussed above demonstrates the effectiveness and feasibility of system (3).Figure 13
The hardware circuit test setup to realize system (3).Figure 14
The screenshots of the digital oscilloscope of the different composite multiscroll attractors with differentR3. (a) Composite 6-scroll with R3 = 2.86 kΩ; (b) composite 4-scroll with R3 = 3.3 kΩ; (c) double-scroll with R3 = 0.5 kΩ; (d) double-scroll with R3 = 4.348 kΩ; (e) composite 3-scroll with R3 = 3.362 kΩ; (f) single-scroll with R3 = 0.263 kΩ.
(a)(b)(c)(d)(e)(f)


## 5. Conclusion

In this paper, a memristor-based chaotic system with abundant dynamic behaviors is reported. Compared with other memristor-based chaotic systems, the new system can display composite one- to six-scroll hidden attractors by changing only one system parameter. And the composite multiscroll attractors are composed of different types of attractors. For example, a composite six-scroll attractor is composed of a Chua’s double-scroll attractor and two jerk double-scroll attractors. The coexistence of symmetric and asymmetric attractors with different numbers of scrolls also reflects the particularity of the system. Furthermore, extreme multistability phenomenon and chaotic transient transition are investigated in the system. Finally, the feasibility of the proposed chaotic system is verified by hardware circuit, and the results match very well with the numerical simulations. Due to the fact that such a new chaotic system constructed in this paper has complex dynamic properties, it could be utilized advantageously in chaos-based engineering applications including image encryption, random bit generation, and chaos-based secure communication, and for future investigation.



---

*Source: 3259468-2020-11-20.xml*