dsp课程设计-fir高通滤波器设计

FIR高通滤波器设计南京师范大学物科院摘要从实现方法方面考虑，将滤波器分为两种，一种是IIR滤波器，另一种是FIR滤波器。FIRDF的最大优点是可以实现线性相位滤波。而IIRDF主要对幅频特性进行逼近，相频特性会存在不同程度的非线性。我们知道，无失真传输与滤波处理的条件是，在信号的有效频谱范围内系统幅频响应应为常数，相频响应为频率的线性函数。另外，FIR是全零点滤波器，硬件和软件实现结构简单，不用考虑稳定性问题。所以，FIRDF是一种很重要的滤波器，在数字信号处理领域得到广泛应用。FIRDF设计方法主要分为两类第一类是基于逼近理想滤波器特性的方法，包括窗函数法、频率采样法和等波纹最佳逼近法；第二类是最优设计法。其中窗函数计法的基本思想是用FIRDF逼近希望的滤波特性。本次设计主要采用窗函数设计法，对理想滤波器进行逼近，从而实现高通滤波器的设计。在MATLAB软件中，有一系列函数用于设计滤波器，应用时十分方便。因此，在本次设计中，滤波器的设计主要采用MATLAB软件，编写适当的程序，得到滤波器的单位脉冲响应。本设计对滤波器的硬件仿真主要使用CCS软件，通过对滤波器的硬件仿真，可以较为真实的看出滤波器的滤波效果。关键字高通、FIRDF、线性相位、HANNING窗、MATLAB、CCS正文1设计目标产生一个多频信号，设计一个高通滤波器消除其中的低频成分，通过CCS的GRAPHVIEW波形和频谱显示，并和MATLAB计算结果比较2设计原理21数字滤波器数字滤波器DIGITALFILTER是由数字乘法器、加法器和延时单元组成的一种装置。其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。由于电子计算机技术和大规模集成电路的发展，数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。数字滤波器广泛用于数字信号处理中，如电视、VCD、音响等。按照滤波电路的工作频带为其命名设截止频率为FP，频率低于FP的信号可以通过，高于FP的信号被衰减的电路称为低通滤波器，频率高于FP的信号可以通过，低于FP的信号被衰减的电路称为高通滤波器；而带通吗，就是频率介于低频段截止频率和高频段截止频率的信号可以通过的电路。22高通滤波器高通滤波器是容许高频信号通过、但减弱（或减少）频率低于截止频率信号通过的滤波器。对于不同滤波器而言,每个频率的信号的减弱程度不同。它有时被称为低频剪切滤波器；在音频应用中也使用低音消除滤波器或者噪声滤波器。高通滤波器与低通滤波器特性恰恰相反。这样的滤波器能够把高频率的声音引导至专用高音喇叭TWEETER,并阻止可能干擾或者损害喇叭的低音信号。使用线圈而不是电容的低通滤波器也可以同时把低频信号引导至低音喇叭WOOFER。高通和低通滤波器也用于数字图像处理中在频域中进行变换。23高通滤波器的分析231高通滤波器的时域分析在时域,信号经过系统的响应YN体现为激励XN跟系统单位抽样响应HN的卷积和YNNHNN1M0HMXNM223。对于长度为N的FIR系统,HN可以看成一个长度为N点的固定窗口,而XN则看成一个队列以齐步走的方式穿过HN窗口,每走一步,位于窗口中的XN部分的点跟HN的对应点的值相乘即加权再求和,所得结果构成此时系统的响应值YN,XN队列每走一步就得到一个响应值YN,即YN是HN对位于其窗口中的XN的加权求和。高通滤波要求HN窗口具有波形锐化作用,即利用HN窗口加权和使得变化快的即高频正弦分量保留理想高通或衰减幅度小实际高通,而变化缓慢即低频的正弦分量正负抵消理想高通或衰减幅度大实际高通。设其中N必须取奇数,/102SIN/102SIN821SSFNFXNX其中FS11KHZ,N为整数,即XN由100HZ的X1N和1KHZ的X2N两种频率的信号组成。高通滤波的目的就是要尽可能地去掉XN中的低频分量X1N,同时尽可能地保留XN中的高频分量X2N。232高通滤波器的频域分析在频域,信号经过系统的响应YN的频谱YEJ体现为激励XN的频谱XEJ跟系统单位抽样响应HN的频谱HEJ即系统的频谱的乘积YEJ|HEJ|XEJ|EJHX,即响应的频谱YEJ的幅值由系统频谱HEJ的幅值对激励频谱XEJ的幅值相乘加权得到,响应的频谱YEJ的幅角由系统的频谱HEJ的幅角跟激励频谱XEJ的幅角相加移相得到122。高通滤波要求系统幅度函数|HJF|对需要保留的高频信号频谱加权权重较大理想时为1,对需要滤除的低频信号频谱加权权重较小理想时为0。其中L为XN的长度L100,0KL1,N为HN的有值长度,M、K均为整数,跟数字频率K相对应的模拟频率为FFSK/LHZ。HN在不同N值时的频谱如图2所示横轴单位为KHZ,当N11时,在FMKHZ即F1KHZ、2KHZ、3KHZ等处为1,而在F2M05KHZ即F15KHZ、35KHZ等处幅度最大。而在F2M05KHZ即F25KHZ、45KHZ等处幅度最小,如图2B当N5时,在FMFS/5处为1,即F22KHZ、33KHZ等处为1,在F4M1FS/10即F33KHZ处幅度最大。而在F4M1FS/10即F55KHZ处幅度最小,如图2C当N21时,在FMFS/21即F524HZ、1047HZ等处为1,在F4M1FS/42处幅度最大,而在F4M1FS/42处幅度最小,如图1图1用不同宽度的HN对XN的滤波在频域上表现如图221所示,图221A为XN的频谱|XJF|,从图2中可以看出XN中含有100HZ和1KHZ两种频率的信号,图1BD为不同长度的HN对同一XN的滤波情况,这跟在时域中分析的结论是一致的。24FIR滤波器从实现方法方面考虑，将滤波器分为两种，一种是IIR滤波器，另一种是FIR滤波器。IIR滤波器保留了模拟滤波器较好的幅度特性，设计简单有效。但这些特性是以牺牲相位特性为代价而获得的，然而现在许多数据传输，图像处理系统都越来越多的要求系统具有线性相位特性。FIRDF的最大优点是可以实现线性相位滤波。而IIRDF主要对幅频特性进行逼近，相频特性会存在不同程度的非线性。我们知道，无失真传输与滤波处理的条件是，在信号的有效频谱范围内系统幅频响应应为常数，相频响应为频率的线性函数。另外，FIR是全零点滤波器，硬件和软件实现结构简单，不用考虑稳定性问题。所以，FIRDF是一种很重要的滤波器，在数字信号处理领域得到广泛应用。FIRDF设计方法主要分为两类第一类是基于逼近理想滤波器特性的方法，包括窗函数法、频率采样法和等波纹最佳逼近法；第二类是最优设计法。其中窗函数计法的基本思想是用FIRDF逼近希望的滤波特性。本次设计主要采用窗函数设计法，对理想滤波器进行逼近，从而实现高通滤波器的设计。25各种窗函数参数各种窗函数的基本参数过渡带宽度B窗函数类型旁瓣峰值DBN近似值精确值阻带最小衰减DBS矩形窗134/N18/N21三角窗258/N61/N25汉宁窗318/N62/N44哈明窗418/N66/N74布莱克曼窗5712/N11/N80HANNING（汉宁）窗的窗函数为12COS5012COS150NNRNNRNNNNNH0510152025303500102030405060708091图2汉宁窗3设计流程31应用MATLAB设计滤波器为了满足高通滤波的效果，设置线性相位高通FIR的参数为通带截止频率，阻带截止频率，通带最大衰减，阻带最RADP2/RADS4/DBP1小衰减。BS40由以上参数，通过查表可得，HANNING窗和HAMING窗都满足条件，我们选择HANNING窗。过渡带宽，HANNING窗的精确过渡带宽为4/PS62/N，所以要求，解之得N25。BN62B输入程序从而得到满足条件的高通FIRDFHN00004,00006,00028,00071,00000,00185,00210,00165,00624,00355,01061,02898,06249,02898,01061,00355,00624,00165,00210,00185,00000,00071,00028,00006,00004其波形如图3所示，频谱如图4，其幅频特性和相频特性如图5。0102030030201001020304050607010203000204060811214图3图400102030405060708091150010005000500PIPHASEDEGREES0010203040506070809110050050PI|HEJ|DB图532用CCS对高通滤波器进行仿真设置输入信号FINSIN2N/30/2COS2N16,使用MATLAB运算得到的高通FIR滤波器进行滤波，滤去低频分量。选择“DEBUG”菜单的“RUN”项，或按F12键运行程序。观察“FIN”、“FOUT”窗口中时域图形；观察滤波效果。鼠标右键单击“INPUT”和“OUTPUT”窗口，选择“PROPERTIES”项，设置“DISPLAYTYPE”为“FFTMAGITUDE”，再单击“OK”按钮结束设置。观察“INPUT”、“OUTPUT”窗口中频域图形；理解滤波效果。4仿真波形41MATLAB对高通FIRDF的设计05101520250302010010203040506070510152025108060402002040608142CCS对高通FIRDF的应用仿真5设计结果由MATLAB和CCS的仿真结果可以看出，此高通FIRDF可以实现高通滤波的功能，在输入为FINSIN2N/30/2COS2N16时，输出波形中已不含W1/30的低频分量，滤波效果良好，输出频谱中低频部分基本没有信号。6心得通过这次课程设计，使我掌握了课程设计的基本思路和方法，掌握了课程设计说明书的基本撰写方法，了解了基本的设计思想和设计方法。同时此次课程设计使我课程的基础知识和基本理论有了更深的理解和掌握，锻炼了我们综合运用所学知识的能力，并在理论分析设计、计算制图运用标准和规范查阅设计手册与资料以及计算机应用能了等方面得到了初步的训练和提高，培养了我们严谨求实的科学态度。这次课程设计使我进一步熟悉了计算工具软件MATLAB并进一步掌握了MATLAB的使用方法对MATLAB语言的发展和特点有了更深的了解，熟悉其工作环境,在其应用方面也有了更深的了解，了解了MATLAB程序书写的一般思路及一般步骤，学会了用MATLAB解决复杂信号处理等问题的方法。在课程设计的这段时间里，我认为收获还是很多的，不但进一步掌握了数字信号处理的基础知识及一门专业仿真软件的基本操作，还提高了自己的设计能力及动手能力，同时对于模拟滤波器来了个系统的总结。更多的是让我看清了自己，明白了凡事需要耐心，实践是检验学习的唯一标准。理论知识的不足在这次课设中表现的很明显。这将有助于我今后的学习，端正自己的学习态度，从而更加努力的学习。参考文献1刘泉阙大顺主编数字信号处理原理与实现北京电子工业出版社，20052程佩青数字信号处理教程（第二版）北京清华大学出版社，20013高西全数字信号处理原理、实现及应用北京电子工业出版社，2010附录MATLAB程序生成高通FIRDFWPPI/2WSPI/4DBWPWSN0CEIL62PI/DBNN0MODN01,2M1NWCWPWS/2/PIHNFIR1N1,WC,HIGH,HANNINGNFWABSFFTHNFIGURE1SUBPLOT1,2,1PLOTM,HNGRIDSUBPLOT1,2,2PLOTM,FWGRIDFIGURE2FREQZHNXLABELOMEGAPIYLABEL|HEJOMEGA|DBSUBPLOT2,1,2XLABELOMEGAPICCS程序INCLUDE“DSP281X_DEVICEH“/DSP281XHEADERFILEINCLUDEFILEINCLUDE“DSP281X_EXAMPLESH“/DSP281XEXAMPLESINCLUDEFILEINCLUDE“F2812AH“INCLUDE“MATHH“DEFINEFIRNUMBER25DEFINESIGNAL1F1000DEFINESIGNAL2F4500DEFINESAMPLEF10000DEFINEPI31415926FLOATINPUTWAVEFLOATFIRFLOATFHNFIRNUMBER00004,00006,00028,00071,00000,00185,00210,00165,00624,00355,01061,02898,06249,02898,01061,00355,00624,00165,00210,00185,00000,00071,00028,00006,00004FLOATFXNFIRNUMBER00FLOATFINPUT,FOUTPUTFLOATFSIGNAL1,FSIGNAL2FLOATFSTEPSIGNAL1,FSTEPSIGNAL2FLOATF2PIINTIFLOATFIN256,FOUT256INTNIN,NOUTMAINVOIDNIN0NOUT0F2PI2PIFSIGNAL100FSIGNAL2PI01FSTEPSIGNAL12PI/30FSTEPSIGNAL22PI16WHILE1FINPUTINPUTWAVEFINNINFINPUTNINNIN25