【正文】 么根據(jù)上面的式子可以推出更新的抽頭權(quán)向量經(jīng)過(guò) n次迭代的最小二乘估計(jì))(nw? 如式（ 319）： )()()()1()()()()()()( 1ndnunPnznPnznPnznnw?????????? （ 319） 然后把 (319）化簡(jiǎn)得出： )()()()1()1()()()1()1()( ndnunPnznPnunknznPnw H ?? ??????? )()()()1()1()()()1()1( 11 ndnunPnznnunknzn H ??? ????????? )()()()1()()()1( ndnunPnwnunknw H ??? ????? （ 320） 將 )()()( nunPnk ? 帶入得： )()()1()]1()()()[()1()(nnknwnwnundnknwnw H??????????????? （ 321） 其中令 )()1()()1()()()( nunwndnwnundn HT ?????? ?? ?? ,其代表先驗(yàn)估計(jì)誤差，)()1( nunwH ?? 代表了 n1 時(shí)抽頭權(quán)值向量最小二乘估計(jì)舊值的期望響應(yīng) d(n)的估值。下圖 31 很好的表示這 RLS算法的這一過(guò)程。 西南科技大學(xué)本科 生畢業(yè)論文 21 橫 向 自 適 應(yīng) 濾 波 器自 適 應(yīng) 權(quán) 值 控 制 機(jī) 制 ∑+輸 入 向 量u ( n )期 望 響 應(yīng) d ( n )輸 出誤 差ξ ( n ) 圖 31 RLS 算法框圖 上面的幾個(gè)式子就組成了完整的 RLS算法 ,式（ 321）表示整個(gè)濾波器的濾波過(guò)程，可以根據(jù)這個(gè)式子對(duì)橫向?yàn)V波器計(jì)算先驗(yàn)估計(jì)誤差 )(n? ，可以根據(jù)此算法的自適應(yīng)過(guò)程在過(guò)去值的基礎(chǔ)上來(lái)增加一個(gè)量來(lái)遞推抽頭權(quán)向量，此時(shí)這個(gè)量和先驗(yàn)估計(jì)誤差 )(n? 復(fù)共軛與時(shí)變的增益向量 k(n)的乘積相等。其中， RLS算法的一個(gè)重要特性是：在每一次迭代中相關(guān)矩陣 )(n? 的逆矩陣是簡(jiǎn)單的標(biāo)向量相除的值。圖 32 給出了整個(gè) RLS 算法的信號(hào)流圖。 K ( n )∑ ∑+ 增 益單 位 負(fù) 反 饋d * ( n )ξ * ( n )ξ * ( n ) k ( n )y * ( n ) 圖 32 RLS 算法信號(hào)流圖 下面簡(jiǎn)單的對(duì) RLS 算法的整個(gè)過(guò)程做一個(gè)整理： 首先就是對(duì)所有的初始只進(jìn)行初始化：令 0)0( ??w ， IP 1)0( ??? ，其中 ? 是表示一個(gè)常數(shù)，當(dāng)信噪比高時(shí)其值取較小的正常數(shù)，相反，當(dāng)信噪比較高時(shí)，其中則取對(duì)較大的正常數(shù)。 然后取得每個(gè)時(shí)刻時(shí)的 d(n)和 u(n) 之后在每個(gè)時(shí)刻都更新增益矢量， )(n? 是一個(gè)中間量，用來(lái)計(jì)算 k(n) )()1()( nunPn ??? )1( ?? nw )()1( nunwH ?? Iz1? )(nuH )(nw? )1( ?? nw 西南科技大學(xué)本科 生畢業(yè)論文 22 )()( )()( nnu nnk H ?? ??? 接下來(lái)就對(duì)濾波器的參數(shù)進(jìn)行更新： )()()1()()()1()()(nnknwnwnunwndn H???????????? 最后每個(gè)時(shí)刻對(duì)逆矩陣進(jìn)行更新： )1()()()1()( 11 ???? ?? nPnunknPnP H?? 以上總結(jié)了 RLS 算法的整個(gè)流程，在清楚了解了整個(gè)過(guò)程后，接下來(lái)將研究如 何應(yīng)用于具體的信道模型中，下圖 33就是基于 RLS 算法的實(shí)驗(yàn)信道模型。 隨 機(jī) 噪 聲發(fā) 聲 器 1信 道∑∑ 隨 機(jī) 噪 聲發(fā) 聲 器 2延 時(shí)基 于 算 R L S 法的 自 適 應(yīng) 橫向 均 衡 器X nν （ n ） + 圖 33 基于 RLS 算法的自適應(yīng)均衡器信道模型框圖 可以看出基于 RLS 算法的自適應(yīng)均衡器的模型與 LMS 算法的有所相同， 系統(tǒng)中 仍然使用兩個(gè)隨機(jī)數(shù)發(fā)聲器，一個(gè) 用來(lái)產(chǎn)生 測(cè)試信號(hào) nx ，是關(guān)于伯努利序列 1??nx ，其值同樣有零均值和單位方差。另一個(gè)表示模擬接收器中加性白噪聲的影響 )(nv ，其值也具有零均值，方差 2v? ，其中方差與信噪比有關(guān)。 兩個(gè)信號(hào)仍然相互獨(dú)立?；赗LS算法的自適應(yīng)均衡器也同樣用來(lái)改變信道畸變的情況，所以這需要用來(lái)上面總結(jié)的 RLS 算法的過(guò)程。這也是與 LMS 算法最根本的不同之處。 經(jīng)過(guò)合理的延遲，隨機(jī)發(fā)聲器 1也會(huì)提供用做訓(xùn)練序列的自適應(yīng)濾波器響應(yīng)。 同樣的首先 把加入到信道的隨機(jī)序列設(shè)為 ??nx ，由伯努利序列組成， 1??nx 且具有零均值和單位方差，令信道的脈沖響應(yīng)可以用升余弦表示 然后用 參數(shù) W是控制均衡器抽頭輸入的相關(guān)矩陣的特征值分布 )(R? 。其中 把均衡器的抽頭系數(shù) M設(shè)為 11， 利 西南科技大學(xué)本科 生畢業(yè)論文 23 用 信道脈沖響應(yīng) nh 的對(duì)稱(chēng)與 均衡器的最優(yōu)抽頭權(quán)值 的對(duì)稱(chēng)實(shí)現(xiàn)兩者 求卷積后得到均衡器的 期望響應(yīng)，通過(guò)選擇橫向?yàn)V波器中均衡器中點(diǎn)的合適延時(shí) RLS 算法可以提供信道響應(yīng)最小相位分量和其值之逆 。在研究 RLS 算法的過(guò)程中為了方便研究將加權(quán)因子 λ的值設(shè)定為 1，并針對(duì)信噪比，特征值擴(kuò)散度等參數(shù)變量來(lái)研究其性能的優(yōu)劣。具體的過(guò)程會(huì)在第四中詳細(xì)分析。 RLS 算法收斂分析 基于以上對(duì) RLS 算法的分析推理，可以想到其方法與 LMS 算法有什么不同，兩種方法其實(shí)最大的性能差別就在收斂性的比較。下面 則需要 對(duì) RLS 算法的收斂性做一個(gè)分析。為了對(duì) RLS 算的收斂性進(jìn)行分析，首先做出三個(gè)假設(shè)，為了方便研究討，在這三種假設(shè) 的情況下λ的值都設(shè)定是 1。 期望響應(yīng) d(n)與抽頭輸入向量 u(n)之間的關(guān)系可以用式（ 322）多重的線(xiàn)性回歸來(lái)描述，其中定義了回歸參數(shù)向量用 0w 表示，測(cè)量噪聲用 )(0ne 表示，其代表的是白噪聲（白噪聲的方差為 20? ，均值是 0，與回歸向量 u(n)無(wú)關(guān)）。 )()()( 0 nenuwnd H ?? （ 322） 輸入向量 u(n)是隨機(jī)生成的，而且 u(n)自相關(guān)函數(shù)是各態(tài)經(jīng)歷的。也就是說(shuō)用時(shí)間平均來(lái)替代集平均，可以把輸入向量 u(n)的集平均相關(guān)矩陣寫(xiě)車(chē)如下： )(1 nnR ?? 當(dāng) nM （ 323） 其中的 )(n? 是關(guān)于 u(n)的時(shí)間平均相關(guān)矩陣，其中 nM 是為了使橫向?yàn)V波器的 每個(gè)抽頭上都能夠有輸入信號(hào)，也就是說(shuō)是將其近似的將隨著時(shí)間增加來(lái)改善其矩陣的值。 加權(quán)誤差向量 )(n? 的波動(dòng)相對(duì)輸入向量 )(nu 的波動(dòng)較慢。因?yàn)榧訖?quán)誤差向量)(n? 是 RLS 算法經(jīng)過(guò) n 次迭代后多次累加的值，那么可以用式（ 324）表示： ?? ?? ???? ni iiknwwn 10 )()()0()()( ??? （ 324） 可以看出 k(i)和 )(i? 都和 u(i)有關(guān)系，但是上式對(duì) )(n? 具有平滑作用，也可以 說(shuō)， RLS 自適應(yīng)濾波器相當(dāng)于是 時(shí)變且低通濾波器?；谝陨系娜齻€(gè)假設(shè)，現(xiàn)在將對(duì)RLS算法的均值收斂性進(jìn)行研究。首先可以先從 )(nw? 開(kāi)始，可以知道： )()()( 1 nznnw ?? ?? nM (325) 其中當(dāng) 1?? 時(shí)可以得出： 西南科技大學(xué)本科 生畢業(yè)論文 24 )0()()()(1 ???? ??niH iuiun （ 326） ???? ni idnunz 1 )()()( （ 327） 然后應(yīng)用上面的假設(shè)一，也就是式（ 322），對(duì)式（ 327）整理得： ???????????????nininiHieiuwwnieiuwiuiunz10000101)()()0()()()()()()( （ 328） 因?yàn)?1?? ，所以可以寫(xiě)成如下形式： ?????????????????????????niniieiunwnwneiunwnwnnnw1010101010101)()()()0()()()()()0()()()()( (329) 然后對(duì)式（ 329）的兩邊求期望，并將上面的兩個(gè)假設(shè)與其整理后可以得到： pnwwRnwwRnwnwE????????????0010010 )0(1)]([ Mn? (330) 其中 p代表的是期望響應(yīng) d(n)與輸入向量 u(n)之間的集平均互相關(guān)向量，從上面的式子可以看出 RLS 算法在均值方面是收斂的。當(dāng)滿(mǎn)足 n 大于 M時(shí)，因?yàn)槌跏蓟?)0(?的值，所以在 )(nw? 是有偏差的，只有當(dāng) n的值無(wú)限大時(shí)，偏差則趨近 0。下面再來(lái)看看 RLS 算 法的均方偏差，首先定義出加權(quán)誤差相關(guān)矩陣為： ]))())(([()]()([)( 00 HH nwwnwwEnnEnK ?? ???? ?? （ 331） 把上面所求的 )(nw? 的值帶入得： ])()()()()()([)( 1 1 0011 ? ?? ? ??? ??? ni nj H jejenjuiunEnK （ 332） 這個(gè)式子中忽略了初始化的影響，在第一個(gè)假設(shè)的情況下，可以知道輸入向量u(n)和 )(1n?? 與測(cè)量噪聲 )(0ne 沒(méi)有關(guān)系，所以將 K(n)用兩個(gè)期望表示為： ? ?)()()()()()()(001 1 11 jeieEnjuiunEnKninjH ?? ??? ?????? ??? ? ? （ 333） 西南科技大學(xué)本科 生畢業(yè)論文 25 在測(cè)量誤差時(shí)使用的是白噪聲，其方差為 20? ，根據(jù)第一個(gè)假設(shè)，可以得出有： ? ??? )()( 00 jeieE 所以加權(quán)誤差相關(guān)矩陣可以寫(xiě)成： ?????? ??? ?? ??niH niuiunEnK11120 )()()()()( ? ? ?? ?)()()()(1201102nEnnnE??????????? （ 335） 最后利用第二個(gè)假設(shè)可以改寫(xiě)為： 1201)( ?? RnnK ? nM （ 336） 最后 就 可以根據(jù)均方偏差的定義，如（ 327），用矩陣求逆算子表示 RLS 算法 的均方偏差，其中用 i? 表示集平均相關(guān)矩陣 R特征值如式（ 328）： )]([)]()([)( nKtrnnEn H ?? ??? （ 337） ??? ?? Mi inRtrnn 120210 11][1)( ???? nM （ 338） 從上面的式子可以得出以下兩點(diǎn)結(jié)論： 均方偏差 )(n? 是通過(guò)最小特征值 min? 的倒數(shù)來(lái)放大的，從更深的層面來(lái)講，就是對(duì)于一節(jié)近似中， RLS 算法對(duì)特征值擴(kuò)散的銘感程度與最小特征值的倒數(shù)是成正比關(guān) 系的，所以通常情況下收斂性回受到病態(tài)最小二乘問(wèn)題影響。 均方偏差 )(n? 與迭代次數(shù) n有關(guān)，兩者之間是一個(gè)線(xiàn)性衰減的關(guān)系，從另一方面來(lái)講就是說(shuō) RLS 算法得到的估計(jì) )(nw? 是隨著時(shí)間線(xiàn)性的按照均方值收斂于多重線(xiàn)性回歸模型的參數(shù)向量 0w 。 本章小結(jié) 本章中導(dǎo)出了遞推的最小二乘 （ RLS）算法，可以看出 RLS 算法與 LMS 算法的區(qū)別就在于 LMS算法的步長(zhǎng)參數(shù) u在 RLS算法中用濾波器的輸入矢量 x(n)的相關(guān)矩陣的逆替代了。另外針對(duì) RLS 算法的收斂性做了分析及推理，可以看出總的來(lái)說(shuō) RLS 算雖然相對(duì) LM