IIR数字滤波器设计(Ⅰ)要点

1、实验题目：iir数字滤波器设计（i ）八实验内容：数字滤波器是对数字信号实现滤波的线性时不变系统。 数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号（数 字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器 本质上是一个完成特定运算的数字计算过程，也可以理解为是一台计 算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信 号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。 时域离散系统的频域特性：,其中f9,力、工（屋中）分别是数字滤波器的输出序列和输入序列的 频域特性（或称为频谱特性），氏是数字滤波器的单位取样响应的 频谱，又称为数字滤波器的频域响应

2、。输入序列的频谱经过滤 波后某（十）h（，,因此，只要按照输入信号频谱的特点和处理信号的 目的，适当选择使得滤波后的星满足设计的要求， 这就是数字滤波器的滤波原理。数字滤波器根据其冲激响应函数的时域特性， 可分为两种，即无 限长冲激响应（iir）数字滤波器和有限长冲激响应（fir）数字滤波器。 iir数字滤波器的特征是，具有无限持续时间冲激响应，需要用递归 模型艮n（h）= z 4*（打一十-来实现，其差分方程为：fcm“=系统函数为：1二一设计iir滤波器的任务就是寻求一个物理上可实现的系统函数h(z),使其频率响应h(z)满足所希望得到的频域指标，即符合给定的 通带截止频率、阻带截止频率、

3、通带衰减系数和阻带衰减系数。设计一个数字巴特沃斯低通滤波器，设计指标如下：vp=0.2n,r=1dbw=0.3h,as=15db米样时间间隔t = 1s。三、实验要求：(1)用单位冲激响应不变变换法进行设计。(2)给出详细的滤波器设计说明书。(3)给出经过运行是正确的程序清单并加上详细的注释。(4)画出所设计滤波器的幅度特性和相位特性。四.程序与实验说明：1 .利用模拟滤波器设计iir数字滤波器方法(1)根据所给出的数字滤波器性能指标计算出相应的模拟滤波 器的设计指标。(2)根据得出的滤波器性能指标设计出相应的模拟滤波器的系统函数h(s)(3)根据得出的模拟滤波器的系统函数 h(s),经某种变

4、换得到对该模拟滤波器相应的数字仿真系统一一数字滤波器。将模拟滤波器转换成数字滤波器的实质是，用一种从s平面到z平面的映射函数将ha(s)转换成h(z)。对这种映射函数的要求是：(1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。(2)数字滤波器的频率响应模仿模拟滤波器的频响， s平面的虚 轴映射z平面的单位圆，相应的频率之间成线性关系。脉冲响应不变 法和双线性变换法都满足如上要求。以低通数字滤波器为例，将设计步骤归纳如下：(1)把数字频率转换为模拟频率：，。(2)计算巴特沃斯模拟滤波器的截止频率和阶数。(3)设计巴特沃斯模拟低通滤波器，给出参数和(此处使用了matlab 的 butt

5、ap(n)函数。(4)把模拟滤波器用单位冲激响应不变变换法转换成数字滤波器 (此处使用了 matlab的residuez函数)。(5)变直接形式为并联形式，并给出结构图。(6)画出幅度特性和相位特性。数字滤波器的设计步骤如图1所示五、程序清单藏程序function iir_imp % impulse invariance transformation% butterworth lowpass filter design %wp=0.2*pi; %数字通带频率(弧度)ws=0.3*pi; %数字阻带频率(弧度)rp=1; %!带波动(db)as=15; 蝌带衰减(db)t=1;麻样周期omega

6、p=wp*t;omegas=ws*t;%数字频率转换为模拟频率ep=sqrt(10八(rp/10)-1); % 通带波动参数ripple=sqrt(1/(1+ep*ep); % 通带波动attn=1/(10八(as/20); % 阻带衰减 % analog butterworth prototype filter calculation cs,ds=afd_butt(omegap,omegas,rp,as); %i+算巴特沃斯数字 滤波器的阶数cs和截止频率dsb,a=imp_invr(cs,ds,t); %用冲激响应不变法将模拟滤波器转化为数 字滤波器，采样频率默认1hz%c,b,a=dir

7、2par(b,a); %直接型转换成并联型figure(1); %图db,mag,pha,grd,w=freqz_m(b,a); %巴特沃斯模拟低通滤波器频率响应subplot(2,2,1);plot(w/pi,mag);title(magnitude response);xlabel(1frequency in pi units);ylabel(7h/);axis(0,1,0,1.1);set(gca,xtickmode,manual,xtick,0,02031);set(gca,ytickmode,manual,ytick,0,attn,ripple,1);grid;subplot(2,2

8、,3);plot(w/pi,db);title(magnitude in db);xlabel(2frequency in pi units);ylabel(decibels);axis(0,1,-40,5);set(gca,xtickmode,manual,xtick,0,0.2,0.3,1);set(gca,ytickmode,manual,ytick,-50,-15,-1,0);grid;set(gca,yticklabelmode,manual,yticklabels,50,15,1,0);subplot(2,2,2);plot(w/pi,pha/pi);title(phase res

9、ponse);xlabel(3frequency in pi units);ylabel(pi units);axis(0,1,-1,1);set(gca,xtickmode,manual,xtick,0,0.2,0.3,1);set(gca,ytickmode,manual,ytick,-1,0,1);grid;subplot(2,2,4);plot(w/pi,grd);title(group delay);xlabel(4frequency in pi units);ylabel(samples);axis(0,1,0,10);set(gca,xtickmode,manual,xtick,

10、0,02031);set(gca,ytickmode,manual,ytick,0:2:10);grid；六、实验结果：n =6butterworth filter order=6omegac =0.7032b =0.0000 0.0006 0.0101 0.0161 0.0041 0.0001a =1.0000-3.36355.0684-4.27592.1066-0.5706 0.0661p = -0.2588 + 0.9659i-0.2588- 0.9659i-0.7071+ 0.7071i-0.7071- 0.7071i-0.9659+ 0.2588i-0.9659- 0.2588ik

11、=10.1946-3.8164bz=-0.00000.00600.10490.06250.0021az=1.0000-2.63404.10062.4437-0.93550.2117c = 口b = 0.40190.15230.3135-0.06930.28460.0829a = 1.00000.51760.99991.00000.51760.99991.00000.51760.9999结构图如下:0.2846. -0.5176f20.0829，-0.9999lz-20.4019l -0.517c*17-1z 0.1523 -0.99991 0.3135-0.5176-： *2 - 0.0693

12、1-0.9999.f，z-2实验结果图如下:figure 1eile edit miew insert tools uesktop 工 indo 附口总国与晅口圜幅度响应magnitude response1 0.89130.177000d20.311 frequency in pi unitsmagnitude in db模值00.20.312frequency in pi units5s-aqobp七、参考文献:helpb ophase response1s 0a00.20.3相位响应j3frequency in pi unitsgroup delay10o 8 6 4 2 0 ss-de

13、em群延迟3 20. o.材：丛玉良：数字信号处理及其 matla映现电子工业出版社主要参考书：程佩青：数字信号处理教程清华大学出版社周 辉：数字信号处理基础及其 matla篁现北京希望电子出版社郭仁剑：matlab7.澈字信号处理 人民邮电出社八、心得体会：1 .数字信号处理以信号与系统、工程数学为基础，要求 学生掌握时域离散信号和系统的基本理论、基本分析方法以及 ffk数字滤波器、谱分析等数字信号处理技术。2 .数字信号处理是一门理论与实践联系紧密的课程，通 过这一周时间的课程设计，使我掌握了数字信号处理技术，提高 了分析问题和解决问题的能力，并通过设计培养了创新意识。利用matla由勺强大运算功能，基于matlab的信号处理工具箱(signal processing toolbox)的数字滤波器设计法可以快速有效地设计由软件组成的常规数字滤波器，设计方便、快捷，大大减 轻了工作量。在设计过程中可以对比滤波器特性。随时更改参数， 以达到滤波器设计的最优化。利用matla股计数字滤波器在数字 通信系统和计算机领域信号处理中，有着