multifractal 多重分形谱的算法的matlab实现

多重分形谱的算法的matlab实现，输入数据为空，请自行选择

function multifractal（）
X=xlsread（‘XXX.xls‘）;%在这里输入数据文件
% 提取数据
[mn] = size（sh）;
m=3800;
A = X（1:m4）;
format long g
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=

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

     文件       2903  2013-04-27 14:34  multifractal1.m

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

                 2903                    1