On Universality of Transition to Chaos Scenario in Nonlinear Systems of Ordinary Differential Equations of Shilnikov’s Type ()
Received 16 March 2016; accepted 20 May 2016; published 23 May 2016

1. Introduction
In recent years, many scientists around the world have tried to analyze and classify the variety of systems with chaotic behavior. In the article [3] , Chen studied the problem of classification of quadratic autonomous dynamic systems of ordinary differential equations, in accordance with the extended Shilnikov’s theorem. Systems of Shilnikov’s type include a lot of well-known chaotic systems, such as the Lorenz system, Chua, Rssler, Chen systems. In that paper [3] , the author classified systems of Shilnikov’s type according to the type of chaotic attractors, existing in these systems: 1) chaos with homoclinic orbit; 2) chaos with heteroclinic orbit; 3) chaos of mixed type with homoclinic and heteroclinic orbits; 4) other types of chaos. The main aim of this work is to present the research of several model systems of ODEs of Shilnikov’s type, mentioned in the above article, conducted in accordance with the trajectory approach of the universal bifurcation theory of Feigenbaum- Sharkovskii-Magnitskii (FSM) theory [1] [2] ; the point is to confirm the existence of the only type of dynamic chaos and universal scenario of transition to chaos. This scenario begins with the Feigenbaum period doubling cascade of bifurcations of some original stable cycle, then continues with the Sharkovskii complete or incomplete subharmonic cascade of bifurcations of stable cycles of arbitrary period up to the cycle of period three and the Magnitskii complete or incomplete homoclinic cascade of bifurcations of stable cycles converging to homoclinic contours of equilibrium points or cycles.
2. Results and Discussion
In this work several chaotic systems from paper [3] are analysed. We use the fourth-order Runge-Kutta method to approximate the solutions of nonlinear autonomous systems of ODEs. Numerical simulations are prefaced by the research of system’s equilibrium points, their type and stability.
2.1. The Sprott(e) System
The first system we study in this paper is a chaotic system Sprott(e) introduced by J. C. Sprott [4] .
(1)
where a is a real parameter (Figure 1).
2.1.1. Analytical Analysis
The system (1) is dissipative everywhere in phase space since
(2)
Position of equilibrium points is defined by conditions:
(3)
,
,
. Thus, the system (1) has one equilibrium point
. Let us find the roots of characteristic polynomial of Jacobi matrix in order to determine its type. We have the following form of linearization matrix of the right-hand side of the system (1):
(4)
![]()
Figure 1. One of the irreglular attractors of the system for
.
For
it reads the form:
(5)
The characteristic equation for the Jacobi matrix at point O looks as follows:
(6)
Thus, we have
,
as roots of the characteristic equation, consequently our equilibrium point O of the system (1) is a center.
2.1.2. Numerical Simulations
The cascade of bifurcations occurs with decreasing a. For
the phase trajectory of the system’s solution appears as a stable cycle, at the value
a period doubling bifurcation of the original stable cycle takes place and for
a double period cycle is observed, for
we found a cycle of period 4, for
we found a cycle of period 8, etc. The Feigenbaum attractor was found at
. Thus, in this system the Feigenbaum period doubling cascade of bifurcations can be observed (Figure 2).
Furthermore, at
a cycle of period 9 was found, for
―a cycle of period 7, for
-a cycle of period 5 (Figure 3). That means, that the scenario of transition to chaos continues with the subharmonic cascade of bifurcations of stable cycles of arbitrary period in accordance with the Sharkovskii’s order. For
a cycle of period 3 was found, its doubling for
and tripling for
. That evidences the existence of the Feigenbaum cascade and complete subharmonic cascade for the cycle of period 3 (Figure 4). No cycles of homoclinic cascade of bifurcations were observed with further decreasing of the parameter a.
2.2. The Sprott(n) System
We consider the other model system-Sprott (n) system [4] .
(7)
where a is a real parameter (Figure 5).
2.2.1. Analytical Analysis
Divergence of the system (7) is negative everywhere in phase space:
(8)
The system has one equilibrium point
.
The Jacobi matrix of this system at point O is following:
(9)
We shall find the eigenvalues of
from the characteristic equation
(10)
![]()
Figure 5. One of the chaotic attractors for
.
Let us take
. We have:
(11)
where
,
.
The equation (11) has one real and two complex conjugate roots if the discriminant
. The discriminant (11) is determined by
(12)
We have:
(13)
Thus, when (13) holds, the original characteristic equation (10) also has one real and two complex conjugate roots.
An equilibrium point is a saddle-focus if the real root is negative and the real part of complex conjugate roots is positive. Let us determine these conditions.
Notice that
, where
are the roots of (10). Consequently, the inequality
reads as
,
reads as
.
We have:
(14)
Thus, we get that under the condition (13) for the existence of one real and two complex conjugate roots and the condition (14) for negativity of a real root and the positivity of complex conjugate roots, the equilibrium point
is an unstable saddle-focus.
2.2.2. Numerical Simulations
For values of the bifurcation parameter
there is a stable cycle in the system. At
the Feigenbaum cascade begins with birth of a double period cycle which exists for
, then a cycle of period 4 for
, a cycle of period 8 for
emerge. At
the period doubling cascade comes to the end with the formation of a Feigenbaum attractor (Figure 6).
Then the scenario of transition to chaos continues with the subharmonic cascade of bifurcations. We found the following cycles: for
a cycle of period 9, for
a cycle of period 7, for
a cycle of period 5, for
a cycle of period 3 (Figure 7). We also found doubling of period 3 cycle for
and its tripling for
, confirming the presence of the complete subharmonic cascade for a cycle of period 3 (Figure 8). For
a homoclinic cycle of period 4 was found (Figure 9). This implies that the scenario further continues with Magnitskii homoclinic cascade of bifurcations of stable cycle with further increase of the parameter a.
2.3. The Burke-Shaw System
The last system we study in this paper is the Burke-Shaw system derived by B. Burke and R. Shaw from the Lorenz system [5] .
(15)
where a and b are real parameters (Figure 10).
2.3.1. Analytical Analysis
Divergence of the system (15) is negative in the area
![]()
Figure 9. The homoclinic cycle of period 4 for
and its doubling for
.
(16)
Taking into consideration the condition (16) we have
,
,
. The system has two equilibrium states
and
, which are symmetric with respect to the z axis.
Eigenvalues of a matrix of linearization are defined at these points by the characteristic equation:
(17)
Furthermore, we shall notice that if
(18)
then
(19)
and
(20)
The characteristic polynomials of linearization matrix at the equilibrium points
and
each have one negative real root and a pair of complex conjugate roots with positive real part, therefore, under these conditions the equilibrium points of the system are unstable saddle-focuses.
2.3.2. Numerical Simutions
We found that shift of dynamics of the system towards more complex solutions (which takes place with the growth of the parameter b) is in agreement with the FSM-scenario. We fix the value of the parameter
.
For
a stable cycle is the only attractor of the system. Next the Feigenbaum period doubling cascade of the original cycle starts: for
a double period cycle, for
a cycle of period 4, at
the Feigenbaum attractor (Figure 11).
Several cycles from Sharkovskii cascade were observed: a cycle of period 9 for
, a cycle of period 7
, a cycle of period 5 for
, a cycle of period 3 for
(Figure 12).
With further increase of the bifurcation parameter we found a cycle of the homoclinic cascade with 4 loops at
and a cycle of the heteroclinic cascade at
(Figure 13).
Thus, the scenario of transition to chaos in the system (15) also occurs in full accordance with the FSM-theory.
![]()
Figure 13. The homoclinic and heteroclinic cycles.
3. Conclusion
The present paper carries out analytical and numerical study of three chaotic systems of Shilnikov’s type referred in Chen’s article [3] , two Sprott’s systems―the Sprott(e) and the Sprott(n) and the Burke-Shaw system. For these systems, bifurcation parameters and types of fixed points were found, as well as the occurrence of initial cycles of the bifurcation cascades and the occurrence of chaotic attractors. With variation of the parameters, we studied phase portraits in this three systems and revealed the same scenario of singular attractors occurrence, which fully corresponded to the mechanism of transition to chaos by Feigenbaum-Sharkovskii-Magnitskii. So, despite classification [3] , all systems seem to fall into one category according to the FSM theory, as chaotic attractors have the same nature.
Acknowledgements
The author is grateful to Prof. N.A. Magnitskii for critical comments concerning analysis of results presented in this paper, to O.I. Ryabkov and D.A. Burov for valuable discussions.