如何使用Matlab制作简单的FIR滤波器？

Question

如何使用Matlab制作简单的低通FIR滤波器(不使用内置功能)？

问题示例:

Implement a FIR LPF with cut-off frequency 250Hz

也可能需要采样频率......

解决方案尝试或我已经知道的:

x = [...] -> input signal
A = 1; -> Since this is FIR
B = [?????]
y = filter(B, A, x) -> Output signal

Afaik,B应包含FIR滤波器的系数.但; 如果我只有截止频率,我该如何计算这些系数？

Answer 1

最简单的是一个"窗口sinc"过滤器:

fs = 44100;
cutoff = 250;
t = -256:256;  % This will be a 513-tap filter
r = 2*cutoff/fs;
B = sinc(r*t).*r .* blackman(length(t))';
freqz(B);

滤波器的长度(参见参考资料t=...)控制过渡带的宽度.cutoff在这种情况下是-6 dB点.blackman是一个流行窗口的名称.您可以查看此 Wikipedia页面,了解有关窗口功能的更多信息.它们在过渡带宽和阻带抑制之间基本上有不同的折衷.