c_cmethod_inf.m

CC_method:CC法

function [Smean,Sdeltmean,Scor,tau,tw]=C_CMethod_inf(X,max_d)
% 用于求延迟时间tau
% X为输入时间序列
% max_d为最大时间延迟
% Smean，Sdeltmean,Scor为返回值
% tau为计算得到的延迟时间
% tw为时间窗口
% zhangli
% 2008-11-30


N=length(X);
Smean=zeros(1,max_d);
Scmean=zeros(1,max_d);
Scor=zeros(1,max_d);
delt=std(X);

% 计算Smean,Sdeltmean,Scor
for t=1:max_d
    S=zeros(4,4);
    Sdelt=zeros(1,4);
    for m=2:5
        for j=1:4
            r=delt*j/2;
            Xdt=disjoint(X,N,t);                         % 将时间序列X分解成t个不相交的时间序列
            Xdt=Xdt';
            s=0;
           for tau=1:t
                N_t=floor(N/t);                          % 分成的子序列长度
                Y=Xdt(:,tau);                            % 每个子序列           
                Cs1(tau)=correlation_integral_inf(Y,N_t,r);% 计算C(1,N/t,r,t)            
                Z=reconstitution(Y,N_t,m,1);             % 相空间重构
                Z=Z';
                M=N_t-(m-1); 
                Cs(tau)=correlation_integral_inf(Z,M,r); % 计算C(m,N/t,r,t)
                s=s+(Cs(tau)-Cs1(tau)^m);                % 对t个不相关的时间序列求和
           end            
            S(m-1,j)=s/tau;            
        end
        Sdelt(m-1)=max(S(m-1,:))-min(S(m-1,:));          % 差量计算
    end
    Smean(t)=mean(mean(S));                              % 计算平均值
    Sdeltmean(t)=mean(Sdelt);                            % 计算平均值
    Scor(t)=abs(Smean(t))+Sdeltmean(t);
end

% 寻找时间延迟tau：即Sdeltmean第一个极小值点对应的t
for i=2:length(Sdeltmean)-1
    if Sdeltmean(i)<Sdeltmean(i-1)&Sdeltmean(i)<Sdeltmean(i+1)
        tau=i;
        break;
    end
end
% 寻找时间窗口tw：即Scor最小值对应的t
for i=1:length(Scor)
    if Scor(i)==min(Scor)
        tw=i;
        break;
    end
end