基于RLS算法自适应滤波器

基于RLS算法自适应滤波器的设计摘要自适应滤波器是统计信号处理的一个重要组成局部。在实际应用中，由于没有充足的信息来设计固定系数的数字滤波器，或者设计规那么会在滤波器正常运行时改变，因此需要研究自适应滤波器。但凡需要处理未知统计环境下运算结果所产生的信号或需要处理非平稳信号时，自适应滤波器可以提供非自适应方法所不可能提供的新的信号处理能力。而且其性能通常远优于用常方法设计的固定滤波器。本文从自适应滤波器研究的意义入手，介绍了自适应滤波器的根本理论思想，具体阐述了自适应滤波器的根本原理、算法及设计方法。自适应滤波器的算法是整个系统的核心。对RLS算法自适应滤波器做了详细的介绍，采用改良的RLS算法设计自适应滤波器，并采用MATLAB进行仿真，通过实验结果来表达该滤波器可以根据信号随时修改滤波参数，到达动态跟踪的效果，使滤波信号更接近于原始信号。关键词：自适应滤波器，RLS算法，噪声消除，FIR第1章绪论1.1课题研究意义和目的滤波技术是信号处理中的一种根本方法和技术，尤其数字滤波技术使用广泛，数字滤波理论的研究及其产品的开发一直受到很多国家的重视。对自适应滤波算法的研究是当今自适应信号处理中最为活泼的研究课题之一。Windrow等于1967年提出的自适应滤波系统的参数能自动的调整而到达最优状况，而且在设计时，只需要很少的或根本不需要任何关于信号与噪声的先验统计知识。这种滤波器的实现差不多像维纳滤波器那样简单，而滤波器性能几乎如卡尔曼滤波器一样好。自适应滤波器与普通滤波器不同，它的冲激响应或滤波参数是随外部环境的变化而变化的,经过一段自动调节的收敛时间到达最正确滤波的要求。自适应滤波器本身有一个重要的自适应算法，这个算法可以根据输入、输出及原参量信号按照一定准那么修改滤波参量，以使它本身能有效的跟踪外部环境的变化。因此，自适应数字系统具有很强的自学习、自跟踪能力和算法的简单易实现性。自适应滤波技术的核心问题是自适应算法的性能问题，提出的自适应算法主要有最小均方〔LMS〕算法、递归最小二乘〔RLS〕算法及相应的改良算法如：归一化〔NLMS〕算法、变步长〔SVSLMS〕算法、递归最小二乘方格形〔RLSL〕算法等。这些算法各有特点，适用于不同的场合。研究自适应算法是自适应滤波器的一个关键内容。递归最小二乘〔RLS〕算法是线性自适应滤波算法中最根本的两类算法之一，由于基于LMS准那么的自适应滤波算法的收敛速度通常较慢，有些在调整过程种的延时也较大。为了克服LMS的算法，我们采用在每个时刻对所有已输入信号重估的平方误差之和最小这样的准那么，即RLS算法。RLS算法复数乘法正比于，使其自适应速度更快。目前应用最多的是系统辨识、回波消除、自适应谱线增强、自适应信道均衡、语音线性预测、自适应天线阵等诸多领域。e(n)1.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〔数字信号处理器〕，但在功能重构，以及应用性修正方面缺乏灵活性；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对自适应滤波器进行了仿真和实现,并分析了该自适应滤波器的性能。第2章自适应滤波器理论根底2.1数字滤波器的根本概念但凡有能力进行信号处理的装置都可以称为滤波器。如果滤波器的输入和输出均为离散信号，称该滤波器为数字滤波器。当滤波器的输出信号为输入端的线性函数时，该滤波器称为线性滤波器，否那么就称为非线性滤波器。一个典型的数字滤波器的框图如图2.1所示。图2.1数字滤波器设输入信号为，输出信号为，该数字滤波器可用以下差分方程来表示：(2-1)式中，称为滤波器系数。当时，上式变为：(2-2)这种滤波器称为全零点滤波器。如果，时，那么称为全极点滤波器或递归滤波器。由上式，可知数字滤波器的传递函数为：(2-3)其单位冲击响应函数为：(2-4)(2-5)如果当<0时，有=0，这样的滤波器系统称之为因果系统。如果冲激响应函数是有限长的，即：(2-6)那么称此滤波器为有限冲激响应FIR(FiniteImpulseResponse)滤波器，否那么，称之为无限冲激响应IIR(InfiniteImpulseResponse)滤波器。如果满足如下条件：(2-7)那么称此滤波器是因果的，并且是稳定的。2.2自适应滤波器的原理自适应滤波器由参数可调的数字滤波器和自适应算法两局部组成〔如图2.2所示〕。参数可调数字滤波器可以是FIR数字滤波器或IIR数字滤波器,也可以是格型数字滤波器。输入信号通过参数可调数字滤波器后产生输出信号(或响应),将其与参考信号〔或称期望响应〕进行比拟，形成误差信号。〔有时还要利用〕通过某种自适应算法对滤波器参数进行调整,最终使的均方值最小。因此,自适应滤波器实际上是一种能够自动调整本身参数的特殊维纳滤波器。在设计时不需要事先知道关于输入信号和噪声的统计特性的知识,它能够在自己的工作过程中逐渐了解,或估计出所需的统计特性,并以此为依据自动调整自己的参数,以到达最正确滤波效果。一旦输入信号的统计特性发生变化,它又能够跟踪这种变化,自动调整参数,使滤波器性能重新到达最正确[1]。图2.2自适应滤波器一般结构2.3自适应滤波器的算法自适应算法主要根据滤波器输入的统计特性进行处理,它可能与滤波器的输入及其他数据有关,因此,存在开环算法和闭环算法。开环算法的控制输出仅取决于滤波器的输入和其他输入数据;闭环的控制输出那么是滤波器输出及其他输入信号的函数。闭环控制利用输出反应,它不但能在滤波器输入信号变化时保持最正确输出,而且还能在某种程度上补偿滤波元件参数的变化和误差及运算误差。但其存在稳定性问题和收敛速度不高的缺点;开环算法的优点是调整速度快,一般无稳定性问题,但通常它要求的计算量大且不能补偿元件参数误差和运算误差,因此,一般采用闭环算法。自适应算法是根据某种判断来设计滤波器的。该算法包括最小均方算法LMS、最小高阶均方算法LMF、最小平方算法OLS和递推最小算法RLS等。由WIDROW和HOFF提出来的最小均方误差算法,因具有计算量小、易于实现、不依赖模型和性能稳健等优点而被广泛采用。在自适应滤波器的实际操作中,一个必须注意的问题是自适应步长的选择。由LMS迭代算法可知,欲使其收敛,那么自适应步长的取值范围为0<<1/,自适应时间常数和步长的关系为=1/()[2]。2.4RLS自适应滤波器基于LMS准那么的自适应滤波算法的收敛速度通常较慢，有些在调整过程种的延时也较大。为了克服LMS的算法，我们采用在每个时刻对所有已输入信号重估的平方误差之和最小这样的准那么，即RLS自适应滤波算法。从物理概念上说，这是一种在现有约束条件下利用了最多可利用信息的原那么。 RLS自适应滤波器的结构框图如图2.3所示：图2.3RLS横向自适应滤波器根本RLS自适应算法所遵循的准那么是确定这样的W，它使=－〔2-8〕的加权平方和：＝〔2-9〕最小。其中：＝〔2-10〕＝〔2-11〕为略小于1的值，，称为遗忘因子。参加这个因资额的物理含义是在权系数W所用到的输入信号中，我们对时间较近的数据加以较大的权来考虑，时间教前的数据其权按指数规律减小。这样可使算法更能反映当前的情况，从而加强对信号统计特性有缓慢变动时的适应性。一般取0.95～0.9995，这种加权的方式为指数权。 为了使加权平方和最小，令：＝0 〔2-12〕即：＝＝－2＝0 〔2-13〕因而有：＝ 〔2-14〕假设令：＝ 〔2-15〕=〔2-16〕那么式2-14可写成：＝ 〔2-17〕如果是非奇异的，那么：＝〔2-18〕上式就是求解滤波器参数的公式，它也是最正确滤波器系数的维纳〔Winear〕方程，只是应注意这里的随时间而改变。但式1-48来调整有两处不变：第一，需要矩阵求逆及矩阵乘法等运算，因而计算量较大；第二，与预测误差之间没建立关系，不能实现直接由预测误差来调整滤波器参数的要求。不过我们注意到预测误差由＝－ 〔2-19〕表示。利用此表达式，可以将式2-27的改写为：＝ ＝＋〔2-20〕将上式带入式2-18，＝＋ ＝＋ 〔2-21〕为了简化第一项，并建立与之间的关系，认为时刻及其以前时刻的滤波器参数相同，即：〔2-22〕那么有：＝＝〔2-23〕为了简化，我们可以认为遗忘因子，这相当于只有本时刻的结果被记忆下来，而以前时刻的结果全部遗忘，于是可写为：＝＝ 〔2-24〕将式2-23和2-24带入2-22，可得：＝＋〔2-25〕式2-36描述了一个滤波器参数受其输入误差控制的自适应滤波算法，被称为递推最小二乘法〔RecursiveLeastSquares〕。为了实现递推计算，还应解决逆矩阵的递推计算问题。这里，我们引入一个著名的引理——矩阵求逆引理，该引理如下：假设是非奇异阵，那么＝ 〔2-26〕这里只要用左乘上式右边，并严正其结果等于单位阵即可。 由的定义式2-12可导出：＝＋ 〔2-27〕利用矩阵求逆引理，得：＝－ 〔2-28〕这样我们可以利用随的迭代式而推得的迭代式。我们将根本RLS的自适应滤波算法综合如下： 初始化步骤：(对于＝0) 令， 运算步骤：对于＝1到所需的终了时刻＝〔2-29〕－ 〔2-30〕=＋ 〔2-31〕从结构来看RLS自适应滤波器的自适应是通过对输入数据进行一定的算法实现的，所以这种结构是“开环〞的。RLS算法中的与LMS算法中的作用相同，但为标量，而那么是随而变的矩阵的逆，这说明不同时刻的每个元素的调整量均随新进的数据的不同步长因子做调整，而不是统一的用同一个因子来调整，这表征了调整的精细性及新信息数据利用的充分性。RLS算法复数乘法正比于，使其自适应速度更快[3]。第3章基于RLS算法自适应滤波器的实现3.1Matlab仿真自适应滤波器在许多场合,一个输入信号往往包含有周期性信号和宽带成分,而周期性信号是期望得到的。如图4.1所示是一个自适应噪声消除滤波器(ANC)的原理图,输入是带有噪声的正弦波,它能够通过自适应调节,别离出信号中所包含的周期性成分和随机成分,从噪声中复原出正弦波。其原理是当周期信号和噪声混合的输入信号被延长一定时间后,其中的周期信号成分是高度相关的,但根据高斯理论的推断,噪声信号是不相关的。于是自适应滤波器就会减小输出信号中噪声的能量,产生周期信号的最正确估计信号。周期信号和噪声都是时变信号,因此滤波器必须根据输入信号的特性适应这种变化,决定权值的选取,最终使得输出信号的能量最低,这样就从某种程度上消除了噪声[4]。图3.1自适应噪声消除滤波器原理框图下面用Matlab中的Simulink工具对自适应滤波器进行模拟仿真。如图3.2所示,是RLS设置参数：图3.2RLS滤波器设置参数即递推最小二乘自适应滤波器的FIR阶数为32，存储指数的权重因子为1.0，filter—taps的初值为0.0，初始输入的估计方差为0.1。Simulink仿真模型界面图，如图3.3所示。图3.3RLS自适应滤波器滤除噪声仿真模型界面图3.2自适应滤波器性能分析观察显示结果,其中TimeScope的显示，如图3.4所示。图3.4FIR阶数位32的波形显示界面图从图3.4中可以看出,第一个显示器中显示的信号为周期信号,代表有用信号;第二个显示器显示的是被噪声干扰后的周期信号;第三个显示器显示的是经过所设计的自适应滤波器后的波形。第一个显示器中也含有滤波后的波形以便与原周期信号进行比拟。第四个是噪声信号。从图3.4中可以分析出,一开始输出信号为0,经过N次迭代后,自适应滤波器慢慢调整权值使输出信号逼近原周期信号,最后与周期信号根本重合。图3.5和图3.6分别是FIR阶数位16和64的仿真输出所示。图3.5阶数16波形显示界面图图3.5阶数64波形显示界面图将图3.3、图3.4及图3.5中的第一个显示器中的波形比拟可以发现,当阶数位16时，自适应速度慢而且滤波效果差，噪声较大，当阶数位64时，自适应的滤波器大幅提高，但滤波效果较阶数位32时就差了许多，因此，为了提高滤波器的滤波效果，在设计RLS自适应滤波器时阶数尽量选择32附近。如图3.7和图3.8分别的滤波后的幅度响应和功率频谱。图3.7幅度响应图3.8功率频谱经过N次的叠加，自适应滤波器的幅度和功率均趋于平稳。通过上述仿真结果可以看出，基于MATLAB的自适应滤波器滤除噪声系统能够有效地从噪声中恢复出原始信号，但在实际应用中，要注意参考信号与噪声信号的相关性，相关性越大，自适应噪声