G-P算法计算关联维数.doc

《G-P算法计算关联维数.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、G-P算法计算关联维数CODE:function[ln_r,ln_C]=G_P(data,N,tau,min_m,max_m,ss)%此程序是用G-P算法计算关联维数Dc%N是时间序列的长度%tau是固定时间间隔%min_m最小的嵌入维数%max_m最大的嵌入维数%ssr的步长form=min_mmax_mY=reconstitution(data,N,m,tau);%重建矢量空间M=N-(m-1)tau;%矢量空间的点数fori=1M-1forj=i+1Md(i,j)=max(abs(Y(,i)-Y(

2、,j)));%计算其余点到点Xi的距离endendmax_d=max(max(d));%是所有点中距离最远的点d(1,1)=max_d;min_d=min(min(d));%是所有点中距离最近的点delt=(max_d-min_d)ss;%是r的步长fork=1ssr=min_d+kdelt;C(k)=correlation_integral(Y,M,r);%计算关联积分函数ln_C(m,k)=log(C(k));%lnC(r)ln_r(m,k)=log(r);%lnrfprintf('%d%d%d%d

3、n',k,ss,m,max_m);endplot(ln_r(m,),ln_C(m,));holdon;endfid=fopen('lnr.txt','w');fprintf(fid,'%6.2f%6.2fn',ln_r);fclose(fid);fid=fopen('lnC.txt','w');fprintf(fid,'%6.2f%6.2fn',ln_C);fclose(fid);