In recent years, we have seen that there is great interest in the subject of hidden and self-Excited attractors, the attractors in chaotic dynamical systems are classified into two categories: dynamical systems with self-Excited attractors and dynamical systems with hidden attractors, as it is said for the systems whose basin of attraction intersects with the neighborhood. The Lorenz system [1] is one of the chaotic systems that have self-Excited attractors, which was proposed in 1963, while the systems in which the basin of attraction does not intersect with the open neighborhood of equilibrium and which does not have fixed points or equilibrium or it has only one stable fixed point or nodal points, so it is called a dynamical system with hidden attractors. Jerk-like system [2] is one of the systems that have hidden attractors. There are many researchers who have provided a number of models in dynamical systems that have these types of attractors, Jafari and Sprott [3] managed to create a list of 17 chaotic models with nonlinear quadratic boundaries that have no equilibrium points, which were classified by Leonov and Kuznetsov [4] [5] as systems with hidden attractors, and Chan and Wang introduced an anarchic system with only one stable point equilibrium [6] through some modifications to Sprott system [7], and Wei [8] designed a fourdimensional system, which turned out to be a hyper-chaotic system without equilibrium points.

Models in dynamical systems have hidden attractors for many applications, especially in the field of engineering designs, because they do not allow for unforeseen accidents [7] such as potential disturbances in the structure of bridges and aircraft wings [9] and others. This paper was arranged as follows: system description, system analysis and finding fixed points, stability analysis of fixed points using stability tests: characteristic roots test, jury test [10] [11], study of the dynamical behavior of the system, finding phase space and bifurcation diagrams [12] [13] and Chaos diagnosing using chaos tests: the Lyapunov exponent test and Lyapunov dimension [14] [15], binary (0 - 1) test [16] [17]. The problem of chaos control attract attention of mathematician, researchers and engineers and there are many practical reasons for controlling chaos one of them chaos causes irregular behavior in nonlinear dynamical systems, therefor chaos control should be eliminated as much as possible or totally suppressed, adaptive control technique [18] [19] were used to control.

Two-dimensional discrete time dynamical system, suggested by Panahi [20] [21] and the define is as follows:-

F = { x t + 1 = x t + y t y t + 1 = 0.1 x t 2 + 0.1 + y t ( − x t − b y t ) + y t (1)

b = 2 (2)

where x t , y t represent variables system and b represent parameter system.

In this section analysis the system was done by means of its fixed points. Let us assume that:-

f 1 ( x t , y t ) = x t f 2 ( x t , y t ) = y t

To find the fixed points of system (1), suppose that:

y t = 0 (3)

0.1 x t 2 + 0.1 + y t ( − x t − b y t ) + y t = y t (4)

By solving Equations (3) and (4), we obtain the following points:-

p 0 = ( i , 0 ) , p 1 = ( − i , 0 )

Since the obtained fixed points are nodal points, we can say that the dynamical system (1) have hidden attractors.

The Jacobian matrix of system (1) obtained as:

J ( x t , y t ) = [ 1 1 0.2 x t − y t − x t − 2 b y t + 1 ] (5)

Proposition (1):- Let

L ( λ ) = a 2 λ 2 + a 1 λ + a 0 = 0 (6)

A characteristic Equation of (5) and | λ i | , i = 1 , 2 are values of the roots of Equation (6), then the following cases are true:

1) If | λ i | < 1 , then the fixed point of system (1) is locally asymptotically stable and is called the sink.

2) If | λ i | > 1 , then the fixed point of system (1) is unstable and is usually called source, but if at least one of the values of the roots of Equation (6) is greater than one, then the fixed point is called Saddle.

3) If | λ i | = 1 , then the fixed point of system (1) is called Non-Hyperbolic point, but if | λ i | ≠ 1 , then the fixed point is called Hyperbolic point.

Proposition (2):-

From the Equation (6), then the jury Table 1 is obtain as:

Such that

b k = | a 0 a n − k a n a k | ,   k = 0 , 1 ,     n = 2 , c k = | b 0 b n − 1 − k b n − 1 b k | ,   k = 0 ,   n = 2

Then the fixed point of system (1) is called stable if the satisfies following conditions:

L ( 1 ) > 0 , ( − 1 ) n L ( − 1 ) > 0 , | a 0 | < a n , | b 0 | > | b n − 1 | , | c 0 | > | c n − 2 | ,

where L is the characteristic equation. Otherwise, the fixed point is called unstable.

In this section the stability analysis of the fixed points of the system (1) is analyzed.

We will test the stability of the fixed points of system (1) by using the test of roots characteristic equation. To test the stability of the fixed point p 0 , we substitute the point p 0 = ( i , 0 ) in Equation (5), we obtain:-

J ( i , 0 ) = [ 1 1 0.2 i − i + 1 ]

Finding the determinant ( Det ( λ I − J ) = 0 ), we get:-

| [ λ 0 0 λ ] − [ 1 1 0.2 i − i + 1 ] | = 0

⇒ λ 2 + ( − 2 + i ) λ + ( 1 − 1.2 i ) = 0 (7)

By analyzing Equation (7), we get:-

λ 1 = 1 , λ 2 = 1

From the Proposition (1) we obtain the following values:-

| λ 1 | = 1 , | λ 2 | = 1

Then the fixed point p 0 = ( i , 0 ) is a non-hyperbolic point, and in the same way point p 1 = ( − i , 0 ) was tested, and its results are shown in Table 2, which turned out to be non-hyperbolic as well.

Note (1)

In this system, we did not use Jury test to study the stability of fixed points of the system, since Jury test depends on the coefficients Equation (7) which represent complex values.

After studying the theoretical side of the system, we will discuss the practical side of the system, and at the beginning we will use Newton Raphson’s numerical method, Newton Raphson’s numerical method was used to generate the best values for the system (1) with the least possible error, as the system parameter was fixed at the value b = 2, the following values were obtained: ( x , y ) = ( 1.8 , − 0.39 ) and an error of (0.0001). Time behavior of system (1) was studied by generating a time series of system states with time (t = 100), shown in Figure 1, which shows the unstable behavior of system states x t , y t with time. The phase space was found for the variables system (1), shown in Figure 2, which shows us the paths for variables system (1) with hidden attractors. The bifurcation diagrams

of the bifurcation parameter of the system (1) were found, the parameter b was fixed at the interval ranging from 1.975 to 2, shown in Figure 3, which shows the bifurcation of parameter b with variable x t , y t indicating chaotic behavior of the system (1).

To study the chaotic behavior of system (1) we will use Lyapunov exponent test, and after applying and calculating Lyapunov exponent of system (1) the following values were obtained:-

L 1 = − 0.250000 , L 2 = 1.853981

And since one of the values of Lyapunov exponent is a positive value, this indicates that system (1) is a chaotic system, and Figure 4 shows us the chaotic behavior of system (1).

To find Lyapunov dimension, we use the following formula:

p + L 1 | L 2 | = 1 + − 0.250000 | 1.853981 | = 1 − 0.134844 = 0.865156 = D L

To study and define the mess more broadly, binary test (0 - 1) was used to determine the chaos of system (1), and by using a mathematical program in MATLAB, the system binary test was applied as a time series of system states and (1000) were generated from the iterations with a fixed value Parameter b = 2 and ( x , y ) = ( 1.8 , − 0.39 ) , and q c ( t ) was calculated against p c ( t ) and shown

in Figure 5, which shows us a behavior similar to Brownian motion, and the mean square of the displacement M c ( t ) was found with time (t = 100) and shown in Figure 6, which shows the linear growth of the average square of the displacement with time, and the mean (k) was found for the growth aligned with K_c with (c) where c ∈ ( 0 , π ) , as it was Getting k = 0.9305 ≅ 1 , which indicates chaotic of the system (1), shown in Figure 7.

To address the chaos of system (1) we will use the adaptive control technique and design an adaptive control law with the unknown parameter b of the system.

By adding the control units to the system, we get:-

{ x t + 1 = x t + y t + u 1 y t + 1 = 0.1 x t 2 + 0.1 + y t ( − x t − b y t ) + y t + u 2 (8)

where u 1 , u 2 are the controllers for the adaptive feeding, defined as follows:-

{ u 1 = − x t − y t − M 1 x t u 2 = − 0.1 x t 2 − 0.1 + x t y t + β y t 2 − y t − M 1 y t (9)

where β is an approximate parameter of b and M 1 , M 2 are positive constants by substituting (9) in (8) we get:

{ x t + 1 = − M 1 x t y t + 1 = ( − b + β ) y t 2 − M 2 y t (10)

Let the error for the discretionary parameter as follows:

e b = − b + β (11)

By substituting (11) into (10):

{ x t + 1 = − M 1 x t y t + 1 = e b y t 2 − M 2 y t (12)

In this paper, the stability and chaos of a two-dimensional discrete time dynamical system with hidden attractors was studied, and the fixed points were found: p 0 = ( i , 0 ) , p 1 = ( − i , 0 ) which shows us that the points are complex. The stability analysis of fixed points was analyzed by using the characteristic equation roots test, which shows us that the points are non-hyperbolic points, so we could not use jury test for the stability, the dynamic behavior of system was analyzed and a time series was generated for each state of system and phase space and bifurcation diagrams of the system were found. Chaos diagnosis in the system by using Lyapunov exponent test and it was obtained: L 1 = − 0.250000 , L 2 = 1.853981 . Lyapunov dimension: D L = 0.865156 , and from the binary test (0 - 1) we get the index value k = 0.9305, which is close to one, all show that the system (1) is chaotic. The adaptive control technique of the chaotic system was carried out. Stability tests: roots of the characteristic equation and jury test for points in the controlling system (10), and Lyapunov exponent test obtaining the values: L 1 = − 0.700000 , L 2 = − 0.41886 , which all show that the system after control is stable and regular.

The authors acknowledge the support university of Mosul and college of computer science and mathematics.

The authors declare no conflicts of interest regarding the publication of this paper.

Aziz, M.M. and Jihad, O.M. (2021) Stability & Chaos Tests of 2D Discrete Time Dynamical System with Hidden Attractors. Open Access Library Journal, 8: e7501. https://doi.org/10.4236/oalib.1107501