function?multifractal（）
X=xlsread（‘XXX.xls‘）;%在這里輸入數(shù)據(jù)文件
%?提取數(shù)據(jù)
[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‘;?????????????????????%?計(jì)算不能被邊長(zhǎng)r整除的余數(shù)
a=L-c‘;?????????????????%?計(jì)算并剔除掉不能被邊長(zhǎng)r整除的原始數(shù)據(jù)
n=length（a）;????????????%?求解格網(wǎng)化邊長(zhǎng)的個(gè)數(shù)，即為?n
TT=[];
j=1;
r=tmin;????????????????%?自定義項(xiàng)，“※-2”
for?i=1:1:n?????????????%?即n=25-10+1，自定義的結(jié)果
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???????????%?自定義項(xiàng)，“※-3”
break;?????????????%?邊長(zhǎng)超過(guò)10+n，超過(guò)初始限制，則程序自動(dòng)終止
end
modify=

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

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

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

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