CC BY
29
УДК 517.929:621.373
ROBUST CHAOS IN AUTONOMOUS TIME-DELAY SYSTEM
D. S. Arzhanukhina1 and S. P. Kuznetsov 2'3
1 Saratov State University 2Kotel'nikov's Institute of Radio-Engineering and Electronics of RAS, Saratov Branch 3Department of Physics and Astronomy, University of Potsdam
We consider an autonomous system constructed as modification of the logistic differential equation with delay that generates successive trains of oscillations with phases evolving according to chaotic maps. The system contains two feedback loops characterized by two generally distinct retarding time parameters. In the case of their equality, chaotic dynamics is associated with the Smale-Williams attractor that corresponds to the double-expanding circle map for the phases of the carrier of the oscillatory trains. Alternatively, at appropriately chosen two different delays attractor is close to torus with Anosov dynamics on it as the phases are governed by the Fibonacci map. In both cases the attractors manifest robustness (absence of regularity windows under variation of parameters) and presumably relate to the class of structurally stable hyperbolic attractors.
Keywords: Attractor, hyperbolic chaos, maps, Anosov dynamics, Arnold cat, Fibonacci map, Smale-Williams attractor.
Introduction
The concept of uniformly hyperbolic chaotic dynamics was advanced in mathematical theory of dynamical systems half a century ago [1-6]. It deals with invariant sets in state space of systems composed exclusively of orbits of saddle type with such a restriction that their stable and unstable manifolds do not touch each other (intersect transversally). The hyperbolic invariant sets are structurally stable. It implies their robustness in respect to (at least small) variation of functions and parameters in the dynamical equations. Henceforth, the hyperbolic dynamics may be regarded as preferable for any possible practical application of chaos and as deserving prior research in a frame of physical and technical disciplines [7, 8].
One important class of the uniform hyperbolicity is Anosov dynamics in systems where the hyperbolic invariant set occupies the whole accessible phase space; usually it is considered in the context of the phase volume preserving maps or flows (conservative dynamics). Alternatively, in the context of dissipative systems, the uniformly hyperbolic attractors were introduced; their mathematical examples are Smale-Williams solenoid,
Plykin type attractors, DA-attractor of Smale1. Although originally they relate to artificial discrete-time systems, analogous attractors can occur as well in Poincare maps associated with continuous time systems governed by differential equations. In the last case it is used to speak about suspension of these attractors provided by the appropriate flows.
Till a short time ago, realistic examples of uniformly hyperbolic dynamics which could relate to real-world systems were poorly presented and discussed in the literature. One example is a hinge mechanism with Anosov dynamics described and studied in [9]. An artificial example of suspension of Plykin-type attractor was constructed in PhD thesis of Hunt [10], but it is surely too complicated to allow a real physical implementation. Possible occurrences of suspension of Plykin-type attractors in neuron model [11] and in modified Lorenz model [12] were discussed, but no convincing data on a level of concrete equations and numerical simulation were provided.
Recently a number of implementable systems with uniformly hyperbolic attractors were advanced and studied due to efforts of Saratov group of nonlinear dynamics [13-20]. It is essential that these examples are constructed on a base of the physical rather than the mathematical toolbox exploiting such entities as oscillators, nonlinear elements, interactions, self-oscillatory and parametric excitation etc.
One particular productive method for design of systems with uniformly hyperbolic attractors is based on using time-delay nonlinear feedback loops [21-24]. Most of the examples suggested were non-autonomous systems functioning in presence of the external periodic driving [21-23]. An exception is an autonomous system with attractor of Smale-Williams type considered in Ref. [24]. The present article is inspired by that work. We construct and study numerically an autonomous time-delay system manifesting different types of uniformly hyperbolic attractors depending on two parameters of time delay.
1. Basic equations
Let us start with logistic delay equation offered in due time in population biology [25]:
r = |[1 - r(t - x1)]r(t). (1)
Here the population is characterized by the positive variable r evolving in time and normalized in such way that the saturation occurs at r = 1. A positive parameter | is the birth rate, ti is the delay time characterizing lag of the effect of saturation. According to Ref. [25], under condition t1 < n/2| the system has a stable stationary state r = 1 and for t > n/2| self-oscillations occur. At large values t1 the generated waveform looks like a periodic sequence of pulses (Fig. 1, a). The period grows with the delay ti as P = (1 + e|Tl )/|. Estimate for the minimal level of the population in between the pulses is rmin = |t exp(-e|Tl + 2|Ti - 1), i.e. it manifests the double exponentially decrease with increase of t1.
Now, let us regard the variable r as squared amplitude of some oscillatory process with frequency wq. For this, we set r = x2 + y2 and require the new variables to satisfy
1The abbreviation DA stands for «derived from Anosov».
the equations
x = -Шоу + (1/2)ц(1 - x2(t - Xi) - y2(t - Ti))x, y = Шох + (1/2)ц(1 - x2(t - Ti) - y2(t - Ti))y.
Solutions of the equations (2) instead of the solitary pulses manifest successive oscillatory trains for x and y while the envelope evolves exactly according to (1) (Fig. 1, b).
Next, let us add terms proportional to x(t — T1)x(t — t2) — y(t — T1)y(t — t2) and x(t — T1)y(t — т2) + x(t — T2)y(t — т1) with small coefficient e in the first and the second equation
x = —Шоу + (1/2)ц(1 — x2(t — Ti) — y2(t — Ti))x+
(3)
+e[x(t - ti)x(t - x2) - y(t - Ti)y(t - t2)], il = mox + (1/2)^(1 - x2(t - Ti) - y2(t - Ti))y+
+e[x(t - Ti)y(t - T2) + x(t - T2)y(t - Ti)],
where we suppose always that the delays satisfy the inequality t2 > Ti. Now, in the case of generation of pulses with low enough level of minimal amplitude between them, just these additional terms will initiate formation of each next pulse of oscillations. Due to this, the phase of the oscillations for the new-born pulse will be determined by the phases of the previous pulses via some mapping as discussed in detail in the following sections. In presence of this initiation of the oscillations the characteristic repetition period of the pulses is reduced comparing to the original system (2).
Formally, the time-delay system is infinite-dimensional. Indeed, an instantaneous state is determined by functions x(t) and y(t) defined on a time interval of length t2. To start solution of the equations (3) the initial conditions have to be chosen as some functions x(t) and y(t) on the time interval [-t2, 0]. In the course of numerical integration they are represented by arrays sampled with time step used in the integration scheme, and the elements of the arrays are replaced step by step by the newly obtained values to be accounted on latter stages of the computations. In other respects the integration schemes
Рис. 1. Self-pulsations in the logistic-delay equation (1) at ^ = 1.6, t1 =2 (a) and oscillatory trains generated by system (2) at ю0 = 2n (b). The period of the pulses is roughly T ~ 16
are analogous to those used for solution of ordinary differential equations. With arbitrary initial conditions the computations are executed for sufficiently long time and then, after reliable arrival at the attractor, processing of the data is performed including plotting waveforms, portraits of attractors, evaluation of the Lyapunov exponents etc.
Now let us turn to some concrete cases of the model (3).
2. Attractor of Smale-Williams type
Let us start with the simplest particular case setting both delay time parameters equal: т2 = ti = т. It leads to the equations
x = -rnoV + (1/2)ц(1 - x2(t - т) - y2(t - т))х + e[x2(t - т) - y2(t - т)],
(4)
y = mox + (1/2)ц(1 - x2(t - т) - y2(t - т))у + 2ex(t - т)y(t - т).
To clarify the principle of functioning of the system suppose that during some current stage of activity the variables behave as x ~ cos(rn0t + ф) and y ~ sin(rn0t + ф); then the additional terms are expressed as
x2(t - т) - y2(t - т) ~ cos2(rn0t + ф) - sin2(rn0t + ф) = cos(2rn0t + 2ф),
(5)
2x(t - т)у (t - т) ~ 2 cos(rn0t + ф) sin(rn0t + ф) = sin(2rn0t + 2ф).
These relations determine actually the stimulating signal for oscillations on the immediate next stage of activity. The phase shift of this signal is transferred to the oscillations arising on this stage. Hence, the phase undergoes the doubling transformation
фга+1 = 2фга + const. (6)
It is the expanding circle map, or the Bernoulli map, which is chaotic and is characterized by the Lyapunov exponent Л = ln 2 & 0.693.
We may consider a Poincare map for the time-delay system that corresponds to transformation of the infinite-dimensional space vector on a time interval from one pulse of oscillations to the next one. Attractor of this Poincare map is supposed to be a kind of Smale-Williams solenoid embedded in the infinite-dimensional state space. Actually, the phase ф plays a role of the angular variable in the Smale-Williams solenoid while in other directions compression of the phase volume will take place. Respectively, attractor of the original autonomous time-delay system (4) is a suspension of the Smale-Williams solenoid in the infinite-dimensional state space.
Figure 2 shows waveforms for dynamical variables x and y illustrating operation of the system (4) at
ю0 = 2п, ^ = 1.6, т = 2, e = 0.05. (7)
In accordance with the above qualitative considerations, the process looks like a sequence of trains of oscillations. Note that the average period of their appearance is less than that in models (1), (2), and in the regime under discussion it is T & 10. Although represented by nearly periodic alternation of excitation and suppression of the oscillations, the process is actually chaotic (in contrast to the regular one observed in Fig. 1, b). Chaos reveals
Рис. 2. Waveforms for the dynamical variables x (black) and y (gray) according to the results of the numerical solution of equation (4) at ^ = 1.6, ю0 = 2я, т = 2, e = 0.05
itself in variations of phases on successive stages of activity. More accurate analysis shows that these phases obey approximately the Bernoulli map. Figure 3 illustrates it by a diagram obtained in the numerical simulation. The phases were determined at moments of maximal amplitude at successive excitation stages as фп = = arg(x - iy) and plotted for sufficiently large number of processed successive oscillatory trains. Surely, the mapping for the phases looks topologically equivalent to the expected Bernoulli map: one complete bypass for the pre-image фп (that is variation by 2n) corresponds to the twofold bypass for the image фга+ь
Fig.4 shows 3D and 2D projections of the attractor from the infinite-dimensional phase space of our time-delay system.
Рис. 3. Diagram illustrating the transformation of phases in the successive stages of the activity plotted according to the results of the numerical solution of equation (4) at ^ = 1.6, ю0 = 2я, т = 2, e = 0.05
Рис. 4. 3D and 2D projections of the attractor from the infinite-dimensional phase space at ю0 = 2п, = 1.6, т = 2, e = 0.05. The first diagram is plotted in the coordinates (x, y, p), where p = \Jx2(t — т) + y2(t — т), and the second on the plane of variables (x, y)
For computation of Lyapunov exponents we use the Benettin algorithm [26, 27] adapted for the time-delayed systems [28, 29]. It is based on simultaneous numerical solutions of the equations (4) and the variation equations
X = -шоу + (1/2)ц[1 - x2(t - т) - y2(t - т)]Х-
-^x[x(t - T)X(t - т) + y(t - T)y(t - т)] + 2e[x(t - T)X(t - т) - y(t - T)y(t - т)],
у = шоХ + (1/2)ц[1 - x2(t - т) - y2(t - т)]у-
-^y[x(t - т)Х^ - т) + y(t - т)у^ - т)] + e[y(t - т)Х^ - т) + x(t - т)у^ - т)].
(8)
Formally, there are an infinite number of Lyapunov exponents in the time-delayed system, but we necessarily restrict ourselves dealing with a finite number of them. To evaluate the larger M exponents we integrate a collection of M sets of equations of form (8). In the computations, the perturbation vectors are represented as finite-dimensional arrays of values X and у sampled on intervals of length т with time step of the integration scheme of the delay equations. The dot product for the vectors involved in normalization and Gram-Schmidt orthogonalization required in the Benettin algorithm is determined as a sum of products of pairs of the elements relating to one and another array.
It is convenient to use normalization of the Lyapunov exponents Xj by the average period of pulse repetition T, namely, to set Л = XT. It makes natural to relate them with the Lyapunov exponent of the Bernoulli map obtained in the qualitative analysis.
According to the computations, in the regime corresponding to parameters (7) the largest four Lyapunov exponents are
Л1 = 0.6964, Л2 = 0.0000,
Л3 = -12.129, Л4 = -15.538, ...
(9)
The estimate of the Kaplane-Yorke dimension for the attractor of the Poincare map (the Smale-Williams solenoid) is D% = 1 + Л1/|Л3| ^ 1.06. Dimension of the attractor in the infinite-dimensional phase space of the autonomous system is Dky = 2+(л1+л2)/|лз| = 1+D°y И 2.06.
Figure 5 shows plots for four Lyapunov exponents of the system (4) versus e, т and ц. Observe that the largest exponent remains almost constant in a wide interval of the parameter variation
Рис. 5. Lyapunov exponents of the system (4) versus the parameters e (a), т (b), ц (c)
and remains close to the expected value ln2 = 0.693.... The second is close to zero (up to numerical inaccuracy). As usual in autonomous systems, it is interpreted as associated with perturbations of infinitesimal shift along the phase trajectory. Other exponents are negative.
Note smooth character of dependences of the Lyapunov exponents on parameters, particularly, there are no drops in the plot of A1, which could correspond to regularity windows intrinsic usually to non-hyperbolic systems. It confirms robust nature of chaos, which is natural in a frame of our qualitative argumentation that it is associated with the uniformly hyperbolic attractor of Smale-Williams type.
3. Attractor associated with Anosov-type dynamics on torus
Now return to the model with two different time delays t1 and t2:
X = -woy + (1/2)|(1 - x2(t - ti) - y2(t - Ti))x+
+e[x(t - Ti)x(t - T2) - y(t - Ti)y(t - T2)],
(10)
y = Wox + (1/2)|(1 - x2(t - Ti) - y2(t - Ti))y+
+e[x(t - Ti)y(t - T2) + x(t - T2)y(t - Ti)]
and select these parameters to get such situation that the delayed signals stimulating excitation of a new (n + 1)-th pulse of oscillations arrive being emitted from the previous two pulses n (delay t1) and n - 1 (delay t2).
Suppose that at the previous two activity stages the phases were q>n and ^n-1, namely, x(t-t1) ~ cos(rn0t+^n), y(t-t1) ~ sin(rn0t+^n), x(t-t2) ~ cos(m0t+^n-1), yx(t - t2) ~ sin(w0t + ^n-1). Then, the additional terms in the right-hand parts are expressed as
x(t - Ti)x(t - T2) - y(t - Ti)y(t - T2) ~ cos(2wot + q>n + tyn_i),
(11)
x(t - Ti)y(t - T2) + x(t - T2)y(t - Ti) ~ sin(2mot + ^n + 9n_i).
As follows, on the next, (n + 1)-th stage of activity the phase will be determined (up to a constant) by the Fibonacci map
^n+i = ^n + ^n-1 + const (mod2n). (12)
Of what kind attractor is in this case? For the infinite-dimensional Poincare map of our time-delay system attractor is an object close geometrically to a two-dimensional torus, and the discrete-time dynamics of the angular variables on this torus obey the map (12). It relates to the class of Anosov maps. Due to structural stability of the Anosov dynamics, one can conjecture that on the attractor in the infinite-dimensional state space of the Poincare map the dynamics remains of the same nature; in other dimensions the phase volume compression occurs that corresponds to the approach of orbits to the attractor. Respectively, attractor in the phase space of the continuous-time system (10) is a suspension of that object, and here (in the autonomous system) we have an additional neutral direction associated with a zero Lyapunov exponent.
Figure 6 shows waveforms for the variables x and y obtained from numerical solution of the equations (10) at
^ = 1.6, e = 0.02, = 2n, xi =2, t2 = 14. (13)
The process looks like a sequence of trains of oscillations for which the phases of the high-frequency carrier vary in a random-like way from one to the next stage of activity. The average period of repetition of the pulses evaluated numerically is roughly T « 10.85.
Figure 7 shows 3D and 2D projections of the attractor from the infinite-dimensional phase space of our time-delay system on the plane of variables (x, y).
Figure 8 illustrates correspondence of the dynamics of the phases to the Fibonacci map. The phases are determined at instants of maximal amplitude in the oscillation trains as ty = arg(x — iy). Diagram (a) presents the data in the coordinates (tyn + tyn-1, and one can see that the dots lie in a strip parallel to the bisector. (Although the strip looks widened, this does not violate the expected topological nature of the mapping.) Diagram (b) shows the 3D diagram for the phases in coordinates (tyn+1, tyn-1), which corresponds visually to the Fibonacci map.
0 22.5 45 67.5 90.0 t
Рис. 6. Waveforms for variables x (black) and y (gray) according to the numerical solution of equation (10) at ц = 1.6, e = 0.02, шо = 2я, тх =2, т2 = 7тi = 14
З.О3.0 -3.0 X 3.0
Рис. 7. 3D and 2D projections of the attractor from the infinite-dimensional phase space at ц = 1.6, e = 0.02, ш0 = 2я, тх = 2, т2 = 7тх = 14. The first diagram is plotted in the coordinates (x,y, p), where p = \Jx2(t — т) + y2(t — т), and the second on the plane of variables (x, y)
ф„+ ф„-1
Рис. 8. Diagrams illustrating dynamics of the phases of the system (10) in successive stages of activity
One more way to illustrate the Anosov dynamics for the phases is similar to the famous picture of Arnold's cat [30, 15, 16]. For each pair of phases фп and фга_1 on the n-th and n -1 stages of activity, we determine either the point with such coordinates fit or not in the picture of the cat face drawn on the plane. If yes, then the point is marked on the graph, and on the following two plots the corresponding points are depicted, respectively after 3 and 6 iteration steps. If the initial point does not fall within the specified cat face area, the dots are not marked, and we carry on with further iterations. Figure 9 shows the pictures obtained by this method comparing with the respective diagram originating from iterations of the Fibonacci map itself.
0 4Vi 2k
Рис. 9. Transformation of the cat face area under successive iterations of the mapping for the phases obtained in computations for the model (10) (a) and for the Fibonacci map (b). (See explanations in the text.)
For computation of the Lyapunov exponents we use the algorithm similar to that described in the previous section. It is based on the joint numerical solutions of the equations (10) and the variation equations
x = -rnoy + (1/2)p,[1- x2(t-Ti)-y2(t-Ti)]x - px(x(t-ti)x(t - ti) + y(t-ti)y(t-ti))+ +e[x(t - Ti)x(t - T2) + x(t - T2)x(t - Ti) - y(t - Ti)y(t - T2) - y(t - r2)y(t - ti)], y= wox+(1/2)^(1 - x2(t - ti)- y2(t-ti))y-^y(x(t-ti)x(t-ti)+ y(t-Ti)y(t-ti))+ +e[y(t - t2)x(t - ti) + y(t - t2)x(t - T2) + x(t - t2)y(t - ti) + y(t - ti)x(t - T2)].
(14)
The perturbation vectors are represented in computations as finite arrays of x and y on time intervals of length t2 sampled with the integration step used for the numerical solution of the equations. We compute four larger exponents and normalize them by the factor T; at parameters (13) they are
Ai = 0.4851, A2 = 0.0003,
A3 = -0.4691, A4 = -0.5404...
(15)
The first and the third Lyapunov exponents are close in magnitude and opposite in sign. Their may be compared with the Lyapunov exponents of the Fibonacci map A = ± ln(1 + V5)/2 = ±0.4812.... The second exponent is zero up to numerical inaccuracy. The fourth and subsequent Lyapunov exponents are all negative, of larger absolute values.
Figure 10 shows plots for the Lyapunov exponents depending on the coupling parameter e (a), parameter of activity ^ (b) and the delay ti (c) at fixed other parameters. As seen, in a fairly wide range of the parameters the largest Lyapunov exponent remains approximately constant, and one of the negative exponents is close in the absolute value to the positive exponent. It indicates persistence of the Anosov dynamics on the attractor.
Viewing the graphics in Fig. 10 one can distinguish two somewhat different
Рис. 10. Lyapunov exponents of the system (10) depending on the parameters e, ^ and т
situations. Namely, in the left parts of the diagrams the first negative exponent is nearly equal to the positive one in absolute value, and the next negative exponent is lesser. As we go to the right some crossover occurs, and the two negative exponents exchange their places. After that, obviously, approach of trajectories to the attractor in time is slower than the convergence and the divergence of orbits on the attractor itself. (Apparently, it implies different nature of the attractors in these two situations, but we leave the comparative analysis of them outside the frame of the present article.)
We outline again the smooth character of dependences of the Lyapunov exponents on parameters that confirms robust nature of chaos. It may be thought that the attractors we deal with in this section are uniformly hyperbolic, close to the so-called «non-strange chaotic attractors» briefly mentioned and discussed in the literature [31-33].
Conclusion
In this paper, we introduce an autonomous system built on the basis of the logistic differential equation with delay. As shown, the dynamics of the phases of generated successive oscillatory trains corresponds to expanding circle map or Anosov map on a torus depending on selection of two retarding time parameters. Numerical simulations show robust chaos generation in these cases, and we conjecture the uniformly hyperbolic nature on the attractors in the infinite-dimensional state space of the time-delay system.
The suggested system and analogous constructions may be realized on a basis of electron devices as generators of robust chaos; physically they may be simpler in implementation comparing, for example, with systems constructed on the basis of coupled oscillators [13-16].
Due to the insensitivity of chaos to variations in parameters, the chaos generators of this kind are of interest from a practical point of view, for example in application for communication, random number generation, cryptographic schemes etc. [7].
D.S.A. acknowledges support of this work by RFBR grant № 14-02-31162, and S.P.K. thanks Prof. A. Pikovsky for discussions.
References
1. Anosov D.V., Gould G.G., Aranson S.K., Grines V.Z., Plykin R.V., Safonov A.V., Sataev E.A., Shlyachkov S.K, Solodov V.V., Starkov A.N., Stepin A.M. Dynamical Systems IX: Dynamical Systems with Hyperbolic Behaviour (Encyclopaedia of Mathematical Sciences) (Vol. 9). Springer, 1995.
2. Smale S. Differentiable dynamical systems // Bull. Amer. Math. Soc. (NS). 1967. Vol. 73. P. 747.
3. Williams R.F. Expanding attractors // Publications mathematiques de l'l.H.E.S. 1974. Vol. 43. P. 169.
4. Afraimovich V. and Hsu S.-B. Lectures on chaotic dynamical systems. AMS/IP Studies in Advanced Mathematics, Vol.28. American Mathematical Society, Providence RI, International Press, Somerville, MA, 2003.
5. Devaney R.L. An Introduction to Chaotic Dynamical Systems. NY: Addison -Wesley, 1989.
6. Shilnikov L. Mathematical problems of nonlinear dynamics: A tutorial // Int. J. of Bif. & Chaos. 1997. Vol. 7, 9. P. 1353.
7. Elhadj Z., Sprott J.C. Robust Chaos and Its Applications. WS, Singapore, 2011.
8. Banerjee S, Yorke J.A., Grebogi C. Robust Chaos // Phys. Rev. Lett. 1998. Vol. 80. P. 3049.
9. Hunt T.J., MacKay R.S. Anosov parameter values for the triple linkage and a physical system with a uniformly chaotic attractor // Nonlinearity. 2003. Vol. 16. P. 1499.
10. Hunt T.J. Low Dimensional Dynamics: Bifurcations of Cantori and Realisations of Uniform Hyperbolicity. PhD Thesis. University of Cambridge, 2000.
11. Belykh V, Belykh I., Mosekilde E. The hyperbolic Plykin attractor can exist in neuron models // Int. J. of Bifurcation and Chaos. 2005. Vol. 15. P. 3567.
12. Morales C.A. Lorenz attractor through saddle-node bifurcations// Ann. de l'Inst. Henri Poincare. 1996. Vol. 13. P. 589.
13. Kuznetsov S.P. Example of a physical system with a hyperbolic attractor of the Smale-Williams type // Phys. Rev. Lett. 2005. Vol. 95. 144101.
14. Kuznetsov S.P. and Seleznev E.P. A strange attractor of the Smale-Williams type in the chaotic dynamics of a physical system // JETP. 2006. Vol. 102. P. 355.
15. Isaeva O.B., Jalnine A.Yu. and Kuznetsov S.P. Arnold's cat map dynamics in a system of coupled nonautonomous van der Pol oscillators // Phys. Rev. E 74, 2006, 046207.
16. Kuznetsov S.P. and Pikovsky A. Autonomous coupled oscillators with hyperbolic strange attractors // Physica. 2007. Vol. D232. P. 87.
17. Kuznetsov S.P. Dynamical chaos and uniformly hyperbolic attractors: from mathematics to physics // Physics-Uspekhi. 2011. Vol. 54, № 2. P. 119.
18. Kuznetsov S.P. Hyperbolic Chaos: A Physicist's View. Higher Education Press: Beijing and Springer-Verlag: Berlin, Heidelberg, 2012.
19. Arzhanukhina D.S., Kuznetsov S.P. A system of three non-autonomous oscillators with hyperbolic chaos. I. The model with dynamics on attractor governed by Arnold's cat map on torus // Applied Nonlinear Dynamics (Saratov). 2012. Vol. 20, № 6. P. 56 (in Russian).
20. Arzhanukhina D.S., Kuznetsov S.P. A system of three non-autonomous oscillators with hyperbolic chaos. II. The model with DA-attractor // Applied Nonlinear Dynamics (Saratov). 2013. Vol. 21, № 2. P. 163 (in Russian).
21. Kuznetsov S.P. and Ponomarenko V.I. Realization of a strange attractor of the Smale-Williams type in a radiotechnical delay-fedback oscillator // Tech. Phys. Lett. 2008. Vol. 34, № 9. P. 771.
22. Kuznetsov S.P. and Pikovsky A.S. Hyperbolic chaos in the phase dynamics of a Q-switched oscillator with delayed nonlinear feedbacks // Europhysics Letters. 2008. Vol. 84. 10013.
23. Kuznetsov A.S., Kuznetsov S.P. Parametric generation of robust chaos with time-delayed feedback and modulated pump source // Communications in Nonlinear Science and Numerical Simulation. 2013. Vol. 18, № 3. P. 728.
24. Kuznetsov S.P. and Pikovsky A. Attractor of Smale-Williams type in an autonomous time-delay system. Preprint nlin. arXiv: 1011.5972, 2010.
25. Fowler A.C. An asymptotic analysis of the delayed logistic equation when the delay is large // IMA Journal of Applied Mathematics. 1982. Vol. 28. P. 41.
26. Benettin G., GalganiL., GiorgilliA., Strelcyn J.-M. Lyapunov characteristic exponents for smooth dynamical systems and for Hamiltonian systems: A method for computing all of them // Meccanica. 1980. Vol. 15. P. 9.
27. Schuster H.G. and Just W. Deterministic chaos: an introduction. Wiley-VCH, 2005.
28. Farmer J.D. Chaotic attractors of an infinite-dimensional dynamical system // Physica D. 1982. Vol. 4. P. 366.
29. Balyakin A.A. and Ryskin N.M. Features of computation of spectrum of Lyapunov exponents in distributed self-oscillatory systems with delayed feedback // Applied Nonlinear Dynamics (Saratov). 2007. Vol. 15, № 6. P. 3. (In Russian.)
30. Arnold V. I. and Avez A. Ergodic Problems in Classical Mechanics. New York, Benjamin, 1968.
31. Farmer J., Ott E., and Yorke J.The dimension of chaotic attractors // Physica D. 1983. Vol. 7. P. 153.
32. Grassberger P. and Procaccia I. Measuring the strangeness of strange attractors // Physica D. 1983. Vol. 9. P. 189.
33. Anishchenko V. S., Astakhov V. V., Neman A. B., Vadivasova T. E., and Schimansky-Geier L. Nonlinear Dynamics of Chaotic and Stochastic Systems: Tutorial and Modern Development. Springer, Berlin, Heidelberg, 2002.
Поступила в редакцию 17.03.2014
ГРУБЫЙ ХАОС В АВТОНОМНОЙ СИСТЕМЕ С ЗАПАЗДЫВАНИЕМ
Д. С. Аржанухина1, С. П. Кузнецов2'3
Саратовский государственный университет им. Н.Г. Чернышевского 2Институт радиотехники и электроники им. В.А. Котельникова РАН, Саратовский филиал 3Департамент физики и астрономии университета Потсдама
Рассматривается автономная система, построенная как модификации логистического дифференциального уравнения с запаздыванием и генерирующая последовательные цуги колебаний с фазой, трансформирующейся в соответствии с хаотическими отображениями. Система содержит две петли обратной связи, характеризующиеся двумя, вообще говоря, разными временами задержки. В случае их равенства хаотическая динамика определяется аттрактором Смейла-Вильямса, который соответствует двукратно растягивающему отображению окружности для фазы несущего сигнала цугов колебаний. При выбранных определенным образом разных временах задержки динамика фаз соответствует отображению Фибоначчи на торе. Таким образом, на аттракторе осуществляется динамика типа Аносова. В обоих случаях аттракторы проявляют грубость (отсутствие окон регулярности при изменении параметров) и, предположительно, относятся к классу структурно-устойчивых гиперболических аттракторов.
Ключевые слова: Аттрактор, гиперболический хаос, отображение, динамика Аносова, кот Арнольда, отображение Фибоначчи, аттрактор Смейла-Вильямса.
Аржанухина Дарья Сергеевна - родилась в Саратове (1988). Окончила факультет нелинейных процессов Саратовского государственного университета (2011). Защитила дипломную работу на тему «Сложная динамика и разрушение гиперболического хаоса в отображении «кот Арнольда» с диссипативным возмущением». В настоящее время является аспирантом базовой кафедры динамических систем факультета нелинейных процессов СГУ. Имеет 4 публикации.
410012 Саратов, ул. Астраханская, 83
Саратовский государственный университет им. Н.Г. Чернышевского E-mail: [email protected]
Кузнецов Сергей Петрович - родился в Москве (1951). Окончил Саратовский государственный университет (1973). Сотрудник Саратовского филиала Института радиотехники и электроники им. В.А. Котельникова РАН (с 1988). Руководит лабораторией теоретической нелинейной динамики. Защитил диссертацию на соискание ученой степени кандидата физико-математических наук в СГУ (1977) и доктора наук (1987) по специальности радиофизика. Профессор базовой кафедры динамических систем СГУ. Автор монографии «Динамический хаос и гиперболические аттракторы: От математики к физике», монографии «Hyperbolic Chaos. A Physicist's View», монографии «Strange Nonchaotic Attractors» (совместно с А. Пиковским и У. Фойдель), а также учебно-научных монографий «Динамический хаос» и «Нелинейные колебания» (в соавторстве с А.П. Кузнецовым и Н.М. Рыскиным). Опубликовал свыше 200 научных статей по нелинейной динамике, радиофизике и электронике. Под его руководством защищены десять кандидатских диссертаций. Лауреат государственной научной стипендии для ученых России (1994-1996), Соросовский доцент (1998), Соросовский профессор (2000, 2001). Член редакционной коллегии журнала «Известия вузов. Прикладная нелинейная динамика».
410019 Саратов, ул. Зеленая, 38
Саратовский филиал ИРЭ имени В.А. Котельникова РАН E-mail: [email protected]
