(LMS算法实验报告.docVIP

基于LMS算法的自适应均衡器设计及MATLAB实现 一.实验目的 了解LMS算法的基本原理； 掌握MATLAB的使用方法； 初步体会分析问题、研究问题的基本步骤和方法，为以后科研积累经验。 二.实验原理 1. LMS算法简介 在移动通信环境中，多径传播效应和频率选择性衰落会导致传输信号失真。失真主要表现为码间干扰,码间干扰是降低数字通信系统性能的一个主要因素。在这样的信道条件下设计实际的数字通信系统以高速传输数据时，往往不能获得足够准确的信道频率响应用于调制和解调器的最佳滤波器的设计。这是因为在每次通信时信道的路由不同,对于这样的信道，要设计最佳固定解调滤波器是不可能的。在这样的情况下，应该采取信道均衡的方式以减小失真。 信道均衡是通信技术和信号处理的基本问题之一，其目的在于克服传送的符号码和符号码之间的相互干扰，这种干扰是因为信道的非理想特性造成的。由于通信信道可能是未知和变化的，就需要自适应的调整均衡器，使得整个传输系统输出的符号码和符号码之间的干扰被消除。信道均衡可以利用发送的训练信号来开始，这称为自动均衡。 在设计自适应均衡器的多种方法中，最小均方自适应算法（LMS）采用梯度搜索法，这使收敛到最优解远比其他算法快，而且该算法原理简单，实施容易，所以目前这一算法已广泛用于计算自适应滤波器的权系数。 LMS算法的原理 （1）自适应滤波原理??? 自适应滤波器的特性变化是由自适应算法通过调整滤波器系数来实现的。一般而言，自适应滤波器由两部分组成，一是滤波器结构，二是调整滤波器系数的自适应算法。自适应滤波器的结构采用FIR或IIR结构均可，由于IIR滤波器存在稳定性问题，因此一般采用FIR滤波器作为自适应滤波器的结构。图1给出了自适应滤波器的一般结构。??? 图1为自适应滤波器结构的一般形式，图中x(n)为输入信号，通过参数可调的数字滤波器后产生输出信号y(n)，将输出信号y(n)与标准信号(或者为期望信号)d(n)进行比较，得到误差信号e(n)。e(n)和x(n)通过自适应算法对滤波器的参数进行调整，调整的目的使得误差信号e(n)最小。 ??? 自适应滤波器设计中最常用的是FIR横向型结构。图2是横向型滤波器的结构示意图。 其中：x(n)为自适应滤波器的输入；w(n)为自适应滤波器的冲激响应：w(n)={w(O)，w(1)，…，w(N-1)}；y(n)为自适应 ??? 自适应滤波算法??? 自适应滤波器除了包括一个按照某种结构设计的滤波器，还有一套自适应的算法。自适应算法是根据某种判断来设计的。自适应滤波器的算法主要是以各种判据条件作为推算基础的。通常有两种判据条件：最小均方误差判据和最小二乘法判据。LMS算法是以最小均方误差为判据的最典型的算法，也是应用最广泛的一种算法。??? 最小均方误差(Least Mean Square，LMS)算法是一种易于实现、性能稳健、应用广泛的算法。所有的滤波器系数调整算法都是设法使y(n)接近d(n)，所不同的只是对于这种接近的评价标准不同。LMS算法的目标是通过调整系数，使输出误差序列e(n)=d(n)-y(n)的均方值最小化，并且根据这个判据来修改权系数，该算法因此而得名。误差序列的均方值又叫“均方误差”(Mean Sqluare Error，MSE)。??? 理想信号d(n)与滤波器输出y(n)之差e(n)的期望值最小，并且根据这个判据来修改权系数wi(n)。由此产生的算法称为LMS。均方误差ε表示为： ??? 对于横向结构的滤波器，代入y(n)的表达式： ??? 其中：R=E[X(n)XT(n)]为N×N的自相关矩阵，它是输入信号采样值间的相关性矩阵。P=E[d(n)X(n)]为N×1互相关矢量，代表理想信号d(n)与输入矢量的相关性。在均方误差ε达到最小时，得到最佳权系数它应满足下式： ??? ??? 这是一个线形方程组，如果R矩阵为满秩的，R-1存在，可得到权系数的最佳值满足：W*=R-1p。用完整的矩阵表示为： ??? ??? 显然φx(m)=E[x(n)x(n-m)]为x(n)的自相关值，φxd(R)=E[x(n)d(n一k)]为x(n)与d(n)互相关值。在有些应用中，把输入信号的采样值分成相同的一段(每段称为一帧)，再求出R，P的估计值得到每帧的最佳权系数。这种方法称为块对块自适应算法。如语音信号的线性预测编码LPC就是把语音信号分成帧进行处理的。R，P的计算，要求出期望值E，在现实运算中不容易实现，为此可通过下式进行估计： ??? 用以上方法获得最佳W*的运算量很大，对于一些在线或实时应用的场合，无法满足其时间要求。大多数场合使用迭代算法，对每次采样值就求出较佳权系数，称为采样值对采样值迭代算法。迭代算法可以避免复杂的R-1和P的运