《数字butterworth滤波器的设计.doc》由会员分享，可在线阅读，更多相关《数字butterworth滤波器的设计.doc（15页珍藏版）》请在三一办公上搜索。
1、摘要巴特沃斯低通滤波器的特点是具有通带内最大平坦的幅度特性，而且随着频率的升高而单调地下降。巴特沃斯滤波器的零点全部在S等于无穷大处,在有限S平面只有极点，因而属于“全极点型”滤波器，本文利用matlab软件设计了低通巴特沃斯模拟滤波器。分别用冲激响应不变法和双线性变换法将模拟滤波器系统函数转换成数字滤波器系统函数，并且比较了两者的优缺点。同时，利用滤波器的设计方法设计一个10阶的Butterworth滤波器。关键词：巴特沃斯滤波器，matlab软件，冲激响应不变法，双线性变换法目录1. MATLAB概述.12. 设计.2 2.1题目一.22.1.1数字滤波器的工作原理.2 2.1.2低通巴特
2、沃斯模拟滤波器设计程序.22.1.3 matlab仿真结果.32.1.4巴特沃斯低通滤波器的低通特性分析.3 2.2题目二.42.2.1冲激响应不变法原理.42.2.2双线性变换法原理.52.2.3matlab设计程序.62.2.4matlab 仿真结果.62.2.5冲激响应不变法的频率混叠失真和优缺点分析.72.2.6双线性变换法存在的非线性频率失真和优缺点分析.8 2.3题目三.82.3.1设计原理分析.82.3.2matlab设计程序.92.3.3 matlab仿真结果.93.设计分析.104.设计体会.115.参考文献.121.MATLAB概述MATLAB软件由美国Math Works
3、公司于1984年推出，经过不断的发展和完善，如今己成为覆盖多个学科的国际公认的最优秀的数值计算仿真软件。MATLAB具备强大的数值计算能力，许多复杂的计算问题只需短短几行代码就可在MATLAB中实现。作为一个跨平台的软件，MATLAB已推出Unix、Windows、Linux和Mac等十多种操作系统下的版本，大大方便了在不同操作系统平台下的研究工作。MATLAB软件具有很强的开放性和适应性。在保持内核不变的情况下，MATLAB可以针对不同的应用学科推出相应的工具箱(toolbox)，目前己经推出了图象处理工具箱、信号处理工具箱、小波工具箱、神经网络工具箱以及通信工具箱等多个学科的专用工具箱，极
4、大地方便了不同学科的研究工作。国内已有越来越多的科研和技术人员认识到MATLAB的强大作用，并在不同的领域内使用MATLAB来快速实现科研构想和提高工作效率。MATLAB提供了20类图像处理函数,涵盖了图像处理的包括近期研究成果在内的几乎所有的技术方法,是学习和研究图像处理的人员难得的宝贵资料和加工工具箱。这些函数按其功能可分为:图像显示;图像文件I/O;图像算术运算;几何变换;图像登记;像素值与统计;图像分析;图像增强;线性滤波;线性二元滤波设计;图像去模糊;图像变换;邻域与块处理;灰度与二值图像的形态学运算;结构元素创建与处理;基于边缘的处理;色彩映射表操作;色彩空间变换;图像类型与类型转
5、换。2.设计2.1题目一：低通巴特沃斯模拟滤波器设计。设计一个低通巴特沃斯模拟滤波器：指标如下：通带截止频率：=3400HZ, 通带最大衰减：=3dB 阻带截至频率：=4000HZ，阻带最小衰减：=40dB2.1.1数字滤波器的工作原理数字滤波器是具有一定传输特性的数字信号处理装置。它的输入和输出均为离散的数字信号，借助数字器件或一定的数值计算方法，对输入信号进行处理，改变输入信号的波形或频谱，达到保留信号中有用成分去除无用成分的目的。如果加上A/D、D/A转换，则可以用于处理模拟信号。2.1.2低通巴特沃斯模拟滤波器设计程序根据设计要求写出如下的matlab实验程序： Wp=2*pi*340
6、0; %通带截止角频率 Ws=2*pi*4000; %阻带截止角频率 Rp=3; %通带最大衰减 Rs=40; %阻带最小衰减 n,Wn = buttord(Wp,Ws,Rp,Rs,s); %求巴特沃斯阶数和3db截止角频率 b,a = butter(n,Wn,s); %求传递函数 z,p,k = butter(n,Wn,s); %求零极点及增益 w=linspace(1,5000)*2*pi; H =freqs(b,a,w); %根据系数向量计算返回模拟滤波器的复频域响应 magH=abs(H); %频率响应的幅度 phaH=unwrap(angle(H); %频率响应的相位(平滑处理) p
7、lot(w/(2*pi),magH); %频率响应的幅度的曲线图 title(巴特沃斯低通滤波器的幅频特性); xlabel(频率/Hz); ylabel(幅度/db) 2.1.3 matlab仿真结果2.1.4巴特沃斯低通滤波器的低通特性分析 巴特沃斯低通滤波器的特点是具有通带内最大平坦的幅度特性，而且随着频率的升高而单调地下降。它的平方幅度响应为其 中 ， n为整数，称为滤波器的阶数，n值越大通带和阻带的近似性越好，过渡带也越陡 。为低通滤波器的截止频率。该滤波器具有一些特殊的性质 ： 对所有的 n ， 都有 当 =0时 ，=1 ； 对所有的 n ， 都有 当 =时 ，=，即在处有3dB的
8、衰减； 是的单调递减函数 ，即不会出现幅度响应的起伏； 当 n时 ， 巴特沃斯滤波器趋向于理想的低通滤波器 ； 在 =0处平方幅度响应的各级导数均存在且等于0， 因此在该点上取得最大值， 且具有最大平坦特性。2.2题目二：模拟低通转换为数字低通滤波器已知一模拟滤波器的系数函数为 分别用冲激响应不变法和双线性变换法将转换成数字滤波器系统函数，并图示和 的幅度相应曲线。分别取采样频率1000Hz和500Hz，分析冲激响应不变法中存在的频率混叠失真和双线性变换法存在的非线性频率失真等。2.2.1冲激响应不变法原理 冲激响应不变法是使数字滤波器的单位冲激响应序列模仿模拟滤波器的单位冲激响应，将模拟滤波
9、器的单位冲激响应加以等间隔抽样，使正好等于的抽样值，即满足：式中：T为抽样周期。 冲激不变法把稳定的转换为稳定的。由此方法可得到一阶系统的最基本的转换关系为： 2.2.2 双线性变换法原理 双线性变换法是从频域出发，使数字滤波器的频率响应与模拟滤波器的频率响应相似的一种变换法。直接使数字滤波器的频率响应，逼近模拟滤波器的频率响应,进而求得H(z)。 先将s平面压缩成s1平面上一个宽度为2/T的水平带状区域，然后通过将这个带状区域映射到整个z平面，实现s平面到z平面的单值映射。引入正切变换将s平面的虚轴变换到s1平面虚轴/T之间 改写成 延拓到整个s平面和s1平面，令，则得 再将s1平面映射到z
10、平面，应用，从而实现了s平面到z平面的单值映射。 或 （简单的代数关系） 2.2.3matlab设计程序clear; %清除变量close all %关闭所有窗口b=1000;a=1,1000; %模拟滤波器分子分母的系数w=0:1000*2*pi; %定义频带宽度hf,w=freqs(b,a,w); %把频域转化到复频域挑选w个频率点来计算频率响应subplot(2,3,1) %分割窗口画图plot(w/2/pi,abs(hf); %画滤波器的幅频特性title(模拟滤波器的幅频特性) %给这个图加个标题grid on %打开网格fs0=1000,500; %定义两个变化频率for m=1:
11、2 fs=fs0(m) d,c=impinvar(b,a,fs) %冲激响应不变法求滤波器系数f,e=bilinear(b,a,fs) %双线性变化法求滤波器系数wd=0:512*pi/512; %频率归一化hw1=freqz(d,c,wd); %求冲激响应不变法的数字滤波器的频率分量hw2=freqz(f,e,wd); %求双线性变化法的数字滤波器的频率分量subplot(2,3,2); %分割窗口画图plot(wd/pi,abs(hw1)/abs(hw1(1); %画滤波器的幅频特性hold on %保持图形不清除grid on %打开网格title(冲激响应不变法); %给这个图加个标题
12、subplot(2,3,3) %分割窗口画图plot(wd/pi,abs(hw2)/abs(hw2(1) %画滤波器的幅频特性hold on; %保持图形不清除title(双线性变化法); %给这个图加个标题end2.2.4matlab 仿真结果VVV2.2.5冲激响应不变法的频率混叠失真和优缺点分析数字滤波器的频率响应与模拟滤波器的频率响应间的关系为： 上式表明，数字滤波器的频率响应是模拟滤波器的周期延拓，根据奈奎斯特抽样定理，只有当模拟滤波器的频率响应是严格限带的，且带限于折叠频率以内时，才能使数字滤波器的频率响应在折叠频率以内，重现模拟滤波器的频率响应而不产生混叠失真。但是，任何一个实际
13、的模拟滤波器响应都不是严格限带的，变换后都会产生周期延拓分量的频谱交叠，即产生频率响应的混叠失真，因此模拟滤波器的频率响应在折叠频率以上衰减越大、越快，变换后频率响应混叠失真就越小。优点：数字滤波器的h(n)能完全模仿模拟滤波器的单位抽样响应，即时域逼近良好。 模拟频率和数字频率之间呈线性关系，频率之间不存在失真，所以一个线性相位的模拟滤波器可以映射成一个线性相位的数字滤波器。缺点：对时域的采样会造成频域的“混叠效应”，故有可能使所设计数字滤波器的频率 响应与原来模拟滤波器的频率响应相差很大。 不适用于设计高通和带阻滤波器。只适用于限带的低通或带通滤波器。2.2.6双线性变换法存在的非线性频率
14、失真和优缺点分析 根据仿真波形图可知，在零频率附近，的频率变换关系接近线性关系外，当增加时，变换关系就是非线性的了，也就是说，与之间存在严重的非线性关系。这就是双线性变换法存在的非线性频率失真。优点：避免了频率响应的混叠失真现象。 在特定数字滤波器和特定模拟滤波器处，频率响应是严格相等的，它可以较准确地控制截止频率的位置。 它是一种简单的代数关系，设计十分方便。缺点：除了零频率附近，与之间严重非线性，即线性相位模拟滤波器变为非线性相位数字滤波器。 要求模拟滤波器的幅频响应为分段常数型，不然会产生畸变。 对于分段常数的滤波器，经双线性变换后，仍得到幅频特性为分段常数的滤波器.但是各个分段边缘的临
15、界频率点产生畸变，这种频率的畸变，可通过频率预畸变来补偿。2.3题目三：设计一个10阶的Butterworth滤波器，通带为100250Hz，采样频率为1000Hz，绘出滤波器的单位脉冲响应。2.3.1设计原理分析 由题目可知，要求设计一个10阶的带通的巴特沃斯滤波器，理论上所用的方法可以是先设计出模拟低通滤波器，再转变成数字带通滤波器。然而，这里我们可以直接利用butter函数直接设计IIR数字带通滤波器。用于设计高通、带通或带阻滤波器，并由参数ftype确定滤波器的形式。当ftype（2）=band (默认值) 时，为带通滤波器，此时Wn=wl wh，则返回的b，a 所构成的滤波器是阶次为
16、2N的带通滤波器，带通范围为：。2.3.2 matlab 设计程序n=10;Wn=100 250/1000;b,a=butter(n,Wn);y,t=impz(b,a,101); %求系统的冲激响应stem(t,y) %绘制二维离散数据的图2.3.3 matlab 仿真结果3.设计分析数字滤波器分为IIR数字滤波器和FIR数字滤波器。IIR滤波器和FIR滤波器的设计方法是不同的。IIR滤波器设计方法有两类，经常采用的第一类设计方法是借助于模拟滤波器的设计方法来进行。这一类方法相对容易一些，因为模拟滤波器设计方法已经很成熟，它不仅有完整的设计公式，还有完善的图表可供查阅，并且还有一些典型的滤波器
17、类型可供设计者使用。另一类设计方法是直接在频域或者时域中进行设计，由于要解联立的方程，设计时则需要计算机做一些辅助设计。而FIR滤波器不能采用先设计模拟滤波器然后在转换为数字滤波器的方法，经常使用的设计方法有窗函数法和频率取样法，还有一种比较有效的方法是切比雪夫等波纹逼近法，需要通过计算机辅助设计来完成。 本次课设要设计的是数字Butterworth滤波器，而设计的内容分为三步：先设计一个低通巴特沃斯模拟滤波器，然后将模拟低通转换为数字低通滤波器，最后完成设计一个10阶的Butterworth带通滤波器。通过这次课程设计就是要我们掌握IIR滤波器第一类设计方法的步骤。 4.设计体会通过此次的课
18、程设计，我对数字通信理论有了更进一步的理解，并且学到了很多关于滤波器设计，及其所涉及到的一些基本函数的用法，还知道了各种模拟滤波器的特性和设计滤波器的一些基本方法与步骤。当然，在设计的过程中，不可避免的遇到了很多问题，刚开始无从下手，因为课本上只是简单的介绍了巴特沃斯滤波器的一些基本特性，以及一些其他类型的滤波器的设计步骤，无法把所学到的理论知识与实际相结合，通过查阅资料与同学的讨论，逐步明确设计方法。但是因为在理论上，将低通转换成高通，一般是变换将高通频率特征转换成低通原型频率特征；而在软件设计中，是对其传递函数进行修改，即变换其z域的表达式，设计初期一直不知道如何将其联系起来。后来发现，其
19、实变换传递函数，也就是变换频率特征，是将变换后的频率代入原低通模型，而后得到高通模型的。 总而言之，通过这次课程设计，无论是对所学课本知识的运用还是对MATLAB软件的操作应用，我都有了很大程度的提高。对数字滤波器的设计流程有了大致的了解，掌握了一些设计滤波器的基本方法，提高了理论用于实践的能力，掌握了更多专业相关的使用知识与技能。5.参考文献1陈怀琛MATLAB在电子信息课程中的应用（第二版）电子工业出版社20062 刘泉数字信号处理原理与实现（第二版）电子工业出版社20083 陈贵明用MATLAB语言处理数字信号与数字图像科学出版社20004 王家文MATLAB7.0图形图像处理国防工业出版社20065 苏金明MATLAB图形图像电子工业出版社20056 朱铭锆DSP 应用系统设计M 电子工业出版社2002.
