用巴特莱特窗函数法设计数字FIR带通滤波器dsp课程设计
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
课程设计课程设计名称:数字信号处理课程设计
专业班级:电信
学生姓名:
学号:
指导教师:乔丽红
课程设计时间: 6.16-6.20
电子信息工程专业课程设计任务书
说明:本表由指导教师填写,由教研室主任审核后下达给选题学生,装订在设计(论文)首页
一需求分析和设计内容
数字信号处理是把许多经典的理论体系作为自己的理论基础,同时又使自己成为一系列新兴学科的理论基础。现如今随着电子设备工作频率范围的不断扩大,电磁干扰也越来越严重,接收机接收到的信号也越来越复杂。为了得到所需要频率的信号,就需要对接收到的信号进行过滤,从而得到所需频率段的信号,这就是滤波器的工作原理。对于传统的滤波器而言,如果滤波器的输入,输出都是离散时间信号,则该滤波器的冲激响应也必然是离散的,这样的滤波器定义为数字滤波器。它通过对采样数据信号进行数学运算来达到频域滤波的目的.
滤波器在功能上可分为四类,即低通(LP)、高通(HP)、带通(BP)、带阻(BS)滤波器等,每种又有模拟滤波器(AF)和数字滤波器(DF)两种形式。对数字滤波器,从实现方法上,具有有限长冲激响应的数字滤波器被称为FIR滤波器,具有无限长冲激响
应的数字滤波器被称为IIR 滤波器。
FIR 数字滤波器的主要优点有:一、具有严格的线性相位特性;二、不存在稳定性问
题;三、可利用DFT 来实现。这些优点使FIR 数字滤波器得到了广泛应用。窗函数法是
一种设计FIR 数字滤波器的基本方法,但它不是最佳设计方法,在满足同样设计指标的
情况下,用这种方法设计出的滤波器的阶数通常偏大。在窗函数法的基础上,以所定义
的逼近误差最小为准则来进行优化设计的算法,由于其中的逼近误差可根据不同的设计
要求进行定义,故此算法适应性强,它既可用于设计选频型滤波器,又适用于非选频型
滤波器的设计。常用的窗函数有矩形窗函数、巴特莱特窗函数、三角窗函数、汉宁(Hann )
窗函数、海明(Hamming )窗函数、布莱克曼(Blackman )窗函数、凯塞(Kaiser )窗函
数等。本设计通过MATLAB 软件对FIR 型滤波器进行理论上的实现,利用巴特莱特窗函数
设计数字FIR 带通滤波器。FIR 系统不像IIR 系统那样易取得较好的通带和阻带衰减特性,
要取得较好的衰减特性,一般要求H (z )阶次要高,也即M 要大。FIR 系统有自己突出的
优点:系统总是稳定的;易实现线性相位;允许设计多通带(或多阻带)滤波器,后两项都
是IIR 系统不易实现的。FIR 数字滤波器的设计方法有多种,如窗函数设计法、频率采样
法和Chebyshev 逼近法等。随着Matlab 软件尤其是Matlab 的信号处理工作箱的不断完
善,不仅数字滤波器的计算机辅助设计有了可能,而且还可以使设计达到最优化。本实验
的数字滤波器的MATLAB 实现是指调用MATLAB 信号处理工具箱函数filter 对给定的输入
信号x(n)进行滤波,得到滤波后的输出信号y(n )。
用巴特莱特窗函数法设计一个数字FIR 带通滤波器,要求通带边界频率为400Hz ,
500Hz ,阻带边界频率为350Hz ,550Hz ,通带最大衰减1dB ,阻带最小衰减40dB ,抽样频
率为2000Hz ,用MATLAB 画出幅频特性,画出并分析滤波器传输函数的零极点;
信号)2sin()2sin()()()(2121t f t f t x t x t x ππ+=+=经过该滤波器,其中=1f 450Hz ,
=2f 600Hz ,滤波器的输出)(t y 是什么?用Matlab 验证你的结论并给出
)(),(),(),(21t y t x t x t x 的图形。
二 设计原理及设计思路
1.设计FIR 数字滤波器的基本方法:
FIR 数字滤波器的系统函数无分母,为11
00()()N N i
n i i i H Z b z h n z ----====∑∑,系统频率响应可写成:1
0()()N jw
jwn n H e h n e --==∑,令()jw H e =()()j w H w e ϕ,H(w)为幅度函数,()w ϕ称为相位函
数。这与模和辐角的表示方法不同,H(w)为可为正可为负的实数,这是为了表达上的方便。如某系统频率响应()jw
e-,如果采用模和幅角的表示方法,sin4w的变
H e=sin4w3j w
号相当于在相位上加上π(因-1= j eπ),从而造成相位曲线的不连贯和表达不方便,用()
H w eϕ则连贯而方便。窗函数法又称傅里叶级数法,其设计是在时域进行的。?函数()j w
一般是无限长且非因果的,设计时需用一个合适的窗函数把它截成有限长的因果序列,使对应的频率响应(的傅里叶变换)尽可能好地逼近理想频率响应。窗函数法的主要缺点是:一、不容易设计预先给定截止频率的滤波器;二、满足同样设计指标的情况下所设计出的滤波器的阶数通常偏大。
一些固定窗函数的特性表
2.FIR数字滤波器设计的基本步骤如下:
(1)确定技术指标
在设计一个滤波器之前,必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中,数字滤波器常被用来实现选频操作。因此,指标的形式一般在频域中给出幅度和相位响应。幅度指标主要以2种方式给出。第一种是绝对指标。他提供对幅度响应函数的要求,一般应用于FIR滤波器的设计。第二种指标是相对指标。他以分贝值的形式给出要求。本文中滤波器的设计就以线性相位FIR滤波器的设计为例。
(2)逼近