CC BY
14
УДК 519.21
Chaotic Dynamics of a Three-dimensional Endomorphism
Hacene Gharout* Nourredine Akroune^
Laboratoire des Mathématiques Appliquées Faculté des Sciences Exactes Universite de Bejaia, Bejaia, 06000
Algerie
Abelkadous Taha*
INSA, University of Toulouse 135 Avenue de Rangueil, Toulouse, 31077
France
Daniele-Fournier Prunaret§
Laboratory of LAAS-CNRS INSA, University of Toulouse 7 Avenue du Colonel Roche, Toulouse, 31077
France
Received 26.02.2018, received in revised form 09.07.2018, accepted 10.09.2018 The present work describes the phases plane bifurcations of some attractors given by a noninvertible three-dimensional map. This study is conducted through the critical manifolds concepts, generalization of critical points and critical lines introduced by Gumowski and Mira [1, 2]. The phase plane shared within two open regions: the first (denoted Z0) each point having no real preimage, and the second (denoted Z2) each point having two real preimages. The regions Z0, Z2 are separated by the critical manifolds, locus of points having two coincident preimages. This requires the visualization of critical manifolds in the three dimensional phases space. And this work also describes the passage of invariant or attractor curves towards weakly chaotic attractors then towards hyper-chaotic attractors via the contact bifurcation through the critical manifolds, which disappear after the contact bifurcation with the its attraction basin boundary.
Keywords: critical manifold, closed invariant curve, weakly chaos, chaos. DOI: 10.17516/1997-1397-2019-12-1-36-50.
1. Introduction and preliminaries
The aim of this work is to study the dynamic of a real three dimensional noninvertible nonlinear map and its phases space. Critical manifolds will be illustrated in the case of three dimensions, concepts which generalize the critical points and critical lines defined studied in dimension one and two [3-9]. The three dimensional noninvertible map is given by:
{Xn+1 = Vn,
Vn+1 = Zn, (1)
Zn+1 = Xn + avn(xn + Zn) + b.
where a and b are two real parameters.
§ [email protected] © Siberian Federal University. All rights reserved
The recurrence T defined in (1), admits two fixed points X1 = (x1,x1 ,x1) and X2 = = (x2, X2, X2), with:
= 1 + - 4b - 8ba =1 - - 4b - 8ba
X1 = 2(1 +2a) ' X2 = 2(1 +2a) '
defined for:
a = ÎT and b ^ 4(20~+Ty if a> ÎT5
1 -1
b > —-- if a < -.
4(2a +1) 2
For a = -21, the map T admits the unique fixed point, X* = (b, b, b).
The classical bifurcations obtained are:
.„„^ of F°u b;furcation ...........
Fold Bifurcation (Node-saddle): The equation of the Fold bifurcation curve, noted A11)0, is
-, with x = --— for a = -—.
4 +8a 2 +4a '2
Flip bifurcation (doubling period): The equation of the bifurcation Flip curve, noted Al, is , -31., -1 b =---a, with x = -.
4 2 2
Neimark-Sacker Bifurcation: The Neimark-Sacker bifurcation occurs when the modulus of the conjugated complex multipliers Si, i = 1, 2 for a stable fixed point or a cycle of order k of T, crosses value 1. The fixed point or the cycle becomes unstable, giving rise to a closed invariant curve which has the same stability as the fixed point or the original cycle. The obtained curves equations of Neimark-Sacker bifurcation are the following:
^ , -2 - 3a - 2a2 - 2a3 + 2^4 + 2a - 7a2 - 4a3 + 2a4 + 2a5 + a6 ri : b=-
r2 : b =
8 + 8a + 2a2
-2 - 3a - 2a2 - 2a3 - 2^4 + 2a - 7a2 - 4a3 + 2a4 + 2a5 + a6
8 + 8a + 2a2
In order to analyse the local stability of the fixed point, we need to evaluate the jacobian matrix at the fixed point and calculate the eigenvalues. Unfortunately the expressions defining the eigenvalues are so complicated that nothing can be said analytically. We can still say something through numerical simulations. After several numerical computations, we have found a generic two-dimensional bifurcation diagram in the (a, b) parameters plane, where the bifurcations curves (Fold (A(1)0, A21)0), Flip (A1) and Neimark-Sacker (r^ r2) bifurcations curves) delimit the zone of stability of the fixed point X2 of the map T in the parameters plane (a, b) (see Fig. 1, where the blue color is the X2 stability area).
b
2. Phases space of the map T
2.1. Critical manifolds of the map T
In the case of the considered transformation T, we note that the critical manifolds are surfaces in the three-dimensional space, and are parallel planes to one of the planes in space (x, y, z), when the parameter a vanishes. The critical curves, denoted EC, are surfaces divide the phases space in areas denoted Zi, i integer, where each Zi corresponds to the set of points of the phase
space that have a preimage of rank i for T. A critical range of rank k +1, denoted ECk is the manifold of rank k using EC0 (see equation (5)), k = 1, 2,... (ECk = Tk(EC0) for all k > 1). Representation of some critical manifolds for different values of a and b are presented in Fig. 2.
Fig. 1. The blue region is the stability zone of the fixed points of T
(a) (b)
Fig. 2. Representation of some critical manifolds for different values of a and b at (a, b) = = (—0.56, —0.75) in the figure (a) and at (a,b) = (—0.1, —0.9758) in the figure (b). Color codes are: EC-1 in brown color, EC in green, EC1 in blue, EC2 in gray, EC3 in cyan, EC5 in black, EC6 in white and EC8 in yellow
The phases space are divided into several regions by the critical manifolds, each of them has a different number of preimages of rank one. EC-1 is the antecedent of rank one of EC0. The equation for the critical manifold EC-1 for the map T satisfies | J(x, y,z) |= 0, where J(x, y, z) is the Jacobian of T at point (x, y, z):
i 0 10 \
J(x, y, z) = I 0 0 1 I . (2)
2x + ay a(x + z) ay The obtained critical manifold EC-1 depends only on y, z and the parameter a, with
x = — ay, y and z e R. (3)
Therefore, EC-t is given by:
EC-t = {( - ay,y,z), y,z G M}.
The manifold EC0 noted EC is defined by:
1 a2y2 4
ECo = T(EC-i) = { (y, z, 4a2y2 + ay( — 2ay + z) + b}, y,z G m}.
(4)
(5)
The critical manifolds or the critical surfaces (cf. [10]), in our cas, are planes for a equal to zero. A graphic illustration of the first eight critical planes is given in Fig. 3, where the critical manifolds are noted PCn instead of ECn, with
PC-i = {(0,y,z),y,z G M}, PCo = T(PC-i) = {(y,z,b), y,z G
(6) (7)
and PCn = Tn+l(PC-t), for n > 0.
Fig. 3. The first eight critical planes of T at a = 0 and b = —0.5
EC is a manifold separating the space R3 into two regions denoted Z0 and Z2. The region Z2 contains all the points that have two distinct antecedents of rank 1, and the region Z0 is formed by points without preimages. The corresponding recurrence is called of type Z0 — Z2 according to the Mira's classification [1].
The map T has two inverse functions corresponding to the two antecedents points of the region Z2, which are defined as follows:
Ti i(x ',y',z') :
T2 i(x \y',z') :
(8)
(9)
where £ =2^a2x'2 — 4ax ' y' + 4z ' — 4b.
These two functions are defined if a?x'2 — ax' y' + z ' — b ^ 0. The condition a2x'2 — ax' y'+ +z ' — b ^ 0 is verified for x' = y, y' = z and z ' ^ 4a?y2 + ay( ——ay + z) + b.
Hence the characterization of both regions Z0 and Z2 is given by:
Z0 = {(x, y, z) e R, a2x2 — 4axy + 4z — 4b < 0}, (10)
Z2 = {(x, y, z) e R, a2x2 — 4axy + 4z — 4b > 0}. (11)
Each of the two inverse determinations is associated with two layers defined in the zone Z2, which join along EC. We also define two unbounded regions R1 and R2 such that EC-1 = R1nR2 and for all X belonging to Z2, T-1(X) e R1 and T-1(X) e R2. When X e EC, we get T-1(X) = T2-1(X) belonging to EC-1 (Fig. 4).
Ti1
Fig. 4. Browsing T and the regions R1 and R2 after application of inverse determinations of T. EC and EC-1 are obtained at a = —0.6 and b = —0.75
2.2. Contact bifurcation
Denote by Q the stable fixed point X2, and r the invariant closed curve, consequent from the destabilization of the fixed point Q e R2, of type focus, via a Neimark-Sacker bifurcation. The curve is such that r C R2 before the bifurcation exposed here, that is r n EC-1 = 0 (see Fig. 5(a)).
The following properties hold (cf. [1] for particular example, [11] for the general case of Z0 — Z2 maps and [3] for two-dimensional noninvertible maps): Property 1. We have the following properties:
1. r n EC-1 = 0, r and its preimage of rank one T-1 (r) = r-1 are separate and do not intersect, that is r n r-1 = 0, with:
T-1(r) = T-1(r) U T2-1(r) = r-1 U r. (12)
2. The inner area bounded by r, denoted a(r) with Q e a(r), is invariant by direct iteration T(a(r)) = a(r); and reverse determination of T defined in a(r) is only given by T2-1. This characterizes the invariance of a(r) and r through the inverse iteration. Therefore, there is a subset U of r wherein T admit a unique inverse defined by T-1 : U —> U. This is the local invariance iteration.
For two values of the parameters a = a1 and b = b1, a contact bifurcation appears when the invariant curve r comes into contact with EC-1 at a point noted u0 (r n EC-1 = {u0}); moreover, the image u1 of u0 by T is a point of tangential contact between r and EC. For small values e1, e2 such that a = a1 + e1 and b = b1 + e2, we have:
r n EC-1 = {u0,v0} and r n R1 = 0, (13)
(a) a = 0.25 (b) a = 0.446
Fig. 5. Before the bifurcation, there is no contact between the closed invariant curve r and preimage of rank one r_i; and after the bifurcation r intersects r_i through EC_1 and comes into tangential contact with EC at b = —0.75 with X0 = (0.1,0.1,0.1)).
which means that r intersects EC_1 in two points u0 and v0; and invariant curve touches the zone R1. This intersection causes the changes of the curve r because the iterates of the rank n for the couple {uo,vo} noted (un,vnjn^1 are points of contact between r and ECn_1. Indeed, r cannot pass through all critical manifolds |Jm=_11 ECi_1 (all critical manifolds delimiting r, except EC_1), which means in the opposite case that a point of Z0 has an antecedent. This bifurcation creates oscillations of r along the ECm. Before the bifurcation, r is invariant by the inverse iteration T-1. After the bifurcation, the arc (u1,v1) of r stays in r, neither by the inverse iteration T_1 but by T1_1. The curve r changes the sheet on u1 and v1 at the contact with EC creating oscillations in the projected three-dimensional space (see Fig. 6).
Fig. 6. After the bifurcation r intersects r_1 in two points (u0,v0) through EC_1 and comes into tangential contact with EC in two points (u1,v1) at a = 0.7 and b = —0.75 with X0 = = (0.1,0.1,0.1)
Property 2. Suppose that r n EC_1 = 0 and r n R1 = 0, then we have the following properties:
1. The area a(T) is no longer invariant by T: T(a(r)) = a(T).
2. r and r_i intersect at least on two points u0 and v0 of EC_i. The part r C R1 is folded and becomes tangent to the critical manifold EC by direct iteration T.
3. r is not invariant by inverse iteration. Indeed, the restricted inverse Tr for the map T in r is defined now by T1_1 and T-1 with respect to the considered part of r and T-1 (r) = r.
4. There exists a neighborhood U(r) of r, for which T do not have a unique inverse from U to U; it is defined by the two inverse determinations T-1 and T—1. This characterizes the local non invariance of T by the inverse iteration.
(a) (b)
Fig. 7. Enlargement of the closed invariant curve at a = 1.01 and b = —0.75 with X0 = = (0.1, 0.1,0.1)
For a = a2 and b = b2, the closed curve r becomes tangent to EC (see Fig. 7) and the some property of tangency occur in all the images of EC (see Fig. 8), with the new change of sheets for r in a point that a neighborhood U verifies the fourth characterization of the Property 2.
Fig. 8. Tangency of the closed invariant curve r (in red color) with the critical surfaces EC0(green), EC1 (blue) and EC2(gray) at a = 1.01 and b = —0.75 with X0 = (0.1,0.1,0.1)
The bifurcation creating self-intersections on the invariant curve r (points of rnr_i G EC_1) causes the transformation of a closed invariant curve into a chaotic attractor. Indeed, for an closed invariant curve, the rotation's number associated is irrational, therefore there exist infinity points of self-intersections reproduced by self-similarity and successive iterations of T, which means that r becomes fractal and its dimension DL in the sense of Lyapounov becomes higher (see Fig. 9(a)). Thus, the formed set is called weakly chaotic and by increasing parameters the complexity of the invariant curve r grows. This results in an enlargement of the weakly chaotic attractor and increased Lyapounov exponents, as well as the associated DL dimension and this generates a chaotic attractor (see Fig. 9).
(c) a=1.14 (d) a=1.14
Fig. 9. Creation of the chaotic attractor at b = —0.75 with X0 = (0.1,0.1,0.1). (a) Chaotic attractor in the space at a = 1.12 and b = —0.75 and the corresponding Lyapounov dimension Dl = 2.12932. (b) The attractor intersection with its antecedent of rank one, through EC_1 and into tangential contact with EC. (c) Chaotic attractor in the space at a = 1.14 and b = —0.75 and Dl = 2.19071). (d) Chaotic attractor limited by the critical surfaces
3. Study of T with b fixed and a variable
As we have already mentioned above, the transformation T defined in (1) is of type Z0 — Z2, has two fixed points Xi and X2, saddle focus and a stable node focus respectively for the values of a and b defined by the Neimark-Sacker, Flip and Fold bifurcation curves. If we vary the parameter
a in the positive direction with b < —0.75, the stable fixed point of map T becomes unstable giving rise to an invariant closed curve via a Neimark-Sacker bifurcation for (a, b) = (0, —0.75). However, if a varies in the negative direction, the birth of two invariant closed orbits is observed. The existence of two dynamics is clearly illustrated by the bifurcation diagram (Fig. 10(a)) for different values of a and b = —0.75 , showing the appearance of double orbits (negative sense of the variations of a) and orbits of order one (positive direction of the variations of a) are illustrated by the projection on the plane of the phases (x,y), as well as the magnification given in Fig. 10(b). For (a, b) = (0, —0.75) the transformation T has a fixed point of order one, which is the only fixed point of order two of T and which belongs to the Neimark-Sacker bifurcation curve (X = (—0.5, —0.5, —0.5) of multipliers S1 = S2 = 0.5 ^ 0.8660254038/, S3 = —1 and of modulus p =1).
■0.6128 0 1.4613
a
(a) Bifurcation diagram (b = —0.85)
(b) Projection in the plane (x, y) (c) Enlargement of the projection
Fig. 10. Bifurcation diagram at b = —0.75 and his projection in the phase plane (x,y)
Indeed, for b = —0.75 and by increasing a, the closed invariant curve r comes into contact with the critical manifolds for the value a = 0.446 (see Fig. 5(a)). Thus, oscillations of the invariant curve appear along the critical manifolds (see Figs. 6 and 7); for a = 1.01 and b = —0.75 the self-intersections generate loops on r, which becomes a weakly chaotic attractor.
By increasing the parameter a, loops are formed and the set of points such that r n r_1 (/ EC_1 grows. The complexity of r increases, and becomes a chaotic attractor (see Fig. 9(a) and 9(b)). We give in Fig. 9(c) the chaotic attractor obtained for a = 1.14 and b = —0.75 and some critical manifolds which delimited and enveloped it (see Fig. 9(d)). The corresponding Lyapounov dimensions DL(a) are DL(1.01) = 1.44203, DL(1.12) = 2.12932 and DL(1.14) = 2.19071. By increasing the parameter a, the attractor comes into contact with the border of its attraction basin which disappears for a = 1.461391 (see Fig. 11).
The two fixed points of the map T2 are stable for b = —0.75 and —0.484 ^ a < 0 becomes
unstable for a < —0.484 and through a Neimark-Sacker bifurcation, to two closed invariant curves r and r2 appears in the zone Z2, such that:
r = ri U r and T-1(r) u T-1(r) = r_ U r^. (14)
For a = —0.5 and b = —0.75, the closed invariant curve r1 comes into contact with EC_1 and intersects its antecedent of rank one r-1 on two points uQ and vQ, once the bifurcation contact occurred; and the invariant curve r2 comes into a tangential contact with EC on two points u1 and v1, which are the images of uQ and vQ by direct application of map T (see Fig. 12(a)):
r n EC_1 = {uQ, vq}. (15)
Varying the parameter a, even in the negative direction the invariant closed curves r1 and r2 change the form (each curve is folded and doubled). From then on, the closed invariant curve r1 comes into contact with r-1 on four points uQ, vQ, rQ and sQ belonging to EC_ 1; and r2 contacts EC in four points u1, v1, r1 and s1 which are the images of the points uQ, vQ, rQ and sQ by T (see Fig. 12(b)).
Fig. 11. The chaotic attractor (obtained at a = 1.14 and b = —0.75 in color green) is going to undergo a contact bifurcation with its attraction basin boundary (in color cyan), which will make it disappear and which occurs at a = 1.461391
For a = —0.556 and b = —0.75, the self-intersections generate loops on r (r = r1 U r2), which becomes a weakly chaotic attractor, with the Lyapounov dimension DL = 2.00322 (see Fig. 13(a)). Furthermore, if we reduce the parameter a, loops are formed and the set of points such that r n r_1 ^ EC_1 grows. Then, the two cyclic attractors expand and merge into a single attractor. The complexity of r and the Lyapounov dimension increase (DL(a): DL(— 0.56) = 2.04589, Dl(—0.6) = 2.65599 and Dl(—0.61) = 3) (see Figs. 13(b) and 13(c)). Hence, the attractor becomes a chaotic. By increasing the parameter a, the attractor touches the border of his attraction basin and disappears for a = —0.612 (see Fig. 14).
4. Study of T with a fixed and b variable
The coexistence of two attractors is observed when the parameter a is fixed and b is varied. The coexistence of two invariant curves for the map T is illustrated by the Fig. 15 with a = —0.1 and b varied (b = —1.25 and b = —1.33). We showing that the two coexisting invariant curves in Fig. 15(a)), becomes a weakly chaotic attractors (Fig. 15(b)) after the bifurcation, before the merge in a single hyper-chaotic attractor at a = —0.1 and b = —1.38 (see Fig. 16(d)).
Fig. 12. Bifurcation contacts of two invariant closed curves with EC—i (b = —0.75, X0 = (0.1,0.1,0.1)). (a) r intersects r— on two points u0 and v0 through EC-i and r2 comes into a tangential contact with EC on two points ui and vi (at a = —0.5 and b = —0.75). (b) Folding and doubling of the two invariant closed curves, ri intersects r-i on four points u0, v0, r0 and so belonging to EC—i and r2 contacts EC in four points ui, vi, ri and si (at a = —0.546 and b = —0.75)
Fig. 13. Evolution of the two closed invariant curves into a chaotic attractor at b = —0.75 with X0 = (0.1,0.1, 0.1)
• Suppose that X0 = (0.1,0.1,0.1); the stable fixed point of the map T becomes unstable and gives rise to an invariant closed curve via a Neimark-Sacker bifurcation (b « —0.7732418029), which turns into a chaotic attractor via the bifurcations of contacts with the critical manifolds and the antecedent of rank one of this attractor. By decreasing the parameter b, oscillations are created on r and Lyapounov dimension DL grows (e.g. DL( —1.08) = 2.00245), and an weakly chaotic attractor appears for b = —1.1 (see Fig. 17).
If we still decreasing the parameter b, the chaotic attractor disappears for b « —1.1058, and an attractor of six parts appears (see Fig. 16(a)). Furthermore, the attractor which is born for b = —1.1058 becomes weakly chaotic and the points number provided from the intersection of this attractor with his antecedent of rank one outside of EC-\ increases, and its complexity also, as well as the Lyapounov dimension (DL(b) : DL( —1.2) = 1.15288, DL( — 1.26) = 1.47679 and Dl ( — 1.3) = 2.08497). The attractor becomes chaotic and Lyapounov dimension DL reaches the
Fig. 14. The chaotic attractor (in color green) is going to undergo a contact bifurcation with its basin boundary (in cyan), which will make it disappear (figure for a = —0.6 and b = —0.75 before the contact bifurcation) and which occurs at a = —0.612
Fig. 15. Coexistence of two attractors. (a) Two invariant curves coexist for the map T at a = —0.1 and b = —1.25. (b) The two invariant curves become chaotic attractors at a = —0.1 and b = —1.33
value of 3 and disappears for b < —1.57879 (see Fig. 16).
• By choosing X0 = (1.5,1, —0.5), the stable cycle with the order two of the map T for — 1.155 < b < —0.7, becomes unstable and gives rise to an invariant curve formed by two closed orbits (see Fig. 18(a)), which evolves to the chaotic attractor previously obtained (see Fig. 16(d)), through the contact bifurcations with the critical manifolds. Here we give some figures which show the qualitative changes of the invariant curve. By decreasing the parameter b, the invariant curve comes into contact with its antecedent of rank one with crossing the EC_i (see Fig. 18(b)), and oscillations are created.
By still decreasing b, a weakly chaotic attractor is formed for b = —1.27, and expands for b = —1.34 and DL( —1.34) = 3 (see Fig. 18(c)); then the two parts of the attractor merge and form the single chaotic attractor obtained at b = —1.38 (see Fig. 16(d)), whose Lyapounov dimension is DL( —1.38) = 3. The hyper chaotic attractor touches the border of his attraction basin and disappears at b = —1.5788. The attraction basins are illustrated for the attractors obtained at b = —1.34 and b = —1.38 with a = —0.1 in Fig. 19.
The figures of attraction basins and the Lyapounov dimension are obtained using the Dynamics logiciel [12].
(c) b = -1.1 (d) b = -1.1
Fig. 16. The second chaotic dynamics for a = —0.1 and X0 = (0.1,0.1,0.1)
(a) b = -1. (b) b = -1.08. (c) b = -1.1.
Fig. 17. The dynamic evolution of a closed invariant curve to chaotic attractor (a = —0.1, X0 = (0.1,0.1, 0.1))
Conclusion
In this work, by setting one of the parameters (a and b) and varying the other parameter of the endomorphism T, the dynamic evolutions of some attractors have been illustrated in the phase space. The passage from invariant curves or attractors to weakly chaotic attractors then to hyper-chaotic attractors via contact bifurcations with critical varieties (critical surfaces) in the three-dimensional phase space, which disappear after the bifurcation of contact with the limit of their basins of attraction. The coexistence of invariant orbits or attractors is observed as well as their dynamic evolutions towards a single chaotic attractor.
(a) b = -1.19. (b) b = -1.27 (c) b = -1.34.
Fig. 18. Evolution of two invariant closed curves to a weakly chaotic attractor, then the two parts of the attractor merge and form a single hyper chaotic attractor (obtained at b = —1.38 and illustrated in Fig. 16(d))); at a = —0.1 with Xq = (1.5,1, —0.5)
(a) (b)
Fig. 19. (a) The attraction basin (in color cyan) of the weakly chaotic attractor obtained at a = —0.1 and b = —1.34 with X0 = (1.5,1, —0.5). (b) The chaotic attractor (at a = —0.1 and b = —1.38 in color green) and its attraction basin (in cyan), which disappears after the contact bifurcation with the its basin boundary which occurs at b = —1.5788
References
[1] I.Gumowski, C.Mira, Dynamique chaotique Transformatios ponctuelles Transition Ordre-Desordre, Cepadues Editions, Toulouse, 1980.
[2] C.Mira, Chaotic Dynamics: From the One-dimensional Endomorphism to the Two-Dimensional Diffeomorphism, World Scientific Publishing, Singapore, 1987.
[3] C.E.Frouzakis, L.Gardini, I.G.Kevrekidis, G.Millerioux, C.Mira, On some proporties of invariant sets of two-dimensional noninvertible maps, International Journal of Bifurcation and Chaos, 7(1997), no. 7, 1167-1194.
[4] C.Mira, A.Agliari, L.Gardini, From the box-within-a-box bifurcation structure to the julia set. Part II: Bifurcation routes to different julia sets from an indirect embedding of a quadratic complex map, International Journal of Bifurcation and Chaos, 19(2009), no. 10, 3235-3282.
[5] G.Millerioux, C.Mira, Noninvertible piecewise linear maps applied to chaos synchronization and secure communications, International Journal of Bifurcation and Chaos, 7(1997), no. 7, 1617-1634.
[6] E.Mosekilde, Y.Maistrenko, D.Postnov, Chaotic Synchronization: Applications to living systems, World Scientific Series on Nonlinear science, Series A, vol. 42, Singapore, 2001.
[7] J.Xu, D.Fournier-Prunaret, A.K.Taha, P.Charge, Chaos generator for secure transmission using a sine map and a RLC series circuit, Science China, 1(2010), no. 53, 129-136.
[8] J.C.Sprott, Simplest dissipative chaotic flow, Physics Letters A, 228(1997), 271-274.
[9] G.M.Zaslavsky, Chaos in Dynamic Systems, Harwood, Amsterdam, 1985.
[10] G.I.Bischi, L.Mroz, H.Hauser, Studying basin bifurcations in nonlinear triopoly games by using 3D visualization, Nonlinear Analysis, Theory, Methods & Applications, 47(2001), no. 8, 5325-5341.
[11] C.Mira, L.Gardini, A.Barugola, J.C.Cathala, Chaotic dynamics in two-dimensional nonin-vertible maps, World Scientific Publishing. Serie A, Singapore, 1996.
[12] H.E.Nusse, J.A.Yorke, Dynamics: Numerical Explorations, Springer-Verlag, New York, 1998.
Хаотическая динамика трехмерного эндоморфизма
Хасен Гаро Ноуредин Акроун
Факультет точных наук Университет Бежа, Бежа, 06000 Алжир
Абелкадоус Таха
Национальный институт прикладных наук Тулузы Рангуил, 135, Тулуза, 31077 Франция
Даниел-Фоурниер Прунаре
Национальный институт прикладных наук Тулузы Полковник Рош-авеню, 7, Тулуза, 31077
Франция
В настоящей работе описываются бифуркации фазовых плоскостей некоторых аттракторов, задаваемых необратимым трехмерным отображением. Это исследование проводится с помощью концепций критических многообразий, обобщения критических точек и критических линий, введенных Гумовским и Мирой [1, 2]. Фазовая плоскость делится в двух открытых областях: первая (обозначается Z0) каждая точка, не имеющая реального прообраза, а вторая (обозначенная Z2) каждая точка имеет два реальных прообраза. Области Z0, Z2 разделены критическими многообразиями, локус точек, имеющих два совпадающих прообраза. Для этого требуется визуализация критических многообразий в пространстве трехмерных фаз. Работа также описывает прохождение инвариантных или аттракторных кривых в сторону слабохаотических аттракторов, а затем к гиперхаотическим аттракторам через контактную бифуркацию, через критические многообразия, которые исчезают после контактной бифуркации с ее границей притяжения.
Ключевые слова: критическое многообразие, замкнутая инвариантная кривая, слабый хаос, хаос.
