Next Article in Journal
Mixed-Integer Linear Programming Model and Heuristic for Short-Term Scheduling of Pressing Process in Multi-Layer Printed Circuit Board Manufacturing
Next Article in Special Issue
A Combined Energy Method for Flutter Instability Analysis of Weakly Damped Panels in Supersonic Airflow
Previous Article in Journal
Multidimensional Fair Fuzzy Equilibrium Evaluation of Housing Expropriation Compensation from the Perspective of Behavioral Preference: A Case Study from China
Previous Article in Special Issue
Performance Analysis on the Use of Oscillating Water Column in Barge-Based Floating Offshore Wind Turbines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hidden Strange Nonchaotic Attractors

by
Marius-F. Danca
1 and
Nikolay Kuznetsov
2,3,*
1
Romanian Romanian Institute of Science and Technology, 400487 Cluj-Napoca, Romania
2
Department of Applied Cybernetics, Saint-Petersburg State University, Peterhof, 198504 Saint-Petersburg, Russia
3
Department of Mathematical Information Technology, University of Jyväskylä, 40014 Jyväskylä, Finland
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(6), 652; https://doi.org/10.3390/math9060652
Submission received: 19 December 2020 / Revised: 4 February 2021 / Accepted: 10 March 2021 / Published: 18 March 2021
(This article belongs to the Special Issue Dynamical Systems and Their Applications Methods)

Abstract

:
In this paper, it is found numerically that the previously found hidden chaotic attractors of the Rabinovich–Fabrikant system actually present the characteristics of strange nonchaotic attractors. For a range of the bifurcation parameter, the hidden attractor is manifestly fractal with aperiodic dynamics, and even the finite-time largest Lyapunov exponent, a measure of trajectory separation with nearby initial conditions, is negative. To verify these characteristics numerically, the finite-time Lyapunov exponents, ‘0-1’ test, power spectra density, and recurrence plot are used. Beside the considered hidden strange nonchaotic attractor, a self-excited chaotic attractor and a quasiperiodic attractor of the Rabinovich–Fabrikant system are comparatively analyzed.

1. Introduction

In the last three decades, strange nonchaotic attractors (SNAs) have attracted much attention from both theoretical and experimental points of view, reportedly arising in many physically relevant situations. Though somewhat exotic, SNAs are not at all rare. In fact, SNAs have been reported in quasiperiodically forced pendulum, quantum particles, quasiperiodic potentials, biological oscillators, quasiperiodically driven Duffing-type oscillators, velocity dependent oscillators, electronic circuits, etc. (see, e.g., [1] and references therein). Additionally, some models are studied in higher dimensions and for both discrete and continuous systems. Commonly existing in quasiperiodically driven nonlinear systems, SNAs are found and studied for non autonomous quasiperiodically forced dynamical systems and can be regarded as “intermediate” structures linking regularity and chaos. However, as this paper shows, there are undriven continuous systems with SNAs, or even with hidden SNAs.
SNAs were introduced in 1984, when Grebogi, Ott, Pelikan, and Yorke constructed dynamical systems with attractors that were strange but not chaotic [2] (see also [3,4,5,6,7]). Results on flows containing SNA were obtained far before than the term SNA was coined (see [8], where some references (in Russian) can be found)). The first experimental observation of an SNA was in a magnetoelastic ribbon [9] (see also [10,11,12]). The ‘strange’ characteristic refers to the fractal, nontrivial, and complicated geometry of the attractor. There are authors who consider SNAs to be “strange” because they are not piecewise differentiable. The relatively new term, SNAs, has still not been precisely formulated and their existence has not been rigorously proven. Therefore, most results on SNAs, such as those in this paper, are based on numerical analysis.
The ‘nonchaotic’ property means that SNAs do not show sensitivity to initial conditions—i.e., with negative Lyapunov exponents (LEs), just like regular systems. Note that some authors consider the nonchaotic property to be defined by nonpositive LEs (see, e.g., [6]). Although it is difficult to distinguish an SNA visually from a strange chaotic attractor, dynamically there are important differences. Indeed, orbits on SNAs and chaotic attractors are non-periodic, but compared to chaotic attractors the trajectories on SNAs are not separate from each other. The dynamical behavior characterizing SNAs occurs in a very narrow range of the control parameter values.
In principle, strange nonchaotic attractors occur in all dissipative dynamical systems that exhibit the period-doubling route to chaos, where the attractors formed at the accumulation points are fractal sets with zero Lyapunov exponents. Such attractors are, however, not physically observable because the set of parameter values for them to arise has a Lebesgue measure of zero in the parameter space [13].
SNAs are typically not robust in the absence of quasiperiodic driving. Therefore, they are not expected to occur naturally, whilst, as shown by Grebogi et al. in [2], if the system is quasiperiodically driven then SNAs can be robust. A robust SNA is typical but not vice versa.
From a computational point of view, attractors can be classified as self-excited (SE) and hidden. While SE attractors can be localized numerically with some standard computational procedures by starting from a point in a neighborhood of unstable equilibrium, the basins of hidden attractors, which can be chaotic (HCAs) or non-chaotic (e.g., stable cycles), are not connected with equilibria. Thus, if the trajectories are attracted by some attractor, the attractor is an SE attractor, while a hidden attractor cannot be excited by the unstable equilibria (if they exist) and trajectories from neighborhoods of unstable equilibria are attracted either by some attractors (stable equilibria or cycles) or diverge to infinity [14,15,16,17,18,19]. Therefore, the HCAs are usually found empirically. In other words, it means to verify if the attraction basins of the HCA do not intersect small neighborhoods of unstable equilibria. However, finding the attraction basins of HCAs analytically remains an open and challenging problem.
The investigation of hidden oscillations actually arose in the second part of Hilbert’s 16th problem and in Bautin’s works, related to nested limit cycles in quadratic systems, which showed the necessity of studying hidden oscillations to solve this problem [15]. HCAs may appear in, e.g., chaotic systems with one or several stable equilibriums or with a line equilibria (infinitely many equilibria), or without equilibrium [20,21,22,23].
In this paper, attractors and transients of the Rabinovich–Fabrikant (RF) are differentiated via their lifetime and it is shown numerically that some HCAs of the RF system present the characteristics of SNAs. For this purpose, the spectrum Λ of the finite-time local Lyapunov exponents (LEs), the ‘0-1’ test (Appendix A), the two-sided Power Spectral Density (PSD) determined for one component, the Correlation Dimension ( D 2 ), and the Recurrence Plot (RP) (Appendix B) are used. Additionally, histograms and the difference between two close trajectories are utilized to differentiate SNAs from chaotic attractors.
This paper is organized as follows: Section 2 presents the Rabinovich–Fabrikant system with some representative attractors and Section 3 shows that the system has hidden strange nonchaotic attractors (HSNAs). The Appendix A, Appendix B, Appendix C present the notions of the ‘0-1’ test and the recurrence plot. Finally, the conclusion in Section 4 ends the paper.

2. Hidden Attractors in RF System

The RF system [24,25], arising from the modulation instability in a nonequilibrium dissipative medium, is modeled by:
x . 1 = x 2 x 3 1 + x 1 2 + a x 1 , x . 2 = x 1 3 x 3 + 1 x 1 2 + a x 2 , x . 3 = 2 x 3 b + x 1 x 2 ,
where a and b are some real positive parameters and b is the bifurcation parameter.
The system has extremely rich dynamics, presenting coexisting attractors, SE chaotic attractors (SECAs), and HCAs [26]. The equilibria of the system are:
X 1 , 2 * b R 1 + 2 b 4 b 3 a , ± b 4 b 3 a R 1 + 2 , a R 1 + R 2 4 b 3 a R 1 + 8 b 6 a , X 3 , 4 * b R 1 2 b 3 a 4 b , ± b 4 b 3 a 2 R 1 , a R 1 R 2 4 b 3 a R 1 8 b + 6 a ,
and X 0 * ( 0 , 0 , 0 ) , with:
R 1 = 3 a 2 4 a b + 4 and R 2 = 4 a b 2 7 a 2 b + 3 a 3 + 2 a .
Due to the general instability of chaotic solutions and due to the particular sensibility of the RF system to numerical integration, the dedicated numerical methods implemented in known software cannot be used here with 100% confidence, especially for long-time integration, because global errors increase with the size of the integration interval [27,28]. To numerically integrate the RF system, in this paper Matlab solver ode45 is used. To note that the system presents some stiffness characteristics (the x 3 component of the solution decreases faster than x 1 , 2 [29]). However, the matlab solvers for stiff systems give, in this case, similar results). Most importantly, the numerical results cannot be improved by decreasing the integration step, since the integration error has an extremum as a function of the integration step-size [30]. Another approach would be higher-accuracy calculations, but this approach has restrictions: on one hand, the way to decrease the error is narrow; on another hand, the number of operations needed, if a very small integration step-size is used, could increase dramatically (see [30] and the references therein). Additionally, note that it is possible to have reliable numerical simulations only in a relatively small time interval [31,32,33]. For example, for the classical Lorenz system, it is a real challenge to get solutions with small errors (in the sense that the largest error between two computed solutions for small but different time-steps has a tolerance of 5 × 10 14 ) for t [ 0 , 100 ] [31].
Therefore, one can never be confident if an apparently chaotic numerical trajectory on larger time intervals represents a real chaotic trajectory or just a chaotic transient. Beside these obstacles in the numerical integration of ODEs of the RF system, the first intensive numerical investigations [25] revealed an extremely strong dependence of the solution on the utilized numerical methods, initial conditions, and the parameter b.
On the other hand, transient chaos appears when trajectories starting from some random initial conditions are chosen near the boundary of the basin of attraction (boundary crisis) or are very close to some bifurcation point. The underlying trajectories behave chaotically, possibly for a relative long time, and then quite abruptly switch to some final periodic state [34,35,36]. Because the lifetime of the transients can be very long (thousands and tens of thousands of numerical iterations), one can never be 100 % confident that an observed chaotic system is not a transient (Sprott’s communication). Therefore, to differentiate transients from attractors of the RF system numerically, a compromise with a sufficiently large time interval is made, with [ 0 , T ] and T = 10 , 000 . For a = 0.1 , b [ 0.2808 , 0.2907 ] and with the initial condition x 0 = [ 0.02 , 0.12 , 0.53 ] , Figure 1 shows the lifetime of the trajectories together with the variation in the spectrum of the LEs, Λ = { λ 1 , λ 2 , λ 3 } , with λ 1 λ 2 > λ 3 . In order to increase T over T = 10 , 000 , o p t i o n s = o d e s e t ( R e l T o l , 1 e 6 , A b s T o l , 1 e 9 ) is used [37]. Note that without using o p t i o n s , some of attractors of the RF system transform into chaotic transients. Thus, in this paper the trajectories which reach a final regular motion before T = 10 , 000 (here the stable equilibria X 1 , 2 * ) with an error of less than 1 × 10 5 are considered transients. The system presents two kinds of transients: spirals which lead to X 1 , 2 * and have negative LEs (black line segments, Figure 1a), and chaotic transients with positive maximal LE (light green line segments, Figure 1a). The spectrum of the local finite-time LEs, Λ = { λ 1 , λ 2 , λ 3 } (Figure 1b), is determined numerically from the system equations with the Wolf algorithm.
Considering the lifetime classification mentioned before, irregular or regular trajectories, which continue their dynamics beyond T = 10 , 000 , are considered in this paper attractors as follows: periodic attractors (see red vertical lines in Figure 1a; stable cycles of the RF system, are not considered in this paper, see [24,25]), quasiperiodic attractors (a quasiperiodic attractor is presented in Figure 2b), hidden chaotic attractors (blue line segments, Figure 1a), and hidden strange nonchaotic attractors (dark green line segments, Figure 1a). Quasiperiodic attractors and irregular transients are merged in the parameter space. There is little distinction between SNAs (in this paper, the attractor H corresponding to b = 0.2875 ) and chaotic attractors (here, the HCA corresponding to b = 0.28708 ), as can be seen in the top of Figure 1. They superficially look quite similar, with their structure remaining actually almost unchanged, and even the maximum LEs have the opposite sign. However, as will be discussed next, there are dynamically significant distinctions (see, e.g., [6]).
Equilibria X 1 , 2 * are stable (stable focus node) for the chosen parameter values, and orbits starting in close neighborhoods will be attracted by X 1 , 2 * .
Equilibria X 3 , 4 * are unstable (attracting saddles) and orbits starting from close neighborhoods or the attraction basin of X 3 , 4 * will be attracted by X 3 , 4 * on the surface determined by the two-dimensional stable manifold; then, after some time, they will exit along the direction of the 1-dimensional unstable manifold.
X 0 * is unstable (repelling focus saddle) for all values of a and b, being globally asymptotically unstable [24], and is a global attractor for the reduced system on x 3 = 0 . Therefore, trajectories starting in the plane x 3 = 0 within neighborhoods of X 0 * will diverge to infinity, while trajectories in neighborhoods of X 0 * with x 3 > 0 are attracted by scrolling and by the stable plane of unstable equilibria X 3 * or X 4 * , after which they are repulsed by the unstable directions to infinity. Additionally, depending the initial conditions on the neighborhood of X 0 * , they can be attracted by X 1 , 2 * along the stable one-dimensional manifold of these equilibria.
Because all equilibria have a pair of complex eigenvalues, all trajectories starting within neighborhoods of equilibria present a scrolling dynamic. Details on equilibria X 0 * and X 1 4 * can be found in [24,25].
Consider next, for a = 0.1 , the attractors corresponding to b = 0.98 , with the initial condition ( 0.158 , 6.188 , 0.977 ) (Figure 2a); b = 0.25675 , with the initial condition ( 0.1 , 0.1 , 0.1 ) , without transients (Figure 2b); and b = 0.2875 , with the initial condition ( 0.2 , 0.38 , 0.001 ) (Figure 2c).
The chaotic attractor corresponding to b = 0.98 (Figure 2a) is SECA, because it can be found starting from initial conditions within neighborhoods of unstable equilibria X 3 , 4 * . To characterize this attractor and trace the chaos, consider Figure 3, where the spectrum of LEs (Figure 3a); the median K (Appendix A) given by the test ‘0-1’ (Figure 3b); the PSD (Figure 3c); the q as a function of p, given by the ‘0-1’ test (Figure 3d); the mean-square displacement M as function of t (Figure 3e); and the RP in Figure 3f (Appendix B) are presented, respectively. The positiveness of the maximal LE, λ 1 = 0.062 , indicates that the attractor is chaotic. This fact is underlined also by the median K 0.98 , close to 1. The broadband of the PSD indicates that the system is chaotic for this value of b. Moreover, the behavior of q as a function of p and the mean-square displacement M are typical for chaotic dynamics (see Appendix A). The RP, with ϵ = 0.95 and t [ 0 , 1000 ] , presents intricate structures that are typical for chaos.
Right after the critical point b = 0.2567 , when λ 1 = λ 2 = 0 and λ 3 = 0.314 (Figure 4a), for b = 0.25675 the obtained attractor is quasiperiodic (Figure 2b). For b [ 0.2567 , 0.2688 ] , the variation in Λ seems to unveil more bifurcations. However, because of the difficulty of analytical study, the bifurcations are not considered here. As for regular motions, K is close to 0, K 0.001 (Figure 4b). The quasiperiodicity is revealed by the PSD (Figure 4c), whose peaks represent the sub-harmonics (multiples) of a main harmonic generated at b = 0.2567 . q as a function of p (Figure 4d) and the evolution of the mean-square displacement M (Figure 4e) show that the attractor is quasiperiodic. Besides the continuous equidistant and parallel off-second diagonal line in the recurrence plot, with ϵ = 0.95 and t [ 0 , 1000 ] (Figure 4f), there exist another set of interrupted short off-second diagonal segments, parallel with the first set of lines (see the larger circled zoom in Figure 4f). Since all these lines are non-equidistant, the dynamics are quasiperiodic. Note that another set of interrupted off-second diagonal lines seem to be born (see the smaller circled zoom in Figure 4f), suggesting a supplementary frequency. However, this frequency cannot be seen in the PSD (Figure 4c).
Consider the attractor corresponding to b = b * = 0.2875 , denoted with H , within a neighborhood of b * , V b * = [ 0.28743 , 0.2876 ] (see Figure 1), obtained with the initial condition x 0 = ( 0.2 , 0.38 , 0.001 ) [37]. Via the mentioned numerical tools, it is shown that the attractor presents the characteristic of SNAs. The attractor H is hidden because its attraction basin has no connection with neighborhoods of unstable equilibria (see [26] for details) and trajectories starting from small neighborhoods of unstable equilibria X 3 , 4 and X 0 * are attracted either by the stable equilibria X 1 , 2 or diverge to infinity, as shown in Figure 2c. For clarity, in Figure 2c only one trajectory from each equilibrium is presented. As can be seen, trajectories starting from a neighborhood of the unstable equilibrium X 0 * tend either to X 1 * (orange plot), to X 2 * (light green plot), or to infinity via the unstable equilibrium X 3 * (black plot) or via the unstable equilibrium X 4 * (grey plot). Trajectories from the neighborhood of X 3 * are attracted either by X 2 * (red plot) or by infinity (brown plot). Trajectories starting from the neighborhood of X 4 * are attracted either by X 1 * (blue plot) or by infinity (brown plot). The attraction basin of H within a parallelipipedic neighborhood of the unstable equilibrium X 0 * is plotted by small red spheres in Figure 5a (note that the system is defined for x 3 > 0 and the axis x 3 is invariant with the reduced equation x 3 ˙ = 2 b x 3 , which has the solution x 3 ( t ) = e 2 b t x 2 ( 0 ) [24]). In Figure 5b, one can see a zoomed cubic area around the unstable equilibrium X 0 * , which underlines the fact that, for the considered resolution, there exists an empty cone-like neighborhood with a radius of about 0.05 , without initial conditions (red points) of H , indicating that the attraction basin of H has no connection with small neighborhoods of the unstable equilibrium X 0 * . From Figure 5c, one can see that beside the red points (initial conditions) leading to the hidden irregular attractor H , there are two other kinds of initial points which lead to the attractive equilibria X 1 , 2 * (blue and green plots) [37].

3. SNAs of the RF System

Next, it is verified that the hidden attractor H verifies the characteristics of SNAs. For this purpose, the ‘0-1’ test, PSD, and RP are used. In Table 1, the characteristics of the considered attractors of the RF system are presented comparatively.
As shown in Figure 6, the LEs are negative for a small but connected interval of b, V b * (Figure 6a); the oscillations of LEs outside V b * indicate presumably attractors’ coexistence; the median K has an intermediate value between 0 and 1, K = 0.62 (Figure 6b); the few peaks in the PSD (F in Figure 6c) not covered by the broadband in the PSD indicate that the attractor is neither quasiperiodic nor chaotic. To note that in [38] it was shown that the PSD of SNAs is discrete, but very dense. Additionally, q as a function of p has a graph between a perfect disc (or circle), for the case of regular behavior, and the Brownian motion (Figure 6d), for chaos; the mean-square displacement M presents an oscillating but bounded increase (Figure 6e); the RP, with ϵ = 0.95 and t [ 0 , 500 ] , presents many non-equidistant off-second diagonal lines parallel with the second diagonal beside several isolated points.
Moreover, the strangeness character, which refers to the nontrivial complicated fractal geometry of the attractor, can be seen from, e.g., the dotted plot in the zoomed area D in Figure 2c, which, besides the non-differentiability of the trajectory, reveals the repeating isolated dust-like points typical for fractals and the correlation dimension D 2 (Appendix C) is D 2 = 1.2 .
Since the system is time-continuous, the time series are generally oversampled. This happens because most numerical integrators have variable step sizes, which can generate extremely large numbers of points within some critical trajectory zones. Therefore, before applying the ‘0-1’ test, the underlying time series has to be resampled. Otherwise, it is possible to obtain small (possible zero) values of K, even if the underlying attractor is chaotic [39]. A simple way to resample in Matlab is to use the function resample.
Because of the negativity of the maximal LE on the HSNA H , trajectories do not diverge by separating from each other and eventually coincide, as presented in Figure 7. Comparatively, the first components x 1 (red plot) and y 1 (blue plot) of two trajectories x and y, starting from two nearby initial conditions x 0 and x 0 + ( 1 × 10 5 , 1 × 10 5 , 0 ) for the SECA corresponding to b = 0.98 (Figure 7a), the quasiperiodic attractor corresponding to b = 0.263 (Figure 7b), and the HSNA H (Figure 7c) are considered. In the case of coupled SNAs, this property characterizes a robust synchronization, in the sense that the non-dependence on initial conditions is expected to synchronize the SNAs quickly, even when they operate with different sets of initial conditions [40].
Figure 7a(i)–c(i) shows overplots of x 1 and y 1 . Details D 1 , D 2 , and D 3 clearly underline the differences between the chaotic attractor, on one side, and the quasiperiodic attractor and H on the other side. Figure 7a(ii)–c(ii) present details of Figure 7a(i)–c(i). Details D 4 , D 5 , and D 6 of the three-dimensional representations are obtained after removing the first 7500 transient points. Meanwhile, Figure 7a shows that the difference between the two chaotic trajectories separate—see Figure 7b,c, with details D 2 and D 3 —suggesting that both quasiperiodic trajectories and both trajectories in the case of HSNA present periodic evolution while remaining close to each other. Note that this property of the H is not so evident as for SNAs in quasiperiodically driven discrete systems.
Because the RPs do not show the differences between the considered attractors clearly enough, in Figure 8 the histograms of the considered attractors are considered. As can be seen, the HSNA H can be viewed as being situated between the quasiperiodic attractors and the chaotic attractors. Hidden SNAs, such as H , also exist for b ( 0.28743 , 0.2876 ) (see Figure 1).

4. Conclusions and Discussion

In this paper, it has been shown numerically that the attractor corresponding to b = 0.2875 presents the characteristics of SNAs. For this purpose, the following tools were utilized: Lyapunov exponents, the ‘0-1’ test, power spectral density, recurrence plots, histograms, and the difference between close trajectories. The considered attractor H has all negative LEs, the median K value is between 0 and 1 ( K 0.62 ), the graph of p and q presents a disk-like structure, the mean-square displacement M has an oscillatory bounded growth, and the recurrence plots present non-equidistant off-second diagonal lines parallel to the second diagonal and also several isolated points. Additionally, contrarily to chaotic attractors, it is verified numerically that the distance between two close trajectories on the HSNA H does not change significantly over time. The fact that LEs look almost similar for the HSNAs and regular attractors could be explained by some inaccuracy in LE calculations for the multi-stable regime observed here. On the other hand, differences between the HSNA and regular attractors of this system are clearly underlined by the 0-1 test, recurrence plot, and power spectral density, which are different.
Even in dynamics do not have explicit quasiperiodic forcing, it is shown numerically that the motion of the RF system could appear on strange nonchaotic attractors. Note that discrete dynamical systems, without a quasiperiodic driving force where the motion is on strange nonchaotic attractors, are considered in [41].
Some open problems deserve further analysis. Thus, the obtained results in this paper regarding the SNA characteristics of H should be considered subject to possible numerical errors due to system’s strong dependence on the numerical integrators and time integration interval [ 0 , T ] , here set to T = 10,000 by using the Matlab o p t i o n s function; strong sensibility on the parameters (a perturbation of order 1 × 10 4 can move the trajectory to different attractors); initial conditions; etc. Probably due to these reasons, the conclusions regarding the SNA characteristics of H are not so evident, just like the cases of quasiperiodic-driven systems. However, the existence of this hidden strange nonchaotic attractor is verified numerically for a connected parameter interval. Therefore, the attractor H would not be just simply a numerical “accident”.
Using different numerical methods to integrate the system, the results could be somewhat different but consistent: there exist small parameter intervals where the RF system has HSNAs.
The usual transition q u a s i p e r i o d i c i t y S N A c h a o s does not seem to be characteristic for this system. Indeed, Figure 1 shows that after the interval of b, which generates HSNAs, it follows a parameter window where, after chaotic (with maximal positive LE) or nonchaotic (with negative maximal LE) transients, the system enters into a regular dynamics domain (stable equilibria X 1 , 2 * ). Note that in quasiperiodically forced dynamical systems, the values of the Lyapunov exponents are negative for both quasiperiodic attractors (tori) and SNAs, and in this case the Lyapunov exponents fail to detect transitions from quasiperiodic dynamics to SNAs [42].
On the other hand, an open problem here is the presumably weak chaos. Weak chaos means that even though two trajectories are initially close, they will diverge to be weaker, and their divergence remains finally bounded. In other words, trajectories are not simply regular nor fully chaotic. The separation of nearby trajectories in the case of weak chaos implies that the corresponding maximal Lyapunov exponent is zero [43,44]. Sometimes, this situation appears for some large parameter intervals (see, e.g., the light green vertical bands in Figure 9a). Therefore, the reliability of the ‘0-1‘ test [45], used for this system, could be interpreted from the perspective of the weak chaos. For b [ 0.255 , b 1 ] [ b 3 , b 4 ] (Figure 9), because the maximal LE, λ 1 , is 0 (Figure 9a), the oscillations of K (Figure 9b or Figure 4b) could be due to the existence of weak chaos. For b [ b 4 , b 5 ] , K increases but cannot reach the value corresponding to chaos, 1, and the maximal LE is only slightly larger than 0. A clear image for b [ b 5 , 0.288 ] can be seen in Figure 6a,b. Actually, even the advantage of the ‘0-1’ test is due to the fact that it is binary: ‘white’ (0) and ‘black’ (2). The ‘gray’-like situations, such as values ’close’ to 1 or 0 for the RF system, require a careful supplementary analysis to clearly state that H is non-chaotic.
Additionally, as Figure 7 show, even when two close trajectories remain at a constant distance to each other, the difference between them seems not to vanish at the end of a relatively large integration interval, unlike in the case of SNAs in quasiperiodically driven discrete systems, where, under some circumstances, trajectories eventually coincide, which represents a characteristic of SNAs: robust synchronization.
The use of the recurrence quantification analysis [46] for the detection of structural changes in the dynamics of complex nonlinear systems could give a supplementary clarification of H .
Searching SNAs in other continuous time systems could offer more exciting results.

Author Contributions

Conceptualization, M.-F.D.; methodology, M.-F.D. and N.K.; software, M.-F.D.; validation M.-F.D. and N.K.; investigation, M.-F.D. and N.K. All authors have read and agreed to the published version of the manuscript.

Funding

The work is supported by the Russian Science Foundation 19-41-02002 and St. Petersburg State University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from authors.

Acknowledgments

Figure 5 was prepared with help from Paul Bourke.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The ‘0-1’ Test

The ’0-1’ test, which has its roots in [47] and was developed in [48] (see also [39] or [42]), is designed to distinguish chaotic behavior from regular behavior in deterministic systems, which could be used to identify SNAs as well. However, the ‘0-1’ test has been used to identify some SNAs (see, e.g., [42]). The test is easy to implement and does not need system equations, with the input being a time series. Consider a dynamical system with a one-dimensional observable data set, constructed from a time series of the underlying system, ϕ ( j ) , j = 1 , 2 , . . . , N . The ‘0-1’ test is based on a theorem, which states that a nonchaotic motion is bounded, while a chaotic dynamic behaves like a Brownian motion [47].
(1)
First, for c [ 0 , 2 π ] , compute the translation variables p and q by [48]:
p ( n ) = j = 1 n ϕ ( j ) cos ( j c ) , q ( n ) = j = 1 n ϕ ( j ) sin ( j c ) ,
for n = 1 , 2 , , N . The choice of c represents an important and sensible algorithm variable (see, e.g., [39], where c is recommended on the interval [ π / 5 , 4 π / 5 ] ).
(2)
To determine the growths of p and q, here we used the mean-square displacement M:
M ( n ) = lim N 1 N j = 1 N [ p ( j + n ) p ( j ) ] 2 + [ q ( j + n ) q ( j ) ] 2 .
where n N (in practice, n = N / 10 would be a good choice).
(3)
The asymptotic growth rate K is defined as follows:
K c = lim n log M ( n ) / log n .
Because of the occurrence of resonances for some isolated values of c (where K c could be larger), the median of the computed values of K c (generally, 100 values), denoted by K, is used, since the median is robust against outliers associated with resonances [48]. If the underlying dynamics are regular (i.e., periodic or quasiperiodic), then K 0 ; if they are chaotic, then K 1 . For SNAs, K takes an intermediate value between 0 and 1.
In Figure A1a, for the cases of the logistic map x n + 1 = r x n ( 1 x n ) and the GOPY map [2] x n + 1 = 2 a tanh ( x n ) cos ( 2 π θ n ) , θ n + 1 = θ n + ω , the plots of q versus p are shown, while in Figure A1b M is presented as a function of n. Figure A1a(i),b(i) represent the regular orbit of the logistic map for r = 3.55 ; Figure A1a(ii),b(ii) present the chaotic orbit for r = 4 , while Figure A1a(iii),b(iii) present the SNA of the GOPY map with a = 1.5 and ω = ( 5 1 ) / 2 .
Figure A1. The elements of the ‘0-1‘ test; (a) presents the plots of q versus p and (b) presents the mean-square displacement M as a function of n. (i) represents the regular orbit of the logistic map x n + 1 = r x n ( 1 x n ) for r = 3.55 ; (ii) presents the chaotic orbit of the logistic map for r = 4 ; (iii) presents an SNA of the GOPY map x n + 1 = 2 a tanh ( x n ) cos ( 2 π θ n ) , θ n + 1 = θ n + ω , with a = 1.5 , and ω = ( 5 1 ) / 2 .
Figure A1. The elements of the ‘0-1‘ test; (a) presents the plots of q versus p and (b) presents the mean-square displacement M as a function of n. (i) represents the regular orbit of the logistic map x n + 1 = r x n ( 1 x n ) for r = 3.55 ; (ii) presents the chaotic orbit of the logistic map for r = 4 ; (iii) presents an SNA of the GOPY map x n + 1 = 2 a tanh ( x n ) cos ( 2 π θ n ) , θ n + 1 = θ n + ω , with a = 1.5 , and ω = ( 5 1 ) / 2 .
Mathematics 09 00652 g0a1

Appendix B. Recurrence Plots

Recurrences in dynamical systems have been analyzed since Poincaré [49], who showed that under certain conditions, orbits of a bounded dynamical system return arbitrarily close to each former points of its own route. However, the time of this return can be arbitrarily long. A tool to visually analyze this recurrence is provided by recurrence plots (RPs), introduced by Eckmann et al. in [50]. Given an autonomous dynamical system and a trajectory { x i } i N , x i R n , the RP is based on the matrix:
R i , j ( ϵ ) = Θ ( ϵ x i x j ) .
Here, ϵ is a predefined threshold, · is a norm (Euclidean or maximum norm), and Θ is the Heaviside function. The ϵ -recurrent (close) points give “1” while separate points give “0”. In a two-dimensional plot, the value is 1 when a point at position j falls into the neighborhood of radius ϵ of point at i, and 0 when the point at j does not enter the neighborhood of points at i; they are depicted as black and white dots, respectively. In this way, a visual representation of the system dynamics is provided. Generally, the threshold ε should be chosen small enough. However, in the case of noise influence, a larger threshold is necessary. At some time moment, the state of the system will recur as close as one wishes ( ϵ threshold) to a former state. RPs can help find transitions. Diagonal lines correspond to trajectories passing in the same region of the phase space at different times. The RPs of periodic trajectories are equidistant parallel to the second diagonal lines (with slope of 45°), indicating some determinism or periodicity (see Figure A2a for the periodic logistic map, with r = 3.4 , and threshold ϵ = 0.05 ). The vertical distances between lines in the RP (perpendicular to the horizontal axis) represents the period. The RP of a quasiperiodic trajectory consists of parallel, off-second diagonal lines with different distances between them, reflecting the existence of different return times. If many intricate structures (clusters) disposed on interrupted parallel lines appear, the RP indicates a chaotic trajectory. Now, the distance between the parallel lines is not constant (see Figure A2b for the chaotic logistic map, with r = 4 , and threshold ϵ = 1.2 ). For details regarding other RP parameters (embedding dimension, time delay, etc.), see [51,52,53]. For the use of the RP to study SNAs in the forced logistic map, see [54].
Figure A2. The recurrence plot applied to the logistic map x n + 1 = r x n ( 1 x n ) ; (a) the recurrence plot of the logistic map with r = 3.4 and threshold ϵ = 0.05 ; (b) the recurrence plot of the logistic map with r = 4 and threshold ϵ = 1.2 .
Figure A2. The recurrence plot applied to the logistic map x n + 1 = r x n ( 1 x n ) ; (a) the recurrence plot of the logistic map with r = 3.4 and threshold ϵ = 0.05 ; (b) the recurrence plot of the logistic map with r = 4 and threshold ϵ = 1.2 .
Mathematics 09 00652 g0a2

Appendix C. Correlation Dimension

Let an d-dimensional embedding space, x 1 , x 2 , . . . , x m a time series and y 1 , y 2 , , y N , with N = m d + 1 , and the delay coordinates y j = ( x j , , x j + d 1 ) , j = 1 , , N . The correlation integral for a given r > 0 is defined as follows:
C ( r ) = 1 k 2 j = 1 k i = 1 , j i k H ( r | x i x j | ) ,
where H is the Heaviside step function and | . | is the Euclidean distance. Here, the H counts the number of pairs ( x i , x j ) , satisfying the condition | x i x j | < r .
If C ( r ) scales like C ( r ) = r D 2 , then:
D 2 = lim r 0 log ( C ( r ) log ( r ) ,
is called the correlation dimension, the straight line of log-log curve of r and C ( r ) .
The G-P algorithm, created by Grassberger and Procaccia [55], is one of the most used algorithms to determine D 2 numerically, and uses the embedding theory together with the idea of phase space reconstruction.

References

  1. Thamilmaran, K.; Senthilkumar, D.V.; Venkatesan, A.; Lakshmanan, M. Experimental realization of strange nonchaotic attractors in a quasiperiodically forced electronic circuit. Phys. Rev. E 2006, 74, 036205. [Google Scholar] [CrossRef] [Green Version]
  2. Grebogi, C.; Ott, E.; Pelikan, S.; Yorke, J.A. Strange Attractors that are not Chaotic. Physica D 1984, 13, 261–268. [Google Scholar] [CrossRef]
  3. Bondeson, A.; Ott, E.; Antonsen, T.M. Quasiperiodically Forced Damped Pendula and Schrödinger Equations with Quasiperiodic Potentials: Implications of Their Equivalence. Phys. Rev. Lett. 1985, 55, 2103. [Google Scholar] [CrossRef] [PubMed]
  4. Broer, H.W.; Takens, F. Mixed spectrum and rotational symmetry. Arch. Rational. Mech. Anal. 1993, 124, 13–42. [Google Scholar] [CrossRef] [Green Version]
  5. Pikovsky, A.S.; Feudel, U. Characterizing strange nonchaotic attractors. Chaos 1995, 5, 253. [Google Scholar] [CrossRef] [PubMed]
  6. Prasad, A.; Negi, S.S.; Ramaswamy, R. Strange nonchaotic attractors. Int. J Bifurc. Chaos 2001, 11, 291–309. [Google Scholar] [CrossRef] [Green Version]
  7. Romeiras, F.J.; Ott, E. Strange Nonchaotic Attractor in a Nonautonomous Oscillatory System with Nonlinear Capacity. Phys. Rev. A 1987, 35, 4404. [Google Scholar] [CrossRef]
  8. Alsedà, L.; Costa, S. On the Definition of Strange Nonchaotic Attractor. Fundam. Math. 2009, 206, 23–39. [Google Scholar] [CrossRef]
  9. Ditto, W.L.; Spano, M.L.; Savage, H.T.; Rauseo, S.M.; Heagy, J.; Ott, E. Experimental observation of a strange nonchaotic attractor. Phys. Rev. Lett. 1990, 65, 533. [Google Scholar] [CrossRef]
  10. Arulgnanam, A.; Prasad, A.; Thamilmaran, K.; Daniel, M. Analytical Study and Experimental Confirmation of SNA through Poincaré Maps in a Quasiperiodically Forced Electronic Circuit. Int. J. Bifurc. Chaos 2015, 25, 1530020. [Google Scholar] [CrossRef]
  11. Bezruchko, B.P.; Kuznetsov, S.P.; Seleznev, Y.P. Experimental observation of dynamics near the torus-doubling terminal critical point. Phys. Rev. E 2000, 62, 7828. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Uenohara, S.; Mitsui, T.; Hirata, Y.; Morie, T.; Horio, Y.; Aihara, K. Experimental distinction between chaotic and strange nonchaotic attractors on the basis of consistency. Chaos 2013, 23, 023110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wang, X.; Zhan, M.; Lai, C.-H.; Lai, Y.-C. Strange Nonchaotic Attractors in Random Dynamical Systems. Phys. Rev. Lett. 2004, 92, 074102. [Google Scholar] [CrossRef] [Green Version]
  14. Deng, Q.; Wang, C. Multi-scroll hidden attractors with two stable equilibrium points. Chaos 2019, 29, 093112. [Google Scholar] [CrossRef]
  15. Leonov, G.; Kuznetsov, N. Hidden attractors in dynamical systems. From hidden oscil-lations in Hilbert-Kolmogorov, Aizerman, and Kalman problems to hidden chaotic attractorsin Chua circuits. Int. J. Bifurc. Chaos 2013, 23, 1330002. [Google Scholar] [CrossRef] [Green Version]
  16. Leonov, G.; Kuznetsov, N.; Mokaev, T. Homoclinic orbits, and self-excited and hidden attractors in a Lorenz-like system describing convective fluid motion. Eur. Phys. J. Special Top. 2015, 224, 1421–1458. [Google Scholar] [CrossRef] [Green Version]
  17. Leonov, G.; Kuznetsov, N.; Vagaitsev, V. Localization of hidden Chua’s attractors. Phys. Lett. A 2011, 375, 2230–2233. [Google Scholar] [CrossRef]
  18. Zhang, X.; Wang, C. Multiscroll Hyperchaotic System with Hidden Attractors and Its Circuit Implementation. Int. J Bifurc. Chaos 2019, 29, 1950117. [Google Scholar] [CrossRef]
  19. Zhou, L.; Wang, C.; Zhou, L. A novel no-equilibrium hyperchaotic multi-wing system via introducing memristor. Int. J. Circ. Theo. App. 2018, 46, 84. [Google Scholar] [CrossRef]
  20. Jafari, S.; Sprott, J.C. Simple chaotic flows with a line equilibrium. Chaos Soliton Fract. 2013, 57, 79–84. [Google Scholar] [CrossRef]
  21. Jafari, S.; Sprott, J.C.; Golpayegani, S.M.R.H. Elementary quadratic chaotic flows with no equilibria. Phys. Lett. A 2013, 377, 699. [Google Scholar] [CrossRef]
  22. Wang, X.; Chen, G. A chaotic system with only stable equilibrium. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1264–1272. [Google Scholar] [CrossRef] [Green Version]
  23. Wei, Z. Dynamical behaviors of a chaotic system with no equilibria. Phys. Lett. A 2011, 376, 102–108. [Google Scholar] [CrossRef]
  24. Danca, M.-F.; Feckan, M.; Kuznetsov, N.; Chen, G. Looking more closely to the Rabinovich-Fabrikant system. Int. J. Bifurc. Chaos 2016, 26, 1650038. [Google Scholar] [CrossRef] [Green Version]
  25. Danca, M.-F.; Kuznetsov, N.; Chen, G. Unusual dynamics and hidden attractors of the Rabinovich-Fabrikant system. Nonlinear Dynam. 2017, 88, 791–805. [Google Scholar] [CrossRef]
  26. Danca, M.-F. Hidden transient chaotic attractors of Rabinovich-Fabrikant system. Nonlinear Dynam. 2016, 86, 1263–1270. [Google Scholar] [CrossRef] [Green Version]
  27. Babuska, I.; Prager, M.; Vitasek, E. Numerical Processes in Differential Equations, SNTL, Praha; John Wiley & Sons (Interscience Publishers): New York, NY, USA, 1966. [Google Scholar]
  28. Teixeira, J.; Reynolds, C.A.; Judd, K. Time-step sensitivity of non-linear atmospheric model. Numerical convergence, truncation error growth and ensemble design. J. Atmos. Sci. 2007, 64, 175–189. [Google Scholar] [CrossRef] [Green Version]
  29. Lambert, J.D. Numerical Methods for Ordinary Differential Systems: The Initial Value Problem; Wiley: New York, NY, USA, 1992. [Google Scholar]
  30. René, L.; Pogonin, V.; Pchelintsev, A. A new accurate numerical method of approximation of chaotic solutions of dynamical model equations with quadratic nonlinearities. Chaos Soliton Fract. 2016, 91, 108–114. [Google Scholar]
  31. Sarra, S.A.; Meador, C. On the numerical solution of chaotic dynamical systemsusing extend precision floating point arithmeticand very high order numerical methods. Nonlinear Anal. Modell. 2011, 16, 340–352. [Google Scholar] [CrossRef] [Green Version]
  32. Wang, P.; Li, J.; Li, Q. Computational uncertainty and the application of a high-performance multiple precision scheme to obtaining the correct reference solution of Lorenz equations. Numer. Algorithms 2012, 59, 147–159. [Google Scholar] [CrossRef]
  33. Liao, S.J.; Wang, P.F. On the mathematically reliable long-term simulation of chaotic solutions of Lorenz equation in the interval [0, 10,000]. Sci. China Phys. Mech. 2014, 57, 330–335. [Google Scholar] [CrossRef] [Green Version]
  34. Grebogi, C.; Ott, E.; Yorke, J.A. Fractal Basin Boundaries, Long-Lived Chaotic Transients, and Unstable-Unstable Pair Bifurcation. Phys. Rev. Lett. 1983, 51, 942. [Google Scholar] [CrossRef]
  35. Grebogi, C.; Ott, E.; Yorke, J.A. Superpersistent chaotic transients. Ergod. Th. Dynam. Syst. 1985, 5, 341–372. [Google Scholar] [CrossRef] [Green Version]
  36. Lai, Y.C.; Tél, T. Transient Chaos: Complex Dynamics on Finite Time Scales; Part of the Applied Mathematical Sciences Book Series (AMS, Volume 173); Springer: New York, NY, USA, 2011. [Google Scholar]
  37. Danca, M.-F.; Bourke, P.; Kuznetsov, N. Graphical Structure of Attraction Basins of Hidden Chaotic Attractors: The Rabinovich-Fabrikant System. Int. J. Bifurc. Chaos 2019, 29, 1930001. [Google Scholar] [CrossRef] [Green Version]
  38. Ding, M.; Grebogi, C.; Ott, E. Evolution of attractors in quasi-periodically forced systems. From quasiperiodic to strange nonchaotic to chaotic. Phys. Rev. A 1989, 39, 2593. [Google Scholar] [CrossRef] [PubMed]
  39. Gottwald, G.A.; Melbourne, I. On the implementation of the 0-1 test for chaos. SIAM J. Appl. Dyn. Syst. 2009, 8, 129–145. [Google Scholar] [CrossRef]
  40. Ramaswamy, R. Synchronization of strange nonchaotic attractors. Phys. Rev. E 1997, 56, 7294. [Google Scholar] [CrossRef] [Green Version]
  41. Negi, S.S.; Ramaswamy, R. A plethora of Strange non-chaotic Attractors. Pramana 2001, 56, 47–56. [Google Scholar] [CrossRef] [Green Version]
  42. Gopal, R.; Venkatesan, A.; Lakshmanan, M. Applicability of 0-1 test for strange nonchaotic attractors. Chaos 2013, 23, 023123. [Google Scholar] [CrossRef] [Green Version]
  43. Klages, R. Weak Chaos, Infinite Ergodic Theory, and Anomalous Dynamics, from Hamiltonian Chaos to Complex Systems; Part of the Nonlinear Systems and Complexity Book Series (NSCH, Volume 5); Springer: New York, NY, USA, 2013; pp. 3–42. [Google Scholar]
  44. Zaslavsky, G.M.; Sagdeev, R.Z.; Usikov, D.A.; Chernikov, A.A. Weak Chaos and Quasi-Regular Patterns; Cambridge Nonlinear Science Series; Cambridge University Press: Cambridge, MA, USA, 2001. [Google Scholar]
  45. Hu, J.; Tung, W.W.; Gao, J.; Cao, Y. Reliability of the 0-1 test for chaos. Phys. Rev. E 2005, 72, 056207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Marwan, N.; Romano, M.; Thiel, M.; Kurths, J. Recurrence Plots for the Analysis of Complex Systems. Phys. Rep. 2007, 438, 237–329. [Google Scholar] [CrossRef]
  47. Nicol, M.; Melbourne, I.; Ashwin, P. Euclidean extensions of dynamical systems. Nonlinearity 2001, 14, 275. [Google Scholar] [CrossRef] [Green Version]
  48. Gottwald, G.A.; Melbourne, I. A new test for chaos in deterministic systems. Proc. R. Soc. London Ser. A 2004, 460, 603–611. [Google Scholar] [CrossRef] [Green Version]
  49. Poincaré, H. Chapitre III. Théorie des Solutions Périodiques. Available online: https://projecteuclid.org/journals/acta-mathematica/volume-13/issue-1-2/Chapitre-III-Th%C3%A9orie-des-solutions-p%C3%A9riodiques/10.1007/BF02392509.full (accessed on 19 November 2020). [CrossRef]
  50. Eckmann, J.P.; Kamphorst, S.O.; Ruelle, D. Recurrence plots of dynamical systems. Europhys Lett. 1987, 4, 973. [Google Scholar] [CrossRef] [Green Version]
  51. Marwan, N.; Donges, J.F.; Zou, Y.; Donner, R.V.; Kurths, J. Complex network approach for recurrence analysis oftime series. Phys. Lett. A 2009, 373, 4246–4254. [Google Scholar] [CrossRef] [Green Version]
  52. Mindlin, G.M.; Gilmore, R. Topological analysis and synthesis of chaotic time series. Physica D 1992, 58, 229–242. [Google Scholar] [CrossRef]
  53. Thiel, M.; Romano, M.C.; Kurths, J.; Meucci, R.; Allaria, E.; Arecchi, T. Influence of observational noise on the recurrence quantification analysis. Physica D 2002, 171, 138–152. [Google Scholar] [CrossRef] [Green Version]
  54. Ngamga, E.J.; Nandi, A.; Ramaswamy, R.; Romano, M.C.; Thiel, M.; Kurths, J. Recurrence analysis of strange nonchaotic dynamics. Phys. Rev. E 2007, 75, 036222. [Google Scholar] [CrossRef] [Green Version]
  55. Grassberger, P.; Procaccia, I. Measuring the strangeness of strange attractors. Physica D 1983, 9, 189–208. [Google Scholar] [CrossRef]
Figure 1. (a) Variation in lifetime intervals [ 0 , T ] of the numerical trajectories; (b) The variation in Λ of the RF system for b [ 0.2808 , 0.2907 ] .
Figure 1. (a) Variation in lifetime intervals [ 0 , T ] of the numerical trajectories; (b) The variation in Λ of the RF system for b [ 0.2808 , 0.2907 ] .
Mathematics 09 00652 g001
Figure 2. Attractors of the RF system. (a) SECA for b = 0.98 . Trajectories from unstable equilibria X 3 , 4 * are attracted by the SECA, while trajectories starting from X 0 * tend to infinity; (b) quasiperiodic attractor for b = 0.25675 after the first transients are removed; (c) HSNA for b = 0.2875 . Beside the attractor H (green plot), trajectories are presented from unstable equilibria X 0 * and X 3 , 4 * , which are attracted either by the stable equilibria X 1 , 2 * or diverge to infinity (see the legend). The dotted plot in the zoomed image D reflects the fractal autosimilariity characteristic of H .
Figure 2. Attractors of the RF system. (a) SECA for b = 0.98 . Trajectories from unstable equilibria X 3 , 4 * are attracted by the SECA, while trajectories starting from X 0 * tend to infinity; (b) quasiperiodic attractor for b = 0.25675 after the first transients are removed; (c) HSNA for b = 0.2875 . Beside the attractor H (green plot), trajectories are presented from unstable equilibria X 0 * and X 3 , 4 * , which are attracted either by the stable equilibria X 1 , 2 * or diverge to infinity (see the legend). The dotted plot in the zoomed image D reflects the fractal autosimilariity characteristic of H .
Mathematics 09 00652 g002
Figure 3. The SECA for b = 0.98 . (a) The variation in the Lyapunov spectrum; (b) the variation in the median K; (c) power spectral density; (d) the graph of p and q; (e) the graph of the mean-square displacement M; (f) the recurrence plot.
Figure 3. The SECA for b = 0.98 . (a) The variation in the Lyapunov spectrum; (b) the variation in the median K; (c) power spectral density; (d) the graph of p and q; (e) the graph of the mean-square displacement M; (f) the recurrence plot.
Mathematics 09 00652 g003
Figure 4. The quasiperiodic attractor for b = 0.25675 . (a) The variation in the Lyapunov spectrum; (b) the variation in the median K; (c) power spectral density; (d) the graph of p and q; (e) the graph of the mean-square displacement M; (f) the recurrence plot and zoomed details.
Figure 4. The quasiperiodic attractor for b = 0.25675 . (a) The variation in the Lyapunov spectrum; (b) the variation in the median K; (c) power spectral density; (d) the graph of p and q; (e) the graph of the mean-square displacement M; (f) the recurrence plot and zoomed details.
Mathematics 09 00652 g004
Figure 5. (a) Three-dimensional parallelipipedic representation of the attraction basin (red plot) of the chaotic attractor H , determined around the unstable equilibrium X 0 * ; (b) zoomed area, revealing the fact that there are no initial points within a small cubic neighborhood of X 0 * ; (c) the same zoomed area presenting, as a supplement, the initial conditions leading to the stable equilibria X 1 , 2 * (blue and red plot).
Figure 5. (a) Three-dimensional parallelipipedic representation of the attraction basin (red plot) of the chaotic attractor H , determined around the unstable equilibrium X 0 * ; (b) zoomed area, revealing the fact that there are no initial points within a small cubic neighborhood of X 0 * ; (c) the same zoomed area presenting, as a supplement, the initial conditions leading to the stable equilibria X 1 , 2 * (blue and red plot).
Mathematics 09 00652 g005
Figure 6. The HSNA H for b = 0.2875 . (a) The variation in the Lyapunov spectrum; (b) the variation in the median K; (c) power spectral density; (d) the graph of p and q; (e) the graph of the mean-square displacement M; (f) the recurrence plot.
Figure 6. The HSNA H for b = 0.2875 . (a) The variation in the Lyapunov spectrum; (b) the variation in the median K; (c) power spectral density; (d) the graph of p and q; (e) the graph of the mean-square displacement M; (f) the recurrence plot.
Mathematics 09 00652 g006
Figure 7. Difference between the first components x 1 and y 1 of two close trajectories x and y. (a) The SECA for b = 0.98 ; (b) the quasiperiodic attractor corresponding to b = 0.25675 ; (c) the HSNA H ; (i) overplot of x 1 and y 1 with zoomed details and three-dimensional representations; (ii) details of (i).
Figure 7. Difference between the first components x 1 and y 1 of two close trajectories x and y. (a) The SECA for b = 0.98 ; (b) the quasiperiodic attractor corresponding to b = 0.25675 ; (c) the HSNA H ; (i) overplot of x 1 and y 1 with zoomed details and three-dimensional representations; (ii) details of (i).
Mathematics 09 00652 g007
Figure 8. Histograms of the component x 2 . (a) The SECA corresponding to b = 0.98 ; (b) the quasiperiodic attractor corresponding to b = 0.25675 ; (c) the HSNA H corresponding to b = 0.2875 .
Figure 8. Histograms of the component x 2 . (a) The SECA corresponding to b = 0.98 ; (b) the quasiperiodic attractor corresponding to b = 0.25675 ; (c) the HSNA H corresponding to b = 0.2875 .
Mathematics 09 00652 g008
Figure 9. (a) LEs of the RF system for b [ 0.255 , 0.289 ] ; (b) median K for b [ 0.255 , 0.289 ] .
Figure 9. (a) LEs of the RF system for b [ 0.255 , 0.289 ] ; (b) median K for b [ 0.255 , 0.289 ] .
Mathematics 09 00652 g009
Table 1. Characterization of the attractors of the RF system.
Table 1. Characterization of the attractors of the RF system.
Motion Λ ‘0-1’ TestPSDRP
(A) chaotic (Figure 2a) λ 1 > 0 > λ 2 , 3 K 1 ;
q , p Brownian;
M linear growth
broadband irregular
structures
(B) quasiperiodic (Figure 2b) 0 = λ 1 λ 2 , 3 K 0 ;
q , p disk;
M oscillatory
many discrete
peaks
non-equidistant
diagonal segments
and structures
(C) SNA (Figure 2c) 0 > λ 1 , 2 , 3 0 < K < 1 ;
q , p disk-like;
M bounded oscillatory growth
broadband with
discrete peaks
characteristics
between (A) and (B)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Danca, M.-F.; Kuznetsov, N. Hidden Strange Nonchaotic Attractors. Mathematics 2021, 9, 652. https://doi.org/10.3390/math9060652

AMA Style

Danca M-F, Kuznetsov N. Hidden Strange Nonchaotic Attractors. Mathematics. 2021; 9(6):652. https://doi.org/10.3390/math9060652

Chicago/Turabian Style

Danca, Marius-F., and Nikolay Kuznetsov. 2021. "Hidden Strange Nonchaotic Attractors" Mathematics 9, no. 6: 652. https://doi.org/10.3390/math9060652

APA Style

Danca, M. -F., & Kuznetsov, N. (2021). Hidden Strange Nonchaotic Attractors. Mathematics, 9(6), 652. https://doi.org/10.3390/math9060652

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