1. uvod u digitalnu sliku - Laboratorija za digitalnu obradu signala
1. uvod u digitalnu sliku - Laboratorija za digitalnu obradu signala
1. uvod u digitalnu sliku - Laboratorija za digitalnu obradu signala
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
praktični dio<br />
Opis funkcije<br />
y = resample(x,p,q) mijenja frekvenciju odmjeravanja tako da je nova frekvencija<br />
umnožak faktora p/q i originalne frekvencije. Pri tome se koriste polifazni filtri, dok su p i q<br />
moraju biti pozitivne cjelobrojne vrijednosti. Dužina izlaznog ni<strong>za</strong> y je jednaka dužini ulaznog<br />
ni<strong>za</strong> pomnoženo sa faktorom p/q a <strong>za</strong>tim <strong>za</strong>okruženo na prvu veću cijelu vrijednost. Ukoliko<br />
je x matrica, ova fukcija radi sa kolonama te matrice. Ova funkcija takođe primjenjuje antialiasing<br />
NF FIR filtar tokom procesa resamplinga. FIR filtar koji se koristi je firls (least<br />
square linear phaseFIR filter odnosno FIR filtar sa linearnom fazom i metodom najmanjih<br />
kvadrata) i to sa Kaiser prozorom.<br />
y = resample(x,p,q,n) koristi n odmjeraka sa obe strane trenutnog odmjerka pri<br />
računanju novih odmjeraka. Dužina FIR filtra je prooporcionalna n , ukoliko koristimo veće<br />
vrijednosti n dobijamo veću tačnost, međutim povećava se vrijeme potrebno <strong>za</strong> računanje.<br />
Početna vrijednost je 10. Ako stavimo da je n=0 dobijamo interpolaciju na osnovu najbližeg<br />
susjeda.<br />
y = resample(x,p,q,n,beta) koristi beta kao parametar pri di<strong>za</strong>jniranju Kaiserovog<br />
prozora koji se koristi pri u toku resamplinga. Početna vrijednost <strong>za</strong> beta je 5.<br />
y = resample(x,p,q,b) filtrira niz x koristeći koeficijente filtra date u vektoru b.<br />
[y,b] = resample(x,p,q) kao rezultat pored izlaznog ni<strong>za</strong> y vraća i vektor b koji<br />
predstavlja koeficijente filtra koji se koristi tokom procesa resamplinga.<br />
� upfirdn<br />
Vrši tri operacije, prvo konverziju naviše, <strong>za</strong>tim primjenjuje FIR filtar i na kraju<br />
konverziju naniže. Kaskada ove tri operacije se primjenjuje kroz tehnike polifaznog filtriranja,<br />
a to je glavna tehnika u multirate filtriranju. Pri tome kvalitet resamplinga <strong>za</strong>visi od kvaliteta<br />
primjenjenog FIR filtra.<br />
Sintaksa<br />
yout = upfirdn(xin,h)<br />
yout = upfirdn(xin,h,p)<br />
yout = upfirdn(xin,h,p,q)<br />
[yout,zf] = upfirdn(xin,h,...,zi)<br />
Opis funkcije<br />
Ova funkcija vrši tri operacije jednu <strong>za</strong> drugom:<br />
- Konverziju naviše ulaznih podataka sa faktorom p (ubacivanje nultih odmjeraka)<br />
- FIR filtriranje ovakvog <strong>signala</strong> sa filtrom čiji je impulsni odziv dat vektorom h<br />
- Konverziju naniže sa faktorom q (odbacivanje odmjeraka)<br />
Ova funkcija je realizovana kao MEX fajl da bi se dobila maksimalna brzina<br />
izvršavanja tako da je potrebno obrađivati samo one odmjerke koji će se naći na izlazu. FIR<br />
filtar koji se koristi je obično NF filtar, i on se mora prethodno projektovati koristeći funkcije<br />
remez ili fir<strong>1.</strong><br />
yout = upfirdn(xin,h) filtrira ulazni signal xin korištenjem FIR filtra sa impulsnim<br />
odzivom datim sa vektorom h . Ukoliko je ulazni signal vektor, on se posmatra kao jedan<br />
signa. Međutim, ako je ulazni signal matrica, svaka kolona se posmatra kao drugi signal i<br />
<strong>za</strong>sebno se filtrira. Ako je h vektor, tada on predstavlja jedan FIR filtar. Ukoliko je on<br />
matrica, svaka kolona predstavlja poseban filtar.<br />
yout = upfirdn(xin,h,p) koristi faktor konverzije naviše p, početna vrijednost mu je <strong>1.</strong><br />
yout = upfirdn(xin,h,p,q) faktor konverzije naniže q početna vrijednost mu je <strong>1.</strong><br />
[yout,zf] = upfirdn(xin,h,...,zi) određuje početne uslove <strong>za</strong> vektor zi. Pri tome dužina<br />
vektora početnih uslova mora biti niža <strong>za</strong> jedan od dužine vektora h.<br />
56