function [tau] = AverageDisplacement(X,maxLags,m)% 平均位移法求混沌時(shí)間序列重構(gòu)的時(shí)間延遲(已知 m,求 tau)% 輸入?yún)?shù)：X         混沌時(shí)間序列%           maxLags　 最大時(shí)間延遲%           m         嵌入維% 輸出參數(shù)：tau　　   時(shí)間延遲　%% 參考文獻(xiàn):呂金虎.混沌時(shí)間序列分析與應(yīng)用.P62%maxLags = maxLags + 1;              % 因?yàn)楹竺嬉鲆粋€(gè)差分，所以這里要加1S_tau = zeros(1,maxLags);for tau = 1:maxLags    xn = PhaSpaRecon(X,tau,m);      % 重構(gòu)相空間    xn_cols = size(xn,2);    temp = zeros(1,xn_cols);    for i = 2:m        temp = temp + (xn(i,:) - xn(1,:)).^2;    end    S_tau(tau) = mean(sqrt(temp));  % tau 所對(duì)應(yīng)的平均位移end% 當(dāng)波形斜率第一次降為初始斜率的 0.4 以下時(shí)的 tau 即為所求 (tau 從 1 開始)slope = diff(S_tau);                % 相鄰 tau 之間的斜率rate = 0.4;gate = slope(1)*rate;temp = find(slope<=gate);if (isempty(temp))    disp('err: max delay time is too small!')    tau = [];else    tau = temp(1);    end