基于自适应滤波对音频信号的处理详解

1、word毕业设计论文 题目基于自适应滤波对音频信号的处理学生 完成日期 2012年5月 基于自适应滤波对音频信号的处理 摘要自适应滤波器是统计信号处理的一个重要组成局部。在信号的传播过程中，通常会受到噪声的干扰，这时候自适应滤波器就表达出了它的重要作用。传统的滤波器只能过滤一种或几种干扰噪声，然而随着信号传输环境的不断变化，传统的滤波器已经不能适应，这就需要一种根据环境的改变，通过自适应算法来改变滤波器的参数和结构，来到达更好滤波效果的滤波器。自适应滤波器是利用前一时刻获得的滤波参数，自动地调节、更新现时刻的滤波参数，以适应信号和噪声未知的统计特性，从而实现最优滤波。本文从自适应滤波器研究的意

2、义入手，介绍了自适应滤波器的根本理论思想，具体阐述了自适应滤波器的根本原理、算法及设计方法。自适应滤波器的算法是整个系统的核心。因LMS算法具有低计算复杂度、在平稳环境中的收敛性好、其均值无偏地收敛到维纳解和利用有限精度实现算法时的稳定性等特性，使LMS算法成为自适应算法中应用最广泛的算法，所以最终采用基于LMS算法设计自适应滤波器。对读取一段wav格式的音频文件采用 MATLAB 进行仿真，通过实验结果来表达该滤波器可以根据信号随时修改滤波参数，到达动态跟踪的效果，使滤波信号更接近于原始信号。关键词：自适应滤波；LMS算法； Matlab ABSTRACT The adaptive filt

3、er is an important part of the digital signal processing. in the spread of the signal process, usually subject to noise interference, and this was reflected when the adaptive filter out of its important role. The traditional filters can only filtration one or more of the noise interference, however,

4、 with the constant change of the signal transmission environment, the traditional filters can not adapt tochanges in the environment which requires a through adaptive algorithm tochange the filter parameters and structure of the device to reach the better the filtering effect of the filter. The adap

5、tive filter is to use a time before get the filter parameters have been automatically adjust and update the current moment of filter parameters, to adapt to the signal and noise statistical properties of the unknown in order to achieve optimum filter. This paper,from the adaptive filter the signific

6、ance of research and introduced its elementary theory, algorithm and design method. The core of the whole system is the auto-adapted filter's algorithm. For LMS algorithm has low computational complexity, in the environment of steady convergence, the mean unbiased to converge to a wiener solutio

7、n and the use of the limited precision of the stability of the algorithm and other characteristics . LMS algorithm as adaptive algorithm in the application of the most a wide range of algorithms.，So Finally the design of adaptive filters based on LMS algorithm.To read a wav format audio files using

8、MATLAB simulation and experimental results to reflectthe filter can be modified at any time according to the signal of the filter parameters to achieve the effect of dynamic tracking, so that the filtered signalcloser to the original signal. Keywords: adaptive filter ; LMS algorithm ; Matlab 目 录第一章

9、绪论6 1.1 引言6 1.2 自适应滤波理论的开展6 1.3 自适应滤波器的根本概念7 1.3.1 自适应滤波器特征及特点7 1.4 自适应滤波器原理8 1.5 课题研究意义和目的9第二章 自适应LMS算法的研究11 2.1自适应LMS算法的开展11 2.1.1 LMS算法的历史11 2.1.2 LMS算法的现状11 2.1.3 LMS算法的开展前景11 2.2 自适应LMS算法的原理12 2.3 自适应LMS算法降噪的原理14第三章 WAV文件格式16 3.1 WAV文件格式介绍16 3.2 PCM数据的存放方式17 3.3 WAV文件格式17第四章 仿真实验20 4.1 MATLAB语言

10、介绍20 4.1.1 MATLAB的特点20 4.1.2 MATLAB的优势20 4.2 MATLAB仿真21 4.2.1 MATLAB程序仿真21 4.2.2 仿真结果22 4.3 LMS自适应算法性能分析24第五章 总结25参考文献26致 谢27 第一章 绪论 1.1 引言 在过去的40年中，信号处理技术获得了巨大开展，数字电路设计技术的进步，是促使人们对数字信号处理领域越来越感兴趣的关键技术原因。相应的数字信号处理系统之所以具有吸引力，是因为它具有可靠性好、精度高、体积小以及灵活性强等特点。 滤波技术是信号处理中的一种根本方法和技术，尤其数字滤波技术使用广泛，数字滤波理论的研究及其产品的

11、开发一直受到很多国家的重视。从总的来说滤波可分为经典滤波和现代滤波。经典滤波要求信号和噪声的统计特性，如维纳滤波和卡尔曼滤波。现代滤波那么不要求信号和噪声的统计特性，如自适应滤波。自适应滤波的原理就是利用前一时刻获得的滤波参数等结果，自动地调节现时刻的滤波参数，从而到达最优化滤波。自适应滤波具有很强的自学习，自跟踪能力，适用于平稳和非平稳随机信号的检测和估计。自适应滤波一般包括3个模块：滤波结构、性能判断和自适应算法。其中，自适应滤波算法一直是人们的研究热点，包括线性自适应算法和非线性自适应算法，非线性自适应算法具有更强的信号处理能力，但计算比拟复杂，实际应用最多的仍然是线性自适应算法。线性自

12、适应算法的种类很多，有LMS自适应滤波算法、R路自适应滤波算法、变换域自适应滤波算法、仿射投影算法、共轭梯度算法等。1.2 自适应滤波理论的开展 早在20世纪40年代，就对平稳随机信号建立了维纳滤波理论。根据有用信号和干扰噪声的统计特性自相关函数或功率谱，以线性最小均方误差估计准那么所设计的最正确滤波器，称为维纳滤波器。这种滤波器能最大程度地滤除干扰噪声，提取有用信号。但是，当输入信号的统计特性偏离设计条件，那么它就不再是最正确的了，这在实际应用中受到了限制。到60年代初，由于空间技术的开展，出现了卡尔曼滤波理论，即利用状态变量模型对非平稳、多输入多输出随机序列作最优估计。现在，卡尔曼滤波器已

13、成功地应用到许多领域，它既可对平稳的和非平稳的随机信号作线性最正确滤波，也可作非线性滤波。实质上，维纳滤波器是卡尔曼滤波器的一个特例。假设设计卡尔曼滤波器时，必须知道产生输入过程的系统的状态方程和测量方程，即要求对信号和噪声的统计特性有先验知识。但在实际中，往往难以预知这些统计特性，因此实现不了真正的最正确滤波。 Widrow B.等于1967 年提出的自适应滤波理论，可使自适应滤波系统的参数自动地调整而到达最正确状况，而且在设计时，只需要很少的或是根本不需要任何关于信号与噪声的先验统计知识。这种滤波器的实现差不多像维纳滤波器那样简单，而滤波性能几乎如卡尔曼滤波器一样好。因此，近十年来，自适应

14、滤波理论的方法得到了迅速开展。1.3 自适应滤波器的根本概念 但凡有能力进行信号处理的装置都可以称为滤波器。在近代电信装备和各类控制系统中，滤波器应用极为广泛。在所有的电子部件中，使用最多，技术最复杂要算滤波器了。滤波器的优劣直接决定产品的优劣，所以，对滤波器的研究和生产历来为各国所重视。滤波器是一种用来消除干扰杂讯的器件，将输入或输出经过过滤而得到纯洁的交流电。您可以通过根本的滤波器积木块二阶通用滤波器传递函数，推导出最通用的滤波器类型：低通、带通、高通、陷波和椭圆型滤波器。传递函数的参数f0、d、hHP、hBP 和hLP，可用来构造所有类型的滤波器。转降频率f0为s项开始占支配作用时的频率

15、。设计者将低于此值的频率看作是低频，而将高于此值的频率看作是高频，并将在此值附近的频率看作是带内频率。阻尼d用于测量滤波器如何从低频率转变至高频率，它是滤波器趋向振荡的一个指标，实际阻尼值从0至2变化。高通系数hHP是对那些高于转降频率的频率起支配作用的分子的系数。带通系数hBP是对那些在转降频率附近的频率起支配作用的分子的系数。低通系数hLP是对那些低于转降频率的频率起支配作用的分子的系数。设计者只需这5个参数即可定义一个滤波器。1.3.1 自适应滤波器特征及特点（1） 滤波器是线性时不变的。（2） 设计过程用到希望的带通、转换波段、带通波纹和阻带衰减。（3） 因为滤波器是频率选择性的，所以

16、当输入信号的各个局部占据不重叠频带时，滤波器工作得最好。例如，它可以轻易别离频谱不重叠的信号和附加噪声。（4） 滤波系数在设计阶段选定，并在滤波器的正常运行中保持不变。然而，在实际应用中有很多问题不能用固定数字滤波器很好地解决，因为我们没有充足的信息去设计固定系数的数字滤波器，或设计规那么会在滤波器正常运行时改变。绝大数这些应用都可以用特殊的智能滤波器，即常说的自适应滤波器来成功解决。自适应滤波器的显著特征是：它在工作过程中不需要用户的干预就能改变响应以改善性能。（5） 滤波结构。这个模块使用输入信号的测量值产生滤波器的输出。如果输出是输入测量值的线性组合，那么这个滤波器就是线性的，否那么称为

17、非线性的。结构有设计者设定，它的参数由自适应算法调整。（6） 性能标准。自适应滤波器的输出和期望的响应当可获得时由 COP 模块处理，并参照特定应用的需要来评估它的质量。（7） 自适应算法。自适应算法使用性能标准的数值或它的函数、输入的测量值和期望值的响应来决定如何修改滤波器的参数，以改善性能。 1.4 自适应滤波器原理 所谓的自适应滤波，就是利用前一刻已获得的滤波参数的结果，自动的调节现时刻的滤波器参数，以适应信号和噪声未知的或随时间变化的统计特性，从而实现最优滤波。自适应滤波器实质上就是一种能调节其自身传输特性以到达最优的维纳滤波器。自适应滤波器不需要关于输入信号的先验知识，计算量小，特别

18、适用于实时处理。图1-1 自适应滤波器原理图图1-1 描述的是一个通用的自适应滤波估计问题，图中离散时间线性系统表示一个可编程滤波器，它的冲击响应为h(n)，或称其为滤波参数。自适应滤波器输出信号为y(n)，所期望的响应信号为d(n)，误差信号e(n)为d(n) 与y(n)之差。这里，期望响应信号d(n) 是根据不同用途来选择的，自适应滤波器的输出信号y(n)是对期望响应信号d(n)进行估计的，滤波参数受误差信号e(n)的控制并自动调整，使y(n)的估计值等于所期望的响应d(n).因此，自适应滤波器与普通滤波器不同，它的冲击响应或滤波参数是随外部环境的变化而变化的，经过一段自动调整的收敛时间到

19、达最正确滤波的要求。但是，自适应滤波器本身有一个重要的自适应算法，这个算法可以根据输入、输出及原参数量值，按照一定准那么改变滤波参量，以使它本身能有效地跟踪外部环境的变化。通常，自适应滤波器是线性的，因而也是一种线性移变滤波器。当然，它可推广到自适应非线性滤波器。 在图1-1中，离散时间线性系统可以分为两类根本结构，其中一类为非递归型横向结构的数字滤波器，它具有有限的记忆，因而称之为有限冲激响应FIR系统，即自适应FIR滤波器。另一类为递归型数字滤波器结构，理论上，它具有无限的记忆，因而称之为无限冲激响应IIR系统，即自适应IIR滤波器。对于上述两类自适应滤波器，还可以根据不同的滤波理论和算法

20、，分为结构不同的自适应滤波器，它们的滤波器性能也不完全相同。 1.5 课题研究意义和目的对自适应滤波算法的研究是当今自适应信号处理中最为活泼的研究课题之一。 Windrow 等于1967年提出的自适应滤波系统的参数能自动的调整而到达最优状况，而且在设计时，只需要很少的或根本不需要任何关于信号与噪声的先验统计知识。这种滤波器的实现差不多像维纳滤波器那样简单，而滤波器性能几乎如卡尔曼滤波器一样好。自适应滤波器与普通滤波器不同，它的冲激响应或滤波参数是随外部环境的变化而变化的，经过一段自动调节的收敛时间到达最正确滤波的要求。自适应滤波器本身有一个重要的自适应算法，这个算法可以根据输入、输出及原参量信

21、号按照一定准那么修改滤波参量，以使它本身能有效的跟踪外部环境的变化。因此，自适应数字系统具有很强的自学习、自跟踪能力和算法的简单易实现性。 自适应滤波技术的核心问题是自适应算法的性能问题，提出的自适应算法主要有最小均方(LMS)算法、递归最小二乘(RLS)算法及相应的改良算法如：归一化(NLMS)算法、变步长(SVSLMS)算法、递归最小二乘方格形(RLSL)算法等。这些算法各有特点，适用于不同的场合。研究自适应算法是自适应滤波器的一个关键内容。最小均方误差(LMS，The least Mean square)算法是线性自适应滤波算法中最根本的两类算法之一，其主要思想 是基于最小均方误差准那么

22、，使滤波器的输出信号与期望输出信号之间的均方误差最小。由于LMS算法简 单有效、鲁棒性好、易于实现，得到了广泛的应用。目前应用最多的是系统辨识、回波消除、自适应谱线增强、自适应信道均衡、语音线性预测、自适应天线阵等诸多领域。 第二章 自适应LMS算法的研究2.1自适应LMS算法的开展2.1.1 LMS算法的历史 1955-1966年期间美国通用公司在研制天线的过程中，为抑制旁瓣，由windows和hoff在60年代初提出了根本LMS算法。随后又开展出了归一化算法和加遗忘因子LMS算法。1977年，makjoul提出了格型滤波器，并由此开展出LMS自适应格型滤波器算法。Herzberg、cohe

23、n和beery提出了延时LMSDLMS算法。2002年，尚勇，吴顺君，项海格提出了并行延时LMS算法。此外，还有复数LMS算法、数据块LMS算法等。2.1.2 LMS算法的现状因LMS算法具有低计算复杂度、在平稳环境中的收敛性好、其均值无偏地收敛到wiener解和利用有限精度实现算法时的稳定性等特性，使LMS算法成为自适应算法中应用最广泛的算法。由于LMS算法的广泛应用，以及在实际条件下，为解决实际问题，基于LMS算法的新LMS类算法不断出现。2.1.3 LMS算法的开展前景 因LMS算法是自适应滤波器中应用最广泛的算法，所以可以说，自适应滤波的开展前景也就是LMS算法的开展前景。它主要包括以

24、下几个方面的应用：1、系统辨识和建模(System Identification and Modeling)。自适应滤波器作为估计未知系统特性的模型。2、自适应信道均衡(Adaptive Channel Equlization)。在数字通信中采用自适应信道均衡器，可以减小传输失真，以及尽可能地利用信道带宽。3、回波消除(Echo Cancellation)。在2线和4线环路 系统中，线路间存在杂散电路耦合，这些杂散导致阻抗不匹配，从而形成了信号的反射，也就是我们在线路两端听到的回声。这种回波能对高速数据传输造成灾难性的后果。回波消除就是预先估计一个回波，然后用返回信号来减此回波，从而到达回波消

25、除的目的。消除心电图中的电源干扰就是它的一个具体应用。4、线性预测编码Linear Predictive Coding。近年来，对语音波形进行编码，它可以大大降低数据传输率。在接收端使用LPC分析得到的参数，通过话音合成器重构话音。合成器实际上是一个离散的随时间变化的时变线性滤波器。时变线性滤波器既当作预测器使用，又当作合成器使用。分析语音波形时作预测器使用，合成语音时作话音生成模型使用。5、自适应波束形成Adaptive Beaamforming。频谱资源越来越紧张，利用现有频谱资源进一步扩展容量成为通信开展的一个重要问题。智能天线技术利用阵列天线替代常规天线，它能够降低系统干扰，提高系统容

26、量和频谱效率，因此智能天线技术受到广泛关注。自适应束波形成通过调节天线各阵元的加权幅度和相位，来改变阵列的方向图，使阵列天线的主瓣对准期望用户，从而提高接收信噪比，满足某一准那么下的最正确接收。在雷达与声纳的波束形成中，自适应滤波器用于波束方向控制，并可在方向图中提供一个零点以便消除不希望的干扰。其应用还有噪声中信号的滤波、跟踪、谱线增强以及预测等。2.2 自适应LMS算法的原理 LMS算法的根本思想：调整滤波器自身参数，使滤波器的输出信号与期望输出信号之间的均方误差最小，这样系统输出为有用信号的最正确估计。LMS算法的结构图如图 2-1 所示。输入信号x依次经过 m-1 个延时单位，在 n

27、时刻构成了一个信号矢量 x(n)。 图2-1 LMS算法结构图输入信号矢量为：加权矢量即滤波器抽头系数矢量为：相应的滤波器的输出为： 1yn相对于期望信号 dn的误差为： 2根据最小均方误差MSE准那么，最正确的滤波器参数 应使得性能函数均方误差 为最小，根据 xn及其过去值 xn-1，xn-2，xn-3，.,xn-m+1去估计 dn，wn的最优估计值 使估计误差的均方值为最小， 满足正那么方程： 3式中为输入信号矢量 xn的自相关函数， 为输入信号矢量 xn与实际需要信号 dn的互相关函数。直接求解该正那么方程的运算量很大，尤其是加权矢量参数个数很多时。LMS算法采用递推收敛到最正确值 的方

28、法来寻求这个最正确的滤波器加权矢量参数 。 LMS算法推导如下：利用最优方法中的最速下降法求最正确权系数向量的近似值。最速下降法，即“下一时刻权系数向量 wn+1应该等于“现时刻权系数向量 wn加上一个负均方误差梯度 的比例项，即： 4为控制收敛速度与稳定性的常数，称之为收敛因子。按照近似方法，直接取作为均方误差 的估计值，即： 5 6 于是可得： 7 其中，为步长，所以为方便起见，仿真时递推公式采用下式： 8 收敛条件为： ，是输入信号自相关矩阵的最大特征值。算法步骤：（1） 算法初始化，W0=0，或者根据对环境的先验知识确定初始权值。（2） 对每一时刻 n=0,1，2，.计算滤波：;误差估

29、计：；权向量更新：其中是用来控制稳定性和收敛速率的步长参数。为确保自适应过程的稳定性，必须满足，其中 为输入功率。综上所述，LMS算法主要步骤为：1数据采集与生成，取得dn和Xn。2自适应滤波算法的参数初始化。3自适应滤波处理。4滤波器系统更新。2.3 自适应LMS算法降噪的原理自适应噪声抵消系统的核心是自适应滤波器，自适应算法对其参数进行控制，以实现最正确滤波，因为这种方法比其他方法多用了一个参考噪声作为辅助输入，从而获得了比拟全面的关于噪声的信息，因而能得到更好的降噪效果，特别是在辅助输入噪声与语音中的噪声完全相关的情况下，自适应噪声抵消法能完全排除噪声的随机性，彻底地抵消语音中的噪声成分

30、，从而无论在信噪比SNR方面还是在语音可懂度方面都能获得较大的提高，自适应噪声抵消器的原理如图 2-2 所示： 图 2-2下面将采用最小均方误差准那么来分析自适应噪声抵消的原理。原始输入为受扰信号；参考输入为，与相关，而与信号sn不相关。原始输入加到自适应滤波器的a 端，参考输入加到自适应滤波器的b 端。图中自适应滤波器AF 接收误差 en的控制，调整权矢量wn，使得它的输出 yn趋于等于 dn中与它相关的 ，于是en作为dn与 yn之差就接近等于信号 sn。证明如下：由图3得： 1 2 两边取数学期望，可得： 3 自适应过程就是自动调节权重 使均方最小的过程，式3中第一项为信号功率，与wn无

31、关。由于信噪不相关，第三项为0。所以，要使 成立，即式3中第二项应最小，即： 4 5 由图 1 得： 6 7 当 被最小化， 也被最小化，en以最小均方差趋于sn，可能的最好情况为：，那么en=sn；最小化输出功率使输出完全没有噪声，是自适应滤波器最理想的情况。因此，自适应滤波器可以用来从噪声中提取信号。 第三章 WAV文件格式3.1 WAV文件格式介绍 WAV文件格式是一种由微软和IBM联合开发的用于音频数字存储的标准，它采用RIFF文件格式结构，非常接近于AIFF和IFF格式。符合RIFF Resource Interchange File Format标准。所有的WAV都有一个文件头，这

32、个文件头有音频流的编码参数。 WAV文件作为多媒体中使用的声波文件格式之一，它是以RIFF格式为标准的。每个WAV文件的头四个字节便是“RIFF。WAV文件由文件头和数据体两大局部组成。其中文件头又分为RIFFWAV文件标识段和声音数据格式说明段两局部。下列图3-1为WAV文件格式说明表。常见的声音文件主要有两种，分别对应于单声道11.025KHz采样率、8Bit的采样值和双声道44.1KHz采样率、16Bit的采样值。采样率是指：声音信号在“模数转换过程中单位时间内采样的次数。采样值是指每一次采样周期内声音模拟信号的积分值。对于单声道声音文件，采样数据为八位的短整数short int 00H

33、-FFH； 而对于双声道立体声声音文件，每次采样数据为一个16位的整数int，高八位和低八位分别代表左右两个声道。WAV文件数据块包含以脉冲编码调制PCM格式表示的样本。WAV文件是由样本组织而成的。在单声道WAV文件中，声道0代表左声道，声道1代表右声道。在多声道WAV文件中，样本是交替出现的。文件头偏移地址字节数数据类型内 容00H4char"RIFF"标志04H4longint 文件长度08H4char"WAV"标志0CH4char"fmt"标志10H4 过渡字节不定14H2int格式类别10H为PCM形式的声音数据)

34、16H2int单声道为1，双声道为2通道数18H2int采样率每秒样本数，表示每个通道的播放速度1CH4long波形音频数据传送速率，其值为通道数×每秒数据位数×每样 本的数据位数8。播放软件利用此值可以估计缓冲区的大小22H2 每样本的数据位数，表示每个声道中各个样本的数据位数。如果有多 个声道，对每个声道而言，样本大小都一样。 24H 4 char 数据标记符data 28H 4 long int 语音数据的长度 图 3-13.2 PCM数据的存放方式PCM数据的存放方式： 8位单声道 0声道 0声道 8位立体声 0声道左 1声道右 0声道左 1声道右 16位

35、单声道 0声道低字节 0声道高字节 0声道低字节 0声道高字节 16位立体声 0声道左低字节 0声道左高字节 1声道右低字节 1声道 右高字节 WAV文件的每个样本值包含在一个整数i中，i的长度为容纳指定样本长度所需 的最小字节数。首先存储低有效字节，表示样本幅度的位放在i的高有效位上， 剩下的位置为0，这样8位和16位的PCM波形样本的数据格式如下所示。样本大小数据格式最大值最小值8位PCMunsigned int225016位PCMint32767 图 3-23.4 WAV文件格式WAVE文件是非常简单的一种RIFF文件，它的格式类型为"WAVE"。RIFF块包含两个子

36、块，这两个子块的ID分别是"fmt"和"data",其中"fmt"子块由结构PCMWAVEFORMAT所组成，其子块的大小就是sizeofof(PCMWAVEFORMAT),数据组成就是PCMWAVEFORMAT结构中的数据。标志符RIFF数据大小格式类型"WAVE""fmt"Sizeof(PCMWAVEFORMAT)PCMWAVEFORMAT"data"声音数据大小声音数据        

37、60;                                       图 3-3 WAVE文件结构 "data"子块包含WAVE文件的数字化波形声音数据，其存放格式依

38、赖于"fmt"子块中wFormatTag成员指定的格式种类，在多声道WAVE文件中，样本是交替出现的。如16bit的单声道WAVE文件和双声道WAVE文件的数据采样格式分别如图四所示：16位单声道：采样一采样二低字节高字节低字节高字节16位双声道：采样一左声道右声道低字节高字节低字节高字节            图 3-4 WAVE文件数据采样格式3.5 WAV文件信息的具体应用 WAV文件中包括了对原始声音的高速率采样，并且以WAVE_PCM_FORMA

39、T脉派编码调变格式，我们可以在VISUALC+程序中实现，在读出WAVEHDR文件头之后，下面就是原始声音的高速率采样信息，我们可以对它作多方面的信息处理。 波形显示我们可以以时域-幅度的方式显示出原始声音的波形，这是最简单同时也是最直接的信息处理方式。在时域范围内，我们可以观察该信号波形是否连续，中间是否有跳变等。 3.5.2 频谱显示 我们可以以频域-幅度的方式显示出原始声音的频谱，在对原始信号经过FFT变换之后，可以得到该信号的频谱，进而得到该信号的能量集中带，分布特征，谱对称系数等等。 3.5.3 用于语音信号识别 讲话者的个体识别是语音信号处理的一个重要内容，但它的一个前提条件是必须

40、提供语音信号的数字波形，通常的方法是将原始的语音信号进行放大、抗混叠滤波、A/D采样、数值编码，最终得到语音信号的数字波形，通常多采用硬件处理，费时费力，如果我们借助非常成熟的声卡技术，将WAV文件翻开，就非常方便地得到语音信号的数字波形，为下一步进行语音信号识别提供良好的前端预处理。 第四章 仿真实验4.1 MATLAB语言介绍当计算涉及矩阵运算或画图时，利用FOTRRAN和C语言等计算机语言进行程序设计是一项很麻烦得工作，不仅需要对所利用的有关算法有深刻的了解，还需要熟练掌握所用语言的语法和编程技巧，MATLAB正是为了免除上述局面产生的，它是集命令翻译，科学计算于一身的一套交互软件系统，

41、在MATLAB下，矩阵的运算变得异常的容易，后来又增添了丰富多彩的图形处理及多媒体功能，使得MATLAB的应用范围越来越广泛。MATLAB具有强大的语言功能，它允许用户以数学形式的语言编写程序，比BASIC、FOTRRAN和C等语言更接近于书写计算公式的思维方式。它的操作和功能函数指令是以平时计算机和数学书上的简单英文单词表达的。由于它在很长一段时间内是用C语言开发的，它的不多的几个程序流控制语句同语言差异甚微，初学者很容易掌握。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创立用户界面、连接其他编程语言的程序等。主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融

42、建模设计与分析等领域。4.1.1MATLAB的特点（1） 其高级语言可用于技术计算。（2） 开发环境可对代码、文件和数据进行管理。（3） 交互式工具可以按迭代的方式探查、设计及求解问题。（4） 数学函数可用于线性代数、统计、傅里叶分析、筛选、优化以及数值积分等。（5） 二维和三维图形函数可用于可视化数据。（6） 各种工具可用于构建自定义的图形用户界面。（7） 各种函数可将基于MATLAB的算法与外部应用程序和语言集成。（8） 不支持大写输入，内核仅仅支持小写。4.1.2 MATLAB的优势（1） 友好的工作平台和编程环境。（2） 简单易用的程序语言。（3） 强大的科学计算机数据处理能力。（4）

43、 出色的图形处理功能。（5） 应用广泛的模块集合工具箱。（6） 实用的程序接口和发布平台。（7） 应用软件开发包括用户界面。4.2 MATLAB仿真 MATLAB程序仿真使用MATLAB编程，采用自适应滤波器技术实现语音去噪过程，程序如下：clear all;clc;hord=20; % - signal and noise generation - s=wavread('e.wav'); %读取wav文件Sound(s); %声音的回放len=length(s);n=1:len; v=(sin(0.6*n)+2*sin(2*pi*n)+sin(0.5*n); %噪声x=(si

44、n(0.6*n)+sin(2*pi*n+0.8)+sin(0.5*n+2*pi); %参考输入d=s+10*v; wavrite(x,15000,16,'d:daizao'); ws=zeros(hord+1,len-hord); errs=zeros(len,1); mus=0.008; %迭代步长 % - LMS noise cancellation - for i=hord+1:len-hord j=i-hord:1:i; x1=x(j); ys(i)=ws(:,i)'*x1; errs(i)=d(i)-ys(i); q(i)=errs(i)-s(i); n1=hord+1;n2=len-hord; Eq(i)=sum(abs(q(i).2); Pq(i)=Eq(i)/(n2-n1); ws(:,i+1)=ws(:,i)+mus*x1*errs(i); end % WAVWRITE(Y,FS,WAVEFILE) wavwrite(errs,15000,'d:shuchu'); figure(1) plot(s),title('original'); figure(2) plot(d),title('no