matlab设计 fir 数 字 滤 波 器.doc

  1. 1、本文档共15页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
FIR 数 字 滤 波 器 专业: 学号: 姓名: 一 课题目的: 学会使用Matlab的各项功能。 学会把自己在课堂上学习的知识运用到实践当中。 了解利用Matlab设计FIR数字滤波器的基本方法。 在课程设计的过程中掌握程序编译及软件设计的基本方法。 提高自己对于新知识的学习能力及进行实际操作的能力。 二 课题要求: 在信号处理过程中所处理的信号往往混有噪音,从接受到的信号中消除或减弱噪音是信号处理过程中十分重要的问题。根据有用信号和噪音的不同特性,提取有用信号的过程称为滤波,实现滤波功能的系统称为滤波器。而数字滤波器又是滤波器中运用极为广泛的一种滤波器。数值滤波技术是数字信号处理的一个重要组成部分,滤波器的设计是信号处理的核心问题之一。FIR数字滤波器在保证幅度特性满足技术要求的同时,很容易做到有严格的线性相位特性。要求通过网络及各种资料解决实际问题设计一个符合要求的FIR数字滤波器。 三 课题内容: 数字滤波器和模拟滤波器有着相同的滤波概念,根据其频率响应特性可分为低通、高通、带通、带阻等类型。与模拟滤波器相比,数字滤波器除了具有数字信号处理固有优点外,还有滤波精度高、稳定性好、灵活性强等优点。在数字信号处理中,由于信号中经常混有各种复杂成分,所以很多信号分析都是基于滤波器而进行的, FIR数字滤波器在数字信号处理中发挥着重要作用,采用Matlab软件对FIR数字滤波器进行仿真设计,简化了设计中繁琐的计算。设计中采用窗函数法,频率采样法和优化设计方法,通过调用Matlab函数设计FIR数字滤波器。绘制出滤波器的特性图。利用所设计的滤波器对多个频带叠加的正弦信号进行处理,对比滤波前后的信号时域和频域图,验证滤波器的效果。最后录制一段语音信号,并对录制的信号进行采样和加噪,绘制出采样后语音信号的时域波形和频谱图,然后用所设计的滤波器对加噪后的信号进行滤波,绘制出滤波后信号的时域波形和频谱,并对滤波前后的信号进行对比,分析信号的变化。 窗函数法的Matlab实现 设计FIR数字滤波器的最简单的方法就是窗函数法,通常也称之为傅里叶级数法,FIR数字滤波器的设计首先给出要求的理想滤波器的 频率响应Hd(ejw),设计一个FIR数字滤波器频率响应Hd(ejw),去逼近理想的滤波响应Hd(ejw)。然而窗函数法设计FIR数字滤波器是在时域进行的,因而必须由理想的频率响应Hd(ejw)推导出对应的单位取样响应hd(n),再设计一个FIR数字滤波器的单位取样响应h(n)去逼近hd(n)。 窗函数主要用来减少序列因截断而产生的Gibbs效应,但当这个窗函数为矩形时,得到的FIR数字滤波器幅频响应会有明显的Gibbs效应,并且任意的增加窗函数的长度,Gibbs效应也不能得到改善。为了克服这种现象,窗函数应该使设计的滤波器具有以下几点: 频率特性的主瓣宽度应该尽量窄,且尽可能地将能量集中在主瓣内; 窗函数频率特性的旁瓣在旁瓣ω趋于π的过程中,其能量迅速减小为零。? 程序中fir1函数的用法:b=fir1(n,Wn,’ftype’,window)? ①n为滤波器的阶数? ②Wn为滤波器的截止频率,它是一个0到1的数。如果Wn是一个含有两个数的向量,则函数返回一个带通滤波器? ③ftype为滤波器的类型,ftype=’high’时,设计的是高通滤波器;ftype=’stop’时,设计的是带阻滤波器;没有此参数时,设计的是低通滤波器? ④window为指定的窗函数,矩形窗为boxcar(n),汉宁窗为hanning(n),海明窗为hamming(n),布莱克曼窗为blackman(n),凯撒窗为kaiser(n,beta),没有此参数时,默认为hamming窗 函数程序如下: ?f1=100; f2=200; %待滤波正弦信号频率? fs=2000; %采样频率? m=(0.3*f1)/(fs/2); %定义过度带宽 ?M=round(8/m); %定义窗函数的长度 ?N=M-1; %定义滤波器的阶数? b=fir1(N,0.5*f2/(fs/2)); %使用fir1函数设计滤波器?%输入的参数分别是滤波器的阶数和截止频率 ?figure(1)? [h,f]=freqz(b,1,512); %滤波器的幅频特性图? ?plot(f*fs/(2*pi),20*log10(abs(h))) %参数分别是频率与幅值? xlabel('频率/赫兹'); ylabel('增益/分贝'); title('滤波器的增益响应'); ?figure(2)? subplot(211

您可能关注的文档

文档评论(0)

pkaokqunw + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档