ThomasLabyrinth¶

Systems / ODEs / Chaotic attractors

Thomas's cyclically symmetric flow at weak damping — a deterministic random walk (labyrinth chaos) across a 3-D lattice of unstable equilibria.

continuous · ODE3 dimensionschaotic

Interactive: drag to rotate

Definition¶

\[ \begin{aligned} \dot{x} &= - a x + b \sin{\left(y \right)} \\ \dot{y} &= - a y + b \sin{\left(z \right)} \\ \dot{z} &= - a z + b \sin{\left(x \right)} \end{aligned} \]

Parameters¶

Symbol	Default	Role
`a`	`0.5`	dissipation (friction) coefficient
`b`	`10`	sinusoidal forcing amplitude

State variables: x, y, z

Properties¶

Lyapunov spectrum

TODO — lyapunov_spectrum failed: AttributeError

Kaplan–Yorke dimension

TODO — requires a numeric Lyapunov spectrum

Divergence ∇·f

$\nabla\!\cdot f = - 3 a$
constant

Equilibria

194 equilibria
0 stable · 194 unstable

Use it in TSDynamics¶

import tsdynamics as ts

system = ts.systems.ThomasLabyrinth()
traj = system.run(final_time=100.0, dt=0.01)
ts.plot(traj)

spectrum = ts.analysis.lyapunov_spectrum(system)
ts.analysis.kaplan_yorke_dimension(spectrum)

Reference¶

Thomas (1999), Int. J. Bifurc. Chaos 9, 1889-1905
doi:10.1142/s0218127499001383

BibTeX

@misc{thomaslabyrinth,
  title = {ThomasLabyrinth system},
  note = {Thomas (1999), Int. J. Bifurc. Chaos 9, 1889-1905},
  doi = {10.1142/s0218127499001383}
}