DSP论文dw.讲义

1、DSP原理及应用课程小论文学 院:XXXXXXXXXXXXXXXXXM班级: 电自普 XXXX学生姓名：XX学 号:XXXXXXXXXX论文题目：基于DSP开发板的语言信号滤波处理完成日期：2013年6月10 日指导教师评语：成绩（五级记分制）：指导教师（签字）:.完成日期：2013 年6月10日指导教师评语成绩（五级记分制）指导教师（签字）数字信号处理技术及其应用，目前正以惊人的速度向前发展着。随着大规模 集成电路的出现和数字部件的成本下降、体积缩小及运算速度提高，数字信号处理的应用日益广泛。目前已制成多种专用数字滤波器，取样率可高达兆赫。高速 专用快速傅里叶变换处理机已有商品出售。简单的数

2、字滤波器已制成集成电路片。目前几乎所有的语音带宽压缩系统都倾于全数字化，因为目前它是最实际可 行的方法。除了专用数字信号处理硬件有所发展之外， 还出现了可编程序的数字 信号处理专用计算机，这种计算机的构造特别适于解决数字信号处理问题。它目 前应用于实时信号处理以及设计、模拟专用数字硬件。在计算量相等的情况下IIR数字滤波器比FIR滤波器的幅频特性优越，频率选择性也好，但是它有着致命的缺点相位特性不好控制，它的相位特性是的 非线性函数。例如双线性变换法产生的IIR滤波器模拟指标的频率与数字化指标 的频率转换关系是，这是使频率产生严重的非线性的原因，这种 W与w的非线性关系，使数字滤波器与模拟滤波

3、器在响应与频率的对应关系上发生了畸变， 如果需要线性相位，就必须用全通网络进行复杂的相位校正但是， 在对程序运行 周期数要求十分严格的DSP处理中加上一个全通均衡器是十分浪费资源的，另外即使加上全通均衡器 对于因果的IIR滤波器仍将得不到线性的相位。在现代 电子系统中如图像处理数据传输等波形传递系统中都越来越多的要求信道具有 线性的相位特性在这方面FIR滤波器具有独到的优点，它可以在幅度特性随意设 计的同时保证精确严格的线性相位。本文主要介绍了利用DSP来实现FIR滤波器的设计，根据FIR滤波器设计的要求，本设计以CCS为开发环境，采用模块化的设计方案来完成设计。本文 详细介绍了硬件电路和软件

4、设计的方法，硬件方面介绍了SEED-DEC2812开发实验箱，软件设计方面包括了程序各模块的流程图和源程序代码，还详细介绍了CCS为开发环境，并详细介绍了仿真的步骤。关键词:MATLAB DSP FIR 数字滤波器1.要求和功能实现1.1目的要求通过课程设计，加深对 DSP芯片TMS320C54X的结构、工作原理的理解， 获得DSP应用技术的实际训练，掌握设计较复杂 DSP系统的基本方法。通过使 用汇编语言编写具有完整功能的图形处理程序或信息系统，使学生加深对所学知 识的理解，进一步巩固汇编语言讲法规则。学会编制结构清晰、风格良好、数据 结构适当的汇编语言程序，从而具备解决综合性实际问题的能力

5、。通过使用Matlab设计相应的FIR低通数字滤波器，得到滤波器H(z)的系数, 然后根据这些系数，编写关于DSP的汇编程序。在CCS软件上对已采集信号进行 处理，最后滤除高频信号，输出我们需要的低频语音信号。效果本文采取比较滤 波前后信号的频谱图。1.2功能实现1.2.1设计滤波器。参数指标：根据被处理对象设计滤波器设计指标。1.2.21.2.31.2.41.2.51.2.6用汇编语言在TMSC54X）编程实现符合要求的滤波器。通过开发板上的语音采集模块输入待处理语音信号（简单介绍）。 滤除语音信号中的高频噪声。显示滤波器的频谱。通过开发板上的语音模块播放处理后的语音信号（简单介绍）。2.系

6、统参数和总体设计方案2.1系统参数一个实际的应用系统中，总存在各种干扰。数字滤波器在语音信号处理、 信 号频谱估计、信号去噪、无线通信中的数字变频以及图像信号等各种信号处理中 都有广泛的应用，数字滤波器也是使用最为广泛的信号处理算法之一。在本设计中，使用MATLAB模拟产生合成信号，然后利用 CCS进行滤波。设定模拟信号的采样频率为8000Hz,。设计一个FIR低通滤波器，其参数为:滤波器名称:FIR低通滤波器采样频率:Fs=8000Hz通带/阻带截止频率:2000Hz 2500Hz通带最大衰减:0.5dB阻带最少衰减:50dB滤波器级数:N=25滤波器系数:由MATLAB根据前述参数求得。1

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

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

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

10、入的信号序列，按规定的算法进行处理，从而得到所期望的 输出序列。一个线性位移不变系统的输出序列y(n)和输入x(n)之间的关系，应满 足常系数线性差分方程，见公式4.1，(431 )N AMy(n)=S biX(ni )S ai y(n-i)n>0i =0i 二N是滤波器的(432 )其中，x(n)为输入序列，y(n)为输出序列，ak和bk为滤波器系数, 阶数。若上式中所有的h均为零，则有FIR滤波器的差分方程为：N -4y(n)=2 akx(n-k)kT对上式进行Z变换得到FIR滤波器的传递函数为:(433 )由上式可以看出，H(z)是Z'的N-1次多项式，它在Z平面内有N-1

11、个零点，同时在原点处有N-1个重极点。N阶滤波器通常采用Nt延迟单元、N个加法器与N+1 个乘法器，取图5中(a)、(b)两种结构。A fflI 龍rfJh t 11<0)图5 FIR滤波器的一般结构通常FIR滤波器的设计方法主要有三种：窗函数法、频率抽样法和切比雪夫 等波纹逼近法。其中窗函数法可以应用比较现成的窗函数， 因而设计简单，在指 标要求不高的场合使用方便灵活。431窗函数法窗函数法也称为傅立叶级数法。理想的数字滤波器频率特性是无法实现的, FIR的设计就是要寻找一个可以得到的频率特性 =来逼近，这相当于用一个可 实现的单位脉冲响应h(n)去逼近一个理想单位脉冲响应。可由理想频

12、率特性 通过傅氏反变换得到一般来说，这样得到的理想单位脉冲响应序列 是个无限长序列，因而是非 因果的。设有一个截止频率为 的理想线性相位低通，延时为T，其频率特性是:0<|q7)< 兌£ I时W得到:一工 <n<这是一个以T为中心偶对称的无限长非因果序列，要想用一个有限长的 因果序列去逼近它，最简单的方法是截取n从0到N-1的一段来表示它，即h(n)= hd (n) (0< n < N -1)；其他 N: h(n)=0。同时，为了保证线性相位，还要满足偶对称 h( n)=h(N-1-n)。这就好像通过 一个窗口观看到的一段hd (n)，因此h(n)

13、就表示成hd (n)和一个窗口函数的乘积， 这样对h(n)的求解就变为h(n) = hd (n) * Wn ，这里的Wn就称为窗口函数,既然 一个频域上的标准的矩形窗口对应于时域是一个无限长的序列，那么在时域上截取一段势必造成频域的矩形窗口的失真。结果就是截取出的信号也相应失真， 为了补偿这种失真，只有改变原来窗口的形状，修正经过时域截取后的窗口失真。窗函数设计方法的基本步骤是：把 Hd(ejw 展成 FS,得 hd(n);对hd(n)自然截短到所需的长度，女口 2M+1;将截短后的hd(n)右移M个采样间隔，得h(n);将h(n)乘以合适的窗口，即得所要滤波器的冲击响应，窗函数以n=M寸 称

14、。利用所求得的单位抽样响应，即可用硬件构成滤波器的转移函数 H(z)，也可 利用h(n)在计算机上用软件来实现滤波。4.3.2频率抽样法窗函数法是从时域出发，用窗函数截取理想的 hd(n )得到h(n)，以此有限长 的h(n)近似hd(n )，这样得到的频率响应H(ej)逼近于理想的频响Hd(ejw)。频率 抽样法是从频率出发， 将给定的理想频响Hd(ejw)加以等间隔抽样。FHd(k)Hd(ejw)“然后以此Hd(k )作为FIR滤波器的频率响应抽样值H(k)，再根据DFT(离散付 氏变换)定义由频域这N个抽样值来唯一确定一个有限长序列 h(n)，同样也可以 算出FIR滤波器的系统函数H(z

15、)及频率响应H(e尬)，可以推出频率响应H(e尬是 频率抽样值H(k)与线性相位因子e咬二F及如下内插函数S( o , k)的线性组 合。-盘 smS他心产” 2Stu -U yj所以，在各频率取样点上，实际滤波器的频响是严格地和所要求的滤波器的 频响一致的，逼近误差为零，但在抽样点之间的频响是各取样点的内插函数的延 伸叠加而成，有一定的逼近误差，误差大小取决于频率响应曲线的圆滑程度和抽 样点的密度 为了减少误差 就要增加抽样点数目即增大采样频率，抽样点之间的 理想频率特性变化越陡，则逼近误差越大，在理想频率特性的不连续点附近会产 生肩峰和纹波。频率抽样法的优点是可以直接在频域设计，适于利用最

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

17、 Approaching Way)一致逼近法(Consistent Approaching Way)切比雪夫最佳一致逼近的基本思想是，对于给定区间a，b上的连续函数f(X )，在所有n次多项式的集合錚n中，寻找一个多项式P(X)，使它在a , b 上对f(x )的偏差和其它一切属于鮎的多项式P(X)对f(x)的偏差相比是最小的，即max0(x ) f(X 审=minmax(p(x) f (x )I切比雪夫逼近理论，这样的多项式是存在的，且是唯一的，并指出了构造这 种最佳一致逼近多项式的方法，就是有名的“交错点组定理”。切比雪夫逼近理论解决了P(X)的存在性、唯一性和如何构造等问题。J.H.Mc

18、Clellan、T.W.Parks、L.R.Rabiner 等人应用切比雪夫逼近理论提出了一 种设计FIR滤波器的计算机辅助算法。这种算法由于是在一致意义上对 H d (ejw 作 最佳逼近，因而获得了较好的通带和阻带性能，并能准确地指定通带和阻带的边 缘。但它的效率依赖于初始极值频率点的估计，且通带和阻带内波纹数较多，这 是Chebyshev方法的两个主要缺点。4.4总体设计步骤1、利用MATLAB来确定FIR滤波器的参数；2、启动CCS,在CCS中建立一个C源文件和一个命令文件，并将这两个文 件添加到工程再编译并装载程序；3、设置波形时域观察窗口，得到滤波前后的波形变化图；4、设置频域观察

19、窗口，得到滤波前后的频谱变化图。4.5利用Matlab求出系数bk在matlab的Comma nd Win doW输入以下语句： c=fir1(25,0.5)按回车键后即可出现系数Q5. FIR滤波器的软件设计及其调试5.1根据语音信号的设计参数确定滤波器的系数 利用MATLAB得的系数为c =Colu mns 1 through 90.0078 -0.0119 -0.01770.0014 -0.0018 -0.00290.00480.02560.0368Colu mns 10 through 180.44840.1452 -0.0821-0.0536 -0.08210.14520.4484-

20、0.05360.0368Colu mns 19 through 260.0256 -0.0177 -0.01190.00780.0048 -0.0029 -0.00180.001445.2汇编源程序采用循环缓冲区实现FIR滤波器的源程序如下:.title "FIR.ASM".mmregs.def start.bss yn,1xnhnin data .usect "in data",1024 outdata .usect "outdata",1024.usect "xn ",25 .usect "hn &q

21、uot;,25.data fir_coff: .word .word .word .word .word .word .word .word .word .word .word .word .word14*32768/10000,-18*32768/10000 -29*32768/10000,48*32768/10000 78*32768/10000,-119*32768/10000 -177*32768/10000,256*32768/10000 368*32768/10000,-536*32768/10000 -821*32768/10000,1452*32768/10000 4484*3

22、2768/10000,4484*32768/10000 1452*32768/10000,-821*32768/10000 -536*32768/10000,368*32768/10000 256*32768/10000,-177*32768/10000 -119*32768/10000,78*32768/10000 48*32768/10000,-29*32768/10000 -18*32768/10000,14*32768/10000.textstart:ssbx frct;stm #hn, ar1 ; rpt #39;mvpd fir_coff,*ar1+ ; stm stm stm s

23、tm rpt使用小数运算系数首地址 将系数移入循环缓冲区#-1,arO #outdata,ar5 #i ndata,ar4 #xn+39,ar2 #39;mvdd *ar4+,*ar2+0% ;将输入数据移入 循环缓冲区stmstmstmstmstm rptb rptz mac sth#xn+39,ar2 ; #hn+39,ar3 #40,bk;#1024,brc;#-1,ar0loo p-1;a,#39;*ar2+0%,*ar3+0%,aa,*(y n);指向第一个输入缓冲区大小40块重复1024次块程序重复大小 计算一个输出保存输出mvkd *(y n) ,*ar5+ mvdd *ar4+

24、,*ar2+0% ; loop: nop b loop.end5.3汇编源程序的连接命令文件在TMS320VC5402DS的硬件系统中，用户可使用的程序存储空间：片内0080H3FFF H片外 48000H4FFFFH用户可使用的数据存储空间：片内 0080H3FFFFH片外 4000H7FFFH 基于 TMS320VC540的资源配置，FIR 滤波器 源程序的连接命令文件如下：MEMORY P AGE 0:PARAM:P AGE 1:SP RAMDARAM:读进一个输入org=1000horg=2060horg=2080hlen=1000hlen=0020h len=9f380h SECTI

25、ONS.text :>.data :>.bss :>P ARAM P ARAM SP RAMxn : alig n (64) >hn : alig n (64) >P AGE 0P AGE 0P AGE 1DARAM PAGE 1DARAM PAGE 1P AGE 1P AGE 1in data :> DARAM outdata :> DARAM6 .FIR滤波器的检验6.1 利用 matlab 导出 “ indata.dat ”文件在matlab中输入下面的程序便可以导出“ indata.dat ”文件:f12=1000; %/Hzf13=16000

26、; %/Hz fs=40000; %/ 采样 HzN=1000;% 数据个数T=1/fs; %采样周期n=0:N;x12=0.7*si n(2* pi*f12* n*T); x13=0.5*si n(2* pi*f13* n*T); x_base=( x12+x13);%寺滤波信号波形figure(1)pl ot(x_base)%寺滤波信号频谱归一化打开文件输出文件头输出figure (2) yff=abs(fft(x_base) df=n *(fs/N) plot(df,yff) xout=x_base/max(x_base);% xto_ccs=ro un d(32767*xout) fi

27、d=fo pen( 'i np ut.dat','w');% fprin tf(fid,'1651 1 0 0 0n');% fprin tf(fid,'%dn',xto_ccs);% fclose(fid);6. 硬件部分简介DSP的开发系统包括硬件和软件两个方面。其中硬件包括一套DSK开发板和仿真器。DSP的开发主要是围绕着这两方面进行的。6.1 DSK简介DSK的全称是DSP Starter Kit ，即DSP初学者套件。DSK为用户提供了一个开发DSP的方便的平台。由于DSK板上提供了基本的硬件支持，结合适当的软件开发工具

28、，简化了开发工作，缩短了开发周期，在DSK板上主要包括如下资源：100 MHz VC5402 DSP芯片64K字外部扩展的SRAM 64K*16256K字的 FLASH存储器 256K*16 6.2仿真器仿真器JTAG测试总线控制器SN74ACT8990 JTAG TB及与主机相连接的并行接口模/数转换器TI TLC320AD50 A/D converter 两个电话接口 (DAA)麦克风/扬声器接口并行口RS-232 串行口扩展子板接口LED*4其中三个可供程序运行时点亮提示DSK板上提供的丰富的硬件资源，使得开发过程中硬件方面可以节省大量的 工作。DSK板上提供了一个八位跳线开关，可以通过

29、它申请外部中断，实现不同 的功能。另外，板上还提供了诸如14针的外部JTAG接头，10针的CPLD勺JTAG 接头，25针的符合IEEE-1284的为JTAG/HPI接入提供的并口等外部接口，通过 JTAG接头可以与计算机内运行的 CCS Code Comp oser Studio进行实时数据交换 RTDX提高编程调试的效率通过JTAG接头与主机进行实时数据交换的过程可见图14示。图14实时数据交换的过程JTAG接口符合IEEE 1149.1标准边界扫描逻辑测试技术规范，这一技术规 范是80年代由联合测试行动组 JTAG Joi nt Test Actio n Group开发的，这个边界扫描测

30、试BST结构提供了有效的测试，弓I线间隔致密的电路板上零部件的能 力，通过BST结构测试引脚连结而不必使用物理测试探针，可以在器件正常工 作时，获功能数据。器件的边界扫描单元能够迫使逻辑追踪引脚信号或是从引脚 或器件核心逻辑信号中捕获数据，强行加入的测试数据串行移入边界扫描单元 捕获的数据串行移出并在器件外部同预期的结果进行比较。通过扩展可以将多个器件连成菊花链，一次进行统一的输入与读出操作。TI提供的增强型JTAG连接，是一种可与任意DSP系统相连的低侵扰式的连 接，使得CCS能够控制程序的执行，实时监视程序运行。通过主机和DSP APIs提供主机和DSP之间的双向实时数据交换，能够使开发者实时连续的观察到 DSP 应用的实际工作方式，仿真接口提供主机一侧的JTAG连接，如TI XSD510。7. 遇到问题和解决方法在调试过程中曾遇到许多的问题，具体总结如下：(1) 采集的语音信号中高频噪音不很明显， 使得最后不能较好的观察到滤波器的性能。解决的方法是：利用matlab产生一个6khz的高频信号与原来的语音信号相 加，这时可以明显的听到语音信号中同时伴随着刺耳的高频噪音。(2) 在选择滤波器的参数时，对具体的音频信息没有整体的了解，使得最后 FIR滤波器不能达到很好的滤波效果。解决的方法是，首先