Loading [MathJax]/extensions/HTML-CSS/handle-floats.js
Filtros Digitais - Filtro FIR Passa-faixa - Janelamento de Hamming, Hann, Blackman e Kaiser em Matlab - Aquela Questão é o que você estava procurando

terça-feira, 18 de setembro de 2018

Textual description of firstImageUrl

Filtros Digitais - Filtro FIR Passa-faixa - Janelamento de Hamming, Hann, Blackman e Kaiser em Matlab

Filtro Digital FIR Passa-faixa com Implementação em MatLab 

Questão: Usando a função fir1 do Matlab, projete um filtro FIR passa-faixa de fase linear com as seguintes especificações: frequências limites da faixa de rejeição de 0,6\pi e 0,8\pi, da faixa de passagem 0,65\pi e 0,75\pi, atenuação máxima na faixa de passagem de 0,2 dB e atenuação mínima na faixa de rejeição de 42 dB. Utilize as janelas de Hamming, Hann, Blackman e Kaiser, e mostre a resposta em magnitude de cada um dos filtros obtidos. Comente os resultados.

Resposta: Conteúdo da resposta.

Código em MatLab para filtro FIR passa-faixa com janelamento Hamming, Hann, Blackman e Kaiser

  1. %% ====== Codigo - Filtro FIR (Janelas) ====== %%
  2. % %
  3. % Aluno: Cleber Gobira %
  4. % Disciplina: Processamento Digital de Sinais %
  5. % %
  6. %% ====== Codigo - Filtro FIR (Janelas) ====== %%
  7.  
  8. FxaR = [0.6 0.8]; %Faixa de Rejeicao
  9. FxaP = [0.65 0.75]; %Faixa de Passagem em radianos
  10. wc = FxaP; %Banda passante
  11.  
  12. fsamp = 1000; %frequencia de amostragem
  13. %Em (Hz). Multiplicado por 1000 por coveniencia
  14. fcuts = [0.3 0.325 0.375 0.4]*1000
  15. mags = [0 1 0]; %magnitude
  16. d1 = 10^(-0.05*42) %Ondulacao da banda de corte
  17. d2 = (10^(0.05*0.2)-1)/(10^(0.05*0.2)+1) %ondulacao da banda de passagem
  18. devs = [d1 d2 d1]; %Definicao do desvio
  19.  
  20. [n,Wn,beta,ftype] = kaiserord(fcuts,mags,devs,fsamp); %para obter a ordem
  21. JanKaiser = fir1(n,Wn,ftype,kaiser(n+1,beta));
  22.  
  23. %% Janela de Hanning
  24. h1 = hann(n+1);
  25. JanHann = fir1(n,Wn,h1);
  26.  
  27. %% janela HAMMING
  28. h2 = hamming(n+1);
  29. JanHamm = fir1(n,Wn,h2);
  30.  
  31. %% janela blackman
  32. b1 = blackman(n+1);
  33. JanBlac = fir1(n,Wn,b1);
  34.  
  35. %% Plotando
  36. plotall = fvtool(JanHamm,1,JanHann,1,JanBlac,1,JanKaiser,1);
  37. legend(plotall,'Hamming','Hanning','Blackman','Kaiser');
  38. title('Resposta em Magnitude (Janela Hamming, Hanning, Blackman, kaiser)');

Clique na imagem para melhor visualização
Filtro passa-faixa FIR com janelamento Kaiser com fase linear.
Figura 1 – Filtro passa-faixa FIR com janelamento Kaiser com fase linear.

Fase linear do filtro passa-faixa FIR com janelamento Kaiser.
Figura 2 – Fase linear do filtro passa-faixa FIR com janelamento Kaiser.

Filtro passa-faixa FIR com janelamento Hamming, Hanning, Blackman e Kaiser, com fase linear.
Figura 3 – Filtro passa-faixa FIR com janelamento Hamming, Hanning, Blackman e Kaiser, com fase linear.

Para reduzir o efeito do fenômeno de Gibbs, deve-se usar janelas com truncamento menos abrupto. Em alguns casos é necessário a utilização de janelas mais sofisticadas, o uso de algumas delas podem ser observadas nas figuras 1 e 3, que são algumas das mais usuais. São poucas as diferenças entre as janelas Hamming, Hanning e Blackman, e essas diferenças estão voltadas a redução do fenômeno de Gibbs.

É importante salientar que entre as janelas citadas a janela de Kaiser é a única onde é possível especificar os parâmetros do filtro. Além disso elá é considerada ótima porque provê um lóbulo principal largo para a dada atenuação da banda de corte, o que implica a mais brusca banda de transição. Na figura 3, pode-se comparar a implementação com diferentes tipos de janela e constatar que quando usado a janela de kaiser é possível conseguir com muito mais precisão o requisitos de projeto do filtro FIR.

REFERÊNCIAS:
OPPENHEIM, Alan V.; SCHAFER, Ronald W. Discrete-time signal processing. Pearson Education, 3rd Edition, 2016.

Tags: filtro digital matlab filtros digitais gibbs janela janelas hamming hanning blackman fir kaiser fase linear janelamento passa-faixa passa faixa alta baixa fir1 

Nenhum comentário:

Postar um comentário

Postagens em Destaque