1、吉林建筑大学电气与电子信息工程学院数字信号处理课程设计报告设计题目：IIR数字滤波器的设计专业班级：信工111学生姓名：马超学号：*指导教师：高晓红王超设计时间：2014.01.062014.01.10教师评语：成绩 评阅教师 日期 目录一、设计目的 1二、设计内容 1三、设计原理 13.1 数字低通滤波器的设计原理 13.2 变换方法的原理 2四、设计步骤 8五、数字低通滤波器 MATLAB 编程及幅频特性曲线 105.1 MATLAB语言编程 105.2 幅频特性曲线 12六、总结 13七、参考文献 13一、设计目的课程设计是理论学习的延伸，是掌握所学知识的一种重要手段，对于贯彻理论联系实
2、际、提高学习质量、塑造自身能力等于有特殊作用。本次课程设计一方面通过MATLAB仿真设计内容，使我们加深对理论知识的理解，同时增强其逻辑思维能力，另一方面对课堂所学理论知识作一个总结和补充。二、设计内容已知通带截止频率fp=0.2kHz,通带最大衰减p=1dB,阻带截止频率fs=0.3kHz，阻带最小衰减s=25dB，T=1ms，按照以上技术要求，用脉冲响应不变法和双线性变换法设计巴特沃斯数字低通滤波器，并观察所设计数字滤波器的幅频特性曲线。三、设计原理3.1数字低通滤波器的设计原理滤波器是自动控制、信号处理和通信领域的重要组成部分,广泛地应用于各种系统中。数字滤波器是指输入、输出均为数字信号
3、，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。与模拟滤波器相比具有很多突出的优点，例如它可以满足滤波器对幅度和相位特性的严格要求，可以避免模拟滤波器所无法克服的电压漂移和噪声问题。设计数字滤波器，首先要按一定规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标。根据转换后的技术指标设计模拟低通滤波器G(s)，再按一定的规则将G(s)转换成H(z)。若是高通、带通或带阻数字滤波器则将它们的技术指标先转化为低通模拟滤波器的技术指标，然后设计出低通G(s)，再将G(s)转换为所需的H(z)将系统函数Ha(s)从s平面转换到z平面的方法有多种，但工程上常用的
4、是脉冲响应不变法和双线性变换法。这两种方法有着各自的优缺点，在设计巴特沃斯数字低通滤波器时需要根据相应要求或想要达成的效果从中选择。3.2变换方法的原理脉冲响应不变法原理：Ha(s) 拉氏逆变换 ha(t) 等间隔采样 ha(nT)=h(n) Z变换 H(z) 脉冲响应不变法是一种将模拟滤波器转化为数字滤波器的基本方法。它利用模拟滤波器理论设计数字滤波器，也就是使数字滤波器能模仿模拟滤波器的特性，这种模仿可从不同的角度出发。脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n) 模仿模拟滤波器的冲击响应ha(t), 使h(n)正好等于ha(t)的采样值，即T为采样周期。
5、利用模拟滤波器理论设计数字滤波器，也就是使数字滤波器能模仿模拟滤波器的特性，这种模仿可从不同的角度出发。脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n) 模仿模拟滤波器的冲击响应ha(t), 使h(n)正好等于ha(t)的采样值，即T为采样周期。如以Ha(s)及H（z）分别表示ha(t)的拉氏变换及h(n)的z变换，即则根据采样序列z变换与模拟信号拉氏变换的关系，得:上式表明，采用脉冲响应不变法将模拟滤波器变换为数字滤波器时，它所完成的S平面到Z平面的变换，正是以前讨论的拉氏变换到Z变换的标准变换关系，即首先对Ha(s)作周期延拓，然后再经过z=est的映射关系映
6、射到Z平面上。z=est的映射关系表明，S平面上每一条宽为2/T的横带部分，都将重叠地映射到Z平面的整个全部平面上。每一横带的左半部分映射到Z平面单位圆以内，每一横带的右半部分映射到Z平面单位圆以外，j轴映射在单位圆上，但j轴上的每一段2/T都对应于绕单位圆一周，如图1所示。图1脉冲响应不变法的映射关系应当指出，z=est的映射关系反映的是Ha(s)的周期延拓与H（z）的关系，而不是Ha(s)本身与H（z）的关系，因此，使用脉冲响应不变法时，从Ha(s)到H(z)并没有一个由S平面到Z平面的简单代数映射关系，即没有一个s=f(z)的代数关系式。另外，数字滤波器的频响也不是简单的重现模拟滤波器的
7、频响，而是模拟滤波器频响的周期延拓，周期为s=2/T=2fs，即如果模拟滤波器的频响带限于折叠频率s/2以内，即|/T这时数字滤波器的频响才能不失真地重现模拟滤波器的频响|但任何一个实际的模拟滤波器，其频响都不可能是真正带限的，因此不可避免地存在频谱的交叠，即混淆，这时，数字滤波器的频响将不同于原模拟滤波器的频响而带有一定的失真。模拟滤波器频响在折叠频率以上衰减越大，失真则越小，这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。虽然脉冲响应不变法能保证S平面与Z平面的极点位置有一一对应的代数关系，但这并不是说整个S平面与Z平面就存在这种一一对应的关系，特别是数字滤波器的零点位置与S平
8、面上的零点就没有一一对应关系，而是随着Ha(s)的极点si与系数Ai的不同而不同。H(ej) 是Ha(j)的周期延拓（周期为fs），因Ha(j)并不是带限，即在超过fs频率部分并不为0，所以就产生了混迭。当为低通或带通滤波器时，fs越大，则Ha(j)的下一周期相隔越远，混迭也就越小。当为带阻或高通滤波器时，Ha(j)在超过fs/2频率部分全为通带，这样就不满足抽样定理，发生了完全的混迭，所以脉冲响应不变法不能设计带阻或高通滤波器。双线性变换法原理：脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整
9、个频率轴上的频率范围压缩到-/T/T之间，再用转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-/T/T一条横带里；第二步再通过标准变换关系将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图2所示。图2双线性变换的映射关系为了将S平面的整个虚轴j压缩到S1平面j1轴上的-/T到/T段上，可以通过以下的正切变换实现（1）式中,T仍是采样间隔。当1由-/T经过0变化到/T时，由-经过0变化到+，也即映射了整个j轴。将式（1）写成将此关系解析延拓到整个S平面和S1平面，令j=s，j1=s1，则得再将S
10、1平面通过以下标准变换关系映射到Z平面z=es1T，从而得到S平面和Z平面的单值映射关系为：（2）（3）式（2）与式（3）是S平面与Z平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双线性变换。式（1）与式（2）的双线性变换符合映射变换应满足的两点要求。首先,把z=ej，可得（4）即S平面的虚轴映射到Z平面的单位圆。其次，将s=+j代入式（4），得因此由此看出，当0时，|z|0时，|z|1。也就是说，S平面的左半平面映射到Z平面的单位圆内，S平面的右半平面映射到Z平面的单位圆外，S平面的虚轴映射到Z平面的单位圆上。因此，稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的
11、。双线性变换法与脉冲响应不变法相比，其主要的优点是避免了频率响应的混叠现象。这是因为S平面与Z平面是单值的一一对应关系。S平面整个j轴单值地对应于Z平面单位圆一周，即频率轴是单值变换关系。这个关系如式（4）所示，重写如下：上式表明，S平面上与Z平面的成非线性的正切关系，如图3所示。由图3看出，在零频率附近，模拟角频率与数字频率之间的变换关系接近于线性关系；但当进一步增加时，增长得越来越慢，最后当时，终止在折叠频率=处，因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象，从而消除了频率混叠现象。但是双线性变换的这个特点是靠频率的严重非线性关系而得到的，如前边式（4）及图3所
12、示。图3双线性变换法的频率变换关系由于这种频率之间的非线性变换关系，就产生了新的问题。首先，一个线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器，不再保持原有的线性相位了；还有一点，就是这种非线性关系要求模拟滤波器的幅频响应必须是分段常数型的，即某一频率段的幅频响应近似等于某一常数（这正是一般典型的低通、高通、带通、带阻型滤波器的响应特性），不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变，如图4所示。图4双线性变换法幅度和相位特性的非线性映射对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但是各个分段边缘的临界频率点产生了畸变，这种频率
13、的畸变，可以通过频率的预畸来加以校正。也就是将临界模拟频率事先加以畸变，然后经变换后正好映射到所需要的数字频率上。由于双线性变换法获得的数字滤波器频率响应特性中不会出现混叠现象，因此可以适用于高通、带通和带阻滤波器的设计。四、设计步骤脉冲响应不变法设计IIR数字低通滤波器步骤：（1）将数字滤波器设计指标转换为相应的模拟滤波器指标。rad/sdBrad/sdB（2）设计相应的模拟滤波器，得到模拟系统函数Ha(s)。根据单调下降要求，选择巴特沃斯滤波器。取N=9（3）按照下面公式，将模拟滤波器系统函数Ha(s)转换成数字滤波器系统函数H(z)。，按照上边的式子求解计算相当复杂，我们可以利用MATL
14、AB软件来简化计算过程，调用MATLAB信号处理工具箱函数进行设计。双线性法设计IIR数字低通滤波器步骤：（1）列出数字低通技术指标。raddBraddB（2）这里T=0.001s，预畸变校正计算相应模拟低通的技术指标为。dBdB（3）设计巴特沃斯低通模拟滤波器。据算阶数N。求出c,，保证阻带技术指标满足要求，通带指标有富余。根据N的值查表得到归一化方程Ga(p)，将p=s/c带入Ga(p)，去归一化得到实际的Ha(s)。（4）用双线性变换法将Ha(s)转换成数字滤波器H(z),即按照上式计算的过程也是非常复杂，我们可以通过MATLAB软件使计算得到简化。根据所给出的技术指标并加以转换，转换之
15、后，调用MATLAB软件，然后应用MATLAB软件设计MATLAB程序来简化计算过程，而且应用MATLAB软件可以得到相应的滤波器图像，从而更直观的观察所设计的滤波器的相关性质，程序中采用双线性变换法，调用MATLAB工具箱函数buttord和butter直接设计数字滤波器五、数字低通滤波器MATLAB编程及幅频特性曲线5.1 MATLAB语言编程（1）模拟滤波器设计程序：Fs=1000;%采样频率为周期倒数Wp=400*pi;Ws=600*pi;%设置归一化通带和阻带截止频率Ap=1;As=25;%设置通带最大和最小衰减N,Wc=buttord(Wp,Ws,Ap,As,s);%调用butto
16、rd函数确定巴特沃斯滤波器阶数B,A=butter(N,Wc,s);%调用butter函数设计巴特沃斯滤波器W=(0:pi:1000*pi);%指定一段频率值hs=freqs(B,A,W);%计算模拟滤波器的幅频响应plot(W/pi/2,abs(hs)/abs(hs(1);%绘出巴特沃斯模拟滤波器的幅频特性曲线gridon;title(巴特沃斯模拟滤波器);xlabel(频率/Hz);ylabel(归一化幅值);(2) 用脉冲响应不变法设计数字滤波器的程序：Fs=1000;%采样频率为周期倒数Wp=400*pi;Ws=600*pi;%设置归一化通带和阻带截止频率Ap=1;As=25;%设置通
17、带最大和最小衰减N,Wc=buttord(Wp,Ws,Ap,As,s);%调用buttord函数确定巴特沃斯滤波器阶数B,A=butter(N,Wc,s);%调用butter函数设计巴特沃斯滤波器W=(0:pi:1000*pi);%指定一段频率值D,C=impinvar(B,A,Fs);%调用脉冲响应不变法W=(0:0.001*pi:pi);Hz=freqz(D,C,W);%返回频率响应plot(W/pi,abs(Hz)/abs(Hz(1);%绘出巴特沃斯数字低通滤波器的幅频特性曲线gridon;title(巴特沃斯数字滤波器（脉冲响应不变法）);xlabel(频率/Hz);ylabel(归一
18、化幅值);(3)用双线性变换法设计数字滤波器的程序：Fs=1000;%采样频率为周期倒数Wp=400*pi;Ws=600*pi;%设置归一化通带和阻带截止频率Ap=1;As=25;%设置通带最大和最小衰减N,Wc=buttord(Wp,Ws,Ap,As,s);%调用buttord函数确定巴特沃斯滤波器阶数B,A=butter(N,Wc,s);%调用butter函数设计巴特沃斯滤波器W=(0:pi:1000*pi);%指定一段频率值D,C=bilinear(B,A,Fs);%调用双线性变换法W=(0:0.001*pi:pi);Hz=freqz(D,C,W);%返回频率响应plot(W/pi,ab
19、s(Hz)/abs(Hz(1);%绘出巴特沃斯数字低通滤波器的幅频特性曲线gridon;title(巴特沃斯数字滤波器（双线性变换法）);xlabel(频率/Hz);ylabel(归一化幅值);5.2 幅频特性曲线六、总结低通滤波器的应用十分广泛，通过MATLAB，可以很容易地设计出巴特沃斯数字低通滤波器，而且通过MATLAB设计滤波器更方便、快捷,节省了很多实践,提高了编程效率,并且参数的修改也十分方便. 还可以进一步进行优化设计。通过这几天的设计过程，让我感触最深的就是MATLAB功能的强大性与掌握编程各种函数和语句的重要性。当然首先要了解所要编程运行的对象的原理。在课程设计的过程中，我深
20、深的感受到我们所掌握的知识还不够，需要学习的东西太多了，真的是学无止境。学习的过程很辛苦，但是我们也学到了很多知识，最后感谢老师给我们这次学习的机会，让我发现自己的不足，让我知道以后应该更投入的学习知识，提高自己的学习能力。七、参考文献1 数字信号处理（第2版），高西全 丁玉美 编著，西安电子科技大学出版社，2001年2MATLAB信号处理详解，陈亚勇主编.，人民邮电出版社，2001年3数字信号处理教程MATLAB释疑与实现陈怀琛编著，电子工业出版社，2004年4信号处理原理及应用，谢平王娜林洪彪编著，机械工业出版社，2008年5 数字信号处理（第3版），高西全丁玉美编著，清华大学出版社，2010年
