Skrypt APSC - MARS

ROZDZIAŁ 7. FILTRY CYFROWE 1317.5.2 Metoda próbkowania w dziedzinie częstotliwościDo projektowania odpowiedzi impulsowej tą metodą przewidziana jest funkcjafir2. Podaje się jej parametry po kolei: rząd filtru, wektor punktów częstotliwościi wektor wartości amplitud w tych punktach. Przykładowy program:f s = 1000; % c z ę s t o t l i w o ś ć próbkowania [ Hz ]f 0 = 400; % c z ę s t o t l i w o ś ć granicznaN = 30; % rządf = [ 0 f 0 / f s f 0 / f s 1 ] ; % wektor punktów c z ę s t o t l i w o ś c im = [ 1 1 0 0 ] ; % wektor amplitudh = f i r 2 (N, f ,m) ;[H,w] = f r e q z (h , 1 , 5 1 2 ) ;p l o t ( f s /2∗w/ pi , 20∗ log10 ( abs (H) ) )g r i d onWektor f może być dłuższy, czyli zawierać więcej punktów. Każdemu z tychpunktów częstotliwości musi odpowiadać wartość amplitudy (wektor m). W tensposób określa się kształt charakterystyki amplitudowo-częstotliwościowej.W przykładzie pokazano również w jaki sposób można przypisać zwracanewartości funkcji freqz do wektorów - odpowiedzi częstotliwościowej (H, zespolone)i częstotliwości (w), aby wykorzystać te dane do dalszych obliczeń lubprezentacji danych. Zadanie.Zaprojektuj filtr o paśmie przepustowym od 300Hz do 2, 5kHz sygnałuspróbkowanego z częstotliwością 8kHz. Porównaj charakterystyki dla kilku wybranychwartości N.7.5.3 Metoda RemezaProjektowanie filtrów metodą Remeza (algorytmem Parks-McClellan’a) zapewniafunkcja firpm. Kolejność i znaczenie parametrów jest bardzo podobne jakw funkcji fir2, przy czym wartości w wektorze częstotliwości i amplitudy podajesię parami tworząc przedziały zmienności charakterystyki amplitudowoczęstotliwościowej.Zależność pomiędzy wektorem częstotliwości i amplitudypokazano na rysunku 7.15 (rysunek zaczerpnięty z pliku pomocy programu Matlab).Zadanie.Zbadaj charakterystykę amplitudowo-częstotliwościową filtru o charakterystycezadanej jak na rysunku 7.15.