C-C法确定相空间重构参数

function [SmeanSdeltmeanScortautw]=C_CMethod（datamax_d）
% 本函数用于求延迟时间tau和时间窗口tw
% data：输入时间序列
% max_d：最大时间延迟
% Smean，SdeltmeanScor为返回值
% tau：计算得到的延迟时间
% tw：时间窗口
N=length（data）;
%时间序列的长度
Smean=zeros（1max_d）;
%初始化矩阵
Scmean=zeros（1max_d）;
Scor=zeros（1max_d）;
sigma=std（data）;
%计算序列的标准差
% 计算SmeanSdeltmeanScor
for t=1:max_d
    S=zeros（44）;
    Sdelt=zeros（14）;
    for m=2:5
        for j=1:4
            r=sigma*j/2;
            Xdt=disjoint（datat）;
            % 将时间序列data分解成t个不相交的时间序列
            s=0;
           for tau=1:t
                N_t=floor（N/t）;
                % 分成的子序列长度
                Y=Xdt（:tau）;
                % 每个子序列
                %计算C（1N/trt）相当于调用Cs1（tau）=correlation_integral1（Yr）            
                Cs1（tau）=0;