一种新的变步长LMS算法及在DSP上实现

liyf

                    　　0 引言
　　最小均方误差(LMS)算法由Widrow和Hoff等人于1960年提出，由于其结构简单，计算量小，稳定性好，易于实现等优点而得到广泛的应用。LMS算法的缺点是收敛速度慢，它克服不了收敛速度和稳态误差这一对固有矛盾：在收敛的前提下，如果步长取较大值，虽然收敛速度能得到提高，但稳态误差会随之增大，反之稳态误差虽然降低但收敛速度就会变慢。为解决这一矛盾，人们提出了许多改进型自适应算法。其中很大一类是变步长LMS算法。文献[4]提出Sigmoid函数变步长LMS算法(SVSLMS)。该算法在初始阶段或未知系统的系数参数发生变化时，其步长较大，从而使该算法有较快的收敛速度；而在算法收敛后，不管主输入端干扰信号e(n)有多大，都保持很小的调整步长，从而获得较小的稳态失调噪声。但Sigmoid函数过于复杂，且在误差e(n)接近零处变化太大，不具有缓慢变化的特性，使得SVSLMS算法在自适应稳态阶段仍有较大的步长变化；文献[5]提出的算法引入了多个调整参数，因而步长因子不易设计和控制；文献[6-8]提出了3种与误差信号成非线性关系的步长设计方法，该类算法具有较好的收敛性能，但3种算法在计算步长因子时，都存在指数运算。在数字信号处理中，进行一次指数运算需要的计算量，相当于进行多次乘法运算的计算量。
　　因此这类算法在实现时，增大了计算复杂度。为克服上述变步长LMS自适应滤波器存在的不足，在此提出了一种新的变步长LMS自适应滤波算法，该算法具有良好的收敛性能，较快的收敛速度，较小的稳态误差．良好的鲁棒性，并且在求变步长因子时计算量较小。
　　1 新的变步长LMS算法分析
　　基本的固定步长LMS算法的迭代公式可以表述为：





　　2 仿真及结果分析
　　下面通过计算机仿真来验证算法的收敛性能。仿真条件为：自适应滤波器的阶数为L=2；未知系统的FIR系数为W=[0，0]T；参考输入信号x(n)是零均值，方差为1的高斯白噪声；v(n)为与x(n)不相关的高斯白噪声。分别做200次独立的仿真，采样点数为1 000，然后求其统计平均，得出学习曲线。
　　图3是α固定，不同β值对应的收敛曲线。随着β值的增大，算法的收敛速度逐渐加快。图4是β保持不变，不同α值对应的收敛曲线，随着α逐渐减小，算法的误差也随之减小，但达到稳态的时间逐渐增加。





　　程序在CCS环境下编译，链接生成公共目标代码文件，再load到DSP中运行。本文采用C语言编程。图6为被噪声污染的输入信号，图7为滤波后的信号。比较两图可以看出，混合信号经过本文算法处理后噪声得到了很好的抑制，而原始信号得到了明显的加强和改善。
　　4 结语
　　本文提出了一种新的变步长LMS算法。该算法通过对最佳Wiener解的分析，并用反正切函数建立了步长因子μ(n)与误差e(n)的非线性函数关系。本算法具有初始阶段和未知系统时变阶段步长自动增大而稳态时步长很小的特点，且具有良好的鲁棒性和较小的计算量。通过在DSP上的实现，说明此算法具有可实行，并且实现方法简单，滤波效果好。（中电网 吴波，李庆梅，刘树东）