Introdução aos filtros de resposta ao impulso infinito (IIR)

Os filtros de Resposta a Impulso Infinito (IIR) são uma classe fundamental de filtros digitais caracterizados por um mecanismo de feedback que permite que a saída dependa tanto de entradas atuais como anteriores, bem como de saídas anteriores. Esta estrutura recursiva permite que os filtros de IIR alcancem transições de frequência acentuadas com muito menos coeficientes do que as suas contrapartes de Resposta a Impulso Finito (FIR). Como resultado, os filtros de IIR são computacionalmente eficientes e amplamente implantados em aplicações de processamento de sinais em tempo real, como equalização de áudio, telecomunicações, instrumentação biomédica e sistemas de controle.

O design dos filtros IIR normalmente começa com um filtro de protótipo analógico – muitas vezes o tipo Butterworth, Chebyshev ou elíptico (Cauer) – que é transformado no domínio digital através de técnicas como a transformada bilinear ou invariância de impulso. Cada protótipo oferece trocas específicas entre flatness da faixa de passagem, atenuação da banda de parada, largura de banda de transição e linearidade de fase. Compreender esses trade-offs é fundamental para selecionar o filtro apropriado para uma determinada aplicação.

Conceitos Principais no Desenho de Filtros IIR

Função de Resposta e Transferência de Frequência

Um filtro IIR é descrito pela sua função de transferência no domínio z:

Os coeficientes (feedback) e (feedforward) determinam a magnitude e a resposta de fase do filtro. A ordem do filtro N (o número de pólos) influencia fortemente a inclinação do roll-off da faixa de passagem para a faixa de paragem.

Especificações do projeto chave

  • Freqüência de borda de banda (]ωp[):] a frequência superior (ou inferior) em que a resposta se encontra dentro de uma tolerância especificada.
  • Frequência de borda de banda de paragem (]ωs):]] a frequência em que a atenuação cumpre um nível mínimo exigido.
  • Ondulação da banda de passa (]Rp[):] o desvio máximo admissível do ganho ideal na faixa de passagem (em dB).
  • A ]s]]]]]]]A[][[[][]][]]a atenuação mínima exigida na parada (em dB).
  • Ordem de filtro (N): determina complexidade computacional e íngreme de transição.

Protótipo analógico e Transformação Digital

A maioria dos métodos de desenho do IIR começa com um filtro analógico de protótipos lowpass que atende às especificações de magnitude. A função de transferência analógica H(s)[ é então convertida para um plano digital H(z)[ usando uma transformação de mapeamento de frequência. O método mais comum para preservar a resposta de magnitude é a transformada bilinear, que mapeia o plano inteiro s[]-z[]-z[-plane via [s = (2/T) * (1 - z^- 1)/(1 + z^-1)[. Pré-warping as frequências críticas é necessário para compensar o mapeamento de frequências não lineares.

Filtros Butterworth: Resposta plana máxima

Características

Os filtros Butterworth são desenhados para fornecer uma resposta de frequência tão plana quanto possível na faixa de passagem, sem ondulações. A resposta de magnitude ao quadrado é dada por:

onde ω c é a frequência de corte (normalmente o ponto -3 dB) e N é a ordem do filtro. Como N aumenta, a transição torna-se mais acentuada, mas a faixa de passagem permanece monotônica.

Vantagens

  • Flat passband: Ideal para aplicações que devem preservar a amplitude do sinal em uma ampla faixa de frequência, como sistemas de áudio de alta fidelidade e monitoramento médico.
  • Resposta de fase previsível: O deslocamento de fase varia linearmente com a frequência dentro da faixa de passagem, minimizando a distorção da forma de onda para sinais transitórios.
  • Equações de projeto simples: As localizações do polo são facilmente computadas em um círculo de raio ω c no plano s, permitindo desenhos de alta ordem simples.

Limitações

  • Rolagem moderna: Alcançar uma transição acentuada requer uma alta ordem de filtro, aumentando o custo computacional e o atraso do grupo.
  • Pobre inclinação de atenuação da banda de parada: A inclinação assintótica é -20N dB/década, que é menos agressiva do que Chebyshev ou filtros elípticos para a mesma ordem.

Passos de Desenho

  1. Especificar borda de banda de passagem, borda de banda de parada, ondulação de banda de passa (tipicamente 0,5-3 dB), e atenuação mínima de banda de parada.
  2. Calcular a ordem de filtragem requerida N utilizando ].
  3. Determinar as localizações do pólo analógico: os pólos são uniformemente espaçados na metade esquerda do círculo de raio ω c.
  4. Converta para filtro digital usando transformada bilinear com pré-warping.
  5. Implementar como uma cascata de seções de segunda ordem (SOS) para estabilidade numérica.

Quadros de design detalhados e fórmulas de forma fechada para coeficientes de Butterworth estão amplamente disponíveis; veja por exemplo a referência autoritária Guia de projeto de filtro IIR da Analog Devices.

Filtros Chebyshev: Rolo Steeper com Ondulações

Tipo I e Tipo II

Os filtros Chebyshev trocam a flatness da faixa de passagem para uma transição mais íngreme. Existem duas variantes:

  • [[FLT: 0]]Chebyshev Tipo I: [FLT: 1]] As ondulações são permitidas na faixa de passagem, enquanto a faixa de parada é monotônica. A resposta de magnitude é [[FLT: 5], onde [[FLT: 2]] T N[[[FLT: 3]] é o polinomial de ordem de Chebyshev [[FLT: 4]] N[[[FLT: 5]] e ε é um fator de ondulação derivado da ondulação de faixa de passo permitida [[FLT: 6]] R p[[[FLT: 7]].
  • Chebyshev Tipo II (inverso Chebyshev): A faixa de passagem é plana, mas as ondas aparecem na faixa de paragem. Este desenho é menos comum porque os nulos da banda de paragem podem ser indesejáveis em algumas aplicações.

Vantagens

  • Roll-off mais rápido do que Butterworth:] Os filtros Chebyshev alcançam a mesma inclinação de transição com uma ordem mais baixa, reduzindo a carga computacional.
  • Bem adequado para especificações onde uma pequena onda de banda de passe é aceitável: Os valores típicos de onda de faixa de passe variam de 0,1 a 3 dB. Muitos sistemas de comunicação toleram essas ondas em troca de filtragem mais nítida.

Limitações

  • Resposta de fase não linear na banda passa: As ondas introduzem distorção de fase, que pode degradar a qualidade dos sinais que dependem da coerência de fase (por exemplo, formas de onda moduladas).
  • ]A magnitude do rípple deve ser cuidadosamente escolhida:A onda mais alta produz um roll-off mais íngreme, mas aumenta a variação do atraso do grupo e o potencial de toque no domínio do tempo.

Abordagem de projecto

  1. Define borda da banda de passagem, borda da banda de parada, ondulação da banda de passa (por exemplo, 1 dB) e atenuação da banda de parada.
  2. Calcular a ordem do filtro: .
  3. Encontre os pólos do protótipo analógico: os pólos estão sobre uma elipse no plano , com o eixo principal ao longo do eixo real para Tipo I.
  4. Mapa para digital usando transformada bilinear com frequência pré-warping.
  5. Decomponha-se em secções SOS para evitar a sensibilidade ao coeficiente.

Para um exemplo prático usando a biblioteca SciPy do Python, consulte A documentação de projeto de filtro IIR do SciPy.

Filtros elípticos: Largura de banda de transição ideal

Características

Os filtros Ellípticos (ou Cauer) conseguem o mais íngreme possível roll-off entre a faixa de passagem e a faixa de parada para uma dada ordem de filtro, permitindo ondulações em ambas as bandas. A resposta ao quadrado de magnitude envolve funções elípticas Jacobianas, tornando o design matematicamente mais envolvido.

A vantagem principal é que para uma ordem fixa N, o filtro elíptico oferece a faixa de transição mais curta. Esta propriedade torna-o o filtro de escolha em aplicações que requerem pontos de corte extremamente afiados, como filtros anti-aliasing em receptores de comunicação ou seleção de canais de banda estreita.

Vantagens

  • Roll-off mais agressivo: Alcança uma dada largura de banda de transição com a menor ordem possível em comparação com Butterworth ou Chebyshev.
  • Comercialidade flexível: Permite especificação independente de ondulação de banda de passagem, atenuação de banda de paragem e relação de transição.

Limitações

  • Rípples em ambas as bandas de passagem e parada: Pode introduzir distorção de sinal inaceitável em aplicações de medição de áudio de alta fidelidade ou precisão.
  • Design complexo: Requer avaliação de funções elípticas e algoritmos especializados. A maioria dos designers confia em software de projeto de filtro (MATLAB, Octave ou ferramentas dedicadas).
  • Sensibilidade de coeficiente maior: A transição estreita pode levar a efeitos de quantização quando implementados em processadores de ponto fixo.

Procedimento de concepção

  1. Especifique todos os quatro parâmetros: ω p, ω s, R p, A s.
  2. Calcular o factor de selectividade ]k = ω p/ω s e o factor de discriminação k1 = ε/
  3. Encontrar a ordem N usando funções integrais elípticas: , onde K[ é a integral elíptica completa do primeiro tipo.
  4. Determinar a localização do pólo e zero das funções elípticas Jacobianas; zeros estão no eixo imaginário na faixa de paragem para filtros elípticos Tipo I.
  5. Aplicar transformada bilinear para obter a função de transferência digital.
  6. Implemente como seções SOS, garantindo que zeros (que são reais para a versão digital) sejam devidamente pareados para evitar o pico de ganho.

Um tratamento abrangente do projeto de filtro elíptico pode ser encontrado na documentação da função elip da MATLAB, que fornece rotinas integradas para gerar coeficientes.

Passos práticos de design comuns a todos os tipos

1. Coleção de especificação

Comece por reunir os requisitos de aplicação: frequência de amostragem, banda de passagem e frequências de stopband (normalizadas para metade da taxa de amostragem), ondulação de banda de passagem admissível e atenuação necessária da banda de paragem. Por exemplo, um filtro anti- aliasing num sistema de áudio de 48 kHz pode exigir uma ondulação de banda de passband inferior a 0,1 dB e uma atenuação de banda de paragem de 80 dB a 24 kHz.

2. Geração de Protótipos Analógicos

Usando ferramentas ou computação manual, gerar o protótipo analógico lowpass que corresponde às especificações de magnitude nas frequências normalizadas. Pré-dobrar as frequências críticas porque a transformada bilinear comprime o eixo de frequência.

3. Conversão Digital

Aplicar a transformada bilinear para obter coeficientes digitais. Para filtros de alta ordem (N > 2), converter a função de transferência em uma cascata de seções de segunda ordem (SOS). Esta decomposição melhora muito a estabilidade numérica em implementações de pontos fixos e flutuantes.

4. Implementação e Verificação

Implemente o filtro no seu ambiente alvo (C, Python, HDL, etc.). Verifique a resposta de frequência usando uma ferramenta de análise espectral e certifique-se de que a ondulação da faixa de passagem, atenuação da banda de parada e largura de banda de transição atendam às especificações originais. Simule com sinais de teste realistas para verificar se há artefatos de domínio temporal, como toque ou sobreposição.

Comparação dos filtros Butterworth, Chebyshev e Elíptica

Property Butterworth Chebyshev (Type I) Elliptic
Passband Response Monotonic, maximally flat Ripples present Ripples present
Stopband Response Monotonic Monotonic Ripples present
Roll-off Sharpness Lowest (N-dependent) Medium Highest
Filter Order for Given Specs Highest required Lower than Butterworth Lowest required
Phase Linearity Best (near linear) Moderate Poor (significant nonlinearity)
Complexity of Design Low Moderate High

Ferramentas e software para projeto de filtro IIR

Os fluxos de trabalho de processamento de sinais digitais modernos dependem de ferramentas de alto nível que automatizam a geração de coeficientes de filtro.

  • [[FLT: 0]]Python com SciPy: [[FLT: 1]] Funções como [[FLT: 8]], [FLT: 9], [FLT: 10] e [[FLT: 11]] retornam diretamente as matrizes SOS dadas frequências normalizadas. Exemplo: [[FLT: 12]].
  • [[FLT: 0]]MATLAB / Octave: Comandos como , , , e numerador de computação e polinômios denominadores. A função ] converte-se em SOS.
  • Pacotes de design de filtro especializados: ferramentas gratuitas e comerciais, como FilterPro (Texas Instruments), Nuhertz Filter Solutions ou AADE Filter Design oferecem especificação gráfica interativa e saída de coeficiente real.

Para uma comparação exaustiva destes ambientes de design, consulte ] um inquérito de 2021 sobre ferramentas de projeto de filtro digital (ResearchGate).

Orientação de Aplicação Real-World

A escolha entre os três tipos de filtro depende das restrições específicas do sistema:

  • Processamento de áudio de alta fidelidade: Use filtros Butterworth para preservar a forma do sinal original com distorção de fase mínima. O leve comprometimento no roll-off é compensado pela qualidade de som transparente.
  • Selecção de canais de telecomunicações: Os filtros elípticos são preferidos quando a rejeição estrita do canal adjacente deve ser alcançada com recursos de hardware limitados.A distorção de fase é tolerável em sinais modulados por amplitude.
  • Acondicionamento de sinal biomédico (por exemplo, ECG):] Os filtros Chebyshev Tipo II oferecem uma faixa de passagem plana para evitar distorcer o complexo QRS, rejeitando ainda a interferência de 50/60 Hz de forma eficaz.Muitos front-ends do ECG usam uma entalhe Chebyshev Tipo II de 4a ordem.
  • Sistemas de controle: Os filtros de Butterworth são frequentemente escolhidos por sua resposta de fase previsível e falta de superação, o que é fundamental para margens de estabilidade do sistema.

Considerações Avançadas

Efeitos Finitos de Comprimento de Palavra

Ao implementar filtros de IIR em processadores de ponto fixo (por exemplo, 16 bits ou 32 bits), a quantização do coeficiente pode mudar de localização do pólo, causando instabilidade ou degradando o desempenho. Usando a implementação de SOS IIR e escalando cuidadosamente as amplitudes do sinal mitiga essas questões. Para filtros elípticos, os zeros podem estar muito próximos do círculo unitário, tornando-os especialmente sensíveis à quantização.

Designs de Passagem Alta, Passagem de Bandas e Paradas de Bandas

Todos os três tipos de protótipos podem ser transformados de lowpass para outras respostas seletivas de frequência através de transformações espectrais (substituindo o ]z[-variável com uma função racional apropriada). Por exemplo, uma transformação lowpass-to-bandpass dobra a ordem do filtro e é facilmente implementada em ferramentas de software. No entanto, as características de fase tornam-se mais complicadas após a transformação.

Comparação com os filtros FIR

Os filtros FIR oferecem estabilidade linear e garantida, mas requerem ordens muito mais elevadas para corresponder ao desempenho de roll-off IIR. Em aplicações onde a fase não é crítica, o menor custo computacional dos filtros IIR, especialmente elípticos, faz deles a escolha preferida.Para sistemas em tempo real com orçamentos de latência rigorosos, os filtros IIIR podem ser implementados com o mínimo de atraso.

Conclusão

A concepção de filtros IIR com características elípticas e Butterworth envolve balancear objetivos concorrentes de flatness da faixa de passagem, íngreme de transição e complexidade de implementação. Os filtros Butterworth fornecem o design mais simples e a resposta de fase mais previsível, tornando-os uma escolha segura para aplicações de áudio e baixa frequência de uso geral. Os filtros Chebyshev oferecem um roll-off melhorado com ondulação controlada, particularmente útil nas comunicações. Os filtros elípticos oferecem a seletividade de frequência mais agressiva para a menor ordem, embora com maior complexidade de design e ondulação de banda.

O engenheiro moderno pode aproveitar bibliotecas de software sofisticadas para rapidamente protótipo desses filtros, deixando tempo para focar em trade-offs de nível de sistema e verificação. Compreender os suportes matemáticos garante que o filtro selecionado atenda tanto as especificações de domínio de frequência quanto as restrições do mundo real do seu sistema de processamento de sinais digitais.

Para mais leituras sobre a teoria do design e derivação destes filtros, consulte O livro online de Julius O. Smith sobre filtros digitais e os textos clássicos de Oppenheim e Schafer.Além disso, o relatório de aplicação do Texas Instruments “Design of IIR Digital Filters” fornece um guia prático orientado para a implementação.