我使用这个过滤器的Python:带通滤波器不尊重截止
def bandpass_firwin(ntaps, lowcut, highcut, fs, window='hamming'):
nyq = 0.5 * fs
taps = firwin(ntaps, [lowcut, highcut], nyq=nyq, pass_zero=False,
window=window, scale=False)
在我的ntaps = 128;低切= 0.7; highcut = 4; fs = 61
我过滤我的信号,其中有61个样本以61 Hz采样(因此它的长度为10秒)。
当我看这已经被过滤通过该带通滤波器的信号的频谱,我看到:
在此光谱中的PEEK是在0.61赫兹。这不在0.7到4赫兹的范围内。
这怎么可能? &我该如何预防?
我无法重现您的结果。我使用你的代码得到了一个非常好的带通滤波器。你怎么生成这个阴谋? – 2013-05-03 09:21:08
哦,这是你的过滤数据。 – 2013-05-03 09:23:47
是的,它已经被卷入。 taps_hamming = bandpass_firwin(ntaps,0.7%,4,FS = FS) Ynew3 = np.convolve(Ynew2,taps_hamming, “相同”) – Ojtwist 2013-05-03 09:25:55