[ Takaisin pääsivulle | Back to main page ]

Esimerkki 9 - vastauksia

Differenssiyhtälö ja impulssivaste:
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])
jossa d[n] on yksikköfunktio delta[n].

Laskentaa varten voi kirjoittaa esimerkiksi Matlab-koodia:
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;

Suotimen lohkokaavion ja amplitudivasteen, eli kompleksiarvoisen funktion
H(w) = (1/7) . (1 + e^-jw + e^-j2w + e^-j3w + e^-4jw + e^-j5w + e^-j6w)
itseisarvon saa näppärästi Suodingeneraattorista osoitteesta http://www.cis.hut.fi/Opinnot/T-61.246/Suodin/generateFilter.php?numerator=0.1429+0.1429+0.1429+0.1429+0.1429+0.1429+0.1429+&denumerator=1&G=&type=dfi&luo=Luo+piirrokset+%2F+Create+figures

[ Takaisin pääsivulle | Back to main page ]

http://www.cis.hut.fi/Opinnot/T-61.246/Demo/esim9vast.shtml
t61246@cis.hut.fi
Friday, 10-Sep-2004 12:50:40 EEST