IIR滤波器的MATLAB实现.pdfVIP

2、IIR滤波器的 MATLAB 实现 2.1 IIR 滤波器的设计方法及原理 IIR数字滤波器是一种离散时间系统，其系统函数为： M b zk k y(z) H(Z) k0  N x(z)  k 1 a z k k1 假设M≤N，当M＞N时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系 统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数a 和b ，它是数 k k 学上的一种逼近问题，即在规定意义上 （通常采用最小均方误差准则）去逼近系 统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近， 就得到数字滤波器。 2.1.1 用脉冲相应不变法设计IIR数字滤波器 利用模拟滤波器来设计数字滤波器，也就是使数字滤波器能模仿模拟滤波器 的特性，这种模仿可以从不同的角度出发。脉冲响应不变法是从滤波器的脉冲响 应出发，使数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应 ha(t)，即将ha(t)进行等间隔采样，使h(n)正好等于ha(t)的采样值，满足 h(n)h (nT) a 式中,T是采样周期。如果令Ha(s)是ha(t)的拉普拉斯变换，H(z)为h(n)的Z变换， 利用采样序列的Z变换与模拟信号的拉普拉斯变换的关系得 1  1   2  X(z)  X (sjk ) X sj k  zesT T a s T a T  (1-1) k k 则可看出，脉冲响应不变法将模拟滤波器的 S 平面变换成数字滤波器的 Z 平面， s z z sT 这个从 到 的变换 =e 是从 S 平面变换到 Z 平面的标准变换关系式。 j 3 / T jIm[z ]  / T o  － 1 o 1 Re[z ] － / T －3 / T Z平 面 S平 面 图 1-1 脉冲响应不变法的映射关系 由（1-1）式，数字滤波器的频率响应和模拟滤波器的频率响应间的关系为    1  2 k j  H(e ) H j  (1-2) T a T  k