[论文解读] Optimum window length of Savitzky-Golay filters with arbitrary order

Mohammad Sadeghi, Fereidoon Behnia|arXiv (Cornell University)|Aug 30, 2018

Image and Signal Denoising Methods参考文献 17被引用 23

一句话总结

本文提出了一种优化算法，用于确定任意阶数的Savitzky-Golay滤波器的最优窗长，通过平衡偏差与方差来最小化均方误差（MSE）。通过使用切比雪夫正交多项式建模滤波器，该方法在多种信号类型下均表现出比现有方法更高的精度和鲁棒性。

ABSTRACT

One of the widely used denoising methods in different domains is the Savitzky-Golay (SG) filter. The SG filter has two design parameters: window length and the filter order. As the length of the window increases, the estimation variance decreases, but the bias error increases at the same time. Mean square error (MSE) measure includes both bias and variance criteria. In this paper, we obtain the optimal window length of an SG filter with arbitrary order which minimizes the MSE. To achieve the optimal window length, we propose an algorithm whose performance is better than the existing methods. In this paper, we follow the viewpoint proposed by Persson and Strang and design the filter on the basis of Chebyshev orthogonal polynomials

研究动机与目标

• 确定任意阶数Savitzky-Golay滤波器的最优窗长，以最小化均方误差（MSE）。
• 解决SG滤波器设计中估计方差与偏差误差之间的权衡问题。
• 开发一种比现有方法更准确、更鲁棒的算法，用于确定最优窗长。
• 通过使用切比雪夫正交多项式重新表述SG滤波器设计，以提高数值稳定性与性能。

提出的方法

• 作者采用Persson和Strang提出的方法框架，将Savitzky-Golay滤波器重新表述为基于切比雪夫权函数的正交多项式。
• 通过将信号投影到对称窗内的正交多项式基上，推导出滤波器系数的解析表达式。
• 通过将均方误差（MSE）表示为信号曲率与噪声功率的函数，确定最优窗长，以最小化MSE。
• 该算法计算不同窗长下的MSE，并选择误差最低的窗长，利用正交多项式展开的性质。
• 通过避免标准单项式基方法中常见的病态多项式基集，确保了数值稳定性。
• 在合成信号和真实世界信号上对所提算法进行了评估，以验证其在最小化MSE方面的优越性。

实验结果

研究问题

• RQ1 对于任意阶数的Savitzky-Golay滤波器，最小化均方误差（MSE）的最优窗长是什么？
• RQ2 与标准多项式基相比，使用切比雪夫正交多项式如何提高SG滤波器设计的精度与稳定性？
• RQ3 所提算法能否在不同信号特征下，优于现有方法以确定最优窗长？
• RQ4 窗长、滤波器阶数与SG滤波中产生的偏差-方差权衡之间存在何种关系？

主要发现

• 所提算法在各种信号类型下，相较于现有方法，实现了更低的均方误差（MSE）。
• 使用切比雪夫正交多项式显著提高了数值稳定性，并降低了条件数，相比标准单项式基多项式拟合。
• 研究表明，最优窗长取决于滤波器阶数和信号曲率，高阶滤波器需要更长的窗长以保持低偏差。
• 该方法提供了MSE关于窗长和滤波器阶数的闭式表达式，实现了精确优化。
• 该算法在具有不同噪声水平和曲率的合成信号上表现出稳健性能，优于传统的启发式窗长选择规则。
• 结果证实，通过该方法最小化MSE可实现比固定或经验选择窗长更优的去噪性能。