Advertisement
Guest User

Untitled

a guest
Sep 30th, 2016
70
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.97 KB | None | 0 0
  1. k = audioread('eq12.wav'); %áudio
  2. fs = 22050;
  3.  
  4. %descomenta os soundsc pra ouvir as coisas
  5.  
  6. noise = randn(length(k),1)*0.2; %ruído
  7. n = [1:length(k)]';
  8. tone = sin(2*pi*0.2*n); % tom
  9.  
  10. q = k %audio + ruido + tom;
  11. %soundsc(q,fs);
  12.  
  13. Q = fft(q,1000); %transformada de fourier do sinal + ruido + tom
  14.  
  15. figure(1);
  16. plot(linspace(0,2*pi,1000)/pi,(abs(Q)));
  17. xlabel('\omega/\pi')
  18. ylabel('|Y(e^{j*\omega})|')
  19. title('Magnitude do sinal + ruido')
  20.  
  21. figure(2);
  22. plot(linspace(0,2*pi,1000)/pi,20*log10(abs(Q)));
  23. xlabel('\omega/\pi')
  24. ylabel('|Y(e^{j*\omega})|(dB)')
  25. title('Magnitude do sinal dB')
  26.  
  27. %% aqui a gente começa a tratar o sinal
  28. hlpf = fir1(300, 0.2); % filtro passa baixa
  29. [H omega] = freqz(hlpf,1,66589); %resposta em frequencia do filtro
  30.  
  31. plot(omega/pi,abs(H));
  32. title('Frequency response of the lowpass filter')
  33. xlabel('\omega/\pi')
  34. ylabel('|H(e^{j*\omega})|')
  35.  
  36. filtered_signal = filter(hlpf,1,q); % filtrando o sinal
  37.  
  38. soundsc(filtered_signal,fs)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement