工程中设计一个数字滤波器

数字设计中的信号处理常需要进行滤波器设计。本文从需求设计出发，分享在工程中的滤波器如何设计。

需求：

如上图所示，只需要调制波信息，theta信息，输入信号中载波信号为sin(wt)，ADC的采样率为载波的8倍，要想保留sin(theta)就需要滤除2w±w'频段，其中theta=(w't)（公式推导省略）因此需要设计低通滤波器。

确定截止频率为w，则2w±w'频段可以被滤除。

（1）借助python求解滤波器因子：

1、使用FIR滤波器实现，确实阶数为64阶，高阶FIR滤波效果较好。

2、截止频率归一化，归一化频率为0.25.

cutoff = fc/nyq  # 截止频率归一化
numtaps = 64  # 阶数
b = signal.firwin(numtaps, cutoff)

采取如上设计，可以看到需求中的2w±w'频段被过滤。

b就是所需要的滤波器参数。对于实现而言只需要滤波器参数即可。

在数字实现侧，只需要将b（浮点型）根据所需精度转化成定点数即可，设计的滤波器增益为1，FIR参数前后对称，采用S1.15定标（S代表有符号，小数点前整数位，小数点后小数位），可以将FIR参数转化为定点数。

（2）定点化处理，通常采用python fxpmath进行定点处理。

from fxpmath import Fxp
RAW = Fxp(None, dtype='S1.15')

以上两步便可得到定点化后的滤波器因子，进行卷积操作即可。