okna, sita in viri

164 19. Sita s FIR
MATLAB 19.3: Izračun koeficientov sita z okenskim postopku pri uporabi Kaiserovega okna.
fs=1000;
% frekvenca vzorčenja
fn=fs/2;
% Nyquistova frekvenca
N=73;
% dolžina sita
beta=5.65;
% parameter valovitosti
fm1=125/fn;
% normalizirana spodnja mejna frekvenca
fm2=275/fn;
% normalizirana zgornja mejna frekvenca
fm=[fm1 fm2];
% vektor mejnih frekvenc
h=fir1(N-1,fm,kaiser(N,beta)); % aproksimacija impulznega odziva
% IZRAČUN IN IZRIS AMPLITUDNE KARAKTERISTIKE
[H,f]=freqz(h,1,fs); % privzame l=512
mag=20*log10(abs(H));
plot(f,mag),grid on
xlabel(’frekvenca [Hz]’)
ylabel(’amplitudni odziv [dB]’)
Koeficienti, ki jih izračuna program 19.3 so zbrani v tabeli 19.7. Amplitudna karakteri-
Tabela 19.7
Koeficienti Kaiserovega pasovnega sita (zgled 19.6.2 na strani 162). Ker je impulzni odziv pozitivno simetričen, je
zapisana le prva polovica koeficientov.
n h[n]
0 −0,0001
1 −0,0004
2 −0,0001
3 −0,0001
4 −0,0007
n h[n]
5 0,0005
6 0,0023
7 0,0008
8 −0,0017
9 −0,0005
n h[n]
10 −0,0005
11 −0,0044
12 −0,0022
13 0,0069
14 0,0066
n h[n]
15 −0,0016
16 0,0000
17 0,0022
18 −0,0117
19 −0,0164
n h[n]
20 0,0069
21 0,0189
22 0,0029
23 0,0044
24 0,0188
n h[n]
25 −0,0125
26 −0,0520
27 −0,0165
28 0,0333
29 0,0104
n h[n]
30 0,0094
31 0,0856
32 −0,0453
33 −0,1665
34 −0,2066
n h[n]
35 0,0891
36 0,2998
stika je na sliki 19.17 na naslednji strani.
♦
Optimalna metoda
Za postopke načrtovanja sit z optimalno metodo nudi program orodje Matlab
v orodnem kovčku Signal Processing Toolbox dve funkciji, ki temeljita na algoritmih,
ki so ji razvili Park in McClellan 12 . Osnova jim je tako imenovani
Remezov izmenjevalni algoritem, ki pomika frekvence, pri katerih ima pogreškovna
funkcija aproksimacije D(ω) − D M (ω) ekstreme. Remez je zanj
12 Glej opombi 8 in 9 na strani 143.
šarko ƒu£ej: Teorija signalov revizija 20040504