%?clc
clear
a=imread（‘C:\0.28%red\N0-red-10X-5500.jpg‘）;
a=rgb2gray（a）;
??x=double（a）;
%???x1=x（65:19265:192）
x1=x/sum（sum（x））;
for?q=-60:2:60
????w=（q+60）/2+1;

????for?k=1:8
????????xe（k）=log（2^k/128）;
????????fun=inline（‘sum（sum（Y））‘）;
?????????b=blkproc（x1[2^k2^k]fun）;
?????????a1=b.^q;
????????b1=sum（sum（a1））;
??????????x12=b;
???[row?vic]=find（x12==0）;
???ss=size（row）;
???for?i=1:ss
????x12（row（i）vic（i））=1;
???end
?????????a2=log（x12）;

?????????a3=a1.*a2;????
?????????b3=sum（sum（a3））;
?????????rq（k）=b3/b1;
?????????fq（k）=b3/b1*q-log（b1）;
????end
????Sxx=sum（xe.^2）-（sum（xe））^2/8;

?????Sxr=xe*rq‘-sum（xe）*sum（rq）/8;
?????Sxf=xe*fq‘-sum（xe）*sum（fq）/8;

?????r（w）=Sxr/Sxx;
?????f（w）=Sxf/Sxx;

???end
???q1=-60:2:60;
%?figureplot（rf）;
q=-60:2:60;
%?figureplot（qr）;%%%%%%%%%%%%%ｑ和阿爾法a
%?figureplot（qf）;%%%%%%%%%%%%%ｑ和f（a）

%%%%%%%%%%%%%%%f形譜寬ｄｅｌｔａa?譜差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形譜寬ｄｅｌｔａa?譜差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形譜寬ｄｅｌｔａa?譜差ｄｅｌｔａf
rmax=max（r）;
rmin=min（r）;
deltar=（rmax-rmin）

%?fmax=max（f）;
%?fmin=min（f）;
%?（fmax-fmin）

ffmin=find（r==rmin）;
fmin=f（ffmin）;

ffmax=find（r==rmax）;
fmax=f（ffmax）;

deltaf=fmin-fmax%%%%%%這個(gè)才是真正的譜差啊啊啊啊啊啊啊
????
%??????amin=find（isnan（hurstH））;
%??????hurstH（hursti）=[];

%%%%%%%%%%%%%%%f形譜寬ｄｅｌｔａa?譜差ｄｅｌｔａf%%%%%%%%%%%%%%%f形譜寬ｄｅｌｔａa?譜差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形譜寬ｄｅｌｔａa?譜差ｄｅｌｔａf%%%%%%%%%%%%%%%f形譜寬ｄｅｌｔａa?譜差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形譜寬ｄｅｌｔａa?譜差ｄｅｌｔａf
figureplot（rf）;
xlabel（‘奇異指數(shù)a‘）ylabel（‘多重分行譜f（a）‘）
%?q=-60:2:60;
%?figureplot（qr）;%%%%%%%%%%%%%ｑ和阿爾法a
%?xlabel（‘權(quán)重因子q‘）ylabel（‘奇異指數(shù)a‘）
%?figureplot（qf）;%%%%%%%%%%%%%ｑ和f（a）
%?xlabel（‘權(quán)重因子q‘）ylabel（‘多重分行譜f（a）‘）

%?ff=find（isnan（f））;
%?f（ff）=[]

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

?????文件???????1914??2012-03-20?10:59??多重分形\fractral.m

?????文件???????4886??2012-10-15?20:05??多重分形\m.asv

?????文件???????4886??2012-10-15?20:06??多重分形\m.m

?????文件???????4212??2012-03-20?15:16??多重分形\mm.m

?????文件??????62913??2012-03-20?09:28??多重分形\N0-red-10X-5500.jpg

?????文件???????2842??2006-12-04?21:21??多重分形\pluslr.m

?????文件???????2233??2012-10-15?20:09??多重分形\sbs.m

?????文件???????3473??2012-10-15?20:10??多重分形\sffs.m

?????文件???????1614??2012-10-15?20:12??多重分形\SFS.m

????..A.SH.?????11264??2012-11-15?09:28??多重分形\Thumbs.db

?????文件???????2790??2012-04-05?20:42??多重分形\多重分形譜算法_1.m

?????文件???????1223??2010-12-29?19:45??多重分形\盒子分形維數(shù).m

?????目錄??????????0??2012-11-15?09:28??多重分形

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

???????????????104250????????????????????13