Henon map
Generates a 2-dimensional time series using the Henon map.
henon( start = runif(min = -0.5, max = 0.5, n = 2), a = 1.4, b = 0.3, n.sample = 5000, n.transient = 500, do.plot = TRUE )
start |
A 2-dimensional vector indicating the starting values for the x and y Henon coordinates. If the starting point is not specified, it is generated randomly. |
a |
The a parameter. Default: 1.4. |
b |
The b parameter. Default: 0.3. |
n.sample |
Length of the generated time series. Default: 5000 samples. |
n.transient |
Number of transient samples that will be discarded. Default: 500 samples. |
do.plot |
Logical value. If TRUE (default value), a plot of the generated Henon system is shown. |
The Henon map is defined as follows:
x[n] = 1 - a*x[n - 1]^2 + y[n - 1]
y[n] = b*x[n - 1].
The default selection for both a and b parameters (a=1.4 and b=0.3) is known to produce a deterministic chaotic time series.
A list with two vectors named x and y containing the x-components and the y-components of the Henon map, respectively.
Some initial values may lead to an unstable system that will tend to infinity.
Constantino A. Garcia
Strogatz, S.: Nonlinear dynamics and chaos: with applications to physics, biology, chemistry and engineering (Studies in Nonlinearity)
## Not run:
henon.map=henon(n.sample = 1000, n.transient=10,do.plot=TRUE,
start=c(-0.006423277,-0.473545134))
# accessing the x coordinate and plotting it
plot(ts(henon.map$x))
## End(Not run)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.