基于DSP的自適應(yīng)濾波器的算法研究

1、基于DSP的自適應(yīng)濾波器的算法研究 作者： 日期：12 個(gè)人收集整理 勿做商業(yè)用途基于DSP的自適應(yīng)濾波器的算法研究梅娟 201310302003摘 要:隨著DSP技術(shù)的高速發(fā)展,人們對(duì)信號(hào)處理的實(shí)時(shí)性、準(zhǔn)確性和靈活性的要求越來越高,DSP技術(shù)在信號(hào)處理中的地位也越來越重要。自適應(yīng)濾波器是一種復(fù)雜的算法,設(shè)計(jì)它是為了在均衡信道,抵消回波，增強(qiáng)譜線,抑制噪聲等方面有所應(yīng)用。而自適應(yīng)濾波器的實(shí)現(xiàn)主要采用最小均方誤差算法完成。自適應(yīng)算法通過調(diào)整濾波器系數(shù)來實(shí)現(xiàn)可以更好地跟蹤信號(hào)的變化,最終實(shí)現(xiàn)自適應(yīng)濾波。關(guān)鍵詞:自適應(yīng)濾波器；數(shù)字信號(hào)處理器；最小均方誤差算法；Matlab仿真Algorithm Re

2、search of Adaptive Filter and DSP Simulation RealizationMei Juan 201310302003Abstract： With the rapid development of DSP technique， there is an increasing requirement of real time processing, accuracy and flexibility in signal processing。 DSP technique plays a more important role in signal processin

3、g。 Adaptive filters have been implemented with some complex algorithms in the past decade。 Designing it is for the balanced channel, counterbalances the echo, the enhancement spectral line， aspects noise elimination has the application and so on. But the realization of adaptive filter mainly uses th

4、e least mean square to complete. The least mean square is realized by adjusting filter and can track signal much better, finally realizes the adaptive filtering。個(gè)人收集整理,勿做商業(yè)用途本文為互聯(lián)網(wǎng)收集，請(qǐng)勿用作商業(yè)用途Keywords： adaptive filter； digital signal processor； least mean square error； matlab simulation0 引言濾波是電子信息處理領(lǐng)

5、域的一種最基本而又極其重要的技術(shù)。在有用信號(hào)的傳輸過程中，通常會(huì)受到噪聲或干擾的污染。利用濾波技術(shù)可以從復(fù)雜的信號(hào)中提取所需要的信號(hào)，同時(shí)抑制噪聲或干擾信號(hào)，以便更有效地利用原始信號(hào)。濾波器實(shí)際上是一種選頻系統(tǒng)，它對(duì)某些頻率的信號(hào)予以很小的衰減，讓該部分信號(hào)順利通過；而對(duì)其他不需要的頻率信號(hào)則予以很大的衰減，盡可能阻止這些信號(hào)通過。在電子系統(tǒng)中濾波器是一種基本的單元電路,使用很多，技術(shù)也較為復(fù)雜，有時(shí)濾波器的優(yōu)劣直接決定產(chǎn)品的性能，所以很多國(guó)家非常重視濾波器的理論研究和產(chǎn)品開發(fā).1 自適應(yīng)濾波器簡(jiǎn)介自適應(yīng)濾波器屬于現(xiàn)代濾波器的范疇,自適應(yīng)濾波器是相對(duì)固定濾波器而言的，固定濾波器屬于經(jīng)典濾波器，

6、它濾波的頻率是固定的,自適應(yīng)濾波器濾波的頻率則是自動(dòng)適應(yīng)輸入信號(hào)而變化的,所以其適用范圍更廣.在沒有任何關(guān)于信號(hào)和噪聲的先驗(yàn)知識(shí)的條件下,自適應(yīng)濾波器利用前一時(shí)刻已獲得的濾波器參數(shù)來自動(dòng)調(diào)節(jié)現(xiàn)時(shí)刻的濾波器參數(shù),以適應(yīng)信號(hào)和噪聲未知或隨機(jī)變化的統(tǒng)計(jì)特性,從而實(shí)現(xiàn)最優(yōu)濾波。所謂自適應(yīng)濾波,就是利用前一時(shí)刻已獲得的濾波器參數(shù)等結(jié)果,自動(dòng)地調(diào)節(jié)現(xiàn)時(shí)刻的濾波器參數(shù),以適應(yīng)信號(hào)和噪聲未知的或隨時(shí)間變化的統(tǒng)計(jì)特性，從而實(shí)現(xiàn)最優(yōu)濾波。自適應(yīng)濾波器實(shí)質(zhì)上就是一種能調(diào)節(jié)其自身傳輸特性以達(dá)到最優(yōu)化的維納濾波器。2 自適應(yīng)濾波原理自適應(yīng)濾波器的特性變化是由自適應(yīng)算法通過調(diào)整濾波器系數(shù)來實(shí)現(xiàn)的。一般而言，自適應(yīng)濾波器由

7、兩部分組成，一是濾波器結(jié)構(gòu),二是調(diào)整濾波器系數(shù)的自適應(yīng)算法。自適應(yīng)濾波器的結(jié)構(gòu)采用FIR或IIR結(jié)構(gòu)均可,由于IIR濾波器存在穩(wěn)定性問題，因此一般采用FIR濾波器作為自適應(yīng)濾波器的結(jié)構(gòu)。圖1給出了自適應(yīng)濾波器的一般結(jié)構(gòu)。圖1為自適應(yīng)濾波器結(jié)構(gòu)的一般形式，圖中x(n）為輸入信號(hào)，通過參數(shù)可調(diào)的數(shù)字濾波器后產(chǎn)生輸出信號(hào)y(n),將輸出信號(hào)y(n)與標(biāo)準(zhǔn)信號(hào)（或者為期望信號(hào))d(n)進(jìn)行比較,得到誤差信號(hào)e（n).e（n)和x（n）通過自適應(yīng)算法對(duì)濾波器的參數(shù)進(jìn)行調(diào)整,調(diào)整的目的使得誤差信號(hào)e（n)最小。圖1 自適應(yīng)濾波器的一般結(jié)構(gòu)自適應(yīng)濾波器設(shè)計(jì)中最常用的是FIR橫向型結(jié)構(gòu)。圖2是橫向型濾波器的結(jié)

8、構(gòu)示意圖。圖2 橫向型濾波器的結(jié)構(gòu)示意圖其中:為自適應(yīng)濾波器的輸入；為自適應(yīng)濾波器的沖激響應(yīng):;為自適應(yīng)濾波器的輸出:.3 自適應(yīng)濾波算法自適應(yīng)濾波器除了包括一個(gè)按照某種結(jié)構(gòu)設(shè)計(jì)的濾波器,還有一套自適應(yīng)的算法。自適應(yīng)算法是根據(jù)某種判斷來設(shè)計(jì)的。自適應(yīng)濾波器的算法主要是以各種判據(jù)條件作為推算基礎(chǔ)的。通常有兩種判據(jù)條件:最小均方誤差判據(jù)和最小二乘法判據(jù)。LMS算法是以最小均方誤差為判據(jù)的最典型的算法,也是應(yīng)用最廣泛的一種算法。最小均方誤差（Least Mean Square, LMS)算法是一種易于實(shí)現(xiàn)、性能穩(wěn)健、應(yīng)用廣泛的算法。所有的濾波器系數(shù)調(diào)整算法都是設(shè)法使接近，所不同的只是對(duì)于這種接近的評(píng)

9、價(jià)標(biāo)準(zhǔn)不同。LMS算法的目標(biāo)是通過調(diào)整系數(shù),使輸出誤差序列的均方值最小化,并且根據(jù)這個(gè)判據(jù)來修改權(quán)系數(shù),該算法因此而得名。誤差序列的均方值又叫/均方誤差0(Mean Square Error, MSE).理想信號(hào)與濾波器輸出之差的期望值最小，并且根據(jù)這個(gè)判據(jù)來修改權(quán)系數(shù)。由此產(chǎn)生的算法稱為L(zhǎng)MS.均方誤差表示為:對(duì)于橫向結(jié)構(gòu)的濾波器，代入的表達(dá)式:其中：為的自相關(guān)矩陣，它是輸入信號(hào)采樣值間的相關(guān)性矩陣。為互相關(guān)矢量，代表理想信號(hào)與輸入矢量的相關(guān)性.在均方誤差E達(dá)到最小時(shí)，得到最佳權(quán)系數(shù)。它應(yīng)滿足下式：=0這是一個(gè)線形方程組,如果R矩陣為滿秩的,存在，可得到權(quán)系數(shù)的最佳值滿足:.用完整的矩陣表示

10、為: 顯然為的自相關(guān)值,為與互相關(guān)值。在有些應(yīng)用中,把輸入信號(hào)的采樣值分成相同的一段（每段稱為一幀)，再求出R，P的估計(jì)值得到每幀的最佳權(quán)系數(shù)。這種方法稱為塊對(duì)塊自適應(yīng)算法.如語(yǔ)音信號(hào)的線性預(yù)測(cè)編碼LPC就是把語(yǔ)音信號(hào)分成幀進(jìn)行處理的。R,P的計(jì)算，要求出期望值E，在現(xiàn)實(shí)運(yùn)算中不容易實(shí)現(xiàn)，為此可通過下式進(jìn)行估計(jì): 用以上方法獲得最佳W*的運(yùn)算量很大，對(duì)于一些在線或?qū)崟r(shí)應(yīng)用的場(chǎng)合，無法滿足其時(shí)間要求。大多數(shù)場(chǎng)合使用迭代算法，對(duì)每次采樣值就求出較佳權(quán)系數(shù)，稱為采樣值對(duì)采樣值迭代算法。迭代算法可以避免復(fù)雜的和P的運(yùn)算,又能實(shí)時(shí)求得近似解，因而切實(shí)可行。LMS算法是以最快下降法為原則的迭代算法,即矢量

11、是矢量按均方誤差性能平面的負(fù)斜率大小調(diào)節(jié)相應(yīng)一個(gè)增量:，這個(gè)是由系統(tǒng)穩(wěn)定性和迭代運(yùn)算收斂速度決定的自適應(yīng)步長(zhǎng)。為次迭代的梯度。對(duì)于LMS算法為下式的斜率：由上式產(chǎn)生了求解最佳權(quán)系數(shù)W的兩種方法，一種是最陡梯度法.其思路為:設(shè)計(jì)初始權(quán)系數(shù),用迭代公式計(jì)算,到與誤差小于規(guī)定范圍。其中計(jì)算可用估計(jì)值表達(dá)式：上式K取值應(yīng)足夠大。如果用瞬時(shí)來代替上面對(duì)的估計(jì)運(yùn)算,就產(chǎn)生了另一種算法-隨機(jī)梯度法,即WidrowHoff的LMS算法。此時(shí)迭代公式為：以后討論的LMS算法都是基于Widrow-Hoff的LMS算法.上式的迭代公式假定濾波器結(jié)構(gòu)為橫向結(jié)構(gòu)。對(duì)于對(duì)稱橫向型結(jié)構(gòu)也可推出類似的迭代公式:4 自適應(yīng)濾波

12、算法的理論仿真使用Matlab編程,采用自適應(yīng)濾波器技術(shù)實(shí)現(xiàn)信噪分離,也就是去噪。程序如下:lear all；clf;signal = cos(2pi*0.0550：10001c）； %產(chǎn)生信號(hào)noise=randn（1,1000); 產(chǎn)生噪聲nfilt=fir1(11,0.4）; 11階低通FIR濾波器fnoise=filter(nfilt，1,noise)； %相關(guān)噪聲數(shù)據(jù)d=signal.c+fnoise； 信號(hào)和噪聲疊加w0 = nfilt。c0.01; %設(shè)置初始化濾波器系數(shù)mu = 0。05; %設(shè)置算法的步長(zhǎng)s=initse(w0,mu)； adaptse函數(shù)的初始化參數(shù)y，e,

13、s = adaptse(noise,d,s）；plot(0:199，signal(1：200),0：199，e（1：200)；程序運(yùn)行的結(jié)果如圖3所示。通過CCS軟件環(huán)境，把濾波程序燒錄到DSP芯片中，在CCS DSK C5000環(huán)境下輸出仿真結(jié)果:輸入信號(hào)為余弦信號(hào)和隨機(jī)噪聲的疊加.程序正確運(yùn)行后,觀察運(yùn)行結(jié)果,得出如圖4,圖5所示的仿真圖。圖3 Matlab 仿真圖圖4 余弦信號(hào)和隨機(jī)噪聲疊加的輸入信號(hào)圖5 輸出信號(hào)5 結(jié)束語(yǔ)通過仿真實(shí)驗(yàn)結(jié)果表明：自適應(yīng)濾波器卻能很好地消除疊加在信號(hào)上的噪聲，雖然也可以用固定濾波器來實(shí)現(xiàn)，但設(shè)計(jì)固定濾波器時(shí)需要預(yù)先知道信號(hào)和噪聲的統(tǒng)計(jì)特性,而自適應(yīng)濾波器則

14、不需要,并且當(dāng)信號(hào)和噪聲的統(tǒng)計(jì)特性發(fā)生變化時(shí)，自適應(yīng)濾波器也能自動(dòng)地調(diào)節(jié)其沖激響應(yīng)特性來適應(yīng)新的情況，因此，自適應(yīng)濾波器具有更加廣闊的應(yīng)用前景。參考文獻(xiàn)1胡廣書.數(shù)字信號(hào)處理理論、算法與實(shí)現(xiàn)M.北京：清華大學(xué)出版社,1998。2王世一.數(shù)字信號(hào)處理M.北京：北京理工大學(xué)出版社,1999.3程衛(wèi)國(guó)，馮峰，姚東,徐昕.Matlab啟用指南M。北京:人民郵電出版社，2000。4鄧重一。濾波器的過去、現(xiàn)在與未來J.世界電子元器件，2003(4)：48-50.5侯志林。過程控制與自動(dòng)化儀表M。北京:機(jī)械工業(yè)出版社,1998。6譚浩強(qiáng)，吳功宜.計(jì)算機(jī)教程M.2版。北京：電子工業(yè)出版社,2001。7于海生。微型計(jì)算機(jī)控制技術(shù)M.北京：清華大學(xué)出版社，1998.8程佩青。數(shù)字信號(hào)處理教程M。北京:清華大學(xué)出版社,2005。9 Simon Haykin。 Adaptive Filter TheoryM.Beijing: Publishing House of Electronics Industry