生成随机信号并使用斯戈莱菲尔特. 指定3的多项式顺序和11的帧长度。绘制原始信号和平滑信号。

顺序=3；framelen=11；lx=34；x=randn（lx，1）；sgf=sgolayfilt（x，顺序，帧）；图（x，“:”)举行在地块（sgf，'.-')传奇(“信号”,“斯戈莱”)

这个斯戈莱菲尔特函数通过将信号与数据中心行进行卷积来执行大部分滤波B，输出sgolay。结果是滤波信号的稳态部分。生成并绘制此部分。

m=（framelen-1）/2；B=sgolay（order，framelen）；stable=conv（x，B（m+1，：），“相同”)；情节（稳定）图例(“信号”,“斯戈莱”,“稳定”)

靠近信号边缘的样本不能放置在对称窗口的中心，必须进行不同的处理。

为了确定启动瞬态，矩阵乘以第一个（framelen-1）/2行B首先框架信号的样本。

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

为了确定终端瞬态，矩阵乘以终端（framelen-1）/2行B通过最后的框架信号的样本。

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

将瞬态和稳态部分连接起来，以生成完整的信号。

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

在最小化中增加权重打破了对称B并且需要额外的步骤才能找到合适的解决方案。

加载采样的语音信号 ${\mathit{F}}_{\mathit{s}}=7418\mathrm{赫兹}$。该文件包含一段女性声音的录音，该声音说“MATLAB®”

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

对长度为21的数据帧应用多项式阶数为9的Savitzky-Golay滤波器平滑信号。绘制原始信号和滤波信号。放大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]）保持从