HyperLorenz¶
Systems / ODEs / Chaotic attractors
A four-dimensional Lorenz extension with a feedback state that lifts the flow into hyperchaos — two positive Lyapunov exponents.
Interactive: drag to rotate
Definition¶
\[
\begin{aligned}
\dot{x} &= - a x + a y + w \\
\dot{y} &= c x - x z - y \\
\dot{z} &= - b z + x y \\
\dot{w} &= d w - x z
\end{aligned}
\]
Parameters¶
| Symbol | Default | Role |
|---|---|---|
a |
10 |
Lorenz-core coefficient |
b |
2.667 |
Lorenz-core coefficient |
c |
28 |
Lorenz-core coefficient |
d |
1.1 |
fourth-variable feedback coefficient |
State variables: x, y, z, w
Properties¶
Lyapunov spectrum
TODO — lyapunov_spectrum failed: AttributeError
Kaplan–Yorke dimension
TODO — requires a numeric Lyapunov spectrum
Divergence ∇·f
$\nabla\!\cdot f = - a - b + d - 1$
constant
constant
Equilibria
1 equilibria
0 stable · 1 unstable
0 stable · 1 unstable
Use it in TSDynamics¶
import tsdynamics as ts
system = ts.systems.HyperLorenz()
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¶
Meier (2003), Presentation of Attractors with Cinema