okna, sita in viri

19.6 Uporaba programa MATLAB pri načrtovanju sit s FIR 173
frekvenca vzorčenja: 2 kHz −−−−−−→ normiranje
2
mejna frekvenca prepustnega pasu: 0,5 kHz −−−−−−→ normiranje
0,5
mejna frekvenca zapornega pasu: 0,6 kHz −−−−−−→ normiranje
0,6
Red sita in njegove koeficiente impulznega odziva izračunamo s programom 19.7.
MATLAB 19.7: Izračun koeficientov nizkopasovnega sita pri uporabi normiranih frekvenc v funkciji remezord
d_p = 3;
% valovitost v dB v prepustnem pasu
d_z = 40;
% valovitost v dB v zapornem pasu
dev = [(10^(d_p/20)-1)/(10^(d_p/20)+1) 10^(-d_z/20)];
f = [0.5 0.6];
% vektor mejnih frekvenc
a = [1 0];
% želene amplitude
[r,fo,ao,w] = remezord(f,a,dev); % računanje reda sita
h = remez(r+1,fo,ao,w);
% računanje koeficientov sita
% upoštevamo red r+1
% IZRAČUN IN IZRIS AMPLITUDNE KARAKTERISTIKE
[H,f]=freqz(h,1,1024);
mag=20*log10(abs(H));
plot(f/pi,mag),grid on;
% f/pi: frekvenca v Hz
xlabel(’frekvenca [Hz]’);
ylabel(’amplitudni odziv [dB]’)
20
0
amplitudni odziv [dB]
−20
−40
−60
−80
Slika 19.22
Amplitudna karakteristika sita iz zgleda 19.6.5.
−100
0 0.2 0.4 0.6 0.8 1
frekvenca [Hz]
Da smo izpolnili zastavljene tolerance, smo v funkciji remez namesto reda r uporabiti
red r+1.
♦
ZGLED 19.6.6 (pasovnoprepustno sito z enakomerno valovitostjo)
Za pasovno prepustno sito z linearno fazno karakteristiko in s parametri
datoteka: signal_C