摘 要： 詳細(xì)討論了高斯濾波器在單片機系統(tǒng)中的快速實現(xiàn)方法，并給出了對于MCS——51系列單片機的具體實現(xiàn)程序，介紹的方法在實時控制、信號檢測與處理方面有很大的實用價值。
關(guān)鍵詞： 濾波器 快速實現(xiàn) 單片機 實時系統(tǒng)
?
濾波器在信號處理、信號檢測、通信領(lǐng)域有非常重要的應(yīng)用，在實時系統(tǒng)中，對濾波器的性能和處理速度有非常嚴(yán)格的要求，特別是快速實時系統(tǒng)中，處理速度至關(guān)重要。目前，為滿足快速處理的需要，用DSP技術(shù)是理想的選擇。但是，目前在實時控制系統(tǒng)中，大多是用單片機實現(xiàn)的，它不僅完成信號的采樣，還需完成信號的處理和控制等功能，如果單片機系統(tǒng)本身可以完成信號的快速處理任務(wù)，將非常方便，我們在一個用MCS－51單片機組成的強噪聲背景下的通信系統(tǒng)中，實現(xiàn)了高斯濾波器的快速實現(xiàn)，滿足了系統(tǒng)的需要。
1 算法原理
高斯濾波器是一個低通濾波器，其方程
，可以證明，高斯濾波器可用均值濾波器多次逼近，一般情況下，大于或等于三次逼近就可近似于高斯濾波器，所以，在設(shè)計高斯濾波器時，可以用設(shè)計均值濾波器逼近代替高斯濾波器。
式中n為當(dāng)前采樣點序號，N為窗口寬度，顯然均值濾波器實際上可由一次加法運算，一次減法運算和一次除法運算完成，而與窗口的寬度無關(guān)，若取窗口寬度N＝2k，則除法運算可用移位來取代。
(1)式中，兩邊同乘N可變?yōu)椋?BR> 
即為了避免除法運算，我們可先采用累加機來代替均值，運算結(jié)果，再除以No上述濾波器結(jié)構(gòu)可用圖1表示。
2 均值濾波器的MCS－51快速實現(xiàn)方法
對于穩(wěn)定信號，上述計算非常簡單，它只需窗口在固定存儲器上滑動即可，對于實時信號，存儲器的內(nèi)容是隨時更新，相當(dāng)于信號經(jīng)過一個固定濾波器。
設(shè)采樣信號循環(huán)存儲在RAM addrl到addrn的存儲空間中，如圖2所示。
顯然，最新采樣值只需取代最早存放在RAM中的數(shù)據(jù)即可，而2式中的累加和實際上就是原累加和加上最新采樣值，減去最早采樣值。我們用一個數(shù)據(jù)指針指向當(dāng)前存放數(shù)據(jù)的存儲單元，為方便編程，把下一個單元的內(nèi)容作為最早采樣值，如圖2中n+1時刻的f(1)，這樣窗口寬度為n+1。
設(shè)R0為當(dāng)前數(shù)據(jù)指針，R7R6存放累加和，當(dāng)前采樣值在A中，則89C51的匯編實現(xiàn)程序為：
mov ＠R0,A,
INC R0
CJNE R0,#addrn,NEXT
mov R0,#addr0
NEXT: add A,R6
mov R6,A
mov A,R7
addc A, #00H
mov R7,A
CLR C
Mov A, R6
SUBB A, ＠R0
mov R6,A
mov A, R7
SUBB A, #00H
mov R7,A
設(shè)窗口寬度為N＝2k，上述累加和只需右移k位即為均值。
顯然，上述算法不管窗口大小如何，其速度都一樣快，所以，本算法對寬窗口濾波更為有效。
上面程序如采用89系列單片機，用24MHz晶振，只需幾個μs，如需速度更快，改用16位單片機，累加和在一個寄存器中，還可減少加法和減法的時間。
作者用上述方法在照明線數(shù)據(jù)通信系統(tǒng)中，對二值信號進(jìn)行實時處理，由于窗口寬度小于256，累加和中只用一個字節(jié)，所以，處理時間只需5μs，令人非常滿意。
本文討論了高斯濾波器在單片機系統(tǒng)中的快速實現(xiàn)。這種方法，特別適用于實時數(shù)據(jù)采集、處理、控制系統(tǒng)中的濾波。實踐表明，這種方法具有很高的實用價值，值得推廣。
參考文獻(xiàn)
1 劉松強.數(shù)字信號處理系統(tǒng)及其應(yīng)用.北京：清華大學(xué)出版社，1996.9
2 孫涵芳.MCS－51/96系統(tǒng)單片機原理及應(yīng)用.北京：北京航空航天大學(xué)出版社，1988.2