信号处理-为什么信号没有以我的截止频率完全滤除？

Question

我需要过滤信号。我想将频率保持在0到51Hz之间。以下是我正在使用的代码，部分代码来自这两个问题（Python: Designing a time-series filter after Fourier analysis，Creating lowpass filter in SciPy - understanding methods and units）：

def butter_lowpass(cutoff, fs, order=5):
    nyq = 0.5 * fs
    normal_cutoff = cutoff / nyq
    b, a = butter(order, normal_cutoff, btype='low', analog=False)
    return b, a

def butter_lowpass_filter(data, cutoff, fs, order=5):
    b, a = butter_lowpass(cutoff, fs, order=order)
    y = lfilter(b, a, data)
    return y

    # y is the original signal
    # getting the unbiased signal
    y = list(np.array(y)-sts.mean(y))
    # saving the original signal
    y_before = y

    # x is the time vector

    # for the spectrum of the original signal
    yps_before = np.abs(np.fft.fft(y_before))**2
    yfreqs_before = np.fft.fftfreq(6000, d = 0.001)
    y_idx_before = np.argsort(yfreqs_before)

    #Filtering
    order = 8
    fs = 1000.0       
    cutoff = 50.0  
    y = butter_lowpass_filter(y, cutoff, fs, order)

    # for the spectrum of the filtered signal
    yps = np.abs(np.fft.fft(y))**2
    yfreqs = np.fft.fftfreq(6000, d = 0.001)
    y_idx = np.argsort(yfreqs)


    fig = plt.figure(figsize=(14,10))
    fig.suptitle(file_name, fontsize=20)
    plt.plot(yfreqs_before[y_idx_before], yps_before[y_idx_before], 'k-', label='original spectrum',linewidth=0.5)
    plt.plot(yfreqs[y_idx], yps[y_idx], 'r-', linewidth=2, label='filtered spectrum')
    plt.xlabel('Frequency [Hz]')
    plt.yscale('log')
    plt.grid()
    plt.legend()
    plt.show()

此代码的结果是经过滤波的信号，但是，这是频谱比较：

如您所见，在100Hz之后，频谱看起来不错，但是，在50Hz和大约100Hz之间仍然存在一个分量。因此，我尝试使用高阶（20）的滤波器，但是作为输出，我得到了一个非常奇怪的光谱：

因此，我知道过滤器不可能而且永远不会是完美的，但是对我来说，这似乎有点过分。以我的经验，我总是能够以截止频率获得一个很好的滤波信号。有什么建议吗？

Answer 1

截止频率通常是传递函数的下降幅度为-6dB的地方。

增加顺序将使过滤器更陡峭，但还会根据过滤器类型添加伪像。通常，您会遇到更大的相位问题（数字问题，相位变化与阶数成正比，波动...）。

在这里，我不知道，直到截止为止，它似乎很好地遵循了原始曲线。

话说回来，20阶滤波器也非常陡峭，并且由于系数中的数字太小而将出现数值问题。通常的技巧是级联2阶滤波器并使全局滤波器遵循相同的曲线（您可以查看Linkwitz-Riley滤波器）。请注意，这些过滤器都是LTI，因此您不能即时修改参数（不是LTI）。