上课教材第6章无限脉冲响应数字滤波器题库.ppt

第6章 无限长单位冲激响应(IIR)滤波器;6.1 引言;一、滤波器的分类; 数字滤波器的理想幅频特性是周期为 的偶函数，只需在区间 上研究即可！ ;一、滤波器的分类;二、滤波器的技术指标;低通滤波器的幅度特性指标;通带最大衰减;三、IIR数字滤波器的技术特点;四、IIR数字滤波器的设计方法概述;该设计方法的实质：是一个S平面到Z平面的复变函数的映射变换 ，即有 注意， 该映射变换必须要满足以下2个条件： （1）S平面虚轴 jΩ必须映射到Z平面的单位圆ejω上 （2）S平面的左半平面必须映射到Z平面单位圆的内部|z|1 下面先讨论模拟滤波器的设计方法，再介绍由模拟滤波器变换为IIR数字滤波器的两种映射变换：冲激响应不变法和双线性变换法 。 ;6.2 模拟滤波器的设计;6.2 模拟滤波器的设计;巴特沃斯（Butterworth）滤波器：具有单调下降的幅度特性 在通带内有最大平坦的幅度特性。 切比雪夫（ Chebyshev ）滤波器：通带或阻带有波动，可以提高选择性。 贝赛尔（ Bessel ）滤波器：通带内有较好的线性相位特性 椭圆（Ellipse0滤波器：选择性相对其它三种是最好的, 但在通带和阻带内均有等波纹的幅频特性。 ;三、模拟低通滤波器的设计;按给定的技术指标，设计出模拟滤波器的系统函数;问题：如何根据Ha(s)Ha(-s) 求 Ha(s) ？;（1）由幅度平方函数求得象限对称的S 平面函数;例2 在MATLAB中设计一个4阶巴特沃斯（Butterworth）模拟低通滤波器，其3dB截止频率为0.3rad/s 。 MATLAB程序如下： [b,a]=butter(4,0.3,low,s); freqs(b,a)； 右图是所设计的 4阶巴特沃斯 （Butterworth） 模拟低通滤波器 的频率响应曲线。 ;例3 在MATLAB中设计一个6阶椭圆低通模拟滤波器，要求通带最大衰减3dB，阻带最小衰减 30dB，其截止频率为100rad/s 。 MATLAB程序如下： [b,a]=ellip(6,3,30,100,low,s); freqs(b,a); 右图是所设计的 6阶椭圆 低通模拟滤波器 的频率响应曲线。;四、 高通模拟滤波器的设计 方法：先设计出低通模拟滤波器，再通过频率变换转换成要求类型(高通,带通等)，而低通是通过幅度平方函数来设计的。;例4 在MATLAB中设计一个24阶贝赛尔（ bessel ）高通模拟滤???器，其3dB通带边界频率为100rad/s 。 MATLAB程序如下： [b,a]=besself(24,100,high,s); freqs(b,a)； 右图是所设计的 24阶贝赛尔（ bessel ） 模拟高通滤波器 的频率响应曲线。;五、带通模拟滤波器的设计 方法：先设计出低通模拟滤波器，再通过频率变换转换成所要设计的带通模拟滤波器。;例5 在MATLAB中设计一个20阶巴特沃斯带通模拟滤波器，其3dB通带频率为100---1000 rad/s 。 MATLAB程序如下： [b,a]=butter(20,[100,1000], s); freqs(b,a)； 右图是所设计的 20阶巴特沃斯 带通模拟滤波器 的频率响应曲线。;六、带阻模拟滤波器的设计 方法：先设计出低通模拟滤波器，再通过频率变换转换成所要设计的带阻模拟滤波器。;例6 在MATLAB中设计一个24阶切比雪夫I型的带阻模拟滤波器，其3dB阻带的上、下边界频率分别为100 rad/s和1000 rad/s 。 通带最大衰减为3dB 。MATLAB程序如下： [b,a]=cheby1(12,3, [100,1000], stop , s); freqs(b,a); 右图是所设计的 24阶切比雪夫I型 带阻模拟滤波器 的频率响应曲线。;6.3 冲激响应不变法;6.3 冲激响应不变法;一、设计原理 具体的转换过程如下： ;是 的单阶极点，得到;例1：;映射关系;, 平面与 平面之间的映射关系;由采样序列的Z变换与原模拟信号的拉氏变换的关系，得到数字滤波器与模拟滤波器的频率响应间的关系为;冲激响应不变法中的频响混叠现象 ;三、冲激响应不变法的特点主要特点与适用范围; ; ;6.4 双线性变换法;一、设计原理 为了克服冲激响应不变法的频率混叠现象，提出了双线性变换法，它采用非线性的频率压缩方法，首先将整个频率轴的频率范围压缩到-π/T～π/T之间，再用z=esT转换到Z平面上，使S平面与Z平面建立起一一对应的单值映射关系。具体的变换规则如下： （1）频率压缩：将整个S平面压缩变换到某个中介S1平面的一横带里（宽度为