第4章___数字滤波器设计.doc

第四章 数字滤波 第二节 数字滤波的设计 通信与电子信息当中，在对信号作分析与处理时，常会用到有用信号叠加无用 噪声的问题。这些噪声信号有的是与信号同时产生的，有的是在传输过程中混入的， 在接收的信号中，必须消除或减弱噪声干扰，这是信号处理中十分重要的问题。根 据有用信号与噪声的不同特性，消除或减弱噪声，提取有用信号的过程就称为滤波。 滤波器的种类很多，实现方法也多种多样，本章利用Matlab来进行数字滤波器的 设计。 数字滤波器是一离散时间系统，它对输入序列x(n)进行加工处理后，输出序列 y(n)，并使y(n)的频谱与x(n)的频谱相比发生某种变化。 由DSP理论得知，无限长冲激响应（IIR）需要递归模型来实现，有限长冲激 响应（FIR）滤波器可以采用递归的方式也可采用非递归的方式实现。本章把FIR 与IIR滤波器分别用Matlab进行分析与设计。数字滤波器的结构参看《数字信号处 理》一书。 数字滤波器的设计一般经过三个步骤： 1．给出所需滤波器的技术指标。 2．设计一个H（Z），使其逼近所需要的技术指标。 3．实现所设计的H（Z）。 4.2.1 IIR数字滤波器设计 设计IIR数字滤波器的任务就是寻求一个因果、物理可实现的系统函数H(z)， 使它的频响H(ejω)满足所希望得到的低通频域指标，即通带衰减Ap、阻带衰减Ar、 通带截频ωp、阻带截频ωr。而其它形式的滤波器由低通的变化得到。 采用间接法设计IIR数字滤波器就是按给定的指标，先设计一个模拟滤波器， 进而通过模拟域与数字域的变换，求得物理可实现的数字滤波器。从模拟滤波器变 换到数字滤波器常用的有：脉冲响应不变法和双线性变换法。IIR滤波器的设计 过程如下 数字频域指标模拟频域指标设计模拟滤波器H(S) 设计数字滤波器H(z) 1. 模拟滤波器简介 模拟滤波器的设计方法已经发展得十分成熟，常用的高性能模拟低通滤波器有 巴特沃斯型、切比雪夫型和椭圆型，而高通、带通、带阻滤波器则可以通过对低通 进行频率变换来求得。必须指出，这三种滤波器都是非线性的相频特性。具体设计 方法参见《数字信号处理》一书。 （1）巴特沃斯低通滤波器的特点是：通、阻带均为单调下降。这种单调下降的特 性使得系统的误差分布不均匀。在设计中，如果在通带满足指标，阻带指标就过于 好了。或者说，在阻带满足指标，通带指标也过于好了。这样不利于以最小的阶数 来满足设计指标。 （2）切比雪夫模拟低通滤波器 分成Ι型和ΙΙ型：Ι型是通带等波纹，阻带单调 下降； ΙΙ型滤波器是通带单调下降，阻带等波纹。切比雪夫的等波纹特性使得 可以用较小的阶数设计出符合指标的滤波器。 （3）椭圆滤波器特点是通、阻带均为等波纹。椭圆滤波器的误差均匀分布在通、 阻带，比较上面几种滤波器，在同样衰耗指标下，椭圆滤波器所要求的阶数最小。 在同样阶数情况下，椭圆滤波器的通带到阻带的变化最陡峭、性能最好。 2．模拟滤波器映射到数字滤波器 从模拟滤波器变换到数字滤波器常用的有：冲激响应不变法和双线性变换法。 （1）冲激响应不变法 冲激响应不变法的基本思路为：设模拟滤波器的单位冲激响应为ha（t），数 字滤波器的单位脉冲响应h(n)，是对ha（t）的均匀取样，即 H(S)h(t) h(n)= h(nT) H(z) 冲激响应不变法由于保持了模拟滤波器的冲激响应在取样时刻的取样值，所以 具有较好的时间特性。但由于从模拟域（s域）变换到数字域（z域）的映射关系 不是一一对应的关系，因而频谱混叠难以避免, 所以在实际应用中有它的局限性。 3． Matlab滤波器设计函数简介 使用MATLAB信号处理工具箱中提供的函数可以很容易地实现IIR滤波器设计。 常用的IIR数字滤波器设计函数有： （1）IIR滤波器阶数（order）选择 buttord - 巴特沃斯（Butterworth）滤波器阶数选择 cheb1ord - 切比雪夫（Chebyshev）I型滤波器阶数选择 cheb2ord - 切比雪夫（Chebyshev）II型滤波器阶数选择 ellipord - 椭圆（Elliptic）滤波器阶数选择 buttord用来确定巴特沃斯滤波器的阶数： [N,Wc] = buttord (Wp,Wr,Ap,Ar,options) Ap,Ar为通带最大衰减和阻带最小衰减，以dB为单位； Wp,Wr为归一化通带截频和阻带截频。 Options=’low’, ’high’, ’bandpass’, ’stop’, ’s’, 分别对应低通、高通、带通、带阻、 模拟滤波器，默认情况下，为低通或带通。对于带通和带阻滤波器Wp=[Wp1,