function?multifractal（）
X=xlsread（‘XXX.xls‘）;%在這里輸入數據文件
%?提取數據
[mn]?=?size（sh）;
m=3800;
A?=?X（1:m4）;
format?long?g
L=length（A）;
i=1;
modify=1;
tmin=2;????????????????%?邊框間距，“※”
tmax=10;
ttmin=-10;
ttmax=10;???????????????%?自定義?q?的范圍
for?r=tmin:1:tmax???????????
c（i1）=mod（Lr）;
i=i+1;
end??????????????????
c‘;?????????????????????%?計算不能被邊長r整除的余數
a=L-c‘;?????????????????%?計算并剔除掉不能被邊長r整除的原始數據
n=length（a）;????????????%?求解格網化邊長的個數，即為?n
TT=[];
j=1;
r=tmin;????????????????%?自定義項，“※-2”
for?i=1:1:n?????????????%?即n=25-10+1，自定義的結果
B=A（1:a（i）1）;
U=reshape（Brlength（B）/r）;
T=mean（U）;
T=T‘.*r^3;????????????
TT（1:?length（T）i）=[T];???
modifying（modify1）=length（T）;
r=r+1;
if?r>=?tmin+n???????????%?自定義項，“※-3”
break;?????????????%?邊長超過10+n，超過初始限制，則程序自動終止
end
modify=

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????2903??2013-04-27?14:34??multifractal1.m

-----------?---------??----------?-----??----

?????????????????2903????????????????????1