【文章內(nèi)容簡(jiǎn)介】 ????Njjjjj wewewee222122? ?() 因?yàn)? jjj XWde T??第三章 自適應(yīng)數(shù)字濾波器 所以 jNjjj Xwewewe?????????????? T21, ?jjj Xe2? ???() () jjjj XeWW ?21 ???FIR濾波器中的第 i個(gè)權(quán)系數(shù)的計(jì)算公式為 Nixeww ijjijij ,3,2,12 ,1 ?????? ?() FIR濾波器中的第 i個(gè)權(quán)系數(shù)的控制電路如圖 , LMS自適應(yīng)濾波器的總框圖如圖 所示 。 第三章 自適應(yīng)數(shù)字濾波器 圖 FIR第 i個(gè)支路的控制電路 z － 1wi( n ＋ 1)＋＋xi( n ) wi( n )2 ?e ( n )控制電路第三章 自適應(yīng)數(shù)字濾波器 LMS算法的加權(quán)系數(shù)按照 ()式進(jìn)行控制 ， 式中加權(quán)矢量的改變量是 2μ ejXj， 梯度的估計(jì)值是 2ejXj。 顯然 ， 這是一個(gè)隨機(jī)變量 ， 這說明 LMS算法的加權(quán)矢量是隨機(jī)變化的 。 因此 ，LMS算法又稱為隨機(jī)梯度法 。 下面對(duì)這種算法的性能進(jìn)行分析 ， 主要分析加權(quán)矢理和性能函數(shù)的平均變化規(guī)律以及它們的隨機(jī)性造成的影響 。 按照 ()式 ， 對(duì)梯度估計(jì)值求統(tǒng)計(jì)平均 ， 得到 jjjj XeEE ????? ][2]?[() 上式說明梯度估計(jì)值是無偏估計(jì)的 ， 梯度的估計(jì)量在理想梯度▽ j附近隨機(jī)變化 ， 權(quán)系數(shù)也是在理想情況下的權(quán)軌跡附近隨機(jī)變化的 。 第三章 自適應(yīng)數(shù)字濾波器 …控制 1控制 2控制 N…x1( n )x2( n )xN( n )w1w2wNy ( n )d ( n )e ( n )＋－…圖 LMS自適應(yīng)濾波器總計(jì)算框圖 第三章 自適應(yīng)數(shù)字濾波器 將誤差公式 ()式代入 ()式，得到 jjjTjjjTjjjjjjdXWXXiWXXdXWW???2]2[][21??????? () 按照 ()式， 對(duì)加權(quán)矢量取統(tǒng)計(jì)平均： *T12][)2(])[(2][][][(2][][2][][WRWERIWERRWEWXXEXdEWEXeEWEWExxjxxjxxdxjjjjjjjjjjj?????????????????() 第三章 自適應(yīng)數(shù)字濾波器 類似于最陡下降法的推導(dǎo) ， 經(jīng)過坐標(biāo)平移和旋轉(zhuǎn) ， 變換到V′ 坐標(biāo)中 。 其公式推導(dǎo)如下 : 令 Vj=WjW* () 那么 E［ Vj］ =E［ Wj］ W* E［ Vj+1］ =E［ Wj+1］ W* () 將上面兩式代入 ()式中，得到 ][]2[][ 1 jxxj VERIVE ????它的遞推解是 0]2[][ VRIVE jxxj ???令 Rxx=QΛ Q T, Λ =QRxxQT () 第三章 自適應(yīng)數(shù)字濾波器 39。],[39。 39。39。239。1T QVVvvvVQV N ??? ?得到 39。039。 ]2[][ VΛI(xiàn)VE jj ???() () 再將 ()、 ()和 ()式代入上式，得到 E［ Wj］ =W*+Q［ I2μΛ］ j Q1(W0W*) () 對(duì)比 ()式 ， 說明 LMS算法加權(quán)矢量的統(tǒng)計(jì)平均值的過渡過程和最陡下降法加權(quán)矢量的過渡過程是一樣的 。 換句話說 ， LMS算法加權(quán)矢量是在最陡下降法加權(quán)矢量附近隨機(jī)變化的 ， 其統(tǒng)計(jì)平均值等于最陡下降法加權(quán)矢量 ， 那么 ， 其收斂條件同樣為 m a x10?? ??() 第三章 自適應(yīng)數(shù)字濾波器 在滿足收斂條件的情況下，才有下式： *][l i m WWEjj ???由于最大的特征值 λ max不可能大于 R的跡 (R的主對(duì)角線元素之和 )， 即 ? ????? )()(tr)()(trm a x 的對(duì)角元素的對(duì)角元素 RRΛΛ?因此收斂條件可以表示為 )(10Rtr?? ?() 第三章 自適應(yīng)數(shù)字濾波器 對(duì)于橫向?yàn)V波器， 式中的跡是 NE［ x2j］，即 N倍的輸入功率， 那么 ][102jxNE?? ?() 實(shí)際中，通常 μ 選得很小，選 ][102jxNE??? ?() 同樣由 ()式，第 i個(gè)分量為 39。0]2[][ ijiji vIvE ????() 第三章 自適應(yīng)數(shù)字濾波器 同樣引入時(shí)常數(shù) τ i, ][][e][21)21l n (139。*39。039。 ijjτjijiiiivQEWWEvvE????????????() () () 同樣，第 i個(gè)權(quán)系數(shù)可以表示成 kτjNkikiji CwwE1* e][ ????() 第三章 自適應(yīng)數(shù)字濾波器 3 . LMS算法性能函數(shù)的過渡過程 ——學(xué)習(xí)過程 由于 LMS算法加權(quán)矢量的平均值的變化規(guī)律與最陡下降法的加權(quán)矢量一樣 ， 可以推想它的均方誤差也會(huì)按照最陡下降的均方誤差變化規(guī)律變化 。 下面進(jìn)行推導(dǎo) 。 按照 ()式， 信號(hào)誤差為 jTjjjTjTjjTjjjjjVXeWWXWXdWXdyde??????????opt**)()(() 第三章 自適應(yīng)數(shù)字濾波器 式中 ， eoptj=djXjTW*， 稱為最佳誤差信號(hào) ， 它對(duì)應(yīng)于最小均方誤差 ， 即 m in22opt ][][ jj eEeE ?按照 ()式寫出均方誤差表示式： ][2][][][ 22 jTjoptjTjjTjoptj VXeEVXXVEeEeE ?????假定 Xj和 Vj不相關(guān)，上式中最后一項(xiàng)為 0，那么 ][m i n jTjjTj VXXVE?? ??第三章 自適應(yīng)數(shù)字濾波器 239。1m i n39。39。m i nm i n])[[][][][][jiNiijTjjxxTjvEVΛEVEVERVE?????????????同樣，假設(shè)加權(quán)系數(shù)變化很小， Vj也變化很小， E［ Vj］ ≈ Vj，這樣： 類似前面的推導(dǎo)，得到 iii m s eNijojiiv????????412e][12239。m i n???? ??() () 對(duì)照最陡下降法性能曲線 ()式 ， LMS均方誤差變化規(guī)律和最陡下降法完全一樣 ， 學(xué)習(xí)曲線同樣近似為幾個(gè)不同時(shí)間常數(shù)的指數(shù)和 。 第三章 自適應(yīng)數(shù)字濾波器 4. 穩(wěn)態(tài)誤差和失調(diào)系數(shù) 由上面分析知道 ， 權(quán)矢量的平均值可以收斂到它的最佳值 ， 但權(quán)矢量變化過程是隨機(jī)的 ， 即使其平均值收斂到最佳值 ， 它仍然按照下式： Wj+1=Wj+2μ ejXj 隨機(jī)地進(jìn)行變化 ， 這樣使權(quán)矢量仍在最佳值附近隨機(jī)變化 ， 但均方誤差將大于最小均方誤差 ， 如圖 所示 。 為此 ， 引入失調(diào)系數(shù) M， M定義為 m inm in??? ??M() 第三章 自適應(yīng)數(shù)字濾波器 加 權(quán) 系 數(shù) 的 變 化均 方 誤 差 的 變 化??m i nvv ＝ w － wo p tj圖 LMS算法穩(wěn)態(tài)誤差 第三章 自適應(yīng)數(shù)字濾波器 可以推出 ［ 5］ 失調(diào)系數(shù)為 ????NiixxRM1]t r[ ???() 或者 M=μ NPin () 式中 ， N是濾波器的階數(shù) ， Pin是輸入信號(hào)功率 。 上式說明 μ 和輸入功率加大都會(huì)增加失調(diào)系數(shù) 。 在保證收斂的情況下 μ 加大 ， 會(huì)提高收斂速度 ， 也說明為了減小失調(diào)系數(shù) ， 應(yīng)該適當(dāng)選擇收斂速度 ， 以保證收斂速度和失調(diào)系數(shù)都滿足要求 。 第三章 自適應(yīng)數(shù)字濾波器 圖 是一個(gè) LMS自適應(yīng)濾波器的計(jì)算機(jī)結(jié)果 ［ 5］ ， 階數(shù) N=5， 其輸入是信號(hào)加白噪聲 ， 輸入信號(hào)功率為 1， 中心頻率是 (fs為采樣頻率 )， 噪聲功率為 , 輸入信號(hào)自相關(guān)函數(shù)的特征值為： 、 、 、 、 ， 權(quán)系數(shù)初始值取 0， μ =。 圖中畫出了一條樣本學(xué)習(xí)曲線和 150條樣本學(xué)習(xí)曲線的平均曲線 。 該圖表明個(gè)別學(xué)習(xí)曲線起伏較大 ，平均學(xué)習(xí)曲線起伏很小 ， 計(jì)算出的維納最小均方誤差為 96， 用 LMS算法得到的穩(wěn)態(tài)誤差大于該值 ， 按 ()式計(jì)算的失調(diào)系數(shù)是 %， 按計(jì)算機(jī)模擬結(jié)果測(cè)得的失調(diào)系數(shù)是 %。 第三章 自適應(yīng)數(shù)字濾波器 圖 LMS算法的學(xué)習(xí)曲線 1 . 5 01 . 2 51 . 0 00 . 7 50 . 5 00 . 2 5均方誤差0 20 40 60 80 1 0 0 1 2 0 1 4 0 1 6 0 1 8 0 2 0 0迭代次數(shù) n?m i n＝ 0 . 7 4 3 9 6 ( 維納解 )個(gè) 別 學(xué) 習(xí)曲 線平 均 學(xué) 習(xí)曲 線( 由 1 5 0 條 個(gè) 別 曲線 平 均 )第三章 自適應(yīng)數(shù)字濾波器 自適應(yīng)格型濾波器 前 、 1. 前向線性預(yù)測(cè)誤差濾波器 為了分析簡(jiǎn)單 ， 假設(shè)信號(hào)屬于實(shí)平穩(wěn)隨機(jī)信號(hào) 。 前向線性預(yù)測(cè)誤差濾波器直接由信號(hào)的線性一步預(yù)測(cè)導(dǎo)出 。 在維納濾波器一章我們已研究了信號(hào)的線性一步預(yù)測(cè)問題 ， 即由 x(n1), x(n2), …, x(np)預(yù)測(cè) x(n)， 其估計(jì)值 x(n)和預(yù)測(cè)誤差 ep(n)用下式表示： ^ ?????pkkp knxanx1, )()(?() ???????pkkpp knxanxnxnxne1, )()()(?)()(第三章 自適應(yīng)數(shù)字濾波器 由于假設(shè)了信號(hào)是實(shí)的 ， 式中預(yù)測(cè)誤差 ep(n)和系數(shù) ap,k均是實(shí)數(shù) 。 ()式表明 是由 n時(shí)刻以前的 p個(gè)數(shù)據(jù) x(n1)、 x(n2)…x(np)得到的估計(jì) ， 因此稱 為前向預(yù)測(cè)誤差 。 將前向預(yù)測(cè)誤差用 表示 ， 上式重寫為 )(? nx)(nep)(nefp?????pkkpfp knxanxne1, )()()(() 對(duì)上式進(jìn)行 Z變換，得到 kpkkpfp zzXazXzE?????1, )()()(() 令 ??????????pkkkppkkkpffpfzazazHzXzEzH0,1,1)()()()(() 第三章 自適應(yīng)數(shù)字濾波器 Hf(z)稱為前向預(yù)測(cè)誤差濾波器的系統(tǒng)函數(shù) 。 前向預(yù)測(cè)誤差濾波器的結(jié)構(gòu)圖如圖 。 圖 前向預(yù)測(cè)誤差濾波器 z － 1ap ,0……z － 1 z － 1ap ,1ap ,2ap , p － 1ap , px ( n ))( ne fp第三章 自適應(yīng)數(shù)字濾波器 用均方誤差最小的準(zhǔn)則求前向預(yù)測(cè)誤差濾波器的最佳系數(shù) ap,k, 0]))([(,2???kpfpaneEk=1, 2, … ， p () 將 ()式代入上式，得到 0)]()([ ??? knxneE fpk=1, 2, 3,… ， p () 上式表明前向預(yù)測(cè)誤差與用于預(yù)測(cè)的數(shù)據(jù)正交 ， 這就是對(duì)于前向預(yù)測(cè)誤差的正交原理 。 按照第二章的推導(dǎo) ， 前向預(yù)測(cè)誤差濾波器的最佳系數(shù) ap,kYuleWalker方程式 ， 重寫如下： 第三章 自適應(yīng)數(shù)字濾波器 ?????????????????pixxipxxppixxipxxirarpkikrakr1,21,)()0(,3,2,10)()(??() 將上式用矩陣方程表示為 ????????????????