Preface


This section gives an example of the chaotic dynamic system known as the Rössler attractor. It is a system of three non-linear ordinary differential equations originally studied by Otto Rössler (born 1940) in the 1970's. These differential equations define a continuous-time dynamical system that exhibits chaotic dynamics associated with the fractal properties of the attractor.

Return to computing page for the first course APMA0330
Return to computing page for the second course APMA0340
Return to Mathematica tutorial for the first course APMA0330
Return to Mathematica tutorial for the second course APMA0340
Return to the main page for the first course APMA0330
Return to the main page for the second course APMA0340
Return to Part I of the course APMA0340
Introduction to Linear Algebra with Mathematica

The Rössler attractor is the attractor for the Rössler system, a system of three non-linear ordinary differential equations originally studied by the German biochemist Otto Eberhard Rössler (born 20 May 1940).

  1. Rössler, O. E. (1976), "An equation for continuous chaos", Physics Letters A, 57 (5): 397--398.
  2. Rössler, O. E. (1979), "An Equation for Hyperchaos", Physics Letters, 71A (2,3): 155--157.
  3. Rössler, Otto E. (1976). "Chaotic behavior in simple reaction system". Zeitschrift für Naturforschung A. 31: 259–264.
  4. Letellier, C.; V. Messager (2010). "Influences on Otto E. Rössler’s earliest paper on chaos". International Journal of Bifurcation & Chaos. 20 (11): 3585–3616.
Rössler was born in Berlin, into an academic family: his father, also named Otto Rössler, was an Austrian Nazi and a scholar of Semitic languages. Rössler was awarded his MD in 1966. After postdoctoral studies at the Max Planck Institute for Behavioral Physiology in Bavaria, and a visiting appointment at the Center for Theoretical Biology at SUNY-Buffalo, in 1969 he became Professor for Theoretical Biochemistry at the University of Tübingen. In 1994, he became Professor of Chemistry by decree. Rössler has authored hundreds of scientific papers in fields as wide-ranging as biogenesis, the origin of language, differentiable automata, chaotic attractors, endophysics, micro relativity, artificial universes, the hypertext encyclopedia, and world-changing technology.

Rössler and his wife Reimara have been involved with a long-running series of disputes with their employer, the University of Tübingen, which they accuse of discrimination and of violations of academic freedom. In 1988, Reimara Rössler, a professor of medicine, was transferred to a different department within the university; in protest, she began working from home. The state of Baden-Württemberg sued her for failure to perform her assigned duties, as a result of which by 1996 she lost her job and was forced to give up a second home to refund her back pay. Meanwhile, in 1993 and 1994, Otto Rössler had been assigned to teach an introductory chemistry course according to the prescribed curriculum for medical students, but insisted instead on teaching his own material. After he was replaced in the course by another lecturer, he continued trying to give the lectures himself, and was removed by police several times. Because of these incidents, in 1995 a state official tried to force Rössler to undergo psychological tests, but after international protests by many academics this plan was dropped. Rössler continued protesting against his and his wife's treatment by the university and in August 2001 he was caught defacing the university auditorium with spray paint in an attempt to draw attention to his protests.

In 1976, Rössler considered the system of ordinary differential equations

\[ \begin{split} \dot{x} &= -y -z , \\ \dot{y} &= x + a\,y , \\ \dot{z} &= b + z \left( x-c \right) , \end{split} \]
with particular values of parameters a = 0.2, b = 0.2, and c = 5.7. However, the values a = 0.1, b = 0.1, and c = 14 are more commonly used.

These Rössler equations are simpler than those Lorenz used since only one nonlinear term appears (the product xz in the third equation). Otto Rössler solved these equations on a computer. When he made a phase space diagram with the variables, he obtained what has become known as the “Rössler attractor.”

Example 1: Consider the Rössler system of equations with parameter values 𝑎 = 0.32, b = 0.3, and c = 5.8. We plot its solution using the following code:

a = 0.32; b = 0.3; c = 5.8;
sol1 = NDSolve[{x'[t] == -(y[t] + z[t]), y'[t] == x[t] + a* y[t],
z'[t] == b + x[t] z[t] - c* z[t], x[0] == 1, y[0] == 2, z[0] == 3}, {x, y, z}, {t, 0, 200}]
ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. sol1], {t, 0, 200}, WorkingPrecision -> 22, PlotRange -> All, PlotTheme -> "Monochrome"]
However, a similar chaotic behavior can be observed with another values of parameters 𝑎 = 0.1, b = 0.1, and c = 18, and initial conditions x0 = 0, y0 = 1, z0 = 0.
a = 0.1; b = 0.1; c = 18; sol2 = NDSolve[{x'[t] == -(y[t] + z[t]), y'[t] == x[t] + a*y[t], z'[t] == b + x[t] z[t] - c*z[t], x[0] == 0, y[0] == 1, z[0] == 0}, {x, y, z}, {t, 0, 200}] ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. sol2], {t, 0, 200}, WorkingPrecision -> 22, PlotRange -> All, PlotTheme -> "Monochrome"]
     
       Rössler attractor for 𝑎 = 0.32, b = 0.3, and c = 5.8            Rössler attractor for 𝑎 = 0.1, b = 0.1, and c = 18

The maxima will occur at points where the derivative is zero and, except in special cases, they will alternate with minima. You can easily detect where the derivative zero using event detection.
Clear[x, y, z, sol, pts];
{{sol}, {pts}} =
Reap[NDSolve[{x'[t] == -(y[t] + z[t]), y'[t] == x[t] + a y[t],
z'[t] == b + x[t] z[t] - c z[t], x[0] == 1, y[0] == 2, z[0] == 3,
WhenEvent[z'[t] == 0, Sow[t]]}, {x, y, z}, {t, 0, 100}]];
z = z /. sol;
maxPts = Last /@ Partition[pts, 2];
Plot[z[t], {t, 0, 100}, PlotRange -> All, Epilog -> Point[{#, z[#]} & /@ maxPts]]

 

Plane Case

Setting z = 0, allows examination of the behavior on the x,y plane

\[ \begin{split} \dot{x} &= -y -z , \\ \dot{y} &= x + a\,y , \end{split} \]

which is a linear system of ODEs. The stability in the x,y plane can then be found by calculating the eigenvalues of the Jacobian \( \begin{bmatrix} 0& -1 \\ 1&a \end{bmatrix} , \) which are \( \left( a \pm \sqrt{a^2 -4} \right) /2 . \) From this, we can see that when 0 < a < 2, the eigenvalues are complex and both have a positive real component, making the origin unstable with an outwards spiral on the x,y plane. Now consider the z-plane behavior within the context of this range for a. So long as x is smaller than c, the c term will keep the orbit close to the x,y plane. As the orbit approaches x greater than c, the z-values begin to climb. As z climbs, though, the -z in the equation for \( \dot{x} \) stops the growth in x.

 

Fixed points

In order to find the fixed points, the three Rössler equations are set to zero and the (x,y,z) coordinates of each fixed point were determined by solving the resulting equations. This yields the general equations of each of the fixed point coordinates:

\[ x= \frac{c\pm \sqrt{c^2 - 4ab}}{2} , \quad y= - \left( \frac{c\pm \sqrt{c^2 - 4ab}}{2a} \right) , \quad z = \frac{c\pm \sqrt{c^2 - 4ab}}{2a} . \]
This in turn can be used to show the actual fixed points for a given set of parameter values:
\[ \left( \frac{c+ \sqrt{c^2 - 4ab}}{2} , \frac{-c- \sqrt{c^2 - 4ab}}{2a} , \frac{c+ \sqrt{c^2 - 4ab}}{2a} \right) \]
and
\[ \left( \frac{c- \sqrt{c^2 - 4ab}}{2} , \frac{-c+ \sqrt{c^2 - 4ab}}{2a} , \frac{c- \sqrt{c^2 - 4ab}}{2a} \right) . \]
As shown in the general plots of the Rössler Attractor above, one of these fixed points resides in the center of the attractor loop. This fixed point is located in the middle of the attractor and is a saddle-focus with an unstable 2D manifold - an unstable spiral mainly in the x,y plane --- when the trajectory settles down onto a chaotic attractor. The other fixed point is outside of the region of the attractor.

 

Topological description

The topology of the Rössler attractor was first described in terms of a paper-sheet model. Typically, the paper-sheet model can be divided in two stripes, one being a "normal band" and one being a Möbius band. These two bands thus define two different topological domains. In the 1990's, such a topological description led to the concept of a branched manifold, also called knot-holder or template. It can be shown that such a paper-sheet model encodes all topological properties of the unstable periodic orbits embedded within the attractor. The Rössler attractor is the most simple chaotic attractor from the topological point of view, that is, it is a simple stretched and folded ribbon.

 

Eigenvalues and eigenvectors

The stability of each of these fixed points can be analyzed by determining their respective eigenvalues and eigenvectors. Beginning with the Jacobian

\[ \begin{pmatrix} 0&-1&-1 \\ 1&a&0 \\ z&0&x-c \end{pmatrix} . \]
The eigenvalues can be determined by solving the following cubic equation:
\[ \lambda^3 - \lambda^2 \left( a+x -c \right) + \lambda \left( ax+1+z -ac \right) +c-x-az =0 . \]
For the centrally located fixed point, Rössler’s original parameter values of a = 0.2, b = 0.2, and c = 5.7 yield eigenvalues of:
\[ \lambda_{1,2} = 0.0971028 \pm {\bf j}\,0.995786 , \quad \lambda_3 = -5.68718 . \]
The magnitude of a negative eigenvalue characterizes the level of attraction along the corresponding eigenvector. Similarly the magnitude of a positive eigenvalue characterizes the level of repulsion along the corresponding eigenvector.

These eigenveluess have several interesting implications. First, the two complex eigenvalues are responsible for the steady outward slide that occurs in the main disk of the attractor. The last negative eigenvalue is attracting along an axis that runs through the center of the manifold and accounts for the z motion that occurs within the attractor. This effect is roughly demonstrated with the figure below.

 

Poincaré map

The Poincaré map is constructed by plotting the value of the function every time it passes through a set plane in a specific direction. An example would be plotting the y,z value every time it passes through the x = 0 plane where x is changing from negative to positive, commonly done when studying the Lorenz attractor. In the case of the Rössler attractor, the x = 0 plane is uninteresting, as the map always crosses the x = 0 plane at z = 0 due to the nature of the Rössler equations. In the x = 0.1 plane for a = 0.1, b = 0.1, c = 14, the Poincaré map shows the upswing in z values as x increases, as is to be expected due to the upswing and twist section of the Rössler plot. The number of points in this specific Poincaré plot is infinite, but when a different c value is used, the number of points can vary. For example, with a c value of 4, there is only one point on the Poincaré map, because the function yields a periodic orbit of period one, or if the c value is set to 12.8, there would be six points corresponding to a period six orbit.

 

Variation of parameters

Rössler attractor's behavior is largely a factor of the values of its constant parameters a, b, and c. In general, varying each parameter has a comparable effect by causing the system to converge toward a periodic orbit, fixed point, or escape towards infinity, however the specific ranges and behaviors induced vary substantially for each parameter. Periodic orbits, or "unit cycles," of the Rössler system are defined by the number of loops around the central point that occur before the loops series begins to repeat itself.

 

Varying a

Here, b is fixed at 0.2, c is fixed at 5.7 and a changes. Numerical examination of the attractor's behavior over changing a suggests it has a disproportional influence over the attractor's behavior. The results of the analysis are:

 

Varying b

Here, a is fixed at 0.2, c is fixed at 5.7 and b changes. As shown in the accompanying diagram, as b approaches 0 the attractor approaches infinity (note the upswing for very small values of b. Comparative to the other parameters, varying b generates a greater range when period-3 and period-6 orbits will occur. In contrast to a and c, higher values of b converge to period-1, not to a chaotic state.

 

Varying c

Here, a = b = 0.1 and c changes. The bifurcation diagram reveals that low values of c are periodic, but quickly become chaotic as c increases. This pattern repeats itself as c increases --- there are sections of periodicity interspersed with periods of chaos, and the trend is towards higher-period orbits as c increases. For example, the period one orbit only appears for values of c around 4 and is never found again in the bifurcation diagram. The same phenomenon is seen with period three; until c = 12, period three orbits can be found, but thereafter, they do not appear.

 

Simulation

Because three simultaneous equations are involved, with one of them being nonlinear, simulating the Rössler equations is not a simple task unless you are an experienced researcher. Glen Kleinschmidt in Australia has an excellent website ( glensstuff.com ) that deals with a variety of interesting topics. One of the things he discusses is simulating the Rössler equations and generating the Rössler attractor using LTspice.

Rössler hyperchaotic equations

In 1979, Otto Rössler developed a set of four equations that described an autonomous “hyperchaotic” system. The word “hyperchaotic” means that the system of equations must be at least four-dimensional (four variables). In addition, there must be two or more positive Lyapunov coefficients and the sum of all the Lyapunov coefficients must be negative.

The hyperchaotic set of equations that Rössler devised was:

\[ \begin{split} \dot{x} &= -y -z , \\ \dot{y} &= x + 0.25\,y +w, \\ \dot{z} &= 3 +xz , \\ \dot{w} &= 0.05\,w -0.5\,z . \end{split} \]
Glen Kleinschmidt also has simulated Rössler’s hyperchaotic system. These files, too, are available on his website.

It is important to note that when simulating or constructing Chau’s circuit, we are working with an actual oscillator electronic circuit. However, when dealing with the Rössler attractor, what we are simulating is an analog computer circuit which solves the Rössler equations.

 

Return to Mathematica page

Return to the main page (APMA0340)
Return to the Part 1 Matrix Algebra
Return to the Part 2 Linear Systems of Ordinary Differential Equations
Return to the Part 3 Non-linear Systems of Ordinary Differential Equations
Return to the Part 4 Numerical Methods
Return to the Part 5 Fourier Series
Return to the Part 6 Partial Differential Equations
Return to the Part 7 Special Functions