双线性变换法IIR数字滤波器设计..docVIP

数字信号处理实验报告 实验名称： 双线性变换法IIR数字滤波器设计 实验时间 2013 年 12 月 19 日 学号： 201111103207 姓名： 张锦龙 成绩： 评语： 一、实验目的 1、掌握用双线性变换法设计低通IIR数字滤波器的基本原理和算法。 2、掌握用双线性变换法设计高通和带通IIR数字滤波器的基本原理和算法。 3、进一步了解数字滤波器和模拟滤波器的频率响应特性二、实验原理与计算方法 １、双线性变换法设计IIR低通数字滤波器的基本原理和算法 双线性变换法设计数字滤波器，采用了二次映射的方法，就是先将整个s平面压缩到s1平面的一个的横形条带范围内，然后再将这个条带映射到z平面上，就能建立s平面到z平面的一一对应关系。对于低通数字滤波器，映射关系为 　　 　 (1) 其中T为抽样周期。 用双线性变换法设计低通IIR数字滤波器的基本步骤，首先根据设计要求确定相应的模拟滤波器的传递函数，再应用(1)式得数字滤波器的传递函数 　　　　　　　　　 　　　　　　　　　　　　　　 (2)通常可以给定的参数为：低通数字滤波器通带边界频率、阻带边界频率和对应的通带衰减函数、阻带衰减函数。s1平面中的模拟角频率与数字角频率的关系为线性关系，在计算模拟滤波器的阶数N、极点si和传递函数之前，应作预畸变处理 　　　 　　　　 (3) 　　模拟滤波器的阶数N、极点si和传递函数的计算方法与冲激响应不变法相同，可以采用Butterworth逼近或Chebyshev逼近。 ２、双线性变换法设计IIR高通带通数字滤波器的基本原理和算法由于双线性变换法获得的数字滤波器频率响应特性中不会出现混叠现象，因此可以适用于高通、带通和带阻滤波器的设计。IIR数字滤波器的设计通常要借助于模拟低通滤波器的设计，给定的对数字滤波器(DF)的技术要求转换为一个低通模拟滤波器(AF)的技术要求，根据这种要求用某种逼近设计出原型的低通模拟滤波器(LP AF)，计算出模拟滤波器的阶数N、极点si和传递函数，再按照双线性变换的变换关系，将模拟滤波器的传递函数转换为数字滤波器的传递函数。表8-1中列出将给定的对数字滤波器(DF)的技术要求直接转换为对一个低通模拟滤波器(AF)的技术要求的频率预校正关系表8-1 双线性变换和频率预校正的计算公式 变换类型 变换关系 频率预校正 备　　　　注 低通变换 其中，T为抽样周期，f为模拟频率 其中，分别为数字带通滤波器通带的上下边界频率数字带滤波器的上下边界频率。 高通变换 带通变换 例：数字高通滤波器的设计 首先将给定的数字高通滤波器的技术指标根据公式转换为模拟低通滤波器的技术指标，利用cheb1ord(Wp,Ws,ap,as,s)函数求出chebyshev模拟低通滤波器的阶数N，再利用cheb1ap(N,ap)函数求出模拟低通滤波器系统函数的零极点，zp2tf(z,p,k)函数将零极点转换为系统函数系数；然后利用lp2hp由模拟低通滤波器的系统函数得到模拟带通滤波器的系统函数，bilinear函数则用于实现双线性变换法由模拟带通滤波器系统函数计算数字数字带通滤波器系统函数的系数。 （2）数字频带变换法 首先将给定的对数字滤波器(DF)的技术要求转换为一个低通模拟滤波器(AF)的技术要求，原型的低通模拟滤波器(LP AF)，滤波器数字滤波器。[bhp,ahp]=zmapping(blp,alp,Nz,Dz)用来实现从数字低通滤波器得到数字高通滤波器的有理函数。 %数字滤波器技术指标 wp=0.2*pi;ws=0.3*pi;Rp=1;As=15; %对应的模拟滤波器技术指标 T=1;Fs=1/T;Wp=(2/T)*tan(wp/2); =(2/T)*tan(ws/2); [cs,cd]=afd_chb1(Wp,Ws,Rp,As); %Chebyshev模拟滤波器 [blp,alp]=bilinear(cs,cd,Fs) %双线性变换 wphp=0.6*pi; %数字高通滤波器截止频率 %低通－高通频带变换 alpha=-(cos((wplp+wphp)/2))/(cos((wplp+wphp)/2)) Nz=-[alpha,1];Dz=[1,alpha]; [bhp,ahp]=zmapping(blp,alp,Nz,Dz) %数字高通滤波器的系统函数系数 （3）IIR数字滤波器的设计 可利用M