354
Views
1
CrossRef citations to date
0
Altmetric
Articles

Markov partition in the attractor of Lozi maps

ORCID Icon
Pages 952-964 | Received 15 Apr 2022, Accepted 20 Feb 2023, Published online: 13 Mar 2023

Abstract

In this paper, we study iterations of two-dimensional maps, in particular iterations of Lozi maps in the region of the parameter space where it has a strange attractor. Using symbolic dynamics techniques for two-dimensional maps, based on the kneading theory of Milnor and Thurston and also in the symbolic dynamic formalism developed by Sousa Ramos, through the kneading sequence for the Lozi maps, we characterize the region in the parameter space that contains the kneading curves and present a method to define a Markov partition for the Lozi attractors. Consequently, the topological entropy for the Lozi map is computed.

Mathematics Subject Classifications:

1. Introduction

Consider the following parametrized family of Lozi maps: (1) Lab(x,y)=(1a|x|+y,bx),(x,y)R2.(1) These piecewise linear plane invertible maps were introduced in 1978 by René Lozi  [Citation5] as a simplified version of the parametrized family of Hénon maps, without loosing the capacity to show chaotic behaviour. Since then, Lozi maps has been an important tool to understand complex behaviour of iterated maps on the plane.

Since it is an invertible map, it is possible to define its inverse, Lab1, given by Lab1(x,y)=(yb,x1+a|yb|).In 1992  [Citation6] Liu et al. showed that for parameters (a,b) satisfying (2) {0<b<1;a>b+1;a>2b/2,(2) every map Lab has indeed a strange attractor. Let L denote the corresponding region of the parameter space; see Figure .

Figure 1. The parameter space region for which Lab has a strange attractor.

Figure 1. The parameter space region for which Lab has a strange attractor.
From this point the Lozi maps considered will have parameter values a and b belonging to L. For simplification purposes, Lab will be denoted by L.

Earlier, Michał Misiurewicz  [Citation7] introduced a region L in R2 satisfying (3) L(L)L,(3) allowing to present the Lozi map attractor, AL, from its successive forward iteration, (4) AL=n=0Ln(L).(4) Considering the fixed point that lies in the first quadrant A=(1a+1b,ba+1b)ifb<a+1,and computing its local stability by evaluating the eigenvalues of the jacobian matrix of the map at the fixed point, in the domain of the parameters for a>|b|+1, the fixed point is a saddle point. Let Ls(A) and Lu(A) be the stable and unstable manifolds of A, respectively.

Definition 1.1

The Lozi map fundamental domain L is the triangle with vertices I, L(I), and L2(I), where (5) I=(2+a+a2+4b2(1+ab),0)(5) is the intersection point of the unstable manifold of the fixed point A, LAu, with the horizontal axis.

Next, some definitions of Lozi maps symbolic dynamics are presented, following [Citation4].

Definition 1.2

Given XR, the itinerary of X is the bi-infinite symbolic sequence (6) it(X)=δ2δ1δ0δ1δ2,(6) where the symbols δn, for nZ, are defined as (7) δn{1,Ln(X)x<0,,Ln(X)x=0,+1,Ln(X)x>0,(7) where Xx corresponds to the x component of X.

By the previous definition, it(X) is a bi-infinite symbolic sequence of {1,,+1}Z.  It is important to distinguish the itinerary for the past and the itinerary for the future of a point X; let δ_(X)=δ2δ1δ0δ1δ2 be an itinerary of X. The symbolic subsequence δ_u(X)=δ2δ1 is the itinerary for the past of X and the symbolic subsequence δ_s(X)=δ0δ1δ2 is the itinerary for the future. Thus, δ_(X)=δ_u(X)δ_s(X).

2. Critical set

Let V be the vertical axis, such that V=V+1V1,where V+1 is the nonnegative part of the vertical axis and V1 is the negative part of the vertical axis.

According to the construction process suggested by Ishii  [Citation4] and Baptista et al. [Citation3], some points belonging to ALVplay a very important role in the study of Lozi's maps dynamics, especially the point L1(I). The orbit associated with the forward itinerary of the point I is the orbit whose itinerary, for the future, corresponds to the maximum symbolic sequence within the set of all admissible symbolic sequences for the future.

Definition 2.1

Given a Lozi map L, its kneading sequence K(L) is the itinerary for the future of I, i.e. K(L)=δ_s(I).

Definition 2.2

Let (a,b)L. Given a Lozi map L, its pruned kneading sequence, kn(L), is the initial subsequence of symbols of K(L) such that the first symbol ★ appears at position n, that is, Ln1(I)x=0.

From the inverse map L1, two maps can be introduced, L+11 and L11, as follows: L+11(x,y)=(yb,x1+ayb)and L11(x,y)=(yb,x1ayb).Using these two maps and a pruned kneading sequence, kn(L), we will define a certain type of line segments, in L, as follows:

Definition 2.3

Let (a,b)L and δ_s=kn(L)=δ0δ1δn2 for n>3.

The line segments Sδiδn2, with 2in2, in L, are defined by Sδn2=Lδn21(LVδn2)L;Sδn3δn2=Lδn31(Sδn2)L;Sδ2δn2=Lδ21(Sδ3δn2)L.

After the definition of the line segments Sδiδn2 with 2in2, for n>3, its possible to define the critical set for L when (a,b)L.

Definition 2.4

Let (a,b)L and δ_s=kn(L). The critical set Cδ_s(L) is defined as the set containing the line segments Sδiδn2, with 2in2, the vertical line segment S=VL and the points L(I) and I.

The following proposition presents the slope of each line segment Sδiδn2, with 2in2 (for n>3).

Proposition 2.5

Let (a,b)L and δ_s=kn(L) =δ0δ1δn2. The line segments Sδiδn2, with 0in2, are given by Sδiδn2=L{(x,y):y=mi(xxi)+yi},where (xi,yi)=Li(I) and mi={aδn2,i=n2,aδi+bmi+1,0i<n2.

Proof.

Let δ_s=kn(L)=δ0δ1δn2 be the pruned kneading sequence of L with dimension n.

Consider the jacobian matrix of L1 at the point XR, defined as follows: DL+11=[01b1ab]ifXy>0and DL11=[01b1ab]ifXy<0,where Xy corresponds to the y component of X.

Let v=(0,1) be the vector of the line segment Vδn2. The vector of a line segment Lδn21(Vδn2)L is given by vδn2=v(DLδn21)T=[01][011baδn2b]=[1babδn2].Depending on the value of δn2, the point Ln1(I) will belong to V+1 or V1. By the definition of Lδn21, Sδn2=Lδn21(Vδn2)L will be a line segment that passes through the point Ln2(I), with vector vδn2=(1b,abδn2) and slope mn2=aδn2.

Let vδn2=(vx,vy). As [vxvy][011baδn3b]=[1bvyvx+abvyδn3],the slope of the line segment Sδn3δn2=Lδn31(Sδn2)L will be mn3=bvxvy+aδn3=aδn3+bmn2.Therefore, for 2in2, the line segment Sδiδn2 that passes through the point (xi,yi)=Li(I) has the vector given by vδiδn2=vδi+1δn2(DLδi1)Tand slope given by mi=aδi+bmi+1.

Proposition 2.6

Let (a,b)L, δ_s=kn(L). The intersection of the line segment Sδiδn2 with the x-axis is the point (ξi,0) where ξi={0ifi=n1.1mi(1ξi+1)if2in2.

Proof.

The line segment S (when i=n1) is the vertical line, and its intersection with the x-axis occurs at the point (ξn1,0)=(0,0). The pre-image of (ξn1,0) is Lδn21(0,0)=(1,0). Let mi be the slope of the line segments Sδiδn2 with 2in2, for n>3. Therefore, the line segment that passes through the points (0,1) and Ln2(I) has slope mn2 and can be given by y=mn2x1. Since this line segment contains the segment Sδn2, the intersection of the line segment Sδn2 with the x-axis occurs at the point (ξn2,0)=(1mn2,0).

Continuing this process, the pre-image of the point (ξn2,0) is the point (0,ξn21). Therefore, the line segment Sδn3δn2 is contained in a line segment that passes through the point (0,ξn21) with slope mn3 and intersects the x-axis at the point (ξn3,0)=(1mn3(1ξn2),0) . In general, the pre-image of a point (ξ,0) is the point (0,ξ1) and consequently the line segment Sδiδn2 with 2in2, for n>3, intersects the x-axis at the point (ξi,0), where ξi=1mi(1ξi+1).

Considering Propositions 2.5 and 2.6, the next result follows.

Proposition 2.7

Let (a,b)L and δ_s=kn(L). The line segments Sδiδn2, with 2in2, for n>3, and the vertical line segment S do not intersect.

2.1. Order relation in the set {1,,+1}N

From [Citation4], the attractor points with the same fixed itinerary for the past ε_u lie in a segment Lu, called the unstable leaf, completely characterized by the parameter values (a,b) and ε_u. Moreover, attractor points with the same fixed itinerary for the future ε_s lie in a segment Ls, called the stable leaf, also completely characterized by (a,b) and ε_s. If a point XAL has a fixed itinerary ε_, X is the intersection of the line segments previously defined, denoted by Lu(X) and Ls(X), respectively.

In order to rewrite the order relation established by Ishii [Citation4] for the set {1,,+1}N, consider the following definition.

Definition 2.8

Let ε_s and δ_s be two symbolic subsequences for the future.

ε_s<sδ_s if one of the following conditions is satisfied:

  1. ε0<δ0;

  2. ε0εi1=δ0δi1 and εiδi,:

  1. the number of +1s in ε0εi1 is even and εi<δi;

  2. the number of +1s in ε0εi1 is odd and εi>δi;

where the order on the symbols is 1<<+1.

If εi=δi for all i0, ε_s=δ_s.

From the previous definition it can be establish that, for any point XL, (8) δ_s(L(I))<sδ_s(X)<sδ_s(I).(8) Considering the set {σi(δ_s)}i=0n1, where σis the shift-operator, the set {ξi}i=2n2{Ix,L(I)x} and denoting by p as a permutation in the set {2,,n1} such that (9) σ(δ_s)<sσp(2)(δ_s)<s<sσp(n1)(δ_s)<sδ_s(9) it follows that (10) L(I)xSσp(2)(δ_s)Sσp(n1)(δ_s)Ix,(10) and where AB means that A is on the left side of B, and consequently (11) L(I)x<ξp(2)<<ξp(n1)<Ix.(11)  Figure  illustrates the critical set, the attractor and the invariant triangle L, considering the pruned kneading sequence δ_s =+11+1+1.

Figure 2. The critical set and the Lozi map attractor for the pruned kneading sequence δs=+11+1+1.

Figure 2. The critical set and the Lozi map attractor for the pruned kneading sequence δs=+1−1+1+1∗.

2.2. Markov partition on AL

Using the elements of the critical set Cδ_s(L), where δ_s=kn(L), it is possible to construct a partition, with n1 regions, for the set L and therefore to AL.

Considering the Definition 2.8 and the order relations given in (Equation9), (Equation10) and (Equation11) we define the n1 regions on the triangle L as follows:

let p be a permutation in the set {2,,n1} such that L(I)xSσp(2)(δ_s)Sσp(n1)(δ_s)Ix.The n1 regions on the triangle L are defined as follows: R1={(x,y)L:(x,y)Sσp(2)(δ_s)}Rk={(x,y)L:Sσp(k)(δ_s)_(x,y)Sσp(k+1)(δ_s)},k=3,,n2Rn1={(x,y)L:Sσp(n1)(δ_s)_(x,y)},where A_(x,y) means that (x,y) is on the right side of A or (x,y) belongs to A and (x,y)A means that (x,y) is on the left side of A.

Considering the previous construction the next result follows:

Theorem 2.9

The set P={R1,,Rn1} is a topological partition of L.

Let X belong to the interior of a set Ri. Suppose that L(Ri)Rj,where Rj is single region or a union of two or more regions of the topological partition P. Assume, without loss of generality, that it is a single region. Thus, L(X)Rj.Moreover, by Ishii [Citation4], the unstable leaf and the stable leaf of X lie in a segment denoted by Lu(X) and Ls(X), respectively, and X is the intersection point of the segments previously defined. Therefore, and due to the Lozi map linearity, we can state that L(Lu(X)Ri)Lu(L(X))Rjand L(Ls(X)Ri)Ls(L(X))Rj.Thus, by Bowen [Citation2], it follows:

Theorem 2.10

The topological partition {R1,,Rn1} provides a Markov partition of L associated to δ_s=kn(L).

3. Transition matrix and topological entropy

Definition 3.1

Let {R1,,Rn1} be a Markov partition of L associated to δ_s=kn(L). The transition matrix over L is a square matrix Aδ_s, of order n1, such that Aδ_s=(aij)={1ifL(Ri)Rj0ifL(Ri)Rj=.

Following [Citation8] and [Citation1], the topological entropy of L can be calculated using the corresponding transition matrix. This result can be stated as follows:

Proposition 3.2

Let δ_s=kn(L)=δ0δ1δn2 and (a,b)L. Let Aδ_s be the transition matrix associated to δ_s. The topological entropy of L is given by htop(L)=log(λmax(Aδ_s)),where λmax(Aδ_s) is the spectral radius of Aδ_s.

Example 3.3

Let consider (a,b)(1.518,0.185728) and the pruned kneading sequence associated δ_s=+11+1+1The itinerary given by the pruned kneading sequence δ_ for (a,b)(1.518,0.185728) is given in Figure .

Figure 3. Itinerary given by the pruned kneading sequence and Lozi attractor for (a,b)(1.518,0.185728).

Figure 3. Itinerary given by the pruned kneading sequence and Lozi attractor for (a,b)≈(1.518,0.185728).

Applying the shift map to the pruned kneading sequence δ_s, it follows: δ_s=+11+1+1σ(δ_s)=1+1+1σ2(δ_s)=+1+1σ3(δ_s)=+1σ4(δ_s)=By Definition 2.8, the following order relation is valid σ(δ_s)<sσ4(δ_s)<sσ2(δ_s)<sσ3(δ_s)<sδ_sand consequently L(I)xSS+1+1S+1Ix.The partition of L is given through the regions () R1={(x,y)L:(x,y)S}R2={(x,y)L:S_(x,y)S+1+1}R3={(x,y)L:S+1+1_(x,y)S+1}R4={(x,y)L:S+1_(x,y)}which can be associated to the following transition matrix Aδ_s given by Aδ_s=[0011000101101000].The transition matrix Aδ_s has maximum eigenvalue λmax1.51288 and the topological entropy is approximately htoplog(1.51288)0.41402.

Example 3.4

Consider (a,b)(1.887,0.0957356) and the pruned kneading sequence associated δ_s=+11111+11.Applying the shift map to δ_s δ_s=+11111+11σ(δ_s)=1111+11σ2(δ_s)=111+11σ3(δ_s)=11+11σ4(δ_s)=1+11σ5(δ_s)=+11σ6(δ_s)=1σ7(δ_s)=that verifies the following order relationship: σ(δ_s)<sσ2(δ_s)<sσ3(δ_s)<sσ6(δ_s)<sσ4(δ_s)<sσ7(δ_s)<sσ5(δ_s)<sδ_sand L(I)xSσ2(δ_s)Sσ3(δ_s)Sσ6(δ_s)Sσ4(δ_s)Sσ7(δ_s)Sσ5(δ_s)Ix.The partition of L is given through the regions () R1={(x,y)L:(x,y)Sσ2(δ_s)},R2={(x,y)L:Sσ2(δ_s)_(x,y)Sσ3(δ_s)},R3={(x,y)L:Sσ3(δ_s)_(x,y)Sσ6(δ_s)},R4={(x,y)L:Sσ6(δ_s)_(x,y)Sσ4(δ_s)},R5={(x,y)L:Sσ4(δ_s)_(x,y)Sσ7(δ_s)},R6={(x,y)L:Sσ7(δ_s)_(x,y)Sσ5(δ_s)},R7={(x,y)L:Sσ5(δ_)_(x,y)},which can be associated to the following transition matrix Aδ_s: Aδ_s=[0100000001100000001000000010000000100011111110000].The transition matrix Aδ_s has a maximum eigenvalue λmax1.9363 and the topological entropy is approximately htop=log(1.9363)0.66078.

Figure 4. Itinerary of point I, the critical set and the regions Rj, j = 1, 2, 3, 4.

Figure 4. Itinerary of point I, the critical set and the regions Rj, j = 1, 2, 3, 4.

3.1. Isentropics curves for Lozi maps

From [Citation3], a pruned kneading sequence corresponds to an algebraic condition on the parameters and, therefore, to a curve on the parameter plane. For example, a simple computation shows that the kneading curve for the pruned kneading sequence δ_s=+11 is described by the following condition: 2+4aa3aba2a2+4b+ba2+4b=0.

Definition 3.5

Given a pruned kneading sequence δ_s of dimension n, its kneading curve Γ(δ_s) is the parameter space curve corresponding to parameter values such that kn(L)=δ_s.

Following the work [Citation3], and considering Proposition 3.2, the topological entropy of L along the curve Γ(δ_s) is equal to log(aδ_s), where (aδ_s,0) is the point of the horizontal axis given by the limit, as b0, of Γ(δ_s) and therefore 0htop(L)log(2).Figure  illustrates the kneading curves (isentropic curves) for pruned kneading sequence up to dimension 7.

Figure 5. Itinerary and the partition associated to the pruned kneading sequence δs=+11111+11  when (a,b)(1.887,0.0957356).

Figure 5. Itinerary and the partition associated to the pruned kneading sequence δs=+1−1−1−1−1+1−1∗  when (a,b)≈(1.887,0.0957356).

Figure 6. Isentropic curves associated to a pruned kneading sequence.

Figure 6. Isentropic curves associated to a pruned kneading sequence.

Acknowledgments

The author expresses his gratitude to Ricardo Severino and Sandra Vinagre for all their support and encouragement in obtaining these results.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Additional information

Funding

This work was partially supported by projects UIBD/00308/2020 and UIDB/04-674/2020, by FCT-Portuguese Foundation for Science and Technology

References

  • V.M. Alekseev and M.V. Yakobson, Symbolic dynamics and hyperbolic dynamic systems, Phys. Rep.75 (1981), pp. 290–325.
  • Rufus Bowen, Markov partitions for axiom A diffeomorphisms, Amer. J. Math. 92 (1970), pp. 725–747.
  • Diogo Baptista, Ricardo Severino, and Sandra Vinagre, Kneading curves for Lozi maps, Grazer Math. Ber. 354 (2009), pp. 6–14.
  • Y. Ishii, Towards a kneading theory for Lozi mappings I: A solution of pruning front conjecture and the first tangency problem, Nonlinearity 10 (1997), pp. 731–747.
  • R. Lozi, Un attracteur étrange du type attracteur de Hénon, J. Physique (Paris) 39 (1978), pp. 9–10.
  • Z. Liu, H. Xie, Z. Zhu, and Q. Lu, The strange attractor of the Lozi mapping, Int. J. Bifur. Chaos Appl. Sci. Engrg. 2 (1992), pp. 831–839.
  • M. Misiurewicz, Strange attractors for the Lozi mappings, in Nonlinear Dynamics, R.G. Helleman, ed., The New York Academy of Sciences, New York, 1980.
  • Peter Walters, An introduction to ergodic theory, in Graduate Texts in Mathematics, Vol. 79, Springer-Verlag, New York, 1982.