自适应滤波器毕业设计论文设计.docx
《自适应滤波器毕业设计论文设计.docx》由会员分享,可在线阅读,更多相关《自适应滤波器毕业设计论文设计.docx(12页珍藏版)》请在冰豆网上搜索。
![自适应滤波器毕业设计论文设计.docx](https://file1.bdocx.com/fileroot1/2022-10/25/ffc12cb9-9547-4c4d-a98e-835a58ded1b3/ffc12cb9-9547-4c4d-a98e-835a58ded1b31.gif)
自适应滤波器毕业设计论文设计
数字信号处理课程要求论文
基于LMS的自适应滤波器设计与应用
学院名称:
专业班级:
学生:
学号:
2013年6月
摘要自适应滤波在统计信号处理领域占有重要地位,自适应滤波算法直接决定着滤波器性能的优劣。
目前针对它的研究是自适应信号处理领域中最为活跃的研究课题之一。
收敛速度快、计算复杂性低、稳健的自适应滤波算法是研究人员不断努力追求的目标。
自适应滤波器是能够根据输入信号自动调整性能进展数字信号处理的数字滤波器。
作为比照,非自适应滤波器有静态的滤波器系数,这些静态系数一起组成传递函数。
研究自适应滤波器可以去除输出信号中噪声和无用信息,得到失真较小或者完全不失真的输出信号。
本文介绍了自适应滤波器的理论根底,重点讲述了自适应滤波器的实现结构,然后重点介绍了一种自适应滤波算法最小均方误差〔LMS〕算法,并对LMS算法性能进展了详细的分析。
最后本文对基于LMS算法自适应滤波器进展MATLAB仿真应用,实验明确:
在自适应信号处理中,自适应滤波信号占有很重要的地位,自适应滤波器应用领域广泛;另外LMS算法有优也有缺点,LMS算法因其鲁棒性强特点而应用于自回归预测器。
关键词:
自适应滤波器,LMS算法,Matlab,仿真
1.引言
滤波技术在当今信息处理领域中有着极其重要的应用。
滤波是从连续的或离散的输入数据中除去噪音和干扰以提取有用信息的过程,相应的装置就称为滤波器。
滤波器实际上是一种选频系统,他对某些频率的信号予以很小的衰减,使该局部信号顺利通过;而对其他不需要的频率信号予以很大的衰减,尽可能阻止这些信号通过。
滤波器研究的一个目的就是:
如何设计和制造最优的(或最优的)滤波器。
Wiener于20世纪40年代提出了最优滤波器的概念,即假定线性滤波器的输入为有用信号和噪音之和,两者均为广义平稳过程且己知他们的二阶统计过程,如此根据最小均方误差准如此(滤波器的输出信号与期望信号之差的均方值最小)求出最优线性滤波器的参数,称之为Wiener滤波器。
同时还发现,在一定条件下,这些最优滤波器与Wiener滤波器是等价的。
然而,由于输入过程取决于外界的信号、干扰环境,这种环境的统计特性常常是未知的、变化的,因而不能满足上述两个要求,设计不出最优滤波器。
这就促使人们开始研究自适应滤波器。
自适应滤波器由可编程滤波器(滤波局部)和自适应算法两局部组成。
可编程滤波器是参数可变的滤波器,自适应算法对其参数进展控制以实现最优工作。
自适应滤波器的参数随着输入信号的变化而变化,因而是非线性和时变的。
2.自适应滤波器的根底理论
所谓自适应滤波,就是利用前一时刻已获得的滤波器参数等结果,自动地调节现时刻的滤波器参数,以适应信号和噪声未知的或随时间变化的统计特性,从而实现最优滤波。
所谓“最优〞是以一定的准如此来衡量的,最常用的两种准如此是最小均方误差准如此和最小二乘准如此。
最小均方误差准如此是使误差的均方值最小,它包含了输入数据的统计特性,准如此将在下面章节中讨论;最小二乘准如此是使误差的平方和最小。
自适应滤波器由数字结构、自适应处理器和自适应算法三局部组成。
数字结构是指自适应滤波器中各组成局部之间的联系。
自适应处理器是前面介绍的数字滤波器〔FIR或IIR〕,所不同的是,这里的数字滤波器是参数可变的。
自适应算法如此用来控制数字滤波器参数的变化。
自适应滤波器可以从不同的角度进展分类,按其自适应算法可以分为LMS自适应滤波器、RLS自适应滤波器等等。
自适应滤波器实际上是一种能够自动调整本身参数的特殊维纳滤波器,在设计时不需要预先知道关十输入信号和噪声的统计特性,它能够在工作过程中逐步了解或估计出所需的统计特性,并以此为依据自动调整自身的参数,以达到最优滤波效果。
图2.1自适应滤波器的原理图
自适应滤波器的特性变化是由自适应算法通过调整可编程滤波器系数来实现的。
图2.1给出了自适应滤波器的一般结构,其中输入信号二x(n)通过可编程滤波器后产生输出信号(或响应)y(n),将其与参考信号d(n)进展比拟,形成误差信号e(n),并以此通过某种自适应算法对滤波器参数进展调整,最终使得e(n)的均方值最小。
利用抽头延迟线做成的横向滤波结构的自适应滤波器,通称为自适应横向滤波器(或自适应FIR滤波器)。
它是研究所有自适应滤波算法的根本结构,由于其结构简单、本钱较低,也是工程领域最常用的一种自适应滤波器。
图2.2单输入自适应滤波器
自适应横向滤波器的结构图如图2.2所示,为可调节抽头权系数表示在n时刻的系数值。
它利用正规直接形式实现全零点传输函数,而不采用反应调节。
权系数的调节过程是首先自动调节滤波器系数的自适应训练步骤,然后利用滤波系数加权延迟抽头上的信号来产生输出信号,将输出信号与期望信号进展比照,所得的误差值通过一定的自适应控制算法再用来调整权值,以保证滤波器处在最优状态,其抽头加权系数集正好等于它的冲激响应,达到实现滤波的目的。
自适应递归滤波器是指零点和极点都能自适应调整的滤波器,它的传输函数中分子和分母通常具有独立的迭代步长因子。
对有些利用横向滤波器实现时需要数百个甚至上千个抽头系数,可以考虑使用自适应递归滤波器。
缺点是递归滤波器要求对极点的稳定性进展监视,而且收敛速度很慢。
采用得最多的自适应递归滤波器结构是如图2-3所示的直接形式结构,自适应滤波器传输函数分母的系数为,传输函数分子的系数为,其中N和M分别是自适应滤波器分母和分子的阶数。
图2.3自适应递归滤波器
3.自适应滤波算法
3.1最小均方误差〔LMS〕算法
LMS算法是基于最小均方误差准如此(MMSE)的维纳滤波器和最陡下降法提出的,依据输入信号在迭代过程中估计梯度矢量,并更新权系数以达到最优的自适应迭代算法。
这算法不需要计算相应的相关函数,也不需要进展矩阵运算。
自适应滤波器最普通的应用就是横向结构。
滤波器的输出信号y(n)是
〔3-1〕
T表示转置矩阵,n是时间指针,N是滤波器次数。
这个例子就是有限脉冲响应滤波器的形式,为x(n)和w(n)两个矩阵卷积。
这种自适应算法使用误差信号
〔3-2〕
为了方便起见,将上述式子表示为向量形式,如此上述式子表示为:
〔3-3〕
误差序列可写为
〔3-4〕
其中d(n)是期望信号,y(n)是滤波器的输出。
使用输入向量x(n)和e(n)来更新自适应滤波器的最小化标准的相关系数。
显然,自适应滤波器控制机理是用误差序列e(n)按照某种准如此和算法对其系数{wi(n)},i=1,2,…,N进展调节的,最终使自适应滤波的目标(代价)函数最小化,达到最优滤波状态。
本节所用的标准是最小均方误差(MSE)即
〔3-5〕
E[]表示算子期望。
假设公式中的y(n)被公式(3-4)取代,式(3-5)就可以表示为
〔3-6〕
是自相关矩阵,是输入信号的自相关矩阵。
是互相关向量,也指出了期望信号和输入信号向量的相互关矢量。
由式(3-6)可见,自适应滤波器的代价函数是延迟线抽头系数的二次函数。
当矩阵R和矢量P己知时,可以由权系数矢量w直接求其解。
最优解最小化MSE,源自解这个公式
〔3-7〕
将式〔3-6〕对w求其偏导数,并令其等于零,假设矩阵R满秩(非奇异),可得代价函数最小的最优滤波系数:
〔3-8〕
这个解称为维纳解,即最优滤波系数值。
因为均方误差(MSE)函数是滤波系数w的二次方程,由此形成一个多维的超抛物面,这好似一个碗状曲面又具有唯一的碗底最小点,通常称之为自适应滤波器的误差性能曲面。
当滤波器工作在平稳随机过程的环境下,这个误差性能曲面就具有固定边缘的恒定形状。
自适应滤波系数的起始值{wi(0)},i=1,2,…,N是任意值,位于误差性能曲面上某一点,经过自适应调节过程,使对应于滤波系数变化的点移动,朝碗底最小点方向移动,最终到达碗底最小点,实现了最优维纳滤波。
对于LMS算法梯度v(n)通过假设平方误差。
2(n)作为公式〔3-7〕的MSE来预测。
因此,梯度预测可以单一化表示为:
〔3-9〕
在实际应用中,2u经常用来代替u。
瞬间梯度预测产生的Widrow一HoffLMS算法,w(n)为自适应滤波器在n时刻的滤波系数或权矢量。
按照最陡下降法调节滤波系数,如此在n+1时刻的滤波系数或权矢量w(n+l)可以用如下简单递归关系来计算:
〔3-10〕
u是自适应步长来控制稳定性和收敛率。
这种瞬时估计是无偏的,因为它的期望值E[]等于最陡下降法的梯度矢量。
以任意初始向量w(0)来开始,向量w(n)集中在最优解决方法w0,假设选u
〔3-11〕
为矩阵R的最大特征值,受限制于
〔3-12〕
Tr[]为指示矩阵的轨迹,是平均输入功率。
对于自适应信号处理应用,最重要的实际考虑是收敛速度,决定滤波器跟踪不稳定型号的能力。
总体来说,权向量要获得收敛只有当最缓慢的权集中一点。
这个最慢的时间
〔3-13〕
这个指出时间连续相反的以u的比例收敛,并且依靠输入矩阵的自相关特征值。
具有全异的特征值,规定时间是受最慢模式的限制。
以梯度预测为根底的自适应导致噪声矩阵的权向量,因此会有性能的损失。
这个自适应处理的噪声导致稳态权向量随意的改变为最适宜的权向量。
稳态权向量的精度通过超额的最小均方误差来测量。
这个LMS算法超过EMS的是
〔3-14〕
是MSE在稳态的最小值。
公式(3-15)和(3-16)产生LMS算法根本协定:
为了在稳态获得高精度(低超自适应滤波算法与应用研究额MSE),需要u的最小值,但是也会降低收敛率。
后面会有进一步关于LMS算法特征的讨论。
对于N维更新u*e(n)是常数,误差信号e(n)乘以u得到u*e(n)。
这个常数首先计算,然后乘以x(n)来更新w(n)。
自适应LMS算法如同最陡下降法,利用时间n=0的滤波系数矢量为任意的起始值w(n),然后开始LMS算法的计算,其步骤如下:
l)由现在时刻n的滤波器滤波系数矢量估值w(n),输入信号矢量x(n)与期望信号d(n),计算误差信号e(n):
〔3-15〕
2)利用递归法计算滤波系数矢量的更新估值。
3)将时间指数n增加1,回到第一步骤,重复上述计算步骤,一直到达稳定状态为止。
由此可见,自适应LMS算法简单,它既不需要计算输入信号的相关函数,又不要求矩阵之逆。
因而得到了广泛的应用。
最小均方差(LMS)算法的性能分析
LMS算法的性能准如此是采用瞬时平方误差性能函数|e(k)|2代替均方误差性能函数E{|e(k)|2},其实质是以当前输出误差、当前参考信号和当前权系数求得下个时刻的权系数。
其输出信号输出误差与权系数的计算公式为:
〔3-16〕
为迭代次数,M为滤波器的阶数。
表示第时刻的输入信号矢量式中,式中,表示参考信号的信号矢量:
〔3-17〕
、分别表示第时刻的输出信号与输出误差,W(k)表示时刻权系数矢量:
〔3-18〕
表示LMS算法步长收敛因子。
自适应滤波器收敛的条件是:
〔3-19〕
其中是输入信号的自相关矩阵R的最大特征值。
的选取必须在收敛速度和失调之间取得较好的折中,既要具有较快的收敛速度,又要使稳态误差最小。
它控制了算法稳定性和自适应速度,如果很小,算法的自适应速度会很慢;如果很大,算法会变得不稳定。
由于LMS算法结构简单、计算量小、稳定性好,因此被广泛应用于系统辨识、信号增强、自适应波束形成、噪声消除以与控制领域等。
关于LMS算法的收敛速度,将讨论两点:
第一,对一个特定的信号环境,收敛速度和步长因子有何关系。
第二,信号环境本身的特性,对收敛速度有何影响。
从收敛速度的角度考虑,步长因子应该尽可能大,再看信号环境,即的特性对算法收敛性能的影响如果当特征值的分布围较大,即最大特征值和最小特征值之比拟大时,公比的取值幅度也将比拟大,算法的总的收敛速度将会变得比拟慢。
4.基于LMS的自适应滤波器应用仿真预测器