医学数字信号处理章滤波器设计.doc

第六章 数字滤波器设计 滤波器在实际信号处理中起到了非常重要的作用。任何检测的信号都含有噪声，而滤波是去除噪声的基本手段。本章介绍IIR滤波器和FIR滤波器的设计。IIR滤波器设计主要内容包括：巴特沃思、切比雪夫模拟低通滤波器设计；脉冲响应不变法和双线性变换法的数字化变换方法；数字高通、带通和带阻滤波器的设计。而FIR滤波器是直接采用的数字式设计方法。针对FIR滤波器特征，首先介绍了其线性相位的实现条件，然后介绍了窗函数法和频率抽样法的设计方法。 6.1 概述 一.滤波器的分类 这里主要讨论经典滤波器的设计。 按功能划分经典滤波器又可分为低通、高通、带通、带阻四种滤波器如图6-1所示. 图6-1 理想低通、高通、带通和带阻滤波器幅度特性 经典滤波器设计从实现方法上分为IIR滤波器和FIR滤波器。它是一个线性时不变离散时间系统，如果滤波器用单位脉冲响应序列h(n)表示，其输入x(n)与输出y(n)之间的关系可以表示为：y(n)=x(n)*h(n) h(n)的Z变换称为系统函数。IIR滤波器和FIR滤波器的系统函数分别是： 二.数字滤波器的性能要求与设计步骤 这种方法在实际中都是用有限精度算法实现的线性非移变离散系统，一般的设计步骤包括： ①根据实际需要确定滤波器的技术指标，例如滤波器频率响应的幅度特性和截止频率等。 ②用一个稳定的因果系统逼近这些指标，具体说，就是有这些指标计算系统函数H（z）。 ③用有限精度的运算实现H（z），包括选择运算结构、进行误差分析和选择存储单元的字长。 第一步与实际应用有关。第三步关于滤波器的结构问题已在上节讨论了，而有关滤波器的误差分析和存储器字长的选择等问题将在后面章节中研究。 一个理想滤波器，要求所在通频带内幅频响应是一常数；相频特性相应为零或是频率的线性函数。但一个实际的滤波器是不可能得到上述幅频和相频响应。以低通滤波器为例，频率响应有通带、过渡带及阻带三个范围。如图6－2所示。 图6－2 逼近理想低通滤波器的容限 其中δ1 ：通带衰减；δ2：阻带衰减；ωp：通带截止频率;ωT：阻带截止频率;ωT－ωP：过渡带 通带最大衰减αp αp越小，通带越平坦 阻带最小衰减αT αT越大，阻带衰减越大 当幅度下降到0.707时，w=wp，此时αp＝3dB，称wp为3dB通带截止频率。 按图6-2的规定给出滤波器的一组技术指标后，下一步就是寻找一个频率响应符合允许指标的离散时间线性系统。这样，滤波器的设计问题便归结为数学逼近问题。显然，对于IIR系统，可以应用有理函数去逼近所希望的频率响应；对于FIR系统，则可用多项式逼近所希望的频率响应。 三.数字滤波器设计方法概述 设计IIR数字滤波器一般有以下两种方法： 1. 利用模拟滤波器的理论来设计:首先根据实际要求设计一个合适的模拟滤波器，然后将它转换成满足给定指标的数字滤波器，这种方法适合于设计幅频特性比较规则的滤波器，例如低通、高通、带通、带阻等。这种设计方法又可分为冲激响应不变法、阶跃响应不变法和双线性变换法等三种。 2.计算机辅助设计：这种方法直接在频域或者时域中进行数字滤波器设计，由于要联立方程，设计时需要计算机作辅助设计。这种方法就是使用最优化技术进行设计。这种设计方法在本书不讨论，可参考其它文献。 6.2 无限冲激响应数字滤波器的设计 一.设计的一般方法 IIR滤波器以模拟低通滤波器为基础的设计方法，为了设计其他的选频滤波器（高通，带通，带阻等），需要对低通滤波器进行频率转换，在设计过程中有两种不同的变换,频带变换和模拟/数字变换。根据这两种变换的先后次序，引出两种设计方法。 如图6－3所示。 图6-3 IIR滤波器的设计流程 二.脉冲响应不变法(Impulse Invariance) 1.变换原理 数字滤波器的单位脉冲响应序列h(n)正好等于模拟滤波器的冲激响应ha(t)的采样值，即h(n)=ha(nT),T为采样周期。如以Ha(s)及H（z）分别表示ha(t)的拉氏变换及h(n)的Z变换，即 数字滤波器的系统函数就是h(n)的Z变换H（z）。 下面我们分析从模拟滤波器到数字滤波器S平面和Z平面之间的映射关系。 设抽样信号 抽样信号的拉氏变换 序列h(n)的z变换 所以 图6-4 S平面与Z平面的映射关系 当σ＝0时，Ω：－Ωs／2～Ωs／2 →｜z|＝1 ∠z:-π～π；即S平面虚轴jΩ：－Ωs／2～Ωs／2映射到Z平面单位圆上。 当σ＜0时,r＜1,｜z|＜1，Ω:-Ωs／2～Ωs／2;即S左半平面的-Ωs／2～Ωs／2带映射到Z平面单位圆内。 当σ＞0时,r＞1,｜z|＞1;即S右半平面的-Ωs／2～Ωs／2带映