【文章內(nèi)容簡(jiǎn)介】 的整個(gè)計(jì)算過(guò)程中非奇異。將上式修正為相關(guān)矩陣的過(guò)程叫做對(duì)角加載。 (532) )(n?河南工業(yè)大學(xué) 橫向?yàn)V波器抽頭輸入與期望響應(yīng)之間的 M １時(shí)間平均互相關(guān)向量 z(n)為 ????niin idiunz1* )()()( ? 它將不受正則化的影響 ， 此處依然假定使用預(yù)加窗法 。 根據(jù)前面討論過(guò)的最小二乘法 ， 可使用代價(jià)函數(shù)獲得最小值的最優(yōu) M １抽頭權(quán)向量 由正則方程定義 。 遞歸最小二乘問(wèn)題的正則方程可用矩陣形式寫(xiě)為： 這里的 Φ(n)和 z(n)分別由式 (532) 和式 (533) 決定。 )(^ nw)()()( nznwn ?? ?(533) (534) 河南工業(yè)大學(xué) 將對(duì)應(yīng)于 i=n的項(xiàng)與式 (532)右邊的求和項(xiàng)分開(kāi)，可寫(xiě)出 根據(jù)定義，上式右邊括號(hào)內(nèi)的表達(dá)式等于相關(guān)矩陣Φ(n1)。于是，可得用于更新抽頭輸入相關(guān)矩陣的遞歸公式 其中 Φ(n1)是相關(guān)矩陣的過(guò)去值，矩陣乘積 在更新過(guò)程中起著“修正”項(xiàng)的作用。注意，上式的遞歸過(guò)程與初始條件無(wú)關(guān)。 ?Φ(n)和 z(n)的遞歸算法 )()( nunu H)()()1()( nununn H????? ?)()(])([)(1111 nunuIiun HninHin ???? ?????? ???? (535) (536) 河南工業(yè)大學(xué) 類(lèi)似地，可以導(dǎo)出抽頭輸入與期望響應(yīng)之間互相關(guān)向量的更新公式 )()()1()( * ndnunznz ??? ? 為了按式 (532)計(jì)算抽頭權(quán)向量 的最小二乘估計(jì)，必須確定相關(guān)矩陣 Φ(n)的逆。然而在實(shí)際中，我們通常盡量避免這樣做，因?yàn)檫@種運(yùn)算非常耗時(shí)，特別是當(dāng)抽頭數(shù) M很大時(shí)。另外，我們希望能夠遞歸計(jì)算 n=1,2,… ∞時(shí)抽頭權(quán)向量 的最小二乘估計(jì)。我們發(fā)現(xiàn)，利用矩陣代數(shù)中矩陣求逆引理，可以實(shí)現(xiàn)上述兩個(gè)目標(biāo)。 )(^ nw)(^ nw(537) 河南工業(yè)大學(xué) 矩陣求逆引理 設(shè) A和 B是兩個(gè) M M正定陣，它們之間的關(guān)系為 其中， D是 N M正定陣， C是 M N矩陣。根據(jù)矩陣求逆引理，可將 A的逆矩陣表示為 該引理在此不做證明，書(shū)中有介紹。在下一節(jié)，我們將說(shuō)明怎樣應(yīng)用矩陣求逆引理，得到計(jì)算抽頭權(quán)向量 最小二乘解的遞歸公式。 HCCDBA 11 ?? ??BCBCCDBCBA HH 11 )( ?? ???)(^ nw(538) (539) 河南工業(yè)大學(xué) 指數(shù)加權(quán)遞歸最小二乘算法 假定相關(guān)矩陣 Φ(n)是非奇異的，因而它可逆。我們對(duì)式 (536) 所表示的遞歸方程應(yīng)用矩陣求逆引理，首先做如下設(shè)定 1)()1()(1????????DnuCnBnA?河南工業(yè)大學(xué) 然后，將這些定義代入矩陣求逆引理，可得計(jì)算相關(guān)矩陣逆陣的遞歸方程如下 為了方便計(jì)算，令 和 )()1()(1)1()()1()1()(11112111nunnunnunnnH ???????????????????????)()( 1 nnP ???)()1(1)()1()(11nunPununPnKH ????????(540) (541) (542) 河南工業(yè)大學(xué) )1()()()1()( 11 ???? ?? nPnunKnPnP H??)()]1()()()1([)()1()()()()1()(1111nunPnunKnPnunPnunKnunPnKHH????????????????用上面的定義，可將式 (540)改寫(xiě)為 M M矩陣 P(n)叫做逆相關(guān)矩陣， M 1向量 K(n)叫做增益向量，后面將會(huì)解釋這樣叫的原因。式 (543)是RLS算法的 Riccati 方程。整理 K(n)表達(dá)式 ，可得 (543) (544) 河南工業(yè)大學(xué) )()( 1 nnP ???)()()( 1 nunnK ???換句話說(shuō)，增益向量 K(n)可以定義為經(jīng)相關(guān)矩陣Φ(n)逆矩陣變換的抽頭輸入向量 u(n). 這一結(jié)論，連同 可以用來(lái)定義增益向量 我們可以將上式簡(jiǎn)化為 K(n)=P(n)u(n) (545) (546) 河南工業(yè)大學(xué) ?抽頭權(quán)向量的時(shí)間更新 下面，我們要導(dǎo)出更新抽頭權(quán)向量最小二乘估計(jì) 的遞歸公式。為此，用式 (532)、式 (537)、式 (541)來(lái)表示 抽頭權(quán)向量 n次迭代時(shí)的最小二乘估計(jì) )(^ nw)()()()1()()()()()()(?*1ndnunPnznPnznPnznnw????????(547) 將式 (547)右邊第一項(xiàng)中 P(n)用式 (543)代替，可得 河南工業(yè)大學(xué) )()()()1()()()1()()()()1()1()()()1()1()()()()1()1()()()1()1()(*^^*11*^ndnunPnwnunKnwndnunPnznnunKnznndnunPnznPnunKnznPnwHHH???????????????????????最后，應(yīng)用 P(n)u(n)等于增益向量 K(n)，可得更新抽頭向量的遞歸方程為 )()()1()]1()()()[()1()(^^^^nnKnwnwnundnKnwnw H???????????(548) (549) 河南工業(yè)大學(xué) 其中 )()1()()1()()()(^^nunwndnwnundnHT?????? ??是一個(gè)先驗(yàn)估計(jì)誤差。內(nèi)積 表示基于 n1時(shí)刻抽頭權(quán)向量最小二乘估計(jì)舊值的期望響應(yīng) d(n)的估值。 根據(jù)調(diào)整抽頭權(quán)向量的表達(dá)式和表示先驗(yàn)估計(jì)誤差的表達(dá)式可用圖 56 (a)所示的框圖表示遞歸最小二乘算法。 )()1(? H nunw ?(550) 河南工業(yè)大學(xué) 一般說(shuō)來(lái)，先驗(yàn)估計(jì)誤差 ξ(n)不同于下式 (551)的后驗(yàn)估計(jì)誤差 其計(jì)算設(shè)計(jì)抽頭權(quán)向量在時(shí)刻 n（當(dāng)前時(shí)刻）的最小二乘估計(jì)。實(shí)際上，我們可以將 ξ(n)視為更新抽頭權(quán)向量之前 e(n)的暫時(shí)值。但要注意的是，在導(dǎo)出式 (549)遞歸算法的最小二乘優(yōu)化問(wèn)題中，我們實(shí)際上是基于 e(n)而不是基于 ξ(n)使代價(jià)函數(shù)ξ(n)最小 。 )()(?)()( nunwndne H?? (551) 河南工業(yè)大學(xué) 輸出 橫向?yàn)V波器 輸入向量 U(n) )()1(^ nunw H ?自適應(yīng)權(quán)值控制機(jī)制 ∑ )1(^ ?nw誤差 ξ(n) 期望響應(yīng) d(n) 圖 56(a)框圖 RLS算法 + ∑ ∑ K(n) Iz1? )(nwH)(nd? )(n?? )()( nkn??)(ny?)(^ nw增益 圖 56(b)信號(hào)流圖 + )1(? ?nw河南工業(yè)大學(xué) ?RLS算法小結(jié) 式 (542) ,(543), (549), (550)組成了 RLS算法，并在表 (51)中總結(jié)。特別要注意的是，式(550)表述了該算法的濾波過(guò)程，據(jù)此激勵(lì)橫向?yàn)V波器以計(jì)算先驗(yàn)估計(jì)誤差 ξ(n) 。式 (549) 描述了算法的自適應(yīng)過(guò)程，據(jù)此可通過(guò)在其過(guò)去值的基礎(chǔ)上增加一個(gè)量來(lái)遞推抽頭權(quán)向量，該量等于先驗(yàn)估計(jì)誤差 ξ(n)復(fù)共軛與時(shí)變?cè)鲆嫦蛄?K(n)的乘積 (“增益向量 ” 由此得名 )。 河南工業(yè)大學(xué) 式 (542) ,(543),使得我們能夠更新增益向量本身。上述 RLS算法的一個(gè)重要特點(diǎn)是，每一次迭代中的相關(guān)矩陣 Φ(n)的逆矩陣為簡(jiǎn)單的標(biāo)量相除所代替。圖 (56a) 給出 RLS算法的框圖，圖 (56b)則是 RLS算法的信號(hào)流圖。 IPw1^)0(0)0(????表 1 RLS算法小結(jié) 算法初始化 ????時(shí)取大的正常數(shù)低時(shí)取小的正常數(shù)高S N RS N R?河南工業(yè)大學(xué) )1n(P)n(u)n(K)1n(P)n(P)n()n(K)1n(w)n(w)n(u)1n(w)n(d)n()n()n(u)n()n(K)n(u)1n(P)n(,2,1nH11^^^HH??????????????????????????和計(jì)算對(duì)每一時(shí)刻， ? 表 1續(xù) 河南工業(yè)大學(xué) 注意，在表 1的總結(jié)中，增益向量 K(n)的計(jì)算分兩步進(jìn)行： ?首先，計(jì)算用 表示的中間量。 ?第二，用 計(jì)算 K(n) 從有限精度運(yùn)算的角度看，分兩步計(jì)算 K(n)比直接用式 (525)計(jì)算 K(n)更可取。 為對(duì) RLS濾波器進(jìn)行初始化，需要指定兩個(gè)量： ?初始權(quán)向量 。習(xí)慣上令 。 ?初始相關(guān)矩陣 Φ(n)。令式 (58)中的 n=0，如果使用預(yù)加窗，可以得到 其中 δ 是正則化參數(shù)。參數(shù) δ的設(shè)定與信噪比有關(guān)，高信噪比時(shí)取小值；低信噪比時(shí)則取較大值。這樣做的合理性可以在正則化的意義得到證明。 )0(^w 0)0(^ ?wI??? )0()(n?)(n?河南工業(yè)大學(xué) 正則化參數(shù)的選擇 在 Moustakides的詳細(xì)研究 (1997)中，評(píng)價(jià)了在平穩(wěn)環(huán)境下 RLS算法的收斂性能， 它有兩個(gè)特殊的可變參數(shù)： ?抽頭輸入數(shù)據(jù)的信噪比 (SNR)，這個(gè)量由流行的運(yùn)行條件決定。 ?正則化參數(shù) δ ，它由設(shè)計(jì)人員控制。 河南工業(yè)大學(xué) 為了總結(jié) Moustakides研究成果的實(shí)驗(yàn)條件，用