IIR数字滤波器设计数字信号处理课程设计毕业设计（论文）.doc

IIR数字滤波器设计（Ⅰ） 数字滤波器简介 在科技发展的前提下，数字滤波器在科研实践中所起到的作用越来越无法取代。数字滤波器是对数字信号处理实现滤波的线性时不变系统。数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号(数字序列)通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为是一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性： ， 其 中,分别是数字滤波器的输出序列和输入序列的频域特性(或称为频谱特性)，是数字滤波器的单位取样响应的频谱，又称为数字滤波器的频域响应。可以看出， 输入序列的频谱经过滤波后变为， 因此，只要按照输入信号频谱的特点和处理信号的目的， 适当选择 ， 使得滤波后的 满足设计的要求。 IIR滤波器简介 数字滤波器在数字信号处理的各种应用中发挥着十分重要的作用。它是通过对采样数据信号进行数学运算处理来达到滤波的目的。其中无限冲击响应数字滤波器也称IIR是采用对离散采样数据作差分方程运算来进行滤波。IIR DF的优点在于可以利用模拟滤波器设计的结果，然后采用双线性变换法或冲激响应不变法将模拟滤波器转换成数字滤波器，而模拟滤波器的设计方便简单、有大量的图表可查。IIR DF的设计还可利用二阶子系统的串联来有效实现。其主要缺点在于Nyquist间隔内无法实现线性相位,若需线性相位，比如在图象处理和数据传输中都要求信道具有线性相位特性，则需要用全通网络进行相位校正。 IIR滤波器的传递函数模型 IIR数字滤波器的基本网络结构可用差分方程、单位脉冲响应及系统函数来描述，其差分方程一般表示为 式(2.1) 式中和分别表示输入和输出信号序列，和是滤波器系数。设式(2.1)中输入信号和输出信号在以前处于零起始状态，则有 式(2.2) 对于线性时不变(LTI)系统来说，上式中和均为常数[2]。 IIR滤波器的状态方程模型 状态方程是描述系统的一种常用的方式。在系统有不可见的状态变量时可使用滤波器状态方程模型。连续LTI系统的状态方程可以写成 式(2.3) 式(2.4) 其中A，B，C，D分别为常数矩阵，在MATLAB中，一般情况下，系统的状态方程可以简记为(A，B，C，D),如果D≡0，则系统的状态方程模型可以简记为(A，B，C)。 对于离散系统来说，状态方程模型可以写成 式(2.5) 式(2.6) 这样的离散时间系统的状态方程模型在MATLAB中写成(A，B，C，D)或(A，B，C)[3]。 IIR滤波器的零极点模型 零极点模型实际是传递函数的另一种表达，它的原理是分别对原系统传递函数的分子和分母进行分解因式，以获得系统的零极点表示形式，对单输入单输出系统来说，零极点模型可以简单的表示为： 式(2.7) 式中，1,2,3…，和，1,2,3…，分别称为系统的零点和极点，它们既可以是实数也可以是复数，而K称为系统的增益。在MATLAB中简记为(Z,P,K)。但为了保证滤波器的稳定性，系统的所有极点都位于左半S平面。如果稳定系统所有的零点都位于左半S平面，则称该系统为最小相位系统，否则称为非最小相位系统。如果系统中的某个零点值恰好等于其中一个极点的值，则它们之间可以对消以直接获得一个完全等效的低价系统。 IIR数字滤波器的实现方法 冲激响应不变法 冲激响应不变法[5]的变换原理是使数字滤波器的单位冲激响应序列模仿模拟滤波器的单位冲激响应序列。将模拟滤波器的单位冲激响应加以等简隔抽样，使正好等于的抽样值，即满足 式(3.1) 其中T是抽样周期。 如果令是的拉普拉斯变换，为的Z变换，根据抽样序列的Z变换和模拟信号的拉普拉斯变换的关系，可得 式(3.2) 可以看出，利用冲激响应不变法将模拟滤波器变换数字滤波器，实际上是先对模拟滤波器的系统函数做周期延拓，再经过的映射变换，从而得到数字滤波器的系统函数。假设在S平面上，S在轴上取值，Z在Z平面内的单位圆周上取值，则可得到数字滤波器的频率响应和模拟滤波器的