Printer Friendly

A Novel Autonomous 5-D Hyperjerk RC Circuit with Hyperbolic Sine Function.

1. Introduction

The study of three-dimensional dynamical systems seems to be mature [1-15]. The interest is now focused on high dimensional systems [16-21]. The reason is that they have been discovered to model natural phenomena more explicitly than the three-dimensional dynamical systems. Recently, hyperjerk systems have received great attention. In 2005, Konstantinos and Sprott proposed a class of chaotic/hyp erchaotic hyperjerk systems [19]. They investigated the dynamical properties of the proposed systems focusing on potential chaotic and hyperchaotic dynamics. They also claimed to have proposed surprisingly simple and most elegant functional forms of hyperjerk systems. A year later, Stefan showed that the concept of hyperjerk systems might serve as an appropriate starting point to study the dynamics of driven or coupled oscillators in a unified way [22]. The five-dimensional (5-D) hyperchaotic Sprott B system has ten terms including two quadratic nonlinear terms and only one control parameter. Hong Mey and collaborators have recently proposed a novel cryptosystem based on 5-D hyperchaotic system [23]. The system is combined with the logistic map for the generation of pseudo random sequences of better properties. In order to analyze the behavior of chaotic systems, some researchers focused on the synchronization and control of the 5-D systems [23-26]. These systems are algebraically simpler than Lorenz and Lorenz-like 5-D systems with twelve terms, three quadratic nonlinear terms including five or six control parameters. Let us note that the practical realization of quadratic nonlinearity found in these works is expensive due to the presence of analogue multipliers. Recently multistability has attracted tremendous research efforts [1, 5,12]. However to the best of the author's knowledge, multistability in 5-D hyperjerk systems is very little documented. Also it is very interesting to question whether there exists a simple chaotic 5-D hyperjerk system capable of multiple coexisting attractors.

This paper investigates the dynamics of a novel 5-D hyperjerk circuit with a very simple nonlinear part (a pair of semiconductor diodes). The new circuit can be regarded as a 5-D version of the jerk circuit previously reported by Kengne and collaborators [27]. It is important to stress that the simplicity of the 5-D autonomous circuit relies on the practical standpoint where the nonlinear part is constructed using semiconductor diodes. Despite the relative simplicity of their electronic circuits, the proposed circuit is characterized by rich and striking nonlinear behaviors such as chaos, anti-monotonicity, multiple coexisting attractors.

The paper is organized as follows: Section 2 deals with the modeling process. The electronic structure of the oscillator is described and a suitable mathematical model is derived to investigate the dynamics of the system. Section 3 is concerned with the numerical analysis. Various bifurcation diagrams combined with their corresponding graphs of the largest Lyapunov exponent are plotted to reveal different transitions to chaos. The occurrence of multiple attractors is also discussed using bifurcation diagrams as arguments. In addition the bubbling phenomenon is presented. Pspice circuit simulations results are carried out in Section 4. Some concluding remarks are presented in Section 5.

2. Presentation of the Novel 5-D Circuit

The electronic circuit of the oscillator under investigation is depicted in Figure 1. It consists of five successive active integrators in a feedback loop. Additionally, another feedback loop involving two of the integrators and a pair of diodes (e.g., type 1N4148) connected in antiparallel is used to synthetize the nonlinear part. The pair of semiconductor diodes is the only element responsible for the complex dynamics of the circuit. For instance, up to seven disconnected attractors have been found depending solely on the system's initial state. Upon applying Kirchhoff's electrical circuit laws and the Shockley diode equation [26] to the circuit of Figure 1, the following set of five coupled first-order differential equations can be derived:

[C.sub.1] [dv.sub.1]/dt = [v.sub.2]/R

[C.sub.2] [dv.sub.2]/dt = [v.sub.3]/R

[C.sub.3] [dv.sub.3]/dt = [v.sub.4]/R

[C.sub.4] [dv.sub.4]/dt = [v.sub.5]/R

[C.sub.5] [dv.sub.4]/dt = [v.sub.5]/[R.sub.[a.sub.0]] - [v.sub.3]/[R.sub.[a.sub.1]] - [v.sub.2]/[R.sub.[a.sub.2]] - [v.sub.1]/[R.sub.aa] - 2[I.sub.s] sinh ([v.sub.4]/[eta][V.sub.T]) (1)

To ease the numerical analysis of the circuit, the following change of variable and parameters is considered:

[v.sub.i] = [x.sub.i][V.sub.ref] (i = 1 ... 5);

t = RC[tau];

[V.sub.ref] = [eta][V.sub.T];

b = R/[R.sub.b];

[a.sub.j] = R/[R.sub.aj] (j = 0 ... 3);

[gamma] = 2R[I.sub.s]/[V.sub.ref] (2)

Thus the dimensionless circuit equations are derived as follows:

[[??].sub.1] = [x.sub.2]

[[??].sub.2] = [x.sub.3]

[[??].sub.3] = [x.sub.4]

[[??].sub.4] = [x.sub.5]

[[??].sub.5] = -[a.sub.0][x.sub.5] - [a.sub.1][x.sub.3] - [a.sub.2][x.sub.2] - [a.sub.3][x.sub.1] - [gamma] sinh ([x.sub.4]) (3)

The simplicity of the new 5-D autonomous system (3) studied in this paper mainly relies on its nonlinear part which is synthetized using only two semiconductor diodes. With reference to some chaotic systems in the literature [16 -21], cubic or quadratic nonlinear product terms are mostly used. Accordingly, the practical construction is more complicated due to the use of multipliers. In addition, the maximum Lyapunov exponent of the studied 5-D chaotic system is bigger than that of some dynamical systems. Hence, it has more complex dynamic behaviors. Table 2 is provided for more illustration.

System (3) is dissipative since its divergence ([nabla]V = -[a.sub.0]) is negative. Consequently, all system orbits will be confined to a specific bounded subset of zero volume in state space and the asymptotic dynamics converges to an attractor. This is true for this particular case. However, the system is dissipative in the sense of Levinson if there exists nonlocal Lyapunov function which defines the absorbing set [23]. In addition, System (3) is invariant under the transformation: ([x.sub.1], [x.sub.2], [x.sub.3], [x.sub.4], [x.sub.5]) [??] (-[x.sub.1], -[x.sub.2], -[x.sub.3], -[x.sub.4], -[x.sub.5]). Consequently, if ([x.sub.1], [x.sub.2], [x.sub.3], [x.sub.4], [x.sub.5]) is a solution of system (3) for a given set of parameters, then (-[x.sub.1], -[x.sub.2], -[x.sub.3], -[x.sub.4], -[x.sub.5]) is also a solution for the same parameter set. This property is responsible for the variety of coexisting attractors in the system. Finally, it is found that the novel oscillator under investigations has only one equilibrium point that is the origin [E.sub.0](0, 0,0, 0).

3. System's Dynamics

3.1. Computational Method. To produce phase portraits, bifurcation diagrams, and Lyapunov spectrum, the system was solved using the classical fourth-order Runge-Kutta algorithm with the time step always [DELTA]t [less than or equal to] 0.005. Two different methods are used to plot the various bifurcation diagrams in order to highlight the sensitivity of the whole system to the tiny changes of its parameters. The first method is by forward and backward continuations, i.e., using a properly continuation package for continuous-time systems [28], which is a standard tool for numerical bifurcation analysis. The second method is by following a specific attractor to obtain its behavior at long term. The later method is used to justify the abundant coexisting attractors in the system (up to seven coexisting attractors).

3.2. Route to Chaos and Antimonotonicity. To reveal the type of transition leading to chaos, a single control parameter (b) was considered to vary in the range 2.4 < [a.sub.2] < 3 while the rest of system's parameters are fixed as follows: [a.sub.0] = 1.75; [a.sub.3] = 1; [a.sub.1] = b = 3; [gamma] = 0.0011. The bifurcation diagram of Figure 2 is obtained by forward continuation. From this diagram it is obvious that the system experiences the classical reverse period doubling bifurcation. It moves from period-1 limit cycle to double band chaotic attractors going through single band chaotic attractors and tiny windows of periodic attractors. Figure 3 (left) illustrates this bifurcation sequence with the numerical phase portraits.

As the system experiences the classical period doubling route to chaos, it is obvious that antimonotonicity can be observed. Represented on Figure 4, this behavior has been searched by varying parameter b in the range 16 to 24 for some discrete values of parameter [a.sub.3]. For instance, if [a.sub.3] = 3 period-2 bubble is observed and each branch develops two branches leading to a stable period-4 bubble for [a.sub.3] = 4. Similarly a period-8 bubble is obtained for [a.sub.3] = 4.2. As the discrete parameter b increases, more bubbles are created until an infinitely tree (chaos) finally occurs when [a.sub.3] = 4.5.

3.3. Multiple Coexisting Attractors. Multistability has been previously revealed in many dynamical systems [29-34], including jerk systems with self-excited attractors [30, 33,34] and hidden attractors [3, 6, 35]. But the models presented up to date achieved at most six coexisting attractors in the phase space. This behavior has not yet been discovered in any 5-D hyperjerk systems. The system under study can challenge this limit by displaying up to seven coexisting attractors depending solely on the system's initial state.

With reference to the bifurcation diagram of Figure 5 the classical forward or backward continuation of parameter b is obtained with the following initial conditions [x.sub.1](0) = 3; [x.sub.2](0) = [x.sub.3](0) = [x.sub.4](0) = [x.sub.5](0) = 0 and [x.sub.1](0) = 1; [x.sub.2](0) = [x.sub.3](0) = [x.sub.4](0) = [x.sub.5](0) = 0. A hysteretic window can be identified leading to the coexistence of four disconnected chaotic and periodic attractors (Figure 6) in the phase space for b = 29.57, [a.sub.0] = [a.sub.2] =[a.sub.3] = 5, [a.sub.1] = 7, [gamma] = 5.4433 x [10.sup.-4].

To observe more than four different attractors in the system, the second method described above (see Section 3.2) is used to plot the bifurcation diagrams of Figures 7 and 12 by varying parameter [a.sub.2]. If we fix systems parameters as [a.sub.2] = 2.458; [a.sub.0] = 1.75; [a.sub.1] = b = 3; [a.sub.3] = 1; [gamma] = 0.0054, the system experiences four different periodic and chaotic attractors (see Figure 8) where the cross section of the basins of attraction of the attractors is also presented. This dynamics is justified using the bifurcation diagrams of Figure 7.

If systems parameters are [a.sub.2] = 3; [a.sub.0] = 1.75; [a.sub.1] = b = 3; [a.sub.3] = 1; [gamma] = 0.0054, the system displays five different period-1 limit cycles (see Figure 9). This dynamics is justified using the bifurcation diagrams of Figure 7.

For [a.sub.2] = 2.71; [a.sub.0] = 1.75; [a.sub.1] = b = 3;[a.sub.3] = 1; [gamma] = 0.0054, the system experiences five different periodic and chaotic attractors (see Figure 10). This dynamics is justified using the bifurcation diagrams of Figure 7.

For [a.sub.2] = 2.8; [a.sub.0] = 1.75; [a.sub.1] = b = 3; [a.sub.3] = 1; [gamma] = 0.0054, the system displays six different periodic and chaotic attractors (see Figure 11). This dynamics is justified using the bifurcation diagrams of Figure 7. The corresponding attraction basin is shown in Figure 7. Let us note that an attractor is hidden if their basin of attraction does not intersect with small neighborhoods of the equilibrium points of system. Given that the single equilibrium point of system (3) E0(0,0,0, 0) intersects with the basin of attraction of the coexisting attractor Figure 7, we can conclude that the said attractor is self-excited (instead of hidden attractors). Some literature can provide systems with hidden attractors [3, 6, 35].

More interestingly, another window of multiple coexisting attractors can be revealed when [a.sub.0] = 1.75; [a.sub.1] = b = 3; [a.sub.3] = 1; [gamma] = 0.0109 and the control parameter [a.sub.2] is varied in a tiny range. The bifurcation diagrams of Figure 12 revealed a window in which up to seven disconnected attractors (see Figure 13) coexist in the phase space depending solely on the system's initial state. Table 2 provides the initial conditions for these coexisting attractors.

This work represents an enriching contribution to the understanding of the nonlinear dynamics of this type of oscillators [36]. However, this striking phenomenon of disconnected coexisting attractors is also reported in other nonlinear dynamic systems such as lazer system [37], chemical reaction [38], and the radio physical system [35]. A special case where infinitely many attractors coexist, also referred to as extreme multistability, is discussed in [26, 39]. The multiplicity of attractors represents an additional type of randomness [40] that is exploited in real applications such as chaos based secret communication, image encryption, and random signal generation as well. However, this type of behavior is not desirable in general, thus the need for control. Detailed analysis on this line is out of the scope of this paper. Also, we suggest the excellent work on control of multistability by [41] to interested readers.

4. Pspice Circuit Simulations

Our motivation in this section is to verify the theoretical/numerical results obtained previously by performing some Pspice based simulations of the circuit. Furthermore, it is important to evaluate the effects of simplifying assumptions (e.g., ideal diode model and ideal op. amplifiers) considered during the mathematical modeling process, on the behavior of a hardware prototype of the 5-D hyperjerk circuit in Pspice. Briefly recall that an interesting aspect of using Pspice is the possibility of setting initial capacitors' voltages and analyzing the corresponding influence on the dynamics of the complete electronic circuit. Thus, the presence of multiple coexisting attractors can be tracked in a straightforward manner.

First to report the reverse period doubling routes to chaos observed during the numerical analysis, the circuit of Figure 1 has been simulated with the following electronic components: R = 10k[ohm], C = 10[eta]F, [R.sub.a0] = 5.71k[ohm], [R.sub.a1] = [R.sub.b] = -3.33k[ohm], [R.sub.a3] = 10k[ohm]. By varying [R.sub.a2], the complete routes to chaos are obtained and depicted in Figure 3-right. For [R.sub.a2] = -15.5k[ohm] a period-1 limit cycle is obtained, for [R.sub.a2] = 16.2k[ohm] a period-2 limit cycle is obtained, and chaotic attractors are obtained for [R.sub.a2] = 17k[ohm] and [R.sub.a2] = 18k[ohm].

Secondly, coexistence of multiple attractors can also be confirmed by Pspice based simulations with the following electronic circuit components: R = 50k[ohm], C = 2[eta]F, [Ra.sub.0] = -28.57k[ohm]; [R.sub.a1] = 16.66k[ohm]; [R.sub.a3] = 50k[ohm]; [R.sub.b] = 16.66k[ohm]. If [R.sub.a2] is fixed to [R.sub.a2] = 17.91k[ohm], four disconnected chaotic and periodic attractors (see Figure 14) coexist depending solely on the system's initial states as indicated in Table 2. For [R.sub.a2] = 15.6k[ohm], five different period-1 limit cycles can be observed (see Figure 15). The corresponding initial states are indicated in Table 2. If [R.sub.a2] is fixed to [R.sub.a2] = 16.5k[ohm], six disconnected chaotic and periodic attractors (see Figure 16) coexist depending solely on the system's initial states as indicated in Table 2. A very good similarity between numerical phase portraits and Pspice simulation results can be observed. However, slight discrepancies that may be attributed to the simplifying assumptions adopted during the modeling process can be noted between the bifurcations points in Pspice compared to the results from the theoretical analysis. It is important to stress that while Pspice software is based on actual circuit components, it still suffers from the discretization and its use can lead to wrong conclusions especially for hidden attractors (the same is true for MATLAB) [42].

5. Conclusion

A novel 5-D hyperjerk circuit with a very simple nonlinear part has been introduced in this work. The circuit is obtained by introducing additional feedback loops in the realization circuit of a jerk system previously reported by J. Kengne and collaborators. The modification yields the simplest 5-D hyperjerk system reported up to date. More interestingly, for some given sets of parameters, the system experiences a plethora of multiple coexisting attractors. For instance, up to seven disconnected attractors coexist in the system depending solely on the initial conditions. To the best of author's knowledge, such dynamics has not yet been reported in any hyperjerk system and thus deserves dissemination. Pspice based simulations were carried out to support the theoretical analysis. A detailed exploration of the parameter space (both experimentally and numerically) in view of revealing hyperchaotic behavior and hidden attractors in system (3) deserves further studies.

https://doi.org/10.1155/2018/1260325

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

[1] D. Dudkowski, S. Jafari, T. Kapitaniak, N. V. Kuznetsov, G. A. Leonov, and A. Prasad, "Hidden attractors in dynamical systems," Physics Reports, vol. 637, pp. 1-50, 2016.

[2] K. Srinivasan, V. K. Chandrasekar, A. Venkatesan, and I. Raja Mohamed, "Duffing-van der Pol oscillator type dynamics in Murali-Lakshmanan-Chua (MLC) circuit," Chaos, Solitons & Fractals, vol. 82, pp. 60-71, 2016.

[3] M. Chen, J. Yu, and B. Bao, "Hidden dynamics and multi-stability in an improved third-order Chua's circuit," The Journal of Engineering, vol. 2015, no. 10, pp. 322-324, 2015.

[4] C. Li, X. Wang, and G. Chen, "Diagnosing multistability by offset boosting," Nonlinear Dynamics, vol. 90, no. 2, pp. 1335-1341, 2017.

[5] S. T. Kingni, S. Jafari, H. Simo, and P. Woafo, "Three-dimensional chaotic autonomous system with only one stable equilibrium: analysis, circuit design, parameter estimation, control, synchronization and its fractional-order form," The European Physical Journal Plus, vol. 129, no. 5, pp. 1-16, 2014.

[6] B.-C. Bao, P. Jiang, Q. Xu, and M. Chen, "Hidden attractors in a practical Chua's circuit based on a modified Chua's diode," IEEE Electronics Letters, vol. 52, no. 1, pp. 23-25,2016.

[7] K. Murali and M. Lakshmanan, "Bifurcation and chaos of the sinusoidally driven Chua's circuit," International Journal of Bifurcation and Chaos, vol. 01, no. 02, pp. 369-384,1991.

[8] X. Wang, A. Akgul, S. Kacar, and V.-T. Pham, "Multimedia Security Application of a Ten-Term Chaotic System without Equilibrium," Complexity, vol. 2017, Article ID 8412093,10 pages, 2017.

[9] B. Bao, Q. Li, N. Wang, and Q. Xu, "Multistability in Chuas circuit with two stable node-foci," Chaos, vol. 26, Article ID 043111, 2016.

[10] Abir Lassoued and Olfa Boubaker, "Dynamic Analysis and Circuit Design of a Novel Hyperchaotic System with Fractional-Order Terms," Complexity, vol. 2017, Article ID 3273408, 10 pages, 2017.

[11] B. Bao, H. Wu, L. Xu, M. Chen, and W. Hu, "Coexistence of Multiple Attractors in an Active Diode Pair Based Chua's Circuit," International Journal of Bifurcation and Chaos, vol. 28, no. 02, Article ID 1850019,13 pages, 2018.

[12] M. Chen, Q. Xu, Y. Lin, and B. Bao, "Multistabiiity induced by two symmetric stable node-foci in modified canonical Chua's circuit," Nonlinear Dynamics, vol. 87, no. 2, pp. 789-802, 2017.

[13] A. Lassoued and O. Boubaker, "On new chaotic and hyperchaotic systems: a literature survey," Lithuanian Association of Nonlinear Analysts. Nonlinear Analysis: Modelling and Control, vol. 21, no. 6, pp. 770-789, 2016.

[14] Y. Zhang and J. Xu, "Classification and Computation of Non-resonant Double Hopf Bifurcations and Solutions in Delayed van der Pol-Duffing System," International Journal of Nonlinear Sciences and Numerical Simulation, vol. 6, no. 1, pp. 67-74,2005.

[15] T. Fozin Fonzin, K. Srinivasan, J. Kengne, and F. Pelap, "Coexisting bifurcations in a memristive hyperchaotic oscillator," AEU-International Journal of Electronics and Communications, vol. 90, pp. 110-122, 2018.

[16] S. Vaidyanathan, C. Volos, and V.-T. Pham, "Hyperchaos, adaptive control and synchronization of a novel 5-D hyperchaotic system with three positive Lyapunov exponents and its SPICE implementation," Archives of Control Sciences, vol. 24, no. 4, pp. 409-446, 2014.

[17] J. Kengne, N. Tsafack, and L. K. Kengne, "Dynamical analysis of a novel single Opamp-based autonomous LC oscillator: anti-monotonicity, chaos, and multiple attractors," International Journal of Dynamics and Control, pp. 1-15, 2018.

[18] C. Li and J. C. Sprott, "Coexisting hidden attractors in a 4-D simplified lorenz system," International Journal of Bifurcation and Chaos, vol. 24, no. 3, Article ID 1450034, pp. 145-161, 2014.

[19] E. Konstantinos and J. C. Sprott, "Chaotic hyperjerk systems," Chaos, Solitons & Fractals, vol. 28, no. 3, pp. 739-746, 2006.

[20] M. Buncha and S. Banlue, "Elementary chaotic snap flows," Chos, Solitons and Fractals, vol. 44, pp. 995-1003, 2011.

[21] O. E. Rossler, "An equation for hyperchaos," Physics Letters A, vol. 71, no. 2-3, pp. 155-157, 1979.

[22] J. Stefan, "On hyperjerky systems," Chos, Solitons and Fractals, vol. 37, pp. 741-747, 2008.

[23] G. A. Leonov, N. V. Kuznetsov, and T. N. Mokaev, "Homoclinic orbits, and self-excited and hidden attractors in a Lorenz-like system describing convective fluid motion," The European Physical Journal Special Topics, vol. 224, no. 8, pp. 1421-1458, 2015.

[24] Q. Jia, "Hyperchaos generated from the Lorenz chaotic system and its control," Physics Letters A, vol. 366, no. 3, pp. 217-222, 2007.

[25] Z. Wei, I. Moroz, J. C. Sprott, A. Akgul, and W. Zhang, "Hidden hyperchaos and electronic circuit application in a 5D self-exciting homopolar disc dynamo," Chaos, vol. 27, 2017.

[26] Z. Wei, K. Rajagopal, W. Zhang, S. T. Kingni, and A. Akgul, "Synchronisation, electronic circuit implementation, and fractional-order analysis of 5D ordinary differential equations with hidden hyperchaotic attractors," Pramana--Journal of Physics, vol. 90, no. 4, 2018.

[27] V.-T. Pham, S. Jafari, C. Volos, A. Giakoumis, S. Vaidyanathan, and T. Kapitaniak, "A chaotic system with equilibria located on the rounded square loop and its circuit implementation," IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 63, no. 9, pp. 878-882, 2016.

[28] M. P. Hanias, G. Giannaris, A. Spyridakis, and A. Rigas, "Time series analysis in chaotic diode resonator circuit," Chaos, Solitons & Fractals, vol. 27, no. 2, pp. 569-573, 2006.

[29] Q. Lai, T. Nestor, J. Kengne, andX.-W. Zhao, "Coexisting attractors and circuit implementation of a new 4D chaotic system with two equilibria," Chaos, Solitons & Fractals, vol. 107, pp. 92-102, 2018.

[30] X. Wang, S. Vaidyanathan, C. Volos, V.-T. Pham, and T. Kapitaniak, "Dynamics, circuit realization, control and synchronization of a hyperchaotic hyperjerk system with coexisting attractors," Nonlinear Dynamics, vol. 89, no. 3, pp. 1673-1687, 2017.

[31] A. N. Pisarchik, R. Jaimes-Reategui, J. R. Villalobos-Salazar, J. H. Garcia-Lopez, and S. Boccaletti, "Synchronization of chaotic systems with coexisting attractors," Physical Review Letters, vol. 96, no. 24, Article ID 244102, 2006.

[32] A. N. Pisarchik and B. K. Goswami, "Annihilation of one of the coexisting attractors in a bistable system," Physical Review Letters, vol. 84, no. 7, pp. 1423-1426, 2000.

[33] J. Kengne, V. R. Signing, J. C. Chedjou, and G. D. Leutcho, "Nonlinear behavior of a novel chaotic jerk system: antimonotonicity, crises, and multiple coexisting attractors," International Journal of Dynamics and Control, vol. 6, no. 2, pp. 318-328, 2017.

[34] V.-T. Pham, S. Vaidyanathan, C. Volos, X. Wang, and D. V. Hoang, "A hyperjerk memristive system with hidden attractors," Studies in Computational Intelligence, vol. 701, pp. 59-80, 2017.

[35] A. P. Kuznetsov, S. P. Kuznetsov, E. Mosekilde, and N. V. Stankevich, "Co-existing hidden attractors in a radio-physical oscillator system," Journal of Physics A: Mathematical and Theoretical, vol. 48, no. 12, Article ID 125101,12 pages, 2015.

[36] J. C. Sprott, "A proposed standard for the publication of new chaotic systems," International Journal of Bifurcation and Chaos, vol. 21, no. 9, pp. 2391-2394, 2011.

[37] C. Masoller, "Coexistence of attractors in a laser diode with optical feedback from a large external cavity," Physical Review A: Atomic, Molecular and Optical Physics, vol. 50, no. 3, pp. 2569-2578, 1994.

[38] A. Massoudi, M. G. Mahjani, and M. Jafarian, "Multiple attractors in Koper-Gaspard model of electrochemical periodic and chaotic oscillations," Journal of Electroanalytical Chemistry, vol. 647, no. 1, pp. 74-86, 2010.

[39] C. Hens, S. K. Dana, and U. Feudel, "Extreme multistability: attractor manipulation and robustness," Chaos, vol. 25, Article ID 053112, 2015.

[40] X. Luo, M. Small, M.-F. Danca, and G. Chen, "On a dynamical system with multiple chaotic attractors," International Journal of Bifurcation and Chaos, vol. 17, no. 9, pp. 3235-3251,2007.

[41] A. N. Pisarchik and U. Feudel, "Control of multistability," Physics Reports, vol. 540, no. 4, pp. 167-218, 2014.

[42] N. V. Kurnetsov, G. A. Leonov, M. V. Yuldashev, and R. V. Yuldashev, "Hidden attractors in dynamical models of phase-locked loop circuits: limitations of simulation in MATLAB and SPICE," Communications in Nonlinear Science and Numerical Simulation, vol. 51, pp. 39-49, 2017.

N. Tsafack (iD) (1,2) and J. Kengne (1)

(1) Research Unit of Laboratory of Automation and Applied Computer (LAIA), Electrical Engineering Department of IUT-FV, University of Dschang, P.O. Box 134, Bandjoun, Cameroon

(2) Research Unit of Laboratory of Condensed Matter, Electronics and Signal Processing (URMACETS) Department of Physics, Faculty of Sciences, University of Dschang, P.O. Box 67, Dschang, Cameroon

Correspondence should be addressed to N. Tsafack; nestor.tsafack@yahoo.fr

Received 2 November 2017; Revised 30 June 2018; Accepted 24 July 2018; Published 3 September 2018

Academic Editor: Jacek Karwowski

Caption: Figure 1: Electronic circuit of the proposed 5-D hyperjerk system.

Caption: Figure 2: Backward continuation of system (3) when decreasing [a.sub.2] from 3 to 2.4 (a) and the corresponding graph of largest Lyapunov exponent ([[lambda].sub.max]) plotted in the range 2.4 [less than or equal to] [a.sub.2] [less than or equal to] 3 (b).

Caption: Figure 3: Numerical phase space trajectories (left) and Pspice based simulation results (right) showing the classical period doubling routes to chaos in the novel 5-D system.

Caption: Figure 4: Bifurcation diagrams showing local maxima of the coordinate [x.sub.3] of the attractor in Poincare cross section in terms of the control parameter [a.sub.2] for remerging Feigenbaum tree (bubbling): period-2 bubble for [a.sub.3] = 3, period-4 bubble for [a.sub.3] = 4, period-8 bubble for [a.sub.3] = 4.2, and full Feigenbaum remerging tree at [a.sub.3] = 4.5.

Caption: Figure 5: Bifurcation diagram for illustrating the coexistence of disconnected chaotic attractors with a pair of period-2 limit cycle. The diagram is plotted by forward or backward continuation of parameter b with the following initial conditions [x.sub.1](0) = 3; [x.sub.2](0) = [x.sub.3](0) = [x.sub.4](0) = [x.sub.5](0) = 0 and [x.sub.1](0) = 1; [x.sub.2](0) = [x.sub.3](0) = [x.sub.4](0) = [x.sub.5](0) = 0.

Caption: Figure 6: Two-dimensional projections ([x.sub.4]-[x.sub.5]) of four coexisting chaotic and periodic attractors for b = 29.57, [a.sub.0] = [a.sub.2] = [a.sub.3], [a.sub.2] = 7, [a.sub.4] = 5.4433 x [10.sup.-4]. Initial conditions are indicated in Table 1.

Caption: Figure 7: Bifurcation diagram for justifying the coexistence of four, five, and six different attractors in the phase space. The diagram in blue and red is plotted by forward and backward continuation while the diagram in green is plotted by following the attractor defined at [a.sub.2] = 2.71 for [x.sub.1] (0) = 3; [x.sub.2] (0) = [x.sub.3] (0) = [x.sub.4] (0) = [x.sub.5] (0) = 0.

Caption: Figure 8: Coexistence of four different attractors (a pair of period-1 limit cycles and two different symmetric attractors) for [a.sub.2] = 2.458 with the rest of system's parameters as follows: b = [a.sub.1] = 3, [a.sub.0] = 1.75, [a.sub.3] = 1, [a.sub.4] = 0.0054. Initial states are given in Table 1.

Caption: Figure 9: Coexistence of five different period-1 limit cycles for [a.sub.2] = 3 with the rest of system's parameters as follows: b = a! = 3, [a.sub.0] = 1.75, [a.sub.3] = 1, [a.sub.4] = 0.0054. Initial states are given in Table 1.

Caption: Figure 10: Two-dimensional projections ([x.sub.1] - [x.sub.2]) of five coexisting attractors for [a.sub.2] = 2.71 (a pair of chaotic attractors, a pair of period-1 limit cycle, and a symmetric chaotic attractor) with the rest of system's parameters as follows: b = [a.sub.1] = 3, [a.sub.0] = 1.75, [a.sub.3] = 1, [a.sub.4] = 0.0054. Initial states are given in Table 1.

Caption: Figure 11: Two-dimensional projections ([x.sub.1] - [x.sub.2]) of six coexisting attractors for [a.sub.2] = 2.8 (two pairs of chaotic attractors and a pair of period-1 limit cycle) with the rest of system's parameters as follows: b = [a.sub.1] = 3, [a.sub.0] = 1.75, [a.sub.3] = 1, [a.sub.4] = 0.0054. Initial states are given in Table 1.

Caption: Figure 12: Bifurcation diagram for illustrating the coexistence of seven different attractors in the phase space. The diagrams are plotted using the same methods as in Figure 7.

Caption: Figure 13: Coexistence of seven disconnected attractors (two pairs of chaotic attractors, a pair of period-1 limit cycle, and a symmetric period-1 limit cycle) for [a.sub.2] = 3 with the rest of system's parameters as follows: b = ax = 3, [a.sub.0] = 1.75, [a.sub.3] = 1, [a.sub.4] = 0.0109. Initial states are given in Table 1.

Caption: Figure 14: Pspice simulation results showing the coexistence of four different attractors for Ra2 = 17.91k [ohm] (a pair of period-1 limit cycles, and two symmetric chaotic attractors). Initial states are indicated in Table 1.

Caption: Figure 15: Pspice simulation results showing the coexistence of five different attractors for Ra2 = 14.65k[ohm] (a pair of chaotic attractors, a pair of period-1 limit cycle, and a symmetric chaotic attractor). Initial states are indicated in Table 1.

Caption: Figure 16: Pspice simulation results showing the coexistence of six different attractors for Ra2 = 16.5k[ohm] (two pairs of chaotic attractors and a pair of period-1 limit cycle). Initial states are indicated in Table 1.
Table 1: Comparative analysis of some dynamical systems by
using the largest Lyapunov exponent ([lambda]max).

References               Dynamical system

[16]               [[??].sub.1] = a([x.sub.2] -
                [x.sub.1]) + [x.sub.4] + [x.sub.5]
                   [[??].sub.2] = [cx.sub.1] -
                  [x.sub.1][x.sub.3] - [x.sub.2]
                [[??].sub.3] = [x.sub.1][x.sub.2]
                           - [bx.sub.3]
                [[??].sub.4] = -[x.sub.1][x.sub.3]
                           + [px.sub.4]
                    [[??].sub.5] = [qx.sub.1]

[17]          [[??].sub.1] = -[x.sub.2] - [x.sub.3]
              [[??].sub.2] = [x.sub.1] + [ax.sub.2]
                           + [x.sub.4]
              [[??].sub.3] = b + [x.sub.1][x.sub.3]
             [[??].sub.4] = -[cx.sub.3] + [dx.sub.4]

[23]         [[??].sub.1] = a([x.sub.2] - [x.sub.1])
                           + [x.sub.4]
              [[??].sub.2] = [cx.sub.1] - [x.sub.1]
                      [x.sub.3] - [x.sub.2]
                [[??].sub.3] = [x.sub.1][x.sub.2]
                           - [bx.sub.3]
                [[??].sub.4] = [x.sub.1][x.sub.3]
                           + [dx.sub.4]

This                 [[??].sub.1] = [x.sub.2]
work                 [[??].sub.2] = [x.sub.3]

                     [[??].sub.3] = [x.sub.4]

                    [[??].sub.4] = [bx.sub.5]

                [[??].sub.5] = -[a.sub.0][x.sub.5]
                 - [a.sub.1][x.sub.3] - [a.sub.2]
                  [x.sub.2] - [a.sub.3][x.sub.1]
                    - [gamma] sinh (x.sub.4])

                                        Maximum Lyapunov
References       Parameters        exponent [[lambda].sub.max]

[16]               a = 10          [[lambda].sub.ma] = 0.4195
                  b = 8/3
              c = 28, p = 1.3
                  q = 2.5

[17]          a = 0.25, b = 3       [[lambda].sub.ma] = 0.112
             c = 0.5, d = 0.05

[23]           a = 10, r = 28      [[lambda].sub.ma] = 0.39854
              b = 8/3, d = 1.3

This              b := 32           [[lambda].sub.ma] = 0.82
work          [a.sub.0] = 5,
               [a.sub.1] = 7
               [a.sub.2] = 5,
              [a.sub.3] = 4.50
             [a.sub.4] = 5.4433
               x [10.sup.-4]

Table 2: Initial conditions for the abundant coexisting attractors.

                                          Corresponding
             Type of     Dimensionless     electronics
Figures    coexistence      parameter      components

6             Four          b = 29.18     [R.sub.b] =
           disconnected                      342.7
             chaotic                         [ohm]
            attractors

              Four
8 & 14     disconnected    [a.sub.2]      [R.sub.a2] =
             chaotic         = 2.458          17.91
           and periodic                     k[ohm]
            attractors

              Five
           disconnected                   [R.sub.a2] =
10           chaotic       [a.sub.2]          14.65
           and periodic      = 2.71         k[ohm]
            attractors

              Five                        [R.sub.a2] =
9 & 15    disconnected     [a.sub.2]          15.6
            periodic           = 3          k[ohm]
            attractors

               Six
11 & 16    disconnected    [a.sub.2]      [R.sub.a2] =
             chaotic          = 2.8           16.5
           and periodic                     k[ohm]
            attractors

              Seven
13        disconnected     [a.sub.2]      [R.sub.a2] =
           chaotic and       = 2.834          35.28
            periodic                        k[ohm]
            attractors

                                  Numerical
                                    initial
                                  conditions
                                 ([x.sub.10],
                                 [x.sub.20],
                                 [x.sub.30],
             Type of             [x.sub.40],
Figures    coexistence           [x.sub.50])

6             Four              (a) (1,0,0,0,0)
           disconnected         (b) (3,0,0,0,0)
             chaotic
            attractors

              Four         (a) ([+ or -] 6,0,0,0,0)
8 & 14     disconnected         (b) (2,0,0,0,0)
             chaotic            (c) (3,0,0,0,0)
           and periodic
            attractors

              Five
           disconnected    (a) ([+ or -] 5,0,0,0,0)
10           chaotic      (b) ([+ or -] 0.6,0,0,0,0)
           and periodic        (c) (5.6,0,0,0,0)
            attractors

              Five         (a) ([+ or -] 5,0,0,0,0)
9 & 15    disconnected    (b) ([+ or -] 0.6,0,0,0,0)
            periodic           (c) (5.6,0,0,0,0)
            attractors

               Six        (a) ([+ or -] 0.45,0,0,0,0)
11 & 16    disconnected   (b) ([+ or -] 0.6,0,0,0,0)
             chaotic      (c) ([+ or -] 0.35,0,0,0,0)
           and periodic
            attractors

              Seven              (a) (1,0,0,0)
13        disconnected      (b) ([+ or -] 3,0,0,0)
           chaotic and      (c) ([+ or -] 2,0,0,0)
            periodic             (d) (4,0,0,0)
            attractors

                                      Pspice
                                      initial
                                    conditions
                                   ([v.sub.10],
                                    [v.sub.20],
                                    [v.sub.30],
             Type of                [v.sub.40],
Figures    coexistence              [v.sub.50])

6             Four                       --
           disconnected
             chaotic
            attractors

              Four          (a-b) ([+ or -] 0.3,0,0,0,0)
8 & 14     disconnected      (c-d) ([+ or -] 1,0,0,0,0)
             chaotic             (e) (0,0,0.2,0,0)
           and periodic
            attractors

              Five
           disconnected   (a-b) ([+ or -] 0.2,0.1,0.1,0.2)
10           chaotic      (c-d) ([+ or -] 0.5,0.1,0.1,0.2)
           and periodic   (c-d) ([+ or -] 0.5,0.1,0.1,0.2)
            attractors

              Five          (a-b) ([+ or -] 0.1,0,0,0,0)
9 & 15    disconnected      (c-d) ([+ or -] 0.2,0,0,0,0)
            periodic        (e-f) ([+ or -] 0.3,0,0,0,0)
            attractors

               Six          (a-b) ([+ or -] 0.1,0,0,0,0)
11 & 16    disconnected     (c-d) ([+ or -] 0.2,0,0,0,0)
             chaotic        (e-f) ([+ or -] 0.3,0,0,0,0)
           and periodic
            attractors

              Seven
13        disconnected                   --
           chaotic and
            periodic
            attractors
COPYRIGHT 2018 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2018 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Tsafack, N.; Kengne, J.
Publication:The Scientific World Journal
Date:Jan 1, 2018
Words:5863
Previous Article:Growth Responses and Photosynthetic Indices of Bamboo Plant (Indocalamus latifolius) under Heavy Metal Stress.
Next Article:Distribution and Seasonal Activity of Phlebotominae Sand Flies in Yazd and Its Outskirts, Center of Iran.

Terms of use | Privacy policy | Copyright © 2020 Farlex, Inc. | Feedback | For webmasters