Lorenz96

Systems · Continuous · Chaotic attractors

Lorenz-96 model on a 1D ring of N weakly coupled scalar variables.

Dimension: variable

Equations

@staticmethod
def _equations(y_sym, t_sym, *, f, N):
    return [
        (y_sym((i + 1) % N) - y_sym((i - 2) % N)) * y_sym((i - 1) % N) - y_sym(i) + f
        for i in range(N)
    ]

Parameters

parameter	default
f	8.0
N	20

Usage

import tsdynamics as ts

sys = ts.Lorenz96()
traj = sys.integrate(final_time=100.0, dt=0.01)

Back to Chaotic attractors