《基于MATLAB的数字带通滤波器课程设计报告》由会员分享，可在线阅读，更多相关《基于MATLAB的数字带通滤波器课程设计报告（7页珍藏版）》请在金锄头文库上搜索。
1、本文格式为Word版，下载可任意编辑基于MATLAB的数字带通滤波器课程设计报告 西安文理学院机械电子工程系 课程设计报告 专业班级 08级电子信息工程1班 题 目 基于MATLAB的数字带通滤波器 学 号 学生姓名 指导教师 2022 年 12 月 西安文理学院机械电子工程系 1 课程设计任务书 学生姓名 _专业班级 _ 学 号 _ 指导教师 _ 职 称 副教授 教研室 电子信息工程 课 程 数 字 信 号 处 理 题目 基于MATLAB的数字带通滤波器设计 任务与要求 设计任务： 要求设计一个IIR带通滤波器，其中通带的中心频率为?po?0.5?，通带的截止频率?p1?0.4?，?p2?0
2、.6?，通带最大衰减?p?3dB；阻带最小衰减?s?15dB，阻带截止频率?s1?0.3?，?s2?0.7?。 设计要求： 1. 根据设计任务要求给出实现方案及实现过程。 2. 给出所实现的滤波器幅频特性及相频特性曲线并加以分析。 3. 论文要求思路明显，布局合理，语言流畅，书写格式符合要求。 开头日期 2022.12.19 完成日期 2022.12.30 2022年 12月 18 日 2 一、设计任务 设计一数字带通滤波器，用IIR来实现，其主要技术指标： 通带边缘频率：wp1=0.4?，wp2=0.6? 通带最大衰减：Ap=3dB 阻带边缘频率：ws1=0.3?，ws2=0.7? 阻带最小
3、衰减：As=15dB 设计总体要求：用MATLAB语言编程举行设计，给出IIR数字滤波器的参数，给出幅度和相位响应曲线，对IIR实现形式和特点等方面举行议论。 二、设计方法 IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR滤波器的设计可以采取在模拟滤波器设计的根基上进一步变换的方法。对比常用的原型滤波器有巴特沃什滤波器（Butterworth）、切比雪夫滤波器（Chebyshev）、椭圆滤波器（Ellipse）和贝塞尔滤波器（Bessel）等。他们有各自的特点，巴特沃什滤波器具有单调下降的幅频特性；切比雪夫滤波器的幅频特性在通带和阻带里有波动，可以提高选择性；贝塞尔滤波器通
4、带内有较好的线性相位特性；椭圆滤波器的选择性最好。本设计IIR数字滤波器采用巴特沃什滤波器3。 设计巴特沃什数字滤波器时，首先应根据参数要求设计出相应的模拟滤波器，其步骤如下： （1）由模拟滤波器的设计指标wp，ws，Ap，As和式（1）确定滤波器阶数N。 100.1Ap?1lg(0.1As)10?1N?2lg(wp)ws （1） （2）由式（2）确定wc。 3 wp(100.1Ap?1)12N?wc?ws(100.1As?1)1 （2） 2N（3）若N是奇数，那么有式（3）和式（4）确定滤波器的系统函数；若N是偶数，那么式（3）和式（5）确定滤波器的系统函数。然后把设计好的模拟Butterw
5、orth滤波器变换成数字滤波器，常采用的方法有脉冲响应不变法和双线性变换法4。 H(s)?HLO(swc) （3） 1(N?1)/21 （4） ?2HLO(s)?(s?1)k?1s?2(sin?k)s?1HLO(s)?N/21 （5） 2k?1s?2(sin?k)s?1巴特沃什滤波器是根据幅频特性在通频内具有最平坦特性而定义的滤波器。一维巴特沃什滤波器的平方幅频特性函数为 H(j?)?211?c2N （6） 可以看出，滤波器的幅频特性随着滤波器阶次N的增加而变得越来越好。在截止频率c处的函数值始终为的处境下，在通带内更多的频带 21区的值接近1；在阻带内函数更急速的趋近0。巴特沃什滤波器的主要
6、特征： （1）对于全体N，Ha(j?)?0?1 （2）对于全体N，Ha(j?)?c? 2（3）Ha(j?)是的单调下降函数。 （4）Ha(j?)随着阶次N的增大而更加接近于梦想滤波器。 最常用于设计IIR数字带通滤波器的方法基于原型模拟传输函数的双线性变换。模拟传输函数通常是：巴特沃什、切比雪夫1型、切比雪夫2 4 22221型和椭圆传输函数。巴特沃什低通传输函数在dc处有一个最大平坦的幅度响应，且随着频率增加幅度响应单调地裁减。本设计IIR数字带通滤波器就采用巴特沃什型。 巴特沃什型数字带通滤波器的设计，MATLAB调用形式如下： N,wn=buttord(wp,ws,Ap,As) num,
7、den=butter(N,wn) 在buttord中，调用参数wp，ws，是数字带通滤波器的归一化通带和阻带截止频率，Ap是通带最大衰减，As是阻带最小衰减。返回参数N是巴特沃什滤波器的阶数，wn是3dB截止频率，函数butter获得数字（den）和（num）是巴特沃什滤波器系统函数H（z）的分母多项式和分子多项式。 三、数字带通滤波器的概括设计及实现 3.1用MATLAB设计的IIR数字带通滤波器的概括程序如下： ws=0.3 0.7; %ws为归一化阻带边界频率 wp=0.4 0.6; %wp为归一化通带边界频率 Ap=3; As=15; %Ap单位为dB的通带波形 %As单位为dB的最小阻带衰减 n1,wn1=buttord(wp,ws,Ap,As);%估计数字滤波器的阶数 %wn1为归一化截止频率 %设计数字滤波器 %计算单位为dB的gain函数 %单位圆上将函数转化为256等份 w=0:pi/255:pi; %计算相位 %频率计算 %计算幅度 num,den=butter(n1,wn1); h=freqz(num,den,w); g=20*log10(abs(h); pha=angle(h); %相角计算 5 %画带通滤波器的幅度曲线图 7
