《数字信号处理》实验指导书 .pdf

《数字信号处理》实验指导书 《数字信号处理》实验指导书 实验 1 数字滤波器的设计 实验序号：1 实验名称：数字滤波器的设计 适用专业：通信工程、电子信息工程 学 时 数：4 学时 一、实验目的 1．掌握双线性变换法设计 IIR 数字滤波器的具体设计方法及其原理和窗函数设计 FIR 滤波器的设计原理与基本方法。 2 ．观察双线性变换的频域特性。熟悉 Butterworth 滤波器的频率特性。 3 ．了解各种不同窗函数对滤波器性能的影响。 4 ．熟悉 Matlab 计算机编程。 二、实验原理 1．用双线性变换法设计 IIR 数字滤波器方法 （1）设计思想：将模拟滤波器转换成数字滤波器的实质是，用一种从 s 平面到 z 平面 的映射函数将Ha(s)转换成 H(z)。对这种映射函数的要求是：(1) 因果稳定的模拟滤波器转 换成数字滤波器，仍是因果稳定的。 (2)数字滤波器的频率响应模仿模拟滤波器的频响，s 平面的虚轴映射 z 平面的单位圆，相应的频率之间成线性关系。脉冲响应不变法和双线性变 换法都满足如上要求。 s 平面与 z 平面之间满足以下映射关系： 1 z  1 s  1 z 1 s 平面的虚轴单值地映射于 z 平面的单位圆上，s 平面的左半平面完全映射到 z 平面的单位 圆内。双线性变换不存在混叠问题。 双线性变换时一种非线性变换  tg ( ω / 2) ，这种非线性引起的幅频特性畸变可通过 预畸而得到校正。 （2 ）以低通数字滤波器为例，将设计步骤归纳如下： ·确定数字滤波器的性能指标：通带临界频率 f 、阻带临界频率 f ；通带内的最大衰 p s 减 A ；阻带内的最小衰减A ； p s · 确定相应的数字角频率， ωp=2πfp；ωs=2πfs； 第 1 页 《数字信号处理》实验指导书   tg ( / 2) ·计算经过预畸的相应模拟低通原型的频率， ； ·根据 Ωp和 Ωs 计算模拟低通原型滤波器的阶数 N，并求得低通原型的传递函数 Ha(s)； ·用上面的双线性变换公式代入 Ha(s)，求出所设计的传递函数 H(z)； ·分析滤波器特性，检查其指标是否满足要求。 2 ．用窗函数设计 FIR 滤波器的基本方法 FIR 滤波器的设计问题在于寻求一系统函数H (z ) ，使其频率响应H ( ej ω ) 逼近滤波器 要求的理想频率响应H d (ej  ) ，其对应的单位脉冲响应hd (n) 。 （1）设计思想：从时域从发，设计h(n) 逼近理想hd (n) 。设理想滤波器H d (ej  ) 的单 位脉冲响应为hd (n) 。以低通线性相位FIR 数字滤波器为例。  H d (ej  )  hd (n)ej n n 1  hd (n)   H d (ej  )ej n d 2  hd (n) 一般是无限长的，且是非因果的，不能直接作为 FIR 滤波器的单位脉冲响应。要 想得到一个因果的有限长的滤波器 h(n) ，最直接的方法是截断h(n)  h (n)w(n) ，即截取为 d 有限长因果序