Hello,
I am trying to use the delay_fnn function in the DynamicalSystems package and I keep getting the following error:
ERROR: MethodError: no method matching StateSpaceSet{2, Float64, SVector{2, Float64}, Nothing}(::Vector{SVector{2, Float64}})
The type StateSpaceSet{2, Float64, SVector{2, Float64}, Nothing} exists, but no method is defined for this combination of argument types when trying to construct it.
Here is the code:
'function lorenz96_rule!(du, u, p, t)
F = p[1]; N = length(u)
# 3 edge cases
du[1] = (u[2] - u[N - 1]) * u[N] - u[1] + F
du[2] = (u[3] - u[N]) * u[1] - u[2] + F
du[N] = (u[1] - u[N - 2]) * u[N - 1] - u[N] + F
# then the general case
for n in 3:(N - 1)
du[n] = (u[n + 1] - u[n - 2]) * u[n - 1] - u[n] + F
end
return nothing # always return nothing for in-place form!
end
N = 6
u0 = range(0.1, 1; length = N)
p0 = [8.0]
lorenz96 = CoupledODEs(lorenz96_rule!, u0, p0)
total_time = 12.5
sampling_time = 0.02
Y, t = trajectory(lorenz96, total_time; Ttr = 2.2, Δt = sampling_time)
z = Y[:, 1]
τ = estimate_delay(z, “mi_min”)
Ed = delay_fnn(z, τ, 2:8)’
Thanks for any help. New to Julia (from Matlab) but familiar with attractor reconstruction.