%本程序用來計算Hurst指數
%y是輸入的信號
%?h?是計算出的?hurst值
function?h=hurst（y）
m=length（y）;
y0=0;
yend=0;
for?i=1:m
????y0=y（1:i）;
????yend（i）=sum（y0）/i;%每個子區間的均值
end


r=0;
for?i=1:m%對所有子區間進行操作
????x=0;
????for?j=1:i%對每個子區間進行操作
????y0=0;
????y2=0;
????y0=y（1:j）;
????y2=y0-yend（i）;%計算離差
????x（j）=sum（y2）;%累積離差
???end