can anyone please suggest matlab code for lyapunov exponent of logistic map
Have a look at https://au.mathworks.com/matlabcentral/fileexchange/63803-lyapunov-exponent-of-logistic-map
Dear Colleague, thank you for your kind response. But unfortunately, it doesn't work well. Actually, this code doesn't get run in the matlab 2022b.
Please suggest how can I run it.
lyapExp = lyapunovExponent( X , fs ) estimates the Lyapunov exponent of the uniformly sampled time-domain signal X using sampling frequency fs . Use lyapunovExponent to characterize the rate of separation of infinitesimally close trajectories in phase space to distinguish different attractors.
Here is a sample Matlab code for calculating the Lyapunov exponent of the logistic map:
Copy code% Define the parameters of the logistic map r = 3.7; x0 = 0.5; % Define the number of iterations and the step size for calculating the derivative n = 10000; h = 1e-6; % Initialize the Lyapunov exponent lyapunov = 0; % Iterate the logistic map and calculate the Lyapunov exponent for i = 1:n % Calculate the derivative of the logistic map x1 = x0 + h; x2 = x0 - h; f1 = r*x1*(1-x1); f2 = r*x2*(1-x2); df = (f1 - f2)/(2*h); % Update the Lyapunov exponent lyapunov = lyapunov + log(abs(df)); % Update the current value of x x0 = r*x0*(1-x0); end % Normalize the Lyapunov exponent lyapunov = lyapunov/n; % Display the Lyapunov exponent disp(['Lyapunov exponent: ', num2str(lyapunov)]);
This code defines the parameters of the logistic map (r and x0) and the number of iterations and step size for calculating the derivative. It then iterates the logistic map and calculates the Lyapunov exponent by repeatedly updating the current value of x and adding the log of the absolute value of the derivative at each iteration. Finally, the Lyapunov exponent is normalized by dividing by the number of iterations and displayed.
- Badges
- Science topic
- Similar topics
- Mathematics
- Statistics
- Statistical Software
- R