自适应滤波器课件.ppt
《自适应滤波器课件.ppt》由会员分享,可在线阅读,更多相关《自适应滤波器课件.ppt(94页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、自适应滤波器第1页,此课件共94页哦3.10 自适应的递归最小二乘方算法3.11 IIR递推结构自适应滤波器的LMS算法3.12 自适应滤波器计算举例3.13 自适应滤波器的数字实现3.14 最小二乘自适应滤波器3.15 最小二乘格形自适应算法3.16 快速横向滤波自适应算法3.17 自适应滤波器的应用第2页,此课件共94页哦设计维纳和卡尔滤波器,要求已知关于信号和噪声统计特性的先验知识。但在许多情况下人们对此并不知道或知道甚少,某些情况下这些统计特性还是时变的。处理上述这类信号需要采用自适应滤波器。自适应信号处理器分为两大类,一类是自适应天线,另一类是自适应滤波器。本章主要讨论自适应滤波器的
2、工作原理、基本原理、重要算法和典型应用。第3页,此课件共94页哦3.1 自适应滤波原理自适应滤波器由参数可调的数字滤波器和自适应算法两部分组成,如下图:d(n)x(n)e(n)y(n)-+参数可调数字Filter自适应算法与参考信号比较经过自适应算法对Filter参数进行调整。自适应算法的原则:最终使e(n)均方值最小!自适应滤波器是一种能自动调整本身参数的特殊维纳滤波器。它在设计时,不需先知道输入信号和噪声的统计特性。它能在自己工作中逐渐学会or估计出所需的统计特性。并以此依据自动调整自己的参数以达到最佳滤波的目的。第4页,此课件共94页哦自适应滤波器常见的例子:1.自适应预测:(可用于语音
3、编码、谱估计、信号白化等)输入信号是s(n),输出响应是预测值期望响应d(n)是n+D时刻的信号值第5页,此课件共94页哦2.自适应建模(a)(b)(a)是正向建模,(b)是逆向建模。在正向建模中,自适应处理器调整自己的权值,使得输出响应y(n)尽可能逼近未知系统的输出d(n)。在逆向建模中,自适应处理器调整自己的权值以成为被建模系统的逆系统,即把被建模系统的输出转换成为输入信号的延时。第6页,此课件共94页哦3.自适应干扰器传感器阵列接收到目标信号,导向延时使其预定观测方向上波束增益最大。固定目标信号滤波器输出为自适应处理器输出是噪声的估计 ,并用来抵消常用于波束形成器。第7页,此课件共94
4、页哦设计自适应滤波器时,首先要确定滤波器的结构(FIR,IIR或格形结构),然后设计自适应算法以调整滤波器参数,其目标是使某一特定的代价函数最小化。(本章选择均方误差为代价函数)第8页,此课件共94页哦3.2 自适应线性组合器自适应线性组合器是一种参数可自适应调整的有限冲激响应数字滤波器,具有非递归结构形式,分析实现简单。在大多数自适应信号处理中得到广泛应用。自适应线性组合器的一般形式:第9页,此课件共94页哦输入信号矢量:输入信号和输出信号之间的关系式为:对单输入情况:对多输入情况:还可表示为:自适应线性组合器按照误差信号均方值(or平均功率)最小的准则(即)来自动调整权矢量,选择什么信号作
5、为参考响应,要根据不同的应用要求来确定。第10页,此课件共94页哦3.3 均方误差性能曲面 由上面三式,得均方误差表示式:将式进一步写成:均方误差 是权矢量W的各分量的二次函数,即若将该式展开,则W各分量只有一次和二次项存在,的图形一定是L+1维空间向中一个中间下凹的超抛物面,有唯一最低点 min,该曲面称为均方误差性能曲面。第11页,此课件共94页哦当输入w(n)只有两个元素时,可得到如图的自适应滤波器:自适应过程,即自动调整权系数w(n),使均方误差达到最小值的过程,相当于沿性能曲面往下搜索达最低点的过程。最常用的搜索方法是梯度法,因此,性能曲面的梯度是一个很重要的概念。第12页,此课件共
6、94页哦均方误差性能曲面的梯度为:,让 得此式表明:当W偏离最佳值W*一个数值 ,将比 大一个数值第13页,此课件共94页哦3.4 二次性能曲面的基本性质平稳随机信号的统计特性是不随时间变化的。因此,其性能曲面在坐标系中是固定不变或“刚性”的。自适应过程就是从性能曲面上某点(初始状态)开始,沿着曲面向下搜索最低点的过程。但对非平稳随机信号来说,这种性能曲面是“晃动的”、“模糊的”自适应过程,不仅要求沿性能曲面向下搜索最低点,而且还对最低点进行跟踪。我们这里只讨论平稳随机过程,且为方便理解,只讨论两个权系数W0和W1的自适应线性组合。此时性能曲面是三维空间 中的一个抛物面。第14页,此课件共94
7、页哦现用一个与 平面平行与其相距 的平面切割该抛物面,交线在 平面上投影是一个椭圆。如图:椭圆中心为 ,它是性能曲面最低点 的投影。如果用若干个与 平面距离不同的平行平面来切割性能曲面,交线投影将是一组中心同在W*的椭圆。它们各与一个确定的 相对应。因此称为等均方误差线or等高线。第15页,此课件共94页哦等高线方程:由 得:常数若将坚持原点平移至 ,得到权偏移矢量全标系 等高线方程:常数 (可由 得到)这是一组同心椭圆,中心位于新坐标原点V=0。将上面讨论推广到L+1个权系数的情况不难想象,等高线将是L+1维空间中的一组同心超椭圆,椭圆中主位于坐标系 的原点。这组同心超椭圆有L+1个主轴,它
8、们也是均方误差 曲面的主轴。第16页,此课件共94页哦是R的特征值矩阵:可由R的特征方程detR-aI=0解出。第17页,此课件共94页哦最终得:由此总结出二次性能曲面的三个基本性质:主轴是R的特征矢量。(1)输入信号自相关矩阵R的特征矢量 确定了性能曲面的主轴;(2)因此它定义的旋转系统 就是椭圆的主轴系统。(3)R的特征值给出了性能曲面沿主轴的二阶导数值。第18页,此课件共94页哦3.5 最陡下降法前面分析知,自适应线性组合器的均方误差性能曲面是权系数的二次函数,但在实际应用中,性能曲面的参数甚至解析式都是未知的。因此,只能由已测数据,采用某种算法对性能曲面自动进行搜索。寻找最低点,从而得
9、到最佳权矢量。牛顿法和最陡下降法是两种著名的方法,牛顿法在数学上有重要意义,但实现很困难。因此,我们只介绍最陡下降法,它在工程上易于实现。最陡下降法是沿性能曲面最陡方向向下搜索曲面最低点。曲面的最陡下降是曲面的负梯度方向。这是一个迭代搜索过程。第19页,此课件共94页哦最陡下降法迭代计算权矢量公式为:是控制搜索步长的参数-称为自适应增益常数或称为收敛因子(常数)。将梯度公式代入上式,得:方程由 计算很困难,一般将w坐标通过平移 坐标,通过旋转到主坐标第20页,此课件共94页哦即:由于它们之间没有耦合,所以可分别由初始权值进行迭代运算求解,可得:第21页,此课件共94页哦即:为确保算法收敛,有
10、,即收敛到V的原点,即W的W*点。因此必须保证可以由给定的 求由R的特征值 在此范围内选取!这样计算仍比较繁锁,可采用直接估计 的方法,让R矩阵的迹:第22页,此课件共94页哦 也可由输入信号取样值进行估计:由于实际自适应F中调整参数是 可将上面结果返回到自然坐标系去,以看清W(n)的自适应调整规律。由 有:利用恒等式 有:从而取第23页,此课件共94页哦3.6学习曲线和收敛速度在自适应调整权系数的过程中,均方误差是迭代n次数的函数,称为学习曲线。这就是最陡下降法学习曲线的表达式,收敛条件约为:第24页,此课件共94页哦收敛速度的快慢,用时间常数来说明:有三个常用时间常数 权系数衰减时间常数:
11、第一就是 ,权系数衰减时间常数:定义:即 衰减为 的 倍时,所经历的迭代次数即为 。通常 其中第25页,此课件共94页哦 学习曲线时间常数其物理意义:衰减为 -的 倍所需迭代次数。3.自适应时间常数 (Samples样本数or取样周期)即:将 迭代次数用其取样间隔来度量。第26页,此课件共94页哦3.7 自适应的最小均方算法最陡下降法,每次迭代都需要知道性能曲面上某点的梯度值,而实际上梯度值只能根据观测数据进行估计。LMS算法是一种很有用的估计梯度的方法。它的突出优点是计算量小,且不脱线计算,只要知道输入信号和参考响应。LMS算法核心思想是用平方误差代替均方误差。LMS算法的基本关系式第27页
12、,此课件共94页哦 LMS算法调整权系数时不需要进行平方运算和统计平均运算,因而实现起来很简单。下图为两个权系数的自适应线性组合器采用LMS算法的计算流程图:d(n)e(n)y(n)+2第28页,此课件共94页哦第29页,此课件共94页哦而最陡下降法梯度公式:说明:LMS算法得到的权矢量的期望值与最陡下降法权矢量本身一样,因此,当收敛条件 或 满足时,随迭代次数趋近于无穷,LMS权矢量的期望值将趋近于最佳权矢量。第30页,此课件共94页哦3.8 权矢量噪声LMS算法之所以简单,主要是因为它对梯度矢量各分量的估计是根据单个数据取样值得到,没有进行平均。也正因此相当于使梯度估计中存在噪声。梯度估计
13、的噪声矢量用N(n)表示,有:协方差:均值为0,协方差=自相关函数,又 不相关第31页,此课件共94页哦将上式变换到主轴坐标系:这是计算噪声的协方差公式。下面考虑自适应调整权矢量这程中,梯度估计噪声的影响,即LMS对权矢量的影响。平移坐标系V:平移坐标系V:由噪声产生第32页,此课件共94页哦用归纳法求上式,即得:另推:这就是LMS算法中梯度估计噪声在稳态权矢量中引起的噪声。稳态解第33页,此课件共94页哦3.9 失调量梯度估计噪声的存在,使得收敛后的稳态权矢量在最佳权矢量附近随机起伏,这意味着稳态均方误差值总是大于最小均方误差 ,且在 附近随机地改变。如下图:将这种偏移量的期望值称为:超量均
14、方误差,用“超量MSE”表示。第34页,此课件共94页哦变换到主轴坐标系,超量 超量将 代入,可得近似公式:超量均方误差是度量这种性能损失的一个量(权系数在最佳解附近随机变动)。另一个度量自适应性能损失的量是失调量。用M表示超量第35页,此课件共94页哦 失调量和收敛速度这两个量要折衷加以考虑,学习曲线时间常数为:根据上式可将R的迹写成:得到:上式说明了失调量、学习曲线时间常数以及权系数的个数三者之间的关系式。在特征值未知的情况下,这个近似式对于设计自适应系统很有用的。在所有特征值相等时:第36页,此课件共94页哦3.10 自适应的递归最小二乘方(RLS)算法1.特征:FIR维纳滤波器的一种时
15、间递归算法 严格以最小二乘方准则为依据的算法 收敛速度快 应用于:快速信道均衡 实时系统辨识 时间序列分析 缺点:每次迭代计算量大,N阶的FIR计算量 曾经一度受到限制,近来有大的发展。2.RLS算法:用二乘方的时间平均的最小化准则取代最 小均方准则,并按时间进行计算。第37页,此课件共94页哦具体说:是要对初始时刻到当前时刻所有误差的平方进行平均并使其最小化,再按照这一准则确定FIR滤波器的权系数矢量W。即其中:期望响应 是N阶FIR的输出响应对于非平稳输入信号,为了能很好的进行跟踪,常引入一个指数加权因子对式进行修正:第38页,此课件共94页哦 遗忘因子(的正数)旧数据的权值按指数规律衰减
16、,越旧对 影响越小,忘性越大。求偏导 这是最小二乘方准则所对应的正交方程。上式经整理后得:或 定义:第39页,此课件共94页哦则标准方程可写为简化形式:该方程的解为:都与时间n有关,所以w是n的函数3.迭代算法将两式写成迭代形式第40页,此课件共94页哦同理:由利用矩阵恒等式则上式可写成:第41页,此课件共94页哦将上两式代入并利用上式得到式中增益误差上式具有卡尔曼滤波器的形式。是根据 及以前所有数据得到的最佳权值,根据它来预测 应该是合理的。可将增益公式写成第42页,此课件共94页哦 上式与LMS算法的差别在于权矢量校正项中出现了因子 。由于R(n)是自相关矩阵 的一种量度,因此,因子 的出
17、现使得RLS算法具有快速收敛的性质。RLS算法步骤如下:(1)在时刻n,已获得w(n-1)、T(n-1)和d(n),x(n)也已存储在滤波器的延时部件中。(2)计算出T(n),w(n),k(n)和e(n|n-1),并得到滤波器输出响应y(n)和误差e(n),即第43页,此课件共94页哦(3)进入第n+1次迭代遗忘因子 的选取 的值对算法影响最大。算法的有效记忆长度用 来度量,定义为 越小对应的 越小,意味着对信号的非平稳跟踪性能越好。太小,会小于信号的每个平稳的有效时间,不能充分利用所有已经获取的取样数据。结果,算出的 将会受到严重噪声影响。对平稳信号,最佳值为1。例子:第44页,此课件共94
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 自适应 滤波器 课件
限制150内