Discrete Dynamics in Nature and Society
Volume 2015 (2015), Article ID 270604, 18 pages
http://dx.doi.org/10.1155/2015/270604

Research Article

Complex Dynamics in Generalized Hénon Map

Meixiang Cai

Institute of Mathematics and Physics, Central South University of Forestry and Technology, Changsha, Hunan 410004, China

Received 19 November 2014; Accepted 10 February 2015

Academic Editor: Viktor Avrutin

Copyright © 2015 Meixiang Cai. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The complex dynamics of generalized Hénon map with nonconstant Jacobian determinant are investigated. The conditions of existence for fold bifurcation, flip bifurcation, and Hopf bifurcation are derived by using center manifold theorem and bifurcation theory and checked up by numerical simulations. Chaos in the sense of Marotto's definition is proved by analytical and numerical methods. The numerical simulations show the consistence with the theoretical analysis and reveal some new complex phenomena which can not be given by theoretical analysis, such as the invariant cycles which are irregular closed graphics, the six and forty-one coexisting invariant cycles, and the two, six, seven, nine, ten, and thirteen coexisting chaotic attractors, and some kinds of strange chaotic attractors.

1. Introduction

The planar map

was first introduced by Hénon [1] as a planar diffeomorphism that imitated essential stretching and folding properties of the Poincaré map of the Lorenz system. This original Hénon map (1) had a strange attractor with fractal structure and had constant Jacobian determinant . Since the late 1970s, Hénon map (1) served as an important but artificial example to illustrate many analytical results and numerical techniques of dynamical theory. For example, Feit [2] introduced the characteristic exponents in order to estimate strange attractors numerically. Marotto [3] proved analytically that the map had a transversal homoclinic orbit, which implied the existence of the chaotic behavior for some parameter values. Curry [4] presented a lot of numerical experiments on Hénon map by using the characteristic exponent, frequency spectrum, and a theorem of Smale [5]. Mora and Viana [6] proposed the more general theory of strange attractors. Sonis [7] presented a detailed description of bifurcation phenomena by using the analysis approach. Cao [8] proved that there exists a set with positive Lebesgue measure, which corresponded to a map possessing a strange attractor. Luo and Guo [9] investigated the complete bifurcation and stability of the stable and unstable periodic solutions and the chaotic layers by introducing the positive and negative iterative mappings. Zhang [10] verified the existence of Wada basin boundaries in a switched Hénon map.

For the measure preserving Hénon map (), there has been some attention. Such as, Brown [11] proved that the measure preserving Hénon map contained an embedded horseshoe for by using geometric methods and a contradiction argument. Kirchgraber and Stoffer [12] proved this Hénon map existing a transversal homoclinic point for a set of parameters which was not small by using shadowing techniques. Jensen [13] proved that the unstable manifold of a hyperbolic fixed point was the iterated limit of a very simple set.

In this paper, we study the following extension of map (1):

where , , and are constants. Following [14], we call this map the generalized Hénon map (GHM), which has nonconstant Jacobian determinant .

Our motivation to study this particular generalized Hénon map is that it appears in the bifurcation analysis of nontransversal homoclinic orbits and heteroclinic cycles and plays an important role in other homoclinic studies. In [15], if a diffeomorphism in has two saddle fixed points connected by two heteroclinic orbits and one of which is nontransversal, then the GHM appears as a rescaled first return map. In [16], If a diffeomorphism in has a codimension 1 homoclinic tangency to a saddle-focus fixed point, the GHM appears when the eigenvalues satisfy some conditions. Moreover, when the so-called effective dimension of the problem in [17] can change, the GHM also can be expected as a rescaled first-return map in other cases of homoclinic and heteroclinic tangencies.

In this paper, the bifurcations and chaos phenomenons in map (2) are investigated. The conditions of existence for fold bifurcation, flip bifurcation, and Hopf bifurcation are derived by using center manifold theorem and bifurcation theory [18]; chaotic behavior in the sense of Marotto’s definition [19] is proved. And numerical simulation results not only show the consistence with the theoretical analysis but also display the new and interesting dynamics behaviors.

This paper is organized as follows. In Section 2, the existence and stability of the fixed points are given. In Section 3, the sufficient conditions of existence for fold bifurcation, flip bifurcation, and Hopf bifurcation are given. In Section 4, we first rigorously prove the existence of chaos in the sense of Marotto’s definition and then give an example to check up the analytic results. The numerical simulations are presented in Section 5 to verify the theoretical analysis and display the new and interesting dynamics. Finally, we give a conclusion in Section 6.

2. Existence and Stability of Fixed Points

The fixed points of map (2) satisfy the following equations:

By a simple analysis about (3), it is easy to obtain the following proposition.

Proposition 1.

(i)If , map (2) has a unique fixed point at ,(ii)if , map (2) has two fixed points at and ,(iii)if , map (2) has three fixed points at , , and , where + +.

The Jacobian matrix of map (2) evaluated at fixed point is given by

The characteristic equation of the Jacobian matrix can be written as

where and .

A simple calculation shows the stability of the fixed points as the following proposition.

Proposition 2. The fixed point of map (2) is stable if one of the following conditions holds:

(i);(ii) and and .

3. Bifurcations

In analysis of the fold bifurcation, flip bifurcation, and Hopf bifurcation of map (2), is used as the bifurcation parameter. By using center manifold theorems and bifurcation theory, the conditions of existences for the bifurcations are given as follows.

Theorem 3. Map (2) undergoes a fold bifurcation at if the following conditions are satisfied at the same time:

Moreover, if (resp., ), the two fixed points bifurcate from for (resp., ) and coalesce as the fixed point at and disappear for (resp., ) (see Appendix A for the proof).

Theorem 4. Map (2) undergoes a flip bifurcation at if the following conditions are satisfied at the same time:

Moreover, if (resp., ), the period-2 points that bifurcate from this fixed point are stable (resp., unstable) (see Appendix B for the proof and expressions of and ).

Theorem 5. Map (2) undergoes a Hopf bifurcation at if the following conditions are satisfied at the same time:

Moreover, if (resp., ) and ; then, an attracting (resp., repelling) invariant closed curve bifurcates from the fixed point for (resp., ) (see Appendix C for the proof and expression of ).

4. Existence of Marotto’s Chaos

In this section, we rigorously prove that map (2) possesses chaotic behavior in the sense of Marotto’s definition in [19].

Suppose that is the fixed point of map (2). We first give the conditions such that the fixed point is a snap-back repeller. The eigenvalues associated with the fixed point are given by

where and .

According to Definition  1 in [19], we begin to find a set in which the norm of conjugate complex eigenvalues exceeds 1 for all . Let , which represents the characteristic equation for any point in . If and , then the characteristic equation have a pair of conjugate complex eigenvalues satisfying .

Let

then,

If , the roots of are , and if , the roots of are .

From the above discussion, we obtain the following results.

Lemma 6. Supposing that , and , we have the following.

(i)If , , , , and or or or , then .(ii)If , , , , and or , then .(iii)If , , , , then .(iv)If , , , , then .(v)If , , , , then .(vi)If , , , , then .(vii)If , , , , then is arbitrary.

Lemma 7. If one of the conditions in Lemma 6 is satisfied and the y-coordinate of fixed point satisfies , then is an expanding fixed point of map (2) in .

Due to Definition 2 of snap-back repeller in [19], we need to find one point such that , , and for some positive integer .

In fact, we have

Now, a map has been constructed to map the point to the fixed point after two iteration if there are solutions different from for (12) and (13). By calculation, the solutions different from for (13) are

Substituting (14) into (12), we have

where .

Next, we expect to find a real root of (15) satisfying . Let , then (15) becomes

Let ; then, (16) becomes

where , , and .

If the real root of (17)  , then (15) has at least a real root .

Obviously, if one of the conditions in Lemma 6 is satisfied and and the solution of (15) satisfies , then is a snap-back repeller in . Thus, we have the following theorem.

Theorem 8. If one of the conditions in Lemma 6 holds and and the solutions and of (14) and (15) satisfy , then the point is a snap-back repeller of map (2), and hence, map (2) is chaotic in the sense of Marotto’s definition, where is given in Lemma 7.

Next, we give specific values of the parameters for illustrating and verifying the conditions in Theorem 8.

Example 9. For , , , , , map (2) has a fixed point , and the eigenvalues of the fixed point are . From Lemma 6, we have , , and , which shows that . We choose ; then, the solution of (15) is and . There exists a point , which is different from and satisfies and . So; is a snap-back repeller.
The bifurcation diagram and the corresponding maximum Lyapunov exponents of map (2) in plane for , , , and are shown in Figures 1(a) and 1(b), respectively. The attractor of map (2) for , , , , and is given in Figure 1(c) (where the maximum Lyapunov exponent ), which is a Maratto’s chaotic attractor.

(a)

(b)

(c)

Figure 1: (a) Bifurcation diagram of map (2) in plane for , , , and . (b) Maximum Lyapunov exponents corresponding to (a). (c) Phase portrait of the chaotic attractor for , , , , and .

5. Numerical Simulations

In this section, we use the bifurcation diagrams, Lyapunov exponents, and phase portraits to illustrate the above analytic results and show the new and more complex dynamical behaviors in the map (2).

5.1. Numerical Simulations for Stability and Codimension One Bifurcations of Fixed Points

In this sbusection, we give the numerical simulations for verifying Propositions 1 and 2 and Theorems 3–5. Four cases are considered as the follows.

Case 1. The bifurcation diagram of fixed points of map (2) in plane for , , , and is given in Figure 2(a). We show that map (2) has one, two, and three fixed points for , , and , respectively. The fold bifurcation (labeled “”) occurs at the fixed point for , and two unstable fixed points bifurcate from this point for , which check the correctness of Theorem 3 (here the coefficient ). The Hopf bifurcations (labeled “” and “”) occur at fixed points and for and , respectively. For example, an attracting invariant cycle bifurcates from the fixed point for , where the coefficients in Theorem 5 are and , which check the correctness of Theorem 5. The stable invariant cycles for and are shown in Figures 2(b) and 2(c), respectively.

(a)

(b)

(c)

Figure 2: (a) Bifurcation diagram of the fixed points of map (2) in plane for , , , and ; (b) and (c): invariant cycles corresponding to Figure 2(a) at and , respectively.

Case 2. The bifurcation diagram of fixed point of map (2) in plane for , , , and is given in Figure 3(a). We show that there is only one fixed point for , and the Hopf bifurcation (labeled “”) occurs at fixed point for , and flip bifurcation (labeled “P-D”) emerges at the fixed point for , which check the correctness of Theorems 5 and 4, respectively. In fact, the fixed point is stable for and loses its stability at flip bifurcation value , an unstable period-2 points bifurcate from the fixed point as increasing (here the coefficients in Theorem 5 are and ).

(a)

(b)

(c)

Figure 3: Bifurcation diagrams of the fixed points of map (2) in plane for (a) , , , and ; (b) , , , and ; (c) , , , and .

Case 3. The bifurcation diagram of fixed point of map (2) in plane for , , , and is given in Figure 3(b). We show that the flip bifurcation (labeled “P-D") occurs at fixed point for , and fold bifurcation (labeled “") emerges from fixed point for .

Case 4. The bifurcation diagram of fixed point of map (2) in plane for , , , and is given in Figure 3(c). We show that the fold bifurcation (labeled “”) emerges from fixed point for , and flip bifurcation (labeled “P-D”) occurs at fixed point for .

5.2. Further Numerical Simulations for Map (2)

In this subsection, we study the complex dynamics changes as the parameters varying with the initial value . The bifurcation parameters are considered in the following five cases.

(a), and , , , .(b), and , , , .(c), and , , , .(d), and , , , .(e), and , , , .

For Case (a). The bifurcation diagram of map (2) in plane and the corresponding maximum Lyapunov exponents are shown in Figures 4(a) and 4(b), respectively. Figures 4(c)–4(e) are the local amplifications of (a). We observe that there is a stable fixed point for ; the fixed point loses its stability as increases; Hopf bifurcation occurs at , and the stable invariant cycle at is shown in Figure 5(a). The region is the invariant cycle region with complex periodic windows, such as period-16, 27 orbits and the invariant cycles at is shown in Figure 5(b). The region is the chaotic region with complex periodic windows, such as period-38 orbit (see Figure 5(c)) and period-11 orbits, and the symmetry breaking of period orbits occurs for the region , and the chaotic attractors at and are shown in the Figures 5(d) and 5(e), respectively. For , the chaotic behavior suddenly converges to period-5 orbits, and the period-5 orbits undergo Hopf bifurcations at separately, and the six coexisting invariant cycles at are shown in Figure 5(f). For the region , we can see the period-doubling bifurcation to chaos and the chaotic behavior suddenly disappears at , and the six coexisting chaotic attractors at , the seven coexisting chaotic attractor at , and the chaotic attractor at are shown in the Figures 5(g)–5(i), respectively.

(a)

(b)

(c)

(d)

(e)

Figure 4: (a) Bifurcation diagram of map (2) in plane for , , , and . (b) Maximum Lyapunov exponents corresponding to (a). (c)–(e). Local amplification of (a) for (c) , (d) , and (e) .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

Figure 5: Phase portraits corresponding to Figure 4: (a) , (b) , (c) , (d) , (e) , (f) , (g) , (h) , and (i) .

For Case (b). The bifurcation diagram of map (2) in plane and the corresponding maximum Lyapunov exponents are shown in Figures 6(a) and 6(b), respectively. In Figure 6, we observe period-doubling bifurcation to chaos and chaotic region with period-6 windows. The two coexisting chaotic attractors at and the chaotic attractor at are shown in Figures 7(a) and 7(c), respectively, and the period-8 orbits at is shown in Figure 7(b).

(a)

(b)

Figure 6: (a) Bifurcation diagram of map (2) in plane for , , , and . (b) Maximum Lyapunov exponents corresponding to (a).

(a)

(b)

(c)

Figure 7: Phase portraits corresponding to Figure 6: (a) , (b) , and (c) .

For Case (c). The bifurcation diagram of map (2) in plane and the corresponding maximum Lyapunov exponents are shown in Figures 8(a) and 8(b), respectively, which show the chaotic region with complex periodic windows. The chaotic attractors at , , and are shown in Figures 9(a)–9(c), respectively.

(a)

(b)

Figure 8: (a) Bifurcation diagram of map (2) in plane for , , , and . (b) Maximum Lyapunov exponents corresponding to (a).

(a)

(b)

(c)

Figure 9: Phase portraits corresponding to Figure 8: (a) , (b) , and (c) .

For Case (d). The bifurcation diagram of map (2) in plane and the corresponding maximum Lyapunov exponents are shown in Figures 10(a) and 10(b), respectively, and the local amplification of (a) for is shown in Figure 10(c). we observe that period-doubling bifurcation to chaos and inverse period-doubling bifurcation to chaos, and the Hopf bifurcation occurs at . The invariant cycle at , the two coexisting chaotic attractors at and the chaotic attractor at are shown in Figures 11(a) and 11(c), respectively.

(a)

(b)

(c)

Figure 10: (a) Bifurcation diagram of map (2) in plane for , , , and . (b) Maximum Lyapunov exponents corresponding to (a). (c) Local amplification of (a).

(a)

(b)

(c)

Figure 11: Phase portraits corresponding to Figure 10: (a) , (b) , and (c) .

For Case (e). The bifurcation diagram of map (2) in plane and the maximum Lyapunov exponents corresponding to (a) are shown in Figures 12(a) and 12(b), respectively. Figures 12(c)–12(h) are the local amplifications of (a). In Figure 12, we observe that the jumping behavior of period-3 orbits to period-1 orbit at , Hopf bifurcation occurs at , and the invariant cycle at is shown in Figure 13(a), and the invariant cycle suddenly converges to period-4 orbits at . The region is the invariant cycle region with complex periodic windows, such as period-21, 17 orbits. The invariant cycle at and the forty-one coexisting invariant cycles at are shown in Figures 13(b) and 13(c), respectively. For , the periodic orbits suddenly become chaotic behaviors. Moreover, the chaotic behaviors suddenly change into a special chaotic behaviors which has 13 small size chaotic attractors as increases. The chaotic attractor at and the thirteen coexisting chaotic attractors at are shown in Figures 13(d) and 13(e), respectively. For , the chaotic behaviors suddenly converge period-13 orbits, which quickly undergo Hopf bifurcation and change into period-13 orbits again. For the region , the period-doubling bifurcation to chaos occurs several times and the chaotic behaviors with complex periodic windows, and the nine coexisting chaotic attractors at , the chaotic attractor at , the ten coexisting chaotic attractors at , and the chaotic attractor at are shown in Figures 13(f)–13(i), respectively.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

Figure 12: (a) Bifurcation diagram of map (2) in plane for , , , and . (b) Maximum Lyapunov exponents corresponding to (a). (c)–(h) Local amplification of (a) for (c) , (d) , (e) , (f) , (g) , and (h) .

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

Figure 13: Phase portraits corresponding to Figure 12: (a) , (b) , (c) , (d) , (e) , (f) , (g) , (h) , and (i) .

6. Conclusion

In this paper, the generalized Hénon map with nonconstant Jacobian determinant is investigated. By using the analytic and numerical methods, the existence conditions of fold bifurcation, flip bifurcation, Hopf bifurcation, and Marotto’s chaos are obtained. Moreover, we discover some interesting dynamical behaviors, such as the invariant cycles which are irregular closed graphics, and the six and forty-one coexisting invariant cycles, and the two, six, seven, nine, ten, and thirteen coexisting chaotic attractors, and some kinds of strange chaotic attractors. The numerical simulations show the consistence with the theoretical analysis and exhibit some specific properties which can not be obtained by theoretical analysis.

Appendices

A. For the Proof of Theorem 3

The characteristic equation associated with the linearization of map (2) about the fixed point is given by

where , .

Supposing that , the eigenvalues of the fixed point are and . For codimension-one bifurcation, it requires , which leads to and .

Let , let , and let , we transform the fixed point of map (2) to the origin and consider the parameter as a new dependent variable, we can write the map (2) in the following form:

where .

We construct an invertible matrix

and use the translation ; then, the map (A.2) becomes

where + + , , and .

By center manifold theory, the stability of near can be determined by studying a one-parameter family of equation on a center manifold, which can be represented as follows:

for and sufficiently small.

Assume that a center manifold is of the form , the center manifold must satisfy

Since only the terms with orders lower than 3 for map (A.4) are concerned, the map restricted to the center manifold directly from (A.4) can be given by

where

Since , , , and , the fixed point is a fold bifurcation point for map (A.4) and the number of fixed points is changed at as .

B. For the Proof of Theorem 4

Supposing that , the eigenvalues of the fixed point are and .

The condition leads to and , and the condition leads to .

Let , let , and let ; we transform the fixed point of map (2) to the origin and consider as a new dependent variable; then, map (2) becomes

where .

We construct an invertible matrix

and use the translation ; map (B.1) becomes

where + + + + ,  , , and .

We again apply the center manifold theorem to determine the nature of the bifurcation of the fixed point at . There exists a center manifold for map (B.3), which can be represented as follows:

for and sufficiently small.

Assuming a center manifold of the form , the center manifold must satisfy

By the approximate calculation for center manifold, the coefficients + , + + , , and the map restricted to the center manifold is given by

where

In order for map (B.6) to undergo a flip bifurcation, it is also required that

C. For the Proof of Theorem 5

The eigenvalues of the characteristic equation at fixed point are , where + + . The eigenvalues are complex conjugate for , which leads to

Assuming that

we have .

We translate the fixed point of map (2) to the origin by using the translation , ; map (2) becomes

For , the eigenvalues of the matrix associated with the linearized map (2) at fixed point are complex conjugate with modulus 1 which are written as

where and .

Under the conditions (C.1)–(C.3), there are and . In addition, if , which leads to and , then we have , .

Let and use the translation ; map (C.4) becomes

where . Let and let ; then, +.

Notice that map (C.6) is exactly in the form on the center manifold, in which the coefficient (see [18]) is given by

where , , , and .

Thus, some complicated calculations give

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The author would like to thank the reviewers for their helpful comments and suggestions. This work was supported by Hunan Provincial Natural Science Foundation of China (no. 13JJ4088 and no. 14JJ3114), the General Project of Hunan Provincial Education Department (no. 14C1191), and the National Natural Science Foundation of China (no. 11301551, no. 11426221, and no. 11426223).

References

M. Hénon, “A two-dimensional mapping with a strange attractor,” Communications in Mathematical Physics, vol. 50, no. 1, pp. 69–77, 1976.

S. D. Feit, “Characteristic exponents and strange attractors,” Communications in Mathematical Physics, vol. 61, no. 3, pp. 249–260, 1978.

F. R. Marotto, “Chaotic behavior in the Hénon mapping,” Communications in Mathematical Physics, vol. 68, no. 2, pp. 187–194, 1979.

J. H. Curry, “On the Hénon transformation,” Communications in Mathematical Physics, vol. 68, no. 2, pp. 129–140, 1979.

S. Smale, “Differentiable dynamical systems,” Bulletin of the American Mathematical Society, vol. 73, pp. 747–817, 1967.

L. Mora and M. Viana, “Abundance of strange attractors,” Acta Mathematica, vol. 171, no. 1, pp. 1–71, 1993.

M. Sonis, “Once more on Hénon map: analysis of bifurcations,” Chaos, Solitons and Fractals, vol. 7, no. 12, pp. 2215–2234, 1996.

Y. Cao, “The nonwandering set of some Hénon map,” Chinese Science Bulletin, vol. 44, no. 7, pp. 590–594, 1999.

A. C. J. Luo and Y. Guo, “Complete bifurcation behaviors of a henon map,” Dynamical Systems: Discontinuity, Stochasticity and Time-Delay, pp. 37–47, 2010.

Y. Zhang, “Switching-induced Wada basin boundaries in the Hénon map,” Nonlinear Dynamics, vol. 73, no. 4, pp. 2221–2229, 2013.

R. Brown, “Horseshoes in the measure-preserving Hénon map,” Ergodic Theory and Dynamical Systems, vol. 15, no. 6, pp. 1045–1059, 1995.

U. Kirchgraber and D. Stoffer, “Transversal homoclinic points of the Hénon map,” Annali di Matematica Pura ed Applicata, vol. 185, no. 5, pp. S187–S204, 2006.

E. Jensen, “A new construction of the unstable manifold for the measure-preserving Hénon map,” Proceedings of the American Mathematical Society, vol. 136, no. 1, pp. 181–192, 2008.

S. V. Gonchenko and V. S. Gonchenko, On Andronov-Hopf Bifurcations of Two-Dimensional Diffeomorphisms with Homoclinic Tangencies, WIAS, Berlin , Germany, 2000.

S. V. Gonchenko, L. P. Shilnikov, and O. V. Stenkin, “On newhouse regions with infinitely many stable and unstable tori,” in Proceedings of the International Conference on Progress in Nonlinear Science, vol. 1, pp. 80–102, RAS Institute of Applied Physics, Nizhny Nogorod, Russia, 2002.

S. V. Gonchenko, D. V. Turaev, and L. P. Shil'nikov, “On dynamic properties of diffeomorphisms with homoclinic tangency,” Journal of Mathematical Sciences, vol. 126, no. 4, pp. 1317–1343, 2005.

D. V. Turaev, “On dimension of non-local bifurcational problems,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 6, no. 5, pp. 919–948, 1996.

J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, Springer, New York, NY, USA, 1983.

F. R. Marotto, “Snap-back repellers imply chaos in ,” Journal of Mathematical Analysis and Applications, vol. 63, no. 1, pp. 199–223, 1978.