信號處理之filter去噪

1、filter概念及函數說明

filter函數是一維的數字濾波器，主要的應用語法以下所示

y=filter（b,a,x）;其中b爲分子係數向量，a爲分母系數向量，x表示輸入序列；

下面採用M點的滑動平均系統去噪實驗說明一下：

N = 201;
n = 0.8*rand(1,N)-0.5;//隨機噪聲
k=0:N-1;
s=2*k.*(0.9.^k)+2.0*cos(0.02*pi*k)+5.0;//信號
x = s + n;//加噪
subplot(2,1,1);
plot(k, n, 'k-', k, s, 'b--', k, x, 'm-');
xlabel('Time index k');
legend('n[k]', 's[k]', ' x[k]');
M = 5;
b = ones(M,1)/M;
a = [1]; 
y = filter(b, a, x);
subplot(2,1,2);
plot(k, s, 'b-', k, y, 'r-');
xlabel('Time index k');
legend('s[k]', 'y[k]');

實驗結果圖

可見對於一些簡單、基礎的信號，用filter濾波器去噪效果較好；可是它容易淹沒再強噪聲中，在此不作詳細說明了。