hurst指数 DFA matlab程序

%Estimation of Hurst parameter
clear all;
Name=‘hb2‘
x=xlsread（strcat（Name‘.xls‘））;  %reand fault recorder
nX = length（x）; 
AV=sum（x）/nX;
y=zeros（nX1）;
for i=1:nX
    y（i）=sum（x（1:i））-i*AV;
end
%xlswrite（‘Yi‘y）
Result=zeros（122）;
FF=zeros（121）;
for kk=24:-2:2       % divided into multispan
    yy=zeros（nX/kkkk）;
    yyy=zeros（nX/kk1）;
    x=1:1:nX/kk;
    for i=1:kk   % check data in each block
        yy（:i）=y（（i-1）*nX/kk+1:i*nX/kk）;   %set data within a span
        p = polyfit（x‘yy（:i）2）;
        f= polyval（px）;
        FF（kk/2）=FF（kk/2）+ sum（（yy（: i） - f（:）） .^ 2）;
    end

    a=mod（nX kk）;      % use information of left data
    if （a）
        yy=[yyyyy];