Van Tijddomein naar Complexe Vlak
De Z-transformatie wijst een rij x_n toe aan een functie X(z) van een complexe variabele:
X(z) = Σ_{n=−∞}^{∞} x_n · z^{−n}
De variabele z parametriseert het complexe vlak. Verschillende regio's van dit vlak komen overeen met verschillende kwalitatieve gedragingen van het filter.
Geometrische Regio's
| Regio | z | Gedrag | ||
|---|---|---|---|---|
| Binnen eenheidscirkel | < 1 | Stabiele polen: afnemende respons | ||
| Eenheidscirkel | = 1 | Frequentieas: z = e^{i2πf} | ||
| Buiten eenheidscirkel | > 1 | Onstabiele polen: groeiende respons |
De eenheidscirkel speelt dezelfde rol in discrete-tijd stabiliteit als de imaginaire as in continue-tijd (Laplace) stabiliteit.
Relatie met Laplace-Transformatie
Voor systemen in continue tijd gebruikt de Laplace-transformatie variabele s. De imaginaire as s = iω is waar de frequentierespons zich bevindt. Stabiliteit: polen moeten Re(s) < 0 hebben (linker halfvlak).
De bilineaire transformatie wijst s → z toe: z = (1 + Ts/2) / (1 − Ts/2). Dit wijst het linker halfvlak toe aan het binnenste van de eenheidscirkel — de geometrische vertaling van 'linker halfvlak stabiel' naar 'binnenste eenheidscirkel stabiel'.
De Bilineaire Transformatie als Conforme Afbeelding
De bilineaire transformatie z = (1 + Ts/2) / (1 − Ts/2) is een Möbius-transformatie — een conforme (hoek-bewarende) afbeelding van het complexe vlak.
De belangrijkste geometrische eigenschappen:
- Wijst s = iω (imaginaire as) toe aan |z| = 1 (eenheidscirkel)
- Wijst Re(s) < 0 (linker halfvlak) toe aan |z| < 1 (binnenste eenheidscirkel)
- Wijst Re(s) > 0 (rechter halfvlak) toe aan |z| > 1 (buitenste eenheidscirkel)
- Frequentievervormning: de afbeelding ω → f is niet-lineair — ω_analog = (2/T)·tan(πf_digital)
Deze vervormning comprimeert hoge frequenties richting het Nyquist-punt. Ontwerpers houden hiermee rekening door de analoge specificatie vooraf te vervormen voordat ze de bilineaire transformatie toepassen.
Butterworth-Polen: Cirkellocusus
Butterworth-filters bereiken een maximaal vlak doorlaatband door analoge polen op een cirkel met straal ω_c in het s-vlak te plaatsen.
Voor een N-de orde Butterworth-filter bevinden de polen zich op:
s_k = ω_c · e^{iπ(2k+N−1)/(2N)} for k = 0, 1, …, N−1
Dit plaatst ze gelijkmatig verdeeld op de linker helft van een cirkel met straal ω_c. (Polen op de rechter helft zouden onstabiel zijn; alleen de polen in het linker halfvlak worden behouden.)
Waarom cirkellocusus → maximaal vlak doorlaatband?
De Butterworth-polynoom |B_N(jω)|² = 1 + (ω/ω_c)^{2N} heeft al zijn polen op |s| = ω_c. De gelijke-straal-beperking betekent dat alle polen gelijkelijk bijdragen aan de grootheid-respons bij ω = ω_c. De maximale-vlakheid-stelling: onder alle N-de orde polynomen met polen op deze cirkel heeft de Butterworth-polynoom de meeste afleidingen gelijk aan nul bij ω = 0.
Chebyshev-Polen: Ellips-Locusus
Chebyshev-polen liggen op een ellips in het s-vlak (geen cirkel). De ellips heeft halflange en halfkorte assen bepaald door de golfregelparameter ε. Het gelijkmatig-golfende doorlaatband ontstaat uit de equioscillatieke eigenschap van Chebyshev-polynomen.
Elliptische Polen: Elliptische-Functie-Locusus
Elliptische (Cauer)-filterpolen liggen ook op een ellips — maar met ZOWEL polen ALS nulpunten die bijdragen aan de frequentierespons. De nulpunten bevinden zich op de imaginaire as (eindige-demping polen in de stopband). De elliptische-functie-afbeelding verdeelt de nulpunten optimaal om gelijkmatig golfing in beide banden tegelijkertijd te bereiken.
Berekenen van Butterworth-Poollocaties
Voor een 4e-orde Butterworth-filter met ω_c = 1 (genormaliseerd), bevinden de polen zich op:
s_k = e^{iπ(2k+3)/8} for k = 0, 1, 2, 3
k=0: s₀ = e^{i3π/8} (in linker halfvlak)
k=1: s₁ = e^{i5π/8} (in linker halfvlak)
k=2: s₂ = e^{i7π/8} (in linker halfvlak)
k=3: s₃ = e^{i9π/8} (in linker halfvlak)
Deze vier polen bevinden zich op gelijke hoekafstand op de eenheidscirkel, allen met negatieve reële delen (linker halfvlak).
Afstand van Polen tot de Eenheidscirkel
Theoretische stabiliteit vereist |p| < 1. In de praktijk ontstaan er twee aanvullende problemen.
Stabiliteitsmarge
De stabiliteitsmarge van een IIR-filter is de minimale afstand van elke pool tot de eenheidscirkel: min_k (1 − |p_k|).
Een pool op |p| = 0,99 is technisch gezien stabiel maar laat slechts 1% marge over. Rekenkundige precisie met eindige precisie (afronding in coëfficiëntrepresentatie & accumulatie van afrondingsfouten) kan effectief polen verplaatsen. Als coëfficiëntenkwantisering een pool verschuift van 0,99 naar 1,001, wordt het filter onstabiel.
Geometrische Gevolg
Polen zeer dicht bij de eenheidscirkel produceren zeer scherpe frequentieresponspunten — resonatoren met smalle bandbreedte. Maar smalle resonatoren vereisen hoge precisie: kleine coëfficiëntenfouten verplaatsen de piekfrequentie aanzienlijk.
De geometrische afweging: piekscherpheid ∝ 1 / (1 − |p|). Als |p| → 1, scherpheid → ∞ maar stabiliteitsmarge → 0 en gevoeligheid voor coëfficiëntenfouten → ∞.
Seconde-Orde Secties
Een IIR-filter met hoge orde, geïmplementeerd als een enkele polynoom, is numeriek gevoelig — afronding van een enkele coëfficiënt kan veel polen verplaatsen. De standaardoplossing: implementeer als een cascade van seconde-orde secties (biquads), elk met slechts één conjugaatpoolpaar en één conjugaatpuntpaar. Fouten in de ene sectie kunnen polen in anderen niet verstoren.