自适应滤波器的设计与实现

1、北京石油化工学院2013届优秀毕业论文自适应滤波器的设计与实现通信工程专业 熊洁指导教师 张宁 教授摘 要: 最小二乘(RLS)法是一种典型的有效的数据处理方法。它使各项实际观测值和计算值之间的差的平方乘以度量其精度的数值以后的和为最小。这就是著名的最小二乘法。递推最小二乘法(RLS)是最小二乘法的一类快速算法。本文基于matlab进行仿真，设计了M文件，研究不同输入和信噪比滤波器的性能。通过对比实验，分析收敛性和滤波效果，确定了滤波器正则化系数的选取。关键词: 自适应滤波器，最小二乘法，matlab仿真1 引言自适应滤波器属于现代滤波的范畴，它是40年代发展起来的自适应信号处理领域的一个重要

2、应用，自适应信号处理主要是研究结构可变或可调整的系统，可以通过自身与外界的接触来改善自身对信号处理的性能，通常这类系统是时变的非线性系统，可以自动适应信号传输的环境和要求，无须详细的知道信号的结构和实际知识，无须精确设计处理系统本身。自适应系统的非线性特性主要是由系统对不同的信号环境实现自身参数的调整来确定的。自适应系统的时变特性主要是由其自适应响应或自适应学习过程来确定的，当自适应过程结束和系统不再进行时，有一类自适应系统可成为线性系统，并称为线性自适应系统，因为这类系统便于设计且易于数学处理，所以实际应用广泛。本文研究的自适应滤波器就是这类滤波器。自适应滤波器是相对固定滤波器而言的，固定滤

3、波器属于经典滤波器，它滤波的频率是固定的，自适应滤波器的频率则是自动适应输入信号而变化的，所以其适用范围更广。在没有任何信号和噪声的先验知识的条件下，自适应滤波器利用前一时刻已获得的滤波器参数来自动调节现时刻的滤波器参数，以适应信号和噪声未知或随机变化的统计特性，从而实现最优滤波。近十几年来，自适应滤波理论和方法得到了迅速的发展，究其原因是因为自适应滤波器相比于其他一般的滤波器在滤波性能、设计实现的难易程度、对外部环境的复杂程度的适应能力和对系统先验统计知识的依赖程度等方面都显现出强大的优势。自适应滤波器具有很强的自学习、自跟踪能力和算法的简单易实现性，它在噪声量化信号的检测增强，噪声干扰的抵

4、消，通信系统的自适应均衡，图象的自适应增强复原以及未知系统的自适应参数辩识等方面都有广泛的应用。本文仅讨论自适应滤波器在噪声干扰的抵消方面的原理、算法及仿真。2 递归最小二乘法所谓的自适应滤波，就是利用前一时刻以获得的滤波器参数的结果，自动的调节现时刻的滤波器参数，以适应信号和噪声未知的或随时间变化的统计特性，从而实现最优滤波。自适应滤波器实质上就是一种能调节其自身传输特性以达到最优的维纳滤波器。自适应滤波器不需要关于输入信号的先验知识，计算量小，特别适用于实时处理。2.1递归最小二乘法原理递归最小二乘估计算法是以最小误差平方和为优化目标，它是在在最小二乘基础上推导出来的。用最小二乘法解决线性

5、滤波问题，这种方法不需要对滤波器输入信号的统计特性进行假设。为了说明最小二乘法的基本思想，假定有一组实数，它们分别取自时刻。要求构造一条曲线，这条曲线能够以某种最优方式拟合这些数据点。现用表示这条曲线与时间的函数关系。根据最小二乘法，“最优”拟合是使与，之差的平方和的代价函数最小。误差函数定义为 公式（1）式中， 为期望响应，是抽头输入向量，为抽头增益向量。递归最小二乘的代价函数定义如下： 公式（2）要完成RLS算法就是要找到抽头向量，使得代价函数最小。RLS算法引入了数加权遗忘因子。该遗忘因子的引入，使RLS算法能够对非平稳信号进行跟踪。粗略的说，的倒数可以用来衡量算法的记忆能力；而的特殊情

6、况，则对应于无限记忆。本文后续实验部分默认。式中是一个正实数，称为正则化参数。将这一项包含在代价函数中，一遍通过平滑作用来稳定递归最小二乘问题的解。将抽头输入向量的时间平均相关矩阵表示为 公式（3）横向滤波器抽头输入与期望响应之间的时间平均互相关向量为 公式（4）通过计算，得出递归最小二乘问题的正则方程可用矩阵形式写为 公式（5）将对应于的项与式（3）右边的求和项分开，可写出 公式（6）根据定义，式（6）右括号内的表达式等于相关矩阵。于是，可使得用于更新抽头输入相关矩阵的递归公式 公式（7）其中是相关矩阵的过去值，矩阵乘积在更新过程中起着修正的作用。式（7）的递归过程与初始条件无关。类似的，可

7、用式（3）导出抽头输入与期望响应之间互相关向量的更新公式 公式（8）为了按式（3-31）计算抽头权向量的最小二乘估计，必须确定相关矩阵的逆。然而在实际中，我们通常尽量避免这样做，因为这种运算非常耗时，特别是当抽头数很大时。为此，我们先引入一个著名的结果矩阵求逆引理。设和是两个正定矩，他们之间的关系为 公式（9）其中，是正定矩，是矩阵。根据矩阵求逆引理，可将的矩阵表示为 公式（10）假定相关矩阵是非奇异的，因而它可逆。我们对式（5）所表示的递归方程应用矩阵求逆引理，首先做如下设定 公式（11）然后将这些定义代入求逆引理，可得计算相关矩阵的递归方程如下 公式（12）为了方便计算，令 公式（13）和

8、 公式（14）用上面的定义，可将式（12）改写为 公式（15） 矩阵叫做逆相关矩阵，向量叫做增益向量。式（15）是RLS算法的Riccati方程。整理式（14），可得 公式（16）即： 公式（17）从式（15）可以看出，式（17）右边最后一行括号里的表达式等于。因此，我们可将式（17）简化为 公式（18）这一结论，连同，可以用来定义增益向量 公式（19）换句话说，增益向量可定义为经相关矩阵逆矩阵变换的抽头输入向量。下面导出更新抽头权向量最小二乘估计的递归公式。为此，用式（2）、式（7）和式（12）来表示抽头权向量n次迭代的最小二乘估计 公式（20）将式（20）右边第一项中用式（15）代替，可得

9、 公式（21）最后，应用等于增益向量，可使得更新抽头权向量的递归方程为 公式（22）其中 公式（23）是一个先验误差。内积表示基于时刻抽头权向量最小二乘估计旧值得期望响应的估值。2.2 RLS算法小结算法初始化 公式（24）对每一时刻，计算 公式（25） 公式（26） 公式（27） 公式（28）3.RLS自适应滤波器M文件的设计与仿真根据前面介绍的LS算法，用MATLAB设计RLS自适应滤波器。本文一共做了三个实验，下面进行仿真与分析。3.1 正弦波去噪实验原始信号为单一正弦波，加入方差为0.8的随机高斯噪声，滤波器输入信号为，输出信号为。运行结果如下： 图1 误差信号输出信号和抽头权值w变化

10、图 图2 误差信号输出信号和抽头权值w变化图为0.8的随机高斯噪声，得到滤波器输入信号，从频谱图可以看出，包含原始信号有很多频率。图3中红色部分为经过RLS滤波器后的输出信号，蓝色部分为误差信号。可以看出误差信号较小。由上可知，滤波效果较好。3.2 滤波器正则化参数的确定RLS算法有两个特殊的可变参数：抽头输入数据的信噪比，这个量由运行的条件决定；正则化参数。以下做了对比实验，在输入为正弦波，控制加性高斯噪声的均值，相同均值的随机的噪声都保持不变。来实现不同信噪比，分析正则化参数的在最佳选择。4.2.1 高信噪比原始输入信号是幅值为1的正弦波，噪声方差为0.1。实验组1：实验组2：运行得到图4

11、。 图3 误差信号输出信号和抽头权值w变化图 图4 误差信号输出信号和抽头权值w变化图实验组2程序与实验组1大部分相同，改变Delta=0.0001运行得到图5。 图5 误差信号输出信号和抽头权值w变化图 图6 误差信号输出信号和抽头权值w变化图当信噪比较高时，实验组1抽头系数在4000开始收敛，实验组2抽头系数在500开始收敛。两组误差信号都很小。不难看出，实验组2收敛速度更快。4.2.2 低信噪比原始输入信号是幅值为1的正弦波，噪声方差为0.8。实验组1：实验组2：实验程序与高信噪比时大部分一样，变动的语句：b=0.8*randn(1,N)。其它部分这里就不再重复。当信噪比较低时，实验组1

12、抽头系数在1800开始收敛，误差信号较小。实验组2抽头系数在500开始收敛，但是误差信号波动较大。综合考虑。实验组1的全局性能较好。4.2.3 结论 综合上述两组对比实验，当抽头输入噪声电平低（即输入信噪比较高）时，RLS算法的正则化系数取较小的正常数收敛速度更快，滤波效果也较好。当抽头输入噪声电平高（即输入信噪比较低）时，RLS算法的正则化系数取较大的正常数收敛速度较快，滤波效果也较好，全局性能更好。图 7 误差信号输出信号和抽头权值w变化图3.3 输入信号不同对滤波效果的影响3.3.1 输入信号为周期信号原始信号输入为周期三角波，噪声方差为2。运行得到图8到图10。 图8 原始信号时域和频

13、域图 图9 滤波器输误差信号和输出信号频域图 图10 滤波器输误差信号和输出信号频域图 图11 滤波器输误差信号和输出信号频域图如图7从时域观察，原始信号峰值为10，噪声峰值为2。滤波之后误差信号值接近于0。从频域上看，接近原始信号频谱。滤波效果较好。3.3.2 输入信号为非周期信号本实验原始输入信号为纯背景音频信号，格式为wav,将其放在matlab文件，存为“2.wav”。所加噪声方差为2。运行得到图11到图14：图 11 滤波器输误差信号和输出信号频域图 图 12 滤波器输误差信号和输出信号频域图 图 13 滤波器输误差信号和输出信号频域图图 14 滤波器输误差信号和输出信号频域图如图1

14、1，从时域分析，背景音乐信号加入相同均值为2噪声，误差信号幅值较三角波实验大。如图14，从频域上看，滤波之后，在100200一级800900之间含有少量噪声频率。滤波效果不如输入为周期信号。4结束语通过对比实验，分析收敛性和滤波效果，确定了滤波器正则化系数的选取，信噪比高正则化系数选取小正常数，反之取大正常数。对不同类型滤波器输入信号独立实验，周期型号和非周期信号的滤波器效果都较好，但周期信号的滤波效果更好。但从程序的运行时间上看，其复杂度高，计算量比较大。 Abstract: Least squares (RLS) method is a kind of typical effective

15、data processing method. It makes the actual observed value multiplied by the square measure of the difference between calculated value and the accuracy of numerical and later to a minimum. This is the famous least-square method. Recursive least squares (RLS) is a fast algorithm of the least-square m

16、ethod. Simulation based on matlab, this paper designs the M file, the different performance of the filter input and SNR. Through contrast experiment, convergence analysis and filtering effect, the filter selection of regularization coefficient is determined. Key words：Adaptive filtering, the least s

17、quare method, emulation verify by MATLAB参 考 文 献1 于洋. 一种改进的LMS自适应滤波器设计D.本科学位论文.北京：北京石油化工学院，20092 丁玉美. 数字信号处理M.西安：西安电子科技大学出版社，20023 王姣. 自适应滤波器的MATLAB仿真D. 本科学位论文.北京：北京石油化工学院，20024 薛小梅. 一种基于LMS算法的自适应滤波器设计D. 本科学位论文.北京：北京石油化工学院，20065 吴燕. 一种自适应滤波器的设计与仿真J.科技广场，2010，第7期：59-616 程玉柱. 自适应滤波器的算法比较研究J.大庆师范学报，2008，第2期：23-267 张立萍. RLS自适应滤波器的Matlabe设计与仿真J.赤峰学院学报，2011，第4期：25-268 蒙俊甫. LMS自适应滤波器的仿真与实现J.设计参考，2009，第9期：67-739 程玉柱. 自适应滤波器的算法比较研究J.大庆师范学报，2008，第2期：2