next up previous contents
Next: Sarkovskii's theorem and the Up: Interval Self-mappings Previous: Down the waterfall

Falling all the way: Cascade diagrams

The fixed points, 2-cycles, and 4-cycles we have seen are all examples of attractors of the forward orbit under f. This concept is the same as we have seen before: the attractor of a sequence of iterates is the set A of all points of accumulation of the orbit. The stable set tex2html_wrap_inline9412 of the attractor A is the set of all tex2html_wrap_inline9416 such that the set of points of accumulation of tex2html_wrap_inline9418 is A. A is a stable attractor if tex2html_wrap_inline9412 is nonempty and open. Here we need to quote a fact gained from deeper analysis:

Theorem: For a fixed tex2html_wrap_inline9426 , there is at most one stable attractor A, and the stable set of A is dense in I.
If the attractor is a cycle, that is the only attracting cycle.

The stable attractor of a dynamical system may be discovered, at least numerically, by simply letting the system run long enough. This is the key to unveiling the most fantastic fractal yet seen in this course. For each tex2html_wrap_inline6624 , we choose a seed tex2html_wrap_inline4950 (it doesn't matter which), and then we compute say 1000 iterates tex2html_wrap_inline9438 . We ignore what happens early on in this sequence; however, once we reach tex2html_wrap_inline9440 , we then record the next 1000 iterates tex2html_wrap_inline8812 , tex2html_wrap_inline9444 . By letting the system run for 1000 iterations, we are guessing that after that much time we are numerically extremely close to the attractor. Then the next 1000 iterates will run over all the points in the attractor. If the attractor is a stable 4-cycle, for instance, the 4 points in that cycle will be recorded 250 times each.

We may now plot the attractors for all tex2html_wrap_inline6624 in a single diagram: the cascade diagram (also called the bifurcation diagram). The horizontal axis measures tex2html_wrap_inline6624 running from 0 to 4; the vertical axis measures x from 0 to 1. The vertical cross-section is the attractor of tex2html_wrap_inline9202 . Generally speaking, a great deal of numerical computation is involved in the creation of this diagram. Fortunately, FRACTINT already has a beautifully optimized version of this program under the fractal type biflambda (the ``bifurcation diagram'' for the ``lambda'' map). The beautiful cascade diagram created by FRACTINT is shown in this figure. To consult the FRACTINT notes on the fractal type biflambda, see the following:

FRACTINT notes on biflambda

   figure2806
Figure 69: The cascade diagram

At tex2html_wrap_inline9454 , we can see the bifurcation as a joining of the curves

align2813

There is also a ``cascade of bifurcations:'' 2-cycle splitting into 4-cycle, 4-cycle splitting into 8-cycle, etc. This can be better examined by using FRACTINT's zooming mechanisms. For a windowed FRACTINT, simply press and hold the left mouse button to drag out a zoom window. Releasing the button creates a larger view of that part of the fractal.

In particular, let's zoom in on the cascade of bifurcations. In the ``Basic Options'' menu of FRACTINT, you may want to set maxiter to 1000, to obtain a more accurate portrayal of the attractors. The result of one zoom is shown in this figure.

   figure2833
Figure 70: The bifurcation diagram up close

Here is another zoom to tex2html_wrap_inline9462 , 0.86<x<0.90.

After all the bifurcations from tex2html_wrap_inline9466 -cycles to tex2html_wrap_inline9468 -cycles, we reach a limiting value of tex2html_wrap_inline6624 where at last there is an infinite attractor, around which the iterates bounce endlessly. The end of the waterfall is commonly considered the beginning of Chaos.

There is clearly a sizeable degree of self-similarity in the waterfall of bifurcations. In fact, careful measurement of the values of tex2html_wrap_inline6624 where the bifurcations occur led Mitchell Feigenbaum to discover an astonishing constant of nature. The table below records the values of tex2html_wrap_inline6624 where bifurcations occur, as measured by FRACTINT. With FRACTINT, it's possible to record just a few digits of accuracy. Set the ``filter'' number of iterations quite high to obtain the best accuracy. Also, in the Windows version of FRACTINT, be sure to activate the coordinate window to be able to read off the coordinates of the bifurcations.

 

tex2html_wrap_inline9476 tex2html_wrap_inline9478 tex2html_wrap_inline9480
tex2html_wrap_inline9482 3.00000
tex2html_wrap_inline9484 3.44949 0.44949
tex2html_wrap_inline9486 3.5441 0.09461 4.75098
tex2html_wrap_inline9488 3.5644 0.0203 4.66059
Table 8: Table of bifurcations: tex2html_wrap_inline4864 and tex2html_wrap_inline4866

 

The ratios of the gaps in the fourth column of Table 8 was discovered by Feigenbaum to have a limit value 4.6692.... Feigenbaum also sought the bifurcation point for other functions that were ``hump-shaped,'' such as tex2html_wrap_inline9496 and the ``roof'' function

displaymath9498

In all cases, the limiting ratio of the gaps between bifurcations turned out to be 4.6692.... The story goes that Feigenbaum was so excited by the universality of his number that he immediately called his mom and related how this number would make him famous. Fascinating historical details may be found in [Gle88]. The proof of one version of the universality of Feigenbaum's number was given in [CE80]. FRACTINT allows the user to vary the function in the bifurcation diagram; try the fractal types bif=sinpi, bifurcation, etc.

Amazingly, the bifurcation diagram shows that there are moments of tranquility beyond the threshold of chaos. We see sizeable ``gaps'' where the attractor again becomes a finite set. If we zoom in very carefully, there is a small region where the attractor has exactly 3 points; that is, a stable 3-cycle occurs. The interval of 3-cycles appears to be tex2html_wrap_inline9502 . Let's try tex2html_wrap_inline4868 to see. The web diagram is shown in this figure.

   figure2868
Figure 71: A periodic point of period 3: tex2html_wrap_inline4868


next up previous contents
Next: Sarkovskii's theorem and the Up: Interval Self-mappings Previous: Down the waterfall

David J. Wright
Mon Aug 19 17:21:15 CDT 1996