largest_lyapunov_exponent_test.m

来自「计算时间序列的最大的里雅普诺夫指数的函数,」· M 代码 · 共 27 行

    clc
    clear all
    close all    
%     N=1024;
%     fs=1000;     %采样频率
%     f=10;      %信号频率
%     dt=1/fs;   %采样周期
%     t=0:dt:(N-1)*dt; %采样时间间隔
%     ff=0:fs/N:(N-1)*fs/N;
%     data=sin(2*pi*f*t);
    lambda = 3.6;
    x = 0.4*ones(1,length(lambda));
    
    k1 = 4000;                   % 前面的迭代点数
    k2 = 1024;                   % 后面的迭代点数
    
    f = zeros(k1+k2,length(lambda));
    for i = 1:k1+k2
        x = lambda .* x .* (1 - x);
        f(i,:) = x;
    end
    f = f(k1+1:end,:);
    X=f';

    le=largest_lyapunov_exponent(X,1000);