语音信号的FIR滤波器处理课程设计_.docx

语音信号的FIR滤波器处理课程设计_此页面能否是列表页或首页？未找到适宜正文内容。语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计摘要随着信息与数字技术的发展，数字信号处理已经成为当今极其重要而学科与技术领域之一。它在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等诸多领域得到了广泛的应用。在数字信号处理的基本方法中，通常会涉及到变换、滤波、频谱分析、调制解调和编码解码等处理。其中滤波是应用非常广泛的一个环节，数字滤波器的理论和相关设计也一直都是人们研究的重点之一。FIR滤波器的是非递归的，稳定性好，精度高；更重要的是，FIR滤波器在知足幅频响应要求的同时，能够获得严格的线性相位特征。因而，它在高保真的信号处理，如数字音频、图像处理、数据传输和生物医学等领域得到广泛应用。CSS集成开发环境使用CCS内置的软件仿真simulator对程序进行编译，调试和运行，主要用于检测目的程序运行的正确性和连贯性,并能通过仿真器与目的板连接，在目的板上实时观察效果。在本次设计中，我们选择的课题是基于DSP的语言信号的FIR滤波处理。首先利用MATLAB进行了仿真，得到滤波前后的时域波形和频谱。然后通过调用MATLAB的分析工具FDATOOL，根据仿真结果导出了滤波器的相关参数，将原始信号数据和滤波器参数输入CCS进行DSP编程。最后在DSP中实现了FIR低通滤波，并通过CCS的频谱分析功能查看了最终DSP的滤波效果。关键词:语音信号处理FIR滤波MATLABCCS语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计目录1.前言(4)2.设计需求(4)3.滤波器设计原理介绍(4)3.1数字滤波器设计原理(4)3.2FIR滤波器的基本构造(5)3.3滤波器的特点(5)3.4窗函数的介绍(6)4.FIR滤波器的设计(7)4.1FIR滤波器设计(7)4.2窗函数设计的基本方法(7)4.3滤波器的1z算法实现(7)5.FIR滤波器的MATLAB设计(8)6.FIR滤波器的DSP实现(10)6.1CCS设计与仿真(10)6.2CCS仿真结果(13)6.调试问题及解决办法(16)7.设计感想(16)参考文献(17)8.附录(17)语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计1前言通过语音传递信息是人类最重要、最有效、最常用、最方便的交换信息的方法，因而，怎样去除混杂在有用信号中的噪声并获得低损甚至无损的有用语音信号成为语音信号处理即用数字信号处理技术和语言学知识对语音信号进行处理并提取有用信息至关重要的问题。但是由于语音识别技术对信号噪声非常敏感，系统本身和器件对系统的处理结果都有着不可忽视的干扰，因而很难实现效果较好的语音信号滤波处理。而DSP由于运算速度快，片上资源丰富和能够实现复杂的线性和非线性算法等特性，在语音信号处理技术方面有明显的优势。在数字信号处理中，滤波占有极其重要的地位。数字滤波是语音信号处理、图像处理、形式识别、频谱分析等应用的基本处理算法。用DSP芯片实现数字滤波除了具有稳定性好、准确度高、不受环境影响等优点外，还具有灵敏性好等特点。本文介绍了数字滤波器的设计基础及用窗函数法设计FIR滤波器的方法，运用MATHLAB语言实现了低通滤波器的设计并用CCS2.0进行观察效果。2设计需求1滤波器的设计的原理介绍2FIR滤波器的设计3FIR滤波器的DSP实现4滤波器为低通滤波器，并用用语音信号通过滤波器后查看滤波效果能否到达设计要求，进行FIR滤波器的检验根据教师的的意见，我们本人拟定了滤波器的指标：Fs=8000hz，Wp=1500Hz，Ws=2000Hz，通带波纹为0.01，阻带波纹为0.1，N=37。3滤波器设计原理介绍3.1数字滤波器的设计原理数字滤波器的设计问题就是寻找一组系数ai和bi，使得其性能在某种意义上逼近所要求的特性。假如在s平面上去逼近，就得到模拟滤波器，假如在z语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计平面上去逼近，则得到数字滤波器。数字滤波是将输入的信号序列，按规定的算法进行处理，进而得到所期望的输出序列。一个线性位移不变系统的输出序列y(n)和输入序列x(n)之间的关系，应知足常系数线性差分方程：2.1.1x(n)为输入序列，y(n)为输出序列，ai、bi为滤波器系数，N为滤波器的阶数。3.2FIR滤波器的基本构造在式2.1.1中，若所有的ai均为0，则得FIR滤波器的差分方程：2.2.1对式2.2.1进行z变换，可得FIR滤波器的传递函数：2.2.2由此可得到FIR滤波器的构造如图1所示。FIR滤波器的单位冲击响应h(n)是一个有限长序列。若h(n)为实数，且知足偶对称或奇对称的条件，即h(n)=h(N-1-n)或h(n)=-h(N-1-n)，则FIR滤波器具有线性相位特性。3.3滤波器的特点)()()(110-=-=ninyainxbnyMiiNii)()(10-=-=Niiinxbny)()()(10-=-=NiiizbzXzYzH语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计优点：1很容易获得严格的线性相位，避免被处理的信号产生相位失真，这一特点在宽频带信号处理、阵列信号处理、数据传输等系统中非常重要；2能够在幅度特性随意设计的同时，保证准确、严格的线性相位，还可得到多带幅频特性；3极点全部在原点永远稳定，无稳定性问题；4任何一个非因果的有限长序列，总能够通过一定的延时，转变为因果序列，所以因果性总是知足；5无反应运算，运算误差小。缺点：1由于无极点，要获得好的过渡带特性，需以较高的阶数为代价；2无法利用模拟滤波器的设计结果，一般无解析设计公式，要借助计算机辅助设计程序完成。3.4窗函数简介数字信号处理的主要数学工具是傅里叶变换。而傅里叶变换是研究整个时间域和频率域的关系。不过，当运用计算机实现工程测试信号处理时，不可能对无限长的信号进行测量和运算，而是取其有限的时间片段进行分析。做法是从信号中截取一个时间片段，然后用观察的信号时间片段进行周期延拓处理，得到虚拟的无限长的信号，然后就能够对信号进行傅里叶变换、相关分析等数学处理。无线长的信号被截断以后，其频谱发生了畸变，原来集中在f(0)处的能量被分散到两个较宽的频带中去了这种现象称之为频谱能量泄漏。为了减少频谱能量泄漏，可采用不同的截取函数对信号进行截短，截断函数称为窗函数，简称为窗。信号截短以后产生的能量泄漏现象是必然的，由于窗函数w(t)是一个频带无限的函数，所以即便原信号x(t)是有限带宽信号，而在截短以后也必然成为无限带宽的函数，即信号在频域的能量与分布被扩展了。又从采样定理可知，无论采样频率多高，只要信号一经截短，就不可避免地引起混叠，因而信号截短必然导致一些误差。语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计泄漏与窗函数频谱的两侧旁瓣有关，假如两侧瓣的高度趋于零，而使能量相对集中在主瓣，就能够较为接近于真实的频谱，为此，在时间域中可采用不同的窗函数来截短信号。4FIR滤波器的设计4.1FIR滤波器设计FIR滤波器的设计问题在于寻求一系统函数，使其应逼近滤波器要求的理想频率响应。4.2窗函数设计的基本方法从时域出发，设计h(n)逼近理想hd(n)，设理想滤波器的单位脉冲响应为hd(n)，则有所求得的一般是无限长的，且是非因果的。要想得到一个因果的有限长的滤波器h(n)，最直接的方法是截断，或者讲用一个窗口函数对进行加窗处理，即。因而，应选择适宜的窗函数。4.3滤波器的1-z算法实现FIR滤波器的输出表达式为式中，为滤波器系数；x(n)表示滤波器在n时刻的输入；y(n)为n时刻的输出。它的基本算法是一种乘法-累加运算，即不断地输入样本x(n)，经过延时后，再进行乘法-累加，最后输出滤波结果y(n)。1线性缓冲区法-=-=10)()(NnnznhzH)(jdeH?-=deeHnhenheHjnjddjnndjd)(21)()()()(nhd)(nhd)(nw)(nhd)()()(nwnhnhd=jezjzHeH=|)()(语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计线性缓冲区法又称延迟线法。其特点：·对于N级的FIR滤波器，在数据存储器中开拓一个N单元的缓冲区滑窗，用来存放最新的N个输入样本；·从最老样本开场取数，每取一个样本后，将此样本向下移位；·读完最后一个样本后，输入最新样本存入缓冲区的顶部。2循环缓冲区法本次设计的FIR滤波器所采用的就是循环缓冲区法。循环缓冲区法的特点如下：·对于N级FIR滤波器，在数据存储器中开拓一个N单元的缓冲区滑窗，用来存放最新的N个输入样本；·从最新样本开场取数；·读完最后一个样本最老样本后，输入最新样本来代替最老样本，而其他数据位置不变；·用片内BK循环缓冲区长度寄存器对缓冲区进行间接寻址，使循环缓冲区地址首尾相邻。5FIR滤波器的matlab实现MATLAB是一种功能强、效率高、便于进行科学和工程计算的交互式软件包，它集数值分析、矩阵运算、信号处理和图形显示于一体，为用户提供了方便、友好的界面环境。MATLAB中的工具箱Toolbox包含了很多实用程序。它提供了多种FIR滤波器设计方法。fir1函数用来设计标准频率响应的基于窗函数的FIR滤波器，可实现加窗线性相位FIR数字滤波器的设计。详细语法如下：b=fir1(n，Wn)b=fir1(n，Wn，ftype)b=fir1(n，Wn，Window)b=fir1(n，Wn，ftype，Window)语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计其中n为滤波器的阶数；Wn为滤波器的截止频率；ftype为用来决定滤波器的类型，当ftype=high时，可设计高通滤波器；当ftype=stop时，可设计带阻滤波器。Window为用来指定滤波器采用的窗函数类型，Window参数可采用的窗口函数有：Boxcar，Hanning，Bartlett，Blackman，Kasier和chebwin等，其默认时为Hamming窗，进而得到滤波器的系数。fir2函数用来设计有任意频率响应的各种加窗FIR滤波器。详细语法如下：b=fir2n，f，mb=fir2n，f，m，Windowb=fir2n，f，m，nptb=fir2n，f，m，npt，Windowb=fir2n，f，m，npt，lapb=fir2n，f，m，nptt，lap，Window其中n为滤波器的阶数；f为频率点矢量；m为幅度点矢量；Window用来指定所使用的窗函数类型，默认值为汉明Hamming窗；npt用来指定fir2函数对频率响应进行内插的点数；lap用来指定fir2函数在重复频率点附近插入的区域大小，进而得到滤波器的系数。详细的代码见附录，输入之后能够导出滤波器的参数。CCS的输入语音信号的辅助产生MATLAB1、理想情况：FIR滤波器的输入数据存储在indata.dat文件中，共300个点，由频率为1000Hz和2500Hz的两个余弦波叠加而成2、为了验证设计的可行性，本身录制一个语音片段，主要信号为低频成分，其中夹杂着高频成分的噪声。在MATLAB中产生一个indata.dat文件供CCS输入，在matlab中输入的代码为：x=wavread('D:matlabworkXINZAO.wav');%语音文件的途径y=fft(x,1024);plot(abs(y(1:1000);%采样100个点x1=round(1024*x);fid=fopen('indata.dat','w');%生成一个dat文件语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计fprintf(fid,'1651100E21000n');fprintf(fid,'%dn',x1);%输出6FIR滤波器的DSP实现CCS是一种针对TMS320系列DSP的集成开发环境，在windows操作系统下，采用图形接口界面，提供环境配置、源文件编辑、程序调试、跟踪和分析等工具，提供了配置、建立、调试、跟踪和分析程序的工具，进而完成编辑、编译、链接、调试和数据分析等工作。6.1CCS设计与仿真语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计6.1.1程序分析FIR滤波器的DSP实现主要有线性缓冲区法和循环缓冲区法两种。线性缓冲区法又称延迟线法。其特点：1对于N级的FIR滤波器，在数据存储器中开拓一个N单元的缓冲区滑窗，用来存放最新的N个输入样本；2从最老样本开场取数，每取一个样本后，将此样本向下移位；3读完最后一个样本后，输入最新样本存入缓冲区的顶部。循环缓冲区法的特点如下：1对于N级FIR滤波器，在数据存储器中开拓一个N单元的缓冲区滑窗，用来存放最新的N个输入样本；语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计2从最新样本开场取数；3读完最后一个样本最老样本后，输入最新样本来代替最老样本，而其他数据位置不变；4用片内BK循环缓冲区长度寄存器对缓冲区进行间接寻址，使循环缓冲区地址首尾相邻。本次设计的FIR滤波器所采用的就是循环缓冲区法。对DSP进行初始化，并设置好存储单元。为数字滤波作准备，将预先设计好的有N个抽头的FIR数字滤波器的冲激响应序列h(n)中的N个数值放入存储单元bn中，从indata数据存储段中读入抽样值，将抽样值逐一放入存储单元xn中。主要循环程序如下：将累加器清零，将存储单元bn与xn值对应相乘并累加到ACC中，循环次将ACC中的值保存输出。再从indata数据存储段中读入一个新的数据，代替原来的最老样本MVDD*AR4+,*AR2+0%,重复循环直到数据全部处理完。最后将处理完的数据输入到outdata数据存储段中。6.1.2数据输入MATLAB辅助DSP实现FIR滤波器滤波,其总体经过是为DSP中编写处理程序提供滤波器系数。在MATLAB中利用滤波器设计、分析工具(FDATOOL),根据指定的滤波器性能快速设计一个FIR,然后把滤波器系数以头文件形式导入CCS中,头文件中包含MATLAB辅助DSP实现FIR数字滤波器的阶数和系数数组,在MATLAB中调试、运行DSP程序并显示、分析处理后的数据。打开CCSIDE中file的data加载初始化data数据，在Address中输入数据存放的起始地址indata，并设置数据栈长度。语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计6.1.3运行程序仿真运行程序，查看输入输出波形，修改相应参数进行调试。在View的Graph中单击Time/frequency出现graphpropertydialog框。将显示类型，图形名称，起始地址，抽样点数，数据类型等分别进行设置，输出各种波形。输出滤波前的信号时域波形。首先，设置如下列图所示，输出滤波前的信号时域波形。图2滤波前语音信号的时域GraphPropertyDialog设置参数图将DisplayType改为FFTMagnitude就能够看到滤波前信号的频谱，将StartAddress改为outdata就能够看到滤波后信号的时域波形和频谱。6.2CCS仿真结果此页面能否是列表页或首页？未找到适宜正文内容。语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计图5滤波后信号时域波形图图6滤波后信号频谱图图3波前的信号，由高频信号和原信号叠加而成的波形较为杂乱，加噪后的信号幅值较大。图5波后的信号，从图上能够看到信号幅度较低，波形较为明晰。比照图4图6能够看到加躁后信号幅值较大，信号的能量集中在高频段，滤波后，高频成分大部分被滤去，信号的能量主要集中在低频段。语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计7调试问题及解决方法1实验时，未做初始化设置直接进入，没有装入gel文件及相关参数，导致输入代码编译链接后不能载入程序。解决方法：初始化设置软件仿真cpu和硬件仿真cpu，并设置其中所用的gel文件，保存并退出启动dsp选择忽略进入编辑编译初始化界面，由软件仿真cpu进入编辑编译工具。2在解决问题1后，载入程序，载入数据，并设置载入数据的初始地址和长度，运行程序，发现程序一直卡死在某个地方，无法显示正确的时域图和频域图。解决方法：换了一台电脑运行，发现问题已被解决，可能原因是本人电脑安装时装少了一些文件。8设计感想刚开场面对大作业的时候，我们都还是觉得很陌生，不知从哪里入手，尽管做过实验。对于本人选的题目，详细需要做什么，能不能做出来，我们并不清楚。但当时间一天天过去后，我们开场认真起来了。从温习课本，再到温习实验，互相讨论编程思路，上网搜索资料，调试实验等，我们都在努力的自学着。首先，是学习上的稳固。一方面，通过温习课本和实验，对DSP的基础知识又有了很大的稳固。其次，通过对用Mathlab实现FIR滤波器的设计，熟悉了matlab软件的一些相关的窗口函数以及相关功能的调用，如如何实现窗函数的调用，如何实现滤波，等等。第二，对FIR的滤波器性能和作用有了更深一层的了解。FIR滤波器的应用特别广泛，当今很多信号处理系统和图像处理系统等都要求信号具有线性相位特性。在这方面，FIR滤波器有特别独特的优点，运用mathlab语言，我们能够很容易的设计出具有严格线性相位的FIR滤波系统，以及比拟容易的实现。此外，我们还知道了FIR滤波器是永远稳定的，由于FIR滤波器的冲激响应是有限长序列，其系统函数为一个多项式，它所含的的极点多为原点。最后，也感谢徐梅宣教师的指导，跑了她办公室两次咨询情况，都很耐心解语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计答，真给了一些写报告或者软件应用的建议。参考文献1邹彦，DSP原理及应用修订版，北京：电子工业出版社，20212王宏.MATLAB6.5及其在信号处理中的应用.M.北京:清华大学出版社,20043高西全,丁玉美,阔永红.数字信号处理-原理、实现及应用.M.北京:电子工业出版社,20064胡广书编著，数字信号处理理论、算法与实现.北京：清华大学出版社，2005；附录1.主程序fir.asm：.title"FIR.ASM".mmregs.defstart.bssyn,1xn.usect"xn",37hn.usect"hn",37indata.usect"indata",300outdata.usect"outdata",256.datafir_coff:.word-10*1024/10000,-36*1024/10000.word0*1024/10000,65*1024/10000.word32*1024/10000,-88*1024/10000.word-93*1024/10000,90*1024/10000.word184*1024/10000,-47*1024/10000.word-297*1024/10000,-71*1024/10000.word417*1024/10000,316*1024/10000.word-524*1024/10000,-848*1024/10000.word598*1024/10000,3108*1024/10000.word4375*1024/10000.word3108*1024/10000,598*1024/10000.word-848*1024/10000,-524*1024/10000语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计.word316*1024/10000,417*1024/10000.word-71*1024/10000,-297*1024/10000.word-47*1024/10000,184*1024/10000.word90*1024/10000,-93*1024/10000.word-88*1024/10000,32*1024/10000.word65*1024/10000,0*1024/10000.word-36*1024/10000,-10*1024/10000.textstart:ssbxfrct;使用小数运算stm#hn,ar1;系数首地址rpt#36;将系数移入循mvpdfir_coff,*ar1+;环缓冲区stm#-1,ar0stm#outdata,ar5stm#indata,ar4stm#xn+36,ar2rpt#36;将输入数据移入mvdd*ar4+,*ar2+0%;循环缓冲区stm#xn+36,ar2;指向第一个输入stm#hn+36,ar3stm#37,bk;缓冲区大小37stm#255,brc;块重复256次stm#-1,ar0rptbloop-1;块程序重复大小rptza,#36;计算一个输出mac*ar2+0%,*ar3+0%,astha,*(yn);保存输出mvkd*(yn),*ar5+mvdd*ar4+,*ar2+0%;读进一个输入loop:nopbloop.end2链接命令文件fird：MEMORYPAGE0:PARAM:org=00E0hlen=1000hPAGE1:SPRAMorg=0060hlen=0020hDARAM:org=0080hlen=1380h语音信号的FIR滤波器处理课程设计语音信号的FIR滤波器处理课程设计SECTIONS.text:>PARAMPAGE0.data:>PARAMPAGE0.bss:>SPRAMPAGE1xn:align(64)>DARAMPAGE1hn:align(64)>DARAMPAGE1indata:>DARAMPAGE1outdata:>DARAMPAGE1