我正在尝试在 SciPy 中设计一个简单的 Butterworth 带通滤波器,但遇到了一些奇怪的结果。
import scipy.signal as signal
import numpy as np
import matplotlib.pyplot as plt
def butter_bandpass(lowcut, highcut, fs, freqs,order=3, label=None):
nyq = 0.5 * fs
low = lowcut / nyq
high = highcut / nyq
sos = signal.butter(order, [low, high], btype='band', output='sos')
w, h = signal.sosfreqz(sos,worN=freqs,whole=True,fs=fs)
return w,h
freqs = 650
for i in np.arange(1,10):
w,h = butter_bandpass(0.01, 0.1, fs=1/0.68, freqs=freqs, order=i)
plt.plot(h)
这给出了奇怪的结果,如下图所示(巴特沃斯滤波器从 1-10 阶)。我认为随着订单的增加,过滤器应该变得越来越矩形?
有谁知道如何在 SciPy 中设计一个简单的巴特沃斯滤波器?