基于MATLAB低通巴特沃斯滤波器仿真.docVIP

基于MATLAB低通巴特沃斯滤波器仿真摘要：本文介绍了低通巴特沃斯滤波器的原理，分析了低通巴特沃斯滤波器的特性 ，然后用MATLAB对低通巴特沃斯滤波器进行了仿真。 关键词：低通巴特沃斯滤波器 MATLAB 仿真 中图分类号：TN911.7 文献标识码：A 文章编号：1007-9416（2013）02-0124-01 1 巴特沃斯滤波器原理 巴特沃斯滤波器是电子滤波器的一种，这种滤波器由英国工程师斯替芬·巴特沃斯（Stephen Butterworth）在1930年提出，它的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。可以构成低通、高通、带通和带阻四种组态，是目前最为流行的一类数字滤波器 ，经过离散化可以作为数字巴特沃思滤波器，较模拟滤波器具有精度高、稳定、灵活、不要求阻抗匹配等众多优点，因而在语音、图像、通信、雷达等众多领域得到了广泛应用，是一种具有最大平坦幅度响应的低通滤波器。 2 MATLAB简介 MATLAB是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案。MATLAB主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。在信号处理中低通巴特沃斯滤波器常用的函数包括： 2.1 buttord函数 [N，wn]=buttord（wp，ws，rp，rs，’s’） 用于计算巴特沃斯数字滤波器的阶数N和3dB截止频率wn。其中，调用参数wp，ws分别为数字滤波器的通带、阻带截止频率的归一化值，要求：0≤wp≤1，0≤ws≤1（1表示数字频率pi）。当ws≤wp时，为高通滤波器；当wp和ws为二元矢量时，为带通或带阻滤波器，这时wn也是二元向量。 rp，rs分别为通带最大衰减和组带最小衰减（dB）。 2.2 butter函数 [B，A]=butter（N，wn，ftype’） 计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式的系数向量B、A。其中，调用参数N和wn分别为巴特沃斯数字滤波器的阶数和3dB截止频率的归一化值（关于pi归一化），一般是调用buttord（1）格式计算N和wn。系数B、A是按照z-1的升幂排列，ftype为滤波器的类型。 N，wn为butter函数的调用参数。 2.3 buttap函数 由[Z，P，K] = buttap（N）函数可设计出N阶巴特沃斯低通滤波器的零、极点。 3 仿真指标 设计一个低通巴特沃斯模拟滤波器，绘制滤波器的的幅频响应及零极点图。指标如下： 通带截止频率：WP=1000HZ，通带最大衰减：RP=3dB 阻带截止频率：Ws=2000HZ，阻带最小衰减：Rs=40 dB 4 程序 subplot（1，2，1）；%两个窗口，幅频图在第一个窗口 wp=1000；ws=2000；rp=3；rs=40；%指标参数设置 [N，wn]=buttord（wp，ws，rp，rs，’s’）%计算巴特沃斯低通滤波器的阶数和3dB截止频率 [B，A]=butter（N，wn，’s’）；%代入N和wn设计低通巴特沃斯模拟滤波器 [Z，P，K]=buttap（N）；%计算滤波器的零、极点 [h，w]=freqs（B，A，1024）；%计算1024点模拟滤波器频率响应h，和对应的频率点w plot（w，20*log10（abs（h）/abs（h（1）））） %画频率响应幅度图 grid； xlabel（’频率Hz’）；ylabel（’幅度（dB）’）； %给x轴和y轴加标注 title（’巴特沃斯幅频响应’） axis（[0，3000， -40，3]）； line（[0，2000]，[-3，-3]）； line（[1000，1000]，[-40，3]）； %绘制巴特沃斯滤波器的极点 subplot（1，2，2） %在第二个窗口画极点图 p=P’；q=Z’； x=max（abs（[p，q]））； x=x+0.1；y=x； axis（[-x，x，-y，y]）； axis（’square’） plot（[-x，x]，[0，0]）；hold on plot（[0，0]，[-y，y]）；hold