基于AVR的數字濾波器濾除工頻干擾的快速算法實現

工頻干擾廣泛存在各種工業現場中，其產生的途徑主要包括輸電饋線、照明設備、發動機以及各種電子儀器設備等。一般可以通過濾波電路消除工頻干擾，但這必將增加硬件結構的復雜程度。實際上，還可以采用數字信號處理的相關算法，通過軟件濾波器濾除工頻干擾。軟件濾波算法的采用，無疑會在簡化電路結構的同時，使系統的硬件資源得到更加充分的利用，并達到降低產品成本的要求。

AVR單片機是Atmel公司生產的8位精簡指令集(RISC)單片機。與同類單片機相比，在運算速度、外設資源、靈活性等方面性能均衡，性價比較高。AVR單片機適合C語言開發，Mega系列AVR單片機還有一個內部硬件乘法器單元。這些特點都為軟件濾波器的實現提供了極大的便利。

1 濾波算法

常用的單片機濾波算法包括中值濾波、均值濾波、滑動平均以及復合濾波算法等。工頻干擾的頻率范圍在50Hz附近，可以采用一個截止頻率遠低于50Hz的低通濾波器來濾除工頻干擾。

假設輸入信號為x(t)，輸出信號為y(t)，則一個RC低通濾波器表達式為：

連續時間信號經過采樣后成為離散時間信號，低通濾波器的表達式也變為：

所以系統的傳遞函數為：

假設采樣頻率Fs=500Hz，a分別取0.8、0.85、0.9、0.95，代入(4)式，利用matlab畫出頻率響應曲線，如圖1。其中50Hz頻率對應的幅度衰減見表1。

2 定點小數表達方式

為了精確構造數字濾波器，經常要用到浮點數據和系數。在進行浮點數乘法運算時，針對AVR單片機設計的C編譯器例如AVR-GCC，需要加入額外的數學庫函數進行編譯。而這會使編譯后程序的代碼量增加、處理時間加長、處理器的開銷也隨之增加。為了更大限度地降低系統開銷，提高程序效率，采用定點小數表示形式進行乘法運算是最佳選擇。

AVR單片機是一種8位精簡指令集(RISC)單片機。其中megaAVR系列內部都帶有一個硬件乘法器，計算一次8位乘8位的定點乘法只需2個時鐘周期。因此采用8位定點采樣數據乘以8位系數的定點乘法方式完成濾波算法是最高效的。

低通濾波處理涉及的運算形式為一個純小數系數和一個已知數據相乘再相加。因此將系數采用定點小數的表示形式，對于提高算法速度是至關重要的。

可以定義一種8位定點小數表示形式——Q8數，其各位權系數如下：

Q8數的表示范圍從0到1-2-8=0.99609375，每兩個數之間的間隔是2-8(0.00390625)，其所能表達的純小數共有28=256個。例如11011000就表示2-1+2-2+2-4+2-5=0.84375，而11011001就是表示2-1+2-2+2-4+2-5+2-8=0.84765625，因此0.84375和0.84765625之間的純小數只能用這兩個數中的一個近似表示了。這對于乘法計算的精度有一定的影響，但是由于濾波公式(3)中的系數a和(1-a)都是常數，在整體性能穩定的情況下，系數微小的不確定性對濾波器整體性能并沒有太大的影響。