下载

2下载券

加入VIP
  • 专属下载特权
  • 现金文档折扣购买
  • VIP免费专区
  • 千万文档免费下载

上传资料

关闭

关闭

关闭

封号提示

内容

首页 基于Hamming窗的FIR数字滤波器设计

基于Hamming窗的FIR数字滤波器设计.doc

基于Hamming窗的FIR数字滤波器设计

距离会产生美只有时间不要太长
2019-05-19 0人阅读 举报 0 0 0 暂无简介

简介:本文档为《基于Hamming窗的FIR数字滤波器设计doc》,可适用于IT/计算机领域

郑州航空工业管理学院《电子信息系统仿真》课程设计级电子信息工程专业    班级题  目基于Hamming窗的FIR数字滤波器设计姓  名      学  号    二О一六年十一月二十五日一、MATLAB软件简介MATLAB是矩阵实验室(MatrixLaboratory)之意。除具备卓越的数值计算能力外它还提供了专业水平的符号计算文字处理可视化建模仿真和实时控制等功能。MATLAB的基本数据单位是矩阵它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完相同的事情简捷得多。MATLAB拥有数百个内部函数的主包和三十几种工具包(Toolbox)工具包又可以分为功能性工具包和学科工具包功能工具包用来扩充MATLAB的符号计算,可视化建模仿真,文字处理及实时控制等功能学科工具包是专业性比较强的工具包,控制工具包,信号处理工具包,通信工具包等都属于此类 开放性使MATLAB广受用户欢迎除内部函数外,所有MATLAB主包文件和各种工具包都是可读可修改的文件,用户通过对源程序的修改或加入自己编写程序构造新的专用工具包。语言特点一种语言之所以能如此迅速地普及显示出如此旺盛的生命力是由于它有着不同于其他语言的特点正如同FORTRAN和C等高级语言使人们摆脱了需要直接对计算机硬件资源进行操作一样被称作为第四代计算机语言的MATLAB利用其丰富的函数资源使编程人员从繁琐的程序代码中解放出来。MATLAB最突出的特点就是简洁。MATLAB用更直观的符合人们思维习惯的代码代替了C和FORTRAN语言的冗长代码。MATLAB给用户带来的是最直观最简洁的程序开发环境。、理论分析.FIR数字滤波器结构特征FIR数字滤波器的冲激响应只能延续一定的时间,N阶因果有限冲激响应滤波器的差分方程表达式为:。其卷积和是一个有限和,可以直接计算,所以其涉及的基本操作就是简单的乘和加,其计算仅需要所给出的初始样本值及所要求的相关样本值。通过精确设计,有限冲激响应滤波器在整个频率范围内均能提供精确的线性相位,由于系统的单位脉冲序列为有限序列,当输入有限时,输出也必然为有限,这样其稳定性总可以独立于滤波器系数之外。因此,在很多情况下,有限冲激响应滤波器成为首选,只要确定能满足要求的转移序列或者脉冲响应的常数,就可以准确地设计出满足要的FIR数字滤波器。通常所采用的设计方法主要有窗函数、频率采样法和等波纹最佳逼近法等,其中窗函数法是从时域进行设计的,其算法简单、物理意义清晰,因此得到了较为广泛的应用。.算法研究为了建立一个具有线性相位和稳定的非递归特性的有限脉冲响应滤波器即FIR滤波器,要考虑两个方面:一是使用有限长的单位取样响应来逼近理想低通二是单位取样响应对(N)对称,保证线性相位根据FIR数字滤波器输入输出关系的差分方程描述,有上式反映了有限列长的单位取样响应。滤波网络传输函数为:利用傅立叶变换得单位取样响应为是截止频率,从而保证线性相位。在窗函数法设计中使用有限列长的逼近用Hamming窗函数将截断进行加权处理后可得。此时频率响应函数为:Hamming窗函数。其中是幅度为,长度为N的矩形序列。.窗函数法的原理可以从时域或频域出发来设计FIR滤波器,从频域出发的方法称为频率采样设计法,从时域出发的方法称为窗函数法,也称为傅立叶级数法,这种方法应用较广泛。窗函数法的设计思想是按照所要求的理想滤波器频率响应,设计一个FIR滤波器,使之频率响应来逼近。因为设计是在时域中进行的,需先由的傅立叶反变换导出序列即由于是矩形频率特性有是一无限长的序列,且是非因果的,而要设计的FIR滤波器的冲激响应序列是有限长的,所以要用有限长的序列来逼近无限长的序列最有效的方法是截断或者说用一个有限长度的窗口函数序列来截取即按照复卷积公式,在时域中的乘积关系可表示成在频域中的周期性卷积关系,即可得所设计的FIR滤波器的频率响应其中,为截断窗函数的频率特性。由此可见,实际的FIR数字滤波器的频率响应逼近理想滤波器频率响应的好坏完全取决于窗函数的频率特性。.窗函数法的设计步骤()给定所要求的频率响应()由的傅立叶反变换求出()根据过渡带宽及阻带最小衰减,选择窗的形状并估计窗口长度N()计算滤波器的单位脉冲响应()验算技术指标是否满足要求。()窗谱主辨尽可能窄,以获得较陡的过渡带()尽量减少窗谱的最大旁瓣的相对幅度,使能量尽量集中于主辨,进而增加阻带的衰减。给定的滤波器指标参数一般为通带截止频率、阻截止频率、实际通带波动和最小阻带衰减。窗函数设计的经验公式为归一化过渡带  滤波器阶数   当≥时当≤≤时三、设计内容.低通滤波器的设计指标。通带截止频率,阻带截止频率,最小阻带衰减,抽样频率为Hz.设计分析由以上参数可得过度带宽trwidth=wswp,wc=(wswp)。考虑选择Hamming窗函数则列长N=*pitrwidth,hd(n)=sin(wc*(n(N)))(pi*(n(N))),h(n)=hd(n)*wh(n)。利用MATLAB编程时,先计算N、trwidth、hd(n)和h(n),并利用MATLAB提供的相关工具计算该低通滤波器的幅度特性值,最终显示Hamming窗频谱、实际单位取样响应频谱、理想单位取样响应频谱和该FIR滤波器幅度特性曲线。.编写Matlab程序()实现FIR滤波器的主程序:wp=*piws=*pi 给出通带和阻带的截止频率trwidth=wswp    得出过渡带宽度N=ceil(*pitrwidth) 滤波器的长度n=::Nwc=(wswp)      计算理想低通截止频率hd=ideallp(wc,N)     计算理想低通单位冲激响应wham=(hamming(N))'  求海明函数h=hd*wham      应用窗函数加权后的有限长单位冲激响应序列db,mag,pha,grd,w=freqzm(h,) 幅检查值响应、相位响应和群延迟响应deltaw=*piAp=(min(db(::wpdeltaw)))    检查通带最大衰减Ar=round(max(db(wsdeltaw::))) 检查阻带最小衰减subplot()stem(n,hd)title('理想单位脉冲响应hd(n)')subplot()plot(n,wham)title('海明窗w(n)')subplot()stem(n,h)title('实际单位脉冲响应h(n)')subplot()plot(wpi,db)title('幅度响应(dB)')n=:Fs=      抽样频率dt=Fst=n*dtTp=*dtf=f=x=sin(*pi*f*t)sin(*pi*f*t)y=filter(wham,,x)figure()subplot(,,)plot(t,x)title('输入信号x(t)')xlabel('ts')ylabel('x(t)')axis(,Tp,min(x),max(x))subplot(,,)plot(t,y)title('输出信号y(t)')xlabel('ts')ylabel('y(t)')axis(,Tp,min(y),max(y))()计算理想低通的单位冲激响应的程序:functionhd=ideallp(wc,N)理想线性相位低通滤波器计算hd=ideallp(wc,n) hd=~N~之间的理想脉冲响应 wc=截止频率(弧度) N=理想线性相位滤波器的长度tao=(N)n=::(N)m=ntaoeps     加一个极小数eps以避免hd的分母为零hd=sin(wc*m)(pi*m)()检查幅值响应、相位响应和群延迟响应functiondb,mag,pha,grd,w=freqzm(b,a)H,w=freqz(b,a,,'whole')H=(H(::))'w=(w(::))'mag=abs(H)db=*log((mageps)max(mag))pha=angle(H)grd=grpdelay(b,a,w)输出结果如下图所示:由上图可以看到:通带纹波为dB最小阻带衰减为dB。满足设计指标的要求。把滤波器的阶数人为的设置成(将N=ceil(*pitrwidth)改为N=)由上图可以看到:通带纹波为dB最小阻带衰减为dB。把滤波器的阶数设置为由上图可以看到:通带纹波为dB最小阻带衰减为dB。不能满足指标要求。综上所述:当滤波器的阶数增加时可以减小通带纹波增加最小阻带衰减。但同时计算量也随之增大成本也要同样增加。用MATLAB中的FilterDesignandAnalysis工具进行同样的设计设置阶数为低通窗口选择hamming归一化频率。 幅频特性曲线相频特性曲线冲击响应零极点图四、总结通过本次课程设计使我更加扎实的掌握了有关数字滤波器方面的知识在这次课程设计中虽然遇到了很多的困难但是经过一次又一次的查阅资料一次一次的认真思考终于找出了原因通过这次课程设计了解知识还是需要通过实践才能真正的掌握住精髓。通过这次学习我不但掌握了FIR数字滤波器窗函数的基本知识及其实际应用的技巧了而且还学会了一些Matlab编程语言对这款软件也有了更深刻的认识得到了很大的收获。我认为本次课程设计对我们来说是很有用的能够提高我们独立去做一个工程的能力体会将来我们要面对的实际问题既巩固了在课堂上学习的知识也在查找资料的过程中也学到了很多课本上没有学习到的知识可以说收获颇丰通过这次课程设计将理论与实践相结合从中提高而来自己独立思考能力和实践动手能力从而不断超越自己学会了坚持不懈坚持到底。这对日后我们走进社会了参加工作非常重要的。

用户评价(0)

关闭

新课改视野下建构高中语文教学实验成果报告(32KB)

抱歉,积分不足下载失败,请稍后再试!

提示

试读已结束,如需要继续阅读或者下载,敬请购买!

评分:

/16

VIP

在线
客服

免费
邮箱

爱问共享资料服务号

扫描关注领取更多福利