Skrypt APSC - MARS

ROZDZIAŁ 7. FILTRY CYFROWE 1291. określenie parametrów filtru: kształtu charakterystyki amplitudowej (iewentualnie fazowej) w dziedzinie częstotliwości, maksymalnych zafalowańw paśmie przepustowym, minimalnego tłumienia w paśmie zaporowym,szerokość pasma przejściowego; wszystkie częstotliwości najlepiejunormować względem częstotliwości próbkowania,2. wyznaczenie odpowiedzi impulsowej,3. wybór okna i wymnożenie jego funkcji z odpowiedzią impulsową (dotyczymetody okien czasowych i próbkowania w dziedzinie częstotliwości),4. sprawdzenie charakterystyki amplitudowo-częstotliwościowej i fazowo-częstotliwościowej.7.5.1 Metoda okien czasowychOdpowiedź impulsową filtru projektowanego tą metodą zwraca funkcja fir1.Jako parametr należy podać rząd filtru n i częstotliwość graniczną Wn unormowanądo połowy częstotliwości próbkowania (częstotliwości Nyquista). JeśliWn jest wektorem 2-elementowym można utworzyć filtr pasmowoprzepustowy.Trzeci parametr może określić typ filtru na dolnoprzepustowy lub pasmowozaporowy.Domyślnie program stosuje okno Hamminga, ale można jako kolejnyparametr funkcji podać inne okno. Dokładny opis znajduje się w dokumentacjiMatlaba (polecenie doc fir1).Do szybkiego sprawdzenia amplitudy i fazy w dziedzinie częstotliwości zaprojektowanegofiltru używamy funkcji freqz. Pierwszym parametrem jest wektorwspółczynników filtru (licznik transmitancji), drugim współczynniki mianownika(dla filtrów FIR wstawiamy 1), trzecim może być liczba punktów tworzącacharakterystykę wyjściową. Wywołując funkcję bez przypisania wartości zwracanychdo wektora tworzony jest wykres.Przedstawmy na prostym przykładzie projekt filtru dolnoprzepustowego o32 współczynnikach i częstotliwości granicznej równej 0,2 częstotliwości próbkowania.h = f i r 1 ( 3 2 , 0 . 2 ) ;f r e q z (h , 1 , 5 1 2 )W projekcie wygodniej, zamiast częstotliwością unormowaną posługiwać sięwartościami bezwzględnymi wyrażonymi w Hercach. Sekwencja poleceń w środowiskuMatlaba może być następująca: