lorenz_system {ARGOS}

Lorenz Chaotic System

Description

Simulates the Lorenz chaotic system with optional noise.

Usage

lorenz_system(n, init_conditions, dt, snr = Inf)

Arguments

n

Number of time points (rounded to the nearest integer).

init_conditions

Initial conditions as a numeric vector of length 3 (X, Y, Z).

dt

Time step between observations.

snr

Signal-to-noise ratio (in dB). Use Inf for no noise.

Details

This function simulates the Lorenz chaotic system with the specified parameters. It uses the specified time step and initial conditions to compute the system's state over time. If a non-Infinite SNR is provided, Gaussian noise is added to the system.

Value

A numeric matrix representing the system's state over time. Each row corresponds to a time point, and each column represents a variable (X, Y, Z).

Examples

# Simulate the Lorenz system with 1000 time points and no noise
data <- lorenz_system(
  n = 1000,
  dt = 0.01,
  init_conditions = c(-8, 7, 27),
  snr = Inf
)