MF_DFA求hurst指数

function F_n=DFA（DATAwin_lengthorder）
      
       N=length（DATA）;   
       n=floor（N/win_length）;
       N1=n*win_length;
       y=zeros（N11）;
       Yn=zeros（N11）;
       
       fitcoef=zeros（norder+1）;
            mean1=mean（DATA（1:N1））;
       for i=1:N1
           y（i）=sum（DATA（1:i）-mean1）;
       end