C-C计算延迟时间和嵌入维数4.docx

《C-C计算延迟时间和嵌入维数4.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、functionlambda_1=lyapunov_wolf(data,N?m.tauT)%该函數用来计算时间序列的最大Lyapunov指方法%m:嵌入维数%tau:时间延迟%data:时间序列%N:时间序列长度%P：时间序列的平均周期，选择演化相点距当前点的位虽差，即若当前相点为I,则演化相点只能在卩-JI>P的相点中搜寻%lambda_l:返回晟犬lyapunov指數值min_point=1;%&&要求呆少搜索到的点数MAX_CISHU=5;%&&聂大増加搜索范围次数%FL1NGHAWI<%求最大、最小和平均相点距离%最大相点距离%最

2、小相点距离max_d=0;min_d=1.0c+10();%相空间重构%重构相空间中柏点的个数avg_dd=0;Y=reconstitution(data,Njn^au);M=N-(m-l)*tau;fi)ri=1:(NI-1)forj=i+1:Md=0;f()rk=1:md=d+(Y(k3)-Y(k,j))*^(k3)-Y(k,j))；精品d=sqrt(d);ifmax_ddmin_d=d;endavg_dd=avg_dd+&end%平均相点距离enddlt_cps=(avg_d-min_d)*0

3、.02;时，对max_cps的放宽幅度min_cps=min_d+dlt_cps/2;max_cps=milled+2*dlt_cps;avg_d=2*avg_dd/1));%若在min.cps—max_cps中找不到演牝相点%演化相点与当前相点距离的最小限%&&演化相点与当前相点距离的最大限从P+1-M-1个相点中找与第一个相点最近的相点位虽(Loc.DIC)71其晟短距离DKDK=1.0c+100;Loc.DK=2;%第i个相点到其最近距离点的距离%第i个相点对应的最近距离点的下精品fori=(P+l)：(M-l)%限制短暂分离，从点P

4、+1开始搜索d=0;f()rk=1:md=d+^(k,i)-Y(k,l))*(Y(k,i)-Y(k,l));endd=sqrt(d);if(d

5、a«sumJmd=0;%存放前i个log2(DK1/PK)的累计和fori=2:(Nkl)%计算演化距离DK1=0;fork=1:mOKI=OKI+(Y(k,i)・Y(k,Loc_DK+l))*(Y(k,i)・Y(k,Loc_DK+l));endOKI=sqrt(DKl);old_Loc_DK=Loc_DK;%保存原最近位宜相点old_PK=PI<;%计算前i个h〉g2(DK1/DK)的累计和以尺保存i点的李氏指数if(DKl〜=0)&(DK〜=0)sumjmd=sum」md+log(DI

6、=sumjmd/(i-l);%以下寻找i点的最短距离：要求距离在指定距离范围内尽量短，与DK1的角度股小point_num=0;%&&在指定距离范围内找到的侯选相点的个数cos_siia=0;%&&夹角余弦的比较初值——要求一定爱锐角zjfwcs=o；%&&增加范围次数while(p()int_num==0)%*搜索相点fbrj=1:(M-1)ifabsQ-i)<=(P-1)%&&候选点距当前点太近，姚过！continue;end%来计算候选点与当前点的距离dncw=0;精品f()rk=1:mdncxv=dncw+(Y(k3)-Y(k,j)

7、)*^(k,i)-Y(k,j))；enddnrw=s屮t(dncyv);if(dnewVmin_cps)

8、(dnew>max_cps)%&&不在距离范围，跳过！continue;end%*计算夹角余弦氏比较DOT=0;f()rk=1:mDOT=D()T+(Y(k,i)・Y(kj))#(Y(hi)・Y(k,old_Loc_DK+l));endCTH=D()T/(dn如DK1);ifacos(CTH)>(3.14151926/4)%&&不畏小于45度的角，跳过！continue;endifCTH>cos.sita%&&新夹角小于过去巳找到的相点

9、的夹角，保留cos_sita=CTH;Loc_DK=j;DK=dnew;endpoinenum=point_num+1;endifp()int_num<=min_poinimax