实验二-IIR数字滤波器的设计.docxVIP

实验二IIR数字滤波器的设计 1、 实验目的 掌握脉冲响应不变法和双线性变换法设计 IIR数字滤波器的具体方法和 原理，熟悉双线性变换法和脉冲响应不变法设计低通、 带通IIR数字滤波器的计 算机编程； 观察双线性变换法和脉冲响应不变法设计的数字滤波器的频域特性，了 解双线性变换法和脉冲响应不变法的特点和区别； 熟悉Butterworth滤波器、Chebyshev滤波器和椭圆滤波器的频率特性。 2、 实验原理与方法 IIR数字滤波器的设计方法可以概括为如图所示， 本实验主要掌握IIR滤波 器的第一种方法，即利用模拟滤波器设计IIR数字滤波器，这是IIR数字滤波器 设计最常用的方法。利用模拟滤波器设计，需要将模拟域的 Ha(s)转换为数字域 H(z)，最常用的转换方法为脉冲响应不变法和双线性变换法。 脉冲响应不变法 用数字滤波器的单位脉冲响应序列 h(n)模仿模拟滤波器的冲激响应ha(t)， 让h(n)正好等于ha(t)的采样值，即 h(n) =ha( nT) 其中T为采样间隔。如果以Hk(s)及H(z)分别表示ha(t)的拉氏变换及h(n) 的Z变换，贝U 1力 2兀 H(z)|zj Ha(s-j k) T T 在MATLAB^，可用函数imp in var实现从模拟滤波器到数字滤波器的脉冲响 应不变映射。 双线性变换法 S平面与z平面之间满足下列映射关系 2 2 1 -z 2 1 -z 或 z二一 s T S平面的虚轴单值地映射于 S平面的虚轴单值地映射于z平面的单位圆上, s平面的左半平面完全映射 到z平面的单位圆内。双线性变换不存在频率混叠问题。 在MATLAB^，可用函数bilinear实现从模拟滤波器到数字滤波器的双线性 变换映射。 双线性变换是一种非线性变换，即I； - 2tan—，这种非线性引起的幅频特 T 2 性畸变可通过预畸变得到校正。 （3）设计步骤 IIR数字滤波器的设计过程中，模拟滤波器的设计是关键。模拟滤波器的设 计一般是采用分布设计的方式，这样设计原理非常清楚，具体步骤如前文所述。 MATLAB信号处理工具箱也提供了模拟滤波器设计的完全工具函数： butter、 chebyl、cheby2、ellip、besself。用户只需一次调用就可完成模拟滤波器的设 计，这样虽简化了模拟滤波器的设计过程，但设计原理却被屏蔽了。 模拟滤波器设计完成之后，利用impinvar或bilinear函数将模拟滤波器映 射为数字滤波器，即完成了所需数字滤波器的设计。 下图给出了实际低通、高通、带通和带阻滤波器的幅频特性和各截止频率的 含义。另外，为了描述过渡带的形状，还引入了通带衰减和阻带衰减的概念。 Wd Wc W咄 W W.i Wo WcK W Wd Wc W咄 W W.i Wo WcK W BPFBEF BPF 图实际滤波器的幅频特性和各截止频率的含义 通带衰减： :-p 二 20log 1 H (e)1 P |H(e?| 阻带衰减：:-s = 20log 1 H (e)1 dBdB在MATLAB言号处理工具箱中，通常用R＞和R dB dB 在MATLAB言号处理工具箱中，通常用 R＞和RS来表示a p和a s。 3、 实验内容 参照教材5.5节所述滤波器设计步骤，利用双线性变换法设计一个 Chebyshev I型数字高通滤波器，观察通带损耗和阻带衰减是否满足要求。已知 滤波器的指标为 fp=0.3kHz, a P=1.2dB, fs=0.2kHz, a s=20dB, T=1ms 已知 fp=0.2kHz, a p=1dB, fs=0.3kHz, a s=25dB, T=1ms 分别用脉冲 响应不变法和双线性变换法设计一个 Butterworth数字低通滤波器，观察所设计 数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两 种方法的优缺点。 设计一个数字带通滤波器,通带范围为 0.25 n ~0.45 n ,通带内最大 衰减3dB, 0.15 n以下和0.55 n以上为阻带,阻带内最小衰减为 15dB,试采用 Butterworth或ellip (椭圆)模拟低通滤波器设计。 利用双线性变换法设计一个带宽为 0.08 n的10阶椭圆带阻滤波器以 滤除数字频率为0.44 n的信号，选择合适的阻带衰减值，画出幅度响应。产生 下面序列的201个样本 x(n) =sin0.44二n , n=0 , 2,…,200 并将它通过这个带阻滤波器进行处理(filter 函数)，讨论所得到的结果。 4、 实验报告 简述实验目的和实验原理。 按实验步骤附上所设计的滤波器传递函数 H(z)及相应的幅频特性曲线，定 性分析所得到的图形，判断设计是否满足要求。 总结脉冲响应不变法和双线性变换法的特点及设计全过程。