《基于LMS最小均方误差法的语音降噪-(DOC 19页)》由会员分享，可在线阅读，更多相关《基于LMS最小均方误差法的语音降噪-(DOC 19页)（21页珍藏版）》请在金锄头文库上搜索。
1、信号处理与分析课程设计报告 第 19 页 共 19 页信号处理与分析课程设计报告 项目名称：基于LMS最小均方误差法的语音降噪姓名：07021102 台斯瑶 07021106 王金泊指导教师：李如玮目录一、课题背景和简介2二、训练目的3三、训练内容3四、最小均方差LMS实现自适应滤波器的方法介绍3五、实验设计及实施过程41、滤波器结构设计52、高斯白噪声的实现方法63、LMS算法的实现8六、实验结果分析9七、从实验中分析LMS算法的缺点11八、实验完整程序12九、参考文献15十、特别鸣谢16一、课题背景和简介本课题是根据电子信息类本科生信号处理和分析课程的学习内容和语音信号处理的实际应用相结合
2、而设计的实践性训练。课程训练以数字信号处理为基础，在掌握基本原理的同时，理解语音信号的相关知识并结合实际应用实现对语音信号的分析和处理。滤波是一种数字信号处理操作，其目的是为了处理某个信号，以便提取出输入信号中所包含的期望信息。在数字信号处理课程中，我们基本掌握了一些线性滤波器的设计方法，有固定的规范可遵循。而在我们的实际生活中，充满了偶然和随机，时不变滤波器已不能够满足更好效果的滤波。在这种情况下，我们就需要自适应滤波器。可以看到，随着数字超大规模集成技术的发展，自适应滤波技术在很多领域得到了广泛应用。最小均方算法是一种搜索算法，他通过对目标函数进行适当的调整，简化了对梯度相量的计算。由于其
3、计算简单性，LMS算法以及其它一些相关算法已广泛应用于自适应滤波的各种应用中。而LMS算法是自适应滤波理论中应用最广泛的算法。这主要归因于其地计算复杂度、在平稳环境中的收敛性、其均值无偏地收敛到维纳解以及利用有限精度算法实现时的稳定特性等等。对LMS最小均方算法的学习，将加深我们对数字信号处理课程的理解，同时对我们今后滤波技术的应用奠定了巩固的基础。二、训练目的1.通过利用c程序实现数字信号处理的相关功能，巩固对信号处理原理知识的理解，提高实际编程和处理数据的综合能力，初步培养在解决信号处理实际应用问题中所应具备的基本素质和要求。2.培养研发能力，通过设计实现不同的信号处理问题，初步掌握在给定
4、条件和功能的情况下，实现合理设计算法结构的能力。3.提高文献整理和资料查询的能力。通过课下对相关语音知识的学习和理解，培养快速解决实际问题的能力，并在文献整理的过程中学会科技文献的写作，提高语言表达能力。三、训练内容根据语音信号的特点，利用不同信噪比的高斯白噪声对语音进行加噪，利用LMS最小均方误差法设计实现自适应滤波器，并讨论语音状态变化下的收敛情况。四、最小均方差LMS实现自适应滤波器的方法介绍自适应滤波器是符合某种准则的最佳滤波器。它是通过对环境进行学习，逐渐达到或逼近最优滤波器。由于学习过程中有“导师”存在，因此它是一种具有监督学习功能的过程。当滤波器的应用环境发生缓慢变化时，相当于滤
5、波器应用于非平稳环境，但环境变化比学习速度更缓慢时，自适应滤波器能够自适应地跟踪这种非平稳变化。开始时，给FIR滤波器赋予任意的初始权系数，在每个时刻，用当前权系数对输入信号进行滤波运算，产生输出信号，输出信号与期望响应的差定义为误差信号，由误差信号与输入信号矢量一起构造一个校正量，自适应地调整权矢量，使误差信号趋于降低的趋势，从而使滤波器逐渐达到或接近最优。LMS最小均方误差的方法是由最速下降法推导而出。最速下降法需要求出其梯度的精确值，要求输入信号和期望信号平稳，且 (Rak=抽头输入向量u(n)与期望响应d(n)的互相关向量；Rxx=抽头输入向量u(n)的相关矩阵；W=抽头权向量)要首先
6、估计 和 ，这给具体实现带来很大困难，因此该算法还不是真正意义的自适应滤波算法，但讨论最陡下降法是有意义，由最陡下降法可以很直观地导出一类自适应滤波算法 - LMS算法。LMS算法的基本思想：调整滤波器自身参数，使滤波器能够自适应地跟踪这种输入信号的变化，实现最优滤波。当横向滤波器运行在实数据的情况下，该算法大体上可描述为：抽头权向量更新值=老的抽头权向量值+学习速率参数*抽头输入向量*误差信号其中误差信号定义为期望向量与抽头输入向量所产生的横向滤波器的实际向量之差设输入信号为u（n）,LMS算法理论推导过程如下：滤波器输出y(n)为： （1）由误差定义得： （2）使用最小均方误差法，得代价函
7、数为均方误差为： （3）式(3)中J 是滤波器的系数k w (k = 0,1,2,)的函数。通过选择最优的系数，使J 达到最小值。定义梯度向量为J ， （4）另外，最陡下降迭代方程为：w(n +1) = w(n) J (n) （5）LMS 是直接利用单次采样数据获得的e2 (n)来代替均方误差J (n)，从而进行梯度估计，每次迭代时计梯度估计为： （6）式（6）代入式（5），得到系数向量自适应迭代法：w(n +1) = w(n) J (n) = w(n) + 2e(n)u(n) （7）式（7）称最小均方自适应算法LMS。五、实验设计及实施过程 自适应滤波器的基本原理：所谓自适应滤波，就是利用前
8、一时刻已经获得的滤波器参数的结果，自动的调节显示课的滤波器的参数，以适应信号和噪声未知的或随时间变化的统计特性，从而实现最优滤波。自适应滤波器实质上就是一种能调节器自身传输特性已达到最优的维纳滤波器。自适应滤波器不需要关于输入信号的先验知识，计算量小，特别是用于实时处理。 由于无法预知信号和噪声的特性，必须设计自适应滤波器，以跟踪信号和噪声的变化。自适应滤波器的特性变化是由自适应算法通过调整滤波器系数来实现的。一般而言，自适应滤波器有两部分组成，已是滤波器的结构，二是调整滤波器系数的自适应算法。本实验的滤波结构采用FIR滤波器设计，自适应算法采用LMS最小均方差算法。总体结构图如图1图 11、
9、滤波器结构设计FIR滤波器原理FIR滤波器的数学表达式为 (8)式中：N为FIR滤波器的抽头数；x(n)为第n时刻的输入样本；h(i)为FIR滤波器第i级抽头系数。在本实验中，即为w(n);其相应的z变换为： （9） 式中：z-i为N-1阶多项式。普通的直接型FIR滤波器结构如图2所示。图2 FIR滤波器直线型结构 在本实验中，抽头系数会根据每一次的输出进行自适应修改。程序实现对输入信号xi和滤波系数wi进行时域的卷积和，求出输出信号y.for(i=0;iBUFLEN;i+) y=y+xi*wi; /实际输出信号的合成 2、高斯白噪声的实现方法高斯白噪声的原理高斯白噪声：具有高斯分布的噪声就成
10、为高斯噪声，而它的功率谱密度又是均匀分布的，则称它为高斯白噪声。在本实验中，我组针对高斯白噪声的特点进行了程序的设计。首先确定的我们需要的白噪声的期望值、方差、和噪声数组的长度，通过一个专门产生白噪声的种子文件进行设计，利用如下公式，最终得到了所需要的白噪声信号。根据定理：设为n 个相互独立的均匀分布的随机数，其期望为，方差为，有中心极限定理可知（如果n个独立随机变量的分布式相同的，并且具有有限的数学期望和方差，当n无穷大时，它们之和的分布趋近于高斯分布；即使n个独立随机变量不是相同分布的，当n无穷大时，如果满足任意一个随机变量都不占优或对和的影响足够小，那么它们之和的分布仍然趋于高斯分布），
11、当n充分大时，有 （10）程序实现方法for(k=0;k=n-1;k+) t=0.0; for(i=1;i=12;i+) /取=1/12次 *r=(*r)*w+v; /求r,s=65536.0; w=2053.0;v=13849.0; *r=*r-m*s; /求整数后余下的数 t=t+(*r)/s; /累加到t中 ak=u+g*(t-6.0); /放到ak中，产生高斯白噪声过程 信噪比的改变 为了测试噪声信号的改变对LMS算法最终效果的影响，我组对噪声信号进行了信号能量强度处理。处理方法即使通过给已产生的信号乘上一个系数，通过系数值的改变来调节信号的能量。系数值得计算方法为 （11）式中，SN
12、R为信噪比，es为原始声音的均值，ns为噪声的均值。程序实现：snrate=sqrt(es/ns/pow(10,SNR/10); /开根号 pow 计算10的SNR/10次方的值，求能量,为了达到信噪比20DB而乘上的系数,SNR为信噪比3、LMS算法的实现 LMS最小均方差算法的原理已在前文给出。本实验就以该原理为基础进行设计。在此不多赘述。需要注意几点： LMS的收敛性LMS算法的收敛因子必须在如下范围内选取： （12）02/(MSmax) 其中Smax是抽头输入u(n)的功率谱密度的最大值，而滤波器长度M为中到大在满足收敛条件下，才有 （13）即，值得合理选取确保了系数向量的平均值接近于最优系数向量W。实际中，通常选得很小，选 （14）在实验中，选定N=10实验程序为：mu=mu+sn*sn;mu=1/(10*(mu); /sn为加噪信号，即为式（14）中x 抽头系数的自适应改变抽头系数的自适应改变是LMS算法中的一个重要部分，没有他的改变，自适应也就无从谈起。我们设计的滤波器权系数控制电路如图3 图3 FIR 滤波器中第个权系数的控制电路图根据权系数的计算方法 （15）我们得到程序设计：for(i=0;iBUFLEN;i+)
