基于matlab的fir数字滤波器设计(多通带,窗函数法).doc

数字信号处理课程设计报告设计名称基于MATLAB的FIR数字滤波器设计彪一、课程设计的目的1、通过课程设计把自己在大学中所学的知识应用到实践当中。2、深入了解利用MATLAB设计FIR数字滤波器的基本方法。3、在课程设计的过程中掌握程序编译及软件设计的基本方法。4、提高自己对于新知识的学习能力及进行实际操作的能力。5、锻炼自己通过网络及各种资料解决实际问题的能力。二、主要设计内容利用窗函数法设计FIR滤波器，绘制出滤波器的特性图。利用所设计的滤波器对多个频带叠加的正弦信号进行处理，对比滤波前后的信号时域和频域图，验证滤波器的效果。三、设计原理FIR滤波器具有严格的相位特性，对于信号处理和数据传输是很重要的。目前FIR滤波器的设计方法主要有三种窗函数法、频率取样法和切比雪夫等波纹逼近的最优化设计方法。常用的是窗函数法和切比雪夫等波纹逼近的最优化设计方法。本实验中的窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求高的时候是比较灵活方便的。如果FIR滤波器的HN为实数,而且满足以下任意条件,滤波器就具有准确的线性相位第一种偶对称，HNHN1N，N1/2第二种奇对称，HNHN1N,N1/2PI/2对称中心在NN1/2处四、设计步骤1设计滤波器2所设计的滤波器对多个频带叠加的正弦信号进行处理3比较滤波前后信号的波形及频谱五、用窗函数设FIR滤波器的基本方法基本思路从时域出发设计HN逼近理想HDN。设理想滤波器的单位响应在时域表达为HDN,则HDN一般是无限长的，且是非因果的，不能直接作为FIR滤波器的单位脉冲响应。要想得到一个因果的有限长的滤波器单位抽样响应HN，最直接的方法是先将HDN往右平移，再迕行截断，即截取为有限长因果序列HNHDNWN，并用合适的窗函数迕行加权作为FIR滤波器的单位脉冲响应。按照线性相位滤波器的要求，线性相位FIR数字低通滤波器的单位抽样响应HN必须是偶对称的。对称中心必须等于滤波器的延时常数，即用矩形窗设计的FIR低通滤波器，所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象，且最大波纹大约为幅度的9，返个现象称为吉布斯（GIBBS）效应。为了消除吉布斯效应，一般采用其他类型的窗函数。MATLAB设计FIR滤波器有多种方法和对应的函数。窗函数设计法不仅在数字滤波器的设计中占有重要的地位，同时可以用于功率谱的估计，从根本上讲，使用窗函数的目的就是消除由无限序列的截短而引起的GIBBS现象所带来的影响。六、程序设计RS001FS200采样频率FCUTS1020405060708090A0,1,0,1,0DEVRSONES1,LENGTHAM,WC,BETA,FTYPEKAISERORDFCUTS,A,DEV,FSM为能够满足设计要求的滤波器的最小阶数，WC为滤波器的截止频率点第一个元素F为待设计滤波器的过渡带的起始点和结束点第二个元素A指定第一个元素频率段的理想幅度值第三个元素DEV中的元素为各通带和阻带内允许的幅度最大误差MMODM,2MWINDOWKAISERM1,BETABFIR1M,WC,FTYPE,WINDOW输入的第一个参数是滤波器的阶数第二个参数是滤波器的截止频率第三个参数是滤波器的类型，STOP为带阻滤波器第四个参数是采用的窗函数H,FFREQZB,1,512滤波器的幅频特性图H,WFREQZB,A,N当N是一个整数时函数返回N点的频率向量和幅频响应向量FIGURE1PLOTFFS/2PI,20LOG10ABSH参数分别是频率与幅值XLABEL频率/赫兹YLABEL增益/分贝TITLE滤波器的增益响应F15F220F330F455F575F695待滤波正弦信号频率T0200/FS定义时间的步长S1SIN2F1PITSIN2F2PITSIN2F3PITSS1SIN2F4PITSIN2F5PITSIN2F6PIT滤波前信号SFFILTERB,1,S使用FILTER函数对信号进行滤波FIGURE2SUBPLOT211PLOTT,S滤波前的信号图像XLABEL时间/秒YLABEL幅度TITLE信号滤波前时域图SUBPLOT212FSFFTS,512AFSABSFSFFS/5120255PLOTF,AFS1256滤波前的信号频域图XLABEL频率/赫兹YLABEL幅度TITLE信号滤波前频域图FIGURE3SUBPLOT211PLOTT,SF滤波后的信号图像XLABEL时间/秒YLABEL幅度TITLE信号滤波后时域图SUBPLOT212FSFFFTSF,512滤波后的信号频域图AFSFABSFSF信号频域图的幅值F0255FS/512频率采样PLOTF,AFSF1256滤波后的信号频域图XLABEL频率/赫兹YLABEL幅度TITLE信号滤波后频域图七、运行结果八、实验心得1、对于MATLAB语句有了更加深刻的理解，也注意到了一些运算符号的使用，例如数组的相乘需用来表示，而一般数字相乘应用。还有当运用数组的法时，必须保持数组是等长的，否则，不能相加。2想要改变图形的尺寸，可调用AXIS函数。3在编程过程中应该注意一些细节问题，例如中英文符号的区别，往往一些错误都是由于粗心而导致的。4设计过程中，学习了许多数字信号处理课程中关于数字滤波器的设计的内容，再通过利用参考文献与网络，完成了用MATLAB进行数字信号处理课程设计。5通过课程设计，加深了对课堂抽象概念的理解,巩固了课堂上所学的理论知识,并能很好地理解与掌握数字信号处理中的基本概念、基本原理、基本分析方法。同时掌握编程方法和解决实际问题的技巧。6与其他高级语言的程序设计相比，MATLAB环境下可以更方便、快捷地设计出具有严格线性相位的FIR滤波器，节省大量的编程时间，提高编程效