现代信号处理论文.doc

《现代信号处理论文.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、递归最小二乘法（RLS）自适应滤波算法摘要所谓自适应，从通俗意义上讲，就是这种滤波器能够根据输入信号统计特性的变化自动调整其结构参数，以满足某种最佳准则的要求。自适应滤波器所采用的最佳准则由最小均方误差准则、最小二乘准则、最大信噪比准则和统计监测准则等。自适应滤波理论和技术是统计信号处理和非平稳随机信号处理的主要内容，它可以在无需先验知识的条件下，通过自学习适应或跟踪外部环境的非平稳随机变化，并最终逼近维纳滤波和卡尔曼滤波的最佳滤波性能。因而，自适应滤波不但可以用来检测确定性信号，而且可以检测平稳的或非平稳的随机信号。自适应技术应用包

2、括自适应谱线检测增强与谱估计方法、自适应噪声抵消技术、自适应均衡技术、只适用阵列处理与波束形成以及自适应神经网络信号处理等内容。关键词：递归最小二乘法；自适应滤波；滤波器设计；自适应算法；1引言滤波可分为经典滤波和现代滤波。经典滤波要求已知信号和噪声的统计特性，如维纳滤波和卡尔曼滤波。现代滤波则不要求己知信号和噪声的统计特性，如自适应滤波。自适应滤波的原理就是利用前一时刻己获得的滤波参数等结果，自动地调节当前时刻的滤波参数，从而达到最优化滤波。自适应滤波一般包括3个模块：滤波结构、性能判据和自适应算法。其中，自适应滤波算法一直是人们的

3、研究热点，包括线性自适应算法和非线性自适应算法，非线性自适应算法具有更强的信号处理能力，但计算比较复杂，实际应用最多的仍然是线性自适应滤波算法，线性自适应滤波算法的种类很多，有LMS自适应滤波算法、R路自适应滤波算法、变换域自适应滤波算法、仿射投影算法、共扼梯度算法等。自适应滤波器主要包括滤波器的结构和自适应算法两部分，这两部分不同的变化与组合，可以导出许多不同形式的自适应滤波器。所谓的自适应滤波，就是利用前一时刻以获得的滤波器参数的结果，自动地调节现时刻的滤波器参数，以适应信号和噪声未知的或随时间变化的统计特性，从而实现最优滤波。自

4、适应滤波器实质上就是一种能调节其自身传输特性以达到最优的维纳滤波器。自适应滤波器不需要关于输入信号的先验知识，计算量小，特别适用于实时处理。2主要算法和程序8RLS算法是最小二乘算法的一类快速算法，观察一个平稳输入信号输入的自适应系数在一点时间内输出误差信号的平均功率（时间平均），将该平均功率是否为最小作为测量自适应系统是否最佳的准则。它包括时间递归最小二乘法（TRLS）与阶数递归最小二乘法（ORLS）两大类。TRLS主要应用于动态系统辨识和在线估计，因此在这里我们采用时间递归的最小二乘法TRLS。通常来说，RLS自适应算法具有快速收

5、敛性，它不像LMS类自适应算法的收敛性对输入信号线管矩阵参数很灵敏，这是低的鲁棒性和较高计算复杂度来得到的。为了提高RLS自适应算法的鲁棒性，目前最新成果是利用QR分解法。QR分解算法由于采用了正交三角化过程，大大提高了自适应算法滤波数值的稳定性，而且具有良好的收敛性和跟踪特性。在数学上实现正交三角化过程的技术有Gram-Schmidt正交化算法与Givens旋转方法。以下是利用RLS算法的实现过程。与一般的最小二乘方法不同，这里考虑一种指数加权的最小二乘方法。使用指数加权的误差平方和作为代价函数，即有估计误差定义为：加权误差平方和的

6、完整表达式为：根据定义，易得到递推估计公式：利用矩阵的求逆公式，可得到逆矩阵的递推公式进一步化简后得综上所述，可得到RLS直接算法：步骤1：初始化：8步骤2：更新：n=1,2,….。其中为一个很小的值。的值越小，相关矩阵初始值在的计算中所占比重越小，若取值过大，就会严重影响RLS算法的收敛速度及收敛结果，这一点是在应用RLS算法时必须注意的。以下是RLS自适应滤波器设计实现过程。1.生成正弦波载波原始信号是一个正弦波波信号，其波形如图1所示，其实现的代码如下。signal=sin(2*pi*0.055*(0:1000-1)');fig

7、ure(1);plot(0:199,signal(1:200),'LineWidth',2);set(gca,'FontSize',30);grid;axis([0200-22]);title('fontsize{30}正弦波载波');8图12.背景噪声模型背景噪声是白噪声，如图2所示。其实现的代码如下。nvar=1.0;%噪声方差noise=randn(1000,1)*nvar;%白噪声figure(2);plot(0:999,noise,'LineWidth',2);set(gca,'FontSize',30);title('

8、fontsize{30}背景噪声');grid;axis([01000-44]);8图23.输入自适应滤波器的信号与噪声输入的噪声是由正弦波信号与白噪声叠加后的结果，其波形如图3所示。其实现的代码如下。nfilt=fir