生成一个随机信号并使其平滑使用sgolayfilt．指定多项式阶为3，帧长为11。绘制原始的和平滑的信号。

订单= 3;framelen = 11;lx = 34;x = randn (lx, 1);山东= sgolayfilt (x,秩序,framelen);情节(x,'：'） 抓住在情节(山东,“。”)传说('信号'，“sgolay”）

的sgolayfilt函数通过将信号与的中心行进行卷积来完成大部分滤波B，输出斯戈拉．其结果就是滤波信号的稳态部分。生成并绘制此部分。

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;plot（cmplt）图例（'信号'，“sgolay”，“稳定”，“完成”） 抓住离开

向最小化添加权重打破对称性B需要额外的步骤才能找到合适的解决方案。

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

负载MTLB.t =（0：长度（mtlb）-1）/ fs;

对长度为21的数据帧应用多项式9阶萨维茨基-戈莱滤波器平滑信号。绘制原始和过滤后的信号。放大0.02秒的间隔。

rd = 9;fl = 21;smtlb = sgolayfilt (mtlb, rd, fl);子图(2,1,1)plot(t,mtlb) axis([0.2 0.22 -3 2]) title(“原始”）网格子图（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])保持离开