logict.m

来自「Chaotic maps. Matlab simulation.」· M 代码 · 共 27 行

%THIS PROGRAM IS WRITTEN FOR IMPLEMENTATION OF CHAOTIC LOGISTIC MAP,ITS
%IMPLEMENTATION,PHASE TRAJECTORY,AUTO CORRELATION.

clc;  close all;  clear all;

A = .5;  B = 4;  phin = -0.35;
phi(1) = ( B*[(A^2) - (phin^2)] ) - A;

for ih = 2:1:6000
  phi(ih) = ( B*[(A^2) - (phi(ih-1)^2)] ) - A;      
end

AST = xcorr(phi,phi);

for tt = 1:1:length(phi)-2    
    XX(tt) =  phi(tt+2);
    YY(tt) =  phi(tt+1);
    ZZ(tt) =  phi(tt);    
end    

figure(1);   plot(phi);             title('\bf LOGISTIC MAP');      xlabel('\bf Time scale');  ylabel('\bf Amplitude fluctuations');
figure(3);   plot(ZZ,YY,'k.');      title('\bf Mapping');           xlabel('\bf X(n)');        ylabel('\bf X(n+1)');
figure(4);   plot(AST);             title('\bf Auto correlation');  xlabel('\bf Time');        ylabel('\bf correlation value');
figure(2);   plot3(ZZ,YY,XX,'r.');  title('\bf Pseudo phase space trajectories');  grid on; 
xlabel('\bf X(n)');  ylabel('\bf X(n+1)');  zlabel('\bf X(n+2)');
%when B = 3.2 P1 when B = 3.5 p2, B = 3.8 CHAOS,3.56 - P3,3.57 More periods