audio wave filter out bandpass 900-1kHz

Question

ting po chun 2023-6-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1987169-audio-wave-filter-out-bandpass-900-1khz

评论： ting po chun 2023-6-23

Hi;

I want to filter the recorded audio file out of the 900-1kHz range and convert it into an impulse response. How can I do it?

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

请先登录，再进行评论。

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

Answer 1

Deep 2023-6-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1987169-audio-wave-filter-out-bandpass-900-1khz#answer_1261284

在 MATLAB Online 中打开

You can read your audio file using the audioread function:

[y, Fs] = audioread('audiofile.wav');

More info: audioread

To filter out the 900-1kHz range, you can use a bandstop filter. MATLAB provides the designfilt function for this purpose:

d = designfilt('bandstopiir', 'FilterOrder', 2, ...
               'HalfPowerFrequency1', 900, 'HalfPowerFrequency2', 1000, ...
               'DesignMethod', 'butter', 'SampleRate', Fs);
filtered_audio = filtfilt(d, y);

More info:

Finally, you can convert the filtered audio into an impulse response using the inverse Fast Fourier Transform, provided by the ifft function:

impulse_response = ifft(filtered_audio);

More info: ifft

This impulse response might not be in a playable range, so you might want to normalize it:

impulse_response = impulse_response / max(abs(impulse_response));

Please explore the documentation links provided to understand more about the functions used.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Deep 2023-6-23

在 MATLAB Online 中打开

Sure, plotting can be done by simply using the plot function.

% Define the sample rate. You can check this in your WAV file's properties.
Fs = 44100;
% Create a time vector
t = (0:length(impulse_response)-1)/Fs;
% Plot the impulse response
plot(t, impulse_response);
xlabel('Time (s)');
ylabel('Amplitude');
title('Impulse Response');