scipy中的带通巴特沃斯滤波器频率
从
我正在按照食谱设计scipy带通滤波器。但是，如果我过多地降低了过滤频率，我最终会在高阶过滤器上产生垃圾。我究竟做错了什么？
from scipy.signal import butter, lfilter
def butter_bandpass(lowcut, highcut, fs, order=5):
nyq = 0.5 * fs
low = lowcut / nyq
high = highcut / nyq
b, a = butter(order, [low, high], btype='band')
return b, a
if __name__ == "__main__":
import numpy as np
import matplotlib.pyplot as plt
from scipy.signal import freqz
# Sample rate and desired cutoff frequencies (in Hz).
fs = 25
# Plot the frequency response for a few different orders.
plt.figure(1)
plt.clf()
for order in [1, 3, 5, 6, 9]:
b, a = butter_bandpass(0.5, 4, fs, order=order)
w, h = freqz(b, a, worN=2000)#np.logspace(-4, 3, 2000))
plt.semilogx((fs * 0.5 / np.pi) * w, abs(h), label="order = %d" % order)
plt.xlabel('Frequency (Hz)')
plt.ylabel('Gain')
plt.grid(True)
plt.legend(loc='best')
plt.figure(2)
plt.clf()
for order in [1, 3, 5, 6, 9]:
b, a = butter_bandpass(0.05, 0.4, fs, order=order)
w, h = freqz(b, a, worN=2000)#np.logspace(-4, 3, 2000))
plt.semilogx((fs * 0.5 / np.pi) * w, abs(h), label="order = %d" % order)
plt.xlabel('Frequency (Hz)')
plt.ylabel('Gain')
plt.grid(True)
plt.legend(loc='best')
plt.show()
更新：已经讨论了该问题，显然已在Scipy 0.14上解决了该问题。但是，在Scipy更新之后，该图看起来仍然非常糟糕。怎么了？
从
显然，该问题是一个已知的错误：
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
编辑于
相关文章
Related 相关文章
- 1
在MATLAB中对图像进行高通巴特沃斯滤波器
- 2
如何使用Matlab绘制巴特沃斯滤波器的频率响应
- 3
带通FIR滤波器
- 4
C/C++ 中的巴特沃斯 BP 滤波，奇怪的频谱
- 5
信号。黄油带通错误：数字滤波器的临界频率必须为0 <Wn <1
- 6
Python中的低通然后逆滤波器
- 7
低通滤波器和采样频率
- 8
Web音频滤波器降低增益而不是滤波频率
- 9
带滤波器的模式匹配
- 10
带firpm的滤波器系数
- 11
设计带阻滤波器（在 SciPy firwin）
- 12
Matlab中的傅立叶变换和LTI滤波器以及频率响应
- 13
如何在MATLAB中绘制数字滤波器的频率响应？
- 14
带通滤波器无法在Matlab中过滤不想要的频率
- 15
如何在MATLAB中绘制数字滤波器的频率响应？
- 16
使用`scipy.signal.butter（）`创建低通滤波器时，“阶次”和“临界频率”是什么
- 17
20hz-20000hz巴特沃思滤波爆炸
- 18
查找带通滤波器的频率响应
- 19
在SciPy中创建低通滤波器-了解方法和单位
- 20
如何在scipy.signal中制作一个低通滤波器？
- 21
C中的FIR滤波器？
- 22
C中的定点FIR滤波器？
- 23
python中的fft带通滤波器
- 24
C ++中的FIR滤波器
- 25
在MATLAB中实现均值滤波器
- 26
dplyr滤波器中的标量比较
- 27
Python中的带通滤波器形状
- 28
OpenCL中的低通滤波器
- 29
R中的输入滤波器
我来说两句