Abstract

There are many recent investigations on chaotic hidden attractors although hyperchaotic hidden attractor systems and their relationships have been less investigated. In this paper, we introduce a hyperchaotic system which can change between hidden attractor and self-excited attractor depending on the values of parameters. Dynamic properties of these systems are investigated. Fractional order models of these systems are derived and their bifurcation with fractional orders is discussed. Field programmable gate array (FPGA) implementations of the systems with their power and resource utilization are presented.

1. Introduction

Many recent works on dynamical systems are categorized into self-excited and hidden attractors [13]. A self-excited attractor has a basin of attraction that is associated with an unstable equilibrium, while a hidden attractor has a basin of attraction that does not intersect with small neighborhoods of any equilibrium points. Hidden attractors are important in most of the engineering problems as they allow chaotic responses [4, 5]. Control of such hidden oscillations is a big challenge because of the multistability nature of the systems [6, 7]. Chaotic attractors are with no equilibrium points [815], with only stable equilibria [1619], and with curves of equilibria [20]. Fractional order with no equilibrium systems with its FPGA implementation has also been reported recently [21, 22].

Recently many researchers have discussed fractional order calculus and its applications [2325]. Fractional order nonlinear systems with different control approaches are investigated [2628]. Fractional order memristor based with no equilibrium chaotic system is proposed by Rajagopal et al. [21, 22]. A novel fractional order with no equilibrium chaotic system is investigated by Li and Chen [29]. Cafagna and Grassi investigated a fractional order hyperchaotic system without equilibrium points [30]. Memristor based fractional order system with a capacitor and an inductor is discussed [31]. Numerical analysis and methods for simulating fractional order nonlinear system are proposed by Petras [32] and MATLAB solutions for fractional order chaotic systems are discussed by Trzaska Zdzislaw [33]. A FPGA implementation of fractional order chaotic system using approximation method is investigated recently for the first time [21, 22].

Recently Jafari et al. announced a 3D chaotic system [34] which can belong to three famous categories of hidden attractors plus systems with self-excited attractors. Motivated by this, in this paper we announce a hyperchaotic chameleon which can be self-excited or a hidden attractor depending on the values of the parameters. This system helps us to better understand the hidden chaotic flows of higher dimension.

2. Novel Chaotic System (NCS)

In this section we introduce a class of novel chaotic and hyperchaotic systems with flux controlled memristor [35, 36], derived from the hyperchaotic system [37] by including parameter which governs the equilibrium of the system and parameter controls the number of Lyapunov exponents in the system, that is, chaotic and hyperchaotic case. The novel hyperchaotic system’s dimensionless state equations arewhere is the memductance of the flux controller memristor where the flux element is defined by the fourth state with and . The parameters and are fixed at and , respectively. We investigate four different choices for the parameters and as in Table 1.

Figures 14 show the 3D phase portraits of the systems , , , and , respectively.

3. Dynamic Properties of the NCS

3.1. Equilibrium Points

The equilibrium points for the NCS can be calculated by equating the state equations to 0. It can be seen that shows two cases of equilibrium points; that is, when the system has origin as the only defined equilibrium point and when the system has no defined equilibrium and hence exhibits hidden attractors. Table 2 shows the equilibrium points for different choices of and . The characteristic equation of and is and the Eigenvalues are ; ; ; and is an unstable focus and thus the two systems are self-excited attractors. As investigated by many researchers [815], chaotic attractors with no equilibrium are hidden thus making and hidden attractors.

3.2. Lyapunov Exponents and Kaplan-Yorke Dimension

Lyapunov exponents of a nonlinear system define the convergence and divergence of the states. The existence of a positive Lyapunov exponent confirms the chaotic behavior of the system [38, 39]. Lyapunov exponents (LEs) are necessary and more convenient for detecting hyperchaos in fractional order hyperchaotic system. A definition of LEs for fractional differential systems was given in [40] based on frequency-domain approximations, but the limitations of frequency-domain approximations are highlighted by Tavazoei and Haeri [41]. Time series based LEs calculation methods like Wolf algorithm [11], Jacobian method [12], and neural network algorithm [13] are popularly known ways of calculating Lyapunov exponents for integer and fractional order systems. Hence we use the Jacobian method to calculate the LEs. Table 3 shows the Lyapunov exponents of the NCS.

3.3. Bifurcation

In this section we derive the bifurcation contours for the NCS. The chaotic behavior of the system largely depends on the parameters , , , and . As discussed in [42], the transient behaviors occurring in memristor based nonlinear systems may result in longer simulation times to reach steady states. Hence we used the ode45 solver for numerical simulations. Four different cases of bifurcations are investigated for the NCS. In case 1 the parameter is varied and the bifurcation of the attractor is investigated. Figure 5 shows the bifurcation plot for . In the second case parameter is varied and the bifurcation of the NCS is studied which is shown in Figure 6. In case 3 the parameter is varied and bifurcation analysis of the system is investigated. Finally bifurcation plot for parameter is derived in case 4. Figures 7 and 8 show the bifurcation contours of and , respectively. As can be seen from the figures the NCS shows strange attractor, hyperchaos, chaos, and quasiperiodic systems. For , , , and the NCS shows strange attractor and the transient behavior of the memristor prevents the bifurcation plots to show the period doubling property even after reducing the transients to 70%. The system exhibits hyperchaotic attractor with hidden oscillations for , , , and . For a small band of with , , and the system shows chaotic oscillations with hidden attractors. A quasichaotic system is seen for , , , and .

3.4. Bicoherence

The motivation to study the bicoherence is twofold. First, the bicoherence can be used to extract information due to deviations from Gaussianity and suppress additive (colored) Gaussian noise. Second, the bicoherence can be used to detect and characterize asymmetric nonlinearity in signals via quadratic phase coupling or identify systems with quadratic nonlinearity. The bicoherence is the third-order spectrum. Whereas the power spectrum is a second-order statistic, formed from , where is the Fourier transform of , the bispectrum is a third-order statistic formed from . The bispectrum is therefore a function of a pair of frequencies . It is also a complex-valued function. The (normalized) square amplitude is called the bicoherence (by analogy with the coherence from the cross-spectrum). The bispectrum is calculated by dividing the time series into segments of length _seg, calculating their Fourier transforms and biperiodogram and then averaging over the ensemble. Although the bicoherence is a function of two frequencies the default output of this function is a one-dimensional output, the bicoherence refined as a function of only the sum of the two frequencies. The autobispectrum of a chaotic system is given by Pezeshki et al. [43]. They derived the autobispectrum with the Fourier coefficients.where is the radian frequency and is the Fourier coefficients of the time series. The normalized magnitude spectrum of the bispectrum known as the squared bicoherence is given bywhere and are the power spectrums at and .

Figures 9 and 10 show the bicoherence contours of the NCS for and and and , respectively. As can be seen from the figures, the NCS shows wider band of power spectrum when and because of the hyperchaotic hidden attractor as when the memristor element introduces a quadratic nonlinearity resulting in the cross-bicoherence. Shades in yellow represent the multifrequency components contributing to the power spectrum. From Figures 9 and 10 the cross-bicoherence is significantly nonzero and nonconstant, indicating a nonlinear relationship between the states. The yellow shades and nonsharpness of the peaks, as well as the presence of structure around the origin in figures (cross-bicoherence), indicate that the nonlinearity between the states , , , and is not of the quadratic nonlinearity and hence may be because of nonlinearity of higher dimensions. The most two dominant frequencies () are taken for deriving the contour of bicoherence. The sampling frequency () is taken as the reference frequency. Direct FFT is used to derive the power spectrum for individual frequencies and Hankel operator is used as the frequency mask. Hanning window is used as the FIR filter to separate the frequencies.

4. Fractional Order NCS (FONCS)

In this section we derive the fractional order model novel chaotic system (FONCS). There are three commonly used definitions of the fractional order differential operator, namely, those of Grunwald–Letnikov, Riemann–Liouville, and Caputo [2325].

In this section, we will study the dynamical behavior of fractional order system derived from the NCS with the Grunwald-Letnikov (GL) definition, which is defined aswhere and are limits of the fractional order equation, is generalized difference, is the step size, and is the fractional order of the differential equation.

For numerical calculations the above equation is modified asTheoretically fractional order differential equations use infinite memory. Hence when we want to numerically calculate or simulate the fractional order equations we have to use finite memory principal, where is the memory length and is the time sampling.The binomial coefficients required for the numerical simulation are calculated asUsing (4)–(6) the fractional order NCS is defined asAs discussed in Section 1, the FONCS also shows chaotic and hyperchaotic systems with no equilibrium and single equilibrium points for a choice of parameter values and as shown in Table 4.

The 3D phase portraits of the FONCS are shown in Figures 1114. The commensurate fractional orders of the system for , , , , , , , and are taken as , , , and , respectively.

4.1. Dynamic Analysis of the FONCS

Most of the dynamic properties of the NCS like the Lyapunov exponents and bifurcation with parameters are preserved in the FONCS [44] if , where The most important analysis of interest when investigating a fractional order system is the bifurcation with fractional order. The largest positive Lyapunov exponents () of the NCS for and appear when against their largest integer order Lyapunov exponents (), largest positive Lyapunov exponents () of the NCS for and appear when against their largest integer order Lyapunov exponents (), largest positive Lyapunov exponents () of the NCS for and appear when against their largest integer order Lyapunov exponents () and largest positive Lyapunov exponent () of the NCS for and appear when against its largest integer order Lyapunov exponent (). It can also be seen that as the fractional order decreases, the FONCS starts losing its largest positive Lyapunov exponent. When the positive Lyapunov exponents of the system become negative and thus the chaotic oscillation in the system disappears. Figures 15 and 16 show the bifurcation of the FONCS for variation in fractional orders for the two distinctive cases and and and .

4.2. Stability Analysis of FONCS
4.2.1. Commensurate Order

For commensurate FONCS of order , the system is stable and exhibits chaotic oscillations if , where is the Jacobian matrix at the equilibrium and are the Eigenvalues of the FONCS, where . As seen from the FONCS, the Eigenvalues should remain in the unstable region and the necessary condition for the FONCS to be stable is . The Eigenvalues of and are ; ; ; and is an unstable focus contributing to the existence of chaotic oscillations.

4.2.2. Incommensurate Order

The necessary condition for the FONCS to exhibit chaotic oscillations in the incommensurate case is , where is the LCM of the fractional orders. If , , , and , then . The characteristic equation of the system evaluated at the equilibrium is and by substituting the values of and the fractional orders, , the characteristic equation is + and is the same for both and . The approximated solution of the characteristic equation is whose argument is zero and which is the minimum argument and hence the stability necessary condition becomes which solves for and hence the FONCS is stable and chaos exists in the incommensurate system.

5. FPGA Implementation of the Fractional Order Novel Cubic Nonlinear Systems

In this section we discuss the implementation of the proposed FONCS in FPGA [21, 22, 38, 4548] using the Xilinx (Vivado) System Generator toolbox in Simulink. The challenge while implementing the systems in FPGAs is to design the fractional order integrator which is not a readily available block in the System Generator [21, 22]. Hence we implement the fractional integrators using the mathematical relation [32, 33] discussed in (4), (5), and (6) and the value of is taken as 0.001 with the initial conditions as described in Table 1 and the commensurate fractional order taken as for FONCS-1, for FONCS-2, for FONCS-3, and for FONCS-4. Figures 17, 18(a), and 18(b) show the Xilinx RTL schematics of the FONCS-1 system implemented in Kintex-7 (device = 7k160t, package = fbg484 S), power utilized by the system, and power utilized for various fractional orders, respectively. Table 5 shows the resources utilized by the FONCS-1 system including the clock frequency. Figure 19 shows the 2D state portraits of the FPGA implemented FONCS-1 system. Figures 20, 21(a), and 21(b) show the Xilinx schematics of the FONCS-2 system implemented in Kintex-7 (Device = 7k160t Package = fbg484 S), power utilized by the system and power utilization for various fractional orders, respectively. Table 6 shows the resources utilized by the FONCS-2 system including the clock frequency. Figure 22 shows the 2D state portraits of the FPGA implemented FONCS-2 system. Figures 23, 24(a), and 24(b) show the Xilinx schematics of the FONCS-3 system implemented in Kintex-7 (device = 7k160t, package = fbg484 S), power utilized by the system, and power utilization for various fractional orders, respectively. Table 7 shows the resources utilized by the FONCS-3 system including the clock frequency. Figure 25 shows the 2D state portraits of the FPGA implemented FONCS-3 system. Figures 26, 27(a), and 27(b) show the Xilinx schematics of the FONCS-4 system implemented in Kintex-7 (device = 7k160t, package = fbg484 S), power utilized by the system, and power utilization for various fractional orders, respectively. Table 8 shows the resources utilized by the FONCS-3 system including the clock frequency. Figure 28 shows the 2D state portraits of the FPGA implemented FONCS-3 system. The sampling rates of the FPGA blocks play a crucial role in the existence of Lyapunov exponents and also increasing the sampling time period in some implementations may lead to a clock frequency mismatch. Maximum power will be consumed by the system when the FONCS shows largest Lyapunov exponents (FONCS-1 , FONCS-2 , and FONCS-3 ). To utilize the power of FPGA, the computation needs to be divided into several independent blocks of threads that can be executed simultaneously [39, 49]. The FONCS power efficiency also depends on the parameter as can be observed from the power efficiency Figures 18, 21, 24, and 27. When the FONCS uses a maximum power of 299 w against 280 w when . This is because of the reason that the system shows their largest positive Lyapunov exponents when . The effect of parameter on the power efficiency is quiet minimal and negligible. The performance on FPGA is directly related to the number of threads and its performances and hence the FONCS are designed as four parallel threads. The fractional order operators are implemented as building blocks and the so-called “frame delay” is not noticeable in the FPGA hardware implementation due to its parallel data structure, unlike a microprocessor-based implementation.

6. Conclusion

This paper introduces a new hyperchaotic system which when changing the values of parameters exhibits self-excited and hidden attractors. Dynamic analyses of the proposed hyperchaotic system are investigated. Fractional order model of the hyperchaotic system is derived and implemented in FPGA. Power efficiency analyses for various fractional orders are derived and it is shown that the system uses maximum power when exhibiting its largest Lyapunov exponent.

Conflicts of Interest

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

Acknowledgments

The authors wish to extend their sincere thanks and gratitude to Dr. Sajad Jafari, Department of Biomedical Engineering, Amirkabir University of Technology, for guiding in achieving crucial results.