ランダムな信号を生成して,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”,“稳定”）

信号エッジの近傍のサンプルは対称ウィンドウの中心に配置できないため,異なる方法で扱わなければなりません。

開始時の過渡状態を求めるため,行列はBの最初の(framelen-1) / 2行を信号の最初のframelenサンプルで乗算します。

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

終了時の過渡状態を求めるため,行列はBの最後の(framelen-1) / 2行を信号の最後の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);子图(2,1,1)plot(t,mtlb) axis([0.2 0.22 -3 2]) title(“原始”图(t,smtlb)轴([0.2 0.22 -3 2])标题(“过滤”网格)

計算を繰り返しますが,今回はカイザーウィンドウを重みベクトルとして使用します。形状係数 $\beta ＝38$を指定します。新しいフィルター処理された信号をプロットします。

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