logistic混沌序列延迟时间计算

airubufeng2010

最近在研究混沌，在用互信息发计算logistic混沌时间序列是发现延迟时间并不是理论值1，而是4，我用的程序如下，请问错在哪里，谢谢！
function  tau=Mutual_Information_main(data)
%互信息法求tau
%data;     % 时间序列，列向量
%max_tau = 100;  % 本程序默认最大时延
%Part = 128;     % 本程序默认box大小

[entropy]=mutual(data,128);
for i = 1:length(entropy)-1           
    if (entropy(i)<=entropy(i+1))
        tau = i-1;            
        break;
    end
end
               
plot(0:length(entropy)-1,entropy)
xlabel('tau');
title('互信息法求tau');