重复1-2-1滤波器的频率响应

时间:2016-08-18 14:44:10

标签: python scipy signal-processing

我正在尝试估计多次应用简单的1-2-1滤波器的效果,并确定残留尺度。具体来说,我正在尝试重现这个情节:

from Small et al., 2013

我使用了scipy.signal.freqz,如下所示

filt = np.array([0.25,0.5,0.25])
w, h=signal.freqz(filt)

我认为对于重复滤波器,我只需要将h自身乘以多次(因为它在频域中,并且滤波只是卷积。)

但是,我不能像在报纸上那样得到同样的情节。我有三个主要问题,

  1. 我认为1-2-1过滤器只是三角形滤波器,还有其他方法可以检查频域的响应吗?

  2. 如何在python中检查重复1-2-1过滤器的频率响应?是不是它本身多次?

  3. 我很难理解freqz输出中的w(归一化频率)单位。有人可以向我解释如何转换为情节中的波数吗?

  4. 谢谢。

1 个答案:

答案 0 :(得分:0)

原来我没错。通过绘制传递函数的绝对值,并将归一化的频率除以2pi,我得到了完全相同的图,并且应用滤波器的多个时间恰好多次对自身的频率响应。

filt = np.array([0.25,0.5,0.25])
w, h=signal.freqz(filt)
plt.plot(w/(2*pi), abs(h**400), label='400 pass')

Comparison between frequency response of repeating 1-2-1 filter