CC方法对一维时间序列的相空间重构

appletree_hello

CC代码是可以对一维的时间序列进行相空间重构吗？ 在下很急 有偿帮助也接受万分感谢

Vickyvictoria

1. function [xn,dn] = PhaSpaRecon(s,tau,m)
   
2. % 混沌序列的相空间重构 (phase space reconstruction)
   
3. % [xn,dn] = PhaSpaRecon(s,tau,m)
   
4. % 输入参数：    s          混沌序列
   
5. %               tau        重构时延
   
6. %               m          重构维数
   
7. % 输出参数：    xn         相空间中的点序列(每一列为相空间中一个点)
   
8. %               dn         一步预测的目标
   
9. 
   
10. len = length(s);
    
11. if (len-1-(m-1)*tau < 1)
    
12.     disp('err: delay time or the embedding dimension is too large!')
    
13.     xn = [];
    
14.     dn = [];
    
15. else
    
16.     xn = zeros(m,len-1-(m-1)*tau);
    
17.     for i = 1:m
    
18.         xn(i,:) = s(1+(i-1)*tau : len-1-(m-i)*tau);   % 相空间重构，每一列为一个点
    
19.     end
    
20.     dn = s(2+(m-1)*tau : end);    % 预测的目标
    
21. end

复制代码