产生一个随机信号，并平滑它使用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{年代}}＝7418\mathrm{赫兹}$．文件中有一段女声的录音说的是“MATLAB®”。

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

通过对长度为21的数据帧应用多项式9阶的Savitzky-Golay滤波器来平滑信号。绘制原始信号和滤波后的信号。放大0.02秒的间隔。

rd = 9;fl = 21;smtlb = sgolayfilt (mtlb, rd, fl);Subplot (2,1,1) plot(t,mtlb)轴([0.2 0.22 -3 2])title(“原始”)网格subplot(2,1,2) plot(t,smtlb)轴([0.2 0.22 -3 2])title(“过滤”网格)

重复计算，但是现在使用Kaiser窗口作为权重向量。指定形状因子 $\beta ＝38$．绘制新的滤波信号。

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