Physics 303, Fall 2014 1 Logistic Map, Feignbaum number and Lyapunov exponent Physics 303, University of New Mexico M. Gold 1 Logistic map Generating the map: an amazing demo http://en.wikipedia.org/wiki/File:LogisticCobwebChaos.gif! Figure 1: “cobweb” plot
29
Embed
1 Logistic map - Physics & Astronomypanda3.phys.unm.edu/nmcpp/gold/phys303/feign.pdfPhysics 303, Fall 2014 4 Table 1: sequence of bifurcation points for period doubling cascade. (period
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Physics 303, Fall 2014 1
Logistic Map, Feignbaum number and Lyapunov exponentPhysics 303, University of New Mexico
M. Gold
1 Logistic map
Generating the map: an amazing demohttp://en.wikipedia.org/wiki/File:LogisticCobwebChaos.gif!
Table 1: sequence of bifurcation points for period doubling cascade. (period doublingvalues from Chaos: An Introduction to Dynamical Systems, Kathleen T. Alligood) Be-yond value of r=3.5699456720 period doubling ends.
Feigenbaum number is limit of δn in the following sequence. Let rn be the n-thbifurcation point,
∆rn = rn − rn−1 and δn = ∆rn/∆rn+1 Feigenbaum number is universal for period-doubling cascade.
Physics 303, Fall 2014 5
Figure 4: Feigenbaum’s calculation used lots of precision!
Physics 303, Fall 2014 6
Feignbaum, Los Alamos Science, 1980
Figure 5: Feigenbaum’s function g1(x)
Physics 303, Fall 2014 7
2 Quantifying Chaos: Lyapunov exponent
Λ = limn→∞
1
n
n−1∑i=0
ln
∣∣∣∣ dfdx |xi
∣∣∣∣
Lyapunov Exponents for the Logistic Map
zoomhorizontal
vertical
reset
detail0.75 0.80 0.85 0.90 0.95 1.00
0.2
0.4
0.6
0.8
1.0Bifurcation Diagram
0.75 0.80 0.85 0.90 0.95 1.000
Lyapunov Exponent
!
This Demonstration plots the orbit diagram of the logistic map xk+1 = 4 l xk H1 - xkL and the corresponding Lyapunov exponents for different ranges of the parameter l.
The Lyapunov exponent is a parameter characterizing the behavior of a dynamical system. It gives the average rate of exponential divergence from nearby initial conditions. The Lyapunov exponent of the logistic map is given by 1
n ⁄k=1n log @4 l H1 - 2 xkLD.
If the Lyapunov exponent is positive, then the system is chaotic; if it is negative, the system will converge to a periodic state; and if it is zero, there is a bifurcation.
By dragging the locator to the left or right or clicking the plot, you can scroll through the whole range of l-values (0.70–1.0), generate the bifurcation diagram, and plot the Lyapunov exponent over that range.
You can zoom to the position of the locator by using the zoom sliders. To replot the graphs at higher zoom scales, use the "detail" button to increase the number of l values and the number of iterations to 5000.
Figure 6: Lyapunov exponent for logistics map, Demo from Wolfram. “Forcing” valueon horizontal axis is our r/4. A value of λ > 0 implies chaos. For super stable fixedpoints, λ→∞
Figure 13: Identifying cancer with fractal dimensional analysis http://www.
sciencedaily.com/releases/2011/04/110418093852.htm “Fractal analysis of im-ages of breast tissue specimens provides a numeric description of tumor growth patternsas a continuous number between 1 and 2. This number, the fractal dimension, is anobjective and reproducible measure of the complexity of the tissue architecture of thebiopsy specimen. The higher the number, the more abnormal the tissue is.”
Figure 14: By C. Fukushima and J. Westerweel, Technical University of Delft, TheNetherlands. This picture is a false-color image of the far-field of a submerged turbulentjet, made visible by means of laser induced fluorescence
Physics 303, Fall 2014 16
5 Fractals in Art
Figure 15: Mount Fuji Seen Below a Wave at Kanagawa, Hokusai
Physics 303, Fall 2014 17
Figure 16: Autumn Rhythm, Jackson Pollack
Fractal? Specifically, can a true Pollack be authenticated by a fractal analysis?
• Yes, Richard P. Taylor, Adam P. Micolich and David Jonas (1999) https://
• Yes, “Authenticating Pollock paintings using fractal geometry”, RP Taylor, RGuzman, TP Martin, GDR Hall, Pattern Recognition Letters, 2007;
“Perceptual and physiological responses to Jackson Pollock’s fractals”, R Taylor,B Spehar, C Hagerhall, P Van Donkelaar, Frontiers in human neuroscience, 2011;
“Order in Pollock’s chaos”, RP Taylor, Scientific American, 2002
Figure 18: Galileo’s original drawing, showing how larger animals’ bones must begreater in diameter compared to their lengths. d2 ∝ L3
Physics 303, Fall 2014 20
Figure 19: ”Scaling of Skeletal Mass to Body Mass in Birds and Mammals”, HenryD. Prange, John F. Anderson and Hermann Rahn The American Naturalist Vol. 113,No. 1 (Jan., 1979), pp. 103-122
Physics 303, Fall 2014 21
• “Scale and Dinension”, LANL preprint Los Alamos Science Summer/Fall 1984:http://panda3.phys.unm.edu/nmcpp/gold/phys303/west-lanl-preprint.pdf
• “ Chaos and fractals in human physiology,” A. L. Goldberger, D. R. Rigney, B.J. West, Sci Am 1990;262:40-49.
• “ A General Model for the Origin of Allometric Scaling Laws in Biology Au-thor(s): Geoffrey B. West, James H. Brown and Brian J. Enquist ,” Science,New Series, Vol. 276, No. 5309 (Apr. 4, 1997), pp. 122-126.
• “The Fourth Dimension of Life: Fractal Geometry and Allometric Scaling ofOrganisms,” Geoffrey B. West [1,2], James H. Brown [3] and Brian J. Enquist [3]Science, New Series, Vol. 284, No. 5420 (Jun. 4, 1999), pp. 1677-1679.
• “The origin of allometric scaling laws in biology from genomes to ecosystems:towards a quantitative unifying theory of biological structure and organization,”Geoffrey B. West [1,2] and James H. Brown [1,3], The Journal of ExperimentalBiology 208, 1575-1592.
1. The Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, NM 87501, USA.
2. Theoretical Division, MS B285, Los Alamos National Laboratory, Los Alamos,NM 87545, USA.
3. Department of Biology, University of New Mexico, Albuquerque, NM 87131,USA.
Abstract, from “The Fourth Dimension of Life” http://www.sciencemag.org/
content/284/5420/1677.abstract
Fractal-like networks effectively endow life with an additional fourth spatial dimen-sion. This is the origin of quarter-power scaling that is so pervasive in biology. Or-ganisms have evolved hierarchical branching networks that terminate in size-invariantunits, such as capillaries, leaves, mitochondria, and oxidase molecules. Natural selec-tion has tended to maximize both metabolic capacity, by maximizing the scaling ofexchange surface areas, and internal efficiency, by minimizing the scaling of transportdistances and times. These design principles are independent of detailed dynamics andexplicit models and should apply to virtually all organisms.
Quoting again, from “The Fourth Dimension of Life:”
We have proposed that the quarter-power allometric scaling laws and otherfeatures of the dynamical behaviour of biological systems reflect the constraintsinherent in the generic properties of these networks. These were postulated tobe: (i) networks are space-filling in order to service all local biologicallyactive subunits; (ii) the terminal units of the network are invariants; and(iii) performance of the network is maximized by minimizing the energyand other quantities required for resource distribution.
These properties of the ‘average idealized organism’ are presumed to beconsequences of natural selection. Thus, the terminal units of the network whereenergy and resources are exchanged (e.g. leaves, capillaries, cells, mitochondria orchloroplasts), are not reconfigured or rescaled as individuals grow from newborn toadult or as new species evolve. In an analogous fashion, buildings are supplied bybranching networks that terminate in invariant terminal units, such as electricaloutlets or water faucets. The third postulate assumes that the continuous feedbackand fine-tuning implicit in natural selection led to ‘optimized’ systems. For example,of the infinitude of space-filling circulatory systems with invariant terminal units thatcould have evolved, those that have survived the process of natural selection,minimize cardiac output. Such minimization principles are very powerful, becausethey lead to ‘equations of motion’ for network dynamics.
Nota Bene– There is a crucial physics difference between pulsatile and non-pulsatile systems.
Physics 303, Fall 2014 25
Figure 22: “space filling fractal”, From Geoffrey B. West, James H. Brown and BrianJ. Enquist Science, New Series, Vol. 276, No. 5309 (Apr. 4, 1997), pp. 122-126
Figure 23: From “The Fourth Dimension of Life”, Geoffrey B. West, James H. Brownand Brian J. Enquist (1999)
Physics 303, Fall 2014 26
Figure 24: From “The Origin of Allometric Scaling:”
Physics 303, Fall 2014 27
Figure 25: From “The Origin of Allometric Scaling:”
Physics 303, Fall 2014 28
Figure 26: From “The Origin of Allometric Scaling:”
Physics 303, Fall 2014 29
Figure 27: From “The Origin of Allometric Scaling:”