Discovering Dynamical Systems Through Experiment and Inquiry
Discovering Dynamical Systems Through Experiment and Inquiry
Dynamical Systems
Through Experiment
and Inquiry
Textbooks in Mathematics
Series editors:
Al Boggess, Kenneth H. Rosen
Advanced Calculus
Theory and Practice, Second Edition
John Srdjan Petrovic
Nonlinear Optimization
Models and Applications
William P. Fox
Linear Algebra
James R. Kirkwood, Bessie H. Kirkwood
Real Analysis
With Proof Strategies
Daniel W. Cunningham
Linear Algebra
What You Need to Know
Hugo J. Woerdman
Reasonable efforts have been made to publish reliable data and information, but the author and pub-
lisher cannot assume responsibility for the validity of all materials or the consequences of their use.
The authors and publishers have attempted to trace the copyright holders of all material reproduced
in this publication and apologize to copyright holders if permission to publish in this form has not
been obtained. If any copyright material has not been acknowledged please write and let us know so
we may rectify in any future reprint.
Except as permitted under U.S. Copyright Law, no part of this book may be reprinted, reproduced,
transmitted, or utilized in any form by any electronic, mechanical, or other means, now known or
hereafter invented, including photocopying, microfilming, and recording, or in any information stor-
age or retrieval system, without written permission from the publishers.
For permission to photocopy or use material electronically from this work, access www.copyright.
com or contact the Copyright Clearance Center, Inc. (CCC), 222 Rosewood Drive, Danvers, MA
01923, 978-750-8400. For works that are not available on CCC please contact mpkbookspermis-
[email protected]
Trademark notice: Product or corporate names may be trademarks or registered trademarks and are
used only for identification and explanation without intent to infringe.
Preface xi
2 Sequences 11
2.1 Introduction to Sequences . . . . . . . . . . . . . . . . . . . . 11
2.2 Convergence of Sequences . . . . . . . . . . . . . . . . . . . . 12
2.3 The Squeeze Theorem . . . . . . . . . . . . . . . . . . . . . . 15
2.4 Arithmetic Limit Theorems . . . . . . . . . . . . . . . . . . . 18
2.5 Bounded and Unbounded Sequences . . . . . . . . . . . . . . 19
2.6 Subsequences . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.7 Liminfs and Limsups . . . . . . . . . . . . . . . . . . . . . . 24
2.8 Cauchy Sequences . . . . . . . . . . . . . . . . . . . . . . . . 29
vii
viii Contents
5 Bifurcations 69
5.1 What is a Bifurcation? . . . . . . . . . . . . . . . . . . . . . 69
5.2 Introduction to Bifurcation Diagrams . . . . . . . . . . . . . 71
5.3 The Tangent Bifurcation . . . . . . . . . . . . . . . . . . . . 72
5.4 The Period Doubling Bifurcation . . . . . . . . . . . . . . . . 76
Bibliography 193
Index 195
Preface
xi
xii Preface
job to ensure that they see those failures as learning opportunities, rather than
disappointments. Be patient, be kind, but hold the students accountable. You
can expect some resistance to the method for a few weeks. We have had
students complain: “He doesn’t teach anything; he just stands there asking
questions.” This is normal, but once a student experiences the joy of self-
discovery, there is no turning back.
If you are new to using IBL methods in the classroom, some help-
ful resources are available from the Academy of Inquiry Based Learn-
ing (http://www.inquirybasedlearning.org/) and the Legacy of R.L. Moore
(http://legacyrlmoore.org/). Both of these organizations offer workshops to
help get instructors involved with IBL. The Academy of Inquiry Based Learn-
ing, in particular has online resources and many active communities for in-
structors new to the method.
We refer to the unique approach to teaching mathematics used in this book
as ECAP for Explore, Conjecture, Apply, and Prove. We believe that ECAP
provides students with a deeper understanding of dynamical systems because
it mimics the actual practice of mathematics. In general, each section begins
with exercises guiding students through explorations of the featured concept in
dynamical systems and concludes with exercises that help the student formally
prove the result.
We chose to put the explorations first (as opposed to end-of-chapter
projects) because in our experience, computer explorations provide students
with an enhanced understanding of the mathematical concepts by allowing
improved visualization that is both dynamic and engaging. Moreover, the
ability to manipulate the equations and/or underlying mathematical aspects
allows students to grasp what aspects are essential for a given phenomenon
and what are not. This exploration, combined with some guided questions,
gives students the confidence to conjecture the formal theorems being formed.
The proof subsection then completes the student’s journey from the discovery
of the mathematics to the final step of proving the validity of the result.
Let’s look at this structure in more detail to understand how this text
might be used during the semester. Most sections of the text are divided into
Explore, Conjecture, Apply, and Prove portions. We find this to be an
effective approach to guide the young mathematician toward the depth of
understanding expected in an IBL classroom.
The Explore sections ask leading questions to introduce a concept and
help students see what may be true. Many of the explorations use nu-
merical activities that we have designed specifically to highlight the con-
cepts and ideas being explored. Each accompanying activity is available on
IBLdynamics.com . There you will find a collection of interactive apps de-
signed to get students thinking about the topics at hand. Students can discover
examples and counterexamples through manipulations built into the software.
This is a unique feature of the text, and one that has proven essential for stu-
dents taking the course in the past.
Preface xiii
On the other hand, we did see a topical theme emerge that we believe
is unique to a dynamical systems text: symbolic dynamics. While symbolic
dynamics is a fairly standard topic in an undergraduate dynamics text, we
have tried to emphasize it in a way that is more detailed and inclusive than is
typically the case. This is partly because the ideas and proofs in this area are
accessible to undergraduates and approachable using the ECAP methodology.
Thus we chose to conclude this book with a chapter dedicated to a general
discussion of the topic.
Finally, we have chosen to include multiple sections on important ideas
from analysis and topology independent from their application to dynamics.
We feel that most students need a refresher (at the least) on these topics
and the study of dynamical systems provides a wonderful motivation and/or
introduction to these ideas.
In piloting this material, we were able to cover the first 10 chapters of this
book in a single semester. If you don’t need to cover the analysis material, you
will most likely be able to cover at least some of the material in Chapters 11–
13. We believe that the pairings of Chapters 11 and 12 or of Chapters 12 and 13
work well. Chapter 11 is not a prerequisite for Chapters 12 and 13 though the
material concerning the Hénon map in Chapter 11 nicely illustrates the theme
of Chapter 12.
Chapter 1 defines discrete dynamical systems and discusses numerical
methods for tracking the orbits of points in a system. It also introduces the
students to some of the tools on IBLdynamics.com . Students can also work
on graphical iteration, and see how dynamics can apply to population models.
In our experience, the upper-level students taking this class have little
experience with the analysis of sequences and other tools from analysis that
arise when studying dynamical systems. Thus, we have chosen to include much
of this material in this book. Chapter 2 provides a fairly complete coverage of
sequences for this reason. You may choose to skip this chapter if that is not
the case with your students.
Chapter 3 defines fixed and periodic points. The number of explorations
on iteration increases here, as we work to build a students’ intuition about
attracting and repelling points. They are challenged to make conjectures about
attracting and repelling fixed points, studying linear, quadratic, and logistic
maps. Population models are revisited and some important theorems about
fixed and periodic points are proven.
Much of Chapter 4 covers analysis topics with the primary focus being the
implicit function theorem and fixed point theorems. The students are given
more opportunity to come up with counterexamples, and the explorations
reflect this. The chapter concludes with a proof of the hyperbolic fixed points
theorem. Section 4.2 can be safely omitted if your students are familiar with
the implicit function theorem.
Bifurcations are covered in Chapter 5. A large number of explorations are
included at the beginning of the chapter to drive students toward conjecturing
Preface xv
Message to Students
We could have structured the content of this text like a typical mathematics
book, but we’ve chosen not to, and for a specific reason. This book is designed
to help you learn to think like a mathematician. We want you to go beyond
memorizing theorems and calculating results. We want you to explore, to
discover, and to prove the results for yourself. Most students are brought up
with the idea that mathematics is entirely about showing what is true. It’s
actually the case that the much more important step is showing when and
why something is true. In this text, we plan for you to achieve exactly that.
Most every section starts with exploration activities. You will be asked to
work things out on paper or use activities at IBLdynamics.com to explore
some topic in dynamical systems. It is expected that you try everything. Take
notes, notice patterns, and don’t be afraid to make mistakes. If you see some-
thing that seems to be generally true, be ready to make that conjecture. Be
prepared to share your insight with the class. And be ready to explain why you
think your insight is true. Your instructor will guide you through the process,
but ultimately, the work must be your own. You will learn that it is acceptable
to fail, and that making a false conjecture, or an error in explanation, is just
part of learning to be a mathematician. You will learn that these failures can
be productive. If we learn from our errors and move forward, an error is not
a failure, it is an opportunity to truly understand something new.
A course like this returns dividends proportional to your efforts. It is much
easier to be given an algorithm, and to execute that algorithm, than to struggle
for the results yourself. It is difficult, but struggling is essential for understand-
ing the mathematics and for personal growth. We sometimes tell our students
to “embrace the stuck” because all new and important ideas originate here.
The joy of discovery, the frustration of exploration, and the satisfaction of
collaboration are all important aspects of mathematics. You will explore. You
will conjecture. You will prove results. And you will apply results to a variety
of related problems. All of these aspects combine when one is creating new
mathematics. It is the intention of the authors that this course will prepare
you, not just for an understanding of dynamical systems, but for exploration
into any area of mathematics you find appealing.
Acknowledgments
We have received the help and support of many. But foremost on that list are
our colleagues in the MCS Department at Gustavus Adolphus College who
have listened to our musings, critiqued our writing, and occasionally covered
Preface xvii
our classes. We would be remiss if we didn’t single out Mike Hvidsten, who
encouraged us to take this project from a collection of classroom activities
and demonstrations to a full-fledged textbook.
We had the help of two wonderful Gustavus students, Emilee Mason and
Zach Dawson. Emilee, who “beta-tested” this project in our Discrete Dynami-
cal Systems course, was invaluable in her proofreading and editing duties. She
was especially helpful in finding all of the places where we omitted commas! If
any are still missing, it is our fault and not hers. Zach Dawson was a lifesaver!
A computer science major, Zach translated all of Tom’s poor Mathematica
code into clean and functioning JavaScript. We would still be struggling with
the basics of JavaScript if it weren’t for Zach.
Tom would also like to thank Francis Su and his colleagues at Harvey
Mudd College. During the infancy of this project, Tom spent a week visiting
Francis picking his brain about writing, inquiry-based learning, and life. As is
often the case, these latter conversations were the most valuable.
Jeff would like to thank Dr. Steven Clontz at University of Southern Al-
abama, and the late Dr. Frank Sturm, for their endless late-night discussions
about the value of inquiry-based learning and the nuts and bolts of how to
make it effective.
Various aspects of this work, including funding for Emilee and Zach, were
possible through a generous gift to Gustavus Adolphus College by the family
of Clifford M. Swanson.
1
An Introduction to Dynamical Systems
{xn }∞
n=0 = {x0 , x1 , . . . }
The space X is called the state space of the dynamical system, and for most of
this book, X will be the real numbers R. What is important here is that it is
very easy to compute the sequences generated by discrete dynamical systems.
Starting with an initial condition x0 , we just apply the function f to compute
x1 . In other words,
f (x0 ) = x1 .
The subscript can be thought of as a discrete time measurement. If time is
measured in years, then x0 is the value of x in year 0, x1 is the value of x in
year 1, and so on.
1
2 An Introduction to Dynamical Systems
x1 = f (x0 ) = f (12) = 6
and
x2 = f (x1 ) = f (6) = 3
and so on. Tables, such as the one inset into figure 1.1, are a natural way
to show sequences generated by iteration. The first column is the iteration
number, and the second column is the sequence of iterates of f with x0 = 12.
The plot in figure 1.1 shows the same iterations but in a graphical format.
Here the horizontal axis is the iteration number (n), and the vertical axis, the
sequence element (xn ).
Notice that each entry in the second column of the table in figure 1.1 is
f applied to the entry immediately above it. This is iteration and the focus
of our study. Given a dynamical system, we would like to describe all of the
possible sequences that are generated by the given function. We call these
sequences orbits. The example in figure 1.1 shows the first six terms of the
orbit of x0 = 12 under iteration by f (x). We might describe this particular
orbit by saying
FIGURE 1.1
The first six iterates of f (x) = x/2 with an initial condition of x0 = 12
displayed as both a table and a graph.
You probably already see that every orbit of this dynamical system, re-
gardless of initial condition, converges to 0 and by doing so you’ve completely
described this dynamical system.
Explore
xn+1 = x2n .
For each of the initial conditions listed below, compute the first 10 iterates of
their orbit. Make both a table and a plot like those in figure 1.1. In each case,
describe the orbit in your own words. Look for commonalities and differences
in the different orbits and try to categorize the different types. Suggestion:
Doing this by hand gets old pretty quick. However, doing this with a spread-
sheet is very easy. After doing the first couple by hand, think about how you
might implement this dynamical system in a spreadsheet program. This can
also be used to create plots such as that show in figure 1.1.
xn+1 = x2n − 1.
4 An Introduction to Dynamical Systems
Use the same initial conditions that you used in exploration 1. Describe each
of the orbits in your own words. Look for commonalities and differences in the
different orbits and try to categorize the different types.
FIGURE 1.2
A plot showing the graphical analysis of xn+1 = x2n with an initial condition
of x0 = 0.9.
Explore
- -
- -
-
-
- -
- -
- -
-
-
- -
FIGURE 1.3
Figures for practicing graphical analysis.
Apply
where
• V (p) is called the viability and measures the percentage of newborns
in the species that survive to reproductive age, and
• λ(p) is the mean reproductive rate and measures the mean number
of newborns per adult.
Application 14. Use the methods described in this chapter to play with
your model. Try a variety of reasonable parameter values for the unassigned
parameters. We suggest using a spreadsheet to do this, but feel free to use
whatever tool you find appropriate. Describe some of the behaviors that your
model produces. Below are a few questions that might be helpful to you as
you do this. Do not limit yourself to these questions. Be creative!
1. Is the behavior that you observe realistic? For example, does it ever
lead to negative population numbers. If it is not realistic, think
about modifying it and repeating this process.
2. Does the population approach some fixed value? Does it oscillate
between multiple values? Does it do something more complicated?
10 An Introduction to Dynamical Systems
Explore
Let’s begin with a fun example, the Fibonacci sequence. The first two
terms of the Fibonacci sequence are F0 = 1 and F1 = 1. Each subsequent
term is determined by the recurrence
Fn = Fn−1 + Fn−2
so that each term is simply the sum of the previous two terms.
Exploration 1. Compute Fn for n = 0, . . . , 10.
It is natural (and useful) to think of a sequence as simply an ordered list of
things (usually numbers). As mathematicians, we need to turn this intuitive
idea into something much more precise.
11
12 Sequences
It is common practice to omit the n = 1 and the ∞ from this notation and
simply write {xn } when referring to a sequence when no confusion will arise.
In this more sophisticated view, a sequence is a function that assigns a
value in the space X to each natural number via n 7→ xn . This gives a list of
numbers with the position of the number in the list corresponding to n. The
space X is the set of all outputs of this function (or codomain) and for us X
will almost always be the real numbers R.
Exploration 2. Compute the first 5 elements of each of the sequences below.
Describe the behavior of each sequence in your own words.
1
1. xn = 1 −
n
2. yn = (−1)n
n
4
3. zn =
5
4
4. w0 = 0, wn+1 = wn
5
n2 + n − 1
5. an =
n+4
2
n +n−1
6. bn =
n2 + 4
lim xn = L.
n→∞
If a sequence does not converge to a limit L, we say that the sequence diverges.
Convergence of Sequences 13
This definition should sound very familiar to you. It is almost exactly the
definition of a limit that you learned in calculus. In particular, it is almost
word for word the definition of a “limit at infinity.” How does definition 2.2
formalize our intuitive concept of convergence? If ε is how closely you want to
approximate the limit L, then N is how far you have to go out in the sequence
to obtain this level of precision.
|A − B|. (2.1)
Read expression 2.1 as “the distance between A and B.” The inequality
|x − 4| < 3
then reads “the distance between x and 4 is less than 3.” When thought
of this way, it is clear that this is the same as
1<x<7
|xn − L| < ε
in definition 2.2 reads “the distance between xn and L is less than ε.”
This naturally leads to
L − ε < xn < L + ε
Explore
Exploration 5. Reread definition 2.2. For the values of ε that you chose in
exploration 4, directly compute the value of N so that
|xn − L| < ε.
How do each of these computed values compare with your answers in Explo-
ration 4?
Apply
It is often the case that the techniques you used in calculus to compute
limits at infinity apply to sequences. That is the case with many of the ex-
ercises in this section. It might be helpful to review this material in general,
and L’Hopital’s rule in particular.
n
4
Application 6. Compute the limit of the sequence zn = and use
5
definition 2.2 to prove that it converges to this limit value.
The Squeeze Theorem 15
2n + 4
Application 8. Compute the limit of the sequence zn = and use
5n + 1
definition 2.2 to prove that it converges to this limit value.
2
Application 9. Compute the limit of the sequence zn = and use
5n + 1
definition 2.2 to prove that it converges to this limit value.
Explore
cn = wan + (1 − w)bn .
What can you say about the relationship of the values cn relative to the values
of an and bn ?
Conjecture
Conjecture 15. Based on your work above, complete the following statement
of the Squeeze Theorem.
Apply
Application 16. Use the Squeeze Theorem to compute the following limits.
cos n
1. lim
n→∞ n
2
n + sin n cos n
2. lim
n→∞ n2
The Squeeze Theorem 17
Application 17. Use the Squeeze Theorem to show that the following se-
quences converge.
cos(n2 )
1. xn =
n4
sin(n)
e
2. xn =
n2
(−1)n
3. xn =
n
Follow the steps below to prove that
Application 18. Let sn = n1/n − 1. To what value should the sequence {sn }
converge in order to prove equation 2.2?
Application 24. Apply the Squeeze Theorem to prove the result in equa-
tion 2.2.
Prove
Proof 25. Use Definition 2.2 to prove the Squeeze Theorem as stated below.
then
lim cn = L.
n→∞
18 Sequences
lim an = A
n→∞
and
lim bn = B.
n→∞
Conjecture
Apply
Application 31. Find the values to which the following sequences converge.
Do these results agree with your conjectures 26 through 30?
3n2 − 1
1. an =
10n + 5n2
Bounded and Unbounded Sequences 19
ln(n + 1)
2. bn =
ln(1 + 4n)
3. cn = e−n
Application 32. Assume {an } and {bn } are both real sequences. Provide
counterexamples to each of the following statements.
1. If limn→∞ an = 0, then limn→∞ (an bn ) = 0.
2. If limn→∞ (an bn ) = 0, then either limn→∞ an = 0 or limn→∞ bn =
0.
3. If limn→∞ (an +bn ) does not exist, then either {an } diverges or {bn }
diverges.
an+1
4. If limn→∞ an = 1, then {an } converges.
Do these problems suggest revising any of your conjectures 26 through 30?
Prove
Definition 2.4 A sequence {xn } is said to have an infinite limit if for every
M > 0 there exists N > 0 such that n > N implies xn > M. We denote this
by
lim xn = ∞.
n→∞
Explore
Exploration 35. List some techniques that you might use to show that a
sequence is monotone?
Conjecture
Apply
n2 + 3
lim = ∞.
n→∞ n − 1
n2 + 3
lim = ∞.
n→∞ n + 1
13n
Application 42. Is the sequence an = monotonic?
n!
Application 43. If a sequence of positive real numbers is strictly decreasing,
must it converge to zero?
Prove
1. S = (0, 1)
2. S = [0, 1]
3. S = {x ∈ (0, 1) | x ∈ Q}
4. S = {x ∈ R | x = 1/n, n ∈ Z}
With this background in mind, we can proceed with the proof of this
theorem. You will prove it for monotone increasing sequences with the proof
for monotone decreasing sequences being similar.
Proof 47. Let {xn }∞n=1 be a monotone increasing sequence bounded above.
Explain why the set {xn } has a supremum u. Let sup{xn } = u.
Our goal is to show that in fact the sequence {xn } converges to u.
Proof 48. Fix ε > 0 and explain why u − ε is not a supremum of {xn }. Use
this to conclude that there exists N > 0 such that if n > N then xn > u − ε.
Proof 49. Conclude the proof by using the result of proof 48 to show that
|xn − u| < ε for all n > N and thus
lim xn = u.
n→∞
2.6 Subsequences
It is sometimes important to consider only a selection of terms in a sequence,
rather than all of the terms. For example, you might want to look at the se-
quence of even indexed terms or the sequence of odd indexed terms because
doing so makes it easy to determine whether the sequence converges or di-
verges. A new sequence taken from a given sequence is called a subsequence
and is formally defined below.
Explore
Exploration 50. Let’s start with a basic example to illustrate this concept.
Consider the sequence
∞
(−1)n
{xn }∞
n=1 = .
n n=1
Prove
Proof 55. Prove that, if a sequence is convergent, its subsequences are also
convergent, and moreover, they converge to the same value.
and
vN = sup {xn : n > N } . (2.4)
Carefully note how both n and N are used in these definitions. For a given
uN or vN , the number N is fixed while n indicates all of the subscripts greater
than this fixed value. Thus, each uN is an infimum of the set that is composed
of all elements of the sequence {xn } that have an index greater than N. You
might say that each uN is an infimum of the “tail” of the original sequence.
Similarly, with vN except that these are supremums.
These uN and vN values form sequences in their own right. We can ask
whether these new sequences converge and, if so, what are the limits. Addi-
tionally, it is natural to ask how these limits do, or do not, relate to the limit
of the original sequence when one exists.
Definition 2.7 The limit infimum of a sequence {xn } (denoted lim inf) is
Similarly,
Explore
(−1)n
xn = 1 +
n
to answer the following questions.
1. Make a list of the sequence elements x1 through x6 . Plot these points
on a graph with n on the horizontal axis and xn on the vertical axis.
2. Let’s start with the infimums. Compute u1 through u6 where
Again, there is a tool on the website that you should use to verify your think-
ing.
Conjecture
Let’s see if we can turn these observations about specific examples into
formal mathematical statements.
Liminfs and Limsups 27
Conjecture 60.
Conjecture 61.
Conjecture 62.
Suppose that the sequence {xn } converges to L. Then the sequence {uN }
and the sequence {vN } .
Conjecture 63.
Suppose that for a sequence {xn } the sequences {uN } and {vN } converge
to the same value L. Then the sequence {xn } .
Apply
Application 65. Find the lim sup and lim inf of the sequence
nπ
an = nsin( 2 ) .
Application 66. Find the lim sup and lim inf of the sequence
nπ
an = sin .
3
Application 67. Find the lim sup and lim inf of the sequence
1
2 n+1 if n is even
an = .
1 if n is odd
28 Sequences
Prove
The examples that you have worked through in this section have led you
to discover the important relationships between the lim sup and lim inf of a
sequence and how these quantities relate to the limit of the given sequence.
In particular, you conjectured the following theorem.
Lemma 2.1 Given any sequence {xn }, the sequence of infimums {un } is
monotone increasing. Similarly, the sequence of supremums {vn } is monotone
decreasing.
Lemma 2.2 If the sequence {xn } is bounded, then the sequence of infimums
{un } converges. Similarly, the sequence of supremums {vn } converges.
lim xn = L (2.7)
n→∞
Proof 71.
1. Prove that for every ε > 0 there exists N0 such that
Proof 72.
1. Prove that for every ε > 0 there exists N1 such that
and explain why this implies that xn > L − ε for all n > N1 .
Proof 73. How should you choose N so that both proof 71 and 72 are
true? Use this to explain why |L − xn | < ε for all n > N . Finish the proof of
Theorem 2.3.
Explore
lim |xn+1 − xn | = 0.
n→∞
|xm − xn | = ε.
Conjecture
Conjecture 82. Suppose that {xn } is a Cauchy sequence on the open interval
(0, 1). Do you think that it must also converge to a number in (0, 1)? Why or
why not?
Cauchy Sequences 31
Apply
P∞ 1
Application 85. Prove that the series n=1 n2 converges by showing that
the sequence of partial sums is Cauchy.
Application 86. Prove that, if {an } is Cauchy, then the sequence {a2n } is
also Cauchy. Show the converse is not true.
Application 87. Suppose {an } is a Cauchy sequence and that there exists
a subsequence of {an } that converges to L. Prove that {an } converges to L
as well.
Prove
We can now work on showing that if our sequence {xn } is Cauchy then it
converges.
Proof 89. Prove that if {xn } is a Cauchy sequence on R, then {xn } is
bounded. Hint: What is the contrapositive statement?
Proof 90. Fix ε > 0 and use the fact that {xn } is a Cauchy sequence to
show that xm + ε is an upper bound for xn for all n > N.
and
lim inf xn ≤ lim sup xn .
Be careful with the subscripts!
Proof 95. Conclude that lim inf xn = lim sup xn and explain why this implies
that the original sequence converges.
3
Fixed Points and Periodic Points
This equation tells us that fixed points can be found algebraically by solving
the equation f (x) = x for x. We call this special value x̃. This is often a
fairly straightforward process. For example, it is easy to see that the fixed
points of f (x) = x2 are x = 0 and x = 1 since these are the only solutions
to x2 = x. Sometimes, we are unable to compute the exact values of fixed
points. Can you algebraically find the fixed points of cos x? This definition
also tells us something important about the graphical nature of fixed points.
It says that fixed points are the points of intersection of the graphs of y =
f (x) and the diagonal line y = x. You saw this when you were doing the
graphical analysis problems in chapter 1. Figure 3.1 illustrates this with the
two functions discussed above.
Once we know the values of a fixed point, the next thing we need to deter-
mine is its relationship to other orbits of the dynamical system. In particular,
do the orbits of nearby starting values converge to the fixed point, go away
from the fixed point, or exhibit some other behavior less easily categorized?
33
34 Fixed Points and Periodic Points
y y
x x
FIGURE 3.1
Fixed points are the intersection points of the graphs of y = f (x) and y = x.
(Left) The function f (x) = x2 has two fixed points. (Right) The function
f (x) = cos x has a single fixed point.
nearby initial conditions move away from the fixed point under iteration. We
will formally define this concept later. Finally, fixed points that are neither
attracting nor repelling are called neutral.
Explore
Exploration 1. Find the values of the fixed point(s) in terms of the param-
eters m and b. Are there any values that either of these parameters cannot
take?
Conjecture
xn+1 = mxn + b.
Apply
In this section, you will model and analyze the population of walleye in
Teelo Lake where there is a strict, fixed limit on the total number of fish
caught per year. Let n measure the number of years since 2020 (i.e., when
n = 10 it is 2030) and let wn denote the walleye population at year n.
Application 5. Model the walleye population assuming that there is no
fishing whatsoever. Assume that the rate of change of the walleye population
from year n to year n + 1 is proportional to the walleye population at year n
with proportionality constant r > 0. Show that
wn+1 = (1 + r)wn .
Application 6. What are the conditions on r that guarantee that the walleye
population will not go extinct in Teelo Lake? Describe the dynamics of the
walleye population in this scenario.
36 Fixed Points and Periodic Points
Application 7. Suppose that exactly h fish are caught per year. Modify
equation 5 to model this situation. In this harvesting model, what is the
equilibrium walleye population and for what values of r and h (if any) is it
attracting?
Application 11. Discuss the pros and cons of these two models.
Prove
In this section, you will prove the theorem about fixed points of linear
dynamical systems from conjecture 4.
Theorem 3.1 Consider the linear dynamical system
xn+1 = mxn + b.
If m 6= 1 then there is a unique fixed point at x = b/(1 − m). If |m| < 1 then
the fixed point is attracting. If |m| > 1 then the fixed point is repelling.
Proof 12. Use mathematical induction to prove that the solution to equa-
tion 3.2 with b = 0 and initial condition x0 is xn = mn x0 .
Proof 13. Prove that the fixed point of equation 3.2 when b = 0 is attracting
if |m| < 1 and repelling if |m| > 1.
Proof 15. Use the results of proof 13 to prove that the fixed point of equa-
tion 3.2 is attracting if |m| < 1 and repelling if |m| > 1.
Fixed Points 37
Explore
Exploration 21. Think about the explorations above and what you observed
about the stability of the fixed points. Reflect on the relationship between the
graph of the dynamical function and the diagonal line. Is there a computable
quantity that determines whether a fixed point is attracting or repelling?
Fixed Points 39
Conjecture
Conjecture 22. Now let’s put this together as a formal mathematical state-
ment.
Apply
Now seems a good point to discuss why the logistic model in equation 3.5
is a reasonable population dynamics model. Recall from section 1.4 that a
general form of many population models is
Application 24. Let’s turn to the viability function and suppose that
p
V (p) = 1 − .
K
1. What is V (0) and V (K)? What is the interpretation of these prop-
erties?
2. More generally, what can you say about the value of V (p) for 0 ≤
p ≤ K? Why is this important?
3. Is V (p) increasing or decreasing? What is the interpretation of this
property?
4. The parameter K is a population size. What does it signify in this
model?
5. If we assume that K = 1, then how can we interpret the variable
p?
40 Fixed Points and Periodic Points
Application 25. Show that the assumptions of the previous two exercises
give the logistic model of equation 3.5.
Now let’s put together the things we’ve learned about fixed point stability
to analyze the fixed points of this model.
Application 26. For equation 3.5, use conjecture 22 to determine the values
of a where the fixed point at x = 0 is attracting. Your answer should be an
interval. Interpret this in terms of population dynamics.
Application 27. For equation 3.5, use conjecture 22 to determine the values
of a where the fixed point at x = (a − 1)/a is attracting. Your answer should
be an interval. Interpret this in terms of population dynamics.
Application 28. Let’s tie things together by considering the fixed point at
x = 0 of fa (x) = a sin x from exploration 20. Determine the values of a where
this fixed point is attracting.
Prove
The following exercises lead you to one proof of this theorem. There are
other approaches to proving this fundamental result and we encourage you to
think about other methods of proof.
In this proof we assume, without loss of generality, that the fixed point is
x = 0. If it is not, we would simply change variables in a manner similar to
proof 14 in section 3.1.1.
Let’s begin by looking at the case where 0 ≤ f 0 (0) < 1.
Proof 29. Show that there exists ε > 0 such that if 0 < x < ε then
0 ≤ f (x) < x.
Proof 30. Show that if 0 < x0 < ε and xn+1 = f (xn ), then the sequence
{xn } is decreasing and bounded below. Cite the reason why this implies that
the sequence converges.
This definition, while a bit complicated sounding, simply says that a fixed
point is repelling if nearby points iterate away from the fixed point. More
precisely, it says that if we start off at any point x0 in some small neighborhood
U of the fixed point, then at some time (the N th iterate) the orbit of x0 leaves
the interval U. Note that it does not say that the orbit leaves and never comes
back. This can happen. We will see instances of this type of behavior later.
Explore
Your goal in this section is almost identical to what you did in section 3.1.2;
to determine conditions that insure that a fixed point is either attracting or
repelling.
Exploration 34. IBLdynamics.com Review the explorations from sec-
tion 3.1.2 and note when each of the fixed points is repelling. On the website
this is broken up into two tools, one for the logistic function, and one for the
sine function.
Conjecture
Conjecture 35. Reflecting on both the explorations and Theorem 3.2 com-
plete the following to conjecture a theorem on repelling fixed points.
42 Fixed Points and Periodic Points
Apply
Application 36. Return to the logistic equation (3.5) and note the interval
of a-values for which the fixed point at x = 0 is attracting. Now determine
the interval of a-values where the fixed point at x = 0 is repelling.
Application 37. Again for the logistic equation, review the interval of values
of a where the fixed point at x = (a − 1)/a is attracting. Determine the
interval(s) of a-values where it is repelling. Restrict your work to the case
where a > 0.
Application 39. If you need any additional practice on finding fixed points
and determining their stability, here are few more problems. Find all of the
fixed points of the following functions and determine whether each is attracting
or repelling.
1. f (x) = x2 − 1/8
2. f (x) = x2 − 2
3. f (x) = ex − 1
4. f (x) = x3
√
5. f (x) = x
Prove
The following exercises lead you to one proof of this theorem. As in the
attracting fixed point proof, we assume without loss of generality that the
fixed point is x = 0 and begin by considering the case where f 0 (0) > 1.
Proof 40. Show that there exists ε > 0 such that if 0 < x < ε, then x < f (x).
Proof 41. Show that if 0 < xn < ε and xn+1 = f (xn ), then xn+1 > xn .
Proof 42. Show that if 0 < xn < ε, then there exists N > 0 such that
xN > ε.
Proof 43. Repeat proofs 40 through 42, modified appropriately, for −ε <
x < 0 to complete the proof.
Proof 44. Now prove the theorem assuming that f 0 (0) < −1 by modifying
the argument of proof 33.
44 Fixed Points and Periodic Points
As you will begin to see in the explorations below, neutral fixed points are
important when thinking about families of dynamical systems. This is because
when a fixed point is neutral, a small change to the function being iterated
can lead to qualitatively different types of dynamics. This change in behavior
is called a bifurcation and will be the subject of chapter 5.
Our immediate goal is not to explore bifurcations, but simply to categorize
some of the dynamics that occur in the neighborhood of a neutral fixed point.
The dynamical ideas presented here are not a complete list of the properties
of neutral fixed points. Instead, we focus on the most common dynamical
behaviors that can occur. These are easily understood through a combination
of calculus ideas and graphical analysis.
Explore
Exploration 47. Pretend that you could reach into the top two plots ((a)
and (b)) and slide the functions up and down. Describe what happens. Pay
particular attention to fixed points. Something different happens to the lower
two plots ((c) and (d)) when you do this? Describe the difference.
Exploration 48. Pretend that you could reach into the lower two plots and
twist them around x = 0 in one direction or the other (by this I mean you can
push one side up and one side down but the graph is pinned to x = 0). Again,
describe what happens. Pay particular attention to fixed points. Something
different happens to the upper two plots when you do this? Describe the
difference.
Fixed Points 45
y y
y y
x x
FIGURE 3.2
Plots of 4 functions that have neutral fixed points.
It will be helpful if we give some names to the dynamics that you observed
in the explorations above. The easiest is captured in the bottom two figures.
We call a fixed point weakly attracting if it is neutral fixed point and there
exists an interval U containing the fixed point such that initial conditions in
U converge to the fixed point under iteration. Weakly repelling is defined
similarly. There is no standard name for the dynamics in the top two figures.
I’m tempted to combine the words “repelling” and “attracting” and call them
retracting points but that has a different meaning. So let’s use the less clever
term of directionally attracting.
Conjecture
46 Fixed Points and Periodic Points
Keeping in mind the examples above, let’s put together all of the conditions
needed for a fixed point to be directionally attracting or weakly attracting.
Conjecture 49. Let’s start with directionally attracting (graphs (a) and (b)
of figure 3.2).
If
• (the condition that implies that x = 0 is a fixed point)
How does the sign of f 00 (0) determine which side is attracting and which side
is repelling?
Conjecture 50. Now let’s turn to the graphs (c) and (d) in figure 3.2.
If
• (the condition that implies that x = 0 is a fixed point)
• (the condition that implies that x = 0 is neutral) and
• (the condition that implies that x = 0 is attracting)
then x = 0 is a weakly attracting fixed point.
• If then x = 0 is a weakly repelling fixed point.
Apply
Prove
Proof 53. Explain why the last condition of Theorem 3.5 (f 000 (x̃) < 0) is
necessary.
We have focused on neutral fixed points where f 0 (x̃) = 1 in the examples
above but neutral fixed points also occur when f 0 (x̃) = −1. This is why we
needed to solve |fa0 (x)| = 1 in the Apply section above.
48 Fixed Points and Periodic Points
f p (x0 ) = x0 . (3.7)
If p is the smallest integer for this equality to hold, then x0 is a prime period
p point. (Note that the adjective prime here does not modify the number p
and thus we are not saying that the number p is prime! A function might have
a prime period 4 orbit for example.) Additionally, if x0 is a period p point,
then the associated periodic orbit is
{x0 , x1 , . . . , xp−1 , x0 , . . . }
Explore
The following questions explore some of the basic properties of periodic orbits.
Exploration 54. Suppose that x0 is a prime period p point. What are all of
its non-prime periods? Explain.
Exploration 56. Suppose that you know that a dynamical system f has a
prime period 2 orbit and you are given a graph of y = f 2 (x). How could you
locate those points using the graph? How do you make sure that they are the
prime period 2 points? How does this concept generalize to prime period p
points?
Periodic Points 49
Conjecture
Conjecture 59. You know that if a differentiable function f has a fixed point
x0 and |f 0 (x0 )| < 1, then x0 is attracting (and if this value is greater than 1,
it is repelling). Now, suppose that f has a prime period p point x0 . Complete
the following theorem that gives conditions on when a periodic point x0 is
attracting or repelling?
Apply
50 Fixed Points and Periodic Points
Application 60. Let f (x) = x2 − 2. What is f 2 (x) and what degree polyno-
mial is it? What is f 3 (x) and what degree polynomial is it? In general, what
is the degree of f n (x)? What does this tell you about the task of computing
periodic points?
Application 61. The function f (x) = x2 − 1 has both fixed points and
period 2 points. Find them and determine their stability.
Prove
Explore
a-value period
1.75 1
3.25 2
3.46 4
3.63 6
3.84 3
TABLE 3.1
For each of these parameter values, complete the following using the peri-
odic orbit tool on the website.
1. Determine whether there are prime period n orbits for n = 1, . . . , 10.
As n gets bigger, it may be difficult to see (and count) all intersec-
tions of the functions graph with the diagonal. Do the best you
can. For each parameter, make a table that summarizes what other
prime periodic orbits exist for that particular value.
2. Which of these parameter values gave you the most periodic orbits?
What was the period of the attracting orbit for that parameter
value?
Conjecture
Prove
Now we turn to proving some of the things that you observed in the ex-
plorations above. Your main tool for these proofs is Theorem 4.1 and using
the fact that f k (f m (x)) = f k+m (x). You first encountered this fact in explo-
ration 6 of chapter 1.
Proof 66. Prove the following. If f is continuous and f has a prime period 2
orbit, then f has a fixed point. (You’ve already done this! In what problem?)
There are two nice corollaries to the theorem you just proved.
Proof 67. Prove the following. If f is continuous and has a prime period 2n
orbit, then f has a period n orbit. Hint: Let G(x) = f n (x). What is f 2n (x)
in terms of G.
Proof 68. Prove the following. If f is continuous and has a prime period 2n
orbit, then f has periodic orbits of period 2k for all k ≤ n. Hint: How can
you use proof 67?
4
Analysis of Fixed Points
53
54 Analysis of Fixed Points
Explore
Both of the theorems in this section concern functions that map a closed
interval [a, b] into itself. In other words, the output interval of the function is
a subset of the input interval to the function. Here is a simple, yet surprisingly
useful, way of visualizing such function that begins with drawing a rectangle
using the directions below.
• Sketch a set of axes and include on it the diagonal line y = x. On the x-axis,
mark the two endpoints of the interval [a, b].
• Use the line y = x to mark the points (a, a) and (b, b). Also, mark the values
a and b on the y-axis.
• Use the points (a, a) and (b, b) to draw a rectangle. When you are finished,
you should have picture similar to figure 4.1.
Exploration 1. Suppose that f maps the interval [a, b] into the interval
[a, b]. What is the relationship between the graph of f (x) and the rectangle
that you constructed?
Exploration 2. Sketch a function f (x) that maps the interval [a, b] into
the interval [a, b] that has exactly one fixed point. Sketch a function f (x) that
maps the interval [a, b] into the interval [a, b] that has multiple fixed points.
Exploration 3. Sketch a function f (x) that maps the interval [a, b] into the
interval [a, b] that has no fixed points. What property must the function f lack
to make this possible?
Fixed Point Existence Theorems 55
y
FIGURE 4.1
The line y = x and the rectangle created using the two corner points (a, a)
and (b, b).
Conjecture
Apply
on the interval [0, 1] for each k = 1, 2, 3, . . . ? Does this mean that there is no
fixed point in the interval [0, 1]? Plot several of these and discuss what you
see and how that does or does not relate to the theorem.
Does this mean that there is no fixed point in the interval [0, 1]? Plot this
function and discuss what you see and how that does or does not relate to the
theorem.
Does this mean that there is no fixed point in the interval [0, 1]? Plot this
function and discuss what you see and how that does or does not relate to the
theorem.
Prove
Theorem 4.1 If f ([a, b]) ⊂ [a, b] and f is continuous, then there exists a
point c ∈ [a, b] such that f (c) = c.
After drawing the pictures, the proof of this statement probably seems so
obvious that it hardly seems worth doing. However, there is a bit of work to
do.
Proof 9. Find your old calculus book and look up the Intermediate Value
Theorem. Use it to prove this theorem.
Fixed Point Existence Theorems 57
Explore
Exploration 11. How many fixed points must a decreasing function that
maps the interval [a, b] into the interval [a, b] have? Repeat the last two items
of exploration 10 in this case.
Conjecture
Conjecture 12. Thinking about the pictures you drew above, complete the
following statement.
Apply
Prove
Below are the exercises that lead you to the proof of this theorem. But
before completing those steps, it is instructive to think about a “non-proof”
that seems easy, but is unfortunately incorrect.
Proof 14. The following “proof” is wrong for at least two reasons! What
are they and why?
Since f ([a, b]) ⊂ [a, b], we know by the Theorem 4.1 that f has a fixed
point x̃ in the interval [a, b]. By assumption, |f 0 (x̃)| < 1. This implies that
the fixed point is attracting. Thus, if x0 ∈ [a, b], then
Proof 16. Look up the Mean Value Theorem in your calculus book. Explain
why that theorem and the derivative condition given here implies that there
exists λ ∈ (0, 1) such that
Proof 17. Now fix an initial condition x0 ∈ [a, b] that is not the fixed point.
Show that the sequence of iterates {f n (x0 )} is a Cauchy sequence (this requires
using the Mean Value Theorem inequality from proof 16 and the geometric
sum formula). Explain why this sequence must converge to a point x0 ∈ [a, b].
Proof 18. Complete the proof by using the continuity of f to show that x0
is a fixed point of f.
The Inverse and Implicit Function Theorems 59
for all points x and y in our metric space. The contraction mapping the-
orem says that if f is a contraction on a complete metric space, then f
has a unique fixed point. The proof of this much more general theorem is
almost identical to the proof that you completed in the exercises above.
Explore
and solve for (f −1 )0 (y). Remember that since f (x) = y we have f −1 (y) = x.
Now, let’s switch gears and think about the idea of invertibility from linear
algebra.
Exploration 25. If A is an n×n matrix, when is A invertible? There are quite
a few equivalent properties so you should list as many as you can remember.
Exploration 28. Compute DF (1, 1) and DF (0, 0). One of these two matrices
is invertible and one is not. Identify which is which.
Exploration 29. What do you think these calculations imply about the
invertibility of the original function F ?
Prove
Unlike other Prove sections, you will not be lead through a proof of this
theorem. Instead, our goal is to simply present to you the theorem and to
take some time to explain it and relate it to the work you did in the Explore
section. We’ll begin with the version of the inverse function theorem that you
learned in calculus that applies to functions of 1 variable.
Like the first version of this theorem, this is a local existence theorem. It
says that if the Jacobian matrix DF of F is invertible at the point x0 , then
the function F is locally invertible. In exploration 25, you were asked to list a
number of conditions that are equivalent to a matrix being invertible. One of
those conditions is that the determinant should be non-zero. So, this theorem
could have read “if det(DF (x0 )) 6= 0 then . . . ”.
This is helpful in understanding why the second theorem is a generalization
of the first. If n = 1, then DF is a 1 × 1 matrix, which is simply a number, and
the determinant of a number is just that number. We know that numbers are
invertible (i.e., have a reciprocal) as long as they are not equal to 0, which is
exactly the condition given here.
Explore
G(x, y) = x2 + y 2 − 1.
Exploration 30. What is {(x, y) | G(x, y) = 0}? Describe this set both in
words and in a picture.
Exploration 32. Can you use this solution to define a function? How?
√ √
Exploration 33.√Does your
√ answer work for the point ( 2/2, 2/2)? What
about the point ( 2/2, − 2/2)?
Exploration 34. Use the function G to compute dy/dx using implicit differ-
entiation. At what values of (x, y) is dy/dx undefined? Identify these points
on your picture from exploration 30.
The Inverse and Implicit Function Theorems 63
Exploration 35. Can you find a single solution to G(x, y) = 0 that works
for all (x, y) in a neighborhood of the point (1, 0)? If so, what is it and what
is the neighborhood? If not, what went wrong?
Exploration 36. Can you find a single solution √ to √G(x, y) = 0 that works
for all (x, y) in a neighborhood of the point ( 2/2, 2/2)? If so, what is it
and what is the neighborhood? If not, what went wrong?
Prove
The inverse function theorem (4.3) tells us when we can solve y = F (x)
for x. The implicit function theorem is very similar. It tells us when we can
solve G(λ, x) = 0 for x. We will usually be interested in cases where λ is a
real number, but the theorem is the same if λ is a list of real numbers. Thus,
we will state the theorem in that setting.
ψ : λ ∈ Rk | |λ − λ0 | < δ → R
such that ψ(λ0 ) = x0 and G(λ, ψ(λ)) = 0 for all λ with |λ − λ0 | < δ.
Again, note that the implicit function theorem is a local existence theorem.
It does not tell us how to solve for x, simply that a solution does exist in a
neighborhood of a given regular point. The function ψ of the theorem is that
solution.
There are a few more things to note about the implicit function theorem.
First, this is actually a special case of the theorem. There is a version of it
that allows x to be of dimension greater than one similar to the second version
of the inverse function theorem. In this more general version, the ability to
64 Analysis of Fixed Points
Explore
and the line y = x but the exact location of these fixed points may change
depending on the value of parameter a. The following exercises explore what
happens to fixed points as a parameter changes. Use the apps on the website
to answer the questions below.
Exploration 39. IBLdynamics.com Before you use any of the website tools,
look at the graph on the site. Use it to estimate the value of the fixed point of
this function and the value of the derivative of the function at the fixed point.
1. According to definition 4.1, is the fixed point of this function hy-
perbolic or non-hyperbolic?
2. As you interact with the graph, does the fixed point move? Do any
new fixed points appear? Is this new fixed point(s) hyperbolic or
non-hyperbolic?
3. Is this system structurally stable?
4. Can you manipulate the graph in such a way that there are multiple
fixed points? Do you have to manipulate the graph a little or a lot?
Let’s repeat this exercise but with a slightly different starting figure.
Exploration 40. IBLdynamics.com Before you use any of the website tools,
look at the graph on the site. Use it to estimate the value of the fixed point of
this function and the value of the derivative of the function at the fixed point.
1. According to definition 4.1, is the fixed point of this function hy-
perbolic or non-hyperbolic?
2. As you interact with the graph, does the fixed point move? Do any
new fixed points appear? Is this new fixed point(s) hyperbolic or
non-hyperbolic?
3. Is this system structurally stable?
4. Can you manipulate the graph in such a way that there are multiple
fixed points? Do you have to manipulate the graph a little or a lot?
The exercises above should have given you a good intuitive understanding
of the relationship between hyperbolic fixed points and structural stability.
Wiggling the graph a little bit (that is changing the parameter slightly) does
not cause a new fixed point to emerge or the existing one to vanish. However,
we need to be more precise than this and that degree of precision is the main
difficulty in discovering the statement of this theorem. One needs to precisely
state the hypotheses and conclusions using formal mathematical language and
notation. The following exercises are designed to walk you through the process
of expressing this concept precisely in a step-by-step fashion. In the final
exercise, you are prompted to put all of this together in a formal statement
of the Hyperbolic Fixed Point Theorem.
66 Analysis of Fixed Points
Exploration 42. Now suppose that along with the conditions above, Fc (x)
is continuous and differentiable with respect to both c and x. Informally, what
can you say about fixed points of Fc (x) for c 6= c0 but nearby? What can you
say about the derivative at the fixed point for c 6= c0 but nearby?
Exploration 43. How can we make the statement “for c = 6 c0 but nearby”
formal using mathematical language and notation? (Hint: Remember, abso-
lute value means distance!)
Exploration 44. How can we express the statement “there is a fixed point
that depends on the value c” using mathematical notation?
Conjecture
Conjecture 45. Use your answers from above to write the Hyperbolic
Fixed Point Theorem.
Apply
xn+1 = ax(1 − x) − hx
Hyperbolic Periodic Points 67
xn+1 = ax(1 − x) − h.
Prove
Proof 49. Use the Implicit Function Theorem to show that you can solve
G(c, x) = 0 for x in terms of c near c = c0 . Call this function γ(c).
Proof 50. Interpret this result in terms of the original family of functions
Fc .
5
Bifurcations
• There is non-zero fixed point as well. How does it depend on the parameter
a? When is it attracting and when is it repelling?
• For what values of a are there period 2 orbits? How many? Attracting or
repelling?
I think you get the basic idea here. We want to compute the values of the
fixed points and periodic points in terms of the parameter. We also want to
determine the stability of these points in terms of the parameter. Ultimately,
69
70 Bifurcations
our goal is to paint a picture (literally) that will describe the dynamics of fa
for every parameter value a. This picture is called a bifurcation diagram and
you will begin to draw one in section 5.2.
But before we do that, we should first think about not how the parameter
affects the dynamics, but instead how the parameter affects the graph of the
function that determines the dynamics.
Explore
Explore
The following exercises will lead you through the construction of a “fixed
point only” bifurcation diagram for the logistic family with the parameter
a > 0. You may want to refer to your answers from exercises 26 and 27 in
subsection 3.1.2.
Exploration 9. Get a sheet of paper and two colored pens (say red and
blue). The colors will be used to differentiate between attracting (red) and
repelling (blue) fixed points.
1. On the sheet of paper, draw a set of axes with the horizontal axis
labelled a (the parameter) and the vertical axis labelled x (the fixed
point). Mark a = 0, 1, 2, 3, 4 on the horizontal axis.
2. You’ve already shown that one of the fixed points of the logistic
function is x = 0. You also showed that this fixed point is attracting
when 0 < a < 1 and repelling when a > 1. Plot this curve (i.e.,
x = 0) on your axes using the blue pen for the attracting region
and the green for the repelling.
72 Bifurcations
3. You also showed that the other fixed point is given by x = (a − 1)/a
and that it is attracting when 1 < a < 3 and repelling for a > 3.
Plot this curve on the same set of axes using the coloring scheme
above.
Exploration 10. IBLdynamics.com Verify that you got the “fixed point
only” bifurcation diagram for the logistic family that is on the website.
Explore
The name “tangent bifurcation” gives a pretty good hint as to what hap-
pens in this kind of bifurcation. The first exercise illustrates the phenomenon.
Exploration 11. IBLdynamics.com There is a tool on the website to ex-
plore the tangent bifurcation. As you manipulate the graph on the site, pay
particular attention to the appearance of fixed points. Adjust the value of the
parameter c to answer the following questions.
1. Identify the value of c where the graph is tangent to the line y = x.
How many fixed points are there at this value of c? What are their
stabilities?
2. How many fixed points are there prior to this value of c? What are
their stabilities?
3. How many fixed points are there after this value of c? What are
their stabilities?
4. Make a rough sketch of the bifurcation diagram or “fixed point
graph.” Remember that the parameter c is the horizontal axis and
the location of the fixed points is on the vertical axis. When you
did this in exploration 9, you used color to indicate stability.
The Tangent Bifurcation 73
Fc0 (x0 ) = x0 .
Exploration 15. The second condition is that at c0 , the graph of y = Fc0 (x)
must be tangent to the line y = x at the fixed point x = x0 . Express this
condition using formal mathematical notation.
74 Bifurcations
Exploration 17. Finally, the fourth condition implies that when the graph
passes through the diagonal, exactly two fixed points appear. An example of
this not occurring was exploration 13. Express this condition using formal
mathematical notation (Hint: Reviewing subsection 3.1.4 might be helpful).
Take another look at the bifurcation diagram that you created in explo-
ration 11. The curve you drew is not a function of the parameter c since it
fails the vertical line test. But if we rotate this picture ninety degrees, it is a
function. In other words, this “fixed point graph” can be viewed as function
whose inputs are fixed point values (x) and whose outputs are the correspond-
ing parameter value (c). This is different than what we did in the hyperbolic
fixed point theorem. Let’s call this function p(x) so that at the bifurcation
value we have p(x0 ) = c0 .
Exploration 18. Look at the fixed point graph again. What can you say
about p0 (x0 ) and p00 (x0 )?
Conjecture
Conjecture 19. Use your answers from above to write the tangent bifur-
cation theorem.
Apply
Prove
Proof 24. Use the implicit function theorem to conclude that there exists a
smooth function p(x) such that G(p(x), x) = 0 near x = 0. Interpret this in
terms of the original function Fc .
Note that you used the implicit function theorem here to solve for the
parameter c in terms of the variable x. This is just the opposite of what was
done in the proof of the hyperbolic fixed point theorem!
Proof 25. Differentiate G(p(x), x) = 0 with respect to x and solve for p0 (x).
Use this to show that p0 (0) = 0.
Proof 26. Differentiate the expression for p0 (x) again with respect to x. Use
this to show that p00 (0) 6= 0.
Proof 27. Explain what the graph of c = p(x) represents in the bifurcation
diagram. What do the derivative conditions above tell you about the shape of
this graph?
FIGURE 5.1
A plot of the graphs y = f3 (x) = 3x(1 − x) and the line y = x. At the fixed
point, x0 = 2/3, f30 (x0 ) = −1. If we increase the parameter through a = 3, the
fixed point remains though its stability changes from attracting to repelling.
Explore
Exploration 30. IBLdynamics.com The tool on the website for this ex-
ploration looks more closely at the function from exploration 29. Explore how
the fixed point moves and the graph itself looks like it is rotating about the
fixed point. Focus on the rotation. Do you see how the graph of y = fa2 (x)
“rotates” through the line y = x as a passes through the bifurcation value of
a = 3?
1. Explain why this property is important for the creation of three
fixed points of fa2 (x) as a passes through the bifurcation value.
2. What can you say about
∂
fa2 (x)
∂x
for x < 3? For x = 3? For x > 3?
3. What does this exploration imply about the value of
∂2
∂ ∂
fa2 (x) = fa2 (x)
∂a ∂x ∂a∂x
evaluated at a = 3, x = 2/3?
The Period Doubling Bifurcation 79
Conjecture
Conjecture 31. Use your answers from above to write the Period Doubling
Bifurcation Theorem.
xn+1 = Fc (xn )
Apply
Prove
80 Bifurcations
xn+1 = Fc (xn )
∂2
Fc2 (x) |(c0 ,x0 ) 6= 0.
•
∂c∂x
Then there exists intervals I 3 x0 and J 3 c0 and a C 1 function p : I → J
with p(x0 ) = c0 such that
2
• Fp(x) (x) = x, but
• Fp(x) (x) 6= x
for all x ∈ I. Moreover, p0 (x0 ) = 0 and p00 (x0 ) 6= 0.
Instead of proving this theorem, let’s spend some time thinking about what
it says and some of the meanings of the hypotheses and the conclusions.
Proof 34. In part 3 of exploration 29 above, you noted from the graphs
that (Fc20 )00 (x0 ) = 0. Use the hypotheses of the Period Doubling Bifurcation
Theorem to prove that this is true in general.
Proof 35. What do you think the purpose of the third bullet point in the
hypothesis of Theorem 5.2 is? In particular, what do you think it says about
how the graph changes with c?
The function p(x) in the conclusion does come from the implicit function
theorem, but not applied directly. That is basically why we are not proving
this theorem here. The last two questions ask you to think about what the
conclusions of this theorem mean.
Proof 36. What do the two bullet points in the conclusion tell us about the
function p(x)? What does p(x) graph in a bifurcation diagram?
Proof 37. What do the derivative conditions on p imply about this graph?
Can you sketch the bifurcation diagram from this?
6
Examples of Global Dynamics
1. Identify the points whose orbits stay bounded (i.e., don’t go off to
infinity).
2. Determine what periodic orbits do and do not exist. If possible,
describe their stability.
3. Develop a method for describing the kinds of non-periodic orbits
that occur.
In short, we would like to catalog all of the different dynamical behaviors
that occur, even though we may not be able to pinpoint the exact orbits that
exhibit these behaviors.
In this section, we will introduce an incredibly helpful tool for describing
and cataloging the global dynamics of some dynamical systems that exhibit
dynamics that initially seem beyond description. This tool boils down to de-
81
82 Examples of Global Dynamics
scribing orbits with infinite sequences of zero’s and one’s. And these sequences
are pretty straightforward to understand.
But before we do that, let’s look back on one of our earlier examples and
paint a global portrait of its dynamics to get a feel for our goal.
Explore
All of the problems in this section are about the logistic family of dynamical
systems
xn+1 = fa (xn ) = axn (1 − xn ).
Our goal here is to describe the fate of every possible initial condition in
R, but we don’t limit our analysis to the interval [0, 1] as we have done for
much of the earlier work. We are not concerned with the population dynamics
interpretation of this dynamical system. We are simply thinking about it as a
mathematical object, and thus are not limited to only positive values of x.
Application 1. IBLdynamics.com Start with the parameter interval 0 <
a < 1. Create a graph of the logistic function y = fa (x) for a value of a between
0 and 1. You can do this using the tool on the website, using graphing software
of your choice, or by hand. Be sure to include the line y = x in your figure.
You probably want a plotting domain of −1 < x < 2 or slightly bigger. Use
this graph to answer the following questions. Some of them review properties
that you have learned previously. What is different here is that this sequence
of exercises guides you in creating a picture that summarizes global dynamics.
For applications 2 through 8, use the graph from application 1.
Exploration 2. What are the fixed points and what are their stabilities?
Exploration 3. Describe the orbits of all initial conditions that are less than
the smallest of the two fixed points.
Exploration 4. Describe the orbits of all initial conditions that are between
the two fixed points.
Exploration 5. On your graph, locate the two points that are eventually
fixed points. There is one that lands on each of the two fixed points.
Exploration 6. Describe the orbits of all initial conditions that are between
the fixed point x = 0 and the eventually fixed point x = 1.
Exploration 7. Describe the orbits of all initial conditions that are between
the eventually fixed point x = 1 and the other eventually fixed point.
The Logistic Map with a = 4 (Part 1) 83
Exploration 8. Describe the orbits of all initial conditions that are greater
than the largest of the two eventually fixed points.
Application 10. You’ve done it! Color a pair of number lines in some way
to summarize the work that resulted from the work you did in applications 1
through 9. Your colors should signify the different types of behaviors that you
observed.
Explore
Some of the work that you did in applications 1 through 9 will be relevant
here. However, much of this is quite different, very interesting, and somewhat
counterintuitive. Most importantly, the dynamics that you will discover here
and the techniques you will use to discover them illustrate some fundamental
concepts in dynamical systems.
Understanding the properties of the graph of y = f4 (x) is essential because
these ultimately determine the dynamics. Figure 6.1 shows the graph of y =
f4 (x). As usual, the line y = x is also shown. There is also something new in
this figure: the lines x = 1 and y = 1 are superimposed in gray. We’ll see why
shortly.
Exploration 11. Prove that the maximum value of f4 (x) is y = 1. Prove that
the fixed points are x = 0 and x = 3/4 and that both of them are repelling.
84 Examples of Global Dynamics
y
FIGURE 6.1
The graph of the logistic function with a = 4.
Exploration 12. Show, using graphical analysis, that the orbits of initial
conditions x0 with x0 < 0 or x0 > 1 tend to negative infinity.
Remember that our ultimate goal is to describe the orbits of all initial
conditions in R that occur in equation 6.1. At this point, we are only left with
the initial conditions in the unit interval.
Exploration 13. Consider an initial condition x0 ∈ [0, 1]. Prove that the
orbit of this point remains in the interval [0, 1] for all time. More precisely,
prove that if x0 ∈ [0, 1], then f4n (x0 ) ∈ [0, 1] for all n ≥ 0.
Exploration 14. Use the graph in figure 6.1 to identify points that are
eventually fixed to x = 0. Identify points that are eventually fixed to x = 3/4?
Exploration 17. Thinking about the graph of y = f4 (x), sketch the graph
of y = f42 (x) on the interval [0, 1] by hand. Note: The formula of f42 (x) will
not help you very much. Instead, think about how the properties of the graph
of y = f4 (x) determine the graph of y = f42 (x).
The Doubling Map 85
Note: The tool on the website for explorations 18 through 23 will plot the
graphs of y = f4n (x) for n = 1 to 7. It will also perform graphical analysis. Use
it to check your answer to exploration 17 and then to answer explorations 18
through 23.
Exploration 18. IBLdynamics.com How many critical points does y =
f4n (x) have? Is there a general formula for the number of critical points in
terms of n? What are the critical values (i.e., y-values of these critical points)?
Describe the orbits of these points under iteration by f4 .
Exploration 20. IBLdynamics.com How many fixed points does f4n (x)
have? Is there a general formula for the number of fixed points of f4n (x) in
terms of n? Describe the orbits of these points under iteration by f4 .
(as well eventually fixed points), and to the non-periodic dynamics observed
in the stair-step diagrams. The first of these is that the function f4 maps the
interval [0, 1] onto the interval [0, 1]. This implies that every orbit remains in
the interval [0, 1] for all n. The second property is that the function is “two-
to-one.” Every point in the interval [0, 1] (except 1/2) has exactly 2 values
that map onto it. We call these two points preimages. These two properties
combine to cause f4n (x) to have 2n − 1 critical points with critical values of
either 0 or 1, and for there to be 2n period-n points for all n > 0.
Another property of f4 makes computing these critical points and periodic
points difficult, if not impossible. That is its non-linearity. While we can use
the quadratic formula to find the fixed points, this tool goes out the window
very quickly since f4n (x) is a polynomial of degree 2n .
Clearly we can’t construct a continuous, linear, two-to-one function since
the graphs of linear functions are lines! However, there is a simple piecewise
linear function that has these important properties. This function is called the
doubling map and is defined by
2x if 0 ≤ x < 1/2
D(x) = 2x (mod 1) = . (6.2)
2x − 1 if 1/2 ≤ x ≤ 1
The graph of the doubling map is shown in figure 6.2. Iterating the doubling
map is easy; simply multiply by 2 each time and if the number ever becomes
greater than 1, ignore the integer part (that is what the notation mod 1
indicates in equation 6.2). So, for example
5 5 5 2 4 1 2
7→ 7→ ≡ 7→ ≡ 7→ . . .
12 6 3 3 3 3 3
So, 5/12 is an eventually period 2 point as it lands on the period 2 orbit of
{1/3, 2/3}.
Explore
Exploration 26. What are the fixed points of D? What are the prime period
2 points? Period 3 points?
The Doubling Map 87
y
FIGURE 6.2
The graph of the doubling map D(x).
Exploration 28. Now let’s think about the graph of D(x) and its iterates.
1. Carefully plot the graph of y = D2 (x) by hand. How many pieces
are there? What are the points of discontinuity? Identify on your
graph the period 2 points of D and eventually fixed points of D.
2. Repeat for y = D3 (x).
3. Do you see a pattern here? What can you say about the graph of y =
Dn (x). How does this translate to periodic points and eventually
fixed points of D?
Conjecture
Conjecture 30. Complete the following theorem about the number of period
n points of the doubling map D.
88 Examples of Global Dynamics
Conjecture 31. Complete the following theorem about the distance between
consecutive period n points of the doubling map D.
Prove
Theorem 6.1 The doubling map D has 2n periodic points for each n > 0.
Theorem 6.2 The distance between consecutive period n points of the dou-
bling map D is less than or equal to 1/2n−1 for each n > 0.
As we will see shortly, the ability to work with geometric series is extremely
important in understanding the complicated dynamics that we are beginning
to explore. If it has been a while since you worked with geometric series in
calculus, I suggest you take some time to go back and review them. One thing
to note is that it does matter where the summation index (k in this case)
starts and ends. In particular, we will often encounter geometric series that
The Doubling Map 89
where each coefficient dk takes an integer value between 0 and 9. Take for ex-
ample the decimal 0.333 . . . . We know this equals 1/3 but why? Well, because
∞ ∞
X 3 X 1 1/10 1
k
= 3 k
= 3 = . (6.5)
10 10 1 − 1/10 3
k=1 k=1
Explore
Let’s now return to the dynamics of the doubling map, but this time
let’s do it not with traditional fractions or decimals, but with their binary
representation.
Exploration 36. Above, we showed that the binary expansion of 1/3 is
.010101 . . . . Write this out as a series and then compute D(1/3) and D2 (1/3).
Does this agree with the work you did in exploration 25?
Conjecture
Conjecture 38.
Conjecture 39. What are the fixed points of the doubling map D in binary?
The period 2 points? The period 3 points?
Conjecture 40. Give several examples (in binary) of eventually fixed points
of D. Similarly, give examples of eventually period 2 points of D.
The Logistic Map with a > 4 (Part 1) 91
Prove
Proof 42. Suppose that a point s = .b1 b2 b3 . . . in binary. Prove that D(s) =
.b2 b3 . . . in binary.
We say that a binary sequence s = .b1 b2 b3 . . . is a periodic binary se-
quence of period n if it repeats blocks of length n. For example
s = .011011011 . . .
is a period sequence of period 3.
Proof 43. Prove that x is a periodic point of the doubling map D if and
only if the binary expansion of x is a periodic binary sequence.
Explore
All of the questions in this section are about the iteration of the logistic
function when a = 4.5 as shown in figure 6.3. As usual, we’ve included the
graph of the line y = x and again have superimposed the lines x = 1 and
y = 1. Printing out at least one copy of figure 6.3 (it is on the website) may
be helpful in answering the questions in this section.
Exploration 44. Show, using graphical analysis, that the orbits of points
with initial conditions less than 0 or greater than 1 go to negative infinity.
92 Examples of Global Dynamics
FIGURE 6.3
The graph of the logistic function with a = 4.5.
Exploration 45. Now, consider the interval [0, 1]. Graphically, identify the
values of x in this interval where f4.5 (x) > 1. What are the orbits of these
points under iteration?
Exploration 46. Identify the two points xL and xR that are eventually fixed
to x = 0 in exactly two iterations? How do these points relate to your answer
from the previous question?
Exploration 47. Graphically, identify the intervals in [0, 1] that get mapped
into the interval (xL , xR ). What are the orbits of these points under iteration?
How many of these intervals are there, and approximately where are they
located?
All of these things and more are true for the examples that you explored in
chapter 6. Moreover, these properties are universal. They are the fundamental
characteristics of what mathematicians refer to as chaotic dynamical systems.
But before we can prove these facts about the examples of chapter 6, we
need to take a brief detour to introduce two essential concepts: the Cantor set
and the shift map. We begin with the Cantor middle-thirds set that exemplifies
the fractal structure that often accompanies chaotic dynamics. Next, we will
introduce a function called the shift map that captures the essence of the
dynamics that you discovered in chapter 6.
Once we understand these two tools, we can then use them to prove that
the dynamical systems of chapter 6 are in fact chaotic. We begin this process
by using the dynamics to understand the structure of the invariant set and
in many cases, this allows us to prove that this set has the structure of the
Cantor set. We then can use this structure to construct a new function that
puts the dynamics of the given dynamical system into correspondence with the
dynamics of the shift map. In this way, we can prove that a given dynamical
system has the properties listed in the bullet points above. This work will be
done in chapter 8.
93
94 The Tools of Global Dynamics
We imagine continuing in this fashion infinitely many times. The result of this
infinite process is the Cantor middle-thirds set K. Formally, we define
∞
\
K= Kn . (7.1)
n=0
This set has some remarkable properties that you will prove below. Probably,
the most confusing result that you will prove is that in one sense the Cantor
set is “infinitely small,” while in another sense, it is “infinitely large.” The
exercises below will make this more precise and lead you through the required
proofs.
Prove
The first thing we need to do is very basic: we need to show that this
infinite process of removing open middle-thirds intervals leaves something in
the set K.
6 ∅.
Proof 1. Prove that K =
Now that we know that K is non-empty, we aim to do two different things
to get a feel for its size. On the one hand, we want to “count” the number of
elements in K. On the other hand, we want to compute the “length” of K.
The Cantor Set 95
Proof 5. Prove, using a geometric series, that the total length of the intervals
removed to construct K sums to 1. What does this imply about the “length”
of K? Hint: Proof 2 suggests how to set up this series.
You might be a little concerned about this last statement since we know that
1 1 0 0
= + 2 + 3 + ··· ∈ K
3 3 3 3
and this is not of the form given in equation 7.2. This is not the problem that
it seems at first.
96 The Tools of Global Dynamics
Proof 10. Explain why proof 9 implies that there are points in K that are
not endpoints of intervals removed.
Proof 11. Explain why proof 9 implies that K is “really large” in the sense
of counting. Discuss why it is remarkable that K is both “really large” in one
sense and “really small” (proof 5) in another. Compare the properties of the
Cantor set to other, more common, subsets of R.
In other words, Σ2 is a space whose points are all infinite sequences of zeros
and ones. We begin by defining a metric (i.e., a distance function) d on Σ2 .
Let s = .s0 s1 s2 . . . and t = .t0 t1 t2 . . . be any two elements of Σ2 . Define a
metric d : Σ2 → R by
∞
X |sk − tk |
d(s, t) = . (7.4)
2k
k=0
Explore
Exploration 14. Let s1 = .11110 and t = .1. Compute d(s1 , t). Let s2 =
.111110 and compute d(s2 , t). Which of s1 and s2 is nearer to t?
Exploration 16. In general terms, when do you think that two points s and
t are “close together”?
98 The Tools of Global Dynamics
Conjecture
d(s, t) ≤ .
Conjecture 18.
1
Let s, t ∈ Σ2 . If d(s, t) < then .
2n
Apply
Application 19. The ideas presented here can be generalized to more than
just infinite sequences of 0s and 1s. Can you define Σ3 , the sequence space on
3 symbols? What is a metric on this space? What about Σn ?
Prove
Proof 24. Prove that d(s, t), as defined in equation 7.4, is a metric by showing
that it has each of the three listed properties.
σ(.s0 s1 s2 . . . ) = .s1 s2 . . . .
The function σ is called the shift map because it simply shifts the sequence
left one space and deletes the first entry. For example,
σ(.011100 . . . ) = .11100 . . . .
Explore
Exploration 28. What are the period 2 points of σ? Period 3 points? Gen-
eralize for period n points?
100 The Tools of Global Dynamics
Exploration 33. Can you think of any other interesting dynamical behavior
that this function might have? Give an example and explain what happens
and why you think that it is interesting.
Exploration 34. For each n, how many points of period n does σ have?
Note: Count all period n points, not just those of prime period n.
Exploration 36. Suppose that I gave you some point s in Σ2 and some
ε > 0. Describe how to find a periodic point p such that d(s, p) < ε?
Prove
While you are familiar with the idea of a continuous function on R, you
may not have seen that the idea of continuity extends to any metric space
X. The definition of continuity below is almost identical to the definition of
continuity that you first learned in calculus. The only difference is that we
replace the absolute value (the distance function on R) with the more general
distance function d.
Definition 7.1 Let X be a metric space with metric d. A function f : X → X
is continuous if for every ε > 0 there exists δ > 0 such that if d(x, y) < δ,
then d(f (x), f (y)) < ε.
To illustrate this definition let’s begin with proofs involving familiar sets.
Proof 38. Show that the rational numbers are dense in the real numbers.
You should already understand the meaning of the first bullet point in the
definition. The other two need definitions and explanations.
Explore
103
104 Examples of Chaos
There are a lot of quantifiers to parse in this definition, and we should take
some time to get a feel for what it says. In short, it says that F displays sensi-
tive dependence on initial conditions if two close points (x and y), eventually
get separated by a distance b under iteration by F. We need to discuss this
more carefully to fully understand all of the nuances of the definition. The
best way to do this is to slowly and carefully analyze each phrase in the defi-
nition and then see how they work together to capture the concept of sensitive
dependence.
First, definition 8.2 says “for every x,” which means that this phenomenon
happens everywhere in X.
Next, it says “there exists a b > 0.” The constant b is a “separation con-
stant” as it gives us a measure of how far apart the orbits of x and y will
become at later time N . This is stated in the conclusion that there exists an
N such that d(F N (x), F N (y)) > b.
Now on to the ε phrase in the definition. This phrase says “for every ε.”
We should think of ε as a measure of initial closeness. This is expressed in the
inequality d(x, y) < ε.
Finally, let’s talk about the y. This time the definition says “there exists
a y,” which implies that the separation might not happen for all points close
to x, but only for at least one point close to x.
Hopefully, this careful description of the definition is helpful. However, it
is probably best to see what happens when we iterate a system that displays
sensitive dependence on initial conditions. The tool on the website for explo-
rations 1 through 5 allows you to explore this phenomenon when iterating the
logistic map with a = 4.
Exploration 1. IBLdynamics.com Using the tool on the website, set ε =
0.0001. Try to find an initial condition where the distance between iterates
remains less than 1/2 for all iterates displayed on the tool. Repeat this with
ε = 0.00001.
F k (U ) ∩ V 6= ∅.
The first topological transitivity tool on the web site uses the logistic function
f4 (x) to illustrate this basic idea. You define an interval I and the tool displays
its image f4 (I).
Exploration 6. IBLdynamics.com Use the tool on the website to choose
the left and right endpoints (xmin and xmax) of an input interval I so that
both endpoints have values less than 1/2 .
1. How are the endpoints of the image interval f4 (I) determined?
2. For a point in the image interval, how many points in the domain
interval I get mapped onto that point?
Answer these questions again when xmin and xmax are both greater than 1/2 .
Exploration 7. IBLdynamics.com Now, take xmin < 1/2 < xmax and again
use the tool on the website to observe f4 (I).
1. How are endpoints of the image interval f4 (I) determined?
106 Examples of Chaos
2. For a point in the image interval, how many points in the domain
interval I get mapped onto that point?
Prove
F n2 (x) ∈ V.
F n2 −n1 (U ) ∩ V 6= ∅.
Prove
d σ N (s) , σ N (t) ≥ 1.
Proof 20. Devise a different method for constructing another point s0 that
also has a dense orbit.
• use this function to show that the dynamics of the original function are
equivalent to the dynamics of the shift map σ. This equivalence is known as
a conjugacy.
Let S be the function that defines the correspondence between X and Σ2 .
Generally, we need S to create a one-to-one correspondence between points of
our dynamic set X and points of our sequence space Σ2 and thus S must be
• one-to-one,
• onto, and
• continuous with a continuous inverse S −1 .
The method used to establish that two different functions possess equiv-
alent dynamics is summarized in the diagram of figure 8.1 where arrows are
used to represent the action of a function. Along the top row, we represent
the dynamical function F that we are studying. Along the bottom row is a
representation of the shift map σ. The homeomorphism S is represented by
the two vertical arrows on the left and right of the diagram. These connect
the underlying phase spaces. In this case, S maps values of X to values of
Σ2 . Because S is a homeomorphism, we say that the spaces X and Σ2 are
topologically equivalent.
Using the homeomorphism S to establish the topological equivalence of the
spaces is only half of the problem. We also need to establish dynamical equiva-
lence of the functions being iterated. This is done by looking at the interaction
between the dynamical functions and the homeomorphism S as represented
in the commutative diagram. The diagram is about function composition and
the easiest way to understand such a diagram is to follow the functions from
the upper left of the diagram to the bottom right of the diagram.
• If we go “right then down,” the diagram tells us to first apply F (top arrow
right) and then apply S (right arrow down). Thus, that pathway is S(F (x)).
• If we go “down then right,” the diagram tells us to first apply S (left ar-
row down) and then apply σ (bottom arrow right). Thus, that pathway is
σ(S(x)).
When this diagram applies, it says that
FIGURE 8.1
A commutative diagram for the topological conjugacy between a function F
and the shift map σ.
Apply
Diagrams such as the one shown in figure 8.1 are known as commuta-
tive diagrams and are used widely in mathematics where the composition
of multiple functions in multiple ways are being compared. We have labeled
the upper and lower topological spaces as X and Σ2 simply because that is
what they will be in our context. These spaces could be R, Rn , C or almost
anything else.
Application 21. You have encountered the concept embodied in the com-
mutative diagram above in linear algebra when you studied the concept of
similar matrices. A matrix A is similar to a matrix B if there exists an in-
vertible matrix P such that
B = P AP −1 . (8.2)
Application 22. What linear algebra concept does the “commutative dia-
gram” below represent?
Return to The Doubling Map 111
A
~v λ~v
P P
D
~e λ~e
Prove
Proof 23. Review section 6.3 on the doubling map D. What dynamical
properties of D suggest that this function is chaotic on [0, 1]?
section 6.4 where you first explored these dynamics. Figure 8.2 shows the
2 3
graphs of y = f4.5 (x), (left) y = f4.5 (x), (center), and y = f4.5 (x), (right) to
help you remember some of the essential properties. In particular, because the
central portion of the graph of y = f4.5 (x) is above the line y = 1, we know
that the orbits of points in the center iterate to negative infinity. Similarly,
there is an interval to the left of the center and one to the right of the center
that get mapped into the center and hence also iterate to negative infinity. In
fact, you showed that for each n there are 2n−1 open intervals such that if x
is in one of these intervals, then
n−1
2
f4.5 (x) > 1.
Hence, these values of x iterate to negative infinity. This leaves very few points
whose orbits remain in the interval [0, 1] for all time.
FIGURE 8.2
n
The graphs of y = f4.5 (x) with n = 1, 2, 3.
Several steps are required to prove that this function is chaotic and it will
be helpful to outline the process before working through the actual proofs.
1. You have already shown that almost every point has an orbit that
goes to negative infinity. However, there is a set of points (that we
will call Λ) whose iterates remain in the unit interval for all time.
Our first step will be to determine the topological structure of this
set Λ.
2. The second step is to develop a method for describing the orbits of
points in Λ. This method (actually a function) is called an itinerary
because it describes the orbit of a point x by telling us where it is
at every time in its orbit.
3. Finally, we will use the itinerary to create a topological conjugacy
between the logistic function restricted to Λ and the shift map σ on
Σ2 . This will prove that the logistic function is chaotic on Λ.
The Logistic Map with a > 4 (Part 2) 113
Prove
and
1
I1 = x ∈ ,1 f (x) ≤ 1 . (8.4)
2
The intervals I0 and I1 are highlighted along the x-axis in figure 8.3. Note that
I0 ∪ I1 is simply the set of all x-values such that f (x) ∈ [0, 1]. In mathematical
notation
I0 ∪ I1 = {x ∈ [0, 1] | f (x) ∈ [0, 1]} .
FIGURE 8.3
The graph of y = f4.5 (x) with the intervals I0 and I1 highlighted.
Proof 31. Prove that I0 has two disjoint subintervals I00 and I01 such that
• if x ∈ I00 then fa (x) ∈ I0 and fa2 (x) ∈ I0 , and
• if x ∈ I01 then fa (x) ∈ I0 and fa2 (x) ∈ I1 .
Identify these intervals on figure 8.3.
114 Examples of Chaos
Proof 34. In general, how many subintervals of type Is0 s1 s2 ...sk are there if
si = 0 or 1 for i = 0, 1, . . . , k?
Proof 35. What kind of set is being constructed using this process?
Now let’s define the invariant set Λ that we are interested in. Let
Proof 40. Show that Λ contains no intervals by following the steps below.
1. Assume, by way of contradiction, that there exists an interval
[α, β] ∈ Λ. Explain why [α, β] is contained in either I0 or I1 .
2. One implication of the mean value theorem is that if m = min|fa0 (x)|
on [α, β], then
|fa (α) − fa (β)| ≥ m |α − β| .
Use this fact to reach a contradiction.
Proof 41. How can we use the subintervals defined in proof 33 to track the
orbit of a point in x ∈ Λ using a sequence of 0’s and 1’s?
The Logistic Map with a = 4 (Part 2) 115
Proof 43. Assume for the moment that S is a homeomorphism (i.e., one-to-
one, onto, and continuous) to answer the following questions. (You will prove
some of these properties shortly).
1. Prove that Λ is uncountable and conclude that Λ must be topolog-
ically equivalent to the Cantor set.
2. Prove that
S (fa (x)) = σ (S (x))
for all x ∈ Λ. Conclude that fa is chaotic on Λ.
As mentioned above, the itinerary function S is a homeomorphism. Veri-
fying some of the required properties is tricky. In the exercises below, we focus
on two properties that are fairly straightforward to prove.
Proof 44. Prove that S is one-to-one. Hint: Assume that it is not. What
must be true about points in the interval connecting two points with identical
itineraries? Use the mean value theorem again.
Proof 45. Prove that S is onto. Hint: Use the itinerary to create a nested
sequence of closed intervals.
Prove
Proof 46. Let h1 (x) = cos(2πx) and q(x) = 2x2 − 1. Show that
Proof 47. Draw the commutative diagram like that shown in figure 8.1 to
illustrate the result of proof 46.
116 Examples of Chaos
Proof 49. Append this conjugacy to the commutative diagram that you
made in proof 47.
Proof 50. Look at your two-level commutative diagram. What is the conju-
gacy function between the doubling map D and the logistic function f4 ?
Proof 51. Briefly explain why this implies that f4 is chaotic on [0, 1].
Proof 52. Recall that a conjugacy function must be one-to-one, onto, and
continuous. The one that you just created violates one of these properties and
that is why it is called a semi-conjugacy. What property is violated and why?
9
From Fixed Points to Chaos
9.1 Introduction
Let’s take a few minutes to review what we know about the logistic family of
functions. In chapter 3, you showed that the logistic family of functions
• has a fixed point at x = 0 and it is attracting if 0 < a < 1, and
• has a fixed point at x = (a − 1)/a and it is attracting if 1 < a < 3.
You summarized this by creating the plot of the value of the attracting fixed
point as a function of the parameter a like the one shown in figure 9.1.
FIGURE 9.1
A fixed point diagram for the logistic family for 0 < a < 3.
117
118 From Fixed Points to Chaos
FIGURE 9.2
A bifurcation diagram for the logistic family for 0 < a < 4.
This basic outline will need some refinement to get an accurate picture, but
it captures the essence of the algorithm. The next few exercises and theorems
will tell you how to make these refinements to display an accurate bifurcation
diagram.
Explore
Exploration 1. Remember that you only want to see the attracting orbits
(when one exists). Modify the plotting of the points in step 2 of the method
outlined above so that you are only plotting approximations of the attracting
orbits and not the iterates that are “far away” from it?
And you always wondered the third derivative might be good for! An impor-
tant fact related to the iteration of functions is that the Schwarzian derivative
remains negative when we compose a function with itself.
Lemma 9.1 If SF < 0, then SF n < 0 for all n ≥ 0.
The other important concept that we need to introduce is that of a basin
of attraction.
Definition 9.2 The basin of attraction of a fixed point x0 of a function F
is the set of all points whose orbit converges to x0 . The immediate basin of
attraction of x0 is the maximal interval containing x0 that lies in the basin
of attraction.
These two definitions and lemma 9.1 are used to prove the following theorem.
Prove
Proof 5. Prove the following for the logistic family of functions. If fa has
an attracting periodic orbit Ω in [0, 1], then f n (1/2) → Ω as n → ∞. Hints:
The points of a period k orbit are fixed points of what function? How many
critical points must this function have?
Proof 6. Modify the algorithm for computing the bifurcation diagram of the
logistic function using this new idea.
Explore
FIGURE 9.3
The bifurcation diagram for the logistic family for 2.95 < a < 3.6.
Exploration 7. What attracting periodic orbit exists for 3.1 < a < 3.4?
Exploration 9. Using figure 9.3, describe what you think happens to the
periods of the attracting orbits as a increases from 3 to about 3.6.
Exploration 12. Based on these explorations, do you think that this cascade
stops at some finite period 2n or not?
Explore
Apply
Application 19. IBLdynamics.com One more time, but now only between
periods 1 and 6.
1 a b A c d B e 6.
(Note the A and B you determined in application 18.)
Application 20. Do you see a pattern? Could you complete one more step
in this process without actually zooming in?
10
Sarkovskii’s Theorem
10.1 Introduction
One of the overarching goals of this book is to describe, as completely as
possible, the entire range of dynamical behavior that can occur in a given
dynamical system. One aspect of this goal is to describe the types of periodic
behavior that can and cannot occur in a given system. In chapter 6, we did
this primarily by using the shift map to show that those dynamical systems
that are conjugate to the shift map on Σ2 have periodic orbits of all periods
and that those orbits are dense in the invariant, chaotic set defined by the
conjugacy.
The goal of this section is a bit different. We aim to determine what prime
periodic orbits must exist if we already know the prime period of one orbit
already. You have already done a little of this work in chapter 3. In that
chapter, you proved that if f is continuous and there exists a period 2 orbit,
then f must also have a fixed point. But what can we say if somehow we know
that f has a period 13 point? Is there a period 12 point? What about a period
14 point?
The answer to this question turns out to be quite remarkable. And what
makes the answer especially satisfying is that its proof uses no high-powered
mathematical machinery. In fact, it relies almost solely on one often overlooked
theorem of calculus: the intermediate value theorem.
125
126 Sarkovskii’s Theorem
Prove
You will prove the intermediate value theorem by completing each of the
proofs in this section. More precisely, you will prove Corollary 10.1, which is
equivalent to theorem 10.1.
Assume that f (a) ≤ 0 ≤ f (b). Let x0 = a and y0 = b. Define
x0 + y0
m0 = .
2
If f (m0 ) = 0, then x = m0 and we are done. If f (m0 ) < 0, let x1 = m0 and
y1 = y0 . If f (m0 ) > 0, let y1 = m0 and x1 = x0 . Continue in this fashion
creating two sequences {xn } and {yn }.
Proof 1. Explain why if either xn = 0 or yn = 0, then we are done. Now,
assume that this does not happen. Prove that f (xn ) < 0 and f (yn ) > 0 for
all n.
Proof 2. Prove that the sequence {xn } is monotone non-decreasing and that
the sequence {yn } is a monotone non-increasing.
Proof 3. Show that sequences {xn } and {yn } are bounded. What does that
imply about each of them?
Prove
Theorem 10.2 If f is continuous and f ([a, b]) ⊃ [a, b], then there exists a
point c ∈ [a, b] such that f (c) = c.
Proof 8. Prove this theorem using corollary 10.1 of the intermediate value
theorem.
Proof 9. Prove that if f is continuous and f has a prime period 2 orbit, then
f has a fixed point.
The following theorem will be quite useful in upcoming work.
Theorem 10.3 If f is continuous, decreasing, and f ([a, b]) ⊃ [a, b], then
there exists a unique point c ∈ [a, b] such that f (c) = c.
Proof 11. Prove this theorem using corollary 10.2 of the intermediate value
theorem.
Explore
The following exercises look at periodic orbits of the logistic family of func-
tions and serve to motivate Sarkovskii’s theorem. Use the tool on the website
to explore each of the following questions. Keep a chart that summarizes what
prime periodic orbits do and do not exist for each given periodic orbit.
Exploration 12. IBLdynamics.com When a = 2.5, the logistic function has
a fixed point. Determine whether it has a prime period n point for n = 1, ..., 8.
Conjecture
The exploration exercises above point toward a more general truth about
the existence of certain kinds of periodic points and how they may or may not
depend on other types of periodic points. Think about those exercises as you
answer the following questions.
Conjecture 21. If you know that f is continuous and has a fixed point,
what other prime periodic points, if any, must it also have?
Conjecture 22. If you know that f is continuous and has a prime period 2
point, what other prime periodic points, if any, must it also have?
Conjecture 23. If you know that f is continuous and has a prime period 3
point, what other prime periodic points, if any, must it also have?
Conjecture 24. If you know that f is continuous and has a prime period 4
point, what other prime periodic points, if any, must it also have?
Conjecture 25. If you know that f is continuous and has a prime period 5
point, what other prime periodic points, if any, must it also have?
Conjecture 26. If you know that f is continuous and has a prime period 6
point, what other prime periodic points, if any, must it also have?
Conjecture 27. If you know that f is continuous and has a prime period 7
point, what other prime periodic points, if any, must it also have?
Conjecture 28. If you know that f is continuous and has a prime period 8
point, what other prime periodic points, if any, must it also have?
130 Sarkovskii’s Theorem
3 . 5 . 7 . ···
3·2 . 5·2 . 7·2 . ···
3 · 22 . 5 · 22 . 7 · 22 . ··· (10.1)
..
.
. . . 23 . 22 . 2 . 1.
Let’s look at this peculiar ordering more carefully. The ordering starts with all
of the odd numbers greater than or equal to 3 listed in increasing order (the
first line of equation 10.1). Next, come the odd numbers greater than or equal
to 3 multiplied by 2 (the second line of equation 10.1) also listed in increasing
order. Then, the odd numbers greater than or equal to 3 multiplied by 4 (the
third line of equation 10.1). This continues through all positive integer powers
of 2. The ordering concludes with the powers of 2 but now in decreasing order
and thus ending in 1.
Explore
Exploration 29. Prove that every natural number n appears once and only
once in equation 10.1.
Exploration 30. Check your answers to the conjectures 21 through 28. Are
they consistent with this theorem?
Apply
Application 31. At the beginning of this chapter, I posed the question “if
f is continuous on R and has a prime period 13 orbit, must it also have a
prime period 12 orbit or a prime period 14 orbit?” What is the answer to this
question?
Prove
�
� � � � �
�
FIGURE 10.1
(Left) The graphs of a function g(x) that has a period 5 orbit but not a period
3 orbit.
Proof 37. Show that there is not a period 3 point in the interval [1, 2] by
finding g 3 ([1, 2]). This should be done using graphical analysis on the interval
[1, 2].
Proof 38. Repeat exploration 37 for the intervals [2, 3] and [4, 5].
Proof 39. For the function g(x) in figure 10.1, prove that g 3 ([3, 4]) = [1, 5].
Explain why this implies that there is a fixed point of g 3 in [3, 4].
Proof 40. Use the analysis of the interval [3, 4] done in exploration 39 to
show that g 3 ([3, 4]) is decreasing. Use that to conclude that the fixed point of
g 3 in [3, 4] is not of prime period 3. What is it?
As you might imagine from the strangeness of the Sarkovskii ordering, the
complete proof of this statement requires considering many different scenarios.
You will tackle a more manageable piece of the proof. You will prove that the
existence of a prime period 3 orbit implies the existence of prime periodic
orbits of all periods.
Sarkovskii’s Theorem 133
Proof 47. Show that f 4 (A3 ) = I1 . Explain why this implies there is a fixed
point off 4 in A3 ?
Proof 48. Why must this fixed point of f 4 be a prime period 4 point of f ?
Proof 49. Discuss how you might generalize the argument of proofs 44
through 48 to prove that f has a prime period n orbit for all n > 3.
11
Dynamical Systems on the Plane
L(x) = Ax (11.1)
is a linear transformation. In section 11.2, you will begin exploring the dynam-
ics of linear transformations on R2 . In this section, we highlight the properties
of linear transformations that will aid you in that work.
The most important properties of the matrix A in the study of dynamical
systems are its eigenvalues and eigenvectors.
Definition 11.1 If A is an n × n matrix, then λ is an eigenvalue of A with
associated eigenvector v if
Note that equation 11.2 is similar to a fixed point equation where the eigen-
vector v is “almost fixed” by the linear transformation L. We say that it is
almost fixed because the eigenvalue λ is essentially a scaling factor telling us
by how much the eigenvector v gets lengthened or shortened. We will need
to make this description more precise because the eigenvalues may or may
not be real numbers. But for now, it begins to explain why eigenvectors are
important in understanding linear dynamical systems.
Because this is not a linear algebra course, we won’t discuss all of the
important properties of eigenvalues and eigenvectors and what they tell us
about the linear transformation. Below is a short summary of the relevant
facts that we will need here.
135
136 Dynamical Systems on the Plane
det(A − λI) = 0
(A − λi I)vi = 0
P −1 AP = D
7a + 14b = 0
0 = 0.
Note that the matrix equation has reduced to essentially one equation and
two unknowns. Thus, we have a free choice of either a or b. If we choose a = 2,
then b = −1 and thus our eigenvector is v1 = (2, −1)T .
There are two things to note here. First, we are not allowed to choose
a = 0 because that forces b = 0 and and then v1 = 0. This is not an eigen-
vector because eigenvectors must be non-zero. Second, you might have made
138 Dynamical Systems on the Plane
a different choice for a than I did. You might have chosen to let a = 1, which
forces b = −1/2 giving v1 = (1, −1/2)T . But if we multiply your eigenvector
by 2 we get mine. This is what is meant when we say that eigenvectors are
unique up to scalar multiplication.
Exploration 1. Find the eigenvector v2 in the example above.
Explore
In the next few exercises, you will look at several examples that illustrate
the most common types of dynamics that can occur in equation 11.3. You
should use what you have learned about linear systems on R to describe the
dynamics of these planar systems.
Exploration 5. Let
.5 0
A= .
0 .2
1. Rewrite the dynamical system xn+1 = Axn as a system of equations
in terms of the coordinates xn and yn .
2. Describe the orbit of the initial condition (x0 , 0) for all x0 ∈ R.
3. Describe the orbit of the initial condition (0, y0 ) for all y0 ∈ R.
4. Use these results to describe the orbit of the initial condition (x0 , y0 )
for all x0 and y0 ∈ R.
5. Draw a picture in the (x, y) plane that illustrates the dynamics of
this system. Your picture should reflect all of the dynamical prop-
erties that you discovered above.
The pictures that you sketched in the three explorations 5 through 7 are
called phase portraits. Like the graphical analysis figures that you used
to study real dynamical systems, these pictures illustrate the dynamics of
the system. However, because we are now working in two dimensions, we are
no longer able to use the graph of the function being iterated to help us
understand the dynamics. The best we can do is to plot the iterates as points
in the (x, y) plane.
These three examples illustrate three fundamental behaviors that occur in
linear systems. The phase portrait of exploration 5 illustrates an attracting
fixed point sometimes called a sink. Exploration 6 is a repelling fixed point
or a source. Finally, Exploration 7 has one attracting direction (the y-axis)
and one repelling direction (the x-axis). A fixed point such as this is called a
saddle.
Exploration 8. Let
−.5 0
A= .
0 .2
140 Dynamical Systems on the Plane
The only difference between this matrix and that given in exploration 5 is that
the first non-zero entry is negative. How does this change the dynamics? In
particular, is the origin still a sink? What would happen if both entries were
negative? What would happen in the other examples if we changed the signs?
In every example above, the axes are invariant in the sense that if an initial
condition is on one of these lines, then the orbit of that point remains on that
line. This motivates the following definition.
Exploration 9. Let
3/2 1
A= .
0 2
Triangular Matrices
You probably noticed that in the previous two exercises, the eigenval-
ues of A are also the diagonal entries of A. Hopefully, you also realized
that this is true because one of the two off-diagonal entries of the matrix
equals zero. Matrices where all of the entries below (or above) the diag-
onal are equal to zero are called triangular matrices. If a matrix A is
triangular, then the eigenvalues of A are the diagonal entries of A.
If the matrix is not triangular, then the eigenvalues are not the diag-
onal entries!
Conjecture
xn+1 = Axn .
• The lines through the origin defined by the eigenvectors v1 and v2 are
.
• If then the origin is a sink.
• If then the origin is a source.
Apply
3 1
/8 /8
Application 12. A = 1 3
/8 /8
5 7
/6 /6
Application 13. A = 7 17
/12 /12
8 −6
Application 14. A =
3 −1
6 −6
Application 15. A =
1 −1
5 −4
/3 /3
Application 16. A = −2 7
/3 /3
Prove
The following exercises will lead you through a proof of theorem 11.1.
Let’s begin by formally defining a “line through the origin defined by a vector
v.” Intuitively, we understand this to mean the line through the origin that
extends in the direction of v. Precisely, a line `(v) defined by a vector v is the
set
`(v) = (x, y) ∈ R2 | (x, y) = tv for some t ∈ R .
Linear Systems with Real Eigenvalues 143
Proof 18. Suppose that |λi | < 1 and that x0 ∈ `(vi ). Prove that
lim xn = 0.
n→∞
This proves that the eigenvectors associated with real eigenvalues less than
one in absolute value define attracting invariant lines. Now let’s do a similar
proof for eigenvalues greater than one in absolute value.
Proof 19. Suppose that |λi | > 1 and that x0 ∈ `(vi ), x0 6= 0. Prove that
lim |xn | = ∞.
n→∞
Proof 20. Explain why these last three proofs prove theorem 11.1.
There is more that we can say about the orbits of initial conditions on an
invariant line. We can, in fact, easily compute the sequence of iterates of these
points.
Proof 21. Suppose that λi and vi are real eigenvalues and eigenvectors of
A and that x0 ∈ `(vi ). Prove that
xn = λni x0 .
Now let’s turn our attention to initial condition that are not on either
invariant line.
Proof 22. Explain why for all x0 ∈ R2 there exists unique constants C1 and
C2 such that
x0 = C1 v1 + C2 v2 .
Proof 24. Prove if |λi | < 1 for i = 1, 2 then for all initial conditions x0
lim |xn | = 0.
n→∞
In other words, prove that if the origin is a sink, then all orbits converge to
it.
144 Dynamical Systems on the Plane
Proof 25. Suppose that |λi | > 1 and that x0 6= 0. Prove that
lim |xn | = ∞.
n→∞
Proof 26. Why can’t we make a similar statement when |λ1 | < 1 < |λ2 |?
Why do you think this is called a saddle?
Complex or Imaginary?
In my experience, students often confuse the terms complex and imag-
inary. A number z = a + ib is a complex number. The value a is called
the real part of z and the value b is referred to as the imaginary part
of z. If z = ib, then we say that z is pure imaginary. In short, 3 + 4i
is a complex number but 5i is pure imaginary.
Explore
that rotates every point in the plane by some fixed angle φ. How do we choose
the constants in equation 11.5 to accomplish this goal? Given that we want
Linear Systems with Complex Eigenvalues 145
to rotate every point, it seems natural to represent our initial points in polar
form so that
xn = r cos θ
(11.5)
yn = r sin θ
for some fixed values of r and θ. We want the point (xn , yn ) to be rotated by
an angle of φ so that means that we need
xn+1 = r cos(θ + φ)
(11.6)
yn+1 = r sin(θ + φ).
Exploration 27. Look up the sum of angles formulas for sin(θ + φ) and
cos(θ + φ) and rewrite equation 11.6 using them.
Exploration 28. Use equation 11.5 to eliminate both r and θ from your
answer to exploration 27 and replace them with xn and yn .
Exploration 30. IBLdynamics.com Use the tool on the website to plot the
dynamics of equation 11.7 with
√
φ = π/4, φ = π/2, φ = π, φ = π/10, φ = π/ 2
and answer the following questions for each.
1. Do the dynamics that you see depend on the initial condition?
2. Is the fixed point (0, 0) a sink, a source, a saddle, or none of these?
3. Are the dynamics periodic or not? If they are what is the period
and how does this relate to φ?
Linear systems like the ones you’ve just explored are called centers. Points
simply rotate around the origin. These orbits do not tend toward the origin
nor away from the origin as n tends to infinity.
Exploration 31. Show that the eigenvalues of the matrix in equation 11.7
are
λ = cos φ ± i sin φ.
The modulus of a complex number z = α + iβ, denoted |z|, is defined by
p
|z| = α2 + β 2
and generalizes the idea of absolute value on R. To see this note that if z =
α + 0i, then |z| = |α|.
Exploration 32. Show that |λ| = 1 for the eigenvalues of exploration 31.
146 Dynamical Systems on the Plane
Let’s build on what you just learned and consider the dynamical system
xn+1 ρ cos φ −ρ sin φ xn
= (11.8)
yn+1 ρ sin φ ρ cos φ yn
where the parameter ρ > 0. Note that when ρ = 1, equation 11.8 reduces to
equation 11.7.
Exploration 33. IBLdynamics.com Use the tool on the website to draw a
variety of phase portraits for systems of the form given in equation 11.8.
Conjecture
Conjecture 35.
xn+1 = Axn
Apply
7 10
/3 /3
Application 37. A = −4 11
/3 /3
2 4
Application 38. A =
−2 −2
−3 −10
/4 /3
Application 39. A = 1 5
/3 /4
3 10
Application 40. A =
−1 −3
Prove
The work that you’ve done in this section points to an interesting fact:
complex eigenvalues correspond to rotational behavior in a dynamical system.
If λ1 and λ2 are complex conjugate pairs of eigenvalues of a dynamical system,
then they can be written so that
λi = ρ(cos φ ± i sin φ)
modulus ρ is the rate of contraction or expansion and if 0 < ρ < 1, then the
iterates tend toward the origin. We thus call this a spiral sink. If ρ > 1, then
the iterates tend away from the origin and we call this a spiral source. We
call the case where ρ = 1 a center.
The exercises in the apply section above show that this relationship be-
tween complex eigenvalues and the dynamic behavior is true even if the matrix
A is not of the form given in equation 11.8. The proof of this fact is not hard
but relies on some linear algebra ideas that are not discussed above and are
generally not covered in an undergraduate linear algebra course. The basic
idea is that a 2 × 2 real matrix with complex eigenvalues is similar to a ma-
trix of the form given in equation 11.8. This similarity is analogous to the
similarity of a matrix A with real distinct eigenvalues to a diagonal matrix D.
All of the results on linear dynamical systems in R2 discussed in this
chapter are summarized in the theorem below.
Theorem 11.2 Consider the linear dynamical system on R2 defined by
xn+1 = Axn
and let A have eigenvalues λ1 and λ2 .
• If |λi | < 1 for i = 1, 2, then (0, 0) is a sink. Moreover, if the imaginary part
of λi is not equal to 0, then (0, 0) is a spiral sink.
• If |λi | > 1 for i = 1, 2, then (0, 0) is a source. Moreover, if the imaginary
part of λi is not equal to 0, then (0, 0) is a spiral source.
• If |λ1 | < 1 < |λ2 |, then (0, 0) is a saddle.
• If |λi | = 1 for i = 1, 2 and the imaginary part of λi is not equal to 0, then
(0, 0) is a center.
Explore
Conjecture
In the predator-prey model of equation 11.12, you found three fixed points.
Two of them were saddles and the third was a spiral sink. Based on this,
you probably realize that all of the fixed point classifications described in
theorem 11.2 for linear systems also occur in nonlinear systems. This should
not be a surprise. We saw that this is the case for non-linear dynamical systems
on R. We showed that fixed points of linear dynamical systems on R of the
form
xn+1 = ax
are attracting if |a| < 1 and are repelling if |a| > 1. For nonlinear systems of
the form
xn+1 = f (xn )
we have an almost identical condition for a fixed point x̃. Except in that set-
ting, it is the value |f 0 (x̃)| that determines whether x̃ is attracting or repelling.
In short, the derivative evaluated at the fixed point determines stability.
The same thing holds for dynamical systems on R2 . The derivative of
F evaluated at the fixed point usually determines the nature of the fixed
point. But what is the derivative of a function F : R2 → R2 ? You learned in
multivariable calculus that the derivative of a function F : R2 → R2 is a 2 × 2
matrix called the Jacobian matrix. This matrix is denoted DF . The Jacobian
matrix consists of four partial derivatives of F and is given by
!
∂f ∂f
∂x ∂y
DF = ∂g ∂g . (11.13)
∂x ∂y
Apply
Application 48. Find the fixed points and compute their stability.
xn+1 = 2xn + 4yn − 23
yn+1 = −2xn − 2yn + 21
Application 49. Find the fixed points and compute their stability.
xn+1 = x2n + yn
yn+1 = 2xn − 3yn
Prove
Proving that the eigenvalues of the Jacobian matrix DF (x̃) determine the
stability of a fixed point is beyond the scope of this text. Below is a statement
of the theorem that justifies the use of this matrix to classify fixed points.
Periodic Points 153
xn+1 = F (xn )
and suppose that F (x̃) = x̃. Let DF (x̃) have eigenvalues λ1 and λ2 .
• If |λi | < 1 for i = 1, 2, then x̃ is a sink. Moreover, if the imaginary part of
λi is not equal to 0, then x̃ is a spiral sink.
F n (x0 ) = x0 .
The orbit of x0 is the set of points {x0 , x1 , . . . , xn−1 } with F (xi ) = xi+1 and
F (xn−1 ) = x0 .
It is also true that the stability of this orbit is determined by the product
of the derivatives evaluated along the periodic orbit. But here we need to
be careful. For functions on the real line, these derivative values are scalars
and scalar multiplication is commutative. Thus, we don’t need to worry about
the ordering of the multiplication. However, this is not true on R2 . Matrix
multiplication is not commutative! If {x0 , x1 , . . . , xn−1 } is a periodic
orbit, then the stability is determined by the eigenvalues of
in this order.
Apply
As you might imagine, equation 11.14 suggests that it is often quite difficult
to compute the stability of a periodic orbit of a planar dynamical system.
However, the predator-prey model that we have been working with does give
us an opportunity to use this idea in a setting that is fairly straightforward.
So, consider equations 11.12 with a = 3.2 and β = 0.5.
154 Dynamical Systems on the Plane
Application 50. Use what you already know about the logistic equation
to show that there exists a period 2 orbit of the form {(x0 , 0), (x1 , 0)} with
0 < x0 < x1 < 1. Explain why this orbit is attracting when restricted to the
x-axis. (You don’t need to compute the exact values of xi .)
Application 52. Show, in general, that the product of two upper triangular
matrices is upper triangular. Observe that the diagonal entries of the product
matrix are the products of the diagonal entries.
xn+1 = 1 − yn + |xn |
(11.15)
yn+1 = xn .
We won’t be studying this system in any great detail, but it is a fun example
that illustrates some of the basic ideas of this section.
Application 55. Show that the point (1, 1) is a fixed point. Compute the
Jacobian matrix DF (1, 1) and find the eigenvalues λi . Show that |λi | = 1.
Application 56. Show that the point (0, 0) lies on a period 6 orbit. Identify
the points of this orbit in figure 11.1.
Application 57. Show that the point (−1, −1) lies on a period 5 orbit.
Identify the points of this orbit in figure 11.1.
Chaos in the Hénon map 155
��� ����������� ���
�
�
�
-� � � � �
-�
FIGURE 11.1
The Gingerbread attractor.
stable and unstable curves are tangent to the stable and unstable lines Es
and Eu defined by the Jacobian matrix DF evaluated at the fixed point. If
you’ve taken a differential equations course, this is identical to what happens
for saddle equilibrium points in that setting.
However, if we look globally, the stable and unstable curves of a fixed point
can interact in a surprising way. In some systems, the stable and unstable
curves can “cross” each other (this can’t happen in differential equations),
and this leads to a remarkable phenomenon known as a homoclinic tangle.
Homoclinic tangles are the hallmark of chaos in planar dynamical systems.
And although we will not discuss homoclinic tangles in any great detail, the
explorations below point at the complexity of this phenomenon in a dynamical
system known as the Hénon map.
Explore
Exploration 61. Show that the Hénon map has a fixed point at approxi-
mately
(x̃, ỹ) = (0.631355, 0.189405).
Exploration 62. Compute the Jacobian matrix of the Hénon map and
evaluate it at (x̃, ỹ). Compute the eigenvalues and show that this fixed point
is a saddle. Note: This problem should not be done by hand.
Recall that a primary reason that the logistic function is chaotic is that it
is a two-to-one function. A consequence of this is that it is not invertible. In
two dimensions, there is more “real estate” that allows a function to smoothly
stretch and fold a region while still being one-to-one. Each of the next two
exercises shows that the Hénon map is invertible and hence one-to-one. Thus,
something quite different is happening here that causes the apparent chaotic
behavior that you observed in exploration 60.
Exploration 63. Our first method of showing that the Hénon map is in-
vertible uses the inverse function theorem (4.4). Show that if b 6= 0, then the
determinant of the Jacobian matrix of the Hénon map is non-zero for all x
and y. Explain why this implies that the Hénon map is invertible.
Chaos in the Hénon map 157
Exploration 64. Let H(x, y) be the Hénon map of equation 11.16. Show
6 0,
that when b =
y/b
H −1 (x, y) = .
x − 1 + ay 2 /b2
The fact that H is invertible and that we can compute H −1 gives us a
method of approximating the stable curve Ws and the unstable curve Wu .
The following outlines a method to approximate Wu .
1. Make a list of initial conditions near the equilibrium point (x̃, ỹ) on
Eu . We know these points are close to the unstable curve Wu since
Wu is tangent to Eu at (x̃, ỹ).
2. For each of these initial conditions, compute a list of iterates.
3. Plot all of these points.
We approximate Ws in an almost identical fashion. The first change is that
instead of taking initial conditions on Eu , we take them from Es . The second
is that instead of iterating H, we iterate H −1 .
� �
��� ���
��� ���
� �
-� -� � � -� -� � �
-��� -���
-��� -���
FIGURE 11.2
Images showing a homoclinic tangle in the Hénon map. The plot on the left
gives a local picture of Ws and Wu . The plot on the right is a more global
picture.
Figure 11.2 shows the result of this method. The plot on the left is the
“local” picture. The parabolic shaped curve is approximately the stable curve
Ws and the shorter, more horizontal curve is approximately Wu . Note that
Ws and Wu intersect at the fixed point. The figure on the right is the “global
picture” and gives a glimpse of the homoclinic tangle. Again, the parabolic
shaped curve is an approximation of Ws . Note that now Wu looks very much
like the Hénon attractor that you observed in exploration 60. To generate
this figure, we just computed a few more iterates of each initial condition
than we did for the figure on the left. In the “global figure,” we see multiple
intersections of Ws and Wu . This is the appearance of the homoclinic tangle
that creates the attractor and the associated chaotic dynamics.
158 Dynamical Systems on the Plane
One final observation: If we take any small box centered at a saddle point,
each iteration stretches that box (in the unstable direction Eu ) and compresses
that box (in the stable direction Es ). The Hénon map adds an additional
aspect to this transformation - a folding. If you look at an image of the Hénon
attractor you will note that, while not quite symmetric with respect to the
x-axis, it does appear to be folded across the x-axis. This leads to one final
question for you.
Exploration 65. What property of the Hénon map of equation 11.16 might
be responsible for this folding across the x-axis?
12
The Smale Horseshoe
159
160 The Smale Horseshoe
If you’ve ever seen taffy being pulled, these processes (especially the second
one) are essentially a one-dimensional version of taffy pulling. But making
taffy is different in one fundamental way: you can’t put two different parts of
the taffy in exactly the same spot! It is not a two-to-one process. And that
suggests a deep question: can a one-to-one function be chaotic? If the function
f : R → R is continuous, then the answer is no. You should be able to prove
this fairly easily.
But as we mentioned in section 11.6, the plane has more real estate and a
function F : R2 → R2 can twist and turn a planar region in a lot of different
ways without mapping two points onto the same point. In this chapter, we
introduce the prototypical example of a chaotic dynamical system that is one-
to-one, onto, and continuous on an invariant set Λ ⊂ R2 .
Explore
( )
FIGURE 12.1
The Smale horseshoe map F . The domain of the function is the stadium-
shaped region S. The image F (S) is shown in gray.
Exploration 5. Show that F has a unique attracting fixed point in the left
semicircular region of S. What theorem did you use here?
Exploration 9. Don’t try and sketch F 3 (S) unless you are very patient.
How many horizontal strips does it have through R? Generalize to F n (S).
These dynamics should remind of the logistic map with a > 4 that you
studied in section 6.4. In that system, you showed the existence of a closed,
invariant set Λ ⊂ [0, 1] by first showing that there was a single open interval
centered around x = 1/2 that left the interval in one iteration. Then you
showed that there were 2 open intervals that mapped into this center interval
and hence left in two iterations. And so on. You showed that what remained
162 The Smale Horseshoe
was a closed subset Λ ⊂ [0, 1] that is a Cantor set and if x ∈ Λ, then fan (x) ∈ Λ
for all n > 0.
The same idea is happening here except instead of open intervals; there
are open vertical strips that get mapped out of the center rectangle R.
Exploration 10. Draw a figure to explain why there are 3 vertical open
strips U` , Uc , and Ur (for left, center, and right) in the rectangle R such that
F (Ui ) 6⊂ R
for i ∈ {`, c, r}. Explain why there are 2 vertical strips V0 and V1 in R such
that F (Vi ) ⊂ R. Label these on your figure.
Exploration 11. Now explain why there are 3 vertical strips in V0 and 3
more in V1 such that if x is in one of these 6 strips, then F 2 (x) 6∈ R. How
many strips are there whose iterates remain in R for 2 iterations?
Exploration 12. How many vertical strips are there in R such that if x is
in one of these strips then F 3 (x) ∈ R? In general, How many vertical strips
are there in R such that if x is in one of these strips then F n (x) ∈ R?
As in section 6.4, this construction continues for all n leaving a set
Exploration 16. If x ∈ Λ, what can you say about both the forward and
backward orbits of x?
More Symbolic Dynamics 163
Explore
infinite series and then use the geometric series formula to evaluate each of
them.
Exploration 18. Show that equation 12.2 can be rewritten as
∞ ∞
X |s−k − t−k | X |sk − tk |
d(s, t) = + .
2k 2k
k=1 k=0
Exploration 20. What is the maximal distance between any two points in
Σ2 ?
Prove
Proof 21. Prove that d, as defined in equation 12.2, is a metric. The prop-
erties of a metric are given in subsection 7.3.1.
Theorem 7.1 provided us with a natural way to see when two one-sided
sequences s and t are near to each other. It stated that if s and t are two
one-sided sequences such that si = ti for i = 0, 1, . . . , n, then d(s, t) ≤ 21n .
This needs to be modified for two-sided sequences.
Proof 22. Prove the following theorem.
Theorem 12.1 Let s, t ∈ Σ2 and suppose that si = ti for i = −n, . . . , n.
Then
1
d(s, t) ≤ n−1 .
2
But, in our new version of Σ2 , we can define a one-to-one version of the shift
map σ : Σ2 → Σ2 by
Notice that the center (i.e., the “dot”) has moved one unit to the right or
equivalently, the sequence has shifted one unit to the left. Hence, the name
shift map.
Explore
Exploration 25. What are the fixed points of σ? What are the period 2
points of σ? In general, what are the period n points of σ?
Because σ is invertible, the dynamical system defined by it is reversible;
we can iterate both forward and backward. This leads to a refined definition
of an orbit.
O+ (x0 ) = {f n (x0 ) : n ≥ 0} .
and
lim σ n (s) = . . . 00.00 . . . .
n→−∞
Prove
Exploration 38. Figure 12.3 shows a rectangle R and its image H 2 (R) by
the Hénon map. Compare figure 12.3 with the Smale horseshoe figure that
you created in exploration 8. What are the similarities and differences?
A Horseshoe in the Hénon Map 167
H(R)
FIGURE 12.2
A region R (dashed box) and its image H(R) via the Hénon map.
H 2 ( R)
FIGURE 12.3
A region R (dashed box) and its image H 2 (R) via the Hénon map.
13
Generalized Symbolic Dynamics
We have used the term open set earlier in the text when we referred to open
intervals or open sets in Rn . The following definitions extends this concept to
general metric spaces.
Definition 13.2 Given a metric space (M, d), an open ball of radius ,
centered at a point x is the set
169
170 Generalized Symbolic Dynamics
A subset U ⊂ M is open if, for every u ∈ U , there exists some , such that
B (u) ⊂ U .
Additionally, both the metric space M and the empty set ∅ are defined to
be open sets.
If we have an open set, we will also need a definition for a closed set.
Definition 13.3 A subset V ⊂ M is closed if its complement is open in M .
In any metric space, there exists sets that are both open and closed
(unimaginatively called “clopen”). Keep that in mind as you go through this
chapter. Open sets and closed sets aren’t mutually exclusive.
Sometimes, you start with an open set, and you want to make it closed.
Definition 13.4 Given an open set U , its closure U is the smallest closed
set containing U .
The one theorem in this section is commonly proven in a course on real
analysis or topology. If you haven’t seen this proof before, it’s worth giving it
an attempt. The fact that this is the only proof assigned here should give you
some insight into its use later in the chapter.
Prove
σ(.s0 s1 s2 . . . ) = .s1 s2 . . . .
Explore
Exploration 4. Prove that the periodic points of the shift map σ are dense
in Σ3 .
The dot to the left of the 0-th term is the origin. In the case that the
alphabet is the set {0, 1, . . . , n − 1}, the space is denoted the full n-shift.
A word or a block in A is any finite list of symbols chosen from A. The
number of symbols in a block u is its length, denoted |u|. Given s ∈ AZ , we
use s[i,j] to represent the block of s from the i-th to the j-th entry of s.
Again, you have encountered this concept previously in this book. Consider
the point s = . . . 01.01 . . . in the two-sided shift space Σ2 of chapter 12.
This sequence is constructed by recursively concatenating the word u = 01.
Additionally, |u| = 2. Of course, this is not the only word in s. It might be
helpful to list all of the words of length 2 and of length 3 that appear in s.
We will often want to consider shifts where certain words or blocks do
not appear. To do that, we first let F be a set of words constructed from
the alphabet A that we call the set of forbidden blocks. We then define
XF ⊂ AZ to be the set of two-sided sequences s ∈ AZ which do not contain
any of the words in F. The set XF is also a shift space.
The shift map from chapter 7 can be extended to shifts on spaces with
forbidden blocks in the natural way, with σ : AZ → AZ given by
The action of the shift map is the same as before. What can be different is
the underlying shift space. Because of this, it is common to use a subscript on
the function name to indicate the underlying shift space.
When we use the term shift dynamical system, we mean a shift space
X and an associated shift map σX . We write this as an ordered pair (X, σX ).
It is sometimes the case that it is easier to describe a shift space by what is
allowed, rather than what is not allowed (i.e., the set of forbidden blocks). Let
X be a subset of a full-shift, and Bn (X) the set of all n-blocks which appear
in any element of X. The language of X is the set
∞
[
B(X) = Bn (X).
n=0
Shift Dynamical Systems 173
.
Our next step is to construct a useful metric on these shift spaces. The gen-
eralization to arbitrary symbols, combined with a restricted language, makes
our original metric in equation 13.2 problematic. So let’s consider a different
metric, ρ : XF × XF → R defined by
2−k
6 t[−k,k] and s[−i,i] = t[−i,i] for all i < k}
if s[−k,k] =
ρ(s, t) =
0 if s = t
(13.4)
In proof 28, you will prove that ρ satisfies the properties of a metric. For now,
note that like all other metrics that we have used thus far, points are “close”
if they agree near the origin (i.e., the decimal point).
It can be confusing when we try to write sequences of points in shift spaces
since each point is also a sequence. Going forward, we’ll be using a subscript
to indicate the coordinate in a point and a superscript to indicate the index
in a sequence of points. Thus, s4 gives the 4th entry of the sequence s, s4
gives the 4th element of the sequence {sk }, and s44 gives the fourth entry of
the 4th element of the sequence {sk }. This is consistent with how we denoted
elements in a shift space in chapter 7, but inconsistent with the notation in
chapter 2. We apparently can’t make shift spaces without breaking some eggs,
so we’ve made a choice of which egg to break.
For example, if we define the sequence {xn } in the one-sided shift on two
symbols as the sequence of points of all terms 0, except for a 1 in coordinate
n, we would have
x0 =.10000 . . .
x1 =.01000 . . .
x2 =.00100 . . .
..
.
Then, we would use the notation x20 = x21 = 0, and x22 = 1. Remember, we
start indexing the sequence at 0, starting to the right of the origin.
Explore
Exploration 9. Consider the full 2-shift. What is the alphabet? What are
the forbidden blocks? Are you restricted in what symbols you use?
Exploration 11. In the shift space {a, b, c}Z , find the distance ρ between
the given points.
1. . . . abbac.abbac . . . and . . . aabac.abbbc . . .
2. . . . abbac.abbac . . . and . . . aabac.abcbb . . .
3. . . . abbac.abbac . . . and . . . abbac.cbbac . . .
Exploration 13. Consider the shift space {0, 1, 2}Z , and the sequence
x0 = . . . 000.000 . . .
x1 = . . . 001.100 . . .
x2 = . . . 0012.2100 . . .
x3 = . . . 0120.0210 . . .
..
.
List x4 through x6 . Does this sequence converge? Why or why not?
Exploration 15. Does the sequence in {0, 1, 2}Z given below converge? Does
it have a convergent subsequence?
x0 = . . . 1010.0101 . . .
x1 = . . . 2121.1212 . . .
x2 = . . . 0202.2020 . . .
x3 = . . . 1010.0101 . . .
..
.
Exploration 18. Look at the shift space {0, 1}Z . Show that the set
{. . . s−2 s−1 .s0 s1 . . . | s−2 = s−1 = s0 = s1 = 0 and si ∈ {0, 1} otherwise}
is both open and closed.
Exploration 19. Suppose you were to change the values on s−2 through s1
in exploration 18. Pick a few different choices for 0 and 1 for each of these
four terms around the origin. Show that the sets are still always open.
Exploration 20. For any number of the sets you described in exploration 19,
show that their unions are also open.
Conjecture
Conjecture 21.
Suppose you have a shift space on n symbols, and a sequence {xn }. Then
sequence {xn } must/must not have a convergent subsequence.
Conjecture 22.
Conjecture 23. What does it mean for a set to be open in a shift space?
Choose a shift we’ve already discussed, describe an open set, and conjecture
a good definition for describing these sets in general.
To check the validity of your conjecture, verify that unions of open sets
are open using your definition. Similarly, verify that finite intersections of
closed sets are closed using your definition. (These are required properties of
a topological space.)
176 Generalized Symbolic Dynamics
Apply
Application 25. List all words in B1 (X)? List all words in B2 (x)? List all
words in B3 (x)?
Apply
Application 27. In this application, you will prove that a shift space X is
compact. The proof uses a variation of the Cantor diagonalization argument
that is commonly used to prove that the real numbers are uncountable. You
may wish to review that proof before working through the steps outlined
below.
1. Assume you have a shift space X with an alphabet A and some
sequence of points {xn }∞n=0 in X. It will be essential to remember
that a shift space has a finite alphabet.
2. Now, consider the set {xn0 }. This is a sequence of characters in
A, made from the 0-th coordinate of every point in the sequence
{xn }∞
n=0 . Explain why there exists a c0 ∈ A and an infinite set
N0 ⊂ N such that xk0 = c0 for all k ∈ N0 .
3. Next, show that there exists some infinite set N1 ⊂ N0 , and some
c−1 , c1 ∈ A, such that, for all k ∈ N1 , the block between index -1
and 1 satisfies
xk[−1,1] = c−1 .c0 c1 .
Prove
Proof 28. Prove that the function ρ defined in equation 13.4 is a metric.
Proof 29. In chapter 7 you proved that the one-sided shift map on Σ2 is
continuous. Using the metric ρ defined in this chapter, prove that the shift
map σ defined on AZ is continuous.
Proof 30. Let X be a shift space, u an element of the language B(X), and
k ∈ Z. A cylinder set in X is defined as the set
We interpret CkX (u) as the set of all points which the block u occurs
starting at index k. We say this cylinder set is centered at u and of radius
k.
1. Prove that a cylinder set in a shift space is open.
2. Prove that a cylinder set in a shift space is closed.
Definition 13.6 If the set F is finite, the shift space XF is a shift of finite
type (SFT).
This allows us to consider only potential blocks of the same length as those
in the set of forbidden blocks.
Explore
Exploration 32. What are the forbidden blocks in a full shift on n symbols?
Exploration 33. Is the empty set, ∅, a shift of finite type? Is it a shift space?
Exploration 35. Let A = {a, b} and F = {bb}. Consider the 2-step shift of
finite type XF .
1. Find all possible blocks of length 2. Since this is a 2-step shift,
remember that all blocks of length greater than 2 must be made up
of these blocks of length 2.
2. Find all of the blocks of length 3 in this space.
3. Find all of the blocks of length 4 in this space.
4. Could you reconstruct this space, changing the number of steps,
where all blocks in F are of length 3? Of length 4? What do you
think would be the advantage or disadvantage of doing so?
Shift Dynamical Systems 179
Exploration 37. Consider the shift space on A = {0, 1}, where for any point
in the space, each pair of 1’s must be separated by an even number of zeros.
So . . . 100100.001 . . . would be allowed, but . . . 10010.001 . . . would not. Is this
a shift of finite type?
Conjecture
Recall the definition of a topological conjugacy from section 8.3 and con-
sider this commutative diagram in figure 13.1. Assume that X is and Y is
a shift space. F1 is the set of forbidden blocks in X, and F2 is the set of
forbidden blocks in Y . We assume S is a conjugacy.
σX
XF1 XF1
S S
σY
YF2 YF2
FIGURE 13.1
A commutative diagram between shift spaces XF1 and YF2 .
Apply
For this reason, the shift space in exploration 35 is usually called the
Fibonacci shift, and the tiling in application 39 is called the Fibonacci
tiling.
Representing Shift Spaces with Graphs 181
The graph G is shown in figure 13.2. In this example, the initial and terminal
functions of edge e0 are i(e0 ) = 0 and t(e0 ) = 1, respectively.
e0
0 1 e2
e1
FIGURE 13.2
Digraph with adjacency matrix 13.5.
So, what does this have to do with dynamical systems and shift maps? We
have a finite set of edges we’re moving along and we can track which edges
we use, moving forward or backward in time. Thus, tracking the sequence of
edges traversed gives a shift space!
182 Generalized Symbolic Dynamics
The last condition checks that the edges in the sequence actually can follow
each other in the directed graph.
Similarly, we can also track the sequence of vertices visited. This is known
as a vertex shift.
Explore
Exploration 40. Draw the graph for each of the following adjacency matri-
ces. Describe the edge shift for each graph. Describe the vertex shift for each
graph.
1.
1 1
1 1
2.
0 1
2 0
3.
1 2 1
1 0 1
1 0 0
Exploration 41. How is the adjacency matrix given in equation 13.5 and
its associated graph in figure 13.2 related to a period 3 orbit of a continuous
dynamical system on R?
Exploration 42. For each edge shift in exploration 40, can you find either
a language or a set of forbidden blocks? What about for the vertex shift?
Exploration 44. For the graph associated with equation 13.5, find the
eigenvalues of the adjacency matrix. What do you suppose is the relation-
ship between these values, the shift space in exploration 35 and the tiling in
application 39?
Conjecture
Conjecture 46.
Conjecture 47. Take a graph from an earlier exercise, and it’s adjacency
matrix A. We know what information about the graph is available from A.
Calculate A2 , A3 , and A4 . What information about the graph is contained in
Ak ? In particular, what do the diagonal entries compute?
Apply
Application 48. In exploration 41, you made the connection between pe-
riod 3 orbits and the adjacency matrix of equation 13.5. How can you use
conjecture 47 to count the number of period k orbits?
Application 49. The trace of a matrix M is the sum of the diagonal entries
of M . For the adjacency matrix of equation 13.5, compute the trace of Ak for
k = 1, 2, . . . , 5. What is the pattern of these numbers? What does that tell
you about periodic orbits?
184 Generalized Symbolic Dynamics
Note that in a higher edge graph, there is at most one edge between any
pair of vertices. We can use this when creating a shift directly from an adja-
cency matrix.
Let’s look at an example to motivate this definition. Start with the graph
G in figure 13.3. This should look familiar.
e0
0 1 e2
e1
FIGURE 13.3
It has vertex set 0 and 1. The first higher edge graph G[1] has a vertex set
made up of paths of length 0. We’ll say a path of length 0 is just a vertex so
that by definition G[1] = G.
Now, construct G[2] . The paths of length 1 in G define the new vertex set.
In other words,
V G[2] = {e0 , e1 , e2 } .
Definition 13.10 tells us how to construct the edges between these vertices.
From the graph of G, we see that after traversing e0 , we could either go on
edge e1 or edge e2 . In other words,
So, we need an edge from node e0 to node e1 and a second edge from node
e0 to node e2 . We label these edges e0 e1 , and e0 e2 . Similarly, we need edges
e1 e0 , e2 e1 , and e2 e2 . This gives the graph G[2] as shown in figure 13.4.
Representing Shift Spaces with Graphs 185
e2 e2
e2
e0 e2 e2 e1
e0 e1
e0 e1
e1 e0
FIGURE 13.4
Second higher edge graph of the graph G in figure 13.2.
Explore
Exploration 50. Using the example from figures 13.3 and 13.4, construct
G[3] . To get you started, the nodes of this graph are the edges of G[2] and so
there are 5 of them.
Exploration 51. Write a few elements from the shift spaces associated with
G, G[2] , and G[3] . Does taking a higher edge graph cause any forbidden words
to appear in elements of the shift space?
Exploration 52. For each graph in exploration 40, find the alphabet of the
associated edge shift. You can use any labeling you like. Find the alphabets
for the 2nd and 3rd higher block codes for each shift you’ve described.
Exploration 53. For each graph in exploration 40, find the adjacency matrix
of the 2nd higher edge graph.
Exploration 55. Draw a directed graph using at least 3 vertices and at least
6 edges. No vertex should be isolated and you should have at least once edge
leaving each vertex. Self-loops are allowed.
In a graph, a cycle is a path that starts and ends at the same vertex.
1. Write the adjacency matrix A for the graph you’ve drawn. Calculate
the trace of A (the sum of the diagonal entries) and count the
number of self-loops you’ve included.
2. Multiply the matrix by itself to get A2 . For each vertex, count the
number cycles using exactly two edges. Compute the trace of A2 .
3. Try this with a few higher powers on the matrix. For each Ak ,
calculate by hand the number of cycles of length k for each vertex,
then calculate the trace of Ak .
4. How does this relate to periodic orbits?
Conjecture
Conjecture 56.
Conjecture 57.
The trace of the k-th power of the adjacency matrix gives you the number
of in the associated edge shift.
Apply
Application 58. Draw the directed graph G for the adjacency matrix
2 3
.
1 1
Label the vertex on the left v and the vertex on the right w. Label the two
self-loops on v as a and b and the three edges from v to w as c, d, and e. Label
the edge from w to v as f and the self-loop on w as g.
Representing Shift Spaces with Graphs 187
1. Introduce a probability function P (G) → [0, 1]. Let P (v) = 1/3 and
P (w) = 2/3.
2. Repeatedly use the random number generator on the website to
choose a value between 1 and 3. This process generates an element
of the shift space associated with this graph. Think of this as moving
through the graph. If you get a 1, you’re at v, if you get a 2 or 3,
you’re at w.
3. Generate a few elements of the vertex shift, using these probabili-
ties to choose which vertex is next in the sequence. Are there any
forbidden words?
4. We next want to assign probabilities to each edge. We first introduce
the idea of conditional probability. Let P (E|V ) be the probability
of taking edge E given that you are at vertex V .
(a) Which conditional probabilities need to be zero?
(b) Choose probabilities for each edge such that
X
P (E|V ) = 1 for all V ∈ V (G).
E∈E(G)
P (E1 E2 . . . En ) = P (i(E1 ))P (E1 |i(E1 ))P (E2 |i(E2 )) · · · P (En |i(En )).
6. Show that the set of possible sequences are exactly the blocks in
the language for the edge shift XG .
What you have constructed is an example of a Markov Chain.
Prove
Proof 60. Prove that, up to renaming of symbols, the 1-step shifts of finite
type are vertex shifts.
188 Generalized Symbolic Dynamics
Proof 62. Prove that if X is an M -step shift of finite type, there is a graph
G such that X [M ] = XbG and X [M +1] = XG . That is, the M -th higher block
code is the same as the vertex shift of G and the M + 1-st higher block code
is the same as the edge shift of G.
Proof 63. Prove that the i, j entry of the k-th power of an adjacency matrix
gives the number of paths involving k edges when moving from vertex i to
vertex j.
Proof 64. Prove that the trace of the k-th power of the adjacency matrix
gives you the number of points of period k in the associated edge shift.
Definition 13.13 Let LP,f be the set of all allowed words for partition P on
(X, f ). The shift space XP,f with language LP,f is the symbolic dynamical
system corresponding to P.
Markov Partitions 189
Explore
Exploration 65. Let M be the unit circle. We think of this as just the
interval [0, 1), but where any function acting on a point in M has its result
taken mod 1 as we did with the doubling map D.
Define the dynamical system φ : M → M by φ(x) = 10x (mod 1).
i i+1
• Show that subdividing M into the intervals ( 10 , 10 ) for i ∈ 0, 1, . . . , 9 is a
topological partition.
• For a given point in M , show that tracking which set in the partition a point
lands under this map gives an element of a one-sided shift. Where do you
encounter an issue?
• Is this a Markov partition?
• Relate this to the shift you described in exploration 8.
Apply
Application 67. If you have not done so, we recommend you review the
materials in section 11.1 prior to working out this application.
Define the torus T2 = R2 /Z2 . You can think of this as the subset of the
plane [0, 1) × [0, 1) where any linear transformation applied to points in the
set is resolved by taking each coordinate modulo 1.
1. Define a linear transformation on R2 with the matrix
1 1
A= .
1 0
Compute the eigenvalues and eigenvectors of this matrix.
2. Let
x
φ(x, y) = A ,
y
with addition modulo 1. Let π : R2 → T2 be given by
A
R2 R2
π π
φ
T2 T2
L2 (0, 1)
P2 (1, 1) L1
P1
P3
(−1, 0)
(0, 1) (1, 1)
P20
P10
P10
P30 P30
P20
(0, 0) (1, 0)
Application 68. Perform the same process as above, using the matrix
2 1
.
1 0
192 Generalized Symbolic Dynamics
Construct the Markov partition and the symbolic representation of the map
on the torus. Show that this map is chaotic.
This is more commonly known as Arnold’s cat map and is an example
of a hyperbolic toral automorphism.
Bibliography
[1] Robert Devaney. A first course in chaotic dynamical systems: Theory and
Experiment. Perseus Books, 1992.
193
Index
195
196 Index