我的test.data文件有两列,时间(x)和rmsd(y)(10000行):
0.0000000 0.0005082
0.0100000 0.0685015
0.0200000 0.0784738
0.0300000 0.0837026
。
。
。
99.9700012 0.1598396
99.9800034 0.1640600
99.9900055 0.1681409
100.0000076 0.1670843
我需要绘制一条旨在产生信号平滑的线,但不使用经典的移动平均线,而是使用Savitzky-Golay算法。我想要这样做是因为,如果移动平均值的窗口大小太短,则噪声不会适当降低:
相反,如果窗口大小太大,则会发生“延迟效应”,并且未正确考虑峰:
在这种情况下,我已经检查了GNUPLOT演示页面,但是找不到任何实现。您可以检查http://www.aip.de/groups/soe/local/numres/bookfpdf/f14-8.pdf,https://c.mql5.com/forextsd/forum/147/sgfilter.pdf和http://195.134.76.37/applets/AppletSmooth/Appl_Smooth2.html以获得有关该算法的详细信息。