python 数字FIR带通滤波器

最新推荐文章于 2024-11-15 15:52:39 发布

2021乐乐

最新推荐文章于 2024-11-15 15:52:39 发布

阅读量2.1k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： python 文章标签：数字信号处理 FIR 巴特沃斯

本文链接：https://blog.youkuaiyun.com/weixin_42528089/article/details/99690963

python 专栏收录该内容

106 篇文章

订阅专栏

本文通过使用Python的SciPy库，展示了如何生成包含10Hz和20Hz频率的正弦波信号，并应用Butterworth带通滤波器进行信号过滤。此外，还介绍了如何利用快速傅里叶变换（FFT）进行频谱分析，可视化处理前后信号的对比。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

from scipy import signal
from scipy.fftpack import fft,ifft
import matplotlib.pyplot as plt
import seaborn
import numpy as np

t = np.linspace(0, 1, 1000, False)  # 1 second
sig = np.sin(2*np.pi*10*t) + np.sin(2*np.pi*20*t)
fig, (ax1, ax2) = plt.subplots(2, 1, sharex=True)
ax1.plot(t, sig)
ax1.set_title('10 Hz and 20 Hz sinusoids')
ax1.axis([0, 1, -2, 2])

sos = signal.butter(10, [15,25], 'bp', fs=1000, output='sos')
filtered = signal.sosfilt(sos, sig)
ax2.plot(t, filtered)
ax2.set_title('After 15 Hz high-pass filter')
ax2.axis([0, 1, -2, 2])
ax2.set_xlabel('Time [seconds]')
plt.tight_layout()
plt.show()