Расчет размерностей Реньи и корреляционного интеграла

MaxVlozhenie = 20;

dc=zeros(l, MaxVlozhenie);

Corrlnt=zeros (10, MaxVlozhenie);

form = 1: MaxVlozhenie

data = zeros (length (signal (:,2)), m);

data (:,1) = detrend (signal (:,2));

for і = 2: m

data (:, i) = circshift (data (:,1), — i+1);

end

data = data (m: length (data),:);

SD = std (data (:, m));

fore=l:10

Epsilon (e) = (l/e)*SD;

for і = 1: length (data)

forj = 1: length (data)

ifi~=j

N=(l/(length (data)* (length (data)-l)));

Corrlnt(e, m) = Corrlnt (e, m) + N*0.5*...

(sign (Epsilon (e)-norm (data (i,:)-data (j,:),!))+!);

end

end

end

end

C=polyfit(log(Epsilon ((Corrlnt(:, m))>0)’),...

log (Corrlnt ((Corrlnt (:, m)) >0, m)), 1);

dc(m)=C(l);

end

figure

plot (de,

hold on

M= (1: MaxVlozhenie);

Cd=polyfit (M, de, 2);

plot (M, Cd (1) *M. V+Cd (2) *M+Cd (3), ’);

% Построение линии тренда и определение корреляционной размерности

106

Вычисление аппроксимированной энтропии

Data=zeros (length (NN (:,2)),3);

data (:, 1) = detrend (NN (:,2));

Сопроводим смещение фазы на т точек data (:,2) = circshift (data (:, 1), 1); data (:,3) = circshift (data (:, 1),2);

%3адаем разбиение ячеек г от СКО

SD = std (data (:, 1)); г = 0.2*SD;

Vs2=zeros (length (data),l); Vs3=zeros (length (data),l); for і = 2: length (data) forj = 1: length (data)

Vs2(i) = Vs2(i) + 0.5* (sign (r-abs (data (i,l)...

  • - data (j,!)))+!)...
  • *0.5* (sign (r-abs (data (i,2)-data (j,2)))+l);

Vs3 (i) = Vs3 (i) + 0.5* (sign (r-abs (data (i,l)...

  • - data (j,!)))+!)...
  • *0.5* (sign (r-abs (data (i,2)-data (j,2)))+l)...
  • *0.5* (sign (r-abs (data (i,3)-data (j,3)))+l); end; end; end;

V2=sum (Vs2); V3=sum (Vs3); ApEn=log (V2./V3);

Приложение 13

Расчет показателя Херста методом накопленной дисперсии

%максималъная длина приращения

ртах=10;

L=length (signal);

for р=1: ртах

for i=l: L-pmax

%набор приращений и их СКО

dX (i)=signal (i+p,2)~ signal (і,2);

S(p)=std(dX);

end

E (p)=log (p);

N(p)=log(S(p));

end

plot (E, N, ’o’);

%коэффициенты MНК-прямой

H=polyfit (E, N,1); %MНК-прямая

h=H (2)+H (1)*E;

hold on

plot (E, h, ’г’); Неполученный показатель Херста

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