Projeto de Filtros Fir

PROJETO DE FILTROS FIR
1. Introduo
FILTRAGEM: alterao das amplitudes das componentes frequenciais de um sinal.
A funo fir1 do Matlab usada no projeto e anlise de filtros digitais de Resposta

Impulsiva Finita (FIR), usando a tcnica de janelamento.
Ambas funes implementam o mtodo clssico de projeto de filtros digitais do tipo

FIR com fase linear baseado na tcnica do janelamento [1]. Ele projeta filtros nas
configuraes padres do tipo passa-baixa, passa-alta, passa-faixa e rejeita-faixa. O
filtro projetado normalizado, ou seja, a Resposta em Frequncia tem magnitude
mxima igual a 0 dB na frequncia central da banda passante.
A funo fir1 apresenta diversas formas de uso (chamada):
1) b = fir1(N,Wn)
2) b = fir1(N,Wn,'tipofiltro')
3) b = fir1(N,Wn,janela)
4) b = fir1(N,Wn,'tipofiltro',janela)
5) b = fir1(...,'normalizao')
Sintaxe 1: b = fir1(n,Wn)
Retorna o vetor b contendo os N+1 coeficientes de um filtro passa-baixa de ordem N.

Trata-se de um filtro de fase linear, baseado na janela de Hamming, com frequncia
normalizada de corte Wn. Os coeficientes do filtro, b, esto ordenados em potncias
decrescentes de z:
b(z) = b(1) + b(2).z-1 + ... + b(n+1).z-n
Wn um nmero entre 0 e 1, onde 1 corresponde frequncia de Nyquist (metade da

frequncia de amostragem, fs/2).
Se Wn um vetor com dois elementos, Wn = [w1 w2], fir1 retorna um filtro
passafaixa com banda passante w1 < w < w2.
Se Wn um vetor multi-elementos, Wn = [w1 w2 w3 w4 w5... wn], fir1 retorna um

filtro multibanda de ordem N com faixas 0 < w < w1, w1 < w < w2, ..., wn < w < 1.
Sintaxe 2: b = fir1(N,Wn,'tipofiltro')
Nesta forma de uso deve-se especificar o tipo de filtro a ser projetado:
'high' para filtro passa alta com frequncia de corte Wn.

'stop' para filtro rejeita faixa, com Wn = [w1 w2] para a faixa de rejeio.
'DC-0' para filtro multibanda, com primeira banda sendo do tipo rejeio.
'DC-1' para filtro multibanda, com primeira banda sendo do tipo passante.
A funo fir1 sempre projeta filtro de ordem par para configuraes passa alta e
rejeita faixa, pois em filtros de ordem mpar, a Resposta em Frequncia na frequncia
de Nyquist nula, o que no apropriado aos filtros passa alta e rejeita faixa.
Sintaxe 3: b = fir1(N,Wn,janela)
Usa a janela especificada no vetor coluna janela para projetar o filtro. Esse vetor deve
ter N+1 elementos. Se nenhuma janela especificada, fir1 usa a janela de Hamming.
Sintaxe 4: b = fir1(n,Wn,'tipofiltro',janela)
Idem s sintaxes 2 e 3 simultaneamente.
Sintaxe 5: b = fir1(...,'normalizao')
Indica se a magnitude do filtro deve ser normalizada, admitindo os valores: 'scale'

(default) se a normalizao for desejada (a magnitude da resposta em frequncia do
filtro na frequncia central da banda passante 0 dB); 'noscale' se a normalizao no
for desejada.
O atraso de grupo do filtro FIR projetado por fir1 n/2 amostras.
Exemplo 1: O arquivo chirp.mat do Matlab contm um sinal y[n] com a gravao do

gorjeio de um pssaro. Esse sinal tem a maior parte de sua potncia acima da frequncia
fs/4 (metade da frequncia de Nyquist). Projete um filtro FIR passa alta de ordem 34
para atenuar as componentes do sinal abaixo de fs/4. Use frequncia de corte 0,48 rad,
e janela de Chebyshev, com atenuao mnima de 30 dB na banda de rejeio.
load chirp % som monofnico, -1 <= y <= +1, fs = 8192 Hz ou amostra/s

amostras = length(y);
w = 0:1/(amostras/2):1*(1-1/(amostras/2)); % freq.Dig.: 0-pi (rad/amostra)
% Filtro Passa Alta, com janela de Chebyshev, ordem 34, wc = 0.48*pi rad/am
% fc = 0.48*fs/2 = 1966 Hz
b = fir1(34,0.48,'high',chebwin(35,30)); % coef.s do filtro FIR desejado
freqz(b,1,512); % resposta em frequncia do filtro digital
f = filter(b,1,y); % filtragem do sinal original y(n)
sound(y) % som original
pause(3)
sound(f) % som filtrado
Y = fft(y); % espectro do sinal original
F = fft(f); % espectro do sinal filtrado
figure; o = length(w);
subplot(2,1,1), plot(w,abs(Y(1:o))); title('Espectro do Sinal Original');
xlabel('Freq. Normalizada (\times\pi rad/amostra)'); ylabel('Magnitude');
subplot(2,1,2), plot(w,abs(F(1:o))); title('Espectro do Sinal Filtrado');
xlabel('Freq. Normalizada (\times\pi rad/amostra)'); ylabel('Magnitude');
Outras funes janelas disponveis no Matlab esto disponveis via funo:
window(@WNAME,N): retorna uma janela de N pontos do tipo especificado pelo

manipulador de funo @WNAME em um vetor coluna. @WNAME pode ser
qualquer nome de funo janela reconhecido pelo Matlab.
@bartlett - Janela de Bartlett
@barthannwin - Janela de Bartlett-Hanning
@blackman - Janela de Blackman
@blackmanharris - Janela Mnima de Blackman-Harris com 4 termos
@bohmanwin - Janela de Bohman
@chebwin - Janela de Chebyshev
@flattopwin - Janela Flat Top
@gausswin - Janela Gaussiana
@hamming - Janela de Hamming
@hann - Janela de Hann
@kaiser - Janela de Kaiser
@nuttallwin - Janela de Nuttall
@parzenwin - Janela de Parzen (de la Valle-Poussin)
@rectwin - Janela Retangular (boxcar)
@taylorwin - Janela de Taylor
@tukeywin - Janela de Tukey
@triang - Janela Triangular
2. Procedimentos Experimentais
1. Para o mesmo sinal usado no Exemplo 1, realize nova filtragem passa alta,
modificando a frequncia de corte do filtro projetado 3.fs/8.
2. Compare os espectros de frequncia dos sinais filtrados no Exemplo 1 e no item
anterior, e identifique as mudanas ocorridas no sinal de udio (ouvindo ambos
os sinais).
3. Projete um filtro FIR passa faixa de ordem 48, com banda passante de 0,65
0,75 (rad/amostra). Repita o item anterior.
3. Questes
1. Descreva a funo fir2 (Matlab).

2. O que se pode dizer do volume de som do udio resultante dos itens do
procedimento 2 e 3? O que pode ser feito para se manter o volume de som do
udio original em ambos os casos?

Projeto de Filtros Fir

Enviado por

Direitos autorais:

Formatos disponíveis

Projeto de Filtros Fir

Enviado por

Dados do documento

Descrição original:

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

Projeto de Filtros Fir

Enviado por

Direitos autorais:

Formatos disponíveis

PROJETO DE FILTROS FIR

FILTRAGEM: alterao das amplitudes das componentes frequenciais de um sinal.

A funo fir1 do Matlab usada no projeto e anlise de filtros digitais de Resposta

Ambas funes implementam o mtodo clssico de projeto de filtros digitais do tipo

A funo fir1 apresenta diversas formas de uso (chamada):

Retorna o vetor b contendo os N+1 coeficientes de um filtro passa-baixa de ordem N.

b(z) = b(1) + b(2).z-1 + ... + b(n+1).z-n

Wn um nmero entre 0 e 1, onde 1 corresponde frequncia de Nyquist (metade da

Se Wn um vetor multi-elementos, Wn = [w1 w2 w3 w4 w5... wn], fir1 retorna um

Nesta forma de uso deve-se especificar o tipo de filtro a ser projetado:

'high' para filtro passa alta com frequncia de corte Wn.

Idem s sintaxes 2 e 3 simultaneamente.

Indica se a magnitude do filtro deve ser normalizada, admitindo os valores: 'scale'

O atraso de grupo do filtro FIR projetado por fir1 n/2 amostras.

Exemplo 1: O arquivo chirp.mat do Matlab contm um sinal y[n] com a gravao do

load chirp % som monofnico, -1 <= y <= +1, fs = 8192 Hz ou amostra/s

Outras funes janelas disponveis no Matlab esto disponveis via funo:

window(@WNAME,N): retorna uma janela de N pontos do tipo especificado pelo

1. Descreva a funo fir2 (Matlab).

Você também pode gostar