How to use a Fir1 filter to make a lowpass filter

Question

Anthony Koning 2022-10-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1821538-how-to-use-a-fir1-filter-to-make-a-lowpass-filter

评论： Star Strider 2022-10-10

Hi, I'm currently trying to make a lowpass filter for a sound clip with the specs of being order 128 and having a curoff frequency of 3575 Hz. Using the tf function does produce a graph, but that graph comes up as blank for any order greater than 70. People are recommending using a fir1, which, while it does produce a graph, though I'm unsure if these are correct given how I'm also getting the following error messages:

"Error using fir1>chkWindow (line 286)

Window length must be the same as the filter length.

Error in fir1>eFir1 (line 161)

wind = chkWindow(L,varargin{windIndex});

Error in fir1 (line 92)

[b,a] = eFir1(varargin{:});

Error in Lowpass (line 25)

blo = fir1(128,0.48,chebwin(35,30));"

My code is as follows:

audioread('beep.wav');

info = audioinfo('beep.wav')

filename = 'beep.wav';

audiowrite(filename,y,Fs);

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

sound(y,Fs)

t = 0:seconds(1/Fs):seconds(info.Duration);

t = t(1:end-1);

figure(1)

plot(t,y)

xlabel('Time')

ylabel('Audio Signal')

title 'Beep Signal'

filename = 'Beep.wav';

audiowrite(filename,y,Fs);

t = (0:length(y)-1)/Fs;

bhi = fir1(34,0.48,'low',chebwin(35,30));

freqz(bhi,1)

blo = fir1(128,0.48,chebwin(35,30));

outlo = filter(blo,1,y);

subplot(2,1,1)

plot(t,y)

title('Original Signal')

ys = ylim;

subplot(2,1,2)

plot(t,outlo)

title('Lowpass Filtered Signal')

xlabel('Time (s)')

ylim(ys)

Any help explaining what is going wrong with this code would be appreciated

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2022-10-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1821538-how-to-use-a-fir1-filter-to-make-a-lowpass-filter#answer_1070698

The filter specification needs to be:

blo = fir1(128,0.48,chebwin(129,30));

because of the filter length.. That will work.

Note that this is correct:

bhi = fir1(34,0.48,'low',chebwin(35,30));

They should both create essentially the same filter, although with different properties since the lengths are different.

(The ‘beep.wav’ file apparently does not exist on the server, or in my offline MATLAB installation, so I cannot run the code.)

.

2 个评论
显示无隐藏无

Anthony Koning 2022-10-10

This sis allow the graph to work, thanks for that. But just out of curiosity, shouldn't a passed audio clip filter look less like this:

and more like this:

Star Strider 2022-10-10

No, it should not.

However if you take the Fourier transform (fft) of the filtered audio clip and divide it element-wise by the Fourier transform of the original, unfiltered audio clip, and the plot the one-sided result of that operation (from 0 Hz to the Nyquist frequency), it should resemble the transfer function in the freqz plot.

The ‘Lowpass Filtered Signal’ plot is in the time domain, and the freqz plot is in the frequency domain.

请先登录，再进行评论。