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‘;?????????????????????%?計算不能被邊長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=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?為任意數，這里取1到n，為n，與?k取值保持一致，q過大，計算機無法
??????????????????????????????????%識別，默認為無窮大，q過小，結果接近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];?????????????%?這里用到兩個循環，即考慮到了冪函數，又需考慮求和
end
j=j+1;
end
XX;?????????????????????%?得到質量分配函數，Xq（ξ），????????????????
X=log（tmin:1:tmax）;
%?X=log（tmax:-1:tmin）;??????????%?此系以前的自定義輸入結果，“※—4”
Y=log（XX）;
%?figure（1）
%?plot（X‘Y‘o-k‘）????????%?至此，計算多重分形譜的第一步，分配函數構建完畢
side_length=?tmin:1:tmax;?????????%?自定義網格邊長，“※—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）;??%?在對數尺度下計算斜率
slope（m1）=b（11）;???????????????????????
m=m+1;
end
slope;?????????????????????????????????%??這里的Slope即為質量指數，τ（q）?
N=polyfit（q‘?slope1）;?
plot（q‘?slope）????%??此步是考察τ（q）-q?之間的關系，
a=diff（slope）./diff（q‘）;???%?第三步計算，diff函數求偏導確實少一列
q=q‘;
f_a=a.*q（1:end-1）-slope（1:end-1）;
%?figure（2）
%?a=sort（a‘ascend‘）;
plot（af_a‘o-k‘）
xlabel（‘α‘‘FontSize‘12）;
ylabel（‘f（α）‘‘FontSize‘12）;
%?polyfit（af_a3）
a=sort（a‘ascend‘）;
%?a+1
%?f_a+1
a????%?奇異性指數

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