数字信号处理尹为民21-62IIRDF设计(一)课件教学.pptVIP

6.2 IIR数字滤波器设计 基本概念 模拟低通滤波器的设计 冲激响应不变法 双线性变换法 设计IIR滤波器的频率变换法 一、 基本概念 1 滤波器的技术指标 1 滤波器的技术指标 2 滤波器的主要参数 3 IIR滤波器设计方法 3 IIR滤波器设计方法 3 IIR滤波器设计方法 二、模拟低通滤波器的设计 1 实际低通滤波器的技术指标 1 实际低通滤波器的技术指标 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 2 幅度平方函数的确定 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 3 由幅度平方函数确定系统函数 4 巴特沃思低通滤波器的实际设计步骤 5 MATLAB实现 由图可以看出：低通滤波器的实际特性除存在通带和阻带外，在通带和阻带之间设置了一个过渡带，而不是突然下降，但在通带内不是完全平直，而是近似于理想的幅度特性，并与理想特性的偏差在规定的范围之内，在阻带内幅度特性也不是零值，而是衰减至规定的偏差范围内，对过渡带内的幅度衰减可根据设计要求进行调整。 常用的模拟原型滤波器有巴特沃思（Butterworth）滤波器、切比雪夫（Chebyshev）滤波器、椭圆（Ellipse）滤波器、贝塞尔（Bessel）滤波器等。这些滤波器都有严格的设计公式，现成的曲线和图表供设计人员使用。这些典型的滤波器各有特点：巴特沃思滤波器具有单调下降的幅频特性；切比雪夫滤波器的幅频特性在通带或者在阻带有波动，可以提高选择性；贝塞尔滤波器通带内有较好的线性相位特性；椭圆滤波器的选择性相对前三种是最好的, 但在通带和阻带内均为等波纹幅频特性。这样根据具体要求可以选用不同类型的滤波器。 Butterworth模拟滤波器幅频的最高值为1，故在画幅频图时没有做归一化操作。 重点：滤波器的技术指标 输入序列的频谱X(ejω)经过滤波后，变为X(ejω)H(ejω)。如果|H(ejω)|的值在某些频率上是比较小的，则输入信号中的这些频率分量在输出信号中将被抑制掉。因此，只要按照输入信号频谱的特点和处理信号的目的，适当选择H(ejω)，使得滤波后的X(ejω)H(ejω)符合人们的要求，这就是数字滤波器的滤波原理。 虽然给出了通带的容限?1及阻带的容限?2，但是，在具体技术指标中往往使用通带允许的最大衰减（纹波）Ap和阻带应达到的最小衰减As描述。 式中，假定|H(ej0)|=1(已被归一化)。例如|H(ejω)|在ωp处满足|H(ejωp)|=0.707，则Ap=3dB；在ωs处满足|H(ejωs)|=0.001，则As=60 dB（参考图6-2）。(注：lg是log10的规范符号表示。) 巴特沃思逼近又称最平幅度逼近。巴特沃思低通滤波器幅度平方函数定义为 式中，N为正整数，代表滤波器的阶数。当Ω=0时，|Ha(j0)|=1； 当Ω=Ωc时，|Ha(jΩc)|=1/ =0.707，20lg|Ha(j0)/Ha(jΩc)|=3 dB， Ωc为3 dB截止频率。当Ω=Ωc时，不管N为多少，所有的特性曲线都通过－3 dB点，或者说衰减为 3 dB。 [例]试设计巴特沃思低通滤波器，要求通带截止频率fp=2kHz，通带衰减 ? p ? 3dB ，阻带截止频率fs=4KHz ，阻带衰减?s ? 15dB。 [解] 代入 s = s / ?c 得到所需的巴特沃思低通滤波器 其中 ?c = 2? f c = 4000? (rad/s) [例]试设计巴特沃思低通滤波器 Ha(s) ，使其满足下列技术指标： [解] 取N = 6 [num,den]=butter(N,wc,s) 确定阶数为N，3-dB截频为wc(radian/s)的Butterworth filter分子和分母多项式。s 表示模拟域。 [N,wc]=buttord(wp,ws,Ap,As,s) 确定模拟Butterworth filter的阶数N和3-dB截频wc。wc是由阻带参数确定的。 s 表示模拟域。 [z,p,k]=buttap(N) 确定N阶归一化的Butterworth filter 的零点、极点和增益(gain) 例：设计满足下列条件的模拟Butterworth低通滤波器 fp=10kHz, fs=15kHz, ?p=1dB, ?s=15dB Wp=2*pi*10000;Ws=3*pi*10000;Ap=1;As=15; [N,Wc]=buttord(Wp,Ws,Ap,A