基于RLS算法自适应滤波器

基于RLS算法自适应滤波器的设计摘要自适应滤波器是统计信号处理的一个重要组成部分。在实际应用中，由于没有充足的信息来设计固定系数的数字滤波器，或者设计规则会在滤波器正常运行时改变，因此需要研究自适应滤波器。凡是需要处理未知统计环境下运算结果所产生的信号或需要处理非平稳信号时，自适应滤波器可以提供非自适应方法所不可能提供的新的信号处理能力。而且其性能通常远优于用常方法设计的固定滤波器。本文从自适应滤波器研究的意义入手，介绍了自适应滤波器的基本理论思想，具体阐述了自适应滤波器的基本原理、算法及设计方法。自适应滤波器的算法是整个系统的核心。对RLS算法自适应滤波器做了详细的介绍，采用改进的RLS算法设计自适应滤波器，并采用MATLAB进行仿真，通过实验结果来体现该滤波器可以根据信号随时修改滤波参数，达到动态跟踪的效果，使滤波信号更接近于原始信号。关键词：自适应滤波器，RLS算法，噪声消除，FIR张少博：基于RLS算法自适应滤波器的设计-1-第1章绪论1.1课题研究意义和目的滤波技术是信号处理中的一种基本方法和技术，尤其数字滤波技术使用广泛，数字滤波理论的研究及其产品的开发一直受到很多国家的重视。对自适应滤波算法的研究是当今自适应信号处理中最为活跃的研究课题之一。Windrow等于1967年提出的自适应滤波系统的参数能自动的调整而达到最优状况，而且在设计时，只需要很少的或根本不需要任何关于信号与噪声的先验统计知识。这种滤波器的实现差不多像维纳滤波器那样简单，而滤波器性能几乎如卡尔曼滤波器一样好。自适应滤波器与普通滤波器不同，它的冲激响应或滤波参数是随外部环境的变化而变化的,经过一段自动调节的收敛时间达到最佳滤波的要求。自适应滤波器本身有一个重要的自适应算法，这个算法可以根据输入、输出及原参量信号按照一定准则修改滤波参量，以使它本身能有效的跟踪外部环境的变化。因此，自适应数字系统具有很强的自学习、自跟踪能力和算法的简单易实现性。自适应滤波技术的核心问题是自适应算法的性能问题，提出的自适应算法主要有最小均方（LMS）算法、递归最小二乘（RLS）算法及相应的改进算法如：归一化（NLMS）算法、变步长（SVSLMS）算法、递归最小二乘方格形（RLSL）算法等。这些算法各有特点，适用于不同的场合。研究自适应算法是自适应滤波器的一个关键内容。递归最小二乘（RLS）算法是线性自适应滤波算法中最基本的两类算法之一，由于基于LMS准则的自适应滤波算法的收敛速度通常较慢，有些在调整过程种的延时也较大。为了克服LMS的算法，我们采用在每个时刻对所有已输入信号重估的平方误差之和最小这样的准则，即RLS算法。RLS算法复数乘法正比于2k，使其自适应速度更快。目前应用最多的是系统辨识、回波消除、自适应谱线增强、自适应信道均衡、语音线性预测、自适应天线阵等诸多领域。1.2国内外研究发展状况自适应滤波的基本理论通过几十年的发展已日趋成熟，近十几年来自适应滤波器的研究主要针对算法与硬件实现。算法研究主要是对算法速度和精度的改张少博：基于RLS算法自适应滤波器的设计-2-进，其方法大都采用软件C、MATLAB等仿真软件对算法的建模和修正。通常，自适应滤波器的硬件实现都是用DSP通用处理器（如TI的TMS320系列）。DSP器件采用改进的哈佛结构，具有独立的程序和数据空间，允许同时存取程序和数据，内置高速的硬件乘法器（MAC），增强的多级流水线。DSP具有的硬件乘法模块（MAC），专用的存储器以及适用于高速数据运行的总线结构，使DSP器件具有高速的数据运算能力。目前，用DSP器件处理数字信号已经成为电子领域的研究热点。在自适应信号处理领域，对于数据处理速度在几兆赫兹以内的，通用DSP器件也是首选。迟男等人在TMS320C32芯片上扩展EPROM和RAM，实现了30阶LMS自适应滤波器，使用的刀D转化器件为AD1674，最高采样频率为l00KHz。陆斌等人采用TMS320C30数字信号处理器与IMSA110专用滤波器并行处理的方法设计出了自适应滤波器并应用于直接序列的扩频接收系统1221。赵慧民等人在TMS320C31上实现了自适应权向量滤波器，完成了信号采样频率为80KHz的自适应滤波。在数据处理速度只要求在几兆赫兹以内的应用场合，这些用DSP实现的自适应滤波器能很好的满足系统实时的需求。在这种需求场合下，DSP具有不可媲美的性价比。但是随着信息化的进程加快和计算机科学与技术、信号处理理论与方法等的迅速发展，需要处理的数据量越来越大，对实时性和精度的要求越来越高。以迅速发展的移动通信技术为例，从IG时代只能传送语音的模拟通信，到2G时代的传送语音和数据的GSM、TDMA与CDMA1595，到2.5G时代传送语音、数据、图片、彩信MMS、简短视频、收发E-mail、网页浏览等的GPRS与CDMA2000lX，到目前正处于研发与测试阶段的能够传送图像、音乐、视频流等多种媒体形式，提供包括网页浏览、电话会议、电子商务的3G通信，以及目前正在研发与憧憬中的能够传送高质量流畅的视频流与多种实时流媒体业务的4G通信。常用的数字系统目标器件除了DSP处理器外还有专用集成电路（ASIC）、专用标准电路模块（ASSP）和现场可编程门阵列（FPGA）。ASIC和ASSP是专门针对完成某种数字信号处理算法的集成电路器件，因此其在性能指标、工作速度、可靠性和成本上优于DSP处理器。其优秀的工作性能主要源于特定的算法全部由ASSP和ASIC中的硬件电路完成。ASSP是半定制集成电路，在许多DSP算法的实现方面都优于DSP（数字信号处理器），但在功能重构，以及应用性修正方面缺张少博：基于RLS算法自适应滤波器的设计-3-乏灵活性；ASIC专用集成电路使用超大规模专用集成电路ASIC的实现方法是实用化的产品唯一可行的方法，只有使用IC，才有高可靠性和可接受的价格及体积功耗等。ASIC虽然有一定的可定制性，但开发周期长，而且有一个最小定制量，在实验室研制开发阶段，开发成本非常高。现代大容量、高速度的FPGA在可重配置的数字信号处理应用领域，特别是对于任务单一、算法复杂的前端数字信号处理运算，有独特的优势。另外，FPGA所具有的大规模并行处理能力和可编程的灵活性使得设计的系统能获得极高的处理性能，并且能够适应日益变化的标准、协议和性能需求。用FPGA实现自适应滤波器，国外起步比较早，发展也非常迅速。HesenerA.于1996年提出了用FPGA实现自适应滤波器的设想，并在FPGA上实现了处理速度可达SM的8阶8位FIR滤波器。WoolfriesN.等人用FPGA实现了自适应栈滤波器，并应用于图象处理。DawoodA.等人用FPGA开发了自适应FIR滤波器并与DSP处理器方案进行了比较研究。国内有一些关于自适应算法硬件实现的研究，但基本是针对自适应滤波器中的算法，如南开大学李国峰的博士论文用VHDL语言描述了正负数的运算问题和浮点数运算问题，完成了基于FIR的LMS自适应滤波器的硬件设计与逻辑综合。国防科学技术大学江和平等人讨论了自适应卡尔曼算法的简化，并完成了FPGA的设计。同济大学梁甲华等人重点讨论了编码方法在FPGA的技术问题。上海交通大学范瑜等人介绍了用VHDL语言实现并行延时LMS算法的自适应数字波束成形器的FPGA设计过程。而针对自适应格型结构采用FPGA硬件实现的文献报导很少，国内中国科学技术大学王显洁等人通过采用流水线结构和运算单元分时复用，提高了运算速度，能够满足实时性预测编码要求。1.3本文研究思路与主要工作在信号的传播路径中以及在信号处理过程中，都会引入噪声。噪声的引入影响了对真实信号的处理，有时候，较强的噪声会“遮盖”了信号，从而难以得到准确、稳定的真实信号。噪声对信号的污染在绝大多数情况下是不可避免的，因而，对噪声的消除和抑制是信号处理中极其重要的工作。在信号传输和处理过程中，最常见的噪声形式为白噪声、带限白噪声、高斯白噪声、瑞利分布白噪声。RLS（RecursiveLeastsquare）自适应滤波器常常用于噪声消除器的构建，本文介绍了自适应滤波器原理，对RLS自适应算法进行分析,最后用MATLAB对自适应滤波器进行了仿真和实现,并分析了该自适应滤波器的性能。张少博：基于RLS算法自适应滤波器的设计-4-第2章自适应滤波器理论基础2.1数字滤波器的基本概念凡是有能力进行信号处理的装置都可以称为滤波器。如果滤波器的输入和输出均为离散信号，称该滤波器为数字滤波器。当滤波器的输出信号为输入端的线性函数时，该滤波器称为线性滤波器，否则就称为非线性滤波器。一个典型的数字滤波器的框图如图2.1所示。图2.1数字滤波器设输入信号为nx，输出信号为ny，该数字滤波器可用以下差分方程来表示：1011)()(MiNiiiinybinxany(2-1)式中ia，ib称为滤波器系数。当0ib时，上式变为：10Miiinxany(2-2)这种滤波器称为全零点滤波器。如果0ia，0ib时，则称为全极点滤波器或递归滤波器。由上式，可知数字滤波器的传递函数为：MiiiMiizbzazH11011(2-3)其单位冲击响应函数为：zHznh1(2-4)张少博：基于RLS算法自适应滤波器的设计-5-iinxihnxnhny(2-5)如果当n0时，有nh=0，这样的滤波器系统称之为因果系统。如果冲激响应函数是有限长的，即：elseNnnhnh,00,(2-6)则称此滤波器为有限冲激响应FIR(FiniteImpulseResponse)滤波器，否则，称之为无限冲激响应IIR(InfiniteImpulseResponse)滤波器。如果nh满足如下条件：Cnhnnhn00,0(2-7)则称此滤波器是因果的，并且是稳定的。2.2自适应滤波器的原理自适应滤波器由参数可调的数字滤波器和自适应算法两部分组成（如图2.2所示）。参数可调数字滤波器可以是FIR数字滤波器或IIR数字滤波器,也可以是格型数字滤波器。输入信号nx通过参数可调数字滤波器后产生输出信号(或响应)ny,将其与参考信号（或称期望响应）nd进行比较，形成误差信号ne。ne（有时还要利用nx）通过某种自适应算法对滤波器参数进行调整,最终使ne的均方值最小。因此,自适应滤波器实际上是一种能够自动调整本身参数的特殊维纳滤波器。在设计时不需要事先知道关于输入信号和噪声的统计特性的知识,它能够在自己的工作过程中逐渐了解,或估计出所需的统计特性,并以此为依据自动调整自己的参数,以达到最佳滤波效果。一旦输入信号的统计特性发生变化,它又能够跟踪这种变化,自动调整参数,使滤波器性能重新达到最佳[1]。张少博：基于RLS算法自适应滤波器的设计-6-图2.2自适应滤波器一般结构2.3自适应滤波器的算法自适应算法主要根据滤波器输入的统计特性进行处理,它可能与滤波器的输入及其他数据有关,因此,存在开环算法和闭环算法。开环算法的控制输出仅取决于滤波器的输入和其他输入数据;闭环的控制输出则是滤波器输出及其他输入信号的函数。闭环控制利用输出反馈,它不但能在滤波器输入信号变化时保持最佳输出,而且还能在某种程度上补偿滤波元件参数的变化和误差及运算误差。但其存在稳定性问题和收敛速度不高的缺点;开环算法的优点是调整速度快,一般无稳定性问题,但通常它要求的计算量大且不能补偿元件参数误差和运算误差,因此,一般采用闭环算法。自适应算法是根据某种判断来设计滤波器的。该算法包括最小均方算法LMS、最小高阶均方算法LMF、最小平方算法OLS和递推最小算法RLS等。由WIDROW和HOFF提出来的最小均方误差算法,因具有计算量小、易于实现、不依赖模型和性能稳健等优点而被广泛采用。在自适应滤波器的实际操作中,一个必须注意的问题是自适应步长的选择。由LMS迭代算法可知,欲使其收敛,则自适