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

Authors: Abir Lassoued; Olfa Boubaker
Journal: Complexity
 (2017)

---

## Abstract

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.


---

## Body

## 1. 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–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–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 Section2, 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.


## 2. Mathematical Model and Basic Properties

Let consider the mathematical model of the novel hyperchaotic system with FO terms expressed by the following differential equations:(1)x˙=y,y˙=z,z˙=-az-by+Gxw˙=ky-hw+Gx,with G(x) being a nonlinear function defined as (2)Gx=-cx2+dxx+mxrx-1,1<r<2,where (a, b, c, d, h, k, m, r) are the system’s parameters and (x, y, z, w) are the state variables. r is a fractional number satisfying 1<r<2. Since r≠1, |x|rx-1 will never be an indeterminate form. The nonlinear function G(x) can be written as follows: (3)Gx=-c-dx2-m-xr-1,ifx<0,0,ifx=0,-c+dx2+mxr-1,ifx>0.System (1) can exhibit chaotic behavior if the general condition of dissipativity is satisfied such as (4)∂x˙∂x+∂y˙∂y+∂z˙∂z+∂w˙∂w=-a-h<0.As long as a+h>0, system (1) is dissipative and it converges to an attractor. Thus, when the parameters (a, b, c, d, m, r, h, k) 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.Figure 1
Projections of the attractor related to the hyperchaotic system (1) onto the spaces (x,y,z), (x,y), (x,z), (y,z), and (w,y).The time series of the state variablesx and w 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 w is extended within [-250,200], unlike the other variables. This point must be considered in practical applications.Figure 2
Time series and sensitive dependence on initial conditions(1,1,1,1) (blue) and (1.001,1,1,1) (red): (a) and (d) variable x; (b) and (c) variable w.
(a)
(b)
(c)
(d)


## 3. Dynamic Analysis

### 3.1. Equilibrium and Stability

The equilibrium points of system (1) are obtained by solving these equations: (5)y=z=w=0,-cx2+dxx+mxrx-1=0.Proposition 1.
(i) Ifx=0, then w=0 and the origin H1=(0,0,0,0) is the first equilibrium of system (1).
(ii) Ifx>0, then w=0 and H2=(((c-d)/m)1/α,0,0,0) is an equilibrium of system (2) where α=r-3.
(iii) Ifx<0, then w=0 and H3=((-(-c-d)/m)1/α,0,0,0) is an equilibrium of system (2) where α=r-3.Proof.
Case (i) is obvious.
For case (ii), we should solve the following equation:-c+d+mxr-3=0 which admits the solution xα=(c-d)/m with α=r-3. The roots of this equation are given by [22] (6)x=c-dm1/αejθ1±2nΠ/αwith n∈N, θ1 being the phase of xα, and α being a fractional number. Notice that the term (c-d)/m is positive when c=-0.11, d=-0.21, and m=6.26. Then, θ1 is equal to zero and we have x=((c-d)/m)1/α.
For case (iii), we should solve the following equation:-c-d-m(-x)r-3=0 which admits the solution (-x)α=Xα=(-c-d)/m with α=r-3. The roots of this equation are given by [22] (7)X=-c-dm1/αejθ2±2nΠ/αwith n∈N, θ2 being the phase of xα, and α being a fractional number. Notice that the term (-c-d)/m is positive when c=-0.11, d=-0.21, and m=6.26. Then, θ2 is equal to zero and we have x=(-(-c-d)/m)1/α.
When the parameters (a, b, c, d, m, r, h, k) are equal to (0.93, 1.11, -0.11, -0.21, 6.26, 1.32, 0.001, 14), system (1) admits three equilibrium points: H1=(0,0,0,0), H2=(11.73,0,0,0), and H3=(-5.87,0,0,0). For the stability analysis, Table 1 gives the Jacobian matrix J and its corresponding eigenvalues calculated for each equilibrium point.Table 1
Stability analysis of system (1).
Equilibrium point
Jacobian matrix
Corresponding eigenvalues
Stability analysis
H1
0
1
0
0
0
0
1
0
0
-
b
-
a
0
0
k
0
-
h
λ
1
=
0  λ2=-0.001  λ3=-0.465+0.945i  λ4=-0.465-0.945i
Stable point
H2
0
1
0
0
0
0
1
0
Δ
1
-
b
-
a
0
Δ
1
k
0
-
h
λ
1
=
-
0.001  λ2=-1.273  λ3=0.171+1.232i  λ4=0.171-1.232i
Unstable point
Δ
1
=
2
d
x
-
2
c
x
+
(
r
-
1
)
m
x
r
-
2
H
3
0
1
0
0
0
0
1
0
Δ
2
-
b
-
a
0
Δ
2
k
0
-
h
λ
1
=
-
0.001  λ2=-2.297  λ3=0.683+1.945i  λ4=0.683-1.945i
Unstable point
Δ
2
=
-
2
d
x
-
2
c
x
-
(
r
-
1
)
m
x
r
-
2


### 3.2. 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 (8)LE1=0.232,LE2=0.020,LE3=0,LE4=-1.169.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.Figure 3
The Lyapunov exponent spectrum of the hyperchaotic system (1).In addition, the corresponding Kaplan-Yorke dimension is(9)DL=3+λ1+λ2+λ3λ4=3.19.


### 3.3. Routes to Chaos

System (1) can display periodic orbits, chaos, and hyperchaos attractors under different conditions. In fact, when the parameter m varies and the parameters (a, b, c, d, r, h, k) are fixed, two Hopf bifurcations are detected as shown in Figure 4. These bifurcations are denoted H in the bifurcation diagram and they appear when m=0.99 and m=2.11, respectively. Each Hopf point is characterized by a 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.Figure 4
Hopf bifurcations.In addition, as the parameter of bifurcationm increases, system (1) undergoes the following routes:(i)
If-1≤m≤0.3, then system (1) exhibits periodic orbit. Figure 5(a) shows this regular attractor with m=-0.5.(ii)
If0.3<m≤2.1, then system (1) converges to a fixed point as shown in Figure 5(b).(iii)
If2.1<m≤3.8, then another periodic orbit is obtained as shown in Figure 5(c) with m=1.(iv)
If3.8<m≤5.5, system (1) exhibits chaotic attractor. Figure 5(d) shows this strange attractor with m=4.5.(v)
If5.5<m<7, then system (1) exhibits hyperchaotic attractor. Figure 5(e) shows this strange attractor with m=6.Figure 5
Different attractors forms exhibited by system (1) when parameter m varies.
(a)
(b)
(c)
(d)
(e)Some typical attractors are tabulated in Table2 according to the parameter m.Table 2
The LEs of some typical attractors of system (1).
m
LE1
LE2
LE3
L
E
4
Attractor
3.6
0
0
-
0.42
0.8
Periodic orbit
3.8
0.05
0
-0.01
-0.92
Chaotic attractor
6
0.18
0.007
0
-1.10
Hyperchaotic attractor
6.5
0.17
0.01
0
-1.16
Hyperchaotic attractor


### 3.4. 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.Table 3
Comparative analysis with related hyperchaotic systems.
Hyperchaotic system
Lyapunov exponents
Kaplan-Yorke dimension
Proposed hyperchaotic system
L
E
1
=
0.231
DKY=3.19
L
E
2
=
0.020
L
E
3
=
0
L
E
4
=
-
1.169
Piecewise linear hyperchaotic circuit [25]
L
E
1
=
0.064
DKY=3.089
L
E
2
=
0.033
L
E
3
=
0
L
E
4
=
-
1.098
Hyperchaotic hyperjerk system [26]
L
E
1
=
0.142
DKY=3.134
L
E
2
=
0.046
L
E
3
=
0
L
E
4
=
-
1.396


## 3.1. Equilibrium and Stability

The equilibrium points of system (1) are obtained by solving these equations: (5)y=z=w=0,-cx2+dxx+mxrx-1=0.Proposition 1.
(i) Ifx=0, then w=0 and the origin H1=(0,0,0,0) is the first equilibrium of system (1).
(ii) Ifx>0, then w=0 and H2=(((c-d)/m)1/α,0,0,0) is an equilibrium of system (2) where α=r-3.
(iii) Ifx<0, then w=0 and H3=((-(-c-d)/m)1/α,0,0,0) is an equilibrium of system (2) where α=r-3.Proof.
Case (i) is obvious.
For case (ii), we should solve the following equation:-c+d+mxr-3=0 which admits the solution xα=(c-d)/m with α=r-3. The roots of this equation are given by [22] (6)x=c-dm1/αejθ1±2nΠ/αwith n∈N, θ1 being the phase of xα, and α being a fractional number. Notice that the term (c-d)/m is positive when c=-0.11, d=-0.21, and m=6.26. Then, θ1 is equal to zero and we have x=((c-d)/m)1/α.
For case (iii), we should solve the following equation:-c-d-m(-x)r-3=0 which admits the solution (-x)α=Xα=(-c-d)/m with α=r-3. The roots of this equation are given by [22] (7)X=-c-dm1/αejθ2±2nΠ/αwith n∈N, θ2 being the phase of xα, and α being a fractional number. Notice that the term (-c-d)/m is positive when c=-0.11, d=-0.21, and m=6.26. Then, θ2 is equal to zero and we have x=(-(-c-d)/m)1/α.
When the parameters (a, b, c, d, m, r, h, k) are equal to (0.93, 1.11, -0.11, -0.21, 6.26, 1.32, 0.001, 14), system (1) admits three equilibrium points: H1=(0,0,0,0), H2=(11.73,0,0,0), and H3=(-5.87,0,0,0). For the stability analysis, Table 1 gives the Jacobian matrix J and its corresponding eigenvalues calculated for each equilibrium point.Table 1
Stability analysis of system (1).
Equilibrium point
Jacobian matrix
Corresponding eigenvalues
Stability analysis
H1
0
1
0
0
0
0
1
0
0
-
b
-
a
0
0
k
0
-
h
λ
1
=
0  λ2=-0.001  λ3=-0.465+0.945i  λ4=-0.465-0.945i
Stable point
H2
0
1
0
0
0
0
1
0
Δ
1
-
b
-
a
0
Δ
1
k
0
-
h
λ
1
=
-
0.001  λ2=-1.273  λ3=0.171+1.232i  λ4=0.171-1.232i
Unstable point
Δ
1
=
2
d
x
-
2
c
x
+
(
r
-
1
)
m
x
r
-
2
H
3
0
1
0
0
0
0
1
0
Δ
2
-
b
-
a
0
Δ
2
k
0
-
h
λ
1
=
-
0.001  λ2=-2.297  λ3=0.683+1.945i  λ4=0.683-1.945i
Unstable point
Δ
2
=
-
2
d
x
-
2
c
x
-
(
r
-
1
)
m
x
r
-
2


## 3.2. 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 (8)LE1=0.232,LE2=0.020,LE3=0,LE4=-1.169.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.Figure 3
The Lyapunov exponent spectrum of the hyperchaotic system (1).In addition, the corresponding Kaplan-Yorke dimension is(9)DL=3+λ1+λ2+λ3λ4=3.19.


## 3.3. Routes to Chaos

System (1) can display periodic orbits, chaos, and hyperchaos attractors under different conditions. In fact, when the parameter m varies and the parameters (a, b, c, d, r, h, k) are fixed, two Hopf bifurcations are detected as shown in Figure 4. These bifurcations are denoted H in the bifurcation diagram and they appear when m=0.99 and m=2.11, respectively. Each Hopf point is characterized by a 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.Figure 4
Hopf bifurcations.In addition, as the parameter of bifurcationm increases, system (1) undergoes the following routes:(i)
If-1≤m≤0.3, then system (1) exhibits periodic orbit. Figure 5(a) shows this regular attractor with m=-0.5.(ii)
If0.3<m≤2.1, then system (1) converges to a fixed point as shown in Figure 5(b).(iii)
If2.1<m≤3.8, then another periodic orbit is obtained as shown in Figure 5(c) with m=1.(iv)
If3.8<m≤5.5, system (1) exhibits chaotic attractor. Figure 5(d) shows this strange attractor with m=4.5.(v)
If5.5<m<7, then system (1) exhibits hyperchaotic attractor. Figure 5(e) shows this strange attractor with m=6.Figure 5
Different attractors forms exhibited by system (1) when parameter m varies.
(a)
(b)
(c)
(d)
(e)Some typical attractors are tabulated in Table2 according to the parameter m.Table 2
The LEs of some typical attractors of system (1).
m
LE1
LE2
LE3
L
E
4
Attractor
3.6
0
0
-
0.42
0.8
Periodic orbit
3.8
0.05
0
-0.01
-0.92
Chaotic attractor
6
0.18
0.007
0
-1.10
Hyperchaotic attractor
6.5
0.17
0.01
0
-1.16
Hyperchaotic attractor


## 3.4. 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.Table 3
Comparative analysis with related hyperchaotic systems.
Hyperchaotic system
Lyapunov exponents
Kaplan-Yorke dimension
Proposed hyperchaotic system
L
E
1
=
0.231
DKY=3.19
L
E
2
=
0.020
L
E
3
=
0
L
E
4
=
-
1.169
Piecewise linear hyperchaotic circuit [25]
L
E
1
=
0.064
DKY=3.089
L
E
2
=
0.033
L
E
3
=
0
L
E
4
=
-
1.098
Hyperchaotic hyperjerk system [26]
L
E
1
=
0.142
DKY=3.134
L
E
2
=
0.046
L
E
3
=
0
L
E
4
=
-
1.396


## 4. 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).


### 4.1. Design of the Analog Circuit with MultiSIM

For the circuit implementation, we choose the particular case study when the system parameterr is fixed to 1.5. Thus, the proposed system will be defined by the following model: (10)x˙=y,y˙=z,z˙=-az-by-cx2+dxx+mxsgnxw˙=ky-hw-cx2+dxx+mxsgnx,where the system parameters (a, b, c, d, m, h, k) are equal to (1, 1, -0.11, -0.21, 5, 0.01,14). System (10) exhibits four LEs such as (11)LE1=0.18,LE2=0.04,LE3=0,LE4=-1.2.Despite the observation of the obtained phase portraits of system (10), we deduce that the maximum value of the signal w can reach the level of 250. Thus, 250V 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 u=x/2, v=y/2, g=z/2, and f=w/160 and then setting the original state variables x, y, z, and w instead of the variable u, v, g, and f, the adjusted system becomes the following one: (12)x˙=y,y˙=z,z˙=az-by-2cx2+2dxx+m2xsgnxw˙=ky-80hw-2cx2+2dxx+m2xsgnx.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.Figure 6
Projections of the attractor related to the adjusted hyperchaotic system (12) onto the spaces (x,w), (x,z), and (y,z).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.Figure 7
Circuit design of the square root function with MultiSIM.For the theoretical study and based on [27], the second voltage source in Figure 7 should be fixed to 2.878V. However, in experimentation applications, we have obtained the root square function by using a stabilized voltage equal to 2.9V 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 720mV and 900mV (≃√0.72 = 0.88), respectively. Thus, the square root function is correctly obtained with 2.9V. In addition, based on MultiSIM results and experimental simulations, if the source voltage is included in [2.7V,3V] 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.Figure 8
Experimental results of the square root function.The corresponding circuit equation of the hyperchaotic system can be described as(13)x˙=1R1C1y,y˙=1R2C2z,z˙=-1R3C3z-1R4C3y+1R5C3x2-1R6C3xx+1R7C3xsgnxw˙=-1R8C4y-1R9C4w+1R5C4x2-1R6C4xx+1R7C4xsgnx.According to system (12) and (13) and design considerations, we fixed the values of the resistances and the capacitors as (14)C1=C2=C3=1nF,C4=0.1nF,R1=R2=100kΩ,R3=R4=10kΩ,R5=4.554kΩ,R6=2.5kΩ,R7=0.2kΩ,R8=7kΩ,R9=1.25kΩ.Finally, the obtained circuit diagram, designed with MultiSIM Software, is provided in Figure 9 where the multiplier is AD633 and the operator amplifier is UA741.Figure 9
Circuit design of the hyperchaotic system with MultiSIM.


### 4.2. Simulation Results

For the oscillator circuit, all active devices (UA741 and AD633) are powered by±15V. 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 Figure10. 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.Figure 10
Simulation results of the hyperchaotic system with MultiSIM Software.


## 4.1. Design of the Analog Circuit with MultiSIM

For the circuit implementation, we choose the particular case study when the system parameterr is fixed to 1.5. Thus, the proposed system will be defined by the following model: (10)x˙=y,y˙=z,z˙=-az-by-cx2+dxx+mxsgnxw˙=ky-hw-cx2+dxx+mxsgnx,where the system parameters (a, b, c, d, m, h, k) are equal to (1, 1, -0.11, -0.21, 5, 0.01,14). System (10) exhibits four LEs such as (11)LE1=0.18,LE2=0.04,LE3=0,LE4=-1.2.Despite the observation of the obtained phase portraits of system (10), we deduce that the maximum value of the signal w can reach the level of 250. Thus, 250V 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 u=x/2, v=y/2, g=z/2, and f=w/160 and then setting the original state variables x, y, z, and w instead of the variable u, v, g, and f, the adjusted system becomes the following one: (12)x˙=y,y˙=z,z˙=az-by-2cx2+2dxx+m2xsgnxw˙=ky-80hw-2cx2+2dxx+m2xsgnx.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.Figure 6
Projections of the attractor related to the adjusted hyperchaotic system (12) onto the spaces (x,w), (x,z), and (y,z).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.Figure 7
Circuit design of the square root function with MultiSIM.For the theoretical study and based on [27], the second voltage source in Figure 7 should be fixed to 2.878V. However, in experimentation applications, we have obtained the root square function by using a stabilized voltage equal to 2.9V 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 720mV and 900mV (≃√0.72 = 0.88), respectively. Thus, the square root function is correctly obtained with 2.9V. In addition, based on MultiSIM results and experimental simulations, if the source voltage is included in [2.7V,3V] 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.Figure 8
Experimental results of the square root function.The corresponding circuit equation of the hyperchaotic system can be described as(13)x˙=1R1C1y,y˙=1R2C2z,z˙=-1R3C3z-1R4C3y+1R5C3x2-1R6C3xx+1R7C3xsgnxw˙=-1R8C4y-1R9C4w+1R5C4x2-1R6C4xx+1R7C4xsgnx.According to system (12) and (13) and design considerations, we fixed the values of the resistances and the capacitors as (14)C1=C2=C3=1nF,C4=0.1nF,R1=R2=100kΩ,R3=R4=10kΩ,R5=4.554kΩ,R6=2.5kΩ,R7=0.2kΩ,R8=7kΩ,R9=1.25kΩ.Finally, the obtained circuit diagram, designed with MultiSIM Software, is provided in Figure 9 where the multiplier is AD633 and the operator amplifier is UA741.Figure 9
Circuit design of the hyperchaotic system with MultiSIM.


## 4.2. Simulation Results

For the oscillator circuit, all active devices (UA741 and AD633) are powered by±15V. 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 Figure10. 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.Figure 10
Simulation results of the hyperchaotic system with MultiSIM Software.


## 5. 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 MultiSIM 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.



---

*Source: 3273408-2017-10-26.xml*