Återkopplingsstrukturen
Ett FIR-filter (finite impulse response) beräknar varje utgångsprove som en viktad summa av nuvarande och tidigare ingångar endast. Ingen återkoppling. Impulsresponsen har begränsad varaktighet.
Ett IIR-filter (infinite impulse response) för tidigare utgångar tillbaka in i beräkningen:
y_n = Σ b_k · x_{n−k} − Σ a_k · y_{n−k}
Den rekursiva termen −Σ a_k · y_{n−k} skapar återkoppling. En enda impuls vid ingången kommer att ekas runt återkopplingsslingan på obestämd tid (avklingar geometriskt om den är stabil).
Varför använda återkoppling?
Ett IIR-filter kan uppnå skarp frekvensselektivitet med mycket färre koefficienter än ett FIR-filter. En IIR med 2 poler kan approximera vad som kräver ett FIR med 50 koefficienter för samma stoppbandsdämpning.
Priset: potentiell instabilitet. Polerna för H(z) bestämmer stabiliteten. Alla poler måste ligga strikt innanför enhetscirkeln.
Hammings återkopplingsduschistoria
Hamming använde en levande personlig historia för att illustrera återkopplingsinstabilitet.
Han stannade på samma hotellrum upprepade gånger eftersom förtrogenhet hjälpte honom att orientera sig när han var trött. Rörmokaren hade installerat rör med stor diameter för varmvatten i duschen. Dessa skapade en betydande fördröjning mellan att vrida på kranen och att känna temperaturförändringen på vattnet.
Varje morgon följde Hamming samma mönster: vattnet för kallt → höj upp varmvattnet → fortfarande kallt → höj upp mer → plötsligt kokande → hoppa ut → sänk ned → upprepa.
Fördröjningen i återkopplingsslingan innebar att hans korrigeringar alltid gick för långt. Han kunde inte anpassa sig till fördröjningen, även efter många upprepningar.
Ingenjörsläxan: instabilitet uppstår antingen från för stor förstärkning i återkopplingsslingan ELLER för stor fördröjning i återkopplingsslingan. Båda manifesteras som samma jaktbeteende. I filterbetingelser: poler på eller utanför enhetscirkeln producerar exakt detta oscillerande eller divergerande svar.
Karakterisering av instabiliteten
Hamming observerade att samma duschistabilitet kunde analyseras på två sätt:
1. Hans respons var för stark (för stor förstärkning i korrigeringsåtgärden).
2. Hans detektion var för försenad (för bråttom att gå in i badkaret innan systemet stabiliserades).
Båda beskrivningarna producerar samma matematiska resultat: återkopplingsslingans pol har flyttat sig utanför enhetscirkeln.
De fyra klassiska familjerna
Analog filterteori utvecklades omkring fyra klassiska designfamiljer, var och en representerar en annan avvägning. Dessa familjer transformeras till diskret tid via den bilineära transformationen eller impulsiv invarians.
Butterworth (maximalt platt)
Passbandrespons: |H(jω)|² = 1 / (1 + (ω/ω_c)^{2N}). Monotont minskande. Ingen krusning i passband eller stoppband. Poler ligger på en cirkel med radie ω_c i s-planet (eller transformerad cirkel i z-planet). Det flataste möjliga passbandet för en given ordning N.
Chebyshev typ I
Lika krusning i passbandet, monotont i stoppbandet. För en given ordning N och krusningsnivå uppnår skarpare avskärning än Butterworth. Poler ligger på en ellips (i s-planet).
Chebyshev typ II
Lika krusning i stoppbandet, monotont i passbandet. Speglad bild av typ I i frekvensdomänen.
Elliptisk (Cauer)
Lika krusning i BÅDA passbandet och stoppbandet. För en given ordning N och krusningsnivåer uppnår den skarpaste möjliga övergången från passband till stoppband. Använder elliptiska funktioner för att placera poler & nollställen optimalt. Hamming: namnet kommer från det faktum att elliptiska funktioner används i härledningen.
Den grundläggande avvägningen
Alla fyra familjerna uppnår samma grundläggande avvägning på olika sätt: högre ordning N ger skarpare övergång. Att tillåta krusning (Chebyshev, elliptisk) uppnår skarpare övergång för samma N. Elliptisk uppnår den absolut skarpaste övergången för vilken given ordning N och krusningsspecifikationer som helst.
Att välja bland filterfamiljer
Valet mellan familjer beror på vad applikationen kan tolerera.
Att ifrågasätta expertpåståndet
Hamming mindes att vissa experter hade hävdat att alla IIR-filter (rekursiva filter) ägde en viss egenskap. Han frågade sig själv om detta verkligen var sant — och hittade ett motexempel.
Hans poäng: experter bär ofta påstår som de absorberade i skolan utan att någonsin återvärdera dem med tanke på nuvarande problem. Om du frågar dig själv om det du får höra verkligen är sant, är det häpnadsväckande hur mycket som faktiskt är, eller ligger nära, falskt, även i väletablerade områden.
Motexemplet var inte den typ av filter du normalt skulle designa, men det visade att påståndet var ytligt. Ett enda motexempel är tillräckligt för att motbevisa ett universellt påstående.
IIR-design i praktiken
Hamming noterade att han självständigt hade utvecklat mycket av IIR-filterteorin medan han löste ett annat problem: att härleda stabila korrektörformler för numeriska vanliga differentialekvationer.
Korrektörformelns form: y_n = Σ a_k · y_{n−k} + Σ b_k · f(y_{n−k})
Återkoppling förekommer både i y-termerna (linjär återkoppling) och i f(y)-termerna (olinjär återkoppling genom differentialekvationen). Stabilitet för IIR-filter är ett specialfall av det mer allmänna problemet med stabilitet för numeriska ODE-integraler.
Att ansluta återkoppling över domäner
Samma matematiska struktur — återkoppling, poler, stabilitetsgräns — förekommer i digitala filter, numeriska ODE-lösare, kontrollsystem, biologiska rytmer, & ekonomiska modeller.
I varje domän: en återkopplingslinga beräknar ett nytt tillstånd från tidigare tillstånd. Stabilitet kräver att återkopplingen inte förstärker störningar på obestämd tid.
Enhetscirkelns stabilitetsgräns i Z-planet motsvarar: den imaginära axeln i Laplace s-planet (kontinuerlig tid), spektralradiusvillkoret ρ(A) < 1 för linjära iterationer, & Lyapunov-exponentvillkoret λ < 0 för olinjära system.