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