DSP课程设计论文.docx

1、DSP课程设计论文DSP小论文基于DSP开发板的语言信号滤波处理系 、 部：电气与信息工程系 学生姓名： 学 号： 指导教师： 职称 讲师 专 业： 通信工程 班 级： 完成时间： 2014-1-1 摘 要数字信号处理技术及其应用，目前正以惊人的速度向前发展着。随着大规模集成电路的出现和数字部件的成本下降、体积缩小及运算速度提高，数字信号处理的应用日益广泛。目前已制成多种专用数字滤波器，取样率可高达兆赫。高速专用快速傅里叶变换处理机已有商品出售。简单的数字滤波器已制成集成电路片。目前几乎所有的语音带宽压缩系统都倾于全数字化，因为目前它是最实际可行的方法。除了专用数字信号处理硬件有所发展之外，还

2、出现了可编程序的数字信号处理专用计算机，这种计算机的构造特别适于解决数字信号处理问题。它目前应用于实时信号处理以及设计、模拟专用数字硬件。在计算量相等的情况下IIR 数字滤波器比FIR滤波器的幅频特性优越，频率选择性也好，但是它有着致命的缺点相位特性不好控制，它的相位特性 是的非线性函数。例如双线性变换法产生的IIR 滤波器模拟指标的频率与数字化指标的频率转换关系是 ，这是使频率产生严重的非线性的原因，这种W与w的非线性关系，使数字滤波器与模拟滤波器在响应与频率的对应关系上发生了畸变，如果需要线性相位，就必须用全通网络进行复杂的相位校正但是，在对程序运行周期数要求十分严格的DSP处理中加上一个

3、全通均衡器是十分浪费资源的，另外即使加上全通均衡器 对于因果的IIR 滤波器 仍将得不到线性的相位。在现代电子系统中如图像处理数据传输等波形传递系统中都越来越多的要求信道具有线性的相位特性在这方面FIR滤波器具有独到的优点，它可以在幅度特性随意设计的同时保证精确严格的线性相位。本文主要介绍了利用DSP来实现FIR滤波器的设计，根据FIR滤波器设计的要求，本设计以CCS为开发环境，采用模块化的设计方案来完成设计。本文详细介绍了硬件电路和软件设计的方法，硬件方面介绍了SEED-DEC2812开发实验箱，软件设计方面包括了程序各模块的流程图和源程序代码，还详细介绍了CCS为开发环境，并详细介绍了仿真

4、的步骤。关键词:MATLAB DSP FIR 数字滤波器目录摘要 11.课程设计的目的要求和功能实现 31.1目的要求 31.2功能实现 32.系统参数和总体设计方案 32.1系统参数 32.2系统总体设计方案流程图 43. DSP及其开发环境 43.1 DSP系统的构成 43.2DSP系统的特点及设计过程 43.3 CCS开发环境 54.数字滤波器的设计原理和方案对比 54.1 IIR滤波器 54.2 FIR滤波器 64.3 FIR滤波器的基本结构 64.4总体设计步骤 94.5利用Matlab求出系数 95.FIR滤波器的软件设计及其调试 95.1根据语音信号的设计参数确定滤波器的系数 9

5、5.2 汇编源程序 105.3 汇编源程序的连接命令文件 116 .FIR滤波器的检验 116.1利用matlab导出“indata.dat”文件 116.2 FIR滤波器检验的验证过程 127.硬件部分简介 168.遇到问题和解决方法 189.心得体会 1810.参考文献 191.课程设计的目的要求和功能实现1.1目的要求通过课程设计，加深对DSP芯片TMS320C54x的结构、工作原理的理解，获得DSP应用技术的实际训练，掌握设计较复杂DSP系统的基本方法。通过使用汇编语言编写具有完整功能的图形处理程序或信息系统，使学生加深对所学知识的理解，进一步巩固汇编语言讲法规则。学会编制结构清晰、风

6、格良好、数据结构适当的汇编语言程序，从而具备解决综合性实际问题的能力。通过使用Matlab设计相应的FIR低通数字滤波器，得到滤波器H(z)的系数，然后根据这些系数，编写关于DSP的汇编程序。在CCS软件上对已采集信号进行处理，最后滤除高频信号，输出我们需要的低频语音信号。效果本文采取比较滤波前后信号的频谱图。1.2功能实现1.2.1设计滤波器。参数指标：根据被处理对象设计滤波器设计指标。1.2.2用汇编语言在TMSc54xx编程实现符合要求的滤波器。1.2.3通过开发板上的语音采集模块输入待处理语音信号(简单介绍)。1.2.4滤除语音信号中的高频噪声。1.2.5显示滤波器的频谱。1.2.6通

7、过开发板上的语音模块播放处理后的语音信号（简单介绍）。2.系统参数和总体设计方案2.1系统参数一个实际的应用系统中，总存在各种干扰。数字滤波器在语音信号处理、信号频谱估计、信号去噪、无线通信中的数字变频以及图像信号等各种信号处理中都有广泛的应用，数字滤波器也是使用最为广泛的信号处理算法之一。在本设计中，使用MATLAB模拟产生合成信号，然后利用CCS进行滤波。设定模拟信号的采样频率为8000Hz，。设计一个FIR低通滤波器，其参数为：滤波器名称： FIR低通滤波器采样频率： Fs=8000Hz通带/阻带截止频率： 2000Hz2500Hz通带最大衰减： 0.5dB阻带最少衰减： 50dB滤波器

8、级数： N=25滤波器系数： 由MATLAB根据前述参数求得。2.2系统总体设计方案流程图3. DSP及其开发环境3.1 DSP系统的构成一个典型的 DSP 系统如图2示。图2 典型的DSP系统图2是一个用DSP做信号处理的典型框图。由于DSP是用来对数字信号进行处理的，所以首先必须将输入的模拟信号变换为数字信号。3.2 DSP系统的特点及设计过程由于数字信号处理系统是以数字信号处理理论为基础，所以具有数字信号处理的全部优点：接口方便；编程方便；具有高速性；稳定性好；精度高；可重复性好；集成方便等等，同时也存在一定的缺点 一般来说DSP的设计过程应遵循一定的设计流程，如图3示。图3 DSP基本

9、设计流程3.3 CCS开发环境CCS 全称是 Code Composer Studio 它提供了基本的代码生成工具，具有一定的调试、分析能力，在CCS下的程序开发过程如图4示。图4 程序开发过程CCS 包括：1、CCS代码生成工具2、CCS集成开发环境IDE3、DSP/BIOS插件程序和API RTDX插件，主机接口和API等。4. 数字滤波器的设计原理和方案对比4.1 IIR滤波器无限长冲激响应系统IIR（Infinite Impulse Response）中，用有理分式表示的系统函数来逼近所需要的频率响应，即其单位冲激响应h(n)是无限长的。IIR滤波器由于吸收了模拟滤波器的结果，有大量的

10、图表可查，可以方便、简单、有效地完成设计，效果很好，但是其相位特性不好控制，必须用全通网络进行复杂的相位较正，才能实现线性相位特性的要求。4.2 FIR滤波器有限长冲激响应系统FIR（Finite Impulse Response）中，用一个有理多项式表示的系统函数去逼近所需要的频率响应，即其单位冲激响应h(n)在有限个n值处不为零。FIR滤波器可在幅度特性随意设计的同时，保证精确、严格的线性相位特性。这在要求相位线性信道的现代电子系统，如图像处理、数据传输等波形传递系统中，是具有很大吸引力的。而且，其单位冲激响应是有限长的，不存在不稳定的因素，并且可用因果系统来实现。4.3 FIR滤波器的基

11、本结构数字滤波是将输入的信号序列，按规定的算法进行处理，从而得到所期望的输出序列。一个线性位移不变系统的输出序列y(n)和输入x(n)之间的关系，应满足常系数线性差分方程，见公式4.1， （4.3.1）其中，x(n)为输入序列，y(n)为输出序列，为滤波器系数，N是滤波器的阶数。若上式中所有的均为零，则有FIR滤波器的差分方程为： （4.3.2）对上式进行Z变换得到FIR滤波器的传递函数为： （4.3.3）由上式可以看出，H(z)是的N-1次多项式，它在z平面内有N-1个零点，同时在原点处有N-1个重极点。N阶滤波器通常采用N个延迟单元、N个加法器与N+1个乘法器，取图5 中(a)、(b)两种

12、结构。图5 FIR滤波器的一般结构通常FIR滤波器的设计方法主要有三种：窗函数法、频率抽样法和切比雪夫等波纹逼近法。其中窗函数法可以应用比较现成的窗函数，因而设计简单，在指标要求不高的场合使用方便灵活。4.3.1 窗函数法窗函数法也称为傅立叶级数法。理想的数字滤波器频率特性 是无法实现的， FIR的设计就是要寻找一个可以得到的频率特性 = 来逼近 ，这相当于用一个可实现的单位脉冲响应 h(n)去逼近一个理想单位脉冲响应 。 可由理想频率特性 通过傅氏反变换得到 一般来说，这样得到的理想单位脉冲响应序列 是个无限长序列，因而是非因果的。设有一个截止频率为 的理想线性相位低通，延时为，其频率特性是

13、： 得到： 这是一个以为中心偶对称的无限长非因果序列，要想用一个有限长的因果序列去逼近它，最简单的方法是截取从0到N-1的一段来表示它，即 h(n)= ；其他N：h(n)=0。同时，为了保证线性相位，还要满足偶对称h(n)=h(N-1-n)。这就好像通过一个窗口观看到的一段，因此h(n)就表示成和一个窗口函数的乘积，这样对h(n)的求解就变为h(n)* ，这里的就称为窗口函数, 既然一个频域上的标准的矩形窗口对应于时域是一个无限长的序列, 那么在时域上截取一段势必造成频域的矩形窗口的失真。结果就是截取出的信号也相应失真，为了补偿这种失真，只有改变原来窗口的形状，修正经过时域截取后的窗口失真。窗

14、函数设计方法的基本步骤是：(1) 把展成FS，得；(2) 对自然截短到所需的长度，如2M+1；(3) 将截短后的右移M个采样间隔，得h(n)；(4) 将h(n)乘以合适的窗口，即得所要滤波器的冲击响应，窗函数以n=M对称。利用所求得的单位抽样响应，即可用硬件构成滤波器的转移函数H(z)，也可利用h(n)在计算机上用软件来实现滤波。4.3.2 频率抽样法窗函数法是从时域出发，用窗函数截取理想的得到h(n)，以此有限长的h(n)近似，这样得到的频率响应逼近于理想的频响。频率抽样法是从频率出发， 将给定的理想频响加以等间隔抽样。 然后以此作为FIR滤波器的频率响应抽样值H(k)，再根据DFT（离散付

15、氏变换）定义由频域这N个抽样值来唯一确定一个有限长序列h(n)，同样也可以算出FIR滤波器的系统函数H(z)及频率响应，可以推出频率响应是频率抽样值H(k)与线性相位因子及如下内插函数S(, k) 的线性组合。所以，在各频率取样点上，实际滤波器的频响是严格地和所要求的滤波器的频响一致的，逼近误差为零，但在抽样点之间的频响是各取样点的内插函数的延伸叠加而成，有一定的逼近误差，误差大小取决于频率响应曲线的圆滑程度和抽样点的密度 为了减少误差 就要增加抽样点数目即增大采样频率，抽样点之间的理想频率特性变化越陡，则逼近误差越大，在理想频率特性的不连续点附近会产生肩峰和纹波。频率抽样法的优点是可以直接在

16、频域设计，适于利用最优化方法，而且这种方法特别适用于窄带选频滤波器， 但频率抽样法的抽样频率只能是2/ N 的整数倍或2/ N 的整数倍加上/ N不能保证截止频率c的准确取值，要实现精确的c就必须取N大，相应的计算量也大。此外，它的阻带最大衰减一般，也只有30-50dB左右，很难满足频域特性要求较高的场合。 4.3.3 Chebyshev逼近法窗函数法和频率采样法设计出的滤波器的频率特性都是在不同意义上对所给理想频率特性的逼近。由数值逼近理论可知，对某个函数f(x)的逼近一般有以下三种方法：插值法(Interpolating Way)最小平方逼近法(Least Square Approachi

17、ng Way)一致逼近法(Consistent Approaching Way)切比雪夫最佳一致逼近的基本思想是，对于给定区间a，b上的连续函数，在所有n次多项式的集合中，寻找一个多项式 p (x)，使它在a，b上对的偏差和其它一切属于的多项式 p(x)对f(x)的偏差相比是最小的，即 切比雪夫逼近理论，这样的多项式是存在的，且是唯一的，并指出了构造这种最佳一致逼近多项式的方法，就是有名的“交错点组定理”。切比雪夫逼近理论解决了p(x)的存在性、唯一性和如何构造等问题。J.H.McClellan、T.W.Parks、L.R.Rabiner 等人应用切比雪夫逼近理论提出了一种设计FIR滤波器的计

18、算机辅助算法。这种算法由于是在一致意义上对作最佳逼近，因而获得了较好的通带和阻带性能，并能准确地指定通带和阻带的边缘。但它的效率依赖于初始极值频率点的估计，且通带和阻带内波纹数较多，这是Chebyshev方法的两个主要缺点。4.4总体设计步骤1、 利用MATLAB来确定FIR滤波器的参数；2、 启动CCS，在CCS中建立一个C源文件和一个命令文件，并将这两个文件添加到工程再编译并装载程序；3、设置波形时域观察窗口，得到滤波前后的波形变化图；4、设置频域观察窗口，得到滤波前后的频谱变化图。4.5利用Matlab求出系数在matlab的Command Window输入以下语句：c=fir1(25,

19、0.5)按回车键后即可出现系数5.FIR滤波器的软件设计及其调试5.1根据语音信号的设计参数确定滤波器的系数利用MATLAB求得的系数为c = Columns 1 through 9 0.0014 -0.0018 -0.0029 0.0048 0.0078 -0.0119 -0.0177 0.0256 0.0368 Columns 10 through 18 -0.0536 -0.0821 0.1452 0.4484 0.4484 0.1452 -0.0821 -0.0536 0.0368Columns 19 through 26 0.0256 -0.0177 -0.0119 0.0078 0

20、.0048 -0.0029 -0.0018 0.001445.2 汇编源程序采用循环缓冲区实现FIR滤波器的源程序如下： .title FIR.ASM .mmregs .def start .bss yn,1 xn .usect xn,25hn .usect hn,25indata .usect indata,1024outdata .usect outdata,1024 .datafir_coff: .word 14*32768/10000,-18*32768/10000 .word -29*32768/10000,48*32768/10000 .word 78*32768/10000,-1

21、19*32768/10000 .word -177*32768/10000,256*32768/10000 .word 368*32768/10000,-536*32768/10000 .word -821*32768/10000,1452*32768/10000 .word 4484*32768/10000,4484*32768/10000 .word 1452*32768/10000,-821*32768/10000 .word -536*32768/10000,368*32768/10000 .word 256*32768/10000,-177*32768/10000 .word -11

22、9*32768/10000,78*32768/10000 .word 48*32768/10000,-29*32768/10000 .word -18*32768/10000,14*32768/10000 .textstart: ssbx frct ;使用小数运算 stm #hn,ar1 ;系数首地址 rpt #39 ;将系数移入循 mvpd fir_coff,*ar1+ ;环缓冲区 stm #-1,ar0 stm #outdata,ar5 stm #indata,ar4 stm #xn+39,ar2 rpt #39 ;将输入数据移入 mvdd *ar4+,*ar2+0% ;循环缓冲区 stm

23、 #xn+39,ar2 ;指向第一个输入 stm #hn+39,ar3 stm #40,bk ;缓冲区大小40 stm #1024,brc ;块重复1024次 stm #-1,ar0 rptb loop-1 ;块程序重复大小 rptz a,#39 ;计算一个输出 mac *ar2+0%,*ar3+0%,a sth a,*(yn) ;保存输出 mvkd *(yn),*ar5+ mvdd *ar4+,*ar2+0% ;读进一个输入loop: nop b loop .end 5.3 汇编源程序的连接命令文件 在TMS320VC5402DSP的硬件系统中，用户可使用的程序存储空间：片内0080H3FF

24、FH,片外48000H4FFFFH;用户可使用的数据存储空间：片内0080H3FFFFH,片外4000H7FFFH。基于TMS320VC5402的资源配置，FIR滤波器源程序的连接命令文件如下：MEMORY PAGE 0: PARAM: org=1000h len=1000h PAGE 1: SPRAM org=2060h len=0020h DARAM: org=2080h len=9f380h SECTIONS .text : PARAM PAGE 0 .data : PARAM PAGE 0 .bss : SPRAM PAGE 1 xn : align (64) DARAM PAGE 1

25、 hn : align (64) DARAM PAGE 1 indata : DARAM PAGE 1 outdata : DARAM PAGE 1 6 .FIR滤波器的检验 6.1 利用matlab导出“indata.dat”文件 在matlab中输入下面的程序便可以导出“indata.dat”文件：f12=1000; %/Hz f13=16000; %/Hz fs=40000; %/采样 Hz N=1000; %数据个数T=1/fs; %采样周期n=0:N; x12=0.7*sin(2*pi*f12*n*T); x13=0.5*sin(2*pi*f13*n*T); x_base=( x12

26、+x13); %待滤波信号波形figure(1) plot(x_base) %待滤波信号频谱figure(2) yff=abs(fft(x_base) df=n*(fs/N) plot(df,yff) xout=x_base/max(x_base);%归一化xto_ccs=round(32767*xout) fid=fopen(input.dat,w);%打开文件fprintf(fid,1651 1 0 0 0n);%输出文件头fprintf(fid,%dn,xto_ccs);%输出fclose(fid);6.2 FIR滤波器检验的验证过程1）进入CCS，编译原文件，生成*.out文件2）选择

27、FileLoad Program打开DlabDebugfir.out3) 选择FileDataLoad，在弹出的对话框中选择本实验文件夹下的indata.dat 文件并打开，在弹出的对话框中输入起始地址indata和长度10244)输入完后点击OK，然后选择Debug Run 5)选择ViewGraphTime/Frequency，通过设置Graph Property Dialog可以查看输入与输出的时域与频域波形,设置的参数和波形截图如下所示：滤波前的语音信号的时域Graph Property Dialog设置及时域波形图分别如图6、7所示：图6 滤波前语音信号的时域Graph Proper

28、ty Dialog设置参数图 图7 滤波前语音信号的时域波形图滤波前的语音信号的频域Graph Property Dialog设置及频域波形图分别如图8、9所示： 图8 滤波前语音信号的频域Graph Property Dialog设置参数图图9 滤波前语音信号的频域波形图滤波后的语音信号的时域Graph Property Dialog设置及时域波形图分别如图10、11所示：图10 滤波后语音信号的时域Graph Property Dialog设置参数图图11 滤波后语音信号的时域波形图滤波后的语音信号的频域Graph Property Dialog设置及频域波形图分别如图12、13所示：图

29、12 滤波后语音信号的频域Graph Property Dialog设置参数图图13 滤波后语音信号的频域波形图结论：比较滤波前和滤波后的语音信号的频域的波形图可以明显看到高频的噪声已经被滤除掉了，证明本组设计的数字FIR滤波器是合理的，它能对所采集的语音信号进行高效滤波，得到稳定的语音信号，滤波效果良好。7.硬件部分简介DSP 的开发系统包括硬件和软件两个方面。其中硬件包括一套DSK开发板和仿真器。 DSP的开发主要是围绕着这两方面进行的。7.1 DSK简介DSK 的全称是DSP Starter Kit，即DSP初学者套件。DSK为用户提供了一个开发 DSP 的方便的平台。由于DSK板上提供

30、了基本的硬件支持，结合适当的软件开发工具，简化了开发工作，缩短了开发周期，在 DSK 板上主要包括如下资源：100 MHz VC5402 DSP芯片64K字外部扩展的SRAM 64K*16256K字的FLASH存储器256K*167.2 仿真器仿真器JTAG测试总线控制器SN74ACT8990 JTAG TBC 及与主机相连接的并行接口模/数转换器 TI TLC320AD50 A/D converter两个电话接口(DAA)麦克风/扬声器接口并行口RS-232串行口扩展子板接口LED*4 其中三个可供程序运行时点亮提示 DSK板上提供的丰富的硬件资源，使得开发过程中硬件方面可以节省大量的工作。DSK板上提供了一个八位跳线开关，可以通过它申请外部中断，实现不同的功能。另外，板上还提供了诸如14 针的外部 JTAG 接头，10针的CPLD的JTAG接头，25针的符合IEEE-1284的为JTAG/HPI接入提供的并口等外部接口，通过JTAG接头可以与计算机内运行的CCS Code Composer Studio进行实时数据交换RTDX提高编程调试的效率 通过JTAG接头与主机进行实时数据交换的过程可见图14