DSPIIR滤波器设计概要.ppt

数字信号处理第八章 IIR滤波器设计 授课教师：胡双红 联系QQ长沙理工大学 计算机与通信工程学院 第一次课 概述 模拟原型滤波器 § 8.1 概述 IIR设计思想 设计步骤 设计指标转换 由幅度平方函数?系统函数 两种基本方法：复值映射 设计步骤： 确定技术要求：由具体应用条件决定 设计模拟低通滤波器 滤波器变换--?数字低通滤波器 频带变换--?其它数字滤波器 主要问题 对IIR数字滤波器的相位特性没有任何控制能力. 因此IIR滤波器设计仅作为幅度设计 模拟滤波器的设计指标 相对线性标尺 低通滤波器在幅度平方响应上的技术指标给出为 由幅度平方函数 确定系统函数 §8.2 原型模拟滤波器特性 IIR滤波器设计方法依赖于已有的模拟滤波器得到数字滤波器，这些模拟滤波器被称为原型滤波器( prototype filters). 实际中广泛应用的三种原型滤波器 Butterworth低通 Chebyshev 低通(Type I and II) Elliptic(椭圆) 低通 Butterworth低通滤波器 表征性质：幅度响应在通带和阻带都是平坦的 定义：幅度平方响应为 Butterworth滤波器性质 幅度响应 |Ha(0)|2 =1, |Ha(jΩc) |2 =0.5, 即对全部 N (在处)都有3dB衰减，称 Ωc为Butterworth低通滤波器的3dB带宽 |Ha(jΩ)|2 随 Ω的增加单调下降 当N→∞时趋近于理想低通滤波器 Ω=0处是最大平坦，所有阶的导数都存在且等于零 确定系统函数Ha(s) |Ha(jΩ)|2 Ω=s/j =Ha(s) Ha(-s)的极点 2N个极点均匀分布在半径为Ωc的圆上，相隔 pi/N N 为奇数时， pk= Ωc ej2pik/N N 为偶数时， pk= Ωc ej(pi/2N+kpi/N) 极点对虚轴(jΩ)对称分布 极点永远不会落在虚轴上，且仅当N为奇数时才会落在实轴上 选取极点全部落在左半平面可得到因果稳定的Ha(s) Butterworth滤波器的系统函数： Butterworth滤波器的设计步骤： N的求解原理 求出归一化系统函数： 归一化N阶巴特沃斯多项式系数 b) 求出极点（左半平面） MATLAB实现 函数 [z,p,k] = buttap(N) 设计一个归一化(Ωc=1)的N阶 Butterworth 模拟原型低通滤波器 Z: 零点; p: 极点; k: 增益值 函数 [b,a] = u_buttap(N,Omegac) 设计非归一化的直接型巴特沃兹模拟原型滤波器. 函数 [C,B,A] = sdir2cas(b,a) 将模拟直接型滤波器转换为模拟级联型滤波器 MATLAB 实现 函数 [b,a] = afd_butt(Wp,Ws,Rp,As) 根据给定技术指标设计模拟低通Butterworth滤波器 函数 [db,mag,pha,w]=freqs_m(b,a,wmax) 以绝对值和相对值dB 尺度计算幅度响应以及相位响应 响应一直计算到最大频率Ωmax为止 模拟滤波器的脉冲响应可用函数impulse计算 调用格式：[ha,x,t]=impulse(b,a) 设计满足下面要求的低通Butterworth滤波器 Ωp=0.2π，Rp=7dB，Ωs=0.3 π ，As=16dB 调用MATLAB函数实现求解 wp=0.2*pi; ws=0.3*pi; Rp=7;As=16; Ripple=10^(-Rp/20); Atten=10^(-As/20); 设计模拟滤波器 [b,a]=afd_butt(wp,ws,Rp,As); *** Butterworth Filter Order = 6 转换为二阶节的级联 [C,B,A]=sdir2cas(b,a) C = 0.4985 B = 0 0 1 A = 1.0000 0.4985 0.2485 0 1.0000 0.4985 Chebyshev低通滤波器 Chebyshev-I型滤波器 在通带内有等波纹响应 Chebyshev-II型滤波器 在阻带内有等波纹响应 Butterworth 滤波器 通带和阻带内都具有单调(monotonic )响应 达到相同的性能指标， Chebyshev滤波器会比Butterworth需要较少的阶次. Chebyshev-I型滤波器幅度平方响应 TN(x) 的关键性质 Chebyshev滤波器的等波纹响应由多项式TN(x)引起 0x1时, TN(x)在–1和 1之间起伏 1xinfinity时, TN(x) 单调增加至