多重分形及其运算方法！

多重分形程序的运行，应用于matlab中，以及其运行结果的运算！

% 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%%%%%%这个才是真正的谱差啊啊啊啊啊啊啊
    
%      amin=find（isnan（hurstH））;
%      hurstH（hursti）=[];

%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf
%%%%%%%%%%%%%%%f形谱宽ｄｅｌｔａa 谱差ｄｅｌｔａf
figureplot（rf）;
xlabel（‘奇异指数a‘）ylabel（‘多重分行谱f（a）‘）
% q=-60:2:60;
% figureplot（qr）;%%%%%%%%%%%%%ｑ和阿尔法a
% xlabel（‘权重因子q‘）ylabel（‘奇异指数a‘）
% figureplot（qf）;%%%%%%%%%%%%%ｑ和f（a）
% xlabel（‘权重因子q‘）ylabel（‘多重分行谱f（a）‘）

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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       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  多重分形\盒子分形维数.m

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

----------- ---------  ---------- -----  ----

               104250                    13