okna, sita in viri

162 19. Sita s FIR
Tabela 19.6
Koeficienti sita d[n] (zgled 19.6.1 na predhodni strani) načrtanega z uporabo pravokotnega okna, koeficienti
Hammingovega sita w[n] in koeficienti sita h[n] načrtanega z uporabo Hammnigovega okna. Velja
h[n] = d[n]·w[n]. Ker je impulzni odziv pozitivno simetričen, je zapisana le prva polovica koeficientov. (koeficienti
so iz starega zgleda)
n d[n] w[n] h[n]
0 0,0072 0,0800 0,0006
1 −0,0127 0,0834 −0,0011
2 0,0078 0,0934 0,0007
3 0,0043 0,1099 0,0005
4 −0,0137 0,1327 −0,0018
5 0,0122 0,1614 0,0020
6 0,0000 0,1957 −0,0000
7 −0,0135 0,2350 −0,0032
8 0,0168 0,2787 0,0047
n d[n] w[n] h[n]
9 −0,0058 0,3262 −0,0019
10 −0,0117 0,3769 −0,0044
11 0,0212 0,4299 0,0091
12 −0,0133 0,4846 −0,0065
13 −0,0076 0,5400 −0,0041
14 0,0252 0,5954 0,0150
15 −0,0234 0,6501 −0,0152
16 0,0000 0,7031 0,0000
17 0,0286 0,7538 0,0216
n d[n] w[n] h[n]
18 −0,0378 0,8013 −0,0303
19 0,0140 0,8450 0,0119
20 0,0311 0,8843 0,0276
21 −0,0635 0,9186 −0,0585
22 0,0467 0,9473 0,0443
23 0,0327 0,9701 0,0318
24 −0,1511 0,9866 −0,1493
25 0,2570 0,9966 0,2566
26 0,6986 1,0000 0,6999
MATLAB 19.2: Izračun koeficientov sita po okenskem postopku. Za primerjavo je dodan še izračun koeficientov
Hammingovega okna in koeficientov sita pri pravokoetnem oknu. Sito aproksimira idealno nizkopasovno sito.
fm=0.7;
N=53;
w_Ham=hamming(N);
h=fir1(N-1,fm,w_Ham);
d=fir1(N-1,fm,boxcar(N));
% normalizirana mejna frekvenca
% dolžina sita
% koeficienti Hammingovega okna
% aproksimacija s Hammingovim oknom
% aproksimacija s pravokotnim oknom
Izračunani koeficienti d[n], w[n] in h[n] so zbrani v tabeli 19.6. Po vzoru programa za
izračun frekvenčne karakteristike (program MATLAB 19.1 na strani 159) še izračunamo
in izrišemo amplitudni karakteristiki za siti z impulznim odzivom d[n] (slika 19.16a) in
h[n] (slika 19.16b).
V program 19.2 bi na začetku lahko vključili še izračun števila koeficientov sita in
izračun normirane mejne frekvence:
N=ceil(3,3/2*0,0625)
fm=2*(1,5+0,25)/5
in si s tem prikrajšali posebno računanje reda sita in normirane frekvence.
♦
ZGLED 19.6.2 (Kaiserovo okno)
Izračunajmo koeficiente sita s FIR z okensko metodo, kjer uporabimo Kaiserovo okno.
Sito naj izpolni naslednje zahteve:
šarko ƒu£ej: Teorija signalov revizija 20040504