反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断累积，有时会产生微弱的寄生振荡。
（3）借助成熟的模拟滤波器的成果：IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不高。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，再通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。
（4）需加相位校准网络：IIR数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准网络。
3.3 FIR数字滤波器概述
FIR数字滤波器（finite impulse response filter）又名“有限脉冲响应数字滤波器”，这类滤波器对于脉冲输入信号的响应最终趋向于0，因此而得名。有限脉冲响应滤波器（FIR filter）的优点：
（1）脉冲响应（impulse response）为有限长：造成当输入数位讯号为有限长的时候，输出数位讯号也为有限长。
（2）比无限脉冲响应滤波器（IIR filter）较容易最佳化（optimize）。
（3）线性相位（linear phase）：造成h(n),是偶对称（even）或奇对称（odd）且有限长。
（4）一定是稳定的（stable）：因为Z转换（Z transform）后所有的极点（pole）都在单位圆内
3.4 设计FIR数字滤波器和IIR数字滤波器比较
不论是IIR滤波器还是FIR滤波器的设计都包括三个步骤： (1) 按照实际任务的要求，确定滤波器的性能指标。
(2) 用一个因果、稳定的离散线性时不变系统的系统函数去逼近这一性能指标。根据不同的要求可以用IIR系统函数，也可以用FIR系统函数去逼近。 (3) 利用有限精度算法实现系统函数，包括结构选择、字长选择等。
但IIR滤波器和FIR滤波器的设计方法完全不同。IIR滤波器设计方法有间接法和直接法，间接法是借助于模拟滤波器的设计进行的。其设计步骤是：先设计过渡模拟滤波器得到系统函数H（s），然后将H（s）按某种方法转换成数字滤波器的系统函数H(z)。FIR滤波器采用间接法，常用的方法有窗函数法、频率采样法和切比雪夫等波纹逼近法。对于线性相位滤波器，经常采用FIR滤波器。
21
3.5 数字滤波器设计的基本思路
数字滤波器的实现有两个关键步骤：一个从数字域到模拟域间的变换，这个变换实现了数字滤波器技术指标到模拟滤波器技术指标的转换，同样也实现了模拟滤波器系统函数到数字滤波器系统函数的转换；另一个是从模拟滤波器技术指标到满足该指标的模拟滤波器的设计。
3.6 设计FIR滤波器
3.6.1 窗函数法及设计步骤
FIR滤波器的设计比较简单，就是要设计一个数字滤波器去逼近一个理想的低通滤波器。通常这个理想的低通滤波器在频域上是一个矩形窗。根据傅里叶变换我们可以知道，此函数在时域上是一个采样函数。通常此函数的表达式为：
sa（n）＝sin（n∩）/n（3-2）
但是这个采样序列是无限的，计算机是无法对它进行计算的。故我们需要对此采样函数进行截断处理。也就是加一个窗函数。就是传说中的加窗。也就是把这个时域采样序列去乘一个窗函数，就把这个无限的时域采样序列截成了有限个序列值。但是加窗后对此采样序列的频域也产生了影响：此时的频域便不在是一个理想的矩形窗，而是成了一个有过渡带，阻带有波动的低通滤波器。通常根据所加的窗函数的不同，对采样信号加窗后，在频域所得的低通滤波器的阻带衰减也不同。通常我们就是根据此阻带衰减去选择一个合适的窗函数。如矩形窗、汉宁窗、汉明窗、BLACKMAN窗、凯撒窗等。 我们介绍用窗函数法设计FIR滤波器的步骤。如下：
（1）根据对阻带衰减及过渡带的指标要求，选择窗函数类型（矩形窗、三角窗、汉宁窗、哈明窗、凯塞窗等），并估计窗口长度N。先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣的窗函数。 （2）构造希望逼近的频率响应函数。 （3）计算h(n).。 （4）加窗得到设计结果。
接下来，我们根据语音信号的特点给出有关滤波器的技术指标： ①低通滤波器的性能指标：
fp=600Hz，fc=800Hz，As=10db ,Ap=1dB ②高通滤波器的性能指标：
fp=3000Hz，fc=2700Hz，As=10dB，Ap=1dB；
22
在Matlab中，可以利用函数fir1设计FIR滤波器，利用函数butter,cheby1和ellip设计IIR滤波器，利用Matlab中的函数freqz画出各步步器的频率响应。 hn=fir1(M，wc，window)，可以指定窗函数向量window。如果缺省window参数，则fir1默认为哈明窗。其中可选的窗函数有Rectangular Barlrtt Hamming Hann Blackman窗，其相应的都有实现函数。
3.6.2源程序与仿真图像
下面我们将给出设计FIR数字滤波器的主要程序和图像。 FIR低通滤波器程序见附录I FIR低通滤波器图像：
图3-1 FIR低通滤波器
FIR高通滤波程序见附录I FIR高通滤波图像：
23
图3-2 FIR高通滤波器
3.7 设计IIR滤波器
3.7.1双线性变换法与设计步骤
无论是数字滤波器还是模拟滤波器，他们技术指标的建立都是以所谓的“固有衰减”参数为参照。以数字滤波器为例，固有衰减参数（关于之旅分量归一化）定义为 A(w)=-20lg|H(ejw)|-20lg|H（ej0）| （式3-1）
目前IIR数字滤波器设计的最通用的方法是借助于模拟滤波器的设计方法。模拟滤波器设计已经有了一套相当成熟的方法，它不但有完整的设计公式，而且还有较为完整的图表供查询，因此，充分利用这些已有的资源将会给数字滤波器的设计带来很大方便。IIR数字滤波器的设计步骤是：
①按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标； ②根据转换后的技术指标设计模拟低通滤波器G?s?；
24
③再按一定规则将G(s)转换成H(z)。若设计的数字滤波器是低通的，那么上述设计工作可以结束，若所设计的是高通，带通或带阻滤波器，那么还有步骤④；
④将高通、带通、或带阻数字l不去的技术指标先转化为低通模拟滤波器的技术指标，然后按照上述步骤②设计出低通G?s?，再将G?s?转换为所需的H(z)。
IIR滤波器设计的一种简单、有效的方法就是双线性变换法。可以借助于模拟滤波器的频率转换设计一个所需类型的过渡模拟滤波器，再经过双线性变换将其转换策划那个所需的数字滤波器。具体设计步骤如下： （1）确定所需类型数字滤波器的技术指标。
（2）将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率，转换公式为Ω=2/T tan(0.5ω)（式3-2）
(3)将相应类型的模拟滤波器技术指标转换成模拟低通滤波器技术指标。 （4）设计模拟低通滤波器。
（5）通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。
（6）采用双线性变换法将相应类型的过渡模拟滤波器转换成所需类型的数字滤波器。 我们知道，脉冲响应不变法的主要缺点是会产生频谱混叠现象，使数字滤波器的频响偏离模拟滤波器的频响特性。为了克服之一缺点，可以采用双线性变换法。 下面我们总结一下利用模拟滤波器设计IIR数字低通滤波器的步骤：
（1）确定数字低通滤波器的技术指标：通带边界频率、通带最大衰减，阻带截止频率、阻带最小衰减。
（2）将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。 （3）按照模拟低通滤波器的技术指标设计及过渡模拟低通滤波器。
（4）用双线性变换法，模拟滤波器系统函数转换成数字低通滤波器系统函数。
函数butter,cheby1和ellip设计IIR滤波器时都是默认的双线性变换法，所以在设计滤波器时只需要代入相应的实现函数即可。
3.7.2源程序与仿真图像
IIR低通滤波器程序见附录I； IIR低通滤波器图像：
25