无限冲激响应IIR滤波器设计.pptVIP

* 非线性关系，但是一对一的转换 * * * Step1. Step2. 频率转换： Step3. Step4. * 所以： So：系数 可省略，双线性Z变换可定义为 ： 这一组定义和前面的定义，对最后的 DF 而言，结果是一样的，区别是中间设计的 AF， 由于缺少了频率定标，将无法给出符合要求的幅频响应。 * 例6.5.1： 设计 IIR LP DF，给定 及 可求出： 【解】 (1) (2) * 求出： 幅频响应曲线见书 P.245 由： * 1．buttord.m 确定 LP DF、或 LP AF的阶次； （1） [N, Wn] = buttord(Wp, Ws, Rp, Rs)； （2）[N, Wn] = buttord(Wp, Ws, Rp, Rs,‘s’)： 与本章内容有关的MATLAB文件 (1)对应 数字滤波器。其中 Wp, Ws分别是通带和阻带的截止频率，其值在 0～1 之间，1对应抽样频率的一半(归一化频率)。对低通和高通，Wp, Ws都是标量，对带通和带阻，Wp, Ws是1×2的向量。Rp, Rs 分别是通带和阻带的衰减(dB)。N是求出的相应低通滤波器的阶次，Wn是求出的3dB频率，它和Wp稍有不同。 (2)对应模拟滤波器，各变量含意和（1）相同，但Wp, Ws及Wn的单位为弧度/秒，它们实际上是频率。 * 2．buttap.m 设计模拟低通(Butt)原型滤波器。 [z, p, k]=buttap(N): N是欲设计的低通原型滤波器的阶次，z, p, k是设计出的极点、零点及增益。 3．lp2lp.m、lp2hp.m、lp2bp.m, lp2bs.m 将模拟低通原型转换为实际的低通、高通、带通及带阻滤波器。 b, a 是AF LP 的分子、分母的系数向量，B, A是转换后的的分子、分母的系数向量；(1)中，Wo是低通或高通滤波器的截止频率； [B, A]=lp2lp(b, a, Wo) [B, A]=lp2hp(b, a, Wo) (1) [B, A]=lp2bp(b, a, Wo , Bw) [B, A]=lp2bs(b, a, Wo , Bw) (2) (2)中Wo是带通或带阻滤波器中心频率，Bw是其带宽。 * 4．bilinear.m ：双线性变换，由模拟滤波器 得到数字滤波器。 [Bz, Az]=bilinear(B, A, Fs) 式中B, A分别是G（s）的分子、分母多项式 的系数向量，Bz, Az分别是H(z)的分子、分 母多项式的系数向量，Fs是抽样频率。 * 5．butter.m 用来直接设计Butterworth数字滤波器，实际上它把 buttord.m, buttap.m, lp2lp.m, bilinear.m等文件都包含了进去，从而使设计过程更简捷。 格式(1)～(3) 用来设计数字滤波器，B，A分别是H(z)的分子、分母多项式的系数向量，Wn是通带截止频率，范围在0～1之间。若Wn是标量，(1)用来设计低通数字滤波器，若Wn是1×2的向量，则(1) 用来设计数字带通滤波器； (2)用来设计数字高通滤波器； (3) 用来设计数字带阻滤波器，显然，这时的Wn是1×2的向量；格式(4) 用来设计模拟滤波器。 [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.7.1(例6.5.1) clear all; fp=100;fs=300;Fs=1000;rp=3;rs=20; wp=2*pi*fp/Fs;ws=2*pi*fs/Fs; Fs=Fs/Fs; % let Fs=1 % Firstly to finish frequency prewarping ; wap=tan(wp/2);was=tan(ws/2); % [n,wn]=buttord(wap,was,rp,rs,s) % Note: s! [z,p,k]=buttap(n); % [bp,ap]=zp2tf(z,p,k) % [bs,as]=lp2lp(bp,ap,wap) % % Note: s=(2/Ts)(z-1)/(z+1);Ts=1,that is 2Fs=1,Fs=0.5; [bz,az]=bilinear(bs,as,Fs/2) % [h,w]=freqz(bz,az,256,Fs*1000); plot(w,abs(h));grid on; 设计 IIR LP DF， * 例6.7.1(例6