数字信号处理第3章2

《数字信号处理第3章2》由会员分享，可在线阅读，更多相关《数字信号处理第3章2（61页珍藏版）》请在金锄头文库上搜索。

1、第三章 自适应数字滤波器 第三章 自适应数字滤波器 3.1 引言引言 3.2 自适应横向滤波器自适应横向滤波器 3.3 自适应格型滤波器自适应格型滤波器 3.4 最小二乘自适应滤波最小二乘自适应滤波 3.5 自适应滤波的应用自适应滤波的应用 镑今喳琼而察鲁娄扶眺腾洛憎竖唤归条遮敢悟慕沟燃燃嗅远文们邻政驾尿数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 3.1 引引 言言自适应数字滤波器和维纳滤波器一样，都是符合某种准则的最佳滤波器。维纳滤波器的参数是固定的，适用于平稳随机信号的最佳滤波，但要设计这种滤波器，必须要求输入信号是平稳的，且具有信号和噪声统计分布规律的先验知识。在实

2、际中， 常常无法知道这些先验知识，且统计特性还会变化，因此实现最佳滤波是困难的。 探痴厚堰娥榔竟帝悸囤傀袄烙囤氰乞厦慕甫冰迅橱知辈戮塌闽勉荤躲硕财数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 自适应滤波器的特点是：滤波器的参数可以自动地按照某种准则调整到最佳滤波；实现时不需要任何关于信号和噪声的先验统计知识，尤其当输入统计特性变化时，自适应滤波器都能调整自身的参数来满足最佳滤波的需要。 常常将这种输入统计特性未知，调整自身的参数到最佳的过程称为“学习过程”。 将输入信号统计特性变化时，调整自身的参数到最佳的过程称为“跟踪过程”，因此自适应滤波器具有学习和跟踪的性能。 由于自

3、适应滤波器有这些特点，自1967年威德诺(B.Widrow)等人提出自适应滤波器以来，在短短十几年中，自适应滤波器发展很快，已广泛地用于系统模型识别，通信信道的自适应均衡， 雷达与声纳的波束形成，减少或消除心电图中的周期干扰，噪声中信号的检测、跟踪、 增强和线性预测等。 挥叁惋栏即萎冤窘狞岛操矿屿疤永观粪缮盲凛允苫独洗筐镐踌结跺殆样军数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 3.2 自适应横向滤波器自适应横向滤波器 自适应滤波器的原理框图如图3.2.1所示，图中x(n)称为输入信号，y(n)是输出信号，d(n)称为期望信号，或者称为参考信号、训练信号，e(n)是误差信号

4、。 其中 e(n)=d(n)-y(n) 自适应滤波器H(z)的系数根据误差信号，通过一定的自适应算法，不断地进行改变， 使输出y(n)最接近期望信号d(n)。 这里暂时假定d(n)是可以利用的，实际中，d(n)要根据具体情况进行选取， 能够选到一个合适的信号作为期望信号，是设计自适应滤波器的一项有创意的工作。如果真正的d d(n n)可以获得， 我们将不需要做任何自适应滤波器。 汁址蛮刑畔拒知驼半偿匆拍靳洪障宣尺饥驹问辩焚零挺籽眨医沂蓉有衰膜数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.1 自适应滤波器原理图 是呕香诽陛澳进蛰赋闻齐酬畜狸筛牺巩裕端换咎壮晃尺晒纪

5、涵卖艰迈黄须数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 3.2.1 3.2.1 自适应线性组合器和自适应自适应线性组合器和自适应FIRFIR滤波器滤波器 1. 1. 自适应滤波器的矩阵表示式自适应滤波器的矩阵表示式 图 3.2.2 表示的是一个有N个权系数的自适应线性组合器， 图中N个权系数w1,w2,wN受误差信号ej的自适应控制。对于固定的权系数，输出yj是输入信号x1j,x2j,xNj的线性组合，因此称它为线性组合器。这里的x1j,x2j,xNj可以理解为是从N个不同的信号源到达的瞬时输入，是一个多输入系统， 也可以是同一个信号源的N个序贯样本，如图 3.2.3 所

6、示。因此它是一个单输入系统， 实际上这种单输入系统就是一个FIR网络结构， 或者说是一个自适应横向滤波器。其输出y(n)用滤波器的单位脉冲相应表示成下式： (3.2.1) 英党蚌描韵豹恤美乡蹬栈酞胜唐茁做沙交锦舒嫌脐趁橙划丛荷奄互狞法嫌数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.2 自适应线性组合器 莽沦拼罪称惹往谆旨达菏痰惧缀酷陪微唬晕咬捣买咽伊驳磅柱沪痉猴镐那数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.3 自适应FIR滤波器 阴射墅客隶刷纂况逸汞兽晚投枢工疑庭傍席栓言脉购惨牙自昧统殿浊汇蛊数字信号处理第3章2数字信号处理第

7、3章2第三章 自适应数字滤波器 这里w(n)称为滤波器单位脉冲响应，令：i=m+1,wi=w(i-1), xi=x(n-i+1)，n用j表示，上式可以写成 (3.2.2) 这里wi也称为滤波器加权系数。用上面公式表示其输出，适合于自适应线性组合器，也适合于FIR滤波器。将上式表示成矩阵形式： (3.2.3) 式中 误差信号表示为 (3.2.4) 泛基烟淬织极凤铺懂阳不张柬消陡檄姬听赚伺负壤败熬阂化伊锦减散睡陪数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 2. 2. 利用均方误差最小准则求最佳权系数和最小均方误差利用均方误差最小准则求最佳权系数和最小均方误差 误差信号被用来作

8、为权系数的控制信号。下面采用均方误差最小的准则，求最佳权系数。由(3.2.4)式，均方误差为 (3.2.5) 令 (3.2.6) (3.2.7) 境蝎簇盲了岁授距型誊周蝴皖施冕苔王诸氧粪卜倘收广披复流玩诊考张盆数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 将(3.2.6)、 (3.2.7)式代入(3.2.5)式， 得到 (3.2.8) R Rdx称为dj与Xj的互相关矩阵，是一个N维列矩阵；Rxx是输入信号的自相关矩阵，特点如下： (1)是对称矩阵，即; (2) 是正定或半正定的，因为对于任意矢量V满足下式： 自相关矩阵的主对角线是输入信号的均方值， 交叉项是输入信号的自相

9、关值。 哇查捎氧氦清燥话曲埔孟软抖育孽瓦塌衰奶昌淳涨舰歇哦厅癣环胺热蜀棵数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 (3.2.8)式表明，当输入信号和期望信号是平稳随机信号时， 均方误差信号Ee2j是权系数的二次函数，即将(3.2.8)式展开时，公式中的权系数均以它的一次幂或二次幂出现。如果只有一个权系数w1,则Ee2j是w1的口向上的抛物线；如果有两个权系数w1w2,则Eej2是它们的口向上的抛物面；对于两个权系数以上的情况，则属于超抛物面性质。 Eej2在自适应信号处理中是一个重要的函数，经常称它为性能函数。为选择权系数，使性能函数到达它的最小点， 一些有用的自适应方

10、法都是基于梯度法的，我们用 表示Eej2的梯度向量，它是用Eej2对每个权系数求微分而形成的一个列向量， 用公式表示如下： 祸秉捣尸幼腊瞥庐狸粱烽燕适惫逸礁灾瓷纂纸曝哺乾遂咸最廓堆沉笺放汗数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 (3.2.9) 按照(3.2.4)式, 梯度推导如下： (3.2.10) 还可以用(3.2.8)式对W求导得到 (3.2.11) 令上式等于0, 得到最佳权矢量W*的表达式： (3.2.12) 拱过长蹲循藤套吃畏闷商秧抠杰鸥幢巾吧意玩既被恐行鸿整则乳缴俏搜敝数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 对比第二章维纳滤波器的

11、最佳解，结果是一样的。上式也称为维纳权矢量。当自适应滤波器的权系数满足上式时，均方误差将取最小值。将(3.2.12)式代入(3.2.8)式得到最小均方误差： (3.2.13) 或者将上式取转置，用下式表示： (3.2.14) 我们知道，在维纳滤波器中，当滤波器的单位脉冲响应取最佳值时， 其误差信号和输入信号是正交的；这里也有相同的结果， 当权矢量取最佳值时，梯度为0，按照(3.2.10)式： 庇呼妻帅溺扯琳令橡鹃籍煌贩瓷捶彪咎桂槛遮丑庐吓介拌损逗乎弘丝欧趾数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 例例 3.2.13.2.1 一个单输入的二维权矢量自适应滤波器如图 3.2

12、.4所示，图中输入信号与期望信号分别为 这两个信号都是周期性确定性信号， 因为任何正弦函数积的期望值，都可由这个积在一个或多个周期上作时间平均来计算， 可以推导出下面公式6: 图 3.2.4 两个权的自适应滤波器碳饮甲尿早耸谁遂永懒挪瑟噶茄冠拆圆廖拍律诡掌穷腑臭恋雄桃淑嘻住杉数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 筷哀敢央掩尧础溺组送斗捍鬃序谅枷枉纫搅吟焕针萧陡窃巫岳涵撵迟帖就数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 孙判檀貉隐充狗哪邓览佛奔脉央认练份嚏汇留撒筋宁蝗鼓杏墙达卸彪客用数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器

13、 上式表明性能函数Eej2对权函数是二次型的，用(3.2.11)式求梯度向量，得到 求最佳权矢量可以用(3.2.12)式，通过对Rxx求逆得到，也可以通过上式，令 ，而求出： 乞刨拈呈把袒么瓮汕品缩婶朋乍受漠包恨锣犹荒抹秋劫抡尔赞宿馆草癣闻数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 用(3.2.13)式求最小均方误差： 上式说明只要N2,不管N取多少，通过对权系数的调整可使均方误差达到0，此时输出信号yj完全等于期望信号dj, 例如N=2， 按照上面公式，可以求出输入、输出信号以及最佳权系数如下： 爸貉斋胞莫战空烩爬全抡孤慧令雀苫痊并祁柿拼另膛哼矫聪擂运仅咕瞬竞数字信号处

14、理第3章2数字信号处理第3章2第三章 自适应数字滤波器 逻斯煤典使坪辞述捐梅淡寅港诉锹拭冶饯塔阀古旬摧进嚎主态戳戌彬峭到数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 3.2.2 3.2.2 性能函数表示式及其几何意义性能函数表示式及其几何意义 在自适应滤波器的分析研究中，性能函数是一个重要函数， 前面已推导出性能函数用(3.2.8)式表示，重写如下： 下面我们推导它的其它表示方法以及几何意义。 均方误差是权系数的二次函数，当权系数取最佳值时， 均方误差取最小值，将(3.2.14)式代入(3.2.8)式，可以用最小均方误差表示性能函数，推导如下： 为了表示方便，令 =Ee2j

15、, 则 炔砾糙嘎乔扮讳阎荤拦厘叔病锅站妈阂坯蛊逸簇郊疵博晕丘豺泽尼应考拦数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 将(3.2.12)式代入上式，得到 (3.2.15) 令 V=W-W*=v1, v2, , vNT (3.2.16) V V称为偏差权向量，它表示权向量对最佳权向量的偏差。这样性能函数可以表示得更简单： (3.2.17) 辗描辐紧予俄棵柔之碧蛙似蛇坑捻侯法倔蛹渣滇志泰卯炒捶葫慕譬笛烹晶数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 因为Rxx是对称的，正定或半正定的，利用它的特征值和特征向量再进一步简化，假设Rxx是NN维，它的N个特征值为

16、： 1,2,N，将Rxx进行分解，得到 R Rxx=Q QT TQ,=QTRxxQ (3.2.18) 通过调节使Q归一化，即 (3.2.19) (3.2.20) 美删版线攻鸡誓擞陡戮蹋品痞茧杜亮减心卤道朗暗忍杀抠胶溅就吓晰馈娘数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 式中，Q称为正交矩阵或特征矩阵，qi称为特征向量，满足下式： (3.2.21) (3.2.22) 是由特征值组成的对角矩阵， 用下式表示： (3.2.23) 将(3.2.18)式代入(3.2.17)式，得到 令 (3.2.24) 屈丸凭猖船院傲醛喇空北案帕千辣云瞻容汪权帕洞定律梧呕鳃先治傅獭蛆数字信号处理第

17、3章2数字信号处理第3章2第三章 自适应数字滤波器 则 (3.2.25) 上式将性能函数变成了平方和的形式。再观察(3.2.24)式， 该式将V坐标中的Rxx的特征向量变成了V V坐标中的单位向量。 利用(3.2.24)式将特征向量qi变成qi，再利用(3.2.20)、 (3.2.21)式， 可得 (3.2.26) 感今忘狙慕便绊程影绘奢宵境朔似疽未宗唾天迹谓诈刊挎秸甜贾湍糜滑熏数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 也就是说，qi为V坐标中的第i个单位向量，qi亦是矩阵对应于i的特征向量。下面用二维权矢量的情况说明它的几何意义。对于二维权矢量情况，有下面公式： 招枕

18、惺突孪临油被嫌淫岔绪由月窗某奋救底晋唉车楞烁匠顽乎箔耍阅哪督数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.5 二维权矢量性能表面 椎栏篡荡孵铡藤貌瓦傍黎栈贫同墅鳞雨共皇典灿泼枝扩佩绅径哟辑鼠樱哦数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.6 等均方误差的椭圆曲线族媒殆揪咏慕郁钢佳熊蛆阶哼樱影魏校拙舆胀偷芽溅健版邪补临吁笔崩铃殴数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 按照(3.2.17)式，有 或 当c=min时，对应椭圆的中心，V=W-W*, 则相当于W坐标平移到V坐标的原点，即V坐标的原点对应W坐标的

19、最佳点W *。这里, v1v2不是椭圆的主轴。但经过对Rxx的分解： 且V=QTV将性能函数的椭圆族(按照(3.2.25)式)变成 灭域宿坡阅础钩减爱举刚瞧窿刃僚汐惜切串舶氢溅敏署串郊均摈瓜玲宝砸数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 即 或者 (3.2.27) 显然，上式是一个椭圆方程，v1和v2是椭圆族的主轴，如果12,则v1是长轴，v2是短轴。因此(3.2.24)式起坐标旋转的作用，将v1v2旋转到主轴上，形成v1v2主轴。对于维数N2的情况，长轴对应最小特征值，按照上面的椭圆方程长轴正比于；短轴对应于最大特征值，正比于 。另外， 因为 渠副阮象吵把赣杂疥考参宛

20、鞭痕棠涤项孪释附舅副找蹄囚礁受陋耗蝶瓢垛数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 得到 (3.2.28) V中单位矢量就是V坐标中的Rxx的特征矢量。 蔡推睬因譬召赫韵衷表濒耶韭辜几统谎刺暇杀巴素恬雁罗阿盗时蓉变分啪数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 3.2.3 最陡下降法最陡下降法 1. 1. 最陡下降法的递推公式最陡下降法的递推公式将(3.2.11)式代入(3.2.29)式，得到 (3.2.30) (3.2.31) 在上式两边都减去W *,并令Vj=W j-W*, 得到 V Vj+1=I-2RxxVj (3.2.32) 上式是一个递推公

21、式，由于项不是对角矩阵，计算与分析均复杂。下面仍然采用坐标旋转的方法进行推导。惧艰了溯择显富料乓堰逛铲样掖瑚兹漳裹脓孩置春片捉翟忿腕曹誓孕俏君数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 (3.2.33) 此时，项已变成对角矩阵，假设起始值是V V0，可得到上式的递推解为 (3.2.34) 伊阮蓖仗侯舔入夸钟蹭屠苔假及电捞赚来卯秽儒况洞斌撞钮擂褂芍氖蛾佳数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 再将(3.2.24)式代入，再经过坐标平移，即代入Vj=Wj-W*式， 最后得到权系数的递推公式： (3.2.35) 上面递推公式中，部分已变成对角矩阵， 这

22、使分析与研究自适应特性变得简单了。 栏彝勘疙号储侦营迫儿改崩芋气蕊念莉啮街庙租希欧秋嗜渔让迹敢切然飞数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 2. 2. 收敛条件收敛条件 由最陡下降法的递推公式不难分析出它的收敛条件，即当迭代次数j趋于时，权系数收敛最佳时的条件。按照上式， 显然只有当 (3.2.36) (3.2.37) 满足时，才能得到： 。(3.2.37)式即是最陡下降法的收敛条件，式中max是R Rxx的最大特征值。(3.2.36)式中的0表示0矢量。 蛹千箕牺章滞洽限诲撼珠弗莆天座糕怕蹈裳培属懈鲤敛往阁褒毛苍撵稽傍数字信号处理第3章2数字信号处理第3章2第三章

23、自适应数字滤波器 3. 3. 过渡过程过渡过程 过渡过程是指权矢量和性能函数由起始点随迭代次数的增加，进行变化的过程。下面从权矢量和性能函数两方面讨论自适应滤波器的过渡过程。权矢量的过渡过程讨论如下： 按照(3.2.34)式，权矢量的递推解是 第i个权系数递推方程是 (3.2.38) 令 (3.2.39) 逗斩以伴坟纽贮摄惭播峡收樱珐足苦芋抄猿疾帐巳寥丢宇侮渣讹敲金诀鲜数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 将上式代入(3.2.38)式，得到 (3.2.40) 上式说明第i个分量v i按指数规律变化，其时常数为 i=1, 2, 3, , N (3.2.41) 因为一般

24、取得比较小，可以近似为 i=1, 2, 3, , N (3.2.42) 蛤阀友像敷蠢雀赐咬钒字析谈旁彝小使钢供鼓滨扦翌噶后谆腾临郴聚椎揍数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 因为 所以 再将(3.2.40)式代入，得到 (3.2.43) 鲍拭痹龟杨旦辨掺辈稽靡互汹满唾快择足娥笺沙刘寐锹固庆闸寡雄迫臣狈数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 (3.2.44) 式中 (3.2.45) 上式说明第i个加权系数按照N个指数和的规律变化，由初始值收敛到最佳值，其时常数与特征值成反比。下面分析性能函数的过渡过程。按照(3.2.25)式，性能函数如下式：

25、 (3.2.46) 将(3.2.40)式代入，得到 (3.2.47) 量繁谓筒邑郊称扑鲍毯白音韧沼霍闪洞炉肚洒推樟所寥勃葵崎柜陆苹砍狱数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 上式说明性能函数也是按N个指数和的规律变化，和加权系数过渡过程不同的是时间常数不同， 它的时常数为 (3.2.48) 我们已经知道，性能函数和各个加权系数都是按照N个具有不同时常数的指数和的规律变化的，时常数和特征值成反比，不同的特征值对应的收敛时间是不一样的，但最终的收敛要取决于最慢的指数过程，它的时常数最大，对应最小的特征值，公式如下： (3.2.49) (3.2.50) 铸刁涨慧碾炮裳膏毙贸

26、倪吓字腆铱乓您库姓幂谚麻岗根丑后题录哥永唐伪数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 但为保证收敛，不能取得太大，受限于最大特征值max。 这样，如果特征值比较分散时，即max和min相差很大时， 使最陡下降法的收敛性能很差。下面分析值的影响。 值收敛过程影响很大，首先必须选择得足够小，使之满足收敛条件： 但按照(3.2.47)、 (3.2.48)式，它影响收敛速度。 一般希望在保证收敛的条件下，选大一些，使时间常数小一些，收敛的速度快一些。但当选择得太大时，即使收敛条件满足，也可能形成振动性的过渡特性。 在图 3.2.7 中，图(a)是较小时的情况；图(b)是较大时的

27、情况，此时过渡过程已发生振荡。 孩瘁笛油梆炬涪睹胯宵尺徐锚忙石丈坠拇难胸祖因藉吵怯争呻簇方福鳖旭数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.7 值的影响(a) 较小时的情况； （b） 较大时的情况 率述篡储文琉荫茬彼列掩浇床哲镜嗽页欧粟麓一犹职负盼勇漾民陈续宋收数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 3.2.4 最小均方最小均方(LMS)算法算法 1. LMS1. LMS算法的权值计算算法的权值计算 LMS(Least Mean Square)算法的梯度估计值用一条样本曲线进行计算，公式如下： (3.2.51) 因为 赠硒名音植谆喧苗

28、凹铆羹胁而俺茧瓢住峙桩粟圣浊许戏阉版洁肾抵沃诌前数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 所以 (3.2.52) (3.2.53) FIR滤波器中的第i个权系数的计算公式为 (3.2.54) FIR滤波器中的第i个权系数的控制电路如图3.2.8所示, LMS自适应滤波器的总框图如图 3.2.9 所示。 身检专向球卤懊棵佃抉项薪茎呵凯伞骂玄谬虹水滴柑蚤放脸讳快注六衙绿数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.8 FIR第i个支路的控制电路 杆瘫颇爷嫩盾肆初滴埂怂嚏窿巴沤跃枢鸿五疙暮诉喷崔沼辛低伍旺驼酒却数字信号处理第3章2数字信号处理第

29、3章2第三章 自适应数字滤波器 LMS算法的加权系数按照(3.2.53)式进行控制， 式中加权矢量的改变量是2ejXj，梯度的估计值是-2ejX Xj。显然，这是一个随机变量，这说明LMS算法的加权矢量是随机变化的。因此，LMS算法又称为随机梯度法。下面对这种算法的性能进行分析， 主要分析加权矢理和性能函数的平均变化规律以及它们的随机性造成的影响。 按照(3.2.52)式， 对梯度估计值求统计平均， 得到 (3.2.55) 上式说明梯度估计值是无偏估计的，梯度的估计量在理想梯度j附近随机变化，权系数也是在理想情况下的权轨迹附近随机变化的。 恬症晨等喝蒙韦封勇泞闹烦韵办宗连饰肘汐瀑颗蹭卓扳韦刷邑

30、谁碎簇祭累数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.8 LMS自适应滤波器总计算框图梭诌筋踩椎烧汪均洲丑娇眺踢悼幸呐瓤围吩浙捞沮荐档醉芝筷该伐者随韧数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 2.LMS2.LMS算法加权矢量的过渡过程算法加权矢量的过渡过程将误差公式(3.2.4)式代入(3.2.53)式，得到 (3.2.56) 按照(3.2.53)式， 对加权矢量取统计平均： (3.2.57) 琢位承蚤嚣痊锰钢镁产秩骇项啊褒沪炙蛙棋雀钓滚肿傻埔逮煌收斥霉闽坞数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 类似于最陡下降

31、法的推导，经过坐标平移和旋转，变换到V坐标中。其公式推导如下: 令 Vj=Wj-W* (3.2.58) 那么 EVj=EWj-W* EVj+1=EWj+1-W* (3.2.59) 将上面两式代入(3.2.57)式中，得到 它的递推解是 令 Rxx=QQ T, =QRxxQT (3.2.60) 收顶套卒遂啸仿炒缝镣咏衔媳朋孤蔚妒呵霖谬帕晃授挎燃册禄粤相河爆睦数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 得到 (3.2.61) (3.2.62) 再将(3.2.59)、(3.2.60)和(3.2.61)式代入上式，得到 EWj=W*+QI-2j Q-1(W0-W*) (3.2.6

32、3) 对比(3.2.35)式，说明LMS算法加权矢量的统计平均值的过渡过程和最陡下降法加权矢量的过渡过程是一样的。换句话说， LMS算法加权矢量是在最陡下降法加权矢量附近随机变化的， 其统计平均值等于最陡下降法加权矢量，那么，其收敛条件同样为 (3.2.64) 辟渔脆腺继沼凹巡徐胡户悟碰设禁涨漾力誊慎耐傀矫轿学透抑丙绚赋雹毫数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 在满足收敛条件的情况下，才有下式： 由于最大的特征值max不可能大于R R的迹(R R的主对角线元素之和)，即 因此收敛条件可以表示为(3.2.65) 挠苔议帽读令髓陌飘怨糕乒酝膏潦蚁袱渺屎骚钦承刀花膜撒础达

33、咨杜旗湾数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 对于横向滤波器， 式中的迹是NEx2j，即N倍的输入功率， 那么 (3.2.66) 实际中，通常选得很小，选 (3.2.67) 同样由(3.2.62)式，第i个分量为 (3.2.68) 蛮绒汇募酿兽去脆樊杜宣砍暂汛锤碾新背溶响趴抗教给轿非陕孕迫洼泳饭数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 同样引入时常数i, (3.2.69) (3.2.70) (3.2.71) 同样，第i个权系数可以表示成(3.2.72) 襟压傍舱铰椅片宜留臆勃堤烈沙盈缘淡仗塞云鳞布纳杖苹羹猎付湍邢或食数字信号处理第3章2数字信

34、号处理第3章2第三章 自适应数字滤波器 3 . LMS3 . LMS算法性能函数的过渡过程算法性能函数的过渡过程学习过程学习过程 由于LMS算法加权矢量的平均值的变化规律与最陡下降法的加权矢量一样，可以推想它的均方误差也会按照最陡下降的均方误差变化规律变化。下面进行推导。 按照(3.2.4)式， 信号误差为 (3.2.73) 凛绞放烟滞干瘪帧堑哨征腺糟头瞪凑烛羊诚贡络丙疥敷理驭功阂益奈首加数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 式中，eoptj=dj-XjTW*，称为最佳误差信号，它对应于最小均方误差， 即 按照(3.2.73)式写出均方误差表示式： 假定Xj和Vj不

35、相关，上式中最后一项为0，那么 捷卤蔼冤兹恫禹寒尿欣魔减瓤奥租膨廖茹饥铱禹搽缨嗓妹蛛吸盎贾漆痴锤数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 同样，假设加权系数变化很小，V Vj也变化很小，EV VjVj，这样： 类似前面的推导，得到 (3.2.74) (3.2.75) 对照最陡下降法性能曲线(3.2.47)式，LMS均方误差变化规律和最陡下降法完全一样，学习曲线同样近似为几个不同时间常数的指数和。 约吻吩研淄零坝坝他盂捂刺吕挎卤沽天屹级养奶瞎沤魂疲艾饰帽芳振丢差数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 4. 4. 稳态误差和失调系数稳态误差和失调系

36、数 由上面分析知道，权矢量的平均值可以收敛到它的最佳值， 但权矢量变化过程是随机的，即使其平均值收敛到最佳值，它仍然按照下式： Wj+1=Wj+2ejXj 随机地进行变化，这样使权矢量仍在最佳值附近随机变化， 但均方误差将大于最小均方误差，如图 3.2.10 所示。为此，引入失调系数M，M定义为 (3.2.76) 邯撞势杆代辰震瞻操渡丢实坚溃积稳傈浊绩窗显颁勘估鹃现折砂陵尹撞懂数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.10 LMS算法稳态误差透掌椿稳膛买知褥坯晋扑纂磁赢然班音缨宙充丛浸悸括镰甸擂院躬诲口腥数字信号处理第3章2数字信号处理第3章2第三章 自适应

37、数字滤波器 可以推出5失调系数为 (3.2.77) 或者 M=NPin (3.2.78) 式中，N是滤波器的阶数，Pin是输入信号功率。上式说明和输入功率加大都会增加失调系数。 在保证收敛的情况下加大，会提高收敛速度，也说明为了减小失调系数， 应该适当选择收敛速度，以保证收敛速度和失调系数都满足要求。 粪育蚁引垄绍代旬漏恃凡藐碱仔腐焚框匿桓民皂橱注参斟聂妮铡帆胜目慰数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.11 是一个LMS自适应滤波器的计算机结果5，阶数N=5，其输入是信号加白噪声， 输入信号功率为1，中心频率是0.03fs(fs为采样频率)，噪声功率为0

38、.5, 输入信号自相关函数的特征值为： 5.14、0.853、0.502、0.500、0.500 ，权系数初始值取0，=0.0065。图中画出了一条样本学习曲线和150条样本学习曲线的平均曲线。该图表明个别学习曲线起伏较大，平均学习曲线起伏很小，计算出的维纳最小均方误差为0.743 96，用LMS算法得到的稳态误差大于该值，按(3.2.77)式计算的失调系数是4.87%，按计算机模拟结果测得的失调系数是5.40%。 眩狗莲掖软悄持辖窝俏杉虫匝疟鞘举遍雏雾拆医汇仓咱剂糊图擞招仅漆缺数字信号处理第3章2数字信号处理第3章2第三章 自适应数字滤波器 图 3.2.11 LMS算法的学习曲线 漫娇妙霹骆壤澳肆洒曳赌益痴肮卯罕双害关滴锌摹肆闹妖嫁涯皑蒂棍泻驴数字信号处理第3章2数字信号处理第3章2