語法規(guī)則

y = sgolayfilt(x,order,framelen)

y = sgolayfilt(x,order,framelen,weights)

y = sgolayfilt(x,order,framelen,weights,dim)

語法描述

y = sgolayfilt(x,order,framelen)：對數(shù)據(jù)向量x使用Savitzky-Golay FIR平滑濾波器。如果x是一個矩陣，則sgolayfilt對每一列進行操作。多項式階次order必須小于框長度framelen，因此framelen必須是奇數(shù)。如果order=framelen-1，則濾波器不進行平滑。

y = sgolayfilt(x,order,framelen,weights)：指定了一個權(quán)重向量weights，它為正實數(shù)，用于最小二乘估計的最小化。如果該值沒有指定，或定義為空[]，它默認為單位矩陣。

y = sgolayfilt(x,order,framelen,weights,dim)：指定了維數(shù)dim。如果dim沒有被指定，則sgolayfilt對第一個不為1的維進行操作；指定為1則對列向量進行操作，指定2則對行向量進行操作。

舉例

穩(wěn)態(tài)和瞬態(tài)Savitzky-Golay濾波器

生成一個隨機信號并使用sgolayfilt平滑。指定多項式階數(shù)為3，框長度為11.畫出原始的和平滑后的信號。

order = 3;

framelen = 11;

lx = 34;

x = randn(lx,1);

sgf = sgolayfilt(x,order,framelen);

plot(x,':')

hold on

plot(sgf,'.-')

legend('signal','sgolay')

sgolayfilt函數(shù)通過對sgolay的輸出B的中間行進行卷積來執(zhí)行大部分濾波。結(jié)果是過濾的信號的穩(wěn)態(tài)部分。生成并畫出這個部分。

m = (framelen-1)/2;

B = sgolay(order,framelen);

steady = conv(x,B(m+1,:),'same');

plot(steady)

legend('signal','sgolay','steady')

靠近信號邊緣的樣本不能放置在對稱窗口的中心，并且必須以不同的方式進行處理。

為了確定啟動瞬態(tài)，矩陣將B的第一個(framelen-1)/2行乘以信號的第一個framelen長度的采樣。

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

為了確定終端瞬態(tài)，矩陣將B的最后的(framelen-1)/2行乘以信號的最后的framelen長度的采樣。

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

連接瞬態(tài)和穩(wěn)態(tài)部分以產(chǎn)生完整的信號。

cmplt = steady;

cmplt(1:m) = ybeg;

cmplt(lx-m+1:lx) = yend;

plot(cmplt)

legend('signal','sgolay','steady','complete')

hold off

添加權(quán)重會破壞B的對稱性，并且需要額外步驟來提供合適的解決方法。

對聲音信號的Savitzky-Golay濾波

用三階Savitzky-Golay平滑mtlb信號，數(shù)據(jù)框長度為41.

load mtlb

smtlb = sgolayfilt(mtlb,3,41);

subplot(2,1,1)

plot(1:2000, mtlb(1:2000))

axis([0 2000 -4 4])

title('mtlb')

grid

subplot(2,1,2)

plot(1:2000,smtlb(1:2000))

axis([0 2000 -4 4])

title('smtlb')

grid

總結(jié)

以上是生活随笔為你收集整理的sgolayfilt函数_Matlab中Savitzky-Golay filtering（最小二乘平滑滤波）函数sgolayfilt的使用方法...的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。