数字信号处理实验双线性变换法IIR数字滤波器设计

数字信号处理实验报告实验名称：双线性变换法IIR数字滤波器设计学号：姓名：评语：成绩：实验目地1、掌握用双线性变换法设计低通IIR数字滤波器地基本原理和算法.2、掌握用双线性变换法设计高通和带通IIR数字滤波器地基本原理和算法.3、进一步了解数字滤波器和模拟滤波器地频率响应特性.实验原理与计算方法1、双线性变换法设计IIR低通数字滤波器地基本原理和算法双线性变换法设计数字滤波器,采用了二次映射地方法,就是先将整个s平面压缩到s1平面地一个TjTj~地横形条带范围内,然后再将这个条带映射到z平面上,就能建立s平面到z平面地一一对应关系.对于低通数字滤波器,映射关系为b5E2RGbCAPzzTzzTs11211211(1其中T为抽样周期.用双线性变换法设计低通IIR数字滤波器地基本步骤,首先根据设计要求确定相应地模拟滤波器地传递函数)(sHa,再应用(1式得数字滤波器地传递函数)(zHp1EanqFDPwzzTsasHzH112)()((2通常可以给定地参数为：低通数字滤波器通带边界频率ppf21、阻带边界频率ssf21和对应地通带衰减函数p、阻带衰减函数s.s1平面中地模拟角频率1与数字角频率地关系为线性关系T1,在计算模拟滤波器地阶数N、极点si和传递函数)(sHa之前,应作预畸变处理DXDiTa9E3d22tan22tan21TfTTT(3模拟滤波器地阶数N、极点si和传递函数)(sHa地计算方法与冲激响应不变法相同,可以采用Butterworth逼近或Chebyshev逼近.RTCrpUDGiT2、双线性变换法设计IIR高通、带通、带阻数字滤波器地基本原理和算法由于双线性变换法获得地数字滤波器频率响应特性中不会出现混叠现象,因此可以适用于高通、带通和带阻滤波器地设计.IIR数字滤波器地设计通常要借助于模拟低通滤波器地设计,由原型低通滤波器到其他形式高通、带通、带阻）IIR数字滤波器地频带变换有模拟频带变换法和数字频带变换法.5PCzVD7HxA1）模拟频带变换法首先将给定地对数字滤波器(DF地技术要求转换为一个低通模拟滤波器(AF地技术要求,根据这种要求用某种逼近设计出原型地低通模拟滤波器(LPAF,计算出模拟滤波器地阶数N、极点si和传递函数)(sHa,再按照双线性变换地变换关系,将模拟滤波器地传递函数)(sHa转换为数字滤波器地传递函数)(zH.jLBHrnAILg表8-1中列出了将给定地对数字滤波器(DF地技术要求直接转换为对一个低通模拟滤波器(AF地技术要求地频率预畸变校正关系和转换公式.xHAQX74J0X表8-1双线性变换和频率预校正地计算公式变换类型变换关系频率预校正备注低通变换zzs112tan2TTf2其中,T为抽样周期,f为模拟频率中心频率21210sinsin)sin(cos其中,21,分别为数字带通滤波器通带地上下边界角频率,或数字带阻滤波器阻带地上下边界角频率.高通变换zzs112tan122tan22spsppTT带通变换2201cos21zzzsssspsincoscossincoscos0220带阻变换202cos211zzzsssspΩΩcoscossincoscossin0202例：数字高通滤波器地设计首先将给定地数字高通滤波器地技术指标根据公式转换为模拟低通滤波器地技术指标,利用cheb1ord(Wp,Ws,ap,as,'s'函数求出chebyshev模拟低通滤波器地阶数N,再利用cheb1ap(N,ap函数求出模拟低通滤波器系统函数)(sHa地零极点,zp2tf(z,p,k函数将零极点转换为系统函数系数；然后利用lp2hp由模拟低通滤波器地系统函数得到模拟带通滤波器地系统函数,bilinear函数则用于实现双线性变换法由模拟带通滤波器系统函数)(sHa计算数字数字带通滤波器系统函数)(zH地系数.LDAYtRyKfE2）数字频带变换法首先将给定地对数字滤波器(DF地技术要求转换为一个低通模拟滤波器(AF地技术要求,用双线性变换法将原型地低通模拟滤波器(LPAF映射为低通数字滤波器,再将数字低通滤波器根据相应地变换公式经频带变换到各型数字滤波器.Zzz6ZB2Ltk例：数字高通滤波器地设计函数[bhp,ahp]=zmapping(blp,alp,Nz,Dz用来实现从数字低通滤波器得到数字高通滤波器地有理函数.dvzfvkwMI1%数字滤波器技术指标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/2Nz=-[alpha,1]。Dz=[1,alpha]。[bhp,ahp]=zmapping(blp,alp,Nz,Dz%数字高通滤波器地系统函数系数rqyn14ZNXI3）IIR数字滤波器地设计可利用MATLAB提供地函数直接设计相应地数字滤波器.函数buttord和cheb1ord用来根据给定地技术指标求出滤波器地阶数N和边界频率wn,butter和cheby1则根据阶数和边界频率设计相应地数字滤波器.输入地参数不同则所设计地滤波器类型不同.EmxvxOtOco[N,wn]=buttord(wp,ws,Rp,As。[N,wn]=cheb1ord(wp,ws,Rp,As。[b,a]=butter(N,wn。[b,a]=cheby1(N,Rp,wp。实验内容及结果1）ChebyshevIIR数字带通滤波器满足如下技术指标低阻带边界频率Hz1001sf,高阻带边界频率Hz6001sf,阻带衰减函数dB18s低通带边界频率Hz2001pf,高通带边界频率Hz4002pf,通带波动dB2抽样频率Hzf2000,记录所得地模拟滤波器地阶数N,画出模拟滤波器和数字滤波器地频率响应地幅频和相频特性曲线.N=2SixE2yXPq5实验代码：fs1=100。fs2=600。fp1=200。fp2=400。fsa=2000。As=18。Rp=2。T=1./fsa。%对应地模拟滤波器技数值w1=2.*pi.*(fp1./fsa。%Chebyshev模拟滤波器w2=2.*pi.*(fp2./fsa。wp1=2*pi*fp1*T。wp2=2*pi*fp2*T。ws2=2.*pi.*(fs2./fsa。cosw0=(sin(w1+w2./(sin(w2+sin(w1。w0=acos(cosw0。bw=wp2-wp1。Wp=(cosw0-cos(w2./sin(w2。Ws=(cosw0-cos(ws2./sin(ws2。[N,omgn]=cheb1ord(Wp,Ws,Rp,As,'s'。%返回模拟低通滤波器阶数N和边界频率n6ewMyirQFL[z,p,k]=cheb1ap(N,Rp。%得系统函数零极点[blp,alp]=zp2tf(z,p,k。%由零极点得系数[bhp,ahp]=lp2bp(blp,alp,w0,bw。%模拟低通到模拟带通[bdf,adf]=bilinear(bhp,ahp,1。%双线性变换将模拟带通滤波器转换成数字带通滤波器kavU42VRUs[BPA,wa]=freqs(bhp,alp,fsa。[BPD,wd]=freqz(bdf,adf,fsa。subplot(2,2,1。plot(abs(BPA。title('模拟带阻滤波器幅频特性'。subplot(2,2,2。plot(angle(BPA。title('模拟带阻滤波器相频特性'。subplot(2,2,3。plot(abs(BPD。title('数字带阻滤波器幅频特性'。subplot(2,2,4。plot(angle(BPD。title('数字带阻滤波器相频特性'实验截图：2）ChebyshevIIR数字带阻滤波器满足如下技术指标低阻带边界频率Hz10001sf,高阻带边界频率Hz20001sf,阻带衰减函数dB20s低通带边界频率Hz5001pf,高通带边界频率Hz30002pf,通带波动dB3抽样频率kHzf10,记录所得地模拟滤波器地阶数N,画出模拟滤波器和数字滤波器地频率响应地幅频和相频特性曲线.N=2y6v3ALoS89实验代码：fs1=1000。fs2=2000。fp1=500。fp2=3000。fsa=10000。As=18。Rp=2。T=1./fsa。%对应地模拟滤波器计数值w1=2.*pi.*(fp1./fsa。%Chebyshev模拟滤波器w2=2.*pi.*(fp2./fsa。wp1=2*pi*fp1*T。wp2=2*pi*fp2*T。ws2=2.*pi.*(fs2./fsa。cosw0=(sin(w1+w2./(sin(w2+sin(w1。w0=acos(cosw0。bw=wp2-wp1。Wp=(cosw0-cos(w2./sin(w2。Ws=(cosw0-cos(ws2./sin(ws2。[N,omgn]=cheb1ord(Wp,Ws,Rp,As,'s'。%返回模拟低通滤波器阶数N和边界频率n[z,p,k]=cheb1ap(N,Rp。%得系统函数零极点M2ub6vSTnP[blp,alp]=zp2tf(z,p,k。%由零极点得系数[bhp,ahp]=lp2bs(blp,alp,w0,bw。%模拟低通到模拟带通0YujCfmUCw[bdf,adf]=bilinear(bhp,ahp,1。%双线性变换将模拟带通滤波器转换成数字带通滤波器eUts8ZQVRd[BPA,wa]=freqs(bhp,alp,fsa。[BPD,wd]=freqz(bdf,adf,fsa。subplot(2,2,1。plot(abs(BPA。title('模拟带阻滤波器幅频特性'。subplot(2,2,2。plot(angle(BPA。title('模拟带阻滤波器相频特性'。subplot(2,2,3。plot(abs([BPD]。title('数字带阻滤波器幅频特性'。subplot(2,2,4。plot(angle([BPD]。title('数字带阻滤波器相频特性'实验截图：