第6章 无限脉冲响应数字滤波器的设计1.pptVIP

第6章 无限脉冲响应数字滤波器的设计 6.1 数字滤波器的基本概念 6.2 模拟滤波器的设计 6.3 用脉冲响应不变法设计IIR数字低通滤波器 6.4 用双线性变换法设计IIR数字 低通滤波器 6.5 数字高通、带通和带阻滤波器的设计 6.6 IIR 数字滤波器的直接设计法 习题 1. 4. 5. 7. 11. IIR和FIR 系统函数分别为： 假设数字滤波器的频率函数H(e jω)用下式表示： 设计步骤：先设计模拟滤波器得到传输函数Ha(s)，然后将Ha(s)按某种方法转换成数字滤波器的系统函数H(z)。 6.2 模拟滤波器的设计 巴特沃斯低通滤波器的幅度平方函数|Ha(jΩ)|2用下式表示： 例6.2.1 已知通带截止频率fp=5kHz，通带最大衰减αp=2dB，阻带截止频率fs=12kHz，阻带最小衰减αs=30dB，按照以上技术指标设计巴特沃斯低通滤波器。 先求3dB截止频率Ωc。按照(6.2.17)式，得到： 3.切比雪夫滤波器的设计方法 图6.2.6示出了阶数N=0,4,5时的切比雪夫多项式特性。 由图可见： (1)切比雪夫多项式的过零点在|x|≤1的范围内； (2)当|x|1时，|CN(x)|≤1,在|x|1范围内具有等波纹性； (3)当|x|1时，CN(x)是双曲线函数，随x单调上升。 按照(6.2.19)式，平方幅度函数与三个参数ε,Ωp和N有关。其中ε与通带内允许的波动大小有关，定义允许的通带波纹δ用下式表示： 设阻带的起始点频率(阻带截止频率)用Ωs表示，在Ωs处的A2(Ωs)用(6.2.19)式确定： 3dB截止频率用Ωc表示， 例6.2.2设计低通切比雪夫滤波器，要求通带截止频率fp=3kHz，通带最大衰减αp=0.1dB，阻带截止频率fs=12kHz，阻带最小衰减αs=60dB。 解 (1) 滤波器的技术要求： (2) 求阶数N和ε： 6.3 用脉冲响应不变法设计IIR数字低通滤波器 设模拟滤波器的传输函数为Ha(s),相应的单位冲激响应是ha(t) ，由拉普拉斯变换， 对ha(t)进行等间隔采样，采样间隔为T，得到： 另外，z=esT是一个周期函数，可写成 如果模拟滤波器二阶基本节的形式为 用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z)。 6.4 用双线性变换法设计IIR数字低通滤波器 (6.4.2) 利用双线性变换法转换，数字滤波器的系统函数H2(z)为 ③设计巴特沃斯低通滤波器。先计算阶数N及3dB截止频率Ωc。 ② 模拟低通的技术指标为 ③ 设计巴特沃斯低通滤波器。 有关的MATLAB文件 1.buttord.m (1)[N,Wn]=buttord(Wp,Ws,Rp,Rs)；数字滤波器，归一化频率（0~1，1对应抽样频率的一半） (2) [N,Wn]=buttord(Wp,Ws,Rp,Rs,’s’)；模拟滤波器，频率单位rad/s; N阶次,Wn——3dB频率 2.buttap.m [z,p,k]=buttap(N) ;原型滤波器G(p),N阶次，z,p,k零、极点、增益 3.lp2lp.m [B,A]=lp2lp(b,a.Wo) 4.lp2hp.m [B,A]=lp2hp(b,a,Wo) 5.lp2bp.m [B,A]=lp2bp(b,a,Wo,Bw) 6.lp2bs.m [B,A]=lp2bs(b,a,Wo,Bw) 7.bilinear.m [Bz,Az]=bilinear(B,A,Fs) 8.butter.m ;直接设计Butterworth数字滤波器 (1)[B,A]=butter(N,Wn); (2)[B,A]=butter(N,Wn,’high’); (3)[B,A]=butter(N,Wn,’stop’); (4)[B,A]=butter(N,Wn,’s’); 9.cheb1ord.m 10.cheb1ap.m 11.cheby1.m 12cheb2ord.m 13.ellipord.m 14.cheb2ap.m 15.ellipap.m 16.besselap.m 17.cheby2.m 18.ellip.m 19.besself