function [tau] = DelayTime_ComplexAutoCorrelation(X,maxLags,m,IsPlot)% (去偏)復自相關法求混沌時間序列重構的時間延遲(已知 m,求 tau)% 輸入參數：X         混沌時間序列%           maxLags　 最大時間延遲%           m         嵌入維% 輸出參數：tau　　   時間延遲　%% 參考文獻:呂金虎.混沌時間序列分析與應用.P63%X_mean = mean(X);C_tau = zeros(1,maxLags);for tau_i = 1:maxLags    xn = PhaSpaRecon(X,tau_i,m);      % 重構相空間    xn_cols = size(xn,2);    temp = zeros(1,xn_cols);    for i = 2:m        temp = temp + (xn(1,:)-X_mean).*(xn(i,:)-X_mean);    end    C_tau(tau_i) = mean(temp);  % tau 所對應的平均位移end% 去偏復自相關函數下降到初始值的 1-1/e 時的 tau 即為所求 (tau 從 1 開始)gate = (1-exp(-1))*C_tau(1);temp = find(C_tau<=gate);if (isempty(temp))    disp('err: max delay time is too small!')    tau = [];else    tau = temp(1);    endif IsPlot    figure;    plot(C_tau)    xlabel('Lag');    title('(去偏)復自相關法');end