%這個方法是利用自相關系數法求解最佳延遲時間
function?Tau=autocorrelation（datatau_max）
%data:輸入的時間序列
%tau_max:最大時間延遲
%Tau:返回所求時間序列的時間延遲
N=length（data）;
%時間序列長度
x_mean=mean（data）;
%時間序列的平均值
data=data-x_mean;
SSd=dot（datadata）;
R_xx=zeros（1tau_max）;
%自相關函數初始化
for?tau=1:tau_max
%計算自相關函數
????for?ii=1:N-tau
????????R_xx（tau）=R_xx（tau）+data（ii）*data（ii+tau）;
????end
????R_xx（tau）=R_xx（tau）/SSd;
end
plot（1:tau_maxR_xx）;
hold?on
%作自相關函數圖
line（[1?tau_max][0?0]）
title（‘自相關法求時間延遲‘）;
ylabel（‘自相關函數‘）;
xlabel（‘時間延遲‘）;
Tau=0;
for?jj=2:tau_max
%求時間序列的時間延遲??
????if?R_xx（jj-1）*R_xx（jj）<=0
???????if?abs（R_xx（jj-1））>abs（R_xx（jj））
???????????Tau=jj;break
???????else
???????????Tau=jj-1;break
???????end
????end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????265429??2012-08-22?22:24??MackeyGlass_t17.txt
?????文件????20323296??2004-10-13?09:14??X02
?????文件????20323296??2004-10-13?09:14??X03
?????文件????20323296??2004-10-13?09:14??X01
?????文件????????6368??2019-11-25?23:28??echo_s_n.m
?????文件?????????833??2019-04-15?16:53??autocorrelation.m