双线性巴特沃兹数字滤波器设计.docx

电气工程学院课程设计任务书课程名称： 数字信号处理课程设计 基层教学单位：仪器科学与工程系 指导教师： 学号100103020168学生姓名王大明（专业）班级 精仪二班设计题目12、 双线性变换法设计数字低通滤波器设计技术参数给定技术指标为：，采样频率。设计要求设计Butterworth低通滤波器，用双线性变换法转换成数字滤波器。(buttord,buttap,butter,bilinear)参考资料数字信号处理方面资料MATLAB方面资料周次前半周后半周应完成内容收集消化资料、学习MATLAB软件，进行相关参数计算编写仿真程序、调试指导教师签字基层教学单位主任签字说明：1、此表一式四份，系、指导教师、学生各一份，报送院教务科一份。 2、学生那份任务书要求装订到课程设计报告前面。电气工程学院 教务科摘 要数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数字处理来达到频域滤波的目的。本文是设计一个数字低通滤波器。根据滤波器的设计思想，通过双线性变换法将巴特沃斯模拟低通滤波器变换到数字低通滤波器，利用MATLAB绘制出数字低通滤波器的系统幅频函数曲线。关键词：数字滤波器；双线性变换法；巴特沃斯；MATLAB目 录1课题描述12设计原理12.1 IIR数字滤波器设计原理12.2巴特沃斯低通滤波器的原理22.3双线性变换法33设计过程64结果分析8总 结11参考文献121课题描述数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。可以设计系统的频率响应，让它满足一定的要求，从而对通过该系统的信号的某些特定的频率成分进行过滤，这就是滤波器的基本原理。如果系统是一个连续系统，则滤波器称为模拟滤波器。如果系统是一个离散系统，则滤波器称为数字滤波器。数字滤波实质上是一种运算过程，实现对信号的运算处理。输入数字信号（数字序列）通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可以理解为是一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性:其中、分别是数字滤波器的输出序列和输入序列的频域特性（或称为频谱特性）, 是数字滤波器的单位取样响应的频谱，又称为数字滤波器的频域响应。输入序列的频谱经过滤波后,因此，只要按照输入信号频谱的特点和处理信号的目的， 适当选择，使得滤波后的满足设计的要求，这就是数字滤波器的滤波原理。2设计原理2.1 IIR数字滤波器设计原理IIR数字滤波器的设计一般是利用目前已经很成熟的模拟滤波器的设计方法来进行设计，通常采用模拟滤波器原型有butterworth函数、chebyshev函数、bessel函数、椭圆滤波器函数等。IIR数字滤波器的设计步骤： 按照一定规则把给定的滤波器技术指标转换为模拟低通滤波器的技术指标； 根据模拟滤波器技术指标设计为响应的模拟低通滤波器； 跟据脉冲响应不变法和双线性不变法把模拟滤波器转换为数字滤波器； 如果要设计的滤波器是高通、带通或带阻滤波器，则首先把它们的技术指标转化为模拟低通滤波器的技术指标，设计为数字低通滤波器，最后通过频率转换的方法来得到所要的滤波器。在MATLAB中，经典法设计数字滤波器主要采用以下步骤：模拟滤波器原型buttap,cheb1ap频率变换模拟离散化bilinear,impinvarIIR数字滤波器图2.1 IIR数字滤波器设计步骤2.2巴特沃斯低通滤波器的原理巴特沃斯滤波器的特点是同频带内的频率响应曲线最为平坦，没有起伏，而在组频带则逐渐下降为零。在振幅的对数对角频率的波特图上，从某一边界见频率开始，振幅随着角频率的增加而逐渐减少，趋向于负无穷大。一阶巴特沃斯滤波器的衰减率为每倍频20分贝，二阶巴特沃斯滤波器的衰减率为每倍频12分贝，三阶的衰减率为每分贝18分贝，如此类推，巴特沃斯滤波器的振幅对角频率单调下降，并且滤波器的结束越高，在组频带振幅衰减速度越快，其他滤波器高阶的振幅对角频率图和低阶数的振幅对角频率有不同的形状。 上述函数的特点是等距离分布在半径为的圆上。因此，极点用下式表示为的表示式：为了使设计公式和图表统一，将频率归一化。巴特沃斯滤波器采用3dB截止频率归一化，归一化后的系统函数为令,称为归一化频率，称为归一化复变量，这样巴特沃斯滤波器的归一化低通原型系统函数为式中，为归一化极点，用下式表示: 2.3双线性变换法双线性变换法是将平面压缩变换到某一中介平面的一条横带里，再通过标准变换关系将此带变换到整个z平面上去，这样就使平面与平面之间建立一一对应的单值关系，消除了多值变换性。为了将平面的轴压缩到平面的轴上的到一段上，可以通过以下的正切变换来实现：这样当由经变化到时，由经过0变化到，也映射到了整个 轴。将这个关系延拓到整个平面和平面，则可以得到再将平面通过标准变换关系映射到平面，即令得到同样对z求解，得到这样的变换叫做双线性变换。为了验证这种映射具有s平面的虚轴映射到z平面单位圆上的特性，考虑 ,，得除了使s平面的虚轴映射到单位圆上之外，s平面的左半部分映射到单位圆的内部，s平面的右半部分映射到单位圆的外部。如图所示图2.2 双线性变化映射关系示意图观察式子，发现的实部为负时，因子的幅度小于1，相当于单位圆的内部。反之，当的实部为负时，该比值的幅度大于1，相当于单位圆的外部。这样就可以看出使用双线性变换可从稳定的模拟滤波器得到稳定的数字滤波器。双线性变换法还避免了使用脉冲响应不变法所遇到的混叠问题，因为它把平面的这个虚轴映射到平面的单位圆上。然而，付出的代价是在频率轴上引入了失真。因此，只有当能容忍或补偿这种失真时，使用双线性变换法设计数字滤波器的方法才是实用的。仅在零频率附近时与之间的频率变换关系接近于线性关系，所产生的数字滤波器的幅频响应相对于原模拟滤波器的幅频响应有畸变。对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但是各分段边缘的临界频率点产生了畸变，这种频率的畸变，可以通过频率的预畸变来加以校正，也就是将临界频率事先加以畸变，然后经变换后正好映射到所需要的频率上。通过的关系变换成一组模拟频率。图2.3 双线性变化法的频率关系为了克服冲击响应不变法产生的频率混叠现象，我们需要使平面与平面建立一一对应的单值关系，即求出，然后将其代入就可以求得，即3设计过程由设计技术参数所给定技术指标计算出归一化系统函数： = N=-lg0.092/lg1.568=5.306 取N=6求得3dB截止频率为，查表得归一化传输函数 ，编写MATLAB程序实现从设计截止频率为的六阶低通巴特沃斯数字滤波器。步骤一：将设计内容题所给归一化巴特沃斯低通滤波器以截止频率为进行去归一化。步骤二：用双线性变化法将低通模拟滤波器变换为低通数字滤波器设计程序如下：clear all; clc; close all fs=1; T=1; N=6;wc=0.7662; omegach=2*tan(wc/2)/T;M=1; N=1,3.8636,7.4637,9.1412,7.4637,3.8636,1h,w=freqs(M,N,512); %模拟滤波器的幅频响应subplot(2,1,1);plot(w,20*log10(abs(h); grid; axis(0,10,-90,0)xlabel(Hz);ylabel(幅度); title(归一化模拟低通滤波器);Ms,Ns=lp2lp(M,N,omegach); %对低通滤波器进行频率变换hs,ws=freqs(Ms,Ns,512); %模拟滤波器的幅频响应subplot(2,1,2);plot(ws,20*log10(abs(hs); grid;axis(0,10,-90,0)xlabel(Hz);ylabel(幅度); title(去归一化模拟低通滤波器);Mz,Nz=bilinear(Ms,Ns,1/T); %对模拟滤波器双线性变换h1,w1=freqz(Mz,Nz); %数字滤波器的幅频响应figureplot(w1/pi,20*log10(abs(h1); grid;xlabel(/);ylabel(幅度(dB); title(数字低通滤波器);axis(0,1,-160,0)运行结果如下图所示：图3.1模拟滤波器的幅频响应图形图3.2低通数字滤波器的幅频响应图形4结果分析比较脉冲响应不变法设计的低通滤波器和双线性法设计的低通滤波器进行比较：优点：是频率坐标变换是线性的，即，如果不考虑频率混叠现象，用这种方法设计的数字滤波器会很好的重现原模拟滤波器的频率特性。另一个优点是数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，时域逼近性好。缺点：是会产生频率混叠现象适合低通、带通滤波器的设计，不适合高通、带阻滤波器的设计。 双线性变换法避免了频率响应的混叠现象。可以设计高通、带阻滤波器。模拟频率与数字频率不再是线性关系，所以一个线性相位模拟滤波器经双线性变换后所得到的数字滤波器不再保持原有的线性相位了。 与脉冲响应不变法相比，双线性变换靠频率的严重非线性关系得到S平面与Z平面的单值一一对应关系，整个轴单值对应于单位圆一周，这个关系就是式 所表示的，其中和为非线性关系。如图中看到，在零频率附近，接近于线性关系，进一步增加时，增长变得缓慢，时 (终止于折叠频率处)，所以双线性变换不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象。双线性变换法的缺点：与的非线性关系，导致数字滤波器的幅频响应相对于模拟滤波器的幅频响应有畸变，(使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变)。例如，一个模拟微分器，它的幅度与频率是线性关系，但通过双线性变换后，就不可能得到数字微分器。若 则另外，一个线性相位的模拟滤波器经双线性变换后，滤波器就不再有线性相位特性。虽然双线性变换有这样的缺点，但它目前仍是使用得最普遍、最有成效的一种设计工具。这是因为大多数滤波器都具有分段常数的频响特性，如低通、高通、带通和带阻等，它们在通带内要求逼近一个衰减为零的常数特性，在阻带部分要求逼近一个衰减为的常数特性，这种特性的滤波器通过双线性变换后，虽然频率发生了非线性变化，但其幅频特性仍保持分段常数的特性。 双线性变换比脉冲响应法的设计计算更直接和简单。由于s与z之间的简单代数关系，所以从模拟传递函数可直接通过代数置换得到数字滤波器的传递函数。置换过程：频响： 这些都比脉冲响应不变法的部分分式分解便捷得多，一般，当着眼于滤波器的时域瞬态响应时，采用脉冲响应不变法较好，而其他情况下，对于IIR的设计，大多采用双线性变换。 总 结本次设计采用了双线性变换法。滤波器设计原理来比较简单，周一周二我查阅了课本和相关书籍，了解了滤波器设计的基本原理。周三周四学习了MATLAB的基本操作和常用函数，了解了其中专门为实现滤波器的函数。然后开始了基本运算和编程，因为使用了MATLAB使得很多原本复杂的运算简单了很多。但中间由于对MATLAB了解的不多，甚至连基本的函数语句都得查书才能知道，这就大大增加了我们设计的时间。有时在调用函数时出现各种错误，得不到正确结果，程序运行时也会有很多的错误。这让我深深体会到，要想正确得出一个波形，需要花很多的时间和很大精力。滤波器的设计难度高，挑战性强，内容多。通过这次设计我对原来课本