% C-C方法重構相空間主文件% 使用平臺 - Matlab6.5 / Matlab7.0% 作者：陸振波，海軍工程大學% 歡迎同行來信交流與合作，更多文章與程序下載請訪問我的個人主頁% 電子郵件：luzhenbo@yahoo.com.cn% 個人主頁：http://luzhenbo.88uu.com.cnclcclearclose all%---------------------------------------------------% 產生 Lorenz 時間序列sigma = 16;          % Lorenz 方程參數 ab = 4;               %                 br = 45.92;           %                 c            y = [1,0,0];        % 起始點 (1 x 3 的行向量)h = 0.01;            % 積分時間步長k1 = 2000;           % 前面的迭代點數k2 = 3000            % 后面的迭代點數z = LorenzData(y,h,k1+k2,sigma,r,b);z = z(k1+1:end,:);X = z(:,1);maxLags = 200;%---------------------------------------------------tic[S_mean,delta_S_mean,S_cor] = CC_luzhenbo(X,maxLags);toc%---------------------------------------------------figure    subplot(311)plot(1:maxLags,S_mean); grid; title('S mean')subplot(312)plot(1:maxLags,delta_S_mean); grid; title('delta S mean')subplot(313)plot(1:maxLags,S_cor); grid; title('S cor')