MFDFA 多重分形程序 MATLAB

一个计算信号多重分形谱的程序，希望对大家有用

function[HHintlogFlogS]=MFDFA（tsorderq）
ts_dim=size（ts）;
if min（ts_dim）>1
    disp（‘输入数据错误！‘）;
    return;
end
if ts_dim（2）>1
    ts=ts‘;
end
ts_m=mean（ts）;
y=cumsum（ts-ts_m）;
N=length（y）;
F=[];
S=[];
for s=5:1:（N/5）
    Fsjj=[];
    Ns=floor（N/s）;
    for j=1:1:（2*Ns）
        if j<=Ns
           t=y（（（j-1）*s+1）:（j*s））;
        else
            t=y（（N-（j-Ns）*s+1）:（N-（j-Ns-1）*s））;
        end
        fitcoef=polyfit（（1:s）‘torder）;
        fitvalue=polyval（fitcoef（1:s）‘）;
        reobject=t-fitvalue;
        Fsj=mean（reobject.^2）;
        Fsjj=[Fsjj;Fsj];
    end
    Fqs=（mean（Fsjj.^（q/2）））^（1/q）;
    F=[F;Fqs];
    S=[S;s];
end
logF=log10（F）;
logS=log10（S）;

len_logS=length（logS）;
X=[ones（len_logS1） logS];
y=logF;
[bbint]=regress（yX0.05）;
H=b（2）;
Hint=bint（2:）;
end

    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        846  2013-07-24 12:52  MFDFA.m

----------- ---------  ---------- -----  ----

                  846                    1