FilterLab šum a odhad stavu

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: okna a lineární fáze

Klouzavý průměr je FIR filtr, ve kterém mají všechny koeficienty stejnou váhu. Když se ta váha smí měnit, dá se charakteristika navrhnout skoro libovolně — a stojí to jen tolik násobení, kolik je koeficientů.

y[k] = Σi=0N−1 h[i]·x[k−i]konečná impulzní odezva

Proč vůbec okno

Ideální dolní propust má obdélníkovou charakteristiku, a její impulzní odezva je funkce sin(x)/x — nekonečně dlouhá. Uříznout ji na N koeficientů je nutné a je to zároveň ta chyba, kterou pak řeší okno: ostré uříznutí se ve frekvenci projeví jako zvlnění a postranní laloky, které při −21 dB ustrnou a víc koeficientů je nezlepší. Tomu se říká Gibbsův jev.

Okno je násobení koeficientů hladkou obálkou, která je na krajích nulová. Laloky spadnou — a zaplatí se to širším přechodovým pásmem.

0.501.02.05.0102050-75-50-250.0−6 dBfcfrekvence [Hz]dBkoeficienty
Okno: Hamming, 31 koeficientů. Nejhorší postranní lalok -52,4 dB, šum × 0,365, zpoždění 150 ms = (N−1)/2 vzorků, na každé frekvenci stejně.
Přepni okno na „žádné" a podívej se, co se stane s postranními laloky: ustrnou na −21 dB a víc koeficientů je nezlepší. Okno je vymění za širší přechod — a to je jediná věc, kterou okno dělá.
Přepni okno na „žádné" a přidávej koeficienty: přechod se zužuje, ale první postranní lalok zůstane trvale na −21 dB. Pak přepni na Blackmanovo a podívej se, co se stalo s oběma čísly.
OknoPostranní lalokyŠířka přechodu
žádné (obdélník)−21 dBnejužší
Hann−31 dB
Hamming−41 dB
Blackman−57 dB

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

Je-li posloupnost koeficientů symetrická, je zpoždění filtru přesně (N−1)/2 vzorku na každé frekvenci. Signál se tedy posune, ale nezdeformuje: složky si zachovají vzájemné časování. To je vlastnost, kterou IIR filtr mít nemůže, a je to důvod, proč se FIR používá tam, kde se z průběhu odečítají časy — v měření, v analýze, ve zpracování zvuku.

Cena je délka. Ostrý přechod potřebuje hodně koeficientů a každý koeficient je půl vzorku zpoždění. FIR s 200 koeficienty při 1 kHz zpozdí o 100 ms. V regulační smyčce je to konec; v offline analýze nikoho netrápí, protože tam se dá filtrovat i pozpátku a zpoždění úplně zrušit.

Šumové zesílení FIR filtru je √(Σh²) — a nezávisí to na tom, jak filtr vypadá ve frekvenci. Dva filtry se stejnou mezí −3 dB můžou propustit znatelně jiný šum. Kalkulačka proto ukazuje potlačení šumu a zpoždění vedle sebe, nikdy jen jedno z nich.

Vzorce v této kapitole

h[i] — okenní sinc
h[i] = w[i]·sin(ω_c(i−M))/(π(i−M)) [—] návrh oknem
t_d — zpoždění symetrického FIR
t_d = (N−1)/2 · Δt [s] lineární fáze
g_n — šumové zesílení FIR
g = √(Σ h[i]²) [—] Parsevalova věta