FourierLab z čeho se to skládá

Tohle je statická kopie kapitoly pro vyhledávače. Interaktivní verze má animované obrázky, kontrolní otázky a tlačítka, která příklad načtou do kalkulačky.

FIR filtr: ideál se musí useknout

Ideální dolní propust má ve frekvenci obdélník: do mezní frekvence propustí všechno beze změny, nad ní nic. Zbývá zjistit, jaká impulsní odezva takové spektrum má — a to je zpětná transformace obdélníku:

h(t) = 2fc · sinc(2fc·t)ideální dolní propust

Jenže ten sinc je nekonečně dlouhý a začíná před nulou. Dvě věci naráz, obě zabíjející: nekonečně mnoho koeficientů se nespočítá, a odezva před příchodem impulsu znamená, že filtr musí znát budoucnost. Ideální filtr tedy neexistuje — ne že by byl drahý, on je nekauzální.

Reálný filtr je useknutý sinc. A jak se usekne, tím se rozhoduje o všem ostatním.

impulsní odezva — useknutý sinc-80-60-40-2000100200300400500Hzfrekvenční charakteristika [dB]Hannkoeficientů41−3 dB naměřeno89,8 Hznepropustné pásmo−62,5 dBzpoždění20,00 msa je stejné provšechny frekvence
Přechod je 61 Hz široký. Dvojnásobek koeficientů by dal 29 Hz — přechod se zužuje nepřímo úměrně počtu koeficientů, ne rychleji.
Nahoře koeficienty (useknutý sinc je vidět), dole výsledná charakteristika. Přidávej koeficienty a sleduj, jak se přechod zužuje.

Co se čím platí

chcizaplatímnaměřeno
ostřejší přechodvíc koeficientů21 → 125 Hz, 161 → 15,6 Hz (8× víc koef., 8× užší)
hlubší útlumširší přechodHann −68 dB, Blackman −80 dB při týchž 51 koeficientech
menší zpožděnímíň koeficientů, horší filtrzpoždění = (N−1)/2 vzorků

Přechod se zužuje nepřímo úměrně počtu koeficientů — osmkrát víc koeficientů dá osmkrát užší přechod, ne šestnáctkrát. To je tvrdá mez a žádný chytřejší návrh ji neobejde; posunout se dá jen konstanta.

Lineární fáze: to, proč se FIR používá

h[n] symetrická ⇒ zpoždění (N−1)/2 stejné pro VŠECHNY frekvencelineární fáze

Symetrická impulsní odezva zpozdí všechny frekvence stejně. To znamená, že tvar signálu zůstane zachovaný — jen přijde později. Filtr, který tuhle vlastnost nemá, deformuje průběh, i když amplitudy nechá na pokoji: hrany se rozmažou a impulsy roztečou. Pro zvuk, obraz i měřicí řetězec je to zásadní, a je to hlavní důvod, proč se FIR volí i tam, kde by IIR filtr stačil s desetinou koeficientů.

Testy v repu ověřují symetrii na 1e-15 a měří jednotkový stejnosměrný přenos přesně na 1 — filtr, který má propustit konstantu, ji nesmí zeslabit.

A na skutečném signálu. Dva tóny, 50 Hz a 300 Hz, dolní propust se 101 koeficienty a mezí 100 Hz: padesátka projde beze změny, tři sta klesne na 0,00 % původní amplitudy. Měřeno spektrem vstupu a výstupu, ne odhadem.

Vzorce v této kapitole

sinc — ideální dolní propust
h(t) = 2fc·sinc(2fc·t) [—] kap. 8
lin. fáze — zpoždění FIR filtru
(N − 1)/2 vzorků, pro všechny frekvence [vzorky] kap. 8