自適應(yīng)濾波器在現(xiàn)代數(shù)字信號處理中的應(yīng)用極為廣泛，分成系統(tǒng)辨識、逆模型、預(yù)測和干擾消除4大類別[1]。對于實(shí)時(shí)性要求很高的場合，往往需要使用硬件設(shè)計(jì)自適應(yīng)濾波器?；?a class="innerlink" href="http://ihrv.cn/tags/FPGA" title="FPGA" target="_blank">FPGA的自適應(yīng)濾波器近年來被不斷地研究和改進(jìn)。

　　使用FPGA中的邏輯單元(LE)很容易實(shí)現(xiàn)乘累加(MAC)運(yùn)算，從而將采用最小均方(LMS)算法的自適應(yīng)濾波器直觀地映射到FPGA上。但是，基于MAC方案的自適應(yīng)濾波器的時(shí)鐘頻率會(huì)很高，功耗會(huì)很大。ALLRED D J等人提出了利用雙LUT結(jié)構(gòu)的DA方案來實(shí)現(xiàn)LMS自適應(yīng)濾波器[2]，雖然提高了吞吐量，降低了功耗，但其控制模塊、地址循環(huán)設(shè)置都比較復(fù)雜；Guo Rui等人提出采用單LUT結(jié)構(gòu)的DA方案來實(shí)現(xiàn)LMS自適應(yīng)濾波器[3]，盡管節(jié)省了近一半的LE，但其處理速度有所下降，并且偏移二進(jìn)制編碼（OBC）的硬件編程頗為復(fù)雜。

　　本文提出了基于雙LUT模塊的特殊DA方案，旨在將LMS順序執(zhí)行算法更高效地映射到FPGA上，充分利用了DA的基本特性以及對LMS算法的改進(jìn)來達(dá)到節(jié)省資源、提高速度的目標(biāo)。

1 自適應(yīng)濾波器的算法

　　設(shè)x(n)和w(n)分別表示自適應(yīng)濾波器的輸入信號和權(quán)值向量，定義：x(n)=[x(n)，x(n-1)，…，x(n-L+1)]T，w(n)=[w0(n)，w1(n)，…，wL-1(n)]T，則LMS算法的基本形式如下。

　　濾波輸出：

　　其中，y(n)是自適應(yīng)濾波器的輸出信號，d(n)是期望信號，e(n)表示誤差信號，是收斂因子，L是自適應(yīng)濾波器的抽頭數(shù)。雖然LMS算法僅涉及乘法、加減和迭代3種運(yùn)算，但它是一種嚴(yán)格順序執(zhí)行的自適應(yīng)算法[1]。倘若利用MAC方案直接將LMS自適應(yīng)濾波器映射到FPGA上，編程雖然簡單，但難以發(fā)揮FPGA并行處理的優(yōu)勢。

　　延遲最小均方（DLMS）算法是對LMS算法的一種改進(jìn)，該算法最大的優(yōu)點(diǎn)就是可以并行執(zhí)行[4]。本文所使用的DLMS算法的權(quán)值更新公式為：

　　而濾波輸出、誤差計(jì)算公式與LMS算法相同。由式(4)知，對第k個(gè)權(quán)值wk(n)進(jìn)行迭代更新時(shí)，必須對相應(yīng)的x(n)延遲k個(gè)采樣周期，k=0，1，…，L-1。圖1給出了DLMS自適應(yīng)濾波器的模塊圖。

　　圖2、圖3分別給出了4抽頭DLMS自適應(yīng)濾波器的抽頭權(quán)值以及輸入、輸出信號的Matlab仿真曲線。由圖2知，自適應(yīng)濾波器的4個(gè)抽頭權(quán)值的總體變化趨勢漸趨平緩。由圖3知，輸出信號相對于輸入信號而言抖動(dòng)幅度逐漸變小，而且愈發(fā)類似于正弦信號。但由于該自適應(yīng)濾波器只有4個(gè)抽頭，輸入信號的信噪比只有3 dB，所以自適應(yīng)濾波器的權(quán)值在520個(gè)采樣點(diǎn)以后才比較平緩，并且輸出信號與期望信號之間仍存在較明顯的差別。若能提高自適應(yīng)濾波器的抽頭數(shù)，或改用信噪比更高的輸入信號，則會(huì)加快抽頭權(quán)值的平緩速度，并且輸出信號將更加接近正弦信號。

2 特殊DA方案

　　在使用DA方案優(yōu)化基于FPGA的權(quán)值固定的經(jīng)典濾波器時(shí)，首先把濾波器的權(quán)值系數(shù)有規(guī)律地存儲在LUT模塊中，然后利用輸入信號的二進(jìn)制形式從最低位到最高位逐一訪問LUT模塊的存儲條目，最后對LUT模塊的各個(gè)輸出量進(jìn)行移位、累加運(yùn)算，得到最終的濾波輸出結(jié)果[5]。而自適應(yīng)濾波器的權(quán)值是不斷更新變化的，必須對存儲在LUT模塊中的權(quán)值進(jìn)行迭代更新，這時(shí)可以考慮使用兩個(gè)LUT模塊聯(lián)合運(yùn)作。

　　設(shè)輸入信號x(n)被量化成B bit二進(jìn)制補(bǔ)碼的形式，最左邊一位是符號位，則：

　　Qj存儲在LUT模塊中時(shí)沒有正負(fù)之別，在LUT模塊輸出端的后面引入一個(gè)符號控制端來判定對Qj做加法還是減法[6]。該LUT模塊定義為核心查找表(K-LUT)。

　　圖4給出了L=4、B=4的DLMS自適應(yīng)濾波器設(shè)計(jì)中K-LUT及其外圍模塊的結(jié)構(gòu)圖。圖中的PSC模塊表示并/串轉(zhuǎn)換器，2-1模塊表示右移器。

　　K-LUT的第r個(gè)條目的更新式如下：

　　該LUT模塊定義為輸入查找表（I-LUT）。當(dāng)檢測到自定義的中間信號clk_data的上升沿時(shí)，新的輸入信號將被讀入，I-LUT的內(nèi)容隨即自動(dòng)更新。

　　對于L=4、B=4的情況，由Quartus II集成開發(fā)軟件的時(shí)序仿真結(jié)果知，clk_data的周期應(yīng)設(shè)置成為時(shí)鐘信號clk周期的6倍。第n時(shí)刻I-LUT的各個(gè)條目存儲內(nèi)容如表1所示。當(dāng)n+1時(shí)刻來臨時(shí)，x(n)、x(n-1)、x(n-2)、x(n-3)自動(dòng)變成x(n+1)、x(n)、x(n-1)、x(n-2)。

　　一旦n時(shí)刻I-LUT的內(nèi)容更新并且y(n)、e(n)、e(n)的計(jì)算都已完成，那么K-LUT的內(nèi)容將被更新。通過設(shè)置合適的clk_data周期來保證K-LUT(n)到K-LUT(n+1)的更新在n+1時(shí)刻到來之前全部完成。式(8)則變?yōu)椋?/p>

　　K-LUT(r)(n+1)=K-LUT(r)(n)+?滋e(n)I-LUT(r)(n)(9)

　　K-LUT的更新步驟如下：(1)讀取K-LUT(n)和I-LUT(n)中具有相同存儲地址的條目；(2)用e(n)乘以I-LUT(n)的輸出；(3)將步驟(2)中的乘積與K-LUT(n)的輸出進(jìn)行求和運(yùn)算；(4)將步驟(3)中的和存儲到K-LUT(n+1)中相同存儲地址的條目。其中?滋和e(n)相乘是通過對e(n)進(jìn)行移位運(yùn)算實(shí)現(xiàn)的，值選成2的冪。

　　圖5為本DA方案的簡化流程圖，可以看出I-LUT的更新由clk_data控制，而濾波運(yùn)算、誤差計(jì)算、e(n)的乘積運(yùn)算都由clk控制。

3 基于FPGA實(shí)例驗(yàn)證

　　由于FPGA中的LUT通常具有4個(gè)輸入端，共24個(gè)存儲條目，故本文以設(shè)計(jì)4階基于FPGA的DLMS自適應(yīng)濾波器為例，這樣使用兩個(gè)LUT模塊就能實(shí)現(xiàn)最基本、最小規(guī)模的自適應(yīng)濾波器。其他的高階設(shè)計(jì)都以4階為基礎(chǔ)。

　　表2給出了分別采用MAC方案和本DA方案設(shè)計(jì)出的DLMS自適應(yīng)濾波器的各項(xiàng)性能指標(biāo)。這里取1/4。必須指出，系統(tǒng)時(shí)鐘clk的最高頻率與系統(tǒng)功耗密切相關(guān)，頻率越高，功耗越大；用于控制輸入信號讀取速度的clk_data的最高頻率表征了DLMS自適應(yīng)濾波器的處理速度。

　　由表2知，本DA方案設(shè)計(jì)出的DLMS自適應(yīng)濾波器不但消耗了更少的LE和寄存器數(shù)量，而且獲得了更低的功耗和更高的處理速度，優(yōu)勢明顯。表3給出了同樣采用本DA方案但值不同時(shí)的DLMS自適應(yīng)濾波器的性能比較。

　　由表3知，取不同值時(shí)采用本DA方案設(shè)計(jì)DLMS自適應(yīng)濾波器所消耗的LE和寄存器數(shù)目是相同的，所以在保證DLMS算法收斂前提下，的取值不會(huì)影響設(shè)計(jì)自適應(yīng)濾波器所消耗的硬件資源，但是自適應(yīng)濾波器的功耗和處理速度會(huì)有少許差別。

　　本文提出的特殊DA方案的基本結(jié)構(gòu)主要包含K-LUT和I-LUT模塊，重點(diǎn)討論了4階DLMS自適應(yīng)濾波器的設(shè)計(jì)方案。通過驗(yàn)證，采用本DA方案設(shè)計(jì)基于FPGA的DLMS自適應(yīng)濾波器相對于MAC方案具有更高的處理速度和更低的資源消耗。本方案推廣到更高階設(shè)計(jì)時(shí)，總體思路不變，但最好要對K-LUT和I-LUT模塊進(jìn)行分解，具體辦法有待進(jìn)一步探究。

參考文獻(xiàn)

　　[1] HAYKIN S.自適應(yīng)濾波器原理(第四版)[M].鄭寶玉，譯.北京：電子工業(yè)出版社，2010.

　　[2] ALLRED D J，HEEJONG Y，KRISHNAN V，et al.A novel high performance distributed arithmetic adaptive filter imple-mentation on FPGA[J].IEEE Signal Processing，2004，5(5)：161-165.

　　[3] Guo Rui，DEBRUNNER L S.A novel adaptive filter imple-mentation scheme using distributed arithmetic[J].IEEE Signals，Systems and Computers，2011，6(11)：160-164.

　　[4] Long Guozhu，Ling Fuyun，PROAKIS J G.The LMS algorithmwith delayed coefficient adaptation[J].IEEE Transactions on Acoustics，Speech and Signal Processing，1989，37(9)：1397-1405.

　　[5] UWE M B.Digital signal processing with field programmablegate arrays[M].Third Edition.Germany Berlin：Springer- Verlag，2007.

　　[6] ALLRED D J，YOO H，KRISHNAN V，et al.LMS adaptive filters using distributed arithmetic for high throughput[J].IEEE Transactions on Circuits and Systems，2005，52(7)：1327-1337.