DSP课程设计 基于DSP开发板的语言信号滤波处理 电信毕业设计.doc

dsp原理与接口技术课程设计报告基于dsp开发板的语言信号滤波处理姓 名： 班 级： 09电信1 学 号： 2009305901 2009305901 2009305901 指导老师： 日期： 2012. xx.xx2012.xx.xx 华南农业大学工程学院 摘 要 通过语音传递信息是人类最重要、最有效、最常用、最方便的交换信息的方法，因此，如何去除混杂在有用信号中的噪声并获得低损甚至无损的有用语音信号成为语音信号处理（即用数字信号处理技术和语言学知识对语音信号进行处理并提取有用信息）至关重要的问题。但是因为语音识别技术对信号噪声非常敏感，系统本身和器件对系统的处理结果都有着不可忽视的干扰，因此很难实现效果较好的语音信号滤波处理。而dsp因为运算速度快，片上资源丰富和能够实现复杂的线性和非线性算法等特性，在语音信号处理技术方面有明显的优势。虽然dsp自身在一定程度上也是语音信号干扰源，但综其上述良好特性，以及干扰源的不可避免性，我们选择用dsp实现语音信号滤波处理。 为了实现基于dsp开发板的语音信号滤波处理，我们对混杂噪声信号的有用信号音频进行时域和频谱分析，得出其频谱特性后，在软件编程方面，我们用汇编语言在ccs中设计出三个滤波器，分别是：低通滤波器、高通滤波器以及带通滤波器，同时实现a/d和d/a转换。在硬件方面，通过接口电路从电脑下载混杂噪声的有用音频，以及将程序下载到芯片中。语音信号在硬件电路芯片上进行adc转换变为数字信号，通过滤波、压缩等变化滤去噪声之后再进行da转换返回模拟语音信号并从硬件电路输出。从而实现dsp的语音信号滤波处理。 关键词: 语音信号处理 滤波 dsp39目 录1 方案比较与选择1.1 方案目标 基于dsp开发板，在ccs中用汇编语言设计一个滤波器，运用该滤波器能将混杂有高频噪声的歌声中的噪声信号去除。滤波的整个过程通过dsp开发板上实现，并且要满足以下要求： 1、根据被处理对象设置滤波器参数；2、该滤波器要用汇编语言实现；3、待处理语音信号采集通过开发板上的信号采集模块输入；4、所设计的滤波器能滤除语音信号中的高频噪声；5、实现滤除噪声功能后，能在电脑上显示滤波器的频谱；6、通过开发板上的语音模块播放处理后的语音信号；我们组将合成信号（混有高频噪声与有用的歌曲信号）通过dsp开发板的耳机输入接口输入tms320vc5416芯片，该音频信号的时域波形和频域波形如图1、图2、图3、图4所示：图1 混合时域波形（1）图2 混合时域波形（2）图3 混合时域波形（3）图4 混合频域波形 通过滤波器对高频噪声的滤除后，理想情况下的有用信号输出时域波形和频域波形如图5、图6所示：图5 理想结果时域图6 理想结果频域而被滤除的高频噪声信号的时域波形和频域波形则分别如图7、图8所示：图7 噪声时域图8 噪声频域1.2方案选择1.2.1 滤波器选择a.iir数字滤波器iir数字滤波器具有无限的脉冲响应，具有如下特点：1. iir数字滤波器的系统函数可以写成封闭函数的形式；2. iir数字滤波器采用递归型结构；3. iir数字滤波器在设计上可借助成熟的模拟滤波器的结果；4. iir数字滤波器需要加相位校准网络； iir数字滤波器的相位特性不好控制，当对相位要求较高的时候，需要外加相位校准网络。由于iir数字滤波器采用递归型结构，结构上带有反馈回路，在运算的舍入处理过程中，会使误差不断累积，致使有时会产生微弱的寄生振荡。但是在设计上，iir数字滤波器能借助成熟的模拟滤波器结果，如巴特沃斯、契比雪夫等，因有现成的设计数据或图表可以查阅，使设计工作量相对比较小，对计算工具的要求也不高。b.fir滤波器fir滤波器为有限长单位冲激响应滤波器，是数字信号处理系统中最基本的原件，它能够保证任意幅频特性的同时又具有严格的线性相频特性，与此同时，它的单位抽样响应是有限长的，因此，fir滤波器是稳定的系统。它具有如下特点：1. fir滤波器系统的单位冲激响应h(n)在有限个n值处不为零；2. fir滤波器的系统函数h（z）在|z|0处收敛，极点全部在z=0处，是因果系统；3. fir滤波器在结构上主要为非递归结构，没有输出到输入的反馈；fir滤波器能够在dsp芯片中实现，dsp芯片中有专用的数字信号处理函数可以调用，因此，在dsp芯片中，fir滤波器的实现相对简单。但是由于程序是顺序执行的，它的速度会受到限制。1.2.2 dsp芯片选择 a.tms320c5402 ti公司第五代16bit定点dsp处理器，它适用于语音通信等实施嵌入应用场合，其特性为：操作速度可以达100mips；先进的多总线结构；整合维特比加速器；数据/程序寻址从、空间为1m16bit，内置4k16bit的rom和16k16bit的ram；内置可编程等待状态发生器、锁相环（pll）时钟产生器、两个多通道缓冲串口、两个16位定时器，与c54x系列的其他芯片相比，c5402具有高性能，低功耗和低价格等特点，它采用6级流水线，而且当rpt（重复指令）有效时，一些多周期的指令就编程单周期的指令；芯片内部ram和rom可以根据pmst寄存器中的ovly和drom位灵活设置，这些都是有利于算法的优化。 b.tms320vc5416 tms320vc5416芯片的运行速度为160m，具有低功耗设计；其工作速度可以高达160mips；它的片内ram存储空间为128k16bit；而其扩张的2路6bita/d接口tlv0832的最大采样速率为20k；具有2路的tlc7528d/a转换，10m/s,8bit；具有uaet串行接口，符合rs232标准；有8m bit的扩展flash，能存储大量固化程序和数据；设计配置有用户可以自定义的开关盒测试指示灯，方便用户进行调试；4组标准化扩展连接器，为用户进行二次开发提供条件；具有能与ieee1149.1相兼容的仅用于测试和仿真的逻辑扫描电路；+5v电源输入，内部为+3.3v、+1.6v电源管理；四层板设计工艺，工作稳定可靠；具有自启动功能设计，可以实现脱机工作；可以选配多张应用接口板，包括usb板。1.3方案选择1.3.1 滤波器选择在很多实际应用中如语音和音频信号处理中，常用数字滤波器来实现选频功能。因此，应将频域中的幅度和相位响应设为选频功能的指标。在通带中，线性相位响应常是被希望的，但精确的线性相位在fir滤波器中可以得到，而在iir滤波器中通带的相位却不可能得到的。在幅度指标的考量方面，iir数字滤波器的设计和模拟滤波器的设计有着紧密的联系，通常是在设计出适当的模拟滤波器后，再通过一定的频带变换把它转换成为所需要的数字iir滤波器。此外，任何数字信号处理系统中也还不可避免地用到模拟滤波器，比如a/d变换器前的抗混叠滤波器及d/a转换后的平缓滤波器，因此模拟滤波器设计也是数字信号处理中应当掌握的技术。从性能上来说，iir数字滤波器传递函数包括零点和极点两组可调因素，对极点的唯一限制是在单位圆内。因此可用较低的阶数获得高的选择性，所用的存储单元少、计算量小、效率高。但是这个高效率是以相位的非线性为代价的。iir数字滤波器的选择性越好，则相位非线性越严重。因为在计算量相等的情况下，虽然iir 数字滤波器比fir滤波器的幅频特性优越，频率选择性也好，但是它有着致命的缺点：相位特性不好控制。它的相位特性是非线性函数。例如双线性变换法产生的iir 滤波器模拟指标的频率与数字化指标的频率转换关系是 ，这是使频率产生严重的非线性的原因，这种w与w的非线性关系，使数字滤波器与模拟滤波器在响应与频率的对应关系上发生了畸变。而fir滤波器传递函数的极点是固定在原点，是不能动的，它只能靠改变零点位置来改变它的性能，所以要达到高的选择性，必须用高的阶数，对于同样的滤波器设计指标，fir滤波器所要求的阶数可能比iir滤波器高5-10倍，结果成本高信号延时也较大，如果按线性相位要求来说，则iir滤波器就必须加全通网络进行相位校正，但是，在对程序运行周期数要求十分严格的dsp处理中加上一个全通均衡器是十分浪费资源的，另外即使加上全通均衡器 对于因果的iir 滤波器仍将得不到线性的相位，反而会大大增加了滤波器的阶数和复杂性。而fir滤波器却可以在幅度特性随意设计的同时保证精确严格的线性相位。fir滤波器的这一特点在现代电子系统中诸如图像处理数据传输等波形传递系统中都越来越多的要求信道具有线性相位特性的情况下，成为了fir滤波器独特的优点。另外，从结构上看，iir滤波器必须采用递归结构来配置极点，并保证极点位置在单位圆内。由于有限字长效应，运算过程中将对系数进行舍入处理，引起极点的偏移，这种情况有时会造成稳定性问题，甚至造成寄生振荡。相反，fir滤波器只要采用非递归结构，不论在理论上还是实际的有限精度运算中都不存在稳定性问题，因此造成的频率特性误差也较小。此外fir滤波器可以采用快速傅立叶变换算法，在相同的阶数条件下运算速度可以快的多。因此，我们小组选择fir滤波器进行后续的语音信号滤波处理。1.3.2 dsp芯片选择 综合分析tms320c5402芯片和tms320vc5416的功能特点，在同样具有低功耗、高稳定性的条件下，tms320c5402芯片具有先进的多总线结构，能整合维特比加速器，同时它采用6级流水线，使工作能同时进行，工作速率提高，其芯片内部ram和rom可以根据pmst寄存器中的ovly和drom位灵活设置等有利于算法的优化，并且价格较低，适合于学生进行开发。但是，tms320vc5416芯片具有更高的运行速度，是目前为止，c54x系列产品中最快的，其速度是c5410的1.6倍，其内部空间更大，是c5410的2倍。同时配置有较大的rom能存储大量的固化程序和数据，有利于运行速度的提升和减少因处理数据量太大而死机的可能，能处理信息量的音频信号。tms320vc5416芯片配置有用户可以自定义的开关盒测试指示灯和4组标准化扩展连接器，为用户调试和进行二次开发提供便利条件，此外，它能够选配诸如usb接口的应用接口板，便于与电脑进行数据交换。 因此，我们小组选择tms320vc5416芯片来进行开发。2 硬件电路概述2.1 dsp开发系统 通用dsp开发系统的三大要素为：1.界面友好、调试方便、具有图形显示功能以及能对程序问题快速定位的开发软件；2.与pc机接口方便、与目标板连接方便、硬件可靠以及通用性好的开发系统；3.具有流行dsp芯片、具有a/d、d/a等丰富外围接口、扩展方便意见实验软件丰富的目标板；此处仅对开发系统以及目标板进行具体描述。2.2 开发系统 在这次语音信号滤波处理设计中，我们小组采用icetek-5100pp系列开发系统，其外观图如图9所示: 图9 icetek-5100pp开发系统外观图其各部件的说明如下：a.仿真头：jtag接口，其各管脚定义如图10所示：图10 jtag管脚定义b.仿真头：mpsd接口，其各管脚定义如图11所示:图11 mpsd仿真接口各管脚定义c.开发系统电源接口d.并口电缆e.开发系统供电指示灯f.正确运行指示灯2.2.1 系统工作环境 icetek-5100pp开发系统对主的硬件要求如下：icetek-5100pp采用并口与计算机相连接，在使用开发系统之前需要正确配置并行端口和开发系统的相应设置。一般的，台式机并行端口有三种工作模式，分别是：spp、epp、ecp；而笔记本电脑则只支持其中的spp和ecp两种工作模式。icetek-5100pp可以再spp和epp模式下工作，所以要把计算机的并行端口配置为spp或epp模式。如果主机支持epp工作模式，则推荐使用epp模式。其中，修改并行端口工作模式需要在计算机的bios中操作。2.2.2 安装步骤 在使用icetek-5100pp系列开发系统之前，需要先安装开发软件，其安装步骤如下：1. 关闭计算机电源；2. +5v电源与开发系统的电源接口相接；3. 仿真头与目标板相接；4. 打开目标板电源，这是开发系统供电指示灯会亮；5. 打开计算机；注意：在安装过程中，要在避免接通计算机电源和目标板电源的情况下插播接口，否则会烧坏计算机并口或开发系统。2.3 系统板 本次设计中，我们将会选用icetekvc5416-a作系统板，其实物图如图12所示：图12 icetekvc5416-a实物图其中各接口具体如下：a. jtag仿真插座，与icetek-5100pp开发板的a接口相连；b. 四位设置拨码开关；c. hpi设置跳线j1、j2；d. u320：四位拨动开关输入；e. +5v电源插座；f. 电源指示；g. 四位用户可控指示灯：d3d6；h. p3扩展插座；i. pi扩展插座；j. 标准rs232接口；k. 语音输入输出接口：mic in、line in、heda phone、line out； l. p2扩展插座；m. p4扩展插座；n. 复位按钮；dsp系统板icetekvc5416-a的模数转换模块：精度为12bit，同相位最多采集路数为6，两路工作时，最高采集数率为450k/路，而6路工作时，则是150k/路，信号输入范围为05v，而信号耦合方式既可以是直流也可以使交流；dsp系统板icetekvc5416-a的数模转换：精度为12bit,路数为4路，其交换数率为100k（10us），输出范围为05v；dsp系统板icetekvc5416-a的周边设计：其主处理器是tms320vc5416，内部存储空间为128k16bit，具有看门狗功能设计。它的上位机接口为rs232串行数据接口，即接口j，同时具有dsp扩展总线。2.4 硬件电路连接1.将dsp开发系统icetek-5100pp的jtag仿真插座连接到实验系统的dsp模块icetekvc5416-a的jtag上，如图13所示：计算机开发系统icetek-5100pp系统板icetekvc5416-aspp/epp/ecpjtag 图13 硬件电路各部件连接说明2.系统上电后，将开发系统icetek-5100pp与计算机相连；3.双击桌面的set up ccs图标，配置并启动ccs软件；3软件编程实现分功能模块写，代码&代码功能说明，以及代码算法思想，和流程图为实现本设计，需在matlab协助中进行软件仿真，仿真成功后，再通过ccs进行硬件仿真、编译并通过串口将程序下载到dsp芯片之中执行。3.1 软件仿真实现3.1.1 matlab简介matlab是美国mathworks公司于1984年正式推出的一套高性能的数值计算和可视化软件，适用于工程应用各个领域的分析设计和复杂计算，是一种功能性很强、效率高、便于进行科学和工程计算的交互式软件包，它集数值分析、矩阵运算、信号处理和图形显示于一体，为用户提供方便和友好的界面环境。 matlab中的工具箱中包含了很多使用程序，如数值分析、矩阵运算、数字信号处理、建模和系统控制等。滤波器的设计就包含在该工具箱的signal中，它提供了很多种fir滤波器设计方法。 matlab因其语言简洁紧凑，使用方便灵活，库函数极其丰富，又具有结构化的控制语句(如for循环，while循环，break语句和if语句)，又有面向对象编程的特性;程序的可移植性很好，基本上不做修改就可以在各种型号的计算机和操作系统上运行，等等优点。因此在各个学科和领域得到了广泛的应用。3.1.2 数字滤波器在matlab的实现 其设计流程如图14所示：初始化产生需要虑波信号 画出原始信号频谱图 生成fir滤波器 进行滤波 画出滤波后信号频谱图 图14 matlab中设计流程图在matlab中，可以利用矩形窗、三角窗、汉宁窗、汉明窗、布莱克曼窗、凯塞窗等设计fir滤波器，在本次课程设计中主要应用汉宁窗设计出fir滤波器。利用matlab中的函数freqz画出各滤波器的频率响应，首先利用数字信号处理里面学过的知识，根据自己选定的参数，用汉宁窗函数法设计fir数字滤波器，得到数字滤波器的参数b，a。其中b为系统函数的分子系数，a为系统函数分母系数。再调用freqz(b,a,512,fs)即可得到该滤波器的频率响应。主程序如下：fpd=1700;fsd=1750;fsu=1850;fpu=1900; % fir滤波器的上下截止频率rp=1;as=36; % 带阻滤波器设计指标fcd=(fpd+fsd)/2;fcu=(fpu+fsu)/2;df=min(fsd-fpd),(fpu-fsu); % 计算上下边带中心频率，和频率间隔wcd=fcd/fs*2*pi;wcu=fcu/fs*2*pi;dw=df/fs*2*pi; % 将hz为单位的模拟频率换算为rad为单位的数字频率wsd=fsd/fs*2*pi;wsu=fsu/fs*2*pi;m=ceil(6.2*pi/dw)+1; % 计算汉宁窗设计该滤波器时需要的阶数n=0:m-1; % 定义时间范围w_ham=hanning(m); % 产生m阶的汉宁窗 hd_bs=ideal_lp(wcd,m)+ideal_lp(pi,m)-ideal_lp(wcu,m); % 调用自编函数计算理想带阻滤波器的脉冲响应h_bs=w_ham.*hd_bs; % 用窗口法计算实际滤波器脉冲响应db,mag,pha,grd,w=freqz_m(h_bs,1); % 调用自编函数计算滤波器的频率特性阻带最大衰减为-75db,fir滤波器的主瓣宽度很小,这样可以使过渡带很陡，旁瓣相对于主瓣也比较小。3.2 硬件仿真实现3.2.1 ccs简介dsp芯片的开发需要一套完整的软、硬件开发工具,可以分为代码生成工具和代码调试工具两类.而ccs(codecomposerstudio)就是ti推出的用于开发其dsp芯片的集成开发环境，它采用windows风格界面，集编辑、编译、链接软件仿真、硬件调试及实时跟踪等功能于一体，极大的方便了dsp程序的设计与开发。ccs集成开发环境中包含simulator(软件仿真器)和emulator(硬件仿真器)两部分。它们使用的是同一个集成开发环境。在本设计中，采用软件仿真器进行程序的调试。3.2.2 matlab与ccs的连接 matlab辅助dsp开发实现的关键是建立matlab与dsp间的连接。以往一般是由开发工具matlab把仿真结果先保存，再调入ccs中，在ccs中的仿真中间结果与matlab的仿真结果进行比较，以此发现dsp程序的不足，这需要反复操作，比较麻烦。 mathworks公司和ti公司共同开发的matlablinkforccs开发工具(ccslink)，实现了在matlab，ticcs开发环境和dsp硬件间的双向连接，开发者可以利用matlab强大的数据处理、分析、可视化功能来处理ccs和目标dsp中的数据，可以大大简化dsp软件开发的分析、调试和验证过程，缩短软件开发周期。matlab可通过3种方式与ccs、目标dsi，进行连接、数据交换。ccslink提供了3种连接对象:与ccs的连接对象可从matlab命令窗运行ccs中的应用程序，向目标dsp的存贮器、寄存器读出/写人数据，检查dsp状态，开始/停止目标dsp中运行的程序。与rtdx(实时数据交换)的连接对象使matlab与目标dsp直接通信，matlab可以实时地向目标dsp取出/发送数据，并不停止dsp中正在执行的程序。嵌入式对象在matlab环境中创建，该对象可代表嵌入在目标c程序中的变量，由其可以直接对嵌入在目标dsp存贮器/寄存器中的变量进行操作。3.2.3 数字滤波器在ccs中的实现其设计流程如图15所示：图15 dsp程序设计流程图3.1音频文件生成y,fs,nbits=wavread(dsp.wav);y=y(:,1);n=length(y);f=fs*(0:511)/1024;t=0:1/22050:(n-1)/22050;au=0.3;d=au*cos(pi*6000*t); %噪声为6khz的余弦信号x2=y+d; %加噪声后的语音信号%sound(x2,fs); %播放加噪声后的语音信号wavwrite(x2,input.wav);n2=length(x2); y1=fft(y,n);y2=fft(x2,1024); figure(1);subplot(2,1,1);plot(y);title(原始的信号波形);subplot(2,1,2);plot(x2);title(加噪后的信号波形); figure(2);plot(f,abs(y2(1:512);title(噪音信号频谱); figure(3);subplot(2,1,1)plot(abs(y1);title(原始信号频谱);subplot(2,1,2)plot(abs(x2);title(加噪后信号频谱);xout=x2/max(x2);% 归一化? xto_ccs=round(32767*xout) ;%clear playsnd;fid=fopen(dspindata.dat,w);% 打开文件fprintf(fid,1651 1 0 0 0n);% 输出文件头fprintf(fid,%dn,xto_ccs);% 输出fclose(fid); 用malab生成混入杂声后的音频文件3.2音频数据读入：stm #temp,ar1 ;temp存输入左声道样本stm #temp+1,ar2;temp+1存输入右声道样本ldm drr10,astl a,*ar1 ;ar1存左声道，ar2存右声道ldm drr20,astl a,*ar2分左右声道分别读入音频样本数据，放入寄存器，等待进行处理3.3判断是否进入fir处理的分支：stm #echo,ar4;是否fir处理的分支cmpm *ar4,#1bc fir,tcb nofir根据寄存器数据内容，判断是否进入fir处理的分支3.4fir处理分支内容（核心部分） fir:ssbx frctstm #temp,ar1;temp存输入左声道样本ld *ar1,a;stm #buffer+1,ar3;stl a,*ar3;将最新样本存入fx1,fx0保留空间,存fir后的输出暂存stm #buffer+25,ar2;将ar2指针指向最高位（最老数据）rptza,#24;24macd*ar2-,coef,astha,*ar2;把fir后的数据放在fx0rsbxfrctld*ar2,astlm a,dxr10;把fir后的数据输出语音端口nopstlm a,dxr20;这里两个声道一起输出同一个数据rsbxfrctnopb l1把读入的声音样本数据在芯片内进行fir处理4结果分析4.1 软件仿真运行结果及分析 为了更好的实现语音信号的滤波处理，我们在matlab上进行软件仿真。在matlab中模拟语音滤波器的功能，将语音信号中的高频噪声去除，还原歌曲音频。 图16 滤波器频率响应图4.2 硬件仿真运行结果及分析 未经滤波的混有噪声的音频频域波形和经过ccs滤波后还原的歌曲音频的时域和频域波形分别如图： 图17 滤波前的音频时域波形图图18 滤波后的音频时域波形图图19 滤波前的音频频域波形图图20 滤波后的音频频域域波形图 比较硬件仿真后得到的还原音频和matlab软件仿真得到的还原音频，可得经ccs硬件仿真后的还原的音频不及matlab中的还原的音频与原来的相近，这是因为数字滤波器的性能主要取决于乘法器的各系数，而这些洗漱在matlab上仿真的时候是没有经过处理的真实数据，因此在matlab上仿真是理想的仿真，因而效果比较好。而在ccs上的各个系数是经过量化处理的，因此系数和原始系数比较会有变化，令在ccs中仿真滤波的结果不如matlab中的好，但是并不影响它的正确性。5问题，调试 1.用ccs3.3进行硬件仿真时，无法打开仿真软件，判断为软硬件不兼容，改用ccs2.3实现，问题解决。2.绘制出滤波后的波形，发现fir滤波器没有滤掉单频噪声。通过自己的仔细检查，是单频噪声的频率改动后，fir滤波器的频率没有改动。所以单频噪声的频率也应该自己先定义，fir滤波器的截止频率应该以单频噪声的频率为中心，这样重新运行后，结果正确。3.采集的语音信号中高频噪音不很明显，使得最后不能较好的观察到滤波器的性能。解决的方法是:利用matlab产生一个6khz的高频信号与原来的语音信号相加，这时可以明显的听到语音信号中同时伴随着刺耳的高频噪音。4.在选择滤波器的参数时,对具体的音频信息没有整体的了解，使得最后fir滤波器不能达到很好的滤波效果。解决的方法是，首先对音频信息进行采集，加上高频噪声后，再整体的观察语音信号的波形、频谱，对高频段进行分析，设计满足要求的滤波器。5.在ccs调试过程中，导入indata.dat文件没有完全体现在波形图和频谱图上。解决的方法是,经过仔细的推敲，原因是导入时没有正确的设定数据的起始地址，再重新设定地址为0x00e5和长度为256后，正确的显示出了对应的波形图和频谱图。6.在编写ccs程序时，没有正确的设定缓冲区长度，引起了对滤波器性能的误判断。解决方法是，重新设定缓冲区的长度为40。6 设计改进 fir数字滤波器的设计有多种方法：有窗函数法，频率采样法，切比雪夫逼近法等。其中，窗函数不仅有截断作用，而且能够起到平滑的功能。加窗虽然会使原信号时域波形发生较大变化，但却更有效地保留了原信号的频率信息，在很多领域得到了应用。对不同类型的信号，因为由不同的窗函数实现的fir滤波器其性能不同，因此在截断处理中不同的滤波器所用的窗函数也不相同。在各种窗函数中，三角窗阻带最小衰减最小，其次是汉宁窗，再次是哈明窗，布莱克曼窗阻带最小衰减最大。主瓣宽度方面，三角形窗，汉宁窗以及哈明窗一样，布莱克曼窗比以三种的主瓣要宽，表现在用布莱克曼窗的滤波器的过渡带比用其他三种窗函数的要宽。主瓣高度与第一旁瓣高度的差值方面，三角窗高度差最小，其次是汉宁窗，再次是哈明窗，blackman窗高度差最大。旁瓣衰减的速率，三角形窗和哈明窗比较缓慢，汉宁窗相对前丽种要快，布莱克曼窗很快。旁瓣幅度汉宁窗和哈明窗较大，三角形窗和布莱克曼窗较小。所设计的fir数字滤波器的阶数n不同，所得到的滤波器的过渡带会不一样，n越小，过渡带越宽。而对于同一种窗，n的值越大主瓣宽度越小，相应地，滤波器的过渡带就越窄，主瓣高度与第一旁瓣高度的差值方面基本不变 而旁瓣的个数越来越多，但旁瓣的宽度随着n的加大而减小。因为hamming窗的fir滤波器使99963的能量集中在主瓣内，旁瓣峰值小于主瓣峰值的1，是各种窗函数中性能较好的一种，因此它被广泛应用。而在四种窗函数中，因为布莱克曼窗的旁瓣峰值最低为一57db可以充分抑制干扰频谱的旁瓣泄露，所以它的干扰抑制性能最为稳定。 针对采用不同的窗函数所造成的频谱泄露问题，则要采用采用hamming窗，因为它可将9999的能量集中在窗谱的主瓣内，有效地消除频谱泄漏，并能在一定程度上减少栅栏效应。7 课程设计心得 在本次课程设计的过程中，虽然遇到很多问题，让组内成员一度灰心丧气，但经过组内各成员同心协力的查找问题，翻阅资料以及讨论解决方案后，成功的完成了课程设计的内容。此过程中，我们不但能够将课堂上学到的知识加以运用，通过实际应用掌握dsp的操作以及原理。同时，也培养了我们克服困难战胜困难的能力，让我们懂得了团队合作的重要性，学会根据组内每个成员自己的特点以及长处不同，针对每个成员对该设计的完成在不同方面的独特作用，充分地了解各成员特性，更好的发挥团队的作用，最大化的提升团队的人员配置，对更好的完成该设计有着不可忽视的作用。 另外，此次设计的过程中，除了增加了对滤波器的设计、实现的认识外，还增加了对matlab和ccs的认识和了解，matlab具有强大的数值计算和图形功能，虽然其执行效率较低，但其编程效率与可读性、可移植性要远高于其他高级语言。它集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个界面友好的用户环境，是一个非常优秀的算法研究与辅助设计的工具。运用它来进行辅助设计，在系统的准确性和高效率方面起到了至关重要的作用。而ccs加速和增强了实时，嵌入信号处理的开发过程，它提供配置、构造、调试、跟踪和分析程序的工具，在基本代码产生工具的基础上增加了调试和实时分析的功能。开发设计人员可在不中断程序运行的情况下检查算法的对错，实现对硬件的实时跟踪调试，从而大大缩减了程序的开发时间。另外，也学会了利用ccslinl命令在matlab 和ccs之间进行通信，提高操作的方便度和速度。 附录主函数：.titlefir.asm.mmregs.defstart.def rxdataleveln.set 25buffer.usectbuffer,leveln+2;25级数组空间，有1个是保留空间in.usectin,8192;输入的数据保存区out.usectout,8192;输出的数据保存区stack.usectstack,400h;堆栈的长度 .bss kong,2 .bss echo,1.bssncount,1.bsstemp,5.bss temp1,1.bss temp2,5;*;* 系统变量定义 *;*cpu_speed_160m.set9007hleft_ptrq.set4000hright_ptr.set8000hflag.set0060hdatabuffsize .set 8192delaytime.set 0 ;6134;*;* mcbsp变量预定义 *;*spsa0.set38hspsd0 .set39hspsa1.set48hspsd1.set49hspcr1.set00hspcr2.set01hrcr1.set02hrcr2.set03hxcr1.set04hxcr2.set05hsrgr1.set06hsrgr2.set07hmcr1.set08hmcr2.set09hrcera.set0ahrcerb.set0bhxcera.set0chxcerb.set0dhpcr.set0ehrcerc.set10hrcerd.set11hxcerc.set12hxcerd.set13hrcere.set14hrcerf.set15hxcere.set16hxcerf.set17hrcerg.set18hrcerh.set19hxcerg.set1ahxcerh.set1bhdrr20.set 20hdrr10 .set 21hdxr20 .set22hdxr10 .set 23hdrr21.set 40hdrr11 .set 41hdxr21 .set42hdxr11 .set 43h;*;* 语音aic23的参数预定义 *;*aic23_reset_reg.set0fhaic23_power_down_ctl.set06haic23_analog_audio_ctl.set04hanapctl_dac.set10haic23_digital_audio_ctl.set05haic23_lt_line_ctl.set00haic23_rt_line_ctl.set01haic23_digital_if_format.set07hdigif_fmt_lrp.set10hdigif_fmt_ms.set40hdigif_fmt_iwl_16.set00hdigif_fmt_for_dsp.set03haic23_sample_rate_ctl.set08hsrc_sr_08.set0chaic23_dig_if_activate.set09hdigifact_act.set01hsrc_bosr.set02hsrc_mo.set01hhechengdigif.set digif_fmt_lrp|digif_fmt_ms|digif_fmt_iwl_16|digif_fmt_for_dsp;*;* fir滤波器的系数 *;*.datacoef:.word20*32768/10000,16*32768/10000;b24,b23.word0*32768/10000,-44*32768/10000;b22,b21.word-117*32768/10000, -181*32768/10000;b20,b19.word-168*32768/10000,0*32768/10000;b18,b17.word359*32768/10000,870*32768/10000;b16,b15.word1415*32768/10000,1835*32768/10000;b14,b13.word1992*32768/10000,1835*32768/10000 ;b12,b11 .word1415*32768/10000,870*32768/10000;b10,b9.word359*32768/10000,0*32768/10000;b8,b7.word-168*32768/10000, -181*32768/10000;b6,b5.word-117*32768/10000,-44*32768/10000;b4,b3.word0*32768/10000,16*32768/10000;b2,b1.word20*32768/10000;b0;*;* 主程序 *;*.textstart:stm#stack+400h,sp;设置堆栈stm#temp,ar1;向3005h的io端口写1st #1,*ar1portw*ar1,3005hstm #0,clkmd;初始化频率rpt #50nopstm#cpu_speed_160m,clkmdrpt #50nop;*stm #spcr1,spsa1;mcbsp1的初始化stm #0,spsd1stm #spcr2,spsa1stm #0,spsd1stm #spcr1,spsa1stm#1800h,spsd1stm #rcr1,spsa1stm#040h,spsd1stm#xcr1,spsa1stm#040h,spsd1stm#srgr2,spsa1stm#2000h,spsd1stm #srgr1,spsa1stm#0ffh,spsd1stm #xcr2,spsa1stm#0,spsd1stm #rcr2,spsa1stm#0,spsd1stm #pcr,spsa1stm#0a08h,spsd1stm#spcr2,spsa1orm#40h,spsd1rpt 100nopstm #spcr1,spsa1orm#01h,spsd1stm #spcr2,spsa1orm #01h,spsd1rpt 100nop;*stm #temp,ar1 ;aic23初始化aic23_write(aic23_reset_reg, 0);st #aic23_reset_reg,*ar1stm #temp+1,ar2st #0,*ar2ld *ar1,9,aadd *ar2,astl a,*ar1stm #spcr2,spsa1testma1:;mvkd spsd1,*ar2ldm spsd1,astl a,*ar2bitf *ar2,#02hbc sendready1,tcb testma1sendready1:;mvkm,*ar1,dxr11;ld *ar1,astlm a,dxr11;stm #temp,ar1 ;aic23_write(aic23_power_down_ctl, 0); st #aic23_power_down_ctl,