function?[ln_rln_C]=G_P（dataNtaumin_mmax_mss）?%?算法2
%?N?是時間序列的長度
%?tau?是固定時間間隔
%?min_m最小的嵌入維數
%?max_m最大的嵌入維數?
%?ssr的步長
data=load（‘E27.txt‘）
N=500
tau=1
min_m=2
max_m=20
ss=15
for?m=min_m?max_m??
????Y=reconstitution（dataNmtau）;%重建矢量空間
max_d=max（max（d））;%是所有點中距離最遠的點
d（11）=max_d;?
??min_d=min（min（d））;%是所有點中距離最近的點??
?delt=（max_d-min_d）*ss;%是r的步長?????
for?k=1ss?
????????r=m