《数字信号处理教学课件》IIR滤波器.ppt

（2）从低通到高通的转换 lp2hp函数可将截止频率为1rad/s的模拟低通滤波器原型变换成截止频率为Wn的高通滤波器。 [bt，at]=lp2hP(b，a，Wn) （3）从低通到带通的转换 lp2bp函数可将截止频率为1rad/s的模拟低通滤波器原型转换成具有指定带宽Bw和中心频率Wn的带通滤波器。 [bt，at]=1p2bp(b，a，Wn，Bw)其中心频率为Wn，带宽为Bw。如果要求的滤波器的低端截止频率为w1，高端截止频率为w2，则可计算出Wn和Bw： Wn=sqrt(w1,w2) Bw=w2-w1 （4）从低通到带阻的转换 lp2bs函数可将截止频率为1rad/s的模拟低通滤波器原型转换成具有指定带宽Bw和中心频率Wn的带阻滤波器。 [bt，at]=lp2bs(b，a，Wn，Bw)其中心频率为Wn，带宽为Bw。Wn和Bw的计算方法与（3）相同。如果要求的滤波器的低端截止频率为w1，高端截止频率为w2，则可计算出Wn和Bw： Wn=sqrt(w1,w2) Bw=w2-w1 3、滤波器离散化 （1）冲激响应不变法 MATLAB工具箱提供了函数impinvar，它采用冲激响应不变法来实现模拟滤波器到数字滤波器的转换。 （2）双线性变换法 为了克服冲激响应不变法产生的频率混叠现象，需要使s平面与z平面建立一一对应的单值映射关系，可采用双线性变换法。MATLAB工具箱提供了函数bilinear()，实现双线性变换。 例7：设计一个巴特沃斯高通数字滤波器，它的抽样频率为Fs为100Hz,通带临界频率为40Hz, rp=0.5dB，rs=40dB，通带上限临界频率为40Hz，阻带下限临界频率为30Hz。 wp=2*pi*40; ws=2*pi*30; rp=0.5; rs=40; %滤波器指标 Fs=100; T=1/Fs;?wp = (2/T)*tan(wp/200) %预畸 ws = (2/T)*tan(ws/200) %预畸[n,wn] = buttord(wp,ws,rp,rs,‘s’) %模拟滤波器设计[z,p,k] = buttap(n);[b1,a1] = zp2tf(z,p,k); [b2,a2] = lp2hp(b1,a1,wn); %频带转换[b,a] = bilinear(b2,a2,Fs) %滤波器离散化[h,w] = freqz(b,a) %画滤波器的频响plot(w*Fs/(2*pi),abs(h)); grid; 例8：试用双线性变换法设计一个带通椭圆数字滤波器,其技术指标为wpl=l0Hz，wsl=9Hz，wp2=20Hz，ws2=21Hz，rp=0.5B，rs=50dB，抽样频率为100Hz。 wp1=2*pi*10; wp2=2*pi*20; ws1=2*pi*9; ws2=2*pi*21;Fs=100; rp=0.5; rs=50; wp1 = (2/T)*tan(wp1*T/2); wp2 = (2/T)*tan(wp2*T/2) ws1 = (2/T)*tan(ws1*T/2); ws2 = (2/T)*tan(ws2*T/2) wp=[wp1 wp2]; ws=[ws1 ws2]; [n,wn] = ellipord(wp,ws,rp,rs,s)Bw=wn(2)-wn(1); Wo=sqrt(wn(2)*wn(1));[z,p,k] = ellipap(n,rp,rs);[b1,a1] = zp2tf(z,p,k)[b2,a2] = lp2bp(b1,a1,Wo,Bw)[b,a] = bilinear(b2,a2,Fs)[H,W] = freqz(b,a);plot(W*Fs/(2*pi),abs(H)); grid; axis([0 50 0 1]); （二）利用MATLAB函数直接设计数字滤波器 1、巴特沃思滤波器 MATLAB提供了函数butter完成巴特沃思滤波器的设计。形式如下： [b,a] = butter(N,wn): 设计N阶低通滤波器，wn为它的3dB截止频率。  [b,a] = butter(N,wn,‘high’): 设计N阶高通滤波器，wn为它的3dB截止频率。 [b,a] = butter(N,wn): 当wn为具有两个元素的矢量wn=[w1 w2]时， 它设计2N 阶带通滤波器，它具有单位为π的3dB通带：w1ww2。 [b,