生成随机信号并使用它使用Sgolayfilt.。指定3的多项式顺序，帧长为11.绘制原始和平滑的信号。

订单= 3;Framelen = 11;lx = 34;x = randn（lx，1）;sgf = sgolayfilt（x，订单，framelen）;plot（x，“:”)举行在情节（SGF，'.-'） 传奇（“信号”那'sgolay'）

这Sgolayfilt.函数通过将信号与中心行卷积来执行大多数过滤B.，输出斯戈拉。结果是滤波信号的稳态部分。生成并绘制这部分。

m =（Framelen-1）/ 2;b = sgolay（顺序，ramelen）;稳定= 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{S.}}=7418\mathrm{赫兹}$。该文件包含一个录制的女性语音，说明“matlab®”这个词。

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

通过将多项式订单9的Savitzky-Golay滤波器应用于长度21的数据帧来平滑信号。绘制原始和过滤的信号。放大0.02-第二间隔。

RD = 9;fl = 21;SMTLB = SGOLAYFILT（MTLB，RD，FL）;子图（2,1,1）曲线（T，MTLB）轴（[0.2 0.22 -3 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]）保持从