Need help to determine stopband

Question

Daemian 2015-3-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/181541-need-help-to-determine-stopband

评论： Star Strider 2015-3-4

采纳的回答： Star Strider

I got a frequency of 6 - 16 Hz to pass through. but how do i determine what is my stopband?

xd = double(myRecording); % Convert to ‘double’

Fs = 8000; % Sampling Frequency

Fn = Fs/2; % Nyquist Frequency

Fpb = [6 16]/Fn; % Passband

Fsb = []/Fn; % Stopband

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2015-3-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/181541-need-help-to-determine-stopband#answer_170083

If you’re designing a bandpass filter, everything within [6 16] Hz (normalised by ‘Fn’ as you correctly compute) is your passband. Everything else is your stopband by definition, so you don’t need to define it specifically. If you need to define passband and stopband ripple (a requirment of buttord and related functions), choose 1 and 10 dB respectively.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Star Strider 2015-3-4

在 MATLAB Online 中打开

The 1 and 10 dB are the passband and stopband ripple, respectively. This is more important in the Chebyshev design than for Butterworth filters, but buttord wants them. Choose whatever values meet your needs. Note that the butter function simply wants the filter order and (for a bandpass filter) passband frequencies. It doesn’t need any other information unless you want to design a highpass or bandstop filter, then it needs a third argument specifying filter type.

I didn’t realise you’re using the buttord function. It does require you to specify the passband and stopband, as well as the ripple for both, so I usually specify the stopbands to be 0.8 and 1.25 times the passband frequencies. I always convert the filter to second-order-section representation for stability.

For your filter, this becomes:

Fs = 8000; % Sampling Frequency
Fn = Fs/2; % Nyquist Frequency
Fpb = [6 16]/Fn; % Passband
Fsb = Fpb.*[0.8  1/0.8]; % Stopband
Rp =  1;
Rs = 10;
[n,Wn] = buttord(Fpb,Fsb,Rp,Rs);
[b,a] = butter(n,Wn);
[sos,g] = tf2sos(b,a);
figure(1)
freqz(sos,1024,Fs)

The freqz call lets you see what your filter looks like. Experiment with the stopband characteristics to get the result you want. You may want to widen them a bit.

Daemian 2015-3-4

Sure thanks a lot :)

Star Strider 2015-3-4

My pleasure!

请先登录，再进行评论。

Need help to determine stopband

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

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Need help to determine stopband

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

采纳的回答

4 个评论 显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

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

4 个评论
显示 2更早的评论隐藏 2更早的评论