Loading

band-pass filtering 

from scipy import signal
# 带通滤波器
# 信号，频率下限，频率上限， 采样率
def bandPass(signals, fre_low, fre_high, fs):
    b, a = signal.butter(8, [2.0 * fre_low / fs, 2.0 * fre_high / fs], 'bandpass', 'sos')
    filtedData = signal.filtfilt(b, a, signals)
    return filtedData

posted @ 2020-02-03 16:30  pgCai  阅读(407)  评论(0)  编辑  收藏  举报