%自相關法求時間延遲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