Posted by: Dave Richeson | January 1, 2010

Geogebra applet for families of discrete dynamical systems

As I mentioned recently, I taught the last two weeks of my colleague’s differential equations course. The topic was discrete dynamical systems. I posted links to a few Geogebra applets that I made, namely, applets for illustrating one-dimensional dynamical systems and an applet to generate cobweb plots for the logistic map.

This is my third and final post on this topic. I created a Geogebra applet that generates cobweb plots for various families of discrete dynamical systems. [Note: for some reason, I can’t get the Javascript to work on Safari on my Mac, but it works with Firefox. If anyone can help me solve the problem I’d be very grateful.]

In class we talked about families of discrete dynamical systems: functions f_k:\mathbb{R}\to\mathbb{R} in which k is a parameter. In particular we focused on bifurcations. From our point of view we say that the family has a bifurcation at a parameter value k_0 if the number or nature (attracting or repelling) of periodic orbits changes as k passes through k_0.

We focused on the following types of bifurcations.

Tangent bifurcation. A pair of fixed points are born—one attracting and one repelling.

Pitchfork bifurcation. A single attracting (or repelling) fixed point splits into three fixed points, one repelling (attracting) and two attracting (repelling).

Period doubling bifurcation. An attracting period n orbit becomes repelling and gives birth to an attracting period 2n orbit. These period doubling bifurcations happen one after another, so we get attracting periodic orbits of period n, then 2n, then 2^2n, then 2^3n,.. etc. This is called a period doubling cascade.

Below is a bifurcation diagram for the logistic map (thanks to Wikipedia), which you can click to enlarge. The vertical axis is the interval [0,1] (the phase space) and the horizontal axis is the parameter, which we call k, but is labeled r. The plot is generated by iterating the dynamical system 100,000 times (with seed x=0.25) for each parameter value, but skipping the first 1,000 iterates. This gives us a way to see the long-term behavior, such as movement toward an attracting periodic orbit or chaotic motion.

You can clearly see several period doubling cascades. The first one, from k=2.4 to 3.58 starts with a fixed point, which doubles to a period 2 orbit, then period 4, and so on. In fact, all of the vertical white bands are period doubling cascades. For example, the large white band at about k=3.85 shows a period doubling cascade with orbits of period 3, 6, 12, 24,…

As I mentioned above, in each case these periodic orbits are attracting, that’s why we can see them. However, there is much more going on in these dynamical systems. For example, consider the logistic map with k=8.4. There is an attracting period 3 orbit. As I explained in my earlier post about Sharkovsky’s Theorem, the existence of a period 3 oribit implies that there are periodic orbits of every period!

By the way, it was neat to see two follow-up posts to my dynamical systems posts from other bloggers: Dan MacKinnon at mathrecreation wrote a curious population model and Mike Croucher at Walking Randomly wrote Using SAGE to investigate the discrete logistic equation.

About these ads

Responses

  1. Very nice set of web diagrams with Geogebra. Thanks.

  2. Dave, do you know how one could use geogebra to plot

    y” / { [ 1 + (y’)^2 ]^1.5 = -M/EI

    M is a function of x, Eand I are constants.
    I had hoped solveode would work but am uncertain how.
    Any help would be appreciated.

    Thanks
    Tim

    • Sorry, I don’t. I haven’t used GeoGebra to solve any DE’s. Their help forums are pretty good. You could ask there.


Categories

Follow

Get every new post delivered to your Inbox.

Join 220 other followers

%d bloggers like this: