matlab一次差分滤波器,差分数字滤波器及其存在原理 - IIR滤波器零相位数字滤波实现及应用...

3 差分数字滤波及其存在的问题

IIR数字滤波器的设计方法主要有双线性Z变换法和冲激响应不变法；FIR数字滤波器的设计方法主要有窗函数法、频率抽样法和切比雪夫逼近法等。有关数字滤波器设计的详细方法。

对于设计好的滤波器，确定出滤波器的系数，通过对差分方程式(1)或式(2)的叠代算法，可以实现IIR滤波器或FIR滤波器的数字滤波。考虑到程序的通用性，这里就按式(1)编写一个差分滤波函数DiffFilter(b，a，x)，调用参数为三个数组，b用于存放滤波器的系数b(r)(r=0，1，…，M)；a用于存放滤波器的系数a(k)(k=0，1，…，N)，且a(0)=1;x用于存放输入信号，返回后用于存放滤波后的输出信号。如果将调用参数a均赋零值，则可实现FIR滤波器的数字滤波。

在此，用一实例进行说明。图2为某型坦克发动机转速1000rpm时，在排气口测得的噪声信号(采样点数1024，采样频率4KHz)。从信号的波形可以看出，它具有一定的周期性，主要是以发动机周期性排气产生的噪声为主。

图2某型坦克发动机噪声信号

采用窗函数(Hanning窗)法设计出某200阶FIR低通滤波器，截止频率为200Hz，分析频率为2KHz，其幅频曲线及相频曲线如图3所示。