duochongfenxing.m

多重分形的MATLAB算法程序 自己做的 大家多指教

% 此函数的功能为求多重分形谱函数，并作出谱函数的图像
% y是存储配分函数的自然对数值
% x是存储q数组的对数的数组
% h1表示了用几个q值作为一组来求a
% h2表示了从a的最小值mina到最大值maxa的变化步长
q=(-4.51:0.0001:4.51);%列出了q值的取值向量数组，可以根据需要取值 
     T=[];           %存储配分函数Tq的值
     aa=[];          %用于存储所有的a值
   X=[];Y=[];
   n=1;
h1=3;
   maxi=fix(length(q)/h1);
   for i=1:maxi                        %先把所有的a都求出来，再求出a的极值
       X=q(h1*(i-1)+1:h1*i);
       Y=T(h1*(i-1)+1:h1*i);
       p=polyfit(X,Y,1);
       aa(n)=p(1);
       n=n+1;
   end
   amin=min(aa);
   amax=max(aa);
   F=[];
   m=1;
h2=0.001;
   for A=amin:h2:amax                  %按对多重分形谱f(a)的离散定义求出f(a)
       b=[];
       l=1;
       for i=1:length(q)
           b(l)=A*q(i)-T(i);
           l=l+1;
       end
       F(m)=min(b);
       m=m+1;
   end

   a00=(amin:h2:amax);
   plot(a00,F,'^b');