引 言
自適應(yīng)濾波器一直是信號處理領(lǐng)域的研究熱點之一,，經(jīng)過多年的發(fā)展，已經(jīng)被廣泛應(yīng)用于數(shù)字通信,、回聲消除,、圖像處理等領(lǐng)域。自適應(yīng)濾波算法的研究始于20世紀(jì)50年代末,，Widrow和Hoff等人最早提出最小均方算法(LMS),。算法由于結(jié)構(gòu)簡單，計算量小,，易于實時處理,，因此在噪聲抵消，譜線增強(qiáng),，系統(tǒng)識別等方面得到了廣泛的應(yīng)用,。為了克服定步長LMS算法中收斂速度、收斂精度及跟蹤速度等對步長大小選取相互矛盾的缺點,，人們提出了許多變步長LMS算法,，但是，當(dāng)輸入信號具有強(qiáng)相關(guān)性時，例如語音信號,，LMS算法及NLMS算法的收斂速度將急劇下降,。因此，本文對進(jìn)入自適應(yīng)濾波器的輸入信號首先進(jìn)行解相關(guān)預(yù)處理,，再對自適應(yīng)濾波器進(jìn)行歸一化LMS,。
很多自適應(yīng)濾波器是在通用DSP處理器上通過編寫軟件代碼實現(xiàn)。這種濾波器在系統(tǒng)實時性要求不高的情況下,，能很好地滿足濾波要求,，但在實時性要求較高的場合(如3G移動通信，網(wǎng)絡(luò)會議等),，這種濾波器在處理速度和抗干擾性能等方面已不能滿足要求,。而用FPGA硬件實現(xiàn)的數(shù)字濾波器就能很好的解決并行性和速度問題,，而且其具有靈活的可配置特性和優(yōu)良的抗干擾能力,，因而自適應(yīng)濾波器的FPGA實現(xiàn)成為近年來數(shù)字濾波器中研究的熱點。Altera公司開發(fā)的DSP Builder可以完成FPGA的DSP系統(tǒng)設(shè)計,，它結(jié)合Matlab／Simulink提供的庫元件和強(qiáng)大的計算仿真能力,，進(jìn)行圖形建模和仿真，還可以完成從Simulink模塊文件向VHDL語言的自動轉(zhuǎn)換,。
1 算法概述
在標(biāo)準(zhǔn)LMS算法中,，其迭代公式如下：
式中：ω(n)為自適應(yīng)濾波器在時刻n的權(quán)矢量；x(n)為時刻n的輸入信號矢量,；d(n)為參考信號,；e(n)是誤差信號；μ是步長迭代因子,，為了增加LMS算法的收斂速度,，文獻(xiàn)[1-4]提出了不同的變步長算法代替定步長因子μ，但是很多變步長因子的計算量大,，存在開方和指數(shù)運算,，故在能滿足要求的情況下，歸一化LMS (NLMS)依然是一個很好的選擇,，即用μn／[xT(n)x(n)]代替定步長因子μ為了控制失調(diào),，在更新公式中引入一個固定的收斂因子μn，這是因為所有的倒數(shù)都是基于平方誤差的瞬時值而不是均方誤差(MSE),。對于NLMS算法,，當(dāng)輸入信號具有強(qiáng)相關(guān)性時，LMS算法及NLMS算法的收斂速度將急劇下降,。因此對進(jìn)入自適應(yīng)濾波器的輸入信號首先進(jìn)行解相關(guān)預(yù)處理,，再對自適應(yīng)濾波器進(jìn)行歸一化LMS，便得解相關(guān)歸一化LMS(DNLMS)即用式(4)代替式(3)：
2 基于FPGA的系統(tǒng)開發(fā)流程
在早期PFGA的DSP系統(tǒng)設(shè)計過程中，DSP開發(fā)者直接使用VHDL或Verilog HDL語音進(jìn)行開發(fā),。這種開發(fā)難度較大,，開發(fā)效率較低，為了克服這些問題,，出現(xiàn)了許多新的基于FPGA的DSP開發(fā)工具,，如Altera公司的DSP Builder、Xilinx公司的System Generator等,。圖1是Altera公司的FPGA系統(tǒng)級開發(fā)流程,。
DSP Builder是Altera公司推出的一個面向DSP開發(fā)的系統(tǒng)級工具。它是作為Matlab／Simuli-nk的一個工具箱出現(xiàn)的,。它從利用Matlab建立DSP電路模型開始,。電路模型可方便地利用Simulink與DSPBuilder中提供的功能模塊和IP核進(jìn)行設(shè)計。電路中的功能模塊和IP的技術(shù)參數(shù),、數(shù)據(jù)格式,、數(shù)據(jù)類型和總線寬度等都能直接設(shè)置。電路模型建立后就能進(jìn)行系統(tǒng)級的仿真,。仿真通過后可以運行Signal Compiler,，將模型轉(zhuǎn)化成RTL級的VHDL代碼和Tcl(工具命令語言)腳本語言。調(diào)用QuartusⅡ進(jìn)行編譯,，QuartusⅡ根據(jù)網(wǎng)表文件及設(shè)置的優(yōu)化約束條件進(jìn)行布線布局和優(yōu)化設(shè)計的適配,，最后生成編程文件和仿真文件。
3 系統(tǒng)建模與仿真
根據(jù)DNLMS算法的思想,，結(jié)合FPGA硬件設(shè)計的特點,。
NLMS算法主要包括延時模塊，乘法模塊,，加法模塊,，除法模塊等。為了讓整個系統(tǒng)看起來清晰,，降低原理圖的復(fù)雜性,，也便于多次應(yīng)用，本文將幾個多次使用的模塊封裝成子模塊,。下面將對各個子模塊作簡單的介紹,。
根據(jù)算法的要求，每次新的數(shù)據(jù)x(n)輸入,，數(shù)據(jù)T,，S，M都需要重新更新一次,，計算量大,，消耗的乘法塊多。為了解決這個問題。根據(jù)T,，S,，M的運算特點，每次更新只需要一個乘法塊和一次加法,，就能完成8次乘法和一次加法的運算,。具體子模塊如圖2所示。
通過圖2的優(yōu)化,，對于8階濾波器,，乘法塊的使用從24個減少到2個，計算速度將得到很大的提高,。ω的迭代運算通過圖3實現(xiàn),。
在頂層模塊的設(shè)計中，加入了PLL模塊和SignalComplier,，根據(jù)各模塊仿真周期,，在PLL模塊中設(shè)置好不同的周期，利用Simulink中的顏色標(biāo)識采樣頻率的高低,。頂層文件的設(shè)計如圖4所示,。
在頂層結(jié)構(gòu)圖中,，子模塊8為方程(9)的實現(xiàn)過程,。
4 計算機(jī)仿真
為了充分模擬相關(guān)性較強(qiáng)的信號，在輸入端加如下的信號源：
x(n)=0.99x(n-1)+v(n)
式中：v(n)為高斯白噪聲,；x(n)的初始值為500,，具體的產(chǎn)生電路如圖5所示。
信號產(chǎn)生波形以及加噪后的仿真波形如圖6所示,。
通過濾波器濾波的結(jié)果如圖7所示,。
在圖7的信號依次為參考信號、加噪聲信號,、誤差信號和濾波信號,。濾波效果良好，對于強(qiáng)相關(guān)的信號,，也有很好的濾波效果,。通過QuartusⅡ編譯，濾器器的速度為35.24 MHz,，編譯后生產(chǎn)如圖8所示的塊,，它作為一個符號模塊可供更高級系統(tǒng)設(shè)計時調(diào)用，通過對該塊的修改和增加濾波器的階數(shù),，可以實現(xiàn)對語音的回聲消除,。
5 結(jié) 語
自適應(yīng)濾波器的研究是近年來研究的熱門方向，在PFGA中實現(xiàn)高速的自適應(yīng)濾波器的設(shè)計更是一個熱點，在此采用DSP Builder系統(tǒng)建模的方法,，在FPGA中實現(xiàn)了歸一化LMS算法,，實驗結(jié)果表明：用DSPBuilder設(shè)計的8階DNLMS算法比用底層的VHDL代碼設(shè)計效率更高，靈活性更大,，速度更快,。同樣比通用的DSP芯片設(shè)計的8階NLMS濾波器處理速度快了20多倍。如果將8階NLMS濾波器擴(kuò)展到512階或者1 024階,，可以很好地應(yīng)用于自適應(yīng)回聲消除等很多自適應(yīng)濾波的領(lǐng)域,。故有一定的參考價值。