Matlab IIR 滤波器

Question

我们的任务是实现一个滤波器（数字高通 IIR 滤波器），以降低 500 Hz 以下的频率并允许更高的频率。 使用 sptool 和椭圆函数，我们构造了以下内容

% Sampling Rate (Hz)
Fs = 46875; 

% High-pass filter
N = 4; % Filter Order
Wn = 500/(46875/2);  % Cutoff Frequency in terms of Passband / Nyquist frequency ratio
Rp = 0.5; % Passband ripple specification
Rs = 20; % Stopband attenuation

[Num,Den] = ellip(N,Rp,Rs,Wn,'high');

secondOrderSection = tf2sos(Num,Den)./2

然后我们使用 L1 范数计算比例因子： 我们需要在 16 位 atmel 微控制器（AC3U3 Xplained）上实现这一点，因此必须以 Q_0.15 形式计算比例因子

% FIRST SECTION
firstScaleFactor = 1/(sum(abs(impz(1,secondOrderSection(1,4:6)))))
sf1_2_Q = round(firstScaleFactor*(2^15))

% SECOND SECTION
secondScaleFactor = 1/(sum(abs(impz(1,secondOrderSection(2,4:6)))))
sf2_2_Q = round(secondScaleFactor*(2^15))

问题：我们的比例因子似乎有点低，第一个比例因子是 77，第二个是 14。 我们的计算是否有任何错误？

Answer 1

德州仪器公司的This application note 包含在定点滤波器实现中避免溢出的技术概述。