Lyapunov Exponent diagram of one variable out of 3D chaotic map
3 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
I am using the code following code to find the lyapunov exponent of variable x:
clear all; clc
rStart = 3.3;
rStep = 0.001;
rEnd = 3.8;
rValues = rStart:rStep:rEnd;
nPoints = length( rValues );
nIterations = 1000; % number of iterations
LE = zeros( 1, nPoints );
x = zeros( 1, nIterations + 1 );
x(1) = 0.1;
y(1) = 0.002;
z(1) = 0.022;
lamda = 3.79;
beta = 0.01;
alpha = 0.0112;
for k = 1:nPoints
sum = 0;
for i = 1:nIterations
x(i+1)=k*(lamda*x(i)*(1-x(i))+beta*((y(i)).^2)*x(i)+alpha*((z(i)).^3)); % compute and save next x
y(i+1)=(lamda*y(i)*(1-y(i))+beta*((z(i)).^2)*y(i)+alpha*((x(i)).^3)); % compute and save next x
z(i+1)=(lamda*z(i)*(1-z(i))+beta*((z(i)).^2)*z(i)+alpha*((y(i)).^3)); % compute and save next x
sum = sum + log( abs( rValues( k ) - 2*rValues( k )*x( i ) ) );
end
LE( k ) = sum / nIterations;
end
figure; plot( rStart:rStep:rEnd, LE, 'k.-' );
But I get a single dot plot. any help!!!
4 Kommentare
Antworten (0)
Siehe auch
Kategorien
Mehr zu Matrix Computations finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!