sgolayfilt

生成一个随机信号和光滑的使用sgolayfilt。指定一个多项式的顺序3和一个框架的长度11。绘制原始和平滑信号。

订单= 3;framelen = 11;lx = 34;x = randn (lx, 1);山东= sgolayfilt (x,秩序,framelen);情节(x,“:”)举行在情节(山东,“。”)传说(“信号”,“sgolay”)

的sgolayfilt函数执行的大部分由卷积滤波信号的中心行B的输出sgolay。结果是稳态过滤信号的一部分。这部分生成和阴谋。

m = (framelen-1) / 2;B = sgolay(顺序,framelen);稳定= conv (x、B (m + 1,:)“相同”);情节(稳定)传说(“信号”,“sgolay”,“稳定”)

样品接近信号边缘不能放置的中心对称的窗口,必须区别对待。

确定启动瞬态,矩阵乘法(framelen-1) / 2行B由第一个framelen样品的信号。

ybeg = B (1: m:) * x (1: framelen);

确定终端瞬态,矩阵乘法决赛(framelen-1) / 2行B通过最后的framelen样品的信号。

yend = B (framelen-m + 1: framelen:) * x (lx-framelen + 1: lx);

连接的瞬态和稳态部分生成完整的信号。

cmplt =稳定;cmplt (1: m) = ybeg;cmplt (lx-m + 1: lx) = yend;情节(cmplt)传说(“信号”,“sgolay”,“稳定”,“完成”)举行从

增加的重量最小化的对称性B一个合适的解决方案,需要额外的步骤。

加载一个语音信号在采样 ${\mathit{F}}_{\mathit{年代}}=\mathrm{7418年}\mathrm{赫兹}$。文件包含一个记录一个女声说“MATLAB®。”这个词

负载mtlbt =(0:长度(mtlb) 1) / Fs;

光滑多项式秩序的信号通过应用Savitzky-Golay滤波器9数据帧的长度21。绘制原始和过滤信号。放大一个0.02秒的时间间隔。

rd = 9;fl = 21;smtlb = sgolayfilt (mtlb, rd, fl);次要情节(2,1,1)情节(t, mtlb)轴([0.2 - 0.22 3 2])标题(“原始”)网格次要情节(2,1,2)情节(t, smtlb)轴([0.2 - 0.22 3 2])标题(“过滤”网格)

重复计算,但现在使用一个Kaiser窗作为权重向量。指定一个形状系数 $\beta =38$。绘制新过滤信号。

kmtlb = sgolayfilt (mtlb, rd, fl, kaiser (fl, 38岁));次要情节(2,1,2)在情节(t, kmtlb)轴([0.2 - 0.22 3 2])从