A New Family of Chaotic Systems with Different Closed Curve Equilibrium

Chaotic systems with hidden attractors, infinite number of equilibrium points and different closed curve equilibrium have received much attention in the past six years. In this work, we introduce a new family of chaotic systems with different closed curve equilibrium. Using the methods of equilibrium points, phase portraits, maximal Lyapunov exponents, Kaplan–Yorke dimension, and eigenvalues, we analyze the dynamical properties of the proposed systems and extend the general knowledge of such systems.


Introduction
Over the past fifty years, investigating chaotic systems has attracted many researchers' attention.Generating chaotic attractors may help one to understand the dynamics of the real system.Recently, new chaotic systems with hidden attractors have been introduced [1,2].An attractor is called self-excited if its basin of attraction intersects any open neighborhood of a stationary state (an equilibrium), otherwise it is called a hidden attractor.It is generally accepted that there are three main families of hidden attractors [1,3,4]: hidden attractors with stable equilibria, hidden attractors without equilibrium and hidden attractors with infinite number of equilibrium points.Some developments include chaotic systems without any equilibrium points [5][6][7][8][9][10], with a single stable equilibrium [11,12], and with uncountable equilibria [13][14][15].However, there is little discussion on systems with different families of hidden attractors [4,16].Undoubtedly, many undiscovered systems surely exist.In this study, the authors present a family of chaotic systems with different closed curve equilibriums, being the promotion of published papers, which may extend the general knowledge of such systems.
Gotthans and Petr žela [14] have introduced a general three-dimensional model generating hidden attractors with infinite equilibria.The model is given as: where x, y and z are state variables, while g 1 (x, y, z) and g 2 (x, y, z) are two nonlinear functions.
Motivated by the works mentioned above, in this paper, we introduce and investigate a new family of chaotic systems with different closed curve equilibriums.Replacing g 1 (x, y, z) with ay + by 2 + xz, and g 2 (x, y, z) with |x| k + |y| k − 1 in (1), i.e., the proposed novel system is obtained as where a and b are both positive parameters, while k is a positive integer.In fact, the new model (3) contains several known systems as special cases.

(Example A)
For the case of k = 2, system (3) corresponds to the following system: Actually, model (4) was studied by Gotthans and Petr žela in [15] and has shown that there exists some chaotic system with a circle of equilibrium points as (x, y, 0) :

(Example B)
For the case of k = 4, system (3) corresponds to the following system: Model ( 5) is exactly the system studied by Pham et al. in [4] and have been shown that there exists a chaotic system with a rounded square equilibrium as (x, y, 0) : Above all, the aim of this study is to introduce a new family of chaotic systems generating hidden attractors with differently shaped equilibrium points, containing systems ( 4) and ( 5) as two particular cases, being the generalization of the existing results, and to broaden the views on the hidden attractors.
For any integer k, by solving the following system we obtain the equilibrium points of system (3).It is easy to see that the proposed system (3) has a family of different curve equilibriums of (x, y, 0) :

Chaotic Behavior of the Proposed System
Consider a general n-dimensional autonomous vector field An equilibrium solution of ( 7) is a point x * ∈ R n such that f (x * ) = 0, i.e., a solution which does not change in time.Other terms substituted for the term equilibrium solution are fixed point or critical point.In order to determine the stability of x * , we must understand the nature of solution near x * .We do this by linearzing (7) at x * , that is, by studying the following linear system where D f = [∂ f i /∂x j ] is the Jacobian matrix of first partial derivatives of the function f = ( f 1 (x 1 , . . ., x n ), f 2 (x 1 , . . ., x n ), . . ., f n (x 1 , . . ., x n )) T (T denotes transpose), and x = x * + ξ.
An equilibrium is called hyperbolic if there are no eigenvalues of D f (x) on the imaginary axis.The following well-known theorem for hyperbolic equilibrium is very crucial in this paper.Theorem 1. [17] Suppose that all of the eigenvalues of D f (x * ) have negative real parts.Then, the equilibrium solution x = x * of the nonlinear vector field is asymptotically stable.
However, it is worth mentioning that the situation is very complicated for the nonhyperbolic equilibrium.Let us recall some definitions that are well known in the literature.Definition 1. [18] An attractor is a minimal closed invariant set that attracts nearby trajectories lying in the domain of stability (or basin of attraction) onto it.

Definition 2. [18]
A chaotic attractor is an attractor that exhibits sensitivity to initial conditions.Definition 3. [18] The spectrum of Lyapunov exponents are quantities that characterize the rate of separation of infinitesimally close trajectories.
Chaos is a phenomenon that is not easily identified or classified.There is also no universally accepted definitions for chaos, but there are some characteristics displayed by the solutions of chaotic systems: sensitivity to initial conditions; long-term nonperiodic behavior and fractal structure.Generally, a system is chaotic if at least one of the Lyapunov exponents is positive.This implies that two trajectories that start initially close to each other on the strange attractor will diverge as time increases.
It is worthwhile to note that system (3) with a family of different curve equilibrium points displays chaotic attractors (see Figure 2) for a = 5, b = 3, and the initial conditions (x(0), y(0), z(0)) = (0, 0, 0).The black closed contour in the middle of each of Figure 2a-f is the plot of the closed curve equilibrium for the different values of k, respectively.Chaotic behavior of the system is generally confirmed by the Lyapunov exponents calculated by the Wolf method and the Kaplan-Yorke dimension D KY (see Table 1 and Figure 3).In simulation, the fourth-order Runge-Kutta integration method has been used with a final time of 10, 000 and the time step equals 0.01.The results are very interesting that, with the increasing values of k, the corresponding largest Lyapunov exponents have a trend of decreasing.Dynamical motion in the close neighborhood of the equilibrium is determined by the eigenvalues and associated eigenspace established along this structure [14,15].The partial derivatives of the state variables are required in estimating the Jacobian matrix.The derivation of |•| can be obtained in several ways.Here, sgn(•) is the signum function.It is known that the signum function is not differentiable at 0.
In the case of the system (3), a state-dependent linearization matrix is easily established as The equilibrium points lie in the plane z = 0, so we have The local behavior along the equilibrium is determined by the eigenvalues, i.e., the roots of the characteristic equation The eigenvalues are given by ( Particularly, for k = 1, we have (See Table 2) ( Particularly, for k = 2, we have (See Table 3) ( In Tables 2 and 3, I, I I, I I I, IV represent the four quadrants, respectively.A pair of purely imaginary eigenvalues represent an unstable center equilibrium.Then, there are concentric periodic orbits around the equilibrium curves.In addition, the pair of purely real eigenvalues implies an unstable saddle.Furthermore, three-dimensional Bogdanov-Takens equilibria (λ 1 = λ 2 = λ 3 = 0) are also presented for the case of k = 1, 2.

Conclusions
The main goal of the current study is to propose a novel and more general family of chaotic systems with different curve equilibriums located on a quadrangle, a circle or other general cases.Furthermore, the chaotic systems here are the extension of existing literature.Dynamics of the new system have been investigated via different tools such as Phase portrait, Kaplan-Yorke dimension, Maximal Lyapunov exponents, and Eigenvalues.The results in this study provide a new and in-depth understanding of such systems.Physical realization of chaotic models have been used in different

6 Figure 2 .
Figure 2. 3D perspective view and plane projections of systems for k from 1 to 6.

Figure 3 .
Figure 3.The graph of the three Lyapunov exponents of system (3) for the value of k from 1 to 10.

Table 2 .
Eigenvalues of the system (3) for k = 1.< y < −0.2324 a pair of purely imaginary IV −0.2324 < y < 0 one positve real, one negative real x