MAT335H1F Lec0101 Burbulla

Similar documents
MAT335H1F Lec0101 Burbulla

Algorithms: Lecture 2

CHAOS AND DYNAMICS KELSEY MACE

2 Problem Set 2 Graphical Analysis

Zoology of Fatou sets

Warm-Up Problem. Please fill out your Teaching Evaluation Survey! Please comment on the warm-up problems if you haven t filled in your survey yet.

Math 354 Transition graphs and subshifts November 26, 2014

We are going to discuss what it means for a sequence to converge in three stages: First, we define what it means for a sequence to converge to zero

Infinity Unit 2: Chaos! Dynamical Systems

Discrete Mathematics and Probability Theory Spring 2014 Anant Sahai Note 20. To Infinity And Beyond: Countability and Computability

Julia Sets and the Mandelbrot Set

MATH 614 Dynamical Systems and Chaos Lecture 2: Periodic points. Hyperbolicity.

Chaos in the Dynamics of the Family of Mappings f c (x) = x 2 x + c

Julia Sets and the Mandelbrot Set

MAT1000 ASSIGNMENT 1. a k 3 k. x =

Mathematics 220 Workshop Cardinality. Some harder problems on cardinality.

MATH 521, WEEK 2: Rational and Real Numbers, Ordered Sets, Countable Sets

UNIVERSITY OF TORONTO Faculty of Arts and Science DECEMBER 2011 EXAMINATIONS. MAT335H1F Solutions Chaos, Fractals and Dynamics Examiner: D.

THE GOLDEN MEAN SHIFT IS THE SET OF 3x + 1 ITINERARIES

Lebesgue Measure and The Cantor Set

What is proof? Lesson 1

Topology Proceedings. COPYRIGHT c by Topology Proceedings. All rights reserved.

Discrete dynamics on the real line

ON SPACE-FILLING CURVES AND THE HAHN-MAZURKIEWICZ THEOREM

Math-2A Lesson 13-3 (Analyzing Functions, Systems of Equations and Inequalities) Which functions are symmetric about the y-axis?

Topological properties of Z p and Q p and Euclidean models

One-to-one functions and onto functions

Chaos and Liapunov exponents

Discrete Mathematics for CS Spring 2007 Luca Trevisan Lecture 27

Real Variables: Solutions to Homework 3

The Dynamics of Two and Three Circle Inversion Daniel M. Look Indiana University of Pennsylvania

A Short Review of Cardinality

CSE 311: Foundations of Computing. Lecture 26: Cardinality

PHY411 Lecture notes Part 5

Math 105A HW 1 Solutions

1 The Cantor Set and the Devil s Staircase

The dynamics of maps close to z 2 + c

MATH31011/MATH41011/MATH61011: FOURIER ANALYSIS AND LEBESGUE INTEGRATION. Chapter 2: Countability and Cantor Sets

v n+1 = v T + (v 0 - v T )exp(-[n +1]/ N )

Math 172 HW 1 Solutions

Stat 451: Solutions to Assignment #1

CITS2211 Discrete Structures (2017) Cardinality and Countability

Math 300: Final Exam Practice Solutions

arxiv: v2 [math.ca] 4 Jun 2017

Date: October 24, 2008, Friday Time: 10:40-12:30. Math 123 Abstract Mathematics I Midterm Exam I Solutions TOTAL

Countability. 1 Motivation. 2 Counting

Section Summary. Sequences. Recurrence Relations. Summations. Examples: Geometric Progression, Arithmetic Progression. Example: Fibonacci Sequence

CHAPTER 8: EXPLORING R

ECS 120 Lesson 18 Decidable Problems, the Halting Problem

Homework 1 Real Analysis

Math 203A - Solution Set 1

The coordinates of the vertex of the corresponding parabola are p, q. If a > 0, the parabola opens upward. If a < 0, the parabola opens downward.

Chapter 20. Countability The rationals and the reals. This chapter covers infinite sets and countability.

WEEK 7 NOTES AND EXERCISES

Outline Inverse of a Relation Properties of Relations. Relations. Alice E. Fischer. April, 2018

SETS AND FUNCTIONS JOSHUA BALLEW

Introductory Analysis I Fall 2014 Homework #5 Solutions

1 Partitions and Equivalence Relations

Solutions to Homework Assignment 2

Discrete Mathematics and Probability Theory Spring 2016 Rao and Walrand Discussion 6B Solution

Extrema and the Extreme Value Theorem

INFINITY: CARDINAL NUMBERS

We first review various rules for easy differentiation of common functions: The same procedure works for a larger number of terms.

MATH 102 INTRODUCTION TO MATHEMATICAL ANALYSIS. 1. Some Fundamentals

McGill University Math 354: Honors Analysis 3

Definition: A "system" of equations is a set or collection of equations that you deal with all together at once.

CHAOTIC UNIMODAL AND BIMODAL MAPS

Math 3361-Modern Algebra Lecture 08 9/26/ Cardinality

Contents. 1 Introduction to Dynamics. 1.1 Examples of Dynamical Systems

Math.3336: Discrete Mathematics. Cardinality of Sets

Iowa State University. Instructor: Alex Roitershtein Summer Exam #1. Solutions. x u = 2 x v

Singular Perturbations in the McMullen Domain

MATH 3300 Test 1. Name: Student Id:

Math 4603: Advanced Calculus I, Summer 2016 University of Minnesota Notes on Cardinality of Sets

Infinite Limits. Infinite Limits. Infinite Limits. Previously, we discussed the limits of rational functions with the indeterminate form 0/0.

consists of two disjoint copies of X n, each scaled down by 1,

Slope Fields: Graphing Solutions Without the Solutions

the time it takes until a radioactive substance undergoes a decay

MATH 1910 Limits Numerically and Graphically Introduction to Limits does not exist DNE DOES does not Finding Limits Numerically

Chapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics

Chapter 8 ~ Quadratic Functions and Equations In this chapter you will study... You can use these skills...

Lecture 3: Sizes of Infinity

Sets are one of the basic building blocks for the types of objects considered in discrete mathematics.

MS 3011 Exercises. December 11, 2013

The Two Faces of Infinity Dr. Bob Gardner Great Ideas in Science (BIOL 3018)

Cardinality of Sets. P. Danziger

LECTURE 22: COUNTABLE AND UNCOUNTABLE SETS

CS 125 Section #10 (Un)decidability and Probability November 1, 2016

Finding numerically Newhouse sinks near a homoclinic tangency and investigation of their chaotic transients. Takayuki Yamaguchi

ADVANCED CALCULUS - MTH433 LECTURE 4 - FINITE AND INFINITE SETS

Chapter One. The Real Number System

Chapter 1A -- Real Numbers. iff. Math Symbols: Sets of Numbers

Contents Ordered Fields... 2 Ordered sets and fields... 2 Construction of the Reals 1: Dedekind Cuts... 2 Metric Spaces... 3

UNM - PNM STATEWIDE MATHEMATICS CONTEST XLII. November 7, 2009 First Round Three Hours

Math 115 Spring 11 Written Homework 10 Solutions

... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré

Extended Essay - Mathematics

Identifying the Graphs of Polynomial Functions

The First Derivative Test

Consider the equation of the quadratic map. x n+1 = a ( x n ) 2. The graph for 0 < a < 2 looks like:

Transcription:

Fall 2011

Q 2 (x) = x 2 2 Q 2 has two repelling fixed points, p = 1 and p + = 2. Moreover, if I = [ p +, p + ] = [ 2, 2], it is easy to check that p I and Q 2 : I I. So for any seed x 0 I, the orbit of x 0 for Q 2 will stay in I, but it will not converge to p = 1, or p + = 2, unless x 0 just happens to be on an orbit that is eventually fixed at one of these two fixed points. As an example, take x 0 = 0.2 and consider n = 200 iterations of x 0 under Q 2. Two graphical representations of this orbit are represented on the next slide. A Chaotic Orbit Not only is the orbit repelled from p = 1, it also seems not to be attracted to any periodic cycle at all; it seems to be completely random, or chaotic.

Periodic Points of Q 2 Paradoxically, Q 2 has infinitely many periodic points, of all possible periods. We can see why by considering the graphs of Q 2 n, for a few values of n. Figure: Q 2 Figure: Q 2 2 Figure: Q 3 2 Figure: Q 4 2 In general, the function Q 2 has at least 2 n periodic points of period n in the interval I. Thus the proliferation of periodic points for Q c that we have already described, for 5/4 < c < 1, seems to increase without bound. But by c = 2 all the periodic points seem to be repelling. We shall investigate this more fully... later.

An Example from the Logistic Family: F 4 (x) = 4x(1 x) F 4 has two repelling fixed points, p = 0 and p = 3 4. Moreover, if J = [0, 1], it is easy to check that F 4 : J J. Again, for any seed x 0 J, the orbit of x 0 for F 4 will stay in J, but it will not converge to p = 0 or p = 3/4 unless x 0 is on an orbit that is eventually fixed at one of these two fixed points. As an example, take x 0 = 0.1 and consider n = 200 iterations of x 0 under F 4. Two graphical representations of this orbit are represented on the next slide. Another Chaotic Orbit Not only is the orbit repelled from both p = 0 and p = 2/3 it also seems not to be attracted to any periodic cycle at all; as with the previous example, it seems to be completely random, or chaotic.

Periodic Points of F 4 As with Q 2, F 4 also has infinitely many periodic points, of all possible periods. By considering the graphs of F4 n, for a few values of n you can see behaviour very similar to the quadratic family. Figure: F 4 Figure: F 2 4 Figure: F 3 4 Figure: F 4 4 Q 2 and F 4 are both quadratic functions, one opens upwards and one downwards. But both share one important feature: every subsequent iteration of Q 2 or F 4 doubles the number of valleys or hills in the graph of Q 2 n or F 4 n on the interval I or J, resp.

The Range of Q c if c < 2. Consider the graph of Q c on [ p +, p + ] for c < 2. The minimum value of Q c is at x = 0; it is y = c. If c < 2, then c < p + since c < p + = 1 + 1 4c 2 2c > 1 + 1 4c 1 2c > 1 4c (1 + 2c) 2 > 1 4c 4c 2 + 4c + 1 > 1 4c c 2 + 2c > 0 c < 2 or c > 0 This means that for c < 2 the image of the interval [ p +, p + ] under Q c is no longer contained in [ p +, p + ]. Let c = 2.5 For the Rest of This Section. Let c = 2.5; then p + = 1 + 11 2 2.16. To the right is the graph of Q 2.5 on the interval [ [ p +, p + ] = 1 + 11, 1 + ] 11. 2 2 Notice that the parabola dips a little below the black square, which is determined by the four points (±p +, ±p + ).

Orbits That Escape to Infinity Q 2.5 (x) = p + x 2 2.5 = 1 + 11 2 x = ± 4 11. 2 These two points are eventually fixed points. If 4 11 4 11 < x 0 <, 2 2 then the orbit of x 0 under Q 2.5 will escape to infinity, since x 1 < p +. Some orbits escape after 2 or 3 iterations: Figure: x 0 = 1.5, n = 2 Figure: x 0 = 1, n = 3

Here s an orbit that escapes after 17 iterations: take x 0 = 1.8. But not all orbits escape to infinity: the orbit of x 0 could be an eventually fixed point, ending up at p + or p. Let I = [ p +, p + ]; let Λ = {x I Q n 2.5(x) I for all n}. Λ is the set of all points in I for which the orbit under Q 2.5 never leaves I. What kind of set is Λ? One way to describe Λ is to describe all the points in I that are not in Λ; that is, describe the set of all points x 0 I such that the orbit of x 0 under Q 2.5 goes to infinity. If the orbit of x I under Q 2.5 goes to infinity then there is a least value of n for which Q 2.5 n (x) / I. Let A n = {x I n is the least value of k such that Q k 2.5(x) / I }. As we have already calculated, A 1 = ( 2.5 p +, ) 2.5 p +. A 1 is an open interval since Q 2.5 sends both endpoints to p +, which is an eventually fixed point. Figure: A 1

We can see the sets A 2 and A 3 by plotting Q 2 2.5 and Q3 2.5 on I. Figure: A 2 Figure: A 3 A 2 is the union of two disjoint open intervals and A 3 is the union of four disjoint open intervals. Remarks About Λ and A n 1. A n is the union of 2 n 1 disjoint open intervals. 2. Q 2.5 (A n ) = A n 1 ; equivalently, A n = Q 1 2.5 (A n 1). 3. n=1a n is the set of all points in I whose orbit under Q 2.5 escapes to infinity. 4. Λ = I n=1a n is the set of all point in I whose orbits remain in I. 5. Λ must be a closed set, since it is the complement of a union of open intervals. 6. Not so obviously: Λ itself contains no open interval. It is a collection of isolated points. That is, Λ is a closed set with total length zero! Such a set is called totally disconnected, or dust.

Λ Contains No Open Interval; Proof by Contradiction Suppose J = (a, b) is an open interval in Λ. J is on the right side of A 1 or on the left side; assume it is on the right side. 0 2.5 p + 2.5 p+ p + a b p + We have 2.5 p + 0.58 < a < b. For x > 0, Q 2.5 is an increasing function, so Q 2.5 (J) = (Q 2.5 (a), Q 2.5 (b)) = (a 2 2.5, b 2 2.5). Thus Q 2.5 (J) is another interval in Λ (why?) with length b 2 a 2 = (b + a)(b a) > 1.16 (b a). Repeat k times: for each k, Q 2.5 k (J) is another interval in Λ with length greater than 1.16 k (b a). For k large enough, 1.16 k (b a) will be bigger than the length of I, contradicting the fact that Λ I. An Algorithmic Construction of Λ (Not to Scale!) 1. Start with I = [ p +, p + ]: 2. Take out A 1 : 3. Take out A 2 : 4. Take out A 3 : Repeat until every open subset A n has been taken out. What s left is Λ. Just such a set was introduced in 1883 by Georg Cantor; it is called Cantor s middle-thirds set, or simply Cantor s set. NOTE: The arguments in this section are actually valid if c (5 + 2 5)/4 2.368, for which c p + 1/2, and for which Q c(x) 1, if x Λ A 1. This last condition is what Devaney uses in his proofs. The results of this section are true for all c < 2, but harder to prove for (5 + 2 5)/4 < c < 2.

Definition of the Cantor Middle-Thirds Set: A Recipe 1. Start with the interval [0, 1]. 2. Remove the middle third (1/3, 2/3), leaving two closed intervals left, [0, 1/3] and [2/3, 1], each of length 1/3. 3. From each of these two intervals remove the open middle thirds, (1/9, 2/9) and (7/9, 8/9), leaving the four closed intervals [0, 1/9], [2/9, 3/9] and [6/9, 7/9], [8/9, 9/9] left, each of length 1/9. 4. Repeat this process: remove the open middle third from each of the previous closed intervals. 5. K is the set of points remaining in [0, 1] in the limit as this process is repeated over and over without end. 6. K φ since it contains the endpoints of every closed interval remaining at each step. Picturing the Cantor Set Here are the first few iterations in the construction of the Cantor Set. Figure: The Cantor Middle-Thirds Set Note that the construction of K is very similar to that of Λ, as described in the last section. In fact, K and Λ have many common properties, which are easier to prove for K, but are also true for Λ.

Properties of the Cantor Set 1. K is a closed subset of [0, 1], because it is the complement of a union of open sets. 2. K is totally disconnected because, as with Λ, you can show that K contains no open intervals. 3. a K if and only if there is a ternary expansion of a such that a = 0.s 1 s 2 s 3 s 4... with s i = 0 or 2. That is, s i 1. This is base 3 arithmetic! 4. K is uncountable. That is, K cannot be put into a one-to-one correspondence with N, the set of natural numbers. Ternary Expansions Let a [0, 1]. In terms of base 3, a can be written as an infinite series s i a = 3 i = s 1 3 + s 2 3 2 + s 3 3 3 +, i=1 with s i {0, 1, 2}. We call 0.s 1 s 2 s 3 s 4... the ternary expansion of a. After all, there is nothing special about decimal expansions; they are base 10. Any base is just as valid as any other. Indeed, in computer science base 2 is very common. In terms of base 2, we could get a binary expansion of a, in which all of its digits would be 0 or 1.

Examples of Ternary Expansions Example 1: The ternary expansion of 1/4 is 0.02020202... since 0 3 + 2 3 2 + 0 3 3 + 2 3 4 + 0 3 5 + 2 3 6 + = 2 9 n=0 n=0 ( ) 1 n = 2 1 9 9 1 1/9 = 1 4. Example 2: The ternary expansion of 5/26 is 0.012012012... since 0 3 + 1 3 2 + 2 3 3 + 0 3 4 + 1 3 5 + 2 3 6 + = 1 ( ) 1 n + 2 ( ) 1 n 9 27 27 27 = 1 9 n=0 1 1 1/27 + 2 1 27 1 1/27 = 3 26 + 2 26 = 5 26. Numbers With Two Ternary Expansions In the decimal system some numbers have two decimal expansions, eg. 1.0000.... = 0.9999.... Some numbers can also have two ternary expansions. For example, since 1 3 = 0.1000... and 1 3 = 0.02222..., 0.0222.... = 0 3 + 2 3 2 + 2 3 3 + 2 3 4 + = 2 9 Similarly, check that 8 9 n=0 = 0.220000... as well as 0.212222.... ( ) 1 n = 2 1 3 9 1 1/3 = 1 3.

Ternary Expansions and The Cantor Set The ternary digits of a [0, 1] tell you at each step which third a is in: s 0 = 0 s 0 = 1 s 0 = 2 s 0 s 1 = 00 01 02 10 11 12 20 21 22 Etc, etc. So a K if all its ternary digits are 2 or 0; a / K if one of its ternary digits is 1. For example 1/4 = 0.020202 K, but 5/26 = 0.012012... / K. But wait: 1/3 = 0.1 K and it has a ternary digit 1. However, as described above, 1/3 can be written as 0.022222.... All in all, we can say: a K if and only if there is a ternary expansion of a such that all of its digits are 0 or 2. The Cantor Set is Uncountable Suppose K is countable. Then there is a bijection Φ : N K. Let k n = Φ(n). List all the elements of K as k 1, k 2, k 3,... : k 1 = 0.s 1 1 s1 2 s1 3 s1 4 s1 5... with s1 i {0, 2} k 2 = 0.s1 2s2 2 s2 3 s2 4 s2 5... with s2 i {0, 2} k 3 = 0.s1 3s3 2 s3 3 s3 4 s3 5... with s3 i {0, 2} k 4 = 0.s1 4s4 2 s4 3 s4 4 s4 5... with s4 i {0, 2} k 5 = 0.s1 5s5 2 s5 3 s5 4 s5 5... with s5 i {0, 2}......... { 0, if s i Define k = 0.s 1 s 2 s 3 s 4 s 5... such that s i = i = 2 2, if si i = 0. Then k K, because all its ternary digits are 0 or 2; but k is not on the above list because it differs with each ternary expansion in at least one digit. So K is uncountable.

What Is An Orbit Diagram? An orbit diagram is an attempt to describe in one picture the dynamics of a family of functions F λ as the parameter λ changes. Values of λ will be on the horizontal axis; values of x n will be on the vertical axis. But not all values of x n : we ignore the first 100 points (say) in the orbit and only plot the values of x n for n > 100. The hope is that if the orbit converges to a point or a periodic cycle, it should have reached it after 100 iterations. So on the orbit diagram we might plot the values of x n, for 100 < n < 300. The choice of x 0 is not random: for x 0 we pick a non-degenerate critical point of F λ, i.e. a point x 0 such that F λ (x 0) = 0 but F λ (x 0) 0. (The reason for this choice of x 0 won t become apparent until Chapter 12.) For the quadratic family, Q c (x) = x 2 + c, the non-denerate critical point is x 0 = 0; for the logistic family, F λ (x) = λx(1 x), the non-degenerate critical point is x 0 = 1/2. The Orbit Diagram for the Quadratic Family Q c So to produce an orbit diagram for the quadratic family we plot, for a given value of c, the values of x n, for 100 < n < 300, in a vertical line. Of course, we always have to pick a finite interval on the vertical axis, so any orbit that escapes to infinity would not show up on the orbit diagram. Based on our previous work we know that for c > 1/4, the orbit of x 0 = 0 goes to infinity, and for c < 2, the orbit of x 0 = 0 also goes to infinity. But for 2 c 1/4, the orbit of x 0 = 0 under Q c is completely contained in the interval [ p +, p + ] [ 2, 2]. Although the production of an orbit diagram is experimental math, our previous work indicates we should expect the orbit to be eventually fixed if 3/4 < c < 1/4, and eventually a 2-cycle if 5/4 < c < 3/4. The orbit diagram is on the next slide....

The Orbit Diagram for Q c (x) = x 2 + c, 2 c < 1/4 Remarks Concerning the Orbit Diagram of Q c 1. As c decreases past 3/4, more period doubling occurs: a 2-cycle is followed by a 4-cycle, which is followed by an 8-cycle, etc. 2. In the period-3 window, cycles of order 6, 12,..., appear as c decreases. 3. The orbit diagram seems self-similar: parts of the diagram resemble the whole diagram. 4. The lines in the orbit diagram that connect periodic cycles do so continuously. 5. For many values of c the orbit of x 0 = 0 under Q c is not attracted to any periodic cycle whatsoever.

The Orbit Diagram for the Logistic Family, F λ Recall (or check) the following features about the logistic family: 1. F λ (x) = λx(1 x) 2. F λ : [0, 1] [0, 1] for 1 λ 4. 3. p = 1 1/λ is an attracting fixed point if 1 < λ < 3. 4. F λ has an attracting 2-cycle if 3 < λ < 1 + 6. λ + 1 ± λ 2 2λ 3 2λ 5. x 0 = 1/2 is a non-degenerate critical point of F λ. The orbit diagram for F λ is on the next slide.... The Orbit Diagram for F λ (x) = λx(1 x), 2.4 λ 4

Geometry Behind Some Features of the Orbit Diagram Consider the following six graphs of Q c, for 2 < c < 1/4. Each graph includes the line y = x and the square determined by the vertices (±p +, ±p + ). Figure: c = 1/4 Figure: c = 0 Figure: c = 3/4 Figure: c = 1 Figure: c = 2 Figure: c = 2.2 The first three diagrams illustrate the tangent bifurcation at c = 1/4, the appearance of an attracting fixed point, and the changing of the fixed point to a repelling fixed point at c = 3/4; the second three diagrams illustrate the cases of an attracting 2-cycle, chaotic orbits, and finally, orbits that can escape to infinity.

Now Consider the Graphs of Q 2 c for 3/4 < c < 1.65 The following three graphs show that in the marked box, lower left, the graphs of Q 2 c for 5/4 < c < 3/4 are similar to the graphs of Q c for 3/4 < c < 1/4. Figure: c = 3/4 Figure: c = 1 Figure: c = 5/4 Figure: c = 1.3 Figure: c = 1.546 Figure: c = 1.65 The graphs of Q 2 c in the indicated boxes resemble very much the last three graphs of Q c, in the previous sequence of graphs. Thus in the orbit diagram of Q c we would expect two bits of the diagram for 1.65 < c < 0.75 to resemble very much the whole orbit diagram for 2 c 0.25, one bit on each part of the 2-cycle.

Period Doubling What the previous graphs illustrate is that just as the fixed point for Q c becomes repelling an attracting 2-cycle appears, namely the fixed points for Q 2 c. Similarly, just as the attracting 2-cycle for Q c becomes repelling, an attracting 4-cycle appears, namely the fixed points of Q 4 c. But the graphs of Q 2 c, Q 4 c, Q 8 c,... are getting steeper and steeper, so the bifurcations occur for values of c closer and closer together. This behaviour is not peculiar to the period-1 window; the same pattern of period doubling occurs in the other windows in the orbit diagram, though some of these windows are very thin compared to the period-1 window. Some Graphs of Q 3 c for 1.8 < c < 1.7 For example, the graphs of Q 3 c pass the line y = x in three different places for 1.8 < c < 1.75. Thus we expect the orbit diagram in the 3-window to resemble three copies of the whole diagram. Figure: c = 1.7 Figure: c = 1.78 Figure: c = 1.8