De Feedback-structuur
Een FIR (eindige impulsrespons) filter berekent elk uitgangsstaal als een gewogen som van huidge & voorbije ingangen alleen. Geen feedback. De impulsrespons heeft eindige duur.
Een IIR (oneindige impulsrespons) filter voert vorige uitgangen terug in de berekening:
y_n = Σ b_k · x_{n−k} − Σ a_k · y_{n−k}
De recursieve term −Σ a_k · y_{n−k} creëert feedback. Een enkele impuls op de ingang zal oneindig om de feedbacklus echoën (geometrisch afnemend als stabiel).
Waarom Feedback Gebruiken?
Een IIR-filter kan scherpe frequentieselectiviteit bereiken met veel minder coëfficiënten dan een FIR-filter. Een 2-pool IIR kan benaderen wat een 50-coëfficiënt FIR voor dezelfde stophogattenuatie vereist.
De prijs: mogelijke instabiliteit. De polen van H(z) bepalen stabiliteit. Alle polen moeten strikt binnen de eenheidscirkel liggen.
Hammings Feedback Doucheverhaal
Hamming gebruikte een levendig persoonlijk verhaal om feedbackinstabiliteit te illustreren.
Hij bleef herhaaldelijk in dezelfde hotelkamer omdat vertrouwdheid hem hielp zich te oriënteren wanneer hij moe was. De loodgieter had dikke warm-waterleidingen in de douche geïnstalleerd. Deze creëerden een aanzienlijke vertraging tussen het aanpassen van de kraan en het voelen van de temperatuurverandering.
Elke ochtend volgde Hamming hetzelfde patroon: water te koud → warm opzetten → nog steeds koud → meer opzetten → plotseling kokend → eruit springen → afzetten → herhalen.
De vertraging in het feedbackpad betekende dat zijn correcties altijd te veel waren. Hij kon zich niet aan de vertraging aanpassen, zelfs na veel herhalingen.
De ingenieurslessen: instabiliteit ontstaat door overmatige versterking in het feedbackpad OF overmatige vertraging in het feedbackpad. Beide tonen zich via dezelfde oscillaties. In filterbegrippen: polen op of buiten de eenheidscirkel produceren precies deze oscillerende of divergerende respons.
Karakterisering van de Instabiliteit
Hamming observeerde dat dezelfde doucheinstabiliteit op twee manieren kon worden geanalyseerd:
1. Zijn reactie was te sterk (overmatige versterking in de correctiewerking).
2. Zijn waarneming was te vertraagd (hij ging te gauw in het bad voordat het systeem zich had gestabiliseerd).
Beide beschrijvingen produceren hetzelfde wiskundige resultaat: de pool van de feedbacklus is buiten de eenheidscirkel verplaatst.
De Vier Klassieke Families
Analoge filterthorie is ontwikkeld rond vier klassieke ontwerp families, die elk een ander afweging vertegenwoordigen. Deze families transformeren naar discrete tijd via de bilineaire transformatie of impulssinvariantie.
Butterworth (Maximaal Vlak)
Doorlaatband-respons: |H(jω)|² = 1 / (1 + (ω/ω_c)^{2N}). Monotoon afnemend. Geen rimpel in doorlaatband of stoptband. Polen liggen op een cirkel met straal ω_c in het s-vlak (of getransformeerde cirkel in het z-vlak). De vlakste mogelijke doorlaatband voor gegeven orde N.
Chebyshev Type I
Gelijke rimpel in de doorlaatband, monotoon in de stoptband. Voor gegeven orde N en rimpelniveau bereikt een scherper afsnijding dan Butterworth. Polen liggen op een ellips (in het s-vlak).
Chebyshev Type II
Gelijke rimpel in de stoptband, monotoon in de doorlaatband. Spiegelbeeld van Type I in het frequentiedomein.
Elliptisch (Cauer)
Gelijke rimpel in ZOWEL doorlaatband als stoptband. Voor gegeven orde N en rimpelniveaus bereikt de scherpst mogelijke overgang van doorlaatband naar stoptband. Gebruikt elliptische functies om polen & nullen optimaal te plaatsen. Hamming: de naam komt van het feit dat elliptische functies in de afleiding worden gebruikt.
De Fundamentele Afweging
Alle vier families bereiken dezelfde basisafweging op verschillende manieren: hogere orde N geeft scherpere overgang. Rimpel toestaan (Chebyshev, elliptisch) bereikt scherpere overgang voor dezelfde N. Elliptisch bereikt de absoluut scherpste overgang voor gegeven N en rimpelspecificaties.
Kiezen Tussen Filterfamilies
De keuze tussen families hangt af van wat de toepassing toestaat.
De Deskundigenclaim in Vraag Stellen
Hamming herinnerde zich dat bepaalde deskundigen hadden geclaimd dat alle IIR (recursieve) filters een bepaalde eigenschap hadden. Hij vroeg zich af of dit echt waar was — en vond een tegenvoorbeeld.
Zijn punt: deskundigen handhaven vaak stellingen die ze op school hebben geleerd, zonder ze ooit in het licht van huidge problemen opnieuw te onderzoeken. Als je jezelf afvraagt of wat je wordt verteld werkelijk waar is, ontdek je dat verbazingwekkend veel daarvan onwaar is, of daaraan grenst, zelfs in een goed ontwikkeld vakgebied.
Het tegenvoorbeeld was niet het soort filter dat je normaal zou ontwerpen, maar het bewees dat de bewering oppervlakkig was. Een enkel tegenvoorbeeld volstaat om een universele bewering te weerleggen.
IIR-ontwerp in de Praktijk
Hamming merkte op dat hij onafhankelijk veel van de IIR-filterthorie had ontwikkeld terwijl hij een ander probleem oploste: het afleiden van stabiele correctieformules voor numerieke gewone differentiaalvergelijkingen.
De correctieformule vorm: y_n = Σ a_k · y_{n−k} + Σ b_k · f(y_{n−k})
Feedback verschijnt zowel in de y-termen (lineaire feedback) als in de f(y)-termen (niet-lineaire feedback door de differentiaalvergelijking). Stabiliteit voor IIR-filters is een speciaal geval van het meer algemene probleem van stabiliteit voor numerieke ODE-integrators.
Feedback Over Domeinen Verbinden
Dezelfde wiskundige structuur — feedback, polen, stabiliteitgrens — verschijnt in digitale filters, numerieke opossers, besturingssystemen, biologische ritmes, & economische modellen.
In elk domein: een feedbacklus berekent een nieuwe status uit vorige statussen. Stabiliteit vereist dat de feedback verstoringen niet onbeperkt versterkt.
De stabiliteitgrens van de eenheidscirkel in het Z-vlak komt overeen met: de imaginaire as in het Laplace s-vlak (continue tijd), de spectraelradiusvoorwaarde ρ(A) < 1 voor lineaire iteraties, & de Lyapunov-exponentvoorwaarde λ < 0 voor niet-lineaire systemen.