Analysis , Synchronization and Circuit Design of a 4 D Hyperchaotic Hyperjerk System

In this work, a 4D hyperchaotic hyperjerk system, with better results for its Lyapunov exponents and Kaplan–Yorke dimension regarding other systems of this family, as well as its circuit implementation, is presented. Hyperchaotic hyperjerk systems depict complex dynamical behavior in a high-dimensional phase space with n ≥ 4, offering robustness against many types of attacks in private communications. For this reason, an adaptive controller in order to achieve global chaos synchronization of coupled 4D hyperchaotic hyperjerk systems with unknown parameters is designed. The adaptive results in this work are proved using Lyapunov stability theory and the effectiveness of the proposed synchronization scheme is confirmed through the simulation results.


Introduction
Nonlinear dissipative dynamical systems characterized by high sensitivity on initial conditions are known today as chaotic systems.The aforementioned feature of sensitivity, which is present in these kinds of dynamical systems, causes the trajectories of two identical chaotic systems to divert exponentially, despite the almost identical initial conditions.
It is known from literature that a hyperchaotic system is defined as a dynamical system with at least two positive Lyapunov exponents.As a consequence, the dynamics of a system like that can expand in several different directions simultaneously.It is well known that the minimum dimension for an autonomous, continuous-time, hyperchaotic system is four.Since the discovery of a first 4D hyperchaotic system by Rössler in 1979 [1], many 4D hyperchaotic systems have been reported in literature, such as the hyperchaotic Lorenz system [2], the hyperchaotic Lü system [3], the hyperchaotic Chen system [4], the hyperchaotic Wang system [5], the hyperchaotic Newton-Leipnik system [6], the hyperchaotic Jia system [7], the hyperchaotic Vaidyanathan system [8], etc.
In mechanics, if the scalar x(t) represents the position of a moving object at time t, then the first derivative, .
x(t), represents the velocity, the second derivative, ..
x(t), represents the acceleration and the third derivative, ... x (t), represents the jerk or jolt [9].In mechanics, a jerk system is described by an explicit third order ordinary differential equation describing the time evolution of a single scalar variable x according to the dynamics: A particularly simple example of a jerk system is the famous Coullet system [10] given by where g(x) is a nonlinear function such as g(x) = b(x 2 − 1), which exhibits chaos for a = 0.6 and b = 0.58.
A generalization of the jerk dynamics is given by the dynamics 1) x dt n−1 , . . ., An ordinary differential equation of the Formula (3) is called a hyperjerk system since it involves time derivatives of a jerk function [11].
Due to the importance of enabling a chaotic or hyperchaotic system to maintain a desirable dynamical behavior, many relative techniques have been proposed describing either a static or a dynamic feedback control, or an open-loop control method.In some contemporary cases, the problem that must be solved is the investigation of a state feedback control law in order to stabilize the system around its unstable equilibrium point [12,13].Some popular techniques in literature are the active control method [14,15], the adaptive control method [16][17][18], the sampled data feedback control method [19,20], the time delay feedback approach [21,22], the backstepping method [23,24], the sliding mode control method [25,26], the fuzzy-model-based control method [27] and the Ott, Grebogy and Yorke (OGY) method [12].
The study of chaos in the last decades had a tremendous impact on the foundations of science and engineering, and one of the most recent exciting developments in this regard is the discovery of chaos synchronization.Therefore, the methods used for controlling nonlinear systems can also be applied for synchronizing two or more coupled chaotic oscillators.Some of the aforementioned methods only apply to chaotic systems while others need several controllers to realize synchronization, focusing on the condition with unknown parameters and disturbances [28,29].Synchronization of coupled chaotic systems is defined as the phenomenon that occurs when a chaotic system drives another chaotic system by adjusting a given property of their motion.Due to the feature of exponential divergence of trajectories of two identical chaotic systems, which start with nearly the same initial conditions, the synchronization of chaotic systems is an interesting research topic in chaos literature.In addition, synchronized chaotic systems can then set up private communication systems with applications in secure communications [30][31][32][33][34], cryptosystems [35,36] and encryption [37,38], offering robustness against various attack methods.
Yamada and Fujisaka [39] conducted the first research on synchronizing a chaotic system, while Pecora and Carroll [40,41] discovered that chaotic systems can be synchronized when an appropriate coupling design is found, such as that the chaotic time evolutions of the two systems become identical.Other approaches regarding the synchronization of chaotic systems are either by the use of state observers or by the use of control laws according to the Lyapunov stability theory.The main types of synchronization are the complete synchronization [42,43], the phase synchronization [44,45], the generalized synchronization [46,47], the antisynchronization (AS) [48,49], the generalized projective synchronization [50] and the hybrid synchronization [51].
Especially in the case of complete synchronization, the main aim is to use the output of the master system in order to control the slave system.In this way, the output of the slave system tracks the output of the master system asymptotically with time.For the aforementioned control methods, adaptive control is the main method, which is used when some or all of the chaotic or hyperchaotic system's parameters are not available for measurement and estimates the uncertain parameters of the systems [16,52].
In this research work, a 4D hyperjerk system with hyperchaotic behavior is studied, by using some of the most used tools of nonlinear dynamics, such as bifurcation diagrams, Lyapunov exponents, phase portraits and Poincaré maps.In addition, the circuit of the proposed dynamical system, which proves its feasibility, is designed.Furthermore, an adaptive synchronization scheme for a system of coupled 4D hyperchaotic hyperjerk systems has been developed.For the obtained, with this method, main adaptive results are proved using Lyapunov stability theory.The effectiveness of the synchronization scheme, which is used in this work, has proved through the simulation results.
The rest of the paper is organized as follows.The next section provides a brief description of the 4D hyperchaotic hyperjerk system.Section 3 presents the analysis of the proposed system's dynamics.The circuit realization of the 4D hyperchaotic hyperjerk system is discussed in Section 4. The adaptive synchronization scheme for a system of two coupled 4D hyperchaotic hyperjerk systems is presented in Section 5. Finally, Section 6 outlines the main points that have been reached with this research study.

Model of the 4D Hyperjerk System
In 2006, Chlouverakis and Sprott [53] discovered a simple hyperchaotic hyperjerk system given by the dynamics which can also be expressed as the following system (5) of differential equations: x 2 = x 3 .
In this work, a 4D hyperjerk system by replacing the absolute nonlinearity (b|x 1 |) with (b|x 2 |) to system (7) and with a slightly different set of values for the system parameters is studied.Thus, the new 4D hyperjerk system is given in system form as: where a, b and c are positive parameters.In this paper, we shall show that system (8) is hyperchaotic when the parameters a, b and c take the values a = 3.8, b = 0.1 and c = 1.5.
For these parameter values, the Lyapunov exponents of the hyperjerk system (8), which are calculated in this work by using Wolf's algorithm [58] are obtained as: L 1 = 0.1909, L 2 = 0.06462, L 3 = 0 and L 4 = −1.81846.It is easily seen that the maximal Lyapunov exponent (MLE) of our novel hyperchaotic hyperjerk system (8) is L 1 = 0.1809, which is greater than the respective MLE of systems ( 5) and (7).In addition, the Kaplan-Yorke dimension of the novel hyperjerk system (8), has been calculated as: D KY = 3.1405.Thus, the Kaplan-Yorke is also greater than the respective D KY of hyperjerk systems ( 5) and (7).This shows that the proposed hyperchaotic 4D hyperjerk system (8) exhibits more complex behavior than systems ( 5) and (7).Furthermore, the proposed in this work system has better results concerning others reported in literature hyperchaotic hyperjerk systems concerning their MLE and only one (see Ref. [59]) has greater Kaplan-Yorke dimensions (see Table 1).Therefore, the high-dimensional phase space of the proposed hyperjerk system combined with its high Kaplan-Yorke dimension, especially in regard to the other reported systems, guarantees the system's complex dynamical behavior, which contributes to the required robustness.

Reported Work
Maximal Lypunov Exponent Kaplan-Yorke Dimension [53] 0.1320 3.1300 [57] 0.1555 3.1171 [59] 0.1448 3.1573 [60] 0.1422 3.1348 [61] 0.0730 3.1300 [62] 0.1250 3.1325 [63] 0.1320 3.1300 This work 0.1809 3.1405 In Figure 1, the hyperchaotic phase portraits of the proposed hyperjerk system (8) are illustrated.In more details, six different phase portraits, which are produced from the combinations of the four variables x 1 -x 4 , show the system's strange attractor for the selected set of parameters and initial conditions.Furthermore, in Figure 2, the Poincaré maps, which are produced by selecting two different planes, display the system (8)'s strange attractor.The study of continuous dynamical systems, like the proposed system (8), through a Poincaré map is one of the most popular topics in nonlinear dynamical analysis.This is done by taking intersections of the system's orbit into the plane x 2 = 0 with dx 3 /dt > 0. Naturally, for a n dimensional attractor, the Poincaré map gives rise to (n − 1) points, which can describe the dynamics of the attractor properly.Thus, this is an artificial way of reducing the map by dropping its dimension.In this direction, the organized set of points in the Poincaré maps of Figure 2 is an indication of system's chaotic behavior, while the discrete number of n points is an indication of periodic state of period-n.Finally, a closed curve in the Poincaré map is an indication of the system's quasiperiodic behavior.

Equilibrium Point Analysis
The 4D hyperchaotic hyperjerk system (8) has only one equilibrium point E(0, 0, 0, 0) and the Jacobian matrix, by linearizing the system at the equilibrium point E, is expressed as: Thus, system (8), at the equilibrium point E, has a characteristic equation, which can described as: By using system (8)'s characteristic Equation (10), we can find the stability of the equilibrium point E. It is obvious that the stability of system's equilibrium point E(0, 0, 0, 0) depends only on the parameters a, b.Therefore, by choosing a = 3.8 and b = 0.1, we find the eigenvalues as: λ 1 = 0.16385 + 1.89874i, λ 2 = 0.16385 − 1.89874i, λ 3 = −0.16385+ 0.498477i and λ 4 = −0.16385− 0.498477i.Based on the theory, according to the aforementioned eigenvalues, equilibrium point E is unstable.

Dissipativity and Invariance
It is also known from the theory that by using the general condition of dissipativity into the proposed system (8), we can find that Therefore, system (8) is obviously dissipative for c > 0. In addition, system (8)'s invariance is indicated by the coordinate transformation (x 1 , x 2 , x 3 , x 4 ) → (−x 1 , −x 2 , −x 3 , x 4 ).Therefore, (−x 1 , −x 2 , −x 3 , x 4 ) is also a solution for the same values of parameters a, b, c and the system can display symmetric attractors.

Analysis of the 4D Hyperjerk Dynamics
Next, the investigation of the effect of parameter a on system's dynamics is presented.Due to the fact that the values of parameters b, c do not affect the system's dynamics but only the size of the attractors, the study of the effect of the parameters b, c has been ignored in this work.The bifurcation diagram and the spectrum of the system's three largest Lyapunov exponents, by varying the value of the bifurcation parameter a in the range from 3.8 to 5, are presented in Figures 3 and 4, respectively.
The bifurcation diagram is a very important tool in nonlinear theory, as it provides an overview of the system's behavior in a range of values of a critical parameter, which is called a bifurcation parameter.These diagrams are produced with repetitive depictions by altering the bifurcation parameter with a small step.Therefore, similar to the Poincaré map, the discrete number of n points in the bifurcation diagram is an indication of periodic state of period-n, while a region with an infinite number of points is an indication of the system's chaotic behavior.For the production of the bifurcation diagram of Figure 3, we have chosen as the trajectory's cutting plane the x 2 = 0 with dx 3 /dt > 0, while the parameter a increases with very small step.In addition, the initial conditions at each iteration are (x 1 (0), x 2 (0), x 3 (0), x 4 (0)) = (0.1, 0.1, 0.1, 0.1).
A route to chaos through a quasiperiodic behavior, when decreasing the value of a, is observed in the bifurcation diagram of Figure 3.In more detail, the proposed system (8) displays periodic behavior for a ≥ 4.21, (period-1 steady state and period-2 steady state), which is confirmed from the phase portraits of Figure 5a,c and Poincaré maps of Figure 5b,d, respectively.Especially in the Poincaré maps, which are produced with the same technique as described for the production of the bifurcation diagram, a discrete number of points showing the period of system (8) are depicted.As the parameter a further decreases, the system through a quasiperiodic behavior, which appears in a very narrow band (4.207 ≤ a < 4.210), is driven to a complex behavior for a < 4.207 (Figure 5e).System (8)'s quasiperiodic behavior is confirmed from the closed curves in the Poincaré map of Figure 5f, according to the theory.In addition, the system's complex behavior for a < 4.207 is mainly hyperchaotic, except for a narrow band in which the system's behavior is chaotic.The diagram of the spectrum of the three largest Lyapunov exponents of Figure 4 confirms the system's hyperchaotic behavior due to the fact that, in the aforementioned region (a < 4.207), the system has two positive Lyapunov exponents.
Computation 2018, 6, x FOR PEER REVIEW 7 of 18 very narrow band (4.207 ≤ a < 4.210), is driven to a complex behavior for a < 4.207 (Figure 5e).System (8)'s quasiperiodic behavior is confirmed from the closed curves in the Poincaré map of Figure 5f, according to the theory.In addition, the system's complex behavior for a < 4.207 is mainly hyperchaotic, except for a narrow band in which the system's behavior is chaotic.The diagram of the spectrum of the three largest Lyapunov exponents of Figure 4 confirms the system's hyperchaotic behavior due to the fact that, in the aforementioned region (a < 4.207), the system has two positive Lyapunov exponents.Computation 2018, 6, x FOR PEER REVIEW 7 of 18 very narrow band (4.207 ≤ a < 4.210), is driven to a complex behavior for a < 4.207 (Figure 5e).System (8)'s quasiperiodic behavior is confirmed from the closed curves in the Poincaré map of Figure 5f, according to the theory.In addition, the system's complex behavior for a < 4.207 is mainly hyperchaotic, except for a narrow band in which the system's behavior is chaotic.The diagram of the spectrum of the three largest Lyapunov exponents of Figure 4 confirms the system's hyperchaotic behavior due to the fact that, in the aforementioned region (a < 4.207), the system has two positive Lyapunov exponents.

Circuit Realization of the Proposed System
The physical realization of theoretical chaotic models through the design of electronic circuits is the classical approach for the verification of their feasibility.Furthermore, the aforementioned approach has been applied in a large number of chaotic systems' engineering applications, such as true and pseudo random bit generation [63][64][65][66][67], chaotic video communication scheme via a wide area network (WAN) remote transmission [68], audio encryption scheme [69], autonomous mobile robots [70], chaotic communication systems [71], image encryption [72][73][74][75], etc.For this reason, analog and digital approaches have been applied to realize chaotic oscillators by using different kinds of electronic devices such as common off-the-shelf electronic components [76,77], integrated circuit technology [78,79], microcontroller [80] or field-programmable gate array (FPGA) [81][82][83].
Therefore, in this section, the circuit implementation of the proposed hyperjerk system (8) is presented.In more detail, Figure 6 depicts the circuit that has been designed by using twenty-five resistors, four capacitors, nine operational amplifiers and three analog multipliers (AD633).Four of the operational amplifiers (U1-U4) are configured as Resistor-Capacitor integrators, which is an important operation required in the derivation of the equivalent electronic circuit.The specific realization of the RC integrators, which are very well known in literature, exploits the properties of the operational amplifiers in the linear region, producing in general a current , where VC is the voltage accross the capacitor.Other three operational amplifiers have been used as inverting amplifiers (U5-U7), which use negative feedback to amplify the input voltage.The last ones (U8, U9) with the two diodes (1N4007) are used for implementing the absolute nonlinearity ( 2x ).The selected configuration for the absolute nonlinearity gives the best results according to literature.Furthermore, the three multipliers (U10-U12) are used for implementing the quintic term c 4 1 4 x x .The low-cost analog multipliers AD633, which are used for this purpose, is the most common approach

Circuit Realization of the Proposed System
The physical realization of theoretical chaotic models through the design of electronic circuits is the classical approach for the verification of their feasibility.Furthermore, the aforementioned approach has been applied in a large number of chaotic systems' engineering applications, such as true and pseudo random bit generation [63][64][65][66][67], chaotic video communication scheme via a wide area network (WAN) remote transmission [68], audio encryption scheme [69], autonomous mobile robots [70], chaotic communication systems [71], image encryption [72][73][74][75], etc.For this reason, analog and digital approaches have been applied to realize chaotic oscillators by using different kinds of electronic devices such as common off-the-shelf electronic components [76,77], integrated circuit technology [78,79], microcontroller [80] or field-programmable gate array (FPGA) [81][82][83].
Therefore, in this section, the circuit implementation of the proposed hyperjerk system (8) is presented.In more detail, Figure 6 depicts the circuit that has been designed by using twenty-five resistors, four capacitors, nine operational amplifiers and three analog multipliers (AD633).Four of the operational amplifiers (U1-U4) are configured as Resistor-Capacitor integrators, which is an important operation required in the derivation of the equivalent electronic circuit.The specific realization of the RC integrators, which are very well known in literature, exploits the properties of the operational amplifiers in the linear region, producing in general a current i = C dV C dt , where V C is the voltage accross the capacitor.Other three operational amplifiers have been used as inverting amplifiers (U5-U7), which use negative feedback to amplify the input voltage.The last ones (U8, U9) with the two diodes (1N4007) are used for implementing the absolute nonlinearity (|x 2 |).The selected configuration for the absolute nonlinearity gives the best results according to literature.Furthermore, the three multipliers (U10-U12) are used for implementing the quintic term c x 4  1 x 4 .The low-cost analog multipliers AD633, which are used for this purpose, is the most common approach for realizing polynomial terms or products, like the quintic term cx 4 x 1 .With the use of Kirchhoff's circuit laws into the proposed circuit of Figure 6, we derive the mathematical model of the hyperjerk system, which is described by the following equations: With the use of Kirchhoff's circuit laws into the proposed circuit of Figure 6, we derive the mathematical model of the hyperjerk system, which is described by the following equations: where x 1 , x 2 , x 3 and x 4 correspond to the voltages on the integrators (U1-U4), respectively, while the power supply is ±15 V. System ( 12) is normalized by using τ = t/RC.It can thus be suggested that system ( 12) is equivalent to system (8 The designed circuit of Figure 6 has been implemented by using the well-known Multisim program.In this way, PSpice results of the phase portraits in various planes are reported in Figure 7.It is easy to see the good agreement between system (8)'s numerical results (Figure 1) and the circuit's simulation results (Figure 7).
where x1, x2, x3 and x4 correspond to the voltages on the integrators (U1-U4), respectively, while the power supply is ±15 V. System ( 12) is normalized by using τ = t/RC.It can thus be suggested that system ( 12) is equivalent to system (8 The designed circuit of Figure 6 has been implemented by using the well-known Multisim program.In this way, PSpice results of the phase portraits in various planes are reported in Figure 7.It is easy to see the good agreement between system (8)'s numerical results (Figure 1) and the circuit's simulation results (Figure 7).(e) (f)

Synchronization Scheme
As it has been mentioned, the synchronization of chaotic systems is an interesting research topic in literature due to its applications in secure communication schemes and cryptosystems.In order to study the possibility of synchronization of two coupled identical 4D hyperchaotic hyperjerk systems with an unknown system's parameters, an adaptive control law has been used.The adaptive control is the main method, which is used when some or all the system parameters are not available for measurement and estimates for the uncertain parameters of the systems [16,52].Therefore, it could be useful in the aforementioned applications.
Therefore, in this work, as the master system, the proposed 4D hypechaotic hyperjerk system (8) is used, while the slave system is given by the following system dynamics: where yi, (i = 1, ..., 4) are the states and ui, (i = 1, ..., 4) are the adaptive controls to be determined.As it is mentioned, the parameters a, b, c, in systems ( 8) and ( 13) are unknown and the design goal is to find adaptive feedback controls ui that uses estimates for the parameters a, b, c, respectively, so as to render the states of the systems (8) and ( 13) fully synchronized asymptotically.The synchronization error between the chaotic systems (8) and ( 13) is defined as: We take the adaptive control laws defined by

Synchronization Scheme
As it has been mentioned, the synchronization of chaotic systems is an interesting research topic in literature due to its applications in secure communication schemes and cryptosystems.In order to study the possibility of synchronization of two coupled identical 4D hyperchaotic hyperjerk systems with an unknown system's parameters, an adaptive control law has been used.The adaptive control is the main method, which is used when some or all the system parameters are not available for measurement and estimates for the uncertain parameters of the systems [16,52].Therefore, it could be useful in the aforementioned applications.
Therefore, in this work, as the master system, the proposed 4D hypechaotic hyperjerk system (8) is used, while the slave system is given by the following system dynamics: where y i , (i = 1, ..., 4) are the states and u i , (i = 1, ..., 4) are the adaptive controls to be determined.As it is mentioned, the parameters a, b, c, in systems ( 8) and ( 13) are unknown and the design goal is to find adaptive feedback controls u i that uses estimates for the parameters a, b, c, respectively, so as to render the states of the systems (8) and ( 13) fully synchronized asymptotically.The synchronization error between the chaotic systems (8) and ( 13) is defined as: Thus, the synchronization error dynamics is given by the following equations: We take the adaptive control laws defined by where k i , (i = 1, ..., 4) are positive gain constants and â, b, ĉ are the parameter update laws.Substituting Label (16) into Label (15), we obtain the closed-loop error dynamics as: As parameters' estimation errors, we used the equations: By differentiating Equations (18) with respect to t, we can obtain . . .
Therefore, by using Equations (18), we can rewrite the closed-loop system (18) as: Next, the following quadratic Lyapunov function is used: V(e, e a , e b , e c ) = By differentiating the quadratic Lyapunov function V of Equation ( 22) along the trajectories of the systems (21) and (20), we obtain the following: Therefore, the following parameter update laws we have taken are: Next, we establish the main results of this section.
Theorem 1.The 4D hyperchaotic hyperjerk systems (8) and ( 13) with unknown parameters are globally and exponentially synchronized for all initial conditions by the adaptive feedback control law ( 16) and the parameter update laws (23), where k i , (i = 1, ..., 4) are positive constants.
Proof.We prove the aforementioned theorem by using the Lyapunov stability theory.For this reason, the quadratic Lyapunov function V of Equation ( 21), which is positive definite on R 7 , is used.Therefore, the time derivative of V, by substituting the parameter update laws ( 23) into ( 22) is obtained as: .
From Equation ( 24), it is clear that .
V is a negative semi-definite function on R 7 .In addition, it is concluded that the synchronization error vector e(t) = (e 1 (t), e 2 (t), e 3 (t), e 4 (t)) and the parameter estimation error (e a (t), e b (t), e c (t)) are globally bounded.
By defining k = min(k 1 , k 2 , k 3 , k 4 ), then it follows from Equation ( Thus, k e(t) By integrating the inequality (26) from 0 to t, we have taken: Therefore, from Equation ( 27), it follows that e ∈ L 2 , while from system (16), it can be deduced that .e ∈ L ∞ .Thus, we can conclude that for all initial conditions e(t) → 0 exponentially as t → ∞, by using the Barbalat's lemma [84].
This completes the proof.As in Section 3, the classical fourth-order Runge-Kutta method is used in order to solve the system of differential Equations (8), ( 13) and ( 23), when the adaptive control laws (16) are applied.In addition, the same parameter values of the 4D hyperjerk systems (8) and ( 13) as in the hyperchaotic case of Section 2 are used, while the gain constants are taken as k i = 10, for i = 1, 2, 3, 4.

USV Symbol Macro(s) Description
Furthermore, the set of initial conditions of the master system (8), is (x 1 (0), x 2 (0), x 3 (0), x 4 (0)) = (0.1, 0.1, 0.1, 0.1), while the set of initial conditions of the slave system (13), is (y 1 (0), y 2 (0), y 3 (0), x 4 (0)) = (-0.1, 0.2, -0.2, 0).In addition, we use the values of â(0) = 3.5, b(0) = 0.15, ĉ(0) = 1.2, as initial conditions of the parameter estimates.The synchronization of the states of the master system (8) and slave system (13) are depicted in Figure 8, while the time-history of the synchronization errors e 1 (t), e 2 (t), e 3 (t), e 4 (t) is depicted in Figure 9.In more details, Figure 8 shows the fast convergence of the respective signals of the two coupled systems by using the proposed synchronization scheme.Furthermore, Figure 9 confirms the feasibility of the synchronization method in the coupled system by confirming that the synchronization errors converge to zero, as it is expected.Therefore, a chaotic complete synchronization has been achieved.
(0.1, 0.1, 0.1, 0.1), while the set of initial conditions of the slave system (13), is (y1(0), y2(0), y3(0), x4(0)) = (-0.1,0.2, -0.2, 0).In addition, we use the values of , as initial conditions of the parameter estimates.The synchronization of the states of the master system (8) and slave system (13) are depicted in Figure 8, while the time-history of the synchronization errors e1(t), e2(t), e3(t), e4(t) is depicted in Figure 9.In more details, Figure 8 shows the fast convergence of the respective signals of the two coupled systems by using the proposed synchronization scheme.Furthermore, Figure 9 confirms the feasibility of the synchronization method in the coupled system by confirming that the synchronization errors converge to zero, as it is expected.Therefore, a chaotic complete synchronization has been achieved.

Discussion
This paper focuses on the study of a 4D hyperchaotic hyperjerk dynamical system (8) with only two nonlinear terms.Recently, Elhadj and Sprott indicated that the hyperjerk form (3) can describe all periodically forced oscillators and many of the coupled oscillators [11].Therefore, the proposed system (8) could be a suitable candidate for use in order to model such kind of oscillators.As proved, the maximal Lyapunov exponent (MLE) and the Kaplan-Yorke dimension of the proposed hyperchaotic system (8) is greater than the respective features of the other hyperjerk systems of the same category, which has been presented in literature.Therefore, the high-dimensional phase space of the hyperjerk system (8) combined with its high Kaplan-Yorke dimension, which guarantees the system's complex dynamical behavior, contributes to the required robustness.In addition, the proposed system (8) is a dissipative system for c > 0 and it has an unstable equilibrium point at E(0, 0, 0, 0).Furthermore, it is invariant under the coordinate transformation (x1, x2, x3, x4) → (−x1, −x2, −x3, x4), which is a very important feature that could drive the system to coexisting attractors.

Discussion
This paper focuses on the study of a 4D hyperchaotic hyperjerk dynamical system (8) with only two nonlinear terms.Recently, Elhadj and Sprott indicated that the hyperjerk form (3) can describe all periodically forced oscillators and many of the coupled oscillators [11].Therefore, the proposed system (8) could be a suitable candidate for use in order to model such kind of oscillators.As proved, the maximal Lyapunov exponent (MLE) and the Kaplan-Yorke dimension of the proposed hyperchaotic system (8) is greater than the respective features of the other hyperjerk systems of the same category, which has been presented in literature.Therefore, the high-dimensional phase space of the hyperjerk system (8) combined with its high Kaplan-Yorke dimension, which guarantees the system's complex dynamical behavior, contributes to the required robustness.In addition, the proposed system (8) is a dissipative system for c > 0 and it has an unstable equilibrium point at E(0, 0, 0, 0).Furthermore, it is invariant under the coordinate transformation (x1, x2, x3, x4) → (−x1, −x2, −x3, x4), which is a very important feature that could drive the system to coexisting attractors.

Discussion
This paper focuses on the study of a 4D hyperchaotic hyperjerk dynamical system (8) with only two nonlinear terms.Recently, Elhadj and Sprott indicated that the hyperjerk form (3) can describe all periodically forced oscillators and many of the coupled oscillators [11].Therefore, the proposed system (8) could be a suitable candidate for use in order to model such kind of oscillators.As proved, the maximal Lyapunov exponent (MLE) and the Kaplan-Yorke dimension of the proposed hyperchaotic system (8) is greater than the respective features of the other hyperjerk systems of the same category, which has been presented in literature.Therefore, the high-dimensional phase space of the hyperjerk system (8) combined with its high Kaplan-Yorke dimension, which guarantees the system's complex dynamical behavior, contributes to the required robustness.In addition, the proposed system (8) is a dissipative system for c > 0 and it has an unstable equilibrium point at E(0, 0, 0, 0).Furthermore, it is invariant under the coordinate transformation (x 1 , x 2 , x 3 , x 4 ) → (−x 1 , −x 2 , −x 3 , x 4 ), which is a very important feature that could drive the system to coexisting attractors.
The dynamical behavior of system (8) has been studied using well-known tools of nonlinear theory, such as bifurcation diagrams, Lyapunov exponents, phase portraits and Poincaré maps.The effect of the parameter a is investigated and the route to chaos through a quasiperiodic behavior, in a narrow band, is driven to hyperchaos, which is maintained in a large range of values of the parameter a.The quasiperiodic behavior has been confirmed by the closed loop, which is produced in system's Poincaré map, for a selected value of the parameter a (a = 4.209) from quasiperiodic's narrow band region.Therefore, despite its simplicity, the system has an enlarged area of hyperchaotic behavior, which is not so usual in other systems of this family.
In order to explore the feasibility of the proposed hyperchaotic system (8), the classical approach, by using the physical realization through the design of an electronic circuit, has been adopted.For this reason, common off-the-shelf electronic components, such as op-amps and multipliers, have been used in this work.PSpice results confirm the good agreement between the circuit's simulation results and numerical results.In more details, the circuit's phase portraits that have been produced by using the PSpice are identical in form and size with the respective system's simulated ones.Therefore, the proposed system could be used as a hyperchaotic generator in many practical applications, such as in encryption and secure communication schemes.
Furthermore, an adaptive synchronization scheme for a system of two coupled 4D hyperchaotic hyperjerk systems has been developed.For the obtained, with this method, the main adaptive results are proved by using Lyapunov stability theory.The effectiveness of the proposed synchronization scheme is proved from the simulation results, in which the fast convergence of the signals of the two coupled systems has been achieved.
One of the possible future research directions is the investigation of the coexistence of multiple attractors, based on system's invariance.In addition, the experimental realization of the circuit that is designed and presented in this work will be done in the near future.Finally, an engineering application, such as an encryption or secure communication, based on the proposed 4D hyperchaotic hyperjerk dynamical system may be designed.

Figure 4 .Figure 3 .
Figure 4. Spectrum of three largest Lyapunov exponents of system (8) versus the parameter a, for b = 0.1 and c = 1.5.

Figure 4 .
Figure 4. Spectrum of three largest Lyapunov exponents of system (8) versus the parameter a, for b = 0.1 and c = 1.5.

Figure 4 .
Figure 4. Spectrum of three largest Lyapunov exponents of system (8) versus the parameter a, for b = 0.1 and c = 1.5.

Figure 6 .
Figure 6.Schematic of the designed circuit for the proposed system (8).

Figure 6 .
Figure 6.Schematic of the designed circuit for the proposed system (8).
with a = R/R a , b = R/R b and c = R/R C .The values of circuit components are R = 100 kΩ, R 1 = 10 kΩ, R 2 = 90 kΩ, R b = 1 MΩ, R c = 66.666 kΩ, and C = 1 nF.In order to change the parameter a, a variable resistor R a can be used.
synchronization error dynamics is given by the following equations:

Table 1 .
Comparison of system (8) with other reported hyperchaotic hyperjerk systems regarding the maximal Lyapunov exponent and the Kaplan-Yorke dimension.