%根據計盒維數原理編寫了求一維曲線分形維數的matlab程序?

function?D=FractalDim（ycellmax）

%求輸入一維信號的計盒分形維數
%y是一維信號
%cellmax:方格子的最大邊長可以取2的偶數次冪次（1248...）取大于數據長度的偶數
%D是y的計盒維數（一般情況下D>=1）D=lim（log（N（e））/log（k/e））
f_v=‘y.txt‘;
y=load（f_v）;
cellmax=16384;
if?cellmaxerror（‘cellmax?must?be?larger?than?input?signal!‘）
end
L=length（y）;%輸入樣點的個數
y_min=min（y）;

%移位操作，將y_min移到坐標0點
y_shift=y-y_min;
%重采樣，使總點數等于cellmax+1
x_ord=[0:L-1]./（L-1）;
xx_ord=[0:cellmax]./（cellmax）;
y_interp=interp1（x_ordy_shiftxx_ord）;
%按比例縮放y，使最大值為2^c
ys_max=max（y_interp）;
factory=cellmax/ys_max;
yy=abs（y_interp*factory）;

t=log2（cellmax）+1;%疊