Monday, June 22, 2015

The Aizawa attractor V

I now have an oscilloscope (a very nice HP 1740A). Here's the Aizawa attractor with computer in FAST mode (i.e. 1 nF integration capacitors) and a 2 second exposure:

Aizawa attractor. 1 V/div both axes; vertical axis is z, horizontal is x.
Note horizontal screen central axis is +1 V.
Aizawa attractor. 1 V/div both axes; vertical axis is y, horizontal is x.
Note horizontal screen central axis is 0 V. 

Thursday, June 18, 2015

The Aizawa attractor IV

Not wishing to underestimate the power of the Dark Side, I've produced some renderings of the Aizawa attractor via MATLAB.

These were obtained using a very simple approach (Euler's method - the (x,y,z) coordinates of the current point are pushed through the differential equations, the resulting new values for the time-derivatives of (x,y,z) then multiplied by a time step, and then added to the current point, to give the next (x,y,z) point, and so on).

The following shows the xz-plane for a time step of 0.0001.


The published examples of the Aizawa attractor (for example, [1]) seem to use 0.01 as the time step - but it's apparent that the size of the time step has (perhaps not surprisingly) some effect on the solution.

The results below are for different time - step sizes (0.01, 0.001, 0.0001 seconds). In all cases the solutions show the first 600 seconds of the trajectory (i.e. for the 0.01 second time step, the equations were iterated 60,000 times, producing 60,000 points, and so on). The parameters are as in the previous posts (α = 0.95, β = 0.7, γ = 0.6, δ = 3.5, ε = 0.25 and ζ = 0.1), with initial conditions also as in previous post (x = 0.1, y = z = 0).

These tentative results tend to suggest that the 0.01 time step is a bit coarse: the 0.001 and 0.0001 cases look similar to each other, but different to the 0.01 case. Interestingly, the central tube-like structure becomes thicker (and more open) when the time step is smaller. Also shown is the attractor in the xy-plane.


The following shows the xy-plane for the 0.0001 case.


Evidently, the tube-like opening, although apparently circular, is not centred on the origin. This appears to agree with results obtained from the analog computer:
xy-plane, 10 minutes' of trajectory

Reference

[1] http://jlswbs.blogspot.co.uk/2011/10/aizawa.html


Tuesday, June 16, 2015

The Aizawa attractor III

Here's my patching for the Aizawa attractor (see previous posts):


This uses all three multipliers, the divider (as a fourth multiplier), and the two squarers. The cubes are obtained directly using the (AD538 based) nth power circuits.


Amazingly, this all worked first time. Still sans oscilloscope - will be interested to see the attractor with the computer in FAST mode. In the meantime, here is the plotter output (1 µF integration capacitors). The plot shows the first ten minutes of the trajectory.

2 V/cm both axes, 10 minutes' worth of trajectory.

Pleasingly, this agrees pretty well with other published images of the Aizawa attractor e.g. [1], although the latter appears to be an inverted mirror image of what I produced - perhaps the axes in [1] are reversed. The image of the attractor given in [2] does have the same handedness as mine.

Clearly, the attractor has some degree of symmetry about the z-axis. In the previous posts, I deduced that the attractor has three equilibrium points, lying along the z-axis:


Now, the first one (z = 1.9432644 V) does appear to correspond to the top of the attractor. But, the other two are well below the attractor. The bottom of the attractor is at about -0.58 V. This discrepancy requires some thought...


Reference

[1] http://www.algosome.com/articles/images/aizawa-attractor.jpg

[2] http://i.imgur.com/wMrmY.jpg

Sunday, June 14, 2015

The Aizawa attractor II

In the previous post, it was apparent (from numerical experiment) that the Aizawa attractor has three purely real equilibrium points along the z-axis, so long as parameter α is greater than a certain value, namely 0.932169751786... (with the other parameter values as listed in the previous post). If α is less than this value, then there is evidently one real equilibrium point, and two complex equilibrium points.

Here we derive an expression for α which gives the demarcation point between all-real and real and two complex equilibrium points.

The two (potentially) complex roots / equilibrium points are given by (see previous post):
Assuming that λ is not itself complex (see below), the above will roots will be purely real if λ2 = α (since this forces the right hand side of the above to be zero). From the previous post we have
Hence we need to solve 
which, by inspection, has the solution

giving the value of 0.932169751786... which agrees with the value determined by numerical experiment.

Note that although the above gives the value for α for which the imaginary parts of the equilibrium points are exactly zero, it is not obvious that values of α greater than this will give purely real equilibrium points, since if α > α0 then λ will itself be complex, because 9γ2 - 4α3 will be negative, and the above argument / derivation of α0 assumes λ to be real.

The Aizawa attractor I

There's a rather beautiful attractor called the Aizawa attractor: a sphere with a tube-like structure penetrating its axis [1] - a good candidate for the analog computer, in anticipation of the arrival of a new (old) oscilloscope from France.

The Aizawa attractor's equations are:


Given that there's precious little literature about this attractor, I have to plough my own furrow; I begin by searching for equilibrium points.

The game seems to be to set the left hand sides of the three equations to zero. Interestingly, by inspection, it's possible to solve the first two equations independently of the third, x = y = 0 being an obvious solution. In fact there are three solutions to the first pair of equations, obtainable, for example via:

Of these solutions, the first two can only yield real equilibrium points if x = y = 0; the third solution yields a complex equilibrium point, unless δ is zero or purely imaginary. If δ is zero then the first two equations would be decoupled. I am not sure what δ being imaginary would mean in terms of the analog computation - perhaps interpretable as a delay?

Given that equilibrium points should be real, since they exist in a real state space, we conclude that x = y = 0 for any such points - i.e. they must lie on the z-axis.

Setting x = y = 0 in the third equation, gives
which has three solutions:

It is not obvious that any of these roots are purely real. But, for certain values of  α and γ it transpires that they are... 

...in all the published cases I can find, the Aizawa attractor uses the following parameter values: α = 0.95, β = 0.7, γ = 0.6, δ = 3.5, ε = 0.25 and ζ = 0.1, and plugging these values into the above equations gives the following (entirely real) roots for z:

Hence, for the parameter values listed, we expect three real equilibrium points, along the attractor's z-axis.

Now, clearly, λ depends solely on α and γ, and, as suspected, it's possible to adjust the value of (for example) α, so as to conjure up some complex values for z. A little experimentation shows that two of the roots become complex if α is reduced slightly, for example α = 0.93, yields:

which suggests just one (real) equilibrium point.

Indeed, it appears that the first root is always real, regardless of the value of α. 

Further experimentation indicates that α  ≈ 0.932169751786 is the point at which the second and third roots start to have a non-zero imaginary part. This value is intriguingly close to the 0.95 value which seems to be the de facto value for α used for the Aizawa attractor. I assume this is not a coincidence, and that the attractor 'needs' the three (real) equilibrium points to exist. It would be interesting to know how the 0.95 value was found.

Vieta's formula for the sum of roots of a polynomial tells us that the sum of the three roots of
(i.e. equilibrium positions along the z-axis) must equal zero (since the quadratic term in z is zero). Hence, if one root (equilibrium position) is real, then the others must either both be real, or both be complex (with equal but opposite imaginary parts). It would appear not to be possible to have just a single complex root (which would mean exactly two real equilibrium positions). 


Reference

[1] www.algosome.com/articles/aizawa-attractor-chaos.html

Tuesday, June 9, 2015

FitzHugh - Nagumo and Fuss - Catalan II

After yet more noodling I have the following result for the exact x1 coordinate of the equilibrium position of the coupled FitzHugh-Nagumo system (see previous posts), this time valid for all values of k:
As in the previous post, C3(n) is the nth Fuss–Catalan number of order 3.

Sunday, June 7, 2015

Coupled FitzHugh - Nagumo system and the Fuss - Catalan numbers

After some mathematical noodling (whilst awaiting a replacement oscilloscope), I've come up with two bizarre expressions for the exact equilibrium position of the coupled FitzHugh-Nagumo attractor (as per previous posts), for the cases in which coupling parameter k = -1 and 0 respectively:

C3(n) is the nth Fuss-Catalan number of order 3, or ternary number [1], the first few terms being 1, 1, 3, 12, 55, 273, 1428, 7752...

Given x1, we can readily obtain y1 from (as per previous post):


Example

Using the above equations for x1, with a = 0.4, b= 0.5, c = 0.01, as per Sprott's figure 6.7, we calculate

x1 = -0.7901341941135... (k = -1 case), and

x1 = -0.8056544203709... (k = 0 case).

These agree with the values for x1 given by equations (7) and (8) in the previous post.

Reference

[1] Sequence A001764, N. J. A. Sloane, A Handbook of Integer Sequences, New York and London: Academic Press (1973). Online at https://oeis.org

Coupled FitzHugh - Nagumo system - equilibrium position

The coupled FitzHugh-Nagumo oscillators produce an attractor (previous post) which has an obvious resting / equilibrium point. The trajectory spends a long time spiralling about this point before making one or more fast excursions, before returning to the equilibrium point, and so on.

The following shows the location of the centre of the equilibrium region for Sprott's figure 6.7 [1]. I've added my estimates as to the coordinates of the equilibrium centre:


Given that the trajectory moves slowly around the equilibrium point, it seems reasonable that one could work out an analytic solution as to the whereabouts of this point. In pursuit of this, I took Sprott's coupled FitzHugh-Nagumo system (equation (6.7)), and set the right hand sides to zero (i.e. the game being to find where the speed of the trajectory is zero):


After a fair bit of algebra, I worked out the following rather compact equations for the coordinates of the equilibrium position (see box below). These equations are valid for small values of c - consistent with Sprott's comment that 'chaotic solutions require that c be rather small' [1]. These equations were derived from a power series development of an exact solution for the coordinates, which was worked out, the derivation of which is given also, below.

It is evident that the y-coordinate does not depend on parameter c, when c is small.

The following tables gives the calculated coordinates of the equilibrium point, for k = -1 and k = 0 (no coupling). All other values are as per Sprott's figure 6.7 (i.e. a = 0.4, b = 0.5, c = 0.01). The first table gives the results obtained using the exact solution, the second gives results obtained using the above, boxed, equations.

Coordinates of equilibrium point calculated using
exact derivation.
Coordinates of equilibrium point calculated using
approximate derivation.
These agree rather nicely with the values I measured directly off Sprott's figure 6.7 - namely -0.78 for x1 and +0.51 for y1.

What's interesting (to me) is that the equilibrium point actually depends on the value of the coupling parameter k. I (perhaps naively) assumed that the equilibrium point for a single decoupled oscillator would be the same as that for the coupled case - but experiments with the analog computer indicated otherwise:

Analog computer results for coupled (k = -1) and uncoupled (k = 0) cases. It is evident
that the equilibrium points have very different y-values.
In particular the y-coordinate of the equilibrium point varies strongly with k.

In contrast, the x-coordinate of the equilibrium position does not depend strongly on the coupling parameter k. To first order, the x-coordinate is dominated by the -a/b term which, in this case is equal to -0.8.

For reasons which remain unexplored (offset voltages?), the analog computer's results are slightly different to Sprott's. My equilibrium positions are shifted with respect to the origin of the plot.

Coordinates of equilibrium point measured from
analog computer output plot.



Derivation of Exact Coordinates of Equilibrium Point

 
Note

The ninth order (degree) polynomial appears to have just one real root. This is equation (7) above.

Reference

[1] Sprott J.C. Elegant Chaos: Algebraically Simple Chaotic Flows, pp. 123 - 125.