语音信号去噪——数字滤波器的设计本科毕业论文设计

皖西学院本科毕业论文设计第1页语音信号去噪数字滤波器的设计摘要在现代各种通信系统中，由于自然界中的各种各样的复杂噪声不免会掺杂在其中，数字信号处理这门经典学科恰好能够解决这个问题，其中最通用的方法就是利用滤波器来滤除这些杂波噪声，FIR数字滤波器就是滤波器设计的基本部分。本论文研究的主要内容就是基于MATLAB软件仿真设计一个数字滤波器，将掺杂在语音信号中的高频噪音消除，在此将分析消除高频噪音前后语音信号的时域及频域特性，对比分析即可验证滤波前后特性差别。在本课题中，将利用简单的窗函数法来设计FIR数字滤波器，通过MATLAB仿真说明所设计滤波器的正确性。通过这次毕业设计，将会进一步理解语音信号原理分析及滤波处理，为更好的设计滤波器打好基础。关键词MATLAB；窗函数法；FIR数字滤波器REMOVENOISEINTHESPEECHSIGNALTHEDESIGNOFDIGITALFILTERABSTRACTINMODERNCOMMUNICATIONSYSTEMS,AVARIETYOFCOMPLEXNOISEMAYMIXINTHENATUREOFSOUNDSTHECLASSICDISCIPLINESOFTHEDIGITALSIGNALPROCESSINGCANSOLVETHISPROBLEM,ONEOFTHEMOSTCOMMONMETHODISTOUSEAFILTERTOFILTERTHOSECLUTTERNOISEFIRDIGITALFILTERISTHEBASICPARTOFFILTERDESIGNTHEMAINRESEARCHCONTENTOFTHISPAPERISBASEDONMATLABSOFTWARESIMULATIONTODESIGNADIGITALFILTER,INWHICHTOCANCELTHEHIGHFREQUENCYNOISEOFTHESPEECHSIGNAL,THENITWILLELIMINATETHEHIGHFREQUENCYNOISEANDTHESPEECHSIGNALSFROMTIMEDOMAINANDFREQUENCYDOMAINCHARACTERISTICSINTHISANALYSISBEFOREANDAFTER,ANDANALYSISTHEDIFFERENCESTESTTHEFILTERINGCHARACTERISTICSINTHISISSUE,USINGASIMPLEWINDOWFUNCTIONMETHODTODESIGNAFIRDIGITALFILTER,MATLABSIMULATIONSHOWSTHECORRECTOFTHEDESIGNEDFILTERTHROUGHTHEDESIGNOFTHISGRADUATIONDESIGN,WEWILLUNDERSTANDTHEPRINCIPLEOFSPEECHSIGNALANALYSISANDFILTERING,ANDLAYTHEFOUNDATIONFORTHEFILTERDESIGNKEYWORDSMATLABWINDOWFUNCTIONMETHODFIRDIGITALFILTER皖西学院本科毕业论文设计第2页目录1引言311课题研究现状312课题研究目的313课题研究内容314MATLAB软件设计平台简介42原始语音信号采集与处理521课题设计步骤及流程图522语音信号处理5221语音信号的采集5222语音信号的时域频谱分析6223语音信号加噪与频谱分析83FIR数字滤波器的设计1031数字滤波器基本概念1032常用窗函数介绍1033FIR数字滤波器概述1034FIR滤波器的窗函数设计1135滤波器的编程实现1336用滤波器对加噪语音信号进行滤波1437回放语音信号164结论18致谢19参考文献20皖西学院本科毕业论文设计第3页1引言11课题研究现状20世纪60年代中期数字信号处理领域形成的诸多富有实践性的的理论和算法，如快速傅立叶变换（FFT）以及各种数字滤波器等是语音信号数字处理的各项理论和技术基础。在70年中后期之后，线性预测技术（LPC）已经用于语音信号的信息压缩和特征提取，并已成为语音信号处理中非常重要的一个工具。80年代语音信号处理技术的重大发展是隐马尔可夫模型（HMM）描述语音信号过程的产生过程。进入上世纪90年代以来，语音信号采集与分析技术在实际应用方面取得了许多突破性的研究进展。近年来，随着科学技术的不断进步，人工神经网络ANN的研究取得了迅速发展，语音信号处理的各项科学研究课题是促进其发展的催生力，同时，它的许多成果也体现在有关语音信号处理的各项技术之中。12课题研究目的语音是人类获取知识和各种各样信息的重要手段和最初来源，人类离不开自然界中各种不同的语音，但在获取语音的过程中，将不可避免的会受到外界环境的干扰和影响，如各种机器的轰鸣声或者自然界太多的电磁噪声干扰等这些不可避免的有害噪声信号都会附加掺杂在语音信号中，这样获取的语音信号将不再是单纯的语音，掺杂的噪声不但降低了语音质量和语音的可懂度，严重时将导致不可预知的不良效果。语音信号处理的好坏将影响并导致语音信号的好坏，只有通过将这些语音信号进行一系列的的数字处理，才能将那些非必要的噪声杂波妥善滤除，得到纯净的单纯的语音信号。现在社会衍生了很多现代的语音通信方式，手机通话、QQ或MSN等这些语音聊天软件以及语音小说等，涉及语音的方方面面已经存在于我们生活中的大部分，因此语音信号去噪处理是具有现实意义的研究课程。13课题研究内容1用MATLAB程序对原始语音信号进行采样，并绘制出采样后语言信号时域波形和频谱图，并针对此图分析语音信号特点。2将干扰噪声加入到已经获取的语音信号中，然后进行频谱分析，并对比未加入噪声的信号，分析差异。3设计FIR数字滤波器，并对被噪声污染的语音信号进行滤波，画出滤波前后信号的时域皖西学院本科毕业论文设计第4页和频谱，对滤波前后有噪音与无噪音的语音信号进行比较，分析信号的变化，从而验证所设计滤波器的滤波效果。4回放语音信号。14MATLAB软件设计平台简介MATLAB是MATRIX和LABORATORY的缩写，是由MATHWORKS公司开发的一套用于科学工程计算的可视化高性能软件，是一种交互式的以矩阵为基本数据结构的系统，应用广泛。MATLAB的主要特点如下1程序的可移植性良好应用于其他程序。2程序限制宽泛，程序设计自由。有大量已经系统定义的函数可直接应用，并且能够用户自定义函数。3语言简洁，使用灵活方便，库函数相当丰富。4源程序向大众开放。用户可灵活的对源文件进行修改以及加入自己的设计语音构成新的工具箱。5最后MATLAB的一个重要特点是功能强大的工具箱。MATLAB包含两个重要的部分核心部分和各种可选的工具箱。皖西学院本科毕业论文设计第5页2原始语音信号采集与处理21课题设计步骤及流程图本课题设计主要是针对一段原始语音信号，加入设计噪声后，用窗函数法设计出的FIR滤波器对加入噪声后的语音信号进行滤波去噪处理，并且分析对比前后时域和频域波形。课题的设计流程图如下图21所示。在WINDOWS下录制一段格式为WAV的语音加入高频噪声信号对语音信号进行频谱分析，绘制出时域和频域图用汉宁窗设计FIR滤波器画出频率响应函数用FIR滤波器对语音信号进行滤波画出语音信号滤波前后频谱波形并进行分析开始结束图21课题设计流程图22语音信号处理221语音信号的采集将话筒输入计算机的语音输入插口上,启动录音机，要求为8000HZ，8位单声道的音频格式，如下图22所示，按下录音按钮,接着对话筒一段话,说完后停止录音,屏幕左侧将显示所录声音的长度。点击放音按钮,可以实现所录音的重现。以文件名“1”保存入CMATLABWORK中。皖西学院本科毕业论文设计第6页图22语音信号的采集222语音信号的时域频谱分析MATLAB软件平台下，利用WAVREAD函数对语音信号进行采样，采集出原始信号波形与频谱，Y,FS,BITSWAVREADBLIP,N1N2，用于读取语音，采样值放在向量Y中，FS表示采样频率HZ，BITS表示采样位数。N1N2表示读取从N1点到N2点的值（若只有一个N的点则表示读取前N点的采样值）。其程序如下Y,FS,BITSWAVREAD1WAVYY,1SIGLENGTHLENGTHYYFFTY,SIGLENGTHPYYYCONJY/SIGLENGTHHALFLENGTHFLOORSIGLENGTH/2FFS0HALFLENGTH/SIGLENGTHFIGUREPLOTF,PYY1HALFLENGTH1XLABELFREQUENCYHZT0SIGLENGTH1/FSFIGUREPLOTT,YXLABELTIMES得到原始语音信号时域波形如下图23所示，频域幅度谱如下图24所示。从图中可以看出语音信号有两个特点在时域内语音信号随着时间的延续而缓慢变化，但在一较短时间内，语音信号基本保持稳定；在频域内语音信号的频谱量主要集中在3003400HZ的范围内，利用这个特点，可以利用一个带通滤波器将此范围内的语音信号频率分量取出，然后按8000HZ的采样频率对语音信号进行采样，就可以取得离散的语音信号。皖西学院本科毕业论文设计第7页图23原始信号时域波形图24原始信号频谱皖西学院本科毕业论文设计第8页223语音信号加噪与频谱分析利用MATLAB程序产生38KHZ的余弦信号噪声加入到语音信号中，模仿语音信号被污染，并对其频谱分析。其主要程序如下FS8000X1WAVREAD1WAVT0LENGTHX11/8000FFS01023/2048AU005DAUCOS2PI3800T噪声为38KHZ的余弦信号X2X1DY1FFTX1,2048Y2FFTX2,2048FIGURE1运行程序后得到加噪后的语言信号波形如下图25所示。图25加噪后的语音信号从上图可以看出，加入38KHZ的噪声信号后，在时域图上与原有信号的时域图有明显差异，在幅度“0”位置处附近多出了高频成分，使加噪后的语音信号显得更加紧凑。PLOTT,X2皖西学院本科毕业论文设计第9页XLABELTIMESYLABEL幅度FIGURE2SUBPLOT2,1,1PLOTF,ABSY111024XLABELHZYLABEL幅度SUBPLOT2,1,2PLOTF,ABSY211024XLABELHZYLABEL幅度SOUNDX2,FS,BITS运行程序后得到原始语音信号和加噪后的语言信号的频谱如下图26所示。图26原始语音信号频谱与加噪后的语音信号频谱比较从上图对比可以看出，加噪后的语音信号表现在频谱图上在38KHZ的位置多出一个高频脉冲成分，表现在回放语音上能听到很刺耳很不舒适的噪音，原有信号听着比较模糊。皖西学院本科毕业论文设计第10页3FIR数字滤波器的设计31数字滤波器基本概念所谓数字滤波器，是指输入、输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的一种数字滤波形式。因此，数字滤波的概念个模拟滤波的相同，只是信号的形式和实现滤波方法不同。由于数字滤波通过数值运算实现滤波，所以其处理精度高、稳定、体积小，重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波功能。32常用窗函数介绍常用的窗函数有矩形窗（RECTANGLEWINDOW）、三角窗（BARTLETTWINDOWS）、汉宁窗（HANNING），又称为升余弦窗、海明窗（HAMMING），又称为改进的升余弦窗，除了以上几种常用窗函数以外，尚有布拉克曼（BLACKMAN）窗、凯塞（KAISER）窗等。对于选择何种窗函数，应充分考虑被处理信号的性质与处理要求等。如果仅要求准确分辨出主瓣频率，而不考虑这种幅值精度，则可选用主瓣宽度较窄从而便于分辨的矩形窗，例如探测物体的自然震动频率等；如果分析的信号是窄带信号，而且具有较强的干扰噪声信号出现的话，则应当选用旁瓣幅度较小的窗函数，如汉宁窗、三角窗等。33FIR数字滤波器概述FIRFINITEIMPULSERESPONSE滤波器又叫有限长单位冲激响应滤波器，是数字信号处理系统中最基本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性，同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。这类滤波器对于脉冲输入信号的响应最终趋向于0，因此而得名。有限脉冲响应滤波器的优点（1）脉冲响应为有限长造成当输入数位信号为有限长的时候，输出数位信号亦为有限序列长度。（2）比无限脉冲响应滤波器较容易实现到最佳状态。（3）线性相位造成HN是偶对称或奇对称且有限长。（4）肯定是稳定的因为Z转换后所有的极点都在单位圆内皖西学院本科毕业论文设计第11页34FIR滤波器的窗函数设计FIR滤波器设计问题在于寻求一系统函数，使其频率响应逼近滤波器要求ZHJWEH的理想频率响应，其对应的单位脉冲响应为。JWDEHNHD用窗函数设计FIR滤波器的的设计思想从时域出发，设计逼近理想。设理HNHD想滤波器的单位脉冲响应为。例如低通线性相位FIR滤波器的理想频率响应与JWDENHD单位脉冲响应分别如下JWNDJWDEHEHNHJNJWDD21NHD一般是无限长的，且是非因果的，不能直接作为FIR滤波器的单位脉冲响应。要想得到一个因果的有限长的滤波器，最直接的方法是截断，即截取NHNWHND为有限长的因果序列，并用合适的窗函数进行加权作为FIR滤波器的单位脉冲响应。按照线性相位滤波器的要求，必须的偶对称的。对称中心应该等于滤波器的时延常数，即21NNHND下面着重介绍用窗函数法设计FIR滤波器的步骤如下1根据对阻带衰减及过渡带的指标要求，选择串窗数类型（矩形窗、三角窗、汉宁窗、哈明窗、凯塞窗等），并估计窗口长度N。先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣的窗函数。然后根据过渡带宽度估计窗口长度N，待求滤波器的过渡带宽度BT近似等于窗函数主瓣宽度，且近似于窗口长度N成反比，NA/BT，A取决于窗口类型。2构造希望逼近的频率响应函数HDEJW。21NJWDGJD皖西学院本科毕业论文设计第12页所谓的“标准窗函数法”，就是选择HDEJW为线性相位理想滤波器，如本论文的低通滤波器，该HDGW应满足CDG013计算HDN。如果给出待求滤波器的截止频率响应函数为HDEJW，那么单位脉冲响应用下式求出WEHNHJNJDD214加窗得到设计结果，验证技术指标是否满足设计要求。D针对该课题用窗函数法设计线性相位FIR数字滤波器的参数如下通带截止频率02CW过渡带宽度40DB具体计算如下1由给定的指标确定窗函数和长度N由于阻带衰减AS40DB，汉明窗和汉宁窗都满足要求，若再考虑从滤波器节数最小的原则出发，这里选用汉宁窗。也可取N21。,208,40W16COS52NRN2确定时延值。0/N3求理想的单位脉冲响应。DENHJNWJWDC21102SIN4求滤波器的单位取样响应HN。NRND210COS50120SI皖西学院本科毕业论文设计第13页35滤波器的编程实现按此要求设计的FIR数字低通滤波器，用MATLAB的程序实现如下DELTW04PIWC02PIAS40NCEIL8PI/DELTW1WINHANNINGNHFIR1N1,WC/PI,WINOMEGALINSPACE0,PI,512MAGFREQZH,1,OMEGAMAGDB20LOG10ABSMAGPLOTOMEGA/PI,MAGDBAXIS011000GRIDXLABEL归一化频率YLABEL幅度/DB此低通滤波器图像如下图31所示。图31低通滤波器的归一化频率皖西学院本科毕业论文设计第14页36用滤波器对加噪语音信号进行滤波上一节利用窗函数法按照论文要求设计了FIR低通滤波器并绘图，观察所设计的滤波器是否能够对采集的一部分语音信号进行相关处理，并将滤波前后的时域波形进行相比较，且对其快速傅里叶变换，即XFFTSIGNAL，其目的是对比前后的频域频谱波形，分析所设计的滤波器能否达到设计要求。在MATLAB程序设计中，FIR滤波器则是利用函数FFTFILT对语音信号进行滤波。程序如下Y,FS,NBITSWAVREAD1WAVNLENGTHYYFFTY,NSOUNDYFIGURE4PLOTYFIGURE5PLOTABSYFP1200FS1100FT8000AS20AP1WP2PIFP/FTWS2PIFS/FTFP2FTTANWP/2FS2FTTANWS/2N,WNBUTTORDWP,WS,AP,AS,SB,ABUTTERN,WN,SNUM,DENBILINEARB,A,1H,WFREQZNUM,DENFIGURE2DAUCOS2PI3800TX2X1DY1FFTX1,2048Y2FFTX2,2048FIGURE3PLOTW800005/PI,ABSHZFILTERNUM,DEN,YSOUNDZMZFIGURE1SUBPLOT2,2,3PLOTABSM,RGRIDSUBPLOT2,2,4PLOTZ,BGRIDSUBPLOT2,2,2PLOTY,BGRIDSUBPLOT2,2,1PLOTABSY,RGRID皖西学院本科毕业论文设计第15页图32频率特性曲线图33滤波前后语音信号频谱的比较皖西学院本科毕业论文设计第16页图34滤波前后的信号波形比较由图32频率特性曲线可得知，该低通滤波器滤除1200HZ以上的高频信号，保留01200HZ以内的低频语音信号，符合设计滤除高频噪音保留低频原始语音信号的特点。由上图33可知，掺有高频噪音的信号经过所设计的低通滤波器之后，保留了原始的低频语音信号，滤除了掺在其中的高频信号，使语音信号听着没有那么的尖锐刺耳了，这说明已经达到了滤除高频噪音信号的目标。分析图34滤波前后的信号波形比较得知，滤波前后语音信号的波形发生了明显改变，滤波前后的信号密度明显减小，这是滤除了高频噪音保留了低频语音信号的结果。37回放语音信号语音信号经过FIR滤波器的滤除噪声的处理，在MATLAB中，函数SOUND可以对声音进行回放。其调用格式SOUNDX，FS，BITS；可以听出来滤波前后的声音所发生的变化，而且声音变得没有加噪时那么刺耳了，比原始语音更加地平滑。用汉宁窗设计FIR滤波器滤掉了在语音中加入的高频的噪声，而且也把原始语音的很小的一部分也滤掉了，余下的语音信号全都是低频语音信号，所以回放语音的时候听起来比以前的更加平滑，说明设计的低通滤波器是成功的。皖西学院本科毕业论文设计第17页4结论人们在获取语音信号的过程中，不可避免的会受到外界环境的干扰和影响，这些干扰和影响不但降低了语音质量和语音的可懂度，严重时将导致不可预知的不良效果。当这些语音信号经过数字处理，用滤波器把噪声杂波滤除，便得到纯净的语音信号。该毕业设计的目的就是设计一个数字滤波器过滤被干扰的语音信号，论文的主要工作有录制一段原始语音信号，然后利用MATLAB程序函数产生一定的噪声信号，使语音信号显得更加浑浊，调用MATLAB程序中SOUND函数功能试听一下语音信号的变化，会发现声音比较刺耳尖锐，同时通过信号的波形与频谱图可以明显看出加入噪声