# Coexistence of Multiple Points, Limit Cycles, and Strange Attractors in a Simple Autonomous Hyperjerk Circuit with Hyperbolic Sine Function

Authors: M. Fouodji Tsotsop; J. Kengne; G. Kenne; Z. Tabekoueng Njitacke
Journal: Complexity
 (2020)

---

## Abstract

In this contribution, a new elegant hyperjerk system with three equilibria and hyperbolic sine nonlinearity is investigated. In contrast to other models of hyperjerk systems where either hidden or self-excited attractors are obtained, the case reported in this work represents a unique one which displays the coexistence of self-excited chaotic attractors and stable fixed points. The dynamic properties of the new system are explored in terms of equilibrium point analyses, symmetry and dissipation, and existence of attractors as well. Common analysis tools (i.e., bifurcation diagram, Lyapunov exponents, and phase portraits) are used to highlight some important phenomena such as period-doubling bifurcation, chaos, periodic windows, and symmetric restoring crises. More interestingly, the system under consideration shows the coexistence of several types of stable states, including the coexistence of two, three, four, six, eight, and ten coexisting attractors. In addition, the system is shown to display antimonotonicity and offset boosting. Laboratory experimental measurements show a very good coherence with the theoretical predictions.


---

## Body

## 1. Introduction

Most authors have been interested in chaotic systems because of their sensitivity to the initial conditions and also to the variation of system parameters. Since the discovery of this phenomenon by Lorenz [1], many classical chaotic systems have emerged. We can mention the Rossler system [2], Chen system [3], Jafari system [4], Pham system [5], and Lü system just to name a few [6]. In the last few years, special attention has been given to “jerk systems” because of their simplicity and complex dynamics [7–12]. From a mathematical point of view, a generalization of the jerk dynamics is usually given in the following form:(1)dnxdtn=fdn−1xdtn−1,…,dxdt,x.Whenn=3, we have d3x/dt3=fd2x/dt2,dx/dt,x, which is called “jerk system” [13]. For n≥4, (1) turns to “hyperjerk system” or “snap system” [14]. In the literature, several authors have studied the latter. Generally, these systems exhibit multistability phenomenon which are the coexistence of multiple attractors solely depending on the initial conditions. These attractors are generally classified into two categories, namely, self-excited and hidden attractors [15–20]. Remember that self-excited attractors exist in systems with unstable equilibrium points [21–23]. In contrast, hidden attractors are characterized by the systems with no equilibrium [24–29], either by a line or a curve of equilibrium points, or system with stable equilibrium points [24, 25]. In addition, hidden attractors have a basin of attraction which does not intersect with the neighborhoods of equilibria.Interested by the self-excited attractors, many authors applied different techniques to hyperjerk systems. Some of these authors introduced different types of nonlinearities. For instance, in 2006, Cklouverakis and Sprott [22] presented a numerical study of a simple subclass of hyperjerk systems and showed that the 4th and 5th order hyperjerk systems developed some simple chaotic behaviors. In 2015, Sundiarapandian and coworkers [23] presented a new hyperchaotic 4-D hyperjerk system by adding a quadratic nonlinearity to the hyperjerk system of Chlouverakis–Sprott system. The authors present some qualitative and quantitative analyses of the new system. In 2017, Daltzis et al. [13] introduced a new hyperjerk system with two nonlinearities (absolute value and quintic term) and showed that the new system can develop hyperchaotic behaviors. Recently, Leutcho et al. [21] presented a new hyperjerk circuit with hyperbolic sine function and demonstrated that the novel proposed system is the unique one which is capable to exhibit the coexistence of nine periodic and chaotic attractors.Motivated by the above mentioned results, we present a new hyperjerk system with nonlinear position feedback involving a hyperbolic sine function. Our circuit is derived from the hyperjerk system proposed by Dalkirian and Sprott [7] by replacing the exponential nonlinearity by the hyperbolic sine function. The striking aspect of the proposed system is its ability to develop the coexistence of up to ten disconnected attractors including periodic, chaotic, and point attractors. The objectives of this work are as follows: (a) to present an analytical study of the proposed hyperjerk system; (b) to highlight regions in which we observe the coexistence of multiple attractors; (c) to point out some striking features like antimonotonicity and offset boosting; and (d) to verify the feasibility of the proposed model through an experimental study.This research is organized as follows. Section2 deals with the modeling process. The electronic consuration of the hyperjerk circuit is presented and the suitable mathematical model is derived to describe the dynamics of the novel hyperjerk, wherein some basic properties of the model are equally presented. In Section 3, the bifurcation structures of the system are investigated numerically. Also, in this section, some tools are used to show multistability observed in the novel system. Section 4 contains experimental study, and at the end of this section, it appears that coherence is observed between the theoretical and experimental analysis. Finally, Section 5 presents conclusion.


## 2. Description and Analysis of the Model

### 2.1. Circuit Description

It is important to know that the new circuit proposed here derives from the hyperjerk system proposed by Dalkirian and Sprott [7]. It is obtained by substituting the exponential nonlinearity by the hyperbolic sine function. Figure 1 represents the schematic diagram of the novel hyperjerk circuit. The circuit consists of four successive integrators associated to several feedback loops. In addition, the nonlinear feedback loop linked with the pair of semiconductor diodes D1,D2 is applied to the first integrator. The symmetrical nature [30] of the system is due to the antiparallel configuration of the diodes. In such type of configuration, the voltage across each diode is equal to the voltage of the resulting two-terminal device, while the current is the addition of the currents flowing through each diode. The symmetrical property of the nonlinearity is necessary for the occurrence of symmetric attractors [30]. We would like to recall that the pair of semiconductor diodes is the only nonlinear element responsible for the chaotic behavior displayed by the whole electronic circuit.Figure 1
Electronic circuit realization of the novel 4-D hyperjerk system with hyperbolic sine nonlinearity. Its simplicity is remarkable. The pair of semiconductor diodes implements the hyperbolic nonlinearity of the model.


### 2.2. State Equations

The following assumptions will be adopted throughout our analysis. Firstly, we considered that capacitors and operational amplifiers are ideal with the latter operating in linear domains. Secondly, the current-voltage characteristic (3) of the pair of semiconductor diodes D1andD2 is obtained from the Shockley diode equation [31, 32] as follows:(2)Id=ID1−ID2=ISexpVd/ηVT−1−ISexp−Vd/ηVT−1=2ISsinhVd/ηVT,where IS, VT=kbT/q, kb, T, q, and η1<η<2 are the intrinsic parameters of the diodes. By applying Kirchhoff’s laws to Figure 1 and considering the above assumptions, it can be shown that the voltages V1, V2, V3, and V4 satisfy the following set of four coupled first-order nonlinear differential equations:(3)C1dV1dt=V2R,C2dV2dt=V3Rm,C3dV3dt=V4Rd,C4dV4dt=V1Rc−V2Rb−V3Re−V4Ra−Id.Applying the following change of variables:(4)t=τRC,Vref=10ηVT,xjVref=Vjj=1,2,3,4,a=R/Ra,b=R/Rb,c=R/Rcd=R/Rd,m=R/Rm,e=R/Re,γ=2RIS/Vref,we get the normalized circuit equations which are expressed by the following smooth nonlinear fourth-order differential equations easy for numerical integration:(5)x˙1=x2,x˙2=mx3,x˙3=dx4,x˙4=cx1−bx2−ex3−ax4−γsinhx1,where the dot represents differentiation concerning the dimensionless time τ. Note that the nonlinear function only depends on the state variable x1 in system (5). γ will be kept constant throughout the numerical analysis: γ=0.0011. Therefore, during the bifurcation analysis of the 4-D system, c is considered like the control parameter (i.e., with respect to Rc). The values of electronic components used for both the numerical and experimental analyses are listed in Table 1. System (5) can be expressed equivalently in the general hyperjerk form as follows:(6)x⃜1=mbcx1−dγsinhx1−mbdx˙1−dex¨1−ax⃛1.Table 1
List of electronic components used during the analyses.
ComponentsValuesCase ACase BRa300Ω−1.7kΩ5.555kΩRb1.666kΩ2.631kΩRc10kΩ4kΩ−6kΩRd2kΩ7.407kΩRe5kΩ673.4ΩRm10kΩ4.545kΩR,R0,R0′′(numerical study)10kΩ10kΩR, R0 (experimental study)10kΩ10kΩR0'′ (experimental study)100kΩ100kΩC1,C2,C3,C410nF10nFD1, D21N41481N4148By observing equation (6), it can be noticed that our model belongs to the wider class of “elegant” hyperjerk dynamical systems defined in [14]. More interestingly, our model (5) represents one of the simplest autonomous 4-D systems reported recently, displaying the coexistence of up to ten fixed points, periodic and chaotic attractors.


### 2.3. Symmetry, Dissipation, and Existence of Attractors

Equation (5) being invariant following the transformation x1,x2,x3,x4⟺−x1,−x2,−x3,−x4, we can conclude that we will have a couple of solutions for a given parameter range. So, if x1,x2,x3,x4 is a solution of our system, then its symmetry −x1,−x2,−x3,−x4 will also be a solution. All these makes it possible to highlight the symmetrical nature of our system. In order to verify the dissipation property of our system, it is necessary to calculate the volume contraction rate (Λ=V−1dV/dt). For every point of space x1,x2,x3,x4T [33, 34], it is given by the following expression:(7)Λ=∂x˙1∂x1+∂x˙2∂x2+∂x˙3∂x3+∂x˙4∂x4=−a<0.The above expression is negative and does not depend on the space coordinates of the system, and thus we can conclude that the introduced system is dissipative.


### 2.4. Fixed Point Analysis

By canceling the right side of equation (5), it is possible to determine the equilibrium points of the system which play a crucial role in the study of the system dynamics. The resolution of equation (8) permits to obtain different equilibrium points of the system.(8)x2=0,x3=0,x4=0,cx1−bx2−ex3−ax4−γsinhx1=0.Note that the pointE00,0,0,0 is a trivial equilibrium point, while E1 and E2 are the solutions of the transcendental equation:(9)cx1−γsinhx1=0.By fixingc=2.442 and maintaining γ at the same previous value, we obtain the other nontrivial equilibrium points E1,E2=±10.79,0,0,0. The stability of the system in the neighborhood of those equilibrium points is studied by resolving the characteristic equation (12), resulting from the below equation:(10)detMJ−λId=0,where Id represents the 4×4 identity matrix and MJ the Jacobian matrix defined as follows:(11)MJ=010000m0000dc−γcoshx10−b−e−a.The characteristic equation obtained isλ4+c3λ3+c2λ2+c1λ+c0=0, with(12)c0=−mdc−γcoshx10,c1=mdb,c2=ed,c3=a.By applying the Lyapunov stability theory, it is shown that at the equilibrium pointE00,0,0,0, c−γ>0⟹c0<0, and thus the equilibrium E0 is unstable since the characteristic equation has coefficients with different signs. In contrast, the stability of the nonzero equilibrium points (E1,2±10.79,0,0,0) depends on the control parameters m and b. By applying the Routh–Hurwitz criterion [32, 33], we have shown that for b<bc1=3.4397 and b>bc2=5.56, the nontrivial equilibrium points are unstable, but for bc1<b<bc2, they are stable. These critical values are obtained from the calculation of the Hopf bifurcation conditions:(13)ωHopf=mdb/a,bc1,c2=12eam±Δ,with Δ=ea/m2+4a2/mdc−cosh10.79.Equation (13) provides the frequency of stable oscillations as well as the critical values of bc1,c2 corresponding to the Hopf bifurcation of the system. From Table 2, it follows that in the regime of (periodic or chaotic) oscillations, the three equilibria are unstable, and thus the system generates self-excited oscillations. For the following parameters a=1.8, d=1.35, e=15, m=3, and c=2.442, the equilibrium point E00,0,0,0 remains unstable for all values of control parameter b. Moreover, for some values of bifurcation parameter b, the nontrivial equilibria have pure imaginary roots, and thus the system presents the Hopf bifurcation. In order to verify the existence of the Hopf bifurcation in the system, eigenvalue locus is plotted. It shows the existence of Hopf bifurcation in the system which is characterized by the intersection of the eigenvalue locus with the imaginary axis. By observing Figures 2(a) and 2(b), we can certify that the new hyperjerk system presents Hopf bifurcation.Table 2
Corresponding eigenvalues of each equilibrium point according to the bifurcation parameterb.
Values of bifurcation parameterbEigenvalues at nontrivial fixedE1,E2, λ1,λ2,λ3,λ4Eigenvalues at the originE00,0,0,0λ1,λ2,λ3,λ43.1−0.9641±3.4906i0.0641±2.7165i (unstable)−0.5719±4.4347i−1.10410.4478 (unstable)3.2−0.9461±3.4712i0.0461±2.7346i (unstable)−0.5618±4.4350i−1.11850.4422 (unstable)3.4397−0.9000±3.4206i±2.7820i (neutral, Hopf bifurcation)−0.5378±4.4358i−1.15360.4292 (unstable)3.5−0.0023±3.5339i−0.8977±2.6358i (stable)−0.5318±4.4361i−1.16250.4260 (unstable)5.55−0.0093±3.5246i−0.8907±2.6459i (stable)−0.3297±4.4591i−1.47570.3351 (unstable)5.56±3.5370i−0.9000±2.6324i (neutral, Hopf bifurcation)−0.3287±4.4593i−1.47720.3347 (unstable)5.59−0.9476±3.2070i0.0476±3.0071i (unstable)−0.3258±4.4598i−1.48190.3336 (unstable)9.120.3784±4.0088i−1.2784±2.0827i (unstable)−0.0077±4.5497i−2.02100.2363 (unstable)Figure 2
Representation of eigenvalue locus in the complex planReλ,Imλ with the following parameter values: (a) a=1.8, b=3.8, c=2.442, d=1.35, e=14.85, and 2≤m≤10; (b) a=1.8, b=3.8, d=1.35, e=14.85, 2≤m≤10, and 2.34≤c≤12.985. The appearance of eigenvalues in complex conjugate pairs justifies the symmetry observed in the system and intersection of the curve with imaginary axis shows the presence of the Hopf bifurcation in the system.
(a)(b)


## 2.1. Circuit Description

It is important to know that the new circuit proposed here derives from the hyperjerk system proposed by Dalkirian and Sprott [7]. It is obtained by substituting the exponential nonlinearity by the hyperbolic sine function. Figure 1 represents the schematic diagram of the novel hyperjerk circuit. The circuit consists of four successive integrators associated to several feedback loops. In addition, the nonlinear feedback loop linked with the pair of semiconductor diodes D1,D2 is applied to the first integrator. The symmetrical nature [30] of the system is due to the antiparallel configuration of the diodes. In such type of configuration, the voltage across each diode is equal to the voltage of the resulting two-terminal device, while the current is the addition of the currents flowing through each diode. The symmetrical property of the nonlinearity is necessary for the occurrence of symmetric attractors [30]. We would like to recall that the pair of semiconductor diodes is the only nonlinear element responsible for the chaotic behavior displayed by the whole electronic circuit.Figure 1
Electronic circuit realization of the novel 4-D hyperjerk system with hyperbolic sine nonlinearity. Its simplicity is remarkable. The pair of semiconductor diodes implements the hyperbolic nonlinearity of the model.


## 2.2. State Equations

The following assumptions will be adopted throughout our analysis. Firstly, we considered that capacitors and operational amplifiers are ideal with the latter operating in linear domains. Secondly, the current-voltage characteristic (3) of the pair of semiconductor diodes D1andD2 is obtained from the Shockley diode equation [31, 32] as follows:(2)Id=ID1−ID2=ISexpVd/ηVT−1−ISexp−Vd/ηVT−1=2ISsinhVd/ηVT,where IS, VT=kbT/q, kb, T, q, and η1<η<2 are the intrinsic parameters of the diodes. By applying Kirchhoff’s laws to Figure 1 and considering the above assumptions, it can be shown that the voltages V1, V2, V3, and V4 satisfy the following set of four coupled first-order nonlinear differential equations:(3)C1dV1dt=V2R,C2dV2dt=V3Rm,C3dV3dt=V4Rd,C4dV4dt=V1Rc−V2Rb−V3Re−V4Ra−Id.Applying the following change of variables:(4)t=τRC,Vref=10ηVT,xjVref=Vjj=1,2,3,4,a=R/Ra,b=R/Rb,c=R/Rcd=R/Rd,m=R/Rm,e=R/Re,γ=2RIS/Vref,we get the normalized circuit equations which are expressed by the following smooth nonlinear fourth-order differential equations easy for numerical integration:(5)x˙1=x2,x˙2=mx3,x˙3=dx4,x˙4=cx1−bx2−ex3−ax4−γsinhx1,where the dot represents differentiation concerning the dimensionless time τ. Note that the nonlinear function only depends on the state variable x1 in system (5). γ will be kept constant throughout the numerical analysis: γ=0.0011. Therefore, during the bifurcation analysis of the 4-D system, c is considered like the control parameter (i.e., with respect to Rc). The values of electronic components used for both the numerical and experimental analyses are listed in Table 1. System (5) can be expressed equivalently in the general hyperjerk form as follows:(6)x⃜1=mbcx1−dγsinhx1−mbdx˙1−dex¨1−ax⃛1.Table 1
List of electronic components used during the analyses.
ComponentsValuesCase ACase BRa300Ω−1.7kΩ5.555kΩRb1.666kΩ2.631kΩRc10kΩ4kΩ−6kΩRd2kΩ7.407kΩRe5kΩ673.4ΩRm10kΩ4.545kΩR,R0,R0′′(numerical study)10kΩ10kΩR, R0 (experimental study)10kΩ10kΩR0'′ (experimental study)100kΩ100kΩC1,C2,C3,C410nF10nFD1, D21N41481N4148By observing equation (6), it can be noticed that our model belongs to the wider class of “elegant” hyperjerk dynamical systems defined in [14]. More interestingly, our model (5) represents one of the simplest autonomous 4-D systems reported recently, displaying the coexistence of up to ten fixed points, periodic and chaotic attractors.


## 2.3. Symmetry, Dissipation, and Existence of Attractors

Equation (5) being invariant following the transformation x1,x2,x3,x4⟺−x1,−x2,−x3,−x4, we can conclude that we will have a couple of solutions for a given parameter range. So, if x1,x2,x3,x4 is a solution of our system, then its symmetry −x1,−x2,−x3,−x4 will also be a solution. All these makes it possible to highlight the symmetrical nature of our system. In order to verify the dissipation property of our system, it is necessary to calculate the volume contraction rate (Λ=V−1dV/dt). For every point of space x1,x2,x3,x4T [33, 34], it is given by the following expression:(7)Λ=∂x˙1∂x1+∂x˙2∂x2+∂x˙3∂x3+∂x˙4∂x4=−a<0.The above expression is negative and does not depend on the space coordinates of the system, and thus we can conclude that the introduced system is dissipative.


## 2.4. Fixed Point Analysis

By canceling the right side of equation (5), it is possible to determine the equilibrium points of the system which play a crucial role in the study of the system dynamics. The resolution of equation (8) permits to obtain different equilibrium points of the system.(8)x2=0,x3=0,x4=0,cx1−bx2−ex3−ax4−γsinhx1=0.Note that the pointE00,0,0,0 is a trivial equilibrium point, while E1 and E2 are the solutions of the transcendental equation:(9)cx1−γsinhx1=0.By fixingc=2.442 and maintaining γ at the same previous value, we obtain the other nontrivial equilibrium points E1,E2=±10.79,0,0,0. The stability of the system in the neighborhood of those equilibrium points is studied by resolving the characteristic equation (12), resulting from the below equation:(10)detMJ−λId=0,where Id represents the 4×4 identity matrix and MJ the Jacobian matrix defined as follows:(11)MJ=010000m0000dc−γcoshx10−b−e−a.The characteristic equation obtained isλ4+c3λ3+c2λ2+c1λ+c0=0, with(12)c0=−mdc−γcoshx10,c1=mdb,c2=ed,c3=a.By applying the Lyapunov stability theory, it is shown that at the equilibrium pointE00,0,0,0, c−γ>0⟹c0<0, and thus the equilibrium E0 is unstable since the characteristic equation has coefficients with different signs. In contrast, the stability of the nonzero equilibrium points (E1,2±10.79,0,0,0) depends on the control parameters m and b. By applying the Routh–Hurwitz criterion [32, 33], we have shown that for b<bc1=3.4397 and b>bc2=5.56, the nontrivial equilibrium points are unstable, but for bc1<b<bc2, they are stable. These critical values are obtained from the calculation of the Hopf bifurcation conditions:(13)ωHopf=mdb/a,bc1,c2=12eam±Δ,with Δ=ea/m2+4a2/mdc−cosh10.79.Equation (13) provides the frequency of stable oscillations as well as the critical values of bc1,c2 corresponding to the Hopf bifurcation of the system. From Table 2, it follows that in the regime of (periodic or chaotic) oscillations, the three equilibria are unstable, and thus the system generates self-excited oscillations. For the following parameters a=1.8, d=1.35, e=15, m=3, and c=2.442, the equilibrium point E00,0,0,0 remains unstable for all values of control parameter b. Moreover, for some values of bifurcation parameter b, the nontrivial equilibria have pure imaginary roots, and thus the system presents the Hopf bifurcation. In order to verify the existence of the Hopf bifurcation in the system, eigenvalue locus is plotted. It shows the existence of Hopf bifurcation in the system which is characterized by the intersection of the eigenvalue locus with the imaginary axis. By observing Figures 2(a) and 2(b), we can certify that the new hyperjerk system presents Hopf bifurcation.Table 2
Corresponding eigenvalues of each equilibrium point according to the bifurcation parameterb.
Values of bifurcation parameterbEigenvalues at nontrivial fixedE1,E2, λ1,λ2,λ3,λ4Eigenvalues at the originE00,0,0,0λ1,λ2,λ3,λ43.1−0.9641±3.4906i0.0641±2.7165i (unstable)−0.5719±4.4347i−1.10410.4478 (unstable)3.2−0.9461±3.4712i0.0461±2.7346i (unstable)−0.5618±4.4350i−1.11850.4422 (unstable)3.4397−0.9000±3.4206i±2.7820i (neutral, Hopf bifurcation)−0.5378±4.4358i−1.15360.4292 (unstable)3.5−0.0023±3.5339i−0.8977±2.6358i (stable)−0.5318±4.4361i−1.16250.4260 (unstable)5.55−0.0093±3.5246i−0.8907±2.6459i (stable)−0.3297±4.4591i−1.47570.3351 (unstable)5.56±3.5370i−0.9000±2.6324i (neutral, Hopf bifurcation)−0.3287±4.4593i−1.47720.3347 (unstable)5.59−0.9476±3.2070i0.0476±3.0071i (unstable)−0.3258±4.4598i−1.48190.3336 (unstable)9.120.3784±4.0088i−1.2784±2.0827i (unstable)−0.0077±4.5497i−2.02100.2363 (unstable)Figure 2
Representation of eigenvalue locus in the complex planReλ,Imλ with the following parameter values: (a) a=1.8, b=3.8, c=2.442, d=1.35, e=14.85, and 2≤m≤10; (b) a=1.8, b=3.8, d=1.35, e=14.85, 2≤m≤10, and 2.34≤c≤12.985. The appearance of eigenvalues in complex conjugate pairs justifies the symmetry observed in the system and intersection of the curve with imaginary axis shows the presence of the Hopf bifurcation in the system.
(a)(b)


## 3. Numerical Computation

### 3.1. Numerical Techniques

System (5) is resolved numerically in order to highlight the rich variety of bifurcation that can be observed in a new hyperjerk system. The dynamic properties of the model were numerically simulated in Turbo Pascal using the fourth-order Runge–Kutta method with a constant time step size of 2×10−3, and parameters are taken in extended precision mode. The transient phase is canceled by integrating the system for a long time. The bifurcation diagram and the Lyapunov exponent are the traditional tools that measure the dependence of the system on the initial conditions as well as the sequence that leads to chaos in the system. The algorithms of Wolf and his collaborators [35] are used for calculating the Lyapunov exponents.


### 3.2. Bifurcation, Chaos in a Novel Hyperjerk Circuit

Different scenarios exhibited by the proposed hyperjerk system are obtained by plotting the bifurcation diagrams. The bifurcation diagram of Figure3(a) is obtained by plotting the local maxima of the variable x1 according to the bifurcation parameter a, the other parameters being fixed at c=1, b=6, d=5, e=2, and m=1. It can be noted that it is a period-doubling route to chaos because the transition from period-1 attractor to double-band chaos is as follows: period-1 ⟶ period-2 ⟶ period-4 ⟶ single-band chaos ⟶ period-5 ⟶ single-band chaos ⟶ double-band chaos. Figure 4 clearly shows the above transition. It is obtained by progressively varying the control parameter. The exact nature of the attractors mentioned above is defined by the graphs of the four largest Lyapunov exponents shown in Figure 3(b). We can observe in Figure 3(b) that periodic attractors are characterized by λ1=0, λ2, λ3, and λ4<0, while chaotic attractors have the following characteristics: λ1>0, λ2=0, λ3, and λ4<0. A perfect coherence is observed between the bifurcation diagram and the corresponding graphs of the four largest Lyapunov exponents. In order to show the complexity of the new hyperjerk circuit, the chaotic attractor has been projected on several planes (Figures 5(a)–5(f)), as well as the Poincaré section (Figure 5(g)). We can observe that the double-band chaos completely changes when moving from one plane to another. For the value of the bifurcation parameter a=17.04, the coexistence of four periodic and chaotic attractors is observed in the novel proposed system (see Figure 6). In order to illustrate the Hopf bifurcation previously proved by theoretical calculations, the bifurcation diagram of Figure 7 has been represented. Stable state is characterized by a fixed point with λmax<0, while oscillatory state is characterized by λmax≥0.Figure 3
(a) Bifurcation diagram showing local maxima of the coordinatex1 versus c and (b, c) the corresponding graphs of four largest Lyapunov exponents plotted in the range 6≤a≤30, with c=1, b=6, d=5, e=2, and m=1 and initial conditions x10,x20,x30,x40=0,0,2.4,0.Figure 4
Numerical phase space trajectories showing routes to chaos in the system when varying the control parameterc: (a) period-1 for a=7, (b) period-2 for a=12.38, (c) period-4 for a=13, (d) period-8 for a=13.49, (e) single-band chaos for a=14.6, (f) period-5 for a=16.08, (g) single -band chaos for a=16.7, (h) single-band chaos for a=20, and (h) double-band chaos for a=30. Initial conditions x10,x20,x30,x40 are 0,0,2.4,0. The others parameters are fixed as follows: c=1, b=6, d=5, e=2, and m=1.
(a)(b)(c)(d)(e)(f)(g)(h)(i)Figure 5
2-D projection of the phase portraits of symmetric double-band chaotic attractors (a-f) of system (3) plotted into planes (x1−x4), (x2−x3), (x2−x4), (x3−x4) and corresponding double-sided Poincaré section (e) in the plane x1=0. Parameters are the same as those in Figure 4.
(a)(b)(c)(d)(e)(f)(g)Figure 6
Coexistence of four attractors (a pair of period-11 limit cycle and a pair of chaotic attractors) witha=17.04, and their corresponding initial conditions are ±4.8,0,0,0 and ±1.2,0,0,0.
(a)(b)(c)(d)Figure 7
(a) Bifurcation diagram showing local maxima of the coordinatex1 versus b and (b) the corresponding graph of largest Lyapunov exponent (λmax) plotted by decreasing b in the range 3.5≤b≤9.12, with a=1.8, m=3, c=2.442, d=1.35, and e=15 and initial conditions x10,x20,x30,x40=0,0,9.333,0.
(a)(b)


### 3.3. Multistability

In this section, we demonstrate the variety of dynamical regimes in the new 4-D system. We show that depending on the values of the system parameters, the system exhibits very rich dynamics and bifurcation scenarios. A multistable system is a system with various coexisting stable states (chaotic, point, and periodic state) under the same system parameters, with different initial conditions. In recent years, the phenomenon of multistability phenomenon has been reported in many nonlinear dynamic systems [13, 36–46].


#### 3.3.1. Coexistence of Attractors with respect to Bifurcation Parameterc

By changing the system parameters and consideringc as bifurcation parameter, we observe a completely different behavior. In addition, a very interesting phenomenon which is the coexistence of multiple attractors appears in the new 4-D hyperjerk. For this phenomenon to be illustrated, the bifurcation diagrams of Figure 8 are plotted using the following method:(i)
The blue diagram is obtained by simultaneously increasing the value of the control parameterc as well as the initial condition x0. At each iteration, we assign to x0 the new value of the control parameter c.(ii)
The red diagram is obtained by incrementingc from its minimum value 2.34 to its maximum value 2.985, with a carefully chosen step. Note that the solutions of the system at each iteration are considered as the initial condition of the next iteration.(iii)
The cyan diagram respects the previous procedure, with the initial condition−10.67,0,0,0, whereas the black diagram follows the same procedure as previously described but the only difference is the decrease of the control parameter c. The initial condition is 10,0,0,0.(iv)
The magenta diagram is obtained by increasing the control parameterc from 2.398 to 2.985, followed by decreasing the bifurcation parameter c from 2.398 to 2.34.Figure 8
symmetric coexisting bifurcation diagrams showing complex structure of the new hyperjerk fora=1.8, b=3.8, d=1.35, e=14.85, and m=2.2. Initial conditions y0=z0=w0=0 and x0 are fixed as follows for obtaining these diagrams: red diagrams for increasing c with x0=5.25, black diagram for decreasing c with x0=−10, cyan diagrams for increasing c with x0=−10.67, magenta diagram for increasing and decreasing, started at c=2.5442 with x0=−4, and blue diagram for starting c with same initial condition x0=c. For more information, see Table 3.
(a)(b)We can observe in Figure8 several windows of coexisting attractors. For more details about the methods used to plot the bifurcation diagrams of Figure 8, see Table 3. The enlarged bifurcation diagram of Figure 9 shows the hysteretic domain, plotted in the range 2.442≤c≤2.488, and the techniques used to plot the diagrams are also presented in Table 3. Figure 9 shows the coexistence of six and eight different limit cycles, chaotic and point attractors. Some sample phase portraits showing the coexistence of six and eight attractors are presented in Figures 10 and 11, respectively. Some basins of attractions showing the initial conditions domains of the coexisting attractors are presented in Figure 12. The coexistence of four attractors is clearly denoted (a pair of periodic attractors (black and yellow) and a pair of chaotic attractors (blue and green)). Note that there is a perfect symmetry between the different cross sections of the competing attractors.Table 3
Techniques used to obtain coexisting bifurcation diagrams and corresponding initial conditions.
Fig.n0Color graphParameter rangeSweeping directionInitial conditionx10,x20,x30,x40Figure8Blue2.34≤c≤2.985Upwardc,0,0,0Red2.34≤c≤2.985Upward5.25,0,0,0Cyan2.34≤c≤2.985Upward−10.67,0,0,0Black2.985≤c≤2.34Downward−10,0,0,0Magenta2.398≤c≤2.985Upward−4,0,0,02.398≤c≤2.34Downward−4,0,0,0Figure9Red2.442≤c≤2.488Upward5,0,0,0Black2.488≤c≤2.442Downward5,0,0,0Cyan2.442≤c≤2.488Upward−10.67,0,0,0Magenta2.442≤c≤2.488Upward12,0,0,0Figure13Blue2.802≤m≤2.865Upward−5.5,0,0,02.738≤m≤2.802DownwardRed2.802≤m≤2.865Upward12.8,0,0,02.738≤m≤2.802DownwardMagenta2.802≤m≤2.865Upward−10.1,0,0,02.738≤m≤2.802DownwardFigure 9
Enlargement of bifurcation diagram of Figure8 plotted in the range 2.442≤c≤2.488 showing the region in which the system exhibits multiple coexisting attractors. The system parameters are the same as the one of Figure 8.Figure 10
Bifurcation like sequence showing local maxima of the coordinatex3 versus the initial condition . plotted in the range 0≤x0≤30, with a=1.8, b=3.8, c=2.5442, d=1.35, and e=14.85. It can be observed the coexistence of ten periodic, chaotic, and point attractors.Figure 11
(a, b) Bifurcation diagrams showing local maxima of the coordinatex1 versus m and (c) the corresponding graph of largest Lyapunov exponent (λmax) plotted in the range 2≤m≤4.18, with a=1.8, b=3.8, c=2.442, and d=1.35.
(a)(b)(c)Figure 12
Enlargement of the bifurcation diagram of Figure11(a) showing the region in which the system exhibits multiple coexisting attractors. This region corresponds to values of m in the range: 2.735≤m≤2.865. Three sets of data are superimposed. For more information about the methods used to plot, see Table 3.The initial conditions of the coexistence of ten attractors exhibited by the proposed hyperjerk system are presented in Figure13. The bifurcation like sequence of Figure 13 shows the variation of x30 in terms of the control parameter c, and the other initial conditions are set to zero (x10=x20=x40=0). Note that chaotic attractors are characterized by an unlimited number of points, while periodic attractors are characterized by a finite number of points. In the same line, the basin of attraction shows the different domains of convergence of similar attractors. The phase portrait of Figure 14 illustrates the coexistence of ten attractors. The initials conditions of the coexisting attractors are given in Table 4.Figure 13
Coexistence of six different attractors (a pair of period-2 limit cycle, a pair of chaotic attractors, and a pair of fixed point attractor) for the following values of system parameters:a=1.8, b=3.8, c=2.442, d=1.35, e=14.85, and m=2.802. Initial conditions x10,x20,x30,x40 are, respectively, ±5.5,0,0,0, ±12.8,0,0,0, and ±10.1,0,0,0.
(a)(b)(c)(d)(e)Figure 14
(a, b) Coexistence of four distinct attractors (a pair of chaotic attractors represented inx1−x2 and x2−x3 plans), (c) a pair of fixed point attractor) and their corresponding basin of attraction. Initial conditions x10,x20,x30,x40 are, respectively,±7.15,0,±8.613,0 and ±12.8,0,0,0. The rest of the parameters area=1.8, b=3.5, c=2.442, d=1.35, e=15, and m=3.
(a)(b)(c)(d)Figure 15
Cross sections of basin of attraction forx30=x40=0, x20=x40=0, x20=x30=0, and x10=x40=0 corresponding to the asymmetric pair of chaotic (blue and green) and period-2 attractors (yellow and black) obtained for c=2.44. Red regions correspond to unbounded motion.
(a)(b)(c)(d)Table 4
Details of the coexistences observed in the novel hyperjerk system.
Fig.n0Type of coexistencesValues of control parameterInitial conditionx10,x20,x30,x40Figure16A symmetric pair of period-2 attractors, a symmetric pair of chaotic attractors, and a pair of fixed point.c=2.44(a)0,0,±44.4,0(b)0,0,±48,0(c)0,0,±26.4,0Figure17A symmetric pair of period-2 attractors, a symmetric pair of chaotic attractors, a symmetric pair of period-9 attractors, and a pair of fixed point.c=2.454(a)0,0,±25.2,0(b)0,0,±48,0(c)0,0,±10.8,0(d)0,0,±26.4,0Figure18A symmetric pair of period-2 attractors, a fixed point, a symmetric pair of period-12 attractors, 2 symmetric pairs of chaotic attractors, and a pair of fixed point.c=2.5442(A”A1,A2“)0,0,±45.6,0(b)0,0,±18,0(c)0,0,±49.2,0(d)0,0,±48,0(e)0,0,±26.4,0


#### 3.3.2. Coexistence of Attractors with respect to Bifurcation Parametersm and b

In order to investigate the sensitivity of the new hyperjerk system in terms of the bifurcation parameterm, the other parameters are fixed as follows: a=1.8, b=3.8, c=2.442, d=1.35, and e=14.85. We found that the novel hyperjerk system can exhibit striking bifurcation sequences when varying the control parameter m in the range 2≤m≤4.18. With reference to Figure 11, the bifurcation diagram in black and the one in blue are obtained by increasing and decreasing the values of the parameter m, while the one in red is obtained by fixing the initial conditions at x10,x20,x30,x40=−4,0,0,0. A window of hysteretic dynamics can be identified in the range 2.6≤m≤2.9. The enlarged bifurcation diagram of Figure 12 clearly illustrates the domain of the coexistence of multiple attractors observed in the new hyperjerk system according to the bifurcation parameter m. Different methods used to plot these bifurcation diagrams are presented in Table 3. Up to six different periodic, chaotic, and point attractors can be obtained by only changing the initial conditions. For instance, sample phase portraits of the coexistence of six distinct attractors are presented in Figure 13.Figure 16
Coexistence of six different attractors (a pair of symmetric period-2 limit cycle, a pair of symmetric period-9 limit cycle, a pair of symmetric chaotic attractors, and a pair of point attractors). The values of others system parameters are fixed as follows:a=1.8, b=3.8, d=1.35, e=14.85, and m=2.2. Initial conditions are given in Table 4.
(a)(b)(c)Figure 17
Coexistence of eight different attractors (a pair of symmetric period-2 limit cycle, a pair of symmetric chaotic attractors, and a pair of point attractors). The parameters are the same as those in Figure16. Initial conditions are given in Table 4.
(a)(b)(c)(d)During the mathematical analyses, it has been shown that the Hopf bifurcation was depending on the control parameterm. The bifurcation diagram of Figure 13(a) clearly illustrates this phenomenon characterized by the stable state followed by the unstable state. Moreover, this control parameter also highlights the coexistence of multiple attractors exhibited by the new 4-D system. By considering the following sets of the parameters: a=1.8, b=3.5, c=2.442, d=1.35, e=15, and m=3, we discover that the new 4-D system displays the coexistence of four distinct chaotic and point attractors. The phase portraits of the Figure 14 and their corresponding cross section of the basin of attraction clearly show the coexistence phenomenon and also give the initial condition domain of each attractor. The green and black domains represent the initial conditions regions of the pair of chaotic attractors, while the yellow domain represents the initial condition regions of the pair of point attractors.Figure 18
Coexistence of ten different attractors (a pair of symmetric period-2 limit cycle, 2 pair of symmetric chaotic attractors, a pair of symmetric period-9 limit cycle, and a pair of point attractors). The parameters are the same as those in Figure16. Initial conditions are given in Table 4.
(a)(b)(c)(d)(e)


### 3.4. Offset Boosting Scenario

Another property of system (3) is the possibility to develop an offset boosting effect. In our model, x1 appears only in the fourth line of equation (3), and thus this variable is a bootable variable [47–52]. Assuming the transformation x1⟶x1+k where k is a constant, equation (5) can be rewritten accordingly as(14)x˙1=x2,x˙2=mx3,x˙3=dx4,x˙4=cx1+k−bx2−ex3−ax4−γsinhx1+k.Figure19 clearly presents offset boosting of the double-band chaotic attractor. The following values of parameter: k=0 (blue), k=10 (red) and k=−10 (green) are used to plot them in x1-x4 and x1-x3 planes.Figure 19
Offset boosting of the double-band chaotic attractor with the following values of parameter:k=0(blue), k=10(red), and k=−10(green). The other parameters are the same as those in Figure 4.
(a)(b)


### 3.5. Antimonotonicity

By decreasing the value of the control parametere15≤e≤20, we can observe the formation and destruction of periodic orbits via reverse period-doubling bifurcation sequences. This interesting phenomenon has been reported in the literature. It is reported in various nonlinear systems such as Duffing oscillator [12], Chua circuit [53], and second-order nonlinear nonautonomous circuit [54, 55]. This phenomenon was reported for the first time in the hyperjerk system by Leutcho et al. [21]. The creation of periodic seas in the parameter space is the necessary requirement for a nonlinear system to experience forward and reverse period-doubling cascade [21]. Sample illustrations are represented in Figure 20, where some bifurcation diagrams are shown. These diagrams are obtained for each discrete value of the control parameter c. In Figure 20, note that for c=2.5442, we have period-2 bubble and for a slight adjustment of the control parameter c, period-4 bubble is observed for c=2.8442, whereas for c=2.9, we have a period-8 bubble. In the same order, chaotic bubbles are formed for c=2.97,andc=2.99. The increase of control parameter c causes the creation of other bubbles, and it finally results in an infinite tree (like chaos).Figure 20
Bifurcation diagrams showing local maxima of the coordinatex1 of the attractor in Poincaré cross section in terms of the control parameter a (bubbling): (a) period-2 bubble for c=2.5442, (b) period-4 bubble for c=2.8442, (c) period-8 bubble for c=2.9, (d) single-band chaos bubble for c=2.93, (e) single-band chaos bubble for c=2.99, and (f) double-band chaos bubble for c=3.0.
(a)(b)(c)(d)(e)(f)


## 3.1. Numerical Techniques

System (5) is resolved numerically in order to highlight the rich variety of bifurcation that can be observed in a new hyperjerk system. The dynamic properties of the model were numerically simulated in Turbo Pascal using the fourth-order Runge–Kutta method with a constant time step size of 2×10−3, and parameters are taken in extended precision mode. The transient phase is canceled by integrating the system for a long time. The bifurcation diagram and the Lyapunov exponent are the traditional tools that measure the dependence of the system on the initial conditions as well as the sequence that leads to chaos in the system. The algorithms of Wolf and his collaborators [35] are used for calculating the Lyapunov exponents.


## 3.2. Bifurcation, Chaos in a Novel Hyperjerk Circuit

Different scenarios exhibited by the proposed hyperjerk system are obtained by plotting the bifurcation diagrams. The bifurcation diagram of Figure3(a) is obtained by plotting the local maxima of the variable x1 according to the bifurcation parameter a, the other parameters being fixed at c=1, b=6, d=5, e=2, and m=1. It can be noted that it is a period-doubling route to chaos because the transition from period-1 attractor to double-band chaos is as follows: period-1 ⟶ period-2 ⟶ period-4 ⟶ single-band chaos ⟶ period-5 ⟶ single-band chaos ⟶ double-band chaos. Figure 4 clearly shows the above transition. It is obtained by progressively varying the control parameter. The exact nature of the attractors mentioned above is defined by the graphs of the four largest Lyapunov exponents shown in Figure 3(b). We can observe in Figure 3(b) that periodic attractors are characterized by λ1=0, λ2, λ3, and λ4<0, while chaotic attractors have the following characteristics: λ1>0, λ2=0, λ3, and λ4<0. A perfect coherence is observed between the bifurcation diagram and the corresponding graphs of the four largest Lyapunov exponents. In order to show the complexity of the new hyperjerk circuit, the chaotic attractor has been projected on several planes (Figures 5(a)–5(f)), as well as the Poincaré section (Figure 5(g)). We can observe that the double-band chaos completely changes when moving from one plane to another. For the value of the bifurcation parameter a=17.04, the coexistence of four periodic and chaotic attractors is observed in the novel proposed system (see Figure 6). In order to illustrate the Hopf bifurcation previously proved by theoretical calculations, the bifurcation diagram of Figure 7 has been represented. Stable state is characterized by a fixed point with λmax<0, while oscillatory state is characterized by λmax≥0.Figure 3
(a) Bifurcation diagram showing local maxima of the coordinatex1 versus c and (b, c) the corresponding graphs of four largest Lyapunov exponents plotted in the range 6≤a≤30, with c=1, b=6, d=5, e=2, and m=1 and initial conditions x10,x20,x30,x40=0,0,2.4,0.Figure 4
Numerical phase space trajectories showing routes to chaos in the system when varying the control parameterc: (a) period-1 for a=7, (b) period-2 for a=12.38, (c) period-4 for a=13, (d) period-8 for a=13.49, (e) single-band chaos for a=14.6, (f) period-5 for a=16.08, (g) single -band chaos for a=16.7, (h) single-band chaos for a=20, and (h) double-band chaos for a=30. Initial conditions x10,x20,x30,x40 are 0,0,2.4,0. The others parameters are fixed as follows: c=1, b=6, d=5, e=2, and m=1.
(a)(b)(c)(d)(e)(f)(g)(h)(i)Figure 5
2-D projection of the phase portraits of symmetric double-band chaotic attractors (a-f) of system (3) plotted into planes (x1−x4), (x2−x3), (x2−x4), (x3−x4) and corresponding double-sided Poincaré section (e) in the plane x1=0. Parameters are the same as those in Figure 4.
(a)(b)(c)(d)(e)(f)(g)Figure 6
Coexistence of four attractors (a pair of period-11 limit cycle and a pair of chaotic attractors) witha=17.04, and their corresponding initial conditions are ±4.8,0,0,0 and ±1.2,0,0,0.
(a)(b)(c)(d)Figure 7
(a) Bifurcation diagram showing local maxima of the coordinatex1 versus b and (b) the corresponding graph of largest Lyapunov exponent (λmax) plotted by decreasing b in the range 3.5≤b≤9.12, with a=1.8, m=3, c=2.442, d=1.35, and e=15 and initial conditions x10,x20,x30,x40=0,0,9.333,0.
(a)(b)


## 3.3. Multistability

In this section, we demonstrate the variety of dynamical regimes in the new 4-D system. We show that depending on the values of the system parameters, the system exhibits very rich dynamics and bifurcation scenarios. A multistable system is a system with various coexisting stable states (chaotic, point, and periodic state) under the same system parameters, with different initial conditions. In recent years, the phenomenon of multistability phenomenon has been reported in many nonlinear dynamic systems [13, 36–46].


### 3.3.1. Coexistence of Attractors with respect to Bifurcation Parameterc

By changing the system parameters and consideringc as bifurcation parameter, we observe a completely different behavior. In addition, a very interesting phenomenon which is the coexistence of multiple attractors appears in the new 4-D hyperjerk. For this phenomenon to be illustrated, the bifurcation diagrams of Figure 8 are plotted using the following method:(i)
The blue diagram is obtained by simultaneously increasing the value of the control parameterc as well as the initial condition x0. At each iteration, we assign to x0 the new value of the control parameter c.(ii)
The red diagram is obtained by incrementingc from its minimum value 2.34 to its maximum value 2.985, with a carefully chosen step. Note that the solutions of the system at each iteration are considered as the initial condition of the next iteration.(iii)
The cyan diagram respects the previous procedure, with the initial condition−10.67,0,0,0, whereas the black diagram follows the same procedure as previously described but the only difference is the decrease of the control parameter c. The initial condition is 10,0,0,0.(iv)
The magenta diagram is obtained by increasing the control parameterc from 2.398 to 2.985, followed by decreasing the bifurcation parameter c from 2.398 to 2.34.Figure 8
symmetric coexisting bifurcation diagrams showing complex structure of the new hyperjerk fora=1.8, b=3.8, d=1.35, e=14.85, and m=2.2. Initial conditions y0=z0=w0=0 and x0 are fixed as follows for obtaining these diagrams: red diagrams for increasing c with x0=5.25, black diagram for decreasing c with x0=−10, cyan diagrams for increasing c with x0=−10.67, magenta diagram for increasing and decreasing, started at c=2.5442 with x0=−4, and blue diagram for starting c with same initial condition x0=c. For more information, see Table 3.
(a)(b)We can observe in Figure8 several windows of coexisting attractors. For more details about the methods used to plot the bifurcation diagrams of Figure 8, see Table 3. The enlarged bifurcation diagram of Figure 9 shows the hysteretic domain, plotted in the range 2.442≤c≤2.488, and the techniques used to plot the diagrams are also presented in Table 3. Figure 9 shows the coexistence of six and eight different limit cycles, chaotic and point attractors. Some sample phase portraits showing the coexistence of six and eight attractors are presented in Figures 10 and 11, respectively. Some basins of attractions showing the initial conditions domains of the coexisting attractors are presented in Figure 12. The coexistence of four attractors is clearly denoted (a pair of periodic attractors (black and yellow) and a pair of chaotic attractors (blue and green)). Note that there is a perfect symmetry between the different cross sections of the competing attractors.Table 3
Techniques used to obtain coexisting bifurcation diagrams and corresponding initial conditions.
Fig.n0Color graphParameter rangeSweeping directionInitial conditionx10,x20,x30,x40Figure8Blue2.34≤c≤2.985Upwardc,0,0,0Red2.34≤c≤2.985Upward5.25,0,0,0Cyan2.34≤c≤2.985Upward−10.67,0,0,0Black2.985≤c≤2.34Downward−10,0,0,0Magenta2.398≤c≤2.985Upward−4,0,0,02.398≤c≤2.34Downward−4,0,0,0Figure9Red2.442≤c≤2.488Upward5,0,0,0Black2.488≤c≤2.442Downward5,0,0,0Cyan2.442≤c≤2.488Upward−10.67,0,0,0Magenta2.442≤c≤2.488Upward12,0,0,0Figure13Blue2.802≤m≤2.865Upward−5.5,0,0,02.738≤m≤2.802DownwardRed2.802≤m≤2.865Upward12.8,0,0,02.738≤m≤2.802DownwardMagenta2.802≤m≤2.865Upward−10.1,0,0,02.738≤m≤2.802DownwardFigure 9
Enlargement of bifurcation diagram of Figure8 plotted in the range 2.442≤c≤2.488 showing the region in which the system exhibits multiple coexisting attractors. The system parameters are the same as the one of Figure 8.Figure 10
Bifurcation like sequence showing local maxima of the coordinatex3 versus the initial condition . plotted in the range 0≤x0≤30, with a=1.8, b=3.8, c=2.5442, d=1.35, and e=14.85. It can be observed the coexistence of ten periodic, chaotic, and point attractors.Figure 11
(a, b) Bifurcation diagrams showing local maxima of the coordinatex1 versus m and (c) the corresponding graph of largest Lyapunov exponent (λmax) plotted in the range 2≤m≤4.18, with a=1.8, b=3.8, c=2.442, and d=1.35.
(a)(b)(c)Figure 12
Enlargement of the bifurcation diagram of Figure11(a) showing the region in which the system exhibits multiple coexisting attractors. This region corresponds to values of m in the range: 2.735≤m≤2.865. Three sets of data are superimposed. For more information about the methods used to plot, see Table 3.The initial conditions of the coexistence of ten attractors exhibited by the proposed hyperjerk system are presented in Figure13. The bifurcation like sequence of Figure 13 shows the variation of x30 in terms of the control parameter c, and the other initial conditions are set to zero (x10=x20=x40=0). Note that chaotic attractors are characterized by an unlimited number of points, while periodic attractors are characterized by a finite number of points. In the same line, the basin of attraction shows the different domains of convergence of similar attractors. The phase portrait of Figure 14 illustrates the coexistence of ten attractors. The initials conditions of the coexisting attractors are given in Table 4.Figure 13
Coexistence of six different attractors (a pair of period-2 limit cycle, a pair of chaotic attractors, and a pair of fixed point attractor) for the following values of system parameters:a=1.8, b=3.8, c=2.442, d=1.35, e=14.85, and m=2.802. Initial conditions x10,x20,x30,x40 are, respectively, ±5.5,0,0,0, ±12.8,0,0,0, and ±10.1,0,0,0.
(a)(b)(c)(d)(e)Figure 14
(a, b) Coexistence of four distinct attractors (a pair of chaotic attractors represented inx1−x2 and x2−x3 plans), (c) a pair of fixed point attractor) and their corresponding basin of attraction. Initial conditions x10,x20,x30,x40 are, respectively,±7.15,0,±8.613,0 and ±12.8,0,0,0. The rest of the parameters area=1.8, b=3.5, c=2.442, d=1.35, e=15, and m=3.
(a)(b)(c)(d)Figure 15
Cross sections of basin of attraction forx30=x40=0, x20=x40=0, x20=x30=0, and x10=x40=0 corresponding to the asymmetric pair of chaotic (blue and green) and period-2 attractors (yellow and black) obtained for c=2.44. Red regions correspond to unbounded motion.
(a)(b)(c)(d)Table 4
Details of the coexistences observed in the novel hyperjerk system.
Fig.n0Type of coexistencesValues of control parameterInitial conditionx10,x20,x30,x40Figure16A symmetric pair of period-2 attractors, a symmetric pair of chaotic attractors, and a pair of fixed point.c=2.44(a)0,0,±44.4,0(b)0,0,±48,0(c)0,0,±26.4,0Figure17A symmetric pair of period-2 attractors, a symmetric pair of chaotic attractors, a symmetric pair of period-9 attractors, and a pair of fixed point.c=2.454(a)0,0,±25.2,0(b)0,0,±48,0(c)0,0,±10.8,0(d)0,0,±26.4,0Figure18A symmetric pair of period-2 attractors, a fixed point, a symmetric pair of period-12 attractors, 2 symmetric pairs of chaotic attractors, and a pair of fixed point.c=2.5442(A”A1,A2“)0,0,±45.6,0(b)0,0,±18,0(c)0,0,±49.2,0(d)0,0,±48,0(e)0,0,±26.4,0


### 3.3.2. Coexistence of Attractors with respect to Bifurcation Parametersm and b

In order to investigate the sensitivity of the new hyperjerk system in terms of the bifurcation parameterm, the other parameters are fixed as follows: a=1.8, b=3.8, c=2.442, d=1.35, and e=14.85. We found that the novel hyperjerk system can exhibit striking bifurcation sequences when varying the control parameter m in the range 2≤m≤4.18. With reference to Figure 11, the bifurcation diagram in black and the one in blue are obtained by increasing and decreasing the values of the parameter m, while the one in red is obtained by fixing the initial conditions at x10,x20,x30,x40=−4,0,0,0. A window of hysteretic dynamics can be identified in the range 2.6≤m≤2.9. The enlarged bifurcation diagram of Figure 12 clearly illustrates the domain of the coexistence of multiple attractors observed in the new hyperjerk system according to the bifurcation parameter m. Different methods used to plot these bifurcation diagrams are presented in Table 3. Up to six different periodic, chaotic, and point attractors can be obtained by only changing the initial conditions. For instance, sample phase portraits of the coexistence of six distinct attractors are presented in Figure 13.Figure 16
Coexistence of six different attractors (a pair of symmetric period-2 limit cycle, a pair of symmetric period-9 limit cycle, a pair of symmetric chaotic attractors, and a pair of point attractors). The values of others system parameters are fixed as follows:a=1.8, b=3.8, d=1.35, e=14.85, and m=2.2. Initial conditions are given in Table 4.
(a)(b)(c)Figure 17
Coexistence of eight different attractors (a pair of symmetric period-2 limit cycle, a pair of symmetric chaotic attractors, and a pair of point attractors). The parameters are the same as those in Figure16. Initial conditions are given in Table 4.
(a)(b)(c)(d)During the mathematical analyses, it has been shown that the Hopf bifurcation was depending on the control parameterm. The bifurcation diagram of Figure 13(a) clearly illustrates this phenomenon characterized by the stable state followed by the unstable state. Moreover, this control parameter also highlights the coexistence of multiple attractors exhibited by the new 4-D system. By considering the following sets of the parameters: a=1.8, b=3.5, c=2.442, d=1.35, e=15, and m=3, we discover that the new 4-D system displays the coexistence of four distinct chaotic and point attractors. The phase portraits of the Figure 14 and their corresponding cross section of the basin of attraction clearly show the coexistence phenomenon and also give the initial condition domain of each attractor. The green and black domains represent the initial conditions regions of the pair of chaotic attractors, while the yellow domain represents the initial condition regions of the pair of point attractors.Figure 18
Coexistence of ten different attractors (a pair of symmetric period-2 limit cycle, 2 pair of symmetric chaotic attractors, a pair of symmetric period-9 limit cycle, and a pair of point attractors). The parameters are the same as those in Figure16. Initial conditions are given in Table 4.
(a)(b)(c)(d)(e)


## 3.3.1. Coexistence of Attractors with respect to Bifurcation Parameterc

By changing the system parameters and consideringc as bifurcation parameter, we observe a completely different behavior. In addition, a very interesting phenomenon which is the coexistence of multiple attractors appears in the new 4-D hyperjerk. For this phenomenon to be illustrated, the bifurcation diagrams of Figure 8 are plotted using the following method:(i)
The blue diagram is obtained by simultaneously increasing the value of the control parameterc as well as the initial condition x0. At each iteration, we assign to x0 the new value of the control parameter c.(ii)
The red diagram is obtained by incrementingc from its minimum value 2.34 to its maximum value 2.985, with a carefully chosen step. Note that the solutions of the system at each iteration are considered as the initial condition of the next iteration.(iii)
The cyan diagram respects the previous procedure, with the initial condition−10.67,0,0,0, whereas the black diagram follows the same procedure as previously described but the only difference is the decrease of the control parameter c. The initial condition is 10,0,0,0.(iv)
The magenta diagram is obtained by increasing the control parameterc from 2.398 to 2.985, followed by decreasing the bifurcation parameter c from 2.398 to 2.34.Figure 8
symmetric coexisting bifurcation diagrams showing complex structure of the new hyperjerk fora=1.8, b=3.8, d=1.35, e=14.85, and m=2.2. Initial conditions y0=z0=w0=0 and x0 are fixed as follows for obtaining these diagrams: red diagrams for increasing c with x0=5.25, black diagram for decreasing c with x0=−10, cyan diagrams for increasing c with x0=−10.67, magenta diagram for increasing and decreasing, started at c=2.5442 with x0=−4, and blue diagram for starting c with same initial condition x0=c. For more information, see Table 3.
(a)(b)We can observe in Figure8 several windows of coexisting attractors. For more details about the methods used to plot the bifurcation diagrams of Figure 8, see Table 3. The enlarged bifurcation diagram of Figure 9 shows the hysteretic domain, plotted in the range 2.442≤c≤2.488, and the techniques used to plot the diagrams are also presented in Table 3. Figure 9 shows the coexistence of six and eight different limit cycles, chaotic and point attractors. Some sample phase portraits showing the coexistence of six and eight attractors are presented in Figures 10 and 11, respectively. Some basins of attractions showing the initial conditions domains of the coexisting attractors are presented in Figure 12. The coexistence of four attractors is clearly denoted (a pair of periodic attractors (black and yellow) and a pair of chaotic attractors (blue and green)). Note that there is a perfect symmetry between the different cross sections of the competing attractors.Table 3
Techniques used to obtain coexisting bifurcation diagrams and corresponding initial conditions.
Fig.n0Color graphParameter rangeSweeping directionInitial conditionx10,x20,x30,x40Figure8Blue2.34≤c≤2.985Upwardc,0,0,0Red2.34≤c≤2.985Upward5.25,0,0,0Cyan2.34≤c≤2.985Upward−10.67,0,0,0Black2.985≤c≤2.34Downward−10,0,0,0Magenta2.398≤c≤2.985Upward−4,0,0,02.398≤c≤2.34Downward−4,0,0,0Figure9Red2.442≤c≤2.488Upward5,0,0,0Black2.488≤c≤2.442Downward5,0,0,0Cyan2.442≤c≤2.488Upward−10.67,0,0,0Magenta2.442≤c≤2.488Upward12,0,0,0Figure13Blue2.802≤m≤2.865Upward−5.5,0,0,02.738≤m≤2.802DownwardRed2.802≤m≤2.865Upward12.8,0,0,02.738≤m≤2.802DownwardMagenta2.802≤m≤2.865Upward−10.1,0,0,02.738≤m≤2.802DownwardFigure 9
Enlargement of bifurcation diagram of Figure8 plotted in the range 2.442≤c≤2.488 showing the region in which the system exhibits multiple coexisting attractors. The system parameters are the same as the one of Figure 8.Figure 10
Bifurcation like sequence showing local maxima of the coordinatex3 versus the initial condition . plotted in the range 0≤x0≤30, with a=1.8, b=3.8, c=2.5442, d=1.35, and e=14.85. It can be observed the coexistence of ten periodic, chaotic, and point attractors.Figure 11
(a, b) Bifurcation diagrams showing local maxima of the coordinatex1 versus m and (c) the corresponding graph of largest Lyapunov exponent (λmax) plotted in the range 2≤m≤4.18, with a=1.8, b=3.8, c=2.442, and d=1.35.
(a)(b)(c)Figure 12
Enlargement of the bifurcation diagram of Figure11(a) showing the region in which the system exhibits multiple coexisting attractors. This region corresponds to values of m in the range: 2.735≤m≤2.865. Three sets of data are superimposed. For more information about the methods used to plot, see Table 3.The initial conditions of the coexistence of ten attractors exhibited by the proposed hyperjerk system are presented in Figure13. The bifurcation like sequence of Figure 13 shows the variation of x30 in terms of the control parameter c, and the other initial conditions are set to zero (x10=x20=x40=0). Note that chaotic attractors are characterized by an unlimited number of points, while periodic attractors are characterized by a finite number of points. In the same line, the basin of attraction shows the different domains of convergence of similar attractors. The phase portrait of Figure 14 illustrates the coexistence of ten attractors. The initials conditions of the coexisting attractors are given in Table 4.Figure 13
Coexistence of six different attractors (a pair of period-2 limit cycle, a pair of chaotic attractors, and a pair of fixed point attractor) for the following values of system parameters:a=1.8, b=3.8, c=2.442, d=1.35, e=14.85, and m=2.802. Initial conditions x10,x20,x30,x40 are, respectively, ±5.5,0,0,0, ±12.8,0,0,0, and ±10.1,0,0,0.
(a)(b)(c)(d)(e)Figure 14
(a, b) Coexistence of four distinct attractors (a pair of chaotic attractors represented inx1−x2 and x2−x3 plans), (c) a pair of fixed point attractor) and their corresponding basin of attraction. Initial conditions x10,x20,x30,x40 are, respectively,±7.15,0,±8.613,0 and ±12.8,0,0,0. The rest of the parameters area=1.8, b=3.5, c=2.442, d=1.35, e=15, and m=3.
(a)(b)(c)(d)Figure 15
Cross sections of basin of attraction forx30=x40=0, x20=x40=0, x20=x30=0, and x10=x40=0 corresponding to the asymmetric pair of chaotic (blue and green) and period-2 attractors (yellow and black) obtained for c=2.44. Red regions correspond to unbounded motion.
(a)(b)(c)(d)Table 4
Details of the coexistences observed in the novel hyperjerk system.
Fig.n0Type of coexistencesValues of control parameterInitial conditionx10,x20,x30,x40Figure16A symmetric pair of period-2 attractors, a symmetric pair of chaotic attractors, and a pair of fixed point.c=2.44(a)0,0,±44.4,0(b)0,0,±48,0(c)0,0,±26.4,0Figure17A symmetric pair of period-2 attractors, a symmetric pair of chaotic attractors, a symmetric pair of period-9 attractors, and a pair of fixed point.c=2.454(a)0,0,±25.2,0(b)0,0,±48,0(c)0,0,±10.8,0(d)0,0,±26.4,0Figure18A symmetric pair of period-2 attractors, a fixed point, a symmetric pair of period-12 attractors, 2 symmetric pairs of chaotic attractors, and a pair of fixed point.c=2.5442(A”A1,A2“)0,0,±45.6,0(b)0,0,±18,0(c)0,0,±49.2,0(d)0,0,±48,0(e)0,0,±26.4,0


## 3.3.2. Coexistence of Attractors with respect to Bifurcation Parametersm and b

In order to investigate the sensitivity of the new hyperjerk system in terms of the bifurcation parameterm, the other parameters are fixed as follows: a=1.8, b=3.8, c=2.442, d=1.35, and e=14.85. We found that the novel hyperjerk system can exhibit striking bifurcation sequences when varying the control parameter m in the range 2≤m≤4.18. With reference to Figure 11, the bifurcation diagram in black and the one in blue are obtained by increasing and decreasing the values of the parameter m, while the one in red is obtained by fixing the initial conditions at x10,x20,x30,x40=−4,0,0,0. A window of hysteretic dynamics can be identified in the range 2.6≤m≤2.9. The enlarged bifurcation diagram of Figure 12 clearly illustrates the domain of the coexistence of multiple attractors observed in the new hyperjerk system according to the bifurcation parameter m. Different methods used to plot these bifurcation diagrams are presented in Table 3. Up to six different periodic, chaotic, and point attractors can be obtained by only changing the initial conditions. For instance, sample phase portraits of the coexistence of six distinct attractors are presented in Figure 13.Figure 16
Coexistence of six different attractors (a pair of symmetric period-2 limit cycle, a pair of symmetric period-9 limit cycle, a pair of symmetric chaotic attractors, and a pair of point attractors). The values of others system parameters are fixed as follows:a=1.8, b=3.8, d=1.35, e=14.85, and m=2.2. Initial conditions are given in Table 4.
(a)(b)(c)Figure 17
Coexistence of eight different attractors (a pair of symmetric period-2 limit cycle, a pair of symmetric chaotic attractors, and a pair of point attractors). The parameters are the same as those in Figure16. Initial conditions are given in Table 4.
(a)(b)(c)(d)During the mathematical analyses, it has been shown that the Hopf bifurcation was depending on the control parameterm. The bifurcation diagram of Figure 13(a) clearly illustrates this phenomenon characterized by the stable state followed by the unstable state. Moreover, this control parameter also highlights the coexistence of multiple attractors exhibited by the new 4-D system. By considering the following sets of the parameters: a=1.8, b=3.5, c=2.442, d=1.35, e=15, and m=3, we discover that the new 4-D system displays the coexistence of four distinct chaotic and point attractors. The phase portraits of the Figure 14 and their corresponding cross section of the basin of attraction clearly show the coexistence phenomenon and also give the initial condition domain of each attractor. The green and black domains represent the initial conditions regions of the pair of chaotic attractors, while the yellow domain represents the initial condition regions of the pair of point attractors.Figure 18
Coexistence of ten different attractors (a pair of symmetric period-2 limit cycle, 2 pair of symmetric chaotic attractors, a pair of symmetric period-9 limit cycle, and a pair of point attractors). The parameters are the same as those in Figure16. Initial conditions are given in Table 4.
(a)(b)(c)(d)(e)


## 3.4. Offset Boosting Scenario

Another property of system (3) is the possibility to develop an offset boosting effect. In our model, x1 appears only in the fourth line of equation (3), and thus this variable is a bootable variable [47–52]. Assuming the transformation x1⟶x1+k where k is a constant, equation (5) can be rewritten accordingly as(14)x˙1=x2,x˙2=mx3,x˙3=dx4,x˙4=cx1+k−bx2−ex3−ax4−γsinhx1+k.Figure19 clearly presents offset boosting of the double-band chaotic attractor. The following values of parameter: k=0 (blue), k=10 (red) and k=−10 (green) are used to plot them in x1-x4 and x1-x3 planes.Figure 19
Offset boosting of the double-band chaotic attractor with the following values of parameter:k=0(blue), k=10(red), and k=−10(green). The other parameters are the same as those in Figure 4.
(a)(b)


## 3.5. Antimonotonicity

By decreasing the value of the control parametere15≤e≤20, we can observe the formation and destruction of periodic orbits via reverse period-doubling bifurcation sequences. This interesting phenomenon has been reported in the literature. It is reported in various nonlinear systems such as Duffing oscillator [12], Chua circuit [53], and second-order nonlinear nonautonomous circuit [54, 55]. This phenomenon was reported for the first time in the hyperjerk system by Leutcho et al. [21]. The creation of periodic seas in the parameter space is the necessary requirement for a nonlinear system to experience forward and reverse period-doubling cascade [21]. Sample illustrations are represented in Figure 20, where some bifurcation diagrams are shown. These diagrams are obtained for each discrete value of the control parameter c. In Figure 20, note that for c=2.5442, we have period-2 bubble and for a slight adjustment of the control parameter c, period-4 bubble is observed for c=2.8442, whereas for c=2.9, we have a period-8 bubble. In the same order, chaotic bubbles are formed for c=2.97,andc=2.99. The increase of control parameter c causes the creation of other bubbles, and it finally results in an infinite tree (like chaos).Figure 20
Bifurcation diagrams showing local maxima of the coordinatex1 of the attractor in Poincaré cross section in terms of the control parameter a (bubbling): (a) period-2 bubble for c=2.5442, (b) period-4 bubble for c=2.8442, (c) period-8 bubble for c=2.9, (d) single-band chaos bubble for c=2.93, (e) single-band chaos bubble for c=2.99, and (f) double-band chaos bubble for c=3.0.
(a)(b)(c)(d)(e)(f)


## 4. Experimental Study

The objective of this section is to confirm the above theoretical results by realizing a laboratory experimental study. For this purpose to be achieved, several approaches have been proposed in the literature to implement chaotic circuit (by using many types of off-the-shelf electronic components [56] or field-programmable gate array (FPGA) technology [57–60] or field-programmable-analog-array (FPAA) technology [61, 62] just to name a few). Only off-the-shelf electronic components (i.e., resistors, capacitors, pair of semiconductor diodes (D1=D2=1N4148), and TL084 operational amplifiers types with a power supply of ±15 VDC) are used to realize the schematic diagram of Figure 1. The following values of electronic circuit components are used during the experimental process: Rc=Rm=R=10kΩ, Rb=1.67kΩ, Re=5kΩ, and Rd=2.941kΩ (for the other parameters, see Table 1 case A). The complete sequence of phase portraits plotted in x1,x4 plan is obtained by adjusting the control resistor Ra in the range 333,33Ω≤Ra≤1.66kΩ. It can be seen in Figure 21 a good coherence between the numerical results (left side) and the experimental ones (right side). By changing the values of electronic components: Re=600Ω,Ra=2.116kΩ, Rb=4.21kΩ, Rm=784Ω, Rd=3.234kΩ, and Rc=5.16kΩ (for the other parameters, see Table 1 case B), the coexistence of attractors emerges. Figure 22 clearly illustrates the coexistence of fixed points, period-2 attractor, and chaotic attractor. Those attractors appear randomly by switching on and off the power supply. We can conclude that the mathematical model proposed in this work perfectly describes the real behavior of the novel hyperjerk circuit.Figure 21
Experimental phase portraits (right column) and corresponding numerical ones (left column) obtained by a direct integration of the system (1) confirming the scenario to chaos in the system for varying Ra (i.e., parameter a): (a) period-1 for Ra=1.428, (b) period-2 for Ra=807, (c) single-band chaos for Ra=684, (d) single–band chaos for Ra=500, and (e) double-band chaos for Ra=333,33. The scales are X=0.2V/div and Y=0.5V/div.
(a)(b)(c)(d)(e)Figure 22
Coexistence of multiples attractors forRc=5.16kΩ. Both periodic and chaotic attractors appear randomly in the experiment when switching on and off the power supply. The scales are X=5V/div and Y=2V/div for all pictures.


## 5. Conclusion

This work has proposed and investigated a new chaotic hyperjerk circuit with three equilibrium points having hyperbolic sine nonlinearity. The chaotic behavior observed in the system is due to the nonlinear component formed by two antiparallel diodes. Classical nonlinear analysis tools have been used to study the complete dynamics of the system. The bifurcation analysis of the new circuit shows that the chaotic double-band attractor arises from the period-doubling scenario followed by the symmetry recovering crisis event. In addition, some properties of the system such as antimonotonicity and offset boosting have been revealed. In particular, various regions in the parameter space in which the system develops the coexistence of up to ten disconnected attractors consisting of stable fixed points, limit cycles, and strange attractors have been reported. The coexistence of periodic, chaotic, and stable fixed points discovered in this work has not yet been reported in a hyperjerk system (at least as simple as the case discussed) and thus merits dissemination. To validate the theoretical study presented in this work, the new chaotic hyperjerk circuit has been realized and used for the investigations. Experimental results agree well with those obtained during the numerical experiment, thus confirming the feasibility of the proposed model. Owing to its extreme simplicity coupled with extremely rich dynamics, the new hyperjerk circuit introduced in this work has potential utility for information encryption as well as for other chaos-based applications [62].



---

*Source: 6182183-2020-07-29.xml*