python - 为什么它在使用 scipy.signal.filtfilt 时返回 nan

Question

我有一个采样率为 160GHz 的日期数组，我希望对数据应用 30kHz-100kHz 带通滤波器。我写如下：

import scipy.signal as dsp
import numpy as np
fs=160e9
data=np.random.rand(int(1e-4*fs))
b, a = dsp.butter(5, [30e3/(fs/2),100e3/(fs/2)],btype='band')
result=dsp.filtfilt(b,a,data)

但是，结果是nan。但是，如果我像这样降低采样率，它会起作用：

import scipy.signal as dsp
import numpy as np
fs=160e9
data=np.random.rand(int(1e-4*fs))
data=data[::10000] # 100 does not work
fs=fs/10000
b, a = dsp.butter(5, [30e3/(fs/2),100e3/(fs/2)],btype='band')
result=dsp.filtfilt(b,a,data)

这样我可以得到正确的结果，但是分析起来很粗糙。因此，当滤波器频率远低于采样率时，巴特沃斯滤波器似乎无法正常工作。我怎么解决这个问题？