数字信号处理实验指导书167.docVIP

《数字信号处理》实验指导书 实验数字滤波器的设计 实验序号：1 实验名称：数字滤波器的设计 1．掌握双线性变换法设计IIR数字滤波器的具体设计方法及其原理和窗函数设计FIR滤波器的设计原理与基本方法。 2．观察双线性变换的频域特性。熟悉Butterworth滤波器的频率特性。 3．了解各种不同窗函数对滤波器性能的影响。 4．熟悉Matlab计算机编程。 二、实验原理 1．用双线性变换法设计IIR数字滤波器方法 （1）设计思想：将模拟滤波器转换成数字滤波器的实质是，用一种从s平面到z平面的映射函数将Ha(s)转换成H(z)。对这种映射函数的要求是：(1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。 (2)数字滤波器的频率响应模仿模拟滤波器的频响，s平面的虚轴映射z平面的单位圆，相应的频率之间成线性关系。脉冲响应不变法和双线性变换法都满足如上要求。 s平面与z平面之间满足以下映射关系： s平面的虚轴单值地映射于z平面的单位圆上，s平面的左半平面完全映射到z平面的单位圆内。双线性变换不存在混叠问题。 双线性变换时一种非线性变换，这种非线性引起的幅频特性畸变可通过预畸而得到校正。 （2）以低通数字滤波器为例，将设计步骤归纳如下： ·确定数字滤波器的性能指标：通带临界频率fp、阻带临界频率fs；通带内的最大衰减Ap；阻带内的最小衰减As； ·?确定相应的数字角频率，ωp=2πfp；ω=2πfs； ·计算经过预畸的相应模拟低通原型的频率，； ·根据Ωp和Ωs计算模拟低通原型滤波器的阶数N，并求得低通原型的传递函数Ha(s)； ·用上面的双线性变换公式代入Ha(s)，求出所设计的传递函数H(z)； ·分析滤波器特性，检查其指标是否满足要求。 2．用窗函数设计FIR滤波器的基本方法 FIR滤波器的设计问题在于寻求一系统函数，使其频率响应逼近滤波器要求的理想频率响应，其对应的单位脉冲响应。 （1）设计思想：从时域从发，设计逼近理想。设理想滤波器的单位脉冲响应为。以低通线性相位FIR数字滤波器为例。 一般是无限长的，且是非因果的，不能直接作为FIR滤波器的单位脉冲响应。要想得到一个因果的有限长的滤波器h(n)，最直接的方法是截断，即截取为有限长因果序列，并用合适的窗函数进行加权作为FIR滤波器的单位脉冲响应。按照线性相位滤波器的要求，h(n)必须是偶对称的。对称中心必须等于滤波器的延时常数，即 用矩形窗设计的FIR低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象，且最大波纹大约为幅度的9%，这个现象称为吉布斯（Gibbs）效应。为了消除吉布斯效应，一般采用其他类型的窗函数。 （2）典型的窗函数 ·矩形窗(Rectangle Window) 其频率响应和幅度响应分别为： ， ·三角形窗(Bartlett Window) 其频率响应为： ·汉宁(Hanning)窗，又称升余弦窗 其频率响应和幅度响应分别为： ·汉明(Hamming)窗，又称改进的升余弦窗 其幅度响应为： ·布莱克曼(Blankman)窗，又称二阶升余弦窗 其幅度响应为： （3）利用窗函数设计FIR滤波器的具体步骤如下： ·按允许的过渡带宽度△ω及阻带衰减AS，选择合适的窗函数，并估计节数N： 其中A由窗函数的类型决定。 ·由给定的滤波器的幅频响应参数求出理想的单位脉冲响应。 ·确定延时值 ·计算滤波器的单位取样响应，。 ·验算技术指标是否满足要求。 三、主要实验仪器及材料 微型计算机、Matlab6.5以上版本的编程环境。 四、实验内容 1．知识准备 在实验前复习数字信号处理理论课中有关滤波器设计的知识，认真阅读本实验的原理部分。熟悉滤波器的特性，掌握设计数字滤波器的具体步骤。 2．编制用双线性变换法设计IIR数字滤波器的程序。采样周期、通带和阻带临界频率以及相应的衰减等参数在程序运行时输入；根据这些输入参数，计算阶数N、传递函数；输出分子分母系数；绘制幅频