混沌时间序列c-c法求时间延迟

clc
clear all
A=load（‘F:/chaos/monthly rainfall number1.txt‘）; %使用相应数据的文件名。
x=A‘;%转化为行向量
%t=input（‘输入时间序列所要分成不相交的子序列的个数‘）;
for t=1:20
    %t1=input（请输入子序列时间延迟ti=:‘）
    t1=1;
    n1=length（x）;
    p=n1/t
    %m0=input（请输入嵌入维数m=;‘）
    St1=0;
    St2=0;
    x1=[];
    s=[];
    for m0=2:5
        for n=1:4
            for i=1:t
                x1（i:）=x（1i:t:n1-t+i）;%将时间序列x分成t个不相交的子序列
                                        %对每个子序列重构相空间的点
                 y=[];
                 for j=1:p-（m0-1）*t1
                     y（j:）=x1（ij:t1:j+（m0-1）*t1）;
                 end
                 %计算各嵌入时间子序列的关联积分cs（m.Ntlt）以及关联维数
                 r=std（x）*n/2;
                 [hl]=size（y）;
                 Cr=0.01;
                 for k=1:h
                     for j=k+1:h
                         d（kj）=r-norm（y（k:）-y（j:））;
                         if d（kj）>0
                             Cr=Cr+1;
                         end
                     end
                 end
                Cr1（i）=