Logistic Map 1 18 2018

Total Page:16

File Type:pdf, Size:1020Kb

Logistic Map 1 18 2018 • CS423/523 website: https://tinyurl.im/0ZCyG – Also available from cs.unm.edu/~melaniem/Courses.html – Subscribe to [email protected] • SFI summer school applications due 1/29 https://www.santafe.edu/engage/learn/schools/sfi- complex-systems-summer-school • Quick Discussion: Mitchell Chapter 1 • CAS in the news • Logistic Map • Project 1, part 1 • Prep with presenters of the Walker paper k is carrying capacity nt= num rabbits at time t xt is fraction of carrying capacity at time t R is intrinsic growth rate: b-d 1 1.0 0.8 0.6 x(t) 0.4 0.2 0 0.0 0 10 20 30 40 50 60 70 80 90 100 t = time 0.0 1.0 Period Doubling Transition to Chaos X-axis: r Y-axis: period of limit cycle • As r increases, between ¼ and ¾: – For any given r, system settles into a limit cycle – Successive period doublings (bifurcations) as r increases – The amount that r increases to get to next period doubling gets smaller and smaller for each new bifurcation. (Feigenbaum’s constant) – At the critical value, the dynamical system falls into essentially an infinite-period limit cycle xt,r ∈ [0,1] x = rx (1 x ) t+1 t − t 1.0 0.9 0.0 0.8 0.65 1.0 0.85 0.90 Feigenbaum’s constant 1 R1 ~ 3.0 (2 = 2 period attractor) 2 R2 ~ 3.44949 (2 = 4 period attractor) 3 R3 ~ 3.54409 (2 = 8 period attractor) R4 ~ 3.564407 . Rinf ~ 3.569946 The rate at which the R values converge is Feigenbaum’s constant ~ 4.6692016 How fast is the next bifurcation relative to the previous bifurcation? From Flake: dk= ak – ak-1 / (ak+1 – ak) where ak is the value of r at which the logistic map bifurcates into a 2k limit cycle Feigenbaum’s constant is the same for all unimodal maps (with a parabolic state space as in Mitchell Fig. 2.4 and Flake 10.2b Characteristics of Chaos • Deterministic • Sensitive • Ergodic • Embedded Characteristics of Chaos • Deterministic: there is structure in the state space. No randomness. Predictable given perfect information. • Sensitive: dependence on initial conditions, perturbations and numerical precision. Butterfly effect. • Ergodic: the state space trajectory will return to all previous local regions. Life repeats itself. • Embedded: unstable limit cycles are embedded in the chaos. • Additional points from Flake: – Stability of fixed points – Pastry: stretch and fold • Short term prediction is possible, information is lost with each recursive application of the function f(x) • The floating point representation needs m + 1 bits to predict values after m recursive iterations of f. • Table 10.1: you can predict when your predictions will fail Logistic Map Exercise xt+1 = R * xt * (1 − xt) X0 = 0.2 or 0.21 R = 2, 3, 3.5 or 3.9 • Do you find: a stable fixed point, periodic cycles or chaos? • Where do you see sensitive dependence on initial conditions? • WIs a population of rabbits growing according to the logistic map a CAS? • Do rabbit populations grow according to the logistic map – What unrealistic simplifications are assumed in the logistic map? – Reductionism is not the only way to simplify Project 1, Part 1 Due Thu Jan 25 on paper in class 1. Plot population vs time for 50 time steps. Demonstrate trajectories with fixed points, periodic cycles and chaotic dynamics. 2. Plot population vs time and demonstrate sensitive dependence on initial conditions. Demonstrate how far into the future you can expect two populations to be correlated for a chaos-generating and a non chaos- generating value of R. 3. Create a bifurcation diagram. You must write code for this yourself (you cannot download it). Label the diagram to show Fiegenbaum’s constant. 4. Create a return map showing xt+1 vs xt for R = 2.75 and for R = 3.85. Create a return map for the population mean of 100 populations with R randomly sampled between 2.5 and 3, and for R between 3.8 and 3.9. 5. Write an introduction, 4 short paragraphs of results and include references following format on the posted template. Provide a caption for each figure. The report is limited to 2 pages. The Project 1 Part 1 report does NOT need to use the template..
Recommended publications
  • An Image Cryptography Using Henon Map and Arnold Cat Map
    International Research Journal of Engineering and Technology (IRJET) e-ISSN: 2395-0056 Volume: 05 Issue: 04 | Apr-2018 www.irjet.net p-ISSN: 2395-0072 An Image Cryptography using Henon Map and Arnold Cat Map. Pranjali Sankhe1, Shruti Pimple2, Surabhi Singh3, Anita Lahane4 1,2,3 UG Student VIII SEM, B.E., Computer Engg., RGIT, Mumbai, India 4Assistant Professor, Department of Computer Engg., RGIT, Mumbai, India ---------------------------------------------------------------------***--------------------------------------------------------------------- Abstract - In this digital world i.e. the transmission of non- 2. METHODOLOGY physical data that has been encoded digitally for the purpose of storage Security is a continuous process via which data can 2.1 HENON MAP be secured from several active and passive attacks. Encryption technique protects the confidentiality of a message or 1. The Henon map is a discrete time dynamic system information which is in the form of multimedia (text, image, introduces by michel henon. and video).In this paper, a new symmetric image encryption 2. The map depends on two parameters, a and b, which algorithm is proposed based on Henon’s chaotic system with for the classical Henon map have values of a = 1.4 and byte sequences applied with a novel approach of pixel shuffling b = 0.3. For the classical values the Henon map is of an image which results in an effective and efficient chaotic. For other values of a and b the map may be encryption of images. The Arnold Cat Map is a discrete system chaotic, intermittent, or converge to a periodic orbit. that stretches and folds its trajectories in phase space. Cryptography is the process of encryption and decryption of 3.
    [Show full text]
  • The Transition Between the Complex Quadratic Map and the Hénon Map
    The Transition Between the Complex Quadratic Map and the Hénon Map by Sarah N. Kabes Technical Report Department of Mathematics and Statistics University of Minnesota Duluth, MN 55812 July 2012 The Transition Between the Complex Quadratic Map and the Hénon map A PROJECT SUBMITTED TO THE FACULTY OF THE GRADUATE SCHOOL OF THE UNIVERSITY OF MINNESOTA BY Sarah Kabes in partial fulfillment of the requirements for the degree of Master of Science July 2012 Acknowledgements Thank you first and foremost to my advisor Dr. Bruce Peckham. Your patience, encouragement, excitement, and support while teaching have made this experience not only possible but enjoyable as well. Additional thanks to my committee members, Dr. Marshall Hampton and Dr. John Pastor for reading and providing suggestions for this project. Furthermore, without the additional assistance from two individuals my project would not be as complete as it is today. Thank you Dr. Harlan Stech for finding the critical value , and thank you Scot Halverson for working with me and the open source code to produce the movie. Of course none of this would be possibly without the continued support of my family and friends. Thank you all for believing in me. i Abstract This paper investigates the transition between two well known dynamical systems in the plane, the complex quadratic map and the Hénon map. Using bifurcation theory, an analysis of the dynamical changes the family of maps undergoes as we follow a “homotopy” from one map to the other is presented. Along with locating common local bifurcations, an additional un-familiar bifurcation at infinity is discussed.
    [Show full text]
  • Chaos: the Mathematics Behind the Butterfly Effect
    Chaos: The Mathematics Behind the Butterfly E↵ect James Manning Advisor: Jan Holly Colby College Mathematics Spring, 2017 1 1. Introduction A butterfly flaps its wings, and a hurricane hits somewhere many miles away. Can these two events possibly be related? This is an adage known to many but understood by few. That fact is based on the difficulty of the mathematics behind the adage. Now, it must be stated that, in fact, the flapping of a butterfly’s wings is not actually known to be the reason for any natural disasters, but the idea of it does get at the driving force of Chaos Theory. The common theme among the two is sensitive dependence on initial conditions. This is an idea that will be revisited later in the paper, because we must first cover the concepts necessary to frame chaos. This paper will explore one, two, and three dimensional systems, maps, bifurcations, limit cycles, attractors, and strange attractors before looking into the mechanics of chaos. Once chaos is introduced, we will look in depth at the Lorenz Equations. 2. One Dimensional Systems We begin our study by looking at nonlinear systems in one dimen- sion. One of the important features of these is the nonlinearity. Non- linearity in an equation evokes behavior that is not easily predicted due to the disproportionate nature of inputs and outputs. Also, the term “system”isoftenamisnomerbecauseitoftenevokestheideaof asystemofequations.Thiswillbethecaseaswemoveourfocuso↵ of one dimension, but for now we do not want to think of a system of equations. In this case, the type of system we want to consider is a first-order system of a single equation.
    [Show full text]
  • Maps, Chaos, and Fractals
    MATH305 Summer Research Project 2006-2007 Maps, Chaos, and Fractals Phillips Williams Department of Mathematics and Statistics University of Canterbury Maps, Chaos, and Fractals Phillipa Williams* MATH305 Mathematics Project University of Canterbury 9 February 2007 Abstract The behaviour and properties of one-dimensional discrete mappings are explored by writing Matlab code to iterate mappings and draw graphs. Fixed points, periodic orbits, and bifurcations are described and chaos is introduced using the logistic map. Symbolic dynamics are used to show that the doubling map and the logistic map have the properties of chaos. The significance of a period-3 orbit is examined and the concept of universality is introduced. Finally the Cantor Set provides a brief example of the use of iterative processes to generate fractals. *supervised by Dr. Alex James, University of Canterbury. 1 Introduction Devaney [1992] describes dynamical systems as "the branch of mathematics that attempts to describe processes in motion)) . Dynamical systems are mathematical models of systems that change with time and can be used to model either discrete or continuous processes. Contin­ uous dynamical systems e.g. mechanical systems, chemical kinetics, or electric circuits can be modeled by differential equations. Discrete dynamical systems are physical systems that involve discrete time intervals, e.g. certain types of population growth, daily fluctuations in the stock market, the spread of cases of infectious diseases, and loans (or deposits) where interest is compounded at fixed intervals. Discrete dynamical systems can be modeled by iterative maps. This project considers one-dimensional discrete dynamical systems. In the first section, the behaviour and properties of one-dimensional maps are examined using both analytical and graphical methods.
    [Show full text]
  • Analysis of the Coupled Logistic Map
    Theoretical Population Biology TP1365 Theoretical Population Biology 54, 1137 (1998) Article No. TP981365 Spatial Structure, Environmental Heterogeneity, and Population Dynamics: Analysis of the Coupled Logistic Map Bruce E. Kendall* Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona 85721, and National Center for Ecological Analysis and Synthesis, University of California, Santa Barbara, California 93106 and Gordon A. Fox- Department of Biology 0116, University of California San Diego, 9500 Gilman Drive, La Jolla, California 92093-0116, and Department of Biology, San Diego State University, San Diego, California 92182 Spatial extent can have two important consequences for population dynamics: It can generate spatial structure, in which individuals interact more intensely with neighbors than with more distant conspecifics, and it allows for environmental heterogeneity, in which habitat quality varies spatially. Studies of these features are difficult to interpret because the models are complex and sometimes idiosyncratic. Here we analyze one of the simplest possible spatial population models, to understand the mathematical basis for the observed patterns: two patches coupled by dispersal, with dynamics in each patch governed by the logistic map. With suitable choices of parameters, this model can represent spatial structure, environmental heterogeneity, or both in combination. We synthesize previous work and new analyses on this model, with two goals: to provide a comprehensive baseline to aid our understanding of more complex spatial models, and to generate predictions about the effects of spatial structure and environmental heterogeneity on population dynamics. Spatial structure alone can generate positive, negative, or zero spatial correlations between patches when dispersal rates are high, medium, or low relative to the complexity of the local dynamics.
    [Show full text]
  • The Logistic Map and Chaos
    The Logistic Map and Chaos Dylana Wilhelm James Madison University April 26, 2016 Dylana Wilhelm Logistic Map April 26, 2016 1 / 18 Logistic Map The logistic map is a first-order difference equation discovered to have complicated dynamics by mathematical biologist Robert May. The general form is given by xn+1 = rxn(1 − xn); where xn is the population of nth generation and r ≥ 0 is the growth rate. Dylana Wilhelm Logistic Map April 26, 2016 2 / 18 Applications of the Logistic Map Generally used in population biology to map the population at any time step to its values at the next time step Additional applications include: Genetics - change in gene frequency Epidemiology - fraction of population infected Economics - relationship between commodity quantity and price Social Sciences - number of people to have heard a rumor Dylana Wilhelm Logistic Map April 26, 2016 3 / 18 Logistic Map Derivation Derived from the logistic difference equation Nn+1 = Nn(r − aNn); by letting x = aN=r Results in simplest non-linear difference equation xn+1 = rxn(1 − xn) Parabola with a maximum value of r=4 at xn = 1=2 For 0 ≤ r ≤ 4, maps 0 ≤ xn ≤ 1 into itself Dylana Wilhelm Logistic Map April 26, 2016 4 / 18 Variation of Growth Rate Limit growth rate to interval 0 ≤ r ≤ 4 Range of behavior as r is varied: Population reaches extinction for r < 1 Non-trivial steady state for 1 < r < 3 Fluctuations in population for r > 3 Dylana Wilhelm Logistic Map April 26, 2016 5 / 18 Variation of Growth Rate Dylana Wilhelm Logistic Map April 26, 2016 6 / 18 Fixed Points of Logistic Map Fixed points satisfy the equation f (x∗) = x∗ = rx∗(1 − x∗); which gives x∗(1 − r + rx∗) = 0: Thus, we have that 1 x∗ = 0 or x∗ = 1 − : r Dylana Wilhelm Logistic Map April 26, 2016 7 / 18 Stability of Fixed Points Stability is given by jf 0(x∗)j < 1.
    [Show full text]
  • Stability and Fractal Patterns of Complex Logistic Map
    BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES • Volume 14, No 3 Sofia • 2014 Print ISSN: 1311-9702; Online ISSN: 1314-4081 DOI: 10.2478/cait-2014-0029 Stability and Fractal Patterns of Complex Logistic Map Bhagwati Prasad, Kuldip Katiyar Department of Mathematics, Jaypee Institute of Information Technology, A-10, Sector-62, Noida, UP-201307 INDIA Emails: [email protected], [email protected] Abstract: The intent of this paper is to study the fractal patterns of one dimensional complex logistic map by finding the optimum values of the control parameter using Ishikawa iterative scheme. The logistic map is shown to have bounded and stable behaviour for larger values of the control parameter. This is well depicted via time series analysis and interesting fractal patterns as well are presented. Keywords: Complex logistic map, Ishikawa iteration, fractals, Julia set. 1. Introduction The name “logistic growth model” is essentially due to Verhulst [20] which he used for the studies on population dynamics. He introduced logistic equation for demographic modelling by extending the Malthus equation of the continuous growth of a population with a view to obtain a stable stationary finite state (see [4]). It took more than hundred years to recognize his founding contributions towards the population dynamics and non-linear sciences. This work received wide attention due to the great implications of the simple looking equation in Chaos theory. In 1963, Edward Lorenz introduced an equivalent version of the logistic model for his famous weather forecast model [13]. R. M a y [15, 16], in 1976, recognized the importance of the logistic model and observed that the continuous time model may not be suitable to reflect the realities in most of the cases and constructed a discrete 14 version of this model.
    [Show full text]
  • Pseudo-Random Number Generator Based on Logistic Chaotic System
    entropy Article Pseudo-Random Number Generator Based on Logistic Chaotic System Luyao Wang and Hai Cheng * Electronic Engineering College, Heilongjiang University, Harbin 150080, China; [email protected] * Correspondence: [email protected] Received: 23 August 2019; Accepted: 27 September 2019; Published: 30 September 2019 Abstract: In recent years, a chaotic system is considered as an important pseudo-random source to pseudo-random number generators (PRNGs). This paper proposes a PRNG based on a modified logistic chaotic system. This chaotic system with fixed system parameters is convergent and its chaotic behavior is analyzed and proved. In order to improve the complexity and randomness of modified PRNGs, the chaotic system parameter denoted by floating point numbers generated by the chaotic system is confused and rearranged to increase its key space and reduce the possibility of an exhaustive attack. It is hard to speculate on the pseudo-random number by chaotic behavior because there is no statistical characteristics and infer the pseudo-random number generated by chaotic behavior. The system parameters of the next chaotic system are related to the chaotic values generated by the previous ones, which makes the PRNG generate enough results. By confusing and rearranging the output sequence, the system parameters of the previous time cannot be gotten from the next time which ensures the security. The analysis shows that the pseudo-random sequence generated by this method has perfect randomness, cryptographic properties and can pass the statistical tests. Keywords: logistic chaotic system; PRNG; Pseudo-random number sequence 1. Introduction With the rapid development of communication technology and the wide use of the Internet and mobile networks, people pay more and more attention to information security.
    [Show full text]
  • Visual Analysis of Nonlinear Dynamical Systems: Chaos, Fractals, Self-Similarity and the Limits of Prediction
    systems Communication Visual Analysis of Nonlinear Dynamical Systems: Chaos, Fractals, Self-Similarity and the Limits of Prediction Geoff Boeing Department of City and Regional Planning, University of California, Berkeley, CA 94720, USA; [email protected]; Tel.: +1-510-642-6000 Academic Editor: Ockie Bosch Received: 7 September 2016; Accepted: 7 November 2016; Published: 13 November 2016 Abstract: Nearly all nontrivial real-world systems are nonlinear dynamical systems. Chaos describes certain nonlinear dynamical systems that have a very sensitive dependence on initial conditions. Chaotic systems are always deterministic and may be very simple, yet they produce completely unpredictable and divergent behavior. Systems of nonlinear equations are difficult to solve analytically, and scientists have relied heavily on visual and qualitative approaches to discover and analyze the dynamics of nonlinearity. Indeed, few fields have drawn as heavily from visualization methods for their seminal innovations: from strange attractors, to bifurcation diagrams, to cobweb plots, to phase diagrams and embedding. Although the social sciences are increasingly studying these types of systems, seminal concepts remain murky or loosely adopted. This article has three aims. First, it argues for several visualization methods to critically analyze and understand the behavior of nonlinear dynamical systems. Second, it uses these visualizations to introduce the foundations of nonlinear dynamics, chaos, fractals, self-similarity and the limits of prediction. Finally, it presents Pynamical, an open-source Python package to easily visualize and explore nonlinear dynamical systems’ behavior. Keywords: visualization; nonlinear dynamics; chaos; fractal; attractor; bifurcation; dynamical systems; prediction; python; logistic map 1. Introduction Chaos theory is a branch of mathematics that deals with nonlinear dynamical systems.
    [Show full text]
  • Coupled Logistic Maps and Non-Linear Differential Equations
    Coupled logistic maps and non-linear differential equations Eytan Katzav1 and Leticia F. Cugliandolo1,2 1Laboratoire de Physique Th´eorique de l’Ecole´ Normale Sup´erieure, 24 rue Lhomond, 75231 Paris Cedex 05, France 2Laboratoire de Physique Th´eorique et Hautes Energies,´ Jussieu, 5`eme ´etage, Tour 25, 4 Place Jussieu, 75252 Paris Cedex 05, France Abstract. We study the continuum space-time limit of a periodic one dimensional array of deterministic logistic maps coupled diffusively. First, we analyse this system in connection with a stochastic one dimensional Kardar-Parisi-Zhang (KPZ) equation for confined surface fluctuations. We compare the large-scale and long-time behaviour of space-time correlations in both systems. The dynamic structure factor of the coupled map lattice (CML) of logistic units in its deep chaotic regime and the usual d = 1 KPZ equation have a similar temporal stretched exponential relaxation. Conversely, the spatial scaling and, in particular, the size dependence are very different due to the intrinsic confinement of the fluctuations in the CML. We discuss the range of values of the non-linear parameter in the logistic map elements and the elastic coefficient coupling neighbours on the ring for which the connection with the KPZ-like equation holds. In the same spirit, we derive a continuum partial differential equation governing the evolution of the Lyapunov vector and we confirm that its space-time behaviour becomes the one of KPZ. Finally, we briefly discuss the interpretation of the continuum limit of the CML as a Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) non-linear diffusion equation with an additional KPZ non-linearity and the possibility of developing travelling wave configurations.
    [Show full text]
  • Symmetry and Self-Similarity with Coupled Logistic Maps W
    Symmetry and Self-Similarity with Coupled Logistic Maps W. Metzler, W. Beau, W. Frees, and A. Ueberla Interdisziplinäre Arbeitsgruppe Mathematisierung, Department of Mathematics, University of Kassel. D-3500 Kassel, F. R. Germany Z. Naturforsch. 42a, 310-318 (1987); received September 24. 1986 The symmetric coupling of two one-dimensional logistic maps displays different types of sym- metric and asymmetric attractors on its route to chaos. Coexisting attractors possess basins in the phase space with Julia-like boundaries. Even classical fractal boundaries occur, when the basins of the coupled map with complex variables are explored. The self-similar features of the objects obtain particular fascination by their symmetries dependent on the map. 1. Introduction 2. The Map Much insight into the dynamics of non-linear It is well known [5] that the iteration scheme systems has been gained through studies of discrete x =x +hx (\-x ), h> 0, (1) mappings. A canonical example of such mappings, k + ] k k k which displays a period-doubling route to chaos, is which is Euler's method for solving the logistic the logistic map [1]. equation The two-dimensional map studied here is a "cross- wise" symmetric coupling of two identical logistic x = x(l- x), maps. Related couplings of nonlinear oscillators may be transformed into have been discussed by Kaneko [2] as well as Hogg & Huberman [3]. They correspondingly found uk + i = ruk( 1 - uk) quasiperiodic behavior with frequency lockings as with well as bifurcations into aperiodic behavior. uk = (h/{l+h))xk, r = 14-/7. (2) The present work should complete their results but sets up priorities on explorations of global Forgetting that (1) is an approximation and treating attraction in the phase space and in the parameter /? as a parameter which need not be small, the space.
    [Show full text]
  • Bug Dynamics: an Example of Chaotic Population Growth and Collapse
    INSTANCES: Incorporating Computational Scientific Thinking Advances into Education & Science Courses Bug Dynamics: An Example of Chaotic Population Growth and Collapse Learning goal: Students will experience some key properties of nonlinear systems, discrete mappings and displaying data, in particular, how to identify bifurcations and chaos. Learning objective Science content objective: • Students will understand how the time dependence of a population is affected by the population’s interactions with its environment, in particular the interdependence of living organisms and their interactions with the components of a biosphere. Science model/computation objectives: • Students will learn how to develop a simple model of a system that can then be computed and used to create simulated bug population data. Characteristics in the simulated data can then be compared to observed phenomena. • Students will develop computational scientific thinking skills by developing a model with successively higher levels of complexity, and by using visualization to uncover simplicity within the complex behaviors found in the model’s output. • Students will learn how mathematical models can be used to describe the relationships among data points. • Students will understand that computer models can be used to explore and develop scientific concepts. Scientific skills objective: • Students will practice the following scientific skills: o Graphing (visualizing) data o Describing trends revealed by data o Proportional and nonproportional (nonlinear) reasoning o Understanding discrete, in contrast to continuous, processes o Test (confirm, modify or reject) hypothesis using analyzed data Activities In this lesson, students will: • Extract data from a table. • Identify key objects/components within a model and study their relationships. • Make statements that describe characteristics of the system.
    [Show full text]