LMS与RLS自适应滤波算法性能比较

应用技术研究作者简介：马文民男山东省济南圣泉集团股份有限公司工程师1LMS与RLS自适应滤波算法性能比较马文民【摘要】：介绍了自适应滤波器去除噪声的原理和从强噪声背景中采用自适应滤波提取有用信号的方法，并对最小均方(LMS,LeastMeanSquares)和递推最小二乘(RLS,RecursiveLeastSquares)两种基本自适应算法进行了算法原理、算法性能分析。计算机模拟仿真结果表明，这两种算法都能通过有效抑制各种干扰来提高强噪声背景中的信号。检测特性相比之下，RLS算法具有良好的收敛性能，除收敛速度快于LMS算法和NLMS算法以及稳定性强外，而且具有更高的起始收敛速率、更小的权噪声和更大的抑噪能力。【关键词】：自适应滤波；原理；算法；仿真引言：自适应滤波是近30年以来发展起来的一种最佳滤波方法。它是在维纳滤波，kalman滤波等线性滤波基础上发展起来的一种最佳滤波方法。由于它具有更强的适应性和更优的滤波性能。从而在工程实际中，尤其在信息处理技术中得到广泛的应用。自适应滤波的研究对象是具有不确定的系统或信息过程。“不确定”是指所研究的处理信息过程及其环境的数学模型不是完全确定的。其中包含一些未知因数和随机因数。任何一个实际的信息过程都具有不同程度的不确定性，这些不确定性有时表现在过程内部，有时表现在过程外部。从过程内部来讲，描述研究对象即信息动态过程的数学模型的结构和参数是我们事先不知道的。作为外部环境对信息过程的影响，可以等效地用扰动来表示，这些扰动通常是不可测的，它们可能是确定的，也可能是随机的。此外一些测量噪音也是以不同的途径影响信息过程。这些扰动和噪声的统计特性常常是未知的。面对这些客观存在的各种不确定性，如何综合处理信息过程，并使某一些指定的性能指标达到最优或近似最优，这就是自适应滤波所要解决的问题。在这几十年里，数字信号处理技术取得了飞速发展，特别是自适应信号处理技术以其计算简单、收敛速度快等许多优点而广泛被使用。它通过使内部参数的最优化来自动改变其特性。自适应滤波算法在统计信号处理的许多应用中都是非常重要的。在工程实际中，经常会遇到强噪声背景中的微弱信号检测问题。例如在超声波无损检测领域，因传输介质的不均匀等因素导致有用信号与高噪声信号迭加在一起。被埋藏在强背景噪声中的有用信号通常微弱而不稳定，而背景噪声往往又是非平稳的和随时间变化的，此时很难用传统方法来解决噪声背景中的信号提取问题。自适应噪声抵消技术是一种有效降噪的方法，当系统能提供良好的参考信号时，可获得很好的提取效果。与传统的平均迭加方法相比采用自适应平均处理方法还能降低样本数量。1自适应滤波器的基本原理所谓的自适应滤波，就是利用前一时刻以获得的滤波器参数的结果，自动的调节现时刻的滤波器参数，以适应信号和噪声未知的或随时间变化的统计特性，从而实现最优滤波。自适应滤波器实质上就是一种能调节其自身传输特性以达到最优的维纳滤波器。自适应滤波器不需要关于输入信号的先验知识，计算量小，特别适用于实时处理。由于无法预先知道信号和噪声的特性或者它们是随时间变化的，仅仅用FIR和IIR两种具有固定滤波系数的滤波器无法实现最优滤波。在这种情况下，必须设计自适应滤波器，以跟踪信号和噪声的变化。自适应滤波器的特性变化是由自适应算法通过调整滤波器系数来实现的。一般而言，自适应滤波器由两部分组成，一是滤波器结构，二是调整滤波器系数的自适应算法。自适应噪声抵消系统的核心是自适应滤波器，自适应算法对其参数进行控制，以实现最佳滤波。不同的自适应滤波器算法，具有不同的收敛速度、稳态失调和算法复杂度。根据自适应算法是否与滤波器输出有关，可将其分成开环算法和闭环算法两类。自适应噪声抵消器中利用了输出反馈，属于闭环算法。其优点是能在滤波器输入变化时保持最佳的输出，而且还能在某种程度上补偿滤波器元件参数的变化和误差以及运算误差。但其缺点是存在稳应用技术研究2定性问题以及收敛速度不高。所以探讨如何提高收敛速度、增强稳定性以满足信号处理的高效性、实时性，一直是人们研究的重点和热点。本文基于自适应噪声抵消对比研究了两类基本的自适应算法，并对它们在分离周期信号和随机噪声中呈现的滤波性能进行了分析。计算机仿真结果表明，RLS算法从背景噪声中提取有用信号的滤波性能明显优于LMS算法。2算法原理图1自适应滤波器原理框图图1给出了用自适应噪声抵消技术来解决噪声背景中的信号提取问题的基本原理。主输入端接收从信号源发来的信号s但是受到噪声源的干扰收到噪声vo。参考输入端的参考信号为vi是一个与有用信号s无关但与vo相关的噪声信号。主输入中含有待抵消的加性噪声，参考输入对准主输入中的噪声vo。利用两输入噪声的相关性和信号与噪声的独立性，使参考输入通过自适应滤波器与主输入中噪声分量逼近并相减，输出误差信号。自适应滤波算法决定滤波器对参考信号v1的处理，使得滤波器的输出尽可能地逼近主输入中的干扰成分。所以，在最佳准则意义下滤波器的输出v逼近vo等效于系统的输出e逼近s。从而在噪声对消器的输出端大大地提高了信噪比。但若参考通道除检测到噪声v1外，还收到信号分量，则自适应滤波器的输出中将包含信号分量，从而使噪声对消效果变坏。因此，为获得良好的噪声对消性能，应使参考通道检测到的信号尽可能小，在信号不可测的噪声环境拾取参考输入信号。3算法：根据自适应算法的优化准则的不同，自适应滤波算法可分为两类最基本的算法：最小均方(LMS)算法和递推最小二乘(RLS)算法。为了解决传统LMS算法存在梯度噪声放大问题，以及为克服常规的固定步长LMS自适应算法在收敛速率、跟踪速率与权失调噪声之间的要求上存在的较大矛盾，许多学者研究出了各种各样的改进型LMS算法，如归一化LMS算法和基于瞬变步长LMS自适应滤波算法以及基于离散小波变换的LMS自适应滤波算法等。a:LMS算法：自适应滤波器在时刻n的向量定义：抽头权向量：011()[(),(),...,()]TMWnbnbnbn参考输入向量：()[(),(1),...,(1)]TXnxnxnxnM()dn是主输入信号，()yn是期望输出值，()en是误差信号，也是系统输出值，M是滤波器长度。由维纳-霍夫方程可知，最小均方误差为：22*min([])[]jjEeEdWP实际上，该方程与维纳滤波器结果完全一样。自适应滤波器与维纳滤波器相比，其差别在于它增加了一个识别控制环节，将输出与期望值进行比较，利用误差()en去控制()Wn，使2[]jEe=最小值，从而得到()Wn的估计*()Wn。根据最优的数学算法最陡下降法，下一个权矢量1()jWn等于现在的权矢量()jWn加一个正比于梯度j的负值变化量，即有：1jjjWW通过梯度下降法：21[]|jjjjWWdEeWWdW推导可知：12jjjjWWeX其中jjjjedWX算法步骤：步骤一：初始化：步骤二：更新：1,2,3,...n滤波：()()()TynWnXn；误差估计：()()()endnyn；权向量更新：*(1)()2()()WnWnenXn；其中是用来控制稳定性和收敛速度的步长参数。为确保自适应过程的稳定性，必须满足应用技术研究302/inMP，其中2[()]inPEXn为输入功率。b:RLS算法：SISO系统动态过程的数学模型：11()()()()()AzzkBzuknk（1）其中()uk，()zk为输入输出量，()nk为噪声。式中11212()1...aannAzazazaz11212()...bbnnBzbzbzbz展开后得到：1212()(1)(2)...(1)(2)...()abnnbzkazkazkabukbukbukn模型（1）可化为最小二乘格式：()()()zkhknk（2）记1212[,,...,,,...,]abnnaaabbb为待估计的参数。()[(1),...,(),(1),...,()]abhkzkzknukukn，对于1,2,...kL（L为数据长度）。方程（2）构成一个线性方程组，写成()()()LLLzkHknk；(1)(2)()LzzZzL，(1)(2)()LhhHhL，(1)(2)()LnnnnL根据最小二乘法一次完成算法，其参数估计为：1ˆ()LLLLLSHHHZ。参数递推估计，每取得一次新的观测数据后，就在前次估计结果的基础上，利用新引入的观测数据对前次估计的结果，根据递推算法进行修正，减少估计误差，从而递推地得出新的参数估计值。这样，随着新观测数据的逐次引入，一次接一次地进行参数估计，直到参数估计值达到满意的精确程度为止。算法步骤：步骤一：初始化(0)0W；1(0)PI，其中I为单位矩阵；步骤二：更新1,2,...n计算更新增益矢量：()(1)()/[()(1)()]TgnPnXnXnPnXn；滤波：()(1)()TynWnXn；误差估计：()()()endnyn；更新权向量：()(1)()()WnWngnen；更新逆矩阵：1()[(1)()()(1)]TPnPngnXnPn；其中，()Pn为自相关矩阵()xxPn的逆矩阵，常数是遗忘因子，且01。总上所述：算法实现的主要步骤为：（1）数据采集与生成，取()dn，()Xn；（2）对参数的初始化；（3）自适应的滤波处理；（4）滤波器系数更新3计算机仿真结果与分析；为了检验两种自适应滤波算法在去噪应用中的滤波性能，下面对LSM算法和RLS算法进行计算机模拟仿真实验。其中采样频率为1000Hz，其算法用MATLAB语言实现。其中图2为幅度为2标准正弦波。图3为幅度为2正弦波叠加带限高斯白噪声的混迭信号，是系统的主输入信号。图4、图5分别为用LMS算法和RLS算法提取得到的正弦信号。表一各自适应滤波各参数设置从图上可以看出，用RLS自适应滤波算法提取得到的正弦信号效果较好。而LMS自适应滤波算法也能将信号提取出来，但是其滤波效果较差，存在没有滤除的随机噪声部分较多。4滤波器性能比较：由于LMS算法只是用以前各时刻的抽头参量等作该时刻数据块估计时的平方误差均方最小的准则，而未用现时刻的抽头参量等来对以往各时刻的数据块作重新估计后的累计平方误差最小的准则，所以LMS算法对非平稳信号的适应性差。RLS算法的基本思想是力图使在每个时刻对所有已输名称N（阶数）μλσLSM80．00026RLS80．990．1应用技术研究4入信号而言重估的平方误差的加权和最小，这使得RLS算法对非平稳信号的适应性要好。与LMS算法相比，RLS算法采用时间平均，因此，所得出的最优滤波器依赖于用于计算平均值的样本数，而LMS算法是基于集平均而设计的，因此稳定环境下LMS算法在不同计算条件下的结果是一致的。在性能方面，RLS的收敛速率比LMS要快得多，因此，RLS在收敛速率方面有很大优势。图6分别为RLS算法和LMS算法在处理过程中的误差曲线，它指出了在迭代过程中的误差减少过程。由图可见，RLS算法在迭代过程中产生的误差明显小于LMS算法。由此可见，RLS在提取信号时，收敛速度快，估计精度高而且稳定性好，可以明显抑制振动加速度收敛过程，故对非平稳信号的适应性强，而LMS算法收敛速度慢，估计精度低而且权系数估计值因瞬时梯度估计围绕精确值波动较大，权噪声大，不稳定。020406080100120140160180200-2-1.5-1-0.500.511.52幅度为2标准正弦波图2幅度为2标准正弦波020406080100120140160180200-4-3-2-101234幅度为2正弦波叠加带限高斯白噪声的混迭信号图3幅度为2正弦波叠加带限高斯白噪声的混迭信号020406080100120140160180200-2.5-2-1.5-1-0.