Tradeoff Sampling rate / bandpass filter design

Question

Peter 2015-8-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/234394-tradeoff-sampling-rate-bandpass-filter-design

评论： Star Strider 2015-8-19

Dear all, I'm usually sampling my data with 20 kHz and want to filter the data with a butterworth filter. My question: if I filter the data from e.g. 30-90 Hz. If I use the simple b,a filter logic

nyquist=10000; 
Wn=[30, 90]/nyquist;   
[b,a]=butter(4,Wn,'bandpass'); 
filtered_signal=filtfilt(b,a,signal);

Is it still a good filter design, or is the order of 4 or even 3 to small to gain a valid filter result? Was I read in the threads, there is a tradeoff between the sampling rate, the filter order and the narrowness of the passband....

Thanks for your experience. Best, Peter

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2015-8-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/234394-tradeoff-sampling-rate-bandpass-filter-design#answer_189645

I suggest using the buttord function to calculate the filter order, and then use the second-order-section implementation in the actual filter operation. The full discussion is here.

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

Star Strider 2015-8-19

在 MATLAB Online 中打开

I’m not sure you need this line:

hd=dfilt.df2sos(sos,g);

I would just use the ‘sos,g’ results in the filter.

Filter design is heuristic. I don’t know what you want to do or what is not working, since you didn’t describe it.

Experiment with different passband and stopband values, since that is the most difficult part of bandpass (or bandstop) filter design. I usually begin specifying the stopband with [0.5, 1/0.5]*Wp (with ‘Wp’ the desired passband in a bandpass filter), and then increase the rolloff to [0.8, 1/0.8]*Wp or so until I get the behaviour I want.

Alternatively, consider a different filter design, for instance Chebyshev or something with a sharper rolloff, if that’s the problem.

请先登录，再进行评论。