lorenzattractor: Simulate the Lorenz Attractor

In crqa: Unidimensional and Multidimensional Methods for Recurrence Quantification Analysis

Simulate the Lorenz Attractor

Description

An implementation of the Lorenz dynamical system, which describes the motion of a possible particle, which will neither converge to a steady state, nor diverge to infinity; but rather stay in a bounded but 'chaotically' defined region, i.e., an attractor.

Usage

lorenzattractor(numsteps, dt, sigma, r, b, plots)

Arguments

numsteps	The number of simulated points
dt	System parameter
sigma	System parameter
r	System parameter
b	System parameter
plots	If TRUE, it plots the Lorenz obtained

Value

It returns a matrix with the 3 dimensions of the Lorenz

Author(s)

Moreno I. Coco ([e-mail])

References

Lorenz, Edward Norton (1963). Deterministic nonperiodic flow. Journal of the Atmospheric Sciences 20(2) 130-141.

Examples

## initialize the parameters
numsteps = 2 ^ 11; dt = .01; sigma = 10; r = 28; b = 8/3;
plots = TRUE

res = lorenzattractor(numsteps, dt, sigma, r, b, plots)

crqa documentation built on Nov. 27, 2023, 5:10 p.m.