Period Three Implies Chaos Tien-Yien Li; James A. Yorke The American Mathematical Monthly, Vol. 82, No.

10. (Dec., 1975), pp. 985-992.

Stable URL: http://links.jstor.org/sici?sici=0002-9890%28197512%2982%3A10%3C985%3APTIC%3E2.0.CO%3B2-H The American Mathematical Monthly is currently published by Mathematical Association of America.

Your use of the JSTOR archive indicates your acceptance of JSTOR's Terms and Conditions of Use, available at http://www.jstor.org/about/terms.html. JSTOR's Terms and Conditions of Use provides, in part, that unless you have obtained prior permission, you may not download an entire issue of a journal or multiple copies of articles, and you may use content in the JSTOR archive only for your personal, non-commercial use. Please contact the publisher regarding any further use of this work. Publisher contact information may be obtained at http://www.jstor.org/journals/maa.html. Each copy of any part of a JSTOR transmission must contain the same copyright notice that appears on the screen or printed page of such transmission.

JSTOR is an independent not-for-profit organization dedicated to and preserving a digital archive of scholarly journals. For more information regarding JSTOR, please contact [email protected].

http://www.jstor.org Fri Apr 13 11:04:22 2007

PERIOD THREE IMPLIES CHAOS

TIEN-YIEN LI AND JAMES A. Y O R K E

1. Introduction. The way phenomena or processes evolve or change in time is often described by differential equations o r difference equations. One of the simplest mathematical situations occurs when the phenomenon can be described by a single number as, for example, when the number of children susceptible to some disease at the beginning of a school year can be estimated purely as a function of the number for the previous year. That is, when the number X,+I at the beginning of the n + 1st year (or time period) can be written

where F maps an interval J into itself. Of course such a model for the year by year progress of the disease would be very simplistic and would contain only a shadow of the more complicated phenomena. For other phenomena this model might be more accurate. This equation has been used successfully to model the distribution of points of impact on a spinning bit for oil well drilling, as mentioned in [8, 111, knowing this distribution is helpful in predicting uneven wear of the bit. For another example, if a population of insects has discrete generations, the size of the n + 1st generation will be a function of the nth. A reasonable model would then be a generalized logistic equation

A related model for insect populations was discussed by Utidain [lo]. See also Oster et a1 [14,15]

,
FIG. 1 . For K = 1, r = 3.9, with x , = .5, the above graph is obtained by iterating Eq. (1.2) 19 times. At right the 20 values are repeated in summary. No value occurs twice. While x, = .975 and x,, = ,973 are close together, the behavior is not periodic with period 8 since x,, = .222.

These models are highly simplified, yet even this apparently simple equation (1.2) may have surprisingly complicated dynamic behavior. See Figure I . We approach these equations with the viewpoint that irregularities and chaotic oscillations of complicated phenomena may sometimes be understood in terms of the simple model, even if that model is not sufficiently sophisticated to allow accurate numerical predictions. Lorenz [1-4] took this point of view in studying turbulent behavior in a fascinating series of papers. H e showed that a certain complicated fluid flow could be modelled

986

TIEN-YIEN LI A N D J. A. YORKE

[December

by such a sequence x, F ( x ) , F Z ( x ) ,. . . , which retained some of the chaotic aspects of the original flow. See Figure 2. In this paper we analyze a situation in which the sequence { F n ( x ) }s non-periodic i and might be called "chaotic." Theorem 1 shows that chaotic behavior for (1.1) will result in any situation in which a "population" of size x can grow for two or more successive generations and then having reached an unsustainable height, a population bust follows to the level x or below.

FIG 2. Lorenz [ I ] studied the equations for a rotating water-filled vessel which is circularly symmetric about its vertical axis. The vessel is heated near the rim and cooled near its center. When the vessel is annular in shape and the rotation rate high, waves develop and alter their shape irregularly. From a simplified set of equations solved numerically, Lorenz let X,, be in essence the maximum kinetic energy of successive waves. Plotting X.,, against X,,, and connecting the points, the above graph is obtained.

In section 3 we give a well-known simple condition which guarantees that a periodic point is stable and then in section 4 we quote a result applicable when F is like the one in Figure 2. It implies that there is an interval J , C J such that for almost every x E J, the set of limit points of the sequence { F n ( x ) } s J,. i A number of questions remain unanswered. For example, is the closure of the periodic points an interval or at least a finite union of intervals? Other questions are mentioned later. Added in proof. May has recently discovered other strong properties of these maps in his independent study of how the behavior changes as a parameter is varied [17].
2. The main theorem. Let F : J + J. For x E J, F O ( x ) enotes x and F n " ( x ) denotes F ( F n ( x ) ) d for n = 0 , 1 , . . . . We will say p is a periodic point with period n if p E J and p = F n ( p ) and p# F k ( p )for 1 5 k < n. We say p is periodic or is a periodicpoint if p is periodic for some n 2 1. We say q is eventually periodic if for some positive integer m, p = F m (q) is periodic. Since F need not be one-to-one, there may be points which are eventually periodic but are not periodic. Our objective is to understand the situations in which iterates of a point are very irregular. A special case of our main result says that if there is a periodic point with period 3, then for each integer n = 1,2,3,. . . , there is a periodic point with period n . Furthermore, there is an uncountable subset of points x in J which are not even "asymptotically periodic."

19751

PERIOD THREE IMPLIES CHAOS

THEOREM Let J be an interval and let F : J J be continuous. Assume there is a point a E J for 1. which the points b = F ( a ) , c = F 2 ( a ) and d = F 3 ( a ) ,satisfy

987

Then TI: for every k

=

1,2,. . . there is a periodic point in J having period k.

Furthermore, T2: there is an uncountable set S C J (containing no periodic points), which satisfies the following conditions:

(A) For evdry p, q E S with p # q, limsup I F " ( p ) - F n ( q ) i > 0

n-=

and

( B ) For every p E S and periodic point q E J,

limsup I F n ( ~ ) F n ( q ) l > 0. n -=

REMARKS Notice that if there is a periodic point with period 3, then the hypothesis of the theorem will be :]atisfied. An example of a function satisfying the hypotheses of the theorem is F ( x ) = rx[l - x l K ] as in (1.2) for r E (3.84,4] with J = [0, K ] and for r > 4, F ( x ) = max{O, r x [ l - x / K ] }with J = [0, K ] . See [ 2 ]for a detailed description of iterates of this function for r E [O,4). The case r = 4 is discussed in [6, 7, 121. While the existence of a point of period 3 implies the existence of one of period 5, the converse is false. (See Appendix 1). We say x E J is asymptotically periodic if .there is a periodic point p for which

It follows from ( B ) that the set S contains no asymptotically periodic points. We remark that it is unknown what the infimum of r is for which the equation (1.2) has points which are not asymptotically periodic.

Proof of Theorem 1. The proof of T1 introduces the main ideas for both T1 and T2. We now give the prodf of TI with necessary lemmas and relegate the tedious proof of T 2 to Appendix 2.
LEMMA Let G : I + R be continuous, where I is an interval. For any compact interval I , C G ( I ) 0. there is a compact interval Q C I such that G ( Q )= I,.

Proof. Let I , = [ G ( p ) ,G ( q ) ] ,where p, q E I. If p < q, let r be the last point of [p,q ] where G ( r )= G ( p ) and let s be the i'rst point after r where G ( s )= G ( q ) . Then G ( [ r ,s ] )= I,. Similar reasoning applies when p > q.
LEMMA Let F : J + J be co,ztinuous and let { I , , } ~ = ,be a sequence of compact intervals with 1. I,, C J and I,,,, C F(I,,) for all n. Then there is a sequence of compact intervals Q, such that Qn+,C Qn C loand F n ( Q , ) = I,, .for n 2 9. For any x E Q = f l Q, we have F n ( x )E I, for all n.

988

TIEN-YIEN L I A N D J. A. YORKE

[December

Proof. Define Qo = I,. Then FO(Qo) lo. If Qn-]has been defined so that F n - l ( Q n - l = I,-,, then = ) I,, C F(I,-]) = Fn(Q,-,). By Lemma 0 applied to G = F n on Q,-] there is a compact interval Q, c Q,-] such that F n ( Q , ) = I,,. This completes the induction. The technique of studying how certain sequences of sets are mapped into or onto each other is often used in studying dynamical systems. For instance, Smale uses this method in his famous "horseshoe example" in which he shows how a homeomorphism on the plane can have infinitely many periodic points [13].

L E M M2. Let G : J+ R be continuous. Let I C J be a compact interval. Assume I C G ( I ) . Then A there is a point p E I such that G ( p )= p.
Proof. Let I = [Po,P I ] .Choose a i ( i = 0 , l ) in I such that G ( a i )= Pi. It follows a,- G ( a o )2 0 and a l - G ( a l )5 0 and so continuity implies G ( P )- P must be 0 for some P in I. Assume d 5 a < b < c as in the theorem. The proof for the case d 2 a > b > c is similar and so is omitted. Write K = [a, b ] and L = [b, c ] . Proof of T I : Let k be a positive integer. For k > 1 let {I,} be the sequence of intervals I,, = L for - 2 and I,-, = K, and define I,, to be periodic inductively, I,,,,, = I, for n = 0 , 1 , 2 , . . . . If = 1, let I, = L for all n. Let Q, be the sets in the proof of Lemma 1. Then notice that Qk C Q0 and F k ( Q k ) Q0 and so by = Lemma 2, G = F k has a fixed point pk in Q,. It is clear that p, cannot have period less than k for F ; otherwise we would need to have Fk-'(p,) = b, contrary to F k + l ( p kE L . The point p, is a periodic ) point of period k for F. n k
= 0 , . . ., k

3. Behavior near a periodic point. For some functions F, the asymptotic behavior of iterates of a point can be understood simply by studying the periodic points. For

a detailed discussion of the points of period 1 and 2 may be found in [I] for a E [0,4] and we now summarize some of those results. For a E [O, 41, F : [O, 11- [O, 11. For a E [O, I], x = 0 is the only point of period 1 ; in fact, for x E [O, I], the sequence Fn(x)+O as n +x. For a E (1, 31, there are two points of period 1, namely 0 and 1 - a - ' , and for x E (0, I), F n ( x ) + l - a - I as n - x . For a > 3 there are also two points of period 2 which we may call p and q and of course F ( p ) = q and F ( q ) = p. For a E ( ? , I + V% = 3.449) and x E (0, I), F Z n( x ) converges to either p or q while F2""(x) converges to the other, except for those x for which there is an n for which F n ( x )equals the point 1 - a-I of period 1. There are only a countable number of such points so that the behavior c of { F n ( x ) } an be understood by studying the periodic points. For a > 1 + V%,there are 4 poilits of period 4 and for a slightly greater than 1 + V%, F 4 " ( x )tends to one of these 4 unless for some n, F n ( x )equals one of the points of period 1 or 2. Therefore we may sfimmarize this situation by saying that each point in [O, 11 is asymptotically periodic. For those values of a for which each point is asymptotically pefiodic, it is sufficient to study only the periodic points and their "stability properties." For any function F a point y E J with period k is said to be asymptotically stable if for some interval I = ( y - 6, y + 6 ) we have

IFk(x)-yl<Ix-yl

for all

xEI.

If F is differentiable at the points y, F ( y ) ,. . .,F k - ' ( y ) ,there is a simple condition that will guarantee this behavior, namely

PERIOD THREE IMPLIES CHAOS

By the chain rule

where y, is the nth iterate, F n ( y ) .Therefore y is asymptotically stable if

This condition of course guarantees nothing about the limiting behavior of points which do not start "near" the periodic point or one of its iterates. The function in Figure 2 which was studied by Lorenz has the opposite behavior, namely, where the derivative exists we have

For such a function every periodic point is "unstable" since for x near a periodic point y of period k, the kth iterate F k ( x )is further from y than x is. T o see this, approximate F k ( x )by

. J Thus for x near y , ) F k ( x ) - y I is approximately Ix - y 1 I ( d l d x ) F k ( y ) l F rom (3.2) I ( d l d x ) F k ( y ) is greater than 1. Therefore F n ( x )is further from y than x is. We do not know when values of a begin to occur for which F in (3.1)has points which are not asymptotically periodic. For a = 3.627, F has a periodic point (which is asymptotically stable) of period 6 (approx. x = .498).This x is therefore a point of period 3 for F Zand so Theorem 1 may be applied to F Z .Since F 2 has points which are not asymptotically periodic, the same is true of F. In order to contrast the situations in this section with other possible situations discussed in the next section, we define the limit set of a point x. The point y is a limit point of a sequence {x,,} c J if there is a subsequence x., converging to y. The limit set L ( x ) is defined to be the set of limit points of { F n ( x ) } I.f x is asymptotically periodic, then L ( x ) is the set { y , F ( y ) ,. . ., ~ ~ ' ( for) some periodic y } point y of period k.
4. Statistical properties of { F a ( x ) } .Theorem 1 establishes the irregularity of the behavior of iterates of points. What is also needed is a description of the regular behavior of the sequence { F n ( x ) w hen F is piecewise continuously differentiable (as is Lorenz's function in Figure 2) and }

>1

where

Jl =

One approach to describing the asymptotic behavior for such functions is to describe L ( x ) , if possible. A second approach, which turns out to be related, is to examine the average behavior of { F n ( x ) } .T he fraction of the iterates {x, . . ., F N - ' ( x ) }of x that are in [ a l ,a z ] will be denoted by 4 ( x , N, a l ,a z ] ) .The limiting fraction will be denoted [

when the limit exists. The subject of ergodic theory, which studies transformation on general spaces,

990

TIEN-YIEN LI A N D J. A. YORKE

[December

motivates the following definition. We say g is the density of x (for F ) if the limiting fraction satisfies

d ( x , [ a , ,a,]) =

1
0I

O2

g ( x ) d x for all a , , a, E J ; a , < a,.

The techniques for the study of densities use non-elementary techniques of measure theory and functional analysis, so that we shall only summarize the results. But their value lies in the fact that for certain F almost all x E J have the same density. Until recently the existence of such densities had not been proved, except for the simplest of functions F. The following result has recently been proved: THEOREM [5]. Let F : J - t J satisfy the following
conditions: 2. 1 ) F is continuous.
2) Except at one point t E J, F is twice continuously diferentiable.
3) F satisfies (4.1).
Then there exists a function g : J--. [O,m),such that for almost all x E J, g is the density o f x. Also for almost all x E J, L ( x ) = { y : g ( y ) > 0 ) which is an interval. Moreover, the set J, = { y : g ( y ) > 0 ) is an interval, and L ( x ) = J, for almost all x. The proof makes use of results in [8]. The problem of computationally finding the density is solved in [9]. A detailed discussion of (3.1) is given.in [16], describing how L ( x ) varies as the parameter a in (3.1) varies between 3.0 and 4.0. A major question left unsolved is whether (for some nice class of functions F ) the existence of a stable periodic point implies that almost every point is asymptotically periodic.
Appendix 1: Period 5 does not imply period 3. In this Appendix we give an example which has a fixed point of period 5 but no fixed point of period 3. Let F : [ l , 5 ]-t [ l ,51, be defined such that F(1) = 3, F(2) = 5, F(3) = 4, F(4) = 2, F(5) = 1 and on each interval [n,n + 11, I 5 n 5 4, assume F is linear. Then

51) Hence, F 3 has no fixed points in [ l , 21. Similarly, F3([2,31) = [3,5]and F3([4, = [ l ,41, so neither of these intervals contains a fixed point of F3. On the other hand,

Hence, F 3 must have a fixed point in [3,4].We shall now demonstrate that the fixed point of F 3 is unique and is also a fixed point of F. Let p E [3,4]be a fixed point of F3. Then F ( p ) E [2,4].If F ( p ) E [ 2 , 3 ] ,t hen F 3 ( p )would be in [ I ,21 which is impossible since then p could not be a fixed point. Hence F ( p ) E [3,4] and F Z ( p ) [ 2 , 4 ] .If F 2 ( p )E [2,3]we would have F 3 ( p )E [ 4 , 5 ] , n impossibility. Hence p, F ( p ) , F 2 ( p ) a are all in [3,4].On the interval [ 3 , 4 ] ,F is defined linearly and so F ( x ) = 10 - 2x. It has a fixed point 1013 and it is easy to see that F 3 has a unique fixed point, which must be 1013. Hence there is no point of period 3.
Appendix 2. Proof of T2 of Theorem 1. Let with

be the set of sequences M

= {M,},=,

of intervals

(A.1)

M,, = K

or M, C L , and if

F(Mn)>Mn+]

M,, = K then

('4.2)

n is the square of an integer and M,,,,, M,,,, C L,

19751

PERIOD THREE IMPLIES CHAOS

99 1

where K = [a,b ] and L = [b, c ] .Of course if n is the square of an integer, then n + 1 and n t 2 are not, so the last requirement in (A.2) is redundant. For M E .Id, let P(M, n ) denote the number of i's in { I , . . ., n } for which M, = K. For each r E (314,l) choose M r = {ML}:=, to be a sequence in ,d such that
n -=

iim P(Mr,n 2 ) / n= r.

is uncountable since M'I # Mr. for r , f r2. For each Let Mo = { M r :r E (314, 1 ) ) C M. Then M r E .llo, by Lemma 1, there exists a point x, with F n ( x , )E MA for all n. Let S = {x,: r E (314, 1)). Then S is also uncountable. For x E S, let P(x, n ) denote the number of i's in ( 1 , . . ., n } for which F 1 ( x ) K. We can never have F k ( x , )= b, because then x, would eventually have period 3, contrary E to (A.2). Consequently P(x,, n ) = P ( M r ,n ) for all n, and so

p(x,) = n - x P(X,, n 2 )= r lim

for all r. We claim that

(A.4) for p, q E S, with p# q, there exist infinitely many n ' s such that F n ( p )E K and F n ( q )E L or vice versa.
We may assume p ( p ) > p(q). Then P(p, n ) - P(q, n)+'m, and so there must be infintely many n's such that F n ( p )E K and F n ( q ) EL. Since F 2 ( b )= d 5 a and F z is continuous, there exists 6 > 0 such that F Z ( x ) ( b + d ) / 2 for all < x E [ b - 6, b ] C K. If p E S and F n ( p )E K, then (A.2) implies F n + ' ( p ) L and Fn+*(p) L. E E Therefore F n ( p )< b - 6. If F n ( q )E L, then F" ( q )2 b so

By claim (A.4), for any p, q E S, p# q, it follows lim sup I F " ( P ) - F n ( q )1 2 6 > 0.

n-=

Hence (2.1) is proved. This technique may be similarly used to prove (B) is satisfied.

Proof o f 2.2. Since F ( b ) = c, F ( c ) = d 5 a, we may choose intervals [ b n ,e n ] ,n = 0,1,2,. . . , such that (a) [ b , c ] = [ b O , c O ] > [ b l , c ' ] > ~ ~ ~ > [ b " , c " ] > ~ ~ . , (b) F ( x ) E ( b ", c ") for all x E ( b " + I , c "+I), (c) F ( b n + ' ) e n , F ( c n + ' )= b". = Let A = r ) ; = , [ b n , c n ] ,b* = infA and c * = sup A, then F ( b * ) = c* and F ( c * ) = b*, because of (c). In order to prove (2.2) we must be more specific in our choice of the sequences M r . In addition to our previous requirements on M E M, we will assume that if Mk = K for both k = n 2 and ( n t I)* then Mk -[ b z n ( 2 1 - ' ) , b *for k = n 2 + ( 2 j - l ) , Mk = [ c * , c 2 " - 2 ' ] k = n 2 + 2 j where j = I ; . . , n . ] for For the remaining k's which are not squares of integers, we assume Mk = L. It is easy to check that these requirements are consistent with ( A . l ) and (A.2), and that we can still choose M r so as to satisfy (A.3). From the fact that p ( x ) may be thought of as the limit of the fraction of n's for which ~ " ' ( xE K, it follows that for any r*, r E (314,l) there exist infinitely many ) n such that M ; = M;' = K for both k = n z and ( n t 1)'. To show (2.2),let x, E S and x,. E S. Since b n + b*, c" + c* as n + x , for any E > O there exists N with ( b n- b * ( < ~ / 2( ,e n- c * ( < ~ / for 2 Then, for any n with n > N a n d M ; = ML* = K for both k = n 2and ( n t I)', we have all n > N.

992

TIEN-YIEN L I A N D J. A. YORKE

with k = n 2 + 1 and ~"'"(x,) and F"~+'(X~.) oth belong to [b2"-I, b*]. Therefore, 1 F"'"(x,)b ~"'"(x,.)) < F . Since there are infinitely many n with this property, liminf,,,/ F n ( x , ) - Fn(x,.)I = 0.
REMARK. theorem can be generalized by assuming that F: J- R without assuming that The F ( J ) c J and we leave this proof to the reader. Of course F ( J ) f l J would be nonempty since it would contain the points a, b, and c, assuming that b, c, and d, are defined.
Research partially supported by National Science Foundation grant GP-31386X
References 1. E. N . Lorenz, The problem of deducing the climate from the governing equations, Tellus, 16 (1%4) 1-1 1. 2. - Deterministic nonperiodic flows, J. Atmospheric Sci., 20 (1%3) 130-141. , 3. - The mechanics of vacillation, J. Atmospheric Sci., 20 (1%3) 448-464. , 4. - The predictability of hydrodynamic flow, Trans. N.Y. Acad. Sci., Ser. 11, 25 (1%3) 409-432. , 5. T . Y. Li and J. A. Yorke, Ergodic transformations from an interval into itself, (submitted for publication). 6. P. R. Stein and S . M. Ulam, Nonlinear transformation studies on electronic computers, Los Alamos Sci. Lab., Los Alamos, New Mexico, 1963. 7. S. M. Ulam, A Collection of Mathematical Problems, Interscience, New York, 1960, p. 150. 8. A. Lasota and J. A. Yorke, On the existence of invariant measures for piecewise monotonic transformations, Trans. Amer. Math. Soc., 186 (1973) 481-488. 9. T. Y. Li, Finite approximation for the Frobenius-Perron operator-A Solution to Ulam's conjecture, (submitted for publication). 10. Syunro Utida, Population fluctuation, an experimental and theoretical approach, Cold Spring Harbor Symposia on Quantitative Biology, 22 (1957) 139-15 1. 11. A. Lasota and P. Rusek, Problems of the stability of the motion in the process of rotary drilling with cogged bits, Archium Gornictua, 15 (1970) 205-216 (Polish with Russian and German Summary). 12. N. Metropolis, M. L . Stein and P. R. Stein, On infinite limit sets for transformations on the unit interval, J. Combinatorial Theory Ser. A 15, (1973) 25-44. 13. S . Smale, Differentiable dynamical systems, Bull. A.M.S., 73 (1967) 747-817 (see 01.5). 14. G. Oster and Y. Takahashi, Models for age specific interactions in a periodic environment, Ecology, in press. 15. D. Auslander, G. Oster and C. Huffaker, Dynamics of interacting populations, a preprint. 16. T. Y. Li and James A. Yorke, The "simplest" dynamics system, (to appear). 17. R. M. May, Biological populations obeying difference equations, stable cycles, and chaos, J. Theor. Biol. (to appear).