# Autonomous Jerk Oscillator with Cosine Hyperbolic Nonlinearity: Analysis, FPGA Implementation, and Synchronization.

1. IntroductionThere has been a significant increase in studies of chaotic oscillators over the past decades because of their complex behaviors and their promising applications [1-9]. There are increasing works on chaos in jerk oscillators [10-12]. It is worth noting that if the scalar x(t) is denoted as a physical variable at time t, the third derivative [??](t) presents the jerk [13]. Different jerk oscillators with chaos were summarized by Sprott [14]. By using Josephson junctions, Yalcin constructed a general jerk circuit with multiscroll and hypercube attractors [15]. Multiscroll attractors in jerk circuits were presented by Ma et al. [16]. An elegant chaotic oscillators based on jerk equation was implemented in a circuit where the nonlinearity was provided by a single diode [17]. A simple chaotic jerk circuit was used in a sound encryption scheme [18]. Especially, a three-dimensional novel jerk chaotic oscillator with two hyperbolic sinusoidal nonlinearities was reported in [19].

Researchers have shown an increased interest in multistability [20-26]. Multistability is associated with the presence of coexisting attractors with different initial conditions [20], Jerk oscillators can exhibit multistability despite their simplicity [27-31], Recently, Njitacke et al. discovered the coexistence of multiple attractors and crisis route to chaos in a chaotic jerk circuit [27] and a simple autonomous jerk oscillator with multiple attractors was introduced in [28], Kengne et al. performed an analysis of a jerk circuit including a pair of semiconductor diodes connected in antiparallel and investigated its multiple attractors [29]. It is interesting that, by applying the diode bridge memristor into the original jerk circuit, authors proposed a novel hybrid diode based jerk circuit [30]. This work revealed that the hybrid diode based jerk circuit exhibits rich dynamical behaviors including multiple coexisting self-excited attractors. Derived from the autonomous jerk circuit, an autonomous memristor-based jerk circuit was constructed [31]. Interestingly, for the same values of system parameters, the coexistence of four different attractors was obtained in such a memristor-based jerk circuit.

Motived by published results related to jerk oscillators, some questions arose as to know, e.g., if a jerk oscillator with cosine hyperbolic nonlinearity can exhibit multistability or how such oscillator can synchronize. The aim of our work is to explore aspects of the unanswered questions. Our paper is structured as follows. Section 2 is devoted to the theoretical analysis of proposed autonomous jerk oscillator with a cosine hyperbolic term. FPGA implementation of proposed autonomous jerk oscillator is presented in Section 3. In Section 4, synchronization of unidirectional coupled jerk oscillators is studied by applying adaptive sliding mode control method. Finally, Section 5 concludes our work.

2. Theoretical Analysis of Proposed Autonomous Jerk Oscillator

In chapter three of "Elegant Chaos: Algebraically Simple flow" book published in 2010, Sprott proposed a list of sixteen autonomous chaotic jerk oscillators with different nonlinearities called memory oscillators [14]. The nonlinearities of memory oscillators include quadratic, cubic, quintic, absolute, maximum, sign, exponential, sine, and tangent hyperbolic functions. Inspired by [14], in this work we introduce an autonomous jerk oscillator with a cosine hyperbolic nonlinearity described by

dx/dt = z, (1a)

dy/dt = ax - y - bz - cosh (x), (1b)

dz/dt = y, (1c)

where x,y,z are state variables of the oscillator and a,b two positive parameters. System (1a), (1b), and (1c) can be converted to a jerk oscillator, as follows:

[d.sup.3]x/d[t.sup.3] = -[d.sup.2]x/d[t.sup.2] - b dx/dt + ax - cosh(x). (2)

System (1a), (1b), and (1c) is dissipative because [[nabla]V = [partial derivative](dx/dt)/[partial derivative]x + [partial derivative](dy/dt)/[partial derivative]y + [partial derivative](dz/dt)/[partial derivative]z = -1. The equilibrium points of system (1a), (1b), and (1c) are obtained by solving dx/dt = 0, dy/dt = 0, dz/dt = 0, which gives

[y.sup.*] = [z.sup.*] = 0, (3a)

a[x.sup.*] - cosh ([x.sup.*]) = 0. (3b)

Equation (3b) cannot be solved analytically. We thus use the Newton-Raphson method to find the value of [x.sup.*]. Depending on the value of parameter a, (3b) presents no roots for 0 [less than or equal to] a < 1.509 and one root ([x.sup.*]) for a [greater than or equal to] 1.509. Therefore, for 0 [less than or equal to] a < 1.509, system (1a), (1b), and (1c) has no equilibrium points while, for a [greater than or equal to] 1.509, it has two equilibrium points [E.sub.1] = ([x.sup.*.sub.1], 0,0) and [E.sub.2] = ([x.sup.*.sub.2],0,0) with [x.sup.*.sub.1] > [x.sup.*.sub.2] (see Figure 1). The characteristic equation evaluated at the equilibrium point E = ([x.sup.*] ,0,0) is

[[lambda].sup.3] + [[lambda].sup.2] + b[lambda] - a + sinh ([x.sup.*]) = 0. (4)

According to the Routh-Hurwitz criterion, the real parts of all the roots [lambda] of (4) are negative if and only if

-a + sinh ([x.sup.*]) > 0 (4b)

a + b - sinh ([x.sup.*]) > 0 (4c)

The stability analysis of equilibrium point [E.sub.1] as function of the parameter a of system (1a), (1b), and (1c) is depicted in Figure 1.

From Figure 1, one can notice that equilibrium point [E.sub.2] is always unstable for 1.65 [less than or equal to] a [less than or equal to] 4, whereas the equilibrium point [E.sub.1] is stable for 1.65 [less than or equal to] a [less than or equal to] 2.269 and unstable for a > 2.269. Since the equilibrium points [E.sub.1] is stable for 1.65 [less than or equal to] a [less than or equal to] 2.269 but it changes stability properties at a [approximately equal to] 2.269, system (1a), (1b), and (1c) can exhibit a Hopf bifurcation from equilibrium point [E.sub.1] when the parameter a varies.

Theorem 1. System (1a), (1b), and (1c) undergoes a Hopf bifurcation at the equilibrium point [E.sub.1] if b > 0 and the parameter a passes through the critical value [a.sub.H] = -b + sinh([x.sup.*.sub.1]).

Proof. Substituting [lambda] = i[omega]([omega] > 0) into (4) and separating real and imaginary parts, we obtain

[omega] = [[omega].sub.0] = [square root of b], (5a)

a = [a.sub.H] = -b + sinh ([x.sup.*.sub.1]). (5b)

Differentiating both sides of (4) with respect to a, we can obtain

3[[lambda].sup.2] d[lambda]/da + 2[lambda] d[lambda]/da + b d[lambda]/da - 1 = 0 (6a)

and

d[lambda]/da = 1/d[[lambda].sup.2] + 2[lambda] + b (6b)

then

[mathematical expression not reproducible]. (7)

All the conditions for Hopf bifurcation to occur are met; therefore system (1a), (1b), and (1c) undergoes a Hopf bifurcation at [E.sub.1] when [a.sub.H] = -b + sinh([x.sup.*.sub.1]) and periodic solutions will exist in a neighbourhood of the point [a.sub.H] (provided that b >0 holds). If b = 3, the critical value is a = [a.sub.H] [approximately equal to] 2.269. For a = 2.2 < [a.sub.H], the trajectories of jerk system (1a), (1b), and (1c) converge to the equilibrium point [E.sub.1] whereas, for a = 2.5 > [a.sub.H], system (1a), (1b), and (1c) exhibits a limit cycle (not shown).

We plot in Figure 2 a two-parameter (a, b) bifurcation diagram depicting the dynamical behaviors of system (1a), (1b), and (1c).

Figure 2 indicated that system (1a), (1b), and (1c) can display unbounded orbits and periodic and chaotic behaviors. In Figure 3, we firstly fix b = 3 and plot the bifurcation diagram with respect to a and the related largest Lyapunov exponent.

The bifurcation diagram of the output x(t) in Figure 2(a) shows that the trajectories of system (1a), (1b), and (1c) converge to the equilibrium point E1 up to a [approximately equal to] 2.269 where a Hopf bifurcation occurs followed by a limit cycle motion (i.e., period-1 oscillations) and period-doubling to chaos interspersed with periodic windows. The largest Lyapunov exponent shown in Figure 3(b) confirms the chaotic behavior found in Figure 3(a). The chaotic behavior is illustrated in Figure 4 for a specific value of parameter a.

The trajectories of chaotic attractor are swirling around one of the two equilibrium points in Figure 4. This is a signature of one-scroll chaotic attractor.

Secondly, for a = 3.3, sample results showing the bifurcation diagram versus b and the related plot of Largest Lyapunov exponent are provided in Figure 5.

By varying parameter b from 2.87 to 3.8, the bifurcation diagram of the output x(t) in Figure 5(a) displays chaotic behavior interspersed with periodic windows. For b > 3.18, period-12-oscillations, period-6-oscillations, and period-3-oscillations are found, respectively, up to b < 3.2012 followed by reverse period-doubling to chaos interspersed with periodic windows. By varying parameter b from 3.8 to 2.87 [see red dots of Figure 5(a)], the output x(t) displays the same dynamical behaviors as shown by black dots of Figure 5(a) in the ranges 2.87 [less than or equal to] b [less than or equal to] 3.182 and 3.2012 [less than or equal to] b [less than or equal to] 3.8, while, in the range 3.182 < b < 3.2012, the output x(t) presents chaotic behavior. By comparing black dots of Figure 5(a) and red dots of Figure 5(a), one can notice that system (1a), (1b), and (1c) displays coexistence of period-6-oscillations and chaotic attractors in the range 3.182 < b < 3.188 and coexistence of period-3-oscillations and chaotic attractors in the range 3.188 [less than or equal to] b < 3.2012. The largest Lyapunov exponent shown in Figure 5(b) confirms the chaotic behavior found in Figure 5(a). Figure 6 depicts the phase portraits of coexisting attractors found in Figure 5(a) for specific values of parameter b.

For b = 3.185, system (1a), (1b), and (1c) can exhibit either one-scroll chaotic attractor or period-6-oscillations depending on initial conditions as shown in Figure 6(a), while, in Figure 6(b) for b = 3.189, system (1a), (1b), and (1c) can exhibit either one-scroll chaotic attractor or period-3oscillations depending on initial conditions. Figure 7 presents the basin of attraction of system (1a), (1b), and (1c) in the plane 7 = 0 for a = 3.3 and b = 3.189.

In Figure 7, red, green, and black regions contain initial conditions that lead to unbounded orbits, periodic and chaotic attractors, respectively. One can see from the basin of attraction that the possibility of occurrence of unbounded orbits is greater than the ones of periodic and chaotic attractors.

3. FPGA Implementation of Proposed Autonomous Jerk Oscillator

There have been many literatures about implementation of chaotic oscillators using FPGA like FPGA based multiscroll attractor discussed in [32-39], digital chaotic oscillators and its FPGA model [36, 40], FPGA based four-wing chaotic attractor [37], and image cryptography with chaotic ciphers and its FPGA implementations and FPGA implemented memristor-based chaotic attractors [41]. Hyperchaotic oscillators with complex mathematical models are also implemented using FPGA and show that the FPGA performs better in power efficiency [42, 43]. More complex synchronization algorithms are also implemented in FPGA such as synchronization of Ikeda oscillators [44], genetically optimized PID and sliding mode control based synchronizations [45, 46], and adaptive sliding mode synchronization [47]. These FPGA realizations are mainly based on hardware-software cosimulations [45-47] and with experimental analysis [44]. FPGA are also used in implementing various chaos-based encryption schemes like chaos-based pseudorandom number generator [48]. In this section we propose to implement the proposed autonomous jerk oscillator using FPGA. The major task in this implementation is to decide the type of numerical method to solve the dynamical system. We choose the forward Euler method [40] to solve the proposed autonomous jerk oscillator and the set of discretized system equations is given as follows:

[x.sub.k+1] = [x.sub.k] + h[[z.sub.k-1]], (8a)

[y.sub.k+1] = [y.sub.k] + h[a[x.sub.k-1] - [y.sub.k-1] - b[z.sub.k-1] - cosh([x.sub.k-1])], (8b)

[z.sub.k+1] = [z.sub.k] + h[[y.sub.k-1]], (8c)

where a, b are the system parameters and h is the step size for the discrete numerical solution. To implement the hyperbolic cosine function, we use the Taylor series expansion with the first three polynomials as

cosh(x) = 1 + [x.sup.2]/2! + [x.sup.4]/4! + [x.sup.6]/6!. (9)

To achieve the accuracy we select the step size h = 0.001 and as per the IEEE754 standards, we use the calculation of 32-bit finite integer instead of floating-point operation. By using the arithmetic operations multiplying, adding or subtracting the system (8a), (8b), and (8c) is implemented and the initial conditions for the next iteration can be obtained. By repeating this process of iteration we can derive the discrete chaotic system. The absolute value operation only needs to set the first bit of the 32-bit integer to 0. The discrete state equations are implemented using hardware and software cosimulations and the needed basic arithmetic operators are implemented using the Xilinx system generator toolkit. We used the Kintex 7 chipset (xc7k160tfbg)) for the cosimulations and Matlab Simulink is used to plot the phase plots. The latency of the arithmetic blocks is kept at 3 and the maximum clock frequency of the FPGA used is 437 MHz. The register-transfer level (RTL) schematic of hyperbolic cosine function given by (9) is not shown while the RTL schematic of proposed autonomous jerk oscillator is shown in Figure 8.

Table 1 shows the resources utilized by the proposed autonomous jerk oscillator for implementing in Kintex-7.

Figure 9 shows the 2D phase portraits of the proposed autonomous jerk oscillator.

Figure 10 shows the coexisting attractors exhibited by the FPGA implemented proposed autonomous jerk oscillator for two different initial conditions.

One-scroll chaotic attractor and coexistence of attractors are clearly seen from Figures 9 and 10, respectively. By comparing Figures 9 and 10 with Figures 9 and 10, it can be noticed that there is a good qualitative agreement between the numerical simulations and FPGA results.

4. Synchronization of Unidirectional Coupled Proposed Autonomous Jerk Oscillators Using Adaptive Sliding Mode Control Method

Here we present synchronization results of coupled chaotic proposed jerk oscillators in order to promote chaos-based synchronization designs of this type of proposed jerk oscillator. Synchronization of chaotic systems has applications in secure communication and cryptography. The highly sensitive nature of chaotic oscillators to initial conditions makes it difficult to synchronize the oscillators with uncertainties and disturbance. Some well-known ways to synchronize chaotic oscillators are using active control method [49, 50] adaptive control method [51, 52], extended back stepping control method [53, 54], sliding mode control method [55, 56], adaptive sliding mode control method [56-59], etc. Sliding mode control method is applied to provide robustness in the face of internal and external disturbances. In section, we use adaptive sliding mode control method for synchronization of unidirectional coupled identical proposed autonomous jerk oscillators using adaptive sliding mode control method.

4.1. Problem Statement. Let the generalized master system be defined as

[??] = f(x) + F(x)c (10)

and the slave system as

[??] = g(y) + G(y)d + u(t) (11)

where f(x),g(y) are n x 1 row vector and F(x),G(y) are m x n matrix elements of the master and slave systems, respectively, c, d are the unknown parameters of the systems, and u(t) is the controller to synchronize the systems. The control objective is to synchronize the slave system with initial condition y(0) with master system of initial condition x(0) such that the synchronization errors (12) approaches zero.

[mathematical expression not reproducible]. (12)

The sliding mode controller design for synchronizing the two systems involves selection of sliding surface for the desired dynamics and designing the reaching law such that any point on the phase space is brought to the sliding surface in the presence of uncertainties.

Let us define the proportional integral sliding surface [60] as

[s.sub.i] = [e.sub.i] + K [integral] [e.sub.i] ([tau]) d[tau] (13)

where K is the proportional constant vector [[k.sub.1],[k.sub.2],[k.sub.3]].

The first derivative of the sliding surface is derived as

[[??].sub.i] = [[??].sub.i] + K[e.sub.i]. (14)

For the existence of the sliding mode, it is necessary and sufficient that the sliding surface and its first derivative should be equal to zero. The error dynamics can be derived as

[??] = g([x.sub.s]) + G([x.sub.s]) d + u(t)-f ([x.sub.m]) - F ([x.sub.m]) c. (15)

In order to avoid the chattering phenomenon caused by discontinuous control signals, the adaptive controller used to synchronize the master and slave system is chosen as

[mathematical expression not reproducible] (16)

where k, [eta], p are positive gain values, [??], [??] are parameter estimates of master and slave systems, and s is the sliding surface.

Using (16) in (15), the error dynamics simplifies to

[mathematical expression not reproducible]. (17)

The stability of the proposed controller can be analyzed using the Lyapunov candidate function:

[mathematical expression not reproducible]. (18)

The dynamics of the Lyapunov candidate function can be derived as follows:

[mathematical expression not reproducible] (19)

Using (17), (16), and (14) in (19),

[mathematical expression not reproducible]. (20)

Let us define the parameter estimate laws as

[??] = [K.sub.c]s.G ([x.sub.s]), (21a)

[??] = [K.sub.d]s.F ([x.sub.m]), (21b)

where [K.sub.c], [K.sub.d] are positive constants.

Using (21a) and (21b) in (20), we can solve the Lyapunov function dynamics as follows:

[??] [less than or equal to] -[eta] [absolute value of ([s.sub.w])] - [rho][s.sup.2], (22)

where [eta] and [rho] are all positive, and [??] is negative definite.

4.2. Numerical Verifications. For numerical validation of the proposed synchronization method, we use the proposed autonomous jerk oscillators as master and slave and apply the adaptive sliding mode control to achieve the synchronization. Let us define the master system as

[[??].sub.m] = [Z.sub.m], (23a)

[[??].sub.m] = a[x.sub.m] - [y.sub.m] - b[z.sub.m] - cosh ([x.sub.m]), (23b)

[[??].sub.m] = [y.sub.m], (23c)

where a, b are the system parameters. The slave system with the adaptive sliding mode controllers ([u.sub.i]) are defined as

[[??].sub.s] = [z.sub.s] + [u.sub.x], (24a)

[[??].sub.s] = a[x.sub.s] - [y.sub.s] - b[z.sub.s] - cosh ([x.sub.s]) + [u.sub.y] (24b)

[[??].sub.s] = [y.sub.s] + [u.sub.z]. (24c)

The parameters of the slave system are assumed to be unknown with parameter estimates [??],[??]. Using the master system (23a), (23b), and (23c) and slave system (24a), (24b), and (24c) in the error dynamics can be derived as follows:

[[??].sub.x] = [e.sub.z] + [u.sub.x], (25a)

[mathematical expression not reproducible] (25b)

[[??].sub.z] = [e.sub.y] + [u.sub.z]. (25c)

Let the adaptive sliding mode controllers be chosen as follows:

[u.sub.x] = - [e.sub.z] - [[gamma].sub.x] tanh ([s.sub.x]) - [[rho].sub.x][s.sub.x] - [k.sub.x][e.sub.x], (26a)

[mathematical expression not reproducible], (26b)

[u.sub.z] = -[e.sub.y] - [[gamma].sub.z] tanh ([s.sub.z]) - [[rho].sub.z][s.sub.z] - [k.sub.z][e.sub.z], (26c)

where [[rho].sub.i] > 0, [[gamma].sub.i] > 0 are the sliding surface gains and [k.sub.i] is the controller gain for i = x,y,z

Using (17) and (18) with (13), the parameter update laws can be defined as follows:

[[??] = -[s.sub.y] [x.sub.m], (27a)

[??] = [s.sub.y][z.sub.m], (27b)

where [??],[??] are the dynamics of the parameter estimates [??],[??]. For numerical simulations, we take the initial conditions of the master system (25a), (25b), and (25c) as [-1,1,2], slave system as [2,2,4], parameter estimates as [??](0) = 2,[??](0) = 5, and the proportional constants K = [1,1,1]. Figure 11 shows the dynamics of the synchronization errors and synchronized states of the master and slave systems and estimated parameters.

When t > 0.3, the synchronization error variables converge to zero with exponentially asymptotical speed [see Figure 11(a)] and thereby guaranteeing the synchronization between the master system (23a), (23b), and (23c) and slave system (24a), (24b), and (24c) [see Figure 11(b)]. The unknown parameters of slave system (24a), (24b), and (24c) are simultaneously successfully estimated to their values [see Figure 11(c)].

5. Conclusion

In this paper, an autonomous jerk oscillator with a cosine hyperbolic nonlinearity and two parameters was proposed and studied. It was demonstrated that the Hopf bifurcation occurs near the equilibrium point as one of jerk oscillator parameter crosses the critical value. The proposed autonomous jerk oscillator exhibits periodic attractors, one-scroll chaotic attractors,, and coexistence between chaotic and periodic attractors. Period-doubling route to chaos and reverse period-doubling route to chaos was found in proposed autonomous jerk oscillator with the variation of each of its two parameters. Further, the proposed autonomous jerk oscillator was implemented using Field Programmable Gate Array in order to show that the proposed autonomous jerk oscillator is hardware realizable. Finally, synchronization of unidirectional coupled identical proposed autonomous jerk oscillators was achieved using adaptive sliding mode control method.

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

Data Availability

No data were used to support this study.

Conflicts of Interest

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

References

[1] E. N. Lorenz, "Deterministic nonperiodic flow," Journal of the Atmospheric Sciences, vol. 20, no. 2, pp. 130-141, 1963.

[2] A. Tamasevicius, A. Cenys, G. Mykolaitis, A. Namajunas, and E. Lindberg, "Hyperchaotic oscillator with gyrators," IEEE Electronics Letters, vol. 33, no. 7, pp. 542-544, 1997.

[3] M. Lakshmanan and D. V. Senthilkumar, Dynamics of Nonlinear Time-Delay Systems, Springer Series in Synergetics, Springer, Berlin, Germany, 2010.

[4] A. Buscarino, L. Fortuna, M. Frasca, and L. V. Gambuzza, "A chaotic circuit based on Hewlett-Packard memristor," Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 22, no. 2, Article ID 023136, 2012.

[5] C.-L. Li, H.-M. Li, W. Li et al., "Dynamics, implementation and stability of a chaotic system with coexistence of hyperbolic and non-hyperbolic equilibria," AEU - International Journal of Electronics and Communications, vol. 84, pp. 199-205, 2018.

[6] D. Lathrop, Physics Today, vol. 68, no. 4, pp. 54-55, 2015.

[7] B. Muthuswamy and S. Banerjee, "Designing Hardware for FPGAs," in A Route to Chaos Using FPGAs, vol. 16 of Emergence, Complexity and Computation, pp. 29-53, Springer International Publishing, Cham, 2015.

[8] V.-T. Pham, S. Jafari, C. Volos, T. Gotthans, X. Wang, and D. V. Hoang, "A chaotic system with rounded square equilibrium and with no-equilibrium," Optik - International Journal for Light and Electron Optics, vol. 130, pp. 365-371, 2017.

[9] F. Nazarimehr, K. Rajagopal, J. Kengne, S. Jafari, and V.-T. Pham, "A new four-dimensional system containing chaotic or hyperchaotic attractors with no equilibrium, a line of equilibria and unstable equilibria," Chaos, Solitons & Fractals, vol. 111, pp. 108-118, 2018.

[10] J. C. Sprott, "Some simple jerk functions," American Journal of Physics, vol. 65, no. 6, pp. 53/-543, 1997

[11] R. Eichhorn, S. J. Linz, and P. Hanggi, "Simple polynomial classes of chaotic jerky dynamics," Chaos, Solitons & Fractals, vol. 13, no. 1, pp. 1-15, 2002.

[12] J.-M. Malasoma, "What is the simplest dissipative chaotic jerk equation which is parity invariant?" Physics Letters A, vol. 264, no. 5, pp. 383-389, 2000.

[13] S. H. Schot, "Jerk: the time rate of change of acceleration," American Journal of Physics, vol. 46, no. 11, pp. 1090-1094, 1978.

[14] J. C. Sprott, Elegant Chaos: Algebraically Simple Chaotic Flows, World Scientific, 2010.

[15] M. E. Yalcin, "Multi-scroll and hypercube attractors from a general jerk circuit using Josephson junctions," Chaos, Solitons & Fractals, vol. 34, no. 5, pp. 1659-1666, 2007.

[16] J. Ma, X. Y. Wu, R. T. Chu, and L. Zhang, "Selection of multiscroll attractors in Jerk circuits and their verification using Pspice," Nonlinear Dynamics, vol. 76, no. 4, pp. 1951-1962, 2014.

[17] J. C. Sprott, "A new chaotic jerk circuit," IEEE Transactions on Circuits and Systems II: Express Briefs, vol. 58, no. 4, pp. 240-243, 2011.

[18] C. Volos, A. Akgul, V.-T. Pham, I. Stouboulos, and I. Kyprianidis, "A simple chaotic circuit with a hyperbolic sine function and its use in a sound encryption scheme," Nonlinear Dynamics, vol. 89, no. 2, pp. 1047-1061, 2017.

[19] S. Vaidyanathan, C. Volos, V.-T. Pham, K. Madhavan, and B. A. Idowu, "Adaptive backstepping control, synchronization and circuit simulation of a 3-D novel jerk chaotic system with two hyperbolic sinusoidal nonlinearities," Archives of Control Sciences, vol. 24, no. 3, pp. 375-403, 2014.

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

[21] 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.

[22] Z. T. Njitacke, J. kengne, H. B. Fotsin, A. Nguomkam Negou, and D. Tchiotsop, Coexistence of multiple attractors and crisis route to chaos in a novel memristive diode bidge-based Jerk circuit, Chaos, Solitons & Fractals, vol. 99, 2017.

[23] B.-C. Bao, Q. Xu, H. Bao, and M. Chen, "Extreme multistability in a memristive circuit," IEEE Electronics Letters, vol. 52, no. 12, pp. 1008-1010, 2016.

[24] B. Bao, T. Jiang, Q. Xu, M. Chen, H. Wu, and Y. Hu, "Coexisting infinitely many attractors in active band-pass filter-based memristive circuit," Nonlinear Dynamics, vol. 86, no. 3, pp. 1711-1723, 2016.

[25] 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.

[26] V.-T. Pham, C. Volos, S. Jafari, and T. Kapitaniak, "Coexistence of hidden chaotic attractors in a novel no-equilibrium system," Nonlinear Dynamics, vol. 87, no. 3, pp. 2001-2010, 2017.

[27] Z. T. Njitacke, J. Kengne, N. A. Nguomkam, T. M. Fouodji, and H. B. Fotsin, "Coexistence of multiple attractorsand crisis route to chaos in a novel chaotic jerk circuit," International Journal of Bifurcation and Chaos, vol. 25, pp. 1550052-1550060, 2015.

[28] J. Kengne, Z. T. Njitacke, and H. B. Fotsin, "Dynamical analysis of a simple autonomous jerk system with multiple attractors," Nonlinear Dynamics, vol. 83, no. 1-2, pp. 751-765, 2016.

[29] J. Kengne, A. Nguomkam Negou, and Z. T. Njitacke, "Antimonotonicity, chaos and multiple attractors in a novel autonomous jerk circuit," International Journal of Bifurcation and Chaos, vol. 27, no. 7,1750100, 20 pages, 2017.

[30] Z. T. Njitacke, J. kengne, and L. Kamdjeu Kengne, "Antimonotonicity, chaos and multiple coexisting attractors in a simple hybrid diode-based jerk circuit," Chaos, Solitons & Fractals, vol. 105, pp. 77-91, 2017.

[31] J. Kengne, A. N. Negou, and D. Tchiotsop, "Antimonotonicity, chaos and multiple attractors in a novel autonomous memristor-based jerk circuit," Nonlinear Dynamics, vol. 88, no. 4, pp. 2589-2608, 2017.

[32] E. Tlelo-Cuautle, A. D. Pano-Azucena, J. J. Rangel-Magdaleno, V. H. Carbajal-Gomez, and G. Rodriguez-Gomez, "Generating a 50-scroll chaotic attractor at 66 MHz by using FPGAs," Nonlinear Dynamics, vol. 85, no. 4, pp. 2143-2157, 2016.

[33] K. Rajagopal, A. Akgul, S. Jafari, A. Karthikeyan, and I. Koyuncu, "Chaotic chameleon: Dynamic analyses, circuit implementation, FPGA design and fractional-order form with basic analyses," Chaos, Solitons & Fractals, vol. 103, pp. 476-487, 2017.

[34] E. Dong, Z. Liang, S. Du, and Z. Chen, "Topological horseshoe analysis on a four-wing chaotic attractor and its FPGA implement," Nonlinear Dynamics, vol. 83, no. 1-2, pp. 623-630, 2016.

[35] S. M. Ismail, L. A. Said, A. A. Rezk et al., "Generalized fractional logistic map encryption system based on FPGA," AEU --International Journal of Electronics and Communications, vol. 80, pp. 114-126, 2017.

[36] E. Tlelo-Cuautle, J. J. Rangel-Magdaleno, A. D. Pano-Azucena, P. J. Obeso-Rodelo, and J. C. Nunez-Perez, "FPGA realization of multi-scroll chaotic oscillators," Communications in Nonlinear Science and Numerical Simulation, vol. 27, no. 1-3, pp. 66-80, 2015.

[37] Y.-M. Xu, L.-D. Wang, and S.-K. Duan, "A memristor-based chaotic system and its field programmable gate array implementation," Acta Physica Sinica, vol. 65, no. 12, Article ID 120503, 2016.

[38] E. Tlelo-Cuautle, L. G. de la Fraga, V.-T. Pham, C. Volos, S. Jafari, and A. D. J. Quintas-Valles, "Dynamics, FPGA realization and application of a chaotic system with an infinite number of equilibrium points," Nonlinear Dynamics, vol. 89, no. 2, pp. 1129-1139, 2017.

[39] Z. Wei, A. Akgul, U. E. Kocamaz, I. Moroz, and W. Zhang, Control, electronic circuit application and fractional-order analysis of hidden chaotic attractors in the self-exciting homopolar disc dynamo, Chaos, Solitons & Fractals, vol. 111, 2018.

[40] D. Valli, B. Muthuswamy, S. Banerjee et al., "Synchronization in coupled Ikeda delay systems," The European Physical Journal Special Topics, vol. 223, no. 8, pp. 1465-1479, 2014.

[41] Karthikeyan Rajagopal, Anitha Karthikeyan, and Prakash Duraisamy, "Hyperchaotic Chameleon: Fractional Order FPGA Implementation," Complexity, vol. 2017, Article ID 8979408, 16 pages, 2017.

[42] Karthikeyan Rajagopal, Laarem Guessas, Sundarapandian Vaidyanathan, Anitha Karthikeyan, and Ashokkumar Srinivasan, "Dynamical Analysis and FPGA Implementation of a Novel Hyperchaotic System and Its Synchronization Using Adaptive Sliding Mode Control and Genetically Optimized PID Control," Mathematical Problems in Engineering, vol. 2017, Article ID 7307452, 14 pages, 2017.

[43] Karthikeyan Rajagopal, Laarem Guessas, Anitha Karthikeyan, Ashokkumar Srinivasan, and Girma Adam, "Fractional Order Memristor No Equilibrium Chaotic System with Its Adaptive Sliding Mode Synchronization and Genetically Optimized Fractional Order PID Synchronization," Complexity, vol. 2017, Article ID 1892618, 19 pages, 2017.

[44] K. Rajagopal, A. Karthikeyan, and A. K. Srinivasan, "FPGA implementation of novel fractional-order chaotic systems with two equilibriums and no equilibrium and its adaptive sliding mode synchronization," Nonlinear Dynamics, vol. 87, no. 4, pp. 2281-2304, 2017.

[45] J. M. Bahi, X. Fang, C. Guyeux, and L. Larger, "FPGA design for pseudorandom number generator based on chaotic iteration used in information hiding application," Applied Mathematics & Information Sciences, vol. 7, no. 6, pp. 2175-2188, 2013.

[46] C. Chen, H. Ma, H. Chen, Y. Meng, and Q. Ding, "FPGA implementation of a UPT chaotic signal generator for image encryption," Pacific Science Review A: Natural Science and Engineering, vol. 17, no. 3, pp. 97-102, 2015.

[47] M. L. Barakat, A. G. Radwan, and K. N. Salama, "Hardware realization of chaos based block cipher for image encryption," in Proceedings of the 2011 23rd International Conference on Microelectronics, 1CM2011, Tunisia, December 2011.

[48] A. M. Atteya and A. H. Madian, "A hybrid Chaos-AES encryption algorithm and its impelmention based on FPGA," in Proceedings of the 2014 12th IEEE International New Circuits and Systems Conference, NEWCAS 2014, pp. 217-220, Canada, June 2014.

[49] S. T. Kingni, V.-T. Pham, S. Jafari, and P. Woafo, "A chaotic system with an infinite number of equilibrium points located on a line and on a hyperbola and its fractional-order form," Chaos, Solitons & Fractals, vol. 99, pp. 209-218, 2017.

[50] B. A. Idowu, U. E. Vincent, and A. N. Njah, "Synchronization of chaos in non-identical parametrically excited systems," Chaos, Solitons & Fractals, vol. 39, no. 5, pp. 2322-2331, 2009.

[51] V. Sundarapandian and R. Karthikeyan, "Anti-synchronization of Lu and Pan chaotic systems by adaptive nonlinear control," European Journal of Scientific Research, vol. 64, no. 1, pp. 94-106, 2011.

[52] V. Sundarapandian and R. Karthikeyan, "Adaptive anti-synchronization of Uncertain Tigan and Li Systems," Journal of Engineering and Applied Sciences, vol. 7, no. 1, pp. 45-52, 2012.

[53] S. S. Majidabad and H. T. Shandiz, "Discrete-time Terminal Sliding Mode Control of Chaotic Lorenz System," Journal of Control and Systems Engineering, vol. 1, no. 1, pp. 1-8, 2013.

[54] S. Vaidyanathan, "A new 3-D jerk chaotic system with two cubic nonlinearities and its adaptive backstepping control," Archives of Control Sciences. Polish Academy of Sciences. Committee of Automatic Control and Robotics, vol. 27(63), no. 3, pp. 409-439, 2017.

[55] O. S. Onma, O. I. Olusola, and A. N. Njah, "Control and Synchronization of Chaotic and Hyperchaotic Lorenz Systems via Extended Backstepping Techniques," Journal of Nonlinear Dynamics, vol. 2014, Article ID 861727, 15 pages, 2014.

[56] B. Wang, Y. Li, and D. L. Zhu, "Simplified sliding mode of a novel class of four-dimensional fractional-order chaos," International Journal of Control and Automation, vol. 8, no. 8, pp. 425-438, 2015.

[57] C. Yin, S. Dadras, S.-M. Zhong, and Y. Chen, "Control of a novel class of fractional-order chaotic systems via adaptive sliding mode control approach," Applied Mathematical Modelling, vol. 37, no. 4, pp. 2469-2483, 2013.

[58] H. Liu and J. Yang, "Sliding-mode synchronization control for uncertain fractional-order chaotic systems with time delay," Entropy, vol. 17, no. 6, pp. 4202-4214, 2015.

[59] S. Wang, Y. G. Yu, and M. Diao, "Hybrid projective synchronization of chaotic fractional order systems with different dimensions," Physica A: Statistical Mechanics and its Applications, vol. 389, no. 21, pp. 4981-4988, 2010.

[60] J. Slotine and W. Li, Applied Nonlinear Control, Prentice Hall, New Jersey, NJ, USA, 1991.

Karthikeyan Rajagopal (iD), (1) Sifeu Takougang Kingni (iD), (2) Gaetan Fautso Kuiate (iD), (3) Victor Kamdoum Tamba, (4) and Viet-Thanh Pham (iD) (5)

(1) Center for Nonlinear Dynamics, Defense University, Ethiopia

(2) Department of Mechanical, Petroleum and Gas Engineering, Faculty of Mines and Petroleum Industries, University of Maroua, P.O. Box 46, Maroua, Cameroon

(3) Department of Physics, Higher Teacher Training College, The University ofBamenda, P. O. Box 39 Bamenda, Cameroon

(4) Laboratory of Electronics and Signal Processing (LETS), Department of Physics, Faculty of Science, University of Dschang, P. O. Box 67, Dschang, Cameroon

(5) School of Electronics and Telecommunications, Hanoi University of Science and Technology, 01 Dai Co Viet, Hanoi, Vietnam

Correspondence should be addressed to Gaetan Fautso Kuiate; fautsoJcuiate@yahoo.com

Received 22 May 2018; Revised 17 July 2018; Accepted 30 July 2018; Published 17 September 2018

Academic Editor: Dimitrios J. Frantzeskakis

Caption: Figure 1: (Color online) Stability diagram of equilibrium points [E.sub.1,2] versus the parameter a for b = 3. Solid black lines indicate the stable branches and the dashed red lines the unstable branches.

Caption: Figure 2: (Color online) Two-parameter (a, b) bifurcation diagram depicting the dynamical behaviors of system (1a), (1b), and (1c). Unbounded orbits are in red; periodic oscillations are in green and chaos is in black.

Caption: Figure 3: The bifurcation diagrams depicting the local maxima (black dots) and local minima (gray dots) of x(t) (a) and the largest Lyapunov exponents (b) versus the parameter a for b = 3.

Caption: Figure 4: 2D phase portrait of system (1a), (1b), and (1c) for a = 3.3 and b = 3. The initial conditions are (x(0), y(0), z(0)) = (1, 0, 0.5).

Caption: Figure 5: (Color online) The bifurcation diagrams depicting the local maxima of x(t) (a) and the largest Lyapunov exponents (b) versus the parameter b for a = 3.3. In the graph (a), two sets of data corresponding, respectively, to increasing (black dot) and decreasing (red dot) values of b are superimposed. Largest Lyapunov exponents are obtained by scanning parameter b upwards (black line) and downwards (red line).

Caption: Figure 6: Phase portraits in the plane (x,z) for specific values of parameter b: (a) b = 3.185 and (b) b = 3.189. The initial conditions are (x(0), y(0), z(0)) = (1.0, 0, 0.5) in left column and (x(0), y(0), z(0)) = (3.0, 0, 0.5) in right column. The other parameter value is a = 3.3.

Caption: Figure 7: (Color online) Cross section of the basin of attraction of system (1a), (1b), and (1c) in the x-z-plane at y = 0 for a = 3.3 and b = 3.189. Initial conditions in the black region lead to unbounded orbits, those in the green region lead to the periodic attractor and those in the red region lead to the strange attractor.

Caption: Figure 8: (Color online) RTL schematic of proposed autonomous jerk oscillator implemented in Kintex 7.

Caption: Figure 9: (Color online) Phase portraits of the FPGA implemented of proposed autonomous jerk oscillator with hardware-software cosimulation for parameter values a = 3.3, b = 3 and initial conditions [1,0,0.5] with h = 0.001.

Caption: Figure 10: (Color online) Phase portraits in the plane (x,z) of the FPGA implemented proposed autonomous jerk oscillator for specific values of parameter b: (a) b = 3.185 and (b) b = 3.189 for the initial conditions (x(0), y(0), z(0)) = (1.0, 0, 0.5) in left column and (x(0), y(0), z(0)) = (3.0, 0, 0.5) in right column with the value of the parameter a = 3.3. The step size is taken as h = 0.001.

Caption: Figure 11: Time history of synchronization error variables (a), synchronized states (b), and estimated parameters (c).

Table 1: Resources consumed by the FPGA implemented of proposed autonomous jerk oscillator. Resource Utilization Available Utilization % LUT 880 101400 0.87 FF 522 202800 0.26 DSP 20 600 3.33 IO 97 285 34.04 BUFG 1 32 3.13

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article; Field Programmable Gate Array |
---|---|

Author: | Rajagopal, Karthikeyan; Kingni, Sifeu Takougang; Kuiate, Gaetan Fautso; Tamba, Victor Kamdoum; Pham, |

Publication: | Advances in Mathematical Physics |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 6170 |

Previous Article: | Optimal Stochastic Control Problem for General Linear Dynamical Systems in Neuroscience. |

Next Article: | Vortex Structures inside Spherical Mesoscopic Superconductor Plus Magnetic Dipole. |

Topics: |