How to create a band pass filter

Question

Anastasiya Kobrina 2015-5-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/214906-how-to-create-a-band-pass-filter

回答： Jeremy 2015-6-18

Hi guys,

I would like to pass my .wav file through a band pass filter. I am not sure what order to use. To give an example- my .wav file contains USVs- aka my band is actually from 20 khz to 100 khz (beyond human hearing).

Sampling frequency is 50 kHz

 fs = 50e3;
 t = linspace(0,1,50e3);
 % 1 kHz and 3 kHz sine waves
 x = cos(2*pi*1e3*t)+0.5*sin(2*pi*3e3*t)+randn(size(t));
 % Lowpass filter everything below 20 kHz
 % Specify the filter
 hlpf = fdesign.lowpass('Fp,Fst,Ap,Ast',20e3,20.1e3,0.5,50,50e3);
 % Design the filter
 D = design(hlpf);
 % apply the filter
 y = filter(D,x);
 subplot(211)
 plot(t(1:1000),x(1:1000)); title('Original Waveform');
 subplot(212)
 plot(t(1:1000),y(388:1000+387)); title('Filtered Waveform');
 figure;
 subplot(211)
 plot(psd(spectrum.periodogram,x,'Fs',fs,'NFFT',length(x)));
 title('Original Signal PSD');
 subplot(212);
 plot(psd(spectrum.periodogram,y,'Fs',fs,'NFFT',length(x)));
 title('Filtered Signal PSD');

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

请先登录，再进行评论。

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

Answer 1

Jeremy 2015-6-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/214906-how-to-create-a-band-pass-filter#answer_183234

you discuss a bandpass filter but you are applying a lowpass filter...

Aside from that, the order really depends on what type of data you have and what you are trying to accomplish, I know this doesn't really help. But you basically want to minimize the filter order that will meet your requirements. The higher the order, the higher the ripple and phase shift, so start low.

You are using a filter definition that defines the requirements and the filter order will be a result of those requirements. With Fp and Fst so close this will be a very high order filter.

You are also not specifying the design method, equiripple will be the default.