Matlab抽取(矢量​​,顺序,'for')内部工作

时间:2013-07-10 16:32:45

标签: c matlab filtering translate

有谁知道在decimate(vector, order, 'fir')中使用什么特定函数来确定滤波器系数?具体来说是fir1Wn = 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并且缺乏文档。

2 个答案:

答案 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