g_p.m

此程序是用G-P算法计算关联维数Dc。 注意：一定要下载混沌工具箱后使用。

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的步长

for m=min_mmax_m
    Y=reconstitution(data,N,m,tau);%重建矢量空间
    M=N-(m-1)tau;%矢量空间的点数
    for i=1M-1
        for j=i+1M
            d(i,j)=max(abs(Y(,i)-Y(,j)));%计算其余点到点Xi的距离           
        end                              
    end
    max_d=max(max(d));%是所有点中距离最远的点
    d(1,1)=max_d;
    min_d=min(min(d));%是所有点中距离最近的点
    delt=(max_d-min_d)ss;%是r的步长
    for k=1ss
        r=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);%lnr
        fprintf('%d%d%d%dn',k,ss,m,max_m);
    end
    plot(ln_r(m,),ln_C(m,));
    hold on;
end
fid=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);