数字滤波器的matlab...,fir数字滤波器的matlab实现.doc

fir数字滤波器的matlab实现.doc

综合性、设计性实验报告姓名贺鹤学号201308002118专业通信工程班级通信1301班实验课程名称FIR数字滤波器的MATLAB实现指导教师及职称李玲香讲师开课学期至学年学期上课时间2015年12月24、31日湖南科技学院教务处编印设计题目实验六FIR滤波器的实现与应用小组合作是()否(√)小组成员1实验目的(1)学会调用MATLAB信号处理工具箱中滤波器设计函数设计各种FIR数字滤波器，学会根据滤波需求确定滤波器指标参数。(2)学习用窗函数法设计FIR数字滤波器的原理及其设计步骤；(3)学习编写数字滤波器的设计程序的方法，并能进行正确编程；根据给定的滤波器指标，给出设计步骤；(4)通过观察滤波器输入输出信号的时域波形及其频谱，建立数字滤波的概念。2实验环境(1)WINDOUSXP或、WINDOUS7之类的操作系统版本；(2)MATLAB703实验相关知识(1)窗函数设计线性相位FIR滤波器步骤①确定数字滤波器的性能要求，临界频率，滤波器单位脉冲响应长度N。}{KW②根据性能要求，合理选择单位脉冲响应HN的奇偶对称性，从而确定理想频率响应的幅频特性和相频特性。JWDEH③求理想单位脉冲响应，在实际计算中，可对采样，并对其求NHDJWDEHIDFT的，用代替。NHM④选择适当的窗函数W(N)，根据求所需设计的FIR滤波器NWHND单位脉冲响应。⑤求，分析其幅频特性，若不满足要求，可适当改变窗函数形式或长度JWDEHN，重复上述设计过程，以得到满意的结果。(2)实验用MATLAB函数可以调用MATLAB工具箱函数FIR1实现本实验所要求的线性相位FIRDF的设计，调用一维快速傅立叶变换函数FFT来计算滤波器的频率响应函数。FIR1是用窗函数法设计线性相位FIRDF的工具箱函数，调用格式如下HNFIR1N,WC,‘FTYPE’,WINDOWFIR1实现线性相位FIR滤波器的标准窗函数法设计。HNFIR1N,WC可得到6DB截止频率为WC的N阶单位脉冲响应HN长度为N1FIR低通滤波器，默认缺省参数WINDOWS选用HAMMIING窗。其单位脉冲响应HN满足线性相位条件HNHN1N其中WC为对Π归一化的数字频率，0≤WC≤1。当WC［WC1,WC2］时，得到的是带通滤波器。HNFIR1N,WC,’FTYPE’可设计高通和带阻滤波器。当FTYPEHIGH时，设计高通FIR滤波器；当FTYPESTOP时，设计带阻FIR滤波器。应当注意，在设计高通和带阻滤波器时，阶数N只能取偶数HN长度N1为奇数。不过，当用户将N设置为奇数时，FIR1会自动对N加1。HNFIR1N,WC,WINDOW可以指定窗函数向量WINDOW。如果缺省WINDOW参数，则FIR1默认为HAMMING窗。可用的其他窗函数有BOXCAR,HANNING,BARTLETT,BLACKMAN,KAISER和CHEBWIN窗。例如HNFIR1N,WC,BARTLETTN1使用BARTLETT窗设计；HNFIR1N,WC,CHEBWINN1,R使用CHEBYSHEV窗设计。HNFIR1N,WC,’FTYPE’,WINDOW通过选择WC、FTYPE和WINDOW参数含义同上，可以设计各种加窗滤波器。4．实验内容设抽样频率为FS1000HZ，已知原信号为XSIN2Π80T2SIN2Π140T，由于某种原因，信号被白噪声污染，实际获得的信号为XNXRANDSIZET，要求设计一个FIR滤波器恢复出原始信号。滤波器设计要求频带HZ0,6575,8595,125135,145155,500幅度01010编程实现要求(1)分析实验内容，设计实验方案；(2)绘制XNXRANDSIZET时域波形，观察输出波形，选择合理的滤波器参数；(3)绘制数字滤波器的频率特性曲线；(4)编程实现恢复的X’信号，并画图，与原信号进行比较。5实验设计方案设计两个带通滤波器将原信号进行两次滤波，恢复原信号。用窗函数法设计FIR带通滤波器设计原理在MATLAB中可直接产生窗函数矩形窗，三角窗，汉宁窗，凯塞窗等，通过调用系统的函数即可实现窗的加载。实际应用中选用传函数大多是他们的折中，凯塞窗可以通过改变参数值来折中选择主瓣宽度和旁瓣衰减，基于此的滤波器使用能力强且比较灵活。窗函数设计法是用一定宽度窗函数截取无线脉冲响应序列活得有现场的脉冲响应序列，设计步骤为(1)窗函数法的设计思想是从时域出发，首先确定理想带通滤波器的频率响应函数(2)通过傅立叶变换获得理想带通滤波器的单位脉冲响应HDN(3)由性能指标确定窗函数WN和窗口长度N(4)求带通滤波器的单位脉冲响应HN，HN即为所设计FIR带通滤波器系数向量(5)检验带通滤波器性能6．MATLAB编程实现编程代码T0000101X1SIN2PI80T2SIN2PI140TX2RANDNSIZETXX1X2FB166FB293FC170FC290FS1000WS12PIFB1/FSWP12PIFC1/FSWS22PIFB2/FSWP22PIFC2/FSN95WN1FB1FC1/FSWN2FB2FC2/FSWNWN1,WN2BFIR1N,WN,BLACKMANN1FIGURE1SUBPLOT2,1,1H1,W1FREQZB,1PLOTW1FS/2PI,20LOG10ABSH1GRIDTITLE FIR带通滤波器的幅频特性 MFILTERB,1,XSUBPLOT2,1,2PLOTT,MTITLE 第一次带通滤波器滤波信号 FB11125FB12155FC11130FC12150WS112PIFB11/FSWP112PIFC11/FSWS122PIFB12/FSWP122PIFC12/FSN195WN11FB11FC11/FSWN12FB12FC12/FSWN1WN11,WN12B1FIR1N1,WN1,BLACKMANN11FIGURE2SUBPLOT2,1,1H11,W11FREQZB1,1PLOTW11FS/2PI,20LOG10ABSH11GRIDTITLE FIR带通滤波器的幅频特性 NFILTERB1,1,XSUBPLOT2,1,2PLOTT,NTITLE 第二次带通滤波器滤波信号 YMNFIGURE3SUBPLOT2,1,1PLOTT,Y, R AXIS0040144GRIDTITLE 滤波后的信号 SUBPLOT2,1,2PLOTT,X1, B GRIDAXIS0040144TITLE 滤波的原信号 结果显示图一图二图三图四分析滤波前后信号的变化，波形图变窄。频谱变化表明信号的低频率段和高频率段被滤除，中间频率段被保留。这说明滤波器设计基本符合要求。实验圆满完成。7．运行结果并回答以下问题如果要求用窗函数法设计带通滤波器，且给定上、下边带截止频率为Ω1和Ω2，试求理想带通的单位脉冲响应。NHD答首先确定理想带通滤波器频率响应，然后通过傅立叶变换获得理想带通滤波器的单位脉冲响应SINWCNA/PINA，取A(N1)NHD/2，N66PI/(W1W2)8实验总结主要对本实验开发过程进行归纳和总结，还应包括在设计过程中所遇到的技术难点及解决方法，尚存在的问题以及进一步开发的见解与建议。在课程设计过程