MF-DFA计算广义Hurst指数

k0=[];
q0=[];
tq=[];
dq=[];
rf=[];
a=2;
frf=[];
p1=xlsread('C:\Users\子彦\Desktop\毕业论文\正文\data\代码\收益序列（周）.xlsx');
y=p1';
for q=-10:a:10
        Fq=[];
        s=[];
        num=[];
        den=[];
        for i=3:length(y)/5       %i=3:118.6
            n=floor(length(y)/i);     %每个区间的长度，向下取整
            err1=[];
            err2=[];
            errq=[];
            errq2=[];
            ss=[];
            ss2=[];
            for j=1:n
                temp=y((1+(j-1)*i):j*i);
                x=1:length(temp);
                p=polyfit(x,temp,1);    %用多项式拟合x,temp的表达式
                yy=polyval(p,x);      %返回n次多项式x处的值
                e=temp-yy;
                er=0;
                for k=1:length(e)