使用 Scipy firwin 在 Python fir1 matlab 中编写

Question

我需要将信号处理 Matlab 包翻译成 Python。我被困在 fir1 函数中。 Matlab中的Fir1描述： b = fir1(n,Wn) 使用汉明窗设计具有线性相位的 n 阶低通、带通或多频带 FIR 滤波器。过滤器类型取决于Wn的元素个数。

我试过了：

import scipy
b = scipy.signal.firwin(11,.5, window = "hamming", pass_zero = False)
-0.00506032, -0.  0.04194288, -0. -0.28848483, 0.49679547, -0.28848483,
 -0. 0.04194288, -0.  -0.00506032

在 Matlab 或 Octave 中：

   fir1(10,.5)
0.00505978  -0.00016286  -0.04194054   0.00066187   0.28847599   0.49581153   0.28847599    
-0.00016286   0.00505978

Answer 1

来自 Matlab 的fir documentation，

如果 Wn 是标量，则 fir1 设计一个 低通或高通滤波器，其截止频率为 Wn...

因此，似乎提供pass_zero = False 可能与Matlab 的结果不符。相反，请尝试pass_zero='lowpass' 进行正确比较。