Оценка показателя Ляпунова

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))0 a(l)=i;

%проверка выполнения условия е-близости 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);

 
Посмотреть оригинал
< Пред   СОДЕРЖАНИЕ   ОРИГИНАЛ   След >