Chaos, Solitons and Fractals
Chaos, Solitons and Fractals
Chaos, Solitons and Fractals
Frontiers
a r t i c l e i n f o a b s t r a c t
Article history: Dynamical systems with megastable properties are very rare in the literature. In this paper, we introduce
Received 4 November 2019 a new two-dimensional megastable dynamical system with a line of equilibria, having an infinite num-
Revised 8 February 2020
ber of stable states. By modifying this new system with temporally-periodic forcing term, a new two-
Accepted 17 February 2020
dimensional non-autonomous nonlinear oscillator capable to generate an infinite number of coexisting
limit cycle attractors, torus attractors and, strange attractors is constructed. The analog implementation
Keywords: of the new megastable oscillator is investigated to further support numerical analyses and henceforth
Forced oscillator validate the mathematical model.
Megastability
Self-excited attractors © 2020 Elsevier Ltd. All rights reserved.
Coexisting attractors
1. Introduction multistable when there are more than two coexisting attractors in
it [29,30]. Bistability represents the situation of two different coex-
Chaotic/hyperchaotic nonlinear dynamical systems have one of isting attractors. In this regard, many other works have been inves-
their main properties to be extreme sensitivity to noise or sensi- tigated in the literature on this interesting topic in the last decade
tive dependence to initial conditions [1–5]. Nonlinear dynamical [31–35]. Indeed, multistability is a very hot and ongoing research
systems have been recently classified into two major groups. That topic. In this scope, there are special cases of multistability which
is self-excited attractors systems and systems with hidden attrac- have attracted much attention in recent years. That is systems with
tors [6–11]. Self-excited attractors have a basin of attraction that extreme multistability (i.e., there is an infinite number of coexist-
is related to an unstable fixed point [3]. In contrast, when the ing attractors, and thus initial states cause bifurcations) [36].
delimitation regions don’t intersect with the small neighborhood Nonlinear oscillators with megastability are another type of ex-
of any equilibrium point [12], then that attractor is hidden. These citing multistable systems found recently in literature [37–41]. This
groups of attractors have been investigated in literature in different latter case has an infinite number of coexisting attractors, but there
systems including Hopfield neural networks [13,14], Chua’s circuit are no such bifurcations in them like systems with extreme multi-
[15], Sprott system [16], jerk [17,18], hyperjerk systems [19–23] and stability [37]. Indeed, the term megastable has been used for sys-
fractional-order chaotic systems [24–27]. In terms of their complex tems with countable infinite attractors [42–46], in contrast to ex-
and extremely rich dynamic, these chaotic systems can generate treme multistability which is related to non-countable infinite at-
different types of stable states for a specific parameter set when tractors [47]. From the point of view of applications, the coexis-
using random initial conditions. This latter characteristic of nonlin- tence of different stable states offers great flexibility in the sys-
ear dynamical systems is called multistability. Recall that multista- tem performance without major parameter changes; that can be
bility is very important in dynamical systems [1,28]. A system is exploited with the right control strategies to induce a definite
switching between different coexisting states [48–52].
Megastable property of the nonlinear dynamical system has
∗
Corresponding author. been investigated in some few cases very recently. In 2017 Sprott
E-mail address: [email protected] (G.D. Leutcho).
https://doi.org/10.1016/j.chaos.2020.109703
0960-0779/© 2020 Elsevier Ltd. All rights reserved.
2 G.D. Leutcho, S. Jafari and I.I. Hamarash et al. / Chaos, Solitons and Fractals 134 (2020) 109703
Fig. 4. A zoomed version of Fig. 3 showing the structure of strange attractor with
Fig. 1. Trajectories showing different coexisting limit cycles of system (2), for 75
omitted transition parts.
different initial conditions (located on a mesh x, y = −150 to x, y = +150 with steps
equal to 4). Each trajectory is for a duration of 20 0 0 sec.
Fig. 5. (a) Bifurcation diagrams of forced system with respect to forcing amplitude A and (b) corresponding maximum values of y versus forcing amplitude A. The initial
conditions are (1, 0) and (10, 0).
for this system in its steady point is an infinite number of stable states around x-axis, which could be
identified under other appropriate initial points.
0 −1
J=
k1 ω1 cos (ω1 x ) + k2 ω2 cos (ω2 x ) − y sin (x ) cos (x )
(3)
0 −1
(x, y ) = (0, 0) k ω + k ω 1
3. The forced chaotic oscillator
−−−−−−−−−−→ 1 1 2 2
Provide that the characteristic equation has coefficient with a With the forcing term introduced in system (5), this latter sys-
different sign [55–57], the equilibrium point O(0, 0) is unstable. For tem can display infinite coexisting chaotic trajectories for an ap-
our simulation, we have selected the following system parameters propriate choice of system parameters and initial conditions as
k1 = 1, k2 = 2, ω1 = √1 , and ω2 = 0.05. Some exciting trajectories well. In this regard, we have selected ω = 0.7and take some dis-
2
in this system are depicted in Fig. 1. This graph represents a phase tinct value of parameter A in the chaotic region in order to re-
portrait in System (2) for 75 distinct initial conditions located on veal these interesting features. In particular, For A = 0.55, system
the x-axis. We can remark that 30 stable states are identified (by (5) can effectively display an infinite number of coexisting chaotic
using those 75 initial conditions). That is, 3 packs, each with 10 and periodic stable states as presented in Fig. 2. These state space
limit cycles are uncovered. Note that, they are just examples in plots are obtained with the same initial conditions employed in
4 G.D. Leutcho, S. Jafari and I.I. Hamarash et al. / Chaos, Solitons and Fractals 134 (2020) 109703
Fig. 6. (a) Phase portrait showing the coexistence of eight stable states with their corresponding time series (b–c) using initial conditions (5, 0), (15, 0), (25, 0), (50, 0), (90,
0), (105, 0), (115, 0), (70, 0) for black, green, cyan, red, yellow, magenta and violet attractor. (d) Represents bifurcation like sequence showing local maxima of the coordinate
y versus x(0) for A = 0.3 and ω = 0.7. (For interpretation of references to color in this figure legend, the reader is referred to the web version of this article)
Fig. 1. One can remark in Fig. 2 that some of the limit cycles which torus attractors can be observed. the results of the investigation of
existed in Fig. 1, have lost their own delimitation regions. Indeed, the new megastable system are also provided in Fig. 7 through the
we have presented in Fig. 3 and Fig. 4, the main pack of coexist- two-parameter Lyapunov exponent diagrams in the same range of
ing solutions around the origin. That is a zoom version of trajecto- parameter space (A, k2 ). These diagrams are very important to bet-
ries shown in Fig. 2. The situation depicted in Fig. 4 is typically for ter localize the regions of multiple coexisting solutions in the sys-
chaotic trajectories. To further highlight these coexisting solutions tem. From the graphs in Fig. 7, it can be seen that for two different
for values of the control parameter A, we provide the bifurcation initial conditions, the coexisting solutions regions are materialized
of Fig. 5(a and b) obtained when monitoring this parameter (i.e., by the simultaneous appearance of different colors on figures in
forcing amplitude) for two distinct initial conditions in the range the same range of parameter space.
A ∈ [0, 0.6].
As shown in Fig. 5, the bifurcation points differ in some range
4. The circuit implementation
of the control parameter for these initial conditions. To analyze the
new megastable system, we also calculate the maximum Lyapunov
The analog implementation of the new megastable oscillator in-
exponents (MLEs) related to the bifurcation diagrams (see Fig. 5(c
vestigated in this work is addressed in this section. The introduced
and d)) using Wolf algorithm [58]. It is easy to observe that on the
two-dimensional system is constructed using electronic component
graph of the Lyapunov exponent diagram, the dynamical behavior
of the Pspice environment [59,60]. This electronic implementation
of the new megastable system alternates between torus, limit cy-
of the model is very singular because the model is build based
cles, and chaotic attractors when A is increased. One can remark
only on linear and trigonometric terms. The circuit of the new
that those regions are torus. They may look like chaos in the bi-
megastable oscillator is designed and shown in Fig. 8. The circuit of
furcation diagram, in the state-space plot, or even in time-series.
the new system is designed using two capacitors (C1 , C2 ), thirteen
However, the largest LE in them is zero (they have two zero LEs
resistors (R1 , ..., R13 ) including, six op-amps TL082CD, one multi-
and one negative LE). Fig. 6 shows some very interesting features
plier which can be implemented practically using AD633JN ver-
in this system. It is a plot of trajectories in the system (5), for
sions of the AD633 four-quadrant voltage multipliers chips used to
8 different initial conditions. The coexistence of limit cycles with
implement the nonlinear terms of our model. The signal (W) at the
G.D. Leutcho, S. Jafari and I.I. Hamarash et al. / Chaos, Solitons and Fractals 134 (2020) 109703 5
5. Conclusion
Fig. 8. The circuit schematic of the novel oscillator with megastable behavior.
6 G.D. Leutcho, S. Jafari and I.I. Hamarash et al. / Chaos, Solitons and Fractals 134 (2020) 109703
Fig. 9. Two dimensional projection of the attractor obtained from circuit realization of the new chaotic circuit and their corresponding time series for three different initial
conditions(5, 0)V, (10, 0)Vand (15, 0)V respectively.
Declaration of Competing Interest [2] Kengne J, Signing VF, Chedjou J, Leutcho G. Nonlinear behavior of a novel
chaotic jerk system: antimonotonicity, crises, and multiple coexisting attrac-
tors. Int J Dyn Control 2018;6:468–85.
The authors declare no conflict of interest. [3] Sprott JC. Elegant chaos: algebraically simple chaotic flows. World Scientific;
2010.
Acknowledgment [4] Wei Z, Moroz I, Sprott J, Akgul A, Zhang W. Hidden hyperchaos and elec-
tronic circuit application in a 5D self-exciting homopolar disc dynamo. Chaos
2017;27:033101.
The authors are grateful to the anonymous reviewers and Dr. [5] Wei Z, Moroz I, Sprott JC, Wang Z, Zhang W. Detecting hidden chaotic regions
Fozin Théophile (University of Buea) for their valuable comments, and complex dynamics in the self-exciting homopolar disc dynamo. Int J Bifurc
which helped improve the content of the present paper. Chaos 2017;27:1730 0 08.
[6] Akgul A, Pehlivan I. A new three-dimensional chaotic system without equilib-
rium points, its dynamical analyses and electronic circuit application. Tech Gaz
Authors contributions 2016;23:209–14.
[7] Bao B, Qian H, Xu Q, Chen M, Wang J, Yu Y. Coexisting behaviors of asymmetric
Authors have contributed equally for this manuscript attractors in hyperbolic-type memristor based Hopfield neural network. Front
Comput Neurosci 2017;11:81.
References [8] Jafari S, Pham V-T, Kapitaniak T. Multiscroll chaotic sea obtained from a simple
3D system without equilibrium. Int J Bifurc Chaos 2016;26:1650031.
[9] Jafari S, Sprott JC, Molaie M. A simple chaotic flow with a plane of equilibria.
[1] Kapitaniak T, Leonov GA. Multistability: uncovering hidden attractors.
Int J Bifurc Chaos 2016;26:1650098.
Springer; 2015.
G.D. Leutcho, S. Jafari and I.I. Hamarash et al. / Chaos, Solitons and Fractals 134 (2020) 109703 7
[10] Pham VT, Volos C, Kapitaniak T, Jafari S, Wang X. Dynamics and circuit [36] Chen M, Sun M, Bao B, Wu H, Xu Q, Wang J. Controlling extreme multista-
of a chaotic system with a curve of equilibrium points. Int J Electron bility of memristor emulator-based dynamical circuit in flux–charge domain.
2018;105:385–97. Nonlinear Dyn 2018;91:1395–412.
[11] Wei Z, Pham V-T, Kapitaniak T, Wang Z. Bifurcation analysis and circuit realiza- [37] He S, Li C, Sun K, Jafari S. Multivariate multiscale complexity analysis of self-
tion for multiple-delayed Wang–Chen system with hidden chaotic attractors. -reproducing chaotic systems. Entropy 2018;20:556.
Nonlinear Dyn 2016;85:1635–50. [38] Li C, Sprott JC. An infinite 3-D quasiperiodic lattice of chaotic attractors. Phys
[12] Bayani A, Rajagopal K, Khalaf AJM, Jafari S, Leutcho G, Kengne J. Dynami- Lett A 2018;382:581–7.
cal analysis of a new multistable chaotic system with hidden attractor: An- [39] Li C, Sprott JC, Hu W, Xu Y. Infinite multistability in a self-reproducing chaotic
timonotonicity, coexisting multiple attractors, and offset boosting. Phys Lett A system. Int J Bifurc Chaos 2017;27:1750160.
2019;383:1450–6. [40] Li C, Sprott JC, Mei Y. An infinite 2-D lattice of strange attractors. Nonlinear
[13] Njitacke Z, Kengne J. Complex dynamics of a 4D Hopfield neural networks Dyn 2017;89:2629–39.
(HNNs) with a nonlinear synaptic weight: Coexistence of multiple attractors [41] Li C, Thio WJ-C, Sprott JC, Iu HH-C, Xu Y. Constructing infinitely many attrac-
and remerging Feigenbaum trees. AEU-Int J Electron Commun 2018;93:242–52. tors in a programmable chaotic circuit. IEEE Access 2018;6:29003–12.
[14] Njitacke Z, Kengne J. Nonlinear dynamics of three-neurons-based Hopfield [42] Sprott JC, Jafari S, Khalaf AJM, Kapitaniak T. Megastability: coexistence of a
neural networks (HNNs): Remerging Feigenbaum trees, coexisting bifurcations countable infinity of nested attractors in a periodically-forced oscillator with
and multiple attractors. J Circuit, Syst Comput 2019;28:1950121. spatially-periodic damping. Eur Phys J Spec Top 2017;226:1979–85.
[15] Bao B, Xu L, Wu Z, Chen M, Wu H. Coexistence of multiple bifurcation [43] Tang Y, Abdolmohammadi HR, Khalaf AJM, Tian Y, Kapitaniak T. Carpet oscilla-
modes in memristive diode-bridge-based canonical Chua’s circuit. Int J Elec- tor: a new megastable nonlinear oscillator with infinite islands of self-excited
tron 2018;105:1159–69. and hidden attractors. Pramana 2018;91:11.
[16] Lai Q, Xu G, Pei H. Analysis and control of multiple attractors in Sprott B sys- [44] Tang Y-X, Khalaf AJM, Rajagopal K, Pham V-T, Jafari S, Tian Y. A new nonlinear
tem. Chaos Solitons Fractals 2019;123:192–200. oscillator with infinite number of coexisting hidden and self-excited attractors.
[17] Kengne J, Negou AN, Tchiotsop D. Antimonotonicity, chaos and multiple at- Chin Phys B 2018;27:040502.
tractors in a novel autonomous memristor-based jerk circuit. Nonlinear Dyn [45] Wang Z, Abdolmohammadi HR, Alsaadi FE, Hayat T, Pham V-T. A new os-
2017;88:2589–608. cillator with infinite coexisting asymmetric attractors. Chaos Solitons Fractals
[18] Kengne J, Njitacke Z, Fotsin H. Dynamical analysis of a simple autonomous jerk 2018;110:252–8.
system with multiple attractors. Nonlinear Dyn 2016;83:751–65. [46] Wei Z, Pham V-T, Khalaf AJM, Kengne J, Jafari S. A modified multistable chaotic
[19] Kengne J, Leutcho GD, Telem ANK. Reversals of period doubling, co- oscillator. Int J Bifurc Chaos 2018;28:1850085.
existing multiple attractors, and offset boosting in a novel memristive [47] Njitacke Z, Kengne J, Tapche RW, Pelap F. Uncertain destination dynam-
diode bridge-based hyperjerk circuit. Analog Integr Circuits Signal Process ics of a novel memristive 4D autonomous system. Chaos Solitons Fractals
2019;101:379–99. 2018;107:177–85.
[20] Leutcho G, Kengne J, Kengne LK. Dynamical analysis of a novel autonomous [48] Canavier CC, Baxter DA, Clark JW, Byrne JH. Control of multistability in ring
4-D hyperjerk circuit with hyperbolic sine nonlinearity: chaos, antimono- circuits of oscillators. Biol Cybern 1999;80:87–102.
tonicity and a plethora of coexisting attractors. Chaos Solitons Fractals [49] Fonzin Fozin T, Kengne R, Kengne J, Srinivasan K, Souffo Tagueu M, Pelap F.
2018;107:67–87. Control of multistability in a self-excited memristive hyperchaotic oscillator.
[21] Leutcho GD, Kengne J. A unique chaotic snap system with a smoothly ad- Int J Bifurc Chaos 2019;29:1950119.
justable symmetry and nonlinearity: chaos, offset-boosting, antimonotonicity, [50] Fonzin Fozin T, Megavarna Ezhilarasu P, Njitacke Tabekoueng Z, Leutcho G,
and coexisting multiple attractors. Chaos Solitons Fractals 2018;113:275–93. Kengne J, Thamilmaran K, et al. On the dynamics of a simplified canonical
[22] Leutcho GD, Kengne J, Kengne R. Remerging Feigenbaum trees, and multiple Chua’s oscillator with smooth hyperbolic sine nonlinearity: hyperchaos, multi-
coexisting bifurcations in a novel hybrid diode-based hyperjerk circuit with stability and multistability control. Chaos 2019;29:113105.
offset boosting. Int J Dyn Control 2019;7:61–82. [51] Fozin TF, Leutcho G, Kouanou AT, Tanekou G, Kengne R, Kengne J, et al. Multi-
[23] Ren S, Panahi S, Rajagopal K, Akgul A, Pham V-T, Jafari S. A new chaotic stability control of hysteresis and parallel bifurcation branches through a linear
flow with hidden attractor: the first hyperjerk system with no equilibrium. augmentation scheme. Zeitschrift für Naturforschung A 2019;75:11–21.
Zeitschrift für Naturforschung A 2018;73:239–49. [52] Pisarchik AN, Kuntsevich BF. Control of multistability in a directly modulated
[24] He S, Sun K, Peng Y. Detecting chaos in fractional-order nonlinear systems us- diode laser. IEEE J Quantum Electron 2002;38:1594–8.
ing the smaller alignment index. Phys Lett A 2019;383:2267–71. [53] Rajagopal K, Singh JP, Roy BK, Karthikeyan A. Dissipative and conservative
[25] Peng D, Sun K, He S, Zhang L, Alamodi AO. Numerical analysis of a simplest chaotic nature of a new quasi-periodically forced oscillator with megastabil-
fractional-order hyperchaotic system. Theor Appl Mech Lett 2019;9:220–8. ity. Chin J Phys 2019;58:263–72.
[26] Peng D, Sun KH, Alamodi AO. Dynamics analysis of fractional-order perma- [54] Kahn PB, Zarmi Y. Nonlinear dynamics: exploration through normal forms.
nent magnet synchronous motor and its DSP implementation. Int J Mod Phys Courier Corporation; 2014.
B 2019;33:1950031. [55] Dolvis L, Vaidyanathan S, Jacques K, Sambas A, Mamat M. A New 4-D Hyper-
[27] Zhang L-M, Sun K-H, Liu W-H, He S-B. A novel color image encryption scheme chaotic System with Four-Scroll Hidden Attractor, Its Properties and Bifurca-
using fractional-order hyperchaotic system and DNA sequence operations. Chin tion Analysis. IOP conference series: materials science and engineering. IOP
Phys B 2017;26:100504. Publishing; 2019.
[28] Pisarchik AN, Feudel U. Control of multistability. Phys Rep 2014;540:167–218. [56] Wei Z, Li Y, Sang B, Liu Y, Zhang W. Complex dynamical behaviors in a 3D
[29] Bao H, Wang N, Bao B, Chen M, Jin P, Wang G. Initial condition-dependent dy- simple chaotic flow with 3D stable or 3D unstable manifolds of a single equi-
namics and transient period in memristor-based hypogenetic jerk system with librium. Int J Bifurc Chaos 2019;29:1950095.
four line equilibria. Commun Nonlinear Sci Numer Simul 2018;57:264–75. [57] Wei Z, Yu P, Zhang W, Yao M. Study of hidden attractors, multiple limit cycles
[30] Lai Q, Chen S. Research on a new 3D autonomous chaotic system with coex- from Hopf bifurcation and boundedness of motion in the generalized hyper-
isting attractors. Optik 2016;127:30 0 0–4. chaotic Rabinovich system. Nonlinear Dyn 2015;82:131–41.
[31] Bao B, Jiang T, Xu Q, Chen M, Wu H, Hu Y. Coexisting infinitely many at- [58] Wolf A, Swift JB, Swinney HL, Vastano JA. Determining Lyapunov exponents
tractors in active band-pass filter-based memristive circuit. Nonlinear Dyn from a time series. Physica D 1985;16:285–317.
2016;86:1711–23. [59] Vaidyanathan S, Dolvis LG, Jacques K, Lien C-H, Sambas A. A new five-dimen-
[32] Fonzin TF, Srinivasan K, Kengne J, Pelap F. Coexisting bifurcations in a mem- sional four-wing hyperchaotic system with hidden attractor, its electronic cir-
ristive hyperchaotic oscillator. AEU-Int J Electron Commun 2018;90:110–22. cuit realisation and synchronisation via integral sliding mode control. Int J
[33] Li C, Sprott JC. Coexisting hidden attractors in a 4-D simplified Lorenz system. Model Ident Control 2019;32:30–45.
Int J Bifurc Chaos 2014;24:1450034. [60] Vaidyanathan S, Pehlivan I, Dolvis LG, Jacques K, Alcin M, Tuna M, et al. A
[34] Negou AN, Kengne J. Dynamic analysis of a unique jerk system with a novel ANN-based four-dimensional two-disk hyperchaotic dynamical system,
smoothly adjustable symmetry and nonlinearity: reversals of period dou- bifurcation analysis, circuit realisation and FPGA-based TRNG implementation.
bling, offset boosting and coexisting bifurcations. AEU-Int J Electron Commun Int J Comput Appl Technol 2020;62:20–35.
2018;90:1–19.
[35] Pham V-T, Akgul A, Volos C, Jafari S, Kapitaniak T. Dynamics and circuit real-
ization of a no-equilibrium chaotic system with a boostable variable. AEU-Int
J Electron Commun 2017;78:134–40.