clc
clear?all
A=load（‘3333333.txt‘）;?%使用相應(yīng)數(shù)據(jù)的文件名。
x=A‘;%轉(zhuǎn)化為行向量
%t=input（‘輸入時(shí)間序列所要分成不相交的子序列的個(gè)數(shù)‘）;
for?t=1:20
????%t1=input（請(qǐng)輸入子序列時(shí)間延遲ti=:‘）
????t1=1;
????n1=length（x）;
????p=n1/t
????%m0=input（請(qǐng)輸入嵌入維數(shù)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）;%將時(shí)間序列x分成t個(gè)不相交的子序列
????????????????????????????????????????%對(duì)每個(gè)子序列重構(gòu)相空間的點(diǎn)
?????????????????y=[];
?????????????????for?j=1:p-（m0-1）*t1
?????????????????????y（j:）=x1（ij:t1:j+（m0-1）*t1）;
?????????????????end
?????????????????%計(jì)算各嵌入時(shí)間子序列的關(guān)聯(lián)積分cs（m.Ntlt）以及關(guān)聯(lián)維數(shù)
?????????????????r=std（x）*n/2;
?????????????????[hl]=size（y）;
?????????????????Cr=0.01;
?????????????????for?k=1:h
?????????????????????for?j=k+1:h
???????????????????