基于MATLAB的数字带通滤波器.docx

西安文理学院机械电子工程系课程设计报告 专业班级 题 目 基于MATLAB的数字带通滤波器学 号 学生姓名指导教师 2011 年 12 月 西安文理学院机械电子工程系 1 PAGE PAGE 10 课程设计任务书 学生姓名 专业班级 学 号 指导教师 职 称 副教授 教研室 电子信息工程 课 程 数 字 信 号 处 理 题目 基于 MATLAB 的数字带通滤波器设计 任务与要求 设计任务： 要求设计一个 IIR 带通滤波器，其中通带的中心频率为? po  ? 0.5? ，通 带的截止频率? ? 0.4? ，? ? 0.6? ，通带最大衰减? ? 3dB ；阻带最小 p1 p 2 p 衰减? ? 15dB ，阻带截止频率? ? 0.3? ， ? ? 0.7? 。 s s1 s 2 设计要求： 根据设计任务要求给出实现方案及实现过程。 给出所实现的滤波器幅频特性及相频特性曲线并加以分析。 论文要求思路清晰，结构合理，语言流畅，书写格式符合要求。 开始日期 2011.12.19 完成日期 2011.12.30 2011 年 12 月 18 日 一、设计任务 设计一数字带通滤波器，用 IIR 来实现，其主要技术指标： 通带边缘频率：wp1=0.4? ，wp2=0.6 ? 通带最大衰减：Ap=3dB 阻带边缘频率：ws1=0.3? ，ws2=0.7? 阻带最小衰减：As=15dB 设计总体要求：用 MATLAB 语言编程进行设计，给出 IIR 数字滤波器的参数，给出幅度和相位响应曲线，对 IIR 实现形式和特点等方面进行讨论。 二、设计方法 IIR 数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR 滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。比较常用的原型滤波器有巴特沃什滤波器（ Butterworth）、切比雪夫滤波器（Chebyshev）、椭圆滤波器（Ellipse）和贝塞尔滤波器（Bessel）等。他们有各自的特点，巴特沃什滤波器具有单调下降的幅频特性；切比雪夫滤波器的幅频特性在通带和阻带里有波动，可以提高选择性；贝塞尔滤波器通带内有较好的线性相位特性；椭圆滤波器的选择性最好。本设计 IIR 数字滤波器采用巴特沃什滤波器[3。] 设计巴特沃什数字滤波器时，首先应根据参数要求设计出相应的模拟滤波器，其步骤如下： 由模拟滤波器的设计指标wp，ws，Ap，As 和式（1）确定滤波器阶数 N。 lg( 10 0 .1 Ap ? 1 ) 2 lg( w p )w sN ? 10 0 2 lg( w p ) w s （1） 由式（2）确定 wc。 (100.1Ap ? 1) 12 N(100.1As ? (100.1Ap ? 1) 12 N (100.1As ? 1) 12 N w （2） 若N 是奇数，则有式（3）和式（4）确定滤波器的系统函数；若N 是偶数，则式（3）和式（5）确定滤波器的系统函数。然后把设计好的模拟 Butterworth 滤波器变换成数字滤波器，常采用的方法有脉冲响应不变法和双线性变换法[4。] H (s) ? H ( s LO ) （3） w c H (s) ? 1 ( N??1) / 2 1 （4） LO (s ? 1)  k ?1 s 2 ? 2(sin? )s ? 1 k H (s) ? N?/ 2 1 （5） LO k ?1 2(sin? )s ? 1 k s s 2 波器。一维巴特沃什滤波器的平方幅频特性函数为 H ( j?) 2 ? ?1 1? ?? c ?2 N  （6） 可以看出，滤波器的幅频特性随着滤波器阶次 N 的增加而变得越来越 好。在截止频率 Ω 1 c 处的函数值始终为 的情况下，在通带内更多的频带 2 区的值接近 1；在阻带内函数更迅速的趋近 0。巴特沃什滤波器的主要特征： 对于所有 N， H 对于所有 N， H ( j?)2 ? 1 a ?? 0 ( j?) 2 ? 1 a ?? ?c 2 H a( j?)2 是 Ω 的单调下降函数。 H a( j?)2 随着阶次 N 的增大而更加接近于理想滤波器。 最常用于设计 IIR 数字带通滤波器的方法基于原型模拟传输函数的双线性变换。模拟传输函数通常是：巴特沃什、切比雪夫 1 型、切比雪夫 2 型和椭圆传输函数。巴特沃什低通传输函数在 dc 处有一个最大平坦的幅度响应，且随着频率增加幅度响应单调地减少。本设计 IIR 数字带通滤波器就采用巴特沃什型。 巴特沃什型数字带通滤波器的设计，MATLAB 调用形式如下： [N,wn]=buttord(wp,ws,Ap,As) [num,den]=butter(N,wn) 在 buttord 中，调用参数 wp，ws