Filtering out high frequency noise of an audio file

Question

Arquelau 2016-12-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/316754-filtering-out-high-frequency-noise-of-an-audio-file

回答： Star Strider 2016-12-13

Hi,

I have an audio that consists of a human voice with a high frequency noise in the background and I want to filter out the noise. Below is the code, in which I use an ideal low pass filter:

rect = @(x) 0.5*(sign(x+0.5) - sign(x-0.5));          % function used in the ideal filter
[audioIN,fs] = audioread('agudo.m4a');                                                
[numSamples, channels] = size(audioIN);                                       
audioIN = audioIN(:,1);  
Ts = 1/fs;                                                     
t = 0:Ts:(numSamples-1)*Ts;                       
fscale = -fs/2:1/max(t):fs/2;
plot(t,audioIN);title('input audio (time domain)');
fftM = fft(audioIN);
fftM = fftshift(fftM);
fftM = abs(fftM);
figure;
plot(fscale,fftM);title('input audio (freq domain)');
filter = rect(fscale/(2*3000));
figure;
plot(fscale,filter);title('filters transfer function');
fftM = fftM.*filter';
figure;
plot(fscale,fftM);title('output audio (freq domain)');
audioOUT = ifft(fftM);
audioOUT = ifftshift(audioOUT);
audioOUT = abs(audioOUT);
audioOUT = audioOUT./max(abs(audioOUT));
figure;
plot(t,audioOUT);title('output audio (time domain)');
filename = 'filtrado.m4a';
audiowrite(filename, audioOUT, fs);

However, in the output, I am not getting what I was expecting (human voice with no noise): I can only hear some random noise (not the high frequency noise of the input audio file).

Can anybody help me? Many thanks.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2016-12-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/316754-filtering-out-high-frequency-noise-of-an-audio-file#answer_247082

Human speech is essentially band-limited to be about 100-6000 Hz. Probably the easiest way to filter your signal is to use bandpass filter with * 50 Hz and 6100 Hz stopband. Use the Signal Processing Toolbox designfilt function to design it. Use the freqz function to be certain it does what you want it to, and the filtfilt function to do the actual filtering.