基于Lorenz时间序列求时间延迟，请指教！

xufinance

与陆振波博士工具箱CC算法采用相同的数据，但结果差异很大，附上源代码，帮我看看，不胜感谢！
陆博士：dela_S_mean第一个极小值 时间延迟tao=10

采用相同的Lorenz序列3000个数据，我的结果：dela_S_mean第一个极小值 时间延迟tao=3

matlab源代码如下：
clear all;
data=load('a2.txt') ;
data=data(:,1)';
%************************************************
N=length(data);
max_d=12;%the maximum value of the time delay
sigma=std(data);%calcute standard deviation s_d
  s_m_t=zeros(5,12);% 存放△S(m,t),m取2到5，最大时间延迟为10
  s1=[0 0 0 0 0 0 0 0 0 0 0 0];%存放S(m,t)之和
  
  for t=1:max_d
     for m=2:5
           s_t2=[0 0 0 0];%存放S(m,rj,t) j=1,2,3,4
        for g=1:4
            r=sigma*g/2;
            p=fix(N/t);
            data_d=zeros(t,p);
            for i=1:t
               for j=1:p        
              data_d(i,j)=data(i+(j-1)*t);
               end
             end %将时间序列分解成t个不相交的时间序列      
            N_d=p;
            s_t3=0;
           for i=1:t
                Y=data_d(i,:);
                C_1(i)=correlation_integralxu(Y,N_d,r);%计算C(1,N_d,r,t)
                X=reconstitutionxu(Y,N_d,m,t);%相空间重构
                N_r=N_d-(m-1)*t;
                C_I(i)=correlation_integralxu(X,N_r,r);%计算C(m,N_r,r,t)
                s_t3=s_t3+(C_I(i)-C_1(i)^m); %求S(m,N,r,t)
            end
            s_t2(g)=s_t3/t;
        end
        s1(t)=s1(t)+sum(s_t2)
        s_m_t(m,t)=max(s_t2)-min(s_t2)
    end
end
for i=1:max_d
    dela_s_mean(i)=sum(s_m_t(:,i))/4;
    s_mean(i)=s1(i)/16;
end

figure;
t=1:max_d;
plot(t,dela_s_mean,'*',t,dela_s_mean);title('delta S mean');


其中两个函数：
function C_I=correlation_integral(X,M,r)
%the function is used to calculate correlation integral
sum_H=0;
for i=1:M-1
      for j=i+1:M
        d=norm((X(:,i)-X(:,j)),inf);%calculat the distances of each two points in matris M with sup-norm
         if (r-d)<0
          sita=0;
        else sita=1;
        end      
        sum_H=sum_H+sita;
    end
end
C_I=2*sum_H/(M*(M-1));%the value of correlation integral

function X=reconstitution(data,N,m,tau)
%该函数用来重构相空间
M=N-(m-1)*tau;%相空间中点的个数
K=zeros(m,M);
for j=1:M           %相空间重构
    for i=1:m
        K(i,j)=data((i-1)*tau+j);
    end
end
X=K;

以上代码也是改于网络源代码，请指教，谢谢！！

心境

您这个程序运行不过去啊