Next Article in Journal
The Phase Space Elementary Cell in Classical and Generalized Statistics
Next Article in Special Issue
Synchronization of a Novel Hyperchaotic Complex-Variable System Based on Finite-Time Stability Theory
Previous Article in Journal
Elastic Properties of CaSiO3 Perovskite from ab initio Molecular Dynamics
Previous Article in Special Issue
Development of Metrics and a Complexity Scale for the Topology of Assembly Supply Chains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relaxation to Fixed Points in the Logistic and Cubic Maps: Analytical and Numerical Investigation

1
Departamento de Física, UNESP, Univ Estadual Paulista Av.24A, 1515, Rio Claro, SP 13506-900, Brazil
2
UNESP, Univ Estadual Paulista, Câmpus São João da Boa Vista, São João da Boa Vista, SP 13874-149, Brazil
3
The Abdus Salam, ICTP, Strada Costiera, 11, Trieste 34151, Italy
*
Author to whom correspondence should be addressed.
Entropy 2013, 15(10), 4310-4318; https://0-doi-org.brum.beds.ac.uk/10.3390/e15104310
Submission received: 14 August 2013 / Revised: 25 September 2013 / Accepted: 1 October 2013 / Published: 14 October 2013
(This article belongs to the Special Issue Dynamical Systems)

Abstract

:
Convergence to a period one fixed point is investigated for both logistic and cubic maps. For the logistic map the relaxation to the fixed point is considered near a transcritical bifurcation while for the cubic map it is near a pitchfork bifurcation. We confirmed that the convergence to the fixed point in both logistic and cubic maps for a region close to the fixed point goes exponentially fast to the fixed point and with a relaxation time described by a power law of exponent - 1 . At the bifurcation point, the exponent is not universal and depends on the type of the bifurcation as well as on the nonlinearity of the map.
Classification:
PACS 05.45.-a; 05.45.Pq; 05.45.Tp

1. Introduction

The motivation for studying discrete mappings comes in the last century from a seminal investigation by May [1], and since then a wide range of applications appeared involving different areas, including biology, physics, chemistry, mathematics, engineering and many others [2,3,4,5,6,7,8,9,10]. Comprehensive discussions on maps can also be found in [11,12,13,14,15,16].
In this paper, we revisit two well known maps namely the logistic map and the cubic map. The nonlinearity of the logistic map is quadratic while for the cubic map, and as the name suggests, is cubic. The two models experience a set of bifurcations reaching the chaos via period-doubling bifurcation and following normal Feingenbaum scaling [17,18]. Before reaching that, the logistic map experiences a transcritical bifurcation where an exchange of stability between fixed points happens and the cubic map experiences a pitchfork bifurcation with a fixed point losing stability, while a twin period one fixed point is born. Our main goal in this paper is to investigate the relaxation to the fixed point around these two bifurcations. We then use a set of numerical simulations and a theoretical investigation to show that, at the bifurcation, a convergence to the fixed point is given by a power law [19] with different exponents for the two bifurcation, while after the bifurcation, the convergence to the fixed point is exponential and with a relaxation time given by a power law with the same exponent for both bifurcations.
The paper is organized as follows. In Section 2 we discuss the mappings considered in this paper and the results obtained. Numerical simulations supporting theoretical findings are given in this Section too. Conclusions are drawn in Section 3.

2. The Mappings and Relaxation to the Fixed Points Investigation

We consider in this section the behavior of the relaxation to the fixed points for two mappings, namely the logistic map given by the expression
x n + 1 = R l x n - R l x n 2
and for the cubic map given by
x n + 1 = R c x n - x n 3
where both R l and R c are control parameters. For our investigations in this paper we consider the ranges R l [ 0 , 4 ] and R c [ 0 , 3 ] . For either R l > 4 and R c > 3 yields the dynamics to go to - and is without interest for us. Figure 1 shows the orbit diagram for the two mappings given by Equations (1) and (2).
As is well known in the literature, the logistic map has two fixed points for R l [ 0 , 3 ] namely
x 1 l = 0
x 2 l = R l - 1 R l
where, according to stability analysis, x 1 l is asymptotically stable for R l [ 0 , 1 ) while x 2 l is asymptotically stable for R l ( 1 , 3 ) for any initial conditions lying x 0 ( 0 , 1 ) . At R l = 1 , the system experiences a transcritical bifurcation and fixed point x 1 l changes stability with x 2 l . For R l = 3 the system exhibits a first period-doubling bifurcation following in a sequence of period-doubling until reaches chaotic behavior. At the bifurcations, the Lyapunov exponents are null given the eigenvalues at the bifurcation points are 1 or - 1 . The sequence of period-doubling follows a Feingenbaum scaling [17,18].
Figure 1. Bifurcation diagrams for: (a) logistic map and; (b) cubic map (for two different initial conditions). The names of some bifurcations are indicated in the figures.
Figure 1. Bifurcation diagrams for: (a) logistic map and; (b) cubic map (for two different initial conditions). The names of some bifurcations are indicated in the figures.
Entropy 15 04310 g001
For the cubic map, there is a fixed point x 1 c = 0 and two period one fixed points x 2 c = R c - 1 and x 3 c = - R c - 1 both are born at R c = 1 . At R c = 1 the system experiences a pitchfork bifurcation when x 1 c loses stability and there is a birth of the x 2 , 3 c fixed points. Each one of them evolve independently suffering period-doubling bifurcations until reach the chaos. Again the sequences of period-doubling are described by the Feingenbaum scaling [17,18]. The two separate chaotic bands are merged only due to a merging chaotic attractors crisis [13,14].
The two regions we are interested in to discuss along this paper correspond to: (i) the region of the transcritical bifurcation in the logistic map and; (ii) the region of the pitchfork bifurcation in the cubic map. Indeed we are seeking to understand and describe how is the relaxation of orbits starting close to the fixed point near both bifurcations. We are then looking to describe the behavior of x approaching to x * denoting the fixed points at R l = 1 and R c = 1 . Exactly like any other variable, x is a function of two entities, i.e., n which is the number of iterations and μ = R l - 1 for the logistic map and μ = R c - 1 for cubic for both R l , c 1 .
Following previous results in the literature [19,20], we start with two hypotheses:
  • For μ = 0 it implies there is an algebraic decay in x so that
    x ( n , μ = 0 ) n β
    where β is a critical exponent and depends on the type of bifurcation.
  • For the parameter μ 0 , we assume the orbit relaxes to the equilibrium exponentially according to
    x ( n , μ ) e - n τ
    where the relaxation time τ has the following form
    τ μ z
    where z is also a critical exponent.
Before showing some theoretical approaches to describe the critical exponents, let us first check what a numerical simulation provides. Figure 2 shows a plot of the convergence to the fixed point considering the logistic map for: (a) μ = 0 and (b) μ 0 .
Figure 2. Convergence to the fixed point for the logistic map considering: (a) μ = 0 where a power law fit furnishes β = - 0.99997 ( 5 ) - 1 and; (b) μ 0 with a slope of z = - 0.994 ( 1 ) - 1 .
Figure 2. Convergence to the fixed point for the logistic map considering: (a) μ = 0 where a power law fit furnishes β = - 0.99997 ( 5 ) - 1 and; (b) μ 0 with a slope of z = - 0.994 ( 1 ) - 1 .
Entropy 15 04310 g002
After doing power law fittings for the two plots of Figure 2 we obtain that β = - 0.99997 ( 5 ) - 1 and z = - 0.994 ( 1 ) - 1 .
On the other hand the convergence to the fixed point for the cubic map is shown in Figure 3 for: (a) μ = 0 leading to a slope of decay given by β = - 0.497 ( 2 ) - 1 / 2 and; (b) μ 0 yielding in a slope of z = - 0.9927 ( 6 ) - 1 .
Figure 3. Convergence to the fixed point for the cubic map considering: (a) μ = 0 where a power law fit furnishes β = - 0.497 ( 2 ) - 1 / 2 and; (b) μ 0 with a slope of z = - 0.9927 ( 6 ) - 1 .
Figure 3. Convergence to the fixed point for the cubic map considering: (a) μ = 0 where a power law fit furnishes β = - 0.497 ( 2 ) - 1 / 2 and; (b) μ 0 with a slope of z = - 0.9927 ( 6 ) - 1 .
Entropy 15 04310 g003
Given the numerical results are now known, we can go ahead with the theoretical argumentation on the characterization of the relaxations. Let us start with the logistic map as example and considering the transcritical bifurcation, i.e., R l = 1 . In this case the mapping is written as
x n + 1 = x n - x n 2
Equation (8) can be rewritten in a more convenient way as
x n + 1 - x n = x n + 1 - x n ( n + 1 ) - n d x d n
that leads to the following approximation
d x d n = - x 2
The approach used in Equation (10) is only valid in the limit of x ( n ) very close to the fixed point. In such a limit, the discrete variables can be treated like a continuous variable, making the derivative possible.
Integrating Equation (10) from both sides leads to
x 0 x 1 x 2 d x = 0 n d n
After doing the integration and rearranging the terms properly we end up with
x ( n ) = 1 n + 1 x 0
As soon as n grows attending to the condition n 1 / x 0 we obtain
x ( n ) n - 1
After a comparison with Equation (5) we see that for the logistic map at R l = 1 the critical exponent β = - 1 , in well agreement to the numerical results presented in Figure 2.
The investigation for R l > 1 is quite similar to the previous case with a minimal detail of subtracting from both sides of Equation (1) a term x n , that leads to
x n + 1 - x n = x n + 1 - x n ( n + 1 ) - n d x d n
yielding at the end with an expression of the type
d x d n = x ( R l - 1 ) - R l x 2
When x is sufficiently close to the fixed point, the second term of Equation (15) which is a quadratic term, becomes rather small as compared to the first one becoming then negligible. Because of this it indeed can be disregarded. Quoting the definition of μ we can rewrite Equation (15) as
d x d n = x μ
which in terms of integral is given as
x 0 x d x x = μ 0 n d n
Doing the integral properly we obtain that
x ( n ) = x 0 e μ n
A comparison with Equation (6) leads us to conclude that the critical exponent z = - 1 , again is in well agreement with the simulation, as confirmed in Figure 2.
Let us now continue the investigation but considering this time the cubic map for R c = 1 . Doing similar procedure as made in Equation (9) we obtain
d x d n = - x 3
After doing the integration and organize the terms properly we obtain that
x ( n ) = 1 2 n + 1 x 0 2 1 2 n - 1 2
for the limit of 2 n 1 / x 0 2 . Comparing the result obtained from Equation (20) with the one presented in Equation (5) we find β = - 1 / 2 which is confirmed by the numerical simulations shown in Figure 3.
Considering the case of R c > 1 but still close to 1, we end up with an expression of the type
d x d n = x μ - x 3
Using similar arguments as used for the logistic map, the cubic term in Equation (21) can be disregarded leading to an identical expression as given by Equation (18). Therefore we conclude that the critical exponent is given by z = - 1 as indeed confirmed by numerical simulations presented in Figure 3.
The two mappings present a bifurcation in R l , c = 1 . At the bifurcation point, both maps exhibit algebraic relaxation to the fixed point but with different critical exponents. For the logistic map β = - 1 while for the cubic map it is given by β = - 1 / 2 . On the other hand, after the bifurcation takes place, the relaxation for the fixed point is given by the same law and with the same critical exponent z = - 1 .
Let us now discuss shortly on the behavior of the entropy for the bifurcations observed at R l , c = 1 . To define a K-entropy, we follow same general discussion as made in [21]. The procedure starts from the evolution in time of a single initial condition converging towards an attractor. The region defining the attractor is therefore covered by a set of discrete cells. An initial condition is started along such cells and its trajectory is followed and marked in the phase space saying what cell is visit at stage n, as for example x ( 0 ) x ( 1 ) x ( 2 ) x ( 3 ) . A second initial condition is started very close to the first one and that may lead to a different sequence of visits. The process continue to a very large number of initial conditions such that an ensemble average on the initial conditions can be made. From the average, a relative number of times a specific sequence of N cells is visited can be defined. Then the entropy S N is given as
S N = < i p ( i ) ln ( p ( i ) ) >
where p ( i ) gives the relative number for the ith sequence and with the summation taken over all possible sequences starting with x ( 0 ) . From this, the K-entropy is then defined as
K = lim N 1 N ( S n - S 0 )
Because the convergence to the attractor at R l , c = 1 is indeed an evolution towards an attracting fixed point, the dynamics is regular. Therefore all sequences starting from the same sufficiently small cell are the same, all orbits follow each other as time passes. This leads to S N = 0 for all the ensemble of N. For a large enough N produces a K = 0 because there is no change in S. This is only observed because the dynamics is regular and no chaos is present for R l , c = 1 .

3. Conclusions

We have considered in this paper the convergence to the fixed point by using two different mappings with different nonlinearities, namely the logistic map with a quadratic nonlinearity and the cubic map, whose nonlinearity is cubic. Both mappings are characterized by a control parameter that induces bifurcations in the system. For the logistic, the convergence to the fixed point at a transcritical bifurcation is given by a power law with exponent β = - 1 . In the cubic map at a pitchfork bifurcation, the decay to the fixed point is also algebraic with slope β = - 1 / 2 . After bifurcation the two systems show an exponential decay to the fixed point whose relaxation time is given by a power law with the same exponent for both maps, namely z = - 1 .

Acknowledgments

JAO sends thanks to CAPES, PROPe/UNESP and CNPq. ERP sends thanks to CAPES. EDL sends thanks to FAPESP (2012/23688-5), CNPq and FUNDUNESP, Brazilian agencies. This research was supported by resources supplied by the Center for Scientific Computing (NCC/GridUNESP) of the São Paulo State University (UNESP).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. May, R.M. Biological populations with non overlapping generations: Stable points, a stable cycles and chaos. Science 1974, 86, 645–647. [Google Scholar] [CrossRef]
  2. Hamacher, K. Dynamical regimes due to technological change in a microeconomical model of production. Chaos 2012, 22, 033149. [Google Scholar] [CrossRef] [PubMed]
  3. McCartney, M. Lyapunov exponents for multi-parameter tent and logistic maps. Chaos 2012, 21, 043104. [Google Scholar] [CrossRef] [PubMed]
  4. Philominathan, P.; Santhiah, M.; Mohamed, I.R.; Murali, K.; Rajasekar, S. Chaotic dynamics of a simple parametrically driven dissipative circuit. Int. J. Bifurc. Chaos 2011, 21, 1927–1933. [Google Scholar] [CrossRef]
  5. Santhiah, M.; Philominathan, P. Statistical dynamics of parametrically perturbed sine-square map. Pramana J. Phys. 2010, 75, 403–414. [Google Scholar] [CrossRef]
  6. Zhang, Y.-G.; Zhang, J.-F.; Ma, Q.; Ma, J.; Wang, Z.-P. Statistical description and forecasting analysis of life system. Int. J. Nonlinear Sci. Numer. Simul. 2010, 11, 157–164. [Google Scholar] [CrossRef]
  7. Hu, W.; Zhao, G.-H.; Zhang, G.; Zhang, J.-Q.; Liu, X.-L. Stabilities and bifurcations of sine dynamic equations on time scale. Acta Phys. Sin. 2012, 17, 170505. [Google Scholar]
  8. Urquizu, M.; Correig, A.M. Fast relaxation transients in a kicked damped oscillator. Chaos, Solitons Fractals 2007, 33, 1292–1306. [Google Scholar] [CrossRef]
  9. Livadiotis, G. Numerical approximation of the percentage of order for one-dimensional maps. Adv. Complex Syst. 2005, 8, 15–32. [Google Scholar] [CrossRef]
  10. Ilhem, D.; Amel, K. One-dimensional and two-dimensional dynamics of cubic maps. Discret. Dyn. Nat. Soc. 2006, 2006, 15840. [Google Scholar] [CrossRef]
  11. Li, T.Y.; Yorke, J.A. Period three implies chaos. Am. Math. Mon. 1975, 82, 985–992. [Google Scholar] [CrossRef]
  12. May, R.M.; Oster, G.A. Bifurcation and dynamical systems in simple ecological models. Am. Nat. 1976, 110, 573–599. [Google Scholar] [CrossRef]
  13. Grebogi, C.; Ott, E.; Yorke, J.A. Chaotic attractors in crisis. Phys. Rev. Lett. 1982, 48, 1507–1510. [Google Scholar] [CrossRef]
  14. Grebogi, C.; Ott, E.; Yorke, J.A. Crises, sudden changes in chaotic attractors, and transient chaos. Physica D 1983, 7, 181–200. [Google Scholar] [CrossRef]
  15. Gallas, J.A.C. Structure of the parameter space of the Hénon map. Phys. Rev. Lett. 1983, 70, 2714–2717. [Google Scholar] [CrossRef] [PubMed]
  16. Collet, P.; Eckmann, J.-P. Iterated Maps on the Interval as Dynamical Systems; Birkhauser: Boston, MA, UA, 1980. [Google Scholar]
  17. Feigenbaum, M.J. Universal metric properties of non-linear transformations. J. of Stat. Phys. 1979, 21, 669–706. [Google Scholar] [CrossRef]
  18. Feigenbaum, M.J. Quantitative universality for a class of non-linear transformations. J. Stat. Phys. 1978, 19, 25–52. [Google Scholar] [CrossRef]
  19. Leonel, E.D.; da Silva, J.K.L.; Kamphorst, S.O. Relaxation and transients in a time-dependent logistic map. Int. J. Bifurc. Chaos 2002, 12, 1667–1674. [Google Scholar] [CrossRef]
  20. Hohenberg, P.C.; Halperin, B.I. Theory of dynamic critical phenomena. Rev. Mod. Phys. 1977, 49, 435–479. [Google Scholar] [CrossRef]
  21. Hilborn, R.C. Chaos and Nonlinear Dynamics: An Introduction for Scientists and Engineers; Oxford University Press: New York, NY, USA, 1994. [Google Scholar]

Share and Cite

MDPI and ACS Style

De Oliveira, J.A.; Papesso, E.R.; Leonel, E.D. Relaxation to Fixed Points in the Logistic and Cubic Maps: Analytical and Numerical Investigation. Entropy 2013, 15, 4310-4318. https://0-doi-org.brum.beds.ac.uk/10.3390/e15104310

AMA Style

De Oliveira JA, Papesso ER, Leonel ED. Relaxation to Fixed Points in the Logistic and Cubic Maps: Analytical and Numerical Investigation. Entropy. 2013; 15(10):4310-4318. https://0-doi-org.brum.beds.ac.uk/10.3390/e15104310

Chicago/Turabian Style

De Oliveira, Juliano A., Edson R. Papesso, and Edson D. Leonel. 2013. "Relaxation to Fixed Points in the Logistic and Cubic Maps: Analytical and Numerical Investigation" Entropy 15, no. 10: 4310-4318. https://0-doi-org.brum.beds.ac.uk/10.3390/e15104310

Article Metrics

Back to TopTop