基于LMS算法的多麦克风降噪模板

1、系主任(或责任教师)签名:2012课程设计任务书指导教师：艾青松工作单位：信息工程学院题目：基于LMS算法的多麦克风降噪初始条件：MATLA颤件数字信号处理基础知识要求完成的主要任务：给定主麦克风录制的受噪声污染的语音信号和参考麦克风录制的噪声，实现语音增强的目标，得到清晰的语音信号。(1)阅读参考资料和文献，明晰算法的计算过程，理解LMSB法基本过程；(2)主麦克风录制的语音信号是LMSprimsp.wav,参考麦克风录制的参考噪声是LMSrefns.wav.用matlab指令读取；(3)根据算法编写相应的MATLA勰序；(4)算法仿真收敛以后，得到增强的语音信号；(5)用matlab指令回

2、放增强后的语音信号；(6)分别对增强前后的语音信号作频谱分析。指导教师签名：2012年月日学生姓名:柯一凡专业班级:通信1102武汉理工大学信息处理课群课程设计报告摘要随着社会工业生产的不断进步，各种噪声污染越来越严重。目前普遍采用的模拟降噪方法已不能满足要求，未来的研究将朝着以数字信号处理器及相关算法为技术支撑的数字降噪技术发展。滤波器设计在数字信号处理中占有极其重要的地位，Matlab功能强大、简单易学、编程效率高，深受广大科技工作者的欢迎。特别是Matlab还具有信号分析工具箱，不需具备很强的编程能力，就可以很方便地进行信号分析、处理和设计。利用MATLAB信号处理工具箱可以快速有效地设

3、计各种数字滤波器。在过去的几十年中，基于多麦克风的噪声消除问题一直是人们关注的课题。而在众多算法中，基于LMS算法更新滤波器权值的广义旁瓣消除器结构应用最为广泛。多麦克风降噪设计主要是通过自适应滤波器来实现的。文中采用LMS算法在MATLAB中实现了自适应滤波器的设计与实现。在MATLAB中建立了数字降噪系统模型，并且针对该模型利用MATLAB语言进行编程，仿真结果表明此设计实现了对信号中混有的环境噪声进行降噪，并且效果远远高于模拟降噪技术。关键词：MATLAB,语音增强，LMS算法，多麦克风武汉理工大学信息处理课群课程设计报告AbstractAbstractWiththesocialprog

4、ressofindustrialproduction,allkindsofnoisepollutionismoreandmoreserious.Thecurrentwidespreaduseofthesimulationofthenoisereductionmethodcannotmeettherequirementhas,andfurtherresearchwilltowardtodigitalsignalprocessorandrelatedalgorithmfortechnicalsupportdigitalnoisereductiontechnologydevelopment.Filt

5、erdesignindigitalsignalprocessingplaysanextremelyimportantrole,.Matlabispowerful,easytolearn,programmingefficiency,whichwaswelcomedbythemajorityofscientists.Matlabalsohasaparticularsignalanalysistoolbox,itneednothavestrongprogrammingskillscanbeeasilysignalanalysis,processinganddesign.UsingMATLABSign

6、alProcessingToolboxcanquicklyandefficientlydesignavarietyofdigitalfilters.Inthepastfewdecades,basedonthenoiseofthemicrophoneisalwaysthepeopletopayattentiontoeliminatethesubject.Andinnumerousalgorithm,LMSalgorithmbasedonthegeneralizedupdatefilterweightvalueside-lobeeliminatemostwidelyusedthestructure

7、.Morenoisereductiondesignmainlythroughthemicrophonetorealizeadaptivefilter.ThisarticleapplyLMSalgorithminMATLABrealizeadaptivefilter,thedesignandimplementation.InMATLABestablisheddigitalnoisereductionsystemmodel,andthemodelforuseofMATLABlanguageprogramming,thesimulationresultsshowthatthedesignedtosi

8、gnalmixinsomeenvironmentalnoisenoisereduction,andtheeffectisfarhigherthananalognoisereductiontechnology.Keywords：MATLAB,Speechenhancement,LMSalgorithm,microphoneII武汉理工大学信息处理课群课程设计报告目录摘要I.AbstractII1前言1.2原始语音信号采集与处理2.2.1语号的采样理论依据22.1.1采样频率采样位数采样定理2.2.2语音信号的采集3.2.3语音信号的时频分析32.4语音信号加噪与频谱分

9、析.5.3基于LMS自适应滤波器的设计7.3.1基本LMS算法7.3.2自适应噪声抵消原理113.3基于最小均方误差准则（LMS）的自适应噪声抵消134程序及运行结果1.54.1程序154.2运行结果1.74.3调试分析1.95心得体会20参考文献21武汉理工大学信息处理课群课程设计报告1前言多麦克风降噪是降噪技术的一个重要应用。我国的降噪技术研究始于80年代初期，采用的手段主要有三种，其中的动态降噪技术（DNR）又可以分为模拟动态降噪技术和数字动态降噪技术。目前国内外解决噪声问题最普遍的方法是采用模拟动态降噪技术，数字降噪技术的研究尚处于初期阶段。数字降噪技术比模拟降噪技术具有更大的优点。模

10、拟降噪技术全采用硬件实施，修改和调试十分困难，对元器件参数的变化也很敏感，技术指标受元器件的误差影响较大，降噪效果不稳定，不利于产品的批量生产。而数字降噪技术由于采用计算机技术实现自适应滤波，通过修改软件算法就可以达到不同的降噪效果，不用更改硬件结构，调试和维修都非常方便；数字降噪技术采用自适应滤波技术，可以实时跟踪噪声的变化进一步进行处理，因此降噪效果较好。另外，数字降噪技术抗干扰能力强，本身具有自恢复能力，并且在整个音频带内降噪比较均衡，而模拟降噪技术偏重于低频段，高频段效果较差。因此降噪技术未来的发展方向是数字降噪技术，以数字信号处理（DSP）及其相关算法为技术支撑的数字降噪技术代表着当

11、今降噪技术的发展。目前市场上的麦克风降噪产品主要是模拟降噪，因此数字降噪的设计在国内属于领先技术。多麦克风数字降噪的系统原理是通过麦克风装置直接检测出噪声信号和音频信号的混合信号，然后将混合信号通过DSP数字降噪模块进行噪声分离并产生降噪信号来抵消噪声，因此人耳就可以只听到较纯净的音频信号而不受环境噪声的干扰。本文采用最小均方误差（LMS）算法，实现了数字降噪DSP中消除噪声的模块自适应滤波器的设计，介绍了其在MATLAB中编程及仿真输出，并通过程序实现了设计。武汉理工大学信息处理课群课程设计报告2原始语音信号采集与处理2.1语音信号的采样理论依据2.1.1采样频率采样频率是指计算机每秒钟采集

12、多少个声音样本，是描述声音文件的音质、音调，衡量声卡、声音文件的质量标准。采样频率越高，即采样的间隔时间越短，则在单位时间内计算机得到的声音样本数据就越多，对声音波形的表示也越精确。采样频率与声音频率之间有一定的关系，根据奎斯特理论，只有采样频率高于声音信号最高频率的两倍时，才能把数字信号表示的声音还原成为原来的声音。这就是说采样频率是衡量声卡采集、记录和还原声音文件的质量标准。2.1.2采样位数采样位数即采样值或取样值，用来衡量声音波动变化的参数，是指声卡在采集和播放声音文件时所使用数字声音信号的二进制位数。采样频率是指录音设备在一秒钟内对声音信号的采样次数，采样频率越高声音的还原就越真实越

13、自然。采样位数和采样率对于音频接口来说是最为重要的两个指标，也是选择音频接口的两个重要标准。无论采样频率如何，理论上来说采样的位数决定了音频数据最大的力度范围。每增加一个采样位数相当于力度范围增加了6dBo采样位数越多则捕捉到的信号越精确。对于采样率来说你可以想象它类似于一个照相机，44.1kHz意味着音频流进入计算机时计算机每秒会对其拍照达441000次。显然采样率越高，计算机摄取的图片越多，对于原始音频的还原也越加精确。2.1.3采样定理在进行模拟/数字信号的转换过程中，当采样频率fsmax大于信号最高频率fmax的2倍时，即：fsmax=2fmax,则采样之后的数字信号完整地保留了原始信

14、号中的信息，一般实际应用中保证采样频率为信号最高频率的510倍，采样定武汉理工大学信息处理课群课程设计报告理又称奈奎斯特定理。1924年奈奎斯特(Nyquist)就推导出：在理想低通信道的最大码元传输速率=2W*10g2N(其中W是理想低通信道的带宽,N是电平强度)。2.2语音信号的采集利用PC机上的声卡和WINDOWS操作系统可以进行数字信号的采集。 将话筒输入计算机的语音输入插口上，启动录音机。按下录音按钮，接着对话筒说话语音信号处理”说完后停止录音，屏幕左侧将显示所录声音的长度。点击放音按钮，可以实现所录音的重现。以文件名OriSound”保存入c:MATLABwork中。 可以看到，

15、文件存储器的后缀默认为.wav， 这是WINDOWS操作系统规定的声音文件存的标准。2.3语音信号的时频分析在MATLAB软件平台下，利用wavread函数对语音信号进行采样，记住采样频率和采样点数，Wavread函数调用格式：y=wavread(file)%读取file所规定的wav文件，返回采样值放在向量y中y,fs,bits=wavread(file)%采样值放在向量y中，fs表示采样频率(hz),bits表示采样位数y=wavread(file,N1,N2)%读取从N1到N2点的采样值放在向量y中对语音信号OriSound.wav进行采样其程序如下：y,fs,nbits=wavread

16、(OriSo谛音信万如载入MATLAB仿真软件平台中画出语音信号的时域波形，再对语音信号进行频谱分析。MATLAB提供了快速傅里叶变换算法FFT计算DFT的函数fft,其调用格式如下：Xk=fft(xn,N)参数xn为被变换的时域序列向量，N是DFT变换区间长度，当N大于xn的长度时，fft函数自动在xn后面补零。，当N小于xn的长度时，fft函数计算xn的武汉理工大学信息处理课群课程设计报告y=wavread(file,N)%读取前N点的采样值放在向量y中(式2-3-1)前N个元素，忽略其后面的元素。在本次设计中，我们利用fft对语音信号进行快速傅里叶变换，就可以得到信号的频谱特性。程序如下

17、：y,fs,nbits=wavread(OriSound)%语音信号的采集sound(y,fs,nbits);%语音信号的播放n=length(y);Y=fft(y,n);%快速傅里叶变换figure;subplot(2,1,1);plot(y);title(原始信号波形，fontweight,bold);axis(7800080000-11);grid;subplot(2,1,2);plot(abs(Y);title(原始信号频谱，fontweight,bold);axis(015000004000);4武汉理工大学信息处理课群课程设计报告2.4语音信号加噪与频谱分析在MATLAB中产生高斯

18、白噪声非常方便，我们可以直接应用两个函数：一个是WGN,另一个是AWGNoWGN用于产生高斯白噪声，AWGN则用于在某一信号中加入高斯白噪声。也可直接用randn函数产生高斯分布序列。在本次设计中， 我们是利用MATLAB中的随机函数(rand或randn)产生噪声加入到语音信号中，模仿语音信号被污染，并对其频谱分析。Randn函数有两种基本调用格式：Randn(n和Randn(m,n),前者产生n刈服从标准高斯分布的随机数矩阵，后者产生mKn的随机数矩阵。在这里，我们选用Randn(m,n)函数c语音信号添加噪声及其频谱分析的主要程序如下：y,fs,nbits=wavread(OriSoun

19、d);sound(y,fs,nbits);n=length(y);Noise=0.2*randn(n,2);s=y+Noise;sound(s);figure;subplot(2,1,1);plot(s);title(加噪语音信号的时域波形,fontweight,bold);axis(7800080000-11);grid;S=fft(s);subplot(2,1,2);plot(abs(S);title(加噪语音信号的频域波形,fontweight,bold);axis(015000004000);grid;武汉理工大学信息处理课群课程设计报告程序结果如下图:加噪语音信号的时域波形x10加噪

20、语音信号的频域波形x10图2-3信号加噪时域波形图与频谱图武汉理工大学信息处理课群课程设计报告3 3 基于 LMSLMS 自适应滤波器的设计在实际应用中，常常无法得到信号和噪声统计特性的先验知识。在这种情况下，自适应滤波技术能够获得极佳的滤波性能，因而具有很好的应用价值。常用的自适应滤波技术有：最小均方(LMS)自适应滤波器、递推最小二乘(RLS)滤波器、格型滤波器和无限冲激响应(IIR)滤波器等。这些自适应滤波技术的应用又包括：自适应噪声抵消、自适应谱线增强和陷波等。LMS自适应滤波器是使滤波器的输出信号与期望响应之间的误差的均方值为最小，因此称为最小均方(LMS)自适应滤波器。3.1基本L

21、MSi：法构成自适应数字滤波器的基本部件是自适应线性组合器，如图3-1的所示,设线性组合器的M个输入为x(k-1)-,x(k-M),其输出y(k)是这些输入加权后的线性组合，即（式3-1-1）定义权向量,W=W1,W2,W3；,WmT,且X(k)=X(k-1)T,X(kM)T)T在图3-1中，令d(k)代表所期望的响应”，并定义误差信号(k)=d(k)y(k)=d(k)-WiX(k-i)i=e武汉理工大学信息处理课群课程设计报告（式3-1-2）（式3-1-3）y(k)=、WiX(k-i)i0图3-1臼适应线性组合器式(3-1-3)写成向量形式e(k)=d(k)WTX(k)=d(k)-XT(k)

22、W(式3-1-4)误差平方为/(k)=d2(k)-2d(k)XT(k)W+WTX(k)XT(k)W上式两边取数学期望后，得均方误差E2(k)=Ed2(k)-2Ed(k)XT(k)W+WTEX(k)XT(k)W(式3-1-5)定义互相关函数行向量Rd：RTd=Ed(k)XT(k)(式3-1-6)和自相关函数矩阵RXX=EX(k)XT(k)(式3-1-7)则均方误差(3-1-5)式可表述为E2(k)=Ed2(k)-2RTdW+WTRXXW(式3-1-8)这表明，均方误差是权系数向量W的二次函数，它是一个中间向上凹的抛物形曲面，是具有唯一最小值的函数。调节权系数使均方误差为最小，相当于沿抛物形曲面下

23、降找最小值。可以用梯度来求该最小值。将式(3-1-8)对权系数W求导数，得到均方误差函数的梯度V(k)=-2RXd+2RxxW(式3-1-9)令V(k)=0,即可求出最佳权系数向量Wopt=RXXRxd(式3-1-10)它恰好是研究Wiener滤波器遇到过的Wiener-Hopf方程。 因此， 最佳权系数向量通常也叫作Wiener权系数向量。将Wopt代入式(3-1-8)得最小均方误差Ee2(k)min=Ed2(k)-RTdWopt(式3-1-11)利用式(3-1-10)求最佳权系数向量的精确解需要知道RXX和Rxd的先验统计知识，而且还需要进行矩阵求逆等运算。WidrowandHoff(19

24、60)提出了一种在这些先验统计知识未知时求Wopt的近似值的方法， 习惯上称为WidrowandHoffLMS算法。 这种算法的根据是最优化方法中的最速下降法。根据最速下降法，“下一时刻”权系数向量W(k+1)应该等于“现时刻”权系数向量W(k)加上一个负均方误差梯度-V(k)的比例项，即武汉理工大学信息处理课群课程设计报告W(k+1)=W(k)-典(k)(式3-1-12)式中，小是一个控制收敛速度与稳定性的常数，称之为收敛因子。不难看出，LMS算法有两个关键：梯度V(k)的计算以及收敛因子小的选择。(一)V(k)的近似计算精确计算梯度V(k)是十分困难的，一种粗略的但是却十分有效的计算V(k

25、)的近似方法是：直接取e2(k)作为均方误差E/(k)的估计值，即V(k)=V2(k)=2(RI(k)(式3-1-13)得到梯度估值V(k)=-2(k)X(k)于是，Widrow-HoffLMS算法最终为W(k+1)=W(k)+2(k)X(k)(式3-1-14)式(3-1-14)的实现方框图如图3-2所示图3-2LMS算法的实现方框图AA卜面分析梯度估值(k)的无偏性。(k)的数学期望为=-2以工火的一=2(旦 a火口.印(幻=(曾(式3-1-15)在上面的推导过程中，利用了d(k)和e(小者皆为标量的事实。在得到最后的结果时，利用了式(3-1-9)。式(3-1-15)表明，,剃雯估值)是无偏

26、估(二)一勺选择对权系数向量更新公式(3-1-14)两边取数学期望，得EW(k+1)=EW(k)+2EE(k)X(k)武汉理工大学信息处理课群课程设计报告=(I-2Rxx)EW(k)+2北(式3-1-16)式中，I为单位矩阵，Rxd=Ed(k)X(k)和Rxx=EX(k)XT(k)。当时，k=0时，EW(1)=(I-2仙Rxx)EW(0)+2Xd对于k=1,利用上式结果，则有EW(2)=(I-2Rxx)EW(1)+2Xd(I-2甲xx)2EW(0)+2立(I-2用xx)iRxdiO起始时，EW(0)=W(0)故重复以上迭代至k+1,则有kEW(k+1)=(I-2甲xx)k，W(0)+2迂(I-

27、2甲xx)iRxd(式3-1-17)i=0由于Rxx是实值的对称阵，我们可以写出其特征值分解式T1Rxx=Q2Q=mQ(式3-1-18)这里，我们利用了正定阵Q的性质QT=Q,，且Adaglk左)是对角阵，其对角元素人是R的特征值。将式(3-1-18)代入式(3-1-19)后得k1k1.1i,、EW(k+1)=(I-2pQ2Q)W+2Z(I-2NQ2Q)Rxd(式3-1-19)i=0注意到以下恒等式及关系式：(3)假定所有的对角元素的值均小于1(这可以通过适当选择以实现)，则ljm(I-2=011-1(4)Rxx=Q2Q将上式代入式(8-1-19),结果有1.11EW(k+1)=Q2QRxd=

28、RXXRxd=Wopt(式3-1-20)由此可见，当迭代次数无限增加时，权系数向量的数学期望值可收敛至Wiener解，其条件是对角阵(I-2仙洲所有对角元素均小于1,即武汉理工大学信息处理课群课程设计报告(1)(I-2pQ3Q,)i=Q(I-2p*Q,limkk一_1i_._1_1Z(I-2qiQ)=ZQ(2NEQi=0if01/一(式3-1-21)max10其中入ma狙RXX的最大特征值。小称为收敛因子，它决定达到式(3-1-20)的速率。事实上，W(k)收敛于W0Pt由比值d=入max/入min定，该比值叫做谱动态范围。大的d值喻示要花费很长的时间才会收敛到最佳权值。克服这一困难的方法之一

29、是产生正交数据。基本LMS自适应算法如下：初始化：W(0)=0;R(0)=I;选择工0-maxFork=1tonfinaldo:W(k)=W(k-1)+2叱x(k-WT(k-1)X(k)X(k)LMS自适应滤波器如图3-3所示：图3-3LMS自适应滤波器3.2自适应噪声抵消原理自适应噪声抵消的目的是要去除主信号中的背景噪声。主信号由有用信号和背景噪声组成，而背景噪声与参考信号中的噪声相关。因此，自适应噪声抵消技术主要依赖于从主信号和噪声中获取参考信号。Widrow和Hoff发展了最小均方误差(LMS)自适应算法和称为自适应线性阈值逻辑单元(ADALINE)的模式识别方法。1965年，基于最小均

30、方误差准则(LMS)的自适应噪声抵消首次得以实现，随后，自适应噪声抵消在信号处理、地11武汉理工大学信息处理课群课程设计报告震和生物医学领域均获得成功应用基于维纳理论的自适应噪声抵消需要无限加权滤波器，以极小化输出误差为了实现维纳滤波方案，必须使用有限加权滤波器。换句话说，自适应滤波器必须假定维纳滤波器是一个有限冲激响应(FIR)滤波器(a)最佳噪声抵消器自适应噪声抵消器图3-4自适应噪声抵消原理方框图如图3-4(a)所示是基于维纳滤波器的自适应噪声抵消原理方框图。主信号由有用信号x(n)和背景噪声v(n)构成，其中s(n)和v(n)不相关。参考信号r(n)可与s(n)AA或v(n)相关。v(

31、n)是背景噪声的最佳估计。v(n)可以通过选择最佳FIR维纳滤波器的最佳加权w(n)计算得出，即AMv(n)=Wwm(n)r(n-m)0M(式3-2-1)i=0其中，M表示滤波器的阶；r(n-m)由延时获得。具有M个权重滤波器的估计误差e(n)由下式定义：Ae(n)=x(n)-v(n)=x(n)-wT(n)r(n)(式3-2-2)由正交原理有，e(n)和r(n)正交。对式(322)两边取平方和数学期望，可得Ee(n)2=Ex(n)2-2PTw+wTRw(式3-2-3)其中，输入信号s(n)和参考矢量r(n)之间的互相关用P表示，即P=Ex(n)匚(n)T(式8-2-4)12武汉理工大学信息处理

32、课群课程设计报告R表示输入自相关矩阵，即R=r(n)r(n)T(式3-2-5)令均方估计误差函数的梯度等于0,可得最佳FIR滤波器(维纳滤波器)权重如下，w=R-P(式3-2-6)实际上，通常P和R的统计量是未知的。然而，用Widrow和Hoff提出的方法迭代求解式(3-2-6)能够克服这一限制。如果参考信号和主信号中的噪声相关，则自适应滤波器将在输出端去除其相关性，具体方法是：从参考信道的噪声中产生一个主信道中背景噪声的估计值v(n),然后从主信道中减去这个估计噪AA声v(n),那么自适应滤波器的输出就是有用信号的估计s(n)0用最速下降法(或梯度下降法)可得到式(3-2-6)的解。 自适应

33、滤波器的加权值W被更新的第(n+1)步迭代式为f(w(n)一w(n+1)=w(n)+Aw(n)=w(n)-N(式3-2-7)二W3.3基于最小均方误差准则(LMS)的自适应噪声抵消如果没有关于参考信号向量P和输入自相关矩阵R的先验信息， 要实现最优滤波器加权是不可能的。因此，Widrow和Hoff提出了另一种可迭代的维纳FIR滤波实现方法。在这种方法中，滤波器的权重可被更新为：w(n+1)=w(n)+2ie(n)(n)(式3-3-1)综上所述，基于最小均方误差准则(LMS)的自适应噪声抵消算法可按以下步骤实现：第一步：设一个初值；wm(0);A第二步：计算自适应FIR滤波器的输出v(n),AM

34、v(n)=Zwm(n)r(n-m)(式3-3-2)i0其中，M表示滤波器的阶第三步：估计当前时刻n的误差e(n),13武汉理工大学信息处理课群课程设计报告e(n)=x(n)-v(n)等(n)(式3-3-3)第四步：用最速下降LMS算法更新滤波器权重：wm(n):Wm(n+1)=Wm(n)+2ie(n)r(m)0mn+1,到下一个时刻，重复以上步骤，直至满足要求为止。收敛参数以必须是正数，并且满足：0仙工(式3-3-5)R0小(式3-3-6)max其中，入maxft示自相关矩阵R的最大特征值。然而，在实际应用中，R的具体值是不知道的，参数小的值也需要试探性地选择。若以取值小，能保证收敛，但需要注

35、意的是，如果取得过小，收敛速度将非常慢；相反，若以取值大，可以提高收敛速度，却是以噪声收敛为代价的。如果参考输入信号r(n)是频率为的正弦信号，自适应滤波器将从主信号中滤除所有的频率为的正弦成分。 基于最小均方误差准则(LMS)的自适应噪声抵消算法的程序见下节。14z=primary(i)-w*(fref(i-Worder+1:i);n2=fref(i-Worder+1:i);w=w+2*u*n2*z;output(i-Worder)=z;end;figure(1);%作图plot(output);%画降噪后的语音波形title(output);wavwrite(output,fs,outyu

36、yin);out,fs,bits=wavread(outyuyin.wav);15武汉理工大学信息处理课群课程设计报告4 4 程序及运行结果4.1程序fs=20000;bits=32;primary,fs,bits=wavread(F:LMSprimsp.wav);%读入主麦克风语音信号sound(primary,fs,bits);fref,fs,bits=wavread(F:LMSrefns.wav);%读入参考语音信号sound(fref,fs,bits);%初始化primary=primary;fref=fref;Worder=10;%滤波器阶数u=0.005;w=zeros(1,Wor

37、der);output=primary;%主语音输出loopsize=max(size(primary);fori=1+Worder:loopsize%写LMS算法公式%生成降噪后的语音武汉理工大学信息处理课群课程设计报告sound(out,fs,bits);figure(2);plot(primary);title(primaryinput);figure(3);plot(fref);title(frefnoise);y1=fft(primary,4096);y2=fft(output,4096);figure(4)subplot(1,2,1);plot(abs(y1);title(输入信号

38、的频谱)axis(0200,040);subplot(1,2,2);plot(abs(y2);title(输出信号的频谱)axis(0200,040)16%画主麦克风语音波形%画参考麦克风语音波形%画输入输出的频谱图武汉理工大学信息处理课群课程设计报告4.2运行结果图4.1主麦克风语音信号17武汉理工大学信息处理课群课程设计报告图4.2参考麦克风语音信号图4.3输出信号18武汉理工大学信息处理课群课程设计报告图4.4输入信号和输出信号的频谱4.3调试分析结果分析：1)通过对滤波前后主麦克风语音信号波形的分析可知，滤波前信号波形受噪声干扰严重，有用信号被噪声覆盖，几乎都是噪声；滤波后滤除了信号中

39、的噪声分量，有用信号被提取出来，起到了信号增强的目的。2)通过对滤波前后主麦克风语音信号频谱的观察分析可知，滤波前后频谱分量明显减弱，滤波后频谱分量较集中，有明显改善。经过以上的降噪处理后，可在Matlab中用函数sound对声音进行回放。其调用格式：sound(primary,fs,bits)和sound(out,fs,bits)。可以察觉滤波前后的声音有明显的变化。19武汉理工大学信息处理课群课程设计报告5 5 心得体会本次课程设计通过给定主麦克风录制的受噪声污染的语音信号和参考麦克风录制的噪声， 实现语音增强的目标，得到清晰的语音信号。运用数字信号处理理论，基于LMST法设计自适应滤波器，并对被噪声污染的语音信号进行滤波，分析滤波前后信号的时域和频域特征，回放语音信号。这次课设对我来说有一定的难度，之前学习过数字信号处理的基础知识，会一些基本的MATLA的句，但仅仅凭借这些知识无法完成本次课设的要求任务。