# A Novel Hypogenetic Chaotic Jerk System: Modeling, Circuit Implementation, and Its Application

Authors: Jiancheng Liu; Karthikeyan Rajagopal; Tengfei Lei; Sezgin Kaçar; Burak Arıcıoğlu; Ünal Çavuşoğlu; Abdullah Hulusi Kökçam; Anitha Karthikeyan
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

When revising the polarity and amplitude information in the feedback, a unique hypogenetic jerk system was obtained which has two controllers to switch the equilibria between stable and unstable. After providing some basic dynamical analysis, an electronic circuit was implemented, and the phase trajectory in the oscilloscope agrees with the numerical simulation. Further exploration shows that this unique chaotic system has superior performance as a random number generator or in voice encryption application.


---

## Body

## 1. Introduction

In the literature, chaos has a wide range of application field. In the recent years, chaos or chaotic systems have been frequently employed in the encryption and random number generation (RNG) studies due to noise-like, aperiodic characteristics of the chaotic systems [1]. In the literature, the encryption studies are not only about image encryption [2–5] but also there are encryption studies about text and other multimedia types like audio and video encryption [6]. It shows that chaos has great performance in the encryption system. In fact, RNG studies in the literature can be categorized into true random number generation (TRNG) and pseudo-random number generation (PRNG) regardless of whether they are chaos-based or not [7, 8]. RNG is usually employed in cryptographic studies for key generation process and performance of the encryption which heavily depends on key randomness. This can be achieved with appropriate chaos-based random number generation process. All these constitute our motivation to develop a new chaotic system.Jerk dynamical systems are realized with a compact electrical circuit structure. Generally, jerk circuits have four connections at the nodex, where the derivative of x is determined by the amplitude and polarity of y, and the amplitude and polarity of x influence the derivative of z. We study those chaotic flows with incomplete information transmission from the node x based on the jerk structure and are therefore named hypogenetic chaotic. Jerk system has the simple structure but can also provide applicable chaotic signal [9–13]. When the feedback information from the other variables is not complete, many of the systems can remain chaos [14–16] even in the jerk structure as hypogenetic flows [14]. Hot research on chaotic systems mainly focuses on chaos control [17–20], Lyapunov exponent calculation and analysis [21, 22], and doubling and growth of attractors [23, 24], and chaotic systems with hidden attractors are also hot research topics because such systems are extremely prone to multistable phenomena, which are a common phenomenon in nature. Considering from equilibrium point, hidden attractors can be mainly of several types, one stable equilibrium[25], a line or plane equilibrium [26, 27], or no equilibrium [28]. Even there are hidden attractors found in some special chaotic systems [29], which have both unstable equilibrium and stable equilibria.Starting from the classic system, hidden attractor is a very important topic. At the same time, linearization has always been a fundamental state. Is there a class of system whose linear linearization satisfies the characteristics of the chaotic hidden attractor (stable equilibrium system)? The circuit of this kind of system is less disturbed because linearization (sign circuit) is stable and similar to the digital circuit.Therefore, we propose a class of hypogenetic chaotic system which has the following unique properties: nonlinearity involved includes amplitude information and polarity information, which are of incomplete feedback; there are two knobs which control the stability of equilibria, by which one can turn off the stable states freely.This paper is organized as follows. In Section2, we give a model description including basic dynamical analysis, and in Section 3, we show the electronic circuit implementation of the system. Further application discussions including random number generator design and voice encryption application are given in Section 4. In the last section, we give a simple conclusion.


## 2. Model Description

In this letter, we announce a new jerk system derived by modifying the jerk system proposed in [14] by replacing the nonlinear terms z2 by z and xz by zsgnx. The proposed system shows a unique property of disorder which has both stable and unstable equilibrium points for different values of parameters.(1)x˙=a1y,y˙=a2z,z˙=a3x+a4z+a5xy+a6zsgnx+a7,where ai for i∈1,7 are the parameters of the system. Jerk systems have the form x˙=yy˙=zz˙=fx,y,z, which can be depicted as x⃛=fx,x˙,x¨ [30]. When a1 = 1 and a2 = 1, system (1) is Jerk system.It is simple to verify that the equilibrium point of system (1) is x=−a7/a3,y=0,andz=0, and the characteristic equation of the system is λ3+a6sgna7/a3λ2+a2a5a7/a3λ−a1a2a3=0. According to the Routh–Hurwitz criterion for the real part of the eigenvalues to be positive, the condition is δ1,δ0>0,δ1δ2>δ0, where δ0=1,δ1=a6sgna7/a3,δ2=a2a5a7/a3,andδ3=−a1a2a3. For the values of parameters a1=1,a2=1, a3=−1,a4=−4,anda6=1 and to discuss the effect of parameters on the type of equilibrium points, we vary the parameters a5,a7 and investigate the type of equilibrium point with the real part of eigenvalues and the Routh–Hurwitz condition δ1δ2>δ0. Figure 1(a) shows the change in the real part of eigenvalues of system (1) with parameter a5, and for a5<1, the real part of eigenvalues become negative making the equilibrium stable. Similarly, Figure 1(b) is change in the real part of eigenvalues with a7, and for a7>−0.6, the equilibrium is stable.Figure 1
Eigenvalues of the system for various values of (a)a5 and (b) a7.
(a)
(b)To analyze the stability of the systems using Routh–Hurwitz (RH) criterion, we investigate the two conditionsδ1,δ0>0andδ1δ2>δ0. For the parameter values a1=1,a2=1, a3=−1,a4=−4,anda6=1, the RH condition modifies to δ0=1,δ1=−a5a7,andδ2=−signa7. We vary the parameters a5 and a7 between 0.7,2 and −1.3,−0.6, respectively, and the condition δ1δ2>δ0 is plotted as shown in Figure 2 which confirms our claim that the equilibrium is stable when a7>−0.6 or a5<1. Both the eigenvalues and RH investigations confirm that system (1) exhibits both stable (nonhyperbolic) and unstable (hyperbolic) equilibrium points. To the best of our knowledge, this feature has not been investigated in the literature.Figure 2
Plots showing the RH conditionδ1δ2>δ0 vs. parameters a5anda7.To show the 2D phase portraits of the stable and unstable system, we fixa7=−1 and plot for a5=0.95 (stable equilibrium) and a5=1.6 (unstable equilibrium) as shown in Figures 3(a) and 3(b), respectively.Figure 3
2D phase portraits of the novel jerk system fora1=1,a2=1,a3=−1,a4=−4, and a6=1 and initial conditions −34−4. Figure 3(a) shows system (1) with a stable equilibrium for a5=0.95anda7=−1, and Figure 3(b) shows system (1) with a unstable equilibrium for a5=1.6anda7=−1.
(a)
(b)The finite-time Lyapunov exponents of system (1) are derived using Wolf’s algorithm [1] for both stable and unstable cases and are given in Table 1 with the Kaplan–Yorke dimension (DKY) calculated for a run time of 20,000 s with initial conditions −34−4 with other parameters except a5,a7 fixed to their respective chaotic values.Table 1
Lyapunov exponents andDKY for different parameters.
Parameters
Lyapunov exponents (LEs)
D
KY
a
5
=
0.95
,
a
7
=
−
1
0.1196
,
0
,
−
1.179
2.101
a
5
=
1.6
,
a
7
=
−
1
0.0727
,
0
,
−
1.375
2.052To further understand the complete dynamical behavior of the system, we derived the bifurcation plots with parametera5 as it governs the equilibrium points of the proposed system. The parameter a5 is varied between 0.8,2, and the local maximum of the state variable z is plotted as shown in Figures 4(a) and 4(b) which show the corresponding LEs of the system. It can be seen from Figure 4(a) that, under the change of parameter a5, the chaos state and periodic state of the system appear alternately, and the way the system changes from the periodic state to chaos is period-doubling bifurcation. By comparing Figures 4(a) and 4(b), the maximum Lyapunov exponent of the system is zero; When the system is chaotic, the largest Lyapunov exponent of the system is positive. When the system is convergent, the largest Lyapunov exponent of the system is negative.Figure 4
(a) Bifurcation of the system witha5; (b) corresponding LEs.
(a)
(b)


## 3. Electronic Circuit Implementation of the System

In this section, system (1) is implemented with an electronic circuit. The circuit designed in PSpice medium and its simulation are carried out. The schematic of the designed circuit is given in Figure 5. Here, we applied the traditional method by the multiplier AD633/AD to carry out the function zsgnx, and the new switch published in [25] was also tried by simulation which shows the same results.Figure 5
The designed circuit for system (1).The electronic circuit consists of OpAmps, multiplier ICs, diodes, resistors, and capacitors. The supply voltages of the active elements are +15 V for positive supply inputs and −15 V for negative supply inputs. The value of the resistors used in the circuit is as follows:R1=R2=R3=400kΩ, R4=R8=R9=R12=R13=R14=R15=R16=100kΩ, R5=42kΩ, R6=40kΩ, and R7=6MΩ. The values of the capacitor in the circuit are as follows: C1=C2=C3=1nF. The given resistor and capacitor values are for the stable equilibrium case when a5=0.95anda7=−1. For the unstable equilibrium case (when a5=1.6anda7=−1), only the value of R5 resistor will be changed. The value of R5 resistor is 25 kΩ.The real-time electronic circuit application is also realized. The results of the real-time electronic circuit are examined on the oscilloscope. 2D phase portraits of system (1) obtained from the oscilloscope are given in Figures 6(a) and 6(b). The phase portrait of the stable equilibrium case a5=0.95anda7=−1 is given in Figure 6(a) while that of the unstable equilibrium case a5=1.6anda7=−1 is given in Figure 6(b). Both of them confirm the theoretical prediction in Figure (3).Figure 6
2D phase portraits, obtained from the oscilloscope, of the novel system fora1=1,a2=1,a3=−1,a4=−4, and a6=1 and initial conditions −34−4 (a) when a5=0.95anda7=−1 and (b) when a5=1.6anda7=−1.
(a)
(b)


## 4. Application Discussion

### 4.1. Random Number Generator Design and NIST Tests

Random number generator (RNG) is one of the most important applications in which chaotic systems are used. RNGs are used in many different engineering fields, especially communication and cryptography. In this section, a new RNG design has been realized by using the new chaotic jerk system proposed in this study. The flow diagram of the method used in the design of RNG is given in Figure7.Figure 7
RNG design stages with the proposed 3D chaotic system.As shown in the flow diagram, the appropriate initial values, the number of bits taken from the state variables, and the step value for the RK4 algorithm are determined with the help of the new chaotic system. Then, the RK4 algorithm is executed to obtain discrete arrays of state variables and convert them into 32 bit binary arrays. Steps to generate the pseudo-random sequences are as follows.The system uses the RK4 algorithm and Matlab software to iterate 2,000 times for simulation.S1: the sequence obtained by the first iteration of the system data = [x (1), y (1), z (1)]; the numerical value is kept to the second decimal place.S2: data = data ∗ 100; data are expressed as a 32 bit binary number:data=b31b30,…,b2b132.S3: Letdata1=b7b6,…,b2b18, data1 for further test use.S4:n = n+1. Then, it is still stated that the sequence to be generated passes through S2–S4 until n = 2001.Random bit arrays are generated for NIST tests by taking the appropriate number of bits from these arrays. If the generated bit array passes these tests successfully, the RNG design is completed and ready for use in practical applications. Otherwise, the procedure is repeated by starting over and setting the initial values, the number of bits taken from the state variables, and the step value for RK4. These operations are carried out until the bit array obtained has passed all of the NIST tests. When the process is complete, an RNG is obtained which successfully passed through the NIST tests, which are the most accepted randomness tests in the literature.The most important randomness test which is accepted internationally is the NIST-800-22 statistical tests. In this test, there are 15 different subtests. An array consisting of at least 1000000 bits must pass all of these tests successfully. That is, thep values obtained at the end of the tests should be equal or greater than 0.001. Table 2 shows the results of the NIST tests of the arrays which are generated by taking 8 bits from the state variables of the new chaotic system in each iteration. These arrays have successfully passed all the tests. Thus, it has been shown that the proposed system can be used conveniently in RNG-based engineering applications. The arrays are obtained in three different ways, in which the bit arrays obtained from all state variables are subjected to XOR operation, the Y and Z state variables are subjected to XOR operation, or using the Z state variable alone.Table 2
3D chaotic system PRNG NIST-800-22 test results.
Statistical tests
p value (X ⊕ Y ⊕ Z)
p value (Y ⊕ Z)
p value (Z)
Result
Frequency (monobit) test
0.373465
0.017981
0.485177
Successful
Block-frequency test
0.932102
0.563486
0.720823
Successful
Cumulative-sum test
0.692331
0.029619
0.813864
Successful
Runs test
0.093974
0.417706
0.661739
Successful
Longest-run test
0.651016
0.544020
0.730122
Successful
Binary matrix rank test
0.489777
0.123925
0.328310
Successful
Discrete Fourier transform test
0.393422
0.291282
0.769023
Successful
Nonoverlapping template test
0.058054
0.045477
0.115145
Successful
Overlapping template test
0.245790
0.504454
0.546871
Successful
Maurer’s universal statistical test
0.550389
0.259675
0.138921
Successful
Approximate entropy test
0.173084
0.067221
0.503813
Successful
Random excursion test (x = −4)
0.617626
0.055541
0.042691
Successful
Random-excursion variant test (x = −9)
0.842508
0.227993
0.065864
Successful
Serial test-1
0.651632
0.454848
0.316808
Successful
Serial test-2
0.675012
0.598698
0.073441
Successful
Linear complexity test
0.073771
0.830488
0.259938
Successful


### 4.2. Voice Encryption Application

In this part of the study, a different application based on the new chaotic system has been realized. In this application, voice encoding and decoding are performed using the RNG designed in the previous section. The block diagram of the application is shown in Figure8.Figure 8
Encryption and decryption process.As can be seen from Figure9, these voice data are recorded with noise, and noise is added in the simulation; the voice data are converted into a binary form and then subjected to XOR operation with a random bit array from the RNG. The encrypted bit array is converted back to the float to get the encrypted voice. In the decoding process, the encrypted voice values are converted into binary form and subjected to XOR operation with bit arrays received from the RNG. The bit array obtained after the XOR operation is converted to a float, and the decoding process is completed.Figure 9
Original, encrypted, and decrypted voice waves, respectively.Figure9 shows the original, encrypted, and decrypted voice plots, respectively. It is seen that there is no similarity between the encrypted voice and the original voice, and that the original voice can be reacquired at the end of the decoding process. The frequency spectrums of the voice signals are shown in Figure 10. When the spectrums are examined, it is seen that the original and decrypted voice spectrums are the same, and the coded spectrum is completely different from the others and has a homogeneous distribution. This indicates that the encryption and decryption processes are quite good.Figure 10
Original, encrypted, and decrypted voice spectrums.


## 4.1. Random Number Generator Design and NIST Tests

Random number generator (RNG) is one of the most important applications in which chaotic systems are used. RNGs are used in many different engineering fields, especially communication and cryptography. In this section, a new RNG design has been realized by using the new chaotic jerk system proposed in this study. The flow diagram of the method used in the design of RNG is given in Figure7.Figure 7
RNG design stages with the proposed 3D chaotic system.As shown in the flow diagram, the appropriate initial values, the number of bits taken from the state variables, and the step value for the RK4 algorithm are determined with the help of the new chaotic system. Then, the RK4 algorithm is executed to obtain discrete arrays of state variables and convert them into 32 bit binary arrays. Steps to generate the pseudo-random sequences are as follows.The system uses the RK4 algorithm and Matlab software to iterate 2,000 times for simulation.S1: the sequence obtained by the first iteration of the system data = [x (1), y (1), z (1)]; the numerical value is kept to the second decimal place.S2: data = data ∗ 100; data are expressed as a 32 bit binary number:data=b31b30,…,b2b132.S3: Letdata1=b7b6,…,b2b18, data1 for further test use.S4:n = n+1. Then, it is still stated that the sequence to be generated passes through S2–S4 until n = 2001.Random bit arrays are generated for NIST tests by taking the appropriate number of bits from these arrays. If the generated bit array passes these tests successfully, the RNG design is completed and ready for use in practical applications. Otherwise, the procedure is repeated by starting over and setting the initial values, the number of bits taken from the state variables, and the step value for RK4. These operations are carried out until the bit array obtained has passed all of the NIST tests. When the process is complete, an RNG is obtained which successfully passed through the NIST tests, which are the most accepted randomness tests in the literature.The most important randomness test which is accepted internationally is the NIST-800-22 statistical tests. In this test, there are 15 different subtests. An array consisting of at least 1000000 bits must pass all of these tests successfully. That is, thep values obtained at the end of the tests should be equal or greater than 0.001. Table 2 shows the results of the NIST tests of the arrays which are generated by taking 8 bits from the state variables of the new chaotic system in each iteration. These arrays have successfully passed all the tests. Thus, it has been shown that the proposed system can be used conveniently in RNG-based engineering applications. The arrays are obtained in three different ways, in which the bit arrays obtained from all state variables are subjected to XOR operation, the Y and Z state variables are subjected to XOR operation, or using the Z state variable alone.Table 2
3D chaotic system PRNG NIST-800-22 test results.
Statistical tests
p value (X ⊕ Y ⊕ Z)
p value (Y ⊕ Z)
p value (Z)
Result
Frequency (monobit) test
0.373465
0.017981
0.485177
Successful
Block-frequency test
0.932102
0.563486
0.720823
Successful
Cumulative-sum test
0.692331
0.029619
0.813864
Successful
Runs test
0.093974
0.417706
0.661739
Successful
Longest-run test
0.651016
0.544020
0.730122
Successful
Binary matrix rank test
0.489777
0.123925
0.328310
Successful
Discrete Fourier transform test
0.393422
0.291282
0.769023
Successful
Nonoverlapping template test
0.058054
0.045477
0.115145
Successful
Overlapping template test
0.245790
0.504454
0.546871
Successful
Maurer’s universal statistical test
0.550389
0.259675
0.138921
Successful
Approximate entropy test
0.173084
0.067221
0.503813
Successful
Random excursion test (x = −4)
0.617626
0.055541
0.042691
Successful
Random-excursion variant test (x = −9)
0.842508
0.227993
0.065864
Successful
Serial test-1
0.651632
0.454848
0.316808
Successful
Serial test-2
0.675012
0.598698
0.073441
Successful
Linear complexity test
0.073771
0.830488
0.259938
Successful


## 4.2. Voice Encryption Application

In this part of the study, a different application based on the new chaotic system has been realized. In this application, voice encoding and decoding are performed using the RNG designed in the previous section. The block diagram of the application is shown in Figure8.Figure 8
Encryption and decryption process.As can be seen from Figure9, these voice data are recorded with noise, and noise is added in the simulation; the voice data are converted into a binary form and then subjected to XOR operation with a random bit array from the RNG. The encrypted bit array is converted back to the float to get the encrypted voice. In the decoding process, the encrypted voice values are converted into binary form and subjected to XOR operation with bit arrays received from the RNG. The bit array obtained after the XOR operation is converted to a float, and the decoding process is completed.Figure 9
Original, encrypted, and decrypted voice waves, respectively.Figure9 shows the original, encrypted, and decrypted voice plots, respectively. It is seen that there is no similarity between the encrypted voice and the original voice, and that the original voice can be reacquired at the end of the decoding process. The frequency spectrums of the voice signals are shown in Figure 10. When the spectrums are examined, it is seen that the original and decrypted voice spectrums are the same, and the coded spectrum is completely different from the others and has a homogeneous distribution. This indicates that the encryption and decryption processes are quite good.Figure 10
Original, encrypted, and decrypted voice spectrums.


## 5. Conclusion and Discussion

Signum function and absolute value function can remove the amplitude or polarity information in the feedback variable, which sometimes still preserves the basic property in a dynamical system. As a new case, such a hypogenetic jerk system is obtained, which has stable or unstable equilibrium points under different parameters. Moreover, physical experiments prove the chaotic oscillation. Further discussion focuses on random number generator, and voice encryption shows that the new derived chaotic system still exhibits great noise-like randomness.



---

*Source: 8083509-2020-05-04.xml*