Horseshoe Chaos

While I struggle with supposedly impossible inversions, I'd like to introduce you all to Smale's horseshoe. This object is where the definition of chaos originates yet it is a relatively simple idea.

You start with a square. Squish is down and pull it out to the sides then bend it at the halfway point so it looks like a horseshoe. Then replace it over the square so that the two lengths form vertical strips on the square. This is the basic transformation of the horseshoe map. It is also invertible. Take the bent strips, rotate back to horizontal, and unbend it then squish it in from the sides and pull from the top and bottom then you get the square again. If you do this inversion again you get strips intersection horizontally.

Repeatedly performing the forward transformations you get more and increasingly thin vertical strips: 2 strips become 4, 4 strips become 8 i.e. each strip gets two thinner strips within it at each iteration. Doing the same backwards, you get lots of increasingly thinner horizontal strips. The points in the strips are the points that remain in the set at that iteration. Most points leave after even just a few iterations in either direction.

If you overlay the two directions, the points at the intersections that remain in the square under all iterations (if you iterate infinitely in either direction). This is the invariant set. This set is a Cantor set (disconnected and unstable fractal).

If you labelled each point with 0 or 1 at each iteration you could describe the positions of all the points in the square e.g. .0101 means a point that is in the 6th vertical strip after four iterations (the right hand strip of the left hand strip of the right hand strip of the left hand strip in the first iteration) or

If you can see that...sorry if it's too small.

Since you can do this in both directions, the whole dynamics of any point in the square can be described by it's bi-infinite sequence of positions.
Knowing this we can construct any orbit we like and we know it will describe at least one point in the square.
So we can make a periodic point: ...010101.0101010... which alternates from the left to the right. Since we have infinite lengths in both direction, we can make (uncountably) infinitely many of these periodic points.
Or we can make a non-periodic point: ...010001101100000.1010100110111... by, for example, 'counting' in binary as above or just by adding the wrong value to a periodic point i.e. ...01010101.0101011...is non-periodic. There are (uncountably) infinite of these too.
Among this set of non-periodic orbits we can find at least one that comes arbitrarily close to the invariant set. This is trivial since you can take just part of the sequence that describes a point in the invariant set and it would come as close as you like!

I hope you are following since I have just demonstrated that Smale's horseshoe perfectly defines chaos.
Chaotic motion must consist of the following:
- infinitely many periodic orbits
- non-periodic orbits
- dense orbits - ones that come arbitrarily close to any point in the set.
There tends also to be the condition of sensitive dependence on initial conditions but this can be shown through the dense orbits - you can take two sequences that are the same for any arbitrary iterations and in the next they can end up in two different strips i.e. they may start very close but can end up very far apart.

Voila!

Strogatz does give a very basic introduction but for a more detailed approach read the first chapter of "Elements of Applied Bifurcation Theory" by Yuri Kuznetsov. Wiki also has its own explanation.

How many drafts...?!

So here we are, nearing the end of my project. Well, actually there are six weeks left til I hand in my thesis but it doesn't feel like it. It feels more like two. :S

I've been asked to do two drafts already. A third one is coming up next week and I'll redo the current one into a fourth. After that they combine into a fifth. Then will follow a sixth from recommended improvements. I'll probably create another couple just for kicks! So that's at leasts one a week...a draft roughly for each month I've worked on it. I'm drowning in them!!

Ah well, at least it means I have shiny pictures for you this week. My draft task this week was to tell the story of one side of my research (decreasing lambda) as purely in images as possible. Obviously the details will be missing but I hope you enjoy the obvious progression and changes in images as lambda varies.
My image report draft is available here.

In other parts of my Masters world I am winding down on lectures and winding up for revision ready for exams in nine weeks. I want to try something new for my revision this year - sort of step it up a level. If anyone has any recommendations that would be very welcome. I have a few ideas but would like to augment these with other influences.

 I also have to produce a poster and a website for my Masters project but I'll share them here first to get opinions! I will of course add a link to this blog in my website (cheeky I know!); maybe they'll give me extra marks for initiative and knowledge sharing never mind the pretty pictures!

Inconvenient Delays

So here we are a month since I last posted. I would explain but that would feel like giving excuses and I'm trying to get out of the habit of feeling guilty for things which were out of my hands. I do however have something for you to read and perhaps enjoy.

My most recent deadline has been for a draft masters thesis for my project. This is available for you to read via this link.

What I have written covers about a third of the area of work I have been doing. It is very minimalist and basic with very little of the style or continuity that I duly put in my work. My introduction and background are perhaps not as full as I would like but maybe you can provide some thoughts on that since it is aimed at bringing people of different ranges of knowledge on the subject up to speed before I share my analysis. If there is any point of logic missing please do let me know.

The version I am sharing is the version I would have like to share with the examiners but my supervisor had other ideas so I actually submitted a stripped down version that is even worse so I will not be sharing that.

One point which you should all love is all the images generated that I have included in the report. In my time I have analysed the Julia set changes, the stability changes, the manifolds and most recently the full bifurcations in one half with the other half due to be completed this week. I have finally reached the crux of my examination of the map which is a strange one. The circle hole in one half has an important bearing on the location of the attractor which "disappears" and yet still points remain in the set - this is only possible with an attractor still present yet at some point there doesn't even exist a higher period attractor, just saddles and repellors. Thus this week my question to answer is whether an attractor still appears to exist for the points under the mapping despite the properties attributed to it or if the boundary of the hole itself becomes the attractor or neither of these things (which would mean either a waste of half a year's work or something totally cool and new).

I'll let you know next week hopefully and also provide you will the full bifurcation diagram.

Nonlinearity gets all the fun

Nonlinearity is inherent in most systems before you make an assumption and reduce the beautiful complexity down to a boring bog standard linear system. There is a section in Strogatz that shows probably the most interesting thing you can do with linear systems (Chapter 5 section 3 page 138 onwards). And I only read it once because linear systems were dull and pretty repetitive long before A level.

Nonlinearity is pretty easy to understand. Take for example square numbers. You can make 4 by either multiplying 2 and 2 or -2 and -2. This means that, under the square operator, 4 has two (well, four really, but two are identical to the other two) possible previous values. This makes squaring a number nonlinear because tracking back you do not know whether two 2s or two -2s made the 4 you see at the moment. This particular brand of nonlinearity is called noninvertibility. You still following me?

Good. Right, to take this further we call the previous values of each number pre-images and the operation, a mapping. Now the quadratic (squaring) mapping that I am examining is:

what happens to z when a=2 and I vary lambda. This mapping does funny things depending on what I set lambda to. If it is between 0 and 1 then it behaves just like the squaring operator - each value of z has two pre-images as I iterate the map but in a circle radius=1-lambda around C the values have no pre-images. Now if I make lambda greater than 1, some values have two pre-images as before but the z values in the circle now have four pre-images.

Still with me? We can take this further but placing this property on curves not just points. A curve is just a line on a graph (an infinite set of points). Instead of looking backwards as we did with the points, we instead go forwards and see how the curve changes. It's first image is a circle, say. This curve then becomes the pre-image of the set of points of the circle under the mapping. And so on, forwards in time. It helps to look at the image below:

You can see the curves collecting on (heading towards) one curve called a manifold (it's attracting). Now consider the relationship of the points and the curves. Because a curve is just a set of points. With noninvertibility inherent in the system, what happens it you take the manifold and reverse the mappings for all the points? Wild Chaos! Tracking the possible points and curves and how and why they interact is the next stage of my project - at least for certain ranges of values for C, a and lambda. This limitation is purely because my project has a time limit of a year.

I am currently working with a PhD student, Stefanie Hittmeyer (who incidentally produced the image above) who is working on the full problem in all its scope. She has been focused on the manifold part of it for the moment while I have been examining the fractal structure side. My work nicely compliments hers. This work is on the fore-front of chaos research. Now it may seem that all this has little value in the real world but the great thing about nonlinear and chaotic systems is that it can accurately mimic real world systems. This is a major advantage over linear systems because when taking a real world problem you have to reduce the picture you are looking at so you can describe it using the simplest tools. With this kind of research, no such reduction is necessary and you can understand the system in it's full complexity. So maybe it'll turn out that neurological signals in the brain will turn out to have similar properties to the mapping that I am examining and my research will contribute to a fuller understanding of how we came to have consciousness and help neurosurgeons get better at fixing the biochemical circuitry that makes us who we are. Or maybe, and more likely since this mapping is a reduction of the five dimensional Lorenz model, it will just help meteorologists get really, really good at predicting weather patterns.

Who knows? It's sort of like the laser - a solution created with no problem that has now become integral to our lifestyle. Just check out how many things it is used in! I hope you managed to follow all that to the end. If you have any specific questions please feel free to email me but I would much rather you follow that urge to take a deeper look yourself (even if that is just reading all the wiki pages). It's way more fun!