如何使用Savitzky Golay滤波器在离散采样的1D信号中找到局部最大值（在采样之间）？

我有一个地震信号y（i）：


这里我找到了一个最大值：i = 152.54，y = 222.29并用红色绘制。

我想自动找到所有的最大值。 SGF的最小值和最大值比其他过滤器要好得多。这听起来很适合我使用。

我找到了一个生成SGF系数的Matlab脚本。
并用它找到了导数的4阶SGF系数。
我编写了一个小的Matlab脚本，该脚本< />通过将信号与导数的
四阶SGF系数进行卷积来找到信号的导数，从而找到样本对（i，i + 1），其中导数更改符号
通过线性插值在i和i + 1之间找到导数的零交叉

脚本：脚本我必须测试导数是否从负变为正，才能使函数获得所需的结果，但这使我感到困惑。

下面是使用此函数在我的信号上找到最大值的结果：


结果看起来不错，但我注意到有些找不到axima：
i = 143.13、190.88、256.97。

这是因为它们接近其他最大值吗？

如何控制最接近的两个最大值？

感谢您的任何回答！