关联维数MATLAB程序

计算关联维数的MATLAB程序，文件名是关联维数的首字母，直接可以用。傻瓜式的。花了1个星期的心血啊。

function [prrpcrp]=glws（xmt） 
%函数名为关联维数的首字母用于单串序列，多串到glsw；
%x为要分析的数据；
%x=xlsread（‘d:\matworks\dbin.xls‘）;
[m1n1]=size（x）;
 n=m1;
[mm1mm]=size（m）;
p=zeros（mm2）; %存放拟合系数的矩阵；
rr=zeros（20mm）;%rr是相当于筛子的那个距离，存放的是对数；
cr=zeros（20mm）;%cr是小于筛子距离的距离个数，存放的是对数；
%prr=zeros（20mm）;%rr是相当于筛子的那个距离，存放的是对数；
%pcr=zeros（20mm）;%cr是小于筛子距离的距离个数，存放的是对数；
scope=zeros（191）;
msr=zeros（191）;
   for k=1:mm
      tt=0;
      nm=n-（m（k）-1）*t;%Nm为列数；
      nr=（nm-1）*nm/2;%Nr为距离的总个数；
      juli=zeros（nr1）;%全部距离搞成一列的长矩阵；
      r=zeros（nmnm）;%各列之间距离矩阵；
      y=zeros（m（k）nm）;%重构相矩阵的值yij；
     for j=1:nm
         for i=1:m（k）
            y（ij）=x（j+（i-1）*t）;
         end
     end
     for i=1:nm-1
         for j=i+1:nm
             for kk=1:m（k）
               r（ij）=r（ij）+（y（kkj）-y（kki））^2;
             end
             r（ij）=sqrt（r（ij））; 
             tt=tt+1;
             juli（tt）=r（ij）;
         end
     end
  %进行r和cr个数的计算；
    rmin=min（juli）;
    rmax=max（juli）;
    for i=1:20  %每次把距离间隔分20分来慢慢加；
       rr（ik）=（rmax-rmin）*（i+1）/21; %距离取法值得研究一下；
       for j=1:nr
           if juli（j）<=rr（ik）
                cr（ik）=cr（ik）+1;
           end
       end
       rr（ik）=log（rr（ik））;
       cr（ik）=log（cr（ik）/nr）;
     end
    %rr=rr‘;
    tt=0;
    for i=1:19
        scope（i）=（cr（i+1k）-cr（ik））/（rr（i+1k）-rr（ik））;%每点的斜率；
        tt=tt+scope（i）;
        plot（iscope（i）‘-bd‘）hold on;
    end
    tt=tt/19;%各相邻点间斜率平均值；
    tshold=（max（scope）-min（scope））/2;%threshold，阈值;
    for i=1:19
          msr（i）=abs（scope（i）-tt）;  %各斜率与平均值的均方根，mean square root;
    end
    tt=0;
    for i=2:18
          if （msr（i-1）>tshold & msr（i+1）>tshold）|（msr（i-1）<0.001 & msr（i+1）<0.001）
              continue
          else
              tt=tt+1;
              prr（tt）=rr（ik）;%符合条件的；
              pcr（tt）=cr（ik）;
          end
    end
     p（k1:2）=polyfit（prrpcr1）;%线性拟合，p为两个数，p1为斜率，p2为截距；
  end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2222  2018-09-15 14:46  glws.m
     文件         184  2018-09-15 14:46  README.md