DSP论文dw.讲义

《DSP原理及应用》课程小论文学院:_XXXXXXXXXXXXXXXXXX_专业班级:电自普XXXX学生姓名:XX学号:XXXXXXXXXX论文题目:基于DSP开发板的语言信号滤波处理完成日期：2013年6月10日指导教师评语:______________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________成绩（五级记分制）:________________指导教师（签字）:________________完成日期：2013年6月10日1指导教师评语:___________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________________成绩（五级记分制）:________________指导教师（签字）:________________摘要数字信号处理技术及其应用，目前正以惊人的速度向前发展着。随着大规模集成电路的出现和数字部件的成本下降、体积缩小及运算速度提高，数字信号处理的应用日益广泛。目前已制成多种专用数字滤波器，取样率可高达兆赫。高速专用快速傅里叶变换处理机已有商品出售。简单的数字滤波器已制成集成电路片。目前几乎所有的语音带宽压缩系统都倾于全数字化，因为目前它是最实际可行的方法。除了专用数字信号处理硬件有所发展之外，还出现了可编程序的数字信号处理专用计算机，这种计算机的构造特别适于解决数字信号处理问题。它目前应用于实时信号处理以及设计、模拟专用数字硬件。在计算量相等的情况下IIR数字滤波器比FIR滤波器的幅频特性优越，频率选择性也好，但是它有着致命的缺点相位特性不好控制，它的相位特性是的非线性函数。例如双线性变换法产生的IIR滤波器模拟指标的频率与数字化指标2的频率转换关系是，这是使频率产生严重的非线性的原因，这种W与w的非线性关系，使数字滤波器与模拟滤波器在响应与频率的对应关系上发生了畸变，如果需要线性相位，就必须用全通网络进行复杂的相位校正但是，在对程序运行周期数要求十分严格的DSP处理中加上一个全通均衡器是十分浪费资源的，另外即使加上全通均衡器对于因果的IIR滤波器仍将得不到线性的相位。在现代电子系统中如图像处理数据传输等波形传递系统中都越来越多的要求信道具有线性的相位特性在这方面FIR滤波器具有独到的优点，它可以在幅度特性随意设计的同时保证精确严格的线性相位。本文主要介绍了利用DSP来实现FIR滤波器的设计，根据FIR滤波器设计的要求，本设计以CCS为开发环境，采用模块化的设计方案来完成设计。本文详细介绍了硬件电路和软件设计的方法，硬件方面介绍了SEED-DEC2812开发实验箱，软件设计方面包括了程序各模块的流程图和源程序代码，还详细介绍了CCS为开发环境，并详细介绍了仿真的步骤。关键词:MATLABDSPFIR数字滤波器1.要求和功能实现1.1目的要求通过课程设计，加深对DSP芯片TMS320C54x的结构、工作原理的理解，获得DSP应用技术的实际训练，掌握设计较复杂DSP系统的基本方法。通过使用汇编语言编写具有完整功能的图形处理程序或信息系统，使学生加深对所学知识的理解，进一步巩固汇编语言讲法规则。学会编制结构清晰、风格良好、数据结构适当的汇编语言程序，从而具备解决综合性实际问题的能力。通过使用Matlab设计相应的FIR低通数字滤波器，得到滤波器H(z)的系数，然后根据这些系数，编写关于DSP的汇编程序。在CCS软件上对已采集信号进行处理，最后滤除高频信号，输出我们需要的低频语音信号。效果本文采取比较滤波前后信号的频谱图。1.2功能实现1.2.1设计滤波器。参数指标：根据被处理对象设计滤波器设计指标。31.2.2用汇编语言在TMSc54xx编程实现符合要求的滤波器。1.2.3通过开发板上的语音采集模块输入待处理语音信号(简单介绍)。1.2.4滤除语音信号中的高频噪声。1.2.5显示滤波器的频谱。1.2.6通过开发板上的语音模块播放处理后的语音信号（简单介绍）。2.系统参数和总体设计方案2.1系统参数一个实际的应用系统中，总存在各种干扰。数字滤波器在语音信号处理、信号频谱估计、信号去噪、无线通信中的数字变频以及图像信号等各种信号处理中都有广泛的应用，数字滤波器也是使用最为广泛的信号处理算法之一。在本设计中，使用MATLAB模拟产生合成信号，然后利用CCS进行滤波。设定模拟信号的采样频率为8000Hz，。设计一个FIR低通滤波器，其参数为：滤波器名称：FIR低通滤波器采样频率：Fs=8000Hz通带/阻带截止频率：2000Hz～2500Hz通带最大衰减：0.5dB阻带最少衰减：50dB滤波器级数：N=25滤波器系数：由MATLAB根据前述参数求得。2.2系统总体设计方案流程图43.DSP及其开发环境3.1DSP系统的构成一个典型的DSP系统如图2示。图2典型的DSP系统图2是一个用DSP做信号处理的典型框图。由于DSP是用来对数字信号进行处理的，所以首先必须将输入的模拟信号变换为数字信号。3.2DSP系统的特点及设计过程由于数字信号处理系统是以数字信号处理理论为基础，所以具有数字信号处理的全部优点：接口方便；编程方便；具有高速性；稳定性好；精度高；可重复性好；集成方便等等，同时也存在一定的缺点一般来说DSP的设计过程应遵循一定的设计流程，如图3示。5图3DSP基本设计流程3.3CCS开发环境CCS全称是CodeComposerStudio它提供了基本的代码生成工具，具有一定的调试、分析能力，在CCS下的程序开发过程如图4示。图4程序开发过程CCS包括：1、CCS代码生成工具2、CCS集成开发环境IDE3、DSP/BIOS插件程序和APIRTDX插件，主机接口和API等。4.数字滤波器的设计原理和方案对比4.1IIR滤波器无限长冲激响应系统IIR（InfiniteImpulseResponse）中，用有理分式表示的系统函数来逼近所需要的频率响应，即其单位冲激响应h(n)是无限长的。DSP应用定义系统性能指标硬件调试选择DSP芯片软件编程硬件设计系统集成软件调试系统测试和调试6IIR滤波器由于吸收了模拟滤波器的结果，有大量的图表可查，可以方便、简单、有效地完成设计，效果很好，但是其相位特性不好控制，必须用全通网络进行复杂的相位较正，才能实现线性相位特性的要求。4.2FIR滤波器有限长冲激响应系统FIR（FiniteImpulseResponse）中，用一个有理多项式表示的系统函数去逼近所需要的频率响应，即其单位冲激响应h(n)在有限个n值处不为零。FIR滤波器可在幅度特性随意设计的同时，保证精确、严格的线性相位特性。这在要求相位线性信道的现代电子系统，如图像处理、数据传输等波形传递系统中，是具有很大吸引力的。而且，其单位冲激响应是有限长的，不存在不稳定的因素，并且可用因果系统来实现。4.3FIR滤波器的基本结构数字滤波是将输入的信号序列，按规定的算法进行处理，从而得到所期望的输出序列。一个线性位移不变系统的输出序列y(n)和输入x(n)之间的关系，应满足常系数线性差分方程，见公式4.1，MiiNiiinyainxbny1100n（4.3.1）其中，x(n)为输入序列，y(n)为输出序列，kkba和为滤波器系数，N是滤波器的阶数。若上式中所有的kb均为零，则有FIR滤波器的差分方程为：10)()(Nkkknxany（4.3.2）对上式进行Z变换得到FIR滤波器的传递函数为：10NikkzbzXzYzH（4.3.3）由上式可以看出，H(z)是1Z的N-1次多项式，它在z平面内有N-1个零点，同时在原点处有N-1个重极点。N阶滤波器通常采用N个延迟单元、N个加法器与N+1个乘法器，取图5中(a)、(b)两种结构。图5FIR滤波器的一般结构通常FIR滤波器的设计方法主要有三种：窗函数法、频率抽样法和切比雪夫7等波纹逼近法。其中窗函数法可以应用比较现成的窗函数，因而设计简单，在指标要求不高的场合使用方便灵活。4.3.1窗函数法窗函数法也称为傅立叶级数法。理想的数字滤波器频率特性是无法实现的，FIR的设计就是要寻找一个可以得到的频率特性=来逼近，这相当于用一个可实现的单位脉冲响应h(n)去逼近一个理想单位脉冲响应。可由理想频率特性通过傅氏反变换得到一般来说，这样得到的理想单位脉冲响应序列是个无限长序列，因而是非因果的。设有一个截止频率为的理想线性相位低通，延时为τ，其频率特性是：得到：这是一个以ｎ＝τ为中心偶对称的无限长非因果序列，要想用一个有限长的因果序列去逼近它，最简单的方法是截取ｎ从0到N-1的一段来表示它，即h(n)=)(nhd)10(Nn；其他N：h(n)=0。同时，为了保证线性相位，还要满足偶对称h(n)=h(N-1-n)。这就好像通过一个窗口观看到的一段)(nhd，因此h(n)就表示成)(nhd和一个窗口函数的乘积，这样对h(n)的求解就变为h(n)＝)(nhd*nW，这里的nW就称为窗口函数,既然一个频域上的标准的矩形窗口对应于时域是一个无限长的序列,那么在时域上截取一段势必造成频域的矩形窗口的失真。结果就是截取出的信号也相应失真，为了补偿这种失真，只有改变原来窗口的形状，修正经过时域截取后的窗口失真。窗函数设计方法的基本步骤是：(1)把jwdeH展成FS，得)(nhd；(2)对)(nhd自然截短到所需的长度，如2M+1；(3)将截短后的)(nhd右移M个采样间隔，得h(n)；(4)将h(n)乘以合适的窗口，即得所要滤波器的冲击响应，窗函数以n=M对8称。利用所求得的单位抽样响应，即可用硬件构成滤波器的转移函数H(z)，也可利用h(n)在计算机上用软件来实现滤波。4.3.2频率抽样法窗函数法是从时域出发，用窗函数截取理想的nhd得到h(n)，以此有限长的h(n)近似nhd，这样得到的频率响应jeH逼近于理想的频响jwdeH。频率抽样法是从频率出发，将给定的理想频响jwdeH加以等间隔抽样。jwdeHkHdNk2然后以此kHd作为FIR滤波器的频率响应抽样值H(k)，再根据DFT（离散付氏变换）定义由频域这N个抽样值来唯一确定一个有限长序列h(n)，同样也可以算出FIR滤波器的系统函数H(z)及频率响应jeH，可以推出频率响应jeH是频率抽样值H(k)与线性相位因子2/1Nje及如下内插函数S(ω,k)的线性组合。所以，在各频率取样点上，实际滤波器的频响是严格地和所要求的滤波器的频响一致的，逼近误差为零，但在抽样点之间的频响是各取样点的内插函数的延伸叠加而成，有一定的逼近误差，误差大小取决于频率响应曲线的圆滑程度和抽样点的密度为了减少误差就要增加抽样点数目即增大采样频率，抽样点之间的理想频率特性变化越陡，则逼近误差越大，在理想频率特性的不连续点附近会产生肩峰和纹波。频率抽样法的优点是可以直接在频域设计，适于利用最优化方法，而且这种方法特别适用于窄带选频滤波器，但频率抽样法的抽样频率只能是2π/N的整数倍或2π/N的整数倍加上π/N不能保证截