e^{i2πf} zeichnet den Einheitskreis nach
Eine komplexe Exponentialfunktion e^{iθ} liegt auf dem Einheitskreis in der komplexen Ebene. Wenn θ zunimmt, rotiert der Punkt gegen den Uhrzeigersinn.
Bei einem digitalen Filter, der zu ganzzahligen Zeiten n = 0, 1, 2, … abgetastet wird, nimmt die Eigenfunktion e^{i2πfn} bei jedem Sample einen Schritt um den Winkel 2πf um den Kreis herum.
Frequenz als Rotationsgeschwindigkeit: f misst, wie viel einer vollständigen Umdrehung pro Sample auftritt.
- f = 0: keine Rotation; Punkt bleibt bei (1, 0)
- f = 1/4: Viertelrotation bei jedem Schritt
- f = 1/2: halbe Rotation bei jedem Schritt (Nyquist-Frequenz)
- f = 1: vollständige Rotation bei jedem Schritt — nicht zu unterscheiden von f = 0
Dieser letzte Punkt enthält die ganze Aliasing-Geschichte geometrisch.
Warum der Einheitskreis
Der Einheitskreis ist die Menge {z : |z| = 1}. Evaluierung der Z-Transformation H(z) auf dem Einheitskreis — setzen z = e^{i2πf} — gibt die Frequenzantwort H(f). Der Einheitskreis ist die Grenze, wo zeitdiskrete Stabilität & Frequenzanalyse aufeinandertreffen.
Winkel & Frequenzen
Jede Frequenz f entspricht einem Winkel θ = 2πf Radianten pro Sample. Der volle Bereich unterschiedlicher Frequenzen umfasst eine vollständige Umdrehung: f ∈ [0, 1) oder äquivalent θ ∈ [0, 2π).
Bei der Nyquist-Frequenz f = 1/2 rückt jedes Sample genau π Radianten vor — eine halbe Umdrehung.
Das geometrische Bild des Aliasing
Der Einheitskreis hat einen Umfang von 2π. Eine vollständige Umdrehung entspricht Frequenz f = 1 (ein vollständiger Zyklus pro Sample). Die unterschiedlichen Frequenzen in einem abgetasteten Signal belegen genau eine Umdrehung.
Was passiert bei f = 1/2 + δ? Die Rotation pro Sample = 2π(1/2 + δ) = π + 2πδ. Nach k Samples ist der Winkel = k(π + 2πδ). Aber Winkel π + 2πδ ist geometrisch identisch mit −π + 2πδ, was der Rotation der Frequenz f = 1/2 − δ entspricht.
Aliasing ist modulare Arithmetik auf dem Kreis. Frequenzen oberhalb der Nyquist-Frequenz umbrechen. Der Kreis hat keine Erinnerung daran, von welcher Richtung sie kamen.
Der Abtastsatz besagt: Bleiben Sie im Halbkreis [0, π). Tasten Sie schnell genug ab, damit Ihr Signal die andere Hälfte nie erreicht. Anti-Aliasing-Filter erzwingen diese Grenze, bevor das Signal den Abtaster erreicht.
Berechnung von Aliasen geometrisch
Der Alias der Frequenz f bei Abtastrate f_s erscheint bei |f − round(f / f_s) · f_s| — die Entfernung zum nächsten Vielfachen von f_s, ausgedrückt als Bruchteil.
Für f_s = 1 (normalisiert): Alias von f = 1 − f für f ∈ (1/2, 1). Dies ist die Spiegelung von f um den Nyquist-Punkt f = 1/2.
Geometrisch: f & 1 − f sitzen an spiegelbildlichen Positionen auf dem Einheitskreis, gleich weit entfernt von der π-Achse.
Magnitudenreaktion als Distanzprodukt
Für eine Übertragungsfunktion H(z) mit Nullstellen z_1, z_2, … und Polen p_1, p_2, …:
|H(f)| = (∏ |e^{i2πf} − z_k|) / (∏ |e^{i2πf} − p_k|)
Dies ist die grafische Methode zum direkten Ablesen der Frequenzantwort aus dem Pol-Nullstellen-Diagramm.
Regeln:
- Eine Nullstelle AUF dem Einheitskreis erzeugt an dieser Frequenz eine perfekte Auslöschung.
- Ein Pol NAHE dem Einheitskreis erzeugt einen Spitzenwert in der Antwort.
- Eine Nullstelle nahe dem Einheitskreis (aber nicht auf ihm) erzeugt eine Absenkung, nicht eine Auslöschung.
- Pole INNERHALB des Einheitskreises halten den Filter stabil.
Die Z-Ebenen-Geometrie kodiert das gesamte Filterverhalten visuell. Ingenieure skizzieren Pol-Nullstellen-Diagramme, bevor sie Koeffizienten berechnen.