用冲激响应不变法设计IIR数字滤波器.ppt

用冲激响应不变法设计IIR数字滤波器 第41组 钱坤 许守都 设计目的 1、掌握构成一个频率响应与给定的滤波特性相接近的模拟滤波器的设计原 ?2、掌握用冲激响应不变法设计IIR数字滤波器的基本原理和算法。? 3、了解数字滤波器和模拟滤波器的频率响应特性，掌握相应的计算方法，分析用冲激响应不变法获得的数字滤波器频率响应特性中出现的混叠现象。 设计原理 １、冲激响应不变法设计IIR数字滤波器的基本原理和算法 采用冲激响应不变法设计数字滤波器，就是使其单位样值响应与相应的模拟滤波器的冲激响应在抽样点处的量值相等，即 (1) 其中T为抽样周期。 因此用冲激响应不变法设计IIR数字滤波器的基本步骤，就是首先根据设计要求确定相应的模拟滤波器的传递函数 Ha(s) ，经Laplace反变换求出冲激响应ha(t)后，对它进行抽样得到的ha(nt)等于数字滤波器的单位样值响应h(n)，再经z变换所得H(z)就是数字滤波器的传递函数。 实验原理 如果模拟滤波器的传递函数Ha(s)的N个极点si都是单极点，则可以将Ha(s)写成部分分式展开的形式 (2) 那么，经Laplace反变换求出的模拟滤波器的冲激响应为 相对应的数字滤波器的单位样值响应为 对上式作z变换，得 (3) 由上面的推导可见，只要模拟滤波器的传递函数Ha(s)的N个极点si都是单极点，当已经求出各个极点值si和部分分式的系数后，则可以从模拟滤波器的传递函数的表达式（2）直接得到数字滤波器的传递函数H(z)的表达式（3）。 实验内容 仿真结果： 2、用脉冲响应不变法设计IIR高通数字滤波器 数字高通的技术指标为：Wp=0.4πrad??????Ap=2dB? T=2s??????????????????Ws=?0.2rad????????????????As=15dB 仿真结果： 3、用脉冲响应不变法设计IIR带通数字滤波器 仿真结果： 4、用脉冲响应不变法设计IIR带阻数字滤波器 仿真结果： 小结： 经观察滤波器幅频特性图得，巴特沃斯滤波器的特点是通带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。?在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。? 低通滤波器容许低频信号通过,?但减弱(或减少)频率高於截止频率的信号的通过。?高通滤波器容许高频信号通过,?但减弱(或减少)频率低于於截止频率的信号的通过。?带通滤波器容许一定频率范围信号通过,?但减弱(或减少)频率低于於下限截止频率和高于上限截止频率的信号的通过。? 带阻滤波器减弱(或减少)一定频率范围信号,?但容许频率低于於下限截止频率和高于上限截止频率的信号的通过。? 对于巴特沃斯滤波器低通和带通效果较好，而高通和帯阻效果较差。 1、用脉冲响应不变法设计IIR低通数字滤波器 ?数字低通的技术指标为：?Wp=0.2πrad?????Ap=1dB? T=2s????Ws=?0.3rad????????????As=15dB ?程序为：? T=2;?????????? ?fs=1/T;???????????????%采样频率为采样周期倒数????????Wp=0.2.*pi;? ?Ws=0.3.*pi;????????%设计归一化通带阻带截止频率????????Ap=1;? ?As=15;????????%设置通带最大最小衰减???? [N,Wc]=buttord(Wp,Ws,Ap,As,s);??????%调用butter函数确定巴特沃斯滤波器阶数?? [B,A]=butter(N,Wc,s);???????????%调用butter函数设计巴特沃斯滤波器?? W=linspace(0,pi,400*pi);??????????????%指定一段频率值? [D,C]=impinvar(B,A,fs);????????????%调用脉冲不变法¨???? Hz=freqz(D,C,W);??????????????????%·返回频率值|??????????????????????? plot(W/pi,abs(Hz)/abs(Hz(1)));????????%绘出巴特沃斯数字低通滤波器的幅频特性曲线?grid?on;?? title(巴特沃斯低通数字滤波器);?xlabel(Frequency/Hz);?ylabel(Magnitude);? 程序为： T=2; %采样周期 fs=1/T; %采样频率 Wp= 0.4.*pi; Ws= 0.2.*pi; %设置归一化通带和阻带截止平率 Ap= 2; As=15; %设置通