自相关法求时间延迟tau

%自相关法求时间延迟tau
%根据自相关函数，参考文献，相空间重构
%方法缺陷：不能够运用到高维

function    tau=zi_xiangguan（data）
%data:输入时间序列
%tau:得到的时间延迟
data=[-0.35
-0.15
0.15
0.09
0.17
0.11
0.05
-0.17
-0.66
-0.2
-0.02
-0.5
-0.62
-0.48
-0.32
-0.29
-0.28
-1.17
-0.05
-0.11
-0.31
-0.49
-0.45
0.12
-0.25
-1.08
-0.46
-1.37
-1.45
-0.5
-1.03
-0.43
-0.35
-0.06
-0.48
-0.32];
tau_init=1;
init_xiangguan=0;
for i=1:length（data）-1
        init_xia