实验七 IIR数字滤波器的MATLAB实现
一、 实验目的:要求掌握IIR数字滤波器的设计原理、设计方法和设计步骤;能根据给定的滤波器指标进行滤波器设计;掌握数字巴特沃斯滤波器、数字切比雪夫滤波器的设计原理和步骤。
二、 实验仪器:电脑一台,MATLAB6.5或更高级版本软件一套。 三、 实验内容:
(一) 实验原理及实例分析
1. IIR数字滤波器的传递函数及特点
设IIR滤波器的输入序列为x(n),则IIR滤波器的输入序列x(n)与输出序列y(n)之间的关系可以用下面的方程式表示:
y(n) bix(n i) ajx(n j)
i 0
j 1
MN
其中,aj和bi是滤波器的系数,其中aj中至少有一个非零。与之相对应的差分方程为:
Y(z)b0 b1z 1 ....bMz MH(Z)
X(z)1 a1z 1 ....aNz N
由传递函数可以发现无限常单位冲激响应滤波器有如下特点: (1) 单位冲激响应h(n)是无限长的。
(2) 系统传递函数H(z)在有限z平面上有极点存在。
(3) 结构上存在着输出到输入的反馈,也就是结构上是递归型的。 2. IIR数字滤波器的设计与实现
IIR数字滤波器的设计有多种方法,如频率变换法、数字域直接设计以及计算辅助设计等。下面只介绍频率变换设计法。首先考虑由模拟低通滤波器到数字低通滤波器的转换,其基本的设计过程如下:
(1) 将数字滤波器的技术指标转换为模拟滤波器的技术指标; (2) 设计模拟滤波器G(S);
(3) 将G(S)转换成数字滤波器H(Z);
在低通滤波器的设计基础上,可以得到数字高通、带通、带阻滤波器的设计流程如下:
(1) 给定数字滤波器的设计要求(高通、带阻、带通); (2) 转换为模拟(高通、带阻、带通)滤波器的技术指标; (3) 转换为模拟低通滤波器的指标;
(4) 设计得到满足第三步要求的低通滤波器传递函数; (5) 通过频率转换得到模拟(高通、带阻、带通)滤波器; (6) 变换为数字(高通、带阻、带通)滤波器。
在MATLAB中设计IIR滤波器的方法及其它们所用到的函数如表6-1所示。