matlab中的数字滤波器

在进行数字信号处理的时候经常需要对时间序列进行滤波。工具箱中常用的函数是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;

结果如下图所示：
原始信号与两种计算方法所得到的铝箔之后的信号：

两种计算方法的差别（实际上是没有差别的）：