okna, sita in viri

170 19. Sita s FIR
V primerih, ko so mejne frekvence blizu ničle ali Nyquistove frekvence,
se lahko zgodi, da funkcija remezord izračuna premajhen reda sita. To
uvidimo s primerjavo doseženih in želenih specifikacij. V takem primeru
moramo izračun koeficientov sita ponoviti pri redu sita r+1 ali r+2.
ZGLED 19.6.4 (nizkopasovno optimalno sito: uporaba dejanskih frekvenc)
Določimo red nizkopasovnega sita, ki še izpolni naslednje zahteve:
mejna frekvenca prepustnega pasu: 0,5 kHz valovitost v prepustnem pasu: 3 dB
mejna frekvenca zapornega pasu: 0,6 kHz slabljenje v zapornem pasu 40 dB
frekvenca vzorčenja:
in izračunajmo koeficiente sita!
2 kHz
REŠITEV: Pri načrtovanju uporabimo funkcijo remezord s sintakso, ki dovoljuje
uporabo dejanskih frekvenc. Ker so koeficienti deviacije podani v decibelih, jih moramo
z (19.85) in (19.87) pretvoriti v absolutno vrednost. To lahko naredimo pomočjo
programa Matlab. Pri tem lahko koeficiente vektorja deviacije izračunamo posebej
(program MATLAB 19.5)
MATLAB 19.5: Izračun vektorja deviacije amplitudne karakteristike
A_p = 3;
% valovitost v dB v prepustnem pasu
A_z = 40;
% valovitost v dB v zapornem pasu
d_p = (10^(A_p/20)-1)/(10^(A_p/20)+1);
d_z = 10^(-A_z/20);
dev = [d_p d_z];
% vektor absolutne valovitosti
ali pa kar v okviru definicije vektorja dev (program MATLAB 19.6 na naslednji strani).
Amplitudno karakteristiko načrtovanega sita izračunamo in narišemo na običajni način.
Če želimo videti še fazno karakteristiko – na njej se sprememba reda sita bolj opazi –
program dopolnimo še s funkcijami za njen izris.
šarko ƒu£ej: Teorija signalov revizija 20040504