巴特沃斯带阻滤波器仿真.docVIP

摘 要 巴特沃斯滤波器是电子滤波器的一种。巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大。 1 理论简单介绍 1.1 MATLAB概述 MATLAB 是一个可视化的计算程序，被广泛地应用在科学运算领域里。它具有功能强大、使用简单等特点，内容包括：数值计算、符号计算、数据拟合、图形图像处理、系统模拟和仿真分析等功能。此外，用MATLAB还可以进行动画设计、有限元分析等。 1.2 MATLAB系统组成 MATLAB系统包含下列五部分： 1）开发环境：这是一组帮助你使用MATLAB的函数和文件的工具和设备。这些工具大部分是图形用户界面。它包括MATLAB桌面和命令窗口，命令历史，和用于查看帮助的浏览器，工作空间，文件和查找路径。 2）MATLAB数学函数库：这里汇集了大量计算的算法，范围从初等函数如：求和，正弦，余弦和复数的算术运算，到复杂的高等函数如：矩阵求逆，矩阵特征值，贝塞尔(Bessel)函数和快速傅立叶变换等。 3）MATLAB语言：这是一种高水平的矩阵/数组语言，含有控制流语句，函数，数据结构，输入/输出，和面向对象编程特征。它允许“小型编程”以迅速创立快速抛弃型程序，以及“大型编程”以创立完整的大型复杂应用程序。 4）句柄制图：这是MATLAB制图系统。它包括高级别的二维、三维数据可视化，图像处理，动画，以及表现图形的命令。它还包括低级别的命令，这使你不但能在MATLAB的应用中建立完整的图形用户界面，而且还能完全定制图形的外观。 5）MATLAB应用程序界面(API)：这是使你编写与MATLAB相合的C或Fortran程序的程序库。它包括从MATLAB中调用程序（动态链接）， 调用MATLAB为计算引擎，和读写MAT-文件的设备。 2 巴特沃斯带阻IIR数字滤波器设计 2.1 数字带阻IIR滤波器设计 IIR数字滤波器在很多领域中有着广阔的应用。与FIR数字滤波器相比,它可以用较低的阶数获得高选择性,所用存储单元少,经济而效率高,在相同门级规模和相同时钟速度下可以提供更好的带外衰减特性。下面介绍一种设计实现IIR数字滤波器的方法：冲击响应不变法。 2.2 设计思路—基于冲激响应不变法的IIR数字滤波器设计 冲激响应不变法的设计原理是利用数字滤波器的单位抽样响应序列H(z)来逼近模拟滤波器的冲激响应g(t)。 按照冲激响应不变法的原理，通过模拟滤波器的系统传递函数G(s)，可以直接求得数字滤波器的系统函数，其转换步骤如: (1) 利用 (可由关系式推导出)，将，转换成，,而，不变； (2) 求解低通模拟滤波器的传递函数； (3) 将模拟滤波器的传递函数转换为数字滤波器的传递函数。 尽管通过冲激响应不变法求取数字滤波器的系统传递函数比较方便，并具有良好的时域逼近特性，但若不是带限的，或是抽样频率不高，那么在中将发生混叠失真，数字滤波器的频率响应不能重现模拟滤波器的频率响应。只有当模拟滤波器的频率响应在超过折叠频率后的衰减很大时，混叠失真才很小，此时采样脉冲响应不变法设计的数字滤波器才能满足设计的要求，这是冲激响应不变法的一个严重的缺点。 2.3 设计指标 抽样频率为,在衰减处的边带频率是，，在衰减处频率是和 2.4 用MTALAB算法设计巴特沃斯带阻IIR数字滤波器 MATLAB仿真程序如下： Fs=10000; wp=[1500 4000]*2/Fs;ws=[2000 3000]*2/Fs; %阻带边界频率 Rp=-2; Rs=-13; %带通衰减和阻带衰减 Nn=512; %显示滤波器的频率特性的数据长度 [N,Wn]=buttord(wp,ws,Rp,Rs); %求的数字滤波器的最小阶数和归一化截止频率 [b,a]=butter(N,Wn,stop); %按最小阶数，截止频率设计数字滤波器 [H,f]=freqz(b,a,Nn,Fs); %绘制滤波器的频率特性 figure(1) subplot(2,1,1),plot(f,20*log10(abs(H))) xlabel(频率/Hz);ylabel(振幅/^o); title(幅度响应曲线（付兴娜）); grid on; subplot(2,1,2),plot(f,180/pi*unwrap(angle(H))) xlabel(频率/Hz);ylabel(振幅/^o); title(相位响应曲线（付兴娜）); grid on; 2.5仿真波形 仿真出的幅度响应及相位响应曲线如下图2.1所示： 图2.1：幅度响应及相位响应曲线 3模拟IIR带阻滤波器的设计 3.1 模拟滤波器的分类 典型模拟滤波器L