Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- k = audioread('eq12.wav'); %áudio
- fs = 22050;
- %descomenta os soundsc pra ouvir as coisas
- noise = randn(length(k),1)*0.2; %ruído
- n = [1:length(k)]';
- tone = sin(2*pi*0.2*n); % tom
- q = k %audio + ruido + tom;
- %soundsc(q,fs);
- Q = fft(q,1000); %transformada de fourier do sinal + ruido + tom
- figure(1);
- plot(linspace(0,2*pi,1000)/pi,(abs(Q)));
- xlabel('\omega/\pi')
- ylabel('|Y(e^{j*\omega})|')
- title('Magnitude do sinal + ruido')
- figure(2);
- plot(linspace(0,2*pi,1000)/pi,20*log10(abs(Q)));
- xlabel('\omega/\pi')
- ylabel('|Y(e^{j*\omega})|(dB)')
- title('Magnitude do sinal dB')
- %% aqui a gente começa a tratar o sinal
- hlpf = fir1(300, 0.2); % filtro passa baixa
- [H omega] = freqz(hlpf,1,66589); %resposta em frequencia do filtro
- plot(omega/pi,abs(H));
- title('Frequency response of the lowpass filter')
- xlabel('\omega/\pi')
- ylabel('|H(e^{j*\omega})|')
- filtered_signal = filter(hlpf,1,q); % filtrando o sinal
- soundsc(filtered_signal,fs)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement