Questions about Butterworth filter

8 次查看(过去 30 天)
Sam
Sam 2014-12-24
I need to apply a Butterworth filter 4th-order with 6Hz cutoff frequency. This is what I get:
[B,A] = BUTTER(4,0.6)
Is this correct? And what does 'B' and 'A' mean? Thanks.

回答(1 个)

Star Strider
Star Strider 2014-12-24
Use lower case letters:
[B,A] = butter(4, 0.6);
‘B’ and ‘A’ are the coefficients of the numerator and denominator coefficients of the filter transfer function, respectively.
The cutoff frequency is normalised by the Nyquist frequency, so a normalised frequency of 0.6 means you have a Nyquist frequency of 10 Hz and a sampling frequency of 20 Hz. If those are not your Nyquist and sampling frequencies, you need to redesign your filter.
I do these steps whenever I design a filter:
Fs = ...; % Sampling Frequency
Fn = Fs/2; % Nyquist Frequency
Fc = 6/Fn; % Passband Frequency (Normalised)
Fz = 7/Fn; % Stopgand Frequency (Normalised)
[n,Wn] = buttord(Fc, Fs, 1, 10); % Order Of Filter
[B,A] = butter(n,Wn); % Designs Lowpass Filter By Default
I always use the freqz function at this stage, to be sure the filter is stable. If it is not, I either redesign it or add these additional steps:
[sos,g] = tf2sos(B,A); % Second-Order-Section Implementation
Then use the filtfilt function to filter your signals, to avoid phase distortion in the filter.
  5 个评论
Sam
Sam 2015-1-3
[n,Wn] = buttord(Fc, Fs, 1, 10);
I think 'Fs' must be replaced by 'Fz'. Am I correct?
Star Strider
Star Strider 2015-1-3
Yes. It is the normalised stopband frequency. I couldn’t test the code (I didn’t have your sampling frequency, Fs) or I’d have caught that.

请先登录,再进行评论。

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by