# Dynamic Analysis and Circuit Implementation of a New 4D Lorenz-Type Hyperchaotic System

Authors: A. Al-khedhairi; A. Elsonbaty; A. H. Abdel Kader; A. A. Elsadany
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

This paper attempts to further extend the results of dynamical analysis carried out on a recent 4D Lorenz-type hyperchaotic system while exploring new analytical results concerns its local and global dynamics. In particular, the equilibrium points of the system along with solution’s continuous dependence on initial conditions are examined. Then, a detailedZ2 symmetrical Bogdanov-Takens bifurcation analysis of the hyperchaotic system is performed. Also, the possible first integrals and global invariant surfaces which exist in system’s phase space are analytically found. Theoretical results reveal the rich dynamics and the complexity of system behavior. Finally, numerical simulations and a proposed circuit implementation of the hyperchaotic system are provided to validate the present analytical study of the system.


---

## Body

## 1. Introduction

Nonlinear dynamics analysis of various phenomena and systems of physics, engineering, biology, chemistry, economy, and industry has attracted a great interest among scientists and considered a very active area of research from 1960s in the last century till now [1–4]. There are two key reasons which interpret this great interest. The first one is that the dynamical systems tools help scientists better comprehend and analyze the varieties of nonlinear characteristics and new phenomena exhibited by systems from different disciplines. In particular, the tools of dynamical systems such as the applied bifurcation theories are successfully employed to investigate the qualitative behaviors of nonlinear systems [5–7]. This includes investigation of equilibrium points and their stability, creation, destruction and stability of periodic orbits, quasiperiodic behavior, homoclinic orbits, creation or destruction of chaotic attractors, and chaos control and synchronization. The second reason is that engineers and scientists can utilize some of the fascinating features of nonlinear dynamical systems in wide range of interesting applications.The spellbinding chaotic dynamics, as an example, is recognized by high sensitivity to initial conditions and positive Lyapunov exponents. The generation of chaos for practical applications can be achieved by exploiting nonlinear electronic circuits, nonlinearities of laser systems with feedbacks, or via digital platforms such as FPGAs and DSPs. The noise-like behavior of chaotic systems, their wide spectrum, and the possibility of attaining chaos synchronization between two chaotic systems render them essential for cutting edge applications related to cryptography and robust physical-layer secure communication systems and ultra-fast physical random bits generation [4–11]. The other applications of dynamical systems methods and chaos theory include financial systems, mathematical biology, nonlinear circuits, nonlinear mechanical systems, plasma physics, chaos control, efficient image encryption, neuroscience research, and geophysics [12–33].The more complex hyperchaotic system possesses at least two positive Lyapunov exponents and it has a phase space of dimension at least four. Clearly, the hyperchaotic systems have more randomness and higher unpredictability than simple chaotic systems. Therefore, hyperchaos is more preferred than simple chaos and its applications have recently become a central topic in research including chaos-based secure communications, image encryption, and cryptographyIn the last two decades, some interesting high-dimensional hyperchaotic systems, in science and engineering, have been explored and their dynamics have been extensively investigated [18, 24, 27, 28, 34, 35]. In fact, it is of great importance from theoretical and practical aspects to explain complicated phenomena and internal structural characteristics of hyperchaotic systems. This research line focuses on applying codimension two or three bifurcation analysis to system under investigation. For example, analysis of Bogdanov-Takens bifurcation, degenerate Hopf bifurcation, and Heteroclinic and Homoclinic bifurcations can be undertaken. Also, it is crucial to determine whether there exists any set of parameters values for which the considered system is integrable and find the corresponding invariant surfaces if exist. As it is practically unattainable goal to apply exhaustive numerical investigation to acquire values of parameters where the hyperchaotic system is integrable, one needs a powerful analytical method that enables achieving this goal easily.Recently, theoretical analysis, based on Fishing Principle, is applied to study some global features of a new four-dimensional Lorenz-type hyperchaotic system [36]. More specifically, conditions for existence of homoclinic orbits are obtained. This work aims at extending the aforementioned work and exploring other aspects of complicated dynamical behaviors of the 4D Lorenz-type hyperchaotic system. Analytical bifurcation structure of the model which includes Z2 symmetrical Bogdanov-Takens, Pitchfork, Andronov-Hopf bifurcation and homoclinic bifurcation are obtained. Existence, uniqueness, and continuous dependence on initial conditions for the solution of 4D hyperchaotic system are all investigated. It is important to ask whether there is any set of parameter values in which system dynamics are regular and the studied system is integrable. The different integrable cases of the system are studied and the closed-form expressions for invariant surfaces corresponding to first integrals of the system are also obtained. Moreover, the hyperchaotic system has been simulated using a proposed electronic circuit realization and numerical simulations are performed confirming the new results of theoretical analysis. Note that attractors found in dynamical system are classified as being self-excited or hidden attractors [37, 38]. The key difference between them is that the self-excited attractor has a basin of attraction excited from unstable equilibria. The hidden attractor, on the other hand, has a basin of attraction with no intersection with neighborhoods of dynamical systems equilibrium points. In this work, we focus on self-excited attractors exist in the hyperchaotic system.The rest of the paper is structured as follows: In Section2, we introduce the 4D Lorenz-type system and discuss the equilibrium points’ existence. A sufficient condition for continuous dependence on initial condition is determined. Phase portraits, bifurcation diagrams, and Lyapunov characteristics spectrum are obtained. In Section 3, the analysis of some possible codimension two bifurcations is performed. It is shown that the 4D Lorenz system undergoes Bogdanov-Takens bifurcation, Andronov-Hopf bifurcation, Pitchfork bifurcation, and homoclinic bifurcation. The integrability analysis of system is investigated in Section 4. A practical application to engineering will be realized by an electronic circuit in Section 5. Finally, Section 6 concludes the paper.


## 2. The 4D Lorenz-Type Hyperchaotic System

The following 4D hyperchaotic system of Lorenz type was presented in [36](1)x˙=ay-x,y˙=cx-dy-xz,z˙=-bz+xy+w,w˙=-rw+kz,where state variables of the system are denoted by x,y,z, and w, the parameters of the system are represented by a,b,c,d,r, and k, and the dot above state variables refers to time derivative of state variables.The following subsections examine the main properties of system (1) and provide elementary dynamical analysis of the model.


### 2.1. Equilibrium Points of the System

The fixed points of system (1) can be obtained when x˙=y˙=z˙=w˙=0, i.e., through solving(2)ay-x=0,cx-dy-xz=0,-bz+xy+w=0,-rw+kz=0.Therefore, system (1) has the following equilibrium solutions:(I) when(d-c)(k-br)/r≤0, the fixed point is given by (3)x=y=z=w=0,(II) when(d-c)(k-br)/r>0, the fixed points are(4)x=y=z=w=0,x=y=±d-ck-brr,z=c-d,w=kc-dr.


### 2.2. Existence and Uniqueness of the Solution

The hyperchaotic system (1) can be put in the form(5)X˙t=ΦXt,t∈0,T,where(6)X=x1x2x3x4,X0=x01x02x03x04,ΦX=ax2-x1,cx1-dx2-x1x3-bx3+x1x2+x4-rx4+kx3,and the initial conditions are given by(7)X0=X0.For the class of continuous functionsF(t)∈C[0,T], we use the following norm in subsequent analysis:(8)F=supt∈0,TFt,while the matrix K=[kij[t]] of continuous functions employs the norm(9)K=∑i,jsupt∈0,Tkijt.It is obvious that 4D system (1) is dissipative if ∇.F(X)=-(a+b+d+r)<0. Now, solution of the system is examined in specific region Γ×J where J=[0,T] and(10)Γ=x1,x2,x3,x4:max⁡x,x2,x3andx4≤M,M>0.Parameter M is utilized to lay a boundary for the phase space region where existence and uniqueness of the solution are investigated.The solution of (5) and (7) can be represented by(11)Xt=X0+∫0tFXsds.The equivalence of the integral equation (11) and system (5)-(7) is obvious. Now, denoting the right hand side of (11) by H(X), then for X1=x11x12x13x14 and X2=x21x22x23x24 we get(12)HX1-HX2=∫0tΦX1s-ΦX2sds,and therefore (13)HX1-HX2≤∫0tΦX1s-ΦX2sds.Then, we obtain (14)HX1-HX2≤Tmax⁡a+c+2M,a+d+M,b+k+M,1+rX1-X2≤KX1-X2,where(15)K=Tmax⁡a+c+2M,a+d+M,b+k+M,1+r>0.Thus, X=H(X), for 0<K<1 as sufficient condition, is a contraction mapping.Theorem 1.
Assume that0<Tmax⁡a+c+2M,a+d+M,b+k+M,1+r<1, then a unique solution of hyperchaotic system (1), starting from initial condition X(0)=X0 in the region Ω×J, exists.


### 2.3. Continuous Dependence on Initial Conditions

The continuous dependence on initial conditions means that solution trajectories of the system which start close to each other still close to each other with evolution of time. This property is contrary to sensitive dependence on initial conditions which specifies chaotic dynamics. The goal of the next analysis is to find the particular parameters and time range where continuous dependence on state variables initial conditions is persevered; i.e., system (1) does not exhibit chaotic dynamics.Assume that there are two points of initial conditions of system (5)-(7); namely, X01 and X02 satisfy (16)X01-X02≤δ.First, suppose that the condition of Theorem 1 holds. Thus (17)X1t=X01+∫0tHX1sds,X2t=X02+∫0tHX2sds,and also we obtain (18)X1-X2≤X01-X02+KX1-X2,and (19)1-KX1-X2≤X01-X02,where (20)0<K<1is defined by (15). Finally, we get (21)X1-X2≤ϵ,where ϵ=δ/(1-K). Thus, we can formulate the following theorem.Theorem 2.
The solution of hyperchaotic system (1) exhibits continuous dependence on initial conditions if system (1) satisfies the condition of Theorem 1. More specifically, ∀ϵ>0∃δ(ϵ)=(1-K)ϵ>0 such that whenever X01-X02≤δ it implies that X1-X2≤ϵ.


## 2.1. Equilibrium Points of the System

The fixed points of system (1) can be obtained when x˙=y˙=z˙=w˙=0, i.e., through solving(2)ay-x=0,cx-dy-xz=0,-bz+xy+w=0,-rw+kz=0.Therefore, system (1) has the following equilibrium solutions:(I) when(d-c)(k-br)/r≤0, the fixed point is given by (3)x=y=z=w=0,(II) when(d-c)(k-br)/r>0, the fixed points are(4)x=y=z=w=0,x=y=±d-ck-brr,z=c-d,w=kc-dr.


## 2.2. Existence and Uniqueness of the Solution

The hyperchaotic system (1) can be put in the form(5)X˙t=ΦXt,t∈0,T,where(6)X=x1x2x3x4,X0=x01x02x03x04,ΦX=ax2-x1,cx1-dx2-x1x3-bx3+x1x2+x4-rx4+kx3,and the initial conditions are given by(7)X0=X0.For the class of continuous functionsF(t)∈C[0,T], we use the following norm in subsequent analysis:(8)F=supt∈0,TFt,while the matrix K=[kij[t]] of continuous functions employs the norm(9)K=∑i,jsupt∈0,Tkijt.It is obvious that 4D system (1) is dissipative if ∇.F(X)=-(a+b+d+r)<0. Now, solution of the system is examined in specific region Γ×J where J=[0,T] and(10)Γ=x1,x2,x3,x4:max⁡x,x2,x3andx4≤M,M>0.Parameter M is utilized to lay a boundary for the phase space region where existence and uniqueness of the solution are investigated.The solution of (5) and (7) can be represented by(11)Xt=X0+∫0tFXsds.The equivalence of the integral equation (11) and system (5)-(7) is obvious. Now, denoting the right hand side of (11) by H(X), then for X1=x11x12x13x14 and X2=x21x22x23x24 we get(12)HX1-HX2=∫0tΦX1s-ΦX2sds,and therefore (13)HX1-HX2≤∫0tΦX1s-ΦX2sds.Then, we obtain (14)HX1-HX2≤Tmax⁡a+c+2M,a+d+M,b+k+M,1+rX1-X2≤KX1-X2,where(15)K=Tmax⁡a+c+2M,a+d+M,b+k+M,1+r>0.Thus, X=H(X), for 0<K<1 as sufficient condition, is a contraction mapping.Theorem 1.
Assume that0<Tmax⁡a+c+2M,a+d+M,b+k+M,1+r<1, then a unique solution of hyperchaotic system (1), starting from initial condition X(0)=X0 in the region Ω×J, exists.


## 2.3. Continuous Dependence on Initial Conditions

The continuous dependence on initial conditions means that solution trajectories of the system which start close to each other still close to each other with evolution of time. This property is contrary to sensitive dependence on initial conditions which specifies chaotic dynamics. The goal of the next analysis is to find the particular parameters and time range where continuous dependence on state variables initial conditions is persevered; i.e., system (1) does not exhibit chaotic dynamics.Assume that there are two points of initial conditions of system (5)-(7); namely, X01 and X02 satisfy (16)X01-X02≤δ.First, suppose that the condition of Theorem 1 holds. Thus (17)X1t=X01+∫0tHX1sds,X2t=X02+∫0tHX2sds,and also we obtain (18)X1-X2≤X01-X02+KX1-X2,and (19)1-KX1-X2≤X01-X02,where (20)0<K<1is defined by (15). Finally, we get (21)X1-X2≤ϵ,where ϵ=δ/(1-K). Thus, we can formulate the following theorem.Theorem 2.
The solution of hyperchaotic system (1) exhibits continuous dependence on initial conditions if system (1) satisfies the condition of Theorem 1. More specifically, ∀ϵ>0∃δ(ϵ)=(1-K)ϵ>0 such that whenever X01-X02≤δ it implies that X1-X2≤ϵ.


## 3.Z2 Symmetry Bogdanov-Takens Bifurcation

Now, we study the case wherec=c∗=-a and d=d∗=-a that implies Jacobian matrix has two real zero eigenvalues, λ1,2=0, and two negative eigenvalues if either b2-2br+4k+r2<b+r or b2-2br+4k+r2<0 evaluated at the origin. The following coordinates’ transformation is applied to (1) in order to put the system in standard form(22)x1x2x3x4=1-1a00100000-b+r-ϖ2k-b+r+b2-2rb+r2+4k2k0011u1u2u3u4,which yields(23)u˙1u˙2u˙3u˙4=ϝ1ϝ2ϝ3ϝ4,where(24)ϝ1=y4y2ϖ2ak-y3y2ϖ2ak-by3y22ak-by4y22ak+ry3y22ak+ry4y22ak+y1y3ϖ2k-y1y4ϖ2k+by1y32k+by1y42k-ry1y32k-ry1y42k+y2,ϝ2=ay1y3ϖ2k-ay1y4ϖ2k+aby1y32k+aby1y42k-ary1y32k-ary1y42k-y2y3ϖ2k+y2y4ϖ2k-by2y32k-by2y42k+ry2y32k+ry2y42k,ϝ3=ky1y2aϖ-b2y34ϖ-b2y44ϖ+bry32ϖ+bry42ϖ-14y3ϖ+14y4ϖ-ky12ϖ-ky3ϖ-ky4ϖ-r2y34ϖ-r2y44ϖ-by32-ry32,ϝ4=-ky1y2aϖ+b2y34ϖ+b2y44ϖ-bry32ϖ-bry42ϖ+ky12ϖ-14y3ϖ+ky3ϖ+r2y34ϖ+14y4ϖ+ky4ϖ+r2y44ϖ-by42-ry42,such that ϖ=b2-2br+4k+r2.The center manifold is assumed in the form of second-order polynomial for sufficiently smallu1 and u2 as(25a)u3=α20u12+α11u1u2+α02u22,(25b)u4=β20u12+β11u1u2+β02u22.then substituting from ((25a)-(25b)) and the first two equations of (23) into the last two equation of (23). By comparing the coefficient of u2iu3j, i+j=2 in both sides of last two equations of (23) after substitution in ((25a)-(25b)), the values of center manifold coefficients are obtained as follows: (26)α20=-2kb2+bϖ-2r+rϖ+r+4k,α11=2k4aϖ+b2+bϖ-2r+rϖ+r+4kab2+bϖ-2r+rϖ+r+4k2,α02=-4kb-r2+4k4a+ϖ+b+rab2+bb-r2+4k-2r+rϖ+r+4k3,β20=-kϖ+b+r2ϖk-br,β11=k-4a+ϖ+b+r2aϖk-br,β02=kaϖ+a-b+r-br+kaϖk-br2. The dynamics of the system on the center manifold is then described by (27)u˙1u˙2=u20+1ΨGu2,u3aGu2,u3,where(28)Ψ=2b24ϖ-br2ϖ+14ϖ+kϖ+r24ϖ+b2+r2,Gu2,u3=2ϖu12u2r5abr-kr2+k2-ϖu13r5br-kr2+k2-ϖu1u22r5a2br-kr2+k2+6ϖu1u22r4abr-kr2+k2+6u1u22r4ar2+k2ϖ-4ϖu12u2r4br-kr2+k2-2ϖu23r4a2br-kr2+k2-4u12u2r4r2+k2ϖ-2u23r4a2r2+k2ϖ+2bu23r3a2r2+k2ϖ+4ϖu1u22r3br-kr2+k2+4u1u22r3r2+k2ϖ+4kϖu12u2r3abr-kr2+k2+4bu12u2r3r2+k2ϖ-2kϖu13r3br-kr2+k2-4ϖu23r3abr-kr2+k2-2kϖu1u22r3a2br-kr2+k2-4u23r3ar2+k2ϖ-6bu1u22r3ar2+k2ϖ+4bu23r2ar2+k2ϖ+2ku23r2a2r2+k2ϖ+6kϖu1u22r2abr-kr2+k2+4ku12u2r2r2+k2ϖ-4kϖu12u2r2br-kr2+k2-2kϖu23r2a2br-kr2+k2-4bu1u22r2r2+k2ϖ-6ku1u22r2ar2+k2ϖ+4kϖu23rabr-kr2+k2+20ku23rar2+k2ϖ+2bku23ra2r2+k2ϖ+2ϖu1u22rbr-k2+2k2ϖu12u2rabr-kr2+k2+4bku12u2rr2+k2ϖ-2ϖu23rabr-k2-k2ϖu13rbr-kr2+k2-4kϖu1u22rbr-kr2+k2-k2ϖu1u22ra2br-kr2+k2-20ku1u22rr2+k2b2-2rb+r2+4k-6bku1u22rar2+k2ϖ+4k2u23a2r2+k2ϖ+4bku1u22r2+k2ϖ+8k2u12u2r2+k2ϖ+-a2r2u13-a2ku13+2ar2u2u12+2aku2u12+2a2u22u1-r2u22u1-ku22u1-2au23a2br-k-u1au1-u22a2-2-r2u23-ku23+ar2u1u22+aku1u22abr-k2-4bku23ar2+k2ϖ-12k2u1u22ar2+k2ϖ.Conditional normal form at bifurcation values of parameters can be obtained by following Kuznetsov approach which implies that the simplified system which represents the conditional normal form at bifurcation value is the following system:(29)v˙1=v2,v˙2=ar1-brv13+ra-3b2+6br-3r2-12-3b2+6br-3r2-121b3r+b2-2r2-1+br1r2+6-1r2-4v12v2The next step in this analysis is to obtain the universal unfolding of the BT bifurcation. In order to put the system in an appropriate form, the bifurcation parameters should be perturbed around bifurcation value, such thatc=-a+Δ1 and d=-a+Δ2. The following system is attained(30)x˙1=a-x1+x2,x˙2=-a+Δ1x2--a+Δ2x2-x1x3,x˙3=x4+x1x2-bx3,x˙4=-rx4+kx3,Δ˙1=0,Δ˙2=0. Now, applying transformation (22) to get(31)u˙1u˙2u˙3u˙4=ϝ~1ϝ~2ϝ~3ϝ~4,(32)Δ˙1=0,Δ˙2=0,where (33)ϝ~1=u4u2ϖ2ak-u3u2ϖ2ak-bu3u22ak-bu4u22ak+ru3u22ak+ru4u22ak-Δ1u2a+u1u3ϖ2k-u1u4ϖ2k+bu1u32k+bu1u42k-ru1u32k-ru1u42k+Δ1u1-Δ2u1+u2,ϝ~2=au1u3ϖ2k-au1u4ϖ2k+abu1u32k+abu1u42k-aru1u32k-aru1u42k+aΔ1u1-aΔ2u1-u2u3ϖ2k+u2u4ϖ2k-bu2u32k-bu2u42k+ru2u32k+ru2u42k-Δ1u2,ϝ~3=ku1u2aϖ-b2u34ϖ-b2u44ϖ+bru32ϖ+bru42ϖ-14u3ϖ+14u4ϖ-ku12ϖ-ku3ϖ-ku4ϖ-r2u34ϖ-r2u44ϖ-bu32-ru32,ϝ~4=-ku1u2aϖ+b2u34ϖ+b2u44ϖ-bru32ϖ-bru42ϖ+ku12ϖ-14u3ϖ+ku3ϖ+r2u34ϖ+14u4ϖ+ku4ϖ+r2u44ϖ-bu42-ru42.The center manifold is assumed as follows:(34)u3=∑i,j,k,l=0i+j+k+l=2α~ijklu1iu2jΔ1kΔ2l,u4=∑i,j,k,l=0i+j+k+l=2β~ijklu1iu2jΔ1kΔ2l,where parameters α~ijkl and β~ijkl can be attained through the same way employed in the first part of this subsection. The following reduced system on center manifold is obtained after some simplifications (35)u˙1u˙2=Au1u2+G~1u1,u2,Δ1,Δ2G~2u1,u2,Δ1,Δ2,A=0100,where (36)
G
~
1
u
1
,
u
2
,
Δ
1
,
Δ
2
=
u
1
2
u
2
b
3
a
b
2
+
ϖ
∗
-
2
r
b
+
r
r
+
ϖ
∗
+
16
2
+
ϖ
∗
u
1
2
u
2
b
2
a
b
2
+
ϖ
∗
-
2
r
b
+
r
r
+
ϖ
∗
+
16
2
-
3
r
u
1
2
u
2
b
2
a
b
2
+
ϖ
∗
-
2
r
b
+
r
r
+
ϖ
∗
+
16
2
+
u
1
3
b
2
b
r
-
4
+
b
2
-
2
r
b
+
r
2
+
16
u
1
2
u
2
b
2
a
ϖ
∗
b
r
-
4
+
4
ϖ
∗
u
1
2
u
2
b
b
2
+
ϖ
∗
-
2
r
b
+
r
r
+
ϖ
∗
+
16
2
+
3
r
2
u
1
2
u
2
b
a
b
2
+
ϖ
∗
-
2
r
b
+
r
r
+
ϖ
∗
+
16
2
+
16
u
1
2
u
2
b
a
b
2
+
ϖ
∗
-
2
r
b
+
r
r
+
ϖ
∗
+
16
2
-
u
1
2
u
2
b
2
a
b
r
-
4
-
b
2
-
2
r
b
+
r
2
+
16
u
1
3
b
2
ϖ
∗
b
r
-
4
+
r
b
2
-
2
r
b
+
r
2
+
16
u
1
2
u
2
2
a
ϖ
∗
b
r
-
4
+
r
u
1
2
u
2
2
a
b
r
-
4
-
r
u
1
3
2
b
r
-
4
-
r
b
2
-
2
r
b
+
r
2
+
16
u
1
3
2
ϖ
∗
b
r
-
4
-
4
ϖ
∗
r
u
1
2
u
2
b
2
+
ϖ
∗
-
2
r
b
+
r
r
+
ϖ
∗
+
16
2
-
16
r
u
1
2
u
2
a
b
2
+
ϖ
∗
-
2
r
b
+
r
r
+
ϖ
∗
+
16
2
,
G
~
2
u
1
,
u
2
,
Δ
1
,
Δ
2
=
-
a
b
u
1
3
b
2
-
2
b
r
+
r
2
+
16
2
ϖ
∗
b
r
-
4
-
a
r
u
1
3
b
2
-
2
b
r
+
r
2
+
16
2
ϖ
∗
b
r
-
4
+
4
a
b
u
1
2
u
2
ϖ
∗
b
2
+
b
ϖ
∗
-
2
r
+
r
ϖ
∗
+
r
+
16
2
+
a
b
u
1
3
2
b
r
-
4
-
a
r
u
1
3
2
b
r
-
4
+
b
u
1
2
u
2
b
2
-
2
b
r
+
r
2
+
16
2
ϖ
∗
b
r
-
4
+
r
u
1
2
u
2
b
2
-
2
b
r
+
r
2
+
16
2
ϖ
∗
b
r
-
4
+
b
2
u
1
2
u
2
ϖ
∗
b
2
+
b
ϖ
∗
-
2
r
+
r
ϖ
∗
+
r
+
16
2
+
16
b
u
1
2
u
2
b
2
+
b
ϖ
∗
-
2
r
+
r
ϖ
∗
+
r
+
16
2
-
16
r
u
1
2
u
2
b
2
+
b
ϖ
∗
-
2
r
+
r
ϖ
∗
+
r
+
16
2
+
b
3
u
1
2
u
2
b
2
+
b
ϖ
∗
-
2
r
+
r
ϖ
∗
+
r
+
16
2
-
b
u
1
2
u
2
2
b
r
-
4
+
r
u
1
2
u
2
2
b
r
-
4
,such that ϖ∗=(b-r)2+16.Now, it necessary to define the next vector valued functions(37)Cs,v,w=∑i,j,k=13∂3G~u1,u2,Δ1,Δ2∂ui∂uj∂uku1=u2=0Δ1=Δ2=0sivjwk,A1s,Λ=∑i=12∑j=12∂2G~u1,u2,Δ1,Δ2∂ui∂Δju2=u3=0Δ1=Δ2=0siΛj,A2s,Λ,Ω=∑i=12∑j,k=12∂2G~u2,u3,Δ1,Δ2∂ui∂Δj∂Δku2=u3=0Δ1=Δ2=0siΛjΩk,where (38)G~u2,u3,Δ1,Δ2=G~1u1,u2,Δ1,Δ2G~2u1,u2,Δ1,Δ2.The exact expression for these vector valued functions is omitted for brevity. Now, following the procedure described in [39], we assume that the following relation holds between unfolding parameters η=(η1, η2)T and perturbation values Δ=(Δ1,Δ2)T(39)Δi=ρiη1+υiη2+γiη12+δiη1η2+ζiη22+Oη3,i=1,2or equivalently(40)ηi=ρ~iΔ1+υ~iΔ2+γ~iΔ12+δ~iΔ1Δ2+ζ~iΔ22+OΔ3,i=1,2.The two real linearly independent generalized eigenvectorsq1 and q2 that satisfy Aq1=0 and Aq2=q1 can be found as (41)q1=10,q2=01.Also, for the transposed matrix AT, there exist eigenvectors p1 and p2 that have the properties (42)ATp2=0,ATp1=p2,p1Tq1=p2Tq2=1,p1Tq2=p2Tq1=0, and they are computed as (43)p1=10,p2=01.After some calculations, the linear part coefficients of (39) can be evaluated as (44)ρ1=1a,ρ2=0,υ1=-1,υ2=-1. Then, the transversality condition is verified by evaluating (45)det⁡1a-10-1=-1a≠0,which implies that transversality condition is fulfilled. The other coefficients of (39) can be evaluated by solving a set of systems of linear equations described in [39] which results in (46)γ1=γ2=0,δ1=δ2=0,ζ1=ζ2=0.By substituting from (39) into (40) and comparing the coefficients of η1iη2j, i+j=1,2, on both sides of (40), the unfolding parameters η1(Δ1,Δ2) and η2(Δ1,Δ2) can be obtained in the following expression (47)η1Δ1,Δ2=aΔ1-aΔ2,η2Δ1,Δ2=-aΔ2.The system representing the universal unfolding of BT bifurcation is thus written as(48)u˙1=u2,u˙2=aΔ1-aΔ2u1-aΔ2u2+ar1-bru13+ra-3b2+6br-3r2-12-3b2+6br-3r2-121b3r+b2-2r2-1+br1r2+6-1r2-4u12u2.Let ϑ30=ar/1-br and ϑ21=ra-3b2+6br-3r2-12-3b2+6br-3r2-12/1b3r+b2-2r2-1+br1r2+6-1r2-4, then apply the following rescaling of coordinates and time to further simplify system (48) (49)u1=s1ϑ30ϑ21ξ1,u2=-s1s2ϑ303/2ϑ212ξ2,t=-ϑ21ϑ30τ. After simplification, we can obtain(50)ξ˙1=ξ2,ξ˙2=μ1ξ1+μ2ξ2+s1ξ13-ξ12ξ2,where derivative is considered with respect to τ.The bifurcation structure of (50) is summarized in the next theorem [39–41].Theorem 3.
The bifurcation diagram of (50) in μ1-μ2 plane involves the following bifurcation curves:
(i)PB={(μ1,μ2):μ1=0,μ2≠0},
(ii)HB1={(μ1,μ2):μ2=0,μ1<0},
(iii)HB2={(μ1,μ2):μ2=μ1+O(μ12),μ1>0},
(iv)HLB={(μ1,μ2):μ2=4/5μ1+O(μ13/2),μ1>0},
(v)SNB={(μ1,μ2):η2=0.752μ1+O(μ13/2),μ1>0}, where μ1=9(aΔ1+a2(Δ12+Δ1Δ2)) and μ2=3(a(Δ1+Δ2)+a2(2Δ12+3Δ1Δ2+Δ22)) where PB represents a pitchfork bifurcation of system (50), HB1 denotes Hopf bifurcation of the origin equilibrium point, HB2 refers to the Hopf bifurcation of the nontrivial equilibrium points, HLB corresponds to Homoclinic connection to the origin, and finally SNB denotes the saddle-node bifurcation of symmetric periodic orbits.The phase portraits of system (1) are obtained via numerical simulations and illustrate examples of different dynamics that can be exhibited by the hyperchaotic system. In particular, Figure 1 verifies the occurrences of Hopf and pitchfork bifurcations of system (2) whereas Figure 2 shows 2D and 3D portraits of hyperchaotic attractor of the system.Figure 1
Periodic orbit of system (1) is created around nontrivial equilibrium points via Hopf bifurcation for a=10; b=7.3; c=-9.8; d=-9.9; k=-9.5; r=0.1 in (a,b), stable equilibrium point at the origin for a = 10; b = 7.3; c = -9.99; d = -9.9; k = -9.5; r = 0.1 is shown in (c), the origin loses its stability, and two stable nontrivial equilibrium points are created via pitchfork bifurcation and shown in (d,e) for a = 7; b = 7.3; c = -7; d = -6.9; k = -9.5; r = 0.1.
(a) (b) (c) (d) (e)Figure 2
Hyperchaotic attractor of system (1) is illustrated via phase portraits on (a) 3D x-y-z space, (b) x-y-w space, (c) 2D x-y plane, (d) 2D x-z plane, (e) 2D y-z plane, and (f) 2D y-w plane at the values of parameters: a=21.7,b=7.3,c=6.6,d=-2,r=0.1, and k=9.5.
(a) (b) (c) (d) (e) (f)Bifurcation diagrams are used to provide useful illustrations for variations of a particular system dynamics with changes in its parameters. Figure3 shows examples of bifurcation diagrams of system (1) with respect to different parameters in the system. Figure 4 depicts the associated Lyapunov exponent plots to the cases presented in Figure 3.Figure 3
Bifurcation diagrams of hyperchaotic system (1) obtained for the following values of parameters: (a) 0≤a≤30,b=7.3,c=6.6,d=-2,r=0.1, and k=9.5. (b) a=21.7,0≤b≤12,c=6.6,d=-2,r=0.1, and k=9.5. (c) a=21.7,b=7.3,0≤c=10,d=-2,r=0.1, and -12≤k≤0.
(a) (b) (c) (d)Figure 4
Lyapunov exponent spectrum corresponding to Figure3.
(a) (b) (c) (d)


## 4. First Integrals of the System

It is well known that, if the analyzed system has chaotic dynamics, the behaviors of the system are irregular behaviors. So, it is impossible to predict the behavior of the system in a long time scale. Therefore, it is important to ask whether there is any set of parameter values in which system dynamics are regular and that the studied system is integrable. In this section, we provide the question of global analytic integrability investigated for hyperchaotic system (1). In particular, we provide conditions of analytic integrability in hyperchaotic system (1) and obtain the possible first integrals of it.Theorem 4.
System (1) has the following first integrals:(1)
I1(x,y,z,w,t)=w2-k(y2+z2) when c=b=d=r=0,(2)
I2(x,y,z,w,t)=e2btbw2-kby2+z2+ckx2 when a=d=r=b,(3)
I3(x,y,z,w,t)=e2bt2c3w2+cbx2-b2y2+z2+c2-2bz+2w+x2 when a=c=d=r=b and k=b2,(4)
I4(x,y,z,w,t)=e2atx2(r-2a)+2az(2a-r)-2aw when k=(2a-b)(2a-r),(5)
I5(x,y,z,w,t)=e2bt2cb(bz-w)-b2y2+z2+w2 when c=d=r=b and k=b2,
wherec2 and c3 are arbitrary constants.Proof.
The functionI(x(t),y(t),z(t),w(t),t) is said to be the first integral of system (1) if it satisfies the condition(51)dIdx=It+x˙Ix+y˙Iy+z˙Iz+w˙Iw=0,Substituting (1) into (51), we obtain(52)It+ay-xIx+cx-dy-xzIy+-bz+xy+wIz+-rw+kzIw=0.Assume that the solution of (52) takes the form(53)Ix,y,z,w,t=Rx,y,z,west.where s is a constant. Substituting (53) into (52), we obtain(54)sR+ay-xRx+cx-dy-xzRy+-bz+xy+wRz+-rw+kzRw=0,Let the solution of (54) take the form(55)Rx,y,z,w=Ax+By+Tz+Qw.Substituting (55) into (54), we obtain(56)sAx+By+Tz+Qw+ay-xA′x+cx-dy-xzB′y+-bz+xy+wT′z+-rw+kzQ′w=0.By derivative (56) twice with respect to w, we obtain(57)kzQ′′′w-rwQ′′′w+s-2rQ′′w=0.Equating the coefficient of z with zero and solving the obtained equation, we obtain(58)Qw=c1+c2w+c3w2,where c1 are arbitrary constant. Substituting (58) into (56), we obtain(59)sTz+xyT′z+ay-xA′x+sAx-bzT′z+B′ycx-dy-xz+sBy+c1s+w-c2r+c2s+2c3kz+T′z+c2kz+w2c3s-2c3r=0.Equating the coefficient of w with zero and solving the obtained equation, we obtain(60)Tz=c2r-sz-c3kz2+c4,where c4 are arbitrary constant. Substituting (60) into (59), we obtain(61)c4s+ay-xA′x+sAx+By+zc2b-ss-r+k+c3kz22b-s+B′ycx-dy+xz-B′y-2c3ky+c1s+c2xyr-s+c3w2s-2r=0Equating the coefficient of xz with zero and solving the obtained equation, we obtain(62)By=c5-c3ky2,where c5 are arbitrary constant. Substituting (62) into (61), we obtain(63)yaA′x+xc2r-s-2cc3k-axA′x+sAx+c2zb-ss-r+k+c3kz22b-s+c1s+c3ky22d-s-2c3rw2+c3sw2+c4s+c5s=0.Equating the coefficient of y with zero and solving the obtained equation, we obtain(64)Ax=12ax22cc3k+c2s-r+c6,where c6 are arbitrary constant. Substituting (64) into (63), we obtain(65)12ax2s-2a2cc3k+c2s-r+c2zb-ss-r+k+c3kz22b-s+sc1+c4+c5+c6+c3ky22d-s+c3w2s-2rEquation (65) satisfies when (1)
s=0,b=c=d=r=0, and c2=0,(2)
s=2b,a=d=r=b,c6=-c1-c4-c5, and c2=0,(3)
s=2b,a=d=r=b,c6=-c1-c4-c5, and k=b2,(4)
s=2a,k=(2a-b)(2a-r),c6=-c1-c4-c5, and c3=0,(5)
s=2b,d=r=b,c6=-c1-c4-c5,k=b2, and c3=-c2/2bc.
So, when we substitute (58), (60), (62), and (64) into (53), we obtain the following first integrals of system (1), respectively: (1)
I1(x,y,z,w,t)=w2-k(y2+z2),(2)
I2(x,y,z,w,t)=e2btbw2-kby2+z2+ckx2,(3)
I3(x,y,z,w,t)=e2bt2c3w2+cbx2-b2y2+z2+c2-2bz+2w+x2,(4)
I4(x,y,z,w,t)=e2atx2(r-2a)+2az(2a-r)-2aw,(5)
I5(x,y,z,w,t)=e2bt2cb(bz-w)-b2y2+z2+w2.
Hence, the proof is completed.


### 4.1. First Integrals Containing the Origin Equilibrium Point

In this subsection, the invariant surfaces in phase space of hyperchaotic system (1) which contain the origin equilibrium point are attained. If we substitute the equilibrium points into the first integrals of the system (1), we obtain(66)I1=w2-ky2+z2=0,(67)I2=y2+z2-cbx2-1kw2=0,(68)I3=2c3w2+cbx2-b2y2+z2+c2-2bz+2w+x2=0,(69)I4=x2r-2a+2az2a-r-2aw=0,(70)I5=2cbbz-w-b2y2+z2+w2=0.Hence, we will classify the solutions (66)-(70) as follows: (1)
For (66), when k>0, (66) becomes circular cone (see Figure 5).(2)
For (67), there are some cases:
(i)
whenc=0 and k>0, (67) becomes circular cone (see Figure 6(a)),(ii)
whenc/b>0, k>0, and w=β, (67) becomes hyperboloid of one sheet (see Figure 6(b)),(iii)
whenc/b<0, k>0, and w=β, (67) becomes ellipsoid (see Figure 6(c)),(iv)
whenc=b, k>0, and w=β, (67) becomes sphere (see Figure 6(d)),(v)
whenc/b>0, k>0, and w=β, (67) becomes hyperboloid of two sheets (see Figure 6(e)).(3)
Letc2=2c3; (68) becomes(71)y2+1b+z2-bc+1b2x2-1b2w+12=0.which is classified as follows:
(i)
whenbc=-1, (71) becomes circular cone (see Figure 7(a)),(ii)
whenbc>-1 and w=β, (71) becomes hyperboloid of one sheet (see Figure 7(b)),(iii)
whenbc<-1 and w=β, (71) becomes ellipsoid (see Figure 7(c)),(iv)
whenc=-b-1/b and w=β, (71) becomes sphere (see Figure 7(d)).(4)
Equation (69) represents parabola (see Figure 8).(5)
Equation (70)can be rewritten in the form(72)y2+z-c2-w-bc2b2=0,which represents circular cone (see Figure9).Figure 5
Plot the equationsI1=0 for k=0.5.Figure 6
Plot of the equationsI2=0 at different values of the parameters: (a) c=0 and k=0.5; (b) c/b=0.2 and β=k, k>0; (c) c/b=-0.5 and β=k, k>0; (d) c/b=-1 and β=k, k>0; (e) c/b=1 and β=-k, k<0.
(a) (b) (c) (d) (e)Figure 7
Plot the equationsI3=0 at different values of the parameters: (a) c=-1 and b=1; (b) c=1,b=2, and β=2; (c) c=2,b=-1, and β=2; (d) c=-2.5,b=2, and β=2.
(a) (b) (c) (d)Figure 8
Plot of the equationI4=0 for r=a=1.Figure 9
Plot of the equationI5=0 for b=c=1.β is a constant.


## 4.1. First Integrals Containing the Origin Equilibrium Point

In this subsection, the invariant surfaces in phase space of hyperchaotic system (1) which contain the origin equilibrium point are attained. If we substitute the equilibrium points into the first integrals of the system (1), we obtain(66)I1=w2-ky2+z2=0,(67)I2=y2+z2-cbx2-1kw2=0,(68)I3=2c3w2+cbx2-b2y2+z2+c2-2bz+2w+x2=0,(69)I4=x2r-2a+2az2a-r-2aw=0,(70)I5=2cbbz-w-b2y2+z2+w2=0.Hence, we will classify the solutions (66)-(70) as follows: (1)
For (66), when k>0, (66) becomes circular cone (see Figure 5).(2)
For (67), there are some cases:
(i)
whenc=0 and k>0, (67) becomes circular cone (see Figure 6(a)),(ii)
whenc/b>0, k>0, and w=β, (67) becomes hyperboloid of one sheet (see Figure 6(b)),(iii)
whenc/b<0, k>0, and w=β, (67) becomes ellipsoid (see Figure 6(c)),(iv)
whenc=b, k>0, and w=β, (67) becomes sphere (see Figure 6(d)),(v)
whenc/b>0, k>0, and w=β, (67) becomes hyperboloid of two sheets (see Figure 6(e)).(3)
Letc2=2c3; (68) becomes(71)y2+1b+z2-bc+1b2x2-1b2w+12=0.which is classified as follows:
(i)
whenbc=-1, (71) becomes circular cone (see Figure 7(a)),(ii)
whenbc>-1 and w=β, (71) becomes hyperboloid of one sheet (see Figure 7(b)),(iii)
whenbc<-1 and w=β, (71) becomes ellipsoid (see Figure 7(c)),(iv)
whenc=-b-1/b and w=β, (71) becomes sphere (see Figure 7(d)).(4)
Equation (69) represents parabola (see Figure 8).(5)
Equation (70)can be rewritten in the form(72)y2+z-c2-w-bc2b2=0,which represents circular cone (see Figure9).Figure 5
Plot the equationsI1=0 for k=0.5.Figure 6
Plot of the equationsI2=0 at different values of the parameters: (a) c=0 and k=0.5; (b) c/b=0.2 and β=k, k>0; (c) c/b=-0.5 and β=k, k>0; (d) c/b=-1 and β=k, k>0; (e) c/b=1 and β=-k, k<0.
(a) (b) (c) (d) (e)Figure 7
Plot the equationsI3=0 at different values of the parameters: (a) c=-1 and b=1; (b) c=1,b=2, and β=2; (c) c=2,b=-1, and β=2; (d) c=-2.5,b=2, and β=2.
(a) (b) (c) (d)Figure 8
Plot of the equationI4=0 for r=a=1.Figure 9
Plot of the equationI5=0 for b=c=1.β is a constant.


## 5. Circuit Implementation

The circuit implementation of hyperchaotic system (1) can be designed using operational amplifier, analog multipliers, resistors, and capacitors [17, 18]. This hyperchaotic nonlinear circuit is equivalent to system (1) and enables the useful experimental study of complex dynamics exhibited by system (1) and verify the results of theoretical analysis carried out upon it. The Multisim package is used to carry out circuit simulations and depict the circuit outputs on the oscilloscope. The nonlinear circuit implementation of hyperchaotic system can be used as a central in a variety of advanced applications, e.g., chaos based physical-layer secure communications and physical random number generators which outperform pseudo random number generators. Figure 10 shows schematics of proposed the circuit realization of hyperchaotic system (1) for the values of parameters given by a=21.7,b=7.3,c=6.6,d=-2,r=0.1, and k=9.5.Figure 10
(a) Circuit implementation of system (1), (b–e) circuit simulations outputs representing x-y,y-z,x-z, and x-w phase portraits of the model at a=21.7;b=7.3;c=6.6;d=-2;r=0.1;k=9.5.
(a) (b) (c) (d) (e)Regarding digital realization of chaotic circuits, it is important to note that dynamics degradation of digital chaos occurs due to impossibility of the chaotic map to reach an ideal chaotic state under the limitations of digital domain with finite-precision. Various methods are proposed to counteract dynamics degradation; see, for example, [42] and references therein.


## 6. Conclusion

In this work, dynamical and bifurcation analysis are applied to explore local and global characteristics of 4D Lorenz-type hyperchaotic system. Through investigation of codimension two local bifurcation of the hyperchaotic system, it is found that the 4D hyperchaotic system undergoes various types of bifurcations including BT bifurcation, Hopf bifurcation, pitchfork bifurcation, and homoclinic bifurcation. Moreover, the paper successfully acquires, for first time, the closed forms of first integrals and regular invariant surfaces in solution space, which pass through the origin equilibrium points, and relates these forms to the parameters of the system. Thus, it predicts some analytical solutions of the system and estimates the associated parameters values. Finally, a circuit realization for 4D hyperchaotic system is proposed as an engineering application and numerical examples are presented to verify theoretical results. Future work can include analytical and numerical investigations of hidden attractors, if exist, in hyperchaotic 4D Lorenz system. Also, other possible first integrals and invariant structures of the state space of the system may be explored in future work.



---

*Source: 6581586-2019-02-06.xml*