Abstract

We announce a new 4D hyperchaotic system with four parameters. The dynamic properties of the proposed hyperchaotic system are studied in detail; the Lyapunov exponents, Kaplan-Yorke dimension, bifurcation, and bicoherence contours of the novel hyperchaotic system are derived. Furthermore, control algorithms are designed for the complete synchronization of the identical hyperchaotic systems with unknown parameters using sliding mode controllers and genetically optimized PID controllers. The stabilities of the controllers and parameter update laws are proved using Lyapunov stability theory. Use of the optimized PID controllers ensures less time of convergence and fast synchronization speed. Finally the proposed novel hyperchaotic system is realized in FPGA.

1. Introduction

The study of chaotic attractors arising in nonlinear dynamical systems has received much attention as this has potential applications in many branches of science and engineering. Chaotic systems are characterized as nonlinear dynamical systems which are irregular, aperiodic, and unpredictable and have sensitive dependence on the initial conditions [1, 2]. Lorenz discovered a chaotic system while he was studying 3D weather model [3]. Some important paradigms of 3D autonomous chaotic systems are Rössler system [4], Chen system [5], Lü system [6], and so forth. Some novel 3D autonomous chaotic systems are Zhu system [7], Vaidyanathan systems [810], Pehlivan system [11], Akgul system [12], Tacha system [13], and so forth.

A hyperchaotic system is mathematically defined as a chaotic system having more than one positive Lyapunov exponent implying that its dynamics are expended in many different directions simultaneously. Thus, a hyperchaotic system has more complex dynamical behaviours than a chaotic system [1, 2]. Rössler was the first to discover a hyperchaotic system [14]. There are many well-known hyperchaotic systems such as hyperchaotic Lorentz system [15], hyperchaotic Chen system [16], hyperchaotic Li system [17], hyperchaotic Vaidyanathan systems [1820], and so forth.

Sliding mode control is one of the best robust control algorithms as discussed in literature [2123]. Most of the works were concentrated towards active sliding mode control assuming that the system parameters are known. Some authors [24, 25] have derived a methodology to synchronize systems assuming the parameters are unknown. Systems with uncertainty and external disturbance is proposed to realize chaos control [25] in which the bounds of the uncertainty and external disturbance are assumed to be unknown. FPGA implementation of the chaotic generators is more useful for real-time applications such as secure communications and cryptography [26]. A new methodology for implementing integer order integrators in FPGA using Xilinx System Generator toolbox was proposed [27].

In this paper, we have derived a novel hyperchaotic system and described the qualitative properties of the system. The Lyapunov exponents of the proposed system are having two positive exponents showing that the system is hyperchaotic. The bifurcation and bicoherence contours of the hyperchaotic system are studied to investigate the parameter dependence and quadratic nonlinearities, respectively.

The synchronization of hyperchaotic systems has many applications like secure communications [28, 29], cryptosystems [30, 31], and so forth. In this work, we also derive new results for the adaptive synchronization of the novel hyperchaotic systems with unknown parameters using sliding mode control and PID control [32] optimized by genetic algorithm [33, 34]. The main adaptive control results are established via Lyapunov stability theory [35]. Finally the novel hyperchaotic system is realized in FPGA using the Xilinx (Vivado) System Generator toolbox.

2. A Novel Hyperchaotic System

In this work, we propose a novel hyperchaotic system described bywhere are the state variables and are positive constant parameters.

In this work, we show that system (1) is hyperchaotic when we take the parameters as

For numerical simulations, we take the initial values of the hyperchaotic system (1) as

Figures 14 show the 3D projections of the novel hyperchaotic system (1) on , and spaces, respectively.

3. Properties of the Novel Hyperchaotic System

In this section, we discuss the qualitative properties of the novel hyperchaotic system (1) introduced in Section 2. We suppose that the parameter values of system (1) are as in the hyperchaotic case for the system parameters (2).

3.1. Dissipativity

In vector notation, the 4D system (1) can be expressed aswhere

Let be any region in with a smooth boundary and also, , where is the flow of the vector field Furthermore, let denote the hypervolume of

By Liouville’s theorem, we have

The divergence of the vector field is easily calculated aswhere

Substituting (7) into (6), we obtain the first-order differential equation

Integrating (8), we obtain the unique solution as

Since , it follows that exponentially as This shows that the novel hyperchaotic system (1) is dissipative. Thus, the system limit sets are ultimately confined into a specific limit set of zero volume, and the asymptotic motion of the novel hyperchaotic system (1) settles onto a strange attractor of the system.

3.2. Equilibrium Points

The equilibrium points of the novel hyperchaotic system (1) are obtained by solving

The parameter values are taken as in the hyperchaotic case (2).

Then the equilibrium points of the novel hyperchaotic system (1) are found as

The Jacobian matrix of the novel hyperchaotic system (1) at any point is found as

The eigenvalues of are easily calculated as

This shows that is a saddle point, which is unstable.

The eigenvalues of are easily calculated as

This shows that is a saddle-focus, which is unstable.

The eigenvalues of are easily calculated as

This shows that is a saddle-focus, which is unstable.

Hence, all the three equilibrium points of the novel hyperchaotic system (1) are unstable.

3.3. Lyapunov Exponents and Kaplan-Yorke Dimension

We take the parameters of the novel 4D system (1) as in the hyperchaotic case (2).

Then the Lyapunov exponents of the novel 4D system (1) are numerically found as

Since there are two positive Lyapunov exponents in (16), it is clear that the 4D system (1) is hyperchaotic.

We note that the sum of the Lyapunov exponents of system (1) is negative. In fact,

This shows that the novel hyperchaotic system (1) is dissipative.

Also, the Kaplan-Yorke dimension of the hyperchaotic system (1) is derived aswhich is fractional.

3.4. Bicoherence

The bicoherence or the normalized bispectrum is a measure of the amount of phase coupling that occurs in a signal or between two signals. Both bicoherence and bispectrum are used to find the influence of a nonlinear system on the joint probability distribution of the system input. Phase coupling is the estimate of the proportion of energy in every possible pair of frequency components . Bicoherence analysis is able to detect coherent signals in extremely noisy data, provided that the coherency remains constant for sufficiently long times, since the noise contribution falls off rapidly with increasing number of segments

The autobispectrum of a chaotic system is given by Pezeshki et al. [36]. They derived the autobispectrum with the Fourier coefficientswhere 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 .

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 statistics, formed from , where is the Fourier transform of , the bispectrum is a third-order statistics 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 M segments of length N_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.

Figures 5 and 6 depict the contours of the bicoherence of the states and for the novel hyperchaotic system (1). Figure 7 shows the bicoherence contours of all the states together. Shades in yellow represent the multifrequency components contributing to the power spectrum. From Figures 5, 6, and 7, the cross-bicoherence is significantly nonzero and nonconstant, indicating a nonlinear relationship between the states. In Figures 5 and 6 (autobicoherence) the yellow shades indicate that the nonlinear relationships of the states and are not narrow-band processes. The yellow shades and nonsharpness of the peaks as well as the presence of structure around the origin in Figure 7 (cross-bicoherence) indicate that the nonlinearity between the states 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.

For the bicoherence plots we used the pseudocode [37] given below.

Analysis of the Input Dataif () thenNot enough data’Is the segment length supplied?if then Is the sampling period supplied?if then

Break the Time Series into M Segments of N Length Eachif () thenFrequency resolutionFrequency array

Calculate the Periodogram for Each Segment Using FFT Command Calculate mean value (DC component) for each segmentSubtract mean (DC component)Calculate the Fourier transform of each segment (row)this gives data = fft(data, 1, dimension = 1,/overwrite)Extract only positive frequencies Calculate the index matrix

Prepare Arrays for Storing Components of Bicoherence where is the complex value, and are the real values

Normalize So That Is in UnitsCalculate the matrix Calculate the matrix

Calculate the Triple Term

Sum over the Ensembles of Each Segment

Define Power Spectrum and its Poisson noise level

Normalize the Sums to Get Means

Calculate the Bicoherence from

Masking the Variables to Avoid/0 Errorsmask = where (, count)Prepare array for 2-dimensional bicoherencenormalize to get bicoherencePlot the results.

3.5. Bifurcation

In order to understand the dynamical behaviour of the MNECS, the bifurcation plots are derived for three cases as follows.

Case 1. Fix and vary between .

Case 2. Fix and vary between .

Case 3. Fix and vary between .

Figures 8(a), 8(b), and 8(c) show the bifurcation plots of the system for the parameters , respectively. As it is observed from Figures 8(a), 8(b), and 8(c), for the values of , the bifurcation plots show denser points confirming the existence of two stable positive Lyapunov exponents (hyperchaos) showing the parameter dependence of system (1).

4. Adaptive Synchronization of the Novel Hyperchaotic Systems Using Sliding Mode Control

In this section, we deploy the adaptive sliding mode control (ASMC) method to derive new results for the complete synchronization of the novel hyperchaotic systems using sliding mode control.

The master system is given as The slave systems with the adaptive sliding mode controllers () are defined asThe parameters of the slave system are assumed to be unknown and hence replaced by their respective parameter estimates .

For the complete synchronization of the master and slave systems, the synchronization errors , where and is defined as The error dynamics are derived by differentiating the synchronization errorsUsing the master system (21) and slave system (22) in the error dynamics (24),

4.1. Control Design

The main control objective of this paper is to design adaptive sliding mode controllers () such that the synchronization errors . Sliding mode control methodology is a simple approach to robust control and good at dealing with dynamic uncertainty [25, 38]. The control design procedure consists of two steps: first constructing a sliding surface which presents the desired dynamics [39] and second selecting a switching control law so as to verify sliding condition.

The integral sliding surface [38] can be defined aswhere and .

The time derivate of the sliding surface is given byFor any initial conditions, the problem of synchronization of master and slave systems is equivalent to that of the synchronization errors remaining on the surface for all . When the system operates in the sliding surface, it satisfies and .

The parameter estimation errors are given byLet the adaptive sliding mode controllers be defined aswhere are the sliding surface gains and is the controller gain for

Let us define the parameter adaptive laws as where are the dynamics of the parameter estimates .

4.2. Stability Analysis of the Controller

Theorem 1. Considering that adaptive sliding mode control input law in (5) is used to control error system in (6) with update laws of parameters in (7), then error system in (6) is globally asymptotically stable.

Proof. To check the stability of the controlled system, let us consider the following Lyapunov candidate function:The derivative of the Lyapunov function is given by Using (27) and (25) in (32),After some mathematical manipulations we arrive atwhere is the sliding reaching law [38], , and are the sliding surface gains.
Introducing the parameter update laws (30) and the adaptive sliding mode controllers (29) in (34), we simplify the Lyapunov first derivative as As can be seen from (35), Lyapunov first derivative is a negative semidefinite as and , hence, by Barbalat’s lemma in Lyapunov stability theory [36, 40], it follows that the synchronization errors converge exponentially to zero as
For numerical simulations, we take the initial conditions of the master system (27) asWe take the initial conditions of the slave system (28) asWe take the initial conditions of the parameter estimates asFigure 9 depicts the time-history of the synchronization errors and Figure 10 shows the synchronized states. Figure 11(a) shows the time-history of the states. Figure 11(b) shows the time-history of adaptive sliding mode (ASM) controllers.

5. PID Controller Based on Genetic Algorithms (GAPID)

Many advanced and intelligent control algorithms have been developed, however, the PID-type controller [16] remains the most popular in industry for the simplicity of control law, few tuning parameters, less time convergence and fast synchronization speed. But finding appropriate parameters values for the PID controller is still a difficult task, so in practice control engineers still often use trial and error for the tuning process. So in practice control engineers still often use trial and error for the tuning process.

PID control consists of three types of control, proportional, integral, and derivative control; they can be implemented in many forms but they are mostly used in feedback loops.

In this paper we use four PID controllers () as given in (39) to synchronize the master system states () with slave system states ().where is the PID controller for is the error signal, and are the proportional, integral, and derivative gains of the PID controller for the synchronization of the master novel hyperchaos system (21) and the slave system (22). Tuning of PID controllers involves the selection of gains for better control performance which is defined with reference to the required performance index.

In this work we will investigate the effectiveness of the PID controller gains when optimized using genetic algorithms (GA). In particular, many PID control based on GA schemes have been successfully applied to the control and the synchronization of chaotic systems [41], GA’s [33, 34] are a stochastic global search method that mimics the process of natural evolution. It starts with no knowledge of the correct solution and depends entirely on responses from its environment and evolution operators (i.e., reproduction, crossover and mutation) to arrive at the best solution. By starting at several independent points (population) and searching in parallel, the algorithm avoids local minima and converging to suboptimal solutions. A genetic algorithm is typically initialized with a random population consisting of between 20–100 individuals. This population is usually represented by a real-valued number or a binary string called a chromosome. How well an individual performs a task is measured and assessed by the objective function. The objective function assigns each individual a corresponding number called its fitness. The fitness of each chromosome is assessed and a survival of the fittest strategy is applied. The magnitude of the error is usually used to assess the fitness of each chromosome. Writing an objective function is the most difficult part of creating a genetic algorithm. An objective function could be created to find a PID controller that gives the smallest overshoot, fastest rise time or quickest settling time, however, in order to combine all of these objectives it was decided to design an objective function that will minimize the error of the controlled system. Each chromosome in the population is passed into the objective function one at a time. The chromosome is then evaluated and assigned a number to represent its fitness, the bigger its number the better its fitness. The genetic algorithm uses the chromosome’s fitness value to create a new population consisting of the fittest members. There are two important performance indices ISE (Integrated Squared Error) and IAE (Integrated Absolute Error) given as

In this paper we use IAE as the objective function and the fitness functions is given as

Using the fitness function (41) and objective function (40) in the optimization toolbox of Matlab with the following options, we derive the optimal values for each of the controllers ().Variable bounds matrix = .Population size = 80; GA. Generally the bigger the population size, the better the final approximation.Number of generations = .Selection function = stochastic uniform.Crossover fraction = .Mutation function = Gaussian.Stopping criteria = error performance criterion.Length of the chromosome = , decimal codage.

With the same parameter values used in (2) and the same the initial conditions in (36) and (37), Table 1 shows the PID gain values after running the genetic algorithm (GA) solver in the optimization tool with the options cited above, we get the best solutions tracked over generations for the complete synchronization of the novel hyperchaotic system:

Figure 12 depicts the time-history of the synchronization errors using genetically optimized PID controllers, Figure 13 shows the synchronized states, and Figure 14 shows the time-history of the PIDs controllers.

As seen from Figures 8 and 11, in the ASMC method of synchronization the errors converges at whereas the GAPID controllers based synchronization the errors converges at . The states synchronization Figures 9 and 12 and the controllers convergence plots shown in Figures 11 and 14 clearly indicates that GAPID controllers are efficient compared to the ASMC controllers. As can be seen from Figures 11 and 14, ASM controllers converging time () is slower than the GAPID controllers (). Also the ASM controllers shoot to peak values of which when realized in electronic circuits need special attenuators to suppress the oscillations whereas in GAPID controllers the shoot values are reduced comparatively.

6. FPGA Implementation of the Novel Hyperchaotic System

In this section we discuss about the implementation of the novel hyperchaotic system in FPGA using the Xilinx (Vivado) System Generator toolbox in Simulink. Firstly we configure the available built in blocks of the System Generator toolbox. The Add/Subblocks are configured with zero latency and 32/16 bit fixed point settings. The output of the block is configured to rounded quantization in order to reduce the bit latency. For the multiplier block a latency of 3 is configured and the other settings are the same as in Add/Subblock. Next we will have to design the integer order integrator which is not a readily available block in the System Generator. Hence we implement the integrators using the mathematical relation and the value of is taken as and the initial conditions are fed in to the forward register. Figure 15 shows the Xilinx RTL schematics of the system using Virtex 7-XC7 chips and Figures 16, 17, 18, and 19 shows the 3D phase portraits of the novel hyperchaotic system using Xilinx System Generator. Comparing Figures 1, 2, 3, and 4 with Figures 16, 17, 18, and 19 one can clearly see that the FPGA implemented novel hyperchaotic system exhibits the same phase portraits for the initial conditions .

7. Conclusions

In this research work, we first described a novel hyperchaotic system and discussed its qualitative properties. Lyapunov exponents are derived to prove that the system is hyperchaotic. Bicoherence plots are drawn to prove the multifrequency terms dependency of the novel hyperchaotic systems. Next, identical hyperchaotic systems are synchronized using adaptive sliding mode control assuming that the system parameters are unknown and also genetically optimized PID controllers. The complete synchronization results derived in this paper have been proved using the Lyapunov stability theory. A comparison between the ASM controllers and GAPID controllers are discussed with various plots. Finally the 4D novel hyperchaotic system is implemented in FPGA implementation to show that the system is hardware realizable.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.