在進行數字訊號處理的時候經常需要對時間序列進行濾波。工具箱中常用的函式是filter

常見形式： y = filter(b,a,x)
其中x是輸入序列，向量b&a是濾波器的引數。y是輸出序列。
下面舉一個例子說明一下實際的計算過程。
比如要對一段語音訊號做預加重處理，其實就是要讓響應的數字序列通過一個高通濾波器。濾波器的傳遞函式模型為——H(z)=1-μ*z^(-1)。實際計算的時候可以這麼處理：y(i) = x(i) - μ*x(i-1)，當index小於1的時候按照 x(index) = 0 處理。
程式碼如下：

[y fs]=audioread('anger_03a01Wa.wav'
 
);
yy = filter([1 -0.9],1,y);
yyy = zeros(length(y),1);
yyy(1) = y(1);
for k = 2:length(y)
    yyy(k) = y(k) - 0.9*y(k-1);
end
%yy = 5*yy;
figure(1)
subplot(3,1,1);plot((1:length(y)),y);title('原始');xlabel('序號');ylabel('幅值'); grid;
subplot(3,1,2);plot((1:length(yy)),yy);title('呼叫濾波器函式');xlabel('序號');ylabel('幅值'
 
); grid;
subplot(3,1,3);plot((1:length(yyy)),yyy);title('按照原理做序列運算');xlabel('序號');ylabel('幅值'); grid;
di = yyy - yy;%兩個結果的差別
figure(2)
plot((1:length(di)),di);title('兩個結果的計算差別');xlabel('序號');ylabel('幅值'); grid;

結果如下圖所示：
原始訊號與兩種計算方法所得到的鋁箔之後的訊號：

兩種計算方法的差別（實際上是沒有差別的）：