Оценка показателя Ляпунова
data=NN (:, 2)/1ООО; %переход к значениям в секундах epsilon=20; %параметр близости epsilon=epsilon/1000; %переход к значениям в секундах tau_max=15; %параметр вложения
N=length (data)-tau max-1;
distanc=cell (N, tau_max+l);
fort=l:N
for tau=O: tau_max
1=1;
clear dist
for i=l: length (data)-tau max
if (abs (data (i)-data (t))
%проверка выполнения условия е-близости dist (l)=abs (data (i+tau)-data (t+tau));
distanc{t, tau+l}=dist;
1=1+1;
end;
end;
close_tr{t, l}=a;
end;
end;
vectau=0: tau max;
vectau=mean (data)*vectau; Ренормировка вектора т fortau=0:tau max
N1=1;
fort=l:N
v_summy=distanc{t, tau+1};
if length (v_summy)>15
%учет только тех данных, у которых не менее 15 е-соседей
staticticka (tau+1, t)=log(sum ((1/length (v_summy))*v_summy));
N1=N1+1;
end;
end;
end;
stat= (1/Nl)*sum (staticticka,2);
hold on
grid on
plot (vectau, stat)
Приложение 10
Вычисление размерностей Реньи методом boxcounting
NN=NN/1000; %переход к секундам
NN (:,2)=NN (:,2)-abs (min (NN (:,2)));
NN (:,2)=NN (:,2)/max (abs (NN(:,2)));
%перенормировка исходного ряда значений
NN (:, 1) =NN (:, 1)+0.001;
%сдвиг исходного ряда аргументов
NN (:, 1) =NN (:, 1)/1000;
%перенормировка исходного ряда аргументов
тахХ=тах (NN); minX=min (NN);
boxsize =linspace (0.01,0.20,20);
%вектор размеров ячеек
boxsize =boxsize’;
counter=zeros (length (boxsize), 1);
Shenon=zeros (length (boxsize), 1);
sD2=zeros (length (boxsize),!);
p=cell (length (boxsize),2);
for i=l: length (boxsize)
counter (i)=0;
box num (i,:)=ceil ((abs (maxX)+abs (minX))/boxsize (i)); доопределение полного количества ячеек для анализа box_left=zeros (box_num (i,l), box_num (і,2));
box_right=zeros (box_num (i,l), box_num (і,2));
box_top=zeros (box_num (i,l), box_num (i,2));
box_bottom=zeros (box_num (i,l), boxjnum (i,2));
n=zeros (box_num (i,l), boxjnum (i,2));
for 1=1: boxjnum (i,2)
forh=l: box_num (i,l) box_left (h, l)= (h-1)* boxsize (i);
box_right (h, I) =h*boxsize (i);
box_top (h, l)= (I)* boxsize (i);
boxjbottom (h, l)= (1-І)*boxsize (i);
if not (isempty (NN (NN (:,2) >box bottom (h, I)...
&NN(:,2)<=box_top (h, 1)&NN (:,l)>box_left (h, I)...
&NN (:, l)<=box_right (h, I),2)))
counter (i)=counter (i)+l;
%количество ячеек, перекрытия BP
end;
ifbox num (i,2)==0
counter=box_num (i,l);
end;
n (h, l)=length ((NN(NN(:,2)>box_bottom (h, I)...
&NN(:,2)<=box_top (h, 1)&NN (:,l)>box_left(h, I)...
&NN(:,l)<=box_right(h, I),2)));
%число точек в каждой ячейке
end;
end;
p{i, l}=n/counter (і);
Соотносительная заселенность ячейки
Norm (i)=sum (sum ((p{i, l})));
%условие нормировки
p{i,l}=p{i,l}/Norm (і);
p{i,2}=log(p{i,l});
Pfi,2} (p{i,2}==-inf)=0;
p{i,2} (isnan (p{i,2}))=0;
Shenon (i)=sum (sum ((p{i,l})* (p{i,2})));
%информация Шенона
sD2 (i)=log (sum (sum ((p{i, 1}). * (p{i, 1}))));
end;
C=polyfit (log (1./boxsize), log (counter), 1); dO=C (1);
C^polyfit (log (boxsize), (Shenon), 1); dl=C (1);
C=polyfit (log (boxsize), (sD2), 1); d2=C (1);