Esimerkki 9 - vastauksia

y[n] = (1/7) . (x[n] + x[n-1] + x[n-2] + x[n-3] + x[n-4] + x[n-5] + x[n-6])
h[n] = (1/7) . (d[n] + d[n-1] + d[n-2] + d[n-3] + d[n-4] + d[n-5] + d[n-6])

n  = [1:500];
x  = cos(2*pi*10/1000*n);
x2 = x + 0.1*randn(1,length(x));  % generoitiin kohinainen kosinisekvenssi (syöte)

N  = 7;                           % MA-N, N-point moving average
y  = zeros(1,length(n));          % alustetaan vaste nollaksi

for k = N : length(x)
  y(k) = sum(x2(k-N+1:k)) / N;    % summataan N kpl ja jaetaan N:llä
end;                              
% HUOM! Matlabissa tämän voi tehdä "kauniimminkin".

plot(n, x2, 'b', n, y, 'r', n, x, 'g');  
% piirretään alkuper. sinisellä,
%  keskiarvoistettu punaisella
%  ja "kohinaton signaali vihreällä
legend({'alkuper. kohinainen','suodatettu','kosini'});
xlabel('aikaindeksit n');
grid on;

H(w) = (1/7) . (1 + e^-jw + e^-j2w + e^-j3w + e^-4jw + e^-j5w + e^-j6w)