A=imread（‘cameraman.tif‘）;
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)超過10+n，超過初始限制，則程序自動(dòng)終止
end
modify=modify+1;
end???????????????
modifying;
TT=?nthroot（TT1）;????????
%?或者不縮小
%?TT;
TT=nonzeros（TT）;??????????
for?cugb=1:1:n?????????????????
modifying_modifying（cugb1）=sum（modifying（1:cugb））;?????%?有影響的新加卷?%?&?*?%?￥?#?@??。?……
end????????????????????????????????????????
modifying_modifying;?????????????
j=1;
???????????????????????%?q?為任意數(shù)，這里取1到n，為n，與?k取值保持一致，q過大，計(jì)算機(jī)無(wú)法
??????????????????????????????????%識(shí)別，默認(rèn)為無(wú)窮大，q過小，結(jié)果接近0，則意義不明確
for?q=ttmin:ttmax???????????%這里取?q=-10:1:10
for?k=1:1:n?
????X=TT（1:modifying_modifying（k1）1）.^q;??????????????
????if?k>1?????????????????
??????X=TT（modifying_modifying（k-11）+1:modifying_modifying（k1）1）.^q;??
????end
t=sum（X）;
XX（kj）=[t];?????????????%?這里用到兩個(gè)循環(huán)，即考慮到了冪函數(shù)，又需考慮求和
end
j=j+1;
end
XX;?????????????????????%?得到質(zhì)量分配函數(shù)，Xq（ξ），????????????????
X=log（tmin:1:tmax）;
%?X=log（tmax:-1:tmin）;??????????%?此系以前的自定義輸入結(jié)果，“※—4”
Y=log（XX）;
%?figure（1）
%?plot（X‘Y‘o-k‘）????????%?至此，計(jì)算多重分形譜的第一步，分配函數(shù)構(gòu)建完畢
side_length=?tmin:1:tmax;?????????%?自定義網(wǎng)格邊長(zhǎng)，“※—5”
side_length=side_length‘;
q=ttmin:ttmax;????????????????????%?q=-5:1:n-10，q=-10:1:10
m=1;
[hahb]=size（XX）;
for?i=1:1:hb
%?XX=XX’;?????????????????????????%?or?not
s=XX（:i）;????????????????????????%?XX
b=polyfit（log（side_length）log（s）1）;??%?在對(duì)數(shù)尺度下計(jì)算斜率
slope（m1）=b（11）;???????????????????????
m=m+1;
end
slope;?????????????????????????????????%??這里的Slope即為質(zhì)量指數(shù)，τ（q）?
N=polyfit（q‘?slope1）;?
plot（q‘?slope）????%??此步是考察τ（q）-q?之間的關(guān)系，
a=diff（slope）./diff（q‘）;???%?第三步計(jì)算，diff函數(shù)求偏導(dǎo)確實(shí)少一列
q=q‘;
f_a=a.*q（1:end-1）-slope（1:end-1）;
%?figure（2）
%?a=sort（a‘a(chǎn)scend‘）;
plot（af_a‘o-k‘）
xlabel（‘α‘‘FontSize‘12）;
ylabel（‘f（α）‘‘FontSize‘12）;
%?polyfit（af_a3）
a=sort（a‘a(chǎn)scend‘）;
%?a+1
%?f_a+1
a????%?奇異性指數(shù)

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-02-15?16:12??duochongfenxing\
?????文件????????2772??2010-12-29?19:40??duochongfenxing\多重分形譜算法_1.m
?????文件????????1223??2010-12-29?19:45??duochongfenxing\盒子分形維數(shù).m