def bandpass2(w, freqlo, freqhi, fs, npass=3):
wn = [2*freqlo/fs, 2*freqhi/fs]
(b,a) = signal.butter(npass, wn, 'bandpass')
w = signal.filtfilt(b, a, w)
return w
微信
公众号
扫码关注公众号