数字信号处理实验冲激响应不变法IIR数字滤波器设计方案.docVIP

个人收集整理 仅供参考学习 个人收集整理 仅供参考学习 PAGE / NUMPAGES 个人收集整理 仅供参考学习 数字信号处理实验报告 实验名称：冲激响应不变法IIR数字滤波器设计 学号：姓名： 评语： 成绩： 实验目地 1、掌握构成一个频率响应与给定地滤波特性相接近地模拟滤波器地设计原理. 2、掌握用冲激响应不变法设计IIR数字滤波器地基本原理和算法. 3、了解数字滤波器和模拟滤波器地频率响应特性,掌握相应地计算方法,分析用冲激响应不变法获得地数字滤波器频率响应特性中出现地混叠现象.b5E2RGbCAP 二、实验原理与计算方法 1、冲激响应不变法设计IIR数字滤波器地基本原理和算法 采用冲激响应不变法设计数字滤波器,就是使其单位样值响应与相应地模拟滤波器地冲激响应在抽样点处地量值相等,即p1EanqFDPw 　　　　 (1) 其中T为抽样周期. 因此用冲激响应不变法设计IIR数字滤波器地基本步骤,就是首先根据设计要求确定相应地模拟滤波器地传递函数,经Laplace反变换求出冲激响应后,对它进行抽样得到地等于数字滤波器地单位样值响应,再经z变换所得就是数字滤波器地传递函数.DXDiTa9E3d 如果模拟滤波器地传递函数地N个极点都是单极点,则可以将写成部分分式展开地形式 　　　　　　　　　　　 (2) 那么,经Laplace反变换求出地模拟滤波器地冲激响应为 相对应地数字滤波器地单位样值响应为 对上式作z变换,得 (3) 由上面地推导可见,只要模拟滤波器地传递函数地N个极点都是单极点,当已经求出各个极点值和部分分式地系数后,则可以从模拟滤波器地传递函数地表达式（2）直接得到数字滤波器地传递函数地表达式（3）.RTCrpUDGiT 2、Butterworth和Chebyshev模拟滤波器地设计方法 （1）Butterworh原型 MATLAB提供了函数[z,p,k]=buttap(N)用来设计N阶归一化地Butterworth模拟低通滤波器；函数[b,a]=u_buttap(N,Omegac)给出未归一化地Butterworth模拟低通滤波器原型；再利用函数[C,B,A]=sdir2cas(b,a)得到级联形式地N阶Butterworth模拟低通滤波器原型.N=3;OmegaC=0.5;5PCzVD7HxA [b,a]=u_buttap(N,Omegac); [C,B,A]=sdir2cas(b,a) （2）按给定技术指标设计Butterworth模拟低通滤波器 函数[b,a]=afd_butt(Wp,Ws,Rp,As)用来实现按给定技术指标设计Butterworth模拟低通滤波器；函数[db,mag,pha,w]=freqs_m(b,a,0.5*pi)绘制模拟滤波器频响特性,得出衰减值,幅频特性,相频特性和自变量w；函数[H,w]=freqs(b,a,0.5*pi)绘制模拟滤波器频响特性,给出复函数H,再由abs(H)和angle(H)求出幅频和相频特性；[ha,x,t]=impulse(b,a)求出模拟滤波器冲激响应.jLBHrnAILg Wp=0.2*pi;Ws=0.3*pi;Rp=7;As=16; [b,a]=afd_butt(Wp,Ws,Rp,As) [db,mag,pha,w]=freqs_m(b,a,0.5*pi) [ha,x,t]=impulse(b,a) （3）数字滤波器设计－Butterworth 函数[b,a]=imp_incr(c,d,T)用来实现冲激响应不变地映射,c、d分别为模拟滤波器系统函数Ha(s)地分子、分母地系数；T为抽样周期；xHAQX74J0X %数字滤波器技术指标 wp=0.2*pi;ws=0.3*pi;Rp=1;As=15; %对应地模拟滤波器技术指标 T=1;Wp=wp/T;Ws=ws/T; [cs,cd]=afd_butt(Wp,Ws,Rp,As); [b,a]=imp_invr(cs,ds,T) %并联形式地系统函数系数 [C,B,A]=dir2par(b,a) 函数[H,w]=freqz(b,a)绘制数字滤波器频响特性,给出复函数H,再由abs(H)和angle(H)求出幅频和相频特性.LDAYtRyKfE （4）Chebyshev原型 MATLAB提供了函数[z,p,k]=cheblap(N,Rp)用来设计N阶通带波动为Rp地归一化Chebyshev模拟低通滤波器；函数[b,a]=u_chb1ap(N,Omegac)给出未归一化地Chebyshev模拟低通滤波器原型.Zzz6ZB2Ltk （5）按给定技术指标设计Chebysheb模拟低通滤波器 函数[b,a]=afd_chb1(Wp,Ws,R