Savitzky-Golay滤波器设计GydF4y2Ba
设计了多项式阶Savitzky-Golay FIR平滑滤波器GydF4y2BaBGydF4y2Ba
=斯戈莱(GydF4y2Ba订单GydF4y2Ba
,GydF4y2Ba框架GydF4y2Ba
)GydF4y2Ba订单GydF4y2Ba
和帧长度GydF4y2Ba框架GydF4y2Ba
.GydF4y2Ba
指定权重向量,GydF4y2BaBGydF4y2Ba
=斯戈莱(GydF4y2Ba订单GydF4y2Ba
,GydF4y2Ba框架GydF4y2Ba
,GydF4y2Ba权重GydF4y2Ba
)GydF4y2Ba权重GydF4y2Ba
,其中包含最小二乘法最小化过程中使用的实数正值权重。GydF4y2Ba
[GydF4y2Ba
返回矩阵GydF4y2BaBGydF4y2Ba
,GydF4y2BaGGydF4y2Ba
) = sgolay (GydF4y2Ba___GydF4y2Ba)GydF4y2BaGGydF4y2Ba
分化的过滤器。可以将这些输出参数与前面的任何输入语法一起使用。GydF4y2Ba
Savitzky-Golay滤波器用于平滑具有较大频率范围的噪声信号。Savitzky-Golay平滑滤波器往往比标准平均FIR滤波器滤除较少的信号高频内容。但是,当噪声水平特别高时,它们在抑制噪声方面不太成功。GydF4y2Ba
一般来说,滤波包括将信号的每个点替换为以该点为中心的移动窗口中包含的信号值的某些组合,前提是附近的点测量的基本值几乎相同。例如,移动平均滤波器将每个数据点替换为周围数据点的局部平均值nts。如果给定的数据点GydF4y2BaKGydF4y2Ba指向左边和GydF4y2BaKGydF4y2Ba指向右侧,窗口总长度为GydF4y2BaLGydF4y2Ba= 2GydF4y2BaKGydF4y2Ba+ 1GydF4y2Ba,则移动平均滤波器进行替换GydF4y2Ba
Savitzky-Golay滤波器通过最小二乘拟合一个GydF4y2BaNGydF4y2Ba通过窗口中的信号值,并将拟合多项式曲线的计算中心点作为新的平滑数据点。对于给定点,GydF4y2BaxGydF4y2BasGydF4y2Ba,GydF4y2Ba
或者,用矩阵来表示,GydF4y2Ba
为了求Savitzky-Golay估计,使用伪逆GydF4y2BaHGydF4y2Ba计算GydF4y2BaA.GydF4y2Ba然后通过GydF4y2BaHGydF4y2Ba:GydF4y2Ba
为了避免病态,GydF4y2BasgolayGydF4y2Ba
使用GydF4y2BaqrGydF4y2Ba
用于计算经济规模分解的函数GydF4y2BaHGydF4y2Ba像GydF4y2BaHGydF4y2Ba=GydF4y2BaQGydF4y2BaRGydF4y2Ba,在这方面GydF4y2BaBGydF4y2Ba=GydF4y2BaQGydF4y2BaQGydF4y2BaTGydF4y2Ba.GydF4y2Ba
有必要进行计算GydF4y2BaBGydF4y2Ba只有一次。大多数信号点的Savitzky-Golay估计是由信号与的中心行卷积得到的GydF4y2BaBGydF4y2Ba. 结果是滤波信号的稳态部分。第一GydF4y2BaKGydF4y2Ba行GydF4y2BaBGydF4y2Ba产生初始瞬态和最终态GydF4y2BaKGydF4y2Ba行GydF4y2BaBGydF4y2Ba产生最终瞬态。请参阅GydF4y2Ba斯戈莱菲尔特GydF4y2Ba
例如,可以通过增加窗口长度来改善噪声抑制,但这会在任何瞬态附近引入类似于吉布斯现象的振铃。GydF4y2Ba
索福克勒斯·奥法尼迪斯GydF4y2Ba信号处理概论GydF4y2Ba恩格尔伍德悬崖,新泽西州:普伦蒂斯大厅,1996年。GydF4y2Ba
[2] 出版社,威廉·H.,索尔·A·特科尔斯基,威廉·T·维特林和布赖恩·P·弗兰纳里。GydF4y2BaC语言中的数值配方:科学计算的艺术GydF4y2Ba纽约:剑桥大学出版社,1992年。GydF4y2Ba
什么GydF4y2Ba是GydF4y2BaSavitzky-Golay过滤器吗?(讲稿)。”GydF4y2BaIEEE信号处理杂志GydF4y2Ba第28卷第4期,2011年7月,第111-117页。GydF4y2Bahttps://doi.org/10.1109/MSP.2011.941097GydF4y2Ba.GydF4y2Ba
过滤器GydF4y2Ba
|GydF4y2Bafir1GydF4y2Ba
|GydF4y2Ba冷杉GydF4y2Ba
|GydF4y2Ba斯戈莱菲尔特GydF4y2Ba