scipy.signal.butter实现带通滤波器

最新推荐文章于 2025-04-12 02:14:37 发布

我有明珠一颗

最新推荐文章于 2025-04-12 02:14:37 发布

阅读量2.1k

点赞数 3

文章标签： scipy numpy python 带通滤波

本文链接：https://blog.csdn.net/m0_37738114/article/details/130372152

版权

该代码定义了一个函数`butter_bandpass_filter`，它利用`scipy.signal`库中的`butter`和`lfilter`函数实现0.5-50Hz的带通滤波。输入数据是一个NumPy数组，表示EEG信号，经过滤波处理后输出过滤后的信号。

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

from scipy.signal import butter, lfilter


# 带通滤波器
def butter_bandpass_filter(data, lowcut, highcut, fs, order):
    low = lowcut * 2 / fs
    high = highcut * 2 / fs
    b, a = butter(order, [low, high], btype='bandpass')
    y = lfilter(b, a, data)
    return y

# 初始信号 data: NumPy Tensor (Sample, EEG Channels, Time Length)
data = np.ndarray([200, 20, 3840])  
data = butter_bandpass_filter(data, 0.5, 50, 128, 4)  # 0.5-50Hz 带通滤波