《基于MATLAB设计巴特沃斯低通滤波器》由会员分享，可在线阅读，更多相关《基于MATLAB设计巴特沃斯低通滤波器（6页珍藏版）》请在金锄头文库上搜索。
1、本文格式为Word版，下载可任意编辑基于MATLAB设计巴特沃斯低通滤波器 基于MATLAB设计巴特沃斯低通滤波器 摘 要： 首先分析了巴特沃斯低通滤波器的特性。然后用MATLAB的信号处理工具箱供给的函 数设计了巴特沃斯低通滤波器，使得巴特沃斯滤波器的设计变得更加简朴、快捷、直观。 巴特沃斯(Butterworth)滤波器是一种具有最大平坦幅度响应的低通滤波器，它在通信领域里已有广应用，在电测中也具有广泛的用途，可以作检测信号的滤波器。MATLAB语言是一种面向科学与工程计算的语言。它编程效率高，测试程序手段丰富，扩展才能强，内涵丰富。它的信号处理工具箱(Signal Processing
2、Toolbox)供给了设计巴特沃斯滤波器的函数，本文充分利用这些函数，举行了巴特沃斯滤波器的程序设计，并将其作为函数文件保存，可便当地举行调用。 1. 巴特沃斯低通滤波器的特性 巴特沃斯低通滤波器的平方幅度响应为： H(j?)2?1?(1?c )2n 其中，n为滤波器的阶数，c为低通滤波器的截止频率。该滤波器具有一些特殊的性质： 对全体的n，都有当 =0时，|H(j0)|2 =1； 对全体的n，都有当=c时，|H(jc)|2 =0.5 ，即在c处有3dB的衰减； |H(j)|2 是的单调递减函数，即不会展现幅度响应的起伏； 当n+时，巴特沃斯滤波器趋向于梦想的低通滤波器； 在=0处平方幅度响应
3、的各级导数均存在且等于0，因此|H(j)|2 在该点上取得最大值，且具有最大平坦特性。 图l表示了2阶、4阶、8阶巴特沃斯低通滤波器的幅频特性。可见阶数n越高，其幅频特性越好，低频检测信号保真度越高。 巴特沃斯与贝塞尔(Besse1)、切比雪夫(Cheby.shev)滤波器的特性差异如图2所示。 从图2可以看出，巴特沃斯滤波器在线性相位、衰减斜率和加载特性三个方面具有特性均衡的优点，因此在实际使用中，巴特沃斯滤波器已被列为首选。 2 .巴特沃斯低通滤波器的MATLAB实现 MATLAB的信号处理工具箱供给了有关巴特沃斯滤波器的函数buttap、buttord、butter。 21 用butta
4、p函数设计归一化的巴特沃斯低通滤波器系数 z，P，k=bttap(n)可设计出n阶巴特沃斯低通滤波器原型，其传递函数为： H(s)?k(s?p(1)(s?p(2)?(s?p(n) 因此实际上Z为空阵。上述零极点形式可以化为: H(s)?b0s?bn?1snnn?1?b1s?b0 其中 b0?0，巴特沃斯滤波器的归一化结果，如表1所示。 22 设计巴特沃斯滤波器的程序实现 buttord函数可在给定滤波器性能的处境下，选择巴特沃斯滤波器的阶数n和截止频率c， 从而可利用butter函数设计巴特沃斯滤波器的传递函数。 n,c=buttord(p,s,Rp,Rs,s)可得到得志性能的模拟巴特沃斯滤波
5、器的最小阶数n及截止频率c， 其中p为通带的拐角频率，s为阻带的拐角频率， p和s的单位均为rad/s；Rs为通带区的最大波动系数，Rp为Rs阻带区的最小衰减系数，Rp和Rs的单位都为dB。 b，a=butter(n，c ,s )可设计截止频率为c的n阶低通模拟巴特沃斯滤波器为: H(s)?B(s)A(s)?b(1)s?b(2)ss?a(2)snnn?1?b(n?1)n?1?a(n?1) 利用buttord函数、buuer函数编制设计巴特沃斯低通滤波器的MATLAB函数文件butterdesignm，其清单如下： FunctionNbutt，Wcbutt，bbutt，abutt=butterd
6、esign(Wp，Rp，Ws，As) Nbutt，Wcbutt=buttord(Wp，Ws，Rp，As，S)； bbutt，abutt=butter(Nbutt，Wcbutt，S)； hbutt，Wbutt=freqs(bbutt，abutt)； subplot(2，1，1)；plot(WbuR，abs(hbutt)； subplot(2，1，2)；plot(Wbutt，angle(hbutt)； 为了使程序更通用，采用以下命令从键盘上输入巴特沃斯滤波器的性能指标。 Wp=input(Please input passband comer frequencyin rads Wp=) Rp=in
7、put(Please input passband ripple in dB Rp= ) Ws=input(Please input stopband comer frequency in rads Ws=) As=input(Please input stopband attenuation in dB As=) 调用butterdesign函数就可设计出得志性能指标的巴特沃斯滤波器，其主程序文件名为buttermainm。下面举例说明巴特沃斯滤波器的计算机辅佐设计过程。 设计一模拟巴特沃斯低通滤波器，其通带的拐角频率p为02，通带区的最大波动系数Rp为7dB，阻带的拐角频率s为03，阻带区
8、的最小衰减系数Rs为16dB。 在MATLAB命令窗口中输入buttermain，根据屏幕上的提示分别输入02*pi，7，03*pi，16，就可得到得志性能指标的巴特沃斯低通滤波器的最小阶数n=3，截止频率c=05122rads，其传递函数为: H(s)?0.1344s?1.0244s?0.5247s?0.1344 32同时得到图3所示的巴特沃斯低通滤波器曲线。 3. 终止语 利用本文的butterdesign函数就可以轻而易举地设计出得志性能指标要求的巴特沃斯滤波器的最小阶数n及截止频率c ，并绘制出巴特沃斯滤波器的频率特性曲线，使得巴特沃斯滤波器的设计变得更加简朴、快捷、直观。 参 考 文 献 1 张殿龙，王福文巴特沃斯低通滤波器在电动机测试中的应用J电测与仪表，2022. 2 王永初Buttworth滤波器在过程操纵中的应用J工业仪表与自动打扮置，1994，6 3 陈思巴特沃斯低通滤波器的简化快速设计J信阳师范学院学报(自然科学版)，1997. 4 楼顺天，李博菡基于MATLAB的系统分析与设计 信号处理M.西安电子科技大学出版社,1998 6
