Die Rückkopplungsstruktur
Ein FIR-Filter (Finite Impulse Response) berechnet jedes Ausgabesample als gewichtete Summe von aktuellen & vorherigen Eingaben nur. Keine Rückkopplung. Die Impulsantwort hat endliche Dauer.
Ein IIR-Filter (Infinite Impulse Response) führt vorherige Ausgaben zurück in die Berechnung:
y_n = Σ b_k · x_{n−k} − Σ a_k · y_{n−k}
Der rekursive Term −Σ a_k · y_{n−k} erzeugt Rückkopplung. Ein einzelner Impuls am Eingang wird auf unbestimmte Zeit um die Rückkopplungsschleife herum widerhallen (geometrisch abfallend, wenn stabil).
Warum Rückkopplung verwenden?
Ein IIR-Filter kann scharfe Frequenzselektivität mit viel weniger Koeffizienten erreichen als ein FIR-Filter. Ein 2-Pol-IIR kann approximieren, was ein 50-Koeffizient-FIR für die gleiche Sperrband-Dämpfung erfordert.
Der Preis: mögliche Instabilität. Die Pole von H(z) bestimmen die Stabilität. Alle Pole müssen streng innerhalb des Einheitskreises liegen.
Hammings Rückkopplungs-Dusch-Geschichte
Hamming nutzte eine lebhafte persönliche Geschichte, um Rückkopplungsinstabilität zu veranschaulichen.
Er blieb wiederholt im selben Hotelzimmer, weil Vertrautheit ihm half, sich zu orientieren, wenn er müde war. Der Klempner hatte Rohre mit großem Durchmesser für heißes Wasser in der Dusche installiert. Diese erzeugten eine erhebliche Verzögerung zwischen der Verstellung des Griffs und dem Spüren der Temperaturänderung.
Jeden Morgen folgte Hamming dem gleichen Muster: Wasser zu kalt → heiß aufdrehen → immer noch kalt → noch mehr aufdrehen → plötzlich kochend → rausspringen → herunterdrehen → wiederholen.
Die Verzögerung in der Rückkopplungsschleife bedeutete, dass seine Korrektionen immer überschossen. Er konnte sich nicht an die Verzögerung anpassen, auch nicht nach vielen Wiederholungen.
Die Ingenieurlektion: Instabilität entsteht entweder durch übermäßige Verstärkung in der Rückkopplungsschleife ODER übermäßige Verzögerung in der Rückkopplungsschleife. Beide manifestieren sich als das gleiche Jagdverhalten. In Filterausdrücken: Pole am oder außerhalb des Einheitskreises produzieren genau diese oszillierende oder divergierende Antwort.
Charakterisierung der Instabilität
Hamming beobachtete, dass die gleiche Dusch-Instabilität auf zwei Wegen analysiert werden konnte:
1. Seine Antwort war zu stark (übermäßige Verstärkung in der Korrekturmaßnahme).
2. Seine Erkennung war zu verzögert (zu voreilig, um die Wanne zu betreten, bevor sich das System beruhigte).
Beide Beschreibungen führen zum gleichen mathematischen Ergebnis: Der Pol der Rückkopplungsschleife hat sich außerhalb des Einheitskreises bewegt.
Die vier klassischen Familien
Die Analogfilter-Theorie entwickelte sich um vier klassische Designfamilien, die jeweils einen anderen Kompromiss darstellen. Diese Familien transformieren in diskrete Zeit über die bilineare Transformation oder Impuls-Äquivalenz.
Butterworth (maximal flach)
Passbandantwort: |H(jω)|² = 1 / (1 + (ω/ω_c)^{2N}). Monoton abnehmend. Keine Welligkeit im Passband oder Sperrband. Pole liegen auf einem Kreis mit Radius ω_c in der s-Ebene (oder transformierter Kreis in der z-Ebene). Das flachste mögliche Passband für eine gegebene Ordnung N.
Chebyshev Typ I
Gleiche Welligkeit im Passband, monoton im Sperrband. Für eine gegebene Ordnung N und Welligkeit-Pegel erreicht einen schärferen Übergang als Butterworth. Pole liegen auf einer Ellipse (in der s-Ebene).
Chebyshev Typ II
Gleiche Welligkeit im Sperrband, monoton im Passband. Spiegelbild von Typ I im Frequenzbereich.
Elliptisch (Cauer)
Gleiche Welligkeit in BEIDEN Passband und Sperrband. Für eine gegebene Ordnung N und Welligkeit-Pegel erreicht den schärfstmöglichen Übergang vom Passband zum Sperrband. Verwendet elliptische Funktionen, um Pole & Nullstellen optimal zu platzieren. Hamming: Der Name kommt daher, dass elliptische Funktionen in der Herleitung verwendet werden.
Der grundlegende Kompromiss
Alle vier Familien erzielen den gleichen grundlegenden Kompromiss unterschiedlich: höhere Ordnung N gibt schärferen Übergang. Die Zulassung von Welligkeit (Chebyshev, elliptisch) erreicht einen schärferen Übergang für die gleiche N. Elliptisch erreicht den absolut schärfsten Übergang für jede gegebene Ordnung N und Welligkeit-Spezifikationen.
Auswahl zwischen Filterfamilien
Die Wahl zwischen Familien hängt davon ab, was die Anwendung toleriert.
Die Expertenbehauptung hinterfragen
Hamming erinnerte sich daran, dass bestimmte Experten behauptet hatten, dass alle IIR-Filter (rekursive) eine bestimmte Eigenschaft besäßen. Er fragte sich selbst, ob dies wirklich wahr war — und fand ein Gegenbeispiel.
Sein Punkt: Experten tragen oft Behauptungen, die sie in der Schule aufgesogen haben, ohne sie jemals im Kontext aktueller Probleme erneut zu prüfen. Wenn Sie sich selbst fragen, ob das, was Ihnen gesagt wird, wirklich wahr ist, ist es erstaunlich, wie viel Sie finden können, das ist, oder sich dem Falschen nähert, sogar in einem gut entwickelten Gebiet.
Das Gegenbeispiel war nicht die Art von Filter, die Sie normalerweise entwerfen würden, aber es bewies die Behauptung als oberflächlich. Ein einzelnes Gegenbeispiel genügt, um eine universelle Behauptung zu widerlegen.
IIR-Design in der Praxis
Hamming vermerkte, dass er viel der IIR-Filtertheorie unabhängig entwickelt hatte, während er ein anderes Problem löste: Ableitung stabiler Korrektorformeln für numerische gewöhnliche Differentialgleichungen.
Die Korrektorformel-Form: y_n = Σ a_k · y_{n−k} + Σ b_k · f(y_{n−k})
Rückkopplung erscheint sowohl in den y-Termen (lineare Rückkopplung) als auch in den f(y)-Termen (nichtlineare Rückkopplung durch die Differentialgleichung). Stabilität für IIR-Filter ist ein Spezialfall des allgemeineren Problems der Stabilität für numerische ODE-Integratoren.
Rückkopplung über Domänen hinweg verbinden
Die gleiche mathematische Struktur — Rückkopplung, Pole, Stabilitätsgrenze — erscheint in digitalen Filtern, numerischen ODE-Lösern, Kontrollsystemen, biologischen Rhythmen, & wirtschaftlichen Modellen.
In jeder Domäne: eine Rückkopplungsschleife berechnet einen neuen Zustand aus vorherigen Zuständen. Stabilität erfordert, dass die Rückkopplung Störungen nicht unbegrenzt verstärkt.
Die Stabilitätsgrenze des Einheitskreises in der Z-Ebene entspricht: der imaginären Achse in der Laplace-s-Ebene (kontinuierliche Zeit), der Spektralradiusbedingung ρ(A) < 1 für lineare Iterationen, & der Lyapunov-Exponenten-Bedingung λ < 0 für nichtlineare Systeme.