我正在尝试估计多次应用简单的1-2-1滤波器的效果,并确定残留尺度。具体来说,我正在尝试重现这个情节:
我使用了scipy.signal.freqz,如下所示
filt = np.array([0.25,0.5,0.25])
w, h=signal.freqz(filt)
我认为对于重复滤波器,我只需要将h自身乘以多次(因为它在频域中,并且滤波只是卷积。)
但是,我不能像在报纸上那样得到同样的情节。我有三个主要问题,
我认为1-2-1过滤器只是三角形滤波器,还有其他方法可以检查频域的响应吗?
如何在python中检查重复1-2-1过滤器的频率响应?是不是它本身多次?
我很难理解freqz输出中的w(归一化频率)单位。有人可以向我解释如何转换为情节中的波数吗?
谢谢。
答案 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