okna, sita in viri

176 19. Sita s FIR
Amplitudna karakteristika, ki jo izračunamo iz izračunanega impulznega odziva,
se v točkah vzorčenja ujema prototipnim sitom. Tega lahko narišemo s
funkcijo plot(f,m).
Funkcija fir2 računa koeficiente sita z inverzno Fourierovo transformacijo.
Če posebej ne navedemo tip okna, fir2 pri njenem računanju uporabi
Hammingovo okno. Če želimo drugo okno, ga navedemo. Na primer
ali direktno, na primer:
window = blackman(r+1)
b = fir2(r,f,m,window)
b = fir2(r,f,m,kaiser(r+1,beta))
Vektor, ki določa okno, mora imeti r+1 koeficientov.
Potek prototipne amplitudne karakteristike, ki jo določata vektorja f in m,
funkcija fir2 upošteva kot oporne točke pri interpolaciji prototipne karakteristike,
ki je določena skupaj v 512 točkah. Če želimo spremeniti to število,
uporabimo parameter npt. Uporabljamo ga lahko brez ali z določitvijo okna:
h = fir2(r,f,m,npt)
h = fir2(r,f,m,npt,window)
Če sta v vektorju f dve zaporedni frekvenci enaki, funkcija fir2 določi
prehodni pas s središčem pri tej frekvenci, v katerem amplitudna karakteristika
naredi gladek prehod med želenima amplitudama (levo in desno limito
prototipne amplitudne karakteristike). Širino prehodnega pasu nadzira parameter
lap. Uporabljamo ga lahko brez ali z določitvijo okna:
b = fir2(n,f,m,npt,lap)
b = fir2(n,f,m,npt,lap,window)
☞
Pri uporabi parametra lap moramo pred njim vedno navesti tudi npt, tudi
če uporabljamo prednastavljeno vrednost! Če lap ni naveden, fir2 uporabi
prednastavljeno vrednost, ki znaša 25.
ZGLED 19.6.7 (Postopek načrtovanja z vzorčenjem frekvenčne karakteristike)
Po postopku vzorčenja frekvenčne karakteristike določimo koeficiente aproksimacije
idealnega nizkopasovnega sita z mejno frekvenco 0,6. Sita bodi 30. reda. Zanj narišimo
prototipno in dejansko amplitudno karakteristiko.
REŠITEV: Potek amplitudne karakteristike idealnega nizkopasovnega sita ima preskok
pri mejni frekvenci 0,6 (normirano). Opišemo jo z:
f = [0 0.6 0.6 1]
m = [1 1 0 0]
šarko ƒu£ej: Teorija signalov revizija 20040504