图像分形的matlab算法,多重分形谱算法

function multifractal(A)

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(i,1)=mod(L,r);

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(B,r,length(B)/r);

T=mean(U);

T=T'.*r^3;

TT(1: length(T),i)=[T];

modifying(modify,1)=length(T);

r=r+1;

if r>= tmin+n