# The Dynamics of a Cubic Nonlinear System with No Equilibrium Point

Authors: J. O. Maaita; Ch. K. Volos; I. M. Kyprianidis; I. N. Stouboulos
Journal: Journal of Nonlinear Dynamics
 (2015)

---

## Abstract

We study the dynamics of a three-dimensional nonlinear system with cubic nonlinearity and no equilibrium points with the use of Poincaré maps, Lyapunov Exponents, and bifurcations diagrams. The system has rich dynamics: chaotic behavior, regular orbits, and 3-tori periodicity. Finally, the proposed system is also reported to verify electronic circuit modeling feasibility.


---

## Body

## 1. Introduction

A lot of work has been done in the field of dynamical system and many systems (Lorenz, Chua, Duffing, Van der pol, Sprott, and many others) have been exhaustively studied. The dynamics of such systems are well known and their properties are used in mechanical and electrical applications and experiments [1–4].In the last two decades a new field of dynamical systems has been “discovered” and attracts the attention of scientists: dynamical systems with no equilibrium points or with conjugate equilibrium points.Equilibrium points are important because their stability determines the dynamics of the system [5–7]. In particular, a stable equilibrium point is a point for which the trajectories around it remain close for small perturbations. On the other hand, an unstable equilibrium point is a point for which the trajectories around it escape even for small perturbations and remove the system from its initial state.Equilibrium points are connected with criteria and theorems that determine the existence of chaotic behavior of a system (Melnikov function, Shilnikov chaos, etc.) [8]. The loss of equilibrium points means that the conventional Shilnikov criteria cannot be applied to prove the chaos in the flow.A dynamical system with no equilibrium points is categorized as chaotic system with hidden attraction because the loss of equilibrium points means that its basin of attraction does not intersect with small neighborhoods of any equilibrium points.Sprott (1994) was the first to introduce a simple flow with no equilibrium points [2]. Since then many researcher have introduced many systems with no equilibrium points or with conjugate equilibrium points [9–15].In this work we study a modified version of the initial Sprott model with a cubic nonlinearity and a constant parameterA.We made a numerical study of the system and used tools such as Poincaré maps, Lyapunov Characteristic Exponents, bifurcations diagrams [16, 17].The system has rich dynamics. In general it has a chaotic behavior but for certain initial values and different values for the parameterA the system may have regular orbits (quasiperiodic or periodic). It is important to note here that for some values of initial conditions we detected transient hyperchaotic behavior of the system.In Section2 we analyze the system and present the behavior of the system for different values of the constant parameter A and different values of the initial conditions. In Section 3 we present an electronic circuit that implements the above nonlinear system and finally conclude in Section 4 of the paper.


## 2. Analysis

We study a nonlinear system with cubic nonlinearity:(1)x˙=y,y˙=-x3-zy,z˙=y2-A,where A>0 is the parameter of the system. As it is obvious, since A>0, the system has no equilibrium point.We used many tools to analyze numerically the above system: Bifurcation diagrams, Poincaré maps (fory=0), and Lyapunov Characteristic Exponents. We simulated the system for many thousands of different initial conditions and different values of the parameter A.The numerical work was done with the help of Mathematica and the programming Languages C and True Basic by using the classical fourth-order Runge-Kutta method.As we see from the bifurcations diagrams (Figure1) the system is in general chaos. This is confirmed by the calculations of the Lyapunov Characteristic Exponents (LCEs) (Figure 2) where we see that for many different initial conditions and different values of the parameter A there are one positive LCE, one negative LCE, and one LCE that equals zero. This confirms the chaotic behavior of the system for these initial conditions.Figure 1
Bifurcation diagrams for different initial conditions taken fory=0 and y′>0.
(a)
x
0
=
0.1, y0=0.0, and z0=1.0
(b)
x
0
=
1.0, y0=10.0, and z0=2.0
(c)
x
0
=
30.1, y0=20.0, and z0=11.0
(d)
x
0
=
1.5, y0=2.0, and z0=1.5Figure 2
Lyapunov Characteristic Exponents for different initial conditions.
(a)
x
0
=
0.1, y0=0.0, z0=1.0, and A=0.6
(b)
x
0
=
1.0, y0=10.0, z0=2.0, and A=20
(c)
x
0
=
3.1, y0=2.0, z0=1.0, and A=5.0
(d)
x
0
=
0.1, y0=0.1, z0=0.2, and A=7.8Furthermore, from the study of the LCEs we detected that for many initial conditions and values ofA the system has a transient hyperchaotic behavior: that is, the system has two positive LCEs (Figure 3).Figure 3
Lyapunov Characteristic Exponents for different initial conditions where transient hyperchaotic behavior is detected.
(a)
x
0
=
0.01, y0=0.01, z0=4.51, and A=2.0
(b)
x
0
=
0.01, y0=1.01, z0=5.01, and A=4In what follows we will present three examples of the dynamics of the system, for different values of the parameterA.First, for the Case Where
A
=
5.16
. We see from the Poincaré section (Figure 4) that the system has both chaotic behavior and regular orbits (quasiperiodic and periodic). Also a paradigm of a 3-tori quasiperiodicity is detected.Figure 4
Poincaré section forA=5.16.For example, for the initial conditions (x0=-0.8, y0=0.0, and z0=1.0), the system has one positive LCE, one negative LCE, and one LCE that equals zero (Figure 5). This confirms the chaotic behavior. For these initial conditions the trajectories in various planes are given by Figure 6.Figure 5
Lyapunov Characteristic Exponents forA=5.16 and initial conditions: x0=-0.8, y0=0.0, and z0=1.0.Figure 6
Trajectories forA=5.16 and initial conditions: x0=-0.8, y0=0.0, and z0=1.0.
(a)
Trajectories on the 3D space
(b)
(
y
(
t
)
,
z
(
t
)
)
(c)
(
x
(
t
)
,
y
(
t
)
)
(d)
(
x
(
t
)
,
z
(
t
)
)For the initial conditions (x0=-0.8, y0=0.0, and z0=1.0) the system has three LCEs that equal zero (Figure 7). That is, the motion is regular and this can be seen in Figure 8.Figure 7
Lyapunov Characteristic Exponents forA=5.16 and initial conditions: x0=-1.6525, y0=0.0, and z0=0.0.Figure 8
Trajectories forA=5.16 and initial conditions: x0=-1.6525; y0=0.0; and z0=0.0
(a)
Trajectories on the 3D space
(b)
(
x
(
t
)
,
y
(
t
)
)
(c)
(
x
(
t
)
,
z
(
t
)
)
(d)
(
y
(
t
)
,
z
(
t
)
)For the Case Where
A
=
0.6
. From the Poincaré sections (Figure 9) we see that the regions of regular orbits are much smaller than the regions of the previous example.Figure 9
Poincaré section forA=0.6.For the initial conditions (x0=-1.2, y0=0.0, and z0=0.0) the system has three LCEs that are equal to zero (Figure 10) and this means that the system has regular orbits. This is also confirmed by the trajectories of the system (Figure 11).Figure 10
Lyapunov Characteristic Exponents forA=0.6 and initial conditions: x0=-1.2, y0=0.0, and z0=0.0.Figure 11
Trajectories forA=0.6 and initial conditions: x0=-1.2, y0=0.0, and z0=0.0.
(a)
Trajectories on the 3D space
(b)
(
y
(
t
)
,
z
(
t
)
)
(c)
(
x
(
t
)
,
z
(
t
)
)
(d)
(
x
(
t
)
,
y
(
t
)
)Also for the initial conditions (x0=2.0, y0=0.0, and z0=2.0) the system has one positive LCE, one negative LCE, and one LCE that equal zero (Figure 12) and the system has a chaotic behavior (Figure 13).Figure 12
Lyapunov Characteristic Exponents forA=0.6 and initial conditions: x0=2.0, y0=0.0, and z0=2.0.Figure 13
Trajectories forA=0.6 and initial conditions: x0=2.0, y0=0.0, and z0=2.0.
(a)
Trajectories on the 3D space
(b)
(
y
(
t
)
,
z
(
t
)
)
(c)
(
x
(
t
)
,
z
(
t
)
)
(d)
(
x
(
t
)
,
y
(
t
)
)For the Case Where
A
=
12.7
. A more complicated structure of the Poincaré section can be seen for the parameter A=12.7 (Figure 14). Beyond the regions of chaotic behavior and regular orbits a paradigm of a 3-tori periodicity also appears.Figure 14
Poincaré section forA=12.7.For the initial conditions (x0=-1.4, y0=0.0, and z0=0.0) the 3-tori periodicity paradigm appears. For this case, as it is shown from the figures of the LCEs and the trajectories (Figures 15 and 16) and as it is expected, the system has regular orbits.Figure 15
Lyapunov Characteristic Exponents forA=12.7 and initial conditions: x0=-1.4, y0=0.0, and z0=0.0.Figure 16
Trajectories forA=12.7 and initial conditions: x0=-1.4, y0=0.0, and z0=0.0.
(a)
Trajectories on the 3D space
(b)
(
y
(
t
)
,
z
(
t
)
)
(c)
(
x
(
t
)
,
z
(
t
)
)
(d)
(
x
(
t
)
,
y
(
t
)
)


## 3. Circuit Realization of the System

Circuital design of chaotic systems plays a crucial role in the field of nonlinear science not only for providing a simple experimental confirmation of phenomena related to nonlinear dynamics but also due to its applications in many engineering approaches, such as secure communication, signal processing, random bit generator, or path planning for autonomous mobile robot [18–23]. In addition, circuital implementation of chaotic systems also provides an effective approach for investigating new dynamics of such theoretical models [24, 25]. For example, time-series of chaotic signals or chaotic attractors can be observed at the oscilloscope easily or experimental bifurcation diagram can be obtained by varying the value of variable resistors [26, 27].In this work, an electronic circuit (Figure17) is introduced for implementing system (1), which is designed by using the operational amplifiers approach [27]. It has three integrators (U1, U5, and U7) and an inverting amplifier (U2), which are implemented with the TL084, as well as four signals multipliers (U3, U4, U6, and U8) by using the AD633. By applying Kirchhoff’s circuit laws, the corresponding circuital equations of designed master circuit can be written as (2)x˙=1R1C1y,y˙=-1R4C2100Vx3-1R5C210Vzy,z˙=1R6C310Vy2-1R7C3VA,where x,y, and z are the voltages in the outputs of the operational amplifiers U2, U5, and U7. Normalizing the differential equations of system (2) by using τ=t/RC we can see that this system is equivalent to the proposed dynamical system (1). The circuit components have been selected as R1=R2=R3=R7=10kΩ, R4=0.1kΩ, R5=R6=1kΩ, C1=C2=C3=10nF and VA adjust the value of the parameter A, while the power supplies of all active devices are ±15VDC.Figure 17
The schematic of the circuit that emulates the proposed dynamical system (1).The designed circuit is implemented in the electronic simulation package Cadence OrCAD and the obtained results are displayed in Figures18 and 19. In more details, these figures depict the simulation phase portraits produced by the OrCAD and the respective ones produced by the system’s arithmetic integration, for A=1 and A=5. The comparison of the chaotic attractors proves that the theoretical attractors are similar with the circuital ones. So, the designed circuit emulates very well the proposed system’s dynamic behavior.Figure 18
Simulation phase portraits produced by the OrCAD and the respective ones produced by the system’s arithmetic integration, forA=1.
(a)
x
-
y plane, OrCAD
(b)
x
-
y plane, numerically
(c)
x
-
z plane, OrCAD
(d)
x
-
z plane, numerically
(e)
y
-
z plane, OrCAD
(f)
y
-
z plane, numericallyFigure 19
Simulation phase portraits produced by the OrCAD and the respective ones produced by the system’s arithmetic integration, forA=5.
(a)
x
-
y plane, OrCAD
(b)
x
-
y plane, numerically
(c)
x
-
z plane, OrCAD
(d)
x
-
z plane, numerically
(e)
y
-
z plane, OrCAD
(f)
y
-
z plane, numerically


## 4. Conclusions

We study a nonlinear system with cubic nonlinearity and no equilibrium point through numerical simulations and confirm that the system has rich dynamics.Specifically, the system has, in general, chaotic behavior. A transient hyperchaotic (two positive LCEs) behavior is also detected.Also, for different initial conditions and different values of the parameterA the system may have regular orbits (periodic and quasiperiodic). Furthermore, examples of a 3-tori periodicity may appear for different values of the constant parameter A and different initial conditions.Finally, the designed nonlinear electronic circuit emulates very well the proposed system’s dynamic behavior.



---

*Source: 257923-2015-09-02.xml*