python - 如何在python中创建带阻滤波器

Question

我有一个 .wav 音频文件，我想做的是过滤某些频率。噪声为 1000Hz，我想创建一个带阻滤波器来过滤 1000Hz 的噪声。我对python很陌生，所以我不知道我是否做得对。我想我做错了，但我不知道使用 signal.butter() 函数的正确方法是什么。

from scipy import signal    
b,a = signal.butter(4, [999,1001], 'bandstop', sampling_rate, output='ba')  
w, h = signal.freqz(b , a)  
plt.plot(w, 20*np.log10(abs(h)))

它生成的图绝对不是带阻滤波器。

标签： pythonsignalssignal-processing

Answer 1

您正在寻找的是一个陷波滤波器。它在 scipy.signal 库中实现：https ://docs.scipy.org/doc/scipy/reference/generated/scipy.signal.iirnotch.html

来自页面：“陷波滤波器是具有窄带宽（高品质因数）的带阻滤波器。它拒绝窄频带并使其余频谱几乎没有变化。”