无限冲激响应(IIR)数字滤波器设计zhjConcise剖析.ppt

求出： 幅频响应曲线见书 P284 由： 给出 数字高通 的技术要求 得到 模拟高通 的技术要求 得到 模拟低通 的技术要求 设计出 得到 模拟高通转移 函数 最后得到 数字高通转移 函数 数字高通滤波器设计步骤 6.6 数字高通, 带通及带阻滤波器的设计 对 带通（BP）、带阻（BS）数字滤波器的设计，只需改变图中 Step2 和 Step4： 带通 带阻 例： 设计一 IIR BP DF,要求： 通带频率范围： 300Hz~ 400Hz ; 阻带频率范围：200Hz、500Hz 要求： 按上述转换办法，可以求出： 1．buttord.m 确定 LP DF、或 LP AF的阶次； （1） [N, Wn] = buttord(Wp, Ws, Rp, Rs)； 与本章内容有关的MATLAB文件 对应 数字滤波器。其中 Wp, Ws分别是通带和阻带的截止频率，其值在 0～1 之间，1对应抽样频率的一半(归一化频率)。对低通和高通，Wp, Ws都是标量，对带通和带阻，Wp, Ws是1×2的向量。Rp, Rs 分别是通带和阻带的衰减(dB)。N是求出的相应低通滤波器的阶次，Wn是求出的3dB频率，它和Wp稍有不同。 （2）[N, Wn] = buttord(Wp, Ws, Rp, Rs,‘s’)： 对应模拟滤波器，式中各个变量的含意和格式（1）相同，但Wp, Ws及Wn的单位为弧度/秒，因此，它们实际上是频率。 2．buttap.m 设计模拟低通(Butt)原型滤波器。 [z, p, k]=buttap(N): N是欲设计的低通原 型滤波器的阶次，z, p, k是设计出的极点、 零点及增益。 [B, A]=lp2lp(b, a, Wo), [B, A]=lp2hp(b, a, Wo) [B, A]=lp2bp(b, a, Wo, Bw), [B, A]=lp2bs(b, a, Wo, Bw) b, a 是AF LP 的分子、分母的系数向量，B, A是转换后的的分子、分母的系数向量；在（1）中，Wo是低通或高通滤波器的截止频率；在（2）中，Wo是带通或带阻滤波器的中心频率，Bw是其带宽。 3．lp2lp.m、lp2hp.m、lp2bp.m, lp2bs.m 将模拟低通原型转换为实际的低通、高通、带通及带阻滤波器。 （1） （2） 4．bilinear.m ：双线性变换，由模拟滤波器 得到数字滤波器。 [Bz, Az]=bilinear(B, A, Fs) 式中B, A分别是G（s）的分子、分母多项式 的系数向量，Bz, Az分别是H(z)的分子、分 母多项式的系数向量，Fs是抽样频率。 5．butter.m 本文件可用来直接设计Butterworth数字滤波器，实际上它把 Buttord buttap、 lP2lp Bilinear 等文件都包含了进去，从而使设计过程更简捷。 格式（1）～（3）用来设计数字滤波器，B，A分别是H（z）的分子、分母多项式的系数向量，Wn是通带截止频率，范围在0～1之间。若Wn是标量，（1）用来设计低通数字滤波器，若Wn是1×2的向量，则（1）用来设计数字带通滤波器；（2）用来设计数字高通滤波器；（3）用来设计数字带阻滤波器，显然，这时的Wn是1×2的向量；格式（4）用来设计模拟滤波器。 （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’)。 6．cheb1ord.m 求Cheb-Ⅰ型滤波器的阶； 7.cheb1ap.m 设计原型低Cheb-I型模拟滤波器； 8．cheby1.m 直接设计数字Cheb-Ⅰ滤波器。 以上三个文件的调用格式和对应的 Butterworth滤波器的文件类似。 9．cheb2ord.m; 10. ellipord.m; 11.cheb2ap.m; 12. ellipap.m; 13.besselap.m; 14