有谁知道在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并且缺乏文档。
答案 0 :(得分:1)
使用Octave的源代码确认过滤器为fir1(30, 1/order)
。
答案 1 :(得分:0)
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