Resposta de Freqüência do Filtro de Corrida Média A resposta de freqüência de um sistema LTI é o DTFT da resposta ao impulso. A resposta ao impulso de uma média de movimento da amostra L é: Como o filtro de média móvel é FIR, a resposta de freqüência reduz à soma finita. pode usar a identidade muito útil para escrever a resposta de freqüência como onde deixamos ae menos jomega. N 0 e M L menos 1. Podemos estar interessados na magnitude dessa função para determinar quais freqüências passam pelo filtro não atenuadas e quais são atenuadas. Abaixo está um gráfico da magnitude desta função para L 4 (vermelho), 8 (verde) e 16 (azul). O eixo horizontal varia de zero a pi radianos por amostra. Observe que nos três casos, a resposta de frequência tem uma característica de baixa passagem. Um componente constante (frequência zero) na entrada passa pelo filtro não atenuado. Certas frequências mais altas, como pi / 2, são completamente eliminadas pelo filtro. No entanto, se a intenção era projetar um filtro lowpass, então não fizemos muito bem. Algumas das freqüências mais altas são atenuadas apenas por um fator de cerca de 1/10 (para a média móvel de 16 pontos) ou 1/3 (para a média móvel de quatro pontos). Nós podemos fazer muito melhor que isso. O gráfico acima foi criado pelo seguinte código Matlab: omega 0: pi / 400: pi H4 (1/4) (1-exp (-iomega4)) ./ (1-exp (-iomega)) H8 (1/8 ) (1-exp (-iomega8)) ./ (1-exp (-iomega)) H16 (1/16) (1-exp (-iomega16)) ./ (1-exp (-iomega)) plot (omega , abs (H4) abs (H8) abs (H16)) eixo (0, pi, 0, 1) Copyright copy 2000- - Universidade da Califórnia, BerkeleyThe Scientist and Engineers Guide to Digital Signal Processing Por Steven W. Smith, Ph. D. Como o nome indica, o filtro de média móvel opera calculando a média de um número de pontos do sinal de entrada para produzir cada ponto no sinal de saída. Na forma de equação, isso está escrito: Onde está o sinal de entrada, é o sinal de saída e M é o número de pontos na média. Por exemplo, em um filtro de média móvel de 5 pontos, o ponto 80 no sinal de saída é dado por: Como alternativa, o grupo de pontos do sinal de entrada pode ser escolhido simetricamente em torno do ponto de saída: . 15-1 de: j 0 a M -1, a: j - (M -1) / 2 a (M -1) / 2. Por exemplo, em um filtro de média móvel de 10 pontos, o índice, j. pode executar de 0 a 11 (média de um lado) ou -5 a 5 (média simétrica). A média simétrica requer que M seja um número ímpar. A programação é um pouco mais fácil com os pontos em apenas um lado, no entanto, isso produz um deslocamento relativo entre os sinais de entrada e saída. Você deve reconhecer que o filtro de média móvel é uma convolução usando um kernel de filtro muito simples. Por exemplo, um filtro de 5 pontos tem o kernel do filtro: 82300, 0, 1/5, 1/5, 1/5, 1/5, 1/5, 0, 08230. Ou seja, o filtro da média móvel é uma convolução do sinal de entrada com um pulso retangular com uma área de um. A Tabela 15-1 mostra um programa para implementar o filtro de média móvel. Movendo o Filtro Médio (filtro MA) Carregando. O filtro de média móvel é um simples filtro FIR de baixa passagem (resposta de impulso finito) comumente usado para suavizar uma matriz de dados / sinal amostrados. Ele toma M amostras de entrada de cada vez e pega a média dessas amostras M e produz um único ponto de saída. É uma estrutura LPF (Low Pass Filter) muito simples que é útil para cientistas e engenheiros filtrarem componentes ruidosos indesejados dos dados pretendidos. À medida que o comprimento do filtro aumenta (o parâmetro M), a suavidade da saída aumenta, enquanto as transições nítidas nos dados são tornadas cada vez mais abruptas. Isto implica que este filtro tem uma excelente resposta no domínio do tempo, mas uma fraca resposta de frequência. O filtro MA executa três funções importantes: 1) Toma M pontos de entrada, calcula a média desses pontos M e produz um único ponto de saída. 2) Devido à computação / cálculos envolvidos. o filtro introduz uma quantidade definida de atraso. 3) O filtro atua como um filtro de baixa passagem (com resposta de domínio de baixa frequência e uma boa resposta de domínio de tempo). Código Matlab: Após o código matlab, simula a resposta no domínio do tempo de um filtro de Média Móvel de ponto M e também plota a resposta de freqüência para vários comprimentos de filtro. Resposta no Domínio do Tempo: No primeiro gráfico, temos a entrada que está entrando no filtro da média móvel. A entrada é ruidosa e nosso objetivo é reduzir o ruído. A figura a seguir é a resposta de saída de um filtro de média móvel de 3 pontos. Pode ser deduzido da figura que o filtro de Média Móvel de 3 pontos não fez muito em filtrar o ruído. Aumentamos os toques de filtro para 51 pontos e podemos ver que o ruído na saída reduziu muito, o que é mostrado na próxima figura. Aumentamos as derivações para 101 e 501 e podemos observar que mesmo que o ruído seja quase zero, as transições são drasticamente atenuadas (observe a inclinação em cada lado do sinal e compare-as com a transição ideal da parede de tijolos em nossa entrada). Resposta de Freqüência: A partir da resposta de freqüência, pode-se afirmar que o roll-off é muito lento e a atenuação da faixa de parada não é boa. Dada essa atenuação de banda de parada, claramente, o filtro de média móvel não pode separar uma banda de freqüências de outra. Como sabemos, um bom desempenho no domínio do tempo resulta em baixo desempenho no domínio da frequência e vice-versa. Em resumo, a média móvel é um filtro de suavização excepcionalmente bom (a ação no domínio do tempo), mas um filtro passa-baixas excepcionalmente ruim (a ação no domínio da frequência) Links Externos: Livros Recomendados: Barra Lateral Principal
No comments:
Post a Comment