lyapunov_kantz_c_main.m

混沌kantz算法的mex实现kantz_lyapunov_c

% 产生 Henon 序列　
%x(n+1) = 1 - a * x(n)^2 + y(n);  
%y(n+1) = b * x(n)
clc
clear
close all

a = 1.4;
b = 0.3;

x0 = 0;
y0 = 0;

N1 = 2000;                   % 前面的迭代点数
N2 = 1000;                   % 后面的迭代点数

z = zeros(N1+N2,2);
for i = 1:N1+N2
    x = 1 - a * x0^2 + y0 ;
    y = b * x0;
    x0 = x;
    y0 = y;
    
    z(i,1) = x;
    z(i,2) = y;
end
z = z(N1+1:end,:);

xn = z(:,1);
Y = z(:,2);

% henon 入口参数%
tau = 1;                        % 时延
m = 3;                          % 嵌入维
P = 100;                        % 序列平均周期
fs = 1;                         % 采样频率
evolvement_length=20;           % 最大离散演化步数
norm=2;

figure;
radius=0.0005;
for m=3:5
    data = PhaSpaRecon(xn,tau,m);    % 每列为一个点
    data = data';
    Y = lyapunov_kantz_c(data,tau,m,evolvement_length,norm,P,fs,radius);
    plot(Y,'g.-');
    hold on;
end

radius=0.002;
for m=3:5
    data = PhaSpaRecon(xn,tau,m);    % 每列为一个点
    data = data';
    Y = lyapunov_kantz_c(data,tau,m,evolvement_length,norm,P,fs,radius);
    plot(Y,'r.-');
    hold on;
end

radius=0.008;
for m=3:5
    data = PhaSpaRecon(xn,tau,m);    % 每列为一个点
    data = data';
    Y = lyapunov_kantz_c(data,tau,m,evolvement_length,norm,P,fs,radius);
    plot(Y,'b.-');
    hold on;
end
xlabel('时间延迟tau');
ylabel('距离增长S(tau)');