% 自相關(guān)法求時(shí)延% 使用平臺 - Matlab6.5 / Matlab7.0% 作者：陸振波，海軍工程大學(xué)% 歡迎同行來信交流與合作，更多文章與程序下載請?jiān)L問我的個(gè)人主頁% 電子郵件：luzhenbo@yahoo.com.cn% 個(gè)人主頁：http://luzhenbo.88uu.com.cnclcclearclose all%---------------------------------------------------% 產(chǎn)生 Lorenz 時(shí)間序列sigma = 16;             % Lorenz 方程參數(shù) a = 16 | 10b = 4;                  %                 b = 4 | 8/3r = 45.92;              %                 c = 45.92 | 28        y = [-1,0,1];       % 起始點(diǎn) (1 x 3 的行向量)h = 0.01;           % 積分時(shí)間步長k1 = 30000;          % 前面的迭代點(diǎn)數(shù)k2 = 10000;          % 后面的迭代點(diǎn)數(shù)z = LorenzData(y,h,k1+k2,sigma,r,b);X = z(k1+1:end,1);%---------------------------------------------------% 自相關(guān)法 (直接求 tau)% 自相關(guān)函數(shù)下降到初始值的 1-1/e 時(shí)的 tau 即為所求 (tau 從 1 開始)maxLags = 100;IsPlot = 1;tau_AutoCorrelation = AutoCorrelation(X,maxLags,IsPlot);