c - Matlab decimate(vector, order, 'fir') 内部工作原理

Question

有谁知道使用什么特定函数来确定滤波器系数decimate(vector, order, 'fir')？具体是fir1与Wn = 0.5? 它说订单是30，所以这是肯定的。我得到了这些系数，但我无法将结果复制到decimate.

C = fir1(30, 0.5)

-1.6994e-03 8.8095e-05 2.9361e-03 ... 2.9361e-03 8.8095e-05 -1.6994e-03

我正在尝试将 MatLab 代码翻译为 C 并遇到缺乏文档的问题。

Answer 1

使用 Octave 的源代码确认过滤器是fir1(30, 1/order).

Answer 2

fir1 用作截止频率，0:1 之间的归一化数字。由于抽取频率的影响是将信号频谱“扩展”了 r 倍，因此 fir1 所做的是确保扩展后信号仍位于归一化频率 0:1 之间。我还尝试了使用 cconv 将滤波器应用于我的信号的代码。与抽取函数比较后，显示出相位延迟。然后我尝试了 filtfilt 并且相位延迟消失了。

lx = length(x);
filt = fir1(30,1/r);
w = filtfilt(filt,1,x);
y = zeros(1,lx/r);
for i=1:lx/r
    y(i) = w(r*i);
end