Do Domínio do Tempo ao Plano Complexo
O Z-transform mapeia uma sequência x_n para uma função X(z) de uma variável complexa:
X(z) = Σ_{n=−∞}^{∞} x_n · z^{−n}
A variável z parametriza o plano complexo. Diferentes regiões deste plano correspondem a diferentes comportamentos qualitativos do filtro.
Regiões Geométricas
| Região | z | Comportamento | ||
|---|---|---|---|---|
| Dentro do círculo unitário | < 1 | Pólos estáveis: resposta decrescente | ||
| Círculo unitário | = 1 | Eixo de frequência: z = e^{i2πf} | ||
| Fora do círculo unitário | > 1 | Pólos instáveis: resposta crescente |
O círculo unitário desempenha o mesmo papel na estabilidade em tempo discreto que o eixo imaginário desempenha na estabilidade em tempo contínuo (Laplace).
Relação com a Transformada de Laplace
Para sistemas em tempo contínuo, a transformada de Laplace usa a variável s. O eixo imaginário s = iω é onde a resposta de frequência reside. Estabilidade: os pólos devem ter Re(s) < 0 (semiplano esquerdo).
A transformação bilinear mapeia s → z: z = (1 + Ts/2) / (1 − Ts/2). Isto mapeia o semiplano esquerdo para o interior do círculo unitário — a tradução geométrica de 'estável no semiplano esquerdo' para 'estável dentro do círculo unitário'.
A Transformação Bilinear como Mapa Conformal
A transformação bilinear z = (1 + Ts/2) / (1 − Ts/2) é uma transformação de Möbius — um mapa conformal (que preserva ângulos) do plano complexo.
Suas propriedades geométricas principais:
- Mapeia s = iω (eixo imaginário) para |z| = 1 (círculo unitário)
- Mapeia Re(s) < 0 (semiplano esquerdo) para |z| < 1 (dentro do círculo unitário)
- Mapeia Re(s) > 0 (semiplano direito) para |z| > 1 (fora do círculo unitário)
- Distorção de frequência: o mapeamento ω → f é não-linear — ω_analog = (2/T)·tan(πf_digital)
Esta distorção comprime frequências altas em direção ao ponto de Nyquist. Os projetistas contabilizam isso pré-deformando a especificação analógica antes de aplicar a transformação bilinear.
Pólos de Butterworth: Locus Circular
Filtros de Butterworth alcançam banda passante maximalmente plana colocando pólos analógicos em um círculo de raio ω_c no plano s.
Para um filtro de Butterworth de N-ésima ordem, os pólos estão em:
s_k = ω_c · e^{iπ(2k+N−1)/(2N)} para k = 0, 1, …, N−1
Isto os coloca igualmente espaçados no semiplano esquerdo de um círculo de raio ω_c. (Pólos na metade direita seriam instáveis; apenas os pólos no semiplano esquerdo são mantidos.)
Por que locus circular → banda passante maximalmente plana?
O polinômio de Butterworth |B_N(jω)|² = 1 + (ω/ω_c)^{2N} tem todos os seus pólos em |s| = ω_c. A restrição de raio igual significa que todos os pólos contribuem igualmente para a resposta de magnitude em ω = ω_c. O teorema de máxima planicidade: entre todos os polinômios de N-ésima ordem com pólos neste círculo, o polinômio de Butterworth tem o máximo de derivadas iguais a zero em ω = 0.
Pólos de Chebyshev: Locus Elíptico
Os pólos de Chebyshev ficam em uma elipse no plano s (não um círculo). A elipse tem eixos semi-maiores e semi-menores determinados pelo parâmetro de ondulação ε. A banda passante de ondulação igual emerge da propriedade de equioscilação dos polinômios de Chebyshev.
Pólos Elípticos: Locus de Função Elíptica
Os pólos de filtro elíptico (Cauer) também ficam em uma elipse — mas com AMBOS os pólos E zeros contribuindo para a resposta de frequência. Os zeros ficam no eixo imaginário (pólos de atenuação finita na faixa de parada). O mapeamento de função elíptica distribui otimamente os zeros para alcançar ondulação igual em ambas as faixas simultaneamente.
Computando Localizações de Pólos de Butterworth
Para um filtro de Butterworth de 4ª ordem com ω_c = 1 (normalizado), os pólos estão em:
s_k = e^{iπ(2k+3)/8} para k = 0, 1, 2, 3
k=0: s₀ = e^{i3π/8} (no semiplano esquerdo)
k=1: s₁ = e^{i5π/8} (no semiplano esquerdo)
k=2: s₂ = e^{i7π/8} (no semiplano esquerdo)
k=3: s₃ = e^{i9π/8} (no semiplano esquerdo)
Estes quatro pólos ficam em espaçamento angular igual no círculo unitário, todos com partes reais negativas (semiplano esquerdo).
Distância de Pólos até o Círculo Unitário
A estabilidade teórica requer |p| < 1. Na prática, duas preocupações adicionais surgem.
Margem de Estabilidade
A margem de estabilidade de um filtro IIR é a distância mínima de qualquer pólo até o círculo unitário: min_k (1 − |p_k|).
Um pólo em |p| = 0.99 é tecnicamente estável mas deixa apenas 1% de margem. A aritmética de precisão finita (arredondamento na representação de coeficientes & acumulação de erros de arredondamento) pode efetivamente mover pólos. Se a quantização de coeficientes desloca um pólo de 0.99 para 1.001, o filtro se torna instável.
Consequência Geométrica
Pólos muito próximos ao círculo unitário produzem picos de resposta de frequência muito agudos — ressonadores de largura de banda estreita. Mas ressonadores estreitos requerem alta precisão: pequenos erros de coeficiente movem a frequência de pico significativamente.
O compromisso geométrico: agudeza de pico ∝ 1 / (1 − |p|). Conforme |p| → 1, agudeza → ∞ mas margem de estabilidade → 0 e sensibilidade aos erros de coeficiente → ∞.
Seções de Segunda Ordem
Um filtro IIR de alta ordem implementado como um polinômio único é numericamente sensível — arredondar um único coeficiente pode mover muitos pólos. A solução padrão: implementar como uma cascata de seções de segunda ordem (biquads), cada uma com apenas um par de pólos conjugados e um par de zeros conjugados. Erros em uma seção não podem perturbar pólos em outras.