function Tsr = Tsrednje_NOVA (kut,kutLMA)
% odreðujemo prvo broj kontrakcija u signalu,a zatim srednje vrijeme jedne
% jedne kontrakcije Tsr [ms]

% TsrLMA je broj uzoraka u postupku usrednjavanja, mora biti neparan broj (odd)
% zasad najbolje poklapanje za LMA=501

Adem=abs(hilbert(kut));
meanAdem= Adem-mean(Adem);

% provodimo DVOSTRUKO FILTRIRANJE da ukupno kašnjenje bude NULA
x1= smooth(meanAdem,kutLMA);
x2= smooth(fliplr(x1'),kutLMA);
x= fliplr(x2');

% srednja vrijednost izglaðenog signala
sr= 0;

% tražimo prvo indekse elemenata koji prolaze kroz srednju vrijednost SR,a
% zatim ukupni broj prolaza dijelimo sa 2 da dobijemo broj kontrakcija
i= [1:length(kut)-1];
broj_prolaza= find((x(i+1)>sr & x(i)<sr) | (x(i+1)<sr & x(i)>sr));
broj_kontrakcija= length(broj_prolaza)/2;

% srednje trajanje kontrakcije
Tsr= length(kut(1:broj_prolaza(end-1)))/(broj_kontrakcija-1);




