Next Article in Journal
Aggregative Game for Distributed Charging Strategy of PEVs in a Smart Charging Station
Next Article in Special Issue
Boundary Value Problems for Fractional Differential Equations of Caputo Type and Ulam Type Stability: Basic Concepts and Study
Previous Article in Journal
Mathematical Modeling of the Tumor–Immune System with Time Delay and Diffusion
Previous Article in Special Issue
Lower and Upper Solution Method for Semilinear, Quasi-Linear and Quadratic Singularly Perturbed Neumann Boundary Value Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Complex Dynamics of Rössler–Nikolov–Clodong O Hyperchaotic System: Analysis and Computations

by
Svetoslav G. Nikolov
1,2,* and
Vassil M. Vassilev
1
1
Institute of Mechanics, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Bl. 4, 1113 Sofia, Bulgaria
2
Department of Mechanics, University of Transport, Geo Milev Str., 158, 1574 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(2), 185; https://doi.org/10.3390/axioms12020185
Submission received: 9 January 2023 / Revised: 3 February 2023 / Accepted: 8 February 2023 / Published: 10 February 2023

Abstract

:
This paper discusses the analysis and computations of chaos–hyperchaos (or vice versa) transition in Rössler–Nikolov–Clodong O (RNC-O) hyperchaotic system. Our work is motivated by our previous analysis of hyperchaotic transitional regimes of RNC-O system and the results recently obtained from another researchers. The analysis and numerical simulations show that chaos–hyperchaos transition in RNC-O system is coupled to change in the equilibria type as one large hyperchaotic attractor occurs. Moreover, we show that for this system, a zero-Hopf bifurcation is not possible. We also consider the cases when the divergence of the system is a constant and detected two families of exact solutions.

1. Introduction

Autonomous nonlinear four-dimensional systems of the form
d x d t = x ˙ = f ( x )     with     f ( x ) = ( f 1 ( x ) , f 2 ( x ) , f 3 ( x ) , f 4 ( x ) ) ;   x = ( x 1 , x 2 , x 3 , x 4 )
can display a rich diversity of periodic, multiple periodic, chaotic and hyperchaotic solutions dependent upon the specific values of one or more bifurcation (control) parameters [1,2,3,4,5]. Such a system is characterized by four Lyapunov exponents λ 1 , λ 2 , λ 3 and λ 4 . If i = 1 4 λ i < 0 , then the system is dissipative and its evolution takes place on an attractor. Now, it is well known that attractors characterized by only one positive Lyapunov exponents are chaotic and they can occur in at least three-dimensional nonlinear systems [5,6,7,8]. On the other hand, for continuous nonlinear four-dimensional dissipative systems, we can observe attractors with two positive Lyapunov exponents [1,9,10]. Historically, from Rössler [11], these systems are for first time named hyperchaotic. For hyperchaotic systems of the form in (1), the possible types of Lyapunov spectra are ( + , + , 0 , ) , ( + , 0 , 0 , ) and ( + , 0 , , ) [1,6,9,11,12].
Recently, a growing interest has been observed in higher-dimensional systems greater than three with chaotic and hyperchaotic behavior [3,7,8,13,14,15,16,17,18,19,20,21,22]. The main reasons are the following: (1) mathematical models of natural phenomena are high dimensional; (2) the possible bifurcations of four-dimensional systems are absent in three-dimensional ones; (3) the characteristic properties (for systems of dimension higher than three) are covering without folds and (4) in high-dimensional systems, a hyperchaotic behavior occurs, which is essentially different from the chaos in low-dimensional models.
Since Meyer et al. [23] and Nikolov [9,24], the transitional processes of hyperchaotic Rössler type systems have been widely considered. Analytical and computational approaches have become increasingly important for nonlinear sciences and are found in more applications in all scientific areas [7,17,22,25,26,27,28]. During the last 10 years, many studies were devoted to the bifurcation scenarios of hyperchaotic systems [13,16,18,19]. This is possible as long as the transition from chaos to hyperchaos (or vice versa) is analytically studied for the case of a small number of parameters (codimension), since in more complex cases, the domain remains widely open.
Recently, in the hyperchaotic Rössler system, main bifurcation scenario leading to a hyperchaos has been demonstrated by Stankevich et al. [13]. The main bifurcation hyperchaotic scenario can be divided into two parts: the first one includes few local and global bifurcations as a homoclinic attractor is created; the second one (which is more complicated) includes cascades of local and global bifurcations as after an external bifurcation, the system becomes hyperchaotic.
Here, we consider the following system:
x ˙ = y z , y ˙ = x + a y + w , z ˙ = b + b 1 x + b 2 y + b 3 z + b 4 w + x z , w ˙ = c z + d w ,
where a = 0.25 , b = 3 , c = 0.5 , d = 0.05 [11], and b 1 , b 2 , b 3 and b 4 are real constants. This system was introduced in [2] as a modification of the original Rössler hyperchaotic system named Rössler–Nikolov–Clodong O hyperchaotic system (RNC-O). In the paper comparing the dynamics of the original and the modified system, it was determined that (1) RNC-O only exhibits hyperchaotic behavior when b 3 is varied; (2) repeated hyperchaos–chaos–hyperchaos transition can be seen; (3) the hyperchaos is possible for a smaller value of information dimension.
In general, system (2) has two equilibrium (fixed) points ( x ¯ 1 , 2 , y ¯ 1 , 2 , z ¯ 1 , 2 , w ¯ 1 , 2 ) , where
x ¯ 1 , 2 = G ± g + G 2 2 d , y ¯ 1 , 2 = G ± g + G 2 2 ( c a d ) , z ¯ 1 , 2 = G ± g + G 2 2 ( c a d ) , w ¯ 1 , 2 = c G ± g + G 2 2 d ( c a d ) ,
and
g = 4 b d ( c a d ) , G = b 4 c b 1 ( c a d ) d ( b 2 b 3 ) ,
which reduce to a single one if
g + G 2 = 0 .
It is easy to check that in the considered case ( a = 0.25 , b = 3 , c = 0.5 , d = 0.05 [11]) g + G 2 > 0 , and hence x ¯ 1 , 2 , y ¯ 1 , 2 , z ¯ 1 , 2 and w ¯ 1 , 2 are all real numbers. Thus, in this case, there are two different equilibrium points. Note that in Figure A1 (see Appendix A), the values of z ¯ 1 , 2 are shown, and those for which the system (2) exhibits hyperchaotic behavior are marked with red color.
In Section 2, we analyze the system (2) in an adiabatic case for its divergence. In Section 3, we investigate equilibria and the coefficient matrix of the linearized system, as we discuss the possible types of transitions to hyperchaos occurring in system (2). In Section 4, numerical examples are presented with some simulation results to illustrate the theoretical results of previous two sections. Finally, we conclude with a summary and some comments in Section 5.

2. Families of Exact Solutions

The divergence of the flow (2) is obtained by
D 4 = x ˙ x + y ˙ y + z ˙ z + w ˙ w = a + d + b 3 + x = 0.3 + b 3 + x .
It can be seen that the RNC-O system has an alternating divergence—see Figure A2 and Figure A3 in Appendix C. The system (2) is dissipative when D 4 < 0 , and an attractor exists if we determined the absorbing domain.
Interesting families of exact solutions to systems of form (2), which are not all constants, arise under the assumption that x = k = c o n s t . Then, x ˙ = 0 , z = y , z ˙ = y ˙ , and Equation (2) transforms into the following overdetermined system:
y ˙ = k + a y + w ,
y ˙ = b b 1 k ( b 2 b 3 k ) y b 4 w
w ˙ = c y + d w
Let b 4 1 and k b 2 b 3 ± a . Then, after subtraction of (7) from (8) and solving for w , one obtains
w = b + ( b 1 + 1 ) k 1 + b 4 b 2 b 3 + a k 1 + b 4 y
Consequently,
w ˙ = b 2 b 3 + a k 1 + b 4 y ˙
and the above system (7)–(9) becomes
y ˙ = b + ( b 1 b 4 ) k 1 + b 4 b 2 b 3 a b 4 k 1 + b 4 y
y ˙ = d b + ( b 1 + 1 ) k b 2 b 3 + a k + ( d ( 1 + b 4 ) c b 2 b 3 + a k ) y
Evidently, Equations (12) and (13) are compatible if and only if
d b + ( b 1 + 1 ) k b 2 b 3 + a k = b + ( b 1 b 4 ) k 1 + b 4
and
d ( 1 + b 4 ) c b 2 b 3 + a k = b 2 b 3 a b 4 k 1 + b 4
Now, it is easy to see that
x = k , y = b + ( b 1 b 4 ) k b 2 b 3 a b 4 k + c 1 e b 2 b 3 a b 4 k 1 + b 4 t , z = b + ( b 1 b 4 ) k b 2 b 3 a b 4 k c 1 e b 2 b 3 a b 4 k 1 + b 4 t , w = a ( b + b 1 k ) ( b 2 b 3 k ) k b 2 b 3 a b 4 k c 1 b 2 b 3 + a k 1 + b 4 e b 2 b 3 a b 4 k 1 + b 4 t ,
where c 1 is an arbitrary constant, is a one-parameter family of solutions of system (2), provided that the compatibility conditions (14) and (15) hold.
Next, let b 4 = 1 . Then, subtracting Equation (8) from Equation (7), we obtain
b + ( b 1 + 1 ) k + ( b 2 b 3 + a k ) y = 0
Hence, b + ( b 1 + 1 ) k = 0 and b 2 b 3 + a k = 0 are the necessary and sufficient conditions for the compatibility of the aforementioned two equations. Suppose now that
k = b 2 b 3 + a , b 1 = 1 b b 2 b 3 + a
Then, solving Equation (7) for w and taking into account the above relations, we obtain
w = y ˙ a y a b 2 + b 3
and
w ˙ = y ¨ a y ˙
Substituting (19) and (20) into (9), we obtain the following equation:
y ¨ ( a + d ) y ˙ + ( c + a d ) y + ( b 2 b 3 + a ) d = 0
whose general solution is
y = A 0 + c 2 e A t + c 3 e B t
where c 2 and c 3 are arbitrary constants, and
A 0 = b 2 b 3 + a c + a d d
A = 1 2 ( a + d ( a d ) 2 4 c ) , B = 1 2 ( a + d + ( a d ) 2 4 c )
Finally, using (19) and (22), we obtain the following two-parameter family of solutions of system (2):
x = k = b 2 b 3 + a , y = A 0 + c 2 e A t + c 3 e B t , z = A 0 c 2 e A t c 3 e B t , w = a ( 1 + A 0 ) b 2 + b 3 + c 2 ( A a ) e A t + c 3 ( B a ) e B t .
Some exact solutions of (2) are shown in Figure A4, Figure A5 and Figure A6 (see Appendix D).

3. Qualitative Behavior for System (2) in Dissipative Case

3.1. Local Analysis

It easy to express (2) in terms of local coordinates x 1 , x 2 , x 3 and x 4 , where x ¯ , y ¯ , z ¯ and w ¯ were defined in (3)
x = x ¯ + x 1 , y = y ¯ + x 2 , z = z ¯ + x 3 , w = w ¯ + x 4
Hence, after some transformations, we determine that
x ˙ 1 = x 2 x 3 , x ˙ 2 = x 1 + a x 2 + x 4 , x ˙ 3 = k 1 x 1 + b 2 x 2 + k 2 x 3 + b 4 x 4 + x 1 x 3 , x ˙ 4 = c x 3 + d x 4 ,
where k 1 = b 1 + z ¯ and k 2 = b 3 + x ¯ . The eigenvalues of the Jacobian matrix of system (6) (for the equilibria (3)) can be computed through the determinant
Ξ ( χ ) = | χ 1 1 0 1 a χ 0 1 k 1 b 2 k 2 χ b 4 0 0 c d χ | = 0 ,
which is equivalent to the characteristic equation
χ 4 + p χ 3 + q χ 2 + r χ + s = 0
where
p = ( a + d + k 2 ) , q = 1 + k 1 + a k 2 + d ( a + k 2 ) + c b 4 , r = b 2 ( 1 + c ) k 2 ( 1 + a d ) k 1 ( a + d ) d a c b 4 , s = k 1 ( a d c ) + d ( k 2 b 2 ) + c b 4 .
For n-dimensional dynamic system, according to [9,29,30,31], in principle, six possible equilibrium types (cases) (with even more sub-cases) depending on the roots of the characteristic equation exist. In four-dimensional systems, there are five distinct types (for t + ) at a hyperbolic fixed (equilibrium) point with even more sub-cases. According to [32], they are
r e a l e i g e n v a l u e : { 1 ) χ 1 , 2 > 0 , χ 3 , 4 < 0 . 2 ) χ 1 > 0 , χ 2 , 3 , 4 < 0 .   O 2 , 2 O 3 , 1
c o m p l e x e i g e n v a l u e s : 3 ) χ 1 , 2 = n 1 ± i m 1 , χ 3 , 4 = n 2 ± i m 2 ; n 1 > 0 , n 2 < 0 , m 1 , m 2 0 .   O 2 , 2
r e a l a n d c o m p l e x e i g e n v a l u e s : | 4 ) χ 1 , 2 = n 1 ± i m 1 , χ 3 , 4 > 0 ; n 1 < 0 , m 1 0 . 5 ) χ 1 , 2 = n 1 ± i m 1 , χ 3 < 0 , χ 4 > 0 ; n 1 < 0 , m 1 0 .   O   ,   2 , 2 O 3 , 1
In all cases when u + v = n , u 0 , and v 0 , the equilibrium state O u , v is (compound) saddle of first/second order. Moreover, for t , we have O 2 , 2 and O 1 , 3 . If u + v n , i.e., there are zero(s) root(s), then compound saddle-focus, compound saddle-knot, etc., take place. For more details, see Case 6 in the Appendix of [9]. Note here that if we have purely imaginary roots and roots whose real parts are of different signs, then the equilibrium state is of a complex saddle type (CS).
Figure 1 shows the possible types (cases) for first ( O 1 ) and second ( O 2 ) equilibria of system (2) depending on the roots of the characteristic Equation (29), when b 1 [ 0.01 , 1.5 ] , b 2 [ 2 , 0.01 ] , b 3 [ 0.01 , 2.5 ] , b 4 [ 0.01 , 0.7 ] and all another parameters values are as those in original Rössler hyperchaotic system. It can be seen that possible types of O 1 and O 2 are one and the same. Interestingly, the dominant type of O 2 is unstable knot ( O 2 0 , 4 ), and the transition between O 2 3 , 1 and O 2 1 , 3 is not typical, as remarked with dotted line. Note that diagrams in Figure 1 are obtained as a result of direct numerical calculations of (29), not shown here. According to [2,9], transition from chaos to hyperchaos (or vice versa) takes place for system (2) when equilibrium points change its type from compound saddle to compound knot (or vice versa).
In order to understand mechanisms of the emergence of hyperchaotic attractors, we will analyze the bifurcation of one-parameter families of the model. If a four-dimensional nonlinear system has an isolated equilibrium with double-zero eigenvalue and a pair of purely imaginary eigenvalues (named zero-Hopf equilibrium), then, under certain conditions, some complicated invariant sets could be bifurcated from a zero-Hopf bifurcation [33]. Thus, the zero-Hopf equilibrium, in some cases, can be generator of a local chaos (see [34] and references therein).
Let us come back to the study of the equilibrium points (3) and the characteristic Equation (29). If we assume that b = c = b 1 = b 2 = b 3 = b 4 = 0 into (3), then we always have an equilibrium point for the hyperchaotic system (2) localized at the origin of its coordinates, i.e., we have E 0 ( 0 , 0 , 0 , 0 ) . In light of this, the following problem arises:
Problem 1.
Is there a one-parameter family of systems of form (2) for which the origin of coordinates (i.e., E 0 ) is a zero-Hopf equilibrium point?
A general answer to this problem is provided below.
Jacobian matrix (28) evaluated at E 0 is
Ξ 0 ( χ ) = | χ 1 1 0 1 a χ 0 1 0 0 χ 0 0 0 0 d χ | = 0 ,
which is equivalent to the characteristic equation
χ 4 + p 0 χ 3 + q 0 χ 2 + r 0 χ + s 0 = 0
where
p 0 = ( a + d ) , q 0 = 1 + a d , r 0 = d , s 0 = 0 .
To satisfy the conditions for E 0 to be an isolated equilibrium with double-zero eigenvalue and a pair of purely imaginary eigenvalues, the coefficients p 0 , q 0 , r 0 and s 0 must be p 0 = 0 , q 0 < 0 , r 0 = 0 and s 0 = 0 , i.e., we must have a = d = 0 . It is easy to see that q 0 = 1 > 0 . Hence, we conclude that the conditions for E 0 to be isolated zero-Hopf equilibrium are not valid, and it follows that the zero-Hopf bifurcation for system (2) is not possible. Then, the following theorem can be formulated:
Theorem 1.
If E 0 is an equilibrium point for the hyperchaotic system (2) localized at the origin of its coordinates, then this equilibrium cannot be of zero-Hopf type. This means that for the system (2), a zero-Hopf bifurcation is not possible.

3.2. Global Behavior

According to [13], the main bifurcation scenarios leading to hyperchaotic behavior can be naturally divided into two parts. The first one includes only a few local and global bifurcations leading to the emergence (in the Poincare map of the system) of a homoclinic attractor—hyperchaotic/chaotic. The second one (which is too complicated) includes (infinite) cascades of local and global bifurcations forming the skeleton of an attractor (with two positive Lyapunov exponents) which becomes hyperchaotic due to an external bifurcation (crisis). As a result, the smaller hyperchaotic attractor is destroyed, as its orbits transfer hyperchaoticity to the larger attractor. In our case, we obtain the same scenario; more precisely, in the second part, the larger attractor (area near O 1 ) starts to become dominant. After the crisis (an external bifurcation for the smaller attractor (area near O 2 )), the smaller attractor is destroyed and finally disappears. In other words, the two attractors merge into one large hyperchaotic attractor [23]. Note that equilibrium state O 2 is mainly a compound unstable knot— O 2 0 , 4 . This type of behavior is demonstrated in Section 4.
We may summarize the (precise) features described in the second part of the scenario leading to hyperchaos: when the dimension of the system is higher than four (all odd dimensions), then the number of smaller hyperchaotic attractors can be two, three or more, and one of them starts becoming dominant, i.e., the shadowland effect takes place. According to [23], for the generalized Rössler system of dimension n = 7 , two smaller and one large attractor exist. The three attractors, their existence, form and location in phase space are the aftermath of the structural symmetry. In the bifurcation diagram (see Figure 12 in [23]), attractor 1, attractor 2, and attractor 3 (from bottom to top) are in different shadings.

4. Numerical Results

Numerical experiments (results) are very useful for creating basic intuitions and suggesting problems for analytical investigations in chaotic dynamic. It is by far truer in the case of nonlinear systems in R 4 , which is the subject of our study.
According to [2,9], when Case 1  b 1 = 0.1 o r 1.2 , b 2 = 0.1 , b 3 = b 4 = 0.1 ; Case 2  b 1 = 0.2 , b 2 = 0.1 , b 3 = b 4 = 0.1 ; and Case 3  b 1 = b 3 = 0.2 , b 2 = 0.3 , b 4 = 0.01 , the system (2) exhibits hyperchaotic behavior. Here, the numerical computations are used to generate pictures of attractors for system (2) and to determined possible equilibria types. We fix the values of a , b , c and d for system (2), and investigate the changes in its behavior with control (bifurcation) parameters b 1 , b 2 , b 3 , and b 4 . The solutions for ( x ( t ) , y ( t ) , z ( t ) , w ( t ) ) have been calculated with initial conditions x 0 = 20 , y 0 = 0 , z 0 = 0 , w 0 = 15 using MATLAB (the MathWorks, Inc., Natick, MA, USA)- ODE45 solver (which is based on an explicit Runge–Kutta (4,5) formula, the Dormand–Prince pair).
In Figure 2, the phase space (in ( x , y , z ) and ( x , y , w ) coordinates) and possible positions of equilibria O 1 (red color point) and O 2 (blue color point) of system (2) in Case 1 are shown. Here, O 1 is compound saddle (i.e., O 1 1 , 3 or O 1 3 , 1 ), and O 2 can be only compound unstable knot (i.e., O 2 0 , 4 ) for b 1 = 0.1 (see Figure 2a,b), and for b 1 = 1.2 (see Figure 2c,d). Note that eigenvalues are provided in Appendix B. These figures suggest that a larger (dominant) attractor occurs near the location (basin of attraction) of O 1 . We see that the qualitative results of Section 3 well predict the behavior of system (2) near the chaos-hyperchaos (or vice versa) transitional point. Interestingly, some of the values of D 4 (for b 1 = 0.1 or b 1 = 1.2 ) in Case 1 can be positive and zero; see Figure A2 in Appendix C.
Figure 3 shows the phase spaces of system (2) in Case 2. The equilibrium points O 1 and O 2 were obtained by using (3). Similar to the results shown in Figure 2, here we see that the motion is also attracted to the first equilibrium point O 1 as a larger (dominant) attractor occurs. The result in Figure 3b clearly shows the way in which the smaller hyperchaotic attractor disappears. In this case, O 1 is compound saddle (i.e., O 1 1 , 3 ), and O 2 is compound unstable knot (i.e., O 2 0 , 4 ); see Appendix B.
Finally, we investigate Case 3 when b 1 = b 3 = 0.2 , b 2 = 0.3 , b 4 = 0.01 . Figure 4 shows the phase spaces of system (2). As can be seen, similar to the previous Figure 2 and Figure 3, the motion is attracted to the first equilibrium point O 1 , which is compound saddle (i.e., O 1 3 , 1 ). Here, O 2 is also compound unstable knot (i.e., O 2 0 , 4 ); see Appendix B.
Comparing Figure 2, Figure 3 and Figure 4, we conclude that our hypothesis appears to be realistic; that is, the hyperchaotic behavior of the system (2) depends on the type of equilibria—compound saddle and compound knot. As a result, a larger (dominant) attractor occurs near the location (basin) of O 1 . Note that this conclusion is in accordance with the scenario described in Section 3.2.

5. Conclusions

The present paper provides some new results and can be seen in part as an extension of earlier works on the RNC-O system (2) in [2,9]. Our main motivation comes from previous theoretical [13,16] and numerical works [9,24] which showed that specific cascades of local and global bifurcations might occur.
According to [35,36], for systems with an alternating divergence, the existence of new very unexpected and interesting dynamical regimes can be seen. To understand the effect of this alternating divergence on the dynamical organization of RNC-O system, in Section 2, we obtain two families of exact solutions of systems (2) under the assumption that its divergence is a constant, i.e., x = k = c o n s t .
As mentioned in Section 3, our goal was to study the equilibria type changes associated to a chaos–hyperchaos transition problem in RNC-O system. Thus, we determine that when equilibrium points change its type from compound saddle to compound knot (or vice versa), the hyperchaos occurs. Moreover, in Section 3.2, we precise the Stankevich’s scenario for hyperchaos transition [13] when one system displays a number of dimensions higher than four, i.e., the so-called shadowland effect takes place. Finally, in Section 3.1, we demonstrate that for RNC-O system, a zero-Hopf bifurcation is not possible. Hyperchaos behavior raises a number of further questions that we leave to future research.

Author Contributions

The initial idea for investigation comes from S.G.N. The model was computed and simulated by all co-authors. Qualitative analysis was performed by S.G.N., and exact solutions were calculated by V.M.V. The manuscript was drafted by all coauthors. All authors have read and agreed to the published version of the manuscript.

Funding

Author Vassil M. Vassilev wishes to acknowledge support for this work in the Bulgarian National Science Fund, grant КП-06-H 22-2.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. z ¯ Coordinates of the First and the Second Equilibrium Points and Values of Lyapunov Exponents

Figure A1. Equilibrium points z ¯ 1 , z ¯ 2 (a) and Lyapunov exponents spectrum (b) when b 2 = 0.1 , b 3 = b 4 = 0.1 and b 1 [ 0.01 , 1.5 ] . All other parameter values are same as those in original Rössler hyperchaotic system; see [11]. Note that red color (see (a)) marks those equilibria for which the system (2) exhibits hyperchaotic behavior.
Figure A1. Equilibrium points z ¯ 1 , z ¯ 2 (a) and Lyapunov exponents spectrum (b) when b 2 = 0.1 , b 3 = b 4 = 0.1 and b 1 [ 0.01 , 1.5 ] . All other parameter values are same as those in original Rössler hyperchaotic system; see [11]. Note that red color (see (a)) marks those equilibria for which the system (2) exhibits hyperchaotic behavior.
Axioms 12 00185 g0a1

Appendix B. Eigenvalues of the First and the Second Equilibrium Points

For system (2), if (1.1) b 1 = b 3 = b 4 = 0.1 , b 2 = 0.1 ; (1.2) b 1 = 1.2 , b 3 = b 4 = 0.1 , b 2 = 0.1 ; (2) b 1 = 0.2 , b 3 = b 4 = 0.1 , b 2 = 0.1 and (3) b 1 = b 3 = 0.2 , b 2 = 0.3 , b 4 = 0.01 , the equilibrium points and their eigenvalues are provided by the following:
Case 1.1
O 1 ( 5.52 , 0.5664 , 0.5664 , 5.66 ) ,     then     ( χ 1 , χ 2 , χ 3 , χ 4 ) = ( 5.29 , 0.031 ± 0.986 i , 0.105 ) ; O 2 ( 5.3 , 0.5433 , 0.5433 , 5.43 ) ,   then   ( χ 1 , χ 2 , χ 3 , χ 4 ) = ( 5.47 , 0.063 ± 0.989 i , 0.1 ) ;
Case 1.2
O 1 ( 2.29 , 0.2346 , 0.2346 , 2.35 )     then     ( χ 1 , χ 2 , χ 3 , χ 4 ) = ( 1.554 , 0.325 ± 1.193 i , 0.317 ) ; O 2 ( 12.79 , 1.312 , 1.312 , 13.12 )     then     ( χ 1 , χ 2 , χ 3 , χ 4 ) = ( 12.89 , 0.118 ± 0.997 i , 0.058 ) ;
Case 2
O 1 ( 5.05 , 0.518 , 0.518 , 5.18 )     then     ( χ 1 , χ 2 , χ 3 , χ 4 ) = ( 4.79 , 0.014 ± 0.988 i , 0.116 ) ; O 2 ( 5.8 , 0.595 , 0.595 , 5.95 )     then     ( χ 1 , χ 2 , χ 3 , χ 4 ) = ( 5.96 , 0.075 ± 0.992 i , 0.092 ) ;
Case 3
O 1 ( 4.775 , 0.4898 , 0.4898 , 4.897 )     then     ( χ 1 , χ 2 , χ 3 , χ 4 ) = ( 4.408 , 0.0011 ± 0.959 i , 0.135 ) ; O 2 ( 6.13 , 0.628 , 0.628 , 6.282 )     then     ( χ 1 , χ 2 , χ 3 , χ 4 ) = ( 6.4 , 0.071 ± 0.0143 i , 0.0824 ) .

Appendix C. Values of D 4 ; see (6)

Figure A2. Values of D4 (see (6)) when (a) a = 0.25, b = 3, c = 0.5, d = 0.05, b1 = 0.1, b2 = −0.1, b3 = b4 = 0.1; (b) a = 0.25, b = 3, c = 0.5, d = 0.05, b1 = 1.2, b2 = −0.1, b3 = b4 = 0.1.
Figure A2. Values of D4 (see (6)) when (a) a = 0.25, b = 3, c = 0.5, d = 0.05, b1 = 0.1, b2 = −0.1, b3 = b4 = 0.1; (b) a = 0.25, b = 3, c = 0.5, d = 0.05, b1 = 1.2, b2 = −0.1, b3 = b4 = 0.1.
Axioms 12 00185 g0a2
Figure A3. Values of D4 (see (6)) when (a) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 ,   b 1 = 0.2 , b 2 = 0.1 , b 3 = b 4 = 0.1 ; (b) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 ,   b 1 = 0.2 , b 2 = 0.3 , b 3 = 0.2 , b 4 = 0.01 .
Figure A3. Values of D4 (see (6)) when (a) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 ,   b 1 = 0.2 , b 2 = 0.1 , b 3 = b 4 = 0.1 ; (b) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 ,   b 1 = 0.2 , b 2 = 0.3 , b 3 = 0.2 , b 4 = 0.01 .
Axioms 12 00185 g0a3

Appendix D. Families of Exact Solutions

Figure A4. Exact solutions of system (2) when b1 = 2.266, b2 = 1, b3 = 2.092, b4 = 2, k = 1, c1 = 1, according to Formula (16).
Figure A4. Exact solutions of system (2) when b1 = 2.266, b2 = 1, b3 = 2.092, b4 = 2, k = 1, c1 = 1, according to Formula (16).
Axioms 12 00185 g0a4
Figure A5. Exact solutions of system (2) when b1 = −1.942, b2 = 1, b3 = 1.064, b4 = −2, k = 1, c1 = 1, according to Formula (16).
Figure A5. Exact solutions of system (2) when b1 = −1.942, b2 = 1, b3 = 1.064, b4 = −2, k = 1, c1 = 1, according to Formula (16).
Axioms 12 00185 g0a5
Figure A6. Exact solutions of system (2) when b1 = 3, b2 = 1, b3 = 2, b4 = −1, k = −0.75, c2 = 1, c3 = 0, according to Formulas (18) and (23)–(25).
Figure A6. Exact solutions of system (2) when b1 = 3, b2 = 1, b3 = 2, b4 = −1, k = −0.75, c2 = 1, c3 = 0, according to Formulas (18) and (23)–(25).
Axioms 12 00185 g0a6

References

  1. Letellier, C.; Rossler, O.E. Hyperchaos. Scholarpedia 2007, 2, 1936. [Google Scholar] [CrossRef]
  2. Nikolov, S.; Clodong, S. Occurrence of regular, chaotic and hyperchaotic behavior in a family of modified Rossler hyperchaotic systems. Chaos Solitons Fractals 2004, 22, 407–431. [Google Scholar] [CrossRef]
  3. Nikolov, S. Estimating of bifurcations and chaotic behavior in a four-dimensional system. J. Calcutta Math. Soc. 2006, 2, 17–28. [Google Scholar]
  4. Panchev, S. Theory of Chaos; Bulgarian Acad. Press: Sofia, Bulgaria, 2001. [Google Scholar]
  5. Peng, J.; Ding, E.; Ding, M.; Yang, W. Synchronizing hyperchaos with a scalar transmitted signal. Phys. Rev. Lett. 1996, 76, 904. [Google Scholar] [CrossRef] [PubMed]
  6. Wolf, A.; Swift, J.; Swinney, H.; Vastano, J. Determining Lyapunov exponents from a time series. Phys. D Nonlinear Phenom. 1985, 16, 285–317. [Google Scholar] [CrossRef]
  7. Jahanshahi, H.; Yousefpour, A.; Wei, Z.; Alcaraz, R.; Bekiros, S. A financial hyperchaotic system with coexisting attractors: Dynamic investigation, entropy analysis, control and synchronization. Chaos Solitons Fractals 2019, 126, 66–77. [Google Scholar] [CrossRef]
  8. Pavlov, A.; Pavlova, O.; Mohammad, Y.; Kurths, J. Characterization of the chaos-hyperchaos transition based on return times. Phys. Rev. E 2015, 91, 022921. [Google Scholar] [CrossRef]
  9. Nikolov, S.; Clodong, S. Hyperchaos-chaos-hyperchaos transition in modified Rossler type systems. Chaos Solitons Fractals 2006, 28, 252–263. [Google Scholar] [CrossRef]
  10. Alligood, K.; Sauer, T.; Yorke, J. Chaos, An Introduction to Dynamical System; Springer: New York, NY, USA, 1996. [Google Scholar]
  11. Rössler, O. An equation for hypechaos. Phys. Lett. A 1979, 71A, 155–157. [Google Scholar] [CrossRef]
  12. Wang, X.; Wang, M. A hyperchaos generalized from Lorenz system. Phys. A Stat. Mech. Its Appl. 2008, 387, 3751–3758. [Google Scholar] [CrossRef]
  13. Stankevich, N.; Kazakov, A.; Gonchenko, S. Scenarios of hyperchaos occurence in 4D Rössler system. Chaos Interdiscip. J. Nonlinear Sci. 2020, 30, 123129. [Google Scholar] [CrossRef]
  14. Starkov, K. On the ultimate dynamics of the four-dimensional Rössler system. Int. J. Bifurc. Chaos 2014, 24, 1450149. [Google Scholar] [CrossRef]
  15. Szczepaniak, A.; Macek, W. Unstable manifolds for the hyperchaotic Rossler system. Phys. Lett. A 2008, 372, 2423–2427. [Google Scholar] [CrossRef]
  16. Barrio, R.; Martinez, M.; Serrano, S.; Wilczak, D. When chaos meets hyperchaos: 4D Rossler model. Phys. Lett. A 2015, 379, 2300–2305. [Google Scholar] [CrossRef]
  17. Kuptsov, P.; Kuznetsov, S. Route to hyperbolic hyperchaos in a nonautonomous time-delay system. Chaos Interdiscip. J. Nonlinear Sci. 2020, 30, 113113. [Google Scholar] [CrossRef] [PubMed]
  18. Stankevich, N.; Volkov, E. Chaos–hyperchaos transition in three identical quorum-sensing mean-field coupled ring oscillators. Chaos Interdiscip. J. Nonlinear Sci. 2021, 31, 103112. [Google Scholar] [CrossRef] [PubMed]
  19. Sataev, I.; Stankevich, N. Cascade of torus birth bifurcations and inverse cascade of Shilnikov attractors merging at the threshold of hyperchaos. Chaos Interdiscip. J. Nonlinear Sci. 2021, 31, 023140. [Google Scholar] [CrossRef]
  20. Li, Q.; Tang, S.; Yang, X. Hyperchaotic set in continuous chaos–hyperchaos transition. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 3718–3734. [Google Scholar] [CrossRef]
  21. Karatetskaia, E.; Shykhmamedov, A.; Kazakov, A. Shilnikov attractors in three-dimensional orientation-reversing maps. Chaos Interdiscip. J. Nonlinear Sci. 2021, 31, 011102. [Google Scholar] [CrossRef]
  22. Lai, Q.; Norouzi, B.; Liu, F. Dynamic analysis, circuit realization, control desing and image encryption application of an extended Lü system with coexisting attractors. Chaos Solitons Fractals 2018, 114, 230–245. [Google Scholar] [CrossRef]
  23. Meyer, T.; Bünner, M.; Kittel, A.; Parisi, J. Hyperchaos in the generalized Rössler system. Phys. Rev. E 1997, 56, 5069–5082. [Google Scholar] [CrossRef] [Green Version]
  24. Nikolov, S. Transitional processes in some modified Rossler type dynamical systems. Comptes Rendus L’academie Bulg. Sci. 2004, 57, 45–52. [Google Scholar]
  25. Schuster, H. Deterministic Chaos: An introduction; Physik-Verlag: Berlin, Germany, 1984. [Google Scholar]
  26. Kapitaniak, T.; Maistrenko, Y.; Popovych, S. Chaos-hyperchaos transition. Phys. Rev. E 2000, 62, 1972–1976. [Google Scholar] [CrossRef] [PubMed]
  27. Hirsch, M.; Smale, S.; Devaney, R. Differential Equations, Dynamical Systems, and An Introduction to Chaos; Academic Press: New York, NY, USA, 2012. [Google Scholar]
  28. Singh, S.; Han, S.; Lee, S.M. Adaptive single input sliding mode control for hybrid-synchronization of uncertain hyperchaotic Lu systems. J. Frankl. Inst. 2021, 358, 7468–7484. [Google Scholar] [CrossRef]
  29. Nemytskii, V.; Stepanov, V. Qualitative Theory of Differential Equations; Princeton University Press: Princeton, NJ, USA, 1960. [Google Scholar]
  30. Mints, R. The character of certain types of complex equilibrium state in n-dimensional spaces. DAN USSR 1962, 147, 31–33. [Google Scholar]
  31. Gavrilov, N. On n-dimensional dynamic systems close to the systems with non-rough homoclinic curve. DAN USSR 1973, 212, 276–279. [Google Scholar]
  32. Wiggins, S. Global Bifurcations and Chaos: Analytical Methods; Springer: New York, NY, USA, 1988. [Google Scholar]
  33. Cid-Montiel, L.; Llibre, J.; Stoica, C. Zero-Hopf bifurcation in a hyperchaotic Lorenz system. Nonlinear Dyn. 2014, 75, 561–566. [Google Scholar] [CrossRef]
  34. Champney, A.; Kirk, V. The entwined wiggling of homoclinic curves emerging from saddle-node/Hopf instabilities. Phys. D Nonlinear Phenom. 2004, 195, 77–105. [Google Scholar] [CrossRef]
  35. Gonchenko, S.; Turaev, D. On three types of dynamics and notation of attractor. Proc. Steklov Inst. Math. 2017, 297, 116–137. [Google Scholar] [CrossRef]
  36. Gonchenko, S.; Gonchenko, A.; Kozlov, A. Three types of attractors and mixed dynamics of nonholonomic models of rigid body motion. Proc. Steklov Inst. Math. 2020, 308, 125–140. [Google Scholar] [CrossRef]
Figure 1. Possible types (cases) and transitions for the first ( O 1 ) and the second ( O 2 ) equilibria of system (2) depending on the roots of the characteristic Equation (29), when b 1 [ 0.01 , 1.5 ] , b 2 [ 2 , 0.01 ] , b 3 [ 0.01 , 2.5 ] , b 4 [ 0.01 , 0.7 ] and all another parameter values are as those in the original Rössler hyperchaotic system, see in [11]. The dotted line remarks that transition between O 2 3 , 1 and O 2 1 , 3 can be seen only in a small number of values for b 1 b 4 and that is why it is not typical.
Figure 1. Possible types (cases) and transitions for the first ( O 1 ) and the second ( O 2 ) equilibria of system (2) depending on the roots of the characteristic Equation (29), when b 1 [ 0.01 , 1.5 ] , b 2 [ 2 , 0.01 ] , b 3 [ 0.01 , 2.5 ] , b 4 [ 0.01 , 0.7 ] and all another parameter values are as those in the original Rössler hyperchaotic system, see in [11]. The dotted line remarks that transition between O 2 3 , 1 and O 2 1 , 3 can be seen only in a small number of values for b 1 b 4 and that is why it is not typical.
Axioms 12 00185 g001
Figure 2. Phase space of system (2) in hyperchaotic regime when (a,b) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 , b 2 = 0.1 , b 3 = 0.1 , b 4 = 0.1 with b 1 = 0.1 ; (c,d) b 1 = 1.2 . The first equilibrium O 1 is marked with red color and the second equilibrium O 2 is marked with blue color.
Figure 2. Phase space of system (2) in hyperchaotic regime when (a,b) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 , b 2 = 0.1 , b 3 = 0.1 , b 4 = 0.1 with b 1 = 0.1 ; (c,d) b 1 = 1.2 . The first equilibrium O 1 is marked with red color and the second equilibrium O 2 is marked with blue color.
Axioms 12 00185 g002
Figure 3. Phase space of system (2) in hyperchaotic regime, when (a,b) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 , b 1 = 0.2 , b 2 = 0.1 , b 3 = b 4 = 0.1 . The first equilibrium O 1 is marked with red color and the second equilibrium O 2 is marked with blue color.
Figure 3. Phase space of system (2) in hyperchaotic regime, when (a,b) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 , b 1 = 0.2 , b 2 = 0.1 , b 3 = b 4 = 0.1 . The first equilibrium O 1 is marked with red color and the second equilibrium O 2 is marked with blue color.
Axioms 12 00185 g003
Figure 4. Phase space of system (2) in hyperchaotic regime, when (a,b) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 , b 1 = 0.2 , b 2 = 0.3 , b 3 = 0.2 , b 4 = 0.01 . The first equilibrium O 1 is marked with red color and the second equilibrium O 2 is marked with blue color.
Figure 4. Phase space of system (2) in hyperchaotic regime, when (a,b) a = 0.25 , b = 3 , c = 0.5 , d = 0.05 , b 1 = 0.2 , b 2 = 0.3 , b 3 = 0.2 , b 4 = 0.01 . The first equilibrium O 1 is marked with red color and the second equilibrium O 2 is marked with blue color.
Axioms 12 00185 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nikolov, S.G.; Vassilev, V.M. Complex Dynamics of Rössler–Nikolov–Clodong O Hyperchaotic System: Analysis and Computations. Axioms 2023, 12, 185. https://doi.org/10.3390/axioms12020185

AMA Style

Nikolov SG, Vassilev VM. Complex Dynamics of Rössler–Nikolov–Clodong O Hyperchaotic System: Analysis and Computations. Axioms. 2023; 12(2):185. https://doi.org/10.3390/axioms12020185

Chicago/Turabian Style

Nikolov, Svetoslav G., and Vassil M. Vassilev. 2023. "Complex Dynamics of Rössler–Nikolov–Clodong O Hyperchaotic System: Analysis and Computations" Axioms 12, no. 2: 185. https://doi.org/10.3390/axioms12020185

APA Style

Nikolov, S. G., & Vassilev, V. M. (2023). Complex Dynamics of Rössler–Nikolov–Clodong O Hyperchaotic System: Analysis and Computations. Axioms, 12(2), 185. https://doi.org/10.3390/axioms12020185

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop