Vom Zeitbereich in den komplexen Raum
Die Z-Transform abbildet eine Folge x_n auf eine Funktion X(z) von einem komplexen Variablen:
X(z) = Σ_{n=−∞}^{∞} x_n · z^{−n}
Die Variable z parametrisiert den komplexen Raum. Verschiedene Regionen dieses Raumes entsprechen unterschiedlichen qualitativen Verhaltens des Filters.
Geometrische Regionen
| Region | |z| | Verhalten |
|--------|-----|---------|
| Innerer Einheitskreis | < 1 | Stabile Pole: abklingender Ansprech |
| Einheitskreis | = 1 | Frequenzachse: z = e^{i2πf} |
| Äusserer Einheitskreis | > 1 | Instabile Pole: wachsendes Verhalten |
Der Einheitskreis spielt bei der diskreten Zeitstabilität die gleiche Rolle wie die Imaginärachse bei der kontinuierlichen Zeit (Laplace) Stabilität.
Beziehung zur Laplace-Transform
Für kontinuierliche Systeme verwendet die Laplace-Transform die Variable s. Die Imaginärachse s = iω ist, wo die Frequenzantwort lebt. Stabilität: Polen müssen haben Re(s) < 0 (linker Halbraum).
Die bilineare Transformiert kartiert s → z: z = (1 + Ts/2) / (1 − Ts/2). Dies kartiert den linken Halbraum in den inneren Einheitskreis - die geometrische Übersetzung von 'linker Halbraum stabil' in 'innerer Einheitskreis stabil.'
Die bilineare Transform als konformes Abbild
Die bilineare Transform z = (1 + Ts/2) / (1 − Ts/2) ist eine Möbius-Transformation - eine konforme (Winkelbeibehaltende) Karte des komplexen Raumes.
Ihre Schlüsselgeometrischen Eigenschaften:
- Kartiert s = iω (Imaginärachse) auf |z| = 1 (Einheitskreis)
- Kartiert Re(s) < 0 (linker Halbraum) auf |z| < 1 (innerer Einheitskreis)
- Kartiert Re(s) > 0 (rechter Halbraum) auf |z| > 1 (äusserer Einheitskreis)
- Frequenzverzerrung: Die Abbildung ω → f ist nichtlineare - ω_analog = (2/T)·tan(πf_digital)
Diese Verzerrung komprimiert hohe Frequenzen in Richtung des Nyquist-Punkts. Entwickler berücksichtigen dies, indem sie die analogen Spezifikation vor der Anwendung des Bilineartransformationen vorverzerrt.
Butterworth-Pole: Kreis Locus
Butterworth-Filter erreichen einen maximal flachen Passband, indem sie analoge Pole auf einem Kreis mit dem Radius ω_c in der s-Ebene platzieren.
Für ein N-tes Ordnungs-Butterworth-Filter befinden sich die Pole an:
s_k = ω_c · e^{iπ(2k+N−1)/(2N)} für k = 0, 1, …, N−1
Dies platziert sie gleichmäßig auf der linken Hälfte eines Kreises mit dem Radius ω_c. (Pole auf der rechten Hälfte wären unstabil; nur die linkshälfteplanen-Pole werden beibehalten.)
Warum kreisförmiger Locus → maximal flacher Passband?
Das Butterworth-Polynom |B_N(jω)|² = 1 + (ω/ω_c)^{2N} hat alle seine Pole bei |s| = ω_c. Die gleichmäßige Radienbeschränkung bedeutet, dass alle Pole gleichmäßig zum Ausgleich der Lautstärke bei ω = ω_c beitragen. Das Maximums-Flachheits-Theorem: Unter allen N-ten Ordnungspolynomen mit Polen auf diesem Kreis hat das Butterworth-Polynom die meisten Ableiter, die in ω = 0 gleich null sind.
Chebyshev-Pole: Ellipse Locus
Chebyshev-Polen lie auf einer Ellipse im s-Ebenen (nicht eine Kreis). Die Ellipse hat Achsenlängen, die durch den Rissparameter ε bestimmt werden. Das gleichmäßige Rauschen im Passband ergibt sich aus der Gleichverteilungseigenschaft von Chebyshev-Polynomen.
Elliptische Polen: Elliptische Funktionen-Basis
Elliptische (Cauer)-Filter-Polen liegen ebenfalls auf einer Ellipse - aber mit BEOHREN Polen UND Nullstellen, die die Frequenzantwort beeinflussen. Die Nullstellen befinden sich auf der imaginären Achse (endliche Abschwächungspolen im Unterband). Die elliptische Funktion-Kartierung verteilt die Nullstellen optimal, um gleichmäßigen Rausch in beiden Bändern gleichzeitig zu erreichen.
Rechnen Sie die Butterworth-Polpositionen aus
Für einen 4. Ordnung Butterworth-Filter mit ω_c = 1 (normalisiert) liegen die Polen bei:
s_k = e^{iπ(2k+3)/8} für k = 0, 1, 2, 3
k=0: s₀ = e^{i3π/8} (in der linken Halbebene)
k=1: s₁ = e^{i5π/8} (in der linken Halbebene)
k=2: s₂ = e^{i7π/8} (in der linken Halbebene)
k=3: s₃ = e^{i9π/8} (in der linken Halbebene)
Diese vier Polen befinden sich gleichmäßig im Winkel auf dem Einheitskreis, und alle haben negative reelle Teile (linkes Halbräumchen).
Entfernung der Pole zur Einheitskreis
Theoretische Stabilität erfordert |p| < 1. In der Praxis ergeben sich zwei zusätzliche Bedenken.
Stabilitätsmarginal
Der Stabilitätsmarginal eines IIR-Filter ist der minimale Abstand von jeder Polstelle zur Einheitskreis: min_k (1 − |p_k|).
Ein Pol bei |p| = 0,99 ist technisch stabil, aber es bleibt nur ein Prozent Marginal. Die endliche Genauigkeit der Arithmetik (Rundung in der Koeffizientenrepräsentation & Akkumulation von Rundungsfehlern) kann Polstellen effektiv verschieben. Wenn die Koeffizientenquantisierung einen Pol von 0,99 auf 1,001 verschiebt, wird das Filter unwirsch.
Geometrische Folge
Pole sehr nahe an der Einheitskreis erzeugen sehr scharfe Frequenzantwortsspitzen - schmale Bandbreitenresonatoren. Aber schmale Resonatoren erfordern hohe Genauigkeit: kleine Koeffizientenfehler verschieben die Spitzenfrequenz signifikant.
Der geometrische Handel: Spitzenpräzision ∝ 1 / (1 − |p|). Als |p| → 1, → Spitzenpräzision ∞, aber Stabilitätsmarginal → 0 und Empfindlichkeit gegenüber Koeffizientenfehlern → ∞.
Zweiten-Ordnungsabschnitte
Ein hochrangiges IIR-Filter, das als einzelnes Polynom implementiert ist, ist numerisch empfindlich - Rundung eines einzelnen Koeffizienten kann viele Pole verschieben. Die Standardlösung: Implementieren Sie es als Kaskade von Zweiten-Ordnungsabschnitten (Biquads), jeder mit nur einem konjugierten Polpaar und einem konjugierten Nullpaar. Fehler in einem Abschnitt können Polstellen in anderen nicht beeinflussen.