python中的savitzky-Golay过滤器，窗口大小错误

Question

我有一个时间序列数据，我想要平滑使用Savitzgy Golay过滤器。根据resaerch的说法：https://arxiv.org/ftp/arxiv/papers/1808/1808.10489.pdf

窗口大小应为n+2，其中n个->多面体。我对我的数据进行了这样的平滑处理，但我没有看到信号的结果。我选错窗户了吗？

我的数据是带有(50858, 2)形状的t和emg列的EMG信号，其中emg是要去噪的值。这是它的头。

​

​

过滤器实现：

Y= data.iloc[:,1].values

Y_filtered= savgol_filter(Y, window_length = 5, polyorder = 3)并绘制它：

plt.subplot(1, 2, 1)
plt.plot(Y[-1000:])
plt.title("EMG with noise")


plt.subplot(1, 2, 2)
plt.plot(Y_filtered[-1000:])
plt.title("SG filter applied ")
plt.tight_layout()
plt.show()

​

​

Answer 1

窗口大小不限于n+2。不过，一定很奇怪。我试过用window_size=21和polyorder=3，它很有效。我认为不太可能在数千个样本信号中看到窗口大小=5的结果，因为您的噪声可能会跨越更多的样本。试着用一个更大的窗口大小。