Filtr cyfrowy (FIR/IIR) — wykres Bodego i charakterystyka częstotliwościowa
Interaktywny projektant filtrów cyfrowych FIR/IIR. Zobacz wykres Bodego (amplituda + faza), rozkład zer i biegunów oraz filtrację sygnału na żywo. Projektuj filtry Butterwortha, Czebyszewa i okienkowane sinc.
🎛️ Filtr cyfrowy (FIR/IIR) — projektowanie DSP
Filtry cyfrowe leżą u podstaw całego nowoczesnego przetwarzania sygnałów — od korekcji dźwięku i rozpoznawania mowy po obrazowanie medyczne i radio definiowane programowo. Ten symulator pozwala interaktywnie projektować filtry FIR i IIR, obserwując w czasie rzeczywistym wykres Bodego (amplituda i faza) oraz rozkład zer i biegunów podczas przesuwania suwaków.
🔬 Co przedstawia
Filtry cyfrowe selektywnie przepuszczają lub blokują zakresy częstotliwości. Filtry FIR (o skończonej odpowiedzi impulsowej) wykorzystują wyłącznie odczepy w przód — są zawsze stabilne i mają dokładnie liniową fazę (co oznacza, że wszystkie częstotliwości są opóźniane jednakowo, zachowując kształt fali). Filtry IIR (o nieskończonej odpowiedzi impulsowej) wykorzystują sprzężenie zwrotne, wymagając mniejszej liczby odczepów do uzyskania tej samej stromości zbocza, ale mogą stać się niestabilne (bieguny muszą pozostać wewnątrz okręgu jednostkowego) i wprowadzają nieliniową fazę. Wykres Bodego pokazuje charakterystykę amplitudową (jak bardzo tłumiona jest każda częstotliwość) i fazową (jak bardzo jest opóźniana). Rozkład zer i biegunów na płaszczyźnie Z pokazuje stabilność: bieguny poza okręgiem jednostkowym → niestabilność.
🎮 Jak korzystać
Wybierz typ filtru (dolnoprzepustowy, górnoprzepustowy, pasmowoprzepustowy) i metodę projektowania (FIR z oknem sinc lub IIR Butterwortha/Czebyszewa). Przesuń suwak częstotliwości granicznej i obserwuj wykres Bodego aktualizujący się w czasie rzeczywistym. Zwiększ rząd filtru, aby uzyskać bardziej strome zbocze (ostrzejsze przejście z pasma przepustowego do zaporowego) kosztem większej liczby obliczeń. Przełącz na Czebyszewa i dostosuj tętnienia — uzyskasz bardziej strome początkowe zbocze w zamian za tętnienia w paśmie przepustowym. Wypróbuj sygnał testowy szumu białego i obserwuj, jak filtr usuwa wysokie lub niskie częstotliwości z podglądu przebiegu.
💡 Czy wiesz, że?
Kodek audio Twojego smartfona stosuje dziesiątki filtrów cyfrowych co milisekundę — korektor, redukcja szumów, usuwanie echa i filtrowanie antyaliasingowe. Kompresja MP3 wykorzystuje 32-pasmowy bank filtrów polifazowych do podziału dźwięku na podpasma częstotliwości. Metoda FIR z oknem sinc została spopularyzowana przez Freda Harrisa w jego przełomowym artykule IEEE z 1978 roku „On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform”. Nowoczesne radia definiowane programowo (SDR) oparte na FPGA uruchamiają filtry FIR z tysiącami odczepów przy częstotliwościach próbkowania rzędu GHz, wykorzystując potokowe jednostki mnożąco-akumulujące (MAC).
O symulatorze filtra cyfrowego (FIR/IIR)
To narzędzie projektuje filtry cyfrowe i pokazuje ich zachowanie w czasie rzeczywistym. Oblicza charakterystykę częstotliwościową, obliczając funkcję przejścia H(z) = B(z)/A(z) wzdłuż okręgu jednostkowego z = e^(jw), rysując wykres Bodego amplitudy (w dB) i fazy (w stopniach), a także rozkład zer i biegunów na płaszczyźnie Z, znajdowany za pomocą solvera pierwiastków Duranda-Kernera. Częstotliwość próbkowania jest ustalona na 44,1 kHz.
Wybierasz typ charakterystyki (dolnoprzepustowa, górnoprzepustowa, pasmowoprzepustowa, pasmowozaporowa) i metodę projektowania: FIR z oknem sinc (idealna odpowiedź impulsowa okienkowana oknem Hamminga, do rzędu 64), IIR Butterwortha lub IIR Czebyszewa (kaskady biquad, ograniczone do rzędu 8). Suwaki ustawiają częstotliwość graniczną, drugą częstotliwość graniczną dla filtrów pasmowych oraz tętnienia w paśmie przepustowym Czebyszewa. Sygnał testowy jest filtrowany na żywo, dzięki czemu widać efekt na przebiegu. Takie filtry są podstawą korektorów audio, usuwania szumów i odbiorników radiowych.
Najczęściej zadawane pytania
Czym jest filtr cyfrowy?
Filtr cyfrowy to algorytm, który selektywnie przepuszcza lub tłumi fragmenty próbkowanego sygnału w zależności od częstotliwości. Działa poprzez łączenie obecnych i przeszłych próbek wejściowych, a w przypadku filtrów ze sprzężeniem zwrotnym także przeszłych próbek wyjściowych, przy użyciu ustalonego zestawu współczynników liczbowych. Wynikiem jest nowy sygnał z wybranymi częstotliwościami wzmocnionymi lub usuniętymi.
Jaka jest różnica między filtrami FIR i IIR?
Filtry FIR (o skończonej odpowiedzi impulsowej) wykorzystują wyłącznie odczepy w przód, więc ich mianownik A(z) wynosi po prostu 1; są zawsze stabilne i mogą mieć dokładnie liniową fazę. Filtry IIR (o nieskończonej odpowiedzi impulsowej) dodają sprzężenie zwrotne poprzez A(z), osiągając zadaną stromość zbocza przy znacznie mniejszej liczbie współczynników, ale mogą stać się niestabilne i mieć nieliniową fazę. Ten symulator oferuje projekty FIR z oknem sinc oraz IIR Butterwortha i Czebyszewa.
Co właściwie robią poszczególne elementy sterujące?
Przyciski typu filtru wybierają filtr dolnoprzepustowy, górnoprzepustowy, pasmowoprzepustowy lub pasmowozaporowy. Przyciski metody projektowania przełączają między FIR z oknem sinc, IIR Butterwortha i IIR Czebyszewa. Suwak rzędu filtru ustawia stromość, suwaki częstotliwości granicznej ustawiają częstotliwości przejścia, a suwak tętnień (tylko dla Czebyszewa) ustawia dopuszczalne tętnienia w paśmie przepustowym w dB. Przyciski sygnału testowego podają na wejście filtru sinusoidę 440 Hz, szum biały lub świergot od 20 Hz do 8 kHz.
Jak wyrażana jest częstotliwość graniczna?
Wewnętrznie częstotliwość graniczna jest wartością znormalizowaną między 0,05 a 0,45, gdzie 0,5 odpowiada częstotliwości Nyquista (połowie częstotliwości próbkowania 44,1 kHz). Wyświetlacz przelicza tę wartość na herce, mnożąc przez 44100, więc znormalizowana wartość 0,10 odpowiada w przybliżeniu 4410 Hz. Ta znormalizowana konwencja jest standardem w DSP, ponieważ matematyka projektowa jest niezależna od rzeczywistej częstotliwości próbkowania.
Jak działa metoda FIR z oknem sinc?
Idealny filtr dolnoprzepustowy ma nieskończenie długą odpowiedź impulsową w kształcie funkcji sinc. Metoda z oknem sinc obcina ją do N+1 odczepów i mnoży przez okno Hamminga (0,54 minus 0,46 cosinus), aby stłumić tętnienia wynikające z gwałtownego obcięcia. Wersje górnoprzepustowa, pasmowoprzepustowa i pasmowozaporowa są budowane poprzez inwersję widmową oraz dodawanie lub odejmowanie prototypów dolnoprzepustowych, co daje filtr o liniowej fazie, zawsze stabilny.
Co mówi mi rozkład zer i biegunów?
Wykres na płaszczyźnie Z pokazuje pierwiastki licznika (zera, okręgi) i mianownika (bieguny, krzyżyki) funkcji przejścia. Zera obniżają charakterystykę amplitudową w pobliskich częstotliwościach; bieguny ją podnoszą. Aby filtr był stabilny, każdy biegun musi leżeć ściśle wewnątrz okręgu jednostkowego. Filtry FIR nie mają biegunów poza początkiem układu współrzędnych, dlatego są bezwarunkowo stabilne.
Dlaczego rząd filtru IIR jest ograniczony do 8, a FIR dozwolony jest do 64?
Filtry IIR osiągają strome zbocze przy bardzo małej liczbie współczynników, więc rząd 8 daje już ostre przejście, a kaskady IIR wysokiego rzędu stają się numerycznie delikatne. Filtry FIR potrzebują znacznie więcej odczepów dla tej samej stromości, ponieważ każdy odczep wnosi tylko niewielki wkład, dlatego narzędzie dopuszcza do 64 odczepów. Przełączenie z FIR na metodę IIR automatycznie ogranicza rząd do 8.
Czym są tętnienia Czebyszewa i po co miałbym ich chcieć?
Filtr Czebyszewa typu I zamienia kontrolowaną ilość tętnień w paśmie przepustowym na bardziej strome przejście niż filtr Butterwortha tego samego rzędu. Suwak tętnień ustawia tę wartość w decybelach, od 0,1 do 3,0 dB. Większe tętnienia kupują ostrzejsze zbocze; projekt Butterwortha daje natomiast maksymalnie płaskie pasmo przepustowe z łagodniejszym zboczem.
Czy ta symulacja jest fizycznie dokładna?
Tak, w granicach narzędzia edukacyjnego. Charakterystyka częstotliwościowa jest dokładną dyskretną transformatą Fouriera współczynników, sekcje biquad wykorzystują standardowe projekty transformaty biliniowej, a filtrowanie wykorzystuje strukturę bezpośrednią II transponowaną. Położenia zer i biegunów pochodzą z prawdziwego solvera pierwiastków wielomianu. Podczas przesuwania suwaków stosowane jest niewielkie wygładzanie wizualne, ale leżące u podstaw liczby są poprawnym DSP.
Gdzie w rzeczywistym świecie stosuje się filtry cyfrowe?
Są wszędzie w nowoczesnej elektronice: korektory i zwrotnice audio, redukcja szumów i usuwanie echa w rozmowach telefonicznych, filtrowanie antyaliasingowe przed konwersją analogowo-cyfrową, oczyszczanie sygnałów EKG i EEG w urządzeniach medycznych oraz wybór kanału w radiach definiowanych programowo. Kodeki takie jak MP3 wykorzystują banki filtrów do podziału dźwięku na podpasma, a szybkie filtry FIR działają na dedykowanym sprzęcie FPGA lub DSP.