Chapter 5: Fractals II Goals: • to Understand Why the Scaling

Total Page:16

File Type:pdf, Size:1020Kb

Chapter 5: Fractals II Goals: • to Understand Why the Scaling Chapter 5: Fractals II In this chapter we address the question of why and have the values that they do and why they are the same for many dierent mapping functions. We will write an equation that embodies the observation of dilation symmetry, and Goals: we will nd that this equation determines the values of the exponents and . To understand why the scaling exhibited by period-doubling bi- In other words, knowing that dilation symmetry exists is enough information furcation sequences of one-dimensional maps is universal, and to to determine what the magnication factors must be. This result explains the calculate the scaling exponents. observation of universality (dierent maps having the same exponents). The To examine other natural processes which give rise to fractal book by Cvitanovic in the reference list is a good resource for further reading structures. about this universality. Several arguments in this chapter are based on this source. This chapter continues our study of dilation-symmetric structures, or frac- We rst review how the iteration of a smooth function leads to n-cycles of tals. These structures are sometimes called “self-similar”. First, we will pursue ever longer length. Any smooth, positive function f(x)=rfo(x) that vanishes our investigation of the period-doubling bifurcation sequence of the logistic map for x = 0 and x = 1 must cross the line f = x when the “amplitude” r exceeds and get some understanding of why the scaling exponents and which char- some positive value. When r is just at this threshold value, f 0 = 1, and as r acterize it are universal. We hope that this will give you some sense of how increases from there, the derivitive at one of the crossing points becomes smaller understanding fractals can give new insight into natural phenomena. Then we than 1, so that the xed point is stable. We follow this point as r continues to will survey a few other systems in nature in which fractals arise. Some of these increases, and see that f 0 becomes more and more negative, eventually reaching other systems are not nearly so well understood, so perhaps you will get some f 0 = 1. At this point the xed point is marginally stable. Any further increase idea of the scope of the question of why fractals are seen in many dierent in r will drive the xed point unstable. To see what happens then, we look at situations. the second iterate f (2)(x). Up to this value of r it has had a stable xed point at the same x value as for f(x). When f 0 becomes -1 at the xed point, f (2)0 A. The universality of the period-doubling bifurcation sequence. 0 0 becomes +1, i.e., f f . It becomes tangent to the f = x line, making the f (2) In the previous chapter and in Required Project 3 we investigate the scaling xed point marginally stable. As r increases from here, this marginally stable behavior of the period-doubling sequence of the logistic map. We saw that the point moves and divides into two, one stable and one unstable. This is the same k sequence of r-values for the orbits of period 2 obeys behavior that f did for smaller r. It must be the same since f (2) is just another smooth function, as f was. The stable xed point represents a two-cycle of f,so r r lim k 1 k 2 = , (2) →∞ that there must be a second stable xed point of f nearby. As r continues to k rk rk1 (2)0 increase, the derivitive f becomes -1 and this xed point becomes unstable, (4) (2) (2) 2k 1 by the same process that happened with f. But at this point f = f (f ) where =4.7...., and that the sequence of values yk f (x) satises (4) (2) 1/2 is just tangent to the f = x line, and f repeats the process that f and f did before it. yk1 yk2 lim , = , Each time the stable xed points switch to a higher iterate of f, it means k→∞ y y k k 1 that f has an n-cycle or period of twice the length as before. For each pe- with =2.5..... These results mean that the bifurcation diagram looks the riod doubling there is a corresponding r value, that we denote rk for the kth (2k) same when it is magnied about the point (x =1/2, r = r∞ =3.57....)bya doubling. Thus between rk and rk+1 f has a stable xed point that aren’t factor in the x-direction and a factor in the r-direction. This mapping stable for lower iterates. By continuing to increase r it appears that we may of the bifurcation diagram onto itself under certain rescaling factors is known increase the length of these cycles as much as we wish, leading to behavior of as dilation symmetry (because the diagram is invariant under dilation of co- arbitrarily great complexity. For the logistic map and many other f functions ordinate scales). In Required Project 3 you show that and are not noticeably the limit cycles become innite for some nite r∞. This complexity as r → r∞ dierent for some dierent choices of the map function. is what we want to probe. Physics 251/CS 279/Math 292 Chapter 5 typeset February 11, 2001 page 1 Physics 251/CS 279/Math 292 Chapter 5 typeset February 11, 2001 page 2 (2k) When we are close to r∞, an arbitrarily slight change in the original func- This means that the function gk dened by a(k)(f (x) x0)=gk(zk) should k tion f can evidently cause a qualitative change in behavior: a period doubling. be smooth, even though f (2 ) itself is not smooth. Figure 5.1 illustrates the We saw above that each period doubling can be understood by looking at a idea using f (2) and f (4). tiny range of x. We expect the same to be true as we approach r∞. A clever observation shows us what x to focus on. We saw that the derivitive of our (2k) iterate f starts out at the marginally stable value of 1 at r = rk. Then it steadily decreases to -1 as we approach r0. At some r value in between, the derivitive must vanish. We denote this r asr ˜k. When the derivitive vanishes at a xed point, this point is called superstable, as noted above. We may deter- mine a superstable xed point x of any iterate by using the Floquet-multiplier (2k) (2k)0 0 0 0 expression for the derivitive of f : f = f (x1)f (x2)...f (x2k ). At a su- perstable point, this multiplier must vanish. This means that one of the f 0 s in the product must vanish. One point of the 2k cycle must be at a maximum or minimum of f(x). These maxima and minima don’t shift as we change r; r is just a multiplicative factor. For example in the logistic map, there is a single maximum located at x =1/2. From now on, we’ll focus on some chosen maximum or minimum of f located at x0 and consider the region around x0. As k →∞the some xed points at the period-doubling rk must lie close to the superstable xed point at the correspondingr ˜k at x0. Limiting dependence on x To see how the instabilities pile up as r → r∞ we k clearly need to know the functional form of f (2 )(x) near x = x0. This functional form cannot be very smooth. We know that each k point in the current 2k-cycle must be a xed point of f (2 ). Thus this function must cross the f = x line 2k times. The number of wiggles must double each time k increases by one. We see this happening as we iterate the logistic map: k each time k increases by one f (2 ) becomes a polynomial of twice as high an order. k We cannot hope to have a manageable limit for f (2 ) as k →∞unless we use a scale of x which expands to compensate for the increasing wiggliness Figure 5.1. Dark line: f (2) light line, f (4) near superstable points. A k k of f (2 ). Accordingly we adopt a rescaled co-ordinate for describing f (2 )(x) section of the f (2) is copied below it. This section is reduced by a (4) called zk,and dened by zk = a(k)(x x0). Evidently a(k) must become huge, factor of 2.74 and inverted and placed near the f curve to show k so that the tiny region of x where f (2 ) looks smooth gets expanded to a nite, that the f (4) curve resembles the f (2) curve except for a change of nonvanishing interval of z. Now when we increase k we can hope to account for horizontal and vertical scale. (2k) the increased wiggliness of f (x) near x0 by the scale factor a(k), so that the (2k) remaining k dependence is smooth. To achieve the hoped-for smooth function, We may use the known behavior of f to infer the behavior of the hy- k k (2 ) we must expand the scale of f (2 ) as well as the scale of x. In words, the pothesized function g. What we know about f is the way it is related to the k+1 behavior we seek is next function in the series,f (2 ): k k+1 k k expanded f (2 )(x)=[smooth function of](expanded x) f (2 )(x)=f (2 )(f (2 )(x)).
Recommended publications
  • Math Morphing Proximate and Evolutionary Mechanisms
    Curriculum Units by Fellows of the Yale-New Haven Teachers Institute 2009 Volume V: Evolutionary Medicine Math Morphing Proximate and Evolutionary Mechanisms Curriculum Unit 09.05.09 by Kenneth William Spinka Introduction Background Essential Questions Lesson Plans Website Student Resources Glossary Of Terms Bibliography Appendix Introduction An important theoretical development was Nikolaas Tinbergen's distinction made originally in ethology between evolutionary and proximate mechanisms; Randolph M. Nesse and George C. Williams summarize its relevance to medicine: All biological traits need two kinds of explanation: proximate and evolutionary. The proximate explanation for a disease describes what is wrong in the bodily mechanism of individuals affected Curriculum Unit 09.05.09 1 of 27 by it. An evolutionary explanation is completely different. Instead of explaining why people are different, it explains why we are all the same in ways that leave us vulnerable to disease. Why do we all have wisdom teeth, an appendix, and cells that if triggered can rampantly multiply out of control? [1] A fractal is generally "a rough or fragmented geometric shape that can be split into parts, each of which is (at least approximately) a reduced-size copy of the whole," a property called self-similarity. The term was coined by Beno?t Mandelbrot in 1975 and was derived from the Latin fractus meaning "broken" or "fractured." A mathematical fractal is based on an equation that undergoes iteration, a form of feedback based on recursion. http://www.kwsi.com/ynhti2009/image01.html A fractal often has the following features: 1. It has a fine structure at arbitrarily small scales.
    [Show full text]
  • Dm/Dt)Vs 2 Fρ = (Duc/Dρ)N = (Mvdv/Dρ)N + (1/2)(Dm/Dρ)V N Y Análogamente La Fuerza Total Sobre M Ahora Será
    0 JUAN RIUS-CAMPS LA IRREVERSIBILIDAD Y EL CAOS 8 DE DICIEMBRE DE 2009 EDICIONES ORDIS 1 2 EDICIONES ORDIS GRAN VÍA DE CARLOS III 59 . 2º . 4ª 08028 BARCELONA JUAN RIUS – CAMPS 8 de Diciembre de 2009 3 4 LA IRREVERSIBILIDAD Y EL CAOS ÍNDICE CAPÍTULO I p. 9 LOS FUNDAMENTOS COSMOLÓGICOS DE LA MECÁNICA Y LAS LEYES FUNDAMENTALES DE LA DINÁMICA INTRODUCCIÓN p. 9 MATERIA Y FORMA p. 10 CAPÍTULO II p. 33 DINÁMICA DE SISTEMAS MACÁNICOS IRREVERSIBLES INTRODUCCIÓN p. 33 SISTEMA DE PUNTOS MATERIALES p. 34 A. FUNDAMENTOS p. 34 B. ESTUDIO ANALÍTICO DE LA EXPRESIÓN DE LA FUERZA EN LA NUEVA DINÁMICA (ND) p. 37 C ENTROPÍA: SEGUNDA LEY FUNDAMENTAL p. 44 D SENTIDO CINEMÁTICO DE LA VELOCIDAD ANGULAR ω* p. 50 5 CAPÍTULO III p. 59 ECUACIONES DE ONDA Y “ECUACIONES DE MAXWELL” A. DETERMINACIÓN DE LA ECUACIÓN DE ONDA CUANDO Up = Up(P, t) p. 59 B. DEDUCCIÓN DE LAS “ECUACIONES DE MAXWELL” p. 67 CAPÍTULO IV p. 74 NOTAS COMPLEMENTARIAS. A. CONDICIÓN DE FUERZA CENTRAL (ND) p. 74 B. CONDICIÓN DE FUERZA CENTRAL (Caso singular de espiral logarítmica) p. 75 C. FUERZA CENTRAL (Expresada en coordenadas polares) p.78 81 CAPÍTULO V p. 81 IRREVERSIBILIDAD Y CAOS CAPÍTULO VI p. 110 PRUEBAS EXPERIMENTALES 6 7 8 CAPÍTULO I LOS FUNDAMENTOS COSMOLÓGICOS DE LA MECÁNICA Y LAS LEYES FUNDAMENTALES DE LA DINÁMICA INTRODUCCIÓN Para expresar nuestro punto de partida se hace preciso afirmar que éste no ha sido mayormente físico sino filosófico, en el sentido de la Metafísica de la Naturaleza o Cosmología, es decir, en la vía abierta desde hace siglos por ARISTÓTELES, continuada por TOMÁS de AQUINO y otros pensadores.
    [Show full text]
  • Birth of a New Class of Period-Doubling Scaling Behavior As a Result of Bifurcation in the Renormalization Equation
    J Stat Phys (2008) 130: 599–616 DOI 10.1007/s10955-007-9442-6 Birth of a New Class of Period-Doubling Scaling Behavior as a Result of Bifurcation in the Renormalization Equation S.P. Kuznetsov · A.A. Mailybaev · I.R. Sataev Received: 2 December 2006 / Accepted: 19 September 2007 / Published online: 24 October 2007 © Springer Science+Business Media, LLC 2007 Abstract It is found that a fixed point of the renormalization group equation corresponding to a system of a unimodal map with extremum of power κ and a map summarizing values of a function of the dynamical variable of the first subsystem, undergoes a bifurcation in the course of increase of κ. It occurs at κc = 1.984396 and results in a birth of the period-2 stationary solution of the RG equation. At κ = 2 this period-2 solution corresponds to the universal period-doubling behavior discovered earlier and denoted as the C-type criticality (Kuznetsov and Sataev in Phys. Lett. A 162:236–242, 1992). By combination of analytical methods and numerical computations we obtain and analyze an asymptotic expansion of the period-2 solution in powers of κ = κ − κc. The developed approach resembles the -expansion in the phase transition theory, in which a “trivial” stationary point of the RG transformation undergoes a bifurcation that gives rise to a new fixed point responsible for the critical behavior with nontrivial critical indices. 1 Introduction One of challenges in the modern nonlinear dynamics concerns elaboration of the renormal- ization group (RG) methods for studying transitions to chaos [1, 2].
    [Show full text]
  • Universality in Transitions to Chaos
    . Chaos: Classical and Quantum I: Deterministic Chaos Predrag Cvitanovic´ – Roberto Artuso – Ronnie Mainieri – Gregor Tanner – Gabor´ Vattay —————————————————————- ChaosBook.org version15.8, Oct 18 2016 printed November 10, 2016 ChaosBook.org comments to: [email protected] Chapter 31 Universality in transitions to chaos When you come to a fork in the road, take it! —Yogi Berra he developments that we shall describe next are one of those pleasing de- monstrations of the unity of physics. The key discovery was made by a T physicist not trained to work on problems of turbulence. In the fall of 1975 Mitchell J. Feigenbaum, an elementary particle theorist, discovered a universal transition to chaos in 1-dimensional unimodal map dynamics. At the time the physical implications of the discovery were nil. During the next few years, howe- ver, numerical and mathematical studies established this universality in a number of realistic models in various physical settings, and in 1980 the universality theory passed its first experimental test. The discovery was that large classes of nonlinear systems exhibit transitions to chaos which are universal and quantitatively measurable. This advance was akin to (and inspired by) earlier advances in the theory of phase transitions; for the first time one could, predict and measure ‘critical exponents’ for turbulence. But the breakthrough consisted not so much in discovering a new set of universal numbers, as in developing a new way to solve strongly nonlinear physical pro- blems. Traditionally, we use regular motions (harmonic oscillators, plane waves, free particles, etc.) as zeroth-order approximations to physical systems, and ac- count for weak nonlinearities perturbational.
    [Show full text]
  • A Complete Bibliography of the Journal of Statistical Physics: 1980–1989
    A Complete Bibliography of the Journal of Statistical Physics: 1980{1989 Nelson H. F. Beebe University of Utah Department of Mathematics, 110 LCB 155 S 1400 E RM 233 Salt Lake City, UT 84112-0090 USA Tel: +1 801 581 5254 FAX: +1 801 581 4148 E-mail: [email protected], [email protected], [email protected] (Internet) WWW URL: http://www.math.utah.edu/~beebe/ 01 March 2019 Version 1.05 Title word cross-reference (1 + 1) [Cor87]. (1=2; 1=2) [Mar87]. (2 + 1) [Cor88]. (d + 1) [Sch86a]. 1 − µjjxjjz [de 88]. 1=2 [TS81, KLS88a]. 1=f [MS83b, NT81]. 1=rn [KH82]. 1=jx − yj2 [ACCN88]. 100 ≤ Γ ≤ 160 [HMV84]. 2n [Kai82]. 3 [BDS86, DC89, Zha89a]. 4 3 4 [Lin87]. 8 [Lin87]. [HP85, deV85]. [HP85, Tou89]. 2 [BR80, BWdF88]. 3 [BWdF88, CR86]. z [BWdF88]. A [Rue86a, Vai89]. AgRBf−g [Spo83]. An (1) 2 [KY88b]. An [KY88b]. B6 [Kra82b]. B7 [Kra82b]. β [GGP86]. β < 8π [LSM89]. β ≥ 2/δ [New87]. d [BL85, O'C84a, OB84, dSVT87, YHM80]. d =2 (1) [GS85]. d = 2 + 1 [BP89]. d =2+ [SH82]. d = 3 [ST81b]. Dn [KY88b]. Dn 4 [KY88b]. δ [OP89]. dv =2∆− γ [ST81b]. g0: φ :d [BK81]. γ [MIH86, SD86]. Γ = 2 [CFS83, FJS83]. H [AC84, BS89b, Del84b, Mar81, Pia87]. i(π/2)kT [Wu86]. K [GGM86, MT84]. L1 [Cer88, Pol89b]. L1 [Ark83]. λ(rφ)4 [GKT84]. ≤ 2 [GM88b, GM96]. M [de 88, KST89]. Zk−1 [May85]. µ [EC81]. N [RG86a, AZ85, Bau83, Bur86, FTZ84, HR83a, Kal81, KD82a, KPS89, KST89, Kur83b, Lav82b, Sch87c, WS82]. n−2 [Sim81a]. O(n) 1 2 4 4 [ABCT86, Dup87, Lar86].
    [Show full text]
  • Universality in Transitions to Chaos
    Chapter 26 Universality in transitions to chaos The developments that we shall describe next are one of those pleasing demonstrations of the unity of physics. The key discovery was made by a physicist not trained to work on problems of turbulence. In the fall of 1975 Mitchell J. Feigenbaum, an elementary particle theorist, discovered a universal transition to chaos in one-dimensional unimodal map dynamics. At the time the physical implications of the discovery were nil. During the next few years, however, numerical and mathematical studies established this universality in a number of realistic models in various physical settings, and in 1980 the universality theory passed its first experimental test. The discovery was that large classes of nonlinear systems exhibit transi- tions to chaos which are universal and quantitatively measurable. This ad- vance was akin to (and inspired by) earlier advances in the theory of phase transitions; for the first time one could, predict and measure “critical ex- ponents” for turbulence. But the breakthrough consisted not so much in discovering a new set of universal numbers, as in developing a new way to solve strongly nonlinear physical problems. Traditionally, we use regular motions (harmonic oscillators, plane waves, free particles, etc.) as zeroth- order approximations to physical systems, and account for weak nonlinear- ities perturbatively. We think of a dynamical system as a smooth system whose evolution we can follow by integrating a set of differential equations. The universality theory tells us that the zeroth-order approximations to strongly nonlinear systems should be quite different. They show an amaz- ingly rich structure which is not at all apparent in their formulation in terms of differential equations; instead, they exhibit self-similar structures which can be encoded by universality equations of a type which we will describe here.
    [Show full text]