Choas Theory: Tent Map (Part 3)

Following from the previous post, we’re going to explore the limitations of iterating the tent map for when \mu = 2. The graph below shows graphically what occurs given a rational number, say, x_0 = 0.4. For clarity, we’ve zoomed in on the point in which the system jumps straight to 0. This is the point in which the binary notation for x is all 0’s.

51479eeaac8bda88aa302df90d8a0ca5.eps

Limitations when using \mu = 2

To see the behaviour of the tent map for all \mu across the range [0,1] of x we plot the bifurcation diagram (seen below).

Bifurcation diagram for the Tent map.

Bifurcation diagram for the Tent map.

Previously we stated that the behaviour of the tent map for \mu < 1 converges to 0 and so it isn’t of any interest to us, and so we’ve plotted the bifurcation for 1 \leq \mu \leq 2. Unlike the Logistic map, the tent map doesn’t follow the period-doubling route to chaos. In fact it can be shown that there are no period-doublings [1].

Tent map histograms: distribution of data

Tent map histograms: distribution of data

Similarly to the logistic map we can represent the distribution of x values and the range on a histogram. Note we cannot use \mu = 2 due to numerical limitations when computing the results as previously discussed. For \mu \leq 1 all the points converge to a fixed point and so nearly all the points are distributed to one value. For \mu > 1 the points converge to a number of fixed points as shown by the histogram for \mu = 1.2. As we can see from the graph above, the distribution for the tent map tends towards a uniform distribution as \mu tends to 2.

Sources:

[1] LAM, L. (1998). Non-Linear Physics for Beginners: Fractals, Chaos, Pattern Formation, Solutions, Cellular Automata and Complex Systems. London, UK. World Scientific Publishing.

Chaos Theory: Tent Map (Part 2)

Continuing from the previous post; the graph below shows the tent map for a range of \mu values, here we’ll use initial condition 0.4.

1

Tent map for a range of values of \mu

For \mu < 1 the system converges to 0 for all initial conditions. If \mu = 1 then all initial conditions less than or equal to \frac{1}{2} are fixed points of the system, otherwise for initial conditions x_0 > \frac{1}{2} they converge to the fixed point 1 - x_0. For example, for x_0 = 0.93 then it converges to the fixed point 0.07 as seen by the red points in the top right plot for \mu = 1 in the figure above.

For \mu > 1 the system has fixed points at 0 and the other at \frac{\mu}{\mu + 1}, we can show this mathematically using the tent map. For when 0 \leq x \leq \frac{1}{2} then the fixed point is when \mu x = x which implies that x = 0 is a fixed point. For when \frac{1}{2} \leq x \leq 1 then the fixed point is when \mu - \mu x = x. Rearranging gives \mu = x + \mu x = (\mu + 1) x. It follows that the fixed point is at \frac{\mu}{\mu + 1}. However, both fixed points are unstable, we can show this mathematically by looking at the gradient of the fixed point. If the gradient is less than one then its said to be stable, if its greater than then its unstable.

Let T(x^*) denote a fixed point of the tent map, then T'(x^*) is the gradient. Given the fixed point 0, then T'(x^*) = T'(0) = \mu. For the fixed point \frac{\mu}{\mu + 1}, then T'(x^*) = T'(\frac{\mu}{\mu + 1}) = \frac{1}{(\mu + 1)^2}. Since \mu > 1 in both cases then both fixed points are unstable. In other words, for a value x near the given fixed point, it will diverge away from it rather than converge. For \mu = 2, the system maps the interval [0,1] onto itself, becoming chaotic.

The dynamics will be aperiodic for initial conditions that are irrational, and periodic when rational. One effect of this is that we cannot run long term simulations of the tent map on a computer because all numbers in a computer are rational. Since x_n is expressed in binary notation, each successive iteration of the tent map will eventually hit 0 since the leftmost bit will always be removed (Ex 1). For irrational numbers there’s an infinite binary expansion and so will never go to 0 (Ex 2) [1].

(Ex 1) \quad x_0 = 0.5322265625 = \frac{1}{2} + \frac{1}{32} + \frac{1}{1024} = 0.1000100001.

(Ex 2) \quad x_0 = \frac{\pi}{10}.

In the next post we will (again) be exploring further the tent map, looking at computational limitations regarding the calculation of the tent map, bifurcation diagrams and histogram.

Sources:

[1] GLEICK, J. (1998). Chaos: the amazing science of the unpredictable. London, UK. Vintage.

Chaos Theory: Tent Map (Part 1)

In previous posts regarding chaos theory we have investigated the logistic map. In this post we’re going to be looking at a similar system known as the tent map; it is also commonly referred to as the triangle map.

The tent map is a recurrence relation, written as:

x_{n+1}=    \begin{cases}    \mu x_{n}, \qquad \qquad 0 \leq x \leq \frac{1}{2}\\    \mu - \mu x_{n}, \qquad \frac{1}{2} \leq x \leq 1 \qquad \qquad \qquad \qquad \qquad \qquad (2)    \end{cases}

for  0 \leq \mu \leq 2 and 0 \leq x \leq 1.

As shown below, the graph of (2) has maximum value \frac{\mu}{2} at x = \frac{1}{2}. The tent map is piecewise linear, because of this characteristic this makes the tent map easier to analyse than the logistic map. However, although the form of the tent map is simple and the equations are linear, for certain parameter values the map can yield complex and chaotic behaviour [1].

Tent map

Tent map

The logistic and tent map are topologically conjugate for r = 4 and \mu = 2 respectively. In other words, the behaviour of the tent map for \mu = 2 is the same as that of the logistic map when r = 4. The proof of this is will be looked at in future posts. First lets look to see if the tent map is sensitive to initial conditions. Here we’ll use similar initial conditions to that found in previous posts for the logistic map; x_{0} = 0.4 and x_{0}* = 0.40001 and \mu = 1.9999

Tent map: Sensitivity to Initial Conditions

Tent map: Sensitivity to Initial Conditions

Similarly to the logistic map, we can see that the tent map is also sensitive to initial conditions. In addition to being sensitive to initial conditions, the tent map is also dependent on its parameter value \mu, ranging from predictable to chaotic behaviour. We will explore the tent map further in the next post.

Sources:

[1] LYNCH, S. (2007). Dynamical Systems with Applications Using Mathematica (eBook). Birkhäuser Basel, 1st edition.

Ajanta Caves

The Ajanta Caves are among the world’s most prestigious religious locations in the world, sacred to Buddhism. Located in India, these caves were carved thousands of years ago from the mountain side; displaying inhuman-like capabilities as the tools and technology we have today would not have been accessible or invented when the caves were formed. The caves include paintings, sculptures and signs of astronomical knowledge and have been subject to wide-spread tourism and speculation as to how they were built.

This slideshow requires JavaScript.

The caves were discovered in 1819 by a British office on a hunting party. Hidden by dense forestation the caves were said to be lost for centuries after the downfall of Buddhism in India in the early millennium. The Ajanta caves were said to be a ancient monastery for Buddhists, including space for living, teaching and worship; this is reflected by the fact that most caves are connected by the exterior.

AjantaChait

Another fascinating feature of these caves are their relationship with the summer and winter solstice’s. The sun shines directly through gaps in the mountain to light up and illuminate significant religious Stupa’s (a place for meditation) at the back of the caves. It’s amazing because precise calculations and tools would have been needed to be able to orient the caves to the solstices since the caves are carved inside the rock. This has led people to believe that alien activity or interaction was involved in the feat of constructing the caves. There’s plenty of information regarding these caves and documentaries surrounding them. Check out the video in the sources below!

Sources:

http://www.youtube.com/watch?v=FHeel1GNFFY

Chaos Theory: Periodic Windows

The bifurcation diagram is much more complex than just a simple division of regions of chaos and periodic behaviour [1]. As previously stated, in this post we will be looking at one of the most iconic features that comes from observing the bifurcation diagram; periodic windows. Amongst all the chaos, a stable period-3 orbit appears at r \approx 3.83 as shown below. These regions (periodic windows) exist for r > r_{\infty}, they are areas whereby the system isn’t chaotic.

Periodic window at approx. r=3.83

First we will explore how the period-3 orbit is created. Note, all periodic windows are created in this manner [2]. Using the logistic equation x_{n+1}=rx_{n}(1-x_{n}) then we can rewrite as x_{n+1} = f(x_n), so x_{n+2} = f(f(x_n)) = f^2(x_n), similarly x_{n+3} = f^3(x_n). This mapping is vital in understanding the birth of the period-3 window. In the figure below we plot the third iterative mapping for r = 3.84 which has 8 fixed points such that f^3(x) = x.

Third iterative map of the logistic map for r=3.84

Third iterative map of the logistic map for r=3.84

Two of the fixed points are actually solutions to f(x)=x, the other six are marked by dots coloured red and black. The red dots correspond to a stable period-3 orbit, alternatively, the black notes correspond to a unstable period-3 orbit. We can show this is true by looking at the gradient at each fixed point. If it’s negative then it’s stable, and if positive then it’s unstable. Decreasing r to r = 3.8 then the graph changes as indicated by the red line below.

Third iterative map of the logistic map for r=3.8, 3.8284, 3.84

Third iterative map for r=3.8, 3.8284, 3.84

Therefore we expect for some value between r = 3.8 and r = 3.84 the graph of f^3(x) becomes tangent to the diagonal. At this value of r then the stable and unstable fixed points unite and annihilate in a tangent bifurcation [2], indicating the birth of the period-3 window. This value of r is approximately represented by the green line above, and is said to be r = 3.8284... [3]. As previously stated we can show that each periodic window is created in this way. In addition, it can be shown that there are periodic windows in all chaotic regions and so is a infinite number of periodic windows [4].

Sources:

[1] GILMORE, R. & LEFRANC, M. (2002). The Topology of Chaos. New York, USA. John Wiley & Sons.

[2] STROGATZ, S.H. (1994). Nonlinear dynamics and chaos. New York, USA. Westview Press, Perseus.

[3] TAKAYASU, H. (1990). Fractals in the physical sciences. Manchester, UK. Manchester University Press.

[4] FELDMAN, D.P. (2012). Chaos and Fractals: An elementary introduction. Oxford, UK. Oxford University Press.

Chaos Theory: Universality – Feigenbaum constant and U-sequence

In this post we’ll be investigating some amazing results gained from non-linear dynamics. First we’ll be looking at the Feigenbaum constant, a quantitative universality followed by the U-sequence, a qualitative universality [1],[2].

       Feigenbaum constant

We can make some extraordinary observations from the bifurcation diagram. The ratio between each bifurcation converges to Feigenbaum’s constant, a universal constant for functions that have a periodic-doubling route to chaos, a fact discovered by Mitchell Feigenbaum in 1975. This constant is approximately 4.669 [2]. We can interpret from this constant that as we approach chaos each periodic region is smaller than the previous region by a factor of 4.669 (the Feigenbaum constant).

Bifurcation diagram showing source of the Feigenbaum constant

Diagram showing source of the Feigenbaum constant

The extraordinary thing about this constant is that its the same for all uni-modal functions, i.e. functions having a quadratic maximum and approach chaos via period doubling. If we call r_n the value of r for which the 2^n period becomes unstable then \delta = \lim_{n\to\infty}\frac{r_{n+1}-r_n}{r_{n+2}-r_{n+1}}=4.669.... The plot above shows the source of the Feigenbaum constant, the red lines indicate where the logistic map bifurcates to a period 2, 4 and 8 orbit respectively.

      U-sequence

The universal sequence, or U-sequence is commonly demonstrated on iterated maps whereby the mathematical model is a uni-modal mapping such as the logistic map [1,p.100]. In 1973, Metropolis, Stein and Stein investigated periods in uni-modal maps, proving that stable periodic solutions appear in an order independent of the map being iterated. In other words the periodic attractors always occur in the same sequence [2]. Up to period 6 the U-sequence is 1,2,4,6,5,3,6,5,6,4,6,5,6 [2]. The values represent the periodic-orbit of the mapping, so the first numbers show that the system goes from a period 1 to a period 2, to a period 4 and so on (see the bifurcation diagram for the logistic map in previous posts to view locations of these periodic orbits). We have highlighted the position where the period-3 window appears, a conspicuous occurrence and feature of the bifurcation diagram; something we will be exploring in the next post. To summarise, these phenomenal results imply that the form of a function isn’t relevant, however, the overall shape is.

Sources:

[1] ZELINKA, I. (2010). Evolutionary Algorithms and Chaotic Systems (eBook). Berlin, Germany. Springer.

[2] STROGATZ, S.H. (1994). Nonlinear dynamics and chaos. New York, USA. Westview Press, Perseus.

The Legend of the S.S. Ourang Medan

We’ve all come across stories of ghost ships and death ships: The Flying Dutchman, Mary Celeste, and so on… However, (in my opinion) none really compare to the nature of the story I’m about to tell you. That is, the legend of the S.S. Ourang Medan.

Image

According to widely circulated reports, in June 1947, multiple ships traversing along the Strait of Melacca picked up on SOS distress signals. The message was as simple as it was disturbing:

“All officers including captain are dead, lying in chartroom and bridge. Possibly whole crew dead”.

What followed was a series of undecipherable morse code and a tomb-like silence after a grim final message:

“I die”.

Image

These distress signals were picked up by two American ships alongside British and Dutch listening posts. The men manning the posts were able to triangulate the position of the signal at the location of the S.S. Ourang Medan. An American Merchant ship, the Silver Star set course towards the ship in response to the urgent SOS. Within hours, the ship caught sight of the incapacitated ship, bobbing up and down hauntingly on the waves of the sea. When nearby, there was no sight of life on deck so the Silver Star decided to dispatch a party to investigate. Little did they know that they were going to confront a living nightmare.

Once on board, the crew quickly realised that the distress signal wasn’t an exaggeration. The deck was littered with corpses; with their eyes and mouths wide open, arms out, faces twisted with the anguish of horror and agony. It is said that even the ships dog was dead, frozen into a haunting grimace. The captain’s remains were discovered on the Bridge, the communications office was still at his post, his fingertips resting on the telegraph. According to reports, all the corpses shared the same wide-eyed expressions on their faces. The cause is still unsolved, reports state that there was no evidence of foul play or damage to the ship – just a ghastly cold chill consuming the atmosphere.

Upon making their way back to the Silver Star, the crew witnessed smoke billowing from the S.S. Ourang, the boarding party were able to make it to safety just before the S.S. Ourang mysteriously exploded. The ship was lost to the sea, and became a part of the realm of maritime legends and mysterious.

Image

As far as what happened to the ship and how it happened, nobody knows. However, if someone does know the truth about it, then there not coming forward. While it may not be as famous as other ghost ship stories that are circulating, it is just as terrifying and mysterious.

Sources:

http://mysteriousuniverse.org/2011/11/death-ship-the-ourang-medan-mystery/

http://www.unexplained-mysteries.com/forum/index.php?showtopic=248895