Received 16 March 2016; accepted 20 May 2016; published 23 May 2016
In recent years, many scientists around the world have tried to analyze and classify the variety of systems with chaotic behavior. In the article  , 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  , 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   ; 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  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  .
where a is a real parameter (Figure 1).
2.1.1. Analytical Analysis
The system (1) is dissipative everywhere in phase space since
Position of equilibrium points is defined by conditions:
, ,. 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):
Figure 1. One of the irreglular attractors of the system for.
For it reads the form:
The characteristic equation for the Jacobi matrix at point O looks as follows:
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.
Figure 2. Cycles of the Feigenbaum cascade accordingly for, , , and the Feigenbaum attractor for.
Figure 3. Cycles of periods 9, 7, 5 of the Sharkovskii subharmonic cascade for the parameter values, ,.
Figure 4. The cycle of period 3 for, its doubling for and its tripling for.
2.2. The Sprott(n) System
We consider the other model system-Sprott (n) system  .
where a is a real parameter (Figure 5).
2.2.1. Analytical Analysis
Divergence of the system (7) is negative everywhere in phase space:
The system has one equilibrium point.
The Jacobi matrix of this system at point O is following:
We shall find the eigenvalues of from the characteristic equation
Figure 5. One of the chaotic attractors for.
Let us take. We have:
The equation (11) has one real and two complex conjugate roots if the discriminant. The discriminant (11) is determined by
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.
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  .
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 6. Cycles of period doubling cascade of bifurcations accordingly for the parameter values, , , ,.
Figure 7. Cycles of the periods 9, 7, 5 for, , correspondingly.
Figure 8. The cycle of period 3 and its doubling and tripling for, and in corresponding order.
Figure 9. The homoclinic cycle of period 4 for and its doubling for.
Taking into consideration the condition (16) we have, ,. The system has two equilibrium states and, which are symmetric with respect to the z axis.
Figure 10. The chaotic attractor for,.
Eigenvalues of a matrix of linearization are defined at these points by the characteristic equation:
Furthermore, we shall notice that if
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
Figure 11. Cycles of period doubling cascade of bifurcations for the parameter values, , , correspondingly.
Figure 12. Cycles of the Sharkovskii subharmonic cascade of bifurcations for, , , in corresponding order.
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.
The present paper carries out analytical and numerical study of three chaotic systems of Shilnikov’s type referred in Chen’s article  , 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  , all systems seem to fall into one category according to the FSM theory, as chaotic attractors have the same nature.
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.