精选实验11_iir滤波器.pptVIP

* 1、简述实验目的及原理。 2、给出上述四个滤波器的频响特性，并对结果加以分析。 3、根据对心电信号进行处理的程序，写出用双线性变换法设计的巴特沃斯低通IIR数字滤波器的系统函数。 五、实验报告及要求 * * * * * * 实验十一 用双线性变换法设计  IIR数字滤波器 * 掌握用双线性变换法设计IIR DF的原理及具体设计方法，熟悉用双线性变换法设计IIR DF的计算机编程。 观察用双线性变换法设计的DF的频响特性，了解双线性变换法的特点。 熟悉用双线性变换法设计BW和CB型DF的全过程。 一、实验目的 * 二、实验设备 2、MATLAB6.5 软件 1、计算机 * 为了克服冲激响应不变法产生的频率混叠现象，这是从S平面到Z平面的标准变换z＝esT的多值对应关系导致的，为了克服这一缺点，产生了双线性变换法。 三、实验原理 1．双线性变换法 * 双线性变换法的映射函数： 分式展开 固定其中一个变量，则另一个是线性的。或者说，这种变换对于s和z是双线性的。 三、实验原理 * （1）把s=σ+jΩ带入上式得： 因此： （2）整个左半平面映射成单位圆的内部，因此这是一个稳定的变换。 （3）虚轴以一一对应的方式映射成单位圆，因此在频域中不会出现混叠。 三、实验原理 * 由于幅度为1，把ω看做Ω的函数求解，得到： 这说明Ω和ω的关系是非线性的，但是没出现混叠。在把Ω变换为ω时产生了非线性畸变。为了克服它带来的问题，通常要使Ω按上式预修正，以抵消畸变的影响。 （4）根据上述变换关系，可以写出： 上式就是模拟滤波器和经采样后的数字滤波器之间的变换关系。 三、实验原理 * 求最小阶数N，ωc Buttord,cheblord Cheb2ord,ellipord 合为一步的设计函数 Butter，cheb1，cheb2，ellip，besself 模拟低通滤波器原型设计 buttap,cheb1ap,cheb2ap, besselap,ellipap函数 频率变换（变为高通、带通、带阻等）lp2lp,lp2hp,lp2bp,lp2bs 模拟-数字 bilinear，impinvar 滤波器 系数B，A 滤波器系数B，A 典型设计 直接设计 设计指标 2. IIR数字滤波器设计流程 设计流程图 三、实验原理 * 利用模拟滤波器设计IIR数字滤波器的步骤 2) 将数字滤波器的技术指标转换成模拟滤波器的技术指标。 确定数字滤波器的技术指标： 或 脉冲响应不变法: 双线性变换法: 通带截止频率ωp(ωl ωu)、阻带截止频率ωs(ωs1 ωs2) 通带衰减αp 、阻带衰减αs。 三、实验原理 * 4) 设计归一化低通滤波器Ga(p) 。 6) 将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字滤波器系统函数H(z)。 脉冲响应不变法或双线性变换法 3) 将模拟滤波器的技术指标转成模拟低通滤波器的技术指标。 频率变换 5) 将Ga(p)转成模拟滤波器Ha(s) 。 去归一化 + 频率变换 三、实验原理 * 巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器 巴特沃斯 切比雪夫I型 切比雪夫II型 椭圆滤波器 通带阻带都单调 通带有波纹 而阻带单调 通带单调， 阻带有波纹 通带阻带 都有波纹 3.典型的模拟滤波器 三、实验原理 * 阶数预测函数： buttord, cheb1ord, cheb2ord, ellipord 直接滤波器设计函数： besself,butter,cheby1, cheby2, ellip 低通模拟原型函数： besselap, buttap, cheb1ap, cheb2ap, ellipap 频率变换函数： lp2bp, lp2bs, lp2hp, lp2lp 滤波器离散函数： bilinear, impinvar 4. Matlab中的IIR滤波器相关函数 三、实验原理 * [n,wn]=buttord(wp,ws,Rp,Rs,’s’); 低通:wpws，其中0≤wp、ws≤1，1对应频率为0.5fs 高通: wpws 带通: wp、ws为二元矢量 带阻: wp、ws为二元矢量 [n,wn]=cheb1ord(wp,ws,Rp,Rs,’s’); [n,wn]=cheb1ord (wp,ws,Rp,Rs); 1)确定阶次与截止频率 [n,wn]=buttord(wp,ws,Rp,Rs); wp与ws的单位为rad, wpws, s指模拟滤波器设计 模拟角频率 w=f/(FS/2) 三、实验原理 * 2)模拟低通原型 [z,p,k]