I'm not quite understand the workflow of this file. Could you explain it step by step? And, what are the value of parameter meanperiod and maxiter? pls reply via my email: firstname.lastname@example.org
Thanks for your help~
hello, am trying to calculate the lyapunov exponents of some meteorological data(daily averages of temp, rainfall etc)collected over 34yrs. i have 365 data points in a yr, is that sufficient? my sample interval is 1day(24hrs) can this work? what do i do? pls reply here and via my emai: email@example.com.
thanks for your efforts.
I have tested my algortihm with well known time series data (for example Lorenz) and i have found LLE for Lorenz data as 1.5013. I used <ln(d)> for the divergence value. Did you try the data with another code and find LLE as 1.5? You can send the values of your parameters and your data via e mail.
I am trying to repeat Rosenstein's results on well known attractors with your lyarosenstein() function. Everything works fine, but the results does not meet the Rosenstein's ones. I construct e.g. Lorenz attractor and obtain data from it with the exact same parameters as Rosenstein did. The plateau and the slope of curve is clearly visible, but the slope is always little bigger than it should (not 1.50 for Lorenz, but e.g. 1.80 etc.). Have you tested your function somehow on those results or do you have idea what might be wrong? I feel there can be a problem with y-axis scope. He uses <ln(d)> and I'm not sure what he means by it.
Many thanks in advance for any help. Its great function anyway.