% Toto je moj vlastny plavajuci priemer
% 
% Verzia 1.0 na prvy pokus 5. 6. 2008 R.B.
%

function y = movingaverage (x, M)

%  x je vstupny vektor, ktory filtrujeme
%  M je rad filtra, kolko clenov priemerujeme

 y=[zeros(1,M-1) x];     % predlzime vektor o M-1 nul na zaciatku
                         % aby sa aj prva hodnota zaratala do priemeru
 for i=1 : length(x)
  y(i) = sum( y( i : (i+M-1) ) );       % spocitavame M-tice
 endfor
 
 y=y/M;                                 % napokon to vydelime

 y=y(1:length(x));                      % a vektor zasa skratime o ten 
                                        % pridavok, lenze na konci

endfunction
