2022年DSP课程设计FIR高通滤波器设计.pdf

DSP课程设计 FIR 高通滤波器设计FIR 高通滤波器设计南京师范大学物科院精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 1 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计摘要从实现方法方面考虑,将滤波器分为两种 ,一种就是IIR 滤波器 ,另一种就是FIR 滤波器。FIRDF 的最大优点就是可以实现线性相位滤波。而 IIRDF 主要对幅频特性进行逼近 ,相频特性会存在不同程度的非线性。我们知道,无失真传输与滤波处理的条件就是 ,在信号的有效频谱范围内系统幅频响应应为常数,相频响应为频率的线性函数。另外 ,FIR 就是全零点滤波器 ,硬件与软件实现结构简单 ,不用考虑稳定性问题。所以,FIRDF 就是一种很重要的滤波器,在数字信号处理领域得到广泛应用。FIRDF 设计方法主要分为两类:第一类就是基于逼近理想滤波器特性的方法,包括窗函数法、 频率采样法与等波纹最佳逼近法;第二类就是最优设计法。 其中窗函数计法的基本思想就是用FIRDF 逼近希望的滤波特性。本次设计主要采用窗函数设计法 ,对理想滤波器进行逼近 ,从而实现高通滤波器的设计。在 MATLAB软件中,有一系列函数用于设计滤波器,应用时十分方便。因此 ,在本次设计中 ,滤波器的设计主要采用MATLAB软件,编写适当的程序 ,得到滤波器的单位脉冲响应。本设计对滤波器的硬件仿真主要使用CCS软件,通过对滤波器的硬件仿真,可以较为真实的瞧出滤波器的滤波效果。关键字 :高通、 FIRDF 、线性相位、 Hanning 窗、MATLAB 、CCS精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计正文1. 设计目标产生一个多频信号 ,设计一个高通滤波器消除其中的低频成分,通过 CCS 的 graph view波形与频谱显示 ,并与 MATLAB 计算结果比较2. 设计原理1.1 数字滤波器数字滤波器 (digital filter) 就是由数字乘法器、加法器与延时单元组成的一种装置。其功能就是对输入离散信号的数字代码进行运算处理,以达到改变信号频谱的目的。由于电子计算机技术与大规模集成电路的发展,数字滤波器已可用计算机软件实现,也可用大规模集成数字硬件实时实现。 数字滤波器广泛用于数字信号处理中,如电视、VCD、音响等。按照滤波电路的工作频带为其命名:设截止频率为 fp,频率低于 fp 的信号可以通过 ,高于 fp 的信号被衰减的电路称为低通滤波器,频率高于 fp 的信号可以通过 ,低于 fp 的信号被衰减的电路称为高通滤波器;而带通不 ,就就是频率介于低频段截止频率与高频段截止频率的信号可以通过的电路。1.2 高通滤波器高通滤波器就是容许高频信号通过、但减弱(或减少 )频率低于截止频率信号通过的滤波器。对于不同滤波器而言,每个频率的信号的减弱程度不同。它有时被称为低频剪切滤波器 ;在音频应用中也使用低音消除滤波器或者噪声滤波器。高通滤波器与低通滤波器特性恰恰相反。这样的滤波器能够把高频率的声音引导至专用高音喇叭(tweeter),并阻止可能干擾或者损害喇叭的低音信号。使用线圈而不就是电容的低通滤波器也可以同时把低频信号引导至低音喇叭(woofer)。 高通与低通滤波器也用于数字图像处理中在频域中进行变换。1.3 高通滤波器的分析1.3.1高通滤波器的时域分析在时域 ,信号经过系统的响应y (n)体现为激励 x(n)跟系统单位抽样响应h(n)的卷积与 y(n)=(n)h(n)=N1m=0h(m)x(n-m)223 。对于长度为 N 的 FIR 系统, h(n)可以瞧成一个长度为N 点的固定窗口 ,而 x(n)则瞧成一个队列以齐步走的方式穿过 h(n)窗口,每走一步 ,位于窗口中的x(n)部分的点跟h(n)的对应点的值相乘(即加权 )再求与 ,所得结果构成此时系统的响应值y(n), x(n)队列每走一步就得到一个响应值 y(n),即 y(n)就是 h(n)对位于其窗口中的x(n)的加权求与。高通滤波要求 h(n)窗口具有波形锐化作用 ,即利用 h(n)窗口加权与使得变化快的 (即高频 )正弦分量保留 (理想高通 )或衰减幅度小 (实际高通 ) ,而变化缓慢 (即低频 )的正弦分量正负抵消 (理想高通 )或衰减幅度大 (实际高通 ) 。设精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 3 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计其中 N 必须取奇数 , )/10002sin(2 .0)/1002sin(8.0)()()(21ssfnfnnxnxnx其中: fs = 11kHz, n为整数 ,即 x ( n)由 100Hz 的x1(n)与 1kHz 的 x2(n)两种频率的信号组成。高通滤波的目的就就是要尽可能地去掉 x(n)中的低频分量 x1(n) ,同时尽可能地保留x(n)中的高频分量 x2(n)。1.3.2高通滤波器的频域分析在频域 ,信号经过系统的响应y(n)的频谱 Y(ej)体现为激励 x(n)的频谱 X (ej)跟系统单位抽样响应h(n)的频谱 H(ej)(即系统的频谱 )的乘积 Y(ej)=|H(ej)|X(ej)|ejH+x,即响应的频谱 Y(ej)的幅值由系统频谱H(ej)的幅值对激励频谱 X(ej)的幅值相乘 (加权)得到,响应的频谱 Y(ej)的幅角由系统的频谱 H(ej)的幅角跟激励频谱X(ej)的幅角相加 (移相)得到122 。高通滤波要求系统幅度函数 |H(jf)|对需要保留的高频信号频谱加权权重较大(理想时为 1) ,对需要滤除的低频信号频谱加权权重较小(理想时为 0) 。其中 L 为 x(n)的长度 (L=100), 0kL- 1,N 为 h(n)的有值长度 ,m、k 均为整数,跟数字频率 k 相对应的模拟频率为f=fsk/L (Hz)。h(n)在不同 N 值时的频谱如图 2 所示(横轴单位为 kHz),当 N=11 时,在 f=m(kHz)即 f=1kHz、2kHz、3kHz 等处为 1,而在 f=(2m-0、 5)(kHz)即 f=1、 5kHz、 3、 5kHz 等处幅度最大。而在 f=(2m+0、5)(kHz)即 f=2、5kHz、4、5kHz 等处幅度最小 ,如图 2(b);当 N=5 时, 在 f=mfs/5 处为 1,即 f=2、2kHz、3、3kHz 等处为 1,在 f=(4m-1)fs/10 即 f=3、3kHz 处幅度最大。而在 f=(4m+1)fs/10 即 f=5、 5kHz 处幅度最小 ,如图 2(c);当 N=21 时,在 f=mfs/21 即f=524Hz、1047Hz 等处为 1, 在 f=(4m-1)fs/42 处幅度最大 ,而在 f=(4m+1)fs/42 处幅度最小 ,如图 1 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 4 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计图 1 用不同宽度的 h(n)对 x(n)的滤波在频域上表现如图2、2、1 所示,图 2、2、1(a)为 x(n)的频谱 |X(jf)| , 从图 2中可以瞧出 x(n)中含有 100Hz与 1kHz 两种频率的信号,图 1(b)(d)为不同长度的 h(n)对同一 x(n)的滤波情况 ,这跟在时域中分析的结论就是一致的。1.4 FIR 滤波器从实现方法方面考虑 ,将滤波器分为两种 ,一种就是 IIR 滤波器 ,另一种就是 FIR 滤波器。IIR 滤波器保留了模拟滤波器较好的幅度特性,设计简单有效。但这些特性就是以牺牲相位特性为代价而获得的,然而现在许多数据传输,图像处理系统都越来越多的要求系统具有线性相位特性。FIRDF 的最大优点就是可以实现线性相位滤波。而IIRDF 主要对幅频特性进行逼近 ,相频特性会存在不同程度的非线性。我们知道,无失真传输与滤波处理的条件就是,在信号的有效频谱范围内系统幅频响应应为常数,相频响应为频率的线性函数。 另外 ,FIR 就是全零点滤波器 ,硬件与软件实现结构简单,不用考虑稳定性问题。所以,FIRDF 就是一种很重要的滤波器 ,在数字信号处理领域得到广泛应用。FIRDF 设计方法主要分为两类 :第一类就是基于逼近理想滤波器特性的方法,包括窗函数法、频率采样法与等波纹最佳逼近法;第二类就是最优设计法。 其中窗函数计法的基本思想就是用 FIRDF 逼近希望的滤波特性。本次设计主要采用窗函数设计法,对理想滤波器进行逼近,从而实现高通滤波器的设计。1.5 各种窗函数参数各种窗函数的基本参数窗函数类型旁瓣峰值n(dB) 过渡带宽度B阻带最小衰减s(dB) 近似值精确值矩形窗-13 4/N 1、8/N -21 三角窗-25 8/N 6、1/N -25 汉宁窗-31 8/N 6、2/N -44 哈明窗-41 8/N 6、6/N -74 布莱克曼窗-57 12/N 11/N -80 Hanning(汉宁)窗的窗函数为 : 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 5 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计)12cos()()(5. 0)()12cos(15. 0)(NnnRnRnRNnnNNNhn0510152025303500.10.20.30.40.50.60.70.80.91图 2 汉宁窗3. 设计流程1.1 应用 MATLAB设计滤波器 : 为了满足高通滤波的效果,设置线性相位高通FIR 的参数为 :通带截止频率radp2/,阻带截止频率rads4/,通带最大衰减dBp1,阻带最小衰减dBs40。由以上参数 ,通过查表可得 ,hanning 窗与 haming 窗都满足条件 ,我们选择 hanning窗。过渡带宽4/psB,hanning 窗的精确过渡带宽为B=6、2/N,所以要求/4/N6.2B,解之得 N=25。输入程序 : 从而得到满足条件的高通FIRDF: h(n)=-0、0004,-0、0006,0、0028,0、0071,-0、0000,-0、0185,-0、0210,0、0165,0、0624, 0、0355,-0、1061,-0、2898,0、6249,-0、2898,-0、1061,0、0355,0、0624,0、0165, 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计-0、0210,0、0185,-0、0000,0、0071,0、0028,-0、0006,-0、0004 其波形如图 3 所示,频谱如图 4,其幅频特性与相频特性如图5。0102030-0.3-0.2-0.100.10.20.30.40.50.60.7010203000.20.40.60.811.21.4图 3 图 4 00.10.20.30.40.50.60.70.80.91-1500-1000-5000500piPhase(degrees)00.10.20.30.40.50.60.70.80.91-100-50050pi|H(ej|dB)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 7 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计图 5 1.2 用 CCS 对高通滤波器进行仿真设置输入信号 fIn=sin(2n/30)/2+cos(2n*1、 6),使用 MATLAB 运算得到的高通 FIR 滤波器进行滤波 ,滤去低频分量。选择“ Debug ” 菜单的“ RUN ” 项,或按 F12 键运行程序。观察“ fin ”、 “fout ”窗口中时域图形 ;观察滤波效果。鼠标右键单击“ Input ”与“Output” 窗口,选择“ Properties ”项,设置“ Display Type” 为“FFT Magitude ”,再单击“ OK ”按钮结束设置。观察“ Input ”、 “Output” 窗口中频域图形 ;理解滤波效果。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计4. 仿真波形1.1 MATLAB对高通 FIRDF 的设计0510152025-0.3-0.2-0.100.10.20.30.40.50.60.70510152025-1-0.8-0.6-0.4-0.200.20.40.60.811.2 CCS 对高通 FIRDF 的应用仿真精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 9 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 10 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计5. 设计结果由 MATLAB 与 CCS 的仿真结果可以瞧出 ,此高通 FIRDF 可以实现高通滤波的功能 ,在输入为 fIn=sin(2n/30)/2+cos(2n*1、6)时,输出波形中已不含w= 1/30 的低频分量 ,滤波效果良好 ,输出频谱中低频部分基本没有信号。6. 心得通过这次课程设计 ,使我掌握了课程设计的基本思路与方法,掌握了课程设计说明书的基本撰写方法 ,了解了基本的设计思想与设计方法。同时此次课程设计使我课程的基础知识与基本理论有了更深的理解与掌握,锻炼了我们综合运用所学知识的能力,并在理论分析设计、计算制图运用标准与规范查阅设计手册与资料以及计算机应用能了等方面得到了初步的训练与提高 ,培养了我们严谨求实的科学态度。这次课程设计使我进一步熟悉了计算工具软件-MATLAB 、 并进一步掌握了 MATLAB的使用方法、对 MATLAB 语言的发展与特点有了更深的了解,熟悉其工作环境 ,在其应用方面也有了更深的了解 ,了解了 MATLAB 程序书写的一般思路及一般步骤,学会了用 MATLAB解决复杂信号处理等问题的方法。在课程设计的这段时间里,我认为收获还就是很多的,不但进一步掌握了数字信号处理的基础知识及一门专业仿真软件的基本操作,还提高了自己的设计能力及动手能力,同时对于模拟滤波器来了个系统的总结。 更多的就是让我瞧清了自己,明白了凡事需要耐心 ,实践就是检验学习的唯一标准。理论知识的不足在这次课设中表现的很明显。这将有助于我今后的学习 ,端正自己的学习态度 ,从而更加努力的学习。参考文献1 刘泉 阙大顺主编、数字信号处理原理与实现、北京: 电子工业出版社 ,2005 2 程佩青、数字信号处理教程 (第二版 ) 、 北京: 清华大学出版社 ,2001 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 11 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计3 高西全 数字信号处理 - 原理、实现及应用、北京: 电子工业出版社 ,2010精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 12 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计附录MATLAB 程序-生成高通 FIRDF wp=pi/2;ws=pi/4;DB=wp-ws;N0=ceil(6、2*pi/DB);N=N0+mod(N0+1,2);m=1:N;wc=(wp+ws)/2/pi;hn=fir1(N-1,wc,high,hanning(N);fw=abs(fft(hn);figure(1)subplot(1,2,1)plot(m,hn);grid;subplot(1,2,2)plot(m,fw);grid;figure(2)freqz(hn) xlabel(omegapi)ylabel(|H(ejomega|dB)subplot(2,1,2)xlabel(omegapi)CCS 程序#include DSP281x_Device、h / DSP281x Header File #include DSP281x_Examples 、h / DSP281x Examples Include File #include f2812a、h #includemath、h #define FIRNUMBER 25 #define SIGNAL1F 1000 #define SIGNAL2F 4500 #define SAMPLEF 10000 #define PI 3、1415926 float InputWave(); float FIR(); float fHnFIRNUMBER= -0、 0004,-0、 0006,0、 0028,0、 0071,-0、 0000,-0、 0185,-0、0210,0、0165,0、0624, 0、 0355,-0、 1061,-0、 2898,0、 6249,-0、 2898,-0、 1061,0、0355,0、0624,0、0165, -0、0210,0、0185,-0、0000,0、0071,0、0028,-0、0006,-0、0004 ; float fXnFIRNUMBER= 0、0 ; 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 13 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计float fInput,fOutput; float fSignal1,fSignal2; float fStepSignal1,fStepSignal2; float f2PI; int i; float fIn256,fOut256; int nIn,nOut; main(void) nIn=0; nOut=0; f2PI=2*PI; fSignal1=0、0; fSignal2=PI*0、1; fStepSignal1=2*PI/30; fStepSignal2=2*PI*1、6; while ( 1 ) fInput=InputWave(); fInnIn=fInput; nIn+; nIn%=256; fOutput=FIR(); fOutnOut=fOutput; nOut+; if ( nOut=256 ) nOut=0; /* 请在此句上设置软件断点*/ float InputWave() for ( i=FIRNUMBER-1;i0;i- ) fXni=fXni-1; fXn0=sin(fSignal1)/2 、0+cos(fSignal2); fSignal1+=fStepSignal1; if ( fSignal1=f2PI ) fSignal1-=f2PI; fSignal2+=fStepSignal2; if ( fSignal2=f2PI ) fSignal2-=f2PI; return(fXn0); float FIR() float fSum; fSum=0; 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 14 页，共 15 页 - - - - - - - - - - DSP课程设计 FIR 高通滤波器设计for ( i=0;iFIRNUMBER;i+ ) fSum+=(fXni*fHni); return(fSum); 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 15 页，共 15 页 - - - - - - - - - -