Dynamic Analysis and Circuit Design of a Novel Hyperchaotic System with Fractional-Order Terms

A novel hyperchaotic system with fractional-order (FO) terms is designed. Its highly complex dynamics are investigated in terms of equilibrium points, Lyapunov spectrum, and attractor forms. It will be shown that the proposed system exhibits larger Lyapunov exponents than related hyperchaotic systems. Finally, to enhance its potential application, a related circuit is designed by using the MultiSIM Software. Simulation results verify the effectiveness of the suggested circuit.


Introduction
Hyperchaos was discovered by Rössler in 1979 [1] and the first hyperchaotic circuit was implemented by Matsumoto in 1986 [2].In these last years, hyperchaotic systems have gained the interest of the scientific community and new systems and circuits are proposed [3][4][5][6][7][8].This great interest can be explained by the aptitude of hyperchaotic systems to generate complex dynamics characterized by more than one positive Lyapunov exponent and attractors deployed in multiple directions.In practical applications and particularly in secure communication, chaotic synchronization has been explored by using electronic circuits, namely, Duffing circuit [9], Chua circuit [10], and Rössler circuit [11].However, for hyperchaotic circuits, many challenging problems are still pending due to their complex behaviors.
On the other hand, several researches have attempted to construct chaotic and hyperchaotic models with simple algebraic structures highly recommended for circuit design.The most famous chaotic one is the Jerk system proposed by sprott, in 1994 [12,13], which contains simple nonlinear terms.However, it is well known that most systems contain conventional nonlinear terms like piecewise linear functions [14][15][16][17], integer order polynomials [8,18], sine functions [19], time delayed functions [20], and switching functions [21].In this framework, fractional-order polynomials could be used to build complex chaotic behaviors and, to the best of our knowledge, they have not been harnessed until now.
The purpose of this paper is to build a novel hyperchaotic system with more complex dynamics than those proposed by related works.Expecting that the PWNL function with FO terms gives us more complex chaotic proprieties than the piecewise linear one, this PWNL function is constructed from absolute functions and FO polynomials.To enhance the potential application of the proposed system, its related circuit is designed afterwards with MultiSIM Software.
The rest of this paper is structured as follows.In Section 2, the mathematical model of the hyperchaotic system is proposed and its basic properties are presented.In Section 3, the dynamic analysis of the novel system is investigated by pointing out its elementary characteristics such as the Lyapunov exponents, the attractor forms, and the equilibrium points.In Section 4, the oscillator circuit of the hyperchaotic system is designed afterwards.

Mathematical Model and Basic Properties
Let consider the mathematical model of the novel hyperchaotic system with FO terms expressed by the following differential equations: with () being a nonlinear function defined as where (, , , , ℎ, , , ) are the system's parameters and (, , , ) are the state variables. is a fractional number satisfying 1 <  < 2. Since  ̸ = 1, ||   −1 will never be an indeterminate form.The nonlinear function () can be written as follows: System (1) can exhibit chaotic behavior if the general condition of dissipativity is satisfied such as As long as  + ℎ > 0, system (1) is dissipative and it converges to an attractor.Thus, when the parameters (, , , , , , ℎ, ) are equal to (0.93, 1.11, −0.11, −0.21, 6.26, 1.32, 0.001, 14) and the initial condition is equal to (1, 1, 1, 1), system (1) generates a strange attractor displayed in Figure 1.This attractor has an asymmetrical form with respect to all the principal axes characterized by two scrolls of different sizes.
The time series of the state variables  and  are described in Figures 2(a) and 2(b).These signals represent the chaotification rates of each variable.On the other hand, system (1) is sensitive to initial conditions as shown in Figures 2(c) and 2(d).Note that the variation range of the variable  is extended within [−250, 200], unlike the other variables.This point must be considered in practical applications.
Proof.Case (i) is obvious.

Lyapunov Exponents Analysis. System (1) exhibits four
Lyapunov exponents (LEs).These LEs are esteemed using the Wolf algorithm [23], as shown in Figure 3 as Since the LE spectrum has two positive Lyapunov exponents; thus system (1) is hyperchaotic. 1 is the largest positive one.This exponent increases the expansion degree of the attractor in the phase space.
In addition, the corresponding Kaplan-Yorke dimension is first Lyapunov coefficient (FLC).A positive FLC indicates the existence of a supercritical Hopf bifurcation, whereas a negative one indicates a subcritical Hopf bifurcation.In system (1), the two points obtained are supercritical Hopf bifurcations.This type of bifurcation indicates that the evolution to chaotic behavior is possible.
Some typical attractors are tabulated in Table 2 according to the parameter .

Comparative Analysis.
Referring to the survey paper [24], the first Lyapunov exponent can be one of the comparative criteria between hyperchaotic systems.Table 3 presents a comparative analysis between system (1) and two related ones, recently proposed in literature.Such a choice is based on the fact that, identical to system (1), the first comparative example contains linear piecewise functions whereas the second one is based on the jerk equation.Based on Table 3, it is clear that system (1) exhibits more complex dynamics.Thus, this confirms the highlight potential applications of noninteger order terms with respect to classical nonlinear terms.

Circuit Design
It is obvious that hardware implementation of chaotic systems is an interesting task in engineering applications, namely, for secure communications and random bits generation.Therefore, the aim of this section is to design an analog circuit that can build hyperchaotic behaviors according to system (1).

Design of the Analog Circuit with MultiSIM.
For the circuit implementation, we choose the particular case study when the system parameter  is fixed to 1.5.Thus, the proposed system will be defined by the following model: LE 3 = 0, Despite the observation of the obtained phase portraits of system (10), we deduce that the maximum value of the signal  can reach the level of 250.Thus, 250 V is a sufficiently high voltage for the common components used in the proposed circuit.Therefore, a linear transformation for system (10) is necessary to decrease the amplitude of the state variables.
Letting  = /2, V = /2,  = /2, and  = /160 and then setting the original state variables , , , and  instead of the variable , V, , and , the adjusted system becomes the following one: The amplitude of the state variables of system (10) has decreased as shown in Figure 6.Moreover, the two systems (10) and ( 12) are equivalent since the linear transformation does not change the physical properties of nonlinear systems.
To design the hyperchaotic circuit of system (12), only common electronic components are used such as resistors, capacitors, diodes, multipliers, and operational amplifiers.In fact, the nonlinear terms of system (12) should be designed first, namely, the quadratic term, the absolute function, the sign function, and the square root function.The quadratic term is implemented with the analog multiplier.The square root element is designed with two operational amplifiers as only active elements [27].The analog circuit of the square root element is provided in Figure 7.
For the theoretical study and based on [27], the second voltage source in Figure 7 should be fixed to 2.878 V.However, in experimentation applications, we have obtained the root square function by using a stabilized voltage equal to 2.9 V as shown in Figure 8.This figure describes two voltages; the first one is a positive source signal and the second one is the output signal of the square root circuit.Based on these results, the observed maximum voltages are equal to 720 mV and 900 mV (≃ √0.72 = 0.88), respectively.Thus, the square root function is correctly obtained with 2.9 V.In addition, based on MultiSIM results and experimental simulations, if the source voltage is included in [2.7 V, 3 V] then system (12) generates strange attractors.To avoid making this paper more cumbersome, details on experiments and experimental results will be soon presented in future works, confirming the MultiSIM results.The corresponding circuit equation of the hyperchaotic system can be described as According to system (12) and ( 13) and design considerations, we fixed the values of the resistances and the capacitors as Finally, the obtained circuit diagram, designed with Multi-SIM Software, is provided in Figure 9 where the multiplier is AD633 and the operator amplifier is UA741.

Simulation Results.
For the oscillator circuit, all active devices (UA741 and AD633) are powered by ±15 V. Several design considerations were taken into account to prevent degrading the hyperchaotic behavior such as the adjustment of the resistors and the capacitors for the integration operations.
The oscilloscope traces of the proposed circuit are shown in Figure 10.Comparing the different hyperchaotic attractors shown in Figures 6 and 10, a good qualitative agreement between the numerical simulations with Matlab and the electrical simulations with MultiSIM Software is observed.In fact, for MultiSIM Software, we have obtained the same attractors forms as those obtained by Matlab simulations.However, in these last attractors, some saturation effects are detected due to the operational amplifiers responses.To avoid making this paper more cumbersome, details on experiments and experimental results will be presented in future works, where saturation effects of amplifiers will be deeply analyzed.

Conclusions
In this paper, a novel hyperchaotic system is proposed by considering fractional-order polynomials.Analytical and numerical results show that this system exhibits more complex behaviors than those proposed by related works.Moreover, its analog circuit is designed and simulated with Mul-tiSIM Software.In future works, experimental realization of the hyperchaotic circuit will be proposed and the saturation effects induced by the operational amplifiers will be analyzed.Thereafter, the proposed circuit will be considered for secure image encryption and decryption applications.

Figure 8 :Figure 9 :
Figure 8: Experimental results of the square root function.

Figure 10 :
Figure 10: Simulation results of the hyperchaotic system with MultiSIM Software.

Table 2 :
The LEs of some typical attractors of system (1).

Table 3 :
Comparative analysis with related hyperchaotic systems.