Od domeny czasu do płaszczyzny zespolonej
Transformata Z odwzorowuje ciąg x_n na funkcję X(z) zmiennej zespolonej:
X(z) = Σ_{n=−∞}^{∞} x_n · z^{−n}
Zmienna z parametryzuje płaszczyznę zespoloną. Różne regiony tej płaszczyzny odpowiadają różnym jakościowym zachowaniom filtru.
Regiony geometryczne
| Rejon | z | Zachowanie | ||
|---|---|---|---|---|
| Wewnątrz koła jednostkowego | < 1 | Bieguny stabilne: odpowiedź zanikająca | ||
| Koło jednostkowe | = 1 | Oś częstości: z = e^{i2πf} | ||
| Poza kołem jednostkowym | > 1 | Bieguny niestabilne: odpowiedź rosnąca |
Koło jednostkowe odgrywa taką samą rolę w stabilności czasu dyskretnego, jak oś urojona w stabilności czasu ciągłego (Laplace'a).
Związek z transformatą Laplace'a
W systemach czasu ciągłego transformata Laplace'a używa zmiennej s. Oś urojona s = iω jest miejscem, gdzie żyje odpowiedź częstościowa. Stabilność: bieguny muszą mieć Re(s) < 0 (lewa półpłaszczyzna).
Transformata biliniowa odwzorowuje s → z: z = (1 + Ts/2) / (1 − Ts/2). To odwzorowuje lewą półpłaszczyznę na wnętrze koła jednostkowego — geometryczne tłumaczenie 'stabilne w lewej półpłaszczyźnie' na 'stabilne wewnątrz koła jednostkowego'.
Transformata biliniowa jako mapa konforemna
Transformata biliniowa z = (1 + Ts/2) / (1 − Ts/2) jest transformacją Möbiusa — konforemną (zachowującą kąty) mapą płaszczyzny zespolonej.
Jej kluczowe właściwości geometryczne:
- Odwzorowuje s = iω (oś urojona) na |z| = 1 (koło jednostkowe)
- Odwzorowuje Re(s) < 0 (lewa półpłaszczyzna) na |z| < 1 (wnętrze koła jednostkowego)
- Odwzorowuje Re(s) > 0 (prawa półpłaszczyzna) na |z| > 1 (zewnętrze koła jednostkowego)
- Zniekształcenie częstości: odwzorowanie ω → f jest nieliniowe — ω_analog = (2/T)·tan(πf_digital)
To zniekształcenie kompresuje wysokie częstości w kierunku punktu Nyquista. Projektanci uwzględniają to poprzez wstępne zniekształcenie specyfikacji analogowej przed zastosowaniem transformaty biliniowej.
Bieguny Butterwortha: locus koła
Filtry Butterwortha osiągają maksymalnie płaskie pasmo przepustowe umieszczając bieguny analogowe na okręgu o promieniu ω_c w płaszczyźnie s.
Dla N-tego rzędu filtru Butterwortha bieguny siedzą na:
s_k = ω_c · e^{iπ(2k+N−1)/(2N)} dla k = 0, 1, …, N−1
To umieszcza je równomiernie rozmieszczone na lewej połowie okręgu o promieniu ω_c. (Bieguny na prawej połowie byłyby niestabilne; zachowywane są tylko bieguny na lewej półpłaszczyźnie.)
Dlaczego locus koła → maksymalnie płaskie pasmo przepustowe?
Wielomian Butterwortha |B_N(jω)|² = 1 + (ω/ω_c)^{2N} ma wszystkie bieguny przy |s| = ω_c. Ograniczenie równego promienia oznacza, że wszystkie bieguny przyczyniają się jednakowo do odpowiedzi amplitudowej przy ω = ω_c. Twierdzenie o maksymalnej płaskości: spośród wszystkich wielomianów N-tego rzędu z biegunami na tym okręgu, wielomian Butterwortha ma najwięcej pochodnych równych zeru przy ω = 0.
Bieguny Czebyszewa: locus elipsy
Bieguny Czebyszewa leżą na elipsie w płaszczyźnie s (nie na okręgu). Elipsa ma półosie główną i małą określone przez parametr tętnień ε. Równomiernie tętniejące pasmo przepustowe wynika z właściwości równooscylacyjnej wielomianów Czebyszewa.
Bieguny eliptyczne: locus funkcji eliptycznej
Bieguny filtru eliptycznego (Cauera) również leżą na elipsie — ale zarówno bieguny JAK I zera przyczyniają się do odpowiedzi częstościowej. Zera siedzą na osi urojonej (skończone bieguny tłumienia w paśmie zatrzymania). Odwzorowanie funkcji eliptycznej optymalnie rozmieszcza zera, aby osiągnąć równomierne tętnienia w obu pasmach jednocześnie.
Obliczanie lokalizacji biegunów Butterwortha
Dla filtru 4-tego rzędu Butterwortha z ω_c = 1 (znormalizowany), bieguny siedzą na:
s_k = e^{iπ(2k+3)/8} dla k = 0, 1, 2, 3
k=0: s₀ = e^{i3π/8} (na lewej półpłaszczyźnie)
k=1: s₁ = e^{i5π/8} (na lewej półpłaszczyźnie)
k=2: s₂ = e^{i7π/8} (na lewej półpłaszczyźnie)
k=3: s₃ = e^{i9π/8} (na lewej półpłaszczyźnie)
Te cztery bieguny siedzą przy równomiernym rozmieszczeniu kątowym na okręgu jednostkowym, wszystkie z ujemnymi częściami rzeczywistymi (lewa półpłaszczyzna).
Odległość od biegunów do koła jednostkowego
Teoretycznie stabilność wymaga |p| < 1. W praktyce pojawiają się dwa dodatkowe obawy.
Margines stabilności
Margines stabilności filtru IIR to minimalna odległość od jakichkolwiek biegunów do koła jednostkowego: min_k (1 − |p_k|).
Biegun przy |p| = 0,99 jest teoretycznie stabilny, ale pozostawia tylko 1% marginesu. Arytmetyka o ograniczonej precyzji (zaokrąglanie w reprezentacji współczynników & akumulacja błędów zaokrąglenia) może efektywnie przesunąć bieguny. Jeśli kwantyzacja współczynników przesunie biegun z 0,99 na 1,001, filtr staje się niestabilny.
Konsekwencja geometryczna
Bieguny bardzo blisko koła jednostkowego dają bardzo ostre szczyty odpowiedzi częstościowej — rezsonatory o wąskim paśmie. Ale wąskie rezsonatory wymagają wysokiej precyzji: małe błędy współczynników znacznie przesuwają częstość szczytu.
Geometryczny kompromis: ostrość szczytu ∝ 1 / (1 − |p|). Gdy |p| → 1, ostrość → ∞ ale margines stabilności → 0 i czułość na błędy współczynników → ∞.
Sekcje drugiego rzędu
Filtr IIR wysokiego rzędu zaimplementowany jako jeden wielomian jest numerycznie wrażliwy — zaokrąglenie jednego współczynnika może przesunąć wiele biegunów. Standardowe rozwiązanie: implementuj jako kaskadę sekcji drugiego rzędu (biquady), każda zawierająca tylko jedną sprzężoną parę biegunów i jedną sprzężoną parę zer. Błędy w jednej sekcji nie mogą zaburzyć biegunów w innych.