中心思想 在基于最小二乘原理对多项式进行拟合时，引入观测值的权值，有明显偏差的信号给予较小的权值，不含明显偏差的信号给予较大的权值，此时最小而成的优化函数有min(VTV) 变为min(VTPV)，其中P是根据某次平差的改正数V基于的权值，在下面的代码中P=1/|V|，迭代求解，知道两次结果的差小于给定的阈值 import numpy as npclass Robust(object):def

概念 照抄百度百科，相关推导wiki可见https://en.wikipedia.org/wiki/Savitzky%E2%80%93Golay_filter英语wiki，部分用户可能不可见，原因大家都懂 Savitzky-Golay滤波拟合法是根据NDVI时间序列曲线的平均趋势，确定合适的滤波参数，用多项式实现滑动窗内的最小二乘拟合；利用Savitzky-Golay滤波方法（基于最小二乘的卷积