e^{i2πf} przebiera okrąg jednostkowy
Funkcja wykładnicza zespolona e^{iθ} lives on the unit circle w płaszczyźnie zespolonej. W miarę wzrostu θ, punkt obraca się przeciwnie do ruchu wskazówek zegara.
Dla filtrujący cyfrowy próbkowany w czasach całkowitych n = 0, 1, 2, …, funkcja własna e^{i2πfn} wykonuje krok kąta 2πf wokół okręgu przy każdej próbce.
Częstotliwość jako szybkość obrotu: f mierzy, jak duża część pełnego obrotu występuje na próbkę.
- f = 0: brak obrotu; punkt pozostaje w (1, 0)
- f = 1/4: obrót o ćwierć na każdym kroku
- f = 1/2: obrót o połowę na każdym kroku (częstotliwość Nyquista)
- f = 1: pełny obrót na każdym kroku — nierozróżnialny od f = 0
Ten ostatni punkt zawiera całą historię aliasingu geometrycznie.
Dlaczego okrąg jednostkowy
Okrąg jednostkowy to zbiór {z : |z| = 1}. Obliczanie transformacji Z H(z) na okręgu jednostkowym — ustawiając z = e^{i2πf} — daje odpowiedź częstotliwościową H(f). Okrąg jednostkowy jest granicą, gdzie spotykają się stabilność dyskretnych czasów & analiza częstotliwości.
Kąty & częstotliwości
Każda częstotliwość f odpowiada kątowi θ = 2πf radiany na próbkę. Pełny zakres różnych częstotliwości obejmuje jeden pełny obrót: f ∈ [0, 1) lub równoważnie θ ∈ [0, 2π).
Przy częstotliwości Nyquista f = 1/2, każda próbka zaawansuje dokładnie π radiany — połowę obrotu.
Geometryczny obraz aliasingu
Okrąg jednostkowy ma obwód 2π. Pełny obrót odpowiada częstotliwości f = 1 (jeden pełny cykl na próbkę). Różne częstotliwości w sygnale próbkowanym zajmują dokładnie jeden obrót.
Co się dzieje przy f = 1/2 + δ? Obrót na próbkę = 2π(1/2 + δ) = π + 2πδ. Po k próbkach, kąt = k(π + 2πδ). Ale kąt π + 2πδ jest geometrycznie identyczny z −π + 2πδ, który odpowiada obrotowi częstotliwości f = 1/2 − δ.
Aliasing jest arytmetyką modularną na okręgu. Częstotliwości powyżej częstotliwości Nyquista zawijają się. Okrąg nie ma pamięci, z którą strony przyszły.
Twierdzenie próbkowania mówi: pozostań w półokręgu [0, π). Próbkuj wystarczająco szybko, aby Twój sygnał nigdy nie dosięgnął drugiej połowy. Filtry anty-aliasingowe wymuszają tę granicę przed sygnałem docierającym do próbnika.
Obliczanie aliasów geometrycznie
Alias częstotliwości f przy szybkości próbkowania f_s pojawia się przy |f − round(f / f_s) · f_s| — odległości do najbliższej wielokrotności f_s, wyrażonej jako ułamek.
Dla f_s = 1 (znormalizowana): alias f = 1 − f dla f ∈ (1/2, 1). To jest odbicie f względem punktu Nyquista f = 1/2.
Geometrycznie: f & 1 − f siedzą w położeniach lustrzanych na okręgu jednostkowym, równo oddalone od osi π.
Odpowiedź amplitudowa jako iloczyn odległości
Dla funkcji przesyłu H(z) z zerami z_1, z_2, … & biegunami p_1, p_2, …:
|H(f)| = (∏ |e^{i2πf} − z_k|) / (∏ |e^{i2πf} − p_k|)
To jest metoda graficzna do bezpośredniego odczytania odpowiedzi częstotliwościowej z wykresu zer i biegunów.
Reguły:
- Zero NA okręgu jednostkowym tworzy doskonałe zero przy tej częstotliwości.
- Biegun BLISKO okręgu jednostkowego tworzy szczyt w odpowiedzi.
- Zero blisko okręgu jednostkowego (ale nie na nim) tworzy wpadlinę, nie zero.
- Bieguny WEWNĄTRZ okręgu jednostkowego utrzymują filtr stabilnym.
Geometria płaszczyzny Z koduje całe zachowanie filtra wizualnie. Inżynierowie szkicują wykresy zer i biegunów przed obliczaniem współczynników.