IIR滤波器的设计方法IIR滤波器的设计方法.doc

实验六 IIR滤波器的设计方法 实验学时：2学时 实验类型：设计、验证 一、实验目的 IR数字滤波器的常用指标和设计过程的理解； 响应不变法和双线性变换法设计IIR滤波器的原理，设计方法，步骤； 现象，比较脉冲响应不变法和双线性变换法设计IIR滤波器的特点 熟悉Matlab信号处理工具箱中的常用函数 二、实验原理与方法 IIR数字滤波器设计步骤： （一）先设计模拟低通原型滤波器 （二）AF数字化为DF 1． 脉冲响应不变法 用数字滤波器的单位脉冲响应序列模仿模拟滤波器的冲激响应,让正好等于的采样值，即，其中为采样间隔，如果以及分别表示的拉式变换及的Z变换，则 2．双线性变换法 ??? S平面与z平面之间满足以下映射关系： （***） s平面的虚轴单值地映射于z平面的单位圆上，s平面的左半平面完全映射到z平面的单位圆内。双线性变换不存在混叠问题。 ??? 双线性变换是一种非线性变换 ，这种非线性引起的幅频特性畸变可通过预畸而得到校正。双线性变换法特别适合用于设计常见的选频性滤波器，下面以双线性变换法设计低通数字滤波器为例，将设计步骤归纳如下： （1）根据需要确定数字滤波器的性能指标：通带截止频率、阻带截止频率、通带最大衰减αp、阻带内的最小衰减αs，采样频率 ； (2)?确定相应的数字角频率，；； (3)?计算经过预畸的相应模拟低通原型的截止频率， ，； (4)根据Ωp，αp和Ωs，αs计算模拟低通原型滤波器的传递函数； (5)用上述（***）表示的双线性变换公式代入，求出； ?分析滤波器特性，检查其指标是否满足要求。 3、IIR低通、高通、带通数字滤波器设计采用双线性原型变换公式： 变换类型 ????????? 变换关系式 ??? 备?? 注 低通 高通 ? 带通 为带通的上下边带临界频率 相关matlab函数介绍 1、Matlab信号处理工具箱中提供了设计巴特沃思模拟滤波器的函数，格式如下： 用于设计巴特沃思模拟滤波器Ha(s)设计的函数 模拟滤波器传输函数表示为： 其中M是向量B的长度，N是向量A的长度，分别是B和A的第k个元素。 （1） 用于根据设计指标计算巴特沃思ALPF幅度平方函数的参数N和3dB截止频率Wc (即本书中的符号)。其中，Wp 和Ws分别是滤波器的通带截止频率和阻止截止频率，单位为rad/s；Rp和Rs分别是通带最大衰减系数和阻带最小衰减系数，单位为dB。 （2） 用于根计算巴特沃思ALPF系统函数Ha(s)中分子和分母多项式系数向量B和A，其中N和分别是滤波器的阶和3dB截止频率，返回向量B和A，其中和分别是上面的表示式中的分母和分子系数。ftype缺省时，设计低通滤波器；ftype=high时，设计高通滤波器；ftype=stop时，设计带阻滤波器，此时为两个元素的向量，且ftype缺省时，设计带通滤波器，带通的频率区间为。S缺省时，设计数字滤波器。 例如：设计一个满足下列指标要求的巴特沃思模拟低通滤波器。指标 Matlab程序： 运行结果如下： N=6 Wc=2.2496 B= 0 0 0 0 0 0 129.5917 A=1.0000 8.6916; 37.7720; 104.0667;191.1447 ;222.5973 ; 129.5917 即： 幅频响应曲线（略）。 2、Matlab中设计数字滤波器的函数都是采用双线性变换法，将模拟滤波器转换为数字滤波器。这些函数及其凋用格式如下（以巴特沃思数字滤波器为例） 数字滤波器的系统函数为： 其中，分别是向量B和A的第k个元素。 （1） 该格式用于计算巴特沃思数字滤波器的阶N和3dB截止频率的归一化值（关于归一化）。调用参数Wp 和Ws分别是数字滤波器的通带截止频率和阻带截止频率的归一化值（关于归一化），要求和，其中1表示数字频率（对应模拟频率，为采样频率）。和分别是通带最大衰减和阻带最小衰减，单位dB。当时，设计高通滤波器；当和是二元向量时，设计带通（）或带阻（）滤波器，这时返回参数也是二元向量。 （2） 该格式用于计算巴特沃思模拟滤波器系统函数中分子和分母多项式系数向量B和A。调用参数N 和WC分别是巴特沃思数字滤波器的阶和3dB截止频率归一化值（关于归一化）。当ftype缺省时，设计低通滤波器；当ftype=high时，设计高通滤波器；ftype=stop时，设计带阻滤波器，此时为二元向量，和分别是带阻滤波器的通带3dB下，上截止频率的归一化值（关于归一化）；为向量，且ftype缺省时，设计带通滤波器，带通的频率区