ssa

1
function [NS,w_avr,w_med,w_iqr] = wav(N,S,W,res,sET)
2
   ET = zeros(N,S);
3
   for j=1:S % öèêë ïî ñåãìåíòàì
4
      for i=1:W
5
         k = (j-1)*res;
6
         ET(i+k,j) = sET(i,j); % ñäâèíóòûé ñåãìåíò ET(:,j)
7
      end
8
   end
9
   TS = zeros(S);
10
   for i=1:N % öèêë ïî îòñ÷åòàì ET
11
      nSi = 0;
12
      for j=1:S % öèêë ïî ñåãìåíòàì ET
13
         if ET(i,j)~=0
14
            nSi      = nSi+1; % òåêóùèé íåíóëåâîé îòñ÷åò ET(i,j)
15
            Smp(nSi) = ET(i,j); 
16
         end
17
      end
18
      NS(i)    = nSi; % êîë-âî ñåãìåíòîâ äëÿ òåêóùåãî i
19
      w_avr(i) = mean(Smp(1:nSi));   % âûáîðî÷íàÿ ñðåäíÿÿ
20
      w_med(i) = median(Smp(1:nSi)); % ìåäèàíà
21
      w_iqr(i) = (prctile(Smp(1:nSi),75)-prctile(Smp(1:nSi),25))/2.0;
22
   end
23
   NS    = NS';
24
   w_avr = w_avr'; % î÷èùåííàÿ ïóëüñîâàÿ âîëíà
25
   w_med = w_med';
26
   w_iqr = w_iqr'; % ïîëîâèíà èíòåðêâàðòèëüíîãî äèàïàçîâà
27
end