文档介绍:第7章自适应滤波器
LMS自适应滤波器的基本原理
Widrow-Hoff LMS算法
自适应滤波器的应用
我们已经研究了维纳滤波器与卡尔曼滤波器。维纳滤波参数是固定的,适用于平稳随机信号最佳滤波。卡尔曼滤波器参数可以是时变的,适用于非平稳随机信号。然而,只有对信号和噪声的统计特性先验已知条件下,这两种滤波器才能获得最优滤波。遗憾的是,在实际应用中,常常无法得到这些统计特性的先验知识。因此,这种情况下,用维纳滤波器和卡尔曼滤波器实现不了最优滤波。然而,自适应滤波能够提供卓越的滤波性能。
所谓的自适应滤波,就是利用前一时刻已获得滤波器参数等结果,自动地调节现时刻的滤波器参数,以适应信号或噪声未知的或随时间变化的统计特性,从而实现最优滤波。设计自适应滤波器时可以不必要求预先知道信号和噪声的自相关函数,而且在滤波过程中信号与噪声的自相关函数即使随时间作慢变化它也能自动适应,自动调节到满足最小均方差的要求。这些都是它突出优点,因此近年来它被广泛的应用于各种信号处理中。
LMS自适应滤波器的基本原理
如图7-1所示,其输入是一随机信号:
其中表示信号的真值, 表示噪声。其输出等于的估计值,用表示。维纳滤波器是具有这样的或的滤波器,它能使与间的均方误差
最小,即
从而达到最好地从噪声中提取信号的目的。而自适应滤波器则能自动调节它的值以满足上述最小均方误差的准则。
如果长为N,则从图7-1可以得到
(7-1)
这里, , 。
由此可见,输出是N个所有过去各输入的线性加权之和,其加权系数就是。在自适应滤波器中这个加权系数常用符号表示,所希望的输出常用表示,并为了书写简化,时间用下标表示,于是式(7-1)成为
(7-2)
由式(7-2)可见,自适应滤波
器可看成是自适应线性组合
器,如图7-2所示。
一般来讲, 其可以是任意一组输入信号,并不一定要求当时, ,…, ,即并不一定要求其各是由同一信号的不同延迟组成延时线抽头形式的所谓横向FIR结构,如图7-3(a)所示。
图7-3(b)是它的原理图。自适应滤波器的要害在于按照
及各值,通过某种算法寻找时的各,从而可自动地调节各值。
自适应滤波器可以用图7-4的简化符号表示。图7-5表示包括个自适应横向滤波器的自适应系统,当所处理的输入信号来自不同信号源时,它实际上等于自适应线性组合器。
利用讨论维纳滤波器时域解时的相同方法可以求得在
时的权系数。将式(7-2)写成矩阵形式有
(7-3)
这里
(7-4)
(7-5)
所以
(7-6)
令
(7-7)
=输入的自相关矩阵(7-8)
于是式(7-6)可以写成
(7-9)
注意,对于平稳输入,式(7-9)的是权矢量的二次方函数,因此~ 是一个凹的抛物体曲面,它具有唯一的极小点。可以用梯度方法沿该曲面调节权矢量的各元素,得到这个均方误差的最小点。
均方误差的梯度(用表示)可以通过将式(7-9),对权矢量的各进行微分得到
(7-10)
置就可得到最佳权矢量,用表示,即
(7-11)
式(7-11)是维纳-霍夫方程的矩阵形式。满足式(7-11)的
即为最佳权矢量或称维纳权矢量。将式(7-11)代回式(7-10),得最小均方误差为
(7-12)
实际上,上述这套方程与维纳滤波器推出的结果完全相同。自适应滤波器与维纳滤波器比较,其差别在于它加了一个识别控制环节,将输出与所希望的值比较,看是否一样。如果有误差,则用去控制
,使为的。因此它的关键在于怎样能简便地寻找,或者说用什么样的算法来求得,最常用的算法是所谓最小均方(Least Mean Square)算法,简称LMS算法。
Widrow-Hoff LMS算法
自适应滤波器的自适应过程的实际目的是要寻求,虽然按式(7-11) 可求得准确的。但需要预先知道相关矩阵P和R。当P和R不能预先获得时,就只能直接用数值计算的方法。在权的数目很大或者输入数据率很高时,这种方法将会遇到计算上的严重困难。这种方法不仅需要计算矩阵的逆,而且还需要测量或估算这么多个自相关和互相关函数才能得到P和R的矩阵元素。不仅如此,当输入信号的统计特性在慢慢地变化时,还必须从头重新计算。由于这些原因,人们宁愿应用另外一些更有使用价值的递推估计的算法。