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）=